基于GM的滄縣地區(qū)淺層地下水控制性管理水位研究

李迎龍+李海明+謝新民+董旺泉+葛亞超+王燁

摘要：由于地下水的超采，已在滄縣地區(qū)引發(fā)一系列環(huán)境地質(zhì)問題。通過確定地下水控制性管理水位，來控制地下水的開采量，實(shí)現(xiàn)滄縣地區(qū)淺層地下水資源的可持續(xù)利用。根據(jù)滄縣地區(qū)淺層水位觀測井的觀測資料，分析了滄縣地區(qū)淺層地下水位動(dòng)態(tài)特征，發(fā)現(xiàn)地下水位呈周期性變化。在充分了解滄縣地區(qū)水文地質(zhì)條件的基礎(chǔ)上，利用GMS軟件建立滄縣地區(qū)地下水水流模型。在現(xiàn)有的開采基礎(chǔ)上，預(yù)測不同頻率降水量情景下滄縣地區(qū)淺層地下水控制性管理水位，并確定其閾值。四種情景下滄縣地區(qū)淺層地下水控制性管理水位閾值分別為：-640～700 m、-430～730 m、-470～710 m、-450～700 m。最終確定滄縣地區(qū)淺層地下水控制性管理水位閾值為-640～730 m。

關(guān)鍵詞：水位動(dòng)態(tài)；數(shù)值模擬；控制性管理水位；閾值；滄縣地區(qū)

中圖分類號：P641文獻(xiàn)標(biāo)識碼：A文章編號：

16721683（2017）06010807

Abstract：Due to the overexploitation of groundwater，a series of environmental and geological problems have been triggered in Cangxian areaThis paper was intended to determine the controlled management level of groundwater as a way to control the exploitation quantity of groundwater and to realize sustainable utilization of the shallow groundwater resources in Cangxian areaBased on the observation data from the shallow water level observation wells in Cangxian area，we analyzed the dynamic characteristics of shallow groundwater level in Cangxian area，and found a periodical change of the groundwater levelBased on a substantial understanding of the hydrogeological conditions in Cangxian area，we developed a groundwater flow model of Cangxian area with the GMS softwareOn the basis of the current exploitation situation，we predicted the controlled management level of shallow groundwater in Cangxian area under 4 different precipitation frequency scenarios，and determined the threshold valueThe threshold values under the 4 scenarios were -640～700 m，-430～730 m，-470～710 m，-450-700 m respectivelyThe controlled management water level threshold of the shallow groundwater in Cangxian area was finalized to be -640～730 m

Key words：dynamics of water level；numerical simulation；controlled management water level；threshold；Cangxian area

地下水是中國北方地區(qū)的重要飲用水源。隨著經(jīng)濟(jì)的不斷發(fā)展，對于地下水的開采也逐年加劇。地下水的超采不僅導(dǎo)致了地下水位的下降，而且引發(fā)了一系列的環(huán)境地質(zhì)問題[1]。針對地下水嚴(yán)重超采的現(xiàn)狀，國家實(shí)行最嚴(yán)格的水資源管理制度[2]，明確提出要“嚴(yán)格地下水管理和保護(hù)。加強(qiáng)地下水動(dòng)態(tài)監(jiān)測，實(shí)行地下水取水總量控制和水位控制”。2015年開始實(shí)施的《河北省地下水管理?xiàng)l例》在全國地方性法規(guī)中首次確立了地下水水量水位“雙控”制度。另一方面，各學(xué)者也積極探討水量和水位控制管理的實(shí)施途徑與方法[34]，特別是閆學(xué)軍等[5]提出并開發(fā)了基于GMS的地下水水位水量 “二元”指標(biāo)管理新模式。第15卷 總第93期·南水北調(diào)與水利科技·2017年12月

李迎龍等·基于GMS的滄縣地區(qū)淺層地下水控制性管理水位研究

