栗現(xiàn)文,周金龍,2

(1.中國地質(zhì)大學(xué)環(huán)境學(xué)院,湖北武漢 430074;2.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆烏魯木齊 830052)

1 高鹽度潛水蒸發(fā)實驗的意義

研究潛水蒸發(fā)機理及其計算方法在干旱區(qū)水資源評價、鹽堿土治理、灌區(qū)改造及生態(tài)需水量計算方面都具有十分重要的意義。在極端干旱區(qū)潛水礦化度常常達到 n×10 g/L,甚至高達 n×100 g/L,在現(xiàn)有的研究中,很少關(guān)注高鹽度潛水蒸發(fā)的機理及高鹽度潛水蒸發(fā)量的精細計算問題。與傳統(tǒng)意義的潛水蒸發(fā)相比,極端干旱區(qū)高鹽度潛水蒸發(fā)條件下,毛管水上升高度、鹽漬臨界深度、土壤水總水勢組成、非飽和導(dǎo)水率與土壤水分特征曲線都將存在一定的差異,將淡水條件下得出的潛水蒸發(fā)有關(guān)參數(shù)直接用于極端干旱區(qū)高鹽度潛水蒸發(fā)的計算或模擬以及調(diào)控鹽荒地水鹽狀況,必將產(chǎn)生偏差。因此,通過蒸滲儀潛水蒸發(fā)模擬實驗及原位鹽荒地潛水蒸發(fā)模擬計算,深入研究極端干旱區(qū)高鹽度潛水蒸發(fā)機理及計算方法,必將有助于加深極端干旱區(qū)潛水蒸發(fā)機理的認識、合理計算鹽生植被生態(tài)需水量、指導(dǎo)鹽荒地開發(fā)與水鹽調(diào)控。

2 國內(nèi)外研究現(xiàn)狀及存在問題

多年來包括土壤物理、水文地質(zhì)、農(nóng)田水利、水文學(xué)及水資源等學(xué)科在內(nèi)的許多專家學(xué)者對干旱區(qū)潛水蒸發(fā)的機理與規(guī)律做了大量的研究工作。國內(nèi) 20世紀 50年代初,為了防治土壤鹽漬化,潛水蒸發(fā)的研究成為當時的熱點。20世紀70年代初,大量開發(fā)利用地下水,為了正確評價地下水資源,潛水蒸發(fā)的研究有了新的進展。由于近年來土壤水動力學(xué)和 SPAC系統(tǒng)研究的不斷深入和發(fā)展,潛水蒸發(fā)研究在以下6個方面都取得了明顯的進展。

在試驗研究的手段方面,負壓計[1]、各種土壤水分測定儀[2]、蒸發(fā)儀[3]等得到廣泛應(yīng)用;吳運卿等(2006)研制了智能化稱重蒸滲儀系統(tǒng)[4];meissner等(2008)研制了室外使用的、在價格上低于傳統(tǒng)鋼質(zhì)或混凝土質(zhì)蒸發(fā)儀的集裝箱式的稱重蒸滲儀[5];Unold等(2008)研究了適用于專門應(yīng)用需要的田間稱重蒸滲儀的模塊化設(shè)計[6]。

在毛管水上升高度分析與測定方面,尹娟(2007)[7]、史文娟等(2007)[8]開展了實驗和研究;Stenitzer等(2007)通過估算土壤轉(zhuǎn)移函數(shù)或測量水力參數(shù),應(yīng)用 SImWASER模擬模型評價潛水毛細上升高度[9]。

在干旱區(qū)潛水蒸發(fā)的影響因素(包氣帶巖性及組合特征、潛水埋深、氣象因素、地溫、植被覆蓋情況等)方面,周金龍等(2002)[10]、Yang等(2002)[11]、周金龍等 (2003)[12]、史文娟等(2006)[13]和 Zhou等(2008)[14]開展了實驗與研究。

在潛水蒸發(fā)量計算經(jīng)驗公式的構(gòu)建和適用性研究方面,國內(nèi)很多學(xué)者[15-17]和少數(shù)國外學(xué)者[18]從不同角度開展了研究。

在潛水蒸發(fā)的數(shù)值模擬模型方面,不少國外學(xué)者[19-20]和國內(nèi)學(xué)者[21]在模型研制和應(yīng)用方面做了大量的工作。值得指出的是,Grunberger等(2008)[20]通過對鹽土土壤剖面氘和 Cl-含量分析,確定土壤蒸發(fā)量。

在潛水蒸發(fā)研究成果的應(yīng)用方面,分別在作物的地下水利用量[22-24]、自然植被生態(tài)需水量計算[25-27]、生態(tài)水位的確定[28]、鹽堿地形成機理分析及其水鹽調(diào)控[29]和巖土工程與環(huán)境工程[30]等方面開展了研究工作。

在上述研究中,均未考慮潛水礦化度對潛水蒸發(fā)的影響,可以視為傳統(tǒng)意義上的潛水蒸發(fā)研究。潛水礦化度(尤其是高鹽度)對土壤水礦化度、土壤總水勢、土壤非飽和滲透率、土壤層地溫、毛管水上升高度、土壤水分測定儀的土中標定方程都可能產(chǎn)生不同程度的影響,從而影響到潛水蒸發(fā)量的大小。潛水蒸發(fā)強烈的地區(qū)主要發(fā)生在干旱區(qū),在極端干旱區(qū)潛水蒸發(fā)強烈的區(qū)域,潛水埋深一般小于 2m,包氣帶巖性一般為粉土、亞黏土,潛水礦化度一般為 3～10g/L,局部為10～30g/L,最大可達 100 g/L以上,高鹽度潛水的水化學(xué)類型一般為 Cl?SO4-Na型。如新疆沖積平原或湖積平原三角洲地帶潛水礦化度由 ng/L增至 n×10 g/L,甚至 300 g/L,高鹽度潛水水化學(xué)類型為 Cl?SO4型或 Cl型[31];在河西走廊、準格爾盆地、塔里木盆地、吐魯番盆地和柴達木盆地均分布有鹽化和濃縮作用的鹽鹵水,其中礦化度一般為 3～10 g/L,有時大于 50 g/L,水化學(xué)類型為以 Cl-Na型為主,河西走廊為 100 g/L,特別是在柴達木盆地潛水礦化度高達 300 g/L以上,除 Cl-Na型外,還出現(xiàn) Cl-mg型鹵水[32]。

我國潛水蒸發(fā)強烈的區(qū)域主要在西北干旱鹽漬土區(qū),包括新疆、青海、甘肅河西走廊和內(nèi)蒙古西部大部分地區(qū),為我國鹽漬土分布最廣的地區(qū),鹽漬化面積約為7000×104hm2。除少數(shù)有灌溉條件的地區(qū)人工栽培各種植被代替了自然鹽生植被,其余絕大部分是尚未開墾的鹽荒地,仍然保持著鹽生自然植被。鹽生植被是鹽漬化地區(qū)植被的主體,應(yīng)當對鹽生植被覆蓋區(qū)的水鹽動態(tài)、潛水蒸發(fā)規(guī)律進行研究[33]。用在淡水條件下完成的潛水蒸發(fā)試驗以及數(shù)值模擬成果計算極端干旱區(qū)高鹽度潛水蒸發(fā)量與鹽生植被生態(tài)需水量以及指導(dǎo)鹽荒地開發(fā)與水鹽調(diào)控,必將產(chǎn)生偏差。

