CN102901985B - 一种适用于起伏地表的深度域层速度修正方法 - Google Patents

Abstract

本发明提供了一种适用于起伏地表的深度域层速度修正方法，属于石油地球物理勘探的地震资料处理领域。所述方法首先对偏移距分段，将输入道集按照偏移距进行分组，并进行起伏地表叠前深度偏移，产生共成像点道集；然后进行起伏地表剩余速度分析，计算高阶剩余曲率谱，并根据高阶剩余曲率谱的局部最大值拾取深度‑剩余曲率对，然后求取速度模型修正量；再对速度模型修正量进行空间插值和平滑，并将速度模型修正量跟初始速度模型相加得到更新后的起伏地表深度域层速度模型；最后利用更新后的起伏地表深度域层速度模型代替初始速度模型，重新进行起伏地表叠前深度偏移和起伏地表剩余速度分析，从而不断提高速度模型的精度。

Description

一种适用于起伏地表的深度域层速度修正方法

技术领域

本发明属于石油地球物理勘探的地震资料处理领域，具体涉及一种适用于起伏地表的深度域层速度修正方法，用于地震勘探中深度域层速度的建模和模型修正更新中。

背景技术

叠前深度偏移是现阶段非常重要的地震资料处理技术之一，是解决复杂地质构造成像的重要方法，它能直接提供深度域的地震构造成像资料，在发现规模油气田中具有举足轻重的作用。在叠前深度偏移中深度域层速度是至关重要的影响因素，速度模型的精度直接决定了成像结果的质量。高精度的层速度来源于对初始层速度的剩余速度分析，所谓的剩余速度分析是指在已有的层速度基础上进一步提高速度精度的方法。

剩余速度分析一般利用偏移后的共成像点道集，叠前深度偏移生成的共成像点道集中蕴含了剩余速度的信息，当深度域层速度不正确时，共成像点道集的同相轴会表现出一定的剩余曲率。基于剩余曲率的剩余速度分析方法可以从中提取这一信息。在水平地表情况下，Al-yahya(1989)给出了地下水平反射层、横向速度均匀假设下的剩余时差曲线公式，Lee和Zhang(1992)给出了带倾角的炮偏移共接收点道集的剩余时差曲线公式。成谷等(2007)将扫描曲线的振幅绝对值作为参数改进了剩余曲率谱的计算方法，提高了剩余曲率谱的分辨率。

复杂地表情况下的勘探问题是近年来研究的重点，探索起伏地表情况下深度域层速度模型的更新技术，对于提高复杂地表情况下叠前深度偏移的成像精度具有重要意义，而现有的基于共成像点道集的剩余曲率速度模型更新方法，是在地表是水平面的假设下推导出来的，没有涉及地表起伏因素。

发明内容

本发明的目的是采用起伏地表剩余速度分析方法提高剩余曲率扫描曲线的拟合精度，采用高阶剩余曲率谱计算方法提高分辨率，采用优化的剩余速度求取方法提高剩余速度精度，进而提供一种适用于起伏地表叠前深度偏移的深度域层速度修正方法。

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

一种适用于起伏地表的深度域层速度修正方法，所述方法首先对偏移距分段，将输入道集按照偏移距进行分组，并进行起伏地表叠前深度偏移，产生共成像点道集；然后进行起伏地表剩余速度分析，计算高阶剩余曲率谱，并根据高阶剩余曲率谱的局部最大值拾取深度-剩余曲率对，然后求取速度模型修正量；再对速度模型修正量进行空间插值和平滑，并将速度模型修正量跟初始速度模型相加得到更新后的起伏地表深度域层速度模型；最后利用更新后的起伏地表深度域层速度模型代替初始速度模型，重新进行起伏地表叠前深度偏移和起伏地表剩余速度分析，从而不断提高速度模型的精度。

所述方法包括以下步骤：

(1)将偏移距分段，并将输入道集按照偏移距进行分组；

(2)利用初始速度模型对各段偏移距道集进行起伏地表叠前深度偏移，产生共成像点道集；产生的共成像点道集按照偏移距排列，每个共成像点道集的道数与步骤(1)中的偏移距分段数相同；

(3)判断共成像点道集是否拉平，如果已经拉平，则转入步骤(7)，如果没有拉平，则转入步骤(4)；

(4)计算起伏地表剩余曲率谱：首先读入每个CDP处偏移产生的共成像点道集，然后利用高阶剩余曲率谱计算公式计算起伏地表剩余曲率谱，并根据局部最大能量值拾取深度-剩余曲率对；

