衡水市咸水入侵問題研究綜述

鄧興 陳康 劉雪友

摘 要：咸水入侵問題是水文地質研究中的一個重要課題。本文以衡水市為例重點介紹了咸水成因、分布以及咸水入侵的機理，綜述了水質變化、水質預測以及咸淡水界面變化的研究現狀，總結出咸水入侵的應對措施，為地下水資源的合理開發利用和淡水資源的保護提供理論依據。

關鍵詞：咸水入侵;地下水水質;咸淡水界面

引言

咸水入侵指的是在人類活動影響下淡水水量和水位的降低，天然狀態下咸淡水之間的平衡被打破，導致咸淡水交界面向淡水一側推進，直到達到新的平衡狀態。

快速發展的社會經濟對于水資源的需求不斷增大，自然界水資源的有限性要求人們必須加強對水資源的管理[1]。衡水市位于華北平原的中部，地表水資源極度匱乏，地下淺部蘊含著大量的咸水，在咸水之下的深層水就是衡水市主要開采的淡水水源。自20世紀70年代以來，衡水市經濟發展迅速，工農業的快速發展導致需水量大量增加，20世紀80年代以后，由于降水量減少，為了滿足發展的需求，衡水市大量超采地下水，導致淺層淡水水質咸化和咸淡水界面下移，已經嚴重威脅到了衡水市的供水安全。因此，衡水市地下水水質以及咸淡水界面的變化成為人們急需研究的重要問題。

1 咸水成因及入侵機理

1.1 咸水成因

衡水市全區除深州北部、安平、饒陽西北部地區有全淡水區外，其它地區均有咸水體分布[2]，咸水體自西北到東南呈逐漸變厚的楔形體分布。咸水成因大致分為2種類型：海相原生咸水，第四紀發生多次海進，海相層中地下水化學成分與現代海水近似，說明海相地層沉積的同時，海水沉積于地下而形成高礦化咸水，在地下水形成過程中，原生海水被逐步改造，形成具有濱海特征的海湖積平原區咸水;陸相次生咸水，該類型的咸水是在氣候、古地理環境等因素影響下的地下水循環過程中形成的[3]。

1.2 咸水入侵機理

1.2.1 淺層咸水入侵淺層淡水

淺層地下水開采區大量開采淺層淡水，周圍咸水水頭高于淡水，打破了咸淡水之間原有的動態平衡，造成咸水區向淡水一側擴散，咸水區面積不斷擴大，形成咸水入侵。根據咸水入侵機理理論分析，衡水地區淺層水咸化現象由于大量開采淺層水而不斷加劇，但實際研究表明淺層咸水有淡化趨勢，并且淡化速度較快。周曉妮等[3]將所取衡水市淺層水樣品礦化度與多種離子化學成分含量間的相關關系進行分析，結合前人所做研究[4]發現1982—2006年衡水地區淺層地下水礦化度不斷減小，得出淺層咸水淡化的結論，結合相關研究通過分析認為，該地區淺層咸水淡化的原因主要是大量開采深層淡水導致咸水體下移，降水入滲填充了淺部空間，造成咸水淡化。

1.2.2 淺層咸水入侵深層淡水

衡水市深層水降落漏斗形成于20世紀70年代，這一時期工、農業發展迅速，用水量需求持續增加，為了滿足發展的需求，衡水市大量開采深層地下水，長期的持續超采導致地下水水位持續下降，深、淺層水之間形成巨大水頭差，為上部淺層咸水越流補給深層淡水提供了水動力條件，最終造成咸水入侵問題[5]。

2 水質及咸淡水界面變化研究

分析前人所做研究并綜合咸水成因以及咸水入侵機理，發現咸水入侵問題的研究主要集中在2個方面：水質變化研究和水質預測及咸淡水界面變化研究。

2.1 水質變化研究

地下水水質變化的研究主要是研究地下水的化學成分在各種因素影響下時間和空間上的變化。對于地下水水質變化研究的方法較多，例如內梅羅指數法、模糊數學、物元可拓法、集對分析法、綜合指數法、單因子評價法、灰色聚類法等。

國外對于地下水水質變化的研究較多。1993年M Loizidou等[6]對雅典地區垃圾填埋場周圍地下水進行取樣監測發現滲濾液對地下水水質影響較大。1999年X.Song等[7]運用化學及同位素技術對Sri Lanka西部地區的水質演化進行研究，發現地下水中主要離子的濃度和電導率的變化與季節的相關性較小。1991年荷蘭建立了國家地下水水質數據庫（INGRES），1991年美國地質調查局開始執行國家水質評價計劃[8]。2011年Polat等[9]對季節性地下水水質變化進行評價。2016年Takeshi Saito等[10]研究了地下溫度變化對地下水質量變化的影響，測試了幾種組分的隨溫度的變化趨勢。

