999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

堰塞湖風險處置研究現狀與展望

2024-05-14 09:14:40周招楊啟貴蔡耀軍盧建華陳朝旭張文傳
人民長江 2024年4期

周招 楊啟貴 蔡耀軍 盧建華 陳朝旭 張文傳

摘要:科學認知堰塞湖基本特征、高效評估堰塞湖風險以及提升堰塞湖應急處置技術是緩和堰塞湖威脅和降低堰塞湖災損的關鍵。從堰塞湖分布、成因、蓄水庫容、生命周期以及堰塞湖災害特點等方面系統分析了堰塞湖基本特征,從定性、定量以及災害鏈角度系統總結了堰塞湖危險性評估方法,并結合堰塞湖潰決機理系統分析當前堰塞湖潰決應急監測困境、引流槽除險技術水平和疏通開挖機械設備現狀?;诋斍把芯康木窒扌裕裏捬呷碧幹醚芯课磥響攸c關注的5個問題:① 堰塞體堆積形態分區機制;② 強非恒定流輸沙理論;③ 大尺度堰塞湖潰決模型試驗;④ 堰塞湖應急處置專用機械設備研制;⑤ 堰塞湖潰決洪峰動態追蹤。研究成果可為高風險堰塞湖風險評估、應急處置及區域防災減災規劃提供有益參考。

關鍵詞:堰塞湖; 潰決洪峰; 風險評估; 應急處置; 實時監測

中圖法分類號: TV143

文獻標志碼: A

DOI:10.16232/j.cnki.1001-4179.2024.04.002

0引 言

堰塞湖是高山峽谷地區因降雨、地震以及火山噴發等動力地質作用造成山體滑坡、坍塌或泥石流堵塞天然河道形成的天然湖泊,堵塞天然河道的固體堆積體則稱為堰塞體。早些年間堰塞湖科學研究較少,部分學者根據堵江成因將堰塞湖命名為地震湖、淤地壩以及滑坡堵江等。魯成公五年“梁山崩,雍河三日不流”便是《春秋》記載的中國早期堰塞湖事件[1-3]。近些年隨著認知水平提升,學者將堵塞河道形成的天然湖泊統一命名為堰塞湖[4-5]。

堰塞體堵塞天然河道,致使堰塞湖出入庫流量極不平衡,引起庫水位持續壅高。尤其是高山峽谷地區大流量河道,極易在短時間內造成堰塞體漫頂潰決,形成破壞威力強悍、災害鏈長遠的非常態潰決洪峰,嚴重威脅下游沿岸地區人民群眾生命財產安全和基礎設施安全,局部地區甚至對生態環境造成不可逆轉損毀[6-8]。1786年四川磨西7.7級地震引發形成摩崗嶺滑坡堰塞湖,漫頂潰決洪水造成沿江大量居民溺斃,死亡人數近10萬[9]。2000年西藏易貢堰塞湖潰決洪峰流量達到12.4萬m3/s,為雅魯藏布江年平均流量的28倍,嚴重損毀沿線公路、光纜以及橋梁等重要基礎設施[10-11]。

近些年受極端氣候影響,堰塞湖總體呈多發、頻發態勢,局部地區甚至反復發生,尤其是長江上游西藏、云南、四川等地區,堰塞湖潰決事件尤為頻繁[12-14]。2008年汶川地震引發多處崩塌、滑坡,形成257處堰塞湖,其中唐家山堰塞湖險情最為嚴重,直接威脅下游江油、綿陽等地130萬人民群眾的生命財產安全[15]。2018年長江上游發生加拉、白格等4次大規模堰塞湖潰決事件,其中“11·3”白格堰塞湖潰決洪峰在下游380 km范圍內均達到超萬年一遇洪水標準,造成沿岸迪慶、麗江及大理等地直接經濟損失高達74億元[16-17]。堰塞湖蓄水能力強、破壞威力大、波及范圍廣,且堰塞湖地處高山峽谷深處,惡劣的水域地理環境及交通通信條件和短暫的應急處置窗口期給堰塞湖防災減災造成極大困難,堰塞湖險情一直是極難控制和防御的自然災害[18-19]。

鑒于高風險堰塞湖潛在的巨大威脅,堰塞湖應急處置及風險防控是當前中國乃至全世界范圍內堰塞湖防災減災領域的重要研究課題。堰塞湖應急處置技術研究能為揭示堰塞湖漫頂潰決機理、人工調控堰塞湖潰決發展過程提供重要理論技術方法,對更大程度避免庫區淹沒損失和潰決洪水損毀具有重要的科學價值和社會意義[20]。因此,本文將結合當前堰塞湖風險處置研究現狀,對堰塞湖基本特征、風險評估和應急處置技術等研究現狀進行評述。

1堰塞湖基本特征

1.1堰塞湖分布

堰塞湖在世界范圍內廣泛分布在高山峽谷地區,尤其是板塊地殼運動劇烈、深切地貌、河流陡降等區域[21-23]。中國的堰塞湖主要集中在青藏高原及周邊秦巴山區地帶,這種特殊的堰塞湖分布主要與中國典型的三大階梯地勢緊密相關,局部地區小規?;禄蚰嗍鞅隳茉陔A梯結合地帶引起大范圍山體崩塌、滑坡,堵塞天然河道形成高風險堰塞湖。

1.2堰塞湖成因

堰塞湖成因主要是由坍塌誘發形成堰塞體的動力地質作用而定。Schuster等[21]對135座堰塞湖動力地質成因進行調查分析,指出降雨、融雪、地震及火山噴發等是誘發形成堰塞湖的主要因素,其中地震是堰塞湖形成的第二大誘發因素,約占總數的37%。石振明等[22]搜集國內外1 298座堰塞湖案例,如圖1所示,指出國外因地震、降雨誘發形成堰塞湖占比分別為23%和35%,相應國內因地震、降雨誘發形成堰塞

湖的比例分別高達53%及38%。國內外堰塞湖成因占比差異主要原因為中國是地震頻發且階梯地勢明顯的國家,單次地震或強降雨便能在階梯地勢結合帶誘發形成多處堰塞湖[23]??傮w而言,降雨及地震是堰塞湖形成的主要誘發因素。

1.3堰塞湖蓄水庫容

堰塞湖蓄水庫容隨機性較強,從幾百m3到上百億m3廣泛分布,其主要與堰塞體垂直高度、崩塌滑坡等地質災害規模、河溝河谷形態以及河道天然流量緊密相關,其中堰塞體作為堵塞天然河道的固體堆積物,其垂直高度直接影響堰塞湖蓄水庫容[24-25]。國內外堰塞湖蓄水庫容分布如圖2所示,堰塞湖蓄水庫容普遍集中在1×105~1×108 m3,部分高風險堰塞湖因堰塞體垂直高度較大及天然河道流量較高,蓄水庫容普遍達到大型水庫規模,甚至達到上百億m3。世界范圍內當前蓄水庫容規模最大的堰塞湖是位于塔吉克斯坦東部帕米爾高原上的薩雷茲堰塞湖,其堰塞體高度達567 m,最大水深505 m,蓄水庫容達170億m3。

1.4堰塞湖生命周期

堰塞湖生命周期是指坍塌滑源體沿滑槽滾落至河道、堵塞河道蓄水成湖直至漫頂潰決的全部過程。受堰塞體物質結構組成、幾何形態以及天然河道下泄流量等諸多因素影響,堰塞湖生命周期長短不一,跨度懸

殊。部分堰塞湖即生即消,即滑源體滾落至河道尚來不及形成堰塞湖便直接垮塌沖走;部分堰塞湖的堰塞體垂直高度大、堆積方量高,局部夾雜大孤石形成穩定支撐框架結構,造成堰塞湖滲漏流量長期出入庫平衡,致使堰塞湖屹立數百年甚至千年依舊穩定[25]。

堰塞湖生命周期直接影響堰塞湖應急處置窗口期,短生命周期堰塞湖將嚴重壓縮應急處置窗口期[26]。柴賀軍等[27]指出堰塞湖在1 d、10 d、6個月及1 a之內自然漫頂潰決的比例依次為20%,50%,80%及90%;石振明等[28]指出堰塞湖在1 h、1 d、1周、1個月、1 a之內自然潰決比例依次為9%,34%,50%,67%,86%。由此可見堰塞湖生命周期普遍較短,超過50%以上的堰塞湖將在1個月之內漫頂潰決,因此堰塞湖應急處置窗口期極為短暫,由于堰塞湖所處惡劣的水域環境及交通通訊障礙等自然環境因素,人工處置堰塞湖險情難度極大[29-31]。

1.5堰塞湖災害特點

堰塞湖庫水位持續壅高回水,在上游庫區造成大量淹沒損失,而一旦漫頂潰決,非常態潰決洪峰勢必肆虐下游沿岸基礎設施。堰塞湖危害主要體現在以下4個方面[3]:

(1) 上游庫區淹沒損失。堰塞體堵塞天然河道致使庫水位持續壅高、回溯,回水長度和淹沒面積持續增大,庫區沿岸農田、道路、工農業基礎設施及城鎮等被嚴重淹沒損毀(“11·3”白格堰塞湖上游庫區淹沒如圖3所示)。

(2) 潰決洪水損毀。堰塞湖非常態潰決洪峰以“刮地三尺、摧枯拉朽”態勢在數百千米范圍內嚴重肆虐沿岸堤防、鄉鎮、農田、生態環境及基礎設施。

(3) 潰決洪水對下游河道的影響。潰決水流流速大、攜沙濃度高、破壞威力強,極易沖毀沿岸河道堤防、造成決口,甚至改道,此外潰決水流裹挾大量泥沙,擁堵河道抬高河床,影響行洪能力。

(4) 殘存堰塞體威脅。殘存堰塞體兩側邊坡陡峭、松散,在強降雨或水流持續沖刷擾動下,極易坍塌、堵塞河道,再次形成堰塞湖。

2堰塞湖風險評估

早期堰塞湖風險評估無明確參考因子,堰塞湖風險等級難以劃分,堰塞湖應急決策難免失準,致使過高或過低評估堰塞湖風險。自2008年唐家山堰塞湖形成以來,國內外眾多學者結合堰塞體危險性和潰決洪水損毀,在堰塞湖風險評估理論和實踐方面均進行了積極探索。隨著不斷深入認識堰塞湖形成-潰決機理,堰塞湖災害系統和災害鏈認知不斷完善,堰塞湖風險評估開始向復雜環境下多致災因子的動態評估方向發展,并逐漸從定性、定量以及災害鏈等方面動態評估堰塞湖風險[32-35]。

2.1定性評估

堰塞湖風險定性評價主要結合堰塞體形成機制、物質組成及幾何形態,判斷堰塞體物質結構抗沖蝕性能,綜合評估堰塞體危險性。定性評估主要特點是依托衛星遙感、無人機航拍以及物探等檢測技術手段,結合同類堰塞湖歷史潰決案例,評估堰塞體危險性[35]。

胡卸文[36]、李勇[37]等利用堰塞體結構特征及堰塞湖蓄水庫容針對汶川地震所引發的堰塞湖進行初步風險分析。何秉順等[38]結合無人機航拍快速生成三維地形技術,短時間內對“5·12”汶川地震所引發的堰塞湖做了初步穩定分析。部分學者[39-40]根據宏觀監控現象及地表位移監測技術對堰塞湖進行風險評估。Graham[41]根據潰決洪水破壞威力、預警信息以及生命損失等因素評估堰塞湖風險。早期堰塞湖定性評估風險因子普遍考慮不周,隨著對堰塞湖認知不斷完善,2021年SL/T 450-2021《堰塞湖風險等級劃分與應急處置技術規范》補充考慮堰塞體結構形態、抗沖蝕特性以及湖區次生災害影響,尤其是生態環境影響,根據堰塞湖庫容、上游來水量、堰塞體物質組成及幾何形態、風險人口、影響城鎮、公共或基礎設施、生態環境等,將堰塞體危險性劃分為極高危險、高危險、中危險和低危險4個等級,有效提升了堰塞湖風險認知水平。

盡管堰塞湖定性評估結果并非完全精確,但在堰塞湖應急處置現場可初步分析堰塞湖風險?,F階段單純的堰塞湖定性風險評估相對甚少,定性評估普遍是定量評價的基礎,是迅速獲取堰塞湖地形信息及水文地理數據評估堰塞湖潰決風險和災損的基礎。

2.2定量評估

定量評估方法是運用大量堰塞湖歷史案例,分析漫頂潰決規律、擬合數學表達式或判別準則,評判堰塞湖風險。定量評估作為堰塞湖風險評估主流方法,根據計算方法特點可分為數理統計法、物理模型試驗法和數值分析法。

