周韓 羅璟 楊沛璋 楊鵬

摘要：

土壤滲透性是影響降雨入滲-徑流過程的關(guān)鍵指標(biāo)，對(duì)準(zhǔn)確開展水利水電工程規(guī)劃設(shè)計(jì)中的行洪論證至關(guān)重要。國道G318川藏段地質(zhì)條件復(fù)雜，土壤類型多樣，生態(tài)環(huán)境敏感脆弱，土壤滲透性差異顯著。為研究影響土壤滲透性能的主控因素，對(duì)沿線取得的土壤樣品進(jìn)行變水頭滲透試驗(yàn)，在顆粒分析的基礎(chǔ)上，將土壤樣品分為粉壤土、砂壤土、壤砂土，研究干密度、不均勻系數(shù)、曲率系數(shù)、植物根系對(duì)土壤滲透性的影響規(guī)律。研究結(jié)果表明：當(dāng)土壤處于相同物性條件下時(shí)，土壤滲透系數(shù)隨砂粒含量增加而增大，即滲透系數(shù)：壤砂土＞砂壤土＞粉壤土；在不含植物根系情況下，土壤滲透系數(shù)均與干密度呈負(fù)相關(guān)，與不均勻系數(shù)和曲率系數(shù)呈正相關(guān)，其中，砂壤土敏感性最高；而含植物根系土壤滲透系數(shù)普遍高于不含植物根系的同類土壤；對(duì)于含植物根系土壤而言，植物根系對(duì)粉壤土滲透系數(shù)的影響最大，并且隨土壤顆粒粒徑增大，植物根系對(duì)土壤滲透系數(shù)影響逐漸降低。

關(guān) 鍵 詞：

滲透系數(shù)； 變水頭試驗(yàn)； 回歸分析； 植物根系； 土壤質(zhì)地分類

中圖法分類號(hào)： TU411.4

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

DOI：10.16232/j.cnki.1001-4179.2023.08.021

0 引 言

青藏高原是中國最大、世界海拔最高的高原，被譽(yù)為“世界屋脊”“第三極”“亞洲水塔”，地質(zhì)條件復(fù)雜，氣候條件多變，生態(tài)環(huán)境敏感脆弱。此外，區(qū)域內(nèi)地形地勢(shì)陡峻，暴雨頻發(fā)，容易在短時(shí)間內(nèi)形成山洪災(zāi)害，并誘發(fā)滑坡、泥石流等地質(zhì)災(zāi)害，嚴(yán)重威脅著人民群眾的生命財(cái)產(chǎn)安全。土壤滲透性表征了土壤的下滲能力，是影響降雨入滲-徑流過程的關(guān)鍵指標(biāo)，開展不同土壤滲透性研究可為準(zhǔn)確評(píng)價(jià)山洪災(zāi)害提供了基礎(chǔ)數(shù)據(jù)支撐，對(duì)保證水利水電工程的安全合理開發(fā)具有重要意義。

自Darcy由水通過飽和砂的試驗(yàn)得出達(dá)西定律和Terzaghi提出Terzaghi公式［1］后，大量有關(guān)土壤滲透系數(shù)的研究開始展開：Taylor［2］，Lambe［3］，Lerourei等［4］將多孔介質(zhì)滲流的毛管模型引入到土壤滲流的研究中，對(duì)不同的土壤進(jìn)行滲透試驗(yàn)，通過數(shù)據(jù)統(tǒng)計(jì)總結(jié)出滲透系數(shù)與孔隙比呈半對(duì)數(shù)關(guān)系。Hansen［5］在達(dá)西定律的基礎(chǔ)上提出了一些更加完善的砂土滲透理論。近些年來，有學(xué)者指出無黏性土（如碎石土、砂土等）的滲透系數(shù)一般與其不均勻系數(shù)和曲率系數(shù)呈正相關(guān)［6-8］。而針對(duì)一些黏性土（如黃土等），通過相關(guān)試驗(yàn)發(fā)現(xiàn)黃土的滲透系數(shù)與干密度總體呈負(fù)相關(guān)，并且體積含水率越高，受干密度的影響越大，不同干密度下壓實(shí)黃土的滲透系數(shù)的差異主要在低吸力段，即滲透系數(shù)陡升段和飽和段［9-14］。在高寒地區(qū)，相同干密度情況下，土壤內(nèi)部植物根系或植物纖維會(huì)提高土壤的滲透能力，滲透系數(shù)隨著植物根系含量增加而增大［15-18］。

