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

    基于DEM方法的滑坡堰塞壩幾何特征分析

    2019-10-21 09:17:21趙高文喬建平姜元俊孟華君王萌楊宗佶
    人民黃河 2019年5期
    關(guān)鍵詞:泥石流滑坡

    趙高文 喬建平 姜元俊 孟華君 王萌 楊宗佶

    摘要:基于離散元法( DEM)對滑坡堰塞壩的幾何特征進行了數(shù)值模擬,考慮溝道斷面形態(tài)、滑坡速度以及溝床坡度三種因素的影響,共27種組合。結(jié)果表明:溝道斷面形態(tài)、滑坡速度對堰塞壩橫向、縱向幾何形態(tài)都有影響,而溝床坡度只對堰塞壩縱向幾何特征產(chǎn)生影響;滑坡速度增大,堰塞壩縱剖面頂寬和底寬的比值減小,堰塞壩的有效厚度減小;溝床坡度增大,堰塞壩下游壩坡內(nèi)角減小,上游壩坡內(nèi)角增大,兩者的比值減小,堰塞壩壩坡傾角的變化規(guī)律與內(nèi)角的變化規(guī)律相反。

    關(guān)鍵詞:滑坡;堰塞壩;泥石流;漫頂破壞;離散元

    中圖分類號:P694;P642.22

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

    doi:10.3969/j.issn.1000- 1379.2019.05 .003

    滑坡堰塞壩屬于堰塞壩的一種,是滑坡造成溝道或河道堵塞而形成的天然壩,該類壩體穩(wěn)定性很差。Costa等[1]統(tǒng)計分析了73座滑坡堰塞壩(堵江型)破壞案例,發(fā)現(xiàn)其中27%在形成后一天內(nèi)破壞、50%在形成后10 d內(nèi)破壞,破壞的主要原因是漫頂侵蝕。柴賀軍等[2]在Costa等的研究基礎(chǔ)上,增加24個破壞案例后進一步指出.21%的滑坡堰塞壩在形成后一天內(nèi)破壞.48%在形成后10 d內(nèi)破壞,78%在形成后半年內(nèi)破壞,88%在形成后一年內(nèi)破壞?;卵呷麎稳舳氯拥佬纬纱笮脱呷赡軙蜎]上游農(nóng)田與民宅,若潰決形成洪水或泥石流,往往給下游帶來巨大損失。例如:1786年6月四川瀘定摩崗嶺滑坡堰塞湖潰決,潰決洪水不僅摧毀了下游一切基礎(chǔ)設(shè)施,還造成了近10萬人死亡[3]:1933年疊溪滑坡堰塞湖潰決,造成下游近2 500人死亡[4]。

    目前,對滑坡堰塞壩幾何特征的研究較少且多以經(jīng)驗性分析為主:Cui等[5]、Xu等[6]、Yin等[7]根據(jù)堰塞壩高度、庫容及壩體材料特性,對汶川地震災(zāi)區(qū)堰塞壩的危險性進行評估,并將其分為高危、危險、一般、低危等類型:王珊珊等[8]和Zhou等[9]基于GIS技術(shù)提取河谷橫斷面特征參數(shù),評價滑坡堵江的可能性;Casagli和Ermini采用地貌指數(shù)法,建立了以堆積指數(shù)Ib、蓄水指數(shù)Ii[10]和無量綱堆積指數(shù)DBI[11]為控制因素的指數(shù)理論,用來評價堰塞壩的穩(wěn)定性;Komp[12]、Nash等[13]、Dong等[14]、Dal Sasso等[15]均采用地貌指數(shù)法進行了相應(yīng)的研究。除經(jīng)驗性分析方法外,數(shù)值模擬是較為常見的研究滑坡堰塞壩幾何形態(tài)的方法:Hun-gr[16]根據(jù)流變理論,建立了滑坡堰塞壩的動力學(xué)數(shù)值模型;Lo等[17]基于離散元法( DEM),反演分析了臺灣小林村滑坡堰塞壩可能的堆積形態(tài):Zhao等[18]采用CFD-DEM耦合的方法模擬滑坡堵江形成堰塞壩的過程。但現(xiàn)有關(guān)于滑坡堰塞壩幾何特征的研究仍然存在以下兩個問題:

    (1)現(xiàn)有經(jīng)驗分析方法主要是根據(jù)已形成堰塞壩的幾何特征,分析堰塞壩未來的發(fā)展趨勢,并不能對堰塞壩的幾何特征作出預(yù)測。有的經(jīng)驗分析準(zhǔn)確性不夠,例如,KOILip[12]采用地貌指數(shù)法分析新西蘭境內(nèi)滑坡堰塞壩穩(wěn)定性時,結(jié)果并不理想。

    (2)數(shù)值模擬研究不成熟。Hungr的研究并不適合滑坡體含水量較低、未發(fā)生液化的情況;Lo等的研究采用了小林村所在區(qū)域的真三維地形,這增加了模擬結(jié)果的真實性和可靠性,但由于地形過于特殊化和具體化,這種研究往往只能用于具體案例的反演分析,研究結(jié)果并非普遍適用,而且真三維地形模擬需要準(zhǔn)確的地形數(shù)據(jù)(包括滑坡前和滑坡后);Dong等[19]指出,準(zhǔn)確的數(shù)字地形數(shù)據(jù)在絕大多數(shù)情況下是無法短時間內(nèi)獲得的:Zhao等的研究主要關(guān)注流固耦合的作用,未考慮地形條件等外在因素的影響,且該研究設(shè)置的離散元顆粒生成方式與野外真實的滑坡堰塞壩形成過程有較大區(qū)別。

    1 堰塞壩幾何特征與潰口的關(guān)系

    幾何特征對滑坡堰塞壩的潰決破壞有重要影響,在地貌指數(shù)法中,滑坡堰塞壩的高度、長度、寬度和體積等變量常常作為堰塞壩穩(wěn)定性分析的關(guān)鍵因素[19]。除了用于堰塞壩的穩(wěn)定性分析外,滑坡堰塞壩的幾何特征還是研究堰塞壩的潰決機理和過程以及開展模型試驗的基礎(chǔ)。由于滑坡堰塞壩的幾何形態(tài)很不規(guī)則,因此不同滑坡堰塞壩漫頂時的潰口位置不同,如著名的疊溪和紅石巖滑坡堰塞壩,如圖1所示。

    絕大多數(shù)滑坡堰塞壩的破壞方式為漫頂侵蝕,由圖1可見,堰塞壩的幾何特征決定了堰塞壩漫頂時的潰口位置,潰口位置又影響漫頂時壩體的侵蝕過程。當(dāng)潰口位于壩體側(cè)面時,壩體受到的側(cè)向侵蝕為單側(cè)侵蝕,當(dāng)潰口位于壩體中間時,側(cè)向侵蝕為雙側(cè)侵蝕。這種因堰塞壩幾何特征引起的侵蝕差異對于滑坡堰塞壩的潰決過程和輸沙過程有重要意義。筆者曾做過驗證性試驗,模擬側(cè)向侵蝕差異對潰壩過程的影響,兩個壩體除初始潰口外,其他條件均相同,結(jié)果如圖2所示,可以看出,雙側(cè)侵蝕的效果明顯強于單側(cè)侵蝕,進而導(dǎo)致潰口的擴展速度和輸沙濃度等參數(shù)也不相同。

    綜上所述,滑坡堰塞壩的幾何特征是堰塞壩的穩(wěn)定性分析、破壞侵蝕過程分析以及開展模型試驗的基礎(chǔ),研究滑坡堰塞壩的幾何特征具有重要意義。

    滑坡堰塞壩的形成過程是一個復(fù)雜且有待深入了解的山坡一溝谷耦合過程[12.20].很多因素都可能對滑坡堰塞壩的幾何特征產(chǎn)生影響,包括滑坡類型、滑坡的運動以及堆積過程等[19]。在滑坡堰塞壩的所有誘發(fā)因素中,地震是最主要的誘發(fā)因素之-[l],汶川地震后,整個災(zāi)區(qū)地層出現(xiàn)松動,形成了上千個不穩(wěn)定斜坡[21],很多不穩(wěn)定斜坡失穩(wěn)后都有可能形成新的滑坡堰塞壩,這給次生地質(zhì)災(zāi)害的發(fā)生提供了條件。如果能預(yù)測這些滑坡堰塞壩的幾何特征,將為災(zāi)害識別和應(yīng)急搶險工作提供有力參考。

    2 研究方法

    采用離散元法( DEM)對滑坡堰塞壩的幾何特征進行初步分析,顆粒接觸采用Hertz-Mindlin with JKR模型。

    離散元法通過力一位移法則計算顆粒間的接觸力,將顆粒運動當(dāng)作顆粒振動進行求解,并根據(jù)牛頓第二定律分析顆粒運動(考慮重力的影響),其控制方程為

    2.2 Hertz-Mindlin with JKR接觸模型

    Hertz-Mindlin with JKR模型是將Hertz-Mindlin無滑動模型與JKR(Johnson-Kendall-Roberts)黏聚力模型[22]相結(jié)合的一種顆粒接觸模型,該模型的優(yōu)點是可以考慮黏聚力的影響,模型中顆粒切向接觸力、法向阻尼力和切向阻尼力的計算方法與Hertz-Mindlin無滑動接觸模型的計算方法相同[23]。但由于考慮了顆粒表面能(y)的影響,因此顆粒的接觸半徑(a)比Hertz-Mindlin無滑動模型中的接觸半徑(ao)要大。法向力( FJKH)和重疊量(δ,可以為負值,表示顆粒未發(fā)生物理接觸)的計算公式為

    當(dāng)表面能γ=0時,Hertz-Mindlin with JKR模型即為Hertz-Mindlin無滑動接觸模型。

    在巖土工程數(shù)值計算領(lǐng)域,Hertz-Mindlin with JKR模型可以考慮土體含水量和黏性顆粒的影響,因此具有較高的適用性。

    3 數(shù)值模擬

    3.1 物理模型

    一般地,滑坡堰塞壩的動力過程包括:坡體失穩(wěn)一加速運動一減速堆積等過程,如圖4(a)所示。由于數(shù)值模擬考慮所有過程的計算量過大,因此有必要進行物理模型概化,本研究不考慮坡體啟動和加速運動過程,通過API(Application Programming Interface)編程直接賦予坡體相應(yīng)的速度,模擬滑坡體以某個速度滑人溝道并堆積成壩,滑坡體的形態(tài)也設(shè)為較規(guī)則的塊體,如圖4(b)所示。

    3.2 影響因素

    滑坡堰塞壩的幾何特征受多種因素影響,筆者在此建議將各類影響因素按能空間尺度和能量尺度進行分類,以便今后開展滑坡堰塞壩動力過程關(guān)鍵因素的提取。

    從空間尺度和能量尺度而言,溝道斷面形態(tài)、滑坡速度(ν)以及溝床坡度(θ)是決定滑坡堰塞壩幾何特征的基本參數(shù),本研究主要考慮這三種因素對滑坡堰塞壩幾何特征的影響。

    3.3 模擬條件

    由于滑坡的規(guī)模變化非常大,從幾萬到幾百萬甚至上千萬立方米,因此本研究采用單位體積分析,滑坡體總體積0.5 m3.考慮到今后開展室內(nèi)模型試驗(滑槽模擬)的需要,此處將滑槽與溝道統(tǒng)稱為幾何體,其材質(zhì)按鋼材建模。滑坡體位于一矩形滑槽內(nèi),尺寸為1.0mxl.0 mx0.6 m(長×寬×高),溝道斷面形態(tài)包括矩形、梯形和V形3種(見圖5),為保證滑坡物質(zhì)能夠完全堵斷溝道,3種溝道的斷面面積與滑槽斷面面積相等( 0.6 m2)。同時,為保證不同溝道內(nèi)的計算結(jié)果具有可比性,溝道高度保持不變(0.6 m),即滑坡體的重力勢能保持不變,滑槽與溝道的夾角保持為30°,溝道縱向長度為5.0 m,滑坡體滑人溝道的角度與滑槽的底部平行。

    滑坡堰塞壩的動力過程實質(zhì)上是能量的釋放和耗散過程,當(dāng)勢能保持不變時,坡體動能是影響滑坡堰塞壩動力過程的根本性因素。本研究設(shè)置的滑坡體滑人溝道的速度為單位速度(1 m/s)的整數(shù)倍,分別為1、2、3 m/s,代表低速、中速、高速(相對于溝道尺寸而言)三種條件,動能比值為1:4:9。溝床坡度設(shè)置為00、50和100三種,其中:00針對底床坡度較小的溝道或河道(統(tǒng)稱溝道),50和100針對底床坡度較大的普通溝道。

    3.4 參數(shù)標(biāo)定

    巖土工程數(shù)值計算對輸人參數(shù)較敏感,為保證模擬結(jié)果的可靠性,必須對輸人參數(shù)進行標(biāo)定。根據(jù)2.2節(jié)的介紹.Hertz-Mindlin with JKR模型中,影響顆粒接觸力的主要因素有表面能、靜摩擦因數(shù)、滾動摩擦因數(shù)、恢復(fù)系數(shù)等。本研究以驗證性試驗(圖2)中的土體材料作為標(biāo)定對象,該土體材料取自都江堰虹口鄉(xiāng)銀洞子溝滑坡堰塞壩,其力學(xué)參數(shù)已通過大型直接剪切試驗測出,黏聚力(c)和內(nèi)摩擦角(φ)分別為6.2kPa和32°。由于土體黏性較低,因此參數(shù)標(biāo)定方法以數(shù)值計算顆粒的自然堆積角(p)等于土體內(nèi)摩擦角時的參數(shù)組合作為數(shù)值計算的輸人參數(shù)。參數(shù)標(biāo)定結(jié)果見表1.所有輸入?yún)?shù)見表2.

    4.1 溝道斷面形態(tài)的影響

    因27種組合的幾何參數(shù)多達243個,此處僅以高速條件(相對于溝道尺寸而言)下,平直溝道內(nèi)滑坡堰塞壩的幾何特征為例進行分析。該條件下不同溝道內(nèi)滑坡堰塞壩的幾何數(shù)據(jù)見表3。

    由表3可見,高速條件下.V形溝道內(nèi)的堰塞壩縱剖面底寬最大,梯形次之,矩形最小;與此規(guī)律相反,堰塞壩的壩坡傾角和內(nèi)角以矩形最大,梯形次之.V形最小。上述規(guī)律說明:V形溝道內(nèi)滑坡堰塞壩向溝道上、下游延伸堆積的范圍最大。同時,由表3還可看出,高速條件下,矩形溝道內(nèi)堰塞壩的壩頂寬度Lt=0,這說明堰塞壩的縱剖面為三角形,但梯形和V形溝道內(nèi)堰塞壩縱剖面均為梯形。低速條件下,三種溝道內(nèi)的堰塞壩縱剖面如圖7所示,全為梯形。由此可見,溝道斷面形態(tài)對堰塞壩的幾何形態(tài)有直接影響。

    4.2 滑坡速度的影響

    此處僅以溝床坡度θ=0°時梯形溝道內(nèi)堰塞壩的幾何特征為例進行分析,不同滑坡速度條件下堰塞壩的幾何數(shù)據(jù)見表4。

    低速條件形成的堰塞壩,hmax與hmin的比值達到4:1,說明堰塞壩的橫斷面存在很大高差,這類堰塞壩比較常見,如Young河滑坡堰塞壩(圖8(a))和都江堰虹口鄉(xiāng)銀洞子溝滑坡堰塞壩(圖8(b))均屬該類型,其中銀洞子溝滑坡堰塞壩剛形成時的高差接近160m。橫斷面存在巨大高差的堰塞壩若發(fā)生漫頂潰決,則其潰口側(cè)壁失穩(wěn)滑動的規(guī)模會越來越大。

    隨著滑坡速度的增大,堰塞壩橫斷面的高差hmax-hmin逐漸減小。由不同速度條件下梯形溝道內(nèi)堰塞壩的橫斷面圖(如圖9所示)可知,隨著滑坡速度增大,滑坡體可以向?qū)Π杜栏叨纬神R鞍狀的斷面形態(tài),這與圖1中的兩個真實案例相吻合。由此可見,滑坡速度對堰塞壩漫頂潰決時的潰口位置有重要影響,低速條件下形成的堰塞壩,初始潰口位于壩體側(cè)面,側(cè)向侵蝕屬單側(cè)侵蝕(如圖l(a)、圖8所示),但當(dāng)滑坡速度增大后,堰塞壩的初始潰口向中部移動,側(cè)向侵蝕變?yōu)殡p側(cè)侵蝕(如圖l(b)所示)。

    隨著滑坡速度的增大,堰塞壩頂寬(Lt)逐漸減小,堰塞壩頂寬與底寬(Lb)的比值亦減小。若以(Lt+Lb)/2表示堰塞壩的有效厚度,則隨著滑坡速度的增大,堰塞壩的有效厚度變小。

    4.3 溝床坡度的影響

    通過對試驗結(jié)果的分析可知,溝床坡度主要對堰塞壩的縱向堆積產(chǎn)生影響。27種組合對應(yīng)的壩坡傾角和內(nèi)角見表5。

    4.3.1 壩坡內(nèi)角(φ'u、φ'd

    根據(jù)表5的統(tǒng)計結(jié)果,可得出以下認識:

    (1)溝床坡度一定時,堰塞壩壩坡內(nèi)角的大小關(guān)系是矩形溝道≥梯形溝道≥V形溝道。V形溝道內(nèi)的堰塞壩壩坡內(nèi)角大小不受滑坡速度的影響。

    (2)整體而言,溝床坡度增大,堰塞壩下游壩坡內(nèi)角(φ'd)減小,上游壩坡內(nèi)角(φ'u)增大,φ'u與φ'd的比值減小,如圖10所示。

    5 結(jié)論

    基于離散元法( DEM)對滑坡堰塞壩的幾何特征進行模擬,采用可以考慮顆粒黏性的Hertz-Mindlinwith JKR模型進行了27種參數(shù)組合的模擬分析,得到以下幾點主要認識。

    (1)滑坡速度和溝道形態(tài)對滑坡堰塞壩的橫向和縱向幾何形態(tài)都有影響,而溝床坡度主要影響滑坡堰塞壩的縱向幾何特征。

    (2)低速條件下形成的堰塞壩在橫向存在較大高差,堰塞壩若發(fā)生漫頂潰決,初始潰口位于堰塞壩側(cè)面,側(cè)向侵蝕為單側(cè)侵蝕。當(dāng)滑坡速度達到高速條件后,滑坡體可以向?qū)Π杜栏叨纬神R鞍狀的橫斷面,這類堰塞壩若出現(xiàn)漫頂潰決,初始潰口位于壩體中部,側(cè)向侵蝕為雙側(cè)侵蝕。

    (3)滑坡速度增大,堰塞壩縱剖面頂寬(Lt)和底寬(Lb)的比值減小,堰塞壩的有效厚度減小。

    (4)溝床坡度增大,堰塞壩下游壩坡內(nèi)角減小,上游壩坡內(nèi)角增大,兩者的比值減小。但是堰塞壩壩坡傾角的變化規(guī)律和內(nèi)角的變化規(guī)律相反。

    參考文獻:

    [1]COSTA J E,SCHUSTER R L.The Formation and Failure ofNatural Dams[J].Ceological Society of America Bulletin,1988, 100(7): 1054-1068.

    [2] 柴賀軍,劉浩吾,劉漢超,等,天然土石壩穩(wěn)定性初步研究[J].地質(zhì)科技情報,2001, 20(1):77-81.

    [3]DAI F C,LEE C F,DENC J H,et al.The 1786 Earthquake-Triggered Landslide Dam and Subsequent Dam-Break Flood onthe Dadu River, Southwestem China [J]. Geomorphology,2005, 65(3):205-221.

    [4] 王蘭生,楊立錚,李天斌,等,四川岷江疊溪較場地震滑坡及環(huán)境保護[J].地質(zhì)災(zāi)害與環(huán)境保護,2000, 11(3):195-199.

    [5]CUI P,ZHUY Y,HAN Y S,et a1.The 12 May WenchuanEarthquake-Induced Landslide Lakes: Distribution and Prelim-inary Risk Evaluation[J]. Landslides, 2009,6(3):209-223.

    [6] XU Q,F(xiàn)AN X M, HUANC R Q,et al.Landslide DamsTriggered by the Wenchuan Earthquake, Sichuan Province,South West China[J].Bulletin of Engineering Ceology andthe Environment, 2009, 68(3): 373-386.

    [7]YIN Y P,WANC F W, SUN P.Landslide Hazards Triggeredby the 2008 Wenchuan Earthquake, Sichuan, China[J]. Land-slides, 2009, 6(2): 139-152.

    [8] 王珊珊,童立強,基于河谷橫剖面形態(tài)特征的滑坡體堵江易發(fā)性評價研究[J].地理與地理信息科學(xué),2016, 32(5):97-102.

    [9]ZHOU L,TAKASHI 0.Longitudinal and Transverse Profilesof Hilly and Mountainous Watersheds in Japan[J].Geomor-phology, 2009, 111(1-2): 17-26.

    [10]CASACLI N,ERMINI L.Ceomorphic Analysis of LandslideDams in the Northem Apennine[J].Japanese Ceomorpho-logical Union, 1999, 20(3):219-249.

    [11]ERMINL,CASACLIN.Prediction of the Behavior of Land-slide Dams Using A Ceomorphologcal Dimensionless Index[J]. Earth Surface Processes and Landforms, 2003, 28(1):31-47.

    [12]KORUP 0.Geomorphometric Characteristics of New ZealandLandslide Dams[J].Engineering Geology, 2004, 73(1):13- 35.

    [13]NASH T, BELL D, DAVIES T, et al. Analysis of the For-mation and Failure of Ram Creek Landslide Dam, South Is-land, New Zealand [ J]. New Zealand Joumal of Ceologyand Ceophysics , 2008, 51( 3) : 187-193.

    [14]Dong J J, Tung Y H, Chen C C, et al. Logistic RegressionModel for Predicting the Failure Probability of A LandslideDam[J]. Engineering Ceology, 2011, 117( 1) : 52-61.

    [15]DAL SASSO S F, SOLE A, PASCALE S, et al. AssessmentMethodology for the Prediction of Landslide Dam Hazard [J] .Natural Hazards and Earth System Sciences, 2014, 14(3) :557-567.

    [16]HUNCR 0. A Model for the Runout Analysis of Rapid FlowSlides, Debris Flows, and Avalanches [ J]. CanadianCeotechnical Journal, 1995 , 32(4) : 610-623.

    [17]10 C M, LIN M L, TANC C L, et al. A Kinematic Modelof the Hsiaolin Landslide Calibrated to the Morphology ofthe Landslide Deposit [ J]. Engineering Geology, 2011,123(1) : 22-39.

    [18]ZHAO T, DAI F, XU N. Coupled DEM-CFD Investigationon the Formation of Landslide Dams in Narrow Rivers [J] .Landslides, 2017, 14(1) : 189-201.

    [19]DONC J J, LAI P J, CHANC C P, et al. Deriving Land-slide Dam Geometry From Remote Sensing Images for theRapid Assessment of Critical Parameters Related to Dam-Breach Hazards [ J] . Landslides, 2014, 11( 1) : 93-105.

    [20]COSTA J E, SCHUSIER R L. Documented Historical LandslideDams from Around the World[ R] . Virginia : US Ceologcal Sur-vey, 1991:91-239.

    [21]殷躍 .汶川八級震地質(zhì)災(zāi)害研究 [J] .工程地質(zhì)學(xué)報,2008, 16(4) : 433-444.

    [22] JOHNSON K L, KENDALL K, ROBERTS A D. Surface Energyand the Contact of Ilastic Solids [ C] //Proceedings of the RoyalSociety of London A: Mathematical, Physical and Engneering Sciences. London : The Royal Society, 1971: 301-313.

    [23] SANTOS K G, CAMPOS A V P, OLIVEIRA O S, et al.Dem Simulations of Dynamic An~e of Repose of AcerolaResidue: A Parametric Study Using a Response SurfaceTechnique [J] . Blucher Chemical Engineering Proceedings ,2015 , 1( 2) : 11326-11333.

    [責(zé)任編輯張帥]

    猜你喜歡
    泥石流滑坡
    水庫滑坡地下水動態(tài)響應(yīng)規(guī)律及浸潤線計算模型——以石榴樹包滑坡為例
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    泥石流
    雜文月刊(2018年21期)2019-01-05 05:55:28
    “民謠泥石流”花粥:唱出自己
    海峽姐妹(2017年6期)2017-06-24 09:37:36
    泥石流
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    基于Fluent的滑坡入水過程數(shù)值模擬
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    機械班長
    奪命泥石流
    亚洲精品国产精品久久久不卡| 国产av精品麻豆| 中亚洲国语对白在线视频| 国产一区二区在线观看av| 午夜免费鲁丝| 久久久久久免费高清国产稀缺| 在线观看一区二区三区激情| 99国产综合亚洲精品| 操出白浆在线播放| 欧美久久黑人一区二区| 欧美精品亚洲一区二区| 国产又爽黄色视频| 满18在线观看网站| www.精华液| 日韩免费高清中文字幕av| 在线观看www视频免费| 午夜福利乱码中文字幕| 国产精品 国内视频| 中文字幕色久视频| 岛国在线观看网站| 日本91视频免费播放| 考比视频在线观看| 汤姆久久久久久久影院中文字幕| 国产一区二区 视频在线| 曰老女人黄片| 日本欧美视频一区| 69av精品久久久久久 | 久久午夜综合久久蜜桃| 97人妻天天添夜夜摸| 国产一级毛片在线| 国产精品久久久久久精品古装| 欧美日韩亚洲高清精品| 国产在线一区二区三区精| 日韩欧美一区二区三区在线观看 | 91大片在线观看| 亚洲成国产人片在线观看| 99久久精品国产亚洲精品| 久久精品国产亚洲av高清一级| 高清视频免费观看一区二区| 新久久久久国产一级毛片| 亚洲少妇的诱惑av| 人人妻人人澡人人爽人人夜夜| 国产91精品成人一区二区三区 | 日韩三级视频一区二区三区| 人人妻人人澡人人爽人人夜夜| 亚洲人成77777在线视频| 视频区欧美日本亚洲| 亚洲国产精品999| 少妇被粗大的猛进出69影院| 亚洲精品粉嫩美女一区| 久久天堂一区二区三区四区| 久久久精品国产亚洲av高清涩受| 欧美精品一区二区免费开放| 精品久久久久久电影网| 极品人妻少妇av视频| av又黄又爽大尺度在线免费看| 三级毛片av免费| 国产成人免费观看mmmm| 国产精品免费大片| 91大片在线观看| 日本av手机在线免费观看| 99国产精品99久久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美人与性动交α欧美软件| 99国产精品一区二区三区| 日韩一区二区三区影片| 国产视频一区二区在线看| 建设人人有责人人尽责人人享有的| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美一区二区三区在线观看 | 嫩草影视91久久| 成年人午夜在线观看视频| 91九色精品人成在线观看| 另类精品久久| 青春草视频在线免费观看| 性高湖久久久久久久久免费观看| 国产av又大| 多毛熟女@视频| 亚洲自偷自拍图片 自拍| 在线观看舔阴道视频| 日韩有码中文字幕| 国产精品成人在线| 老汉色av国产亚洲站长工具| cao死你这个sao货| 狠狠狠狠99中文字幕| 国产精品一区二区在线观看99| 俄罗斯特黄特色一大片| 国产免费现黄频在线看| 99国产精品免费福利视频| 国产在线视频一区二区| 日韩欧美一区视频在线观看| 国产av一区二区精品久久| 考比视频在线观看| 在线观看免费视频网站a站| 亚洲欧美色中文字幕在线| 久久精品国产a三级三级三级| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美一区视频在线观看| 亚洲,欧美精品.| 精品少妇黑人巨大在线播放| 免费高清在线观看日韩| 精品国产一区二区三区久久久樱花| 色精品久久人妻99蜜桃| 亚洲国产欧美一区二区综合| 国产精品九九99| 一级片免费观看大全| 色婷婷av一区二区三区视频| av在线播放精品| 欧美日韩精品网址| 每晚都被弄得嗷嗷叫到高潮| 国产高清视频在线播放一区 | 90打野战视频偷拍视频| 欧美xxⅹ黑人| 母亲3免费完整高清在线观看| 久久天躁狠狠躁夜夜2o2o| 十八禁人妻一区二区| 欧美人与性动交α欧美精品济南到| 18禁黄网站禁片午夜丰满| 亚洲综合色网址| 国产片内射在线| 精品高清国产在线一区| 国产精品99久久99久久久不卡| 搡老熟女国产l中国老女人| 18禁观看日本| 亚洲精品中文字幕一二三四区 | 在线 av 中文字幕| 亚洲视频免费观看视频| 中文字幕人妻熟女乱码| 一进一出抽搐动态| 久久久国产成人免费| 美女扒开内裤让男人捅视频| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品一区三区| 精品国产国语对白av| 亚洲情色 制服丝袜| 亚洲va日本ⅴa欧美va伊人久久 | 免费在线观看完整版高清| 免费在线观看完整版高清| 国产精品av久久久久免费| 天天操日日干夜夜撸| 亚洲一区二区三区欧美精品| 女警被强在线播放| 夫妻午夜视频| 丝袜美足系列| 午夜两性在线视频| 999久久久国产精品视频| 国产在线观看jvid| 亚洲精华国产精华精| 欧美久久黑人一区二区| 免费黄频网站在线观看国产| 91国产中文字幕| 一级毛片女人18水好多| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 一区二区三区四区激情视频| 精品一区二区三卡| 69精品国产乱码久久久| 国产精品熟女久久久久浪| 在线永久观看黄色视频| 成人黄色视频免费在线看| 男男h啪啪无遮挡| 亚洲成人免费av在线播放| 大型av网站在线播放| 成人影院久久| 精品人妻一区二区三区麻豆| 这个男人来自地球电影免费观看| 在线观看舔阴道视频| 日韩精品免费视频一区二区三区| 黄片大片在线免费观看| 精品久久久精品久久久| 亚洲第一青青草原| 久久久久久亚洲精品国产蜜桃av| 少妇 在线观看| 精品人妻熟女毛片av久久网站| 精品国产一区二区久久| 99热全是精品| 叶爱在线成人免费视频播放| 淫妇啪啪啪对白视频 | av超薄肉色丝袜交足视频| 久久久久久久精品精品| 国产一区有黄有色的免费视频| 欧美精品一区二区大全| 中文字幕av电影在线播放| 在线观看人妻少妇| 久9热在线精品视频| 一本久久精品| www.999成人在线观看| 精品久久久久久电影网| 久久久久国产精品人妻一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 在线看a的网站| 搡老乐熟女国产| 纯流量卡能插随身wifi吗| 免费女性裸体啪啪无遮挡网站| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 大香蕉久久成人网| 热99国产精品久久久久久7| 男人添女人高潮全过程视频| 男女下面插进去视频免费观看| 最新在线观看一区二区三区| 首页视频小说图片口味搜索| 老司机福利观看| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区在线观看国产| 成人影院久久| 黑人巨大精品欧美一区二区mp4| 亚洲成人免费av在线播放| 日韩,欧美,国产一区二区三区| 男女之事视频高清在线观看| 操出白浆在线播放| svipshipincom国产片| 真人做人爱边吃奶动态| 欧美 日韩 精品 国产| 亚洲成人国产一区在线观看| 满18在线观看网站| 首页视频小说图片口味搜索| 美女福利国产在线| 97精品久久久久久久久久精品| 丝袜喷水一区| 国产三级黄色录像| 男女边摸边吃奶| 一进一出抽搐动态| 高清欧美精品videossex| h视频一区二区三区| 电影成人av| 高清在线国产一区| 超碰97精品在线观看| 老汉色av国产亚洲站长工具| 中亚洲国语对白在线视频| 一级片'在线观看视频| 五月天丁香电影| 搡老岳熟女国产| 19禁男女啪啪无遮挡网站| 欧美日韩亚洲综合一区二区三区_| 69精品国产乱码久久久| av又黄又爽大尺度在线免费看| 夜夜骑夜夜射夜夜干| 在线 av 中文字幕| 嫩草影视91久久| 免费黄频网站在线观看国产| 日韩欧美一区视频在线观看| 色94色欧美一区二区| 国产成人av教育| 日本av手机在线免费观看| 国产日韩欧美视频二区| 欧美在线黄色| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费午夜福利视频| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 在线观看人妻少妇| 免费高清在线观看日韩| 三级毛片av免费| 中文字幕av电影在线播放| 欧美 日韩 精品 国产| 99国产精品免费福利视频| 欧美日韩成人在线一区二区| 婷婷丁香在线五月| 在线十欧美十亚洲十日本专区| 少妇精品久久久久久久| 精品福利永久在线观看| 日本撒尿小便嘘嘘汇集6| 国产成人a∨麻豆精品| 999精品在线视频| 无遮挡黄片免费观看| 国产1区2区3区精品| 悠悠久久av| 少妇裸体淫交视频免费看高清 | 亚洲 国产 在线| 国产男女超爽视频在线观看| 首页视频小说图片口味搜索| 青春草亚洲视频在线观看| 少妇猛男粗大的猛烈进出视频| 精品少妇一区二区三区视频日本电影| 男女无遮挡免费网站观看| 免费观看av网站的网址| 欧美黄色淫秽网站| 我要看黄色一级片免费的| 日本欧美视频一区| 久久国产精品大桥未久av| 久久久国产成人免费| 久久这里只有精品19| 窝窝影院91人妻| 又紧又爽又黄一区二区| 日韩欧美国产一区二区入口| 黄片小视频在线播放| 亚洲欧美成人综合另类久久久| 亚洲欧美清纯卡通| 两个人免费观看高清视频| 亚洲五月婷婷丁香| 亚洲精品av麻豆狂野| 97人妻天天添夜夜摸| 久久精品熟女亚洲av麻豆精品| 亚洲一码二码三码区别大吗| 99久久精品国产亚洲精品| av电影中文网址| 日本91视频免费播放| 精品国产乱码久久久久久男人| 亚洲免费av在线视频| 黄网站色视频无遮挡免费观看| 俄罗斯特黄特色一大片| 精品卡一卡二卡四卡免费| 午夜福利一区二区在线看| 国产一区二区三区av在线| 午夜成年电影在线免费观看| 亚洲人成电影观看| 麻豆av在线久日| 老司机午夜十八禁免费视频| 亚洲美女黄色视频免费看| cao死你这个sao货| 久久久久久久精品精品| 在线观看免费视频网站a站| 国产欧美日韩综合在线一区二区| 欧美精品亚洲一区二区| av有码第一页| 日韩大片免费观看网站| 久久久精品区二区三区| 美女高潮到喷水免费观看| 国产亚洲欧美精品永久| 91字幕亚洲| 日韩人妻精品一区2区三区| 国产人伦9x9x在线观看| 久久精品国产a三级三级三级| av网站免费在线观看视频| 日日夜夜操网爽| 国产在线观看jvid| 汤姆久久久久久久影院中文字幕| 美女福利国产在线| 午夜免费观看性视频| 日韩人妻精品一区2区三区| 午夜免费鲁丝| 亚洲av成人不卡在线观看播放网| 欧美成人性av电影在线观看| 国产精品一区二区免费欧美| 亚洲美女黄片视频| 国产亚洲av嫩草精品影院| 熟妇人妻久久中文字幕3abv| 十八禁网站免费在线| 老司机福利观看| 欧美一区二区精品小视频在线| 最近最新免费中文字幕在线| 亚洲国产欧美网| 一进一出好大好爽视频| 亚洲成人久久爱视频| 777久久人妻少妇嫩草av网站| 亚洲av日韩精品久久久久久密| 国内久久婷婷六月综合欲色啪| 嫩草影视91久久| 午夜福利在线观看吧| 免费观看精品视频网站| 中文字幕人成人乱码亚洲影| 狂野欧美白嫩少妇大欣赏| 国产精品1区2区在线观看.| 亚洲美女视频黄频| 一进一出抽搐动态| 久热爱精品视频在线9| 九九热线精品视视频播放| 中文字幕最新亚洲高清| 97人妻精品一区二区三区麻豆| 男插女下体视频免费在线播放| 免费在线观看日本一区| 夜夜看夜夜爽夜夜摸| 亚洲av美国av| www.www免费av| 人人妻人人看人人澡| 男女下面进入的视频免费午夜| 亚洲国产欧美一区二区综合| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 欧美3d第一页| 亚洲国产高清在线一区二区三| or卡值多少钱| 老司机午夜福利在线观看视频| 精品一区二区三区视频在线观看免费| 日韩欧美免费精品| www国产在线视频色| 国产伦在线观看视频一区| 久久精品aⅴ一区二区三区四区| 99riav亚洲国产免费| 日本在线视频免费播放| 精品福利观看| 国产精品免费一区二区三区在线| 久久九九热精品免费| 亚洲一区二区三区不卡视频| 搡老妇女老女人老熟妇| 成年版毛片免费区| 亚洲精品久久国产高清桃花| 夜夜看夜夜爽夜夜摸| 禁无遮挡网站| 日韩欧美 国产精品| 99精品久久久久人妻精品| 午夜视频精品福利| x7x7x7水蜜桃| 亚洲自拍偷在线| 性欧美人与动物交配| 少妇人妻一区二区三区视频| 欧美一级毛片孕妇| av欧美777| √禁漫天堂资源中文www| 神马国产精品三级电影在线观看 | 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 午夜视频精品福利| 90打野战视频偷拍视频| 国产又黄又爽又无遮挡在线| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 曰老女人黄片| 国产高清视频在线观看网站| 国产区一区二久久| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 国内精品久久久久精免费| 亚洲成人中文字幕在线播放| 日韩有码中文字幕| 国产视频一区二区在线看| 狠狠狠狠99中文字幕| 精品久久久久久久人妻蜜臀av| 美女黄网站色视频| 蜜桃久久精品国产亚洲av| 91字幕亚洲| 国产av又大| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 欧美日本视频| 88av欧美| 亚洲天堂国产精品一区在线| 午夜a级毛片| 老司机午夜福利在线观看视频| 黄色 视频免费看| 18禁裸乳无遮挡免费网站照片| 亚洲午夜精品一区,二区,三区| 精品久久久久久久毛片微露脸| www.精华液| 人妻夜夜爽99麻豆av| www.999成人在线观看| 久久亚洲精品不卡| 制服诱惑二区| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 无限看片的www在线观看| 99国产综合亚洲精品| 精品高清国产在线一区| 男女床上黄色一级片免费看| 国产aⅴ精品一区二区三区波| 巨乳人妻的诱惑在线观看| 免费在线观看成人毛片| 国产爱豆传媒在线观看 | 国产99久久九九免费精品| 久久精品夜夜夜夜夜久久蜜豆 | 国内精品一区二区在线观看| 90打野战视频偷拍视频| 精品国产乱子伦一区二区三区| 亚洲av美国av| 国产单亲对白刺激| 午夜福利在线在线| 真人一进一出gif抽搐免费| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久午夜综合久久蜜桃| 婷婷六月久久综合丁香| 国产亚洲欧美98| 国产精品综合久久久久久久免费| 俺也久久电影网| 亚洲国产精品合色在线| 国产精品精品国产色婷婷| 亚洲在线自拍视频| 国产一区在线观看成人免费| 久久九九热精品免费| 黑人巨大精品欧美一区二区mp4| 亚洲一码二码三码区别大吗| 一本综合久久免费| 国产精品自产拍在线观看55亚洲| 日日爽夜夜爽网站| 岛国视频午夜一区免费看| 夜夜躁狠狠躁天天躁| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品一区二区www| 亚洲精品中文字幕一二三四区| 精品不卡国产一区二区三区| 18禁黄网站禁片午夜丰满| 少妇人妻一区二区三区视频| 国产人伦9x9x在线观看| 视频区欧美日本亚洲| 国产精品乱码一区二三区的特点| 男插女下体视频免费在线播放| 男人舔女人的私密视频| 可以免费在线观看a视频的电影网站| 日韩精品免费视频一区二区三区| 夜夜躁狠狠躁天天躁| 嫩草影院精品99| 老司机在亚洲福利影院| 亚洲一卡2卡3卡4卡5卡精品中文| 国产野战对白在线观看| 国产高清videossex| 国产黄a三级三级三级人| 国产成+人综合+亚洲专区| 亚洲黑人精品在线| 三级国产精品欧美在线观看 | 国产激情偷乱视频一区二区| 两个人视频免费观看高清| 日本黄色视频三级网站网址| 老司机午夜十八禁免费视频| 国产精品香港三级国产av潘金莲| 日韩大尺度精品在线看网址| 我要搜黄色片| 麻豆成人av在线观看| 一进一出抽搐动态| 一区福利在线观看| 人妻丰满熟妇av一区二区三区| 欧美日韩亚洲综合一区二区三区_| 黄色成人免费大全| 国产麻豆成人av免费视频| 午夜福利欧美成人| 脱女人内裤的视频| 日日干狠狠操夜夜爽| 在线观看免费视频日本深夜| 久久精品成人免费网站| 777久久人妻少妇嫩草av网站| cao死你这个sao货| 国产黄片美女视频| 18美女黄网站色大片免费观看| 麻豆成人av在线观看| 99久久国产精品久久久| 在线观看一区二区三区| 国产精品电影一区二区三区| 最近最新中文字幕大全免费视频| 麻豆久久精品国产亚洲av| 最新在线观看一区二区三区| 亚洲激情在线av| 天堂av国产一区二区熟女人妻 | 亚洲成av人片在线播放无| av超薄肉色丝袜交足视频| 欧美三级亚洲精品| a在线观看视频网站| 亚洲一区二区三区色噜噜| 黄色片一级片一级黄色片| 久久精品夜夜夜夜夜久久蜜豆 | 宅男免费午夜| 女人爽到高潮嗷嗷叫在线视频| 久久精品成人免费网站| av天堂在线播放| 亚洲欧洲精品一区二区精品久久久| 亚洲无线在线观看| 又紧又爽又黄一区二区| aaaaa片日本免费| 午夜免费观看网址| 99久久99久久久精品蜜桃| 久久久久国内视频| 午夜影院日韩av| 欧美久久黑人一区二区| 国产黄色小视频在线观看| 亚洲中文字幕日韩| 国产单亲对白刺激| 国产av不卡久久| 欧美激情久久久久久爽电影| 精品熟女少妇八av免费久了| 一区福利在线观看| 中文字幕高清在线视频| 欧美成人一区二区免费高清观看 | 久久精品综合一区二区三区| 国产精品1区2区在线观看.| 美女 人体艺术 gogo| 午夜福利欧美成人| 99国产精品99久久久久| 久久精品成人免费网站| 2021天堂中文幕一二区在线观| 天堂影院成人在线观看| 叶爱在线成人免费视频播放| 动漫黄色视频在线观看| 国产成人欧美在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品合色在线| 亚洲成人久久爱视频| 国产av一区二区精品久久| 9191精品国产免费久久| 日韩高清综合在线| 亚洲国产精品999在线| 国产成人av激情在线播放| 国产成人精品无人区| 大型黄色视频在线免费观看| 正在播放国产对白刺激| 欧美日韩乱码在线| 每晚都被弄得嗷嗷叫到高潮| √禁漫天堂资源中文www| 舔av片在线| 亚洲精品在线美女| 校园春色视频在线观看| 国产三级黄色录像| 国产精品国产高清国产av| 99久久国产精品久久久| 国产私拍福利视频在线观看| 99国产精品一区二区三区| 在线观看日韩欧美| 色综合亚洲欧美另类图片| 国产伦一二天堂av在线观看| 麻豆国产av国片精品| 久久这里只有精品中国| 韩国av一区二区三区四区| 久久午夜亚洲精品久久| av欧美777| 91九色精品人成在线观看| 亚洲色图 男人天堂 中文字幕| 99re在线观看精品视频| 久久久久久久久中文| 久久久久久久久中文| 十八禁人妻一区二区| 丁香欧美五月| 国产精品电影一区二区三区| 免费看十八禁软件| 久久人人精品亚洲av| 欧美最黄视频在线播放免费| 久久精品91无色码中文字幕| 欧美成人午夜精品|