地表水與地下水相互作用研究進(jìn)展

朱金峰,劉悅憶,章樹安,鄭 航

地表水與地下水相互作用研究進(jìn)展

朱金峰1,2,劉悅憶2,3,章樹安1,鄭 航2,4*

(1.水利部水文局,北京 100053；2.清華大學(xué)水利水電工程系,水沙科學(xué)與水利水電工程國家重點(diǎn)實(shí)驗(yàn)室,北京 100084；3.南水北調(diào)東線總公司,北京 100088；4.東莞理工學(xué)院,東莞 523015)

地表水與地下水的水量轉(zhuǎn)化是水文循環(huán)的重要過程,地表水與地下水的相互作用及轉(zhuǎn)化關(guān)系是水文地質(zhì)等領(lǐng)域研究的熱點(diǎn)和難點(diǎn).為分析和掌握其規(guī)律特征,更好的進(jìn)行水資源分析評價、物質(zhì)遷移和能量變化的相關(guān)研究,本文對國內(nèi)外近20年的地表地下水循環(huán)及其相互作用的相關(guān)研究進(jìn)行了系統(tǒng)的分析和探討.通過梳理自然因子和人類活動影響下地表水與地下水相互作用的研究熱點(diǎn)及前沿問題,總結(jié)了該研究領(lǐng)域發(fā)展的主要研究方法與關(guān)鍵模擬技術(shù).通過有關(guān)研究方法的分析對比,當(dāng)前開始更多結(jié)合新技術(shù)對不同實(shí)驗(yàn)方法、解析法和數(shù)值模型的綜合研究,耦合自然因子和人類活動影響的流域地表水和地下水相互作用的研究成為研究的前沿和熱點(diǎn).

地表水與地下水；水循環(huán)；模型技術(shù)；水資源管理

地表水和地下水是相互關(guān)聯(lián)的水文連續(xù)體,存在于山區(qū)巖層、河流系統(tǒng)、海岸和喀斯特等地形環(huán)境中,是流域水文循環(huán)和水量收支計(jì)算必須考慮的要素[1-2].由于地表水和地下水的水量、熱量和物質(zhì)交換維持著河流生態(tài)系統(tǒng)的基本功能,對于流域的水資源管理保護(hù)、水污染防治和生態(tài)健康具有重要意義[3-4].水質(zhì)方面,該過程影響著水化學(xué)成分的分布和演變規(guī)律[5];水量方面,地下水是一些流域水文循環(huán)和水資源轉(zhuǎn)化的主導(dǎo)因素,尤其對一些干旱區(qū)的生態(tài)植被狀況影響顯著.在降水稀少地區(qū),河床滲漏對地下水補(bǔ)給占有較大比例,枯水期地下水含水層以基流的形式向河道排泄,從而保證河道流量和維持生態(tài)系統(tǒng)[6-7].隨著人類社會和經(jīng)濟(jì)活動的發(fā)展,許多河流水利工程的修建和運(yùn)行,地下水開采和灌溉渠系管網(wǎng)系統(tǒng)的完善,在解決流域水資源時空分布不均、水資源供需矛盾和水資源合理配置等方面發(fā)揮著巨大作用,同時干擾了地表水和地下水系統(tǒng),灌溉引水和地下水開采改變了地表地下水的循環(huán),導(dǎo)致河道徑流減少和生態(tài)退化等生態(tài)環(huán)境問題[8-9].然而,歷史上水資源管理更多著眼于獨(dú)立的地表水系統(tǒng)或地下水資源的開發(fā)利用,地表水與地下水的交換轉(zhuǎn)化常被忽視.污染物的運(yùn)移和擴(kuò)散也受到河流-含水層相互作用的密切影響,考慮地表地下水的水量交換、揭示流域水循環(huán)的演變特征是水資源科學(xué)管理的重要內(nèi)容[10].

由于地表水和地下水相互作用涉及的領(lǐng)域和學(xué)科比較廣泛,近些年關(guān)于其規(guī)律特征的系統(tǒng)研究較少,對相關(guān)的研究成果、研究內(nèi)容和技術(shù)方法缺乏系統(tǒng)的梳理;雖然有關(guān)的理論研究已經(jīng)比較完善,部分新的研究方法和技術(shù)還有待發(fā)展,相關(guān)的基礎(chǔ)科學(xué)研究與水資源保護(hù)管理的實(shí)踐應(yīng)用也需要更好的結(jié)合.首先,從研究內(nèi)容上進(jìn)行揭示,該領(lǐng)域的研究熱點(diǎn)主要有地表水和地下水循環(huán)的時空過程、河床下層的“底流區(qū)”的水量和物質(zhì)能量交換、不同自然因素和人類活動對河流-含水層水量交換的定量影響評價等;其次,從研究方法上總結(jié)該領(lǐng)域的幾類研究方法(實(shí)驗(yàn)觀測和數(shù)學(xué)模型)和模型技術(shù)(解析法和數(shù)值模型),分析了其主要特點(diǎn)和發(fā)展應(yīng)用情況;最后,通過該領(lǐng)域有關(guān)研究內(nèi)容和研究方法的分析對比,提出了當(dāng)前地表水與地下水研究的主要問題和挑戰(zhàn),以及未來該領(lǐng)域研究的發(fā)展方向.

1 地表水-地下水循環(huán)轉(zhuǎn)化規(guī)律及其影響因素

1.1 地表水與地下水循環(huán)形式

地表水與地下水相互作用的核心是地表水和地下水循環(huán),準(zhǔn)確掌握其機(jī)理和過程是支撐流域水資源綜合利用和保護(hù)的基礎(chǔ)條件.地表地下水循環(huán)同水資源管理密切相關(guān),它們的定量研究是流域尺度水量收支的組成部分,也是地表水和地下水聯(lián)合調(diào)度的重要內(nèi)容[11-12].從20世紀(jì)90年代,水文地質(zhì)學(xué)家開始更多關(guān)注近河道的地表地下水交換,相關(guān)研究在近20年持續(xù)增長[13].一般而言,地形高的地帶是地下水補(bǔ)給區(qū),而地形低洼區(qū)產(chǎn)生地下水排泄[14].地表水和地下水補(bǔ)給-排泄過程機(jī)理涉及多維度多因素、空間變異性和水力特性、不同尺度和學(xué)科的綜合影響,存在許多復(fù)雜的水文過程,本質(zhì)上都是通過地表水入滲和地下水排泄實(shí)現(xiàn)的,表現(xiàn)為降水補(bǔ)給、河道基流、地表水體滲漏(湖泊、河道、渠道等)等水循環(huán)現(xiàn)象[15].其中,地表水與地下水的轉(zhuǎn)化包括湖泊與地下水、地下水與濕地/海岸、河流-含水層系統(tǒng)、泉水排泄等幾種形式,以及人類活動引發(fā)的地表水與地下水的轉(zhuǎn)化[8,16].

河流-含水層系統(tǒng)的相互作用是流域地表地下水循環(huán)中的核心內(nèi)容.作為沖積含水層地下水排泄的主要過程,河流不僅是地下水流出和匯入的地方,也是濱河帶和河流生態(tài)系統(tǒng)的重要組成部分[8,17].相關(guān)定量研究始于一百多年前Boussinesq對河流與連續(xù)沖積含水層作用規(guī)律的探討[1].近些年關(guān)于河流范圍內(nèi)地表水和地下水相互影響的研究,主要是為河流系統(tǒng)生態(tài)結(jié)構(gòu)評價、河流恢復(fù)和濱河綠洲濕地管理等服務(wù)[18-19].在沖積平原地區(qū),受地形和水文特性因素影響,河流同地下含水層的聯(lián)系表現(xiàn)為獲取型、損耗型、穿透流動和平行流動幾種形式,并且在時間和空間上動態(tài)變化[4,8,20].

河流-含水層系統(tǒng)水循環(huán)研究的熱點(diǎn)和前沿是確定河流-含水層水量交換的時空分布,它們受自然和人類等影響因子的控制.在系統(tǒng)內(nèi)部,沖積平原的水流、運(yùn)移和地表地下水交換主要由三類因素控制:含水層和河床特性、水力梯度、河流的位置和結(jié)構(gòu)等,它們共同決定著地表地下水交換的水量、方向和空間分布等.在系統(tǒng)外部,氣候因子和人類水資源開發(fā)活動則間接影響著該系統(tǒng)演化[14].

1.2 自然影響因子

在河流-含水層系統(tǒng)的相關(guān)研究中,綜合分析影響其水循環(huán)的自然因子,包括流域地質(zhì)地貌、河床特性和水文情勢等,它們在系統(tǒng)內(nèi)部驅(qū)動著地表水-地下水循環(huán),改變著地表水與地下水的轉(zhuǎn)化,有關(guān)的定量研究一直是熱點(diǎn)領(lǐng)域.

(1)流域地質(zhì)地貌:對于地下水含水層而言,從河流與地下水關(guān)系的宏觀角度,流域的地質(zhì)地貌因素控制著其相互關(guān)系和演化的基本格局,如許多天然河流從上游河谷深切、中游河流非完整切割含水層、下游平原河谷開闊等呈現(xiàn)出顯著不同的河流-含水層系統(tǒng)演化過程[21].對于宏觀不同的地形地貌和地質(zhì)構(gòu)造等因素對河流-含水層系統(tǒng)、水量轉(zhuǎn)化的影響,王文科等以黃河流域?yàn)槔?首次系統(tǒng)深入研究了河水與地下水轉(zhuǎn)化的8種模式類型,對進(jìn)一步研究河流-含水層的相互關(guān)系和定量模擬分析具有重要意義[22].在一定的地貌條件下,宏觀的地質(zhì)特性表現(xiàn)為含水層的滲透系數(shù)、給水度和地表高程等微觀描述,它們在面狀空間范圍影響著地表地下水循環(huán),是決定流域尺度地下水運(yùn)動和地表水文過程的重要因子[23-24].

不同的顆粒組成、巖性礦物和成巖等決定了含水層具有非均質(zhì)性,表現(xiàn)為影響地下水流的水力參數(shù)情況,非均質(zhì)程度會影響到地下水水位、河流-含水層的連通狀態(tài)[25-26].在一些地表水和地下水交換強(qiáng)烈的區(qū)域,含水層非均質(zhì)性被認(rèn)為是重要的不確定因素[27-28].具體來說,地下水運(yùn)動很大程度上受到值的大小及非均勻分布控制,而給水度的影響作用較小[14,29],值范圍跨度對河流-含水層水量交換的影響也較為明顯[30].通過模型的模擬分析,能夠定量研究含水層非均質(zhì)性對虧損型河流向含水層補(bǔ)給量(大小及空間分布)的影響[26,31].在數(shù)學(xué)意義上,方差和相關(guān)長度能夠刻畫空間分布的非均質(zhì)性,河流-含水層水量轉(zhuǎn)化的不確定性隨非均質(zhì)性的增加而增大[26,32].

(2)河床特性:對于河道而言,在河流-含水層系統(tǒng)中,河床特性主要是指河床巖性相關(guān)的物理參數(shù).在河道范圍內(nèi),地表水與地下水的流動通過河床下層的“底流區(qū)”產(chǎn)生復(fù)雜的作用,維持著底棲動物、微生物生存和物理化學(xué)、水文地質(zhì)等復(fù)雜過程,近些年許多野外試驗(yàn)和實(shí)驗(yàn)室研究的熱點(diǎn)[33-34].有學(xué)者通過定量研究模擬包氣帶水流和污染物的運(yùn)移對地下水的影響,地下水沿河床排泄的強(qiáng)度影響著過渡區(qū)的物質(zhì)、水量交換和面積,也影響著地表河道的水文徑流[35-36].

河流-含水層的水量主要通過河床進(jìn)行轉(zhuǎn)化,河床介質(zhì)和滲透速率直接影響到其相互作用,如有研究對天山北麓瑪納斯河從上游到下游入滲速率的變化研究,表明其與河床巖性由粗變細(xì)一致[21].許多定量的研究和分析表明,河床水力滲透系數(shù)、曼寧糙率系數(shù)、坡度和河流寬度均對河流和地下水循環(huán)具有重要影響[37-39].河床的非均質(zhì)性及多期沉積物也會影響到地表水與地下水交換的時空分布、地下水運(yùn)動和溶質(zhì)運(yùn)移情況等,在相關(guān)的研究模擬中考慮較少.

(3)水文情勢:河道徑流和地下水位的變化會直接改變河流-含水層的補(bǔ)排關(guān)系[21].如河流水位和地下水位相互關(guān)系的季節(jié)性變化,不但影響水力坡度的大小,有時甚至能夠影響到水量交換的方向,引起地表地下水循環(huán)復(fù)雜的更替演化[22].河流水位對河流-含水層水量交換影響可能不大,但是隨著河道徑流情勢和地下水位的相對變化,含水層水位和河流間的水力梯度在不同位置會發(fā)生季節(jié)性改變,從而直接影響到近河床的地下水運(yùn)動、河水下滲和地下水排泄量[28-29].

1.3 人類活動影響

在水資源緊缺的干旱半干旱流域,人類活動對水資源系統(tǒng)的影響更為明顯,引起了一些流域地表水和地下水交換量的下降,改變了水文循環(huán)過程[40].人類在社會經(jīng)濟(jì)發(fā)展中通過水資源開發(fā)利用,干擾到天然狀態(tài)的地表地下水循環(huán),進(jìn)而改變著河流與地下含水層的補(bǔ)排關(guān)系.例如,河道引水會減少河道徑流、河流向地下水的補(bǔ)給;依賴地下水灌溉的灌區(qū),地下水開采、渠系滲漏和田間入滲等過程間接改變了地表水和地下水的運(yùn)動,部分地區(qū)的渠系滲漏和田間入滲已成為地下水的主要補(bǔ)給形式.特別是在一些干旱地區(qū),人類活動引起的土地利用和渠系滲漏的變化對地表地下水的時空分布影響強(qiáng)烈[41].地下水開采被認(rèn)為是影響地下水位最敏感的因子,它通過含水層補(bǔ)給量增加、原有自然排泄量減少或含水層儲量減小等途徑進(jìn)行平衡,稱為對其他水源的“捕捉”,引起河道徑流和地下水蒸發(fā)減少等現(xiàn)象,從而影響到長期的生態(tài)水文演變[42].

2 主要研究方法

為改善河道徑流和維持地下水位,合理地模擬地表水-地下水循環(huán)對于水文過程演變分析和水量合理調(diào)度十分重要[43-45].隨著社會經(jīng)濟(jì)的發(fā)展和人類活動的進(jìn)一步增強(qiáng),在水資源的開發(fā)利用與調(diào)度配置中開展地表水與地下水的相互作用研究十分必要.對于河流-含水層系統(tǒng),地表水和地下水水量交換的研究方法有直接水量測量、間接實(shí)驗(yàn)法、水量平衡法和基于達(dá)西定律的研究方法[4,13,46-47].主要研究方法的比較見表1.

表1 地表水和地下水相互作用主要研究方法對比表

2.1 直接水量測量

利用滲漏測量儀對含水層-河床界面的點(diǎn)尺度水量通量進(jìn)行觀測,能夠確定河流滲漏和地下水排泄量.但是河床平均滲漏量的確定需要很多觀測點(diǎn)的數(shù)據(jù),測得的也是點(diǎn)尺度數(shù)據(jù),較大尺度的數(shù)據(jù)獲取較困難,在實(shí)際的研究應(yīng)用相對較少[4].

2.2 間接實(shí)驗(yàn)法

利用熱力學(xué)或同位素化學(xué)等方法進(jìn)行地表水和地下水交換水量的推算[2,48].由于地下水的水溫是相對恒定的,而河流水溫一般具有季節(jié)性變化規(guī)律,其水量交換過程可通過相關(guān)的溫度觀測和熱力學(xué)方法進(jìn)行建模推算.Conant[49]最早指出多測點(diǎn)短期的河床溫度觀測同地下水排泄的空間分布密切相關(guān),建立了溫度分布和河床垂向水量的簡單數(shù)學(xué)模型,但是其地下水排泄量的計(jì)算仍是基于達(dá)西定律和實(shí)測水頭梯度,在應(yīng)用中也需要很多假設(shè)條件.Hatch等[50]提出了采用時間序列熱力學(xué)數(shù)據(jù)確定河床滲漏量的方法,該法具有較少的數(shù)據(jù)需求和較快的速度,能夠分析復(fù)雜地表水和地下水交換的時空變化動態(tài).基于實(shí)測溫度對地表水和地下水交換量的研究仍局限于均一性和邊界簡化地區(qū),且相關(guān)的熱力學(xué)參數(shù)依賴于估計(jì)值,對地表水量的變化響應(yīng)也難以考慮[51].

水循環(huán)過程中的水化學(xué)成分一定程度上記錄著水分運(yùn)移、轉(zhuǎn)化的歷史,水化學(xué)分析方法也是了解地下水與地表水相互作用的有效方法之一[52].通過水化學(xué)分析識別河流不同地帶地表地下水的相互轉(zhuǎn)化關(guān)系,能夠給出不同水源的組成比例及循環(huán)轉(zhuǎn)化特征[53-54],如采用水文地質(zhì)化學(xué)分析方法研究地下水補(bǔ)給的水源和機(jī)制,分析潛水補(bǔ)給的主要來源[55].結(jié)合同位素技術(shù)進(jìn)行量化分析,能夠?qū)Φ叵滤a(bǔ)給源、變化特征,以及人類活動和氣候變化均對地表地下水情勢的影響進(jìn)行定量評價[56],如通過利用穩(wěn)定同位素追蹤水流運(yùn)動研究河流水量的來源組成[3],利用氚模型模擬數(shù)十年尺度的地表水和地下水交換,估算多年平均的地下水補(bǔ)給和排泄量[15].

2.3 水量平衡法

水量平衡法又分為多斷面流量觀測、河道徑流分割和物質(zhì)追蹤等,借助離子追蹤、同位素技術(shù)等確定地表水和地下水交界處的水量交換.多斷面流量觀測通過不同斷面的流量觀測確定河道的凈流入或凈流出[48].河道徑流分割假定基流完全代表地下水排泄量,在蒸散發(fā)量等水均衡項(xiàng)資料缺乏時具有較大不確定性[57];物質(zhì)追蹤根據(jù)不同水源的濃度不同確定河水和地下水的貢獻(xiàn)量,詳見文獻(xiàn)[4].另外,還有基于長系列地表水和地下水觀測資料,利用水文系統(tǒng)分析法對地表水和地下水相互轉(zhuǎn)化規(guī)律進(jìn)行研究[58],以及依據(jù)已有的地表水(流量)和地下水(水位)觀測數(shù)據(jù),結(jié)合調(diào)研和用水?dāng)?shù)據(jù)的整理成果,對流域內(nèi)水量的收支平衡進(jìn)行計(jì)算分析,推算地表水與地下水轉(zhuǎn)化的規(guī)律[59-60].

3 主要模擬技術(shù)

這些實(shí)驗(yàn)觀測的研究方法常用于模型建立和結(jié)果評價的基礎(chǔ)依據(jù),但是其定性分析較多,且傳統(tǒng)的觀測手段需要大量的人力和物力,實(shí)際研究中難以真實(shí)刻畫地表地下水轉(zhuǎn)化的物理過程,隨著許多新技術(shù)和儀器的發(fā)展應(yīng)用,新的觀測手段能夠更好的提高實(shí)驗(yàn)觀測法的水平和應(yīng)用能力.其中,基于達(dá)西定律的研究方法應(yīng)用非常廣泛,它一般結(jié)合壓力水頭觀測和解析法、數(shù)值模型法等技術(shù)進(jìn)行模擬研究[13,46].

3.1 解析法

在數(shù)值模型出現(xiàn)以前,受計(jì)算能力的限制,河流-含水層相互的水量轉(zhuǎn)化更多的是依靠地下水運(yùn)動的相關(guān)解析模型進(jìn)行計(jì)算,包括求解一維和二維的Boussinesq方程,并一直沿用和發(fā)展[61-62].該方法主要是在自由表面邊界條件假定下,借助拉普拉斯變換或降解法得到均勻各向同性條件的地下水結(jié)果,并利用地下水位觀測值進(jìn)行驗(yàn)證.解析法根據(jù)河流和含水層類型(河流切割程度、水力聯(lián)系等),將河流-含水層的解析模型分為三種類型[46].

該方法多用于簡單或假設(shè)邊界條件的河流-含水層系統(tǒng)等特定條件下的求解,對于河流水管理者的情景分析和短期的水位預(yù)報(bào)比較有效和方便[11,63-64].但是難以代表流域水文的真實(shí)物理過程,也無法體現(xiàn)相關(guān)水文參數(shù)的空間異質(zhì)性,有學(xué)者發(fā)現(xiàn)它會高估河流的滲漏量,并且由于解析法采用過多的簡化條件而存在較大誤差,使得解析解模型在實(shí)際應(yīng)用時存在許多限制.由于解析法能夠便于分析不同邊界條件和影響因子下徑流及地下水的響應(yīng),其結(jié)果可用于一些數(shù)值模型的驗(yàn)證,通過兩種方法的結(jié)合應(yīng)用能夠提高地表水和地下水運(yùn)動的模擬能力和可靠性.

3.2 數(shù)值模型法

隨著計(jì)算機(jī)運(yùn)算能力的提高,利用數(shù)值模型成為研究流域尺度地表地下水循環(huán)的重要方法,它能夠快速處理更為復(fù)雜和現(xiàn)實(shí)的水文地質(zhì)條件,比解析法更具有優(yōu)勢和廣泛的應(yīng)用[42].20世紀(jì)60年代,計(jì)算機(jī)的應(yīng)用和數(shù)值模型的出現(xiàn)成為水文地質(zhì)定量研究進(jìn)展的重要突破.在70年代前,地下水?dāng)?shù)值模型基本是穩(wěn)定流計(jì)算,隨后出現(xiàn)了基于有限元計(jì)算的綜合的地表地下水模型和瞬變流的模擬模型[65].用于模擬河流-含水層系統(tǒng)的數(shù)值模型于80年代開始建立并迅速發(fā)展,早期主要用于水量管理和聯(lián)合調(diào)度情景的評價,但是實(shí)際受到調(diào)查資料和計(jì)算能力限制導(dǎo)致時空離散比較粗糙,河流地表水運(yùn)動、河床和含水層等物理?xiàng)l件均存在很大簡化[13].隨著計(jì)算機(jī)運(yùn)算和存儲能力的提高,數(shù)值模型方法不再受研究尺度和復(fù)雜邊界條件限制,在物理過程和參數(shù)處理方面具有較高靈活性,成為研究地表地下水循環(huán)、相互作用及水資源管理的主要手段[66-67].

綜合地表地下水循環(huán)的數(shù)值模型在流域水循環(huán)和水資源分析評價的研究中應(yīng)用廣泛,它們能分布式地模擬流域地表水文過程及地下水的運(yùn)動,以及地表水與地下水的交換.地表水模型是利用相關(guān)經(jīng)驗(yàn)公式、明渠圣維南方程等建立的,地下水模型一般是是根據(jù)達(dá)西定律和水量平衡原理進(jìn)行計(jì)算.關(guān)于地表水流動、地下水流運(yùn)動及其耦合方法等方面的系統(tǒng)研究,Spanoudaki 等[68]對已有的地表地下水綜合模型進(jìn)行了系統(tǒng)總結(jié).地表地下水循環(huán)的耦合模型從地下水結(jié)構(gòu)上,分為單層平面模型、雙層準(zhǔn)三維模型和多層三維模型[69];從耦合方法上,不同方式對地表水和地下水模型進(jìn)行耦合或連接,分為完全耦合,迭代耦合和松散耦合[70].

在物理模擬中,河流-含水層水量轉(zhuǎn)化依賴河流水位和地下水位這兩個變量,它們一般分別由河道和地下水運(yùn)動方程計(jì)算得到,也是地表水和地下水模擬模型進(jìn)行連接和耦合計(jì)算的關(guān)鍵要素[71].目前,許多用于河流-含水層系統(tǒng)研究的數(shù)值模型主要是基于美國地質(zhì)調(diào)查局開發(fā)的MODFLOW[10,29,72],以及將其與其它成熟的地表水模擬模型進(jìn)行耦合,通過GIS平臺建立流域地表水與地下水計(jì)算單元的轉(zhuǎn)換關(guān)系,有較好的應(yīng)用[43,73].其他一些成熟的數(shù)值模型在地表水和地下水耦合模擬中應(yīng)用,如利用ParFlow開展的沖積地區(qū)的河流-含水層系統(tǒng)水量交換的空間特征和時間動態(tài)進(jìn)行了研究,結(jié)合地質(zhì)統(tǒng)計(jì)、高斯模擬和蒙特卡洛方法模擬不同尺度非均質(zhì)性的影響[74];利用HydroGeoSphere模擬不同含水層非均質(zhì)情景的地下水流動影響[32];利用GSFLOW能夠較好地模擬復(fù)雜地下水與地表水的相互作用過程,分析不同水文參數(shù)對于流域水資源的影響[75].

