CN1773310A - 转换波延迟时静校正方法 - Google Patents

Abstract

转换波静校正方法发明是一套完整的转换波静校正方法。它采用纵波激发三分量接收，对Z分量和X分量对应的检波点进行互相关，并把检波点互相关结果叠加得到互相关叠加剖面。根据该剖面的峰值位置确定两分量的初至差，由两分量的初至差计算转换波延迟时。然后根据转换波延迟时建立横波近地表模型，最终求取精确的检波点横波基准面静校正量。获得的静校正量是校正到基准面的结果。

Description

转换波延迟时静校正方法

技术领域

本发明涉及石油地球物理勘探数据处理技术，是一种利用上行透射转换波与纵波在风化层引起的初至时差求取转换波延迟时，建立横波近地表模型，计算横波静校正量的方法。

背景技术

转换波地震勘探是一种特殊的勘探方法。它通过井中爆炸或地面可控震源激发纵波，纵波在地层中向下传播，到达界面后，一部分能量转换为横波(SV波和SH波)。纵波、SV波和SH波三个分量分别向上传播到达地表。地表用三分量检波器进行接收，以获取完整的地震波场信息。

与纵波(P波)不同，横波(S波)不受近地表潜水面起伏的影响。这主要是由于流体饱和的沉积物比干燥沉积物具有更大的胀缩力，剪切模量受水饱和的影响很小。上行S波进入潜水层后，速度变化很小，继续保持低速；而上行P波的低速层终止于潜水面。因此，P波和S波的近地表低速层厚度不一致，使得纵波和横波的静校正量之间几乎没有关系。用纵波静校正量线性近似横波静校正量通常是不正确的。同一位置上的横波静校正量通常比纵波大2～10倍。

现有的《利用上行透射转换波计算横波静校正的方法》发明专利假设上行纵波与横波在低速层中为近垂直路径传播，因此利用转换波初至与纵波初至的时差，加上纵波的静校正量，作为转换波的静校正量，获得了较好的效果。但实际上，上行波与垂直方向存在一定夹角，并且时差与转换波延迟时不同，直接利用相对时差作为相对的静校正量，仍存在一定的静校正量问题。

另外，许多静校正算子同时求解炮点和检波点静校正量。但对于P-SV波而言，只需求解检波点的静校正量。这就要求抑制炮点静校正量求解。因此，必须采用有针对性的静校正方法。

发明内容

本发明利用转换波与纵波的初至时差，提出转换波延迟时的概念。然后根据转换波延迟时建立横波表层模型，从而计算横波的基准面静校正量。

针对背景技术中存在的问题，本发明根据纵波和转换波射线路径的关系，导出转换波延迟时、纵波和转换波初至时差和纵波延迟时的关系。根据这一关系，首先求取转换波延迟时，再建立横波近地表模型，最终求取横波基准面静校正量。具体步骤如下：

1、计算转换波延迟时(图1)

对于两层结构，纵波速度分别为VP1和VP2，横波速度为VS1，A和A′为激发点，D为接收点。

由A点激发，D点接收，纵波折射波路径为ABCD，α为临界角，AB段和CD段的传播速度为VP1，BC段传播速度为VP2；转换波沿界面以VP2速度滑行，到达E点产生转换波，沿射线ED以速度VS1到达D点。

根据折射理论： $\sin \mathrm{\α}=\frac{V_{P1}}{V_{P2}},\sin \mathrm{\β}=\frac{V_{S1}}{V_{P2}}$

由于VP1＞VS1，所以α＞β，转换点E比C点更靠近D点，如图1。

从折射路径可以看出，对于任何一个激发点A或A′，纵波与转换波到接收点D的时间差Δt相同。根据射线路径可以得到：

Δt＝[T_AB+T_BE+T_ED]-[T_AB+T_BC+T_CD]

＝[T_A′B′+T_B′E+T_ED]-[T_A′B′+T_B′C+T_CD]                          (1)

＝T_CE+T_ED-T_CD

转换波与纵波时差Δt可转化为如下的形式：

Δt＝T_CE+T_ED-T_CD＝T_ED-T_EF+T_EF+T_CE-T_CD

＝(T_ED-T_EF)-(T_CD-T_CF)                                              (2)

令：d_P＝T_CD-T_CF，d_S＝T_ED-T_EF

由此得到转换波延迟时的两个部分组成：转换波与纵波的初至时差和纵波延迟时，两者的和即为转换波延迟时：

d_S＝Δt+d_P                                                         (3)

其中，ds定义为转换波延迟时，dp为纵波延迟时，Δt为转换波与纵波的初至时差。

2、建立横波近地表模型，确定厚度h或表层横波速度VS1

根据表层调查可以获得表层横波速度，内插后得到全区的表层横波速度；或者根据纵波速度纵横比参数，间接得到表层横波速度。根据转换波延迟时、表层横波速度和纵波折射波速度，可计算表层的厚度，公式如下：

$h=\frac{V_{S1}{d}_S}{\sqrt{1-{\left(\frac{V_{S1}}{V_{P2}}\right)}^2}}---\left(4\right)$

反之，由于纵波和转换波发生在同一界面，因此可以从纵波数据计算横波近地表模型的厚度，然后根据转换波延迟时等参数计算近地表模型的表层横波速度，公式如下：

$V_{S1}=\frac{h{V}_{P2}}{\sqrt{h^2+V_{P2}^2{d}_s^2}}---\left(5\right)$

3、计算横波静校正量

根据近地表模型，计算横波基准面静校正量。计算公式如下：

$T=-(\underset{i=0}{\overset{n}{\mathrm{\Σ}}}\frac{h_i}{V_i}-\frac{H_d-H_g}{V_s})+\mathrm{\τ}---\left(6\right)$

其中，hi为近地表模型各层的厚度，Vi为各层的横波速度，Hd为基准面高程，Hg为最深层界面高程，Vs为横波替换速度，τ为炮点地震波从井底到地面的垂直传播时间。

本发明的优点在于：

本发明是一种转换波地震勘探的完整的转换波静校正方法。它根据纵波延迟时和初至时差计算转换波延迟时，根据转换波延迟时建立横波近地表模型，最终计算横波静校正量。该方法提供的静校正量是校正到基准面的精确结果。

图2是转换波地震数据静校正前后的叠加剖面对比。可以看出横波静校正问题得到较好解决，同相轴连续性明显提高，扭曲的同相轴变得光滑和连续，尤其对于静校正问题比较突出的左半部，效果显著。

附图说明

图1纵波折射波与转换波折射波射线路径示意图

图2静校正前和静校正后检波点叠加剖面对比图

图3Z分量与X分量互相关叠加剖面

图4横波近地表模型

具体实施方案

实施例1

由于共检波点道集上转换波与纵波初至时差为常数，通过检波点X分量与Z分量(X分量初至是上行透射转换横波的能量，Z分量初至是上行纵波折射的能量)在给定时窗范围内进行互相关，并把检波点互相关结果叠加得到互相关叠加剖面，根据能量峰值确定转换波与纵波的初至时差Δt(图3)。由Δt与对应检波点纵波延迟时dp之和计算该检波点的转换波延迟时ds。

根据横波表层调查已知表层横波速度时，根据(4)式计算低速带厚度，建立横波近地表模型(图4)。然后根据(6)式计算横波的基准面静校正量。

实施例2

由于共检波点道集上转换波与纵波初至时差为常数，通过检波点X分量与Z分量在给定时窗范围内进行互相关，并把检波点互相关结果叠加得到互相关叠加剖面，根据能量峰值确定转换波与纵波的初至时差Δt。由Δt与对应检波点纵波延迟时dp之和计算该检波点的转换波延迟时ds。

由纵波数据计算近地表模型厚度，根据(5)式计算表层横波速度，建立横波近地表模型。然后根据(6)式计算横波的基准面静校正量。

Claims (5)

1.一种转换波延迟时静校正方法，它是地震勘探中利用上行透射转换波计算横波静校正量的方法，其特征在于：根据转换波与纵波的初至时差(Δt)和纵波延迟时(dp)的和计算转换波延迟时(ds)：d_S＝Δt+d_P，然后根据转换波延迟时建立横波表层模型，最后计算横波基准面静校正量。

2.如权利要求1所述的一种转换波延迟时静校正方法，其特征在于：利用共检波点道集上行透射转换横波的初至(X分量)和上行纵波折射的初至(Z分量)，在给定时窗范围内进行互相关，并把检波点互相关结果叠加得到互相关叠加剖面，根据能量峰值确定转换波与纵波的初至时差Δt。

3.如权利要求1或2所述的一种转换波延迟时静校正方法，其特征在于：所述的横波表层模型中，根据表层调查可以获得表层横波速度VS1，内插后得到全区的表层横波速度；或者根据纵波速度纵横比参数，间接得到表层横波速度；然后根据转换波延迟时ds、表层横波速度VS1和纵波折射波速度VP2，计算近地表模型的厚度h，计算公式如下：

$h=\frac{V_{S1}{d}_S}{\sqrt{1-{\left(\frac{V_{S1}}{V_{P2}}\right)}^2}}$

4.如权利要求1或2所述的一种转换波延迟时静校正方法，其特征在于：所述的横波表层模型中，根据纵波数据计算横波近地表模型的厚度h，然后根据转换波延迟时ds和纵波折射波速度VP2，计算近地表模型的表层横波速度VS1，计算公式如下：

$V_{S1}=\frac{{\mathrm{hV}}_{P2}}{\sqrt{h^2+V_{P2}^2{D}_s^2}}$

5.如权利要求1-4所述的一种转换波延迟时静校正方法，其特征在于：所述的横波基准面静校正量计算公式如下：

$T=-(\underset{i=0}{\overset{n}{\mathrm{\Σ}}}\frac{h_i}{V_i}-\frac{H_d-H_g}{V_S})+\mathrm{\τ}$

其中，hi为近地表模型各层的厚度，Vi为各层的横波速度，Hd为基准面高程，Hg为最深层界面高程，Vs为横波替换速度，τ为炮点地震波从井底到地面的垂直传播时间。

Images