In geophysics (primarily in oil-and-gas exploration/development), seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion may be pre- or post-stack, deterministic, random or geostatistical; it typically includes other reservoir measurements such as well logs and cores.[1]
Geophysicists routinely perform seismic surveys 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 side-lobe 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:
pre-stack or post-stack
seismic resolution or well-log 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 advances 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 field data and well data, where well data serves to add the high 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 band-limited, reducing resolution and quality. To extend the frequency band available, low-frequency data is derived from log data, pre-stack 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 limited-trend 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).
Post-stack seismic resolution inversion
An example of a post-stack seismic resolution inversion technique is the constrained sparse-spike 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 p-wave 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 pseudo-acoustic 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 high-quality acoustic impedance volumes from full or post-stack seismic data: full-bandwidth impedance, bandlimited Impedance, reflectivity model, and low-frequency 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 mixed-norm approach and allow a weighting factor between minimizing the sparsity of the solution and minimizing the misfit of the residual traces.
Pre-stack seismic resolution inversion
Pre-stack inversion is often used when post-stack inversion fails to sufficiently differentiate geologic features with similar P-impedance signatures.[6] Simultaneous inversion solves for S-impedance and density, in addition to P-impedance. While many geologic features can express similar P-impedance characteristics, few will share combined P-impedance and S-impedance 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 P-impedance alone or whether S-impedance is also required. This will dictate whether a pre- or post-stack inversion is needed.
Simultaneous Inversion (SI) is a pre-stack method that uses multiple offset or angle seismic sub-stacks and their associated wavelets as input; it generates P-impedance, S-impedance 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 offset-dependent phase, bandwidth, tuning and NMO stretch effects.[7]
The inversion algorithm works by first estimating angle-dependent P-wave reflectivities for the input-partial stacks. Next, these are used with the full Zoeppritz equations (or approximations, such as Aki–Richards, for some algorithms) to find band-limited elastic reflectivities. These are in turn merged with their low-frequency counterparts from the model and integrated to elastic properties. This approximate result is then improved in a final inversion for P-impedance, S-impedance 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 angle-dependent well log reflection coefficients in the region of interest to the corresponding offset stack at the well locations. Reflection coefficients are calculated from P-sonic, S-sonic 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 offset-dependent bandwidth, scaling and tuning effects. A near-stack wavelet can be used as the starting point for estimating the far-angle (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 angle-stacked seismic data sets and generates three elastic parameter volumes as output.
The resulting elastic parameters are real-rock properties that can be directly related to reservoir properties. The more advanced algorithms use the full Knott–Zoeppritz equations and there is full allowance for amplitude and phase variations with offset. This is done by deriving unique wavelets for each input-partial stack. The elastic parameters themselves can be directly constrained during the seismic inversion and rock-physics relationships can be applied, constraining pairs of elastic parameters to each other. Final elastic-parameter models optimally reproduce the input seismic, as this is part of the seismic inversion optimization.
Post stack geostatistical inversion
Geostatistical inversion integrates high resolution well data with low resolution 3-D seismic, and provides a model with high vertical detail near and away from well control. This generates reservoir models with geologically-plausible shapes, and provides a clear quantification of uncertainty to assess risk. Highly detailed petrophysical models are generated, ready for input to reservoir-flow 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 one-step 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 co-simulate 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.
Pre-stack log-detail 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 (P-impedance, S-impedance 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 leading-edge geostatistical techniques, including Markov chain Monte Carlo (MCMC) sampling and pluri-Gaussian 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 multiple-angle-stack 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 rock-physics 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 p-sonic, s-sonic and density logs tied to the appropriate events on the seismic data. Wavelets are extracted individually for each well. A final "multi-well" 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 3-D bodies of hydrocarbon-bearing units are captured with their corresponding rock properties, and the uncertainty in reservoir size and properties is quantified.
In petroleum exploration and development, formation evaluation is used to determine the ability of a borehole to produce petroleum. Essentially, it is the process of "recognizing a commercial well when you drill one".
In the oil and gas industry, the term wireline usually refers to the use of multi-conductor, single conductor or slickline cable, or "wireline", as a conveyance for the acquisition of subsurface petrophysical and geophysical data and the delivery of well construction services such as pipe recovery, perforating, plug setting and well cleaning and fishing. The subsurface geophysical and petrophysical information results in the description and analysis of subsurface geology, reservoir properties and production characteristics.
Reflection seismology is a method of exploration geophysics that uses the principles of seismology to estimate the properties of the Earth's subsurface from reflected seismic waves. The method requires a controlled seismic source of energy, such as dynamite or Tovex blast, a specialized air gun or a seismic vibrator. Reflection seismology is similar to sonar and echolocation.
Well logging, also known as borehole logging is the practice of making a detailed record of the geologic formations penetrated by a borehole. The log may be based either on visual inspection of samples brought to the surface or on physical measurements made by instruments lowered into the hole. Some types of geophysical well logs can be done during any phase of a well's history: drilling, completing, producing, or abandoning. Well logging is performed in boreholes drilled for the oil and gas, groundwater, mineral and geothermal exploration, as well as part of environmental and geotechnical studies.
In geophysics and reflection seismology, amplitude versus offset (AVO) or amplitude variation with offset is the general term for referring to the dependency of the seismic attribute, amplitude, with the distance between the source and receiver. AVO analysis is a technique that geophysicists can execute on seismic data to determine a rock's fluid content, porosity, density or seismic velocity, shear wave information, fluid indicators.
Exploration geophysics is an applied branch of geophysics and economic geology, which uses physical methods at the surface of the Earth, such as seismic, gravitational, magnetic, electrical and electromagnetic, to measure the physical properties of the subsurface, along with the anomalies in those properties. It is most often used to detect or infer the presence and position of economically useful geological deposits, such as ore minerals; fossil fuels and other hydrocarbons; geothermal reservoirs; and groundwater reservoirs. It can also be used to detect the presence of unexploded ordnance.
Seismic anisotropy is the directional dependence of the velocity of seismic waves in a medium (rock) within the Earth.
Petrophysics is the study of physical and chemical rock properties and their interactions with fluids.
Seismic migration is the process by which seismic events are geometrically re-located in either space or time to the location the event occurred in the subsurface rather than the location that it was recorded at the surface, thereby creating a more accurate image of the subsurface. This process is necessary to overcome the limitations of geophysical methods imposed by areas of complex geology, such as: faults, salt bodies, folding, etc.
A Stoneley wave is a boundary wave that typically propagates along a solid-solid interface. When found at a liquid-solid interface, this wave is also referred to as a Scholte wave. The wave is of maximum intensity at the interface and decreases exponentially away from it. It is named after the British seismologist Dr. Robert Stoneley (1894–1976), a lecturer in the University of Leeds, who discovered it on October 1, 1924.
Heavy oil production is a developing technology for extracting heavy oil in industrial quantities. Estimated reserves of heavy oil are over 6 trillion barrels, three times that of conventional oil and gas.
Michael Schoenberg (1939–2008) was an American theoretical geophysicist noted for his fundamental contributions to the understanding of anisotropy in the real earth and its application to the determination of texture, fracture porosity, and flow properties of reservoir rocks.
A synthetic seismogram is the result of forward modelling the seismic response of an input earth model, which is defined in terms of 1D, 2D or 3D variations in physical properties. In hydrocarbon exploration this is used to provide a 'tie' between changes in rock properties in a borehole and seismic reflection data at the same location. It can also be used either to test possible interpretation models for 2D and 3D seismic data or to model the response of the predicted geology as an aid to planning a seismic reflection survey. In the processing of wide-angle reflection and refraction (WARR) data, synthetic seismograms are used to further constrain the results of seismic tomography. In earthquake seismology, synthetic seismograms are used either to match the predicted effects of a particular earthquake source fault model with observed seismometer records or to help constrain the Earth's velocity structure. Synthetic seismograms are generated using specialized geophysical software.
In the oil and gas industry, reservoir modeling involves the construction of a computer model of a petroleum reservoir, for the purposes of improving estimation of reserves and making decisions regarding the development of the field, predicting future production, placing additional wells and evaluating alternative reservoir management scenarios.
In reflection seismology, a seismic attribute is a quantity extracted or derived from seismic data that can be analysed in order to enhance information that might be more subtle in a traditional seismic image, leading to a better geological or geophysical interpretation of the data. Examples of seismic attributes can include measured time, amplitude, frequency and attenuation, in addition to combinations of these. Most seismic attributes are post-stack, but those that use CMP gathers, such as amplitude versus offset (AVO), must be analysed pre-stack. They can be measured along a single seismic trace or across multiple traces within a defined window.
Seismic inversion involves the set of methods which seismologists use to infer properties through physical measurements. Surface-wave inversion is the method by which elastic properties, density, and thickness of layers in the subsurface are obtained through analysis of surface-wave dispersion. The entire inversion process requires the gathering of seismic data, the creation of dispersion curves, and finally the inference of subsurface properties.
In reflection seismology, the anelastic attenuation factor, often expressed as seismic quality factor or Q, quantifies the effects of anelastic attenuation on the seismic wavelet caused by fluid movement and grain boundary friction. As a seismic wave propagates through a medium, the elastic energy associated with the wave is gradually absorbed by the medium, eventually ending up as heat energy. This is known as absorption and will eventually cause the total disappearance of the seismic wave.
Inverse modeling is a mathematical technique where the objective is to determine the physical properties of the subsurface of an earth region that has produced a given seismogram. Cooke and Schneider (1983) defined it as calculation of the earth's structure and physical parameters from some set of observed seismic data. The underlying assumption in this method is that the collected seismic data are from an earth structure that matches the cross-section computed from the inversion algorithm. Some common earth properties that are inverted for include acoustic velocity, formation and fluid densities, acoustic impedance, Poisson's ratio, formation compressibility, shear rigidity, porosity, and fluid saturation.
Seismic data acquisition is the first of the three distinct stages of seismic exploration, the other two being seismic data processing and seismic interpretation. Seismic acquisition requires the use of a seismic source at specified locations for a seismic survey, and the energy that travels within the subsurface as seismic waves generated by the source gets recorded at specified locations on the surface by what is known as receivers.
The Marmousi model is a complex 2D structural model and its seismic response devised by the Institut Français du Petrole, with strong horizontal and vertical velocity changes. Since its inception in 1988, Marmousi has become a standard to compare depth-migration and velocity determination models and to this day remains one of the most published geophysical data sets.
↑ Pendrel, J., "Seismic Inversion—A Critical Tool in Reservoir Characterization", Scandinavian Oil-Gas Magazine, No. 5/6, 2006, pp. 19–22 .
1 2 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 Seismic-Derived Acoustic Impedance Data", The Leading Edge, March 2000, pp. 242–256.
↑ 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, pp.213–214.
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., Torres-Verdin, C., Chesters, W., Kvien, K., Globe, M., "Joint Stochastic Inversion of Petrophysical Logs and 3D Pre-Stack Seismic Data to Assess the Spatial Continuity of Fluid Units Away from Wells: Application to a Gulf-of-Mexico 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 Poro-elastic 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, pp.5–11.
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, Thirty-First Annual Convention and Exhibition, May 2007.
Russell, B., Hampson, D., "The Old and the New in Seismic Inversion", CSEG Recorder, December 2006, pp.5–11.
Stephen, K., MacBeth, C., "Reducing Reservoir Prediction Uncertainty by Updating a Stochastic Model Using Seismic History Matching", SPE Reservoir Evaluation & Engineering, December 2008.
Vargas-Meleza, L., Megchun, J., Vazquez, G., "Petrophysical Properties Estimation by Integrating AVO, Seismic Inversion and Multiattribute Analysis in a 3-D 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, 799–813, December 2006.
Watson, I., Lines, L., "Seismic Inversion at Pike’s Peak, Saskatchewan", CREWES Research Report, Volume 12, 2000.
Whitfield, J., "The Relation of Net Pay to Amplitude Versus Offset Gradients: A Gulf of Mexico Case Study", University of Houston Master's Thesis, 1993.
Zou, Y., Bentley, L., Lines, L., "Integration of Reservoir Simulation with Time-Lapse Seismic Modeling", 2004 CSEG National Convention.
This page is based on this Wikipedia article Text is available under the CC BY-SA 4.0 license; additional terms may apply. Images, videos and audio are available under their respective licenses.
