劉莉莉，喬 偉，孫青言

(1.山東省水利科學(xué)研究院，山東 濟(jì)南 250013；2.濱州市水利局，山東 濱州 256603；3.中國水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室，北京 100038)

土壤水分入滲是陸面水文過程的重要組成部分，降水入滲過程及入滲量的變化直接影響地面的產(chǎn)流過程和產(chǎn)流量，因此入滲機(jī)理的理解與刻畫是研究流域(區(qū)域)水文循環(huán)的關(guān)鍵內(nèi)容之一[1]。水文模型是描述水文循環(huán)機(jī)理的重要工具，而降水(灌溉)入滲模擬則是水文模型的核心功能之一[2]。早期降水入滲數(shù)學(xué)模型通常是建立在均質(zhì)非凍土壤基質(zhì)之上，對模型的假設(shè)約束較多，限制了模型的應(yīng)用范圍。如Green-Ampt模型在寒區(qū)應(yīng)用時往往需要考慮土壤凍融循環(huán)造成的土壤水力條件變化[3]。寒區(qū)溫度低，土壤水易凝固，形成季節(jié)性凍土，甚至多年凍土。凍土的凍融不但會改變土壤的強(qiáng)度，導(dǎo)致構(gòu)筑其上的建筑設(shè)施發(fā)生變形、塌陷和破壞[4]，而且會改變土壤的物理性質(zhì)和水力學(xué)性質(zhì)，使以土壤為傳遞基質(zhì)的水分運(yùn)移發(fā)生變化[5]。其中，凍土直接改變了土壤的滲透性，使降水、融雪、灌溉等地表來水的入滲受到影響，進(jìn)而影響寒區(qū)的水文循環(huán)過程。因此，對寒區(qū)凍土凍融循環(huán)規(guī)律及其對入滲的影響機(jī)理進(jìn)行深入研究具有重要的科學(xué)意義。

中國寒區(qū)面積約占國土面積的43.5%左右[6]，包含了高山、高原、荒漠、草原、森林,以及農(nóng)田等多種景觀類型。這些地區(qū)受氣溫波動的影響，出現(xiàn)土壤層凍結(jié)與融化交替的現(xiàn)象，使水文過程中的入滲機(jī)制出現(xiàn)了不同于長年無凍土地區(qū)的特殊性。在利用水文模型研究寒區(qū)水文循環(huán)的過程中，如果忽略了凍土的影響，可能會使產(chǎn)流、入滲、地下水補(bǔ)給等水文循環(huán)過程的路徑和通量模擬產(chǎn)生較大誤差；在基于區(qū)域(流域)水文循環(huán)轉(zhuǎn)化的水資源評價中造成地表、地下水資源量核算誤差。如凍土層阻滯了融雪下滲補(bǔ)給地下水，使這部分水分以地表徑流和壤中流的形式匯入江河，形成地表水，增加了地表水資源量[7]。因此，建立具有寒區(qū)特點(diǎn)的土壤水分入滲理論和模型是研究寒區(qū)水文循環(huán)機(jī)制、完善寒區(qū)水文水資源理論、改進(jìn)相關(guān)應(yīng)用實(shí)踐的基礎(chǔ)。

對于寒區(qū)凍土凍融循環(huán)條件下的水分入滲機(jī)理研究，均是通過大量的室內(nèi)、室外科學(xué)實(shí)驗(yàn)，從把握土壤的凍融循環(huán)規(guī)律入手，揭示凍土在水分結(jié)晶和融化影響下的水力學(xué)性質(zhì)和結(jié)構(gòu)特性變化，在此基礎(chǔ)上改進(jìn)原有土壤入滲方程，或者創(chuàng)建新的適用于寒區(qū)凍土特點(diǎn)的土壤入滲方程，為更精確地刻畫寒區(qū)產(chǎn)流(入滲)這一關(guān)鍵水文環(huán)節(jié)提供理論、方法和工具上的支持。

1 寒區(qū)土壤水凍融機(jī)制