綜上所述，國內(nèi)外學(xué)者針對(duì)土壤滲透性影響因素進(jìn)行了大量試驗(yàn)研究，但缺乏從土壤質(zhì)地分類角度出發(fā)對(duì)土壤滲透性開展系統(tǒng)研究，難以基于衛(wèi)星土壤數(shù)據(jù)快速獲取區(qū)域滲透性，制約水文模型參數(shù)量化。國道G318川藏段跨度長，穿越高原草甸區(qū)、亞高山森林、高寒草甸等生態(tài)敏感地區(qū)，地質(zhì)環(huán)境復(fù)雜，氣候惡劣溫差較大，土壤類型多樣，土壤滲透性差異顯著。本文選取沿線14個(gè)位置，43個(gè)原狀樣進(jìn)行顆粒分析和變水頭滲透試驗(yàn)，研究干密度、不均勻系數(shù)、曲率系數(shù)、植物根系對(duì)土壤滲透性的影響規(guī)律，揭示土壤類型與滲透特性之間的關(guān)系，為土壤滲透性區(qū)域化方法提供新思路，為建立該地區(qū)水文模型提供基礎(chǔ)數(shù)據(jù)支撐。

1 室內(nèi)試驗(yàn)

1.1 試樣性質(zhì)

為探究取土點(diǎn)淺表層土滲透特性開展室內(nèi)變水頭試驗(yàn)，考慮取樣的場(chǎng)地條件以及時(shí)間條件，取樣點(diǎn)均勻分布在取土點(diǎn)上游，共對(duì)14處取土點(diǎn)流域進(jìn)行淺表層土取樣，共計(jì)43個(gè)原狀土試樣，現(xiàn)場(chǎng)取樣點(diǎn)位置如圖1所示（衛(wèi)星土壤數(shù)據(jù)來源于國家科技資源共享服務(wù)平臺(tái)-國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心-土壤分中心）。

根據(jù)衛(wèi)星土壤數(shù)據(jù)結(jié)合現(xiàn)場(chǎng)取土點(diǎn)位置，初步根據(jù)圖1將取得的土樣劃分為如表1所列的類別。

1.2 室內(nèi)滲透試驗(yàn)及篩分試驗(yàn)

依據(jù)在實(shí)驗(yàn)室內(nèi)測(cè)定滲透系數(shù)的相關(guān)規(guī)范及試驗(yàn)原理，將試驗(yàn)分為常水頭和變水頭兩種，前者適用于滲透性大的粗粒砂、土（即K=10-3～10-2cm/s），變水頭試驗(yàn)適用于滲透系數(shù)小的細(xì)粒砂、土（即K=10-6～10-3cm/s），由表1可得試樣基本為細(xì)粒土，因此采用變水頭滲透試驗(yàn)方法。

變水頭試驗(yàn)裝置如圖2所示，本次室內(nèi)試驗(yàn)采用TST-55型滲透儀，滲透用水從直立水頭滴管流入試驗(yàn)儀器，自下而上穿過土樣。試驗(yàn)過程中，將水頭滴管中液面加至設(shè)計(jì)高度，開動(dòng)秒表，記錄起始水頭H1。經(jīng)時(shí)間t2-t1后，再記錄終止水頭H2，根據(jù)GB/T 50123-1999《土工試驗(yàn)方法標(biāo)準(zhǔn)》，變水頭滲透系數(shù)應(yīng)按下式計(jì)算：

KT=2.3aLA（t2-t1）lgH1H2（1）

式中：kT為水溫為T時(shí)試樣滲透系數(shù)；a為變水頭管的斷面面積，m2；L為試樣高度，cm；A為試樣的斷面面積，m2；t1和t2分別為測(cè)讀水頭的起始和終止時(shí)間，s；H1和H2分別為起始和終止水頭高度，cm；2.3為ln和log的變換因數(shù)。

