[发明专利]一种检测地震信号的控制方法及装置

权利要求书

1.一种检测地震信号的控制方法，其特征在于，包括：

获取地震波的N个样点的离散序列{x_i}，当前时窗长度为M；

计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和；将时窗向前移动一个样点，计算新时窗内样点的各值，是在上一个时窗的值基础上，减去上一个时窗第1个样点的值，并加上当前时窗最后一个样点的值，依次类推，得到所有时窗N个样点的峭度和偏斜度的数值；并通过得到的所有时窗N个样点的峭度和偏斜度的数值得到峭度曲线和偏斜度曲线；

通过得到的地震波的峭度曲线和偏斜度曲线，判断地震波到时的位置。

2.如权利要求1所述的方法，其特征在于，

所述得到所有时窗N个样点的峭度K的数值的步骤，包括：

计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和，其中，将当前时窗长度为M的时窗内离散序列的算术和、平方和、立方和、4次方和分别记为sm_j,j=1,2,3,4，公式如下所示：

sm1=Σi=1Mxi,]]>

sm2=Σi=1Mxi2,]]>

sm3=Σi=1Mxi3,]]>

sm4=Σi=1Mxi4;]]>

然后通过公式

K=1(M-1)·σx(sm4-4x‾·sm3+6x‾2·sm2-3Mx‾4)]]>得到当前时窗长度为M的峭度K，

其中为：

x‾=1Msm1,]]>

标准差σx=sqrt(1MΣi=1M(xi-x‾)2)=sqrt(1M(sm2-M·x‾2));]]>

将时窗向前移动一个样点，计算新时窗内的各值，是在上一个时窗的值基础上，减去上一个时窗第1个样点的值f，并加上当前时窗最后一个样点的值r，公式如下所示：

sm1n=sm1+(r-f),]]>

sm2n=sm2+(r2-f2),]]>

sm3n=sm3+(r3-f3),]]>

sm4n=sm4+(r4-f4),]]>

得到新时窗的后，再计算新时窗的峭度K的数值，依次类推，

直到得到所有时窗的峭度K的数值。