Separation of P- and S-wavefields from wide-angle multicomponent OBC data for a basalt model

Wide-angle multicomponent ocean-bottom cable (OBC) data should further enhance sub-basalt imaging by using both compressional and converted shear wavefields. The first step in analysing multicomponent OBC data is to decompose the recorded wavefields into pure P- and pure S-wavefields, and extract the upgoing P- and S-waves. This paper presents a new scheme to separate P- and S-wavefields from wide-angle multicomponent OBC data in the τ–p domain. By considering plane-wave components with a known horizontal slowness, the P- and S-wavefields are separated into the directions of observed P- and S-wave oscillations using the horizontal and vertical components of the data. The upgoing P- and S-waves are then extracted from the separated P- and S-wavefields. The parameters used in the separation are the seismic wave velocities and the density at the receiver location, which can be estimated from the first reflection phase observed on the horizontal and vertical components. Numerical tests on synthetic data for a plane-layered model show good performance and demonstrate the accuracy of the scheme. Separation of wavefields from a basalt model is performed using synthetic wide-angle multicomponent OBC data. The results show that both near-offset and wide-angle reflections and conversions from within and below basalt layers are enhanced and clearly identified on the separated wavefields.     

Elastic reverse-time migration based on amplitude-preserving P- and S-wave separation

Imaging the PP- and PS-wave for the elastic vector wave reverse-time migration requires separating the P- and S-waves during the wave field extrapolation. The amplitude and phase of the P- and S-waves are distorted when divergence and curl operators are used to separate the P- and S-waves. We present a P- and S-wave amplitude-preserving separation algorithm for the elastic wavefield extrapolation. First, we add the P-wave pressure and P-wave vibration velocity equation to the conventional elastic wave equation to decompose the P- and S-wave vectors. Then, we synthesize the scalar P- and S-wave from the vector Pand S-wave to obtain the scalar P- and S-wave. The amplitude-preserved separated P- and S-waves are imaged based on the vector wave reverse-time migration (RTM). This method ensures that the amplitude and phase of the separated P- and S-wave remain unchanged compared with the divergence and curl operators. In addition, after decomposition, the P-wave pressure and vibration velocity can be used to suppress the interlayer reflection noise and to correct the S-wave polarity. This improves the image quality of P- and S-wave in multicomponent seismic data and the true-amplitude elastic reverse time migration used in prestack inversion.     

基于纵横波解耦的三维弹性波逆时偏移

纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性.

     

低信噪比海底多分量地震数据波场分离方法

以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量地震资料都无法基于现有的海底波场分离方法与流程取得合理的结果.本文以海底声波场与弹性波场分离基本原理为基础,通过对方法流程的修正,摆脱常规流程对中小偏移距直达波信号的依赖性.借助模拟数据实验讨论了波场分离对海底介质参数、噪声的敏感性.结合东海YQ探区海底多分量地震资料上-下行P/S波分离及其叠前深度偏移处理,验证了本文方法流程的可行性.

     

ELASTIC EXTRAPOLATION OF PRIMARY SEISMIC P- AND S-WAVES1

The elastic Kirchhoff-Helmholtz integral expresses the components of the monochromatic displacement vector at any point A in terms of the displacement field and the stress field at any closed surface surrounding A. By introducing Green's functions for P- and S-waves, the elastic Kirchhoff-Helmholtz integral is modified such that it expresses either the P-wave or the S-wave at A in terms of the elastic wavefield at the closed surface. This modified elastic Kirchhoff-Helmholtz integral is transformed into one-way elastic Rayleigh-type integrals for forward extrapolation of downgoing and upgoing P- and S-waves. We also derive one-way elastic Rayleigh-type integrals for inverse extrapolation of downgoing and upgoing P- and S-waves. The one-way elastic extrapolation operators derived in this paper are the basis for a new prestack migration scheme for elastic data.     

Multiwave velocity analysis based on Gaussian beam prestack depth migration

