李杰彪，蘇銳，田霄，周志超

（核工業(yè)北京地質(zhì)研究院，中核高放廢物地質(zhì)處置評(píng)價(jià)技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京 100029）

干旱-半干旱地區(qū)地下水補(bǔ)給方法研究綜述

李杰彪，蘇銳，田霄，周志超

（核工業(yè)北京地質(zhì)研究院，中核高放廢物地質(zhì)處置評(píng)價(jià)技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京 100029）

地下水補(bǔ)給是水文地質(zhì)研究中的一個(gè)重要方面?？偨Y(jié)了干旱-半干旱地區(qū)地下水補(bǔ)給研究的4種方法：①直接測(cè)定方法；②物理方法；③同位素方法；④數(shù)值模擬方法。并對(duì)每種方法的原理、研究現(xiàn)狀、優(yōu)缺點(diǎn)以及適用性等幾個(gè)方面進(jìn)行論述，旨在為水文地質(zhì)方法學(xué)研究提供參考。

地下水補(bǔ)給；水文地質(zhì)；干旱-半干旱；方法學(xué)

地下水已成為制約經(jīng)濟(jì)發(fā)展的重要因素之一，在干旱-半干旱地區(qū)顯得尤為突出［1］。一個(gè)地區(qū)地下水的合理開發(fā)和科學(xué)管理首先要評(píng)價(jià)地下水的補(bǔ)給條件?！暗叵滤a(bǔ)給是指含水層或含水系統(tǒng)從外界獲得水量的過程。補(bǔ)給的研究包括：補(bǔ)給來源、補(bǔ)給條件以及補(bǔ)給量［2］”。地下水補(bǔ)給是地下水資源評(píng)價(jià)的重要組成部分，是研究水文循環(huán)過程的基礎(chǔ)，是分析區(qū)域地下水分布規(guī)律、進(jìn)行地下水流數(shù)值模擬不可缺少的參數(shù)。因此，在高放廢物地質(zhì)處置庫場址預(yù)選和評(píng)價(jià)中地下水補(bǔ)給研究是一個(gè)重要內(nèi)容。

1 干旱-半干旱地區(qū)地下水補(bǔ)給研究方法介紹

干旱-半干旱地區(qū)地下水補(bǔ)給的研究對(duì)象為地表水、包氣帶和飽水帶。目前，國內(nèi)、外干旱-半干旱地區(qū)地下水補(bǔ)給研究方法主要有：①直接測(cè)定方法；②物理方法（零通量面方法、達(dá)西方法和地下水位動(dòng)態(tài)方法）；③同位素方法（歷史示蹤劑、環(huán)境示蹤劑和人工示蹤劑）；④數(shù)值模擬方法。干旱-半干旱地區(qū)地下水補(bǔ)給研究方法見表1［3］。

1.1 直接測(cè)定方法

直接測(cè)定方法又稱為地中滲透儀方法，就是用埋在土中的地中滲透儀直接測(cè)定，通過測(cè)量土柱下滲量確定地下水的補(bǔ)給量。

一般認(rèn)為，最早的地中滲透儀是1976年在英國和瑞士同時(shí)開始應(yīng)用的。國外對(duì)地中滲透儀的建造及其在地下水補(bǔ)給中的應(yīng)用都做了深入的研究。Aboukhaled（1981）對(duì)不同類型地中滲透儀的設(shè)計(jì)、建造以及優(yōu)、缺點(diǎn)等進(jìn)行了系統(tǒng)的總結(jié)［4］。Gee等（1992）收集了美國漢福德市20多年的地中滲透儀觀測(cè)數(shù)據(jù)來研究該區(qū)地下水補(bǔ)給的空間變異性。我國這方面起步較晚，自20世紀(jì)60年代才開始運(yùn)用地中滲透儀開展地下水入滲補(bǔ)給和蒸發(fā)方面的研究。吳金泉等（1994）在河北某試驗(yàn)場應(yīng)用地中滲透儀對(duì)大氣降水入滲補(bǔ)給機(jī)理進(jìn)行了分析［5］。郭會(huì)榮等（2007）通過在某均衡試驗(yàn)場地中的滲透儀上開展土壤水流穿透試驗(yàn)，分析灌溉水對(duì)地下水的補(bǔ)給方式。

