CN105676286B - 一种用于地震预警系统的实时地震震级估计方法 - Google Patents
一种用于地震预警系统的实时地震震级估计方法 Download PDFInfo
- Publication number
- CN105676286B CN105676286B CN201610049554.XA CN201610049554A CN105676286B CN 105676286 B CN105676286 B CN 105676286B CN 201610049554 A CN201610049554 A CN 201610049554A CN 105676286 B CN105676286 B CN 105676286B
- Authority
- CN
- China
- Prior art keywords
- mrow
- earthquake
- magnitude
- seismic
- moment
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 230000005855 radiation Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims 1
- 230000001133 acceleration Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009466 transformation 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/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
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
本发明公开了一种用于地震预警系统的实时地震震级估计方法,利用地震矩求出的矩震级,即:检测到信号后每1秒估算出相对应的矩震级,随着时间增长至20s左右得出与实际震级非常接近的估计震级。本发明利用地震矩从而求出的矩震级能有效克服中大型地震(M>6)中出现的震级饱和问题,而且地震矩与地震动中零频时的振幅Ω0有明确联系,具有明确的物理意义。本发明方法具有实时性,检测到信号后每间隔1s均能估算出相对应的矩震级,随着时间增长至20s左右就能得出与实际震级非常接近的估计震级,能有效应用于地震预警系统。
Description
技术领域
本发明涉及一种新型的用于地震预警系统的地震震级的估计方法,尤其涉及一种用于克服中大型地震动中震级饱和问题,估计实时震级的谱反演方法。
背景技术
近年来,世界各地大
地震日益频繁,甚至会引发海啸、滑坡、崩塌、地裂缝等次生灾害。在我国,历史上也发生过很多破坏性很大的地震,以08年的四川大地震为甚,由此可见,地震灾害对人们生命财产安全造成了极大的威胁。
地震预警系统,是指在地震发生后,利用地震波传播速度小于电波传播速度的特点,提前对地震波尚未到达的地方进行预警。如果能够利用地震测站获取的地震信息,对地震震级进行快速评估,就可以在地震波到达之前发出预警。
地震震级是通过仪器给出地震大小的一种量度,我国目前普遍采用的方法是由仪器观测同时用地震波波形通过统计得到的经验公式来确定震级的大小,该方法有三个主要的缺点。
(1)震级饱和问题。从上述描述中可以看出,不同震级代表了地震时在不同频域中地震辐射波能的大小,ML、MB、MS这三种主要震级代表地震时约0.1-20s周期范围地震的能量大小,对于大于20s的地震波能量就反映不出了。因此确定上述震级对应的地震仪只能反映一定频带范围内的地震,而在此范围之外,由于仪器性能限制,无法测出地震发生时所辐射的周期大于20s的地震波(特别是大震级7.0面波震级以后),所确定的震级并不能反映实际地震的大小,即出现了“震级饱和”现象。
(2)物理意义不清晰。目前采用的震级估计方法仅由仪器指针摆动的幅值确定,没有明确的物理意义,无法深刻揭示地震震级与频度之间的关系。
(3)不具有实时性。目前采用的震级估计方法耗时较长,对于一个地震分析往往需要1-3分钟才能较为准确的得出其震级,无法实际应用于地震预警系统。
发明内容
本发明的目的就在于为了解决上述问题而提供一种用于克服中大型地震动中震级饱和问题,估计实时震级且应用于地震预警系统的方法。
为了达到上述目的,本发明采用了以下技术方案:
一种用于地震预警系统的实时地震震级估计方法,利用地震矩求出的矩震级,即:检测到信号后每间隔1s估算出相对应的矩震级,随着时间增长至20s左右得出与实际震级非常接近的估计震级,当地震在A点发生后,附近的B测站的操作包括以下步骤:
1)当第一个测站检测到地震动信号的第一秒时间内,采样频率为100Hz则有100个样点,进行滤波以及去零漂处理后可以得出该地震动第一秒信号的时间-加速度曲线,进行两次积分处理即可得出该地震动第一秒信号的时间-位移曲线,再进行快速傅里叶变换可以得出在频域下的频率-位移幅值(f-Ω)曲线,该曲线在对数坐标系下有如下关系:
(1)
式中,Ω0为地震谱零频时的振幅;fc为拐角频率;n和γ均为衰减常数;Q为独立频率因子;t为地震波传播时间。
在频率-位移幅值曲线对数坐标系下,对于小于拐角频率的频率值,其对应的位移幅值几乎趋于一个平台(见图3),该平台与y轴的交点即零频时的振幅值Ω0。
而零频时的振幅值Ω0与地震矩存在如下关系:
(2)
式中M0为地震矩(N·m);ρr为台站当地土密度(kg/m3);ρs为震源地土密度(kg/m3);r为震源距(km);cr为台站处剪切波速(m/s);cs为震源处剪切波速(m/s);Ω0为地震谱零频时的振幅;Rc(=0.63)为平均辐射常数;Fs(=2)为自由液面反射常数。
得出地震矩之后利用公式(3)估算出相应的矩震级Mw,马上以矩震级Mw形式向相关机构发出地震预警信号;
logM0=1.5Mw+9.05 (3)
式中Mw即为矩震级;M0为地震矩;
2)余下的n秒钟内,在每1秒重复上述步骤操作得出实时的地震震级,持续地发出地震预警信号。
本发明的有益效果在于:
本发明利用地震矩从而求出的矩震级能有效克服中大型地震(M>6)中出现的震级饱和问题,而且地震矩与地震动中零频时的振幅Ω0有明确联系,具有明确的物理意义。而且该方法具有实时性,检测到信号后每间隔1s均能估算出相对应的矩震级,随着时间增长至20s左右就能得出与实际震级非常接近的估计震级,能有效应用于地震预警系统。
附图说明
图1是本发明具体实施方法的操作流程示意图;
图2是地震预警系统示意图;
图3是5个测站对于某地震在双对数坐标系下的频率-位移幅值。
图4是2011.3.11东日本大地震FKS016测站的整个地震记录;
图5是图1滤波及去零漂后检测到P波信号后20s的地震加速度数据。
图6是图2中加速度经过两次积分得出的地震位移数据。
图7是图3中位移数据进行傅里叶变换在双对数坐标系下的频率-位移幅值。
图8是以2011.03.11东日本大地震M9.0一个测站为例利用零频时的振幅求出的矩震级及各参数列表。
具体实施方式
下面结合附图对本发明作进一步具体描述。
本发明是一种新型的用于克服中大型地震动中震级饱和问题,估计实时震级并应用于地震预警系统的谱反演方法,具体操作流程如图1所示。
实施例1:一种新型的用于地震预警系统的实时地震震级估计方法,以2011年东日本大地震为例,包括以下步骤:
步骤1:如图2,当地震在A点发生后,附近的B测站将连续接收到地震动信号,该信号为加速度信号。
步骤2:在B测站,每间隔一秒对该地震动信号进行分析,经过滤波和去零漂处理后作出时间-加速度曲线,如图5。
步骤3:对时间-加速度曲线进行两次积分处理即可得出该信号的时间-位移曲线,如图6。
步骤4:对时间-位移曲线进行快速傅里叶变换得出在频域下的频率-位移幅值(f-Ω)曲线,如图7。
步骤5:频率-位移曲线在双对数坐标系下有公式(1)的关系,利用公式(1)可求出零频时的振幅值Ω0。
步骤6:利用公式(2)可用零频时的振幅值Ω0求出相应的地震矩M0。
步骤7:得出地震矩后,利用公式(3)即可估算出相应的矩震级Mw,马上向C、D、E、F各个地方发出地震预警信号。
步骤8:每间隔1秒重复上述步骤操作,即可得出实时的地震震级,持续地发出地震预警信号。若利用多个测站的信号可以求出其平均值,能更加接近实际震级,如图8。
本发明的有益效果在于:
本发明利用地震矩从而求出的矩震级能有效克服中大型地震(M>6)中出现的震级饱和问题,而且地震矩与地震动中零频时的振幅Ω0有明确联系,具有明确的物理意义。而且该方法具有实时性,检测到信号后每间隔1s均能估算出相对应的矩震级,随着时间增长至20s左右就能得出与实际震级非常接近的估计震级,能有效应用于地震预警系统。
Claims (3)
1.一种用于地震预警系统的实时地震震级估计方法,利用地震矩求出的矩震级,即:检测到信号后每间隔1s估算出相对应的矩震级,随着时间增长至20s左右得出与实际震级非常接近的估计震级,当地震在A点发生后,附近的B测站的操作包括以下步骤:
1)当检测到地震动信号的第一秒时间内,进行滤波以及去零漂处理得出该地震动第一秒信号的时间-加速度曲线;进行两次积分处理即可得出该地震动第一秒信号的时间-位移曲线,再进行快速傅里叶变换可以得出在频域下的频率-位移幅值(f-Ω)曲线,该曲线在对数坐标系下有如下关系:
<mrow>
<mi>&Omega;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>&Omega;</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>f</mi>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&gamma;</mi>
<mi>n</mi>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mi>&gamma;</mi>
</mfrac>
</msup>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
<mi>Q</mi>
</mfrac>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Ω0为地震谱零频时的振幅;fc为拐角频率;n和γ均为衰减常数;Q为独立频率因子;t为地震波传播时间;
因零频时的振幅值Ω0与地震矩存在如下关系:
<mrow>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msubsup>
<mi>&pi;&rho;</mi>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msubsup>
<mi>&rho;</mi>
<mi>s</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msubsup>
<mi>rc</mi>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msubsup>
<mi>c</mi>
<mi>s</mi>
<mrow>
<mn>5</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msub>
<mi>&Omega;</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<msub>
<mi>F</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中M0为地震矩(N·m);ρr为台站当地土密度(kg/m3);ρs为震源地土密度(kg/m3);r为震源距(km);cr为台站处剪切波速(m/s);cs为震源处剪切波速(m/s);Ω0为地震谱零频时的振幅;Rc为平均辐射常数,取值为0.63;Fs为自由液面反射常数,取值为2;
得出地震矩之后利用公式(3)估算出相应的矩震级Mw,马上以矩震级Mw形式向相关机构发出地震预警信号;
logM0=1.5Mw+9.05 (3)
式中Mw即为矩震级;M0为地震矩;
2)余下的n秒钟内,在每1秒重复上述步骤操作得出实时的地震震级,持续地发出地震预警信号。
2.根据权利要求1所述的用于地震预警系统的实时地震震级估计方法,其特征在于,所述步骤2)中n的取值根据实际地震的震级升高而升高,n不超过20。
3.根据权利要求1所述的用于地震预警系统的实时地震震级估计方法,其特征在于,当采用多个测站的信号时,以多个测站的矩震级Mw的平均值发出地震预警信号,能更加接近实际震级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610049554.XA CN105676286B (zh) | 2016-01-25 | 2016-01-25 | 一种用于地震预警系统的实时地震震级估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610049554.XA CN105676286B (zh) | 2016-01-25 | 2016-01-25 | 一种用于地震预警系统的实时地震震级估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676286A CN105676286A (zh) | 2016-06-15 |
CN105676286B true CN105676286B (zh) | 2017-11-14 |
Family
ID=56302557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610049554.XA Expired - Fee Related CN105676286B (zh) | 2016-01-25 | 2016-01-25 | 一种用于地震预警系统的实时地震震级估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105676286B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106770656B (zh) * | 2016-12-21 | 2023-02-21 | 西南交通大学 | 一种用于模型试验中土体动剪切模量的测试装置 |
CN107798826B (zh) * | 2017-09-28 | 2020-04-24 | 北京交通大学 | 一种用于地震预警的震级预测方法 |
CN108037529B (zh) * | 2017-12-02 | 2019-10-11 | 桂林理工大学 | 一种基于初始振动信号的地震事件快速识别方法 |
CN109188524A (zh) * | 2018-10-26 | 2019-01-11 | 辽宁工程技术大学 | 考虑介质品质因子的地震中场地反应计算方法 |
CN110657935B (zh) * | 2019-09-24 | 2021-02-12 | 浙江大学 | 一种地震波加速度积分处理方法及系统 |
CN110727026B (zh) * | 2019-10-10 | 2021-05-28 | 太原理工大学 | 一种基于地震-海啸连续作用的非弹性位移比谱模型 |
CN110824560A (zh) * | 2019-11-11 | 2020-02-21 | 宜宾凯翼汽车有限公司 | 一种基于车联网的地震监测系统及方法 |
CN113126154A (zh) * | 2019-12-30 | 2021-07-16 | 李智 | 一种地震预报的面波频谱分析方法 |
CN112379423A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 微地震事件矩震级的确定方法及装置 |
CN113053064B (zh) * | 2021-03-15 | 2022-11-04 | 中国地震局地球物理研究所 | 一种面向区域性地震预警系统评测方法及评测系统 |
CN114355438B (zh) * | 2022-01-11 | 2022-08-09 | 西南交通大学 | 地震震级的估算方法、装置、设备及可读存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4385948B2 (ja) * | 2005-01-11 | 2009-12-16 | 大成建設株式会社 | リアルタイム地震情報を利用したリアルタイム地震応答波形推定方法 |
CN102628955A (zh) * | 2012-04-24 | 2012-08-08 | 成都高新减灾研究所 | 一种地震预警震级获取方法 |
CN102636806A (zh) * | 2012-05-10 | 2012-08-15 | 武汉大学 | 一种利用高频gnss数据确定地震震级的方法 |
JP5313706B2 (ja) * | 2009-01-27 | 2013-10-09 | 大成建設株式会社 | 地震動波形推定方法 |
CN104730577A (zh) * | 2015-03-04 | 2015-06-24 | 曾雄飞 | 一种基于共振包监测数据的地震预测方法 |
-
2016
- 2016-01-25 CN CN201610049554.XA patent/CN105676286B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4385948B2 (ja) * | 2005-01-11 | 2009-12-16 | 大成建設株式会社 | リアルタイム地震情報を利用したリアルタイム地震応答波形推定方法 |
JP5313706B2 (ja) * | 2009-01-27 | 2013-10-09 | 大成建設株式会社 | 地震動波形推定方法 |
CN102628955A (zh) * | 2012-04-24 | 2012-08-08 | 成都高新减灾研究所 | 一种地震预警震级获取方法 |
CN102636806A (zh) * | 2012-05-10 | 2012-08-15 | 武汉大学 | 一种利用高频gnss数据确定地震震级的方法 |
CN104730577A (zh) * | 2015-03-04 | 2015-06-24 | 曾雄飞 | 一种基于共振包监测数据的地震预测方法 |
Non-Patent Citations (2)
Title |
---|
地震的震级;陈运泰 等;《地震地磁观测与研究》;20041230;第25卷(第6期);第1-12页 * |
地震预警震级确定方法研究;金星 等;《地震学报》;20120930;第34卷(第5期);第593-610页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105676286A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676286B (zh) | 一种用于地震预警系统的实时地震震级估计方法 | |
Walsh et al. | Silver and Chan revisited | |
CN104007176B (zh) | 一种复杂岩土工程介质的全波场检测系统及方法 | |
CN102680575B (zh) | 一种复杂岩土介质的冲击映像方法及系统 | |
CN106052837B (zh) | 一种用于高速铁路地震预警中列车振动噪声识别方法 | |
Nanjo et al. | Seismicity prior to the 2016 Kumamoto earthquakes | |
Heidari | Quick estimation of the magnitude and epicentral distance using the P wave for earthquakes in Iran | |
CN102565855A (zh) | 油田压裂地面微地震数据处理方法 | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
Johnson et al. | Statistical comparison of methods for estimating sediment thickness from horizontal-to-vertical spectral ratio (HVSR) seismic methods: An example from Tylerville, Connecticut, USA | |
CN105204065A (zh) | 一种拾取初至波的方法和装置 | |
US20120186342A1 (en) | Method for Detecting Underground Tunnels | |
Zhan et al. | Accelerogram-based method for quick assessment of liquefaction occurrence | |
US20160245938A1 (en) | Method and device for increasing frequency of seismic digital signal | |
CN105954790B (zh) | 一种用于地震预警系统的快速地震震源距估计方法 | |
Kuo et al. | Identification of co-seismic ground motion due to fracturing and impact of the Tsaoling landslide, Taiwan | |
JP2022087776A (ja) | 地震波を時間領域に分解して地震前兆現象を捉える地震予測時間領域解析システム | |
CN112230275A (zh) | 地震波形的识别方法、装置及电子设备 | |
Özcebe et al. | On the comparison of 3D, 2D, and 1D numerical approaches to predict seismic site amplification: The case of Norcia basin during the M6. 5 2016 October 30 earthquake | |
Kong et al. | Evaluation of earthquake signal characteristics for early warning | |
KR101386255B1 (ko) | 지진측정장치 및 지진측정방법 | |
Bora et al. | Shear-Wave Splitting and Crustal Anisotropy in the Shillong–Mikir Plateau of Northeast India | |
CN115170428A (zh) | 一种声波远探测成像图的降噪方法 | |
Xu et al. | Application of MEMS data to fast inversion of rupture process: Tests with recordings from the IRREEW network | |
Zhao et al. | Comparison of Near‐Fault Velocity Pulse‐Like Ground Motions from the 2018 Mw 6.4 Hualien, Taiwan, Earthquake with the Next Generation Attenuation (NGA)‐West2 Ground‐Motion Models and Directivity Models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171114 Termination date: 20210125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |