• <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)行提高客戶黏性淺析
    人人妻人人澡欧美一区二区| 日韩强制内射视频| 亚洲av第一区精品v没综合| 91在线精品国自产拍蜜月| av视频在线观看入口| 亚洲五月天丁香| 春色校园在线视频观看| 免费在线观看成人毛片| 美女高潮的动态| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲网站| 午夜福利成人在线免费观看| 又爽又黄a免费视频| 久久久久免费精品人妻一区二区| av.在线天堂| 天天躁夜夜躁狠狠久久av| 亚洲国产精品sss在线观看| 日韩欧美一区二区三区在线观看| 床上黄色一级片| 日韩强制内射视频| 午夜老司机福利剧场| 亚洲最大成人手机在线| 久久精品久久久久久久性| 久久精品国产鲁丝片午夜精品| 国产精品电影一区二区三区| 麻豆成人av视频| 亚洲成人av在线免费| 日韩人妻高清精品专区| 精华霜和精华液先用哪个| 99久久中文字幕三级久久日本| 中国国产av一级| 亚洲婷婷狠狠爱综合网| 人妻系列 视频| 久久精品国产鲁丝片午夜精品| 亚洲人成网站在线观看播放| 在线a可以看的网站| 日本免费一区二区三区高清不卡| 热99在线观看视频| 有码 亚洲区| 精品免费久久久久久久清纯| 在线免费十八禁| 中出人妻视频一区二区| 啦啦啦韩国在线观看视频| 观看免费一级毛片| 亚洲在久久综合| 99视频精品全部免费 在线| 此物有八面人人有两片| 亚洲18禁久久av| 中国国产av一级| 日本欧美国产在线视频| 久99久视频精品免费| 精品99又大又爽又粗少妇毛片| 免费观看精品视频网站| 国产精品1区2区在线观看.| 亚洲成人av在线免费| 波多野结衣高清无吗| 菩萨蛮人人尽说江南好唐韦庄 | 好男人在线观看高清免费视频| 久久精品夜色国产| 日韩欧美国产在线观看| 国产伦一二天堂av在线观看| 日韩人妻高清精品专区| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久av| 日本爱情动作片www.在线观看| 国产老妇伦熟女老妇高清| 男人狂女人下面高潮的视频| 精品久久久久久久人妻蜜臀av| 欧美不卡视频在线免费观看| 国产精品福利在线免费观看| 国产一区二区在线av高清观看| 亚洲国产精品成人综合色| 国产精品一区www在线观看| 国产91av在线免费观看| 国产成人午夜福利电影在线观看| 国产老妇女一区| ponron亚洲| 亚洲成人中文字幕在线播放| 亚洲国产精品成人综合色| 午夜激情欧美在线| 中文在线观看免费www的网站| 欧美最黄视频在线播放免费| 黄色欧美视频在线观看| 婷婷精品国产亚洲av| 午夜福利在线在线| 天堂√8在线中文| 亚洲国产精品成人综合色| 欧美xxxx性猛交bbbb| 亚洲人成网站在线播| 国产亚洲精品av在线| 成人特级黄色片久久久久久久| 国产午夜精品久久久久久一区二区三区| 亚洲18禁久久av| 欧美zozozo另类| 最近2019中文字幕mv第一页| 午夜久久久久精精品| 欧美日韩国产亚洲二区| 国产爱豆传媒在线观看| 国产精品免费一区二区三区在线| 日本av手机在线免费观看| 国产片特级美女逼逼视频| 特大巨黑吊av在线直播| 蜜桃亚洲精品一区二区三区| 18+在线观看网站| 成人特级黄色片久久久久久久| 日本色播在线视频| 我要搜黄色片| 日韩一区二区三区影片| 国产在线男女| 人妻夜夜爽99麻豆av| 精品一区二区三区人妻视频| 久久久久久久久大av| 18禁在线播放成人免费| 非洲黑人性xxxx精品又粗又长| 国内久久婷婷六月综合欲色啪| 欧美高清成人免费视频www| 国产亚洲精品久久久久久毛片| 欧美日本亚洲视频在线播放| 亚洲精品色激情综合| 欧美激情在线99| 中文字幕精品亚洲无线码一区| 国产乱人偷精品视频| 欧美日本视频| 身体一侧抽搐| 亚洲一级一片aⅴ在线观看| 免费大片18禁| 久久久久网色| 国产乱人偷精品视频| 黄色配什么色好看| 久久99蜜桃精品久久| 久久久久久久久久黄片| 日韩高清综合在线| 欧美成人免费av一区二区三区| 日韩欧美在线乱码| 亚洲精品亚洲一区二区| 国产成年人精品一区二区| 激情 狠狠 欧美| 国产精品嫩草影院av在线观看| 亚洲图色成人| 婷婷亚洲欧美| 国产午夜福利久久久久久| 日韩精品有码人妻一区| 狠狠狠狠99中文字幕| 午夜爱爱视频在线播放| 成人一区二区视频在线观看| .国产精品久久| av免费观看日本| 久久久成人免费电影| 青青草视频在线视频观看| 久久久久网色| 国产日韩欧美在线精品| 老熟妇乱子伦视频在线观看| 少妇人妻一区二区三区视频| or卡值多少钱| 欧美色欧美亚洲另类二区| 精品久久久久久久久久免费视频| 看十八女毛片水多多多| 欧美日韩精品成人综合77777| 99九九线精品视频在线观看视频| 成人性生交大片免费视频hd| 免费观看精品视频网站| 亚洲人成网站在线播放欧美日韩| 草草在线视频免费看| 听说在线观看完整版免费高清| 日韩亚洲欧美综合| 成人综合一区亚洲| 欧美区成人在线视频| 男人舔女人下体高潮全视频| 最近中文字幕高清免费大全6| 欧美xxxx性猛交bbbb| 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡| 深夜精品福利| 插阴视频在线观看视频| 国产大屁股一区二区在线视频| 欧美激情久久久久久爽电影| 亚洲精品日韩av片在线观看| 久久精品影院6| 在线观看午夜福利视频| 国产亚洲av片在线观看秒播厂 | 中文字幕免费在线视频6| 99热这里只有是精品50| 午夜福利在线观看免费完整高清在 | 久久人人精品亚洲av| 丰满乱子伦码专区| 人妻系列 视频| 网址你懂的国产日韩在线| 在线观看一区二区三区| 国产亚洲精品久久久com| 国产色爽女视频免费观看| 少妇熟女aⅴ在线视频| 国产高清不卡午夜福利| 99热网站在线观看| 一本久久中文字幕| 最近最新中文字幕大全电影3| 亚洲人成网站在线播放欧美日韩| 免费电影在线观看免费观看| 最近中文字幕高清免费大全6| 欧美3d第一页| 久久久久久久久久成人| 婷婷色综合大香蕉| 毛片一级片免费看久久久久| 亚洲国产色片| 久久这里有精品视频免费| 国产女主播在线喷水免费视频网站 | 久久精品国产自在天天线| 欧美激情在线99| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品合色在线| 午夜a级毛片| 91久久精品国产一区二区成人| 久久久精品大字幕| 最新中文字幕久久久久| 色视频www国产| 欧美区成人在线视频| 中国国产av一级| 乱码一卡2卡4卡精品| av在线观看视频网站免费| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 成人国产麻豆网| 国产成人aa在线观看| 黄色一级大片看看| 精品久久国产蜜桃| 男女那种视频在线观看| 国产不卡一卡二| 成人午夜高清在线视频| 精品日产1卡2卡| 人妻制服诱惑在线中文字幕| 99久久中文字幕三级久久日本| 日韩亚洲欧美综合| 久久亚洲国产成人精品v| 2022亚洲国产成人精品| 中文字幕熟女人妻在线| 久久精品久久久久久久性| 国产成人一区二区在线| 亚洲精品日韩在线中文字幕 | 午夜久久久久精精品| 变态另类成人亚洲欧美熟女| 亚洲自偷自拍三级| 可以在线观看的亚洲视频| 亚洲不卡免费看| 伦理电影大哥的女人| 久久国产乱子免费精品| 成人午夜高清在线视频| 直男gayav资源| 欧美不卡视频在线免费观看| 丰满人妻一区二区三区视频av| 男人的好看免费观看在线视频| 免费看av在线观看网站| 亚洲天堂国产精品一区在线| 亚洲无线在线观看| 最近视频中文字幕2019在线8| 中国国产av一级| 毛片一级片免费看久久久久| 99久国产av精品国产电影| 可以在线观看的亚洲视频| 中文字幕制服av| 可以在线观看毛片的网站| 亚洲最大成人中文| 边亲边吃奶的免费视频| 日韩av不卡免费在线播放| 26uuu在线亚洲综合色| 又粗又硬又长又爽又黄的视频 | 三级国产精品欧美在线观看| 亚洲精品成人久久久久久| 亚洲av成人av| 亚洲国产欧美在线一区| 嫩草影院精品99| 国产精品嫩草影院av在线观看| 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| av在线播放精品| 国产精品一二三区在线看| www.av在线官网国产| 村上凉子中文字幕在线| 毛片一级片免费看久久久久| 亚洲图色成人| av.在线天堂| 亚洲国产欧美在线一区| 三级毛片av免费| 色5月婷婷丁香| 日韩欧美精品v在线| 国产av不卡久久| 亚洲精品国产成人久久av| 一级黄片播放器| 亚洲成人久久爱视频| 天堂网av新在线| 三级男女做爰猛烈吃奶摸视频| 一进一出抽搐gif免费好疼| 91麻豆精品激情在线观看国产| 日本av手机在线免费观看| 免费观看的影片在线观看| 亚洲av一区综合| 午夜福利视频1000在线观看| 桃色一区二区三区在线观看| 联通29元200g的流量卡| 少妇熟女aⅴ在线视频| 神马国产精品三级电影在线观看| 国产精品久久久久久亚洲av鲁大| 久久久久久国产a免费观看| 国产成人影院久久av| 精品人妻熟女av久视频| 内地一区二区视频在线| a级一级毛片免费在线观看| h日本视频在线播放| 国内精品美女久久久久久| 午夜a级毛片| 亚洲av成人av| 男女视频在线观看网站免费| 欧美色欧美亚洲另类二区| 一级黄片播放器| 观看免费一级毛片| 亚洲国产精品合色在线| 在线免费观看的www视频| 亚洲精品久久国产高清桃花| 国产成人午夜福利电影在线观看| 哪个播放器可以免费观看大片| 18禁在线播放成人免费| 级片在线观看| 精品人妻偷拍中文字幕| 99九九线精品视频在线观看视频| 观看免费一级毛片| 自拍偷自拍亚洲精品老妇| 亚洲五月天丁香| 日韩一区二区视频免费看| 亚洲国产精品sss在线观看| 一本久久精品| 国产伦理片在线播放av一区 | 日韩精品有码人妻一区| 老司机影院成人| 在线观看66精品国产| 国产高清三级在线| 啦啦啦韩国在线观看视频| 亚洲激情五月婷婷啪啪| 国内精品宾馆在线| 免费观看在线日韩| 男人和女人高潮做爰伦理| 日本成人三级电影网站| 嫩草影院入口| 黄色日韩在线| 午夜福利在线观看吧| 国产精品.久久久| 日韩成人av中文字幕在线观看| 国产亚洲精品av在线| 91久久精品国产一区二区三区| 啦啦啦观看免费观看视频高清| 97超视频在线观看视频| 十八禁国产超污无遮挡网站| 国产高清视频在线观看网站| 啦啦啦韩国在线观看视频| 五月玫瑰六月丁香| 99精品在免费线老司机午夜| 国产一区二区亚洲精品在线观看| av专区在线播放| 国产亚洲91精品色在线| 国产伦精品一区二区三区四那| 少妇的逼好多水| 99热精品在线国产| a级毛片免费高清观看在线播放| 日韩一区二区视频免费看| 变态另类丝袜制服| 国产午夜精品一二区理论片| 黄色视频,在线免费观看| 日韩在线高清观看一区二区三区| 国产视频首页在线观看| 国产伦一二天堂av在线观看| 国产高清有码在线观看视频| 黄色一级大片看看| 好男人视频免费观看在线| av天堂中文字幕网| 亚洲精品久久国产高清桃花| 国产极品精品免费视频能看的| 99久国产av精品国产电影| 国产极品精品免费视频能看的| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 老司机影院成人| 悠悠久久av| 成人性生交大片免费视频hd| 亚洲精品久久久久久婷婷小说 | 三级毛片av免费| 此物有八面人人有两片| 干丝袜人妻中文字幕| 精品人妻视频免费看| 中出人妻视频一区二区| 热99在线观看视频| 丝袜喷水一区| 亚洲内射少妇av| 免费一级毛片在线播放高清视频| 亚洲精品国产成人久久av| 在线国产一区二区在线| 一级二级三级毛片免费看| 国产午夜精品论理片| 国产黄色小视频在线观看| 黄色配什么色好看| 一边摸一边抽搐一进一小说| 日本黄色视频三级网站网址| 精品人妻熟女av久视频| 欧美日本亚洲视频在线播放| 亚洲av电影不卡..在线观看| 26uuu在线亚洲综合色| 国产久久久一区二区三区| 亚洲欧美成人综合另类久久久 | 成人特级黄色片久久久久久久| 国内精品宾馆在线| 国产高潮美女av| 在线免费十八禁| 亚洲精品亚洲一区二区| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 国产精品永久免费网站| 美女黄网站色视频| av专区在线播放| 久久久精品欧美日韩精品| 欧美+日韩+精品| 一夜夜www| 亚洲色图av天堂| 性色avwww在线观看| 永久网站在线| 99热6这里只有精品| av天堂在线播放| 美女高潮的动态| 看非洲黑人一级黄片| 国内久久婷婷六月综合欲色啪| 丰满人妻一区二区三区视频av| 欧美一级a爱片免费观看看| 色哟哟·www| 久久久久久久久中文| 乱人视频在线观看| 国产精品一区二区三区四区免费观看| 好男人视频免费观看在线| 成年女人看的毛片在线观看| 午夜激情欧美在线| 五月伊人婷婷丁香| 国产乱人偷精品视频| 亚洲18禁久久av| 波多野结衣巨乳人妻| 黄色视频,在线免费观看| 亚洲最大成人中文| 久久久久九九精品影院| 91午夜精品亚洲一区二区三区| 日本黄色片子视频| 最新中文字幕久久久久| 亚洲在线观看片| 91久久精品国产一区二区三区| 久久亚洲国产成人精品v| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 人妻少妇偷人精品九色| 久久国产乱子免费精品| 乱码一卡2卡4卡精品| 亚洲内射少妇av| 尾随美女入室| 少妇猛男粗大的猛烈进出视频 | 免费看av在线观看网站| 久久久成人免费电影| 久久精品影院6| av在线老鸭窝| av又黄又爽大尺度在线免费看 | 亚洲经典国产精华液单| 99在线视频只有这里精品首页| 国产av在哪里看| 精品久久久久久成人av| АⅤ资源中文在线天堂| 中国美白少妇内射xxxbb| 亚洲aⅴ乱码一区二区在线播放| 国产白丝娇喘喷水9色精品| 一本久久精品| 亚洲精品日韩在线中文字幕 | 亚洲熟妇中文字幕五十中出| 久久久久免费精品人妻一区二区| 亚洲高清免费不卡视频| 久久久久久九九精品二区国产| 免费观看的影片在线观看| 亚洲人成网站在线播放欧美日韩| 日产精品乱码卡一卡2卡三| 免费观看人在逋| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 亚洲国产精品久久男人天堂| 亚洲欧美日韩无卡精品| 精品无人区乱码1区二区| 亚洲精品亚洲一区二区| 国产真实乱freesex| 国内精品美女久久久久久| 亚洲精品国产av成人精品| 成人亚洲欧美一区二区av| 国产成人a∨麻豆精品| 黄色一级大片看看| 91在线精品国自产拍蜜月| 搞女人的毛片| 一级二级三级毛片免费看| 久久久久免费精品人妻一区二区| 乱人视频在线观看| 亚洲精华国产精华液的使用体验 | 亚洲久久久久久中文字幕| 色噜噜av男人的天堂激情| 亚洲av不卡在线观看| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 国产一区二区三区av在线 | 国产真实伦视频高清在线观看| 桃色一区二区三区在线观看| av在线观看视频网站免费| 人妻系列 视频| 亚洲欧美精品专区久久| 国产 一区 欧美 日韩| 久久久精品欧美日韩精品| 亚洲av二区三区四区| 久99久视频精品免费| 26uuu在线亚洲综合色| 中文字幕久久专区| 一个人免费在线观看电影| 国产高清不卡午夜福利| 精品无人区乱码1区二区| 久久精品国产亚洲av香蕉五月| 国产精品野战在线观看| 99久久无色码亚洲精品果冻| 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 99在线人妻在线中文字幕| 特级一级黄色大片| 免费大片18禁| 熟女电影av网| 欧美丝袜亚洲另类| 亚洲欧洲国产日韩| 亚洲丝袜综合中文字幕| 免费搜索国产男女视频| 大型黄色视频在线免费观看| 亚洲图色成人| 国产精品日韩av在线免费观看| 久久鲁丝午夜福利片| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 国产探花在线观看一区二区| 女人十人毛片免费观看3o分钟| 边亲边吃奶的免费视频| 成年av动漫网址| 国国产精品蜜臀av免费| 伦理电影大哥的女人| 国产成人aa在线观看| 婷婷六月久久综合丁香| 免费av观看视频| 性色avwww在线观看| 免费看a级黄色片| 久久久a久久爽久久v久久| 青青草视频在线视频观看| 午夜福利在线观看免费完整高清在 | 中文亚洲av片在线观看爽| 99热全是精品| 亚洲一级一片aⅴ在线观看| 特级一级黄色大片| 国产片特级美女逼逼视频| 在线观看一区二区三区| 少妇的逼水好多| 能在线免费看毛片的网站| 自拍偷自拍亚洲精品老妇| 国产在线精品亚洲第一网站| 综合色丁香网| 一区二区三区高清视频在线| 久久精品国产亚洲av天美| 中文字幕熟女人妻在线| 久久精品91蜜桃| 国产精品蜜桃在线观看 | 成年av动漫网址| 男女边吃奶边做爰视频| 三级男女做爰猛烈吃奶摸视频| 99热这里只有是精品在线观看| 黑人高潮一二区| 黄色配什么色好看| 最近的中文字幕免费完整| 国产蜜桃级精品一区二区三区| 在线国产一区二区在线| 国产精华一区二区三区| eeuss影院久久| 久久这里有精品视频免费| 一区二区三区免费毛片| 亚洲真实伦在线观看| 级片在线观看| 青春草亚洲视频在线观看| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 97在线视频观看| 插逼视频在线观看| 中出人妻视频一区二区| 国产一区二区在线av高清观看| 国产毛片a区久久久久| 99热6这里只有精品| 我的女老师完整版在线观看| 亚洲国产日韩欧美精品在线观看| 国模一区二区三区四区视频| 国产精品福利在线免费观看| 亚洲在线观看片| 狠狠狠狠99中文字幕| 日本免费一区二区三区高清不卡| 国产在视频线在精品| 青青草视频在线视频观看| 欧美精品一区二区大全| 亚洲av电影不卡..在线观看| 激情 狠狠 欧美| 天堂中文最新版在线下载 | 2022亚洲国产成人精品| 天堂网av新在线| 亚洲五月天丁香| 深夜a级毛片| 欧美变态另类bdsm刘玥| 禁无遮挡网站| 亚洲欧美成人综合另类久久久 | 亚洲精品日韩av片在线观看| 麻豆一二三区av精品| 国产亚洲av嫩草精品影院| 日韩精品有码人妻一区| 91狼人影院| 精品人妻熟女av久视频|