CN108562933A - 一种融合多类数据源的震级快速估计方法 - Google Patents
一种融合多类数据源的震级快速估计方法 Download PDFInfo
- Publication number
- CN108562933A CN108562933A CN201810359680.4A CN201810359680A CN108562933A CN 108562933 A CN108562933 A CN 108562933A CN 201810359680 A CN201810359680 A CN 201810359680A CN 108562933 A CN108562933 A CN 108562933A
- Authority
- CN
- China
- Prior art keywords
- magnitude
- earthquake
- gps
- formula
- station
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000004927 fusion Effects 0.000 title claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 6
- 241001269238 Data Species 0.000 claims description 4
- 230000001960 triggered effect Effects 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 7
- 230000002045 lasting effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000116 mitigating effect Effects 0.000 description 2
- 235000001167 Melaleuca cajuputi Nutrition 0.000 description 1
- 244000304222 Melaleuca cajuputi Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 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/01—Measuring or predicting earthquakes
-
- 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
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/6161—Seismic or acoustic, e.g. land or sea measurements
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
本发明公开了一种融合多类数据源的震级快速估计方法。针对多种预警震级计算方法适用范围的不同,提出联合利用高频GPS、强震动及测震多类型数据快速测定震级的方法。解决了仅使用单一方法测定震级时容易出现误报的问题,以及利用P波初始记录测定震级时出现震级的饱和问题,实现地震震级的动态、持续测定,最大程度的提高了震级估计的准确性,能够应用于地震预警系统,从而有效的降低系统的误报率。
Description
技术领域
本发明涉及防震减灾技术领域,特别是一种融合多类数据源的震级快速估计方法。
背景技术
我国是世界上多地震国家之一,近几年强震频发,2008年5月12日汶川8.0级地震、2010年4月14日玉树7.1级地震和2013年4月20日芦山7.0级地震的发生,对当地的经济社会发展产生了巨大的影响,对人民生命和财产造成了严重的损失。
地震预报没有实现重大突破的今天,地震预警系统是世界上公认的能够有效减轻地震灾害的新技术手段之一。地震预警系统利用地震波传播速度远小于电磁波的原理,通过快速分析地震监测台网实时数据,在具有破坏性的S波到达预警目标区之前发布地震预警信息,提醒人们紧急避险,同时对重要工程设施进行紧急处置,从而达到减轻地震灾害的目的。
震级(Ms)的快速估计是地震预警系统中最关键的部分,同时也是最困难的技术环节。预警震级估计的准确性直接影响到极震区及预警区的烈度预测,从而关乎是否对预警区发布预警信息的策略选择,是地震预警系统产生误报和漏报的重要因素。然而,由于预警信息发布的实时性要求,预警震级的计算仅能使用最先触发的几个台站初至波数据,传统常规的利用多个台站的面波信息计算震级的方法并不适用于地震预警系统,因此震级的快速估计也成为地震预警系统中最困难的技术环节。
目前,国际上比较常用的测定实时震级的方法主要分为两种,即利用测震或强震动台站记录的P波前几秒数据的卓越周期计算震级的τc方法和量取P波前几秒数据最大幅值计算震级的Pd方法。两种方法均可以在地震破裂过程没有结束之前,快速估计整个地震的震级规模。从τc和Pd方法适用范围来看,独立使用一种方法测定震级存在一定的局限性。
(1)τc方法更适用于震级偏高的地震事件,对于中强以下的地震(Ms≤4.5)不敏感。Pd方法受到震中距的影响,对于部分震中台站较少记录的地震事件,测定的震级信度较低。
(2)对于计算的较大地震(Ms≥6.5),两种方法都会出现震级饱和的现象,计算得到的震级偏低,识别率不稳定,容易出现误报的问题。
发明内容
为了解决上述问题,本发明提供一种融合多类数据源的震级快速估计方法。
为实现上述目的,本发明提供如下技术方案:一种融合多类数据源的震级快速估计方法,利用测震、强震动及高频GPS实时数据,当检测到第一个站点地震事件后,进入以下地震震级测定步骤:
1、阈值触发判断
首先对测震及强震动台站产出实时数据流进行地震事件触发判断,在满足三个触发台站的空间集中性后,对测震与强震动垂直向数据分别进行一次和二次积分得到位移记录,并利用Butterworth高通滤波器对数据进行滤波处理,检测触发的三个台站中是否有一个台站同时满足τc和Pd阈值(对于Ms=3.0级地震,设定τc和Pd的阈值分别为0.2s和5μm),τc为记录数据的卓越周期,Pd为记录数据的最大幅值,τc计算公式见式(1):
式中,u和分别代表高通滤波后台站垂直向地动位移和速度;
如果有一个台站达到阈值条件,将锁定此事件为一个需要发布地震预警的地震事件,随后计算各个触发台站的τc和Pd;
2、逐秒震级计算
当锁定一个地震事件后,对触发台站的实时数据流中的每一秒的数据包进行分析,计算τc和Pd,去除超出一倍标准差的异常数据,计算预测震级,是根据触发各个台站τc的平均值计算得到的震级,计算公式见式(2):
Mpd是根据震中距80km以内的所有触发台站的最大Pd值计算得到的震级,计算公式见式(3):
式中,Δ为震中距;
取和Mpd均值作为实时震级Mre;
3、GPS震级校正
当Mre≥6.5时,搜索震中附近是否有GPS台站触发,如有GPS台站触发,获取GPS台站的最大地面位移,由震中距得量规函数值,计算各个台站的GPS震级MGPS,MGPS计算公式见式(4):
MGPS=Log(A/T)+1.66Log(Δ)+3.3 (4)
式中,A为绝对位移的最大振幅,T为最大振幅对应的周期,Δ为震中距;
取各个台站MGPS均值对原发布震级进行更新校正;
4、震级最终确定
当所有触发台站的数据长度满3s,将计算得到的各个GPS台站MGPS均值作为最终震级Mre发布。
优选地,式(2)与τc的统计关系由中国大陆247条地震记录计算得到。
优选地,式(3)与Pd的统计关系由中国大陆247条地震记录计算得到。
与现有技术相比,本发明有益效果:
本发明针对多种预警震级计算方法适用范围的不同,提出联合利用高频GPS、强震动及测震多类型数据快速测定震级的方法。解决了仅使用单一方法测定震级时容易出现误报的问题,以及利用P波初始记录测定震级时出现震级的饱和问题,实现地震震级的动态、持续测定,最大程度的提高了震级估计的准确性,能够应用于地震预警系统,从而有效的降低系统的误报率。
附图说明
图1是本发明震级快速估计方法流程图。
图2是2013年4月20日四川芦山7.0级地震震级快速估计过程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供一种融合多类数据源的震级快速估计方法,包括以下步骤:
步骤1:阈值触发判断
首先对测震及强震动台站产出实时数据流进行地震事件触发判断,在满足三个触发台站的空间集中性后,对测震与强震动垂直向数据分别进行一次和二次积分得到位移记录,并利用Butterworth高通滤波器对数据进行滤波处理,检测触发的三个台站中是否有一个台站同时满足τc和Pd阈值(对于Ms=3.0级地震,设定τc和Pd的阈值分别为0.2s和5μm);
τc为记录数据的卓越周期,Pd为记录数据的最大幅值,τc计算公式见式(1):
式中,u和分别代表高通滤波后台站垂直向地动位移和速度;
如果至少有一个台站达到阈值条件,将锁定此事件为一个需要发布地震预警的地震事件,随后计算各个触发台站的τc和Pd;
步骤2:逐秒震级计算
当锁定一个地震事件后,对触发台站的实时数据流中的每一秒的数据包进行分析,计算全部触发台站的τc,去除超出一倍标准差的异常数据,取τc平均值,计算预测震级计算公式见式(2):
计算震中距80km以内的所有触发台站的最大Pd值,计算预测震级Mpd,计算公式见式(3):
式中,Δ为震中距;
取和Mpd均值作为实时震级Mre;
步骤3:GPS震级校正
当Mre≥6.5时,搜索震中附近是否有GPS台站触发,如有GPS台站触发,获取GPS台站的最大地面位移,由震中距得量规函数值,计算各个台站的GPS震级MGPS,MGPS计算公式见式(4):
MGPS=Log(A/T)+1.66Log(Δ)+3.3 (4)
式中,A为绝对位移的最大振幅,T为最大振幅对应的周期,Δ为震中距;
取各个台站MGPS均值对原发布震级进行更新校正;
步骤4:震级最终确定
当所有触发台站的数据长度满3s,将计算得到的各个GPS台站MGPS均值作为最终震级Mre发布。
图2是2013年4月20日四川芦山7.0级地震震级快速估计过程图,图中4行数据分别为距离震中最近的宝兴台测震仪垂直向速度记录、距离震中最近的宝兴台强震仪垂直向加速度记录、距离震中最近的芦山站GPS三分项记录以及实时震级变化过程。图中可以清晰的看出,实时震级在测震与强震台站触发后,开始逐秒计算,计算结果为与的均值。第一秒计算震级为6.2,第二秒计算震级为7.6,第三秒计算震级为7.5。此时,震中附近GPS台站震级计算完成,计算震级为6.9。第四秒将作为最终震级进行更新修正。
显然,本领域的技术人员可以对发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种融合多类数据源的震级快速估计方法,其特征在于,利用测震、强震动及高频GPS实时数据,当检测到第一个站点地震事件后,进入以下地震震级测定步骤:
1)、阈值触发判断
首先对测震及强震动台站产出实时数据流进行地震事件触发判断,在满足三个触发台站的空间集中性后,对测震与强震动垂直向数据分别进行一次和二次积分得到位移记录,并利用Butterworth高通滤波器对数据进行滤波处理,检测触发的三个台站中是否有一个台站同时满足τc和Pd阈值,τc为记录数据的卓越周期,Pd为记录数据的最大幅值,τc计算公式见式(1):
式中,u和分别代表高通滤波后台站垂直向地动位移和速度;
如果有一个台站达到阈值条件,将锁定此事件为一个需要发布地震预警的地震事件,随后计算各个触发台站的τc和Pd;
2)、逐秒震级计算
当锁定一个地震事件后,对触发台站的实时数据流中的每一秒的数据包进行分析,计算τc和Pd,去除超出一倍标准差的异常数据,计算预测震级。是根据触发各个台站τc的平均值计算得到的震级,计算公式见式(2):
Mpd是根据震中距80km以内的所有触发台站的最大Pd值计算得到的震级,计算公式见式(3):
式中,Δ为震中距;
取和Mpd均值作为实时震级Mre;
3)、GPS震级校正
当Mre≥6.5时,搜索震中附近是否有GPS台站触发,如有GPS台站触发,获取GPS台站的最大地面位移,由震中距得量规函数值,计算各个台站的GPS震级MGPS,MGPS计算公式见式(4):
MGPS=Log(A/T)+1.66Log(Δ)+3.3 (4)
式中,A为绝对位移的最大振幅,T为最大振幅对应的周期,Δ为震中距;
取各个台站MGPS均值对原发布震级进行更新校正;
4)、震级最终确定
当所有触发台站的数据长度满3s,将计算得到的各个GPS台站MGPS均值作为最终震级Mre发布。
2.根据权利要求1所述的震级快速估计方法,其特征在于,所述τc和Pd的阈值设定为MS=3.0级地震时的τc和Pd值,分别为0.2s和5μm。
3.根据权利要求1所述的震级快速估计方法,其特征在于,所述式(2)中与τc的统计关系由中国大陆247条地震记录计算得到。
4.根据权利要求1所述的震级快速估计方法,其特征在于,所述式(3)中与Pd的统计关系由中国大陆247条地震记录计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359680.4A CN108562933B (zh) | 2018-04-20 | 2018-04-20 | 一种融合多类数据源的震级快速估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359680.4A CN108562933B (zh) | 2018-04-20 | 2018-04-20 | 一种融合多类数据源的震级快速估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108562933A true CN108562933A (zh) | 2018-09-21 |
CN108562933B CN108562933B (zh) | 2019-08-16 |
Family
ID=63535902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810359680.4A Active CN108562933B (zh) | 2018-04-20 | 2018-04-20 | 一种融合多类数据源的震级快速估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108562933B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507728A (zh) * | 2018-11-08 | 2019-03-22 | 长沙矿山研究院有限责任公司 | 一种基于微震监测的地压灾害预警方法 |
CN109581479A (zh) * | 2018-12-10 | 2019-04-05 | 南京云创大数据科技股份有限公司 | 一种地震预警信息的处理方法 |
CN111552217A (zh) * | 2020-05-26 | 2020-08-18 | 河南安达盛电子技术有限公司 | 一种地震台站运维系统及方法 |
CN116699691A (zh) * | 2023-06-14 | 2023-09-05 | 中国石油大学(华东) | 一种基于强震动pgv的震级快速估算方法、系统及终端 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002277557A (ja) * | 2001-03-21 | 2002-09-25 | Railway Technical Res Inst | 震央距離及びマグニチュード推定方法とそのための装置 |
JP2006170739A (ja) * | 2004-12-15 | 2006-06-29 | Kajima Corp | 緊急地震速報を用いた地震防災システム |
JP2010091480A (ja) * | 2008-10-09 | 2010-04-22 | Ohbayashi Corp | 長周期地震動検知器、長周期地震動検知方法 |
CN102736101A (zh) * | 2012-06-19 | 2012-10-17 | 成都高新减灾研究所 | 一种分布式地震预警信息处理方法及系统 |
CN104459784A (zh) * | 2014-12-11 | 2015-03-25 | 中国科学院地质与地球物理研究所 | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 |
-
2018
- 2018-04-20 CN CN201810359680.4A patent/CN108562933B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002277557A (ja) * | 2001-03-21 | 2002-09-25 | Railway Technical Res Inst | 震央距離及びマグニチュード推定方法とそのための装置 |
JP2006170739A (ja) * | 2004-12-15 | 2006-06-29 | Kajima Corp | 緊急地震速報を用いた地震防災システム |
JP2010091480A (ja) * | 2008-10-09 | 2010-04-22 | Ohbayashi Corp | 長周期地震動検知器、長周期地震動検知方法 |
CN102736101A (zh) * | 2012-06-19 | 2012-10-17 | 成都高新减灾研究所 | 一种分布式地震预警信息处理方法及系统 |
CN104459784A (zh) * | 2014-12-11 | 2015-03-25 | 中国科学院地质与地球物理研究所 | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 |
Non-Patent Citations (5)
Title |
---|
GEOFFREY BLEWITT,ET AL.: "Rapid determination of earthquake magnitude using GPS for tsunami warning systems", 《GEOPHYSICAL RESEARCH LETTERS》 * |
YIH-MIN WU,ET AL.: "Development of an Earthquake Early Warning System Using Real-Time Strong Motion Signals", 《SENSORS》 * |
刘希强,等: "地震预警系统中进一步提高三台震中定位精度的方法研究", 《地震研究》 * |
刘希强,等: "山东省地震动强度(烈度)实时速报系统研究", 《中国地球物理2010》 * |
曲保安,等: "用P波前3秒估算预警震级的源误差分析及消除", 《中国地球科学联合学术年会 2014》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507728A (zh) * | 2018-11-08 | 2019-03-22 | 长沙矿山研究院有限责任公司 | 一种基于微震监测的地压灾害预警方法 |
CN109581479A (zh) * | 2018-12-10 | 2019-04-05 | 南京云创大数据科技股份有限公司 | 一种地震预警信息的处理方法 |
CN109581479B (zh) * | 2018-12-10 | 2020-11-06 | 南京云创大数据科技股份有限公司 | 一种地震预警信息的处理方法 |
CN111552217A (zh) * | 2020-05-26 | 2020-08-18 | 河南安达盛电子技术有限公司 | 一种地震台站运维系统及方法 |
CN116699691A (zh) * | 2023-06-14 | 2023-09-05 | 中国石油大学(华东) | 一种基于强震动pgv的震级快速估算方法、系统及终端 |
CN116699691B (zh) * | 2023-06-14 | 2024-11-05 | 中国石油大学(华东) | 一种基于强震动pgv的震级快速估算方法、系统及终端 |
Also Published As
Publication number | Publication date |
---|---|
CN108562933B (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108562933B (zh) | 一种融合多类数据源的震级快速估计方法 | |
US11783691B1 (en) | Beidou communication-based integrated monitoring system for slope surface crack | |
US7855640B2 (en) | Method and system of determining alarm conditions | |
KR102468888B1 (ko) | IoT 기반 건물 상태 모니터링 시스템 | |
CN102385787B (zh) | 一种区域地震监测网的预警方法 | |
CN116129612B (zh) | 一种基于物联网的分布式滑坡监测预警系统 | |
CN114495429B (zh) | 一种崩塌灾害早期预警方法及系统 | |
KR101763337B1 (ko) | 진동 가속도와 변위 계측 기반 재난 경보 시스템 및 방법 | |
Nakamura | A new concept for the earthquake vulnerability estimation and its application to the early warning system | |
JP2019067217A (ja) | 安全度判定プログラム、安全度判定方法および安全度判定装置 | |
JP2023550091A (ja) | 分散型光ファイバセンシングによる振動の垂直距離予測 | |
CN118011269A (zh) | 一种基于北斗的智能接地线系统及方法 | |
JP3296764B2 (ja) | 光学式振り子地震センサとそれを用いた防災警報装置及び防災警報判定法 | |
CN116989679A (zh) | 一种基于图像处理的高速公路高边坡位移集成监测方法 | |
Bressan et al. | Detecting the 11 March 2011 Tohoku tsunami arrival on sea-level records in the Pacific Ocean: application and performance of the Tsunami Early Detection Algorithm (TEDA) | |
CN117894157B (zh) | 一种基于大数据的地震灾害预警系统 | |
CN207050704U (zh) | 一种可用于建筑结构形变监测的位移监测系统 | |
CN109139019A (zh) | 潜埋暗挖隧道下穿地铁既有线路的监测方法及监测装置 | |
JP2004060311A (ja) | 地すべり監視方法、およびシステム | |
CN106157539A (zh) | 地震临震报警系统 | |
JP6326155B1 (ja) | 地盤傾斜計による地盤傾斜確認方法 | |
US20230296798A1 (en) | Rock fall analyser | |
CN115953886A (zh) | 灾害预警方法、装置、系统和电子设备 | |
CN116241327A (zh) | 一种隧道掌子面稳定性监测方法及系统 | |
CN107101610A (zh) | 可用于建筑结构形变监测的位移监测系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |