地表水地下水聯(lián)合耦合模擬進展

張永雷+許玉鳳+潘網(wǎng)生

摘要 地表水與地下水相互作用的耦合模擬在一個地區(qū)水資源的開發(fā)利用與管理中是十分重要的。本文梳理并分析了地表水與地下水耦合模擬中的主要問題、主要方法以及地表水地下水聯(lián)合調(diào)度的主要策略，認為分布式水文模型對較長時空跨度的地表水與地下水水文循環(huán)進行耦合是未來水文模擬的重要方向，以流域為尺度的地表水地下水聯(lián)合模擬調(diào)度將成為提高水資管理水平、保障水資源安全的重要手段，隨著人類活動的不斷加大，流域整體水循環(huán)模擬必須綜合考慮人類活動的影響。

關(guān)鍵詞 地表水；地下水；耦合模擬；聯(lián)合調(diào)度；數(shù)學模型

中圖分類號 P641 文獻標識碼 A 文章編號 1007-5739（2017）24-0153-03

Summary on Joint Coupled Simulation of Surface Water and Groundwater

ZHANG Yong-lei XU Yu-feng PAN Wang-sheng

（School of Tourism and Resources Environment，Qiannan Normal University for Nationalities，Duyun Guizhou 558000）

Abstract The coupling simulation of the interaction between surface water and groundwater is very important in the development，utilization and management of water resources in a region.This paper combed and analyzed the main problems，methods，and strategies in both the coupling simulation and joint scheduling of surface water and groundwater，considering that the use of distributed hydrological model to simulate the hydrological cycle in long span of time and space between the surface water and the groundwater was an important direction in the future；To jointly simulate and deploy the surface water and the groundwater under watershed scale would become a key method in improving the management level and ensuring the safety of water resources；With the increasing human activities，its impact to the overall water cycle simulation of the basin must be taken into account.

Key words surface water；groundwater；coupling simulation；joint scheduling；mathematical model

水是生命之源，水資源短缺已成為全球性的問題。據(jù)統(tǒng)計，在地球上總水儲量中，海水和其他咸水占97.47%，淡水僅占2.53%，其中人類真正可以利用的地表水和淺層地下水僅占全球總水儲量的0.3%[1]。我國人口眾多，水資源短缺問題尤為嚴峻，水資源緊缺的城市已超過550個，2014年人均水資源不足2 000 m3，遼河、黃河、長江、淮河等7大水系水污染嚴重，人均淡水資源量只相當于世界人均淡水資源量的1/4[2-3]。在嚴峻的水資源情勢下，提高水資源利用率和有效治理水污染已成為整個社會的共識，受到人們尤其是學者和政府機構(gòu)的高度重視。在我國，可以利用的淡水資源主要有淺層地下水和以河流水、湖泊水為代表的地表水，而河流、湖泊等水體的污染減少了可利用水資源數(shù)量，加劇了水資源緊缺，已成為制約經(jīng)濟和社會發(fā)展的重要瓶頸并曾引起嚴重的生態(tài)和社會問題[4]；隨著工業(yè)化和城鎮(zhèn)化的深入發(fā)展，地表水污染向地下延伸成為另一個顯著特征，正在逐步危及地下水利用的可持續(xù)性[5]。

要提高水資源利用的可持續(xù)性，把地表水和地下水納入統(tǒng)一的系統(tǒng)進行統(tǒng)籌管理是可行方法之一，也是未來提高水資源可持續(xù)利用的必由之路[6]。因此，了解并追蹤地表水和地下水耦合模擬和調(diào)度的國內(nèi)外研究現(xiàn)狀和研究熱點，對于科學把握我國未來水資源情勢的變化、有效制定相關(guān)政策措施以保證水資源的安全與可持續(xù)利用具有重要意義，同時對于遏制一個地區(qū)地下水位下降和生態(tài)退化也有重要的現(xiàn)實意義[7]。

1 地表水地下水耦合模擬的主要問題

在以流域為尺度的水文循環(huán)中，地表水和地下水密切關(guān)聯(lián)的統(tǒng)一整體，地表水和地下水的循環(huán)轉(zhuǎn)化關(guān)系是探究流域尺度水循環(huán)和計算水量收支必須首先明確的問題[8]。研究地表水和地下水的循環(huán)與轉(zhuǎn)化關(guān)系以及聯(lián)合調(diào)度是對流域水資源進行綜合利用的重要一環(huán)[9]。無論國內(nèi)還是國外，早期的研究往往把地表水與地下水分別作為一個獨立的系統(tǒng)進行研究，割裂了地表水和地下水的聯(lián)系，導致地表水與地下水的交換被忽視，對流域水資源的綜合利用產(chǎn)生了不利的影響[10-11]。從20世紀70年代開始，在制定流域水資源利用規(guī)劃時，研究人員始終將地表水和地下水納入統(tǒng)一的整體進行考慮[12]，近些年對河流水與沿岸淺層地下水的相互影響模擬較多[13]，人類活動影響下流域整體水循環(huán)的模擬研究相對較少，水循環(huán)的演變和水資源的調(diào)度管理缺乏長時空尺度的綜合模擬，這在人水矛盾突出、生態(tài)環(huán)境脆弱的干旱區(qū)內(nèi)陸河流域表現(xiàn)得尤為明顯[14]。endprint

隨著溫室氣體排放所帶來的一系列升溫效應，氣候變暖對一個地區(qū)的水文水資源的影響也愈加引起眾多學者的重視，與之相關(guān)的研究如冰川對升溫的響應、氣候變化對河流—綠洲生態(tài)系統(tǒng)的影響等被置于越來越重要的地位并引起政府的高度重視。在干旱區(qū)地內(nèi)河流域，靠河流維系的綠洲水文過程與生態(tài)過程成為研究的重點和熱點，河流改道或斷流對生態(tài)與地下水位埋深的影響，以及地下水位與沿河生態(tài)對于河道輸水的響應成為研究的核心問題[15-17]，氣候變暖導致極端水文事件的增多與流域內(nèi)地表水地下水管理分割、統(tǒng)籌協(xié)調(diào)不足、開發(fā)利用失控等問題相疊加，使流域水資源情勢比較嚴峻，加之河流水在山前沖積扇和沖積平原上向地下水的轉(zhuǎn)化過程對綠洲生態(tài)環(huán)境影響顯著，一些學者由此提出了建立地表水與地下水統(tǒng)籌調(diào)度模式、建立地表水地下水統(tǒng)一管理機構(gòu)的建議并被政府采納[18-19]。在農(nóng)灌區(qū)，由于化肥的大量使用，引河水灌溉的余水含有大量礦物溶質(zhì)，通過田間入滲反補地下水造成地下水礦化度增加、水質(zhì)惡化，水資源利用的可持續(xù)性受到干擾[20]。

與資源性缺水情況明顯的干旱區(qū)不同，濕潤地區(qū)降水量較豐富，資源性缺水狀況不甚明顯，水資源問題主要是以水體污染、水資源浪費以及管理不當所導致的可利用水水資源減少，即水資源的結(jié)構(gòu)性緊缺[18，21]。在一些平原地區(qū)，人類活動改變了地表水和地下水的穩(wěn)定氫氧同位素的比例，導致河流水水化學類型在人類活動影響下由上游的Na-HCO3型和Ca-Mg-HCO3型向中下游的Ca（Mg）-Cl（SO4）型演化[22]。河水中的污染物質(zhì)通過下滲與側(cè)向入滲進入傍河淺層地下水并在地下水中累積，引起傍河地下水較為嚴重的污染，加大了水污染治理的難度[23]，這種現(xiàn)象在其他流域也具有一定的普遍性。人口增加和農(nóng)業(yè)發(fā)展使用水量成倍增加導致地表水資源出現(xiàn)短缺，人們大量開發(fā)利用地下水，由于缺少長期的和嚴密的規(guī)劃，地下水超采嚴重，一些地區(qū)出現(xiàn)水位下降、地下水漏斗和地面下沉等問題，造成較大的物質(zhì)財富損失，一些學者認為這些地區(qū)必須統(tǒng)籌管理和調(diào)度地表水與地下水才能最大程度發(fā)揮水資源潛力，緩解用水壓力，保障用水安全[24]。

當一個地區(qū)的水資源不足時，跨地區(qū)或跨流域調(diào)水是解決方法之一，但效果有限，已有的事實表明，在地勢比較低平的地區(qū)，外來調(diào)水可以季節(jié)性地改善溝渠水和溝渠附近淺層地下水的水質(zhì)，而對深層地下水和坑塘水水質(zhì)則無明顯的改善作用，而且調(diào)水對溝渠水水質(zhì)的改善僅體現(xiàn)在調(diào)水時節(jié)，對淺層地下水水質(zhì)的改善則明顯滯后[25]。更令人憂心的是，不少地區(qū)與之有關(guān)的科學研究明顯滯后，水資源開發(fā)利用管理術(shù)推廣和升級依舊困難，在一些巖溶地區(qū)，相關(guān)的研究主要集中在水資源安全度評價、水體礦化度、城市地表水地下水交換等方面，對于河流與淺層地下水相互作用、地表水地下水的聯(lián)合調(diào)度等方面的研究較少[26-29]。

2 地表水地下水耦合模擬方法及趨勢

2000年之后，采用數(shù)學模型對地表水地下水的循環(huán)轉(zhuǎn)化關(guān)系進行耦合模擬已成為普遍的做法[30-33]。目前，在研究河流與沿岸淺層地下水相互關(guān)系方面，基于達西定律的研究方法應用最廣，該方法包括壓力水頭觀測、解析解法和數(shù)值模型法等[30]。此外，直接水量測量、間接實驗法、水量平衡法等方法也有比較廣泛的應用。在以流域尺度的水循環(huán)研究中，綜合地表地下水的數(shù)學模型受到青睞，這類模型能分布式地模擬流域地表水文過程及地下水的運動，模擬結(jié)果與實際觀測耦合較好。數(shù)學模型與軟件平臺的結(jié)合可以實現(xiàn)較為復雜的耦合模擬，比如，Larocque等[34]利用HydroGeoSphere模擬了不同情景下含水層的地下水流動，Spanoudaki[8]則建立了耦合三維地表水和三維地下水的IRENE模型，采用半顯示和顯示有限差分求解地表水方程和地下水方程。數(shù)學模型可以與解析解法結(jié)合使用并通過試驗觀測數(shù)據(jù)對模型進行驗證[35-36]。目前，數(shù)學模型在精度、穩(wěn)定或收斂性問題等方面也存在一些問題[37-38]，然而無論是從區(qū)域尺度還是流域尺度研究河流與淺層地下水的循環(huán)與轉(zhuǎn)化，數(shù)學模型均是最主要的選擇對象[30，39]。

國內(nèi)對研究方法的采用多借鑒和改造國外方法，其中示蹤理論法的原理以及傳統(tǒng)示蹤法與溫度示蹤法在地表水地下水相互作用的識別和量化方面應用較廣[40]，根據(jù)驅(qū)動力—壓力—狀態(tài)—影響—響應概念查模型構(gòu)建水資源安全評價框架，運用多級模糊綜合評判法和熵權(quán)法可以定量計算水資源安全度綜合指數(shù)，并初步分析水資源安全度的影響機制[27]；在地表水—地下水的相互轉(zhuǎn)化方面，基于同位素水文學的研究方法在研究河流水與地下水的轉(zhuǎn)化關(guān)系方面有較多的應用[41]，通過水樣品采取分析測定，采用端元法可以定量計算地表水地下水相互轉(zhuǎn)化的比例[22]；基于排隊理論建立灌區(qū)渠系優(yōu)化配水模型可以比較真實地反映用于灌溉的地下水水量隨時間變化的不確定性[42]，改進的分布式水文模型（SWAT）模型與模塊化地下水動力模型（MODFLOW）相結(jié)合可以實現(xiàn)灌區(qū)地表水—地下水分布式模型的耦合，不僅可準確模擬地表水水文過程、地下水動力過程和地下水動態(tài)變化過程，還可用于預測一個地區(qū)地表水和地下水轉(zhuǎn)化規(guī)律的改變和構(gòu)建以流域為尺度的地表水與地下水耦合模型[24，43-44]，而將水量分配效益模型、地下水位控制模型、輸沙效果模型三者相結(jié)合，并通過分層序列法將模型逐步優(yōu)化，則可以得到綜合考慮輸沙效果及地下水水位控制的水量優(yōu)化配置方案[45]。

3 地表水地下水聯(lián)合調(diào)度的主要策略

在水資源調(diào)度方面，進行流域地表水和地下水的聯(lián)合調(diào)度成為國內(nèi)外水資源管理的普遍的重要策略[14，25，46-47]。利用GAMS平臺建立以流域為尺度的地表地下水聯(lián)合調(diào)度優(yōu)化模型，可以評估在人類活動影響下水庫水位的變化對地下水位動態(tài)變化的影響；對于受人類控制較為強烈的河流如干旱區(qū)內(nèi)陸河，生態(tài)輸水的目的是改善河流的水文情勢和下游退化的生態(tài)環(huán)境并在農(nóng)業(yè)用水和生態(tài)用水之間取得協(xié)調(diào)，包含水文—生態(tài)—社會—經(jīng)濟等因素的綜合水資源模型可以用于制定流域水量分配計劃[48]，基于邊際效益相等原則建立的水量優(yōu)化配置模型，能夠較好地反映生態(tài)系統(tǒng)對不同生態(tài)輸水情景的響應，為調(diào)節(jié)生態(tài)輸水時空分配提供參照[49]。通過求解綜合水文—經(jīng)濟—政策的地表水—地下水優(yōu)化調(diào)度模型，可以得到不同生態(tài)流量約束下的社會經(jīng)濟效益的損益表，進而協(xié)調(diào)生態(tài)與經(jīng)濟的平衡發(fā)展。endprint

國內(nèi)學者針對不同地區(qū)的現(xiàn)實情況也提出了較為系統(tǒng)完備的地表水地下水聯(lián)合調(diào)度方案：鄧銘江等[47]提出氣候干旱、降水稀少的地區(qū)可采用建設山間凹陷、山前凹陷和深覆蓋河谷3種地下水庫相結(jié)合的方式，地表水地下水聯(lián)合調(diào)度，以地下水庫調(diào)蓄為主的開發(fā)模式；雍 會等[19]認為完善流域與區(qū)域的地下水資源管理法規(guī)體系完善源流與干流地下水資源管理體制也是統(tǒng)籌管理地表水與地下水的重要舉措；孔曉樂等[25]通過調(diào)研發(fā)現(xiàn)華北平原等半濕潤平原地區(qū)采用調(diào)水、淺層地下水、坑塘水混合灌溉可以滿足農(nóng)業(yè)的用水需求，并使淺層地下水水質(zhì)得到一定必善，深層地下水水位得到一定恢復；陳志兵等[50]和羅 美等[26]提出在西南喀斯特巖溶地貌區(qū)，地表破碎導致水資源蓄積困難，地下水的開發(fā)利用可采取打井、建造小水柜水池以及修筑部分水利工程設的方式，但是地下水的開發(fā)必須把其所在的打井深度、流量、地下地質(zhì)與當?shù)厝丝?、?jīng)濟、生態(tài)環(huán)境等緊密結(jié)合起來才能較好地統(tǒng)籌地表水地下水的聯(lián)合調(diào)度；井柳新等[51]在分析了我國有關(guān)地表水地下水相互轉(zhuǎn)化、污染物在地表水—地下水系統(tǒng)中的運移等方面的研究后認為，要統(tǒng)籌規(guī)劃地表水—地下水的污染防治，建立健全水環(huán)境監(jiān)測體系，深化地表水—地下水水質(zhì)水量相互轉(zhuǎn)化研究，構(gòu)建地表水—地下水水源地聯(lián)合風險預警平臺，加大地表水—地下水污染防治科研投入，多管齊下，才能保障飲水安全。

4 結(jié)論

（1）較短時空尺度的水循環(huán)模擬往往無法反映長期水循環(huán)規(guī)律，長期規(guī)律運行下的水循環(huán)演變規(guī)律需要在較長時空尺度下進行綜合模擬才具有意義。隨著人類活動的不斷加強，人類活動對流域整體水循環(huán)規(guī)律的影響是不容忽視的，因而流域尺度的整體水循環(huán)模擬必須綜合考慮人類活動的影響。

（2）地表水地下水聯(lián)合模擬調(diào)度決策離不開數(shù)學模型的支持，分布式流域水文模型對一個地區(qū)地表水和地下水相互作用的耦合模擬具有較高的精度，采分布式水文模型對較長時空跨度的地表水與地下水水文循環(huán)進行耦合應是未來水文模擬的重要方向。

（3）以流域為尺度的地表水地下水聯(lián)合模擬調(diào)度可以為水資源的合理開發(fā)利用提供決策支持，也將為治理水污染提供新的思路和方法，隨著水資源緊缺程度的加深和水安全問題的突出，地表水地下水聯(lián)合模擬調(diào)度將成為提高水資管理水平、保障水資源安全的重要手段。

（4）總體來講，國外多以模型的應用和優(yōu)化為地表水地下水的可持續(xù)利用和高效管理提供決策支持，國內(nèi)則以工程措施和法律法規(guī)及信息平臺建設為地表水地下水的聯(lián)合調(diào)度與高效管理提供參考。

5 參考文獻

[1] 顧大釗，張勇，曹志國.我國煤炭開采水資源保護利用技術(shù)研究進展[J].煤炭科學技術(shù)，2016（1）：1-7.

[2] 徐紅甫.我國缺水地區(qū)中小城市水資源可持續(xù)開發(fā)的現(xiàn)狀與對策[J].水利發(fā)展研究，2007（11）：48-50.

[3] 陳文科.我國重點流域地區(qū)的水污染趨向與治水治污路徑[J].江漢論壇，2010（2）：5-14.

[4] 趙艷民，秦延文，鄭丙輝，等.突發(fā)性水污染事故應急健康風險評價[J].中國環(huán)境科學， 2014（5）：1328-1335.

[5] 張新鈺，辛寶東，王曉紅，等.我國地下水污染研究進展[J].地球與環(huán)境，2011（3）：415-422.

[6] 孫棟元，伊力哈木，馮省利，等.干旱內(nèi)陸河流域地表水地下水聯(lián)合調(diào)度研究進展[J].地理科學進展，2009（2）：167-173.

[7] 王喜華.三江平原地下水：地表水聯(lián)合模擬與調(diào)控研究[D].北京：中國科學院，2015.

[8] SPANOUDAKI K，STAMOU AI，NANOU-GIANNAROU A.Development and verification of a 3-D integrated surface water-groundwater model[J].Journal of Hydrology，2009，375（3/4）：410-427.

[9] SOPHOCLEOUS M.Interactions between groundwater and surface water：the state of the science[J].Hydrogeology Journal，2002，10（1）：52-67.

[10] 王禹.礦坑地下水流量計算[J].有色金屬（冶煉部分），1954（7）：6-10.

[11] 汪東云.地下水過量開采與人工補給[J].建筑知識，1981（4）：42-43.

[12] FRANK G，陳恩健.肥料使用與地表水和地下水污染的關(guān)系[J].土壤學進展，1979（3）：26-37.

