堰塞體狀態(tài)相關(guān)剪脹理論與壩體潰決演化規(guī)律研究構(gòu)想

蔡正銀，鐘啟明，何 寧，夏云峰，朱俊高

(1.南京水利科學(xué)研究院 巖土工程研究所，江蘇 南京 210024；2.水利部水庫大壩安全重點實驗室，江蘇 南京 210029；3.南京水利科學(xué)研究院 河流海岸研究所，江蘇 南京 210024；4.河海大學(xué) 土木與交通學(xué)院，江蘇 南京 210098)

堰塞體是一種特殊的天然土石壩，一般由于地震和降雨引起的滑坡與泥石流阻塞河道，造成上游段壅水而形成[1-2]。堰塞體堵江形成堰塞湖作為重大水旱自然災(zāi)害，具有集雨面積廣、蓄水量大等特點，全球范圍內(nèi)有文獻記載的堰塞體案例有1 700余處，中國有800余處[3]。近年來，受地形地貌、地質(zhì)構(gòu)造及氣象水文等條件綜合作用，中國由于滑坡形成的堰塞體呈多發(fā)、頻發(fā)態(tài)勢。

堰塞體一般由崩滑土石料快速堆積而成，沒有經(jīng)過充分壓實，結(jié)構(gòu)較為松垮，組成物質(zhì)雜亂，局部存在由大顆粒骨架組成的高滲透區(qū)域，滲流和力學(xué)穩(wěn)定性較差[4]，且缺乏必要的洪水溢流設(shè)施，容易發(fā)生潰決造成嚴重的洪水災(zāi)害，對下游公眾生命財產(chǎn)和基礎(chǔ)設(shè)施構(gòu)成巨大威脅。2000年，易貢堰塞體潰決，21 億m3洪水下泄，導(dǎo)致中國墨脫、波密、林芝等縣90余鄉(xiāng)、近萬人受災(zāi)，印度布拉馬普特拉河沿岸7個邦94人死亡，250萬人無家可歸[5]；2008年，“5?12”汶川地震形成了257處滑坡堰塞體[6]，其中，唐家山是集雨面積最廣、蓄水量最大、威脅最嚴重的堰塞體，在人工干預(yù)下于2008年6月7日應(yīng)急泄流，共轉(zhuǎn)移下游風(fēng)險人口約25萬人，所幸未造成人員傷亡（圖1）[7]；2018年10—11月，中國金沙江和雅魯藏布江各接連發(fā)生2次滑坡事件，形成了白格[8]和加拉[9]堰塞體，并在短期內(nèi)發(fā)生潰決，對人民群眾生命財產(chǎn)安全構(gòu)成巨大威脅（圖2）。

圖1 唐家山堰塞體潰決Fig. 1 Breaching of Tangjiashan landslide dam

圖2 “11?3”白格堰塞體潰決Fig. 2 Breaching of “11?3” Baige landslide dam

Costa和Schuster[1]對全球73座堰塞體的壽命統(tǒng)計發(fā)現(xiàn)：85%的堰塞體壽命小于1 a，潰決模式主要是水流漫頂沖刷、滲透破壞或邊坡失穩(wěn)，其中，93%為水流漫頂沖刷潰決，4%為滲透破壞潰決，3%為邊坡失穩(wěn)潰決。Peng和Zhang[10]、Shen等[11]分別通過對全球204座和352座堰塞體的壽命統(tǒng)計也得出了類似結(jié)論。堰塞體的潰決雖然會帶來巨大的災(zāi)難，但有的堰塞體形成后結(jié)構(gòu)穩(wěn)定，使得堰塞湖得以保留，至今仍未潰決，如：位于塔吉克斯坦因1911年地震形成的Usoi堰塞體，其體積約為2 200 萬m3，高度約600 m，堰塞湖庫容約170 億m3，是目前世界上現(xiàn)存庫容最大的堰塞湖[12]；其他典型案例如中國重慶小南海堰塞體和瑞士Klontalersee堰塞體[1]。

中國歷來高度重視自然災(zāi)害防治，《中共中央關(guān)于制定國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和二〇三五年遠景目標的建議》對完善國家應(yīng)急管理體系，提高防災(zāi)、減災(zāi)、抗災(zāi)、救災(zāi)能力，防范和化解極端自然災(zāi)害事件系統(tǒng)性風(fēng)險提出更高要求。堰塞體致災(zāi)風(fēng)險大，破壞力強，嚴重威脅人們的生命財產(chǎn)和公共安全，而現(xiàn)有的理論和方法還不能很好地反映堰塞體的性態(tài)演化規(guī)律及潰決過程，亟需加強復(fù)雜條件下堰塞體的安全評價方法、災(zāi)變機理、災(zāi)害預(yù)測的基礎(chǔ)理論與關(guān)鍵技術(shù)研究，為應(yīng)急搶險和防災(zāi)減災(zāi)提供科學(xué)依據(jù)和技術(shù)支撐。

1 國內(nèi)外研究進展

1.1 土體剪脹理論與本構(gòu)模型

剪脹是指土體受剪引起的體積膨脹或收縮現(xiàn)象，是粒狀材料區(qū)別于普通固體材料的一個顯著特性。剪脹被定義為塑性體積應(yīng)變增量與塑性偏應(yīng)變增量之比。1962年，Rowe[13]提出經(jīng)典的應(yīng)力剪脹理論，其剪脹方程為：

式中，dεv為體積應(yīng)變增量，dε1為大主應(yīng)變增量，σ1為大主應(yīng)力，σ3為小主應(yīng)力，φcv為土體的極限摩擦角。

由式（1）看出，應(yīng)力剪脹理論中剪脹只與應(yīng)力水平有關(guān)，與材料的狀態(tài)（級配、密度等）無關(guān)[14]。對于正常固結(jié)黏土，由于應(yīng)力水平與土體密度之間存在一一對應(yīng)的關(guān)系，故經(jīng)典的應(yīng)力剪脹理論中實際上隱含了剪脹與密度的相關(guān)性。應(yīng)力剪脹理論描述黏性土的剪脹與應(yīng)力水平之間的關(guān)系，取得了巨大的成功，很多著名的黏性土本構(gòu)模型都是采用應(yīng)力剪脹理論，如劍橋模型及修正劍橋模型的剪脹方程，分別見式（2）和（3）：

對于無黏性土，相同應(yīng)力水平條件下，密度可能不同。正是由于經(jīng)典應(yīng)力剪脹理論對于材料內(nèi)部狀態(tài)參量的忽視，使該理論應(yīng)用到砂土或堆石料時出現(xiàn)較大誤差[14]。

砂土剪脹理論研究方面，20世紀80年代開始，國內(nèi)外學(xué)者將孔隙比、砂土的狀態(tài)參量或狀態(tài)指數(shù)引入剪脹方程，建立了一系列狀態(tài)相關(guān)剪脹方程，其中：最具代表性的是Li等[15]將同時反映砂土的應(yīng)力水平與砂土密度的狀態(tài)參量引入剪脹方程，提出砂土狀態(tài)相關(guān)剪脹理論；其后，蔡正銀等[16-18]開展了砂土剪切帶的形成過程、砂土的變形特性與臨界狀態(tài)和砂土的漸進破壞模擬等方面的研究工作，進一步完善了砂土的狀態(tài)相關(guān)剪脹理論。

堆石料剪脹理論研究方面，隨著中國一大批高土石壩工程開工建設(shè)，作為大壩主要填筑材料的堆石料的剪脹理論的合理性引起了學(xué)術(shù)界的廣泛關(guān)注。與砂土相比，堆石料粒徑大、易破碎，工程特性也不盡相同，將砂土的剪脹理論直接用于堆石料的剪脹變形特性存在一定的局限性。學(xué)者們針對不同初始密度的堆石料進行了一系列大型三軸試驗，研究了密度和圍壓對堆石料變形特性的影響[19]，探討了堆石料強度和變形特性[20]、堆石料臨界狀態(tài)[21]和剪脹方程[22]，但堆石料的相關(guān)剪脹理論還很不成熟。

本構(gòu)模型方面，目前應(yīng)力變形分析中應(yīng)用最廣泛的主要有“鄧肯-張”非線性彈性模型、“南水”彈塑性模型，以及其他一些彈塑性本構(gòu)模型。其中：“鄧肯-張”模型基于廣義虎克定律，不能反映土體的應(yīng)變軟化特性、剪應(yīng)力引起的體脹或體縮現(xiàn)象；“南水”模型通過引入體積屈服面和剪切屈服面，可較好地反映堆石料的卸荷體縮特性和各向異性。但上述模型僅從數(shù)學(xué)角度解釋了表觀現(xiàn)象，未考慮材料初始狀態(tài)（如級配、密度）的影響，復(fù)雜應(yīng)力路徑下應(yīng)力-應(yīng)變關(guān)系的適用性也未得到充分驗證。

與筑壩堆石料相比，堰塞體材料的級配更為寬泛，不同部位的密度和級配差異較大。隨著堰塞湖水位的抬升和其他外荷載的作用，以及材料自身的蠕變，材料密度的改變將貫穿堰塞體的整個生命周期，更增加了其變形特性研究的復(fù)雜性。目前，針對堰塞體材料是否存在臨界狀態(tài)的相關(guān)研究很少，在剪脹方程中需引入何種狀態(tài)參量才能合理反映材料的剪脹性尚不明晰；常用于土石壩數(shù)值分析的本構(gòu)模型都沒有考慮土的密度的影響，更沒有考慮級配的影響?；谝陨戏治?，弄清堰塞體材料變形及強度特性隨應(yīng)力、密度和級配等因素的變化規(guī)律，建立能同時考慮應(yīng)力水平與土體狀態(tài)（密度、級配）影響的堰塞體材料剪脹理論與本構(gòu)模型，并應(yīng)用于堰塞體的應(yīng)力-變形及穩(wěn)定性計算，具有重要的學(xué)術(shù)價值和現(xiàn)實意義。

1.2 土石材料沖蝕特性

