程大偉， 陳茜，安鵬，郭鴻 ，鄭睿

(1. 長安大學(xué) a.環(huán)境科學(xué)與工程學(xué)院;b.旱區(qū)地下水文與生態(tài)效應(yīng)教育部重點實驗室;c.地質(zhì)工程與測繪學(xué)院，西安 710000； 2. 陜西理工學(xué)院 土木工程與建筑學(xué)院，陜西 漢中723001)

?

一維穩(wěn)態(tài)流非飽和土滲透系數(shù)垂直分布模型及其線性簡化

程大偉1a，b， 陳茜2，安鵬1c，郭鴻2，鄭睿1a

(1. 長安大學(xué) a.環(huán)境科學(xué)與工程學(xué)院;b.旱區(qū)地下水文與生態(tài)效應(yīng)教育部重點實驗室;c.地質(zhì)工程與測繪學(xué)院，西安 710000； 2. 陜西理工學(xué)院 土木工程與建筑學(xué)院，陜西 漢中723001)

摘要：滲透系數(shù)是控制地下水流動的重要參數(shù)，對滲透系數(shù)的空間分布規(guī)律進行研究具有重要的意義。基于Gardner模型獲取了一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直分布模型，該模型用指數(shù)函數(shù)描述，受飽和滲透系數(shù)和無量綱的深度與流動率等因素的控制；該模型表明一維穩(wěn)態(tài)流條件下均質(zhì)典型土類的滲透系數(shù)其沿垂直方向變化趨勢主要受比流量與飽和滲流系數(shù)的負數(shù)值二者之間的相對大小影響。分別采用泰勒級數(shù)方法和以地下水位處及地表處的滲透系數(shù)作為控制條件方法對一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直分布模型進行線性近似簡化。采用泰勒級數(shù)方法獲取的簡化模型其計算誤差隨無量綱的深度增大而增大。簡化后的模型具有形式簡單、參數(shù)少等特點。通過算例對比簡化模型與原模型的差異，計算結(jié)果表明：采用以地下水位處及地表處的滲透系數(shù)作為控制條件的方法進行線性近似簡化的模型計算誤差比采用泰勒級數(shù)方法獲得的線性近似簡化模型的計算誤差小。

關(guān)鍵詞：一維穩(wěn)態(tài)流;滲透系數(shù);垂直分布;簡化模型

非飽和土的滲透系數(shù)是控制土層中水分和污染物遷移的重要參數(shù)之一[1]。試驗獲取滲透系數(shù)的方法分為穩(wěn)態(tài)方法和非穩(wěn)態(tài)方法的瞬時截面法[2]。滲透試驗的測試結(jié)果表明，滲透系數(shù)對飽和度的依賴性較大[3]，在低飽和度時，滲透系數(shù)的測試非常困難。采用經(jīng)驗公式和理論預(yù)測的方式確定非飽和土的滲透系數(shù)往往是較為理想的選擇。為此，許多學(xué)者對此進行過探討，并提出了一些非飽和土滲透系數(shù)的經(jīng)驗公式或是嘗試經(jīng)驗預(yù)測非飽和土的滲透性函數(shù)[4-7]。

滲透系數(shù)的空間分布與許多巖土工程問題有關(guān)，如入滲或蒸發(fā)條件下所引起的巖土質(zhì)邊坡穩(wěn)定[8-10]以及堤壩穩(wěn)定性分析[11]、污染物入滲運移[12-13]等。本文擬以滲透函數(shù)Gardner模型為基礎(chǔ)獲得一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直方向分布模型，并對該模型進行簡化，使其簡潔易用。

1一維穩(wěn)態(tài)流均質(zhì)非飽和土滲透系數(shù)沿垂直分布模型

很多模型可以用來研究滲透系數(shù)對基質(zhì)吸力的依賴特征。其中Gardner模型[14]廣泛應(yīng)用于求解非飽和流動的解析解，該模型可表達為：

(1)

式中：k為滲透系數(shù)；ks為飽和土滲透系數(shù)，對于均質(zhì)土ks為一單值，與位置無關(guān)；α為土水特征曲線擬合參數(shù)，又稱為孔徑分布參數(shù)，kPa-1；(ua-uw)為基質(zhì)吸力。

在穩(wěn)定狀態(tài)下，垂直非飽和流動可用達西定律進行描述，并依照慣例規(guī)定水流向下運動為負號，向上運動為正號，因此，垂直方向的比流量為

