[发明专利]一种黄土塬地震资料井控中长波长静校正方法

说明书

技术领域

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

背景技术

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

发明内容

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

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

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

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

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

所述步骤(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)井控中长波长静校正函数拟合；

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

通过选择拟合函数s(i,j)＝a(xx(i,j)-x₀)³+b(xx(i,j)-x₀)+c (2)

其中：s(i,j)为点(i,j)中长波长静校正量，

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

x₀为平滑后的地表高程的算术平均值，

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

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

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

根据拟合函数s(i,j)＝a(xx(i,j)-x₀)³+b(xx(i,j)-x₀)+c计算得到中长波长静校正量。

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

附图说明

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

图2是井控深度图。

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

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

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

具体实施方式

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

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

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

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

k＝Max_offset/Dcmp_inline，

l＝Max_offset/Dcmp_xline，

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