鄭永通，全建軍，龔 薇，劉愫昀，王紹然，賴福霖

（1. 福建省地震局 龍巖地震監(jiān)測中心站, 福建 龍巖 364000；2. 福建省地震局 南平地震監(jiān)測中心站, 福建 南平 353000）

0 引言

b值是用來衡量某地區(qū)地震活動水平的一種標(biāo)志，是反映一個地區(qū)承受平均應(yīng)力和接近強度極限程度的一個參數(shù)[1]。高國英認(rèn)為大震前b值異常時間長短、異常量值大小與未來強震震級有關(guān)[2]；王筱榮等認(rèn)為b值判定預(yù)備震源1年內(nèi)是否發(fā)震正確率為75%，2年內(nèi)發(fā)震正確率為92%[3]；鄭福婉等認(rèn)為震前b值存在著不同程度的異常，其變化方式有兩種：一種是震前升高而后下降發(fā)震，一種是震前b值下降后回升時發(fā)震，多數(shù)震例屬于后者[4]。對于某個區(qū)域內(nèi)已發(fā)生的地震，通過b值空間掃描，分析該區(qū)域低b值或高b值異常，來判斷未來該區(qū)域是否存在發(fā)生強震或中強地震的可能提供依據(jù)。本文著重研究福建及其鄰近沿海地區(qū)發(fā)生中強地震前的b值空間掃描圖像特征，探討b值空間掃描異常區(qū)域與中強震之間的關(guān)系，為預(yù)測地震提供參考借鑒。

1 福建內(nèi)陸地質(zhì)構(gòu)造和地震活動概況

福建內(nèi)陸主要有3條NE-NNE向地震構(gòu)造帶，分別是：邵武-河源斷裂帶、政和-海豐斷裂帶、長樂-詔安斷裂帶；NW-NNW向主要有4條斷裂帶，分別是：上杭-云霄斷裂帶、九龍江斷裂帶、永安-晉江斷裂帶、閩江斷裂帶，兩組斷裂帶構(gòu)成福建內(nèi)陸縱橫交錯的地震地質(zhì)構(gòu)造背景。據(jù)史料記載，歷史上福建內(nèi)陸僅在泉州發(fā)生一次8.0級大地震（1604年），此次地震破壞較重，分布范圍廣，有感范圍最遠(yuǎn)達(dá)1 000 km。近幾十年來，福建內(nèi)陸中強地震時有發(fā)生，是東南沿海中強震較為活躍的地區(qū)。

福建內(nèi)陸自1990年至今共發(fā)生11次ML≥4.5地震，其中最大地震為1997年5月31日發(fā)生在永安小陶的ML5.3地震。從地震發(fā)生時所處的斷裂帶來看，中強震主要發(fā)生在政和-海豐斷裂帶，共發(fā)生7次ML≥4.5地震；其次是長樂-詔安斷裂帶，共發(fā)生4次ML≥4.5地震（圖1）。

圖1 福建內(nèi)陸主要活動斷裂帶及地震活動特點

2 b值計算

2.1 b值的定義

Gutenberg et al[5]在研究世界地震活動時，根據(jù)全球各大地震區(qū)6級以上地震數(shù)目統(tǒng)計結(jié)果，提出地震震級與頻度的關(guān)系：

式（1）亦稱古登堡公式或G-R關(guān)系式。式中：N為震級M以上地震的頻度；a和b為統(tǒng)計常數(shù)，a表示區(qū)域內(nèi)的地震活動水平，b值代表一個區(qū)域內(nèi)不同大小地震頻度的比例關(guān)系。b值的計算主要有線性最小二乘法、極大似然法、非線性最小二乘法、矩估計方法等4種。就計算誤差而言，極大似然法最小[6]。因此，本文b值計算方法選取極大似然法。

2.2 起算震級的選取

對范圍較大的空間區(qū)域進(jìn)行b值空間掃描計算時，應(yīng)根據(jù)臺網(wǎng)監(jiān)測能力，并參考震級-頻度關(guān)系式確定整個研究區(qū)可使用的震級下限[7]。由于1998年之前，測震記錄信號為模擬信號，地震監(jiān)測能力比較弱，根據(jù)LgN-M圖的擬合結(jié)果來看（圖2a），起算震級可選為ML1.9。實行數(shù)字化改造后，地震監(jiān)測能力有一定的提高。根據(jù)LgN-M圖的擬合結(jié)果來看（圖2b），起算震級可選為ML1.7。綜合以上分析結(jié)果認(rèn)為，起始震級應(yīng)選取ML≥1.9。本文選取ML=2.0為起算震級，對福建內(nèi)陸進(jìn)行b值空間掃描。

圖2 福建內(nèi)陸最小完整性震級的確定

2.3 參數(shù)設(shè)定

本文研究區(qū)域為福建內(nèi)陸（115°～121°E，23°～29°N），地震目錄選取福建省地震局監(jiān)測中心提供的1970年1月1日至2021年4月30日ML≥2.0地震目錄，窗長設(shè)為12個月，時間步長設(shè)為1個月；緯度窗長為1.0°、步長為0.5°；經(jīng)度窗長為1.0°、步長為0.5°，分別對福建內(nèi)陸自1970年1月至今發(fā)生的11次中強震震前作b值空間掃描圖像。

2.4 福建內(nèi)陸中強震前b值空間掃描分析

從福建內(nèi)陸1990年以來ML≥4.5地震（表1）可知，福建內(nèi)陸中強震主要集中在1992—1999年、2007—2008年、2013年3個時段。其中1992—1999年共發(fā)生中強地震4次，2007—2008年共發(fā)生中強地震4次，2013年共發(fā)生中強地震3次。下面分別對以上幾個時間段發(fā)生的中強震震前進(jìn)行b值空間掃描分析，探討中強震前b值空間掃描異常特征。

