The Technologic Research and Application of Frequency-Divided Inversion in SL Oilfield

2014 ◽  
Vol 1010-1012 ◽  
pp. 1750-1753
Author(s):  
Peng Yu

Differing from the current inversion methods, the frequency-divided inversion is depend on well logging and seismic data, using the methods of support vector computer to study the amplitude response in different detective frequencies (amplitude versus frequency, simply called AVF), taking AVF as independent information of inversion and creating the non-linear mapping relationship between wave impedance curves of logging data and seismic waveform, fully using the information of low, medium and high frequencies in seismic data to reduce the inversion uncertainty and get a higher resolution result. Frequency-divided inversion is also a high resolution seismic inversion method which no needs wavelet extraction and prior model.

Author(s):  
Michael Gineste ◽  
Jo Eidsvik

AbstractAn ensemble-based method for seismic inversion to estimate elastic attributes is considered, namely the iterative ensemble Kalman smoother. The main focus of this work is the challenge associated with ensemble-based inversion of seismic waveform data. The amount of seismic data is large and, depending on ensemble size, it cannot be processed in a single batch. Instead a solution strategy of partitioning the data recordings in time windows and processing these sequentially is suggested. This work demonstrates how this partitioning can be done adaptively, with a focus on reliable and efficient estimation. The adaptivity relies on an analysis of the update direction used in the iterative procedure, and an interpretation of contributions from prior and likelihood to this update. The idea is that these must balance; if the prior dominates, the estimation process is inefficient while the estimation is likely to overfit and diverge if data dominates. Two approaches to meet this balance are formulated and evaluated. One is based on an interpretation of eigenvalue distributions and how this enters and affects weighting of prior and likelihood contributions. The other is based on balancing the norm magnitude of prior and likelihood vector components in the update. Only the latter is found to sufficiently regularize the data window. Although no guarantees for avoiding ensemble divergence are provided in the paper, the results of the adaptive procedure indicate that robust estimation performance can be achieved for ensemble-based inversion of seismic waveform data.


2017 ◽  
Vol 5 (4) ◽  
pp. T523-T530
Author(s):  
Ehsan Zabihi Naeini ◽  
Mark Sams

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.


Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. O21-O37 ◽  
Author(s):  
Dario Grana ◽  
Ernesto Della Rossa

A joint estimation of petrophysical properties is proposed that combines statistical rock physics and Bayesian seismic inversion. Because elastic attributes are correlated with petrophysical variables (effective porosity, clay content, and water saturation) and this physical link is associated with uncertainties, the petrophysical-properties estimation from seismic data can be seen as a Bayesian inversion problem. The purpose of this work was to develop a strategy for estimating the probability distributions of petrophysical parameters and litho-fluid classes from seismics. Estimation of reservoir properties and the associated uncertainty was performed in three steps: linearized seismic inversion to estimate the probabilities of elastic parameters, probabilistic upscaling to include the scale-changes effect, and petrophysical inversion to estimate the probabilities of petrophysical variables andlitho-fluid classes. Rock-physics equations provide the linkbetween reservoir properties and velocities, and linearized seismic modeling connects velocities and density to seismic amplitude. A full Bayesian approach was adopted to propagate uncertainty from seismics to petrophysics in an integrated framework that takes into account different sources of uncertainty: heterogeneity of the real data, approximation of physical models, measurement errors, and scale changes. The method has been tested, as a feasibility step, on real well data and synthetic seismic data to show reliable propagation of the uncertainty through the three different steps and to compare two statistical approaches: parametric and nonparametric. Application to a real reservoir study (including data from two wells and partially stacked seismic volumes) has provided as a main result the probability densities of petrophysical properties and litho-fluid classes. It demonstrated the applicability of the proposed inversion method.


2021 ◽  
pp. 1-97
Author(s):  
Lingxiao Jia ◽  
Subhashis Mallick ◽  
Cheng Wang

The choice of an initial model for seismic waveform inversion is important. In matured exploration areas with adequate well control, we can generate a suitable initial model using well information. However, in new areas where well control is sparse or unavailable, such an initial model is compromised and/or biased by the regions with more well controls. Even in matured exploration areas, if we use time-lapse seismic data to predict dynamic reservoir properties, an initial model, that we obtain from the existing preproduction wells could be incorrect. In this work, we outline a new methodology and workflow for a nonlinear prestack isotropic elastic waveform inversion. We call this method a data driven inversion, meaning that we derive the initial model entirely from the seismic data without using any well information. By assuming a locally horizonal stratification for every common midpoint and starting from the interval P-wave velocity, estimated entirely from seismic data, our method generates pseudo wells by running a two-pass one-dimensional isotropic elastic prestack waveform inversion that uses the reflectivity method for forward modeling and genetic algorithm for optimization. We then use the estimated pseudo wells to build the initial model for seismic inversion. By applying this methodology to real seismic data from two different geological settings, we demonstrate the usefulness of our method. We believe that our new method is potentially applicable for subsurface characterization in areas where well information is sparse or unavailable. Additional research is however necessary to improve the compute-efficiency of the methodology.


Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. B1-B12 ◽  
Author(s):  
Josiane Pafeng ◽  
Subhashis Mallick ◽  
Hema Sharma