土壤水凍融機(jī)制是指土壤水在各種內(nèi)部和外部因素的影響下凍結(jié)和融化的動態(tài)特征，其中內(nèi)部和外部因素包括氣溫、土壤溫度、土壤含水量、土壤類型、積雪覆蓋、植物殘余覆蓋、植被類型，以及地下水等，正是由于這些因素與土壤水的相互作用機(jī)制，使得寒區(qū)凍土的凍融過程在不同的區(qū)域、時期和下墊面條件下出現(xiàn)了不同的動態(tài)特征。

溫度是寒區(qū)土壤水凍融的直接控制因素。氣溫與土壤溫度并不一致，以表層土壤溫度與氣溫的比較為例。氣溫在波動中逐漸降低，降到一定程度土壤開始從地表凍結(jié)，直到氣溫出現(xiàn)上升的趨勢，表層土壤開始融化。表層土壤從凍結(jié)到融化之間的時期內(nèi)，表層土壤溫度波動幅度明顯小于氣溫波動幅度[8]。土壤凍結(jié)的溫度因土壤類型的不同有所差異。細(xì)沙型土壤的凍結(jié)溫度在-0.26～0 ℃之間，粉土的凍結(jié)溫度在-0.42～0 ℃之間，黏土的凍結(jié)溫度則在-0.47～0 ℃之間，可見凍結(jié)溫度有隨土壤粒徑減小而降低的趨勢[9]。凍土的厚度、空間分布范圍均與溫度存在密切關(guān)系[10-12]。

土壤水分的凝固和融化是土壤凍融的本質(zhì)[7]。土壤凍融是指土壤中水分的凍結(jié)與融化，因此與土壤水含量密切相關(guān)。土壤凍融隨季節(jié)變化的范圍又稱為活動層，是土壤層物質(zhì)和能量交換最為活躍的區(qū)域，是凍土層最大季節(jié)融化厚度[13]。根據(jù)大量的實(shí)驗(yàn)研究，活動層厚度與土壤表層含水量存在較強(qiáng)的負(fù)線性相關(guān)性。另外，土壤含水量對凍土凍結(jié)鋒面的移動速率具有重要影響。吳禮舟等[14]根據(jù)凍結(jié)鋒面的移動規(guī)律，建立了凍脹量隨時間的關(guān)系式，并通過算例發(fā)現(xiàn)含水量和干密度對凍土溫度分布產(chǎn)生不小的影響，從而得出“含水量越大，凍結(jié)鋒面移動越快”的結(jié)論。

植物殘余對土表的覆蓋會影響土壤的凍結(jié)，實(shí)質(zhì)上是通過改變土壤的微氣候來改變土壤凍結(jié)時間、速度、深度等。以玉米收割后秸稈殘余的處理為例，秸稈殘莖站立、殘葉松散平鋪地面，對土壤具有很好的保溫作用，相對于裸土地，可有效減小凍土層的厚度，使凍土層提前融化[15]。地表覆蓋通過影響凍土的凍融狀態(tài)改變土壤的滲透性質(zhì)，對土壤墑情也起到一定的保持作用[16-18]。另外，與植物殘余對凍土的影響類似，積雪覆蓋也對土壤凍融過程產(chǎn)生影響。大量實(shí)驗(yàn)表明，地表積雪可以很好地保持土壤溫度，改變地表的熱量平衡模式[19]，從而減緩了土壤凍結(jié)速度，減小了土壤的凍深，增加了土壤水分向冷端的遷移量[20-21]。

淺層地下水(潛水)對上層包氣帶土壤具有補(bǔ)給作用，與包氣帶凍土層存在密切的相互作用關(guān)系。受土壤凍融的影響，不同地下水位對土壤水的補(bǔ)給強(qiáng)度不同。室內(nèi)試驗(yàn)研究發(fā)現(xiàn)，地下水埋深存在一個閾值，小于該閾值時，地下水對土壤水的累積補(bǔ)給量隨埋深增加而增加，大于該閾值時，累積補(bǔ)給量則隨埋深增加趨于減少[22]?？梢姡叵滤坏淖兓瘜ν寥篮康淖兓嬖谳^大影響，進(jìn)而改變土壤的凍融規(guī)律。雷志棟等[23]利用水熱耦合遷移模型研究了地下水淺埋深條件下凍土的凍融循環(huán)規(guī)律：凍土在快速凍結(jié)階段凍結(jié)速率隨凍結(jié)深度增加而減??；凍結(jié)前、凍結(jié)和凍融過程中均存在土壤含水量的變化，其中，土壤凍結(jié)過程中含水量的增量與土壤凍結(jié)速度之間呈雙曲型相關(guān)關(guān)系。同時，土壤的凍融循環(huán)對地下水埋深也存在一定影響。地下水位在凍結(jié)期降低，在融化期升高，融通后迅速抬高；土壤凍結(jié)過程中地下水初始埋深與地下水最大埋深存在一定的相關(guān)關(guān)系[24]。地下水埋深與凍土層的相互作用還影響了土壤水分入滲的特性：凍融土壤的相對穩(wěn)滲率隨地下水埋深的減小而減??；凍融期間地下水埋深小的土壤入滲能力始終小于地下水埋深大的土壤入滲能力。地下水埋深對凍融土壤水分入滲能力的影響通過其對地表土壤含水量的影響而實(shí)現(xiàn)[25]。

