[发明专利]一种三维VSP数据初至波拾取方法

权利要求书

1.一种三维VSP数据初至波拾取方法，其特征在于：具体步骤包括：

1)用地震勘探方法取得三维三分量VSP数据集，进行观测系统定义、去除单频干扰、振幅恢复、静校正处理；

2)抽取若干炮不同井源距的VSP数据，初步拾取其初至，形成每个观测点深度的初至时间与井源距关系表；

3)对每一道数据，利用初至井源距关系表采用三次样条插值出此道的先导初至t_p；

4)采用一个大于零的时窗Δt，截取(t_p-Δt)～(t_p+Δt)时间段的数据；

5)对截取的数据进行三分量合成模数道处理，公式为：

x_i、y_i、z_i表示三分量VSP数据第i个样点振幅值，m_i表示模数道第i个样点振幅值；

6)根据精度要求确定新的采样率，采用三次样条算法按新采样率把模数道数据重采样；

7)采用以下公式计算模数道数据每个样点的相邻振幅差的乘积值dm_i，dm_i＝(m_i-m_i-1)·(m_i+1-m_i)，m_i表示模数道第i个样点振幅值，dm_i小于零的值为波峰，并找出其中的最大波峰值A；

8)用小于1的比例因子α，计算每一个波峰与A的比值，确定比值大于α的波峰，其中的第一个波峰是初至后的第一个波峰，这个波峰值计为B；

9)用小于1的比例因子β，计算B之前的数据与B的比值，确定比值小于β的距离B最近的一个样点，对应的时间为拾取初至；

10)重复(3)～(9)至全部三维VSP数据拾取；

11)在炮集内利用相邻深度初至时差辅助手动拾取初至；

12)对全部数据按井源距分组，在每个井源距内计算相邻深度初至时差dt_j，公式为：dt_j＝(t_j-t_j-1)/(d_j-d_j-1)，t_j表示第j道初至，d_j表示第j道深度值；对每一个深度的时差值序列按大小排序，找到中间值meddt，用允许时差范围ΔT₂判断，找出dt＞meddt+ΔT₂或dt＜meddt-ΔT₂的初至，返回到步骤11)重新拾取初至；

13)对全部数据按炮检水平方位进行分组，在每个炮检水平方位组内把初至按井源距从小到大排序，对每一个深度初至计算后一炮初至与前一炮的初至差dt_ij，公式为：dt_ij＝(t_ij-t_(i-1)j)，t_ij表示第j深度道第i炮初至。找出dt_ij＜0的初至，返回到步骤11重新拾取初至；

14)对全部数据按接收点深度进行分组，在每个接收点深度内有地理坐标y、x及初至t三列数据，利用基于三角形的线性插值算法对(y，x，t)进行二元函数曲面拟合，形成规则网格数据(Y，X，T)，对(Y，X，T)进行多道二维中值滤波处理得到(Y，X，Tm)，计算dT＝|T-Tm|，用允许时差范围ΔT₃判断，找出dT＞ΔT₃的初至，返回到步骤11重新拾取初至。

2.根据权利要求1所述的方法，其特征在于：步骤2)所述的井源距VSP数据包括最小井源距和最大井源距数据。

3.根据权利要求1所述的方法，其特征在于：步骤6)新的采样率为0.01-0.1ms。