蘇 珊，郭培蘭

（廣西壯族自治區(qū)地震局，南寧 530000）

0 引言

2008 年至2021 年，廣西及鄰區(qū)先后發(fā)生平果M4.5、蒼梧M5.4、北流M5.4 和靖西M5.2 等13 次4級以上地震，地震活動顯著增強，多數(shù)地震沿桂西地區(qū)的北西—南東向斷裂呈帶狀分布，桂東南的北東—南西向斷裂發(fā)生了兩次5 級以上地震（圖1）。廣西及鄰區(qū)地處地震活躍的川滇地區(qū)與相對穩(wěn)定的華南地塊過渡地帶，區(qū)內(nèi)地震分布具有明顯的分區(qū)特征，整體為西強東弱。既往研究表明，區(qū)域應(yīng)力環(huán)境與地震發(fā)生機理和孕震過程等密切相關(guān)，其變化趨勢是地震危險性研判的重要指標，亟待進一步研究。

圖1 2008年以來廣西及鄰區(qū)M4級以上地震分布Fig.1 Distribution of earthquakes withM≥4 in Guangxi and its adjacent areas since 2008

1944 年 Gutenberg 和 Richter[1]提 出 的震 級 — 頻度公式（G-R關(guān)系）lgN=a-b×M，其中b值時地震學中最重要的地震活動性參數(shù)之一。實驗室?guī)r石聲發(fā)射微震實驗結(jié)果表明，應(yīng)力與b值之間存在一個負相關(guān)關(guān)系，當應(yīng)力增強時，b值降低，而當應(yīng)力減弱時，b值升高，這一關(guān)系在真實的地震活動中也得到體現(xiàn)，b值也因此被作為地下應(yīng)力指示計[2-6]。大震分發(fā)生與地下應(yīng)力的集中增強密切相關(guān)，而b值與應(yīng)力的負相關(guān)關(guān)系使它的時空變化特征對于地震危險性評估和地震預(yù)報有著重要意義。在實際應(yīng)用中，b值的時空變化與大震的發(fā)生關(guān)系往往更為復雜。對川滇1981年以來19次M6.5地震與b值的時空分布關(guān)系研究結(jié)果表明，川滇地區(qū)b值空間分布特征與強震活動之間聯(lián)系并不顯著，但震前b值的時間分布呈現(xiàn)系統(tǒng)地下降，58%的震源區(qū)的b值在強震前出現(xiàn)顯著下降趨勢[7]。陳學忠等[8-9]對 1976 年唐山M7.8 和 2019 年長寧M7.8 地震前b值時空變化特征進行研究，結(jié)果表明這兩次地震前出現(xiàn)顯著的b值下降現(xiàn)象。史海霞等[10]在對2008 年汶川M8.0 地震前b值時空變化特征的研究中發(fā)現(xiàn)，震前b值從2002 年始至地震前呈現(xiàn)出一個長期趨勢性減小，而后出現(xiàn)持續(xù)約兩年的穩(wěn)定并略有上升趨勢，在地震前約半年，出現(xiàn)快速、顯著的下降，這個結(jié)果與與Nanjo 等[11]報道的日本東北9.0 級地震前b變化特征非常相似。本文將系統(tǒng)的研究2008 年以來廣西及鄰區(qū)b值的時空變化，探索其與區(qū)域中強震孕育之間的聯(lián)系，為地震風險評估提供參考。

1 資料和研究方法

1.1 研究資料的選取

本文將廣西及鄰區(qū)（104°～113°E，20°～27°N）作為研究區(qū)，廣西地區(qū)2008 年1 月廣西數(shù)字測震臺網(wǎng)竣工投入正式運行[12]，區(qū)域監(jiān)測能力得到較大提升，因此選取研究區(qū)內(nèi)2008年1月至2022年3月地震目錄作為計算b值的基礎(chǔ)數(shù)據(jù)。根據(jù)在此期間研究區(qū)域內(nèi)的的地震震級—頻度（G-R關(guān)系）分布圖（圖2），震級在ML1.5～4.3 的范圍內(nèi)G-R關(guān)系線性擬合效果最好，確定最小完備性震級MC=ML1.5，震級上限為ML4.3。

圖2 2008年1月至2022年3月廣西及鄰區(qū)地震震級-頻度（G-R）關(guān)系Fig.2 The G-R relationship between magnitude and frequency of earthquakes in Guangxi and its adjacent areas from January 2008 to March 2022

1.2 b值分析方法

本文采用基于格點搜索的b值分析方法程序?qū)ρ芯繀^(qū)進行b值時空分析[13]。該程序選用最大似然法[14]來計算b值：

b值95%的置信概率為：

在計算b值隨時間的變化時，采用基于地震事件數(shù)量選取窗長和步長的方法，即選取等事件數(shù)量作為計算b值的窗長，并將窗口以等事件數(shù)量的增量進行滑動，使得各個時間節(jié)點上用于計算b值的地震事件數(shù)量是均勻的，確保計算結(jié)果不受樣本量變化的影響。

