劉 莉， 張永浩， 劉 艷， 石建芳

(西北核技術(shù)研究所，陜西西安 710024)

淺談非飽和帶補(bǔ)給測量方法

劉 莉， 張永浩， 劉 艷， 石建芳

(西北核技術(shù)研究所，陜西西安 710024)

非飽和帶補(bǔ)給測量一直是水文地質(zhì)研究的一個難點(diǎn)問題。近年來，由于對非飽和帶污染物的運(yùn)移以及干旱半干旱地區(qū)環(huán)境問題關(guān)注度的提升，使得其測量技術(shù)有了很大發(fā)展。介紹了物理方法、示蹤劑方法和數(shù)值模擬方法的應(yīng)用條件以及選擇原則。明確補(bǔ)給研究背景，全面了解場地特點(diǎn)。并要考慮測量方法的適用范圍、測量周期及其它一些因素。同時由于各種方法均具有局限性，因此在實踐中要綜合應(yīng)用多種方法，以增加計算結(jié)果的可靠性。

水文地質(zhì);非飽和帶;補(bǔ)給;測量方法

劉莉，張永浩，劉艷，等.淺談非飽和帶補(bǔ)給測量方法［J］.東華理工大學(xué)學(xué)報：自然科學(xué)版，35(3)：276-280.

Liu Li，Zhang Yong-hao，Liu Yan，et al.Preliminary discussion on the recharge estimated methods in the unsaturated zone［J］.Journal of East China Institute of Technology(Natural Science)，35(3)：276-280.

補(bǔ)給測量的方法有多種，研究人員將其分為:直接方法和間接方法;物理、化學(xué)和同位素方法;含水層響應(yīng)和徑流流入/流出的方法;飽和帶非飽和帶方法;地下水、土壤水或地下水-土壤水?dāng)?shù)值模擬及場地尺度或流域范圍內(nèi)的水文平衡模擬等方法(Sophocleous，2004)。非飽和帶是介于地面與潛水面之間的復(fù)雜三相體系，是連接潛水與大氣降水的重要通道，土壤水文活動活躍。與地下水和地表水補(bǔ)給測量相比，非飽和帶受氣候和外界因素影響強(qiáng)烈，因此其補(bǔ)給測量的難度更大。綜合上述方法并結(jié)合非飽和帶補(bǔ)給研究，可將其補(bǔ)給測量方法分為物理方法、示蹤劑法和數(shù)值模擬方法(Scanlon et al.，2002;汪丙國，2008;Allison et al.，1994)。

近年來，由于人們對非飽和帶污染物的運(yùn)移以及廢物處置問題關(guān)注度的提升，使得其測量技術(shù)有了很大發(fā)展。目前，非飽和帶補(bǔ)給測量技術(shù)主要應(yīng)用于干旱半干旱地區(qū)。如何選擇合適的測量方法，需要考慮多種因素，透徹理解各種方法的原理和適用范圍。

1 物理方法

1.1 水均衡法

水均衡法是根據(jù)質(zhì)量守恒定律，對土壤層(或非飽和帶)寫出水量平衡方程，只要弄清計算區(qū)域邊界范圍內(nèi)外的水分交換量，取得水量平衡分量的測定值，就可以得到所求補(bǔ)給量的值，Scanlon等(1997)列出了水均衡法計算公式，Kumar(2003)給出了其實施詳細(xì)步驟。該方法的優(yōu)點(diǎn)是限制條件少，適用范圍大，但其精度取決于各平衡分量的測定，其中有效降水量、土壤水蓄水量和蒸發(fā)量較難確定，特別是在干旱地區(qū)潛在蒸發(fā)量往往等于甚至遠(yuǎn)遠(yuǎn)超過降雨量(Allison et al.，1994)。

式中R為補(bǔ)給量(mm/a);P為降水量(mm/a);ET為蒸散發(fā)量(mm/a);R0為地表徑流(mm/a);ΔS為土壤蓄水量的變化(mm/a)。

1.2 蒸滲儀法

