2022年9月5日瀘定地震震級的計算與對比分析

肖 勇 王旭亮 劉哲函 李 健 王曉明

（中國北京100085 禁核試北京國家數據中心）

0 引言

震級是地震的基本參數之一，用來表示地震的強度（陳運泰等，2004）。震級的測定是地震學研究的基本內容之一，可用于地震活動性、區域構造等的研究，根據震級還可以近似得到其他地震參數。此外，震級還是地震災后應急與救援決策的重要技術依據，快速準確地測定震級，有利于科學合理實施應急救援，減少人員傷亡和經濟損失（劉瑞豐，2018）。在禁核試核查領域，地震震級是估算地下核試驗爆炸當量的關鍵數據。

Richter 等（1935）最早引入了地方性震級ML的概念，其定義為地震信號幅值的對數加上與震中距相關的量規函數。在此基礎上，地震學家發展出了體波震級mb、面波震級MS、矩震級MW等震級概念。雖然這些震級的基本定義是一致的，但不同機構在計算不同震級時數據處理的方式及公式參數并不盡相同，如對于mb震級，全面禁止核試驗條約組織（CTBTO）國際數據中心（IDC）在測量幅值時濾波頻帶為0.8—4.5 Hz，而美國國家地質調查局（USGS）則采用2 個通帶分別為1.05—2.65 Hz、0.5—6.5 Hz 的頻帶（靳平，2022）；此外，這2 個機構采用的量規函數也不相同。這些計算上的不同可能會導致不同機構發布的震級結果出現較大差異。對于2022 年9 月5 日四川瀘定地震，IDC 地震事件公報中此次事件的震級為mb5.4，中國地震臺網中心（CENC）發布的M震級為6.8，USGS 發布的震級為mb6.6。禁核試北京國家數據中心（NDC）也監測到了此次地震，給出的震級為mb5.6。由此可見，相比于USGS 和CENC，IDC 和NDC 的體波震級mb明顯偏低，其中，NDC采用了與IDC 相同的震級算法，但所用臺站不同，因此震級結果較一致，但略有偏差。

不同地震監測機構對地震震級的定義和計算方法的不同，導致對同一事件不同機構發布的震級也不盡相同。本文系統梳理了國內外各主要地震監測機構的震級計算方法，考慮到實際中常用震級為體波震級mb，重點對不同機構發布的mb震級測定方法進行總結，對比IDC、USGS 發布的歷史地震事件mb震級的差異，對IDC 發布的mb震級偏低的原因進行分析，并以四川瀘定地震為例，基于NDC 現有地震監測臺網計算不同機構發布的mb震級并與其他機構發布的震級進行了對比分析。

1 不同震級計算方法比較

1.1 中國地震局震級

我國最早實施的震級國家標準是GB 17740—1999《地震震級的規定》（國家質量技術監督局，2004），在此基礎上規定了地方性震級ML、水平向面波震級MS、垂向面波震級MS7、短周期體波震級mb、中長期體波震級mB等5 種震級。隨著地震觀測技術的發展，震級的測定方法也發生了變化，如原國標震級定義均基于模擬記錄，而目前地震觀測數據基本均為數字記錄。此外，矩震級也逐漸被國際主流地震監測機構所使用。2013 年，國際地震學與地球內部物理協會（IASPEI）發布了新的震級標準，我國與國際主流地震監測機構之間在震級測定方面的差異逐漸體現。為此，中國地震局制定了新的震級國家標準B17740—2007《地震震級的規定》（中華人民共和國國家質量監督檢驗檢疫總局等，2017），并于2017 年發布。新國標中關于ML、MS、mb的定義與之前一致，但取消了mB和MS7，并增加了MS(BB)、MB(BB)、MW，這樣既保持了震級測定的連續性，又保持了與國際機構的一致性。

新國標中測定短周期體波震級mb時，將垂向的寬頻帶記錄仿真成DD-1 短周期地震儀記錄，進而測量P 波波列位移的最大值，并按下式計算震級

其中，A為質點位移最大值，單位為μm；T為對應周期；Δ為震中距；h為震源深度；Q為中國地震局量規函數。其中，DD-1 儀器響應函數為

體波震級mb是根據特定頻段的地震波振幅測定的，在對大地震的震級測量時會出現震級飽和現象，而矩震級MW則不存在震級飽和問題（陳運泰等，2018）。根據矩震級定義（Kanamori，1977），標量地震矩與矩震級間關系如下

其中，標量地震矩M0是衡量地震大小的物理量，單位為牛頓 ·米（N ·m），是由地震斷層的面積A、斷層的平均滑動量（平均錯距）D、斷層附近介質的剪切模量μ三者的乘積來定義的，即M0=μDA。

1.2 美國地質調查局震級

USGS 采用2 個通帶分別為1.05—2.65 Hz、0.5—6.5 Hz 的2 階Butterworth 帶通濾波器對寬頻帶速度記錄濾波，計算最大峰峰值的1/2（單位為μm）及相應周期T，然后根據測量的周期大小對2 個帶通濾波器的衰減作補償。得到A、T后，利用式（1）計算mb，其中，Q采用Granville 等（2022）提出的量規函數Q(GR)。

1.3 CTBT 國際數據中心震級

IDC 采用通帶為0.8—4.5 Hz 的3 階Butterworth 帶通濾波器對原始數據作正向、反向濾波，然后從P 波波前0.5 s 開始，取長度6 s 的數據測量峰峰值（單位nm）及相應周期T，之后利用式（1）計算mb，并采用Veith 等（1972）給出的量規函數Q(VC)。

1.4 IDC 震級與USGS 震級的歷史地震事件結果比較