地下水位的變化是地下水資源量變化的直觀反映，將其與地下水取水總量同時(shí)作為地下水的控制指標(biāo)相對于僅把取水總量作為地下水的管理控制指標(biāo)更加合理。近幾年國外多用建模的思路來對地下水資源管理進(jìn)行研究[6]。Yeh W W G[7]綜述了地下水建模與管理的優(yōu)化方法，Theodossiou N P等[8]研究了非線性模擬和優(yōu)化模型在地下水含水層管理中的應(yīng)用。Gabor Bekesi等[9]根據(jù)地下水位對開采的響應(yīng)，限制性的分配利用地下水資源。Mohamed EI Alfy [10]利用地下水?dāng)?shù)值模擬探究地下水開采的最優(yōu)方案。Gaur S等[11]研究了人工神經(jīng)網(wǎng)絡(luò)和粒子群算法在地下水管理中的應(yīng)用。國內(nèi)對于地下水資源管理的研究集中于地下水控制性管理水位及管理模型的實(shí)用性。謝新民等[12]對地下水控制性關(guān)鍵水位進(jìn)行了研究，提出不同類型地下水紅藍(lán)線水位的分析成果。方樟等[13]以河南省安陽市為例，利用地下水流數(shù)值模擬的方法研究了平原區(qū)地下水控制性管理水位。吳鳴等[14]、李平等[15]、周念清等[16]分別探索了遺傳算法在地下水管理中的應(yīng)用。endprint

本文基于GMS軟件建立滄縣地區(qū)地下水水流模型，研究不同水平年情景下，滄縣地區(qū)淺層地下水控制性管理水位及其閾值，以便為滄縣地區(qū)淺層地下水的合理開采提供監(jiān)管依據(jù)。

1滄縣地區(qū)水文地質(zhì)條件

11研究區(qū)概況

滄縣位于華北平原東部，屬海河流域南系下游[17]，全境介于北緯38°05′至38°33′，東經(jīng)116°27′至117°09′之間。滄縣地勢低平，由西南向東北緩緩傾斜。該地區(qū)屬溫帶半濕潤大陸性季風(fēng)氣候，滄縣多年平均（1981-2005）蒸發(fā)量1 02890 mm，多年平均降雨量55330 mm，降水量年際變化較大，年內(nèi)分配不均勻，水資源先天不足[1819]。滄縣地下水主要賦存于第四系松散巖層的砂土和亞沙土空隙中，其水文地質(zhì)剖面圖見圖1，多年平均（1980-2005）地下水資源量為192億m3，地下水資源嚴(yán)重短缺。滄州市地勢自西南向東北略傾斜，全區(qū)屬?zèng)_積、湖積、海積地形地貌，該區(qū)基底構(gòu)造地處華北斷陷，地質(zhì)構(gòu)造較復(fù)雜。由于連續(xù)超采已形成一系列環(huán)境地質(zhì)問題，如形成地下水降落漏斗、地面沉降等[2021]。

12研究區(qū)范圍及行政分區(qū)

整個(gè)研究區(qū)（即滄縣地區(qū)）包括滄縣及滄州市區(qū)。為了方便研究，且根據(jù)河北省滄州地區(qū)地下水管理的需求，在縣級行政分區(qū)與水文地質(zhì)分區(qū)的基礎(chǔ)上，劃分縣級行政區(qū)地下水管理亞區(qū)或單元，并以此作為縣級行政區(qū)地下水管理的最小基本單元。共分為4個(gè)區(qū)：運(yùn)西Ⅰ區(qū)、運(yùn)西Ⅱ區(qū)、運(yùn)東Ⅲ區(qū)、運(yùn)東Ⅳ區(qū)，其分區(qū)見圖2。

13含水層空間分布

研究區(qū)地下水為第四系孔隙水，分為四個(gè)含水組。第一含水組為淺層地下水，底板埋深約為4000 m。運(yùn)西地區(qū)大部分為淡水，運(yùn)東地區(qū)大部分為咸水；運(yùn)西全區(qū)及運(yùn)東Ⅲ區(qū)含水層巖性主要以粉砂為主，運(yùn)東Ⅳ區(qū)主要以黏砂為主。第二含水組為弱承壓水，運(yùn)西上部為咸水，下部為淡水，到運(yùn)東逐漸變成全咸水；含水層巖性以粉細(xì)砂為主，見有細(xì)砂。第三含水組為深層承壓水，底板埋深22000～42000 m；含水層巖性以細(xì)砂為主，見有中細(xì)砂和粉砂。第四含水組系湖積形成，含水層巖性主要以細(xì)砂為主，見有粉細(xì)砂。

