TTI介质的交错网格伪P波正演方法

研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果.     

半径-斜率域加权反假频地震数据重建

为减小地震数据缺失给地震后续处理工作带来的影响,需要对地震数据进行插值重建.针对反假频插值重建这个难点问题,进行了相关研究,并由此提出了一种改进的R-P（半径-斜率）域加权反假频地震数据插值重建方法.该方法将F-K（频率-波数）谱变换到R-P域,在R-P域设计一个权函数并将其作用于每次的迭代插值过程.通过模型数据和实际数据的测试,证明了该方法具有较好的反假频插值重建能力.     

零偏VSP反演Q值CFS方法及影响因素研究

当地震波在地层中传播时,能量的衰减是传播介质非弹性性质的反映.这种介质所固有的衰减特性通常用品质因子(Quality Factor, Q)来描述.本文基于地震波在传播过程中主频降低的衰减特性,利用质心频率偏移法(Centroid Frequency Shift,简称 CFS)对零偏VSP正演模型进行Q值计算,分析薄层、上下行波场、频带宽度、反演波速、震源位置等因素的影响.模型结果显示:CFS法比频谱比法、振幅衰减法能更加准确识别薄层界面,可准确反演出厚度10m的地层;上行波场的加入降低反演结果准确度,尤其对高Q值层;当有效波频带增宽时,高频部分衰减明显,Q值反演结果接近理论值;反演波速误差对反演准确性影响很大.最后,利用CFS对实际井旁地面地震资料进行Q值反演,反Q滤波后其主频由25Hz提高到35Hz,分辨率得到改善,证明了CFS方法的实用性.     

数据自相关多次波偏移成像

在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义.     

叠前逆时偏移影响因素分析

反射地震勘探中的偏移成像技术是获取地下介质构造形态最有效的手段之一.在叠前深度域偏移方法中,目前工业界采用的方法包括基于射线理论的波动方程积分解法和基于波动理论的微分波动方程单程波解法,这两类方法难以处理地震波横向速度变化剧烈的高陡倾角构造成像问题.近年来勘探地震学研究领域发展起来的叠前逆时偏移采用了双程波求解微分波动方程的算法,这种方法具有相位准确、不受介质横向速度变化和高陡倾角构造的影响、成像精度高、可以利用回转波正确成像等优点,从理论上弥补了当前工业界常规地震偏移所面临的成像缺陷.然而,叠前逆时偏移成像方法从理论走向实用尚需解决如下问题:计算速度和数据存储空间的节省、初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除等等.对于计算速度和存储量大的问题,随着计算机硬件的快速发展,将会不断得到改善,同时可以采取一些计算技术和存储策略来加以缓解.本文主要针对初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除这些因素,利用简单模型进行了分析.对于反射波造成的传播路径上的假像,给出了一种振幅补偿滤波方法.对勘探地球物理学界给出的SEG/EAGE二维盐丘模型、Marmousi模型和本研究设计的崎岖海底模型进行了叠前逆时偏移成像,均取得了较好的成像效果.     

波路径偏移压制层间多次波的理论与应用

消除层间多次波是地震勘探资料处理研究领域的难题,尤其对于实际资料的处理,到目前为止还很难找到一种完全有效的方法. 本文给出了仅对一次波成像既波路径偏移方法压制层间多次波方法,在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角,由此计算的角度作为射线追踪的初始角度,计算地震波射线的传播路径. 结合由程函方程计算的走时表,判断偏移范围是反射波还是多次波. 在前期偏移过程压制多次波的理论研究基础上,本文主要研究波路径偏移消除多次波的应用部分. 为了进一步说明效果的有效性,计算了在单炮和共成像点道集压制层间多次波,给出了实际资料的压制多次波的偏移结果.     

各向同性介质长偏移距地震同相轴动校正

传统二阶动校正方法基于较小最大偏移距与目标层深度比和地震波沿直线传播假设,进行长偏移距地震资料处理时,这些假设不再成立.高阶项动校正公式能提高长偏移动校正精度,文中对几种典型的高阶项动校正方法进行了比较,并提出了优化四阶、优化六阶动校正方法.模型计算表明,高阶项动校正方法能取得较常规动校正方法好的动校正结果,但并非阶数越高动校正精度就越高;在纵向速度变化剧烈时,高阶动校正或优化高阶动校正方法一般不能适用于最大偏移距与目标层深度大于3.5的地震反射同相轴,优化四阶和优化六阶动校正公式由于考虑了无穷大偏移距的影响,具有更稳定、更加精确动校正效果,适合于实际的各向同性长偏移距地震资料处理.     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分,但该方法会因为网格离散化而产生数值频散,降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散,然而纵、横波速度取值差异较大,在弹性波有限差分正演模拟中,在满足纵波最大速度确定的稳定性条件下,浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散,提高数值模拟精度,本文首先采用优化差分网格系数降低数值频散,然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制,但由于弹性波场中纵、横波速度差异明显,横波波场频散明显强于纵波,为了压制横波波场的数值频散,往往需要选取较大的频散压制参数,但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法,对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正,可以有效压制数值频散,削弱纵波FCT产生的假象.通过理论分析和数值算例发现,本文方法能有效削弱弹性波场有限差分数值...     

基于纵横波解耦的三维弹性波逆时偏移

纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性.     

Hausdorff分数维识别地震道初至走时

地震波初至走时的识别在地震勘探、人工地震层析成像以及全球地震层析成像方法研究中起重要作用．初至走时拾取的精度在很大的程度上影响地震层析成像及演的精度．本研究以提高地震波初至走时拾取的精度及定量化程度为目标,利用计算地震道时间序列分数维的方法,实现了地震波初至走时的自动拾取．本文以分形理论为基础,进行了地震道时间序列Hausdorff分数维的计算．计算结果表明地震道时间序列的分数维在初至到达前后具有不同的数值,其变化点能够定量指示出初至走时的位置．本文还给出了利用该方法对实测数据进行初至走时拾取的实例．