基于深度卷積神經(jīng)網(wǎng)絡(luò)的2022年日本福島7.4級(jí)地震震級(jí)估算

曾奎原，朱景寶，宋晉東，李山有

（1.中國(guó)地震局工程力學(xué)研究所地震工程與工程振動(dòng)重點(diǎn)實(shí)驗(yàn)室，黑龍江 哈爾濱 150080；2.地震災(zāi)害防治應(yīng)急管理部重點(diǎn)實(shí)驗(yàn)室，黑龍江 哈爾濱 150080）

引言

地震預(yù)警可以有效的減輕地震破壞所帶來(lái)的人員傷亡和財(cái)產(chǎn)損失，地震預(yù)警通過(guò)P波和S波的波速不同的原理（即P波波速小于S波波速），利用震源周?chē)荚O(shè)的地震監(jiān)測(cè)臺(tái)站獲得的初至P波攜帶的信息估算地震基本參數(shù)和破壞區(qū)域，搶在破壞性地震波S波到來(lái)之前，向公眾發(fā)出警報(bào)，提醒其采取緊急避險(xiǎn)措施，為人們提供幾秒甚至數(shù)十秒的避險(xiǎn)時(shí)間［1－2］，進(jìn)而減少地震所造成的損失。美國(guó)、日本、意大利、墨西哥等［3－17］地震活動(dòng)多發(fā)地帶已經(jīng)建立了較為完善的地震預(yù)警系統(tǒng)。同時(shí)，中國(guó)地震預(yù)警系統(tǒng)［18－19］也在不斷進(jìn)行驗(yàn)證測(cè)試與更新完善。

對(duì)于地震預(yù)警系統(tǒng)而言，震級(jí)估算的結(jié)果決定著地震破壞區(qū)估算的準(zhǔn)確性和預(yù)警信息發(fā)布的可靠性，因此精確并且快速地估算出震級(jí)，是地震預(yù)警中尤為重要的一個(gè)環(huán)節(jié)。震級(jí)估算的傳統(tǒng)方法主要通過(guò)P波初至后的單一特征參數(shù)和震級(jí)的經(jīng)驗(yàn)關(guān)系來(lái)估算出震級(jí)［20］。常用的特征參數(shù)包括幅值類(lèi)參數(shù)［21－22］、周期類(lèi)參數(shù)［5，13－14］和能量類(lèi)參數(shù)［9，23－24］。雖然這上述幾類(lèi)特征參數(shù)已經(jīng)被前人研究所驗(yàn)證和震級(jí)存在著一定的相關(guān)性，但是由于單獨(dú)采用某一類(lèi)特征參數(shù)進(jìn)行震級(jí)估算時(shí)，其所攜帶的與震級(jí)相關(guān)的信息較為單一，所以在估算震級(jí)的時(shí)候容易產(chǎn)生較大誤差，且會(huì)發(fā)生小震高估、大震低估的現(xiàn)象。近年來(lái)，提取多種特征參數(shù)、利用人工智能方法建立震級(jí)估算模型快速發(fā)展，極大的減小了震級(jí)估算的誤差［25－27］，緩解了小震高估與大震低估的現(xiàn)象。

日本當(dāng)?shù)貢r(shí)間2022年3月16日23時(shí)36分，日本福島縣近海發(fā)生MJMA7.3級(jí)地震，震中位于北緯37°41.8′，東經(jīng)141°37.3′，震源深度60 km，最大觀測(cè)烈度（日本烈度）為6強(qiáng)。次日在日本氣象廳發(fā)布的震源參數(shù)報(bào)告中將此次地震的震級(jí)由MJMA7.3級(jí)提升至MJMA7.4級(jí)，震源深度變?yōu)?7 km。據(jù)相關(guān)媒體報(bào)道，在此次地震中，福島縣和宮城縣是地震活動(dòng)發(fā)生的主要地區(qū)，有較為強(qiáng)烈的震感，此次地震造成了1人不幸遇難，超80人受傷。在福島縣等地，公路路面發(fā)生了龜裂和扭曲、建筑物窗玻璃碎裂，部分地區(qū)道路交通臨時(shí)封閉。此外，地震還造成福島、宮城以及東京等多個(gè)地區(qū)地鐵路和地下鐵交通臨時(shí)中斷。