(5)利用剩余速度求取公式对每一个共成像点道集的深度-剩余曲率对计算深度域层速度修正量；

(6)对步骤(5)获得的所有深度域层速度修正量进行空间插值和平滑处理，获得整个空间的速度模型修正量，并将该速度模型修正量与初始速度模型相加，得到更新后的速度模型；转入步骤(2)，同时用更新后的速度模型代替步骤(2)中的初始速度模型；

(7)输出偏移剖面。

在步骤(1)中根据偏移孔径将偏移距分为20-35段，每段包含至少一道；各偏移距段取等间距排列；偏移距的段数和间隔在偏移过程中给出。

所述步骤(4)中的要计算的CDP位置是人为给定的。

所述更新后的速度模型可直接用于地震处理中的叠前深度偏移偏移，产生的共成像点道集又可以不断更新地震勘探区的深度层速度模型。

所述高阶剩余曲率谱计算公式为：

$\mathrm{Spe}(z,\mathrm{\γ})=\frac{{\left[\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}U(z_i,i)\right]}^{2n}}{\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}{U}^{2n}(z_i,i)}---\left(13\right)$

其中，U(z_i，i)表示当前γ值对应的剩余曲率扫描曲线所经过共成像点道集网格(z_i，i)处的振幅值，i是共成像点道集内的道数，代表了偏移距信息，z_i是当前γ值对应的剩余曲率扫描曲线在第i道的偏移深度，由起伏地表剩余速度分析公式求得，Ntr表示一个共成像点道集的总道数；n大于1即为高阶谱计算公式，n等于1为常规公式；

所述起伏地表剩余速度分析公式为：

$z_m=({\mathrm{dz}}_1+{\mathrm{dz}}_2)/2+\sqrt{\frac{L_{\mathrm{offset}}^2}{4}({\mathrm{\γ}}^2-1)+\frac{{(2z-{\mathrm{dz}}_1-{\mathrm{dz}}_2)}^2}{4}{\mathrm{\γ}}^2}---\left(10\right)$

其中，dz₁表示地表CDP点C与炮点S的高程差，dz₂表示地表CDP点C与检波点R的高程差，z表示地表CDP点C与地下反射点P之间的距离，z_m是深度域层速度为v_m时对应的偏移深度，L_offset是炮点和检波点之间的偏移距，γ＝v_m/v，表示叠前深度偏移所用的深度域层速度与准确层速度之间的比值，即剩余曲率值。

求取所述速度模型修正量的公式为：

${\mathrm{deltV}}_1={\stackrel{\‾}{v}}_1/{\mathrm{\γ}}_1-{\stackrel{\‾}{v}}_1---\left(15\right)$

${\mathrm{deltV}}_i=\frac{n_i-n_{i-1}}{\frac{n_i}{{\stackrel{\‾}{v}}_i/{\mathrm{\γ}}_i}-\frac{n_{i-1}}{{\stackrel{\‾}{v}}_{i-1}/{\mathrm{\γ}}_{i-1}}}-\frac{n_i-n_{i-1}}{\frac{n_i}{{\stackrel{\‾}{v}}_i}-\frac{n_{i-1}}{{\stackrel{\‾}{v}}_{i-1}}}---\left(16\right)$

(15)式表示第一个拾取点处的速度修正量，(16)式表示当i大于1时，第i个拾取点处的速度修正量；

其中， 表示第i个拾取点处的平均速度，n_i代表从地表到第i个拾取点处的深度样点数；

γ＝v_m/v，表示叠前深度偏移所用的深度域层速度与准确层速度之间的比值，即剩余曲率值。

与现有技术相比，本发明的有益效果是：

(1)本发明利用精度较高的起伏地表剩余速度分析公式和剩余曲率谱计算公式进行剩余速度谱的计算，该公式结合了地表炮检点高程，能够处理起伏地表叠前深度偏移产生的共成像点道集；同时可以向下兼容平地表公式；成像道集同相轴拉平时，成像深度为真实的深度。

(2)应用本发明可以得到精度更高的深度域起伏地表剩余速度模型和起伏地表速度更新模型，应用更新的深度域速度模型进行叠前深度偏移可以得到更加准确的构造地震成像结果。

(3)本发明对剩余速度修正量的插值和平滑等措施都为得到准确度较高的剩余速度打下良好的基础。

附图说明

图1是起伏地表射线路径示意图。

图2-1是利用本发明方法获得的起伏地表剩余曲率扫描曲线，其中， 表示剩余曲率值0.8 表示剩余曲率值1.0 表示剩余曲率值1.2-表示地表。

图2-2是水平地表剩余曲率扫描曲线，其中， 表示剩余曲率值0.8表示剩余曲率值1.0 表示剩余曲率值1.2-表示地表。

图3-1是利用常规剩余曲率公式获得的计算结果图。

图3-2是利用本发明方法中的优化剩余曲率公式获得的计算结果图，n取3。

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

图5是本发明实施例中初始深度域层速度模型。

图6-1是本发明实施例中初始速度模型起伏地表叠前深度偏移剖面。

图6-2是本发明实施例中初始速度模型起伏地表叠前深度偏移共成像点道集。

图7是本发明实施例中初始模型共成像点道集的剩余曲率谱。

图8是本发明实施例中深度域层速度模型修正量。

图9是本发明实施例中更新后的深度域层速度模型。

图10-1是本发明实施例中更新后速度模型起伏地表叠前深度偏移剖面。

图10-2是本发明实施例中更新后速度模型起伏地表叠前深度偏移共成像点道集。

具体实施方式

下面结合附图对本发明作进一步详细描述：一种适用于起伏地表的深度域层速度修正方法，所述方法首先对偏移距分段，将输入道集按照偏移距进行分组，并进行起伏地表叠前深度偏移，产生共成像点道集；然后进行起伏地表剩余速度分析，计算高阶剩余曲率谱，并根据高阶剩余曲率谱的局部最大值拾取深度-剩余曲率对，然后求取速度模型修正量；再对速度模型修正量进行空间插值和平滑，并将速度模型修正量跟初始速度模型相加得到更新后的起伏地表深度域层速度模型；最后利用更新后的起伏地表深度域层速度模型代替初始速度模型，重新进行起伏地表叠前深度偏移和起伏地表剩余速度分析，从而不断提高速度模型的精度。所述初始速度模型是已有模型，是通过对常规的时间域均方根速度做时深转换得到的。

本发明中所用的公式如下：

(1)起伏地表剩余速度分析公式的推导

起伏地表情况下射线路径如图1所示，在起伏地表情况下，旅行时t_r可以表示为：

t_r＝(SP+PR)/v (1)

其中S是激发点(炮点)，R是接收点(检波点)，P是地下反射点，SP、PR表示射线路径，v是层速度，C点是地表CDP点。

炮点和检波点之间的偏移距L_offset可以表示为：

L_offset＝a+b (2)

其中a与b分别表示炮点半偏移距和检波点半偏移距。

旅行时可以进一步表示为：

$t_r=\frac{1}{v}(\sqrt{a^2+{(z-{\mathrm{dz}}_1)}^2}+\sqrt{b^2+{(z-{\mathrm{dz}}_2)}^2})---\left(3\right)$

其中dz₁表示地表CDP点C与炮点S的高程差，dz₂表示地表CDP点C与检波点R的高程差，z表示地表CDP点C与地下反射点P之间的距离。

在地层水平的情况下，由斯奈尔定律可知，∠SPC＝∠CPR，从而可以得到下式：

$\frac{a}{z-{\mathrm{dz}}_1}=\frac{b}{z-{\mathrm{dz}}_2}---\left(4\right)$

利用(2)式、(4)式求取a、b的表达式并带入(3)式，旅行时可以进一步表示为：

$t_r=\frac{1}{v}(2z-{\mathrm{dz}}_1-{\mathrm{dz}}_2)\sqrt{{\left(\frac{L_{\mathrm{offset}}}{2z-{\mathrm{dz}}_1-{\mathrm{dz}}_2}\right)}^2+1}---\left(5\right)$

我们令 那么(5)式可化简为：

$t_r=\sqrt{{\left(\frac{L_{\mathrm{offset}}}{v}\right)}^2+{\left(\frac{2\hat{Z}}{v}\right)}^2}---\left(6\right)$

当深度域层速度为v_m时，旅行时为：

$t_{\mathrm{rm}}=\sqrt{{\left(\frac{L_{\mathrm{offset}}}{v_m}\right)}^2+{\left(\frac{2{\hat{Z}}_m}{v_m}\right)}^2}---\left(7\right)$

其中 z_m是深度域层速度为v_m时对应的偏移深度。

因为旅行时t_r或t_rm都是地表检波点观测得到的测量时间值，是同一个值，所以在偏移中保持不变，改变的只是成像深度，利用旅行时相等，得到如下关系式：

