球面波的反射P波AVO分析

本文对球面波反射P波反射系数计算公式进行了推导,根据推导公式计算出第I类AVO和第Ⅲ类AVO介质中反射系数及相位随入射角的变化,并与基于平面波的Zoeppritz方程计算的结果进行对比.结果显示,对于第I类AVO介质,球面波反射系数在临界角附近较基于平面波的Zoeppritz方程计算结果更精确;同时受界面深度的影响较大,随深度的增加,球面波AVO趋势接近平面波AVO.最后,通过数值模拟技术,对模型数据和实际数据进行模拟,对平面波AVO道集和球面波AVO道集进行了验证分析.浅层的第I类AVO现象,在近临界角和超临界角处折射引起反射系数与相位的变化较大,对于第Ⅲ类AVO现象,由于不存在临界角问题,球面波模拟结果与基于平面波的Zoeppritz计算结果差别较小.上述计算分析,可为实际资料的大偏移距道集的AVO分析提供理论基础.     

Weak‐anisotropy approximation for P‐wave reflection coefficient at the boundary between two tilted transversely isotropic media

Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy.     

Elastic full waveform inversion of multi‐component OBC seismic data

Elastic full waveform inversion of seismic reflection data represents a data‐driven form of analysis leading to quantification of sub‐surface parameters in depth. In previous studies attention has been given to P‐wave data recorded in the marine environment, using either acoustic or elastic inversion schemes. In this paper we exploit both P‐waves and mode‐converted S‐waves in the marine environment in the inversion for both P‐ and S‐wave velocities by using wide‐angle, multi‐component, ocean‐bottom cable seismic data. An elastic waveform inversion scheme operating in the time domain was used, allowing accurate modelling of the full wavefield, including the elastic amplitude variation with offset response of reflected arrivals and mode‐converted events. A series of one‐ and two‐dimensional synthetic examples are presented, demonstrating the ability to invert for and thereby to quantify both P‐ and S‐wave velocities for different velocity models. In particular, for more realistic low velocity models, including a typically soft seabed, an effective strategy for inversion is proposed to exploit both P‐ and mode‐converted PS‐waves. Whilst P‐wave events are exploited for inversion for P‐wave velocity, examples show the contribution of both P‐ and PS‐waves to the successful recovery of S‐wave velocity.     

Reflection and transmission of inhomogeneous waves in a composite porous solid saturated by two immiscible fluids

This paper is concerned with reflection and transmission of a plane, elastic, and inhomogeneous wave striking obliquely at some discontinuity inside a porous medium composed of two distinct solids and saturated by two immiscible fluids. It is found that four P‐ and two SV‐waves are reflected, whereas four P‐ and two SV‐waves are transmitted at the interface. All reflected and transmitted waves are inhomogeneous in nature and specified with different directions of propagation and attenuation vectors. An expression for the Umov–Poynting energy flux vector is derived for the system. Continuity of energy flux along normal to the interface gives 12 required boundary conditions. Expressions of amplitude and energy ratios of various reflected and transmitted waves are derived. Variations in amplitude and energy coefficients of reflected and transmitted waves with angle of incidence are numerically studied for a porous matrix composed of shaley sandstone and clay, saturated with water and oil. The effects of change in oil saturation and volume fraction of clay are also observed on amplitude ratios. Numerical simulation reveals that the change in sign in the difference of capillary pressure across the interface causes jump in the values of amplitude ratios of all waves.     

TTI介质各向异性参数多波反演与PS波AVO分析

把遗传算法引入到了TTI介质AVO信息反演各向异性参数的过程中,依据TTI介质PP波、PS波反射系数公式,建立Thomsen参数和TTI介质对称轴倾角、方位角的目标函数,分别通过PP波和PS波的反射系数反演出了各向异性参数和对称轴倾角、方位角等信息.文中对反演结果的精确度和稳定性进行了分析,发现PS波的反演结果优于PP波反演结果;对称轴倾角的反演准确性明显优于对称轴方位角.本文通过模型正演合理解释了这一现象的原因.最后,本文通过对PS波AVO梯度的研究,提出了利用PS波振幅定性分析TTI介质对称轴倾角的方法.     