蒸滲儀是一種埋設(shè)在田間或溫室內(nèi)地下裝滿土壤的大型儀器，其土壤表面或裸露或種植各種作物，在限定三維邊界條件下測定水體轉(zhuǎn)化的裝置。根據(jù)水均衡原理對土壤水收支的每一分量精確測量，以獲得裸土或作物的騰發(fā)量、潛在騰發(fā)量以及深層滲漏量。一般分為非稱重式和稱重式土壤蒸滲儀兩類，非稱重式蒸滲儀可用于測量降雨和蓄水量，稱重式土壤蒸滲儀則可精確測量蒸散發(fā)(姜峻等，2008;陳建耀等，1999;陳建峰，2009)。隨著電子技術(shù)、核技術(shù)、土壤物理學(xué)等領(lǐng)域的進(jìn)展，蒸滲儀在形式、結(jié)構(gòu)、精度和功能方面不斷改進(jìn)，其精度高，不僅用于農(nóng)田水利，而在環(huán)境工程等領(lǐng)域得到新的應(yīng)用。使用蒸滲儀的缺點(diǎn)在于其設(shè)備安裝復(fù)雜、成本高、對環(huán)境有擾動，只能限于監(jiān)測期內(nèi)的補(bǔ)給量計算。

1.3 零通量面法

零通量面是因為非飽和帶中往往存在著水分上行運(yùn)動和下行運(yùn)動的界面(蒸散發(fā)和排泄)，其垂直水頭壓力梯度為零，以零通量面稱之。零通量面土壤水收支變化可采用負(fù)壓計、中子水份儀及TDR等儀器測量或計算(何其華等，2003;王永義，2006)，適合全年土壤含水率變化較大、地下水埋深遠(yuǎn)大于零通量面的地區(qū)。其局限性在于降雨量比較大或氣候比較寒冷的地區(qū)，零通量面存在的時間較短(葛帆等，2004)，只能在一年中某一時間內(nèi)使用;另外不能用于土壤水運(yùn)動貫穿整個剖面的情況，因為下行水濕潤鋒會模糊零通量面。與蒸滲儀一樣，零通量面法也只能計算監(jiān)測期內(nèi)的補(bǔ)給。楊建鋒等(1997)采用零通量面法在沈陽地區(qū)測定了地下水入滲補(bǔ)給量;范琦等(2006)在河北欒城地區(qū)應(yīng)用零通量面法進(jìn)行包氣帶水分運(yùn)移監(jiān)測，并計算了地下水補(bǔ)給。

1.4 達(dá)西定律法

非飽和帶補(bǔ)給也可用達(dá)西定律來計算，該法要求計算垂直水力梯度和非飽和水力傳導(dǎo)率，計算公式見式2。該方法的顯著特點(diǎn)是可以全年應(yīng)用，能在較長時間范圍內(nèi)對單點(diǎn)補(bǔ)給進(jìn)行測量;當(dāng)在深滲流區(qū)應(yīng)用時，其測量值代表的區(qū)域范圍更大。采用達(dá)西定律和飽和—非飽和帶模型模擬時，都會因為水力傳導(dǎo)率測量準(zhǔn)確性、數(shù)據(jù)適用性(如實驗室數(shù)據(jù)與野外數(shù)據(jù))及其空間變異性導(dǎo)致誤差(Weinzettel et al.，2002)，并且當(dāng)水力傳導(dǎo)率與體積含水率非線性相關(guān)時，會導(dǎo)致誤差放大，這些不確定因素綜合起來會導(dǎo)致測量結(jié)果發(fā)生量級變化。李明香等(2000)在中輻院山西太原試驗場進(jìn)行了黃土包氣帶水分運(yùn)移研究，應(yīng)用達(dá)西定律計算出包氣帶補(bǔ)給;Palolo等(2002)在阿根廷某平原地區(qū)應(yīng)用達(dá)西定律和水均衡法測定了包氣帶入滲補(bǔ)給，兩種方法測得的結(jié)果相近。

式中 K(θ)為水力傳導(dǎo)率(cm/s);H為總水頭(cm);h為基巖壓力水頭(cm);z為位置水頭(cm)。

2 示蹤劑方法

目前常用兩類示蹤劑:一是人工示蹤劑，包括溴化物、氚和染色劑等;二是環(huán)境示蹤劑，包括天然環(huán)境中存在的或者由過去人類活動、事件產(chǎn)生的物質(zhì)，如氫氧同位素、大氣層核試驗產(chǎn)生的T/36Cl、工業(yè)或農(nóng)業(yè)污染源如溴化物、硝酸鹽等。這些標(biāo)志物的出現(xiàn)，表明該處在某時間段內(nèi)接受了補(bǔ)給，為計算近期補(bǔ)給提供了定量依據(jù)。通常采用的方法是開挖探槽或剖面對示蹤劑的遷移進(jìn)行觀察和取樣，分三類計算方法:⑴根據(jù)示蹤劑峰值的位置計算;⑵根據(jù)土壤中示蹤劑剖面的形態(tài)來計算;⑶根據(jù)示蹤劑輸入總量及剖面累積量來計算(Sophocleous，2004)。如Norris等(1987)在美國Yucca山附近應(yīng)用核爆產(chǎn)生的36Cl峰值計算包氣帶大氣降水入滲速率;張之淦等(1990)在山西應(yīng)用環(huán)境氚峰值位置研究黃土包氣帶水入滲補(bǔ)給量;王福剛等(2007)在新鄉(xiāng)應(yīng)用D、18O同位素峰值位移法計算了通過包氣帶的大氣降水入滲補(bǔ)給量。近年來出現(xiàn)了氯離子示蹤法，即根據(jù)降水輸入氯離子濃度與土壤水中氯離子濃度之比值反求降水的入滲補(bǔ)給量，在實際應(yīng)用中又分為氯離子質(zhì)量平衡法和氯離子累積法(Sophocleous，2004;陳植華等，1996)。如 Scanlon(1991)應(yīng)用氯質(zhì)量平衡法測定了某廢物處置場包氣帶水通量;馬金珠等(2004)在巴丹吉林沙漠采用氯質(zhì)量平衡法結(jié)合氫氧同位素法測定了地下水多年平均補(bǔ)給率;聶振龍等(2011)利用氯質(zhì)量平衡法計算了張掖盆地降水入滲速率。

在補(bǔ)給量較小的地區(qū)，由于缺少其它合適的方法，多采用示蹤劑法(萬軍偉，2003;汪丙國等，2006)。示蹤劑方法的計算結(jié)果及解譯受其入滲機(jī)理的影響，比如示蹤劑的二次輸入及混合問題、優(yōu)先流問題等，導(dǎo)致計算結(jié)果偏高或者偏低，甚至無法應(yīng)用(Scanlon，2000)。并且由于示蹤劑是間接計算補(bǔ)給，不能提供相關(guān)補(bǔ)給細(xì)節(jié)。

3 數(shù)值模擬方法

數(shù)值模擬法就是在計算機(jī)上采用離散化的方法求解數(shù)學(xué)模型，它的解是數(shù)學(xué)模型的近似解。由于這種方法能較好的反映復(fù)雜地質(zhì)條件下地下水流的狀態(tài)，具有較高的仿真度，因此發(fā)展較快(郭曉東等，2010)。近年來由于計算機(jī)技術(shù)的發(fā)展，使得模擬長期補(bǔ)給成為可能，對于非飽和多孔介質(zhì)中水流運(yùn)動，主要以Richards方程為基礎(chǔ)(楊建鋒等，2005)，廣泛采用的程序有 HYDRUS-1D，HYDRUS-2D，VS2DT，UNSAT-H，F(xiàn)EFLOW 軟件等(表 1)(Fayer，2000;Gassman et al.，2007;王旭東等，2004)，但同時還須結(jié)合實驗研究和具體地質(zhì)條件的研究(芮孝芳等，2004;薛禹群，2010)，比如結(jié)合蒸滲儀、示蹤劑、含水率和溫度等數(shù)據(jù)(Scanlon，1997)。數(shù)值模擬還可用來進(jìn)行水流過程以及參數(shù)敏感性分析。其缺點(diǎn)在于模擬所需數(shù)據(jù)量大，并且結(jié)果受各輸入?yún)?shù)的影響。數(shù)值模擬方法對時間尺度沒有限制，但是氣候數(shù)據(jù)的時效性將其測量時間限制在 100 a左右(Scanlon，1997)。韓錦萍(2005)、范磊(2008)采用HYDRUS-1D程序?qū)Χ鯛柖嗨古璧氐陌鼩鈳诌\(yùn)移及降水入滲補(bǔ)給進(jìn)行了數(shù)值模擬，作為實驗結(jié)果的驗證手段。

