CN108089225A - 一种基于单台站初至p波的震级快速估算方法 - Google Patents
一种基于单台站初至p波的震级快速估算方法 Download PDFInfo
- Publication number
- CN108089225A CN108089225A CN201711255119.3A CN201711255119A CN108089225A CN 108089225 A CN108089225 A CN 108089225A CN 201711255119 A CN201711255119 A CN 201711255119A CN 108089225 A CN108089225 A CN 108089225A
- Authority
- CN
- China
- Prior art keywords
- data
- formula
- magnitude
- earthquake
- maximum
- 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 21
- 230000033001 locomotion Effects 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 5
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 238000003786 synthesis reaction Methods 0.000 abstract description 2
- 230000007812 deficiency Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000011435 rock Substances 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/24—Recording seismic data
-
- 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
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
一种基于单台站初至P波的震级快速估算方法,通过对检测到的初始1.5秒P波信号进行时域和频域分析,得到由时域幅值参数和频域周期参数组成的综合多种参数的震级估算公式,利用日本强震动记录和差分进化算法回归震级估算公式中各参数的系数,给出震级的最终估算公式。本发明方法仅需单台站初至1.5秒的P波数据估算震级,具有良好的准确性、时效性和易用性,适用于地震监测和地震预警。
Description
技术领域
本发明涉及一种基于单台站初至P波的震级快速估算方法,主要用于利用初始的地震P波数据估算震级。
背景技术
地震灾害可造成人员伤亡、毁坏人工结构物、又可引发火灾、洪水、泥石流、海啸、核泄漏等次生灾害,严重危害社会经济发展。尽管地震预报至今仍是世界性的难题,但地震预警作为震灾预防的新技术,在很多地震多发国家和地区引起重视并得到应用。地震预警技术是指在地震发生后,在破坏性的地震波到达设防地区前,快速确定地震参数(震级、震中距、震源深度等),对设防地区发出报警的一种技术。地震预警技术分为多台站地震预警和单台站地震预警,多台站地震预警是利用多个监测台站确定地震参数,可用数据多,准确性较高,但时效性差;单台站地震预警是利用单个监测台站确定地震参数,可用数据少,准确性较差,但时效性好。单台地震预警技术能够为设防地区提供更长的预警时间,因此如何提高单台站地震预警的准确性,是地震预警技术的一个重要研究内容。震级的准确性至关重要,决定了地震影响的大小和范围,当前地震预警中震级估算普遍采用Pd法、τc法、τp法等,这些方法的不足在于:1)小地震的震级估算结果偏大,大地震的震级估算结果震级饱和,准确性低;2)需要利用初至3秒以上的P波数据进行计算,在时效性上仍难以满足实时的地震预警***需求。因此,发展一种准确、高效、易用的单台站震级快速估算方法成为地震预警的一个亟需解决的难题。
发明内容
针对震级估算方法的前述现状和不足,本发明提供了一种基于单台站初至P波的震级快速估算方法,通过对检测到的初始1.5秒P波信号进行时域和频域分析,得到由时域幅值参数和频域周期参数组成的综合多种参数的震级估算公式,利用日本强震动记录和差分进化算法回归震级估算公式中各参数的系数,给出震级的最终估算公式。本发明方法仅需单台站初至1.5秒的P波数据估算震级,具有良好的准确性、时效性和易用性,适用于地震监测和地震预警。
本发明为实现上述目的所采用的技术方案是:
1)计算初始1.5秒P波数据的时域参数:求出三个方向加速度数据的合成加速度(平方和开根号)的最大值PGA;求出竖向加速度的最大值ZPGA;对竖向加速度数据进行积分得到速度数据和位移数据,求出速度最大值ZPGV和位移最大值ZPGD;对竖向加速度数据的绝对值进行积分得到累积绝对速度数据,求出最大值CAV。
2)计算初始1.5秒P波数据的频域参数:按公式(1)计算最大卓越周期TPmax;按公式(2)计算平均卓越周期Tc。
公式(1)和公式(2)中,v是速度数据,u是位移数据,t数据时程时间,i是离散数据点序号,N离散数据点总数,T数据总时间长。
3)通过1)和2)的所有求参数,建立震级公式(3)。
公式(4)中,Mag是待估算震级,a1,a2,a3,a4,a5,a6,a7,a8,a9是待回归系数。
5)利用日本NIED强震动记录(K-NET,Kik-net:震中距20km~200km,震级大于3.5级,震源深度小于100km,P波到时明确的记录)的初始1.5秒P波数据和差分进化算法对公式(3)进行系数回归,得到震级估算经验公式(4)
本发明方法的优点:
1)本发明方法综合多种时域参数和频域参数建立的震级估算公式,兼顾了多方面信息,并利用日本NIED的高质量强震动数据和差分进化算法回归系数,可大幅提高震震级估算的准确性高;
2)本发明方法使震级估算时间由3秒降低到1.5秒,满足实时或秒级的地震数据处理***的要求,时效性好;
3)本发明方法原理简单,易于编程实现,计算效率高,可广泛地应用地震监测和地震预警的各种软硬件***中,易用性好。
附图说明
图1为一种基于单台站初至P波的震级快速估算方法。
具体实施方式
实施例:
下面结合实施例对本发明做进一步说明,但本发明并不限于以下实施例。
本发明所提供的基于单台站初至P波的震级快速估算方法的具体实施流程如下:
1)三分向(东西方向、南北方向和竖直方向)加速度数据采集仪器检测到P波信号后,缓冲1.5秒P波数据,对1.5秒数据进行基线处理(基线处理方法为:1.5秒P波数据减去P波到时之前几秒内噪声数据的平均值。)。
2)计算时域参数:计算三个方向加速度数据的合成加速度(平方和开根号)的最大值PGA;计算竖向加速度的最大值ZPGA;对竖向加速度数据进行积分求出速度数据,并求出最大值ZPGV;对速度数据再次积分求出位移数据,并求出最大值ZPGD;对竖向加速度数据的绝对值进行积分得到累积绝对速度数据,求出最大值CAV。
3)计算频域参数:按公式(1)计算最大卓越周期TPmax;按公式(2)计算平均卓越周期Tc。
4)计算震级:按公式(4)计算震级。
Claims (1)
1.一种基于初始振动信号的地震事件快速识别方法,其特征在于具体步骤为:
1)计算初始1.5秒P波数据的时域参数:求出三个方向加速度数据的合成加速度(平方和开根号)的最大值PGA;求出竖向加速度的最大值ZPGA;对竖向加速度数据进行积分得到速度数据和位移数据,求出速度最大值ZPGV和位移最大值ZPGD;对竖向加速度数据的绝对值进行积分得到累积绝对速度数据,求出最大值CAV。
2)计算初始1.5秒P波数据的频域参数:按公式(1)计算最大卓越周期TPmax;按公式(2)计算平均卓越周期Tc。
公式(1)和公式(2)中,v是速度数据,u是位移数据,t数据时程时间,i是离散数据点序号,N离散数据点总数,T数据总时间长。
3)通过1)和2)的所有求参数,建立震级公式(3)。
公式(4)中,Mag是待估算震级,a1,a2,a3,a4,a5,a6,a7,a8,a9是待回归系数。
5)利用日本NIED强震动记录(K-NET,Kik-net:震中距20km~200km,震级大于3.5级,震源深度小于100km,P波到时明确的记录)的初始1.5秒P波数据和差分进化算法对公式(3)进行系数回归,得到震级估算经验公式(4)
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711255119.3A CN108089225B (zh) | 2017-12-02 | 2017-12-02 | 一种基于单台站初至p波的震级快速估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711255119.3A CN108089225B (zh) | 2017-12-02 | 2017-12-02 | 一种基于单台站初至p波的震级快速估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108089225A true CN108089225A (zh) | 2018-05-29 |
CN108089225B CN108089225B (zh) | 2019-08-13 |
Family
ID=62173436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711255119.3A Active CN108089225B (zh) | 2017-12-02 | 2017-12-02 | 一种基于单台站初至p波的震级快速估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108089225B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108957537A (zh) * | 2018-06-20 | 2018-12-07 | 西安理工大学 | 倾斜地震波作用下场地卓越周期的计算方法 |
CN113009566A (zh) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
CN114355438A (zh) * | 2022-01-11 | 2022-04-15 | 西南交通大学 | 地震震级的估算方法、装置、设备及可读存储介质 |
CN115032693A (zh) * | 2022-07-11 | 2022-09-09 | 中国地震局地球物理研究所 | 一种强震的前震自动识别方法及装置 |
IT202100021416A1 (it) * | 2021-08-06 | 2023-02-06 | Cavagna Group Spa | Sistema di prevenzione |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040260475A1 (en) * | 2002-04-10 | 2004-12-23 | Cherry J. Theodore | Systems for predicting earthquakes and methods of employing such systems |
US20080183393A1 (en) * | 2007-01-25 | 2008-07-31 | Rundle John B | Method and system for forecasting earthquakes |
CN101967860A (zh) * | 2010-10-19 | 2011-02-09 | 哈尔滨工业大学 | 最不利设计地震动记录的确定方法 |
CN102628955A (zh) * | 2012-04-24 | 2012-08-08 | 成都高新减灾研究所 | 一种地震预警震级获取方法 |
CN105528520A (zh) * | 2015-12-09 | 2016-04-27 | 哈尔滨工业大学 | 一种基于主余震序列的损伤指数评估方法 |
-
2017
- 2017-12-02 CN CN201711255119.3A patent/CN108089225B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040260475A1 (en) * | 2002-04-10 | 2004-12-23 | Cherry J. Theodore | Systems for predicting earthquakes and methods of employing such systems |
US20080183393A1 (en) * | 2007-01-25 | 2008-07-31 | Rundle John B | Method and system for forecasting earthquakes |
CN101967860A (zh) * | 2010-10-19 | 2011-02-09 | 哈尔滨工业大学 | 最不利设计地震动记录的确定方法 |
CN102628955A (zh) * | 2012-04-24 | 2012-08-08 | 成都高新减灾研究所 | 一种地震预警震级获取方法 |
CN105528520A (zh) * | 2015-12-09 | 2016-04-27 | 哈尔滨工业大学 | 一种基于主余震序列的损伤指数评估方法 |
Non-Patent Citations (2)
Title |
---|
罗光财,等: "PGV/PGA和PGD/PGA随震级和震中距变化的研究", 《西北地震学报》 * |
齐娟,等: "竖向地震标准设计反应谱曲线研究", 《工程抗震与加固改造》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108957537A (zh) * | 2018-06-20 | 2018-12-07 | 西安理工大学 | 倾斜地震波作用下场地卓越周期的计算方法 |
CN113009566A (zh) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
CN113009566B (zh) * | 2021-03-19 | 2024-01-30 | 中国地震局地球物理研究所 | 现地地震动预测模型及其构建方法 |
IT202100021416A1 (it) * | 2021-08-06 | 2023-02-06 | Cavagna Group Spa | Sistema di prevenzione |
CN114355438A (zh) * | 2022-01-11 | 2022-04-15 | 西南交通大学 | 地震震级的估算方法、装置、设备及可读存储介质 |
CN115032693A (zh) * | 2022-07-11 | 2022-09-09 | 中国地震局地球物理研究所 | 一种强震的前震自动识别方法及装置 |
CN115032693B (zh) * | 2022-07-11 | 2023-02-03 | 中国地震局地球物理研究所 | 一种强震的前震自动识别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108089225B (zh) | 2019-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108089225B (zh) | 一种基于单台站初至p波的震级快速估算方法 | |
Wu | Progress on development of an earthquake early warning system using low-cost sensors | |
Husen et al. | Induced seismicity during the construction of the Gotthard Base Tunnel, Switzerland: hypocenter locations and source dimensions | |
CN104765064B (zh) | 一种微地震干涉成像的方法 | |
Cornou et al. | Contribution of dense array analysis to the identification and quantification of basin-edge-induced waves, Part II: Application to Grenoble basin (French Alps) | |
CN102628955B (zh) | 一种地震预警震级获取方法 | |
CN104732728B (zh) | 一种智能终端预警地震*** | |
CN106886044A (zh) | 一种基于剪切波与Akaike信息准则的微地震初至拾取方法 | |
CN114152980B (zh) | 一种快速自动化产出震源机制解的方法与装置 | |
Crowell et al. | Hypothetical real‐time GNSS modeling of the 2016 M w 7.8 Kaikōura earthquake: Perspectives from ground motion and tsunami inundation prediction | |
TW201819955A (zh) | 自動化校正地盤特性之現地型地震預警系統及相關方法 | |
CN108051853B (zh) | 一种基于单台站初至p波的震中距快速估算方法 | |
CN115826042A (zh) | 一种边云端结合的分布式地震数据处理方法与装置 | |
Walter et al. | Calving event detection by observation of seiche effects on the Greenland fjords | |
Tolstoy et al. | Hydroacoustic constraints on the rupture duration, length, and speed of the great Sumatra-Andaman earthquake | |
Nazaruddin et al. | Intraplate earthquake occurrence and distribution in Peninsular Malaysia over the past 100 years | |
Moya et al. | Comparison of coseismic displacement obtained from GEONET and seismic networks | |
CN110786863B (zh) | 一种基于移动设备的行人步态检测方法 | |
Li et al. | Adaptive denoising approach for high-rate GNSS seismic waveform preservation: Application to the 2010 EI Mayor-Cucapah earthquake and 2012 Brawley seismic swarm | |
CN111405469A (zh) | 基于手机移动传感网络的矿震监测***及群智定位方法 | |
Suriñach et al. | A template to obtain information on gravitational mass movements from the spectrograms of the seismic signals generated | |
Rajabi et al. | Real time earthquake prediction using cross-correlation analysis & transfer function model | |
Lin et al. | A fast magnitude estimation for the 2011 Mw 9.0 great Tohoku earthquake | |
Wu et al. | Progress on earthquake rapid reporting and early warning systems in Taiwan | |
Groos | Broadband seismic noise: classification and Green's function estimation |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211116 Address after: 518053 tower 413, block B, Tianjian technology building, Shenyun West Second Road, Gaofa community, Shahe street, Nanshan District, Shenzhen, Guangdong Patentee after: Zhongzhen Kejian (Shenzhen) Holding Co., Ltd Address before: 541004 No. 12, Jiangan Road, Qixing District, Guilin City, Guangxi Zhuang Autonomous Region Patentee before: Guilin University of Technology |