總之，影響土壤凍融循環(huán)規(guī)律的任何因素都是通過改變兩個基本條件來實(shí)現(xiàn)的：土壤溫度和土壤含水量。氣溫變化、植被類型、植物殘余和積雪覆蓋、太陽輻射、氣候變化等都是通過改變土壤溫度而影響土壤的凍融規(guī)律；降水、灌溉、融雪、地下水變化等則通過改變土壤的含水量動態(tài)分布影響凍土的分布格局。不同土壤類型實(shí)質(zhì)上也是由于土壤熱量傳遞性質(zhì)的差異和持水能力的不同具有不同的土壤溫度和含水量，因此產(chǎn)生不同的凍融規(guī)律。

2 土壤凍融對水力學(xué)性質(zhì)的影響

在土壤凍融循環(huán)規(guī)律及其影響因素明確的基礎(chǔ)上進(jìn)行凍土入滲機(jī)理的探索，往往通過建立土壤水入滲方程的方式實(shí)現(xiàn)。入滲方程是水文模型的重要組成部分[26-28]。然而，由于土壤結(jié)構(gòu)復(fù)雜，各向異性，且土壤水滲流過程中水分狀態(tài)和水力學(xué)性質(zhì)是動態(tài)變化的，從而增大了入滲模型化的難度[29]，而凍土的凍融循環(huán)過程則進(jìn)一步加劇了這一難度[30-31]。土壤水入滲在凍土和非凍土內(nèi)部的水力學(xué)性質(zhì)差異巨大，水分在兩種狀態(tài)土壤介質(zhì)中的滲透過程往往會導(dǎo)致數(shù)值模型不能收斂、模型假設(shè)不適用等問題?？傊c非凍土入滲方程相比，面向凍土的入滲方程發(fā)展更加滯后[32]。

一般土壤入滲方程都可由若干基本的土壤水力學(xué)性質(zhì)參數(shù)及其函數(shù)關(guān)系(土壤水分特征曲線)確定(如土壤含水量，水力傳導(dǎo)度，土壤水勢等參數(shù)及其函數(shù)關(guān)系)。這些參數(shù)及其關(guān)系通常與土壤類型、體積密度、粒徑分布等土壤屬性有關(guān)[33-34]。對于凍土凍融循環(huán)過程中的土壤入滲機(jī)理，應(yīng)考慮土壤中液態(tài)水的含量和固態(tài)水對水力傳導(dǎo)度的阻抗[30-31]。

2.1 土壤水力學(xué)性質(zhì)參數(shù)

由于凍土入滲模型研究的相對滯后性，導(dǎo)致刻畫土壤水入滲的大部分參數(shù)取值均來自非凍土的室內(nèi)外實(shí)測，因此，這些參數(shù)不能直接用在凍土入滲的模擬中。水文模型引用較多的幾個參數(shù)取值，如Clapp和Hornberger[33]利用經(jīng)驗(yàn)?zāi)Ｐ陀嬎愕?1種無機(jī)土壤和Letts等[35]給出的3種有機(jī)土壤的土壤孔隙度、土壤水勢、土壤水力傳導(dǎo)度等參數(shù)，土壤樣本均來自非凍土區(qū)，這些參數(shù)值不適用于凍土區(qū)入滲機(jī)理的刻畫。

