The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,w...The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.展开更多
The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors...The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors,vibration of axisymmetric infinite membrane and displacement intensity factors which all involve this type of integration.However,traditional numerical integration algorithms cannot be used due to the high oscillation characteristics of the Bessel function,so it is particularly important to propose a high precision and efficient numerical algorithm for calculating the integral of high oscillation.In this paper,the improved Gaver-Stehfest(G-S)inverse Laplace transform method for arbitrary real-order Bessel function integration is presented by using the asymptotic characteristics of the Bessel function and the accumulation of integration,and the optimized G-S coefficients are given.The effectiveness of the algorithm is verified by numerical examples.Compared with the linear transformation accelerated convergence algorithm,it shows that the G-S inverse Laplace transform method is suitable for arbitrary real order Hankel transform,and the time consumption is relatively stable and short,which provides a reliable calculation method for the study of electromagnetic mechanics,wave propagation,and fracture dynamics.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant No.42104131)the Natural Science Foundation of Sichuan Province of China(Grant No.2022NSFSC1140)Open Fund(PLC20211101)of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation
文摘The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.
基金Supported by the National Natural Science Foundation of China(42064004,12062022,11762017,11762016)
文摘The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors,vibration of axisymmetric infinite membrane and displacement intensity factors which all involve this type of integration.However,traditional numerical integration algorithms cannot be used due to the high oscillation characteristics of the Bessel function,so it is particularly important to propose a high precision and efficient numerical algorithm for calculating the integral of high oscillation.In this paper,the improved Gaver-Stehfest(G-S)inverse Laplace transform method for arbitrary real-order Bessel function integration is presented by using the asymptotic characteristics of the Bessel function and the accumulation of integration,and the optimized G-S coefficients are given.The effectiveness of the algorithm is verified by numerical examples.Compared with the linear transformation accelerated convergence algorithm,it shows that the G-S inverse Laplace transform method is suitable for arbitrary real order Hankel transform,and the time consumption is relatively stable and short,which provides a reliable calculation method for the study of electromagnetic mechanics,wave propagation,and fracture dynamics.
