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

    信息量與多模型耦合的滑坡易發(fā)性評價研究

    2021-12-09 02:43:24儀政宋琨黃海峰李輝
    人民長江 2021年10期
    關(guān)鍵詞:三峽庫區(qū)

    儀政 宋琨 黃海峰 李輝

    摘要:滑坡易發(fā)性評價是識別滑坡災(zāi)害危險性“靶區(qū)”的基礎(chǔ)和關(guān)鍵。三峽庫區(qū)秭歸縣的侏羅系地層區(qū)是滑坡災(zāi)害易發(fā)區(qū),選取坡度、地表切割深度、曲率、距水系與路網(wǎng)距離等8個影響因子,通過信息量模型(I)與人工神經(jīng)網(wǎng)絡(luò)(ANN)、隨機森林(RF)、XGBboost、支持向量機(SVM)等4種模型耦合,進行了滑坡易發(fā)性評價。通過受試者工作特征(ROC)曲線對比分析發(fā)現(xiàn):信息量和支持向量機耦合模型的ROC曲線AUC值最大(0.848),滑坡易發(fā)性的分區(qū)圖顯示高易發(fā)區(qū)主要集中在水系和路網(wǎng)兩側(cè),其結(jié)果與實際滑坡分布基本一致。研究成果可為區(qū)域滑坡易發(fā)性評價和滑坡災(zāi)害的防治提供參考。

    關(guān)鍵詞:滑坡易發(fā)性評價; 滑坡易發(fā)性分區(qū); 信息量模型; 多模型耦合; 侏羅系地層; 三峽庫區(qū)

    中圖法分類號: P642

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2021.10.022

    0引 言

    滑坡災(zāi)害的頻繁發(fā)生會對山區(qū)居民的生命財產(chǎn)安全造成巨大威脅。為降低滑坡災(zāi)害的損失,需要在滑坡易發(fā)“靶區(qū)”識別的基礎(chǔ)上進行合理的工程防治或監(jiān)測預(yù)警。而滑坡易發(fā)“靶區(qū)”的識別主要是基于滑坡易發(fā)性的評價結(jié)果而進行的。因此,滑坡易發(fā)性的評價成為山區(qū)地質(zhì)災(zāi)害防治的基礎(chǔ)性工作。

    滑坡易發(fā)性評價的方法主要有:信息量法[1]、確定性系數(shù)[2-3]、邏輯回歸[2-3]、層次分析法[4]、熵值法[5]、模糊度[6]、支持向量機[7]、人工神經(jīng)網(wǎng)絡(luò)[7]、隨機森林[8]、決策樹[9]、XGBoost[10]等。目前國內(nèi)外學者對此進行了較多研究,如Devkota[2]等運用確定性因子、熵值法、邏輯回歸模型對喜馬拉雅山區(qū)滑坡進行了易發(fā)性評價,并分析了各模型的準確度;喬德京等[7]將隨機森林和確定系數(shù)耦合,進行滑坡易發(fā)性評價,發(fā)現(xiàn)耦合后的模型精度更高;張玘愷[3]等將信息量法、確定性系數(shù)、邏輯回歸相耦合,對九寨溝縣滑坡進行易發(fā)性評價,并指出耦合后的模型要比單獨信息量模型的評價準確高。由于傳統(tǒng)信息量法無法體現(xiàn)各因子的權(quán)重,因此,將傳統(tǒng)信息量和確定性系數(shù)模型及其他模型耦合[3-8],可賦予信息量和確定系數(shù)相應(yīng)的權(quán)重,使得滑坡易發(fā)性評價獲得更高的準確度和可靠度。

    三峽庫區(qū)秭歸縣的侏羅系地層區(qū)(占秭歸縣土地面積的22%)是滑坡地質(zhì)災(zāi)害的多發(fā)區(qū)。李長冬[11]等以三峽庫區(qū)侏羅系軟硬相間地層作為研究對象,并對三峽庫區(qū)侏羅系軟硬相間地層中抗滑樁布設(shè)和結(jié)構(gòu)優(yōu)化進行了研究。據(jù)統(tǒng)計,秭歸縣39%的滑坡災(zāi)害發(fā)生在該地層區(qū)。為評價該區(qū)域滑坡的易發(fā)性,采用信息量模型(I)與人工神經(jīng)網(wǎng)絡(luò)(ANN)、隨機森林(RF)、XGBoost、支持向量機(SVM)等4種模型耦合,并對比各耦合模型的準確度尋找最優(yōu)模型,以為其他區(qū)域滑坡易發(fā)性評價提供借鑒。

    1研究方法

    1.1信息量法

    信息量法是以信息論為基礎(chǔ)的統(tǒng)計方法。信息量法把影響區(qū)域穩(wěn)定性因素的數(shù)值轉(zhuǎn)化為能夠反映出對區(qū)域穩(wěn)定影響程度的信息量值,信息量越大,滑坡發(fā)生的可能性越大,反之亦然。信息量模型計算公式如下[1]:

    1.3隨機森林

    隨機森林屬于一種特殊的自助抽樣集成(Bagging)方法,并以決策樹作為Bagging中的模型[12]。該模型對大小為N的原始樣本進行有放回的隨機抽樣,得到m個大小為N的新訓練集,之后對每個訓練集進行決策樹建模,然后對這些決策樹進行組合,最終采用投票的方式以少數(shù)服從多數(shù)的原則來得到最終的分類和預(yù)測結(jié)果。隨機森林模型不僅對數(shù)據(jù)選取隨機化,且對特征的選取進行了隨機化,通過平均決策樹極大減少了模型的過擬合。

    1.4XGBoost

    XGBoost模型是改進的梯度提升樹(GBDT)算法,模型使用牛頓法求解損失函數(shù)的極值,將損失函數(shù)通過泰勒形式展開到二階,另外損失函數(shù)中加入了正則化項。訓練時的目標函數(shù)由兩部分組成,第一部分為梯度提升算法損失,第二部分為正則化項[13]。

    1.5SVM(支持向量機)

    支持向量機(SVM)是Vapnik基于統(tǒng)計學理論提出的一種機器學習方法,通過尋求結(jié)構(gòu)化風險最小化從而提高學習機的泛化能力[14]。為了減小歸納誤差,本文選取徑向基核函數(shù)作為核函數(shù)來進行滑坡易發(fā)性評價。

    徑向基核函數(shù)為

    1.6信息量與多模型耦合模型

    傳統(tǒng)的信息量模型雖能客觀地獲取到各個因子對滑坡災(zāi)害易發(fā)性的影響程度,但是沒有確定各影響因子的權(quán)重,而ANN、RF、XGBoost、SVM模型的可以確定各影響因子的權(quán)重。因此為了更準確地預(yù)測滑坡易發(fā)性地區(qū)區(qū)劃,本文將信息量模型和ANN、RF、XGBoost、SVM模型耦合,進而得到滑坡易發(fā)性區(qū)劃圖。

    2研究區(qū)概況及評價流程

    2.1研究區(qū)概況

    三峽庫區(qū)秭歸縣是滑坡災(zāi)害頻發(fā)的重點區(qū)域。近幾年曾發(fā)生多起大型滑坡災(zāi)害事件,如2003年7月13日的沙鎮(zhèn)溪鎮(zhèn)千將坪滑坡、2014年9月2日的沙鎮(zhèn)溪鎮(zhèn)杉樹槽滑坡、2017年10月29日的歸州鎮(zhèn)鹽關(guān)滑坡和2019年12月5日的泄灘鄉(xiāng)卡門子灣滑坡等。這些新發(fā)生的滑坡災(zāi)害所處地層均為侏羅系地層。因此,將秭歸縣侏羅系地層區(qū)作為研究區(qū)進行滑坡災(zāi)害易發(fā)性評價。

    研究區(qū)地理坐標為110°27′45″E~110°46′53″E,31°11′47″N~30°54′12″N,國土面積532 km2。地層主要為侏羅系的沙溪廟組、聶家山組和香溪組。巖性主要為紫紅色泥巖夾石英砂巖、灰綠色粉砂質(zhì)泥巖、粉砂巖夾長石砂巖和炭質(zhì)頁巖,以及中厚層砂屑灰?guī)r和泥質(zhì)白云巖等[15]。據(jù)調(diào)查數(shù)據(jù)統(tǒng)計可知,區(qū)域內(nèi)共有滑坡災(zāi)害點183處(見圖1),占秭歸縣全境滑坡總數(shù)(458處)的39%,多數(shù)分布在河流兩側(cè)。

    2.2數(shù)據(jù)來源和滑坡易發(fā)性評價流程

    本文研究所需要的數(shù)據(jù)包括遙感影像數(shù)據(jù)、滑坡數(shù)據(jù)、河流路網(wǎng)矢量數(shù)據(jù)等。高程數(shù)據(jù)采用分辨率為30m ASTER GDEM數(shù)據(jù),遙感數(shù)據(jù)使用Landsat4-5 TM衛(wèi)星數(shù)字產(chǎn)品,二者均從地理空間數(shù)據(jù)云(http:∥www.gscloud.cn)下載獲得。滑坡數(shù)據(jù)通過秭歸縣地質(zhì)災(zāi)害點排查資料獲取。本文采用分辨率50 m×50 m的柵格作為滑坡易發(fā)性評價的基本單元,使用ArcGIS將整個研究區(qū)劃分為213 785個柵格單元。首先利用ArcGIS軟件提取各評價因子圖層,然后依據(jù)研究區(qū)內(nèi)滑坡影響因子的分析結(jié)果,按照每個因子的不同類別及對滑坡的影響程度進行因子分級劃分,再將分級后的因子圖層和滑坡災(zāi)害點圖層屬性連接,計算各分級因子信息量;之后為各分級圖層賦予信息量數(shù)據(jù),進而得到各因子信息量圖層。使用GIS中柵格轉(zhuǎn)點工具將研究區(qū)的柵格單元轉(zhuǎn)換成點數(shù)據(jù),之后使用多值提取至點工具與各因子信息量柵格圖層疊加分析,提取到整個研究區(qū)信息量屬性數(shù)據(jù),建立滑坡易發(fā)性數(shù)據(jù)庫。將183個災(zāi)害點和等數(shù)量的非災(zāi)害點作為總樣本,從中隨機選取70%的樣本點作為訓練樣本,導入到各模型當中進行訓練;剩余的30%作為測試樣本,運用ROC曲線來檢驗?zāi)P头诸惖男Ч??;乱装l(fā)性評價流程如圖2所示。

    3指標因子構(gòu)建

    指標因子的選取影響著滑坡易發(fā)性評價的精度。通過對研究區(qū)現(xiàn)場地質(zhì)調(diào)查以及災(zāi)害點的時空分布規(guī)律,選取了以下8種環(huán)境指標因子:坡度、地表切割深度、曲率、地形濕度指數(shù)(TWI)、歸一化植被覆蓋指數(shù)(NDVI)、距水系與路網(wǎng)距離、土地分類。其中,坡體的坡度直接影響到坡體內(nèi)部的應(yīng)力分布狀態(tài),控制坡體的穩(wěn)定性;地表切割深度直觀地反映了地表被侵蝕切割的情況;曲率將影響著坡面物質(zhì)的搬運和沉積;歸一化植被覆蓋指數(shù)反映了地表植被覆蓋情況,植被頻繁的地區(qū)水土保持越好,有利于增強坡體的穩(wěn)定性。研究區(qū)位于三峽庫區(qū),水庫和河流旁的坡體坡腳常年受到?jīng)_刷掏蝕,加速了坡體的失穩(wěn),因此選取距水系距離及地形濕度指數(shù)作為環(huán)境指標因子。土地分類和離公路距離可反映人類活動的區(qū)域,交通建設(shè)削坡改變了坡體的自然地貌,增大了災(zāi)害發(fā)生的概率。

    通過整理收集的資料(高程數(shù)字模型DEM,水系圖,路網(wǎng)圖,遙感影像圖,滑坡點分布圖),利用ArcGIS空間分析功能和柵格計算器處理DEM數(shù)據(jù)得到坡度、地形濕度指數(shù)、地表切割深度、曲率因子圖層。對水系和路網(wǎng)進行緩沖區(qū)分析,緩沖距離均設(shè)置為0~300 m,300~600 m,600~900 m,>900 m。把獲得Landsat4-5 TM遙感影像數(shù)據(jù)導入到ENVI軟件中,利用波段運算得到研究區(qū)歸一化植被覆蓋指數(shù),之后使用監(jiān)督分類工具以最大似然法進行對遙感影像進行土地分類。最終得到的各評價因子圖層如圖3所示。

    4滑坡易發(fā)性評價

    4.1基于ArcGIS信息量計算

    根據(jù)災(zāi)害點和各因子圖層疊加分析,依據(jù)信息量模型計算出不同因子等級下的信息量值(見表1),然后將信息量值賦予各因子分級圖層中,利用ArcGIS柵格計算器將各因子信息量圖層疊加,得到研究區(qū)的滑坡總信息量。

    4.2基于多模型與信息量耦合的滑坡易發(fā)性評價

    選取含有信息量的災(zāi)害點和等數(shù)量的非災(zāi)害點作為總樣本,隨機抽取70%總樣本作為訓練樣本,30%總樣本作為測試數(shù)據(jù)。將樣本數(shù)據(jù)導入到SPSS Modeler以及Python進行訓練,得到ANN,RF,XGBboost,SVM易發(fā)性評價模型,然后將整個研究區(qū)信息量總數(shù)據(jù)代入模型中進行計算,得到研究區(qū)的滑坡易發(fā)性指數(shù)。將獲得的易發(fā)性指數(shù)導入ArcGIS當中,利用自然斷點法劃分5個等級,最終得到多種模型與信息量耦合的研究區(qū)滑坡易發(fā)性評價圖(見圖4)。

    4.3精度分析

    為評價各模型預(yù)測的精度,本文采用受試者工作特征曲線(ROC)對模型進行精度評價。ROC曲線越靠近左上方,表明預(yù)測結(jié)果的準確度越高,曲線與X坐標軸圍成的面積AUC值通常來衡量預(yù)測結(jié)果的精確度,AUC值越大,表明預(yù)測預(yù)測精度越高。ROC曲線縱軸代表真陽性率(敏感性),在滑坡災(zāi)害易發(fā)性評價中代表研究區(qū)中各易發(fā)性等級內(nèi)真實發(fā)生滑坡災(zāi)害的累加百分比;其橫軸代表假陽性率(特異性),在滑坡災(zāi)害易發(fā)性評價中代表研究區(qū)內(nèi)各易發(fā)性等級面積累加百分比[3]。從ROC曲線(見圖5)可知,信息量模型,I-RF模型,I-XGBoost模型,I-ANN模型,I-SVM模型的AUC值分別為0.653,0.784,0.813,0.825,0.848。結(jié)果表明單獨的信息量模型預(yù)測效果遠不如耦合模型的準確度,且I-SVM耦合模型精度最高,耦合效果最好。

    5結(jié) 論

    (1) 本文以三峽庫區(qū)秭歸縣侏羅系地層區(qū)為研究區(qū)域,選取了坡度、曲率、地形濕度指數(shù)(TWI)、歸一化植被覆蓋指數(shù)(NDVI)、地表切割深度、土地分類、水系距離、路網(wǎng)距離共8個影響因子建立滑坡災(zāi)害易發(fā)性評價指標體系,運用單一信息量模型以及信息量模型與多模型耦合模型,建立了滑坡易發(fā)性區(qū)劃圖。單一信息量模型預(yù)測精度不如信息量和機器學習耦合模型,各模型中水系距離和道路距離均占比較大。

    (2) 根據(jù)滑坡易發(fā)性區(qū)劃圖可知,滑坡極高和高易發(fā)區(qū)主要分布在水系和路網(wǎng)的兩側(cè),其結(jié)果與滑坡災(zāi)害點分布基本一致。

    (3) 不同信息量和機器學習耦合模型結(jié)果對比發(fā)現(xiàn),信息量與支持向量機耦合模型(I-SVM模型)的ROC曲線其AUC值最大,為0.848,表明信息量與支持向量機耦合模型的預(yù)測效果最佳。

    參考文獻:

    [1]陳立華,李立豐,吳福,等.基于GIS與信息量法的北流市地質(zhì)災(zāi)害易發(fā)性評價[J].地球與環(huán)境,2020,48(4):471-479.

    [2]DEVKOTA K C,REGMI A D,POURGHASEMI H R,et al.Landslide susceptibility mapping using certainty factor,index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya[J].Natural Hazards,2013,65(1):135-165.

    [3]張玘愷,凌斯祥,李曉寧,等.九寨溝縣滑坡災(zāi)害易發(fā)性快速評估模型對比研究[J].巖石力學與工程學報,2020,39(8):1595-1610.

    [4]許嘉慧,孫德亮.基于GIS與改進層次分析法的奉節(jié)縣滑坡易發(fā)性區(qū)劃[J].重慶師范大學學報(自然科學版),2020,37(2):36-44.

    [5]韓玲,張庭瑜,張恒.基于IOE和SVM模型的府谷鎮(zhèn)滑坡易發(fā)性分區(qū)[J].水土保持研究,2019,26(3):367-372.

    [6]MENG Q,MIAO F,ZHEN J,et al.GIS-based landslide susceptibility mapping with logistic regression,analytical hierarchy process,and combined fuzzy and support vector machine methods:a case study from Wolong Giant Panda Natural Reserve,China[J].Bulletin of Engineering Geology and the Environment,2016,75(3):923-944.

    [7]夏輝,殷坤龍,梁鑫,等.基于SVM-ANN模型的滑坡易發(fā)性評價:以三峽庫區(qū)巫山縣為例[J].中國地質(zhì)災(zāi)害與防治學報,2018,29(5):13-19.

    [8]喬德京,王念秦,郭有金,等.加權(quán)確定性系數(shù)模型的滑坡易發(fā)性評價[J].西安科技大學學報,2020,40(2):259-267.

    [9]楊永剛,殷坤龍,趙海燕,等.基于C5.0決策樹-快速聚類模型的萬州區(qū)庫岸段鄉(xiāng)鎮(zhèn)滑坡易發(fā)性區(qū)劃[J].地質(zhì)科技情報,2019,38(6):189-197.

    [10]林報嘉,劉曉東,楊川,等.XGBoost機器學習模型與GIS技術(shù)結(jié)合的公路崩塌災(zāi)害易發(fā)性研究[J].公路,2020,65(7):20-26.

    [11]李長冬.三峽庫區(qū)侏羅系地層滑坡中抗滑樁布設(shè)和結(jié)構(gòu)優(yōu)化研究[M].武漢:中國地質(zhì)大學出版社,2015.

    [12]劉堅,李樹林,陳濤,等.基于優(yōu)化隨機森林模型的滑坡易發(fā)性評價[J].武漢大學學報(信息科學版),2018,43(7):1085-1091.

    [13]雷明.機器學習[M].北京:清華大學出版社,2019.

    [14]VAPNIK V.The nature of statistical learning theory[M].Boston:Springer Science &Bussiness Media,2013.

    [15]彭令,牛瑞卿,趙艷南,等.區(qū)域滑坡災(zāi)害風險評估:以長江三峽庫區(qū)秭歸縣為例[J].吉林大學學報(地球科學版),2013,43(3):891-901.

    (編輯:胡旭東)

    Abstract:Landslide susceptibility evaluation is the basis and key to identify landslide hazard target area.The Jurassic stratum area in Zigui County,the Three Gorges Reservoir area,is a landslide hazard prone area.Eight influencing factors including slope,surface cutting depth,curvature,and distance from water system to road network were selected.The information model(I)was coupled with four common machine learning models,namely artificial neural network(ANN),random forest(RF),XGBboostand support vector machine(SVM)to carry out the evaluation of landslide susceptibility.Through a comparative analysis on the receiver operating characteristic curve,it was found that the AUC valuein the ROC curve of the model coupled with information valuemodel and SVMmodel was the largest(0.848).The zoning map of landslide susceptibility showed that the high-prone areas were concentrated on both sides of the water system and the road network,and the results were basically consistent with the actual landslide distribution.The research results can provide reference for regional landslide susceptibility evaluation and landslide disaster prevention.

    Key words:landslide susceptibility evaluation;landslide susceptibility zoning;information valuemodel;multi-model coupling;Jurassic stratum;Three Gorges Reservoir area

    猜你喜歡
    三峽庫區(qū)
    三峽庫區(qū)萬家壩滑坡變形區(qū)穩(wěn)定性復(fù)核研究
    三峽庫區(qū)藥用資源重慶市重點實驗室
    三峽庫區(qū)滑坡土體抗剪強度參數(shù)相關(guān)性構(gòu)建的最優(yōu)Copula函數(shù)
    三峽庫區(qū)不同稻田分布格局下農(nóng)業(yè)小流域徑流磷排放特征
    昭君今若在,定驚故里殊 三峽庫區(qū)興山縣移民搬遷側(cè)記
    中國三峽(2016年5期)2017-01-15 13:58:50
    三峽庫區(qū)產(chǎn)業(yè)培育及結(jié)構(gòu)調(diào)整的思考
    三峽庫區(qū)紫色土坡耕地氮磷徑流特征研究
    長江三峽庫區(qū)農(nóng)村移民增收的設(shè)想
    學習月刊(2015年14期)2015-07-09 03:37:52
    三峽庫區(qū)秭歸縣郭家壩碼頭滑坡特征及穩(wěn)定性評價
    河南科技(2014年15期)2014-02-27 14:12:46
    三峽庫區(qū)消落帶生態(tài)庫岸整治工程設(shè)計概述
    亚洲欧美一区二区三区黑人| 好男人电影高清在线观看| 男人舔女人的私密视频| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区91| 久9热在线精品视频| 法律面前人人平等表现在哪些方面| 91在线观看av| 老熟妇仑乱视频hdxx| 久久精品aⅴ一区二区三区四区| 国产一区二区在线av高清观看| 国模一区二区三区四区视频 | 免费在线观看成人毛片| 精品国产乱子伦一区二区三区| 1024视频免费在线观看| 成人三级做爰电影| 久久精品亚洲精品国产色婷小说| 神马国产精品三级电影在线观看 | 最新美女视频免费是黄的| 757午夜福利合集在线观看| 午夜福利高清视频| 国产一区二区在线av高清观看| 桃红色精品国产亚洲av| 久久久久久久久免费视频了| 一区二区三区激情视频| 亚洲男人天堂网一区| 少妇裸体淫交视频免费看高清 | 一边摸一边做爽爽视频免费| 午夜免费观看网址| 精品久久久久久久末码| 久久久久国产一级毛片高清牌| 日韩欧美在线乱码| 午夜福利在线在线| or卡值多少钱| 女警被强在线播放| 成人精品一区二区免费| 男人舔奶头视频| 麻豆国产97在线/欧美 | 日本黄色视频三级网站网址| 国产私拍福利视频在线观看| 午夜免费激情av| 日韩欧美在线二视频| 久久精品国产综合久久久| 国产高清视频在线播放一区| 九色成人免费人妻av| 午夜免费激情av| 国产精品 欧美亚洲| 亚洲自偷自拍图片 自拍| 国产精品久久久久久亚洲av鲁大| 嫩草影院精品99| 99在线视频只有这里精品首页| xxx96com| 亚洲精品在线美女| 免费高清视频大片| 一级作爱视频免费观看| 看片在线看免费视频| 亚洲精品一卡2卡三卡4卡5卡| 久久国产乱子伦精品免费另类| 午夜福利免费观看在线| 在线观看免费午夜福利视频| 看片在线看免费视频| 国产高清视频在线播放一区| 啦啦啦免费观看视频1| 国产亚洲精品第一综合不卡| 身体一侧抽搐| 欧美日韩精品网址| 黄色成人免费大全| 欧美性长视频在线观看| 国产亚洲精品第一综合不卡| 亚洲一码二码三码区别大吗| 亚洲欧美精品综合久久99| 国产亚洲精品综合一区在线观看 | 成人国语在线视频| 日韩精品中文字幕看吧| 国产精品香港三级国产av潘金莲| 欧美在线一区亚洲| 日韩欧美在线二视频| 色老头精品视频在线观看| 亚洲国产精品999在线| 一个人免费在线观看的高清视频| 国产精品亚洲av一区麻豆| 国产亚洲精品综合一区在线观看 | 婷婷丁香在线五月| 老司机在亚洲福利影院| 欧美午夜高清在线| 成人手机av| 亚洲天堂国产精品一区在线| 中文亚洲av片在线观看爽| 精品国产美女av久久久久小说| 亚洲黑人精品在线| 中文在线观看免费www的网站 | 久99久视频精品免费| 窝窝影院91人妻| 午夜视频精品福利| 国产精品自产拍在线观看55亚洲| 成熟少妇高潮喷水视频| 一区二区三区激情视频| 久久久国产精品麻豆| 三级毛片av免费| 天堂av国产一区二区熟女人妻 | 十八禁人妻一区二区| www.精华液| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出| av福利片在线| 欧美日韩国产亚洲二区| 久久精品aⅴ一区二区三区四区| 亚洲五月天丁香| 日本三级黄在线观看| 国产精品一区二区精品视频观看| www日本黄色视频网| 亚洲片人在线观看| 亚洲,欧美精品.| 久久久水蜜桃国产精品网| 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 天堂av国产一区二区熟女人妻 | 在线观看午夜福利视频| 怎么达到女性高潮| 此物有八面人人有两片| 亚洲黑人精品在线| 日韩大尺度精品在线看网址| 男人舔奶头视频| 俄罗斯特黄特色一大片| 蜜桃久久精品国产亚洲av| 欧美精品亚洲一区二区| 美女高潮喷水抽搐中文字幕| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区av网在线观看| 级片在线观看| 99精品欧美一区二区三区四区| 免费在线观看完整版高清| 色综合站精品国产| 久久久国产成人精品二区| 免费在线观看视频国产中文字幕亚洲| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av| 国产成人av教育| 午夜福利在线观看吧| 中文字幕精品亚洲无线码一区| 成人国语在线视频| 久久 成人 亚洲| 免费在线观看影片大全网站| 日韩欧美三级三区| 最近视频中文字幕2019在线8| 亚洲成人久久爱视频| 怎么达到女性高潮| 少妇粗大呻吟视频| 九色成人免费人妻av| 欧美色视频一区免费| 级片在线观看| 婷婷精品国产亚洲av| 亚洲乱码一区二区免费版| 久久久久久国产a免费观看| 欧美人与性动交α欧美精品济南到| 国产人伦9x9x在线观看| 在线观看免费视频日本深夜| 国产精品久久久人人做人人爽| 岛国在线免费视频观看| 亚洲欧美激情综合另类| 看免费av毛片| 亚洲精品在线美女| 亚洲国产看品久久| 日韩欧美免费精品| 久久久久性生活片| 日日夜夜操网爽| 中文亚洲av片在线观看爽| 在线观看一区二区三区| 嫩草影院精品99| 亚洲五月婷婷丁香| 国产精品 欧美亚洲| 国产熟女xx| 亚洲成人免费电影在线观看| 美女大奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 高清毛片免费观看视频网站| 欧美另类亚洲清纯唯美| 国产伦一二天堂av在线观看| 国产高清视频在线播放一区| 香蕉av资源在线| 欧美黑人巨大hd| 又爽又黄无遮挡网站| 国产亚洲精品综合一区在线观看 | 日韩欧美三级三区| 制服丝袜大香蕉在线| 国产成人影院久久av| 首页视频小说图片口味搜索| 高潮久久久久久久久久久不卡| 亚洲va日本ⅴa欧美va伊人久久| 三级男女做爰猛烈吃奶摸视频| 91字幕亚洲| 午夜a级毛片| 18禁美女被吸乳视频| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 老司机在亚洲福利影院| 免费看a级黄色片| 国产免费av片在线观看野外av| 香蕉丝袜av| 成人18禁在线播放| 亚洲自偷自拍图片 自拍| 国产真人三级小视频在线观看| 亚洲国产欧美一区二区综合| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 好男人在线观看高清免费视频| 搡老熟女国产l中国老女人| 日本五十路高清| 久久久久久国产a免费观看| 亚洲成av人片免费观看| 国产精品亚洲av一区麻豆| 日日摸夜夜添夜夜添小说| 小说图片视频综合网站| 亚洲美女黄片视频| 午夜免费观看网址| 黄色视频,在线免费观看| 真人一进一出gif抽搐免费| 级片在线观看| 久久香蕉精品热| 亚洲美女黄片视频| 久久精品国产清高在天天线| 俄罗斯特黄特色一大片| a级毛片在线看网站| 亚洲五月婷婷丁香| 国产日本99.免费观看| 亚洲欧美精品综合久久99| 两个人看的免费小视频| 不卡av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 中文字幕久久专区| 99久久精品热视频| 老司机午夜福利在线观看视频| 亚洲电影在线观看av| 香蕉丝袜av| 精品第一国产精品| 国产精品久久久av美女十八| 两性夫妻黄色片| 国产成+人综合+亚洲专区| 宅男免费午夜| 18禁黄网站禁片午夜丰满| 亚洲国产精品久久男人天堂| 国产在线精品亚洲第一网站| 法律面前人人平等表现在哪些方面| www国产在线视频色| 亚洲精华国产精华精| 日韩av在线大香蕉| 亚洲一码二码三码区别大吗| 亚洲av第一区精品v没综合| 午夜福利在线观看吧| 国产成人aa在线观看| 欧美大码av| 欧美绝顶高潮抽搐喷水| 日本五十路高清| 在线观看午夜福利视频| 亚洲国产精品合色在线| 国产精品爽爽va在线观看网站| www.精华液| 国产精品久久久久久亚洲av鲁大| 丝袜人妻中文字幕| 岛国在线观看网站| 亚洲熟妇熟女久久| 久久九九热精品免费| 日韩av在线大香蕉| 国产成人av教育| 久久精品亚洲精品国产色婷小说| 老汉色∧v一级毛片| 男女做爰动态图高潮gif福利片| 制服丝袜大香蕉在线| 熟女少妇亚洲综合色aaa.| 男人舔奶头视频| 高清毛片免费观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩中文字幕国产精品一区二区三区| 久久久精品大字幕| 99久久精品热视频| 色av中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 午夜精品在线福利| 久久香蕉精品热| 18禁美女被吸乳视频| 丁香六月欧美| 88av欧美| 国产精品永久免费网站| av在线播放免费不卡| 久久久国产成人免费| 国产精品99久久99久久久不卡| 成在线人永久免费视频| 精品久久久久久久人妻蜜臀av| 久久婷婷人人爽人人干人人爱| 亚洲成人久久爱视频| 国产精品1区2区在线观看.| 亚洲人成77777在线视频| 成人18禁高潮啪啪吃奶动态图| 男女午夜视频在线观看| 久99久视频精品免费| 国产视频一区二区在线看| 亚洲精品中文字幕一二三四区| 一级毛片精品| 久久精品aⅴ一区二区三区四区| 中文字幕精品亚洲无线码一区| 久久香蕉精品热| 免费看美女性在线毛片视频| 一级毛片高清免费大全| 国产午夜精品久久久久久| 国产又黄又爽又无遮挡在线| 午夜免费观看网址| 亚洲男人的天堂狠狠| 在线视频色国产色| 狂野欧美白嫩少妇大欣赏| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清 | 日韩 欧美 亚洲 中文字幕| 99久久国产精品久久久| 亚洲精品中文字幕在线视频| 十八禁人妻一区二区| 欧美精品啪啪一区二区三区| 久久人妻福利社区极品人妻图片| 男男h啪啪无遮挡| 日韩精品免费视频一区二区三区| 国产精品亚洲av一区麻豆| 国产精品综合久久久久久久免费| 日韩av在线大香蕉| 欧美性长视频在线观看| 丁香欧美五月| 欧美成人性av电影在线观看| 看片在线看免费视频| 亚洲精品色激情综合| 午夜免费激情av| 人人妻人人看人人澡| 国产精品免费一区二区三区在线| 久久精品国产亚洲av高清一级| videosex国产| 在线视频色国产色| 一二三四社区在线视频社区8| 久久这里只有精品19| 91在线观看av| 在线视频色国产色| 黄频高清免费视频| 熟女少妇亚洲综合色aaa.| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 免费人成视频x8x8入口观看| 欧美3d第一页| 精品午夜福利视频在线观看一区| 亚洲国产日韩欧美精品在线观看 | 亚洲精华国产精华精| e午夜精品久久久久久久| 久久久精品大字幕| 精品国产超薄肉色丝袜足j| 草草在线视频免费看| 色播亚洲综合网| 又爽又黄无遮挡网站| 国产伦在线观看视频一区| 国产av在哪里看| 叶爱在线成人免费视频播放| xxx96com| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 欧美日韩亚洲综合一区二区三区_| 村上凉子中文字幕在线| 国内毛片毛片毛片毛片毛片| 中亚洲国语对白在线视频| 热99re8久久精品国产| 日日摸夜夜添夜夜添小说| 亚洲av成人一区二区三| 亚洲国产高清在线一区二区三| 午夜免费观看网址| 国产亚洲av嫩草精品影院| 国产97色在线日韩免费| 国产成人av激情在线播放| 色老头精品视频在线观看| 欧美激情久久久久久爽电影| 好看av亚洲va欧美ⅴa在| 亚洲国产精品合色在线| 国产亚洲精品av在线| 999久久久精品免费观看国产| 免费看a级黄色片| 亚洲黑人精品在线| 欧美成人免费av一区二区三区| 在线a可以看的网站| 少妇人妻一区二区三区视频| 久久久国产精品麻豆| 精品熟女少妇八av免费久了| 日本在线视频免费播放| 亚洲电影在线观看av| 亚洲片人在线观看| 男女床上黄色一级片免费看| 91国产中文字幕| 日本一本二区三区精品| 欧美精品啪啪一区二区三区| 男插女下体视频免费在线播放| 国产精品久久视频播放| 国产精品自产拍在线观看55亚洲| 国产单亲对白刺激| 亚洲一区中文字幕在线| 男女午夜视频在线观看| 99精品久久久久人妻精品| 久久久久久久久免费视频了| 在线视频色国产色| 非洲黑人性xxxx精品又粗又长| 无人区码免费观看不卡| 午夜免费激情av| 亚洲国产精品久久男人天堂| 国产精品亚洲一级av第二区| 久久精品综合一区二区三区| 在线播放国产精品三级| 国内精品一区二区在线观看| 精品久久久久久久人妻蜜臀av| 午夜两性在线视频| 草草在线视频免费看| 久久精品成人免费网站| xxx96com| 亚洲真实伦在线观看| cao死你这个sao货| 欧美成狂野欧美在线观看| 欧美一级a爱片免费观看看 | 在线永久观看黄色视频| 午夜精品久久久久久毛片777| 特级一级黄色大片| 亚洲美女视频黄频| 国产成人精品无人区| 国产高清视频在线观看网站| 成人手机av| 成人精品一区二区免费| 日本 av在线| 精品免费久久久久久久清纯| 亚洲成a人片在线一区二区| 欧美日韩瑟瑟在线播放| 丰满的人妻完整版| 久9热在线精品视频| 无限看片的www在线观看| 国产野战对白在线观看| 精品久久久久久久末码| 曰老女人黄片| 国产精品久久久久久精品电影| 亚洲成a人片在线一区二区| 国产亚洲精品av在线| 日本五十路高清| av视频在线观看入口| 国产v大片淫在线免费观看| 此物有八面人人有两片| 欧美人与性动交α欧美精品济南到| 国产精品98久久久久久宅男小说| 日本一本二区三区精品| 美女大奶头视频| 国产真实乱freesex| 99久久久亚洲精品蜜臀av| a级毛片在线看网站| www日本在线高清视频| 少妇被粗大的猛进出69影院| 国产aⅴ精品一区二区三区波| 成人国产一区最新在线观看| 欧美成狂野欧美在线观看| 两个人的视频大全免费| 国产成人影院久久av| 欧美性长视频在线观看| 成在线人永久免费视频| 一区二区三区国产精品乱码| 婷婷精品国产亚洲av| 男女视频在线观看网站免费 | av中文乱码字幕在线| 91麻豆av在线| 亚洲欧美日韩无卡精品| 午夜福利成人在线免费观看| 国产真人三级小视频在线观看| 欧美黄色淫秽网站| 久久久久亚洲av毛片大全| 国产精华一区二区三区| 悠悠久久av| 嫁个100分男人电影在线观看| 亚洲av五月六月丁香网| 国产黄色小视频在线观看| 中文资源天堂在线| 极品教师在线免费播放| 一级毛片高清免费大全| 欧美日韩中文字幕国产精品一区二区三区| 欧美大码av| 一进一出好大好爽视频| 亚洲精品av麻豆狂野| 亚洲精品美女久久av网站| 国产精品野战在线观看| 在线观看一区二区三区| 亚洲av电影不卡..在线观看| 久久久久久免费高清国产稀缺| 精品久久久久久久末码| 午夜成年电影在线免费观看| 国产不卡一卡二| 亚洲一区高清亚洲精品| 成人三级黄色视频| 国产精品一区二区免费欧美| 国产亚洲精品久久久久久毛片| bbb黄色大片| 天天躁夜夜躁狠狠躁躁| 久久午夜亚洲精品久久| 日韩成人在线观看一区二区三区| 精品一区二区三区四区五区乱码| 国产不卡一卡二| 欧美国产日韩亚洲一区| 久久天堂一区二区三区四区| 高清在线国产一区| av在线播放免费不卡| 美女扒开内裤让男人捅视频| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 91av网站免费观看| 欧美一级a爱片免费观看看 | 国产精品永久免费网站| 亚洲精品美女久久久久99蜜臀| 久久 成人 亚洲| 两个人视频免费观看高清| 国产亚洲精品综合一区在线观看 | 色老头精品视频在线观看| 51午夜福利影视在线观看| 国产伦一二天堂av在线观看| 90打野战视频偷拍视频| 久久国产精品人妻蜜桃| 精品国产美女av久久久久小说| 国产成人啪精品午夜网站| 欧美成人午夜精品| 国产精品久久视频播放| 天天躁夜夜躁狠狠躁躁| 成人av在线播放网站| 老熟妇仑乱视频hdxx| 欧美日韩一级在线毛片| 国产99白浆流出| 给我免费播放毛片高清在线观看| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡免费网站照片| 亚洲一区二区三区不卡视频| 国产欧美日韩一区二区三| 久久天堂一区二区三区四区| 国产精品久久久久久人妻精品电影| 国内精品久久久久精免费| 美女午夜性视频免费| 日韩有码中文字幕| 啦啦啦观看免费观看视频高清| 日韩精品青青久久久久久| 久久香蕉国产精品| 国产久久久一区二区三区| 国产精品免费视频内射| 亚洲av电影不卡..在线观看| 国产又色又爽无遮挡免费看| 欧美黄色淫秽网站| 精品人妻1区二区| 一本精品99久久精品77| 色噜噜av男人的天堂激情| 国产激情欧美一区二区| 国产欧美日韩一区二区三| 精品久久久久久久久久久久久| 草草在线视频免费看| 亚洲天堂国产精品一区在线| 欧美黑人欧美精品刺激| 精品国产乱子伦一区二区三区| АⅤ资源中文在线天堂| 看黄色毛片网站| 国产三级在线视频| 国产69精品久久久久777片 | www.熟女人妻精品国产| 国产精品久久久久久久电影 | 亚洲av日韩精品久久久久久密| 久久草成人影院| 欧美一级a爱片免费观看看 | 精品国产亚洲在线| 男人的好看免费观看在线视频 | 亚洲五月婷婷丁香| 国产精品久久久久久久电影 | 婷婷精品国产亚洲av在线| 亚洲成人国产一区在线观看| 欧美大码av| 在线观看舔阴道视频| 久久精品综合一区二区三区| 亚洲国产日韩欧美精品在线观看 | 全区人妻精品视频| 亚洲精华国产精华精| 女人爽到高潮嗷嗷叫在线视频| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 欧美在线一区亚洲| 曰老女人黄片| 亚洲成人中文字幕在线播放| 午夜福利成人在线免费观看| 国产激情欧美一区二区| 麻豆av在线久日| 亚洲 欧美 日韩 在线 免费| 在线观看一区二区三区| 制服人妻中文乱码| 婷婷精品国产亚洲av在线| 国产97色在线日韩免费| 久久精品影院6| 日本在线视频免费播放| 精品久久久久久,| 老司机午夜福利在线观看视频| www日本黄色视频网| 国产精品久久电影中文字幕| 国产在线观看jvid| 午夜福利成人在线免费观看| 国产精品免费一区二区三区在线| 亚洲最大成人中文| 亚洲国产精品成人综合色| 日本五十路高清| 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 久久久久久亚洲精品国产蜜桃av| 国产野战对白在线观看| 欧美中文日本在线观看视频| 床上黄色一级片| netflix在线观看网站| 十八禁人妻一区二区| 成人永久免费在线观看视频| 欧美日韩一级在线毛片| 日韩成人在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 成人午夜高清在线视频|