計算b值在空間上的變化特征時，首先將研究區(qū)劃分為一定大小的格點，以格點為中心選取固定大小的正方形區(qū)域作為它的空間窗，計算每個空間窗內(nèi)b值隨時間的變化曲線，而后，挑選出研究區(qū)內(nèi)b值顯著下降的格點，并將其作為異常點，多個異常點集中的區(qū)域作為b值異常區(qū)。選取這個異常區(qū)內(nèi)的地震計算b值隨時間的變化，從而可以得到整個b值異常區(qū)的范圍及其b值隨時間變化的特征。本文將研究區(qū)劃分為0.1°×0.1°的格點，選取的空間窗大小為2°×2°。

對于b值的時間變化曲線b（t），b值的相對變化表示為：

其中b1是時間t1的b值，而b2是時間t2的。通過計算各節(jié)點b值的相對變化，可以得到b值相對變化的空間分布。

2 結(jié)果

本文將研究區(qū)劃分為0.1°×0.1°的格點，選取的空間窗大小為2°×2°，每個格點以100個地震計算b值，每3個地震向前滑動計算時間變化曲線。Lei[15]研究發(fā)現(xiàn)，統(tǒng)計窗口中事件數(shù)量越大，b值曲線越平滑，事件數(shù)量越少，b值演化的分辨率越高。經(jīng)過試算后，最終確定100 個地震事件為窗長、3 個事件為步長，在保證b值變化信息不丟失的情況下盡量平滑曲線，b值置信概率為95%。在計算過程中，如果格點內(nèi)得到b值數(shù)量較少，則其時間曲線上節(jié)點的均勻性可能較低，無法反應(yīng)出可靠的b值變化信息[8]。我們設(shè)定時間曲線上最少數(shù)據(jù)點數(shù)為10 個，最終共得到研究區(qū)內(nèi)4014 個格點的b值隨時間變化曲線。圖3 展示了2008 年以來研究區(qū)內(nèi)b值變化的空間分布情況：廣西中西部地區(qū)存在一個顯著的b值下降區(qū)域，覆蓋靖西—崇左、百色—合浦、巴馬—博白三條北西向斷裂和桂林—南寧、防城—靈山兩條南東向斷裂，也是這五條斷裂的交匯區(qū)域。該b值下降區(qū)域先后發(fā)生了5 次M4 以上地震，包括 2019 年 11 月 25 日靖西M5.2 和 2021年8月4日德保M4.8地震；位于廣西東南部粵桂交界區(qū)域的合浦—北流斷裂附近存在一個b值上升區(qū)域，于2019年10月12日發(fā)生北流M5.2地震。

圖3 2008-01—2022-03廣西地區(qū)b值空間分布Fig.3 Spatial distribution of b value in Guangxi from January 2008 to March 2022

為分析研究區(qū)內(nèi)b值顯著變化區(qū)域的時間曲線，本文篩選了b值下降超過10%和上升超過30%的區(qū)域（圖4（a），（c）），分別計算b值在這兩個區(qū)域的時間變化情況。受區(qū)域范圍和地震數(shù)量的影響，b值下降區(qū)域以300 個地震為時間窗10 個地震滑動計算b值（圖4（b）），上升區(qū)域以100 個地震為時間窗5 個地震滑動計算b值（圖4（d））。最終結(jié)果顯示，b值下降區(qū)域在2019 年靖西M5.2 地震之前，存在一個近4 年的長時間b值下降（圖4（d）方框內(nèi)），2014 年 11 月至 2018 年 8 月b值從 1.17 下降到了0.8，而后b值出現(xiàn)顯著的回升，持續(xù)六個月至 2019 年 4 月，2019 年 5 月開始b值再次下降直至發(fā)震。靖西5.2 級地震后，b值短暫回升轉(zhuǎn)平之后發(fā)生2021 年德保M4.8 地震。在粵桂交界地區(qū)，2019年北流5.2級地震前b值處于從低值上升階段，震后而后出現(xiàn)短暫的b值上升速率增加，約一年后恢復震前上升速率。

圖4 2008-01—2022-03廣西及鄰區(qū)b值空間分布和時間變化情況Fig.4 Spatial distribution and temporal variation of b value in Guangxi and its adjacent areas from January 2008 to March 2022

3 討論

