Seismic inversion of a carbonate buildup: A case study

2017 ◽  
Vol 5 (4) ◽  
pp. T641-T652 ◽  
Author(s):  
Mark Sams ◽  
Paul Begg ◽  
Timur Manapov
Keyword(s):  
Seismic Data ◽  
Probability Distributions ◽  
Low Frequency ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Accurate Estimation ◽  
Amplitude Analysis ◽  
Simultaneous Inversion ◽  
Band Limited ◽  
Using Data

The information within seismic data is band limited and angle limited. Together with the particular physics and geology of carbonate rocks, this imposes limitations on how accurately we can predict the presence of hydrocarbons in carbonates, map the top carbonate, and characterize the porosity distribution through seismic amplitude analysis. Using data for a carbonate reef from the Nam Con Son Basin, Vietnam, the expectations based on rock-physics analysis are that the presence of gas can be predicted only when the porosity at the top of the carbonate is extremely high ([Formula: see text]), but that a fluid contact is unlikely to be observed in the background of significant porosity variations. Mapping the top of the carbonate (except when the top carbonate porosities are low) or a fluid contact requires accurate estimates of changes in [Formula: see text]. The seismic data do not independently support such an accurate estimation of sharp changes in [Formula: see text]. The standard approach of introducing low-frequency models and applying rock-physics constraints during a simultaneous inversion does not resolve the problems: The results are heavily biased by the well control and the initial interpretation of the top carbonate and fluid contact. A facies-based inversion in which the elastic properties are restricted to values consistent with the facies predicted to be present removes the well bias, but it does not completely obviate the need for a reasonably accurate initial interpretation in terms of prior facies probability distributions. Prestack inversion improves the quality of the facies predictions compared with a poststack inversion.

Reservoir Prediction Under Control of Sedimentary Facies

2017 ◽  
Vol 25 (03) ◽  
pp. 1750022
Author(s):  
Xiuwei Yang ◽  
Peimin Zhu
Keyword(s):  
Seismic Data ◽  
Low Frequency ◽  
Prediction Method ◽  
Rock Physics ◽  
Real Data ◽  
Sedimentary Facies ◽  
Seismic Inversion ◽  
Rock Properties ◽  
Reservoir Prediction ◽  
Frequency Model

Acoustic impedance (AI) from seismic inversion can indicate rock properties and can be used, when combined with rock physics, to predict reservoir parameters, such as porosity. Solutions to seismic inversion problem are almost nonunique due to the limited bandwidth of seismic data. Additional constraints from well log data and geology are needed to arrive at a reasonable solution. In this paper, sedimentary facies is used to reduce the uncertainty in inversion and rock physics modeling; the results not only agree with seismic data, but also conform to geology. A reservoir prediction method, which incorporates seismic data, well logs, rock physics and sedimentary facies, is proposed. AI was first derived by constrained sparse spike inversion (CSSI) using a sedimentary facies dependent low-frequency model, and then was transformed to reservoir parameters by sequential simulation, statistical rock physics and [Formula: see text]-model. Two numerical experiments using synthetic model and real data indicated that the sedimentary facies information may help to obtain a more reasonable prediction.

Mapping lithofacies and pore‐fluid probabilities in a North Sea reservoir: Seismic inversions and statistical rock physics

Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 988-1001 ◽  
Author(s):  
T. Mukerji ◽  
A. Jørstad ◽  
P. Avseth ◽  
G. Mavko ◽  
J. R. Granli
Keyword(s):  
North Sea ◽  
Seismic Data ◽  
Acoustic Impedance ◽  
Probability Distributions ◽  
Rock Physics ◽  
Normal Incidence ◽  
Seismic Inversion ◽  
Training Data ◽  
Well Log ◽  
Statistical Rock Physics