(1) 數理統計法。

數理統計法即是根據堰塞湖基礎資料,考慮堰塞湖水動力學指標、地貌學指標以及物質結構組成等要素,采用邏輯回歸方法,提出數學表達式或判別準則。1999年Casagli等[42]選用堰塞體體積和堰塞湖流域面積等2個參數提出堆積指標法(DBI)評估堰塞體穩定性;Ermini[43]、Dong[44]等根據堰塞體體積、流域面積、堰塞體高度、河道坡度、堰塞體物質結構組成以及回水長度等因子,提出堰塞湖風險評估模塊化方法。但數理統計法依賴堰塞湖歷史案例樣本以及風險因子合理性,不同樣本及不同風險因子堰塞湖風險定量評估稍有差異。

(2) 物理模型試驗法。

部分國內學者基于不同致災因子利用水槽物理模型試驗還原堰塞體潰決失穩過程,評價堰塞體穩定性[45-48]。盡管水槽物理模型能直觀地反映堰塞湖潰決發展過程、評估堰塞體穩定性,但當前堰塞湖水槽物理模型普遍是參照人工土石壩,堰塞體寬緩、矮胖等幾何特征以及寬級配物質組成等特性普遍被忽略[49]。蔡耀軍[20]、楊啟貴[50]等結合唐家山、白格等高風險堰塞湖指出堰塞湖堰塞體物質組成相差顯著,顆粒粒徑從毫米級砂石料跨越至數米乃至數十米孤石,同等高度堰塞體順水向長度遠大于垂直高度,呈寬緩斷面結構形態。而常規堰塞湖物理模型幾何垂直高度有限,漫頂水流極易引起堰塞體坍塌瞬潰,致使難以出現溯源陡坎現象,或因蓄水庫容不足造成難以持續沖刷堰塞體,致使殘留堰塞體過高,引發水力失真現象。蔡耀軍等[20]指出為盡可能避免縮尺效應影響,堰塞體頂部縱向長度需達到2~3倍垂直高度,垂直高度至少需要達到1 m以上規模。

(3) 數值分析法。

數值分析法是通過堰塞體物質結構組成、幾何形態及力學特性判斷分析堰塞體穩定性,分析堰塞體坍塌變形過程中應力應變關系。堰塞體穩定性數值分析法主要有極限分析法、有限差分法、有限單元法、離散單元法和非連續變形分析法等確定性分析方法[24],其中離散單元法和非連續變形分析法在堰塞體穩定性評估中應用較少。Mizuyama等[51]利用簡化畢肖普法考慮堰塞體形狀、物質結構組成及滲漏量等因素分析堰塞體穩定性;Hu等[52]利用瑞典條分法考慮不同水位及地震荷載等因素分析堰塞體穩定性。近些年學者嘗試利用模糊層次分析法、熵值法、模糊綜合分析法等不確定分析法開展堰塞體危險性評估研究,但數學模型理論與工程實際結合還存在較大差距[53-56]??傮w而言,堰塞體數值分析方法普遍選用人工土石壩模擬方法評估自身穩定性,但較少考慮堰塞體寬級配及寬緩幾何形態特點。

2.3災害鏈評價方法

持續壅高庫水位致使河道兩岸山體孔隙水壓力上升,極易誘發二次滑坡、形成大落差涌浪;漫頂潰決洪水嚴重肆虐沿岸堤防、生態環境及基礎設施,尤其是下游地區存在梯級堰塞湖,漫頂水流極易如“多米諾效骨牌”呈放大效應形成非常態潰決洪水,更進一步增大潰決洪水損毀。

當前堰塞湖風險評估風險因子有限,普遍局限于堰塞體自身穩定性及潰決洪水損毀,對堰塞湖災害系統及災害鏈普遍缺乏有效認知[57]。Niu等[58]首次研究梯級堰塞湖潰決機制,提出上游堰塞湖漫頂潰決時漫頂潰決洪峰直接影響下游堰塞湖潰決模式:低潰決洪峰逐漸沖蝕堰塞體,引起漸進式漫頂潰決,高潰決洪峰則快速沖刷掏蝕堰塞體,引發堰塞體坍塌瞬潰。Zhou等[59]研究了梯級堰塞湖動態潰決演化過程,指出梯形堰塞湖潰決洪峰呈疊加放大效應,潰決洪峰遠超單座堰塞湖最大潰決洪峰。Shi等[60]指出梯級堰塞湖潰決顯著改變河道匯水速率。當前梯級堰塞湖潰決研究尚處于初步階段,小尺度物理模型中潰決洪水放大效應難以體現,未來可考慮利用大中型水庫建模,更易于揭示潰決洪峰擴大效應。

3堰塞湖應急處置技術

針對堰塞湖險情,科研學者采用堰塞湖水文及地理信息應急監測和堰塞體快速疏通開挖等人工應急處置技術,快速識別堰塞湖水文、水情、庫區地形等相關信息,監測堰塞湖庫水位及潰決洪峰變化,評估堰塞湖風險,劃定堰塞湖上下游影響區域,轉移人民群眾。其中引流槽除險技術是當前應用最為廣泛、效果最為明顯的堰塞體快速疏通開挖應急處置技術。但堰塞湖應急處置技術普遍需根據堰塞湖基本特征及結合堰塞湖潰決發展過程特點,采取相應人工措施干預堰塞湖潰決發展過程,緩和堰塞湖險情[61-64]。

3.1堰塞湖潰決機理研究

非恒定攜沙水流動態沖刷掏蝕堰塞體,致使潰口呈間歇性非線性坍塌變形,潰決機理極為復雜,洞悉查明堰塞湖潰決機理是堰塞湖應急處置的基礎和先決條件。國內外學者在堰塞湖參數模型、物理模型及數值模擬等方面開展了大量堰塞湖潰決試驗研究,但因堰塞湖潰決過程涉及學科理論過多、影響因子復雜,當前潰決機理仍然不明,給堰塞湖應急處置造成巨大困境[65]。