文中主要采用朱景寶等［28］建立的深度卷積神經(jīng)網(wǎng)絡(luò)DCNN－M模型，分別使用P波到達(dá)后1～40 s不同時(shí)間窗下的特征參數(shù)［29］作為輸入，對(duì)2022年3月16日發(fā)生在日本福島縣近海的地震進(jìn)行單臺(tái)震級(jí)估算，并采用多臺(tái)平均方法實(shí)時(shí)估算事件震級(jí)，驗(yàn)證該模型對(duì)此次日本福島縣近海地震震級(jí)估算的可行性。

1 方法

1.1 DCNN－M模型

文中所使用的模型是朱景寶等通過(guò)選用日本K－net強(qiáng)震臺(tái)網(wǎng)中心中所記錄的MJMA3～7.5地震事件并基于深度卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練而得到的。且在其研究中表明，對(duì)于所采用的地震事件，在P波初至后3 s的時(shí)間窗內(nèi)，該模型得到的震級(jí)估算結(jié)果相比于傳統(tǒng)的單參數(shù)方法τc和Pd方法有著更高準(zhǔn)確性，且小震高估、大震低估的現(xiàn)象得到了明顯的改善。

DCNN－M模型的基本網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示，其中主要包括輸入層、4個(gè)卷積層、4個(gè)BN層、4個(gè)池化層以及3個(gè)全連接層和輸出層，且在每個(gè)BN層和池化層之間都有插入一個(gè)激活函數(shù)用于對(duì)卷積層所輸出的結(jié)果做非線性映射，以提高模型的非線性表達(dá)能力。

圖1 DCNN－M模型網(wǎng)絡(luò)結(jié)構(gòu)Fig.1 The network structure of DCNN－M model

1.2 特征參數(shù)

DCNN－M模型所選取的是與震級(jí)相關(guān)的4類(lèi)共12個(gè)特征參數(shù)，即幅值參數(shù)（峰值位移Pd、峰值速度Pv、峰值加速度Pa）、能量參數(shù)（速度平方積分IV2、累積絕對(duì)速度CAV、累積能量變化率DI）、周期參數(shù)（特征周期τc、峰值比Tva、構(gòu)造參數(shù)）和衍生參數(shù)（豎向累積絕對(duì)位移cvad、豎向累積絕對(duì)速度cvav、豎向累積絕對(duì)加速度cvaa）作為輸入。同時(shí)，為了消除作為模型輸入的幅值類(lèi)參數(shù)、能量類(lèi)參數(shù)以及衍生類(lèi)參數(shù)因距離衰減而對(duì)震級(jí)估計(jì)產(chǎn)生的影響，需要對(duì)這些參數(shù)進(jìn)行震源距的校正，統(tǒng)一校正到參考震源距10 km［30－31］。

1.3 實(shí)時(shí)多臺(tái)平均

基于地震預(yù)警的時(shí)效性分析，文中通過(guò)該模型對(duì)P波最先觸發(fā)的10個(gè)臺(tái)站進(jìn)行震級(jí)估算，并將實(shí)時(shí)多臺(tái)平均的結(jié)果作為實(shí)時(shí)的震級(jí)估計(jì)。實(shí)時(shí)多臺(tái)震級(jí)平均計(jì)算公式如下：

式中：Mi為第i秒時(shí)的平均估算震級(jí)；是第i秒時(shí)單個(gè)臺(tái)站的估算震級(jí)；N是第i秒時(shí)參與平均計(jì)算的臺(tái)站數(shù)量。

2 強(qiáng)震數(shù)據(jù)及處理

