王延貴，齊梅蘭，金亞昆

（1.國際泥沙研究培訓(xùn)中心，北京 100048；2.北京交通大學(xué)土木建筑工程學(xué)院，北京 100044）

河道岸灘穩(wěn)定性綜合評價方法

王延貴1，齊梅蘭2，金亞昆2

（1.國際泥沙研究培訓(xùn)中心，北京 100048；2.北京交通大學(xué)土木建筑工程學(xué)院，北京 100044）

為了評價河岸邊界條件和河流動力條件對岸灘穩(wěn)定性的影響，以河道岸灘崩塌影響因子層次結(jié)構(gòu)模型和權(quán)重系數(shù)為基礎(chǔ)，將定性指標(biāo)定量化和定量指標(biāo)標(biāo)準(zhǔn)化，構(gòu)造了枯水期和洪水期岸灘穩(wěn)定綜合評價函數(shù)，提出了考慮多種因素影響的岸灘穩(wěn)定性綜合評價方法。對擬定的3種不同土壤性質(zhì)、邊界條件和水動力條件的河段岸灘穩(wěn)定性評價結(jié)果表明，岸灘穩(wěn)定性綜合評價過程和評價方法是合理的。

河道；岸灘穩(wěn)定性；崩岸；層次分析法；綜合評價方法

在天然沖積河流演變過程中，崩岸是一種普遍存在的自然現(xiàn)象，國內(nèi)外主要河流都存在嚴(yán)重的崩岸問題［1-5］，如中國的長江中下游和黃河下游河道、美國的Ohio河、西非的Niger河和澳大利亞的Lockyer Valley。崩岸往往會給河流防洪、岸邊建筑物安全等造成嚴(yán)重威脅，特別是近些年來隨著大型水利工程的建設(shè)或運行，河勢將發(fā)生明顯變化，崩岸強(qiáng)度會有所增強(qiáng)，如三峽水庫蓄水運用以來，長江中下游河道發(fā)生崩岸655次，崩岸強(qiáng)度有增大的趨勢［5］。因此，開展河流岸灘穩(wěn)定性監(jiān)測與綜合評估就顯得極為重要，通過監(jiān)測和評估可以及時掌控河流岸灘安全狀態(tài)的變化，從而采取相應(yīng)的防護(hù)措施。

不同河段岸灘水流邊界條件存在差別，對應(yīng)的岸灘穩(wěn)定性和崩塌機(jī)理也不一樣。有許多學(xué)者就岸灘穩(wěn)定性和崩塌機(jī)理進(jìn)行了研究［3-4，7-8］，還有一些學(xué)者利用數(shù)學(xué)方法對岸灘安全性進(jìn)行了分析［9-11］。王坤［9］運用模糊決策理論對不同斷面崩岸危險程度進(jìn)行了評價；廖小永等［10］采用模糊統(tǒng)計聚類理論，建立了以岸灘安全程度為聚類標(biāo)準(zhǔn)的隸屬函數(shù)，判斷得到石首河段不同斷面岸灘的安全程度；筆者［11-12］也利用層次分析法對崩岸主要影響因素的作用與權(quán)重進(jìn)行了分析。本文在前期研究的基礎(chǔ)上，建立河道岸灘安全性評估的層次結(jié)構(gòu)模型，考慮專家建議推求各影響指標(biāo)的權(quán)重，進(jìn)而通過模糊數(shù)學(xué)方法將影響指標(biāo)進(jìn)行量化處理，構(gòu)造岸灘穩(wěn)定綜合評價函數(shù)，并對擬定的不同河段穩(wěn)定性進(jìn)行評估。

1　岸灘穩(wěn)定性層次結(jié)構(gòu)模型與影響指標(biāo)權(quán)重

1.1岸灘穩(wěn)定性的層次結(jié)構(gòu)模型

對河道岸灘安全性進(jìn)行綜合評價時，選取岸灘穩(wěn)定性作為目標(biāo)，列為第一層，即目標(biāo)層。河流崩岸是來水來沙條件和河岸邊界條件等諸多因素共同作用的結(jié)果，故將崩岸影響因素分為河岸邊界條件和河流動力條件兩大類，列為第二層。河岸邊界條件包括土壤特性和岸坡形態(tài)，主要的影響因子有土壤黏聚力、內(nèi)摩擦角、干密度及邊坡形態(tài)。河流動力條件主要為來水來沙條件，而來水來沙條件對河道的作用具體表現(xiàn)為河道沖刷，包括河岸側(cè)向沖刷和河床垂向沖刷兩個方面。另外，河流的水位變化及洪水浸泡也對崩岸有重要的影響，而洪水浸泡主要是通過對河岸浸泡時間長短來影響崩岸。因此，河岸側(cè)向沖刷、河床垂向沖刷、水位變化、河岸浸泡時間4個因子作為河流動力條件中的分析對象。層次分析結(jié)構(gòu)模型如圖1所示。