20世紀90年代以來國內對于水質變化也做了大量的研究。沈照理[11]在“水文地球化學基礎”中對地下水化學成分的特征進行了詳細講述。張宗祜等[12]描述了華北平原地下水環境的演化與發展。邱漢學等[13]敘述了海水入侵區地下水化學成分的形成作用。高太忠等[14]通過取樣對比分析發現武邑縣圈頭地區淺層淡水有被咸化的現象，并對20世紀80—90年代的水位埋深進行分析，得出河北平原地下水位持續下降、水資源逐年衰竭的結論，提出實現水資源和社會經濟可持續發展的對策。趙勇勝等[15]將水質預警模型與GIS技術相結合，建立地下水水質預警系統，并且應用于吉林省西部平原地區地下水水質預警。謝菲等[16]分析了涇惠渠灌區地下水基本特征、地下水化學成分的形成原因以及不同灌溉水源區地下水的暫時性變化規律。欒鳳嬌等[17]采用單因子評價法對葉爾羌河流域平原區農村地下水水質進行評價，評價結果符合研究區的實際情況。

對于水質變化的研究國內外學者都非常重視，在研究中運用了多種理論方法，為水質變化研究做出了大量的貢獻。

2.2 水質預測及咸淡水界面變化研究

20世紀初國外就已經利用數學模型研究地下水問題，對于水質研究的數學模型應用開始于20世紀60年代。1976年前蘇聯科學家作出了對孔隙介質中水動力彌散研究的詳細綜述[18]。1972年J.Fried[19]通過對經典模型和水動力彌散方程的研究，提出了新的水動力彌散方程。Birkholzer等[20]利用美國地質調查局NWIS數據庫建立了含水層地球化學模型。Ben V.Furlong等[21]利用區域地下水流動模型預測區域井水水質分布。2016年Sakizadeh等[22]采用神經網絡、神經網絡集成神經網絡和貝葉斯正則化神經網絡對地下水水質指數進行了預測，并對3種預測方法的性能進行了對比。近年來隨著溶質運移模型的廣泛研究，獲得的成果越來越多。

國內對于地下水水質模型的研究開始于20世紀80年代。吳吉春等[23]首次建立了越流含水層系統的地下水污染數學模型，并且導出了包括越流項和井流項的溶質運移方程。1982年鄧聚龍教授[24]提出了灰色系統理論，實現了對地下水水質的科學預測，其中的GM（1，1）模型是常用的預測模型。孫鳴等[25]詳細敘述了GM （1，1）模型的建模、精度檢驗和模型修正過程，并以大武水源地為例預測了水質變化趨勢。呂海濤等[26]指出了應用灰色模型時注意的問題。近年來人工神經網絡的預測應用發展迅速，劉國東等[27]以成都市金堂縣東風水庫水質資料為例，建立了地面水水質綜合評價的BP人工神經和Hopfield網絡模型，結果發現Hopfield網絡更優于BP人工神經網絡。王秉忱等[28]對地下水污染與水質模擬實驗的基本理論作出敘述，以濟寧市區地下水污染發展為例建立地下水動力彌散方程并展示了詳細的求解過程，指明求解過程只適用于以彌散為主的問題。郭永海等[29]對河北平原地下水環境進行模擬，結果表明咸水入侵深層淡水并不是簡單地機械混合，其中發生了多種物質的溶解、沉淀、陽離子交換等水文地球化學作用。吳吉春等[30]通過野外實地調查，建立山西柳林泉裂隙發育區三維可混溶溶質運移模型，將模擬結果與實際情況對比，證實三維模型的建立是合理的。苗添升等[31]為了減小參數取值的偶然性，提高模擬結果的合理性，使模擬結果更具參考價值，使用蒙特卡洛方法對地下水數值模擬模型進行不確定性分析。趙春蘭等[32]利用Modflow和MT3D軟件模擬垃圾滲濾液對地下水水質的影響。杜超等[33]采用VB語言編程技術建立了下遼河平原地下水水質預測模型。王現國等[34]利用PHREEQC軟件對三門峽盆地地下水化學成分進行了水文地球化學模擬，對地下水化學環境演化過程定量評價。張將偉等[35]以河谷地區為原型建立地表水地下水的水流水質聯合模擬模型，應用HydroGeoSphere軟件對構建的模型進行求解，對未來水流和總氮的變化特征進行預測。武雅潔等[36]利用OpenGeoSys軟件建立二維飽和多孔介質變密度地下水流的溶質運移模型，對地下截滲墻影響下的咸水入侵問題進行了數值模擬。隨著近年來計算機模擬軟件的推廣，地下水水質模擬預測的技術日趨成熟，準確度也越來越高。

