引言

日本是第一个建立地震预警系统的国家，日本预警系统在2011年3月11日日本发生9.0级特大地震时发挥了作用，东京地区民众在地震S波最大震动到达之前收到地震报警信息，得以避险(陈会忠等，2011)。随着经济的快速发展和社会的不断进步，社会公众对地震速报提出了越来越高的要求。地震预警具有极其重要的社会效益，它不仅可以减轻人员伤亡和降低次生地震灾害的发生，还可为震后紧急救援和抢修提供第一手的信息(刘希强等，2009)。快速确定地震参数是地震预警系统的关键技术环节之一，关系到预警时间的长短。现场地震预警系统利用P波进行预警可以获得更长的预警时间，即在预警的目标区安装地震计，利用P波比S波传播速度快的原理，用P波到达后数秒的数据来确定地震参数，对S波还没到达的地区通过电磁波发出预警，从而达到预警的目的(周彦文等，2010)。目前国外已经涌现出大量关于地震事件自动实时检测、P波初动拾取、震相类型识别方面的研究方法和成果(Allen，1978，1982;Bear，Kradolfer，1987;Sleeman，Van，1999;Bai，Ken-nett，2000;Diehletal，2009;Kuperkochetal，2010)。国内学者提出了基于分形分维(常旭，刘伊克，2002)、自回归模型(Zhangetal.，2003;王海军等，2003;杨配新等，2004;王继等，2006)、人工神经网络法(Wang，Teng，1995;张范民，李清河，1998;李慧婷，黄文辉，2000;王娟等，2004)、波形变化增长法(朱元清等，2002)、小波变换的主成分分析法(刘希强等，2000)、定尺度比方法(刘希强等，2002)、地震事件自动检测的EFGLP方法(刘希强等，2009)、初至波震相自动精确识别的TOC-AIC方法(刘希强等，2009)、基于幅值和频率的P波识别(毛燕等，2011)等来识别震相，取得了较好的效果。Kuperkoch等(2010)根据高阶统计量(HOS)中的峰度(kurtosis)和偏斜度(skewness)与Akaike信息准则提出了一种地震事件和震相的自动检测和识别算法，并将该算法与Allen(1978)、Bear和Kradolfer(1987)的检测算法对自动获得P波震相读取值进行了比较研究。研究表明，高阶统计量方法显示出了更佳的识别效果，尤其峰度识别效果最佳，且分别使用偏斜度和峰度对3000多个人工测量的P波震相的到时识别其平均偏差分别为(0.38±0.75)s和(0.26±0.64)s。笔者使用峰度方法对山东地震台网记录的单台垂向记录进行了区域地震事件自动检测和震相识别探索研究，同时应用AIC方法对峰度方法所确定的包含地震P波在内的时间窗内记录进行自动分析，探索进一步提高P波震相初至识别的技术和方法。地震波形资料的椭圆型带通滤波方法椭圆型带通滤波器属于无限脉冲响应滤波器，它在有限频率范围内存在传输零点和极点;其通带和阻带都具有等波纹特性，通带、阻带逼近特性良好;对于同样的性能要求，该滤波器比巴特沃斯、切比雪夫这两种滤波器所需用的阶数都低，而且它的过渡带比较窄。图1a给出了阶数为7的椭圆形带通滤波器，其基本参数如下:通带波动系数为0.1，阻带衰减系数为30、数据采样率为100p/s、通带低端截止频率5.5Hz、通带高端截止频率7.5Hz、带通阻带低端截止频率为5Hz、带通阻带高端截止频率为8Hz。从幅频特性可以看出，椭圆型带通滤波器具有通带、阻带逼近特性良好的特点。图1b给出的是由1.5、6.5和10.5Hz正弦波组成的合成信号，其中6.5Hz的信号起始点在2s。图1c给出了根据上述滤波器进行数据滤波的结果，可以看出，滤波后结果与仿真信号相比，在振幅和相位两方面非常一致。在处理区域地震事件(震中距小于120km)的数据时，考虑到其主频分布范围，对椭圆型滤波器的设计参数如下:通带低端截止频率2Hz、通带高端截止频率15Hz、带通阻带低端截止频率为1.5Hz、带通阻带高端截止频率为16Hz。

