CN108089225A - 一种基于单台站初至p波的震级快速估算方法 - Google Patents

一种基于单台站初至p波的震级快速估算方法 Download PDF

Info

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
Application number
CN201711255119.3A
Other languages
English (en)
Other versions
CN108089225B (zh
Inventor
王延伟
兰景岩
曹振中
黄林
胡明祎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongzhen Kejian (Shenzhen) Holding Co., Ltd
Original Assignee
Guilin University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guilin University of Technology filed Critical Guilin University of Technology
Priority to CN201711255119.3A priority Critical patent/CN108089225B/zh
Publication of CN108089225A publication Critical patent/CN108089225A/zh
Application granted granted Critical
Publication of CN108089225B publication Critical patent/CN108089225B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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波的震级快速估算方法，主要用于利用初始的地震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)
。
CN201711255119.3A 2017-12-02 2017-12-02 一种基于单台站初至p波的震级快速估算方法 Active CN108089225B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 哈尔滨工业大学 一种基于主余震序列的损伤指数评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
罗光财，等: "PGV/PGA和PGD/PGA随震级和震中距变化的研究", 《西北地震学报》 *
齐娟，等: "竖向地震标准设计反应谱曲线研究", 《工程抗震与加固改造》 *

Cited By (7)

* Cited by examiner, † Cited by third party
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波的震级快速估算方法
Duputel et al. Real-time W phase inversion during the 2011 off the Pacific coast of Tohoku Earthquake
Husen et al. Induced seismicity during the construction of the Gotthard Base Tunnel, Switzerland: hypocenter locations and source dimensions
CN102628955B (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)
CN104765064B (zh) 一种微地震干涉成像的方法
CN106886044B (zh) 一种基于剪切波与Akaike信息准则的微地震初至拾取方法
Wu et al. A test of earthquake early warning system using low cost accelerometer in Hualien, Taiwan
CN103926623A (zh) 一种压制逆时偏移低频噪音的方法
CN104732728A (zh) 一种智能终端预警地震系统
TW201819955A (zh) 自動化校正地盤特性之現地型地震預警系統及相關方法
CN108051853B (zh) 一种基于单台站初至p波的震中距快速估算方法
CN114152980A (zh) 一种快速自动化产出震源机制解的方法与装置
Tolstoy et al. Hydroacoustic constraints on the rupture duration, length, and speed of the great Sumatra-Andaman earthquake
Sanderson et al. A pilot experiment on infrasonic lahar detection at Mount Adams, Cascades: ambient infrasound and wind‐noise characterization at a quiescent stratovolcano
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
CN111405469B (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
CN115826042A (zh) 一种边云端结合的分布式地震数据处理方法与装置
CN110786863A (zh) 一种基于移动设备的行人步态检测方法
Suriñach et al. A template to obtain information on gravitational mass movements from the spectrograms of the seismic signals generated
US20040098198A1 (en) Method of using ocean acoustic sensors for coastal earthquake prediction
Rajabi et al. Real time earthquake prediction using cross-correlation analysis & transfer function model
Hloupis et al. Wavelet-based methods for rapid calculations of magnitude and epicentral distance: an application to earthquake early warning system

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