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

    金沙江白格滑坡殘留體穩(wěn)定性及堵江風(fēng)險(xiǎn)分析

    2023-05-21 00:16:45唐岳灝姜清輝
    水利水電快報(bào) 2023年5期

    唐岳灝 姜清輝

    摘要:

    為評估2018年金沙江白格村兩次特大型滑坡堰塞堵江事件后所形成潛在不穩(wěn)定塊體(殘留體)的穩(wěn)定性及堵江風(fēng)險(xiǎn),在應(yīng)用嚴(yán)格三維極限平衡分析白格滑坡殘留體穩(wěn)定性的基礎(chǔ)上,結(jié)合應(yīng)急處置方案,基于SPH-DEM流固耦合模型,對殘留體滑坡范圍和堆積體厚度進(jìn)行了預(yù)測,定量評估崩滑堵江風(fēng)險(xiǎn)。結(jié)果表明:① K1-Ⅰ和K2-Ⅰ滑塊處于臨界狀態(tài),失穩(wěn)風(fēng)險(xiǎn)較大,存在再次發(fā)生滑坡堵江的可能;② 河道清淤方法可降低堰塞體整體高度,增大滑坡物質(zhì)堆積空間,有效降低再次滑坡堵江的風(fēng)險(xiǎn)。研究成果可為高位滑坡堵江風(fēng)險(xiǎn)預(yù)測提供一定參考。

    關(guān)鍵詞:

    白格殘留體; 滑坡堵江; SPH-DEM; 數(shù)值模擬; 風(fēng)險(xiǎn)評估

    中圖法分類號:P642.22

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

    DOI:10.15974/j.cnki.slsdkb.2023.05.006

    文章編號:1006-0081(2023)05-0038-07

    0 引 言

    2018年10月10日和11月3日,在西藏自治區(qū)江達(dá)縣波羅鄉(xiāng)白格村與四川省白玉縣絨蓋鄉(xiāng)則巴村交界處的金沙江西藏岸,先后兩次發(fā)生大規(guī)模高位滑坡[1-4],堵塞金沙江,形成堰塞湖?;卵呷蜎]了上游的村莊和各種生產(chǎn)生活設(shè)施,同時(shí)潰壩洪水嚴(yán)重沖毀滑坡下游的村莊、農(nóng)田和公路、橋梁等基礎(chǔ)設(shè)施[5-6]。白格滑坡經(jīng)兩次滑動后,松散物質(zhì)已大為減少,但滑坡體對山體擾動作用巨大(主要表現(xiàn)為拉、拽、刮、鏟、刷),滑坡三面由于臨空卸荷作用不斷加強(qiáng),變形跡象十分明顯,后緣及兩側(cè)發(fā)育多條深大裂縫,存在大量滑坡殘留體。無人機(jī)航空攝影測量及工程地質(zhì)測繪、鉆探、物探等詳細(xì)勘察工作的成果表明:前兩次白格滑坡周界后緣有3個(gè)主要?dú)埩糇冃误w,其滑坡周界明顯,已經(jīng)發(fā)生了明顯的下滑、錯(cuò)動,仍具備進(jìn)一步下滑的趨勢。同時(shí),地表宏觀變形跡象和監(jiān)測成果表明殘留體一直存在前緣溜滑和蠕滑變形,其穩(wěn)定安全問題不容樂觀。因此,開展白格滑坡殘留體的防災(zāi)減災(zāi)研究,分析殘留體的穩(wěn)定性和變形趨勢,可為滑坡應(yīng)急處置提供支撐,并為白格滑坡第三次堵江的可能性、危害性預(yù)測提供基本數(shù)據(jù)。

    目前,已有一些學(xué)者利用數(shù)值方法對白格殘留體堵江風(fēng)險(xiǎn)進(jìn)行風(fēng)險(xiǎn)評價(jià)。曹水合等[7]利用數(shù)值模擬軟件Massflow對白格殘留體堵江范圍和高度進(jìn)行了風(fēng)險(xiǎn)預(yù)測研究,趙程等[8]利用MassMov2D碎屑流模擬軟件對白格滑坡滑源區(qū)的3處潛在不穩(wěn)定巖體進(jìn)行了預(yù)測,周禮等[9]利用PFC3D軟件對白格滑坡滑源區(qū)殘留不穩(wěn)定部分可能失穩(wěn)的運(yùn)動路徑和堆積范圍進(jìn)行預(yù)測。蔡耀軍等[10]利用PFC3D軟件對殘留體不同失穩(wěn)規(guī)模進(jìn)行了堰塞體堆積形態(tài)預(yù)測。盡管這些數(shù)值仿真取得了一定進(jìn)展,但也存在一些不足:這些軟件大多使用固體力學(xué)的方法模擬滑坡動力學(xué)過程,無法處理滑坡堵江產(chǎn)生的流固耦合問題;有的將滑坡過程簡化為流體流動過程,導(dǎo)致模擬結(jié)果與實(shí)際情況不完全對應(yīng)。此外,上述大部分軟件都是基于國外商用程序,后處理的二次開發(fā)功能比較有限,在精細(xì)化展示堆積體的三維幾何形態(tài)方面缺乏可擴(kuò)展性。

    本文首先采用嚴(yán)格三維極限平衡分析方法分析計(jì)算滑體安全系數(shù),同時(shí)確定滑體主滑方向,建立三維滑坡模型。然后使用基于拉格朗日力學(xué)體系的離散元(DEM)方法和光滑粒子流體動力學(xué)(SPH)方法完全自主開發(fā)的SPH-DEM耦合程序——利用DEM模擬大變形的滑坡動力過程、利用SPH模擬碎屑體入江產(chǎn)生的水動力學(xué)問題,以更好地反映滑坡(固體)與水體(流體)間的相互作用。通過SPH-DEM耦合程序可以實(shí)現(xiàn)真實(shí)場景下的滑坡堵江動力學(xué)過程模擬,對強(qiáng)變形區(qū)里不穩(wěn)定體崩落后的發(fā)展過程和堆積形態(tài)進(jìn)行預(yù)測,定量評估崩落后再次堵江的風(fēng)險(xiǎn),為潰壩及洪水災(zāi)害分析做參考。

    1 滑坡體特征

    白格滑坡地質(zhì)剖面如圖1所示。根據(jù)勘察,按殘留體變形特征與空間位置關(guān)系,將殘留體劃分為3個(gè)區(qū)塊(圖2):滑坡后緣為K1,滑坡左側(cè)(金沙江下游側(cè))為K2,滑坡右側(cè)(金沙江上游)為K3。K1,K2和K3都是由滑坡牽引卸荷形成。各殘留體面積較大,根據(jù)不同變形邊界特征又可分為多個(gè)子塊,其中,K1殘留體細(xì)分為K1-Ⅰ,K1-Ⅱ和K1-Ⅲ,K2殘留體細(xì)分為K2-Ⅰ,K2-Ⅱ。

    根據(jù)《西藏自治區(qū)昌都市江達(dá)縣波羅鄉(xiāng)金沙江白格滑坡應(yīng)急勘查報(bào)告》:K1-Ⅰ和K1-Ⅲ殘留體潛在滑帶為元古界雄松群千枚巖組,巖性以絹云母石英千枚巖為主,為深灰色、灰黑色,變余泥質(zhì)結(jié)構(gòu),千枚狀構(gòu)造,礦物成分以絹云母為主,含少量石墨、石英;K1-Ⅱ潛在滑帶為片麻巖與侵入蛇紋巖接觸面,巖性以全/強(qiáng)風(fēng)化蛇綠巖為主,深灰綠色,風(fēng)化呈淺綠-淺白色,碎塊石土狀;K2-Ⅰ殘留體潛在滑帶為元古界雄松群千枚巖組,巖性以絹云母石英千枚巖為主,深灰色、灰黑色,變余泥質(zhì)結(jié)構(gòu),千枚狀構(gòu)造,礦物成分以絹云母為主,含少量石墨、石英;K2-Ⅱ殘留體潛在滑面不明顯;K3潛在滑帶為元古界雄松群片麻巖強(qiáng)風(fēng)化層。典型巖石樣品如圖3所示。

    2 滑坡后緣殘留體穩(wěn)定性評價(jià)

    2.1 計(jì)算模型和計(jì)算方法

    根據(jù)白格滑坡殘留體的地形等高線、滑坡邊界以及地質(zhì)鉆孔、深部測斜儀監(jiān)測確定的滑面深度,建立白格滑坡殘留體(K1-Ⅰ,K1-Ⅱ,K1,K2-Ⅰ,K2-Ⅱ,K3)的三維穩(wěn)定分析計(jì)算模型。計(jì)算方法采用滑源區(qū)嚴(yán)格三維極限平衡分析法,在對現(xiàn)有的邊坡穩(wěn)定性定量評價(jià)方法進(jìn)行總結(jié)分析的基礎(chǔ)上,克服了現(xiàn)有滑坡穩(wěn)定性三維分析需要預(yù)先給定主滑方向以及滑體只滿足部分平衡條件的局限性,適用于多結(jié)構(gòu)面組合形成、具有復(fù)雜空間形態(tài)滑裂面的滑坡嚴(yán)格三維穩(wěn)定分析理論,可以同時(shí)計(jì)算得到滑體的安全系數(shù)和主滑方向,且嚴(yán)格滿足空間力系的6個(gè)平衡方程,能夠?qū)崿F(xiàn)對邊坡/滑坡在自然、降雨、地震不同工況條件下的穩(wěn)定狀態(tài)或加固后斜坡穩(wěn)定性的準(zhǔn)確評估。圖4為計(jì)算所得的K1-Ⅰ、K1-Ⅱ和K2-Ⅰ滑面。

    2.2 計(jì)算工況和參數(shù)取值

    根據(jù)滑坡殘留體的實(shí)際情況,可能作用在殘留體上的荷載有:坡體自重、水壓力和地震力。① 坡體自重。在天然狀態(tài)下,坡體自重按天然重度計(jì)算。據(jù)當(dāng)?shù)厮臍庀筚Y料及斜坡巖性組成的實(shí)際情況,在連續(xù)降雨或暴雨條件下,降水入滲深度按全入滲計(jì),入滲范圍內(nèi)的土體重度按飽和重度計(jì)。② 水壓力。地下水產(chǎn)生的荷載主要是降雨誘發(fā)坡面暫態(tài)飽和區(qū)產(chǎn)生的地下水壓力等。③ 地震力。工區(qū)地震設(shè)防烈度為7度,設(shè)計(jì)地震分組為第二組,設(shè)計(jì)基本地震加速度值為0.20 g 。計(jì)算工況考慮3種:正常工況(自重荷載)、降雨工況(自重荷載+降雨水壓力)、地震工況(自重荷載+地震荷載)。殘留體滑帶力學(xué)計(jì)算參數(shù)主要依據(jù)巖土力學(xué)試驗(yàn)和工程經(jīng)驗(yàn)選取,如表1所示。

    2.3 穩(wěn)定性評估與穩(wěn)定狀態(tài)分級

    采用三維嚴(yán)格極限平衡法計(jì)算得到的白格滑坡后緣殘留體K1-Ⅰ,K1-Ⅱ,K1,K2-Ⅰ,K2-Ⅱ和K3在不同工況條件下穩(wěn)定性安全系數(shù)以及各個(gè)滑塊的體積、主滑方向和穩(wěn)定狀態(tài),如表2所示,其中主滑方向角為與正北方向夾角,順時(shí)針為正。穩(wěn)定狀態(tài)根據(jù)《青藏高原重大滑坡風(fēng)險(xiǎn)防控指南》中滑坡穩(wěn)定狀態(tài)分為不穩(wěn)定 (K<1.0);臨界穩(wěn)定(1.0≤K≤1.05);弱穩(wěn)定(1.05

    根據(jù)殘留體在正常、降雨、地震工況下的穩(wěn)定性分區(qū)、穩(wěn)定狀態(tài)分級和相應(yīng)的預(yù)警等級,得到殘留體的風(fēng)險(xiǎn)等級(圖2),為下一步進(jìn)行滑坡殘留體應(yīng)急處置和堵江風(fēng)險(xiǎn)評估提供了直觀的參考。圖5為計(jì)算得到的各滑塊主滑方向與地表外觀測點(diǎn)位移矢量的對比,從中可以看出,主滑方向與外觀測點(diǎn)變形方向吻合,進(jìn)一步論證了所使用的三維嚴(yán)格分析方法的有效性和可靠性。

    基于對滑源區(qū)殘留體的穩(wěn)定性評估,可以看出K1-Ⅰ和K2-Ⅰ滑塊處于臨界狀態(tài),失穩(wěn)風(fēng)險(xiǎn)較大,因此,首先考慮的失穩(wěn)模式為K1-Ⅰ+K2-Ⅰ發(fā)生崩滑(模式1);其次考慮處于臨界狀態(tài)和弱穩(wěn)定狀態(tài)的滑塊發(fā)生失穩(wěn),即K1-Ⅰ+K2-Ⅰ+K1-Ⅱ+K3發(fā)生崩滑(模式2)。

    3 滑坡體應(yīng)急處置

    由于白格滑坡后緣殘留體K1-Ⅰ和K2-Ⅰ處于臨界穩(wěn)定狀態(tài),存在再次下滑和堵江風(fēng)險(xiǎn),因此,需要采取應(yīng)急處置措施。對白格滑坡后緣殘留變形體的處置原則采用滑源區(qū)—運(yùn)移區(qū)—堆積區(qū)全過程處理中的“顧兩頭”,分別為滑源區(qū)的處置和金沙江河道堰塞體的處置?;磪^(qū)應(yīng)急處置采用的方案是削方減載,主要對K1-Ⅰ殘留體進(jìn)行削坡,其目的是提高殘留體的穩(wěn)定性,降低滑坡風(fēng)險(xiǎn),如圖6所示。削坡坡比根據(jù)GB 50330-2013《建筑邊坡工程技術(shù)規(guī)范》坡率法,按碎石土邊坡坡率允許值,結(jié)合深大裂縫分布位置,確定為1∶1.7~1∶1.3,分1~8級,坡高8 m,馬道寬5~6 m;總開挖體積約48萬m3。堰塞堆積區(qū)采用的方案是河道清淤,主要針對四川省一側(cè)開展,主要目的是降低堰塞體整體高度,提高水流下泄能力,增大滑坡物質(zhì)堆積空間,降低再次滑坡堵江的風(fēng)險(xiǎn)。具體方案為將堰塞體殘?bào)w中靠近河道側(cè)不小于40.0 m范圍內(nèi)、2 920 m高程以上的殘?bào)w全部清除。開挖設(shè)計(jì)方案共有4級馬道,馬道寬度0.2 m,放坡坡比1∶2.0,每一級臺階高度20 m,總開挖體積約245萬m3,如圖7所示。

    4 滑坡殘留體堵江評估

    4.1 光滑粒子流體動力學(xué)(SPH)

    4.4 模型計(jì)算參數(shù)選擇

    從DEM建模角度,需要將滑坡碎屑體概化為等直徑的小球。由于潛在滑坡體的體積超過500萬m3(模式2),如果采用小顆粒構(gòu)建滑坡體會導(dǎo)致顆粒數(shù)量巨大,嚴(yán)重影響DEM和SPH模型的計(jì)算效率,但過大的顆粒會影響模型計(jì)算精度。經(jīng)反復(fù)驗(yàn)算,最終確定采用直徑2.4 m的球形顆粒。模型計(jì)算參數(shù)見表3。

    4.5 堵江風(fēng)險(xiǎn)評估

    采用滑坡啟動—運(yùn)移—堵江災(zāi)害鏈全過程模擬的SPH-DEM耦合方法,對白格滑坡后緣殘留體不同塊體組合形成的失穩(wěn)模式在失穩(wěn)觸發(fā)后的動力學(xué)過程及其入江后與水流相互作用的過程進(jìn)行了模擬。根據(jù)滑坡最后形成的堰塞堆積體形態(tài)對滑坡堵江的危害性進(jìn)行了評估。計(jì)算工況:① K1-Ⅰ+K2-Ⅰ(模式1);② K1-Ⅰ+K2-Ⅰ+K1-Ⅱ+K3(模式2)。圖8~9為兩種失穩(wěn)模式最后形成的堰塞體堆積體厚度和橫剖面、縱剖面。為了對比分析,進(jìn)行源區(qū)和河道工程應(yīng)急處置后,相應(yīng)失穩(wěn)模式造成的堰塞體堆積體厚度和橫剖面如圖10~11所示。表4為應(yīng)急處置前后滑坡堰塞壩幾何特征統(tǒng)計(jì)。

    根據(jù)圖8~11和表4分析所得結(jié)果如下。

    (1) 應(yīng)急處置前,兩種失穩(wěn)模式堆積體方量分別為440萬m3和527萬m3,最后堆積形成的堰塞壩高分別為60.2 m和63.0 m。相對模式1,模式2的滑塊方量增加了87萬m3,但最后形成堰塞壩的高度僅增加了約3 m。主要原因是模式2中的K1-Ⅱ滑塊和K3滑塊主要位于上游側(cè)坡,根據(jù)其失穩(wěn)觸發(fā)后運(yùn)移軌跡,大部分滑坡物質(zhì)堆積到靠金沙江上游側(cè)河道,因此,模式2形成堰塞壩的寬度(順河向)與模式1相比增加了446 m,但壩高增加較少。

    (2) 采用應(yīng)急處置措施后,金沙江河道清淤245萬m3,為后緣殘留體失穩(wěn)后形成的滑坡壩提供了堆積空間,因此,模式1和模式2形成的堰塞壩高度分別為47.0 m和50.5 m,相對應(yīng)急處置前壩高減少了約13 m,減少了堰塞湖庫容,有效降低了滑坡堰塞湖的危害性。

    (3) 根據(jù)計(jì)算結(jié)果,按照最危險(xiǎn)的情況(模式2),白格后緣殘留體處于臨界穩(wěn)定和弱穩(wěn)定的子滑塊觸發(fā)失穩(wěn)后,其堆積形成的最大壩高為50.5 m,與白格“10·10”滑坡堰塞壩高程大致相同。考慮滑坡區(qū)交通條件大為改善、人工干預(yù)簡單易行等因素,裂縫區(qū)滑坡堵江及其災(zāi)害鏈風(fēng)險(xiǎn)處于可控范圍。

    5 結(jié) 論

    (1) 本文采用嚴(yán)格三維穩(wěn)定性分析方法對白格后緣殘留體6個(gè)滑塊(K1-Ⅰ、K1-Ⅱ、K1、K2-Ⅰ、K2-Ⅱ、K3)在不同工況條件下的穩(wěn)定性進(jìn)行評估,得到滑塊的安全系數(shù)和主滑方向,計(jì)算結(jié)果與現(xiàn)場宏觀定性分析和外觀監(jiān)測得到的各子區(qū)域變形方向吻合。根據(jù)滑塊穩(wěn)定性分析結(jié)果,對白格后緣殘留體進(jìn)行了穩(wěn)定性分區(qū)、分級,并給出了分級監(jiān)測預(yù)警建議。

    (2) 針對堆積區(qū)的河道清淤應(yīng)急處置,采用滑坡災(zāi)害鏈全過程模擬SPH-DEM方法,對后緣殘留體不同失穩(wěn)模式可能導(dǎo)致的滑坡堵江效果進(jìn)行了分析。數(shù)值模擬結(jié)果表明:金沙江河道清淤為后緣殘留體失穩(wěn)后形成的滑坡壩提供了堆積空間,減少了堰塞湖庫容,有效降低了滑坡堰塞湖的危害性。

    (3) 通過定量分析白格殘留體滑坡穩(wěn)定性與堵江風(fēng)險(xiǎn),驗(yàn)證了本文提出的方法可以用于碎屑流堵江的動力學(xué)過程及風(fēng)險(xiǎn)評價(jià)研究。

    參考文獻(xiàn):

    [1] 王立朝,溫銘生,馮振,等.中國西藏金沙江白格滑坡災(zāi)害研究[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),2019(1):1-9.

    [2] 馮文凱,張國強(qiáng),白慧林,等.金沙江“10·11”白格特大型滑坡形成機(jī)制及發(fā)展趨勢初步分析[J].工程地質(zhì)學(xué)報(bào),2019(2):415-425

    [3] 鄧建輝,高云建,余志球,等.堰塞金沙江上游的白格滑坡形成機(jī)制與過程分析[J].工程科學(xué)與技術(shù),2019,51(1):9-16.

    [4] 王塞,鄧建輝,陳菲,等.“10·10”西藏白格滑坡運(yùn)動特征反演分析[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào),2022,54(2):162-170.

    [5] 羅再均,次仁卓瑪,德吉白姆.2018年西藏堰塞湖搶險(xiǎn)救災(zāi)實(shí)踐及思考[J].水利水電快報(bào),2019,40(3):18-21.

    [6] 謝忱,陳骎,侯奇東,等.金沙江白格堰塞壩自然泄流沖刷潰決過程數(shù)值模擬[J].人民長江,2021,52(8):22-29.

    [7] 曹水合,吳新明,鐘東.基于Massflow的金沙江白格滑坡數(shù)值模擬及堵江風(fēng)險(xiǎn)預(yù)測[J].地質(zhì)災(zāi)害與環(huán)境保護(hù),2021(4):3-7.

    [8] 趙程,范宣梅,楊帆,等.金沙江白格滑坡運(yùn)動過程分析及潛在不穩(wěn)定巖體預(yù)測[J].科學(xué)技術(shù)與工程,2020,20(10):3860-3867.

    [9] 周禮,范宣梅,許強(qiáng),等.金沙江白格滑坡運(yùn)動過程特征數(shù)值模擬與危險(xiǎn)性預(yù)測研究[J].工程地質(zhì)學(xué)報(bào),2019,27(6):1395-1404.

    [10] 蔡耀軍,欒約生,楊啟貴,等.金沙江白格堰塞體結(jié)構(gòu)形態(tài)與潰決特征研究[J].人民長江,2019(3):15-22.

    [11] TANG Y,JIANG Q,ZHOU C.A Lagrangian-based SPH-DEM model for fluid-solid interaction with free surface flow in two dimensions[J].Applied Mathematical Modelling,2018,62(10):436-460.

    [12] 徐文杰.滑坡涌浪流–固耦合分析方法與應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2020,39(7):1420-1433.

    (編輯:江 燾,高小雲(yún))

    Abstract:

    In order to assess the stability and risk assessmentof river blockage for potentially unstable blocks (residual bodies) formed after two super-large landslides and river blockage occurred in Baige village,Jinsha River in 2018,the strict three-dimensional limit equilibrium method was applied to analyze the stability status of the residual bodies.Combing with the emergency response plans,the SPH-DEM fluid-solid coupling model was deployed to predict the deposit range and thickness,and to assess the risk of river blockage quantitatively.The results showed that:① The K1-I and K2-I blocks were in a critical state,with relatively large instability,leading to a high possibility of blocking Jinsha River again.② The clearance of river could decrease the overall height of the landslide damand increase the accumulation space for landslide debris consequently,which could effectively reduce the risk of river blockage.The research result could provide references for predicting the risk of landslide-caused river blockage.

    Key words:

    Baige residual body; river blockage; SPH-DEM; numerical simulation; risk assessment

    噜噜噜噜噜久久久久久91| 欧美成人一区二区免费高清观看| 亚洲人与动物交配视频| 日韩视频在线欧美| 亚洲成人av在线免费| 国产亚洲5aaaaa淫片| 久久久久视频综合| 国产老妇伦熟女老妇高清| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 国产爱豆传媒在线观看| 午夜激情久久久久久久| 国产精品一区www在线观看| 国产乱人偷精品视频| 国产在线免费精品| 最新中文字幕久久久久| 久久99热这里只频精品6学生| 亚洲精品成人av观看孕妇| 大又大粗又爽又黄少妇毛片口| 人体艺术视频欧美日本| 久久久久视频综合| 久久久久精品性色| 成人午夜精彩视频在线观看| 亚洲精品自拍成人| 欧美性感艳星| 男人狂女人下面高潮的视频| 国产精品一区二区三区四区免费观看| 国内揄拍国产精品人妻在线| av国产久精品久网站免费入址| 亚洲成人一二三区av| 18+在线观看网站| 亚洲av欧美aⅴ国产| 亚洲av福利一区| 日本一二三区视频观看| 国产淫语在线视频| 男女边摸边吃奶| 男女边吃奶边做爰视频| 亚洲欧美日韩另类电影网站 | 欧美最新免费一区二区三区| 我的女老师完整版在线观看| 久久国内精品自在自线图片| 一级毛片aaaaaa免费看小| 制服丝袜香蕉在线| 国产精品三级大全| 人妻系列 视频| 成年女人在线观看亚洲视频| 国产极品天堂在线| 搡老乐熟女国产| 免费黄色在线免费观看| 国产成人一区二区在线| 少妇丰满av| 久久 成人 亚洲| 欧美 日韩 精品 国产| 国产探花极品一区二区| 亚洲美女视频黄频| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 成人毛片60女人毛片免费| 777米奇影视久久| 老司机影院成人| 国产精品麻豆人妻色哟哟久久| 精品国产露脸久久av麻豆| 亚洲国产精品专区欧美| 亚洲国产最新在线播放| 大又大粗又爽又黄少妇毛片口| 日韩欧美 国产精品| 午夜激情久久久久久久| 天堂俺去俺来也www色官网| 欧美高清成人免费视频www| 日日啪夜夜撸| 激情 狠狠 欧美| 三级国产精品片| 欧美日本视频| 国产黄频视频在线观看| 美女福利国产在线 | 国产极品天堂在线| 在线观看免费高清a一片| 内地一区二区视频在线| 永久网站在线| 丰满乱子伦码专区| 如何舔出高潮| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 国产一区二区在线观看日韩| 国产深夜福利视频在线观看| 亚洲av不卡在线观看| 久久久久久人妻| 亚洲国产毛片av蜜桃av| 人妻 亚洲 视频| 国产成人午夜福利电影在线观看| 亚洲无线观看免费| 精华霜和精华液先用哪个| av在线播放精品| 香蕉精品网在线| 国产成人精品福利久久| 欧美精品国产亚洲| 麻豆精品久久久久久蜜桃| 尾随美女入室| 国产精品成人在线| 超碰av人人做人人爽久久| 亚洲av综合色区一区| 亚洲欧美日韩东京热| 最新中文字幕久久久久| 日韩免费高清中文字幕av| 在线观看一区二区三区激情| 成人毛片a级毛片在线播放| 欧美另类一区| 欧美极品一区二区三区四区| 精品熟女少妇av免费看| 久久久久人妻精品一区果冻| 高清在线视频一区二区三区| 国产黄色免费在线视频| 久久国产亚洲av麻豆专区| av黄色大香蕉| 22中文网久久字幕| 91在线精品国自产拍蜜月| av线在线观看网站| 在线观看av片永久免费下载| 美女国产视频在线观看| 免费黄频网站在线观看国产| 伊人久久国产一区二区| 日韩国内少妇激情av| 一级毛片黄色毛片免费观看视频| 亚洲国产av新网站| 高清不卡的av网站| 亚洲人与动物交配视频| 精品久久久久久久久av| 亚洲成人中文字幕在线播放| 国产视频首页在线观看| 一本一本综合久久| 亚洲不卡免费看| 高清视频免费观看一区二区| 欧美精品人与动牲交sv欧美| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 亚洲欧美中文字幕日韩二区| 久久久a久久爽久久v久久| 国产精品.久久久| av播播在线观看一区| av一本久久久久| 国产精品久久久久久精品电影小说 | 亚洲av电影在线观看一区二区三区| 日韩亚洲欧美综合| 如何舔出高潮| tube8黄色片| 人人妻人人看人人澡| 99精国产麻豆久久婷婷| 午夜福利在线在线| 中国三级夫妇交换| 免费观看无遮挡的男女| 亚洲精品456在线播放app| 熟女av电影| av在线观看视频网站免费| 色视频在线一区二区三区| 成年免费大片在线观看| 日韩在线高清观看一区二区三区| 黄片无遮挡物在线观看| 免费久久久久久久精品成人欧美视频 | 欧美日韩在线观看h| 黑人猛操日本美女一级片| 我要看黄色一级片免费的| 免费黄色在线免费观看| 中文字幕久久专区| 欧美日韩综合久久久久久| 国产精品免费大片| 亚洲人成网站在线播| 高清午夜精品一区二区三区| 青青草视频在线视频观看| 九九爱精品视频在线观看| 精品熟女少妇av免费看| 亚洲精品久久午夜乱码| 亚洲国产毛片av蜜桃av| 免费高清在线观看视频在线观看| 亚洲国产最新在线播放| 99re6热这里在线精品视频| 久久久久久久久大av| 丝袜脚勾引网站| 国产精品99久久久久久久久| 国产亚洲91精品色在线| 中文字幕制服av| 午夜福利在线观看免费完整高清在| 我要看日韩黄色一级片| 国产日韩欧美在线精品| 久久国产乱子免费精品| 国产片特级美女逼逼视频| 人妻一区二区av| 色婷婷av一区二区三区视频| 欧美丝袜亚洲另类| 你懂的网址亚洲精品在线观看| 丰满乱子伦码专区| 在线观看人妻少妇| 少妇人妻一区二区三区视频| 国产精品99久久99久久久不卡 | 欧美老熟妇乱子伦牲交| 各种免费的搞黄视频| 久久精品国产亚洲网站| 精品熟女少妇av免费看| 91精品国产国语对白视频| 亚洲第一av免费看| 搡女人真爽免费视频火全软件| 久久久久精品性色| 人体艺术视频欧美日本| 国产黄片视频在线免费观看| 日韩视频在线欧美| 亚洲国产精品专区欧美| 极品教师在线视频| 男女啪啪激烈高潮av片| 夜夜骑夜夜射夜夜干| 欧美老熟妇乱子伦牲交| 一级毛片 在线播放| 久久久久国产精品人妻一区二区| 日韩电影二区| 免费高清在线观看视频在线观看| 亚洲综合精品二区| 五月开心婷婷网| 亚洲国产欧美人成| 精品亚洲成国产av| 欧美精品一区二区免费开放| 国产 精品1| 国产黄片美女视频| 美女xxoo啪啪120秒动态图| 高清在线视频一区二区三区| 国产永久视频网站| 国模一区二区三区四区视频| 伦精品一区二区三区| 国产精品av视频在线免费观看| 国产精品一区www在线观看| 国产精品一区二区三区四区免费观看| 丰满迷人的少妇在线观看| av在线app专区| 国产爱豆传媒在线观看| 春色校园在线视频观看| 99国产精品免费福利视频| 啦啦啦在线观看免费高清www| 国产v大片淫在线免费观看| 亚洲av日韩在线播放| 一区二区三区四区激情视频| 亚洲无线观看免费| 成人特级av手机在线观看| 国产91av在线免费观看| 国产一区二区三区av在线| 热99国产精品久久久久久7| 不卡视频在线观看欧美| 国产亚洲欧美精品永久| 国产乱人偷精品视频| 日本与韩国留学比较| 男人和女人高潮做爰伦理| 日本黄色日本黄色录像| 国产男女内射视频| videos熟女内射| 亚洲精品日韩av片在线观看| 亚洲精品第二区| 亚洲国产高清在线一区二区三| 久久精品久久精品一区二区三区| 亚洲精品一二三| 五月玫瑰六月丁香| 在线观看人妻少妇| 亚洲国产精品专区欧美| 欧美高清性xxxxhd video| 亚洲精品自拍成人| 国产精品久久久久久精品电影小说 | 亚洲欧洲日产国产| 亚洲国产毛片av蜜桃av| 亚洲精品aⅴ在线观看| 十分钟在线观看高清视频www | 最近中文字幕高清免费大全6| 欧美成人午夜免费资源| 成人免费观看视频高清| 一级毛片 在线播放| 久久人人爽人人片av| 亚洲成色77777| 亚洲精品国产成人久久av| 欧美丝袜亚洲另类| 久久国产精品大桥未久av | 精品一区在线观看国产| 麻豆成人午夜福利视频| 亚洲欧美成人精品一区二区| 国产精品成人在线| 久久人人爽人人片av| 亚洲欧美一区二区三区黑人 | 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 久久精品夜色国产| 身体一侧抽搐| 热99国产精品久久久久久7| 亚洲成色77777| 超碰av人人做人人爽久久| 成人一区二区视频在线观看| 国产精品人妻久久久影院| 精品99又大又爽又粗少妇毛片| 久久99精品国语久久久| 少妇 在线观看| 身体一侧抽搐| 毛片女人毛片| 久久久久久久久久久免费av| 久久国产亚洲av麻豆专区| 亚洲性久久影院| 97在线视频观看| 国产 一区 欧美 日韩| 日韩中字成人| 特大巨黑吊av在线直播| 哪个播放器可以免费观看大片| 少妇人妻精品综合一区二区| 在线观看av片永久免费下载| 简卡轻食公司| 男的添女的下面高潮视频| 亚洲精品国产av成人精品| 国产白丝娇喘喷水9色精品| 深爱激情五月婷婷| 欧美日本视频| 欧美性感艳星| 人妻一区二区av| 91久久精品国产一区二区三区| 美女中出高潮动态图| 七月丁香在线播放| 伊人久久国产一区二区| 亚洲av.av天堂| 中文字幕免费在线视频6| 亚洲熟女精品中文字幕| 免费看av在线观看网站| 国产精品成人在线| 成人二区视频| 日本午夜av视频| 国产精品国产三级国产专区5o| 精品人妻偷拍中文字幕| 乱系列少妇在线播放| 欧美三级亚洲精品| 美女cb高潮喷水在线观看| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 亚洲国产高清在线一区二区三| 欧美成人一区二区免费高清观看| 人体艺术视频欧美日本| 我要看黄色一级片免费的| 欧美日韩综合久久久久久| 六月丁香七月| 美女脱内裤让男人舔精品视频| 亚洲高清免费不卡视频| 大陆偷拍与自拍| 久久午夜福利片| h日本视频在线播放| av国产精品久久久久影院| 人人妻人人爽人人添夜夜欢视频 | av国产久精品久网站免费入址| 成人免费观看视频高清| videossex国产| 亚洲成色77777| 久久韩国三级中文字幕| 黄色日韩在线| 97热精品久久久久久| 久久女婷五月综合色啪小说| 18禁裸乳无遮挡免费网站照片| 亚洲成色77777| 毛片一级片免费看久久久久| 久久女婷五月综合色啪小说| 能在线免费看毛片的网站| 亚洲经典国产精华液单| 一区二区三区乱码不卡18| 一个人看的www免费观看视频| 欧美zozozo另类| 亚洲最大成人中文| 高清不卡的av网站| 一个人看的www免费观看视频| 在线观看免费高清a一片| 伦理电影大哥的女人| 久久 成人 亚洲| 性高湖久久久久久久久免费观看| 欧美xxxx黑人xx丫x性爽| h日本视频在线播放| 嫩草影院入口| 肉色欧美久久久久久久蜜桃| 男人和女人高潮做爰伦理| 三级国产精品欧美在线观看| 国产精品一区www在线观看| av网站免费在线观看视频| av女优亚洲男人天堂| 在线看a的网站| av女优亚洲男人天堂| 91久久精品电影网| 99热全是精品| 哪个播放器可以免费观看大片| 亚洲国产毛片av蜜桃av| 久久久久视频综合| 精品一品国产午夜福利视频| 91久久精品国产一区二区三区| 蜜桃在线观看..| 搡老乐熟女国产| 尾随美女入室| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说| 99久久精品国产国产毛片| 永久网站在线| 亚洲人与动物交配视频| 中文欧美无线码| 热re99久久精品国产66热6| 国产午夜精品久久久久久一区二区三区| 免费看av在线观看网站| 视频区图区小说| 国产深夜福利视频在线观看| 各种免费的搞黄视频| 国产色婷婷99| 欧美少妇被猛烈插入视频| 国产成人午夜福利电影在线观看| 天美传媒精品一区二区| 女人久久www免费人成看片| videossex国产| 免费观看在线日韩| 麻豆成人午夜福利视频| 日本黄大片高清| 国产 精品1| 婷婷色麻豆天堂久久| 青春草国产在线视频| 少妇人妻一区二区三区视频| 国产深夜福利视频在线观看| 精品午夜福利在线看| 在线观看免费视频网站a站| 91久久精品国产一区二区三区| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 久久国产精品大桥未久av | 久久国产精品男人的天堂亚洲 | 亚洲中文av在线| 亚洲成人手机| 亚洲国产高清在线一区二区三| 精品一区二区三区视频在线| 亚洲真实伦在线观看| 一级av片app| 日韩不卡一区二区三区视频在线| 久久久a久久爽久久v久久| 久久久国产一区二区| 久久久欧美国产精品| 欧美日韩视频高清一区二区三区二| 黄色配什么色好看| 欧美成人一区二区免费高清观看| 两个人的视频大全免费| 下体分泌物呈黄色| av播播在线观看一区| 精品99又大又爽又粗少妇毛片| 久久久久久久久大av| 中文字幕久久专区| 欧美一区二区亚洲| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 91午夜精品亚洲一区二区三区| 色视频www国产| 免费播放大片免费观看视频在线观看| 一个人免费看片子| 丰满乱子伦码专区| 久久久久国产网址| 国产亚洲精品久久久com| 纵有疾风起免费观看全集完整版| 99热全是精品| 成人国产麻豆网| 国产 一区 欧美 日韩| 一级a做视频免费观看| 亚洲国产精品999| 卡戴珊不雅视频在线播放| 熟妇人妻不卡中文字幕| 久久精品国产自在天天线| 亚洲国产成人一精品久久久| 各种免费的搞黄视频| 国产淫语在线视频| 午夜免费鲁丝| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 街头女战士在线观看网站| 九九爱精品视频在线观看| 成人国产av品久久久| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 精品久久久噜噜| 伦理电影免费视频| 欧美激情国产日韩精品一区| 午夜福利影视在线免费观看| 亚洲av国产av综合av卡| 内地一区二区视频在线| 日韩av免费高清视频| 久久久久久久久久成人| 午夜福利影视在线免费观看| freevideosex欧美| 免费黄网站久久成人精品| 爱豆传媒免费全集在线观看| 亚洲精华国产精华液的使用体验| 欧美成人a在线观看| 1000部很黄的大片| 热re99久久精品国产66热6| www.色视频.com| 中文字幕久久专区| 久久99热这里只频精品6学生| 网址你懂的国产日韩在线| 亚洲综合精品二区| 欧美亚洲 丝袜 人妻 在线| 熟妇人妻不卡中文字幕| 日本猛色少妇xxxxx猛交久久| 国产伦理片在线播放av一区| 亚洲国产最新在线播放| freevideosex欧美| 亚洲国产精品国产精品| 精品久久国产蜜桃| 色综合色国产| 色婷婷久久久亚洲欧美| 成人黄色视频免费在线看| 免费黄频网站在线观看国产| 欧美性感艳星| 一级毛片黄色毛片免费观看视频| 国产视频内射| 视频区图区小说| 精品视频人人做人人爽| 在线观看美女被高潮喷水网站| 成人亚洲欧美一区二区av| 在线免费十八禁| 精华霜和精华液先用哪个| 又黄又爽又刺激的免费视频.| av一本久久久久| 亚洲经典国产精华液单| 狂野欧美白嫩少妇大欣赏| 在线观看一区二区三区| 成人一区二区视频在线观看| 久久6这里有精品| 日本黄大片高清| 一级av片app| 日韩中文字幕视频在线看片 | 王馨瑶露胸无遮挡在线观看| 日韩av不卡免费在线播放| 亚洲人成网站在线播| 亚洲av福利一区| 伊人久久精品亚洲午夜| 热99国产精品久久久久久7| 国产一区二区三区av在线| 亚洲国产毛片av蜜桃av| 日本av手机在线免费观看| 亚洲欧美精品自产自拍| 如何舔出高潮| 免费久久久久久久精品成人欧美视频 | 久久99蜜桃精品久久| 精品亚洲成国产av| 超碰av人人做人人爽久久| av.在线天堂| 啦啦啦啦在线视频资源| 亚州av有码| 精品少妇黑人巨大在线播放| 国产老妇伦熟女老妇高清| 3wmmmm亚洲av在线观看| 久久这里有精品视频免费| 久久6这里有精品| 中文天堂在线官网| 高清毛片免费看| 男女边吃奶边做爰视频| .国产精品久久| 亚洲精品成人av观看孕妇| 国产成人免费观看mmmm| 亚洲,欧美,日韩| av不卡在线播放| 美女主播在线视频| 日韩欧美一区视频在线观看 | 中文在线观看免费www的网站| 99九九线精品视频在线观看视频| 免费观看av网站的网址| 国产精品久久久久久久久免| 日本猛色少妇xxxxx猛交久久| 亚洲色图av天堂| 少妇人妻精品综合一区二区| 亚洲av男天堂| 亚洲久久久国产精品| 热99国产精品久久久久久7| 国产一区二区在线观看日韩| av视频免费观看在线观看| 好男人视频免费观看在线| av在线观看视频网站免费| 一级爰片在线观看| 日本爱情动作片www.在线观看| 日韩欧美 国产精品| 蜜桃久久精品国产亚洲av| 丰满乱子伦码专区| 国产av码专区亚洲av| 免费看日本二区| 色哟哟·www| 久久久久久久久久人人人人人人| 色视频www国产| 国产成人精品婷婷| 亚洲成色77777| 久久99热6这里只有精品| 777米奇影视久久| 成人免费观看视频高清| a级毛色黄片| 啦啦啦在线观看免费高清www| 女性生殖器流出的白浆| 日韩一区二区三区影片| 在线观看三级黄色| 国产黄片美女视频| 建设人人有责人人尽责人人享有的 | 美女cb高潮喷水在线观看| 亚洲精品456在线播放app| 麻豆国产97在线/欧美| 欧美最新免费一区二区三区| 国产69精品久久久久777片| 在线观看三级黄色| 亚洲成色77777| 免费看av在线观看网站| 免费久久久久久久精品成人欧美视频 | 老熟女久久久| 成人亚洲欧美一区二区av| 欧美日韩国产mv在线观看视频 | 黄色怎么调成土黄色| 欧美日韩综合久久久久久| 狂野欧美白嫩少妇大欣赏| 五月伊人婷婷丁香| 久久久久久久久久久丰满| 欧美少妇被猛烈插入视频| 午夜免费观看性视频| av国产精品久久久久影院| 99视频精品全部免费 在线| 韩国av在线不卡| 色婷婷av一区二区三区视频| 人人妻人人看人人澡|