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

    變量取值方式對崩塌易發(fā)性的影響

    2022-01-17 06:37:06陳麗霞
    地理空間信息 2021年12期
    關鍵詞:源區(qū)信息量易發(fā)

    趙 玉,陳麗霞*,付 圣

    (1.中國地質(zhì)大學地球物理與空間信息學院,湖北 武漢 430074)

    高陡斜坡上的巖土體突然脫離母巖并堆積到坡腳的崩塌現(xiàn)象,極易造成人員傷亡與財產(chǎn)損失,尤其是交通沿線的高陡邊坡問題最為突出[1]。在防控崩塌時,對崩塌源區(qū)的提前判識是首要任務,其結果的準確性直接影響落石的運動軌跡預測和風險的量化判定[2]。

    識別崩塌源區(qū)的常用方法有野外地質(zhì)判斷、統(tǒng)計分析或智能模型分析等[3-8]。其中,邏輯回歸模型(LRM)是一種廣泛用于斜坡地質(zhì)災害易發(fā)性分析或空間預測的統(tǒng)計分析方法[3-7]。該模型能揭示因變量和多個互不相關的自變量之間的多元回歸關系[5],LRM 的建模過程本身具有挑選變量的優(yōu)點[6]。隨著應用需求的不斷提高,目前LRM 已從二元發(fā)展到多元,再到與其他模型聯(lián)合使用的階段[7]。

    在預測模型的使用中,無論何種方式,都面臨輸入變量取值方式的問題。對于崩塌,其分析因素包括地形地貌、地層巖性、地質(zhì)構造、人類工程活動等。這些數(shù)據(jù)以兩類方式存在:連續(xù)型數(shù)據(jù)(如坡度值)和分類型數(shù)據(jù)(如巖性類型)。連續(xù)變量的取值方式有變量的原值、變量的歸一化值和變量的標準化值。其中,變量的歸一化值是連續(xù)變量原值與評價單元內(nèi)的最大和最小值的差值之比獲得[8];變量的標準化值是將變量值標準化到-1 到+1 之間[9]。分類變量的取值類型有屬性,分級值、模型處理分級值、數(shù)學公式處理分級值和權重分級值,包括根據(jù)因子本身屬性的分級或者使用某種模型獲得的因子對災害的權重或貢獻程度分級。因子屬性分級值(如坡度因子),以10°為間隔從小到大分為1 到9 的序號值[10-11],結合模型處理后的分級值,例如分級后的因子經(jīng)過信息量模型處理后獲得的信息量值[11];結合數(shù)學公式處理后的分級值,例如結合加權頻率比模型[12]進行歸一化的分級值;權重分級值,如根據(jù)專家經(jīng)驗或數(shù)學模型對不同因子賦權重,也可以對不同因子的不同級別賦權重,最終獲得分級值[13]。可見,變量的取值方式多樣,但其對災害評價結果是否有影響以及影響程度如何,需要進行對比驗證。

    鑒于此,本文以湖北省利川龍駒壩G318 沿線崩塌為例,使用評價因子的3 種取值方式(原始值、信息量值和信息量排序值),采用LRM 評價崩塌源區(qū)識別的差異性,最終選用最優(yōu)取值方案實現(xiàn)崩塌源的識別。

    1 研究方法

    在資料收集、遙感解譯和野外調(diào)查的基礎上,利用統(tǒng)計分析方法從地質(zhì)數(shù)據(jù)中提取崩塌的主要影響因子,深入研究崩塌的發(fā)育規(guī)律。因子包括坡度、坡向、曲率、粗糙度、巖性、斜坡結構、節(jié)理密度、水系線密度、歸一化植被指數(shù)(NDVI)、土地利用類型。通過統(tǒng)計崩塌集中分布的坡度,初步篩選出崩塌源范圍。之后對崩塌形成內(nèi)在條件和外在誘因上進行變量提取,采用3 種變量輸入方式(原始值、信息量值和信息量排序值),在ArcGIS 內(nèi)實施崩塌源區(qū)易發(fā)性評價,用5 類易發(fā)性等級(極高、高、中、低和極低)表征崩塌源區(qū)發(fā)生災害的可能性大小。最后分別從易發(fā)區(qū)分布空間特征、極高和高易發(fā)區(qū)的面積比例和評價精度三方面,判斷3 種結果的崩塌源識別能力,技術路線見圖1。

    圖1 技術路線圖

    3 種變量取值方案為:

    方案(1)通過ArcGIS 提取各個因子的原始值。例如坡度因子以連續(xù)型坡度值賦給研究區(qū)的柵格單元;斜坡結構因子則以分類型數(shù)據(jù)進行賦值。連續(xù)型數(shù)據(jù)和分類型數(shù)據(jù)分別保持原值。

    方案(2)通過信息量模型[10]獲取各個因子的信息量值。信息量值為分類型,各因子以離散形式輸入。

    方案(3)是在(2)的基礎上,對各因子的各級別的信息量值由大到小進行重要性由高到低的排序后獲取其排序值,排序值為分類型,各因子的各級別以離散的編號值輸入。

    1.1 信息量模型原理

    信息量模型最早由殷坤龍[10]提出并用于地質(zhì)災害易發(fā)性預測,其計算公式為:

    式中,I為單元信息量值;Ii為因素xi對地質(zhì)災害所提供的信息量;Si為因素xi所占單元總面積;Si0為因素xi單元中發(fā)生地質(zhì)災害的面積之和;A0為研究區(qū)地質(zhì)災害的單元面積之和;A為研究區(qū)總單元面積之和。

    1.2 LRM 原理

    LRM 是目前廣泛用于地質(zhì)災害易發(fā)性評價的數(shù)學統(tǒng)計模型[3,4,7],LRM 中表示崩塌的因變量Y是一個二分類變量,其取值Y=1 和Y=0,分別代表發(fā)生崩塌和未發(fā)生崩塌。影響Y取值的n個自變量因子分別為X1,X2,…,Xn,在n個自變量作用下崩塌發(fā)生的條件概率為P=P(Y=1|X1,X2,…,Xn),則LRM 可表示為:

    式中,zi為中間變量參數(shù),無實際意義;a0為回歸常數(shù);ai為第i個變量的回歸系數(shù)(i=1,2,…,n);Xij為第i號單元中第j個變量的取值,存在崩塌取1,否則取0;Pi為第i號單元內(nèi)崩塌發(fā)生概率回歸值(i=1,2,…,n)。

    1.3 評價精度分析

    采用ROC 曲線[14]對3 種方案的LRM 結果進行檢驗。曲線下的面積AUC 用來表示預測準確性,AUC越高,預測精度值越高。

    2 研究區(qū)及數(shù)據(jù)

    2.1 研究區(qū)背景

    研究區(qū)位處重慶市萬州區(qū)和湖北省利川市,包含磨刀溪中下游及其支流龍駒河(圖2a),坐標范圍為108°30′E~108°45′E,30°30′N~30°40′N。地處四川盆地東緣,以江津至奉節(jié)沿齊岳山一線為界,平行展布北東向窄條狀中低山脈、地勢相對低凹的寬闊臺地狀山地和平緩丘陵(圖2b)。

    圖2 研究區(qū)地理位置圖和地質(zhì)構造圖

    區(qū)內(nèi)地層除背斜核部出露中生代三疊系外,其余均為侏羅系,主要為侏羅系中統(tǒng)上沙溪廟組(J2s上部以泥巖為主,中部粉砂質(zhì)泥巖,下部厚層長石砂巖)、下沙溪廟組(J2xs泥巖、泥巖夾長石石英砂巖)和新田溝組(J2x頂部砂質(zhì)泥巖,中部夾鈣質(zhì)砂巖),侏羅系下統(tǒng)自流井組(J1~2z頁巖夾石英砂巖)、珍珠沖組(J1z中厚層石英砂巖、頁巖)、以及三疊系上統(tǒng)須家河組(T3xj厚層長石石英砂巖夾頁巖或煤層)和三疊系中統(tǒng)巴東組(T2b中厚~厚層含粉砂質(zhì)灰?guī)r夾鈣質(zhì)頁巖)。

    研究區(qū)崩塌共105 處,面積累計39.45×104 m2,體積累計321.89×104 m3。受地形地貌和巖土工程性質(zhì)的影響,崩塌在上沙溪廟組密集發(fā)育,在其他地層相對較弱(圖3a)。這種近水平地層內(nèi)軟硬巖互層的組合極易產(chǎn)生崩塌,嚴重威脅道路G318 和x553 的安全(圖3b,3c)。

    圖3 崩塌分布圖及典型崩塌

    2.2 研究區(qū)數(shù)據(jù)源

    基礎數(shù)據(jù):①Pleiades 影像數(shù)據(jù)(分辨率為0.5 m,影像時間:2 014.9.22)和高分一號數(shù)據(jù)(分辨率為1 m,影像時間:2015-03-30);②地質(zhì)圖(1∶10 000);③野外踏勘崩塌數(shù)據(jù)。

    2.3 崩塌源區(qū)初步篩選

    崩塌特征受高陡地形控制,為了鎖定崩塌源區(qū),利用最小坡度對區(qū)內(nèi)柵格進行初步篩選。采用5°間隔將坡度分成16 類,通過統(tǒng)計各級內(nèi)崩塌的比例,確定出其分布的最小坡度值。崩塌位于10°~70°內(nèi);在坡度達到25°時,崩塌數(shù)量陡增; 在40°時達到峰值(圖4)。因此,選取25°為最小坡度。

    圖4 崩塌源區(qū)與坡度范圍對應圖

    2.4 因子預處理

    通過對歷史崩塌的發(fā)育特點與規(guī)律分析,選擇10 個因子進行評價,各因子的取值情況見表1。

    表1 崩塌源區(qū)評價變量分級表

    1)地形地貌因子。坡度是控制崩塌重要因子,在較陡的斜坡帶,應力集中,結構面發(fā)育時易崩塌,且?guī)r石風化強度大,其屬于連續(xù)型,通過等間距劃分(圖5a)。坡向控制不同坡面所受到的光照強度,其為連續(xù)型,通過等間距劃分(圖5b)。曲率反映斜坡坡面形態(tài),其在空間幾何上是斜坡表面上某個點的切線方向角所對弧長的轉(zhuǎn)動率(圖5c)。坡面形態(tài)反映斜坡的應力狀態(tài)以及對地表水的匯聚能力,粗糙度反映坡體表面的起伏變化與侵蝕程度,定義為曲面面積與投影面積之比[16]。在GIS 中用自然間斷法分級(圖5d)。

    圖5 崩塌源區(qū)評價因子圖

    2)工程地質(zhì)因子。巖性分為三類:軟硬巖互層(J2s,J2xs)、軟巖夾硬巖(J1~2z,J2x)和硬巖含軟弱夾層(T3xj,T2b,J1z),屬于分類型(圖5e)。斜坡結構反映結構面與臨空面的空間關系,表征危巖體掉落概率大小。根據(jù)斜坡劃分標準[17]分為5 類:順向坡[0°~30°)、斜順坡[30°~60°)、橫向坡[60°~120°)、逆斜坡[120°~150°)及逆向坡[150°~180°],其屬于分類型(圖5f),節(jié)理是崩塌重要控制因素,通過野外結構面測量獲得,其為連續(xù)型,通過自然間斷法分級(圖5g)。

    3)其他誘發(fā)因子。沖溝密度體現(xiàn)臨空條件或反映水的侵蝕能力,屬于連續(xù)型(圖5h)。NDVI 反映地表植被的覆蓋狀況,良好的坡體穩(wěn)定性更好。通過近紅外和紅外波段計算得到,屬于連續(xù)型,通過等間距劃分(圖5i)。土地利用類型反映人類活動對斜坡的擾動情況。利用高分一號數(shù)據(jù)并采用波譜角分類(SAM)方法進行地物識別,得到土地利用類型圖,其為分類型(圖5j)。

    3 結 果

    3.1 崩塌源區(qū)空間分布特征

    3 種變量取值方式在識別崩塌源極高易發(fā)區(qū)的能力上有差異。經(jīng)過與信息量模型耦合后(方案2 中信息量值和方案3 中的信息量排序值)的評價,識別極高易發(fā)等級的崩塌源區(qū)的能力較原值的輸入方式強(表2)。在原值方案下,研究區(qū)只有30%面積被判定為極高易發(fā)區(qū);而在信息量值或者排序值的方案下,極高易發(fā)區(qū)的識別比例上升10%。

    表2 不同變量取值方式下崩塌源易發(fā)區(qū)面積對比

    圖6a、6b、6c 結果表明,變量輸入方式不會影響高易發(fā)區(qū)的分布情況,地質(zhì)條件才是決定崩塌的重要因素。無論是哪一種輸入方式,其崩塌源高易發(fā)區(qū)分布特征相同,即南北高,中部低。實質(zhì)上,受地質(zhì)構造和巖性的控制,主要分布在趕場向斜的北翼和馬頭場向斜的南翼(圖3)的上沙溪廟組地層內(nèi)。

    3.2 崩塌源區(qū)識別精度檢驗與對比

    變量的3 種取值方式在崩塌源識別精度上有差異(圖6d)。通過比較ROC 曲線下AUC 值發(fā)現(xiàn),采用信息量排序值作為變量的輸入值得到的評價結果精度最高,達到0.855,其次是信息量值(0.832)和原值(0.821)的輸入方案。進一步表明,在使用信息量模型與LRM 耦合評價崩塌源區(qū)易發(fā)性時,信息量排序值結果最優(yōu)。從輸入變量的特點分析其原因是:信息量值比原值數(shù)據(jù)更多包含了因子變量的權重值,而排序值相當于對信息量值大的因子賦予了二次權重,LRM 回歸系數(shù)相當于三次權重,多權重的計算促使評價精度最優(yōu)。該研究結果與樊芷吟[10]所得結論一致,即將信息量模型與LRM 耦合時,結果最優(yōu)。

    圖6 研究區(qū)崩塌源區(qū)預測結果圖和評價精度曲線圖

    4 結 語

    針對變量存在多種輸入方式從而影響崩塌源區(qū)識別精度的問題,以龍駒壩境內(nèi)崩塌為例,討論了變量的不同輸入方式對崩塌源區(qū)預測的影響,所得結論如下:

    1)崩塌源的分布受地形和地質(zhì)條件控制,邏輯回歸模型中變量的取值方式不會改變其總體空間分布特點,說明在評價地質(zhì)災害易發(fā)性工作中,災害發(fā)育的地形、地質(zhì)條件比模型重要。

    2)3 種變量取值方式在崩塌源識別精度上有差異,以信息量排序值為取值方式能得到最高評價精度,同時,與信息量模型進行耦合能優(yōu)化邏輯回歸模型的預測結果,這與前人研究結論一致。

    3)邏輯回歸模型只是數(shù)學統(tǒng)計模型中的代表模型,變量取值方式對于其他預測模型是否有同樣的精度提升的效果還需進一步研究。

    4)案例研究區(qū)崩塌源總體處于坡度高于25°的斜坡帶,各方案結果均表明極高和高易發(fā)區(qū)分布于趕場向斜的北翼和馬頭場向斜的南翼的上沙溪廟組地層內(nèi)。

    猜你喜歡
    源區(qū)信息量易發(fā)
    機用鎳鈦銼在乳磨牙根管治療中的應用
    貴州省地質(zhì)災害易發(fā)分區(qū)圖
    大眾科學(2022年5期)2022-05-18 13:24:20
    夏季羊易發(fā)疾病及防治方法
    冬小麥蒸散源區(qū)代表性分析
    冬季雞腸炎易發(fā) 科學防治有方法
    基于信息理論的交通信息量度量
    渭河源區(qū)徑流量變化特征及趨勢分析
    如何增加地方電視臺時政新聞的信息量
    新聞傳播(2016年11期)2016-07-10 12:04:01
    基于多尺度互信息量的數(shù)字視頻幀篡改檢測
    計算機工程(2015年4期)2015-07-05 08:29:20
    基于聯(lián)合熵和交互信息量的視頻篡改檢測
    国产免费av片在线观看野外av| 久久精品国产综合久久久| 极品人妻少妇av视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲美女黄色视频免费看| tocl精华| 在线观看免费午夜福利视频| 制服诱惑二区| 亚洲色图 男人天堂 中文字幕| 制服诱惑二区| 国产极品粉嫩免费观看在线| 久久久国产成人免费| 男女高潮啪啪啪动态图| 亚洲精品第二区| 精品一区二区三卡| 亚洲欧美一区二区三区黑人| 一级片'在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区三区在线| 十八禁网站网址无遮挡| 在线观看免费午夜福利视频| 午夜日韩欧美国产| 韩国精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩一区二区精品| av在线播放精品| 国产精品国产三级国产专区5o| 人妻 亚洲 视频| 男女免费视频国产| 大片电影免费在线观看免费| 777久久人妻少妇嫩草av网站| 午夜影院在线不卡| 男女免费视频国产| 2018国产大陆天天弄谢| 久热爱精品视频在线9| 男人添女人高潮全过程视频| 久久久久久久久久久久大奶| av天堂久久9| 一级黄色大片毛片| 深夜精品福利| 欧美人与性动交α欧美精品济南到| 亚洲情色 制服丝袜| 麻豆乱淫一区二区| 久久久久精品人妻al黑| 久久这里只有精品19| 丁香六月欧美| 久久久国产欧美日韩av| 最近中文字幕2019免费版| 三级毛片av免费| 99国产精品99久久久久| 精品亚洲成国产av| 亚洲一码二码三码区别大吗| 伦理电影免费视频| 国产精品免费视频内射| 国产欧美日韩精品亚洲av| 成年人黄色毛片网站| 亚洲情色 制服丝袜| 妹子高潮喷水视频| 国产亚洲欧美在线一区二区| 亚洲国产av影院在线观看| 美女福利国产在线| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久久久久久大奶| 91老司机精品| 欧美日韩黄片免| 日本精品一区二区三区蜜桃| 国产成人一区二区三区免费视频网站| 国产成人a∨麻豆精品| 色视频在线一区二区三区| 国产精品一区二区精品视频观看| 久久中文看片网| 一区二区三区精品91| 久久国产精品人妻蜜桃| 在线亚洲精品国产二区图片欧美| 精品国产乱码久久久久久男人| 亚洲伊人久久精品综合| 9191精品国产免费久久| 美女中出高潮动态图| 男女下面插进去视频免费观看| 高清av免费在线| a 毛片基地| 久久精品国产亚洲av香蕉五月 | 男人添女人高潮全过程视频| 黄片大片在线免费观看| 91精品伊人久久大香线蕉| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 久久这里只有精品19| 一本综合久久免费| 两个人免费观看高清视频| 18在线观看网站| 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 国产免费现黄频在线看| 成人国语在线视频| 亚洲成人手机| 黄色视频在线播放观看不卡| 咕卡用的链子| 午夜激情av网站| 免费在线观看影片大全网站| 日韩欧美国产一区二区入口| 欧美精品高潮呻吟av久久| 亚洲精品一二三| 亚洲av日韩精品久久久久久密| 亚洲一码二码三码区别大吗| 久9热在线精品视频| av又黄又爽大尺度在线免费看| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三区av网在线观看 | 国产精品一区二区在线不卡| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 一边摸一边抽搐一进一出视频| 久久狼人影院| 欧美成狂野欧美在线观看| 人人妻人人添人人爽欧美一区卜| 国产精品影院久久| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久 | 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 日韩视频在线欧美| 如日韩欧美国产精品一区二区三区| 一区二区三区四区激情视频| 天天影视国产精品| 亚洲天堂av无毛| 99国产精品99久久久久| 精品一品国产午夜福利视频| 99精品久久久久人妻精品| 在线观看舔阴道视频| 天天影视国产精品| 欧美人与性动交α欧美精品济南到| 热re99久久国产66热| 午夜福利视频在线观看免费| 日韩大码丰满熟妇| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三卡| 精品国产乱码久久久久久小说| 91九色精品人成在线观看| 久久国产精品男人的天堂亚洲| 国产黄色免费在线视频| 男男h啪啪无遮挡| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到| 黑人操中国人逼视频| 天天影视国产精品| 99香蕉大伊视频| 天天躁日日躁夜夜躁夜夜| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 中文字幕精品免费在线观看视频| 两个人免费观看高清视频| 亚洲天堂av无毛| 天天躁日日躁夜夜躁夜夜| 久久亚洲精品不卡| 在线观看免费午夜福利视频| 一二三四在线观看免费中文在| www.999成人在线观看| 黄色片一级片一级黄色片| 我要看黄色一级片免费的| 另类亚洲欧美激情| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 日韩欧美一区视频在线观看| 午夜精品久久久久久毛片777| 在线观看舔阴道视频| 久久久精品免费免费高清| 久久九九热精品免费| 在线观看免费高清a一片| 亚洲国产毛片av蜜桃av| 国产熟女午夜一区二区三区| 欧美日韩精品网址| 亚洲av电影在线进入| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 久久国产精品大桥未久av| 国产成人啪精品午夜网站| 国产精品久久久久久精品古装| 国产成人精品在线电影| 正在播放国产对白刺激| 波多野结衣av一区二区av| 亚洲一区二区三区欧美精品| 真人做人爱边吃奶动态| h视频一区二区三区| 两性夫妻黄色片| 亚洲欧美清纯卡通| 精品人妻熟女毛片av久久网站| 亚洲av男天堂| 国产亚洲一区二区精品| 亚洲一码二码三码区别大吗| 一区二区三区四区激情视频| 国产欧美日韩一区二区三 | 国产老妇伦熟女老妇高清| 一二三四在线观看免费中文在| xxxhd国产人妻xxx| 精品一区在线观看国产| 久久精品亚洲熟妇少妇任你| 亚洲av成人不卡在线观看播放网 | 亚洲第一欧美日韩一区二区三区 | 成人国产av品久久久| 欧美中文综合在线视频| 91精品伊人久久大香线蕉| 国产一级毛片在线| 精品少妇黑人巨大在线播放| 亚洲精华国产精华精| 大码成人一级视频| 天天躁日日躁夜夜躁夜夜| 十八禁高潮呻吟视频| 女性生殖器流出的白浆| 亚洲色图 男人天堂 中文字幕| 青青草视频在线视频观看| 午夜激情av网站| 日韩精品免费视频一区二区三区| 国产成人精品久久二区二区91| 国产男人的电影天堂91| 欧美人与性动交α欧美软件| 99re6热这里在线精品视频| 各种免费的搞黄视频| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 亚洲av电影在线观看一区二区三区| 乱人伦中国视频| 中亚洲国语对白在线视频| 99久久综合免费| 两个人免费观看高清视频| 亚洲国产av影院在线观看| 99热全是精品| 亚洲一码二码三码区别大吗| av电影中文网址| 欧美激情 高清一区二区三区| 成人国产av品久久久| 国产淫语在线视频| 日本vs欧美在线观看视频| 亚洲五月色婷婷综合| 国产成人精品在线电影| 看免费av毛片| 免费在线观看完整版高清| 亚洲第一欧美日韩一区二区三区 | 国产日韩欧美在线精品| 亚洲国产av新网站| 国产一级毛片在线| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 18禁国产床啪视频网站| 大码成人一级视频| 精品亚洲成国产av| 菩萨蛮人人尽说江南好唐韦庄| 精品少妇一区二区三区视频日本电影| 99国产精品一区二区蜜桃av | 国产精品一区二区免费欧美 | 天堂8中文在线网| 天天躁狠狠躁夜夜躁狠狠躁| 性少妇av在线| 在线亚洲精品国产二区图片欧美| 久久九九热精品免费| 一边摸一边抽搐一进一出视频| 90打野战视频偷拍视频| 一级毛片精品| 91精品三级在线观看| www.精华液| 99国产极品粉嫩在线观看| 亚洲人成电影观看| 国产精品久久久久成人av| 一二三四社区在线视频社区8| 777米奇影视久久| 国产成人欧美| 久久性视频一级片| 麻豆乱淫一区二区| 亚洲精品自拍成人| 美女国产高潮福利片在线看| 午夜福利在线免费观看网站| 国产1区2区3区精品| av视频免费观看在线观看| av片东京热男人的天堂| 男人舔女人的私密视频| 亚洲色图 男人天堂 中文字幕| 97人妻天天添夜夜摸| 日韩视频一区二区在线观看| 看免费av毛片| 91精品三级在线观看| 欧美另类一区| 欧美精品高潮呻吟av久久| 99久久国产精品久久久| 纯流量卡能插随身wifi吗| 国产成人影院久久av| 成年人午夜在线观看视频| 黄色视频,在线免费观看| 一区二区三区四区激情视频| 99精品欧美一区二区三区四区| 2018国产大陆天天弄谢| 男人添女人高潮全过程视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲专区中文字幕在线| 99热国产这里只有精品6| 97人妻天天添夜夜摸| 午夜福利,免费看| 亚洲精品久久成人aⅴ小说| 久久亚洲国产成人精品v| 十八禁人妻一区二区| 日本av手机在线免费观看| 999久久久精品免费观看国产| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 中文字幕精品免费在线观看视频| 欧美黑人精品巨大| 欧美亚洲 丝袜 人妻 在线| 日日爽夜夜爽网站| 一个人免费看片子| 国产成人精品久久二区二区91| 男女高潮啪啪啪动态图| 婷婷色av中文字幕| 亚洲av日韩精品久久久久久密| 各种免费的搞黄视频| 日韩一区二区三区影片| 捣出白浆h1v1| 亚洲国产精品一区三区| av天堂在线播放| 久久精品成人免费网站| 热99re8久久精品国产| 午夜激情久久久久久久| 下体分泌物呈黄色| 国产亚洲av片在线观看秒播厂| 日本黄色日本黄色录像| av福利片在线| av国产精品久久久久影院| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 日本猛色少妇xxxxx猛交久久| 免费日韩欧美在线观看| 国产免费av片在线观看野外av| 制服诱惑二区| 天天躁日日躁夜夜躁夜夜| 国产精品国产av在线观看| 99精品久久久久人妻精品| 五月天丁香电影| 国产免费av片在线观看野外av| 欧美日韩一级在线毛片| www日本在线高清视频| 亚洲国产精品999| 别揉我奶头~嗯~啊~动态视频 | 欧美变态另类bdsm刘玥| 免费观看a级毛片全部| 国产在视频线精品| 久久毛片免费看一区二区三区| 亚洲伊人久久精品综合| 十八禁高潮呻吟视频| 日本一区二区免费在线视频| 在线观看舔阴道视频| 少妇被粗大的猛进出69影院| 久久久久精品国产欧美久久久 | 18禁裸乳无遮挡动漫免费视频| 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 亚洲伊人色综图| 女警被强在线播放| 男人爽女人下面视频在线观看| 搡老岳熟女国产| 国产亚洲av高清不卡| 久久毛片免费看一区二区三区| 亚洲一区二区三区欧美精品| 国产av国产精品国产| 伊人久久大香线蕉亚洲五| 精品福利永久在线观看| 久久99热这里只频精品6学生| 波多野结衣av一区二区av| 国产成人影院久久av| 人妻人人澡人人爽人人| 欧美黄色片欧美黄色片| 亚洲av男天堂| av免费在线观看网站| 少妇被粗大的猛进出69影院| 精品久久久久久电影网| 国产成人精品久久二区二区免费| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 老司机在亚洲福利影院| 亚洲九九香蕉| 黄片播放在线免费| 丝袜喷水一区| 91国产中文字幕| 热99国产精品久久久久久7| 一二三四社区在线视频社区8| av又黄又爽大尺度在线免费看| 婷婷色av中文字幕| 中文字幕人妻丝袜制服| 久久综合国产亚洲精品| 亚洲欧美色中文字幕在线| 男女免费视频国产| 久久久久精品国产欧美久久久 | 性色av一级| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 人妻人人澡人人爽人人| 18禁国产床啪视频网站| 一区二区三区精品91| 色视频在线一区二区三区| 日韩欧美国产一区二区入口| 黄色视频,在线免费观看| 免费不卡黄色视频| 无限看片的www在线观看| 久久精品aⅴ一区二区三区四区| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频| 精品人妻1区二区| 一级片免费观看大全| 精品第一国产精品| 成人av一区二区三区在线看 | 国产成人免费观看mmmm| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 久久性视频一级片| 丝袜美足系列| 青春草视频在线免费观看| 欧美成狂野欧美在线观看| 欧美日韩国产mv在线观看视频| 欧美在线黄色| 中文字幕另类日韩欧美亚洲嫩草| √禁漫天堂资源中文www| 丰满饥渴人妻一区二区三| 两个人免费观看高清视频| 69精品国产乱码久久久| 亚洲欧美清纯卡通| 欧美老熟妇乱子伦牲交| a级毛片在线看网站| 久久热在线av| 国产日韩欧美亚洲二区| 18在线观看网站| 精品国产国语对白av| 考比视频在线观看| 亚洲天堂av无毛| 色播在线永久视频| 免费黄频网站在线观看国产| 亚洲av国产av综合av卡| 一本色道久久久久久精品综合| 午夜福利影视在线免费观看| 国产精品久久久人人做人人爽| 成人手机av| 91精品伊人久久大香线蕉| 美女午夜性视频免费| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 亚洲国产欧美一区二区综合| 久久久久久久精品精品| 国产在线观看jvid| 后天国语完整版免费观看| 不卡av一区二区三区| 日本黄色日本黄色录像| 在线av久久热| 亚洲熟女毛片儿| 亚洲性夜色夜夜综合| 搡老乐熟女国产| 人人妻,人人澡人人爽秒播| 最近最新免费中文字幕在线| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区| 99香蕉大伊视频| 下体分泌物呈黄色| 国产精品九九99| 亚洲av电影在线观看一区二区三区| 久久精品熟女亚洲av麻豆精品| 日韩中文字幕视频在线看片| 99国产极品粉嫩在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 韩国高清视频一区二区三区| 中亚洲国语对白在线视频| 久久精品亚洲熟妇少妇任你| 久久国产精品大桥未久av| 精品国内亚洲2022精品成人 | 国产麻豆69| 亚洲三区欧美一区| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 男人添女人高潮全过程视频| 精品一区二区三区四区五区乱码| 在线观看www视频免费| 久久精品亚洲av国产电影网| 一本久久精品| 亚洲视频免费观看视频| 免费不卡黄色视频| 国产视频一区二区在线看| 国产精品自产拍在线观看55亚洲 | 日韩一卡2卡3卡4卡2021年| 亚洲五月色婷婷综合| 男女免费视频国产| 亚洲成人手机| 一区福利在线观看| 免费人妻精品一区二区三区视频| 精品久久久精品久久久| 捣出白浆h1v1| 欧美精品av麻豆av| 欧美精品高潮呻吟av久久| 日韩有码中文字幕| 亚洲精品第二区| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 91字幕亚洲| 一区二区三区精品91| 欧美性长视频在线观看| 日本黄色日本黄色录像| 黑人猛操日本美女一级片| 精品亚洲成a人片在线观看| 国产熟女午夜一区二区三区| 老司机影院成人| 日韩,欧美,国产一区二区三区| 大码成人一级视频| 韩国精品一区二区三区| 国产黄频视频在线观看| 成年人午夜在线观看视频| 性高湖久久久久久久久免费观看| 多毛熟女@视频| 啦啦啦免费观看视频1| 日韩制服骚丝袜av| 久久久精品区二区三区| 国产精品久久久久久人妻精品电影 | 亚洲精品美女久久av网站| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 欧美人与性动交α欧美精品济南到| 一级黄色大片毛片| 男女下面插进去视频免费观看| 午夜福利视频精品| 宅男免费午夜| 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 99久久国产精品久久久| 亚洲国产欧美日韩在线播放| 国产黄频视频在线观看| 欧美精品亚洲一区二区| 在线看a的网站| 国产成人精品久久二区二区免费| 国产在线视频一区二区| 伦理电影免费视频| 午夜免费观看性视频| 欧美一级毛片孕妇| videos熟女内射| 精品久久久久久电影网| 亚洲全国av大片| 久久香蕉激情| 婷婷成人精品国产| 久久久久国内视频| 久久久久国产精品人妻一区二区| 满18在线观看网站| 丰满少妇做爰视频| 十八禁网站免费在线| 2018国产大陆天天弄谢| 飞空精品影院首页| 在线精品无人区一区二区三| 欧美另类亚洲清纯唯美| 午夜精品国产一区二区电影| 国产免费福利视频在线观看| 日日摸夜夜添夜夜添小说| 欧美黑人欧美精品刺激| 91大片在线观看| 久久毛片免费看一区二区三区| 国产男女内射视频| 男人添女人高潮全过程视频| 在线av久久热| av在线老鸭窝| av片东京热男人的天堂| 精品亚洲乱码少妇综合久久| 亚洲欧美精品自产自拍| 99国产精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 丰满少妇做爰视频| 一个人免费在线观看的高清视频 | 久久久久久免费高清国产稀缺| 18在线观看网站| 亚洲精品在线美女| 日韩有码中文字幕| 丁香六月天网| √禁漫天堂资源中文www| 人人妻人人爽人人添夜夜欢视频| 9色porny在线观看| 亚洲中文av在线| 咕卡用的链子| 欧美在线一区亚洲| 亚洲熟女毛片儿| 亚洲人成电影免费在线| 涩涩av久久男人的天堂| 一级片'在线观看视频| kizo精华| 日日爽夜夜爽网站| 黄色视频,在线免费观看| 高清欧美精品videossex| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 美女大奶头黄色视频| 国产日韩欧美视频二区| 国产av又大| 丝瓜视频免费看黄片| 亚洲成国产人片在线观看| 成年动漫av网址| 91大片在线观看| 美女高潮喷水抽搐中文字幕| 黄片播放在线免费| 久久久久精品人妻al黑| 亚洲精品av麻豆狂野| 亚洲av电影在线进入| 一本大道久久a久久精品| 51午夜福利影视在线观看| 亚洲午夜精品一区,二区,三区| av超薄肉色丝袜交足视频| 一边摸一边做爽爽视频免费| 欧美亚洲 丝袜 人妻 在线| 美女视频免费永久观看网站| 最黄视频免费看| 亚洲av成人不卡在线观看播放网 | 飞空精品影院首页| 视频区图区小说| 午夜福利视频在线观看免费| 国产亚洲欧美精品永久| 国产精品久久久人人做人人爽| 黄色a级毛片大全视频|