山區(qū)庫岸塌岸預(yù)測方法綜述

張 梁，王俊杰，，閻宗嶺

(1.重慶交通大學(xué)水利水運工程教育部重點實驗室，重慶400074;2.國家山區(qū)公路工程技術(shù)研究中心，重慶400067)

隨著我國經(jīng)濟的快速發(fā)展和西部大開發(fā)戰(zhàn)略的全面實施，許多水庫正處于論證、設(shè)計或建設(shè)階段，大量內(nèi)河港口、碼頭正在或?qū)⒁ㄔO(shè)，其中大多數(shù)位于水電資源非常豐富的西南、中部山區(qū)，由此而引起的庫區(qū)庫岸再造及塌岸情況的預(yù)測就顯得尤為重要。

水庫蓄水及運行過程中，庫岸所處的地質(zhì)環(huán)境將發(fā)生顯著改變，自然平衡條件遭到破壞，從而引起岸坡變形失穩(wěn)，庫岸線也逐漸后退，直至達到新的平衡狀態(tài)為止，這一過程稱為庫岸再造。庫岸再造是一個十分復(fù)雜的動力地質(zhì)過程，受岸坡物質(zhì)組成、結(jié)構(gòu)特征、形態(tài)及水流等多因素控制，塌岸過程復(fù)雜，尚無法精確地通過數(shù)學(xué)計算式來表達。

目前，庫岸塌岸預(yù)測方法有以佐洛塔廖夫［1－2］為代表的條件類比圖解法、以康德拉捷夫為代表的數(shù)學(xué)分析法［3］、平衡剖面法［4］、動力法［5］、徐瑞春的塌岸預(yù)測圖解的若干修正［6］、超前信息法等。常用的方法以圖解法為主，包括卡丘金法［7－8］、佐洛塔廖夫法［1－2］、兩段法［9］、岸坡結(jié)構(gòu)法［10］4 種。

范云沖等［7］在系統(tǒng)分析庫岸再造相關(guān)因素的基礎(chǔ)上，簡要地介紹了水庫塌岸寬度預(yù)測的不同方法，并對這些方法進行了評價。馬淑芝等［4］在劃分地質(zhì)結(jié)構(gòu)單元的基礎(chǔ)上，用穩(wěn)態(tài)坡形類比法對巫山新城庫岸再造進行預(yù)測。王躍敏等［9］通過外福鐵路線水口水庫庫岸的塌岸觀測研究，提出適合我國南方山區(qū)峽谷型水庫塌岸的兩段預(yù)測法。李彥軍［11］從巖土統(tǒng)計參數(shù)分析入手，引入蒙特卡羅法在庫岸預(yù)測參數(shù)方面的應(yīng)用，使巖土體參數(shù)的獲取簡便和可操作性，具有一定的實用價值。劉天翔［10］、湯明高［12－13］、許強等［14－15］對三峽庫區(qū)庫岸劃分了不同塌岸類型，可用岸坡結(jié)構(gòu)法對沖(磨)蝕型和坍(崩)塌型塌岸進行預(yù)測，用極限平衡搜索預(yù)測法和Flac3D數(shù)值模擬預(yù)測法對滑移型的塌岸進行預(yù)測。

由于我國的地質(zhì)地貌的關(guān)系及國家經(jīng)濟發(fā)展的需要，大多數(shù)水庫修建在西南、中部山區(qū)，選擇合適的預(yù)測方法迫在眉睫，筆者從探討山區(qū)水庫塌岸過程、影響因素兩方面出發(fā)，對各種水庫塌岸預(yù)測方法進行比較，分析各方法的適用性。

1 山區(qū)庫岸塌岸過程及影響因素

1.1 山區(qū)庫岸塌岸過程

