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

    基于FLO-2D數(shù)值模擬的泥石流流動與堆積影響因素研究*

    2016-07-06 01:15:49梁鴻熙
    工程地質(zhì)學(xué)報 2016年2期
    關(guān)鍵詞:屈服應(yīng)力黏性泥石流

    梁鴻熙 尚 敏 徐 鑫

    (三峽大學(xué)土木與建筑學(xué)院宜昌443002)

    基于FLO-2D數(shù)值模擬的泥石流流動與堆積影響因素研究*

    梁鴻熙尚敏徐鑫

    (三峽大學(xué)土木與建筑學(xué)院宜昌443002)

    摘要泥石流和普通的山體滑坡不同,它是發(fā)生在山區(qū)的一種含有砂土和石塊的暫時性水流,具有賓漢體的性質(zhì)和運(yùn)動阻塞雙重特性,國內(nèi)外眾多學(xué)者在流動學(xué)實(shí)驗(yàn)及流動學(xué)模型的基礎(chǔ)上,對數(shù)值模擬方法進(jìn)行了大量研究,通過泥石流的災(zāi)害范圍、運(yùn)動速度、運(yùn)動時間與實(shí)際情況進(jìn)行對比,開發(fā)出適用于泥石流預(yù)報的數(shù)值模擬程序,F(xiàn)LO-2D就是其中之一。工程實(shí)踐證明,在對泥石流組成物質(zhì)的流動特性研究基礎(chǔ)上,對泥石流移動及沉積特性進(jìn)一步研究是非常重要的。本文運(yùn)用FLO-2D流動模型數(shù)值模擬方法對泥石流的移動特性進(jìn)行模擬,以此分析泥石流流動及堆積特性與黏性系數(shù)和屈服應(yīng)力的關(guān)系。通過分析得到:隨著黏性系數(shù)的增大,泥石流流速呈非線性減少; 屈服應(yīng)力作為影響泥石流發(fā)生及停止的因素,隨著屈服應(yīng)力的增大,流深也隨著增大,但在泥石流移動特性中,屈服應(yīng)力對泥石流流速的影響并不明顯; 隨著黏性系數(shù)的增加,泥石流沖擊力減小,而隨著屈服應(yīng)力的增加,泥石流的沖擊力則增大。

    關(guān)鍵詞數(shù)值模擬FLO-2D黏性系數(shù)屈服應(yīng)力

    0引言

    泥石流和普通的山體滑坡不同,它是發(fā)生在山區(qū)的一種包含泥砂、塊石的固液兩相流體,其介于塊體重力運(yùn)動和水流等液體水力運(yùn)動之間,呈黏性層流或稀性紊流狀態(tài)。據(jù)有關(guān)學(xué)者的大量研究表明,泥石流具有賓漢體的性質(zhì)和運(yùn)動阻塞雙重特性。O’Brienetal.(1985)、Iversonetal.(1997)、Imranetal.(2001)在對泥石流進(jìn)行分析時,根據(jù)泥石流所具有的移動特性,把它作為一種流體來進(jìn)行分析研究。由于工程地質(zhì)、氣象等要素的綜合作用影響,使得泥石流運(yùn)動模型的研究變得困難。同時,由于泥石流運(yùn)動速度快和移動距離長及受災(zāi)范圍廣,使得泥石流的預(yù)報工作難度也非常大。但是,泥石流往往突然爆發(fā),來勢兇猛,而且其能量巨大,破壞性強(qiáng),還具有運(yùn)動過程的直進(jìn)性、陣流性、非定常性、不均勻性和侵蝕性等特點(diǎn)。有時形成危害較大的陣流,嚴(yán)重威脅人民群眾的生命財(cái)產(chǎn)安全。為了使泥石流災(zāi)害區(qū)的損失最小化,所以對災(zāi)害區(qū)進(jìn)行提前災(zāi)害預(yù)報是十分必要的。由于泥石流流動與堆積特性直接決定泥石流的破壞及影響范圍,也是影響預(yù)測預(yù)報的主要因素。因此,對兩者的研究對泥石流的預(yù)測預(yù)報是至關(guān)重要的,有較強(qiáng)的實(shí)際應(yīng)用價值。

    近年來,泥石流危險性評價得以發(fā)展,研究成果的實(shí)踐意義也得到很大提高。研究泥石流的主要方法有野外調(diào)查、數(shù)值模擬以及實(shí)驗(yàn)研究等多種方法。由于數(shù)值模擬方法的成本低,而且重復(fù)性強(qiáng),在近年來被廣泛運(yùn)用。在數(shù)值模擬方面,尤其對泥石流主要致災(zāi)區(qū)——堆積區(qū)的研究較多。其中,日本學(xué)者石川芳治等(1991)通過水工模型試驗(yàn)用數(shù)值模擬研究了泥石流的堆積泛濫范圍,羅元華(1998)考慮侵蝕因素的控制方程,進(jìn)行泥石流堆積數(shù)值模擬。還有其他眾多學(xué)者在流動學(xué)實(shí)驗(yàn)及流動學(xué)模型的基礎(chǔ)上,張鵬等(2014)、Wuetal.(2013)、Linetal.(2005)、Bertoloetal.(2005)運(yùn)用FLO-2D模型較好地模擬了小流域泥石流運(yùn)動淤積過程。Johnsonetal.(1970)提出利用賓漢體黏性流解釋泥石流運(yùn)動過程,建立泥石流運(yùn)動控制方程,計(jì)算出泥石流最大流速的分布。Marcel(2006)等采用有限差分方法對泥石流溝谷進(jìn)行數(shù)值模擬確定了泥石流的危險范圍,并在此基礎(chǔ)上進(jìn)行了危險性分區(qū)。唐川等(1993)、唐川(1994a,b,c)運(yùn)用剖開算子有限差分法對泥石流進(jìn)行計(jì)算模擬,得出泥深及流速的分布,對泥石流危險度進(jìn)行評價。

    在眾多學(xué)者研究的基礎(chǔ)上,本文使用FLO-2D軟件加上雨量資料、數(shù)值地形資料與相關(guān)參數(shù)模擬泥石流的流動情況,以2011年湖北黃石富水河地區(qū)青石溝所發(fā)生泥石流為數(shù)值模擬原型,配合現(xiàn)有的地質(zhì)調(diào)查資料,來模擬泥石流的黏性流動參數(shù)及屈服應(yīng)力變化情況,分析泥石流的流動及沉積特性。

    1研究區(qū)概況

    青石溝位于陽新沿鎮(zhèn),屬于富水河流域,流域面積5310km2,屬北亞熱帶氣候區(qū),干流長196km,河道總落差613m。流域內(nèi)5km以上支流有130條,其中一級支流30條,二級支流52條,地下水豐富,但地下水類型與分布規(guī)律呈現(xiàn)多樣性。研究區(qū)內(nèi)溝谷呈南北向展布,形態(tài)多呈“V”字形,地勢北高南低,相對高差198m,為泥石流形成的典型地形。根據(jù)氣象資料顯示,該地區(qū)降雨量豐富,年均降雨量1389.6mm。年均降雨日147個,夏季最多, 4~7月平均降雨量739.9mm,雨量多,強(qiáng)度大,常造成洪澇災(zāi)害,此氣象條件極有利于泥石流的形成。在2011年7月12日,當(dāng)日累計(jì)降雨量超過250mm,當(dāng)天傍晚6點(diǎn)32分左右爆發(fā)了泥石流,方量約為1.5×104m3,為小型泥石流,整個流徑約為556m,但并未造成人員及財(cái)產(chǎn)損失。

    根據(jù)現(xiàn)場調(diào)查,青石溝泥石流形成區(qū)內(nèi)為三面環(huán)山、一面出口的圓形寬闊地段,縱長約200m,周圍山坡陡峭多為30°~60°的陡坡,經(jīng)測算匯水面積約10000m2。陡坡坡體巖體破碎,植被覆蓋率較低; 斜坡常被沖溝切割強(qiáng)烈,形成較多不穩(wěn)定斜坡,崩塌與滑坡時有發(fā)生。加之該地區(qū)出露的巖層主要為震旦系、奧陶系灰?guī)r、砂巖與泥巖互層,巖層巖性較軟,易受風(fēng)化侵蝕,形成第四系崩坡積松散堆積物,構(gòu)成了泥石流物質(zhì)來源,為泥石流的爆發(fā)提供了豐富的物質(zhì)基礎(chǔ)。流通區(qū)縱長約為200m,多為狹窄而深切的峽谷或沖溝,谷壁陡峻而縱坡降大,平均坡度為25°,泥深為2~3m,由于所夾挾物質(zhì)多為松散堆積物,并未對原溝谷產(chǎn)生較強(qiáng)烈的下切作用。堆積區(qū)縱長約為156m,位于青石溝的溝口處,地形較為開闊、平緩,平均坡度為10°,堆積厚度為5~8m,形成面積為2000m2的堆積區(qū)。

    2FLO-2D數(shù)值模擬

    2.1概要

    FLO-2D是美國科羅拉多州立大學(xué)的O’Brienetal.(1993)用有限差分法來求解運(yùn)動控制的程序,它可以在工程設(shè)計(jì)、洪水災(zāi)害管理、泥石流等問題上得到很好運(yùn)用。為了求取流動深度及流量,我們采用的數(shù)值地形為規(guī)則網(wǎng)格單元,并且每個網(wǎng)格都給定固定單一的高程和粗糙系數(shù),然后輸入各項(xiàng)參數(shù)進(jìn)行計(jì)算,在FLO-2D程序中,有以下幾項(xiàng)假設(shè)及限制:

    (1)差分時間間隔內(nèi)為穩(wěn)定流,

    (2)水道斷面和粗糙度是規(guī)則的,

    (3)凈水壓力分布,

    (4)每一個網(wǎng)格有單一的高程及N值,

    (5)無法模擬水道的刷深現(xiàn)象,

    (6)無法模擬上游崩塌可能對泥石流造成的影響,

    (7)假定泥石流流通區(qū)為定溝床模式,故無法模擬溝道的侵蝕現(xiàn)象。

    2.2建模

    FLO-2D的連續(xù)方程(1)及運(yùn)動方程(2)、(3)如下所示:

    (1)

    (2)

    (3)

    上式中,式(1)控制泥石流或洪水以及質(zhì)量守恒,其中,t為時間;i為降雨強(qiáng)度;u為X方向的平均流速;v為Y方向的平均流速;h為流動深度。

    式(2)與式(3)為力平衡的動量方程,動力波模式為動量方程主體,擴(kuò)散波模式為省略式中最后3項(xiàng),運(yùn)動波模式則是省略式中的壓力梯度。FLO-2D提供動力波模式及擴(kuò)散波模式來模擬洪水與泥石流問題。式中,g為重力加速度;Sfx及Sfy為摩擦坡降;Sox及Soy為底床坡降。

    針對泥石流的流變方程,O’Brienetal.(1985)設(shè)計(jì)了一套關(guān)于泥石流的計(jì)算方程式,如式(4),式(5)所示:

    (4)

    (5)

    上式中,τy為屈服應(yīng)力;η為黏性系數(shù);C為慣切剪應(yīng)力;K為層流阻滯系數(shù),取值范圍24~50000;ρ為泥石流密度;g為重力加速度;n為曼寧n值,參考FLO-2D手冊;h為流深;V為平均流速。

    又有式(6)~式(7)如下:

    (6)

    (7)

    式中,Cv為泥石流體積濃度,其中經(jīng)驗(yàn)系數(shù)α1、α2、β1、β2取值通過查閱相關(guān)FLO-2D手冊資料確定,其取值如表2所示。

    表1 不同沉積物濃度流動特性表

    表2 FLO-2D 輸入?yún)?shù)值

    FLO-2D中按體積濃度和流動特性分為4個等級,分別是滑坡,泥石流,泥流,洪水(表1)。研究中,根據(jù)對現(xiàn)場泥石流的調(diào)查情況,取表中的0.45作為體積濃度。

    2.3輸入變量

    本文以2011年富水河流域青石溝泥石流數(shù)據(jù)為基礎(chǔ),選取參數(shù)進(jìn)行計(jì)算。由于泥石流發(fā)生機(jī)制的復(fù)雜性,我們選取泥石流災(zāi)害地區(qū)流域面積比較大的青石溝流域,同時對該流域的資料進(jìn)行整理利用。根據(jù)實(shí)測情況,該區(qū)斷面是由形成區(qū)坡度為38°、流通區(qū)坡度15°、堆積區(qū)坡度6°所構(gòu)成的長556m,高198m斜面(圖1),在資料使用中,為了便于掌握泥石流堆積特性,距離為556m以后的視為平地。

    圖1 流域橫截面簡圖Fig. 1 Basin cross section

    流量按照實(shí)際的流域面積和降雨強(qiáng)度來計(jì)算,泥石流發(fā)生時當(dāng)日累計(jì)雨量超過250mm,降雨強(qiáng)度105mm·h-1。通過計(jì)算得到,流域最大流量為1.89m3·s-1。

    為了獲得移動特性和堆積特性隨著黏性系數(shù)和屈服應(yīng)力的變化關(guān)系明顯的規(guī)律,研究中各取標(biāo)準(zhǔn)黏性系數(shù)和屈服應(yīng)力值的10倍、100倍的差異值來進(jìn)行模擬分析。表2為本次模擬所涉及的主要參數(shù)取值,各參數(shù)的取值采取現(xiàn)場原位試驗(yàn)、室內(nèi)試驗(yàn)與經(jīng)驗(yàn)參數(shù)相綜合的辦法。并且,為了分析隨著黏性系數(shù)及屈服應(yīng)力的變化泥石流流動和沉積變化特點(diǎn),黏性系數(shù)依次為η0、10η0、100η0,屈服應(yīng)力依次為τy0、10τy0、100τy0,流入流量,體積濃度、運(yùn)移時間等其他條件保持不變。

    3分析結(jié)果

    3.1泥石流移動及沉積特性與黏性系數(shù)的關(guān)系

    從圖2 中可以看出,泥石流移動速度更依賴于斜坡坡角。同時,隨著黏性系數(shù)的增加,移動速度呈減小的趨勢,當(dāng)η0較小時,不同斷面傾斜角的移動速度呈現(xiàn)較大的差異,而隨著η0的增大,不同斷面傾斜角的移動速度差異慢慢減小。

    圖2 泥石流流速與黏性系數(shù)關(guān)系圖Fig. 2 Velocity of debris flow according to viscosity

    圖3為泥石流流深與黏性系數(shù)變化關(guān)系圖,泥石流的流深并沒有隨著傾斜角度的不同而呈現(xiàn)較大的差異。黏性系數(shù)越大,粒子間的凝聚力增大,流速越來越小,但是流深卻不斷變大。

    圖3 泥石流流深與黏性系數(shù)關(guān)系圖Fig. 3 Flow depth of debris flow according to viscosity

    圖4為泥石流流動范圍隨著黏性系數(shù)變化結(jié)果圖,其流動范圍結(jié)果圖是根據(jù)泥石流發(fā)生當(dāng)日,根據(jù)降雨情況及其他資料,通過數(shù)值分析得到的。分析表明黏性系數(shù)越大,流動距離和擴(kuò)散范圍則越小,當(dāng)黏性系數(shù)為100η0時,在堆積區(qū)坡角為6°的坡面中,其活動范圍在沒有達(dá)到平地的時候已經(jīng)停止,這是由于黏性增加,泥石流粒子間的凝聚力增加,泥石流的流動阻力增大。

    表3為FLO-2D模擬得出的泥石流沉積距離、

    平均沉積寬度、最大沉積寬度和堆積面積變化情況。隨著黏性系數(shù)的變大,泥石流沉積距離、沉積寬度,沉積面積逐漸減小。隨著黏性系數(shù)的增加,移動速度、移動距離、擴(kuò)散范圍所表現(xiàn)出來的反比例現(xiàn)象,是由于泥石流本身移動特性和沉積特性所影響的。

    表3 泥石流運(yùn)動特征與黏性系數(shù)關(guān)系表

    圖4 不同黏度的泥石流擴(kuò)散范圍Fig. 4 Spreading range of debris flow according to viscositya. η0; b. 10η0; c. 100η0

    隨著泥石流的流動,為了模擬分析其危險范圍,需要對泥石流的沖擊力進(jìn)行分析。FLO-2D中所用到的泥石流沖擊力計(jì)算公式如(8)~式(10),它受到流體密度,沖擊角度等多種因素的綜合影響。

    (8)

    k=1.261eCw

    (9)

    F=Pih

    (10)

    式中,Pi為沖擊壓力,通過流體密度ρf、系數(shù)k、網(wǎng)格中最大移動速度V來計(jì)算。單位面積沖擊力F為Pi所對應(yīng)的網(wǎng)格最大流深h的乘積。即沖擊力與最大移動速度的平方和流深呈正比例關(guān)系。

    圖5是利用不同黏性系數(shù)時堆積部的沖擊力和表3數(shù)據(jù)所繪制的,結(jié)合圖4 的計(jì)算結(jié)果可知,黏性系數(shù)越大,沖擊力則減小,這是由于移動速度快速減小的原因所導(dǎo)致。

    圖5 沖擊力、沉積區(qū)與黏性系數(shù)關(guān)系圖Fig. 5 Impact force and depositional area according to viscosity

    3.2泥石流移動及堆積特性與屈服應(yīng)力的關(guān)系

    圖6為形成區(qū)、流通區(qū)和堆積區(qū)中泥石流流速隨屈服應(yīng)力變化關(guān)系圖。盡管屈服應(yīng)力越大,但是移動速度的變化量并不是很明顯,而在屈服應(yīng)力相同的情況下,傾角越大,其流速也越大。所以,相對于屈服應(yīng)力而言,傾角對流速的變化影響更大。

    圖6 泥石流流速與屈服應(yīng)力關(guān)系圖Fig. 6 Velocity of debris flow according to yield stress

    圖7 泥石流流深與屈服應(yīng)力關(guān)系圖Fig. 7 Flow depth of debris flow according to yield stress

    圖9 不同屈服應(yīng)力下泥石流的擴(kuò)散范圍Fig. 9 Spreading range of debris flow according to yield stressa. τy0; b. 10τy0; c. 100τy0

    圖7為不同屈服應(yīng)力狀態(tài)下,流深的變化圖,隨著屈服應(yīng)力的增大,流深也變大。特別是在堆積區(qū)傾角為6°時的斜面內(nèi),隨著屈服應(yīng)力的增大,流深變大的尤為明顯,但是通過與圖3 對比可以得知,流深的變化更依賴于黏性系數(shù)的變化。

    圖8為堆積區(qū)不同屈服應(yīng)力下,沖擊力和堆積面積的變化圖。雖然沖擊力與黏性系數(shù)呈反比例關(guān)系,但是,隨著屈服應(yīng)力的增大,沖擊力也隨著增大。

    圖8 沖擊力、沉積面積與屈服應(yīng)力關(guān)系圖Fig. 8 Impact force and depositional area according to yield stress

    圖9為FLO-2D模擬下,不同屈服應(yīng)力狀態(tài),泥石流移動范圍變化結(jié)果圖,圖9 中,隨著屈服應(yīng)力的增大,傾斜面的運(yùn)動并沒有出現(xiàn)太大的變化,但是堆積區(qū)的擴(kuò)散范圍變小。屈服應(yīng)力愈大,發(fā)生初始流動的剪切應(yīng)力也越大,這使得泥石流在流動過程中會更快停止運(yùn)動。隨著屈服應(yīng)力的變化,泥石流的沉積特性呈現(xiàn)更大的差異性。由此可知,在影響泥石流發(fā)生和停止的因素中,與屈服應(yīng)力對泥石流流動特性影響相比,屈服應(yīng)力對堆積特性的影響更大。表4是通過圖9 得到堆積區(qū)移動特性表,隨著屈服應(yīng)力的變大,堆積距離、堆積寬度、堆積面積不斷減小。

    表4 泥石流運(yùn)動特征與屈服應(yīng)力關(guān)系表

    4結(jié)論

    本論文中,利用所收集到的湖北省黃石富水河流域青石溝泥石流的數(shù)據(jù),借助FLO-2D數(shù)值模擬分析程序,對泥石流的移動特性進(jìn)行模擬,以此分析泥石流流動及堆積特性與黏性系數(shù)和屈服應(yīng)力的關(guān)系,所獲得的結(jié)論如下:

    (1)隨著黏性系數(shù)的增大,泥石流流速呈非線性減少,流經(jīng)斜面坡角越大,其變化越大。即:流速同時受黏性系數(shù)與流經(jīng)斜面坡角的影響。斜面坡角對流深的影響不大,但是隨著黏性系數(shù)的增大,流深也隨著增加。黏性系數(shù)對泥石流的移動和堆積特性有比較大的影響,與泥石流流速、移動距離、擴(kuò)散范圍呈反比例關(guān)系。

    (2)屈服應(yīng)力作為影響泥石流發(fā)生及停止的因素,隨著屈服應(yīng)力的增大,流深也隨著增大,在泥石流移動特性中,屈服應(yīng)力對泥石流流速的影響并不明顯。隨著屈服應(yīng)力的增加,堆積距離,堆積寬度,堆積面積則減小。由此可知,相比屈服應(yīng)力對移動特性的影響,屈服應(yīng)力對堆積特性的影響更大。

    (3)隨著黏性系數(shù)的增加,泥石流沖擊力減小,而隨著屈服應(yīng)力的增加,泥石流的沖擊力則增大。

    參考文獻(xiàn)

    BertoloLP,WieczorekGF. 2005.CalibrationofnumericalmodelsforsmalldebrisflowsinYosemiteValley,California,USA[J].NaturalHazardsandEarthSystemSciences,5(6): 993~1001.

    HuelimannM,CoponsR,AltimirJ. 2006.DetaileddebrisflowhazardassessmentinAndora:AMultidisciplinarApproach[J].Geomorphology,78(3-4): 359~372.

    ImranJ,HarffP,ParkerG. 2001,Anumericalmodelofsubmarinedebrisflowwithgraphicaluserinterface,Computers&Geosciences,27(6): 717~729.

    Iverson,RM,RiedME,LaHusenRG. 1997.Debris-flowmobilizationfromlandslides[J].AnnualReviewofEarthandPlanetarySciences,25: 85~138.

    JohnsonAM,RahnPH. 1970.Mobilizationofdebrisflows[J].ZeitschriftfurGeomorphologie,3(9): 168~186.

    LinML,WangKL,HuangJJ. 2005.Debrisflowrunoffsimulationandverification-casestudyofChen-You-Lanwatershed[J].NaturalHazardsandEarthSystemSciences5(3): 439~445.

    LuoYH. 1998.Methodologicalresearchonnumericalsimulationanddisasterriskassessmentofdebrisflowdepositionfan[D].Wuhan:ChinaUniversityofGeosciences.

    O’BrienJS,JulienPY,F(xiàn)ullertonWT. 1993.Two-dimensionalwaterfloodandmudfloodsimulation.JournalofHydraulicEngineeringASCE,119(2): 244~260.

    O’BrienJS,JulienPY. 1985.Physicalpropertiesandmechanicsofhyper-concentratedsedimentflows,Proceedingsofthespecialtyconferenceondelineationoflandslide,F(xiàn)lashfloodanddebrisflowhazardinUtah[M].UtahStateUniversity:Utah.260~279.

    O’Brien,JS,Julien,PY. 1988. “LaboratoryanalysisofMudflowproperties”.[J]JournalofHydraulicEngineering,114(8): 877~887.

    TangC,LiuXL. 1993.Dynamicsimulationandpredictionmodelofdebrisflow[J].JournalofSoilandWaterConservation,5(3): 37~40.

    TangC.1994a.Numericalsimulationofdebrisflowinundationonthealluvialfansanditspredictionmodeloftheriskareas[J].JournalofsoilandWaterconservation,8(1): 45~50.

    TangC.1994b.Anumericalsimulationstudyontheriskassessmentofdebrisflowaccumulationfan[J].JournalofCatastrophology,9(4): 7~13.

    TangC.1994c.Discussiononnumericalsimulationmethodofplanetwo-dimensionaldebrisflow[J].HydrogeologyandEnginee-ringGeology,21(5): 9~12.

    WuYH,LiuKF,ChenYC, 2013.ComparisonbetweenFLO-2DandDebris-2Dontheapplicationofassessmentofgranulardebrisflowhazardswithcasestudy.[J].JournalofMountainScience,10(2): 293~304.

    ZhangP,MaJZ,ShuHP,WangG. 2014.NumericalsimulationoferosionanddepositiondebrisflowbasedonFLO-2DModel[J].JournalofLanzhouUniversity(NaturalSciences),50(3): 363~368.

    羅元華. 1998. 泥石流堆積數(shù)值模擬及泥石流災(zāi)害風(fēng)險評估方法研究[D]. 武漢:中國地質(zhì)大學(xué).

    石川芳治,水山高久,井戶清尾. 1991. 堆積扇上泥石流堆積泛濫機(jī)理[M]. 泥石流及洪水災(zāi)害防御國際學(xué)術(shù)討論會文集. 成都: 27~31.

    唐川,劉希林. 1993. 泥石流動力堆積模擬和危險范圍預(yù)測模型[J]. 水土保持學(xué)報,5(3): 37~40.

    唐川. 1994. 泥石流堆積泛濫過程的數(shù)值模擬及其危險范圍預(yù)測模型的研究[J]. 水土保持學(xué)報,8(1): 45~50.

    唐川. 1994. 泥石流堆積扇危險度分區(qū)評價的數(shù)值模擬研究[J]. 災(zāi)害學(xué),9(4): 7~13.

    唐川. 1994. 平面二維泥石流數(shù)值模擬方法的探討[J]. 水文地質(zhì)工程地質(zhì).21(5): 9~12.

    張鵬,馬金珠,舒和平,等, 2014. 基于FLO-2D模型的泥石流運(yùn)動沖淤數(shù)值模擬[J]. 蘭州大學(xué)學(xué)報(自然科學(xué)版).50(3): 363~368.

    RESEARCHONTHEINFLUENCEFACTORSOFFLOWANDDEPOSITIONOFDEBRISFLOWBASEDONTHEFLO-2DSIMULATION

    LIANGHongxiSHANGMinXUXin

    (College of Civil Engineering & Architecture,China Three Gorges University,Yichang 443002)

    AbstractDebris flow is different from the common landslide. It is a kind of temporary water with sand and stones occurred in a mountainous area. As researches shown, the debris flow is characterized with the nature of Bingham Body and moving jams. Many scholars had made many researches of numerical simulation method, which based on rheology experiment and model. And compared the scope, movement speed and movement time of the debris flow with the actual situation, the numerical simulation programs were developed for the forecast of debris flow. The FLO-2D is one of the software. As engineering practice proved, movement and sedimentary characteristics of debris flow is very important for further study based on the flow characteristics of research on debris flow substance composition. In order to analyze the relationship between debris flow movement-accumulation characteristics and viscosity-yield stress, we simulated the moving characteristics of debris flow by FLO-2D model. Through the analysis, we conclude that debris flow velocity is reduced with the increase of the viscosity, and the reduction is nonlinear. With the increase of yield stress, that is one of the impact factors about debris flow occurrence and stop, flow depth increases. But among the debris flow movement characteristics, the effect between yield stress and debris flow velocity is not obvious. With the increase of the viscosity, the impact force of debris flow decreases. Meanwhile, with the increase of the yield stress, the impact force of debris flow is increased.

    Key wordsNumerical simulation, FLO-2D,Viscosity, Yield stress

    DOI:10.13544/j.cnki.jeg.2016.02.008

    * 收稿日期:2015-01-30; 收到修改稿日期: 2015-05-05.

    基金項(xiàng)目:梁鴻熙(1991-),男,碩士生,土木與建筑學(xué)院地質(zhì)工程專業(yè). Email: 421486282@qq.com.

    第一作者簡介:尚敏(1977-),男,博士,副教授,主要從事地質(zhì)災(zāi)害機(jī)理與防治研究. Email:summing@126.com

    中圖分類號:P642.23

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

    猜你喜歡
    屈服應(yīng)力黏性泥石流
    潤滑劑對磁流變液屈服應(yīng)力的影響
    復(fù)雜流體的屈服應(yīng)力及其測定與應(yīng)用
    中國制筆(2020年2期)2020-07-03 09:20:52
    富硒產(chǎn)業(yè)需要強(qiáng)化“黏性”——安康能否玩轉(zhuǎn)“硒+”
    如何運(yùn)用播音主持技巧增強(qiáng)受眾黏性
    傳媒評論(2019年4期)2019-07-13 05:49:28
    泥石流
    雜文月刊(2018年21期)2019-01-05 05:55:28
    玩油灰黏性物成網(wǎng)紅
    華人時刊(2017年17期)2017-11-09 03:12:03
    “民謠泥石流”花粥:唱出自己
    海峽姐妹(2017年6期)2017-06-24 09:37:36
    泥石流
    鈣基潤滑脂替代鋰基潤滑脂可行性研究
    基層農(nóng)行提高客戶黏性淺析
    日本撒尿小便嘘嘘汇集6| 黄色a级毛片大全视频| 脱女人内裤的视频| 免费一级毛片在线播放高清视频| 国产亚洲欧美精品永久| 成年免费大片在线观看| 老司机福利观看| 久久精品国产99精品国产亚洲性色| 日本一区二区免费在线视频| 亚洲欧美激情综合另类| www国产在线视频色| 欧洲精品卡2卡3卡4卡5卡区| 久久热在线av| 午夜免费鲁丝| 欧美最黄视频在线播放免费| 国产精品野战在线观看| av电影中文网址| 777久久人妻少妇嫩草av网站| 黑人巨大精品欧美一区二区mp4| 亚洲一区二区三区色噜噜| 国产视频一区二区在线看| 麻豆久久精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 欧美日韩中文字幕国产精品一区二区三区| 国产极品粉嫩免费观看在线| 国内精品久久久久久久电影| 亚洲成av人片免费观看| 母亲3免费完整高清在线观看| 午夜免费激情av| 日本在线视频免费播放| 欧美精品啪啪一区二区三区| av在线播放免费不卡| 亚洲性夜色夜夜综合| 亚洲人成电影免费在线| 国产激情欧美一区二区| 男人舔女人的私密视频| 欧美性猛交黑人性爽| 1024手机看黄色片| 高潮久久久久久久久久久不卡| 日本成人三级电影网站| 黄色女人牲交| 99riav亚洲国产免费| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区三| 99国产精品一区二区蜜桃av| www.精华液| 欧美精品啪啪一区二区三区| 久久婷婷人人爽人人干人人爱| 狠狠狠狠99中文字幕| www国产在线视频色| 美国免费a级毛片| 国产av在哪里看| 国产精品爽爽va在线观看网站 | 级片在线观看| 国产午夜福利久久久久久| 国产私拍福利视频在线观看| 国产亚洲av嫩草精品影院| 黄色丝袜av网址大全| 亚洲欧美精品综合一区二区三区| 久久精品国产清高在天天线| 午夜精品在线福利| 国产激情欧美一区二区| 91成人精品电影| 一个人免费在线观看的高清视频| 日韩大码丰满熟妇| 99久久无色码亚洲精品果冻| 欧美日韩乱码在线| av视频在线观看入口| 欧美黄色淫秽网站| 亚洲av片天天在线观看| 在线观看66精品国产| 亚洲成人精品中文字幕电影| 深夜精品福利| 免费看十八禁软件| 国产成人一区二区三区免费视频网站| 法律面前人人平等表现在哪些方面| 国产1区2区3区精品| 欧美不卡视频在线免费观看 | 色婷婷久久久亚洲欧美| av免费在线观看网站| 亚洲精品在线观看二区| 久久婷婷人人爽人人干人人爱| 国产精品av久久久久免费| 免费女性裸体啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 制服诱惑二区| 亚洲一区高清亚洲精品| 国产久久久一区二区三区| 国产免费男女视频| 欧美日韩福利视频一区二区| 他把我摸到了高潮在线观看| 夜夜看夜夜爽夜夜摸| 亚洲avbb在线观看| 亚洲一码二码三码区别大吗| 99精品欧美一区二区三区四区| 精品人妻1区二区| netflix在线观看网站| 黄色丝袜av网址大全| 欧美日韩一级在线毛片| 成在线人永久免费视频| 啦啦啦 在线观看视频| 午夜福利在线观看吧| 午夜福利视频1000在线观看| 国内精品久久久久久久电影| 可以在线观看的亚洲视频| 观看免费一级毛片| 日本熟妇午夜| 成人av一区二区三区在线看| 18禁黄网站禁片免费观看直播| 在线观看一区二区三区| 欧美成人性av电影在线观看| 精品欧美国产一区二区三| www国产在线视频色| 日韩欧美在线二视频| 超碰成人久久| 熟女电影av网| 在线视频色国产色| 亚洲精品美女久久久久99蜜臀| 亚洲中文av在线| 久久香蕉激情| 国产片内射在线| 黄色成人免费大全| 午夜福利在线在线| 精品国产美女av久久久久小说| 国产三级黄色录像| 丝袜人妻中文字幕| 在线观看日韩欧美| 老司机福利观看| 国产高清videossex| 高清在线国产一区| 精品高清国产在线一区| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 成人国产一区最新在线观看| 亚洲激情在线av| 精品久久久久久久毛片微露脸| 91国产中文字幕| 欧美黄色淫秽网站| 可以在线观看毛片的网站| 国语自产精品视频在线第100页| 老汉色∧v一级毛片| 韩国精品一区二区三区| 女性生殖器流出的白浆| 国产精华一区二区三区| 99热6这里只有精品| 成人亚洲精品av一区二区| 好看av亚洲va欧美ⅴa在| 亚洲专区国产一区二区| av在线播放免费不卡| 亚洲九九香蕉| 一级作爱视频免费观看| 色综合欧美亚洲国产小说| 国产成人av激情在线播放| 亚洲国产欧美日韩在线播放| 国内揄拍国产精品人妻在线 | 淫秽高清视频在线观看| 亚洲国产精品成人综合色| 久久久久九九精品影院| 十分钟在线观看高清视频www| 女警被强在线播放| 成年人黄色毛片网站| 成人手机av| 国产精品 国内视频| 免费观看人在逋| 国产精品日韩av在线免费观看| 亚洲男人天堂网一区| 色综合亚洲欧美另类图片| 美女午夜性视频免费| 成人三级黄色视频| 婷婷精品国产亚洲av在线| 白带黄色成豆腐渣| 色播在线永久视频| 搡老熟女国产l中国老女人| 黄色丝袜av网址大全| 国内精品久久久久精免费| 精品国内亚洲2022精品成人| 一级毛片女人18水好多| 怎么达到女性高潮| 久久精品91蜜桃| 两个人免费观看高清视频| 1024视频免费在线观看| 日韩欧美国产在线观看| 动漫黄色视频在线观看| 久久久久久久久久黄片| 视频在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久| 国产视频一区二区在线看| 国产av一区二区精品久久| 亚洲精品久久国产高清桃花| 亚洲欧洲精品一区二区精品久久久| 一二三四在线观看免费中文在| 一区二区日韩欧美中文字幕| 国产高清有码在线观看视频 | 日本 av在线| 精品卡一卡二卡四卡免费| 9191精品国产免费久久| 欧美日韩黄片免| 男女床上黄色一级片免费看| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 国产精品久久久av美女十八| 嫩草影院精品99| 桃色一区二区三区在线观看| 亚洲专区字幕在线| 国产av在哪里看| av视频在线观看入口| 午夜成年电影在线免费观看| 欧美三级亚洲精品| 亚洲精品美女久久av网站| 午夜激情福利司机影院| 欧美大码av| 亚洲五月色婷婷综合| 91麻豆精品激情在线观看国产| 久久亚洲真实| 国产v大片淫在线免费观看| 香蕉久久夜色| 波多野结衣高清无吗| a级毛片a级免费在线| 一级作爱视频免费观看| 精品一区二区三区av网在线观看| 亚洲 欧美 日韩 在线 免费| 美女大奶头视频| 成人亚洲精品一区在线观看| 国产真实乱freesex| 高清在线国产一区| 欧美不卡视频在线免费观看 | 亚洲片人在线观看| 亚洲av第一区精品v没综合| 成年免费大片在线观看| 黑人操中国人逼视频| 波多野结衣高清无吗| 成人亚洲精品一区在线观看| 国产国语露脸激情在线看| 可以在线观看毛片的网站| 少妇粗大呻吟视频| 久久久国产精品麻豆| 非洲黑人性xxxx精品又粗又长| 两性午夜刺激爽爽歪歪视频在线观看 | 91大片在线观看| 此物有八面人人有两片| 1024手机看黄色片| 国产亚洲精品久久久久久毛片| 欧美性猛交黑人性爽| 级片在线观看| 国内揄拍国产精品人妻在线 | 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 亚洲av片天天在线观看| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 日本五十路高清| 桃红色精品国产亚洲av| 日韩国内少妇激情av| 99热只有精品国产| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 午夜福利视频1000在线观看| 757午夜福利合集在线观看| 国产精品,欧美在线| 午夜老司机福利片| 久久国产精品影院| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 国内少妇人妻偷人精品xxx网站 | 久久久久久免费高清国产稀缺| 国产午夜精品久久久久久| 国产亚洲欧美98| 亚洲熟女毛片儿| 国产精品爽爽va在线观看网站 | 色婷婷久久久亚洲欧美| 欧美日韩中文字幕国产精品一区二区三区| 日本三级黄在线观看| 18禁观看日本| 禁无遮挡网站| 国产精品国产高清国产av| 免费在线观看完整版高清| 久久久国产成人免费| 亚洲国产精品合色在线| 久久人妻福利社区极品人妻图片| 国产人伦9x9x在线观看| 国产精品久久久av美女十八| 99久久国产精品久久久| 亚洲精品在线美女| 国产精品久久久久久人妻精品电影| xxx96com| 日本五十路高清| 国产成人影院久久av| 久久国产亚洲av麻豆专区| 国产高清激情床上av| 热99re8久久精品国产| 精品福利观看| 午夜亚洲福利在线播放| 中文字幕久久专区| 一区福利在线观看| 精品一区二区三区四区五区乱码| 国产午夜福利久久久久久| 午夜亚洲福利在线播放| 欧美最黄视频在线播放免费| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 免费看十八禁软件| 久久久久国产精品人妻aⅴ院| 一二三四社区在线视频社区8| 成人三级做爰电影| av免费在线观看网站| 中文字幕人成人乱码亚洲影| 国产成人一区二区三区免费视频网站| 免费无遮挡裸体视频| 黄色视频不卡| 欧美另类亚洲清纯唯美| 不卡一级毛片| 一边摸一边抽搐一进一小说| 精品福利观看| 视频区欧美日本亚洲| 久久草成人影院| www.自偷自拍.com| 午夜精品久久久久久毛片777| 一二三四在线观看免费中文在| 精品免费久久久久久久清纯| 搞女人的毛片| 国产一区二区三区在线臀色熟女| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av香蕉五月| 日韩有码中文字幕| 91九色精品人成在线观看| 黄色女人牲交| 免费在线观看成人毛片| 国产精品一区二区精品视频观看| 免费观看人在逋| 亚洲av成人一区二区三| 成人亚洲精品一区在线观看| 人人妻人人看人人澡| 亚洲成av片中文字幕在线观看| 正在播放国产对白刺激| 他把我摸到了高潮在线观看| 级片在线观看| 村上凉子中文字幕在线| 最近最新免费中文字幕在线| 久久香蕉激情| 99riav亚洲国产免费| 欧美中文日本在线观看视频| 欧美日韩精品网址| www.999成人在线观看| 十分钟在线观看高清视频www| 久久久久国产精品人妻aⅴ院| 亚洲第一青青草原| 国产色视频综合| 亚洲午夜理论影院| 日韩欧美一区视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲专区字幕在线| 1024香蕉在线观看| 中文字幕精品亚洲无线码一区 | ponron亚洲| 欧美日韩乱码在线| 国产aⅴ精品一区二区三区波| 男人舔奶头视频| 欧美黑人巨大hd| 啦啦啦韩国在线观看视频| 精华霜和精华液先用哪个| 欧美乱妇无乱码| 无限看片的www在线观看| 国产一区二区三区在线臀色熟女| 中文资源天堂在线| 黄频高清免费视频| 激情在线观看视频在线高清| 一个人观看的视频www高清免费观看 | cao死你这个sao货| 国产精品1区2区在线观看.| 国产成人精品久久二区二区91| 99国产精品一区二区三区| av在线天堂中文字幕| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 精品卡一卡二卡四卡免费| 淫秽高清视频在线观看| 国产伦一二天堂av在线观看| 免费女性裸体啪啪无遮挡网站| 搞女人的毛片| 国产精品久久久久久人妻精品电影| 亚洲午夜精品一区,二区,三区| 亚洲七黄色美女视频| 在线十欧美十亚洲十日本专区| а√天堂www在线а√下载| 欧美国产精品va在线观看不卡| 日本免费a在线| 国产伦一二天堂av在线观看| 成人午夜高清在线视频 | 精品久久久久久久末码| 午夜两性在线视频| 一本一本综合久久| 12—13女人毛片做爰片一| 制服诱惑二区| 欧美三级亚洲精品| 欧美日韩福利视频一区二区| 人人妻,人人澡人人爽秒播| 97超级碰碰碰精品色视频在线观看| www.自偷自拍.com| 女人被狂操c到高潮| 国产野战对白在线观看| 人妻久久中文字幕网| 日韩 欧美 亚洲 中文字幕| 精品日产1卡2卡| 天天躁狠狠躁夜夜躁狠狠躁| 老司机福利观看| 欧美黑人精品巨大| 日韩欧美国产在线观看| 亚洲av成人av| 亚洲av片天天在线观看| 51午夜福利影视在线观看| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频| 婷婷精品国产亚洲av| 又大又爽又粗| 麻豆久久精品国产亚洲av| av福利片在线| 最好的美女福利视频网| 欧美另类亚洲清纯唯美| 亚洲第一青青草原| 国产精品野战在线观看| 久久人妻福利社区极品人妻图片| 欧美不卡视频在线免费观看 | 国产aⅴ精品一区二区三区波| 久久久久久久久久黄片| 国产1区2区3区精品| 黑丝袜美女国产一区| 免费在线观看视频国产中文字幕亚洲| 一个人观看的视频www高清免费观看 | 别揉我奶头~嗯~啊~动态视频| 亚洲 国产 在线| 很黄的视频免费| 欧美av亚洲av综合av国产av| 亚洲男人的天堂狠狠| 久久精品国产亚洲av香蕉五月| 观看免费一级毛片| 国产亚洲欧美在线一区二区| 禁无遮挡网站| 少妇粗大呻吟视频| 黄色片一级片一级黄色片| 午夜免费成人在线视频| 女人被狂操c到高潮| 在线观看一区二区三区| 久久人妻av系列| 免费高清视频大片| 午夜日韩欧美国产| 亚洲精品中文字幕一二三四区| 国产精品99久久99久久久不卡| 一本久久中文字幕| 伊人久久大香线蕉亚洲五| 首页视频小说图片口味搜索| 久久99热这里只有精品18| 黄色a级毛片大全视频| videosex国产| 久久久久国产精品人妻aⅴ院| 视频在线观看一区二区三区| 午夜视频精品福利| 桃色一区二区三区在线观看| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 天堂动漫精品| 精品卡一卡二卡四卡免费| 18禁美女被吸乳视频| 麻豆成人av在线观看| 欧美人与性动交α欧美精品济南到| 欧美又色又爽又黄视频| 午夜福利成人在线免费观看| 亚洲男人天堂网一区| 国产1区2区3区精品| 午夜亚洲福利在线播放| 亚洲人成网站在线播放欧美日韩| 欧美日韩一级在线毛片| 欧美av亚洲av综合av国产av| 欧美大码av| 亚洲中文av在线| 搞女人的毛片| 桃红色精品国产亚洲av| 国产精品影院久久| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出| 国产精品亚洲美女久久久| 亚洲精华国产精华精| 老司机靠b影院| 99久久精品国产亚洲精品| 巨乳人妻的诱惑在线观看| 国产又黄又爽又无遮挡在线| 90打野战视频偷拍视频| 国产黄片美女视频| 中文字幕人妻熟女乱码| 丁香欧美五月| 国产aⅴ精品一区二区三区波| 琪琪午夜伦伦电影理论片6080| 日本熟妇午夜| 亚洲精品久久国产高清桃花| 久久精品人妻少妇| 亚洲国产中文字幕在线视频| 亚洲精品美女久久久久99蜜臀| 91九色精品人成在线观看| 亚洲五月天丁香| 999精品在线视频| av片东京热男人的天堂| 国产主播在线观看一区二区| 人妻丰满熟妇av一区二区三区| 嫩草影院精品99| 亚洲一区中文字幕在线| 一本一本综合久久| 在线视频色国产色| 久久久久久久久免费视频了| 日韩视频一区二区在线观看| 一本综合久久免费| 国产成人影院久久av| 精品电影一区二区在线| 色综合欧美亚洲国产小说| 日本精品一区二区三区蜜桃| 亚洲久久久国产精品| 亚洲av成人一区二区三| 午夜福利在线观看吧| 大型黄色视频在线免费观看| 人妻久久中文字幕网| 亚洲精品美女久久久久99蜜臀| 成人永久免费在线观看视频| 99riav亚洲国产免费| 人人妻,人人澡人人爽秒播| 色尼玛亚洲综合影院| 国产精品日韩av在线免费观看| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| 国产三级在线视频| 日韩欧美在线二视频| 国产三级黄色录像| 欧美乱色亚洲激情| 成人手机av| 最好的美女福利视频网| 国产精品永久免费网站| 亚洲人成77777在线视频| 在线观看舔阴道视频| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 级片在线观看| 黄色女人牲交| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| xxx96com| 国产一卡二卡三卡精品| 黄片播放在线免费| 黄色视频不卡| 好看av亚洲va欧美ⅴa在| 久久人妻福利社区极品人妻图片| av电影中文网址| av视频在线观看入口| 亚洲专区国产一区二区| 久久午夜亚洲精品久久| 亚洲熟女毛片儿| bbb黄色大片| 成人三级黄色视频| 亚洲国产精品久久男人天堂| 99热只有精品国产| 精品久久久久久久人妻蜜臀av| videosex国产| 久久精品国产综合久久久| 国产一区二区三区在线臀色熟女| 国产免费男女视频| 欧美zozozo另类| 亚洲三区欧美一区| 桃红色精品国产亚洲av| 18禁观看日本| 欧美日韩福利视频一区二区| 中文亚洲av片在线观看爽| 狠狠狠狠99中文字幕| 久久精品亚洲精品国产色婷小说| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 1024手机看黄色片| 午夜影院日韩av| 欧美另类亚洲清纯唯美| 少妇被粗大的猛进出69影院| 欧美一级a爱片免费观看看 | 神马国产精品三级电影在线观看 | 国产99久久九九免费精品| a级毛片a级免费在线| 久久久久国产精品人妻aⅴ院| 一区二区三区精品91| 国产精品影院久久| 久久伊人香网站| 天堂√8在线中文| 韩国av一区二区三区四区| 啦啦啦 在线观看视频| 女人被狂操c到高潮| 国产精品久久久人人做人人爽| 黄色a级毛片大全视频| 1024香蕉在线观看| 亚洲国产欧美网| 一个人观看的视频www高清免费观看 | 久久99热这里只有精品18| 午夜免费鲁丝| 免费看美女性在线毛片视频| 亚洲专区国产一区二区| 大香蕉久久成人网| 国产精品影院久久| 亚洲第一电影网av| 久久久国产精品麻豆| 色老头精品视频在线观看| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 黄频高清免费视频| 日本成人三级电影网站| 久久香蕉激情| 美女大奶头视频| 成人亚洲精品av一区二区| 国产成人系列免费观看| 成年免费大片在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 久久中文字幕一级| 一级黄色大片毛片|