地震转换波对叠前同时反演的影响研究

叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠.     

基于Zoeppritz偏导方程精确解的地层密度多角度反演

基于Zoeppritz方程对介质密度偏导数所建立的偏导方程的精确解,构造了多角度反演地层介质密度的反演方程,在偏导数求解过程中考虑了介质密度对波速度的影响因素,并由此实现了利用反射系数梯度精确解计算地层密度的多角度联合反演.通过数值算例考察了计算方法,结果显示:反演方法对层状地层模型不论反射波是否存在相干现象均获得了较好的反演结果,反演迭代10次后计算结果的最大相对误差能够收敛到1%之内;随着反演角度的增加地层介质密度反演的精度逐步提高,反演具有自动校正能力,有快的计算速度.本方法克服了传统AVO(Amplitude Versus Offset)基于Zoeppritz方程近似所遇到的困难,不受反演角度大小及反射界面对波反射强弱的限制,为地层介质密度的多角度包括大角度反演提供了一种新的快速有效的计算方法.     

基于Russell近似的纵横波联合反演方法研究

PP波和PS波联合反演方法作为有效的地震技术,比单纯纵波反演精度要高,能够提高地震储层识别的精度.以Russell近似理论为基础,推导了新的转换波AVO近似公式,双层模型界面的反射特征数值模拟显示,新公式具有较高的近似精度,且具备直接反演流体因子f、剪切模量μ和密度ρ等参数的优势,有效避免间接反演带来的误差.结合纵横波联合反演理论,提出了基于贝叶斯理论的新型联合反演算法.在实际应用中,对纵波和转换波角道集进行同相轴匹配处理,综合利用纵波和转换波资料携带的信息,实现基于Russell近似的多波联合反演.模型数据和实际资料测试结果表明,反演结果与真实值或测井结果匹配度较高,证实该方法真实有效.     

Joint PP‐ and PSV‐wave amplitudes versus offset and azimuth inversion for fracture compliances in horizontal transversely isotropic media

We propose a robust approach for the joint inversion of PP‐ and PSV‐wave angle gathers along different azimuths for the elastic properties of the homogeneous isotropic host rock and excess compliances due to the presence of fractures. Motivated by the expression of fluid content indicator in fractured reservoirs and the sensitivity of Lamé impedances to fluid type, we derive PP‐ and PSV‐wave reflection coefficients in terms of Lamé impedances, density, and fracture compliances for an interface separating two horizontal transversely isotropic media. Following a Bayesian framework, we construct an objective function that includes initial models. We employ the iteratively reweighted least‐squares algorithm to solve the inversion problem to estimate unknown parameters (i.e., Lamé impedances, density, and fracture compliances) from PP‐ and PSV‐wave angle gathers along different azimuths. Synthetic tests reveal that the unknown parameters estimated using the joint inversion approach match true values better than those estimated using a PP‐wave amplitude inversion only. A real data test indicates that reasonable results for subsurface fracture detection are obtained from the joint inversion approach.     

An Accurate Method of Computing the Gradient of Seismic Wave Reflection Coefficients (SWRCs) for the Inversion of Stratum Parameters

The optimization inversion method based on derivatives is an important inversion technique in seismic data processing, where the key problem is how to compute the Jacobian matrix. The computational precision of the Jacobian matrix directly influences the success of the optimization inversion method. Currently, most of the AVO (amplitude versus offset) inversions are based on approximate expressions for the Zoeppritz equations to obtain the derivatives of the seismic wave reflection coefficients (SWRCs) with respect to the stratum parameters. As a result, the computational precision and range of applications of these AVO inversions are restricted undesirably. In order to improve the computational precision and to extend the range of applications of AVO inversions, the partial derivative equations of the Zoeppritz equations are established, with respect to the ratios of wave velocities and medium densities. By solving the partial derivative equations of the Zoeppritz equations accurately, we obtained the partial derivative of SWRCs with respect to the ratios of seismic wave velocities and medium densities. With the help of the chain rule for derivatives, the gradient of the SWRCs can be accurately computed. To better understand the behavior of the gradient of the SWRCs, we plotted the partial derivative curves of the SWRCs, analyzed the characteristics of these curves, and gained some new insight into the derivatives. Because only a linear system of equations is solved in our method without adding any new restrictions, the new computational method has both high precision and a quick running speed; it is not only suitable for small incident angles and weak reflection seismic waves but also for large incident angles and strong reflection seismic waves. With the theoretical foundations established in the article, we can further study inversion problems for layered stratum structures and we can further improve the computational speed and precision of the inversions.     