Applying seismic inversion to estimate subsurface elastic earth properties for reservoir characterization is a challenge in exploration seismology. In recent years, waveform-based seismic inversions have gained popularity, but due to high computational costs, their applications are limited, and amplitude-variation-with-offset/angle inversion is still the current state-of-the-art. We have developed a genetic-algorithm-based prestack seismic waveform inversion methodology. By parallelizing at multiple levels and assuming a locally 1D structure such that forward computation of wave equation synthetics is computationally efficient, this method is capable of inverting 3D prestack seismic data on parallel computers. Applying this inversion to a real prestack seismic data volume from the Rock Springs Uplift (RSU) located in Wyoming, USA, we determined that our method is capable of inverting the data in a reasonable runtime and producing much higher quality results than amplitude-variation-with-offset/angle inversion. Because the primary purpose for seismic data acquisition at the RSU was to characterize the subsurface for potential targets for carbon dioxide sequestration, we also identified and analyzed some potential primary and secondary storage formations and their associated sealing lithologies from our inversion results.


2017 ◽  
Vol 5 (1) ◽  
pp. T1-T9 ◽  
Author(s):  
Rui Zhang ◽  
Kui Zhang ◽  
Jude E. Alekhue

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.


2021 ◽  
pp. 1-54
Author(s):  
Song Pei ◽  
Xingyao Yin ◽  
Zhaoyun Zong ◽  
Kun Li

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.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. WA185-WA200
Author(s):  
Yuqing Chen ◽  
Gerard T. Schuster

We present a wave-equation inversion method that inverts skeletonized seismic data for the subsurface velocity model. The skeletonized representation of the seismic traces consists of the low-rank latent-space variables predicted by a well-trained autoencoder neural network. The input to the autoencoder consists of seismic traces, and the implicit function theorem is used to determine the Fréchet derivative, i.e., the perturbation of the skeletonized data with respect to the velocity perturbation. The gradient is computed by migrating the shifted observed traces weighted by the skeletonized data residual, and the final velocity model is the one that best predicts the observed latent-space parameters. We denote this as inversion by Newtonian machine learning (NML) because it inverts for the model parameters by combining the forward and backward modeling of Newtonian wave propagation with the dimensional reduction capability of machine learning. Empirical results suggest that inversion by NML can sometimes mitigate the cycle-skipping problem of conventional full-waveform inversion (FWI). Numerical tests with synthetic and field data demonstrate the success of NML inversion in recovering a low-wavenumber approximation to the subsurface velocity model. The advantage of this method over other skeletonized data methods is that no manual picking of important features is required because the skeletal data are automatically selected by the autoencoder. The disadvantage is that the inverted velocity model has less resolution compared with the FWI result, but it can serve as a good initial model for FWI. Our most significant contribution is that we provide a general framework for using wave-equation inversion to invert skeletal data generated by any type of neural network. In other words, we have combined the deterministic modeling of Newtonian physics and the pattern matching capabilities of machine learning to invert seismic data by NML.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R11-R19 ◽  
Author(s):  
Ronghuo Dai ◽  
Cheng Yin ◽  
Nueraili Zaman ◽  
Fanchang Zhang

Poststack seismic impedance inversion is an effective approach for reservoir prediction. Due to the sensitivity to noise and the oscillation near the bed boundary, Gaussian distribution constrained seismic inversion is unfavorable to delineate the subtle-reservoir and small-scale geologic features. To overcome this shortcoming, we have developed a new method that incorporates a priori knowledge in the seismic inversion through a preconditioning impedance model using the adaptive edge-preserving smoothing (Ad-EPS) filter. The Ad-EPS filter preconditioned impedance model for a blocky solution makes the formation interfaces and geologic edges more precise and sharper in the inverted impedance results and keeps the inversion procedure robust even if random noise exists in the seismic data. Furthermore, compared with the conventional EPS filter, the Ad-EPS filter is able to resolve thick and thin geologic features through window size scanning, which is used to find the best-fitting window size for each sample to be filtered. The results of numerical examples and real seismic data test indicate that our inversion method can suppress noise to obtain a “blocky” inversion result and preserve small geologic features.


1981 ◽  
Vol 21 (1) ◽  
pp. 155
Author(s):  
D. B. Hays ◽  
J. Wardell

The G-LOG process is a method of seismic inversion which provides direct estimates of subsurface acoustic impedance from wavelet process stacked or migrated data. The fundamentals and characteristics of the inversion method will be discussed and examples of its use on Australian seismic data will be presented.G-LOG functions are derived by an iterative subsurface modelling technique based on a rigorous inversion of one- dimensional wave equation. This process finds the acoustic impedance model, or log, whose resulting wave-equation- consistent synthetic seismogram best matches the input seismic data in a least mean squared error sense. Multiple reflections are included in the synthetic seismogram, so that they become useful information in the determination of the log.Interval velocity logs are derived from the acoustic impedance logs. The results can be displayed in various forms, including detailed velocity logs, and colour-coded log 'sections' to match with the seismic section. Several examples of such results are presented.The G-LOG process is a revolutionary technique of subsurface modelling, and the logs it provides are strong indicators of subsurface lithology and will be an important tool in the evaluation and re-evaluation of potential hydrocarbon-bearing prospects.*Trademark of G.S.I.


Sign in / Sign up

Export Citation Format

Share Document