CN103293554A - 一种基于地震数据的静校正方法 - Google Patents

Abstract

本发明公开了一种基于地震数据的静校正方法，其具体是基于炮点井深建立表层模型的静校正方法。该方法利用地震资料采集时的炮点井深资料，插值得到所有检波点处的井深，该井深值为检波点处的虚拟的炮点井深。根据炮点和检波点高程，在给定基准面和替换速度的前提下，可计算炮点静校正量和检波点静校正量。其中检波点静校正量由两部分求和得到，一是检波点到潜水面静校正，二是潜水面到最终基准面静校正。应用炮点静校正量和检波点静校正量进行基准面校正，消除地表起伏和潜水面横向变化不均匀对地震记录产生的时差。通过本发明处理后的地震资料，可解决微测井资料少，信噪比低的地震资料的静校正问题。提高地震资料成像的信噪比和分辨率，有利于构造的准确落实。

Description

一种基于地震数据的静校正方法

技术领域

本发明涉及地球物理勘探领域，具体涉及一种基于地震数据的静校正方法，该方法特别适用为应用炮点井深建立表层模型的静校正方法。

背景技术

静校正是地震资料处理中的一个重要环节，静校正方法有高程静校正、折射波静校正、微测井静校正、组合静校正方法等。高程静校正对地表条件的假设是单一介质，该方法在实际应用中精度较低；折射波静校正虽然精度较高，但是该方法在地表高程变化较大的地震勘探区的应用，易出现空间低频假构造现象；微测井静校正直接应用微测井资料，可高精度求取井点处表层速度和厚度资料，但是空间上的插值精度直接影响该方法的静校正效果；目前工业界广泛应用的是组合静校正方法，即微测井静校正低频与折射波静校正高频进行组合的静校正方法。在表层调查点(微测井)较密的情况下，依据表层的微测井资料可较为准确的求取表层模型结构，进而可精确求取低频(长波长)静校正量，再和折射波静校正的高频量(短波长量)进行组合后，可得到精度很高的静校正量，应用该组合静校正量后，既可保持成像的宏观构造，又可提高地震资料的信噪比和分辨率。但是对于微测井资料稀疏或缺少的地区，目前的组合静校正方法就无法有效应用。

发明内容

为了解决现有技术的问题，本发明提供一种基于地震数据的静校正方法，本发明应用炮点井深建立表层模型的静校正方法，在微测井资料缺少的情况下，可求取准确的静校正量，提高地震资料的信噪比和分辨率。

依据本发明的技术方案，适用于基于炮点井深建立表层模型的静校正的基于地震数据的静校正方法，其采用以下步骤：

1)表层资料的获取，应用测量仪记录野外炮点的井深H_s0和高程H_s、检波点的高程H_r，炮点x坐标x_s，炮点y坐标y_s，检波点x坐标x_r，检波点y坐标y_r，并记录在文本文件上；

2)给定固定基准面高程H和替换速度值V；

3)计算任一检波点处虚拟的炮点井深H_r0，假设该检波点的坐标为(x_r，y_r)，该检波点的前一个炮点坐标为(x_s1，y_s1)，炮点井深为H_s01，后一个炮点坐标为(x_s2，y_s2)，炮点井深为H_s02，则该检波点处虚拟的炮点井深的计算公式可表示为

$H_{r0}=H_{s01}+\frac{\sqrt{{(y_r-y_{s1})}^2+{(x_r-x_{s1})}^2}}{\sqrt{{(y_{s2}-y_{s1})}^2+{(x_{s2}-x_{s1})}^2}}(H_{s02}-H_{s01});$

4)根据小折射资料的直达波时距曲线求取研究区低速层的速度V₀，

5)计算炮点到基准面的校正量t_s，

6)计算检波点到潜水面的静校正量t_r0，

7)计算潜水面到最终基准面静校正量t_r，

8)计算检波点到基准面的静校正量t_R，t_R＝t_r0+t_r；

9)基于步骤5)中获取的基准面的校正量t_s和步骤8)中获取的检波点到基准面的静校正量t_R，实现对地震资料的静校正。

其中，步骤9)为将t_S和t_R直接应用到地震资料中，实现静校正计算；或者分解得到t_S和t_R的低频成分，将得到的t_S和t_R的低频成分和折射波静校正的高频成分进行组合应用。

在上述技术方案中，步骤1)中所述的炮点井深H_s0的单位为m，炮点高程H_s的单位为m，检波点高程H_r的单位为m；步骤2)固定基准面高程H的单位为m，替换速度V的单位为m/s；步骤4)偏移距x的单位为m，旅行时t的单位为s，低速层速度V₀的单位为m/s；炮点基准面校正量t_s的单位为秒(s)；检波点到潜水面静校正量t_r0的单位为秒(s)；潜水面到最终基准面静校正量t_r的单位为秒(s)；检波点静校正量t_R的单位为秒(s)。

优选地，步骤3)的检波点处虚拟的炮点井深H_r0可应用该点处前后的炮点空间坐标和炮点井深通过线性插值求取。

优选地，步骤4)的低速层的速度V₀可应用表层调查的小折射方法求取，在小折射资料的时距曲线上求取直达波时距曲线的斜率(偏移距对时间导数)，该值即为低速层的速度值。

使用本发明的技术方案，应用炮点井深建立表层模型的静校正，重要的是静校正是地震资料处理的基础工作，在微测等表层调查资料缺少时，应用本发明方法进行静校正量的求取有效提高了地震资料的信噪比和分辨率，为勘探部署提供强有力的技术支撑。

附图说明

图1为本发明方法求取静校正量的原理示意图。

图2松辽盆地北部西坡地区东吐莫平洋工区未做静校正的地震剖面图。

图3松辽盆地北部西坡地区东吐莫平洋工区应用折射波静校正后的地震剖面图。

图4松辽盆地北部西坡地区东吐莫平洋工区应用本发明的方法和折射波组合静校正后的地震剖面图。

图5松辽盆地北部西坡地区巴彦查干工区110线应用折射波静校正的地震剖面。

图6松辽盆地北部西坡地区巴彦查干工区110线应用本发明的方法和折射波组合静校正后的地震剖面。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提供了基于地震数据的静校正方法，其适用于基于炮点井深建立表层模型的静校正，该方法采用以下步骤：

1)表层资料的获取，应用测量仪记录野外炮点的井深H_s0和高程H_s，检波点的高程H_r，并记录在文本文件上，同时记录测量点的大地坐标等参数；

2)根据工区情况和以往的处理经验，给定固定基准面高程H(如大庆油田给定H＝120)和替换速度值V(如大庆油田给定V＝1600m/s)；

3)计算任一检波点处虚拟的炮点井深 $H_{r0}=H_{s01}+\frac{\sqrt{{(y_r-y_{s1})}^2+{(x_r-x_{s1})}^2}}{\sqrt{{(y_{s2}-y_{s1})}^2+{(x_{s2}-x_{s1})}^2}}(H_{s02}-H_{s01});$

4)根据小折射资料求取研究区低速层的速度V₀；

5)计算炮点到基准面的校正量

6)计算检波点到潜水面的静校正量

7)计算潜水面到最终基准面静校正量

8)计算检波点到基准面的静校正量t_R＝t_r0+t_r；

9)静校正量的应用，一种方法是将t_S和t_R直接应用到地震资料中，实现静校正计算；另一种方法是分解得到t_S和t_R的低频成分，和折射波静校正的高频进行组合应用。

优选地，步骤3)的任一检波点处虚拟炮点井深可通过该点前后两个炮点的空间坐标和炮点井深通过空间显性插值求取。

优选地，步骤4)的低速层的速度V₀可应用表层调查的小折射方法求取。

下面结合附图，说明对松辽盆地西坡地区的二维地震资料进行静校正处理的具体实施方法：

1)表层资料的获取，应用测量仪记录野外炮点的井深H_s0和高程H_s，检波点的高程H_r，并记录在文本文件上，同时记录测量点的大地坐标等参数；

2)根据工区情况和以往的处理经验，给定固定基准面高程H＝120m和替换速度值V＝1600m/s；

3)根据炮点的空间坐标和炮点井深计算所有检波点处虚拟炮点井深H_r0；

4)根据小折射资料求取研究区低速层的速度V₀；

5)计算炮点到基准面的校正量t_S；

6)计算检波点到基准面的静校正量t_R；

7)分解本发明方法的静校正量的低频和折射波静校正的高频进行组合应用。

上述步骤1)至步骤6)所计算的基本原理图如图1所示。

参考图2，图2为松辽盆地北部西坡地区东吐莫平洋工区未做静校正的地震剖面图。由图可以看出该剖面的信噪比和分辨率都很低，特别是剖面中部几乎看不出同相。

参考图3，图3为松辽盆地北部西坡地区东吐莫平洋工区应用折射波静校正后的地震剖面图。从图可以看出，在应用静校正后，地震剖面的中部同相轴较为清晰，地震资料的信噪比和分辨率都有一定程度的提高。