部分學者[66-67]利用堰塞湖歷史案例推導潰口深度、潰口寬度及峰值流量等參數模型,但參數模型需要大量堰塞湖實測資料率定,有限的原型觀測案例降低了參數模型的可靠性。Zhong[48]、Begam[67]等采用數值模擬分析潰決過程參數敏感性以及潰決洪峰關鍵影響因子,但數值模型普遍假定潰口最終形態并利用寬頂堰公式計算潰口流量,而潰口形態呈非線性變化特征,與寬頂堰相距甚遠,且最終形態普遍難以假定。此外數值模擬潰口侵蝕率、擴展系數等參數受物質結構組成影響顯著,與實際潰決過程亦存在偏差。

周招等[68-69]開展室內物理模型試驗(圖4),將堰塞湖潰決發展過程依次劃分為尾部下切階段、陡坎溯源階段、全斷面下切階段以及沖淤平衡-恢復穩定階段。蔡耀軍等[20]認為堰塞體坍塌潰決主要集中在陡坎溯源階段及全斷面下切階段。楊興國[70]、張倩[71]、楊陽[72]等指出堰塞體漫頂潰決主要由潰決水流橫向展寬及縱向下切沖刷堰塞體引起坍塌失穩所致,但二者主導作用目前仍存在爭議:部分科研學者認為堰塞體坍塌主要是高速潰決水流剪切沖刷所致,部分學者則認為堰塞體坍塌是高速潰決水流掏蝕潰口底坡,致使潰口邊坡失穩所致,更有學者認為堰塞體坍塌是二者聯合相互作用所致[20]。

盡管物理模型試驗能有效模擬堰塞湖潰決發展過程,但當前幾乎所有物理模型普遍并未完全遵循重力相似準則選取寬級配顆粒材料堆筑堰塞體,部分細顆粒砂石料被統一概化,造成模型材料抗沖蝕性能與原型材料相差明顯。此外局部物理模型幾何尺度有限,垂直高度及蓄水庫容規模不足,致使漫頂水流積蓄勢能及持續沖刷能力不足,引起潰決發展過程失真等現象[73-75]。

3.2堰塞湖應急監測

堰塞湖通常地處交通險阻、人跡罕至之處,惡劣水域環境及水文、地質和物質組成等基礎資料的匱乏給人工監測堰塞湖險情提出了巨大挑戰,尤其是在有限應急處置窗口期內。

自2008年唐家山堰塞湖形成之后,學者廣泛開展堰塞湖技術研究。全球定位系統和非接觸式監測技術發展迅速,堰塞湖信息獲取逐步向快速感知、三維可視化以及潰口動態監測等方向發展[76-78]。

唐家山堰塞湖應急處置采用航空攝影測量和衛星遙感測量獲得堰塞湖周邊地理和地形地質資料;白格堰塞湖搶險中采取星載InSAR監測、無人機航拍等技術手段對堰塞體周邊不穩定體進行跟蹤監測。盡管星載InSAR技術具有全天候、全天時、覆蓋范圍廣等優點,但其觀測結果為地表真實形變在雷達視線上的投影,難以反映地物的真實變形特征及其演化規律[78-80]。

堰塞體三維可視化普遍利用高密度電法物理探測設備感知獲取堰塞體物質結構信息,并利用搭載集雙目相機、可見光相機、激光測距機、紅外熱像儀于一體的五目協同無人機快速查看堰塞體幾何形態,劃分堰塞體表觀顆粒粒徑并監測堰塞體邊坡變形[81-82]。陸二男等[83]在唐家山堰塞湖應急處置過程中采用高密度電法對引流槽進行了物探測試解譯;吳學明等[84]在紅石巖堰塞湖應急處置中采用被動源面波法、瞬變電磁法對堰塞體和古滑坡區開展探測。但堰塞湖三維可視化技術難以準確獲取堰塞體散粒料堆積分區,其對堰塞體內部物質結構探測深度及精度有限,且堰塞體變形及物質結構探測解譯速度有待提高。

針對堰塞湖潰口快速動態演變及現場惡劣的環境,至今全球范圍內均無人工或自動化野外監測堰塞湖潰口形態變化的實例?,F階段僅有的免棱鏡全站儀和手持激光測距儀等非接觸式監測設備僅適用于常規堤防潰口;常規聲吶探測技術因潰決水流攜沙濃度較高及沖刷掏蝕能力較強,難以固定于堰塞湖引流槽內監測潰口變形[85]。現階段主要是利用無人機航拍監測潰口展寬,但潰口縱向下切監測目前尚處于空白領域。近些年蔡耀軍等[20]提出堰塞湖引流槽沖刷變形自動化動態監測設備,如圖5所示,一定程度上可以動態監測引流槽橫向展寬及縱向下切變形,但動態監測設備依然面臨被潰決水流沖走而漂浮或懸浮于潰決水流中的技術挑戰。

3.3堰塞湖引流槽除險技術

堰塞湖引流槽除險技術基本原理是在堰塞體頂部縱向開挖引流槽,引導漫頂水流漸進式沖刷掏蝕堰塞體,平穩下泄堰塞湖囤蓄洪水。為提升引流槽泄流效率,周招等[86]在傳統梯形斷面引流槽基礎上優化引流槽橫、縱斷面結構型式,提出復式及陡坎式引流槽,并利用野外模型試驗證實復式及陡坎式引流槽能有效提升引流槽泄流效率,降低堰塞湖最大壅高水位及潰決洪峰流量。

盡管引流槽除險技術能一定程度降低堰塞體過水高程、削減潰決洪峰,但部分高風險堰塞湖蓄水庫容規模龐大,堰塞體物質結構疏松,持續壅高漫頂水流極易引起引流槽迅速變形,致使堰塞體坍塌失控,甚至瞬潰,形成難以控制的潰決洪峰[88]。

針對當前堰塞湖潰決發展過程難以控制的工程難題,現階段的鋼筋防護網、巨石鎖口和人工結構體等引流槽防護措施[87]能在尾部下切階段減緩潰決水流沖刷掏蝕引流槽,但在全斷面下切階段普遍因自身柔韌性不足而難以隨引流槽坍塌變形,調整自身結構型式,護坡固壁持續能力不足,調控堰塞湖潰決發展過程作用有限,甚至攔阻引流槽正常泄流,壅高庫水位,激起反噬作用。周招[86]、蔡耀軍[20]等提出引流槽控潰削峰設想,即提升尾部下切階段引流槽泄流效率,減緩全斷面下切階段堰塞體坍塌,并利用室外模型試驗證實掛壁式引流槽柔性防護措施能隨引流槽邊坡坍塌持續調整自身結構型式,減緩潰決水流橫向展寬及縱向下切掏蝕引流槽,延緩堰塞體坍塌,實現降低堰塞湖最大壅高水位及削減潰決洪峰(圖6)。但引流槽柔性防護措施布置形式極易影響削峰效果,鋪設不當甚至會攔阻引流槽泄流,激起反噬作用。

引流槽除險技術能有效降低堰塞湖最大壅高水位,削減潰決洪峰,但受堰塞湖潰決機理不明等諸多因素限制,當前堰塞湖引流槽除險技術更多依靠科研學者工程經驗開展。

3.4堰塞體快速疏通用開挖機械設備

近些年眾多堰塞湖險情得到成功處置,但當前并無堰塞體快速疏通開挖專用機械設備及應急處置專業隊伍。現階段堰塞湖應急處置主要是依靠中國安能建設集團有限公司(前身武警水電部隊)及常規基建施工隊伍駕駛挖掘機和自卸式重型卡車等機械設備快速進場,開挖輸運堰塞體。但堰塞湖水陸交通環境惡劣,常規大型機械設備因尺寸及重量等原因難以快速進場,如2008年唐家山堰塞湖應急部門不得已臨時開辟空中通道、租用米26直升機(最大載重20 t、主旋翼直徑32 m)吊運重型挖掘機進場施工作業。

武警消防部門常用橡皮艇(載重1 t,最多可乘載7人)臨時搭載科研設備及人員登頂堰塞體,實地調查堰塞湖第一手資料[88]。中國安能建設集團有限公司研制動力浮橋(容許通行荷載80 t)、HZQL26應急機械化橋(容許通行荷載60 t、跨度26.6 m)和臨時硬質路面機動鋪設系統(容許通行荷載60 t),開辟堰塞湖水上通道,有效解決軟質地面常規重型機械設備難以快速進場難題[89-91]。北京中瀚環球真空流體科技有限責任公司利用真空虹吸原理研發大管徑、大流量虹吸泄流裝備(最大流量5 m3/s),可有效解決大流量虹吸裝備脫流問題。但其跨越寬度有限,僅有8 m,僅適用于中小型堰塞湖應急排水,面對高風險堰塞湖(庫容普遍達到億m3量級)作用甚微,只能緩和壅水速度[92]。中國葛洲壩集團勘測設計有限公司研制了裝配式模塊化土方連續開挖輸運機械設備(模塊單元最大重12 t),如圖7所示,滿足常規直升飛機吊運要求,取得堰塞體開挖轉運一體化突破,顯著提升了堰塞體疏通開挖施工效率[2]。

上述堰塞湖應急處置機械設備在部分堰塞湖應急處置現場取得了良好應用效果,但總體而言當前堰塞湖應急處置針對性機械裝備仍然偏少,機械設備投送能力仍然偏弱,以至于在窗口期短暫的應急處置現場顯得力不從心,裝配式、輕量化工程機械設備(吊車、推土機、挖掘機、鑿碎機等)仍有待進一步開發研究。

4結論與展望

堰塞湖作為高山峽谷地區一種常見自然災害,極易在短時間內漫頂潰決,形成非常態潰決洪水,嚴重威脅沿岸人民群眾生命財產及重要設施安全。本文結合國內外堰塞湖案例及唐家山和白格等大型高風險堰塞湖實踐經驗,系統總結堰塞湖成因、蓄水庫容、生命周期等基本特征,分析了堰塞湖危險性定性、定量評估方法,闡釋了堰塞湖潰決機理、引流槽除險和疏通開挖機械設備現狀,可為后續堰塞湖風險處置提供技術指導。受堰塞湖潰決復雜水沙耦合過程影響,盡管現階段學者取得了一系列的成果,潰決機理認知水平有所提升,但堰塞湖風險處置技術研究仍存在如下科學難題:

(1) 堰塞體堆積形態分區機制。堰塞體堆積形態及物質結構分區是堰塞湖潰決過程模擬和預報的先決基本參數,但堰塞湖應急處置窗口期短暫、水域環境惡劣、堰塞體幾何規模龐大,常規物探設備難以快速獲取堰塞體物質結構組成及分區。堰塞體砂石料堆積形態分區機制是滑源體失穩-啟動-運動-停歇-堆積過程的能量耗散和轉換機制,是堰塞體物質結構組成和滑源體堆積過程的內在聯系,現階段亟需建立快速判別堰塞體物質結構組成的方法。

(2) 強非恒定流輸沙理論?,F有泥沙輸移理論普遍是在恒定均勻流條件下建立的砂石料沖刷搬運規律。堰塞湖潰決水流及潰口形態均隨潰決發展過程呈強非恒定流、非均勻流特征,泥沙輸移理論普遍未考慮與潰決形態變形相關的砂石料起動流速、推移質運動和懸移質運動規律。堰塞體松散堆積體無輸移搬運歷史觀測資料,現階段的泥沙輸移理論外延應用難免存在誤差。

(3) 大尺度堰塞湖潰決模型試驗。當前堰塞湖潰決數值模擬人為假定因素過多,難以真實反映潰口展寬及下切非線性變形過程,而堰塞湖潰決物理模型尺度有限,蓄水庫容及潰決水流流速不足,堰塞體幾何形態與砂石料顆粒級配相似性(幾何相似、重力相似等)難以兼容。因此極有必要開展符合模型相似律的大尺度物理模型試驗,分析研究潰決水流沖刷搬運機制及潰口展寬下切變化過程。

(4) 堰塞湖應急處置專用機械設備研制。堰塞湖事發突然、水域環境惡劣、交通險阻,應急處置窗口期短暫,常規堰塞體疏通開挖等重型機械設備因尺寸或重量等原因難以快速進場。因此極有必要開展堰塞湖應急處置保通進場及疏通開挖等裝配式、輕量化機械設備研制,便于以水陸空等方式快速進場,及時疏通開挖堰塞體,有效提升作業效率。

(5) 堰塞湖潰決洪峰動態追蹤。堰塞湖潰決洪水瞬息萬變,破壞能力強,常規水文站極易因水位變化跨度大、頻率快而難以實時監測潰決洪峰,甚至出現水文站淹沒損毀的情況。因此極有必要研發適應能力強、響應時間短的非接觸式潰決洪峰動態追蹤技術,以便快速評估堰塞湖傳播影響范圍。

參考文獻:

[1]劉寧,楊啟貴,陳祖煜.堰塞湖風險處置[M].武漢:長江出版社,2016.

[2]蔡耀軍,李建清,彭文祥,等.堰塞湖應急處置[M].武漢:長江出版社,2021.

[3]劉寧,陳尊蘭,崔鵬.堰塞湖及其風險處置[M].北京:科學出版社,2013.

[4]蔡耀軍,金新平,楊啟貴,等.金沙江白格堰塞湖風險統合管理與應急處置[M].武漢:長江出版社,2020.

[5]蔡耀軍,欒約生,彭文祥,等.堰塞湖典型案例分析成因、致災、應急處置[M].武漢:長江出版社,2023.

[6]JOHANNES T W.Case history and hazard analysis of two lake-damming landslides in the Himalayas[J].Journal of Asian Earth Sciences,1998,16(2):323-331.

[7]LI M H,HSU M H,HSIEH L S,et al.Inundation potentials analysis for Tsao-Ling landslide lake formed by Chi-Chi earthquake in Taiwan[J].Natural Hazards,2002,25:289-303.

[8]SHANG Y Y,YANG Z F,LI L H,et al.A super-large landslide in Tibet in 2000:background,occurrence,disaster,and origin[J].Geomorphology,2003,54:225-243.

[9]周洪福,韋玉婷,王運生,等.1786年磨西地震觸發的摩崗嶺滑坡演化過程與成因機理[J].成都理工大學學報(自然科學版),2017,44(6):649-658.

[10]邢愛國,徐娜娜,宋新遠.易貢滑坡堰塞湖潰壩洪水分析[J].工程地質學報,2010,18(1):78-83.

[11]DELANEY K B,EVANS S G.The 2000 Yigong landslide (Tibetan Plateau),rockslide dammed lake and outburst flood:review,remote sensing analysis,and process modelling[J].Geomorphology,2015,246:377-393.

[12]LIU N,ZHANG J X,LIN W,et al.Draining Tangjiashan Barrier Lake after Wenchuan Earthquake and the flood propagation after dam break[J].Science China Series E:Technological Sciences,2009,52(4):801-809

[13]LI J,CAO Z X,CUI Y F,et al.Hydro-sediment-morphodynamic processes of the baige landslide-induced barrier Lake,Jinsha River,China[J].Journal of Hydrology,2021,596:126134.

[14]ZHANG Z,HE S,LIU W,et al.Source characteristics and dynamics of the October 2018 Baige landslide revealed by broadband seismograms[J].Landslides,2019,16(4):777-785.

[15]崔鵬,韓用順,陳曉清.汶川地震堰塞湖分布規律與風險評估[J].四川大學學報(工程科學版),2009,41(3):35-42.

[16]ZHANG L M,XIAO T,HE J,et al.Erosion-based analysis of breaching of Baige landslide dams on the Jinsha River,China,in 2018[J].Landslides,2019,16(10):1965-1979.

[17]CHEN C,ZHANG L,XIAO T,et al.Barrier lake bursting and flood routing in the Yarlung Tsangpo Grand Canyon in October 2018[J].Journal of Hydrology,2020,583:124603.

[18]蔡耀軍,周招,楊興國,等.堰塞湖風險評估快速檢測與應急搶險技術和裝備研發[J].巖土工程學報,2022,44(7):1266-1280.

[19]周興波,周建平,杜效鵠,等.白格堰塞湖搶險處置應急管理與經驗啟示[J].水利學報,2021,52(10):1229-1239.

[20]蔡耀軍,楊興國,周招,等.基于物理模擬的堰塞湖溢流潰決機理[J].工程科學與技術,2023,55(1):150-160.

[21]SCHUSTER R L,COSTA J E.The formation and failure of natural dams[J].Geological Scoeity of America Bulletin,1988,100(7):1054-1068.

[22]石振明,馬小龍,彭銘,等.基于大型數據庫的堰塞壩特征統計分析與潰決參數快速評估模型[J].巖石力學與工程學報,2014,33(9):1780-1790.

[23]孔紀名,阿發友,鄧宏艷,等.基于滑坡成因的汶川地震堰塞湖分類及典型實例分析[J].四川大學學報(工程科學版),2010,42(5):44-51.

[24]年廷凱,吳昊,陳光齊,等.堰塞壩穩定性評價方法及災害鏈效應研究進展[J].巖石力學與工程學報,2018,37(8):1796-1812.

[25]單熠博,陳生水,鐘啟明.堰塞體穩定性快速評價方法研究[J].巖石力學與工程學報,2020,39(9):1847-1859.

[26]劉寧.堰塞湖應急處置實踐與認識[J].水科學進展,2010,21(4):541-549.

[27]柴賀軍,劉漢超,張倬元.中國滑坡堵江的類型及其特點[J].成都理工學院學報,1998(3):60-65.

[28]石振明,周明俊,彭銘,等.崩滑型堰塞壩漫頂潰決機制及潰壩洪水研究進展[J].巖石力學與工程學報,2021,40(11):2173-2188.

[29]ALICE T,AGNES H,JEAN-PHILIPPE M,et al.Seismic monitoring of soft-rock landslides:the Super-Sauze and Valoria case studies[J].Geophysical Journal International,2013,3:1-22.

[30]LIU Z P,GUO X L,ZHOU X B,et al.Cascading dam breach process simulation using a coupled modeling platform[J].Science China(Technological Sciences),2019,62(8):1455-1466.

[31]STUMPF A,MALET J P,DELACOURT C.Correlation of satellite image time-series for the detection and monitoring of slow-moving landslides[J].Remote Sensing of Environment,2017,189:40-55.

[32]崔鵬,郭劍.溝谷災害鏈演化模式與風險防控對策[J].工程科學與技術,2021,53(3):5-18.

[33]杜鎮瀚,鐘啟明,董海洲,等.堰塞湖風險評估研究綜述[J].河海大學學報(自然科學版),2022,50(4):13-25.

[34]石振明,熊永峰,彭銘,等.堰塞湖潰壩快速定量風險評估方法:以2014年魯甸地震形成的紅石巖堰塞湖為例[J].水利學報,2016,47(6):742-751.

[35]鐘啟明,陳小康,梅勝堯,等.滑坡堰塞湖潰決風險與過程研究進展[J].水科學進展,2022,33(4):659-670.

[36]胡卸文,韓玫,梁敬軒,等.汶川地震災區泥石流若干關鍵問題[J].西南交通大學學報,2016,51(2):331-340.

[37]李勇,黃潤秋,DENSMORE A,等.汶川8.0級地震的基本特征及其研究進展[J].四川大學學報(工程科學版),2009,41(3):7-25.

[38]何秉順,丁留謙,王玉杰,等.四川安縣肖家橋堰塞湖穩定性初步評估[J].巖石力學與工程學報,2009,28(增2):3626-3631.

[39]孫黎明.高山峽谷區滑坡堰塞體快速感知與模擬計算方法研究:以白格堰塞湖為例[J].水利水電技術(中英文),2021,52(7):44-52.

[40]WANG W,YIN K,CHEN G,et al.Practical application of the coupled DDA-SPH method in dynamic modeling for the formation of landslide dam[J].Landslides,2019,16(5):1021-1032.

[41]GRAHAM W J.A procedure for estimating loss of life caused by dam failure[J].Sedimentation & River Hydraulics,1999,6(5):1-43.

[42]CASAGLI N,ERMINI L.Geomorphic analysis of landslide dams in the Northern Apennine[J].Chikei,1999,20(3):219-249.

[43]ERMINI L,CASAGLI N.Prediction of the behaviour of landslide dams using a geomorphological dimensionless index[J].Earth Surface Processes & Landforms,2003,28(1):31-47.

[44]DONG J J,LI Y S,KUO C Y,et al.The formation and breach of a short-lived landslide dam at Hsiaolin village,Taiwan—Part I:post-event reconstruction of dam geometry[J].Engineering Geology,2011,123(1/2):40-59.

[45]OUYANG C Y,AN H C,ZHOU S,et al.Insights from the failure and dynamic characteristics of two sequential landslides at Baige village along the Jinsha River,China[J].Landslides,2019,16(3):1397-1414.

[46]PENG M,ZHANG L M.Analysis of human risks due to dam break floods-part 1:a new model based on Bayesian networks[J].Natural Hazards,2012,64(1):903-933.

[47]DANKA J,ZHANG L M.Dike failure mechanisms and breaching parameters[J].Journal of Geotechnical and Geoenvironmental Engineering,2015,141(9):04015039.

[48]ZHONG Q M,CHEN S S,WANG L,et al.Back analysis of breaching process of Baige landslide dam[J].Landslides,2020,17:1681-1692.

[49]PENG M,ZHANG L M.Breaching parameters of landslide dams[J].Landslides,2012,9(1):13-31.

[50]楊啟貴,楊文俊.5·12唐家山堰塞湖應急處置[M].武漢:長江出版社,2020.

[51]MIZUYAMA T,SATOHUKA Y,OGAWA K,et al. Estimating the outflow discharge rate from landslide dam outbursts[M]∥Disaster Mitigation of Debris Flows Slope Failures and Landslides.Tokyo:UniversalAcademic Press,2006(1):365-377.

[52]HU X W,LUO G,LU X P,et al.Analysis on dam-breaking mode of Tangjiashan barrier dam in Beichuan County [J].Journal of Mountain Science,2011,8(2):354-362.

[53]BAYER B,SIMONI A,MULAS M,et al.Deformation responses of slow moving landslides to seasonal rainfall in the Northern Apennines,measured by InSAR[J].Geomorphology,2018,308:293-306.

[54]FAN X,XU Q,ALONSO-RODRIGUEZ A,et al.Successive landsliding and damming of the Jinsha River in eastern Tibet,China:prime investigation,early warning,and emergency response[J].Landslides,2019,16(5):1003-1020.

[55]PENG M,ZHANG L M.Dynamic decision making for dam-break emergency management-part 1:theoretical framework[J].Natural Hazards and Earth System Science,2013,13:425-437.

[56]SHI Z M,XIONG X,PENG M,et al.Risk assessment and mitigation for the Hongshiyan landslide dam triggered by the 2014 Ludian earthquake in Yunnan,China[J].Landslides,2017,14(1):269-285.

[57]程尊蘭,崔鵬,李泳,等.滑坡、泥石流堰塞湖災害主要的成災特點與減災對策[J].山地學報,2008(6):733-738.

[58]NIU Z P,XU W L,LI N W,et al.Experimental investigation of the failure of cascade landslide dams[J].Journal of Hydrodynamics,Ser B,2012,24(3):430-441.

[59]ZHOU G D,CUI P,CHEN H Y,et al.Experimental study on cascading landslide dam failures by upstream flows[J].Landslides,2013,10(5):633-643.

[60]SHI Z M,GUAN S G,PENG M,et al.Cascading breaching of the Tangjiashan landslide dam and two smaller downstream landslide dams[J].Engineering Geology,2015,193:445-458.

[61]SHAN T,ZHAO J.A coupled CFD-DEM analysis of granular flow impacting on a water reservoir[J].Acta Mechanica,2014,225(8):2449-2470.

[62]LIU W,HE S.Dynamic simulation of a mountain disaster chain:landslides,barrier lakes,and outburst floods[J].Natural Hazards,2018,90(2):757-775.

[63]CARRIVICK J L.Dam break-outburst flood propagation and transient hydraulics:a geosciences perspective[J].Journal of Hydrology,2010,380(3/4):338-355.

[64]趙天龍,陳生水,王俊杰,等.堰塞壩漫頂潰壩離心模型試驗研究[J].巖土工程學報,2016,38(11):1965-1972.

[65]ZHAO T,DAI F,XU N W.Coupled DEM-CFD investigation on the formation of landslide dams in narrow rivers [J].Landslides,2017,14(1):189-201.

[66]CAO Z X,YUE Z Y,PENDER G.Flood hydraulics due to cascade landslide dam failure[J].Journal of Flood Risk Manage,2011,4(2):104-114.

[67]BEGAM S,SEN D,DEY S.Moraine dam breach and glacial lake outburst flood generation by physical and numerical models[J].Journal of Hydrology,2018,563:694-710.

[68]周招,蔡耀軍,彭文祥,等.高危堰塞湖引流槽結構形式優化試驗研究[J].人民長江,2023,54(2):200-205,219.

[69]周招,蔡耀軍,李建清,等.堰塞湖潰決過程人工干預試驗研究[J].水電能源科學,2022,40(10):109-114.

[70]楊興國,李海波,廖海梅,等.滑坡-堰塞湖應急處理與引流泄水優化技術[J].工程科學與技術,2018,50(3):95-104.

[71]張婧,曹叔尤,楊奉廣.堰塞壩局部潰決輸沙率試驗研究[J].泥沙研究,2011(6):32-36.

[72]楊陽,曹叔尤.堰塞壩漫頂潰決與演變水槽試驗指標初探[J].四川大學學報(工程科學版),2015,47(2):1-7.

[73]王光謙,王永強,劉磊,等.堰塞壩及其潰決模擬研究評述[J].人民黃河,2015,37(9):1-7.

[74]HU P,CAO Z X,PENDER G,et al.Numerical modelling of riverbed grain size stratigraphic evolution[J].International Journal of Sediment Research,2014,29:329-343.

[75]XUY,ZHANG L M.Breaching parameters for earth and rockfill dams[J].Journal of Geotechnical and Geoenvironmental Engineering,2009,135(12):1957-1970.

[76]黃艷,馬強,吳家陽,等.堰塞湖信息獲取與潰壩洪水預測[J].人民長江,2019,50(4):12-19,52.

[77]熊瑩,周波,鄧山.堰塞湖水文應急監測方案研究與實踐:以金沙江白格堰塞湖為例[J].人民長江,2021,52(增1):73-76,84.

[78]ZHOU W,LI S,ZHOU Z,et al.InSAR observation and numerical modeling of the earth-dam displacement of Shuibuya dam (China)[J].Remote Sensing,2016,8(10):877.

[79]GUO R,LI S M,CHEN Y N,et al.Identification and monitoring landslides in Longitudinal Range-Gorge Region with InSAR fusion integrated visibility analysis[J].Landslides,2021,18:551-568.

[80]LIU B,HE K,HAN M,et al.Application of UAV and GB-SAR in Mechanism Research and Monitoring of Zhonghaicun Landslide in Southwest China[J].Remote Sensing,2021,13:1653.

[81]陳越,許惠平,陳華根.地理信息系統中矢量數據的三維可視化技術[J].同濟大學學報(自然科學版),2004,32(9):1197-1201.

[82]李云,范子武,吳時強,等.大型行蓄洪區洪水演進數值模擬與三維可視化技術[J].水利學報,2005,36(10):1158-1164.

[83]陸二男,張建清,胡志虎.唐家山堰塞湖導流渠地層結構的地球物理探測[J].人民長江,2008,39(22):55-56.

[84]吳學明,高才坤,王俊,等.綜合物探方法在紅石巖堰塞體探測中的應用研究[J].物探化探計算技術,2018,40(3):364-370.

[85]香天元.先進水文儀器在堰塞湖除險中的應用[J].中國水利,2008(13):65-67.

[86]周招,蔡耀軍,李建清,等.堰塞湖柔性防護措施控潰削峰技術研究[J].人民長江,2023,54(1):220-226.

[87]陳曉清,趙萬玉,高全,等.堰塞湖潰決人工結構體滯洪效果實驗研究[J].山地學報,2011,29(2):217-225.

[88]郭欽元.水域救援常用船艇對比研究[J].中國應急救援,2022(5):47-52.

[89]周志東,高靜,林振奎.堰塞湖應急處置施工技術概述[J].水利水電快報,2019,40(3):22-25.

[90]陳伯智.動力舟橋在特大洪澇災害搶險中的應用研究[J].四川水利,2022(增1):16-20.

[91]張金美,張宏建,曾智剛.制式應急橋梁在搶險救援中的應用研究[J].中國應急救援,2016(5):19-22.

[92]翁怡萌,蔡耀軍,王小波,等.堰塞湖大流量虹吸泄流關鍵技術與裝備研發[J].中國新技術新產品,2022(10):99-10.

(編輯:郭甜甜)

主站蜘蛛池模板: 中文无码毛片又爽又刺激| 国产精品9| 无码不卡的中文字幕视频| 91久久偷偷做嫩草影院| 欧美啪啪一区| 中文字幕久久波多野结衣| 男人的天堂久久精品激情| 另类专区亚洲| 久青草免费在线视频| 婷婷丁香在线观看| 国产亚洲欧美在线视频| 亚洲国产看片基地久久1024| 亚洲免费黄色网| A级毛片无码久久精品免费| 伊人婷婷色香五月综合缴缴情| 久久天天躁狠狠躁夜夜2020一| 亚洲VA中文字幕| 丁香婷婷综合激情| 91外围女在线观看| 中文字幕永久在线观看| 久久福利网| 久久久精品久久久久三级| 亚洲无线一二三四区男男| 欧美成人精品欧美一级乱黄| 国产成人h在线观看网站站| 国产美女自慰在线观看| 在线看片中文字幕| 秋霞国产在线| 五月天婷婷网亚洲综合在线| 亚洲欧美人成电影在线观看| 特级aaaaaaaaa毛片免费视频| 一区二区影院| 曰韩免费无码AV一区二区| 欧美中出一区二区| 亚洲人成网址| 国产精品乱偷免费视频| 露脸国产精品自产在线播| 国产九九精品视频| 国产成人91精品| 色噜噜综合网| 国产内射在线观看| 国产亚洲精品97在线观看| 免费Aⅴ片在线观看蜜芽Tⅴ | 狠狠操夜夜爽| 亚洲无线一二三四区男男| 在线观看免费国产| 婷婷综合色| 亚洲日韩Av中文字幕无码| 91最新精品视频发布页| 天堂亚洲网| 欧美自拍另类欧美综合图区| 久久网欧美| h网址在线观看| 视频二区国产精品职场同事| 免费看美女自慰的网站| 国产成人一区免费观看| 欧美在线黄| 欧美不卡视频在线| 久久国产精品麻豆系列| 真实国产乱子伦视频| 精品福利网| 欧美一区精品| 国产91视频免费观看| 亚洲最大在线观看| 亚洲Av综合日韩精品久久久| 在线精品视频成人网| 国产激情影院| 少妇高潮惨叫久久久久久| 成人午夜视频网站| 亚洲福利一区二区三区| 任我操在线视频| 欧美在线综合视频| 国产午夜小视频| 亚洲伊人久久精品影院| 国产偷国产偷在线高清| 天天色综合4| 萌白酱国产一区二区| 18禁高潮出水呻吟娇喘蜜芽| 亚洲高清中文字幕在线看不卡| 99久久精品国产麻豆婷婷| 四虎国产在线观看| 麻豆国产在线观看一区二区 |