Anisotropic Anelastic Impedance Inversion for Attenuation Factor and Weaknesses Combining Newton and Bayesian Markov Chain Monte Carlo Algorithms

Huaizhen Chen,Jian Han,Jianhua Geng,Kristopher A. Innanen
DOI: https://doi.org/10.1007/s10712-022-09727-y
2022-01-01
Surveys in Geophysics
Abstract:The problem of frequency-dependent seismic inversion for estimating attenuation and anisotropy based on rock physics effective model is reviewed and analyzed. Based on an extended periodically layered medium model, frequency-dependent stiffness parameters of anisotropic anelastic media as a function of attenuation factor and weaknesses are presented. For the case of an interface separating two anisotropic anelastic media, a new parameterized P -to- P reflection coefficient and anisotropic anelastic impedance 𝒜_EI are proposed. Based on the anisotropic 𝒜_EI , an inversion workflow of employing frequency components of seismic amplitudes to estimate attenuation factor and weaknesses for anisotropic anelastic media is presented. In the first stage inversion, the anisotropic 𝒜_EI is calculated using frequency components of partially incidence-angle-stacked seismic data; and the second-stage inversion for unknown attenuation factor and weaknesses is implemented combining a Newton method and Bayesian Markov chain Monte Carlo algorithm by calculating first- and second-order derivatives of anisotropic 𝒜_EI with respect to unknown parameters. Tests on synthetic seismic gathers confirm the accuracy of new parameterized reflection coefficient, and the unknown parameters are estimated reliably even in the case of seismic data of signal-to-noise ratio of 2 being employed. Applying the inversion approach and workflow to a field data set acquired over a shale reservoir, reliable attenuation factors and weaknesses that match well log curves are obtained, which are preserved as indicators for identifying potential hydrocarbon-bearing cracked reservoirs.
What problem does this paper attempt to address?