堰塞體的潰決絕大多數(shù)由水流沖蝕導(dǎo)致，因此研究土石材料的沖蝕特性至關(guān)重要。土石材料的沖蝕過程實質(zhì)上是水流與土石顆粒之間的耦合作用過程，其沖蝕特性一般通過起動流速（或臨界剪應(yīng)力）和沖蝕速率表征。目前，常用的測定土石材料沖蝕的方法包括旋轉(zhuǎn)圓柱試驗、孔洞沖蝕試驗、噴射沖蝕試驗、水槽試驗等（圖3）。

圖3 常用土石材料沖蝕特性測試方法Fig. 3 Common test methods for erosion characteristics of soil and rock materials

旋轉(zhuǎn)圓柱試驗裝置是通過旋轉(zhuǎn)外層透明的圓柱體帶動水流產(chǎn)生作用于圓柱體試樣的剪應(yīng)力測量材料沖蝕特性。1962年，Moore和Masch[23]首次采用繞圓柱體土試驗的旋轉(zhuǎn)水流提供剪切力來測量表面沖蝕。隨后，學(xué)者們對旋轉(zhuǎn)圓柱測量設(shè)備進行了改進，與原來的設(shè)計相比，新設(shè)備能夠測定重塑土和原狀土[24]。旋轉(zhuǎn)圓柱測量裝置的局限性在于不能測算試樣在飽和狀態(tài)時的沖蝕率，且沖刷掉的土顆粒仍在圓柱體內(nèi)而影響測試效果。

孔洞沖蝕試驗裝置是通過在土體試樣內(nèi)部開鑿孔洞，測量水流沖蝕后孔洞內(nèi)徑的發(fā)展確定土體內(nèi)部沖蝕特征。2004年，Wan和Fell[25]首次采用孔洞試驗確定了土體各參數(shù)與沖蝕率的關(guān)系。但沖蝕率結(jié)果可能受試驗過程中孔洞的堵塞、內(nèi)徑變化不均勻等因素影響而相差很遠。

噴射沖蝕試驗裝置通過對土體試驗表面進行射流沖刷，確定土體的抗沖蝕能力。1991年，Hanson[26]提出材料抵抗力的概念，并建議了噴射指數(shù)；Chang等[27]應(yīng)用此方法現(xiàn)場測定了汶川地震造成的2座堰塞體的土體沖蝕特性。但由于堰塞體材料的寬級配特征和不均勻分布特性，測試結(jié)果的隨機性較大，對堰塞體材料沖蝕特性進行全面評價還存在困難。

綜上可知，旋轉(zhuǎn)圓柱試驗、孔洞和噴射沖蝕試驗一般用于粒徑均勻的原狀土或重塑土沖蝕特性的測量，對于寬級配的堰塞體材料，一般需采用水槽模型試驗研究其沖蝕特性。寬級配土體的起動特性較均勻土體復(fù)雜得多，主要是粗細顆粒間受到顯著的蔭暴和床沙粗化作用。1950年，Einstein[28]指出了采用均勻沙代替非均勻沙研究沖蝕特性可能導(dǎo)致的誤差。數(shù)十年來，國內(nèi)外學(xué)者圍繞非均勻顆粒起動特性和運動規(guī)律開展了廣泛研究。

對于非均勻顆粒起動問題：Paintal[29]提出暴露度的概念研究大小顆粒之間的相互作用；韓其為[30]基于床面顆粒相對位置，提出含義清晰且便于應(yīng)用的絕對暴露度和相對暴露度概念；其后，學(xué)者們又開展一系列相關(guān)研究，在一定程度上揭示了堰塞體顆粒的分布特征及其起動特性[31]。但近底水流流速、顆粒非均勻性及隨機分布規(guī)律對寬級配材料起動特性的影響仍值得深入研究。

對于非均勻顆粒的運動規(guī)律，根據(jù)顆粒在水流中的運動特征，一般將其分為推移質(zhì)和懸移質(zhì)進行研究。非均勻推移質(zhì)沖蝕率計算方法可分為3類：直接分組法[32]、修正剪應(yīng)力法[33]和級配法[34]。對于懸移質(zhì)的運動規(guī)律，潰壩水流的高含砂特性使其具有與普通挾沙水流不同的運動特性和沖淤規(guī)律，國內(nèi)外學(xué)者在懸移質(zhì)濃度分布[35]、水流挾砂力[36]等方面取得了系列的研究成果，但高速水流作用下的寬級配顆粒懸浮、擴散與交換機理尚需進一步研究。

總之，在土石材料沖蝕特性研究方面，國內(nèi)外已經(jīng)積累了許多寶貴的經(jīng)驗，但試驗條件與堰塞體實際沖蝕情況仍有較大差異，如：水流的強非恒定，土石料的寬級配、非均勻特征，并且缺乏與潰口變形直接相關(guān)的沖蝕特性模擬方法。

1.3 堰塞體潰決機理

近年來，國內(nèi)外學(xué)者基于不同尺度的水工模型試驗，開展了一系列堰塞體潰決機理和潰決過程模型試驗。小尺度模型試驗方面：Zhou等[37]基于4組壩高0.7 m的模型試驗，將堰塞體的潰決過程分為3個階段：溯源沖蝕→加速沖蝕→沖蝕衰退；Zhu等[38]基于12組壩高0.3 m的模型試驗，將堰塞體的潰決過程分為4個階段：初始沖蝕→溯源沖蝕→加速沖蝕→河道再平衡。大尺度模型試驗方面：Zhang等[39]基于3組壩高2.5 m的模型試驗，將堰塞體的潰決過程分為3個階段：溯源沖蝕→加速沖蝕→沖蝕衰退；Takayama等[40]基于2組壩高1.0 m的模型試驗，將堰塞體的潰決過程分為2個階段：漸進沖蝕→漫頂沖蝕?？傮w而言，盡管對堰塞體潰決過程有不同的劃分標準，國內(nèi)外的研究成果整體上體現(xiàn)了對潰決過程的一致認識。另外，國內(nèi)外學(xué)者圍繞堰塞體形態(tài)（坡比[41]、高度[42]）、材料物理力學(xué)特性（顆粒級配[43]、初始含水率[44]、密實度[45]）、河道水動力條件（河道坡比[46]、上游來流量[47]、涌浪[48]）、初始潰口形態(tài)[49]等影響堰塞體潰決過程的因素開展試驗，通過對潰口演化規(guī)律和潰口流量過程的分析，獲取了很多關(guān)于堰塞體潰決機理的有益認識。

隨著試驗技術(shù)的不斷提高、模型尺寸的不斷增大，測試精度有所提高，但小尺度模型試驗（也包括現(xiàn)場大尺度模型試驗）在堰塞體應(yīng)力水平、顆粒粒徑大小方面與原型相比仍然有較大差別，試驗中的縮尺效應(yīng)將直接影響研究成果的準確性，所得結(jié)論能否真實反映原型堰塞體的潰決機理仍值得商榷。值得一提的是，為有效解決小尺度潰壩試驗?zāi)Ｐ团c原型應(yīng)力水平相差過大，大尺度潰壩試驗場地難尋、耗時長、成本高、隨壩高增加風(fēng)險難以控制等問題，南京水利科學(xué)研究院利用離心機高速旋轉(zhuǎn)形成的超重力場具有的“時空放大”效應(yīng)的原理，成功研發(fā)一套基于NHRI-400gt離心機的土石壩潰壩離心模型試驗系統(tǒng)[50]，如圖4所示，創(chuàng)建了潰壩離心模型試驗方法，可在較短時間內(nèi)重現(xiàn)各類100 m級土石材料壩的潰決過程，為正確揭示堰塞體的潰決機理提供了科學(xué)先進的技術(shù)手段。

圖4 NHRI-400gt土石壩潰壩離心模型試驗系統(tǒng)Fig. 4 NHRI-400gt centrifugal model test system for earth-rock dam breach

1.4 堰塞體潰決過程數(shù)值模擬

數(shù)值模擬是潰壩過程預(yù)測的重要手段，潰壩數(shù)學(xué)模型一般分為3類[51]：參數(shù)模型、基于潰決機理的簡化數(shù)學(xué)模型、基于潰決機理的精細化數(shù)學(xué)模型。

參數(shù)模型大多針對潰壩案例進行數(shù)據(jù)統(tǒng)計，采用回歸算法或機器學(xué)習(xí)方法獲得潰壩相關(guān)參數(shù)；由于其公式簡單、計算快捷，也常用于潰壩致災(zāi)后果的快速評價。目前的參數(shù)模型多用來描述土石壩的潰決，專門針對堰塞體潰決的參數(shù)模型較少[3]。1985年，美國學(xué)者Costa[52]基于10座堰塞體潰決案例，最早提出3個預(yù)測潰口峰值流量的回歸方程。早期的參數(shù)模型僅簡單考慮了堰塞體的形態(tài)特征和堰塞湖水動力條件，且這些模型中只給出了潰口峰值流量的表達式，未提供其他潰決參數(shù)的表達式。近年來，Peng和Zhang[10]、石振明等[53]分別提出潰決參數(shù)快速預(yù)測模型，可綜合考慮堰塞體形態(tài)特征、堰塞湖水動力條件及堰塞體材料的沖蝕特性，模型可以預(yù)測潰口峰值流量、潰口最終尺寸（頂寬、底寬、深度）及潰壩歷時，取得了良好的模擬效果；但參數(shù)模型無法提供潰口演化過程和潰壩洪水流量過程線。

對于基于潰決機理的簡化數(shù)學(xué)模型，目前大多借鑒均質(zhì)壩模型，但由于堰塞體與均質(zhì)壩在材料特性和結(jié)構(gòu)特點上存在巨大差異，其潰決機理也有明顯不同，因此模擬結(jié)果的可靠性存疑。近年來，學(xué)者們也開發(fā)了一些可反映堰塞體潰決機理的簡化模型，一般選擇寬頂堰流公式模擬潰口流量，采用基于剪應(yīng)力的不同沖蝕公式模擬潰口的縱向下切和橫向擴展，建立不同的方法考慮堰塞體縱斷面的溯源沖蝕，并采用極限平衡法模擬潰口邊坡失穩(wěn)，比較有代表性的模型有香港科技大學(xué)的DABA[54]、中國水利水電科學(xué)研究院DB-IWHR[55]、南京水利科學(xué)研究院DB-NHRI[56]。這些模型已應(yīng)用于國內(nèi)易貢、唐家山、小崗劍、白格、加拉等堰塞體潰決案例的反演分析或潰決參數(shù)預(yù)測。總的來說，簡化數(shù)學(xué)模型的優(yōu)勢在于一定程度上考慮了堰塞體的潰決機理，且計算速度較快，在數(shù)值模擬中應(yīng)用最為廣泛；但模型一般假設(shè)潰口的形狀和演化過程，且無法真正模擬堰塞體潰決過程中的水土耦合作用。

