本发明属于核爆炸地震事件识别，具体涉及一种识别核爆炸地震事件的方法。
背景技术：
：核爆炸地震事件的识别是各国发展禁核试核查国家技术手段的核心需求，现有的地震数据处理系统采取的事件识别技术是事件深度、体波震级与面波震级比率和区域判据高频p/s幅值比判据等三个判据，而且只给出天然地震的事件分类。根据全面禁止核试验条约组织国际数据中心判据筛选的统计结果表明，这些判据仅筛选出了约27％的天然地震事件，未被筛选出的事件占事件总数约三分之二，未筛选出的事件中是否含有核爆炸地震事件是各国所关注的重点，目前监测系统中运行的事件识别判据及识别技术无法满足我国对核爆炸地震识别的需求。技术实现要素：本发明的目的是提供一种识别核爆炸地震事件的方法，解决了核爆炸地震事件监测中的事件识别困难的问题。本发明的技术方案为：一种识别核爆炸地震事件的方法，其特征在于：包括以下步骤：s1、根据特定地区范围及地震监测台网选取合适的地震监测台站作为核爆炸地震事件识别台站；s2、根据选取的台站垂向通道记录到地震事件的初至波时刻t，截取台站垂向通道一定时间段的时间序列数据x(t)；s3、根据天然地震和核爆炸地震在时域、频域上特征的差异以及在时频域的综合特征差异，选取波形复杂度判据、谱比值判据和综合比值判据作为基础判据：s4、利用已知性质的核爆炸地震和天然地震，将波形复杂度判据、谱比值判据和综合比值判据的计算数值分别代入到各自的人工神经网络模型中训练，得到波形复杂度判据、谱比值判据和综合比值判据各自对应的天然地震概率、爆炸地震概率和不确定概率；s5、对各判据计算的天然地震、爆炸地震和不确定的概率进行基于d-s证据理论的证据组合，最终得出各判据联合后的概率，通过比较各类事件概率的大小，确定事件的类型。s1具体包括以下步骤：s1.1、计算地震监测台网中第i个台站到特定地区范围内的最大距离di：di＝max(distance(stai,range))(1)式中stai为第i个台站的经纬度坐标，range为特定地区经纬度范围，distance为台站和特定地区之间距离计算函数，max为取最大值函数；s1.2、选取距离在给定距离d内的台站作为核爆炸地震事件识别台站。s2中，所述时间段为[t-t1，t+t2]；t1、t2为数据的时间段长度，根据地震波列的持续时间设定。s3具体包括以下步骤：s3.1、根据天然地震和核爆炸地震在时域上特征的差异，计算波形复杂度判据数值c：其中x(t)为s2中截取的时间序列数据，l为短时时间窗长度；h为长时时间窗长度；s3.2、根据天然地震和核爆炸地震在频域上特征的差异，计算波形谱比值判据数值sr：其中x(f)为时间序列数据x(t)的振幅谱，l1和h1是低频的起、止频率，l2和h2是高频的起、止频率。s3.3、根据天然地震和核爆炸地震在时频域的综合特征差异，计算综合比值判据数值src：src＝sr/c。s5中，根据d-s证据理论的计算公式，首先得出其中两个判据联合后的天然地震概率、爆炸地震概率和不确定概率，然后将联合后的各概率与第三个判据进行联合，得出最终的天然地震概率、爆炸地震概率和不确定概率。所述d-s证据理论的计算公式如下式所示：其中为冲突因子，反映了证据的冲突程度，m()为基本概率赋值函数；a、b分别为用于联合计算的两个判据，z为两个判据的相交结果，包括天然地震、爆炸地震和不确定三种结果。a∩b＝φ的组合情形包括：a判据地震概率×b判据爆炸概率、b判据地震概率×a判据爆炸概率。天然地震的组合情形包括：a判据地震概率×b判据地震概率、a判据地震概率×b判据不确定概率、a判据地震概率×b判据不确定概率；爆炸地震的组合情形包括：a判据爆炸概率×b判据爆炸概率、a判据爆炸概率×b判据不确定概率、a判据爆炸概率×b判据不确定概率；不确定的组合情形包括：a判据不确定概率×b判据不确定概率。有益效果：本发明通过采用d-s证据理论的证据组合方法，能够正确、有效的给出核爆炸地震事件的识别结果，而且在单个判据失效的条件下，本发明仍然可以给出正确的识别结论。具体实施方式为使本发明的目的、内容和优点更加清楚，下面结合实施例，对本发明的具体实施方式作进一步详细描述。本发明提供一种识别核爆炸地震事件的方法，其特殊之处在于：其包括以下步骤：s1、根据特定地区范围及地震监测台网选取合适的地震监测台站作为核爆炸地震事件识别台站，特定地区为地震事件发生的地区；具体步骤如下：s1.1、计算地震监测台网中第i个台站到特定地区范围内的最大距离di：di＝max(distance(stai,range))(1)式中stai为第i个台站的经纬度坐标，range为特定地区经纬度范围，distance为台站和特定地区之间距离计算函数，max为取最大值函数。s1.2、选取距离在给定距离d内的台站作为核爆炸地震事件识别台站：di<d(2)s2、根据选取的台站垂向通道记录到地震事件的初至波时刻t，截取台站垂向通道时间段为[t-t1，t+t2]的时间序列数据x(t)；t1、t2为数据的时间段长度，根据地震波列的持续时间设定；s3、选取波形复杂度判据、谱比值判据和综合比值判据作为基础判据：s3.1、根据天然地震和核爆炸地震在时域上特征的差异，计算波形复杂度判据数值c：其中x(t)为s2中截取的时间序列数据，l为短时时间窗长度，单位为秒；h为长时时间窗长，单位为秒。s3.2、根据天然地震和核爆炸地震在频域上特征的差异，计算波形谱比值判据数值sr：其中x(f)为时间序列数据x(t)的振幅谱，l1和h1是低频的起、止频率，l2和h2是高频的起、止频率。s3.3、根据天然地震和核爆炸地震在时频域的综合特征差异，计算综合比值判据数值src：src＝sr/c(5)其中，sr为s3.2中计算的谱比值判据数值，c为s3.1中的计算的波波形复杂度判据数值。s4、利用已知性质的核爆炸地震和天然地震对波形复杂度判据、谱比值判据和综合比值判据进行人工神经网络训练，确定各判据使用的人工神经网络模型及参数，各判据所使用的人工神经网络输入为判据计算数值，分别输出各判据的天然地震概率pc_eq、psr_eq、psrc_eq、爆炸地震概率pc_ex、psr_ex、psrc_ex和不确定概率pc_un、psr_un、psrc_un。s5、由d-s证据理论的计算公式得出各判据联合后的概率，包括天然地震概率、爆炸地震概率和不确定概率，如下公式(6)所示：其中为冲突因子，反映了证据的冲突程度，m为基本概率赋值函数；a、b分别为用于联合计算的两个判据，z为两个判据的相交结果；s5.1、将波形复杂度判据与谱比值判据进行证据组合，形成组合1，见表1：表1组合1证据组合方式由表1可知，公式(6)中k值计算如下：k＝pc_eq×psr_ex+pc_ex×psr_eq组合1条件下给出的地震概率pzh1_eq、爆炸概率pzh1_ex和不确定概率pzh1_un计算如下：pzh1_eq＝(pc_eq×psr_eq+pc_eq×psr_un+psr_eq×pc_un)/(1-k)pzh1_ex＝(pc_ex×psr_ex+pc_ex×psr_un+psr_ex×pc_un)/(1-k)pzh1_un＝pc_un×psr_un/(1-k)s5.2、将组合1的结果与综合比值判据的证据组合形成组合2，如表2所示：表2组合1结果与综合比值判据的证据组合由表2可知，公式(6)中k值计算如下：k＝pzh1_eq×psrc_ex+pzh1_ex×psrc_eq组合2条件下给出的地震概率pzh2_eq、爆炸概率pzh2_ex和不确定概率pzh2_un计算如下：pzh2_eq＝(pzh1_eq×psrc_eq+pzh1_eq×psrc_un+psrc_eq×pzh1_un)/(1-k)pzh2_ex＝(pzh1_ex×psrc_ex+pzh1_ex×psrc_un+psrc_ex×pzh1_un)/(1-k)pzh2_un＝pzh2_un×psrc_un/(1-k)由组合2给出的事件识别地震概率pzh2_eq、爆炸概率pzh2_ex和不确定概率pzh2_un即为基于d-s证据理论的最终综合判定结果，通过比较各类事件概率的大小，由此给出事件识别结论。实施例：按照上述方法，对某地区的核爆炸地震和天然地震事件进行了识别测试。先根据事件位置及地震监测台网选取了某地震台站作为核爆炸地震事件识别台站；由台站记录到的两次事件的初至波到时时刻截取了垂向通道的时间序列数据；计算波形复杂度判据、谱比值判据和综合比值判据等判据的计算结果，并分别代入到各自的人工神经网络模型中计算天然地震、爆炸地震和不确定的概率；对各判据计算的天然地震、爆炸地震和不确定的概率进行基于d-s证据理论的证据组合，最终给出了事件的识别结论。表3为判据联合方法对天然地震事件的识别结果，表4为判据联合方法对核爆炸地震事件的识别结果，可以看到，在单个判据失效的条件下，判据联合方法仍然给出了正确的识别结论。表3对某地区天然地震事件的识别结果判据类型地震概率爆炸概率不确定概率结论波形复杂度69.85％8.38％21.77％天然地震谱比值0.00％51.26％48.74％爆炸综合比值75.22％1.77％23.01％天然地震联合结果84.72％10.28％4.99％天然地震表4对某地区核爆炸地震事件的识别结果判据类型地震概率爆炸概率不确定概率结论波形复杂度73.13％7.67％19.20％天然地震谱比值3.36％79.58％17.06％爆炸综合比值6.69％67.24％26.07％爆炸联合结果18.01％79.10％2.89％爆炸以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。当前第1页12