Seismic inversion
Seismic inversion, in Geophysics (primarily Oil and Gas exploration/development), is the process of transforming seismic reflection data into a quantitative rockproperty description of a reservoir. Seismic inversion may be preor poststack, deterministic, random or geostatistical, and typically includes other reservoir measurements such as well logs and cores.
Contents
Introduction
Seismic surveys are routinely performed to gather information about the geology of an oil or gas field. These surveys record sound waves which have traveled through the layers of rock and fluid in the earth. The amplitude and frequency of these waves can be estimated so that any sidelobe and tuning effects^{[2]} introduced by the wavelet may be removed.
Seismic data may be inspected and interpreted on its own without inversion, but this does not provide the most detailed view of the subsurface and can be misleading under certain conditions. Because of its efficiency and quality, most oil and gas companies now use seismic inversion to increase the resolution and reliability of the data and to improve estimation of rock properties including porosity and net pay.^{[3]}
There are many different techniques used in seismic inversion.^{[4]} These can be roughly grouped into two categories: prestack or poststack, and seismic resolution or welllog resolution. The combination of these categories yields four technical approaches to the inversion problem, and the selection of a specific technique depends on the desired objective and the characteristics of the subsurface rocks. Although the order presented reflects the advancement in inversion techniques over the past 20 years, each grouping still has valid uses in particular projects or as part of a larger workflow.
Wavelet estimation
All modern seismic inversion methods require seismic data and a wavelet estimated from the data. Typically, a reflection coefficient series from a well within the boundaries of the seismic survey is used to estimate the wavelet phase and frequency. Accurate wavelet estimation is critical to the success of any seismic inversion. The inferred shape of the seismic wavelet may strongly influence the seismic inversion results and, thus, subsequent assessments of the reservoir quality.
Wavelet amplitude and phase spectra are estimated statistically from either the seismic data alone or from a combination of seismic data and well control using wells with available sonic and density curves. After the seismic wavelet is estimated, it is used to estimate seismic reflection coefficients in the seismic inversion.
When the estimated (constant) phase of the statistical wavelet is consistent with the final result, the wavelet estimation converges more quickly than when starting with a zero phase assumption. Minor edits and "stretch and squeeze" may be applied to the well to better align the events. Accurate wavelet estimation requires the accurate tie of the impedance log to the seismic. Errors in well tie can result in phase or frequency artifacts in the wavelet estimation. Once the wavelet is identified, seismic inversion computes a synthetic log for every seismic trace. To ensure quality, the inversion result is convolved with the wavelet to produce synthetic seismic traces which are compared to the original seismic.^{[4]}
Components of inversion
Inversion includes both seismic data and well data, where well data serves to add the low frequency below the seismic band and to constrain the inversion. Well logs are first conditioned and edited to ensure there is a suitable relationship between impedance logs and the desired properties. The logs are then converted to time, filtered to approximate the seismic bandwidth and edited for borehole effects, balanced and classified by quality.
Seismic data is bandlimited, reducing resolution and quality. To extend the frequency band available, lowfrequency data is derived from log data, prestack depth or time migrated velocities and/or a regional gradient.^{[5]} High frequency can be derived from well control or geostatistical analysis.
Initial inversions are often run with relaxed constraints, starting with the seismic and then adding limitedtrend data from the wells. This provides a rough overview of the reservoir in an unbiased manner. It is critical at this point to evaluate the accuracy of the tie between the inversion results and the wells, and between the original seismic data and the derived synthetics. It is also important to ensure that the wavelet matches the phase and frequency of seismic data.
Without a wavelet, the solution is not unique. Deterministic inversions address this problem by constraining the answer in some way, usually to well log data. Stochastic inversions address this problem by generating a range of plausible solutions, which can then be narrowed through testing for best fit against various measurements (including production data).
Poststack seismic resolution inversion
An example of a poststack seismic resolution inversion technique is the Constrained SparseSpike Inversion (CSSI). This assumes a limited number of reflection coefficients, with larger amplitude. The inversion results in acoustic impedance (AI), which is the product of rock density and pwave velocity. Unlike seismic reflection data (which is an interface property) AI is a rock property. The model generated is of higher quality, and does not suffer from tuning and interference caused by the wavelet.
CSSI transforms seismic data to a pseudoacoustic impedance log at every trace. Acoustic impedance is used to produce more accurate and detailed structural and stratigraphic interpretations than can be obtained from seismic (or seismic attribute) interpretation. In many geological environments acoustic impedance has a strong relationship to petrophysical properties such as porosity, lithology, and fluid saturation.
A good (CSSI) algorithm will produce four highquality acoustic impedance volumes from full or poststack seismic data: fullbandwidth impedance, bandlimited Impedance, reflectivity model, and lowfrequency component. Each of these components can be inspected for its contribution to the solution and to check the results for quality. To further adapt the algorithm mathematics to the behavior of real rocks in the subsurface, some CSSI algorithms use a mixednorm approach and allow a weighting factor between minimizing the sparsity of the solution and minimizing the misfit of the residual traces.
Prestack seismic resolution inversion
Prestack inversion is often used when poststack inversion fails to sufficiently differentiate geologic features with similar Pimpedance signatures.^{[6]} Simultaneous inversion solves for Simpedance and density, in addition to Pimpedance. While many geologic features can express similar Pimpedance characteristics, few will share combined Pimpedance and Simpedance traits (allowing improved separation and clarity). Often a feasibility study using the wells logs will indicate whether separation of the desired lithotype can be achieved with Pimpedance alone or whether Simpedance is also required. This will dictate whether a pre or poststack inversion is needed.
Simultaneous Inversion (SI) is a prestack method that uses multiple offset or angle seismic substacks and their associated wavelets as input; it generates Pimpedance, Simpedance and density as outputs (although the density output resolution is rarely as high as the impedances). This helps improve discrimination between lithology, porosity and fluid effects. For each input partial stack, a unique wavelet is estimated. All models, partial stacks and wavelets are input to a single inversion algorithm —enabling inversion to effectively compensate for offsetdependent phase, bandwidth, tuning and NMO stretch effects.^{[7]}
The inversion algorithm works by first estimating angledependent Pwave reflectivities for the inputpartial stacks. Next, these are used with the full Zoeppritz equations (or approximations, such as AkiRichards, for some algorithms) to find bandlimited elastic reflectivities. These are in turn merged with their lowfrequency counterparts from the model and integrated to elastic properties. This approximate result is then improved in a final inversion for Pimpedance, Simpedance and density, subject to various hard and soft constraints. One constraint can control the relation between density and compressional velocity; this is necessary when the range of angles is not great enough to be diagnostic of density.
An important part in the inversion procedure is the estimation of the seismic wavelets. This is accomplished by computing a filter that best shapes the angledependent well log reflection coefficients in the region of interest to the corresponding offset stack at the well locations. Reflection coefficients are calculated from Psonic, Ssonic and density logs using the Zoeppritz equations. The wavelets, with amplitudes representative of each offset stack, are input directly into the inversion algorithm. Since a different wavelet is computed for each offset volume, compensation is automatically done for offsetdependent bandwidth, scaling and tuning effects. A nearstack wavelet can be used as the starting point for estimating the farangle (or offset) wavelet.
No prior knowledge of the elastic parameters and density beyond the solution space defined by any hard constraints is provided at the well locations. This makes comparison of the filtered well logs and the inversion outputs at these locations a natural quality control. The lowest frequencies from the inversion are replaced with information from the geologic model since they are poorly constrained by the seismic data. When applied in global mode a spatial control term is added to the objective function and large subsets of traces are inverted simultaneously. The simultaneous inversion algorithm takes multiple anglestacked seismic data sets and generates three elastic parameter volumes as output.
The resulting elastic parameters are realrock properties that can be directly related to reservoir properties. The more advanced algorithms use the full KnottZoeppritz equations and there is full allowance for amplitude and phase variations with offset. This is done by deriving unique wavelets for each inputpartial stack. The elastic parameters themselves can be directly constrained during the seismic inversion and rockphysics relationships can be applied, constraining pairs of elastic parameters to each other. Final elasticparameter models optimally reproduce the input seismic, as this is part of the seismic inversion optimization.
Poststack log detail inversion
Geostatistical inversion integrates high resolution well data with low resolution 3D seismic, and provides a model with high vertical detail near and away from well control. This generates reservoir models with geologicallyplausible shapes, and provides a clear quantification of uncertainty to assess risk. Highlydetailed petrophysical models are generated, ready for input to reservoirflow simulation.
Geostatistics differs from statistics in that it recognizes that only certain outcomes are geologically plausible. Geostatistical inversion integrates data from many sources and creates models that have greater resolution than the original seismic, match known geological patterns, and can be used for risk assessment and reduction.
Seismic, well logs and other input data are each represented as a probability density function (PDF), which provides a geostatistical description based on histograms and variograms. Together these define the chances of a particular value at a particular location, and the expected geological scale and composition throughout the modeled area.
Unlike conventional inversion and geomodeling algorithms, geostatistical inversion takes a onestep approach, solving for impedance and discrete property types or lithofacies at the same time. Taking this approach speeds the process and improves accuracy.
Individual PDFs are merged using bayesian inference techniques, resulting in a posterior PDF conditioned to the whole data set. The algorithm determines the weighting of each data source, eliminating potential bias. The posterior PDF is then input to a Markov Chain Monte Carlo algorithm to generate realistic models of impedance and lithofacies, which are then used to cosimulate rock properties such as porosity. These processes are typically iterated until a model emerges that matches all information. Even with the best model, some uncertainty remains. Uncertainty can be estimated using random seeds to generate a range of realizations. This is especially useful when dealing with parameters that are sensitive to change; an analysis of this sort enables greater understanding of development risk.
Prestack logdetail inversion
Amplitude versus offset (AVO) (AVA) geostatistical inversion incorporates simultaneous AVO (AVA) inversion into the geostatistical inversion algorithm so high resolution, geostatistics, and AVO may be attained in a single method. The output model (realizations) are consistent with well log information, AVO seismic data, and honor rock property relationships found in the wells. The algorithm also simultaneously produces elastic properties (Pimpedance, Simpedance and density) and lithology volumes, instead of sequentially solving for lithology first and then populating the cell with impedance and density values. Because all output models match all input data, uncertainty can be quantitatively assessed to determine the range of reservoir possibilities within the constraining data.
AVA geostatistical inversion software uses leadingedge geostatistical techniques, including Markov Chain Monte Carlo (MCMC) sampling and pluriGaussian lithology modeling. It is thus possible to exploit "informational synergies" to retrieve details that deterministic inversion techniques blur out or omit. As a result, geoscientists are more successful in reconstructing both the overall structure and the fine details of the reservoir. The use of multipleanglestack seismic volumes in AVA geostatistical inversion enables further evaluation of elastic rock properties and probable lithology or seismic facies and fluid distributions with greater accuracy.
The process begins with a detailed petrophysical analysis and well log calibration. The calibration process replaces unreliable and missing sonic and density measurements with synthesized values from calibrated petrophysical and rockphysics models. Well log information is used in the inversion process to derive wavelets, supply the low frequency component not present in the seismic data, and to verify and analyze the final results. Next, horizon and log data are used to construct the stratigraphic framework for the statistical information to build the models. In this way, the log data is only used for generating statistics within similar rock types within the stratigraphic layers of the earth.
Wavelet analysis is conducted by extracting a filter from each of the seismic volumes using the well elastic (angle or offset) impedance as the desired output. The quality of the inversion result is dependent upon the extracted seismic wavelets. This requires accurate psonic, ssonic and density logs tied to the appropriate events on the seismic data. Wavelets are extracted individually for each well. A final "multiwell" wavelet is then extracted for each volume using the best individual well ties and used as input to the inversion.
Histograms and variograms are generated for each stratigraphic layer and lithology, and preliminary simulations are run on small areas. The AVA geostatistical inversion is then run to generate the desired number of realizations, which match all the input data. The results are quality controlled by direct comparison of the inverted rock property volumes against the well logs. Further QC involves review by a multidisciplinary team of all input parameters and the results of the simulation. Analysis of multiple realizations produces mean (P50) property cubes or maps. Most often these are lithology or seismic facies cubes and predicted lithology or facies probabilities, but other outputs are also possible. Selected lithology and facies cubes are also generated for P15 and P85 probabilities (for example). Reservoir 3D bodies of hydrocarbonbearing units are captured with their corresponding rock properties, and the uncertainty in reservoir size and properties is quantified.
See also
References
 ↑ J. Pendrel, Scandinavian Oil Gas Magazine, No. 5/6 2006
 ↑ Oilfield glossary Retrieved 20110603.
 ↑ Pendrel, J., "Seismic Inversion—A Critical Tool in Reservoir Characterization", Scandinavian OilGas Magazine, No. 5/6, 2006, p. 1922.
 ↑ ^{4.0} ^{4.1} Sen, M. K., "Seismic Inversion", Society of Petroleum Engineers, 2006.
 ↑ Latimer, R., Davison, R., van Riel, P., "An Interpreter’s Guide to Understanding and Working with SeismicDerived Acoustic Impedance Data", The Leading Edge, March 2000, p. 242256.
 ↑ Pendrel, J., "Seismic Inversion—The Best Tool for Reservoir Characterization", CSEG Recorder.
 ↑ Pendrel, J., Dickson, T., "Simultaneous AVO Inversion to P Impedance and Vp/Vs", SEG.