近年來，隨著計算流體力學(xué)和泥沙科學(xué)的發(fā)展，產(chǎn)生了一系列以水動力學(xué)和平衡/非平衡輸沙理論為基礎(chǔ)，基于靜水壓力分布和淺水假設(shè)的1維、2維和3維潰壩過程數(shù)學(xué)模型，此類模型可稱為基于潰決機理的精細化數(shù)學(xué)模型。該類模型主要包括清水或渾水的連續(xù)性方程、動量守恒方程和潰口底床變化方程，有些包含潰口邊坡穩(wěn)定性判別方程，可精細地模擬土石材料壩的潰決過程。依據(jù)模型中選用沖蝕公式的類型，可將精細化模型分為4類[57]：平衡輸沙模型、非平衡輸沙模型、兩相流模型和兩層流模型。平衡輸沙模型[58]一般采用清水淺水方程描述水流運動，采用Exner方程描述潰口底床的變化過程，沖蝕過程以平衡輸沙為主，主要采用不同的推移質(zhì)（少量包含懸移質(zhì)）公式直接計算壩料的沖蝕率。非平衡輸沙模型[59]中假定過壩水流為渾水，采用含水土混合物密度的淺水方程描述潰壩水流運動；假定被沖蝕的壩料以推移質(zhì)和懸移質(zhì)形式運動，推移質(zhì)與懸移質(zhì)交互過程由經(jīng)驗公式確定，將實際沖蝕率折算成渾水的體積濃度，通過渾水流量與體積濃度計算沖蝕率。兩相流模型[60]假定固相顆粒在液相自由水流的驅(qū)動下發(fā)生移動，固相顆粒的體積濃度相對較低；固相和液相分別基于連續(xù)介質(zhì)假設(shè)建立模型，采用連續(xù)性方程和動量守恒方程描述其運動過程，其中，固相表現(xiàn)出類似顆粒流的特性。兩層流模型[61]將潰口底床上方的水流假定為兩層，分別為上方的清水層和中間的推移質(zhì)水流層，每層都有各自的深度和濃度，清水層與推移質(zhì)水流層之間存在清水交換，而推移質(zhì)水流層與潰口底床之間存在土體交換，交換量由經(jīng)驗公式確定，對清水層和含沙水層分別建立連續(xù)性方程與動量守恒方程得到兩層流模型。上述4類模型通常采用近似黎曼解法和全變差遞減法（TVD）等激波捕捉方法，采用有限體積法、光滑顆粒流體動力學(xué)法等數(shù)值模擬方法對控制方程進行求解。此類模型可實現(xiàn)對堰塞體潰決過程的精細化模擬，但目前多用于顆粒較為均勻的壩體材料。由于精細化模擬方法可考慮潰壩過程中復(fù)雜的水土交互作用，并可模擬復(fù)雜邊界條件和潰口演化規(guī)律，是潰壩數(shù)值模擬的發(fā)展方向。

1.5 堰塞體穩(wěn)定性評價方法

堰塞體穩(wěn)定性評價是開展相關(guān)風(fēng)險評估和應(yīng)急響應(yīng)的重要內(nèi)容，而堰塞體穩(wěn)定性又與其堆積形態(tài)、物質(zhì)組成和結(jié)構(gòu)特征、水動力條件及次生地質(zhì)災(zāi)害密切相關(guān)，按特性可分為定性和定量評價方法。

定性評價方法主要包括工程類比法和歷史分析法。前者是根據(jù)形成條件和地質(zhì)條件與同類堰塞體進行類比分析；后者是針對某一堰塞體的形成歷史和發(fā)育過程進行推測分析，以此評價堰塞體當(dāng)前的穩(wěn)定性情況。定性評價中，常借助InSAR監(jiān)測技術(shù)、無人機航拍、地面變形監(jiān)測等技術(shù)手段，主要以堰塞體的形成機制、物質(zhì)組成和結(jié)構(gòu)特征為基礎(chǔ)，判斷堰塞體的抗沖蝕性能，最終評價堰塞體的整體穩(wěn)定性[62]。定性評價的特點是不通過數(shù)學(xué)計算，而是利用收集的資料，通過類比分析得出結(jié)論，綜合考慮了多方面因素的影響，其分析結(jié)果可作為堰塞體穩(wěn)定性定量評價的基礎(chǔ)。

