Depth-domain seismic reflectivity inversion with compressed sensing technique

2017 ◽  
Vol 5 (1) ◽  
pp. T1-T9 ◽  
Author(s):  
Rui Zhang ◽  
Kui Zhang ◽  
Jude E. Alekhue
Keyword(s):  
Compressed Sensing ◽  
Time Domain ◽  
Seismic Data ◽  
Reservoir Characterization ◽  
Synthetic Data ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Depth Migration ◽  
Band Limited ◽  
The Time Domain

More and more seismic surveys produce 3D seismic images in the depth domain by using prestack depth migration methods, which can present a direct subsurface structure in the depth domain rather than in the time domain. This leads to the increasing need for applications of seismic inversion on the depth-imaged seismic data for reservoir characterization. To address this issue, we have developed a depth-domain seismic inversion method by using the compressed sensing technique with output of reflectivity and band-limited impedance without conversion to the time domain. The formulations of the seismic inversion in the depth domain are similar to time-domain methods, but they implement all the elements in depth domain, for example, a depth-domain seismic well tie. The developed method was first tested on synthetic data, showing great improvement of the resolution on inverted reflectivity. We later applied the method on a depth-migrated field data with well-log data validated, showing a great fit between them and also improved resolution on the inversion results, which demonstrates the feasibility and reliability of the proposed method on depth-domain seismic data.

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.

High resolution fixed-point seismic inversion

2021 ◽  
pp. 1-54
Author(s):  
Song Pei ◽  
Xingyao Yin ◽  
Zhaoyun Zong ◽  
Kun Li
Keyword(s):  
Fixed Point ◽  
Objective Function ◽  
Seismic Data ◽  
Target Function ◽  
Recursive Formula ◽  
Frequency Component ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Good Convergence ◽  
Band Limited

Resolution improvement always presents the crucial task in geological inversion. Band-limited characteristics of seismic data and noise make seismic inversion complicated. Specifically, geological inversion suffers from the deficiency of both low- and high-frequency components. We propose the fixed-point seismic inversion method to alleviate these issues. The problem of solving objective function is transformed into the problem of finding the fixed-point of objective function. Concretely, a recursive formula between seismic signal and reflection coefficient is established, which is characterized by good convergence and verified by model examples. The error between the model value and the inverted value is reduced to around zero after few iterations. The model examples show that in either case, that is, the seismic traces are noise-free or with a little noise, the model value can almost be duplicated. Even if the seismic trace is accompanied by the moderate noise, the optimal inverted results can still be obtained with the proposed method. The initial model constraint is further introduced into the objective function to increase the low-frequency component of the inverted results by adding prior information into the target function. The singular value decomposition (SVD) method is applied to the inversion framework, thus making a high improvement of anti-noise ability. At last, the synthetic models and seismic data are investigated following the proposed method. The inverted results obtained from the fixed-point seismic inversion are compared with those obtained from the conventional seismic inversion, and it is found that the former has a higher resolution than the latter.

Time-domain inversion of GPR data containing attenuation due to conductive losses

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. K103-K109 ◽  
Author(s):  
Qingyun Di ◽  
Meigen Zhang ◽  
Maioyue Wang
Keyword(s):  
Time Domain ◽  
Ground Penetrating Radar ◽  
Seismic Data ◽  
Random Noise ◽  
Synthetic Data ◽  
The Time Domain ◽  
Ground Penetrating ◽  
Domain Inversion ◽  
Inversion Techniques ◽  
And Inversion

Many seismic data processing and inversion techniques have been applied to ground-penetrating radar (GPR) data without including the wave field attenuation caused by conductive ground. Neglecting this attenuation often reduces inversion resolution. This paper introduces a GPR inversion technique that accounts for the effects of attenuation. The inversion is formulated in the time domain with the synthetic GPR waveforms calculated by a finite-element method (FEM). The Jacobian matrix can be computed efficiently with the same FEM forward modeling procedure. Synthetic data tests show that the inversion can generate high-resolution subsurface velocity profiles even with data containing strong random noise. The inversion can resolve small objects not readily visible in the waveforms. Further, the inversion yields a dielectric constant that can help to determine the types of material filling underground cavities.

