Real-time estimation of rupture length of earthquake source based on records of strong motion array
-
摘要: 本文提出一种通过小孔径台阵波束形成技术快速估测震源破裂长度的方法，同时提出了调整各台站相位的方法，即在每个台站的延时无法确定的情况下，近似地得到每个台站之间相位差的方法。在此基础上，利用自贡强震动观测台阵在汶川地震中得到的主震地震动记录进行了破裂长度的离线快速估测。结果表明，该方法的时效性较好，在地震动加速度峰值到来前约10 s时，便得到了断层破裂长度（约120 km）和相应破裂走向的估测结果，虽然此时的估测结果与最终估测结果（约230 km）有较大的偏差，但从实现快速地震预警的角度来看，还是有较好的实际意义。通过与前人的研究结果进行对比，结果显示该方法低估了东北方向的破裂长度，于是探讨了该现象产生的原因及可能的解决方案。最后对这种算法的计算速度进行了测试分析，结果表明通过并行计算该方法能够满足实时断层破裂长度估测的要求。Abstract: We develop a fast rupture length estimation method by using small-aperture seismic arrays beam-forming technology, and proposed a method to adjust the phase of each station. By the method the phase difference of each station can be calculated approximately, when delay of each station cannot be determined. The result of Wenchuan earthquake outline test shows that the time efficiency of this method is good, for example, at about 10 s before the arrival of peak ground acceleration, the estimation result of rupture length is about 120 km. Although this estimation has a large deviation compared with the final estimation results （about 230 km）; however, from the point of view of fast earthquake warning, it still has good practical significance. By comparing with the previous research results, it is indicated that this method underestimates the rupture length in northeast direction, and the causes and possible solutions are discussed. Finally, the computational speed of the algorithm is studied, the results show that the method can meet the requirement of real-time fault rupture length estimation by parallel computing.
-
引言
随着我国地震观测台站数量的不断增加，区域地震监测能力的逐渐增强，可监测震级的下限不断降低。密集地震观测台网产出高质量地震观测数据的同时，也对现有业务系统提出了更多的新需求，其中从实时产出的地震观测数据中准确判别和拾取不同类型的震相是最基础的、最迫切的需求之一。这些震相信息为更好地认识地壳介质结构、了解发震断层、获取震源信息、进行地震预警服务等提供数据支持。
为了更快更准确地从连续地震记录波形中判别和拾取出地震震相，许多学者都开展了系列研究，并提出了很多实用的算法。这些震相拾取算法可归纳为能量准则（Allen，1982；Baer，Kradolfer，1987；Lomax et al，2012；马强等，2013）、极化分析（Vidale，1986；Ruud，Husebye，1992；Amoroso et al，2012）、模糊分析（Chu，Mendel，1994）、人工神经网络（Dai，MacBeth，1995）、高阶统计量（Saragiotis et al，2002；Küperkoch et al，2012）、小波变换（Anant，Dowla，1997；刘希强等，1998，2000）及基于机器学习（Kong et al，2019；Perol et al，2018；Ross et al，2018a，b；Zhu et al，2018；李安等，2020）等几大类相关算法。
相较于初动清晰、信噪比较高的P波震相，S波震相由于受到P波尾波和（或）其它转换波震相的干扰，信噪比通常较P波要低得多，因此S波震相拾取结果的准确性普遍低于P波震相。即使是有经验的分析人员，仍无法保证S波震相拾取的准确性。现有的S波震相自动拾取算法大都是基于极化特征分析的方法（Flinn，1965；Vidale，1986；Jurkevics，1988；Cichowicz，1993），这些算法利用三分向地震波形记录，基于P波震相与S波震相的不同偏振特征，通过计算质点运动的偏振度、线性度等参数，寻找特征函数突变发生的点，进而判定S波震相到时的位置。此类方法大都适用于单台记录，算法有效性依赖于其有效剥离多种震相的能力，如将剪切波与面波分离等。Lois等（2013）即利用一个基于特征值分析的特征函数，提出了一种S波自动拾取时域算法，该算法对于P波、S波的极化方向不作任何假设，是一种近乎“零参数”的时域实时算法。本文则在该算法的基础上，进一步优化相关参数，并利用福建地震台网2015—2018年期间的9 855条三分向波形记录分析进行验证。
1. 方法
在实时地震数据处理中，通常采用滑动窗的方式将连续观测数据离散为固定长度的数据，而后再采用相应离线数据处理算法进行分析。同样，对于一段已知P波到时位置的波形记录，本文在P波到时后采用滑动窗将记录划分，并根据一定的滑动步长进行滑动处理，然后对各小窗的三分向波形记录求取特征值，进而可获得特征值时程，本文将该特征值时程作为S波识别特征函数。其中，由三分向波形求取特征值可采用数学分析中协方差矩阵分解的方式得到，特征值时程的物理含义是确定质点运动偏振方向。假定由三分向波形数据计算得到的协方差矩阵为
$${\boldsymbol{C}} {\text{＝}} \left[ {\begin{array}{*{20}{c}} {{\rm{cov}} {\text{（}}x{\text{，}}x{\text{）}}}&{{\rm{cov}} {\text{（}}x{\text{，}}y{\text{）}}}&{{\rm{cov}} {\text{（}}x{\text{，}}z{\text{）}}} \\ {{\rm{cov}} {\text{（}}y{\text{，}}x{\text{）}}}&{{\rm{cov}} {\text{（}}y{\text{，}}y{\text{）}}}&{{\rm{cov}} {\text{（}}y{\text{，}}z{\text{）}}} \\ {{\rm{cov}} {\text{（}}z{\text{，}}x{\text{）}}}&{{\rm{cov}} {\text{（}}z{\text{，}}y{\text{）}}}&{{\rm{cov}} {\text{（}}z{\text{，}}z{\text{）}}} \end{array}} \right]{\text{，}}$$ (1) 式中，
${\rm{cov}} {\text{（}}x{\text{，}}y{\text{）}} = \Big[\sum\limits_{i {\text{＝}} 1}^L {{\text{（}}{x_i} {\text{－}} \overline x{\text{）}}{\text{（}}{y_i} {\text{－}} \overline y{\text{）}}}\Big]\Big/ L$ ，L为滑动窗长度，且${\rm{cov}} {\text{（}}x{\text{，}}y{\text{）}} {\text{＝}} {\rm{cov}} {\text{（}}y{\text{，}}x{\text{）}}$ ，${\rm{cov}} {\text{（}}x{\text{，}}z{\text{）}} {\text{＝}} $ $ {\rm{cov}} {\text{（}}z{\text{，}}x{\text{）}}$ ，${\rm{cov}} {\text{（}}y{\text{，}}z{\text{）}} {\text{＝}} {\rm{cov}} {\text{（}}z{\text{，}}y{\text{）}}$ 。将式（1）分解即可分别得到三个特征值${\lambda _1}{\text{，}}{\lambda _2}{\text{，}}{\lambda _3}$ （${\lambda _1} {\text{＞}} {\lambda _2} {\text{＞}} {\lambda _3}$ ）及相应的三个特征向量${ {\boldsymbol{u}}_{{1}}}{\text{，}}{ {\boldsymbol{u}}_{{2}}}{\text{，}}{ {\boldsymbol{u}}_{{3}}}$ 。对于连续波形记录，即可得到特征值时程${\lambda _1}{\text{（}}t{\text{）}}{\text{，}}{\lambda _2}{\text{（}}t{\text{）}}{\text{，}} $ $ {\lambda _3}{\text{（}}t{\text{）}}$ 。据此也可分别计算得到信号偏振线性度、质点运动平面度、P波传播方位角、线性运动入射角等参数（Jurkevics，1988；马强，2008）。由于
${\lambda _1}{\text{（}}t{\text{）}}$ 对信号传播方向上的能量突变更为敏感（图1d），因此本文仅选择其作为特征函数。为了更好地体现信号中能量的微弱变化并突显信号突变位置，本文求取${\lambda _1}{\text{（}}t{\text{）}}$ 的平方根$f{\text{（}}t{\text{）}} {\text{＝}} \sqrt {{\lambda _1}{\text{（}}t{\text{）}}}$ （图1e），放大${\lambda _1}{\text{（}}t{\text{）}}$ 中低值成分的同时压低${\lambda _1}{\text{（}}t{\text{）}}$ 中高值成分进而减小其动态范围，实现对特征函数的压缩。尤其对于信噪比较低的事件波形，采用此方式的特征函数凸显信号中的能量突变位置更为有效。随后，选取P波初至到S尾波段对$f{\text{（}}t{\text{）}}$ 求取峰度系数$K{\text{（}}t{\text{）}} {\text{＝}} $ $ kur[f{\text{（}}t{\text{）}}]$ ，即可从特征函数中将S波到时位置进一步凸显出来（图1f）。应用该时间窗是为了避免将后续事件的初至震相错误识别为本次事件S波震相，但对于密集事件序列，该时间窗的选取比较困难。由图1f可见，在S波到时附近$K{\text{（}}t{\text{）}}$ 呈现一个显著的陡坎，但S波到时点位于陡坎起始最低处，因此在实际应用中很难准确拾取，需要进一步处理，即对$K{\text{（}}t{\text{）}}$ 求取一阶差分$\Delta K{\text{（}}t{\text{）}}$ ，此时S波到时与$\Delta K{\text{（}}t{\text{）}}$ 峰值位置相一致，S波到时位置将更容易判别，据此即可对S波粗到时进行拾取。最后，应用自回归赤池信息准则（Akaike information criteria，缩写为AIC）（Akaike，1974）在所拾取S波粗到时位置附近精确判别（前后各取0.3 s），即可获得准确可靠的S波到时信息。![]()
图 1 本文S波拾取方法算例示意图事件发震时刻为2015-01-02 00:33:56.89，震中距为3.8 km的GTSK台站记录，窗长为0.2 s。图（a−c）为三分向波形记录；图（d）为特征值时程；图（e）为特征函数f （t）；图（f）为峰度系数时程K（t）；图（g）为求取差分后特征函数ΔK（t）Figure 1. An example of S phase picking by using the algorithm of this studyThe origin time of the event is 00:33:56.89 on 2 January 2015. The wavforms were recorded by the station GTSK with epicentral distance 3.8 km，and window length is taken as 0.2 s. Figs. （a） to （c） are three-component seismic records；Fig. （d） shows the history of three eigenvalues；Fig. （e） is the characteristic function f （t）；Fig. （f） shows time history of Kurtosis coefficient K（t），and Fig. （g） shows the differential characteristic function ΔK（t）滑动时窗长度是本文算法中需要确定的重要参数。若选择较短的时窗，该算法对于特征函数中微小的改变过于敏感，拾取到的S波到时位置将超前；若选择过长的时窗，算法将缺乏必要的敏感度，拾取到的S波到时位置将显著滞后。实时地震数据处理中，算法的鲁棒性是一个重要的指标，选择某一固定时间窗长度显然无法满足所有事件的需求。为此，本文采用多窗口综合加权方法，窗口长度为0.2—1.4 s，间隔0.2 s，共7个。首先在各窗口中拾取S波到时Si，并计算S震相信噪比Ri，然后依据信噪比综合加权给出S波到时位置，即
$${S_{\rm{f}}} {\text{＝}} \dfrac{{\displaystyle\sum\limits_{i {\text{＝}} 1}^7 {{S_i}{R_i}} }}{{\displaystyle\sum\limits_{i {\text{＝}} 1}^7 {{R_i}} }}{\text{.}}$$ (2) 2. 数据
本文研究中收集了福建台网2015—2018年记录的1 862次地方震事件，震级范围为ML−0.5—4.2，仅挑选其中震中距小于100 km的9 855条三分向地震动记录计算分析，在此震中距范围内，初至S波震相均为Sg。本文以人工地震编目提供的Sg震相拾取结果为参考，对比采用本文方法的拾取震相到时拾取位置，以验证本文算法的准确性及实用性。本研究中使用的台站观测记录以仙游震群和台湾海峡南部震群为主，这两个震群的观测数据量约占所有分析记录1/3左右，各台站记录数量空间分布如图2所示。本文所用台站记录在震中距100 km内基本呈均匀分布，震中距分布如图3所示。
![]()
图 2 本文所用台站记录随震中距和震级分布Figure 2. Epicentral distance and magnitude distribution for earthquake records used in this study3. 结果与讨论
应用本文方法得到的S波到时拾取结果与人工编目S震相拾取结果间的偏差统计情况如图4所示。结果显示，相较于人工拾取结果，本文方法的S波震相到时拾取平均偏差为（−0.003±1.34） s，其中拾取偏差小于0.2 s的记录所占比例为63.9%，拾取偏差小于0.5 s的记录所占比例为79.6%，即应用本文方法能够较准确地对大部分S波震相到时进行拾取。但从图3中也可以明显看出，应用本文方法的S波震相到时拾取结果有超前于人工拾取结果的趋势，因而造成图4b中统计直方图的偏态分布。本文将测试数据进一步划分为仙游震群、台湾海峡震群及其它事件三类，每类分别包括1 276条、2 330条和6 249条记录，并分析不同类别的S波到时拾取结果精度，结果如图5所示。由该图可见，不同的类别事件应用本文方法的S波拾取结果偏差分布较为一致，均表现出超前于人工拾取结果的趋势，三类事件的S波拾取偏差分别为（−0.018±1.68） s，（−0.053±1.56） s和（0.011±1.17） s。综上可见，本文所研究的相关算法具有较好的适用性，应用于不同区域的事件时均表现稳定。
![]()
图 4 应用本文方法的S波震相到时拾取偏差统计图（a）为各记录S震相拾取偏差散点图，图（b）和（c）为S震相拾取偏差统计直方图Figure 4. Statistics of S phase arrival time pick errors by using our methodFig. （a） gives S phase picking error for each record，Figs. （b） and （c） show the histogram of the phase picking error![]()
图 5 分区域S波震相到时拾取偏差结果统计（a） 仙游震群；（b） 台湾海峡南部震群；（c） 其它事件Figure 5. Statistics of S phase arrival time picking error in different regions（a） Xianyou sequence；（b） Southern Taiwan Strait sequence；（c） Other events此外，由图4还可以看出，仍有约11.2%的记录的S波震相到时拾取偏差超过±1.0 s，其中约4.1%的记录的S波震相到时拾取偏差超过±2.0 s，图6分别列举出三个拾取偏差较大的记录情况。进一步分析后，可将引起拾取偏差过大的原因归纳为三个方面。首先，信噪比较低是造成S波震相拾取偏差较大的显著原因。如图6a中所示波形记录，由于S波震相受到前序P波尾波的污染，信噪比较低，即使对于有经验的分析人员仍较难准确拾取，这也是通用S波震相到时拾取算法准确度不高的主要困难所在。本文采用式（3）截取人工拾取S波到时前后各2 s窗长内的波形记录并计算均方根进而得到各记录的信噪比，对于采样率为100 的记录：
![]()
图 6 S波震相到时拾取偏差较大的波形记录（a） 低信噪比记录结果；（b） 图（a）中记录1—20 Hz带通滤波后的结果；（c） 多事件叠加记录结果；（d） EW向异常记录结果Figure 6. Waveforms with large S phase arrival time picking error（a） The result for a low SNR record；（b） The result for the record in Fig. （a） after 1−20 Hz band-pass filtering；（c） The result for a multi-event record；（d） The result for an abnormal record in EW component$${\rm{SNR}} {\text{＝}} {\dfrac{{\sqrt {\dfrac{1}{n}\displaystyle\sum\limits_{i {\text{＝}} s}^{s {\text{＋}} n} {{x_i}} } } }{ {\sqrt {\dfrac{1}{n}\displaystyle\sum\limits_{i {\text{＝}} s {\text{－}} n}^s {{x_i}} } }}}{\text{，}}$$ (3) 式中s为人工拾取S波到时位置，n为数据长度（取200）.
随后，对406条S波到时拾取偏差超过±2.0 s的记录的信噪比进行统计，并与其它记录的信噪比情况进行对比，如图7所示。可以直观地看出低信噪比（＜2.0）记录在拾取偏差超过±2.0 s的记录中占有最大比例（图7a），与之相应，对于其余拾取误差较小的记录，高信噪比记录所占的比例较高（图7b）。图7以信噪比2.0为界，分别给出了低信噪比和高信噪比情况下S波震相拾取误差分布情况。由该结果可见：对于低信噪比波形记录，S波震相拾取误差超过1.0 s的比例为19%，显著多于高信噪比波形记录中的比例（9%）；且高信噪比波形记录中S波震相拾取误差小于0.2 s的比例（70%）更明显高于低信噪比波形记录中的比例（52%）。由此，作者认为观测记录信噪比是显著影响S波震相拾取精度的最重要因素。表1分别列出了针对不同信噪比记录，应用本文方法统计得到的S波震相拾取精度结果。可见，随着记录信噪比的提高，S波拾取残差呈平均值减小、标准差降低的趋势，与图7、图8所示结果一致。此外，从应用本文方法拾取的S－P到时差与人工拾取S－P到时差的对比（图9a）中也可以看出，偏差较大（如＞±2.0 s）记录的信噪比大多较低，而且随着震中距增加（S－P到时差增大），S波拾取结果的偏差分布也随之更为离散，显然这也与记录信噪比逐渐降低有关。同时，本文也对比了不同震群事件的S波拾取精度（图9b-d），本文算法在不同震群中的表现也较为一致，其中：仙游震群（图9b）中震中距较小的记录占比较高，算法对于此部分记录（＜5.0 s）也显示了更佳的性能，显然其信噪比也更高；同理海峡震群（图9c）和其它事件（图9d）中大震中距记录占比也相对更高，S波拾取残差随信噪比降低也变得更为离散，再次说明了记录的信噪比在应用本文S波拾取算法中的重要性。
![]()
图 7 地震记录信噪比的分布（a） 拾取偏差大于2.0 s的406条记录；（b） 拾取偏差小于2.0 s的9449条记录Figure 7. SNR distribution of seismic records（a） 406 records with S phase picking errors larger than ±2.0 s；（b） 9449 records with S phase picking errors less than ±2.0 s表 1 不同信噪比记录的S波到时拾取偏差统计Table 1. Statistic on S phase picking error for different SNR recordsSNR 记录数量 偏差均值/s 偏差中值/s 偏差标准差/s ＜2 2 953 0.135 0.040 1.92 2—5 6 416 −0.051 0.050 0.99 ≥5 486 −0.202 0.040 0.96 ![]()
图 8 不同信噪比记录S波震相拾取偏差分布（a） 信噪比小于2.0的记录；（b） 信噪比大于2.0的记录Figure 8. Pi chart of S phase detection error for different SNR records（a） Records with SNR less than 2.0；（b） Records with SNR larger than 2.0![]()
图 9 应用本文方法所得的S－P到时差与人工拾取结果对比及其信噪比分布（a） 所有记录；（b） 仙游震群；（c） 台湾海峡震群；（d） 其它事件Figure 9. Comparison of the arrival time differences for S−P by our method with those by manual picking and SNR distributions（a） All records；（b） Xianyou sequence；（c） Taiwan Strait sequence；（d） Other events若对事件波形进行带通滤波（如1—20 Hz），则可在一定程度上提高记录质量，也将有利于提高S波震相的拾取精度。图6b展示了图6a所示波形记录经1—20 Hz零相移带通滤波后再次应用本文方法的S波震相拾取结果。图10则分别为该记录滤波前后应用本文方法的具体识别效果，其中各子图所展示内容与图1一致。如前文已提及，不同滑动窗长下，S波到时拾取位置也会有所差异，因而本文采用式（2）对不同窗长下的S波到到时拾取位置进行加权综合，而图1和图10仅为某一窗长（分别为0.2 s和1.0 s）下的拾取结果，与图6中所示S波到时综合加权判别结果（虚线）会有所差异。对比可见，由于波形记录质量显著提高，S波震相拾取也更为准确，相较于人工拾取结果的偏差仅为0.36 s。需要说明的是，本文图4中所示震相拾取偏差均基于原始波形，未经带通滤波。鉴于上述分析，可推测出，若能够通过带通滤波提高部分低信噪比记录的信噪比水平，则震相拾取的偏差将得以改善。作者也将收集相关资料，继续深化该部分的研究。
![]()
图 10 图6a中信噪比较低记录带通滤波前（a）、后（b）的S波到时拾取各子图意思同图1，事件发震时间为2015-01-08 05:48:44.56，台站ZPCH的震中距为99.8 km，窗长为1.0 sFigure 10. S phase picking for a low SNR record in Fig. 6a before （a） and after （b） band-pass filteringEach subfigure has the same meaning as Fig. 1. The event occurred at 05:48:44.56 on 8 January 2015，which was recorded by the station ZPCH with epicentral distance 99.8 km. And time window length is taken as 1.0 s针对地震序列（尤其是大震后的密集余震序列）的处理能力也是考验震相拾取算法实用性的一个重要指标。鉴于序列震事件记录的复杂性及多样性，现有常规处理算法（如STA/LTA，AR-AIC，Z-detect等）均不能较为有效地处理和应对，本研究所涉及相关算法本质上也属于这类算法范畴。如前文所述，该算法不仅依赖于P波到时位置，而且仅能在P波到时之后的一段时间窗内识别唯一一个S波到时位置。显然，该算法基本不具备处理更为复杂的情况下序列地震的能力。图6c展示了一个密集序列的波形记录，在30 s内记录了多个小震事件，由于记录过于复杂，人工分析也仅标注了其中主要的S波震相到时（图6c中所标识Sg震相），而本文方法则只拾取到了序列中的一个可能的S波震相。由于本文方法尚不具备连续拾取后续震相的功能，因此本文方法在此类序列地震记录中的处理结果欠佳。在实际业务系统应用中，可适当提高S波尾波段阈值（即缩短事件持续时间），部分程度上可弥补在序列事件处理中的不足。
此外，由于波形记录本身存在的质量缺陷也是造成S波震相拾取偏差过大的显著影响因素。本文中所讨论相关算法的有效性和可靠性均依赖特征值分析结果，而一旦所分析的波形存在异常，则将直接影响特征值的分解结果，进而造成S波到时拾取结果的较大偏差。由图6d可见，由于EW向分量波形显著异常（该分向无输出），特征值计算错误，最终导致S波拾取结果出现较大偏差。尽管此类问题在实际业务应用中无法完全避免，但通过对记录波形质量的分析监控，能够有效地避免将异常波形纳入处理，从而提高S波到时拾取结果的准确性。
4. 结论
本文基于特征值分解算法，研究了一种可用于地方震S波震相实时提取的实用化算法。该算法对P、S波偏振方向不作任何假设，通过应用多个滑动时间窗即可较准确地拾取出S波震相到时，避免由于窗口长度选择不合理而造成的拾取偏差，满足地震信号的实时处理需求。应用福建地震台网地震观测的数据测试结果表明，该算法S波震相拾取精度高，适用于实时地震信息处理系统，但尚不具备密集序列连续处理能力。此外，由于应用该算法前需要首先已知P波到时信息，因此S波震相到时拾取的可靠性一定程度上也依赖于P波震相到时拾取的准确程度。但相较于S波震相，P波震相初动更为明显，因而更易于拾取，拾取精度更有保障。本文仅以福建地震台网100 km范围内的地方震事件为主，所拾取震相也主要为Sg震相，重点对该算法的适用性展开分析讨论，因而所得相关结论也可能不够全面。作者将继续收集相关波形数据，以更充分地论证该算法的适用性及其在实际数据处理系统中应用的可靠性。
随着地震观测台站数量的不断增加及观测数据质量的逐步提升，近年来深度学习等人工智能算法在地震数据处理领域内的应用也越来越多。尤其是在震相拾取方面，已经有多个学者提出了高性能、高可信的相关震相拾取算法，如GPD算法（Ross et al，2018a）、U-net算法（赵明等，2019）、PhaseNet算法（Zhu，Beroza，2018）等。作者也将应用这类算法开展相关研究，并将其与传统方法进行整合，以期更快、更准、更全地产出震相到时信息。
审稿专家为本文提出了建设性的修改意见，作者在此表示感谢。
-
![]()
图 2 自贡强震台阵记录到的汶川地震三分量加速度时程
Figure 2. The three-component acceleration time histories of Wenchuan earthquake recorded by Zigong strong motion array
![]()
图 3 强震动记录的筛选过程示意图
（a） 1号台站前60 s三分量强震动记录及其S波标志震相；（b） 本文采用的汶川地区地下速度结构；（c） 自贡台阵记录到的汶川地震各台站加速度时程之间的相关性
Figure 3. Diagram of the screening process for strong ground motion records
（a） The pre 60 s three-component strong motion records of the station No.1 and the seismic phase of the S wave；（b） Underground velocity structure of seismic wave in Wenchuan area；（c） Correlation between acceleration time histories of each station recorded by Zigong array during Wenchuan earthquake
![]()
图 4 经过震相对齐、滤波和相位调整后的地震动加速度时程
Figure 4. The time histories of ground motion after seismic phase alignment，filtering and station delay processing
![]()
图 5 离线条件下汶川地震瞬时主能量源实时估测结果（a）及同时刻2号参考台站的垂向加速度时程（b）
图（a）中圆点的位置代表这一时刻的主能量源位置，圆点的大小代表这一时刻相关系数和的极值，虚线表示两阶段估测的破裂长度范围
Figure 5. The outline result of the real-time estimation of the main energy source of Wenchuan earthquake （a） and the vertical acceleration record of the reference station No.2 at the same time （b）
In Fig.（a），the dot position represents the location of main energy source at the moment，the dot size represents the maximum size of correlation coefficient at the moment，and the dotted lines represent range of the rupture
-
阿布都瓦里斯•阿布都瓦衣提. 2013. 基于台阵观测资料的地震动空间变化特性研究[D]. 北京: 中国地震局地球物理研究所: 19–23. Abduwahit A. 2013. Study on Spatial Variation of Ground Motions Based on Array Records[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 19–23 （in Chinese）.
金星, 张红才, 李军, 韦永祥, 马强. 2012. 地震预警连续定位方法研究[J]. 地球物理学报, 55（3）: 925-936. Jin X, Zhang H C, Li J, Wei Y X, Ma Q. 2012. Research on continuous location method used in earthquake early warning system[J]. Chinese Journal of Geophysics, 55（3）: 925-936 （in Chinese）
李小军. 2001. 对近年大震震害现象与工程地震问题研究的思考[J]. 国际地震动态, （8）: 26-32. Li X J. 2001. Comments on phenomena of great earthquake damages and engineering seismological researches[J]. Recent Developments in World Seismology, （8）: 26-32 （in Chinese）.
刘瑞丰, 陈翔, 沈道康, 郑秀芬, 杨辉, 陈宏峰. 2014. 宽频带数字地震记录震相分析[M]. 北京: 地震出版社: 100–113. Liu R F, Chen X, Shen D K, Zheng X F, Yang H, Chen H F. 2014. Phase Analysis on the Broad Frequency Band Digital Seismic Recording[M]. Beijing: Seismological Press: 100–113 （in Chinese）.
马强. 2008. 地震预警技术研究及应用[D]. 哈尔滨: 中国地震局工程力学研究所: 34–40. Ma Q. 2008. Study and Application on Earthquake Early Warning[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 34–40 （in Chinese）.
唐晖, 李小军, 李亚琦. 2012. 自贡西山公园山脊地形的场地效应分析[J]. 振动与冲击, 31（8）: 74-79. Tang H, Li X J, Li Y Q. 2012. Site effect of topography on ground motions of Xishan park of Zigong city[J]. Journal of Vibration and Shock, 31（8）: 74-79 （in Chinese）.
吴建平, 黄媛, 张天中, 明跃红, 房立华. 2009. 汶川MS8.0级地震余震分布及周边区域P波三维速度结构研究[J]. 地球物理学报, 52（2）: 320-328. Wu J P, Huang Y, Zhang T Z, Ming Y H, Fang L H. 2009. Aftershock distribution of the MS8.0 Wenchuan earthquake and three dimensional P-wave velocity structure in and around source region[J]. Chinese Journal of Geophysics, 52（2）: 320–328 （in Chinese）.
杨宇, 李小军, 贺秋梅. 2011. 自贡西山公园山脊场地地形和土层效应数值模拟[J]. 震灾防御技术, 6（4）: 436-447. Yang Y, Li X J, He Q M. 2011. Numerical simulation for site effect of ridge terrain and overlaying soil in Zigong Xishan park[J]. Technology for Earthquake Disaster Prevention, 6（4）: 436-447 （in Chinese）.
张小飞, 陈华伟, 仇小峰, 汪飞, 王成华, 张弓, 徐大专. 2015. 阵列信号处理及MATLAB实现[M]. 北京: 电子工业出版社: 50–55. Zhang X F, Chen H W, Qiu X F, Wang F, Wang C H, Zhang G, Xu D Z. 2015. Array Signal Processing and MATLAB Implemen-tation[M]. Beijing: Publishing House of Electronics Industry: 50–55 （in Chinese）.
张勇, 冯万鹏, 许力生, 周成虎, 陈运泰. 2008. 2008年汶川大地震的时空破裂过程[J]. 中国科学: D辑, 38（10）: 1186-1194. Zhang Y, Feng W P, Xu L S, Zhou C H, Chen Y T. 2009. Spatio-temporal rupture process of the 2008 great Wenchuan earthquake[J]. Science in China Series D: Earth Sciences, 2009, 52（2）: 145-154.
Akaike H. 1973. Information theory and an extension of the maximum likelihood principle[G]//Selected Papers of Hirotugu Akaike. New York: Springer: 199–213.
Allen R M, Kanamori H. 2003. The potential for earthquake early warning in southern California[J]. Science, 300（5620）: 786-789.
Allen R M, Brown H, Hellweg M, Khainovski O, Lombard P, Neuhauser D. 2009. Real-time earthquake detection and hazard assessment by ElarmS across California[J]. Geophys Res Lett, 36（5）: L00B08.
Allen R V. 1978. Automatic earthquake recognition and timing from single traces[J]. Bull Seismol Soc Am, 68（5）: 1521-1532.
Allen R V. 1982. Automatic phase pickers: Their present use and future prospects[J]. Bull Seismol Soc Am, 72（6B）: S225-S242.
Borcea L, Papanicolaou G, Tsogka C, Berryman J. 2002. Imaging and time reversal in random media[J]. Inverse Probl, 2002, 18（5）: 1247-1279.
Böse M. 2006. Earthquake Early Warning for Istanbul Using Artificial Neural Networks[D]. Karlsruhe: Karlsruhe University: 19–24.
Böse M, Hauksson E, Solanki K, Kanamori H, Heaton T H. 2009a. Real-time testing of the on-site warning algorithm in southern California and its performance during the July 29 2008 MW5.4 Chino Hills earthquake[J]. Geophys Res Lett, 36（5）: L00B03.
Böse M, Sokolov V, Wenzel F. 2009b. Shake map methodology for intermediate-depth Vrancea （Romania） earthquakes[J]. Earthq Spectra, 25（3）: 497-514.
Böse M, Heaton T H, Hauksson E. 2012. Real-time finite fault rupture detector （FinDer） for large earthquakes[J]. Geophys J Int, 2012, 191（2）: 803-812.
Chen D Y. 2015. Development and Study of Earthworm Platform for Earthquake Early Warning in Taiwan[D]. Taipei: College of Science, Taiwan University: 43–47.
Espinosa-Aranda J M, Jiménez A, Ibarrola G, Alcantar F, Aguilar A, Inostroza M, Maldonado S. 1995. Mexico city seismic alert system[J]. Seismol Res Lett, 66（6）: 42-53.
Espinosa-Aranda J M, Cuellar A, Garcia A, Ibarrola G, Islas R, Maldonado S, Rodriguez F H. 2009. Evolution of the Mexican seismic alert system （SASMEX）[J]. Seismol Res Lett, 80（5）: 694-706.
Fletcher J B, Spudich P, Baker L M. 2006. Rupture propagation of the 2004 Parkfield, California, earthquake from observations at the UPSAR[J]. Bull Seismol Soc Am, 96（4B）: S129-S142.
Horiuchi S, Negishi H, Abe K, Kamimura A, Fujinawa Y. 2005. An automatic processing system for broadcasting earthquake alarms[J]. Bull Seismol Soc Am, 95（2）: 708-718.
Hoshiba M, Iwakiri K, Hayashimoto N, Shimoyama T, Hirano K, Yamada Y, Ishigaki Y, Kikuta H. 2011. Outline of the 2011 off the Pacific coast of Tohoku earthquake （Mw 9.0）: Earthquake early warning and observed seismic intensity[J]. Earth Planets Space, 63（7）: 547-551.
Kodera Y, Saitou J, Hayashimoto N, Adachi S, Morimoto M, Nishimae Y, Hoshiba M. 2016. Earthquake early warning for the 2016 Kumamoto earthquake: Performance evaluation of the current system and the next-generation methods of the Japan Meteorological Agency[J]. Earth Planets Space, 68（1）: 202-213.
Meng L, Ampuero J P, Sladen A, Rendon H. 2012. High-resolution backprojection at regional distance: Application to the Haiti M7.0 earthquake and comparisons with finite source studies[J]. J Geophys Res, 117（B4）: B04313.
Meng L, Allen R M, Ampuero J P. 2014. Application of seismic array processing to earthquake early warning[J]. Bull Seismol Soc Am, 104（5）: 2553-2561.
Nakamura Y. 1988. On the urgent earthquake detection and alarm system （UrEDAS）[C]//Proceedings of 9th World Conference on Earthquake Engineering （Vol. VII）. Tokyo-Kyoto, Japan: Japan Association for Earthquake Disaster Prevention: 673.
Odaka T. 2003. A new method of quickly estimating epicentral distance and magnitude from a single seismic record[J]. Bull Seismol Soc Am, 93（1）: 526-532.
Pazos A, Romeu N, Lozano L, Colom Y, López M, Goula X, Jara J A, Cantavella J V, Zollo A, Hanka W, Carrilho F. 2015. A regional approach for earthquake early warning in South West Iberia: A feasibility study. Bull Seismol Soc Am, 105（2A）: 560-567.
Peng H S, Wu Z L, Wu Y M, Yu S M, Zhang D N, Huang W H. 2011. Developing a prototype earthquake early warning system in the Beijing capital region[J]. Seismol Res Lett, 82（3）: 394-403.
Pittore M, Bindi D, Stankiewicz J, Oth A, Wieland M, Boxberger T, Parolai S. 2014. Toward a loss-driven earthquake early warning and rapid response system for Kyrgyzstan （Central Asia）[J]. Seismol Res Lett, 85（6）: 1328-1340.
Vallée M, Landès M, Shapiro N M, Klinger Y. 2008. The 14 November 2001 Kokoxili （Tibet） earthquake: High-frequency seismic radiation originating from the transitions between sub-Rayleigh and supershear rupture velocity regimes[J]. J Geophys Res, 113（B7）: B07305.
Wells D L, Coppersmith K J. 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement[J]. Bull Seismol Soc Am, 84（4）: 974–1002.
Wu Y M, Chung J K, Shin T C, Hsiao N C, Tsai Y B, Lee W H K, Teng T L. 1999. Development of an integrated earthquake early warning system in Taiwan: Case for the Hualien area earthquakes[J]. Terr Atmos Oceanic Sci, 10（4）: 719-736.
Wu Y M, Chen D Y, Lin T L, Hsieh C Y, Chin T L, Chang W Y, Li W S, Ker S H. 2013. A high-density seismic network for earthquake early warning in Taiwan based on low cost sensors[J]. Seismol Res Lett, 84（6）: 1048-1054.
Zhang Y, Wang R J, Zschau J, Chen Y T, Parolai S, Dahm T. 2014. Automatic imaging of earthquake rupture processes by iterative deconvolution and stacking of high-rate GPS and strong motion seismograms[J]. J Geophys Res, 119（7）: 5633-5650.
Zollo A, Lancieri M, Nielsen S. 2006. Earthquake magnitude estimation from peak amplitudes of very early seismic signals on strong motion records[J]. Geophys Res Lett, 33（23）: 430-452.
Zollo A, Amoroso O, Lancieri M, Wu Y M, Kanamori H. 2010. A threshold-based earthquake early warning using dense accelerometer networks[J]. Geophys J Int, 183（2）: 963-974.
-
期刊类型引用(1)
1. 王雪，李欣刚，闫科宇，杨梦雪. 空间扭索面曲塔斜拉桥振动破损有限元仿真. 计算机仿真. 2024(09): 101-104+204 . 
百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 1390
- HTML全文浏览量: 820
- PDF下载量: 71
- 被引次数: 2
