劉穎超，邢國(guó)章，劉 凱，謝振華

（1.北京市水文地質(zhì)工程地質(zhì)大隊(duì)，北京 100195；2.北京市地質(zhì)工程勘察院，北京 100048）

北京典型地區(qū)灌溉條件下包氣帶水分運(yùn)移特征研究

劉穎超1，邢國(guó)章2，劉 凱1，謝振華1

（1.北京市水文地質(zhì)工程地質(zhì)大隊(duì)，北京 100195；2.北京市地質(zhì)工程勘察院，北京 100048）

以張家灣包氣帶水鹽運(yùn)移試驗(yàn)場(chǎng)灌溉試驗(yàn)資料為依據(jù)，運(yùn)用土壤水分能量觀點(diǎn)，分析了灌溉條件下包氣帶水分運(yùn)移特征。試驗(yàn)結(jié)果表明，在一定灌溉量下包氣帶水分運(yùn)移具有明顯的分帶特征以及灌溉水明顯運(yùn)移到達(dá)的深度；在一定灌溉量條件下，灌溉后地下水位上升微小，灌溉補(bǔ)給地下水量較小，并通過非飽和水流達(dá)西定律計(jì)算了向下入滲補(bǔ)給量。

包氣帶；含水量；土水勢(shì)；水分運(yùn)移；補(bǔ)給量

0 引言

包氣帶是地下水與大氣水、地表水相互聯(lián)系的紐帶，在廣大的平原地區(qū)，降雨、灌溉水的入滲補(bǔ)給是淺層地下水的主要補(bǔ)給來源，地下水獲得垂向補(bǔ)給必須通過包氣帶來實(shí)現(xiàn)。目前隨著地下水運(yùn)動(dòng)的研究與發(fā)展，地下水運(yùn)動(dòng)的基本理論也基本成熟。目前，在確定降雨和灌溉水入滲補(bǔ)給量的計(jì)算方面，主要采用降雨入滲補(bǔ)給系數(shù)法和灌溉水入滲補(bǔ)給系數(shù)法［1～8］。由于沒有考慮降雨、灌溉水進(jìn)入到非飽和帶土壤后的運(yùn)移機(jī)理及過程，以及補(bǔ)給系數(shù)的準(zhǔn)確性，致使地下水補(bǔ)給量的計(jì)算精度不高。近些年來，非飽和包氣帶水分運(yùn)移特征研究成為熱點(diǎn)，各地相繼建成了包氣帶水鹽運(yùn)移試驗(yàn)場(chǎng)。針對(duì)不同地區(qū)、不同巖性包氣帶進(jìn)行了各種原位試驗(yàn)研究。李曼等對(duì)太行山前平原田間土壤水分運(yùn)移機(jī)理研究［9］；李明香等對(duì)黃土包氣帶水分運(yùn)移特征進(jìn)行了研究［10］；李訓(xùn)華，商潔對(duì)沙漠地區(qū)包氣帶水分運(yùn)移特征進(jìn)行了研究［11～12］；趙貴章，段鵬鄂爾多斯盆地風(fēng)沙灘地區(qū)包氣帶水分運(yùn)移特征進(jìn)行了研究［13～14］。研究非飽和帶土壤水分運(yùn)移規(guī)律，可以準(zhǔn)確的確定降雨、灌溉對(duì)地下水的補(bǔ)給量，以此提高地下水資源量的計(jì)算精度，滿足工農(nóng)業(yè)生產(chǎn)發(fā)展的需要。本文依據(jù)張家灣包氣帶原位水鹽運(yùn)移試驗(yàn)場(chǎng)進(jìn)行的灌溉試驗(yàn)，運(yùn)用土壤水能量理論分析了灌溉條件下包氣帶水分運(yùn)移的特征。

1 試驗(yàn)場(chǎng)概況

包氣帶水鹽運(yùn)移試驗(yàn)場(chǎng)位于北京通州區(qū)張家灣鎮(zhèn)國(guó)土資源部平原區(qū)地下水—北京野外基地內(nèi)。該試驗(yàn)場(chǎng)地處潮白河沖洪積扇的中下部，第四系主要為河流沖洪積作用形成，厚度在150m左右，代表著北京平原區(qū)沖洪積扇中下部大部分區(qū)域的表層地質(zhì)環(huán)境。包氣帶巖性以砂質(zhì)粉土、粉砂夾粘質(zhì)粉土為主。試驗(yàn)場(chǎng)多年降雨量為556.2mm，多年地下水位埋深3～4m。

試驗(yàn)場(chǎng)南北長(zhǎng)17m，東西寬12m，試驗(yàn)場(chǎng)內(nèi)修建地下觀測(cè)室一處，在觀測(cè)室北、東、南方向分別開挖3個(gè)試驗(yàn)儀器埋設(shè)槽，根據(jù)地層巖性在不同深度埋設(shè)負(fù)壓計(jì)，試驗(yàn)場(chǎng)內(nèi)設(shè)有簡(jiǎn)易氣象站、水位及含水量觀測(cè)系統(tǒng)。試驗(yàn)場(chǎng)平面布置圖見圖1，試驗(yàn)場(chǎng)內(nèi)地層及負(fù)壓計(jì)埋設(shè)情況見圖2。

2 試驗(yàn)方法

本次灌溉試驗(yàn)選取南槽作為灌溉區(qū)，灌溉面積9.97 m2，灌水量1500L，試驗(yàn)水源選用涼水河河水。試驗(yàn)過程中通過WM-1型負(fù)壓計(jì)系統(tǒng)，監(jiān)測(cè)各土體深度的基質(zhì)勢(shì)；運(yùn)用the TRIME Data Pilot system監(jiān)測(cè)各深度土體含水量；同時(shí)監(jiān)測(cè)試驗(yàn)場(chǎng)內(nèi)地下水位和氣象因素。灌溉開始前對(duì)地下水位及各深度土體基質(zhì)勢(shì)、含水量進(jìn)行觀測(cè)作為背景值。灌溉開始后基質(zhì)勢(shì)每10分鐘觀測(cè)1次，7小時(shí)后每30分鐘觀測(cè)一次，10小時(shí)后每1小時(shí)觀測(cè)一次，15小時(shí)后每3小時(shí)觀測(cè)一次，灌溉7天后改為日常觀測(cè)每天2次。含水量的觀測(cè)，灌溉當(dāng)天每30分鐘觀測(cè)一次，灌溉第2天后每3小時(shí)觀測(cè)1次，灌溉7天后改為日常觀測(cè)每天1次。水位和氣象每天觀測(cè)一次。

圖1 試驗(yàn)場(chǎng)平面布置圖

圖2 試驗(yàn)場(chǎng)地層及負(fù)壓計(jì)埋設(shè)圖

3 結(jié)果與分析

灌溉后灌溉水滲入地下，儲(chǔ)存在包氣帶中，并逐漸由上至下改變土壤的勢(shì)能狀態(tài)，一次灌溉能夠真正補(bǔ)給地下水取決于多方面的因素，主要與灌溉前土壤水勢(shì)能剖面有關(guān)，還跟灌溉量與灌溉時(shí)間等有關(guān)。灌溉水向下入滲的驅(qū)動(dòng)力為土水勢(shì)，土壤水總是從高土水勢(shì)處向低土水勢(shì)處運(yùn)移，測(cè)定了土壤中連續(xù)兩點(diǎn)的土水勢(shì)以后，即可確定土壤水的運(yùn)移方向，即土壤水在水勢(shì)梯度大于零時(shí)，水分向下運(yùn)移；當(dāng)水勢(shì)梯度小于零，水分向上運(yùn)移，水勢(shì)梯度等于零時(shí)，該處的水分通量為零，稱該處的水平面為零通量面［15］，以此判斷灌溉水的運(yùn)移情況。通過對(duì)含水量和基質(zhì)勢(shì)的監(jiān)測(cè)和計(jì)算，繪制各深度土體含水量、土水勢(shì)隨時(shí)間、深度變化曲線，分析灌溉過程中土壤水分運(yùn)移特征。土壤水分運(yùn)動(dòng)呈非飽和流狀態(tài)，可近似看作一維垂向水流運(yùn)動(dòng)，遵循非飽和達(dá)西定律，通過非飽和達(dá)西定律計(jì)算灌溉向下入滲補(bǔ)給量。

圖3 灌溉當(dāng)天包氣帶含水量隨深度變化曲線

3.1 含水量變化情況

從含水量隨深度變化曲線上可以看出從地表至地表以下120 cm土體含水量逐漸減?。▓D3、圖4），受為粉土粉砂互層（130～138cm）、169～180cm粉質(zhì)粘土層（169～180cm）阻水作用，地表以下120～200cm范圍內(nèi)土體含水量出現(xiàn)增加，220cm以下土體為毛細(xì)水作用影響帶，含水量有一定波動(dòng)。灌溉后土壤剖面上含水量分布發(fā)生明顯變化，其中0～60cm范圍內(nèi)土體吸收水分后，含水量出現(xiàn)明顯增大，變幅自地表依次減小，剖面上遇到相對(duì)阻水巖層后含水量會(huì)明顯增加（圖5）。

圖4 灌溉后包氣帶含水量隨深度變化曲線

圖5 包氣帶各深度土體含水量隨時(shí)間變化情況

3.2 土水勢(shì)變化情況

灌溉試驗(yàn)于5月6日上午9:00開始，灌溉過程中保持水面穩(wěn)定在8cm。灌溉前從土水勢(shì)隨深度變化圖上可以看出，由于蒸發(fā)強(qiáng)烈，水分處于虧缺狀態(tài)，土水勢(shì)值遠(yuǎn)離飽和水勢(shì)線(重力勢(shì))，土壤吸水能力增強(qiáng)、給水能力減弱，灌溉入滲補(bǔ)給地下水的條件最為困難。包氣帶土壤水運(yùn)動(dòng)為蒸發(fā)-入滲型，并在80cm處形成零通量面，80cm以上土體土壤水處于蒸發(fā)狀態(tài)，向上運(yùn)移；80cm以下至220cm包氣帶水處于入滲狀態(tài)，向下運(yùn)移。灌溉開始后，在入滲剖面上隨著時(shí)間推移灌溉水逐漸入滲深入包氣帶，較干燥的土壤得以充分吸水，各深度土體土水勢(shì)逐步增大接近飽和水勢(shì)線，當(dāng)土體吸持一定水分后，土體進(jìn)入過水狀態(tài)，土水勢(shì)值則繼續(xù)變大。若入滲水量有限，這時(shí)的土體含水率和土水勢(shì)都達(dá)到該次灌溉入滲過程的極值，然后土體進(jìn)入脫水狀態(tài)，土水勢(shì)值逐漸減少，慢慢遠(yuǎn)離飽和土水勢(shì)線。隨著入滲水流下移深度的增大，濕潤(rùn)鋒的水平軸距變短而垂向距變長(zhǎng)，同時(shí)濕潤(rùn)鋒到達(dá)各深度的時(shí)間間隔逐漸變大，即濕潤(rùn)鋒的運(yùn)移逐漸減慢。在灌溉當(dāng)天土水勢(shì)明顯影響深度為140cm，隨著埋深增加，各層土水勢(shì)變化幅度逐漸減少，140cm以下各深度土水勢(shì)沒有明顯變化（圖6），即灌溉當(dāng)天土壤水分以入滲型緩慢運(yùn)移到140cm。從灌溉后各深度隨時(shí)間變化曲線中可以看出為0～120cm各層土水勢(shì)隨時(shí)間推移逐漸緩慢降低；140～180cm各層土水勢(shì)隨時(shí)間推移先緩慢增大，之后長(zhǎng)時(shí)間維持穩(wěn)定，180cm以下土各層水勢(shì)維持穩(wěn)定（圖7），故此次灌溉水分運(yùn)移明顯影響深度為180cm。

圖6 灌溉當(dāng)天各包氣帶土體土水勢(shì)隨時(shí)間變化曲線

圖7 灌溉后各包氣帶土體土水勢(shì)隨時(shí)間變化曲線

從土水勢(shì)隨深度深度圖上看（圖8、圖9），灌溉后由于蒸發(fā)強(qiáng)烈土壤水分從5月6日－5月7日的入滲型迅速過渡到5月8日后的蒸發(fā)-入滲型，且其形成的零通量面隨著時(shí)間推移緩慢運(yùn)移到90cm。

