喬 岡，王文科

1.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，西安 7100542.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，西安 710054

西北干旱內(nèi)陸盆地區(qū)裸土蒸發(fā)強度

喬 岡1，2，王文科1

1.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，西安 710054
2.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，西安 710054

鑒于包氣帶蒸發(fā)的水分是從地表散失，故查明裸土蒸發(fā)規(guī)律對研究包氣帶水分蒸發(fā)、潛水蒸發(fā)問題具有重要意義。為了查明西北干旱內(nèi)陸盆地區(qū)裸土蒸發(fā)強度規(guī)律，通過開展裸土含水率觀測、包氣帶負(fù)壓觀測等野外原位試驗，結(jié)合氣象資料，對包氣帶剖面流場特征進行分析，采用有限差分法對裸土蒸發(fā)強度進行計算。結(jié)果表明：依據(jù)包氣帶水分在垂向上遷移的不同狀態(tài)，將裸土區(qū)包氣帶剖面劃分為水分向上傳輸帶和水分向下遷移帶2個帶；西北內(nèi)陸干旱盆地裸土包氣帶極限蒸發(fā)深度為1.2m；受外界蒸發(fā)能力和包氣帶向上輸水能力交替控制作用，盆地內(nèi)裸土蒸發(fā)出現(xiàn)在4月下半月至9月下半月，其中4月下半月至5月裸土蒸發(fā)強度最大，平均蒸發(fā)強度達0.56mm/d；裸土蒸發(fā)強度平均僅占水面蒸發(fā)強度的3.5%，最大達6.7%。

裸土；蒸發(fā)強度；總水勢；包氣帶；干旱內(nèi)陸盆地

0 引言

西北干旱內(nèi)陸盆地降水稀少、蒸發(fā)強烈，蒸發(fā)排泄已成為包氣帶水和潛水主要的天然排泄途徑，而裸土層作為大氣水分和土壤水分交換的界面，對包氣帶水分以及地下水的蒸發(fā)具有重要作用。盆地內(nèi)以黏性土、黏土為主要巖性類型的包氣帶是否存在極限蒸發(fā)深度，其值是多少以及裸土蒸發(fā)強度有多大等問題，已備受關(guān)注［1－5］。近年來，胡光成［6］、孫艷等［7］對裸土蒸發(fā)強度進行了研究。但是，從裸土蒸發(fā)機理的角度研究裸土蒸發(fā)強度仍然是薄弱環(huán)節(jié)。地處準(zhǔn)噶爾盆地南緣之天山北麓區(qū)是我國西北典型的干旱內(nèi)陸盆地之一，地處亞歐大陸腹地，遠離海洋，屬典型的大陸干旱性氣候，區(qū)內(nèi)土壤鹽漬化、荒漠化嚴(yán)重，生態(tài)環(huán)境極其脆弱。因此，研究包氣帶水分蒸發(fā)機理對在西北干旱內(nèi)陸盆地正確理解裸土蒸發(fā)過程、掌握包氣帶水動態(tài)規(guī)律以及在地下水資源評價中合理確定潛水蒸發(fā)量具有重要的指導(dǎo)意義。

1 研究方法

初步研究［8］表明，準(zhǔn)噶爾盆地地下水循環(huán)主要以垂向交換為主，受內(nèi)陸盆地區(qū)干旱氣候的影響，包氣帶水分蒸發(fā)、潛水蒸發(fā)的水分都是從地表散失。因此，本文主要采用監(jiān)測裸土區(qū)含水率、監(jiān)測包氣帶剖面負(fù)壓的方法，結(jié)合氣象要素的分析，來達到裸土蒸發(fā)強度研究的目的。為了便于查明監(jiān)測數(shù)據(jù)的變化規(guī)律，本次研究將逐日監(jiān)測數(shù)據(jù)按每半個月取平均值處理后再進行分析。

研究區(qū)每年10月初至翌年3月底為土壤冰凍期，在此期間，包氣帶水分凍結(jié)，裸土和地下水的蒸發(fā)作用很微弱。因此，原位試驗從4月初開始至9月底結(jié)束，為期6個月，共計180d，獲取裸土含水率監(jiān)測數(shù)據(jù)170個，包氣帶負(fù)壓監(jiān)測數(shù)據(jù)1 920個。

1.1 氣象要素觀測

裸土含水率的變化規(guī)律主要受氣象要素、土質(zhì)的影響［9］，因此，試驗區(qū)對日降水量、日蒸發(fā)量、平均氣溫、最高氣溫、最低氣溫、相對濕度、最小濕度、凍土深度以及包氣帶埋深0.00、0.10、0.15、0.20和0.40m位置的地溫等14個氣象要素指標(biāo)進行觀測，觀測頻率1次/d；采用篩析、密度計相結(jié)合的方法對試驗區(qū)包氣帶巖性進行顆粒分析定名，以便于對計算結(jié)果進行綜合分析。

1.2 裸土含水率測定

由于裸土含水率的變化規(guī)律綜合體現(xiàn)了包氣帶水分運移受外界影響的程度，故在原位試驗場的裸土區(qū)，負(fù)壓計探頭對應(yīng)地表的中心位置取裸土樣進行測定，取樣的頻率為1次/5d，如遇下雨天加密為1次/d，連續(xù)觀測5d，與負(fù)壓計讀數(shù)的觀測時間、觀測周期保持一致。

1.3 剖面負(fù)壓觀測

在包氣帶剖面不同埋深設(shè)置負(fù)壓觀測系統(tǒng)，依據(jù)每1個層段至少設(shè)置1支負(fù)壓計、若層較厚時適當(dāng)加密的布設(shè)原則，在包氣帶剖面上總計安裝負(fù)壓計5支，根據(jù)試坑剖面巖性顆粒分析結(jié)果以及結(jié)構(gòu)特征布設(shè)裸土區(qū)負(fù)壓計位置（圖1）。

圖1 負(fù)壓監(jiān)測位置Fig.1 Pressure monitoring location

負(fù)壓觀測頻率1次/5d，若遇雨天每天加密觀測1次，延續(xù)至5d。運用包氣帶剖面負(fù)壓變化可以定性判斷包氣帶水分的運動方向，并且可以作為裸土蒸發(fā)強度計算結(jié)果驗證的依據(jù)。

1.4 數(shù)值計算

在地表裸露的條件下，太陽總輻射、氣溫、空氣飽和差等綜合作用使得包氣帶水分散失到大氣中，從而影響剖面負(fù)壓水頭的分布。因此，可借助裸土區(qū)實測裸土含水率（計算區(qū)域的上邊界條件）和實測剖面負(fù)壓水頭資料（計算區(qū)域的初始條件和下邊界條件），采用有限差分法計算原理［10］確定一維包氣帶水分運動條件下裸土蒸發(fā)強度。

依據(jù)達西定律及有限差分法，流經(jīng)第n個單元的水分通量qzi由下式計算得到：

式中：Kn、Kn－1分別為第n個、第n－1個單元格中心節(jié)點位置處的非飽和滲透系數(shù)，由Van Geunchten［11］于1980年提出的經(jīng)驗公式計算；hn、hn－1分別為第n個、第n－1個單元格中心節(jié)點負(fù)壓水頭；zn、zn－1分別為第n個、第n－1個單元格中心節(jié)點埋深。

一維包氣帶水分通量計算差分格式見圖2。

圖2 一維包氣帶水分通量計算差分格式示意圖Fig.2 Differential scheme schematic diagram of the vadose zone 1Dwater flux computer

2 結(jié)果及討論

2.1 裸土含水率

通過對觀測數(shù)據(jù)的分析可知，年內(nèi)裸土含水率變化較大，具體表現(xiàn)為：1）9月上半月至翌年4月上半月裸土含水率總體呈增加趨勢，4月上半月裸土含水率達到最大值，達到25.0%。其原因是：受外界氣象要素變化的影響，隨著降水量的增加和地表氣溫的降低，近地表的氣態(tài)水轉(zhuǎn)化為液態(tài)水，故裸土含水率呈增加趨勢。2）4月下半月至8月下半月裸土含水率總體呈下降趨勢，含水率為4.8%～23.0%，9月份上半月裸土含水率最小，僅為4.8%。其原因是：隨著地表氣溫升高，近地表的液態(tài)水轉(zhuǎn)化為氣態(tài)水，進而使裸土含水率呈減小趨勢（圖3）。

圖3 裸土含水率、裸土地溫隨時間變化規(guī)律Fig.3 Changes curve of the bare soil moisture content and ground temperature with time

以上分析反映出包氣帶裸土含水率與氣象要素具有一定相關(guān)性。經(jīng)對實測資料與氣象要素統(tǒng)計分析，發(fā)現(xiàn)裸土含水率與0m處大于0℃的裸土地溫呈負(fù)指數(shù)相關(guān)關(guān)系（圖4）。

圖4 裸土含水率與地溫相關(guān)性曲線Fig.4 Relationship between the bare soil moisture content and ground temperature

2.2 剖面負(fù)壓

包氣帶水分在水勢的作用下，從高水勢狀態(tài)向低水勢狀態(tài)遷移。據(jù)此，根據(jù)包氣帶剖面負(fù)壓監(jiān)測數(shù)據(jù)，計算繪制了不同時段包氣帶剖面總水勢等值線圖（圖5）。依據(jù)包氣帶水分遷移的不同狀態(tài)，可將裸土區(qū)包氣帶剖面劃分為上下2個帶，分別是水分向上傳輸帶和水分向下遷移帶。由于2個帶之間的分界面上沒有水分交換，其水分通量為0，因此，該分界面則是零通量面。它可以綜合反映包氣帶水分在垂向上的流動特征，其在裸土區(qū)垂向上的位置稱為包氣帶水分蒸發(fā)深度，該深度的最大值即是極限蒸發(fā)深度。從圖5可以看出，原位試驗場裸土包氣帶極限蒸發(fā)深度為1.2m。

圖5 裸土含水率、裸土地溫及包氣帶總水勢等值線分布圖Fig.5 Distribution map among the bare soil moisture content，ground temperature and the water head isoline

2.2.1 水分向上傳輸帶

水分向上傳輸帶主要分布在零通量面以上的區(qū)域，該區(qū)包氣帶水分受外界蒸發(fā)能力及包氣帶向上輸水能力的綜合作用，水分呈向上傳輸狀態(tài)，到地表蒸發(fā)散失到大氣中，轉(zhuǎn)化為大氣水。

2.2.2 水分向下遷移帶

水分向下遷移帶主要分布在零通量面以下的區(qū)域，該區(qū)包氣帶水分不受外界蒸發(fā)能力及包氣帶向上輸水能力的作用，在自重的驅(qū)動下，水分呈向下遷移的狀態(tài)，最終到達潛水面，補給地下水。

2.3 裸土蒸發(fā)強度

水由裸露的土面直接汽化散失到空氣中稱為土面蒸發(fā)。單位時間從單位面積土表面蒸發(fā)損失的水量稱為土面蒸發(fā)強度（mm/d）。由于包氣帶蒸發(fā)的水分是從地表散失，因此為了研究包氣帶水分蒸發(fā)和潛水蒸發(fā)問題，首先必須了解裸土蒸發(fā)規(guī)律。應(yīng)用前述的包氣帶水分運移數(shù)值模型，將原位試驗點包氣帶剖面第一個節(jié)點（即地表面）通量作為實際裸土蒸發(fā)強度，所計算不同時間的裸土蒸發(fā)強度同降水量、平均氣溫的相關(guān)性曲線見圖6。

圖6 裸土蒸發(fā)強度與日降水量、平均氣溫關(guān)系圖Fig.6 Relation curve of bare soil evaporation intensity and rainfall，air temperature

由圖6可見，裸土蒸發(fā)強度與降水量、平均氣溫之間存在著密切的關(guān)系：1）4月上半月裸土蒸發(fā)微弱，降水量較大，以入滲為主；4月下半月至5月上半月受氣溫逐日升高和降水頻率增大等因素控制，裸土蒸發(fā)強度亦開始逐漸增大，且與氣溫之間存在2～3d的滯后。其原因是：這段時間內(nèi)降水量增加及前期凍土和地表積雪消融，使裸土及包氣帶含水量增加，包氣帶向上輸水能力大；而同時氣溫呈波動狀急劇升高，溫度急劇升高使裸土中更多的水分汽化散失到空氣中，裸土蒸發(fā)強度增大，此時，裸土蒸發(fā)強度主要受外界蒸發(fā)能力的控制。2）5月下半月至9月下半月，氣溫持續(xù)升高，但降水量少，地表土層和包氣帶受降水量小、含水率低、向上輸水能力弱的限制，裸土蒸發(fā)強度相對減小，變化比較平穩(wěn)。只有在較大降水強度發(fā)生2～3d后，裸土蒸發(fā)強度才出現(xiàn)明顯的增大；隨后受氣溫高等因素影響，裸土蒸發(fā)強度經(jīng)歷快速減小到逐漸緩慢減小的過程，直到下一次較大降水強度的出現(xiàn)，又重復(fù)上述過程。這種過程在這段時間內(nèi)呈交替發(fā)生，此時裸土蒸發(fā)強度主要受包氣帶向上輸水能力控制。3）9月下半月以后，氣溫持續(xù)下降，降水增加，裸土蒸發(fā)強度很小，地表土層呈入滲狀態(tài)，此時，裸土蒸發(fā)強度又開始轉(zhuǎn)變?yōu)槭芡饨缯舭l(fā)能力控制。

