楊雅慧 劉洪良 張展偉 張慶揚 昝書林

1 河北紅山巨厚沉積與地震災(zāi)害國家野外科學(xué)觀測研究站，河北省邢臺市隆堯縣山口鎮(zhèn)，055350 2 河北省地震局保定地震監(jiān)測中心站，河北省保定市天威中路206號，071000

河北測區(qū)流動重力測網(wǎng)自1968年開始觀測，最早由邢臺地區(qū)4個閉合環(huán)組成，之后陸續(xù)建立華北火車測網(wǎng)、文霸汽車測網(wǎng)、邯安測網(wǎng)，2009年火車測網(wǎng)停測，文霸、邯安、邢臺測網(wǎng)組合成一個橫跨河北及周邊省市的大型重力測網(wǎng)并與全國測網(wǎng)連接，構(gòu)成大華北流動重力測網(wǎng)的一部分，并于2010年開始正式測量[1-2]。近年來，通過本測區(qū)觀測結(jié)果發(fā)現(xiàn)，河北東南部太行山以東區(qū)域重力場正異常累積明顯。為進一步了解本測區(qū)重力場長期動態(tài)變化機制及影響因素，為地震前兆信息的辨別及中長期地震預(yù)報提供理論依據(jù)，本文基于河北測區(qū)2015～2020年積累的12期相對重力復(fù)測資料，并搜集周邊相鄰測區(qū)的相對重力復(fù)測數(shù)據(jù)，同時以區(qū)域內(nèi)白家疃、太原、泰安、薊縣等絕對點同期重力點值作為控制，得到高精度的重力平差結(jié)果，獲得重力場動態(tài)變化圖像；然后采用多項式擬合方法計算重力場的平均年變率，結(jié)合已有數(shù)據(jù)及研究成果，對本測區(qū)重力場變化機制進行分析。

1 觀測數(shù)據(jù)與處理

河北流動重力測網(wǎng)在2015～2020年每年3月、8月各進行一期觀測工作，每期觀測用時約為45 d。使用2臺CG-5型相對重力儀進行觀測，測量精度優(yōu)于5 μGal，讀數(shù)分辨率為1 μGal[3]。

整理河北及周邊測區(qū)重力資料及絕對點信息，使用LGADJ處理軟件進行絕對重力點控制下的平差處理，得到河北及部分周邊地區(qū)的平差結(jié)果。參與計算的各期測點總數(shù)為285～550個，測段數(shù)為1 274～2 642條，平差后平均點值精度優(yōu)于16 μGal，聯(lián)測精度優(yōu)于17 μGal(表1)。

表1 重力數(shù)據(jù)處理精度Tab.1 Processing accuracy of the gravity data

參與計算的河北及周邊地區(qū)重力資料的整體測網(wǎng)分布如圖1所示。其他測區(qū)使用儀器多為CG-5型相對重力儀及LCR-G型相對重力儀。

圖1 2015年河北及周邊地區(qū)重力測網(wǎng)分布Fig.1 Gravity network distribution of Hebei and surrounding areas in 2015

2 重力場動態(tài)特征

重力場的動態(tài)變化可以反映重力場在研究時段內(nèi)的持續(xù)演變情況，由上述重力平差結(jié)果可得到重力場累積動態(tài)變化圖(圖2)。可以看出，2015年重力場累積變化較小，大部分區(qū)域呈現(xiàn)低負(fù)值，僅在石家莊西部、邯鄲北部、天津東北部出現(xiàn)小范圍的低正值，整體重力在-20～10 μGal之間變化。在2015～2016年的重力累積圖像中，零值線呈NE-SW向分布，其中，薊縣北部出現(xiàn)輕微負(fù)異常累積，最高達(dá)-40 μGal，并在之后的2017年間持續(xù)累積，至2018年底負(fù)異常消失，之后變化平穩(wěn)，未出現(xiàn)異常高值。另一方面，位于河北東南部、石家莊以東區(qū)域的天津-邯鄲一線自2015年起出現(xiàn)正異常并持續(xù)累積，至2017年高正值異常明顯，并逐漸向北東方向延伸，2020年底形成走向近NE-SW的橢圓形高正值異常區(qū)。該高值區(qū)西部以太行山為界，北部至燕山一帶，東部與山東地區(qū)重力累積高值區(qū)相連接；零值線沿太行山山前斷裂帶展布，高梯度帶明顯，累積重力場正負(fù)差異最高達(dá)120 μGal，區(qū)域內(nèi)重力場異常變化顯著。2015～2020年，在異常區(qū)高梯度附近相繼發(fā)生臨城MS3.7、永清MS4.3、豐南MS4.5、古冶MS5.1地震，重力異常仍在持續(xù)。其他地區(qū)重力場累積變化平緩，無明顯異常。