直接測(cè)定方法的優(yōu)點(diǎn)是測(cè)量精度高，能夠直接測(cè)得地下水的補(bǔ)給并且可以長期監(jiān)測(cè)。該方法通常作為估計(jì)地下水補(bǔ)給量的基準(zhǔn)。存在的不足主要有：①不適用于植被根系埋藏過深的地區(qū)；②代表性的問題；③建造和維護(hù)都需要很大的人力、物力和財(cái)力；④從建造完成到得到真實(shí)可靠數(shù)據(jù)，需要較長的時(shí)間。由于這些方面的原因使得其推廣受到一定程度的限制。

1.2 物理方法

1.2.1 零通量面方法

零通量面（又稱為ZFP）是指土壤中通過該面的水分通量為零的平面或曲面。零通量面方法“就是利用土壤的天然含水率和水勢(shì)剖面資料，分析水勢(shì)運(yùn)動(dòng)方向，找出零通量面所在的位置，計(jì)算地下水補(bǔ)給量［6］”。

Richards（1956）等人首先應(yīng)用該方法研究土壤水分通量的變化并推導(dǎo)出土壤非飽和導(dǎo)水率，隨后該方法被廣泛應(yīng)用在干旱-半干旱地區(qū)的地下水補(bǔ)給研究中。Sharma等（1991）運(yùn)用ZFP方法評(píng)價(jià)澳大利亞西部半干旱地區(qū)包氣帶中不同深度的地下水補(bǔ)給強(qiáng)度。Flint A L等（2000）用ZFP方法研究美國高放廢物處置庫場址尤卡山地區(qū)，對(duì)階地、山脊、邊坡和洼溝4類地貌區(qū)用了99個(gè)中子測(cè)孔，估算大氣降水入滲補(bǔ)給。我國于1982年由原地礦部水文地質(zhì)與工程地質(zhì)研究所從國外引進(jìn)了該方法，隨后該方法得到了廣泛的應(yīng)用。雷志棟等（1988）利用定位通量法計(jì)算潛水入滲補(bǔ)給量和土壤水分蒸發(fā)量［7］。邱景唐（1992）對(duì)零通量面的類型及其發(fā)生、遷移和消失的規(guī)律進(jìn)行了綜合分析，并提出了影響零通量面發(fā)生變化的主要因素［8］。

零通量面方法的主要優(yōu)點(diǎn)有：①可通過土壤水分直接測(cè)得下滲補(bǔ)給量；②能夠追蹤連續(xù)的補(bǔ)給率及其時(shí)空變化。該方法最大的缺點(diǎn)就是有可能出現(xiàn)零通量面消失的情況，如出現(xiàn)單純的下滲剖面時(shí)，由于濕鋒面的不斷向下推進(jìn)，零通量面可能消失。當(dāng)這種情況發(fā)生時(shí)通常采用達(dá)西方法或者水均衡方法。另外，該方法的費(fèi)用相對(duì)比較昂貴且數(shù)據(jù)采集量很大。該方法適合應(yīng)用在土壤含水量變化大、水位埋深比較大的地區(qū)。

1.2.2 達(dá)西方法

達(dá)西方法的優(yōu)點(diǎn)是其計(jì)算所需要的參數(shù)都可以通過野外觀測(cè)獲得，并且其既適用于飽水帶也適用于包氣帶。其難點(diǎn)在于確定滲透系數(shù)的空間分布和獲得可靠的水力梯度值。和零通量面方法相比，達(dá)西方法的好處就是可以常年使用，不存在“失效期”的問題。

1.2.3 地下水位動(dòng)態(tài)方法