Reliably predicting lithologic and saturation heterogeneities is one of the key problems in reservoir characterization. In this study, we show how statistical rock physics techniques combined with seismic information can be used to classify reservoir lithologies and pore fluids. One of the innovations was to use a seismic impedance attribute (related to the [Formula: see text] ratio) that incorporates far‐offset data, but at the same time can be practically obtained using normal incidence inversion algorithms. The methods were applied to a North Sea turbidite system. We incorporated well log measurements with calibration from core data to estimate the near‐offset and far‐offset reflectivity and impedance attributes. Multivariate probability distributions were estimated from the data to identify the attribute clusters and their separability for different facies and fluid saturations. A training data was set up using Monte Carlo simulations based on the well log—derived probability distributions. Fluid substitution by Gassmann’s equation was used to extend the training data, thus accounting for pore fluid conditions not encountered in the well. Seismic inversion of near‐offset and far‐offset stacks gave us two 3‐D cubes of impedance attributes in the interwell region. The near‐offset stack approximates a zero‐offset section, giving an estimate of the normal incidence acoustic impedance. The far offset stack gives an estimate of a [Formula: see text]‐related elastic impedance attribute that is equivalent to the acoustic impedance for non‐normal incidence. These impedance attributes obtained from seismic inversion were then used with the training probability distribution functions to predict the probability of occurrence of the different lithofacies in the interwell region. Statistical classification techniques, as well as geostatistical indicator simulations were applied on the 3‐D seismic data cube. A Markov‐Bayes technique was used to update the probabilities obtained from the seismic data by taking into account the spatial correlation as estimated from the facies indicator variograms. The final results are spatial 3‐D maps of not only the most likely facies and pore fluids, but also their occurrence probabilities. A key ingredient in this study was the exploitation of physically based seismic‐to‐reservoir property transforms optimally combined with statistical techniques.

Incorporation of geology with rock physics enables subsalt poststack inversion: A case study in the Gulf of Mexico

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB53-WB65 ◽  
Author(s):  
Huyen Bui ◽  
Jennifer Graham ◽  
Shantanu Kumar Singh ◽  
Fred Snyder ◽  
Martiris Smith
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Low Frequency ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Reservoir Properties ◽  
Blind Test ◽  
Frequency Model ◽  
Layer Sequence

One of the main goals of seismic inversion is to obtain high-resolution relative and absolute impedance for reservoir properties prediction. We aim to study whether the results from seismic inversion of subsalt data are sufficiently robust for reliable reservoir characterization. Approximately [Formula: see text] of poststack, wide-azimuth, anisotropic (vertical transverse isotropic) wave-equation migration seismic data from 50 Outer Continental Shelf blocks in the Green Canyon area of the Gulf of Mexico were inverted in this study. A total of four subsalt wells and four subsalt seismic interpreted horizons were used in the inversion process, and one of the wells was used for a blind test. Our poststack inversion method used an iterative discrete spike inversion method, based on the combination of space-adaptive wavelet processing to invert for relative acoustic impedance. Next, the dips were estimated from seismic data and converted to a horizon-like layer sequence field that was used as one of the inputs into the low-frequency model. The background model was generated by incorporating the well velocities, seismic velocity, seismic interpreted horizons, and the previously derived layer sequence field in the low-frequency model. Then, the relative acoustic impedance volume was scaled by adding the low-frequency model to match the calculated acoustic impedance logs from the wells for absolute acoustic impedance. Finally, the geological information and rock physics data were incorporated into the reservoir properties assessment for sand/shale prediction in two main target reservoirs in the Miocene and Wilcox formations. Overall, the poststack inversion results and the sand/shale prediction showed good ties at the well locations. This was clearly demonstrated in the blind test well. Hence, incorporating rock physics and geology enables poststack inversion in subsalt areas.

Integrated VTI model building with seismic data, geologic information, and rock-physics modeling — Part 1: Theory and synthetic test

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C177-C191 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols
Keyword(s):  
Seismic Data ◽  
Seismic Anisotropy ◽  
Model Building ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Additional Information ◽  
Synthetic Test ◽  
Rock Physics Modeling ◽  
Physics Modeling ◽  
Geologic Information

Seismic anisotropy plays an important role in structural imaging and lithologic interpretation. However, anisotropic model building is a challenging underdetermined inverse problem. It is well-understood that single component pressure wave seismic data recorded on the upper surface are insufficient to resolve a unique solution for velocity and anisotropy parameters. To overcome the limitations of seismic data, we have developed an integrated model building scheme based on Bayesian inference to consider seismic data, geologic information, and rock-physics knowledge simultaneously. We have performed the prestack seismic inversion using wave-equation migration velocity analysis (WEMVA) for vertical transverse isotropic (VTI) models. This image-space method enabled automatic geologic interpretation. We have integrated the geologic information as spatial model correlations, applied on each parameter individually. We integrate the rock-physics information as lithologic model correlations, bringing additional information, so that the parameters weakly constrained by seismic are updated as well as the strongly constrained parameters. The constraints provided by the additional information help the inversion converge faster, mitigate the ambiguities among the parameters, and yield VTI models that were consistent with the underlying geologic and lithologic assumptions. We have developed the theoretical framework for the proposed integrated WEMVA for VTI models and determined the added information contained in the regularization terms, especially the rock-physics constraints.

