引言

我国是世界上岩溶最为发育的国家之一，岩溶面积达363万平方千米，主要分布在黔、桂、滇东、川南、湘西等省区，其中分布最广的是贵州省，约占全省面积的73%，其次是广西壮族自治区，约占全区面积的34% （覃子建等，1996）。近年来随着我国经济的飞速发展，工农业、居民生活用水迅速增加，大量地下水被抽取，导致地下水位不断下降，塌陷事件也逐年增多，岩溶塌陷已成为我国六大地质灾害之一（赵博超等，2015）。然而，天然地震与岩溶塌陷的波形具有较高的相似性，分析人员依靠经验根据波形特征进行二者的识别费时较多，缺乏时效性。因此亟需研究一种能够及时准确地识别天然地震与岩溶塌陷事件的方法，建立准确的事件目录，为快速开展塌陷灾害防治工作提供重要的科技支撑，为地方经济发展、国土资源规划提供理论基础。

近年来国内外研究人员针对地震事件类型的识别开展了广泛而深入的研究，总结了许多时域和频域方面的识别判据。时域的判据主要包括初动方向、震源深度、体波震级与面波震级比、P波初动振幅与P波和S波最大振幅比、P波与S波最大振幅比、P波和S波最大振幅与尾波持续时间比、体波和面波震相特征、P波视速度差异等（吴宝峰，闫峰，2009；Pasyanos，2010；Taylor，2011；刘莎等，2012；王婷婷，边银菊，2015）；而频域的判据主要有能量频谱分布、拐角频率、倒谱曲线、瞬态谱等（张萍等，2005；Yılmaz et al，2013；靳玉贞等，2015；Kortström et al，2016；叶庆东等，2017）。这些方法各有特色，均取得了很好的效果，对于建立准确的地震事件目录具有重要意义；但是这些方法多是单独使用时域特征或频域特征来进行识别，而且地质条件、仪器参数、波形传播路径等差异使得其普适性不足。因此，许多研究人员使用在时域和频域均具有很高分辨率的小波包变换方法来进行特征判据提取。和雪松等（2006）通过对地震和矿震波形进行小波包变换，再计算二者的奇异值，从而达到识别地震的目的；王婷婷等（2014）使用小波包对地震和爆破波形进行分解，根据P波和S波的发育频段不同，使用不同分解系数中的P波和S波能量进行对比，获得了地震和爆破的识别判据；Beccar-Varela等（2016）通过小波变换对美国亚利桑那州M3.0—3.5地震和爆破事件进行分解，利用分解得出的时频能量特征进行地震事件的识别。尽管上述利用小波包分解提取判据的方法能够有效地进行地震识别，但是这些方法主要是针对地震与爆破或矿震，至于地震与岩溶塌陷事件识别的研究则较少。此外，上述方法尽管提取了地震波形在时频域的特征，但均未涉及地震波形分形维数的特征。

为此，本文拟根据天然地震和岩溶塌陷在不同小波包分解频带内存在的差异，结合分形维数计算方法，对比各小波包分解系数的维数差异，以期提出一种能够体现出事件波形时频域特征和分形维数特征并且能够快速有效地识别天然地震与岩溶塌陷的新方法。

1 资料

本文从广西地震台网中心2008—2015年记录到的塌陷和地震事件中，选取记录台站较多、波形较为清晰的广西境内的39个事件，包括17次塌陷事件和22次浅源地震事件，震源深度为5—10 km，震级范围为ML1.0—2.9，所有事件的类型均已由广西地震台网中心确认。从这39个事件中，提取出400条记录清晰的事件波形记录，其中天然地震事件和岩溶塌陷事件的波形各200条，最大、最小震中距分别为237 km和25 km。所选地震、塌陷事件和记录台站的分布如图1所示。

生活中类似的人性之恶与无耻并不罕见。2011年下半年，因为所在大学接受教育部本科合格评估，我对学生的出勤抓得比较严，一有迟到、早退、旷课的行为都登记在册，学生上课讲小话什么的，也会遭到我的批评。没想到期末教学检查时，居然有学生利用匿名向老师书面提意见的机会将我说得一无是处，甚至还捏造了我在课堂上的“违规行为”。好在我所在的文学院主持正义，当年将我推荐为学校的优岗，并公之于众，才将个别学生的这种无耻之举压下去。

本文所选用事件的波形记录的采样率均为100 sps，波形记录的频率范围为0—50 Hz；所有事件只使用垂直向波形进行计算，事件波形的截取宽度均为120 s（即12 000个点），S波的初动位置尽量位于截取窗口的中部。

图1 事件和地震台站位置示意图 Fig. 1 Location of earthquakes （green circles），karst collapses （red circles） and seismic stations （triangles） in this study

2 识别算法设计及基本原理

2.1 识别算法设计

一般来说，一套完整的自动识别算法主要由两个模块组成：其一为特征提取模块，用于从目标数据中提取出识别所需的特征参数；其二为识别判定模块，利用已提取出的特征参数对目标类型进行判定。本文在设计两个模块时考虑到：① 本文所使用的特征提取算法应能体现出天然地震和岩溶塌陷的震源机制在时频域的特点，同时还要能描述出两种事件在不同频带下波形的复杂性；② 由于测震仪器、地质构造、传播路径等差异，每一组事件所提取的特征数据可能存在一定差异，因此识别算法需要具有一定的鲁棒性和容错性，同时从现实应用的角度出发，算法还需具有较快的运算速度。

但该技术也存在着一定的缺点，其能耗较高，在灌溉的过程中必须借助灌溉设备才能达到更好的喷射效果，因此在采用该技术的过程中应建设灌溉机的专用通道，出现大量浪费农田面积的现象[4]。

图2 本文算法流程图 Fig. 2 The flow chart of algorithm in this study

2.2 小波包变换

小波变换是一种重要的时频域波形分析方法，在数字信号处理、图像处理、声纹识别等方面均有广泛应用。而小波包变换是对小波变换的进一步扩展，弥补了小波变换对高频细节特征提取的不足。小波包变换方法是通过多组低通、高通正交滤波器将信号的各频段不断细分，以此来提高信号在整个频带内的分辨能力，进而获得信号中所包含的更多细节。也正因为此，小波包变换也被称作信号时频域的显微镜。

图3 波形分形盒维数示意图 Fig. 3 Schematic diagram of fractal box dimension

2.3 分形维数

分形维数是分形理论中用于描述曲线或者图形不规则程度的重要特征量，但是，在实际使用中往往难以直接计算得到不规则分形体的维数，而只能计算其近似值。因此，本文从快速计算的角度出发，选用计盒维数来计算维数。

考察组一行详细了解了十三师工会网上服务职工的现状和需求，先后前往黄田农场、火箭农场，实地考察了哈密喜折惠干鲜果销售有限公司、哈密果然良品电子商务有限公司的经营状况和电商平台的运行情况，并提出了意见建议。

式中，n为样本输入数据的数量，m为隐含层节点数，gi（x）为神经网络的径向基函数。

鉴于此，本文使用小波包分形方法来提取天然地震和岩溶塌陷的特征，再通过径向基函数(radial basis function，简写为RBF)神经网络进行事件分类的识别判定。小波包分形方法利用了小波包在时频域下高分辨的特点，并通过分形理论对不同频带下波形信号的复杂度进行描述，从而达到兼顾地震信号的细节分析与整体描述合理性的目的（曹茂森等，2005；谢全民等，2011；朱权洁等，2012）。而RBF神经网络除了神经网络普遍具有的自适应性强、能够实现各种非线性映射、有较高的鲁棒性和容错性、运算速度快等特点外，还具备结构简单、收敛快等优势。所以本文使用小波包分形作为特征提取模块，RBF神经网络作为识别判定模块。首先，对事件波形进行归一化处理，然后利用小波包分形提取出事件波形特征；其次，利用已知事件组成训练集完成对RBF神经网络的训练，并将训练好的神经网络用于对未知事件的类型进行识别。图2为本文算法的流程框图。

3 特征提取

RBF神经网络是一种经典的前馈3层神经网络，根据人脑皮层神经网络结构提出，具备以任意精度逼近任意一个连续函数的能力，由输入层、隐含层、输出层组成，其中输入层将输入数据与网络连接起来，隐含层对输入层传递过来的数据进行非线性变换，输出层则负责对隐含层进行加权后输出（Haykin，2001）。图8给出了单输出的RBF神经网络结构图，由图可知整个网络的输出yn可以表示为

图4 地震（a）、岩溶塌陷（b）事件的波形（上）及其时频谱（下） Fig. 4 Waveforms （upper） and spectrograms （lower） of the earthquake （a） and the karst collapse （b）

3.1 事件波形的小波包分解

神经网络的训练和测试由训练集和测试集完成。本文按照发震时间的先后顺序，将前10次地震事件的100条波形数据和前9个岩溶塌陷事件的100条波形数据组成训练集，用于神经网络训练；测试集则是由余下12次地震事件的100条波形数据和8次岩溶塌陷事件的100条波形数据组成，用于验证分类能力。之后分别计算测试集和训练集的地震和岩溶塌陷事件的小波包分形维数，计算结果列于表1。

图5 2010年1月25日ML2.9地震事件的第四层小波包分解 Fig. 5 The fourth-layer wavelet packet decomposition of the ML2.9 earthquake on 25 January 2010

图6 2015年9月9日ML1.8岩溶塌陷事件的第四层小波包分解 Fig. 6 The fourth-layer wavelet packet decomposition of the ML1.8 karst collapse on 9 September 2015

对本文所选的400条事件波形记录进行db4小波包的4层分解，得到16个分解系数波形，对应的频带宽度分别为0—3.125 Hz，3.125—6.250 Hz，…，46.875—50 Hz。图5和图6分别给出了河池市天峨县的地震和塌陷事件的小波包分解结果。可以看出，虽然地震与塌陷事件的优势频率不同，但是在各频段内均有震动波形的记录。因此，为了获取事件在不同频段内的波形特征，提取特征时需对各频段的数据进行计算。

图7 地震与岩溶坍塌事件的小波包分形维数计算结果对比 Fig. 7 Comparison of wavelet packet fractal dimension of earthquake with that of collapse

3.2 小波包系数的分形维数计算

本文采用计盒维数法计算分形维数（张济忠，1995），分别计算各个事件的小波包分解系数，最终获得由16个维数值组成的一维向量，此即为本文所要提取的地震和岩溶塌陷事件的特征。下面将以这16个维数值为输入向量，进行RBF人工神经网络训练和事件类型识别。图7给出了地震和岩溶塌陷事件小波包分解系数的分形维数计算结果，可以看出，岩溶坍塌事件在不同频段内波形的维数高于地震事件波形的维数。

4 RBF 神经网络

由于地震与岩溶塌陷的震源性质和传播路径有所不同，台站记录到两种事件的波形在时频谱中呈现不同的特征。因此，通过时频域分析来进行地震和塌陷事件的特征提取具有可行性。图4为归一化处理后的地震和塌陷事件的时频能量分布图，所示的两事件是由贞丰地震台记录到的地震事件（2010-01-25，ML2.9，震中距Δ=128 km）和塌陷事件（2015-09-09，ML1.8，Δ=124 km），两个事件的震中均位于河池市天峨县。可以看出，地震的优势频率处于0—20 Hz范围内，而塌陷的优势频率处于0—10 Hz范围内。

何良诸十多年前，就领教了赵集的憨劲。那天，何良诸来到北大坎煤矿生活区，临街的小勺酒店，飞出猜拳喝令声。何良诸探头往里瞅，一个女人招呼道：“进来呀。”店内光线暗，酒气、汗腥、臭脚丫子味，呛鼻子，何良诸“啊乞”，打个喷嚏。

