[发明专利]一种改进偏斜度的自动拾取S波震相方法

摘要

本发明提出一种改进偏斜度的自动拾取S波震相方法，流程包括：分别输入地震东西向、南北向地震记录；计算特征函数，即地震动水平向合振幅；计算特征函数的改进偏斜度；计算改进的偏斜度相对变化率绝对值序列；拾取改进的偏斜度相对变化率绝对值序列的最大值，该最大值对应的时刻就是S波初至时刻。本发明采用改进偏斜度，其主要部分与标准差的高次方成反比，突出了标准差对描述涨落的作用；改进方法标准差高次方的使用，使得改进偏斜度计算区间进入地震信号区时单调变化，也使得地震信号与非地震信号的改进偏斜度的差异明显，进而通过改进偏斜度相对变化率确定震相到时，使得拾取震相到时的误差小；本发明不进行任何滤波处理，不会使波形失真。

主权项

1.一种改进偏斜度的自动拾取S波震相方法，其特征在于，具体流程如下：步骤1：分别输入地震东西向、南北向地震记录；步骤2：根据东西、南北向地震记录，用公式(1)计算特征函数xi(t)，即地震动水平向合振幅，其中，xEW(t)表示东西向地震记录，xNS(t)表示南北向地震记录：步骤3：在整个地震记录时段内，用公式(2)、(3)和(4)计算特征函数改进的偏斜度K(t)：其中，M是记录点数，即记录时长；为在采样时段M内记录数据平均值，σ(t)为采样时段M内的记录数据x_i(t)标准差，n为正整数；步骤4：在整个地震记录时段内，用公式(5)计算改进偏斜度相对变化率绝对值，得到序列{ri(t)}：其中，Ki‑1(t‑1)为前一时刻改进偏斜度，Ki(t)为当前时刻改进偏斜度；步骤5：用公式(6)拾取改进偏斜度相对变化率绝对值序列{ri(t)}的最大值rmax，其对应的时刻就是S波初至时刻：rmax＝max{ri(t)}。           (6)