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

    考慮巖體劣化的三峽庫(kù)區(qū)順層巖質(zhì)邊坡破壞分析

    2022-10-06 03:12:42盼,黃林,秦臻,閆國(guó)強(qiáng),白
    人民長(zhǎng)江 2022年9期
    關(guān)鍵詞:順層青石滑體

    秦 盼 盼,黃 波 林,秦 臻,閆 國(guó) 強(qiáng),白 林 豐

    (1.三峽大學(xué) 防災(zāi)減災(zāi)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 宜昌 443002; 2.三峽大學(xué) 湖北長(zhǎng)江三峽滑坡國(guó)家野外科學(xué)觀測(cè)研究站,湖北 宜昌 443002; 3.中國(guó)地質(zhì)大學(xué)(武漢) 工程學(xué)院,湖北 武漢 430074)

    0 引 言

    三峽庫(kù)區(qū)廣泛發(fā)育順層基巖滑坡,對(duì)沿岸居民和長(zhǎng)江航道帶來(lái)巨大的威脅。例如,2003年7月13日發(fā)生的秭歸縣千將坪滑坡,體積約為2 400萬(wàn)m3,產(chǎn)生了約39 m高的涌浪,造成24人死亡或失蹤[1-2]。同時(shí),順層巖質(zhì)滑坡特殊的地質(zhì)結(jié)構(gòu),導(dǎo)致水庫(kù)蓄水對(duì)順層滑坡的穩(wěn)定性影響強(qiáng)烈。

    Yin等認(rèn)為三峽水庫(kù)誘發(fā)滑坡可分為3個(gè)階段,不同階段其主要誘發(fā)因素不同,分別是水位變化疊加降雨和巖體劣化[3]。從順層滑坡典型案例來(lái)看,千將坪滑坡是2003年三峽庫(kù)區(qū)首次蓄水達(dá)到135.00 m誘發(fā)形成的,杉樹(shù)槽滑坡是在2014年9月1日水位大幅上漲和暴雨期間發(fā)生的[4],卡門子灣滑坡是2019年12月10日發(fā)生,由阻滑段巖體強(qiáng)度劣化所致[3]。

    在誘發(fā)滑坡的前兩個(gè)階段,對(duì)三峽庫(kù)區(qū)順層巖質(zhì)岸坡研究較多,取得了大量認(rèn)識(shí)。例如,在水位變動(dòng)下順層滑坡研究方面,劉才華等研究指出庫(kù)水位上升會(huì)使巖土體水致弱化,滑帶抗剪強(qiáng)度急劇降低,坡體抗滑力降低,影響岸坡穩(wěn)定性[5];鄧華鋒等研究指出在長(zhǎng)期水-巖循環(huán)作用下的力學(xué)參數(shù)弱化,會(huì)使穩(wěn)定滑坡向不穩(wěn)定方向轉(zhuǎn)變[6];Jian等研究發(fā)現(xiàn)千將坪滑坡失穩(wěn)的主要誘發(fā)因素是庫(kù)水位上升和持續(xù)降雨[7];代貞偉的研究表明坡體軟硬相間順層單斜結(jié)構(gòu)是藕塘滑坡形成的內(nèi)在因素,降雨及水庫(kù)蓄水是導(dǎo)致藕塘滑坡復(fù)活的主要誘因[8]。但關(guān)于第3個(gè)階段庫(kù)水位長(zhǎng)期波動(dòng)條件對(duì)順層岸坡變形破壞的影響研究較少。Tang等通過(guò)對(duì)黃土坡滑坡進(jìn)行長(zhǎng)期現(xiàn)場(chǎng)原位監(jiān)測(cè),發(fā)現(xiàn)坡體位移曲線中的“跳躍點(diǎn)”與暴雨和庫(kù)區(qū)水位快速下降的時(shí)間相吻合[9]。黃波林等研究發(fā)現(xiàn)巖溶岸坡巖體劣化早期,以破壞淺表、岸坡后退為主,引起的主要破壞模式為崩塌、滑移和傾倒3種類型[10]。劉廣寧等也發(fā)現(xiàn)受庫(kù)水位波動(dòng)影響,坡體前緣受江水侵蝕,坡腳部位的軟弱夾層向內(nèi)凹陷,巖體會(huì)發(fā)生滑移、錯(cuò)動(dòng)、坍塌等破壞模式[11]。

    2008年以來(lái),三峽庫(kù)區(qū)已經(jīng)經(jīng)歷了12個(gè)周期的145.00~175.00 m水位變動(dòng)。三峽庫(kù)區(qū)巫峽段內(nèi)青石-抱龍一帶順層消落帶巖體強(qiáng)烈劣化。本文以青石6號(hào)坡為研究對(duì)象,利用連續(xù)-非連續(xù)數(shù)值模擬,分析斜坡裂縫擴(kuò)展、變形破壞機(jī)理和破壞運(yùn)動(dòng)過(guò)程。然后結(jié)合室內(nèi)飽和-風(fēng)干循環(huán)試驗(yàn)結(jié)果,討論了巖體劣化對(duì)岸坡變形破壞過(guò)程所產(chǎn)生的影響。研究結(jié)果可為三峽庫(kù)區(qū)順層巖質(zhì)岸坡的防災(zāi)減災(zāi)提供技術(shù)支撐。

    1 青石6號(hào)坡工程地質(zhì)概況

    青石6號(hào)坡位于三峽庫(kù)區(qū)巫峽段右岸,距上游巫山縣城16 km。坡體在高程170 m處,地形起伏較大,將坡體分為上下兩部分。坡體上部呈馬鞍狀,凸起端長(zhǎng)度達(dá)218 m,坡體下部坡面角度逐漸轉(zhuǎn)緩,形成陡緩平臺(tái)。坡體兩側(cè)沖溝為上下游邊界,前緣入江,受庫(kù)水位波動(dòng)影響(見(jiàn)圖1)。水位為145.00 m時(shí),斜坡縱長(zhǎng)約410 m,寬約為170 m,總體坡向345°。坡體表層出露三疊系下統(tǒng)嘉陵江組3段(T1j3)白云巖夾泥質(zhì)灰?guī)r,巖層產(chǎn)狀為340°~350°∠40°~50°(見(jiàn)圖2)。坡體由多個(gè)巖層組成,且?guī)r層傾角略大于坡角,為典型的順層巖質(zhì)岸坡。研究區(qū)位于楠木園背斜北翼、官渡口向斜南翼,地層穩(wěn)定,無(wú)活動(dòng)斷裂通過(guò)。

    岸坡消落帶范圍內(nèi)巖體發(fā)育有兩組產(chǎn)狀分別為34°∠70°和270°∠80°的優(yōu)勢(shì)結(jié)構(gòu)面,平均間距約為2.8 m,僅在坡體表層發(fā)育。巖體結(jié)構(gòu)面與層面相互切割坡面巖體,形成碎塊狀板狀六面體巖體。巖層面間以鈣質(zhì)膠結(jié)物填充,少量巖層面處于張開(kāi)狀態(tài)。在庫(kù)水長(zhǎng)期淘蝕作用下,巖層節(jié)理的張開(kāi)程度逐漸擴(kuò)大。在上覆巖體的重力擠壓作用下,消落帶范圍巖體的破碎程度不斷增加,巖體強(qiáng)度不斷下降。

    經(jīng)長(zhǎng)期變形演化后,坡體上出現(xiàn)明顯的變形破壞跡象。坡體高程165 m處,間續(xù)發(fā)育有多處波狀彎曲現(xiàn)象,從坡體上游邊界間斷延伸至坡體下游邊界。在坡體兩側(cè)沖溝內(nèi)觀察,發(fā)現(xiàn)波狀彎曲由表層向內(nèi)部延伸有一定深度(見(jiàn)圖3)。波狀彎曲撓曲段由于彎折角度過(guò)大,新生大量縱向裂縫。在對(duì)青石岸坡范圍內(nèi)的21個(gè)斜坡進(jìn)行調(diào)查后發(fā)現(xiàn),在互不相鄰的9個(gè)斜坡上,均發(fā)育有彎曲程度不同的波狀彎曲。青石岸坡上所出現(xiàn)空間分布不連續(xù)、彎曲程度大小不一的波狀彎曲現(xiàn)象說(shuō)明,在青石岸坡上出現(xiàn)的波狀現(xiàn)象并不是由地質(zhì)營(yíng)力所致,而是由坡體上部重力作用所致。

    2 青石6號(hào)坡變形破壞全過(guò)程分析

    受自重和庫(kù)水位波動(dòng)的影響,青石6號(hào)坡坡面巖體較為破碎,質(zhì)量較差,具有明顯變形跡象。不同破壞模式下,坡體運(yùn)動(dòng)方式不同,產(chǎn)生的災(zāi)害大小也有所差異。本文選用Munjiza提出的FEM/DEM耦合分析方法,對(duì)青石6號(hào)坡進(jìn)行數(shù)值計(jì)算,預(yù)測(cè)其變形破壞全過(guò)程。

    2.1 FEM/DEM原理

    本文采用文獻(xiàn)[12]中的方法,該方法結(jié)合了FEM(有限元)和DEM(離散元)優(yōu)點(diǎn),能夠模擬出層狀巖體中裂縫擴(kuò)展和破碎過(guò)程,得到巖質(zhì)邊坡漸進(jìn)式損傷和破壞過(guò)程[12-13]。

    FEM/DEM耦合分析的基本思路是將整個(gè)坡體劃分為三角形有限元網(wǎng)格,網(wǎng)格之間由4個(gè)節(jié)點(diǎn)無(wú)厚度的節(jié)理單元相連[12](見(jiàn)圖4),形成完整的坡體結(jié)構(gòu)。通過(guò)在有限元網(wǎng)格中插入離散元網(wǎng)格表示滑體區(qū)域,允許離散元網(wǎng)格發(fā)生轉(zhuǎn)動(dòng)與位移。模型中的破壞準(zhǔn)則包括抗剪強(qiáng)度、抗拉強(qiáng)度和斷裂能。

    對(duì)坡體模擬過(guò)程中,受重力作用影響,連接離散元網(wǎng)格的節(jié)理單元發(fā)生錯(cuò)動(dòng)時(shí)滑體僅沿著滑動(dòng)面產(chǎn)生蠕動(dòng)變形;進(jìn)一步演化時(shí),節(jié)理單元發(fā)生分離,離散元網(wǎng)格相互張開(kāi),滑體開(kāi)始出現(xiàn)破裂。通過(guò)對(duì)坡體連續(xù)-非連續(xù)的分析[13-14],實(shí)現(xiàn)了變形破壞全過(guò)程模擬。上述模型的主要特點(diǎn)為既可對(duì)坡體初始階段所發(fā)生的小變形進(jìn)行模擬,也可對(duì)坡體在破壞過(guò)程中所發(fā)生的大規(guī)模破壞、運(yùn)動(dòng)巖體相互分離進(jìn)行模擬,從而反映出邊坡巖體變形破裂這一復(fù)雜巖體力學(xué)行為。

    在巖質(zhì)邊坡分析方面,F(xiàn)EM/DEM耦合分析方法可以模擬坡體中裂縫擴(kuò)展和巖體破碎過(guò)程,完整展示坡體的整個(gè)運(yùn)動(dòng)過(guò)程,對(duì)分析巖質(zhì)邊坡具有獨(dú)特的優(yōu)勢(shì)[15-16]。國(guó)內(nèi)外學(xué)者在巖體裂紋擴(kuò)展[17-18]、洞室開(kāi)挖[19]、巖質(zhì)邊坡[20]、危巖體分析[21]等方面均取得了較好成果,驗(yàn)證了FEM/DEM方法的可行性和普適性。但由于FEM/DEM后處理程序缺少對(duì)網(wǎng)格位移量的追蹤,不能得到坡體運(yùn)動(dòng)速度和位移量,無(wú)法對(duì)坡體變形情況定量分析。

    2.2 模型建立

    青石6號(hào)坡數(shù)值計(jì)算模型長(zhǎng)377 m,高297 m,由基巖巖體與滑體構(gòu)成(見(jiàn)圖5)。整個(gè)計(jì)算模型由2 000個(gè)均勻大小的計(jì)算網(wǎng)格組成,網(wǎng)格之間由節(jié)點(diǎn)相連。利用有限元網(wǎng)格單元模擬基巖巖體,離散元網(wǎng)格單元模擬滑體區(qū)域,節(jié)理單元模擬層面。從坡體結(jié)構(gòu)來(lái)看,青石6號(hào)坡坡體中部出現(xiàn)了明顯的凸起段,在上部巖體的重力作用下,巖層順著層面滑動(dòng)。由于坡體前緣坡面角度逐漸轉(zhuǎn)緩,巖層滑移受阻,使巖層發(fā)生波狀彎曲。故將坡體凸起段劃分為滑體區(qū)域,其余區(qū)域在模擬過(guò)程中不發(fā)生變化,將剩余區(qū)域劃分為基巖巖體。

    計(jì)算模型中所需的巖體強(qiáng)度參數(shù)如表1所列,各部分取值先結(jié)合現(xiàn)場(chǎng)實(shí)測(cè)和經(jīng)驗(yàn)進(jìn)行綜合取值,再通過(guò)對(duì)比坡體實(shí)際變形破壞情況和數(shù)值模擬結(jié)果進(jìn)行不斷調(diào)整優(yōu)化。表中基巖巖體和滑體區(qū)域主要由灰?guī)r巖體組成,強(qiáng)度較高?;瑒?dòng)面位于巖層層面之間,強(qiáng)度主要由填充其中的泥質(zhì)填充物維持。經(jīng)長(zhǎng)期庫(kù)水掏蝕,泥質(zhì)填充物之間的膠結(jié)能力越來(lái)越差,松散程度不斷增大,致使滑動(dòng)面的巖土體強(qiáng)度大大降低。所取得滑動(dòng)面內(nèi)摩擦角遠(yuǎn)低于滑動(dòng)面傾角,使得上覆巖體具備沿滑移面向下滑移的條件。

    表1 計(jì)算模型中巖體強(qiáng)度參數(shù)Tab.1 Rock strength parameters of calculation model

    在計(jì)算過(guò)程中,使用張拉破壞斷裂能和剪切破壞斷裂能對(duì)巖質(zhì)邊坡破壞前和破壞過(guò)程中所出現(xiàn)的巖體破裂、損傷和軟化現(xiàn)象進(jìn)行模擬[22],通過(guò)不斷計(jì)算節(jié)點(diǎn)位置得到巖質(zhì)邊坡的運(yùn)動(dòng)后狀態(tài)。對(duì)模型底部和左側(cè)邊界節(jié)點(diǎn)/單元施加固定約束,固定坡體位置。對(duì)坡體內(nèi)部節(jié)點(diǎn)/單元施加重力,模擬坡體在重力作用下變形破壞過(guò)程。計(jì)算模型中的最大不平衡力與初始最大不平衡力比值低于10-5或離散網(wǎng)格的物理位移趨近于一個(gè)常數(shù)不再變化便視為模型已達(dá)到應(yīng)力平衡,單元網(wǎng)格不再發(fā)生變化[23]。

    2.3 數(shù)值模擬結(jié)果分析

    本文中的工況僅考慮坡體在重力作用下的變形演化過(guò)程,并未涉及其他因素?;跀?shù)值模擬結(jié)果,可以將其變形演化過(guò)程劃分為裂縫新生、擴(kuò)展、彎曲-隆起、裂縫貫通、大規(guī)模滑移等5個(gè)階段。初始階段,坡體結(jié)構(gòu)完整,無(wú)網(wǎng)格單元張開(kāi)(見(jiàn)圖6),以連續(xù)分析為主?;鶐r巖體中節(jié)理單元強(qiáng)度較大,在接下來(lái)的演化過(guò)程中處于穩(wěn)定狀態(tài),計(jì)算全過(guò)程為連續(xù)分析。

    裂縫新生階段:坡體中部巖體(主動(dòng)傳力區(qū))在重力作用下內(nèi)部應(yīng)力逐步調(diào)整。坡體中下部巖體(被動(dòng)擠壓區(qū))發(fā)生蠕滑,滑體發(fā)生剪張破壞,剪力向臨空面方向作用,形成張拉裂縫(見(jiàn)圖7)。此階段,僅坡體中下部巖體(被動(dòng)擠壓區(qū))的節(jié)理單元斷裂,少量網(wǎng)格相互分離,滑體以連續(xù)分析為主,占坡體變形破壞總時(shí)步的33%。

    裂縫擴(kuò)展階段:在坡體中部巖體重力的持續(xù)作用下,坡體中應(yīng)力不斷調(diào)整,坡體中上部巖體(主動(dòng)傳力區(qū))繼續(xù)蠕變,中下部巖體(被動(dòng)擠壓區(qū))中的新生裂縫由內(nèi)部向兩端不斷張開(kāi)擴(kuò)展,向上延伸長(zhǎng)度增加(見(jiàn)圖8)。此階段,被動(dòng)擠壓區(qū)中節(jié)理單元斷裂數(shù)量不斷增加,使得三角形單元張開(kāi)度逐漸增加?;w下部開(kāi)始轉(zhuǎn)為非連續(xù)分析,上部仍保持連續(xù)介質(zhì)。這一階段占坡體變形破壞總時(shí)步的23%。

    彎曲-隆起階段:隨著作用時(shí)間的增加,被動(dòng)擠壓區(qū)垂直于層面的壓應(yīng)力遠(yuǎn)小于順層壓應(yīng)力,由于兩者壓力差過(guò)大,致使坡面輕微彎曲隆起。此階段中,在坡體內(nèi)部形成明顯的應(yīng)力集中區(qū)[24]和“X”形節(jié)理。彎曲-隆起區(qū)域位于坡面彎曲起伏轉(zhuǎn)折段,使得滑移段巖體下滑受阻,有利于此現(xiàn)象出現(xiàn)。隨著巖層彎曲隆起程度不斷增大,當(dāng)彎折角度過(guò)大時(shí),彎曲部位開(kāi)始出現(xiàn)縱向裂縫,巖體逐漸碎裂化發(fā)生松動(dòng),在彎曲部位前緣出現(xiàn)巖體崩落(見(jiàn)圖9)。此階段,被動(dòng)擠壓區(qū)中節(jié)理單元的斷裂數(shù)量持續(xù)增加,滑體與基巖巖體發(fā)生錯(cuò)動(dòng)進(jìn)而發(fā)生節(jié)理單元變形,滑體以非連續(xù)分析為主,占坡體變形破壞總時(shí)步的33%。

    裂縫貫通階段:隨著坡體中部巖體(主動(dòng)傳力區(qū))蠕滑位移的持續(xù)增加,對(duì)坡體中下部巖體(被動(dòng)擠壓區(qū))的擠壓程度不斷增大,彎曲-隆起程度增加。隆起部位內(nèi)的巖體裂縫相互貫通形成完整剪出面,坡體向臨空面已產(chǎn)生一定位移(見(jiàn)圖10)。此階段,滑體中三角形網(wǎng)格的張開(kāi)度不斷擴(kuò)大,使得相鄰張開(kāi)區(qū)域相連,滑體以非連續(xù)分析為主,占變形破壞總時(shí)步的11%。

    滑移破壞階段:坡體中下部巖體(被動(dòng)擠壓區(qū))阻抗能力達(dá)到極限,滑移面貫通,形成完整剪出口。坡體順著坡面方向發(fā)生“滑移-剪斷”引起大規(guī)模滑移,最終堆積在坡腳(見(jiàn)圖11)。此階段,滑體中的三角形單元相互分離逐漸碎片化,以非連續(xù)分析為主。

    從以上分析可知,青石6號(hào)坡裂縫主要是由于巖層相互擠壓形成?;w對(duì)阻滑端巖體形成擠壓,巖層發(fā)生彎曲-隆起,隨著彎曲程度的增加,巖層破碎程度增大,逐漸形成完整的滑移面,沿著剪出口發(fā)生大規(guī)?;?。青石6號(hào)坡數(shù)值分析過(guò)程中所出現(xiàn)的滑移-彎曲現(xiàn)象,符合現(xiàn)階段坡體表面所出現(xiàn)的變形跡象,可以用于預(yù)測(cè)斜坡的破壞過(guò)程。對(duì)于青石6號(hào)坡這類庫(kù)水型岸坡,水位變動(dòng)對(duì)消落帶巖體所產(chǎn)生的損傷劣化作用,在岸坡的變形演化過(guò)程中起到了極為重要的作用[25]。如進(jìn)一步結(jié)合庫(kù)水位變動(dòng)對(duì)消落帶巖體強(qiáng)度和結(jié)構(gòu)所產(chǎn)生的不良影響,討論不同劣化情況下岸坡變形破壞模式如何變化,所得結(jié)果會(huì)更貼近庫(kù)區(qū)的實(shí)際工況,適用范圍更加廣泛。

    3 庫(kù)水位長(zhǎng)期變動(dòng)對(duì)巖體劣化的影響

    3.1 庫(kù)水對(duì)岸坡的損傷劣化作用

    自三峽水庫(kù)開(kāi)始正式運(yùn)行后,庫(kù)水位經(jīng)周期性調(diào)度后,消落帶巖體長(zhǎng)期處于“飽和-風(fēng)干”條件,水巖相互作用會(huì)造成巖體強(qiáng)度劣化。劉新榮等通過(guò)干濕循環(huán)試驗(yàn)?zāi)M砂巖在庫(kù)水位漲落情況下,砂巖抗剪強(qiáng)度劣化規(guī)律[26-27]。鄧華鋒等發(fā)現(xiàn)經(jīng)循環(huán)加載作用后,內(nèi)部損傷的砂巖巖樣在浸泡-風(fēng)干作用下抗剪、抗壓強(qiáng)度劣化效應(yīng)明顯,具有明顯的時(shí)間效應(yīng)和非均勻性[28]。以上說(shuō)明庫(kù)水位周期性升降會(huì)使巖體強(qiáng)度逐步劣化,但是對(duì)岸坡穩(wěn)定性和變形破壞過(guò)程會(huì)產(chǎn)生何種影響,仍然值得探討。

    本文選取青石-抱龍段4個(gè)不同區(qū)域的灰?guī)r(T1j3)巖樣,按照試驗(yàn)規(guī)程要求進(jìn)行“浸泡-風(fēng)干”循環(huán)試驗(yàn)[29]。其中試樣1,2位于岸坡消落帶之上,試樣3,4位于消落帶范圍內(nèi),已經(jīng)歷過(guò)多次水位變動(dòng),由此試樣3,4的初始單軸抗壓強(qiáng)度低于試樣1,2的初始單軸抗壓強(qiáng)度。每循環(huán)5個(gè)周期對(duì)巖樣進(jìn)行力學(xué)強(qiáng)度測(cè)定,測(cè)定經(jīng)過(guò)0,5,10,15,20,25,30,35,40,45,50次循環(huán)試驗(yàn)后灰?guī)r(T1j3)巖樣的飽和單軸抗壓強(qiáng)度,得到在不同循環(huán)周期時(shí)灰?guī)r單軸抗壓強(qiáng)度下降曲線(見(jiàn)圖12)。經(jīng)過(guò)50次循環(huán)試驗(yàn)后,4個(gè)試樣飽和單軸抗壓強(qiáng)度的總下降率約為18%~26%,單次循環(huán)平均下降率約為0.37%~0.51%(見(jiàn)圖13),總體相差不大。

    受尺寸效應(yīng)的影響,室內(nèi)巖樣的下降速率并不能表征岸坡巖體下降速率。閆國(guó)強(qiáng)等根據(jù)灰?guī)r的劣化損傷特征,引入核衰變劣化控制方程[30]:

    N(t)=N0e-λt

    (1)

    式中:N0為初始物理強(qiáng)度參數(shù),λ為劣化常數(shù),t為循環(huán)周期。

    式(1)可以較好地描述灰?guī)r巖體物理力學(xué)參數(shù)的衰減過(guò)程。本文利用核衰變劣化控制方程對(duì)上述結(jié)果進(jìn)行修正,得到在不同循環(huán)次數(shù)下灰?guī)r的飽和單軸抗壓強(qiáng)度,如表2所列。

    由表2可知,4個(gè)試樣飽和單軸抗壓強(qiáng)度的總下降率近似相等,在50次循環(huán)周期后約為29.88%。此下降率明顯大于室內(nèi)干濕循環(huán)作用下的下降率結(jié)果,考慮到了巖體的實(shí)際物理力學(xué)屬性,更加貼近實(shí)際結(jié)果。

    表2 采用劣化控制方程得到的灰?guī)r飽和單軸抗壓強(qiáng)度Tab.2 Saturated uniaxial compressive strength of limestone by deterioration control equation

    由此將消落帶巖體(坡面高程145~175 m內(nèi))單個(gè)循環(huán)周期強(qiáng)度下降率近似取值為采用劣化控制方程下灰?guī)r飽和單軸抗壓強(qiáng)度平均下降率,約為0.6%。在數(shù)值模擬計(jì)算時(shí),對(duì)消落帶范圍巖體黏聚力、內(nèi)摩擦角與滑動(dòng)面黏聚力以0.6%的下降率折減。統(tǒng)計(jì)在不同循環(huán)周期下,巖體強(qiáng)度下降率逐漸增加時(shí),坡體失穩(wěn)所需時(shí)步(坡體內(nèi)部裂縫相互貫通時(shí)坡體失穩(wěn))和變形破壞模式是否發(fā)生變化。

    經(jīng)計(jì)算發(fā)現(xiàn),隨著消落帶范圍巖體強(qiáng)度下降率增加,坡體失穩(wěn)所需時(shí)步逐漸減少,變形破壞模式并不會(huì)發(fā)生變化(見(jiàn)表3)。消落帶巖體強(qiáng)度降低,坡體阻滑段巖體阻抗能力下降,是造成坡體失穩(wěn)速度加快的主要原因,坡體出現(xiàn)彎曲-隆起的時(shí)間也逐步提前。消落帶巖體強(qiáng)度下降程度的高低,僅會(huì)影響岸坡變形演化速度,并不會(huì)改變坡體變形破壞模式,仍以滑移-彎曲為主。此劣化方程未考慮到巖體結(jié)構(gòu)面對(duì)巖體物理力學(xué)強(qiáng)度的影響,會(huì)導(dǎo)致所取下降率依然偏小,坡體失穩(wěn)所需時(shí)步增大。如在后續(xù)的研究中考慮在水位波動(dòng)條件下,岸坡巖體結(jié)構(gòu)面劣化效應(yīng)對(duì)巖體強(qiáng)度下降率的影響,會(huì)更加貼近實(shí)際工況。

    表3 不同循環(huán)周期下坡體失穩(wěn)所需時(shí)步變化Tab.3 Time step of slope instability under different cycle periods

    3.2 巖體劣化對(duì)岸坡變形破壞模式的影響

    巖體強(qiáng)度下降僅是消落帶巖體劣化的表現(xiàn)形式之一,更為常見(jiàn)的是對(duì)岸坡結(jié)構(gòu)的破壞。在巖層面發(fā)育的順層岸坡上,經(jīng)長(zhǎng)年累月的庫(kù)水掏蝕,逐步發(fā)育有大量沿層面與垂直于層面的節(jié)理裂隙切割巖體,使得坡面逐漸碎塊化,并被侵蝕掏空形成凹腔。將消落帶范圍巖體掏蝕深度定為0.97,2.68,4.49,6.74 m,并作為數(shù)值模擬的4種工況(見(jiàn)圖14)。得到不同巖體結(jié)構(gòu)劣化結(jié)果下,對(duì)坡體變形失穩(wěn)所產(chǎn)生的影響。

    當(dāng)掏蝕深度達(dá)0.97 m時(shí),在消落帶區(qū)域形成凹腔,使得坡體阻滑段巖體的阻抗能力急劇下降,失穩(wěn)所需時(shí)步大幅度縮減。由于掏蝕區(qū)域較少,并未影響坡體變形失穩(wěn)模式。當(dāng)掏蝕深度達(dá)2.68 m時(shí),形成的凹腔區(qū)域擴(kuò)大,阻滑段巖體變薄?;w前緣臨空,順著層面直接滑移,破壞模式由滑移-彎曲轉(zhuǎn)變?yōu)榛?見(jiàn)圖15)。當(dāng)淘蝕深度分別為4.49,6.74 m時(shí),凹腔區(qū)域不斷擴(kuò)大,阻滑段巖體逐步消失,滑體前緣臨空程度增大,坡體更易直接滑移,失穩(wěn)所需時(shí)步進(jìn)一步縮減(見(jiàn)表4)。

    表4 不同劣化程度下坡體失穩(wěn)模式變化情況Tab.4 Change of slope instability mode under different deterioration degree

    對(duì)青石6號(hào)坡而言,由于其阻滑段巖體處于三峽庫(kù)區(qū)消落帶范圍內(nèi),長(zhǎng)期庫(kù)水位波動(dòng)的影響會(huì)對(duì)阻滑段巖體強(qiáng)度和結(jié)構(gòu)產(chǎn)生劣化效應(yīng),影響坡體穩(wěn)定性。從上述分析可知,消落帶巖體強(qiáng)度的劣化并不會(huì)改變坡體的變形失穩(wěn)模式;而坡體結(jié)構(gòu)對(duì)劣化程度的響應(yīng)達(dá)到臨界點(diǎn)時(shí),坡體變形破壞模式會(huì)發(fā)生變化。相比岸坡其他區(qū)域,消落帶區(qū)域巖體所處應(yīng)力環(huán)境更復(fù)雜,劣化效應(yīng)極為明顯,影響到未來(lái)庫(kù)岸邊坡變形演化方式。黃波林等也發(fā)現(xiàn):“水-巖-岸”互饋?zhàn)饔茫涌炝藢用尕炌?;同時(shí)由于趨表-后退效應(yīng),不斷下切的消落帶巖體,改變了消落帶巖體結(jié)構(gòu),兩者共同作用加快了岸坡變化速度[25]。王健等也發(fā)現(xiàn)坡腳巖體在不同下切深度下,巖質(zhì)順層岸坡的變形破壞模式會(huì)發(fā)生相互轉(zhuǎn)換[31],印證了本文通過(guò)降低消落帶巖體強(qiáng)度和破壞巖體結(jié)構(gòu),模擬在長(zhǎng)期庫(kù)水波動(dòng)條件下岸坡變形破壞過(guò)程的可行性。

    4 結(jié)論與建議

    以FEM/DEM耦合分析方法作為青石6號(hào)坡數(shù)值模擬手段,結(jié)合室內(nèi)飽和-風(fēng)干試驗(yàn)結(jié)果,對(duì)坡體變形破壞機(jī)理、庫(kù)水對(duì)岸坡?lián)p傷劣化作用與變形破壞模式的影響等方面進(jìn)行研究分析,得出以下結(jié)論:

    (1) 以FEM/DEM耦合分析方法預(yù)測(cè)得到了青石6號(hào)坡變形破壞全過(guò)程,將其劃分為裂縫萌生、擴(kuò)展、彎曲-隆起、裂縫貫通、大規(guī)?;频?個(gè)階段,坡體變形破壞模式為彎曲-隆起。模擬結(jié)果中出現(xiàn)的彎曲-隆起現(xiàn)象與坡體目前出現(xiàn)的變形跡象相對(duì)應(yīng),證明了模擬結(jié)果的準(zhǔn)確性和可靠性。

    (2) 根據(jù)室內(nèi)“飽和-風(fēng)干”試驗(yàn),得到灰?guī)r(T1j3)巖樣飽和單軸抗壓強(qiáng)度下降速率,以此對(duì)消落帶巖體進(jìn)行強(qiáng)度折減。計(jì)算表明,巖體強(qiáng)度的下降僅會(huì)加快岸坡變形演化速度,并不會(huì)影響坡體變形破壞模式。

    (3) 隨著坡體結(jié)構(gòu)劣化程度的增加,致使阻滑段巖體逐漸變薄,阻滑能力降低,滑體前緣臨空程度不斷增大,使得青石6號(hào)坡變形破壞模式由滑移-剪斷向順層滑移轉(zhuǎn)變。

    目前青石6號(hào)坡正處于蠕變階段,但是隨著庫(kù)水位對(duì)消落帶巖體的長(zhǎng)期劣化作用,會(huì)使得岸坡巖體的破碎程度逐漸增加,加快坡體的變形演化速度。整個(gè)青石-抱龍段存在類似的20個(gè)順層巖質(zhì)斜坡,應(yīng)當(dāng)對(duì)其進(jìn)行長(zhǎng)期的監(jiān)測(cè)。針對(duì)各個(gè)斜坡目前所處階段,進(jìn)行相應(yīng)的監(jiān)測(cè)或治理。否則一旦發(fā)生失穩(wěn)破壞,滑體入江會(huì)帶來(lái)嚴(yán)重的涌浪災(zāi)害,對(duì)長(zhǎng)江航道與沿岸居民帶來(lái)嚴(yán)重的威脅。

    猜你喜歡
    順層青石滑體
    珍瓏·青石子度假莊園
    小青石和水滴
    輕質(zhì)高強(qiáng)堇青石多孔陶瓷的制備與表征
    基于三維數(shù)值模擬的含軟弱夾層順層巖質(zhì)邊坡開(kāi)挖穩(wěn)定性研究
    預(yù)應(yīng)力錨索在公路順層巖質(zhì)邊坡中的應(yīng)用
    紅砂巖順層邊坡監(jiān)測(cè)及變形破壞探究
    滑坡碎屑流顆粒分選效應(yīng)的數(shù)值模擬
    立式旋壓機(jī)縱向進(jìn)給機(jī)構(gòu)液壓配重設(shè)計(jì)
    青石崢嶸
    寶藏(2018年1期)2018-01-31 02:05:09
    萬(wàn)梁高速某滑坡降雨入滲穩(wěn)定性及處治技術(shù)研究*
    久久九九热精品免费| 2021天堂中文幕一二区在线观| 亚洲成人久久性| 成人国产一区最新在线观看| 亚洲中文字幕一区二区三区有码在线看 | 男女那种视频在线观看| 亚洲国产中文字幕在线视频| 国产v大片淫在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲成a人片在线一区二区| 国产成人欧美在线观看| 久久午夜综合久久蜜桃| 中文字幕av在线有码专区| 亚洲精品一卡2卡三卡4卡5卡| 色哟哟哟哟哟哟| 国内精品美女久久久久久| av天堂在线播放| 亚洲真实伦在线观看| 国产精品久久久av美女十八| 男人舔奶头视频| 我的老师免费观看完整版| 欧美绝顶高潮抽搐喷水| 国产成+人综合+亚洲专区| 波多野结衣巨乳人妻| 亚洲国产日韩欧美精品在线观看 | 国产高清三级在线| 国产激情偷乱视频一区二区| 99国产精品一区二区三区| 99久久无色码亚洲精品果冻| 午夜福利免费观看在线| 日本 欧美在线| 午夜福利视频1000在线观看| xxxwww97欧美| 天堂影院成人在线观看| 韩国av一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久九九精品二区国产| 2021天堂中文幕一二区在线观| 亚洲18禁久久av| 免费高清视频大片| 我要搜黄色片| 国产成人精品久久二区二区91| 国产精华一区二区三区| 很黄的视频免费| 成熟少妇高潮喷水视频| 人人妻人人看人人澡| 亚洲av五月六月丁香网| 亚洲精品一卡2卡三卡4卡5卡| 人妻夜夜爽99麻豆av| 亚洲精品一区av在线观看| 久久精品国产亚洲av香蕉五月| 搡老岳熟女国产| 午夜福利在线在线| 日韩免费av在线播放| 黑人巨大精品欧美一区二区mp4| 午夜福利高清视频| 最新中文字幕久久久久 | 九九久久精品国产亚洲av麻豆 | 中亚洲国语对白在线视频| 俺也久久电影网| e午夜精品久久久久久久| 国内久久婷婷六月综合欲色啪| 最新美女视频免费是黄的| 国产伦精品一区二区三区四那| 国产伦在线观看视频一区| h日本视频在线播放| 亚洲成人免费电影在线观看| ponron亚洲| 美女午夜性视频免费| 国产精品一区二区三区四区久久| 午夜精品在线福利| 亚洲在线观看片| 人妻夜夜爽99麻豆av| 床上黄色一级片| 久久久国产成人免费| 欧美最黄视频在线播放免费| 国产伦在线观看视频一区| 波多野结衣高清无吗| 亚洲av电影在线进入| aaaaa片日本免费| 亚洲成av人片在线播放无| 好男人电影高清在线观看| 欧美黑人欧美精品刺激| 亚洲欧美激情综合另类| 色老头精品视频在线观看| 男人舔女人下体高潮全视频| 两性夫妻黄色片| 级片在线观看| 好男人电影高清在线观看| 天天添夜夜摸| 给我免费播放毛片高清在线观看| 午夜精品久久久久久毛片777| 在线看三级毛片| 国产97色在线日韩免费| 黄色女人牲交| 久久婷婷人人爽人人干人人爱| 亚洲激情在线av| 成年女人永久免费观看视频| 俄罗斯特黄特色一大片| 国产黄色小视频在线观看| 成人无遮挡网站| 99热6这里只有精品| 亚洲国产中文字幕在线视频| 亚洲国产日韩欧美精品在线观看 | 久久久国产成人精品二区| 日本黄色片子视频| 国产麻豆成人av免费视频| 亚洲人与动物交配视频| 日本撒尿小便嘘嘘汇集6| 精品久久蜜臀av无| 久久久精品大字幕| 国产伦在线观看视频一区| 小说图片视频综合网站| 婷婷精品国产亚洲av在线| 丰满人妻一区二区三区视频av | 成人av在线播放网站| 国产激情欧美一区二区| 2021天堂中文幕一二区在线观| 亚洲七黄色美女视频| 精品久久久久久成人av| 欧美日韩黄片免| 女生性感内裤真人,穿戴方法视频| 中文字幕精品亚洲无线码一区| 亚洲一区二区三区不卡视频| 免费电影在线观看免费观看| 麻豆成人午夜福利视频| 99国产精品99久久久久| 精品国产三级普通话版| 亚洲五月婷婷丁香| 国产探花在线观看一区二区| 欧美av亚洲av综合av国产av| 亚洲精品久久国产高清桃花| 麻豆一二三区av精品| 在线a可以看的网站| 日韩欧美国产一区二区入口| 免费电影在线观看免费观看| 天天添夜夜摸| 国产精品爽爽va在线观看网站| 欧美成人性av电影在线观看| 两个人视频免费观看高清| 丝袜人妻中文字幕| 午夜激情福利司机影院| 亚洲在线自拍视频| 亚洲欧洲精品一区二区精品久久久| 99久久综合精品五月天人人| 最近最新中文字幕大全电影3| 国产精品电影一区二区三区| 亚洲成人免费电影在线观看| 一a级毛片在线观看| 波多野结衣高清作品| 男女视频在线观看网站免费| 狂野欧美白嫩少妇大欣赏| 91av网一区二区| 好男人电影高清在线观看| 丰满人妻一区二区三区视频av | 亚洲18禁久久av| 女生性感内裤真人,穿戴方法视频| 亚洲激情在线av| 99riav亚洲国产免费| 成人18禁在线播放| av天堂中文字幕网| 久久精品国产亚洲av香蕉五月| 日本免费a在线| 国产成人一区二区三区免费视频网站| 欧美最黄视频在线播放免费| 老司机在亚洲福利影院| 亚洲欧美日韩无卡精品| 天堂√8在线中文| 欧美xxxx黑人xx丫x性爽| 首页视频小说图片口味搜索| 少妇的逼水好多| 国产精品久久久久久亚洲av鲁大| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 一二三四社区在线视频社区8| aaaaa片日本免费| 在线观看66精品国产| 床上黄色一级片| 久久精品国产99精品国产亚洲性色| 一个人看的www免费观看视频| 亚洲国产欧洲综合997久久,| 亚洲精品乱码久久久v下载方式 | 久久精品国产清高在天天线| 欧美3d第一页| 国产一区二区三区视频了| 国产精品av久久久久免费| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人澡欧美一区二区| 久久久色成人| 一二三四在线观看免费中文在| 麻豆国产97在线/欧美| 国产成人欧美在线观看| 九九在线视频观看精品| 性色av乱码一区二区三区2| 久久热在线av| 亚洲国产看品久久| 久久久色成人| 一级黄色大片毛片| 日本a在线网址| 国产精品一区二区免费欧美| 无遮挡黄片免费观看| www.www免费av| 久久久久久久久中文| 精品国产三级普通话版| 99热6这里只有精品| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 亚洲av成人精品一区久久| 亚洲欧美日韩无卡精品| 亚洲人与动物交配视频| 国产精品亚洲一级av第二区| 久久久精品欧美日韩精品| 日韩欧美国产在线观看| 亚洲成人中文字幕在线播放| 中文字幕精品亚洲无线码一区| av天堂中文字幕网| 国内毛片毛片毛片毛片毛片| 免费观看精品视频网站| 国产精品一及| 精品99又大又爽又粗少妇毛片 | 成人鲁丝片一二三区免费| 一区二区三区激情视频| 母亲3免费完整高清在线观看| 国内精品久久久久精免费| 美女午夜性视频免费| 在线观看日韩欧美| 国产亚洲精品久久久久久毛片| 天天添夜夜摸| 亚洲av成人精品一区久久| 啦啦啦免费观看视频1| 欧美日韩一级在线毛片| 熟女电影av网| 国产成人精品久久二区二区91| 成人欧美大片| 亚洲精品色激情综合| 欧美中文日本在线观看视频| 又粗又爽又猛毛片免费看| 中亚洲国语对白在线视频| 国产亚洲精品av在线| 久久精品国产清高在天天线| 欧美日韩亚洲国产一区二区在线观看| 欧美日本亚洲视频在线播放| 久久久色成人| 成人特级黄色片久久久久久久| 久久久精品欧美日韩精品| 国产高清三级在线| 久久婷婷人人爽人人干人人爱| 老司机午夜十八禁免费视频| 18禁观看日本| 欧美成人一区二区免费高清观看 | 99久久无色码亚洲精品果冻| 一个人看视频在线观看www免费 | 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 舔av片在线| 夜夜躁狠狠躁天天躁| 最新中文字幕久久久久 | 我的老师免费观看完整版| 母亲3免费完整高清在线观看| 欧美高清成人免费视频www| 99热只有精品国产| 日韩欧美国产在线观看| 婷婷六月久久综合丁香| 精品午夜福利视频在线观看一区| 久久国产乱子伦精品免费另类| 男插女下体视频免费在线播放| 757午夜福利合集在线观看| 国产精品九九99| 国产人伦9x9x在线观看| 精品国产三级普通话版| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| 午夜精品在线福利| h日本视频在线播放| 香蕉丝袜av| 国产伦精品一区二区三区四那| 久久婷婷人人爽人人干人人爱| 欧美3d第一页| 午夜福利在线在线| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 国产成人精品久久二区二区91| netflix在线观看网站| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 国产三级在线视频| 老司机午夜十八禁免费视频| 日韩欧美国产一区二区入口| 国产成人av教育| 日本黄大片高清| 麻豆av在线久日| 亚洲va日本ⅴa欧美va伊人久久| 欧美一区二区国产精品久久精品| 亚洲 欧美一区二区三区| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 搡老岳熟女国产| 国内精品久久久久久久电影| 神马国产精品三级电影在线观看| 国产蜜桃级精品一区二区三区| 在线免费观看不下载黄p国产 | 久久午夜亚洲精品久久| 国产精品久久久久久久电影 | 成人三级做爰电影| 欧美av亚洲av综合av国产av| 久久久久性生活片| 免费在线观看影片大全网站| tocl精华| 国产午夜福利久久久久久| 午夜福利在线观看吧| 特级一级黄色大片| 欧美性猛交╳xxx乱大交人| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 欧美大码av| 一级a爱片免费观看的视频| 一级作爱视频免费观看| 精品一区二区三区视频在线观看免费| 一本久久中文字幕| 免费无遮挡裸体视频| 日本a在线网址| 好看av亚洲va欧美ⅴa在| 久久久国产精品麻豆| 国内精品久久久久精免费| 亚洲,欧美精品.| av黄色大香蕉| 18禁国产床啪视频网站| 国产高清三级在线| 亚洲18禁久久av| 久久国产精品人妻蜜桃| 99精品欧美一区二区三区四区| 精品久久久久久久久久久久久| 国产又色又爽无遮挡免费看| 特大巨黑吊av在线直播| 五月伊人婷婷丁香| 99久国产av精品| 亚洲自拍偷在线| 亚洲欧美激情综合另类| 欧美日本亚洲视频在线播放| 国产成人欧美在线观看| 欧美在线黄色| 欧美成人免费av一区二区三区| 男女午夜视频在线观看| 国产成人欧美在线观看| 99re在线观看精品视频| 欧美中文综合在线视频| 一进一出好大好爽视频| 国产高清三级在线| 一二三四社区在线视频社区8| 日本一本二区三区精品| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产欧美网| 天堂√8在线中文| 日韩大尺度精品在线看网址| 国产一区二区激情短视频| 久久久精品大字幕| 19禁男女啪啪无遮挡网站| 午夜精品一区二区三区免费看| 岛国视频午夜一区免费看| 久久午夜综合久久蜜桃| 亚洲熟妇中文字幕五十中出| 一级毛片精品| 一个人免费在线观看的高清视频| 国产精品亚洲美女久久久| 欧美性猛交╳xxx乱大交人| 亚洲成av人片在线播放无| 中文字幕熟女人妻在线| 一个人看视频在线观看www免费 | 国产av一区在线观看免费| 国产主播在线观看一区二区| 久久午夜综合久久蜜桃| 长腿黑丝高跟| 99国产精品99久久久久| 99国产综合亚洲精品| 国产亚洲精品一区二区www| 国产一区二区三区在线臀色熟女| 欧美日韩国产亚洲二区| 日本a在线网址| 久久久久国产精品人妻aⅴ院| 51午夜福利影视在线观看| 免费看光身美女| 91在线精品国自产拍蜜月 | 亚洲 欧美 日韩 在线 免费| 精品久久久久久久毛片微露脸| 最新美女视频免费是黄的| 18禁黄网站禁片免费观看直播| 熟女人妻精品中文字幕| 伊人久久大香线蕉亚洲五| 午夜激情欧美在线| 亚洲美女黄片视频| 婷婷丁香在线五月| 在线永久观看黄色视频| 天堂网av新在线| 变态另类丝袜制服| 亚洲国产日韩欧美精品在线观看 | 成年女人毛片免费观看观看9| 在线视频色国产色| 亚洲国产精品久久男人天堂| 欧美黄色片欧美黄色片| 免费一级毛片在线播放高清视频| 成在线人永久免费视频| 99国产精品一区二区三区| 欧美在线一区亚洲| 国产精品亚洲一级av第二区| 国产淫片久久久久久久久 | 亚洲av成人一区二区三| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 欧美zozozo另类| 国语自产精品视频在线第100页| 老鸭窝网址在线观看| 精品不卡国产一区二区三区| 99riav亚洲国产免费| 亚洲中文字幕一区二区三区有码在线看 | 欧美高清成人免费视频www| 久久中文字幕人妻熟女| 色老头精品视频在线观看| 欧美xxxx黑人xx丫x性爽| 日本黄色视频三级网站网址| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 狂野欧美激情性xxxx| 一区福利在线观看| 亚洲在线自拍视频| 精品久久久久久久毛片微露脸| 国产免费男女视频| 老司机午夜十八禁免费视频| 午夜福利在线在线| 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 国内毛片毛片毛片毛片毛片| 国产精品一区二区精品视频观看| 国产亚洲av嫩草精品影院| 一区二区三区高清视频在线| 国产日本99.免费观看| 亚洲自拍偷在线| 熟女少妇亚洲综合色aaa.| 99久久精品热视频| 精品国内亚洲2022精品成人| 精品国产乱码久久久久久男人| 亚洲午夜精品一区,二区,三区| 国产精品 国内视频| 亚洲精品在线美女| 黄色片一级片一级黄色片| 亚洲精品在线观看二区| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 99热精品在线国产| 国产精品女同一区二区软件 | 黄片大片在线免费观看| 久久精品人妻少妇| 婷婷丁香在线五月| 亚洲av美国av| 日韩欧美在线乱码| 99国产精品一区二区蜜桃av| 亚洲中文字幕一区二区三区有码在线看 | 婷婷亚洲欧美| or卡值多少钱| 免费无遮挡裸体视频| 午夜福利视频1000在线观看| 国产精品影院久久| 午夜成年电影在线免费观看| 欧美又色又爽又黄视频| 欧美另类亚洲清纯唯美| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 午夜福利免费观看在线| 成年版毛片免费区| 欧美色视频一区免费| 久久久久亚洲av毛片大全| 校园春色视频在线观看| 日本精品一区二区三区蜜桃| 国产亚洲精品av在线| 男女视频在线观看网站免费| 精品不卡国产一区二区三区| 俺也久久电影网| 成人国产一区最新在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲欧美日韩东京热| 日本免费一区二区三区高清不卡| 在线观看免费午夜福利视频| 97碰自拍视频| av片东京热男人的天堂| 亚洲国产精品sss在线观看| 狂野欧美激情性xxxx| 18禁黄网站禁片免费观看直播| 99久久综合精品五月天人人| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看| 男人舔女人下体高潮全视频| 国产激情久久老熟女| 男人舔女人下体高潮全视频| av片东京热男人的天堂| av中文乱码字幕在线| av片东京热男人的天堂| 99re在线观看精品视频| 午夜成年电影在线免费观看| 男女下面进入的视频免费午夜| 一个人免费在线观看电影 | 国产欧美日韩精品亚洲av| 成人欧美大片| 久久久久久大精品| 国产激情偷乱视频一区二区| 国产精品久久久久久精品电影| 国产乱人视频| 国产成人一区二区三区免费视频网站| 免费看a级黄色片| 欧美日韩亚洲国产一区二区在线观看| 国内揄拍国产精品人妻在线| 久久久久精品国产欧美久久久| 日本 av在线| 亚洲天堂国产精品一区在线| 哪里可以看免费的av片| 成人欧美大片| 亚洲性夜色夜夜综合| 亚洲va日本ⅴa欧美va伊人久久| 一个人观看的视频www高清免费观看 | 成年版毛片免费区| 99riav亚洲国产免费| 欧美三级亚洲精品| 亚洲欧美日韩高清专用| 国产一区二区三区视频了| 久久这里只有精品19| 在线视频色国产色| 欧美中文综合在线视频| 国产高清激情床上av| 制服人妻中文乱码| 色综合站精品国产| 国产亚洲精品综合一区在线观看| 性色av乱码一区二区三区2| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 久久精品人妻少妇| 婷婷丁香在线五月| 国产91精品成人一区二区三区| 丁香六月欧美| 一个人免费在线观看电影 | 日韩av在线大香蕉| 日本黄色视频三级网站网址| 国产美女午夜福利| 亚洲成人中文字幕在线播放| 一本精品99久久精品77| 国产成人福利小说| 久久午夜综合久久蜜桃| 久久久久久久午夜电影| 日韩国内少妇激情av| 又大又爽又粗| 成人av一区二区三区在线看| 国产97色在线日韩免费| 久久这里只有精品中国| 脱女人内裤的视频| 久久久久久大精品| 欧美日韩瑟瑟在线播放| 黄片大片在线免费观看| 91字幕亚洲| 少妇的丰满在线观看| 他把我摸到了高潮在线观看| 国产高清三级在线| 精品久久久久久久毛片微露脸| 激情在线观看视频在线高清| 18美女黄网站色大片免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产成人aa在线观看| 午夜免费成人在线视频| 国产v大片淫在线免费观看| 久久天堂一区二区三区四区| 国产激情欧美一区二区| 免费观看人在逋| 在线观看66精品国产| 桃红色精品国产亚洲av| 欧美又色又爽又黄视频| 法律面前人人平等表现在哪些方面| 长腿黑丝高跟| 国产免费男女视频| 亚洲国产欧美一区二区综合| 欧美日韩黄片免| 欧美丝袜亚洲另类 | 久久精品影院6| 亚洲av片天天在线观看| 亚洲黑人精品在线| 亚洲电影在线观看av| 国产日本99.免费观看| 最新中文字幕久久久久 | 亚洲成人精品中文字幕电影| 午夜成年电影在线免费观看| 十八禁网站免费在线| 国产亚洲精品av在线| 国产精品 国内视频| 大型黄色视频在线免费观看| 色综合欧美亚洲国产小说| 亚洲激情在线av| 美女被艹到高潮喷水动态| 嫩草影院精品99| 成年版毛片免费区| 国产真实乱freesex| 精品久久久久久,| 色尼玛亚洲综合影院| 亚洲精品色激情综合| 韩国av一区二区三区四区| 亚洲国产欧美一区二区综合| 国产午夜福利久久久久久| 国产日本99.免费观看| 成人三级黄色视频| 中文在线观看免费www的网站| 色综合站精品国产| h日本视频在线播放| 一进一出抽搐动态| 精华霜和精华液先用哪个| 亚洲在线观看片| 51午夜福利影视在线观看| 不卡一级毛片| 欧美中文日本在线观看视频|