張娜　龔燕民　王艷

摘要河北省易縣地震臺 DSQ 型水管傾斜儀 NE 向2022年2月14日以來觀測數(shù)據(jù)曲線出現(xiàn)趨勢下降 E 傾變化，年變觀測數(shù)據(jù)曲線不同于往年。對該異常進行了分析，并結(jié)合三維集中荷載模型模擬計算結(jié)果表明：2022年2月14日—4月6 日旺隆水庫共泄流總庫容量達到279.5萬 m3，造成傾斜變化理論值為1.27×10?7，高于傾斜干擾的限定指標1.455×10?8一個量級，推測旺隆水庫泄流荷載是造成易縣臺水管儀 NE 分量加速異常變化的主要因素。

關(guān)鍵詞水管傾斜儀；荷載；水庫；異常

中圖分類號： P315.63文獻標識碼： A文章編號：2096-7780（2023）05-0228-06

doi：10.19987/j.dzkxjz.2022-146

Analysis of water tube tiltmeter observation affected by reservoir load

Zhang Na1），Gong Yanmin2），Wang Yan1）

1） Baoding Seismic Monitoring Center Station， Hebei Baoding 071051， China 2） Yixian Seismic Station， Hebei Yixian 071000， China

AbstractThe NE component curve of DSQ water tube tiltmeter at Yixian seismic station in Hebei Province had a downward trend since February 14，2022. The annual change curve was different from that of previous years. Combined with the simulation results of three-dimensional concentrated load model，it is shown that the total reservoir capacity of Wanglong reservoir reached 279.50×104 m3 from February 14 to April 6，2022. The theoretical value of tilt change was 1.27×10?7，which is higher than the limit index of tilt interference of 1.455×10?8. It was inferred as the main factor causing the abnormal change of NE component acceleration of water tube instrument at Yixian station.

Keywordswater tube tiltmeter; load; reservoir; abnormal

引言

河北省易縣地震臺（以下稱易縣臺）隸屬于保定地震監(jiān)測中心站，始建于1974年，位于易縣縣城西北方向龍里華村，距縣城15 km。臺站位于太行山隆起帶與華北平原拗陷區(qū)的分界部位，臺址處于太行山北段低山丘陵地區(qū)，東鄰紫荊關(guān)斷裂帶。臺基巖性為中生代侵入閃長巖花崗基底巖與震旦系沉積巖和閃長巖小侵入體的接合部。地形變觀測方法包括地傾斜、地應變、臺站跨斷層測量等。易縣臺是形變觀測綜合的國家基本臺。

易縣臺 D SQ 水管傾斜儀（以下簡稱水管儀）1995年正式投入觀測，2001年7月依托首都圈改造項目對觀測系統(tǒng)通過“九五”并網(wǎng)改造專項進行了改造，2014年10月由武漢地震科學儀器研究院有限公司將儀器升級改造成“十五”儀器。易縣臺水管儀共有3個觀測分量，其中 NS 分量方位角為0°00.00′，基線長10.15 m；EW 分量方位角為263°44.52′，基線長15.76 m；NE 分量方位角為233°46.63′，基線長15.30 m。很多學者對易縣臺水管儀觀測數(shù)據(jù)進行了分析與研究，1998年張北6.2級地震前易縣臺水管儀 NS 分量出現(xiàn)異常，在震前5天異常變化明顯[1]；2004年5月山西平遙4.0級地震前1—2天易縣臺水管儀 NE 向出現(xiàn)明顯的固體潮畸變[2]；2008年汶川8.0級大地震前水管儀 EW 分量曲線在2008年4—5月間整體呈現(xiàn)上升—轉(zhuǎn)折下行—急速回升形態(tài)，地面由 E 傾轉(zhuǎn)為 W 傾[3]。自2001年正式觀測以來，水管儀 NE 分量觀測數(shù)據(jù)連續(xù)，記錄曲線清晰， M2波潮汐因子、觀測精度較高，同震響應良好，每年年初 NE 分量觀測曲線呈 N 傾變化，5月轉(zhuǎn)為 E 傾變化，7、8月受降雨干擾影響觀測曲線出現(xiàn)不同程度的“鼓包”現(xiàn)象，特別是2016年和2021年尤為突出，直至11月后轉(zhuǎn)為 N 傾變化。然而易縣臺水管儀 NE 分量自2022年2月14日起觀測數(shù)據(jù)曲線持續(xù)出現(xiàn)趨勢下降 E 傾變化，年變觀測曲線周期變化不同于往年，本文針對易縣臺水管儀 NE 分量加速 E 傾短臨異常變化進行一系列的核實調(diào)查及跟蹤分析。

1 觀測山洞與周圍環(huán)境調(diào)查

易縣臺觀測山洞（圖1）建于1950年，人工開鑿，位于易縣臺臺址正南800 m 華蓋山山腰，屬于軍隊備用山洞，高程100 m。山洞大約呈北東東走向，山洞洞深長約200 m，共有7個硐室，在距洞口150余米處分岔為2條路，兩條路夾角約為90° ， 在山洞除離硐室最近的 NWS 硐室外，其他6個硐室都是與華蓋山山脈走向基本一致，山洞的洞寬不等，約在3.0—3.5 m 之間，高度在1.8—2.0 m 之間。山洞頂部覆蓋層厚約40 m，洞溫在（13.0±0.3）℃ ， 保溫效果好。山洞巖性底部為閃長巖，頂部為中厚層夾薄層石灰?guī)r和頁巖，臺址山洞巖石完整、致密，山洞洞口不在風口上，無滑坡及泥石流發(fā)生的可能，符合《地震臺站建設(shè)規(guī)范地形變臺站第1部分：洞室地傾斜和地應變臺站》的要求。

易縣臺山洞周圍的環(huán)境干擾源比較復雜，包括采石場、抽水井、水泥廠及水庫。其中巨石礦業(yè)有限公司位于山洞東偏北3.3 km，2016年至今主要以采石為主，2017年4—5月因大量采石荷載變化致使易縣臺水管儀 NS 分量出現(xiàn)短期加速 S 傾異常變化，后恢復正常；山洞周圍6 km 內(nèi)有3口正常使用的抽水井，綜合易縣臺歷年各類形變觀測曲線分析可知，抽水干擾對鉆孔應變儀器觀測數(shù)據(jù)影響較大，對傾斜觀測儀器觀測數(shù)據(jù)影響較小[4]；易縣臺站西南方向9 km 處的京蘭水泥廠，2009年運行至今，無異常出現(xiàn)；易縣臺臺站周圍共有2座水庫，一座是處在易縣臺山洞西北方向5.06 km 的旺隆水庫，另一座是處在易縣臺山洞東北方向7.62 km 的黃蒿水庫，異常核實過程中發(fā)現(xiàn)：旺隆水庫自2022年2月14日起因河北省“十四五”時期水庫除險加固和運行管護工作要求進行除險加固工程，每日持續(xù)泄流，泄流的水主要自北易水河分支流向梁各莊，在百益泉分支流向下游，匯入護城河（圖2）。直到4月6 日全部泄流完成，4月10日開始動工修防滲網(wǎng)。此次泄流總量達到279.50萬 m3，泄流途徑地距易縣臺山洞大約2 km，泄流造成易縣臺山洞西北荷載減少，使得該區(qū)域的應力場發(fā)生變化。據(jù)張磊等[5]研究表明，當水庫載荷突然發(fā)生大幅變化時，破壞了水庫與觀測點之間的相對平衡，造成觀測點一定程度的位移變化，從而影響該點的地傾斜觀測。為進一步驗證此次水管儀加速 E 傾變化是否是水庫泄流所致，將采用荷載模型進行定量計算分析。

2 水管儀 NE 分量異常及水庫泄流影響

2.1 異常概述

2022年2月14日易縣臺水管儀 NE 分量觀測數(shù)據(jù)曲線持續(xù)出現(xiàn)趨勢下降 E 傾變化，年變觀測曲線不同于往年。2022年2月14日—5月7 日水管儀 NE 分量觀測數(shù)據(jù)異常變化幅度達到0.29″，同硐室易縣臺伸縮儀 NE 分量觀測數(shù)據(jù)同一時段內(nèi)出現(xiàn)明顯的壓性變化，異常變化幅度比水管儀 NE 分量變化幅度?。▓D3）。2022年5月6—7 日依據(jù)中國地震局監(jiān)測預報司編著的《地震前兆異常落實工作指南》進行異常核實排查：①通過對觀測儀器、數(shù)采、通訊、外線路等設(shè)備進行逐一檢查，均未出現(xiàn)異常；2022年5月7 日進行水管儀儀器標定，結(jié)果為0.59109×10?3″，符合規(guī)范要求，標定結(jié)果確認水管儀觀測儀器運行正常；②通過對山洞硐室內(nèi)壁、地面的檢查，未發(fā)現(xiàn)積水、漏水現(xiàn)象；③通過對氣象因素對比分析，未出現(xiàn)影響因素；④未有人員在異常出現(xiàn)時對水管儀進行破壞、維修、調(diào)零，排除人為干擾；⑤現(xiàn)場環(huán)境干擾核實中，未發(fā)現(xiàn)新增水庫、水井，只有旺隆水庫2—4月放水泄流造成的區(qū)域荷載變化，其他均未出現(xiàn)任何變化。因此，我們只需針對水庫泄流荷載進行異常分析。

2.2 旺隆水庫泄流影響分析

2.2.1 水庫泄流庫容情況

旺隆水庫是一座中型水庫，地處易縣縣城西13 km，位于河北省易縣北易水上游旺隆溝。1958年興建，1960年竣工?？刂屏饔蛎娣e37 km2，總庫容1275萬 m3，有大壩、放水洞、溢洪道等工程。本流域50%年徑流量510萬 m3，原興利庫容486萬 m3，1990年在溢洪道興建橡膠壩后，興利庫容達886萬 m3[6]。

2.2.2 水庫泄流日庫容對水管儀 NE 分量觀測影響的定性分析

旺隆水庫2022年2月14日起開始日泄流，時間持續(xù)52天，最多日泄流容量達29萬 m3，速度達到3 m3/s。2022年2月14—28日，15天內(nèi)泄流總?cè)萘扛哌_158.48萬 m3，3月1 日—4月6 日，37天泄流總?cè)萘窟_到121.02萬 m3。將日泄流容量與水管儀 NE 分量日變化值進行對比（圖4），可以清晰地看出，2月17日和18日旺隆水庫日泄流容量達到29萬 m3和23.32萬 m3時，對應的水管儀 NE 分量日變化幅度卻只有?4.08×10?3″和?2.84×10?3″；而當3月12日水管儀 NE 分量日變化幅度為?10.38×10?3″時，而當日旺隆水庫卻未進行泄流。因此，日庫容變化總量與水管 NE 分量數(shù)據(jù)變化關(guān)系一致性較差。但是，旺隆水庫2022年2月起52天加速泄流庫容總量達到300萬 m3，荷載變化或許是引起水管儀 NE 分量加速 E 傾的主要原因，考慮通過荷載模擬定量分析。

2.2.3 水庫泄流引起水管儀 NE 分量加速 E 傾的理論分析

水庫卸載都是由于水位變化產(chǎn)生的，可以用三維集中荷載模型進行分析[7]。作為一級近似，可用三維無限半空間邊界上作用集中荷載模型研究這種荷載 F 變化，其位移解：

為計算荷載對傾斜觀測影響，只需對 x 水平方向求導即可，即：

