[发明专利]拾取地震波时距曲线的方法及装置

权利要求书

1.一种拾取地震波时距曲线的方法，其特征在于，所述方法包括：

设定包含时距曲线的时窗与时窗内的时间样点数，获取每个时窗内的地震记录；

将获取的地震记录中信噪比最高的地震记录对应的地震道确定为参考地震道；

确定目标地震道，将每个时窗内获取的目标地震道的地震记录与相同时窗内获取的参考地震道的地震记录进行同相叠加，得到叠加结果，所述目标地震道为拾取时距曲线的地震道；

根据所述叠加结果，确定目标地震道的地震波时距曲线相对于参考地震道的地震波时距曲线的时移量；

根据参考地震道的地震波时距曲线以及目标地震道的地震波时距曲线相对于参考地震道的地震波时距曲线的时移量，确定目标地震道的地震波时距曲线。

2.根据权利要求1所述的方法，其特征在于，所述将采集的地震记录中信噪比最高的地震记录对应的地震道确定为参考地震道，包括：

确定每一道地震记录分别与所有道地震记录之间的互相关函数；

构造地震记录互相关函数矩阵，所述互相关函数矩阵中每一列为一道地震记录与地震记录中的每一道地震记录之间的互相关函数；

确定互相关函数矩阵中每一列互相关函数的平均相似系数；

将所述平均相似系数最高的一列确定为参考列，将所述参考列对应的地震道确定为参考地震道。

3.根据权利要求2所述的方法，其特征在于，所述确定每一道地震记录分别与所有道地震记录之间的互相关函数，包括：

利用公式计算每一道地震记录分别与所有道地震记录之间的互相关函数；

其中，k和l为地震道的序号，t为时间样点，Φ_k,l(τ)为第k道与第l道相关的互相关函数，S_k(t)为第k道地震道在时间样点t采集的地震记录，S_l(t)为第l道地震道在时间样点t采集的地震记录，τ为S_l(t)的延迟时间，k＝1,2,…,n，l＝1,2,…,n，n为设定时窗内的地震道数，t＝1,2,…,m，m为时间样点的数量。

4.根据权利要求3所述的方法，其特征在于，确定互相关函数矩阵中每一列互相关函数的平均相似系数，包括：

利用公式计算互相关函数矩阵中每一列互相关函数的平均相似系数；

其中，S_k为第k道的平均相似系数，τ₀为使S_k达到最大值时的延迟时间。

5.根据权利要求3或4所述的方法，其特征在于，在将所述平均相似系数最高的一列确定为参考列之后，所述方法还包括：

按照每列对应的平均相似系数由高到低的顺序，对互相关函数矩阵中的列进行重新排序，得到排序后的互相关函数矩阵，所述排序后的互相关函数矩阵中的第1列为参考列。

6.根据权利要求5所述的方法，其特征在于，所述将每个时窗内采集的目标地震道的地震记录与相同时窗内采集的参考地震道的地震记录进行同相叠加，得到叠加结果，包括：

利用公式计算每列对应的二阶互相关函数；

其中，ΦΦ_1,l(δ)为排序后的互相关函数矩阵中，第l道对应的互相关函数与第1列对应行上的互相关函数的二阶互相关函数，δ为第l道二阶互相关函数的延迟时间，Φ'_1,k(τ)为排序后的互相关函数矩阵中，第k道与第1列对应的地震道相关的互相关函数，Φ'_l,k(τ+δ)为排序后的互相关函数矩阵中，第k道与第l道相关的互相关函数。

7.根据权利要求6所述的方法，其特征在于，根据所述叠加结果，确定目标地震道的时距曲线相对于参考地震道的时距曲线的时移量，包括：

构造函数F(δ)＝ΦΦ_1,l(δ)，确定使函数F(δ)取最大值时δ的值δ₀；

将δ₀确定为第l道地震道的时距曲线相对于参考地震道的时距曲线的时移量。