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.     

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

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

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

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

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

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

有效应力变化对致密砂岩孔隙结构及弹性波响应的影响规律

致密砂岩气藏具有裂缝发育和有效应力高的特征,研究不同有效压力下孔、裂隙介质地震波传播特征,有利于地震解释与地下储层的识别.但是前人的研究较少考虑岩石内部微观孔隙结构特征与孔隙、裂隙间流体流动的关系.本文首先通过选取四川盆地典型致密砂岩岩样,在不同有效压力下对岩石样本进行超声波实验测量.然后基于实验测得的纵、横波速度进行裂隙参数反演,得到不同有效压力下致密砂岩样本的裂隙孔隙度.再将裂隙孔隙度和样本岩石物理参数代入双重孔隙介质模型,模拟得到不同有效压力下饱水致密砂岩样本纵横波速度频散和衰减的变化规律.结果表明模型预测的速度频散曲线与纵波速度实验测量结果能够较好的吻合.最后统计分析了致密砂岩裂隙参数,得到了致密砂岩储层裂隙参数随有效压力及孔隙度变化特征.依据实际岩石物理参数建立模型,其裂隙参数三维拟合结果能够较好描述致密砂岩裂隙结构与孔隙度、应力的关联,可为实际地震勘探中预测储层裂缝性质提供基础依据.     

DECOMPOSITION OF MULTICOMPONENT SEISMIC DATA INTO PRIMARY P- AND S-WAVE RESPONSES1

Inversion of multicomponent seismic data can be subdivided in three main processes: (1) Surface-related preprocessing (decomposition of the multicomponent data into ‘primary’ P-and S-wave responses). (2) Prestack migration of the primary P- and S-wave responses, yielding the (angle-dependent) P-P, P-S, S-P and S-S reflectivity of the subsurface. (3) Target-related post-processing (transformation of the reflectivity into the rock and pore parameters in the target). This paper deals with the theoretical aspects of surface-related preprocessing. In a multicomponent seismic data set the P- and S-wave responses of the subsurface are distorted by two main causes: (1) The seismic vibrators always radiate a mixture of P- and S-waves into the subsurface. Similarly, the geophones always measure a mixture of P- and S-waves. (2) The free surface reflects any upgoing wave fully back into the subsurface. This gives rise to strong multiple reflections, including conversions. Therefore, surface-related preprocessing consists of two steps: (1)Decomposition of the multicomponent data (pseudo P- and S-wave responses) into true P- and S-wave responses. In practice this procedure involves (a) decomposition per common shot record of the particle velocity vector into scalar upgoing P- and S-waves, followed by (b) decomposition per common receiver record of the traction vector into scalar downgoing P- and S-waves. (2) Elimination of the surface-related multiple reflections and conversions. In this procedure the free surface is replaced by a reflection-free surface. The effect is that we obtain ‘primary’ P-and S-wave responses, that contain internal multiples only. An interesting aspect of the procedure is that no knowledge of the subsurface is required. In fact, the subsurface may have any degree of complexity. Both the decomposition step and the multiple elimination step are fully determined by the medium parameters at the free surface only. After surface-related preprocessing, the scalar P- and S-wave responses can be further processed independently by existing scalar algorithms.     

小近震三分向地震图S波记录振幅包络形状及Q值测定

本文采用三层地壳介质模型,讨论了地方台网记录的小近震三分向地震图S波段振幅包络形状问题.认为S波包络的规则形状主要是由于强反射界面上的一次和多次全反射波形成的.分析了地壳分层结构,震源辐射特性,介质衰减对n次全反射波振幅的影响.在此基础上,从理论上给出了从直达S波到尾波之间的S波段记录部分振幅包络形状函数和S波最大振幅.利用包络形状函数方法测定了宁夏地区介质品质因素Q值.      

Linearized elastic parameter sections1

Contrasts in elastic parameters can be estimated directly from seismic data using offset-dependent information in the PP reflection coefficient. A linear approximation to the PP reflection coefficient including three coefficients is fitted to the data, and relative contrasts in various elastic parameters are obtained from linear combinations of the estimated coefficients. Linearized elastic parameter sections for the contrasts in P-wave impedance, P-wave velocity, density, plane-wave modulus, and the change in bulk modulus and shear modulus normalized with the plane-wave modulus are estimated. If the average P- to S-wave velocity ratio is known, linearized parameter sections including the contrast in the average P- to S-wave velocity ratio and a fluid factor section can be computed. Applied to synthetic data, visual comparison of the estimated and true elastic parameter sections agree qualitatively, and the results are confirmed by an analysis of the standard deviation of the estimated parameters. The parameter sections obtained by inversion of a shallow seismic anomaly in the Barents Sea are promising, but the reliability is uncertain because neither well data nor regional trends are available.     

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

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

Pure P- and S-wave equations in transversely isotropic media

Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media.     

基于非饱和Biot-Stoll模型的海底沉积物介质声频散特性研究

