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

    基于斜率模型的突發(fā)型黃土滑坡失穩(wěn)時(shí)間預(yù)測

    2021-07-23 06:13:38方汕澳修德皓趙寬耀李志剛
    水文地質(zhì)工程地質(zhì) 2021年4期
    關(guān)鍵詞:倒數(shù)切線黃土

    方汕澳,許 強(qiáng),修德皓,趙寬耀,李志剛,蒲 楓

    (成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610059)

    我國西北黃土地區(qū)降水稀少、氣候干旱,臨河地區(qū)多采用提水漫灌等灌溉方式,導(dǎo)致地下水位逐漸升高,滑坡等地質(zhì)災(zāi)害頻發(fā)[1]。由灌溉誘發(fā)的滑坡加速階段歷時(shí)短、突發(fā)性強(qiáng)、危害性大[2],因此,開展針對此類突發(fā)型黃土滑坡的監(jiān)測預(yù)警預(yù)報(bào)工作對減少人員傷亡和降低經(jīng)濟(jì)損失具有重要意義。

    已有多位學(xué)者針對突發(fā)型黃土滑坡監(jiān)測預(yù)警預(yù)報(bào)展開了研究。亓星等[3]通過研發(fā)具有智能變頻功能的滑坡位移計(jì)成功獲取了黑方臺(tái)地區(qū)突發(fā)型黃土滑坡完整變形數(shù)據(jù);趙超英等[4]、史緒國等[5]多位學(xué)者利用InSAR技術(shù)實(shí)現(xiàn)了黃土滑坡的早期識(shí)別和時(shí)序監(jiān)測,提出黃土滑坡監(jiān)測新方法;亓星等[6?7]先后通過速度倒數(shù)法和改進(jìn)的Saito模型,對突發(fā)型黃土滑坡失穩(wěn)時(shí)間預(yù)測進(jìn)行了研究;許強(qiáng)等[8]通過建立以變形速率閾值和改進(jìn)切線角為預(yù)警指標(biāo)的四級預(yù)警判據(jù),成功實(shí)現(xiàn)了對甘肅黑方臺(tái)突發(fā)型黃土滑坡的預(yù)警,為解決突發(fā)型黃土滑坡“什么時(shí)間發(fā)生”的難題提供了有效途徑。

    失穩(wěn)時(shí)間預(yù)測是監(jiān)測預(yù)警預(yù)報(bào)工作的重要組成部分[9],先后經(jīng)歷了現(xiàn)象預(yù)測[10?11]、確定性預(yù)測、統(tǒng)計(jì)分析預(yù)測[12?15]、非線性預(yù)測[16?18]、確定性預(yù)測[19]等幾個(gè)階段。確定性預(yù)測主要是基于滑坡運(yùn)動(dòng)參數(shù),通過物理或數(shù)學(xué)方法將室內(nèi)試驗(yàn)數(shù)據(jù)或滑坡監(jiān)測數(shù)據(jù)進(jìn)行推導(dǎo),并不斷總結(jié)改進(jìn)得到的預(yù)測方法,能夠反映滑坡的物理實(shí)質(zhì)[19]。Saito等[20]率先通過大量室內(nèi)試驗(yàn)總結(jié)出理想狀態(tài)下土質(zhì)邊坡失穩(wěn)時(shí)間預(yù)測經(jīng)驗(yàn)公式,并提出了基于蠕變第三階段的圖解法;Fukuzono[21]在室內(nèi)試驗(yàn)的基礎(chǔ)上進(jìn)一步總結(jié)出恒定荷載作用下滑坡臨滑階段位移速度倒數(shù)與時(shí)間的關(guān)系,提出速度倒數(shù)法預(yù)測模型;Mufundirwa等[22]在前人研究的基礎(chǔ)上提出斜率模型(又稱SLOPE模型、SLO模型)來表示巖石蠕變末期應(yīng)變發(fā)散特征,并應(yīng)用到滑坡實(shí)例和巖石蠕變試驗(yàn)中。在這之后有多位學(xué)者[23?24]基于以上預(yù)測模型進(jìn)行了總結(jié)改進(jìn)。

    上述預(yù)測方法多應(yīng)用于室內(nèi)試驗(yàn)或巖質(zhì)滑坡,關(guān)于自然條件下的黃土滑坡的研究很少,特別是突發(fā)型黃土滑坡。本文以2019年黑方臺(tái)發(fā)生的突發(fā)型黃土滑坡為研究對象,通過改進(jìn)的切線角預(yù)警模型確定滑坡變形階段,基于斜率模型提出以切線角為指標(biāo)的簡化計(jì)算方法進(jìn)行滑坡失穩(wěn)時(shí)間預(yù)測研究,以期為突發(fā)型黃土滑坡的監(jiān)測預(yù)警和應(yīng)急管理提供一定的建議。

    1 研究區(qū)概況

    黑方臺(tái)地區(qū)位于甘肅省臨夏州永靖縣內(nèi),地貌以黃土臺(tái)塬、河谷地貌為主。整個(gè)臺(tái)塬由黑臺(tái)和方臺(tái)組成,長7.7 km,寬2.5 km,面積為13.7 km2,地勢西高東低[25?26]。受常年大水漫灌影響,臺(tái)塬地下水位以0.3 m/a左右的速率上升,塬邊滑坡頻發(fā)[8]。

    據(jù)調(diào)查統(tǒng)計(jì),黑方臺(tái)地區(qū)發(fā)育滑坡多達(dá)77處,按發(fā)生位置的不同可將該區(qū)域滑坡劃分為新塬段、黨川段、黃茨段、焦家崖段、焦家段和陳家段共六個(gè)區(qū)段(圖1)。從滑坡成災(zāi)模式來看,黑方臺(tái)地區(qū)發(fā)生的滑坡可劃分為黃土基巖型、黃土泥流型、滑移崩塌型和靜態(tài)液化型滑坡等4類,其中滑移崩塌型和靜態(tài)液化型滑坡具有明顯的突發(fā)性特征[27]。2019年成功預(yù)警了4起突發(fā)型滑坡,主要分布在北側(cè)陳家段和南側(cè)黨川段,其中滑移崩塌型滑坡有陳家6#、黨川6#和黨川7#,靜態(tài)液化型滑坡有黨川4#,具體滑坡預(yù)警信息如表1所示。

    圖1 黑方臺(tái)滑坡分布Fig.1 Location of the Heifangtai landslides

    表1 2019年成功預(yù)警的滑坡Table 1 Basic characteristics of landslides predicted in 2019

    2 研究方法

    2.1 監(jiān)測數(shù)據(jù)的獲取與處理

    本文采用具有智能變頻功能的滑坡位移計(jì)進(jìn)行位移監(jiān)測,該儀器能獲取突發(fā)型黃土滑坡完整變形數(shù)據(jù)。結(jié)合黑方臺(tái)突發(fā)型黃土滑坡的發(fā)育規(guī)律和特征進(jìn)行監(jiān)測儀器的布設(shè)(圖2):對于靜態(tài)液化型滑坡,監(jiān)測儀器布設(shè)在地下水滲水強(qiáng)烈的滑坡后緣附近;對于滑移崩塌型滑坡,監(jiān)測儀器主要集中在多級裂縫發(fā)育的滑坡后緣附近。多次成功預(yù)警的案例也驗(yàn)證了上述布設(shè)原則的可行性[8]。許強(qiáng)[28]認(rèn)為當(dāng)滑坡進(jìn)入加速變形階段,滑坡區(qū)會(huì)以滑動(dòng)塊體形式整體滑動(dòng),各監(jiān)測點(diǎn)的位移矢量方向?qū)②呌谝恢?,選取加速變形階段不同監(jiān)測點(diǎn)的監(jiān)測數(shù)據(jù)進(jìn)行預(yù)警預(yù)報(bào)的差別不大。斜率模型是基于加速變形階段的監(jiān)測數(shù)據(jù)進(jìn)行預(yù)測計(jì)算,故本文所選取的監(jiān)測數(shù)據(jù)能夠反映滑坡變形趨勢,達(dá)到有效預(yù)警預(yù)報(bào)的目的。

    圖2 滑坡監(jiān)測點(diǎn)布設(shè)示意圖Fig.2 Monitoring systems of the landslides

    采用移動(dòng)平均法對監(jiān)測數(shù)據(jù)進(jìn)行濾波處理,以降低儀器噪聲和周圍環(huán)境對監(jiān)測數(shù)據(jù)的影響,反映滑坡變形趨勢。計(jì)算公式如下:

    式中:Si—i時(shí)刻滑坡的累計(jì)位移/mm;

    N—移動(dòng)平均項(xiàng)數(shù);

    —移動(dòng)平均處理后的位移/mm。

    處理后的數(shù)據(jù)預(yù)測標(biāo)準(zhǔn)誤差大小與N的取值呈正線性相關(guān),N的取值越小誤差越小,在這里N取24[8]。

    2.2 滑坡變形階段確定

    通過坐標(biāo)變化進(jìn)行無量綱處理得到的改進(jìn)切線角模型能夠使滑坡變形曲線在相同切線角特征下進(jìn)行比較[29],因此,采用該模型確定滑坡變形階段,計(jì)算公式如下:

    式中:Si—i時(shí)刻滑坡的累計(jì)位移/mm;

    Ti、Ti?1—變換后與時(shí)間相同量綱的縱坐標(biāo)值,兩者相差一個(gè)監(jiān)測周期,如1天、1 周等;

    —滑坡等速變形階段的速率/(mm·d?1);

    αi—改進(jìn)切線角/(°);

    ti、ti?1—監(jiān)測時(shí)刻,兩者相差一個(gè)監(jiān)測周期。

    當(dāng)切線角大于45°時(shí),可認(rèn)為滑坡進(jìn)入初加速變形階段;當(dāng)切線角大于80°時(shí),滑坡進(jìn)入中加速階段;當(dāng)切線角大于85°時(shí),滑坡進(jìn)入加加速階段;切線角越接近90°,滑坡越接近發(fā)生。突發(fā)型黃土滑坡前期變形緩慢且不穩(wěn)定,切線角存在較大幅度的波動(dòng),通過計(jì)算發(fā)現(xiàn),當(dāng)切線角達(dá)到60°后,波動(dòng)性減小且變化趨勢較明顯,亓星等[7]也指出在滑坡初加速變形階段(45°~80°)的前期,距離滑坡實(shí)際發(fā)生時(shí)間還較長,較小的波動(dòng)都會(huì)導(dǎo)致預(yù)測結(jié)果產(chǎn)生較大的差異。故本文以60°切線角為起點(diǎn),每5°為間隔,將滑坡加速階段后期劃分為60°~65°、65°~70°、70°~75°、75°~80°、80°~85°、85°~90°共6 個(gè)變形階段。

    2.3 斜率模型介紹

    SLO失穩(wěn)時(shí)間預(yù)測模型源于Fukui等[30]根據(jù)“加速蠕變理論”總結(jié)出的巖石加速蠕變階段應(yīng)變與時(shí)間關(guān)系,見式(4):

    式中:ε—應(yīng)變;

    t—時(shí)間/d;

    tf—失穩(wěn)時(shí)間/d;

    A、B—常數(shù)。

    Mufundirwa等[22]在式(4)的基礎(chǔ)上令滑坡累計(jì)位移s代替應(yīng)變?chǔ)牛瑫r(shí)式(4)兩側(cè)對時(shí)間t進(jìn)行微分,得到式(5),重新整理后得到式(6)。

    式中:v—滑坡位移速度/(mm·d?1)。

    將式(6)中的v·t視作y,v視作x,tf視作斜率,則式(6)轉(zhuǎn)化為以tf為斜率的線性方程,求失穩(wěn)時(shí)間tf也轉(zhuǎn)變?yōu)榍笤摲匠痰男甭蔾,即求該線性方程的斜率(SLOPE),即SLO預(yù)測模型。該模型實(shí)質(zhì)上是速度倒數(shù)模型在線性條件下的變形,前提假設(shè)為滑坡加速階段的速度倒數(shù)變化符合線性趨勢。與其他預(yù)測模型相比,該模型將預(yù)測失穩(wěn)時(shí)間轉(zhuǎn)化為斜率進(jìn)行求解,不需要進(jìn)行復(fù)雜的擬合和參數(shù)求解,在進(jìn)行突發(fā)型黃土滑坡失穩(wěn)時(shí)間預(yù)測方面具有一定的優(yōu)勢。

    2.4 預(yù)測計(jì)算數(shù)據(jù)選取

    斜率模型是基于滑坡加速蠕變階段提出的,為保證模型適應(yīng)性和預(yù)測結(jié)果可靠性,計(jì)算數(shù)據(jù)應(yīng)從加速變形階段選取。計(jì)算截點(diǎn)的選擇對及時(shí)判斷滑坡危險(xiǎn)性意義重大,計(jì)算截點(diǎn)過早可能會(huì)導(dǎo)致時(shí)間預(yù)報(bào)偏差過大,過晚可能會(huì)導(dǎo)致應(yīng)急處理時(shí)間不足。大量成功預(yù)警的滑坡案例表明,當(dāng)切線角達(dá)到85°時(shí),滑坡進(jìn)入臨滑階段,失穩(wěn)概率極大,此時(shí)發(fā)布預(yù)警預(yù)報(bào)信息能夠有效地指導(dǎo)災(zāi)情應(yīng)急處理工作,降低滑坡造成的經(jīng)濟(jì)損失和人員傷亡。因此,為了更好地評估預(yù)測結(jié)果,本文在進(jìn)行滑坡失穩(wěn)時(shí)間預(yù)測計(jì)算時(shí),以各加速變形階段第一組數(shù)據(jù)為計(jì)算起點(diǎn),以切線角達(dá)到85°為計(jì)算截點(diǎn),采用累計(jì)計(jì)算的方式進(jìn)行預(yù)測計(jì)算(圖3)。

    圖3 累計(jì)計(jì)算示意圖Fig.3 Schematic diagram of cumulative calculation

    但研究發(fā)現(xiàn),由于測量誤差、環(huán)境變化、滑坡局部滑動(dòng)和數(shù)據(jù)處理手段的局限性等多種因素影響,通過監(jiān)測儀器獲得的滑坡變形曲線并不穩(wěn)定,新數(shù)據(jù)的增加會(huì)對預(yù)測結(jié)果產(chǎn)生較大的影響,所以本文在累計(jì)計(jì)算的基礎(chǔ)上提取關(guān)鍵變形階段的數(shù)據(jù)進(jìn)行簡化計(jì)算,只選取滑坡切線角在60 °、65°、70°、75°、80°、85°時(shí)對應(yīng)的監(jiān)測數(shù)據(jù)進(jìn)行累計(jì)計(jì)算,這既能保證計(jì)算數(shù)據(jù)能反映滑坡變形趨勢,又能降低變形波動(dòng)對預(yù)測結(jié)果的影響。

    2.5 預(yù)測結(jié)果評價(jià)

    本文主要采用預(yù)測壽命圖和預(yù)測結(jié)果均方根誤差(RMSE)從定性和定量的角度進(jìn)行模型預(yù)測結(jié)果評價(jià)。

    預(yù)測壽命指的是預(yù)測失穩(wěn)時(shí)間與當(dāng)前預(yù)測時(shí)間節(jié)點(diǎn)的差值,即tfp?tm。預(yù)測壽命圖以預(yù)測壽命為y軸,時(shí)間為x軸,能夠反映預(yù)測結(jié)果提前預(yù)測、滯后預(yù)測、預(yù)測失效和預(yù)測壽命穩(wěn)定性等情況(圖4)。

    圖4 預(yù)測壽命圖(據(jù)文獻(xiàn)[22],有刪改)Fig.4 Predicted life expectancy(modified from Ref.[22])

    采用均方根誤差(RMSE)用來評價(jià)模型預(yù)測精度,RMSE的值越小,預(yù)測精度越高,其表達(dá)式如下:

    式中:yMSE—均方根誤差值;

    xi—滑坡實(shí)際壽命/h;

    N—預(yù)測總數(shù)。

    3 研究結(jié)果

    3.1 滑移崩塌型滑坡

    滑移崩塌型滑坡從發(fā)育到發(fā)生歷時(shí)較短,具有明顯的突發(fā)性特征,但其運(yùn)動(dòng)距離和體積較小?;掳l(fā)育與灌溉水關(guān)系不大,發(fā)育過程主要是坡體在自重和降雨條件下向臨空面發(fā)生不均勻沉降和變形,在坡體后緣發(fā)育有多級裂縫和錯(cuò)臺(tái)。以黨川6#滑坡的預(yù)測結(jié)果為例進(jìn)行分析。該滑坡失穩(wěn)時(shí)間為2019年3月26日凌晨5點(diǎn)。通過改進(jìn)的切線角模型得到滑坡變形曲線(圖5),切線角在凌晨4點(diǎn)20分達(dá)到85°,進(jìn)入加速階段,此時(shí)距滑坡失穩(wěn)還有40 min。

    圖5 黨川6#滑坡變形曲線Fig.5 Deformation curve of the Dangchuan 6# landslide

    黨川6#滑坡各變形階段的預(yù)測壽命如圖6所示,簡化累計(jì)計(jì)算得到的各變形階段預(yù)測壽命波動(dòng)性較小,變化趨勢較平穩(wěn),且與實(shí)際壽命線變化趨勢相近。而累計(jì)計(jì)算得到的各變形階段預(yù)測壽命波動(dòng)性較大,在切線角達(dá)到80°后,預(yù)測壽命變化趨勢才趨于平穩(wěn)。通過計(jì)算兩種方式得到的滑坡各變形階段預(yù)測結(jié)果均方根誤差(表2)可以看出,簡化累計(jì)計(jì)算得到各變形階段預(yù)測結(jié)果均方根誤差均小于0.5 h,預(yù)測精度較高。

    圖6 黨川6#滑坡預(yù)測壽命Fig.6 Predicted life of the Dangchuan 6# landslide

    表2 黨川6#滑坡預(yù)測結(jié)果均方根誤差Table 2 RMSE of the predicted results of the Dangchuan 6#landslide

    通過對比另外兩起滑移崩塌型滑坡(陳家6#和黨川7#)的預(yù)測結(jié)果也可以看出(圖7)。兩種計(jì)算方式得到滑坡各變形階段的預(yù)測壽命變化趨勢相近,累計(jì)計(jì)算得到的各變形階段預(yù)測壽命波動(dòng)性較大,在切線角達(dá)到80°后才趨于穩(wěn)定;而簡化累計(jì)計(jì)算得到的各變形階段預(yù)測壽命波動(dòng)性較小,變化趨勢較平穩(wěn),相應(yīng)的預(yù)測結(jié)果均方根誤差也明顯小于累計(jì)計(jì)算方式,預(yù)測精度較高(表3)。其中,黨川7#滑坡的預(yù)測結(jié)果的均方根誤差較大,但考慮到預(yù)測失穩(wěn)時(shí)間多早于實(shí)際失穩(wěn)時(shí)間,且隨著起算角度的增加,預(yù)測誤差逐漸減小,所以斜率模型在黨川7#滑坡的預(yù)測結(jié)果較“安全”??傮w來看,斜率模型在滑移崩塌型滑坡中的預(yù)測效果較好。

    表3 陳家6#和黨川7#滑坡預(yù)測結(jié)果的均方根誤差Table 3 RMSE of the predicted results of the Chenjia 6# and Dangchuan 7# landslide

    圖7 陳家6#和黨川7#滑坡預(yù)測壽命Fig.7 Predicted life expectancy of Chenjia 6# and Dangchuan 7#

    3.2 靜態(tài)液化型滑坡

    靜態(tài)液化型滑坡整體呈座椅狀,凹向臺(tái)塬內(nèi)部,底部有大量地下水滲出,后壁陡峭?;碌陌l(fā)育一方面是由于該地區(qū)黃土顆粒以粉粒為主,黏粒含量較少,擁有易于產(chǎn)生靜態(tài)液化的顆粒級配;另一方面,該地區(qū)多年大水漫灌的灌溉方式導(dǎo)致地下水位不斷抬升,由于黃土層底部粉質(zhì)黏土層的隔水性,塬邊黃土層底部易形成飽水層,形成軟基效應(yīng),這成為靜態(tài)液化產(chǎn)生的先決條件[31]?;蠡⌒伟疾鄣孛伯a(chǎn)生的匯水作用和坡腳滑坡體堆積造成的地下水雍高也容易導(dǎo)致該類滑坡呈漸進(jìn)后退式發(fā)育。黨川4#滑坡就是一起典型的靜態(tài)液化型滑坡,滑坡變形曲線如圖8所示。滑坡前期變形緩慢且不明顯,位移和變形速率沒有明顯變化,切線角在失穩(wěn)前18 min才達(dá)到85°,突發(fā)性極強(qiáng)。

    圖8 黨川4#滑坡變形曲線Fig.8 Deformation curve of the Dangchuan 4# landslide

    從黨川4#滑坡各變形階段的預(yù)測壽命圖(圖9)可以看出,兩種計(jì)算方式得到的各變形階段預(yù)測壽命變化趨勢相近,且預(yù)測壽命多大于實(shí)際滑坡壽命,說明預(yù)測失穩(wěn)時(shí)間多晚于實(shí)際失穩(wěn)時(shí)間。在切線角達(dá)到80°后,預(yù)測壽命開始向?qū)嶋H壽命線收斂。

    圖9 黨川4#滑坡預(yù)測壽命Fig.9 Predicted life expectancy of the Dangchuan 4# landslide

    對比兩模型得到的各變形階段預(yù)測結(jié)果均方根誤差(表4)可以看出,簡化累計(jì)計(jì)算方式的預(yù)測精度略高于累計(jì)計(jì)算方式,但兩種計(jì)算方式的預(yù)測誤差均較大。從80°切線角起算得到的預(yù)測結(jié)果均方根誤差較小,在切線角達(dá)到80°之后,各變形階段的預(yù)測壽命開始向?qū)嶋H壽命線收斂??傮w來看,斜率模型在靜態(tài)液化型滑坡中的預(yù)測效果較差。

    表4 黨川4#滑坡預(yù)測結(jié)果均方根誤差Table 4 Predicted RMSE of the Dangchuan 4# landslide

    3.3 預(yù)測結(jié)果分析

    3.3.1 速度倒數(shù)變化趨勢的影響

    各變形階段預(yù)測結(jié)果提前或滯后預(yù)測概率的大小與速度倒數(shù)曲線的凹凸性有關(guān)。以黨川4#和黨川7#滑坡為例進(jìn)行分析,在切線角達(dá)到70°之后,黨川4#滑坡各變形階段得到的預(yù)測失穩(wěn)時(shí)間多晚于實(shí)際失穩(wěn)時(shí)間,而黨川7#各變形階段得到的預(yù)測失穩(wěn)時(shí)間均早于實(shí)際失穩(wěn)時(shí)間。對比70°切線角之后兩滑坡各變形階段速度倒數(shù)曲線圖(圖10)可以看出,黨川4#滑坡速度倒數(shù)曲線呈“凸”型,黨川7#滑坡速度倒數(shù)曲線呈“凹”型,由于斜率模型是速度倒數(shù)法在線性條件下的變形,前提假設(shè)為滑坡加速階段的速度倒數(shù)變化符合線性趨勢,即速度倒數(shù)曲線為一條直線,該直線與時(shí)間軸的交點(diǎn)即為失穩(wěn)時(shí)間點(diǎn)[9]。對于黨川4#這類速度倒數(shù)圖像呈“凸”型的滑坡,其前期速度倒數(shù)變化符合線性趨勢,但臨滑階段速度迅速增加,速度倒數(shù)迅速向時(shí)間軸收斂,偏離了前期線性趨勢,前期線性趨勢線與時(shí)間軸的交點(diǎn)多滯后于實(shí)際失穩(wěn)點(diǎn),故滯后預(yù)測概率大。而對于黨川7#這類速度倒數(shù)圖像呈“凹”型的滑坡,其前期速度倒數(shù)的線性趨勢與時(shí)間軸的交點(diǎn)多早于實(shí)際失穩(wěn)點(diǎn),故提前預(yù)測概率大。

    圖10 滑坡速度倒數(shù)曲線圖Fig.10 Inverse velocity of landslides

    3.3.2 滑坡成災(zāi)模式的影響

    滑坡成災(zāi)模式不同是造成不同滑坡預(yù)測結(jié)果差異的根本原因。靜態(tài)液化型滑坡前期變形緩慢且不明顯,主要表現(xiàn)為下伏軟基在重力和地下水作用下向臨空方向緩慢蠕動(dòng),變形量大于1 mm/d,坡體后緣發(fā)育有裂縫,寬度大于3 mm,坡腳局部滲水;隨著變形不斷發(fā)展,上覆黃土裂縫向坡體內(nèi)部發(fā)育,形成潛在滑動(dòng)面,下伏飽和土體結(jié)構(gòu)發(fā)生改變,坡腳有大量地下水滲出,且有顆粒物質(zhì)被帶出,存在地下水局部雍高現(xiàn)象;當(dāng)滑坡進(jìn)入臨滑階段,由變形產(chǎn)生的孔隙水壓力來不及消散,逐漸增大,并形成超孔隙水壓力,底部飽和黃土層產(chǎn)生液化,坡體發(fā)生整體潰散,移動(dòng)速度驟升[32?33]?;鲁蔀?zāi)模式如圖11所示。所以,受靜態(tài)液化的影響,該類滑坡臨滑階段速度驟升,速度倒數(shù)圖像近乎與時(shí)間軸垂直,嚴(yán)重偏離了前期線性變形趨勢,導(dǎo)致其預(yù)測結(jié)果偏差大,模型預(yù)測精度偏低。

    圖11 靜態(tài)液化型滑坡成災(zāi)模式圖(據(jù)文獻(xiàn)[27],有刪改)Fig.11 The disaster-causing mode of loess flow landslides(modified from Ref.[27])

    而滑移崩塌型滑坡的變形過程受地下水影響小,滑坡變形主要是坡體在自重作用下向臨空面發(fā)生移動(dòng),上覆黃土體發(fā)生不均勻沉降和變形,隨著變形的不斷發(fā)展,裂縫繼續(xù)向黃土內(nèi)部擴(kuò)展延伸,形成潛在滑動(dòng)面,錯(cuò)臺(tái)沿平行臨空面方向多級發(fā)育,其成災(zāi)模式如圖12所示。受上覆黃土不均勻沉降和裂縫間歇錯(cuò)動(dòng)的影響,滑移崩塌型滑坡在變形過程中也會(huì)出現(xiàn)速度的波動(dòng),預(yù)測壽命也會(huì)隨速度變化出現(xiàn)相應(yīng)的波動(dòng)。

    圖12 滑移崩塌型滑坡成災(zāi)模式圖(據(jù)文獻(xiàn)[27],有刪改)Fig.12 Disaster-causing mode of loess fall landslides(modified from Ref.[27])

    4 討論

    本文以2019年黑方臺(tái)地區(qū)發(fā)生的4起突發(fā)型黃土滑坡為例,驗(yàn)證了斜率模型在該類滑坡失穩(wěn)時(shí)間預(yù)測方面的可行性,但模型預(yù)測準(zhǔn)確性和適應(yīng)性還需進(jìn)行驗(yàn)證總結(jié)。對于加速階段歷時(shí)極短的黃土滑坡,例如靜態(tài)液化型滑坡,該模型還是有很大的局限性,僅依靠位移一項(xiàng)監(jiān)測指標(biāo)很難實(shí)現(xiàn)滑坡失穩(wěn)臨界點(diǎn)的準(zhǔn)確預(yù)測。有關(guān)這類強(qiáng)突發(fā)性滑坡的預(yù)警預(yù)報(bào)工作應(yīng)在考慮滑坡成因機(jī)制的基礎(chǔ)上,結(jié)合位移、傾角、地下水位等多種指標(biāo)進(jìn)行綜合預(yù)警預(yù)測。本文提出了以切線角為指標(biāo)的簡化累計(jì)計(jì)算,降低了數(shù)據(jù)波動(dòng)對預(yù)測精度的影響,根據(jù)滑坡不同變形階段的變形趨勢選擇更優(yōu)的數(shù)據(jù)處理方法也是今后預(yù)警預(yù)測工作開展的重點(diǎn)。

    5 結(jié)論

    (1)SLO失穩(wěn)時(shí)間預(yù)測模型參數(shù)少,無需進(jìn)行復(fù)雜的擬合,計(jì)算效率高。通過2019年黑方臺(tái)發(fā)生的滑坡實(shí)例驗(yàn)證了該模型在突發(fā)型黃土滑坡失穩(wěn)時(shí)間預(yù)測方面具有一定的可行性。

    (2)由于觀測不確定性和處理方法局限性的影響,將所有的監(jiān)測數(shù)據(jù)進(jìn)行累計(jì)計(jì)算并不一定能取得很好預(yù)測效果。本文在累計(jì)計(jì)算的基礎(chǔ)上提出以切線角為劃分指標(biāo)的簡化累計(jì)計(jì)算,選取指定切線角對應(yīng)的監(jiān)測數(shù)據(jù)進(jìn)行簡化累計(jì)計(jì)算有助于提高模型預(yù)測精度。

    (3)簡化累計(jì)計(jì)算得到的滑坡各變形階段預(yù)測壽命變化趨勢較穩(wěn)定,而累計(jì)計(jì)算方式得到的滑坡各變形階段預(yù)測壽命的變化趨勢波動(dòng)較大,在切線角達(dá)到80°后才趨于穩(wěn)定。簡化累計(jì)計(jì)算得到的滑坡各變形階段預(yù)測結(jié)果均方根誤差小,預(yù)測精度高,從切線角80°起算得到的預(yù)測精度最高。

    (4)預(yù)測結(jié)果提前或滯后預(yù)測概率的大小與速度倒數(shù)曲線的凹凸性有關(guān)。速度倒數(shù)變化趨勢呈“凹”型,預(yù)測時(shí)間多早于實(shí)際失穩(wěn)時(shí)間;速度倒數(shù)變化趨勢呈“凸”型,預(yù)測時(shí)間多晚于實(shí)際失穩(wěn)時(shí)間;速度倒數(shù)變化趨勢符合線性變化趨勢時(shí),模型預(yù)測精度較高。

    (5)滑坡成災(zāi)模式不同是影響各滑坡預(yù)測結(jié)果不同的根本原因,靜態(tài)液化型滑坡能在極短時(shí)間內(nèi)迅速失穩(wěn),速度變化嚴(yán)重偏離了前期變形趨勢,模型預(yù)測結(jié)果偏差大。滑移崩塌型滑坡變形趨勢相對較穩(wěn)定,但也會(huì)因速度波動(dòng)產(chǎn)生預(yù)測偏差。

    猜你喜歡
    倒數(shù)切線黃土
    圓錐曲線的切線方程及其推廣的結(jié)論
    各路創(chuàng)新人才涌向“黃土高坡”
    黃土成金
    只要有信心 黃土變成金
    驚喜倒數(shù)日歷
    切線在手,函數(shù)無憂
    《劉文西:繪不盡是黃土情》
    過圓錐曲線上一點(diǎn)作切線的新方法
    巧用倒數(shù) 求異創(chuàng)新
    讀寫算(下)(2015年11期)2015-11-07 07:21:13
    欧美日韩亚洲综合一区二区三区_| 国产一区二区三区在线臀色熟女 | 少妇 在线观看| 欧美激情久久久久久爽电影 | 午夜福利,免费看| 1024香蕉在线观看| 午夜福利一区二区在线看| 日本五十路高清| 国产日韩欧美亚洲二区| 两个人免费观看高清视频| 午夜免费观看性视频| 亚洲国产精品999| 捣出白浆h1v1| 国产三级黄色录像| 黄色视频,在线免费观看| 午夜福利免费观看在线| 国产成人精品久久二区二区91| 国产成人影院久久av| tocl精华| 亚洲综合色网址| 成人影院久久| 国产成人精品久久二区二区免费| 成人影院久久| 美国免费a级毛片| 久久久久网色| 黄色 视频免费看| 丝瓜视频免费看黄片| 精品亚洲成国产av| 欧美黑人欧美精品刺激| 欧美黑人欧美精品刺激| 少妇裸体淫交视频免费看高清 | 国产成人影院久久av| 如日韩欧美国产精品一区二区三区| 日本wwww免费看| 免费不卡黄色视频| 久久热在线av| 午夜日韩欧美国产| 亚洲精品一卡2卡三卡4卡5卡 | 在线十欧美十亚洲十日本专区| 精品福利观看| 好男人电影高清在线观看| 亚洲成人免费电影在线观看| 国产亚洲精品第一综合不卡| 国产亚洲精品久久久久5区| 亚洲av男天堂| 搡老熟女国产l中国老女人| 久久人人97超碰香蕉20202| 亚洲久久久国产精品| 国产成人影院久久av| 久久久久精品国产欧美久久久 | 日韩欧美一区视频在线观看| 久久免费观看电影| 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 亚洲成av片中文字幕在线观看| 欧美亚洲 丝袜 人妻 在线| 男女高潮啪啪啪动态图| 丰满饥渴人妻一区二区三| 精品视频人人做人人爽| 在线精品无人区一区二区三| 精品乱码久久久久久99久播| 久久狼人影院| 女人久久www免费人成看片| 黄片大片在线免费观看| 久久99热这里只频精品6学生| 妹子高潮喷水视频| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 男女之事视频高清在线观看| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 国产野战对白在线观看| 99国产综合亚洲精品| 免费不卡黄色视频| 下体分泌物呈黄色| 亚洲第一青青草原| 久久精品久久久久久噜噜老黄| 窝窝影院91人妻| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 成人国产av品久久久| 亚洲精品一卡2卡三卡4卡5卡 | 精品欧美一区二区三区在线| 97人妻天天添夜夜摸| 日韩电影二区| 亚洲美女黄色视频免费看| 青青草视频在线视频观看| 国产无遮挡羞羞视频在线观看| 一本一本久久a久久精品综合妖精| 青草久久国产| 99国产精品99久久久久| 啦啦啦免费观看视频1| 51午夜福利影视在线观看| 精品久久蜜臀av无| 久久久久国产精品人妻一区二区| 日本黄色日本黄色录像| 夜夜夜夜夜久久久久| av电影中文网址| 一个人免费看片子| 1024视频免费在线观看| 中文字幕色久视频| 一级黄色大片毛片| 国产一区二区三区综合在线观看| 日本撒尿小便嘘嘘汇集6| 午夜91福利影院| 无遮挡黄片免费观看| 亚洲欧美一区二区三区黑人| 亚洲av电影在线观看一区二区三区| 韩国高清视频一区二区三区| 高清视频免费观看一区二区| 三级毛片av免费| 亚洲精品国产av成人精品| 男人添女人高潮全过程视频| 免费人妻精品一区二区三区视频| 91精品国产国语对白视频| 男人舔女人的私密视频| 啪啪无遮挡十八禁网站| 欧美精品av麻豆av| 999久久久国产精品视频| 一个人免费在线观看的高清视频 | 丝瓜视频免费看黄片| 中文字幕人妻丝袜制服| 午夜久久久在线观看| 黄网站色视频无遮挡免费观看| 在线精品无人区一区二区三| 日本撒尿小便嘘嘘汇集6| 老司机亚洲免费影院| 亚洲精品久久成人aⅴ小说| 麻豆av在线久日| 久久久久久久精品精品| 天天躁狠狠躁夜夜躁狠狠躁| 汤姆久久久久久久影院中文字幕| 精品少妇久久久久久888优播| 亚洲全国av大片| 99精国产麻豆久久婷婷| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 在线看a的网站| www.熟女人妻精品国产| 亚洲专区中文字幕在线| 亚洲一码二码三码区别大吗| 精品第一国产精品| 99国产极品粉嫩在线观看| 狠狠婷婷综合久久久久久88av| 午夜激情久久久久久久| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 欧美日韩亚洲国产一区二区在线观看 | av网站免费在线观看视频| 欧美大码av| 三上悠亚av全集在线观看| 黄色视频,在线免费观看| 日本猛色少妇xxxxx猛交久久| 亚洲三区欧美一区| 视频区欧美日本亚洲| 亚洲成国产人片在线观看| 下体分泌物呈黄色| 男女床上黄色一级片免费看| 精品久久蜜臀av无| 每晚都被弄得嗷嗷叫到高潮| 欧美xxⅹ黑人| 高潮久久久久久久久久久不卡| 成人国产一区最新在线观看| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 韩国精品一区二区三区| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看 | 欧美变态另类bdsm刘玥| 精品熟女少妇八av免费久了| 一边摸一边抽搐一进一出视频| 国产一卡二卡三卡精品| 天堂中文最新版在线下载| 亚洲精品乱久久久久久| 久久久精品区二区三区| 国产一区二区三区在线臀色熟女 | 精品乱码久久久久久99久播| 少妇的丰满在线观看| 精品福利观看| 国产免费现黄频在线看| 久久ye,这里只有精品| 纯流量卡能插随身wifi吗| 亚洲av电影在线进入| 日韩中文字幕视频在线看片| 一边摸一边做爽爽视频免费| 久久女婷五月综合色啪小说| √禁漫天堂资源中文www| 丝袜在线中文字幕| 一级,二级,三级黄色视频| 久久国产亚洲av麻豆专区| 搡老岳熟女国产| 欧美精品一区二区免费开放| 色视频在线一区二区三区| 九色亚洲精品在线播放| 婷婷成人精品国产| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 高清av免费在线| 老汉色∧v一级毛片| 丝袜美腿诱惑在线| 50天的宝宝边吃奶边哭怎么回事| 国产精品香港三级国产av潘金莲| 热99re8久久精品国产| 男人爽女人下面视频在线观看| 日本精品一区二区三区蜜桃| 久久久国产一区二区| 亚洲精品中文字幕在线视频| 欧美在线一区亚洲| 大片免费播放器 马上看| 日韩欧美一区二区三区在线观看 | 国产亚洲欧美在线一区二区| 欧美日韩亚洲国产一区二区在线观看 | 女人精品久久久久毛片| 老汉色∧v一级毛片| 亚洲七黄色美女视频| 午夜免费鲁丝| 日韩精品免费视频一区二区三区| 精品久久久久久电影网| 亚洲熟女毛片儿| 黄色视频,在线免费观看| 中国国产av一级| videos熟女内射| 汤姆久久久久久久影院中文字幕| 欧美激情久久久久久爽电影 | 少妇人妻久久综合中文| 激情视频va一区二区三区| 国产亚洲精品久久久久5区| 黄色a级毛片大全视频| 久久国产精品男人的天堂亚洲| 欧美精品一区二区免费开放| 一级黄色大片毛片| 国产免费福利视频在线观看| 91精品伊人久久大香线蕉| 国产精品国产av在线观看| 一区二区av电影网| 亚洲精品中文字幕在线视频| 久久精品成人免费网站| 搡老熟女国产l中国老女人| 美女高潮喷水抽搐中文字幕| 国产精品欧美亚洲77777| 成人国产一区最新在线观看| 久9热在线精品视频| 岛国在线观看网站| 一个人免费在线观看的高清视频 | h视频一区二区三区| 国产欧美日韩综合在线一区二区| 99热全是精品| 国产精品久久久久久人妻精品电影 | 日韩免费高清中文字幕av| 美女高潮喷水抽搐中文字幕| 国产人伦9x9x在线观看| 男女无遮挡免费网站观看| 18禁黄网站禁片午夜丰满| 日韩 亚洲 欧美在线| 老司机影院毛片| 黄色怎么调成土黄色| 美女中出高潮动态图| 免费在线观看黄色视频的| 黑人巨大精品欧美一区二区mp4| 天天添夜夜摸| 日日夜夜操网爽| 欧美变态另类bdsm刘玥| 午夜激情久久久久久久| 国产精品影院久久| 99精品久久久久人妻精品| 九色亚洲精品在线播放| 欧美在线一区亚洲| 欧美国产精品va在线观看不卡| 性高湖久久久久久久久免费观看| 91老司机精品| 蜜桃在线观看..| 日韩中文字幕视频在线看片| 男女午夜视频在线观看| 国产极品粉嫩免费观看在线| 啦啦啦啦在线视频资源| 视频在线观看一区二区三区| 亚洲一区中文字幕在线| 性色av乱码一区二区三区2| 日日夜夜操网爽| 久久亚洲国产成人精品v| 可以免费在线观看a视频的电影网站| 狠狠精品人妻久久久久久综合| 伊人亚洲综合成人网| 亚洲精品国产av成人精品| 国产在线免费精品| 国精品久久久久久国模美| 热99久久久久精品小说推荐| 久久久精品94久久精品| 久久久精品免费免费高清| 18禁裸乳无遮挡动漫免费视频| 男人舔女人的私密视频| 女警被强在线播放| 久久久久久久国产电影| 亚洲色图 男人天堂 中文字幕| 老熟妇乱子伦视频在线观看 | 99国产精品99久久久久| 国产三级黄色录像| 黑人猛操日本美女一级片| 亚洲欧美激情在线| 日韩人妻精品一区2区三区| 色视频在线一区二区三区| 日韩欧美国产一区二区入口| 亚洲熟女精品中文字幕| 丰满饥渴人妻一区二区三| 一级片免费观看大全| 欧美性长视频在线观看| 国内毛片毛片毛片毛片毛片| 国产成人a∨麻豆精品| 欧美另类亚洲清纯唯美| 国产成人免费观看mmmm| 18禁观看日本| 久久久国产成人免费| 国产亚洲一区二区精品| 两性夫妻黄色片| 亚洲五月色婷婷综合| 男女免费视频国产| 国产在视频线精品| 99精国产麻豆久久婷婷| 亚洲中文日韩欧美视频| 性高湖久久久久久久久免费观看| 美女福利国产在线| 欧美日本中文国产一区发布| 久久国产精品人妻蜜桃| a级片在线免费高清观看视频| 久久久国产一区二区| 国产一区有黄有色的免费视频| 国产成人一区二区三区免费视频网站| 999久久久国产精品视频| 男人操女人黄网站| 蜜桃国产av成人99| 香蕉国产在线看| 亚洲专区国产一区二区| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品电影小说| 久久中文字幕一级| 免费在线观看视频国产中文字幕亚洲 | 可以免费在线观看a视频的电影网站| 国产精品影院久久| 国产区一区二久久| 丝瓜视频免费看黄片| 亚洲人成电影免费在线| 久久久久久人人人人人| 人人妻人人澡人人看| 老汉色∧v一级毛片| 成年人黄色毛片网站| 日韩中文字幕欧美一区二区| 日本av免费视频播放| 日韩一区二区三区影片| 欧美精品一区二区免费开放| 亚洲国产欧美日韩在线播放| 久久香蕉激情| 自线自在国产av| 别揉我奶头~嗯~啊~动态视频 | 狂野欧美激情性xxxx| 丝袜喷水一区| 狂野欧美激情性xxxx| 亚洲国产毛片av蜜桃av| 大码成人一级视频| 午夜免费观看性视频| 法律面前人人平等表现在哪些方面 | 精品国产乱码久久久久久小说| 欧美日韩亚洲高清精品| 在线观看免费午夜福利视频| 成年av动漫网址| 操出白浆在线播放| 精品人妻1区二区| 一区二区av电影网| 在线观看人妻少妇| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人系列免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 巨乳人妻的诱惑在线观看| 国产男人的电影天堂91| avwww免费| 日本猛色少妇xxxxx猛交久久| 午夜免费观看性视频| 美女扒开内裤让男人捅视频| 精品视频人人做人人爽| 亚洲精品国产一区二区精华液| 精品一品国产午夜福利视频| 男女国产视频网站| 亚洲成av片中文字幕在线观看| 精品国产乱子伦一区二区三区 | 无限看片的www在线观看| 国产精品av久久久久免费| 一区二区三区乱码不卡18| 黑丝袜美女国产一区| 国产成人免费无遮挡视频| www.av在线官网国产| 真人做人爱边吃奶动态| 日韩视频一区二区在线观看| 亚洲精品国产av成人精品| 欧美xxⅹ黑人| 亚洲欧美精品综合一区二区三区| 国产精品成人在线| 亚洲国产欧美日韩在线播放| 日韩制服骚丝袜av| 波多野结衣av一区二区av| 大码成人一级视频| 国产免费av片在线观看野外av| 制服人妻中文乱码| 日本黄色日本黄色录像| netflix在线观看网站| 免费高清在线观看日韩| 欧美精品啪啪一区二区三区 | 一区二区三区乱码不卡18| 免费观看av网站的网址| 男人爽女人下面视频在线观看| av超薄肉色丝袜交足视频| 麻豆乱淫一区二区| 99香蕉大伊视频| 午夜精品国产一区二区电影| 一进一出抽搐动态| 国产亚洲精品第一综合不卡| 免费久久久久久久精品成人欧美视频| 欧美激情 高清一区二区三区| 91国产中文字幕| 国产精品一二三区在线看| 欧美日韩视频精品一区| 91字幕亚洲| 日韩一卡2卡3卡4卡2021年| 免费黄频网站在线观看国产| 国产亚洲午夜精品一区二区久久| 宅男免费午夜| av在线播放精品| 五月开心婷婷网| 人妻人人澡人人爽人人| 亚洲av欧美aⅴ国产| 久久久久国产精品人妻一区二区| 国产欧美日韩一区二区精品| 亚洲av成人不卡在线观看播放网 | 国产精品国产av在线观看| 国产精品一区二区免费欧美 | 男人爽女人下面视频在线观看| 久久天躁狠狠躁夜夜2o2o| 国产日韩欧美视频二区| 日本一区二区免费在线视频| videos熟女内射| 欧美另类亚洲清纯唯美| 欧美黄色片欧美黄色片| 热re99久久国产66热| 女人久久www免费人成看片| 一本久久精品| 母亲3免费完整高清在线观看| 精品少妇一区二区三区视频日本电影| av福利片在线| 在线观看人妻少妇| 国产片内射在线| 捣出白浆h1v1| 亚洲欧美清纯卡通| 久久久久国内视频| 亚洲专区中文字幕在线| 精品国产乱码久久久久久小说| av在线app专区| 国产男女内射视频| 亚洲av美国av| 国产成人精品久久二区二区免费| 国产成人精品无人区| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 国产精品影院久久| 欧美中文综合在线视频| 18禁国产床啪视频网站| 久久久久久久精品精品| 精品国产乱码久久久久久男人| 亚洲成人国产一区在线观看| 精品国内亚洲2022精品成人 | 少妇精品久久久久久久| 免费av中文字幕在线| 深夜精品福利| 黄色怎么调成土黄色| 黑丝袜美女国产一区| 亚洲国产精品成人久久小说| 国产精品成人在线| av在线app专区| 国产精品欧美亚洲77777| 精品人妻熟女毛片av久久网站| 咕卡用的链子| 免费看十八禁软件| 伊人久久大香线蕉亚洲五| 一本色道久久久久久精品综合| 午夜免费观看性视频| 一级毛片女人18水好多| 天堂中文最新版在线下载| 精品少妇一区二区三区视频日本电影| 人人妻人人澡人人看| 国产真人三级小视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 精品乱码久久久久久99久播| 久久精品国产综合久久久| 老鸭窝网址在线观看| 12—13女人毛片做爰片一| 夜夜夜夜夜久久久久| 国产免费av片在线观看野外av| 亚洲 欧美一区二区三区| 国产伦人伦偷精品视频| 亚洲精品成人av观看孕妇| 免费在线观看黄色视频的| 啦啦啦免费观看视频1| 久久久久久久久免费视频了| 精品国产乱码久久久久久男人| 黑人巨大精品欧美一区二区mp4| 视频区图区小说| 国产精品香港三级国产av潘金莲| 国产99久久九九免费精品| 性色av乱码一区二区三区2| 韩国高清视频一区二区三区| 十分钟在线观看高清视频www| 久久天堂一区二区三区四区| 午夜日韩欧美国产| 一级,二级,三级黄色视频| 狠狠精品人妻久久久久久综合| 国产亚洲欧美精品永久| 在线天堂中文资源库| www日本在线高清视频| 亚洲五月婷婷丁香| 亚洲成人免费av在线播放| 成人影院久久| 美女脱内裤让男人舔精品视频| 国产精品99久久99久久久不卡| 国产欧美日韩综合在线一区二区| 精品一品国产午夜福利视频| 国产免费视频播放在线视频| 日韩人妻精品一区2区三区| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲国产一区二区在线观看 | 丰满少妇做爰视频| 国产成人影院久久av| 日本五十路高清| 国产高清国产精品国产三级| 色94色欧美一区二区| 国产精品久久久久久精品电影小说| 中文字幕人妻熟女乱码| 欧美97在线视频| 国产av国产精品国产| 亚洲精品中文字幕一二三四区 | 久久久久国内视频| 国产精品国产三级国产专区5o| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 成人影院久久| 捣出白浆h1v1| 大片免费播放器 马上看| 男女床上黄色一级片免费看| 性色av乱码一区二区三区2| 9191精品国产免费久久| 美女脱内裤让男人舔精品视频| 亚洲中文日韩欧美视频| 亚洲性夜色夜夜综合| 91精品伊人久久大香线蕉| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 岛国毛片在线播放| 午夜福利视频精品| 俄罗斯特黄特色一大片| 在线 av 中文字幕| 美女大奶头黄色视频| 国产日韩欧美视频二区| 青草久久国产| 999久久久国产精品视频| 操美女的视频在线观看| 国产精品麻豆人妻色哟哟久久| 美女扒开内裤让男人捅视频| 国产欧美日韩一区二区三 | 我要看黄色一级片免费的| 人成视频在线观看免费观看| 97精品久久久久久久久久精品| 91国产中文字幕| 自线自在国产av| 可以免费在线观看a视频的电影网站| avwww免费| 丰满迷人的少妇在线观看| 久久久久精品人妻al黑| 亚洲中文日韩欧美视频| 国产亚洲精品久久久久5区| 秋霞在线观看毛片| 欧美日韩福利视频一区二区| 国产日韩欧美亚洲二区| 多毛熟女@视频| 国产免费现黄频在线看| 一区二区三区激情视频| 可以免费在线观看a视频的电影网站| 啦啦啦 在线观看视频| 亚洲专区国产一区二区| 首页视频小说图片口味搜索| 日韩,欧美,国产一区二区三区| 手机成人av网站| 少妇人妻久久综合中文| 99re6热这里在线精品视频| 9热在线视频观看99| 国产成人av教育| 亚洲,欧美精品.| 老司机深夜福利视频在线观看 | 精品一区二区三卡| 国产免费视频播放在线视频| 国产成人影院久久av| 日本欧美视频一区| 99国产精品99久久久久| 国产精品偷伦视频观看了| 男女无遮挡免费网站观看| 久久精品国产综合久久久| 后天国语完整版免费观看| 80岁老熟妇乱子伦牲交| 我要看黄色一级片免费的| 99热国产这里只有精品6| 国产精品国产av在线观看| 不卡av一区二区三区| 热99re8久久精品国产| 男女床上黄色一级片免费看| 精品福利观看| 18禁裸乳无遮挡动漫免费视频| 国产一区二区激情短视频 | 两人在一起打扑克的视频|