山西靜樂觀測井水位干擾分析

劉俊芳呂 芳郭 宇高文玉郭寶仁

1）中國山西035403定襄地震臺

2）中國太原030021山西省地震局

3）中國山西030025太原大陸裂谷動力學國家野外觀測站

山西靜樂觀測井水位干擾分析

劉俊芳1），3）呂 芳2），3）郭 宇1），3）高文玉1），3）郭寶仁1），3）

1）中國山西035403定襄地震臺

2）中國太原030021山西省地震局

3）中國山西030025太原大陸裂谷動力學國家野外觀測站

山西靜樂觀測井具有含水層系統對地殼應力—應變響應靈敏的優勢，其水位變化對應地震，可間接反映地殼應力—應變特征。2015年7月以來，靜樂井水位快速上升，對觀測環境、儀器及人為干擾進行分析，并進行同層地下水開采及抽水試驗，分析該井水位異常原因，結果發現，降雨量、抽水試驗對靜樂井水位變化影響不大，而生產用同層地下水開采量減小是該井水位變化的主要原因。

水位；干擾分析；地下水開采；調查；異常

0 引言

地下流體觀測可以獲得地球深部信息，地下水位異常可能是最直接的宏觀短臨地震前兆，地下水位觀測已廣泛用于地震監測，是開展地震監測的重要手段。降水及人為開采造成水位升降，并因此影響水位觀測是不可避免的。調查研究觀測井水位異常及干擾因素具有積極意義。多個研究者對山西地下流體觀測井水位異常進行研究，孟彩菊（2010）認為，水位變化影響短水準高差，二者相關性及年最大變化量的異常與山西地震帶內發生的地震有較好對應。沈曉松等（2012）研究太原井水位異常，認為與煤窯關閉、汾河清水復流工程、關井壓采、黃河水大量引用具有直接關系。楊曉芳等（2013）研究認為，地震波引起含水層介質及孔隙壓力變化，導致觀測井內上下層水體進入井孔，水溫表現為臺階變化。近年來山西靜樂井水位變化明顯，文中對引起井水位變化的各種干擾因素進行分析，探討水位變化與地下水開采的關系。

1 觀測井概況

山西靜樂井位于晉西北汾河上游，地質構造處于呂梁山斷塊隆起區、靜樂盆地東緣、東碾河斷裂帶上，出露基巖為寒武系、奧陶系碳酸鹽巖系，上覆第四系砂礫石層。靜樂井深362.92 m，觀測層為中奧陶灰巖裂隙承壓水，含水層巖石溶洞發育，井區處于強徑流區，地下水補給為來自北、南山區基巖裸露區的大氣降雨補給和東碾河河水沿斷裂破碎帶入滲補給。

多年觀測資料表明，靜樂井水位干擾因素主要有觀測系統、自然環境和人為干擾。其中：觀測系統干擾包括儀器故障、電源干擾、路由器死機、通信單元故障等；自然環境干擾包括雷電、刮風和降雨干擾等；人為干擾包括模擬換紙干擾和校測、儀器更換等。

2 井水位正常動態特征

2.1 趨勢變化特征

從1983—2015年靜樂井水位日均值曲線（圖1）可以看出，該井水位長趨勢變化呈現4個階段：1983—1998年為高值期，1999—2006年為低值期，2007—2011年為高值期，2012—2015年為低值期。2012年以來的趨勢性轉折變化主要由周圍工業用井抽水干擾影響所致。

圖1　靜樂水位長趨勢曲線Fig.1 Long-term trend curve of Jingle well water level

2.2 年動態特征

多年的觀測資料顯示，靜樂井水位具有較規律的年變動態特征，表現為兩峰兩谷型，高水位期在每年5月、9月，低水位期為2月、7月。文中選擇井水位趨勢變化的4個階段的部分年份進行分析，靜樂井水位受氣壓、固體潮和降雨影響。由圖2可知，降雨量大的年份靜樂井水位變幅大，可見降雨是影響該井水位動態變化的主要因素。

圖2　靜樂水位日均值年變化Fig.2 Annual change of daily average water level of Jingle well

3 井水位異常變化

3.1 井水位快速上升

靜樂井水位2007年7月23日—10月3日累計上升1.89 m，之后持續高值（圖1）。2013年—2015年6月受生產用同層地下水開采影響，水位處于低值，2015年7月水位突然快速上升，且已超過多年平均值（圖2），是否為地震前兆異常亦或受同層地下水開采影響所致，在此進行分析。

3.2 井水位地下應力響應分析排除前兆異常

靜樂井水位可記錄國內外6.5級以上地震的同震響應。水位對地震的響應為震蕩型水震波，震蕩幅度大小和周期長短與井—含水層系統特征有關。水震波主要由地震面波引起，其響應幅度與震級正相關（圖3）。

