scholarly journals Impact of poroelastic effects on the inversion of fracture properties from amplitude variation with offset and azimuth data in horizontal transversely isotropic media

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. N27-N39
Author(s):  
Nicolás D. Barbosa ◽  
Corinna Köpke ◽  
Eva Caspari ◽  
J. Germán Rubino ◽  
James Irving ◽  
...  

The identification and characterization of fractures is an important objective in many areas of earth and environmental sciences. Amplitude variation with offset and azimuth (AVOAz) analysis of seismic reflection data is a key method for achieving these tasks. Theoretical and experimental studies have shown that the presence of pore fluids together with the strong mechanical contrast between the fractures and their embedding background give rise to wave-induced fluid flow (WIFF) effects. This implies that the effective stiffness tensor of a fluid-saturated fractured rock defining its seismic response becomes viscoelastic and frequency-dependent. In spite of this, AVOAz analysis typically relies on end-member-type elastic stiffness models that either assume a relaxed (i.e., equilibrated) or unrelaxed (i.e., unequilibrated) state of the wave-induced fluid pressure in the rock. In general, however, neither the appropriateness of the chosen model nor the associated errors in the inversion process are known. To investigate this topic, we have considered a poroelastic medium containing parallel vertical fractures and generate synthetic seismic AVOAz data using the classic Rüger approximations for PP-wave reflection coefficients in horizontally transversely isotropic media. A Markov chain Monte Carlo method is used to perform a Bayesian inversion of the synthetic seismic AVOAz data. We quantify the influence of WIFF effects on the AVOAz inversion results when elastic relaxed and unrelaxed models are used as forward solvers of inversion schemes to estimate the fracture volume fraction, the elastic moduli, and the porosity of the background rock, as well as the overall weakness of the medium due to the presence of fractures. Our results indicate that, when dealing with single-frequency data, relaxed elastic models provide biased but overall better inversion results than unrelaxed ones, for which some fracture parameters cannot be resolved. Improved inversion performance is achieved when using frequency-dependent data, which illustrates the importance of accounting for poroelastic effects.

2019 ◽  
Vol 220 (2) ◽  
pp. 839-855
Author(s):  
Da Shuai ◽  
Alexey Stovas

SUMMARY We develop a method to compute frequency-dependent kinematic parameters for an effective orthorhombic (ORT) medium. In order to investigate the influence of fracture weaknesses on the kinematic parameters, the effective ORT medium is composed based on the linear slip theory and derived by applying the limited Baker–Campbell–Hausdorff series. The frequency-dependent kinematic parameters including vertical velocity, two normal moveout velocities defined in vertical symmetry planes, and three anelliptic parameters (two of them are defined in vertical symmetry plane and one parameter is the cross-term one). We also investigate the influence of volume fraction, frequency, velocity ratio and fracture weaknesses on the effective kinematic parameters.


2020 ◽  
Vol 68 (8) ◽  
pp. 2471-2493
Author(s):  
Amir Abbas Babasafari ◽  
Deva Prasad Ghosh ◽  
Ahmed Mohammed Ahmed Salim ◽  
Masoumeh Kordi

Geophysics ◽  
1994 ◽  
Vol 59 (5) ◽  
pp. 775-781 ◽  
Author(s):  
J. P. Blangy

The amplitude variation with offset (AVO) responses of elastic transversely isotropic media are sensitive to contrasts in both of Thomsen’s anisotropic parameters δ and ε. The equation describing P-P reflections indicates that the smaller the contrasts in isotropic properties (compressional velocity, shear velocity, and density) and the larger the contrasts in δ and ε across an interface of reflection, the greater the effects of anisotropy on the AVO signature. Contrasts in δ are most important under small‐to‐medium angles of incidence as previously described by Banik (1987), while contrasts in ε can have a strong influence on amplitudes for the larger angles of incidence commonly encountered in exploration (20 degrees and beyond). Consequently, using Rutherford and Williams’ AVO classification of isotropic gas sands, type I gas sands overlain by a transversely isotropic (TI) shale exhibit a larger decrease in AVO than if the shale had been isotropic, and type III gas sands overlain by a transversely isotropic (TI) shale exhibit a larger increase in AVO than if the shale had been isotropic. Furthermore, it is possible for a “type III” isotropic water sand to exhibit an “unexpected) increase in AVO if the overlying shale is sufficiently anisotropic. More quantitative AVO interpretations in TI media require considerations of viscoelastic TI in addition to elastic TI and lead to complicated integrated earth models. However, when elastic and viscoelastic TI have the same axis of symmetry in a shale overlying an isotropic sand, both elastic and viscoelastic TI affect the overall AVO response in the same direction by constructively increasing/decreasing the isotropic component of the AVO response. Continued efforts in this area will lead to more realistic reservoir models and hopefully answer some of the unexplained pitfalls in AVO interpretation.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB33-WB53 ◽  
Author(s):  
Milana Ayzenberg ◽  
Ilya Tsvankin ◽  
Arkady Aizenberg ◽  
Bjørn Ursin

