CN101609165B - 一种基于起伏地表的折射波静校正方法 - Google Patents

Abstract

本发明是物探地震资料处理过程中静校正的基于起伏地表的折射波静校正方法。在炮集上拾取来自同一折射层的初至时间，根据第一个折射层的旅行时公式，计算出炮点风化层厚度、风化层速度、第一个折射层速度、入射角和风化层的地层倾角模型参数，得到的单个折射层的近地表模型，计算静校正量；根据第二个折射层的折射波旅行时公式，求得第二个折射层的近地表模型，根据得到的两个折射层的近地表模型，计算静校正量。本发明直接求取风化层速度，用得到的近地表模型计算静校正量，对于近地表结构复杂的地区，可利用近地表模型作为初始模型进行层析反演，提高层析反演的运行效率。

Description

一种基于起伏地表的折射波静校正方法

技术领域

本发明涉及地球物理勘探技术，属于地震资料处理过程中静校正的技术范畴，是一种基于起伏地表的折射波静校正方法。

背景技术

目前，地震资料处理过程中静校正的方法很多，基本可以分为两类：一类是折射波法，另一类是层析反演法。

折射波法包括截距时间法、延迟时法、减去法、加减法、互换法、广义互换法和分解法等，这些方法都是根据折射初至时间和折射波旅行时公式直接计算高速层速度和延迟时的，利用延迟时求取风化层厚度时需要已知风化层的速度，一般风化层的速度需要应用小折射和微测井的资料求得。对于起伏地表的情况，一般都采用插值的方式，这样大大降低了近地表模型的精度，无法准确的得到风化层的速度。

另一类为层析反演方法，使用层析反演方法的关键问题有两个：一个是如何提高运行效率，另一个是提高反演法结果的精度。当初始模型与真实模型相近时，不但迭代次数少效率高，且反演结果的精度也高；如果初始模型与真实模型相差较远，不但迭代次数多效率低，且反演结果可能收敛到一个局部极值，反演结果的精度也较差。因此一个好的初始速度模型非常重要。如果不了解勘探地区信息，一般选择常速模型作为初始模型，但需要的迭代次数多计算时间较长并且反演结果的可靠性差并且运行速度很慢。

发明内容

本发明目的在于提供一种针对起伏地表小折射施工，运行效率高的基于起伏地表的折射波静校正方法。

本发明具体实施为如下步骤：

1)在起伏地表采集地震数据；

步骤1所述的地表采集是中间放炮，每炮都应该具有正、负排列。

步骤1所述的地表采集是先采集叠前炮集地震记录，然后拾取初至时间，并把初至时间记录到道头中。

2)用通常的方法在炮集上拾取来自同一折射层的初至时间；

3)对于第一个折射层的初至时间，求解折射波旅行时公式

步骤3)所述的计算折射波旅行公式为：

$T_{m,n}=\frac{{2h}_{m0}}{v_0}\cos i+\frac{H_n-H_m}{v_0}\cos (i-{\mathrm{\θ}}_1)+\frac{X_{\mathrm{mn}}}{v_0}\sin (i-{\mathrm{\θ}}_1)---\left(1\right)$

式中T_m，n表示炮点m、检波点n的折射波旅行时，h_m0为炮点m的风化层厚度，i为第一个折射界面入射角 $\sin i=\frac{v_0}{v_1},$ H_n，H_m表示检波点、炮点高程，X_mn表示炮检距，θ₁是折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为折射层速度；

步骤3)所述折射波旅行时公式解为(即根据初至时间和旅行时公式求解以下参数)：

$p=\frac{\sin (i-{\mathrm{\θ}}_1)}{v_0},$ $q=\frac{\cos (i-{\mathrm{\θ}}_1)}{v_0},$ ${\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

$\stackrel{\‾}{p}=\frac{\sin (i+{\mathrm{\θ}}_1)}{v_0},$ $\stackrel{\‾}{q}=\frac{\cos (i+{\mathrm{\θ}}_1)}{v_0},$ ${\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

4)将步骤3)中得到的解换算成炮点风化层厚度h_m0、风化层速度v₀、第一个折射层速度v₁、入射角i和风化层的地层倾角θ₁模型参数；

步骤4)所述的换算方法如下：

入射角 $i=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right),$