表1 1990年以來ML≥4.5地震

2.4.1 1992年11月26日龍巖5.1級地震前b值空間掃描分析

從b值空間掃描結(jié)果來看，1992年11月26日龍巖發(fā)生ML5.1地震前，從1990年9月開始在震中及其附近出現(xiàn)b值異常，異常值為2.5；之后b值開始下降至1.5左右，異常持續(xù)到1991年10月后結(jié)束，異常持續(xù)時間為14個月；在異常結(jié)束12個月后發(fā)生地震（圖3）。

圖3 1990年9月至1992年11月b值空間掃描圖像

2.4.2 1994年5月24日龍巖4.8級地震前b值空間掃描分析

1994年5月24日龍巖4.8級地震前，從1993年1月開始，在震中附近出現(xiàn)b值異常，異常值為1.0，1993年5月異常值開始上升，至1993年6月，震中區(qū)b值異常值達(dá)到最大為1.4；之后開始下降，1993年7月，異常值下降至1.0；1993年8月開始，震中附近異常值消失，在異常結(jié)束9個月后發(fā)震（圖4）。

圖4 1992年12月至1994年5月b值空間掃描圖像

2.4.3 1997年5月31日永安ML5.3地震前b值空間掃描分析

1997年5月31日永安ML5.3地震前，從1996年4月開始在震中附近出現(xiàn)b值異常，異常值為1.2左右；從1997年1月開始，b值出現(xiàn)下降，1997年2月b值下降至0.6，持續(xù)到發(fā)震，異常持續(xù)時間為14個月（圖5）。

圖5 1996 年 5 月至 1997 年 5 月b 值空間掃描圖像

2.4.4 1999年2月13日永安ML4.6地震前b值空間掃描分析

1999年2月13日永安ML4.6地震前，從1997年6月開始在震中附近出現(xiàn)b值空間掃描異常，異常值為0.6左右；從1998年3月開始b值異常值上升至0.9；持續(xù)至1999年2月異常結(jié)束，異常持續(xù)時間為21個月（圖6）。

圖6 1997年6月至1999年2月b值空間掃描圖像。

2.4.5 2007—2008年，永春、順昌、古田、長泰中強震前b值空間掃描分析

2007年8月29日永春ML4.6地震前，從2006年5月開始在震中附近出現(xiàn)b值異常，異常值為0.9；之后異常開始上升，2007年5月b值達(dá)到最大為1.5；2007年6月b值下降，震前震中異常值為0.6。2008年3月6日，古田ML4.8地震前，從2007年10月開始，在震中附近開始出現(xiàn)b值異常，異常值為0.5；異常一直持續(xù)到發(fā)震，異常持續(xù)時間為6個月，震前震中b值異常值為0.5～0.6。2008年7月5日長泰ML4.7地震前，從2007年8月開始在震中附近出現(xiàn)b值異常，異常值為0.4；之后b值開始緩慢上升，2007年12月之后，異常值穩(wěn)定在0.8～0.9之間；一直持續(xù)到發(fā)震，異常持續(xù)時間為12個月（圖7）。順昌地震震前未出現(xiàn)b值異常。

圖7 2006年5月至2008年7月b值空間掃描圖像

2.4.6 2013年仙游幾次中強震前b值空間掃描分析

2013年仙游發(fā)生3次中強地震前，從2011年12月開始在震中附近開始出現(xiàn)b值異常，異常值為0.8；之后b值有升有降，異常值范圍在0.6～1.0之間。在仙游陸續(xù)發(fā)生3次中強震期間，異常未結(jié)束（圖8）。2013年8—10月在仙游連續(xù)發(fā)生3次中強震期間，b值變小，異常值僅為0.3～0.4。

圖8 2011年12月至2013年10月b值空間掃描圖像

從福建內(nèi)陸中強震前b值異常情況來看（表2），除2007年3月13日順昌ML4.9地震震前未出現(xiàn)b值空間掃描異常外，其余中強地震震前在震中或附近均出現(xiàn)b值異常，且多數(shù)中強震發(fā)生后，b值異常仍然存在。另外， 2013年8月—10月仙游連續(xù)發(fā)生3次中強地震前，異常持續(xù)時間也相對較長。

表2 福建內(nèi)陸中強震前b值異常情況表

3 結(jié)論

通過對1990年至今發(fā)生在福建內(nèi)陸的中強震震前b值空間掃描分析，得出以下結(jié)論。

1）不同的起始震級計算出的b值結(jié)果差別較大，在計算b值時，要考慮研究區(qū)域的臺網(wǎng)監(jiān)測能力，選取適合本區(qū)域的起算震級來計算b值。

2）1990年至今，發(fā)生在福建內(nèi)陸的中強震震前絕大多數(shù)地震存在b值空間掃描異常，且大部分異常一直持續(xù)到發(fā)震仍然存在。震前b值異常變化有3種形式，即上升、下降或上升-下降交替變化。異常持續(xù)時間往往與異常范圍1～2年內(nèi)發(fā)生中強震的次數(shù)有關(guān)，若異常范圍內(nèi)發(fā)生的中強震次數(shù)越多，則異常持續(xù)時間也越長。

3）當(dāng)某一區(qū)域出現(xiàn)b值空間掃描異常時，要注意該區(qū)有可能發(fā)生中強地震。由于有些中強震發(fā)生前b值異常已結(jié)束，有些中強震發(fā)生前b值異常仍然存在，所以對發(fā)震時間的預(yù)測，要參考其他前兆儀器的觀測資料分析結(jié)果來進(jìn)行綜合判定。