波动方程的高阶广义屏叠前深度偏移

不同于常规广义屏传播算子的推导中使用散射理论,本文利用单平方根算子的渐近展开,推导出了单程波方程广义屏传播算子的高阶表达式.高阶广义屏传播算子不仅可提高常规广义屏传播算子的计算精度,而且还能改善广义屏传播算子对速度强横向变化介质的适应性.把高阶广义屏传播算子应用于波动方程叠前深度偏移,可得到比常规广义屏传播算子更好的效果.高阶广义屏传播算子的阶数越高,计算精度越高,但计算量也越多.以SEG EAGE二维盐丘模型数据的波动方程叠前深度偏移为例,二阶广义屏传播算子相对于常规（一阶）广义屏传播算子增加了30%的计算量.高阶广义屏传播算子是常规广义屏传播算子理论的发展和完善.     

基于单平方根算子的叠前最小二乘分步傅里叶偏移方法

基于反演思想的最小二乘类偏移方法,具有较常规方法更好的成像保幅性和振幅均衡性,逐渐成为研究的热点.本文将单平方根分步傅里叶反偏移算子和偏移算子引入最小二乘框架,实现了基于单平方根算子的叠前最小二乘分步傅里叶偏移(SSR-LSSSF)算法.同时,为了抑制离散傅里叶变换引起的周期性伪波场,在偏移算法中引入了余弦衰减边界.在实现优化算法的基础上,通过对复杂盐丘模型的成像试处理及与其他方法成像对比可知:(1)本文发展的SSR-LSSSF偏移方法具有较好的迭代收敛特性,成像具有较好的保幅性,成像结果逐渐逼近真实反射系数.(2)通过与RTM、LSRTM及传统SSF偏移方法成像结果对比,表明本文方法在成像精度和偏移噪声压制方面具有一定的优势.     

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

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

面炮叠前深度偏移与控制照明技术

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

基于粒子群优化算法的高阶广义屏偏移成像（英文）

频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。     

基于Chebyshev多项式的非对称走时Kirchhoff叠前时间偏移角道集求取

地震射线走时的求取方法是叠前时间偏移研究的核心问题之一,也是影响计算时间域角道集角度精确性的关键因素之一.本文基于Kirchhoff叠前时间偏移,应用第一类切比雪夫多项式,对弯曲射线对称走时加以改进,引进非对称项,优化后得到切比雪夫非对称走时方程,与高精度走时进行比较和误差分析,再将该走时求取方法应用于时间域角道集的求取中,得到地下较真实的入射角.通过模型计算和实际地震资料处理证明,此种非对称走时及其角道集的求取方法具有精度高、计算量少的优点.     

单程波李代数深度积分的精度分析和算法改进

本文在利用拟微分算子和象征理论研究单程波传播规律的基础上,详细分析了单平方根算子李代数深度积分存在的问题,利用改进的算法对单平方根算子进行大步长延拓,计算其李代数积分和指数映射,得到地震波走时的解析表达式.对方法改进前后李代数积分象征中系数校正情况和线性横向变速介质中走时的计算误差进行了对比,数值算例表明,改进后的方法收敛更快,减少了计算时间且提高了计算精度.     

稳定的保幅高阶广义屏地震偏移成像方法研究

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

共炮检距道集波动方程保幅叠前深度偏移方法

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

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

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