地震波动方程方向行波波场分离正演数值模拟与逆时成像

为了进一步提高对地震波传播规律的认识,将波印廷矢量引入到二维地震波动方程方向行波波场分离正演数值模拟中。 根据地震波波印廷矢量的波场数值特征,实现了对全地震波场中左行波、右行波、上行波和下行波波场的自动识别与分离。以均匀介质模型、倾斜界面模型以及 Marmousi 模型为例,开展了相应的数值模拟实验与逆时成像处理。 计算结果表明,该方法准确有效,能够实现任意时刻波场快照中方向行波的波场分离,并合成分别由左行波、右行波、上行波和下行波形成的波场快照与数值模拟记录。该方法简单易行,计算量较小,对实际地震资料中方向行波波场的识别、分离、成像及验证具有一定的应用价值。     

各向异性弹性介质方向行波波场分离正演数值模拟

为了进一步提高对弹性波波场传播规律的认识,将波印廷矢量应用于多波多分量各向异性介质弹 性波波动方程方向行波波场分离正演数值模拟中。 根据弹性波波印廷矢量的波场数值特征,在多分量弹 性波正演数值模拟过程中,实现了上行波、下行波、左行波和右行波的方向行波波场分离。 以均匀各向异 性弹性介质模型、倾斜界面模型和 Marmousi 模型为例,开展了相应的方向行波波场分离数值模拟实验。 结果表明,这种方法计算量小,算法简单,能够准确实现波场快照和数值模拟记录的方向行波波场分离。 因此,在多波多分量弹性波资料的地震波场模拟分析和成像方面具有一定的应用价值。     

地震波逆时偏移中不同域共成像点道集偏移噪声分析

地震波逆时成像方法通常输出2种共成像点道集：一种是共成像点偏移距道集,另一种是共成像点角度道集。开展了针对地震波逆时偏移的2种共成像点道集波场特征及其差异的数值实验。以倾斜界面模型和复杂的Marmousi模型为例,合成了2种共成像点道集,其中角度道集是在逆时延拓过程中采用波印廷矢量制作的。计算结果表明,在速度模型准确的情况下,2种共成像点道集均可拉平,其中偏移噪声在偏移距道集上的分布规律较差,不利于后续的精细处理,而在角度道集上的偏移噪声主要集中在90°附近的高角度区域．对小角度范围内的地震道进行叠加时,逆时成像剖面的信噪比和地层的刻画能力均可得到有效提高,同时低频噪声能量与界面的反射系数存在正相关性。     

弹性波联合叠前逆时偏移数值试验

与常规纵波地震资料相比,多波多分量地震资料含有更为丰富的弹性波场信息,因此,研究多波多分量弹性波场的逆时成像方法可以更好地挖掘其潜在的应用价值。从非均匀各向异性介质弹性波波动方程出发,尝试从弹性波正演过程来计算纵波/准纵波初至走时,并以此走时作为弹性波联合叠前逆时成像值的筛选条件;同时给出了具体的逆时成像原理和计算步骤,并引入基于散度和旋度的纯纵、横波的波场分离算子,最终实现了弹性波波场中混合波场、纯纵波和纯横波(转换波)波场的联合逆时成像。均匀弹性介质模型算例验证了所提出的走时计算方法的准确性和有效性;各向异性弹性介质模型的数值模拟试验结果表明,所给出的弹性波联合叠前逆时成像方法能够考虑介质岩性和构造特征的复杂变化,获得符合理论和实际地质情况的成像结果。     

TI介质中qP波方程逆时偏移技术的研究现状与展望

理论上说,TI介质逆时偏移技术更有利于解决各向异性介质的地震波成像问题,工业实践也证明了此项技术在解决复杂地层地震波成像问题方面的巨大潜力。总结了TI介质中qP波逆时偏移技术的研究现状,重点分析了qP波方程的建立、qP波方程的延拓算法和qP波成像方法这3个逆时偏移关键技术的实现思路和存在的主要问题,并在此基础上探讨了今后该领域的研究重点。分析认为,现有技术主要存在以下问题:①已有的TI介质qP波方程具有较高的运动学精度,但对方程的动力学精度缺乏深入研究;②qP波方程炮点波场重构中的随机边界技术存在不稳定隐患,业界在提升稳定性方面缺乏有效的解决方案;③二阶或高阶qP波方程难以准确求取qP波的坡印廷矢量,且现有的坡印廷矢量技术无法解决多组波场问题;④子波拉伸校正缺乏针对性技术等。因此,今后TI介质中qP波方程逆时偏移领域的研究重点应该包括qP波方程的动力学精度分析多参数随机边界设置技术和qP波传播方向的准确求取等方面。     

