郭學(xué)潮
(山東省鄄城縣水務(wù)局,山東 鄄城 274600)
岸坡穩(wěn)定可以有效保障汛期城鎮(zhèn)安全,防止洪水泛濫。楊文東等[1]使用云模型選取相關(guān)因素劃分權(quán)重對(duì)邊坡穩(wěn)定性進(jìn)行預(yù)測(cè)評(píng)價(jià)研究,該方法可以較為快速、準(zhǔn)確地獲取邊坡穩(wěn)定性情況,適用性較好。邢志龍[2]以青走道水庫(kù)庫(kù)岸邊坡為例,通過(guò)室內(nèi)試驗(yàn)和理論計(jì)算方法,預(yù)測(cè)蓄水后岸坡塌方量。李鵬岳等[3]使用可拓學(xué)模型結(jié)合層次分析法獲取指標(biāo)權(quán)重,綜合評(píng)價(jià)長(zhǎng)壽至豐都段庫(kù)岸穩(wěn)定性。徐鎮(zhèn)凱等[4]認(rèn)為云模型在評(píng)價(jià)邊坡穩(wěn)定性中的應(yīng)用快捷、方便、準(zhǔn)確。許波等[5]以孟底溝水庫(kù)為例,確定影響因素,使用層次分析法進(jìn)行單體庫(kù)岸邊坡穩(wěn)定性評(píng)價(jià)。劉明發(fā)[6]利用Geo-studio研究不同地下水位下邊坡穩(wěn)定性情況,結(jié)果表明地下水是影響邊坡穩(wěn)定性的一個(gè)重要因素。朱文煒等[7]、熊爽等[8]分別使用不同的數(shù)值模擬方法,研究邊坡穩(wěn)定性情況,認(rèn)為巖土體強(qiáng)度是影響邊坡穩(wěn)定性的重要因素。
夏季雨水多發(fā),河道水位變化劇烈,岸坡侵蝕、沖刷現(xiàn)象嚴(yán)重,岸坡穩(wěn)定性必將急劇變化。因此,快速獲取現(xiàn)狀岸坡穩(wěn)定性等級(jí)的現(xiàn)實(shí)意義重大。
岸坡穩(wěn)定性影響因素評(píng)價(jià)指標(biāo)選取應(yīng)當(dāng)全面、客觀。綜合分析,岸坡自身特征、地下水位特征是影響其穩(wěn)定性的主要因素。因此,選取巖土體強(qiáng)度參數(shù)(內(nèi)聚力、內(nèi)摩擦角)、水位(岸坡內(nèi)、外水位)、坡型(坡高、坡比)以及外部荷載因素共7個(gè)定量因素建立評(píng)價(jià)指標(biāo)體系。參考前人研究,可將岸坡穩(wěn)定性分為3個(gè)等級(jí),Ⅰ-Ⅲ級(jí)表示穩(wěn)定性逐漸提高。指標(biāo)體系見(jiàn)圖1。
圖1 評(píng)價(jià)指標(biāo)體系圖
表1 指標(biāo)評(píng)價(jià)等級(jí)表
李德毅等[9]在1995年提出云模型的概念,現(xiàn)在應(yīng)用于風(fēng)險(xiǎn)評(píng)價(jià)、算法改進(jìn)、數(shù)據(jù)挖掘、智能控制等各個(gè)行業(yè),效果顯著。云模型針對(duì)岸坡穩(wěn)定性評(píng)價(jià)中的不確定問(wèn)題,有一定解決問(wèn)題的優(yōu)勢(shì),可獲取不確定性的度量值。
權(quán)重表示某一指標(biāo)對(duì)評(píng)價(jià)結(jié)果影響的重要程度。使用改進(jìn)層次法和熵權(quán)法組合賦權(quán)獲取各指標(biāo)權(quán)重。
2.2.1 改進(jìn)層次分析法
1979年,Satty[10]首次提出層次分析法用于解決評(píng)價(jià)指標(biāo)權(quán)重計(jì)算的問(wèn)題。但人為主觀因素對(duì)權(quán)重值影響較大,因此采用粗糙集[11]改進(jìn)層次分析法,獲取客觀權(quán)重,提高評(píng)價(jià)結(jié)果的準(zhǔn)確性。
計(jì)算過(guò)程如下:
σCB(A)=γC(B)-γC-A(B)
(1)
(2)
式中:σ為不同的集合;γ為集合中的元素;Δσij為各指標(biāo)的相對(duì)重要程度。
2.2.2 熵值法
河道岸坡的穩(wěn)定性具有時(shí)空差異性,是一個(gè)隨時(shí)間、空間變化的結(jié)果。熵權(quán)法使用理論計(jì)算的方法獲取權(quán)重,避免了人為因素?cái)_動(dòng),可獲取客觀權(quán)重。
(3)
(4)
式中:yji為各指標(biāo)實(shí)際值;p為選取指標(biāo)總數(shù)。
2.2.3 組合權(quán)重計(jì)算
(5)
組合權(quán)重為wz:
(6)
式中:α,β為分配系數(shù),滿足以下條件:
(7)
α+β=1
(8)
各指標(biāo)權(quán)重見(jiàn)表2。
表2 評(píng)價(jià)指標(biāo)權(quán)重表
2.2.4 不同指標(biāo)間的相互作用分析
當(dāng)河道水位上漲時(shí),岸坡巖土體類型決定了河道(坡外)水位補(bǔ)給岸坡地下(坡內(nèi))水位的速率,同時(shí)坡外水位也決定了坡內(nèi)水位最終能達(dá)到的水位高度。反之,若是坡內(nèi)水位補(bǔ)給坡外水位時(shí),坡內(nèi)水位的高度決定了坡外水位的高度。兩者之間的轉(zhuǎn)化是一個(gè)動(dòng)態(tài)過(guò)程,此時(shí)岸坡的穩(wěn)定性也是一個(gè)動(dòng)態(tài)變化過(guò)程。
同時(shí)水位的變化,造成坡體內(nèi)部巖土體含水率產(chǎn)生變化,含水率的變化造成巖土體強(qiáng)度(黏聚力和內(nèi)摩擦角)變化。
在選取評(píng)價(jià)指標(biāo)的條件上,李健等[12]提出了云參數(shù)(Ex,En,He)的求解公式:
Ex=(Bmin+Bmax)/2
En=(Bmax-Bmin)/6
He=k
(9)
式中:Bmax,Bmin分別為某等級(jí)標(biāo)準(zhǔn)的最大和最小邊界;k為常數(shù)。
為了評(píng)價(jià)結(jié)果清晰可辨別,作出以下修改,見(jiàn)表3。
表3 云模型的數(shù)字特征表
使用云模型確定參數(shù)Ex,En,He,生成云模型見(jiàn)圖2。
綜合確定度U計(jì)算方法如下:
(10)
式中:δi為確定度;ωz為指標(biāo)組合權(quán)重,岸坡穩(wěn)定性等級(jí)L為:
L=max{U1,U2,…,Un}
(11)
圖2 各評(píng)價(jià)指標(biāo)云圖
以山東菏澤市某河流岸坡為例,使用本文研究方法評(píng)價(jià)其穩(wěn)定性,并與專家現(xiàn)場(chǎng)評(píng)價(jià)結(jié)果進(jìn)行對(duì)比,驗(yàn)證本文所提出方法的準(zhǔn)確性。隨著雨季到來(lái),河道水位上漲,河水侵蝕、沖刷造成河道岸坡(未治理段)形態(tài)改變,同時(shí)坡體內(nèi)部地下水位發(fā)生改變,造成岸坡穩(wěn)定性急劇變化。通過(guò)現(xiàn)場(chǎng)調(diào)查以及室內(nèi)力學(xué)實(shí)驗(yàn),岸坡主要巖土體類型為堆積層,力學(xué)試驗(yàn)結(jié)果表明內(nèi)聚力為22 kPa,內(nèi)摩擦角為18°,由于天然狀態(tài)下含水率較高,飽和狀態(tài)下強(qiáng)度參數(shù)變化較小。其中一段河道岸坡P1坡比1.75,坡高3.0 m,工程區(qū)人類工程活動(dòng)較少,坡頂荷載選擇為0 kPa,岸坡內(nèi)、外水位分別為1.5和3 m(平均高度)。
以內(nèi)聚力為22 kPa為例,使用云發(fā)生器確定該指標(biāo)隸屬于各級(jí)別的確定度:U(Ⅰ)=0.000 0,U(Ⅱ)=0.009 6,U(Ⅲ)=0.345 5,所以指標(biāo)內(nèi)聚力為22 kPa屬于Ⅲ級(jí)。綜合計(jì)算5處岸坡P1-P5的穩(wěn)定性評(píng)價(jià)結(jié)果,評(píng)價(jià)結(jié)果見(jiàn)表4。
表4 評(píng)價(jià)結(jié)果對(duì)比表
由表4可知,P1-P5的穩(wěn)定性等級(jí)依次為Ⅱ、Ⅲ、Ⅰ、Ⅰ、Ⅲ;組合賦權(quán)-云模型獲取岸坡穩(wěn)定性等級(jí)與現(xiàn)場(chǎng)專家評(píng)價(jià)結(jié)果一致,具有一定的適應(yīng)性。
1) 通過(guò)組合賦權(quán)獲取選取7個(gè)指標(biāo)的權(quán)重,巖土體強(qiáng)度參數(shù)為影響河道岸坡穩(wěn)定性的主要因素。
2) 通過(guò)云模型確定P1-P5的穩(wěn)定性依次為Ⅱ、Ⅲ、Ⅰ、Ⅰ、Ⅲ,與專家評(píng)價(jià)結(jié)果一致,組合賦權(quán)-云模型可快速評(píng)價(jià)河道岸坡的穩(wěn)定性現(xiàn)狀,且準(zhǔn)確性較高,具有工程實(shí)踐意義。