计盒维数 DB（F）定义为（张济忠，1995）：对于需要计算维数的目标图形F，构造足够多的边长均为δ的正方形盒子将其完全覆盖（图3），所用到的盒子数量为Nδ（F），那么F的计盒维数DB（F）则为当δ趋于0时Nδ（F）增加的对数速率，可以表示为

本文所使用的RBF神经网络结构为：① 输入层：神经网络的输入为小波包分形得出的16个特征值，所以输入层节点数为16，节点的权值均为1；② 隐含层：径向基函数选择应用广泛的高斯函数，其表达式为

式中，Ci为第i个基函数的中心向量，为标准差，为所选中心之间的最大距离，隐含层节点的数量及其到输出层的权值需经训练得出；③ 输出层：神经网络用于事件类型的判定，故为单输出。

图8 RBF神经网络结构图 Fig. 8 Structure of RBF neural networks

RBF神经网络是否优秀，关键在于隐含层的训练，即基函数中心的选取、节点数量的选择以及节点权值的训练。本文使用matlab神经网络工具箱来进行隐含层的训练，目标均方差设为0.01，高斯函数扩展步长为2.0。训练的具体步骤为：① 建立一个0神经元的网络；② 寻找最大方差样本；③ 增加一个神经元，并使用最大方差的输入样本作为基函数中心；④ 重新设计隐含层权值，使用能够实现最小输出误差的值作为隐含层权值；⑤ 检查输出的均方差是否小于所设定的目标值，若满足该条件则网络训练完毕，否则返回第二步，从剩余样本中继续寻找最大方差样本。

在笔者看来，传记史学应该是传主、传记史学作品、传记作者和作品受众为研究对象，归纳总结传记形态及内容的科学理论与实践，探索解决传记史学作品创作中的实际问题，从而揭示传记史学一般发展规律的一门历史学边缘学科。传记史学具有历史性、阶段性、特殊性、现实性等特征。

5 识别能力检验

使用小波包进行变换前首先需要选取小波包的基函数。常用的基函数有Haar，Daubechies，Mexican Hat，Morlet，Meyer，Symlet，Coiflet，双正交小波等，不同的基函数性质不同。考虑到基函数需具备快速计算的能力，以便增强算法的时效性，同时需要具备精确重构的性质，以便能够完整地反映原始波形的频率特性，为此，本文选取具有紧支撑正交性的Daubechies作为基函数。考虑到需要避免因支撑区间太大而产生边界问题或者因区间太小而使得能量无法集中，选择区间长度为7的4阶Daubechies小波，简写为db4。

表1 地震和岩溶塌陷事件的小波包分形维数均值和方差 Table 1 Mean and variance of fractal dimension of earthquake and collapse events

?

表2 小波包分形维数判据的识别结果 Table 2 The recognition result of the criteria based on wavelet packet fractal

?

从表1可以看出，岩溶塌陷事件各小波包分解系数计算出的分形维数值均明显大于地震事件对应分解系数的维数值。使用训练集中200条波形小波包分形特征向量进行RBF神经网络的训练，训练完成后，分别使用训练集和测试集的波形进行分类检验。训练集和测试集的识别结果如表2所示，可以看出：对于训练集，本文的识别方法可以达到100%的识别率；而对于测试集，地震和岩溶塌陷的识别率可分别达到87%和92%，能够有效地识别出天然地震与岩溶塌陷事件。由表2还可以看出，本文错误识别的共有9个地震事件13条波形记录和5个岩溶塌陷事件的8条波形，其相关信息列于表3。通过分析错误识别的事件波形，本文认为导致错误识别的原因在于：地脉动等背景噪声幅值较大，导致事件波形信号不明显；事件波形存在零漂现象；事件波形中叠加了少量干扰，出现了多个波形信号；事件波形的毛刺多，信号不够平滑；事件波形中有断记，部分波形为直线。

表3 错误识别事件的相关信息 Table 3 The directory of false recognition

?

