构造单程波算子的优化方法及其在偏移中的应用

波动方程叠前深度偏移具有处理复杂构造的特点,高精度偏移方法是勘探地球物理研究方向之一,其单程波算子的构造是波动方程叠前深度偏移的重要内容.本文着重对混合域中单程波算子构造的两种优化方法进行分析,讨论他们与其他单程波算子构造方法之间相互关系,并用上述方法计算变速介质的脉冲响应及偏移Marmousi模型数据,以此说明该类方法的偏移精度.     

频率-空间域有限差分法叠前深度偏移

为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果.     

面向目标的合成曲面波叠前深度偏移及其应用

面向目标的控制照明叠前深度偏移属于面炮偏移方法的一种特殊情况.通过对目标区域波场照射方向的控制,实现波场在目的层的最佳入射,提高目标区域成像质量.将目的层定义的目标合成算子外推到地表形成地表合成算子,应用地表合成算子进行波面源和记录的合成,最后通过地表合成算子的一系列旋转,实现旋转控制照明叠前深度偏移.本文系统阐述了面向目标的合成曲面波控制照明叠前深度偏移方法的基本原理和实现方法,通过数值模型试算验证了该方法的有效性,最后将其应用到某探区实际资料的处理中,并进行了实际资料的应用效果分析,总结了其在实际资料应用中的特点.     

地表旋转控制照明叠前深度偏移

基于炮集数据的常规叠前深度偏移波场外推的计算量非常大,控制照明叠前深度偏移是基于快速的组合炮技术,通过在目的层上控制震源波场,高效实现对地质目标的高质量成像. 组合炮所需的合成算子,由在目的层定义的震源波场反向外推到地表建立. 本文提出一种更快速的地表旋转控制照明方法,直接在地表对震源波场进行面向目标的照明控制,实现目标区域高质量成像. 该方法只需对某一核合成算子进行简单的几何旋转建立合成算子函数,避免了复杂的波场外推. 这种新的控制照明方法与快速的一般组合炮偏移效率相同,但对复杂目标区域的成像能力显著增强. 应用本文方法对Marmousi模型进行计算,得到深层复杂构造非常理想的成像结果.     

地震照明叠前深度偏移方法综述

对地震照明叠前深度偏移的基本概念、实现方法进行了分类和阐述.地震照明叠前深度偏移是通过对震源和炮记录进行合理的选择和合成,从而进行地震照明成像的一种有效方法.其可以分为平面波偏移和小束波偏移.理论模型的处理效果证明地震照明叠前深度偏移成像技术有很高的计算效率,并且还可以提高地下特定目标的成像质量.     

起伏地表条件下基于复Pade逼近的叠前深度偏移

叠前深度偏移是解决复杂地表和复杂构造地震成像的有效技术,而波场“直接下延”法实现了复杂地表条件下的地震成像.基于上述成果,结合高精度的波场延拓算子,本文提出了一种新的叠前深度偏移方法,这种方法是在波场延拓时,对声波方程中的平方根项进行复Pade逼近,通过推导得到基于复Pade逼近的傅里叶有限差分算子,结合波场“直接下延”法,实现了起伏地表条件下的叠前深度偏移,该算法减少了偏移噪音,从而得到准确、稳定的偏移成像结果.通过理论模型试算和实际资料试处理,验证了该方法的有效性.     

特征高斯波包叠前深度偏移方法

高斯波包（Gaussian packet）传播算子可在局部时空域高效地计算局部波包的传播.高斯波包叠前深度偏移的基础是在Gabor变换域描述观测数据,再利用高斯波包传播算子计算炮点波场和检波点波场,两者相关即可得到偏移结果.利用炮道集的局部τ-p特征在Gabor变换域表达观测数据,可以仅关注部分高斯波包框架函数上的数据投影,这样既实现了波场的压缩存储,同时可利用高斯波包传播算子反传框架函数以实现整个炮道集的快速反传.这些综合了观测数据局部τ-p特征的高斯波包函数称为特征高斯波包（characteristic Gaussian packet,CGP）,相应的波场反传称为特征高斯波包反传.理论及数值分析证明了上述特征高斯波包反传方法是有效且快速的.炮点正传波场也利用高斯波包传播算子模拟.利用互相关成像条件可实现特征高斯波包叠前深度偏移（characteristic Gaussian packet pre-stack depth migration,CGPM）.由于高斯波包传播算子描述了局部方向及局部空间的波场,所以CGPM可以自然地提取角度域成像道集（ADCIG）,并易于实现面向目标叠前深度偏移,从而作为偏移引擎为偏移速度分析（MVA）服务.数值实验证明了CGPM和面向目标CGPM的有效性和实用性.     

适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.     

面炮成像、控制照明与AVA道集

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点.面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题.本文采用面炮成像技术进行叠前深度偏移,通过对面炮震源下行波场的质量控制和优选射线个数和范围,以达到最佳的成像效果.采用控制照明技术,较大地提高了目标地层的成像精度.与此同时,得到振幅随入射角变化(AVA)道集,有利于叠前振幅解释和储层岩性预测.数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算.     

非稳态相移法叠前深度偏移

介绍一种能够适应介质速度横向变化的非稳态相移算子及其叠前深度偏移方法．为了克服常规相移偏移算法中要求速度横向不变的缺点,出现了基于非稳态滤波器理论的非稳态相移算子,即PSPI算子、NSPS算子和SNPS算子,其中SNPS算子是将前二者结合起来的一种对称的非稳态相移算子,它比前二者具有更高的精度和稳定性．为了提高运算速度,基于非稳态相移算子的叠前深度偏移算法采取了分片均匀近似的策略,Marmousi模型的叠前深度偏移结果证明了该算法的可行性和有效性。     

基于Chebyshev多项式的弯曲射线Kirchhoff叠前时间偏移

弯曲射线Kirchhoff叠前时间偏移被认为是一种综合了叠前时间偏移效率和叠前深度偏移精度的方法.本文以保精度地减少高阶Kirchhoff叠前时间偏移走时计算量为目标,在分析了Chebyshev正交多项式性质的基础上,建立了Chebyshev多项式约简系数表,进而用模拟退火法对转换系数进行分段优化,从而实现了在大炮检距...     

逆时偏移去噪方法研究进展

逆时偏移(RTM)作为精度较高的深度域成像方法之一越来越受到人们的重视..基于双程波动方程的逆时偏移具有原理简单、无倾角限制、可以用于复杂区域和陡倾角构造成像等优点.本文概述了国内外逆时偏移的发展历程和现状,介绍了应用随机边界条件后逆时偏移基本原理以及存在的问题并分析了成像噪声的产生机制,重点给出成像后压制逆时偏移低频噪声的几种方法及特征并给出总体认识.     

保幅Kirchhoff弯曲射线叠前时间偏移技术在天然气水合物三维地震资料处理中的应用

天然气水合物的富集往往与断裂、底辟及泥火山等构造有关,这就要求地震成像要精确,而针对水合物的地震处理又要以保真保幅为前提,因此快速高效而又有较高成像质量的的保幅Kirchhoff弯曲射线叠前时间偏移技术被广泛应用于三维水合物资料处理中。与直射线Kirchhoff叠前时间偏移技术相比,弯曲射线Kirchhoff叠前时间偏移同样具有快速高效的特点,同时成像精度在一定程度上可媲美叠前深度偏移。在实际资料的应用中可发现,基于保幅Kirchhoff弯曲射线叠前时间偏移技术处理的地震剖面可精确地刻画气体通道,有利于天然气水合物富集区的识别。     

τ-P域地震资料三维叠前深度偏移及其实用算法

本文提出了在τ-P域实现三维叠前深度偏移的计算方法，该方法把地震道集由时-空域变换到τ-P域，然后对共P数据体进行相移或相移加插值偏移，将偏移后的各共P数据体叠加得到叠前深度偏移结果．该算法能够实现陡倾地层的正确归位，降低偏移处理维数，减少计算机内存需求，易于并行处理，但是计算量和所花费的计算机时间仍然很大．为此我们又提出了针对目的层的分时（层）偏移方法，由于只考虑目的层，可节省大量的计算时间．文中以二维实际地震资料及三维数值模型资料的处理说明了上述方法的可行性与有效性．     

微机群并行实现Marmousi模型叠前深度偏移

波动方程叠前深度偏移等地震资料处理需要相当大规模的计算资源条件。巨型并行机因其计算能力强而备受工业界的青睐。但其价格的昂贵也限制了使用范围。近年来国际上出现的微机群是一种大规模并行计算的廉价实现方案，本文即是讨论利用微机群来实现Marmousi模型叠前深度偏移所涉及的数值实现。为比较起见，我们选择国际上大规模分布式并行计算机代表－IBM SP－2以进行计算性能的比较。     

基于共聚焦点道集的叠前深度偏移

共聚焦点(CFP)偏移技术是一种基于等时原理,将Kirchhoff积分法的一步偏移分两步聚焦（即激发聚焦和检波聚焦）来完成的叠前地震成像方法.该方法借助于逆时聚集算子和共聚焦点道集来实现叠前偏移成像.基于共聚焦点道集的叠前深度偏移是把基于共炮集的深度偏移的算法引入到CFP技术上来,基于波场延拓的理论来实现偏移成像,该方法首先生成共聚焦点道集,然后基于面炮合成的理论合成聚焦震源,最终通过相关成像来实现叠前偏移成像.该方法选取较少的聚焦点就可以实现对于地下构造的偏移成像,和炮域波动方程偏移相比,其计算效率得到了提高.通过模型试算和实际资料的试处理,验证了该方法在实现叠前深度偏移成像上的正确性和有效性.     

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

层析成像反演速度建模方法研究

本文应用层析成像反演方法优化叠前深度偏移速度模型。阐述叠前深度偏移速度建模流程,详细介绍利用两种不同层析反演方法来优化深度域的层速度模型,对同一实际地震资料速度模型建立的应用效果表明,在资料信噪比较高的地区能提供准确解释层位信息的前提下,模型层析反演优化后的速度模型精度要比网格层析反演优化后的速度模型更高,且稳定性强。