地下水位動(dòng)態(tài)方法是一種傳統(tǒng)的以飽水帶為研究對(duì)象，基于水均衡原理確定地下水補(bǔ)給量的方法，該方法假設(shè)地下水位的上升全部是由于補(bǔ)給引起的。Meinzer（1923）、Meinzer和Stearns（1929）最早應(yīng)用該方法估計(jì)地下水的補(bǔ)給量，隨后很多學(xué)者對(duì)該方法進(jìn)行了更深入的研究。Baumann（1952）發(fā)現(xiàn)地下水位隨補(bǔ)給的變化依賴于入滲率和入滲歷時(shí)，及其他一些因素如流域的形狀與大小、含水層特性等等。Sophocleous（1991）提出了基于土壤水平衡的復(fù)合地下水位動(dòng)態(tài)方法，使估算誤差大大減小。Richard W等（2002）對(duì)該方法進(jìn)行了詳細(xì)描述，并討論了該方法在裂隙基巖地區(qū)中的應(yīng)用［10］。

這種方法的優(yōu)點(diǎn)是不需要查明地下水的補(bǔ)給機(jī)理，使用方便，參數(shù)獲取相對(duì)簡單。由于這些優(yōu)點(diǎn)使得該方法在干旱-半干旱地區(qū)得到了廣泛應(yīng)用。該方法的局限性體現(xiàn)在：①最好應(yīng)用在短期內(nèi)地下水水位變化明顯的地區(qū)；②代表性的問題，即使在同一個(gè)地方，不同觀測(cè)井的地下水位也可能由于地質(zhì)條件的不同而有較大的差異；③該方法不適用于穩(wěn)定流條件下；④很難確定給水度。

1.3 同位素方法

20 世紀(jì)60年代開始，同位素方法才相繼應(yīng)用于地下水補(bǔ)給研究，由于同位素方法的優(yōu)越性使其成為干旱-半干旱地區(qū)研究地下水補(bǔ)給最為普遍的一種方法。Allison（1988）詳細(xì)介紹了同位素方法在地下水補(bǔ)給研究中的應(yīng)用［11］。目前，用于研究地下水補(bǔ)給的同位素示蹤劑可以分為3種：歷史示蹤劑、環(huán)境示蹤劑和人工示蹤劑。

1.3.1 歷史示蹤劑

歷史示蹤劑是指歷史時(shí)期內(nèi)環(huán)境中產(chǎn)生的某些化學(xué)物質(zhì)，或者是因?yàn)槿祟惢顒?dòng)生成的物質(zhì)，或者是核試驗(yàn)產(chǎn)生的放射性物質(zhì)。目前用于研究地下水補(bǔ)給的主要是由核爆試驗(yàn)產(chǎn)生的放射性同位素3H和36Cl。

Zimmermann等（1967）首先應(yīng)用土壤中核爆3H的峰值位移，對(duì)地下水的補(bǔ)給和蒸發(fā)進(jìn)行了研究。20世紀(jì)70～80年代，該方法得到了普遍的應(yīng)用。Dincer（1974）等根據(jù)鉆取土心3H的峰值確定了沙特阿拉伯達(dá)赫納沙丘地區(qū)的補(bǔ)給量為23mm/a。我國學(xué)者也做了大量研究，但主要都集中在黃土包氣帶中。張之淦等（1990）將核爆3H應(yīng)用于山西黃土包氣帶水分運(yùn)移及降水入滲研究［12］。武清華等（2000）利用土壤包氣層實(shí)驗(yàn)，研究了3H示蹤劑在黃土包氣層中的遷移規(guī)律。Lin和Wei（2001）于1998年再次對(duì)該處包氣帶進(jìn)行了環(huán)境同位素剖面研究，進(jìn)一步證實(shí)了前者的成果。

由于核爆試驗(yàn)的停止和自身的衰減（半衰期只有12.32 a），目前環(huán)境中的3H含量已經(jīng)很低且3H不易被測(cè)定使得該方法的應(yīng)用受到一定限制。只有在滲透性差、補(bǔ)給強(qiáng)度低或者包氣帶厚度比較大的地方才適用。