Prestack depth migration of multicomponent seismic data improves the imaging accuracy of subsurface complex geological structures. An accurate velocity field is critical to accurate imaging. Gaussian beam migration was used to perform multicomponent migration velocity analysis of PP- and PS-waves. First, PP- and PS-wave Gaussian beam prestack depth migration algorithms that operate on common-offset gathers are presented to extract offset-domain common-image gathers of PP- and PS-waves. Second, based on the residual moveout equation, the migration velocity fields of P- and S-waves are updated. Depth matching is used to ensure that the depth of the target layers in the PP- and PS-wave migration profiles are consistent, and high-precision P- and S-wave velocities are obtained. Finally, synthetic and field seismic data suggest that the method can be used effectively in multiwave migration velocity analysis.     

陆上纵波源VSP资料中的纯横波

常规陆上VSP(Vertical Seismic Profiling)勘探普遍采用纵波震源激发,三分量检波器接收,主要利用的是纵波和转换横波信息。已有的研究表明,炸药震源在井下激发、可控震源在地面垂向振动,均会产生较强的纯纵波和一定强度的纯横波;泊松比差别较大的分界面有利于形成较强的透射转换横波。本文通过对激发形成的纯横波和下行转换形成的横波进行对比分析,认为纯横波的主频往往低于纯纵波的主频,而下行转换横波的主频通常接近纵波的主频。本文分别对两个陆上纵波源零偏和非零偏VSP资料进行分析,结果表明这些资料中普遍存在纯横波,只是横波的强弱存在不同程度的变化。利用纵波源零偏VSP资料,可以获得横波速度信。最后对VSP纵波和横波联合应用前景进行了分析,应该充分利用纵波源VSP资料中的横波信息。     

海底天然气水合物OBS多分量地震正演模拟（英文）

根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。     

基于快速非支配排序遗传算法的VTI介质多分量叠前联合反演

在VTI介质中,由于引入了各向异性参数使得多分量多参数地震反演问题的非线性程度显著增加,因此采用传统的权重加权法构建单目标函数进行反演得到的反演结果往往并不理想.本文以反射率法为基础,结合快速非支配排序遗传算法研究了一种VTI介质的多分量叠前联合反演方法.该方法以反射率法为正演方程,应用互相关原理构建PP波和PSV波的多目标函数,进而采用快速非支配排序遗传算法全局寻优获得VTI介质的厚度、纵横波速度、密度和各向异性参数等多个参数.在正演的过程中,反射率法可以考虑几何扩散、吸收衰减、透射损失、多次波以及纵横波旅行时不匹配等地震波传播效应,更能精确地描述地震波在地下地层中的真实传播情况;在反演的过程中,快速非支配排序遗传算法可以在不引入权重系数的条件下同时优化多个目标函数,获得联合反演问题的Pareto最优解,既不添加权重系数影响又充分利用多分量地震数据.模型测试结果验证了该反演方法的有效性和可行性.

     

Exact elastic impedance tensor for isotropic media

The existing expressions of elastic impedance,as the generalized form of acoustic impedance,represent the resistance of subsurface media to seismic waves of non-normal incidence,and thus include information on the shear-wave velocity.In this sense,conventional elastic impedance is an attribute of the seismic reflection and not an intrinsic physical property of the subsurface media.The derivation of these expressions shares the approximations made for reflectivity,such as weak impedance contrast andisotropic or weakly anisotropic media,which limits the accuracy of reflectivity reconstruction and seismic inversion.In this paper,we derive exact elastic impedance tensors of seismic P-and S-waves for isotropic media based on the stress-velocity law.Each componentof the impedance tensor represents a unique mechanical property of the medium.Approximations of P-wave elastic impedance tensor components are discussed for seismic inversion and interpretation.Application to synthetic data and real data shows the accuracy and robust interpretation capability of the derived elastic impedance in lithology characterizations.     

2D and 3D amplitude-preserving elastic reverse time migration based on the vector-decomposed P- and S-wave records