14淺層地下水補(bǔ)、徑、排條件

滄縣地區(qū)淺層地下水以孔隙水為主。地下水的補(bǔ)給以降水入滲和地表入滲為主，其次為地下水回灌及側(cè)向徑流補(bǔ)給；主要排泄方式為潛水蒸發(fā)和人工開采，天然狀態(tài)下地下水的流向與地形傾斜相一致，即由西南流向東北，由于開采影響，地下水流向開采區(qū)。

15地下水位動(dòng)態(tài)特征

151地下水位年際動(dòng)態(tài)分析

滄縣地區(qū)的淺層觀測井共有11眼，其分布見圖3。滄縣地區(qū)地下水位年際變化有一定的周期性，可以分為兩個(gè)階段。一是前一年的6月到9月，地下水位有一個(gè)短暫的上升期；二是從9月到下一年的6月，地下水位處于下降期。2012年淺層地下水年平均水位為421 m，比2011年上升136 m；2012年末平均水位467 m，比2011年同期上升了056 m。2013年淺層地下水年平均水位為478 m，比2012年上升了057 m；2013年末平均水位483 m，比2012年同期上升了016 m，淺層地下水位動(dòng)態(tài)變化曲線見圖4。

152地下水位年內(nèi)動(dòng)態(tài)分析

滄縣地區(qū)地下水位年內(nèi)變化均有一定的周期性，可以分為三個(gè)階段。一是水位下降期，從1月到6月水位逐漸下降，6月水位最低。二是水位上升期，從6月到9月水位逐漸升高，9月水位最高。三是水位下降期，從9月到12月水位再次下降。水位上升的幅度大于水位下降的幅度，水位逐步恢復(fù)。

153地下水位變化與降水量關(guān)系

根據(jù)資料的完整性，選取東關(guān)、張韓店和捷地三個(gè)觀測孔來研究地下水位變化與降水量的關(guān)系，其水位及降水資料從2011年6月-2013年12月。地下水位的變化與降水有關(guān)，降水量越大，水位變化幅度也越大，但具有一定的滯后性，滯后期約為1個(gè)月。其關(guān)系見圖5。

2地下水?dāng)?shù)值模型的建立

根據(jù)滄縣地區(qū)水文地質(zhì)條件，結(jié)合地下水長期觀測資料及地下水實(shí)際開采情況等，利用三維地下水流數(shù)值模擬軟件GMS，將滄縣地區(qū)潛水層概化為一個(gè)具有統(tǒng)一水力聯(lián)系的地下水系統(tǒng)，將該系統(tǒng)概化為參數(shù)隨空間變化、非均質(zhì)的非穩(wěn)定流地下水系統(tǒng)。

21水文地質(zhì)概念模型

211模型范圍及邊界條件

整個(gè)研究區(qū)包括滄縣全境及整個(gè)滄州市區(qū)，邊界條件具體如下：（1） 研究區(qū)AG、BC、EF段有河流經(jīng)過，因此將其概化為定水頭邊界；（2）AB、CD、FG段，從2012年地下水位等值線分析，該邊界與地下水位等值線基本垂直，可視為隔水邊界；（3）DE段與研究區(qū)外有水量交換，因此將其概化為流量邊界。模型范圍及地下水位等值線見圖6。

212水文地質(zhì)參數(shù)分區(qū)及初值的確定

本次用于地下水流模型的水文地質(zhì)參數(shù)主要分為兩類：（1）用于計(jì)算地下水的源匯項(xiàng)參數(shù)，如降水入滲補(bǔ)給系數(shù)、灌溉田間入滲補(bǔ)給系數(shù)、潛水蒸發(fā)系數(shù)等；（2）含水層的水文地質(zhì)參數(shù)，主要包括潛水含水層的滲透系數(shù)（k）、給水度（μ），滲透系數(shù)的參數(shù)分區(qū)如圖7、圖8所示。根據(jù)滄縣地下水資源評價(jià)確定水文地質(zhì)參數(shù)初值。