圖2 重力場累積動態(tài)變化Fig.2 Cumulative dynamic changes of gravity field

3 平均年變率變化

平均年變率可以剔除重力場的短期、局部影響因素，從而反映出由于地下活動或地下物質(zhì)密度變化引起的重力場長時間變化趨勢。為保證數(shù)據(jù)結(jié)果的可靠性，在進行重力場年變率相關(guān)計算之前，首先對參與計算的觀測數(shù)據(jù)進行篩選，刪除其中測量期次少于6期的測點數(shù)據(jù)，避免部分測點由于數(shù)據(jù)量少造成誤差。

假設(shè)重力觀測數(shù)據(jù)g(φ,λ,Δt)包含長期變化、周期性變化以及噪聲信息，則g(φ,λ,Δt)可寫為[4-6]：

g(φ,λ,Δt)=A+BΔt+

(1)

式中，B為重力場變化速率，Δt為時間差，i=1、2分別表示年周期、半年周期項，Ci和Di為振幅。由于本文使用流動觀測重力數(shù)據(jù)，觀測周期長，周期變化不明顯，故在實際計算中僅保留長期趨勢項進行擬合。

最終共有478個重力測點參與計算，其中河北地區(qū)有181個。部分高值區(qū)重力測點擬合情況如圖3所示，2015～2020年重力平均年變率情況如圖4所示。

圖3 重力測點擬合結(jié)果Fig.3 Gravity measuring points fitting results

圖4 重力平均年變率Fig.4 Average annual gravity change rate

圖3(a)～3(d)依次按測點地理位置由北至南排序，平均年變率范圍為4.57～26.90 μGal/a。圖3(b)中測點擬合相關(guān)系數(shù)最小，其余3個測點R2均在0.5以上，擬合效果較好。重力值均呈現(xiàn)整體上升的累積態(tài)勢。

從圖4可以看出，河北北部、西部區(qū)域重力平均年變率較小，多在-10～10 μGal/a區(qū)間內(nèi)變化，R2<0.5居多，重力場多為季節(jié)性等短期變化，長期趨勢變化無明顯異常；而東南部天津-邯鄲一線平均年變率高值顯著，最大值可達(dá)27 μGal/a，重力持續(xù)正向累積。高值區(qū)域整體走向近NE-SW向，與累積重力變化圖中高值區(qū)域形態(tài)及位置大致相符，但整體范圍相對縮小。區(qū)域內(nèi)年變率擬合相關(guān)系數(shù)R2≥0.7居多，擬合結(jié)果較好，這些區(qū)域一直處于較規(guī)律的變動狀態(tài)。

為排除測量期間測點發(fā)生變動對重力平差結(jié)果的影響，對河北東南部重力平均年變率高值區(qū)域內(nèi)重力測點逐一進行點位變動情況核實。經(jīng)統(tǒng)計，共核實重力測點103個，在2015～2020年測量期間點位發(fā)生變動的測點為6個。基于其發(fā)生變動前的部分觀測數(shù)據(jù)進行擬合，然后與全部數(shù)據(jù)擬合結(jié)果對比，發(fā)現(xiàn)兩者重力平均年變率差值在0.2～5.8 μGal/a之間。

根據(jù)已有研究，河北地區(qū)重力變化主要受地表沉降及地下水流失的影響[7-8]。地表沉降會使區(qū)域重力值上升，地下水流失則會使區(qū)域重力值下降。2015～2020年河北東南部的天津-邯鄲一線重力累積變化及平均年變率均呈現(xiàn)明顯高值，表明地表沉降的影響覆蓋了地下水流失的影響，過度使用地下水的情況可能有所改善，與以往研究結(jié)果相符[5-7]。

