单程波动方程叠前深度偏移并行算法 被引量：4

1

作者 单联瑜 徐兆涛 《石油地球物理勘探》 EI CSCD 北大核心 2008年第6期636-640,606+742,共7页

根据单程波动方程波场外推原理,文中总结了共炮道集数据混合域波场外推与成像算法流程,该流程主要包括地震数据预处理、波动方程速度分析和波场外推成像3部分;鉴于单程波动方程波场外推算法和数据流具有非常适合粗粒度的任务级并行计算... 根据单程波动方程波场外推原理,文中总结了共炮道集数据混合域波场外推与成像算法流程,该流程主要包括地震数据预处理、波动方程速度分析和波场外推成像3部分;鉴于单程波动方程波场外推算法和数据流具有非常适合粗粒度的任务级并行计算,处理器间的通讯较少,各个计算节点上可能有多个处理任务,且各节点的处理任务不等,负载极不平衡等特点,本文采用了主从模式动态负载平衡并行算法,该算法可以充分利用多个处理节点的空闲资源,提高作业的计算效率;将混合域单程波动方程波场外推算法与主从模式动态负载平衡并行算法相结合,实现了软件级断点恢复功能。应用国产"Stseis叠前深度成像处理软件系统"对胜利探区的实际资料进行了生产性应用处理,取得了较好的叠前成像结果。 展开更多

关键词 叠前深度偏移 单程波动方程 负载平衡 并行算法 断点恢复 成像

在线阅读 下载PDF 职称材料

一种单程波动方程的构造方法与应用

2

作者 柯本喜 戈良玉 +1 位作者 严昌言 范祯祥 《石油地球物理勘探》 EI CSCD 北大核心 1992年第3期318-330,共13页

本文通过引入一个与波的入射角度有关的参量,给出一种单程波动方程的构造方法。该方法的特点是有两个方向入射的平面波可满足所构造的方程,并可将该方程分裂为15°与45°形式的低阶方程。在判断波的传播方向时,只需判别出近法... 本文通过引入一个与波的入射角度有关的参量,给出一种单程波动方程的构造方法。该方法的特点是有两个方向入射的平面波可满足所构造的方程,并可将该方程分裂为15°与45°形式的低阶方程。在判断波的传播方向时,只需判别出近法线入射和近掠射两种波,给出θ=0,π/4、π/2三个值,即可有效地解决对入射角瞬时判别的不准确性及所引起的边界不适定性,同时又可有效地解决边界反射问题。文中以平面波源制作的叠加剖面和点源炮集记录的边界处理为例,说明该方法的处理效果是良好的。最后给出了构造的单程波动方程在开发常规WIMIG偏移模块中的应用。 展开更多

关键词 边界处理 偏移 单程波动方程

在线阅读 下载PDF 职称材料

数学检波器与波动方程地震叠前正演 被引量：13

3

作者 贺振华 胡光岷 黄德济 《成都理工大学学报（自然科学版）》 CAS CSCD 北大核心 2004年第6期675-678,共4页

地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义。叠前正演通常用射线追踪法和全程波动方程法来实现。这些方法在计算效率和记录信噪比等方面存在许多固有的问题。作者通过提出数学检波器(MG)的概念和方法,使利用单... 地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义。叠前正演通常用射线追踪法和全程波动方程法来实现。这些方法在计算效率和记录信噪比等方面存在许多固有的问题。作者通过提出数学检波器(MG)的概念和方法,使利用单程波动方程实现叠前正演成为可能。文中给出了用单程波动方程法对一个复杂地质模型形成8个共炮点地震记录的例子。在这些共炮记录中,来自复杂地质模型的反射波和绕射波十分清晰。记录质量好且计算效率高。更重要的是,利用作者研制的新方法可获得没有多次波、直达波和其它干扰的高信/噪比地震叠前正演记录。 展开更多

关键词 数学检波器 叠前正演模拟 单程波动方程 复杂地质模型

在线阅读 下载PDF 职称材料

波动方程保幅叠前深度偏移与AVO响应 被引量：3

4

作者 刘定进 印兴耀 +2 位作者 陆树勤 杨瑞娟 王丽 《中国石油大学学报（自然科学版）》 EI CAS CSCD 北大核心 2009年第4期45-51,共7页