风化层速度为 $v_0=\frac{1}{\sqrt{p^2+q^2}},$

第一个折射层速度 $v_1=\frac{v_0}{\sin i},$

炮点风化层厚度 $h_{m0}=\frac{{\mathrm{\τ}}_m{v}_0}{\cos i},$

风化层的地层倾角 ${\mathrm{\θ}}_1=i-\mathrm{acr}\sin {\mathrm{pv}}_0=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right)-\mathrm{acr}\sin {\mathrm{pv}}_0;$

式中h_m0为炮点m的风化层厚度，i为第一个折射界面入射角 $\sin i=\frac{v_0}{v_1},$ θ₁是折射面视倾角，v₀为风化层速度，v₁为折射层速度；

5)在共检波点集合上重复步骤3)和步骤4)得检波点处的风化层厚度h_n0；

6)根据得到的单个折射层的近地表模型，计算静校正量；

7)对于第二个折射层的初至时间，求第二个折射层的折射波旅行时：

步骤7)所述的折射波旅行时用以下公式；

$T_{m,n}-T_{m,k}=\frac{h_{n0}-h_{k0}}{v_0}\cos \mathrm{\β}+\frac{(X_{\mathrm{mn}}-X_{\mathrm{mk}})\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$ (2’)

$+\frac{(H_n-H_k)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$

式中T_m，n表示炮点m、检波点n的折射波旅行时，T_m，k表示炮点m、检波点k的折射波旅行时，h_n0为检波点n处的风化层厚度，h_k0为检波点k处的风化层厚度，i为第二个折射界面的入射角 $\sin i=\frac{v_1}{v_2},$ H_n，H_k表示检波点n、检波点k的高程，X_mn表示炮点m到检波点n的炮检距，X_mk表示炮点m到检波点k的炮检距，θ₁是第一个折射面视倾角，θ₂是第二个折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n，k的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为第一个折射层速度，v₂为第二个折射层速度；

8)利用步骤7)求得的值，并结合步骤4)第一个折射层速度v₁，可得到第二个折射层速度v₂及第二个折射界面的倾角θ₂，再用步骤4)中得到的炮点m处的风化层厚度h_m0和检波点n处的风化层厚度h_n0，代入下公式，求得炮点与第一折射面垂足到第二折射面的距离h_m1；

$T_{m,n}=\frac{h_{m0}\cos \mathrm{\α}}{v_0}+\frac{h_{n0}\cos \mathrm{\β}}{v_0}+\frac{2h_{m1}\cos i}{v_1}+\frac{X_{\mathrm{mn}}\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$ (2)

$+\frac{(H_n-H_m)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$

式中T_m，n表示炮点m、检波点n的折射波旅行时，h_m0为炮点m处的风化层厚度，h_n0为检波点n处的风化层厚度，i为第二个折射界面的入射角 $\sin i=\frac{v_1}{v_2},$ H_m，H_n表示炮点m和检波点n的高程，X_mn表示炮点m到检波点n的炮检距，θ₁是第一个折射面视倾角，θ₂是第二个折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为第一个折射层速度，α为炮点第一折射界面入射角，β为检波点第一折射界面出射角，θ₂为第二个折射面视倾角；

9)在共检波点集合上重复步骤7)和步骤8)处理，得检波点处的第一个折射层厚度h_n1；

10)根据得到的两个折射层的近地表模型，计算静校正量。

本发明的优点是可以直接求取风化层速度，对于一般近地表结构不是很复杂的地区，直接用本发明得到的近地表模型计算静校正量，对于近地表结构比较复杂的地区，可以利用本发明得到的近地表模型作为初始模型进行层析反演，提高了层析反演的运行效率。

附图说明

图1是本发明步骤6)所得到的单个折射层的速度模型；

图2是应用图1速度模型计算的静校正量得到的叠后剖面；

图3是本发明步骤10)所得到的两个折射层的速度模型；

图4是应用图3速度模型计算的静校正量得到的叠后剖面。

具体实施方案

以下结合附图详细说明本发明。

本发明先采集叠前炮集记录，用常规方法拾取初至时间，通过初至时间计算出风化层的速度和厚度，以及两个折射层的速度和厚度。根据所求得的这些参数组成的近地表模型，进行静校正处理。

具体步骤如下：