22地下水流數(shù)值模型

221數(shù)學(xué)模型

根據(jù)研究區(qū)水文地質(zhì)條件，將研究區(qū)內(nèi)水文地質(zhì)概念模型概化為水平占主導(dǎo)，垂向運(yùn)動(dòng)較弱、參數(shù)隨空間變化，非均質(zhì)、非穩(wěn)定地下水流系統(tǒng)。

222地下水流數(shù)值模型的建立

（1） 空間剖分。

研究區(qū)面積為1 46950 km2，長為5589 km，寬4480 km。，進(jìn)行網(wǎng)格剖分，共剖分為10 000個(gè)網(wǎng)格（100行×100列），其中活動(dòng)單元格5 017個(gè)，見圖9。

（2） 定解條件。

初始水位：以2012年1月地下水位資料為基礎(chǔ)，采用反距離加權(quán)插值法獲取含水層初始水位。endprint

（3） 源匯項(xiàng)處理。

主要的補(bǔ)給項(xiàng)有：降水入滲補(bǔ)給、側(cè)向補(bǔ)給、灌溉補(bǔ)給等；主要的排泄項(xiàng)有：人工開采、蒸發(fā)及側(cè)向排泄等。采用達(dá)西定律確定邊界側(cè)向補(bǔ)給、排泄量；根據(jù)降水量乘以降水入滲補(bǔ)給系數(shù)，得到降水入滲補(bǔ)給量；灌溉入滲補(bǔ)給量、蒸發(fā)排泄量的確定方法與降水入滲補(bǔ)給量確定方法相同；開采量以井的形式輸入模型。

223模型的識別和驗(yàn)證

本次模擬采用2012年1月1日至2012年12月31日的實(shí)測地下水位數(shù)據(jù)對模型進(jìn)行識別，采用2013年1月1日至2013年12月31日實(shí)測地下水位數(shù)據(jù)對模型進(jìn)行驗(yàn)證。在研究區(qū)內(nèi)選取11個(gè)觀測孔，對模型進(jìn)行識別及驗(yàn)證，發(fā)現(xiàn)計(jì)算值與觀測值在Y=X直線附近，說明計(jì)算值與觀測值相當(dāng)，觀測孔擬合結(jié)果見圖10。

根據(jù)觀測孔計(jì)算值與觀測值的擬合結(jié)果可知，有超過80%的觀測孔誤差范圍在05 m之內(nèi)，說明該模型的精度較高，可以用來預(yù)測。

3各分區(qū)控制性管理水位及閾值

31不同降水典型水平年的確定

利用滄縣地區(qū)1987-2013年共34年資料進(jìn)行頻率分析，采用經(jīng)驗(yàn)頻率線求得不同典型年的降水量頻率。在降水保證率為20%時(shí)視為豐水年，代表年份為2010年；在降水保證率為50%時(shí)視為平水年，代表年份為1998年；在降水保證率為75%時(shí)視為枯水年，代表年份為2007年。

32壓采

為減緩滄縣地區(qū)地下水的過度開采，滄縣采取逐步壓采地下水的措施。要求2015年地下水的開采量控制在077億m3，2020年地下水的開采量控制在062億m3，2030年地下水的開采量減少到042億m3。由于僅有壓采總量的資料，本次預(yù)測根據(jù)2013年各鄉(xiāng)鎮(zhèn)開采量，按比例將壓采后開采量分配到各鄉(xiāng)鎮(zhèn)，最終分配到各開采井。

33不同情景設(shè)計(jì)

利用已建立的地下水流模型，以2014年1月地下水流場作為初始流場，預(yù)測期為2015年1月1日到2030年1月1日。通過對不同情景下的地下水流場的預(yù)測，得到不同情景下代表性監(jiān)測點(diǎn)的地下水位，作為該情景下地下水控制性管理水位。

