• <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)的遺傳算法的模糊聚類算法
    国产v大片淫在线免费观看| a级毛片a级免费在线| 日韩中文字幕欧美一区二区| 亚洲精品456在线播放app | 中文字幕高清在线视频| 精品一区二区三区视频在线 | ponron亚洲| 中文字幕高清在线视频| 亚洲av电影不卡..在线观看| 国产成人精品久久二区二区91| 免费电影在线观看免费观看| 1000部很黄的大片| 中文字幕av在线有码专区| 18禁黄网站禁片午夜丰满| 天天添夜夜摸| 午夜免费成人在线视频| 婷婷亚洲欧美| 亚洲人成网站高清观看| 久久久久性生活片| 欧美日韩乱码在线| 亚洲成人中文字幕在线播放| 国产精品一区二区精品视频观看| 中文字幕久久专区| 欧美日韩一级在线毛片| 天天添夜夜摸| 国产私拍福利视频在线观看| 日本免费一区二区三区高清不卡| 亚洲av成人一区二区三| www国产在线视频色| 又大又爽又粗| 国产一区在线观看成人免费| 中文字幕久久专区| 亚洲人与动物交配视频| 日本五十路高清| 免费在线观看日本一区| 毛片女人毛片| 露出奶头的视频| 男人舔女人下体高潮全视频| 女同久久另类99精品国产91| 久久久久免费精品人妻一区二区| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久末码| 久久99热这里只有精品18| 精品国产超薄肉色丝袜足j| 日韩欧美在线乱码| 久久国产精品人妻蜜桃| 国内精品久久久久久久电影| 这个男人来自地球电影免费观看| 日本免费a在线| 91麻豆精品激情在线观看国产| 国内少妇人妻偷人精品xxx网站 | 日本与韩国留学比较| 久久久久国产精品人妻aⅴ院| 伦理电影免费视频| 精品免费久久久久久久清纯| 亚洲中文字幕一区二区三区有码在线看 | 特大巨黑吊av在线直播| 成人无遮挡网站| 舔av片在线| 日韩 欧美 亚洲 中文字幕| 好男人电影高清在线观看| 手机成人av网站| 日本一本二区三区精品| 男人的好看免费观看在线视频| 精品国产乱码久久久久久男人| 在线永久观看黄色视频| 男人舔奶头视频| 白带黄色成豆腐渣| 中文字幕高清在线视频| 国产毛片a区久久久久| 脱女人内裤的视频| 国产单亲对白刺激| 后天国语完整版免费观看| 日韩av在线大香蕉| 88av欧美| 国产精品久久久av美女十八| 日本黄大片高清| www日本在线高清视频| 欧美大码av| 99热6这里只有精品| 亚洲欧美一区二区三区黑人| 免费高清视频大片| 精品不卡国产一区二区三区| 九色国产91popny在线| 国产日本99.免费观看| 久久午夜亚洲精品久久| 欧美在线一区亚洲| 视频区欧美日本亚洲| 欧美午夜高清在线| 欧美在线黄色| 男人舔女人的私密视频| 黑人巨大精品欧美一区二区mp4| 日本成人三级电影网站| 欧美不卡视频在线免费观看| 91九色精品人成在线观看| 人妻丰满熟妇av一区二区三区| 91麻豆av在线| 久99久视频精品免费| 午夜精品在线福利| www.www免费av| 一边摸一边抽搐一进一小说| 亚洲欧美日韩卡通动漫| 久久精品影院6| 国产伦精品一区二区三区视频9 | 成人性生交大片免费视频hd| 国产精品久久久久久人妻精品电影| 国产精品 国内视频| 国产亚洲欧美在线一区二区| 国产精品,欧美在线| 国产免费av片在线观看野外av| 国产男靠女视频免费网站| 国语自产精品视频在线第100页| 日韩欧美在线乱码| 一区二区三区国产精品乱码| 久久久久精品国产欧美久久久| 久久精品国产99精品国产亚洲性色| 一级作爱视频免费观看| 两个人视频免费观看高清| 国产午夜精品久久久久久| 男女做爰动态图高潮gif福利片| 亚洲精品在线观看二区| 精品国产亚洲在线| 在线十欧美十亚洲十日本专区| 99久久99久久久精品蜜桃| 啦啦啦韩国在线观看视频| 国产高清videossex| 午夜成年电影在线免费观看| 欧美最黄视频在线播放免费| 国产黄片美女视频| 一本综合久久免费| 国产成人精品久久二区二区免费| 亚洲av中文字字幕乱码综合| 成人欧美大片| 天堂网av新在线| 丰满人妻熟妇乱又伦精品不卡| 精品99又大又爽又粗少妇毛片 | 免费观看人在逋| 亚洲中文av在线| 变态另类丝袜制服| 欧美日韩乱码在线| 我要搜黄色片| 精品国产乱码久久久久久男人| 日本成人三级电影网站| 国产伦一二天堂av在线观看| 好看av亚洲va欧美ⅴa在| 搞女人的毛片| 天天一区二区日本电影三级| 成人国产综合亚洲| 人人妻人人看人人澡| 亚洲国产中文字幕在线视频| 好男人在线观看高清免费视频| 亚洲欧美日韩卡通动漫| 成在线人永久免费视频| 男女床上黄色一级片免费看| 神马国产精品三级电影在线观看| 黄色成人免费大全| 色老头精品视频在线观看| 99热精品在线国产| 成人亚洲精品av一区二区| 亚洲成av人片免费观看| 男女午夜视频在线观看| 又黄又粗又硬又大视频| 亚洲专区字幕在线| 日日摸夜夜添夜夜添小说| 级片在线观看| 国产探花在线观看一区二区| 亚洲av电影不卡..在线观看| 国产成人av教育| 久久天躁狠狠躁夜夜2o2o| 免费大片18禁| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av香蕉五月| 男女午夜视频在线观看| 国产亚洲精品久久久久久毛片| 久久精品国产综合久久久| 久久久国产成人精品二区| 久久久国产成人免费| 狠狠狠狠99中文字幕| 久久久久久久精品吃奶| 欧美成人性av电影在线观看| 国产伦人伦偷精品视频| 88av欧美| 国产午夜精品久久久久久| 91麻豆精品激情在线观看国产| av国产免费在线观看| 亚洲乱码一区二区免费版| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 国产亚洲精品久久久久久毛片| 久久精品国产综合久久久| 国产亚洲精品一区二区www| 一个人看的www免费观看视频| 成人国产综合亚洲| 亚洲国产高清在线一区二区三| 亚洲中文字幕日韩| 国产精品99久久99久久久不卡| 国产激情偷乱视频一区二区| 午夜日韩欧美国产| 欧美精品啪啪一区二区三区| 成人国产综合亚洲| 网址你懂的国产日韩在线| 国产精品98久久久久久宅男小说| 动漫黄色视频在线观看| 欧美一区二区精品小视频在线| 黄色 视频免费看| 又粗又爽又猛毛片免费看| 久久人人精品亚洲av| 黄色女人牲交| 久久久久久久久久黄片| 国产三级黄色录像| av视频在线观看入口| 免费观看的影片在线观看| 精品久久久久久久久久免费视频| 日韩欧美 国产精品| 国产午夜精品久久久久久| 精品无人区乱码1区二区| 搞女人的毛片| 女人高潮潮喷娇喘18禁视频| www.www免费av| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 一个人观看的视频www高清免费观看 | 国产麻豆成人av免费视频| 91九色精品人成在线观看| 亚洲熟女毛片儿| 无遮挡黄片免费观看| 久久精品aⅴ一区二区三区四区| 国产精品,欧美在线| 又大又爽又粗| 国产精品电影一区二区三区| 一级a爱片免费观看的视频| 成年人黄色毛片网站| 可以在线观看毛片的网站| 欧美又色又爽又黄视频| 欧美黑人巨大hd| 麻豆一二三区av精品| 97碰自拍视频| 真人一进一出gif抽搐免费| 国产私拍福利视频在线观看| 亚洲国产色片| 三级男女做爰猛烈吃奶摸视频| 琪琪午夜伦伦电影理论片6080| 性欧美人与动物交配| 999精品在线视频| 最近最新中文字幕大全免费视频| 国产欧美日韩精品一区二区| 国产成+人综合+亚洲专区| 一级a爱片免费观看的视频| 欧美日韩福利视频一区二区| 亚洲精品在线美女| 18美女黄网站色大片免费观看| 人人妻人人看人人澡| 久久中文看片网| 少妇的丰满在线观看| 三级国产精品欧美在线观看 | 国产精品一区二区免费欧美| 日韩 欧美 亚洲 中文字幕| 黄色片一级片一级黄色片| 久久久久久久久中文| 久久婷婷人人爽人人干人人爱| xxxwww97欧美| 熟女电影av网| 午夜福利在线观看吧| 欧美大码av| 2021天堂中文幕一二区在线观| 亚洲国产欧美一区二区综合| 欧美绝顶高潮抽搐喷水| 男女做爰动态图高潮gif福利片| 日韩高清综合在线| 欧美激情久久久久久爽电影| 国产亚洲欧美在线一区二区| 亚洲欧美日韩高清在线视频| 18禁黄网站禁片免费观看直播| 国产美女午夜福利| 精品福利观看| АⅤ资源中文在线天堂| 老司机福利观看| 999久久久精品免费观看国产| 国产成人精品无人区| 欧美一级a爱片免费观看看| 嫩草影视91久久| 18禁裸乳无遮挡免费网站照片| 国产激情欧美一区二区| 午夜a级毛片| 国产精华一区二区三区| 一二三四在线观看免费中文在| 床上黄色一级片| 亚洲精品456在线播放app | 久久精品综合一区二区三区| 成人亚洲精品av一区二区| 色哟哟哟哟哟哟| 久久精品91无色码中文字幕| 国产av麻豆久久久久久久| 巨乳人妻的诱惑在线观看| 欧美乱妇无乱码| 一级毛片高清免费大全| 成年女人永久免费观看视频| 91麻豆av在线| 日本免费一区二区三区高清不卡| 久久亚洲真实| 亚洲欧美精品综合久久99| 精品久久久久久久人妻蜜臀av| 国产精品99久久久久久久久| 久久精品人妻少妇| 日韩欧美三级三区| 69av精品久久久久久| 三级毛片av免费| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 国产精品久久久久久久电影 | 久久精品亚洲精品国产色婷小说| 91字幕亚洲| 18禁黄网站禁片免费观看直播| 长腿黑丝高跟| 国产成人精品无人区| 亚洲一区高清亚洲精品| xxx96com| 手机成人av网站| 欧美又色又爽又黄视频| 久久久久久大精品| 久久热在线av| 在线视频色国产色| 女警被强在线播放| 日本黄色视频三级网站网址| 禁无遮挡网站| 国产亚洲精品一区二区www| 国产伦精品一区二区三区视频9 | 黑人操中国人逼视频| 日韩国内少妇激情av| 又粗又爽又猛毛片免费看| 1024香蕉在线观看| 18禁黄网站禁片午夜丰满| 在线观看66精品国产| 欧美性猛交╳xxx乱大交人| 久久热在线av| 午夜福利18| 给我免费播放毛片高清在线观看| 精品国产乱子伦一区二区三区| 成人国产一区最新在线观看| 狠狠狠狠99中文字幕| 夜夜夜夜夜久久久久| 国产成人av激情在线播放| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 男人的好看免费观看在线视频| 精品午夜福利视频在线观看一区| 又黄又爽又免费观看的视频| 久久香蕉国产精品| e午夜精品久久久久久久| 黑人操中国人逼视频| 午夜免费成人在线视频| 美女高潮喷水抽搐中文字幕| 久久精品91无色码中文字幕| 国产亚洲精品一区二区www| 丰满人妻一区二区三区视频av | 一个人免费在线观看电影 | 日韩欧美三级三区| 俺也久久电影网| 国产又黄又爽又无遮挡在线| 小说图片视频综合网站| 亚洲国产精品sss在线观看| 国产高清有码在线观看视频| 狂野欧美白嫩少妇大欣赏| 亚洲第一电影网av| 亚洲中文字幕日韩| 五月玫瑰六月丁香| 亚洲精品久久国产高清桃花| 欧美一级a爱片免费观看看| 美女午夜性视频免费| 亚洲国产色片| 老司机深夜福利视频在线观看| 国产精品精品国产色婷婷| 日韩欧美国产在线观看| 搡老岳熟女国产| 国产麻豆成人av免费视频| 在线观看免费视频日本深夜| 亚洲成av人片免费观看| 在线a可以看的网站| 国产精品女同一区二区软件 | 亚洲欧美精品综合一区二区三区| 精品免费久久久久久久清纯| 精品国产乱子伦一区二区三区| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久 | 欧美又色又爽又黄视频| 日本 欧美在线| 欧美丝袜亚洲另类 | 亚洲国产看品久久| 久久久精品欧美日韩精品| 狂野欧美白嫩少妇大欣赏| 国产熟女xx| 亚洲精品美女久久av网站| 99热这里只有精品一区 | 欧美午夜高清在线| a在线观看视频网站| 美女扒开内裤让男人捅视频| 欧美激情久久久久久爽电影| 国产黄a三级三级三级人| 亚洲男人的天堂狠狠| 午夜福利18| 热99re8久久精品国产| 不卡一级毛片| 午夜激情欧美在线| 国产精品综合久久久久久久免费| 亚洲五月婷婷丁香| 亚洲色图av天堂| 国产精品自产拍在线观看55亚洲| 日本一本二区三区精品| 国产一区二区在线观看日韩 | 国产人伦9x9x在线观看| 欧美乱色亚洲激情| 国产一区二区在线观看日韩 | 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 法律面前人人平等表现在哪些方面| 亚洲中文日韩欧美视频| 黑人巨大精品欧美一区二区mp4| 亚洲乱码一区二区免费版| 天堂√8在线中文| 最近最新免费中文字幕在线| tocl精华| 国产精品永久免费网站| 又黄又粗又硬又大视频| 免费av不卡在线播放| www日本黄色视频网| 制服丝袜大香蕉在线| 一本一本综合久久| 国产精品 欧美亚洲| 国内精品一区二区在线观看| 亚洲美女黄片视频| 久久久水蜜桃国产精品网| 久久久精品欧美日韩精品| 欧美日韩黄片免| 一进一出抽搐gif免费好疼| 亚洲精品色激情综合| 免费观看的影片在线观看| 国产成人精品久久二区二区免费| 伊人久久大香线蕉亚洲五| 哪里可以看免费的av片| 日本 欧美在线| 级片在线观看| av在线蜜桃| 成人特级av手机在线观看| 18禁黄网站禁片午夜丰满| 韩国av一区二区三区四区| 国产一级毛片七仙女欲春2| 久久婷婷人人爽人人干人人爱| 亚洲七黄色美女视频| 免费看日本二区| 日本熟妇午夜| 国产不卡一卡二| 亚洲精品中文字幕一二三四区| 日本黄色片子视频| 久久精品国产综合久久久| 久久精品人妻少妇| 一二三四社区在线视频社区8| 欧美乱妇无乱码| 看免费av毛片| 麻豆久久精品国产亚洲av| 亚洲18禁久久av| 桃红色精品国产亚洲av| 真实男女啪啪啪动态图| av福利片在线观看| 欧美av亚洲av综合av国产av| 国产亚洲欧美在线一区二区| 亚洲欧美激情综合另类| 久久精品国产综合久久久| 欧美日本亚洲视频在线播放| 男女之事视频高清在线观看| 观看免费一级毛片| 日韩精品青青久久久久久| 国产黄色小视频在线观看| 国产亚洲欧美98| 色老头精品视频在线观看| 精品电影一区二区在线| 美女大奶头视频| 久久久久久大精品| 两人在一起打扑克的视频| 搞女人的毛片| 少妇人妻一区二区三区视频| 亚洲精品中文字幕一二三四区| 成人亚洲精品av一区二区| 天堂av国产一区二区熟女人妻| 国内少妇人妻偷人精品xxx网站 | 欧美成人一区二区免费高清观看 | 天堂动漫精品| 99热6这里只有精品| 不卡av一区二区三区| 老汉色av国产亚洲站长工具| 全区人妻精品视频| 免费看美女性在线毛片视频| 午夜成年电影在线免费观看| 女同久久另类99精品国产91| 国产精品一及| 一本一本综合久久| 日韩欧美在线乱码| 91av网一区二区| or卡值多少钱| 成年免费大片在线观看| 曰老女人黄片| 免费电影在线观看免费观看| 香蕉国产在线看| 99热精品在线国产| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放| 特大巨黑吊av在线直播| 国产精品99久久99久久久不卡| 深夜精品福利| 亚洲成av人片免费观看| 国产伦精品一区二区三区四那| 又粗又爽又猛毛片免费看| 黄频高清免费视频| 亚洲人成电影免费在线| 欧美精品啪啪一区二区三区| 草草在线视频免费看| 男女视频在线观看网站免费| 中文在线观看免费www的网站| 国产精品久久久久久亚洲av鲁大| 亚洲欧美日韩卡通动漫| 中文字幕久久专区| 一区二区三区高清视频在线| 欧美黄色片欧美黄色片| 真人一进一出gif抽搐免费| 黑人操中国人逼视频| av国产免费在线观看| 国产成年人精品一区二区| 俄罗斯特黄特色一大片| 狂野欧美激情性xxxx| 国产精品久久久久久久电影 | 在线国产一区二区在线| 五月玫瑰六月丁香| 亚洲成人久久爱视频| 91av网站免费观看| 精品久久久久久久久久免费视频| 一区福利在线观看| 老司机午夜十八禁免费视频| 国内精品美女久久久久久| www.www免费av| 老司机在亚洲福利影院| 男女之事视频高清在线观看| 欧美乱码精品一区二区三区| 亚洲熟妇熟女久久| 精品久久久久久久久久免费视频| 国产精品久久久人人做人人爽| 日本撒尿小便嘘嘘汇集6| 国产三级中文精品| 两个人视频免费观看高清| 久久久精品大字幕| 亚洲成人精品中文字幕电影| 很黄的视频免费| 三级国产精品欧美在线观看 | 欧美丝袜亚洲另类 | 亚洲欧美日韩东京热| 久久久久久大精品| 亚洲在线观看片| 亚洲激情在线av| 成熟少妇高潮喷水视频| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 一区二区三区激情视频| 99久久无色码亚洲精品果冻| 超碰成人久久| 黄色视频,在线免费观看| 日日摸夜夜添夜夜添小说| 久久久久久久久久黄片| 亚洲国产精品久久男人天堂| 成人三级黄色视频| 中文字幕精品亚洲无线码一区| 身体一侧抽搐| 午夜免费激情av| 又黄又爽又免费观看的视频| 观看美女的网站| 欧美日本视频| 九九热线精品视视频播放| 国内精品久久久久精免费| 色吧在线观看| 日本三级黄在线观看| 国产精品香港三级国产av潘金莲| 性欧美人与动物交配| 最好的美女福利视频网| 在线免费观看不下载黄p国产 | 国产1区2区3区精品| 五月玫瑰六月丁香| 国产综合懂色| 免费av毛片视频| 欧美极品一区二区三区四区| 日韩精品青青久久久久久| 在线观看日韩欧美| 色老头精品视频在线观看| 88av欧美| 熟女少妇亚洲综合色aaa.| 亚洲国产精品999在线| 日本一二三区视频观看| 超碰成人久久| 变态另类丝袜制服| 看免费av毛片| 国产一区二区在线观看日韩 | 日本黄大片高清| 欧美成人一区二区免费高清观看 | 亚洲自偷自拍图片 自拍| 免费观看人在逋| 成人一区二区视频在线观看| 午夜a级毛片| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久亚洲av鲁大| 两性夫妻黄色片| 成人高潮视频无遮挡免费网站| 变态另类丝袜制服| 欧美最黄视频在线播放免费| 美女午夜性视频免费| 精品无人区乱码1区二区| 成年女人看的毛片在线观看| 久久性视频一级片| 老司机在亚洲福利影院|