統計了2021 年1 月1 日至5 月31 日IDC、USGS 地震事件公報結果，其中，IDC地震事件數18 396 個，USGS 地震事件數6 718 個，匹配地震事件數3 198 個，2 個機構mb震級對比如圖1 所示，震級偏差分布如圖2 所示。由圖1、2 可見，相比USGS 震級，IDC 的mb震級整體明顯偏小，平均偏差約0.56 個震級單位。

圖1 mb(IDC)與mb(USGS)的mb 震級對比Fig.1 Comparision of magnitude (mb) between IDC and USGS

圖2 IDC 與USGS 的mb 震級偏差分布Fig.2 Distribution of magnitude(mb) deviations between IDC and USGS

從IDC、USGS 的震級定義來看，其不同之處為：①測量的幅值不同，前者是單位為nm 的峰—峰值，后者是單位為μm 的峰—峰值的1/2，幅值測量方式的不同使得在測量同一信號時，IDC 測量的幅值對數lgA會比USGS 測量的lgA約大3.3 個震級單位；②2 個機構采用的量規函數也不相同，通過對比《新地震觀測手冊》（Bormann，2022）中2 個機構的量規函數可知，相同震中距和深度條件下，USGS 的Q(GR)約比IDC 的Q(VC)大3—4個單位，即USGS 對幅值測量的偏差進行了校正；③周期T的不同，由于測量頻帶不同，周期T也會有細微的不同。IDC、USGS 的mb震級偏差既受幅值測量方法不同的影響，又受量規函數不同的影響，二者間約0.5 級的震級偏差可以認為是由震級計算方式的差異所帶來的系統性偏差。此外，由于DD-1 的儀器響應函數在大于0.8 Hz 時幾乎是平坦的（接近1），故USGS、CENC 關于震級mb的幅值測量方式基本一致，只是采用的量規函數不同。根據新國標給出的量規函數，在相同震中距和深度時與Q(GR)相差不大，因此USGS和CENC 對此次瀘定地震給出的mb震級較一致。

2 瀘定地震震級計算結果分析

為對比不同機構的震級計算結果，NDC 結合自建的數據處理系統開發了一套震級計算程序，該程序可以同時計算IDC、USGS、國標（CENC）的震級，實現了在同一地震臺網條件下對不同機構發布震級的比對。

2.1 體波震級計算結果

對于此次瀘定地震，NDC 通過申請獲取的記錄到此次事件的臺站分布如圖3 所示。

圖3 NDC 通過申請獲取的監測到瀘定地震的臺站分布紅色五角星為瀘定地震震中Fig.3 Distribution map of stations monitoring Luding earthquake

基于獲取的臺站數據和震級計算程序，利用不同機構震級計算公式及參數的震級計算結果如圖4 所示。由圖4 可見，IDC 大部分臺站的震級相比于其他機構的明顯偏低。我們進一步統計了不同機構發布的震級及NDC 震級程序的計算結果（表1）（由于CENC 震級發布規則，此次事件未發布mb震級）。由表1 可見，NDC震級計算程序給出的結果與各機構發布的震級大體一致，但由于NDC 所用臺站不同，NDC 計算的USGS 震級與該機構發布的震級間有一定差異；NDC 臺網與IDC 臺網所包含的地震臺站有一定的重合，因此基于IDC 震級公式得到的震級與IDC 基本一致；由表1 還可見，USGS 的mb震級明顯大于IDC、NDC 的，這主要是由于對于大地震，USGS 測量P 波到時后60 s 內信號，而IDC、NDC 只用6 s 的信號。

表1 不同震級公式計算結果比較Table 1 The comparison of mb magnitude from different institutions

圖4 基于NDC 臺網計算的不同機構mb 震級Fig.4 Magnitude(mb) calculated of different institutions for NDC network

2.2 矩震級計算結果

利用gCAP 方法（Zhu et al，2013）反演了此次瀘定地震的震源機制和矩震級，反演時取Pz、Sh 的濾波頻段分別為0.02—0.10 Hz、0.02—0.05 Hz，得到震源機制為走向155°、傾角74°、滑動角-3°，矩震級6.60（圖5）。由圖5 可見，臺站的觀測波形與理論波形擬合程度較好，反演結果相對穩定。

圖5 瀘定地震震源機制反演結果Fig.5 Focal mechanism of the Luding earthquake

對于此次瀘定地震事件，USGS 利用W-phase 計算的Mww為6.64，中國地震局地球物理研究所張喆等計算的Mw為6.61，IDC 不發布Mw；NDC 震級計算程序給出的Mw為6.60，與USGS、CENC 給出的矩震級基本一致。

3 討論

以四川瀘定地震為例，梳理了國內外主要震級發布機構的mb震級定義，分析了不同機構之間mb震級存在偏差的原因。同時，對IDC、USGS 發布的歷史地震事件mb震級進行了統計對比，得到了2 個機構mb震級之間的系統偏差，并從震級定義上對這種系統偏差進行了分析。結果表明，不同機構之間mb震級的系統偏差主要是由于幅值測量方式和量規函數的不同所致，2 種因素綜合作用導致了最終的系統偏差。

此外，基于NDC 的震級計算程序和獲取的臺站數據計算了USGS、IDC、國標震級，發現同一震級計算公式在不同的臺網條件下震級計算結果大體一致，僅略有差異。該結果表明，即使震級定義一致，所用地震臺網不同也會造成最終的震級存在偏差，即事件震級計算結果會受到觀測臺網的影響。

綜上所述，不同機構mb震級存在偏差的根本原因是震級定義的不同，同時最終震級結果還會受觀測臺網的影響。因此，不同機構之間震級存在偏差是客觀的，只要每個機構的震級能保持一致性和連續性，該震級就是合理的。在計算特定機構的震級時，應采取相對應的量規函數對震級進行校正，從而避免系統性的錯誤。