王志濤，繳錫云，韓紅亮，虞曉彬

(1.河海大學水文水資源與水利工程科學國家重點實驗室，江蘇南京 210098;2.河海大學水利水電學院，江蘇南京 210098)

隨著計算機技術(shù)的發(fā)展，采用數(shù)值計算方法來模擬土壤水分以及溶質(zhì)運移越來越受到人們的重視，并逐漸成為研究土壤水分運動及溶質(zhì)運移的主要方法［1］。然而，在進行模擬之前，必須先確定模型中的相關(guān)參數(shù)，土壤水分特征曲線就是其中必不可少的參數(shù)之一［2-4］，其準確程度直接影響到模型的精度。描述土壤水分特征曲線的模型很多［5-6］，應用較多的是Van Genuchten模型［7］(簡稱VG模型)。

在應用模型進行數(shù)值模擬時，測定、分析模型中諸多參數(shù)的變化對模擬結(jié)果的影響程度，判斷各個參數(shù)對模擬結(jié)果的重要性，這就是參數(shù)的敏感性分析。目前，模型參數(shù)的敏感性分析大多是針對水文模型而言的［8-9］，對土壤水分運動模型參數(shù)敏感性的研究較少。聶衛(wèi)波等［10-11］對一維垂直入滲過程的數(shù)值模擬進行了驗證。畢經(jīng)偉等［12］、孟麗江［13］以底部通量作為敏感性分析指標，得出飽和水力傳導度是最敏感的參數(shù)、飽和含水率次之。梁冰等［14］分析了VG模型中的α和n這2個經(jīng)驗參數(shù)對土壤重力釋水和入滲過程的影響，認為α較小時不易于土壤釋水過程的進行，但有利于水分入滲;n較小時不利于重力釋水的進行，對于入滲過程的影響則較為復雜。

本文以HYDRUS-1D模型［15］模擬土壤垂直一維入滲，采用單因素局部擾動分析方法，分別對VG模型中各個參數(shù)進行擾動(擾動幅度為±20%)，研究參數(shù)擾動后對濕潤鋒運移距離、累計入滲量的影響，旨在找到對土壤入滲影響顯著的VG模型參數(shù)，以便在實際應用時準確地確定敏感性較強的參數(shù)，避免在模擬過程中由于參數(shù)不精確所帶來的誤差。

1 入滲模型

1.1 VG模型及試驗土壤參數(shù)

VG 模型表達式［7］如下:

式中:h——壓力水頭，cm;θr——殘余含水率，cm3/cm3;θs——飽和含水率，cm3/cm3;α——土壤特征曲線的經(jīng)驗參數(shù)，cm-1;n——土壤特征曲線的經(jīng)驗參數(shù)，無量綱;Ks——飽和導水率，cm/min;K(h)——非飽和導水率，cm/min;m——經(jīng)驗系數(shù)，m=1-1/n(n ＞1)。

試驗用土取自河北省滄州市中國科學院南皮生態(tài)農(nóng)業(yè)試驗站。利用SEDIMA4-12粒度分析儀(德國)進行顆粒分析，得出該土樣為粉壤土(美國制分類標準)，密度為1.40 g/cm3，其砂粒、粉粒、黏粒的質(zhì)量分數(shù)分別為11.84%，66.21%，21.95%。根據(jù)土壤的機械組成及土壤密度，應用HYDRUS-1D模型中內(nèi)置的RETC模塊預測 VG 模型參數(shù)［15-17］，得出 θr=0.0743cm3/cm3，θs=0.4312cm3/cm3，α =0.0056cm-1，n=1.6345，Ks=0.0107 cm/min。

1.2 邊界條件及初始條件

在HYDRUS-1D模型中，土壤質(zhì)地、密度是通過VG模型中的參數(shù)取值體現(xiàn)的。土壤質(zhì)地、密度對土壤入滲過程的入滲速率、累計入滲量、土壤濕潤鋒運移距離有較大影響;入滲水頭對入滲速率、累計入滲量及土壤濕潤鋒運移距離影響微弱;初始含水率對入滲速率、累計入滲量的影響也較小，但是對濕潤鋒運移距離有較明顯的影響［10-11］。初始條件和邊界條件不會影響VG模型參數(shù)的敏感性，故在模擬時土壤初始體積含水率取為恒定值(0.15 cm3/cm3);上邊界取為恒定的壓力水頭(h=2.5 cm);下邊界取為自由排水界面。模擬土柱高度為100 cm，模擬時間為750 min。入滲過程中，濕潤鋒始終未到達下邊界，可視為半無限土體，因此下邊界壓力水頭為一恒定壓力水頭。上、下邊界條件分別為

2 入滲過程對VG模型參數(shù)的敏感性分析