水庫蓄水后，庫岸巖土體在地表水及地下水的浸泡和波浪的磨蝕作用下，發(fā)生軟化、崩解等現(xiàn)象，庫岸巖土體開始坍塌，使岸線逐漸向邊岸方向后移，促使磨蝕淺灘不斷的增長，波浪的回流則把破壞的物質(zhì)搬運、分選并堆積下來，促使堆積淺灘不斷的向水庫方向增長。水位變幅帶的岸坡再造是一個磨蝕、搬運、堆積綜合作用的過程(圖1)［16］。由于堆積淺灘的形成使近岸的坡角變緩，增加了波浪的破碎距離，磨蝕作用也相應(yīng)的減弱。同時由于堆積淺灘的形成，使庫岸再造的起點向水庫方向移動，造成水庫再造線向庫內(nèi)的平移，從而也減小了庫岸的再造量。

圖1 土質(zhì)岸坡坍塌過程Fig.1 Collapse process of soil bank failure

對于典型的高山峽谷河流，其岸坡陡峻，在水位變幅帶內(nèi)很難形成堆積部分，其淺灘部分全部由磨蝕淺灘組成。在此種情況下岸坡的磨蝕作用較形成堆積灘時要劇烈，更有利于岸坡的破壞。經(jīng)過上述過程的重復(fù)作用，當(dāng)岸坡的坡度和其相應(yīng)的長度達到某一特定值時，在現(xiàn)有水文、水動力條件不變的前提下，岸坡的磨蝕作用力與其抵抗作用的能力相平衡時，岸坡的再造進入相對穩(wěn)定狀態(tài)，水庫塌岸也就進入平衡剖面的形成階段。

平衡剖面中這種跟一定的地層巖性條件與水文條件相適應(yīng)的淺灘極限坡角稱為最終淺灘坡角。平衡剖面的結(jié)構(gòu)由3部分形成:①水上岸坡;②淺灘;③淺灘前緣斜坡。

1.2 山區(qū)庫岸塌岸影響因素

對于山區(qū)水庫而言，起主導(dǎo)作用的是現(xiàn)代地質(zhì)作用，岸坡的破壞主要體現(xiàn)在巖土體在庫水作用下物理力學(xué)性質(zhì)的變化，靜止庫水位時的浮托力，水位升降時的滲流力及地下水的潛蝕、溶蝕作用等因素［17－22］對岸坡穩(wěn)定性的影響。

導(dǎo)致水庫塌岸的因素可歸為內(nèi)因和外因［21－29］。地形地貌、地層巖性及地質(zhì)構(gòu)造條件、風(fēng)化程度等是促使庫岸再造的內(nèi)因;庫水位升降、水文地質(zhì)條件、波浪作用、水位變動作用、地震力作用、降雨、人類活動等是影響庫岸再造的外因。內(nèi)因是庫岸再造的根本原因，外因是庫岸再造的重要影響因素，庫岸再造過程是各種內(nèi)因和外因共同作用的結(jié)果。

2 山區(qū)土質(zhì)庫岸塌岸預(yù)測方法

水庫塌岸預(yù)測理論最早來源于前蘇聯(lián)。1935年，前蘇聯(lián)學(xué)者薩瓦連斯基首次提出了水庫塌岸問題，認為水庫塌岸是由于河道壅水后，水位抬高，吃水線與岸坡相接觸，岸壁的天然平衡條件遭到破壞而引起的。同時指出波浪是造成水庫塌岸的主要因素之一。薩瓦連斯基開啟了水庫塌岸預(yù)測研究的先河，從而帶動了一大批學(xué)者進行這方面的研究，也因此而產(chǎn)生了一系列對現(xiàn)今塌岸預(yù)測研究仍具有深遠意義的成果。

1937年，什利亞莫夫開始研究水庫塌岸問題，他認為水庫塌岸的特征決定于3組因素:①水文因素;②地質(zhì)地貌因素;③其他因素。

在20世紀40～50年代，前蘇聯(lián)薩瓦連斯基、卡丘金、佐洛塔寥夫等研究了蘇聯(lián)的水庫塌岸問題，提出的塌岸范圍預(yù)測的基本計算方法和圖解法。近年來，有關(guān)水庫塌岸預(yù)測方法的研究仍在進行?，F(xiàn)階段，預(yù)測水庫塌岸或水庫邊岸再造范圍和規(guī)模的方法可分為如下幾種。

