許亞吉 龐衛(wèi)東

摘要：介紹了烈度速報子系統(tǒng)的產品產出流程、云南地區(qū)臺站接入情況、儀器烈度算法，并以云南漾濞MS6.4地震為例，從烈度速報子系統(tǒng)的輸入數(shù)據(jù)、產出情況、儀器烈度與宏觀調查烈度的對比介紹該系統(tǒng)的應用情況。該系統(tǒng)產出的震中附近儀器烈度與宏觀調查烈度在空間走向上基本吻合，二者在Ⅵ～Ⅷ度區(qū)儀器烈度與宏觀調查有一定程度的對應，但各烈度區(qū)涉及的分布范圍有所不同，儀器烈度能較為客觀地反映震區(qū)災害情況。

關鍵詞：漾濞MS6.4地震;烈度速報子系統(tǒng);儀器烈度

中圖分類號：P315.6?? 文獻標識碼：A?? 文章編號：1000-0666（2021）03-0407-07

0 引言

地震烈度是衡量地震對地表和工程建筑結構的影響和破壞程度的重要參數(shù)，通常由地震工作者在震后數(shù)天通過分析處理實地調查資料得到。儀器地震烈度是指利用儀器記錄到的地震動觀測數(shù)據(jù)，通過峰值加速度、峰值速度等地震動參數(shù)與地震烈度之間的統(tǒng)計關系估算出臺站附近的破壞情況，在震后數(shù)分鐘或數(shù)十分鐘便可得到，其快速、直觀的特點有助于抗震救災工作的及時展開。

目前國內外使用較多的儀器地震烈度算法有：美國Shakemap方法綜合利用加速度峰值、速度峰值計算儀器烈度（Wald et al，1999）;日本氣象廳（JMA）利用儀器記錄到的三分向記錄，經傅立葉變換、濾波等方法處理后，合成有效峰值加速度進而計算儀器烈度（Shabestari，Yamazaki，2001;Karim，Yamazaki，2002;解全才等，2020）。袁一凡（1998）提出地震烈度模糊評定方法，綜合利用峰值加速度、傅立葉譜卓越頻率、相對持時（地震動三要素）和相應頻率點對應的反應譜值確定強震臺的儀器烈度。崔建文等（2008）利用隨機有限斷層法合成觀測點的地震動，然后采用模糊評定方法計算各觀測點的儀器烈度。王玉石等（2008）利用譜烈度與地震烈度的相關性來計算強震動記錄對應的儀器烈度。金星等（2013）參考日本氣象廳儀器烈度的算法，充分考慮我國的房屋建筑類型和建筑材料的差異、主要建筑的自振周期，設計了純幅值的帶通濾波器，用該濾波器對三分量的加速度記錄濾波后合成有效峰值加速度。李亮等（2018）在研究了美國地質調查局和日本氣象廳的儀器烈度算法的基礎上，對地震動數(shù)據(jù)記錄進行濾波處理，之后將三分量的加速度和速度進行合成后的峰值代入經驗公式求出儀器烈度值。我國先后頒布實施《中國地震烈度表》（GB/T 17742—2008）、《地震儀器烈度表》（DB35/T 1308—2012）、《儀器地震烈度計算規(guī)程》（試行）（2015）、《中國地震烈度表》（GB/T 17742—2020）對地震動記錄的選取、處理、記錄轉換、數(shù)字濾波、記錄合成、計算合成峰值加速度和峰值速度及儀器烈度的方法做了明確規(guī)定。部分地區(qū)結合自身實際，又出臺了各自的地方標準，例如福建地區(qū)就有自己的地震儀器烈度表。

日本、美國、中國臺灣等國家和地區(qū)均建設了地震烈度速報系統(tǒng)或城市地震烈度速報臺網。隨著國家地震烈度速報與預警工程云南子項目的推進，云南省地震局于2020年7月部署由福建省地震局技術團隊研發(fā)的預警項目定制軟件——烈度速報子系統(tǒng)。試運行1年，該系統(tǒng)對責任區(qū)內M≥3.0地震均有響應，產出了儀器烈度圖、峰值加速度等值線圖、峰值速度等值線圖等一系列烈度速報產品。本文從烈度速報產品的產出流程、臺站接入情況、儀器烈度算法對該系統(tǒng)進行介紹，并以2021年5月21日云南漾濞MS6.4地震為例，從烈度速報子系統(tǒng)的輸入數(shù)據(jù)、產出情況、儀器烈度與宏觀調查烈度的對比等方面介紹其實際應用情況。

1 地震烈度速報子系統(tǒng)

