任意广角方程逆时偏移及其在煤田资料中的应用

利用任意广角方程逆时延拓的高阶有限差分格式实现了纵波资料的叠前逆时深度偏移,给出了稳定性条件及吸收边界条件;采用归一化互相关成像条件,在成像条件的计算中摒弃双程波方程,采用单程波方程对检波点逆时延拓波场和炮点正向延拓波场进行计算,得出了最佳归一化互相关成像条件.该算法在煤田资料偏移处理中,对较深处的构造清晰成像,取得了较好的应用效果.     

弹性波叠前逆时深度偏移技术及其应用

研究了二维各向同性介质中两分量记录的基于弹性波标量波动方程的叠前逆时深度偏移方法.从弹性波波动方程出发,推导各向同性介质中弹性波逆时深度偏移技术中相关算法,并探讨其中有关问题的解决方案,实现二维多分量资料的叠前逆时深度偏移,最后进行模型试算和实际资料试处理.结果表明,叠前逆时深度偏移方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的多分量地震资料处理.     

隧道反射波超前探测有限差分正演模拟与偏移处理

针对地下工程领域隧道不良地质灾害超前探测问题,构建断层、倾斜软弱夹层等地质灾害体的隧道介质模型,研究隧道反射波超前探测的波场传播机理与偏移成像问题。利用一阶速度-应力弹性波方程和高阶交错网格差分计算方法,导出隧道工程反射波探测数值模拟的差分计算格式、数值稳定性条件以及边界条件,对上述模型进行正演模拟,得到相应模型不同时刻的波场快照与地震合成记录;通过对这些模拟结果进行分析,验证高阶交错网格差分法在隧道波场正演中的有效性。在导出弹性波方程逆时延拓交错网格差分计算格式基础上,通过求解程函数方程,获取激发时间成像条件,进而得到隧道全波场逆时偏移算法及偏移处理程序,并应用于隧道波场偏移成像中。研究结果表明:采用高阶交错网格有限差分逆时偏移算法能使反射波归位,绕射波收敛,从而大大提高了隧道工程反射波探测的分辨率,为隧道超前探测资料的解释提供理论依据。     

各向异性介质弹性波场正演及偏移

大量数据表明地下岩层的地震各向异性是普遍存在的.然而在常规处理中,通常视地下介质为各向同性介质,这必然导致一定的误差.作者从弹性波动方程出发,用虚谱法正演模拟了横向各向同性介质中弹性波场,用弹性波多分量联合逆时深度偏移法,并在计算中很好地解决了边界反射问题,实现了多波多分量资料的偏移,获得了一些与各向异性介质有关的波场信息和认识.这些信息和认识有助于复杂地层成像和岩性参数反演.     

三维VSP波动方程保幅叠前深度偏移方法

VSP偏移成像技术在VSP资料处理中具有重要意义。为了对地下复杂构造进行精确成像,常要求对VSP数据进行三维叠前偏移处理。传统的单程波方程延拓计算时只能保证相位信息的正确性,没有对振幅做任何处理,仍存在保幅性较差的问题,不利于将其结果运用到岩性分析中去。本文以保幅型单程波偏移算子为基础,研究了三维VSP单程波动方程保幅偏移成像方法。并在实际资料中,采用了保幅型三维地震道插值方法,对非规则网格炮点采集的三维VSP数据进行插值,保证了波动方程保幅偏移的效果。     

基于波场外推的3D VSP波动方程深度偏移成像方法研究——以胜利垦71-检41井为例

针对3DVSP成像方法中的保幅性差、对复杂构造的归位不准、分辨率低等问题,提出基于波场外推的波动方程深度偏移成像方法。该方法以三维声波方程为出发点,对有限差分法进行误差补偿,多方面确保成像的稳定性,并充分考虑到波场在地下传播过程中的绕射和折射效应,避免了Kirchhoff方法中的高频近似假设和在复杂构造区的焦散及多路径现象。实现过程是将三维VSP数据按井中检波器的位置抽成共检波点道集,对共波点道集进行三维叠前深度偏移,对偏移后的数据进行同深度叠加。该成像方法在胜利垦71井区3DVSP实际资料处理中取得了比较好的成像效果。实际处理结果证明：该方法正确、可行、实用。     

波动方程保幅叠前深度偏移与AVO响应

从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的单程波动方程;按照叠前深度偏移方法中的波场双重向下延拓理论实现共炮域公式到CMP(共中心点)域公式的转换,推导出由双平方根方程定义的保幅单程波动方程,从而将振幅误差补偿作为偏移的一部分在偏移过程中实现,使得地震成像在给出正确位置的同时也给出真实振幅.在此基础上获得的CRP(共反射点)道集在保证真实振幅的前提下,有效解决了绕射波干涉、反射点弥散、与倾角有关的NMO等问题,使AVO响应更加清晰,提高了AVO资料的分析精度.     

基于惠更斯原理的波动方程共炮点道集地震正演

地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义.作者利用下行波将震源场延拓到地质介质任意深度,实现单程波动方程模拟叠前正演记录.用单程波动方程法对一个盐丘地质模型生成129个共炮点地震记录,并对所有共炮记录进行叠前深度偏移.在共炮记录中,来自盐丘地质模型的反射波和绕射波十分清晰.叠前深度偏移剖面和盐丘地质模型完全吻合.该方法生成记录质量好且计算效率高,并且叠前正演记录没有多次波、直达波和其他干扰波.     

全波动方程有限差分法倾角时差校正

首先从倾角时差极正（DMO）微笑响应方程出发,导出了共偏移距地震剖面DMO波动方程,在浮动坐标系中,该方程化为一个只与炮检距有关而与速度无关的二阶偏微分方程,用时间步长上的分裂方法对此方程进行分裂,获得了以地面观测波场为这值,最大记录时间波场为初值的波动方程波场外推法DMO的定解问题.由于时间步长上的分裂方法没有任何近似,所以理论上讲能适应任何倾角反射波的成像,同时由于采用浮动坐标变换,波场只需延拓到t＝τ,处即可以成像,这样和其他波场外推法DMO相比计算量也减少了一半.最后给出了此方法子波响应的计算结果并就计算精度做了定性分析     

虚谱法无反射速度-应力方程地震数值模拟

基于单程波波场廷拓算子的零偏移距地震数值模拟方法难以适应复杂地质条件下地震波速度的空间剧烈变化,且其最大成像角度有不同程度的限制.本文依据地震反射的形成机理,重置模型空问的密度分布,使得模型空间内纵波阻抗为常量以压制层间多次波;采用虚谱法可得到高精度的地震波场空间导数,因此可有效抑制空间数值频散,同时提高了数值模拟的垂向分辨率.数值模拟结果表明,得到的零偏移距地震记录绕射波场完整,偏移归位效果优于频率-波数域或频率-空间域的单程波模拟结果.无反射一阶应力-速度方程零偏移距数值模拟可有效压制多次反射波,在较大网格间距时,虚谱法仍可有效抑制数值频散.