1)在起伏地表采集地震数据，本发明在每公里大约有30米左右的高程差为起伏地表。本发明要求中间放炮，每炮都应该具有正、负排列。

2)在炮集上拾取初至时间，分析来自同一折射层的初至时间，一般来自同一折射层的初至时间通常呈线性关系。

对于近地表结构比较简单的情况，可以直接把初至时间认定为来自于一个折射层；当近地表结构复杂的时候，可以把初至分成来自两个折射层的初至时间。

对于只有一个折射层的近地表结构，仅使用第一个折射层公式求解，得到仅为单个折射界面的近地表结构；如果有两个折射层的近地表结构，分别用第一个折射层公式和第二个折射层公式求解，得到具有两个折射界面的近地表结构。

3)对于第一个折射层的初至时间，本发明用公式：

$T_{m,n}=\frac{{2h}_{m0}}{v_0}\cos i+\frac{H_n-H_m}{v_0}\cos (i-{\mathrm{\θ}}_1)+\frac{X_{\mathrm{mn}}}{v_0}\sin (i-{\mathrm{\θ}}_1)---\left(1\right)$

进行求解。

具体的求解过程为：公式(1)中已知为X_mn、H_m和H_n，未知数为h_m0、v₀、i和θ₁，对于这个公式，本发明设未知数为：

$p=\frac{\sin (i-{\mathrm{\θ}}_1)}{v_0},$ $q=\frac{\cos (i-{\mathrm{\θ}}_1)}{v_0},$ ${\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

对于负排列，地层倾角方向相反，公式(1)变为：

$T_{m,n}=\frac{{2h}_{m0}}{v_0}\cos i+\frac{H_n-H_m}{v_0}\cos (i+{\mathrm{\θ}}_1)+\frac{X_{\mathrm{mn}}}{v_0}\sin (i+{\mathrm{\θ}}_1)---\left(1^{,}\right)$

设负排列的未知数为：

$\stackrel{\‾}{p}=\frac{\sin (i+{\mathrm{\θ}}_1)}{v_0},$ $\stackrel{\‾}{q}=\frac{\cos (i+{\mathrm{\θ}}_1)}{v_0},$ ${\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

本发明用最小二乘法解这两个方程。

最小二乘法的实现过程(以正排列为例)：

T_m，n＝2τ_m+(H_n-H_m)q+X_mnp

最小二乘法的目标函数 $Q=\underset{n=1}{\overset{N}{\mathrm{\Σ}}}[T_{m,n}-2{\mathrm{\τ}}_m-(H_n-H_m)q-X_{\mathrm{mn}}p{]}^2\→\min$

N的取值通常根据排列中所含的同一折射层的初至时间来选取，N值越大拟和的效果越好，一般为6-20之间。

对Q求偏导展开，得到：

$\left[\begin{array}{ccc}2N& 2\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}& 2\underset{n=1}{\overset{N}{\mathrm{\Σ}}}(H_n-H_m)\\ 2\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}& \underset{n=1}{\overset{N}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}^2& \underset{n=1}{\overset{N}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}(H_n-H_m)\\ 2\underset{n=1}{\overset{n}{\mathrm{\Σ}}}(H_n-H_m)& \underset{n=1}{\overset{n}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}(H_n-H_m)& \underset{n=1}{\overset{N}{\mathrm{\Σ}}}{(H_n-H_m)}^2\end{array}\right]\left[\begin{array}{c}{\mathrm{\τ}}_m\\ p\\ q\end{array}\right]=\left[\begin{array}{c}2\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{T}_{m,n}\\ \underset{n=1}{\overset{N}{\mathrm{\Σ}}}{X}_{\mathrm{mn}}{T}_{m,n}\\ \underset{n=1}{\overset{N}{\mathrm{\Σ}}}(H_s-H_r)T_{m,n}\end{array}\right]$

解这个方程式，得到最佳的τ_m、p和q值。

同样的处理过程，可以得到负排列的最佳参数τ_m、p和q值，取两个τ_m的均值作为τ_m的值

4)用得到的τ_m、p、q、p、q值，求出未知数为h_m0、v₀、v₁、i和θ₁等参数。

求解的过程如下：