從灌溉后土水勢(shì)變幅可以看出(表1)，灌溉后土壤水分運(yùn)移明顯分為3個(gè)帶：0～70cm土體變幅大于5cm Hg，土水勢(shì)、水勢(shì)梯度隨時(shí)間變化快，灌溉水運(yùn)移明顯，為水分運(yùn)移的強(qiáng)交替帶；70～180cm范圍內(nèi)土體土水勢(shì)變幅為0～5cm Hg，土水勢(shì)、水勢(shì)梯度隨時(shí)間變化較大，灌溉水運(yùn)移也比較明顯，為水分運(yùn)移弱交替帶；180～220cm范圍內(nèi)土體土水勢(shì)變幅小于1cm Hg，土水勢(shì)、水勢(shì)梯度隨時(shí)間變化很小，基本處于穩(wěn)定狀態(tài)，灌溉水運(yùn)移不太明顯，為水分運(yùn)移基本穩(wěn)定帶；220～330cm為毛管水帶。

圖8 灌溉當(dāng)天包氣帶土水勢(shì)隨深度變化曲線

圖9 灌溉后包氣帶土水勢(shì)隨深度變化曲線

3.3 入滲補(bǔ)給量計(jì)算

從地下水位隨時(shí)間變化圖上看（圖10），灌溉后引起地下水位的微小上升，而后又慢慢下降，表明灌溉補(bǔ)給地下水量微小，而5月18日16.1mm的降雨引起地下水位的明顯抬升，表明大面積的降雨通過地層優(yōu)先流或其它途徑能夠?qū)Φ叵滤倪M(jìn)行快速補(bǔ)給。

根據(jù)土水勢(shì)變化分析，本次試驗(yàn)明顯入滲影響深度為180cm，由于169～200cm為粘質(zhì)粉土層，入滲剖面上土體粒徑由粗變細(xì)，對(duì)土的下滲具有很強(qiáng)的阻滯作用，入滲能力變小，但180～220cm土體間存在土水勢(shì)梯度，土壤水在水勢(shì)梯度的作用下還能緩慢下滲，根據(jù)非飽和水流達(dá)西定律（公式1）計(jì)算灌溉當(dāng)月通過該層的水量。

表1 各深度土體土水勢(shì)變化情況表

式中q—土壤水分通量

k(θ)—以土壤含水率 表示的非飽和土壤導(dǎo)水率

k(ψm)—以土壤基質(zhì)勢(shì) 表示的非飽和土壤導(dǎo)水率

根據(jù)非飽和達(dá)西定律計(jì)算5月1日－31日土壤水由180cm到220cm通過粘質(zhì)粉土層向下入滲量約為1.6mm。該粘質(zhì)粉土層位于土壤剖面下部，以下都為粉砂土，處于毛細(xì)水上升帶附近，通過該層的水量可近似的估計(jì)為該時(shí)間段內(nèi)灌溉水入滲補(bǔ)給量。

圖10 地下水位隨時(shí)間變化曲線

4 結(jié)論

試驗(yàn)場(chǎng)地灌溉條件下水分運(yùn)移特征具有如下顯著的特點(diǎn)：

（1）灌溉后包氣帶水分運(yùn)移經(jīng)歷了蒸發(fā)–入滲型、入滲型、蒸發(fā)–入滲型3種運(yùn)動(dòng)狀態(tài)的轉(zhuǎn)變。

（2）根據(jù)灌溉后土水勢(shì)變幅，試驗(yàn)場(chǎng)內(nèi)土體水分運(yùn)移明顯分為3個(gè)帶，即強(qiáng)交替帶、弱交替帶、基本穩(wěn)定帶。

（3）通過土水勢(shì)隨時(shí)間的變化，確定本次灌溉試驗(yàn)水分明顯運(yùn)移至180cm。灌溉后引起試驗(yàn)場(chǎng)內(nèi)地下水位微小上升，補(bǔ)給地下水量小，而5月18日降雨后，引起了試驗(yàn)場(chǎng)內(nèi)地下水位的快速上升，表明試驗(yàn)場(chǎng)內(nèi)灌溉或降雨都能補(bǔ)給地下水。

（4）當(dāng)包氣帶中存在粘質(zhì)粉土層或粉質(zhì)粘土層時(shí)時(shí)對(duì)水分向下入滲具有很強(qiáng)的阻滯作用，但土體間仍存在土水勢(shì)梯度，通過非飽和達(dá)西定律計(jì)算得出灌溉當(dāng)月入滲補(bǔ)給地下量約為1.6mm。

