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

    基于地表變形數(shù)據(jù)的潛在滑坡識(shí)別研究

    2020-02-29 10:39:27李萌彭思佳白艷萍黃兆歡
    科技創(chuàng)新與應(yīng)用 2020年6期
    關(guān)鍵詞:試驗(yàn)區(qū)滑坡聚類

    李萌 彭思佳 白艷萍 黃兆歡

    摘? 要:潛在滑坡識(shí)別對(duì)地質(zhì)災(zāi)害高發(fā)區(qū)的滑坡研究有重要意義。文章以中巴經(jīng)濟(jì)走廊災(zāi)害多發(fā)區(qū)某一試驗(yàn)區(qū)為例,結(jié)合由2014-2019的Sentinel1-A雷達(dá)數(shù)據(jù)計(jì)算出的地表變形數(shù)據(jù),選取試驗(yàn)區(qū)地形、水文、地質(zhì)、聚類共8種特征變量基于隨機(jī)森林模型進(jìn)行滑坡識(shí)別,滑坡識(shí)別結(jié)果總體精度為88.62%,Kappa系數(shù)為0.78。結(jié)果表明,結(jié)合地表變形數(shù)據(jù)的隨機(jī)森林模型可以較好的識(shí)別潛在滑坡,聚類特征對(duì)于識(shí)別結(jié)果有重要貢獻(xiàn)。

    關(guān)鍵詞:隨機(jī)森林;中巴經(jīng)濟(jì)走廊;地質(zhì)災(zāi)害;特征變量

    中圖分類號(hào):P642.22 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-2945(2020)06-0009-05

    Abstract: The identification of potential landslides is of great significance for landslide research in high-risk areas of geological hazards. Taking a test area in a disaster-prone area of the China-Pakistan Economic Corridor as an example, combined with the surface deformation data set calculated from the Sentinel1-A radar data of 2014-2019, eight characteristic variables of terrain, hydrology, geology, and clustering in the test area were selected. Landslide identification was performed using a random forest model. The overall accuracy of the landslide identification result was 88.62%, and the Kappa coefficient was 0.78. The results show that the random forest model combined with the surface deformation data set can better identify the distribution of potential landslides, and the clustering features make an important contribution to the recognition results.

    Keywords: random forest; China-Pakistan Economic Corridor; geological disaster; characteristic variable

    引言

    滑坡是一種多發(fā)常見、分布廣、范圍大、破壞性嚴(yán)重的地質(zhì)災(zāi)害之一,嚴(yán)重威脅著人民生命財(cái)產(chǎn)安全[1]。傳統(tǒng)的識(shí)別方法周期長(zhǎng)、工作量大、危險(xiǎn)性高,而遙感技術(shù)具有時(shí)間快、成本低等優(yōu)勢(shì),將其應(yīng)用到滑坡的識(shí)別中對(duì)于滑坡發(fā)生的可能性預(yù)測(cè),助力防災(zāi)減災(zāi)救災(zāi)等綜合工作具有較大的價(jià)值。

    近年來,隨著傳感器的多元化,眾多學(xué)者使用不同數(shù)據(jù)源和方法進(jìn)行滑坡識(shí)別的研究。張群等借助Geoeye-1衛(wèi)星影像提取光譜、空間、紋理及鄰域特征,基于面向?qū)ο蠓椒ㄌ崛「拭C黑方臺(tái)黃土滑坡區(qū)域[2]。徐喬等使用高分辨率QuickBird 影像,提取光譜、紋理及形狀特征識(shí)別出疑似滑坡區(qū)[3]。余坤勇等利用ALOS 多光譜數(shù)據(jù)預(yù)測(cè)山體滑坡空間分布趨勢(shì),結(jié)果表明隨機(jī)森林模型預(yù)測(cè)精度高于邏輯回歸模型,達(dá)到90.8%[4]。Chaoying Z等使用一種典型的小基線子集(SBAS)干涉合成孔徑雷達(dá)(InSAR)技術(shù)分析烏東德水庫滑坡的時(shí)空變形特征[5]。Aggeliki K等基于Sentinel-1,Sentinel-2和Landsat-8多源遙感數(shù)據(jù)對(duì)馬其頓西北部某煤礦區(qū)進(jìn)行滑坡分布識(shí)別和繪制[6]?;伦R(shí)別的一個(gè)主流思路就是使用機(jī)器學(xué)習(xí)的方法,包括支持向量機(jī)、神經(jīng)網(wǎng)絡(luò)、隨機(jī)森林等[7-9],還有學(xué)者將其它方法應(yīng)用到滑坡識(shí)別,如主成分變換、物探法[10-11]等。

    本文借助地表變形數(shù)據(jù),即利用PS變形點(diǎn)進(jìn)行空間自相關(guān)分析,計(jì)算出平均速度的聚類特征,再結(jié)合地形、水文、地質(zhì)共8種特征變量,基于隨機(jī)森林算法,對(duì)中巴經(jīng)濟(jì)走廊某試驗(yàn)區(qū)進(jìn)行潛在滑坡識(shí)別和研究,并分析不同特征變量的重要性。

    1 研究區(qū)概況

    中巴經(jīng)濟(jì)走廊是中國(guó)和巴基斯坦兩國(guó)之間交通、能源、海洋等領(lǐng)域交流與合作的戰(zhàn)略通道,是我國(guó)“一帶一路”重點(diǎn)工程之一,走廊內(nèi)地質(zhì)環(huán)境極為復(fù)雜、自然災(zāi)害活躍,山體滑坡、落石、塌方、地震等地質(zhì)災(zāi)害時(shí)常發(fā)生,給人們的生命財(cái)產(chǎn)安全和社會(huì)經(jīng)濟(jì)發(fā)展帶來嚴(yán)重的威脅,因此對(duì)走廊內(nèi)的潛在滑坡識(shí)別很有必要和意義。

    2 數(shù)據(jù)源與處理

    2.1 雷達(dá)影像數(shù)據(jù)

    Sentinel-1衛(wèi)星是歐洲航天局哥白尼計(jì)劃(GMES)中的地球觀測(cè)衛(wèi)星,由兩顆衛(wèi)星組成,載有C波段合成孔徑雷達(dá),衛(wèi)星數(shù)據(jù)不易受天氣環(huán)境的影響,可提供不論晝夜和云雨霧等狀況下的連續(xù)圖像。從歐空局?jǐn)?shù)據(jù)中心(https:∥scihub.Copernicus.eu/dhus/#/home)獲取試驗(yàn)區(qū)2014年10月31日到2019年10月29日共122景Sentinel-1A IW SLC數(shù)據(jù),用于生成試驗(yàn)區(qū)的地表變形數(shù)據(jù)。獲取的Sentinel-1影像基本參數(shù)如表1所示。

    基于PS-InSAR技術(shù),輸入研究區(qū)Sentinel-1A影像,進(jìn)行裁剪配準(zhǔn)、相位估計(jì)、差分干涉、相位解纏等一系列處理步驟,計(jì)算得到試驗(yàn)區(qū)時(shí)間序列上的地表變形數(shù)據(jù)即PS變形點(diǎn),試驗(yàn)區(qū)位置及PS變形點(diǎn)如圖1 所示,不同顏色的PS變形點(diǎn)顯示了不同區(qū)間的平均形變速率,可以看出試驗(yàn)區(qū)內(nèi)沿著溝谷分布的區(qū)域較為穩(wěn)定,但是仍然存在下降速率較大的PS點(diǎn)。

    2.2 其它數(shù)據(jù)

    地形因子對(duì)滑坡的識(shí)別有很大的影響,獲取空間分辨率為30米的ASTER GDEM高程數(shù)據(jù),作為地形特征變量的提取。在ArcGIS中對(duì)DEM 數(shù)據(jù)進(jìn)行裁剪,提取坡度、坡向、曲率、水系。試驗(yàn)區(qū)的地質(zhì)構(gòu)造斷對(duì)滑坡也具有一定影響,從巴基斯坦地質(zhì)中心獲取的覆蓋研究區(qū)范圍斷層空間分布數(shù)據(jù),作為滑坡識(shí)別的又一特征變量。

    3 研究方法

    3.1 隨機(jī)森林模型

    隨機(jī)森林算法是由Breiman于2001年提出的一種非參數(shù)新型機(jī)器學(xué)習(xí)方法,即在變量(列)的使用和數(shù)據(jù)(行)的使用上進(jìn)行隨機(jī)化,生成很多分類樹,再匯總分類樹的結(jié)果,隨機(jī)森林在運(yùn)算量沒有顯著提高的前提下提高了預(yù)測(cè)精度,該方法分類準(zhǔn)確性高、抗噪能力較強(qiáng)、泛化性能較好,已經(jīng)廣泛應(yīng)用于分類和識(shí)別的相關(guān)研究中。

    3.2 樣本數(shù)據(jù)選取

    樣本數(shù)據(jù)來自于Google earth高清影像的目視解譯。結(jié)合試驗(yàn)區(qū)部分實(shí)地考察照片以及由Sentinel-1A計(jì)算出的PS變形點(diǎn)的平均速度分布范圍勾畫滑坡區(qū)域,生成滑坡區(qū)域內(nèi)共682個(gè)的滑坡點(diǎn),再利用ArcGIS 10. 2軟件在已知勾畫的滑坡區(qū)之外隨機(jī)生成相等數(shù)目的穩(wěn)定點(diǎn)。將獲取的682個(gè)滑坡點(diǎn)和682個(gè)穩(wěn)定點(diǎn),按照比例7:3隨機(jī)分配分為訓(xùn)練樣本集和驗(yàn)證樣本集,分別用于隨機(jī)森林識(shí)別滑坡的模型建立和精度驗(yàn)證,圖2 是試驗(yàn)區(qū)位置及選取的樣本空間分布,圖3 是在Google earth高清影像輔助樣本解譯的PS變形點(diǎn)和實(shí)地考察照片。

    3.3 特征變量選取

    滑坡是自然界多種因子綜合作用的結(jié)果,各個(gè)因子之間往往有著相互關(guān)聯(lián)的因果關(guān)系。滑坡發(fā)生與否,首先受到自身地理環(huán)境的影響,包括地形地貌、土壤環(huán)境和植被信息等。一般情況下,在地勢(shì)平坦區(qū)不易發(fā)生滑坡災(zāi)害,在地形起伏大的區(qū)域會(huì)隨著地震、降雨等事件而發(fā)生滑坡,且滑坡發(fā)生的概率會(huì)隨著坡度和高程的增加而增加。同時(shí),構(gòu)造面又為降雨等水流進(jìn)入斜坡提供了通道,故各種節(jié)理、裂隙、層面、斷層發(fā)育的斜坡,特別是當(dāng)平行和垂直斜坡的陡傾角構(gòu)造面及順坡緩傾的構(gòu)造面發(fā)育時(shí),最易發(fā)生滑坡。另外,離水系越近的土壤,其含水性越大,而含水性越大往往有利于滑坡的發(fā)生,因而到水系距離也是滑坡發(fā)生的一個(gè)影響因子。

    空間自相關(guān)是指一些變量在同一個(gè)分布區(qū)內(nèi)的觀測(cè)數(shù)據(jù)之間潛在的相互依賴性,如果這個(gè)變量的空間分布存在規(guī)律性,這個(gè)變量是空間自相關(guān)的。ArcGIS中的空間自相關(guān)分析,可以計(jì)算某個(gè)變量在某個(gè)位置上的值受相鄰位置上該變量值的影響程度,并對(duì)該變量空間分布的獨(dú)立性和隨機(jī)性進(jìn)行檢驗(yàn),從而表現(xiàn)變量的聚類特征。本文基于試驗(yàn)區(qū)的地表變形數(shù)據(jù)即PS變形點(diǎn),以平均速率作為變量,基于 Getis和Moran's I 兩種模型進(jìn)行統(tǒng)計(jì)分析,分別得出兩種模型的Gi*指數(shù)和Morans I指數(shù),公式如下。

    公式(1)中,Wij(d)為在d距離內(nèi)的空間相鄰權(quán)重矩陣,同樣地,若i與j相鄰,該Wij(d)為1,不相鄰為0。公式(2)中:N為樣本數(shù),Xi變量X在i位置的值,Xj變量X在j位置的值,Wij表示空間權(quán)重。綜上考慮,最終選取試驗(yàn)區(qū)地形、水文、地質(zhì)、聚類共8種特征變量(表2),提取結(jié)果如圖4所示。

    4 結(jié)果和分析

    4.1 分類與精度評(píng)價(jià)

    輸入8個(gè)特征變量以及訓(xùn)練樣本集,在隨機(jī)森林模型中,人為調(diào)整樹的T個(gè)數(shù)和變量個(gè)數(shù)m,計(jì)算RFC模型,設(shè)置T從1-100逐一變化對(duì)融合特征進(jìn)行滑坡識(shí)別,并統(tǒng)計(jì)總體精度,技術(shù)路線如圖5所示。

    通過隨機(jī)森林算法的特征變量和訓(xùn)練樣本的輸入和建模,得到識(shí)別滑坡的結(jié)果,并計(jì)算了滑坡識(shí)別結(jié)果的混淆矩陣,結(jié)果顯示滑坡識(shí)別的總體精度為88.62%,Kappa精度為0.78,用戶精度和生產(chǎn)精度分別為89.18%和87.87%(表3),總體來看識(shí)別效果良好。

    4.2 重要性評(píng)估

    隨機(jī)森林模型可以對(duì)輸入的特征變量進(jìn)行重要性排序,如圖6所示,可以看出不同特征變量在此次滑坡識(shí)別實(shí)驗(yàn)中的對(duì)識(shí)別結(jié)果的貢獻(xiàn)。

    由圖6可以看出對(duì)滑坡識(shí)別貢獻(xiàn)最大的是Gi*指數(shù),即熱(冷)點(diǎn)分析中的 Getis模型的聚類特征。地形特征中高程的貢獻(xiàn)最大,而坡向和曲率貢獻(xiàn)最小,坡度和坡向的評(píng)分接近,遠(yuǎn)低于高程,由此可以推測(cè)在海拔越高的地方發(fā)生滑坡的可能性越大。斷層因子的評(píng)分較高,在構(gòu)造方面,例如斷層發(fā)育的斜坡,當(dāng)平行和垂直斜坡的陡傾角構(gòu)造面及順坡緩傾的構(gòu)造面發(fā)育時(shí),最易發(fā)生滑坡,斷層評(píng)分也已說明研究區(qū)內(nèi)的滑坡易發(fā)區(qū)受斷層的影響因素較大。

    通過熱(冷)點(diǎn)分析計(jì)算的Gi*指數(shù)貢獻(xiàn)最大,主要是因?yàn)镻S變形點(diǎn)的平均速度的聚類特征,平均速度為負(fù),表示該點(diǎn)地面下沉,但是一個(gè)或者兩個(gè)的下沉點(diǎn)往往不一定代表滑坡點(diǎn),但是當(dāng)下沉點(diǎn)成一簇一簇聚集起來分布,則該區(qū)域最有可能發(fā)生滑坡,如特征變量的重要性排序所示,Gi*指數(shù)的貢獻(xiàn)遠(yuǎn)大于其它7個(gè)特征變量。

    4.3 數(shù)據(jù)對(duì)比驗(yàn)證

    通過分類后處理步驟設(shè)置閾值,消除分類結(jié)果中的微小細(xì)碎斑塊,得到圖7所示黑色不規(guī)則多邊形的滑坡區(qū)域,統(tǒng)計(jì)樣本點(diǎn)中分為滑坡點(diǎn)的穩(wěn)定點(diǎn)和分為穩(wěn)定點(diǎn)的滑坡點(diǎn),可以看出大多滑坡點(diǎn)都落在黑色不規(guī)則多邊形內(nèi),即滑坡點(diǎn)落在滑坡區(qū),而只有少數(shù)穩(wěn)定點(diǎn)誤分到滑坡點(diǎn)。出現(xiàn)誤分的原因可能有以下原因:樣本點(diǎn)本身存在的誤差;隨機(jī)森林模型的算法誤差;滑坡的發(fā)生機(jī)理和誘發(fā)因子復(fù)雜,特征選擇不夠全面。以上都有可能造成分類結(jié)果精度的誤差。

    通過將PS變形點(diǎn)覆蓋到滑坡識(shí)別結(jié)果中進(jìn)行對(duì)比(圖8),可以看出在空間趨勢(shì)上,PS變形點(diǎn)的下沉速度較大的點(diǎn)大多落在黑色不規(guī)則多邊形內(nèi),即下降位移較大的點(diǎn)落在識(shí)別的滑坡區(qū)里,通過PS變形點(diǎn)和識(shí)別結(jié)果的對(duì)比,可以看出地表變形數(shù)據(jù)和滑坡識(shí)別區(qū)域的對(duì)應(yīng)效果較好。

    5 結(jié)束語

    本文結(jié)合地表變形數(shù)據(jù)中的PS變形點(diǎn),通過空間自相關(guān)分析,使用 Getis和Moran's I兩種模型計(jì)算了平均速度的聚類特征,再結(jié)合地形、水文、地質(zhì)特征變量,對(duì)中巴經(jīng)濟(jì)走廊某試驗(yàn)區(qū)進(jìn)行潛在滑坡識(shí)別,識(shí)別結(jié)果總體精度為88.62%,Kappa系數(shù)為0.78。

    特征重要性排序中Gi*指數(shù)評(píng)分最高,表明下沉點(diǎn)成團(tuán)聚集的地方,最有可能發(fā)生滑坡。地形特征中高程的貢獻(xiàn)大于坡向和曲率,由此可以推測(cè)在海拔越高的地方發(fā)生滑坡的可能性越大。坡向和曲率和貢獻(xiàn)重要性較小,不是滑坡發(fā)生的主導(dǎo)因素。

    降水量尤其暴雨也是導(dǎo)致滑坡可能發(fā)生的重要因素,然而由于研究區(qū)降水?dāng)?shù)據(jù)缺乏,因此只考慮了地表變形數(shù)據(jù)、地形、水文、地質(zhì)等特征,而人類活動(dòng)等其它潛在因子也對(duì)滑坡產(chǎn)生影響,這可能對(duì)模型的預(yù)測(cè)結(jié)果造成一定的誤差。

    參考文獻(xiàn):

    [1]趙寶強(qiáng),韓守富,白艷萍,等.時(shí)序InSAR技術(shù)在大型滑坡監(jiān)測(cè)中的應(yīng)用[J].科技創(chuàng)新與應(yīng)用,2019(1):21-24.

    [2]張群,趙超英.基于面向?qū)ο蟮母叻诌b感數(shù)據(jù)甘肅黑方臺(tái)黃土滑坡半自動(dòng)識(shí)別[J].災(zāi)害學(xué),2017,32(3):210-215.

    [3]徐喬,孟凡利,余紹淮.基于面向?qū)ο箅S機(jī)森林分類模型的滑坡遙感解譯[J].中外公路,2019,39(3):30-33.

    [4]余坤勇,姚雄,邱祈榮,等.基于隨機(jī)森林模型的山體滑坡空間預(yù)測(cè)研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(10):338-345.

    [5]Chaoying Z,Ya K,Qin Z,et al.Landslide Identification and Monitoring along the Jinsha River Catchment (Wudongde Reservoir Area),China,Using the InSAR Method[J].Remote sensing,2018, 10, 993.

    [6]Aggeliki K,Konstantinos N.Landslide mapping using optical and radar data a case study from Aminteo Western Macedonia Greece[J].European Journal of Remote Sensing,2019,2279-7254.

    [7]李松,鄧寶昆,徐紅勤,等.地震型滑坡災(zāi)害遙感快速識(shí)別方法研究[J].遙感信息,2015,30(4):25-28.

    [8]李勛,楊環(huán),殷宗敏,等.基于DEM和遙感影像的區(qū)域黃土滑坡體識(shí)別方法研究[J].地理與地理信息科學(xué),2017,33(4):86-92.

    [9]姚雄,余坤勇,劉健,等.基于隨機(jī)森林模型的降水誘發(fā)山體滑坡空間預(yù)測(cè)技術(shù)[J].福建農(nóng)林大學(xué)學(xué)報(bào):自然科學(xué)版,2016,45(2):219-227.

    [10]李成龍,張景發(fā).基于主成分分析的遙感震害變化檢測(cè)方法與應(yīng)用[J].地震,2013,33(2):103-108.

    [11]郭瑞,李素敏,陳婭男,等.基于 SBAS-InSAR 的礦區(qū)采空區(qū)潛在滑坡綜合識(shí)別方法[J].地球信息科學(xué)學(xué)報(bào),2019,21(7):1109-1120.

    猜你喜歡
    試驗(yàn)區(qū)滑坡聚類
    滑坡推力隱式解與顯式解對(duì)比分析——以河北某膨脹土滑坡為例
    18個(gè)自貿(mào)試驗(yàn)區(qū)
    商周刊(2019年19期)2019-10-17 03:23:44
    各試驗(yàn)區(qū)先行先試 探索創(chuàng)造模范做法
    我國(guó)自由貿(mào)易試驗(yàn)區(qū)發(fā)展現(xiàn)狀以及未來展望
    基于DBSACN聚類算法的XML文檔聚類
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    4個(gè)自貿(mào)試驗(yàn)區(qū)總體方案實(shí)施率超過90%
    大社會(huì)(2016年6期)2016-05-04 03:42:08
    基于Fluent的滑坡入水過程數(shù)值模擬
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    基于改進(jìn)的遺傳算法的模糊聚類算法
    2021天堂中文幕一二区在线观| 日韩av在线大香蕉| 男人和女人高潮做爰伦理| 一卡2卡三卡四卡精品乱码亚洲| 欧美高清成人免费视频www| 久久精品91蜜桃| 国产高清激情床上av| 嫩草影院新地址| 国产人妻一区二区三区在| 国产aⅴ精品一区二区三区波| 国产毛片a区久久久久| 午夜影院日韩av| 淫秽高清视频在线观看| 熟女人妻精品中文字幕| 国产在线精品亚洲第一网站| 国产一级毛片七仙女欲春2| 一区福利在线观看| 免费人成视频x8x8入口观看| 蜜桃久久精品国产亚洲av| 国产免费男女视频| 亚洲精品成人久久久久久| 99国产综合亚洲精品| 精品一区二区三区视频在线| 可以在线观看毛片的网站| 永久网站在线| 香蕉av资源在线| 啪啪无遮挡十八禁网站| 尤物成人国产欧美一区二区三区| 小说图片视频综合网站| 精品乱码久久久久久99久播| 最近在线观看免费完整版| 亚洲天堂国产精品一区在线| 亚洲人成网站在线播放欧美日韩| 宅男免费午夜| 在现免费观看毛片| 十八禁国产超污无遮挡网站| 国产精品一区二区性色av| 亚洲av免费在线观看| 99riav亚洲国产免费| 黄色丝袜av网址大全| 一个人看的www免费观看视频| 一级a爱片免费观看的视频| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 天堂动漫精品| 成年人黄色毛片网站| av在线老鸭窝| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 91久久精品国产一区二区成人| 精品无人区乱码1区二区| 成人一区二区视频在线观看| 99国产极品粉嫩在线观看| 国产精品爽爽va在线观看网站| 一进一出抽搐动态| 国产午夜精品论理片| 午夜福利18| 黄色视频,在线免费观看| 国产中年淑女户外野战色| 波多野结衣巨乳人妻| 久久热精品热| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 亚洲av二区三区四区| 草草在线视频免费看| 日韩av在线大香蕉| 亚洲在线自拍视频| 欧美一区二区精品小视频在线| 国产三级黄色录像| 欧美色欧美亚洲另类二区| 日日摸夜夜添夜夜添av毛片 | 亚洲成人久久性| 亚洲熟妇中文字幕五十中出| 非洲黑人性xxxx精品又粗又长| 美女被艹到高潮喷水动态| 国产高潮美女av| 成年女人永久免费观看视频| 亚洲一区高清亚洲精品| 国产黄片美女视频| 日韩欧美一区二区三区在线观看| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 久久精品人妻少妇| 搡老妇女老女人老熟妇| 极品教师在线视频| 亚洲专区国产一区二区| 国产成年人精品一区二区| 一边摸一边抽搐一进一小说| 国产色婷婷99| 色哟哟哟哟哟哟| 99在线人妻在线中文字幕| 亚洲国产色片| 91狼人影院| 精品人妻视频免费看| 我要搜黄色片| 天堂影院成人在线观看| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区三区四区久久| 麻豆成人午夜福利视频| 色5月婷婷丁香| 日韩欧美国产一区二区入口| 51午夜福利影视在线观看| 亚洲精品456在线播放app | av中文乱码字幕在线| 我要看日韩黄色一级片| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 精品久久国产蜜桃| 欧美黄色淫秽网站| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看的www视频| 午夜精品久久久久久毛片777| 日日夜夜操网爽| 波多野结衣高清无吗| 国产精品女同一区二区软件 | 人人妻人人看人人澡| 亚洲成人精品中文字幕电影| 一个人观看的视频www高清免费观看| 一本精品99久久精品77| 国产爱豆传媒在线观看| 一区二区三区四区激情视频 | 久久精品综合一区二区三区| 日本黄色片子视频| 久久香蕉精品热| 亚洲av成人av| 免费看日本二区| 三级男女做爰猛烈吃奶摸视频| 高清日韩中文字幕在线| 午夜激情欧美在线| 变态另类成人亚洲欧美熟女| 淫秽高清视频在线观看| 在线播放国产精品三级| 波多野结衣高清无吗| 91在线精品国自产拍蜜月| 中出人妻视频一区二区| 精品一区二区三区人妻视频| 成人特级av手机在线观看| 欧美3d第一页| 全区人妻精品视频| or卡值多少钱| 国产又黄又爽又无遮挡在线| 日韩欧美国产一区二区入口| 久久精品综合一区二区三区| 免费在线观看亚洲国产| 亚洲,欧美,日韩| 亚洲一区二区三区色噜噜| 一个人免费在线观看的高清视频| 别揉我奶头 嗯啊视频| 国产精品一区二区性色av| 嫩草影院精品99| 蜜桃亚洲精品一区二区三区| 我的老师免费观看完整版| 嫩草影院精品99| 国产精品久久电影中文字幕| 国产免费一级a男人的天堂| 乱人视频在线观看| 又黄又爽又刺激的免费视频.| 嫩草影视91久久| 日本a在线网址| 中文在线观看免费www的网站| 不卡一级毛片| 岛国在线免费视频观看| 成人一区二区视频在线观看| 国产高清激情床上av| a在线观看视频网站| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 国产精品嫩草影院av在线观看 | 亚洲色图av天堂| 一二三四社区在线视频社区8| 狂野欧美白嫩少妇大欣赏| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产清高在天天线| 麻豆av噜噜一区二区三区| 又黄又爽又免费观看的视频| 深夜精品福利| 欧美极品一区二区三区四区| 精品人妻偷拍中文字幕| 天天躁日日操中文字幕| netflix在线观看网站| 五月伊人婷婷丁香| 精品欧美国产一区二区三| 9191精品国产免费久久| 国内揄拍国产精品人妻在线| 简卡轻食公司| 国产亚洲精品久久久久久毛片| 琪琪午夜伦伦电影理论片6080| 老司机福利观看| 1000部很黄的大片| 18美女黄网站色大片免费观看| 国产欧美日韩精品一区二区| 国产精品av视频在线免费观看| 1024手机看黄色片| h日本视频在线播放| 国产三级黄色录像| 黄色丝袜av网址大全| 久久精品人妻少妇| 青草久久国产| 国产毛片a区久久久久| 天堂√8在线中文| 如何舔出高潮| 亚洲av成人不卡在线观看播放网| 久久性视频一级片| 亚洲成av人片在线播放无| 亚洲综合色惰| 91麻豆av在线| 国产高潮美女av| 91久久精品国产一区二区成人| 国产精品电影一区二区三区| 高清日韩中文字幕在线| 欧美极品一区二区三区四区| 成人毛片a级毛片在线播放| 日日干狠狠操夜夜爽| 欧美潮喷喷水| 色吧在线观看| 特大巨黑吊av在线直播| av在线观看视频网站免费| 十八禁人妻一区二区| 在线播放无遮挡| 国产精品一及| 特大巨黑吊av在线直播| av在线观看视频网站免费| 亚洲欧美清纯卡通| 国产白丝娇喘喷水9色精品| 国产一区二区亚洲精品在线观看| 亚洲人成网站高清观看| 亚洲美女视频黄频| 在线看三级毛片| 国产精品亚洲一级av第二区| 国产黄色小视频在线观看| 久久久精品大字幕| 黄色女人牲交| 欧美成人免费av一区二区三区| 两个人的视频大全免费| 十八禁人妻一区二区| 欧美+日韩+精品| 长腿黑丝高跟| 久久人妻av系列| 亚洲国产精品sss在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| 69人妻影院| 窝窝影院91人妻| 久久久成人免费电影| 久久久久性生活片| 精品人妻1区二区| 成人无遮挡网站| 亚洲一区二区三区不卡视频| 色哟哟·www| 色尼玛亚洲综合影院| 国产精品影院久久| 精品熟女少妇八av免费久了| 日本黄大片高清| 日本 欧美在线| 色5月婷婷丁香| 国产国拍精品亚洲av在线观看| 精品久久久久久成人av| 五月伊人婷婷丁香| 99热这里只有是精品在线观看 | 国产视频内射| 国产一区二区三区视频了| 久久久久久久精品吃奶| 特级一级黄色大片| 黄色视频,在线免费观看| 人妻丰满熟妇av一区二区三区| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 国产精品久久视频播放| 免费无遮挡裸体视频| 偷拍熟女少妇极品色| 天堂影院成人在线观看| 五月伊人婷婷丁香| 老司机午夜福利在线观看视频| 日本免费一区二区三区高清不卡| 18禁裸乳无遮挡免费网站照片| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美日韩高清在线视频| 男女做爰动态图高潮gif福利片| 在线看三级毛片| 成人av一区二区三区在线看| 青草久久国产| 日本黄色视频三级网站网址| 一区二区三区免费毛片| av福利片在线观看| 色播亚洲综合网| 国产熟女xx| 欧美一区二区国产精品久久精品| 窝窝影院91人妻| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 久久精品国产自在天天线| 桃色一区二区三区在线观看| 此物有八面人人有两片| 国产美女午夜福利| 日本免费一区二区三区高清不卡| 国产人妻一区二区三区在| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久人妻精品电影| 国产精品不卡视频一区二区 | 又黄又爽又免费观看的视频| 日日干狠狠操夜夜爽| 91九色精品人成在线观看| 免费高清视频大片| 99热精品在线国产| 亚洲内射少妇av| av在线老鸭窝| 五月玫瑰六月丁香| 淫妇啪啪啪对白视频| 久久这里只有精品中国| 亚洲男人的天堂狠狠| 亚洲自拍偷在线| 非洲黑人性xxxx精品又粗又长| 成人av一区二区三区在线看| 国产人妻一区二区三区在| 午夜福利在线观看免费完整高清在 | 97超视频在线观看视频| 日韩欧美国产一区二区入口| 欧美一区二区国产精品久久精品| x7x7x7水蜜桃| 日韩精品中文字幕看吧| 97人妻精品一区二区三区麻豆| 免费观看的影片在线观看| 国产大屁股一区二区在线视频| 国产免费一级a男人的天堂| 99在线人妻在线中文字幕| 国产视频内射| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 国产精品伦人一区二区| 小说图片视频综合网站| 久久精品综合一区二区三区| 亚洲美女搞黄在线观看 | 亚洲成a人片在线一区二区| 激情在线观看视频在线高清| 少妇的逼好多水| 亚洲第一欧美日韩一区二区三区| 日本黄色视频三级网站网址| 18禁在线播放成人免费| av黄色大香蕉| 婷婷精品国产亚洲av在线| 日韩高清综合在线| 亚洲无线观看免费| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| 国产乱人视频| 国内精品久久久久久久电影| 亚洲精品一区av在线观看| x7x7x7水蜜桃| 无遮挡黄片免费观看| 久久99热6这里只有精品| 日韩免费av在线播放| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区性色av| 在现免费观看毛片| 久久热精品热| 国产视频内射| 亚洲精品亚洲一区二区| 天天躁日日操中文字幕| 久久久久久国产a免费观看| 日本 欧美在线| 极品教师在线视频| 国产亚洲精品久久久久久毛片| 性色av乱码一区二区三区2| 国产精品一区二区免费欧美| 一级a爱片免费观看的视频| 国产成人av教育| 一本综合久久免费| 久久精品国产亚洲av香蕉五月| bbb黄色大片| 亚洲午夜理论影院| 成人毛片a级毛片在线播放| 天天躁日日操中文字幕| 成人精品一区二区免费| 婷婷六月久久综合丁香| 熟女人妻精品中文字幕| 亚洲经典国产精华液单 | 中出人妻视频一区二区| 中文亚洲av片在线观看爽| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 成人国产综合亚洲| 国产精品一区二区三区四区久久| 老司机午夜十八禁免费视频| 白带黄色成豆腐渣| 麻豆国产97在线/欧美| 欧美丝袜亚洲另类 | 看十八女毛片水多多多| 少妇丰满av| 久久午夜亚洲精品久久| 免费av观看视频| 亚洲成人精品中文字幕电影| 久久久久国内视频| 中文字幕久久专区| 色在线成人网| 免费搜索国产男女视频| 男女那种视频在线观看| 99国产极品粉嫩在线观看| 老鸭窝网址在线观看| 国产成+人综合+亚洲专区| 亚洲av日韩精品久久久久久密| 亚洲av不卡在线观看| 在线播放国产精品三级| 久久久久免费精品人妻一区二区| 欧美bdsm另类| 99在线视频只有这里精品首页| av专区在线播放| 日韩欧美三级三区| 国产日本99.免费观看| 欧美激情久久久久久爽电影| 嫩草影院精品99| 成熟少妇高潮喷水视频| 三级国产精品欧美在线观看| 毛片一级片免费看久久久久 | 99热这里只有精品一区| 美女大奶头视频| 好男人在线观看高清免费视频| 中文字幕高清在线视频| 国产主播在线观看一区二区| 国产v大片淫在线免费观看| 老司机深夜福利视频在线观看| 色吧在线观看| 成人av在线播放网站| 婷婷亚洲欧美| 麻豆国产97在线/欧美| 欧美最黄视频在线播放免费| 久久久成人免费电影| 久久99热6这里只有精品| 国产精品影院久久| 午夜精品在线福利| 老司机深夜福利视频在线观看| 成人特级av手机在线观看| 亚洲久久久久久中文字幕| 在线观看66精品国产| 亚洲美女视频黄频| 99热精品在线国产| 亚洲国产精品sss在线观看| 禁无遮挡网站| 亚洲午夜理论影院| 亚洲片人在线观看| 无人区码免费观看不卡| 成人av一区二区三区在线看| 亚洲av熟女| 欧美日韩国产亚洲二区| 亚洲国产欧洲综合997久久,| 特大巨黑吊av在线直播| 亚洲 国产 在线| 亚洲欧美清纯卡通| 亚洲无线观看免费| 色哟哟·www| 日本三级黄在线观看| av国产免费在线观看| 日本黄色片子视频| 国内精品一区二区在线观看| 欧美黑人欧美精品刺激| 中文在线观看免费www的网站| 亚洲自拍偷在线| 国产69精品久久久久777片| 日本成人三级电影网站| 久久久精品欧美日韩精品| 亚洲精品456在线播放app | 一个人看的www免费观看视频| 欧美午夜高清在线| 成年女人毛片免费观看观看9| 中文字幕av成人在线电影| 久久久色成人| 亚洲av五月六月丁香网| 51午夜福利影视在线观看| 美女被艹到高潮喷水动态| 美女大奶头视频| 最新在线观看一区二区三区| 国产av在哪里看| 18禁裸乳无遮挡免费网站照片| 亚洲中文日韩欧美视频| 男女床上黄色一级片免费看| 精品久久久久久成人av| 亚洲乱码一区二区免费版| 久久久国产成人精品二区| 日韩成人在线观看一区二区三区| 午夜激情福利司机影院| 在线天堂最新版资源| 国产高清有码在线观看视频| 亚洲最大成人中文| 亚洲自偷自拍三级| 成年免费大片在线观看| 波野结衣二区三区在线| 国产精品免费一区二区三区在线| 中文字幕精品亚洲无线码一区| 亚洲中文字幕日韩| 少妇的逼好多水| 欧美日韩亚洲国产一区二区在线观看| 动漫黄色视频在线观看| 亚洲狠狠婷婷综合久久图片| 丁香六月欧美| 亚洲av成人av| 中文字幕高清在线视频| av在线天堂中文字幕| 亚洲av中文字字幕乱码综合| 日韩欧美三级三区| 国产成人a区在线观看| 久久人人爽人人爽人人片va | 亚洲精品亚洲一区二区| 成年人黄色毛片网站| 三级国产精品欧美在线观看| 在现免费观看毛片| 国产乱人伦免费视频| 听说在线观看完整版免费高清| 欧美极品一区二区三区四区| 国产野战对白在线观看| 99久久精品国产亚洲精品| 日韩欧美 国产精品| 性色avwww在线观看| 国产精品亚洲一级av第二区| or卡值多少钱| 中文字幕高清在线视频| 亚洲av免费在线观看| 成人鲁丝片一二三区免费| 又紧又爽又黄一区二区| 波多野结衣高清作品| 一进一出抽搐动态| 亚洲专区国产一区二区| 欧美性感艳星| 亚洲欧美日韩高清在线视频| eeuss影院久久| 一进一出抽搐动态| 成人毛片a级毛片在线播放| 亚洲人与动物交配视频| 大型黄色视频在线免费观看| 一级黄片播放器| 精品人妻视频免费看| 美女高潮的动态| 国产免费一级a男人的天堂| 男人和女人高潮做爰伦理| 三级国产精品欧美在线观看| 午夜免费激情av| 夜夜看夜夜爽夜夜摸| 亚洲成av人片免费观看| 精品一区二区三区人妻视频| 精品一区二区三区视频在线| 99久久久亚洲精品蜜臀av| 婷婷精品国产亚洲av在线| 国产日本99.免费观看| 日韩成人在线观看一区二区三区| 男人狂女人下面高潮的视频| 久久精品国产自在天天线| 国产精品久久电影中文字幕| 亚洲激情在线av| 九色国产91popny在线| 波多野结衣高清作品| 麻豆一二三区av精品| 国产精品自产拍在线观看55亚洲| 搞女人的毛片| 老鸭窝网址在线观看| 午夜两性在线视频| 99精品久久久久人妻精品| 日韩欧美国产一区二区入口| 免费看日本二区| 久久这里只有精品中国| 又黄又爽又刺激的免费视频.| 亚洲国产欧洲综合997久久,| 国产伦精品一区二区三区四那| 欧美性猛交黑人性爽| 亚洲色图av天堂| 亚洲人成伊人成综合网2020| 黄片小视频在线播放| 天天一区二区日本电影三级| 精品熟女少妇八av免费久了| 日韩欧美一区二区三区在线观看| 欧美日韩国产亚洲二区| 亚洲av免费高清在线观看| www日本黄色视频网| 天堂av国产一区二区熟女人妻| 中国美女看黄片| 久久国产乱子伦精品免费另类| 免费在线观看日本一区| 夜夜爽天天搞| av在线观看视频网站免费| 久久久精品大字幕| 啪啪无遮挡十八禁网站| 亚洲精品一区av在线观看| 国产色婷婷99| 嫩草影院入口| 国产野战对白在线观看| 脱女人内裤的视频| 嫩草影院入口| 久久人人精品亚洲av| 99精品在免费线老司机午夜| 午夜福利在线观看吧| 无遮挡黄片免费观看| 国产av麻豆久久久久久久| 亚洲精品久久国产高清桃花| 午夜福利18| 欧美一级a爱片免费观看看| 久久亚洲真实| 亚洲美女黄片视频| 啦啦啦观看免费观看视频高清| 精品久久国产蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 久久久久九九精品影院| 精品久久久久久久久亚洲 | 99热精品在线国产| 午夜福利免费观看在线| 国内精品一区二区在线观看| 丁香欧美五月| 高清毛片免费观看视频网站| 久久欧美精品欧美久久欧美| 内地一区二区视频在线| 久久人人爽人人爽人人片va | 我的老师免费观看完整版| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 级片在线观看|