Soppe等(2003)[34]用稱重式蒸滲儀研究了紅花對鹽度為14ds/m的地下水的利用情況。Jorenush等(2003)[35]應(yīng)用改進的 TSA m模型模擬研究了灌溉和非灌溉條件下,不同潛水位埋深(0.3～1.20m)和不同鹽度(0.5～13.0 dS/m)條件下的毛管水上升高度和土壤鹽度的變化規(guī)律。Weisbrod等(2004)[36]對5mol/LNaNO3溶液在變飽和孔隙介質(zhì)中的遷移進行了實驗研究。Northey等(2006)[37]應(yīng)用蒸滲儀定量化研究了潛水位和鹽度存在短期變化條件下淺層地下水的毛管上升水流和潛在鹽化風險。Tikhonravova(2007)[38]研究了壤土的熱擴散率與土壤含水量和含鹽量的關(guān)系。劉廣明等(2003)[39]建立了地下水累積蒸發(fā)量與地下水埋深、地下水礦化度的關(guān)系。秦耀東(2003)[40]認為在研究植物吸收土壤水和土壤蒸發(fā)時需要考慮溶質(zhì)勢的作用,并介紹了依據(jù)溶液質(zhì)量濃度和溶質(zhì)的摩爾質(zhì)量計算溶質(zhì)勢的 Van′t Hoff方程。胡守忠等(2006)[41]開展了鹽漬化地區(qū)土壤水分運動參數(shù)的測定及標定的研究工作,建立了土壤基質(zhì)勢和飽和度之間的數(shù)學(xué)模型。在上述研究中,有關(guān)學(xué)者已經(jīng)注意到了礦化度對潛水蒸發(fā)的影響,但僅限于灌區(qū)潛水礦化度的范圍(一般小于10g/L),未涉及潛水礦化度 10～300 g/L的情況,且以室內(nèi)模擬實驗為主。因此,上述研究成果仍無法揭示極端干旱區(qū)潛水礦化度一般為 10～300 g/L的鹽荒地和鹽生植被覆蓋區(qū)的潛水蒸發(fā)和水鹽運移規(guī)律。因未進行淡水(小于1 g/L)、咸水(1～10 g/L)、鹽水(10～50 g/L)和鹵水(50～300 g/L)蒸發(fā)的系統(tǒng)對比試驗,因此,無法回答諸如礦化度對毛管水上升高度及潛水蒸發(fā)的影響程度如何?高鹽度潛水蒸發(fā)量與潛水埋深、包氣帶巖性、潛水礦化度間的關(guān)系如何?高鹽度潛水分布區(qū)鹽生植被生態(tài)需水量如何計算?如何在高鹽度條件下進行潛水 -土壤 -植被 -大氣連續(xù)體水鹽運移的數(shù)值模擬?等問題。上述問題的解決,對于提高極端干旱區(qū)潛水蒸發(fā)量、鹽荒地天然植被生態(tài)需水量的計算精度以及指導(dǎo)宜農(nóng)鹽荒地開發(fā)與水鹽調(diào)控具有重要的應(yīng)用價值,也有利于深化極端干旱區(qū)潛水蒸發(fā)機理的認識,即具有重要的科學(xué)意義。

3 值得進一步研究的問題及擬采取的技術(shù)路線

分析國內(nèi)外研究現(xiàn)狀,極端干旱區(qū)高鹽度(10～300 g/L)潛水蒸發(fā)的實驗研究幾乎是空白,機理研究尚處于起步階段,有以下幾個方面的問題值得深入研究。

3.1 礦化度對毛細水上升高度及潛水蒸發(fā)影響的研究

可利用小型蒸滲儀開展不同土質(zhì)、不同礦化度潛水條件下毛管水上升高度的測定工作。建立不同土質(zhì)毛管水上升高度與礦化度間的統(tǒng)計模型,為高鹽度地區(qū)水鹽調(diào)控和生態(tài)水位的合理確定提供依據(jù);開展不同土質(zhì)、不同礦化度潛水蒸發(fā)的對比實驗,掌握礦化度對潛水蒸發(fā)的影響程度,建立不同土質(zhì)潛水蒸發(fā)量與潛水位埋深、潛水礦化度間的統(tǒng)計模型。

3.2 高鹽度潛水蒸發(fā)條件下的土壤水動力學(xué)參數(shù)率定(測定)和數(shù)值模擬研究

可利用大型蒸滲儀開展高鹽度潛水蒸發(fā)條件下的土壤總水勢、土壤水分特征曲線和中子儀標定的試驗(測定)研究工作,掌握潛水礦化度對土壤總水勢、土壤水分特征曲線和中子儀標定方程的可能影響及其機理。在上述研究工作基礎(chǔ)上,應(yīng)用 HYDRUS-1D模型等開展不同土質(zhì)、不同礦化度條件下非飽和水動力參數(shù)的率定和不同土質(zhì)(組合)、不同礦化度條件下潛水蒸發(fā)的數(shù)值模擬計算,并與實測潛水蒸發(fā)量進行對比。

3.3 野外原位鹽荒地潛水-土壤水水 -鹽運移規(guī)律與鹽荒地水鹽調(diào)控研究