2.1 圖解法

利用現(xiàn)階段不同巖土體水下穩(wěn)定坡角、水位變幅帶坡角和水上穩(wěn)定坡角，與將來水庫蓄水后不同庫水條件下的庫岸岸坡類比，通過各土層不同的穩(wěn)定坡角繪制庫岸剖面線來預(yù)測庫岸塌岸距離;根據(jù)水庫現(xiàn)今的庫岸岸坡剖面，通過一定的參數(shù)來計算并繪制待預(yù)測庫岸的塌岸剖面從而進行塌岸預(yù)測。

目前，主要的圖解預(yù)測方法包括工程地質(zhì)類比法、卡丘金法和卓洛塔寥夫法3種。

2.1.1 工程地質(zhì)類比圖解法［5－6，30］

由于現(xiàn)階段天然河道的平均枯水位、江水漲幅帶、平均洪水位分別與水庫運行期低水位、調(diào)節(jié)水位(即水位變動帶)、最高設(shè)計水位存在可類比性，因此可以通過地質(zhì)調(diào)查，并統(tǒng)計現(xiàn)天然河道的平均枯水位以下、江水漲幅帶以及平均洪水位以上三帶內(nèi)不同巖土體的穩(wěn)態(tài)坡角。以此作為該巖土層在不同庫水位條件下的穩(wěn)定坡角，進而類比圖解水庫蓄水運行時的庫岸再造范圍。根據(jù)實測岸坡剖面，自現(xiàn)江河枯水位起，首尾相連依次繪出在不同庫水位條件下相應(yīng)巖土層的穩(wěn)定坡角，并以各段穩(wěn)定坡角連線代表最終庫岸再造邊界線，進而量取庫岸再造的最終寬度與高程。

圖解中根據(jù)地質(zhì)測繪與勘探資料，現(xiàn)場調(diào)查統(tǒng)計不同巖土體在天然河道的平均枯水位以下、江水漲幅帶以及平均洪水位以上三帶內(nèi)巖土體的穩(wěn)定坡角。采用調(diào)查統(tǒng)計的數(shù)據(jù)，按公式(1)計算各類巖土層在不同庫水條件下的穩(wěn)定坡角。

式中:α為一個統(tǒng)計范圍內(nèi)該巖土層的穩(wěn)定坡角;αi為單個統(tǒng)計點該巖土層的坡角;Li為單個統(tǒng)計順坡向之間的平面距離。

因枯水位以下巖土層穩(wěn)態(tài)坡角通常無法量取，可將江水漲幅帶穩(wěn)態(tài)坡角按0.8的系數(shù)折減而得。

然后根據(jù)各巖土層自然岸坡坡度統(tǒng)計值與前述類比原則，得出各巖土層在不同庫水位狀態(tài)下的穩(wěn)定坡角建議值，最后采用圖解法求得塌岸范圍。

類比圖解法是一種普遍適用的方法，可直接用于均質(zhì)土質(zhì)岸坡，但不能用于預(yù)測岸坡結(jié)構(gòu)復(fù)雜的庫岸岸坡。

2.1.2 卡丘金法［7－8，31］

E.r.卡丘金在前人研究的資料的基礎(chǔ)上，根據(jù)多年的實測資料于1949年提出岸坡最終塌岸預(yù)測計算公式［3］:式中:S為最終塌岸寬度，m;N為堆積系數(shù)(與土石顆粒成分有關(guān)的系數(shù)，土石顆粒愈粗，就愈加易于形成水下堆積岸坡，所以按卡丘金提供的經(jīng)驗數(shù)據(jù)，砂土的N值為0.5，亞黏土的為0.6，黏土的為1.0;當(dāng)原始岸坡較陡，庫水水深較大時，難于形成水下堆積階地，此時N實際應(yīng)等于1);A為水位變化幅度，設(shè)計高水位與設(shè)計低水位差值，m;hp為波浪影響深度，m(設(shè)計低水位以下波浪影響深度一般取1～2倍浪高:如果浪高取0.5 m時，浪高影響深度取1 m);hσ為正常高水位以上岸坡的高度，m;α為最終磨蝕淺灘坡角;β為水上岸坡的穩(wěn)定坡角;γ為原始岸坡坡角;hb為波浪爬升高度，設(shè)計高水位以上浪爬高度可按式(3)計算:

式中:k為被沖蝕的岸坡表面糙度系數(shù)，一般砂質(zhì)岸坡k取0.55～0.75，礫石質(zhì)岸坡k=0.85～0.9，混凝土k=1，拋石k=0.775;坡度α可參照河谷邊岸平水位處河濱淺灘坡角值，當(dāng)己知作用于該岸坡地帶波浪波高和組成岸坡土石顆粒成分時，也可根據(jù)各種顆粒成分沉積物的水下岸坡坡度與浪高的關(guān)系圖解確定;h為波浪波高，根據(jù)各庫區(qū)實際情況取值。

卡丘金法的實質(zhì)是依據(jù)實測的洪枯水位變幅帶各類巖性岸坡長期穩(wěn)定坡角，根據(jù)幾何關(guān)系用圖解法求解岸坡最終塌岸預(yù)測寬度，其精度取決于計算參數(shù)的選定?？紤]到計算參數(shù)大多選自經(jīng)驗值，因此，在實際預(yù)測時必須對類似水動力條件和類似巖土體條件下的已有岸坡塌岸進行觀測，以獲得相應(yīng)的較為可靠的計算參數(shù)。有時參考水庫蓄水前的洪、枯水位變幅帶岸坡形態(tài)數(shù)據(jù)來計算，也具有較好的預(yù)測效果。

該方法適用于小型水庫以及由松散沉積層——黃土、砂、砂質(zhì)黏土與黏土所覆蓋的具有不高的岸坡地帶及平原地區(qū)水庫的塌岸預(yù)測;在山區(qū)河道型水庫的塌岸預(yù)測中，適用于水庫的中游和上游地帶，而預(yù)測水庫下游的結(jié)果往往與實際相差甚遠，實際塌岸寬度要比預(yù)測的小得多，卡丘金法的預(yù)測結(jié)果偏于安全。

2.1.3 佐洛塔寥夫法［1－2，32］

此法為前蘇聯(lián)學(xué)者佐洛塔寥夫于1955年提出，通過圖解法進行岸坡最終塌岸寬度預(yù)測。其原理為:對于水庫庫岸再造，波浪作用是主要的，庫岸再造后的岸坡可分為淺灘外緣陡坡、堆積淺灘、沖蝕淺灘、爬升帶斜坡以及水上岸坡帶五段，通過作圖得到上述5段岸坡，即為庫岸再造的最終岸坡。

對于松散堆積岸坡(殘坡積、崩坡積、滑坡堆積以及人工棄渣岸坡)，大型水庫的中、下游地段，一般采用佐洛塔寥夫提出的圖解法。該法認為:水庫中下游地段，水深較大，水面較廣，波高增加，對庫岸的破壞，波浪作用是主要的。

佐洛塔寥夫圖解法考慮了水下堆積淺灘的影響，將下部起點設(shè)在堆積淺灘臺坎前緣，并且假設(shè)以該點的深度在正常高水位以下浪高h處作為預(yù)測10 a期間的坍岸剖面，以該點深度位于保證率5%枯水位以下3 h(黏土類邊岸)或2 h(砂土邊岸)深處作為確定最終坍岸剖面的起點。其作圖方法與卡丘金方法相同，并且還分別考慮了下游寬闊部位和上游狹窄部位兩種情況，這種方法在實際運用中必須查明有多少比例的沖蝕土可組成堆積淺灘，因而實際運用較為復(fù)雜。

2.2 以康德拉捷夫為代表的數(shù)學(xué)分析法

康德拉捷夫［3］在1953年提出了一種塌岸寬度的計算分析方法，與卡丘金和卓洛塔廖夫的方法不同的是他認為淺灘表面的外形輪廓線不是直線和折線型，而是拋物線型。這種方法一般適用于巖性條件單一，由砂土和粉砂土等非黏性上組成的非滑動斜坡，并都不考慮泥沙縱向移動的情況。他建議用下面的公式計算任一時間的塌岸寬度:

式中:l為水庫蓄水任一年的塌岸寬度，m;Wt為水庫蓄水任一年的塌岸體積，m3;L為水庫蓄水最終的塌岸寬度，m;W0為水庫蓄水最終的塌岸體積，m3;V0為水庫蓄水第1 a的塌岸速度，m/a;t為計算年限，a;Wl為水庫蓄水第t年的塌岸體積，m3。

該塌岸公式用來預(yù)測塌岸寬度，其難點在于各參數(shù)的確定。

2.3 平衡剖面法

在水庫風(fēng)(波)浪和船行波的長期作用下，岸坡斷面將逐漸調(diào)整至平衡位置，形成平衡斷面。因此可以根據(jù)水庫運行性質(zhì)、波浪作用規(guī)律，以及岸坡巖土體工程地質(zhì)特征，運用水力學(xué)、泥沙運動學(xué)等理論以及實際觀察數(shù)據(jù)，可以建立基于經(jīng)驗的數(shù)學(xué)模型，用于預(yù)測此平衡斷面，從而獲得水庫塌岸的空間規(guī)模。

平衡剖面法［4］需要根據(jù)觀測和試驗數(shù)據(jù)，總結(jié)分析出各種水位變幅帶之間、波浪作用帶之間的穩(wěn)定岸坡坡角與波浪要素間的關(guān)系曲線，以供繪制平衡斷面。該方法要求太高，實用性差。

2.4 動力法

動力法［5］的預(yù)測依據(jù)是塌岸量與波能和巖石抗沖刷強度之間的“關(guān)系方程”:

式中:Q為庫岸在單位寬度內(nèi)被沖刷的巖土體的體積，m3/m;E為波浪作用于單位寬度庫岸的動能，kg·m;Kp為巖土體的抗沖刷系數(shù)，m3/t;t為水庫運營年限，a;b為經(jīng)驗常數(shù)，取決于濱岸淺灘中堆積部分寬度，變幅為0.45～0.95。

該法有一定物理依據(jù)，但“關(guān)系方程”的建立同時也需要一定量的觀測樣本。在海洋工程科學(xué)領(lǐng)域該方法已得到一定程度的應(yīng)用，在水庫較少使用。

2.5 關(guān)于紅層塌岸預(yù)測的若干修正

由紅層構(gòu)成的岸坡比由堅硬巖體構(gòu)成的岸坡穩(wěn)定性差，但又遠較由松軟土體構(gòu)成的岸坡穩(wěn)定性好。所以，預(yù)測紅層岸坡垮塌的方法不能照搬固結(jié)巖體或者松散土體岸坡垮塌的預(yù)測方法。

徐瑞春［6］在多年對四川盆地紅層河谷邊坡研究的基礎(chǔ)上，對紅層塌岸預(yù)測方法進行了一些有益的探索，提出了紅層塌岸預(yù)測圖解方法的若干修正。

該法屬于經(jīng)驗圖解法，其有效性有待進一步驗證。

2.6 兩段法

該方法有王躍敏等［9］提出。其原理為:預(yù)測塌岸線由水下穩(wěn)定岸坡線和水上穩(wěn)定岸坡線的連線組成時，水下穩(wěn)定岸坡線由原河道多年最高洪水位及水下穩(wěn)定坡角確定，水上穩(wěn)定岸坡線由設(shè)計洪水位和毛細水上升高度及水上穩(wěn)定坡角確定。

該法適用于我國南方山區(qū)的峽谷型水庫，庫面較窄，風(fēng)浪作用較小，岸坡地層為黏性土、砂性土、碎石類土、棄碴及巖石的全風(fēng)化地層，有較完整的水文氣象資料等?！皟啥畏ā币彩且环N類比法，只是在確定水上穩(wěn)定坡度角和水下穩(wěn)定坡度角所選用的參數(shù)不同。在使用“兩段法”預(yù)測坍岸寬度的同時，還需用卡丘金法進行比較，全面驗證其合理性，科學(xué)地取舍。也就說其適用性還有待于進行大量的應(yīng)用驗證。

