In VTI media,the conventional inversion methods based on the existing approximation formulas are difficult to accurately estimate the anisotropic parameters of reservoirs,even more so for unconventional reservoirs wit...In VTI media,the conventional inversion methods based on the existing approximation formulas are difficult to accurately estimate the anisotropic parameters of reservoirs,even more so for unconventional reservoirs with strong seismic anisotropy.Theoretically,the above problems can be solved by utilizing the exact reflection coefficients equations.However,their complicated expression increases the difficulty in calculating the Jacobian matrix when applying them to the Bayesian deterministic inversion.Therefore,the new reduced approximation equations starting from the exact equations are derived here by linearizing the slowness expressions.The relatively simple form and satisfactory calculation accuracy make the reduced equations easy to apply for inversion while ensuring the accuracy of the inversion results.In addition,the blockiness constraint,which follows the differentiable Laplace distribution,is added to the prior model to improve contrasts between layers.Then,the concept of GLI and an iterative reweighted least-squares algorithm is combined to solve the objective function.Lastly,we obtain the iterative solution expression of the elastic parameters and anisotropy parameters and achieve nonlinear AVA inversion based on the reduced equations.The test results of synthetic data and field data show that the proposed method can accurately obtain the VTI parameters from prestack AVA seismic data.展开更多
Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion...Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion on the propagation of plane waves is given.Elastic plane waves are characterized by slowness vectors, wave vectors, polarization vectors and group velocity vectors, etc. The results obtained are presented in a coordinate-free form due to the introduction of the crystal axis' orieniation vector which specifies the anisotropy of the media. Therefore, the results are the most general and convenient for further application展开更多
Fourier transform method is used to obtain an approximate solution of Green's tensor to homogeneous and transversely isotropic media like unidirectional fiber re-inforced composites and austenitic stainless steel...Fourier transform method is used to obtain an approximate solution of Green's tensor to homogeneous and transversely isotropic media like unidirectional fiber re-inforced composites and austenitic stainless steel materials in order to provide the theoretical basis for the scattering problems. A comparison to homogeneously isotropic media is presented and a brief discussion of the main features of the solution is given展开更多
Two-dimensional scalar equation for the displacement of steady cross-plane shear (SH) waves in homogeneous and transversely isotropic media like unidirectional fibrous com-posites is given. Then, thrbugh a simple coor...Two-dimensional scalar equation for the displacement of steady cross-plane shear (SH) waves in homogeneous and transversely isotropic media like unidirectional fibrous com-posites is given. Then, thrbugh a simple coordinate system transform, the scalar equation is standardized into a Helmholtz equation. Corresponding integral equations are derived for the scattering problems and boundary element method (BEM) is used to calculate the scattered fields of arbitrarily shaped obstacles with both soft and rigid boudary conditions numerically.A discussion is given on the numerical results which is mainly focused on the influence of the a-nisotropy of the media to the directivity of the scattered fields by circular cylindrical voids.展开更多
The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of ...The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot's wave equations. The method of research on non-axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.展开更多
Dynamic behavior of single pile embedded in transversely isotropic layered media is investigated using the finite element method combined with dynamic stiffness matrices of the soil derived from Green's function f...Dynamic behavior of single pile embedded in transversely isotropic layered media is investigated using the finite element method combined with dynamic stiffness matrices of the soil derived from Green's function for ring loads. The influence of soil anisotropy on the dynamic behavior of piles is examined through a series of parametric studies.展开更多
The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled ...The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.展开更多
In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on...In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on homogenous and isotropic materials.However,soil and rock are often modeled in heterogeneous and anisotropic media because of the inherent anisotropy of sedimentary rock masses and their stratified structure.In recent decades,coupled hydro-mechanical(HM)interactions in isotropic porous media have been studied;however,the behavior of transversely isotropic porous media is rarely considered.In addition,it is difficult for commercial software such as Plaxis and Flac3D to express complex rock formation where the anisotropy of the material and the associated cracks and fractures could be assembled into a single model.In this study,a finite element implementation using Differential Equation Analysis Library(DEAL.II),an open-source library of finite element codes,was developed to model the fully coupled HM behavior of transversely isotropic porous media.The proposed implementation can be applied to both isotropic and transversely isotropic porous media based on Biot’s theory.The developed code can be used to model poroelastic media with(1)equations of linear elasticity for the solid matrix and(2)diffusion equations for fluid flow based on mass and linear-momentum conservation laws.We verified the performance and accuracy of the code through two examples,i.e.,Mandel’s problem with a compared analytical solution and a tunnel excavation process with the Flac3D software.On the basis of these numerical applications,we present the code to model the behavior of various geo-structures such as tunnels and pile–soil interactions with anisotropic materials.展开更多
文摘In VTI media,the conventional inversion methods based on the existing approximation formulas are difficult to accurately estimate the anisotropic parameters of reservoirs,even more so for unconventional reservoirs with strong seismic anisotropy.Theoretically,the above problems can be solved by utilizing the exact reflection coefficients equations.However,their complicated expression increases the difficulty in calculating the Jacobian matrix when applying them to the Bayesian deterministic inversion.Therefore,the new reduced approximation equations starting from the exact equations are derived here by linearizing the slowness expressions.The relatively simple form and satisfactory calculation accuracy make the reduced equations easy to apply for inversion while ensuring the accuracy of the inversion results.In addition,the blockiness constraint,which follows the differentiable Laplace distribution,is added to the prior model to improve contrasts between layers.Then,the concept of GLI and an iterative reweighted least-squares algorithm is combined to solve the objective function.Lastly,we obtain the iterative solution expression of the elastic parameters and anisotropy parameters and achieve nonlinear AVA inversion based on the reduced equations.The test results of synthetic data and field data show that the proposed method can accurately obtain the VTI parameters from prestack AVA seismic data.
文摘Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion on the propagation of plane waves is given.Elastic plane waves are characterized by slowness vectors, wave vectors, polarization vectors and group velocity vectors, etc. The results obtained are presented in a coordinate-free form due to the introduction of the crystal axis' orieniation vector which specifies the anisotropy of the media. Therefore, the results are the most general and convenient for further application
文摘Fourier transform method is used to obtain an approximate solution of Green's tensor to homogeneous and transversely isotropic media like unidirectional fiber re-inforced composites and austenitic stainless steel materials in order to provide the theoretical basis for the scattering problems. A comparison to homogeneously isotropic media is presented and a brief discussion of the main features of the solution is given
文摘Two-dimensional scalar equation for the displacement of steady cross-plane shear (SH) waves in homogeneous and transversely isotropic media like unidirectional fibrous com-posites is given. Then, thrbugh a simple coordinate system transform, the scalar equation is standardized into a Helmholtz equation. Corresponding integral equations are derived for the scattering problems and boundary element method (BEM) is used to calculate the scattered fields of arbitrarily shaped obstacles with both soft and rigid boudary conditions numerically.A discussion is given on the numerical results which is mainly focused on the influence of the a-nisotropy of the media to the directivity of the scattered fields by circular cylindrical voids.
文摘The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot's wave equations. The method of research on non-axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.
文摘Dynamic behavior of single pile embedded in transversely isotropic layered media is investigated using the finite element method combined with dynamic stiffness matrices of the soil derived from Green's function for ring loads. The influence of soil anisotropy on the dynamic behavior of piles is examined through a series of parametric studies.
基金financially supported by the National Key R&D Program of China (No. 2019YFC0605503)the Major Scientific and Technological Projects of CNPC (No. ZD2019-183-003)the National Natural Science Foundation of China (No. 41922028,41874149)。
文摘The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.
文摘In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on homogenous and isotropic materials.However,soil and rock are often modeled in heterogeneous and anisotropic media because of the inherent anisotropy of sedimentary rock masses and their stratified structure.In recent decades,coupled hydro-mechanical(HM)interactions in isotropic porous media have been studied;however,the behavior of transversely isotropic porous media is rarely considered.In addition,it is difficult for commercial software such as Plaxis and Flac3D to express complex rock formation where the anisotropy of the material and the associated cracks and fractures could be assembled into a single model.In this study,a finite element implementation using Differential Equation Analysis Library(DEAL.II),an open-source library of finite element codes,was developed to model the fully coupled HM behavior of transversely isotropic porous media.The proposed implementation can be applied to both isotropic and transversely isotropic porous media based on Biot’s theory.The developed code can be used to model poroelastic media with(1)equations of linear elasticity for the solid matrix and(2)diffusion equations for fluid flow based on mass and linear-momentum conservation laws.We verified the performance and accuracy of the code through two examples,i.e.,Mandel’s problem with a compared analytical solution and a tunnel excavation process with the Flac3D software.On the basis of these numerical applications,we present the code to model the behavior of various geo-structures such as tunnels and pile–soil interactions with anisotropic materials.