文中所采用的數(shù)據(jù)是由日本防災(zāi)科學(xué)技術(shù)研究所K－net強(qiáng)震觀測(cè)臺(tái)網(wǎng)在本次地震中所獲取的加速度記錄，共576組（1 728條）強(qiáng)震記錄，同時(shí)對(duì)這些強(qiáng)震記錄做如下處理：

（1）首先采用馬強(qiáng)等［32］提出的P波到時(shí)自動(dòng)撿拾方法對(duì)強(qiáng)震加速度記錄進(jìn)行P波一次撿拾，然后對(duì)一次自動(dòng)撿拾的P波記錄繼續(xù)2次人工識(shí)別確認(rèn)；

（2）隨后對(duì)撿拾后的強(qiáng)震動(dòng)加速度記錄先進(jìn)行一次積分得到速度記錄，然后對(duì)積分后得到的速度記錄再一次積分得到位移記錄。最后為了消除積分帶來(lái)的低頻漂移的影響，對(duì)得到的記錄用4階0.075 Hz巴特沃斯高通濾波器進(jìn)行濾波［33－34］。

通過(guò)以上處理并篩選后得到的強(qiáng)震記錄共374組（1 122條）。

圖2展示了所有臺(tái)站記錄的三分向PGA與震中距的關(guān)系，圖中紅圈標(biāo)注的為首臺(tái)FKS005和加速度記錄最大臺(tái)站FKS002的PGA，藍(lán)圈標(biāo)注的為剩下臺(tái)站PGA。圖1和圖2表明：此次地震中臺(tái)站所記錄到的峰值加速度較大值主要集中在震源周?chē)?，且隨著震中距的增加，峰值加速度呈現(xiàn)不斷衰減的趨勢(shì)。

圖2 三分向PGA和震中距的關(guān)系Fig.2 The relationship of three-component PGA and epicenter distance

圖3則展示了本次地震事件中的強(qiáng)震記錄數(shù)量與震中距的關(guān)系，其中震中距的范圍在63.27～761.36 km，臺(tái)站主要集中分布在200～350 km的震中距范圍內(nèi)，有14個(gè)臺(tái)站位于100 km震中距以?xún)?nèi)，其中距離震中最近的臺(tái)站為FKS005，震中距為63.27 km，該臺(tái)站記錄的最大加速度幅值為NS向的607.817 8 cm/s2，而獲得最大加速度幅值記錄的臺(tái)站是FKS002，震中距為97.95 km，最大加速度幅值為EW向的－751.501 cm/s2。

圖3 震中距和強(qiáng)震動(dòng)記錄數(shù)量的關(guān)系Fig.3 The relationship of epicenter distance and strong-motion records numbers

圖4（a）和（b）分別展示了這次地震事件中首個(gè)觸發(fā)的臺(tái)站FKS005和所記錄到加速度最大的臺(tái)站FKS002的UD、EW、NS方向的加速度時(shí)程波形圖，峰值加速度以及P波到時(shí)標(biāo)注。

圖4 臺(tái)站FKS005和臺(tái)站FKS002三分向加速度記錄Fig.4 Three-component acceleration records of station FKS005 and station FKS002

3 結(jié)果

3.1 DCNN－M模型結(jié)果