裸土蒸發(fā)強度受外界蒸發(fā)能力和包氣帶向上輸水能力雙重作用控制，不同時間、不同強度的降水對裸土蒸發(fā)強度影響不同：4月上半月至5月下半月小于5mm的降水量就可對裸土蒸發(fā)強度產(chǎn)生較明顯的影響；而5月下半月至9月下半月只有大于10 mm的降水量才對裸土蒸發(fā)強度產(chǎn)生較明顯的影響，小于5mm的降水量對裸土蒸發(fā)強度影響不大。

經(jīng)對計算結(jié)果的統(tǒng)計表明，區(qū)內(nèi)裸土蒸發(fā)出現(xiàn)在4月下半月至9月下半月，其中：4月下半月至5月裸土蒸發(fā)強度最大，平均蒸發(fā)強度0.46～0.56 mm/d；6月至7月裸土平均蒸發(fā)強度0.19～0.32 mm/d；8月上半月裸土平均蒸發(fā)強度為0.29mm/d；8月下半月至9月底裸土蒸發(fā)強度呈減小趨勢，至9月底平均蒸發(fā)強度只有0.13mm/d（圖7）。裸土蒸發(fā)強度平均占水面蒸發(fā)強度的3.5%，最大可達到6.7%，出現(xiàn)在5月份（圖8）。

3 結(jié)論

1）依據(jù)包氣帶水分在垂向上遷移的不同狀態(tài)，可將裸土區(qū)包氣帶剖面劃分為水分向上傳輸帶和水分向下遷移帶2個帶；受干旱氣候的影響，對于黏性土質(zhì)而言，裸土蒸發(fā)的極限深度為1.2m，且出現(xiàn)在裸土含水率最小的8月下半月。

圖7 裸土蒸發(fā)強度圖Fig.7 Topsoil average evaporation intensity

圖8 裸土蒸發(fā)強度占水面蒸發(fā)強度的比例Fig.8 Percentage of topsoil to water surface average evaporation intensity

2）在西北干旱內(nèi)陸盆地，裸土蒸發(fā)主要發(fā)生在4月下半月至9月下半月，其中4月下半月至5月裸土蒸發(fā)強度最大，且裸土蒸發(fā)強度主要受外界蒸發(fā)能力的限制；6月至9月下半月裸土蒸發(fā)強度逐漸減小，且裸土蒸發(fā)強度主要受包氣帶向上輸水能力制約。

3）裸土蒸發(fā)強度平均占水面蒸發(fā)強度的3.5%，最大達6.7%。

（References）：

［1］束龍倉，荊艷東，黃修東，等.改進的無作物潛水蒸發(fā)經(jīng)驗公式［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2012，42（6）：1859-1865.

Su Longcang，Jing Yandong，Huang Xiudong，et al.On the Development of Improve Empirical Formulas for Calculating the Phreatic Water Evaporation for Bare Land［J］.Journal of Jilin University：Earth Science Edition，2012，42（6）：1859-1865.

［2］霍思遠，靳孟貴，梁杏.包氣帶弱滲透性黏土透鏡體對降雨人滲補給影響的數(shù)值模擬［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2013，43（5）：1579-1587.

Huo Siyuan，Jin Menggui，Liang Xing.Impacts of Low-Permeability Clay Lens in Vadose Zone onto Rainfall Infiltration and Groundwater Recharge Using Numerical Simulation of Variably Saturated Flow［J］.Journal of Jilin University ：Earth Science Edition，2013，43（5）：1579-1587.

［3］高曉飛，王曉嵐.微型蒸發(fā)器口徑影響土壤蒸發(fā)測量值的試驗研究［J］.灌溉排水學(xué)報，2011，30（1）：1-4.

Gao Xiaofei， Wang Xiaolan. Effect of Micro-Lysimeters on Soil Evaporation Measurement［J］.Journal of Irrigation and Drainage，2011，30（1）：1-4.

［4］王文科，楊澤元，程東會，等.面向生態(tài)的干旱半干旱地區(qū)區(qū)域地下水資源評價的方法體系［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2011，41（1）：159-167.

Wang Wenke，Yang Zeyuan，Cheng Donghui，et al.Method of Ecology-Oriented Groundwater Resource Assessment in Arid and Semi-Arid Area［J］.Journal of Jilin University：Earth Science Edition，2011，41（1）：159-167.

［5］胡順軍，雷加強，徐新文，等.塔克拉瑪干沙漠風(fēng)沙土潛水極限蒸發(fā)強度的理論分析［J］.科學(xué)通報，2008，53（增刊2）：107-111.

