共炮检距二维弹性波地震剖面的偏移方法

共均匀各向同性介质的分界面上入射纵波时，除产生纵波反射Ｐ－Ｐ波外，还会产生转换横波的反射Ｐ－ＳＶ波。同样，当入射波为ＳＶ横波时，在分界面上会同时产生ＳＶ－ＳＶ反射波和ＳＶ－Ｐ反射转换波。对这种二维弹性波地震记录如何进行偏移是探讨的问题。为了能够处理实际地震数据，本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。     

共炮检距二维弹性波地震剖面的偏移方法

在均匀各向同性介质的分界面上入射纵波时，除产生纵波反射Ｐ－Ｐ波外，还会产生转换横波的反射Ｐ－ＳＶ波。同样，当入射波为ＳＶ横波时，在分界面上会同时产生ＳＶ－ＳＶ反射波和ＳＶ－Ｐ反射转换波，对这种二维弹性波地震记录如何进行偏移是正在探讨的问题。为了能够处理实际地震数据，本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。具体做法是，首先用弹性波方程将两分量的地震记录分解为Ｐ－Ｐ（或ＳＶ－ＳＶ）波和Ｐ－ＳＶ（或ＳＶ－Ｐ）波，然后用我们已导出的适合于处理共炮检距剖面的叠前偏移方法对它们进行处理。用模拟数据进行的试验说明方法正确，效果良好。     

多波高斯束叠前深度偏移速度分析(英文)

多分量地震资料叠前深度偏移技术可以对地下复杂地质构造进行更准确的成像,精确成像的前提是获取准确的纵横波偏移速度。本文采用高斯束偏移方法对多波地震数据进行偏移速度分析,首先分别给出纵波和转换波共偏移距域高斯束叠前深度偏移方法原理,在此基础上抽取纵波和转换波偏移距域共成像点道集;然后根据共成像点道集拉平准则,分别对纵波和横波速度进行更新;当两种波成像深度不一致时,对纵波和转换波成像剖面进行深度匹配,完成高精度的纵横波偏移速度分析。模型数据和实际资料试算表明,该方法是一种有效的多波偏移速度分析方法。     

基于一阶速度-胀缩-旋转方程的纵横波地震偏移角度域共成像点道集提取方法（英文）

角度域共成像点道集(angle-domain common-image gathers,ADCIGs)是多波AVA反演和偏移速度分析的基础数据。纵、横波坡印廷矢量的准确求取是获得多分量地震资料ADCIGs的前提。由常规速度-应力弹性波方程得到的坡印廷矢量代表了纵、横两种波的混合波场的传播方向,而非纯纵波或纯横波的传播方向。一阶速度-胀缩-旋转弹性波方程能够求取纯纵波或纯横波的坡印廷矢量,获得单一类型波的传播方向信息。本文基于一阶速度-胀缩-旋转弹性波方程逆时偏移实现纵、横波ADCIGs的求取,具体思路为:利用胀缩振动速度矢量和胀缩标量分别求取炮点和接收点波场中纯纵波的坡印廷矢量,利用剪切振动速度矢量和旋转矢量求取纯横波的坡印廷矢量,取炮点波场中各网格点纵波能量最大时刻的纵波坡印廷矢量为入射纵波传播方向矢量,取接收点波场中各网格点纵(横)波能量最大时刻的纵(横)波坡印廷矢量为反射纵波(反射转换横波)传播方向矢量,据此求取入射角并最终得到纵、横波角道集。模型试算表明:本文方法可以准确求取界面处入射波和反射波的传播方向矢量,得到地震波入射角,实现纵、横波角道集的高精度提取。     

垂向不均匀介质中波传播特点的研究

本文由波传播问题的变分原理导出了垂向不均匀介质中的弹性波波动方程及纵波和横波波动方程。在不同的假设条件下得到了纵、横波波动方程的几种简化形式,其中包括地震勘探中常用的变系数波动方程。利用有限单元法,求得了垂向不均匀介质中波动方程的数值解。通过数值结果对纵、横波之间的耦合程度以及介质的不均匀性对波的衰减作用进行了讨论。     