1.1 烈度速報產品產出流程

烈度速報子系統(tǒng)通過核心處理模塊從數(shù)據(jù)流服務器接收實時波形數(shù)據(jù)，在接收到來自EQIM平臺的地震三要素信息后進行觸發(fā)計算，得到各觀測臺站儀器烈度數(shù)據(jù)。儀器烈度繪圖服務接收核心處理程序的計算結果后，調用GMT進行繪圖，生成儀器烈度分布圖、峰值加速度等值線圖、峰值速度等值線圖、加速度反應譜等值圖等數(shù)據(jù)產品，并通過消息中間件等對外服務。同時烈度速報客戶端接收計算結果，在客戶端進行展示。烈度速報產品產出流程如圖1所示。

1.2 臺站接入情況

目前烈度速報子系統(tǒng)接入的云南地區(qū)臺站為測震臺站和預警一般站，其中測震臺站68個、一般站1 230個，臺站平均臺間距約為15 m，臺站分布如圖2所示。臺站觀測數(shù)據(jù)通過公用網絡，實時將數(shù)據(jù)傳回中心。測震臺站配備寬頻帶和甚寬頻帶地震儀，記錄數(shù)據(jù)為地震動速度，采樣率為100 sps，使用GPS進行授時。一般站設備為泰德TMA-33、瑞琪VH-GL-LDY01、港震GL-P2B、臺灣Palert Advance四種烈度儀，主要技術參數(shù)見表1，其傳感器為微機械加速度傳感器（MEMS加速度傳感器），記錄數(shù)據(jù)為地震動加速度，采樣率為100 sps，使用NTP服務（網絡時間協(xié)議）進行授時。MEMS加速度傳感器是微電子技術與微機械工程結合發(fā)展的一種慣性器件，與傳統(tǒng)加速度計相比體積更小、成本更低，利于大面積布設。在云南地震局預警項目攻堅團隊和云南鐵塔公司的共同努力下，云南地區(qū)一般站運行狀態(tài)良好，臺站運行率保持在97%以上，數(shù)據(jù)傳輸延時在2 s以內的臺站達到96%以上。

1.3 儀器烈度計算方法

使用臺站實測記錄，在進行基線校正、記錄轉換、數(shù)字濾波后，計算合成峰值加速度和合成 峰值速度，根據(jù)儀器烈度計算的國家標準（或地方標準）得到單個測點的儀器烈度值。利用震中附近臺站資料擬合烈度衰減關系，臺站的選擇遵循以下條件：M<5.5地震選取震中距在100 km內臺站，5.5≤M<6.5地震選取震中距150 km內臺站，M≥6.5地震選取震中距200 km內臺站。利用得到的衰減關系插值結合臺站實際烈度得出儀器烈度分布圖（金星等，2010，2013;張紅才，王士成，2008）。臺站對最近的插值網格的權重影響較大，對于臺站較稀疏的地區(qū)，以衰減關系為主;對于臺站較密集地區(qū)（平均臺間距在12 km以內），使用各測點的儀器烈度已基本可以得到烈度分布情況（王士成等，2017;徐欽等，2019）。單個測點的儀器烈度值在作為輸入數(shù)據(jù)計算震中附近烈度時，需對臺站的可信度做篩選，對可信度較低的臺站采取權重為零處理。篩選條件為：①利用臺站波形記錄的信噪比和功率譜密度（PSD）對異常臺站進行篩查，剔除信噪比低、功率譜密度曲線異常臺站;②根據(jù)震中附近烈度衰減關系，用實際計算值與理論計算值做比較，剔除烈度偏差超過2度的臺站。

2 烈度速報子系統(tǒng)在漾濞MS6.4震中的檢驗

2.1 產出情況

2021年5月21日21時48分云南漾濞發(fā)生MS6.4地震，22：00：48系統(tǒng)產出第一張烈度分布圖（圖3a），隨后陸續(xù)產出震中附近儀器烈度分布圖（3b）、峰值加速度等值線圖（圖4a）、峰值速度等值線圖（圖4b）、加速度反應譜等值圖（圖5）、2度以上鄉(xiāng)鎮(zhèn)烈度報告、2度以上縣市烈度報告。震后16 min產出所有烈度速報產品。

