A data-driven amplitude variation with offset inversion method via learned dictionaries and sparse representation

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. R725-R748 ◽  
Author(s):  
Bin She ◽  
Yaojun Wang ◽  
Jiandong Liang ◽  
Zhining Liu ◽  
Chengyun Song ◽  
...  

Amplitude variation with offset (AVO) inversion is a typical ill-posed inverse problem. To obtain a stable and unique solution, regularization techniques relying on mathematical models from prior information are commonly used in conventional AVO inversion methods (hence the name model-driven methods). Due to the difference between prior information and the actual geology, these methods often have difficulty achieving satisfactory accuracy and resolution. We have developed a novel data-driven inversion method for the AVO inversion problem. This method can effectively extract useful knowledge from well-log data, including sparse dictionaries of elastic parameters and sparse representation of subsurface model parameters. Lateral continuity of subsurface geology allows for the approximation of model parameters for a work area using the learned dictionaries. Instead of particular mathematical models, a sparse representation is used to constrain the inverse problem. Because no assumption is made about the model parameters, we consider this a data-driven method. The general process of the algorithm is as follows: (1) using well-log data as the training samples to learn the sparse dictionary of each elastic parameter, (2) imposing a sparse representation constraint on the objective function, making the elastic parameters be sparsely represented over the learned dictionary, and (3) solving the objective function by applying a coordinate-descent algorithm. Tests on several synthetic examples and field data demonstrate that our algorithm is effective in improving the resolution and accuracy of solutions and is adaptable to various geologies.

2016 ◽  
Vol 4 (4) ◽  
pp. T613-T625 ◽  
Author(s):  
Qizhen Du ◽  
Bo Zhang ◽  
Xianjun Meng ◽  
Chengfeng Guo ◽  
Gang Chen ◽  
...  

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reflectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reflectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.


Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. R245-R255 ◽  
Author(s):  
Ali Gholami ◽  
Hossein S. Aghamiry ◽  
Mostafa Abbasi

The inversion of prestack seismic data using amplitude variation with offset (AVO) has received increased attention in the past few decades because of its key role in estimating reservoir properties. AVO is mainly governed by the Zoeppritz equations, but traditional inversion techniques are based on various linear or quasilinear approximations to these nonlinear equations. We have developed an efficient algorithm for nonlinear AVO inversion of precritical reflections using the exact Zoeppritz equations in multichannel and multi-interface form for simultaneous estimation of the P-wave velocity, S-wave velocity, and density. The total variation constraint is used to overcome the ill-posedness while solving the forward nonlinear model and to preserve the sharpness of the interfaces in the parameter space. The optimization is based on a combination of Levenberg’s algorithm and the split Bregman iterative scheme, in which we have to refine the data and model parameters at each iteration. We refine the data via the original nonlinear equations, but we use the traditional cost-effective linearized AVO inversion to construct the Jacobian matrix and update the model. Numerical experiments show that this new iterative procedure is convergent and converges to a solution of the nonlinear problem. We determine the performance and optimality of our nonlinear inversion algorithm with various simulated and field seismic data sets.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. R669-R679 ◽  
Author(s):  
Gang Chen ◽  
Xiaojun Wang ◽  
Baocheng Wu ◽  
Hongyan Qi ◽  
Muming Xia