可選擇典型的高鹽度潛水地段,測定模型所需要的鹽運移參數(shù),利用 3.2率定的非飽和水動力參數(shù),建立野外原位鹽荒地潛水 -土壤水水、鹽運移數(shù)值模擬模型,借助數(shù)值模擬手段研究潛水礦化度對鹽荒地植物蒸騰量、土壤剖面積鹽狀況的影響,為提出鹽荒地水鹽調(diào)控方案提供技術(shù)支撐。

極端干旱區(qū)高鹽度(10～300 g/L)潛水蒸發(fā)實驗研究擬采取的技術(shù)路線如下圖 1所示。

圖 1高鹽度(10～300 g/L)潛水蒸發(fā)實驗研究技術(shù)路線框圖

[1]孔凡哲,王曉贊.利用土壤水吸力計算潛水蒸發(fā)初探[J].水文,1997,(3):44-47.

[2]高照陽,張紅梅,常明勛,等.國內(nèi)外土壤水分監(jiān)測技術(shù)[J].節(jié)水灌溉,2004,(2):28-29.

[3]Chuck Young,W esley Wallender,Gerrit Schoups,et al.modeling shallow water table evaporation in irrigated regions[J].Irrig Drainage Syst,2007,(21):119-132.

[4]吳運卿,羅金耀,王富慶.智能化稱重式蒸滲儀系統(tǒng)的研制與實現(xiàn)[J].實驗室研究與探索,2006,25(4):432-434,438.

[5]Ralph meissner,Holger Rupp,manfred Seyfarth.Advances in out door lysi meter techniques[J].Water Air Soil Pollut:Focus(2008)8:217-225.

[6]G.von Unold,J.Fank.modular design of field lysimeters for specific application needs[J].Water Air Soil Pollut:Focus(2008)8:233-242.

[7]尹娟,費良軍,程東娟.均質(zhì)土壤毛管水上升特性室內(nèi)試驗研究[J].農(nóng)業(yè)工程學(xué)報,2007,23(6):91-94.