文中根據(jù)日本氣象廳提供的實(shí)時(shí)定位結(jié)果，將幅值類(lèi)參數(shù)、能量類(lèi)參數(shù)和衍生類(lèi)參數(shù)進(jìn)行震源距校正后，聯(lián)合頻率類(lèi)參數(shù)作為深度卷積神經(jīng)網(wǎng)絡(luò)震級(jí)估算模型DCNN－M的輸入進(jìn)行地震預(yù)警震級(jí)估算。圖6展示了DCNN－M模型在首臺(tái)觸發(fā)后1～40 s內(nèi)對(duì)2022年3月16日日本福島縣近海MJMA7.4地震的離線數(shù)值模擬得到的震級(jí)估計(jì)結(jié)果。其中紅色實(shí)心圓代表的是DCNN－M模型震級(jí)估計(jì)的結(jié)果，黑色虛線代表此次地震的實(shí)際震級(jí)，藍(lán)色虛線代表±1級(jí)誤差，紫色虛線代表±0.5級(jí)誤差。表1給出了日本氣象廳緊急地震速報(bào)詳情給出的速報(bào)結(jié)果，以及對(duì)應(yīng)于DCNN－M模型給出的震級(jí)預(yù)測(cè)結(jié)果。當(dāng)估計(jì)的地震震級(jí)超過(guò)特定閾值（MJMA3.5）時(shí)，日本氣象廳就會(huì)向高級(jí)用戶(hù)發(fā)布預(yù)警信息［35］。由于本次地震日本氣象廳給出的震級(jí)結(jié)果存在錯(cuò)誤，都為1.0級(jí)，因此無(wú)法將DCNN－M模型的結(jié)果與日本氣象廳的震級(jí)結(jié)果作出直接一對(duì)一比較，因此文中以日本氣象廳緊急地震速報(bào)每一報(bào)發(fā)布時(shí)刻為對(duì)照，只分析DCNN－M模型在日本氣象廳緊急地震速報(bào)的每一報(bào)時(shí)刻震級(jí)估算結(jié)果。

表1 日本氣象廳和DCNN－M模型震級(jí)估算結(jié)果Table 1 The magnitude estimation results of Japan Meteorological Agency and the DCNN－M model

結(jié)合圖5和表1，DCNN－M模型在首臺(tái)觸發(fā)后1.6 s（即：日本氣象廳預(yù)警速報(bào)第1報(bào)），給出的震級(jí)估算結(jié)果為5.85級(jí)，且滿足日本氣象廳預(yù)警發(fā)布的閾值（MJMA3.5）。同時(shí)，隨著首臺(tái)觸發(fā)后時(shí)間的不斷增加，觸發(fā)的臺(tái)站數(shù)目也在增多，DCNN－M模型給出的震級(jí)估算結(jié)果也在迅速趨近于實(shí)際震級(jí)。在首臺(tái)觸發(fā)后4.2 s（即：日本氣象廳預(yù)警速報(bào)第3報(bào)）時(shí)，震級(jí)估計(jì)達(dá)到了6.47級(jí)，震級(jí)估計(jì)誤差小于1個(gè)震級(jí)單位；在首臺(tái)觸發(fā)后11 s，DCNN－M模型估算的結(jié)果就達(dá)到了6.94級(jí)，與實(shí)際震級(jí)的誤差已經(jīng)小于0.5個(gè)震級(jí)單位，并且在首臺(tái)觸發(fā)后23.4 s（即：日本氣象廳預(yù)警速報(bào)第21報(bào)）給出的震級(jí)估算結(jié)果為7.37級(jí)，與實(shí)際震級(jí)的誤差僅僅有0.03個(gè)震級(jí)單位。可以看出：對(duì)于此次地震，在日本氣象廳預(yù)警速報(bào)第1報(bào)時(shí)，DCNN－M模型震級(jí)估計(jì)滿足閾值；在首臺(tái)觸發(fā)初期，DCNN－M模型可以得到魯棒的震級(jí)估算，且隨著時(shí)間的增加，預(yù)測(cè)震級(jí)逐漸接近實(shí)際震級(jí)。

圖5 DCNN－M模型多臺(tái)平均震級(jí)估算結(jié)果Fig.5 Estimation results of multiple mean magnitudes in DCNN－M model

從表1中也可以得到：在日本氣象廳預(yù)警速報(bào)第3報(bào)時(shí)，該模型所給出的震級(jí)估算結(jié)果與實(shí)際震級(jí)的誤差已經(jīng)小于1個(gè)震級(jí)單位；在首臺(tái)觸發(fā)后11 s，日本氣象廳未有速報(bào)發(fā)布，但是依據(jù)DCNN－M模型可以得到6.94級(jí)，且預(yù)測(cè)震級(jí)與實(shí)際震級(jí)的誤差小于0.5個(gè)震級(jí)單位。

