全建軍 ， 賴見深 ， 陳珊樺 ， 鄭永通 ， 陳美梅 ， 林慧卿

（1.福建省地震局永安地震臺(tái)， 福建 永安 366000； 2.福建省地震局東山地震臺(tái)， 福建 邵武 354000；3.福建省地震局泉州地震臺(tái)， 福建 泉州 364000； 4.福建省地震局龍巖地震臺(tái)， 福建 龍巖 363000）

0 引言

地震是地殼長時(shí)間構(gòu)造運(yùn)動(dòng)的生成物， 當(dāng)?shù)卣鹪杏l(fā)展時(shí)， 地殼發(fā)生變形是最直接、 最密切的伴隨現(xiàn)象， 大地出現(xiàn)變形是構(gòu)造運(yùn)動(dòng)在地表的反應(yīng)。 由于地震波的刺激， 地震發(fā)生過程中， 在觀測(cè)點(diǎn)通過數(shù)字化前兆形變監(jiān)測(cè)儀器， 能觀測(cè)到傾斜、 應(yīng)變及應(yīng)力變化情況， 這就是形變波， 用以區(qū)分地震波[1]。 地震引起的共震形變波包含地震孕育以及發(fā)動(dòng)過程中的多數(shù)信息， 利用這些同震信息能估算部分地震參數(shù)值， 通過共震形變波物理方面的深入研究， 能更好的了解短臨地震前兆的動(dòng)力學(xué)特點(diǎn)。

數(shù)字化體應(yīng)變前兆觀測(cè)使同震形變記錄成為可能， 小陶臺(tái)采用的是TJ-Ⅱ型鉆孔體應(yīng)變儀， 該儀器被用來監(jiān)測(cè)地殼應(yīng)變信息以及探索地球物理學(xué)與地球動(dòng)力學(xué)。 鉆孔體應(yīng)變儀一般固定在深井中， 其通過專用的水泥當(dāng)做耦合材料， 具備很好的運(yùn)行平穩(wěn)性以及監(jiān)測(cè)準(zhǔn)確性[2]。 它可以清楚地觀測(cè)到固體潮汐改變， 也可觀測(cè)到地震發(fā)生過程中地殼出現(xiàn)的形變， 尤其適用于監(jiān)測(cè)固體潮汐以及地殼形變的持續(xù)改變， 可為研究地震發(fā)展過程中地殼形變改變方式收集資料， 也可為研究地球彈性變化給予關(guān)鍵參數(shù)依據(jù)[3]。 利用TJ-Ⅱ鉆孔體應(yīng)變儀對(duì)地球內(nèi)部的形變變化實(shí)施高精準(zhǔn)監(jiān)測(cè)， 進(jìn)而了解地殼前兆形變的長、 中、 短期、 臨震與震后改變的時(shí)空排列和調(diào)整過程， 可為研究地震學(xué)與地區(qū)防災(zāi)項(xiàng)目的開展給予基礎(chǔ)性分析資料。

在地震前兆異常探索尚未成熟的今天， 前兆觀測(cè)設(shè)備同震與映震能力的探索已得到眾多科研人員重視。 如陸明勇、 呂品姬、 關(guān)華、 焦成麗等[4-7]， 均對(duì)形變監(jiān)測(cè)設(shè)備的映震強(qiáng)弱做了深度分析， 特別是北京殼力研究所匯編的《應(yīng)變儀同震變化研究分析報(bào)告》， 它組織全國鉆孔應(yīng)變觀測(cè)站點(diǎn)， 對(duì)形變?cè)O(shè)備同震與映震強(qiáng)弱進(jìn)行較完整的歸納和比較分析， 取得了顯著效果。 本文對(duì)小陶臺(tái)鉆孔體應(yīng)變儀自投入運(yùn)行以來的觀測(cè)資料進(jìn)行研究； 統(tǒng)計(jì)分析所記錄的地震， 總結(jié)出同震形變波形特征； 并通過震級(jí)與震中距關(guān)系、 最大震幅與震級(jí)關(guān)系兩方面的研究， 量化評(píng)估小陶臺(tái)鉆孔體應(yīng)變儀的映震能力， 可為形變臺(tái)站場(chǎng)地監(jiān)測(cè)效能以及儀器布設(shè)密度提供依據(jù)。

1 井孔概況及儀器安裝