[8]史文娟,沈 冰,汪志榮,等.層狀土壤毛管水最大上升高度分析[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(1):94-97.

[9]E.Stenitzer,H.Diestel,Th.Zenker,etal.Assessmentof capillary rise fro m shallow groundwater by the simulationmodel SImWASER using either estimated pedotransfer functions ormeasured hydraulic para meters[J].Water Resour manage,2007,21:1567-1584.

[10]周金龍,虎膽?吐馬爾白,董新光,等.新疆平原區(qū)大氣降水 -灌溉水 -土壤水 -地下水水量轉(zhuǎn)化關(guān)系實驗研究[M].烏魯木齊:新疆科技出版社,2002,57-79.

[11]mingdi Yang,Ernest K.Yanful.Water balance during evaporation and drainage in cover soils under differentwater table conditions[J].Advances in Environ mental Research,2002,6:505-521.

[12]周金龍,董新光,王斌.新疆平原區(qū)潛水蒸發(fā)研究[J].工程勘察,2003,(5):23-27.

[13]史文娟,沈冰,汪志榮,等.層狀土壤中砂層層位對潛水蒸發(fā)的影響[J].干旱區(qū)地理,2006,(2):282-286.

[14]ZHOU Hong-W ei,LISheng-Yu,SUN Shu-guo,et al.Effects of natural coverson soil evaporation of the shelterbelt along the Tari m Desert Highway[J].Chinese Science Bulletin,2008,53(Supp.II):137-145.

[15]HU Shun-jun,TIAN Chang-yan,SONG Yu-dong.model for calculating phreatic water evaporation on bare and Tamarix-vegetated lands[J].Chinese Science Bulletin,2006,51(Supp.Ⅰ):43-50.

[16]J.H.Si,Q.Feng,X.Y.Zhang,etal.Growing season evapotranspiration from Tamarix ramosissimastands under extremearid conditions in northwest China[J].Environ Geol(2005)48:861-870.

[17]劉淼,郝振純,王加虎,等.淮北平原裸地潛水蒸發(fā)計算公式比較[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2008,(38)4:649-653.

[18]Laura K.Lautz.Estimating groundwater evapotranspiration ratesusing diurnal water-table fluctuations in a sem i-arid riparian zone[J].Hydrogeology Journal,2008,16:483-497.

[19]Chuck Young,Wesley Wallender,Gerrit Schoups,et al.modeling shallow water table evaporation in irrigated regions[J].Irrig D rainage Syst,2007,21:119-132.

[20]Olivier Grunberger,Peggy macaigne,Jean-Luc m ichelot,et al.Salt crust development in paddy fields owing to soil evaporation and drainage:Contribution of chloride and deuteriu m profile analysis[J].Journal of Hydrology,2008,(348):110-123.

[21]雷志棟,尚松浩,楊詩秀,等.土壤凍結(jié)過程中潛水蒸發(fā)規(guī)律的模擬研究[J].水利學(xué)報,1999,(6):6-10.

[22]A.R.Sepaskhah,A.Kanooni,m.Ghasem i.Esti mating water table contribution to corn and sorghu m water use[J].Agricultural water manage ment,2003,58:70-79.

[23]Caecelia A.Hurst,Peter J.Thorburn,David Lockington,et al.Sugarcane w ater use from shallow water tables:implications for i mproving irrigation water use efficiency[J].Agricultural water management,2004,65:1-19.

[24]王曉紅.潛水蒸發(fā)及作物的地下水利用量估算方法的研究[D].武漢:武漢大學(xué),2004.

[25]王根緒,程國棟.干旱內(nèi)陸流域生態(tài)需水量及其估算[J].中國沙漠,2002,22(2):129-134.

[26]張麗,董增川,趙斌.干旱區(qū)天然植被生態(tài)需水量計算方法[J].水科學(xué)進展,2003,14(6):745-748.

[27]周金龍,劉豐,侯江濤,等.新疆平原非灌區(qū)自然植被生態(tài)耗水量的計算[J].干旱區(qū)資源與環(huán)境,2006,20(4):162-165.

[28]樊自立,陳亞寧,李和平,等.中國西北干旱區(qū)生態(tài)地下水埋深適宜深度的確定[J].干旱區(qū)資源與環(huán)境,2008,22(2):1-5.

[29]m.Ibrakhimov,A.Khamzina,I.Forkutsa,et al.Groundwater table and salinity:spatial and te mporal distribution and influence on soil salinization in Khorez m region(Uzbekistan,A ral Sea Basin)[J].Irrig Drainage Syst,2007,21:219-236.

[30]Yu-Jun Cui,Jorge G.Zornberg.W ater balance and evapotranspiration monitoring in geotechnical and geoenviron mental engineering[J].Geotechnical and Geological Engineering,2008,26(6):783-798.

[31]中國科學(xué)院新疆綜合考察隊,中國科學(xué)院地質(zhì)研究所,中國科學(xué)院新疆分院.新疆地下水[M].北京:科學(xué)出版社,1965,55-57.

[32]中國科學(xué)院《中國自然地理》編輯委員會.中國自然地理?地下水[M].北京:科學(xué)出版社,1981,69,72-75.

[33]徐恒剛主編.中國鹽生植被及鹽漬化生態(tài)治理[m].北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2004,1-36.

[34]R.W.O.Soppe,J.E.Ayars.Characterizing ground water use by safflower using weighing lysimeters[J].AgriculturalWatermanage ment,2003,60:59-71.

[35]m.H.Jorenush,A.R.Sepaskhah.modelling capillary rise and soil salinity for shallow saline water table under irrigated and nonirrigated conditions[J].Agricultural Water management,2003,61:125-141.

[36]Noam W eisb rod,m ichael R.Niemet,mark L.Rockhold,et al.m igration of saline solutions in variably saturated porousmedia[J].Journal of Conta m inant Hydrology,2004,72:109-133.

[37]J.E.Northey,E.W.Christen,J.E.Ayars,etal.Occurrence and measurement of salinity stratification in shallow groundwater in the murru mbidgee irrigation area,south-eastern Australia[J].Agricu lturalW atermanagement,2006,81:23-40.

[38]P.I.Tikhonravova.Effectof thewater contenton the thermal diffusivity of clay loams with different degrees of salinization in the Transvolga region[J].Eurasian Soil Science,2007,40(1):47-50.

[39]劉廣明,楊勁松.地下水作用條件下土壤積鹽規(guī)律研究[J].土壤學(xué)報,2003,(1):8-12.

[40]秦耀東.土壤物理學(xué)[M].北京:高等教育出版社,2003:44-47.

[41]胡守忠,喬冬梅,史海濱.鹽漬化地區(qū)土壤水分運動參數(shù)的測定及標定[J].灌溉排水學(xué)報,2006,(2):30-33,38.