利用水文地球化學和水文物理學方法分析定襄七巖泉水氡異常

劉俊芳 李惠玲 范雪芳 郭寶仁 郭 宇 高文玉 劉金柱

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

2）中國山西 034200 忻州地震監測中心站

0 引言

對于地震地下流體觀測，可利用水文地球化學方法進行數據分析，如：胡小靜等（2018）利用水文地球化學方法，研究云南江川漁村井的地下水循環特征和補給來源；田雷等（2018）通過分析滇17 井水化學測項，發現觀測數據變化由地下應力改變所致；王云等（2018）對位于滇東南楔形構造區的觀測點進行水化學和氣體化學特征分析，判斷不同泉點的水巖作用、循環特征、逸出氣成因；楊靜等（2019）利用水文地球化學方法，分析了夏縣溫泉水補給來源、溫泉類型。

山西定襄七巖泉水氡觀測始于1981 年6 月，日常采樣為混合泉水，長期觀測的背景資料有待完善。2008 年以來，該泉水氡測值呈上升趨勢，在水源點、采樣點及其周圍采集水樣并分析離子組分和同位素特征，判斷水質類型，分析水—巖平衡狀態，并根據氫氧同位素特征判斷泉水補給來源，利用水文地球化學和水文物理學分析方法，分析該泉水氡測值異常上升原因及降水補給來源。

1 取樣環境

定襄七巖泉為承壓上升泉，屬寒武系灰巖喀斯特裂隙溶洞水。泉點地處定襄縣茶房口村以南約500 m 的七巖山山谷西側，流體取樣點靠近系舟山北麓斷裂，附近石灰巖、變質巖等裂隙節理及溶洞較為發育，其中賦存的裂隙溶洞水在斷層附近呈群泉溢出。泉口處原有6—7 個出水泉眼，被沖積物所覆蓋，水氡觀測日常采樣為混合泉水，水溫約10.0 ℃，高于當地年平均氣溫（約8.9 ℃）。由地貌、含水層及水溫條件判斷，地下水循環深度相對較大。

因水氡觀測房建造已久，取水點上游幾十米處陶瓷水管破裂，水流暴露在外。2014年8 月，改造為PVC 引水管道，之后水流量增大，水氡測值下降12 Bq/L（改造前后水氡采樣點在同一位置，水樣為混合泉水）。異常核實認為，氡值下降由觀測房內取樣出水口密封不當所致。此后，水氡觀測動態特征趨于穩定。

2 水氡異常及分析

2008 年以來，七巖泉水氡測值呈上升趨勢變化，2019 年4 月，在泉水上游到下游布設定襄1#—5#采樣點（圖1），依次進行現場取樣，并將樣品送原中國地震局地殼應力研究所地殼動力學重點實驗室測試水化學離子成分和氫氧同位素成分，并采用水文地球化學和物理學方法進行分析，判斷泉水補給來源，核實異常原因。

圖1 現場采樣點分布Fig.1 Distribution of sampling points on site

2.1 水文地球化學分析

采用水化學類型、離子成分、水—巖反應、氫氧同位素等方法，對七巖泉水樣進行水文地球化學分析，判斷七巖泉補給來源，分析水氡高值異常原因。

2.1.1 水巖化學平衡特征。水巖平衡分析用于體系開放與封閉的判斷，時間及運移過程的判斷（伍劍波，2013）。地下水中含F-、Cl-、NO3-、Na+、K+、Ca2+、Mg2+、SO42+等離子，承載著流體來源、運移和巖石圈信息（張磊等，2016）。而水的電導率與水中所含的溶質和雜質呈正相關性（陳愛華，2013）。

（1）采樣點水樣離子成分。據有資質的檢測機構對送檢水樣的測試結果，定襄七巖泉5 個采樣點水樣離子成分主要有F-、Cl-、NO3-、Na+、K+、Ca2+、Mg2+、SO42+，見圖2。

