• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于動力過程的冰湖潰決洪水侵蝕演化特征研究

    2023-10-05 14:17:44任思宇翟勝強葉曉華
    冰川凍土 2023年4期
    關鍵詞:支溝冰湖洪峰流量

    任思宇, 姜 亮, 翟勝強, 葉曉華, 郭 江

    (1. 四川省地質環(huán)境調查研究中心,四川 成都 610036; 2. 四川省地質災害防治工程技術研究中心,四川 成都 610036)

    0 引言

    受全球氣候變暖的影響,喜馬拉雅山地區(qū)冰川急劇退縮,冰川消融與后退不僅會影響冰川徑流變化[1],而且會導致由冰磧物攔擋的冰川湖的形成[2],由此引發(fā)的冰湖潰決災害也隨之增加[3]。冰湖潰決洪水具有突發(fā)性強、規(guī)模大、破壞力強和危害范圍廣等特點[4]。潰決洪水在演進的過程中,由于海拔落差大,潰決洪水對溝床沉積物進行沖刷侵蝕,進而形成泥石流。在泥石流形成、運動和堆積的過程中又會激發(fā)其他次生災害,從而形成災害鏈,因此致災能力非常強[5]。波曲河流域內“極高”危險的冰湖有10 個,“高”危險的有22 個[6],其下游為聶拉木縣及連接中國和尼泊爾的重要通道,該區(qū)域受冰湖影響危險性很高,波曲河流域中下游左側的一級支流章藏布于1964 年、1981 年和1983 年先后三次暴發(fā)大規(guī)模冰湖潰決災害[7]。其中1981 年的次仁瑪錯冰湖潰決最為嚴重[8],并引發(fā)泥石流、滑坡等次生災害。潰決洪水沖毀下游50 km 范圍內的中尼公路和交通設施,誘發(fā)了多處古滑坡復活,并沖毀曲鄉(xiāng)附近的建筑及尼泊爾境內的遜科西水電站,導致尼泊爾200 人死亡[9]。近年來流域氣候由干冷向濕熱轉化,章藏布冰川融水不斷補給次仁瑪錯冰湖,導致冰湖面積不斷擴張,除此之外章藏布冰川存在多組裂隙,融水匯入裂隙,增加了冰崩冰滑坡的風險,因此對冰湖潰決洪水防災減災志在必行。

    因此,要定量評估冰湖潰決洪水的危險性,就需要對其演化過程及動力特征進行量化分析,從而制定合理有效的防災減災措施[10]。近年來,有關冰湖潰決的研究正逐步從以分析冰湖潰決災害案例中的泥石流成因[11]、冰湖潰決模式[12]、機理分析[13]等為主的定量分析階段向以物理實驗[14]、數(shù)值模擬[8]為主的冰湖潰決洪水災害定量預測階段發(fā)展。

    目前國內外學者對冰湖潰決洪水的模擬演進進行了大量的研究,如楊瑞敏等[15]利用BREACH 模型和SMPDBK 模型估算和模擬來古冰湖潰決洪水并對其進行災害預警分析,BREACH 模型可以模擬潰決時壩體潰口特征及潰口處洪峰流量,SMPDBK 模型可以模擬潰決洪水的洪峰演進及洪水深度,但不能反映洪水演進過程中洪水水淹范圍、洪水速度,及洪水對溝床侵蝕等特征,在潰決洪水災害評估中受到了極大限制;舒有鋒[16]運用FLDWAV模型模擬了桑旺錯潰決洪水,模型可以考慮多種潰決洪水的啟動方式和壩體潰決破壞過程,但模型不能揭示潰決洪水在演進過程中流深的變化及洪水溝床的侵蝕特征;Sattar 等[17]利用HEC-RAS 一維和二維模型對南洛納克湖進行潰決模擬,模型可以考慮不同潰壩模式,可以采用一維模型模擬潰口處流量隨時間的變化趨勢,并將一維模型計算的結果導入二維模型,模擬潰決洪水的演進過程,模型能反映潰決洪水演進過程流深和流速變化,能較好地對冰湖潰決洪水災害進行定量分析,但模型難以反映潰決洪水侵蝕演化動力過程。歐陽朝軍采用連續(xù)介質開發(fā)的Massflow 數(shù)值平臺可以根據(jù)用戶需求對本構模型進行二次開發(fā),用戶可根據(jù)自身需求獲得相應模擬參數(shù),在滑坡、泥石流、洪水領域得到了廣泛應用。

    本文基于地表過程數(shù)值仿真軟件Massflow,根據(jù)實地調查和冰磧土力學參數(shù)對潰決洪水模型進行二次開發(fā),考慮了冰湖潰決洪水在演進的過程中對溝床及溝岸松散沉積物及坡積物的侵蝕。模型將泥沙輸移模型與水動力模型進行耦合,使模型更符合冰湖潰決洪水動力演化特征研究。通過建立冰湖潰決洪水動力模型,模擬反演次仁瑪錯冰湖潰決洪水的演進過程,著重分析潰決洪水在下游基礎設施處的洪峰流量、流速、流深、侵蝕、沉積及潰決洪水的傳播等特征,將模擬結果與實地調查結果進行對比分析,驗證了模型的可行性與適用性。對冰湖潰決洪水危險性進行定量評估,為冰湖潰決洪水防災減災提供了理論依據(jù)。

    1 研究區(qū)域概況

    1.1 地質背景

    次仁瑪錯冰湖(86° 03′35″~86° 4′23″ E,28°03′28″~28°4′10″ N)位于波曲河流域章藏布支溝(85°44′06″~85°59′12″ E,28°06′30″~28°19′18″ N)內。冰湖從西向東延伸,位于支溝的一個小冰川侵蝕的支流溝中,北部僅3.97 km2。在這個小支流中,由于冰川周期性退縮和發(fā)展形成了次仁瑪錯冰湖。冰湖距離中尼友誼橋約24 km,溝道平均縱比降為77.3‰。支溝干流長8.5 km,海拔從6 109 m 到3 168 m,溝道平均縱比降為176.0‰,匯水面積約50.5 km2。支溝內巖性主要為云母石英片巖,黑云母花崗巖片麻巖和第四紀沖積土。支溝內有六條主要支流,河床下切深度4~10 m。在雨季(6—9月),章藏布支溝干流流量為14 m3·s-1。圖1 顯示了次仁瑪錯冰湖和章藏布支溝的位置及物源條件,圖2為潛在冰崩冰滑坡航拍圖片。

    圖1 次仁瑪錯冰湖位置及章藏布支溝地層巖性及物源條件Fig. 1 The location of Cirenmaco glacial lake (a) and the formation lithology (b),and materials conditions of Zhangzangbu branch ditch (c)

    圖2 潛在冰崩冰滑坡Fig. 2 Potential ice avalanche and landslide

    研究區(qū)域分布兩條斷層如圖1 所示,北部斷層穿過章藏布。在這樣的地質背景下,第四紀層廣泛分布于溝谷,因此很容易發(fā)生滑坡和塌陷。此外,冰磧臺地分布在溝谷兩岸,冰磧土分布于冰舌前端,冰磧土顆粒粗,孔隙率高,顆粒級配級不均勻(圖3)。次仁瑪錯冰湖是由冰川作用形成的冰磧湖,下游由疏松的第四紀沉積土所覆蓋。

    圖3 章藏布冰磧土顆粒級配Fig. 3 Grain gradation of moraine in Zhangzangbu branch ditch

    章藏布支溝具有豐富的物源,為泥石流的形成和發(fā)育提供了有利條件,除此之外溝道兩岸分布眾多潛在滑坡,如707滑坡、樟木隧道滑坡和樟木口岸滑坡群,自尼泊爾4·25(2015 年4 月25 日)地震以來,古滑坡變形速率加快,次仁瑪錯冰湖一旦潰決,其導致的直接災害及次生災害將嚴重威脅著樟木口岸的正常運營(圖4~7)。

    圖4 707滑坡Fig. 4 Landslide of 707

    圖5 樟木隧道滑坡Fig. 5 Landslide of Zhangmu tunnel

    圖6 章藏布支溝物源Fig. 6 Materials of Zhangzangbu branch ditch

    圖7 樟木口岸滑坡群Fig. 7 Landslide group of Zhangmu port

    1.2 氣象條件

    聶拉木氣象站位于章藏布支溝溝口北部10 km2處,本文統(tǒng)計了氣象站1981—2017 年的年均降雨、年均氣溫和年均濕度數(shù)據(jù)如圖8 和圖9 所示,流域內呈年均氣溫逐漸升高,降雨逐漸減小的趨勢,氣象條件由濕冷轉為干熱,冰川融水加劇,冰湖面積加速擴張,增加了冰湖潰決洪水的風險。

    圖8 聶拉木氣象站1981—2017年的年均降雨和年均氣溫數(shù)據(jù)Fig. 8 Average annual rainfall and average annual temperature data of Nyalam Meteorological Station from 1981 to 2017

    圖9 聶拉木氣象站1981—2017年年均濕度數(shù)據(jù)Fig. 9 Annual average humidity data of Nyalam Meteorological Station from 1981 to 2017

    1.3 水文條件

    波曲流域氣候跨越了亞熱帶、暖溫帶、寒溫帶、高山苔原帶和雪山冰漠帶。受喜馬拉雅山脈地殼隆升的影響,河谷深切、河道彎曲,呈現(xiàn)典型的山區(qū)河流特點:河床坡降大,平均縱坡降達81‰。中、上游形成相對切深達1 000~1 500 m 的高山深谷,多呈“V”侵蝕型谷,一般寬在60~80 m;而下游峽谷深切,相對切深達1 500~3 000 m,其橫斷面多呈“V”型,兩岸懸崖峭壁,谷坡坡度大于35°,極有利于雨洪的集流和冰雪運移。波曲河流域年平均降水量582.9 mm,最大日降水量107.6 mm,雨季主要集中在6—8 月,降雪時間長達6個月以上。洪水期為5月中旬至9月中旬,7—8月為洪峰期,主河年均徑流流量31.7 m3·s-1,平均徑流深度約500 mm,平均徑流模數(shù)15.9 L·s-1·km-2,出境年徑流量1×108m3(友誼橋),主河和各主要支溝的徑流變幅較小。

    2 冰湖潰決洪水動力演化過程的模型構建

    2.1 模型構建

    模型采用全局笛卡爾坐標。模型分為上層運動流流體和底部可侵蝕沉積物,上層自由表面定義為Z1top點,上層和下層之間的邊界為Z1bot,h1=Z1top-Z1bot為流動深度。(u1,v1,w1)表示三個獨立坐標軸上深度平均的速度分量[18]。流動層質量和動量守恒方程可以表示如下[19]:

    式中:τ1zxbot和τ1zybot為上部流動層在底部靜態(tài)層上施加的基礎剪應力;E1bot為流體的侵蝕率;u1(Z1bot)和v1(Z1bot)為進入流動層材料的邊界速度。u1(Z1bot)和v1(Z1bot)用平均速度分量u1和v1表示。質點跳躍條件ρ1E1bot=ρ2E1top,其中E1top下層頂部的侵蝕率,ρ2為下層的密度[20]。

    對方程(1)~(3)進行推導,其表達式為:

    式中:E和D為侵蝕和沉積公式;P為溝床沉積物孔隙度。侵蝕和沉積方向垂直于基底表面,當采用第1 層的基礎剪應力τ1zxbot和τ1zybot時,在動量守恒方程右邊存在動量交換項ρ2u1(E-D)和ρ2v1(E-D)和Ouyang[18]。

    方程式的左邊去除ρ1,則質量和動量方程式為:

    侵蝕沉積方程通過刪減冗余項,得到如下方程:

    考慮到實際情況的復雜性,可將溝床沉積物的漲縮特征以及輕微影響項進行忽略,將溝床演化方程簡化為:

    式中:E和D為侵蝕和沉積速率;p為河床沉積物孔隙度。

    式中:α為經驗參數(shù),為近床濃度和深度平均濃度之間的差異,深度平均濃度由min[m,(1-p)/c];c1是流動中的固體濃度;m是經驗指數(shù);ws0是單個顆粒在平靜水中的沉降速度。它可以表示為[21]:

    式中:υ為水的運動黏度;d為沉積物粒徑;s=ρs/ρw-1;g是重力加速度,取9.8 m·s-2;ce為推移質泥沙輸移能力,可表示為:

    式中:θc為決定沉積物運動開始的希爾茲參數(shù)的臨界值[22];θ=(sgd)為屏蔽參數(shù);u*為摩擦速度。剪應力τ1zxbot和τ1zybot表示為[18]:

    其中n是曼寧粗糙度。

    2.2 求解方法

    將控制方程(8)至(17)進一步量化為向量模型:

    采用MacCormack-TVD 有限差分法來求解以上方程。與解傳統(tǒng)的淺水方程一樣,首先利用算子分裂法將方程(18)分解為兩個一維方程,表示如下:

    然后,(n+ 1)Δt時刻由以下方程求解:

    執(zhí)行預測步驟p:

    執(zhí)行修正步驟c:

    執(zhí)行平均步驟:

    式中:

    函數(shù)G的表達式為:

    函數(shù)φ(x)代表最小流體通量限制函數(shù)表達式為:

    式中:x代表自變量,變量C的表達式示為:

    式中:Cr為當?shù)氐膸炖蕯?shù)(Courant number)表達式如下:

    2.3 模型參數(shù)選取

    DEM 數(shù)據(jù)源自ALOS 衛(wèi)星數(shù)據(jù)。采用原始地形精度12.5 m×12.5 m的網格進行計算。通過實地調查和遙感影像(高分二號)對河道進行分析,模擬范圍河床糙率系數(shù)為0.07。冰磧土作為特殊的巖土材料,其性質與溝床沉積物區(qū)別大[23],由于潰決洪水侵蝕夾帶物源主要為冰磧土,因此我們采集了次仁瑪錯冰湖下游冰磧土體,進行土工實驗(圖3),由于流域范圍較大,不同河段溝床堆積碎屑顆粒級配相差巨大,難以在模型中實現(xiàn)區(qū)別對待,因此本文采用章藏布流域分布最廣泛的冰磧土作為模型的特征參數(shù),冰磧物容重為1.6~2.0 g·cm-3,潰決洪水容重1.0 g·cm-3。

    表1 次仁瑪錯冰湖潰決洪水數(shù)值模型參數(shù)Table 1 Numerical model parameters of the outburst flood of the Cirenmaco glacial lake

    3 模型驗證

    3.1 次仁瑪錯冰湖潰決誘因及危害

    次仁瑪錯冰湖曾在1964 年、1981 年和1983 發(fā)生過多次潰決,其中規(guī)模最大的一次為1981 年,冰川融水加速湖水位上升,導致冰磧壩體發(fā)生管涌潰決,潰決歷時60 min,最大潰決流量為16 000 m3·s-1,潰決洪水在壩體上沖刷了一條長50 m,底寬40~60 m 的沖溝,潰決洪水在章藏布支溝演進的過程中沖刷侵蝕約2×106m3固體物質,其中1.6×106m3的固體物質形成泥石流沖入波曲河,成為這次泥石流的主要物源。

    潰決洪水在演進過程中,對河床及溝岸殘坡積物進行沖刷側蝕,導致大規(guī)模的崩塌和滑坡,潰決洪水演進至友誼橋附近時,由于其運動受阻而雍高達20 m,比洪峰水位高約5 m,沖毀橋梁及東岸附近的全部建筑物,并引起東岸坡積、殘積層大規(guī)模的崩塌和滑坡[24]。泥石流過后,友誼橋上游段河谷中留下較大范圍的泥石流沉積。自友誼橋至水電站約8.3 km 的河段,除沿途繼續(xù)強烈的側向侵蝕外,大量泥石流物質在這里沉積,河床顯著升高,次仁瑪錯的潰決給下游造成了巨大的經濟損失。

    3.2 次仁瑪錯冰湖潰決洪水災害反演結果

    徐道明[24]對1981 年次仁瑪錯冰湖進行了實地調查,推測了潰口處的流量隨時間的演化特征和洪水在演進過程中洪峰流量隨距離的變化特征,因此在反演過程中以徐道明實地調查數(shù)據(jù)為依據(jù),采用在潰口處添加流量曲線作為模擬的啟動條件。

    潰決洪水在下游的演進過程中呈衰減趨勢,如圖10~11所示為實測洪峰流量、流深、流速與模擬洪峰流量、流深、流速隨距離的演化過程對比圖,模擬結果與實測結果誤差范圍小于20%,模擬結果與實測結果擬合度較好。

    圖10 模擬洪峰流量與實測流量Fig. 10 Simulated peak flow and measured flow

    圖11 模擬流深、流速與實測結果對比Fig. 11 Comparison of simulated flow depth and velocity with measured results

    4 預測分析

    次仁瑪錯冰湖從(0.10±0.08) km2(1988年)擴大至(0.35±0.04) km2(2021年),最大水深為(115±2) m,庫容為1.8×107m3,通過歷史遙感影像觀測,冰湖有明顯的擴張趨勢,流域氣象條件變化顯著,由濕冷轉向干熱(圖8~9),冰川融水流入冰川冰裂隙中,一方面導致冰川強度降低,另一方面起到潤滑作用,冰崩冰滑坡發(fā)生風險逐漸增加。

    由于冰湖在潰決的過程中,洪水流態(tài)及洪水對壩體的沖刷機理極為復雜。模型難以考慮潰口的演化過程,然而潰口假設法與實際情況偏差較大,得出的結果誤差率高。通過大量的文獻查閱,基于前人對大量冰湖潰決案例調查數(shù)據(jù)的基礎上,本文推導得到經驗公式法,擬合潰口處的流量曲線,以此作為冰湖潰決洪水啟動的參考方法。經驗公式法需要的參數(shù)較少,便于獲取,通過經驗公式(表2)計算潰決洪水的洪峰流量,潰口寬度及潰決歷時,擬合最危險工況(全潰)時的流量曲線,如圖12 所示,其中V為冰湖庫容,Hw為冰湖水深,K0為經驗參數(shù)取值為1。

    表2 冰湖潰決洪峰流量、潰口寬度和潰決歷時經驗公式Table 2 Empirical formula of glacial lake outburst peak discharge, outburst width and outburst duration

    圖12 潰口流量曲線Fig. 12 Flow curve at the breach

    圖13 為全潰條件下不同剖面點處的模擬洪峰流量[圖13(a)]、流深和流速[圖13(b)]。潰決洪水對潰口(No.1)進行下切侵蝕,侵蝕深度21~29 m(勘查最大侵蝕深度28 m),洪水模擬最大深度達16 m,洪峰流量時刻流速11.4 m·s-1,冰磧壩址(No.2)處地形開闊,洪水在該處發(fā)散,流速15 m·s-1、流深13 m,洪水裹挾的大顆粒土體在壩址處堆積,形成沖積扇,洪峰流量從1.80×104m3·s-1降至1.64×104m3·s-1。No.7~No.12 段,河床受冰川融水沖刷而下切侵蝕,河床及溝岸巖體破碎,河床縱比降較大,潰決洪水在該段對河床基底及坡岸進行強烈的沖刷侵蝕,潰決洪水逐漸演化為稀性泥石流,洪峰流量不斷增大。并在章藏布支溝溝口(No.12)和下游形成11 m高泥石流堰塞壩短暫堵塞主河。No.13~No.21 段,稀性泥石流逐漸進入峽谷,地形相對較緩,流深有所增加,流速有所放緩,洪峰流量過水斷面急劇降低,洪峰流量在演進過程中不斷衰減,在No.21~No.22段,泥石流對滑坡群坡腳進行侵蝕夾帶,洪峰流量有微弱的增長,下游溝床相對開闊,泥石流逐漸沉積,洪峰流量逐漸衰減,泥石流演進至水電站時洪峰流量為9 520 m3·s-1,對下游仍具有較強的破壞性。

    圖13 全潰條件下不同剖面點處的模擬洪峰流量(a)、流深和流速(b)Fig. 13 Simulated peak flow (a), flow depth and velocity (b) at different profile points under the condition of total collapse

    圖14 潰決洪水演進過程沿線侵蝕深度(a)和沉積深度(b)Fig. 14 Erosion depth (a)and siltation depth (b)along the evolution of the outburst flood

    707 老橋(No. 11)距離潰口6 440 m,潰決洪水到達該位置用時3 300 s,洪峰到達老橋用時4 500 s,洪峰流量為1.41×104m3·s-1,洪峰到達老橋時流深16.3 m,流速13 m·s-1,稀性泥石流對溝道左岸707滑坡坡腳進行沖刷側蝕[圖15(c)],最大侵蝕深度約為9 m,威脅707 滑坡穩(wěn)定性。稀性泥石流沖出章藏布支溝,沖擊波曲河右岸,對右岸土體進行沖刷侵蝕,侵蝕深度約9 m[圖15(c)],隨后稀性泥石流在章藏布支溝溝口下游沉積(No.12),形成泥石流堰塞壩[圖15(d)],短暫堵塞波曲河。對比No.11~No.12段的模擬結果和實地調查結果,模擬結果能較為準確地揭示潰決洪水流深、流速及溝床及溝岸松散沉積物的沖淤特征。因此能較為全面地對冰湖潰決洪水造成的危害進行定量評估。

    圖15 章藏布支溝溝口處洪峰流量時刻流深(a)、流速(b)和最終侵蝕深度(c)和沉積深度(d)Fig. 15 Flow depth (a), flow velocity (b), final erosion depth (c) and silting depth (d)at peak discharge time at the mouth of Zhangzangbu branch ditch

    滑坡群(No. 21~No. 23)位于樟木口岸附近,受尼泊爾2015 年4 月25 日8.1 級地震影響,古滑坡群復活,通過模擬結果可知洪峰到達該位置用時3.4 h,由圖16(a)和圖16(b)可以得出洪峰到達滑坡群時的最大流深25 m,最大流速為16.8 m·s-1,潰決洪水在該段對滑坡群坡腳進行沖刷掏蝕,最大侵蝕深度達12.0 m。流體中固體顆粒沉積于溝床,沉積深度高達11.0 m。潰決洪水對滑坡群坡腳沖刷側蝕,引起潛在滑坡群坡腳大規(guī)模崩塌,可能增加滑坡群變形速率或失穩(wěn)概率進而誘發(fā)更大規(guī)模的次生災害。

    圖16 滑坡群洪峰流量時刻流深(a)、流速(b)和最終侵蝕深度(c)和沉積深度(d)Fig. 16 Flow depth (a), flow velocity(b), final erosion depth (c) and siltation depth (d)of landslide group at flood peak discharge time

    水電站(No. 24)位于尼泊爾境內,距離潰口25.0 km,洪峰到達該位置用時4.8 h,洪峰流量為9 520 m3·s-1,流深17 m,流速8 m·s-1,圖17為水電站位置洪峰流量時刻的流體深度和流體速度及最終的侵蝕深度和淤埋深度。水電站高度約為15 m,稀性泥石流淤滿水電站,沉積深度最大約17 m,導致水電站失效。稀性泥石流漫過水電站繼續(xù)向下游演進,對下游基礎設施仍具有較高的破壞性。

    圖17 水電站洪峰流量時刻流深(a)、流速(b)和最終侵蝕深度(c)和沉積深度(d)Fig. 17 Flow depth (a), flow velocity (b), final erosion depth (c) and siltation depth (d)at peak discharge time of hydropower station

    5 討論

    冰湖潰決洪水在演化過程中由清水流轉換為高含沙洪水,再到泥石流[28]。前人采用BASEMENT[29],HEC-RAS[8],F(xiàn)LO-2D[30]等模型模擬高原冰湖潰決事件。然而,這些模型無法進行冰湖潰決洪水對沿途基底侵蝕夾帶的模擬。研究發(fā)現(xiàn)冰湖潰決誘發(fā)的泥石流的峰值流量比僅計算洪水的結果要高得多[28]。MASSFLOW 模型則可以模擬潰決洪水對基底的侵蝕夾帶效應,進而模擬冰湖潰決泥石流形成演化的動力過程,模型基于深度積分的二維數(shù)值模型,可根據(jù)實際需求對模型進行改進,進而考慮流體對基底的侵蝕效應[18]。除此之外,高精度的DEM 對于準確模擬冰湖潰決泥石流至關重要,它可以更加真實地反映潰決洪水在溝道的侵蝕和沉積狀況,由于模擬范圍較大,普通計算機難以實現(xiàn)高精度模擬,通過模擬結果與實地調查對比分析12.5 m 精度地形用于流域風險評估和危險性分析可以被接受。

    冰湖潰決事件模擬中的輸入?yún)?shù)的選擇具有極大的挑戰(zhàn)性。相對于經驗模型的參數(shù)校準,物理模型則利用質量和動量守恒等物理原理表示過程。然而,即使是復雜的物理模型在應用于真實事件時也有很大的不確定性,因此模型校準是必不可少的[31]。在現(xiàn)實應用中,受限冰湖潰決事件高質量的現(xiàn)場數(shù)據(jù)獲取,有時必須做出相應的妥協(xié)。模型校準的關鍵問題在于選擇合適的調整參數(shù)和定義一個能被接受的參數(shù)調整范圍[28],而敏感性分析可以識別模型校準中使用的敏感輸入?yún)?shù),這需要建模者證明標準參數(shù)值的偏差。

    在模型敏感性分析中,通過調整曼寧系數(shù)和潰決水量來觀察模擬結果的差異,其中曼寧系數(shù)控制著流體的紊流運動。敏感性分析表明,曼寧系數(shù)變化,泥石流淹沒面積、流深、流速等參數(shù)變化較小,在宏觀上這些誤差可以被接受。但流體在演化過程中的時間具有明顯差異,且泥石流侵蝕距離也顯著延長,但在泥石流平均流速上差異很?。?0]。同時,潰決水量的增加將加快了泥石流演進的速度。較大的曼寧系數(shù)將導致泥石流到達時間延遲,并增加泥石流深度和侵蝕距離,這與前人研究成果一致[32]。分析結果對區(qū)域風險評估、國土空間規(guī)劃及防災減災工作的開展具有一定的參考意義,但由于冰湖潰決洪水演進過程中的復雜性質和流動行為可能發(fā)生的變化,在解釋模型結果時仍然需要謹慎處理。

    6 結論

    通過資料收集、遙感解譯、野外調查等方法對西藏波曲河流域次仁瑪錯冰湖潰決洪水動力過程進行研究,對1981 年洪災進行模擬反演,并對次仁瑪錯冰湖再次潰決后的洪水演進動力過程進行預測,得到以下結論:

    (1)對次仁瑪錯1981年冰湖潰決洪水進行模擬反演,對比洪峰流量、流深、流速等特征參數(shù),得出模擬結果與實測結果具有良好的一致性,驗證了模型的適用性和可行性。

    (2)次仁瑪錯冰湖面積從自1988年以來急速擴張,最大水深為(115±2) m,庫容為(1.8×107) m3,流域氣候由濕冷轉向干熱,冰川融水流入冰川冰裂隙中,增加冰崩冰滑坡發(fā)生風險,進而導致冰湖潰決風險增加。

    (3)次仁瑪錯冰湖一旦再次潰決,洪水對下游冰磧物及松散坡積物和沉積物進行侵蝕夾帶,逐漸演化為稀性泥石流,對下游潛在滑坡坡腳進行沖刷側蝕,加速潛在滑坡變形。對下游友誼橋構成巨大威脅,稀性泥石流在水電處進行沉積,導致水電站失效。

    (4)潰決洪水在演進的過程中(No. 7~No. 12)對冰磧土及溝床沉積物進行沖刷侵蝕,洪水逐漸演化為泥石流,洪峰流量逐漸增強。稀性泥石流在(No.12~No.20)段,一方面對溝床及岸坡松散沉積物進行沖刷侵蝕,另一方面稀性泥石流在溝床開闊處形成沉積,基本上為沖淤平衡。稀性泥石流到達潛在滑坡群后(No.20~No.21),對潛在滑坡群坡腳進行沖刷側蝕,洪峰流量有微弱增長。泥石流到達水電站處,洪峰流量仍然有9 520 m3·s-1,對下游仍具有較強的危害性。

    (5)通過對次仁瑪錯冰湖潰決洪水的動力演化過程分析,潰決洪水對章藏布支溝松散物質進行侵蝕夾帶,逐漸演化為稀性泥石流,對下游707滑坡和樟木口岸滑坡群坡腳進行沖刷侵蝕,易誘發(fā)滑坡變形,導致次生災害的發(fā)生,除此之外,潰決洪水引發(fā)的泥石流將對中尼友誼橋和下游水電站的安全運營構成威脅。

    猜你喜歡
    支溝冰湖洪峰流量
    基于歷史邊界的喜馬拉雅山脈冰湖提取方法對比研究
    冰川凍土(2023年4期)2023-10-05 14:18:04
    按支溝穴緩解習慣性便秘
    點按支溝穴治岔氣
    跑步岔氣,點按支溝穴
    老友(2020年11期)2020-11-28 07:40:58
    退耕還林工程對渭河洪峰流量的影響
    冰湖奇觀
    空中之家(2017年11期)2017-11-28 05:28:04
    佛岡縣潖江流域年洪峰流量P-Ⅲ分布參數(shù)估算
    可可西里冰湖旁的白色帳篷
    中國三峽(2016年5期)2017-01-15 13:58:45
    大南川流域設計洪峰流量計算分析
    某特小流域設計洪峰流量計算分析
    他把我摸到了高潮在线观看| 亚洲全国av大片| 亚洲成人免费电影在线观看| 国产91精品成人一区二区三区| 国产一区二区激情短视频| 一区二区三区精品91| 久久久久国产一级毛片高清牌| 色播在线永久视频| 亚洲少妇的诱惑av| 欧美日韩乱码在线| 日本vs欧美在线观看视频| 亚洲在线自拍视频| 亚洲国产欧美日韩在线播放| 欧美黄色片欧美黄色片| 超色免费av| 久久草成人影院| 国产熟女xx| 久久欧美精品欧美久久欧美| 成年人免费黄色播放视频| 久久性视频一级片| 一级毛片精品| 韩国精品一区二区三区| 欧美日本中文国产一区发布| 国产av又大| 精品第一国产精品| 国产亚洲精品久久久久久毛片| 日韩欧美一区视频在线观看| 电影成人av| 欧美成人免费av一区二区三区| 亚洲一区二区三区不卡视频| 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 欧美色视频一区免费| 成人18禁高潮啪啪吃奶动态图| 久热爱精品视频在线9| 日本 av在线| 国产成人精品无人区| 狂野欧美激情性xxxx| 女人爽到高潮嗷嗷叫在线视频| 久久欧美精品欧美久久欧美| 琪琪午夜伦伦电影理论片6080| 国产欧美日韩一区二区精品| 日日夜夜操网爽| 亚洲五月色婷婷综合| 日本 av在线| 啪啪无遮挡十八禁网站| 在线播放国产精品三级| 一区二区三区激情视频| 久久精品影院6| 欧美成狂野欧美在线观看| 黑人猛操日本美女一级片| 国产成+人综合+亚洲专区| 91在线观看av| 大陆偷拍与自拍| av在线天堂中文字幕 | 在线看a的网站| 精品久久久久久久毛片微露脸| 久久性视频一级片| 欧美中文日本在线观看视频| 中文字幕高清在线视频| 亚洲人成网站在线播放欧美日韩| www.自偷自拍.com| 国产三级在线视频| 欧美激情极品国产一区二区三区| 18禁美女被吸乳视频| 男人操女人黄网站| 18禁黄网站禁片午夜丰满| 精品欧美一区二区三区在线| 欧美久久黑人一区二区| 成人精品一区二区免费| 两性午夜刺激爽爽歪歪视频在线观看 | 又大又爽又粗| 久久精品成人免费网站| 两性夫妻黄色片| 性少妇av在线| 国产97色在线日韩免费| svipshipincom国产片| 亚洲在线自拍视频| 一级毛片高清免费大全| 日韩国内少妇激情av| 美女国产高潮福利片在线看| 黄网站色视频无遮挡免费观看| 999久久久精品免费观看国产| 大型黄色视频在线免费观看| 国产精品久久久av美女十八| 国产又爽黄色视频| 黄片小视频在线播放| 午夜福利在线免费观看网站| av有码第一页| 两个人免费观看高清视频| 在线观看一区二区三区| 嫩草影视91久久| 午夜福利一区二区在线看| 满18在线观看网站| 人人澡人人妻人| 丝袜在线中文字幕| 欧美激情高清一区二区三区| 男女高潮啪啪啪动态图| 老熟妇仑乱视频hdxx| 五月开心婷婷网| 国产精品秋霞免费鲁丝片| 午夜91福利影院| 亚洲成人免费av在线播放| а√天堂www在线а√下载| 校园春色视频在线观看| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 国产欧美日韩一区二区精品| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精华国产精华精| 成人18禁高潮啪啪吃奶动态图| 久久婷婷成人综合色麻豆| 国产精品九九99| 热99re8久久精品国产| 欧美乱码精品一区二区三区| 亚洲激情在线av| tocl精华| 国产亚洲精品一区二区www| 在线观看一区二区三区激情| 一级片免费观看大全| 国产人伦9x9x在线观看| 韩国精品一区二区三区| 国产精品久久久久久人妻精品电影| 免费久久久久久久精品成人欧美视频| 日本黄色视频三级网站网址| 久久久精品欧美日韩精品| 狂野欧美激情性xxxx| 亚洲欧美精品综合一区二区三区| 国产一区二区在线av高清观看| 极品人妻少妇av视频| 高清欧美精品videossex| 午夜a级毛片| 在线观看免费日韩欧美大片| 国产成人免费无遮挡视频| 国产野战对白在线观看| 日韩大码丰满熟妇| 日日爽夜夜爽网站| 在线观看免费视频网站a站| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 最好的美女福利视频网| 18禁裸乳无遮挡免费网站照片 | 在线观看免费午夜福利视频| 国产亚洲精品一区二区www| 999久久久国产精品视频| 天堂影院成人在线观看| 亚洲色图 男人天堂 中文字幕| 欧美成狂野欧美在线观看| 亚洲欧洲精品一区二区精品久久久| 老司机深夜福利视频在线观看| 午夜两性在线视频| 国产成人系列免费观看| 另类亚洲欧美激情| 国内毛片毛片毛片毛片毛片| 国产熟女午夜一区二区三区| 免费在线观看视频国产中文字幕亚洲| www.999成人在线观看| 日韩一卡2卡3卡4卡2021年| 岛国在线观看网站| 久久久久久久午夜电影 | 99久久久亚洲精品蜜臀av| 午夜福利在线免费观看网站| 亚洲国产中文字幕在线视频| 妹子高潮喷水视频| 18禁美女被吸乳视频| 制服人妻中文乱码| 新久久久久国产一级毛片| 亚洲国产欧美一区二区综合| 亚洲av成人一区二区三| 亚洲精品一二三| 国产精品美女特级片免费视频播放器 | 女性被躁到高潮视频| 亚洲精品一二三| 精品国产一区二区久久| 在线观看66精品国产| 国产一区二区激情短视频| 成年女人毛片免费观看观看9| 国产精品亚洲一级av第二区| 国产有黄有色有爽视频| 国产成人av激情在线播放| 日本五十路高清| 老鸭窝网址在线观看| 免费在线观看日本一区| 男女做爰动态图高潮gif福利片 | 法律面前人人平等表现在哪些方面| 欧美成人午夜精品| 天堂中文最新版在线下载| 两个人看的免费小视频| 天天躁狠狠躁夜夜躁狠狠躁| a在线观看视频网站| 可以在线观看毛片的网站| 久久人人精品亚洲av| 麻豆av在线久日| 在线观看免费高清a一片| 在线观看一区二区三区激情| 国产精品美女特级片免费视频播放器 | 亚洲 欧美 日韩 在线 免费| 国内毛片毛片毛片毛片毛片| 黄色怎么调成土黄色| 国产三级在线视频| 一区在线观看完整版| 又黄又粗又硬又大视频| 国产一区二区激情短视频| 一级片'在线观看视频| 热re99久久精品国产66热6| 一个人免费在线观看的高清视频| 丰满迷人的少妇在线观看| 免费看a级黄色片| 中国美女看黄片| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| 亚洲国产精品一区二区三区在线| 美女福利国产在线| 精品一品国产午夜福利视频| 免费在线观看黄色视频的| 欧美不卡视频在线免费观看 | 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影 | 久久香蕉国产精品| 成人影院久久| 涩涩av久久男人的天堂| 日韩欧美三级三区| 老司机福利观看| 每晚都被弄得嗷嗷叫到高潮| 久久精品91无色码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲美女黄片视频| 午夜老司机福利片| 麻豆久久精品国产亚洲av | 老司机在亚洲福利影院| 亚洲九九香蕉| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 国产国语露脸激情在线看| 欧美激情久久久久久爽电影 | 国产国语露脸激情在线看| 一二三四在线观看免费中文在| 黑丝袜美女国产一区| 亚洲男人的天堂狠狠| 国产深夜福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美人与性动交α欧美软件| 午夜福利在线免费观看网站| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999在线| 欧美精品一区二区免费开放| 激情在线观看视频在线高清| 一夜夜www| 精品国内亚洲2022精品成人| 宅男免费午夜| 欧洲精品卡2卡3卡4卡5卡区| 欧美在线黄色| 97超级碰碰碰精品色视频在线观看| 动漫黄色视频在线观看| 狠狠狠狠99中文字幕| 国产1区2区3区精品| 国产主播在线观看一区二区| av片东京热男人的天堂| 好男人电影高清在线观看| 久久久久九九精品影院| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 精品国产乱子伦一区二区三区| 国产不卡一卡二| x7x7x7水蜜桃| 日韩免费av在线播放| 久久久久久久精品吃奶| 久久久久精品国产欧美久久久| 国产精品亚洲av一区麻豆| 亚洲一卡2卡3卡4卡5卡精品中文| 校园春色视频在线观看| 欧美最黄视频在线播放免费 | 亚洲国产精品sss在线观看 | www.www免费av| 久久热在线av| 在线观看免费视频日本深夜| 国产蜜桃级精品一区二区三区| 国产黄色免费在线视频| 亚洲一区高清亚洲精品| 看片在线看免费视频| 老熟妇乱子伦视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品二区激情视频| 黄频高清免费视频| 精品久久久久久电影网| 久久人妻av系列| 亚洲精品国产区一区二| 色综合站精品国产| 黄色片一级片一级黄色片| 亚洲国产欧美网| 9热在线视频观看99| 又紧又爽又黄一区二区| 两个人免费观看高清视频| 日韩大码丰满熟妇| 夜夜躁狠狠躁天天躁| 午夜免费成人在线视频| 久久99一区二区三区| 中文字幕最新亚洲高清| 身体一侧抽搐| 成在线人永久免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产乱人伦免费视频| 两个人免费观看高清视频| 香蕉国产在线看| 美女福利国产在线| 亚洲美女黄片视频| 国产精品 国内视频| xxx96com| 妹子高潮喷水视频| 两性夫妻黄色片| 成年版毛片免费区| 精品国产乱码久久久久久男人| 搡老乐熟女国产| 男人舔女人下体高潮全视频| 人妻久久中文字幕网| 最好的美女福利视频网| 亚洲欧美激情综合另类| 久久人人97超碰香蕉20202| 亚洲色图 男人天堂 中文字幕| 国产视频一区二区在线看| 精品免费久久久久久久清纯| 91精品国产国语对白视频| 18禁国产床啪视频网站| www.熟女人妻精品国产| 另类亚洲欧美激情| 成人18禁在线播放| 久久精品国产清高在天天线| 久久久国产一区二区| 午夜免费成人在线视频| 三级毛片av免费| 精品国产超薄肉色丝袜足j| 麻豆av在线久日| 97人妻天天添夜夜摸| 在线看a的网站| 国产精品久久久久成人av| av在线天堂中文字幕 | 免费av毛片视频| 欧美日韩亚洲国产一区二区在线观看| 一a级毛片在线观看| 99久久综合精品五月天人人| 男人舔女人下体高潮全视频| 日本三级黄在线观看| 在线观看免费视频网站a站| 人妻久久中文字幕网| 国产又爽黄色视频| 亚洲男人的天堂狠狠| 国产不卡一卡二| 婷婷丁香在线五月| 男女下面进入的视频免费午夜 | 欧美在线一区亚洲| 精品福利观看| 国产精品 国内视频| 老司机在亚洲福利影院| 久久精品成人免费网站| 动漫黄色视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 男人舔女人的私密视频| 无人区码免费观看不卡| 新久久久久国产一级毛片| aaaaa片日本免费| 婷婷精品国产亚洲av在线| x7x7x7水蜜桃| 欧美av亚洲av综合av国产av| 国产亚洲欧美精品永久| 国产高清激情床上av| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费观看网址| 午夜免费鲁丝| 免费高清视频大片| 美女福利国产在线| 精品国产乱码久久久久久男人| 美国免费a级毛片| 精品国产乱码久久久久久男人| 午夜福利影视在线免费观看| 日日干狠狠操夜夜爽| www日本在线高清视频| 亚洲欧美一区二区三区黑人| 欧美日韩av久久| 女警被强在线播放| xxxhd国产人妻xxx| 看黄色毛片网站| 一区福利在线观看| 97超级碰碰碰精品色视频在线观看| 国产成人av激情在线播放| 久久久国产欧美日韩av| 99热国产这里只有精品6| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 久久久久久久久免费视频了| 一区二区日韩欧美中文字幕| 极品人妻少妇av视频| 一级黄色大片毛片| 嫩草影院精品99| 久久国产精品人妻蜜桃| 亚洲精品久久成人aⅴ小说| 99热国产这里只有精品6| 亚洲精品在线美女| 欧美成人性av电影在线观看| 久久久久久久久中文| 国产免费现黄频在线看| 91精品国产国语对白视频| 99在线人妻在线中文字幕| 黑人猛操日本美女一级片| 日本欧美视频一区| 桃色一区二区三区在线观看| 日本免费一区二区三区高清不卡 | 美女高潮到喷水免费观看| 国产亚洲欧美98| 麻豆久久精品国产亚洲av | 又黄又爽又免费观看的视频| 长腿黑丝高跟| 免费久久久久久久精品成人欧美视频| 伊人久久大香线蕉亚洲五| 国产黄色免费在线视频| 欧美日韩国产mv在线观看视频| 午夜福利一区二区在线看| 久久草成人影院| 精品电影一区二区在线| 久久久久国内视频| 久久精品亚洲精品国产色婷小说| 日韩欧美三级三区| 熟女少妇亚洲综合色aaa.| 国产深夜福利视频在线观看| 大香蕉久久成人网| 国产精品乱码一区二三区的特点 | 成人三级做爰电影| 亚洲激情在线av| 精品国产乱子伦一区二区三区| 国产精品日韩av在线免费观看 | 久久国产乱子伦精品免费另类| 成熟少妇高潮喷水视频| 免费在线观看完整版高清| 日日爽夜夜爽网站| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 在线观看免费日韩欧美大片| 嫩草影视91久久| 亚洲成国产人片在线观看| 精品国产超薄肉色丝袜足j| 亚洲精品久久午夜乱码| 日本精品一区二区三区蜜桃| 国产成人免费无遮挡视频| 国产一区二区在线av高清观看| 热re99久久国产66热| 国产亚洲精品久久久久5区| 老司机在亚洲福利影院| 久久午夜亚洲精品久久| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久电影网| 国产精品香港三级国产av潘金莲| 级片在线观看| 99在线人妻在线中文字幕| 久久精品国产综合久久久| 亚洲欧美一区二区三区黑人| 亚洲欧美日韩另类电影网站| 99精品久久久久人妻精品| 国产三级在线视频| 多毛熟女@视频| 欧美黄色片欧美黄色片| 精品高清国产在线一区| 天堂中文最新版在线下载| 水蜜桃什么品种好| 久久午夜亚洲精品久久| 久久人人97超碰香蕉20202| 国产单亲对白刺激| 精品国产一区二区三区四区第35| 91大片在线观看| 琪琪午夜伦伦电影理论片6080| av电影中文网址| 久久久久久久精品吃奶| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 巨乳人妻的诱惑在线观看| 久久久久精品国产欧美久久久| 人人妻人人爽人人添夜夜欢视频| 精品高清国产在线一区| 大型av网站在线播放| 侵犯人妻中文字幕一二三四区| 久久久国产成人免费| 国产精品影院久久| 国产亚洲欧美精品永久| 国产欧美日韩精品亚洲av| 精品久久久精品久久久| 中亚洲国语对白在线视频| 99久久久亚洲精品蜜臀av| 女人被狂操c到高潮| 丰满的人妻完整版| 久久久久久久久中文| 精品福利永久在线观看| 国产高清视频在线播放一区| 精品午夜福利视频在线观看一区| 亚洲久久久国产精品| 国产熟女xx| 一级,二级,三级黄色视频| 老熟妇仑乱视频hdxx| 精品久久蜜臀av无| www.999成人在线观看| 日韩免费高清中文字幕av| 一边摸一边抽搐一进一小说| 亚洲av成人不卡在线观看播放网| 国产成人啪精品午夜网站| 在线天堂中文资源库| 亚洲专区字幕在线| 又大又爽又粗| 国产男靠女视频免费网站| www.熟女人妻精品国产| 露出奶头的视频| 亚洲色图av天堂| 久久香蕉激情| 国产不卡一卡二| 国产99久久九九免费精品| 精品免费久久久久久久清纯| 欧美日韩黄片免| 欧美丝袜亚洲另类 | 免费高清视频大片| 欧美人与性动交α欧美软件| 中文欧美无线码| 村上凉子中文字幕在线| 精品第一国产精品| 日韩欧美免费精品| 久久久久久免费高清国产稀缺| 欧美丝袜亚洲另类 | 99精品久久久久人妻精品| 日本欧美视频一区| 脱女人内裤的视频| 人成视频在线观看免费观看| 精品电影一区二区在线| 午夜成年电影在线免费观看| 色综合站精品国产| 亚洲专区国产一区二区| 久久这里只有精品19| 他把我摸到了高潮在线观看| 91在线观看av| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 日韩大尺度精品在线看网址 | 午夜久久久在线观看| av有码第一页| 变态另类成人亚洲欧美熟女 | 高清欧美精品videossex| 中文欧美无线码| 国产区一区二久久| 91麻豆精品激情在线观看国产 | 国产黄色免费在线视频| 男女午夜视频在线观看| 人人澡人人妻人| 51午夜福利影视在线观看| 老司机靠b影院| 欧美另类亚洲清纯唯美| 久久久国产精品麻豆| 又紧又爽又黄一区二区| 桃色一区二区三区在线观看| 99精品欧美一区二区三区四区| 欧美性长视频在线观看| 亚洲人成电影免费在线| 国产免费av片在线观看野外av| 看免费av毛片| 男女下面插进去视频免费观看| 国产精品综合久久久久久久免费 | 国产1区2区3区精品| 日韩人妻精品一区2区三区| 成人国产一区最新在线观看| 国产高清国产精品国产三级| 中文字幕另类日韩欧美亚洲嫩草| 久久婷婷成人综合色麻豆| 亚洲精品国产一区二区精华液| 亚洲av日韩精品久久久久久密| 在线观看66精品国产| 夜夜夜夜夜久久久久| 亚洲成av片中文字幕在线观看| 女性被躁到高潮视频| 无人区码免费观看不卡| 久久中文看片网| 亚洲va日本ⅴa欧美va伊人久久| 99香蕉大伊视频| 丁香欧美五月| 美女福利国产在线| 久久久久九九精品影院| 国产xxxxx性猛交| 日日爽夜夜爽网站| 又大又爽又粗| 天天影视国产精品| 国产1区2区3区精品| 丝袜人妻中文字幕| 亚洲第一av免费看| 一区二区三区国产精品乱码| 久久精品成人免费网站| 亚洲国产精品999在线| 日日摸夜夜添夜夜添小说| www.www免费av| 一级a爱片免费观看的视频| 一级毛片精品| 国产黄色免费在线视频| 国产在线观看jvid| 成人精品一区二区免费| 国产精华一区二区三区| 国产精品香港三级国产av潘金莲| 国产亚洲精品综合一区在线观看 | 亚洲国产精品sss在线观看 | 男女床上黄色一级片免费看| 国产精品国产高清国产av| 亚洲欧美激情综合另类| 亚洲av成人av| 变态另类成人亚洲欧美熟女 | 亚洲一区中文字幕在线| 精品久久久久久电影网| 亚洲精品美女久久久久99蜜臀| 亚洲精品一卡2卡三卡4卡5卡| 桃色一区二区三区在线观看| 在线天堂中文资源库| 国产精品亚洲一级av第二区| 婷婷丁香在线五月|