1.1 井孔概況

小陶臺(tái)鉆井孔臺(tái)址位于永安市小陶鎮(zhèn)第二中學(xué)內(nèi)， 117.131°E， 25.729°N， 高程： 65 m， 東邊是新華夏系政和—海豐斷裂帶， 西邊是邵武—河源斷裂帶， 南北向各是NW 方向的弱一等斷裂：永安—晉江以及九龍江斷裂所控制（圖1）。 井孔所處地表為第四系土層所覆蓋， 厚度約10 m， 表現(xiàn)為軟塑—可塑狀態(tài)。 第四系覆蓋層下分布地層主要巖性為花崗巖， 風(fēng)化較弱， 巖石強(qiáng)度大、 結(jié)構(gòu)致密、 巖石顆粒硬度大， 局部巖石硅化富含石英，具備強(qiáng)研磨性， 最高研磨性能達(dá)10 級(jí)。 鉆井頂層巖土的能鉆程度是Ⅰ至Ⅵ度， 下層能鉆程度是Ⅸ至Ⅺ度。 小陶鉆孔體應(yīng)變觀測(cè)井處在小陶中學(xué)西邊10～12 m， 觀測(cè)井南邊200 m 位置有一小溪流過， 該溪流速緩慢且水位常年較低， 即使降雨頻繁的4 至6 月， 也未看到溪流給小陶臺(tái)形變監(jiān)測(cè)帶來影響。

圖1 福建小陶臺(tái)鉆孔體應(yīng)變所處地質(zhì)斷層分布圖Fig.1 Geological fault distribution of the borehole body strain at Xiaotao Station

2 鉆孔體應(yīng)變儀觀測(cè)資料質(zhì)量分析

地震前兆觀測(cè)中， 判斷觀測(cè)資料質(zhì)量的關(guān)鍵依據(jù)主要有三點(diǎn)， 分別是完整性、 穩(wěn)定性以及可靠性。 只有前兆資料達(dá)到相關(guān)監(jiān)測(cè)規(guī)范標(biāo)準(zhǔn)后，該資料才具有真正的使用意義， 才可用來判斷前兆資料的日常動(dòng)態(tài)變化， 才可正確反映前兆監(jiān)測(cè)點(diǎn)的受力狀態(tài)。 因此， 對(duì)小陶臺(tái)鉆孔體應(yīng)變儀前兆數(shù)據(jù)的完整性、 穩(wěn)定性及可靠性進(jìn)行分析， 非常有必要。

2.1 資料完整性分析

小陶地震臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀正式運(yùn)行后，日觀測(cè)曲線平滑無頓挫， 正式觀測(cè)的15 d 時(shí)間里，未出現(xiàn)開閥情況。 觀測(cè)儀器所監(jiān)測(cè)的固體潮變化完整清晰， 日變化觀測(cè)曲線所包含的噪聲較小，每日的周期改變性較顯著， 具備較好的規(guī)律特點(diǎn)，與理論固體潮表現(xiàn)形式擁有較高的同步性 （見圖2）。 通過 （1）式， 對(duì)小陶臺(tái)鉆孔體應(yīng)變儀 2016 年 1月至2020 年5 月觀測(cè)資料的完整率進(jìn)行統(tǒng)計(jì)計(jì)算， 其數(shù)據(jù)完整率達(dá)99.95%。

2.2 資料穩(wěn)定性分析

通過觀測(cè)數(shù)據(jù)的年變幅度與年零漂兩個(gè)指標(biāo)，可以對(duì)觀測(cè)儀器的穩(wěn)定性進(jìn)行檢測(cè)。 由于鉆孔體應(yīng)變儀的傳感探頭埋深較深， 而且加入的膨脹水泥已和井中的基巖耦合成一塊。 所以從理論上分析， 儀器探頭不會(huì)因地表溫度改變而影響到年變形態(tài)， 正常年變幅度會(huì)位于一個(gè)波動(dòng)較小且較穩(wěn)定的范圍之內(nèi)[8]。 日常觀測(cè)中， 判斷監(jiān)測(cè)設(shè)備和墩基是否具有較高穩(wěn)定性， 一般是利用數(shù)據(jù)的年度零漂值， 如果觀測(cè)設(shè)備的年度漂移較低， 墩基通常比較穩(wěn)固， 而且具有較好的耦合性[9]， 觀測(cè)到的前兆資料準(zhǔn)確性也較高。 圖3 是小陶臺(tái)鉆孔體應(yīng)變儀2016 年1 月至9 月的整點(diǎn)值變化曲線圖， 從圖3 可以看出小陶臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀觀測(cè)值月潮汐形態(tài)清晰， 一階差分去除趨勢(shì)變化后， 大小潮清晰明顯， 月變幅在100E-9以內(nèi)， 趨于穩(wěn)定。