圖1　岸灘穩(wěn)定性評價層次結(jié)構(gòu)模型

1.2不同時期岸灘穩(wěn)定性影響因子權(quán)重

由于河道水位高低變化與岸灘浸泡狀況有很大的差異，河岸穩(wěn)定性也將隨之發(fā)生變化。在汛期，河道水位較高，且變化幅度較大，對岸灘浸泡作用明顯，河道岸灘崩岸往往會發(fā)生在汛末水位迅速降低過程中或之后。長江中下游崩岸90%以上發(fā)生在枯水期或汛后，美國密西西比河下游大多數(shù)崩岸也發(fā)生在汛后，尤以大水年之后崩岸明顯。在枯水期，河道水位較低，河道側(cè)向沖刷明顯，特別是彎道凹岸崩塌嚴(yán)重。顯然，枯水期和洪水期河流岸灘崩塌過程中各影響因子的作用是不一樣的，相應(yīng)的權(quán)重也隨之發(fā)生改變，進(jìn)而影響岸灘安全性的綜合評價。因此，利用層次分析法分別對洪水期和枯水期河流崩岸的影響因子進(jìn)行分析，對各層次中諸影響因子兩兩比較，考慮它們之間的相對重要性，列出判斷矩陣，進(jìn)行岸灘安全性影響分析，求得枯水期和洪水期條件下岸灘穩(wěn)定性各影響因子的權(quán)重系數(shù)［11-14］，如表1所示，相應(yīng)也就得到對應(yīng)枯水期和洪水期各影響因子的權(quán)重矩陣W1、W2：

表1　河道岸灘穩(wěn)定性影響因子權(quán)重［12］

2　岸灘穩(wěn)定性影響因子的量化分析

2.1主要影響因子的量化分析方法

由于岸灘崩塌影響因子的量綱是不同的，很難有效地利用多目標(biāo)函數(shù)評價不同條件下岸灘的穩(wěn)定性，因此，對岸灘影響指標(biāo)進(jìn)行無量綱化處理是非常必要的。河道岸灘安全性影響指標(biāo)可分為定量指標(biāo)和定性指標(biāo)，如與土壤特性相關(guān)的黏聚力、內(nèi)摩擦角、干密度等皆為定量指標(biāo)，具有準(zhǔn)確的數(shù)值；河道沖刷（側(cè)向、垂向）、水位變化和河岸浸泡時間等雖然為定量指標(biāo)，但需要給出進(jìn)一步的說明。如：側(cè)向沖刷寬度為某河段在某一時期內(nèi)側(cè)向沖刷的距離，垂向沖刷深度為某河段在某一時期內(nèi)河床垂向沖刷的深度；水位變化則是某河段在某一時期內(nèi)水位變化的幅度；河岸浸泡時間為某河段在某一時期內(nèi)河岸持續(xù)浸泡的時長。

對于定量指標(biāo)，某些指標(biāo)值越大，岸灘越安全，為正向指標(biāo)，如土壤黏性系數(shù)、土壤內(nèi)摩擦角；反之，為負(fù)向指標(biāo)，如土壤干密度、河岸側(cè)向沖刷、河床垂向沖刷、河道水位變化、河岸浸泡時間。為便于分析，本文統(tǒng)一約定指標(biāo)值大者為優(yōu)（即正向指標(biāo)），對指標(biāo)值小為優(yōu)者（即負(fù)向指標(biāo)）可利用指標(biāo)同趨勢化方法進(jìn)行處理，形成統(tǒng)一的正向指標(biāo)。所謂指標(biāo)同趨勢化就是將指標(biāo)通過整理變換［15］，使所有指標(biāo)轉(zhuǎn)化為同一方向，具體可采用指標(biāo)轉(zhuǎn)置的方法處理，計算公式為

式中：E′i為處理后的正向指標(biāo)值；Emax為負(fù)向指標(biāo)最大值；Emin為負(fù)向指標(biāo)最小值；Ei為原負(fù)向指標(biāo)值。

指標(biāo)進(jìn)行同趨勢化處理后，為了消除各指標(biāo)量綱的影響，需要對指標(biāo)進(jìn)行規(guī)范化或無量綱化處理，具體方法就是用同一指標(biāo)數(shù)列中的最大值去除數(shù)列中的每一個指標(biāo)值，即

式中：Eii為處理后的無量綱指標(biāo)值；Emax為正向指標(biāo)序列中的最大值；Ei為正向序列中第i個指標(biāo)。

岸灘邊坡形態(tài)作為一個定性指標(biāo)，只有定性描述或簡單的等級劃分，而沒有定量數(shù)值。為了進(jìn)一步探討河岸形態(tài)對河岸穩(wěn)定的影響，通過分析不同類型岸灘邊坡形態(tài)（圖2）的穩(wěn)定狀況，引入岸灘邊坡形態(tài)坡度參數(shù)［16］：

式中：K為岸灘邊坡形態(tài)坡度；K0為直線岸坡的坡度；Δu、Δd分別為岸坡上部凸出或凹進(jìn)的厚度（凹進(jìn)取負(fù)值），即圖2中的EF和HI；Su、Sd分別為岸坡上部凸出或凹進(jìn)部分面積（凹進(jìn)取負(fù)值）；L為直線岸坡的坡長，即圖2中的AB段；a為岸坡坡度有效系數(shù)。一般說來，岸灘邊坡形態(tài)坡度小于穩(wěn)定坡度時，岸灘處于穩(wěn)定狀態(tài)；當(dāng)岸灘邊坡形態(tài)坡度介于崩塌坡度和穩(wěn)定坡度之間時，岸灘處于不穩(wěn)定狀態(tài)；當(dāng)岸灘邊坡形態(tài)坡度大于崩塌坡度時，岸灘崩塌。

圖2　岸灘邊坡形態(tài)類型

岸灘邊坡形態(tài)這類指標(biāo)現(xiàn)階段還很難用準(zhǔn)確的數(shù)值表示，仍通過模糊評判法對定性指標(biāo)進(jìn)行量化處理，具體方法是根據(jù)相關(guān)專家對所給定的指標(biāo)按規(guī)定的評語進(jìn)行評判，并賦予指標(biāo)隸屬度。指標(biāo)評語集V為

賦予指標(biāo)評語的隸屬度集U為

根據(jù)專家對各個指標(biāo)不同方案下所給的評語，按其標(biāo)準(zhǔn)隸屬度進(jìn)行平均，取平均值作為該指標(biāo)的隸屬度。

2.2岸灘穩(wěn)定性影響因子的指標(biāo)量化

為便于分析比較和實際應(yīng)用，結(jié)合河流岸灘邊坡形態(tài)、土壤特性和水流條件的實際情況，設(shè)計了河段土體性質(zhì)和水位變化差異較大的3個不同河段進(jìn)行穩(wěn)定性評估，3個河段分別稱為河段1、河段2和河段3。擬定3種河段在不同時期河岸土壤性質(zhì)、河流動力條件和邊界條件所對應(yīng)的指標(biāo)，如表2所示。

在不同河段河道岸灘安全性綜合評價過程中，通過對定性指標(biāo)和定量指標(biāo)進(jìn)行量化處理，得到層次結(jié)構(gòu)模型中底層各影響因子指標(biāo)的量化值：

式中X1、X2分別為枯水期和洪水期不同河段河道岸灘各影響因子的量化指標(biāo)矩陣。

表2　3個河段枯水期和洪水期影響指標(biāo)

3　河道岸灘穩(wěn)定性綜合評價與分析

河道岸灘穩(wěn)定性綜合評價就是通過岸灘穩(wěn)定性評價層次結(jié)構(gòu)模型和權(quán)重系數(shù)，利用模糊數(shù)學(xué)方法將定性指標(biāo)定量化和定量指標(biāo)標(biāo)準(zhǔn)化，進(jìn)而構(gòu)造枯水期和洪水期岸灘穩(wěn)定綜合評價函數(shù)，以評價河道岸灘穩(wěn)定性。

在確定了岸灘底層各影響指標(biāo)的量化值后，根據(jù)各影響因子指標(biāo)的權(quán)重系數(shù)，可以構(gòu)造河道岸灘穩(wěn)定綜合評價函數(shù)，即

式中：G為河道岸灘穩(wěn)定綜合評價指標(biāo)函數(shù)值，G越大，河段穩(wěn)定性越好；xi為底層各影響因子的標(biāo)準(zhǔn)量化值；wi為底層各影響因子對岸灘安全性作用的權(quán)重；n為底層影響因子的數(shù)量。