情景一：選取對供水最不利的降水量連續(xù)特枯水年組33240 mm，把對應(yīng)的2002年降水量逐月資料帶入輸入模型計(jì)算。地下水在2015-2030年按照滄縣壓采方案逐步減少開采量。在各開采量的基礎(chǔ)上預(yù)測2015-2030年地下水位的變化情況。

情景二：選取降水量連續(xù)枯水年組44820 mm，把對應(yīng)的2007年降水量逐月資料帶入模型計(jì)算。地下水在2015-2030年按照滄縣壓采方案逐步減少開采量。在各開采量的基礎(chǔ)上預(yù)測2015-2030年地下水位的變化情況。

情景三：考慮地下水以豐補(bǔ)歉的調(diào)節(jié)功能，降水量取多年平均值53090 mm（平水年），把對應(yīng)的1998年降水量逐月資料帶入模型計(jì)算。地下水在2015-2030年按照滄縣壓采方案逐步減少開采量。在各開采量的基礎(chǔ)上預(yù)測2015-2030年地下水位的變化情況。

情景四：連續(xù)豐水年組63940 mm，把對應(yīng)的2010年降水量逐月資料帶入模型計(jì)算。地下水在2015-2030年按照滄縣壓采方案逐步減少開采量。在各開采量的基礎(chǔ)上預(yù)測2015-2030年地下水位的變化情況。

34各分區(qū)地下水位閾值確定

（1）情景一。

研究區(qū)淺層地下水控制性管理水位整體上從西向東呈現(xiàn)出先降后升再降的趨勢。淺層地下水控制性管理水位的最大值與最小值均出現(xiàn)在運(yùn)東Ⅳ區(qū)，其值分別為700 m和-640 m。在黃遞鋪鎮(zhèn)形成了局部的降落漏斗，見圖11。從地下水控制性管理水位閾值區(qū)間來看，情景一下滄縣地區(qū)淺層地下水位閾值區(qū)間范圍在-640～700 m之間。

（2）情景二。

研究區(qū)淺層地下水控制性管理水位整體上從西向東呈現(xiàn)出先降后升再降的趨勢。淺層地下水控制性管理水位的最大值為730 m，出現(xiàn)在運(yùn)東Ⅳ區(qū)。最小值為-430 m，出現(xiàn)在運(yùn)西Ⅱ區(qū)。從地下水控制性管理水位閾值區(qū)間來看，情景二下滄縣地區(qū)淺層地下水位閾值區(qū)間范圍在-430～730 m之間。

（3）情景三。

研究區(qū)淺層地下水控制性管理水位整體上從西向東呈現(xiàn)出先降后升再降的趨勢。 淺層地下水控制性管理水位的最大值為710 m，出現(xiàn)在運(yùn)東Ⅳ區(qū)。最小值為-470 m，出現(xiàn)在研究區(qū)運(yùn)西Ⅱ區(qū)。從地下水控制性管理水位閾值區(qū)間來看，情景三下滄縣地區(qū)淺層地下水位閾值區(qū)間范圍在-470～710 m之間。

（4）情景四。

研究區(qū)淺層地下水控制性管理水位整體上從西向東呈現(xiàn)出先降后升再降的趨勢。淺層地下水控制性管理水位的最大值為700 m，出現(xiàn)在運(yùn)東Ⅳ區(qū)。最小值為-450 m，出現(xiàn)在研究區(qū)運(yùn)東Ⅲ區(qū)。從地下水控制性管理水位閾值區(qū)間來看，情景四下滄縣地區(qū)淺層地下水位閾值區(qū)間范圍在-450～700 m之間。

4結(jié)果與討論

不同情景可以確定各自的地下水控制性管理水位，將所有情景整合在一起構(gòu)成一個(gè)由最小值到最大值的區(qū)間即為地下水位閾值。滄縣地區(qū)地下水控制性管理水位閾值為-640～730 m，各行政分區(qū)淺層地下水位閾值及水位埋深閾值見表1。