由于 $p=\frac{\sin (i-{\mathrm{\θ}}_1)}{v_0},$ $q=\frac{\cos (i-{\mathrm{\θ}}_1)}{v_0},$ ${\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0},$ $\stackrel{\‾}{p}=\frac{\sin (i+{\mathrm{\θ}}_1)}{v_0},$ $\stackrel{\‾}{q}=\frac{\cos (i+{\mathrm{\θ}}_1)}{v_0}$ 已经由步骤3)已经求出，可以求出 $v_0=\frac{1}{\sqrt{p^2+q^2}}$ 或者 $v_0=\frac{1}{\sqrt{{\stackrel{\‾}{p}}^2+{\stackrel{\‾}{q}}^2}},$ $i=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right),$ ${\mathrm{\θ}}_1=i-\mathrm{acr}\sin {\mathrm{pv}}_0=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right)-\mathrm{acr}\sin {\mathrm{pv}}_0,$ $v_1=\frac{v_0}{\sin i},$ $h_{m0}=\frac{{\mathrm{\τ}}_m{v}_0}{\cos i}.$

5)在共检波点集合上做同样处理(重复步骤3)和步骤4))，可得检波点处的风化层厚度h_n0

6)根据步骤1)-步骤5)得到的单个折射层的近地表模型计算静校正量，进行静校正处理。

静校正的实现过程为：一般是根据所得到的近地表模型，沿风化层底界面，把风化层剥掉(即利用近地表模型的风化层厚度和风化层速度计算出地震波在风化层中的传播时间)，然后用替换速度(一般给定恒定的替换速度，通常为高速顶界面速度的平均值)填充到统一基准面(即填充的过程为用统一基准面到风化层底界面的厚度和替换速度计算出射线的传播时间)。静校正的过程是减去风化层的传播时间，然后加上填充过程中射线传播时间。

7)对于第二个折射层的初至时间，本发明用公式：

$T_{m,n}=\frac{h_{m0}\cos \mathrm{\α}}{v_0}+\frac{h_{n0}\cos \mathrm{\β}}{v_0}+\frac{2h_{m1}\cos i}{v_1}+\frac{X_{\mathrm{mn}}\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$ (2)

$+\frac{(H_n-H_m)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$

来求解，第二个折射层还存在如下关系

$\frac{\sin \mathrm{\α}}{v_0}=\frac{\sin (i+{\mathrm{\θ}}_2-{\mathrm{\θ}}_1)}{v_1},$ $\frac{\sin \mathrm{\β}}{v_0}=\frac{\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$

本发明利用不同道之间相减消去前两项：

$T_{m,n}-T_{m,k}=\frac{h_{n0}-h_{k0}}{v_0}\cos \mathrm{\β}+\frac{(X_{\mathrm{mn}}-X_{\mathrm{mk}})\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$ (2’)

$+\frac{(H_n-H_k)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}$

用步骤3)所叙述的最小二乘法求解方程(2’)，求得

$p=\frac{\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1},$ $q=\frac{\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1},$ $\stackrel{\‾}{p}=\frac{\cos {\mathrm{\θ}}_1\sin (i+{\mathrm{\θ}}_2-{\mathrm{\θ}}_1)}{v_1},$

$\stackrel{\‾}{q}=\frac{\sin {\mathrm{\θ}}_1\sin (i+{\mathrm{\θ}}_2-{\mathrm{\θ}}_1)}{v_1},$ $r=\frac{\cos \mathrm{\β}}{v_0},$ $\stackrel{\‾}{r}=\frac{\cos \mathrm{\α}}{v_0}.$

8)用得到的r、p、q、r、p和q值，并结合第一折射层已求得的v₁，可求得v₂及θ₂，再用第一折射层的h_m0和h_n0，可从(2)式中求得h_m1。

求解的过程如下： ${\mathrm{\θ}}_1=\mathrm{arctg}\frac{p}{q}$ 或 ${\mathrm{\θ}}_1=\mathrm{arctg}\frac{\stackrel{\‾}{p}}{\stackrel{\‾}{q}},$ $v_0=\frac{1}{\sqrt{p^2+q^2+r^2}}$ 或者 $v_0=\frac{1}{\sqrt{{\stackrel{\‾}{p}}^2+{\stackrel{\‾}{q}}^2+{\stackrel{\‾}{r}}^2}},$ $\mathrm{\α}=\arccos \frac{\stackrel{\‾}{r}}{\sqrt{{\stackrel{\‾}{p}}^2+{\stackrel{\‾}{q}}^2+{\stackrel{\‾}{r}}^2}},$ $\mathrm{\β}=\arccos \frac{r}{\sqrt{p^2+q^2+r^2}},$ $i=\frac{1}{2}(\arcsin \frac{{\mathrm{qv}}_1}{\sin {\mathrm{\θ}}_1}+\arcsin \frac{{\stackrel{\‾}{q}v}_1}{\sin {\mathrm{\θ}}_1}),$ ${\mathrm{\θ}}_2=\frac{1}{2}(\arcsin \frac{{\mathrm{qv}}_1}{\sin {\mathrm{\θ}}_1}-\arcsin \frac{\stackrel{\‾}{q}{v}_1}{\sin {\mathrm{\θ}}_1})+{\mathrm{\θ}}_1,$ $v_2=\frac{v_1}{\sin i},$ 把所求得的参数代入公式(2)，便可以得到炮点与第一折射面垂足到第二折射面的距离h_m1。

9)在共检波点集合上做同样处理(重复步骤7)和步骤8))，可得检波点处的第一个折射层厚度h_n1。

