1.一种基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,包括以下步骤:
A、对地震信号叠前道集进行预处理,选取叠前道集的中间道作为基准道,采用DTW计算方法计算待校准道与基准道的最佳偏移量;
B、根据步骤A中待校准道与基准道的最佳偏移量,计算待校准道与基准道的最佳匹配点振幅误差;再采用LCMA计算方法求取最佳匹配点振幅误差中的关键段;
C、将步骤B中关键段的平均偏移量作为该关键段的中心点偏移量,并采用三次样条插值方法求取待校准道中各点的偏移量,根据待校准道中各点的偏移量移动各点的位置,实现叠前道集拉平。
2.如权利要求1所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述步骤A中采用DTW计算方法计算待校准道与基准道的最佳偏移量的目标函数具体表示为:
其中,u为最佳偏移量,N为每一道信号的长度,l为偏移范围向量,e为校准对齐误差矩阵,i∈{1,2...N}。
3.如权利要求2所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述最佳偏移量的目标函数求解方法包括以下分步骤:
S1、根据迭代公式求取最小累加距离d,再按照最短距离目标函数从0到N遍历l得到最短距离DL;所述迭代公式具体表示为:
d[0,l]=e[0,l]
所述最短距离目标函数具体表示为:
S2、根据步骤S1中最短距离对应点的偏移量u[N-1],采用回溯法递归求取最佳偏移量u[i-1];所述回溯法递归公式具体表示为:
u[N-1]=arg min d[N-1,l]
r=N-1,N-2,…,1。
4.如权利要求3所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述步骤B中计算待校准道与基准道的最佳匹配点振幅误差的计算公式具体为:
其中,Ω为最佳匹配点振幅误差,f为基准道,g为待校准道。
5.如权利要求4所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述步骤B中采用LCMA计算方法求取最佳匹配点振幅误差中的关键段的计算公式具体为:
ζk=LCMA(Ω,L)
其中,ζk为求取的第k个关键段,L为关键段长度。
6.如权利要求5所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述步骤C中将关键段的平均偏移量作为该关键段的中心点偏移量具体表示为:
7.如权利要求6所述的基于Segmental DTW的地震信号叠前道集拉平方法,其特征在于,所述步骤C中三次样条插值方法的插值函数具体表示为:
其中,Dnew为原始叠前道集拉平后的数据,Xnew为叠前道集偏移矩阵,X为叠前道集中某一道的自变量,D为叠前道集中某一道的函数值,f为三次样条插值函数。