$\sqrt{{\left(\frac{L_{\mathrm{offset}}}{v_m}\right)}^2+{\left(\frac{2{\hat{Z}}_m}{v_m}\right)}^2}=\sqrt{{\left(\frac{L_{\mathrm{offset}}}{v}\right)}^2+{\left(\frac{2\hat{Z}}{v}\right)}^2}---\left(8\right)$

化简得到：

${\hat{Z}}_m=\sqrt{\frac{L_{\mathrm{offset}}^2}{4}({\mathrm{\γ}}^2-1)+{\hat{Z}}^2{\mathrm{\γ}}^2}---\left(9\right)$

其中，γ＝v_m/v，表示叠前深度偏移所用的深度域层速度与准确层速度之间的比值，也称为剩余曲率值。

将 带入(9)式，从而得到最终的起伏地表剩余速度分析公式：

$z_m=({\mathrm{dz}}_1+{\mathrm{dz}}_2)/2+\sqrt{\frac{L_{\mathrm{offset}}^2}{4}({\mathrm{\γ}}^2-1)+\frac{{(2z-{\mathrm{dz}}_1-{\mathrm{dz}}_2)}^2}{4}{\mathrm{\γ}}^2}---\left(10\right)$

根据(10)式我们就可以对扫描区间内的每一个γ值，利用已知的偏移距信息和高程信息求取不同偏移距对应的偏移深度z_m。一个γ值的所有偏移深度构成剩余曲率扫描曲线，三个γ值的剩余曲率扫描曲线形态如图2-1所示。

当地表水平时，dz₁＝dz₂＝0，(10)式退化为平地表的(11)式即：

$z_m=\sqrt{\frac{L_{\mathrm{offset}}^2}{4}({\mathrm{\γ}}^2-1)+{\mathrm{\γ}}^2{z}^2}---\left(11\right)$

这说明本发明的起伏地表剩余速度分析公式可以兼容平地表剩余速度分析公式。(11)式三个γ值的剩余曲率扫描曲线形态如图2-2所示。通过对比图2-1和图2-2可以看出，两者之间有明显的差别，在复杂地表情况下应该用(10)式进行曲率扫描，(10)式就是起伏地表剩余速度分析公式。

(2)剩余曲率谱的计算

为了获得剩余速度值，我们首先需要计算剩余曲率谱，然后拾取剩余速度谱对应的局部极值，得到‘深度-剩余曲率对’。对于剩余曲率谱的计算，Al-Yahya(1989)提出如下公式计算剩余曲率扫描谱：

$\mathrm{Spe}(z,\mathrm{\γ})=\frac{{\left[\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}U(z_i,i)\right]}^2}{\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}{U}^2(z_i,i)}---\left(12\right)$

其中U(z_i，i)表示当前γ值对应的剩余曲率扫描曲线所经过共成像点道集网格(z_i，i)处的振幅值，i是道集内的道数，代表了偏移距信息，z_i是当前γ值对应的剩余曲率扫描曲线在第i道的偏移深度，由(10)式决定，Ntr表示一个共成像点道集的道数。

式(12)的分辨率比较低，实用中不容易拾取‘深度-剩余曲率对’，为了进一步提高剩余曲率谱的精度和分辨率，本发明对(12)式进行了改进，研究了高阶剩余曲率谱计算公式，具体形式如下：

$\mathrm{Spe}(z,\mathrm{\γ})=\frac{{\left[\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}U(z_i,i)\right]}^{2n}}{\underset{i=1}{\overset{\mathrm{Ntr}}{\mathrm{\Σ}}}{U}^{2n}(z_i,i)}---\left(13\right)$

当n取1时就是(12)式，当n取2，3或其它自然数时，就是高阶剩余曲率谱计算公式，计算的曲率谱收敛程度会更高，分辨率也更高，从而更有利于精确地拾取最佳的‘深度-剩余曲率对’。当n取1和3时剩余曲率谱的计算结果分别如图3-1，图3-2所示，通过对比可以看出当n取3时，剩余速度谱的分辨率得到了明显提高，因此局部极值的自动拾取精度、效率会更高。(13)式就是本文中的高阶剩余曲率谱计算公式，n取3。

通过对图3-2所示的剩余曲率谱按照局部能量最大值进行自动拾取，可以获得该共成像点道集对应的一系列‘深度-剩余曲率对’，即(z_n，γ_n)，n代表拾取的剩余曲率值的个数，z_n代表每个拾取剩余曲率值点的深度。利用(z_n，γ_n)就可以计算当前道的速度模型修正量。