2区域地震事件触发判断的峰度方法

随机变量X的数学期望表示为式中，p(x)是随机变量X的概率密度。利用期望函数，随机变量X的k阶统计量mk可表示为式中，当k=2时表示X的方差，m2表示X的取值与数学期望的偏离程度。若X比较集中，则m2较小;反之，若取值比较分散，则m2较大。因此，m2是刻画X取值分散程度的一个量，它是衡量X取值分散程度的一个尺度。峰度是指频数分布的集中程度，也就是分布曲线的尖峭程度。随机变量的峰度的表达式为峰度K是四阶中心动差m4除以标准差四次方。将四阶动差m4除以标准差四次方，同样也是为了得到相对数，以便进行不同频数分配的峰度比较。当峰度K＞3时，表示频数分布比正态分布更集中，分布呈尖峰状态，平均数的代表性更大;K＜3时表示频数分布比正态分布更分散，分布呈平坦峰，平均数的代表性较小。

3初至波震相自动精确识别的Kurto-sis-AIC方法

常用的震相自动识别技术是AR-AIC方法，其基本原理是假定地震波记录可分为若干个局部稳态的过程，每个过程满足自回归模型，并认为新的地震波震相到达前后的信号是两个不同的稳态过程(Sleeman，Van，1999)。设长度为N的地震波记录包含有背景信号和地震波信号，通过滑动变量k的不断改变，得到长度分别为［1，k－M］和［k－M+1，N－M－k］数据段的最大似然函数的最大值，进而得到AIC函数:式中，k表示滑动变量(下同)，M表示AR拟合过程的阶数，σ21，max和σ22，max表示最大似然函数的最大值，C2是常数。AR过程的阶数需要通过实际记录并经多次试验计算得到，AR系数根据Yule-Walker方程确定(Haykin，1996)。当所选数据中包含有背景噪声和地震波信号时，AIC函数的最小值对应的时间就是地震震相初至。Maeda(1985)提出了无需计算AR系数，直接根据地震图记录而得到AIC函数的方法，称为VAR-AIC方法。AIC函数的最小值对应的时间就是地震震相初至。AIC函数表示为式中，var(x［1，k］)和var(x［k+1，N］)分别表示两个时间段内数据的方差。VAR-AIC方法属于二阶统计量方法。二阶统计量(如二阶矩、相关函数、功率谱密度函数等)方法需要假设系统是最小相位系统，噪声为高斯白噪声，然而实际信号常常不满足这个假设。高阶统计量与二阶统计量相比具有适于检测和描述系统的非线性，在低信噪比条件下，可有效抑制加性高斯有色噪声的影响，检测和识别较弱的地震信号。用峰度进行精细震相识别的具体思路是用两个时间段数据的峰度作为特征函数代替公式(5)中的方差，本文称为Kurtosis-AIC方法。Kur-tosis-AIC函数表示为动确定的地震事件触发判断结果为参考，选择一个包含有噪声和地震信号的时间窗，Kurtosis-AIC函数的最小值对应时刻即为震相初至。

4资料和结果分析

