呂媛媛

（杭州交通投資建設(shè)管理集團有限公司，浙江 杭州 310024）

0 引 言

高速公路建設(shè)對我國社會經(jīng)濟發(fā)展起著至關(guān)重要的作用，是推進國家交通網(wǎng)絡(luò)戰(zhàn)略的關(guān)鍵基礎(chǔ)。然而，在大力發(fā)展公路工程建設(shè)的同時也將不可避免地會對周圍環(huán)境產(chǎn)生一定程度的擾動[1]。近年來，由于高速公路的擴張修建而導(dǎo)致邊坡大量挖方填方的現(xiàn)象屢見不鮮，嚴(yán)重改變了自然環(huán)境的生態(tài)面貌，大大削弱了巖土體的抗剪強度，最終導(dǎo)致地質(zhì)災(zāi)害頻頻發(fā)生。因此，為了增強路堤邊坡穩(wěn)定性并且達到防止水土流失、恢復(fù)生態(tài)的目的[2]，高速公路邊坡生物防護技術(shù)得到了廣泛關(guān)注。其中，植物護坡對于生態(tài)文明建設(shè)更是具有重要意義。

目前，已有學(xué)者通過研究發(fā)現(xiàn)植物根系對預(yù)防邊坡淺層滑坡具有重要作用[3]。在植物根系力學(xué)性質(zhì)方面，程洪等[4]采用抗剪強度試驗測定根土復(fù)合體，結(jié)果表明植物根系能有效增強土體的抗剪強度。而稽曉雷等[5]利用數(shù)值模擬方法，發(fā)現(xiàn)根土復(fù)合體的強度隨著狗牙根的根系數(shù)量的增多、根系長度增大而增大。在模型實驗方面，吳宏偉[6]利用根土相互作用模型箱研究了大氣-植被-土體相互作用系統(tǒng)，結(jié)果表明植物的蒸騰作用能顯著提高邊坡的安全系數(shù)；徐文秀等[7]通過植物護坡模型試驗發(fā)現(xiàn)草本植物的根、莖、葉的綜合作用對增強邊坡抗沖刷能力具有顯著效果。

然而，上述學(xué)者的研究僅僅針對某一種植物根系對土體的加固作用，尚未分析不同品種植物根系在不同種植密度下對土體的加筋效果。本文針對千黃高速庫灣高填路基段開展研究，該庫彎路基段由于為一臨湖邊坡，冬夏季水位變化高差大，在進行植物生態(tài)護岸時需優(yōu)先考慮植物在耐淹狀態(tài)下的生長狀態(tài)，為此本文首先通過多種常規(guī)護坡植物的水淹脅迫試驗（百喜草、香根草、狗牙根、結(jié)縷草、高羊茅、黑麥草和果嶺草），篩選出3種不同根系形態(tài)的植物作為研究對象（百喜草、香根草和狗牙根），并通過室內(nèi)試驗獲得植物根系力學(xué)性質(zhì)數(shù)據(jù)為有限元模擬的參數(shù)提供依據(jù)，進一步建立二維平面應(yīng)變模型，分析了所篩選耐淹植物在降雨條件下不同植物種植方案對邊坡的加固作用。

1 植物預(yù)選試驗

為了挑選出3種能夠快速生長、品質(zhì)良好且根系形態(tài)不同的草本植物，首先本文進行了植物預(yù)選試驗，篩選出適合用于該地區(qū)固堤護坡的植物種類。根據(jù)千島湖庫區(qū)野外調(diào)查情況，選取該區(qū)域常見的7種草本植物作為供試植物：結(jié)縷草、果嶺草、黑麥草、高羊茅、狗牙根、百喜草和香根草，進行耐淹性能的試驗研究。

1.1 水淹脅迫對種子萌發(fā)率的影響

本試驗主要對每種植物進行兩種處理：無脅迫處理（對照組）和水淹脅迫處理（試驗組）。其中，試驗組采用室內(nèi)水淹處理種子，水淹時將種子置于水桶內(nèi)，水淹高度為10 cm，每24 h換一次水，水淹時長240 h，水淹結(jié)束后馬上進行萌發(fā)實驗。

萌發(fā)實驗時，首先使用10%的次氯酸鈉溶液對種子進行消毒，然后將對照組和試驗組的種子浸泡10 h。浸泡結(jié)束后，在培養(yǎng)皿底墊一層濾紙，每個濾紙上放置 20粒植物種子，每日在濾紙上均勻噴灑蒸餾水使其保持濕潤并置于植物培養(yǎng)箱中，每日記錄每個培養(yǎng)皿中種子萌發(fā)的數(shù)目。

