2022年青海門源M S6.9地震震源機制解＊

韓立波

（中國地震局地球物理研究所，北京100081）

引言

據中國地震臺網正式測定，2022年1月8日1時45分，青海海北州門源縣發生了MS6.9地震（37.77°N，101.26°E），震源深度10 km。這次地震是該地區近年來最大的一次地震，包括寧夏銀川市、中衛市及甘肅蘭州市在內的多個城市震感強烈，該地區歷史上強震頻發，青藏高原東北緣歷史上發生過多次7級以上的地震，其中1920年以來發生過3次7.5級以上的大地震，包括1920年海原M8.5地震、1927年古浪M8.0地震和1932年昌馬M7.6地震，它們分別沿海原斷裂、冷龍嶺斷裂和昌馬斷裂產生了顯著的地震地表破裂帶。新中國成立以來，在這一地區附近曾發生過1954年山丹M7.2地震，1986年門源M6.4地震，2016年門源再次發生M6.4地震，該地區的地震危險性也引起社會公眾和地震工作者的廣泛關注。

在區域構造位置上，此次門源MS6.9地震發生在青藏高原東北緣冷龍嶺斷裂、托萊山斷裂和肅南—祁連斷裂的階區部位，構造較為復雜（圖1）。本文利用CAP方法反演了主震震源機制解和震源深度，依據這些結果討論了該地區的活動構造特征，并指出其對于認識該地區地震危險性的潛在意義。

圖1 區域構造背景圖Fig.1 Tectonic settingsof the study region

1 方法

1.1 震源機制解反演

本文利用CAP方法[2-3]進行震源機制解反演，其主要思想是利用近震數據，把寬頻帶數字波形記錄分為體波部分（Pnl）和面波部分，分別計算它們的理論地震圖和實際觀測波形的目標誤差函數，在給定參數空間中進行網格搜索，同時反演震源機制解和震源深度。

考慮到因幾何擴散產生的衰減對波形的影響，使用經震中距矯正后的絕對誤差值作為目標誤差函數，定義為：

式中，r為震中距，r0為選定的參考震中距，p為比例因子，用以保障震中距r處的權重與r0相當。基于所定義的目標誤差函數，采用網格搜索的方法在M0、θ?φ、δ、λ以及震源深度空間進行搜索，得到最佳的震源機制解、矩震級和震源深度。

前人的研究結果表明，CAP方法在反演震源機制解和確定震源深度上優勢較為明顯[4-5]。

1.2 理論地震圖計算

在計算理論地震圖時，這里采用了目前廣泛使用的頻率—波數（F-K）法[6]。F-K法適用于水平分層地殼模型，通過對頻率和波數分別進行積分，采用傳播矩陣計算地震的全波場位移分布，能夠計算各種頻率下包含了體波和面波波形等成分的全部波形。速度結構模型使用了Crust 2.0（http://igppweb.ucsd.edu/~gabi/crust2.html）。

1.3 數據資料處理

本文使用了國家數字測震臺網數據備份中心[7]提供的青海、甘肅等區域臺網的寬頻帶數字地震波形記錄。按照方位角覆蓋及信噪比的要求，挑選出寬頻帶地震臺數據進行處理，圖1給出了臺站的分布。

對挑選出的寬頻帶數據首先去除儀器響應，并旋轉至大圓路徑。然后對波形中Pnl部分和面波部分使用帶寬為0.02—0.05 Hz的帶通濾波器進行濾波。這些頻帶在美國南加州震源機制解研究中采用，適合三維結構并不十分復雜的情形；而Pnl和面波的相對權重為2:1，一些研究表明這樣的權重可以較好兼顧Pnl和面波的優點。相應地，采用相同的濾波參數，對計算得到的理論地震圖進行濾波。

2 結果

2.1 震源持續時間反演結果

反演過程中考慮到主震震級較大，部分震中距較小臺站存在顯著的限幅，數據處理過程中仔細去除了限幅臺站。為了確定本次根據公式（1）給出的誤差目標函數，在參數全空間范圍內搜索最佳震源機制解、震源深度和矩震級。對于中強地震，由于其震源較為復雜，其震源破裂往往持續一定時間，稱為震源持續時間，與破裂尺度和破裂速度相關。為了獲得最佳震源持續時間，我們還搜索了不同震源持續時間的擬合誤差。圖2給出了反演擬合誤差隨著震源持續時間變化的結果，在震源持續時間過大或者過小時，誤差均顯著增加，當震源持續時間為7.3 s時，擬合誤差最小。

圖2 誤差隨著震源持續時間變化搜索結果圖Fig.2 Waveform fit errors as function of source duration time in 2022 Menyuan M S6.9 event moment-tensor inversion

2.2 震源機制解反演結果

我們選用震源持續時間為7.3 s，圖3給出了觀測波形與理論波形擬合誤差隨深度的分布，當地震震源深度為3 km時，擬合誤差達到最小，據此得到最佳地震震源深度為3 km，此時對應的雙力偶解即為最佳雙力偶解。其中，主震最佳雙力偶解為節面Ⅰ：走向191°、傾角62°、滑動角173°；節面Ⅱ：走向284°、傾角82°、滑動角21°。此外，反演得到的矩震級為Mw6.7。由圖3還可看出，隨著震源深度的改變，震源機制解變化并不顯著，說明反演得到的震源機制解較為穩定。