A depth variant seismic wavelet extraction method for inversion of poststack depth-domain seismic data

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. R569-R579 ◽  
Author(s):  
Rui Zhang ◽  
Zhiwen Deng
Keyword(s):  
Seismic Data ◽  
Extraction Method ◽  
Reservoir Characterization ◽  
Real Data ◽  
Velocity Model ◽  
Seismic Inversion ◽  
Data Depth ◽  
Spectral Variation ◽  
Seismic Wavelet ◽  
The Time Domain

Prestack depth seismic imaging is increasingly being used in industry, which has also led to an increasing need for its inversion results, such as acoustic impedance (AI), for reservoir characterization. Conventional seismic inversion methods for reservoir characterization are usually implemented in the time domain. A depth-time conversion would be required before inversion of depth-domain seismic data, which would depend on an accurate velocity model and a fine time-depth conversion algorithm. Thus, it could be beneficial that we can directly invert the depth migrated seismic data. Depth-domain seismic data could indicate a strong nonstationarity, such as spectral variation, which makes it difficult to use a constant wavelet for direct inversion in depth. To address this issue, we have developed a new wavelet extraction method by using a depth-wavenumber decomposition technique, which can generate depth variant wavelets to accommodate the nonstationarity of the depth-domain seismic data. The synthetic and real data applications have been used to test the effectiveness of our method. The directly inverted depth-domain AI indicates a good correlation with well-log data and a strong potential for reservoir characterization.

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.

Applications of quantitative prestack seismic analysis to unconventional resource play characterization in the Permian/Delaware Basin

2019 ◽  
Vol 38 (2) ◽  
pp. 106-115 ◽  
Author(s):  
Phuong Hoang ◽  
Arcangelo Sena ◽  
Benjamin Lascaud
Keyword(s):  
Hydraulic Fracturing ◽  
Seismic Data ◽  
Reservoir Characterization ◽  
Seismic Analysis ◽  
Seismic Inversion ◽  
Reservoir Properties ◽  
Yield Attributes ◽  
Horizontal Drilling ◽  
Delaware Basin ◽  
Unconventional Resource

The characterization of shale plays involves an understanding of tectonic history, geologic settings, reservoir properties, and the in-situ stresses of the potential producing zones in the subsurface. The associated hydrocarbons are generally recovered by horizontal drilling and hydraulic fracturing. Historically, seismic data have been used mainly for structural interpretation of the shale reservoirs. A primary benefit of surface seismic has been the ability to locate and avoid drilling into shallow carbonate karsting zones, salt structures, and basement-related major faults which adversely affect the ability to drill and complete the well effectively. More recent advances in prestack seismic data analysis yield attributes that appear to be correlated to formation lithology, rock strength, and stress fields. From these, we may infer preferential drilling locations or sweet spots. Knowledge and proper utilization of these attributes may prove valuable in the optimization of drilling and completion activities. In recent years, geophysical data have played an increasing role in supporting well planning, hydraulic fracturing, well stacking, and spacing. We have implemented an integrated workflow combining prestack seismic inversion and multiattribute analysis, microseismic data, well-log data, and geologic modeling to demonstrate key applications of quantitative seismic analysis utilized in developing ConocoPhillips' acreage in the Delaware Basin located in Texas. These applications range from reservoir characterization to well planning/execution, stacking/spacing optimization, and saltwater disposal. We show that multidisciplinary technology integration is the key for success in unconventional play exploration and development.

Case study: Seismic resolution and reservoir characterization of thin sands using multiattribute analysis and bandwidth extension in the Daqing field, China

2020 ◽  
Vol 8 (1) ◽  
pp. T89-T102
Author(s):  
David Mora ◽  
John Castagna ◽  
Ramses Meza ◽  
Shumin Chen ◽  
Renqi Jiang
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Statistical Significance ◽  
Physical Property ◽  
Seismic Attributes ◽  
Oil Field ◽  
Thin Layers ◽  
Inversion Method ◽  
Bandwidth Extension ◽  
Multiattribute Analysis