[13] KAZEZYELMAZ-ALHAN CM，MEDINA MA.The effect of surface/gro-und water interactions on wetland sites with different characteristics[J].Desalination，2008，226（1/2/3）：298-305.

[14] 朱金峰，王忠靜，鄭航，等.黑河流域中下游全境地表-地下水耦合模型與應用[J].中國環(huán)境科學，2015（9）：2820-2826.

[15] 姚俊強，楊青，陳亞寧，等.西北干旱區(qū)氣候變化及其對生態(tài)環(huán)境影響[J].生態(tài)學雜志， 2013（5）：1283-1291.

[16] 陳亞寧.氣候變化對西北干旱區(qū)水循環(huán)影響機理與水資源安全研究[J].中國基礎(chǔ)科學，2015（2）：15-21.

[17] 陳亞寧，楊青，羅毅，等.西北干旱區(qū)水資源問題研究思考[J].干旱區(qū)地理，2012（1）：1-9.endprint

[18] 陳亞寧，李稚，范煜婷，等.西北干旱區(qū)氣候變化對水文水資源影響研究進展[J].地理學報，2014（9）：1295-1304.

[19] 雍會，張風麗，張曉莉，等.干旱區(qū)塔里木河流域地下水資源統(tǒng)籌管理研究[J].中國農(nóng)業(yè)資源與區(qū)劃，2016（1）：90-94.

[20] 高業(yè)新，王貴玲，劉花臺，等.石羊河流域的水化學特征及其地表水與地下水的相互轉(zhuǎn)化[J].干旱區(qū)資源與環(huán)境，2006（6）：84-88.

[21] 駱東玲，陸培東，喻國華.東部平原城市水污染現(xiàn)狀及河網(wǎng)結(jié)構(gòu)對水污染影響[J].環(huán)境工程，2011（增刊1）：59-62.

[22] 張兵，宋獻方，張應華，等.第二松花江流域地表水與地下水相互關(guān)系[J].水科學進展 2014（3）：336-347.

[23] 張龍江，朱維斌，陶然，等.濉河符離集段地表水污染與傍河淺層地下水相互關(guān)系的研究[J].地下水，2006（5）：44-46.

[24] 王中根，朱新軍，李尉，等.海河流域地表水與地下水耦合模擬[J].地理科學進展，2011（11）：1345-1353.

[25] 孔曉樂，王仕琴，劉丙霞，等.外來調(diào)水對華北低平原區(qū)地表水和地下水水化學特征的影響：以河北省南皮縣為例[J].中國生態(tài)農(nóng)業(yè)學報，2016（8）：1135-1144.

[26] 羅美，王鎮(zhèn)，郭明仁.貴州普定喀斯特三水（雨水、地表水、地下水）資源開發(fā)利用與合理配置技術(shù)[J].安徽農(nóng)業(yè)科學，2014（20）：6662-6664.

[27] 張鳳太，易武英.典型喀斯特峰叢洼地地區(qū)水資源安全度評價：以貴州平塘縣為例[J].貴州科學，2016（5）：39-43.

[28] 周冰洋，趙兵，張曉詩，等.黔北大方地區(qū)地下水及地表水化學類型及特征分析[J].四川建材，2016（1）：243-245.

[29] 郎赟超，劉叢強，H S，等.貴陽地表水—地下水的硫和氯同位素組成特征及其污染物示蹤意義[J].地球科學進展，2008（2）：151-159.

[30] FLECKENSTEIN J H，KRAUSE S，HANNAH D M，et al.Groundwat-ersurface water interactions：New methods and models to improve unde-rstanding of processes and dynamics[J].Advances in Water Resources，2010，33（11）：1291-1295.

[31] 王蕊，王中根，夏軍.地表水和地下水耦合模型研究進展[J].地理科學進展，2008（4）：37-41.

[32] 胡立堂，王忠靜，趙建世，等.地表水和地下水相互作用及集成模型研究[J].水利學報，2007（1）：54-59.