由圖2 可知：采樣點定襄3#各離子成分較其他采樣點均偏高，其中Cl-含量較高，表明3#泉水流程較遠，含溶質和雜質多；定襄泉1#、2#、4#、5#采樣點水樣化學組分接近，3#采樣點水樣離子成分明顯不同，其中前者水質表現為Ca-Mg-HCO3型，3#采樣點為Ca-Mg-HCO3-SO4型，表明3#與其他采樣點水樣補給來源不同。

圖2 定襄七巖泉不同采樣點水樣離子成分Fig.2 The composition of different sampling points in Dingxiang Qiyan spring

（2）采樣點水質分析。從piper 三線圖、durov 圖和Na-K-Mg 三角圖進行水質分析，定襄泉各采樣點水樣水質分析結果見圖3（圖中各化學元素均未添加離子符號）。由圖3（a）、（b）可見，定襄七巖泉為Ca-Mg-HCO3型水，表明地表水或降水對泉水影響較大。由圖3（c）可見，定襄七巖泉泉水屬于“未成熟水”，且為淺層地下水，主要接受大氣降水補給，循環周期相對較快，水—巖之間尚未達到離子平衡狀態，溶解作用仍在進行。

圖3 定襄泉及其周圍水質（a）piper 三線分析；（b）durov 分析；（c）Na-K-Mg 三角分析Fig.3 Water quality of Dingxiang Qiyan spring and its surroundings

2.1.2 氫氧同位素特征。氫氧同位素方法是通過大氣降水，確定補給水源的水文地球化學方法。經實驗室測試，定襄七巖泉δD 范圍在-73.57‰—-65.85‰，δ18O 范圍在-10.61‰—-9.50‰，結果見圖4。由上述分析可知，七巖泉水Ca2+占比較高，屬Ca-Mg-HCO3型低礦化度水質，是典型的大氣水與巖石之間第一階段作用的反映，顯示了淺層水文循環特征。由圖4 可見：①七巖泉水氫氧同位素比值沿全球大氣降水線分布，結合賈振興等（2015）的研究結果，其比值亦靠近太原地區大氣降水線。泉水補給源為大氣降水的淺層地下水，水—巖反應程度較低，循環深度不大；3#與其他采樣點水源不同，體現在：3#采樣點出現在2 條大氣降水線上，為大氣降水，而其他采樣點位于2 條大氣降水線附近，為大氣降水經巖石徑流積累水。

圖4 定襄泉氫氧同位素分析Fig.4 The analysis of hydrogen and oxygen isotope of Dingxiang spring

2.2 水文物理學分析

采用降雨干擾分析與地震活動分析的方法，對定襄七巖泉水氡進行高值異常分析，研判降雨對水氡動態的影響程度、引起水氡高值趨勢性異常的可能性及水氡高值異常是否對地震活動具有指示意義。

2.2.1 降雨干擾。定襄七巖泉為上升泉，降水量增大，泉水流量增加，水氡測值則出現明顯下降變化，分析認為氡值下降與降水量及泉水流量增大有直接關系。以2016 年夏季降水干擾為例，探討降雨量對定襄七巖泉氡含量的影響。據調查，2016 年7 月—9 月定襄當地降雨量累計314.7 mm，七巖泉水流量增大，氡含量短期內出現大幅度下降，降幅大于50%（圖5）。實驗檢測雨水中氡的含量為0，可見降雨會降低測點泉水樣品氡含量。整理統計定襄當地2010—2019 年降雨量及對應的七巖泉氡值變化，結果見表1。

圖5 2016 年7—9 月定襄水氡、流量及降雨量日值曲線Fig.5 The curves for data analysis of water radon from July to September in 2016

表1 定襄臺降雨量與水氡測值變化對比Table 1 Comparison of rainfall and water radon of Dingxiang Seismic Station