36Cl也是研究地下水補(bǔ)給的優(yōu)良示蹤劑，其主要來源于核試驗(yàn)，并且36Cl在土壤水和地下水中的對(duì)流、彌散、擴(kuò)散過程中其守恒性質(zhì)優(yōu)于3H。3H在蒸散發(fā)氣態(tài)下會(huì)有所損失，而36Cl在正常條件下完全不揮發(fā)，只因水分蒸發(fā)損失而增加其在水中的濃度，在包氣帶中發(fā)生積累。Marianne Guerin（2001）應(yīng)用3H 和36Cl示蹤劑對(duì)美國高放廢物處置場址尤卡山地區(qū)包氣帶中的水分運(yùn)移情況進(jìn)行了研究［13］。

1.3.2 環(huán)境示蹤劑

環(huán)境示蹤劑方法是研究干旱、半干旱地區(qū)地下水補(bǔ)給方面最為廣泛的一種方法。該方法具有直接參與地下水循環(huán)、不污染環(huán)境、成本低廉、不需要大量觀測(cè)資料以及能夠追蹤長時(shí)間尺度地下水運(yùn)動(dòng)等優(yōu)點(diǎn)。環(huán)境示蹤劑既適用于包氣帶也適用于飽水帶。常用的環(huán)境示蹤劑有Cl-、穩(wěn)定同位素（如D、18O）等。

包氣帶Cl-示蹤的特性最早是由澳大利亞水文地質(zhì)學(xué)家Anderson于1945年提出的。Eriksson等（1969）應(yīng)用Cl-示蹤的特性研究以色列海岸平原地區(qū)的地下水補(bǔ)給，由于其忽略海水入侵的影響，因此其準(zhǔn)確性未能使人信服。Allison等（1978）經(jīng)過修正提出了利用包氣帶Cl-計(jì)算地下水補(bǔ)給量的氯質(zhì)量平衡法（CMB方法）和氯離子剖面法（CPM方法）。Allison（1988）提出了利用Cl-濃度剖面計(jì)算入滲補(bǔ)給量的另一種方法：累積法。Subyani和Zekai（2006）在研究沙特阿拉伯地區(qū)的地下水補(bǔ)給時(shí)將攝動(dòng)方法引入到CMB方法中，完成了對(duì)CMB方法的改進(jìn)，提高了計(jì)算的精度［14］。我國學(xué)者也對(duì)CMB方法的應(yīng)用做了深入的研究。陳植華等（1996）對(duì)Cl-示蹤法在干旱-半干旱地區(qū)的應(yīng)用作了詳細(xì)介紹。陳宗宇等（2001）利用CMB方法探討識(shí)別包氣帶剖面所代表的古水文-氣候記錄的方法，恢復(fù)研究區(qū)的古水文補(bǔ)給歷史［15］。馬金珠等（2004）利用CMB方法研究巴丹吉林沙漠地區(qū)的800 a以來的地下水補(bǔ)給量及其所反映的氣候變化特征［16］。黃天明等（2010）利用改進(jìn)的CMB方法研究黃土高原地區(qū)土地利用變化對(duì)地下水補(bǔ)給的影響［17］。聶振龍等（2011）應(yīng)用CMB方法求得張掖盆地地下水位埋深大于5m時(shí)仍存在降水入滲補(bǔ)給。

Cl-示蹤劑應(yīng)用的局限性主要依賴以下幾個(gè)方面：①氯離子輸入來源的多樣性；②氯離子背景值的不確定性；③降水入滲補(bǔ)給方式的不確定性；④植被的吸附作用［18］。