圖3 門源M S6.9地震矩張量反演中波形擬合誤差隨深度變化Fig.3 Waveform fit errors as function of depth in 2022 Menyuan M S6.9 event moment-tensor inversion

圖4給出了本次地震理論波形與觀測波形的擬合情況，大部分臺站各震相擬合效果較好。由于數據記錄質量差或區域地殼速度結構復雜等原因，簡單的一維速度結構模型在反演中不能完全適用，部分臺站或部分分量的波形擬合情況可能較差。為避免此類數據參與計算影響結果的可靠性，反演中未予使用。同時部分臺陣的面波也由于限幅未參與反演。

圖4 門源M S6.9地震矩張量反演理論地震波形（紅色）與實際觀測地震波形（黑色）波形圖。下方第一行數字為各段理論地震波形相對實際觀測波形的移動時間，正值表示理論波形相對觀測波形超前。第二行數字為理論波形與觀測波形的相關系數（百分比）。波形圖左側字母為臺站，其下數字分別為臺站震中距（km）以及理論地震圖相對實際觀測數據整體移動時間。震源球上黑色區域代表壓縮區，白色代表拉張區，紅色十字符號代表臺站。震源球采用下半球投影Fig.4 Comparison between synthetic（red）and observed （black）seismograms of 2022 Menyuan M S6.9 event.The numbers on the lower left side of the seismograms are the time shifts（upper）and cross-correlation coefficient in percent （lower）.Positive time shiftsmean that the observed data have been delayed.Theletters on the left side arestations，below which thenumbers are epicentral distances and time shifts between synthetics and the observed data.The black color in beach-ball denotes compression，while white color isextension.The red crosses are stations.Lower hemisphere projection is used

3 討論與結論

本文利用CAP方法反演了2022年青海門源MS6.9地震震源機制解，結果表明，本次地震為走滑型，這一結果與其他結果接近。Fan等[8]對主震后三天的余震進行了精定位，結果顯示余震序列分東西兩段，西段走向近90°，東段走向133°左右，推測余震序列西段位于托萊山斷裂，東段位于冷龍嶺斷裂，本次地震發生在托萊山斷裂和冷龍嶺斷裂交匯部位。前人認為冷龍嶺斷裂帶全新世活動強烈，晚第四紀主要表現為左旋走滑運動[9-13]，震源機制解的性質與冷龍嶺斷層活動性質吻合。震源機制解的節面Ⅰ走向也與余震序列的東段走向接近，因此，本次地震為冷龍嶺斷裂活動所致。值得注意的是，該地區1986年8月26日和2016年1月21日發生了2次M6.4地震，根據前人研究，1986年和2016年2次地震震源機制解為逆沖型震源機制解（https://www.globalcmt.org/）[14]。本次地震與1986年和2016年2次地震震源機制解差異較大，可能是因為本次地震發生在冷龍嶺斷裂上，而前2次M6.4地震發生在冷龍嶺北側斷裂上[15-16]，冷龍嶺斷裂以走滑運動為主，冷龍嶺北側斷裂是冷龍嶺斷裂西北端一條伴生斷裂，在地表表現為一條晚第四紀變形不明顯的逆斷裂，這也說明了該地區構造運動復雜。

本次地震震源深度為3 km左右，相對較淺。李振洪等[17]采用InSAR技術獲取了同震地表形變場，并基于分布式滑動模型反演了地震斷層面上的滑動分布，發現本次地震發震斷層的最大滑動值達到3.5 m，處于地下4 km左右。據現場考察，本次地震地表破裂帶可能長達22 km[18]，這對于如此震級的地震來說，是極為顯著的。相對而言，例如2013年蘆山M7.0地震矩震級為MW6.7，2017年九寨溝地震矩震級也為MW6.7，均未發現顯著的地表破裂。動力學模擬表明，地表破裂模式受多種因素影響，例如受震源位置影響，也就是受斷層帶上非均勻的應力場分布影響[19]，因此，這些地震盡管震級接近，但發生在不同的震源位置，且本次地震震源深度也較淺，其斷層面上應力場分布也不同，這些均可能造成地表破裂差異較大。

冷龍嶺、金強河、毛毛山、老虎山斷裂所在的區段是祁連—海原斷裂帶的地震空區（天祝地震空區）[20]，該地區的地震危險性長期受到關注。朱琳等[21]基于分層黏彈性流變模型計算了青藏高原北部1900年以來的強震對祁連—海原斷裂帶的庫侖應力加載，結果顯示，祁連—海原斷裂帶西段木里江倉斷裂和托萊山斷裂以及中段的金強河—老虎山斷裂應力增強顯著，最大庫侖應力加載可達1 MPa以上。因此，該區域未來地震危險性高，值得進一步關注。

致謝

中國地震局地球物理研究所“國家數字測震臺網數據備份中心”為本研究提供了地震波形數據。