由表1 可知，受連續降雨影響，水氡測值先上升變化，當降雨量達到一定閾值時出現下降變化。其中，歷年降雨量累計均值為116.09 mm，當累計連續降雨量超過一定閾值（116 mm）時，水氡測值明顯下降。多年觀測發現，降雨僅造成氡值的短期變化，不會改變其年動態趨勢，由此排除降雨引起水氡高值趨勢性異常的可能。

2.2.2 地震活動。多年觀測發現，定襄七巖泉水氡具有映震異常，始測以來共出現4 組高值異常，前3 組高值異常與華北地區5 級以上地震有較好的對應關系，如：高文玉（2006）分析發現，定襄水氡高值異常與華北地震活動存在同源關系，反映了區域應力場的活動性；范雪芳等（2010）探討了定襄水氡、靜樂水位的異常特征，認為二者在華北和山西北部中強地震前均具有明顯的中期和中短期異常變化；呂芳等（2013）根據靜樂井的井—含水層系統對地殼應力—應變的響應能力較山西地區其他觀測井高及震例研究結果，認為靜樂井水位高值異常幅度與對應地震的震級及發震距離成正比，2007 年以來高值異常與觀測環境、降雨等無明顯關系，可能與區域應力場活動增強有關，其高值異常對山西北部中等地震具有指示意義；劉俊芳等（2016）分析認為，靜樂井能較為靈敏地反映地下應力變化特征，2015 年7 月以來靜樂井水位快速上升受降雨量、抽水試驗的影響較小，主要與生產用水的同層地下水開采量減少有關。根據歷史震例，結合定襄和靜樂當地降雨量統計結果，對比分析1983—2019 年定襄泉水氡和靜樂井水位的高值異常與華北地區中強地震的對應關系，結果見圖6、表2。

由圖6、表2 可知，定襄水氡相對高值異常狀態與靜樂井水位高值異常同步。水氡和水位歷年動態顯示，降雨干擾不會造成高值趨勢性異常；定襄水氡與靜樂水位變化動態相似，其高值異常與華北及晉冀蒙交界區中強地震有較好的對應關系。因此，2008 年至今，定襄水氡高值異?？赡芊从沉藚^域地殼應力變化。因泉點改造，2014 年七巖泉采樣點出水口裸露，造成水氡逸散，測值降低（郭宇等，2016）。

表2 定襄泉水氡、靜樂井水位觀測映震統計Table 2 Statistics of earthquake reflections of the radon observation in Dingxiang spring and water level observation in Jingle well

圖6 靜樂井水位、定襄泉水氡及其對應降雨月均值曲線對比Fig.6 Comparison of water level of Jingle well,water radon of Dingxiang spring,and their corresponding monthly mean value rainfall curves

3 結果與討論

對于定襄七巖泉2008 年以來出現的水氡測值上升趨勢變化，應用水文地球化學和水文物理學進行異常分析，研究認為：①水文地球化學分析：七巖泉泉水屬于“未成熟水”，主要補給源為大氣降水，循環周期相對較快，為大氣成因水，其水—巖反應程度較低，循環深度不大；在水化學類型、離子成分、水質分析、大氣降水線位置分布、采樣位置等方面，3#采樣點水源為大氣降水，其他采樣點水樣為大氣降水經巖石徑流積累的水；②水文物理學方法分析：降雨僅對水氡短臨動態造成干擾，年變動態不受影響；當累計連續降雨量超過一定閾值（116 mm）后，七巖泉水氡測值明顯下降。定襄水氡與靜樂水位變化動態相似，其高值異常與華北及晉冀蒙交界區中強地震有較好的對應關系。

通過綜合分析，發現定襄七巖泉水氡處于相對高值狀態與靜樂井水位2017 年9 月以來高值異常顯著同步，表明七巖泉水氡高值異常可能是區域構造應力變化的反映。定襄七巖泉水氡與靜樂井水位準同步高值現象是否反映了二者具有地下流體的同源關系，以及該準同步現象是否反映區域應力可能達到新的動態平衡，有待積累資料進一步研究。