3.2 單臺(tái)結(jié)果與多臺(tái)結(jié)果對(duì)比

如圖6所示，給出了參與平均的單臺(tái)震級(jí)估算的結(jié)果。結(jié)合圖5和圖6可以得到，個(gè)別單臺(tái)估算的結(jié)果在首臺(tái)觸發(fā)初期估算的結(jié)果相比于多臺(tái)平均結(jié)果給出的震級(jí)要明顯過(guò)低，如果僅僅用單臺(tái)的估算結(jié)果，就容易導(dǎo)致估算結(jié)果趨近于實(shí)際震級(jí)的時(shí)間要增加，不利于地震預(yù)警所注重的時(shí)效性，因此用多臺(tái)估算結(jié)果參與平均能夠在首臺(tái)觸發(fā)初期快速的使得震級(jí)結(jié)果趨近于實(shí)際震級(jí)結(jié)果。同時(shí)，隨著觸發(fā)時(shí)間的增加，單臺(tái)估算的結(jié)果也在趨近于實(shí)際震級(jí)，但是部分臺(tái)站震級(jí)估算的結(jié)果超過(guò)了實(shí)際的震級(jí)，出現(xiàn)了高估的現(xiàn)象，以臺(tái)站FKS006尤為明顯，而反觀多臺(tái)平均的結(jié)果，并未出現(xiàn)高估的現(xiàn)象，且估算結(jié)果一直在趨近實(shí)際震級(jí)，并未超出。這是因?yàn)槎嗯_(tái)參與平均將部分高估的結(jié)果給降低了，使其能更好的向?qū)嶋H結(jié)果趨近且不至于產(chǎn)生與實(shí)際結(jié)果有較大的誤差。

圖6 單臺(tái)震級(jí)估算結(jié)果Fig.6 Magnitude estimation of single-station

3.3 DCNN－M模型與τc和Pd方法震級(jí)估算結(jié)果對(duì)比

如圖7所示，圖中分別給出了DCNN－M模型、τc和Pd方法對(duì)本次地震事件的多臺(tái)平均震級(jí)估算結(jié)果。從圖中可以看出，τc方法所預(yù)測(cè)的震級(jí)結(jié)果明顯較差，在首臺(tái)觸發(fā)初期震級(jí)估算結(jié)果過(guò)低，且隨著時(shí)間窗的增加，估算值有所改善但仍與實(shí)際震級(jí)有著很大的誤差，而Pd方法的震級(jí)估算結(jié)果比τc方法的估算結(jié)果有了很大程度上的改善，且在首臺(tái)觸發(fā)初期震級(jí)估算的結(jié)果甚至優(yōu)于文中所研究的DCNN－M模型給出的震級(jí)估算結(jié)果，但是隨著時(shí)間窗的增加，震級(jí)高估的現(xiàn)象愈發(fā)明顯，也未得到降低。傳統(tǒng)的單參數(shù)τc和Pd方法在預(yù)測(cè)震級(jí)時(shí)，分別出現(xiàn)了一些諸如震級(jí)低估高估，首臺(tái)觸發(fā)初期未能給出較好的預(yù)測(cè)結(jié)果等現(xiàn)象，而DCNN－M模型在這些方面都有了明顯的改善：首先在首臺(tái)觸發(fā)初期，給震級(jí)估算誤差就已經(jīng)滿足日本氣象廳預(yù)警信息發(fā)布的震級(jí)閾值（MJMA3.5），而τc方法在2 s給出的結(jié)果為MJMA2.95未能滿足閾值。其次，隨著首臺(tái)觸發(fā)時(shí)間的增加，DCNN－M模型的預(yù)測(cè)震級(jí)迅速趨近于實(shí)際震級(jí)，而Pd方法在首臺(tái)觸發(fā)后14 s的預(yù)測(cè)結(jié)果就已經(jīng)超出了實(shí)際震級(jí)，且在隨后高估的現(xiàn)象越發(fā)明顯。

