一种黄土塬地震资料井控中长波长静校正方法

摘要

本发明提供了一种黄土塬地震资料井控中长波长静校正方法，属于石油地球物理勘探领域。本方法包括：(1)输入资料处理范围的地表高程文件，包括X、Y坐标及高程信息；(2)计算高程平滑参数k、l，然后对输入的地表高程文件进行平滑；(3)井控中长波长静校正函数拟合；(4)计算中长波长静校正量。本发明能够通过少量井的控制有效计算黄土塬地震资料中长波长静校正问题。

权利要求

1.一种黄土塬地震资料井控中长波长静校正方法，其特征在于：所述方法包括：(1)输入资料处理范围的地表高程文件，包括X、Y坐标及高程信息；

(2)计算高程平滑参数k、l，然后对输入的地表高程文件进行平滑；

(3)井控中长波长静校正函数拟合；

(4)计算中长波长静校正量；

其中，所述步骤(2)是这样实现的：①计算高程平滑参数k、l：

k＝Max_offset/Dcmp_inline，l＝Max_offset/Dcmp_xline，其中：k为沿测线方向的平滑点数，l为垂直测线方向平滑点数，

Max_offset为地震采集最大排列长度，Dcmp_inline为沿测线方向道间距离，Dcmp_xline为垂直测线方向道间距离；

②对输入的地表高程文件进行平滑：利用下式对地表高程文件进行平滑：其中：xx(i,j)为平滑后的点(i,j)的地表高程；

x(m,n)为点(m,n)的地表高程；

其中，所述步骤(3)是这样实现的：通过选择拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c   (2)其中：s(i,j)为点(i,j)中长波长静校正量，xx(i,j)为平滑后点(i,j)的地表高程，x0为平滑后的地表高程的算术平均值，a、b、c分别为相关系数，

通过3口井以上的数据，迭代求解得到系数a、b、c，然后剔除误差较多的井，重新计算求解系数a、b、c，最后分别求取a、b、c平均值。

2.根据权利要求1所述的黄土塬地震资料井控中长波长静校正方法，其特征在于：所述步骤(4)是这样实现的：根据拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c计算得到中长波长静校正量。

说明书

一种黄土塬地震资料井控中长波长静校正方法

技术领域

[0001] 本发明属于石油地球物理勘探领域，具体涉及一种黄土塬地震资料井控中长波长静校正方法。

背景技术

[0002] 鄂尔多斯盆地南部黄土塬地区地表条件沟壑纵横，地表起伏剧烈，地表巨厚的黄土产生严重的静校正问题，虽然在地震资料处理过程中采用一次静校正(层析静校正)和剩余静校正多次叠代，较好的解决了中波长和短波长静校正问题，但仍然存在较为严重的长波长静校正问题，处理成果存在与地表高程相关的假构造，影响了地震资料的构造解释。如何消除剩余中长波长静校正问题是黄土塬地震资料处理的一大难点，目前还没有有效的方法解决这个问题。

发明内容

[0003] 本发明的目的在于解决上述现有技术中存在的难题，提供一种黄土塬地震资料井控中长波长静校正方法，利用地表高程和关键控制井钻井分层数据，通过钻井分层和地震剖面同相轴时差，通过对地表高程进行平滑，然后通过控制井位置标志层时差，拟合中长波长静校正量函数，计算全区中长波长静校正量。消除与地表高程相关的假构造。

[0004] 本发明是通过以下技术方案实现的：

[0005] 一种黄土塬地震资料井控中长波长静校正方法，包括：

[0006] (1)输入资料处理范围的地表高程文件，包括X、Y坐标及高程信息；

[0007] (2)计算高程平滑参数k、l，然后对输入地表高程文件进行平滑；

[0008] 所述步骤(2)是这样实现的：

[0009] ③计算高程平滑参数k、l：

[0010] k＝Max_offset/Dcmp_inline，

[0011] l＝Max_offset/Dcmp_xline，

[0012] 其中：k为沿测线方向的平滑点数，

[0013] l为垂直测线方向平滑点数，

[0014] Max_offset为地震采集最大排列长度，

[0015] Dcmp_inline为沿测线方向道间距离

[0016] Dcmp_xline为垂直测线方向道间距离

[0017] ④对输入的地表高程文件进行平滑：

[0018] 利用下式对地表高程文件进行平滑：

[0019]

[0020] 其中：xx(i,j)为平滑后的点(i,j)的地表高程；

[0021] x(m,n)为点(m,n)的地表高程；

[0022] (3)井控中长波长静校正函数拟合；

[0023] 所述步骤(3)是这样实现的：

[0024] 通过选择拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c   (2)[0025] 其中：s(i,j)为点(i,j)中长波长静校正量，

[0026] xx(i,j)为平滑后点(i,j)的地表高程，

[0027] x0为平滑后的地表高程的算术平均值，

[0028] a、b、c分别为相关系数，

[0029] 通过3口井以上的数据，迭代求解得到系数a、b、c，然后剔除误差较多的井，重新计算求解系数a、b、c，最后分别求取a、b、c平均值。

[0030] (4)计算中长波长静校正量。

[0031] 根据拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c计算得到中长波长静校正量。

