基于扩展伪多道匹配的保幅型多次波压制方法

波动方程预测减去法压制多次波主要分为两步:多次波模型预测和多次波匹配相减。两步骤同时决定着最终多次波压制效果的好坏。本文在现有SRME方法基础上进行了改进。为了实现保幅条件下的多次波压制,除在匹配过程中采用扩展伪多道匹配方法之外,提出了三点改进:①按多次波阶次进行预测和匹配;②匹配前有效能量保持;③匹配过程中迭代应用短滤波器。通过Sigsbee2B模型及实际资料试算,表明本文方法能够在较好保持有效能量的前提下得到更为理想的多次波压制剖面。     

多轴卷积完全匹配层吸收边界条件

在通常情况下,完全匹配层（PML）吸收边界条件较其他吸收边界条件具有更优越的吸收性能,已广泛用于地震波正演模拟。但传统的PML吸收边界条件也存在一定缺陷,如不分裂卷积完全匹配层（C-PML）吸收边界在某些介质中不稳定,分裂的多轴完全匹配层（M-PML）吸收边界的吸收效果不好。为此,本文提出了多轴卷积完全匹配层（MC-PML）吸收边界条件,既可提高C-PML吸收边界的稳定性,又可改善M-PML吸收边界的吸收效率,且采用不分裂算法实现,可有效降低计算成本。正演模拟结果表明,MC-PML吸收边界条件具有更高的计算稳定性和较好的吸收效果。     

TTI介质一阶qP波稳定方程波场数值模拟及逆时偏移

相比于各向异性介质弹性波波动方程,利用声学近似的各向异性介质qP波方程进行波场数值模拟及逆时偏移更具优势。常规的声学近似方法往往会造成非均匀TTI介质中倾角剧变区域出现数值不稳定。为此,基于精确的TTI介质qP-qSV波耦合频散关系,首先引入一个各向异性控制参数σ,推导了新的TTI介质二阶qP波稳定方程,并通过引入波场的伪速度分量,将其转换为等价的一阶应力-速度形式波动方程。然后,利用优化的最小二乘交错网格高阶有限差分(LS-SGFD)方法数值求解TTI介质一阶qP波稳定方程,构建波场延拓算子,实现了精确的各向异性介质波场模拟及逆时偏移成像。模型试算结果表明,TTI介质一阶qP波方程能够稳定地模拟qP波的波场传播特征,利用优化的LS-SGFD方法能够有效地提高波场模拟的精度,进一步可以改善偏移成像质量。     

构造倾角约束的有效能量叠加成像方法及应用

基于反射波、绕射波以及野外资料中的噪声在倾角域共成像点道集(Dip-angle Domain Common Image Gathers,DDCIGs)中的响应特征,本文实现了一种基于构造倾角约束的有效能量叠加成像方法。通过高斯束叠前深度偏移提取倾角域共成像点道集,在DDCIGs中识别有效的反射和绕射能量及低信噪比数据引起的噪声,并结合偏移剖面和探区地质信息划分出叠加区块,进行基于构造倾角约束的成像处理。Marmousi模型测试以及实际资料偏移结果表明:本文方法可以有效压制偏移噪声,提高剖面信噪比及增强剖面同相轴的连续性,具有很好的人机交互功能,便于在实际资料处理及参数选取中进行质量控制。     

有效邻域波场近似框架下三维起伏地表高斯束地震正演模拟

在传统高斯束正演模拟的基础上,针对三维起伏地表模型使用有效邻域波场近似导出了三维起伏地表高斯束正演公式,并给出了笛卡尔坐标系下相应的实现算法.对典型的三维起伏地表模型进行试算,结果表明:本文方法可以较好地刻画由起伏地表引起的波场畸变和由复杂构造引起的同相轴干涉现象,充分体现了高斯束在处理起伏地表条件的优越性;试算得到的射线路径及单炮记录证实了本方法的正确性以及对三维复杂地表模型较好的适应性.     

高斯波束角道集的旅行时层析速度分析

