基于矢量波场逆时偏移的保幅角道集提取方法（英文）

将偏移后的炮域偏移距道集转换为角度域共成像道集(ADCIGs)可为偏移速度分析(MVA)和叠前反演提供输入道集,并且ADCIGs是理论上没有假象的叠前反演道集,也是目前公认的精度最高的叠前反演道集。本文研究了基于矢量波场逆时偏移的弹性波保幅ADCIGs的提取方法,以保幅弹性波逆时偏移方程为基础,其核心是求取不同震源位置的纵、横波场在地下各成像点的入射角,对于转换波勘探,二者共享一个入射角,即震源纵波入射角。根据几何关系,震源纵波波场的传播角、构造的局部地层倾角之差为震源纵波入射角,震源纵波波场的传播角利用解耦后纵波场的极化向量得到,构造的局部地层倾角利用偏移叠加剖面的复波数得到。对纵、横波的共炮点偏移道集按入射角重新排列即可得到各自的ADCIGs。文中利用水平层状介质模型、倾斜层状介质模型、Marmousi-Ⅱ弹性波部分模型和实测资料验证了算法的有效性,计算结果表明,本文方法计算的纵、横波角度具有较高的精度,提取的角道集具有较好的振幅保真性,能够为MVA和叠前反演提供可靠的输入道集。     

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

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

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

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

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

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

裸眼井中软地层多极子源激发的弹性波的传播

利用全波场和分波场算法对裸眼井软地层中由单极子源和多极子源（偶极子源和四极子源）激发的弹性波进行了数值研究,给出了声波全波列理论波形、纯纵波和纯横渡理论波形及纵横波的频率传输特性曲线.结果表明,纯纵波和纯横波是无频散几何衰减波;除了在v_s>（1-ρ₁/ρ₂）^1/2v₁的井眼中由单极子源激发的横波无传播特征频率外,纵横波都具有各自的传播特征频率,纵波的传播特征频率有一系列离散值,横波的传播特征频率只有一个且等于非对称模式波的截止频率;多极子源激发的纵波传播特征频率的最小值大于横波的传播特征频率,随着极子数增多,二者的特征频率值皆增高,差别也越大,因此,当声源发射主频较低,接近横波传播特征频率而远离纵波传播特征频率时,低频多极子源将压制纵波并激发出较强的横波,分析结果还表明,通常所说的低频多极子源激发出的横波波群实际上是纯横波与截止频率附近的非对称模式波叠加的结果,在其它条件不变时,二老对横波波群相对贡献的大小取决于声源激发主频的高低及其带宽的大小.     

基于能量密度的自解耦互相关成像条件

基于能量密度构建的弹性波能量成像条件,可以准确提取弹性波场的能量信息,并可以有效压制背向散射.但是该成像条件得到的成像结果将所有弹性波场能量信息耦合在一起,难以区分纯波模式的能量信息.为此,我们从势能及动能两方面将能量场解耦,得到纯纵波（PP）、纯横波（SS）、转换波（C）能量场.根据能量守恒定律,基于Helmholtz原理对位移场的空间导数进行因式分解,将弹性势能自解耦得到纯纵波、纯横波及转换波势能;引入体应力构建得到P波及S波速度,将弹性动能分解为纯纵波、纯横波和转换波动能.从而,构建基于能量密度自解耦得到纯波能量互相关成像条件.数值实验表明自解耦成像条件可实现弹性波场能量的分解,并得到背向散射压制、振幅有效保持的PP波、SS波和C波成像结果.     

联合纵波确定转换波接收点大的静校正量

在转换波地震勘探数据处理中,接收点的横波静校正一直是较难解决的问题.因为横波速度基本不受潜水面的影响,横波低速带厚且不均匀,速度又比较低,通常会造成转换波接收点的横波静校正量大且横向变化剧烈,与纵波静校正没有直接相关性.针对这种情况,我们采用改进的纵波构造时间控制方法求取转换波接收点的静校正,在转换波的共接收点叠加剖面上追踪拾取多组同相轴,联合纵波对应层位的地下构造时间确定接收点大的横波静校正量.通过不同的时移将几个相邻的层位合成一个反射层位,解决同相轴不能连续追踪的难题,同时利用浅、中和深层多组层位对一个接收点的静校正量求其平均,降低拾取误差对静校正量的影响,兼顾地下浅、中和深层的静校正效果.实际数据处理结果表明,该静校正方法能稳定的解决转换波静校正问题,取得了较好的应用效果.     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.     

部分叠加资料纵横波速度测井约束反演

本文提出一种角度部分叠加资料同时反演纵波和横波速度的方法,以角度部分叠加资料为基础,利用地震波振幅随入射角变化与弹性参数间的数学关系,基于非线性最优化理论,通过将最小平方问题转化为大型带状矩阵的求解问题,采用逐道外推技术依次求得该角度剖面的每个点纵波和横波速度值,利用地震波速度在各向异性介质中与入射角的关系,进一步得到地震波零入射角纵、横波速度,从而可以求得纵横波速度比及泊松比,为地震资料的岩性及含油气性解释提供了丰富的参数信息,为保证反演的稳定性,降低多解性,采用测井、叠后构造层位解释结果的共同约束,并做好叠前去噪以及振幅保持处理.     

疏松砂岩气藏地震AVO属性敏感性分析与评价

