吸收边界条件在探地雷达偏移处理中的应用

在雷达波偏移处理中。通常将边界条件设为零，但是这会在边界处产生由非介质因素引起的反射现象，使得剖面的解释精度降低。因此，为了消除反射边界的影响，提出将吸收边界条件引入到逆时偏移中，并将其运用到实测雷达波剖面上。结果表明。该方法很好地消除了人为边界，精度也得到了明显提高。     

二维矢量地震波场的叠前逆时深度偏移

从矢量波动理论出发，导出了二维弹性波逆时传播的高阶差分格式，实现了弹性波在数值空间中的逆时延拓，采用逆时格式差分求解程函方程，得到网格空间中各点的直达波旅行时，以此作为弹性波逆时偏移的成像条件。实现了二维多波多波多分量资料的叠前逆时深度偏移，数值试验得到了满意结果。     

各向异性介质中弹性波多分量联合逆时深度偏移

当今石油工业界使用的大部分偏移算子都是在假设地球为均匀各向同性介质的前提下经过一定的近似得到的,其应用有一定的局限性。因为经过大量的研究已经证实地壳和上地幔是各向异性介质,不考虑介质各向异性的偏移算子必然带来一些在反射点归位方面的不可估计的错误。从弹性波动方程出发,推导出了适用于各向异性介质中进行弹性波多分量联合偏移处理的偏移算子,并用数据实例进行证明。     

基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

声波方程隐式高阶有限差分数值模拟

推导了声波方程空间二阶导数的隐式求解公式及差分系数的求解方法,讨论了该方法的数值频散特征。利用该方法分别对均匀介质及Marmousi模型进行了数值模拟,将其结果与传统的显式差分格式的模拟结果进行了对比分析。结果表明：该方法较传统的显式求解方法具有更低的数值频散、更高的计算精度。     

任意广角方程逆时偏移的脉冲响应及模型试算

采用高阶差分法求解任意广角声波方程,推导了各向同性介质中任意广角声波方程叠前逆时深度偏移的差分算子,给出了适用于任意广角声波方程逆时延拓的吸收边界条件和其计算方法,在此基础上实现了任意广角声波方程的高精度逆时延拓。分析了不同模型条件下任意广角声波方程叠前逆时深度偏移的脉冲响应,分析结果表明,任意广角方程逆时偏移算子可压制层间反射波,减少偏移噪声,理论模型的成像结果表明该方法能够对复杂构造模型精确成像。     

逆时偏移中的归一化加权互相关成像条件及应用

成像条件是决定声波方程叠前逆时深度偏移效果的重要因素。常规的归一化互相关成像条件能补偿深部地层的反射波能量,但无法压制地震波逆时偏移过程中产生的偏移噪声。针对归一化互相关成像条件的主要缺陷,从二维声波方程出发,在交错网格空间中推导了声波方程逆时延拓的高阶有限差分格式和PML吸收边界条件,进而在对常规互相关成像条件分析的基础上对其进行改进,给出了一种加权互相关成像条件,即在成像过程中利用炮检波场的夹角进行加权校正,以实现只利用炮点下行波场与接收点上行波场进行成像的目标。模型数据和实测资料偏移结果表明,归一化加权互相关成像条件成像效果明显优于常规归一化互相关成像条件。     

含直立裂隙介质的弹性波动方程正演模拟

从具有水平对称轴的横向各向同性（HTI）介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。     

几种叠前深度偏移技术效果的对比

在回顾叠前深度偏移发展历程的基础上,概述了当前几类叠前深度偏移方法:克希霍夫积分法、波动方程法、逆时偏移。在分析叠前深度偏移的波场成像方法、适用性、近地表条件敏感性及存在问题的基础上,对比偏移剖面,比较不同偏移方法的成像效果,指出叠前深度偏移的发展方向,为以后的叠前深度偏移的研究提供参考建议。     

有限差分法探地雷达波动方程偏移成像

在探地雷达（Ground-Penetrating Radar,简称GPR）剖面中，由于绕射波的存在，使得资料的处理解释十分困难，其结果的准确性与真实度也会降低。针对这一问题，作者提出了探地雷达有限差分波动方程偏移法，首先进行了理论模型的实验分析，在此基础上我们对实测GPR剖面资料进行分析处理，取得了较好的成果。     

横向各向同性介质中的弹性波方程逆时偏移及其成像条件

作者在本文中导出了在横向各向同性介质情况下，二维弹性波逆时传播的高阶差分格式，实现了弹性波在数值空间中的逆时延拓，从snell定量出发，采用插值法求取网格空间中各点的直达波旅行时，以此作为弹性波逆时偏移的成像条件，实现多波多分量资料的逆时偏移，数值试验表明，逆时偏移能够使地表接收到的波场准确归位，提高资料处理的精度。     

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

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

转换波波动方程叠前深度偏移

克希霍夫叠前时间偏移被广泛运用于转换波地震成像领域,该方法具有计算效率高、速度模型易于获取的优点。但是,纵横波速度差使纵波和转换波的成像在时间域不匹配,且克希霍夫偏移也无法较好地处理多次到达和保幅等问题。这里提出一种基于波动方程的转换波叠前深度偏移方法,以改善克希霍夫时间偏移的不足。     

VSP数据的逆时偏移

据逆时偏移的基本思想,首次引用傅氏变换和差分法相结合的方法(FDM)数值求解全声波方程,实现对垂直地震剖面(VSP)记录数据的波场外推。同时,利用射线方程实现对震源的向下延拓。从而解决了非零井源距VSP数据的偏移。 所用方法对地下界面倾角较大、速度横向变化没有限制,因而适应于地下构造比较复杂的地质情况。此外,还可用于地面地震记录的迭前深度偏移。     

合成平面波数据的Rytov近似叠前深度偏移

作者把Berkhout(1992)提出的合成平面波技术应用于线性Rytov近似叠前深度偏移，提高了波动方程叠前深度偏移的计算效率。通过在Marmousi模型数据试验，我们认为在达到基本相同成像效果的条件下，利用合成平面波技术的波动方程叠前深度偏移相对于基本共炮集数据的波动方程叠前深度偏移可提高计算效率约10倍。     

油气储集层突变理论识别技术及其应用

本文论述了突变论中的尖点突变模型及其表达式,将地震信号化成尖点突变的标准形式,在时间域或频率域计算突变参数,并将突变参数与Ｂｕｒｇ谱及油气综合预测参数联合应用。实际应用例子表明,突变理论识别技术具有广阔的应用前景。     

叠前深度偏移在复杂构造模型成像中的应用

研究了叠前深度偏移在复杂构造模型成像中的应用。主要对一个复杂的地质模型所正演出的地震数据资料进行试验性处理,并对该试验性处理结果进行归纳分析,得出速度~深度模型的准确性决定于对偏移结果的质量以及根据地下构造选取偏移孔径的原则等定性认识和结论。通过这些认识和结论,对实际地震资料的处理有一定的指导意义。