CN105572738A - 一种采用单个台站检测特定地区核爆炸地震事件的方法 - Google Patents
一种采用单个台站检测特定地区核爆炸地震事件的方法 Download PDFInfo
- Publication number
- CN105572738A CN105572738A CN201610064683.6A CN201610064683A CN105572738A CN 105572738 A CN105572738 A CN 105572738A CN 201610064683 A CN201610064683 A CN 201610064683A CN 105572738 A CN105572738 A CN 105572738A
- Authority
- CN
- China
- Prior art keywords
- station
- data
- time
- nuclear explosion
- particular locality
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6167—Nuclear
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/65—Source localisation, e.g. faults, hypocenters or reservoirs
Abstract
本发明公开了一种采用单个台站检测特定地区核爆炸地震事件的方法,基于特定地区历史核爆炸地震事件在台站监测记录的信号特征,先截取了历史核爆炸地震事件在选定台站上的信号记录,并进行了时频域频谱分析确定了滤波频带,在选定的台站的实时数据检测中发现检测信号后,截取了实时波形数据与历史数据进行信号相关计算,判定该特定地区发生了核爆炸地震。本发明解决了特定地区核爆炸地震事件不能使用单台进行检测的问题,能够对特定地区的天然地震事件进行过滤,实现针对特定地区核爆炸地震事件的单台监测报警,完善了对于特定地区核爆炸事件的监测能力。
Description
技术领域
本发明属于地震事件检测领域,具体涉及一种采用单个台站检测特定地区核爆炸地震事件的方法。
背景技术
地震监测是全面禁止核试验条约中规定用于监测地下核试验产生的核爆炸地震事件主要技术,全面禁止核试验条约组织国际数据中心和的地震数据处理系统采取的报警方法是基于台网监测的报警方法,在汇集了台网台站数据后进行台站信号特征提取,使用全球格点搜索方法进行事件关联形成自动处理事件列表,然后对事件列表中的事件信息进行判断,符合报警条件时触发报警。目前全球核试验地震监测台网较为稀疏,核试验也在向小型化的趋势发展,对于较小震级的核爆炸地震事件很有可能只有个别的台站能够记录到,满足不了系统形成检测事件所要求的三个或三个以上监测台站记录到的条件,导致无法对此类核爆炸事件进行报警。
发明内容
本发明的目的是提供一种采用单个台站检测特定地区核爆炸地震事件的方法,解决了现有地震检测技术不能使用单台进行报警的技术难点,进一步完善了对于特定地区核爆炸事件的检测能力。
本发明的技术方案为:一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于,包括以下步骤:
S1、选取一个地震监测台站作为信号相关报警台站;
S2、读取特定地区历史核爆炸地震事件在选定台站上的监测数据,使用公式(5)进行时频域频谱分析得到信号时频域谱图:
S=spectrogram(x,window,noverlap,nfft,fs)(5)
式中x为监测台站记录的数据,window为滑动窗口内的数据点数,noverlap为计算过程中滑动窗口重叠的数据点数,nfft为快速傅里叶变换的数据长度,fs为监测台站的采样率,spectrogram为时频域频谱分析函数;
观察信号时频域谱图确认信号时间,截取该时间前后指定时间长度L的数据作为信号相关处理的模板数据H;根据信号时频域谱图确定能量集中分布频带F;
S3、获取监测台站实时监测数据,计算信号检测时刻短时绝对平均值STAj和长时绝对平均值LTAj:
式中Ak表示是采样点序号k的数值,j表示信号检测时刻采样点序号,n表示短时绝对平均值时间窗口数据点数,m表示长时绝对平均值时间窗口数据点数;
S4、比较台站数据的STAj/LTAj结果与设定的信号检测阈值Th大小,当 时,
根据比较结果和台站数据信息给出台站的信号检测点时间Tdet。
S5、截取该台站信号检测点时间Tdet前后指定时间长度为L的数据W,结合步骤S2得到的模板数据H及能量集中分布频带F,同时在该集中分布频带F进行滤波后按照公式(10)逐点错动计算两列数据W和H的互相关系数ρxy:
式中x(l)和y(l)分别是数据W和模板数据H的时间序列,n表示数据序列的个数,和分别是数据W和模板数据H的n个数据点的均值;
S6、选取逐点错动计算得到的互相关系数ρxy序列中的最大值作为两段数据的最终相关系数结果,当最终相关系数大于等于设定的相关系数阈值时,即判定该特定地区发生核爆炸地震;
S7、当最终相关系数小于设定的相关系数阈值时,重复进行步骤S3-S6,实现对特定地区单个台站的核爆炸地震信号进行检测。
有益效果:本发明可以对特定地区的天然地震事件进行了过滤,实现了针对特定地区核爆炸地震事件的单台检测并报警,完善了对于特定地区核爆炸事件的监测能力。
附图说明
图1为本发明对某特定地区天然地震信号相关检测结果;
图2为本发明对某特定地区核爆炸地震信号相关检测结果。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面对本发明的具体实施方式作进一步详细描述。
一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于,包括以下步骤:
S1、选取一个合适的地震监测台站作为信号相关报警台站,具体步骤如下:
S1.1、按照公式(1)计算第i个台站到特定地区范围内的最大距离Di:
Di=max(distance(stai,range))(1)
式中stai为第i个台站的经纬度坐标,range为特定地区经纬度范围,distance为台站和特定地区之间距离计算函数,max为取最大值函数;
S1.2、选取距离在给定距离D内的所有台站,即每个台站的距离符合以下要求:
Di<D(2)
S1.3、统计给定时间段范围内特定地区范围内发生的事件总数为N,1.2)步骤中选取的台站中第i个台站检测到的事件数为Ki,其中1<Ki≤N,按照公式(3)计算第i个台站的监测灵敏度Ei:
S1.4、选取灵敏度在给定阈值E内的一个台站,符合以下要求,作为信号相关报警台站:Ei>E(4)
S2、读取特定地区历史核爆炸地震事件在选定台站上的监测数据,使用公式(5)进行时频域频谱分析得到信号时频域谱图,确认信号时间和能量集中分布频带F:S=spectrogram(x,window,noverlap,nfft,fs)(5)
式中x为监测台站记录的数据,window为滑动窗口内的数据点数,noverlap为计算过程中滑动窗口重叠的数据点数,nfft为快速傅里叶变换的数据长度,fs为监测台站的采样率,spectrogram为时频域频谱分析函数;
观察信号时频域谱图确认信号时间,截取该时间前后指定时间长度L(该长度可以为10s-20s)的数据作为信号相关处理的模板数据H;根据信号时频域谱图确定能量集中分布频带F;
S3、使用连续数据传输协议获取监测台站实时监测数据,通过解析帧数据类型、帧数据起始时间、帧数据结束时间获取台站监测数据信息。
S4、根据实时获取的台站监测数据,计算信号检测时刻短时绝对平均值STAj和长时绝对平均值LTAj:
式中Ak表示是采样点序号k的数值,j表示信号检测时刻采样点序号,n表示短时绝对平均值时间窗口数据点数,m表示长时绝对平均值时间窗口数据点数;n、m为根据采样率来设定的固定值;
S5、比较台站数据的STA/LTA结果与设定的信号检测阈值Th大小,当
根据比较结果和台站数据信息给出台站的信号检测点时间Tdet。
S6、截取该台站信号检测点时间Tdet前后指定时间长度为L的数据W,结合步骤S2得到的模板数据H及能量集中分布频带F,同时在该集中分布频带F进行滤波后按照公式(10)逐点错动计算两列数据W和H的互相关系数ρxy:
式中x(l)和y(l)分别是数据W和模板数据H的数据序列,d表示数据序列的个数,和分别是x(l)和y(l)的均值;
S7、选取逐点错动计算得到的互相关系数ρxy序列中的最大值作为两段数据的最终相关系数结果,当最终相关系数大于等于设定的相关系数阈值时,即判定该特定地区发生核爆炸地震,则启动报警;
S8、当最终相关系数小于设定的相关系数阈值时,重复进行步骤S3-S7,实现对特定地区单个台站的核爆炸地震信号进行检测。
由于特定地区一般设定为某国家的核试验场,这些核试验场一般不处于地震带上,极少会发生天然地震,因此可以推断出该特定地区发生的地震事件极可能为核爆炸地震事件。
采用本发明的方法对某特定地区的核爆炸地震事件进行了信号相关报警测试。先截取了历史核爆炸地震事件在选定台站上的信号记录,并进行了时频域频谱分析确定了滤波频带,在选定的台站的实时数据检测中发现检测信号后,截取了实时波形数据与历史数据进行信号相关计算并给出了事件报警。通过对比图1和图2的天然地震信号和核爆炸地震信号相关计算结果,采用本方法得到的天然地震事件最大相关系数为0.11,核爆炸地震事件最大相关系数为0.85。可以看出,利用本发明可以对特定地区的天然地震事件进行了过滤,实现了针对特定地区核爆炸地震事件的单台检测并报警,完善了对于特定地区核爆炸事件的监测能力。
Claims (4)
1.一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于,包括以下步骤:
S1、选取一个地震监测台站作为信号相关报警台站;
S2、读取特定地区历史核爆炸地震事件在选定台站上的监测数据,使用公式(5)进行时频域频谱分析得到信号时频域谱图:
S=spectrogram(x,window,noverlap,nfft,fs)(5)
式中x为监测台站记录的数据,window为滑动窗口内的数据点数,noverlap为计算过程中滑动窗口重叠的数据点数,nfft为快速傅里叶变换的数据长度,fs为监测台站的采样率,spectrogram为时频域频谱分析函数;
观察信号时频域谱图确认信号时间,截取该时间前后指定时间长度L的数据作为信号相关处理的模板数据H;根据信号时频域谱图确定能量集中分布频带F;
S3、获取监测台站实时监测数据,计算信号检测时刻短时绝对平均值STAj和长时绝对平均值LTAj:
式中Ak表示是采样点序号k的数值,j表示信号检测时刻采样点序号,n表示短时绝对平均值时间窗口数据点数,m表示长时绝对平均值时间窗口数据点数;
S4、比较台站数据的STAj/LTAj结果与设定的信号检测阈值Th大小,当
根据比较结果和台站数据信息给出台站的信号检测点时间Tdet;
S5、截取该台站信号检测点时间Tdet前后指定时间长度为L的数据W,结合步骤S2得到的模板数据H及能量集中分布频带F,同时在该集中分布频带F进行滤波后按照公式(10)逐点错动计算两列数据W和H的互相关系数ρxy:
式中x(l)和y(l)分别是数据W和模板数据H的数据序列,d表示数据序列的个数,和分别是x(l)和y(l)的均值;
S6、选取逐点错动计算得到的互相关系数ρxy序列中的最大值作为两段数据的最终相关系数结果,当最终相关系数大于等于设定的相关系数阈值时,即判定该特定地区发生核爆炸地震事件;
S7、当最终相关系数小于设定的相关系数阈值时,重复进行步骤S3-S6,实现对特定地区单个台站的核爆炸地震信号进行检测。
2.根据权利要求1所述的一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于,步骤S1具体包括以下步骤:
S1.1、按照公式(1)计算第i个台站到特定地区范围内的最大距离Di:
Di=max(distance(stai,range))(1)
式中stai为第i个台站的经纬度坐标,range为特定地区经纬度范围,distance为台站和特定地区之间距离计算函数,max为取最大值函数;
S1.2、选取距离在给定距离D内的所有台站,即每个台站的距离符合以下要求:
Di<D(2)
S1.3、统计给定时间段范围内特定地区范围内发生的事件总数为N,1.2)步骤中选取的台站中第i个台站检测到的事件数为Ki,其中1<Ki≤N,按照公式(3)计算第i个台站的监测灵敏度Ei:
S1.4、选取灵敏度在给定阈值E内的一个台站,符合以下要求,作为信号相关报警台站:
Ei>E(4)。
3.根据权利要求2所述的一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于:指定时间长度L为10s-20s。
4.根据权利要求1-3任一项所述的一种采用单个台站检测特定地区核爆炸地震事件的方法,其特征在于:所述特定地区为某国家的核试验场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610064683.6A CN105572738B (zh) | 2016-01-29 | 2016-01-29 | 一种采用单个台站检测特定地区核爆炸地震事件的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610064683.6A CN105572738B (zh) | 2016-01-29 | 2016-01-29 | 一种采用单个台站检测特定地区核爆炸地震事件的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105572738A true CN105572738A (zh) | 2016-05-11 |
CN105572738B CN105572738B (zh) | 2016-11-30 |
Family
ID=55883079
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610064683.6A Active CN105572738B (zh) | 2016-01-29 | 2016-01-29 | 一种采用单个台站检测特定地区核爆炸地震事件的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105572738B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107180512A (zh) * | 2017-06-22 | 2017-09-19 | 禁核试北京国家数据中心 | 一种特定地区地震事件的报警方法 |
CN107272061A (zh) * | 2017-06-29 | 2017-10-20 | 禁核试北京国家数据中心 | 一种次声信号与地震事件的自动关联方法 |
CN107290787A (zh) * | 2017-06-29 | 2017-10-24 | 禁核试北京国家数据中心 | 一种地震次声同址台站的监测信号关联方法 |
CN107300715A (zh) * | 2017-06-22 | 2017-10-27 | 禁核试北京国家数据中心 | 一种识别核爆炸地震事件的方法 |
CN109709598A (zh) * | 2018-12-11 | 2019-05-03 | 中国铁道科学研究院集团有限公司 | 一种基于伪源衰减法的地震误报判识方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033844A (zh) * | 2012-12-12 | 2013-04-10 | 中国地震局地震研究所 | 单站地震p波检测预警仪 |
CN104459784A (zh) * | 2014-12-11 | 2015-03-25 | 中国科学院地质与地球物理研究所 | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 |
-
2016
- 2016-01-29 CN CN201610064683.6A patent/CN105572738B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033844A (zh) * | 2012-12-12 | 2013-04-10 | 中国地震局地震研究所 | 单站地震p波检测预警仪 |
CN104459784A (zh) * | 2014-12-11 | 2015-03-25 | 中国科学院地质与地球物理研究所 | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 |
Non-Patent Citations (5)
Title |
---|
AGIUS R. MATTHEW 等: "A Single-Station Automated Earthquake Location System at Wied Dalam Station, Malta", 《SEISMOLOGICAL RESEARCH LETTERS》 * |
BOBROV 等: "Perspective of Cross-Correlation in Seismic Monitoring at the International Data Centre", 《PURE AND APPLIED GEOPHYSICS》 * |
FORHLICH C. 等: "地震的单台定位:为深入研究而进行的考察和讨论", 《世界地震译丛》 * |
PULLIAM J. 等: "单台定位的精度和不确定性", 《世界地震译丛》 * |
刘贤伦: "单台识别核爆的方法初探", 《地震地磁观测与研究》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107180512A (zh) * | 2017-06-22 | 2017-09-19 | 禁核试北京国家数据中心 | 一种特定地区地震事件的报警方法 |
CN107300715A (zh) * | 2017-06-22 | 2017-10-27 | 禁核试北京国家数据中心 | 一种识别核爆炸地震事件的方法 |
CN107300715B (zh) * | 2017-06-22 | 2018-12-11 | 禁核试北京国家数据中心 | 一种识别核爆炸地震事件的方法 |
CN107272061A (zh) * | 2017-06-29 | 2017-10-20 | 禁核试北京国家数据中心 | 一种次声信号与地震事件的自动关联方法 |
CN107290787A (zh) * | 2017-06-29 | 2017-10-24 | 禁核试北京国家数据中心 | 一种地震次声同址台站的监测信号关联方法 |
CN107290787B (zh) * | 2017-06-29 | 2018-12-11 | 禁核试北京国家数据中心 | 一种地震次声同址台站的监测信号关联方法 |
CN107272061B (zh) * | 2017-06-29 | 2019-02-05 | 禁核试北京国家数据中心 | 一种次声信号与地震事件的自动关联方法 |
CN109709598A (zh) * | 2018-12-11 | 2019-05-03 | 中国铁道科学研究院集团有限公司 | 一种基于伪源衰减法的地震误报判识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105572738B (zh) | 2016-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105572738A (zh) | 一种采用单个台站检测特定地区核爆炸地震事件的方法 | |
CN102606891B (zh) | 漏水检测仪及检测系统、及检测方法 | |
CN102338830B (zh) | 闪电探测定位系统及方法 | |
CN102563361A (zh) | 基于常规数据和音波信号的输气管道泄漏检测定位装置及方法 | |
JP2013525790A (ja) | 配管漏れをチェックするデータ処理方法およびシステム | |
CN105676287B (zh) | 一种检测特定地区核爆炸地震事件的方法 | |
CN104373820B (zh) | 降低管道泄漏监测误报率的方法 | |
CN106382981B (zh) | 一种单站次声波信号识别提取方法 | |
CN105676286A (zh) | 一种用于地震预警系统的实时地震震级估计方法 | |
CN105467428A (zh) | 一种基于短时能量检测和频谱特征分析的地震波预警方法 | |
Liu et al. | A sound monitoring system for prevention of underground pipeline damage caused by construction | |
Vogel et al. | Brief Communication: Likelihood of societal preparedness for global change: trend detection | |
CN103499831A (zh) | 一种地震数据监控系统 | |
CN105067101A (zh) | 振源识别的基于振动信号的基音频率特征的提取方法 | |
Liguori et al. | Estimation of the minimum measurement time interval in acoustic noise | |
CN104113872A (zh) | 一种数据业务监控方法及系统 | |
CN110448849A (zh) | 一种基于机器学习的消火栓压力监测系统及方法 | |
CN202093101U (zh) | 闪电探测定位系统 | |
Roarty et al. | Automated quality control of High Frequency radar data | |
CN104048165B (zh) | 管道泄漏诊断的方法 | |
CN105738950B (zh) | 一种针对特定地区的指向性聚束检测方法 | |
CN104105045A (zh) | 响度检测方法及系统 | |
CN113050158A (zh) | 近场微震信号波形的分析方法、装置、设备及存储介质 | |
CN106322123B (zh) | 管道测漏方法、系统及其采样器和主机 | |
Picozzi et al. | Trans-national earthquake early warning (EEW) in north-eastern Italy, Slovenia and Austria: first experience with PRESTo at the CE 3 RN network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |