地震预警系统中利用Pd估测PGV的改进方法

方案：

lgR≤0.86+0.17Mw（1）

依据以上条件确定的强地震动记录包括：2016年日本熊本7.3级地震主余震强地震动记录152组、2018年汶川8.0级地震主余震强地震动记录99组、2013年芦山7.0级地震主余震强地震动记录98组，1999年中国台湾集集7.6级地震主震强地震动记录128组，地震动记录信息如图1所示。

数据处理过程：对获得的加速度时程进行积分得到速度时程;采用Wu等（2005，2007）的方法，使用一个0.075 Hz的高通Butterworth滤波器对积分得到的速度、位移时程进行滤波;采用马强（2008）使用的STA/LTA方法对竖向强地震动记录位移时程进行P波的自动触发，并通过人工手段标记P波以获得Pd;采取矢量合成的方法对2个水平向速度时程计算，得到水平向的峰值速度。强地震动记录为非持续触发的记录，某些记录偶尔有丢头的情况，在这些记录之前增加0点，对于仍存在异常的数据，进行手动清除。

2 PGV-Pd关系统计分析

利用汶川地震余震、芦山地震主余震及日本熊本地震主余震事件的314组强地震动记录，统计Pd与PGV之间的关系（简称统计关系A）。由于高Pd记录较少、低Pd记录较多，在进行统计分析时需要提高高Pd记录的权重，因此将314组记录按Pd最大值和最小值平均分为5部分，将每个数据点的权重设置为该数据点所处部分数据数量的倒数，使每部分的权重相同、归一，并使用最小二乘法对带权重的PGV-Pd关系进行拟合。最终得到统计关系A，如图2所示。统计结果A反映出2个特征：①汶川地震余震、芦山地震主余震、熊本地震主余震记录均较均匀的分散于统计关系A两侧;②从统计结果的标准差来看，随着时间窗的增加，统计结果并未有明显的提升，因此从地震预警追求更早估测地震动参数的需求来说，采用3 s时窗的Pd观测值来估测PGV是比较合理的。

将统计关系A与汶川地震主震、集集地震主震强地震动记录进行对比，结果如图3所示，可以观察到2个明显的特征：①通过汶川地震主震、集集地震主震强地震动记录得到的散点大部分落在统计关系A的上方，即存在明显的PGV低估现象;②随着时间窗的增加，以统计关系A估测汶川地震主震、集集地震主震PGV的结果有明显提升的趋势，即PGV低估的现象减弱。

图2，3展示出来的特征表明，在原地地震预警系统中统计关系A能够较好地估测PGV，但当加入汶川地震、集集地震这样的大地震的记录时，对统计关系A的影响非常强烈，且存在着明显的低估效应。这可能有2方面的原因：①大地震的破裂过程复杂、破裂尺度大，且近场地震动特性与较小的地震存在差异，PGV-Pd关系在近场与小地震差异也会较大;②大地震能量释放较多、传播的距离较远，从图1b的数据分布来看，50 km范围内的4，5，6级地震分布相对均匀，而震中距大于50 km的记录多来自汶川地震和集集地震的主震。刘辰等（2019）研究发现，不同震中距条件下，PGV-Pd统计关系存在明显差异，将大地震、大震中距的记录用于传统的PGV-Pd统计不利于得到可信的统计关系。虽然依靠增加时間窗能够明显地降低加入大地震进行统计得到结果的差异（图2，3），但是与之相对应的有效预警时间将会明显地减少，对于地震预警系统分秒必争的需求来说并不是一个理想的选择。

鉴于以上原因，本文提出利用Pd估测PGV的改进方法：当Pd观测值较小时，以传统的PGV-Pd统计关系进行PGV估测;当Pd观测值大于一定阈值时给出高和一般2种PGV估测方法供用户和自动控制系统参考。

第一种方法，因为集集地震主震强地震动记录3 s时窗下Pd的最小值约为0.1 cm，所以设置分界值x0=0.1 cm，然后根据统计关系A确定分界点Q0（x0，y0），利用集集地震和汶川地震的主震记录进行过定点Q0的线性拟合，得到用于估测高PGV的统计关系B。一般PGV估测结果仍采用统计关系A得到，具体如图4所示。第二种方法，根据需求任意设计分界值，根据统计关系A得到斜率k，利用集集地震和汶川地震主震记录拟合平行于统计关系A的、用于估测高PGV结果的统计关系C。一般PGV估测结果仍采用统计关系A得到，如图5所示。统计关系A，B，C的详细参数见表1。

图4，5所示的方案具有3个特点：①由于采用过定点的拟合方式，所以统计关系B与A之间保证了连续性;②随着时间窗的增加，无论采用统计关系B还是C，2种PGV估测方法均趋于一致;③统计关系B的标准差比较高，且随着时窗的增加而减小，而统计关系C的标准差不仅较为稳定，且较小。

3 讨论

将本文方法与其他学者采用类似方法统计得到的PGV-Pd关系进行对比（图6），结果表明：本文方法和Peng等（2014）的研究结果非常接近，这一方面可能是因为汶川余震、熊本地震主余震和芦山地震主余震遵循着相似的统计关系，另一方面Peng等（2014）的研究中采用的296组记录中仅包含了4条主震记录，在建立统计关系时主震记录所占的权重比较小;Zollo等（2010）在研究中采用了意大利、日本和中国台湾省的数据，其中以意大利的数据为主，中国台湾省的数据所占比重较低，所得统计结果的斜率稍高于统计关系A;Wu和Kanamori（2005）在研究中仅采用中国台湾省的数据，统计关系明显高于其它研究的结果，因此我们推测地震预警系统PGV-Pd统计关系在不同的地区存在差异。

无论PGV-Pd统计关系是否存在地区差异，该关系对于大地震常常会低估PGV是一个不争的事实，在本文、Peng等（2014）以及Wu和Kanamori（2005）的研究中都存在着这样的现象。受到强地震动记录数量、质量的限制，现有的地震预警系统PGV-Pd关系统计研究中，从具有复杂破裂过程的汶川地震和集集地震中获得的强地震动记录数量较少，从而导致在计算中无法直接进行拟合，而如果单纯依靠增加这些强地震动记录的权重，则得到的统计结果可能会低估大地震的PGV或高估一般性地震的PGV。而本文则是通过大地震的少数记录，结合统计关系A得到统计关系B和C，当Pd达到一定阈值时给出高、一般2种PGV估测建议值，当Pd达到一定阈值时，建议参考这2种方法对PGV进行估测。统计关系B作为一种提供高PGV估测值的方法，和用于一般PGV值估测的统计关系A之间保持了连续性，但其斜率太高导致PGV随Pd的增加上升过快。统计关系C虽然和统计关系A之间未能保持连续，但统计关系C对数残差最接近正态分布，且与前人的研究较为一致，因此，并未出现PGV随Pd上升过快的现象。因此建议以统计关系C来提供高PGV的估测值。

为了便于比较本文设计的高PGV估测方法对数残差的频度分布，将3种方法（A，A&B，A&C）下汶川地震和集集地震主震强地震动记录的对数残差直方图示于图7。从图7a可以看出，统计关系A的PGV估测结果和观测值对数残差的优势分布偏离0值，即PGV估测值和观测值存在系统偏差，PGV存在低估现象。图7b，c表明统计关系B和C的对数残差结果优于统计关系A，其中优度最高的是统计关系C。

those by Relation A（a），Relation B（b），Relation C（c）

4 結论

本文筛选熊本地震主余震、汶川地震主余震、芦山地震主余震、集集地震主震强地震动记录，探讨了以Pd估测PGV的地震预警改进方法，提出了在Pd观测值高于一定阈值后提供高和一般2种PGV估测结果的方法，主要结论如下：

（1）利用我国汶川地震余震、芦山地震主余震、及日本熊本地震主余震事件的314组强地震动记录进行PGV-Pd关系的统计，得到以3～10 s为时间窗、以Pd连续估测一般PGV的统计关系A;利用汶川地震主震、集集地震主震强地震动记录，结合统计关系A，得到2种用于估测高PGV的统计关系B和C。

（2）通过和前人研究进行对比认为，统计方法B虽然在应用过程中和统计关系A保持了连续性，但由于斜率过高导致其和前人的研究结果偏差较大，且随着Pd的增加，PGV测值上升太快;而统计关系C与前人研究的一致性较好。

（3）分析了3种方法（A，A&B，A&C）的对数残差，发现统计关系A存在明显的优势分布偏离现象，再次证明了统计关系A在用于大地震时PGV的低估现象;统计关系B和C的对数残差较接近正态分布，统计关系C最接近正态分布。

（4）建议在地震预警系统中应用PGV-Pd方法时，当Pd达到一定阈值时，利用高、一般2种方法对PGV进行估测，即以统计关系A来提供一般PGV估测值，以统计关系C来提供高PGV的估测值。

由于早期地震信息的不足，实现准确的地震预警仍然是一件具有挑战性的工作，尤其是当大地震发生时，如何通过地震预警系统避免大地震复杂性带来的影响，降低漏报、低报的可能性，仍需要不断研究探讨。本文提出了一种新的思路，在地震发生时给出高、一般2种PGV估测值，但用户和自动判定系统采取怎样的手段或者如何从更多的途径去获得信息用于决定采用高PGV还是一般PGV的估测结果仍是一项需要研究的课题。

参考文献：

金星，张红才，李军，等.2012.地震预警震级确定方法研究[J].地震学报，34（5）：563-610.

刘辰，李小军，景冰冰，等.2019.基于地震动记录的距离分段地震预警PGV-Pd关系研究[J].地球物理学报，62（4）：1413-1426.

马强，李水龙，李山有，等.2014.不同地震动参数与地震烈度的相关性分析[J].地震工程与工程振动，34（4）：83-92.

马强.2008.地震预警技术研究及应用[D].哈尔滨：中国地震局工程力学研究所.

杨黎薇，林国良，邱志刚，等.2017.强震预警中P波到时STA/LTA和贝叶斯BIC双步骤捡拾研究[J].地震研究，40（4）：619-628.

杨黎薇，林国良，邱志刚，等.2018.基于人工神经元网络和多特征参数的预警震级估算[J].地震研究，41（2）：302-310.

Allen R M，Kanamori H.2003.The Potential for Earthquake Early Warning in Southern California[J].Science，300（5620）：786-789.

Nakamura Y.1988.On the Urgent Earthquake Detection and Alarm System（UrEDAS）[C].Tokyo-Kyoto：Proceedings of Ninth World Conference on Earthquake Engineering.

Peng C Y，Yang J S，Xue B，et al.2014.Exploring the Feasibility of Earthquake Early Warning Using Records of the 2008 Wenchuan Earthquake and Its Aftershocks[J].Soil Dynamics and Earthquake Engineering，57：86-93.

Wald D J，Quitoriano V，Heaton T H，et al.1999.Relationships between Peak Ground Acceleration，Peak Ground Velocity，and Modified Mercalli Intensity in California[J].Earthquake Spectra，15：557-564.

Wu Y M，Kanamori H，Allen R M，et al.2007.Determination of Earthquake Early Warning Parameters，τc and Pd，for Southern California[J].Geophysical Journal International，170：711-717.

Wu Y M，Kanamori H.2005.Rapid Assessment of Damaging Potential of Earthquakes in Taiwan from the Beginning of P Waves[J].Bullet Seismology Society of America，95（3）：1181-1185.

Wu Y M，Kanamori H.2008.Development of an Earthquake Early Warning System Using Real Time Strong Motion Signals[J].Sensors，（8）：1-9.

Zollo A，Amoroso O，Lancieri M，et al.2010.A Threshold-Based Earthquake Early Warning Using Dense Accelerometer Networks[J].Geophysical Journal International，183（2）：963-974.

Improved Method of Estimating PGV Using Pd for Earthquake Early Warning

LIU Chen1，3，LI Xiaojun2，3，JING Bingbing2，YIN Jianming4，GUO Xing5，ZHOU Yue3

（1.Postdoctoral workstation of China Reinsurance（Group）Corporation，Beijing 100033，China）

（2.The College of Architecture and Civil Engineering，Beijing University of Technology，Beijing 100124，China）

（3.Institute of Geophysics，China Earthquake Administration，Beijing 100081，China）

（4.China Property and Casualty Reinsurance Company LTD，Beijing 100033，China）

（5.Nuclear and Radiation Safety Center，Beijing 100082，China）

Abstract

Referring to PGV underestimation of great earthquake by PGV-Pd relation during earthquake early warning，we proposed an improved method of estimating PGV using Pd.Based on 314 strong ground motion records of the Wenchuan aftershocks，the Lushan earthquake and its aftershocks，the Kumamoto earthquake and its aftershocks，a 3s-10s time window of the relationship between displacement amplitude Pd and peak velocity PGV of earthquake is analyzed.Comparing the analysis results with 163 sets of strong ground motion records of the Chi-Chi earthquake and the Wenchuan earthquake，we presents a method for PGV estimation by Pd.When Pd is greater than a set threshold，two kinds of estimating results，high and general，are given for reference to users and automatic judging systems，and two schemes are designed to give the overestimation results.Through comparison among our relations and previous studies，as well as residual analysis of our relations，we come to the conclusion of a high and general estimation suggestion.

Keywords：Earthquake early warning system;Pd method;PGV estimation;Ground motion estimation