参考图4，图4为松辽盆地北部西坡地区东吐莫平洋工区应用本发明的方法和折射波组合静校正后的地震剖面图。该图4和图3的信噪比和分辨率几乎相当，但是同相轴所反映的宏观构造和图3不一样，显示本发明的方法的低频成分的静校正量可以较为精确描述宏观构造的变化。也就是说图2未进行静校正的剖面和图3折射波静校正的剖面存在假构造，应用本发明方法重新处理后，有效消除了静校正不准的影响，有利于构造的准确落实。

参考图5，图5为松辽盆地北部西坡地区巴彦查干工区110线应用折射波静校正的地震剖面。可以看出黑色椭圆内的断层成像不太清晰，同时蓝色椭圆内的断层和同相轴成像效果不理想。

参考图6，图6松辽盆地北部西坡地区巴彦查干工区110线应用本发明的方法和折射波组合静校正后的地震剖面。和图5对比分析，显示应用本方法的方法处理后的地震剖面断点、断面清晰，能够满足断层的准确识别，且新处理成果剖面各反射层波组特征明显。

本发明主要实现缺少微测井资料的表层模型建立的方法，假设炮点深度都在潜水面(高速度层顶界)以下，和炮点相比，检波点大都埋在地表很浅的范围内，所以可根据炮点井深求取检波点位置的虚拟井深，这样就可以根据炮点井深度和虚拟井深的资料求取表层模型，在给定固定基准面和替换速度后，可实现基于表层模型的静校正计算。在后续的静校正应用中，可以单独应用炮点井深静校正量，也可以应用炮点井深静校正量的低频和折射波静校正量的高频进行组合的方法。

显然上述实施例仅为清楚的说明本发明所做的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上，还可以做出其他不同形式的变化或变动，这里无需也无法对所有实施方式予以穷举。由此所引申的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (6)

1.一种基于地震数据的静校正方法，其适用于基于炮点井深建立表层模型的静校正，该方法采用以下步骤：

1)表层资料的获取，应用测量仪记录野外炮点的井深H_s0和高程H_s、检波点的高程H_r，炮点x坐标x_s，炮点y坐标y_s，检波点x坐标x_r，检波点y坐标y_r，并记录在文本文件上；

2)给定固定基准面高程H和替换速度值V；

3)计算任一检波点处虚拟的炮点井深H_r0，假设该检波点的坐标为(x_r，y_r)，该检波点的前一个炮点坐标为(x_s1，y_s1)，炮点井深为H_s01，后一个炮点坐标为(x_s2，y_s2)，炮点井深为H_s02，则该检波点处虚拟的炮点井深的计算公式可表示为

$H_{r0}=H_{s01}+\frac{\sqrt{{(y_r-y_{s1})}^2+{(x_r-x_{s1})}^2}}{\sqrt{{(y_{s2}-y_{s1})}^2+{(x_{s2}-x_{s1})}^2}}(H_{s02}-H_{s01});$

4)根据小折射资料的直达波时距曲线求取研究区低速层的速度V₀，

5)计算炮点到基准面的校正量t_s，

6)计算检波点到潜水面的静校正量t_r0，

7)计算潜水面到最终基准面静校正量t_r，

8)计算检波点到基准面的静校正量t_R，t_R＝t_r0+t_r；

9)基于步骤5)中获取的基准面的校正量t_s和步骤8)中获取的检波点到基准面的静校正量t_R，实现对地震资料的静校正。

2.根据权利要求1所述的基于地震数据的静校正方法，其中步骤9)为将t_S和t_R直接应用到地震资料中，实现静校正计算。

3.根据权利要求1所述的基于地震数据的静校正方法，其中步骤9)为分解得到t_S和t_R的低频成分，将得到的t_S和t_R的低频成分和折射波静校正的高频成分进行组合应用。

4.根据权利要求1所述的基于地震数据的静校正方法，其特征在于，步骤1)中所述的炮点井深H_s0的单位为m，炮点高程H_s的单位为m，检波点高程H_r的单位为m；步骤2)固定基准面高程H的单位为m，替换速度V的单位为m/s；步骤4)偏移距x的单位为m，旅行时t的单位为s，低速层速度V₀的单位为m/s；炮点基准面校正量t_s的单位为秒(s)；检波点到潜水面静校正量t_r0的单位为秒(s)；潜水面到最终基准面静校正量t_r的单位为秒(s)；检波点静校正量t_R的单位为秒(s)。

5.根据权利要求2或权利要求3所述的基于地震数据的静校正方法，其特征在于，所述步骤3)的检波点处虚拟的炮点井深H_r0可应用该点处前后的炮点空间坐标和炮点井深通过线性插值求取。

6.根据权利要求2或权利要求3所述的基于地震数据的静校正方法，其特征在于，所述步骤4)的低速层的速度V₀可应用表层调查的小折射方法求取。

Images