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

鄧興 陳康 劉雪友

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

關(guān)鍵詞：咸水入侵;地下水水質(zhì);咸淡水界面

引言

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

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

1 咸水成因及入侵機理

1.1 咸水成因

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

1.2 咸水入侵機理

1.2.1 淺層咸水入侵淺層淡水

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

1.2.2 淺層咸水入侵深層淡水

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

2 水質(zhì)及咸淡水界面變化研究

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

2.1 水質(zhì)變化研究

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

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

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

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

2.2 水質(zhì)預(yù)測及咸淡水界面變化研究

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

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

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

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

3 咸水入侵問題的應(yīng)對措施

3.1 減少地下水的開采

根據(jù)咸水入侵的機理，防止咸水入侵問題最直接的方式是減少對地下水地開采，通過合理規(guī)劃地下水的開采量，可以有效地控制地下水位，減小咸淡水之間的水頭差。以研究所得的結(jié)論作為科學(xué)依據(jù)，為當(dāng)?shù)厮Y源的開發(fā)利用提出合理化建議。

3.2 開發(fā)利用咸水、微咸水

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

3.3 建立健全咸水入侵監(jiān)測網(wǎng)和檢測系統(tǒng)

在全市范圍內(nèi)建立咸水入侵三維監(jiān)測網(wǎng)，對地下水水質(zhì)進行動態(tài)監(jiān)測，利用數(shù)值模擬技術(shù)建立咸水入侵模型，預(yù)測咸水的發(fā)展趨勢。

參考文獻

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

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

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

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

[5]白春艷， 周金龍. 中鹽度地下水灌溉試驗研究現(xiàn)狀及展望[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]周仰效，李文鵬. 地下水水質(zhì)監(jiān)測與評價[J]. 水文地質(zhì)工程地質(zhì)， 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]沈照理. 水文地球化學(xué)基礎(chǔ)[M]. 地質(zhì)出版社， 1993.

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

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

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

[15]洪梅， 趙勇勝， 張博. 地下水水質(zhì)預(yù)警信息系統(tǒng)研究[J]. 吉林大學(xué)學(xué)報（地）， 2002， 32（04）：364-368.

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

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

[18]王秉忱， 楊天行， 王寶金. 地下水污染地下水水質(zhì)模擬方法[M]. 北京師范學(xué)院出版社， 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]吳吉春， 薛禹群， 謝春紅. 越流含水層系統(tǒng)中的溶質(zhì)運移方程[J]. 水文地質(zhì)工程地質(zhì)， 1998（01）：30-31.

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

[25]孫鳴， 徐品. 灰色模型在地下水水質(zhì)預(yù)測中的應(yīng)用[J]. 河北地質(zhì)大學(xué)學(xué)報， 1994（02）：192-197.

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

[27]劉國東， 黃川友， 丁晶. 水質(zhì)綜合評價的人工神經(jīng)網(wǎng)絡(luò)模型[J]. 中國環(huán)境科學(xué)， 1998， 18（06）：514-517.

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

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

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

[31]苗添升， 盧文喜， 歐陽琦， 等. 地下水?dāng)?shù)值模擬的不確定性分析在水質(zhì)預(yù)測中的應(yīng)用[J]. 水電能源科學(xué)， 2016（08）：20-23.

[32]趙春蘭， 凌成鵬， 吳勇， 等. 垃圾滲濾液對地下水水質(zhì)影響的數(shù)值模擬預(yù)測——以冕寧縣漫水灣生活垃圾填埋場為例[J]. 環(huán)境工程， 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]王現(xiàn)國， 楊國華， 谷芳瑩， 等. 三門峽盆地地下水化學(xué)成分演化機理與模擬[J]. 人民黃河， 2018， 40（06）：82-86.

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

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

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

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

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

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

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

[42]牟純?nèi)澹?張建平. 河北省典型區(qū)咸淡水界面下移現(xiàn)狀及其入侵機制分析[J]. 水科學(xué)與工程技術(shù)， 2002（01）：37-39.

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

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

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

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

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

[48]孫曉明， 徐建國， 楊齊青， 等. 環(huán)渤海地區(qū)海（咸）水入侵特征與防治對策[J]. 地質(zhì)調(diào)查與研究， 2006， 29（03）：203-211.

[49]劉小京. 環(huán)渤海缺水區(qū)鹽堿地改良利用技術(shù)研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報， 2018， 26（10）：102-108.