圖2 福建小陶臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀與固體潮理論值對(duì)比曲線圖Fig.2 Comparison curve between theoretical values of TJ- II borehole body strain gauge and solid tide at Xiaotao station

圖3 福建小陶臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀月整點(diǎn)值（a）及一階差分曲線圖（b）Fig.3 Monthly point value（a）and the first-order difference curve（b）of TJ-Ⅱ borehole body strain gauge at Xiaotao Station

2.3 觀測(cè)資料可靠性分析

對(duì)固體潮汐形變進(jìn)行分析、 檢測(cè)觀測(cè)儀器與觀測(cè)資料的內(nèi)在精度， 通常使用Venedikov 調(diào)和分析法[10]。 本文使用 Venedikov 調(diào)和分析法， 對(duì)小陶臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀 2016 年1 月至 2020 年 5 月的觀測(cè)資料進(jìn)行逐年分析， 并通過分析結(jié)果來追蹤潮汐參數(shù)伴隨時(shí)間的改變。 表1 列出了2016—2020 年按照年序列的潮汐參數(shù)， 可以看出小陶臺(tái)觀測(cè)資料滿足形變學(xué)科提出的 α≥2.0 和 δα/α≤0.05 的內(nèi)在精準(zhǔn)度指標(biāo)。

3 地震波記錄及同震變化分析

地殼在長期構(gòu)造運(yùn)動(dòng)中發(fā)生形變， 不斷累積應(yīng)力， 當(dāng)應(yīng)力超過巖石承受極限時(shí)， 發(fā)生斷裂，并產(chǎn)生錯(cuò)動(dòng)， 應(yīng)力迅速釋放形成構(gòu)造地震。 因此，構(gòu)造應(yīng)力與斷層滑動(dòng)互成因果， 前者的超限促進(jìn)后者加速， 后者反過來釋放前者， 使地殼巖層受力重新達(dá)到平衡穩(wěn)定。 TJ-Ⅱ型鉆孔體應(yīng)變儀不僅具有較好高頻特性與耐振動(dòng)性， 而且能將動(dòng)態(tài)消息真實(shí)記錄下， 它可以較好地感應(yīng)到應(yīng)力觸發(fā)時(shí)產(chǎn)生的應(yīng)變波動(dòng)。 用鉆孔體應(yīng)變儀進(jìn)行地殼形變動(dòng)態(tài)記錄時(shí)， 其具有較強(qiáng)的地震波記錄能力， 相當(dāng)于 “速度型” 無擺地震計(jì)[11]。 邱澤華等指出， 構(gòu)造地震發(fā)生時(shí)， 有可能使一些在 “應(yīng)力影” 區(qū)域內(nèi)的斷層釋放應(yīng)力后更趨于穩(wěn)定， 也可能因 “應(yīng)力觸發(fā)” 作用造成另一些斷層上的應(yīng)力增大， 進(jìn)而產(chǎn)生滑動(dòng)。 因此， 近震體波傳播會(huì)引起這些斷層上測(cè)點(diǎn)的形變儀器瞬間振動(dòng)[12]， 在曲線上表現(xiàn)為直接上升、 下降或突跳， 稱為階躍； 而遠(yuǎn)距離地震面波傳到測(cè)點(diǎn)時(shí)， 如果頻率與鉆孔體應(yīng)變儀的自振頻率相同則會(huì)引起共振， 就能記錄到脈沖或地震震蕩波形態(tài)的信號(hào)。

表1 小陶臺(tái)鉆孔體應(yīng)變儀觀測(cè)資料調(diào)和分析逐年統(tǒng)計(jì)表Table 1 Yearly statistics of harmonic analysis of borehole body strain data at Xiaotao Station

3.1 地震波記錄情況

