Based on the long-wavelength approximation,a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the ...Based on the long-wavelength approximation,a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study,we propose a reliable method for predicting the Thomsen's weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen's weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen's weak anisotropic parameters and fracture weaknesses.Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth(AVAZ) inversion for Thomsen's weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares(IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen's weak anisotropic parameters and fracture weaknesses.展开更多
Sedimentary strata typically exhibit the characteristics of transverse isotropy(VTI)with a vertical axis of symmetry.However,fractures in sedimentary strata tend to produce anisotropic closure due to horizontal in sit...Sedimentary strata typically exhibit the characteristics of transverse isotropy(VTI)with a vertical axis of symmetry.However,fractures in sedimentary strata tend to produce anisotropic closure due to horizontal in situ stress,resulting in pronounced orthorhombic anisotropy in VTI media under such stress conditions and influencing the propagation behavior of seismic waves.Previous studies have primarily focused on the elastic wave velocity anisotropy induced by applied stress in isotropic background media,neglecting the impact of VTI background media on the anisotropy induced by horizontal in situ stress and the response characteristics of seismic wave propagation.To address these gaps,we first establish the effective elastic stiffness tensor of VTI media under horizontal in situ stress using nonlinear acoustoelastic theory.Then,we derive the accurate and linearized approximate equations for P-wave seismic reflectivity of VTI media under horizontal in situ stress,based on wave equations and scattering theory,respectively.Finally,we compare and analyze the characteristics of orthorhombic anisotropic seismic response induced by horizontal in situ stress at various types of elastic reflection interfaces.Our results demonstrate that the linearized approximation of the seismic reflection response characteristics closely aligns with the accurate equations under conditions of small stress below 10 MPa,effectively capturing the azimuth-dependent orthorhombic anisotropy induced by horizontal in situ stress in VTI media.The results of this study also provide a novel theoretical approach and valuable insights into the seismic prediction of in situ stress.展开更多
文摘Based on the long-wavelength approximation,a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study,we propose a reliable method for predicting the Thomsen's weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen's weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen's weak anisotropic parameters and fracture weaknesses.Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth(AVAZ) inversion for Thomsen's weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares(IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen's weak anisotropic parameters and fracture weaknesses.
基金supported by the National Natural Science Foundation of China(Grant Nos.42130810,42004107)the Science and Technology Innovation Program of Hunan Province(Grant No.2022RC1238)。
文摘Sedimentary strata typically exhibit the characteristics of transverse isotropy(VTI)with a vertical axis of symmetry.However,fractures in sedimentary strata tend to produce anisotropic closure due to horizontal in situ stress,resulting in pronounced orthorhombic anisotropy in VTI media under such stress conditions and influencing the propagation behavior of seismic waves.Previous studies have primarily focused on the elastic wave velocity anisotropy induced by applied stress in isotropic background media,neglecting the impact of VTI background media on the anisotropy induced by horizontal in situ stress and the response characteristics of seismic wave propagation.To address these gaps,we first establish the effective elastic stiffness tensor of VTI media under horizontal in situ stress using nonlinear acoustoelastic theory.Then,we derive the accurate and linearized approximate equations for P-wave seismic reflectivity of VTI media under horizontal in situ stress,based on wave equations and scattering theory,respectively.Finally,we compare and analyze the characteristics of orthorhombic anisotropic seismic response induced by horizontal in situ stress at various types of elastic reflection interfaces.Our results demonstrate that the linearized approximation of the seismic reflection response characteristics closely aligns with the accurate equations under conditions of small stress below 10 MPa,effectively capturing the azimuth-dependent orthorhombic anisotropy induced by horizontal in situ stress in VTI media.The results of this study also provide a novel theoretical approach and valuable insights into the seismic prediction of in situ stress.
