Lame Constant: Lambda and Mu Elastic Moduli, Geomechanics, and Montney Frac Modeling
The Lame constants are the two independent elastic parameters that fully describe a linear, isotropic, homogeneous solid in the framework of classical elasticity, named for French mathematician and engineer Gabriel Lame (1795 to 1870). They are conventionally written as lambda (λ, the first Lame parameter, sometimes called Lame's first constant or the modulus of dilatation) and mu (μ, the second Lame parameter, identical to the shear modulus G). Together λ and μ relate the stress tensor σ_ij to the strain tensor ε_ij through the constitutive equation σ_ij = λ · ε_kk · δ_ij + 2μ · ε_ij, where ε_kk is the trace of the strain tensor (volumetric strain) and δ_ij is the Kronecker delta. From these two constants every other isotropic elastic modulus can be derived: Young's modulus E = μ(3λ + 2μ)/(λ + μ), Poisson's ratio ν = λ/(2(λ + μ)), bulk modulus K = λ + (2/3)μ, and the P-wave modulus M = λ + 2μ. The second-form expression quoted in the original SLB-style definition (λ = K − (2/3)μ) is simply the rearrangement of K = λ + (2/3)μ to solve for λ given K and μ measured directly from log or core data. In petroleum geomechanics the Lame constants are foundational to hydraulic-fracture modeling, sanding prediction, wellbore stability analysis, 4D seismic time-lapse interpretation, and pore-pressure prediction in deep WCSB plays. Modern crosswell-tomography and full-waveform sonic logs (Schlumberger Sonic Scanner, Halliburton XBAT) record compressional (Vp) and shear (Vs) velocities, which together with density (ρ_b from RHOB log) yield μ = ρ_b · Vs² and λ = ρ_b · (Vp² − 2 · Vs²). In Montney and Duvernay siltstone-shale completions, typical values are λ of roughly 12 to 22 GPa and μ of 8 to 16 GPa, giving Young's modulus of 25 to 45 GPa and Poisson's ratio of 0.22 to 0.28; these values feed directly into Mangrove, GOHFER, ResFrac, and Petrel-RE frac simulators to predict fracture half-length, height growth into the bounding Doig phosphate, and proppant-placement efficiency. AER Directive 083 (Hydraulic Fracturing Subsurface Integrity) and BC OGC's IB 2018-08 both require documented geomechanical inputs, including Lame-derived stress profiles, in completion permits for wells within 200 m of an offset wellbore. Operators including Canadian Natural Resources Limited, ARC Resources, and Tourmaline Oil publish Lame-based stress models in their NI 51-101 reserves disclosures and use them to optimize multi-well pad spacing and fracture-stage clustering across the WCSB unconventional fairway.
Key Takeaways
- Two independent elastic constants: Lambda (λ) and mu (μ) together fully describe a linear isotropic elastic solid. Every other elastic modulus (Young's E, Poisson's ν, bulk K, P-wave M) is derivable from the pair. Mu (μ) is identical to the shear modulus G; lambda (λ) has no direct physical analog but governs the coupling between volumetric strain and normal stress.
- Computed from sonic and density logs: Modern dipole sonic logs deliver compressional (Vp) and shear (Vs) slownesses, and bulk-density logs deliver ρ_b. Then μ = ρ_b · Vs² and λ = ρ_b · (Vp² − 2 · Vs²), with values in pascals when ρ_b is in kg/m³ and velocities in m/s. Routine WCSB practice computes λ and μ at every 0.15 m sample point along a logged interval.
- Frac-design input: Hydraulic-fracture simulators (Mangrove, GOHFER, ResFrac, Petrel-RE) require Lame-derived stress profiles to predict fracture geometry, height growth, and proppant placement. In Montney and Duvernay completions, λ of 12 to 22 GPa and μ of 8 to 16 GPa give Young's modulus 25 to 45 GPa and Poisson's ratio 0.22 to 0.28, the working stress envelope for stage design across the WCSB unconventional fairway.
- LMR attribute interpretation: Pre-stack seismic inversion outputs λρ (lambda-rho) and μρ (mu-rho) attributes (the Goodway, Chen, and Downton 1997 LMR method), which separate brittle, low-Poisson "sweet spot" zones from ductile shale baffles. In the Montney, μρ greater than 25 GPa·g/cc and λρ less than 30 GPa·g/cc routinely flag the most productive landing intervals across the Karr-Kakwa and Pipestone fairways.
- Regulatory frame: AER Directive 083 and BC OGC IB 2018-08 require geomechanical documentation in completion permits for wells within 200 m of offset wellbores. Lame-based stress profiles, in-situ stress estimates, and frac-barrier predictions are core deliverables, and inadequate documentation can delay frac permits by 4 to 12 weeks, with associated standby and rig-rate costs of CAD $25,000 to $80,000 per day.
Computing Lambda and Mu from a Dipole Sonic Log
On a Montney horizontal logged with a Schlumberger Sonic Scanner, the petrophysicist exports compressional slowness DTC (typically 65 to 80 µs/ft in the Montney C interval) and shear slowness DTS (110 to 145 µs/ft) alongside bulk density RHOB (2.55 to 2.68 g/cc). Converting slowness to velocity (V = 304,800/DT in m/s, with DT in µs/ft) yields Vp of 3,810 to 4,690 m/s and Vs of 2,100 to 2,770 m/s. With ρ_b of 2,620 kg/m³, μ comes to 11.5 to 20.1 GPa and λ to 13.8 to 24.5 GPa. These per-foot values feed the geomechanical earth model.
Why the Two Constants Cannot Be Reduced to One
An isotropic elastic solid has two independent elastic moduli; no single number captures its mechanical response under arbitrary loading. A material with high μ (stiff in shear) but low λ deforms very differently under triaxial stress than one with the reverse, even at identical Young's modulus. In a Duvernay completion, a low-λ, high-μ interval is brittle and fractures readily, while a high-λ, low-μ shale is ductile, absorbs fracture energy, and acts as a frac barrier. The two-constant nature of isotropic elasticity is, in this sense, exactly what makes the LMR pre-stack inversion useful for landing-zone selection.
Fast Facts
Gabriel Lame published his Lecons sur la theorie mathematique de l'elasticite des corps solides in Paris in 1852, formalizing the two-constant framework that has carried unchanged for 172 years into modern petroleum geomechanics. Lame himself never saw a hydraulic-fracture simulator, but his constants determine whether a CAD $14 million Montney horizontal frac job places proppant productively in the C interval or wastes 30 to 50 percent of its energy growing height into the Doig phosphate, a daily reality for every WCSB completions engineer in 2026.
Related Terms
The Lame constant μ is identical to the shear modulus, the resistance of a solid to shape change under shear stress. Lame's λ feeds directly into the bulk modulus K = λ + (2/3)μ, the resistance to volumetric compression. Together they derive Poisson's ratio, the ratio of transverse to axial strain, which controls fracture geometry and proppant placement in unconventional completions. The Lame parameters are the primary inputs to a geomechanical model used in WCSB Montney, Duvernay, and Cardium completions to predict in-situ stress, frac height, and lateral spacing.
Montney Frac-Stage Design Scenario
A completions engineer at a WCSB Montney operator is designing a 60-stage, 3,200 m lateral frac in the Pipestone fairway at a budgeted CAD $14.2 million completed-well cost. Dipole-sonic-derived Lame constants across the C interval show μ averaging 14.8 GPa and λ averaging 17.2 GPa, yielding Young's modulus of 36 GPa and Poisson's ratio of 0.24 (clearly brittle and frac-favourable). The overlying Doig phosphate shows μ of 8.5 GPa and λ of 23 GPa, Poisson's ratio of 0.32 (ductile, a good upward frac barrier). The Mangrove model predicts 92 m frac half-lengths and 38 m height growth, contained below the Doig.
The completion is executed to design and the well delivers a first-six-month cumulative of 142,000 boe versus a 125,000 boe type-curve expectation, a 14 percent outperformance attributed directly to the Lame-constant-driven landing-zone selection and stage-spacing design, and the same workflow is rolled out across the operator's CAD $260 million 2027 Pipestone development program.