小陶地震臺(tái)TJ-Ⅱ鉆孔體應(yīng)變儀自運(yùn)行以來，從統(tǒng)計(jì)數(shù)據(jù)及儀器地震波形圖上， 可看出映震能力較強(qiáng)， 多次記錄下強(qiáng)遠(yuǎn)震發(fā)生時(shí)的明顯同震響應(yīng)過程， 具有顯著振幅效應(yīng)以及豐富的應(yīng)變頻率，它能很好地反映出每個(gè)大地震發(fā)生時(shí)導(dǎo)致的地下應(yīng)變變化過程。

本文統(tǒng)計(jì) 2016 年 1 月 1 日至 2020 年 5 月 31日， 全球共發(fā)生7 級(jí)以上地震53 次， 主要分布在環(huán)太平洋板塊及板塊交界地區(qū)； 全國發(fā)生6 級(jí)以上地震22 次， 主要分布在青藏塊體、 新疆和臺(tái)灣地區(qū)； 華南地區(qū)未記錄到5 級(jí)以上地震； 臺(tái)灣地區(qū)發(fā)生5 級(jí)以上地震37 次。 震時(shí)小陶鉆孔體應(yīng)變儀共記錄到271 次全球強(qiáng)震同震響應(yīng)， 其中臺(tái)灣發(fā)生的37 次5.0 級(jí)以上地震全部記錄到； 對(duì)全球出現(xiàn)的7 級(jí)以上強(qiáng)震均能清楚地記錄下地震波整個(gè)過程， 且地震波波形振幅較大， 波形較完整；對(duì)全國出現(xiàn)的6 級(jí)以上地震波也具有很好記錄效果， 顯示出較強(qiáng)的地震響應(yīng)能力。

3.2 同震變化分析

地震分析預(yù)報(bào)時(shí)， 僅研究地震前的細(xì)微變化還不行， 首先應(yīng)熟悉地震發(fā)生的本質(zhì)， 尤其是地震同震響應(yīng)， 這是探索地震本質(zhì)的最關(guān)鍵數(shù)據(jù)[13]。小陶臺(tái)鉆孔體應(yīng)變儀記錄到的同震響應(yīng)形態(tài)以脈沖波、 地震震蕩波為主， 同時(shí)同震記錄的各種地震波也十分豐富、 時(shí)間持久、 形態(tài)完整。 當(dāng)最大振幅到達(dá)后， 后續(xù)應(yīng)變振幅隨時(shí)間推移而出現(xiàn)有規(guī)律衰減， 但各個(gè)地震的振幅存在不同， 部分呈現(xiàn)的是日變幅度的多倍， 也有部分被干擾覆蓋而未呈現(xiàn)。 有些地震同震波形表現(xiàn)為單純的壓性或張性， 還有的表現(xiàn)為對(duì)稱性， 這與地震發(fā)生地點(diǎn)有一定關(guān)系。 同震響應(yīng)振幅的大小不僅和地震震級(jí)大小有關(guān)， 還與震中距緊密聯(lián)系。 通常情況下，地震震級(jí)較大、 震中距較近， 那么同震振幅就比較大， 反之就越??； 但在震級(jí)相同情況下， 振幅與距離的關(guān)系還受其他因素影響。

2016 年 02 月 06 日臺(tái)灣高雄市（北緯 22.94°， 東經(jīng) 120.54°， 震中距 465 km）發(fā)生了 6.7 級(jí)地震，小陶鉆孔體應(yīng)變儀清晰地記錄到該震的地震波（圖4a）， 地震發(fā)生時(shí)刻 03 時(shí) 57 分， 地震波到達(dá)時(shí)刻04 時(shí)00 分， 高雄地震距離較近， 以快速變化為主， 變化形態(tài)為突降型， 下降幅度為 179.5×10-9；2017 年 9 月 8 日墨西哥 （北緯 15.05°， 西經(jīng) 93.90°，震中距14 388 km）發(fā)生了8.2 級(jí)地震， 同樣小陶鉆孔體應(yīng)變儀也清晰記錄到該震的地震波（圖4b），地震發(fā)生時(shí)刻12 時(shí)49 分， 地震波到達(dá)時(shí)刻13 時(shí)08 分， 墨西哥地震較遠(yuǎn)， 變化形態(tài)為震蕩型， 地震波最大振幅為 58.8×10-9， 振蕩持續(xù)時(shí)間 134 min。

