基于互相关成像条件的保幅高斯束逆时偏移

常规的逆时偏移算法都是基于双程波动方程的数值解法,实际上利用Kirchhoff积分进行波场延拓,也可以实现逆时偏移。以Kirchhoff积分为基础,采用高斯束叠加积分表示的格林函数来表示正反向延拓波场,将正反向延拓波场进行保幅的互相关得到地下成像。保幅高斯束逆时偏移方法结合了射线类偏移方法的高计算效率和逆时波动方程偏移的高精度,能很好地处理焦散点、陡倾角等问题,并且具有面向目标成像的能力,并消除了偏移过程中的几何扩散效应和入射角变化引起的振幅误差。     

基于修正伪Helmholtz分解算子的横向各向同性介质纵横波波场分离

纵横波波场分离是弹性波逆时偏移成像中的关键步骤,可以有效消除纵横波串扰和成像畸变,提高成像精度。在各向异性介质中,常用的波场分离法有非平稳滤波器或者低秩近似的方法,但由于使用多次傅里叶变换,导致波场分离计算成本较高。借鉴基于Helmhotz算子构建各向同性介质纵横波解耦方程的思路,提出消除振幅畸变的修正伪Helmholtz分解算子,通过待定系数方法求解横向各向同性介质中P波和S波2种模式的表达式,并将其变换为一阶解耦的伪弹性波方程,实现时空域纵横波波场解耦。通过简单模型的波场分离测试,得到分离后的P波与S波波场,验证了该方法的波场分离的有效性。此外,将解耦得到的矢量P波与矢量S波应用到弹性逆时偏移之中,利用矢量点积互相关成像条件得到清晰的弹性波逆时偏移成像结果,说明本方法在复杂介质中有着较好的适用性,同时可以有效地应用到VTI介质弹性逆时偏移之中。     

弹性波叠前逆时偏移

利用弹性波一阶速度—应力方程,采用交错网格高阶有限差分法进行震源波场的模拟和接收波场的逆时外推,并结合了能流密度矢量的互相关成像条件,实现了弹性波叠前逆时偏移,且能很好的消除叠前逆时偏移中的假象干扰.笔者通过一个二层模型和一个金属矿模型进行了叠前逆时偏移成像,结果证明所述方法可以有效地消除假象干扰并能对复杂陡倾角界面进行清晰地成像.     

基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

三维混合延拓一步法波动方程正演模拟

提出三维相位移加有限差分混合延拓法,以此为基础,实现速度纵横向变化,构造任意复杂的三维介质一步法正演模拟,为三维复杂波场的地震,地质解释提供了精确,实用的方法。对于二维介质,相位移延拓方法对陡倾斜地层成像精度高,稳定性好,计算快,但实现的速度的横向变化困难,而频率、空间域的有限差分法算法简单,稳定性好,能适应速度的纵横向任意变化,但偏移陡倾地层存在很多问题,用快速的４５°有限差分对相位移延拓作补充     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

地震波逆时偏移中的层位校正与去噪方法

波动方程逆时偏移是目前比较精确的基于波动理论的深度域偏移成像方法,但逆时偏移成像条件的运用会产生低频噪音,影响成像质量及层位的准确性。依据低频噪声的产生机理,运用交错网格有限差分法求取各向同性介质中的地震波场,通过对比炮点与检波点波场在不同频率条件下的互相关成像特征,分析其子波频率与周期差异(或时移)条件下逆时偏移成像对层位解释精度与低频噪音产生的影响,提出了改进的时移逆时偏移成像层位校正方法和低频噪声压制策略。数据算例表明,改进后的时移逆时偏移成像有效的压制了低频噪声,成像层位准确。     

基于互相关成像条件的VTI介质保幅高斯束逆时偏移

保幅高斯束逆时偏移是一种兼具高斯束偏移灵活、高效和逆时偏移高精度特性的深度域成像算法,不仅能够消除偏移过程中由于入射角变化和几何扩散引起的振幅误差,而且具有面向目标成像的能力。其以Kirchhoff积分为基础,采用高斯束的加权积分来构建格林函数,进而实现正反向波场的延拓,最后通过相应的成像条件进行成像。这里从基于相速度的各向异性射线追踪算法出发,采用互相关成像条件,将保幅高斯束逆时偏移推广到VTI介质。随后通过模型试算对方法的正确性和有效性进行了验证。     