Characterization of a carbonate geothermal reservoir using rock-physics-guided deep neural networks

2021 ◽  
Vol 40 (10) ◽  
pp. 751-758
Author(s):  
Fabien Allo ◽  
Jean-Philippe Coulon ◽  
Jean-Luc Formento ◽  
Romain Reboul ◽  
Laure Capar ◽  
...  
Keyword(s):  
Neural Networks ◽  
Seismic Data ◽  
Deep Neural Networks ◽  
Synthetic Data ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Geothermal Reservoir ◽  
Rock Properties ◽  
Hybrid Techniques ◽  
Geothermal Wells

Deep neural networks (DNNs) have the potential to streamline the integration of seismic data for reservoir characterization by providing estimates of rock properties that are directly interpretable by geologists and reservoir engineers instead of elastic attributes like most standard seismic inversion methods. However, they have yet to be applied widely in the energy industry because training DNNs requires a large amount of labeled data that is rarely available. Training set augmentation, routinely used in other scientific fields such as image recognition, can address this issue and open the door to DNNs for geophysical applications. Although this approach has been explored in the past, creating realistic synthetic well and seismic data representative of the variable geology of a reservoir remains challenging. Recently introduced theory-guided techniques can help achieve this goal. A key step in these hybrid techniques is the use of theoretical rock-physics models to derive elastic pseudologs from variations of existing petrophysical logs. Rock-physics theories are already commonly relied on to generalize and extrapolate the relationship between rock and elastic properties. Therefore, they are a useful tool to generate a large catalog of alternative pseudologs representing realistic geologic variations away from the existing well locations. While not directly driven by rock physics, neural networks trained on such synthetic catalogs extract the intrinsic rock-physics relationships and are therefore capable of directly estimating rock properties from seismic amplitudes. Neural networks trained on purely synthetic data are applied to a set of 2D poststack seismic lines to characterize a geothermal reservoir located in the Dogger Formation northeast of Paris, France. The goal of the study is to determine the extent of porous and permeable layers encountered at existing geothermal wells and ultimately guide the location and design of future geothermal wells in the area.

Improving seismic QP estimation using rock-physics constraints

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. MR187-MR198 ◽  
Author(s):  
Yi Shen ◽  
Jack Dvorkin ◽  
Yunyue Li
Keyword(s):  
Seismic Data ◽  
Closed Form Solution ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Analytical Relation ◽  
Gas Reservoir ◽  
Seismic Section ◽  
Data Set ◽  
Attenuation Model ◽  
Text Model

Our goal is to accurately estimate attenuation from seismic data using model regularization in the seismic inversion workflow. One way to achieve this goal is by finding an analytical relation linking [Formula: see text] to [Formula: see text]. We derive an approximate closed-form solution relating [Formula: see text] to [Formula: see text] using rock-physics modeling. This relation is tested on well data from a clean clastic gas reservoir, of which the [Formula: see text] values are computed from the log data. Next, we create a 2D synthetic gas-reservoir section populated with [Formula: see text] and [Formula: see text] and generate respective synthetic seismograms. Now, the goal is to invert this synthetic seismic section for [Formula: see text]. If we use standard seismic inversion based solely on seismic data, the inverted attenuation model has low resolution and incorrect positioning, and it is distorted. However, adding our relation between velocity and attenuation, we obtain an attenuation model very close to the original section. This method is tested on a 2D field seismic data set from Gulf of Mexico. The resulting [Formula: see text] model matches the geologic shape of an absorption body interpreted from the seismic section. Using this [Formula: see text] model in seismic migration, we make the seismic events below the high-absorption layer clearly visible, with improved frequency content and coherency of the events.

Low frequency component of seismic data estimation and its application in seismic inversion

2020 ◽  
Author(s):  
Ding Jicai ◽  
Zhao Xiaolong ◽  
Jiang Xiudi ◽  
Wang Yandong ◽  
Huang Xiaogang ◽  
...  
Keyword(s):  
Seismic Data ◽  
Low Frequency ◽  
Frequency Component ◽  
Seismic Inversion ◽  
Data Estimation

Accuracy of wavelets, seismic inversion, and thin-bed resolution