已有研究表明，土壤的滲透系數(shù)與其不均勻系數(shù)及曲率系數(shù)擁有較好的相關(guān)性，所以在滲透試驗(yàn)完成后，對(duì)試樣烘干并碾碎，使用細(xì)粒土標(biāo)準(zhǔn)篩網(wǎng)對(duì)重塑樣進(jìn)行篩分，同時(shí)篩分出土樣中的植物根系并稱重。

2 試驗(yàn)結(jié)果及分析

2.1 試驗(yàn)結(jié)果

依據(jù)美國土壤質(zhì)地分類中砂粒組分含量分類標(biāo)準(zhǔn)，將篩分完成后的土樣進(jìn)行分類，如表2所列。

將室內(nèi)篩分試驗(yàn)的結(jié)果與衛(wèi)星土壤數(shù)據(jù)進(jìn)行對(duì)比，發(fā)現(xiàn)篩分試驗(yàn)所獲得的土壤質(zhì)地分類與衛(wèi)星土壤數(shù)據(jù)質(zhì)地分類有所出入。詳細(xì)情況見表3，同時(shí)為較好描述衛(wèi)星數(shù)據(jù)質(zhì)地分類與篩分試驗(yàn)質(zhì)地分類差異，定義：

β=衛(wèi)星數(shù)據(jù)土壤質(zhì)地分類個(gè)數(shù)室內(nèi)試驗(yàn)土壤質(zhì)地分類個(gè)數(shù)×100%（2）

其中，β為衛(wèi)星數(shù)據(jù)質(zhì)地分類與室內(nèi)試驗(yàn)的質(zhì)地分類重合度。

試驗(yàn)結(jié)果對(duì)于所有土壤質(zhì)地分類β為79.07%，其中粉壤土β最好，達(dá)100%；砂壤土β次之，為78.57%；壤砂土β最差，為45.45%。

由此可以說明：對(duì)于粉壤土而言，衛(wèi)星土壤數(shù)據(jù)與實(shí)際現(xiàn)場(chǎng)數(shù)據(jù)重合度較高，能夠通過衛(wèi)星土壤數(shù)據(jù)對(duì)現(xiàn)場(chǎng)區(qū)域的粉壤土進(jìn)行初步劃分。

此次滲透試驗(yàn)結(jié)果的處理當(dāng)中，根據(jù)同一個(gè)取土點(diǎn)同一類土進(jìn)行滲透系數(shù)的求平均值處理，同時(shí)去除掉取土點(diǎn)8、9這些只含一個(gè)土樣的取土點(diǎn)，以期減少試驗(yàn)當(dāng)中的誤差。同時(shí)一個(gè)試樣中是否含有植物根系對(duì)滲透系數(shù)有較大影響，依據(jù)此次試驗(yàn)結(jié)果分析認(rèn)為土壤中植物根系質(zhì)量占單位體積土壤總質(zhì)量1%以上的試樣其滲透系數(shù)被植物根系影響。

此次滲透試驗(yàn)結(jié)果如表4所列。

2.2 干密度對(duì)土壤滲透性的影響

土壤干密度指從原位取出的土在經(jīng)過孔烘箱烘干至完全脫水后的密度，即所取土壤的固體顆粒質(zhì)量與總體積的比值。本文試驗(yàn)所用土壤干密度多處于1.2～1.6 g/cm3區(qū)間之內(nèi)，土壤的干密度越大，單位體積土壤內(nèi)的固體顆粒越多，孔隙及水流通道越少，土壤更加致密，難以滲透。由圖3可知滲透試驗(yàn)所用粉壤土、砂壤土、壤砂土在干密度上升時(shí)，其對(duì)應(yīng)的滲透系數(shù)呈現(xiàn)下降趨勢(shì)，一般地，含有植物根系的試樣其滲透系數(shù)要比未含有植物根系的試樣大。同時(shí)，圖3（a）中兩條回歸線間距較大，說明植物根系對(duì)粉壤土滲透性影響較大，而圖3（b）中兩條回歸線間距很小，幾乎重合，說明隨著土壤中大粒徑砂礫含量的增加，植物根系對(duì)土壤滲透系數(shù)的影響逐漸降低。