本文综合高斯波束叠前深度偏移准确、灵活、高效及其角道集能够直接反映成像点偏移深度与局部入射角之间的关系,避免多波至带来的假象干扰等优势,提出了基于高斯波束角道集的旅行时层析速度分析方法,推导了灵敏度矩阵和角道集剩余曲率之间的定量关系式,并利用自动拟合策略求取剩余曲率。模型和实际资料试算验证了该方法具有较高的反演精度和计算效率。     

基于成像域的多次波压制

常规多次波压制方法大都是在数据域基于滤波理论或者波动理论实现的,这些方法都不能有效地兼顾多次波压制精度和计算效率。为此,本文将波动理论与滤波理论相结合,基于成像角度域共成像点道集(ADCIGs)实现多次波压制。首先对成像域多次波映射进行了讨论,通过推导可知,成像域ADCIGs中多次波和一次波剩余时差不同,这一结论即为在成像域应用滤波方法压制多次波的理论基础;考虑到基于DSR方程的叠前深度偏移可以获得更多的角度信息,然后基于DSR方程提取ADCIGs,再引入适用于ADCIGs的正切平方估计高分辨率拉东变换分离一次波和多次波,实现多次波压制。由于本文方法将波动理论与滤波理论相结合,因此能同时兼顾多次波压制精度和计算效率,模型试算和实际资料处理验证了方法的有效性。     

一阶弹性波方程的变网格高阶有限差分数值模拟

 使用可变网格的有限差分法进行地震模拟有许多独特的优点，主要表现为对地质模型的离散化更为合理，在低速带和复杂构造区域，可将局部网格划分得相对精细些，不仅提高了模拟精度，消除了因采样不足导致的频散现象，而且可以减少计算机内存需求，保持模型计算的灵活性。本文提出一种新的基于高阶交错网格技术的弹性波数值模拟方法，通过改变网格的空间步长实现了局部网格加密技术，弥补了常规网格的缺陷和不足。试算结果表明，本文提出的算法稳定性较好，且能够提高模拟精度，减少计算时间，提高计算效率。     

双相黏弹VTI介质一阶速度—应力方程正演模拟及双程波照明研究

地层岩石具有各向异性、黏弹性及双相性等特性,因此建立更加精确的地下介质模型,研究地震波传播规律,对认识复杂介质中地震波的传播特征和实际资料的解释有着重要的实际意义。为此,基于黏弹性广义标准线性体（GSLS）模型,首先推导了双相黏弹VTI介质的一阶速度—应力方程并进行正演模拟,与单相弹性各向同性介质相比,双相黏弹VTI介质中地震波能量明显衰减,波前面变成椭圆,且出现慢纵波。对二维双相黏弹VTI介质进行地震双程波照明分析,结果表明介质的黏弹性使地震波能量衰减,且介质各向异性导致地震波的传播速度随方向变化,双相介质慢纵波的出现使地震波能量分布复杂化。最后,根据照明结果对观测系统进行了优化。     

基于波动方程的聚焦点控制照明叠前深度偏移

 基于波动方程的聚焦点控制照明叠前深度偏移技术借助于差分计算，把速度、密度等介质参数的影响体现在差分计算的矩阵方程中，能够自动适应速度场的任意变化，快速傅里叶变换的使用也加速了波场延拓的计算速度。因此此法兼具有限差分偏移方法和傅里叶偏移方法的优点，既可适应速度场的剧烈变化，又可保证对陡倾地层的成像效果，是目前针对复杂构造最有效的成像方法之一。对于单个聚焦点及其周围的成像步骤为：①采用矩形网格情况下绕射走时的有限差分计算方法生成合成算子；②应用合成算子来合成面炮震源和面波记录；③对合成的面炮震源和面波记录做傅里叶有限差分法波动方程叠前深度偏移，得到该聚焦点及其附近区域的成像结果。按照上述成像步骤，将震源波场和炮集记录依据相应的外推公式进行延拓，最终应用成像条件求取成像值。在地质目标处选取多个聚焦点，可以得到面向目标的控制照明偏移成像，在多个层位上选取多个聚焦点进行控制照明叠前深度偏移，可以得到整块区域的成像。通过对Marmousi模型的试算，取得了较好的效果。