Estimating the fluid property factor and density from amplitude-variation-with-offset (AVO) inversion is important for fluid identification and reservoir characterization. The fluid property factor can distinguish pore fluid in the reservoir and the density estimate aids in evaluating reservoir characteristics. However, if the scaling factor of the fluid property factor (the dry-rock [Formula: see text] ratio) is chosen inappropriately, the fluid property factor is not only related to the pore fluid, but it also contains a contribution from the rock skeleton. On the other hand, even if the angle gathers include large angles (offsets), a three-parameter AVO inversion struggles to estimate an accurate density term without additional constraints. Thus, we have developed an equation to compute the dry-rock [Formula: see text] ratio using only the P- and S-wave velocities and density of the saturated rock from well-logging data. This decouples the fluid property factor from lithology. We also developed a new inversion method to estimate the fluid property factor and density parameters, which takes full advantage of the high stability of a two-parameter AVO inversion. By testing on a portion of the Marmousi 2 model, we find that the fluid property factor calculated by the dry-rock [Formula: see text] ratio obtained by our method relates to the pore-fluid property. Simultaneously, we test the AVO inversion method for estimating the fluid property factor and density parameters on synthetic data and analyze the feasibility and stability of the inversion. A field-data example indicates that the fluid property factor obtained by our method distinguishes the oil-charged sand channels and the water-wet sand channel from the well logs.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. N31-N50 ◽  
Author(s):  
Jun Lu ◽  
Yun Wang ◽  
Jingyi Chen ◽  
Ying An

With the increase in exploration target complexity, more parameters are required to describe subsurface properties, particularly for finely stratified reservoirs with vertical transverse isotropic (VTI) features. We have developed an anisotropic amplitude variation with offset (AVO) inversion method using joint PP and PS seismic data for VTI media. Dealing with local minimum solutions is critical when using anisotropic AVO inversion because more parameters are expected to be derived. To enhance the inversion results, we adopt a hierarchical inversion strategy to solve the local minimum solution problem in the Gauss-Newton method. We perform the isotropic and anisotropic AVO inversions in two stages; however, we only use the inversion results from the first stage to form search windows for constraining the inversion in the second stage. To improve the efficiency of our method, we built stop conditions using Euclidean distance similarities to control iteration of the anisotropic AVO inversion in noisy situations. In addition, we evaluate a time-aligned amplitude variation with angle gather generation approach for our anisotropic AVO inversion using anisotropic prestack time migration. We test the proposed method on synthetic data in ideal and noisy situations, and find that the anisotropic AVO inversion method yields reasonable inversion results. Moreover, we apply our method to field data to show that it can be used to successfully identify complex lithologic and fluid information regarding fine layers in reservoirs.


2020 ◽  
Vol 8 (1) ◽  
pp. SA63-SA72
Author(s):  
Wu Haibo ◽  
Cheng Yan ◽  
Zhang Pingsong ◽  
Dong Shouhua ◽  
Huang Yaping

The brittleness index (BI) is an important parameter for coal-bed methane (CBM) reservoir fracturing characterization. Most published studies have relied on petrophysical and well-log data to estimate the geomechanical properties of reservoir rocks. The major drawback of such methods is the lack of control away from well locations. Therefore, we have developed a method of combining BI calculation from well logs with that inverted from 3D seismic data to overcome the limitation. A real example is given here to indicate the workflow. A traditional amplitude-variation-with-offset (AVO) inversion was conducted first. BI for the CBM reservoir was then calculated from the Lamé constants inverted from prestack seismic data through a traditional AVO inversion method. We build an initial low-frequency model based on the well-log data. Comparison of the seismic inverted BI at the target reservoir and BI extracted from the well-log data showed satisfactory results. This method has been proved to be efficient and effective enough at identifying BI sweet spots in CBM reservoirs.


2017 ◽  
Vol 5 (3) ◽  
pp. SL57-SL67 ◽  
Author(s):  
Guangsen Cheng ◽  
Xingyao Yin ◽  
Zhaoyun Zong

Prestack seismic inversion is widely used in fluid indication and reservoir prediction. Compared with linear inversion, nonlinear inversion is more precise and can be applied to high-contrast situations. The inversion results can be affected by the parameters’ sensitivity, so the parameterization of nonlinear equations is very significant. Considering the poor nonlinear amplitude-variation-with-offset (AVO) inversion results of impedance and velocity parameters, we adjust the parameters of the nonlinear equation, avoid the inaccuracy caused by parameters sensitivity and get the ideal nonlinear AVO inversion results of the Lamé parameters. The feasibility and stability of the nonlinear equation based on the Lamé parameters and method are verified by the model and the real data examples. The resolution and the lateral continuity of nonlinear inversion results are better compared with the linear inversion results.