Extended reflectivity method for modelling the propagation of diffusive–viscous wave in dip‐layered media

The reflectivity method plays an important role in seismic modelling. It has been used to model different types of waves propagating in elastic and anelastic media. The diffusive–viscous wave equation was proposed to investigate the relationship between frequency dependence of reflections and fluid saturation. It is also used to describe the attenuation property of seismic wave in a fluid‐saturated medium. The attenuation of diffusive–viscous wave is mainly characterised by the effective attenuation parameters in the equation. Thus, it is essential to obtain those parameters and further characterise the features of the diffusive–viscous wave. In this work, we use inversion method to obtain the effective attenuation parameters through quality factor to investigate the characteristics of diffusive–viscous wave by comparing with those of the viscoacoustic wave. Then, the reflection/transmission coefficients in a dip plane‐layered medium are studied through coordinate transform and plane‐wave theory. Consequently, the reflectivity method is extended to compute seismograms of diffusive–viscous wave in a dip plane multi‐layered medium. Finally, we present two models to simulate the propagation of diffusive–viscous wave in a dip plane multi‐layered medium by comparing the results with those in a viscoacoustic medium. The numerical results demonstrate the validity of our extension of reflectivity method to the diffusive–viscous medium. The numerical examples in both time domain and time–frequency domain show that the reflections from a dip plane interface have significant phase shift and amplitude change compared with the results of horizontal plane interface due to the differences in reflection/transmission coefficients. Moreover, the modelling results show strong attenuation and phase shift in the diffusive–viscous wave compared to those of the viscoacoustic wave.     

P-SV波AVO方法研究进展

为有效地利用AVO(amplitude versus offset)信息来反演岩性参数和预测油气储层，详细介绍了P—SV转换波AVO方法的研究现状和主要进展，并针对几种具有代表性的方法的基本思路、方法特点及参数反演等方面进行了对比和评述。同时，给出了两个不同形式的近似公式，并讨论了利用这两种近似公式进行AVO分析的思路和参数反演的方法，通过对比和讨论可以得出：不同近似的主要目的是为了提高近似精度并体现不同的岩性参数对反射系数的影响、敏感程度及其在参数反演中的意义；根据弹性模量和波速之间的关系，目前所有的PP波和P—SV波反射系数近似公式都可以统一表示成射线参数幂级数的形式，这些结论对于AVO理论研究和参数反演都具有重要意义。     

AN ATTEMPT AT THE INVERSION OF REFLECTION DATA*

The main problem in seismic prospecting is to infer from the observed reflection response the distribution of density and seismic velocity with depth. This process is generally called the inversion of the reflection data. For plane waves propagating through plane parallel stratification, it can be shown that at any depth the ratio between the amplitude of the transmitted and reflected wave satisfies the Riccati equation. Based on this equation we have formulated an iterative inversion method, which is found to be suitable for numerical computations. We have applied this method on synthetic reflection data, and found that it provides a very fast and accurate inversion.     

The form of the reflection operator and extended image gathers for a dipping plane reflector in two space dimensions