Hu Shunjun，Lei Jiaqiang，Xu Xinwen，et al.Theory Analysis of the Limit Evaporates Strength in Taklimakan Desert Aeolian Sandy Soil［J］.Sciences Bulletin，2008，53（Sup.2）：107-111.

［6］胡光成.銀川平原地表蒸發(fā)量的估算及其在生態(tài)水文地質(zhì)中的應(yīng)用［D］.北京：中國地質(zhì)大學(xué)，2010.

Hu Guangcheng.Estimation of Land Surface Evapotranspiration in Yinchuan Plain and Its Application in Eco-Hyduogeology［D］.Beijing：China University of Geosciences，2010.

［7］孫艷，王益權(quán)，劉軍.緊實脅迫條件下土壤水分的蒸發(fā)及數(shù)學(xué)模擬［J］.西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版，2009，37（11）：219-223.

Sun Yan， Wang Yiquan，Liu Jun.Soil Water Evaporation and Its Maths Simulation Under Compaction Stress Condition［J］.Journal of Northwest A ＆F University：Natural Science Edition，2009，37（11）：219-223.

［8］鄧銘江，王世江，董新光，等.新疆水資源及可持續(xù)利用［M］.北京：中國水利水電出版社，2005.

Deng Mingjiang，Wang Shijiang，Dong Xinguang，et al.Water Resources in Xinjiang and It’s Sustainable Utilization［M］.Beijing：China Waterpower Press，2005.

［9］朱學(xué)愚，錢孝星.地下水水文學(xué)［M］.北京：中國環(huán)境科學(xué)出版社，2005.

Zhu Xueyu，Qian Xiaoxing.Groundwater Hydrology［M］.Beijing：China Environmental Science Press，2005.

［10］李俊亭.地下水流數(shù)值模擬［M］.北京：地質(zhì)出版社，1988.

Li Junting.Groundwater Numerical Simulation［M］.Beijing：Geological Publishing House，1988.

［11］Van Genuchten M Th.A Closed-Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soil［J］.Soil Sci Soc Amer J，1980，44：892-898.

Evaporation Intensity of Bare Soil in Northwest Arid Inland Basin

Qiao Gang1，2，Wang Wenke1

1.SchoolofEnvironmentScienceandEngineering，Chang’anUniversity，Xi’an710054，China
2.Xi’anCenterofGeologicalSurvey，CGS，Xi’an710054，China

Owing to the vadose zone evaporation water from the ground surface，finding out the bare soil evaporation rule has important meaning to study some vadose zone water and underwater evaporation problems.To find out the rule of bare soil evaporation intensity in northwest arid inland basin，we analyze the vadose zone cross section water flow features with weather data and calculate the bare soil evaporation intensity adopting the finite-difference method through developing some bare soil water content and vadose zone pressure monitoring site experiments.Express as a result：According to vertical transportation different status of some vadose zone water，the cross section in the bare soil area vadose zone is divided into the vadose zone water upward to transfer zone and the vadose zone water downward to move zone.The extreme limit evaporation depth is 1.2min the bare soil zone of Northwest arid inland basin.Subjecting to the external space weather evaporation ability and the vadose zone water upward to transfer ability，the bare soil evaporation inside the basin appears in the second half month of April to the second half month of September，among which the bare soil evaporation quantity in the second half month of April to May is the biggest，and the evaporation intensity reaches to 0.56mm/d on the average.The percentage of the bare soil evaporation intensity only equivalent to surface evaporation intensity on the average of 3.5%，as biggest as 6.7%.

bare soil；evaporation intensity；total water potential；vadose zone；arid inland basin

10.13278/j.cnki.jjuese.201404211

P641.131

A

喬岡，王文科.西北干旱內(nèi)陸盆地區(qū)裸土蒸發(fā)強度.吉林大學(xué)學(xué)報：地球科學(xué)版，2014，44（4）：1327-1332.

10.13278/j.cnki.jjuese.201404211.

Qiao Gang，Wang Wenke.Evaporation Intensity of Bare Soil in Northwest Arid Inland Basin.Journal of Jilin University：Earth Science Edition，2014，44（4）：1327-1332.doi：10.13278/j.cnki.jjuese.201404211.

2013-10-22

教育部科學(xué)技術(shù)研究重大項目（308021）；教育部科學(xué)技術(shù)研究重點項目（02078）

喬岡（1980—，男，博士研究生，助理研究員，主要從事干旱區(qū)水資源合理開發(fā)利用與生態(tài)環(huán)境保護方面的研究，E-mail：qgcyboy＠163.com。