許健生 李 麗 隗永剛 張 旭

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

0 引言

根據全球數字地震儀臺網（GSN）的測定結果，2016 年11 月13 日11 時02 分56.0 秒（GMT），在新西蘭南島（42.76°S，173.08°E，h=23.0 km）發(fā)生一次MW7.8 地震，震后發(fā)生浪高達3 m 的海嘯（Hamling et al，2017）。在震中距1.9°—144.5°范圍內，GSN 有100 個臺的寬頻帶（BH）和甚寬頻帶（VH）數字地震儀較好記錄到該地震。

李普春等（2013）在實驗室進行的巖石粘滑錯動實驗表明，粘滑錯動過程一般可分為預滑、粘滑和止滑3 個階段，且往往不是一次、單點錯動，而是由多次、多點粘滑錯動組成，表現出在斷層的不同部位、不同時間、多點、多次粘滑錯動的特征。在天然地震觀測中，破裂前的預滑過程激發(fā)出的信號較弱，一般不易被觀測到。但破裂后的同震粘滑和止滑過程激發(fā)出的信號較強，容易被觀測到。分析此次MW7.8 地震數字波形記錄后發(fā)現：在100個臺的長周期（LP）記錄上，在與P 震相同時到達的Xs1 震相之后，可進一步識別出子粘滑錯動震相Xs2、Xs3、Xs4 和粘滑錯動的峰值震相Xsm。在粘滑錯動達到峰值后還可識別出止滑過程激發(fā)出的長周期止滑面波震相XsQ（勒夫型）和XsR（瑞雷型）。其中，距震中最近的SNZO 地震臺（41.31°S，174.70°E，Δ=1.9°）記錄到的XsQ 和XsR 震相，周期長約600 s。對此次新西蘭南島MW7.8 地震的同震粘滑震相進行分析，以期在地震記錄圖上識別出Xs、XsQ 和XsR 震相，為認識同震粘滑錯動過程及預判地震次生災害提供幫助。

1 數據

本文所用資料來自全球地震臺網（GSN）100 個臺的寬頻帶（BH）和甚寬頻帶（VH）速度型數字記錄，拾震器型號多為STS-1 或STS-2。在震中距100°以內，用CDSN 地震臺站分析軟件將BH 通道的記錄仿真成短周期（SP）后讀取P 震相數據。在震中距100°—144.5°間，仿真成長周期記錄（LP）后讀取Pdif 震相到時。全部的P+Xs1、Xs2、Xs3 和Xs4 震相數據均為仿真成長周期記錄（LP）后讀取，在VH 通道讀取了全部的XsQ 和XsR震相數據。

地震參數取自美國地質調查局（USGS）國家地震信息中心（NEIC）發(fā)布的地震目錄，發(fā)震時間和地震圖上標注的時間均為世界標準時（UTC）。震中位置和GSN 臺站位置見圖1。

圖1 GSN 臺站位置和震中分布紅圓點為震中，藍三角為臺站Fig.1 Distribution of earthquake epicenter and GSN stations

2 Xs、Xsm、XsQ 和XsR 震相的記錄特征

在寬頻帶（BH）長周期（LP）和甚寬頻帶（VH）的數字波形記錄上，Xs、Xsm、XsQ 和XsR 震相有以下記錄特征。P+Xs1、Xs2、Xs3、Xs4 和Xsm 震相典型記錄（VH）見圖2，XsQ 和XsR 震相典型記錄見圖3。

圖2 P+Xs1、Xs2、Xs3、Xs4 和Xsm 震相典型記錄Fig.2 The typical recordings of P+Xs1,Xs2,Xs3,Xs4 and Xsm seismic phases

圖3 XsQ 和XsR 震相典型記錄Fig.3 The typical recordings of XsQ and XsR seismic phases

（1）由圖2 可見，Xs1 震相與P 震相同時到達，Xs2、Xs3 和Xs4 震相是破裂后平均約7.8 s、21.6 s 和 61.1 s 出現的長周期震相。P+Xs1 震相后平均約80.0 s 的Xsm 是Xs 的峰值震相。在長周期記錄上，Xsm 周期范圍在15.2—34.6 s，Xsm 震相平均周期約25.4 s。雖然P+Xs1、Xs2、Xs3 和Xs4 震相的振幅隨震中距增大而逐漸衰減，但此次地震，在Δ＜144.5°時都能很好地觀測到P+Xs1、Xs2、Xs3 和Xs4 震相。由圖2 可見，Xs2、Xs3和Xs4 震相與P+Xs1 震相之間的到時差是一定值，不隨震中距增大而變化。

（2）由圖3 可見，在震中距1.9°—144.5°間，低頻且大振幅的勒夫型面波XsQ 和瑞雷型面波XsR，不僅是整個地震記錄中周期最長，振幅最大的震相，而且是面波波列中速度較快的震相，其速度明顯大于在地殼中傳播的瑞雷型面波Rm。XsQ震相的周期范圍在36.7—600 s，平均周期約62.5 s，XsR 震相的周期范圍在32.6—600 s，平均周期約55.4 s。

3 Xs、XsQ 和XsR 震相的物理機制

早在20 世紀60 至70 年代，美國學者Brace 和Byerlee 就針對Reid 的“彈性回跳”學說與觀測事實相悖的問題，把“粘滑”作為一種地震機制提了出來，指出彈性回跳的物理機制可用摩擦滑動過程中出現的不平穩(wěn)滑動（即粘滑）來解釋（Brace et al，1966；Byerlee et al，1975；吳曉媧等，2016）。并根據巖石的壓力實驗結果指出，主破裂過程可由數次粘滑過程構成。

李普春等（2013）在實驗室進行的巖石粘滑錯動實驗表明，粘滑錯動過程一般可分為預滑、粘滑和止滑3 個階段，且錯動過程往往不是一次、單點錯動，而是由多次、多點粘滑錯動組成，表現出在斷層的不同部位、不同時間，多點、多次粘滑錯動的特征。

由GSN 100 個臺的寬頻帶（BH）和甚寬頻帶（VH）數字地震儀的記錄可見，Xs1 震相是與P 震相同時到達的震相，其生成機制可能是斷層面上的凹凸體開始破裂后，斷層由“粘結”的靜態(tài)轉變?yōu)椤盎瑒印钡膭討B(tài)，在粘滑錯動時輻射出的長周期波動。粘滑錯動是一個在構造力作用下不斷克服斷層面上磨擦阻力的過程。在粘滑錯動過程開始后，粘滑錯動的速度和位移有一個由慢到快、由小到大的過程。當構造力較小時斷層錯動速度慢，位移小。而當構造力可克服較大磨擦阻力時，斷層粘滑錯動速度加快，位移增大。當遇到另一凹凸體或磨擦阻力較大時又停滯或降速。當構造力再加大到可克服磨擦阻力時，斷層又向前滑動。如此，便形成了多次、多點的子滑動事件。Xsm 震相則是整個粘滑錯動過程中，粘滑錯動速度和位移達到峰值時的震相。在粘滑錯動達到峰值后，斷層粘滑錯動開始減速，進入止滑階段。李普春等（2013）的實驗結果表明，止滑過程可以激發(fā)長周期、大振幅的止滑波動震相。這種波動輻射出了長周期的勒夫型面波XsQ 和瑞雷型面波XsR。

4 結果討論

由時差值與震中距的關系[圖4(a)]可知，Xs2、Xs3、Xs4 和Xsm 與P+Xs1 震相之間的到時差與震中距無關，基本為一定值。它們的平均差值分別為7.8 s、21.6 s、61.1 s 和80.0 s?？梢姶舜蔚卣鸬恼郴e動過程主要由4 個子粘滑錯動事件構成。此次地震發(fā)生后，張勇等（2016）用遠場地震波形記錄反演了主震的震源破裂過程，得到了主震破裂過程中靜態(tài)滑動量分布隨時間變化的結果，見圖4(b)中a、b、c、d 矩形框。將本文給出的4 個子粘滑錯動事件的起始時間點與張勇等（2016）給出的靜態(tài)滑動量隨時間變化過程相對照可見，本文給出的粘滑錯動幅度隨時間的變化與張勇等（2016）給出的靜態(tài)滑動量分布隨時間的變化在幅度和時間點上比較吻合，說明粘滑錯動幅度隨時間的變化與彈性破裂后的靜態(tài)滑動量隨時間的變化是同步發(fā)生的。在0—7.8 s、7.8—21.6 s、21.6—61.1 s 和61.1—80.0 s 四個時段同步變化并逐漸增強，在P+Xs1 震相后約80.0 s 達到峰值。最終，與彈性破裂同步形成4 次同震粘滑錯動?？梢?，彈性破裂不是此次MW7.8 地震能量釋放的唯一形式，在彈性破裂釋放能量的同時，斷層粘滑錯動也在釋放能量。此次MW7.8 地震在彈性破裂開始后伴隨著4 次較大的粘滑錯動，是斷層粘滑錯動和彈性破裂共同作用的結果，是一種“粘滑錯動+彈性破裂”的機制。

圖4 時差值與震中距的關系（a）和靜態(tài)滑動量分布（b）Fig.4 The relationship between epicenter distance and different time(a) and the distribution of static slip quantity (b)