10)根据得到的两个倾斜界面的近地表模型，做静校正处理。静校正过程如步骤6)所述。

本发明具体实施例：

1)采集叠前炮集地震记录，总炮数为392炮。

2)拾取初至时间。

3)用第一个折射层的初至时间，利用公式(1)求解未知数τ_m、p、q、p、q。

4)利用步骤3)求得的未知数，换算风化层和第一个折射层的模型参数，包括厚度、速度以及地层倾角等等。

5)在共检波点集合上做同样处理(重复步骤3)和步骤4))，可得检波点处的厚度。

6)根据步骤1)-步骤5)得到的单个折射层的近地表模型计算静校正量，应用到地震数据中，得到叠加剖面。图1是所得到的速度模型，图2是根据图1速度模型计算静校正量，然后把静校正量应用到地震记录上并叠加得到的地震剖面。

7)用第二个折射层的初至时间，利用公式(2’)和第一个折射层的模型参数，求解未知数r、p、q、r、p和q。

8)利用步骤7)求得的未知数，并结合第一折射层的近地表模型参数，求取第二个折射层的模型参数，包括厚度、速度以及地层倾角等等。

9)在共检波点集合上做同样处理(重复步骤7)和步骤8))，可得检波点的厚度。

10)用得到的两个倾斜界面的近地表速度模型，计算静校正量，应用到地震数据中，得到叠加剖面。图3是所得到的速度模型，图4是根据图3速度模型计算静校正量，然后把静校正量应用到地震记录上并叠加得到的地震剖面。

Claims (3)

1.一种基于起伏地表的折射波静校正方法，其特征是采用如下步骤：

1)在起伏地表采集地震数据；

2)用通常的方法在炮集上拾取来自同一折射层的初至时间；

3)对于第一个折射层的初至时间，求解折射波旅行时公式；

$T_{m,n}=\frac{{2h}_{m0}}{v_0}\cos i+\frac{H_n-H_m}{v_0}\cos (i-{\mathrm{\θ}}_1)+\frac{X_{\mathrm{mn}}}{v_0}\sin (i-{\mathrm{\θ}}_1)---\left(1\right)$

根据初至时间和旅行时公式求解以下参数：

$p=\frac{\sin (i-{\mathrm{\θ}}_1)}{v_0},q=\frac{\cos (i-{\mathrm{\θ}}_1)}{v_0},{\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

$\stackrel{\‾}{p}=\frac{\sin (i+{\mathrm{\θ}}_1)}{v_0},\stackrel{\‾}{q}=\frac{\cos (i+{\mathrm{\θ}}_1)}{v_0},{\mathrm{\τ}}_m=\frac{h_{m0}\cos i}{v_0}$

式中T_m，n表示炮点m、检波点n的折射波旅行时，h_m0为炮点m的风化层厚度，i为第一个折射界面入射角

H_n，H_m表示检波点、炮点高程，X_mn表示炮检距，θ₁是第一个折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为第一个折射层速度；

4)将步骤3)中得到的解换算成炮点风化层厚度h_m0、风化层速度v₀、第一个折射层速度v₁、入射角i和风化层的第一个折射面视倾角θ₁模型参数；