(2)

式中：q為比流量；γw為水的容重；y為垂直方向距地下水位處的距離，取地下水位處y=0，地下水位以上為正。

利用邊界條件，在地下水位處(y=0)基質(zhì)吸力為0，并結(jié)合式(1)，可得基質(zhì)吸力沿y方向分布模型[15]

(3)

式中：e是自然常數(shù)。

將式(3)代入式(1)，即可獲得一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直分布模型，如式(4)所示。

(4)

式中：kwy為沿垂直方向(y方向)分布的滲透系數(shù)。

式(4)表明，一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向(y方向)分布受飽和滲透系數(shù)ks和無量綱的深度αγwy與流動率q/ks等因素的控制。

進一步地，可以對入滲或蒸發(fā)條件下均質(zhì)典型土類(砂土、粉土和粘土)的滲透系數(shù)kwy沿垂直方向(y方向)分布規(guī)律的變化趨勢進行討論。

如前所述，對于均質(zhì)土其飽和滲透系數(shù)ks是一單值，即與位置無關(guān)；q是比流量，在一維穩(wěn)態(tài)流條件下其值與位置無關(guān)；α為土水特征曲線擬合參數(shù)，γw為水的容重，此二者也與位置無關(guān)。

對式(4)求取一階導(dǎo)數(shù)，有

(5)

式(5)的駐點應(yīng)滿足

(6)

對于砂土、粉土和粘土的典型水文參數(shù)如表1[15]所示。

表1　砂土、粉土和粘土的典型水文參數(shù)

如表1所示，典型土類的土水特征曲線擬合參數(shù)α和飽和土滲透系數(shù)ks均大于0。同時,考慮到e-αγwy恒大于零，因此,式(5)的駐點存在的條件是

(7)

對于式(3)而言，比流量q與滲透系數(shù)ks之間存在如式(8)所示關(guān)系。

(8)

當(dāng)q≥-ks時，有。

(9)

式(7)～(9)表明：對于一維穩(wěn)態(tài)流條件下均質(zhì)典型土類的滲透系數(shù)kwy而言，其沿垂直方向(y方向)變化趨勢主要受比流量q與飽和滲流系數(shù)ks的負數(shù)值之間的相對大小影響，而與位置無關(guān)。值得注意的是，按照慣例規(guī)定水流向下運動取負號，水流運動向上運動取正號，同時考慮到典型土類的飽和土滲透系數(shù)ks其值大于0， 因此,對于一維穩(wěn)態(tài)流均質(zhì)非飽和土可以得到如下結(jié)論：

1)若為蒸發(fā)條件，恒有0

2)若為入滲條件，當(dāng)滿足-ks

3)若為入滲條件，當(dāng)滿足q=-ks，滲透系數(shù)kwy沿垂直方向(y方向)保持不變，即為飽和滲透系數(shù)ks。

2一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直分布模型的簡化

一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向分布，如式(4)所示，為非線性函數(shù)，涉及參數(shù)包括飽和土滲透系數(shù)ks、土水特征曲線擬合參數(shù)α和比流量q，以下對式(4)進行近似簡化。

對式(4)利用泰勒級數(shù)展開，可得

(9)

略去二階及以上的高階項，則有：

(10)

式(10)表明，一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向分布可近似用線性函數(shù)表達，且始終過(0，ks)點。但簡化模型式(10)的誤差與無量綱的深度αγwy有關(guān)，隨無量綱的深度αγwy增大，簡化模型的誤差增大。其次，經(jīng)簡化后的表達式形式上雖更為簡潔，但涉及參數(shù)個數(shù)并未減小。

在保證近似后表達式的線性化和過(0，ks)點的特點的前提下，為獲取誤差較小且參數(shù)更少的簡化模型,取地表處y=H，設(shè)該處的滲透系數(shù)為kwH，則滲透系數(shù)kwy沿垂直方向分布必然經(jīng)過(H，kwH)點，故一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向分布還可以近似表達為

(11)

式(11)所涉及的參數(shù)僅有兩個：飽和滲透系數(shù)ks和地表處土的滲透系數(shù)kwH。該簡化模型與原模型相比具有形式簡單、參數(shù)少、參數(shù)易測量等特點。

若把地表處(y=H)的基質(zhì)吸力記為(ua-uw)H,并代入Gardner模型，則有

(12)

將式(12)代入式(11),有

(13)

