徐芳芳 朱成林 賈 媛 李 杰 周長(zhǎng)禮 周均太

1 山東省地震局，濟(jì)南市文化東路20號(hào)，250014

鉆孔應(yīng)變觀測(cè)是研究地殼運(yùn)動(dòng)的手段之一，主要分為體積式和分量式(包括RZB-2和YRY-4分量式鉆孔應(yīng)變)兩種，觀測(cè)分辨率達(dá)10-10，可以記錄到各種應(yīng)變潮汐信息。鉆孔應(yīng)變觀測(cè)是通過(guò)測(cè)量鉆孔孔徑變化來(lái)直接或間接計(jì)算地殼應(yīng)力-應(yīng)變狀態(tài)相對(duì)變化的，可記錄到與地震的孕育及發(fā)生直接相關(guān)的應(yīng)力-應(yīng)變信息。目前，國(guó)內(nèi)許多專(zhuān)家已在相關(guān)領(lǐng)域開(kāi)展了研究工作[1-2]，部分學(xué)者對(duì)客觀評(píng)定鉆孔應(yīng)變四分量觀測(cè)資料質(zhì)量的方法進(jìn)行了分析[3-5]，為科學(xué)應(yīng)用四分量鉆孔應(yīng)變觀測(cè)資料奠定了基礎(chǔ)。本文選取鄒城地震臺(tái)RZB-2分量式鉆孔應(yīng)變2020-01～2021-07更換探頭前后的數(shù)據(jù)資料，利用相對(duì)標(biāo)定、四分量自檢分析、M2波潮汐因子相對(duì)中誤差計(jì)算及與同臺(tái)站TJ-2體應(yīng)變觀測(cè)資料進(jìn)行相關(guān)性對(duì)比等方法，計(jì)算更換探頭前后儀器的內(nèi)精度、自檢情況及M2波潮汐因子精度，以綜合分析分量式鉆孔應(yīng)變觀測(cè)系統(tǒng)的穩(wěn)定性和觀測(cè)數(shù)據(jù)的可靠性。

1 臺(tái)站及儀器概況

鄒城地震臺(tái)位于鄒城市吉玉山東南麓，海拔138 m，地處尼山斷裂與嶧山斷裂之間，地勢(shì)北高南低，地表出露巖性為花崗巖，地下水埋深4～5 m。臺(tái)站周?chē)鸁o(wú)大型工業(yè)廠房施工，周邊有幾處作坊使用地下水，對(duì)臺(tái)站儀器產(chǎn)生一定影響。RZB-2分量式鉆孔應(yīng)變安裝于2017-11-29，鉆孔開(kāi)口孔徑為146 mm，終孔孔徑為130 mm，孔井基巖完整，探頭深度81 m，應(yīng)變探頭內(nèi)應(yīng)變傳感器由上到下依次呈順時(shí)針45°排列，NS分量方位角為NW5°，EW分量方位角為SW5°，NE分量方位角為NE40°，NW分量方位角為NW40°，于2019-06接入前兆數(shù)據(jù)庫(kù)。因EW分量故障，2020-10-29更換探頭并重新安裝，新探頭深度73 m，方位角不變，采樣率均為秒采樣。自運(yùn)行以來(lái)，儀器觀測(cè)數(shù)據(jù)固體潮清晰，有年變趨勢(shì)。TJ-2體應(yīng)變安裝于2015-12-12，鉆孔深度71 m，探頭在井下68 m處，水位探頭在水下8 m處，2017-02正式加入中國(guó)地震局前兆臺(tái)網(wǎng)管理系統(tǒng)，觀測(cè)數(shù)據(jù)質(zhì)量較好，2臺(tái)儀器相距5 m。

2 評(píng)價(jià)指標(biāo)

2.1 相對(duì)標(biāo)定自檢內(nèi)精度

