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

    基于精細DEM 的崩塌滑坡災害識別及主控因素分析
    ——以雅魯藏布江縫合帶加查—朗縣段為例*

    2019-10-26 07:07:26王瑞琪王學良劉海洋孫娟娟王新輝
    工程地質(zhì)學報 2019年5期
    關鍵詞:面密度坡向柵格

    王瑞琪 王學良 劉海洋 孫娟娟 王新輝 張 蘇

    ( ①中國科學院地質(zhì)與地球物理研究所,中國科學院頁巖氣與地質(zhì)工程重點實驗室 北京 100029)

    ( ②中國科學院大學 北京 100049)

    ( ③中國科學院地球科學研究院 北京 100029)

    ( ④北京洛斯達數(shù)字遙感技術有限公司 北京 100120)

    ( ⑤國網(wǎng)經(jīng)濟技術研究院有限公司 北京 102209)

    0 引 言

    雅魯藏布江縫合帶加查—朗縣段位于青藏高原東南部,崩塌、滑坡等地質(zhì)災害十分發(fā)育,伴隨著青藏高原地區(qū)人類工程的日益增多,該區(qū)域的地質(zhì)災害識別、發(fā)育特征與成因機理研究迫在眉睫,以求為后續(xù)的風險評估工作打下基礎。

    近年來,針對雅魯藏布江縫合帶的工程地質(zhì)條件及地質(zhì)災害分布情況,有許多學者展開了調(diào)查研究。如張路青等( 2004) 對川藏公路南線八宿—林芝段滾石災害進行了工程地質(zhì)調(diào)查,給出了初步的防護建議; 袁廣祥等( 2010) 對雅魯藏布江大拐彎北部川藏公路地質(zhì)災害發(fā)育與分布開展了研究; 陳洪凱等( 2011) 運用層次分析法分析了川藏公路沿線的地質(zhì)災害的危險性,并建立了相應的評價模型和評價數(shù)據(jù)庫; 吳瑞安等( 2017) 對川藏鐵路加查—朗縣段幾處重大滑坡災害進行了詳細闡述; 李秀珍等( 2017) 結(jié)合現(xiàn)場調(diào)查和遙感解譯等手段,研究了川藏交通廊道康定至林芝段滑坡崩塌的空間分布特征和危害方式。限于數(shù)據(jù)精度和工作程度等方面的原因,如何在研究區(qū)開展災害面線上的精細識別和主控因素分析尚顯不足。本文利用機載雷達獲取了精細DEM( 10 m 分辨率) ,UAV 獲取高清航攝數(shù)據(jù),在已有的崩塌滑坡災害遙感目視解譯識別方法的基礎上,提出了區(qū)域尺度上基于統(tǒng)計學方法的崩塌滑坡災害的精細識別方法; 開展了研究區(qū)崩塌滑坡主控因素分析,尤其是各地形地貌因子對崩塌滑坡的影響程度,并進行了可靠性驗證。

    1 地質(zhì)背景條件

    加查—朗縣段位于藏東南地區(qū),雅魯藏布江中上游地帶,區(qū)內(nèi)河流下切與側(cè)蝕作用強烈,山高谷深,斜坡陡立,地形起伏大,坡體植被發(fā)育較少,巖體裸露。結(jié)合區(qū)域地質(zhì)資料,區(qū)內(nèi)分布有3 條與雅魯藏布江縫合帶大致平行的大型伴生斷裂( 圖1) ,新構造運動強烈,斷裂帶附近巖體破碎。區(qū)內(nèi)巖性組成較為簡單,主要出露千枚巖、礫巖和花崗巖。研究區(qū)主要位于雅魯藏布江兩岸,高程在3060~5330 m 之間,主要集中在3100~4400 m 的高程范圍。在內(nèi)外動力地質(zhì)作用下,研究區(qū)容易發(fā)生崩塌、滑坡等地質(zhì)災害。

    2 研究區(qū)地質(zhì)災害高精度識別

    2.1 崩塌、滑坡點識別

    研究區(qū)的崩塌、滑坡災害的識別技術主要采取遙感技術,包括衛(wèi)星遙感和無人機遙感技術( Mckean et al.,1991; Mantovani et al.,1996; Akgun et al.,2012; Scaioni et al.,2014; 黃海峰等,2017; 許強等,2018; 唐堯等,2019) 。數(shù)據(jù)來源主要為機載雷達獲取的精細DEM 數(shù)據(jù)( 10 m 精度) 、高清航拍數(shù)據(jù),以結(jié)合衛(wèi)星遙感影像作為輔助數(shù)據(jù)。根據(jù)王學良( 2018) ,本文采用的歷史崩塌滑坡識別方法步驟如下:( 1) 建立研究區(qū)的崩塌滑坡解譯標志。主要依據(jù)高精度DEM 顯示的歷史崩塌、滑坡形態(tài),輔以遙感影像的災害特征,建立相應的解譯標志。( 2) 利用ArcGIS 平臺,遵循建立的解譯標志,進行研究區(qū)的歷史崩塌滑坡識別。( 3) 現(xiàn)場驗證與成果完善,采用無人機和傳統(tǒng)工程地質(zhì)調(diào)查手段,進行現(xiàn)場驗證,驗證數(shù)為26 處,準確數(shù)為24 處,遙感解譯準確率可達92.3%。( 4) 對解譯的歷史崩塌滑坡進行細化圈定,基于10 m 精度DEM,利用統(tǒng)計學方法,分析崩塌滑坡與各因素的相關性,從而形成精細地形數(shù)據(jù)的崩塌滑坡識別方法。

    對研究區(qū)滑坡、崩塌災害進行遙感解譯,共計解譯崩塌41 處,滑坡92 處,分布情況如圖1 所示。

    2.2 崩塌、滑坡規(guī)模識別

    基于ArcGIS 平臺,圈定崩塌滑坡的范圍之后,可以在屬性表中自動計算各災害點的面積,為獲取各災害點的規(guī)模,還需獲取崩塌滑坡的厚度。利用高精度DEM 數(shù)據(jù),對這些崩塌滑坡截取剖面。滑坡深度的估算主要基于現(xiàn)場調(diào)查,依據(jù)滑坡前緣高差和后緣陡坎的高差,通過在GIS 中截取滑坡剖面并獲取其陡坎高差實現(xiàn),而崩塌堆積物主要呈錐形,其平均厚度取最大厚度值的二分之一。

    圖1 研究區(qū)崩塌滑坡識別結(jié)果Fig. 1 Identification results of collapse and landslides in the study area

    研究區(qū)內(nèi)部分滑坡由于風化作用,后緣陡坎已經(jīng)漸趨平緩,后緣不存在明顯陡坎,而能夠在坡腳形成堆積體的大型崩塌較少,因此選取其中50 處具有明顯陡坎的滑坡以及6 處具有錐形堆積體的大型崩塌,共計56 處崩塌滑坡災害點進行統(tǒng)計。分析發(fā)現(xiàn)災害面積主要集中在105~106im2,占災害總數(shù)的66.04%。將面積與厚度進行乘法計算后,發(fā)現(xiàn)災害體積主要集中在105~107im3,占災害總數(shù)的84.91%。

    對這56 處崩塌滑坡的面積與規(guī)模進行分析,發(fā)現(xiàn)地質(zhì)災害體積與面積大致呈現(xiàn)指數(shù)關系,Larsen et al. ( 2010) 認為,滑坡( 包括巖質(zhì)滑坡和土質(zhì)滑坡) 的體積和面積滿足V = αAλ,對統(tǒng)計結(jié)果進行擬合( 圖2) ,根據(jù)擬合結(jié)果,α = 5.6367,λ = 1.154,R2= 0.7558,即V = 5.636 7A1.154,擬合效果較好,反映了數(shù)據(jù)的準確性和規(guī)律性,并可以通過此公式估算區(qū)內(nèi)崩塌滑坡體積。

    圖2 研究區(qū)地質(zhì)災害體積與面積關系擬合Fig. 2 Fitting curve of geological hazard volume and area

    3 研究區(qū)崩塌滑坡發(fā)育主控因素分析

    3.1 地層巖性

    崩塌主要集中發(fā)育在花崗巖、礫巖等硬脆巖體之中( Wang et al.,2012,2014) 研究區(qū)的41 處崩塌均發(fā)生在硬脆的礫巖和花崗巖區(qū)域,花崗巖位于高海拔位置,坡度通常較大,受到內(nèi)外動力地質(zhì)作用,極易發(fā)生崩塌,礫巖巖體由于坡度較陡也常伴有崩塌發(fā)生,多為中型崩塌。千枚巖具有千枚構造,在構造作用下,巖體破碎,且坡度通常較緩,因而92 處滑坡均在千枚巖中,以中型滑坡為主,邊界通常為硬巖和軟巖的接觸界面( 圖3) 。

    圖3 研究區(qū)特征剖面( 扎西林村)Fig. 3 Typical section of the study area( Zhaxilin Village)

    3.2 地質(zhì)構造

    區(qū)內(nèi)伴生雅魯藏布江縫合帶的斷裂,走向呈現(xiàn)近東西向,其傾向大致為180°~220°,傾角在65°左右,斷裂帶對區(qū)內(nèi)地質(zhì)災害發(fā)育具有控制作用。

    3.2.1 斷層

    結(jié)合現(xiàn)場對斷層位置的進一步確認和修正,對92 處滑坡與斷裂的相對位置進行統(tǒng)計分析發(fā)現(xiàn),以斷層為前邊界的牽引式滑坡數(shù)量最多,有46 處,占滑坡總數(shù)的50%; 以斷層為后邊界的推動式滑坡最少,為13 處,約占滑坡總數(shù)的14.1%; 滑坡體前部和后部均有斷層穿過的混合式滑坡為33 處,占滑坡總數(shù)的35.9%。對滑坡的滑坡方向與斷裂的走向進行統(tǒng)計分析,發(fā)現(xiàn)有78 處滑坡的滑動方向垂直于斷裂走向,約占滑坡總數(shù)的84.8%,如扎西林村滑坡、拉崗村滑坡等。在斷裂帶的影響下,斷裂附近崩塌分布較為密集,且普遍規(guī)模較大。

    3.2.2 巖體結(jié)構

    研究區(qū)內(nèi)主要發(fā)育滑移式崩塌、傾倒式崩塌和拉裂式崩塌。以加查鎮(zhèn)某崩塌點為例,分析該類崩塌的破壞模式和控制特征。崩塌點的破壞模式分析主要依靠無人機航攝技術,利用小型多旋翼無人機獲取巖體的結(jié)構面信息( 劉海洋等,2017) 。工作區(qū)內(nèi)高差大,為保證地面分辨率能達到設計要求,采取旁向65%重疊度,航向75%重疊度,航高大致保持在300 m。航空影像數(shù)據(jù)主要包括正射影像數(shù)據(jù)和三維模型影像數(shù)據(jù)。

    三維模型的生產(chǎn)主要通過影像密集匹配,三角網(wǎng)構建、自動賦予紋理等步驟,最終得到三維模型。應用AgisoftPhotoscan 建模軟件對三維模型進行制作加工,處理成果影像,完成三維實景模型制作( 圖4a) 。選取巖體結(jié)構明顯的部分坡體進行結(jié)構面提取分析,將點云數(shù)據(jù)以文本文件的形式,從Photoscan 中導入Coltop3D( 圖4b) ,利用Coltop3D軟件提取巖體結(jié)構產(chǎn)狀信息。

    圖4 加查鎮(zhèn)某崩塌三維點云模型Fig. 4 A 3D point cloud model of a collapse in Jiacha Town

    此崩塌點以礫巖為主,礫石成分為花崗質(zhì),主要發(fā)育3 組結(jié)構面,其中J1為礫巖層面。利用Coltop軟件進行產(chǎn)狀統(tǒng)計,共統(tǒng)計57 條結(jié)構面。J1的產(chǎn)狀統(tǒng)計中,最大值為199°∠46°,最小值為180°∠35°,均值為190°∠40°; J2的產(chǎn)狀統(tǒng)計中,最大值為231°∠80°,最小值為208°∠59°,均值為221°∠74°; J3的產(chǎn)狀統(tǒng)計中,最大值為289°∠90°,最小值為263°∠70°,均值為271°∠87°; 取均值作為這3 組結(jié)構面產(chǎn)狀,即:J1=190°∠40°,J2=221°∠74°,J3=271°∠87°,坡面產(chǎn)狀為203°∠50°。

    統(tǒng)計3 組結(jié)構面間距,其中J1結(jié)構面間距在3.2~20.1 m 之間,在5 ~6 m 間距范圍內(nèi)的結(jié)構面較多,此外,間距范圍在3 ~9 m 各間距范圍內(nèi)的結(jié)構面間距分布較均勻,均值為5.5 m,間距較小; J2結(jié)構面間距在1.3 ~10.7 m 之間,主要集中在4 ~5 m 之間,均值為8.5 m,間距較大; J3結(jié)構面間距在1.7~72.1 m 之間,主要集中在2~3 m 范圍之間,均值為18.0 m,間距大。就斜坡整體而言,各組結(jié)構面間距較大,切割較為強烈,巖體較破碎,容易發(fā)生大塊石的崩落。

    圖5 主要結(jié)構面的穩(wěn)定性分析Fig. 5 Stability analysis of main structural planes

    利用這57 條結(jié)構面的產(chǎn)狀和坡面產(chǎn)狀繪制赤平極射投影圖( 圖5) 。由于上半球投影較為常用,因而采用此投影方法。其中,綠色線段DO,EO,F(xiàn)O的方向分別為J1、J2、J3傾向,坡面傾向為MO,與J1產(chǎn)狀相近。J1、J2、J3分別相交于A、B、C 3 點,連接3點與圓心O 點,形成的紅色線段即為交線AO,BO和CO。根據(jù)李智毅等( 1994) ,J2與J3的交線BO傾向與坡面相近,但傾角大于坡角,因而不易沿J2、J3交線滑動; J1、J2交線AO 和J1、J3的交線CO 傾角均小于坡角,AO 位于J1、J2傾向線DO 和EO 的夾角之外,因而楔形體沿著J1發(fā)生滑動,J2僅起到切割作用; CO 在J1、J3傾向線DO 和FO 的夾角之外,因而楔形體沿J1發(fā)生滑動,J3同樣僅起到切割作用。綜合上述分析,結(jié)構面切割形成的楔形體主要沿J1發(fā)生滑移式崩塌。

    3.3 地形地貌

    分析崩塌滑坡識別結(jié)果可以發(fā)現(xiàn),滑坡主要發(fā)育在雅魯藏布江南岸以及北岸坡體的中下部,而崩塌主要發(fā)生在北岸坡體的中上部??刂票浪?、滑坡發(fā)生的地形地貌因素主要為坡度、坡向和高程,為準確確定崩塌、滑坡受坡度、坡向和高程的控制情況,對于崩塌物源區(qū)的2390 個柵格和滑坡面內(nèi)的311 792 個柵格進行統(tǒng)計分析。利用ArcGIS 軟件的轉(zhuǎn)換工具和提取工具,將這些柵格的坡度、坡向以及高程信息提取至點文件。將坡度數(shù)據(jù)按5°間隔,劃分為18 個等級; 將坡向按照20°間隔,劃分為18 個等級; 將高程按照100 m 間隔,劃分為24 個等級。

    按照這個分類對于全區(qū)內(nèi)坡度、坡向和高程數(shù)據(jù)進行重分類,并統(tǒng)計各數(shù)據(jù)范圍內(nèi)的面積,以崩塌和滑坡面范圍內(nèi)各坡度、坡向、高程內(nèi)的柵格面積與區(qū)內(nèi)各坡度、坡向、高程范圍內(nèi)的柵格面積作比,作為崩塌和滑坡的面密度,由于每個柵格面積相同,即可轉(zhuǎn)化為柵格個數(shù)的比值。

    分析發(fā)現(xiàn),在0°~80°的坡度范圍內(nèi),崩塌面密度隨坡度上升而大致呈現(xiàn)增大趨勢,而滑坡面密度在20°~25°坡度范圍內(nèi)之間呈現(xiàn)一個峰值,而后隨坡度上升呈現(xiàn)下降趨勢( 圖6) 。不同坡向的崩塌面密度變化趨勢與崩塌數(shù)量一致,主要分布在140°~260°的坡向范圍內(nèi); 滑坡面密度與坡向關系較小,不同坡度范圍內(nèi)滑坡分布較為離散( 圖7) 。崩塌面密度隨高程上升而大致呈現(xiàn)增大趨勢,主要分布在3800~4900 m 高程范圍內(nèi); 滑坡在3200~4000 m 高程范圍內(nèi)的分布情況較為分散,在3800 m 以下和4000 m 以上坡度范圍內(nèi)滑坡面密度較低( 圖8) 。

    圖6 災害與坡度關系圖Fig. 6 Relationship between hazards and slope

    圖7 災害與坡向關系圖Fig. 7 Relationship between hazards and slope direction

    圖8 災害與高程關系圖Fig. 8 Relationship between hazards and elevation

    為研究各因子對崩塌滑坡發(fā)生的影響程度,應用曲線下面積法( AUC,Area Under the Curve) 。以各因子不同等級柵格累計比值為橫軸,災害( 崩塌或滑坡) 柵格累計比值為縱軸,面密度從高到低排列,用每組的面密度值高低表征發(fā)生崩塌或滑坡可能性的大小。利用Origin2019 專業(yè)繪圖軟件,對兩組數(shù)據(jù)繪制折線圖,即為受試者工作特征曲線( ROC曲線,Receiver-operating Characteristic Curve) 。

    對于各ROC 曲線求其AUC 值,分析發(fā)現(xiàn),針對坡度而言,崩塌的AUC 值為0.914,略大于滑坡的AUC 值( 0.835) ,表明坡度對于崩塌的控制作用明顯大于滑坡; 就坡向和高程而言,崩塌的AUC( 0.787 和0.848) 值均遠大于滑坡AUC 值( 0.528和0.682) ,表明坡向和高程對于崩塌具有一定的控制作用,而對滑坡的影響更小。因此,可以認為研究區(qū)內(nèi)的崩塌由坡度、坡向和高程共同控制,其中坡度為主控因素; 而滑坡與坡度的相關性較大,高程和坡向?qū)碌陌l(fā)育具有一定的影響,但影響較小。

    4 結(jié) 論

    (1) 利用遙感影像數(shù)據(jù)與10 m 精度無人機航攝數(shù)據(jù),對崩塌滑坡進行識別,共計解譯崩塌41 處,滑坡92 處,在構造作用和巖性控制下,呈現(xiàn)條帶狀分布。崩塌主要集中在斷層帶附近的礫巖、花崗巖巖體,以滑移式崩塌為主,而滑坡主要發(fā)育在千枚巖區(qū)域,以牽引型為主,84.8%的滑坡滑動方向與斷裂走向大致垂直。

    (2) 基于10 m 精度的DEM 數(shù)據(jù),通過截取崩塌滑坡剖面獲取其厚度( 或深度) ,發(fā)現(xiàn)區(qū)內(nèi)崩塌滑坡的體積和面積滿足公式V=5.636 7A1.154,擬合效果較好,可為崩塌滑坡災害的規(guī)律識別提供參考。

    (3) 從地層巖性、地質(zhì)構造、地形地貌3 個方面分析崩塌的主控因素,表明崩塌的主控因素較多,地層巖性決定了研究區(qū)崩塌和滑坡的發(fā)育類型,地形中的坡度、坡向和高程因素均對崩塌的形成重要影響作用,其中,坡度與崩塌發(fā)生的相關性最大,滑坡的發(fā)育受斷層的控制作用明顯,地形地貌中的坡度因子對滑坡影響作用較大,而坡向和高程對于滑坡的貢獻較小。

    猜你喜歡
    面密度坡向柵格
    基于預測PI的電極面密度控制系統(tǒng)
    基于預測PI的電極面密度控制系統(tǒng)
    化工機械(2025年1期)2025-01-31 00:00:00
    一種適用于復雜加載狀態(tài)下微噴射物質(zhì)面密度測量的Asay 膜方法
    爆炸與沖擊(2024年3期)2024-11-01 00:00:00
    基于鄰域柵格篩選的點云邊緣點提取方法*
    DEM地表坡向變率的向量幾何計算法
    測繪學報(2019年11期)2019-11-20 01:31:42
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    車身空腔用膨脹膠塊聲學性能影響因素分析
    粘接(2014年9期)2014-04-29 20:57:28
    基于CVT排布的非周期柵格密度加權陣設計
    雷達學報(2014年4期)2014-04-23 07:43:13
    動態(tài)柵格劃分的光線追蹤場景繪制
    国内少妇人妻偷人精品xxx网站| 最近的中文字幕免费完整| 一级av片app| 一级毛片aaaaaa免费看小| 伦理电影大哥的女人| 亚洲真实伦在线观看| 777米奇影视久久| 成人无遮挡网站| 中文字幕亚洲精品专区| 免费看光身美女| 国产亚洲5aaaaa淫片| 夫妻午夜视频| 秋霞在线观看毛片| 久久毛片免费看一区二区三区| 男女免费视频国产| 色哟哟·www| 国内精品宾馆在线| 又爽又黄a免费视频| 国产黄片美女视频| 亚洲欧洲国产日韩| 一级毛片我不卡| 直男gayav资源| 内射极品少妇av片p| 2022亚洲国产成人精品| 视频中文字幕在线观看| 国产精品蜜桃在线观看| 亚洲国产日韩一区二区| 免费播放大片免费观看视频在线观看| 91在线精品国自产拍蜜月| 性色av一级| kizo精华| 免费观看a级毛片全部| 偷拍熟女少妇极品色| 午夜激情福利司机影院| 91狼人影院| 欧美亚洲 丝袜 人妻 在线| 精品一区在线观看国产| 欧美3d第一页| 看十八女毛片水多多多| 最近最新中文字幕免费大全7| videossex国产| 老熟女久久久| 午夜福利在线观看免费完整高清在| 久久久久久久久大av| 亚洲色图综合在线观看| 精品99又大又爽又粗少妇毛片| 国产在线视频一区二区| 视频区图区小说| 国精品久久久久久国模美| 男人狂女人下面高潮的视频| 一个人看的www免费观看视频| 成人国产麻豆网| 高清不卡的av网站| av在线app专区| 老司机影院毛片| 日韩不卡一区二区三区视频在线| 免费观看的影片在线观看| 中文在线观看免费www的网站| 久久国产精品大桥未久av | 偷拍熟女少妇极品色| 99热这里只有是精品在线观看| 日韩成人伦理影院| 久久精品久久久久久噜噜老黄| 久久97久久精品| 极品教师在线视频| 国产高清不卡午夜福利| 亚洲av福利一区| 一本一本综合久久| 精品一品国产午夜福利视频| av在线观看视频网站免费| 舔av片在线| 妹子高潮喷水视频| 欧美最新免费一区二区三区| av免费在线看不卡| 下体分泌物呈黄色| 亚洲在久久综合| 亚洲熟女精品中文字幕| 久久午夜福利片| av黄色大香蕉| av网站免费在线观看视频| 亚洲内射少妇av| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添av毛片| 久久久久精品久久久久真实原创| 午夜老司机福利剧场| 国产成人a区在线观看| 日本黄大片高清| 2018国产大陆天天弄谢| 色视频www国产| 国产探花极品一区二区| 亚洲av.av天堂| 美女cb高潮喷水在线观看| 国产一区有黄有色的免费视频| 欧美xxⅹ黑人| 超碰av人人做人人爽久久| 婷婷色av中文字幕| 在线天堂最新版资源| 国产成人91sexporn| 中国国产av一级| 国内精品宾馆在线| 精品亚洲成a人片在线观看 | 性色av一级| 精品少妇黑人巨大在线播放| 3wmmmm亚洲av在线观看| 国产精品久久久久成人av| 亚洲丝袜综合中文字幕| 五月伊人婷婷丁香| av不卡在线播放| 联通29元200g的流量卡| 亚洲精品亚洲一区二区| 国产精品久久久久久av不卡| 成人综合一区亚洲| 欧美日韩视频精品一区| 能在线免费看毛片的网站| tube8黄色片| 大话2 男鬼变身卡| 国产精品精品国产色婷婷| 色婷婷av一区二区三区视频| 麻豆成人午夜福利视频| 国产精品99久久久久久久久| 亚洲国产精品专区欧美| a 毛片基地| 日韩成人伦理影院| 中国三级夫妇交换| 国语对白做爰xxxⅹ性视频网站| 热re99久久精品国产66热6| 纯流量卡能插随身wifi吗| 亚洲成人一二三区av| 欧美3d第一页| 极品教师在线视频| 在线观看免费视频网站a站| 一个人看视频在线观看www免费| 国产熟女欧美一区二区| 国产 一区精品| 一级毛片电影观看| 高清欧美精品videossex| 国产v大片淫在线免费观看| 80岁老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 在线观看免费日韩欧美大片 | 十八禁网站网址无遮挡 | av国产免费在线观看| 在线观看国产h片| av女优亚洲男人天堂| 国产欧美日韩精品一区二区| 在线观看免费日韩欧美大片 | 在线观看美女被高潮喷水网站| 中文字幕免费在线视频6| 欧美区成人在线视频| 久久精品久久久久久噜噜老黄| 日韩三级伦理在线观看| 亚洲无线观看免费| 国产在线免费精品| 夫妻性生交免费视频一级片| 男人舔奶头视频| 精品午夜福利在线看| 最近中文字幕高清免费大全6| 五月天丁香电影| 成年女人在线观看亚洲视频| 精品久久久久久久久亚洲| 一级毛片我不卡| 久久久精品免费免费高清| 久久99蜜桃精品久久| 卡戴珊不雅视频在线播放| av在线老鸭窝| 亚洲天堂av无毛| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 国产69精品久久久久777片| 亚洲av二区三区四区| 大片电影免费在线观看免费| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 久久影院123| 国产在线一区二区三区精| 成人毛片60女人毛片免费| 欧美日韩视频高清一区二区三区二| 99久久人妻综合| 国产高潮美女av| 久久青草综合色| 国产 一区 欧美 日韩| 久久热精品热| av国产久精品久网站免费入址| 免费人成在线观看视频色| 免费观看a级毛片全部| 午夜福利影视在线免费观看| 亚洲欧美成人精品一区二区| 在线观看三级黄色| 欧美极品一区二区三区四区| 午夜日本视频在线| 91精品伊人久久大香线蕉| 亚洲三级黄色毛片| 国产精品秋霞免费鲁丝片| 狂野欧美激情性bbbbbb| 网址你懂的国产日韩在线| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 中文欧美无线码| 97热精品久久久久久| 国产精品国产av在线观看| 少妇丰满av| 丝袜脚勾引网站| 久久鲁丝午夜福利片| 亚州av有码| 国产亚洲一区二区精品| 午夜福利网站1000一区二区三区| 午夜福利视频精品| 欧美+日韩+精品| 七月丁香在线播放| 特大巨黑吊av在线直播| 在线看a的网站| 久久久久久久精品精品| 成人午夜精彩视频在线观看| 国产精品一及| 久久久久久久大尺度免费视频| 日本-黄色视频高清免费观看| 美女xxoo啪啪120秒动态图| 又爽又黄a免费视频| 国产亚洲5aaaaa淫片| 国产在视频线精品| 黄色一级大片看看| 久久影院123| 丝袜喷水一区| 日本av免费视频播放| 国产亚洲精品久久久com| 尾随美女入室| 一级a做视频免费观看| 国产91av在线免费观看| 汤姆久久久久久久影院中文字幕| 久久久色成人| 99精国产麻豆久久婷婷| 啦啦啦视频在线资源免费观看| 男人添女人高潮全过程视频| www.色视频.com| 秋霞在线观看毛片| av网站免费在线观看视频| 好男人视频免费观看在线| 人人妻人人看人人澡| 色吧在线观看| 欧美高清性xxxxhd video| 亚洲av国产av综合av卡| 亚洲国产色片| 亚洲国产高清在线一区二区三| 黄色怎么调成土黄色| 中文字幕久久专区| 99久久中文字幕三级久久日本| 国产无遮挡羞羞视频在线观看| 国产免费福利视频在线观看| av国产久精品久网站免费入址| 日韩成人伦理影院| 丰满乱子伦码专区| 亚洲三级黄色毛片| 爱豆传媒免费全集在线观看| 国产久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区 | 国产一区有黄有色的免费视频| 婷婷色综合大香蕉| 汤姆久久久久久久影院中文字幕| 好男人视频免费观看在线| 国产成人午夜福利电影在线观看| 亚州av有码| 精品99又大又爽又粗少妇毛片| freevideosex欧美| a 毛片基地| 欧美区成人在线视频| 亚洲精品久久午夜乱码| 国产伦理片在线播放av一区| 久久久国产一区二区| 国产精品久久久久久久久免| 日日啪夜夜撸| 欧美精品亚洲一区二区| 最近中文字幕2019免费版| 午夜免费鲁丝| 中国美白少妇内射xxxbb| 777米奇影视久久| 久久精品国产a三级三级三级| 精品人妻视频免费看| 日韩 亚洲 欧美在线| 免费观看av网站的网址| 丰满少妇做爰视频| 一个人看的www免费观看视频| 青青草视频在线视频观看| 国产淫语在线视频| xxx大片免费视频| 99久久人妻综合| 亚洲精品日韩av片在线观看| 三级国产精品片| 男的添女的下面高潮视频| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 国产黄片视频在线免费观看| 尾随美女入室| 欧美性感艳星| 人妻一区二区av| 男女边吃奶边做爰视频| 男女无遮挡免费网站观看| 久久久久久伊人网av| 91精品国产九色| 男女边摸边吃奶| 中文乱码字字幕精品一区二区三区| 欧美日韩综合久久久久久| 亚洲,一卡二卡三卡| 18+在线观看网站| av又黄又爽大尺度在线免费看| 国产91av在线免费观看| 久久精品久久久久久噜噜老黄| 国产无遮挡羞羞视频在线观看| 天天躁夜夜躁狠狠久久av| 亚洲欧美清纯卡通| av.在线天堂| 欧美日韩国产mv在线观看视频 | 亚洲四区av| www.av在线官网国产| 99九九线精品视频在线观看视频| 尾随美女入室| 国产综合精华液| 日韩在线高清观看一区二区三区| 激情 狠狠 欧美| a级一级毛片免费在线观看| 亚洲一级一片aⅴ在线观看| 国产在线一区二区三区精| 99热国产这里只有精品6| 性色av一级| 亚洲精品日韩在线中文字幕| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 在现免费观看毛片| 制服丝袜香蕉在线| 久久久久网色| 交换朋友夫妻互换小说| 麻豆成人午夜福利视频| 久久久久久九九精品二区国产| 日韩亚洲欧美综合| 亚洲av欧美aⅴ国产| 99久久精品国产国产毛片| 九九在线视频观看精品| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 一区二区三区四区激情视频| 黄色配什么色好看| 亚洲欧美日韩卡通动漫| 三级经典国产精品| 久久午夜福利片| 黄色日韩在线| 午夜老司机福利剧场| 老司机影院毛片| 国产成人a∨麻豆精品| 成人高潮视频无遮挡免费网站| 深爱激情五月婷婷| 亚洲美女搞黄在线观看| 欧美精品一区二区大全| 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 全区人妻精品视频| 一级a做视频免费观看| 男人狂女人下面高潮的视频| 九九在线视频观看精品| 搡老乐熟女国产| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 国产免费一区二区三区四区乱码| 99久久精品一区二区三区| 国内揄拍国产精品人妻在线| 国产片特级美女逼逼视频| 黄色配什么色好看| 麻豆国产97在线/欧美| 蜜桃久久精品国产亚洲av| 男女啪啪激烈高潮av片| av卡一久久| 中国美白少妇内射xxxbb| 国产一区亚洲一区在线观看| 国产欧美日韩一区二区三区在线 | 天堂俺去俺来也www色官网| 亚洲国产高清在线一区二区三| 久久精品久久精品一区二区三区| 欧美97在线视频| 亚洲伊人久久精品综合| 涩涩av久久男人的天堂| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 麻豆国产97在线/欧美| 亚洲av中文字字幕乱码综合| 国产视频首页在线观看| 精品少妇黑人巨大在线播放| 毛片一级片免费看久久久久| 国产高清不卡午夜福利| 男人和女人高潮做爰伦理| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 91午夜精品亚洲一区二区三区| 国产中年淑女户外野战色| 亚洲精品视频女| 99热全是精品| 亚洲精品日韩av片在线观看| av线在线观看网站| 亚洲精品乱久久久久久| 国产伦精品一区二区三区四那| 男女边吃奶边做爰视频| 高清视频免费观看一区二区| av在线播放精品| 五月开心婷婷网| 国产免费一区二区三区四区乱码| 亚洲欧美一区二区三区黑人 | 久热这里只有精品99| 伦理电影免费视频| 人人妻人人爽人人添夜夜欢视频 | 国产老妇伦熟女老妇高清| 我的老师免费观看完整版| 熟女人妻精品中文字幕| 欧美精品一区二区免费开放| 亚洲不卡免费看| 七月丁香在线播放| h视频一区二区三区| 黑人猛操日本美女一级片| 中国美白少妇内射xxxbb| 亚洲精品久久午夜乱码| 老司机影院毛片| 毛片女人毛片| 2022亚洲国产成人精品| 午夜精品国产一区二区电影| 国产精品国产三级国产av玫瑰| 在线观看一区二区三区激情| 男女边吃奶边做爰视频| 亚洲在久久综合| 天天躁日日操中文字幕| 免费观看在线日韩| 亚洲人与动物交配视频| 亚洲成人中文字幕在线播放| 中文字幕人妻熟人妻熟丝袜美| 爱豆传媒免费全集在线观看| 午夜免费男女啪啪视频观看| 亚洲,一卡二卡三卡| 99精国产麻豆久久婷婷| 国产 精品1| 亚洲精品久久久久久婷婷小说| 免费看av在线观看网站| 国产成人91sexporn| 日韩视频在线欧美| 亚洲av电影在线观看一区二区三区| 日韩欧美一区视频在线观看 | 日韩视频在线欧美| 日韩伦理黄色片| 久久99热6这里只有精品| 天堂中文最新版在线下载| 女人十人毛片免费观看3o分钟| 午夜免费男女啪啪视频观看| 少妇人妻 视频| 建设人人有责人人尽责人人享有的 | 午夜福利在线在线| 中文字幕制服av| 成人国产av品久久久| 国产精品久久久久久久久免| 高清日韩中文字幕在线| 日韩欧美一区视频在线观看 | 国产精品蜜桃在线观看| 丰满人妻一区二区三区视频av| 日韩电影二区| 观看美女的网站| 一个人看的www免费观看视频| 一本一本综合久久| 国产精品成人在线| 尤物成人国产欧美一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲在久久综合| 日韩欧美精品免费久久| 日韩av在线免费看完整版不卡| 嘟嘟电影网在线观看| 色婷婷av一区二区三区视频| a级一级毛片免费在线观看| 日本av手机在线免费观看| 波野结衣二区三区在线| 成年人午夜在线观看视频| 亚洲精品,欧美精品| 午夜精品国产一区二区电影| 自拍欧美九色日韩亚洲蝌蚪91 | 日本av免费视频播放| 亚洲国产精品国产精品| 国产亚洲欧美精品永久| 大陆偷拍与自拍| 99热网站在线观看| 中文字幕免费在线视频6| 精华霜和精华液先用哪个| 色视频www国产| 国产精品国产三级专区第一集| 简卡轻食公司| 在线看a的网站| 久久久色成人| 国产成人精品久久久久久| 免费观看的影片在线观看| 你懂的网址亚洲精品在线观看| 美女内射精品一级片tv| 国产精品一区二区在线观看99| av在线播放精品| 日韩,欧美,国产一区二区三区| 久久久久久久精品精品| 肉色欧美久久久久久久蜜桃| 人妻系列 视频| 免费播放大片免费观看视频在线观看| 日韩电影二区| 亚洲,一卡二卡三卡| 亚洲欧美一区二区三区黑人 | 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 欧美精品人与动牲交sv欧美| 国产精品蜜桃在线观看| 久久人人爽人人爽人人片va| 国产在线视频一区二区| 五月天丁香电影| 精品人妻视频免费看| 99久久人妻综合| 一级毛片 在线播放| 国产成人精品一,二区| 亚洲av福利一区| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 丝袜喷水一区| 久久久成人免费电影| 免费黄网站久久成人精品| 亚洲av成人精品一二三区| 免费少妇av软件| 校园人妻丝袜中文字幕| 三级国产精品欧美在线观看| 一级av片app| 亚洲自偷自拍三级| 高清视频免费观看一区二区| 国产亚洲精品久久久com| 国产午夜精品久久久久久一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 激情 狠狠 欧美| 欧美日韩综合久久久久久| 少妇熟女欧美另类| 国产成人91sexporn| 亚洲精华国产精华液的使用体验| 日韩欧美精品免费久久| 亚洲精品,欧美精品| 老女人水多毛片| 久久久久人妻精品一区果冻| av福利片在线观看| 少妇 在线观看| 美女国产视频在线观看| 色5月婷婷丁香| 亚洲国产最新在线播放| 国产黄频视频在线观看| 在线观看人妻少妇| 观看美女的网站| 人妻夜夜爽99麻豆av| 中文字幕亚洲精品专区| 国国产精品蜜臀av免费| 性高湖久久久久久久久免费观看| 久久久久国产精品人妻一区二区| 亚洲国产精品专区欧美| 亚洲国产精品999| 婷婷色综合www| 亚洲熟女精品中文字幕| 成年女人在线观看亚洲视频| 毛片一级片免费看久久久久| 婷婷色麻豆天堂久久| 在线观看人妻少妇| 校园人妻丝袜中文字幕| 91久久精品电影网| 国产爽快片一区二区三区| 精品久久国产蜜桃| 国产高清不卡午夜福利| 久久精品久久久久久久性| 亚洲自偷自拍三级| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 国产精品秋霞免费鲁丝片| 国产精品不卡视频一区二区| 毛片女人毛片| 精品国产乱码久久久久久小说| av专区在线播放| 欧美成人一区二区免费高清观看| 免费人妻精品一区二区三区视频| 国产成人精品一,二区| 啦啦啦视频在线资源免费观看| 日韩av在线免费看完整版不卡| 精品亚洲成a人片在线观看 | 全区人妻精品视频| 国产精品久久久久久久久免| 亚洲国产欧美人成| 99精国产麻豆久久婷婷| 亚洲精品456在线播放app| 麻豆国产97在线/欧美| 亚洲成人av在线免费| 亚洲av综合色区一区| 国产国拍精品亚洲av在线观看| 亚洲av二区三区四区| 久久99热6这里只有精品| 国产亚洲91精品色在线| 人人妻人人澡人人爽人人夜夜| 免费观看的影片在线观看| 亚洲激情五月婷婷啪啪| 久热久热在线精品观看| 国产伦精品一区二区三区视频9| 国产色爽女视频免费观看| 国产男女内射视频| 少妇的逼水好多| 久久人人爽av亚洲精品天堂 | 久久国产亚洲av麻豆专区| 亚洲第一区二区三区不卡| 秋霞在线观看毛片| 视频区图区小说| 精品一区二区三卡| 国产综合精华液| 成人亚洲欧美一区二区av| 国产精品一区www在线观看| 中文欧美无线码| 国产大屁股一区二区在线视频| 国国产精品蜜臀av免费| 日本黄色日本黄色录像| 中文字幕精品免费在线观看视频 |