The Daqing field, located in the Songliao Basin in northeastern China, is the largest oil field in China. Most production in the Daqing field comes from seismically thin sand bodies with thicknesses between 1 and 15 m. Thus, it is not usually possible to resolve Daqing reservoirs using only conventional seismic data. We have evaluated the effectiveness of seismic multiattribute analysis of bandwidth extended data in resolving and making inferences about these thin layers. Multiattribute analysis uses statistical methods or neural networks to find relationships between well data and seismic attributes to predict some physical property of the earth. This multiattribute analysis was applied separately to conventional seismic data and seismic data that were spectrally broadened using sparse-layer inversion because this inversion method usually increases the vertical resolution of the seismic. Porosity volumes were generated using target porosity logs and conventional seismic attributes, and isofrequency volumes were obtained by spectral decomposition. The resulting resolution, statistical significance, and accuracy in the determination of layer properties were higher for the predictions made using the spectrally broadened volume.

Estimating interannual variability arising from weather events

2001 ◽  
Vol 38 (A) ◽  
pp. 274-288 ◽  
Author(s):  
Xiaogu Zheng ◽  
James Renwick
Keyword(s):  
Frequency Domain ◽  
Interannual Variability ◽  
Time Domain ◽  
Synthetic Data ◽  
Estimation Procedure ◽  
Frequency Domain Method ◽  
Domain Method ◽  
Weather Events ◽  
Long Time ◽  
The Time Domain

The advantages and limitations of frequency domain and time domain methods for estimating the interannual variability arising from day-to-day weather events are summarized. A modification of the time domain method is developed and its application in examining a precondition for the frequency domain method is demonstrated. A combined estimation procedure is proposed: it takes advantage of the strengths of both methods. The estimation procedures are tested with sets of synthetic data and are applied to long time series of three meteorological parameters. The impacts of the different methods on tests of potential long-range predictability for seasonal means are also discussed.

Colored and linear inversions to relative acoustic impedance

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. N15-N27 ◽  
Author(s):  
Carlos A. M. Assis ◽  
Henrique B. Santos ◽  
Jörg Schleicher
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Computational Cost ◽  
Synthetic Data ◽  
Seismic Source ◽  
Alternative Methods ◽  
Well Log ◽  
Linear Inversion ◽  
Amplitude Inversion ◽  
Band Limited

Acoustic impedance (AI) is a widely used seismic attribute in stratigraphic interpretation. Because of the frequency-band-limited nature of seismic data, seismic amplitude inversion cannot determine AI itself, but it can only provide an estimate of its variations, the relative AI (RAI). We have revisited and compared two alternative methods to transform stacked seismic data into RAI. One is colored inversion (CI), which requires well-log information, and the other is linear inversion (LI), which requires knowledge of the seismic source wavelet. We start by formulating the two approaches in a theoretically comparable manner. This allows us to conclude that both procedures are theoretically equivalent. We proceed to check whether the use of the CI results as the initial solution for LI can improve the RAI estimation. In our experiments, combining CI and LI cannot provide superior RAI results to those produced by each approach applied individually. Then, we analyze the LI performance with two distinct solvers for the associated linear system. Moreover, we investigate the sensitivity of both methods regarding the frequency content present in synthetic data. The numerical tests using the Marmousi2 model demonstrate that the CI and LI techniques can provide an RAI estimate of similar accuracy. A field-data example confirms the analysis using synthetic-data experiments. Our investigations confirm the theoretical and practical similarities of CI and LI regardless of the numerical strategy used in LI. An important result of our tests is that an increase in the low-frequency gap in the data leads to slightly deteriorated CI quality. In this case, LI required more iterations for the conjugate-gradient least-squares solver, but the final results were not much affected. Both methodologies provided interesting RAI profiles compared with well-log data, at low computational cost and with a simple parameterization.

Sign in / Sign up

Export Citation Format

Share Document