首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 570 毫秒
1.
单程波李代数深度积分的精度分析和算法改进   总被引:1,自引:1,他引:0       下载免费PDF全文
张廉萍  刘洪  李幼铭 《地球物理学报》2010,53(11):2739-2746
本文在利用拟微分算子和象征理论研究单程波传播规律的基础上,详细分析了单平方根算子李代数深度积分存在的问题,利用改进的算法对单平方根算子进行大步长延拓,计算其李代数积分和指数映射,得到地震波走时的解析表达式.对方法改进前后李代数积分象征中系数校正情况和线性横向变速介质中走时的计算误差进行了对比,数值算例表明,改进后的方法收敛更快,减少了计算时间且提高了计算精度.  相似文献   

2.
适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG/EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。  相似文献   

3.
三维各向异性介质中的波动方程叠前深度偏移方法   总被引:1,自引:0,他引:1       下载免费PDF全文
刘礼农  张剑锋 《地球物理学报》2011,54(11):2906-2915
基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀...  相似文献   

4.
双平方根单程波动方程叠前τ偏移方法   总被引:1,自引:0,他引:1       下载免费PDF全文
本文将常规双平方根(DSR)单程波动方程从深度域变换到双程垂直走时(τ)域,由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正,可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件,提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难,本文采用稳相近似,在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明,DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.  相似文献   

5.
本文将常规双平方根(DSR)单程波动方程从深度域变换到双程垂直走时(τ)域,由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正,可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件,提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难,本文采用稳相近似,在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明,DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.  相似文献   

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

7.
共炮检距道集波动方程保幅叠前深度偏移方法   总被引:9,自引:9,他引:0       下载免费PDF全文
本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括:(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率-波数域和频率-空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.  相似文献   

8.
柯璇  石颖 《地球物理学报》2017,60(11):4468-4479
通常工业界实现逆时偏移算法时采用有限差分数值方法模拟地震波场,波场模拟常常受稳定性条件限制,且易产生数值频散,成像精度降低.本文引入了一步法波场延拓方法,首先构建声波传播算子,借助Chebyshev多项式和Jacobi-Anger展开式近似传播算子中的e指数项,进而实现波场递推,该方法时间步长的选取不受稳定性条件限制而且不存在空间频散现象.本文将一步法波场延拓方法用于逆时偏移成像的波场模拟,并提出双缓冲区存储策略,在不增加计算量的前提下,大幅降低了逆时偏移方法的波场存储量.波场模拟和逆时偏移成像测试表明,本文提出的一步法波场延拓方法模拟地震波场精度高,消除了频散影响,可在较大时间步长的情况下实现高精度波场模拟;提出的基于一步法波场延拓的逆时偏移方法成像质量好;基于双缓冲区存储策略的逆时偏移成像方法存储成本低.  相似文献   

9.
单程波算子积分解的象征表示   总被引:7,自引:5,他引:7       下载免费PDF全文
单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法,他们代表了算子的不同表示方法,当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式,给出了算子象征的代数运算的头几阶表达式,这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式,对改善Kirchhoff积分的聚焦性能有重要意义.  相似文献   

10.
本文介绍的偏移成像方法为递归F-K法、相移插值法及串联偏移技术的有效结合,简称CFPI法.该方法中,波场外推采用相移插值法中递归向下延拓方式,而每一延拓层间的成像则用文中给出的常速F-K加剩余速度F-K偏移,通过与横向速度变化有关的插值来实现,该过程用FFT完成.该方法可采用较大的延拓步长速归向下外推与成像.理论分析及实际计算结果比较表明,该算法效率远高于相移插值法偏移.  相似文献   

11.
稳定的保幅高阶广义屏地震偏移成像方法研究   总被引:1,自引:0,他引:1       下载免费PDF全文
以先进的波动理论为基础的波动方程保幅地震偏移成像是在给出正确位置的同时也给出真实振幅的一种特殊完善.作者从保幅单程波动方程的非稳态相移公式出发,基于反问题求解中常用的摄动理论,利用单平方根算子的渐进展开,从而推导出保幅叠前深度偏移方程的高阶广义屏形式;针对散射波场计算项对于横向变速介质的不稳定性,通过数学近似提出一个有效提高稳定性的策略,应用到波场递归外推过程中,从而得到一种稳定的保幅高阶广义屏叠前深度偏移算子.理论模型试算和实际资料处理表明,该方法不但可以更精确地使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,使AVO响应更加清晰,提高了AVO资料的分析精度.  相似文献   

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

13.
波动方程深度偏移的频率相关变步长延拓方法   总被引:9,自引:1,他引:8       下载免费PDF全文
发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据(Marmousi模型)和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法.  相似文献   

14.
An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.  相似文献   

15.
在合成炮叠前深度偏移的基础上,提出了一种高效的合成炮叠前深度偏移方法,即在螺旋坐标下用因子分解进行波场外推的混合法. 用因子分解进行波场外推分为因果过程和反因果过程两个显式求解过程. 这种螺旋坐标下的显式求解过程,提高了波场外推的效率. 根据相位编码原理,对多个射线参数的合成波场进行了编码叠加,基于射线参数实现了MPI并行计算,进一步提高了计算效率. 在推导了有关公式并进行定量分析之后,对Marmousi复杂模型进行了计算,并作了比较,结果表明本文方法具有精度高和速度快的特点,可用于实际资料的计算.  相似文献   

16.
随着能源和资源勘查开采工作的深入,地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一,能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓,由于传统的地震波场延拓技术往往基于水平地表条件,相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法:基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法,采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造,获得了相应起伏地形条件下的复杂几何模型,开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性,可适应不同类型起伏地表采集的地震资料,显示出该方法在地震勘探领域的良好应用前景.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号