[0032] 与现有技术相比，本发明的有益效果是：本发明解决了与地表地形相关的假构造问题，效果非常明显。本发明能够通过少量井的控制有效计算黄土塬地震资料中长波长静校正问题。

附图说明

[0033] 图1是存在剩余中长波长静校正量的等T0图。

[0034] 图2是井控深度图。

[0035] 图3是井控剩余中长波长静校正处理后等T0图。

[0036] 图4是井控剩余中长波长静校正处理前后叠加效果对比。

[0037] 图5是本发明方法的步骤框图。

具体实施方式

[0038] 下面结合附图对本发明作进一步详细描述：

[0039] 剩余中长波长静校正量与地表高程具有相关性，本发明研制了一种井控中长波长静校正技术来解决黄土塬地区剩余中长波长静校正问题。

[0040] 由于三维地震资料采集，采用多条排列中间激发观测系统，受到排列长度的影响，存在长波长静校正问题，采用层析静校正处理技术，反演近地表模型计算一次静校正量，反演近地表模型时一般采用2500m以内的排列初至时间进行计算，以提高叠加成像效果为评判依据，中长波长静校正问题难以解决。

[0041] (1)计算高程平滑参数k、l

[0042] k＝Max_offset/Dcmp_inline，

[0043] l＝Max_offset/Dcmp_xline，

[0044] 其中：k为沿测线方向的平滑点数，

[0045] l为垂直测线方向平滑点数，

[0046] Max_offset为地震采集最大排列长度，

[0047] Dcmp_inline为沿测线方向道间距离

[0048] Dcmp_xline为垂直测线方向道间距离

[0049] (2)地表高程平滑

[0050] 通过研究，分析了剩余中长波长静校正问题和地表高程趋势有相关性，因此采用平滑地表高程，得到地表趋势，计算分析中长波长静校正量的合理性。

[0051]

[0052] 其中：xx(i,j)为平滑后点(i,j)的地表高程

[0053] x(m,n)为点(m,n)地表高程

[0054] k为沿测线方向的平滑点数

[0055] l为垂直测线方向平滑点数

[0056] (3)井控中长波长静校正函数拟合

[0057] 通过选择拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c   (2)[0058] 其中：s(i,j)为点(i,j)中长波长静校正量

[0059] xx(i,j)为平滑后点(i,j)的地表高程

[0060] x0为平滑后高程的算术的平均值

[0061] a、b、c分别为相关系数

[0062] 通过3口井以上的数据，就可以迭代求解系数a、b、c

[0063] 3口井以上的数据，可以分别迭代求解系数a、b、c后，剔除误差较多的井，重新计算求解系数a、b、c，最后分别求取a、b、c平均值。

[0064] (4)计算中长波长静校正量

[0065] 根据拟合函数s(i,j)＝a(xx(i,j)-x0)3+b(xx(i,j)-x0)+c计算中长波长静校正量。

[0066] (5)技术实现的步骤

[0067] 如图5所示，本发明完整的技术实现步骤如下：

[0068] ①输入需要处理的地表高程文件，包括X、Y坐标及高程信息。

[0069] ②按照以下公式计算平滑参数：

[0070] k＝Max_offset/Dcmp_inline

[0071] l＝Max_offset/Dcmp_xline

[0072] 其中：Max_offset为地震采集最大排列长度

[0073] k为沿测线方向的平滑点数

[0074] l为垂直测线方向平滑点数

[0075] Dcmp_inline为沿测线方向道间距离

[0076] Dcmp_xline为垂直测线方向道间距离

[0077] ③按照上文发明内容(2)中的公式(1)方法，对输入地表高程文件进行平滑；

[0078] ④在地表高程平滑的基础上，按照上文发明内容(3)以及公式(2)中的方法计算出相关系数a、b、c。

[0079] ⑤根据发明内容(4)计算得到中长波长静校正量。

[0080] 本发明是为了解决剩余中长波长静校正问题，这些问题反映在与地表高程具有相关性。本发明只要输入原始的地表高程文件，指定对应的平滑半径，程序就能自动输出高程平滑文件。再输入控制井位置及时差文件，就可以自动输出中长波长静校正量文件。

[0081] 图1是剩余中长波长静校正前的处理成果，得到的目的层等T0图，图2是利用钻井分层数据控制的深度图，图3是利用井控约束中长波长静校正处理后得到的目的层等T0图。对比图1与图2，构造差异较大，存在与地表具有相关性的假鼻隆构造，经过处理解释多次结合分析，认为还存在剩余中长波长静校正问题。对比图2和图3，等T0图与钻井分层数据控制的深度图趋势、形态较为一致，说明该方法较好的解决了黄土塬覆盖区严重的静校正问题。

[0082] 图4上部所示是地表高程，中部、下部所示分别是剩余中长波长静校正处理前后的叠加剖面。通过对比分析，可以看到井控剩余中长波长静校正处理后，剖面上对应高程较低位置的表现为“下凹”的假构造得到消除。

[0083] 本发明针对黄土塬地区存在地较为严重的中长波长静校正问题，通过输入地表高程

[0084] 和控制井位置的时间误差差，计算工区内的剩余中长波长静校正量，实际资料表明本发明较好的解决了与地表地形相关的假构造问题，效果非常明显。

[0085] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。