1.2 水淹脅迫對植物生理特性的影響

本試驗選用10 cm×10 cm×10 cm的盆缽，在每個盆缽中添加等量泥炭土，每盆缽中移入 1 cm高的相同植物幼苗20棵，并于溫室內(nèi)定期澆水養(yǎng)護。待植物生長一個月后，對于無脅迫對照處理組的植物維持原來生長狀態(tài)；對于水淹脅迫處理組，采用雙套盆法，取一個大水箱，水箱中注入水，將水淹脅迫處理組的植物連同盆缽一起放入其中（如圖 1所示），調(diào)整水位至將植物完全淹沒，連續(xù)淹沒168 h后取出植物盆缽。水淹植物取出之后，在正常環(huán)境下放置3 d，3 d后檢測7種植物的生理指標(biāo)。

圖1 水淹處理Fig.1 Flooding treatment

1.3 植物耐淹試驗效果分析

本研究進行了不同植物的水淹脅迫處理，通過檢測4種生理指標(biāo)（葉綠素含量、SOD活性、丙二醛含量和脯氨酸含量）和萌發(fā)率，以耐淹性為主要篩選、評價指標(biāo)，性狀、生長特性為次要指標(biāo)進行比較。

本研究發(fā)現(xiàn)水淹提高了植物的 SOD活性，按照變化大小依次為：百喜草、香根草、狗牙根、高羊茅、結(jié)縷草、黑麥草、果嶺草。其中變化率較大的狗牙根、百喜草和香根草的變化率在80%～90%之間，變化率較小的果嶺草和黑麥草的變化率在40%左右，如圖2所示。

圖2 不同水淹脅迫對植物葉綠素含量的影響Fig.2 Effects of different waterlogging stress on chlorophyll content in plants

植物中脯氨酸含量除了作為植物細胞質(zhì)內(nèi)滲透調(diào)節(jié)物質(zhì)外，還在穩(wěn)定生物大分子結(jié)構(gòu)、降低細胞酸性、解除氨毒以及作為能量庫調(diào)節(jié)細胞氧化還原勢等方面起重要作用。本研究發(fā)現(xiàn)水淹提高了植物的脯氨酸含量，按照變化大小依次為：香根草、狗牙根、百喜草、高羊茅、果嶺草、黑麥草、結(jié)縷草。其中變化率較大的狗牙根和香根草的變化率在100%～120%之間，變化率較小的黑麥草和結(jié)縷草的變化率在70%左右，如圖3所示。

圖3 不同水淹脅迫對脯氨酸含量的影響Fig.3 Effects of different waterlogging stress on proline content

通過整理 SOD活性、丙二醛含量以及萌發(fā)率試驗結(jié)果，然后使用隸屬函數(shù)法對7種植物進行了耐淹性的四因素綜合評價，計算出了7種植物耐淹性綜合評價指數(shù)以及耐淹性能的綜合排序，見表1。

表1 7種植物在水淹處理下因子得分、綜合評價值（D）耐淹性排序Table 1 Factor scores and comprehensive evaluation value(D) ranking of submergence tolerance of 7 plants under waterlogging treatment

通過表1綜合分析顯示狗牙根、香根草和百喜草的耐淹性最強，其次為結(jié)縷草、高羊茅、黑麥草、果嶺草。因此，綜合考慮各方面因素，初步確定狗牙根、百喜草和香根草作為本研究對象。

2 有限元模型建立

植物根系對邊坡的穩(wěn)定性具有重要的加固作用，是增強路堤邊坡穩(wěn)定性并且達到防止水土流失、恢復(fù)生態(tài)環(huán)境的關(guān)鍵方法。本節(jié)將參考千黃高速項目路基邊坡，以狗牙根、百喜草和香根草作為護坡植物，結(jié)合室內(nèi)試驗獲得的植物根系和土體物理力學(xué)數(shù)據(jù)，建立二維模型研究降雨條件下不同植物種植方案對邊坡的加固作用。

2.1 模型幾何尺寸和土體材料參數(shù)

本研究采用 ABAQUS有限元軟件建立二維平面應(yīng)變數(shù)值模型，坡體的幾何參數(shù)參照了千黃高速項目淳安段的路基邊坡橫斷面圖，土體采用常規(guī)理想彈塑性本構(gòu)模型，即Mohr-Coulomb的破壞準(zhǔn)則，強度指標(biāo)采用固結(jié)不排水剪切強度，土坡示意圖如圖 4所示。邊坡模型的左右邊界約束其水平位移，模型的下邊界約束其水平和豎向位移。土體的材料參數(shù)根據(jù)室內(nèi)土工試驗獲得，如表2所示。

