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

权利要求书

1.一种三维密度变速地震反射空校方法，其特征在于采用以下技术步骤实现：

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

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

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

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

5)对公共区域内的各地震测点下的各反射层进行密度校正；

6)根据新的层速度按照下式计算反射层厚度NewThick_i：

NewThick_i＝T_i*NewV_i/2000；

按照下式计算反射层界面埋深NewD_i：

NewD_i＝(NewThick₁+NewThick₂+…+NewThick_i)；

按照下式计算反射层的新的平均速度NewVa：

NewVa＝NewD_i*2/T_i；

7)按坐标得到每一地震层校正前后的反射层速度、平均速度和埋深，获得到经过三维密度变速地震反射空校的结果。

2.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤1)所述的三维剩余密度数据体为坐标X、坐标Y、坐标Z、剩余密度Δσ，坐标Z为地下半空间某一剖分长方体单元的中心海拔，剩余密度Δσ为该剖分单元的相对正负值，单位为g/cm³。

3.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤1)所述的地震速度数据体为坐标X、坐标Y、反射层名LName、反射层的反射时间T，单位为毫秒、反射层的层速度V、平均速度Va、反射层相对于统一基准面Height的底界埋深D。

4.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤2)所述的地震速度数据体计算是确定数据体的坐标点数，每个坐标点的坐标及该位置下面的地震反射层数目，循环地震速度数据体。

5.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤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。

6.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤3)所述的计算是首先确定数据体的记录数和求取数据体在空间三个方向上的循环数，再循环读取重力三维剩余密度体数据。

7.权利要求1的三维密度变速地震反射空校方法，其特征在于步骤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层剖分长方体单元的中心海拔。