非零井源距VSP多分量地震资料逆时偏移

上行波叠前深度偏移是提高非零井源距VSP地震资料成像精度的主要途径之一。从二维各向同性介质中的弹性波动方程出发,在交错网格空间中推导了弹性波场逆时延拓的高精度有限差分算子,给出了弹性波场逆时延拓的最佳匹配层吸收边界条件,从时间一致性准则出发,通过求解程函方程获取VSP地震资料逆时偏移的成像条件,在此基础上实现了非零井源距VSP多分量地震资料联合逆时深度偏移。理论模型与实际资料均得到了高质量的成像结果。     

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

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

基于归一化互相关成像条件的GPR逆时偏移成像

针对基于零时刻成像条件的探地雷达(GPR)逆时偏移精度低、难以对复杂结构进行高精度成像的缺点,将归一化互相关成像条件应用于GPR叠前逆时偏移成像。从二维GPR电磁波方程出发,采用基于完全匹配层(PML)边界条件的时域有限单元法(FETD)模拟电磁波正向和逆向传播,采用归一化互相关成像条件获取叠前逆时偏移的偏移结果,将空间高通滤波用于压制互相关过程中产生的低频噪声,然后编制相应的GPR叠前逆时偏移程序。在此基础上,建立2个复杂的GPR模型,利用基于归一化互相关成像条件的GPR逆时偏移程序进行计算,并与基于零时刻成像条件的GPR逆时偏移剖面进行对比。研究结果表明:与基于零时刻成像条件的GPR逆时偏移剖面相比,基于归一化互相关成像条件的GPR叠前逆时偏移剖面能更清晰地反映异常体空间形态和内部结构信息,其分辨率和成像精度更高。     

起伏地表条件下波动方程法共聚焦点成像技术

为解决起伏地表条件下的偏移成像,引入了共聚焦点成像技术.采用合成聚焦算子和面炮合成的共聚焦点方法和波场逐步累加的直接下延法对起伏地表条件下的模型资料和实际资料进行偏移成像研究.结果表明:该方法克服了起伏地表的影响,直接实现了基于起伏地表的叠前深度成像;采用波场延拓的方法实现了起伏地表条件下的单点聚焦成像和多点聚焦成像;在一定程度上解决了起伏地表的静校正问题,并且计算效率得到提高,在选取较少聚焦点的情况下,就可以达到和共炮偏移一样的成像效果.     

时移成像条件下的波动方程叠前深度偏移

为提高炮域叠前深度偏移的计算效率,引入时移成像条件进行叠前深度偏移.基于频率空间域有限差分波场延拓算子进行延拓,采用时移成像条件来求取大延拓步长间的层位成像值,成像条件虽然增加一次逆傅里叶变换,但对计算时间的影响却不大.脉冲响应和Marmousi模型的测试结果表明,所提出的方法在得到与小步长效果相当的成像结果的同时,又可以成倍地提高计算效率.     

多次波信息偏移成像技术在火山机构精准刻画中的应用

为提高渤海M油田小尺度火成岩的成像质量,对基于单程波理论的多次波偏移方法进行研究。多次波偏移方法采用地震波场一次波与多次波的全波场信息,利用单程波理论对多次波进行互相关成像,能够提高火成岩的照明度及成像分辨率,同时提高了计算效率。由火成岩模型测试及实际数据应用证明,该方法能够改善火山通道的成像质量,提高小尺度火成岩的信噪比和分辨率,使火成岩同相轴刻画得更为清晰。利用该成像结果进行精细地震资料解释能够大幅度提高解释精度。     

基于双平方根算子叠前深度偏移成像技术的结构损伤成像与探测

介绍一种基于双平方根算子的叠前深度偏移方法用于结构损伤探测.首先推导上行波传播时间(走时)的全微分表达式,推导在频率波数域和频率空间域中,用源点和接收点坐标表示的双平方根算子波动方程.然后,用频率空间域的有限差分法及混合法进行双平方根算子的波场外推,推导关于炮点和检波点的差分方程,以及双平方根算子混合法波场外推公式.最后,用Marmousi数据体对结构损伤进行偏移成像.根据双程波动方程分析该方法的数值频散性、奇异性以及波场、位移、振幅、能量、动应力集中系数等力学特征,计算得到波场快照图和剖面图,达到对结构损伤成像的目的.数值算例证明了这种方法的适用性.     

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

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

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

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

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

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