表2 土體材料參數(shù)Table 2 Soil material parameters

圖4 土坡示意圖Fig.4 Schematic diagram of soil slope

根據(jù)復(fù)合材料的觀點和方法，在不考慮界面單元的情況下，以土體與植物根系的物理力學(xué)參數(shù)為依據(jù)，將植物根系作為外荷載考慮，直接作用在土體單元上，通過公式將根土復(fù)合體轉(zhuǎn)化為等價的均質(zhì)材料，利用均質(zhì)材料的有限元方法進行求解，分析香根草根系、百喜草根系和狗牙根百喜草混種等3種情況下植物根系對土坡的加固作用。其中，有限元的網(wǎng)格圖劃分如圖5所示，對坡體采用四結(jié)點平面應(yīng)變四邊形單元、雙線性位移。

圖5 網(wǎng)格劃分圖Fig.5 Grid generation diagram

為了有效地反映根系的淺層固坡效果，本研究通過改變坡面邊界條件的孔壓和滲流速度來分別模擬地下水分布和降雨條件。本研究選擇在坡面加入降雨入滲條件，通過對比素土坡和不同分布密度的根土復(fù)合體土坡的應(yīng)力及位移等變量，分析植物根系對土坡穩(wěn)定性的加固效果。

2.2 不同植物根系二維模型及材料參數(shù)

（1）香根草

香根草的根系形態(tài)復(fù)雜，通過數(shù)理統(tǒng)計分析的方法可以大致了解草根的分布規(guī)律。香根草根系能穿透堅硬紅黏土，并穿透礫石和巖層之間薄弱的地方，可長至2～3 m深，最深達到5～6 m，并且有耐淹的特性，是典型的護坡植物。在本研究計算中，由于相對于根系的長度，香根草根系的直徑變化量較小，因此植物根系的直徑可以認為是不變的。根據(jù)香根草根須分布情況，并采用二維的桿單元將其簡化成 4枝須根和 1枝主根，其中主根長度為760 mm，另外兩種須根長為200 mm和370 mm，且與主根夾角分別為 33°和 25°。通過室內(nèi)抗拉試驗可知，在受拉時植物根系呈現(xiàn)出先受拉變形后驟然破壞的現(xiàn)象，根系處于塑性變形的應(yīng)變范圍較窄[8]，因此在材料參數(shù)定義中近似將根系材料設(shè)置成只能夠受拉力的彈性材料。并在網(wǎng)格劃分中采用2節(jié)點桿單元T2D2，方便計算，在簡化后的香根草二維模型如圖6所示。在ABAQUS有限元軟件材料模塊中輸入材料橫截面面積，以實現(xiàn)模擬植物根系的平均直徑，該平均值通過植物根系直徑的統(tǒng)計值確定，由于本文采用二維平面應(yīng)變模型，未計入根系三維尺寸的影響，對其輸入值參考表3。

圖6 香根草根系及其二維模型Fig.6 Vetiver grass root system and its two-dimensional model

（2）百喜草

通過對百喜草進行調(diào)研，可知百喜草的根系主要是須根根系，主根不明顯，具有粗壯、木質(zhì)、多節(jié)的根狀莖。因此，可將百喜草根系簡化成5枝須根，長度為400 mm，根系之間的夾角為25°。如圖7所示。

圖7 百喜草根系及其二維模型Fig.7 Paspalum natatum root system and its two-dimensional model

（3）狗牙根

根據(jù)狗牙根根系須根數(shù)量的分布情況（如圖 8所示），將狗牙根根系簡化成3枝須根，根須的長度為200 mm，根系之間的夾角為45°。

圖8 狗牙根根系及二維模型Fig.8 Bermuda grass root system and its two-dimensional model

（4）根系力學(xué)參數(shù)

植物根系的物理力學(xué)參數(shù)通過查閱相關(guān)文獻資料[9-11]，取值如表3所示。

表3 植物根系力學(xué)參數(shù)Table 3 Mechanical parameters of plant roots

（5）根土接觸定義

根土復(fù)合體采用嵌入式接觸，以根系為嵌入實體，土體為被嵌入實體，由于植物根系的自重相對較輕，在計算過程中可以忽略根系的自重。

