• <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久久男人的天堂| 亚洲人成电影免费在线| 色尼玛亚洲综合影院| 99久久99久久久精品蜜桃| 国产视频一区二区在线看| 美女大奶头视频| 久久久久久国产a免费观看| 嫁个100分男人电影在线观看| 久久精品国产99精品国产亚洲性色 | 欧美乱码精品一区二区三区| 人人妻人人澡欧美一区二区 | 中文字幕精品免费在线观看视频| 国产在线观看jvid| 两个人视频免费观看高清| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区av网在线观看| 国产精品亚洲一级av第二区| 黄片小视频在线播放| 亚洲精华国产精华精| 亚洲欧美精品综合久久99| 香蕉久久夜色| 亚洲第一av免费看| 两个人视频免费观看高清| 少妇的丰满在线观看| 黑人操中国人逼视频| 999久久久国产精品视频| 亚洲成av片中文字幕在线观看| 99riav亚洲国产免费| 午夜久久久在线观看| 天天一区二区日本电影三级 | √禁漫天堂资源中文www| 免费在线观看视频国产中文字幕亚洲| 午夜精品国产一区二区电影| 亚洲熟妇中文字幕五十中出| 日韩高清综合在线| 成人18禁在线播放| 老熟妇乱子伦视频在线观看| 丁香欧美五月| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 日韩国内少妇激情av| 涩涩av久久男人的天堂| 亚洲性夜色夜夜综合| 99久久国产精品久久久| 看免费av毛片| 精品一区二区三区视频在线观看免费| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 桃色一区二区三区在线观看| 丝袜美足系列| 亚洲专区中文字幕在线| 在线天堂中文资源库| 制服诱惑二区| 欧美久久黑人一区二区| 中文亚洲av片在线观看爽| 国产精品久久久人人做人人爽| 1024视频免费在线观看| 日本在线视频免费播放| 亚洲午夜精品一区,二区,三区| 国产又色又爽无遮挡免费看| 欧美日韩亚洲综合一区二区三区_| 国产极品粉嫩免费观看在线| 亚洲 欧美一区二区三区| 麻豆久久精品国产亚洲av| 99精品欧美一区二区三区四区| 男人的好看免费观看在线视频 | 成人亚洲精品av一区二区| 99国产精品一区二区三区| 精品久久久久久,| 精品国产乱码久久久久久男人| 国产97色在线日韩免费| 咕卡用的链子| 国产三级在线视频| 亚洲九九香蕉| 亚洲一区高清亚洲精品| 久久精品亚洲精品国产色婷小说| 在线视频色国产色| 日本a在线网址| 婷婷丁香在线五月| 老司机在亚洲福利影院| 两个人免费观看高清视频| 亚洲精华国产精华精| 午夜久久久久精精品| 夜夜看夜夜爽夜夜摸| 黄片小视频在线播放| 亚洲精品久久国产高清桃花| 欧美色欧美亚洲另类二区 | www.自偷自拍.com| 国产主播在线观看一区二区| av视频在线观看入口| 欧美成人免费av一区二区三区| 在线十欧美十亚洲十日本专区| 国产xxxxx性猛交| cao死你这个sao货| 久久久久精品国产欧美久久久| 麻豆成人av在线观看| 欧美老熟妇乱子伦牲交| 久久性视频一级片| 久久青草综合色| 国产精品综合久久久久久久免费 | 久久国产精品人妻蜜桃| 妹子高潮喷水视频| 啦啦啦 在线观看视频| bbb黄色大片| 99久久99久久久精品蜜桃| 日韩成人在线观看一区二区三区| 搡老妇女老女人老熟妇| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 亚洲第一青青草原| 午夜精品在线福利| 99国产精品99久久久久| 婷婷精品国产亚洲av在线| 亚洲av电影在线进入| 夜夜夜夜夜久久久久| 国产片内射在线| 曰老女人黄片| 人人妻人人澡欧美一区二区 | 又黄又粗又硬又大视频| 男女下面进入的视频免费午夜 | 很黄的视频免费| 咕卡用的链子| 久久国产精品男人的天堂亚洲| 一级,二级,三级黄色视频| 久久人妻av系列| 悠悠久久av| 精品久久久久久久毛片微露脸| 欧美黄色淫秽网站| 9热在线视频观看99| 免费av毛片视频| 国产熟女午夜一区二区三区| 真人做人爱边吃奶动态| 男女午夜视频在线观看| 免费人成视频x8x8入口观看| 国产人伦9x9x在线观看| 国产精品综合久久久久久久免费 | 日日摸夜夜添夜夜添小说| 精品人妻1区二区| 男女下面插进去视频免费观看| 久久久久精品国产欧美久久久| 黄色 视频免费看| 淫秽高清视频在线观看| 国产成人欧美在线观看| 老熟妇乱子伦视频在线观看| 黄片小视频在线播放| 黄色片一级片一级黄色片| 亚洲一区二区三区不卡视频| 午夜a级毛片| 1024视频免费在线观看| 国产99久久九九免费精品| 乱人伦中国视频| 最近最新中文字幕大全免费视频| 亚洲第一电影网av| 欧美亚洲日本最大视频资源| 露出奶头的视频| 精品欧美一区二区三区在线| 精品少妇一区二区三区视频日本电影| 91精品三级在线观看| 性少妇av在线| 777久久人妻少妇嫩草av网站| 亚洲av美国av| 国产成人一区二区三区免费视频网站| 国产精品乱码一区二三区的特点 | 国产成人精品久久二区二区免费| 搡老岳熟女国产| 男人的好看免费观看在线视频 | 日韩大尺度精品在线看网址 | 高清在线国产一区| 免费在线观看影片大全网站| 99国产精品免费福利视频| 免费不卡黄色视频| 9热在线视频观看99| 十八禁网站免费在线| e午夜精品久久久久久久| 老熟妇乱子伦视频在线观看| 天堂√8在线中文| 国产亚洲欧美98| 又黄又粗又硬又大视频| 亚洲激情在线av| 在线观看免费午夜福利视频| 超碰成人久久| 91麻豆精品激情在线观看国产| 三级毛片av免费| 手机成人av网站| 天堂动漫精品| 欧美丝袜亚洲另类 | 99久久国产精品久久久| 亚洲一区中文字幕在线| 亚洲久久久国产精品| 最新在线观看一区二区三区| 一级a爱视频在线免费观看| 亚洲熟妇熟女久久| 欧美中文日本在线观看视频| 精品少妇一区二区三区视频日本电影| 国产精品日韩av在线免费观看 | 国产一区二区三区综合在线观看| 女人被狂操c到高潮| 精品不卡国产一区二区三区| 久久久久久国产a免费观看| 午夜福利成人在线免费观看| 不卡一级毛片| 久久人妻福利社区极品人妻图片| 正在播放国产对白刺激| 久久欧美精品欧美久久欧美| 18美女黄网站色大片免费观看| 黄色毛片三级朝国网站| av欧美777| 久久久国产成人免费| 不卡av一区二区三区| 岛国视频午夜一区免费看| 成人国语在线视频| 亚洲精品中文字幕在线视频| 两个人视频免费观看高清| 免费在线观看日本一区| 免费少妇av软件| 好男人在线观看高清免费视频 | 久久精品aⅴ一区二区三区四区| 免费在线观看黄色视频的| 色老头精品视频在线观看| 老汉色av国产亚洲站长工具| 久久精品国产99精品国产亚洲性色 | 亚洲欧美日韩高清在线视频| 女性被躁到高潮视频| 国产精品免费一区二区三区在线| 日本撒尿小便嘘嘘汇集6| 午夜两性在线视频| 国产精品av久久久久免费| 天堂影院成人在线观看| 亚洲精品在线观看二区| 十八禁人妻一区二区| 咕卡用的链子| 深夜精品福利| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 香蕉丝袜av| 九色国产91popny在线| 国产野战对白在线观看| 成年人黄色毛片网站| 日韩国内少妇激情av| 亚洲熟妇熟女久久| 国产精品 欧美亚洲| 自线自在国产av| 亚洲五月天丁香| 久久久久久久久中文| 国产精品久久久人人做人人爽| 日本撒尿小便嘘嘘汇集6| 动漫黄色视频在线观看| 国产伦人伦偷精品视频| 日日爽夜夜爽网站| 欧美黑人精品巨大| 免费不卡黄色视频| 妹子高潮喷水视频| 免费人成视频x8x8入口观看| 亚洲国产欧美网| 国产精品免费一区二区三区在线| 可以免费在线观看a视频的电影网站| 亚洲中文字幕一区二区三区有码在线看 | 精品一区二区三区av网在线观看| 99re在线观看精品视频| 国产国语露脸激情在线看| 少妇 在线观看| 精品国产乱码久久久久久男人| 久久青草综合色| 少妇被粗大的猛进出69影院| 日本撒尿小便嘘嘘汇集6| 999久久久精品免费观看国产| а√天堂www在线а√下载| 国产精品电影一区二区三区| 国产av一区二区精品久久| 国产区一区二久久| 国产成人欧美| av网站免费在线观看视频| 国产成人精品久久二区二区免费| 深夜精品福利| 色尼玛亚洲综合影院| 亚洲欧美日韩无卡精品| 别揉我奶头~嗯~啊~动态视频| 在线观看66精品国产| 国产伦人伦偷精品视频| 一区二区三区国产精品乱码| 国产91精品成人一区二区三区| 久久精品国产99精品国产亚洲性色 | 国产精品九九99| 最新美女视频免费是黄的| 99re在线观看精品视频| aaaaa片日本免费| 国产午夜福利久久久久久| 亚洲国产欧美一区二区综合| 久久精品亚洲精品国产色婷小说| 国产精品九九99| xxx96com| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 午夜福利高清视频| 桃红色精品国产亚洲av| 又黄又爽又免费观看的视频| or卡值多少钱| 国产激情欧美一区二区| 人人妻人人爽人人添夜夜欢视频| 亚洲avbb在线观看| 亚洲少妇的诱惑av| 久久香蕉激情| 亚洲欧美激情在线| 最近最新免费中文字幕在线| 搞女人的毛片| 天天添夜夜摸| 啦啦啦韩国在线观看视频| 他把我摸到了高潮在线观看| 欧美另类亚洲清纯唯美| 亚洲欧美日韩无卡精品| 亚洲欧美激情综合另类| 天天躁夜夜躁狠狠躁躁| 欧美乱妇无乱码| 国产亚洲av高清不卡| 色播亚洲综合网| 美国免费a级毛片| videosex国产| 中文字幕久久专区| 18禁国产床啪视频网站| 色尼玛亚洲综合影院| 国产精品香港三级国产av潘金莲| 桃色一区二区三区在线观看| 嫩草影视91久久| 亚洲色图 男人天堂 中文字幕| 一级黄色大片毛片| 露出奶头的视频| 窝窝影院91人妻| 99国产精品免费福利视频| 久久影院123| 在线观看免费视频网站a站| 最新在线观看一区二区三区| 婷婷精品国产亚洲av在线| bbb黄色大片| 老司机靠b影院| 午夜免费激情av| 少妇 在线观看| 夜夜爽天天搞| 亚洲av五月六月丁香网| 国产精品爽爽va在线观看网站 | 亚洲最大成人中文| 欧美日韩乱码在线| 亚洲国产日韩欧美精品在线观看 | 国产亚洲欧美98| 神马国产精品三级电影在线观看 | 亚洲成av人片免费观看| 亚洲中文字幕日韩| 中亚洲国语对白在线视频| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久久5区| 国产亚洲精品综合一区在线观看 | 99热只有精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 最好的美女福利视频网| 久久久精品国产亚洲av高清涩受| 国产高清视频在线播放一区| 久久久国产成人免费| 身体一侧抽搐| 国产av又大| 欧美日韩黄片免| 午夜福利视频1000在线观看 | 极品人妻少妇av视频| 韩国精品一区二区三区| 很黄的视频免费| 国产精品影院久久| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 中文字幕人妻熟女乱码| 12—13女人毛片做爰片一| 国产欧美日韩一区二区精品| 悠悠久久av| 99香蕉大伊视频| 欧美亚洲日本最大视频资源| 麻豆国产av国片精品| 每晚都被弄得嗷嗷叫到高潮| 欧美人与性动交α欧美精品济南到| 激情视频va一区二区三区| 午夜免费激情av| 精品国产亚洲在线| 真人一进一出gif抽搐免费| 在线观看www视频免费| 午夜免费激情av| 精品国产超薄肉色丝袜足j| 亚洲五月婷婷丁香| 亚洲电影在线观看av| 十八禁人妻一区二区| 黄网站色视频无遮挡免费观看| 嫁个100分男人电影在线观看| 欧美精品啪啪一区二区三区| av有码第一页| 色综合欧美亚洲国产小说| 黄网站色视频无遮挡免费观看| 成在线人永久免费视频| 人妻丰满熟妇av一区二区三区| 国产91精品成人一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲激情在线av| 九色国产91popny在线| 中文字幕人妻丝袜一区二区| 精品国产乱子伦一区二区三区| 99国产精品免费福利视频| 欧美最黄视频在线播放免费| av有码第一页| 咕卡用的链子| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 久久久久国内视频| 身体一侧抽搐| 人人妻人人澡人人看| 97碰自拍视频| 欧美丝袜亚洲另类 | 日日夜夜操网爽| 亚洲中文av在线| 亚洲中文日韩欧美视频| 嫩草影院精品99| 一本综合久久免费| а√天堂www在线а√下载| 激情在线观看视频在线高清| 视频在线观看一区二区三区| 日韩中文字幕欧美一区二区| 后天国语完整版免费观看| 99久久久亚洲精品蜜臀av| 人人妻人人澡人人看| 国产三级在线视频| 成人三级做爰电影| 欧美黑人精品巨大| 国产在线精品亚洲第一网站| 亚洲精品美女久久久久99蜜臀| 欧美最黄视频在线播放免费| 亚洲人成77777在线视频| 韩国av一区二区三区四区| 一夜夜www| 九色亚洲精品在线播放| a在线观看视频网站| 悠悠久久av| 亚洲人成伊人成综合网2020| 欧美日本亚洲视频在线播放| 可以免费在线观看a视频的电影网站| 日韩精品免费视频一区二区三区| 国产片内射在线| 亚洲av成人av| 一区二区三区国产精品乱码| 久久人妻熟女aⅴ| av免费在线观看网站| a在线观看视频网站| 午夜视频精品福利| 国产亚洲欧美98| 99在线视频只有这里精品首页| 亚洲最大成人中文| 国产精品自产拍在线观看55亚洲| 人妻丰满熟妇av一区二区三区| 国产真人三级小视频在线观看| 久久久精品国产亚洲av高清涩受| 亚洲成人免费电影在线观看| 亚洲av熟女| 美国免费a级毛片| 久久精品人人爽人人爽视色| 久久人人爽av亚洲精品天堂| 两个人免费观看高清视频| 色综合欧美亚洲国产小说| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 99久久久亚洲精品蜜臀av| 如日韩欧美国产精品一区二区三区| 女同久久另类99精品国产91| 高清黄色对白视频在线免费看| 又黄又粗又硬又大视频| 天天添夜夜摸| 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 亚洲天堂国产精品一区在线| 午夜a级毛片| 婷婷丁香在线五月| 波多野结衣av一区二区av| 久久精品国产99精品国产亚洲性色 | 88av欧美| 亚洲av成人一区二区三| 日韩av在线大香蕉| 国产一级毛片七仙女欲春2 | 亚洲一区二区三区不卡视频| 一边摸一边抽搐一进一小说| 天天躁夜夜躁狠狠躁躁| 国产亚洲欧美精品永久| 国产欧美日韩精品亚洲av| 咕卡用的链子| 在线观看午夜福利视频| 国产欧美日韩一区二区三区在线| 少妇的丰满在线观看| 宅男免费午夜| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 丝袜在线中文字幕| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| 啦啦啦韩国在线观看视频| 日本免费a在线| 国产精品免费视频内射| bbb黄色大片| 男人舔女人下体高潮全视频| 亚洲色图综合在线观看| av在线天堂中文字幕| 免费无遮挡裸体视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美日韩无卡精品| 亚洲情色 制服丝袜| 久久久国产精品麻豆| 女性被躁到高潮视频| 精品欧美国产一区二区三| 亚洲精品在线观看二区| 黄色丝袜av网址大全| 19禁男女啪啪无遮挡网站| 国产成+人综合+亚洲专区| 18美女黄网站色大片免费观看| 大型黄色视频在线免费观看| 国产成人啪精品午夜网站| 亚洲精品美女久久久久99蜜臀| 一本大道久久a久久精品| 男人的好看免费观看在线视频 | 日本免费一区二区三区高清不卡 | 欧美黄色片欧美黄色片| 国产精品亚洲美女久久久| 久久精品成人免费网站| 亚洲av成人不卡在线观看播放网| 日本在线视频免费播放| 午夜福利免费观看在线| 97人妻天天添夜夜摸| 女警被强在线播放| av欧美777| 丁香六月欧美| 久久亚洲真实| 97碰自拍视频| 日韩欧美国产在线观看| 国产精品精品国产色婷婷| 国产精品亚洲一级av第二区| 老司机午夜十八禁免费视频| 国产私拍福利视频在线观看| 国产1区2区3区精品| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 99久久综合精品五月天人人| 天堂√8在线中文| 女性被躁到高潮视频| 99香蕉大伊视频| 成人三级黄色视频| 日韩欧美在线二视频| 午夜免费观看网址| 免费高清视频大片| 九色国产91popny在线| 午夜免费观看网址| 人人澡人人妻人| 日本精品一区二区三区蜜桃| 日韩欧美在线二视频| 国产精品香港三级国产av潘金莲| 久久草成人影院| 亚洲自拍偷在线| 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 窝窝影院91人妻| 国产精品 国内视频| 久久久久亚洲av毛片大全| 丰满人妻熟妇乱又伦精品不卡| 久久久国产精品麻豆| 亚洲最大成人中文| 久久久久国内视频| 麻豆久久精品国产亚洲av| 国产欧美日韩一区二区精品| 成人特级黄色片久久久久久久| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| 国产精品久久久久久人妻精品电影| 亚洲一区二区三区不卡视频| 黄色视频,在线免费观看| 97人妻天天添夜夜摸| 免费在线观看亚洲国产| 亚洲男人天堂网一区| 亚洲国产精品合色在线| 欧美日韩福利视频一区二区| 国产主播在线观看一区二区| 精品国产国语对白av| 91国产中文字幕| 日韩 欧美 亚洲 中文字幕| 丁香六月欧美| 丝袜人妻中文字幕| 久久久国产欧美日韩av| 99精品在免费线老司机午夜| 亚洲精华国产精华精| 精品电影一区二区在线| 亚洲一区高清亚洲精品| 久久精品aⅴ一区二区三区四区| 无人区码免费观看不卡| 黑人欧美特级aaaaaa片| 亚洲久久久国产精品| 看免费av毛片| 国产精品爽爽va在线观看网站 | 一级毛片高清免费大全| 视频在线观看一区二区三区| 在线视频色国产色| 免费一级毛片在线播放高清视频 | 欧美+亚洲+日韩+国产| 最近最新中文字幕大全电影3 | 欧美人与性动交α欧美精品济南到| 亚洲片人在线观看| 美国免费a级毛片| av福利片在线| 在线永久观看黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | 丁香欧美五月| 国产精品一区二区免费欧美| 亚洲黑人精品在线| 久久人人精品亚洲av| 啦啦啦观看免费观看视频高清 | 国产高清视频在线播放一区| 亚洲欧洲精品一区二区精品久久久| 欧美激情高清一区二区三区| 日韩免费av在线播放| 级片在线观看|