Elastic Constants

Key Takeaways

• Young's modulus and Poisson's ratio together determine hydraulic fracture width and closure stress — in hydraulic fracturing, the net pressure in the fracture (fracturing pressure minus closure stress) drives fracture opening, and the fracture width response to net pressure is directly proportional to 1/E (higher Young's modulus means stiffer rock that opens less per unit of net pressure); fracture width governs proppant transport (narrower fractures restrict proppant entry), fracture conductivity (wider proppant packs have higher permeability), and screenout risk (narrow fractures bridge more readily); Poisson's ratio controls the relationship between vertical and horizontal stresses through the equation sigma_h = (ν/(1-ν)) × sigma_v (for uniaxial strain conditions) — higher ν means higher horizontal stress relative to vertical stress, affecting both the closure stress and the fracture orientation; in heterogeneous formations with alternating stiff (high E) and compliant (low E) layers, fracture height growth is controlled by the elastic contrast between layers — fractures tend to be contained within stiff layers because they open less per unit net pressure and the stress intensity factor at the fracture tip is lower for the same fracture dimensions in stiffer rock.
• Elastic constants from sonic logs require core calibration to be used quantitatively in geomechanical models — dynamic elastic constants from sonic logs (Young's modulus_dynamic and Poisson's ratio_dynamic) are consistently higher than the static values measured in laboratory tests because: dynamic measurements at sonic frequencies (hundreds of Hz) are essentially undrained (pore fluid does not have time to redistribute, making the rock appear stiffer than the drained static value); dynamic measurements at micro-strain amplitudes have all microcracks closed by confining stress, while static measurements at larger strains open these cracks (making the rock appear softer in static tests); empirical relationships between dynamic and static elastic constants (E_static = a × E_dynamic + b, where a and b are regression coefficients specific to each lithology type) are established by measuring both dynamic and static elastic constants on the same core samples from the same well, then applying these corrections to the full sonic log profile to obtain continuous static elastic constants; without this core calibration, using dynamic sonic log elastic constants directly in fracture design or wellbore stability models will systematically overestimate rock stiffness.
• Anisotropic elastic constants describe formations where properties differ in different directions — real sedimentary rocks are often not isotropic (having the same elastic constants in all directions) but transversely isotropic (TI), with five independent elastic constants rather than two: the TI model has different stiffness in the direction parallel to bedding versus perpendicular to bedding, reflecting the layered fabric of the sediment; shale is the most common and strongly anisotropic sedimentary rock, with P-wave velocity parallel to bedding typically 15-30% higher than perpendicular to bedding, and correspondingly different Young's modulus values in the two directions; for wellbore stability analysis in highly deviated wells drilled through anisotropic shale, isotropic elastic constants give incorrect predictions of wellbore failure orientation and mud weight requirements because the failure criterion depends on the orientation-dependent stress and strength distribution around the wellbore; cross-dipole acoustic logging that measures shear wave velocity in two perpendicular directions provides the additional data needed to characterize TI anisotropy and improve geomechanical model accuracy.
• Elastic constants change with effective stress, temperature, and saturation, affecting logging versus reservoir conditions — in-situ elastic constants differ from laboratory measurements on core samples because: effective stress at reservoir depth may differ from the applied confining pressure in laboratory tests; reservoir temperature affects clay mineral stiffness and fluid properties that contribute to the elastic response; the pore fluid (oil, gas, or brine) affects the bulk modulus through Gassmann's equation, so a gas-saturated formation has lower bulk modulus than the same rock brine-saturated; Gassmann fluid substitution calculations allow petrophysicists to predict how elastic constants and sonic velocities change when reservoir fluid is substituted (as occurs during production when oil is replaced by water or when gas is produced from a gas cap) — essential for 4D seismic feasibility studies that require knowing whether the time-lapse seismic response will be detectable given the expected fluid saturation changes during production.
• Reservoir compaction is governed by the elastic (and plastic) constants of the reservoir rock and overburden — when reservoir pressure declines during production, effective stress on the reservoir rock increases and the rock compresses elastically (reversibly), reducing pore volume and driving some additional hydrocarbon production (compaction drive); the uniaxial compressibility (Cm = 1/E × (1+ν)(1-2ν)/(1-ν)) derived from elastic constants relates the vertical effective stress increase to the vertical strain in the reservoir, allowing calculation of reservoir compaction and the resulting surface subsidence; in reservoirs with low Young's modulus (soft formations such as chalk or unconsolidated sand), compaction can be dramatic — the Ekofisk chalk field in the Norwegian North Sea compacted several meters over decades of production, causing significant seafloor subsidence that required raising the platform decks; in stiffer formations, compaction is smaller but still contributes to production and must be included in material balance and reserve calculations; distinguishing elastic (recoverable) from plastic (permanent) compaction is essential for predicting whether compaction damage to the reservoir rock will occur during pressure depletion.