式中，F(xiàn) 為水庫泄流造成的荷載集中力，可根據(jù) F=mg 得出 F=2.79×1010 N（水庫變化量 m=2.79×109 kg， g≈10 N/kg），x、y、z 為易縣臺山洞到旺隆水庫邊緣的方位距離（x=5 km，y=0，z≈100 m）；G、λ為拉梅常數(shù)，可根據(jù)公式 G=E/[2×（1+ν）]得出，其中水的泊松比ν=0.5，水的體壓縮楊氏模量 E ＝2.2×109 Pa[8]，由于水是液體只有一個拉梅常數(shù) G，得出 G=0.7×109 Pa。拉梅常數(shù)λ我們定義為0 Pa ，由公式（2）計算得到ex=1.27×10?7。依據(jù)邱澤華[7]參考傾斜干擾的限定指標1.455×10?8可知，此次水庫泄流造成的荷載影響已經(jīng)超過近1個量級。因此，易縣臺水管儀 NE 分量加速 E 傾異常變化主要與水庫泄流造成荷載變化有關(guān)。

2.2.4 同硐室內(nèi)地傾斜實際變化與傾斜荷載理論值的對比

2022年2月14日—4月6 日旺隆水庫泄流庫容量為279.5萬 m3，依據(jù)公式（2）計算得到傾斜變化理論值為26.19×10?3 s，即旺隆水庫2022年2—4月造成易縣臺山洞傾斜儀器水平方向變化量為26.19 ms。結(jié)合易縣臺山洞2022年2月14日—4月6 日垂直擺、水管儀觀測數(shù)據(jù)曲線變化（圖5）發(fā)現(xiàn)：①垂直擺傾斜儀2022年2月14日—4月6 日 NS 分量觀測數(shù)據(jù)曲線出現(xiàn)N 傾，變化幅度為33.70 ms；EW 分量觀測數(shù)據(jù)曲線出現(xiàn) E 傾，變化幅度為12.00 ms；②水管儀2022年2月14日—4月6 日 NS 分量觀測數(shù)據(jù)曲線出現(xiàn)S 傾，變化幅度為661.52 ms（3月19日調(diào)零干擾影響，變化幅度增大）；EW 分量觀測數(shù)據(jù)曲線出現(xiàn) E 傾，變化幅度為65.53 ms；NE 分量觀測數(shù)據(jù)曲線出現(xiàn)加速 E 傾，變化幅度為155.29 ms；③垂直擺兩分量觀測數(shù)據(jù)變化幅度與旺隆水庫泄流荷載造成的傾斜理論變化值（26.19 ms）相差較少；④垂直擺 EW 分量觀測數(shù)據(jù)曲線與水管儀 EW 分量、NE 分量觀測數(shù)據(jù)曲線形態(tài)變化一致，都為 E 傾（表1）。

3 結(jié)語

易縣臺水管儀 NE 分量觀測數(shù)據(jù)2022年2月起加速 E 傾異常變化，經(jīng)過異常核實排查與論證，我們可以看出：

（1）2022年水管 NE 分量傾斜觀測數(shù)據(jù)發(fā)生了 E 傾異常變化，依據(jù)《地震前兆異常落實工作指南》異常核實要求，通過對觀測儀器、環(huán)境干擾、網(wǎng)絡(luò)通訊、人為原因等因素核實，判定得出：此次水管儀E傾異常變化主要與旺隆水庫泄流造成荷載影響有關(guān)，這是因為水庫載荷突然發(fā)生大幅變化，破壞了水庫與觀測點之間的相對平衡，造成觀測點一定程度的位移變化，從而影響到地傾斜結(jié)果；

（2）利用邱澤華[7]三維集中荷載模型進行模擬計算得出：2022年2月14日—4月6 日旺隆水庫共泄流總庫容量達到279.5萬 m3，造成地傾斜變化模擬理論值1.27×10?7，高于傾斜干擾的限定指標1.455×10?8一個量級，此是造成易縣臺水管儀 NE 分量觀測數(shù)據(jù)加速異常變化的主要因素；