烈度速報子系統(tǒng)在接收到來自EQIM平臺的地震信息后，選用位于（23.67°～28.17°N、97.37°～102.17°E）范圍內的臺站數(shù)據(jù)進行處理，通過上述的篩選方法最終選取476個臺站的數(shù)據(jù)作為輸入數(shù)據(jù)集。選取震中距在150 km范圍內的臺站擬合烈度衰減關系，根據(jù)得到的衰減關系插值并結合臺站實際計算結果得出儀器烈度分布圖。對于此次地震，由于系統(tǒng)中沒有較為詳細的發(fā)震斷層信息和震源過程信息，故選用基于點源模型的衰減關系進行插值計算，參照EPRI（2003）推薦的方法將衰減關系式中的震中距轉換為等價斷層距（張紅才，王士成，2016）。地震動峰值（PGA、PGV）等值線圖通過擬合地震動衰減關系進行插值結果與臺站實測值得到。目前擬合地震動衰減關系時，臺站場地均選擇為Ⅰ類場地。后續(xù)隨著軟件本地化——臺站場地效應工作的持續(xù)開展，場地條件將按照臺站實際場地選擇。

此次地震最高儀器烈度為Ⅷ度，主要涉及漾濞縣蒼山西鎮(zhèn)、太平鄉(xiāng)2個鄉(xiāng)（鎮(zhèn)）。儀器烈度Ⅶ度區(qū)主要涉及漾濞縣漾江鎮(zhèn)、瓦廠鄉(xiāng)、平坡鎮(zhèn)、龍?zhí)多l(xiāng)、富恒鄉(xiāng)、順濞鎮(zhèn)，永平縣北斗彝族鄉(xiāng)，巍山縣紫金鄉(xiāng)以及大理市太邑鄉(xiāng)9個鄉(xiāng)（鎮(zhèn)）。儀器烈度Ⅵ度區(qū)主要涉及漾濞縣、洱源縣、永平縣、彌渡縣、巍山縣、賓川縣、大理市、云龍縣、保山隆陽區(qū)、昌寧縣的42個鄉(xiāng)（鎮(zhèn)）。

2.2 儀器烈度與宏觀調查烈度對比

圖6為云南省地震局現(xiàn)場工作隊通過宏觀調查繪制的漾濞MS6.4地震烈度分布圖。地震烈度Ⅷ度區(qū)主要涉及漾濞縣蒼山西鎮(zhèn)、漾江鎮(zhèn)、太平鄉(xiāng)3個鄉(xiāng)鎮(zhèn)。地震烈度Ⅶ度區(qū)主要涉及漾濞縣蒼山西鎮(zhèn)、漾江鎮(zhèn)、富恒鄉(xiāng)、太平鄉(xiāng)、順濞鎮(zhèn)、平坡鎮(zhèn)，云龍縣團結鄉(xiāng)，洱源縣西山鄉(xiāng)、煉鐵鄉(xiāng)，巍山縣紫金鄉(xiāng)和大理市太邑鄉(xiāng)等11個鄉(xiāng)（鎮(zhèn)）。Ⅵ度區(qū)主要涉及漾濞縣、洱源縣、云龍縣、永平縣、巍山縣、大理市的42個鄉(xiāng)（鎮(zhèn)）和2個街道辦事處。圖3b震中附近儀器烈度在空間上呈北西向分布，與宏觀調查烈度走向基本吻合。儀器烈度分區(qū)與各宏觀調查烈度分區(qū)具有一定程度的對應關系，但在面積和分布范圍上有所不同。這主要是因為儀器烈度計算的是地震動的物理參數(shù)，直接由臺站的實測記錄計算得到;宏觀調查烈度是通過房屋破壞程度等受災情況得到，受地震動強度、房屋結構等因素影響，故二者不可能完全一致。

表2給出了位于宏觀烈度Ⅵ度區(qū)以上的28個臺站記錄的儀器烈度。將各臺站儀器烈度與宏觀調查烈度進行對比可以看出，Ⅵ度區(qū)內儀器烈度值具有一定的離散性，與宏觀調查烈度值偏差在±1度以內的占66.7%，Ⅶ～Ⅷ度區(qū)內儀器烈度值與宏觀調查烈度值偏差均在±1度以內。如果儀器烈度值與宏觀調查烈度值偏差在±1度以內是可接受的，對儀器烈度值四舍五入取整后，二者有一定程度的對應，能較為客觀地反映震區(qū)災害情況。

3 結論

