Enabling isotropic reflectivity modeling and inversion in anisotropic media

2021 ◽  
Vol 40 (4) ◽  
pp. 267-276
Author(s):  
Peter Mesdag ◽  
Leonardo Quevedo ◽  
Cătălin Tănase
Keyword(s):  
Synthetic Data ◽  
Anisotropic Media ◽  
Forward Modeling ◽  
Seismic Inversion ◽  
Stratified Medium ◽  
Effective Parameters ◽  
Elastic Parameters ◽  
Pp Wave ◽  
Anisotropic Elastic ◽  
And Inversion

Exploration and development of unconventional reservoirs, where fractures and in-situ stresses play a key role, call for improved characterization workflows. Here, we expand on a previously proposed method that makes use of standard isotropic modeling and inversion techniques in anisotropic media. Based on approximations for PP-wave reflection coefficients in orthorhombic media, we build a set of transforms that map the isotropic elastic parameters used in prestack inversion into effective anisotropic elastic parameters. When used in isotropic forward modeling and inversion, these effective parameters accurately mimic the anisotropic reflectivity behavior of the seismic data, thus closing the loop between well-log data and seismic inversion results in the anisotropic case. We show that modeling and inversion of orthorhombic anisotropic media can be achieved by superimposing effective elastic parameters describing the behavior of a horizontally stratified medium and a set of parallel vertical fractures. The process of sequential forward modeling and postinversion analysis is exemplified using synthetic data.

Prestack Gaussian-beam depth migration in anisotropic media

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. S133-S138 ◽  
Author(s):  
Tianfei Zhu ◽  
Samuel H. Gray ◽  
Daoliu Wang
Keyword(s):  
Gaussian Beam ◽  
Ray Tracing ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Elastic Parameters ◽  
Depth Migration ◽  
Kirchhoff Migration ◽  
Beam Depth ◽  
Dynamic Ray Tracing

Gaussian-beam depth migration is a useful alternative to Kirchhoff and wave-equation migrations. It overcomes the limitations of Kirchhoff migration in imaging multipathing arrivals, while retaining its efficiency and its capability of imaging steep dips with turning waves. Extension of this migration method to anisotropic media has, however, been hampered by the difficulties in traditional kinematic and dynamic ray-tracing systems in inhomogeneous, anisotropic media. Formulated in terms of elastic parameters, the traditional anisotropic ray-tracing systems aredifficult to implement and inefficient for computation, especially for the dynamic ray-tracing system. They may also result inambiguity in specifying elastic parameters for a given medium.To overcome these difficulties, we have reformulated the ray-tracing systems in terms of phase velocity.These reformulated systems are simple and especially useful for general transversely isotropic and weak orthorhombic media, because the phase velocities for these two types of media can be computed with simple analytic expressions. These two types of media also represent the majority of anisotropy observed in sedimentary rocks. Based on these newly developed ray-tracing systems, we have extended prestack Gaussian-beam depth migration to general transversely isotropic media. Test results with synthetic data show that our anisotropic, prestack Gaussian-beam migration is accurate and efficient. It produces images superior to those generated by anisotropic, prestack Kirchhoff migration.

Forward modeling and inversion of tensor CSAMT in 3D anisotropic media

2017 ◽  
Vol 14 (4) ◽  
pp. 590-605 ◽  
Author(s):  
Tao Wang ◽  
Kun-Peng Wang ◽  
Han-Dong Tan
Keyword(s):  
Anisotropic Media ◽  
Forward Modeling ◽  
And Inversion

scholarly journals DETERMINATION OF AZIMUTHAL ANISOTROPIC MEDIA ELASTIC PARAMETERS FROM MULTIWAVE AVOA DATA BY NONLINEAR OPTIMIZATION METHOD

2019 ◽  
pp. 14-26 ◽  
Author(s):  
T. V. Nefedkina ◽  
P. A. Lykhin ◽  
G. A. Dugarov
Keyword(s):  
Signal To Noise Ratio ◽  
Anisotropic Media ◽  
Optimization Method ◽  
Parameter Estimates ◽  
Reflection Coefficients ◽  
Elastic Parameters ◽  
Ray Method ◽  
Pp Wave ◽  
Hti Medium

In this paper, we investigate optimization algorithm of joint nonlinear AVOA inversion of PP+PS reflections in anisotropic media. Algorithm is based on the exact solution for PP and PS waves reflection coefficients in anisotropic HTI medium. The PP and PS wave’s reflections from the top of the anisotropic layer are examined. We use synthetic seismograms generated by ray method for the algorithm testing. We show that joint compressional and converted wave’s inversion allows increasing the robustness of the method and the accuracy of medium-parameter estimates. Coefficients of anisotropy are determined with better accuracy if signal-to-noise ratio is bigger than 5 for PP wave and bigger than 2 for PS wave.

Azimuthal amplitude variation with offset parameterization and inversion for fracture weaknesses in tilted transversely isotropic media