基于保幅波场分离的弹性波逆时偏移方法研究（英文）

欲实现基于弹性波方程的矢量波场逆时偏移纵、横波独立成像,必须在波场延拓过程中实现纵、横波场的分离,散度和旋度算子分离的纵、横波出现振幅与相位的畸变,导致输出成像结果的振幅失真。本文提出一种在弹性波场延拓过程中实现纵、横波保幅分离的方法,在传统的弹性波方程中加入纵波压力、纵波振动速度和横波振动速度方程,实现纵横波的矢量分解,再对分解后的矢量纵波和矢量横波做标量化合成得到保幅分离的纵、横波场,对保幅分离的纵、横波场应用成像条件,然后实现矢量波场逆时偏移的保幅纵横波成像。该方法可以保证分离后纵、横波的振幅与相位不变;同时,分解后的纵波压力和纵波振动速度可用于层间反射噪音压制和横波极性校正,提高多分量地震资料联合逆时偏移的纵、横波成像质量,从而实现保幅弹性波逆时偏移的目的,为叠前深度剖面应用于叠前反演工作奠定基础。     

用非零偏P-P波、P-SV波计算零偏反射纵波和反射横波

在P波入射时，地下同一点产生的P－P波和P－SV波是入射角、界面两侧纵波阻抗和横波阻抗的函数，并且两者是耦合的，利用两者的耦合关系，可在入射角确定的情况下提取与入射角无关的零偏纵波反射和横波反射。本文对此进行了理论推导，结论是：入射角为零时的纵波与横波反射系数是入射角非零时的P－P波和P－SV波反射系数的线性组合，进而可以用纵横波反射计算泊松比反向间面，把多波叠前反演转化到叠后进行，通过合成弹性参数反射剖面，直接反演弹性参数等，都很有意义。文中给出了有关的公式，并对具体实现方法和有关参数进行了讨论。     

利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

矢量分离纵横波场的弹性波逆时偏移

弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像.     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分,但该方法会因为网格离散化而产生数值频散,降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散,然而纵、横波速度取值差异较大,在弹性波有限差分正演模拟中,在满足纵波最大速度确定的稳定性条件下,浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散,提高数值模拟精度,本文首先采用优化差分网格系数降低数值频散,然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制,但由于弹性波场中纵、横波速度差异明显,横波波场频散明显强于纵波,为了压制横波波场的数值频散,往往需要选取较大的频散压制参数,但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法,对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正,可以有效压制数值频散,削弱纵波FCT产生的假象.通过理论分析和数值算例发现,本文方法能有效削弱弹性波场有限差分数值...     

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.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

角度域弹性波Kirchhoff叠前深度偏移速度分析方法

为提高地震成像结果的准确性并真实反映实际地震波场在介质中的传播特性,应该充分利用多分量地震数据的矢量特征进行弹性波成像,其中,最为棘手的问题是纵横波偏移速度场的确定,为此,本文提出了直接利用多分量地震数据进行弹性波角度域偏移速度分析的方法.基于空移成像条件的弹性波Kirchhoff偏移方程提取了弹性波局部偏移距域共成像...     

Sign-change correction for prestack migration of P-S converted wave reflections1

The polarization direction or 'sign’ of reflected converted P–S waves depends upon the angle of incidence of the incident P-wave. Sign reversal due to reversal of the angle of incidence is often encountered and is an impediment to P–S wave processing and imaging, because when P–S events or P-S migrated images with mixed signs are stacked, destructive interference occurs. We have created and demonstrated a means of correcting for this reversal. To do this, a P-wave angle of incidence is calculated for every point in the image space. This is done by calculating a P–S reflected waveform for every point, by extrapolating the reflected S-wavefield backwards from the receiver line, and then cross-correlating this waveform with the S-wave reflections observed at the receiver line. A multiplier, (sgn α) is assigned to each point in the image space, where α is the angle of incidence of the P-wave. The multiplier was applied to a set of prestack reverse time migration images derived from a cross-borehole physical elastic model data set. The improvement in the stacked image when the sign correction is applied is spectacular. The P-S image quality is comparable to, or better than, stacked migrated P-P images. The method appears to be applicable to all reflection modes and to all recording geometries.     