為進一步研究地下水對重力場的影響，由北至南選取北京南部的河間觀測井(井1)、石家莊觀測井(井2)以及石家莊南部永年觀測井(井3)的水位數(shù)據(jù)(具體井位置見圖1)與觀測井最近的重力測點(測點1、2、3)數(shù)據(jù)進行對比，結(jié)果如圖5所示。

從圖5(a)中看出，在2015～2020年3口觀測井的水位均呈不同程度的下降趨勢，造成重力值的下降；而圖5(b)中重力點值呈現(xiàn)上下波動，與水位的持續(xù)下降無明顯對應(yīng)。故短期來看，地下水流失對重力的影響相對被覆蓋。在長期變化趨勢方面，除測點2呈現(xiàn)輕微下降趨勢外，測點1、3均為整體上升趨勢。由于河北地區(qū)地勢多平坦，對于觀測井附近一定區(qū)域，可以按照無限平面近似公式計算其水位變化對重力場的影響[9]：

圖5 水位和重力測點變化Fig.5 Variation of groundwater and gravity measuring points

Δg=2πGφ(H2-H1)=2πGφΔh

(2)

式中，G=6.672 59×10-11Nm2/kg2，φ為給水密度，Δh為水位變化。根據(jù)已有研究，觀測井1的φ取0.25 g/cm3[9]，觀測井2和3的φ均取0.15 g/cm3[10]。

若不考慮其他影響因素，將重力平均年變率全部歸于地下水流失及地表沉降的影響，則3個觀測井的計算結(jié)果如表2所示。

表2 水位變化及地表沉降影響Tab.2 The effects of groundwater variation and surface subsidence

由表2可見，除井1水位變化對重力平均年變率的影響稍大外，其他2口井對重力平均年變率的影響均相對較??；地表沉降對重力平均年變率的影響明顯高于地下水流失。根據(jù)已有地表沉降速率研究結(jié)果，天津-邯鄲一帶重力場變化高值區(qū)多位于沉降中-強發(fā)育地區(qū)，沉降速率在30 mm/a以上，其中強發(fā)育地區(qū)達(dá)到50 mm/a以上[8]，僅地表沉降則至少可以引起重力場10 μGal/a的變化，強發(fā)育地區(qū)至少可以引起重力場17 μGal/a的變化。而重力平均年變率計算結(jié)果中，該部分區(qū)域重力變化多在15～30 μGal/a，兩者具有較好的一致性。

4 結(jié) 語

1)為提高數(shù)據(jù)處理的精度，除本測區(qū)觀測資料外，另外添加周邊測區(qū)相對重力觀測資料及絕對重力值進行整體平差控制，使最終點值精度控制在16 μGal以內(nèi)，數(shù)據(jù)結(jié)果更加精準(zhǔn)；在計算年變率時，首先對測點進行篩選，使結(jié)果能更好地反映重力場的真實變化。

2)基于長時間的觀測結(jié)果，可以分析得出重力場長期的背景趨勢變化情況。在重力累積動態(tài)變化圖中，河北東南部天津-邯鄲一線重力場正異常持續(xù)累積，高值區(qū)逐漸擴大，最終在2020年底形成走向近NE-SW的橢圓形高值異常區(qū)；北至燕山附近，西至太行山，正負(fù)差異最高可達(dá)120 μGal。在重力平均年變率結(jié)果中，正高值區(qū)域分布與重力累積動態(tài)變化結(jié)果基本一致，且該區(qū)域擬合相關(guān)系數(shù)較高。說明除去季節(jié)性變化等短期、局部影響因素外，重力變化呈現(xiàn)長期規(guī)律上升態(tài)勢。

3)綜合上述結(jié)論及河北地區(qū)高值區(qū)域內(nèi)選取的部分水位數(shù)據(jù)、相關(guān)地表沉降研究結(jié)果得出，相對于地表沉降因素，重力場受地下水流失的影響較小，綜合表現(xiàn)為地表沉降引起的正異常累積，沉降發(fā)育嚴(yán)重區(qū)域影響至少可達(dá)17 μGal/a。