式(13)是采用地表處基質(zhì)吸力記為(ua-uw)H和飽和滲透系數(shù)ks為參數(shù)的一維穩(wěn)態(tài)流非飽和土滲透系數(shù)垂直分布的線性近似簡化模型。

對于不同的非飽和土其土性、土的密度、濕度和結(jié)構(gòu)的變化均會引起滲透系數(shù)發(fā)生變化。當(dāng)采用Gardner模型描述滲透性函數(shù)時，通過孔徑分布參數(shù)α反映不同土性、土的密度和結(jié)構(gòu)的變化帶來的影響，因而簡化模型式(13)與原模型式(4)能夠反映土性、土的密度和結(jié)構(gòu)的變化的影響。至于在一維穩(wěn)態(tài)流條件下均質(zhì)非飽和土土層內(nèi)其濕度變化對非飽和土滲透系數(shù)影響則是通過垂直方向距地下水位處的距離y的變化間接反映出來。原因在于，一維穩(wěn)態(tài)流條件下均質(zhì)非飽和土層內(nèi)基質(zhì)吸力隨深度單調(diào)變化(一一對應(yīng)關(guān)系)[14],而在吸濕或脫濕條件下基質(zhì)吸力與含水率亦呈單調(diào)變化(一一對應(yīng)關(guān)系)[2]，因而，在一維穩(wěn)態(tài)流條件下均質(zhì)非飽和土層內(nèi)深度與含水率之間具有一一對應(yīng)的關(guān)系。在簡化模型式(13)與原模型式(4)中正是利用這一關(guān)系通過垂直方向距地下水位處的距離y的不同計算并反映土層內(nèi)不同深度處濕度對非飽和土滲透系數(shù)的影響。

3算例

利用前節(jié)所述方法，分別利用滲透系數(shù)沿垂直方向分布的簡化模型和原模型計算粉土和粘土兩種典型土類滲透系數(shù)并進行對比分析。計算土層是厚度為10 m的均質(zhì)土層，土層中的地下水位為y=0 m位置處，地表處為y=10 m。所用水文參數(shù)引自文獻[15],如表2所示。

表2　水文參數(shù)

注：表中q為正值表示蒸發(fā)。

采用簡化模型和原模型分別對粉土、粘土的滲透系數(shù)進行計算，計算結(jié)果如圖1所示。圖中R表示式(4)計算結(jié)果，即原模型計算結(jié)果；R1表示式(10)計算結(jié)果，即對原模型利用泰勒級數(shù)進行線性近似簡化后的計算結(jié)果；R2表示式(11)計算結(jié)果，即以地下水位處和地表處的滲透系數(shù)作為控制條件進行線性近似簡化后的計算結(jié)果。由圖1可見，本算例中的粉土、粘土的滲透系數(shù)，其沿垂直方向分布趨勢具有相似性，即自地下水位處至地表處滲透系數(shù)均呈現(xiàn)出遞減的趨勢，其原因在于本算例中兩種土類的比流量q=-3.14×10-8m/s均比相應(yīng)土類的飽和滲透系數(shù)的負值要大(粉土-ks=-1.0×10-7m/s，粘土-ks=-5.0×10-8m/s)。從圖1中還可以看出，式(10)和式(11)的計算結(jié)果與式(4)計算結(jié)果在變化趨勢上相似，但是式(11)的總體誤差水平更小。

圖1　粉土、粘土的滲透系數(shù)計算結(jié)果對比Fig.1　the calculation comparing between linear model and original model of silt and

圖2　粉土、粘土兩種模型計算值相對誤差Fig.2　the relative error between linear model and original model of silt and clay

為進一步說明上述計算誤差差異情況，圖2給出了簡化模型和原模型計算出的滲透系數(shù)之間的相對誤差隨計算深度的分布規(guī)律。圖2中E1表示式(10)計算結(jié)果與式(4)計算結(jié)果之間的相對誤差隨垂直方向距離y的變化趨勢，E2表示式(11)計算結(jié)果與式(4)計算結(jié)果之間的相對誤差隨垂直方向距離y的變化趨勢。由圖2可見，粉土和粘土的滲透系數(shù)采用式(10)計算時所產(chǎn)生的相對誤差隨垂直方向距離y的增大而增大，最大誤差對于粉土為-42.65%，對于粘土為-4.46%，負號表示計算結(jié)果比式(4)要??；式(11)計算時所產(chǎn)生的相對誤差隨垂直方向距離y的增大呈先增大后遞減的趨勢，最大誤差對于粉土為7.01%，對于粘土為0.96%。這說明一維穩(wěn)態(tài)流均質(zhì)非飽和土滲透系數(shù)kwy沿垂直方向分布可以近似用線性簡化模型描述，但從計算結(jié)果的相對誤差來看，采用式(11)，即以地下水位處和地表處的滲透系數(shù)作為控制條件進行線性近似簡化的模型，計算結(jié)果產(chǎn)生的誤差要小。