波场延拓反Q 滤波的正则化方法

基于波场向下延拓的反Q 滤波是补偿大地衰减效应的一种常用方法,但其振幅补偿的不稳定性限制了该方法在中、深层地震资料处理中的应用。针对于此,提出波场延拓反Q 滤波的正则化方法,该方法不仅较好地解决了地震波能量补偿的不稳定问题,同时对经过品质因子较低地层的深层地震波有着较好的能量补偿效果。模型试验和实例应用均表明: 与常规稳定的反Q 滤波方法相比,基于正则化法的反Q 滤波不仅具有较强的抗干扰能力和较好的振幅补偿能力,而且有效地提高了中、深层地震资料的分辨能力。     

基于P＋S波震源的弹性波叠前逆时偏移方法

基于弹性波动方程,推导了二维非均匀介质条件下纵横波解耦且无转换波的一阶速度一应力弹性波方程,使得弹性波在延拓过程中实现矢量偏移保证偏移质量,且不产生其他类型的转换干扰波。从激发时间成像条件以及反射系数理论出发,并与纵横波解耦且无转换波的一阶速度一应力弹性波方程结合,给出了适用于P＋s波震源模拟数据的成像条件以及偏移方法。文中给出了两个典型模型的理论记录的偏移结果,试验表明,笔者提出的方法能够对P＋s波震源的产生的波场准确归位,提高资料处理精度。     

地震波逆时偏移中两种成像条件应用效果对比

成像条件是实现地震波场准确成像的关键环节,而优化逆时成像条件在改善逆时成像应用效果中的作用尤为明显。为此,从零井源距和非零井源距脉冲响应出发,系统对比了常规相关法逆时成像条件和基于行波分离法逆时成像条件的逆时偏移脉冲响应波场差异,分析了逆时偏移低波数噪声的成因机理及波场规律,同时以复杂断陷理论模型和SZ工区实际三维地震资料为例,系统对比了两种逆时成像条件的成像应用效果。研究结果表明,采用基于行波分离法逆时成像条件的成像结果能够有效衰减逆时偏移低波数背景噪声,恢复掩盖在这种背景噪声之下的有效地震成像信息,地层细节刻画更加清晰,这可为当前复杂波场和复杂构造高精度地震成像提供方法指导。     

基于谱元法的各向异性黏弹性介质地震波场模拟和分析

基于谱元法计算理论,通过坐标旋转算法实现任意各向异性黏弹性介质中的地震波场计算,研究复合介质中地震波场特征。根据计算结果,qSV波表现为剪切滞弹性形变,qP波表现为膨胀滞弹性形变;品质因子不仅影响地震波振幅,还影响其传播频率;黏弹性介质吸收高频快,在一定频率范围内,黏弹性介质没有频散现象。     

地震波的有限差分法正演模拟

基于波动方程的地震波正演模拟能够准确揭示波在介质中传播的振幅、频率与相位变化,真实地反映出波的动力学特征。有限差分法因其精度高、计算效率高而成为波场数值模拟的一个重要方法。根据地震波的近似方程———声波方程,采用二阶有限差分网格数值算法,利用透明边界条件,设计典型地质模型,进行地面波场的地震记录模拟。通过稳定性条件分析、震源子波优选以及频散抑制,成功合成了倾斜地层模型与正断层模型的地面共炮点波场记录,直观地显示出直达波、反射波、绕射波等丰富的波场信息,获取各模型的波场快照。结果表明,有限差分法结合透明边界条件对于获取介质丰富的波场信息具有很大的优势。     

双复杂介质条件下的反Q滤波偏移延拓算子研究

地下介质的非弹性吸收衰减效应造成地震资料中深层信号频带窄、能量弱、成像困难,因此需要对地层的非弹性吸收衰减进行补偿。在已有的补偿方法基础上,基于单程波方程,采用与品质因子和频率相关的复速度代替介质真实速度,将粘滞传播算子在频率—空间域采用连分式展开,得到相移项和补偿项,针对连分式的误差,在频率—波数域进行补偿校正,并用逐步累加法解决起伏地表问题。该方法提高了成像精度,有效补偿了由于地层非弹性吸收造成的振幅和频率损失。模型试算和实际资料处理表明,该方法具有一定的适应性和实用性。     