Author(s):  
Liping Niu ◽  
Jianhua Geng ◽  
Xinming Wu ◽  
Luanxiao Zhao ◽  
Hong Zhang

Abstract Linearised approximations of the Zoeppritz equation are widely used for amplitude versus offset (AVO) forward modelling and prestack seismic inversion. However, current linearised approximations typically exhibit large errors for strong-contrast media and large incidence angles (even for those less than the critical angle), leading to poor quality estimation of elastic parameters. We therefore develop a data-driven method to improve the linearised AVO approximation, based on the concept that the inaccuracy of the conventional linearised AVO forward operator is solely responsible for generating residuals of the PP reflection coefficients calculated by the conventional linearised approximation compared with using the Zoeppritz equation. Therefore, if true model parameters from well-logging data are known, the residuals of the PP reflection coefficients can be used to estimate linear modifications of the inaccurate linearised AVO forward operators at well sites. We apply estimated linear modifications to original inaccurately linearised forward operators to obtain an improved linear AVO operator (ILAO). Because ILAOs can only be estimated at the well sites by the data-driven method, we construct spatially variant ILAOs using structure-guided interpolation for those at a distance from the well sites. Numerical example results reveal that the improved linear AVO approximation is more accurate than the Aki–Richards approximation for strong-contrast media and large incident angles. Moreover, the accuracy and resolution of the inverted P- and S-wave velocities and density using ILAOs are higher than those using the conventional Aki–Richards operators. Finally, application to the field data successfully demonstrates the feasibility and effectiveness of the improved linearised AVO inversion method.


Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1426-1436 ◽  
Author(s):  
Wojciech Dȩbski ◽  
Albert Tarantola

Seismic amplitude variation with offset data contain information on the elastic parameters of geological layers. As the general solution of the inverse problem consists of a probability over the space of all possible earth models, we look at the probabilities obtained using amplitude variation with offset (AVO) data for different choices of elastic parameters. A proper analysis of the information in the data requires a nontrivial definition of the probability defining the state of total ignorance on different elastic parameters (seismic velocities, Lamé’s parameters, etc.). We conclude that mass density, seismic impedance, and Poisson’s ratio constitute the best resolved parameter set when inverting seismic amplitude variation with offset data.


2021 ◽  
Vol 40 (9) ◽  
pp. 646-654
Author(s):  
Henning Hoeber

When inversions use incorrectly specified models, the estimated least-squares model parameters are biased. Their expected values are not the true underlying quantitative parameters being estimated. This means the least-squares model parameters cannot be compared to the equivalent values from forward modeling. In addition, the bias propagates into other quantities, such as elastic reflectivities in amplitude variation with offset (AVO) analysis. I give an outline of the framework to analyze bias, provided by the theory of omitted variable bias (OVB). I use OVB to calculate exactly the bias due to model misspecification in linearized isotropic two-term AVO. The resulting equations can be used to forward model unbiased AVO quantities, using the least-squares fit results, the weights given by OVB analysis, and the omitted variables. I show how uncertainty due to bias propagates into derived quantities, such as the χ-angle and elastic reflectivity expressions. The result can be used to build tables of unique relative rock property relationships for any AVO model, which replace the unbiased, forward-model results.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. E49-E55 ◽  
Author(s):  
Jonathan E. Downton ◽  
Charles Ursenbach

Contrary to popular belief, a linearized approximation of the Zoeppritz equations may be used to estimate the reflection coefficient for angles of incidence up to and beyond the critical angle. These supercritical reflection coefficients are complex, implying a phase variation with offset in addition to amplitude variation with offset (AVO). This linearized approximation is then used as the basis for an AVO waveform inversion. By incorporating this new approximation, wider offset and angle data may be incorporated in the AVO inversion, helping to stabilize the problem and leading to more accurate estimates of reflectivity, including density reflectivity.


Sign in / Sign up

Export Citation Format

Share Document