同種土壤相同含水量條件下，凍結(jié)后的有效孔隙小于未凍結(jié)土壤的有效孔隙，因?yàn)橥寥揽紫吨械牟糠只蛘呷克中纬杀?，阻塞了土壤孔隙，降低了液態(tài)水的儲量和傳導(dǎo)度。凍土中水力傳導(dǎo)度對溫度的變化極為敏感，在凍土凍結(jié)極小的溫度區(qū)間內(nèi)，水力傳導(dǎo)度能出現(xiàn)數(shù)量級上的變化。如Burt和Williams[9]觀測到細(xì)沙凍土溫度從-0.26 ℃上升到0 ℃時，水力傳導(dǎo)度增大了近8個數(shù)量級。由此可見，土壤凍結(jié)狀態(tài)與非凍結(jié)狀態(tài)下，水力傳導(dǎo)度存在質(zhì)的變化。

2.2 土壤水分特征曲線

水文模型中比較常用的3組參數(shù)關(guān)系方程如式(1)～式(3)所示：

(1)

(2)

(3)

式中：ψ為土壤基質(zhì)勢，m；ψ0為飽和土壤基質(zhì)勢，m；θl為土壤體積含水率，m3/m3；θS為飽和土壤體積含水率，m3/m3；θr為土壤殘留含水率，m3/m3；Se為有效飽和含水率，m3/m3；b為經(jīng)驗(yàn)系數(shù)；K為不飽和水力傳導(dǎo)度，m/s；KS為飽和水力傳導(dǎo)度，m/s；α為常量，1/m；λ、m和n為經(jīng)驗(yàn)系數(shù)。其中，式(1)由Clapp和Hornberger[33]提供，式(2)由Brooks和Corey[36]提供，式(3)由Van Genuchten[37]提供。

上述土壤參數(shù)關(guān)系方程用于凍土?xí)r需要加以修正，以體現(xiàn)土壤中冰晶對入滲的影響。通常的做法是：(1)用土壤液體體積含水率代替土壤總體積含水率[3]。(2)飽和水力傳導(dǎo)度考慮一個阻抗因子[38]。(3)用飽和土壤液體體積含水率代替飽和土壤總體積含水率[39]。

3 入滲模型研究進(jìn)展

入滲模型用于刻畫降水、灌溉、融雪等水分向土壤入滲的過程，該過程是土壤水分運(yùn)動的重要形式之一。土壤入滲過程可以由土壤的入滲強(qiáng)度描述，而入滲強(qiáng)度取決于水分供應(yīng)條件(降水、灌溉強(qiáng)度、灌溉方式等)和土壤滲水性能(與土壤初始濕度、基質(zhì)勢、土壤質(zhì)地、結(jié)構(gòu)等有關(guān))。適用于各種條件的入滲模型很多，但大致分為三類：理論模型、半經(jīng)驗(yàn)?zāi)Ｐ秃徒?jīng)驗(yàn)?zāi)Ｐ蚚40]。

3.1 經(jīng)典入滲模型

Green-Ampt模型和Philip模型是兩個應(yīng)用比較廣泛的理論模型。Green-Ampt模型又稱為活塞模型，是建立在毛管理論基礎(chǔ)之上的具有一定物理基礎(chǔ)的簡化模型[41]。該模型假設(shè)土壤為均質(zhì)且厚度足夠大，土壤初始含水率在剖面上分布均勻，地表積水且深度恒定，濕潤鋒面水平分布且飽和含水率、初始含水率在鋒面上下截然分開[42]。

Green-Ampt模型常用的基本形式如式(4)所示：

(4)

式中：finf,t為t時刻的入滲速率,mm/h；KS為飽和水力傳導(dǎo)度,mm/h；ψwf為濕潤峰處的土壤基質(zhì)勢,mm；Δθv為濕潤峰兩端的土壤含水率差值；Fint,t為t時刻累積的入滲量,mm。