Plane-wave reflection coefficients (PWRCs) are routinely used in amplitude-variation-with-offset analysis and for generating boundary data in Kirchhoff modeling. However, the geometrical-seismics approximation based on PWRCs becomes inadequate in describing reflected wavefields at near- and postcritical incidence angles. Also, PWRCs are derived for plane interfaces and break down in the presence of significant reflector curvature. Here, we discuss effective reflection coefficients (ERCs) designed to overcome the limitations of PWRCs for multicomponent data from heterogeneous anisotropic media. We represent the reflected wavefield in the immediate vicinity of a curved interface by a generalized plane-wave decomposition, which approximately reduces to the conventional Weyl-type integral computed for an apparent source location. The ERC then is obtainedas the ratio of the reflected and incident wavefields at each point of the interface. To conduct diffraction modeling, we combine ERCs with the tip-wave superposition method (TWSM), extended to elastic media. This methodology is implemented for curved interfaces that separate an isotropic incidence half-space and a transversely isotropic (TI) medium with the symmetry axis orthogonal to the reflector. If the interface is plane, ERCs generally are close to the exact solution, sensitive to the anisotropy parameters and source-receiver geometry. Numerical tests demonstrate that the difference between ERCs and PWRCs for typical TI models can be significant, especially at low frequencies and in the postcritical domain. For curved interfaces, ERCs provide a practical approximate tool to compute the reflected wavefield. We analyze the dependence of ERCs on reflector shape and demonstrate their advantages over PWRCs in 3D diffraction modeling of PP and PS reflection data.


2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo

ABSTRACT. Reverse time migration (RTM) is one of the most powerful methods used to generate images of the subsurface. The RTM was proposed in the early 1980s, but only recently it has been routinely used in exploratory projects involving complex geology – Brazilian pre-salt, for example. Because the method uses the two-way wave equation, RTM is able to correctly image any kind of geological environment (simple or complex), including those with anisotropy. On the other hand, RTM is computationally expensive and requires the use of computer clusters. This paper proposes to investigate the influence of anisotropy on seismic imaging through the application of RTM for tilted transversely isotropic (TTI) media in pre-stack synthetic data. This work presents in detail how to implement RTM for TTI media, addressing the main issues and specific details, e.g., the computational resources required. A couple of simple models results are presented, including the application to a BP TTI 2007 benchmark model.Keywords: finite differences, wave numerical modeling, seismic anisotropy. Migração reversa no tempo em meios transversalmente isotrópicos inclinadosRESUMO. A migração reversa no tempo (RTM) é um dos mais poderosos métodos utilizados para gerar imagens da subsuperfície. A RTM foi proposta no início da década de 80, mas apenas recentemente tem sido rotineiramente utilizada em projetos exploratórios envolvendo geologia complexa, em especial no pré-sal brasileiro. Por ser um método que utiliza a equação completa da onda, qualquer configuração do meio geológico pode ser corretamente tratada, em especial na presença de anisotropia. Por outro lado, a RTM é dispendiosa computacionalmente e requer o uso de clusters de computadores por parte da indústria. Este artigo apresenta em detalhes uma implementação da RTM para meios transversalmente isotrópicos inclinados (TTI), abordando as principais dificuldades na sua implementação, além dos recursos computacionais exigidos. O algoritmo desenvolvido é aplicado a casos simples e a um benchmark padrão, conhecido como BP TTI 2007.Palavras-chave: diferenças finitas, modelagem numérica de ondas, anisotropia sísmica.


1995 ◽  
Vol 26 (2-3) ◽  
pp. 431-436 ◽  
Author(s):  
Patrick N.(Jr). Okoye ◽  
N. F. Uren ◽  
W. Waluyo

2020 ◽  
Vol 41 (4) ◽  
pp. 805-833 ◽  
Author(s):  
Jidong Yang ◽  
Biaolong Hua ◽  
Paul Williamson ◽  
Hejun Zhu ◽  
George McMechan ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document