圖5 給出當震中距范圍為1.9°—144.5°時，（Pdif）+Xs1、X s 2、X s 3、X s 4、X s m和PP 震相的走時關系。P、PP 震相的走時由《IASEP1991 地震走時表 》（h＜35 km）（肯尼特，1993）給出，Pdif 震相走時由《Pdif 震相走時表》給出（許健生等，2013）。Xs1、Xs2、Xs3 和Xsm 震相的走時曲線基本平行，其走時與P 震相走時差是一定值，不隨震中距變化而變化，僅與粘滑錯動持續(xù)時間有關。

圖5 P(Pdif)+Xs1、Xs2、Xs3、Xs4、Xsm、和PP 震相的走時關系Fig.5 The relationship between the travel times of P(Pdif)+Xs1,Xs2,Xs3,Xs4,Xsm,and PP seismic phases

圖6 給出震中距1.9°—144.5°，P（Pdif）、PP、XsQ、XsR 和Rm 震相的走時關系。Rm震相走時由《JB 地震走時表》（Jeffreys et al，1958）給出，XsQ、XsR 是本文用全球地震臺網（GSN）100 個臺的甚寬頻帶（VH）數字波形記錄得出的結果。XsQ、XsR 的走時與震中距的擬合關系為

圖6 XsQ 和XsR 震相與P 震相的走時關系Fig.6 The relationship between the travel time of XsQ,XsR,and P seismic phases

式中，T的單位為s，Δ的單位為（°）。XsQ 和XsR 不僅是整個地震記錄中周期最長，振幅最大的震相，而且是面波波列中速度較快的震相，其速度明顯大于在地殼中傳播的面波Rm。

圖7 給出了震中距1.9°—144.5°，Xsm 與XsQ 和XsR 震相的走時差關系。Xsm 與XsQ 和XsR 震相的走時差在震中距為0 時，趨近于0。說明長周期的勒夫型面波XsQ 和瑞雷型面波XsR 應是粘滑錯動達到峰值時激發(fā)的面波。

圖7 Xsm 與XsQ、XsR 震相的走時差關系Fig.7 The relationship between the travel time differences of Xsm,XsQ,and XsR seismic phases

已有觀測結果表明，在2008 年5 月12 日中國汶川MW7.9 地震時，同震粘滑錯動在震中附近的ENH 地震臺（Δ=5.3°）激發(fā)出了周期約115.3 s 的XsQ 波動和周期約328.6 s 的XsR 波動；在2011 年3 月11 日日本MW9.1 地震時，同震粘滑錯動在震中附近的ERM 地震臺（Δ=3.8°）激發(fā)出周期約553.8 s 的XsQ 波動，在MAJO 地震臺（Δ=3.8°）激發(fā)出周期約448.8 s 的XsR 波動（許健生等，2019）。發(fā)生在內陸的中國汶川MW7.9 地震引發(fā)了大面積山體滑坡、泥石流等次生災害；發(fā)生在海里的日本MW9.1 地震引發(fā)了罕見的海嘯等次生災害；同樣，在本次新西蘭南島近海岸MW7.8 地震時，同震粘滑錯動在震中附近的SNZO 地震臺（Δ=1.9°）激發(fā)出周期約600 s 的XsQ 和XsR 波動，引發(fā)浪高達3 m 的海嘯等次生災害。據此推測，這種周期長達數百秒的同震粘滑錯動容易引發(fā)山體滑坡、泥石流和海嘯等次生災害。該新發(fā)現和初步推測尚需更多震例加以驗證。若證實這種規(guī)律的存在，則有助于在震后短時間內預判地震次生災害。

5 結論

根據上述理論研究、實驗和觀測結果，得出以下幾點認識：①P+Xs1、Xs2、Xs3 和Xs4 震相是彈性破裂開始后的同震粘滑錯動震相，Xsm 是粘滑錯動的峰值震相，標志著粘滑錯動過程中能量釋放的峰值點；②XsQ 和XsR 震相是粘滑錯動達到峰值后，止滑過程激發(fā)的長周期面波震相；③彈性破裂不是MW7.8 地震能量釋放的唯一形式，彈性破裂開始后有4 次粘滑錯動，MW7.8 地震是粘滑錯動和彈性破裂共同作用的結果，是一種“粘滑錯動+彈性破裂”的機制。

從地震記錄圖上識別出Xs 震相有助于認識同震粘滑錯動過程，識別出XsQ 和XsR 震相有助于預判地震次生災害。

數據轉換和震相標注均由黑龍江省地震局和躍時高級工程師提供的CDSN 臺站分析軟件完成，在此表示感謝。