鄒福建，張 楚，李洋洋，楊 寧

(1. 徐州市水利建筑設(shè)計(jì)研究院，江蘇 徐州 221000；2. 江蘇省水利勘測(cè)設(shè)計(jì)研究院有限公司徐州分公司，江蘇 揚(yáng)州 221000)

黃河故道的土堤問(wèn)題是多方面的。毀滅性的洪水和季節(jié)性過(guò)多的降雨正在加速黃河故道土堤的破壞過(guò)程，每年都會(huì)增加對(duì)農(nóng)業(yè)和基礎(chǔ)設(shè)施造成巨大的破壞可能性。在過(guò)去的幾十年里，黃河故道的土堤每年都面臨侵蝕、決口等問(wèn)題。河堤主要的損壞原因可以分解為人類(lèi)活動(dòng)、溢流、侵蝕、滲漏和地基滑動(dòng)[1- 3]。此外，施工期間的監(jiān)督不足導(dǎo)致土方工程質(zhì)量差，使用了不合適的施工材料、土壤結(jié)塊、土塊粉碎不徹底、土壤壓實(shí)不足、或表土層鋪設(shè)不足、使用劣質(zhì)材料、維護(hù)不足、河流遷移和公眾砍伐[4- 7]。在諸多原因中，設(shè)計(jì)方法和施工程序不當(dāng)是主要原因路堤破壞最重要的原因之一。土堤的穩(wěn)定性受河流或水庫(kù)附近水位升降過(guò)程中發(fā)生的滲流影響[8- 10]。然而，要解釋特定問(wèn)題的失效現(xiàn)象，必須對(duì)材料行為有清楚的了解。為了獲得適當(dāng)?shù)睦斫?，一些研究人員已經(jīng)嘗試描述潛在的大規(guī)模破壞問(wèn)題(如滑坡或含有有機(jī)物)的土壤的巖土性質(zhì)[11- 13]。在本項(xiàng)研究中，黃河故道的堤防破壞和河岸侵蝕問(wèn)題已經(jīng)從穩(wěn)定性和巖土工程特性的因素方面進(jìn)行了研究。簡(jiǎn)言之，本研究旨在調(diào)查黃河故道上段和黃河故道下段堤防材料的物理和機(jī)械性能，確定黃河故道河岸材料的巖土性質(zhì)，并闡明徐州范圍內(nèi)的黃河故道的河岸破壞機(jī)制，評(píng)估現(xiàn)有設(shè)計(jì)，通過(guò)一個(gè)實(shí)例分析路堤穩(wěn)定性的方法。

土壤樣本取自黃河故道上段右岸堤防的破碎部分和黃河故道下段侵蝕河岸。實(shí)驗(yàn)室試驗(yàn)在上海海洋大學(xué)。試驗(yàn)包括粒度分析、顆粒密度、液限和塑限、壓實(shí)特性、固結(jié)、滲透性、直剪試驗(yàn)和無(wú)側(cè)限抗壓強(qiáng)度試驗(yàn)。在試驗(yàn)結(jié)果的基礎(chǔ)上，對(duì)土壤進(jìn)行了分類(lèi)。采用落差法測(cè)定了水的滲透系數(shù)滲透性，對(duì)不同含水率的試樣進(jìn)行了無(wú)側(cè)限抗壓強(qiáng)度試驗(yàn)。

1 材料和方法

基本物理性質(zhì)參數(shù)見(jiàn)表1。試驗(yàn)結(jié)果表明，土壤的塑限和液限之間存在顯著差異。路堤土的塑性極限是指非塑性(NP)，其中河岸土的塑性極限為27%。從粒徑分布分析可知，路堤土的最大粒徑為425μm，有效粒徑為7.6μm，如圖1所示。據(jù)此，采用工程分類(lèi)法將土壤劃分為細(xì)砂(SF)類(lèi)。同樣，河岸土被劃分為低液限(ML)的粉土。

表1 土壤物理成分

圖1 土壤粒徑分布曲線圖

根據(jù)擊實(shí)試驗(yàn)結(jié)果，在最佳含水量(Wopt)為21.2%時(shí)，河堤土(ES)和河岸土(RS)的最大干密度分別為1.59g/cm3和1.72g/cm3，如圖2所示。河堤土的滲透性和強(qiáng)度值隨含水量的增加而變化，如圖3所示。在含水量為24%時(shí)，滲透系數(shù)最小值為1.29×10-5cm/s，超過(guò)最佳含水量，接近液限(w=25.8%)。干側(cè)單位含水率變化時(shí)，滲透系數(shù)的變化幾乎比最小滲透系數(shù)大一個(gè)數(shù)量級(jí)。

