地下水與地表水相互作用的識別和量化方法研究

摘要結合國內外研究現狀，將地下水與地表水相互作用的識別和量化方法歸納總結為直接測量法、調查分析法、模擬計算法和環境示蹤法4類，并對各類方法的優缺點和適用范圍進行了比較，最后對其未來發展趨勢作出了展望。

關鍵詞地下水與地表水相互作用；滲透儀；流量過程線；耦合模型；環境示蹤劑

中圖分類號S181.3文獻標識碼A文章編號0517-6611（2015）24-203-03

在自然界中，地下水與地表水是一個不可分割的整體，其中任何一方的開發和污染都將會給另一方造成不可忽視的影響[1]。例如，地下水的過度開采可能造成地表徑流的衰減甚至干涸，地下水中營養物質的排泄會導致地表水體富營養化，海水入侵引起地下水源鹽度升高。如果水資源的管理政策中不考慮地下水與地表水的相互作用，水資源匱乏、水質惡化和生態環境破壞等一系列水問題將得不到有效的解決[2]。對此，聯合國教科文組織（UNESCO）、國際水文科學協會（IAHS）和國際水文地質學家協會（IAH）等組織和機構，都把地下水與地表水的相互作用作為目前研究的重要熱點和前沿課題[3]，以期為準確評價水資源，合理規劃管理水資源，有效保護水資源和生態環境系統提供技術支撐。

關于地下水與地表水的相互作用研究，早在1877年，Boussinesq就對河流與連續沖擊含水層作用的規律進行探討。而直到19世紀60年代，人們對水體富營養化、酸雨以及生態系統功能退化等問題的關注，才引發了地下水與湖水、河水、濕地以及海水等地表水相互作用的大量研究[4]。至今，有關兩水相互作用的研究對象幾乎涵蓋了沖積平原、基巖山區、喀斯特巖溶區以及濱海地區等所有的地下水和地表水，研究內容涉及了地下水與地表水的聯合調度規劃[5]、地下水與地表水的綜合污染防治[6]以及水生生態系統的環境健康[7]等多個方面。但是，如何準確識別和量化地下水與地表水的相互作用一直是一個非常復雜和困難的問題，它對分析流域水循環機制、準確評價水資源量和控制流域水污染具有重要的意義。筆者對地下水與地表水相互作用的識別和量化方法進行總結歸納，將目前國內外研究中所用到的方法分為直接測量法、調查分析法、模擬計算法和環境示蹤法，并指出了其未來發展的方向。

1直接測量法

地下水與地表水水量和溶質交換量的直接測量工具是滲透儀。滲透儀最早由Israelsen和Reeve于1944年設計用于測量灌溉渠中地表水的滲透損失量。之后，Lee[8]用一個208 L的鐵筒接一個塑料收集袋子制成了最為經典的滲透儀，即將一個開口的圓筒倒扣在河床沉積物中，用袋子收集地下進入沉積物的水，或者在袋子中裝入已知體積的水測量地表損失的水。基于此，為了避免收集袋子受水流壓迫影響，研究者們將收集的袋子改為排水管，并安裝熱脈沖裝置、超聲波裝置或者電磁表等設計出不同類型的自動滲透儀來記錄水流速度[9-11]。目前，滲透儀這種直接測量法已被廣泛應用于湖、河、水庫、海洋以及河流與地下水的定量研究中[12-13]。

盡管這種用來直接量化地下水和地表水交換水量的滲透儀原理簡單，經濟實用，但是其操作非常麻煩。一般來說，為了獲得有代表性的平均滲透通量，往往需要進行多個點的測量，比較耗時費力，且在測量時間段內滲透通量的變化過程也無法掌握。特別是通過滲透儀測量的滲透通量可能不完全是地下水排泄造成的，還可能包括河流沉積物中排泄的水等[14]。因此，滲透儀的測試結果可能存在一定的誤差。

2調查分析法

調查分析包括長期的流量資料調查分析（即流量過程線分析）和某一均衡期的河流水均衡項調查計算。

2.1流量過程線分析

河流的流量過程曲線是指河流流量隨時間變化的曲線，該流量包括河流的降水徑流量和基流量（一般是指地下水補給量）兩個部分。在地下水與地表水相互作用的研究中，研究者通過圖解法、滑動最小值法、HYSPEP和數字濾波法等方法將流量過程曲線中的基流量分割出來，它屬于基流分割法的一種[15]。例如，林學鈺等[16]通過流量過程線分析方法查明了黃河流域地下水對黃河的貢獻。伍立群等[17]將基流分割法應用于云南地下水資源量的評價，提高了山丘區地下水資源評價的精度。

但是，由于流量過程線分析的方法將壤中流和河岸儲量都歸為地下水對地表水的基流量[18]，而且要求長序列的河流流量數據，因此它在計算地下水對地表水的補給量中應用較少。

2.2水均衡計算

水均衡法是研究地下水與地表水相互作用較為簡單的方法。若以地表水或地下水為研究主體，將一條河流分成若干段或盡量選取一個具有隔水邊界的完整水文地質單元，查明該地區某一時間段內進入和流出的均衡項，即可計算出地下水與地表水的水量交換量。其中，涉及河流流量測量的方法有水位-流量法、體積法、速度-面積法、坡度-面積法、稀釋法、三角堰法等[19]。如肖長來等[20]以地下水為研究主體，用水均衡法研究了洮兒河扇形地地下水與地表水資源的轉化規律。

水均衡法雖然在確定地下水與地表水大致的補排關系，粗略估算兩水交互量的研究中快速方便，但因為其各均衡項存在多方面的不確定性，且容易忽略短期內地下水與地表水交互量的微小變化，會導致計算結果不可信。

3模擬計算法

模擬計算法是基于達西定律的一種常用量化方法，它是將兩水之間的界面看作地下透水層的一部分，通過野外勘察獲得地下水和地表水的水位、河床沉積物的滲透系數，用達西定律計算確定地下水排泄到地表水或地表水補給地下水的水量。

Rushton等[21]曾對地下水與地表水之間交互帶的達西定律水量計算提出改進的線性、非線性和綜合模型，認為地下水與地表水水位差較大時，達西線性模型的計算結果偏大，地下水與地表水的交換量應該與水頭差之間存在線性和非線性關系。但是，由于難以獲取參數，該結論并沒有得到廣泛應用。

迄今為止，基于達西定律的地下水與地表水耦合模型一直是地下水與地表水相互作用的重要研究方向。胡立堂等[22]將現有的水流耦合模型分為嵌入地表水模塊的地下水模型、嵌入地下水模塊的地表水模型和地下水與地表水模塊兼容模型3種類型。

嵌入地表水模塊的地下水模型和嵌入地下水模塊的地表水模型是最早用到的耦合模型，這兩類模型不考慮地下水或地表水水位的時空變化，只是將其簡單概化為邊界條件或源匯項輸入地下水或地表水模型。如Mcdonal和Harbaugh按照簡單的達西公式將地表水模塊（Stream模塊）添加到美國地質調查局開發的MOFLOW地下水模擬軟件中。席丹等[23]利用SWAT模型的水文響應單元（HRU）和MODFLOW模型的有限差分網格（CELL），構建地下水與地表水模型的交互面，將SWAT模型計算的地下水補給量和潛水蒸發量分別輸入MODFLOW模型的地下水補給（RCH）模塊和潛水蒸發（EVT）模塊中，并應用于西安市黑河流域的地下水模擬計算，取得了很好的效果。

地下水與地表水模塊兼容模型是通過交換通量將成熟的地下水方程與地表水方程聯系起來，考慮水頭隨時間的變化，它較嵌入型模型更能真實地反映實際情況，如MODBRANCH、MIKE-SHE、SWATMOD和MODHMS、HydroGeoSphere和GSFOW等耦合模擬軟件中的模型[24]。龍玉橋[25]應用HydroGeoShpere軟件，將二維圣維南方程組和三維Richard方程通過交互帶水量交換量耦合模擬計算了遼寧白石水庫的最小下泄量。

基于達西定律的方法對于認識地下水系統內部結構非常重要，它是理解地下水與地表水相互作用機制，確定兩水交互量的有效手段。但是，資料短缺情況下的模擬和地下水系統的空間非均質性處理是該方法需要解決的難點。

4環境示蹤法

天然水的化學成分在一定程度上記載著水體形成和運移歷史的信息。在地下水與地表水的相互作用中，兩水獨具特色的溶解組分在彼此的水力聯系下發生了特征性的變化，同時也留下了相互作用的印記。因此，將地下水與地表水中有明顯差異的物質作為環境示蹤劑來追蹤水流運動，是識別和量化兩水相互作用非常有效的手段之一。目前，研究中用到的示蹤劑包括水環境參數、常見溶解性組分、同位素及其他物質。