天然地震事件和岩溶塌陷事件的波形差异的本质是两种事件产生的震源机制不同。天然地震从本质上来讲地壳介质在构造应力作用下的破裂过程，通常震源深度较深，破裂过程较长，同时地震波经地下介质多次反射、折射、散射、叠加，导致地震仪所记录到的地震波频率范围更宽，尾波持续时间更久；而岩溶塌陷则发生在地表或者地下几十米处，是由于塌陷体从周边岩体上脱落并撞击底部岩体而产生的瞬时压缩源，能量释放时间短，波的传播路径较简单，未经过多层介质干扰，频带范围窄，尾波持续时间短（靳玉贞等，2015）。

三是交通运输主管部门应搭建统一的信息平台，实现航道主管部门建设的电子航道图与海事AIS数据，以及船闸运行单位、过闸船舶之间的信息互通共享。

为了能够得到量化的地震波形小波包分解的特征，本文使用分形维数来计算特征参数。地震波形分形维数的大小能够在一定程度上反映出岩石所受加载应力的差异。当岩石受到不同大小应力作用时，所产生的瞬态弹性波（称为声发射技术）的分形维数变化规律也验证了这一点。尹贤刚等（2009）运用分形理论方法研究了不同加载应力比阶段岩石声发射的分形维数，结果表明加载应力达到岩石峰值应力的40%左右时，分形维数即可达到最大值，之后逐步减小直至发生破裂；张黎明等（2015）对大理岩卸围压破坏过程的研究结果表明，应力达到峰值前，声发射分形维数会大幅度降低。由此可见，由于岩溶塌陷事件震动所产生的加载力并未达到岩石破裂峰值，而天然地震事件的加载力已经达到破裂峰值，所以岩溶塌陷事件波形的分形维数要大于地震事件的分形维数，这可以作为识别两种事件的判据。

6 讨论与结论

2） 时频谱图和小波包分解系数对比显示，天然地震波形的优势频段为0—20 Hz，岩溶塌陷的优势频段为0—10 Hz。尽管两种事件的优势频率不同，但是二者所产生的波形在各个频段上均有记录，因此进行时频分析时不能只计算优势频率，应将全频段考虑在内。

1） 小波包分形维数识别方法对于训练集的识别能力能够达到100%，对于测试集的平均识别能力达89.5%，说明该方法能够有效地识别天然地震与岩溶塌陷事件。而且，一个地震事件通常会被多个台站记录到，综合多个台站的识别结果，可提高事件识别的准确率。

本文利用基于小波包的分形维数和RBF神经网络识别技术，对2008—2015年广西境内的17次岩溶塌陷事件和22次浅源地震事件进行了识别研究和识别效果的检验，得到以下几点认识：

网络信息平台也是地方政府提供政务服务的有效载体。依托网络开展 “互联网+”政务服务，可以整合各类政务服务事项和业务办理等信息，构建全流程一体化服务平台，通过网上大厅、办事窗口、移动客户端、自助终端等多种形式为自然人和法人等提供一站式办理的政务服务，实现 “优化服务流程，创新服务方式，推进数据共享，打通信息孤岛，推行公开透明服务，降低制度性交易成本，让企业和群众少跑腿、好办事、不添堵”[1]。2016年来，国务院陆续出台专项政策，布局 “互联网+”政务服务，各地积极响应，快速建设，取得了良好效果。

3） 地震事件的各小波包分形维数均小于塌陷事件的小波包分形维数，这可能与岩石受到的应力大小有关，岩石受力越接近临界状态，所产生震动的分形维数越小。

在实际工作中，可将本文获得的量化指标作为快速自动化识别地震与岩溶塌陷事件的判定依据之一。鉴于地震破裂过程、传播路径等因素的复杂性，地震事件性质的判定还需结合近台记录P波初动符号、P波与S波振幅比、现场地震地质调查、地震宏观现象等资料综合判定。

不同地区的地质背景、塌陷方式、传播路径等方面差异较大，波形记录的特征也有所不同。因此，本研究获得的指标仅具有一定参考价值，不一定适用于其它地区。本文研究样本为已知类型事件，下一步将收集更多的事件样本，进行更广泛的识别训练和研究，不断地完善和验证判定指标的普适性。