（3）通過對易縣臺山洞內(nèi)水管儀、垂直擺兩套地傾斜儀器2022年2月14日—4月6 日觀測數(shù)據(jù)實際異常幅度變化情況與水庫泄流造成荷載傾斜理論值進行對比分析得出：①垂直擺傾斜儀兩分量觀測數(shù)據(jù)變化幅度與旺隆水庫泄流荷載造成的傾斜理論變化值最接近；②垂直擺 EW 分量觀測數(shù)據(jù)曲線與水管儀 EW 分量、NE 分量觀測數(shù)據(jù)曲線異常形態(tài)變化一致，都為 E 傾，這是因為旺隆水庫位于易縣臺形變山洞西北方向，應力主方向為東西向，因而此次泄流造成荷載應力對水管儀 EW 分量、NE 分量及垂直擺 EW 分量觀測數(shù)據(jù)影響較大。

依據(jù)以上異常分析，我們得出旺隆水庫泄流荷載是造成易縣臺2022年2月水管儀 NE 分量觀測數(shù)據(jù)加速 E 傾的主要因素。

參考文獻

［1］陳開坊，李占錄.張北6.2級強震前易縣臺地形變臨震異常[J].華北地震科學，1999，17（4）：70-71

ChenKF， LiZL. ImpendingearthquakeanomalyoftopographicdeformationatYixianstationbeforeZhangbeiM6.2 strong earthquake[J]. North China Earthquake Sciences，1999，17（4）：70-71

［2］陳建國，張娜，王秀文，等.形變數(shù)字化資料對山西地區(qū)地震前兆反應研究[J].山西地震，2009（1）：1-8

Chen J G，Zhang N，Wang X W， et al. Study on the response of digital deformation data to earthquake precursor of Shanxi[J]. Earthquake Research in Shanxi，2009（1）：1-8

［3］蘇樹朋，李津津，張勤，等.易縣地震臺形變資料對汶川8.0級大地震的反應[J].華北地震科學，2009，27（2）：54-58

Su S P，Li J J，Zhang Q，et al. The response of crustal deformation data in Yixian seismic station to Wenchuan M8.0 earthquake[J]. North China Earthquake Sciences，2009，27（2）：54-58

［4］張娜，龔燕民，李博.洞體應變觀測中荷載干擾實例分析[J].地震地磁觀測與研究，2017，38（4）：153-157

Zhang N，GongY M，Li B. Case analysis of load interference of cave strain observation[J]. Seismological and Geomagnetic Observation and Research，2017，38（4）：153-157

［5］張磊，邵永新，曹井全，等.薊縣臺水管傾斜儀 NS 向加速北傾異常分析[J].大地測量與地球動力學，2016，36（4）：372-376

Zhang L，Shao Y X，Cao J Q，et al. Analysis of acceleration incline to north direction of water tube tilt-meter at JixianSeismostation[J]. Journal of Geodesy and Geodynamics，2016，36（4）：372-376

［6］趙連勛.旺隆水庫成了聚寶盆[J].河北水利，1994（4）：22-23

Zhao L X. Wanglong Reservoir become a cornucopia[J]. Hebei Water Resources，1994（4）：22-23

［7］邱澤華.傾斜觀測載荷干擾源最小“安靜”距離理論分析[J].地震，2004，24（4）：80-86

Qiu Z H. Theoretical study on the least“ Quiet” distance between borehole tilt observatory site and surficial mechanical disturbance sources[J]. Earthquake，2004，24（4）：80-86

［8］劉瀾波，鄭香媛.井水固體潮分析結(jié)果及其在地震預報中的應用[J].地震，1985（1）：7-12

Liu L B，Zheng X Y. The analysis results of water-well earth tides and their applications to earthquake prediction[J]. Earthquake，1985（1）：7-12