圖2 河堤土壤和河岸土壤的壓實(shí)曲線

圖3 河堤土壤透水性和強(qiáng)度隨含水量變化圖

2 結(jié)果與討論

現(xiàn)實(shí)情況是，沙袋僅在洪水緊急情況下放置。但是，沙袋也在其他時(shí)間保護(hù)河岸坡面，防止雨水和徑流沖刷細(xì)顆粒物，同樣十分重要。在其他一些地方，預(yù)制混凝土塊被代替沙袋。但是，預(yù)制混凝土塊相對(duì)于沙袋非常昂貴，所以在某種程度上是使用范圍有限的。沙子是最便宜且容易獲得的材料。因此，沙袋的河岸、河堤應(yīng)用可能是堤壩防護(hù)的一種切實(shí)可行的方法。除沙袋外，薄層水泥復(fù)合材料可作為護(hù)坡措施。但是，要找到一個(gè)有效的、可持續(xù)的邊坡防護(hù)體系，還需要綜合研究其對(duì)邊坡穩(wěn)定性的影響。含水量對(duì)河岸土巖土性質(zhì)的影響如圖4所示。當(dāng)飽和度超過(guò)90%時(shí)，滲透率為3.5×10-7cm/s。通過(guò)直剪試驗(yàn)，得到了抗剪強(qiáng)度參數(shù)黏聚力(c)和內(nèi)摩擦角(?)分別為153kPa和22°。

圖4 河岸土壤透水性和強(qiáng)度隨含水量變化圖

當(dāng)含水量為12.8%時(shí)，最大抗壓強(qiáng)度為280kPa，低于最佳含水量(16.45%)。但隨著含水量的增加，強(qiáng)度迅速下降，在含水量為22%時(shí)達(dá)到150kPa左右。在峰值后，隨著含水量增加10%，強(qiáng)度降低了46%。

圖5 河岸裂縫

根據(jù)現(xiàn)場(chǎng)觀測(cè)以及樣本材料的試驗(yàn)結(jié)果，確定了研究區(qū)域黃河故道下段的塌岸機(jī)理。破壞與岸線后幾厘米至幾十厘米處的張力裂縫形成有關(guān)，如圖5所示。也有報(bào)道稱(chēng)，在洪水退去期間，當(dāng)河流流速似乎很高時(shí)，會(huì)發(fā)生嚴(yán)重的塌岸。這種類(lèi)型的破壞可以描述為板狀旋轉(zhuǎn)破壞，它包括2種機(jī)制，例如具有拉伸裂紋的平面破壞和涉及具有相同幾何形狀的破壞塊的傾倒。從巖土工程的角度來(lái)看，岸坡土的強(qiáng)度隨著含水率的增加而迅速降低。同時(shí)，滲透性也變得非常低(<1×10-6cm/s)，最終增加了材料的重量并引發(fā)了堤壩的大規(guī)模破壞，如圖6所示。

此外，在洪水結(jié)束時(shí)河水水位下降時(shí)，由于河岸材料的低滲透性，河岸地下水位下降的速度低于地表水位。這一現(xiàn)象造成了河岸水壓力與坡面地表水圍壓的不平衡，也導(dǎo)致了岸坡的突然大規(guī)模破壞。在這種情況下，為了采取必要的措施保護(hù)研究區(qū)的河岸，很有必要找出臨界岸坡和岸坡高度的極限值。河岸穩(wěn)定性分析的極限平衡法可以遵循Darby和Thorne(1996)提出的方法。這是一種成熟的河岸穩(wěn)定性分析方法，適用于沿平面破壞面發(fā)生破壞的陡峭、粘性、非層狀河岸。在這種方法中，還包括孔隙水和靜水圍壓，其中破壞面不受約束地穿過(guò)堤腳。

對(duì)徐州地區(qū)的黃河故道堤防進(jìn)行了案例研究。根據(jù)通識(shí)，將路堤中的含水層自由表面視為假定的飽和線，該飽和線由土壤類(lèi)型決定。本研究采用有限元法進(jìn)行滲流分析，合理確定自由面。采用極限平衡法對(duì)邊坡穩(wěn)定性進(jìn)行了分析，并與已有的計(jì)算值進(jìn)行了比較。