圖3　靜樂水位同震震蕩曲線（a）2015年4月25日尼泊爾8.1級地震；（b）2015年7月3日新疆和田6.5級地震；（c）山西朔州水位（2015年1月—5月）與華北地震對應關系；（d）靜樂水位（2015年4月1日—18日）與華北地震對應關系Fig.3 The magnitudes of earthquakes concussion type graph of Jingle well water level

根據呂芳（2015）等的研究，靜樂井對地殼應力—應變的響應能力取決于井—含水層系統的特征參數，如含水層厚度（d）、滲透系數（K）、導水系數（T）、井孔半徑（rw）等。根據Cooper等（1965）、張昭棟等（1999）的研究，井—含水層系統對地殼應力—應變的響應能力與導水系數、滲透參數、固有周期、有效水柱高度成正比，與井孔半徑成反比，與含水層巖性關系密切，一般灰巖含水層觀測井對地殼應力—應變的響應能力大于其他巖性的含水層。靜樂井觀測層屬于中奧陶灰巖裂隙承壓水，地殼應力—應變響應能力大于其他巖性的含水層。因此，靜樂井能較為靈敏地反應地下應力變化的特征。

由圖3（c）可見，2015年4月15日15：39內蒙古阿拉善左旗MS6.1地震前（4月1日—15日），朔州觀測井（距震中約701.04 km）水位持續下降約0.414 m，變幅3.54%，達到異常變化范圍；震后水位持續下降，至4月24日降到最低，共下降0.772 m，變幅6.2%，之后逐漸回升至正常狀態。靜樂井在此次MS6.1地震震中東南方向，距震中約657.71 km。由圖3（d）可見，靜樂井水位在該地震前（4月12日—15日15：38）上升0.053 m，震后（4月15日15：43—16日0：51）急速下降0.049 m，變化幅度均為0.5%。靜樂井水位2015年7—9月快速上升約0.93 m，但無地震對應。

綜上所述，靜樂井水位2015年7月快速上升，不是地震前兆異常，應為干擾所致。

4 井水位干擾分析

靜樂井水位呈季節性趨勢變化，引起水位測值產生變化原因多樣，結合井水位動態特征，從觀測環境、觀測儀器及人為干擾、生產用同層地下水開采及相關抽水試驗等，分析2015年7月靜樂井水位快速上升干擾因素。

4.1 觀測環境

（1）靜樂井水位短期動態變化與觀測環境間接相關，距觀測井約200 m的東碾河河水雨季上漲，枯水季節河水流量減小，水位回落，河水漲落對水位變化具有荷載作用。

（2）靜樂井水位受降雨影響，隨降雨補給趨勢上升。李自紅等（2011，2015）和張淑亮等（2009）分析認為，上升幅度與雨量正相關，由井—含水層特征及井區環境條件分析，降雨是導致水位快速上升的主要原因。而靜樂井水位與降雨之間，并不是每次強降雨都會造成水位大幅度上升。水位上升伴隨強降雨，只有在應力場快速積累過程中，降雨的調制作用才會造成水位高值異常（圖4）。2015年降雨量與往年同期相比無明顯變化，因此，7月以來水位的快速上升受降雨量影響的可能性較小。通過分析降雨、水溫與水位的變化（圖4），排除降雨對靜樂井水位的影響。

圖4　靜樂井水位、水溫日均值與降雨量對比Fig.4 The comparison of daily average of water level and water temperature at Jingle well with rainfall

4.2 觀測儀器及人為干擾

靜樂水位出現不明原因變化，于2015年5月通過同層觀測數字水位儀（SWY-Ⅱ）對比觀測，發現兩套儀器的水位曲線變化趨勢一致（圖5），排除儀器及人為干擾。

圖5　靜樂井同層水位同步觀測儀器對比Fig.5 The comparison of synchronized observation instruments with the same layer water level of Jingle well

4.3同層地下水開采

周圍開采環境對觀測井的影響，間接反映了井水位的動態變化。用于生產的山西天柱山化工廠2個同層開采井在測點允許范圍內［圖6（a）］，1號井位于靜樂井上游1.5 km處，2號井在2 km處。1號井夜間抽水，2號井全天抽水，影響了該地區地下水分布，改變了靜樂井水位的年變動態。根據靜樂井水位資料，2012年以來水位持續低值；2013年—2015年7月受化工廠抽水干擾，每年1—2月水位上升，其他時段水位處于起伏變化態勢，整體處于低值異常。

當1號井和2號井抽水作業時，靜樂井水位下降，且與水溫的臺階變化時間相對應［圖6（b）］。當同層開采井大量抽水時，導致不同層位地下水進行交換，打破水動態平衡狀態。受不同層位冷熱水混入等影響，觀測層水溫發生變化。靜樂井水溫隨井深增加而下降，水位下降或上升，導致地下水循環速度加快，引起水溫升高或下降。觀測井上游的同層地下水開采，抽走一部分低溫地下水，水循環要達到平衡，部分地表水通過滲流補給，含水層溫度通過熱傳導，溫度上升；遠處深層水的補給，溫度下降。因此，水位下降，水溫則呈現“上升—下降”起伏變化，屬于非構造原因引起。

圖6　靜樂井水位受抽水干擾影響曲線（a）測點環境；（b）2015年6月1日—15日觀測井水位、水溫分鐘值曲線Fig.6 The pumping interference of Jingle well water level

經調查，山西天柱山化工有限公司生產用井2012年9月至2014年8月為試生產階段，生產用水以地表水為主，地下水間斷使用；2014年9月以后為正式生產階段，生產用水以地表水為主，地下水為輔，地表水水量不足或渾濁，使用地下水；2015年8月中旬，地表水經改造滿足生產用水，地下水作為生活用水。2014年8月至2015年10月27日，山西天柱山化工有限公司生產用井抽水統計見表1。

表1　靜樂井測點周圍地下水生產用井抽水統計Table 1 The conditions of pumping groundwater with the production wells surrounding Jingle well

由表1可知，2015年7月以來，地下水抽水減少約6 000 m3，靜樂井水位快速上升，水溫也同步上升（圖4）。進一步說明，工業用井抽水對靜樂井水位存在干擾。

4.4 抽水試驗

對化工廠進行二級水源地影響半徑測定，對實驗用井G1井于2015年5月26日至8月8日進行抽水試驗，該井位于靜樂井下游偏西約1.8 km的水源地［圖6（a）］。在試驗過程中，G1井抽水量較大（表2），但因G1井所處位置，且在2015年7—8月G1井抽水期間，靜樂井水位呈上升趨勢（圖2），認為該井抽水對靜樂井水位變化影響不大。

表2　靜樂井測點周圍地下水抽水統計Table 2 The conditions of pumping groundwater surrounding Jingle well

5 討論

山西靜樂井因同層井工業用水的抽水干擾，導致地下水位分布發生變化，水位觀測資料的正常年變動態發生改變，降雨量、抽水試驗對靜樂井水位變化影響不大，而生產用同層地下水開采量減小是該井水位變化的主要原因。

李自紅，范雪芳，呂芳，等.山西省2011年度地震趨勢研究報告［R］.山西省地震局，2010：36-37.

李自紅，范雪芳，呂芳，等.山西省2015年度會商報告［R］.山西省地震局，2014：32-38.

呂芳，張淑亮，胡玉良，等.山西靜樂井水位異常調查與分析［J］.山西地震，2013，3：28-31.

孟彩菊.太原地震臺短水準受水位影響分析及異常識別［J］.地震地磁觀測與研究，2010，31（1）：77-83.

沈曉松，靳玉貞，孟彩菊，等.太原地震臺水位異常分析［J］.地震地磁觀測與研究，2012，33（5/6）：178-182.

楊曉芳，崔勇，張濤，等.新疆地下流體動態響應特征研究［J］.山西地震，2013，（2）：22-28.

張淑亮，范雪芳，李斌.山西靜樂、孝義井水位異常調查報告［R］.山西省地震局，2009：109-114.

張昭棟，遲鎮樂，陳會民，等.水井含水層系統與水位觀測系統對固體潮與地震波的響應［J］.地震學報，1999，10（2）：171-181.

Cooper H H.The response of well-aquifer systems to seismic waves ［J］.Geophys Res，1965，70（16）：3 915-3 926.

Abstract

The observation well of Jingle in Shanxi Province has an advantage which aquifer systems are of sensitive response to the crustal stress-strain.The corresponding to the earthquake of Jingle well water level can reflected the crustal stress-strain characteristics.The abnormality cause of the well water level is found through the analysis of the observational environment，equipment，human disturbance and the same layer of groundwater mining and pumping test.The final results indicated that the rapid rise of water level of Jingle well since July 2015 is mainly caused by the reduction of groundwater exploitation in the well from the same aquifer.

Analysis of water level observation interference of Jingle well in Shanxi Province

Liu Junfang1），3），Lv Fang2），3），Guo Yu1），3），Gao Wenyu1），3）and Guo Baoren1），3）
1） Dingxiang Seismic Station，Shanxi Province 035403，China
2） Earthquake Administration of Shanxi Province，Taiyuan 030021，China
3） National Continental Rift Valley Dynamics Observatory of Taiyuan，Shanxi Province 030025，China

water level，interference analysis，groundwater exploitation，investigation，anomaly

10.3969/j.issn.1003-3246.2016.03.005

劉俊芳（1986—），女，工程師，主要從事臺站前兆監測分析工作