克希霍夫法VSP多波联合成像

VSP 转换波跟VSP 纵波或常规地面转换波相比,具有较高的分辨率和信噪比,但传统的VSP成像方法只利用了反射P波信息,而把转换波（反射S波、透射S波）以及透射P波当作影响成像质量的噪音.本文给出了一种VSP共炮点道集多分量地震资料克希霍夫法偏移成像的方法.本方法充分利用了多波（反射P波、反射S波、透射P波、透射S波）信息,根据转换点处四种波同时起跳,能量叠加最大的原理,从接收点分别用向绕射点延拓它们的能量,并将其叠加起来,求得的和最大的一点即反射点.通过模型试算和实际资料处理表明,此法成像精度高,信噪比高,且有利于改善剖面的频率特性.     

弹性波有限元逆时偏移技术研究

改善弹性波有限元逆时偏移剖面的质量是目前该方法面临的难题，为此，首先在有限元方程中，合理地加入边界力项，并找到准确的成像条件，以及对炮集资料用惠更斯（Huygens）原理消除直达波和面波干扰等措施，提高了复杂含油气结构的偏移剖面质量．本文重点研究单分量资料弹性波偏移的可靠性，利用单分量、抽炮集资料，模拟实际资料处理典型含油气构造，得到了高质量的自激自收和炮集资料的偏移剖面．     

弹性波克希霍夫积分偏移法

P、S波同时偏移法是为了适应复杂构造地区勘探情况而开始研究的一种偏移方法。在前人工作的基础上,应用弹性波克希霍夫积分偏移的基本原理,改进了纵波（PP）及转换波（PS）的同时克希霍夫积分偏移方法,给出了一种近于实用的计算机流程和程序;并用其分别对射线理论合成反射波资料、有限元模拟地震波传播资料进行了检验,证明此种方法是有效的。在解决复杂构造地区的地震波成象方面,该方法具有很大的潜力。     

2-D Frequency-domain Waveform Inversion of Coupled Acoustic-Elastic Media with an Irregular Interface

In order to correctly interpret marine exploration data, which contain many elastic signals such as S waves, surface waves and converted waves, we have developed both a frequency-domain modeling algorithm for acoustic-elastic coupled media with an irregular interface, and the corresponding waveform inversion algorithm. By applying the continuity condition between acoustic (fluid) and elastic (solid) media, wave propagation can be properly simulated throughout the coupled domain. The arbitrary interface is represented by tessellating square and triangular finite elements. Although the resulting complex impedance matrix generated by finite element methods for the acoustic-elastic coupled wave equation is asymmetric, we can exploit the usual back-propagation algorithm used in the frequency domain through modern sparse matrix technology. By running numerical experiments on a synthetic model, we demonstrate that our inversion algorithm can successfully recover P- and S-wave velocity and density models from marine exploration data (pressure data only).     

双相横向各向同性介质分界面上弹性波反射与透射问题研究

从双相横向各向同性介质弹性波波动方程出发 ,利用边界上的 4个连续性条件 ,计算双相横向各向同性介质分界面上弹性波反射和透射系数 .计算表明 ,快纵波在双相横向各向同性介质分界面上 ,要产生反射快纵波、反射转换 qSV波、反射转换慢纵波和透射快纵波、透射转换qSV波、透射转换慢纵波 .反射转换慢纵波振幅和透射转换慢纵波振幅均较小 .频率、耗散和各向异性大小影响着快纵波反射系数的大小 .