2017 ◽  
Vol 5 (4) ◽  
pp. T523-T530
Author(s):  
Ehsan Zabihi Naeini ◽  
Mark Sams
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Low Frequency ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Frequency Content ◽  
Frequency Model ◽  
North West ◽  
The North ◽  
Well Data

Broadband reprocessed seismic data from the North West Shelf of Australia were inverted using wavelets estimated with a conventional approach. The inversion method applied was a facies-based inversion, in which the low-frequency model is a product of the inversion process itself, constrained by facies-dependent input trends, the resultant facies distribution, and the match to the seismic. The results identified the presence of a gas reservoir that had recently been confirmed through drilling. The reservoir is thin, with up to 15 ms of maximum thickness. The bandwidth of the seismic data is approximately 5–70 Hz, and the well data used to extract the wavelet used in the inversion are only 400 ms long. As such, there was little control on the lowest frequencies of the wavelet. Different wavelets were subsequently estimated using a variety of new techniques that attempt to address the limitations of short well-log segments and low-frequency seismic. The revised inversion showed greater gas-sand continuity and an extension of the reservoir at one flank. Noise-free synthetic examples indicate that thin-bed delineation can depend on the accuracy of the low-frequency content of the wavelets used for inversion. Underestimation of the low-frequency contents can result in missing thin beds, whereas underestimation of high frequencies can introduce false thin beds. Therefore, it is very important to correctly capture the full frequency content of the seismic data in terms of the amplitude and phase spectra of the estimated wavelets, which subsequently leads to a more accurate thin-bed reservoir characterization through inversion.

scholarly journals Low-frequency seismic: the next revolution in resolution

2014 ◽  
Vol 54 (1) ◽  
pp. 69
Author(s):  
Andrew Long ◽  
Cyrille Reiser
Keyword(s):  
New Technologies ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Seismic Inversion ◽  
Reservoir Properties ◽  
Well Control ◽  
Velocity Models ◽  
The Earth ◽  
Air Gun ◽  
Band Limited

Ultra-low seismic frequencies less than about 7 Hz cannot be produced by conventional air gun arrays, for any configuration and for any towing depth. There is a profound difference between improving low-frequency recovery by removing source and receiver ghosts (achievable) and improving low-frequency injection on the source side (an unrealised dream). If 1–7 Hz amplitudes could be usefully injected into the earth, it would be possible to facilitate much sharper seismic representation of geological contacts and internal features, and seismic inversion would yield robust and precise predictions of reservoir properties—without well control. The net result is fewer exploration and appraisal wells, greatly reduced exploration and development risks, and optimised recoverable reserves. Furthermore, an emerging seismic pursuit known as full waveform inversion (FWI) makes the bold promise that raw seismic field gathers can be directly used to invert for the highest achievable velocity models, almost without any human intervention. These models will bypass the traditional lack of low-frequency information in band-limited seismic data, and facilitate the aforementioned ambition of seismic inversion without well control. FWI, however, is confronted by the paradox that ultra-low-frequency seismic gathers are the necessary input for stable results. This paper describes new technologies that may enable the injection of strong 2–7 Hz amplitudes into the earth, and explains in simple terms how FWI can already be pursued as a robust complement to the prediction of accurate reservoir properties. The low-frequency revolution is already here.

Simultaneous inversion of prestack seismic data for rock properties using simulated annealing

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1877-1885 ◽  
Author(s):  
Xin‐Quan Ma
Keyword(s):  
Simulated Annealing ◽  
Seismic Data ◽  
Low Frequency ◽  
Optimization Procedure ◽  
P Wave ◽  
Rock Properties ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
Simultaneous Inversion ◽  
Reflection Seismic

A new prestack inversion algorithm has been developed to simultaneously estimate acoustic and shear impedances from P‐wave reflection seismic data. The algorithm uses a global optimization procedure in the form of simulated annealing. The goal of optimization is to find a global minimum of the objective function, which includes the misfit between synthetic and observed prestack seismic data. During the iterative inversion process, the acoustic and shear impedance models are randomly perturbed, and the synthetic seismic data are calculated and compared with the observed seismic data. To increase stability, constraints have been built into the inversion algorithm, using the low‐frequency impedance and background Vs/Vp models. The inversion method has been successfully applied to synthetic and field data examples to produce acoustic and shear impedances comparable to log data of similar bandwidth. The estimated acoustic and shear impedances can be combined to derive other elastic parameters, which may be used for identifying of lithology and fluid content of reservoirs.

Sign in / Sign up

Export Citation Format

Share Document