本文对50次地震的198个单台垂直向记录进行了分析与研究，涉及的台站数量有53个(图2)。所选取单台记录的最大震中距为110.0km、最小为6.0km，震中距小于100km的单台记录有190个，占总数的96%。50次地震中的最小震级为0.6级，最大震级为4.2级。采用逐点滑动的方式对数据的峰度进行实时计算和分析，计算时间窗长为10s，步长为1个采样点。为了总结地震事件触发阈值的分布规律，选择大于噪声段最大峰度值作为每次地震事件触发的最低阈值，峰度超过阈值所对应的时间初定为P波震相初至。然后应用时间段的峰度数据进行Kurtosis-AIC方法分析，其最小值对应时刻即为精细识别的P波震相。图3给出了梁山台记录的2009年6月15日22时27分范县2.4级地震垂直向波形、事件触发时刻及应用峰度和Kurtosis-AIC方法对初至震相到时的拾取结果。当选取峰度阈值等于3时，由峰度方法得到的地震事件触发时刻为49.69s，应用Kurtosis-AIC方法得到的P波震相初至结果为49.68s，与人机交互拾取结果49.64s相比，震相识别的误差为0.04s，识别结果非常理想。图4给出了数据滤波前后信噪比的比较及信噪比随震中距的变化。从图4可以看出，由于对信号进行了椭圆型带通滤波，有效提高了信噪比。图4a中数据滤波前信噪比的变化范围为－3.8～47.6dB，均值为8.7dB，标准差为9.4dB，滤波后信噪比的变化范围为5.1～49.8dB，均值为22.0dB，标准差为9.0dB，可见对观测数据进行滤波后信噪比得到了明显提高。图4b表明对发生在山东地区2级左右小地震进行滤波处理后的信噪比随震中距的变化并未出现明显衰减的特点。图5给出了地震事件触发判断的峰度最低阈值和最大峰度值分布图。从图5a中可以看出，在对数据样本处理时，保证地震事件峰度触发阈值变化范围为2.5～8，均值为3.35，标准差为0.70。与之相对应，图5b中的峰度值的变化范围为4.2～691.9，其中位于4.2～8.0区间有3次，均值为137.0，标准差为118.5。为了保证地震事件既不漏报又不误报，需要在积累的大量样本中选择1个经验阈值。如果本文选定峰度阈值为8，研究样本中有3个单台记录出现漏报，这3次漏报地震的震级都为1.3级，漏报率为1.5%，误报率为0。上述出现漏报资料的信噪比都不高，分别是3.9dB、1.08dB和－1.7dB。如果设置峰度阈值为8，并且地方震的单台垂向记录数据满足信噪比大于4dB，则较大地震不会出现漏报和误报。图6给出了Kurtosis-AIC震相自动识别结果与人机交互识别结果对比的误差分布及误差分布与观测数据信噪比的关系。由图6a可以看出，震相识别误差范围为－0.19～0.5s，平均绝对误差为0.09s，标准差为0.08s。198个单台垂直向记录的处理结果中有127个误差在0.1s内，占总数的64.14%;有174个误差在0.17s内，占88%;有184个误差在0.2s内，占93%。从图6a和图6b对比可以看出，震相识别误差较大的样本相对应的垂向信噪比多数比较低，反应了在人工分析中同样面临的当P波信号较弱时震相初至较难识别这一难题。但总体上Kurtosis-AIC方法对较弱的P波地震信号，也具有较好的分辨能力，达到了比较高的精度。

5结论

本文利用山东地震台网记录的单台垂向记录和峰度方法进行了区域地震事件自动检测和震相识别方法的探索研究，为了减小噪声干扰又不失地方震记录信号，本文设计了保持信号的振幅和相位畸变最小的地震波形资料的椭圆型带通滤波方法。通过研究取得了以下结果:(1)椭圆型滤波器提高了含噪声P波信号的信噪比。(2)提出了一种基于直达P波信号进行区域地震事件实时检测的峰度方法，并用该方法研究得到了山东地震台网记录判断地震事件触发的经验阈值。表明选择合适的峰度经验阈值可同时有效降低地震事件的漏报率和误报率。(3)提出了峰度与AIC方法相结合的精细震相识别新方法对198个单台垂向记录的分析研究表明，与人工识别震相到时的结果相比，根据Kurto-sis-AIC震相自动识别方法得到的震相到时的平均绝对值误差为0.09±0.08s，精度非常高。(4)本文研究样本仍偏少，经验阈值和识别精度有待进一步检验。

特别申明：本站部分文章由网络收集整理.如需使用，请与原作者联系,如有侵权，请联系我们删除,谢谢！