［1］韓占成,溫 浩.滄州市年降雨人滲補(bǔ)給系數(shù)分析［J］.河北水利科技,1997,18(3):11～13.

［2］張榮旺.阿拉善盟賀蘭山西麓水資源調(diào)查研究區(qū)降水入滲補(bǔ)給系數(shù)α的確定［J］.內(nèi)蒙古水利,2009,(2):50～51.

［3］劉延璽,朝倫巴根,馬 龍等.規(guī)定時(shí)段降雨入滲補(bǔ)給系數(shù)的確定［J］.內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào),2002,23(2):44～48.

［4］牛振紅.降水入滲補(bǔ)給系數(shù)的實(shí)驗(yàn)研究與分析計(jì)算［J］.地下水,2003,25(3):152～154.

［5］王言思.內(nèi)蒙孿井灌區(qū)地下水補(bǔ)給的研究［D］.青島.中國(guó)海洋大學(xué),2009.

［6］周金龍,王 斌.新疆平原區(qū)灌溉水入滲補(bǔ)給地下水試驗(yàn)分析［J］.地下水,2002,24(4):202～203.

［7］謝永玉.淮北平原降雨入滲補(bǔ)給系數(shù)隨地下水埋深變化特征［J］.地下水,2012(1):9～11.

［8］李 娟,葉 立,張維江.寧夏中部干旱風(fēng)沙區(qū)降雨入滲補(bǔ)給系數(shù)的估算［J］.水土保持研究,2009(1):22～26.

［9］李 曼,何雨江,藺文靜等.太行山前平原田間土壤水分運(yùn)移機(jī)理研究［J］.干旱區(qū)資源與環(huán)境,2014,28(3):101～106.［10］李明香,馬炳輝,劉錦峰等.黃土包氣帶水分運(yùn)移的現(xiàn)場(chǎng)研究［J］.輻射防護(hù),2000,20(1-2):91～99.

［11］李訓(xùn)華.沙漠包氣帶水分運(yùn)移機(jī)理研究［J］.煤田地質(zhì)與勘探,2004,32(5):38～39.

［12］商 潔.巴丹吉林沙漠腹地包氣帶水分運(yùn)移研究［D］.北京.中國(guó)地質(zhì)大學(xué)（北京）,2014.

［13］趙貴章.鄂爾多斯盆地風(fēng)沙灘地區(qū)包氣帶水—地下水轉(zhuǎn)化機(jī)理研究［D］.西安.長(zhǎng)安大學(xué)，2011.

［14］段 鵬,王文科,藏 鵬.鄂爾多斯盆地風(fēng)沙灘地區(qū)包氣帶水分運(yùn)移特征研究［J］.地下水,2013,35(1):4～6.

［15］荊恩春,費(fèi)瑾,張孝和等.土壤水分通量法試驗(yàn)研究［M］.北京:地震出版社.1994.

Research on the Characteristics of Water M ovement under Irrigation Condition in Vadose Zone in the Typical Area of Beijing

LIU Yingchao1, XING Guozhang2,LIU Kai1,XIE Zhenhua1

（Beijing Institu te of Hydrogeo logy and Engineering Geo logy, Beijing 100195；
Beijing Institu te of Geo logical Engineering, Beijing 100048）

Base on the irrigation experiment of moisture and salt content m igration of the vadose zones in-situ in the Zhangjiawan area, the paper analyzes the m igration characteristics of moisture movement of unsaturated zones from the view of moisture energy. It shows that moisture m igration of vadose zones can be roughly divided into three zones and the process of groundwater recharge under irrigation, the deep of the moisture m igration. The ground water level rises very low under the condition of irrigation and the irrigation recharge quantity was very little, which was calculated by unsaturated fl ow Darcy's law.

Vadose zone；Water content；Soil water potential；Moisture movement；Recharge quantity

P641.131

A

1007-1903（2015）01-0020-05

劉穎超（1980–），男，工程師，主要從事水文地質(zhì)、工程地質(zhì)、地?zé)岬刭|(zhì)方面研究。E-mail:lyc95588@163.com