2.7 岸坡結(jié)構(gòu)法

許強等［10］通過三峽庫區(qū)數(shù)百段庫岸的塌岸地質(zhì)現(xiàn)場調(diào)查和預(yù)測分析的基礎(chǔ)上，提出了岸坡結(jié)構(gòu)法，并給出了詳細的圖解方法和預(yù)測步驟。

岸坡結(jié)構(gòu)法的主要原理就是根據(jù)岸坡上各種不同物質(zhì)的水下堆積坡角、沖磨蝕角(水位變動帶內(nèi)坡角)、水上穩(wěn)定坡角和水庫的設(shè)計低水位、設(shè)計高水位來進行預(yù)測，也是一種圖解法和類比法。庫岸結(jié)構(gòu)法采用3個坡度，主要是針對水庫調(diào)度，并沒有考慮波浪因素。該法綜合考慮了水庫的各水位之間存在的聯(lián)合作用，適合于三峽水庫這種山區(qū)水庫，水位變化主要源于水庫的調(diào)度。適用于沖磨蝕型和坍塌型塌岸。

2.8 其它方法

除了上述目前較常使用的幾種塌岸預(yù)測方法外，日本京都大學(xué)的 Nagata N等［33］(2000)還開展了運用數(shù)值模擬方法進行了塌岸速率、塌岸范圍預(yù)測的嘗試，但由于塌岸形成機制的復(fù)雜，以及影響塌岸的水文動態(tài)都很難模擬，因此，這方面的研究還處在實驗階段。Elei，Sebnem and Paul A.(2002)等［34］對美國喬治亞州與南卡羅萊納州交界的Savannah河上的Hartwell Lake水庫塌岸進行了研究，并針對黏性土庫岸塌岸的長期觀測，從能量法的角度提出一預(yù)測長期塌岸速率的表達式:

式中:ER為塌岸速率;kn為比例系數(shù);Hb為近岸碎浪的波高。

該方法在Hartwell Lake水庫塌岸的預(yù)測中得到了較好的應(yīng)用，但該方法其實質(zhì)仍然是能量法，其中的塌岸預(yù)測參數(shù)具有地域性的限制，并且還必須要有長期的觀測資料做支撐，故在山區(qū)庫岸塌岸預(yù)測中有一定的時間局限性及處理上的復(fù)雜性。

根據(jù)官廳水庫塌岸預(yù)測，在前蘇聯(lián)學(xué)者提出的方法的基礎(chǔ)上，孫廣忠［35］提出了一種水庫下游區(qū)的塌岸預(yù)測方法。他認為淺灘外形具有指數(shù)曲線性質(zhì)，預(yù)測時關(guān)鍵是要弄清楚在塌岸過程中庫岸剖面外形的變化及堆積系數(shù)，確定庫岸剖面為任一形狀時的位置，從而可以求得塌岸帶的寬度及塌岸量。孫廣忠提出的方法在官廳水庫的塌岸預(yù)測中得到了很好的應(yīng)用。

馬淑芝［4］、范云沖等［7］根據(jù)張咸恭等提出的穩(wěn)定坡形、坡角地質(zhì)類比法提出了穩(wěn)態(tài)坡形類比法，它是根據(jù)工程地質(zhì)類比法的原則，對不同工程性質(zhì)的岸坡選擇不同的穩(wěn)態(tài)坡形，充分考慮了岸坡巖體工程性質(zhì)的差異，使預(yù)測結(jié)果更趨于合理、準確。穩(wěn)態(tài)坡形類比法在三峽庫區(qū)移民遷建城鎮(zhèn)庫岸再造預(yù)測中得到了較好的應(yīng)用，可以進一步推廣到土質(zhì)岸坡庫岸再造的預(yù)測中。

