[发明专利]三维密度变速地震反射空校方法

说明书

技术领域

本发明属于地球物理勘探方法，涉及重力-地震勘探的数据处理技术，具体是一种三维密度变速地震反射空校方法。

背景技术

地震勘探是利用地下地层层速度差异引起相应的反射波特征。现行地震勘探中反射层速度观测和研究是基于地震勘探方法建立的，即通过地震反射能量速度扫描等技术实现，因而在地震资料品质差的地区，无法进行有效的速度扫描，观测困难，甚至造成速度失真。在一些山区或构造复杂地区，地震勘探的反射层的速度受种种因素的影响，精度降低，导致地质模型的畸变，其结果与实际的差异较大，无法满足石油勘探对复杂地区的需要。

重力勘探是油气勘探中的另一种方法，它是针对地下地层密度差异引起对应的重力异常进行测量和分析，用于研究地质构造。

目前采用的三维密度分布是重力勘探获得的三维反演结果，比以往的二维反演技术更进一步。但在以往的地震速度勘探中很少利用重力勘探获得的三维反演结果，也没有利用三维密度对地震的反射层速度进行校正，而提高重力资料利用价值和提高地震勘探成功率的技术。

发明内容

本发明的目的在于提供一种重力与地震联合勘探，提高反射层速度的精度，解决复杂区实际地质问题的三维密度变速地震反射空校方法。

本发明采用以下技术步骤实现：

1)野外采集的重力数据和地震数据，分别处理得到三维剩余密度数据体和地震速度数据体；

步骤1)所述的三维剩余密度数据体为坐标X、坐标Y、坐标Z、剩余密度Δσ，坐标Z为地下半空间某一剖分长方体单元的中心海拔，剩余密度Δσ为该剖分单元的相对正负值，单位为g/cm³。

步骤1)所述的地震速度数据体为坐标X、坐标Y、反射层名LName、反射层的反射时间T，单位为毫秒、反射层的层速度V、平均速度Va、反射层相对于统一基准面Height的底界埋深D。

2)利用地震速度数据体数据计算各反射层的厚度和底界海拔埋深；

步骤2)所述的地震速度数据体计算是确定数据体的坐标点数，每个坐标点的坐标及该位置下面的地震反射层数目，循环地震速度数据体。

步骤2)按照下式计算各层的厚度Thick_i：

Thick_i＝T_i*V_i/2000；

按照下式计算各层的厚度相加得到各层的深度D_i：

D_i＝(Thick₁+Thick₂+…+Thick_i)；

按照下式计算各层的海拔埋深D_i：

D_i＝Height-D_i；

式中：地震各层的反射时间为T_i，层速度为V_i，统一基准面为Height。

3)利用重力三维剩余密度体数据计算厚度和底界深度；

步骤3)所述的计算是首先确定数据体的记录数和求取数据体在空间三个方向上的循环数，再循环读取重力三维剩余密度体数据。

步骤3)所述的厚度为剖分单元的厚度GThick，底界深度为该剖分单的底界海拔GD，按照下式计算：

GThick₁＝Z₁-Z₂，GD₁＝Z₁-(Z₁-Z₂)/2；GThick_i＝2*(GD_i-1-Z_i)，GD_i＝Z_i-GThick_i/2

式中Z₁，Z₂，……，Z_i，分别为三维剩余密度数据体中地下第一，第二，……，第i层剖分长方体单元的中心海拔。

4)对公共区域的重力三维剩余密度数据进行插值；

步骤4)所述的插值是依据地震速度数据体的坐标位置，对重力三维剩余密度数据体进行插值，计算出地震数据体所在位置的三维剩余密度。