• <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)險分析

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

    唐岳灝 姜清輝

    摘要:

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

    關(guān)鍵詞:

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

    中圖法分類號: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è)施,同時潰壩洪水嚴(yán)重沖毀滑坡下游的村莊、農(nóng)田和公路、橋梁等基礎(chǔ)設(shè)施[5-6]。白格滑坡經(jīng)兩次滑動后,松散物質(zhì)已大為減少,但滑坡體對山體擾動作用巨大(主要表現(xiàn)為拉、拽、刮、鏟、刷),滑坡三面由于臨空卸荷作用不斷加強,變形跡象十分明顯,后緣及兩側(cè)發(fā)育多條深大裂縫,存在大量滑坡殘留體。無人機航空攝影測量及工程地質(zhì)測繪、鉆探、物探等詳細(xì)勘察工作的成果表明:前兩次白格滑坡周界后緣有3個主要殘留變形體,其滑坡周界明顯,已經(jīng)發(fā)生了明顯的下滑、錯動,仍具備進一步下滑的趨勢。同時,地表宏觀變形跡象和監(jiān)測成果表明殘留體一直存在前緣溜滑和蠕滑變形,其穩(wěn)定安全問題不容樂觀。因此,開展白格滑坡殘留體的防災(zāi)減災(zāi)研究,分析殘留體的穩(wěn)定性和變形趨勢,可為滑坡應(yīng)急處置提供支撐,并為白格滑坡第三次堵江的可能性、危害性預(yù)測提供基本數(shù)據(jù)。

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

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

    1 滑坡體特征

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

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

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

    2.1 計算模型和計算方法

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

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

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

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

    采用三維嚴(yán)格極限平衡法計算得到的白格滑坡后緣殘留體K1-Ⅰ,K1-Ⅱ,K1,K2-Ⅰ,K2-Ⅱ和K3在不同工況條件下穩(wěn)定性安全系數(shù)以及各個滑塊的體積、主滑方向和穩(wěn)定狀態(tài),如表2所示,其中主滑方向角為與正北方向夾角,順時針為正。穩(wěn)定狀態(tài)根據(jù)《青藏高原重大滑坡風(fēng)險防控指南》中滑坡穩(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)險等級(圖2),為下一步進行滑坡殘留體應(yīng)急處置和堵江風(fēng)險評估提供了直觀的參考。圖5為計算得到的各滑塊主滑方向與地表外觀測點位移矢量的對比,從中可以看出,主滑方向與外觀測點變形方向吻合,進一步論證了所使用的三維嚴(yán)格分析方法的有效性和可靠性。

    基于對滑源區(qū)殘留體的穩(wěn)定性評估,可以看出K1-Ⅰ和K2-Ⅰ滑塊處于臨界狀態(tài),失穩(wěn)風(fēng)險較大,因此,首先考慮的失穩(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)險,因此,需要采取應(yīng)急處置措施。對白格滑坡后緣殘留變形體的處置原則采用滑源區(qū)—運移區(qū)—堆積區(qū)全過程處理中的“顧兩頭”,分別為滑源區(qū)的處置和金沙江河道堰塞體的處置。滑源區(qū)應(yīng)急處置采用的方案是削方減載,主要對K1-Ⅰ殘留體進行削坡,其目的是提高殘留體的穩(wěn)定性,降低滑坡風(fēng)險,如圖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)險。具體方案為將堰塞體殘體中靠近河道側(cè)不小于40.0 m范圍內(nèi)、2 920 m高程以上的殘體全部清除。開挖設(shè)計方案共有4級馬道,馬道寬度0.2 m,放坡坡比1∶2.0,每一級臺階高度20 m,總開挖體積約245萬m3,如圖7所示。

    4 滑坡殘留體堵江評估

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

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

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

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

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

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

    5 結(jié) 論

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

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

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

    參考文獻(xiàn):

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

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

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

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

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

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

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

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

    [9] 周禮,范宣梅,許強,等.金沙江白格滑坡運動過程特征數(shù)值模擬與危險性預(yù)測研究[J].工程地質(zhì)學(xué)報,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é)報,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

    亚洲美女黄片视频| 午夜福利在线观看免费完整高清在 | 中文字幕人成人乱码亚洲影| 婷婷亚洲欧美| 午夜激情福利司机影院| 观看美女的网站| 亚洲va日本ⅴa欧美va伊人久久| 国产成人av教育| 久久天堂一区二区三区四区| 亚洲人成网站高清观看| 在线观看午夜福利视频| 欧美在线黄色| 国产高清视频在线播放一区| 亚洲专区字幕在线| 天天躁日日操中文字幕| 五月伊人婷婷丁香| 一区二区三区高清视频在线| 国产精品一区二区精品视频观看| 国产精品 欧美亚洲| 成人18禁在线播放| 国产精品久久电影中文字幕| 精品一区二区三区四区五区乱码| 久9热在线精品视频| 91在线精品国自产拍蜜月 | 大型黄色视频在线免费观看| 婷婷精品国产亚洲av| 亚洲专区字幕在线| 九色国产91popny在线| 两人在一起打扑克的视频| 久久久久国内视频| a级毛片a级免费在线| 长腿黑丝高跟| 亚洲欧美日韩无卡精品| 成人三级黄色视频| 在线观看免费视频日本深夜| 亚洲午夜理论影院| 看黄色毛片网站| 啦啦啦观看免费观看视频高清| 日本撒尿小便嘘嘘汇集6| 亚洲五月婷婷丁香| 精品免费久久久久久久清纯| 九九热线精品视视频播放| 国产精品久久久人人做人人爽| 免费在线观看影片大全网站| 亚洲精品一卡2卡三卡4卡5卡| 9191精品国产免费久久| 国产精品 欧美亚洲| 色视频www国产| 18禁裸乳无遮挡免费网站照片| 欧美在线一区亚洲| 欧美日韩乱码在线| 叶爱在线成人免费视频播放| 国产成人精品久久二区二区91| 天堂√8在线中文| 国产成年人精品一区二区| av中文乱码字幕在线| 午夜福利在线在线| av天堂中文字幕网| 亚洲欧美精品综合一区二区三区| 午夜精品一区二区三区免费看| 久久久国产成人精品二区| 国产精品精品国产色婷婷| 每晚都被弄得嗷嗷叫到高潮| 香蕉久久夜色| 99热精品在线国产| 久久国产精品人妻蜜桃| 特大巨黑吊av在线直播| 精品久久蜜臀av无| 两性夫妻黄色片| cao死你这个sao货| 欧美性猛交╳xxx乱大交人| 国产免费av片在线观看野外av| 亚洲成人精品中文字幕电影| 一个人看的www免费观看视频| 啦啦啦观看免费观看视频高清| 久久国产精品人妻蜜桃| 亚洲欧美激情综合另类| 免费搜索国产男女视频| 欧美一区二区国产精品久久精品| 男女做爰动态图高潮gif福利片| 国产野战对白在线观看| 999久久久精品免费观看国产| 又黄又爽又免费观看的视频| 九色国产91popny在线| 99久久综合精品五月天人人| 欧美日韩福利视频一区二区| 久久午夜亚洲精品久久| 一个人免费在线观看电影 | 国产成人一区二区三区免费视频网站| 国产真实乱freesex| 亚洲av日韩精品久久久久久密| 亚洲精品色激情综合| 欧美在线黄色| 亚洲真实伦在线观看| 麻豆成人av在线观看| 亚洲成av人片在线播放无| 美女被艹到高潮喷水动态| 日韩欧美国产在线观看| 国产三级中文精品| 看黄色毛片网站| 日韩大尺度精品在线看网址| 高清在线国产一区| 性欧美人与动物交配| 欧美丝袜亚洲另类 | 国产三级黄色录像| 亚洲精品国产精品久久久不卡| 午夜免费观看网址| 91av网一区二区| 亚洲国产精品999在线| 91字幕亚洲| 2021天堂中文幕一二区在线观| 2021天堂中文幕一二区在线观| 成人av一区二区三区在线看| 国产精品久久久久久久电影 | 成人精品一区二区免费| 日韩欧美在线乱码| 美女被艹到高潮喷水动态| 成人永久免费在线观看视频| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区免费观看 | 又紧又爽又黄一区二区| 午夜福利高清视频| 夜夜看夜夜爽夜夜摸| 国产爱豆传媒在线观看| 男人舔女人的私密视频| av天堂在线播放| 91字幕亚洲| 日本五十路高清| 久久精品夜夜夜夜夜久久蜜豆| 人人妻,人人澡人人爽秒播| 亚洲专区中文字幕在线| 我的老师免费观看完整版| 日韩欧美 国产精品| 色哟哟哟哟哟哟| 小说图片视频综合网站| 亚洲熟女毛片儿| 国产精品98久久久久久宅男小说| 日本成人三级电影网站| 亚洲欧美精品综合一区二区三区| 日本黄大片高清| av天堂中文字幕网| 日韩欧美国产一区二区入口| 亚洲欧美日韩无卡精品| 亚洲av成人不卡在线观看播放网| 麻豆久久精品国产亚洲av| 999精品在线视频| 欧美日韩综合久久久久久 | 色视频www国产| 色尼玛亚洲综合影院| 99久国产av精品| 老司机午夜福利在线观看视频| 中文字幕久久专区| 亚洲国产欧美网| 亚洲av成人一区二区三| 国产亚洲精品综合一区在线观看| 亚洲av成人一区二区三| 无限看片的www在线观看| 免费看光身美女| 国产成+人综合+亚洲专区| 欧美日本视频| 亚洲国产日韩欧美精品在线观看 | 免费观看精品视频网站| 韩国av一区二区三区四区| 久久久久久久久免费视频了| 看免费av毛片| 日韩大尺度精品在线看网址| 老司机深夜福利视频在线观看| 国产精品亚洲美女久久久| 超碰成人久久| 亚洲一区二区三区不卡视频| h日本视频在线播放| 黄色日韩在线| 国产午夜福利久久久久久| 亚洲一区二区三区色噜噜| 国产野战对白在线观看| 国产精品一区二区三区四区久久| 成年女人永久免费观看视频| 国产一区二区激情短视频| 一夜夜www| 欧美日韩福利视频一区二区| 黄色片一级片一级黄色片| 国产单亲对白刺激| 蜜桃久久精品国产亚洲av| 99热这里只有是精品50| 成人欧美大片| 好看av亚洲va欧美ⅴa在| 真实男女啪啪啪动态图| 午夜激情福利司机影院| 桃色一区二区三区在线观看| 欧美高清成人免费视频www| 又紧又爽又黄一区二区| 日本黄色视频三级网站网址| 少妇人妻一区二区三区视频| 男女下面进入的视频免费午夜| 丰满的人妻完整版| 国产一区二区三区在线臀色熟女| 免费观看人在逋| 亚洲黑人精品在线| 久久久久九九精品影院| bbb黄色大片| 成年女人毛片免费观看观看9| 欧美国产日韩亚洲一区| 高清在线国产一区| 搡老熟女国产l中国老女人| 不卡av一区二区三区| 99视频精品全部免费 在线 | 国产成人精品无人区| 日本一二三区视频观看| 亚洲国产色片| 欧美国产日韩亚洲一区| 国产在线精品亚洲第一网站| 午夜日韩欧美国产| av福利片在线观看| 色哟哟哟哟哟哟| 久久这里只有精品中国| 亚洲精品色激情综合| 精品久久久久久久末码| 国产欧美日韩一区二区精品| 精品久久蜜臀av无| 一本精品99久久精品77| 亚洲无线观看免费| 亚洲午夜精品一区,二区,三区| 日本黄色视频三级网站网址| 嫩草影视91久久| 精品乱码久久久久久99久播| 两个人视频免费观看高清| 亚洲五月婷婷丁香| 色av中文字幕| 99热这里只有是精品50| 手机成人av网站| 国产三级中文精品| 婷婷六月久久综合丁香| 久久中文字幕人妻熟女| 黄色片一级片一级黄色片| 国产精品电影一区二区三区| 亚洲国产精品999在线| 免费搜索国产男女视频| av黄色大香蕉| 亚洲午夜精品一区,二区,三区| 他把我摸到了高潮在线观看| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华精| 美女 人体艺术 gogo| 法律面前人人平等表现在哪些方面| 超碰成人久久| 国产99白浆流出| 黄色女人牲交| 成年女人永久免费观看视频| 在线观看日韩欧美| 中文字幕最新亚洲高清| 亚洲午夜精品一区,二区,三区| 精品电影一区二区在线| 人人妻人人澡欧美一区二区| 欧美黑人欧美精品刺激| 一个人免费在线观看电影 | 夜夜看夜夜爽夜夜摸| 97人妻精品一区二区三区麻豆| 丰满人妻一区二区三区视频av | 91字幕亚洲| 男人舔女人下体高潮全视频| www.999成人在线观看| a级毛片在线看网站| 99久久国产精品久久久| 亚洲av成人精品一区久久| 又黄又粗又硬又大视频| 女人高潮潮喷娇喘18禁视频| 丰满人妻一区二区三区视频av | 国产精品久久久久久亚洲av鲁大| 久久99热这里只有精品18| 噜噜噜噜噜久久久久久91| 国产成人影院久久av| 99久国产av精品| 亚洲午夜理论影院| 一进一出抽搐gif免费好疼| 国产精品香港三级国产av潘金莲| 极品教师在线免费播放| 亚洲色图 男人天堂 中文字幕| 久久久国产成人精品二区| 亚洲午夜精品一区,二区,三区| 亚洲国产精品sss在线观看| 亚洲国产欧洲综合997久久,| 男人的好看免费观看在线视频| 国产精品亚洲av一区麻豆| 黑人操中国人逼视频| 精品一区二区三区视频在线观看免费| 在线永久观看黄色视频| 日日干狠狠操夜夜爽| 亚洲精品456在线播放app | 美女cb高潮喷水在线观看 | 亚洲国产精品999在线| 床上黄色一级片| 一本久久中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲国产日韩欧美精品在线观看 | 女人高潮潮喷娇喘18禁视频| 免费大片18禁| 国产av一区在线观看免费| 不卡一级毛片| aaaaa片日本免费| 熟妇人妻久久中文字幕3abv| 国产免费av片在线观看野外av| 免费观看的影片在线观看| 一区二区三区国产精品乱码| 亚洲国产精品999在线| 成人无遮挡网站| 国产精品一区二区免费欧美| 999久久久国产精品视频| 国产高清有码在线观看视频| 国产人伦9x9x在线观看| 成年女人看的毛片在线观看| 国产精品九九99| 三级毛片av免费| 中文资源天堂在线| 色在线成人网| 制服丝袜大香蕉在线| 国产日本99.免费观看| 天堂av国产一区二区熟女人妻| 久久精品91无色码中文字幕| 美女高潮的动态| 中文字幕高清在线视频| 午夜福利高清视频| 岛国视频午夜一区免费看| 久久久色成人| 中文字幕av在线有码专区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久久久人妻精品电影| 欧美中文综合在线视频| 国产精品亚洲一级av第二区| 久久久水蜜桃国产精品网| 久久久久性生活片| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆 | 色综合欧美亚洲国产小说| 老司机在亚洲福利影院| 国产精品久久久人人做人人爽| 午夜a级毛片| a级毛片在线看网站| 成人亚洲精品av一区二区| 又黄又爽又免费观看的视频| 国产精品久久久久久精品电影| 高清毛片免费观看视频网站| 99热精品在线国产| 欧美日本亚洲视频在线播放| 午夜福利在线观看吧| 中文字幕久久专区| 人人妻,人人澡人人爽秒播| 欧美在线一区亚洲| 久久人妻av系列| 成人av一区二区三区在线看| 久久久精品欧美日韩精品| 桃色一区二区三区在线观看| 一二三四在线观看免费中文在| 亚洲18禁久久av| 狂野欧美白嫩少妇大欣赏| www日本在线高清视频| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 亚洲一区二区三区不卡视频| 国产伦精品一区二区三区视频9 | 国产99白浆流出| 黄片小视频在线播放| 成人无遮挡网站| 精品国内亚洲2022精品成人| 国产av在哪里看| 久久久久久久精品吃奶| 精品久久久久久久久久免费视频| 99国产极品粉嫩在线观看| 91在线精品国自产拍蜜月 | 色视频www国产| 黑人巨大精品欧美一区二区mp4| 亚洲精品美女久久av网站| 操出白浆在线播放| 欧美精品啪啪一区二区三区| 久久国产乱子伦精品免费另类| 99在线人妻在线中文字幕| 欧美另类亚洲清纯唯美| 国产亚洲av高清不卡| 国产精品一及| 色精品久久人妻99蜜桃| 黄色日韩在线| 国产探花在线观看一区二区| 国产精品久久久av美女十八| 欧美色欧美亚洲另类二区| 国产免费男女视频| 久久久久亚洲av毛片大全| 亚洲精品中文字幕一二三四区| 成人性生交大片免费视频hd| 别揉我奶头~嗯~啊~动态视频| 亚洲真实伦在线观看| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 色综合婷婷激情| 村上凉子中文字幕在线| 欧美日韩乱码在线| 1024手机看黄色片| 九色国产91popny在线| 国产乱人伦免费视频| 午夜a级毛片| 性色avwww在线观看| 亚洲在线观看片| 亚洲av成人不卡在线观看播放网| 黄频高清免费视频| 法律面前人人平等表现在哪些方面| 成人高潮视频无遮挡免费网站| 嫩草影视91久久| 亚洲av日韩精品久久久久久密| 欧美3d第一页| 欧美精品啪啪一区二区三区| 亚洲成人久久性| 国产伦精品一区二区三区视频9 | 亚洲人成网站在线播放欧美日韩| 人人妻,人人澡人人爽秒播| 他把我摸到了高潮在线观看| 国内毛片毛片毛片毛片毛片| 18禁裸乳无遮挡免费网站照片| 黄色片一级片一级黄色片| 久久久久精品国产欧美久久久| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 叶爱在线成人免费视频播放| 操出白浆在线播放| 久久久久久久久免费视频了| 波多野结衣高清作品| 夜夜躁狠狠躁天天躁| 天堂动漫精品| 国产精品野战在线观看| 久久久精品欧美日韩精品| 欧美日本视频| 性色av乱码一区二区三区2| 两性午夜刺激爽爽歪歪视频在线观看| 91在线精品国自产拍蜜月 | 一个人看视频在线观看www免费 | 日本黄色片子视频| 亚洲成人中文字幕在线播放| 色尼玛亚洲综合影院| 丰满的人妻完整版| 亚洲av成人精品一区久久| 五月玫瑰六月丁香| 天堂av国产一区二区熟女人妻| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人一区二区三| 日韩欧美免费精品| 国产精品一区二区三区四区免费观看 | 一个人免费在线观看电影 | 亚洲五月婷婷丁香| 1000部很黄的大片| 三级毛片av免费| 亚洲精品在线观看二区| 成人三级做爰电影| 免费av毛片视频| 久久热在线av| 国产午夜精品久久久久久| 99国产精品99久久久久| 极品教师在线免费播放| 国产精品乱码一区二三区的特点| 天天躁日日操中文字幕| 全区人妻精品视频| 国产真实乱freesex| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 五月伊人婷婷丁香| 99久久成人亚洲精品观看| 欧美黄色片欧美黄色片| 成人欧美大片| 叶爱在线成人免费视频播放| 亚洲色图av天堂| 亚洲精品久久国产高清桃花| 国内精品美女久久久久久| 在线观看免费午夜福利视频| 人人妻人人澡欧美一区二区| 成人亚洲精品av一区二区| 婷婷六月久久综合丁香| 国产成人av激情在线播放| 国产三级黄色录像| 成人永久免费在线观看视频| 母亲3免费完整高清在线观看| 国产亚洲精品久久久久久毛片| tocl精华| av天堂在线播放| 久久久久国产一级毛片高清牌| 精品国产三级普通话版| 久久热在线av| 中文字幕熟女人妻在线| 又黄又爽又免费观看的视频| 国产亚洲精品久久久com| ponron亚洲| 99久久精品国产亚洲精品| 精华霜和精华液先用哪个| svipshipincom国产片| 88av欧美| 日韩欧美在线乱码| 亚洲欧美日韩卡通动漫| 国产一区二区三区视频了| 国产麻豆成人av免费视频| av中文乱码字幕在线| 色综合欧美亚洲国产小说| 国产免费av片在线观看野外av| 99热这里只有是精品50| 午夜福利免费观看在线| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 欧美+亚洲+日韩+国产| 精品日产1卡2卡| 国产私拍福利视频在线观看| 国产激情偷乱视频一区二区| 搡老妇女老女人老熟妇| 九九在线视频观看精品| 国内精品久久久久久久电影| 在线观看一区二区三区| 精品久久久久久久久久免费视频| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 男人和女人高潮做爰伦理| 国产精品久久久久久亚洲av鲁大| 成人亚洲精品av一区二区| 中文字幕人妻丝袜一区二区| 欧美日本视频| 亚洲精品久久国产高清桃花| 欧美又色又爽又黄视频| 亚洲男人的天堂狠狠| 99精品欧美一区二区三区四区| 色综合婷婷激情| 欧美黑人巨大hd| 黄色 视频免费看| 最新在线观看一区二区三区| 日韩欧美一区二区三区在线观看| 99精品在免费线老司机午夜| 搡老妇女老女人老熟妇| 精品一区二区三区视频在线观看免费| aaaaa片日本免费| 精品国产三级普通话版| 99久久精品一区二区三区| 日韩欧美免费精品| 动漫黄色视频在线观看| 国产单亲对白刺激| 丰满人妻一区二区三区视频av | 激情在线观看视频在线高清| 香蕉国产在线看| 精品久久久久久,| 亚洲欧洲精品一区二区精品久久久| 19禁男女啪啪无遮挡网站| 欧美绝顶高潮抽搐喷水| 日本一二三区视频观看| 久久久久亚洲av毛片大全| 国内久久婷婷六月综合欲色啪| 色综合站精品国产| 一二三四在线观看免费中文在| 欧美日韩一级在线毛片| 巨乳人妻的诱惑在线观看| 国产精品一区二区免费欧美| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合一区二区三区| 久99久视频精品免费| 一区福利在线观看| 国产精品电影一区二区三区| 狠狠狠狠99中文字幕| 国产极品精品免费视频能看的| 成年免费大片在线观看| 亚洲激情在线av| 亚洲av片天天在线观看| 最近最新免费中文字幕在线| 18禁观看日本| 亚洲国产中文字幕在线视频| 欧美性猛交黑人性爽| 老司机在亚洲福利影院| 99热6这里只有精品| 九九久久精品国产亚洲av麻豆 | 久久中文看片网| 国产亚洲精品av在线| 欧美成人一区二区免费高清观看 | 国产久久久一区二区三区| 久久久国产精品麻豆| 一个人看的www免费观看视频| 可以在线观看的亚洲视频| 男女下面进入的视频免费午夜| 真人做人爱边吃奶动态| 欧美黑人巨大hd| 精品人妻1区二区| 天堂影院成人在线观看| 精品国产乱码久久久久久男人| 亚洲一区二区三区色噜噜| 人妻久久中文字幕网| 哪里可以看免费的av片| h日本视频在线播放| 狂野欧美白嫩少妇大欣赏| 免费av毛片视频| 国产成+人综合+亚洲专区| 久久性视频一级片| 国产真人三级小视频在线观看| 国产1区2区3区精品| 日本 欧美在线| 人人妻人人澡欧美一区二区| 国产成人av激情在线播放| cao死你这个sao货| 亚洲专区中文字幕在线| svipshipincom国产片| av在线天堂中文字幕| 日本黄大片高清| 色综合欧美亚洲国产小说| 色综合站精品国产| 午夜两性在线视频| 中出人妻视频一区二区| 床上黄色一级片| 婷婷精品国产亚洲av| 亚洲欧美激情综合另类| 国产高清视频在线观看网站| 久久久久国产精品人妻aⅴ院| 欧美日韩一级在线毛片| 欧美午夜高清在线| 香蕉久久夜色| 麻豆成人午夜福利视频| 国产精品久久久人人做人人爽| 国产亚洲精品av在线| 1024手机看黄色片| 久久久成人免费电影|