已有數(shù)值模型對河流-含水層的交換項(xiàng)和河道的概化通常比較簡化,如MODFLOW采用了矩形網(wǎng)格差分法,地表水流近似處理為一維明渠恒定流;對于地表水和地下水的交換水量當(dāng)作源匯項(xiàng)概化處理.武強(qiáng)等(2005)建立了一維明渠非恒定漸變流偏微分描述的地表河網(wǎng)和“擬三維”非穩(wěn)定流地下水模擬的耦合模型SGS[76]; Spanoudaki等[68]首次建立了耦合三維地表水和三維地下水的復(fù)雜模型(IRENE),采用半顯示和顯示有限差分求解地表水方程和地下水方程,能夠更為精確細(xì)致地模擬地表水文過程、地下水運(yùn)動及其相互作用.可見,隨著計(jì)算能力和新技術(shù)(如GIS和RS)的發(fā)展,數(shù)值模型作為強(qiáng)大的工具手段廣泛用于地表水和地下水模型的耦合研究,許多分布式的、多維度的河流-含水層水文循環(huán)過程的模擬也更多依賴數(shù)值模型進(jìn)行分析研究[77-79].

這些數(shù)值模型技術(shù)通?？山Y(jié)合解析法、其他數(shù)值模型或?qū)嶒?yàn)觀測數(shù)據(jù)進(jìn)行相互驗(yàn)證分析,提高模型的可靠性[80-81].在缺少準(zhǔn)確的參數(shù)分布和估計(jì)時,可通過敏感性分析和隨機(jī)模擬等方法研究河流-含水層的相互作用規(guī)律[82].雖然數(shù)值模型在精度、穩(wěn)定性或收斂性上可能存在一些問題,目前利用數(shù)值模型研究從微觀尺度、區(qū)域尺度到流域尺度的河流-含水層的循環(huán)與相互轉(zhuǎn)化規(guī)律仍是主要的研究方法,結(jié)合高精度的時空實(shí)驗(yàn)數(shù)據(jù)和數(shù)值物理模型研究地表水-地下水的機(jī)理過程,以及地表水-地下水耦合模型受非均質(zhì)性的影響研究仍面臨挑戰(zhàn).

4 結(jié)論與展望

4.1 結(jié)論

地表水和地下水有關(guān)的監(jiān)測和實(shí)驗(yàn)是地表水與地下水相互作用研究的重要基礎(chǔ),采用直接水量測量、間接實(shí)驗(yàn)法、水量平衡法和基于達(dá)西定律的研究方法仍在結(jié)合新技術(shù)手段不斷發(fā)展.當(dāng)前開始更多結(jié)合新技術(shù)對不同實(shí)驗(yàn)方法、解析法和數(shù)值模型的綜合研究,耦合自然因子和人類活動影響的流域地表水和地下水相互作用的研究成為研究的前沿和熱點(diǎn).美國地質(zhì)調(diào)查局的懷俄明蒙特拿水科學(xué)中心開展了Fish河流區(qū)域的地表水-地下水作用監(jiān)測項(xiàng)目,通過在河流水文監(jiān)測信息納入附近實(shí)時的地下水位和水溫監(jiān)測信息,擴(kuò)大監(jiān)測功能,系統(tǒng)采集長系列的數(shù)據(jù)資料,研究氣候變化、灌溉和水文徑流等的響應(yīng)關(guān)系. 地下水信息的監(jiān)測日益受到重視,我國從2015年開始了國家地下水監(jiān)測工程項(xiàng)目,由水利部和國土資源部組織實(shí)施,通過完善地下水監(jiān)測網(wǎng)絡(luò),基本覆蓋全國平原、盆地和巖溶山區(qū)350萬km2的地下水動態(tài)有效監(jiān)控,其中部分泉流量等的監(jiān)測信息,將來可為地表水和地下水相互作用的研究提供基礎(chǔ)支撐.

4.2 展望

研究內(nèi)容上,人類活動和地表地下水循環(huán)強(qiáng)烈的流域,對地表地下水循環(huán)的演變和水資源的調(diào)度管理有待更長時空尺度的綜合模擬.已有研究更多著眼于自然因子影響的模擬分析,對人類通過不同的水資源開發(fā)利用方式進(jìn)行的河流水量調(diào)配、地下水開采和水資源調(diào)度管理,需要研究綜合地表水的循環(huán)利用、地下水系統(tǒng)的演變和人類活動的相互影響等,以及不同要素對河流-含水層影響的定量評價方法.

技術(shù)方法上,綜合新技術(shù)應(yīng)用的實(shí)驗(yàn)觀測能增加對含水層地質(zhì)參數(shù)的理解,需要通過不同信息源數(shù)據(jù)的融合和分析研究,提高對地表水與地下水相互作用規(guī)律的認(rèn)識.地表水和地下水運(yùn)動的模擬模型通常需要建立在較大的空間尺度,受傳統(tǒng)人工監(jiān)測手段的影響,有限的觀測數(shù)據(jù)還難以用于模型的研究尺度.在較大尺度的地下水模擬中通過現(xiàn)場觀測和地質(zhì)調(diào)查數(shù)據(jù)得到參數(shù)分布仍比較困難且存在很大不確定性,將來結(jié)合國家地下水監(jiān)測工程項(xiàng)目的運(yùn)行和地下水監(jiān)測技術(shù)的提高,能更好的滿足不同時空尺度地下水觀測資料的需求.

模擬概化中,完善模型的物理機(jī)理和提高模型的模擬預(yù)測能力是重要的發(fā)展方向.由于數(shù)學(xué)解析方法受到復(fù)雜物理過程和基礎(chǔ)理論的限制發(fā)展較緩,需要研究同求解能力相適應(yīng)的概化方法,同時建立與數(shù)值模型相互驗(yàn)證結(jié)合的解析方法也具有重要前景,方便應(yīng)用于地表水-地下水轉(zhuǎn)化的預(yù)測分析.物理概化中,已有研究關(guān)于不同河流-含水層模式、邊界條件(非飽和帶地下水同地表水、飽和地下水的耦合)研究相對薄弱,其模擬還需要更細(xì)致考慮不同關(guān)系類型和非飽和帶地下水.

模型參數(shù)上,如何合理改善模型參數(shù)和不確定性研究具有重要意義.在不同的尺度上,數(shù)值模型從概念性和理論性出發(fā)廣泛用于河流-地下水系統(tǒng)的模擬分析,但目前研究多是根據(jù)水文地質(zhì)特性將研究區(qū)或河流劃分為不同單元,各單元具有相同的參數(shù)值.河床和含水層參數(shù)大小及分布的空間變化,使得確定性數(shù)值模型模擬的地表地下水循環(huán)具有較大不確定性.地表水-地下水的相互轉(zhuǎn)化能夠反映地表河道徑流和地下水流場的動力特征,而非均質(zhì)性對其影響的不確定性研究仍顯不足,發(fā)展隨機(jī)數(shù)值模型也是將來重要的研究方向.

[1] Winter T C. Recent advances in understanding the interaction of groundwater and surface water [J]. Reviews of Geophysics, 1995, 33(S2):985-994.

[2] Sophocleous M. Interactions between groundwater and surface water: the state of the science [J]. Hydrogeology Journal, 2002, 10(1):52-67.

[3] Lambs L. Interactions between groundwater and surface water at river banks and the confluence of rivers [J]. Journal of Hydrology, 2004,288(3):312-326.

[4] Kalbus E, Reinstorf F, Schirmer M. Measuring methods for groundwater–surface water interactions: a review [J]. Hydrology and Earth System Sciences, 2006,10(6):873-887.

[5] 張亞麗,張依章,張 遠(yuǎn),等.渾河流域地表水和地下水氮污染特征研[J]. 中國環(huán)境科學(xué), 2014,34(1):170-177.

[6] Hussein M, Schwartz F W. Modeling of flow and contaminant transport in coupled stream–aquifer systems [J]. Journal of contaminant hydrology, 2003,65(1):41-64.

[7] 武 強(qiáng),徐軍祥,張自忠,等.地表河網(wǎng)-地下水流系統(tǒng)耦合模擬Ⅱ:應(yīng)用實(shí)例[J]. 水利學(xué)報(bào), 2005,(6):754-758.

[8] Woessner W W. Stream and fluvial plain ground water interactions: rescaling hydrogeologic thought [J]. Ground Water, 2000,38(3):423-429.

[9] 陳永金,李衛(wèi)紅,陳亞寧,等.塔里木河流域綜合治理的生態(tài)效應(yīng)[J]. 中國環(huán)境科學(xué), 2007,27(1):24-28.

[10] 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):497-510.

[11] Hantush M M. Modeling stream–aquifer interactions with linear response functions [J]. Journal of Hydrology, 2005,311(1):59-79.