四分量應(yīng)變自洽好壞關(guān)系到觀測(cè)數(shù)據(jù)的可信度。儀器探頭安裝最重要的條件是孔井基巖要完整，當(dāng)孔井基巖完整時(shí)，儀器觀測(cè)系統(tǒng)基本符合自洽方程。而在有完整基巖的孔井中，也可能因?yàn)樵`敏度發(fā)生變化，造成與實(shí)驗(yàn)室靈敏度不同，導(dǎo)致觀測(cè)偏離自洽。邱澤華等[3]根據(jù)四分量鉆孔應(yīng)變觀測(cè)的特點(diǎn)，提出四分量鉆孔應(yīng)變觀測(cè)的實(shí)地相對(duì)標(biāo)定方法，唐磊等[4-5]則在此基礎(chǔ)上改進(jìn)了相對(duì)標(biāo)定自檢內(nèi)精度法，進(jìn)一步保證客觀地評(píng)價(jià)數(shù)據(jù)質(zhì)量。相對(duì)標(biāo)定自檢內(nèi)精度法的具體計(jì)算原理見(jiàn)文獻(xiàn)[4]。

根據(jù)月數(shù)據(jù)，選用整點(diǎn)值計(jì)算得到2020-01～2021-07鄒城臺(tái)RZB-2分量式鉆孔應(yīng)變觀測(cè)資料的相對(duì)標(biāo)定結(jié)果及年平均相對(duì)標(biāo)定結(jié)果，具體見(jiàn)表1。通過(guò)對(duì)比探頭更換前(2020-01～09)和更換后(2020-12～2021-07，因2020-11儀器零漂較大不參與計(jì)算)的結(jié)果發(fā)現(xiàn)，更換探頭后鉆孔應(yīng)變各分量每月K值明顯增大，偏差S明顯減小，自檢內(nèi)精度指標(biāo)α數(shù)值大幅度降低，表明更換探頭后觀測(cè)數(shù)據(jù)的內(nèi)精度明顯增高。

表1 鄒城臺(tái)RZB-2鉆孔應(yīng)變2020-01～2021-07觀測(cè)資料相對(duì)標(biāo)定結(jié)果

鉆孔施工方、儀器安裝專(zhuān)家及臺(tái)站技術(shù)人員對(duì)更換探頭前觀測(cè)數(shù)據(jù)內(nèi)精度較低的情況共同進(jìn)行了分析。首先，由于鉆井工藝問(wèn)題，在探頭安裝并灌注水泥后，泥漿無(wú)法完全固結(jié)，導(dǎo)致探頭與圍巖耦合效果不好[2,6]。而更換探頭時(shí)，泥漿已經(jīng)過(guò)長(zhǎng)時(shí)間沉淀，探頭與圍巖耦合效果較好，因此內(nèi)精度較高。由于臺(tái)站周邊存在抽水干擾，需結(jié)合四分量自檢和M2潮汐因子相對(duì)中誤差對(duì)數(shù)據(jù)質(zhì)量進(jìn)行進(jìn)一步分析。

2.2 數(shù)據(jù)自檢分析

鉆孔應(yīng)變各觀測(cè)時(shí)刻滿足雙襯套模型觀測(cè)方程組[7]：

(1)

式中，Si(i=1,2,3,4)為分量式鉆孔應(yīng)變的4個(gè)觀測(cè)分量，ε1為最大主應(yīng)變，ε2為最小主應(yīng)變，φ為最大主應(yīng)變方向，θ為第1個(gè)分量的方位角，A、B為與鉆孔各層材料的彈性參數(shù)和半徑相關(guān)的耦合系數(shù)。據(jù)此方程組可以得到：

S1+S3=S2+S4=2A(ε1+ε2)

(2)

式(2)即為四分量鉆孔應(yīng)變觀測(cè)特有的自恰方程，可以計(jì)算面應(yīng)變，從而驗(yàn)證鉆孔應(yīng)變觀測(cè)結(jié)果的可靠性和儀器數(shù)據(jù)的穩(wěn)定性[8]。采用整點(diǎn)值計(jì)算2組面應(yīng)變的相關(guān)系數(shù)，相關(guān)系數(shù)越接近1表明觀測(cè)系統(tǒng)越穩(wěn)定。為對(duì)比降雨的影響，采用鉆孔水位數(shù)據(jù)進(jìn)行對(duì)比分析(圖1)，計(jì)算結(jié)果表明，RZB-2鉆孔應(yīng)變?cè)诟鼡Q探頭前2組面應(yīng)變相關(guān)性較低(圖2(a))，更換探頭后2組面應(yīng)變相關(guān)系數(shù)較高，接近1，具有很好的同步特征(圖2(b))，說(shuō)明鄒城臺(tái)RZB-2分量式鉆孔應(yīng)變?cè)诟鼡Q探頭后觀測(cè)系統(tǒng)穩(wěn)定，數(shù)據(jù)符合四分量鉆孔應(yīng)變觀測(cè)原理。

圖1 鄒城臺(tái)RZB-2鉆孔應(yīng)變數(shù)據(jù)自檢分析Fig.1 Surface strain sliding correlation analysis of RZB-2 borehole strain at Zoucheng station

圖2 鄒城臺(tái)RZB-2鉆孔應(yīng)變2組面應(yīng)變滑動(dòng)相關(guān)性分析Fig.2 Observation system self-checking analysis of RZB-2 borehole strain at Zoucheng station

2.3 潮汐因子相對(duì)中誤差分析

應(yīng)變固體潮M2波潮汐因子的內(nèi)精度是指M2波潮汐因子與M2波潮汐因子中誤差的比值，即潮汐因子相對(duì)中誤差[6]?？捎迷撝笜?biāo)數(shù)值大小來(lái)說(shuō)明觀測(cè)固體潮數(shù)據(jù)的可靠性，相對(duì)中誤差(內(nèi)精度)數(shù)值越小，說(shuō)明觀測(cè)值越可靠。選取鄒城臺(tái)

RZB-2分量式鉆孔應(yīng)變和TJ-2體應(yīng)變2020-01～2021-07整點(diǎn)值觀測(cè)資料，利用維尼迪柯夫調(diào)和分析方法進(jìn)行逐月分析，計(jì)算2套鉆孔應(yīng)變M2波潮汐因子相對(duì)中誤差。由圖3可以看出，鄒城臺(tái)RZB-2鉆孔應(yīng)變潮汐因子內(nèi)精度在2020-01～08較差，平均為0.016 7，TJ-2體應(yīng)變內(nèi)精度較高，平均為0.008 2；2021-01以來(lái)，鉆孔應(yīng)變各分量數(shù)據(jù)內(nèi)精度較差，平均為0.073 7，體應(yīng)變平均內(nèi)精度為0.008 4。經(jīng)現(xiàn)場(chǎng)核實(shí)發(fā)現(xiàn)，2020-12在距離臺(tái)站東側(cè)約200 m處新增一作坊，不定時(shí)通過(guò)30 m深抽水井進(jìn)行抽水，導(dǎo)致干擾明顯，在2021-01及春季3～5月尤其嚴(yán)重，EW分量曲線干擾顯著。根據(jù)TJ-2體應(yīng)變和RZB-2鉆孔應(yīng)變面應(yīng)變相關(guān)性分析(圖4)可知，雖然鉆孔應(yīng)變潮汐因子的內(nèi)精度較差，但2套儀器的面應(yīng)變相關(guān)性較高，說(shuō)明觀測(cè)系統(tǒng)的運(yùn)行是穩(wěn)定可靠的；而TJ-2體應(yīng)變幾乎不受抽水干擾，RZB-2鉆孔應(yīng)變受抽水干擾嚴(yán)重，可能由于TJ-2體應(yīng)變與RZB-2鉆孔應(yīng)變觀測(cè)的分辨率相差2個(gè)量級(jí)，且RZB-2鉆孔應(yīng)變?cè)隗w應(yīng)變東側(cè)5 m處，距離抽水井更近。另外，鉆井之間可能存在一定的裂隙，導(dǎo)致儀器對(duì)干擾的響應(yīng)程度不一致。

圖3 鄒城臺(tái)RZB-2鉆孔應(yīng)變更換探頭前后觀測(cè)精度曲線Fig.3 Observation precision curves before and after change sensor probe of RZB-2 borehole strain at Zoucheng station

2.4 相關(guān)性分析

鄒城臺(tái)TJ-2體應(yīng)變2017年正式加入中國(guó)地震前兆臺(tái)網(wǎng)以來(lái)，觀測(cè)資料精度較高。RZB-2分量式鉆孔應(yīng)變和TJ-2體應(yīng)變的探頭安裝在同一深度，相距僅5 m。利用窗長(zhǎng)為20 d的整點(diǎn)值進(jìn)行日相關(guān)性滑動(dòng)分析，為方便比較，將2套儀器數(shù)據(jù)均換算至10-8量級(jí)。更換探頭前，RZB-2鉆孔應(yīng)變的2組面應(yīng)變和TJ-2體應(yīng)變的面應(yīng)變?nèi)障嚓P(guān)系數(shù)相對(duì)較低且不穩(wěn)定，多集中在0.8～0.95范圍內(nèi)(圖4(a))；更換探頭后，2套儀器的面應(yīng)變?nèi)障嚓P(guān)系數(shù)接近1(圖4(b))。將2套儀器的整點(diǎn)值時(shí)間序列進(jìn)行歸一化處理，圖5中紅色虛線為同步變化產(chǎn)生的線性相關(guān)性，可以看出，更換探頭前儀器的面應(yīng)變一致性較差，更換探頭后面應(yīng)變的一致性較好。綜合分析發(fā)現(xiàn)，更換探頭前儀器觀測(cè)結(jié)果的一致性、內(nèi)精度及自檢情況均較差，說(shuō)明RZB-2分量式鉆孔應(yīng)變的觀測(cè)系統(tǒng)不穩(wěn)定，觀測(cè)數(shù)據(jù)可靠性較差；更換探頭后儀器的觀測(cè)系統(tǒng)正常，2套儀器的觀測(cè)結(jié)果一致性好，數(shù)據(jù)可靠性高。

圖4 鄒城臺(tái)RZB-2鉆孔應(yīng)變更換探頭前后與TJ-2體應(yīng)變面應(yīng)變滑動(dòng)相關(guān)性分析Fig.4 The sliding correlation analysis of surface strain before and after change sensor probe of RZB-2 borehole strain with TJ-2 volumetric strain at Zoucheng station

圖5 鄒城臺(tái)RZB-2鉆孔應(yīng)變更換探頭前后與TJ-2體應(yīng)變散點(diǎn)分布Fig.5 Scatter distribution of correlated time series of surface strain before and after change sensor probe of RZB-2 borehole strain with TJ-2 volumetric strain at Zoucheng station

3 結(jié) 語(yǔ)

鄒城臺(tái)RZB-2分量式鉆孔應(yīng)變第1次安裝后內(nèi)精度較差，2組面應(yīng)變不符合四分量自檢方程，面應(yīng)變相關(guān)性較差，主要原因是探頭與圍巖耦合不好；第2次更換探頭后，RZB-2鉆孔應(yīng)變的內(nèi)精度較高，2組面應(yīng)變符合四分量自檢方程，面應(yīng)變相關(guān)系數(shù)接近1，自檢曲線非常好，與體應(yīng)變面應(yīng)變的相關(guān)性也接近1，說(shuō)明RZB-2鉆孔應(yīng)變的觀測(cè)系統(tǒng)正常穩(wěn)定，數(shù)據(jù)可靠。2021年臺(tái)站周邊新增一工廠導(dǎo)致抽水干擾增多，數(shù)據(jù)潮汐因子內(nèi)精度較差，平均值為0.07左右，嚴(yán)重影響觀測(cè)數(shù)據(jù)質(zhì)量，建議盡快采取相關(guān)環(huán)境保護(hù)措施，減小干擾。根據(jù)分析結(jié)果可知，日常工作中僅憑潮汐因子內(nèi)精度或相對(duì)標(biāo)定自檢內(nèi)精度來(lái)評(píng)價(jià)儀器的工作情況及數(shù)據(jù)質(zhì)量是不全面的，應(yīng)綜合考慮各種參數(shù)。