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

程大偉 陳茜 安鵬 郭鴻 鄭睿

摘要：滲透系數(shù)是控制地下水流動的重要參數(shù)，對滲透系數(shù)的空間分布規(guī)律進行研究具有重要的意義?；贕ardner模型獲取了一維穩(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ù)；垂直分布；簡化模型

中圖分類號：TU435

文獻標志碼：A 文章編號：1674-4764（2016）02-0029-06

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 Taylors series and permeability coefficient of groundwater level and the earths surface as control conditions. The error of linear model with the method of Taylors 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 earths surface as control conditions is more accurate than the linear model with the method of Taylors series.

Keywords：unidimensional steady flow； permeability coefficient； vertical distribution； simplified model

非飽和土的滲透系數(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ù)

對于不同的非飽和土其土性、土的密度、濕度和結(jié)構(gòu)的變化均會引起滲透系數(shù)發(fā)生變化。當采用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所示。

采用簡化模型和原模型分別對粉土、粘土的滲透系數(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）的總體誤差水平更小。

為進一步說明上述計算誤差差異情況，圖2給出了簡化模型和原模型計算出的滲透系數(shù)之間的相對誤差隨計算深度的分布規(guī)律。圖2中E1表示式（10）計算結(jié)果與式（4）計算結(jié)果之間的相對誤差隨垂直方向距離y的變化趨勢，E2表示式（11）計算結(jié)果與式（4）計算結(jié)果之間的相對誤差隨垂直方向距離y的變化趨勢。由圖2可見，粉土和粘土的滲透系數(shù)采用式（10）計算時所產(chǎn)生的相對誤差隨垂直方向距離y的增大而增大，最大誤差對于粉土為-4265%，對于粘土為-4.46%，負號表示計算結(jié)果比式（4）要??；式（11）計算時所產(chǎn)生的相對誤差隨垂直方向距離y的增大呈先增大后遞減的趨勢，最大誤差性簡化模型描述，但從計算結(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，F(xiàn)REDLUND 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.Rainfalls 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.

（編輯 王秀玲）