利用完全Zoeppritz方程分析了疏松砂岩气藏反射系数对密度比、纵波速度比、横波速度比、泊松比和含气饱和度的敏感性.通过求取反射系数对各弹性参数的偏导得到反射系数对弹性参数的敏感性,通过求取反射系数对含气饱和度的完全导数得到对含气饱和度的敏感性.对于疏松砂岩气藏界面,P-P波AVO(振幅随偏移距变化)当入射角大于40°对纵波速度比、横波速度比在低含气饱和度时敏感,对泊松比的敏感性在气层饱和度较高和较低时差异明显,对含气饱和度的敏感性在低饱和度时敏感,且随入射角增加而增大;P-SV波AVO对泊松比的敏感性在含气饱和度较高和较低时差异明显.研究表明,利用反射系数对弹性参数和含气饱和度的敏感性分析可以帮助识别低饱和气藏.     

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

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

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

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

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

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

基于应力偏量的横波纯应力逆时偏移成像方法

基于弹性波解耦延拓方程的波场分离方法不仅可以得到解耦的纵、横波质点振动速度场,还可以获得纵、横波应力场.针对横波纯应力场在利用单一分量成像时不具有明确物理意义的问题,本文将应力偏张量引入到横波应力场中,基于应力偏量第二不变量构建得到横波应力不变量并将其用于应力场的逆时偏移成像中,获得了可完整表征横波应力场的成像结果.模型试算表明,本文构建的横波应力不变量可以有效利用横波应力张量中的波场信息,并得到准确的弹性逆时偏移成像结果.     

厚风化层覆盖区转换波静校正方法

P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.     

共转换点道集S变换谱相似性的多波地震数据横波波场提取技术

横波速度动校正后的共转换点(CCP)道集内,同时刻的各道横波信号S变换(ST)谱与其叠加道ST谱具有相似关系.因此,可基于这种相似关系设计自适应滤波器来提取多波地震数据中的横波波场.首先对共中心点(CMP)道集应用纵波速度动校正并在各道减去叠加道来去除数据中的纵波波场;然后在CCP道集应用横波速度动校正,将地震道振幅水平调整至叠加道振幅水平并做S变换,以叠加道ST谱为参考对地震道ST谱进行自适应滤波,去除数据中的残余纵波和噪声;最后,将滤波结果的振幅水平恢复至滤波前振幅水平.理论和实际数据试算表明,本文方法可有效提取多波地震数据中的横波波场,为多波多分量横波数据处理提供新思路.     

求解弹性波解耦方程的一种优化拟解析方法

常规伪谱方法二阶时间差分格式时间精度较低,且对于大步长时间采样间隔,常规伪谱方法不稳定.拟解析方法对于速度变化剧烈的模型,在时间和空间上均有较大误差.本文提出了一种基于解耦的二阶位移弹性波方程波场模拟及矢量波场分解的优化拟解析方法,将归一化的拟拉普拉斯算子分别应用于P波和S波波场延拓,延拓矢量波场的同时,可分解并延拓纯纵波和纯横波波场.利用弹性波优化拟微分算子表示拟拉普拉斯算子,该拟微分算子不仅包括原始微分算子的谱估计,而且还包含一个时间补偿项,其可在波数-空间域精确地补偿波动方程在时间方向上采用二阶有限差分引起的误差.利用低秩分解近似求解弹性波优化拟微分算子,可有效提高计算效率.2D均匀模型、层状模型以及部分盐丘模型数值正演模拟结果表明:相比较于常规的伪谱法和拟解析法,本文方法在时间和空间上均有很高的精度,并且稳定性条件比较宽松.     

琼东南海域水合物OBS资料成像处理关键技术

海底地震仪(OBS)资料同时具备纵波和转换横波的优势,但其接收站点少、接收间距大的特点导致其处理方式相对常规海洋地震资料而言具有极大的挑战.对本文主要针对于OBS资料处理的关键技术:利用直达波进行OBS重定位处理,保证后续处理的正确性;利用水检资料和陆检资料进行波场分离,获得上行波场和下行波场;利用下行波进行P波镜像叠前时间偏移处理;对转换横波利用欧拉重构进行重定向处理及极化旋转处理、进行转换横波叠前时间偏移处理.最终获得海底成像清楚、海底以下成像分辨率高、偏移归位准确的OBS纵剖面及相应的转换横波剖面.同时具备的纵、横波剖面利于解释专家对水合物成矿区域进行综合分析.     

层状弹性介质的波阻抗反演

讨论层状弹性介质的纵波阻抗和横波阻抗的同时反演问题，由纵波线源激发产生的弹性波场，通过倾斜叠加可以得到其平面波分解，分解出来的平面波场对应于以不同倾角斜入射激发所产生的波场，它们满足由二维弹性波基本方程派生出的一维偏微分方程组。利用在地面同时接收的位移速度的水平和垂直分量，可以由每个特定角度斜入射产生的平面波场所满足的一维方程组反演层状弹性介质的两个参数，在已知密度情况下，可以反演纵波阻抗和横波阻     

水饱和裂纹对地壳岩样中地震波速及各向异性的影响

选择４种地壳岩石样品,经干燥或水饱和处理后在不同围压条件下测量了在其中传播的纵、横波的速度及其各向异性．在大气压条件下低孔隙度(＜１％岩样中,水饱和样品中的纵波速度明显地比干燥样品中的高,但横波速度的差别不大．因为在低孔隙度岩样中纵波速度对孔隙流体的反应比横波速度敏感,可以用泊松比的变化来反映随着围压的增加晶粒间流体对弹性波传播特性的影响．根据实验数据,按Ｏ’Ｃｏｎｎｅｌｌ模型分别计算了干燥和水饱和岩样中的裂纹密度,与通过实测体应变曲线得到的裂纹孔隙度十分吻合．利用横波的速度和偏振特性可以推断岩样中定向排列微裂纹的空间取向情况．研究表明,同时测量在岩样中传播的纵、横波的速度,通过Ｖｐ／Ｖｓ比值可以给出有关颗粒边界流体的证据,也可以估计岩样中的裂纹密度．