研究區(qū)淺層地下水控制性管理水位整體上自西向東呈現(xiàn)出先降后升再降的趨勢。相對來說運(yùn)西Ⅱ區(qū)與運(yùn)東Ⅳ區(qū)的水位比運(yùn)西Ⅰ區(qū)和運(yùn)東Ⅲ區(qū)低。運(yùn)西Ⅱ區(qū)水位較低，主要原因是淺層地下水的開采主要集中在該區(qū)域。

在壓采總量確定的情況下，最終確定各行政分區(qū)地下水位管理閾值，該閾值可以用來監(jiān)督管理滄縣地區(qū)地下水開采情況。如果地下水位超出該閾值說明地下水超采，需采取相應(yīng)措施嚴(yán)格限制地下水的開采。

5結(jié)論

地下水位的變化具有周期性，其周期與降水的周期一致，且具有滯后性，滯后期月為1個(gè)月。滄縣地區(qū)淺層地下水控制性管理水位閾值為-470～730 m。對于地下水位的監(jiān)測不足，缺乏長期的水位觀測資料，建議增加地下水位監(jiān)測點(diǎn)，對地下水進(jìn)行長期的監(jiān)測，尤其是降落漏斗區(qū)。未來研究的方向集中于地下水控制性管理體系的進(jìn)一步完善。endprint

參考文獻(xiàn)（References）：

[1]王仕琴，宋獻(xiàn)方，王勤學(xué)，等華北平原淺層地下水水位動(dòng)態(tài)變化[J]地理學(xué)報(bào)，2008，63（5）：435445（WANG S Q，SONG X F，WANG Q X，et alDynamic features of shallow groundwater in North China Plain[J]Acta Geographica Sinica，2008，63（5）：435445（in Chinese）） DOI：103321jissn：03755444200805002

[2]左其亭，張保祥，王宗志，等2011年中央一號文件對水科學(xué)研究的啟示與討論[J]南水北調(diào)與水利科技，2011，9（05）：6873（ZUO Q T，ZHANG B X，WANG Z Z，et alRevelation and discussion on water science research from the "2011 No1 Document by the Central Committee of the CPC"[J]SouthtoNorth Water Transfers and Water Science & Technology，2011，9（5）：6873（in Chinese）） DOI：103724SPJ1201201105068

[3]葉勇，謝新民地下水控制性水位管理分區(qū)研究[J]黑龍江水專學(xué)報(bào)，2009，36（1）：116119（YE Y，XIE X MStudy on the management subarea based on groundwater dominate water level[J]Journal of Heilongjiang Hydraulic Engineering College，2009，36（1）：116119（in Chinese）） DOI：103969jissn2095008X200901033

[4]張遠(yuǎn)東，王策地下水取用水總量與水位雙重控制芻議[J]中國水利，2014（9）：79（ZHANG Y D，WANG CDual control on total volume and level of groundwater consumption[J]China Water Resources，2014（9）：79（in Chinese））

[5]閆學(xué)軍，周亞萍，張偉，等地下水開發(fā)利用水位水量“二元”指標(biāo)管理模式研究[J]河北工業(yè)大學(xué)學(xué)報(bào)，2012，41（2）：6568（YAN X J，ZHOU Y P，ZHANG W，et alStudy on the management modes of groundwater levelquality two indicators in exploitation[J]Journal of Hebei University of Technology，2012，41（2）：6568（in Chinese）） DOI：103969jissn10072373201202015

[6]SINGH AGroundwater resources management through the applications of simulation modeling：A review[J]Science of the Total Environment，2014，499：414423DOI：101016jscitotenv201405048

[7]YEH W W GOptimization methods for groundwater modeling and management[J]Hydrogeology Journal，2015，23（6）：10511065

[8]THEODOSSIOU N PApplication of nonlinear simulation and optimisation models in groundwater aquifer management[J]Water Resources Management，2004，18（2）：125141

[9]BDKESI G，McGuire M，MOILER DGroundwater allocation using a groundwater level response management methodGnangara groundwater system，Western Australia[J]Water Resources Management，2009，23（9）：16651683