从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的单程波动方程;按照叠前深度偏移方法中的波场双重向下延拓理论实现共炮域公式到CMP(共中心点)域公式的转换,推... 从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的单程波动方程;按照叠前深度偏移方法中的波场双重向下延拓理论实现共炮域公式到CMP(共中心点)域公式的转换,推导出由双平方根方程定义的保幅单程波动方程,从而将振幅误差补偿作为偏移的一部分在偏移过程中实现,使得地震成像在给出正确位置的同时也给出真实振幅。在此基础上获得的CRP(共反射点)道集在保证真实振幅的前提下,有效解决了绕射波干涉、反射点弥散、与倾角有关的NMO等问题,使AVO响应更加清晰,提高了AVO资料的分析精度。 展开更多

关键词 叠前深度偏移 保幅偏移 单程波动方程 双平方根 AVO

在线阅读 下载PDF 职称材料

基于惠更斯原理的波动方程共炮点道集地震正演 被引量：7

5

作者 吴永国 贺振华 黄德济 《成都理工大学学报（自然科学版）》 CAS CSCD 北大核心 2008年第3期274-278,共5页

地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义。作者利用下行波将震源场延拓到地质介质任意深度,实现单程波动方程模拟叠前正演记录。用单程波动方程法对一个盐丘地质模型生成129个共炮点地震记录,并对所有共炮记... 地震非零炮检距正演模拟技术在叠前地震资料处理和解释中有重要意义。作者利用下行波将震源场延拓到地质介质任意深度,实现单程波动方程模拟叠前正演记录。用单程波动方程法对一个盐丘地质模型生成129个共炮点地震记录,并对所有共炮记录进行叠前深度偏移。在共炮记录中,来自盐丘地质模型的反射波和绕射波十分清晰。叠前深度偏移剖面和盐丘地质模型完全吻合。该方法生成记录质量好且计算效率高,并且叠前正演记录没有多次波、直达波和其他干扰波。 展开更多

关键词 下行波 叠前正演 单程波动方程 波场延拓

在线阅读 下载PDF 职称材料

粘弹VTI介质单程波正演模拟与照明分析 被引量：6

6

作者 陈雪 韩立国 +1 位作者 杨贺龙 段超然 《石油物探》 EI CSCD 北大核心 2015年第6期643-651,共9页

实际地下介质不仅表现各向异性的特征,同时还表现粘弹性的特征,对入射波和反射波的传播及能量分布造成一定的影响。为了更加准确地分析地震波在实际地下介质中的传播规律,利用改进的单程波动方程数值模拟方法,对粘弹VTI介质模型进行qP... 实际地下介质不仅表现各向异性的特征,同时还表现粘弹性的特征,对入射波和反射波的传播及能量分布造成一定的影响。为了更加准确地分析地震波在实际地下介质中的传播规律,利用改进的单程波动方程数值模拟方法,对粘弹VTI介质模型进行qP波正演模拟以及点震源和不同角度入射平面波震源的照明度分析,并对正演及照明结果进行了对比分析。通过两组模型试验,总结了地震波在粘弹VTI介质中与在弹性各向同性介质中传播的不同之处:一方面,介质存在粘弹性会对在其中传播的地震波能量造成吸收;另一方面,介质还存在各向异性,导致地震波在介质中传播速度随方向变化,对地震波的能量分布造成影响。 展开更多

关键词 粘弹性 各向异性 单程波动方程 正演模拟 照明分析

在线阅读 下载PDF 职称材料

基于耦合反射/透射系数单程波传播算子的地震波模拟研究 被引量：7

7

作者 孙伟家 符力耘 姚振兴 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第10期2558-2565,共8页

单程波近似实际土是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,... 单程波近似实际土是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成:分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明:在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 展开更多

关键词 单程波动积分方程 分离变量逼近 单程波近似 反射/透射系数 地震模拟

在线阅读 下载PDF 职称材料

适合起伏地表的炮点波场延拓正演模拟 被引量：2

8

作者 杨海生 刘云生 《石油地球物理勘探》 EI CSCD 北大核心 2009年第2期146-151,共6页