波动方程地震定向照明分析

 本文从地震波能量理论出发，提出了基于玻印亭矢量的地震波场方向性分解方法和地震波定向照明度计算方法。根据声波方程高精度正演模拟，计算出声波波场的能量密度矢量和能量传播方向，并通过角域滤波器实现对波场的方向性分解；进而基于方向性分解后的波场，可以计算出震源定向照明度和源—检定向照明度。数值模型试验结果表明：该法计算精度高且速度快，具有非常良好的实用性；定向照明可以用于定量分析地震采集参数对成像质量的影响；目标层CRP点照明角分布可用于地震分辨率分析。     

地震波叠前逆时偏移脉冲响应研究与应用

为了更好地应用地震波叠前逆时偏移成像技术,开展了逆时偏移脉冲响应研究.从有偏移距层状介质模型的脉冲响应出发,首先探讨了常规相关法逆时成像条件低频逆时噪声的成因机理,分析了上、下行波分离逆时成像条件在低频噪声压制中的优势,对比了两种逆时成像条件在偏移速度平滑条件下的脉冲响应特征;同时分析了不同偏移网格参数、速度分布等情况下地震波数值频散对逆时偏移脉冲响应精度的影响.计算结果表明,地震波波场延拓过程中相对波阻抗差界面的上行反射波场是逆时背景噪声的主要成因;速度平滑能够有效降低低频背景噪声,但以牺牲成像位置准确性为代价;地震波数值频散是影响逆时成像精度的一个重要因素,逆时偏移参数的选择必须严格依据数值频散关系.与此同时,还通过模型验证了逆时偏移中震源和检波点位置具有可互换性,为共检波点道集的逆时偏移处理提供了方法指导.通过逆时偏移参数的优选,在实际陆上地震资料叠前逆时偏移中取得了较好的成像效果.     

多分量波场的矢量法叠前深度偏移技术

本文研究了二维各向同性介质中两分量记录的矢量法叠前深度偏移方法———叠前逆时深度偏移方法。首先从弹性波波动方程出发,在交错网格空间中推导了各向同性介质中弹性波逆时延拓的高阶有限差分格式;然后从程函方程出发,采用逆时差分格式求解地下各点的地震波走时。为确保算法满足地震波传播的因果性条件及其对复杂模型的适应性,用扩展波阵面算法追踪波前并搜索全局极小,以上述方法的计算结果作为弹性波逆时偏移的成像条件,实现二维多分量资料的叠前逆时深度偏移。模型试算和实际单炮记录试验表明,叠前逆时深度偏移方法考虑了地震波的矢量特征,是一种有效的矢量波场处理技术。     

各向异性弹性波高阶有限差分法叠前逆时深度偏移

采用高阶交错网格有限差分法算子构建了各向异性介质多分量弹性波动方程的高精度正演模拟和叠前逆时深度偏移的数值离散方程,在正演过程中采用最大绝对振幅能量法实现各向异性介质多分量初至旅行时的计算,并以此作为多分量双程弹性波叠前逆时成像条件,同时给出计算所需的稳定性条件、多分量叠前逆时成像的基本原理、吸收边界条件等,采用凹陷模型合成了共炮点道集,并进行多分量双程各向异性弹性波叠前逆时成像研究,并与对正演炮集进行多分量各向同性叠前逆时成像结果进行对比。计算结果表明,考虑了各向异性弹性波场的矢量特性,能够更为准确地实现叠前多分量弹性波场的成像问题,使地质层位中的断层、断点等复杂目标成像更加清晰准确,并且偏移成像精度较高,因此开展各向异性叠前逆时成像可为当前高精度岩性地震勘探提供方法指导。     

双程声波方程叠前逆时深度偏移的成像条件

 成像条件是决定波动方程叠前逆时深度偏移效果的重要因素之一。本文从构造成像角度研究了不同的成像条件计算方法在双程声波方程叠前逆时深度偏移中的应用效果。模型试验表明,不论采用何种成像条件,基于双程波方程的叠前逆时深度偏移均会产生较强的层间反射,该反射会恶化浅层成像结果。应用波场归一化互相关成像条件时,层间反射的影响相对较弱;基于波场互相关的逆时偏移成像条件对深部构造的成像能力总体上优于求解程函方程得到的成像条件;而归一化互相关成像条件能够在对地震波进行偏移成像的同时补偿深层能量,在相同条件下,它对深部地层的成像能力更强。