對于枯水期和洪水期的岸灘，已知不同時期各影響因子的權(quán)重矩陣W1和W2，以及不同河段影響因子的量化指標(biāo)矩陣X1和X2，便可求得枯水期和洪水期河段的河道岸灘穩(wěn)定綜合評價函數(shù)Gk和Gh，即

利用公式（11）和（12）計算枯水期和洪水期3種河段的河道岸灘穩(wěn)定綜合評價函數(shù)值，結(jié)果如表3所示。

表3　不同河段河道岸灘穩(wěn)定綜合評價函數(shù)值

結(jié)合表1和表2的實際情況，從表3可以得出：

a.枯水期河段1、河段2和河段3的岸灘穩(wěn)定綜合評價函數(shù)值分別為0.830、0.695和0.620，表明河段1最安全，河段2次之，河段3最不安全。3個河段的土壤性質(zhì)差別較大，河段1土壤性質(zhì)對應(yīng)的主要影響因子對岸灘的穩(wěn)定性最有利，河段3土壤性質(zhì)對應(yīng)的影響因子對岸灘穩(wěn)定性不利，雖然河段1的沖刷較為嚴(yán)重，其岸灘穩(wěn)定性排序仍為第一位。

b.洪水期河段1、河段2和河段3的岸灘穩(wěn)定綜合評價函數(shù)值分別為0.934、0.736和0.548，表明河段1最安全，河段2次之，河段3最不安全。同樣，雖然河段1的沖刷、水位變化和河岸浸泡時間都較大，但仍不能超過河岸土壤性質(zhì)的作用，因此河段1的安全性仍然最好；河段3由于按土壤性質(zhì)的安全性較弱，再加上水位變化和河岸浸泡時間的影響都較大，使河段3的穩(wěn)定性最差。

在利用該方法評價河道岸灘穩(wěn)定性時，需要獲得河岸邊界條件和河流動力條件中的8個影響因子的指標(biāo)值。在實際河流中，多數(shù)情況沒有如此全面的岸灘實測資料，獲得諸多實測資料仍有很大的困難。因此，根據(jù)河道岸灘邊界水流動力特點，本文設(shè)計了3種河道岸灘，并擬定了相應(yīng)影響因子的指標(biāo)值，目的是介紹河道岸灘穩(wěn)定性評價的步驟和過程。目前，河道岸灘穩(wěn)定性綜合評價方法還沒有針對實際河流岸灘利用實測資料進(jìn)行的評價和檢驗，也就是說，該方法若要應(yīng)用于實際，還有許多工作需要深入研究，這也是我們努力的方向之一。

4　結(jié) 語

結(jié)合不同時期影響河道岸灘穩(wěn)定性的主要因素及其差異，利用層次分析法研究了枯水期和洪水期各主要影響因子的作用，根據(jù)岸灘穩(wěn)定性影響因子權(quán)重和無量綱化指標(biāo)的特點，構(gòu)造了枯水期和洪水期河道岸灘穩(wěn)定綜合評價函數(shù)，并對設(shè)計的3種不同岸灘河段的穩(wěn)定性進(jìn)行了評價。鑒于河道岸灘穩(wěn)定性綜合評價方法還沒有應(yīng)用于實際的河道岸灘穩(wěn)定性評價，以后還有許多工作需要深入研究。

［1］長江水利委員會.長江中下游護(hù)岸工程40年［C］//長江中下游護(hù)岸工程論文集（4）.武漢：水利部長江水利委員會，1990：15.

［2］趙業(yè)安，周文浩，費祥俊，等.黃河下游河道演變基本規(guī)律［M］.鄭州：黃河水利出版社，1998.

［3］ABAM T K S.Factors affecting distribution of instability of river banks in the Niger delta［J］.Engineering Geology，1993，35：123-133.

［4］HAGERTY D J，SPOOR M F，ULLRICH C R.Bank failure and erosion on the Ohio River［J］.Engineering Geology，1981，17：141-158.

［5］THOMPSON C，CROKE J，GROVE J，et al.Spatiotemporal changes in river bank mass failures in the LockyerValleyQueensland，Australia［J］. Geomorphology，2013，19（1）：129-141.

［6］中華人民共和國水利部.中國河流泥沙公報2011［M］.北京：中國水利水電出版社，2012.