[12] 王文科,楊澤元,程東會,等.面向生態(tài)的干旱半干旱地區(qū)區(qū)域地下水資源評價的方法體系[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2011,(1):159-167.

[13] Fleckenstein J H, Krause S, Hannah D M, et al. Groundwater- surface water interactions: new methods and models to improve understanding of processes and dynamics [J]. Advances in Water Resources, 2010,33(11):1291-1295.

[14] Winter T C. Relation of streams, lakes, and wetlands to groundwater flow systems [J]. Hydrogeology Journal, 1999,7(1): 28-45.

[15] Harvey J W, Newlin J T, Krupa S L. Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA [J]. Journal of Hydrology, 2006, 320(3):400-420.

[16] Restrepo J I, Montoya A M, Obeysekera J. A wetland simulation module for the MODFLOW ground water model [J]. Ground Water, 1998,36(5):764-770.

[17] Rodríguez L B, Cello P A, Vionnet C A. Modeling stream-aquifer interactions in a shallow aquifer, Choele Choel Island, Patagonia, Argentina [J]. Hydrogeology Journal, 2006,14(4):591-602.

[18] Kazezyelmaz-Alhan C M, Medina M A. The effect of surface/ ground water interactions on wetland sites with different characteristics [J]. Desalination, 2008,226(1):298-305.

[19] Kazezyelmaz-Alhan C M, Medina M A, Richardson C J. A wetland hydrology and water quality model incorporating surface water/groundwater interactions [J]. Water Resources Research, 2007,43(4):W04434.

[20] Osman Y Z, Bruen M P. Modelling stream–aquifer seepage in an alluvial aquifer: an improved loosing-stream package for MODFLOW [J]. Journal of Hydrology, 2002,264(1):69-86.

[21] 王文科,李俊亭,王 釗,等.河流與地下水關(guān)系的演化及若干科學(xué)問題[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2007,(2):231-238.

[22] 王文科,孔金玲,段 磊,等.黃河流域河水與地下水轉(zhuǎn)化關(guān)系研究[J]. 中國科學(xué)E輯:技術(shù)科學(xué), 2004,(S1):23-33.

[23] Sophocleous M, Perkins S P. Calibrated models as management tools for stream-aquifer systems: the case of central Kansas, USA [J]. Journal of Hydrology, 1993,152(1):31-56.

[24] Vereecken H, D?ring U, Hardelauf H, et al. Analysis of solute transport in a heterogeneous aquifer: the Krauthausen field experiment [J]. Journal of contaminant Hydrology, 2000,45(3): 329-358.Cardenas M B.

[25] Kalbus E, Schmidt C, Molson J W, et al. Influence of aquifer and streambed heterogeneity on the distribution of groundwater discharge [J]. Hydrology and Earth System Sciences, 2009,13(1): 69-77.

[26] Bruen M P, Osman Y Z. Sensitivity of stream–aquifer seepage to spatial variability of the saturated hydraulic conductivity of the aquifer [J]. Journal of Hydrology, 2004,293(1):289-302.

[27] Kollet S J, Zlotnik V A. Stream depletion predictions using pumping test data from a heterogeneous stream–aquifer system (a case study from the Great Plains, USA) [J]. Journal of Hydrology, 2003,281(1):96-114.

[28] Derx J, Blaschke A P, Bl?schl G. Three-dimensional flow patterns at the river–aquifer interface—a case study at the Danube [J]. Advances in Water Resources, 2010,33(11):1375-1387.

[29] Wroblicky G J, Campana M E, Valett H M, et al. Seasonal variation in surface-subsurface water exchange and lateral hyporheic area of two stream-aquifer systems [J]. Water Resources Research, 1998,34(3):317-328.

[30] 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): W01506.

[31] Fleckenstein J H, Niswonger R G, Fogg G E. River-aquifer interactions, geologic heterogeneity, and low-flow management [J]. Ground Water, 2006,44(6):837-852.

[32] Larocque M, Cook P G, Haaken K, et al. Estimating flow using tracers and hydraulics in synthetic heterogeneous aquifers [J]. Ground water, 2009,47(6):786-796.

[33] Cardenas M B, Wilson J L, Zlotnik V A. Impact of heterogeneity, bed forms, and stream curvature on subchannel hyporheic exchange [J]. Water Resources Research, 2004,40(8):W08307.

[34] Salehin M, Packman A I, Paradis M. Hyporheic exchange with heterogeneous streambeds: Laboratory experiments and modeling [J]. Water Resources Research, 2004,40(11):W11504.

[35] 李 瑋,何江濤,劉麗雅,等.Hydus-1D軟件在地下水污染風(fēng)險(xiǎn)評價中的應(yīng)用[J]. 中國環(huán)境科學(xué), 2013,33(4):639-647.

[36] Stream-aquifer interactions and hyporheic exchange in gaining and losing sinuous streams [J]. Water Resources Research, 2009, 45(6):W06429.

[37] Chen X. Streambed hydraulic conductivity for rivers in south- central Nebraska 1 [J]. JAWRA Journal of the American Water Resources Association, 2004,40(3):561-573.

[38] Ryan R J, Boufadel M C. Evaluation of streambed hydraulic conductivity heterogeneity in an urban watershed [J]. Stochastic Environmental Research and Risk Assessment, 2007,21(4):309- 316.

[39] Genereux D P, Leahy S, Mitasova H, et al. Spatial and temporal variability of streambed hydraulic conductivity in West Bear Creek, North Carolina, USA [J]. Journal of Hydrology, 2008, 358(3):332-353.

[40] Webb R H, Leake S A. Ground-water surface-water interactions and long-term change in riverine riparian vegetation in the southwestern United States [J]. Journal of Hydrology, 2006, 320(3):302-323.

[41] Wang Genxu, Yang Lingyuan, Chen Ling, et al. Impacts of land use changes on groundwater resources in the Heihe River Basin [J]. Journal of Geographical Sciences, 2005,15(4):405-414.

[42] Chen X, Shu L. Stream-Aquifer Interactions: Evaluation of Depletion Volume and Residual Effects from Ground Water Pumping [J]. Ground Water, 2002,40(3):284-290.

[43] 朱金峰,王忠靜,鄭 航,等.黑河流域中下游全境地表-地下水耦合模型與應(yīng)用 [J]. 中國環(huán)境科學(xué), 2015,35(9):2820–2826.

[44] Christensen S, Rasmussen K R, Moller K. Prediction of regional ground water flow to streams [J]. Ground Water, 1998,36(2):351- 360.

[45] Newman B D, Vivoni E R, Groffman A R. Surface water–groundwater interactions in semiarid drainages of the American southwest [J]. Hydrological Processes, 2006,20(15): 3371-3394.

[46] 王才川,王文科,張 淵,等.河流—地下水水流模型研究進(jìn)展[J]. 地下水, 2010,(6):4-7.

[47] Srivastava K, Serrano S E, Workman S R. Stochastic modeling of transient stream–aquifer interaction with the nonlinear Boussinesq equation [J]. Journal of Hydrology, 2006,328(3): 538-547.

[48] Hatch C E, Fisher A T, Ruehl C R, et al. Spatial and temporal variations in streambed hydraulic conductivity quantified with time-series thermal methods [J]. Journal of Hydrology, 2010, 389(3):276-288.

[49] Conant B. Delineating and quantifying ground water discharge zones using streambed temperatures [J]. Groundwater, 2004, 42(2):243-257.

[50] Hatch C E, Fisher A T, Revenaugh J S, et al. Quantifying surface water–groundwater interactions using time series analysis of streambed thermal records: Method development [J]. Water Resources Research, 2006,42(10).

[51] Schmidt C, Bayer-Raich M, Schirmer M. Characterization of spatial heterogeneity of groundwater-stream water interactions using multiple depth streambed temperature measurements at the reach scale [J]. Hydrology and Earth System Sciences Discussions, 2006,3(4):1419-1446.

[52] 張 弛.埃及西奈地區(qū)地下水污染的水化學(xué)模擬及評估[J]. 中國環(huán)境科學(xué), 2014,34(5):1315.

[53] 聶振龍,陳宗宇,程旭學(xué),等.黑河干流淺層地下水與地表水相互轉(zhuǎn)化的水化學(xué)特征[J]. 吉林大學(xué)學(xué)報(bào)(地理科學(xué)版), 2005, 35(1):48-53.

[54] 張光輝,劉少玉,張翠云,等.黑河流域地下水循環(huán)演化規(guī)律研究[J]. 中國地質(zhì), 2004,31(3):289-293.

[55] Wang Ping, Yu Jingjie, Zhang Yichi et al. Groundwater recharge and hydrogeochemical evolution in the Ejina Basin, northwest China [J]. Journal of Hydrology, 2013,476:72-86.

[56] 張光輝,聶振龍,張翠云,等.黑河流域走廊平原地下水補(bǔ)給變異特征與機(jī)制[J]. 水利學(xué)報(bào), 2005,36(6):715-720.

[57] Moon S K, Woo N C, Lee K S. Statistical analysis of hydrographs and water-table fluctuation to estimate groundwater recharge [J]. Journal of Hydrology, 2004,292(1):198-209.

[58] 胡興林,蘭永超,王 靜,等.黑河中游盆地水資源轉(zhuǎn)化規(guī)律研究[J]. 地下水, 2008,30(2):23-26.

[59] 胡廣錄,趙文智.恢復(fù)生態(tài)地下水位的需水量及恢復(fù)方案研究[J]. 干旱區(qū)研究, 2009,26(1):94-101.

[60] 劉少玉,盧耀如,程旭學(xué),等.黑河中下游盆地地下水系統(tǒng)與水資源開發(fā)的資源環(huán)境效應(yīng)[J]. 地理學(xué)與國土研究, 2002, 18(4):90-96.

[61] Workman S R, Serrano S E, Liberty K. Development and application of an analytical model of stream/aquifer interaction [J]. Journal of Hydrology, 1997,200(1):149-163.

[62] Serrano S E, Workman S R, Srivastava K, et al. Models of nonlinear stream aquifer transients [J]. Journal of Hydrology, 2007,336(1):199-205.

[63] Moench A F, Barlow P M. Aquifer response to stream-stage and recharge variations. I. Analytical step-response functions [J]. Journal of Hydrology, 2000,230(3):192-210.

[64] Wurster F C, Cooper D J, Sanford W E. Stream/aquifer interactions at Great Sand Dunes National Monument, Colorado: influences on interdunal wetland disappearance [J]. Journal of Hydrology, 2003,271(1):77-100.

[65] LaBolle E M, Ahmed A A, Fogg G E. Review of the Integrated Groundwater and Surface-Water Model (IGSM) [J]. Ground Water, 2003,41(2):238-246.

[66] Langevin C, Swain E, Wolfert M. Simulation of integrated surface-water/ground-water flow and salinity for a coastal wetland and adjacent estuary [J]. Journal of Hydrology, 2005, 314(1):212-234.

[67] Sarris T S, Paleologos E K. Numerical investigation of the anisotropic hydraulic conductivity behavior in heterogeneous porous media [J]. Stochastic Environmental Research and Risk Assessment, 2004,18(3):188-197.

[68] Spanoudaki K, Stamou A I, Nanou-Giannarou A. Development and verification of a 3-D integrated surface water–groundwater model [J]. Journal of Hydrology, 2009,375(3):410-427.

[69] 王旭升,周 劍.黑河流域地下水流數(shù)值模擬的研究進(jìn)展[J]. 工程勘察, 2009,(9):35-38.

[70] 王 蕊,王中根,夏 軍.地表水和地下水耦合模型研究進(jìn)展[J]. 地理科學(xué)進(jìn)展, 2008,27(4):37-41.

[71] Gunduz O, Aral M M. River networks and groundwater flow: a simultaneous solution of a coupled system [J]. Journal of Hydrology, 2005,301(1):216-234.

[72] McDonald M G, Harbaugh A W. The history of MODFLOW [J]. Ground Water, 2003,41(2):280-283.

[73] 初京剛,張 弛,周惠成. SWAT與MODFLOW模型耦合的接口及框架結(jié)構(gòu)研究及應(yīng)用[J]. 地理科學(xué)進(jìn)展, 2011,(3):335-342.

[74] Frei S, Fleckenstein J H, Kollet S J, et al. Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model [J]. Journal of Hydrology, 2009,375(3):383- 393.

[75] 張浩佳,吳劍鋒,林 錦,等.基于GSFLOW的地下水-地表水耦合模擬與分析[J]. 工程勘察, 2015,(5):34-38.

[76] 武 強(qiáng),孔慶友,張自忠,等.地表河網(wǎng)-地下水流系統(tǒng)耦合模擬Ⅰ:模型[J]. 水利學(xué)報(bào), 2005,(5):588-592+597.

[77] 胡立堂,王忠靜,趙建世,等.地表水和地下水相互作用及集成模型研究[J]. 水利學(xué)報(bào), 2007,38(1):54-59.

[78] Panday S, Huyakorn P S. A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow [J]. Advances in Water Resources, 2004,27(4):361-382.

[79] Sulis M, Meyerhoff S B, Paniconi C, et al. A comparison of two physics-based numerical models for simulating surface water–groundwater interactions [J]. Advances in water resources, 2010,33(4):456-467.

[80] 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.

[81] 凌敏華,陳 喜,程勤波,等.地表水文過程與地下水動力過程耦合模擬及應(yīng)用[J]. 水文, 2011,6:8-13.

[82] 束龍倉,魯程鵬,李 偉.考慮參數(shù)不確定性的地表水與地下水交換量的計(jì)算方法[J]. 水文地質(zhì)工程地質(zhì), 2008,(5):68-71.

Review on the research of surface water and groundwater interactions.

ZHU Jin-feng1,2, LIU Yue-yi2,3, ZHANG Shu-an1, ZHENG Hang2,4*

(1.Bureau of Hydrology, Ministry of Water Resources, Beijing 100053, China；2.State Key Laboratory of Hydroscience and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China；3.Eastern Route Corporation of South-to-North Water Diversion, Beijing 100088, China；4.Dongguan University of Technology, Dongguan 523015)., 2017,37(8)：3002~3010

The surface water and groundwater transformation is a critical issue within hydrology circle, and their interactions are focused and difficult in the fields of hydrology, geology and environment. To analyze the exchange features between surface water and groundwater and improve researches on water resources analysis and evaluation, material movement, and energy distribution, the researches on surface water-groundwater cycle and their interactions in recent 20years were systematically analyzed and investigated. By summarizing the related research fields and core problems under natural factors and human activities, it concluded the commonly used research methods and modeling technologies. Through the comparison of related studies, it showed that comprehensive study methods (different experiment methods, analytical method and numerical method) combined with new technologies were paid more attention, and surface water and groundwater interactions research coupling natural factors and human activities became frontiers and major focuses.

surface water and groundwater；water cycle；modeling technology；water resources management

X52

A

1000-6923(2017)08-3002-09

朱金峰(1987-),男,山東濰坊人,博士,主要從事水文水資源研究.發(fā)表論文7篇.

2017-01-10

國家自然科學(xué)基金資助項(xiàng)目(51479089,51323014)

* 責(zé)任作者, 講師, zhenghang@tsinghua.edu.cn