(3)速度模型修正量的求取

根据(z_n，γ_n)中的z_n深度信息和深度间隔可以获得当前的深度样点数n_i，然后利用深度域初始偏移速度计算平均速度，公式如下

${\stackrel{\‾}{v}}_i=n_i/\left(\underset{j=1}{\overset{n_i}{\mathrm{\Σ}}}\frac{1}{v_j}\right)---\left(14\right)$

表示第i个拾取点处的平均速度，n_i代表从地表到z_n处的深度样点数，v_j是初始偏移速度模型第j个点处的速度值。

那么，第一个拾取点处的速度修正量可以表示为：

$\mathrm{delt}{V}_1={\stackrel{\‾}{v}}_1/{\mathrm{\γ}}_1-{\stackrel{\‾}{v}}_1---\left(15\right)$

当i大于1时，第i个拾取点处的速度修正量可以表示为：

${\mathrm{deltV}}_i=\frac{n_i-n_{i-1}}{\frac{n_i}{{\stackrel{\‾}{v}}_i/{\mathrm{\γ}}_i}-\frac{n_{i-1}}{{\stackrel{\‾}{v}}_{i-1}/{\mathrm{\γ}}_{i-1}}}-\frac{n_i-n_{i-1}}{\frac{n_i}{{\stackrel{\‾}{v}}_i}-\frac{n_{i-1}}{{\stackrel{\‾}{v}}_{i-1}}}---\left(16\right)$

对deltV进行普通的空间差值(如反距离加强插值)和平滑处理(如高斯平滑)，就可以得到整个模型空间的速度修正量，将该速度模型修正量与初始速度模型相加就得到更新后的速度模型。(15)和(16)式就是优化的速度模型修正量的求取公式。

如图4所示，本发明方法的步骤包括：

(1)将偏移距分段，并将输入道集按照偏移距进行分组。一般根据偏移孔径将偏移距分为20-35段，每段包含若干输入地震道。各偏移距段一般取等间距排列。偏移距段数和间隔一般在偏移过程中给出；

(2)产生共成像点道集。利用初始速度模型对各段偏移距输入道集进行起伏地表叠前深度偏移，产生共成像点道集。此处一般用起伏地表Kirchhoff叠前深度偏移方法进行偏移，产生的共成像点道集按照偏移距排列，每个共成像点道集的道数与步骤(1)中的偏移距分段数相同(偏移过程是每段偏移距的输入道集产生共成像点道集的一道)；

(3)判断所有共成像点道集是否拉平。此步一般是处理人员浏览道集，不需要计算，如果已经拉平，则转入步骤(7)，如果没有拉平，则转入步骤(4)；

(4)计算剩余曲率谱。首先读入每个需要计算的共成像点道集(所有共成像点或等间隔抽取的共成像点道集)，然后利用起伏地表剩余速度分析公式(10)和高阶剩余曲率谱计算公式(13)计算剩余曲率谱，最后根据局部最大能量值拾取‘深度-剩余曲率对’。

(5)求取速度修正量。对每一个共成像点道集拾取的‘深度-剩余曲率对’利用速度修正量公式(15)、(16)计算深度层速度修正量。对所有已拾取的共成像点进行循环处理；

(6)更新速度模型。对步骤(5)求取的速度修正量进行空间插值和平滑，获得整个模型空间的速度修正量，并将该修正量与初始速度模型(已有模型，对常规的时间域均方根速度做时深转换得到的)相加，得到更新后的速度模型。转入步骤(2)，同时用更新后的速度模型代替步骤(2)中的初始速度模型；

(7)输出偏移剖面。

所述更新后的速度模型可直接用于地震处理中的叠前深度偏移偏移，偏移产生的共成像点道集又可以进一步求取剩余速度修正量并更新速度模型，从而形成一个完整的起伏地表剩余速度分析流程。

下面通过一个实施例来说明本发明方法的效果。