考慮到庫岸邊坡工程的復(fù)雜性和非線性，王征亮等(2005)［36］應(yīng)用神經(jīng)網(wǎng)絡(luò)和可拓學(xué)等方法進行了庫岸邊坡穩(wěn)定性評價和塌岸預(yù)測，將非線性分析方法引入塌岸研究，值得進一步探索。

3 結(jié)語

1)總結(jié)的塌岸預(yù)測方法中，有些在具體的庫岸塌岸預(yù)測實踐中效果較好，對實際有很好的指導(dǎo)意義;有些方法還值得進一步的探索與研究。

2)由于山區(qū)庫岸地質(zhì)條件的復(fù)雜多變，在工程具體應(yīng)用中，仍要堅持理論與實踐相結(jié)合的原則，根據(jù)實際選用合理的方法手段。

3)塌岸的預(yù)測是一個眾所周知的世界性難題，至今為止，雖然取得了令人欣喜的成果，但要想及時精確地預(yù)測，現(xiàn)在還做不到，因地質(zhì)條件的復(fù)雜多變都會對庫岸預(yù)測結(jié)果帶來影響。不同的地質(zhì)條件和塌岸機理，必須選取不同的預(yù)測模型和預(yù)測參數(shù)，進而選擇不同的庫岸塌岸預(yù)測方法。所以對庫岸塌岸的預(yù)測，還期望更多的科技工作者熱切地關(guān)注它和重視它，并繼續(xù)對其進行更深入更系統(tǒng)地研究和探索。

［1］田一德，汪小蓮.三峽水庫庫岸崩滑體處理問題初探［J］.人民長江，1999，30(11):6 －7.

［2］中科院地質(zhì)研究所.水利水電工程地質(zhì)［M］.北京:科學(xué)出版社，1974.

［3］張悼元，王士天，王蘭生.工程地質(zhì)分析原理［M］.北京:地質(zhì)出版社，1997.

［4］馬淑芝，賈洪彪，唐輝明.利用穩(wěn)態(tài)坡形類比法預(yù)測基巖岸坡的庫岸再造［J］.中國地質(zhì)大學(xué)學(xué)報:地球科學(xué)，2002，27(2):231－234.

［5］［前蘇聯(lián)］В.Д.洛姆塔澤.工程動力地質(zhì)學(xué)［M］.李生林，劉蕙蘭，譯.北京:地質(zhì)出版社，1985:3.

［6］徐瑞春.紅層與大壩［M］.武漢:中國地質(zhì)大學(xué)出版社，2003.

［7］范云沖，張友誼，胡卸文.庫岸再造預(yù)測方法及其評價［J］.四川水力發(fā)電，2002，21(4):69 －71，83.

［8］張奇華，丁秀麗，張杰，等.三峽庫區(qū)奉節(jié)河段庫岸蓄水再造研究［J］.巖石力學(xué)與工程學(xué)報，2002 21(7):1007－1012.

［9］王躍敏，唐敬華，凌建明.水庫塌岸預(yù)測方法研究［J］.巖土工程學(xué)報，2000，22(5):569 －571.

［10］劉天翔，許強，黃潤秋，等.三峽庫區(qū)塌岸預(yù)測評價方法初步研究［J］.成都理工大學(xué)學(xué)報:自然科學(xué)版，2006，33(1):77 －83.

［11］李彥軍.蒙特卡羅法在水庫塌岸預(yù)測中的應(yīng)用［J］.水土保持研究，2006，13(6):16 －17.

［12］湯明高，許強，黃潤秋.三峽庫區(qū)典型塌岸模式研究［J］.工程地質(zhì)學(xué)報，2006，14(2):172 －177.

［13］湯明高，許強，黃潤秋，等.三峽庫區(qū)塌岸預(yù)測參數(shù)及其影響因素研究［J］.成都理工大學(xué)學(xué)報:自然科學(xué)版，2006，33(5):460－464.

［14］黃潤秋，許強，李渝生，等.三峽庫區(qū)塌岸模式、塌岸預(yù)測參數(shù)與塌岸范圍預(yù)測評價研究報告［R］.成都:成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家專業(yè)實驗室，2005.