邱澤華等曾指出， 觀測(cè)點(diǎn)應(yīng)力情況受震源處斷層錯(cuò)動(dòng)及測(cè)點(diǎn)構(gòu)造穩(wěn)定性的共同影響， 也就是如果測(cè)點(diǎn)構(gòu)造較為穩(wěn)定， 那么震源的應(yīng)力釋放不一定會(huì)使測(cè)點(diǎn)應(yīng)力狀態(tài)發(fā)生變化。 小陶臺(tái)鉆孔體應(yīng)變儀同震波僅有兩種形態(tài)： 近震、 小震為脈沖波， 遠(yuǎn)震、 強(qiáng)震為地震震蕩波， 沒有其他臺(tái)站所觀測(cè)到的階躍變化， 這可能因?yàn)樾√站@孔巖性致密， 無明顯裂隙發(fā)育， 所處的政和-海豐斷裂帶近年構(gòu)造活動(dòng)也較為平靜穩(wěn)定， 所以測(cè)點(diǎn)地下應(yīng)力場(chǎng)不易受影響而出現(xiàn)階躍變化。

圖4 小陶鉆孔體應(yīng)變儀記錄臺(tái)灣高雄6.7、墨西哥8.2 級(jí)地震的同震效應(yīng)曲線Fig.4 Coseismic effect curve recorded by Xiaotao body strain gauge during Gaoxiong 6.7 and Mosco 8.2 earthquakes

4 同震響應(yīng)分析

4.1 方法原理

我們所說的記震能力， 指觀測(cè)儀器能夠捕捉到的相同震中距條件下震級(jí)響應(yīng)下限， 或各檔震級(jí)中所能記錄到對(duì)應(yīng)的最遠(yuǎn)震中距。 從理論上來說， 相同震級(jí)， 震中距越小， 臺(tái)站觀測(cè)儀器能記錄到同震響應(yīng)概率越大； 而同等距離下， 觀測(cè)儀器記錄到同震變化的概率也隨震級(jí)的增大而增大。

我們通過分析鉆孔體應(yīng)變儀的觀測(cè)曲線， 整理出同震響應(yīng)目錄表， 并根據(jù)該目錄表與同期實(shí)際地震事件進(jìn)行對(duì)比， 使用震中距（D）和地震震級(jí)（MS）繪制出目標(biāo)臺(tái)站震級(jí)—震中距的對(duì)數(shù)關(guān)系圖，設(shè)定直線：

可將臺(tái)站觀測(cè)到的震例與未觀測(cè)到的地震事件， 在震級(jí)（MS）一震中距關(guān)系圖中分隔開， 則可根據(jù)式（1） 推導(dǎo)出不同震級(jí)的響應(yīng)震中距離。

4.2 資料選取原則

小陶臺(tái)鉆孔體應(yīng)變儀正式觀測(cè)于2016 年， 至今已積累大量典型震例。 本文統(tǒng)計(jì)的地震為2016—2020 年小陶臺(tái)鉆孔體應(yīng)變儀記錄到的地震，地震名稱、 發(fā)震時(shí)刻、 震級(jí)、 經(jīng)緯度均采用 “中國地震臺(tái)網(wǎng)目錄” 提供的參數(shù)。 震級(jí)統(tǒng)一采用面波震級(jí)MS， 沒有使用面波震級(jí)MS的， 則通過下式進(jìn)行相應(yīng)轉(zhuǎn)換：

首先， 統(tǒng)一將國家臺(tái)網(wǎng)速報(bào)目錄中的發(fā)震時(shí)間轉(zhuǎn)換為北京時(shí)間（+8）， 經(jīng)緯度也按正負(fù)表示（東經(jīng)、 北緯為正， 西經(jīng)、 南緯為負(fù)）， 震級(jí)使用 MS，統(tǒng)計(jì)國外 MS≥6 級(jí)； 國內(nèi) MS≥5 級(jí)； 福建省內(nèi)四級(jí)以上的地震， 總計(jì)582 次地震， 其中有同震反映的271 次， 占總數(shù)的46.56%。

其次， 利用地心余緯度公式， 計(jì)算各個(gè)地震事件到小陶臺(tái)（25.73°N， 117.13°E）的距離， 計(jì)算結(jié)果如表2：