國內對于咸水入侵問題的研究始于20世紀80年代后期，通過觀察咸淡水界面位置的改變可以看出咸水入侵的趨勢，是對咸水入侵最直觀的體現。薛禹群，吳吉春等對這一方面研究較多。薛禹群等[37]在研究咸水入侵問題時，認為對水頭和鹽分在垂向上的動態變化分析是十分重要的，最好的辦法是建立三維監測網，進行長期觀測，通過設置長期監測點可以更好地觀察咸水體的運移趨勢，為咸水入侵的動態監測提供基礎資料。明木河等[38，39]對衡水市咸淡水界面的變化做了大量研究。出于對咸淡水界面運移預測中突變界面模型和動力彌散界面模型等確定性模型方法局限性的考慮，潘世兵等[40]運用灰色建模原理提出一種非確定性模型方法-灰色突變模型方法，并應用于廣饒縣井灌區咸水入侵模擬預測，模擬結果表明灰色突變模型方法具有簡潔有效的特點，但是這種模型方法的應用必須要具備系列長期的監測數據，并且只適用于淺層咸水入侵趨勢的研究，其中山東省廣饒縣水利局針對咸淡水混合區布置了咸水入侵監測網，為咸淡水界面運移規律分析和預測提供了基礎數據和資料。王永昌等[41]認為，氣候和過量開采地下水是導致滄州地區咸水入侵的主要原因，通過搜集河間市300諸多深機井資料，將各測井的咸水底界埋深值及成井時間與背景值（1967年）對照比較，證實了20世紀60—90年代近30a的時間內河間市咸水底界面全市平均下移10.69m。牟純儒等[42]對河北省典型區咸淡水界面下移現狀和機制進行了分析。張亞哲等[43]將礦化度的大小作為咸淡水分界的指標，詳細論述河北平原地下咸水的分布及儲量，并且利用多年的動態監測結果證實咸淡水界面呈下移趨勢，提出咸水的開發利用應遵循“以用為主，以用促改”的策略。張素娥等[44]利用典型井孔礦化度的動態監測值證實了衡水市在淺層咸水分布區均有不同程度的咸水入侵深層淡水現象，并指出超采深層地下水是導致咸水入侵深層淡水的主要原因，通過統計物探測井資料計算出了咸淡水界面下移的多年平均值。費宇紅、張兆吉、宋海波等[45]對華北平原咸水入侵問題研究較多，通過對地下水水質動態監測、機電測井、成井物探等資料的分析，證實了咸淡水界面下移的事實，并提出地下水開采是造成咸淡水界面下移的主要原因。孫曉敏等[46]利用取自研究區的原狀土進行室內化學滲透模擬試驗，實驗結果表明化學滲透效應在一定程度上會影響咸水下移的速度。周曉妮等[47]對咸淡水界面的空間分布特征及其演變的過程進行了野外觀測。

大量學者針對咸水入侵問題從入侵機理、入侵速度、咸淡水界面運移預測、水質變化預測、微咸水、咸水的開發利用等諸多方面進行了研究，研究成果對于防治咸水入侵問題具有重要意義。

3 咸水入侵問題的應對措施

3.1 減少地下水的開采

根據咸水入侵的機理，防止咸水入侵問題最直接的方式是減少對地下水地開采，通過合理規劃地下水的開采量，可以有效地控制地下水位，減小咸淡水之間的水頭差。以研究所得的結論作為科學依據，為當地水資源的開發利用提出合理化建議。

3.2 開發利用咸水、微咸水

前人對于華北平原咸水灌溉利用問題已經做了較多研究，研究結果表明：對于礦化度大于3.0g/L的地區，實行咸淡水混合灌溉，通過咸淡水的混合降低灌溉水的含鹽量。在礦化度小于3.0g/L的地區，可以直接進行灌溉[48，49]。

3.3 建立健全咸水入侵監測網和檢測系統