刘名卓[9]认为网络课程中既包括教学内容和学习资源，也包括教学活动。“集教师角色、教材角色和教学媒体为一体。”这是从内容上对网络课程进行定义的。

参 考 文 献

曹茂森，任青文，王怀洪. 2005. 基于小波与分形理论的地震异常检测[J]. 地球物理学报，48（3）：672–679.

Cao M S，Ren Q W，Wang H H. 2005. A method of detecting seismic singularities using combined wavelet with fractal[J].Chinese Journal of Geophysics，48（3）：672–679 （in Chinese）.

和雪松，李世愚，沈萍，冯全雄. 2006. 用小波包识别地震和矿震[J]. 中国地震，22（4）：425–434.

He X S，Li S Y，Shen P，Feng Q X. 2006. A wavelet packet approach to wave classification of earthquakes and mining shocks[J]. Earthquake Research in China，22（4）：425–434 （in Chinese）.

靳玉贞，林木金，范晓瑜，刘晓萍，何佳，杨世英，孟彩菊. 2015. 山西地区爆破、塌陷（矿震）特殊地震动特征识别[J]. 地震地磁观测与研究，36（3）：63–66.

Jin Y Z，Lin M J，Fan X Y，Liu X P，He J，Yang S Y，Meng C J. 2015. Feature recognition of explosion，collapse （mine earthquake） to specific ground motion in Shanxi region[J]. Seismological and Geomagnetic Observation and Research，36（3）：63–66 （in Chinese）.

刘莎，杨建思，田宝峰，郑钰，姜旭东，徐志强. 2012. 首都圈地区爆破、矿塌和天然地震的识别研究[J]. 地震学报，34（2）：202–214.

Liu S，Yang J S，Tian B F，Zheng Y，Jiang X D，Xu Z Q. 2012. Discrimination between explosions，mine collapses and earthquakes in capital region of China[J]. Acta Seismologica Sinica，34（2）：202–214 （in Chinese）.

覃子建，罗远模，王继安. 1996. 喀斯特地震研究[J]. 贵州师范大学学报：自然科学版，14（4）：6–16.

Qin Z J，Luo Y M，Wang J A. 1996. The research on earthquakes in karst areas[J]. Journal of Guizhou Normal University：Natural Science，14（4）：6–16 （in Chinese）.

王婷婷，边银菊，张博. 2014. 地震与爆破的小波包识别判据研究[J]. 地震学报，36（2）：220–232.

Wang T T，Bian Y J，Zhang B. 2014. Recognition criteria of earthquakes and explosions based on wavelet packet analysis[J].Acta Seismologica Sinica，36（2）：220–232 （in Chinese）.

王婷婷，边银菊. 2015. 振幅衰减特性在地震与爆破识别中的应用[J]. 地震学报，37（1）：169–179.

Wang T T，Bian Y J. 2015. Amplitude attenuation and its application to earthquake and explosion discrimination[J]. Acta Seismologica Sinica，37（1）：169–179 （in Chinese）.

吴宝峰，闫峰. 2009. P波视速度方法识别爆破群的初步研究[J]. 地震地磁观与研究，30（2）：74–76.

Wu B F，Yan F. 2009. Preliminary study on identification of blasts cluster using the P-wave apparent velocity[J]. Seismological and Geomagnetic Observation and Research，30（2）：74–76 （in Chinese）.

谢全民，龙源，钟明寿，刘好全，周翔. 2011. 小波包与分形组合技术在爆破振动信号分析中的应用研究[J]. 振动与冲击，30（1）：11–15.

Xie Q M，Long Y，Zhong M S，Liu H Q，Zhou X. 2011. Application of wavelet packet and fractal combination technology in blasting vibration signal analysis[J]. Journal of Vibration and Shock，30（1）：11–15 （in Chinese）.