Geophysics ◽  
2020 ◽  
Vol 86 (1) ◽  
pp. C1-C18
Author(s):  
Xinpeng Pan ◽  
Lin Li ◽  
Shunxin Zhou ◽  
Guangzhi Zhang ◽  
Jianxin Liu
Keyword(s):  
Tilt Angle ◽  
Stiffness Matrix ◽  
Incidence Angle ◽  
Fractured Reservoirs ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Pp Wave ◽  
And Inversion

The characterization of fracture-induced tilted transverse isotropy (TTI) seems to be more suitable to actual scenarios of geophysical exploration for fractured reservoirs. Fracture weaknesses enable us to describe fracture-induced anisotropy. With the incident and reflected PP-wave in TTI media, we have adopted a robust method of azimuthal amplitude variation with offset (AVO) parameterization and inversion for fracture weaknesses in a fracture-induced reservoir with TTI symmetry. Combining the linear-slip model with the Bond transformation, we have derived the stiffness matrix of a dipping-fracture-induced TTI medium characterized by normal and tangential fracture weaknesses and a tilt angle. Integrating the first-order perturbations in the stiffness matrix of a TTI medium and scattering theory, we adopt a method of azimuthal AVO parameterization for PP-wave reflection coefficient for the case of a weak-contrast interface separating two homogeneous weakly anisotropic TTI layers. We then adopt an iterative inversion method by using the partially incidence-angle-stacked seismic data with different azimuths to estimate the fracture weaknesses of a TTI medium when the tilt angle is estimated based on the image well logs prior to the seismic inversion. Synthetic examples confirm that the fracture weaknesses of a TTI medium are stably estimated from the azimuthal seismic reflected amplitudes for the case of moderate noise. A field data example demonstrates that geologically reasonable results of fracture weaknesses can be determined when the tilt angle is treated as the prior information. We determine that the azimuthal AVO inversion approach provides an available tool for fracture prediction in a dipping-fracture-induced TTI reservoir.

Ensemble-based seismic inversion for a stratified medium

Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. R29-R39 ◽  
Author(s):  
Michael Gineste ◽  
Jo Eidsvik ◽  
York Zheng
Keyword(s):  
Waveform Inversion ◽  
Seismic Inversion ◽  
Stratified Medium ◽  
Elastic Parameters ◽  
Waveform Data ◽  
Optimization Task ◽  
Estimation Uncertainty ◽  
Seismic Waveform ◽  
Probabilistic Setting ◽  
Iterative Ensemble Smoother

Seismic waveform inversion is a nontrivial optimization task, which is often complicated by the nonlinear relationship between the elastic attributes of interest and the large amount of data obtained in seismic experiments. Quantifying the solution uncertainty can be even more challenging, and it requires considering the problem in a probabilistic setting. Consequently, the seismic inverse problem is placed in a Bayesian framework, using a sequential filtering approach to invert for the elastic parameters. The method uses an iterative ensemble smoother to estimate the subsurface parameters, and from the ensemble, a notion of estimation uncertainty is readily available. The ensemble implicitly linearizes the relation between the parameters and the observed waveform data; hence, it requires no tangent linear model. The approach is based on sequential conditioning on partitions of the whole data record (1) to regularize the inversion path and effectively drive the estimation process in a top-down manner and (2) to circumvent a consequence of the ensemble reduced rank approximation. The method is exemplified on a synthetic case, inverting for elastic parameters in a 1D medium using a seismic shot record. Our results indicate that the iterative ensemble method is applicable to seismic waveform inversion and that the ensemble representation indeed indicates estimation uncertainty.

2.5-D forward modeling and inversion of time domain IP method

2017 ◽  
Vol 70 (0) ◽  
pp. 69-79
Author(s):  
Hideki Mizunaga ◽  
Kiyotaka Ishinaga
Keyword(s):  
Time Domain ◽  
Forward Modeling ◽  
And Inversion

Extension of forward modeling phase-screen code in isotropic and anisotropic media up to critical angle

Geophysics ◽  
2007 ◽  
Vol 72 (5) ◽  
pp. SM107-SM114 ◽  
Author(s):  
James C. White ◽  
Richard W. Hobbs
Keyword(s):  
Critical Angle ◽  
Model Space ◽  
Anisotropic Media ◽  
Forward Modeling ◽  
Phase Screen ◽  
Computationally Efficient ◽  
Reflection And Transmission ◽  
Transmission Coefficients ◽  
Phase Corrections ◽  
Converted Phases

The computationally efficient phase-screen forward modeling technique is extended to allow investigation of nonnormal raypaths. The code is developed to accommodate all diffracted and converted phases up to critical angle, building on a geometric construction method. The new approach relies upon prescanning the model space to assess the complexity of each screen. The propagating wavefields are then divided as a function of horizontal wavenumber, and each subset is transformed to the spatial domain separately, carrying with it angular information. This allows both locally accurate 3D phase corrections and Zoeppritz reflection and transmission coefficients to be applied. The phase-screen code is further developed to handle simple anisotropic media. During phase-screen modeling, propagation is undertaken in the wavenumber domain where exact expressions for anisotropic phase velocities are available. Traveltimes and amplitude effects from a range of anisotropic shales are computed and compared with previous published results.

Sign in / Sign up

Export Citation Format

Share Document