一种针对特定地区地震事件的定位方法技术

本发明专利技术涉及一种针对特定地区地震事件的定位方法，首先选取合适的地震台站组建快速监测地震台网；计算地震台网中的每个台站初至其P波震相的平均速度；将所监测的特定地区区域进行格点划分，根据台站的初至P波震相检测时间和某格点推算事件的发震时间；计算检测到事件的各台站对于推算的事件发震时间的均值和方差；然后选取方差最小的格点作为事件的震中位置，其对应的均值即为事件发震时间，从而对地震事件进行定位。本发明专利技术在特定地区重复地震事件发生的条件下，提出了区域精细格点搜索定位方法，流程简洁，计算量小，便于实现且定位准确。解决了特定地区地震事件快速监测定位的技术难题。

【技术实现步骤摘要】
一种针对特定地区地震事件的定位方法
本专利技术属于地震探测领域，具体涉及一种针对特定地区地震事件的定位方法。
技术介绍
地震监测是全面禁止核试验条约中规定用于监测地下核试验产生的核爆炸地震事件主要技术，现有的地震数据处理系统采取的定位方法是全球格点搜索(GlobalAssociation，简称GA)，通过提取到时最早的三个台站所检测的P震相和S震相的到时差进行画圆交叉确定初始定位区域，并根据最小二乘原理选取初始格点位置，再根据其他台站的到时信息对格点进行迭代，最终得到收敛的定位结果。GA算法原理复杂，对计算资源要求较高，需要高性能的服务器提供计算平台，对于特定地区小范围的事件定位精度较差，不适用于对于特定地区地震事件的快速监测定位。
技术实现思路
本专利技术的目的是提供一种针对特定地区地震事件的定位方法，在特定地区重复地震事件发生的条件下，解决了特定地区地震事件无法快速监测定位的技术难题。本专利技术的技术方案为：一种针对特定地区地震事件的定位方法，其特征在于：包括以下步骤：S1、根据特定地区区域范围内发生的历史地震事件及关联的地震台站，选取合适的地震台站组建快速监测地震台网；S2、根据历史地震事件信息及台站的P波震相检测信息，计算地震台网中的每个台站初至其P波震相的平均速度Vi；S3、将所监测的特定地区区域进行格点划分，得到特定地区区域格点数列；S4、在对特定地区地震事件的实时监测中，在符合震相走时差报警的情况下，对于每个台站i，假设格点数列中的某个格点是产生该信号的地震事件震中位置，根据台站i的初至P波震相检测时间deti和该格点推算事件的发震时间；S5、计算检测到事件的各台站对于推算的事件发震时间的均值和方差；S6、选取方差最小的格点作为事件的震中位置，其对应的均值即为事件发震时间，从而对地震事件进行定位。有益效果：本专利技术在特定地区重复地震事件发生的条件下，提出了区域精细格点搜索定位方法，流程简洁，计算量小，便于实现且定位准确。解决了特定地区重复地震事件无法快速监测定位的技术难题。具体实施方式为使本专利技术的目的、内容和优点更加清楚，下面结合实施例，对本专利技术的具体实施方式作进一步详细描述。本专利技术提供一种针对特定地区地震事件的定位方法，其特征在于：包括以下步骤：S1、选取合适的地震台站组建快速监测地震台网，具体步骤如下：S1.1、计算第i个台站到特定地区范围内的最大距离Di：Di＝max(DISTANCE(stai,range))(1)式中stai为第i个台站的经纬度坐标，range为特定地区经纬度范围，DISTANCE为台站和特定地区之间距离计算函数，max为取最大值函数。特定地区可以根据监测需求自行设定，为重复地震事件发生的地区。S1.2、选取距离在给定距离D内的台站：Di<D(2)S1.3、统计给定时间段范围内特定地区范围内发生的事件总数为N，S1.2中选取的台站中第i个台站检测到的事件数为Ki，其中1<Ki≤N，计算台站的监测灵敏度Ei：S1.4、选取灵敏度在给定阈值E内的台站组成地震台网：Ei>E(4)S2、计算地震台网中的每个台站到特定地区初至P波震相的平均速度Vi，具体步骤如下：式中M为台站i记录到的特定地区历史地震事件的个数，stai为第i个台站的经纬度坐标，Eventj为第j个历史地震事件震中的经纬度坐标，tij为台站i记录到第j个历史地震事件的初至P波震相时间，Tj为第j个历史地震事件的发震时间；S3、将所监测的特定地区区域格点化，计算特定地区区域格点数列，其中某个格点(nx，ny)的位置计算公式如下：式中N1为特定地区经度范围的格点个数，N2为特定地区纬度范围的格点个数，nx为经度格点序列，取值为1～N1；ny为纬度格点序列，取值为1～N2；LONmin和LONmax分别为经度的最大值和最小值，LATmin和LATmax分别为纬度的最大值和最小值。S4、在对特定地区地震事件的实时监测中，对于每个台站i，假设格点数列中的某个格点(nx，ny)是产生地震事件信号的震中位置，根据台站i的初至P波震相检测时间deti和该格点(nx，ny)推算事件的发震时间S5、计算检测到事件的各台站对于推算的事件发震时间的均值和方差S6、选取方差最小的格点作为事件的震中位置，其对应的均值即为事件发震时间。实施例：按照上述方法，对某特定地区的重复地震事件进行了定位测试。先根据特定地区区域范围在历史地震事件数据库中搜索该地震发生的历史地震事件及关联的地震台站，根据上述选台方法组建了对特定地区的地震监测台网；由历史地震事件信息及台站的P波震相检测信息计算各台站到该地区的P波平均速度；对划定的特定地区区域进行格点划分，得到该地区的格点文件；在对特定地区地震事件的实时监测中，在符合震相走时差报警的情况下，根据各台站的P波信号的检测时间、P波平均速度以及区域格点文件对各假设事件格点的发震时间进行了回推，选取了各台站对假设格点事件发震时间方差最小的格点作为了震中位置，并由发震时间均值给出了事件的发震时间，该定位结果与通过其它定位方法的结果一致。以上所述仅是本专利技术的优选实施方式，应当指出，对于本
的普通技术人员来说，在不脱离本专利技术技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本专利技术的保护范围。本文档来自技高网...