VG模型參數(shù)θr，θs，α，n等擾動顯著地影響土壤水分特征曲線的形狀。θr和α的擾動均對VG模型曲線的高吸力段(低含水率)影響顯著，在中、低吸力段(中、高含水率)擾動曲線逐漸與標準曲線趨于一致;θs的擾動對VG模型曲線高、中、低吸力段均有影響;n的擾動對VG模型曲線的影響最顯著，且擾動曲線在飽和含水率處與基準曲線相交。Ks的擾動雖然不影響土壤水分特征曲線的形狀，但是Ks決定著土壤導水能力的強弱，對土壤水分運動的影響不容忽視。文獻［12-13］指出，Ks和θs是影響底部通量的關(guān)鍵性因素。

2.1 垂直濕潤距離計算值

θr，θs和α的擾動與垂直濕潤距離計算值Rz呈負相關(guān)關(guān)系，n和Ks與Rz呈正相關(guān)關(guān)系。n的擾動對Rz的影響最大，尤其是當n減小時，Rz明顯減少，這是由于n減小后相同體積含水率下的土壤吸力明顯增大，入滲速率減小，土壤水分運動減慢。此外，θs減小時入滲速率下降，但是土壤可含蓄的水量減小，因此Rz增加。VG模型各參數(shù)擾動對Rz的影響見圖1。

圖1 VG模型參數(shù)擾動對R z的影響Fig.1 Influence of perturbation of VG model parameters on R z

Rz與入滲歷時t可采用單參數(shù)冪函數(shù)擬合，對于粉壤土，采用式(4)擬合［11］:

式中m為擬合參數(shù)。式(4)擬合的結(jié)果見表1。

表1 VG模型參數(shù)擾動對R z，m的影響Table 1 Influence of perturbation of VG model parameters on R z and m

θs，n和Ks的擾動對Rz及m的影響較大，且Rz和m的擾動變幅均超過了10.00%;α的擾動影響較小，Rz和m的擾動變幅最大值分別為8.11%和8.40%;θr擾動影響最小，Rz和m的擾動變幅不足3.00%;n的負擾動影響最大，Rz和m的擾動變幅分別為-31.76%，-33.19%。

2.2 累計入滲量計算值

θs，n和Ks的擾動均與累計入滲量計算值I呈正相關(guān)關(guān)系，且n的擾動對I的影響最大;α的擾動與I呈負相關(guān)關(guān)系;θr的擾動基本不影響I。各參數(shù)的擾動對I的影響見圖2。

為了定量分析VG模型各參數(shù)的擾動對I影響，采用Philip入滲公式進行擬合:

式中S為滲吸系數(shù)，cm/min0.5。VG模型各參數(shù)擾動對I及Philip公式中滲吸系數(shù)的影響結(jié)果見表2。

圖2 VG模型參數(shù)擾動對累計入滲量計算值的影響Fig.2 Influence of perturbation of VG model parameters on calculated values of cumulative infiltration

表2 VG模型參數(shù)擾動對I及S的影響Table 2 Influence of perturbation of VG model parameters on I and S

結(jié)果表明，θs，n和Ks的擾動對I及Philip公式中S的影響較大，I和S的變幅均在10.00%以上;α的擾動對I及S的影響較小，I和S的變幅在-6.69% ～9.16%之間;θr的擾動影響最小，I和S的變幅均低于0.10%，可忽略不計。由模擬結(jié)束時刻(750 min)的I可知，n的正、負擾動對其影響最大，I的變幅分別為17.63%和-26.33%;其次是θs，在其擾動下I的變幅分別為11.56%和-13.32%;Ks擾動下，I的變幅為11.47%和-12.77%。VG模型各參數(shù)n，θs，Ks，α，θr的擾動對S的影響依次減弱。n的負擾動對S的影響最大，S的變幅為-27.97%。

3 結(jié) 論

a.n，θs和Ks的擾動對Rz和I的影響均較大，模擬結(jié)束時Rz和I的擾動影響均不低于10%，且對Rz，I模型參數(shù)的擾動影響也超過了10%;α的擾動影響次之;θr的擾動影響最小。

b.n值的負擾動對Rz和I的影響最大，且n偏小時對模擬結(jié)果帶來的誤差較大。因此，在實際應用中應保證VG模型參數(shù)n，θs，Ks等的精度。

［1］雷志棟，楊師秀，謝森傳.土壤水動力學［M］.北京:清華大學出版社，1988.

［2］徐紹輝，劉建立.土壤水力性質(zhì)確定方法研究進展［J］.水科學進展，2003，14(4):394-401.(XU Shaohui，LIU Jianli.Advances in approaches for determining unsaturated soil hydraulic properties［J］.Advances in Water Science，2003，14(4):394-401.(in Chinese))

［3］劉建立，徐紹輝，劉慧.估計土壤水分特征曲線的間接方法研究進展［J］.水利學報，2004，35(2):68-76.(LIU Jianli，XU Shaohui，LIU Hui.A review of development in estimating soil water retention characteristics from soil data ［J］.Journal of Hydraulic Engineering，2004，35(2):68-76.(in Chinese))

［4］李峰，繳錫云，李盼盼，等.田間土壤水分特征曲線參數(shù)反演［J］.河海大學學報:自然科學版，2009，37(4):373-377.(LI Feng，JIAO Xiyun，LI Panpan，et al.Parametric inversion of soil water characteristic curves of farmland［J］.Journal of Hohai University:Natural Sciences，2009，37(4):373-377.(in Chinese))

［5］呂殿青，邵明安.非飽和土壤水力參數(shù)的模型及確定方法［J］.應用生態(tài)學報，2004，15(1):163-166.(LYU Dianqing，SHAO Ming'an.Models and determining methods of unsaturated soil hydraulic parameters［J］.Chinese Journal of Applied Ecology，2004，15(1):163-166.(in Chinese))

［6］盧小慧，靳孟貴，汪丙國.欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗站土壤水分特征曲線分析［J］.中國農(nóng)村水利水電，2006(12):30-32.(LU Xiaohui，JIN Menggui，WANGBingguo.Discussion on the soil water characteristic curve of the agricultural eco-system experiment station in Luancheng，Hebei Province［J］.China Rural Water and Hydropower，2006(12):30-32.(in Chinese))

［7］GENUCHTEN M T.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils［J］.Soil Sci Soc Am J，1980，44(5):892-898.

［8］劉毅，陳吉寧，杜鵬飛.環(huán)境模型參數(shù)識別與不確定性分析［J］.環(huán)境科學，2002，23(6):6-10.(LIU Yi，CHEN Jining，DU Pengfei.Parameters identification and uncertainty analysis for environmental model［J］.Environmental Science，2002，23(6):6-10.(in Chinese))

［9］田雨，雷曉輝，蔣云鐘，等.水文模型參數(shù)敏感性分析方法研究評述［J］.水文，2010，30(4):9-12.(TIAN Yu，LEI Xiaohui，JIANG Yunzhong，et al.Comment on parameter sensitivity analysis of hydrological model［J］.Journal of China Hydrology，2010，30(4):9-12.(in Chinese))

［10］聶衛(wèi)波.畦溝灌溉水流運動模型與數(shù)值模擬研究［D］.楊凌:西北農(nóng)林科技大學，2009.

［11］聶衛(wèi)波，馬孝義，王術(shù)禮.一維土壤入滲特性數(shù)值模擬［J］.灌溉排水學報，2009，28(3):53-57.(NIE Weibo，MA Xiaoyi，WANG Shuli.Numerical simulation of one-dimensional soil infiltration characteristic ［J］.Journal of Irrigation and Drainage，2009，28(3):53-57.(in Chinese))

［12］畢經(jīng)偉，張佳寶，陳效民，等.農(nóng)田土壤中土壤水滲漏與硝態(tài)氮淋失的模擬研究［J］.灌溉排水學報，2003，22(6):23-26.(BI Jingwei，ZHANG Jiabao，CHEN Xiaomin，et al.Modeling of soil water erainage and nitrate leaching in arable soils［J］.Journal of Irrigation and Drainage，2003，22(6):23-26.(in Chinese))

［13］孟江麗.水鹽和排水系統(tǒng)的模擬研究［D］.烏魯木齊:新疆農(nóng)業(yè)大學，2004.

［14］梁冰，姜利國.非飽和帶水分運動模型的Picard解法及參數(shù)靈敏度分析［J］.水資源與水工程學報，2009，20(5):23-27.(LIANG Bing，JIANG Liguo.Picard solution of moisture movement model in unsaturated zone and parameter sensitivity analysis［J］.Journal of Water Resources＆ Water Engineering，2009，20(5):23-27.(in Chinese))

［15］?IMUNEK J，?EJNA M，SATIO H，et al.The HYDRUS-1D software package for simulating the one-dimensional movement of water，heat，and multiple solutes in variably-saturated media［M］.Version 4.California:U.S.Salinity Laboratory，2008.

［16］SCHAAP M G，LEIJ F J.Using neural networks to predict soil water retention and soil hydraulic conductivity［J］.Soil and Tillage Research，1998，47(1/2):37-42.

［17］魏義長，劉作新，康玲玲.遼西淋溶褐土土壤水動力學參數(shù)的推導及驗證［J］.水利學報，2004，35(3):81-86.(WEI Yichang，LIU Zuoxin，KANG Lingling.Derivation and verification of soil hydrodynamic parameters in cinnamon soil［J］.Journal of Hydraulic Engineering，2004，35(3):81-86.(in Chinese))