首先给定一个初始的起伏地表速度模型，其构造如图5所示，对其进行起伏地表叠前深度偏移，产生的偏移剖面如图6-1所示，共成像点道集如图6-2所示。可以看出，速度不准确时，偏移剖面质量比较差，反射层的位置不正确(正确位置应该是3000米处)，道集没有拉平，道集的形态受地表影响改变明显，不再是规则的曲线。对所有共成像点道集进行起伏地表剩余速度分析，剩余曲率谱计算并拾取‘深度-剩余曲率对’。计算得到剩余曲率谱如图7所示，可以看出起伏地表剩余曲率谱的局部极值非常明显。速度模型修正量如图8所示，更新后的速度模型如图9所示。利用更新后的速度模型重新进行起伏地表叠前深度偏移处理，产生的深度偏移剖面如图10-1所示，共成像点道集如图10-2所示，可以看出，道集被很好地拉平，反射层位置正确，剖面质量得到很大提高，验证了该方法的正确性。

本发明的方法涉及到起伏地表偏移产生的共成像点道集和炮检点高程，共成像点道集由起伏地表叠前深度偏移产生，共成像点道集的偏移距参数由人工给定；炮检点高程由观测系统给出，在野外测量时就已确定。因此，在地下反射层倾角比较小的情况下，可以直接使用本发明进行起伏地表情况下的剩余速度分析，计算剩余曲率谱以及速度模型修正量，能较好的解决深度域起伏地表速度模型修正问题，并且兼容水平地表情况下的速度模型修正，没有增加内存资源和计算量，计算精度和效率都比较高。

上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。

Claims (3)

1.一种适用于起伏地表的深度域层速度修正方法，其特征在于：所述方法包括：

(1)将偏移距分段，并将输入道集按照偏移距进行分组；

(2)利用初始速度模型对各段偏移距道集进行起伏地表叠前深度偏移，产生共成像点道集；产生的共成像点道集按照偏移距排列，每个共成像点道集的道数与步骤(1)中的偏移距分段数相同；

(3)判断共成像点道集是否拉平，如果已经拉平，则转入步骤(7)，如果没有拉平，则转入步骤(4)；

(4)计算起伏地表剩余曲率谱：首先读入每个地表CDP点处偏移产生的共成像点道集，然后利用高阶剩余曲率谱计算公式计算起伏地表剩余曲率谱，并根据局部最大能量值拾取深度-剩余曲率对；

(5)利用剩余速度求取公式对每一个共成像点道集的深度-剩余曲率对计算深度域层速度修正量；

(6)对步骤(5)获得的所有深度域层速度修正量进行空间插值和平滑处理，获得整个空间的速度模型修正量，并将该速度模型修正量与初始速度模型相加，得到更新后的速度模型；转入步骤(2)，同时用更新后的速度模型代替步骤(2)中的初始速度模型；

(7)输出偏移剖面；

所述高阶剩余曲率谱计算公式为：

$Spe(z_m,\mathrm{\γ})=\frac{{\[\underset{i=1}{\overset{Ntr}{\Σ}}U(z_{mi},i)\]}^{2n}}{\underset{i=1}{\overset{Ntr}{\Σ}}{U}^{2n}(z_{mi},i)};$

所述起伏地表剩余速度分析公式为：

$z_{mi}=({\mathrm{dz}}_{1i}+{\mathrm{dz}}_{2i})\mathrm{/2+}\sqrt{\frac{L_{offseti}^2}{4}({\mathrm{\γ}}^2-1)+\frac{{(2z_i-{\mathrm{dz}}_{1i}-{\mathrm{dz}}_{2i})}^2}{4}{\mathrm{\γ}}^2};$

其中，z_m是深度域层速度为v_m时对应的偏移深度，γ＝v_m/v，表示叠前深度偏移所用的深度域层速度与准确层速度之间的比值，即剩余曲率值，U(z_mi，i)表示当前γ值对应的剩余曲率扫描曲线所经过共成像点道集网格(z_mi，i)处的振幅值，i是共成像点道集内的道数，代表了偏移距信息，z_mi是当前γ值对应的剩余曲率扫描曲线在第i道的偏移深度，由起伏地表剩余速度分析公式求得，Ntr表示一个共成像点道集的总道数，n为大于1的自然数；

dz_1i表示地表CDP点与炮点的高程差，dz_2i表示地表CDP点与检波点的高程差，z_i表示地表CDP点与地下反射点之间的距离，L_offseti是炮点和检波点之间的偏移距。

2.根据权利要求1所述的适用于起伏地表的深度域层速度修正方法，其特征在于：在步骤(1)中根据偏移孔径将偏移距分为20-35段，每段包含至少一道；各偏移距段取等间距排列；偏移距段数和间隔在偏移过程中给出。

3.根据权利要求1所述的适用于起伏地表的深度域层速度修正方法，其特征在于：所述步骤(4)中的地表CDP点位置是人为给定的。