Seismic inversion (post-stacking vs pre-stacking)
Seismic inversion, most of the time used in oil and gas industry, is the process used to reconstruct earth properties. It combines seismic and well data to predict rock properties (lithology, fluid content, porosity) across a survey. These rock properties can be used to identify hydrocarbon and reservoir.
- 1 Description
- 2 Post-stacking inversion
- 3 Pre-stacking inversion
- 4 Differences and Similarities
- 5 References
- 6 External link
Importance of impedance attribute
Many of these rocks properties can be identified by using well log data (gamma ray, water saturation, shale volume) or seismic data (more difficult to obtain). However, with seismic inversion we can obtain other rocks properties such as impedance and its attributes (P-impedance, S-impedance, Poisson’s Ratio, Vp/Vs, Lambda*Rho, Mu*Rho). Those properties are linked to fluid content, porosity, lithology. For example, for a given lithology, if we know the P-impedance, we can predict its porosity. From this relationship, by combining impedance at the well and the impedance calculated from seismic data, we can predict hydrocarbon across the survey.
How to calculate impedance
P-impedance= density *P-velocity
Note that all other attributes can be calculated from the impedance.
From well log data we used density and velocity data to get the impedance and the Poisson’s ratio. From Seismic data, we obtain impedance by using seismic inversion which converts seismic from a boundary property to a layer property. Indeed, well data show us properties of the rock’s layer and seismic data show us information about the boundary between rock layers. In seismic inversion, the seismic data is changed into impedance.
How seismic inversion works?
Seismic Inversion removes the imprints of the wavelet in the seismic data by deconvolution and then by converting the result into impedance.
Seismic inversion can be divided into two techniques which are post-stacking and pre-stacking inversion.
The first technique, Post-Stacking inversion is the most common approach used for inversion. This technique transforms a single seismic information volume into an acoustic impedance volume by using seismic data, well data, and basics knowledge in stratigraphy for interpretation. By removing the wavelet from seismic data, we assist at the creation of a high resolution image of the subsurface.
The post-stacking methods are colored inversion, model-based inversion, sparse spike inversion and Band-limited impedance inversion.
Colored Inversion (CI) is a process to approximately match or convolve the amplitude spectrum of a seismic with the acoustic impedance spectra issued from well data. It calculates its own inverted operator based on the frequency of seismic and well data. This method is very easy to use, fast, and strong even if there is noise. It can also be applied with or without well data and background model. The inconvenient is that the seismic input must be at zero-phase and it considers the reflectivity data of one well as the data of all the region.
In this method, the seismic trace (initial model) is convoluted with a wavelet to get a synthetic seismic trace. Then the impedance is confronted to many iterations until the difference between the inverted trace and the initial trace is reduced to a limit value.
Seismic trace = (wavelet*Reflectivity) +noise
[[File:Mbi.png|thumb|center|Model-based inversion example. It starts with an initial model and ends with a residual model. The initial model is process by iteration. Credit: Veeken 
Sparse Spike Inversion
Sparse spike inversion assumes the seismic trace models the subsurface reflectivity by using small number of acoustic impedance interface. It presents two techniques: linear programming and maximum likelihood. Note that it final purpose is to obtain high resolution impedance.
It is an algorithm which extracts an estimate of the reflectivity by using frequency domain constraints to recover the high frequencies of the seismic spectrum. This reflectivity is by the end incorporated under the initial model, where the sparse reflectivity is created. This method is to restore the impedance by minimizing error between the modeled and the initial trace.
The particularity of this algorithm is to disturb the reflectivity from the seismic trace. Assuming the wavelet is known, its job will be to add reflections coefficients until a best match is found.
Band-Limited Impedance inversion
Band-limited impedance inversion (BLI) changes post stack seismic data into impedance, density and P-wave velocity. This method is defined by the relationship between the seismic trace and seismic impedance.
To obtain multiple impedance attributes we do pre-stacking. Pre-stacking inversion changes seismic into P-impedance, S-impedance, and density by integrating well and seismic data.
The two techniques used on pre-stacking are: simultaneous and elastic inversion. Note that those technics required global wavelets and background model.
Simultaneous pre-stacking inversion is defined by :
Aki-Richards equation 
Elastic Impedance 
Elastic inversion model:
Differences and Similarities
- Farfour, M., Yoon, W. and Kim, J. (2015). Seismic attributes and acoustic impedance inversion in interpretation of complex hydrocarbon reservoirs. Journal of Applied Geophysics, 114, pp.68-80.
- Petrologic.dmt.de. (2018). Post-Stack Inversion: Petrologic Geophysical Services GmbH. [online] Available at: http://petrologic.dmt.de/services/post-stack-inversion.html [Accessed 24 Apr. 2018].
- Veeken, P. C. H., and Da Silva, M., 2004, Seismic inversion methods and some of their constraints: First Break, 22, 47-70, accessed October 28, 2017; https://www.researchgate.net/publication/277392423_Seismic_Inversion_Methods_and_some_of_their_constraints
- Maurya, S.P., and Sarkar, P., 2016, Comparison of Post stack Seismic Inversion Methods:A case study from Blackfoot Field, Canada: International Journal of Scientific & Engineering Research, 7, 1091-1101, accessed October 28, 2017; https://www.ijser.org/onlineResearchPaperViewer.aspx?Comparison-of-Post-stack-Seismic-Inversion-Methods-A-case-study-from-Blackfoot-Field,-Canada.pdf
- Petrologic.dmt.de. (2018). Pre-Stack Inversion: Petrologic Geophysical Services GmbH. [online] Available at: http://petrologic.dmt.de/services/pre-stack-inversion.html [Accessed 24 Apr. 2018].
- Doyen, P. M. (1988). Porosity from seismic data: A geostatistical approach. Geophysics, 53(10):1263-1275. https://library.seg.org/doi/pdf/10.1190/1.1442404
- Dufour, J., Squires, J., Goodway, W. N., Edmunds, A., and Shook, I. (2002). Integrated geological and geophysical interpretation case study, and lame rock parameter extractions using avo analysis on the blackfoot 3c-3d seismic data, southern alberta, canada. Geophysics, 67(1):27-37.https://library.seg.org/doi/pdf/10.1190/1.1451319
- Li, Q., 2002, Sparse Spike Inversion and the Resolution Limit: CSEG, accessed October 28, 2017; https://cseg.ca/assets/files/resources/abstracts/2002/Li_Q_Sparse_Spike_Inversion_MOD-1.pdf
- Russell, Brian H., 1988, Introduction to Seismic Inversion Methods: Society of Exploration Geophysicists, accessed October 28, 2017; https://library.seg.org/doi/book/10.1190/1.9781560802303