表1 部分?jǐn)?shù)值模擬軟件簡介Table 1 Introductionof the numerical simulation software

4 測量方法應(yīng)用原則

毫無疑問，補(bǔ)給測量方法的選擇與環(huán)境因素緊密相關(guān)，包括地理環(huán)境、氣候、土壤地質(zhì)及水文地質(zhì)條件、植被、土地使用情況等。除此之外，主要關(guān)注以下三個方面:

⑴補(bǔ)給研究背景。包括研究目的、補(bǔ)給源、補(bǔ)給機(jī)理、補(bǔ)給量及地下水埋深等，并通過對具有相似氣候、地理環(huán)境的場地的對比分析，預(yù)測補(bǔ)給量，以選擇恰當(dāng)?shù)姆椒?。如在水資源評價對非飽和帶水分遷移轉(zhuǎn)化規(guī)律研究中，其時空尺度通常較大;示蹤劑方法有時是唯一能提供完整、長期補(bǔ)給信息的方法，例如采用氯離子示蹤法;較厚的非飽和帶會導(dǎo)致某些示蹤劑如CFCs濃度及比值的改變(萬軍偉等，2003);在應(yīng)用環(huán)境示蹤劑時，要求補(bǔ)給量能使示蹤劑既能穿過根系層但又不會遷移太深，反映了非飽和帶厚度和示蹤劑峰值定位對該方法的限制。

⑵測量方法適用范圍。包括計算方法所適用的時空尺度、測量精度、可靠性等，圖1為各補(bǔ)給測量方法綜合比較(Scanlon et al.，1997，2006)。

⑶其它因素。主要是測量時間及費(fèi)用方面的考慮，如水均衡法和蒸滲儀需要長期監(jiān)測，而其它方法所需時間則較短;就研究成本而言，化學(xué)和同位素示蹤劑方法的采樣和分析都較為復(fù)雜繁瑣，甚至需要質(zhì)譜儀和液閃儀等高新技術(shù)和設(shè)備(萬軍偉等，2003)，但采樣通常為一次性，與長期持續(xù)監(jiān)測所需要的設(shè)備以及數(shù)據(jù)收集、分析測定所耗費(fèi)的資源相比，其效率也許更高。

5 結(jié)語

補(bǔ)給測量首先需要全面了解場地特點(diǎn)，建立正確的概念模型，針對場地特點(diǎn)和研究目的選擇恰當(dāng)?shù)姆椒?。水均衡法較為適合濕潤地區(qū)，該法明顯受植被和土地使用情況影響;蒸滲儀法的監(jiān)測周期長、費(fèi)用高，對環(huán)境有擾動;零通量面法則與補(bǔ)給量大小及零通量面的存在相關(guān);達(dá)西定律法在干旱地區(qū)應(yīng)用時，土壤含水率及非飽和傳導(dǎo)率對其影響較大;示蹤劑法技術(shù)應(yīng)用較廣，是干旱地區(qū)其它方法難以應(yīng)用時的最佳選擇;數(shù)值模擬方法則受概念模型及各輸入?yún)?shù)精度的限制?？傮w來說，非飽和帶補(bǔ)給通常是對單點(diǎn)或者局部區(qū)域進(jìn)行測量，但如果將區(qū)域內(nèi)多點(diǎn)地形數(shù)據(jù)及其補(bǔ)給計算結(jié)果結(jié)合起來，就可推及較大區(qū)域;由于每種技術(shù)方法都有局限性，所以要綜合應(yīng)用多種方法，反復(fù)迭代、逐步優(yōu)化和相互補(bǔ)充。

陳建峰.2009.淺析包氣帶土壤水變化特征［J］.地下水，31(1):53-54.

圖1 補(bǔ)給測量方法比較Fig.1 Comparison of the recharge estimate methods(Modified from Scanlon et al.)

陳建耀，劉昌明.1999.利用大型蒸滲儀模擬土壤、植物、大氣連續(xù)體水分蒸散［J］.應(yīng)用生態(tài)學(xué)報，10(1):45-48.

陳植華，徐恒力.1996.確定干旱-半干旱地區(qū)降水入滲補(bǔ)給量的新方法——氯離子示蹤法［J］.地質(zhì)科技情報，15(3):87-91.

范磊.2008.降水入滲過程中包氣帶水運(yùn)移轉(zhuǎn)化機(jī)理研究［D］.西安:長安大學(xué)圖書館.

范琦，王驥等.2006.包氣帶增厚條件下地下水補(bǔ)給規(guī)律研究［J］.水文地質(zhì)工程地質(zhì)，(3):21-24.

葛帆，王釗.2004.蒸滲儀及其應(yīng)用現(xiàn)狀［J］.節(jié)水灌溉，(2):30-35.

郭曉東，田輝，張梅桂，等.2010.我國地下水?dāng)?shù)值模擬軟件應(yīng)用進(jìn)展［J］.地下水，32(7):5-7.

韓錦萍.2005.包氣帶水分運(yùn)移的數(shù)值模擬分析［D］.西安:長安大學(xué)圖書館.

何其華，何永華，包維楷.2003.干旱半干旱區(qū)山地土壤水分動態(tài)變化［J］.山地學(xué)報，21(2):149-156.

姜峻，都全勝，趙軍，等.2008.稱重式蒸滲儀系統(tǒng)改進(jìn)及在農(nóng)田蒸散研究中的應(yīng)用［J］.水土保持通報，28(6):67-72.

李明香，馬炳輝.2000.黃土包氣帶水分運(yùn)移的現(xiàn)場研究［J］.輻射防護(hù)，20(1-2):91-100.

馬金珠，李丁，李相虎，等.2004.巴丹吉林沙漠包氣帶Cl－示蹤與氣候記錄研究［J］.中國沙漠，24(6):674-679.

聶振龍，連英立.2011.利用包氣帶環(huán)境示蹤劑評估張掖盆地降水入滲速率［J］.地球?qū)W報，32(1):117-122.

芮孝芳，黃國如.2004.分布式水文模型的現(xiàn)狀與未來［J］.水利水電科技進(jìn)展，(2):55-58.

萬軍偉.2003.同位素水文學(xué)理論與實踐［M］.北京:中國地質(zhì)大學(xué)出版社.

汪丙國，靳孟貴，王文峰，等.2006.氯離子示蹤法在河北平原地下水垂向入滲補(bǔ)給量評價中的應(yīng)用［J］.節(jié)水灌溉，(3):16-20.

汪丙國.2008.地下水補(bǔ)給評價方法研究——以華北平原為例［D］.武漢:中國地質(zhì)大學(xué):10-18.

王福剛，廖資生.2007.應(yīng)用D、18O同位素峰值位移法求解大氣降水入滲補(bǔ)給量［J］.吉林大學(xué)學(xué)報，37(2):284-334.

王旭東，蔣云鐘，趙紅莉，等.2004.分布式水文模擬模型在流域水資源管理中的應(yīng)用［J］.南水北調(diào)與水利科技，(1):21-25.

王永義.2006.非飽和帶水研究的重要性及其研究方法［J］.西北水利發(fā)電，22(5):28-30.

薛禹群.2010.中國地下水?dāng)?shù)值模擬的現(xiàn)狀與展望［J］.高校地質(zhì)學(xué)報，16(1):1-6.

楊建鋒，李同斌，等.1997.零通量面法在沈陽地區(qū)地下水資源評價中的應(yīng)用［J］.世界地質(zhì)，16(3):55-60.

楊建鋒，萬書勤，鄧偉，等.2005.地下水淺埋條件下包氣帶水和溶質(zhì)運(yùn)移數(shù)值模擬研究述評［J］.農(nóng)業(yè)工程學(xué)報，21(6):158-165.

張之淦，劉芳珍.1990.應(yīng)用環(huán)境氚研究黃土包氣帶水分運(yùn)移及入滲補(bǔ)給量［J］.水文地質(zhì)工程地質(zhì)，(3):5-7.

Allison G B，Gee et G W.1994.Vadose-Zone Techniques for estimating groundwater recharge in arid and semiarid regions［J］.Soil Science.Soc.AM，58:6-14.

Fayer M J.2000.Unsaturated soil water and heat flow model theory，user manual，and examples［OL］. http://221.7.11.12:81/1Q2W3E4R5T6Y7U8I9O0P1Z2X3C4V5B/www.pnl.gov/main/publications/external/technical_reports/PNNL-13249.pdf.

Gassman P W，Reyes M R，et al.2007.The soil and water assessment tool:historical development applications，and future research directions［J］.Transactions of the ASABE，50(4):1211-1250.

Kumar C P.2003.Estimation of groundwater recharge using soil moisture balance approach ［J］.Journal of soil and water conservation，Soil Conservation Society of India，2(1):53-58.

Norris A E，Wolfsberg K.1987.Infiltration at Yucca Mountain，Nevada traced by Cl-36 ［J］.Nuclear instruments and methods in physics research B29:376-379.

Palolo A W，Eduardo J U et al.2002.Groundwater recharge estimations from studies of the unsatured zone［J］.Groundwater and human development:731-738.

Scanlon B R，Healy R W.2002.Choosing appropriate techniques for quantifying groundwater recharge［J］.Hydrogeology Journal，(10):18-39.

Scanlon B R，Keese K E，et al.2006.Global synthesis of groundwater recharge in semiarid and arid regions［J］.Hydrological processes，(20):3335-3370.

Scanlon B R，Scott W T.1997.Hydrologic issues in arid，unsaturated systems and implications for contaminant transport［J］.Review of Geophysics，35(4):461-490.

Scanlon B R.1991.Evaluation of moisture flux from chloride data in desert soils［J］.Hydrol，128:137-156.

Scanlon B R.2000.Uncertainties in estimating water fluxes and residence times using environmental tracers in an arid unsaturated zone［J］.Water resources research，36(2):395-409.

Sophocleous M.2004.Groundwater recharge［EB/OL］.http://www.eolss.net.

Weinzettel P A，Usunoff E J.2002.Groundwater recharge estimations from studies of the unsaturated zone［J］.Groundwater and human development:731-738.

Preliminary Discussion on the Recharge Estimated Methods in the Unsaturated Zone

LIU Li， ZHANG Yong-hao， LIU Yan， SHI Jian-fang
(Northwest Institute of Nuclear Technology，Xi'an，SX 710024，China)

To estimate the recharge in the unsaturated zone is a difficult problem in the hydrogeological fields.In recent years，the estimate methods develop greatly because the problems of contaminant migration and waste disposal in unsaturated zone have attracted more attention.The physical，tracers and numerical methods respectively are introduced，and some advices on choosing the appreciated methods are given.First of all，the background of the study is important.The applicability and reliability of the methods must be considered as well as some other factors.For each method has its limitations，multiple techniques should be used to increase reliability of results in practice.

hydrogeology;unsaturated zone;recharge;estimate methods

P64.131

A

1674-3504(2012)03-0276-05

10.3969/j.issn.1674-3504.2012.03.012

2011-12-16 責(zé)任編輯:吳志猛

劉 莉(1973—)，女，工程師，主要從事環(huán)境地質(zhì)研究。E-mail:liuliml@sohu.com