定量評價方法可概括為統(tǒng)計學(xué)法、物理模擬法和數(shù)值模擬法等3種方法。統(tǒng)計法基于收集獲取的大量已潰和未潰堰塞體的資料，從中選取評價指標，提出判別準則，并據(jù)此快速評價堰塞體的穩(wěn)定性。1999年，Casagli和Ermini[63]基于70座堰塞體案例，選取堰塞體體積和堰塞湖流域面積2個參數(shù)，最早提出堆積指標法（BI）評價堰塞體的穩(wěn)定性。其后，各國學(xué)者基于不同的案例提出一系列堰塞體穩(wěn)定性的快速評價方法，如基于多元回歸的方法[64]和基于邏輯回歸的方法[65]。這些評價方法主要圍繞堰塞體的形態(tài)參數(shù)和物質(zhì)組成，以及堰塞湖的水動力條件選取不同的評價指標。該類方法的優(yōu)點是避免復(fù)雜的計算，缺點是未考慮堰塞體內(nèi)部的應(yīng)力-應(yīng)變關(guān)系及其演化規(guī)律，可用于堰塞體穩(wěn)定性初步評價。物理模擬法是評價某一特定堰塞體最直觀的方法。在搜集堰塞體的地質(zhì)資料之后，利用相似的物理材料制成模型模擬原型堰塞體的變形失穩(wěn)過程，可通過施加不同的外荷載，模擬地震[66]、涌浪[67]等不同因素作用下堰塞體的穩(wěn)定性。該方法避免了現(xiàn)場試驗操作周期長且易受到外界多種因素干擾等缺點，也優(yōu)化了利用統(tǒng)計學(xué)方法得出的結(jié)果。但由于堰塞體內(nèi)部結(jié)構(gòu)復(fù)雜，物理模型不能完全模擬其真實狀態(tài)，再加上模擬試驗的尺寸效應(yīng)等影響，其結(jié)果往往是不全面的。數(shù)值模擬法能夠考慮堰塞體的力學(xué)特性和不同結(jié)構(gòu)面的影響，分析堰塞體的變形過程，對其穩(wěn)定性做出定量評價。其中，常用的應(yīng)力-應(yīng)變分析又分為連續(xù)介質(zhì)分析和不連續(xù)介質(zhì)分析。連續(xù)介質(zhì)分析包括有限單元法[68]和邊界單元法[69]兩大類；不連續(xù)介質(zhì)分析方法引入不連續(xù)分塊剛體模型，主要包括極限平衡法[70]、離散單元法[71]、不連續(xù)變形分析[72]。上述方法均為確定性方法，由于材料、結(jié)構(gòu)、荷載等參數(shù)存在不確定性，采用傳統(tǒng)的確定性壩體穩(wěn)定性評價方法的研究方向大多局限在形態(tài)、湖水位和物質(zhì)組成等方面，沒有充分考慮堰塞體材料的隨機性和狀態(tài)相關(guān)性、結(jié)構(gòu)空間變異性和荷載不確定性等問題。近年來，學(xué)者逐漸認識到巖土工程和水利水電工程中存在大量材料參數(shù)不確定和荷載不確定的問題，有學(xué)者將不確定性方法，如可靠度理論、模糊數(shù)學(xué)、灰色理論等數(shù)學(xué)方法應(yīng)用在結(jié)構(gòu)穩(wěn)定性評價中[73]，取得了較好的成效，但對于堰塞體的穩(wěn)定性采用不確定性方法的報道較少。

綜上所述，由于堰塞體復(fù)雜的外部形態(tài)、內(nèi)部材料和結(jié)構(gòu)特征，目前國內(nèi)外對于堰塞體狀態(tài)相關(guān)剪脹理論與壩體潰決演化規(guī)律的研究工作仍處于起步階段。有必要圍繞滑坡堰塞體安全評價方法與災(zāi)變機理開展深入研究，針對高邊坡、堰塞體及原始河道構(gòu)成的系統(tǒng)，充分考慮堰塞體的形態(tài)特征、物質(zhì)組成、結(jié)構(gòu)特點和堰塞湖水動力條件，揭示堰塞體多體、多過程、多尺度相互作用機理，闡明堰塞體性態(tài)演變過程與多過程作用的關(guān)聯(lián)機制，揭示堰塞體長期性態(tài)演變與潰決機理。提出能考慮流固耦合的堰塞體滲流、變形、穩(wěn)定和潰決過程大規(guī)模一體化數(shù)值模擬方法與穩(wěn)定性評價方法。為提升滑坡堰塞體應(yīng)急搶險決策水平和防災(zāi)減災(zāi)救災(zāi)能力，保障人民生命財產(chǎn)安全和國家社會經(jīng)濟可持續(xù)發(fā)展提供理論與技術(shù)支撐。

2 研究內(nèi)容

2.1 關(guān)鍵科學(xué)問題

本文提出主要探究的3個科學(xué)問題：1）寬級配堰塞體材料的剪脹規(guī)律；2）寬級配堰塞體材料的非恒定動態(tài)沖蝕特性；3）堰塞體漫頂潰決與滲透破壞機理。