由于現(xiàn)場(chǎng)土壤干密度較其他參數(shù)更易獲取，故使用干密度獲取擬合經(jīng)驗(yàn)公式，如表5所列。

2.3 不均勻系數(shù)對(duì)土壤滲透性的影響

不均勻系數(shù)Cu反映了土壤顆粒組成離散程度，其數(shù)值越大，則表示土顆粒分布愈不均勻，粗顆粒與細(xì)顆粒的大小相差較為懸殊。而對(duì)于Cu較小的土，其顆粒粒徑區(qū)間變化范圍較窄，粒徑區(qū)間單一，大顆粒間的孔隙缺乏較小顆粒的填充，進(jìn)而在土壤內(nèi)部形成較大的孔隙，水流滲透相對(duì)較簡(jiǎn)單。

本次試驗(yàn)土樣的不均勻系數(shù)多位于10～30之間。對(duì)于粒徑級(jí)配曲線連續(xù)的土，不均勻系數(shù)越大，表示土中含有的不同粒徑級(jí)的粗細(xì)顆粒越多，粒組變化范圍越大。由圖4可知當(dāng)粉壤土，砂壤土，壤砂土的不均勻系數(shù)上升時(shí)，其對(duì)應(yīng)的滲透系數(shù)呈現(xiàn)上升趨勢(shì)，在相同不均勻系數(shù)情況下，3類土的滲透系數(shù)大小順序?yàn)椋喝郎巴粒旧叭劳粒痉廴劳痢?/p>

同時(shí)，含植物根系的大部分試樣其滲透系數(shù)大于未含植物根系的試樣。對(duì)于含植物根系土壤而言，植物根系對(duì)粉壤土滲透系數(shù)的影響最大，并且隨土壤顆粒粒徑增大，植物根系對(duì)土壤滲透系數(shù)影響逐漸降低。

2.4 曲率系數(shù)對(duì)土壤滲透性的影響

曲率系數(shù)是反映土顆粒徑分布曲線形態(tài)的系數(shù)，定義為顆粒級(jí)配曲線上，累積顆粒含量為30%的粒徑d30的平方與累積顆粒含量為60%的粒徑d60和有效粒徑d10乘積的比值。該系數(shù)是描述級(jí)配曲線平滑程度的指標(biāo)，一般與不均勻系數(shù)配合一起確定土的級(jí)配的好壞，用來判斷土的工程性質(zhì)。

一般地，土的曲率系數(shù)與滲透系數(shù)呈現(xiàn)正相關(guān)的相關(guān)關(guān)系。對(duì)于此次試驗(yàn)3類土，其不均勻系數(shù)與滲透系數(shù)的關(guān)系如圖5所示。同時(shí)也可得，隨著土壤中大粒徑砂礫含量的增加，植物根系對(duì)土壤滲透系數(shù)的影響逐漸降低。

2.5 不同類型土壤滲透性差異

由圖6可知在相同干密度、不均勻系數(shù)、曲率系數(shù)的情況下，滲透系數(shù)：壤砂土＞砂壤土＞粉壤土。

2.6 滲透性影響因素耦合分析

大量研究表明，土壤的滲透特性是在諸多外在和內(nèi)在因素耦合作用下形成的［20］。其中，內(nèi)在因素主要指土壤的理化性質(zhì)，例如土壤干密度、粒徑級(jí)配、不均勻系數(shù)、曲率系數(shù)、孔隙度、礦物組成成分等，各因素之間互相影響，相互耦合，共同決定土壤的理化性質(zhì)，并進(jìn)一步影響土壤的滲透特性。例如，在其他情況一致的情況下，土壤的不均勻系數(shù)越大，級(jí)配越不良，孔隙度越大，干密度越小，則滲透系數(shù)越大。