滲流分析的水力模型[14- 15]是模擬非飽和流過(guò)程中使用最廣泛的模型。因此，在本研究中，它被用來(lái)表示導(dǎo)水率和壓頭之間的關(guān)系：

(1)

式中，kh—任何壓頭下的滲透系數(shù)，m/s；ks—飽和滲透系數(shù)，m/s；h—土壤水頭；n—曲線形狀參數(shù)。

(2)

(3)

這些水力參數(shù)(ks，α和n)是通過(guò)使用數(shù)據(jù)庫(kù)程序ROSETTA(Schaap et al.，2001)預(yù)測(cè)的，該程序使用土壤質(zhì)地?cái)?shù)據(jù)來(lái)預(yù)測(cè)參數(shù)值。黃河故道堤防土壤類(lèi)型為砂質(zhì)黏土，其水力參數(shù)ks，α和n分別為3.33×10-7m/s、2.7和1.23。利用這些水力參數(shù)繪制了土壤水分保持曲線，如圖6所示。

圖6 保水曲線

在有限元滲流分析中，各種單元的邊界條件和水力性質(zhì)相同。本次分析采用四邊形單元和三角形單元對(duì)區(qū)域進(jìn)行離散。用不同的元素種類(lèi)和數(shù)量進(jìn)行了4次試驗(yàn)，比較了自由表面位置的變化。前2個(gè)分析考慮了3個(gè)節(jié)點(diǎn)三角形單元，而另外2個(gè)則考慮了4個(gè)節(jié)點(diǎn)四邊形單元。每個(gè)分析設(shè)置域的邊界條件如圖7所示。此外，還計(jì)算了不同離散區(qū)域的截面中部總流量，以便于比較。滲透分析總結(jié)見(jiàn)表2。

圖7 有限元單元和邊界條件

表2 有限元計(jì)算滲漏總結(jié)

在計(jì)算孔隙水壓力分布時(shí)，不需要精確地確定滲透系數(shù)。因此，近似滲透率函數(shù)足以用于分析目的。通過(guò)中間部分的流量略有不同，見(jiàn)表2。結(jié)果表明，對(duì)于任何類(lèi)型和任何數(shù)量的單元，水位自由面的位置幾乎相同。圖8為滲流分析結(jié)果。

圖8 自由面和滲漏面位置

根據(jù)安全系數(shù)(FS)值檢查黃河故道堤防的穩(wěn)定性。在有限元滲流分析的基礎(chǔ)上，對(duì)假定的飽和線(1∶5.45)和有限元滲流分析預(yù)測(cè)的自由面進(jìn)行極限平衡邊坡穩(wěn)定性分析。土壤的飽和和非飽和單位重量分別為120 pcf(18.8kN·m-3)和115 pcf(18.1kN·m-3)。分別采用莫爾-庫(kù)侖強(qiáng)度準(zhǔn)則，抗剪強(qiáng)度參數(shù)c=100psf(4.78kPa)，摩擦角φ=18°。用隨機(jī)生成圓形表面的技術(shù)，從100個(gè)試驗(yàn)表面中找出一個(gè)臨界破壞面。

趾部各破壞面角度下限分別設(shè)為-5°和-45°(-ve為逆時(shí)針?lè)较?。表3為極限平衡分析總結(jié)。

表3 邊坡穩(wěn)定分析結(jié)果

3 結(jié)語(yǔ)

結(jié)果表明，在假定飽和線的情況下，安全系數(shù)被高估了22%～24%。此外，假設(shè)的飽和線不滿足二維均勻各向同性介質(zhì)流動(dòng)的Laplace方程。因此，滲流分析是解決堤防滲流問(wèn)題的必要條件，也是確定堤防安全值和安全設(shè)計(jì)的可靠因素。路堤材料的巖土工程性質(zhì)需要通過(guò)使用添加劑或加固材料(如水泥土、天然纖維或土工合成纖維等)來(lái)改善。還需要通過(guò)使用土工袋、水泥復(fù)合材料和鋼筋等來(lái)保護(hù)邊坡。還建議在邊坡穩(wěn)定性分析之前進(jìn)行滲流分析，以獲得設(shè)計(jì)穩(wěn)定路堤時(shí)更可靠的安全系數(shù)。