一种二维迴折射线积分法叠前深度偏移方法技术

一种二维迴折射线积分法叠前深度偏移方法，步骤为1)采用地震方法得到成像空间对应的深度域速度场；2)根据速度场划分矩形网格；3)在已经划分的矩形网格上进行射线追踪，而且对射线角度大于90度的射线继续追踪；记录射线在空间位置走过的时间，并将射线没有到达的地方进行插值得到充满空间的旅行时表；4)应用旅行时对叠前炮集数据进行叠前深度偏移处理，得到共成像点道集；5)将共成像点道集数据进行叠加得到该工区深度域的构造成像结果，并将结果显示为地层剖面图像。该方法射线追踪的精度高，为迴折射线旅行时产生创造条件，提高旅行时计算精度。

【技术实现步骤摘要】
【国外来华专利技术】
本专利技术涉及反射波地震数据处理过程中的二维叠前深度偏移构造成像技术，具体 的说是一种精确的二维迴折射线积分法叠前深度偏移构造成像方法。
技术介绍
随着油田勘探程度的提高，复杂地区和复杂油气藏的油气勘探逐步成为勘探的主 要目标。根据地表复杂度和地下构造复杂度可以分为以下四类：地表简单地下构造简单的 情况用叠后时间偏移就可以得到很好的成像效果；地表复杂地下构造简单的情况用叠前时 间偏移可以得到很好的成像效果；地表简单地下构造复杂的情况用叠后深度偏移可以得到 很好的成像效果；地表复杂地下构造复杂的情况用叠前深度偏移可以得到很好的成像效 果。在国内西北部有许多地区既是复杂地表又是地下构造复杂的情况，对于这种情况必须 利用叠前深度偏移加以解决。 叠前深度偏移单程波方法可以分为三大类：积分法、差分法和傅立叶变换法。三类 方法各有优缺点，其中积分法的优点是计算效率高、输入输出灵活和陡倾角成像效果好，缺 点是整体成像信噪比不如其它两类方法高；差分法的优点是成像精度高，缺点是计算效率 低；傅立叶变换法介于以上两种方法之间。需要注意的是叠前深度偏移的计算量是非常大 的，大的非常可怕。如果计算机的运算速度足够快，差分法叠前深度偏移是最好的选择，然 而目前计算机运算速度还不足以对大块的三维工区进行差分法叠前深度偏移，目前，国内 外进行大块工区叠前深度偏移的绝对主力是积分法，至少在未来的几年内是不会改变的。 常规的单程波叠前深度偏移方法能够得到地下复杂构造较好的成像结果，这里常 规的单程波叠前深度偏移方法指的是利用下行波场外推成像。但对于逆掩褶皱、逆冲断层 和倒转盐体等非常复杂的构造，侧翼成像还存在问题，尤其是大于90度的侧翼部分不能很 好的成像。解决这类问题需要利用迴折波成像，在高频近似情况下为迴折射线。所谓的迴 折波是没有经过任何地质体的反射而波变换传播方向，由向下传播变为向上传播再到达地 表的波（如图1 (a)所示）。正是迴折波可以到达盐丘侧翼大于90度的部分，再反射被地 面接受，如果能很好的利用迴折波，倒转盐体侧翼大于90度部分的构造成像就不再是问题 了。改善上述特定地质体的构造成像质量，需要对常规单程波叠前深度偏移方法进行改进。 将单程波有限差分法和傅立叶变换法改造成带有耦合项单程波方程算法可以对 构造倾角大于90度的构造成像（主要针对盐丘侧翼），当然这种方法的计算量是非常大的， 对于大块三维数据进行叠前深度偏移还不现实。如果能利用积分法叠前深度偏移计算效率 高的特点，将其改造成具有迴折波成像能力的技术对复杂区大规模叠前深度偏移的推广应 用有非常大的意义。目前，在国内还没有出现波积分法对构造倾角大于90度的构造进行成 像的技术，所以开展积分法迴折射线叠前深度偏移研宄势在必行，本专利技术在二维积分法叠 前深度偏移的基础上，开发出具有迴折波成像能力的积分法叠前深度偏移算法。
技术实现思路
本专利技术提供。 具体具体步骤包括： 1)采用地震的方法取得成像空间对应的深度域速度场； 2)根据速度场划分矩形网格； 3)在已经划分的矩形网格上进行射线追踪，而且对射线角度大于90度的射线继 续追踪；记录射线在空间位置走过的时间，并将射线没有到达的地方进行插值得到充满空 间的旅行时； 4)应用旅行时对叠前炮集数据进行叠前深度偏移处理，得到共成像点道集； 5)将共成像点道集数据进行叠加即可得到该工区深度域的构造成像结果，并将结 果显示为地层剖面图像。 本专利技术的二维迴折射线积分法叠前深度偏移方法，其【具体实施方式】为： 利用下面所示的射线追踪公式在已经划分的矩形网格上进行射线追踪： Kx4=f(x+At,x(t) +AtKx3,Px(t) +AtLx3) Kz4=f(x+At,z(t) +AtKz3,Pz(t) +AtLz3) Lx4= g〇+A t，x〇) + A tKx3，Px〇) + A t Lx3) Lz4=g〇+A t，z(t) + A tKz3，Pz〇) + A t Lz3) 其中，其中x、z是位置坐标分量，v是速度矢量，关于x、z的函数。匕、匕是慢度矢 量在x、z方向的分量，也是关于x、z的函数。 本专利技术能更好的参与旅行时空间插值，采用高阶龙格库塔公式使射线追踪的精度 更高，同时在单条射线追踪的过程中对射线方向不做限制，使迴折射线能很好的产生，本发 明在较大的射线阴影区不进行旅行时插值的方式代替阴影区填充的方式以提高旅行时计 算精度。【附图说明】 附图1表示迴折射线形成以及迴折射线反射示意图； 附图2表示V(z)介质纵波速度模型，即V(Z)介质的速度场； 附图3表示V(z)介质中某一点的射线追踪结果图； 附图4表示脉冲响应偏移测试结果，其中，附图4(a)为常规积分法脉冲响应，（b) 为迴折波偏移的脉冲响应测试结果； 附图5表示典型的倒转盐体模型纵波速度场； 附图6表示常规积分法叠前深度偏移结果； 附图7表示迴折射线积分法叠前深度偏移结果； 附图8表示常规积分法和迴折射线偏移结果之差； 附图9表示某地区实际资料的积分法叠前深度偏移结果，其中附图9 (a)表示常规 积分法偏移结果，（b)表示迴折射线偏移结果。【具体实施方式】 如图1所示，本专利技术的二维迴折射线叠前深度偏移方法基于射线方程，通过四阶 龙格库塔积分公式进行求解，得到每一条射线的射线路径和旅行时。传统的射线追踪方法 是逐个波前面向前追踪，并对每一条射线的传播方向做了不能超过90度的限制。当前第1页1 2 本文档来自技高网...

【技术保护点】
一种二维迴折射线积分法叠前深度偏移方法，其特征在于：1)采用地震方法得到成像空间对应的深度域速度场；2)根据速度场划分矩形网格；3)在已经划分的矩形网格上进行射线追踪，而且对射线角度大于90度的射线继续追踪；记录射线在空间位置走过的时间，并将射线没有到达的地方进行插值得到充满空间的旅行时表；4)应用旅行时对叠前炮集数据进行叠前深度偏移处理，得到共成像点道集；5)将共成像点道集数据进行叠加得到该工区深度域的构造成像结果，并将结果显示为地层剖面图像。

【技术特征摘要】
【国外来华专利技术】

【专利技术属性】
技术研发人员：杨顺伟，
申请(专利权)人：杨顺伟，
类型：发明
国别省市：广东;44