式（4）中， Δ 為臺(tái)站震中距， λe為震中經(jīng)度， λs為臺(tái)站經(jīng)度， Φe為震中地心余緯度， Φs為臺(tái)站地心余緯度。 地心余緯度為

式（5）中， ψ 是地球中心緯度， 同地理緯度 φ 之間關(guān)系如下：

式（6）中， f 是地球的扁平率， 其 f=1/298.25≈0.003，如果忽略不計(jì)， 則 ψ=Φ， 即可將 ψ 近似看成 Φ"，于是

式（7）、（8）中， Φe為震中地理緯度， Φs為臺(tái)站地理緯度， 則

表2 小陶臺(tái)鉆孔體應(yīng)變儀地震記錄表Table 2 Earthquake catalog recorded by borehole body strain gauge at Xiaotao Station

若臺(tái)站震中距Δ 用長度D 可表示為

式（10）中， R 是地球平均半徑（6371 km）， Δ 是震中距圓心角（弧度）。

4.3 記震能力線

形變觀測(cè)儀器對(duì)地震記震能力大小， 一般通過震級(jí)和震中距相關(guān)性來表示。 以同震變化的震中距對(duì)數(shù)lgD（以10 為底）為橫坐標(biāo)， 以地震MS震級(jí)為縱坐標(biāo)， 繪制小陶臺(tái)鉆孔體應(yīng)變儀記震能力圖。由圖5 可見， 有同震響應(yīng)記錄（黑色實(shí)心圓點(diǎn)）與沒有同震響應(yīng)記錄（黑色實(shí)心三角）的地震分布顯示在圖中， 在分布圖中可直觀以一條直線分割圓點(diǎn)集中區(qū)與三角集中區(qū)。 直線上部分布基本為圓點(diǎn)， 即有同震響應(yīng)記錄點(diǎn)， 直線下部則基本為三角， 即無同震響應(yīng)記錄點(diǎn)， 其中震級(jí)5～6 級(jí)的響應(yīng)地震分布震中距對(duì)數(shù) 2.70～3.11 （501～1288 km）范圍內(nèi)， 震級(jí)6～7 級(jí)的響應(yīng)地震分布震中距對(duì)數(shù)3.09～3.63（1230～4265 km）范圍內(nèi)。

該直線即為小陶臺(tái)鉆孔體應(yīng)變儀記震能力線，在直線上任取兩點(diǎn)， 代入直線方程：

式（11）中， lgD1是首點(diǎn)橫向坐標(biāo)（即震中距對(duì)數(shù)），M1是首點(diǎn)縱向坐標(biāo) （即震級(jí)）， lgD2是次點(diǎn)橫向坐標(biāo)（即震中距對(duì)數(shù)）， M2是次點(diǎn)縱向坐標(biāo)（ 震級(jí)），可推導(dǎo)出小陶臺(tái)鉆孔體應(yīng)變儀映震能力數(shù)學(xué)表達(dá)式為：

式（12）中 MS為震級(jí)， D 為震中距離。

由式（12）計(jì)算出小陶臺(tái)鉆孔體應(yīng)變儀地震響應(yīng)范圍為： 震級(jí) 5.5 級(jí)、 震中距 480 km； 震級(jí) 6級(jí)、 震中距 984 km； 震級(jí) 7.0 級(jí)、 震中距 4130 km； 震級(jí) 7.5 級(jí)、 震中距 8472 km； 由于地球半徑為6371 km， 所以7.5 級(jí)以上地震， 小陶臺(tái)鉆孔體應(yīng)變儀基本就有映震反映， 該結(jié)果與圖5 顯示相符。

4.4 最大振幅和震級(jí)相關(guān)性

如果地震出現(xiàn)在相同位置， 我們可把它們的震中距當(dāng)成一樣， 分析最大振幅和震級(jí)二者之間的相關(guān)性。 本文選取2016 至2020 年日本發(fā)生的近海地震， 挑選了發(fā)震時(shí)刻、 同震波形清晰明顯的地震， 除去外界影響較嚴(yán)重部分的觀測(cè)資料，從而實(shí)現(xiàn)最大振幅的選取。 我們選擇原則是震中距相似而震級(jí)各異， 將MS震級(jí)設(shè)成橫向坐標(biāo)， 將最大振幅設(shè)成縱向坐標(biāo)， 擬合出小陶臺(tái)鉆孔體應(yīng)變儀最大振幅與震級(jí)關(guān)系圖（圖6）。