對應(yīng)這3個科學(xué)問題的科學(xué)技術(shù)挑戰(zhàn)及重要性的內(nèi)涵分別為：1）由于滑坡體物質(zhì)來源和運動形式的差異，堰塞體呈現(xiàn)結(jié)構(gòu)性特征，物質(zhì)組成表現(xiàn)出寬級配特性，堰塞體材料表現(xiàn)出顯著的狀態(tài)相關(guān)性，常規(guī)剪脹理論形成的壩體應(yīng)力變形分析方法已無法滿足模擬的需求；亟需揭示應(yīng)力水平、相對密度、顆粒級配等因素影響下寬級配堰塞體材料的變形特性及臨界狀態(tài)變化規(guī)律，提出堰塞體材料的狀態(tài)相關(guān)剪脹方程。2）目前的沖蝕公式大多未考慮水流含砂濃度變化對其沖蝕特性的影響，無法合理反映堰塞體材料的寬級配特征，高速挾砂水流作用下潰口動邊界的演化規(guī)律仍需深入研究，有必要探明寬級配堰塞體材料在挾砂水流作用下的輸運規(guī)律，以及強非恒定流條件下的水沙耦合作用機制，揭示堰塞體材料在非恒定流作用下的動態(tài)沖蝕機理。3）無論何種致災(zāi)因子，堰塞體的最終潰決模式均表現(xiàn)為漫頂或滲透破壞，目前堰塞體潰決初期的滲流侵蝕、溯源沖刷和后期高流挾砂水流作用下沿程侵蝕的機理和轉(zhuǎn)化規(guī)律認識不清，堰塞體潰決不同階段的歷時與影響因素尚不明晰；亟需揭示堰塞體沖蝕模式的內(nèi)在機制，闡明不同潰決模式作用下，材料組成、結(jié)構(gòu)和形態(tài)及外在影響因子對堰塞體潰口演化機制與洪水流量過程的影響。

2.2 重點研究內(nèi)容

圍繞提升堰塞體安全評價與災(zāi)變模擬水平，需重點突破的方向展開3個層面的研究。機理層面：堰塞體外觀形態(tài)、內(nèi)部結(jié)構(gòu)和材料宏觀力學(xué)特性及其時空變異規(guī)律；堰塞體全生命周期性能演化規(guī)律；非恒定流作用下寬級配堰塞體材料的沖蝕特性；不同致災(zāi)因子作用下堰塞體的潰決機理。理論與方法層面：與堰塞體結(jié)構(gòu)和成分相適應(yīng)的物理力學(xué)性質(zhì)測試方法；堰塞體材料狀態(tài)相關(guān)剪脹理論；堰塞體材料彈塑性本構(gòu)模型與參數(shù)確定方法；不同荷載作用下堰塞體的穩(wěn)定性評價方法；貫穿堰塞體全生命周期的安全評價與災(zāi)變模擬方法。應(yīng)用層面：以唐家山、紅石巖、白格等已潰與未潰的典型堰塞體為示范工程和現(xiàn)場試驗基地，驗證并完善堰塞體穩(wěn)定性評價和災(zāi)變機理的理論體系和方法。研究內(nèi)容包括：

1）堰塞體材料狀態(tài)相關(guān)剪脹理論和穩(wěn)定分析方法：①堰塞體材料特性室內(nèi)試驗縮尺方法；②狀態(tài)相關(guān)的堰塞體材料強度變形特性與長期流變特性；③考慮狀態(tài)相關(guān)剪脹理論的堰塞體材料本構(gòu)模型與壩體穩(wěn)定分析方法。

2）挾砂水流作用下堰塞體材料的沖蝕機理：①堰塞體顆粒運動特性及試驗?zāi)M技術(shù)；②堰塞體材料非恒定流動態(tài)沖蝕特性；③堰塞體漫頂與滲透破壞潰口演化規(guī)律。

3）考慮流固耦合的堰塞壩潰決過程模擬方法：①寬級配堰塞體材料沖蝕方程；②堰塞體漫頂潰壩與滲透破壞過程數(shù)值分析方法；③堰塞體潰決過程3維重建。

4）堰塞體穩(wěn)定性評價方法：①堰塞體空間結(jié)構(gòu)識別及材料參數(shù)特性檢測方法；②不同外荷載作用下堰塞體穩(wěn)定性分析方法；③堰塞體穩(wěn)定性可靠度評價方法；④堰塞體滲流、變形、穩(wěn)定和潰決過程大規(guī)模一體化數(shù)值模擬方法。

3 研究方法與技術(shù)路線

圍繞構(gòu)建堰塞體全生命周期安全評價和災(zāi)變模擬的理論體系與方法這一核心問題，從基本理論與模擬方法，關(guān)鍵科學(xué)問題，機理揭示、理論創(chuàng)新與技術(shù)創(chuàng)新，方法集成，科學(xué)體系5個層次開展研究工作。采用多學(xué)科交叉和多手段融合的方法，基于傳統(tǒng)理論，借助各種物理和數(shù)值模擬方法，圍繞堰塞體安全評價、災(zāi)變機理與成災(zāi)過程中面臨的三大關(guān)鍵科學(xué)問題開展深入系統(tǒng)的研究。揭示堰塞體外觀形態(tài)、細觀結(jié)構(gòu)和材料特性時空變異規(guī)律，堰塞體全生命周期性能演化規(guī)律，寬級配堰塞體材料起動與沉積特征及堰塞體漸進破壞潰決機理；提出堰塞體材料狀態(tài)相關(guān)剪脹理論、堰塞體材料彈塑性本構(gòu)理論、堰塞體穩(wěn)定性可靠度評價理論、堰塞體材料非線性沖蝕模擬理論、堰塞體漸進潰決模擬理論等系列創(chuàng)新理論；實現(xiàn)堰塞體空間結(jié)構(gòu)識別與物料特性勘測技術(shù)、寬級配材料縮尺試驗技術(shù)、大型三軸試驗高精度體變測試技術(shù)、寬級配堰塞體潰決離心模型試驗技術(shù)等試驗和測試技術(shù)創(chuàng)新；通過模型試驗縮尺方法、流固耦合模擬方法、顯式與隱式數(shù)值算法、高性能軟件集成方法等各類方法的集成，形成堰塞體全生命周期安全評價與災(zāi)變模擬理論體系與方法，并在已潰與未潰的典型堰塞體案例中驗證和應(yīng)用。所采用的研究方法如圖5所示。

