[發明專利]一種檢測地震信號的控制方法及裝置

有效

• 申請號: 201210204753.5
• 申請日: 2012-06-15
• 公開（公告）號: CN102707312A
• 公開（公告）日: 2012-10-03
• 發明（設計）人: 劉勁松
• 申請（專利權）人: 中國科學院地質與地球物理研究所
• 主分類號: G01V1/22
• 分類號: G01V1/22;G01V1/28
• 代理公司: 北京國昊天誠知識產權代理有限公司 11315
• 代理人: 許志勇

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的數值。