本文提出一种波形相对保真、运算效率更高的适合起伏地表的炮点波场延拓正演模拟方法,炮点波场按波的传播方向延拓,运算中使用了适合起伏地表的波场延拓起始面技术,波形相对保真,无需做横向叠加。其基本思路是:①根据惠更斯原理,按地震... 本文提出一种波形相对保真、运算效率更高的适合起伏地表的炮点波场延拓正演模拟方法,炮点波场按波的传播方向延拓,运算中使用了适合起伏地表的波场延拓起始面技术,波形相对保真,无需做横向叠加。其基本思路是:①根据惠更斯原理,按地震波的传播方向,将炮点波场从起伏地表炮点深度处开始逐层向下延拓,与界面反射系数相关,直至最深处;②逐层向上延拓,并加上向下延拓时产生的新反射波炮点波场,直至起伏地表排列处,将炮点波场逐个与其检波点波场相关,经反傅里叶变换后可以获得一个排列的地震合成记录;③重复上述步骤,得到整个工区的所有炮记录。用文中方法对山峰体模型和高速透镜体模型的理论试算表明,地震合成记录的反射波和绕射波十分清晰,没有多次波、直达波和其他干扰波,信噪比高,所有炮记录的叠前深度偏移成像结果与地质模型完全吻合。 展开更多

关键词 起伏地表 炮点 波场延拓 正演模拟 单程波动方程傅里叶变换 叠前深度偏移

在线阅读 下载PDF 职称材料

基于有理式近似的高精度波场传播算子

9

作者 白敏 陈小宏 +1 位作者 熊晓军 吴娟 《石油地球物理勘探》 EI CSCD 北大核心 2013年第2期176-180,332+155,共5页

以地震波单程波传播方程为基础,利用有理式展开方法,推导出一种新的高精度波场传播算子。该算子波场延拓满足分裂步框架,采用快速傅里叶变换(FFT)进行计算。频散曲线分析以及数值试验的对比研究表明:本文方法的波场成像能力明显优于裂... 以地震波单程波传播方程为基础,利用有理式展开方法,推导出一种新的高精度波场传播算子。该算子波场延拓满足分裂步框架,采用快速傅里叶变换(FFT)进行计算。频散曲线分析以及数值试验的对比研究表明:本文方法的波场成像能力明显优于裂步傅里叶法(SSF),可以达到与傅里叶有限差分法(FFD)类似的波场成像效果。计算效率方面,其波场延拓只需三次FFT计算和一次插值,仅比SSF多了一次FFT和插值,但比FFD(计算量相当于五次FFT)少了两次FFT。这种算子不仅精度高,稳定性好,而且计算效率高,具有重要的理论意义和实际应用价值。 展开更多

关键词 单程波动方程 波场延拓 有理式近似 地震偏移 精度分析

在线阅读 下载PDF 职称材料

时-空局域化地震波传播方法:Dreamlet叠前深度偏移 被引量：3

10

作者 吴帮玉 吴如山 高静怀 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第9期3105-3114,共10页