圖5 研究方法框架Fig. 5 Research methods and framework

秉承由物理機理到數(shù)學(xué)模型、從交叉到融合的研究方法，采用原位測試（外觀識別、空間結(jié)構(gòu)、物料組成、巖性特征）、單元試驗（組構(gòu)成分、物理力學(xué)特性）、大型物理模型試驗（水工模型、離心模型）、理論分析和數(shù)值模擬等多種研究手段，在堰塞體結(jié)構(gòu)與材料特性時空變異規(guī)律、復(fù)雜環(huán)境下堰塞體孕災(zāi)過程模擬、不同致災(zāi)因子作用下堰塞體潰決過程模擬、堰塞體安全評價與災(zāi)變過程模擬一體化系統(tǒng)等4個方面開展深入研究。同時，集成堰塞體特性識別技術(shù)、堰塞體工作性態(tài)演化模擬技術(shù)、堰塞體災(zāi)變過程模擬技術(shù)等貫穿堰塞體形成—孕災(zāi)—致災(zāi)過程的成套模擬方法，構(gòu)建以堰塞體全生命周期滲流、變形、穩(wěn)定和潰決過程一體化平臺為核心內(nèi)容的安全評價與災(zāi)變模擬系統(tǒng)，為中國堰塞體防災(zāi)減災(zāi)決策提供理論基礎(chǔ)與科學(xué)依據(jù)，采用的技術(shù)路線如圖6所示。

圖6 技術(shù)路線Fig. 6 Technology road

4 主要創(chuàng)新之處

課題組擬基于空天地一體化的勘察技術(shù)，并通過試驗技術(shù)的改進創(chuàng)新，采用理論研究和數(shù)值分析等方法，在如下4個方面取得突破：

1）建立合理反映堰塞體材料物理力學(xué)性質(zhì)、適應(yīng)復(fù)雜應(yīng)力路徑、統(tǒng)一模擬堰塞體應(yīng)力應(yīng)變特性的狀態(tài)相關(guān)剪脹理論與彈塑性本構(gòu)理論。

2）基于堰塞體彈塑性本構(gòu)模型與多相多場耦合理論，采用可靠度隨機有限元方法，建立考慮滲流、變形、穩(wěn)定的堰塞體性能演化規(guī)律的安全評價理論。

3）建立正確反映堰塞體寬級配材料沖蝕特性和不同致災(zāi)因子作用下的潰決機理，合理預(yù)測其潰口演化規(guī)律和洪水流量過程的堰塞體潰決過程模擬理論。

4）集成應(yīng)用GIS技術(shù)、數(shù)據(jù)庫管理技術(shù)、多維數(shù)值仿真模擬技術(shù)，開發(fā)高效的數(shù)值仿真平臺，實現(xiàn)考慮流固耦合的堰塞體滲流、變形、穩(wěn)定和潰決過程的大規(guī)模一體化數(shù)值模擬。

5 結(jié)論與展望

針對堰塞體的形態(tài)特征、物質(zhì)組成、結(jié)構(gòu)特點和堰塞湖水動力條件，擬通過空天地一體化的測量分析技術(shù)、大型三軸試驗、大型水槽模型試驗、離心模型試驗、理論分析和數(shù)值仿真等手段，研究寬級配堰塞體材料的剪脹規(guī)律，建立堰塞體材料的狀態(tài)相關(guān)本構(gòu)模型，基于可靠度計算方法，對堰塞體的穩(wěn)定性進行評價；揭示堰塞體材料非恒定流動態(tài)沖蝕特性和不同潰決模式下堰塞體的潰口演化規(guī)律，提出考慮流固耦合的堰塞體滲流、變形、穩(wěn)定和潰決過程大規(guī)模一體化數(shù)值模擬方法，為堰塞體的應(yīng)急處治提供理論與技術(shù)支撐。

擬解決的問題主要針對堰塞體形成后性態(tài)的時空演化規(guī)律，以及不同致災(zāi)因子作用下的潰決機理和潰決過程。目前，對于堰塞體的滑坡堵江過程及其誘發(fā)因素，以及堰塞體潰決后的泥石流演進致災(zāi)過程的研究仍面臨重大的科技挑戰(zhàn)。為進一步完善堰塞體的防災(zāi)減災(zāi)體系，未來應(yīng)深入開展“山體滑坡—堰塞體堵江—堰塞湖潰決—泥石流演進—災(zāi)害損失評估”災(zāi)害鏈發(fā)展演化過程的機理與模擬研究。