為了增強(qiáng)Green-Ampt模型在不同條件下的適用性，許多學(xué)者對模型進(jìn)行了改進(jìn)，使對入滲的刻畫更符合實(shí)際[29,43-44]。土壤具有很強(qiáng)的空間變異性，各個方向可能存在不同的質(zhì)地、結(jié)構(gòu)、初始含水率等，造成土壤水的入滲并非如Green-Ampt模型假設(shè)的理想狀態(tài)。為了體現(xiàn)土壤垂直方向上的非均勻性，面向分層土壤的Green-Ampt模型及其應(yīng)用不斷出現(xiàn)[45-47]。土壤水分入滲過程中，濕潤鋒后土壤層并非完全飽和[48-49]，因此水力傳導(dǎo)度在剖面上的分布也必然存在差異[50-51]，基于這一觀點(diǎn)的模型改進(jìn)也不斷涌現(xiàn)。地表積水深度和降水、灌溉強(qiáng)度也對土壤入滲過程產(chǎn)生影響，不同積水條件下的修正模型相繼出現(xiàn)[43,52-53]。Philip入滲方程如式(5)所示[54-56]：

I=S·t0.5+A·t

(5)

土壤入滲理論模型是土壤水分運(yùn)動模型(如Richards模型)在特定土壤水分條件和水分供應(yīng)模式下的解析解模型，盡管受很多假設(shè)條件的限制，但仍是水文模型中使用最多的模型，因?yàn)檫@類模型具有較強(qiáng)的物理基礎(chǔ)，且模型參數(shù)可以通過實(shí)地監(jiān)測、土壤類型相關(guān)性或者土壤轉(zhuǎn)換函數(shù)獲得，實(shí)用性強(qiáng)[33-34]。

在生產(chǎn)實(shí)踐中應(yīng)用比較廣泛的經(jīng)驗(yàn)?zāi)Ｐ陀蠯ostiakov公式和Horton公式。其中Kostiakov公式如式(6)所示[58]：

i=i1·t-α

(6)

式中：i為入滲速率,cm/min；i1為起始入滲速率,cm/min,與土壤質(zhì)地和初始土壤含水率有關(guān)；t為時間，h;α為經(jīng)驗(yàn)指數(shù)。該模型在計算短歷時土壤入滲或者初期土壤入滲時結(jié)果較為精確，隨著時間的延長，模型結(jié)果趨于0，與實(shí)際情況不符[41]。

Horton公式如式(7)所示[59]：

i=if+(i0-if)e-β t

(7)

式中：if為穩(wěn)定入滲速率,mm/min；i0為初始入滲速率,mm/min；β為常數(shù)，反映土壤的特性。

3.2 寒區(qū)入滲模型

寒區(qū)由于凍土的凍融使得入滲具有不同于非凍土入滲的特點(diǎn)，需要建立適合寒區(qū)的入滲模型或者在原有模型的基礎(chǔ)上進(jìn)行改進(jìn)。Granger等[60]在15 a雪地水文和5 a凍土入滲實(shí)驗(yàn)(加拿大大草原棕壤和暗棕壤區(qū))的基礎(chǔ)上，建立了融雪入滲、雪水當(dāng)量和土壤含水率之間的數(shù)學(xué)關(guān)系，用于估算整個融雪期的入滲量，公式如式(8)所示：

INF=5(1-θp)SWE0.584

(8)

式中：INF為融雪入滲量，mm；SWE為雪水當(dāng)量，即雪融化后的水量，mm；θp為0～30 cm土壤層的平均含水率，cm3/cm3。該公式是在大量的實(shí)地實(shí)驗(yàn)中總結(jié)出的經(jīng)驗(yàn)公式，與前文的入滲公式最大的區(qū)別在于考慮了融雪的影響。

Zhao和Gray[39]應(yīng)用數(shù)值模擬方法建立了預(yù)測凍土入滲的半經(jīng)驗(yàn)參數(shù)模型，用于刻畫中細(xì)土壤中準(zhǔn)穩(wěn)態(tài)流隨時間變化的累積入滲量，公式如式(9)所示：

[(273.15-TI/273.15)]-0.45t0.44

(9)

式中：C1為土壤的比熱容，J/(kg·℃)；S0為入滲過程中表層土壤飽和度，m3/m3；SI為0～30 cm土層預(yù)融化孔隙飽和度，m3/m3；TI為0～30 cm土層預(yù)融化平均土壤溫度，℃；t為入滲時間,h。該模型將入滲與初始土壤飽和度、表層土壤飽和度、土壤水力傳導(dǎo)度、初始土壤溫度和入滲時間聯(lián)系起來，與實(shí)測數(shù)據(jù)擬合良好。

