Lithostratigraphic Inversion: AVO Analysis, Rock Properties, and Pore Fluid Discrimination from Seismic
Lithostratigraphic inversion is a quantitative seismic interpretation technique that uses the variation in reflected seismic amplitude as a function of source-receiver offset (amplitude-versus-offset or AVO behaviour) to recover estimates of subsurface lithology, rock-physics properties, and pore-fluid distribution for individual rock layers rather than simply imaging reflector geometry. The method extends conventional seismic inversion, which recovers acoustic impedance from stacked seismic data, by including the angle-dependent reflection coefficients described by the Zoeppritz equations or their Aki-Richards linearization, enabling separation of P-wave impedance, S-wave impedance, and density into discrete rock-property estimates that can then be linked to porosity, mineralogy, and pore-fluid saturation through calibrated rock-physics models. In the Western Canadian Sedimentary Basin lithostratigraphic inversion has become a foundational element of unconventional resource appraisal in the Montney and Duvernay shales, where the technique is used to map sweet spots characterized by high brittleness, low clay volume, and elevated hydrocarbon saturation that cannot be discriminated using conventional post-stack inversion alone. The data inputs are pre-stack gathered seismic (typically processed to common-image-point gathers preserving relative amplitudes across the available offset or angle range), well log data calibrated to rock-physics templates, and a low-frequency model that constrains the absolute impedance level outside the seismic bandwidth. Pre-stack processing must preserve true amplitude relationships, requiring careful attention to surface-consistent amplitude balancing, residual moveout correction, and angle-mute design. Outputs include 3D volumes of P-impedance, S-impedance, density, Vp/Vs ratio, and rock-physics-derived attributes such as lambda-rho, mu-rho, Young's modulus, and Poisson's ratio that map to engineering parameters required for hydraulic fracture design and completions optimization under AER Directive 083. Typical 3D lithostratigraphic inversion projects for a Montney pad development covering 50 to 200 square kilometres run $250,000 to $750,000 CAD inclusive of pre-stack processing, well-tie work, low-frequency model building, inversion, and rock-physics analysis. The technique was pioneered in the 1990s through SPE and SEG publications by researchers including Bill Goodway, Brian Russell, and Dan Hampson, and has since matured into commercial software platforms including CGG Jason, Hampson-Russell, Schlumberger ISIS, and DUG Insight that are routinely deployed by WCSB operators including Tourmaline, ARC Resources, ConocoPhillips, and Ovintiv.
Key Takeaways
- Mathematical basis: Lithostratigraphic inversion solves the Zoeppritz equations or their Aki-Richards linearization in reverse, taking observed reflection amplitudes at multiple offsets or angles and inverting them to estimate P-impedance, S-impedance, and density at each layer interface. Modern workflows use Bayesian or deterministic inversion algorithms with well-log conditioning and stochastic rock-physics templates to produce calibrated property volumes with quantified uncertainty.
- Pre-stack data requirement: Conventional post-stack inversion can recover only acoustic impedance and works on stacked seismic data. Lithostratigraphic inversion requires offset or angle gathers preserved through processing, with careful relative-amplitude preservation. Pre-stack data volumes are 10 to 30 times larger than post-stack data, requiring substantial compute resources for processing and inversion runs that can take days to weeks on modern HPC clusters.
- Rock-physics calibration: The inversion output is only as good as the rock-physics model linking elastic properties to lithology, porosity, and pore-fluid saturation. Calibration uses sonic and density logs from offset wells, with the Hashin-Shtrikman, Voigt-Reuss-Hill, or self-consistent approximations defining how mineralogy and porosity combine to predict elastic moduli. Gassmann fluid substitution then converts dry-rock moduli to saturated values for any assumed pore-fluid mixture.
- Montney application: WCSB Montney appraisal commonly uses lithostratigraphic inversion to map sweet spots within the Upper, Middle, and Lower Montney informal members. Brittleness ratio (Young's modulus divided by Poisson's ratio) and lambda-rho versus mu-rho cross-plots discriminate quartz-rich brittle layers (preferred frac targets) from clay-rich ductile layers (frac-barrier candidates). The map products then guide horizontal well landing and frac stage design across 4,500 to 6,000 metre laterals.
- Cost and turnaround: A typical 3D lithostratigraphic inversion project over 100 square kilometres of Montney acreage runs $400,000 to $600,000 CAD across 12 to 16 weeks. Inputs are pre-existing 3D seismic (already shot at $30,000 to $60,000 CAD per square km for new acquisition), 4 to 12 calibration wells with sonic and density log suites, and a starting rock-physics model. Outputs feed directly into well-design software like Petrel or Kingdom for landing and frac plan optimization.
AVO Class Discrimination and Hydrocarbon Indicators
The five Rutherford-Williams AVO classes (Class I through IV plus Class III conformable-sand) describe how reflection amplitude varies with offset for different rock and fluid combinations. Class III anomalies (negative intercept becoming more negative with offset) are the classical bright-spot DHI signature for gas sands; Class II (small intercept with strong gradient) characterizes many WCSB Montney and Duvernay reservoirs where the elastic contrast against the host shale is subtle. Lithostratigraphic inversion converts these AVO classes into quantitative impedance and density volumes, enabling direct mapping of probable hydrocarbon-bearing intervals rather than the qualitative bright-spot interpretation that AVO offers alone.
Low-Frequency Model and Well-Tie Discipline
Seismic data lacks frequencies below approximately 5 to 8 Hz, but absolute impedance requires the full bandwidth from 0 Hz up to the seismic Nyquist limit. The missing low-frequency component is supplied by a model built from interpreted horizons and well-log filter-back impedance trends. Errors in the low-frequency model produce errors in absolute impedance that propagate directly into rock-property estimates. Best-practice WCSB workflows use 4 to 12 calibration wells distributed across the survey area, with each well tied to seismic using a synthetic seismogram derived from sonic and density logs and a well-calibrated wavelet estimated from the seismic-to-well comparison.
Fast Facts
The most cited paper in WCSB lithostratigraphic inversion practice is Bill Goodway's 1997 SEG abstract introducing the lambda-rho and mu-rho cross-plot for fluid and lithology discrimination, which was developed using data from the Alberta Plains and became known as the LMR method (Lambda-Mu-Rho). LMR cross-plot interpretation is now embedded in nearly every commercial inversion package and remains the de facto standard for Montney and Duvernay sweet-spot mapping more than 25 years after its initial publication at the SEG Annual Meeting in Dallas.
Related Terms
Lithostratigraphic inversion is closely related to the broader category of seismic inversion, which includes post-stack acoustic-impedance inversion as a simpler subset. The discipline depends entirely on accurate rock physics models that translate elastic properties into lithology and fluid estimates, and the output volumes are routinely cross-validated against well log measurements at calibration wells. The technique is also closely tied to AVO analysis, since AVO is the underlying phenomenon that pre-stack inversion exploits to separate P-impedance, S-impedance, and density from observed reflection amplitudes.
Real-World WCSB Scenario: Duvernay Sweet-Spot Mapping
An operator developing Duvernay acreage near Kaybob, Alberta commissioned a lithostratigraphic inversion of an existing 180 square kilometre 3D seismic survey in 2022 to refine sweet-spot mapping ahead of a 14-well pad development. Pre-stack gathers were re-processed for amplitude preservation, 8 offset wells with full triple-combo log suites were used to build a calibrated rock-physics model spanning the Duvernay, Ireton, and Leduc intervals, and the inversion ran on commercial software over 14 weeks for total cost of $520,000 CAD. Outputs included P-impedance, S-impedance, density, brittleness index, and lambda-rho-mu-rho cross-plots integrated into the geological model.
The resulting maps shifted the planned lateral landing of 9 of 14 wells by 4 to 11 metres vertically and 80 to 220 metres laterally relative to the pre-inversion plan, targeting brittle high-organic zones identified in the inversion volume. Post-drill production data from the first 6 wells showed average IP-30 production 18 percent above the offset average, a payback that recovered the inversion investment in 6 weeks of incremental production.