外在因素中，生長在土壤中的植被根系是影響土壤滲透特性的主要因素［21-22］。植物根系在土壤中穿插固結(jié)土壤，在根系生長的過程中，土壤中形成新的孔隙，連接形成通道，改變土壤結(jié)構(gòu)。同時(shí)，根系分泌物可以增加土壤中微型濾水團(tuán)聚體［23］的含量，阻止土顆粒的搬運(yùn)，從而降低土壤的干密度。在多因素耦合作用下，干密度降低，土壤孔隙度升高，粒度的分布范圍變大，引起土壤不均勻系數(shù)、曲率系數(shù)升高，進(jìn)而提高土壤的滲透性能。

3 結(jié) 論

本文對(duì)國道G318川藏段沿線14個(gè)取樣點(diǎn)的43個(gè)原狀土樣進(jìn)行顆粒分析和變水頭滲透試驗(yàn)，按照美國土壤質(zhì)地分類標(biāo)準(zhǔn)將其劃分為粉壤土、砂壤土、壤砂土，研究了干密度、不均勻系數(shù)、曲率系數(shù)、植物根系對(duì)土壤滲透性的影響規(guī)律，得到以下結(jié)論：

（1） 國道G318沿線土壤類型分布中粉壤土衛(wèi)星土壤數(shù)據(jù)與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)高度重合，可通過衛(wèi)星土壤數(shù)據(jù)對(duì)粉壤土進(jìn)行初步區(qū)劃。

（2） 在不含植物根系情況下，土壤滲透系數(shù)均與干密度呈負(fù)相關(guān)，與不均勻系數(shù)和曲率系數(shù)呈正相關(guān)，其中，砂壤土敏感性最高；隨著土壤砂礫含量的增加，干密度對(duì)土壤滲透系數(shù)的影響呈現(xiàn)上升的趨勢(shì)。

（3） 對(duì)于含植物根系土壤而言，植物根系對(duì)粉壤土滲透系數(shù)的影響最大，并且隨土壤顆粒粒徑增大，植物根系對(duì)土壤滲透系數(shù)影響逐漸降低。

（4） 當(dāng)處于相同物性條件下時(shí)，即相同干密度或者不均勻系數(shù)或者曲率系數(shù)條件下，土壤滲透系數(shù)隨砂粒含量增加而增大，即滲透系數(shù)：壤砂土＞砂壤土＞粉壤土。

參考文獻(xiàn)：

［1］ 太沙基，帕克.工程實(shí)用土力學(xué)［M］.蔣彭年，譯.北京：水利電力出版社，1960.

［2］ TAYLOR D W.Fundamentals of soil mechanics［M］.New York：John Wiley & Sons，Inc，1948.

［3］ LAMBE T W，WHITMAN S E.Soil mechanics［M］.New York：John Wiley＆Sons，Inc，1969.

［4］ LEROUEIL S.BOUCLIN G.TAVENAS F.Permeability anisotropy of natural clays as a functionof strain［J］.Canadian Geotechnical Journal，1990，27（5）：568-579.

［5］ HANSEN D.Discussion of “on the use of the Kozeny-Carman equation to predict the hydraulic conductivity of soils［J］.Canadian Geotechnical Journal，2004，41（5）：990-993.

［6］ 張國棟，鄧全勝，李泯蒂，等.干濕循環(huán)對(duì)水庫滑坡滲透系數(shù)影響的試驗(yàn)研究［J］.人民長江，2016，47（22）：79-82.

［7］ 朱崇輝，劉俊民，王增紅.無粘性粗粒土的滲透試驗(yàn)研究［J］.人民長江，2005，36（11）：55-57.

［8］ 蘇立君，張宜健，王鐵行.不同粒徑級(jí)砂土滲透特性試驗(yàn)研究［J］.巖土力學(xué)，2014，35（5）：1289-1294.

［9］ 王鐵行，盧靖，張建鋒.考慮干密度影響的人工壓實(shí)非飽和黃土滲透系數(shù)的試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2006（11）：2364-2368.

［10］ 鄒錫云，許強(qiáng)，趙寬耀，等.浸水作用下黑方臺(tái)黃土滲透特性變化微觀研究［J］.人民長江，2020，51（6）：166-171.