［7］王延貴.河流岸灘挫落崩塌機(jī)理及其分析模式［J］.水利水電科技進(jìn)展，2013，33（5）：21-25.（WANG Yangui. Study on mechanism and analysis mode of sinking failure for theriverbank［J］.AdvancesinScienceand Technology of Water Resources，2013，33（5）：21-25.（in Chinese））

［8］張芳枝，陳曉平.河流沖刷作用下堤岸穩(wěn)定性研究進(jìn)展［J］.水利水電科技進(jìn)展，2009，29（4）：84-88.（ZHANG Fangzhi，CHEN Xiaoping.Research progress on stability analysis of embankment under effects of river scouring［J］.Advances in Science and Technology of Water Resources，2009，29（4）：84-88.（in Chinese））

［9］王坤.模糊決策理論在崩岸研究中的應(yīng)用［J］.長江科學(xué)院院報，2005，22（8）：12-15.（WANGKun. Application of fuzzy decision-making theory in bank-failure study［J］.Journal of Yangtze River Scientific Research Institute，2005，22（8）：12-15.（in Chinese））

［10］廖小永，王坤.模糊統(tǒng)計聚類理論在崩岸問題中的應(yīng)用研究［J］.長江科學(xué)院院報，2007，24（4）：5-8.（LIAO Xiaoyong，WANG Kun.Study on application of fuzzy statistic clustering theory in bank-failure［J］.Journal of Yangtze River Scientific Research Institute，2007，24（4）：5-8.（in Chinese））

［11］WANG Yangui，JIN Yakun，SHI Hongling.Study on roles of influence factors in river bank failure［C］//Proceedings of35thIAHRWorldCongress.Beijing：Tsinghua University Press，2013：213.

［12］王延貴，金亞昆.模糊層次分析在河道岸灘穩(wěn)定性評價中的應(yīng)用［J］.浙江水利科技，2014，42（5）：38-41.（WANG Yangui，JIN Yakun.Application of analytic hierarchy process in evaluation of river bank stability［J］. Zhejiang Hydrotechnics，2014，42（5）：38-41.（in Chinese））

［13］SAATY T L.The analytical hierarchy process［M］.New York：McGraw Hill，1980.

［14］汪浩，馬達(dá).層次分析標(biāo)度評價與新標(biāo)度方法［J］.系統(tǒng)工程理論與實踐，1993（9）：25-26.（WANG Hao，MA Da.Scale evaluation and new scale methods［J］.Systems Engineering-Theory&Practice，1993（9）：25-26.（in Chinese））

［15］張慶華，白玉慧，倪紅珍.節(jié)水灌溉方式的優(yōu)化選擇［J］.水利學(xué)報，2002，33（1）：47-51.（ZHANG Qinghua，BAI Yuhui，NI Hongzhen.Optimal choice of water saving irrigation mode［J］.Journal of Hydraulic Engineering，2002，33（1）：47-51.（in Chinese））

［16］王延貴，匡尚富.河岸崩塌類型與崩塌模式的研究［J］.泥沙研究，2014（1）：13-20.（WANG Yangui，KUANG Shangfu.Study on types and collapse modes of the bank failures［J］.Joural of Sediment Research，2014（1）：13-20.（in Chinese））

Comprehensive evaluation method for river bank stability

WANG Yangui1，QI Meilan2，JIN Yakun2
（1.International Research and Training Center on Erosion and Sedimentation（IRTCES），Beijing 100048，China；2.School of Civil Engineering，Beijing Jiaotong University，Beijing 100044，China）

In order to evaluate the influences of the bank boundary condition and the river dynamic condition on bank stability，through quantification of qualitative indices and standardization of quantitative indices，comprehensive evaluation functions for bank stability during the dry period and the flood period were created based on the hierarchical framework model and weight coefficients of the factors of bank failures，and a comprehensive evaluation method for bank stability，considering many factors，was put forward.The results of stability evaluation of three reach banks with different soil properties，bank boundary conditions，and river dynamic conditions show that the comprehensive evaluation function and method are suitable for bank stability evaluation.

river；bank stability；bank failure；analytic hierarchy process；comprehensive evaluation method

TV143+.3

A

10067647（2016）05005505

10.3880/j.issn.10067647.2016.05.010

國家自然科學(xué)基金（511792080）；中國水利水電科學(xué)研究院科研專項（沙基本科研1438）；“十二五”國家科技支撐計劃（2012BAB04B02-2）

王延貴（1963—），男，教授級高級工程師，博士，主要從事水力學(xué)及河流動力學(xué)研究。E-mail：wangyg@iwhr.com

（20150628 編輯：熊水斌）