［15］許強，劉天翔，湯明高，等.三峽庫區(qū)塌岸預(yù)測新方法——岸坡結(jié)構(gòu)法［J］.水文地質(zhì)工程地質(zhì)，2007(3):110－115.

［16］重慶市公路局，重慶交通大學(xué)，重慶交通科研設(shè)計院，等.三峽庫區(qū)蓄水初期公路病害防治對策研究總報告［R］.重慶:重慶交通科研設(shè)計院，2005.

［17］宋岳，段世委，陳書文.官廳水庫塌岸影響因素分析［J］.水利水電工程設(shè)計，2004，23(1):34 －37.

［18］李永樂.三門峽水庫庫岸坍塌成因分析與防治措施研究［J］.水土保持學(xué)報，2003，17(6):129 －132.

［19］毛昶熙.滲流計算分析與控制［M］.北京:中國水利水電出版社，2003.

［20］高小育，廖紅建，丁春華.滲流對土質(zhì)邊坡穩(wěn)定性的影響［J］.巖土力學(xué)，2004，25(1):69 －72.

［21］左建，郭成久.水利工程地質(zhì)［M］.北京:中國水利水電出版社，2004:163.

［22］劉紅星，夏金梧，王小波，等.長江中下游岸坡變形破壞的主要型式及處理［J］.人民長江，2002，33(6):8 －10，27.

［23］崔中興，宋克強，仵彥卿，等.中小型水庫庫岸失穩(wěn)機理分析及其穩(wěn)定性評價方法［J］.西安理工大學(xué)學(xué)報，2000，16(1):37－41.

［24］包太，劉新榮，稅月.水位下降卸荷誘發(fā)庫岸邊坡快速失穩(wěn)機理分析［J］.水文地質(zhì)工程地質(zhì)，2004(5):7－11.

［25］吳俊杰，王成華，李廣信.非飽和土基質(zhì)吸力對邊坡穩(wěn)定的影響［J］.巖土力學(xué)，2004，25(5):732 －736，744.

［26］Hawkins R D.Sensitivity of sandstone strength and deformability to changes in moisture content［J］.Quarterly Journal of Engineering Geology，1992，25:115－130.

［27］張文杰，陳云敏，凌道盛.庫岸邊坡滲流及穩(wěn)定性分析［J］.水利學(xué)報，2005，36(12):1510 －1516.

［28］鄭穎人，時衛(wèi)民，孔位學(xué).庫水位下降時滲透力及地下水浸潤線的計算［J］.巖石力學(xué)與工程學(xué)報，2004，23(18):3203－3210.

［29］巫黎明，馮衛(wèi)兵.反弧型挑浪嘴式防洪墻對削減波浪爬高的試驗研究［J］.電力勘測設(shè)計，2004(1):39－43.

［30］周代榮，張丙先.工程地質(zhì)類比法在水庫塌岸預(yù)測中的應(yīng)用［J］.廣東水利水電，2007，4:13 －15.

［31］彭輝，劉德富，童廣勤.重慶市巫溪縣王家河庫岸139段塌岸預(yù)測研究［J］.水土保持研究，2006，13(4):107－110，113.

［32］曾顥.三峽庫區(qū)興山縣平邑口庫岸塌岸預(yù)測［J］.三峽大學(xué)學(xué)報:自然科學(xué)版，2008，30(1):27 －30.

［33］Nagata N，Hosoda T，Muramoto Y.Numerieal analiysis of river channel proeesses with bank erosion［J］.Journal of Hydraulic Engineering，2000，126(4):243－252.

［34］Elei S，Paul A.Reservoir Shore Line Erosion and Sediment Deposition with Cohesive Sediments［R］.Clemson，South Carolina:South Carolina Water Resources Center，Clemson University，2002.

［35］孫廣忠.水庫坍岸研究［M］.北京:水利電力出版社，1958.

［36］王征亮.三峽庫區(qū)長壽區(qū)庫岸塌岸預(yù)測的可拓學(xué)研究［D］.長春:吉林大學(xué)，2005.