［11］ 梁燕，邢鮮麗，李同錄，等.晚更新世黃土滲透性的各向異性及其機(jī)制研究［J］.巖土力學(xué)，2012，33（5）：1313-1318.

［12］ 趙彥旭，張虎元，呂擎峰，等.壓實(shí)黃土非飽和滲透系數(shù)試驗(yàn)研究［J］.巖土力學(xué)，2010，31（6）：1809-1812.

［13］ 趙樂，倪萬魁，張鎮(zhèn)飛，等.不同干密度壓實(shí)黃土垂直積水入滲特性［J］.科學(xué)技術(shù)與工程，2021，21（33）：14310-14317.

［14］ 李華，李同錄，張亞國，等.不同干密度壓實(shí)黃土的非飽和滲透性曲線特征及其與孔隙分布的關(guān)系［J］.水利學(xué)報(bào)，2020，51（8）：979-986.

［15］ 范昊明，李貴圓，周麗麗，等.凍融作用對(duì)草甸土物理力學(xué)性質(zhì)的影響［J］.沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，2011，42（1）：114-117.

［16］ 李勇，徐曉琴，朱顯謨，等.黃土高原植物根系強(qiáng)化土壤滲透力的有效性［J］.科學(xué)通報(bào)，1992（4）：366-369.

［17］ 高金花，李明，陳卓，等.河床基底重構(gòu)復(fù)合土控滲試驗(yàn)研究［J］.人民黃河，2022，44（8）：136-140.

［18］ 黃鋼，鄭明新，董旭，等.草灌植物根系對(duì)煤系土加固效果研究［J］.北京交通大學(xué)學(xué)報(bào)，2020，44（4）：49-56.

［19］ 吳克寧，趙瑞.土壤質(zhì)地分類及其在我國應(yīng)用探討［J］.土壤學(xué)報(bào)，2019，56（1）：227-241.

［20］ 吳宏偉.大氣-植被-土體相互作用：理論與機(jī)理［J］.巖土工程學(xué)報(bào)，2017，39（1）：1-47.

［21］ 鄒全，王桂堯，張文豪，等.百喜草根系對(duì)土壤滲透性能的影響［J］.公路與汽運(yùn)，2022（3）：60-64.

［22］ 桑凱新，胡淦林，黃超，等.黃河河岸帶5種植物類型根系結(jié)構(gòu)特征對(duì)土壤滲透性的影響［J］.中國水土保持科學(xué)，2020，18（5）：1-8.

［23］ 張偉偉，江朝華，程星，等.草本植物根系對(duì)黃河故道區(qū)非飽和土特性的影響［J］.水利水電科技進(jìn)展，2017，37（1）：73-78.

（編輯：黃文晉）

Abstract：

Soil permeability is a key factor affecting the rainfall infiltration-runoff process，which is very important to accurately carrying out the flood discharge demonstration in the planning and design of hydropower and water conservancy projects.The geological conditions along the "G318" Sichuan-Tibet section are complex，the soil types are diverse，the ecological environment is sensitive and fragile，and the soil permeability is significantly different.To study the main influence factors affecting soil permeability，the variable-head permeability test was carried out on soil samples along the route.Based on particle size，the soil samples were divided into silty loam，sandy loam，and loamy sandy soil according to the classification standard of sand component content in the soil texture classification.The effects of dry density，non-uniformity coefficient，curvature coefficient，and plant roots on soil permeability were studied.The results show that when the soil is in the same physical properties，the soil permeability coefficient increases with the increase of sand content，the permeability coefficient is in an order of loamy sand>sandy loam>silty loam.In the condition of no plant roots，the soil permeability coefficient is negatively correlated with dry density，and positively correlated with non-uniformity coefficient and curvature coefficient，among which the sensitivity of sandy loam is the highest.The permeability coefficient of soil with plant roots is generally higher than that of the same soil without plant roots.The influence of plant roots on the permeability coefficient of silty loam soil is the largest，and the influence of plant roots on the permeability coefficient of soil gradually decreases with the increase of soil particle size.

Key words：

permeability coefficient；varying water-head test；regression analysis；plant roots；soil texture classification