在全市范圍內建立咸水入侵三維監測網，對地下水水質進行動態監測，利用數值模擬技術建立咸水入侵模型，預測咸水的發展趨勢。

參考文獻

[1] 潘俊， 劉洋， 常玉輝， 等. 下遼河平原南部第三系含水層咸水入侵的調控[J]. 沈陽建筑大學學報（自然科學版）， 2013， 29（04）：710-714.

[2]河北省地礦局第三水文工程地質大隊地質環境監測站. 河北省衡水市地質環境監測報告（2006—2010年）[R]. 2011.

[3]周曉妮， 劉少玉， 王哲， 等. 華北平原典型區淺層地下水化學特征及可利用性分析——以衡水為例[J]. 水科學與工程技術， 2008（02）：56-59.

[4]張宗祜， 沈照理， 薛禹群. 華北平原地下水環境演化[M]. 地質出版社， 2000.

[5]白春艷， 周金龍. 中鹽度地下水灌溉試驗研究現狀及展望[J]. 地下水， 2013， 35（01）：69-71.

[6]Loizidou M， Kapetanios E G. Effect of leachate from landfills on underground water quality[J]. Science of the Total Environment， 1993， 128（01）：69.

[7]Conceptual model of the evolution of groundwater quality at the wet zone in Sri Lanka[J]. Environmental Geology， 1999， 39（02）：149-164.

[8]周仰效，李文鵬. 地下水水質監測與評價[J]. 水文地質工程地質， 2008， 35（01）：1-11.

[9]ElI A， Polat R. Assessment of the statistical significance of seasonal groundwater quality change in a karstic aquifer system near Izmir-Turkey[J]. Environmental Monitoring & Assessment， 2011， 172（1-4）：445-462.

[10]Masoud A A， Koike K， Mashaly H A， et al. Spatio-temporal trends and change factors of groundwater quality in an arid area with peat rich aquifers： Emergence of water environmental problems in Tanta District， Egypt[J]. Journal of Arid Environments， 2016（124）：360-376.

[11]沈照理. 水文地球化學基礎[M]. 地質出版社， 1993.

[12]張宗祜， 施德鴻， 沈照理， 等. 人類活動影響下華北平原地下水環境的演化與發展[J]. 地球學報， 1997， 18（04）：337-344.

[13]邱漢學， 劉貫群. 海水入侵區地下水化學成分的形成作用[J]. 工程勘察， 1999（01）：43-46.

[14]高太忠， 羅人明， 黃群賢， 等. 河北平原地下水開發中的問題及可持續利用[J]. 地質災害與環境保護， 2002， 13（03）：6-9.

[15]洪梅， 趙勇勝， 張博. 地下水水質預警信息系統研究[J]. 吉林大學學報（地）， 2002， 32（04）：364-368.

[16]謝菲， 蔡煥杰， 趙春曉， 等. 涇惠渠灌區地下水化學特征及其對不同水源灌溉的響應[J]. 灌溉排水學報， 2017（04）：79-84.

[17]欒風嬌， 周金龍， 曾妍妍， 等. 葉爾羌河流域平原區農村地下水水質安全評價[J]. 人民黃河， 2017（06）：68-71.

[18]王秉忱， 楊天行， 王寶金. 地下水污染地下水水質模擬方法[M]. 北京師范學院出版社， 1985.

[19]Fried J J. Groundwater Pollution： Theory， Methodology， Modelling and Practical Rules[M]. Amsterdam： Elsevier Scientific Publishing Company， 1975.

[20]Birkholzer J， Apps J A， Zheng L， et al. Prediction of Groundwater Quality Changes in Response to CO2 Leakage from Deep Geological Storage[C]// Agu Fall Meeting. AGU Fall Meeting Abstracts， 2008：1887-1894.

[21]Furlong B V， Riley M S， Herbert A W， et al. Using regional groundwater flow models for prediction of regional wellwater quality distributions[J]. Journal of Hydrology （Amsterdam）， 2011， 398（1-2）：1-16.

[22]Sakizadeh， Mohamad. Artificial intelligence for the prediction of water quality index in groundwater systems[J]. Modeling Earth Systems and Environment， 2016， 2（01）：8.

[23]吳吉春， 薛禹群， 謝春紅. 越流含水層系統中的溶質運移方程[J]. 水文地質工程地質， 1998（01）：30-31.

[24]鄧聚龍.灰色控制理論[M].華中理工大學出版社， 1985.

[25]孫鳴， 徐品. 灰色模型在地下水水質預測中的應用[J]. 河北地質大學學報， 1994（02）：192-197.