孟昭彤等[16]根據(jù)以往震列總結(jié)出了4 種大震前b值隨時間變化的類型：①在大地震發(fā)生之前的幾年時間，b值升高，主震發(fā)生之后，b值下降；②在大地震發(fā)生之前幾年，b值降低，主震發(fā)生之后，b值升高；③大震發(fā)生前b值會明顯下降，并在回升過程中發(fā)生大震；④在大地震發(fā)生前，b值升高，但在主震前短時間內(nèi)有下降的前兆。顯然圖4（b）中展示的b值變化情況與以上4 種類型接皆不相符，靖西地震發(fā)生在長時間b值下降導致的低b值背景下，發(fā)震前出現(xiàn)顯著的回升后快速降低轉(zhuǎn)折。這種b值變化特征與2008 年汶川8.0 級地震[10]和2011 年日本東北9.0 級地震非常相似。一般認為，b值與應(yīng)力呈反比，地下高應(yīng)力狀態(tài)更有利于地震的發(fā)生，大震之前低b值的特性已經(jīng)得到普遍驗證[17-20]。廣西中西部地區(qū)長時間的b值下降反映了區(qū)域內(nèi)應(yīng)力的持續(xù)性增強，而GNSS 基線結(jié)果顯示，華南塊體西部2011 至今持續(xù)受到川滇地塊的順時針壓扭作用，為區(qū)域應(yīng)力增強提供了動力來源[21]。馬瑾等[22]通過對地震前亞失穩(wěn)狀態(tài)的研究，將粘滑型地震孕育分為平穩(wěn)變形，偏離線性，亞失穩(wěn)，失穩(wěn)和失穩(wěn)后5個階段，其中偏離線性階段是應(yīng)力釋放的開始，在此階段處于由應(yīng)力積累向應(yīng)力釋放的過渡階段，應(yīng)力釋放尚不占優(yōu)勢，亞失穩(wěn)階段應(yīng)力釋放逐步占據(jù)優(yōu)勢。這或許可以解釋靖西地震前一年左右出現(xiàn)的b值上升再下降的擾動變化現(xiàn)象，即b值上升意味著地震孕育處于偏離線性階段，在這個階段應(yīng)力累積已經(jīng)結(jié)束而釋放并未完全開始，造成了階段性的區(qū)域地震活動性減弱，b值上升。而后地震孕育進入亞失穩(wěn)階段，應(yīng)力釋放逐漸顯著，b值也同步下降，直至失穩(wěn)，地震發(fā)生。圖4（d）中展示的b值上升區(qū)于2019 年10月12日發(fā)生北流5.2級地震，地震前后，除了震后余震造成的相對短時間的b值擾動，其上升的整體態(tài)勢不變，這不同于研究領(lǐng)域?qū)τ谡鹎癰值下降的一般認識。1996 年麗江M7 地震震前同樣具有高b值動態(tài)演化特點，同時震中及震后余震區(qū)位于高波速比異常區(qū)[23]，其所反映出的孕震模式被稱為堅固體孕震模式，即非均勻地殼中高速擾動體（堅固體）的存在是局部地區(qū)應(yīng)力高度集中的根本條件，當高速擾動體被低速擾動區(qū)包圍時，即震源體的剛度比四周的環(huán)境剛度大時可以滿足失穩(wěn)條件導致突然破裂發(fā)生大地震；否則可以導致緩慢破裂而不發(fā)生大地震[24]，區(qū)域出現(xiàn)高b值特征。黃強強等[25]針對北流地震震源區(qū)短周期密集臺陣背景噪聲成像的研究結(jié)果表明，震源區(qū)存在NNWSSE 向的高速異常，而震中上方為上凸形態(tài)的高速異常，這些高速異?？赡転橹猩◢弾r基，這種局部異常高速體可能是華南“馬賽克”結(jié)構(gòu)中微陸塊在上地殼中的延伸，上地殼結(jié)構(gòu)的不均勻性和大尺度構(gòu)造相互作用，會在局部區(qū)域?qū)е聭?yīng)力積累，從而引發(fā)中強地震，這種地震成因機理與堅固體孕震模式非常類似，推測為北流地震震源區(qū)及周邊b值持續(xù)上升的根本原因。

事實上靖西5.2 級地震前b值變化特征并不孤立，考慮b值地下應(yīng)力計的性質(zhì)，其所反應(yīng)的區(qū)域應(yīng)力累積、釋放等與地震孕震相關(guān)內(nèi)在機制還有待深入研究。但不能否認，在b值長期下降的背景下其值出現(xiàn)先上升后下降的變化形態(tài)是具有明顯的前兆意義的，需要持續(xù)關(guān)注，同時這種b值變化特征對于中強地震的普適性還需要深入研究。北流地震前后區(qū)域內(nèi)b值始終升高的現(xiàn)象，與地震孕震模式密切相關(guān)，震源區(qū)地下高速異常表明其具有明顯的堅固體孕震特征。

4 結(jié)論

本文基于廣西臺網(wǎng)提供的地震目錄，利用基于格點搜索的b值分析方法程序，對研究區(qū)進行b值時空分析，主要認識如下：

（1）華南西部存在兩個b值顯著變化區(qū)域，分別為廣西中西部b值顯著下降區(qū)和粵桂交界b值顯著上升區(qū)，這兩個區(qū)域內(nèi)分別發(fā)生的2019 年靖西5.2及地震和2019年北流5.2級地震。

（2）廣西中西部地區(qū)在長時間b值下降導致的低b值背景下，出現(xiàn)顯著的回升后快速降低轉(zhuǎn)折，而后發(fā)生靖西5.2 及地震。這種b值變化特征與2008 年汶川 8.0 級地震和 2001 年日本東北 9.0 級地震震前區(qū)域b值變化非常相似，具有明顯前兆意義。

（3）粵桂交界區(qū)域 2019 年 10 月發(fā)生北流 5.2 級地震，發(fā)震前后b值上升態(tài)勢不變，結(jié)合地下速度結(jié)構(gòu)分析，認為其具有明顯的堅固體孕震特征，從而導致b值升高異常。