圖7 DCNN－M、τc和Pd方法多臺(tái)平均震級(jí)估算結(jié)果Fig.7 Estimation results of multiple mean magnitudes in DCNN－M model、τc and Pd method

4 結(jié)論與討論

基于朱景寶等［29］使用日本K－net強(qiáng)震數(shù)據(jù)構(gòu)建的DCNN－M模型，本研究使用日本2022年3月16日福島縣近海MJMA7.4地震在P波到達(dá)后1～40 s不同時(shí)間窗下的特征參數(shù)作為DCNN－M模型的輸入，對(duì)此次地震進(jìn)行震級(jí)估算。通過(guò)分析單臺(tái)DCNN－M模型離線數(shù)值模擬的結(jié)果和對(duì)其結(jié)果進(jìn)行多臺(tái)平均，探索在該模型下對(duì)此次地震事件震級(jí)估算的可行性，并給出以下結(jié)論：

（1）在首臺(tái)觸發(fā)后1.6 s（即：日本氣象廳預(yù)警速報(bào)第1報(bào)），DCNN－M模型給出的震級(jí)估算結(jié)果為5.8級(jí)，且滿足日本氣象廳預(yù)警信息發(fā)布的震級(jí)閾值（MJMA3.5）；

（2）在首臺(tái)觸發(fā)后4.2 s（即：日本氣象廳預(yù)警速報(bào)第3報(bào)），DCNN－M模型給出了6.47級(jí)的估算結(jié)果，震級(jí)估算誤差小于1；

（3）在首臺(tái)觸發(fā)后11 s，DCNN－M模型的震級(jí)估算結(jié)果為6.94級(jí)，與實(shí)際震級(jí)偏差縮小到可容忍的0.5個(gè)震級(jí)誤差范圍內(nèi)；

（4）在首臺(tái)觸發(fā)后23.4 s（即：日本氣象廳預(yù)警速報(bào)第21報(bào)），DCNN－M模型的震級(jí)估算結(jié)果為7.37級(jí)，與實(shí)際震級(jí)僅有0.03個(gè)震級(jí)的誤差，并且在隨后的估算結(jié)果也與實(shí)際震級(jí)的誤差維持在0.1～0.2個(gè)震級(jí)范圍內(nèi)；

（5）通過(guò)與τc和Pd方法的多臺(tái)平均結(jié)果進(jìn)行對(duì)比，DCNN－M模型所給出的震級(jí)估算結(jié)果優(yōu)化了這2種傳統(tǒng)單參數(shù)方法中出現(xiàn)的震級(jí)低估高估，首臺(tái)觸發(fā)初期震級(jí)估計(jì)誤差較大現(xiàn)象。

基于地震預(yù)警的時(shí)效性分析，DCNN－M模型對(duì)本次日本福島縣近海地震，在首臺(tái)觸發(fā)后1.6 s，震級(jí)估算誤差滿足日本氣象廳預(yù)警信息發(fā)布的震級(jí)閾值（MJMA3.5）；在首臺(tái)觸發(fā)后4.2 s時(shí)，震級(jí)估算誤差小于1；在首臺(tái)觸發(fā)后11 s時(shí)，震級(jí)估算誤差小于0.5。這也表明：對(duì)于此次地震，在首臺(tái)觸發(fā)初期，DCNN－M模型可以得到魯棒的震級(jí)估計(jì)。同時(shí)，隨著時(shí)間窗的增加，震級(jí)估計(jì)誤差逐漸減小。

致謝：日本防災(zāi)科學(xué)技術(shù)研究所為本研究提供數(shù)據(jù)支持（https：//www.kyoshin.bosai.go.jp/kyoshin/quake/index_en.html）。