[26]呂海濤， 趙林明. 應用灰色模型（GM）時注意的一些問題[J]. 華北水利水電大學學報（自然科學版）， 2000， 21（01）：76-78.

[27]劉國東， 黃川友， 丁晶. 水質綜合評價的人工神經網絡模型[J]. 中國環境科學， 1998， 18（06）：514-517.

[28]王秉忱， 楊天行， 侯印偉， 等. 地下水質模擬試驗研究的基本理論與計算方法——以濟寧市區地下水污染發展趨勢預測的地下水質模擬為例[J]. 吉林大學學報（地球科學版）， 1981（03）：106-124.

[29]郭永海， 沈照理， 鐘佐燊. 河北平原地下水化學環境演化的地球化學模擬[J]. 中國科學（D輯：地球科學）， 1997（04）：360-365.

[30]吳吉春， 張政治. 山西柳林泉裂隙發育區溶質運移三維數值模擬[J]. 南京大學學報（自然科學）， 2000， 36（06）：728-734.

[31]苗添升， 盧文喜， 歐陽琦， 等. 地下水數值模擬的不確定性分析在水質預測中的應用[J]. 水電能源科學， 2016（08）：20-23.

[32]趙春蘭， 凌成鵬， 吳勇， 等. 垃圾滲濾液對地下水水質影響的數值模擬預測——以冕寧縣漫水灣生活垃圾填埋場為例[J]. 環境工程， 2017（02）：163-167.

[33]杜超， 肖長來， 呂軍， 等. Development and Application of Real Time Groundwater Quality Prediction System： An Example in Lower Liaohe River Plain[J]. Journal of Jilin University， 2014， 44（05）：1625-1632.

[34]王現國， 楊國華， 谷芳瑩， 等. 三門峽盆地地下水化學成分演化機理與模擬[J]. 人民黃河， 2018， 40（06）：82-86.

[35]張將偉， 盧文喜， 陳末， 等. 基于HydroGeoSphere的河谷地區地表水地下水水流水質聯合模擬[J]. 中國農村水利水電， 2018（01）：29-32.

[36]武雅潔， 馮峰， 雷鑫. 地下截滲墻影響下的咸水入侵規律研究[J]. 中國海洋大學學報（自然科學版）， 2018（08）：131-138.

[37]薛禹群， 吳吉春， 謝春紅， 等. 萊州灣沿岸海水入侵與咸水入侵研究[J]. 科學通報， 1997（22）：2360-2368.

[38]明木和. 衡水深層淡水咸化問題的探討[J]. 水文地質工程地質， 1986（05）：46-48.

[39]明木和， 沈珍瑤. 咸水底界變化規律初探[J]. 物探與化探， 1992（05）：68-73.

[40]潘世兵， 張福存， 安永會， 等. 咸水入侵趨勢預測的非確定模型研究[J]. 勘察科學技術， 1999（06）：3-6.

[41]王永昌. 滄州區域咸水入侵初探[J]. 地下水， 1995（02）：78-80.

[42]牟純儒， 張建平. 河北省典型區咸淡水界面下移現狀及其入侵機制分析[J]. 水科學與工程技術， 2002（01）：37-39.

[43]張亞哲， 申建梅， 王瑩， 等. 河北平原地下（微）咸水的分布特征及開發利用[J]. 農業資源與環境學報， 2009， 26（06）：29-33.

[44]張素娥， 寇秋煥， 李志軍. 河北省衡水市地質災害分布特征及其成因分析[J]. 地質災害與環境保護， 2006， 17（04）：8-12.

[45]費宇紅， 張兆吉， 宋海波， 等. 華北平原地下咸水垂向變化及機理探討[J]. 水資源保護， 2009， 25（06）：21-23.

[46]孫曉敏， 吳劍鋒. 華北平原咸-淡含水層間粘性土層的化學滲透效應研究[J]. 南京大學學報（自然科學）， 2016， 52（03）：421-428.

[47]周曉妮， 石建省， 馬榮， 等. 咸淡水界面空間分布特征演變的野外觀測[J]. 中國科學院大學學報， 2017， 34（05）：618-624.

[48]孫曉明， 徐建國， 楊齊青， 等. 環渤海地區海（咸）水入侵特征與防治對策[J]. 地質調查與研究， 2006， 29（03）：203-211.

[49]劉小京. 環渤海缺水區鹽堿地改良利用技術研究[J]. 中國生態農業學報， 2018， 26（10）：102-108.