入射角 $i=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right),$

风化层速度为 $v_0=\frac{1}{\sqrt{p^2+q^2}},$

第一个折射层速度 $v_1=\frac{v_0}{\sin i},$

炮点风化层厚度 $h_{m0}=\frac{{\mathrm{\τ}}_m{v}_0}{\cos i},$

第一个折射面视倾角 ${\mathrm{\θ}}_1=i-\mathrm{acr}\sin {\mathrm{pv}}_0=\mathrm{arctg}\left(\frac{{\mathrm{pv}}_0+\stackrel{\‾}{p}{v}_0}{{\mathrm{qv}}_0+\stackrel{\‾}{q}{v}_0}\right)-\mathrm{acr}\sin {\mathrm{pv}}_0;$

式中h_m0为炮点m的风化层厚度，i为第一个折射界面入射角

θ₁是第一个折射面视倾角，v₀为风化层速度，v₁为第一个折射层速度；

5)在共检波点集合上重复步骤3)和步骤4)得检波点处的风化层厚度h_n0；

6)根据得到的单个折射层的近地表模型，计算静校正量；

7)对于第二个折射层的初至时间，求第二个折射层的折射波旅行时：

$\begin{array}{c}{T}_{m,n}-T_{m,k}=\frac{h_{n0}-h_{k0}}{v_0}\cos \mathrm{\β}+\frac{(x_{\mathrm{mn}}-X_{\mathrm{mk}})\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}\\ +\frac{(H_n-H_k)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}\end{array}---\left(2^{,}\right)$

式中T_m，n表示炮点m、检波点n的折射波旅行时，T_m，k表示炮点m、检波点k的折射波旅行时，h_n0为检波点n处的风化层厚度，h_k0为检波点k处的风化层厚度，i为第二个折射界面的入射角

H_n，H_k表示检波点n、检波点k的高程，X_mn表示炮点m到检波点n的炮检距，X_mk表示炮点m到检波点k的炮检距，θ₁是第一个折射面视倾角，θ₂是第二个折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n，k的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为第一个折射层速度，v₂为第二个折射层速度；

8)利用步骤7)求得的值，并结合步骤4)第一个折射层速度v₁，可得到第二个折射层速度v₂及第二个折射界面视倾角θ₂，再用步骤4)中得到的炮点m处的风化层厚度h_m0和检波点n处的风化层厚度h_n0，代入下公式，求得炮点与第一折射面垂足到第二折射面的距离h_m1；

$\begin{array}{c}{T}_{m,n}=\frac{h_{m0}\cos \mathrm{\α}}{v_0}+\frac{h_{n0}\cos \mathrm{\β}}{v_0}+\frac{{2h}_{m1}\cos i}{v_1}+\frac{x_{\mathrm{mn}}\cos {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}\\ +\frac{(H_n-H_m)\sin {\mathrm{\θ}}_1\sin (i-{\mathrm{\θ}}_2+{\mathrm{\θ}}_1)}{v_1}\end{array}---\left(2\right)$

式中T_m，n表示炮点m、检波点n的折射波旅行时，h_m0为炮点m处的风化层厚度，h_n0为检波点n处的风化层厚度，i为第二个折射界面的入射角

H_m，H_n表示炮点m和检波点n的高程，X_mn表示炮点m到检波点n的炮检距，θ₁是第一个折射面视倾角，θ₂是第二个折射面视倾角，其正负按以炮点到检波点方向的水平线为坐标轴，顺时针为负，逆时针为正，n的范围由来自同一折射层的炮检距范围确定，v₀为风化层速度，v₁为第一个折射层速度，α为炮点第一折射界面入射角，β为检波点第一折射界面出射角；

9)在共检波点集合上重复步骤7)和步骤8)处理，得检波点处的第一个折射层厚度h_n1；

10)根据得到的两个折射层的近地表模型，计算静校正量。

2.根据权利要求1所述的基于起伏地表的折射波静校正方法，其特征在于步骤1)所述的地表采集是中间放炮，每炮都应该具有正、负排列。

3.根据权利要求1所述的基于起伏地表的折射波静校正方法，其特征在于步骤1)所述的地表采集是先采集叠前炮集地震记录，然后拾取初至时间，并把初至时间记录到道头中。

Images