4結(jié)論

非飽和含水層滲透系數(shù)的空間分布規(guī)律一直是研究的熱點和難點問題。本文基于Gardner模型獲取了一維穩(wěn)態(tài)流非飽和土滲透系數(shù)沿垂直分布模型，并分別采用泰勒級數(shù)方法和以地下水位處和地表處的滲透系數(shù)作為控制條件方法對模型進行線性簡化。具體結(jié)論如下：

1)一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向(y方向)分布規(guī)律可用指數(shù)函數(shù)描述，受飽和滲透系數(shù)ks和無量綱的深度αγwy與流動率q/ks等因素的控制。

2)對于一維穩(wěn)態(tài)流條件下均質(zhì)典型土類的滲透系數(shù)kwy沿垂直方向(y方向)單調(diào)遞減。

3)一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向分布可用泰勒級數(shù)展開，可近似簡化為線性模型表達，該線性簡化模型的誤差與無量綱的深度αγwy有關(guān)，隨αγwy增大，該線性簡化模型的誤差增大；

4)一維穩(wěn)態(tài)流非飽和土滲透系數(shù)kwy沿垂直方向分布還可用以地下水位處和地表處的滲透系數(shù)作為控制條件進行線性近似簡化，簡化模型與原模型相比具有形式簡單、參數(shù)少、參數(shù)易測量等特點。

5)從算例的計算結(jié)果來看，采用以地下水位處和地表處的滲透系數(shù)作為控制條件進行線性近似簡化的模型計算誤差要比采用泰勒級數(shù)方式獲得的線性近似簡化模型的計算誤差要小。

參考文獻：

[1] 葉為民,錢麗鑫,白云,等.由土-水特征曲線預(yù)測上海非飽和軟土滲透系數(shù)[J].巖土工程學(xué)報,2005,27(11):1262-1265.

YE W M,QIAN L X,BAI Y,et al.Predicting coefficient of permeability from soil-water characteristic curve for shanghai soft soil [J].Chinese Journal of Geotechnical Engineering,2005,27(11):1262-1265. (in Chinese)

[2] FREDLUND D G,RAHARDIO H.非飽和土土力學(xué)[M].陳仲頤，譯.北京:中國建筑工業(yè)出版社,1997.

FREDLUND D G,RAHARDIO H.Soil mechanics for unsaturated soils [M].Trans.Chen Z Y.Beijing:China Building Industry Press,1997. (in Chinese)

[3] 孫大松，劉鵬，夏小和,等.非飽和土的滲透系數(shù)[J].水利學(xué)報,2004,(3):71-75.

SUN D S,LIU P,XIA X H,et al.Permeability coefficient of unsaturated soils [J].Shuili Xuebao,2004,(3):71-75. (in Chinese)

[4] CAMPBELL J D.Pore pressures and volume changes in unsaturated soils [D].Urbana-Champaign University of Illinois,1973.

[5] MUALEM Y.Hydraulic conductivity of unsaturated soils: prediction and formulas [C]//Methods of Soil Analysis.Madison,Wis:American Society of Agronomy,Part 1,1986:799-823.

[6] VANAPALLI S K,FREDLUND D G,PUFAHL D E.The influence of soil structure and stress history on the soil-water characteristics of a compacted till [J].Geotechnique,1999,49(2):143-159.

[7] GARDNER W R.Some steady state solutions of the unsaturated moisture flow equation with application to evaporation from a water table [J].Soil Science,1958,85(4):228-232.

[8] 孔郁斐,宋二祥,楊軍,等.降雨入滲對非飽和土邊坡穩(wěn)定性的影響[J].土木建筑與環(huán)境工程,2013,35(6):16-21.

KONG Y F,SONG E X,YANG J,et al.Rainfall’s effect on the stability of unsaturated slopes [J].Journal of Civil，Architectural & Environmental Engineering,2013,35(6):16-21. (in Chinese)