The elastic reverse time migration approach based on the vector-wavefield decomposition generally uses the scalar product imaging condition to image the multicomponent seismic data. However, the resulting images contain the crosstalk artefacts and the polarity reversal problems, which are caused by the nonphysical wave modes and the angle-dependent reduction of image amplitudes, respectively. To overcome these two problems, we develop an amplitude-preserving elastic reverse time migration approach based on the vector-decomposed P- and S-wave seismic records. This approach includes two key points. The first is that we employ the vector-decomposed P- and S-wave multicomponent records to independently reconstruct the PP and PS reflection images to mitigate the crosstalk artefacts. The second is that we propose two schemes in addressing the issue of polarity reversal problem in the conventional PP image. One solution is to adopt the angle-dependent equation. Another one is to reconstruct an amplitude-preserving PP image with a separated scalar P-wave particle velocity, which has a clear physical meaning. Numerical examples using two-dimensional and three-dimensional models demonstrate that the proposed elastic reverse time migration approach can provide the images with better amplitude-preserving performance and fewer crosstalk artefacts, compared with the conventional elastic reverse time migration approach based on the scalar product imaging condition.     

基于改进反射系数近似方程的纵横波阻抗同步反演

基于地震波反射系数近似公式的叠前反演是油气勘探的重要工具.本文在已有研究的基础上,推导了一个改进的射线参数域地震纵波反射系数近似方程.该方程建立了地震纵波反射系数与纵波阻抗和横波阻抗的非线性关系,在中、小角度的范围内较现有的反射系数线性近似公式精度更高.另外,由于该方程仅包含纵波和横波阻抗反射系数项,因此基于新方程的反演能够有效地降低同步反演纵波速度、横波速度、密度三个参数的不适定性.在此基础上,结合广义线性反演法（GLI）理论和贝叶斯理论,相应地发展了一种叠前地震同步反演方法.模型测试和实际资料的应用表明,基于新方程的反演方法能够利用有限角度（偏移距）的数据稳定地反演纵波和横波阻抗,由于在反演过程中,不需要假设纵横波速度为常数,因此该方法还能有效地提高反演结果的精度.

     

SEPARATION OF P- AND S-WAVES*

In combined P- and SV-observations with vertical and in-line horizontal geophones, recognition and separation of the individual P- and S-arrivals may form a considerable problem. Due to the geophone receiving characteristics, P- and S-waves are observed on both horizontal and vertical geophones. A new method of separating P- and S-waves is presented that inverts the geophone records for the receiving characteristics, thereby separating the two wavetypes. The filter coefficients are determined from the near-surface P- and S-wave velocity and from the geometry of the geophone groups. The P/S-separation filter was tested on synthetic data with good results. Even PS- and SP-converted reflections, which are completely coinciding for horizontal reflectors, could be well separated. The filter appears to be stable and to have low sensitivity to noise in the data or to errors in the near-surface P- and S-wave velocity. The filter is superior to other separation methods, which were applied to the same dataset for comparison. This method can also be applied to VSP data.     

Extraction of amplitude-preserving angle gathers based on vector wavefield reverse-time migration

Angle-domain common-image gathers (ADCIGs) transformed from the shotdomain common-offset gathers are input to migration velocity analysis (MVA) and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reversetime migration for calculating the incidence angle of P-and S-waves at each image point and for different source locations. The P-and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors’ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration (RTM) images. The ADCIGs of P-and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P-and S-wave ADCIGs of the elastic-wave reverse-time migration, the P-and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion.     

Detection of gas and water using HHT by analyzing P- and S-wave attenuation in tight sandstone gas reservoirs

A direct detection of hydrocarbons is used by connecting increased attenuation of seismic waves with oil and gas fields. This study analyzes the seismic attenuation of P- and S-waves in one tight sandstone gas reservoir and attempts to give the quantitative distinguishing results of gas and water by the characteristics of the seismic attenuation of P- and S-waves. The Hilbert–Huang Transform (HHT) is used to better measure attenuation associated with gas saturation. A formation absorption section is defined to compute the values of attenuation using the common frequency sections obtained by the HHT method. Values of attenuation have been extracted from three seismic sections intersecting three different wells: one gas-saturated well, one fully water-saturated well, and one gas- and water- saturated well. For the seismic data from the Sulige gas field located in northwest Ordos Basin, China, we observed that in the gas-saturated media the S-wave attenuation was very low and much lower than the P-wave attenuation. In the fully water-saturated media the S-wave attenuation was higher than the P-wave attenuation. We suggest that the joint application of P- and S-wave attenuation can improve the direct detection between gas and water in seismic sections. This study is hoped to be useful in seismic exploration as an aid for distinguishing gas and water from gas- and water-bearing formations.     