上述類型的寒區(qū)土壤入滲模型均是經(jīng)驗(yàn)、半經(jīng)驗(yàn)?zāi)Ｐ?，在研究?dāng)?shù)氐耐寥廊霛B時具有較好的應(yīng)用效果，但這些模型似乎都未考慮凍土層對入滲機(jī)制的影響，缺乏土壤凍融循環(huán)的物理基礎(chǔ)。

陸面參數(shù)模式SWAP(Soil-Water-Atmosphere-Plants)中涉及的土壤入滲模型考慮了熱量傳遞引起的土壤凍融對入滲的影響，建立了土壤水-冰相互轉(zhuǎn)化過程中的土壤水入滲方程如式(10)所示[61]：

(10)

式中：I為入滲速率；Ms now為融雪量，mm；Ks為非凍土飽和水力傳導(dǎo)度，m/s；θ為凍土層濕潤區(qū)液態(tài)水含水率，m3/m3；θs為土壤孔隙度，m3/m3；θi為土壤中冰所占孔隙率，m3/m3；θl為土壤靜態(tài)水含水率，m3/m3。該模型也是半經(jīng)驗(yàn)?zāi)Ｐ?，但是考慮了土壤凍融的影響，更符合凍土入滲的機(jī)理。

專門用于刻畫凍土入滲的模型相對較少，一般都是通過大量實(shí)驗(yàn)建立的經(jīng)驗(yàn)、半經(jīng)驗(yàn)?zāi)Ｐ?。理論模型方面，往往直接采用非凍土的入滲模型，不對模型的結(jié)構(gòu)進(jìn)行改進(jìn)，只是對模型中的水力學(xué)參數(shù)進(jìn)行凍土特性的修正。如SWAP中的土壤水分運(yùn)動部分，基本方程仍然是基于達(dá)西定律和質(zhì)量守恒定律的Richards方程，只是對水力傳導(dǎo)度進(jìn)行了修正，將土壤水分凍結(jié)對土壤水運(yùn)動的影響考慮在內(nèi)[62]。

4 結(jié) 論

本文分析了寒區(qū)土壤凍融循環(huán)規(guī)律及其影響因素，已有的研究文獻(xiàn)顯示，影響土壤凍融循環(huán)規(guī)律的任何因素都是通過改變土壤溫度和土壤含水量兩個基本條件來實(shí)現(xiàn)的。討論了土壤凍融對土壤水力學(xué)性質(zhì)的影響，包括土壤水力學(xué)性質(zhì)參數(shù)和土壤水分特征曲線，這些參數(shù)用于凍土?xí)r需要加以修正，以體現(xiàn)土壤中冰晶對入滲的影響。介紹了若干經(jīng)典土壤水入滲模型和寒區(qū)經(jīng)驗(yàn)、半經(jīng)驗(yàn)土壤入滲模型，指出寒區(qū)由于凍土的凍融使得入滲具有不同于非凍土入滲的特點(diǎn)，需要建立適合寒區(qū)的入滲模型或者在原有模型的基礎(chǔ)上進(jìn)行改進(jìn)。

生產(chǎn)實(shí)踐與水文模型發(fā)展的需求推動入滲理論與模型不斷發(fā)展。但是寒區(qū)凍土由于其復(fù)雜的土壤水凍融過程，使得入滲模式產(chǎn)生不同于非凍土的特殊性，需要建立新的入滲模型，或者在原有入滲模型的基礎(chǔ)上進(jìn)行模型結(jié)構(gòu)的再設(shè)計、模型參數(shù)的重修正。這些模型上的創(chuàng)新與改進(jìn)必須基于寒區(qū)凍土凍融機(jī)理和理論的不斷完善，包括凍土凍融機(jī)理、凍土水分運(yùn)動規(guī)律、凍融機(jī)制下土壤水力學(xué)性質(zhì)的變化等。在大量總結(jié)歷史文獻(xiàn)和開展室內(nèi)外實(shí)驗(yàn)的基礎(chǔ)上，進(jìn)行凍土入滲理論和模型的深入研究，對于寒區(qū)水文學(xué)理論發(fā)展和生產(chǎn)實(shí)踐均有重要科學(xué)意義。