利用含气非饱和Biot-Stoll模型研究了声波在海底表层沉积物介质中的传播,讨论了骨架耗散、含气饱和度对快纵波、慢纵波和横波速度和衰减的影响,并与Biot模型的结果进行了对比.研究结果表明:孔隙流体黏滞耗散与骨架耗散共同影响声波传播速度和衰减,低频情况下骨架耗散引起的衰减占主要地位,高频情况下骨架耗散引起的衰减较小;少量气体(<1%)的引入显著改变了快纵波速度,气体含量的变化对快纵波衰减影响很大,低频情况下气体对慢纵波速度的影响不大,而对横波速度的影响较大,气体含量的变化对慢纵波和横波衰减影响较小.利用超声波测量系统测量了一例杭州湾海底沉积物样品的纵波速度和衰减,当含气量趋近0%时,Biot-Stoll模型预测的纵波速度和实验测量结果较为一致.     

应用三分量浅层地震反射方法探测隐伏活动断裂

多波多分量地震勘探方法在油气田勘探开发中日益受到重视,但在工程和活动断裂探测中过去还没有开展多波多分量地震勘探的先例.本文简要介绍了在北京黄庄—高丽营断裂上开展的三分量浅层地震反射试验数据采集和资料处理方法,以及取得的初步结果.试验结果表明,采用纵波震源或横波震源激发、三分量接收的浅层地震勘探方法,可获得较高信噪比的浅层地震纵波、横波和转换波剖面,特别是水平分量采集的横波和转换波可提供很有意义的地下结构和构造信息.综合利用纵波、横波和转换波剖面特征进行的地质分层和断层解释结果取得了与跨断层的钻孔联合地质剖面结果较好的一致性.     

Ultrasonic wave propagation in dry and brine-saturated sandstones as a function of effective stress: laboratory measurements and modelling1

Compressional and shear-wave velocities have been measured and a novel approach using digital processing employed to study wave attenuation for brine- and gas- saturated sandstones, over a range of effective stresses from 5 to 60 MPa. Also measured were the complex conductivity in the brine-saturated state and permeability in the gas-saturated state over the same range of stresses as for the velocity measurements. Broadband ultrasonic pulses of P- and orthogonally polarized S-waves in the frequency range 0.3–0.8 MHz are transmitted through the specimen to be characterized for comparison with a reference (aluminium) having low attenuation. The attenuation is calculated in terms of the quality factor Q from the Fourier spectral ratios, using the frequency spectral ratios technique. The corrections necessary for the effects of diffraction due to the finite size of the ultrasonic transducers have been carried out for the case of measurements under lower confining stress. To interpret the laboratory measured velocity and attenuation data under the physical conditions of this study and to estimate the effects of pore structure, numerical modelling of velocities and attenuation as functions of the confining stress have been performed, based on the MIT model. Theoretical models based on several hypothesized attenuation mechanisms are considered in relation to laboratory data of the effects of confining pressure, fluid saturation and pore structure on attenuation. Numerical calculations using these models with the experimental data indicate that friction on thin cracks and grain boundaries is the dominant attenuation mechanism for dry and brine-saturated sandstones at low effective stresses for the frequencies tested. However, for brine-saturated sandstones at moderately high effective stresses, fluid flow could play a more important role in ultrasonic S-wave attenuation, depending on the pore structure of the sample.     

南海北部神狐海域天然气水合物分解的测井异常

南海北部神狐海域GMGS-1钻探揭示SH3井天然气水合物位于稳定带上部,厚度约为10 m.氯离子异常计算的水合物饱和度最高达26%,高水合物饱和度层出现高电阻率和低纵波速度.为分析该低纵波速度异常,本文基于简化的三相介质理论计算了饱和水纵波速度,在深度195 m附近,测量的纵波速度小于饱和水纵波速度.利用阿尔奇公式,基于原位温度、盐度、密度孔隙度和测量的电阻率,利用交会分析确定了该井的阿尔奇常数为a=1.1和m=2.3.基于该参数,利用阿尔奇方程计算的水合物饱和度占孔隙空间5%~20%,局部地层水合物饱和度达26.8%,在垂向上分布不均匀.由于钻探可能导致水合物发生分解而产生游离气,原位游离气和水合物分解产生的气体都能造成低纵波速度异常.由于地震资料采集在测井之前完成,利用不同速度制作合成地震记录并与地震资料进行对比,能够确定水合物稳定带上部的低速异常形成原因.     

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.     

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.     

三维地震与OBS联合勘探揭示的神狐海域含水合物地层声波速度特征

以三维高分辨地震与海底高频地震仪（OBS）联合勘探数据为基础,获得海底之下沉积层的地震反射成像剖面及多波信息,并以此确定研究区含天然气水合物沉积层的纵、横波速度的变化特征.根据走时反演获得的横波速度与纵波速度对比分析发现,研究区海底之下500 m深度范围内的某些沉积层具有较高的纵横波速度,这一纵波速度升高区域与水合物稳定带对应,而纵波速度下降并且横波速度变化较小的区域,可能与游离气的存在相关.游离气的可能存在与基于这一区域2007年钻探测井结果的普遍认识不完全相符.     

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

以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量...     

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

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