CN105572738B - 一种采用单个台站检测特定地区核***地震事件的方法 - Google Patents
一种采用单个台站检测特定地区核***地震事件的方法 Download PDFInfo
- Publication number
- CN105572738B CN105572738B CN201610064683.6A CN201610064683A CN105572738B CN 105572738 B CN105572738 B CN 105572738B CN 201610064683 A CN201610064683 A CN 201610064683A CN 105572738 B CN105572738 B CN 105572738B
- Authority
- CN
- China
- Prior art keywords
- station
- data
- particular locality
- time
- signal
- 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.)
- Active
Links
- 238000004880 explosion Methods 0.000 title claims abstract description 38
- 238000001514 detection method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012544 monitoring process Methods 0.000 claims abstract description 31
- 238000012937 correction Methods 0.000 claims abstract description 12
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 11
- 238000012790 confirmation Methods 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 3
- 230000000052 comparative effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
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. for interpretation or for event detection
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
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,步骤S1.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 CN105572738A (zh) | 2016-05-11 |
CN105572738B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107300715B (zh) * | 2017-06-22 | 2018-12-11 | 禁核试北京国家数据中心 | 一种识别核***地震事件的方法 |
CN107180512B (zh) * | 2017-06-22 | 2018-12-11 | 禁核试北京国家数据中心 | 一种特定地区地震事件的报警方法 |
CN107272061B (zh) * | 2017-06-29 | 2019-02-05 | 禁核试北京国家数据中心 | 一种次声信号与地震事件的自动关联方法 |
CN107290787B (zh) * | 2017-06-29 | 2018-12-11 | 禁核试北京国家数据中心 | 一种地震次声同址台站的监测信号关联方法 |
CN109709598B (zh) * | 2018-12-11 | 2020-12-11 | 中国铁道科学研究院集团有限公司 | 一种基于伪源衰减法的地震误报判识方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033844A (zh) * | 2012-12-12 | 2013-04-10 | 中国地震局地震研究所 | 单站地震p波检测预警仪 |
CN104459784B (zh) * | 2014-12-11 | 2017-05-24 | 中国科学院地质与地球物理研究所 | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 |
-
2016
- 2016-01-29 CN CN201610064683.6A patent/CN105572738B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105572738A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105572738B (zh) | 一种采用单个台站检测特定地区核***地震事件的方法 | |
CN105676287B (zh) | 一种检测特定地区核***地震事件的方法 | |
CN107276707B (zh) | 一种基于多属性分析的黑广播自动化分析方法及智能化识别装置 | |
CN111060965B (zh) | 一种基于卷积神经网络的地震震相拾取及事件检测方法 | |
US10026418B2 (en) | Abnormal frame detection method and apparatus | |
CN106683687B (zh) | 异常声音的分类方法和装置 | |
CN105676286B (zh) | 一种用于地震预警***的实时地震震级估计方法 | |
CN109409308A (zh) | 一种基于鸟类鸣声的鸟类物种识别的方法 | |
CN106382981B (zh) | 一种单站次声波信号识别提取方法 | |
CN101430882A (zh) | 一种抑制风噪声的方法及装置 | |
CN105716707B (zh) | 一种超低频次声异常信号判别方法 | |
CN105357228A (zh) | 一种基于动态阈值的突发流量检测方法 | |
CN106205606A (zh) | 一种基于语音识别的动态定位监控方法及*** | |
CN105067101A (zh) | 振源识别的基于振动信号的基音频率特征的提取方法 | |
CN106023499A (zh) | 一种光纤安防信号双重识别方法及*** | |
CN107993304A (zh) | 一种数据处理方法、装置、设备及计算机可读存储介质 | |
CN105467428A (zh) | 一种基于短时能量检测和频谱特征分析的地震波预警方法 | |
CN109116196A (zh) | 一种电力电缆故障放电声智能识别方法 | |
CN105738950B (zh) | 一种针对特定地区的指向性聚束检测方法 | |
Lara et al. | A real-time microearthquakes-detector based on voice activity detection and endpoint detection: An approach to Cotopaxi Volcano | |
CN107300715B (zh) | 一种识别核***地震事件的方法 | |
CN104597486A (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 | |
Zhuang et al. | Preliminary analysis of observations on the ultra-low frequency electric field in the Beijing region | |
CN108037529A (zh) | 一种基于初始振动信号的地震事件快速识别方法 |
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 |