提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(drea... 提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度. 展开更多

关键词 Dreamlet 时空局域化 多维局部余弦变换 单程波动方程 叠前深度偏移

在线阅读 下载PDF 职称材料

复杂海底自由表面多次波预测方法 被引量：14

11

作者 黎孝璋 邓勇 +2 位作者 赫建伟 任婷 顾汉明 《石油地球物理勘探》 EI CSCD 北大核心 2020年第1期64-70,6,共8页

针对崎岖海底多次波路径复杂、绕射波及绕射多次波发育等情形,自由表面多次波压制技术(SRME)预测多次波模型存在数据信噪比低、动力学和运动学信息欠准确等不利因素,提出了基于单程波及双程波波动方程延拓的多次波预测方法。即基于波动... 针对崎岖海底多次波路径复杂、绕射波及绕射多次波发育等情形,自由表面多次波压制技术(SRME)预测多次波模型存在数据信噪比低、动力学和运动学信息欠准确等不利因素,提出了基于单程波及双程波波动方程延拓的多次波预测方法。即基于波动方程并联合SRME技术,利用岩石物理模型对地震记录进行波场延拓,得到预测的多次波模型;同时采用自适应匹配相减方法提高复杂海底多次波压制效果。理论模型及实际数据试验表明:相对于单一的SRME多次波预测法,所提方法预测的多次波模型具有更高信噪比、与实际多次波在大多情况下的运动学和动力学特征吻合度更高。该方法理论上需已知准确岩石物理模型,由于海洋地震资料多次波主要是自由表面相关多次波,因此在仅已知海水及海底岩石物理参数情况下,该方法也能预测大部分与自由表面相关的多次波。 展开更多

关键词 复杂海底 多次波预测 SRME 双程波动方程 单程波动方程

在线阅读 下载PDF 职称材料

矩形网格三维显式差分偏移 被引量：1

12

作者 谢毅 Notfors Carl 《石油物探》 EI CSCD 2007年第6期588-593,共6页

由Hale提出的基于McClellan变换的三维显式差分偏移算法能够对速度横向变化剧烈的陡倾角复杂介质偏移成像,该算法计算效率高,但只适用于方形网格上采集的资料。对于矩形网格的地震记录,通常需将其先插值到方形网格上,然后再应用基于McCl... 由Hale提出的基于McClellan变换的三维显式差分偏移算法能够对速度横向变化剧烈的陡倾角复杂介质偏移成像,该算法计算效率高,但只适用于方形网格上采集的资料。对于矩形网格的地震记录,通常需将其先插值到方形网格上,然后再应用基于McClellan变换的三维显式差分偏移算法进行偏移成像,这样既降低了计算效率,又增加了对内存的需求。为此,提出了一种修正McClellan变换的三维显式差分偏移算法,即将插值嵌于外推算子中,无需对数据重采样,可以直接对矩形网格资料进行偏移成像。理论分析与数值实验结果表明,应用修正方法偏移成像的结果与基于数据插值的显式偏移算法结果相吻合,但计算量显著降低。 展开更多

关键词 显式差分 三维偏移 契比雪夫递归滤波器 单程波动方程 McClellan变换

在线阅读 下载PDF 职称材料

退化的Fourier偏移算子及其在复杂断块成像中的应用 被引量：19

13

作者 符力耘 孙伟家 李东平 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2007年第4期1241-1250,共10页

波动方程宽角抛物逼近得到的通常是非常系数的单程波传播算子,其系数是速度横向变化的函数,因此需要利用有限差分(FD)进行数值实施.通过对Lippmann-Schwinger单程波动积分方程的退化核逼近,本文研究了一类宽角退化算子的偏移成像.这种... 波动方程宽角抛物逼近得到的通常是非常系数的单程波传播算子,其系数是速度横向变化的函数,因此需要利用有限差分(FD)进行数值实施.通过对Lippmann-Schwinger单程波动积分方程的退化核逼近,本文研究了一类宽角退化算子的偏移成像.这种退化偏移算子只用快速Fourier变换进行波场延拓,将常规的Fourier分裂步地震偏移方法(SSF)推广适应强速度横向变化介质和大角度传播波场.退化的Fourier偏移算子通过在两个分裂步项之间作波数域线性插值来实现波场延拓,每延拓一层需要比常规的SSF地震偏移方法多一次快速Fourier变换(FFT).通过SEG/EAGE盐丘模型和实际地震资料的应用表明,退化Fourier偏移算子能很好地对盐下的陡倾角断层和实际地震剖面上的复杂小断块和大断裂地质构造成像. 展开更多

关键词 单程Lippmann—Schwinger波动积分方程 积分方程退化核逼近 退化算子 FOURIER 地震偏移 复杂断块成像

在线阅读 下载PDF 职称材料

GeoEast-Lightning叠前深度偏移软件在准噶尔盆地南缘的应用

14

作者 马晶晶 王艺豪 +2 位作者 罗勇 潘龙 李晓峰 《石油地球物理勘探》 EI CSCD 北大核心 2022年第S01期78-82,10,共6页

在简要介绍GeoEast-Lightning叠前深度偏移软件的Kirchhoff叠前深度(Q)偏移、高斯射束叠前深度偏移、单程波波动方程叠前深度(Q)偏移方法的基础上,对比分析了不同叠前深度偏移方法在准噶尔盆地南缘呼图壁三维地震工区的应用效果,以期为... 在简要介绍GeoEast-Lightning叠前深度偏移软件的Kirchhoff叠前深度(Q)偏移、高斯射束叠前深度偏移、单程波波动方程叠前深度(Q)偏移方法的基础上,对比分析了不同叠前深度偏移方法在准噶尔盆地南缘呼图壁三维地震工区的应用效果,以期为类似工区地震资料处理提供借鉴。结果表明:①Kirchhoff叠前深度(Q)偏移技术在实际生产中普遍应用,但它对复杂构造的成像能力及成像精度有限;②高斯束深度偏移克服积分法偏移精度不高的问题,保留了积分法偏移高效、灵活的特点,相比波动方程偏移计算效率大大提高;③单程波深度(Q)偏移方法耗时短、效率高,适应大多数具有一定倾角的地层成像。 展开更多

关键词 叠前深度偏移 KIRCHHOFF积分法 高斯射束偏移 单程波波动方程

在线阅读 下载PDF 职称材料