频率域单程方程三维叠前正演中的波场替换

利用定位原理进行三维叠前正演,仅靠单程声波方程的延拓即可形成复杂地质条件下的正演记录,运算效率高无多次反射出现。但速度模型中对于过大或过小的速度常使波场在延拓时有频率和波数成分的损失,这种损失引起波场畸变并严重影响正演记录的质量。本文采用波场替换技术,在每次延拓时用零相位子波形成一个同时保留了原波场的动力学特征与空间位置信息的新波场,然后用新波场取代原畸变波场并参与正演记录合成。在频率域单程方程三维叠前正演中特别适合做波场替换,因为震源脉冲与接收响应(均可视为点脉冲)的延拓波场在形态上相似,而且在延拓波场中每道记录之振幅包络只有一个极值点。理论试算结果表明,采用波场替换技术后最终正演记录质量比替换前有了很大改善。      

叠前弹性波逆时深度偏移及波场分离技术探讨

这里探讨了Sun和McMechan提出的叠前弹性波标量逆时深度偏移方法,即在地表附近,对地表接收到的弹性波波场分量进行波场逆时延拓,然后分别对在波场延拓过程中通过波场分离得到的纯纵、横波反射波场进行偏移。在模型中,每个网格点的纵、横波成像时间,为震源到该网格点的纵波初至时间。利用纵波速度模型和声波传播方程的有限差分解,对反射纵波进行逆时延拓和成像,利用横波速度模型和声波传播方程的有限差分解,对转换横波进行逆时延拓和成像,并在偏移前进行极性校正。然后,在Sun和McMechan偏移成像方法的基础上,提出了直接偏移成像方法,即对地表接收到的弹性波波场分量在整个模型范围内进行波场逆时延拓。在波场延拓过程中,对符合成像条件的网格点进行波场分离、偏移成像和转换波极性校正。最后,通过数值模拟对二种偏移成像方法的结果进行了比较。     

桩井爆破振动在邻近边坡产生的特殊效应

为研究桩井爆破地震波在邻近自然边坡传播规律,对桩井开挖爆破振动速度进行了现场测试。在邻近桩井的台阶坡脚及边缘各布置一个测点。测试结果表明,位于台阶边缘测点的径向、垂向振动速度及主频率较位于坡脚处的测点存在增大趋势,出现了坡面效应;切向振动速度及主频率存在减小的趋势,说明爆破地震波的速度及频率变化具有一致性,且爆破地震波的放大效应具有方向性,以垂向为主。坡面效应合理地解释了爆破地震波在台阶边缘的振速放大效应,同时通过“鞭梢效应”对爆破地震波频率的分析,验证了坡面效应的合理性。     

基于波场延拓的反Q滤波方法比较

波场延拓法是最常用的反Q滤波方法,但其振幅因子的不稳定性致使深层的高频信息放大,从而影响反Q滤波的效果。本文通过推导证明了完全稳定振幅补偿只和Glim(增益限制)有关难以提高深层地震波分辨率,系统分析了波场延拓反Q滤波算法中不同振幅补偿算子的滤波特性,增益控制法对某一频率之后既不放大也不压制,正则化滤波和迭代滤波法不仅具有较强的抗噪声干扰能力,而且对深层地震波具有较好的振幅补偿能力。通过对理论模型和井间地震实际资料应用不同振幅补偿算子,发现正则化振幅补偿因子和迭代法振幅补偿因子不仅对高频干扰有较强的压制效果,而且对于地震记录的中深层部分和低Q能量衰减严重的地震记录,都有较好的补偿效果。     

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

地震勘探基于波场分离的逆时偏移去噪方法

基于双程波动方程的逆时偏移（RTM）具有原理简单、无倾角限制、可以适用于速度任意变化的模型,并具有对各类波正确成像等优点,是当今最好的偏移方法之一.本文阐述了逆时偏移存在的缺陷并分析了低频噪声的产生机制,重点探讨了使用波场分离的逆时偏移去噪方法,并通过模型的试算证明了该方法的有效性.