4.1水環境參數

研究發現，地下水與地表水差異明顯的水環境參數包括pH、溶解氧、氧化還原電位、電導率、堿度和溫度等。如Edet等[26]通過尼日利亞東南部Calabar地區連續12個月的地下水與河水溫度、EC、pH、Eh和DO的測定，發現該區地下水和河水的溫度、EC、Eh和DO都存在顯著性差異，其中地下水中只有DO受季節和潮汐的顯著性影響，而河水中除了DO受季節和潮汐的共同影響外，溫度和Eh值還會隨季節的變化而明顯變化。宋獻方等[27]利用地下水和地表水電導率的差別，根據流域內不同水體的電導率空間分布趨勢推斷了兩水的補排關系。Rodgers等[28]將堿度與二氧化硅作為天然示蹤劑，識別了蘇格蘭Cairngorms地區Feshie河地下水與地表水的水力聯系。

值得說明的是，從20世紀80年代，溫度示蹤劑開始被逐漸引入不同類型的水生環境研究，直到最近，地表水和沉積物溫度的監測才成為識別和量化地下水與地表水水力聯系的有效手段[29-30]。用耦合熱傳導對流方程來計算地下水與地表水之間水流速度在國外已被廣泛應用[31]。

4.2常見溶解性組分

用來示蹤地下水與地表水相互作用的常見溶解性組分包括主要離子和可溶性SiO2等。

郎赟超[32]通過比較貴陽不同季節地下水與地表水中K+、Na+、Ca2+、Mg2+、Cl-、SO2-4、HCO-3、NO-3、SiO2和TDS的變化，發現地下水對地表水的水化學組成具有明顯的響應關系，說明喀斯特地區地下水與地表水之間的水化學聯系密切，水流交換比較活躍。Moore等[33]使用多元統計的方法如主成分分析（PCA）、因子分析和聚類分析，對大量的水化學數據進行了處理，識別了控制地下水和地表水水化學成分的主要因子，指出了不同水樣點之間的關系。在這些主要離子中，Cl-作為保守性離子，可通過質量平衡方程定量評價地下水與地表水之間的混合比例[34]。

自然界水中的可溶性SiO2是巖石風化和水循環條件的有效指示劑，它可為地下水與地表水的相互作用提供可靠的依據[35]。由于降水中SiO2的含量較低，一般為1 mg/L以下，而地下水中由于周圍硅酸鹽礦物的水解會存在大量的可溶性SiO2。因此，在地下水輸入地表水時，地表水中SiO2含量則會增加；反之，如果地表水滲入地下水，則地下水中SiO2含量會減小。Gao等[36]對地下水中的可溶SiO2做了等值線，通過靠近地表水的地下水SiO2濃度高低來判別地下水與地表水之間的補給關系。

4.3同位素及其他物質

與水環境參數和常見溶解性固體一樣，不同水體的同位素組成受不同環境分餾作用的影響也有所不同。

如地表水中2H和18O因強烈的蒸發作用有明顯富集，Kumar等[37]利用這兩個穩定同位素辨識了地下水和地表水的水力聯系。由于放射性同位素222Rn從巖石或沉積物中釋放出來隨地下水流遷移排泄于地表時，因周圍分壓力的迅速減小造成地表水中222Rn濃度低于地下水。袁曉婕等[38]應用222Rn評價了膠州灣的海底地下水排泄。而鍶不發生動力學分餾，87Sr/86Sr值只隨著巖石礦物中Rb/Sr值和水巖作用的時間而發生變化。Négrel等[39]將穩定同位素δD和δ18O值與87Sr/86Sr值相結合，分析了法國盧瓦爾河附近不同層位地下水與地表水的復雜補排關系。此外，18O/16O、碳和鐳等的同位素以及人工合成化合物SF6和CFC（氟利昂）也有應用[40]。

總體來說，流域水體中的水化學成分來源多樣，發生的化學反應較為復雜，研究者們不會使用單獨一種示蹤劑來識別地下水與地表水的相互作用，而往往將各種水文地球化學成分相互比對，全面掌握地下水與地表水的水化學聯系，并由幾種敏感的環境示蹤劑相互驗證，共同確定兩水的相互轉化關系，甚至通過不同補給源的水化學成分分析，采用端元混合法（EMMA）得到地下水與地表水的混合比例。

5結論與展望

近年來，地下水與地表水相互作用的識別和量化技術手段不斷提高，但各方法都有一定的優缺點，適用的范圍不同。

①直接測量法原理簡單、成本低，但操作耗時費力，誤差較大，適用于小尺度范圍（10-2～1 m）；

②調查分析法方便、快速，但不確定因素多，難以估算短期變化量，結果可能不可信，適用于大尺度范圍（10～103 m）；

③模擬計算法具有可視化、仿真性，但耦合模擬存在難點，需要的數據量大，空間非均質性難以處理，操作復雜，適用于大、中、小尺度范圍（10-2～103 m）；

④環境示蹤法原理簡單、實用有效，但精度不高，溫度和同位素示蹤成本高，同位素采樣要求高，干擾因素多，適用于中尺度范圍（1～102 m）。

其中，溫度示蹤和模擬計算是目前地下水與地表水相互作用識別和量化方法應用的熱點方向。但是，由于地下水溫度的測量費效比（費用/效果）不高，耦合模擬的仿真性在實際情況下限制較多。因此，還需要不斷改進或開發出更為簡單有效且適用性廣的方法，來推動地下水與地表水相互作用的研究。

43卷24期鄭潔瓊地下水與地表水相互作用的識別和量化方法研究

參考文獻

[1] SOPHOCLEOUS M.Interaction between groundwater and surface water：the state of the science[J].Journal of Hydrology，2002，10：52-67.

[2] WINTER T C，HARVEY J W，FRANK O L，et al.Ground water and surface water：A single resource[M].USGS Circular 1139，1998.

[3] 胡俊鋒，王金生，滕彥國.地下水與河水相互作用的研究進展[J].水文地質工程地質，2004，31（1）：108-113.

[4] WINTER T C，ROSENBERRY D O.The interaction of ground water with prairie pothole wetlands in the Cottonwood Lake area，eastcentral North Dakota，1979-1990[J].Wetlands，1995，15（3）：193-211.

[5] 孫棟元，伊力哈木，馮省利，等.干旱內陸河流域地表水地下水聯合調度研究進展[J].地理科學進展，2009，28（2）：167-173.

[6] HILDEBRANDT A，GUILLAMN M，LACORTE S，et al.Impact of pesticides used in agriculture and vineyards to surface and groundwater quality（North Spain）[J].Water Research，2008，42（13）：3315-3326.

[7] HAYASHI M，ROSENBERRY D O.Effects of ground water exchange on the hydrology and ecology of surface water[J].Ground Water，2002，40（3）：309-316.

[8] LEE D R.Device for measuring seepage flux in lakes and estuaries[J].Limnol Oceanogr，1977，22（1）：140-147.

[9] KRUPA S L，BELANGER T V，HECK H H，et al.Krupaseep-the next generation seepage meter[J].Journal of Coastal Research，1998，26：210-213.

[10] PAULSEN R J，SMITH C F，O'ROURKE D，et al.Development and evaluation of an ultrasonic groundwater seepage meter[J].Ground Water，2001，39（6）：904-911.

[11] ROSENBERRY D O，MORIN R H.Use of an electromagnetic seepage meter to investigate temporal variability in lake seepage[J].Ground Water，2004，42（1）：68-77.

[12] LANDON M K，RUS D L，HARVEY F E.Comparison of instream methods for measuring hydraulic conductivity in sandy streambeds[J].Ground Water，2001，39（6）：870-885.

[13] 袁瑞強.黃河三角洲淺層地下水特征及其向海輸送[D].青島：中國海洋大學，2006.

[14] KALBUS E，REINSTORF F，SCHIRMER M.Measuring methods for groundwatersurface water interactions：A review[J].Hydrology and Earth System Sciences，2006，10：873-887.

[15] 陳利群，劉昌明，李發東.基流研究綜述[J].地理科學進展，2006，25（1）：1-15.

[16] 林學鈺，廖資生，錢云平，等.基流分割法在黃河流域地下水研究中的應用[J].吉林大學學報：地球科學版，2009，39（6）：959-967.

[17] 伍立群，代興蘭.河川基流分割法在山丘區地下水資源量評價中的運用[J].中國農村水利水電，2005（1）：35-38.

[18] BRODIE R S，HOSTETLER S.A review of techniques for analysing baseflow from stream hydrographs[C]//Proccedings of the NZHS-IAH-NZSSS 2005 Conference.Auckland，New Zealand，2005.

[19] HAUER F R，LAMBERTI G A.Methods in stream ecology[M].San Diego：Academic Press，1996.

[20] 肖長來，張力春，方樟，等.洮兒河扇形地地表水與地下水資源的轉化關系[J].吉林大學學報：地球科學版，2006，36（2）：234-239.

[21] RUSHTON K R，TOMLINSON L M.Possible mechanisms for leakage between aquifers and rivers[J].J Hydrol，1979，40：49-65.

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

[23] 席丹，張世艷，何慶龍，等.基于GPR技術的地表水地下水的耦合計算及應用[J].水資源研究，2014（3）：298-306.

[24] 徐力剛，張奇，左海軍.地表水地下水的交互與耦合模擬研究現狀與進展[J].水資源保護，2009，25（5）：82-102.

[25] 龍玉橋.HydroGeoSphere在白石水庫最小下泄流量研究中的應用[D].長春：吉林大學，2008.

[26] EDET A，WORDEN R H.Monitoring of the physical parameters and evaluation of the chemical composition of river and groundwater in Calabar（Southeastern Nigeria）[J].Environ Monit Assess，2009，157：243-258.

[27] 宋獻方，劉相超，夏軍，等.基于環境同位素技術的懷沙河流域地表水和地下水轉化關系研究[J].中國科學D輯：地球科學，2007，37（1）：102-110.

[28] RODGERS P，SOULSBY C，PETRY J，et al.Groundwater-surface-water interactions in a braided river：A tracerbased assessment[J].Hydrol Process，2004，18（7）：1315-1332.

[29] SILLIMAN S E，BOOTH D F.Analysis of timeseries measurements of sediment temperature for identification of gaining vs.losing portions of Juday Creek，Indiana[J].Journal of Hydrology，1993，146：131-148.

[30] 劉傳琨，胡玥，劉杰，等.基于溫度信息的地表-地下水交互機制研究進展[J].水文地質工程地質，2014，41（5）：5-10.

[31] TANIGUCHI M，BURNETT W C，SMITH C R，et al.Spatial and temporal distributions of submarine discharge rates obtained from vavious types of seepage meters at a site in the Northeastern Gulf of Mexico[J].Biogeochemistry，2003，66：35-53.

[32] 郎赟超.喀斯特地下水文系統物質循環的地球化學特征[D].貴陽：中國科學院地球化學研究所，2005.

[33] MOORE P J，MARTIN J B，SCREATON E J.Geochemical and statistical evidence of recharge，mixing，and controls of spring discharge in an eogenetic karst aquifer[J].Journal of Hydrology，2009，376：443-455.

[34] MULLANEY J R，LORENZ D L，ARNTSON A.Chloride in groundwater and surface water in areas underlain by the glacial aquifer system，Northern United States[R].USGS，2009.

[35] ASANO Y，OHTE N，UCHIDA T.Sources of weatheringderived solutes in two granitic catchments with contrasting forest growth[J].Hydrol Proc，2004，18：651-666.

[36] GAO X B，WANG Y X，WU P L，et al.Trace elements and environmental isotopes as tracers of surface watergroundwater interaction：A case study at Xin'an karst water system，Shanxi Pronvince，Northern China[J].Environ Earth Sci，2010，59：1223-1234.

[37] KUMAR U S，SHARMA S，NAVADA S V.Recent studies on surface watergroundwater relationships at hydroprojects in India using environmental isotopes[J].Hydrol Process，2008，22：4543-4553.

[38] 袁曉婕，郭占榮，馬志勇，等.基于222Rn質量平衡模型的膠州灣海底地下水排泄[J].地球學報，2015，36（2）：237-244.

[39] NGREL P H，PETELETGIRAUD E，BARBIER J.Surface watergroundwater interactions in an alluvial plain：Chemical and isotopic systematic[J].Journal of Hydrology，2003，277：248-267.

[40] COOK P G，FAVREAU G，DIGHTON J C，et al.Determining natural groundwater influx to a tropical river using radon，chlorofluorocarbons and ionic environmental tracers[J].Journal of Hydrology，2003，277：74-88.