[发明专利]一种用于地震预警系统的实时地震震级估计方法

说明书

技术领域

本发明涉及一种新型的用于地震预警系统的地震震级的估计方法，尤其涉及一种用于克服中大型地震动中震级饱和问题，估计实时震级的谱反演方法。

背景技术

近年来，世界各地大

地震日益频繁，甚至会引发海啸、滑坡、崩塌、地裂缝等次生灾害。在我国，历史上也发生过很多破坏性很大的地震，以08年的四川大地震为甚，由此可见，地震灾害对人们生命财产安全造成了极大的威胁。

地震预警系统，是指在地震发生后，利用地震波传播速度小于电波传播速度的特点，提前对地震波尚未到达的地方进行预警。如果能够利用地震测站获取的地震信息，对地震震级进行快速评估，就可以在地震波到达之前发出预警。

地震震级是通过仪器给出地震大小的一种量度，我国目前普遍采用的方法是由仪器观测同时用地震波波形通过统计得到的经验公式来确定震级的大小，该方法有三个主要的缺点。

(1)震级饱和问题。从上述描述中可以看出，不同震级代表了地震时在不同频域中地震辐射波能的大小，M_L、M_B、M_S这三种主要震级代表地震时约0.1-20s周期范围地震的能量大小，对于大于20s的地震波能量就反映不出了。因此确定上述震级对应的地震仪只能反映一定频带范围内的地震，而在此范围之外，由于仪器性能限制，无法测出地震发生时所辐射的周期大于20s的地震波(特别是大震级7.0面波震级以后)，所确定的震级并不能反映实际地震的大小，即出现了“震级饱和”现象。

(2)物理意义不清晰。目前采用的震级估计方法仅由仪器指针摆动的幅值确定，没有明确的物理意义，无法深刻揭示地震震级与频度之间的关系。

(3)不具有实时性。目前采用的震级估计方法耗时较长，对于一个地震分析往往需要1-3分钟才能较为准确的得出其震级，无法实际应用于地震预警系统。

发明内容

本发明的目的就在于为了解决上述问题而提供一种用于克服中大型地震动中震级饱和问题，估计实时震级且应用于地震预警系统的方法。

为了达到上述目的，本发明采用了以下技术方案：

一种用于地震预警系统的实时地震震级估计方法，利用地震矩求出的矩震级，即：检测到信号后每间隔1s估算出相对应的矩震级，随着时间增长至20s左右得出与实际震级非常接近的估计震级，当地震在A点发生后，附近的B测站的操作包括以下步骤：

1)当第一个测站检测到地震动信号的第一秒时间内，采样频率为100Hz则有100个样点，进行滤波以及去零漂处理后可以得出该地震动第一秒信号的时间-加速度曲线，进行两次积分处理即可得出该地震动第一秒信号的时间-位移曲线，再进行快速傅里叶变换可以得出在频域下的频率-位移幅值(f-Ω)曲线，该曲线在对数坐标系下有如下关系：

(1)

式中，Ω₀为地震谱零频时的振幅；f_c为拐角频率；n和γ均为衰减常数；Q为独立频率因子；t为地震波传播时间。

在频率-位移幅值曲线对数坐标系下，对于小于拐角频率的频率值，其对应的位移幅值几乎趋于一个平台(见图3)，该平台与y轴的交点即零频时的振幅值Ω₀。

而零频时的振幅值Ω₀与地震矩存在如下关系：

(2)

式中M₀为地震矩(N·m)；ρ_r为台站当地土密度(kg/m³)；ρ_s为震源地土密度(kg/m³)；r为震源距(km)；c_r为台站处剪切波速(m/s)；c_s为震源处剪切波速(m/s)；Ω₀为地震谱零频时的振幅；R_c(＝0.63)为平均辐射常数；F_s(＝2)为自由液面反射常数。

得出地震矩之后利用公式(3)估算出相应的矩震级M_w，马上以矩震级M_w形式向相关机构发出地震预警信号；

logM₀＝1.5M_w+9.05(3)

式中M_w即为矩震级；M₀为地震矩；

2)余下的n秒钟内，在每1秒重复上述步骤操作得出实时的地震震级，持续地发出地震预警信号。

本发明的有益效果在于：