穩(wěn)定同位素D、18O也常用于大氣降水入滲補(bǔ)給研究。Thoma等（1978）利用大氣降水中18O的季節(jié)性變化求解以色列的Pelat沙丘地區(qū)的降水入滲量，并與3H求出的結(jié)果進(jìn)行對(duì)比，發(fā)現(xiàn)結(jié)果基本一致［19］。我國學(xué)者也開展了這方面的研究：李發(fā)東等（2007）應(yīng)用D、18O示蹤劑計(jì)算太行山地區(qū)的入滲補(bǔ)給量［20］；王福剛等（2007）應(yīng)用D、18O示蹤劑求解河南新鄉(xiāng)地區(qū)的大氣降水入滲補(bǔ)給量，并通過野外地中滲透儀方法和降水入滲系數(shù)方法證明了該方法的可靠性。

1.3.3 人工示蹤劑

人工示蹤劑是指人工注入包氣帶中的示蹤劑。一般被投放到地表或地下一定深度，根據(jù)示蹤劑峰值的變化來推求入滲補(bǔ)給量。常用的主要有Br-、3H和染色劑。其中染色劑經(jīng)常用來評(píng)價(jià)優(yōu)先流。

國內(nèi)、外對(duì)人工示蹤劑在地下水補(bǔ)給中的應(yīng)用都開展了大量的研究。Rice等（1986）應(yīng)用Br-示蹤劑研究裸地在灌溉和降水條件下的補(bǔ)給強(qiáng)度。Rangarajan等（2000）應(yīng)用3H作為人工示蹤劑選取25個(gè)典型代表區(qū)確定印度地區(qū)的排泄量為24～198 mm/a。陳宗宇等（2003）應(yīng)用3H示蹤劑研究正定試驗(yàn)場的入滲補(bǔ)給強(qiáng)度為0.095 m/a。汪丙國等（2008）運(yùn)用人工示蹤劑Br-、3H確定河北平原地區(qū)的入滲補(bǔ)給強(qiáng)度為127.8mm/a。

Sr同位素也可以用于估計(jì)地下水補(bǔ)給，其代表的時(shí)間尺度從數(shù)百年到數(shù)千年，該方法在干旱-半干旱地區(qū)特別適用。國外對(duì)該方法進(jìn)行了一定的探討，Maher等（2004）應(yīng)用Sr同位素研究美國華盛頓州漢福德地區(qū)包氣帶中的滲透速率。我國關(guān)于Sr同位素的研究成果很多，但是，關(guān)于Sr同位素研究地下水補(bǔ)給方面的報(bào)道幾乎沒有。因此，這一方法在我國具有很大的發(fā)展前景［21］。

示蹤劑的選擇應(yīng)遵循的原則有：①在包氣帶中不易被吸收、不發(fā)生化學(xué)變化；②在環(huán)境中背景值低；③測(cè)試和操作較簡單等。人工示蹤劑方法具有定時(shí)、定位和定量投放的優(yōu)點(diǎn)，且不需要測(cè)定水文地質(zhì)參數(shù)，花費(fèi)少。該方法的缺點(diǎn)是其代表的時(shí)空尺度有限，不能克服地層在空間上的變化，某些示蹤劑易污染環(huán)境，并且該方法由于忽視優(yōu)先流的補(bǔ)給容易使其評(píng)價(jià)量偏小。

1.4 數(shù)值模擬方法

隨著計(jì)算機(jī)技術(shù)的發(fā)展，數(shù)值模擬方法已成為計(jì)算地下水補(bǔ)給量的有效方法之一。目前用于計(jì)算地下水補(bǔ)給的模型主要有：①包氣帶水均衡模型（土壤水均衡模型、基于Richards方程的模型）；②流域模型；③地下水流動(dòng)模型；④流域模型和地下水流動(dòng)模型相結(jié)合的方法［22］。常用的用于包氣帶水均衡模型的軟件主要有：VADOSE/W、HYDRUS、WHI UnSatSuite和TOUGH2等。

國內(nèi)、外對(duì)數(shù)值模擬方法從理論到應(yīng)用都開展了大量研究。Ward Sanford（2002）詳細(xì)敘述了在地下水流模型的校正和反演過程中利用水力梯度、滲透系數(shù)和其他參數(shù)等信息預(yù)測(cè)地下水補(bǔ)給率［23］。O’Reilly（2004）利用土壤水均衡模型估計(jì)美國佛羅里達(dá)州中部地區(qū)的地下水補(bǔ)給量。我國學(xué)者雷志棟等（1982）在有限元分析基礎(chǔ)上編制了適合于DJ—130機(jī)的非飽和一維流動(dòng)BASIC計(jì)算程序。楊建鋒（2005）對(duì)包氣帶水和溶質(zhì)運(yùn)移數(shù)值模型進(jìn)行了詳細(xì)的介紹。

數(shù)值模擬理論上適用于各種條件的模擬，并且其不受時(shí)間尺度和空間尺度的限制，利用該方法還可以預(yù)測(cè)地下水補(bǔ)給隨氣候、土地利用以及其他因素的變化。數(shù)值模擬最大的難點(diǎn)在于精確獲取其中的參數(shù)。

2 高放廢物地質(zhì)處置中的地下水補(bǔ)給研究

高放廢物地質(zhì)處置的基本概念是“把玻璃固化后的核廢物裝入容器，再將其置入地下深部的基巖洞室中，同時(shí)以膨潤土緩沖回填材料和密封材料回填，即，設(shè)置多重屏障以阻止放射性物質(zhì)釋放、遷移和返回生物圈”［24］。因此，放射性廢物處置庫系統(tǒng)是一個(gè)能將廢物同人類和環(huán)境隔離開來的天然和工程屏障系統(tǒng)。在天然系統(tǒng)內(nèi)，地下水的運(yùn)移是放射性核素到達(dá)人類環(huán)境中最有可能的途徑。因而，在選擇放射性廢物地質(zhì)處置場址并進(jìn)行特性評(píng)價(jià)時(shí)，水文地質(zhì)是必須考慮的一個(gè)重要因素。

然而，考慮隔離的時(shí)間段超過數(shù)千年甚至上萬年，國際上普遍使用數(shù)學(xué)模型來研究整個(gè)水文地質(zhì)系統(tǒng)，而地下水補(bǔ)給量是模型其中的源匯項(xiàng)，也是構(gòu)成處置庫水文地質(zhì)特性評(píng)價(jià)的數(shù)據(jù)基礎(chǔ)［25-26］。目前，美國、瑞典和德國等一些有核國家都在進(jìn)行放射性廢物處置庫場址研究。這些國家考慮的放射性廢物處置的地質(zhì)環(huán)境不同，因而，水文地質(zhì)系統(tǒng)及對(duì)未來行為預(yù)測(cè)中存在的不確定性也不同，但候選處置場址的特性評(píng)價(jià)中需要了解的水文地質(zhì)系統(tǒng)特性的信息是相似的，采用的途徑和研究的方法也是相似的。Flint（2002）等在研究美國高放廢物處置場址尤卡山地區(qū)的入滲補(bǔ)給量時(shí)，綜合應(yīng)用了前述的物理方法、經(jīng)驗(yàn)值方法、環(huán)境示蹤劑方法以及數(shù)值模擬方法，獲得了可靠數(shù)據(jù)［27］。

我國高放廢物地質(zhì)處置庫的選址工作，起始于1985年，工作重點(diǎn)在甘肅北山地區(qū)，其氣候干旱，年降雨量僅為幾十毫米，地表水和地下水都十分貧乏［28］。對(duì)于預(yù)選場址而言具有可利用數(shù)據(jù)少、水文地質(zhì)條件復(fù)雜以及要求精度高等特點(diǎn)。筆者認(rèn)為，在我國高放廢物處置庫選址過程中，對(duì)于地下水的入滲補(bǔ)給研究，可以借鑒國外經(jīng)驗(yàn)，把同位素方法、地下水位動(dòng)態(tài)方法、達(dá)西方法和數(shù)值模擬方法相結(jié)合，以獲取可靠數(shù)據(jù)。同時(shí)還應(yīng)嘗試?yán)眯录夹g(shù)、新方法，如：遙感技術(shù)、地球物理技術(shù)和惰性氣體方法等，從不同角度揭示影響地下水入滲補(bǔ)給的因素，為后續(xù)建立大尺度概念模型提供科學(xué)的參考依據(jù)。

3 結(jié)語

通過調(diào)研大量國內(nèi)、外資料，總結(jié)出干旱-半干旱地區(qū)有關(guān)地下水補(bǔ)給研究的主要方法：①直接測(cè)定方法；②物理方法；③同位素方法；④數(shù)值模擬方法。由于每種地下水補(bǔ)給研究方法都有各自相應(yīng)的優(yōu)、缺點(diǎn)以及適用性，加上干旱-半干旱地區(qū)水文地質(zhì)條件的復(fù)雜性，因此，在選取地下水補(bǔ)給研究方法時(shí)需要考慮的原則主要有：研究目的、補(bǔ)給類型/機(jī)理、可利用/可獲取的數(shù)據(jù)、準(zhǔn)確度、時(shí)間/空間以及經(jīng)費(fèi)情況等。例如：在人力、物力、財(cái)力都比較充足且精度要求比較高時(shí)，可以采用地中滲透儀方法；在地下水位埋深大、土壤含水量變化大的地區(qū)可以采用零通量面方法；在經(jīng)費(fèi)較少、地下水補(bǔ)給機(jī)理相對(duì)復(fù)雜尤其是在存在優(yōu)先流的情況下，可以考慮地下水位動(dòng)態(tài)方法。在地質(zhì)條件復(fù)雜、研究程度低、經(jīng)費(fèi)較少時(shí)應(yīng)優(yōu)先選用同位素方法；若要估算大尺度的地下水補(bǔ)給，應(yīng)采用以飽水帶為研究對(duì)象的方法：物理方法、同位素方法和數(shù)值模擬方法；如果要預(yù)測(cè)萬年以后大氣降水入滲補(bǔ)給特征，則必須采用數(shù)值模擬方法。

［1］BridgetR.Scanlon，Kelley E.Keese，Alan L.Flint，et al.Global synthesis of groundwater recharge in semiaridand arid regions［J］.Hydrological Processes，2006:3 335-3 370.

［2］王大純，張人權(quán)，史毅虹，等.水文地質(zhì)學(xué)基礎(chǔ)［M］.北京：地質(zhì)出版社，2006.

［3］Bridget R.Scanlon，Richard W，Healy Peter G Cook.Choosing appropriate techniques forquantifying groundwater recharge［J］.Hydrogeology Journal，2002，10:18-39.

［4］Aboukhaied A，lfaro A A，Smith M.Lysimeters Food and Agriculture organization of the United Nations ［J］.FAO Irrigation and Drainage Paper，1981:39.

［5］Jinquan Wu，Renduo Zhang，Jinzhong Yang. Estimating infiltration recharge using a response functionmodel［J］.Journal of Hydrology，1997，198: 124-139.

［6］顧慰祖.同位素水文學(xué)［M］.北京：科學(xué)出版社，2011.

［7］雷志棟，楊詩秀，謝森傳.田間土壤水量平衡與定位通量法的應(yīng)用［J］.水利學(xué)報(bào)，1988，5:1-7.

［8］邱景唐.非飽和土壤水零通量面的研究［J］.水利學(xué)報(bào)，1992，5:27-32.

［9］Sibanda T，Nonner JC，Uhlenbrook S.Comparison of groundwater recharge estimation methods for the semi-arid Nyamandhlovu area，Zimbabwe［J］. Hydrogeology Journal，2009，17:1 427-1 441.

［10］Richard W Healy，Peter G Cook.Using groundwater levels to estimate recharge［J］.Hydrogeology Journal，2002，10:91-109.

［11］Allison，G B.A review of some the physical，chemical，and isotopic techniques available for estimating groundwater recharge［J］.Estimation of Natural Groundwater Recharge，1987:49-72.

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