【技术保护点】
一种针对特定地区地震事件的定位方法，其特征在于：包括以下步骤：S1、根据特定地区区域范围内发生的历史地震事件及关联的地震台站，选取合适的地震台站组建快速监测地震台网；S2、根据历史地震事件信息及台站的P波震相检测信息，计算地震台网中的每个台站初至其P波震相的平均速度Vi；S3、将所监测的特定地区区域进行格点划分，得到特定地区区域格点数列；S4、在对特定地区地震事件的实时监测中，在符合震相走时差报警的情况下，对于每个台站i，假设格点数列中的某个格点是产生该信号的地震事件震中位置，根据台站i的初至P波震相检测时间deti和该格点推算事件的发震时间；S5、计算检测到事件的各台站对于推算的事件发震时间的均值和方差；S6、选取方差最小的格点作为事件的震中位置，其对应的均值即为事件发震时间，从而对地震事件进行定位。

【技术特征摘要】
1.一种针对特定地区地震事件的定位方法，其特征在于：包括以下步骤：S1、根据特定地区区域范围内发生的历史地震事件及关联的地震台站，选取合适的地震台站组建快速监测地震台网；S2、根据历史地震事件信息及台站的P波震相检测信息，计算地震台网中的每个台站初至其P波震相的平均速度Vi；S3、将所监测的特定地区区域进行格点划分，得到特定地区区域格点数列；S4、在对特定地区地震事件的实时监测中，在符合震相走时差报警的情况下，对于每个台站i，假设格点数列中的某个格点是产生该信号的地震事件震中位置，根据台站i的初至P波震相检测时间deti和该格点推算事件的发震时间；S5、计算检测到事件的各台站对于推算的事件发震时间的均值和方差；S6、选取方差最小的格点作为事件的震中位置，其对应的均值即为事件发震时间，从而对地震事件进行定位。2.根据权利要求1所述的一种针对特定地区地震事件的定位方法，其特征在于：S1具体包括以下步骤：S1.1、计算第i个台站到特定地区范围内的最大距离Di：Di＝max(DISTANCE(stai,range))；式中stai为第i个台站的经纬度坐标，range为特定地区经纬度范围，DISTANCE为台站和特定地区之间距离计算函数，max为取最大值函数；S1.2、选取距离在一定距离D内的台站：Di<D；S1.3、统计给定时间段范围内特定地区范围内发生的事件总数为N，S1.2中选取的台站中第i个台站检测到的事件数为Ki，其中1<Ki≤N，计算台站的监测灵敏度Ei：

【专利技术属性】
技术研发人员：刘哲函，王晓明，唐伟，刘俊民，王燕，李靓，王媛，
申请(专利权)人：禁核试北京国家数据中心，
类型：发明
国别省市：北京,11