Further reading
 Caulfield, C., Feroci, M., Yakiwchuk, K. "Seismic Inversion for Horizontal Well Planning in Western Saskatchewan", Evolving Geophysics Through Innovation, p. 213214.
 Chakrabarty, C., Fossey, J., Renard, G., Gadelle, C. "SAGD Process in the East Senlac Field: From Reservoir Characterization to Field Application", No. 1998.192.
 Contreras, A., TorresVerdin, C., Chesters, W., Kvien, K., Globe, M., "Joint Stochastic Inversion of Petrophysical Logs and 3D PreStack Seismic Data to Assess the Spatial Continuity of Fluid Units Away from Wells: Application to a GulfofMexico Deepwater Hydrocarbon Reservoir", SPWLA 46th Annual Logging Symposium, June 26–29, 2005.
 De Barros, Dietrich, M., "Full Waveform Inversion of Shot Gathers in Terms of Poroelastic Parameters", EAGE, London, June 2007.
 Deutsch, C., Geostatistical Reservoir Modeling, New York: Oxford University Press, 2002, 376 pages.
 Francis, A., "Limitations of Deterministic and Advantages of Stochastic Seismic Inversion", CSEG Recorder, February 2005, p. 511.
 Hasanusi, D., Adhitiawan, E., Baasir, A., Lisapaly, L., van Eykenhof, R., "Seismic Inversion as an Exciting Tool to Delineate Facies Distribution in Tiaka Carbonate Reservoirs, Sulawesi – Indonesia", Proceedings, Indonesian Petroleum Association, ThirtyFirst Annual Convention and Exhibition, May 2007.
 Russell, B., Hampson, D., "The Old and the New in Seismic Inversion", CSEG Recorder, December 2006, p. 511.
 Stephen, K., MacBeth, C., "Reducing Reservoir Prediction Uncertainty by Updating a Stochastic Model Using Seismic History Matching", SPE Reservoir Evaluation & Engineering, December 2008.
 VargasMeleza, L., Megchun, J., Vazquez, G., "Petrophysical Properties Estimation by Integrating AVO, Seismic Inversion and Multiattribute Analysis in a 3D Volume of Playuela, Veracruz", AAPG International Conference: October 24–27, 2004, Cancun, Mexico.
 Wang, X., Wu, S., Xu, N., Zhang, G., "Estimation of Gas Hydrate Saturation Using Constrained Sparse Spike Inversion: Case Study from the Northern South China Sea", Terr. Atmos. Ocean. Sci., Vol. 17, No. 4, 799813, December 2006.
 Watson, I., Lines, L., "Seismic Inversion at Pike’s Peak, Saskatchewan", CREWES Research Report, Volume 12, 2000.
 Zou, Y., Bentley, L., Lines, L., "Integration of Reservoir Simulation with TimeLapse Seismic Modeling", 2004 CSEG National Convention.
External links
 American Association of Petroleum Geologists
 Australian Society of Exploration Geophysicists (ASEG)
 Earthworks Seismic Inversion company
 Understanding Stochastic Inversion
 FugroJason, Jason Geoscience Workbench (JGW)
 CGGVERITAS, Simultaneous Elastic Inversion
 European Association of Geoscientists and Engineers
 Petroleum Exploration Society of Australia (PESA)
 Society of Exploration Geophysicists
 Society of Petroleum Engineers
 Society of Petrophysicists and Well Log Analysts (SPWLA)
 The University of Texas at Austin Petroleum and Geoscience Engineering Reading Room
 Geological Survey Publications Warehouse