[9] 于玉貞,林鴻州,李榮建,等.非穩(wěn)定滲流條件下非飽和土邊坡穩(wěn)定分析[J].巖土力學(xué),2008,29(11):2892-2898.

YU Y Z,LIN H Z,LI R J,et al.Stability analysis of unsaturated soil slope under transient seepage flow state [J].Rock and Soil Mechanics,2008,29(11):2892-2898. (in Chinese)

[10] LING H,LING H I.Centrifuge model simulations of rainfall-induced slope instability [J].Journal of Geotechnical and Geo-environmental Engineering,2012,138(9):1151-1157.

[11] 李湛，欒茂田，劉占閣,等.滲流作用下邊坡穩(wěn)定性分析的強度折減彈塑性有限元法[J].水利學(xué)報,2006,37(5):554-559.

LI Z,LUAN M T,LIU Z G,et al.Elasto-plastic FEM for embankment stability analysis under the condition of seepage based on the technique of shear strength reduction [J].Journal of Hydraulic Engineering,2006,37(5):554-559. (in Chinese)

[12] ERTAN D,CORAPCIOGLU M Y,TUNCAY K.Modeling of settlement in saturated and unsaturated municipal landfills [J].International Journal of Geomechanics,ASCE,2006:269-278.

[13] 趙穎,梁冰,薛強,等.地表水入滲對垃圾填埋場水質(zhì)水量影響的數(shù)值模擬分析[J].巖土力學(xué),2010,31(7):2295-2302.

ZHAO Y,LIANG B,XUE Q,et al.Numerical simulation analysis of effect of surface water infiltration on water quality and quantity in landfills [J].Rock and Soil Mechanics,2010,31(7):2295-2302. (in Chinese)

[14] CHEN X H,SONG J X,WANG W K.Spatial variability of specific yield and vertical hydraulic conductivity in a highly permeable alluvial aquifer [J].Journal of Hydrology,2010,388(3/4):379-388.

[15] NING L,WILLIAM J L.非飽和土力學(xué)[M].韋昌富.譯.北京:高等教育出版社,2012.

NING L,WILLIAM J L.Unsaturated soil mechanics [M].Trans.Wei C F.Beijing:Higher Education Press,2012.

(編輯王秀玲)

Permeability coefficient vertical distribution model and it’s linear simplified models of unsaturated soil under unidimensional steady flow condition

Cheng Dawei1a，b,Chen Xi2,An Peng1c,Guo Hong2，Zheng Rui1a

(1a. School of Environmental Science and Engineering; 1b. Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education; 1c. School of Geology Engineering and Geomatics, Chang’an University, Xi′an 710000, P. R. China; 2. School of Civil Engineering and Architecture, Shaanxi University of Techinology, Hanzhong 723001,Shaanxi, P. R. China)

Abstract:The permeability coefficient is important to investigate permeability coefficient spatial distribution of groundwater flow. Based on Gardner model, the permeability coefficient vertical distribution model is established, expressed as exponential, which is controlled by saturation permeability coefficient, dimensionless depth and dimensionless flow rate. The model shows that the permeability coefficient vertical distribution for homogeneous typical soils is influenced by the difference between flow rate and aturation permeability coefficient. And then, simplifies the model is simplified as linear model based on Taylor’s series and permeability coefficient of groundwater level and the earth’s surface as control conditions. The error of linear model with the method of Taylor’s series increases when dimensionless depth increases. The calculation results of the simplified model and the original model are compared . The result shows that the linear model with the method of permeability coefficient of groundwater level and the earth’s surface as control conditions is more accurate than the linear model with the method of Taylor’s series.

Keywords:unidimensional steady flow; permeability coefficient; vertical distribution; simplified model

doi:10.11835/j.issn.1674-4764.2016.02.004

收稿日期：2015-10-23

基金項目：國家自然科學(xué)基金(41230314)；中央高?；究蒲袠I(yè)務(wù)費(310829151076)

作者簡介：程大偉(1984-) ,男,博士,主要從事非飽和土相關(guān)理論研究，(E-mail)250446719@qq.com。

中圖分類號：TU435

文獻標(biāo)志碼：A

文章編號：1674-4764(2016)02-0029-06

Received:2015-10-23

Foundation item：National Natural Science Foundation of China (No. 41230314); Basic Scientific Research Foundation of Central University (No. 310829151076)

Author brief：Cheng Dawei(1984-),PhD,main research interest: the theory of unsaturated soil,(E-mail)250446719@qq.com.