Mechanical properties are used to design hydraulic fracture stimulation programs in oil and gas wells, and in the design of mines and gas storage caverns. In this situation, the mechanical properties are derived in the laboratory or from well log analysis, calibrated to the lab results.
The main purpose for finding these attributes is to distinguish reservoir quality rock from non-reservoir. The ultimate goal is to determine porosity, lithology, and fluid type by "reverse-engineering" the seismic attributes. The process is sometimes called "quantitative seismic interpretation". In high porosity areas such as the tar sands, and in high contrast areas such as gas filled carbonates,, modest success has been achieved, usually after several iterative calibrations to log and lab data. Something can be determined in almost all reservoirs, but how "quantitative" it is may not be known.
There are many other types of
seismic attributes related to the signal frequency, amplitude, and
phase, as well as spatial attributes that infer geological structure
and stratigraphy, such as dip angle, dip azimuth, continuity,
thickness, and a hundred other factors. While logs may be used to
calibrate or interpret some of these attributes, they are not
discussed further here. Dynamic elastic constants can also be determined in the laboratory using high frequency acoustic pulses on core samples. Static elastic constants are derived in the laboratory from tri-axial stress-strain measurements (non-destructive) or the chevron notch test (destructive).
Elastic constants
are needed by five distinct disciplines in the petroleum industry: The elastic constants of rocks are defined by the Wood-Biot-Gassmann Equations. The equations can be transformed to derive rock properties from log data. If crossed dipole sonic data is available, anisotropic stress can be noticed by differences in the X and Y axis displays of both the compressional and shear travel times. When this occurs, all the elastic constants can be computed for both the minimum and maximum stress directions. This requires the original log to be correctly oriented with directional information, and may require extra processing in the service company computer center.
Elasticity is a property of matter,
which causes it to resist deformation in volume or shape.
Hooke's Law, describing the behavior of elastic materials,
states that within elastic limits, the resulting strain is
proportional to the applied stress. Stress is the external
force applied per unit area (pressure), and strain is the fractional
distortion which results because of the acting force. Different types of deformation can result, depending upon the mode of the acting force. The three elastic moduli are: Young's Modulus
Y (also abbreviated E in various literature), Bulk Modulus
Kc, Shear Modulus
N, (also abbreviated as
Poisson's Ratio
PR (also abbreviated
All of these
moduli can be derived directly from well logs and indirectly from
seismic attributes:
Some people prefer different abbreviations: Mu or In the seismic industry, it is common to think in terms of velocity and acoustic impedance in addition to the more classical mechanical properties described above.
The compressional to shear velocity ratio is a good
lithology indicator:
Acoustic impedance:
Where:
Velocity of sound, density, and elastic properties of rocks are
intimately connected by a series of transforms. Knowledge of any two
of these properties means all the others can be calculated.
The velocity of longitudinal (compressional) waves in solids can
be predicted from the following two equations. Where: The
transverse (shear) wave velocity is defined by the following two
equations: Where: To translate these formulae into metric, convert density into gm/cc, velocity to Km/sec and the various moduli to megabars, and change the constant terms to 1.0. To convert moduli in megabars to psi, multiply by 6.89 * 10^-6. To convert megabars to Kilopascals, multiply by 10^4. The elastic constants K, N, Y and PR are often known, and many values are listed in handbooks. Identities exist which show that knowledge of any two constants infers knowledge of the other two. This in turn, infers knowledge of velocity. These identities follow.
3: K = L + 2 * N / 3 4: K = Y * N / (3 * (3 * N - Y)) 5: K = L * (1 + PR) / (3 * PR) 6: K = S * (2 * (1 + PR)) / (3 * (1 - 2 * PR)) 7: K = Y / (3 * (1 - 2 * PR)) 8: K = DENS * (Vp ^ 2 - 4 / 3 * Vs ^ 2)
9: Y = N * (3 * L + 2 * N) / (L + N) 10: Y = 9 * K * (K - L) / (3 * K - L) 11: Y = 9 * K * L / (3 * K + L) 12: Y = L * (1 + PR) * (1 - 2 * PR) / PR 13: Y = 2 * N * (1 + PR) 14: Y = 3 * K * (1 - 2 * PR) 15: Y = ((9 * DENS * R3 ^ 2 * R2 ^ 2) / (3* R2 ^ 2 + 1)) Where:
18: L = K - 2 * N / 3 19: L = N * (Y - 2 * N) / (3 * N - Y) 20: L = 3 * K * (3 * K - Y) / (N * K - Y) 21: L = N * (2 * PR / (1 - 2 * PR)) 22: L = 3 * K * (PR / (1 - PR)) 23: L = Y * PR / ((1 + PR) * (1 - 2 * PR)) 24: L = DENS * (Vp^2 - 2 * Vs ^ 2)
25: PR = L / 2 * (L + N) 26: PR = L / (3 * K - L) 27: PR = (3 * K - 2 * N) / (2 * (3 * K + N)) 28: PR = (Y / (2 * N)) - 1 29: PR = (3 * K - Y) / (6 * K) 30: PR = ((R1^2 - 2) / (R1^2 - 1) / 2) 31: PR = ((3 * (R2^2) - 2) / (3 * (R2^2) + 1) / (3 * (R3^2) + 1) / 2) Where:
33: DENS = (L + 2 * N) / (Vp ^ 2) 34: DENS = (3 * K - 2) / (Vp ^ 2) 35: DENS = (K + 4 * N / 3) / (Vp ^ 2) 36: DENS = N * (4 * N - Y) / (3 * N - Y) / (Vp ^ 2) 37: DENS = 3 * K * (3 * K + Y) / (9 * K - Y) / (Vp ^ 2) 38: DENS = L * ((1 - PR) / PR) / (Vp ^ 2) 39: DENS = N * (2 - 2 * PR) / (1 - 2 * PR) / (Vp ^ 2) 40: DENS = 3 * K * (1 - PR) / (1 + PR) / (Vp ^ 2) 41: DENS = Y * (1 - PR) / ((1 + PR) * ( 1 - 2 * PR)) / (Vp ^ 2) 42: DENS = 3 * ( K - L) / 2 / (Vs ^ 2) 43: DENS = 3 * K * Y / (9 * K - Y) / (Vs ^ 2) 44: DENS = L * ((1 - 2 * PR) / (2 * PR) / Vs ^ 2) 45: DENS = 3 * K * (1 - 2 * PR) / (2 + 2 * PR) / (Vs ^ 2) 46: DENS = Y / (2 + 2 * PR) / (Vs ^ 2) Such relationships are used to reconstruct density logs in bad hole conditions by using sonic log data and assumed values for Poisson's ratio. PR is often a function of shale volume and lithology, which can be determined in zones where hole condition is good. Where:
The
first and simplest relates compressibility (which is the inverse
of the bulk modulus K) and pressure: Where: The constants a and b, for particular solids can be found in the Handbook of Physical Constants. For
example assume the following measured values on a limestone sample: K
= Y / 3 * (1 - 2 * P) = 11.42*10^6 / 3 (1 - 2 * 0.32) = 10.6 *
10^6 psi
1. VOIGHT 48: a = (C11 + C22 + C33) * 4.83*10^6 49: b = (C23 + C31 + C12) * 4.83*10^6 50: c = (C44 + C55 + C66) * 4.83*10^6 51: K = (a + 2 * b) / 3 52: N = (a - b + 3 * c) / 5 2.
REUSS Where: The Cij and Sij values are obtained from the tables in The Handbook of Physical Constants. Other coefficients for the aggregate may be obtained from K and N, by use of the relationships between the various elastic constants given earlier. Examples of these two methods are also shown in the Handbook of Physical Constants. For
many rocks, elastic constants are known, whereas velocity is unknown.
This is especially true when the effects of pressure and temperature
are being considered. It is also clear that given a reasonable
set of elastic constants and either a velocity or density log,
the other log can be constructed with confidence. This is particularly
useful in seismography. Note that the sonic velocity log as a
rule, measures the travel time associated with the longitudinal
or compressional wave. Therefore, the appropriate equations should
be used for log interpretation work. |
||

Page Views ---- Since 01 Jan 2015
Copyright 2023 by Accessible Petrophysics Ltd. CPH Logo, "CPH", "CPH Gold Member", "CPH Platinum Member", "Crain's Rules", "Meta/Log", "Computer-Ready-Math", "Petro/Fusion Scripts" are Trademarks of the Author |