圖5 小陶臺(tái)鉆孔體應(yīng)變儀記震能力線示意圖Fig.5 Schematic diagram of seismic recording capacity line of borehole body strain gauge at Xiaotao Station

圖6 小陶臺(tái)鉆孔體應(yīng)變儀震級(jí)與最大振幅關(guān)系圖Fig.6 The relationship between the magnitude and the maximum amplitude of borehole body strain gauge at Xiaotao Station

本文最大振幅與震級(jí)相關(guān)性分析采用的最大地震是6.7 級(jí)， 由圖6 能發(fā)現(xiàn)， 統(tǒng)計(jì)資料大部分位于同一曲線附近。 如果地震超過6 級(jí)， 那么最高振幅將伴隨地震震級(jí)增大而變大， 與指數(shù)關(guān)系相契合。 日本2016 年4 月16 日的7.3 級(jí)地震最大振幅為143×l0-9， 因震級(jí)和振幅太大， 本次統(tǒng)計(jì)予以舍棄。

鉆孔體應(yīng)變儀記錄的地震震級(jí)同最大振幅基本呈冪指數(shù)相關(guān)性， 伴隨地震震級(jí)增大， 最大振幅增長速率也提升。 從圖6 可看出， 即使人工剔除偏差較大的地震， 可部分點(diǎn)依然非常離散， 從而得知震中距與最高振幅之間， 并非全部是反例關(guān)系。 通過相關(guān)分析可知， 形變觀測(cè)儀器記錄的最高振幅， 并非是地震波表現(xiàn)的最高振幅值。 由于鉆孔體應(yīng)變儀采樣率不高， 為分鐘值采樣， 導(dǎo)致最高振幅會(huì)出現(xiàn)遺漏， 還有可能因?yàn)橹芷诒鎰e不準(zhǔn)或者錯(cuò)誤， 導(dǎo)致地震最高振幅的記錄出現(xiàn)較大偏離。

5 結(jié)論

（1）小陶臺(tái)鉆孔應(yīng)變儀觀測(cè)資料具有較高連續(xù)性與穩(wěn)定性。 在日常觀測(cè)中， 可記錄到清晰、 完整的固體潮曲線， 并且潮汐因子 “內(nèi)精度” 指標(biāo)遠(yuǎn)超過形變臺(tái)網(wǎng)規(guī)范要求， 可為預(yù)報(bào)分析、 研究所用。

（2）小陶臺(tái)鉆孔體應(yīng)變儀記錄到的同震響應(yīng)形態(tài)主要有兩類： 脈沖波、 地震波。 距離較遠(yuǎn)的大震面波發(fā)育， P 波也會(huì)存在部分記錄， 波形特點(diǎn)大部分顯示為： 曲線上下震動(dòng)變化， 后期振幅慢慢變小并逐漸恢復(fù)至正常日變。 據(jù)有關(guān)科技文獻(xiàn)記載， 鉆孔體應(yīng)變儀可采集到同震響應(yīng)與震后曲線發(fā)生的階躍變化， 但小陶臺(tái)至今尚未采集到曲線出現(xiàn)階躍的數(shù)據(jù)， 說明該測(cè)點(diǎn)構(gòu)造較為穩(wěn)定， 不易受到震源應(yīng)力釋放的影響。

（3）小陶臺(tái)鉆孔體應(yīng)變儀觀測(cè)時(shí)間長， 積累震例多， 大震響應(yīng)率為46.56%， 地震響應(yīng)范圍為： 震級(jí) 5.5 級(jí)、 震中距 480 km； 震級(jí) 6 級(jí)、 震中距 984 km； 震級(jí) 7.0 級(jí)、 震中距 4130 km； 震級(jí) 7.5 級(jí)、震中距8472 km； 由于地球半徑為6371 km， 所以7.5 級(jí)以上地震基本都有映震反映。

（4）鉆孔體應(yīng)變儀對(duì)同震應(yīng)變信息的捕捉能力與響應(yīng)形態(tài)較好， 干擾信息對(duì)其影響較小， 僅與地震震級(jí)、 震中距直接相關(guān)， 記震能力線與其他方法相比， 不僅能直觀反映測(cè)點(diǎn)構(gòu)造情況， 也可量化顯示鉆孔體應(yīng)變儀的映震能力。