[10][ZK（#]El Alfy MNumerical groundwater modelling as an effective tool for management of water resources in arid areas[J]Hydrological Sciences Journal，2014，59（6）：12591274

[11]GAUR S，CH S，GRAILLOT D，et alApplication of artificial neural networks and particle swarm optimization for the management of groundwater resources[J]Water Resources Management，2013，27（3）：927941

[12]謝新民，柴福鑫，顏勇，等地下水控制性關(guān)鍵水位研究初探[J]地下水，2007，29（6）：4750，64（XIE X M，CHAI F X，YAN Y，et alPreliminary study on critical depth of ground water table[J]Ground Water，2007，29（6）：4750，64（in Chinese）） DOI：103969jissn10041184200706015endprint

[13]方樟，謝新民，馬喆，等河南省安陽市平原區(qū)地下水控制性管理水位研究[J]水利學(xué)報(bào)，2014，45（10）：12051213（FANG Z，XIE X M，MA Z，et alConfirmation of controlled groundwater management level and demonstration study in Anyang City，Henan Province[J]Journal of Hydraulic Engineering，2014，45（10）：12051213（in Chinese）） DOI：1013243jcnkislxb201410009

[14]吳鳴，吳劍鋒，施小清，等基于諧振子遺傳算法的高效地下水優(yōu)化管理模型[J]吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2015，45（5）：14851492（WU M，WU J F，SHI X Q，et alNew harmonic oscillator genetic algorithm for efficient groundwater optimization and management[J]Journal of Jilin University （Earth Science Edition），2015，45（5）：14851492（in Chinese）） DOI：1013278jcnkijjuese201505203

[15]李平，盧文喜，靳孟貴改進(jìn)的模擬退火遺傳算法在地下水管理中的應(yīng)用[J]水文地質(zhì)工程地質(zhì)，2011，38（3）：913（LI P，LU W X，JIN M GApplication of improved simulated annealing genetic algorithm to groundwater management[J]Hydrogeology and Engineering Geology，2011，38（3）：913（in Chinese）） DOI：103969jissn10003665201103002

[16]周念清，陳劍橋，江思珉基于遺傳算法的模式搜索法求解地下水管理模型[J]勘察科學(xué)技術(shù)，2011（1）：1821，24（ZHOU N Q，CHEN J Qi，JIANG S MPattern search method for solving groundwater management model based genetic algorithm[J]Site Investigation Science and Technology，2011（1）：1821，24（in Chinese）） DOI：103969jissn10013946201101005

[17]呂慶玉滄州市水資源多年特征分析及利用現(xiàn)狀[J]地下水，2014，36（4）：157160（LYU Q YMultiyear characteristics analysis and current utilization situation of water resources in Cangzhou city[J]Ground Water，2014，36（4）：157160（in Chinese）） DOI：103969jissn10041184201404064

[18]卞海彬滄州水資源現(xiàn)狀及預(yù)測分析[J]水利科技與經(jīng)濟(jì)，2011，17（2）：5758（BIAN H BCurrent situation and forecast of water resources in Cangzhou[J]Water Conservancy Science and Technology and Economy，2011，17（2）：5758（in Chinese）） DOI：103969jissn10067175201102022

[19]張可義，張傳河，張雷滄州市水文特性分析[J]水利科技與經(jīng)濟(jì)，2010，16（11）：12111212（ZHANG K Y，ZHANG C H，ZHANG LAnalysis of hydrological characteristics of Cangzhou City[J]Water Conservancy Science and Technology and Economy，2010，16（11）：12111212（in Chinese）） DOI：103969jissn10067175201011005

[20]馮書倉滄州市地下水現(xiàn)狀及地下水監(jiān)測工作研究[J]中國水利，2013（S2）：3940（FENG S CCurrent situation of groundwater and research of its inspection work in Cangzhou City[J]China Water Resources，2013（S2）：3940（in Chinese）） DOI：103969jissn100011232013z1015

[21]龐炳義，王艷麗滄州市水資源現(xiàn)狀與開發(fā)利用對策淺析[J]地下水，2014，36（3）：4546（PANG B Y，WANG Y LAnalysis on current situation and development of water resources in Cangzhou City[J]Ground Water，2014，36（3）：4546（in Chinese））endprint