[33] 殷禹宇，胡友彪，劉啟蒙，等.地表水與地下水相互作用研究進展[J].綠色科技，2016（4）：50-52.

[34] LAROCQUE M，COOK P G，HAAKEN K，et al.Estimating flow using tr-acers and hydraulics in synthetic heterogeneous aquifers[J].Ground water，2009，47（6）：786-796.

[35] KOLLET S J，MAXWELL R M.Integrated surface-groundwater flow modeling：A free-surface overland flow boundary condition in a parallel groundwater flow model[J].Advances in Water Resources，2006，29（7）：945-958.

[36] CHEN X.Hydrologic connections of a stream-aquifer-vegetation zone in south-central Platte River valley，Nebraska[J].Journal of Hydrology，2007，333（2/3/4）：554-568.

[37] WERNER A D，GALLAGHER M R，WEEKS S W.Regional-scale，fully coupled modelling of stream-aquifer interaction in a tropical catchment[J].Journal of Hydrology 2006，328（3/4）：497-510.

[38] ENGDAHL N B，VOGLER E T，WEISSMANN G S.Evaluation of aquifer heterogeneity effects on river flow loss using a transition probability framework[J].Water Resources Research，2010，46（1）：183-186.

[39] FLECKENSTEIN J H，NISWONGER R G，F(xiàn)OGG G E.River-aquifer interactions，geologic heterogeneity，and low-flow management[J].Gr-ound Water，2006，44（6）：837-852.

[40] 謝建文，鄭曉靜.地下水與地表水相互作用的識別和量化方法研究[J].海峽科技與產(chǎn)業(yè)，2016（8）：128-129.endprint

[41] 宋獻方，劉相超，夏軍，等.基于環(huán)境同位素技術(shù)的懷沙河流域地表水和地下水轉(zhuǎn)化關(guān)系研究[J].中國科學（D輯：地球科學），2007（1）：102-110.

[42] 楊改強，郭萍，李睿環(huán)，等.基于排隊理論的灌區(qū)渠系地表水及地下

水優(yōu)化配置模型[J].農(nóng)業(yè)工程學報，2016（6）：115-120.

[43] 白樂，李懷恩，何宏謀，等.煤礦開采區(qū)地表水-地下水耦合模擬[J].煤炭學報，2015（4）：931-937.

[44] 劉路廣，崔遠來.灌區(qū)地表水—地下水耦合模型的構(gòu)建[J].水利學報，2012（7）：826-833.

[45] 孔珂，徐晶，王昕，等.小開河灌區(qū)地表水沙及地下水聯(lián)合優(yōu)化配置模型[J].中國農(nóng)村水利水電，2016（1）：71-74.

[46] KANG S，SU X，TONG L，et al.The impacts of human activities on the water-land environment of the Shiyang River basin，an arid region in northwest China[J].Hydrological Sciences Journal，2004，49（3）：413-427.

[47] 鄧銘江，裴建生，王智.干旱區(qū)內(nèi)陸河流域地貌單元特征及地下水儲水構(gòu)造[J].水利學報，2006（11）：1360-1366.

[48] WELSH W D，VAZE J，DUTTA D，RASSAM D，et al：An integrated mo-delling framework for regulated river systems[J].Environmental Modelling & Software，2013，39：81-102.

[49] GORLA L，PREONA P.On quantifying ecologically sustainable flow rel-eases in a diverted river reach[J].Journal of Hydrology，2013，489：98-107.

[50] 陳志兵，潘曉東，唐建生.貴州省重點巖溶流域烏江下游段地下水賦存特征及開發(fā)利用方式探討[J].安徽農(nóng)業(yè)科學，2010（32）：18292-18293.

[51] 井柳新，孫愿平，孫宏亮，等.中國地表水-地下水污染協(xié)同管理控制模式初探[J].環(huán)境污染與防治，2016（3）：95-98.endprint