預警工程定制軟件——烈度速報子系統(tǒng)在云南部署近一年時間，對責任區(qū)內M≥3.0地震均有響應。通過漾濞MS6.4地震，該系統(tǒng)經受了一次很好的實戰(zhàn)檢驗。本文從輸入數(shù)據(jù)、烈度速報產出情況、儀器烈度與宏觀調查烈度的對比情況詳細介紹了烈度速報子系統(tǒng)在此次地震的應用。該系統(tǒng)在震后16 min快速產出了烈度速報產品。對比儀器烈度與宏觀調查烈度發(fā)現(xiàn)震中附近二者在空間走向上基本吻合，在Ⅵ～Ⅷ度區(qū)有一定程度的對應，但涉及的面積和分布范圍有所不同。總體上看，儀器烈度能較為客觀地反映震區(qū)災害情況。該系統(tǒng)產出的結果在地震應急初期可為抗震救援決策提供一定程度的科學依據(jù)，對地震現(xiàn)場工作隊開展現(xiàn)場災害調查有一定程度的指導意義。從圖件的實用性上來看，該系統(tǒng)產出的震中附近烈度分布圖未能包含震中附近全部Ⅵ區(qū)，圖幅自適應方面還有進一步優(yōu)化的空間。

感謝國家地震烈度速報與預警工程云南子項目攻堅團隊、云南省地震局漾濞地震現(xiàn)場工作隊及福建省地震局王士成在本文成稿過程中給予的幫助。兩位審稿專家對本文的完善提出了非常寶貴的意見，在此一并表示感謝。

參考文獻：

崔建文，盧大偉，高東，等.2008.基于合成地震動的震區(qū)烈度劃分[J].地震研究，31（4）：388-393.

金星，張紅才，李軍，等.2013.地震儀器烈度標準初步研究[J].地球物理學進展，28（5）：2336-2351.

金星，張紅才，韋永祥，等.2010.基于地震監(jiān)測臺網資料近實時插值計算震動圖的初步研究[J].防災減災學報，26（1）：1-11.

李亮，李山有，紀忠華，等.2018.儀器烈度計算方法研究[J].震災防御技術，13（4）：801-809.

王士成，金星，張紅才，等.2017.臺網密度對地震烈度速報的影響研究[J].地震工程與工程振動，37（6）：162-168.

王玉石，周正華，王偉.2008.基于假設檢驗的地震動強度（烈度）速報方法[J].地震工程與工程振動，28（5）：49-54.

解全才，馬 強，張景發(fā).2020.2016年日本熊本MW7.0地震強震動記錄特征分析[J].地震研究，43（1）：125-133.

徐欽，陳琳榮，張宇，等.2019.吉林松原MS5.7、MS5.1地震強震動記錄特征分析[J].地震研究，42（4）：627-634.

袁一凡.1998.由地震動三要素確定地震動強度（烈度）的研究[R].哈爾濱：國家地震局工程力學研究所，1-42.

張紅才，王士成.2016.地震烈度速報產品可靠性分析[J].地震工程與工程振動，36（5）：65-75.

張紅才.2008.基于地震監(jiān)測臺網資料的震動圖及震動烈度研究[D].哈爾濱：中國地震局工程力學研究所.

中國地震局.2015.儀器地震烈度計算規(guī)程（試行）[S].

Wald D J，Vincent Q，Heaton T H，et al.1999.Relationship between peak ground acceleration，peak ground velocity，and Modified Mercalli Intensity in California[J].Earthquake Spectra，15（3）：557-564.

Karim K R，Yamazaki F.2002.Correlation of JMA instrumental seismic intensity with strong motion parameters[J].Earthquake Engineering and Structural Dynamics，31（5）：1191-1212.

Shabestari K T，Yamazaki F.2001.A proposal of instrumental seismic intensity scale compatible with MMI evaluated from three-component acceleration records[J].Earthquake Spectra，17（4）：711-723.

GB/T 17742—2008，中國地震烈度表[S].

GB/T 17742—2020，中國地震烈度表[S].

DB35/T 1308—2012，地震儀器烈度表[S].

Application of the Rapid Intensity Reporting Subsystemto the 2021 Yangbi，Yunnan MS6.4 Earthquake

XU Yaji，PANG Weidong

（Yunnan Earthquake Agency，Kunming 650224，Yunnan，China）

Abstract

In this paper，we firstly introduced the Rapid Intensity Reporting Subsystem about its products output process，linked-up stations and algorithm of the instrumental intensity.Then we described the subsystem about its data input，intensity output，and the comparison between the instrumental intensity and the macroscopic intensity when the subsystem is used to handle the MS6.4 earthquake.We found that the spatial trend of the instrumental intensity was roughly consistent with the macroscopic intensity in the vicinity of the epicenter.The instrumental intensity and the macroscopic intensity in Ⅵ-Ⅷ areas to some extent corresponded to each other.In general，the instrumental intensity objectively reflected the damage in the disaster area.

Keywords：the Yangbi MS6.4 earthquake;the Rapid Intensity Reporting Subsystem;instrumental intensity