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

    基于Flow-R和FLO-2D的溝谷型泥石流危險(xiǎn)性評價(jià)

    2022-06-17 00:53:21林文周偉李婧楊莉唐澤
    人民長江 2022年5期
    關(guān)鍵詞:危險(xiǎn)區(qū)危險(xiǎn)性泥石流

    林文 周偉 李婧 楊莉 唐澤

    摘要:泥石流是中國西南山區(qū)頻發(fā)的地質(zhì)災(zāi)害之一,具有突發(fā)性、流速快以及破壞力強(qiáng)等特點(diǎn),往往造成巨大的人員傷亡和經(jīng)濟(jì)損失。地質(zhì)災(zāi)害危險(xiǎn)性評估是防災(zāi)減災(zāi)管理和防治環(huán)節(jié)中的有效措施之一,用于合理量化泥石流運(yùn)動堆積線路沿程泥石流災(zāi)害危險(xiǎn)性的空間分布特征。以四川省都江堰市龍池鎮(zhèn)麻柳溝為研究對象,選取坡度、平面曲率和上坡匯水面積作為溝谷型泥石流啟動區(qū)識別因子,基于Flow-R模型識別出的泥石流啟動點(diǎn),耦合FLO-2D模型模擬麻柳溝泥石流運(yùn)動堆積過程,將泥石流堆積深度、運(yùn)動速度和暴發(fā)頻率相結(jié)合分析泥石流的危險(xiǎn)性,并采用混淆矩陣對模擬結(jié)果進(jìn)行了精度評估。結(jié)果表明:基于Flow-R和FLO-2D耦合模型的麻柳溝泥石流的運(yùn)動堆積模擬結(jié)果與實(shí)際情況吻合度較高。模擬泥石流堆積區(qū)面積2.19萬m2,堆積體積4.92萬m3,泥石流危險(xiǎn)性由溝道中心向兩側(cè)逐漸降低,高危險(xiǎn)區(qū)占比42.3%,中危險(xiǎn)區(qū)占比23.1%,低危險(xiǎn)區(qū)占比34.6%。麻柳溝形成區(qū)、流通區(qū)溝道中散落堆積物體積大約為30萬m3,麻柳溝泥石流仍處于高度活躍階段,危險(xiǎn)性較高,今后在降雨和水動力充分的條件下大量松散物質(zhì)可能將被重新激活,一旦發(fā)生泥石流將有可能造成重大地質(zhì)災(zāi)害。

    關(guān) 鍵 詞:泥石流; 危險(xiǎn)性評價(jià); 泥石流啟動點(diǎn); Flow-R; FLO-2D; 都江堰市

    中圖法分類號: P642

    文獻(xiàn)標(biāo)志碼: A

    DOI:10.16232/j.cnki.1001-4179.2022.05.023

    0 引 言

    泥石流是中國西南山區(qū)頻發(fā)的地質(zhì)災(zāi)害之一,具有突發(fā)性、流速快以及破壞力強(qiáng)等特點(diǎn),往往會造成巨大的人員傷亡和經(jīng)濟(jì)損失。泥石流危險(xiǎn)性評價(jià)是災(zāi)害學(xué)重點(diǎn)研究的課題之一。泥石流活動規(guī)模是泥石流危險(xiǎn)性評價(jià)的一個(gè)重要指標(biāo)[1],它取決于控制流量觸發(fā)的運(yùn)動特征、泥石流沿程侵蝕和泥石流峰值流量等因素。FLO-2D是一種被廣泛用于模擬洪水的軟件,近年來也被用來模擬溝谷型泥石流的運(yùn)動和堆積,便于對溝谷型泥石流進(jìn)行危險(xiǎn)性評價(jià)。阮德修等[2]將FLO-2D與3DMine進(jìn)行耦合,模擬了尾礦庫潰壩時(shí)的泥石流流動過程,并預(yù)測了尾礦庫潰壩和半潰壩情況下的致災(zāi)程度。Chang等[3]通過實(shí)地調(diào)查,分析了汶川地震災(zāi)區(qū)泥石流的形成條件,運(yùn)用FLO-2D模型,模擬了泥石流運(yùn)動及其危害范圍。張建石[4]通過統(tǒng)計(jì)分析,得出震后多年間的物源演變規(guī)律及特征,并利用FLO-2D模型模擬在不同降水頻率下肖家溝泥石流的沖出特征。Nocentini[5]將FLO-2D模型與DAN-W模型進(jìn)行耦合,預(yù)測了研究區(qū)潛在泥石流的危險(xiǎn)性。黃靖玲等[6]將HEC-HMS模型與FLO-2D模型耦合,對馬貴河流域進(jìn)行了危險(xiǎn)性評估。如果要獲得較好的FLO-2D模型模擬效果,關(guān)鍵在于能找出相對準(zhǔn)確的泥石流啟動點(diǎn)。這就需要結(jié)合實(shí)地情況和專業(yè)經(jīng)驗(yàn),還應(yīng)考慮到一定的客觀因素。為了解決泥石流啟動點(diǎn)識別的問題,本文引入了一種識別溝谷型泥石流啟動點(diǎn)的軟件Flow-R。該軟件是一種基于GIS的重力災(zāi)害評價(jià)模型,可高效快速識別溝谷型泥石流啟動點(diǎn)[7]。

    四川省龍池鎮(zhèn)龍溪河流域自“5·12”地震后,流域內(nèi)松散物源量增加,地質(zhì)災(zāi)害活動頻繁,大量研究表明,泥石流災(zāi)害已成為汶川地震重災(zāi)區(qū)龍池鎮(zhèn)最大的地質(zhì)安全隱患[8]。2010年8月13日,流域在強(qiáng)降水條件下暴發(fā)泥石流災(zāi)害,造成巨大損失。本文以麻柳溝為研究對象,將FLO-2D模型與Flow-R模型耦合使用,來識別泥石流大多啟動點(diǎn),模擬泥石流淹沒區(qū)范圍,計(jì)算泥石流運(yùn)動的最大動量和最大堆積深度[9];同時(shí),結(jié)合流速、泥深與降水暴發(fā)頻率來探究麻柳溝的危險(xiǎn)性。研究成果可為麻柳溝的泥石流災(zāi)害的危險(xiǎn)分析、預(yù)警及發(fā)展規(guī)劃提供科學(xué)的參考依據(jù)。

    1 研究區(qū)概況

    麻柳溝位于四川省都江堰市龍池鎮(zhèn),溝口距離龍池鎮(zhèn)約0.65 km。溝口地理坐標(biāo)為103°33′20″E,31°03′55.11″N(見圖1)。麻柳溝流域面積約為0.93 km2,流域相對高差為913 m,主溝長度為2.3 km,平均縱坡降401‰。巨大的相對高差對流域內(nèi)的松散堆積物而言具有足夠的勢能條件,使得該溝易形成危害極大的高動能泥石流。

    研究區(qū)地處龍門山斷裂帶。2008年5月12日汶川縣映秀鎮(zhèn)發(fā)生8.0級大地震,其主發(fā)震斷裂映秀-北川斷裂橫貫研究區(qū),導(dǎo)致研究區(qū)內(nèi)巖層破碎嚴(yán)重,成為典型的構(gòu)造不穩(wěn)定區(qū)域。震后流域內(nèi)發(fā)生大量崩塌滑坡,為暴雨泥石流提供了充足的松散物質(zhì)來源。

    研究區(qū)屬于中亞熱帶濕潤季風(fēng)氣候區(qū),日照時(shí)間較少,降水充沛,陰雨天氣頻繁。多年平均降水量為1 134.8 mm(1987~2008年),5~9月降水量占全年降水量的80%。常年充分的降水使坡面松散堆積物直接運(yùn)移到溝道堆積,形成泥石流溝道物源。2010年8月13日龍池鎮(zhèn)發(fā)生強(qiáng)降水事件,最大1 h降水量達(dá)75.0 mm,連續(xù)3 h降水量達(dá)到150.0 mm,導(dǎo)致泥石流暴發(fā)。

    2 研究方法

    Flow-R模型將坡度、平面曲率、上坡匯水面積作為溝谷型泥石流啟動點(diǎn)的識別因子,然后將其轉(zhuǎn)化為Flow-R模型能識別的ASCII文件,通過對不同的因子設(shè)置閾值,識別出研究區(qū)的啟動點(diǎn)。

    基于DEM數(shù)據(jù)轉(zhuǎn)化為FLO-2D能識別的ASCII文件,結(jié)合現(xiàn)場實(shí)地調(diào)查進(jìn)行網(wǎng)格劃分。調(diào)查發(fā)現(xiàn)研究區(qū)泥石流溝道流動寬度大于10 m,所以將本次數(shù)值模擬網(wǎng)格劃分為10 m的網(wǎng)格滿足本次研究的需要,對每個(gè)網(wǎng)格賦值唯一的高程值、粗糙系數(shù)等基礎(chǔ)地形因子,將Flow-R識別的啟動點(diǎn)導(dǎo)入網(wǎng)格,輸入啟動點(diǎn)泥石流流量及流變參數(shù),模擬麻柳溝“8·13”特大泥石流運(yùn)動堆積過程,構(gòu)建混淆矩陣對運(yùn)動堆積范圍進(jìn)行精度驗(yàn)證。

    危險(xiǎn)性評價(jià)建模是將耦合模型結(jié)果獲取的淹沒區(qū)流速、泥深和降水頻率相結(jié)合,將泥石流溝劃分為高危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)和低危險(xiǎn)區(qū)。工作流程如圖2所示。

    2.1 泥石流啟動點(diǎn)識別

    坡度(地形條件)、松散物源可用性(物源條件)、輸水量(水源條件)是影響泥石流形成的關(guān)鍵因素,當(dāng)這些關(guān)鍵因素滿足一定條件時(shí)即可誘發(fā)泥石流[10]。Flow-R模型考慮了泥石流形成的三大影響因素,通過設(shè)置各個(gè)因子泥石流啟動識別閾值(詳見3.1節(jié)),將同時(shí)滿足閾值的柵格單元劃分為“啟動區(qū)”,將至少有一次不滿足閾值的柵格單元劃分為“不是啟動區(qū)”。

    2.2 泥石流運(yùn)動流通堆積模擬

    FLO-2D模型用于模擬泥石流的運(yùn)動過程和跳動特征[11]。模型以網(wǎng)格為單元計(jì)算泥石流流動,在軟件劃分出大小相同的規(guī)則網(wǎng)格后,將高程、粗糙系數(shù)等數(shù)值賦予給各網(wǎng)格。二維流動是通過運(yùn)動方程的數(shù)值積分和流體體積守恒來實(shí)現(xiàn)的,即網(wǎng)格中固體物質(zhì)和水的比值不斷變化,而總體積變化,且在模擬過程中,當(dāng)?shù)蜐舛攘黧w進(jìn)入網(wǎng)格遇到已停淤的高濃度流體時(shí),匯流處流體發(fā)生混合并繼續(xù)流動,控制方程如下。

    (1) 連續(xù)方程。

    ht+(uh)x+(vh)y=l(1)

    式中:l表示水力坡降;h表示泥深;v表示泥石流在y方向的速度;u表示泥石流在x方向的速度。

    (2) 運(yùn)動方程。

    Sfx=Scx-hx-VxgVxx-VygVyy-1gVxt(2)

    Sfy=Scy-hy-VygVyy-VygVyx-1gVyt(3)

    式中:Sfx和Sfy分別表示x、y方向的摩擦坡降;Scx和Scy分別表示x、y方向的底床坡降。

    (3) 流變方程。

    τ=τc+τmc+τt+τv+τd(4)

    式中:τc表示黏性屈服應(yīng)力;τv表示黏性剪應(yīng)力;τmc表示摩爾-庫倫剪應(yīng)力;τt表示紊流剪應(yīng)力;τd表示擴(kuò)散剪應(yīng)力。

    2.3 泥石流危險(xiǎn)性評價(jià)

    對單溝泥石流危險(xiǎn)性評價(jià),根據(jù)唐川等的研究,結(jié)合泥石流運(yùn)動速度和泥石流深度對泥石流強(qiáng)度進(jìn)行劃分[12](見表1)。在此基礎(chǔ)上,Chang等[3]以泥石流強(qiáng)度與暴發(fā)頻率相結(jié)合的方式進(jìn)行泥石流危險(xiǎn)性評價(jià),將研究區(qū)危險(xiǎn)性分為高、中、低3個(gè)等級(見圖3)。本文將耦合模型模擬結(jié)果結(jié)合泥石流運(yùn)動深度、速度以及降水頻率,對研究區(qū)的泥石流危險(xiǎn)性進(jìn)行評價(jià)。

    3 結(jié)果分析

    3.1 泥石流啟動點(diǎn)識別

    (1) 坡度。

    坡度是泥石流形成的主要因素。大多數(shù)地區(qū)發(fā)生泥石流的地形坡度大于15°[6],Ortigao等[13]認(rèn)為,泥石流啟動時(shí)的地形坡度介于20°~25°之間。麻柳溝坡度主要分布在20°~50°,占整個(gè)流域的86.41%。結(jié)合麻柳溝地形坡度特征,將泥石流啟動的地形坡度閾值定為20°。

    (2) 平面曲率。

    麻柳溝泥石流是典型的溝谷型泥石流,泥石流的物源一般分布在曲率為凹形的區(qū)域。平面曲率表示地形的凹凸變化,可用于識別泥石流溝道等負(fù)地形,間接體現(xiàn)泥石流的物源?,F(xiàn)有研究表明,泥石流啟動時(shí)的平面曲率閾值大多在-0.5/100 m-1到-2/100 m-1之間[14]。Delmonaco等[15]在平面曲率為-1/100 m-1的凹地形中識別出泥石流物源。麻柳溝負(fù)地形平面曲率統(tǒng)計(jì)結(jié)果表明,該溝平面曲率大于-3/100 m-1的比例達(dá)到84%。本文將研究區(qū)平面曲率閾值定為-3/100 m-1。

    (3) 匯水面積。

    上坡匯水面積,即匯流累積量,表示DEM中流經(jīng)每個(gè)柵格的流水累積量,旨在識別研究區(qū)內(nèi)任何活躍的河流、溝渠或暗流。Park等通過觀測實(shí)際發(fā)生的泥石流物源,構(gòu)建了泥石流源區(qū)坡度與匯水面積的閾值預(yù)測模型[16]。該模型如下:

    θ=15e32.34A-160.85(5)

    式中:A表示匯水面積,θ表示坡度。參考公式,通過設(shè)定100,500 m2和1 000 m2匯水面積閾值提取水系,結(jié)果表明1 000 m2閾值最符合實(shí)際情況。對麻柳溝流域進(jìn)行了匯水面積統(tǒng)計(jì),結(jié)果顯示,該溝78%的匯水面積在0~1 000 m2范圍內(nèi),因此,將研究區(qū)匯水面積率閾值定為1 000 m2。

    由圖5啟動點(diǎn)模擬結(jié)果可知:麻柳溝泥石流啟動區(qū)主要分布在溝道及溝道兩側(cè),啟動點(diǎn)坡度分布在20°~45°之間,其中,超過一半的啟動點(diǎn)分布于30°~40°之間,啟動點(diǎn)高程主要分布在1 100~1 300 m。溝頂細(xì)小匯水溝道為泥石流提供了豐富的水源,促使下部溝道沉積物轉(zhuǎn)化為泥石流。

    3.2 泥石流運(yùn)動流通堆積模擬

    參考水文地質(zhì)手冊,結(jié)合實(shí)地考察,泥石流啟動點(diǎn)主要分布在楊家溝和磯子溝。采用雨洪法分別計(jì)算了楊家溝和磯子溝處啟動點(diǎn)的洪峰流量。泥石流在運(yùn)動過程中挾帶有大量的松散物質(zhì),參考FLO-2D使用手冊[11],泥石流中水和沉積物總混合物的體積可以通過將水體積乘以放大系數(shù)(BF)來確定。放大系數(shù) BF 與泥石流的體積濃度密切相關(guān)[10],采用野外配漿實(shí)驗(yàn)獲取體積濃度為0.58,計(jì)算出麻柳溝啟動點(diǎn)流量。對洪峰流量采用簡單概化的五邊形方法來求取泥石流峰值流量過程線[17]。

    泥石流的屈服強(qiáng)度τy及黏滯參數(shù)η與堆積范圍相關(guān)[18]。τy、n和η可由下式計(jì)算:

    n=0.33C-0.15VeCv-0.15lnh(6)

    η=α1eβ1Cv(7)

    τy=α2eβ2Cv(8)

    式中:α1、α2、β1和β2為經(jīng)驗(yàn)系數(shù);Cv為泥石流體積濃度;h為泥深;n為曼寧系數(shù)。曼寧系數(shù)常被用來確定水和泥石流固體物質(zhì)之間的接觸關(guān)系。

    對麻柳溝進(jìn)行現(xiàn)場調(diào)查發(fā)現(xiàn),溝道內(nèi)植被稀疏,而且固體堆積物平均厚度在10 cm以上。曼寧系數(shù)可以間接表現(xiàn)流體的流變特征,根據(jù)公式(6)計(jì)算出的麻柳溝曼寧系數(shù)為0.12。層流阻滯系數(shù)K、α和β均為經(jīng)驗(yàn)系數(shù),一般通過查閱 FLO-2D 使用手冊或進(jìn)行專項(xiàng)試驗(yàn)獲取。本次研究,根據(jù)前人對經(jīng)驗(yàn)系數(shù)的取值[19]以及FLO-2D使用手冊[11]中的參數(shù)值列表,獲得了麻柳溝FLO-2D模擬的主要參數(shù)(見表3)。

    麻柳溝數(shù)值模擬結(jié)果如圖6所示。由圖6(a)可知:溝口的堆積面積為2.19×104 m2,占總流域面積的2.3%,最大堆積厚度為4.5 m,平均堆積深度為2.25 m,沖出總量為4.92×104 m3。研究區(qū)最大堆積深度達(dá)5.90 m,最小堆積深度為0.04 m,溝道中心泥深最大,向溝道兩側(cè)泥深減少,溝道中和堆積扇前段泥深較大,速度較快,溝道兩側(cè)和堆積扇翼端堆積深度較小,運(yùn)動速度較慢。

    利用混淆矩陣對模擬結(jié)果進(jìn)行精度評價(jià),并采用準(zhǔn)確率(A)和敏感度(S)指標(biāo)對模擬結(jié)果進(jìn)行評價(jià),準(zhǔn)確率和敏感度越高,表明模擬結(jié)果越好[20]。

    A=TP+TNN(9)

    S=TPTP+FN(10)

    式中:TP為實(shí)際淹沒區(qū)與模擬淹沒區(qū)均覆蓋的柵格數(shù);TN為實(shí)際淹沒區(qū)與模擬淹沒區(qū)均未覆蓋的柵格數(shù);FN為實(shí)際淹沒區(qū)覆蓋而模擬淹沒區(qū)未覆蓋的柵格數(shù);N為研究區(qū)柵格總數(shù)。根據(jù)麻柳溝2011年4月26日谷歌影像提取麻柳溝泥石流淹沒區(qū)范圍(見圖6(b)),根據(jù)公式(9)和公式(10)計(jì)算出的模擬結(jié)果準(zhǔn)確率為91.0%,敏感度為72.3%。表明耦合模型應(yīng)用于麻柳溝泥石流的運(yùn)動堆積模擬的結(jié)果與實(shí)際情況吻合度較高。

    3.3 泥石流危險(xiǎn)性評價(jià)

    模擬還原了麻柳溝“8·13”泥石流運(yùn)動堆積過程。由麻柳溝泥石流危險(xiǎn)性分區(qū)圖(見圖7)可知:高危險(xiǎn)區(qū)主要分布在溝道中心,面積約為58 300 m2,占危險(xiǎn)區(qū)總面積的42.3%。由溝道中心向兩側(cè)泥石流危險(xiǎn)性逐漸降低,中危險(xiǎn)區(qū)面積為31 200 m2,占總面積的23.1%;低危險(xiǎn)區(qū)面積為45 400 m2,占總面積的34.6%。高危險(xiǎn)區(qū)域分布范圍較大,低危險(xiǎn)和中危險(xiǎn)分布范圍相對較小。溝道兩側(cè)和堆積扇翼端主要為低危險(xiǎn)區(qū)域;泥石流形成區(qū)和堆積區(qū)邊緣主要為中危險(xiǎn)區(qū);高危險(xiǎn)區(qū)則主要集中在溝道中和堆積扇前段。在形成區(qū)和流通區(qū)的松散堆積體總量約30萬m3,在極端降水天氣下極易發(fā)生次生泥石流災(zāi)害,應(yīng)該加強(qiáng)防災(zāi)減災(zāi)工作和泥石流風(fēng)險(xiǎn)管理工作。

    4 結(jié) 論

    本文以麻柳溝為研究對象,基于Flow-R模型識別泥石流啟動點(diǎn)后耦合FLO-2D模型,用以模擬泥石流運(yùn)動;采用速度、泥石流堆積深度和降水暴發(fā)頻率相結(jié)合的方法,對單溝泥石流危險(xiǎn)性進(jìn)行了評價(jià),可以得出如下結(jié)論。

    (1) 麻柳溝泥石流啟動區(qū)絕大部分布在溝道及溝道兩側(cè),啟動點(diǎn)坡度分布在20°~45°之間,其中超過一半的啟動點(diǎn)分布于30°~40°之間,高程主要分布在1 100~1 300 m范圍內(nèi)。

    (2) 模擬泥石流堆積區(qū)面積為2.19萬m2,堆積體積為4.92萬m3。計(jì)算形成區(qū)、流通區(qū)溝道中散落堆積物體積大約為30萬m3,耦合模型應(yīng)用于麻柳溝模擬的準(zhǔn)確率為91.0%,敏感度為72.3%,模擬精度較高,可以為土地資源規(guī)劃、災(zāi)害防治設(shè)計(jì)提供參考依據(jù)。

    (3) 麻柳溝高危險(xiǎn)區(qū)分布范圍較大,占危險(xiǎn)區(qū)面積的42.3%,溝道中心危險(xiǎn)性較高,沿溝道中心向兩側(cè)危險(xiǎn)性降低,高危險(xiǎn)區(qū)主要集中在溝道中和堆積扇前沿,低危險(xiǎn)區(qū)主要分布于溝道兩側(cè)和堆積扇翼端,泥石流上溝段和堆積區(qū)邊緣主要為中危險(xiǎn)區(qū)。

    綜上所述,基于Flow-R和FLO-2D耦合模型將泥石流源區(qū)和泥石流運(yùn)動相結(jié)合,可以得出泥石流運(yùn)動堆積線路沿程災(zāi)害危險(xiǎn)性空間分布特征。根據(jù)危險(xiǎn)分布特征圖可知,麻柳溝泥石流仍處于高度活躍狀態(tài),今后在降水和水動力充分的條件下,大量的松散物質(zhì)可能被重新激活,一旦發(fā)生泥石流將可造成重大的地質(zhì)災(zāi)害。本文研究成果可為泥石流災(zāi)害監(jiān)測預(yù)警提供建設(shè)性考慮。本文研究未考慮物源量、巖性等其他因子對泥石流啟動點(diǎn)識別的影響,以及未考慮泥石流在運(yùn)動過程中側(cè)蝕和下蝕的影響,使得數(shù)值模擬的結(jié)果稍顯偏小,且較為保守,在今后的研究中需要對其進(jìn)行補(bǔ)充和完善。

    參考文獻(xiàn):

    [1] LUNA Q B,REMATRE A,VAN ASCH T W J,et al.Analysis of debris flow behavior with a one dimensional run-out model incorporating entrainment[J].Engineering Geology,2012,128:63-75.

    [2] 阮德修,胡建華,周科平,等.基于FLO-2D與3DMine耦合的尾礦庫潰壩災(zāi)害模擬[J].中國安全科學(xué)學(xué)報(bào),2012,22(8):150-156.

    [3] CHANG M,TANG C,VAN ASCH T W J,et al.Hazard assessment of debris flows inthe Wenchuan earthquake-stricken area,southwest China[J].Landslides,2017,14(5):1783-1792.

    [4] 張建石.汶川縣肖家溝泥石流物源演變及沖出規(guī)模研究[J].人民長江,2020,51(8):37-43.

    [5] NOCENTINI M,TOFANI V,GIGLI G,et al.Modeling debris flows in volcanic terrains for hazard mapping:the case study of Ischia Island(Italy)[J].Landslides,2015,12(5):831-846.

    [6] 黃靖玲,徐粒,周華真,等.基于歷史災(zāi)害情景再現(xiàn)的小流域山洪風(fēng)險(xiǎn)研究[J].災(zāi)害學(xué),2020,35(3):194-200.

    [7] HORTON P,JABOYEDOFF M,RUDAZ B,et al.Flow-R,a model for susceptibility mapping ofdebris flows and other gravitational hazards at a regionalscale[J].Natural Hazards and Earth System Sciences,2013,13(4):869-885.

    [8] 許強(qiáng).四川省8.13特大泥石流災(zāi)害特點(diǎn)、成因與啟示[J].工程地質(zhì)學(xué)報(bào),2010,18(5):596-608.

    [9] CHEN H,DADSON S,CHI Y.Recent rainfall-induced landslides and debris flows in northern Taiwan[J].Geomorphology,2018,77:112-125.

    [10] RICKENMANN D,ZIMMERMANN M.The 1987 debris flows in Switzerland:documentation and analysis[J].Geomorphology,1993,8(2-3):175-189.

    [11] FLO-2D Software Inc.FLO-2D,2-Dimensional Flood Routine Model Manual[Z].State of Arizona:FLD-2D Software,Inc.,2004.

    [12] 唐川,周鉅乾,朱靜,等.泥石流堆積扇危險(xiǎn)度分區(qū)評價(jià)的數(shù)值模擬研究[J].災(zāi)害學(xué),1994,9(4):7-13.

    [13] ORTIGAO J A R,KANJI M A.Landslide classification and risk management[M].Berlin:Springer,2004.

    [14] BLAHUT J,WESTEN C,STERLACCHINI S.Analysis of landslide inventories for accurate prediction of debris-flow source areas[J].Geomorphology,2010,119(1-2):36-51

    [15] DELMONACO G,LEONI G,MARGOTTINI C,et al.Large scale debris-flow hazard assessment:a geotechnical approach and GIS modeling[J].Natural Hazards and Earth System Sciences,2003,3:443-455.

    [16] PARK D W,LEE SR,VASU N N,et al.Coupled model for simulation of landslides and debris flows at local scale[J].Natural Hazards,2016,81:1653-1682.

    [17] 夏添.震區(qū)泥石流危險(xiǎn)性評價(jià)及預(yù)警減災(zāi)系統(tǒng)研究[D].成都:成都理工大學(xué),2013.

    [18] O’BRIEN J S,JULIEN P,F(xiàn)ULLERTON W.Two-dimensional water flood andmudflow simulation[J].Journal of Hydraulic Engineering,1993,119:244-261.

    [19] 程思.都江堰市龍溪河流域震后多溝同發(fā)泥石流危險(xiǎn)性及易損性研究[D].成都:成都理工大學(xué),2015.

    [20] CARRANZA E,CASTRO O T.Predicting lahar-inundation zones:case study in West Mount Pinatubo,Philippines[J].Natural Hazards,2006,37:331-372.

    (編輯:趙秋云)

    Gully debris flow hazard assessment based on Flow-R and FLO-2D coupling models

    LIN Wen,ZHOU Wei,LI Jing,YANG Li,TANG Ze

    (State Key Laboratory of Geohazard Prevention and Geoenvironment Protection,Chengdu University of Technology,Chengdu 610059,China)

    Abstract:

    Debris flow is one of the frequent geological disasters in southwest? mountainous areas of China.It is characterized by sudden occurrence,fast velocity and strong destructive force,which often causes huge casualties and economic losses.Risk assessment of geological hazards is one of the effective measures in disaster prevention and mitigation management and prevention,which can reasonably quantify the spatial distribution characteristics of debris flow hazard along the debris flow movement and accumulation route.Taking Maliu gully in Longchi Town,Dujiangyan City,Sichuan Province as the research object,slope,plane curvature and upslope catchment area were selected as the identification factors for the initiation area of gully debris flow.Based on the starting point of debris flow identified by Flow-R model,the FLO-2D model was coupled to simulate the movement and accumulation process of debris flow in Maliu gully.The debris flow hazard can be determined by combination of accumulation depth,velocity and outbreak frequency of debris flow.In addition,the accuracy of the simulation results was assessed by confusion matrix.The results showed that the simulation results of Maliu gully debris flow by Flow-R and FLO-2D coupling models was in good agreement with the actual conditions.The simulated debris flow accumulation area is 2.19×104m2,and the accumulation volume is 4.92×104m3.The debris flow risk gradually decreases from the center of the channel to both sides.High risk area accounts for 42.3%,medium risk area accounts for 23.1%,and low risk area accounts for 34.6%.In formation area and circulation area of Maliu gully,the volume of scattered accumulation materials in the channel is about 300000 m3.The Maliu gully debris flow is still in a highly active stage with high risk.In the future,under the condition of sufficient rainfall and hydrodynamic force,a large quantity of loose materials may be reactivated,which may cause major geological disasters once it occurs.

    Key words:

    debris flow;hazard assessment;starting point of debris flow;Flow-R;FLO-2D;Dujiangyan City

    猜你喜歡
    危險(xiǎn)區(qū)危險(xiǎn)性泥石流
    安徽省山洪危險(xiǎn)區(qū)動態(tài)化管理技術(shù)研究
    大科技(2022年28期)2022-07-25 08:02:24
    O-3-氯-2-丙烯基羥胺熱危險(xiǎn)性及其淬滅研究
    危險(xiǎn)性感
    輸氣站場危險(xiǎn)性分析
    基于AHP對電站鍋爐進(jìn)行危險(xiǎn)性分析
    泥石流
    雜文月刊(2018年21期)2019-01-05 05:55:28
    自升式鉆井平臺(JU2000E)的電氣防爆問題淺析
    “民謠泥石流”花粥:唱出自己
    海峽姐妹(2017年6期)2017-06-24 09:37:36
    泥石流
    機(jī)械班長
    国产精品无大码| 国产综合精华液| 一本一本综合久久| 国产91av在线免费观看| 人妻一区二区av| 亚洲国产高清在线一区二区三| 国产精品福利在线免费观看| av国产免费在线观看| 日本av免费视频播放| 日本午夜av视频| 99久久精品热视频| 亚洲综合色惰| 我的女老师完整版在线观看| 欧美精品国产亚洲| 在线观看国产h片| 久久久色成人| 男女无遮挡免费网站观看| 人妻制服诱惑在线中文字幕| 人人妻人人澡人人爽人人夜夜| 全区人妻精品视频| 你懂的网址亚洲精品在线观看| 岛国毛片在线播放| 欧美区成人在线视频| 男女边摸边吃奶| 高清av免费在线| 插阴视频在线观看视频| 亚洲精品乱码久久久v下载方式| 国产探花极品一区二区| 3wmmmm亚洲av在线观看| 亚洲色图综合在线观看| 一级黄片播放器| 日韩欧美一区视频在线观看 | 不卡视频在线观看欧美| 欧美bdsm另类| 99久久综合免费| 久久久久久久大尺度免费视频| 日日摸夜夜添夜夜添av毛片| 只有这里有精品99| 免费大片18禁| 免费看av在线观看网站| 看免费成人av毛片| 在线观看一区二区三区| 日日摸夜夜添夜夜爱| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 美女主播在线视频| 国产av国产精品国产| 国产无遮挡羞羞视频在线观看| 蜜桃在线观看..| 亚洲av中文av极速乱| 2018国产大陆天天弄谢| 亚洲国产精品成人久久小说| 免费看不卡的av| 日韩av在线免费看完整版不卡| 国产精品三级大全| 最近手机中文字幕大全| 日本欧美视频一区| 亚洲精品亚洲一区二区| 在线观看国产h片| 麻豆精品久久久久久蜜桃| 国产成人午夜福利电影在线观看| 汤姆久久久久久久影院中文字幕| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品古装| 色5月婷婷丁香| 免费人妻精品一区二区三区视频| 最近中文字幕高清免费大全6| 欧美丝袜亚洲另类| 亚洲欧美日韩东京热| 欧美bdsm另类| 天天躁日日操中文字幕| 国产日韩欧美在线精品| 亚洲av国产av综合av卡| 久久99热这里只频精品6学生| 婷婷色av中文字幕| .国产精品久久| 国产精品秋霞免费鲁丝片| 午夜福利视频精品| 亚洲精品一区蜜桃| 大陆偷拍与自拍| 日韩精品有码人妻一区| 久久久国产一区二区| 一级二级三级毛片免费看| 精品久久久久久久末码| 亚洲人与动物交配视频| 成年免费大片在线观看| 久久久久久久大尺度免费视频| 高清毛片免费看| 在线观看人妻少妇| 成人亚洲精品一区在线观看 | 久久久久国产精品人妻一区二区| 亚洲在久久综合| av在线观看视频网站免费| 哪个播放器可以免费观看大片| 丰满人妻一区二区三区视频av| 久久久精品94久久精品| 久久97久久精品| 国产精品久久久久久av不卡| 毛片女人毛片| 亚洲不卡免费看| 亚洲经典国产精华液单| 一边亲一边摸免费视频| 91精品国产九色| 国产精品一区二区三区四区免费观看| 少妇人妻 视频| 国产精品成人在线| 免费大片18禁| 老女人水多毛片| 欧美日韩综合久久久久久| 成人国产av品久久久| 美女国产视频在线观看| 人妻少妇偷人精品九色| 看免费成人av毛片| 亚洲精品色激情综合| 中文乱码字字幕精品一区二区三区| 80岁老熟妇乱子伦牲交| 日本午夜av视频| 久久午夜福利片| 精品久久久噜噜| 22中文网久久字幕| 深爱激情五月婷婷| 尤物成人国产欧美一区二区三区| 国产乱人视频| 久久国产亚洲av麻豆专区| 高清毛片免费看| 精品熟女少妇av免费看| 免费大片18禁| av线在线观看网站| 久久久亚洲精品成人影院| 免费大片18禁| 精品久久久噜噜| 久久久午夜欧美精品| 少妇裸体淫交视频免费看高清| 特大巨黑吊av在线直播| 免费观看在线日韩| 久久99热这里只频精品6学生| 国产在线一区二区三区精| 寂寞人妻少妇视频99o| 国产国拍精品亚洲av在线观看| 大片免费播放器 马上看| 国产黄色视频一区二区在线观看| 三级国产精品欧美在线观看| 九九久久精品国产亚洲av麻豆| videossex国产| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 观看av在线不卡| 青青草视频在线视频观看| 日韩一区二区三区影片| 成人一区二区视频在线观看| 精品国产露脸久久av麻豆| 美女主播在线视频| 久久午夜福利片| 亚洲av成人精品一区久久| 国产女主播在线喷水免费视频网站| 777米奇影视久久| 在现免费观看毛片| 尤物成人国产欧美一区二区三区| 99视频精品全部免费 在线| 亚洲精品色激情综合| 简卡轻食公司| 成人美女网站在线观看视频| 三级国产精品欧美在线观看| 一级毛片电影观看| 久久精品久久久久久久性| 国产精品不卡视频一区二区| 中文欧美无线码| 久久青草综合色| 中文字幕亚洲精品专区| 午夜视频国产福利| 婷婷色麻豆天堂久久| 精品一区在线观看国产| 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 国产精品欧美亚洲77777| 国产精品99久久99久久久不卡 | 老师上课跳d突然被开到最大视频| 国产在线男女| 久久女婷五月综合色啪小说| 永久网站在线| 国产精品一及| 国模一区二区三区四区视频| 国产精品久久久久成人av| 久久久久人妻精品一区果冻| 91精品伊人久久大香线蕉| 18禁动态无遮挡网站| 久久精品国产a三级三级三级| 久久久久久久久久成人| tube8黄色片| 蜜桃在线观看..| 黄色日韩在线| 女人久久www免费人成看片| 亚洲怡红院男人天堂| 亚洲综合精品二区| 精品视频人人做人人爽| 2022亚洲国产成人精品| 伦理电影免费视频| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 日韩人妻高清精品专区| 日韩中文字幕视频在线看片 | 免费大片黄手机在线观看| 日韩中字成人| 亚洲色图综合在线观看| 亚洲美女搞黄在线观看| 插逼视频在线观看| 亚洲综合色惰| 亚洲国产最新在线播放| 国产日韩欧美亚洲二区| 国产有黄有色有爽视频| 亚洲国产精品专区欧美| 高清不卡的av网站| 国产91av在线免费观看| 99视频精品全部免费 在线| 亚州av有码| 水蜜桃什么品种好| 中文字幕久久专区| 插逼视频在线观看| 日韩大片免费观看网站| 大香蕉97超碰在线| 另类亚洲欧美激情| 国产黄片美女视频| 寂寞人妻少妇视频99o| 午夜激情久久久久久久| 热re99久久精品国产66热6| 性高湖久久久久久久久免费观看| 亚州av有码| 午夜激情福利司机影院| 欧美国产精品一级二级三级 | 热99国产精品久久久久久7| 成人美女网站在线观看视频| 亚洲欧美精品自产自拍| 噜噜噜噜噜久久久久久91| 中文字幕精品免费在线观看视频 | 男人爽女人下面视频在线观看| 国产 一区 欧美 日韩| 蜜桃在线观看..| 天堂中文最新版在线下载| 婷婷色综合大香蕉| 日日撸夜夜添| 青春草亚洲视频在线观看| 能在线免费看毛片的网站| 91精品国产九色| 少妇的逼水好多| 九九在线视频观看精品| 久久国产乱子免费精品| 日韩av免费高清视频| 卡戴珊不雅视频在线播放| 国产淫片久久久久久久久| 国产91av在线免费观看| 亚洲丝袜综合中文字幕| 国产在线免费精品| 1000部很黄的大片| 深夜a级毛片| 国产男女内射视频| 亚洲成人手机| 久久久久性生活片| 国产精品麻豆人妻色哟哟久久| 欧美+日韩+精品| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 狂野欧美激情性xxxx在线观看| 黄色视频在线播放观看不卡| 成人综合一区亚洲| 国产白丝娇喘喷水9色精品| 自拍偷自拍亚洲精品老妇| av又黄又爽大尺度在线免费看| 国精品久久久久久国模美| 夜夜骑夜夜射夜夜干| 性色av一级| 又黄又爽又刺激的免费视频.| 91狼人影院| 精品人妻偷拍中文字幕| 国产一区亚洲一区在线观看| 午夜福利视频精品| 男人狂女人下面高潮的视频| 人人妻人人看人人澡| 午夜福利在线观看免费完整高清在| 最黄视频免费看| 国产亚洲最大av| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| 亚洲精品一二三| 永久免费av网站大全| 菩萨蛮人人尽说江南好唐韦庄| 在线观看国产h片| 国产成人aa在线观看| 一个人免费看片子| 色综合色国产| 国产一区二区三区av在线| 国产精品三级大全| 五月天丁香电影| 又黄又爽又刺激的免费视频.| 亚洲,欧美,日韩| 两个人的视频大全免费| 色婷婷av一区二区三区视频| 人人妻人人添人人爽欧美一区卜 | 高清在线视频一区二区三区| 日韩成人伦理影院| 精品人妻偷拍中文字幕| 2018国产大陆天天弄谢| 99久久精品国产国产毛片| 久久精品国产亚洲av天美| 十八禁网站网址无遮挡 | 亚洲av二区三区四区| 国产 一区精品| 久久久精品94久久精品| 麻豆成人午夜福利视频| 国产日韩欧美亚洲二区| 亚洲av二区三区四区| 精品国产乱码久久久久久小说| 深爱激情五月婷婷| 久久精品国产亚洲网站| av专区在线播放| 成人毛片60女人毛片免费| 国产精品av视频在线免费观看| 欧美zozozo另类| 99热这里只有是精品在线观看| 国产在线男女| 22中文网久久字幕| 日本猛色少妇xxxxx猛交久久| 婷婷色综合www| 国产在线一区二区三区精| 少妇高潮的动态图| 日韩伦理黄色片| 国产精品女同一区二区软件| 人妻系列 视频| 小蜜桃在线观看免费完整版高清| 男女免费视频国产| 午夜激情久久久久久久| 久久久久网色| 亚洲av成人精品一区久久| 精品熟女少妇av免费看| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 老熟女久久久| 黑丝袜美女国产一区| 国产成人91sexporn| 亚洲高清免费不卡视频| 久久久久久久久久久丰满| 国产精品一区二区三区四区免费观看| 亚洲精品久久久久久婷婷小说| 多毛熟女@视频| 国产黄色视频一区二区在线观看| 我的老师免费观看完整版| 各种免费的搞黄视频| 毛片女人毛片| 成人综合一区亚洲| 亚洲av欧美aⅴ国产| 国产成人精品婷婷| 我要看黄色一级片免费的| 精品人妻一区二区三区麻豆| 精品国产一区二区三区久久久樱花 | 国产精品蜜桃在线观看| 99国产精品免费福利视频| 亚洲成人中文字幕在线播放| 最近中文字幕高清免费大全6| 一级av片app| 日韩,欧美,国产一区二区三区| 99热国产这里只有精品6| 91精品国产国语对白视频| 亚洲不卡免费看| 人体艺术视频欧美日本| 中国三级夫妇交换| 国产老妇伦熟女老妇高清| 久久久精品免费免费高清| 亚洲精品自拍成人| 国产精品女同一区二区软件| 极品教师在线视频| h日本视频在线播放| 久久久成人免费电影| 国产在视频线精品| 国产淫语在线视频| 亚洲av中文av极速乱| av卡一久久| 亚洲精华国产精华液的使用体验| 草草在线视频免费看| 国产精品一区www在线观看| 成人亚洲精品一区在线观看 | 美女高潮的动态| 国产免费福利视频在线观看| 亚洲欧美日韩卡通动漫| 大又大粗又爽又黄少妇毛片口| 日韩伦理黄色片| 久久久久国产精品人妻一区二区| 成人亚洲精品一区在线观看 | 国产精品久久久久久精品古装| 色网站视频免费| 高清av免费在线| 99热这里只有是精品在线观看| 久久精品熟女亚洲av麻豆精品| 汤姆久久久久久久影院中文字幕| 国产亚洲最大av| 少妇 在线观看| 老司机影院毛片| 精品一区在线观看国产| 亚洲人成网站在线播| 22中文网久久字幕| 国产视频首页在线观看| 亚洲成人一二三区av| 色综合色国产| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| 男女无遮挡免费网站观看| 亚洲精品日韩在线中文字幕| 超碰av人人做人人爽久久| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人 | 赤兔流量卡办理| 黄片无遮挡物在线观看| 91精品国产国语对白视频| 又爽又黄a免费视频| 午夜福利高清视频| 久久影院123| 在线亚洲精品国产二区图片欧美 | 久久国产亚洲av麻豆专区| 观看免费一级毛片| 欧美精品人与动牲交sv欧美| 视频中文字幕在线观看| av黄色大香蕉| 边亲边吃奶的免费视频| 亚洲国产精品一区三区| 99热这里只有精品一区| 国产免费一级a男人的天堂| 午夜免费鲁丝| 免费黄频网站在线观看国产| 国产乱人视频| 一区二区三区四区激情视频| 尤物成人国产欧美一区二区三区| 日本-黄色视频高清免费观看| av一本久久久久| 特大巨黑吊av在线直播| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 国产v大片淫在线免费观看| 亚洲电影在线观看av| 成人18禁高潮啪啪吃奶动态图 | 免费少妇av软件| 婷婷色麻豆天堂久久| 欧美成人a在线观看| 小蜜桃在线观看免费完整版高清| 最黄视频免费看| 在现免费观看毛片| 欧美+日韩+精品| www.色视频.com| 国产中年淑女户外野战色| 成人免费观看视频高清| 有码 亚洲区| 亚洲人成网站在线播| 日韩欧美 国产精品| 亚洲国产最新在线播放| av福利片在线观看| 亚洲av电影在线观看一区二区三区| 观看免费一级毛片| 久久精品久久久久久噜噜老黄| 久久精品国产亚洲网站| 亚洲av日韩在线播放| videossex国产| 国产成人freesex在线| 国产亚洲最大av| 国产免费又黄又爽又色| 少妇人妻 视频| 在线亚洲精品国产二区图片欧美 | 国产精品国产三级国产专区5o| 九九在线视频观看精品| 搡老乐熟女国产| videos熟女内射| 国产午夜精品久久久久久一区二区三区| 亚洲欧美清纯卡通| 黄色配什么色好看| 97在线视频观看| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 免费观看在线日韩| 激情五月婷婷亚洲| 另类亚洲欧美激情| 我要看日韩黄色一级片| 丝瓜视频免费看黄片| 国产精品一及| 中文字幕制服av| 亚洲av日韩在线播放| 久久97久久精品| 国产精品99久久久久久久久| 中国国产av一级| 亚洲国产毛片av蜜桃av| 国产午夜精品久久久久久一区二区三区| 我要看黄色一级片免费的| 亚洲在久久综合| 亚洲欧洲日产国产| 亚洲国产成人一精品久久久| 22中文网久久字幕| av国产久精品久网站免费入址| 国产精品人妻久久久久久| 色婷婷av一区二区三区视频| 色视频www国产| 永久免费av网站大全| 午夜福利视频精品| kizo精华| 亚洲伊人久久精品综合| 亚洲欧美清纯卡通| 少妇精品久久久久久久| 亚洲欧美日韩无卡精品| 久久久久精品性色| 久久女婷五月综合色啪小说| 丰满迷人的少妇在线观看| 久久久久久久久久久丰满| 国产精品福利在线免费观看| 中文乱码字字幕精品一区二区三区| 伦理电影大哥的女人| av一本久久久久| 涩涩av久久男人的天堂| 亚洲美女黄色视频免费看| 在现免费观看毛片| av福利片在线观看| 我要看日韩黄色一级片| a级毛色黄片| 国产有黄有色有爽视频| 国产精品爽爽va在线观看网站| 国产亚洲一区二区精品| 欧美日韩一区二区视频在线观看视频在线| 亚洲成色77777| 天美传媒精品一区二区| 秋霞伦理黄片| 亚洲av男天堂| 国产成人a区在线观看| 国产男人的电影天堂91| 国产av一区二区精品久久 | 十分钟在线观看高清视频www | 91精品国产九色| 少妇人妻久久综合中文| 日本黄色片子视频| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线 | 日韩欧美一区视频在线观看 | 多毛熟女@视频| 97超碰精品成人国产| 欧美性感艳星| 亚洲在久久综合| 久久青草综合色| 免费av不卡在线播放| 2021少妇久久久久久久久久久| 国产女主播在线喷水免费视频网站| 亚洲四区av| 高清黄色对白视频在线免费看 | 久久久亚洲精品成人影院| 2021少妇久久久久久久久久久| 日韩欧美 国产精品| 亚洲精品第二区| 人妻夜夜爽99麻豆av| 国产成人aa在线观看| 精品人妻视频免费看| 99热这里只有是精品在线观看| 欧美日韩在线观看h| 91精品国产九色| 精品国产三级普通话版| 欧美成人午夜免费资源| 视频区图区小说| 精品人妻一区二区三区麻豆| 亚洲自偷自拍三级| 人妻 亚洲 视频| 综合色丁香网| 夫妻性生交免费视频一级片| 日韩免费高清中文字幕av| 男的添女的下面高潮视频| 人妻系列 视频| 精品国产乱码久久久久久小说| 99久久人妻综合| 久久久国产一区二区| 国产成人一区二区在线| 欧美zozozo另类| 成年美女黄网站色视频大全免费 | av又黄又爽大尺度在线免费看| 久热这里只有精品99| 国产精品精品国产色婷婷| 亚洲色图av天堂| 久久久色成人| 只有这里有精品99| 丝瓜视频免费看黄片| 久久久国产一区二区| 国产伦在线观看视频一区| 久久国产亚洲av麻豆专区| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 国产亚洲5aaaaa淫片| 欧美成人午夜免费资源| 国产成人精品婷婷| 麻豆成人av视频| 国语对白做爰xxxⅹ性视频网站| 中文字幕av成人在线电影| 婷婷色av中文字幕| 国产精品一及| 一级毛片aaaaaa免费看小| 一级二级三级毛片免费看| 久久久久久久久久久免费av| 女性被躁到高潮视频| 国产真实伦视频高清在线观看| 国产深夜福利视频在线观看| 日韩av免费高清视频| 亚洲在久久综合| 美女cb高潮喷水在线观看| 亚洲精品aⅴ在线观看| 一级黄片播放器| 看非洲黑人一级黄片| 乱码一卡2卡4卡精品| 成人18禁高潮啪啪吃奶动态图 | 最后的刺客免费高清国语| 中文字幕av成人在线电影| 九九爱精品视频在线观看| 少妇的逼好多水| 国产精品一区二区三区四区免费观看| 国产老妇伦熟女老妇高清| av国产免费在线观看| 男人添女人高潮全过程视频| 国产黄片视频在线免费观看|