叶庆东，丁志峰，王生文，余大新，郑晨. 2017. 汶川地震科学钻探3号井孔附近微震震源参数确定及其意义[J]. 地球物理学报，60（7）：2716–2732.

Ye Q D，Ding Z F，Wang S W，Yu D X，Zheng C. 2017. Determining the source parameters of the microearthquakes near the third borehole of the Wenchuan Earthquake Fault Scientific Drilling （WFSD-3） and its implications[J]. Chinese Journal of Geophysics，60（7）：2716–2732 （in Chinese）.

尹贤刚，李庶林，唐海燕，裴建良. 2009. 岩石破坏声发射平静期及其分形特征研究[J]. 岩石力学与工程学报，28（增刊2）：3383–3390.

Yin X G，Li S L，Tang H Y，Pei J L. 2009. Study on quiet period and its fractal characteristics of rock failure acoustic emission[J]. Chinese Journal of Rock Mechanics and Engineering，28（S2）：3383–3390 （in Chinese）.

张济忠. 1995. 分形[M]. 北京：清华大学出版社：123.

Zhang J Z. 1995. Fractal[M]. Beijing：Tsinghua University Press：123 （in Chinese）.

张黎明，任明远，马绍琼，王在泉，王金良. 2015. 大理岩卸围压破坏全过程的声发射及分形特征[J]. 岩石力学与工程学报，34（增刊1）：2862–2867.

Zhang L M，Ren M Y，Ma S Q，Wang Z Q，Wang J L. 2015. Acoustic emission and fractal characteristics of marble during unloading failure process[J]. Chinese Journal of Rock Mechanics and Engineering，34（S1）：2862–2867 （in Chinese）.

张萍，蒋秀琴，苗春兰，邹向荣. 2005. 爆破、矿震与地震的波谱差异[J]. 地震地磁观测与研究，26（3）：24–34.

Zhang P，Jiang X Q，Miao C L，Zou X R. 2005. The differences of wave spectrum among explosion，mine tremor and earthquake[J]. Seismological and Geomagnetic Observation and Research，26（3）：24–34 （in Chinese）.

赵博超，朱蓓，王弘元，赖柄霖. 2015. 浅谈岩溶塌陷的影响因素与模型研究[J]. 中国岩溶，34（5）：515–521.

Zhao B C，Zhu B，Wang H Y，Lai B L. 2015. Influence factors and mathematical models of karst collapses[J]. Carsologica Sinica，34（5）：515–521 （in Chinese）.

朱权洁，姜福兴，尹永明，于正兴，温经林. 2012. 基于小波分形特征与模式识别的矿山微震波形识别研究[J]. 岩土工程学报，34（11）：2036–2042.

Zhu Q J，Jiang F X，Yin Y M，Yu Z X，Wen J L. 2012. Classification of mine microseismic events based on wavelet-fractal method and pattern recognition[J]. Chinese Journal of Geotechnical Engineering，34（11）：2036–2042 （in Chinese）.

Beccar-Varela M P，Gonzalez-Huizar H，Mariani M C，Tweneboah O K. 2016. Use of wavelets techniques to discriminate between explosions and natural earthquakes[J]. Phys A Stat Mech Appl，457：42–51.

Haykin S. 2001. Neural Networks：A Comprehensive Foundation[M]. Beijing：Tsinghua University Press：278.

Kortström J，Uski M，Tiira T. 2016. Automatic classification of seismic events within a regional seismograph network[J]. Comput Geosci，87：22–30.

Pasyanos M E. 2010. A general method to estimate earthquake moment and magnitude using regional phase amplitudes[J]. Bull Seismol Soc Am，100（4）：1724–1732.

Taylor S R. 2011. Statistical discriminants from two-dimensional grids of regional P/S spectral ratios[J]. Bull Seismol Soc Am，101（4）：1584–1589.

Yılmaz Ş，Bayrak Y，Çınar H. 2013. Discrimination of earthquakes and quarry blasts in the eastern Black Sea region of Turkey[J]. J Seismol，17（2）：721–734.