本文在研究香根草根系分布密度時，采用1、3、5、8、10株/m等5種情況進行計算分析；狗牙根系的分布密度采用1、4、8、12、15株/m等5種情況進行分析；當(dāng)狗牙根與百喜草混種時，根系的混合比例取 1∶1，其分布密度采用 1、5、10、15株/m等4種分布密度。

3 計算結(jié)果分析

通過有限元數(shù)值模擬計算，發(fā)現(xiàn)在邊坡表面植入植物根系之后，由于土體和植物根系之間的彈性模量差異較大，因此會對邊坡表面土層形成一定程度上的加筋作用。當(dāng)植物根系的分布密度逐步增大時，坡體的水平位移、豎向位移以及總位移都呈現(xiàn)逐漸減小的趨勢，這也說明了植物根系的植入對防止邊坡淺層滑坡具有重要的作用。表 4～6分別對應(yīng)為香根草、百喜草以及狗牙根與百喜草混種工況下的水平位移、豎向位移和總位移的計算結(jié)果。

表4 香根草分布密度對邊坡表層位移的影響Table 4 Effect of density distribution of Vetiver grass on surface displacement of slope

表5 百喜草分布密度對邊坡表層位移的影響Table 5 Effect of density distribution of Paspalum natatum on surface displacement of slope

表6 狗牙根與百喜草1∶1混種時分布密度對邊坡表層位移的影響Table 6 Effect of density distribution on surface displacement of slope with 1∶1 mixture of Bermuda grass and Paspalum natatum

根據(jù)上述計算結(jié)果，將不同種植方案下的水平位移、豎向位移和總位移轉(zhuǎn)換為位移減小量，可得到不同種植密度以及不同植物種植方案與邊坡位移減小量的關(guān)系（如圖9所示）。由圖可知，由于植物根系的加筋作用，使得邊坡表層土體的水平位移、豎向位移以及總位移都明顯減小。其中，植物根系的植入對坡體的水平位移影響最大，其次是對總位移的影響，影響最小的是豎向位移，尤其是植物根系的分布密度對邊坡位移的影響。對比不同植物種植方案以及植物分布密度來看，可以發(fā)現(xiàn)不同分布密度對香根草和百喜草表現(xiàn)出近似線性的提升效果。而當(dāng)分布密度由5株/m轉(zhuǎn)至15株/m過程時，狗牙根與百喜草混種方案所對應(yīng)的位移減小量突然增加而后保持穩(wěn)定線性增長，該現(xiàn)象可以說明當(dāng)該方案分布密度為5～8株/m時，對邊坡加固的效果尤為顯著，混種方案進一步增加種植密度效果有限。

圖9 香根草、百喜草、狗牙根與百喜草混種分布密度與位移減小量關(guān)系圖Fig.9 Relationship between density distribution and displacement reduction of Vetiver grass, Paspalum natatum and Bermuda grass mixed with Paspalum natatum

對比不同植物種類對邊坡土體的加筋效果，可得到如下結(jié)果（見圖10）。

圖10 分布密度和水平位移減小量關(guān)系圖Fig.10 Relation between distribution density and horizontal displacement reduction

由圖10可以明顯看出，在同樣的分布密度下，相較于狗牙根與百喜草混種的種植方案，香根草的根系對坡面土體的加筋效果最好，而百喜草根系對抵抗邊坡位移的效果稍顯薄弱。加之，當(dāng)隨著分布密度的逐漸增加，香根草所對應(yīng)的位移減小量與狗牙根相比差距也逐漸增大。由此可以得出，在本研究中香根草根系的加筋作用最好，究其原因，系香根草根系以主根為主、須根次之，相較于其他兩種植物方案更加深入邊坡土體，對邊坡的加固效果愈好。并由此得出不同方案加固效果排序：香根草方案＞狗牙根與百喜草混種方案＞百喜草方案。

4 結(jié) 論

根據(jù)以上綜合分析，可得出以下結(jié)論：

（1）通過試驗研究狗牙根、香根草和百喜草在水淹處理下種子萌發(fā)率以及植物生理特性，發(fā)現(xiàn)香根草的耐淹性最強，可作為庫灣區(qū)路堤邊坡綠化工程首選植物。

（2）植物根系的加入對邊坡水平位移的影響最大，其次是總位移，影響最小的是豎向位移。

（3）香根草對邊坡的加固作用最好，其次是狗牙根與百喜草混種，而百喜草的加固作用不如前面兩者。

（4）從植物根系形態(tài)方面可知，香根草發(fā)達的主根相較于茂密的須根更能深入邊坡土體，起到更好的護坡作用。