STATIC CORRECTIONS FOR SHEAR WAVE SECTIONS*

The computation of static corrections requires information about subsurface velocities. This information can be obtained by different methods: surface wave analysis, short refraction lines, downhole times, uphole times and first arrivals from seismograms. For pure shear waves generated by SH sources the analysis of first arrivals from seismograms combined, if necessary, with short refraction lines has proved to be most accurate and economic. A comparison of first-arrival plots from P- and S-wave surveys of the same line measured in areas of unconsolidated sediments in northern Germany illustrates the characteristic differences between the two velocity models. P-waves show a marked velocity increase at the water table from about 600 to 1800 m/s. S-wave velocities of the same strata increase gradually from about 100 to 400 m/s. As a consequence, S-wave models are vertically and laterally more complex and, in general, show no significant velocity increase at a defined boundary as P-wave models do. Therefore, other suitable correction levels with specific velocities must be chosen. A comparison of “tgd-corrections” (correction time between geophone position and datum level) for P- and S-waves in areas of unconsolidated sediments shows that their ratio is different from the P-/S-velocity ratio for the respective correction level because of the greater depth of the S-wave refractor. Therefore, P- and S-waves are influenced by different near-surface anomalies, and time corrections calculated for both wave types are largely independent.     

基于散度和旋度纵横波分离方法的改进

纵、横波的分离是多波多分量地震资料处理中很重要的一步,其分离结果直接影响到后续数据处理的质量.各向同性介质中纵波为无旋场,横波为无散场,因此可以在频率-波数域利用散度和旋度算子对地震记录进行纵、横波分离,但是此处理过程必须知道地表处的纵、横波速度.本文给出了一种估算地表纵、横波速度的方法,可以在纵、横波速度值未知的情况下,将其估算出来.针对弹性波场进行散度和旋度运算时,纵、横波的相位和振幅比发生改变的问题,本文给出了相位和纵、横波振幅比的校正方法.     

地震各向异性流体检测技术在胜利油田垦71井区的应用(英文)

根据Chapman理论模型,在各向异性介质(如HTI介质)中,当入射角在0-45。范围内,慢横波会发生较大的衰减和频散,且对流体粘度敏感,而P波和快横波则比较小。对于沿裂隙法向传播的慢横波,其振幅受流体影响很大。因此,在P波响应对流体不敏感的情况下,可利用慢横波来获得裂隙型油气藏的流体信息。本文分析了胜利油田垦71地区三维三分量地震数据,检测出的慢横波振幅和旅行时异常与该区的测井资料十分吻合。分析结果还发现,与含油区相比,含水区会产生更高的横波分裂。在含水区,慢横波振幅会产生明显变化,而在含油区则几乎没有变化。     

区域地震信号自动识别方法及应用(英文)

地震信号的实时、自动、准确识别对于地震自动速报和地震预警十分重要。仿真信号试验分析表明,观测数据的四阶统计量函数(BKCF)对信号与噪声在能量和(或)频率方面的微弱差异变化具有较高的分辨能力。以此为基础,本文提出了一种新的自动探测区域地震事件的方法和测定直达波震相到时的BKCF-AIC方法。为了进一步提高波震相到时测定的精度,本文首先对指定时段的P-波记录进行偏振特性分析,其次对含有P波的S波记录进行偏振滤波处理,再次应用上述方法测定震相到时。与传统算法相比,基于山东测震台网记录的区域地震震例分析结果表明,使用本文提出的方法能够大幅度降低地震事件误检、漏检率,进一步提高了震相识别精度。     

高铁地震震源子波时间函数及验证

基于纵横波产生机理并结合高铁列车及铁轨结构,本文给出了高铁列车产生地震纵波及横波震源（点源及等效源）子波的时间函数;通过实际高铁地震数据及合成高铁地震数据分析,验证了高铁地震震源子波时间函数的合理性;在高铁列车匀速、加速及减速运行情况下,利用高铁地震点源子波合成记录分析了高铁地震震源子波的频谱特征,为实际高铁地震数据处理解释奠定了良好的基础.