［13］Marianne Guerin.Tritium and36Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain［J］.Journal of Contaminant Hydrology，2001:257-288.

［14］Ali Subyani，Zekai Sen.Refined chloride massbalancemethod and itsapplication in SaudiArabia［J］. HydrologicalProcesses，2006，20:4 373-4 380.

［15］陳宗宇，畢二平，聶振龍，等.包氣帶剖面中古水文-氣候信息的初步研究［J］.地球?qū)W報(bào)，2001，22 （4）:335-339.

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

［17］Tianming Huang，Zhonghe Pang.Estimating groundwater recharge following land-use changeusing chloridemassbalance ofsoilprofiles：A case study at Guyuan and Xifeng in the Loess Plateau of China［J］. Hydrogeology Journal，2011，19:177-186.

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

［19］Thoma G，Esser N，Sonntag C，et al.New technique of insitu soil moisture sampling for environmental isotope analysis applied at Pilot sand dune near Bordeaux［C］//Isotope hydrology 1978 （Proc Symp Neuherberg，1978）.Vienna:IAEA，1978：191-204.

［20］Fadong Li，Xianfang Song，Changyuan Tang.Tracing infiltration and recharge using stable isotope in Taihang MT.，North China［J］.Environ Geol，2007，53:687-696.

［21］馬傳明，劉存富，周愛國，等.同位素水文學(xué)新技術(shù)新方法［M］.武漢：中國地質(zhì)大學(xué)出版社，2010:12.

［22］Richard W.Healy.Estimating groundwater recharge ［M］.Cambridge:Printed in the United Kingdom at the University Press，2010.

［23］Ward Sanford.Recharge and groundwater models：An overview［J］.Hydrogeology Journal，2002，10: 110-120.

［24］王駒，范顯華，徐國慶，等.中國高放廢物地質(zhì)處置十年進(jìn)展［M］.北京：原子能出版社，2004.

［25］郭永海，王駒.高放廢物地質(zhì)處置中的地質(zhì)、水文地質(zhì)、地球化學(xué)關(guān)鍵科學(xué)問題［J］.巖石力學(xué)與工程學(xué)報(bào)，2007，26（2）:3 926-3 931.

［26］王駒，陳偉明，蘇銳，等.高放廢物地質(zhì)處置及其若干關(guān)鍵科學(xué)問題［J］.巖石力學(xué)與工程學(xué)報(bào)，25（4）:801-812.

［27］Flint A L，F(xiàn)lint L E，Kwicklis E M，et al. Estimating recharge at Yucca Mountain，Nevada，USA：Comparison of methods［J］.Hydrogeology Journal，2002，10:180-204.

［28］王駒，蘇銳，陳偉明，等.中國高放廢物深地質(zhì)處置［J］.巖石力學(xué)與工程學(xué)報(bào).2006，25（4）: 649-658.

Reviews on study methods of groundwater recharge in arid and sem i-arid regions

LIJie-biao，SU Rui，TIAN Xiao，ZHOU Zhi-chao
（CNNC Key Laboratory on Geological Disposal of High-level RadioactiveWaste，Beijing Research Institute of Uranium Geology，Beijing 100029，China）

Groundwater recharge is an important aspect in the study of hydrogeology.Four study methods of groundwater recharge in arid and semiarid regions were summarized in this paper:they are ①direct measurement methods；②physical methods；③isotopic methods；④numerical modeling methods.The principle，research situation，merits and demerits as well as application were discussed for eachmethod，in order to provide reference for the research of hydrogeologymethodology.

groundwater recharge；hydrogeology；arid and semi-arid regions；methodology

P641

A

1672-0636（2013）03-0168-06

2013-04-08；

2013-05-10

李杰彪（1987—），男，河北冀州人，碩士研究生，水文地質(zhì)專業(yè)，主要從事高放廢物地質(zhì)處置方面的工作。E-mail：hgylijiebiao@126.com

10.3969/j.issn.1672-0636.2013.03.008