一种基于切比雪夫展开的广义屏叠前深度偏移方法技术

技术编号:7050157 阅读:298 留言:0更新日期:2012-04-11 18:40
为了提高强横向速度变化介质中波传播的描述精度,本发明专利技术提出了利用切比雪夫多项式展开求取高阶广义屏算子,从而改进广义屏叠前深度偏移方法,实现陡倾构造及横向速度变化剧烈地区的成像。

【技术实现步骤摘要】

本专利技术涉及地震资料处理领域,特别是涉及叠前深度偏移方法。
技术介绍
当前地震勘探的重点是构造复杂或速度变化剧烈的地区。在这些地区采用常规的偏移处理手段不能得到精确的地下构造成像,叠前深度偏移能适应构造复杂或速度变化剧烈区的精确成像,是复杂构造区精确成像最有效的方法。深度偏移技术始于70年代,80年代在理论和方法上有了较大的发展,90年代被广泛应用。叠前深度偏移中最主要的是波场外推算子的计算,目前常用计算波场外推算子的方法有基于射线理论的kirchhoff积分法、相位移法、空间一频率域有限差分法、分步 Fourier法、Fourier有限差分法及广义屏法。20世纪90年代以来,广泛应用于声学、电磁学等领域的波场传播相位屏法引入到反射地震学中,用于地震反射波场的波动方程偏移成像和波场模拟。相位屏法对高频或窄角度传播的波场描述比较准确,适用于弱散射和地下介质空间弱横向变化,同时还具有极高的计算效率和对介质性质空间变化的自适应性。但地震波的主频较低,在对空间强横向变化介质中地震波场传播进行描述时,该法会出现较大的误差。为了准确的描述强横向变化介质中地震波的传播,人们在相位屏理论的基础上,发展了广义屏法。广义屏法不仅保持了裂步Rmrier法中相空间坐标变量间的分离避免了物理本空间的坐标变量与其对偶空间坐标变量间的耦合,使之具有高效与易于实现的特点,而且对于宽角度传播的地震波具有相当高的描述精度,对速度强横向变化介质也具有很好的适应性。当前的波场传播中的散射波场的Born近似和Rytov近似的广义屏算子都是一阶近似传播算子。上述波场传播算子对强横向变化介质中波传播的描述精度较低。为了提高波场传播的偏移角度,本专利技术通过切比雪夫多项式对单程波方程的单平方频散方程进行多项式展开,推导出广义屏传播算子的高阶表达式。
技术实现思路
为更有效地克服相位屏叠前深度偏移中存在的上述缺陷,本专利技术目的是通过对频散方程进行高阶展开,从而提高对地震波的描述精度。本专利技术的目的是提供一种在提高广义屏波场传播算子阶次的基础上,来提高对强横向变化介质中波传播的描述精度。该方法步骤如下首先,介质速度为V的精确上行波的频散关系式为本文档来自技高网
...

【技术保护点】
1.一种基于切比雪夫展开的广义屏叠前深度偏移方法,该方法包括:步骤1对频率—波数域的上行波的精确方程进行多项式展开,利用切比雪夫多项式展开求取展开多项式的系数,使其展开的多项式式与精确的上行波的频散方程的逼近程度最高,求取展开多项式的各项系数;步骤2 求出切比雪夫展开多项式系数,如式(1)所示:    (1)式(1)中:,为角频率(Hz),为方向的波数(/m), 为介质速度(m/s),为介质速度(m/s),为展开的阶数,为切比雪夫系数,为切比雪夫多项式;步骤3 将式(1)代入波场外推方程并利用指数的一阶近似,得切比雪夫多项式展开的波场外推方程;步骤4 读入二维SEG/EGEA模型数据,对数据进行傅里叶分析;步骤5 在频率—波数域中,在每一个延拓步长上,采用切比雪夫多项式展开的高阶广义屏算子向下延拓波场;步骤6延拓过程中下一个深度的波场为上一个深度波场延拓后的结果,将延拓后的结果进行反傅里叶变换到频率—空间域进行叠加成像,输出成像结果。

【技术特征摘要】
1.一种基于切比雪夫展开的广义屏叠前深度偏移方法,该方法包括步骤1对频率一波数域的上行波的精确方程进行多项式展开,利用切比雪夫多项式展开求取展开多项式的系数,使其展开的多项式式与精确的上行波的频散方程的逼近程度最高,求取展开多项式的各项系数;步骤2求出切比雪夫展开多项式系数,如式(1)所示/(J^ = 1.00000154 0J5167624c 01401117V + 0.0575672ε5 + 0.069I4592r4- ⑴-0.1964S4&C5 -0.17172 2&fii- m1 I I式(1)中1=),为角频率(Hz), 为方向的波数(/m), 为介质速、C V 极ZC度(m/S),v为介质速度(m/s),;i为展开的阶数,Cf为切比雪夫系数,:ζ为切...

【专利技术属性】
技术研发人员:罗仁泽黄元溢
申请(专利权)人:西南石油大学
类型:发明
国别省市:90

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1