微屈多次反射与倾斜反射层

当进行地震反射采集的时候,无论是一次反射还是多次反射,在共深度点道集记录中均是一簇双曲线。这些双曲线在时间上的极小值对应于检波点与炮点相重合的那一点(即炮检距SGD=0处)。然而我们发现了这个规律的例外情形。如果反射界面之一是倾斜的,一部分微屈多次反射的时距曲线的双曲线顶点就会出现在炮检距SGD不为零的地方。而且当射线路径至少涉及一个倾斜界面时,本来在所有界面都是水平界面情形下具有相同旅行时间的数组微屈多次反射便分裂成互不相同的同相轴。每一波组具有不同的迭加视速度,有些比不是微屈多次反射的那些多次波的迭加速度低,有些则比它们高。一般来说,微屈多次波的射线离开震源时的射线角度与到达检波点时的射线角度不尽相同。我们导出了密度、层速度均不相等的一套地层的表达微屈多次反射时距关系的方程式。我们还讨论了具有不同层速度的模型,给出了由计算机产生的共深度点数据以及迭加之后的数据。由于每组微屈多次反射有几个不同的速度,因而使得用迭加来消除多次反射的方法受到影响。看来,当地下界面中包含倾斜反射层时,微屈多次反射可能是一个重要的噪声源。     

双检理论研究及合成处理

本文深入探讨了水、陆检波器的接收机制及波场动力学特征,认识到水层鸣震是由下行鬼波(虚反射)和上行微曲多次波组成。水检和陆检对于一次波、鬼波和微屈多次波具有不同的响应特征:对于一次波,双检的响应表现为同向;对于鬼波,双检的响应表现为反向;对于微曲多次波,双检的响应表现为同向。因此,本文通过水检和陆检同时采集的数据叠加去除鬼波;在此基础上,进一步求取海底反射系数,从而消除上行微屈多次波。本文方法在实际OBC资料处理中取得了良好应用效果。     

南黄海盆地滩浅海区多次波组合压制技术的应用

针对南黄海盆地滩浅海区地震资料中的鸣震干扰、全程多次波和层间多次波的特点,根据各类压制多次波方法的实用条件及其优缺点,提出了多次波压制策略:通过双检合并处理技术压制滩浅海区的鬼波和微屈多次波,通过预测反褶积压制残余鸣震及层间多次波;通过高分辨率Radon变换压制全程多次波和层间多次波;在叠前CRP道集上应用高分辨率Radon变换、内切和中值滤波压制小时差层间多次波。应用效果验证了该组合压制技术的有效性。     

三维层状介质中基于走时梯度的多次波射线追踪

本文结合分区多步思想改进初至波走时梯度射线追踪算法(MTG),使之能适用于三维层状介质中多次波射线追踪计算。在计算过程中按照射线类型和速度分界面将三维模型分为不同的计算区域,根据多次波的传播顺序从震源开始在对应的区域内采用网格剖分线性插值法计算波前走时场。然后从检波点开始根据多次波的类型在对应的计算区域内沿走时梯度方向分步追踪多次波射线路径。在计算多次波波前走时场的过程中,为提高计算效率加入了逐次网格剖分技术,使分界面上波前新的震源位置更加精确;在多次波射线追踪的过程中根据走时梯度方向追踪射线路径,在速度分界面区域内采用局部追踪策略求取下一个追踪点以提高多次波的射线精度。本文在多层层状介质模型及扩展到三维的Marmousi模型中进行试验,计算结果表明,改进的MTG算法结合分区多步技术计算多次波时具有常规MTG法计算初至波时的优点,计算精度高、算法简单、计算效率高,验证了MTG算法在复杂的三维模型中计算多次波的有效性。     

利用射线对复杂构造可定量成像吗?

只利用单次波(初至波或最强波)时,基于射线的叠前深度偏移不能对复杂构造定量成像。本文探索是否在基于射线保幅叠前深度偏移中考虑多次波可改善复杂介质的定量成像。 起初用于处理单次波的渐近Bron射线偏移反演,通过考虑所有震源和检波器射线路径的相互作用就可扩展到处理多次波。借助基于波前构造技术的动力学射线追踪来计算重叠射线场中的多次波。 采用复杂的Marmousi模型,我们证实即使模型复杂的深部Bron射线反演能提供一个可靠的相对波阻抗扰动的定量估计,而当只利用单次波(初至波或最强彼)时,Bron射线反演的振幅会被严重低估。 2D研究表明利用射线对复杂构造可进行定量成像。未来的研究需要在作3D之前,应对多次波Bron射线偏移或反演进行优化。     

利用射线对复杂构造可定量成像吗？

只利用单次波（初次波或最强波）时，基于射线的叠前深度偏移不能对复杂构造定量成像。本文探索是否在基于射线保幅叠前深度偏移中考虑多次波可改善复杂介质的定量成像。起初用于处理单次波的渐近Bron 射线偏移反演，通过考虑所有震源和检波器射线路径的相互作用就可扩展到处理多次波。借助基于波前构造技术的动力学射线追踪来计算重叠射线场中的多次波。采用复杂的Marmousi模型，我们证实即使模型复杂的深剖Bron射线反演能提供一个可靠的相对波阻抗扰动的定量估计，而 当只利用单次波（初至波或最强波）时，Bron射线反演的振幅会被严重 低估。2D研究表明利用射线对复杂构造可进行定量成像。未来的研究需要在作3D之前，应对多次波Bron射线偏移或反演进行优化。     

用x-t域预测反褶积法压制长周期多次波

用常规反褶积压制周期长于100毫秒的多次波时,有两种系统误差。一种误差是时窗效应,它是由根据数据的有限部分计算实际自相关函数而引起的。另一种误差是周期性假设所引起的,这种假设只有在一维介质、零偏条件下才能成立。误差的幅度随多次波周期长度的增加而增加。本文描述并证实了一种能克服常规一维法中两种系统误差的预测反褶积方程的精确二维解。该法可用于简单波场系统或捕波系统以及复杂或微屈多次波中。该方法并不要求将时窗设计成算子大小的6～10倍,在较大的传播角内均很精确。从经典预测意义上看,该法所用的公式是严密的。     

三江盆地低信噪比资料的多次波处理

三江盆地多次波形成的主要因素是区内存在一套反射很强的平行不整合地层(穆棱组),并且还发育大量高速煤层。本文从时空域、拉冬域、频率域及自相关等多个方面对多次波的地球物理特征进行了较详细的分析,发现了多次波在本区存在的规律性,并应用拉冬域、频率域正、反变换及自适应多次波切除等方法对多次波进行了分离和研究。表明拉冬域分离效果最好。在此工作的基础上,开发了一套多次波去除技术,由于该技术是在频率域应用拉冬正变换技术,提高了多次波的聚焦性能,并将滤波函数与相干技术相结合,从而实现自适应去除多次波的目的,而且稳定性很好。     

海上多次波的交互模拟

通过射线追踪计算模拟海底多次波、层间次波及全程多次波。计算模拟在工作站上通过交互菜单管理。针对南海珠江口盆地南部多次波严重地区的地质模型,进行了模拟试验。根据与实际的炮集记录和零偏移记录的对比可以判定测区深部的反射为有效波,而不是多次波。     

起伏层状TI介质中多次波射线追踪

现有的TI介质中的射线追踪算法很少涉及多次波的追踪计算,加之大多数算法又是建立在规则网格（或单元）模型划分的基础上,不利于追踪不规则界面下的多次波射线路径和相应走时。本文在模型参数化时引入不规则网格（或单元）划分,在主节点上对TI介质的6个参数（五个弹性参数,外加介质对称轴倾角）进行采样,计算得到三种类型 （qP,qSV,qSH） 波的群速度,通过单元主节点的群速度值线性插值得到次级节点上的群速度值。然后采用网格（或单元）波前扩展的方法进行分区波前扫描（射线追踪）,实现了起伏层状TI介质中多次（透射、反射、及转换）波的追踪计算。均匀各向异性介质中的计算结果表明,数值解与解析解吻合,验证了算法的有效性;模型数值模拟结果表明,该算法可以有效地模拟起伏地表和速度间断面,具有较高的计算精度和较快的CPU运行时间,能够适用于任意倾角的TI介质中的多次波射线追踪。