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

    基于I法與CF法的榕江縣地質(zhì)災(zāi)害易發(fā)性評價

    2022-10-13 11:03:28李微洪托徐世光陳明超秦煜杰曾磊
    關(guān)鍵詞:榕江縣確定性信息量

    李微,洪托,徐世光,陳明超,秦煜杰,曾磊

    (1.昆明理工大學(xué),昆明 650093;2.云南地礦工程勘察集團公司,昆明 650041;3.貴州省地質(zhì)環(huán)境監(jiān)測院,貴陽 550081;4.重慶市生態(tài)環(huán)境科學(xué)研究院,重慶 401147)

    1 引言

    地質(zhì)災(zāi)害是人類生存發(fā)展、工程建設(shè)中的重要問題[1]。當(dāng)前地質(zhì)災(zāi)害防治普遍實行治理為主,預(yù)防為輔,但地質(zhì)災(zāi)害易發(fā)性評價能從空間上給出地質(zhì)災(zāi)害的分布及發(fā)生的概率,如此大大提高了地質(zhì)災(zāi)害預(yù)防的可行性,降低了預(yù)警的高成本,將更有限的財力物力放在危險性較高影響較大的地質(zhì)災(zāi)害上,實現(xiàn)有的放矢的地質(zhì)災(zāi)害防治,對國土開發(fā)具重要的作用[2]。

    目前常用的地質(zhì)災(zāi)害評價方法包括定性評價方法(滑坡編錄方法和知識驅(qū)動方法)和定量評價方法(數(shù)據(jù)驅(qū)動方法和物理驅(qū)動方法)[3]。定量評價方法又可分為:頻率比[4]、邏輯回歸[5]、信息量法[[6]、層次分析[7]、確定性系數(shù)[8]、人工神經(jīng)網(wǎng)絡(luò)[9]和支持向量機[6]等模型。

    本文采用的方法即就是定量評價方法下的信息量模型及確定性系數(shù)模型進行分別評價,再檢驗對比兩種評價模型的評價結(jié)果,選用精度較大的結(jié)果。

    2 研究區(qū)概況

    2.1 研究區(qū)地質(zhì)背景

    本次研究區(qū)為貴州省黔東南州榕江縣,總面積約為3 315.8 km2,位于東經(jīng)108°04′~108°44′,北緯25°26′~26°28′之間;榕江縣距黔東南首府凱里市173 km、省府貴陽市367 km;榕江縣地處貴州高原向廣西丘陵過渡的邊緣地帶,主要為中低山、丘陵、盆地地貌;年代地層為寒武系、震旦系、二疊系、新近系及第四系;地質(zhì)構(gòu)造由水尾、色同向斜和雞冠嶺背斜三大褶皺帶和24條斷層組成;境內(nèi)有河流67條,總長1 303 km;地形西北高東南低,最高海拔1 704 m,最低海拔213 m。氣候?qū)賮啛釒е衼啔夂騾^(qū)[10]。

    2.2 地質(zhì)災(zāi)害類型及發(fā)育規(guī)律

    根據(jù)野外調(diào)查資料顯示,榕江縣現(xiàn)存隱患點特大型不穩(wěn)定斜坡1個,大型5個,包含3個滑坡、1個泥石流、1個不穩(wěn)定斜坡;中型地質(zhì)災(zāi)害24個,為滑坡、崩塌、不穩(wěn)定斜坡,小型地質(zhì)災(zāi)害若干。研究區(qū)現(xiàn)狀地質(zhì)災(zāi)害分布圖如圖1所示。通過對這些災(zāi)點進行詳細分析,可歸納出該區(qū)地質(zhì)災(zāi)害發(fā)生的大致規(guī)律。

    (1) 大型滑坡、崩塌及泥石流等主要發(fā)生在高程及地形起伏度較大的地區(qū)。

    (2) 滑坡及泥石流集中發(fā)生時間段往往為4~10月雨季,受雨水沖刷及地下水作用控制。

    (3) 不穩(wěn)定斜坡隱患點形成往往與人為工程活動相關(guān)聯(lián),往往在進行道路施工或房建施工時形成,若不及時治理,將對工程造成影響。

    (4) 研究區(qū)地層巖性出露以二疊系灰?guī)r為主,碳酸鹽巖類巖溶水也是造成滑坡發(fā)生的重要原因。

    1.第四系;2.新近系;3.二疊系;4.寒武系;5.震旦系;6.道路;7.榕江地質(zhì)構(gòu)造;8.水系;9.災(zāi)害點圖1 現(xiàn)狀地質(zhì)災(zāi)害分布圖

    3 模型概述

    3.1 數(shù)據(jù)來源

    本文進行易發(fā)性評價的數(shù)據(jù)來源主要包括:(1)地理空間數(shù)據(jù)云30 m分辨率DEM數(shù)據(jù);(2)榕江縣1∶50 000地質(zhì)圖;(3)1∶100 000交通圖;(4)野外調(diào)查。

    3.2 評價模型

    結(jié)合實際情況,本文將選用信息量模型及確定性系數(shù)模型進行地質(zhì)災(zāi)害易發(fā)性評價,兩種方法均屬于數(shù)據(jù)型驅(qū)動法,信息量的概念由Shannon在1948年首次提出,1985年左右才被應(yīng)用到地質(zhì)災(zāi)害危險性評價領(lǐng)域,此后發(fā)展迅猛,并深入到理論經(jīng)濟學(xué)、水利工程等多個學(xué)科;確定性系數(shù)最早在1975年由Bchanan等提出,如今已經(jīng)成為地質(zhì)災(zāi)害區(qū)域危險性評價的常用模型[11],確定性系數(shù)模型不僅在本領(lǐng)域應(yīng)用廣泛,而且已深入到計算機科學(xué)與技術(shù)、大氣科學(xué)等多個學(xué)科。

    信息量模型以統(tǒng)計理論為基礎(chǔ),物理意義明確,流程清晰,操作簡單[11];確定性系數(shù)模型可計算災(zāi)害數(shù)量為0的分類,可達到數(shù)據(jù)標(biāo)準(zhǔn)化的效果。

    3.2.1 信息量模型概述

    信息量法(Information Value, I),公式如下[6]:

    (1)

    式中,n為影響因子數(shù)量;I(xi,H)為因子xi影響地質(zhì)災(zāi)害發(fā)生(H)的信息量;xi為評價單元內(nèi)所取的因子等級;Ni為分布在特定評價因子xi內(nèi)的地質(zhì)災(zāi)害面積;N為地質(zhì)災(zāi)害的總面積;Si為含有評價因子xi的面積;S為研究區(qū)總面積[6]。

    信息量總值I越大,說明地質(zhì)災(zāi)害越容易發(fā)生。

    3.2.2 確定性系數(shù)模型概述

    確定性系數(shù)模型(Certainty Factor, CF),其計算公式為:

    (2)

    式中,PPa為地質(zhì)災(zāi)害在評價因子a類中發(fā)生的條件概率;PPs表示地質(zhì)災(zāi)害在整個研究區(qū)中發(fā)生的先驗概率[8];CF的值域為[-1,1];當(dāng)CF為正值時,地質(zhì)災(zāi)害發(fā)生概率增長,越接近1越易于發(fā)生;CF為負值時,地質(zhì)災(zāi)害發(fā)生概率下降,越接近-1越不易于發(fā)生;CF=0時則不確定災(zāi)害是否發(fā)生[12]。

    4 評價因子分析

    榕江縣主要災(zāi)種有滑坡、崩塌、泥石流,隱患災(zāi)點有不穩(wěn)定斜坡。根據(jù)榕江縣地質(zhì)環(huán)境條件、地質(zhì)災(zāi)害成因及分布特征[13],將選取地層巖性、地質(zhì)構(gòu)造、地形起伏、高程、坡度、坡向、距水系距離、距道路距離等8個因子作為本次地質(zhì)災(zāi)害易發(fā)性評價的評價因子,將評價因子用細化等間距的分級方式,結(jié)合區(qū)間內(nèi)的災(zāi)害密度進行分級,分級等級選取不同間隔進行多次嘗試,最終將總體規(guī)律最明顯的等級進行選用。各因子分級如下表1所示。

    高程、坡度、坡向及地形起伏度歸屬于地形地貌類因子。高程因子評價結(jié)果如圖2(a)和表1所示,小于500 m及500~700 m的范圍內(nèi)高程因子的頻率比值大于1,說明在此區(qū)間內(nèi)地質(zhì)災(zāi)害易發(fā)生。坡度因子評價結(jié)果如圖 2(b)和表1所示,坡度20°~40°時,頻率比值均大于1,且頻率比值隨坡度越來越大,說明該區(qū)地質(zhì)災(zāi)害主要發(fā)生在高坡度的邊坡上。坡向因子如圖 2(c)和表1所示,坡向為東南、南向和西南時,頻率比均大于1,表明地質(zhì)災(zāi)害主要發(fā)生在東南-西南方向,位于該方向的邊坡更易發(fā)生地質(zhì)災(zāi)害。由圖2(d)和表1可知,地形起伏度在15~45 m范圍內(nèi)的頻率比值均大于1,此時地質(zhì)災(zāi)害容易發(fā)生。

    地層巖性和地質(zhì)構(gòu)造因子是造成地質(zhì)災(zāi)害發(fā)生的重要原因,它們均屬于災(zāi)害發(fā)生的內(nèi)部條件。榕江縣出露5套地層從圖 2(e)和表1中可知,寒武系、二疊系和第四系的頻率比值均大于1,在這些地層中,地質(zhì)災(zāi)害易于發(fā)生。地質(zhì)構(gòu)造是通過改變坡體結(jié)構(gòu)進而促使災(zāi)害發(fā)生,在此以斷層對災(zāi)害發(fā)生的影響作為地質(zhì)構(gòu)造評價因子,當(dāng)斷層越接近坡體,則坡體結(jié)構(gòu)越破碎,地質(zhì)災(zāi)害發(fā)生頻率則越高。將全縣斷層構(gòu)造線按1 km步長緩沖,根據(jù)圖2(f)和表1可知,距離斷層0~900 m處頻率比值大于1,從整體趨勢而言,距離斷層越遠,頻率比值越小,表明了距離斷層越近,地質(zhì)災(zāi)害越易于發(fā)生。

    地質(zhì)老話說“十水九滑”,說明了水對于地質(zhì)災(zāi)害發(fā)生的必要性。地表水系統(tǒng)的沖刷與切割給地質(zhì)災(zāi)害提供了巨大的誘因[1],故依據(jù)榕江縣地下水系統(tǒng)的規(guī)模、分布狀況,提取流域面積超過10 km2的小流域作為水網(wǎng)線,按200 m步長加以緩沖。據(jù)圖2(g)及表1可知,距離水系600 m內(nèi)的區(qū)域,頻率比值大于1,且0~200 m 范圍內(nèi)頻率比值最大,整體來說,距離水系越遠,頻率比值越小,表明距離水系越近,地質(zhì)災(zāi)害越易于發(fā)生。

    人類工程活動對地質(zhì)災(zāi)害發(fā)生起到觸發(fā)作用,許多人工形成的邊坡往往發(fā)展成為較大的滑坡,威脅工程安全。本文主要考慮人類工程活動的公路交通建設(shè)對地質(zhì)災(zāi)害發(fā)生的影響,所以將據(jù)道路距離作為評價因子進行分析。根據(jù)圖2(h)和表1可知,距離公路距離在400 m內(nèi)時,頻率比值遠大1。

    表1 榕江縣地質(zhì)災(zāi)害評價因子指標(biāo)分級

    5 地質(zhì)災(zāi)害易發(fā)性評價結(jié)果

    5.1 信息量法(I)評價結(jié)果

    本文兩種評價模型均利用ArcGIS軟件來進行評價因子指標(biāo)信息量值的疊加,疊加后的數(shù)據(jù)可生成榕江縣地質(zhì)災(zāi)害信息量值分布圖,利用自然間斷法,將榕江縣的地質(zhì)災(zāi)害易發(fā)性分為極高易發(fā)區(qū)、高易發(fā)區(qū)、中易發(fā)區(qū)和低易發(fā)區(qū)等4個易發(fā)性等級區(qū)間[13、14],經(jīng)綜合分析概化,即得到榕江縣地質(zhì)災(zāi)害易發(fā)性分區(qū)成果圖[15]。結(jié)果如圖3(a)和表1、表2所示。

    圖2 地貌類影響因子分級

    從表2可知,從低易發(fā)區(qū)到極高易發(fā)區(qū),對應(yīng)的頻率比值持續(xù)增大的趨勢,極高易發(fā)區(qū)和高易發(fā)區(qū)頻率比值之和達總頻率比值的89.5%,這說明信息量(I)模型評價榕江縣地質(zhì)災(zāi)害易發(fā)性是有效可行的。由圖3(a)及表2可知,在紅色極高易發(fā)區(qū)地質(zhì)災(zāi)害面積達41.06%,橙色高易發(fā)區(qū)地質(zhì)災(zāi)害面積達27.15%,黃綠色中易發(fā)區(qū)地質(zhì)災(zāi)害面積為26.49%,綠色低易發(fā)區(qū)地質(zhì)災(zāi)害面積僅5.298%。

    5.2 確定性系數(shù)(CF)評價結(jié)果

    榕江縣地質(zhì)災(zāi)害確定性系數(shù)值分布圖的獲取方法和信息量值分布圖的獲取方法相同,將獲取的榕江縣地質(zhì)災(zāi)害確定性系數(shù)值分布圖劃分為極高易發(fā)區(qū)、高易發(fā)區(qū)、中易發(fā)區(qū)和低易發(fā)區(qū)等4個易發(fā)性等級區(qū)間,結(jié)果如圖3(b)和表1、表2所示。

    圖3 地質(zhì)災(zāi)害易發(fā)性分區(qū)

    表2 頻率比值表

    從表2可知,從低易發(fā)區(qū)到極高易發(fā)區(qū),對應(yīng)的頻率比值持續(xù)增大的趨勢,極高易發(fā)區(qū)和高易發(fā)區(qū)頻率比值之和達總頻率比值的89.8%,這說明確定性系數(shù)(CF)模型評價榕江縣地質(zhì)災(zāi)害易發(fā)性也是有效可行的。由圖3(b)及表2可知,在紅色極高易發(fā)區(qū)地質(zhì)災(zāi)害面積達41.72%,橙色高易發(fā)區(qū)地質(zhì)災(zāi)害面積達39.73%,黃綠色中易發(fā)區(qū)地質(zhì)災(zāi)害面積為15.23%,綠色低易發(fā)區(qū)地質(zhì)災(zāi)害面積僅3.31%。

    6 對比分析

    6.1 頻率比檢驗對比

    FR法是建立在地質(zhì)條件相似的地區(qū),孕育地質(zhì)災(zāi)害的概率也相似的假設(shè)之上。FR值可以定量表示環(huán)境因子各屬性區(qū)間對地質(zhì)災(zāi)害發(fā)生的相對影響程度。公式[16]如下:

    (3)

    式中,li是某個環(huán)境因子的第i個屬性區(qū)間內(nèi)的地質(zhì)災(zāi)害面積;L是研究區(qū)內(nèi)地質(zhì)災(zāi)害的總面積,si是環(huán)境因子第i個屬性區(qū)間的面積;S是研究區(qū)總面積[16]。

    Fr大于1表明該環(huán)境因子屬性區(qū)間利于地質(zhì)災(zāi)害發(fā)育,值越大表示對地質(zhì)災(zāi)害發(fā)育的影響也越大;反之,則不利于地質(zhì)災(zāi)害發(fā)育[16]。

    由上表2可看出,兩種評價模型的易發(fā)性由低到高時,頻率比值均逐漸增大,符合評價規(guī)律,當(dāng)頻率比值大于1時,地質(zhì)災(zāi)害發(fā)生概率逐步增大。

    6.2 ROC曲線檢驗對比

    ROC 曲線(receiver operating characteristic curve)是驗證地質(zhì)災(zāi)害易發(fā)性評價精度的手段[17]。ROC 曲線縱軸代表敏感度,代表區(qū)內(nèi)某易發(fā)等級發(fā)生災(zāi)害的累計百分比;其橫軸代表特異度,代表區(qū)內(nèi)各易發(fā)性等級面積累加百分比[12]。

    將ROC曲線下面積AUC作為一個客觀定量的評價指標(biāo)來衡量模型預(yù)測的準(zhǔn)確程度能更加直觀地表現(xiàn)出評價結(jié)果。AUC值介于 0~1 范圍,越接近 1 模型預(yù)測準(zhǔn)確性越高[12]。

    本文2種模型ROC檢驗結(jié)果如圖4所示,CF模型、I模型的AUC值分別為 0.860及0.774。2種模型的AUC值均大于0.75,兩模型均能精確地評價榕江縣地質(zhì)災(zāi)害易發(fā)性,且CF模型較I模型更精確。

    圖4 ROC曲線圖

    7 結(jié)論

    (1) 基于ArcGIS 技術(shù)采用信息量模型(I)、確定性系數(shù)模型(CF)開展地質(zhì)災(zāi)害易發(fā)性評價,I模型和CF模型的ROC精度分別達0.774和0.860,說明該易發(fā)性分區(qū)與研究區(qū)實際情況比較一致。

    (2) 本文用以上兩種模型得出的榕江縣地區(qū)滑坡災(zāi)害易發(fā)性面積大致相同,滑坡易發(fā)性等級的頻率比值由低易發(fā)區(qū)到極高易發(fā)區(qū)均明顯上升,且極高易發(fā)區(qū)與高易發(fā)區(qū)的頻率比值占總體頻率比值均達到了89%,說明兩種評估模型都能夠客觀正確地評估榕江縣地區(qū)地質(zhì)災(zāi)害的易發(fā)性。

    (3) 在進行比選后,將精度較大的CF模型評價結(jié)果作為本次研究成果,在榕江縣,地質(zhì)災(zāi)害極高易發(fā)區(qū)占9.289%,高易發(fā)區(qū)占37.236%,中易發(fā)區(qū)占40.26%,低易發(fā)區(qū)占13.216%。在國土空間利用時應(yīng)該避開地質(zhì)災(zāi)害極高易發(fā)區(qū)及高易發(fā)區(qū),在中易發(fā)區(qū)及低易發(fā)區(qū)應(yīng)當(dāng)綜合考慮建筑所需安全類別來進行合理利用。

    猜你喜歡
    榕江縣確定性信息量
    論中國訓(xùn)詁學(xué)與經(jīng)典闡釋的確定性
    論法律解釋的確定性
    法律方法(2022年1期)2022-07-21 09:18:56
    含混還是明證:梅洛-龐蒂論確定性
    橙香滿園
    盛裝的苗女
    民族文化有傳人
    基于信息理論的交通信息量度量
    榕江冬日臍橙黃
    如何增加地方電視臺時政新聞的信息量
    新聞傳播(2016年11期)2016-07-10 12:04:01
    法律確定性的統(tǒng)合理性根據(jù)與法治實施
    女人精品久久久久毛片| 久久久久久久久久久丰满| 嫩草影院入口| 久久97久久精品| 中文字幕av电影在线播放| 黄色毛片三级朝国网站 | 九九爱精品视频在线观看| 久久久久久久大尺度免费视频| 一二三四中文在线观看免费高清| 日韩三级伦理在线观看| 午夜日本视频在线| 欧美激情极品国产一区二区三区 | 十八禁高潮呻吟视频 | 精品酒店卫生间| 久久毛片免费看一区二区三区| 美女视频免费永久观看网站| 免费观看无遮挡的男女| 精品久久国产蜜桃| 26uuu在线亚洲综合色| 免费人妻精品一区二区三区视频| 亚洲第一区二区三区不卡| 久久女婷五月综合色啪小说| 亚洲av成人精品一区久久| 国产日韩欧美视频二区| 啦啦啦视频在线资源免费观看| 日本av免费视频播放| 亚洲成人av在线免费| 国产午夜精品一二区理论片| 亚洲美女搞黄在线观看| a 毛片基地| 啦啦啦中文免费视频观看日本| www.色视频.com| 九九在线视频观看精品| 午夜免费男女啪啪视频观看| 亚洲精品国产av蜜桃| 国产淫片久久久久久久久| 久久午夜福利片| 一级毛片久久久久久久久女| 中文精品一卡2卡3卡4更新| 亚洲欧美成人综合另类久久久| av.在线天堂| 亚洲在久久综合| 国产淫语在线视频| 中文字幕人妻熟人妻熟丝袜美| 国产视频内射| 久久国内精品自在自线图片| 嫩草影院新地址| 丰满乱子伦码专区| 色婷婷av一区二区三区视频| 亚洲精品久久久久久婷婷小说| 欧美最新免费一区二区三区| av视频免费观看在线观看| 国产深夜福利视频在线观看| 亚洲精品第二区| 高清欧美精品videossex| 久久99热这里只频精品6学生| 自拍欧美九色日韩亚洲蝌蚪91 | 在现免费观看毛片| 91久久精品国产一区二区成人| 国产亚洲午夜精品一区二区久久| 视频区图区小说| 日韩一本色道免费dvd| av天堂中文字幕网| 国产伦在线观看视频一区| tube8黄色片| 久久久久久久大尺度免费视频| 少妇人妻精品综合一区二区| 日韩免费高清中文字幕av| 校园人妻丝袜中文字幕| 亚洲精品国产色婷婷电影| 桃花免费在线播放| 午夜av观看不卡| 纵有疾风起免费观看全集完整版| 大香蕉久久网| 人体艺术视频欧美日本| 成人毛片60女人毛片免费| 高清午夜精品一区二区三区| 久久久欧美国产精品| 久久久久久久亚洲中文字幕| 纯流量卡能插随身wifi吗| 国产成人aa在线观看| 丰满少妇做爰视频| 亚洲,一卡二卡三卡| 国产综合精华液| 午夜免费鲁丝| 国产精品久久久久久av不卡| 好男人视频免费观看在线| 王馨瑶露胸无遮挡在线观看| 99九九线精品视频在线观看视频| 一本—道久久a久久精品蜜桃钙片| 日韩av不卡免费在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 国产在视频线精品| 国产亚洲精品久久久com| 精品一区二区三区视频在线| 日韩一区二区三区影片| 免费黄色在线免费观看| 精品久久久精品久久久| 国产精品嫩草影院av在线观看| 亚洲人与动物交配视频| 狂野欧美激情性bbbbbb| 高清在线视频一区二区三区| 免费黄色在线免费观看| 水蜜桃什么品种好| 亚洲图色成人| 免费人成在线观看视频色| 精品视频人人做人人爽| 99热这里只有是精品在线观看| 亚洲欧美精品专区久久| av免费在线看不卡| 国产精品免费大片| 日本爱情动作片www.在线观看| 久久午夜福利片| 久久国产亚洲av麻豆专区| 久久午夜福利片| 亚洲va在线va天堂va国产| 久久国内精品自在自线图片| 免费在线观看成人毛片| 内地一区二区视频在线| 乱码一卡2卡4卡精品| 赤兔流量卡办理| 国产伦在线观看视频一区| 最近中文字幕2019免费版| 久久久久久久久久久久大奶| 国产真实伦视频高清在线观看| 男女边摸边吃奶| 欧美精品一区二区大全| 精品熟女少妇av免费看| 99精国产麻豆久久婷婷| 成年人免费黄色播放视频 | 色94色欧美一区二区| 国产伦在线观看视频一区| 女人久久www免费人成看片| 亚洲欧美一区二区三区黑人 | 午夜福利,免费看| 中文字幕久久专区| 乱人伦中国视频| 亚洲av不卡在线观看| 最近的中文字幕免费完整| 简卡轻食公司| 丰满饥渴人妻一区二区三| 久久99热6这里只有精品| 丝袜在线中文字幕| 丝袜在线中文字幕| 国产av一区二区精品久久| 国产一区二区三区综合在线观看 | 精华霜和精华液先用哪个| av在线观看视频网站免费| 亚洲精品,欧美精品| 亚洲精品久久久久久婷婷小说| 亚洲真实伦在线观看| 深夜a级毛片| 亚洲欧美成人综合另类久久久| 菩萨蛮人人尽说江南好唐韦庄| 国产成人91sexporn| 99热网站在线观看| 国产成人免费观看mmmm| 久久久久久久久大av| 亚洲精品乱久久久久久| 偷拍熟女少妇极品色| 久久青草综合色| 美女内射精品一级片tv| 色视频在线一区二区三区| 青青草视频在线视频观看| 国产av国产精品国产| 免费观看的影片在线观看| 亚洲精品第二区| 中文字幕亚洲精品专区| 一个人免费看片子| 午夜视频国产福利| 国产精品国产av在线观看| 亚洲,欧美,日韩| 老司机影院毛片| 国产熟女欧美一区二区| 22中文网久久字幕| 免费播放大片免费观看视频在线观看| 亚洲精品日韩av片在线观看| 如日韩欧美国产精品一区二区三区 | 欧美日韩一区二区视频在线观看视频在线| 男女免费视频国产| 久久免费观看电影| 国产极品粉嫩免费观看在线 | 亚洲国产精品999| 国产伦理片在线播放av一区| 亚洲av二区三区四区| 亚洲av二区三区四区| 成人亚洲欧美一区二区av| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 亚洲综合色惰| 97在线视频观看| 黄色毛片三级朝国网站 | 老女人水多毛片| 中文字幕人妻熟人妻熟丝袜美| 99热这里只有是精品50| 欧美一级a爱片免费观看看| 国产又色又爽无遮挡免| 看免费成人av毛片| 日韩人妻高清精品专区| av福利片在线| 午夜激情久久久久久久| 成人影院久久| 亚洲高清免费不卡视频| 精品一区二区三区视频在线| 亚洲欧美中文字幕日韩二区| 亚洲国产精品成人久久小说| 国产精品一二三区在线看| 亚洲av综合色区一区| 晚上一个人看的免费电影| 久久精品国产鲁丝片午夜精品| 国产色婷婷99| 亚洲精华国产精华液的使用体验| 国产白丝娇喘喷水9色精品| 日本-黄色视频高清免费观看| 免费观看a级毛片全部| 一级毛片电影观看| av国产久精品久网站免费入址| 欧美xxxx性猛交bbbb| 亚洲不卡免费看| 91精品一卡2卡3卡4卡| 国产高清不卡午夜福利| 中文欧美无线码| 免费观看a级毛片全部| 亚洲欧美精品专区久久| 99久久精品热视频| 大片电影免费在线观看免费| 青青草视频在线视频观看| 午夜91福利影院| 久久久久精品性色| 妹子高潮喷水视频| 一级爰片在线观看| 欧美人与善性xxx| 久久人妻熟女aⅴ| 国产极品天堂在线| 国产国拍精品亚洲av在线观看| 噜噜噜噜噜久久久久久91| 一区二区三区四区激情视频| 一级毛片我不卡| 一本久久精品| 成人二区视频| 欧美成人精品欧美一级黄| 欧美 日韩 精品 国产| 十八禁高潮呻吟视频 | 99热国产这里只有精品6| 啦啦啦视频在线资源免费观看| 下体分泌物呈黄色| 欧美 日韩 精品 国产| 99热这里只有是精品50| 国产精品.久久久| 一区在线观看完整版| 亚洲美女视频黄频| 噜噜噜噜噜久久久久久91| 国产免费一区二区三区四区乱码| 亚洲欧美日韩另类电影网站| 晚上一个人看的免费电影| 中文字幕久久专区| av有码第一页| 日韩一本色道免费dvd| 3wmmmm亚洲av在线观看| 99久久中文字幕三级久久日本| 在线观看国产h片| 欧美日韩国产mv在线观看视频| a级毛色黄片| 午夜影院在线不卡| 亚洲国产日韩一区二区| 午夜福利在线观看免费完整高清在| 男人舔奶头视频| 亚洲欧洲国产日韩| 激情五月婷婷亚洲| 99九九在线精品视频 | 亚洲欧美一区二区三区黑人 | 精品人妻偷拍中文字幕| 午夜av观看不卡| 啦啦啦在线观看免费高清www| 亚洲va在线va天堂va国产| 久久精品国产a三级三级三级| 九草在线视频观看| 亚洲精品日韩在线中文字幕| 亚洲av.av天堂| 男男h啪啪无遮挡| 亚洲精品中文字幕在线视频 | 99视频精品全部免费 在线| 国产一区二区三区av在线| 精品少妇内射三级| 在线亚洲精品国产二区图片欧美 | 夜夜爽夜夜爽视频| 九九在线视频观看精品| 国产精品秋霞免费鲁丝片| 青青草视频在线视频观看| 欧美丝袜亚洲另类| 国产精品人妻久久久久久| 美女脱内裤让男人舔精品视频| 3wmmmm亚洲av在线观看| 国产精品女同一区二区软件| 91精品国产九色| www.av在线官网国产| 只有这里有精品99| √禁漫天堂资源中文www| 久久久久精品性色| 乱码一卡2卡4卡精品| 亚洲经典国产精华液单| av在线观看视频网站免费| 久久久国产一区二区| 成人无遮挡网站| 国产91av在线免费观看| 五月开心婷婷网| 欧美日韩精品成人综合77777| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| 伦精品一区二区三区| 亚洲欧美精品专区久久| 99九九线精品视频在线观看视频| 香蕉精品网在线| 国产在线视频一区二区| 91aial.com中文字幕在线观看| 一级毛片久久久久久久久女| 国产深夜福利视频在线观看| 全区人妻精品视频| 精品一区二区免费观看| 男人和女人高潮做爰伦理| 大又大粗又爽又黄少妇毛片口| 欧美精品一区二区大全| 91久久精品国产一区二区三区| 少妇熟女欧美另类| 免费少妇av软件| 欧美激情极品国产一区二区三区 | 久久精品夜色国产| 国产精品不卡视频一区二区| 纵有疾风起免费观看全集完整版| 99热网站在线观看| 国产成人a∨麻豆精品| 纯流量卡能插随身wifi吗| 伦精品一区二区三区| 日日啪夜夜撸| 欧美少妇被猛烈插入视频| 欧美日韩av久久| 亚洲精品一二三| 免费高清在线观看视频在线观看| 精品午夜福利在线看| 99热全是精品| av视频免费观看在线观看| 纵有疾风起免费观看全集完整版| 少妇的逼好多水| 亚洲国产精品成人久久小说| 亚洲精华国产精华液的使用体验| 又大又黄又爽视频免费| 久久 成人 亚洲| 少妇精品久久久久久久| 日韩av在线免费看完整版不卡| 美女xxoo啪啪120秒动态图| 熟妇人妻不卡中文字幕| 十八禁网站网址无遮挡 | 曰老女人黄片| 亚洲久久久国产精品| 蜜臀久久99精品久久宅男| 一本色道久久久久久精品综合| 国产精品久久久久成人av| 国产深夜福利视频在线观看| 精品国产一区二区三区久久久樱花| 亚洲一级一片aⅴ在线观看| 高清欧美精品videossex| 男人添女人高潮全过程视频| 边亲边吃奶的免费视频| 亚洲欧洲精品一区二区精品久久久 | 丰满少妇做爰视频| 黑人巨大精品欧美一区二区蜜桃 | 精品少妇内射三级| 国产一级毛片在线| 黄色欧美视频在线观看| av在线观看视频网站免费| 18禁裸乳无遮挡动漫免费视频| 久久99精品国语久久久| 少妇的逼好多水| 丰满饥渴人妻一区二区三| 岛国毛片在线播放| 中文精品一卡2卡3卡4更新| 一级黄片播放器| 人妻制服诱惑在线中文字幕| 亚洲综合色惰| 亚洲国产色片| 亚洲无线观看免费| 少妇人妻久久综合中文| 国产美女午夜福利| 伦精品一区二区三区| 日本色播在线视频| av卡一久久| 亚洲三级黄色毛片| 亚洲精品视频女| 综合色丁香网| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 大片电影免费在线观看免费| 在线观看www视频免费| 九九在线视频观看精品| 久久ye,这里只有精品| 99热这里只有精品一区| 国产精品一区二区性色av| 一本色道久久久久久精品综合| 国产美女午夜福利| 高清视频免费观看一区二区| 少妇人妻 视频| 一本一本综合久久| 免费看不卡的av| 亚洲欧美成人综合另类久久久| 男男h啪啪无遮挡| 国产日韩欧美在线精品| 亚洲激情五月婷婷啪啪| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 男女啪啪激烈高潮av片| 国产黄片视频在线免费观看| 极品少妇高潮喷水抽搐| 美女大奶头黄色视频| 在线观看免费高清a一片| 美女中出高潮动态图| 亚洲精品视频女| 国产av国产精品国产| 99热这里只有是精品50| 国产精品偷伦视频观看了| av国产久精品久网站免费入址| 日本91视频免费播放| 国产黄频视频在线观看| 能在线免费看毛片的网站| 少妇的逼水好多| 亚洲综合色惰| 国产成人aa在线观看| 免费看光身美女| 18禁裸乳无遮挡动漫免费视频| 秋霞在线观看毛片| 色网站视频免费| 国产伦精品一区二区三区视频9| 一本一本综合久久| 欧美激情国产日韩精品一区| 日韩大片免费观看网站| 国产成人aa在线观看| 三级经典国产精品| 青青草视频在线视频观看| 日本欧美视频一区| 日韩精品有码人妻一区| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| 最新的欧美精品一区二区| 嫩草影院入口| 国产伦理片在线播放av一区| 丰满人妻一区二区三区视频av| 国产精品久久久久久久久免| 国产成人精品无人区| 好男人视频免费观看在线| 久久99热6这里只有精品| 中文资源天堂在线| 一边亲一边摸免费视频| 亚洲不卡免费看| 自拍偷自拍亚洲精品老妇| 一级毛片我不卡| 色5月婷婷丁香| 国产亚洲91精品色在线| av线在线观看网站| 亚洲精品日韩av片在线观看| 少妇熟女欧美另类| 日本黄大片高清| 国产免费一区二区三区四区乱码| 亚洲精品一二三| 国产精品人妻久久久影院| 久久精品国产亚洲av天美| 美女大奶头黄色视频| 亚洲,一卡二卡三卡| 成年女人在线观看亚洲视频| 国产一区二区在线观看日韩| 99久国产av精品国产电影| 波野结衣二区三区在线| 色视频在线一区二区三区| 欧美精品高潮呻吟av久久| 国产精品成人在线| av女优亚洲男人天堂| 女性生殖器流出的白浆| 国产国拍精品亚洲av在线观看| 亚洲av综合色区一区| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 欧美97在线视频| 少妇被粗大猛烈的视频| 国产黄色视频一区二区在线观看| 成年美女黄网站色视频大全免费 | 只有这里有精品99| 国语对白做爰xxxⅹ性视频网站| 亚洲精华国产精华液的使用体验| 99久久中文字幕三级久久日本| 国产乱来视频区| 精品熟女少妇av免费看| av在线播放精品| 国产在线男女| 一区二区av电影网| 国产69精品久久久久777片| 永久免费av网站大全| 免费观看无遮挡的男女| 久久这里有精品视频免费| 国产一区二区在线观看av| 色视频www国产| 亚洲欧美精品自产自拍| 日本爱情动作片www.在线观看| 精品一区二区三卡| 精品国产乱码久久久久久小说| 黄色视频在线播放观看不卡| 亚洲经典国产精华液单| 亚洲欧洲精品一区二区精品久久久 | 日韩 亚洲 欧美在线| 麻豆成人av视频| 婷婷色综合大香蕉| 久久人人爽人人片av| 免费观看在线日韩| 亚洲成人手机| 国产精品不卡视频一区二区| 一级片'在线观看视频| 亚洲图色成人| 亚洲天堂av无毛| 国产日韩欧美视频二区| 纵有疾风起免费观看全集完整版| 六月丁香七月| 国产成人freesex在线| 国产亚洲最大av| 久久婷婷青草| 女人精品久久久久毛片| 国产一级毛片在线| 最新的欧美精品一区二区| 国产精品久久久久成人av| 性高湖久久久久久久久免费观看| 丰满乱子伦码专区| 一级爰片在线观看| 亚洲av成人精品一二三区| 性色av一级| freevideosex欧美| 精品人妻熟女av久视频| .国产精品久久| 美女福利国产在线| 欧美激情国产日韩精品一区| 91久久精品国产一区二区三区| 亚洲三级黄色毛片| 亚洲精品乱久久久久久| 蜜桃久久精品国产亚洲av| 自拍偷自拍亚洲精品老妇| 成年人午夜在线观看视频| 国产精品一区二区在线观看99| 嫩草影院新地址| 精品人妻熟女毛片av久久网站| 国产亚洲最大av| 美女中出高潮动态图| 亚洲精品乱久久久久久| videos熟女内射| 精品国产乱码久久久久久小说| 性色av一级| 国产黄频视频在线观看| 午夜激情福利司机影院| av网站免费在线观看视频| 丰满少妇做爰视频| 精品国产露脸久久av麻豆| 亚洲色图综合在线观看| 国产探花极品一区二区| 亚洲精品国产色婷婷电影| 美女中出高潮动态图| 99久久中文字幕三级久久日本| 97在线人人人人妻| 十分钟在线观看高清视频www | 两个人免费观看高清视频 | 国产乱人偷精品视频| 亚洲综合色惰| 亚洲精品国产av蜜桃| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区国产| 久久人妻熟女aⅴ| 日韩欧美精品免费久久| 九草在线视频观看| 亚洲欧美日韩卡通动漫| 国产成人freesex在线| 日日摸夜夜添夜夜爱| 国产成人freesex在线| 久久6这里有精品| 久久鲁丝午夜福利片| 久久青草综合色| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 伊人久久精品亚洲午夜| 中文字幕av电影在线播放| 女性生殖器流出的白浆| 亚洲,一卡二卡三卡| 黄色配什么色好看| 日本av手机在线免费观看| 亚洲性久久影院| 成人漫画全彩无遮挡| 在线精品无人区一区二区三| 国产精品女同一区二区软件| 色5月婷婷丁香| 91在线精品国自产拍蜜月| 少妇高潮的动态图| 国产免费一区二区三区四区乱码| 在线天堂最新版资源| 精品久久久精品久久久| 天堂中文最新版在线下载| 国产精品熟女久久久久浪| 中文天堂在线官网| 韩国av在线不卡| 3wmmmm亚洲av在线观看| 午夜精品国产一区二区电影| 少妇被粗大猛烈的视频| 老熟女久久久| 日本vs欧美在线观看视频 | 国产精品成人在线| 亚洲av综合色区一区| 美女主播在线视频| 中文字幕制服av| 国产极品粉嫩免费观看在线 | 国产亚洲最大av| 亚洲精品一二三| 国产综合精华液|