The reflection operator for a simple flat‐lying interface can be thought of as the set of all its plane‐wave reflection coefficients or as the set of virtual surveys with sources and receivers along the interface. When there is dip, however, it is necessary to include the varying effects of propagation between the virtual‐survey level and the interface. Hence, step one in this paper is to derive the reflection operator for a dipping plane interface as observed at a datum level some distance away. The key assumption is that the aperture at the datum level is sufficient to characterize the reflector properties around a particular point. This translates into an assumption that the dip is moderate, though no explicit small‐angle approximation is required. The second step is to find the apparent reflection operator that would relate data that have been extrapolated from the datum towards and possibly beyond the reflector using an assumed migration velocity. This apparent reflection operator is closely related to extended common‐image gathers. The apparent reflection operator may be analysed asymptotically in terms of rays and other signals, shedding light on the structure of extended image gathers. In keeping with the virtual‐survey idea, the results are considered in a subsurface space‐time or slowness‐time domain at various extrapolation levels around the interface. An important distinction is drawn between using subsurface midpoint‐offset coordinates and the wavefield coordinates of the incident and reflected waves. The latter reveal more clearly the effects of dip, because they lead to a more asymmetric apparent reflection operator. Properties such as an up‐dip shift of a traveltime minimum and its associated curvature theoretically provide information about the reflector location and dip and the migration‐velocity error. The space‐time form of the reflection operator can be highly intricate around the offset‐time origin and it was described for a simple flat interface in a background paper. To avoid a layer of mathematics, the reflection‐operator formulas presented here are in the intermediate space‐frequency domain. They are analysed by considering their stationary‐phase and branch‐point high‐frequency contributions. There is no Born‐like assumption of weak reflector contrast and so wide‐angle, total reflection and head‐wave effects are included. Snell’s law is an explicit part of the theory. It is hoped that the work will therefore be a step towards the goal of unifying amplitude‐versus‐offset, imaging and waveform inversion.     

声波在两种多孔介质界面上的反射和透射

本文导出了声波在两种多孔介质界面上反射、透射的一般计算公式.作为例子,数值计算了P₁波入射于界面时,P₁、P₂和S波的反射、透射系数与声波频率、入射角等量之间的关系.结果表明,各种模式波的反射、透射系数与入射角、多孔介质性质有关,在Biot特征频率附近与频率有关,并用界面两侧的法向能流相等验证了结果的正确性.若把多孔介质当作均匀固体处理,将会得到显著不同的结果.     

基于耦合反射/透射系数单程波传播算子的地震波模拟研究

单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合.     

AVO investigations of shallow marine sediments

Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m³, quality factor for P‐wave Q_P=15±200, quality factor for S‐wave Q_S=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m³,Q_P=790±660,Q_S=25±10.     

电场垂直入射面的非均匀电磁波在导电媒质界面的类全反射横向偏移

利用电磁波在导电媒质界面反射系数的附加相角，导出了电场垂直于入射面的非均匀电磁波在导电介质界面的类全反射横向偏移，并对横向偏移进行了相关计算，绘出了横向偏移随入射角的变化曲线. 结果表明：入射角在相移常数临界角、衰减常数临界角和90°处，横向偏移曲线存在三个间断点. 当入射角等于这三个角时，横向偏移为无穷大，即电磁波将沿界面传播；当入射角在这三个角附近时，横向偏移变得非常大.     

Multi-parameter nonlinear inversion with exact reflection coefficient equation

Amplitude variation with amplitude or angle (AVO/AVA) inversion has been widely utilized in exploration geophysics to estimate the formation of elastic parameters underground. However, conventional AVO/AVA inversion approaches are based on different approximate equations of Zoeppritz equations under various hypotheses, such as limited incident angles or weak property contrast, which reduces their prediction precision theoretically. This study combines the exact P-wave Zoeppritz equation with a nonlinear direct inversion algorithm to estimate the six parameters imbedded in the exact equation simultaneously. A more direct and explicit expression of the Zoeppritz equation is discussed in the case of P-wave exploration, under which condition the incident longitudinal wave produces the reflected longitudinal (P–P) wave and upgoing converted shear (P–SV) wave. Utilizing this equation as the forward solver, a nonlinear direct inversion method is introduced to implement the direct inversion of the six parameters including P-wave velocities, S-wave velocities, and densities in the upper and lower media around an interface, respectively. This nonlinear algorithm is able to estimate the inverse of the nonlinear function in terms of model parameters directly rather than in a conventional optimization way. Model tests illustrate that the nonlinear direct inversion method shows great potential to estimate multiple parameters with the exact Zoeppritz equation.