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

    多因子協(xié)同作用下的云南省山火風險評估

    2022-04-28 12:14:50文剛高孟平高振宇周仿榮潘浩李志棟張濤羅繼強劉星星
    廣東電力 2022年4期
    關鍵詞:火點山火植被

    文剛,高孟平,高振宇,周仿榮,潘浩,李志棟,張濤,羅繼強,劉星星

    (1.電力遙感技術聯(lián)合實驗室(云南電網有限責任公司電力科學研究院),云南 昆明 650217;2.云南電網有限責任公司,云南 昆明 650200;3.北京空間飛行器總體設計部,北京 100094;4.廣東工業(yè)大學,廣東 廣州 510006)

    云南省的大部分架空輸電線路處于環(huán)境復雜的崇山峻嶺中,輸電線路遇風很容易擺動,與樹木或山體坡面頻繁接觸,引起短路放電造成山火[1]。一旦有山火發(fā)生,火勢很容易在架空線路走廊周邊蔓延,使輸電線路暫停運作,同時由于山火持續(xù)時間一般較長,導致絕大部分重合閘失敗,進而造成停電事故。當山火蔓延影響到電網區(qū)域中多條輸電線路時,很可能會造成大范圍停電事故,嚴重影響電網的安全穩(wěn)定運行[2]。為了盡可能減少山火隱患,防止電網的安全穩(wěn)定運行受到山火威脅,有必要開展山火風險評估的研究。

    山火的發(fā)生受植被類型、植被狀態(tài)(含水率)、氣象因素、地形、人為因素等影響,山火風險預測系統(tǒng)是受多種誘導因子作用的復雜系統(tǒng)。目前,山火評估預測的方法主要分為傳統(tǒng)方法和模擬法。

    傳統(tǒng)方法包括相關關系法、統(tǒng)計分析法、以火報火法。H.t.Gisborne在1936年提出了多因子預測火險等級的相關關系預測方法[3]。統(tǒng)計分析法是通過收集歷史火災相關氣象資料資料,統(tǒng)計分析歷史火災的天氣狀況、時間、地形、次數等氣象因子,確定火災與各個氣象因子的相關性,總結火災發(fā)生的規(guī)律。但是氣象因子對火災發(fā)生的影響機制十分復雜,同時火災發(fā)生也與人類活動有一定關系,因而統(tǒng)計分析法對火災的評估預測可靠性比較低[4]。以火報火法通過建立不同森林可燃物的點火燃燒實驗與氣象要素之間的關系來預測火災的危險性,具有代表性的是澳大利亞曾進行點火試驗800次,建立了森林火災預報模型,該系統(tǒng)由McArthur設計,不僅能預測林火的可能性,還能輸出定量的火行為參數,例如林火蔓延速度、火險強度、飛火距離等[5]。該方法預測火災的危險性較為準確,但需要進行大量的野外試驗,建立模型需要花費大量的時間。

    隨著計算機技術的迅猛發(fā)展及其在火災評估中的應用,模擬法利用衛(wèi)星監(jiān)測技術和地理信息系統(tǒng)(geographic information system,GIS)技術,使山火預測更為準確。其中Logistic回歸是應用最為廣泛的模型。文獻[6]基于Logistic回歸模型對浙江省2000—2016年林火發(fā)生的主要影響因素進行分析,確定驅動因子,建立林火預測模型并劃分火險等級,模型總體預測準確率達到79.1%。文獻[7]以江西省為研究對象,利用Logistic回歸模型對高程、坡度、距道路最近距離、可燃物含水率、陸地表面溫度、歸一化植被指數、全局植被濕度指數進行訓練,接受者操作特征(receiver operating characteristic,ROC)曲線下面積(area under ROC curve,AUC)值達到0.757。文獻[8]基于衛(wèi)星火點數據和氣象數據,用Logistic回歸模型對浙江省林火變化趨勢和影響因子進行研究。文獻[9]運用加權Logistic回歸建立廣東省輸電走廊林火預警模型。另外,人工神經網絡(artificial neural network,ANN)也在該領域得到應用。文獻[10]收集歷史氣候數據,采用BP神經網絡進行預測。文獻[11]采用中分辨率成像光譜儀(moderate-resolution imaging spectroradiometer,MODIS)數據,以西班牙加利西亞為研究對象,比較ANN和Logistic回歸模型在森林火災評估中的優(yōu)劣,結果顯示二者的評估精度接近。

    ANN本質是一種擬合技術,由于山火訓練數據樣本只有0和1,所以模型評估結果大多落在0和1附近,而落在0.4~0.6之間的值很少[12]。Logistic回歸模型在山火評估領域中的應用更加廣泛,因此本文基于Logistic回歸模型選取更加全面的山火評價因子對山火進行評估,并進行Person相關性分析,剔除冗余的評價因子。最后,將本文所提協(xié)同多因子山火風險評估方法與文獻[7,13-14]的方法對比,以驗證所提方法的優(yōu)越性。

    1 研究區(qū)域和數據準備

    1.1 研究區(qū)域

    云南省位于中國西南地區(qū),介于北緯21°8′~29°15′、東經97°31′~106°11′之間。地形以山地、高原為主,山地面積占全省總面積的88.64%,森林覆蓋率達到65.04%。云南氣候基本屬于亞熱帶和熱帶季風氣候,這種氣候年溫差小、日溫差大,降水在空間和時間上分布也極不均勻。云南省的森林面積分布廣泛,而且氣候很容易引起山火,導致省內山火頻繁,給電網設備等人民財產造成了重大損失。本文以云南省為研究區(qū)域,進行山火風險評估研究。

    1.2 數據準備

    本文通過文獻調研,選取直接或間接的山火影響因子,數據集屬性見表1。

    表1 數據集屬性Tab.1 Data set attributes

    本文從地形、氣象、植被和人為因素4個方面選取山火影響因子,各個因子的詳細說明如下。

    a)歷史火點數據來源于NASA網站MODIS傳感器偵測到的火點(https://firms.modaps.eosdis.nasa.gov/map/),根據云南省每天的歷史火點提取火點位置的17個山火評價因子,同時隨機選擇非火點,并提取非火點位置的17個山火評價因子。

    b)地表溫度數據來源于3級水平MODIS地表溫度數據MOD11A1,此數據是具有1 km空間分辨率的正弦投影數據。地表溫度高的區(qū)域土壤水分少,植被干燥,容易發(fā)生山火;反之,地表溫度低的區(qū)域土壤水分相對多,不容易發(fā)生山火。

    c)歸一化紅外指數(NDII7)來源于MODIS表面反射率數據系列,選取500 m空間分辨率、8 d時間分辨率的3級正弦投影數據。該數據具有7個波段,通過波段2和波段7計算得到NDII7遙感影像。NDII7小的區(qū)域植被含水率低,容易引起火災;反之,NDII7大的區(qū)域植被含水率高,不容易引起火災。文獻[15]通過實驗和科學研究發(fā)現(xiàn),NDII7較小的干燥植被比NDII7較大的新鮮植被更容易發(fā)生燃燒。因此,本文用NDII7來表征可燃物的含水率。

    d)歸一化指標指數(NDVI)采用全球MODIS植被指數數據提供的時空連續(xù)的植被監(jiān)測數據[16-17]。采用正弦投影和16 d合成的500 m空間分辨率的植被數據,用于檢測植被狀態(tài)和土地覆蓋率變化情況。

    e)高程數據描述一個地方的海拔高度,不同高程對山火發(fā)生的影響程度不一樣。高海拔的地方溫度比較低,空氣濕度大,發(fā)生山火的可能性比較低;反之,低海拔的地方溫度比較高,空氣濕度小,發(fā)生山火的可能性大[18]。

    f)坡度描述地面的傾斜程度,受重力等因素的影響,坡度大的地方植被的保水能力小,水土流失比較嚴重,更容易發(fā)生山火。坡向描述一個地方的朝向,不同坡向決定著對陽光照射的接受程度,接受陽光多的坡向發(fā)生山火的可能性更高。根據高程,利用ARCGIS軟件的3D分析功能可得到研究區(qū)域的坡度和坡向。

    g)距道路距離和距離居民的距離反映人為因素對山火的影響關系。距離道路和居民地近的地方,人類活動頻繁,容易引起山火;距離水體近的區(qū)域土壤濕度高,不容易引發(fā)山火[19]。

    h)不同的土地類型對山火發(fā)生的影響程度不同,森林、灌木叢和覆蓋度高的草原更容易引發(fā)山火。

    i)氣象對于山火的發(fā)生有重要的影響,選取降水量、空氣濕度、平均溫度、最小溫度、最大溫度、最大陣風風速、最大陣風風向7個氣象因子。氣象要素的精度對于山火評估有重要影響,文獻[20]使用基于ANUSPLIN的樣條插值方法對全國667個氣象站數據進行插值,利用1 667個加密站點對氣象插值結果進行準確性檢驗,并與反向距離權重法和普通克里金插值方法進行實驗對比,結果表明基于ANUSPLIN的樣條插值方法誤差最小。因此,本文選擇該方法對研究區(qū)域的125個氣象站獲得的氣象資料進行插值,得到整個研究區(qū)域的氣象數據。

    2 評估模型與評價指標

    本文提出的山火評估流程如圖1所示,評估模型為Logistic回歸模型,系統(tǒng)性能指標為AUC。

    圖1 山火評估流程Fig.1 Wildfire assessment flow chart

    2.1 評估模型

    Logistic回歸模型是線性回歸分析模型,常用于數據分析、金融、醫(yī)學等領域。通過Logistic回歸分析可以得到自變量的權重,輸出0~1之間概率值,模型根據結果是否大于0.5分類為0或1。本文基于山火影響因子,通過Logistic回歸模型計算山火風險概率,計算公式為:

    (1)

    z=β0+β1x1+β2x2+…+βnxn.

    (2)

    式(1)、(2)中:Pz為發(fā)生山火風險概率;β0為模型中常數,表示方程的截距;n為自變量的個數,這里表示山火影響因子個數;xi(i=1,2,…,n)為山火影響因子種類,βi為模型中與xi對應的變量參數。

    2.2 評價指標AUC

    ROC曲線可用于揭示敏感性(真正率)與特異性(假正率)之間的相互關系。通過設置連續(xù)變量的閾值,計算出敏感性a和b,以a為縱坐標、以(1-b)為橫坐標繪制而成的曲線即為ROC曲線。AUC可以評估預測模型的質量,其值介于0.5~1:當AUC為0.5時,說明模型的預測效果很差,相當于隨機分類;當AUC接近1時,說明模型的預測效果很好。

    3 實驗與討論

    3.1 相關性分析

    相關性分析是研究模型自變量之間相關程度的分析,衡量2個自變量因素的相關密切程度,目的在于發(fā)掘自變量間的相關特性,是描述客觀事物相互間關系的密切程度并用適當的統(tǒng)計指標表示出來的過程。

    利用統(tǒng)計產品與服務解決方案(statistical product service solutions,SPSS)軟件分析火點評價因子的Person相關性,Person系數接近1時表示有強正相關,接近-1時表示有強負相關。地表溫度與平均氣溫的相關性系數為0.607,平均氣溫與最低氣溫、最高氣溫的相關性系數分別為0.931、0.943,最低氣溫與最高氣溫的相關性系數為0.822。一般認為,相關性系數絕對值達到0.6即為高度相關[12],可利用逐步回歸剔除相關性系數高的其中一個因子。本文通過相關性分析,剔除平均氣溫、最低氣溫2個因子。

    3.2 結果和分析

    3.2.1 評價指標比較

    本文提取MODIS和NPP信息中2018年和2019年所有遙感偵測到的植被大火火點的17個特征,整理得到29 442個樣本火點,同時選取相同數量的非火點樣本,最終數據集包含58 884個樣本,其中70%的樣本用于訓練集,30%的樣本用于測試集。

    用高程、坡度、坡向、植被覆蓋類型、距最近道路距離,距最近水體距離、距最近居民點距離、地表溫度、可燃物含水率、歸一化植被指數、降水量、最高氣溫、平均相對濕度、最大陣風風速、最大陣風風向15個因子來評估云南省山火,與文獻[7,13-14]的山火評估方法進行實驗對比,結果見表2。

    表2 幾種山火評估方法AUC的對比Tab.3 AUC comparisons of different wildfire assessment methods

    文獻[13]采用高程、坡度、坡向、植被覆蓋類型、距最近道路距離、距最近水體距離、距最近居民點距離7個評價因子來評估山火。本文方法計算得到的AUC達到0.80,與文獻[13]方法相比提升了0.16,說明地表溫度、可燃物含水率、歸一化植被指數、降水量、最高氣溫、平均相對濕度、最大陣風風速、最大陣風風向等因子能夠顯著提升評估效果。

    文獻[14]采用高程、坡度、坡向、植被覆蓋類型、距最近道路距離,距最近水體距離、距最近居民點距離、地表溫度、可燃物含水率、歸一化植被指數10個評價因子來評估山火。與文獻[14]方法相比,本文方法計算得到的AUC提升了0.10,說明降水量、最高氣溫、平均相對濕度、最大陣風風速、最大陣風風向是評估山火的有效指標。

    對比文獻[13]和文獻[14],文獻[14]方法考慮地表溫度、植被含水率和歸一化植被指數這3個影響因子,其AUC相比文獻[13]提升了0.06,說明這3個遙感特征用于評估云南省山火并不會出現(xiàn)冗余特征的現(xiàn)象。

    文獻[7]在較大的時間尺度和省級空間尺度上進行森林火災風險分區(qū)研究,基于Logistic回歸建立火災風險評估模型,實驗結果顯示火災分布與火災影響因素之間具有很好的擬合度。但是文獻[7]只采用海拔、土地利用類型、土地表面溫度、歸一化植被指數和全局植被水分指數對火災進行預測評估。本文雖然缺少全局植被水分指數數據,但是引入了歸一化紅外指數來表征植被水分狀態(tài)。與文獻[7]方法相比,本文方法計算得到的AUC提升了0.11,可見氣象要素對于云南省山火評估有著重要作用。

    測試集混淆矩陣見表3。從混淆矩陣可以得到:本文方法非火點的錯誤率為0.21,火點的錯誤率為0.19;文獻[13]方法非火點的錯誤率為0.34,火點的錯誤率為0.38;文獻[14]方法非火點的錯誤率為0.28,火點的錯誤率為0.33;文獻[7]方法非火點的錯誤率為0.28,火點的錯誤率為0.34。從預測準確性來看,本文選取的評價因子更具有優(yōu)勢。

    表3 測試集混淆矩陣Tab.4 Test set confusion matrix

    3.2.2 風險圖比較

    選取不同評價因子,制作2018年1月17日的山火風險概率圖,并將概率按等區(qū)間劃分的原則劃分為極低、低、中等、高、極高5個等級,如圖2所示。0~0.2為山火發(fā)生概率極低;0.2~0.4為山火發(fā)生概率低;0.4~0.6為山火發(fā)生概率中;0.6~0.8為山火發(fā)生概率高;0.8~1為山火發(fā)生概率極高。從圖2來看,雖然大部分火點落在高風險的紅色區(qū)域和黃色區(qū)域,但是文獻[13]、文獻[14]和文獻[7]的山火風險概率圖中,云南省大部分區(qū)域都預測為紅色和黃色的高風險區(qū)域,與實際情況不符合。根據本文選取的15個評價因子預測得到的云南省山火風險概率圖中,大部分區(qū)域為低風險區(qū)域,而且大部分火點也處于山火發(fā)生概率高和極高的高風險區(qū)域,與實際情況相符。因此,可以采用本文的15個評價因子評估云南省的山火風險狀況,結合相應的消防設施資源,做好撲救山火的準備。

    圖2 選取不同評價因子時,云南省2018年1月17日的山火風險概率圖Fig.2 Wildfire risk probability in January 17, 2018 in Yunnan when selecting different evaluation factors

    4 結束語

    云南省的許多輸電線路工程需要穿越山嶺地帶,輸電線路容易受到山火的影響,引發(fā)電力事故,因此有必要進行山火風險評估的研究。為了有效進行山火風險評估,本文基于Logistic回歸模型,從地形、人為因素、遙感和氣象4個方面選取山火評價因子,選取高程、坡度、坡向、植被覆蓋類型、距最近道路距離,距最近水體距離、距最近居民點距離、地表溫度、可燃物含水率、歸一化植被指數、降水量、平均氣溫、最低氣溫、最高氣溫、平均相對濕度、最大陣風風速、最大陣風風向17個山火影響因子,并通過相關性分析剔除平均氣溫和最低氣溫2個冗余因子。實驗結果顯示,所選15個因子能全面、準確地評估云南省的山火狀況。本文選取的15個山火因子實質上是針對山火的空間規(guī)律進行評估,未來可以研究從時間序列的角度出發(fā)進行山火評估。

    猜你喜歡
    火點山火植被
    基于植被復綠技術的孔植試驗及應用
    河北地質(2022年2期)2022-08-22 06:24:04
    亞像元火點對紅外預警衛(wèi)星的輻射干擾特性
    紅外技術(2021年1期)2021-01-29 01:41:54
    人身上有5個祛火點
    美國加州各地山火肆虐 已造成6名消防員死亡
    點煙頌
    鹿鳴(2018年1期)2018-01-30 12:05:42
    綠色植被在溯溪旅游中的應用
    江蘇省海門市如何實現(xiàn)連續(xù)4年秸稈焚燒“零火點”?
    山火難熄
    基于原生植被的長山群島植被退化分析
    基于NDVI的鶴壁市植被覆蓋動態(tài)分析研究
    河南科技(2014年4期)2014-02-27 14:07:25
    精品人妻偷拍中文字幕| 国产一区二区三区av在线| 久久精品aⅴ一区二区三区四区 | 欧美成人午夜免费资源| 91在线精品国自产拍蜜月| 国产欧美日韩综合在线一区二区| 国产精品久久久久久久久免| 中文字幕最新亚洲高清| 午夜日本视频在线| 精品人妻在线不人妻| 亚洲精品国产av成人精品| 久久这里只有精品19| 婷婷色麻豆天堂久久| 久久精品人人爽人人爽视色| 美女中出高潮动态图| 纯流量卡能插随身wifi吗| 男人爽女人下面视频在线观看| 黄片播放在线免费| 国产 一区精品| 日韩熟女老妇一区二区性免费视频| 国产日韩欧美亚洲二区| 丝袜脚勾引网站| 一本大道久久a久久精品| 搡老乐熟女国产| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩一区二区视频在线观看视频在线| 欧美人与善性xxx| 国产伦理片在线播放av一区| 欧美97在线视频| 飞空精品影院首页| 9191精品国产免费久久| 黑人欧美特级aaaaaa片| av卡一久久| 中文字幕av电影在线播放| 极品人妻少妇av视频| 18在线观看网站| 欧美日韩一区二区视频在线观看视频在线| 国产一区有黄有色的免费视频| 国产欧美日韩一区二区三区在线| 欧美xxxx性猛交bbbb| 亚洲高清免费不卡视频| 美女内射精品一级片tv| 日本av手机在线免费观看| 精品久久国产蜜桃| 精品国产一区二区久久| 亚洲美女黄色视频免费看| 视频在线观看一区二区三区| 久久女婷五月综合色啪小说| 秋霞伦理黄片| 精品亚洲成a人片在线观看| 欧美成人午夜精品| 97精品久久久久久久久久精品| 久久99热这里只频精品6学生| 少妇的丰满在线观看| 只有这里有精品99| 女人久久www免费人成看片| 亚洲丝袜综合中文字幕| 欧美 日韩 精品 国产| 久久久久精品性色| 亚洲 欧美一区二区三区| 日本猛色少妇xxxxx猛交久久| 日韩不卡一区二区三区视频在线| 男女边摸边吃奶| 最近中文字幕2019免费版| 免费看不卡的av| 日本91视频免费播放| 韩国高清视频一区二区三区| 午夜福利影视在线免费观看| 久久热在线av| 少妇的丰满在线观看| 中文字幕人妻熟女乱码| 亚洲精品视频女| 国产熟女欧美一区二区| a级片在线免费高清观看视频| 一个人免费看片子| 黄色 视频免费看| 久久久久国产精品人妻一区二区| 大片免费播放器 马上看| 国产精品一区二区在线不卡| 97在线视频观看| 草草在线视频免费看| 午夜日本视频在线| 狠狠婷婷综合久久久久久88av| 天天躁夜夜躁狠狠躁躁| 黑人巨大精品欧美一区二区蜜桃 | 一区二区日韩欧美中文字幕 | 在线精品无人区一区二区三| a级片在线免费高清观看视频| 欧美97在线视频| 国产又爽黄色视频| 国产色爽女视频免费观看| 久久精品国产鲁丝片午夜精品| 欧美成人精品欧美一级黄| 狂野欧美激情性bbbbbb| 国产精品99久久99久久久不卡 | 久久精品国产亚洲av涩爱| 久久婷婷青草| 国产色爽女视频免费观看| 国产深夜福利视频在线观看| 久久精品久久久久久久性| 曰老女人黄片| 日韩欧美精品免费久久| 久久国产亚洲av麻豆专区| 视频中文字幕在线观看| 丰满乱子伦码专区| 欧美精品国产亚洲| 亚洲国产精品专区欧美| 日本午夜av视频| 久久99精品国语久久久| 在线看a的网站| 99视频精品全部免费 在线| 欧美激情 高清一区二区三区| 考比视频在线观看| 欧美精品高潮呻吟av久久| 极品少妇高潮喷水抽搐| 日韩一区二区视频免费看| 免费人成在线观看视频色| 少妇精品久久久久久久| 国产高清三级在线| 自线自在国产av| 久久精品久久久久久久性| 久久精品夜色国产| 妹子高潮喷水视频| 亚洲伊人色综图| 亚洲少妇的诱惑av| 熟妇人妻不卡中文字幕| 九九爱精品视频在线观看| 女人被躁到高潮嗷嗷叫费观| 麻豆乱淫一区二区| 啦啦啦啦在线视频资源| 久热久热在线精品观看| 成人二区视频| 999精品在线视频| 免费观看性生交大片5| 丰满乱子伦码专区| 精品午夜福利在线看| 免费在线观看完整版高清| 亚洲成色77777| 午夜免费男女啪啪视频观看| 国产淫语在线视频| 乱码一卡2卡4卡精品| 下体分泌物呈黄色| 亚洲综合精品二区| videossex国产| 18禁在线无遮挡免费观看视频| 99九九在线精品视频| av在线播放精品| 久久 成人 亚洲| 高清不卡的av网站| 国产精品久久久久久av不卡| 好男人视频免费观看在线| 国产成人a∨麻豆精品| 亚洲婷婷狠狠爱综合网| 久久精品国产综合久久久 | 人妻 亚洲 视频| 午夜免费观看性视频| 中国美白少妇内射xxxbb| 久久久久精品人妻al黑| 久久精品人人爽人人爽视色| 美女脱内裤让男人舔精品视频| 中文字幕人妻丝袜制服| 成人毛片60女人毛片免费| 永久免费av网站大全| 日韩一本色道免费dvd| 国产极品天堂在线| 91aial.com中文字幕在线观看| 国产黄色免费在线视频| 一级片免费观看大全| 日产精品乱码卡一卡2卡三| 爱豆传媒免费全集在线观看| 人人妻人人添人人爽欧美一区卜| 99九九在线精品视频| 亚洲欧美一区二区三区国产| 久热久热在线精品观看| 国产精品蜜桃在线观看| 97在线人人人人妻| h视频一区二区三区| 伦精品一区二区三区| 免费大片黄手机在线观看| 亚洲经典国产精华液单| 免费观看在线日韩| freevideosex欧美| 国产黄频视频在线观看| 少妇被粗大猛烈的视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av.av天堂| 国产激情久久老熟女| 久久99热这里只频精品6学生| 久久 成人 亚洲| 欧美人与性动交α欧美软件 | 亚洲成av片中文字幕在线观看 | 少妇人妻精品综合一区二区| 侵犯人妻中文字幕一二三四区| 最近的中文字幕免费完整| 欧美丝袜亚洲另类| 丝瓜视频免费看黄片| 又黄又爽又刺激的免费视频.| 国产淫语在线视频| 欧美+日韩+精品| 日本-黄色视频高清免费观看| 久久免费观看电影| 超色免费av| 十分钟在线观看高清视频www| 国产片内射在线| 国国产精品蜜臀av免费| 中文字幕av电影在线播放| 丝袜美足系列| 久久国内精品自在自线图片| a级毛片黄视频| 免费看av在线观看网站| 久久久国产欧美日韩av| 日本91视频免费播放| 丝袜喷水一区| kizo精华| 赤兔流量卡办理| 91精品三级在线观看| 满18在线观看网站| 这个男人来自地球电影免费观看 | 亚洲国产最新在线播放| 久久精品夜色国产| 中文字幕免费在线视频6| 亚洲欧美色中文字幕在线| 亚洲,欧美,日韩| 欧美xxⅹ黑人| 国产精品无大码| 99久久中文字幕三级久久日本| 有码 亚洲区| 国精品久久久久久国模美| 在线 av 中文字幕| 飞空精品影院首页| 亚洲国产毛片av蜜桃av| 欧美+日韩+精品| 国产精品久久久久久久电影| 中文字幕人妻丝袜制服| 中文字幕人妻熟女乱码| 日日撸夜夜添| 国产亚洲午夜精品一区二区久久| 亚洲国产最新在线播放| 九草在线视频观看| 性高湖久久久久久久久免费观看| 香蕉精品网在线| 五月开心婷婷网| 亚洲av欧美aⅴ国产| 国产成人精品婷婷| 2021少妇久久久久久久久久久| 精品99又大又爽又粗少妇毛片| 看免费成人av毛片| 久久这里只有精品19| 国产av国产精品国产| 欧美精品一区二区免费开放| 亚洲av国产av综合av卡| 亚洲av福利一区| 99九九在线精品视频| 国产成人欧美| 久久这里有精品视频免费| 女性生殖器流出的白浆| 欧美精品一区二区免费开放| 国产国语露脸激情在线看| 国产男女内射视频| 美国免费a级毛片| 久久久久精品性色| 黄色 视频免费看| 美女中出高潮动态图| 国产成人精品福利久久| 国产免费一级a男人的天堂| 国产熟女午夜一区二区三区| 成人国语在线视频| 人人澡人人妻人| 午夜福利影视在线免费观看| 国产一级毛片在线| 国产色婷婷99| 精品熟女少妇av免费看| 免费播放大片免费观看视频在线观看| 9色porny在线观看| 2021少妇久久久久久久久久久| 日日啪夜夜爽| 十八禁网站网址无遮挡| 亚洲av中文av极速乱| 免费久久久久久久精品成人欧美视频 | 高清毛片免费看| 宅男免费午夜| 国产免费一级a男人的天堂| av福利片在线| 中国国产av一级| 日本欧美国产在线视频| 中文字幕免费在线视频6| 亚洲人成网站在线观看播放| 制服丝袜香蕉在线| 女人精品久久久久毛片| 国产1区2区3区精品| 国产精品一区www在线观看| 97在线人人人人妻| 久久国产精品大桥未久av| 欧美人与性动交α欧美精品济南到 | 国产乱人偷精品视频| 尾随美女入室| 亚洲精品aⅴ在线观看| 国产精品无大码| 亚洲情色 制服丝袜| 人妻少妇偷人精品九色| 国产深夜福利视频在线观看| 国产精品国产av在线观看| 国产白丝娇喘喷水9色精品| 制服诱惑二区| 高清av免费在线| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩另类电影网站| 国产深夜福利视频在线观看| 欧美变态另类bdsm刘玥| 22中文网久久字幕| 菩萨蛮人人尽说江南好唐韦庄| 久久影院123| 亚洲四区av| 免费黄网站久久成人精品| 80岁老熟妇乱子伦牲交| 久久精品人人爽人人爽视色| 亚洲美女黄色视频免费看| 国产男女内射视频| 春色校园在线视频观看| 国产爽快片一区二区三区| 91精品三级在线观看| 免费人妻精品一区二区三区视频| 国产视频首页在线观看| 国产精品一区www在线观看| 亚洲精华国产精华液的使用体验| 国产熟女欧美一区二区| 夫妻午夜视频| 水蜜桃什么品种好| a级片在线免费高清观看视频| 22中文网久久字幕| 国产精品三级大全| 美女国产高潮福利片在线看| 久久久a久久爽久久v久久| 99精国产麻豆久久婷婷| 中文字幕人妻熟女乱码| 欧美日韩视频高清一区二区三区二| 最近2019中文字幕mv第一页| 国产免费又黄又爽又色| 国产高清不卡午夜福利| 日本午夜av视频| 亚洲色图综合在线观看| 男人添女人高潮全过程视频| 日本色播在线视频| 丝袜美足系列| 久久这里有精品视频免费| 两性夫妻黄色片 | 国产av国产精品国产| 亚洲伊人色综图| 男女边摸边吃奶| 亚洲 欧美一区二区三区| 亚洲三级黄色毛片| 久久精品国产自在天天线| 我要看黄色一级片免费的| 深夜精品福利| 精品一区二区免费观看| 妹子高潮喷水视频| 女人久久www免费人成看片| 一级爰片在线观看| 香蕉国产在线看| 亚洲国产最新在线播放| 秋霞伦理黄片| 少妇人妻久久综合中文| 久久99热6这里只有精品| 久久av网站| 日本wwww免费看| 建设人人有责人人尽责人人享有的| 少妇高潮的动态图| 亚洲一级一片aⅴ在线观看| 亚洲人与动物交配视频| 欧美人与善性xxx| 哪个播放器可以免费观看大片| 午夜免费鲁丝| 国产色婷婷99| 99热国产这里只有精品6| 9热在线视频观看99| 少妇被粗大的猛进出69影院 | 精品久久蜜臀av无| 高清欧美精品videossex| 国产色婷婷99| 久久99一区二区三区| 欧美精品国产亚洲| 最近2019中文字幕mv第一页| 日韩欧美一区视频在线观看| 伊人亚洲综合成人网| 精品第一国产精品| 欧美成人午夜精品| 亚洲精品美女久久久久99蜜臀 | 亚洲国产精品一区二区三区在线| 国产一区二区在线观看av| 精品少妇久久久久久888优播| 丰满饥渴人妻一区二区三| 少妇被粗大的猛进出69影院 | 在线观看免费日韩欧美大片| 国产国拍精品亚洲av在线观看| 夜夜骑夜夜射夜夜干| 国产欧美日韩一区二区三区在线| 精品99又大又爽又粗少妇毛片| 国产精品不卡视频一区二区| 亚洲成av片中文字幕在线观看 | 国产乱人偷精品视频| 免费大片黄手机在线观看| 哪个播放器可以免费观看大片| 亚洲国产精品一区三区| 97精品久久久久久久久久精品| 亚洲精品久久成人aⅴ小说| 亚洲av日韩在线播放| 亚洲av福利一区| 五月天丁香电影| 精品亚洲乱码少妇综合久久| 国产欧美亚洲国产| 国产一区二区在线观看日韩| 亚洲精品久久午夜乱码| 亚洲av综合色区一区| 久久亚洲国产成人精品v| 两性夫妻黄色片 | 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡 | 国产日韩欧美视频二区| 久久久欧美国产精品| 亚洲美女黄色视频免费看| 亚洲精品久久午夜乱码| 欧美成人午夜免费资源| 亚洲av电影在线进入| av片东京热男人的天堂| 18禁在线无遮挡免费观看视频| 女的被弄到高潮叫床怎么办| 成年人免费黄色播放视频| 美女主播在线视频| 亚洲精品av麻豆狂野| videos熟女内射| 成年美女黄网站色视频大全免费| 在线观看三级黄色| 下体分泌物呈黄色| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 黑人欧美特级aaaaaa片| 韩国av在线不卡| av有码第一页| 看非洲黑人一级黄片| 精品国产国语对白av| 丁香六月天网| 欧美日韩综合久久久久久| 日本wwww免费看| 三级国产精品片| 亚洲第一区二区三区不卡| 嫩草影院入口| 高清视频免费观看一区二区| 中文精品一卡2卡3卡4更新| 18禁观看日本| 国产精品国产三级专区第一集| 大片免费播放器 马上看| 国产不卡av网站在线观看| 人人妻人人澡人人看| 久久精品久久精品一区二区三区| 国产69精品久久久久777片| 精品人妻熟女毛片av久久网站| 肉色欧美久久久久久久蜜桃| 哪个播放器可以免费观看大片| 欧美最新免费一区二区三区| 婷婷色麻豆天堂久久| 久久ye,这里只有精品| 国产亚洲一区二区精品| 捣出白浆h1v1| 日韩在线高清观看一区二区三区| 男人舔女人的私密视频| 亚洲在久久综合| 男人操女人黄网站| 狂野欧美激情性xxxx在线观看| 一本大道久久a久久精品| av线在线观看网站| 99国产精品免费福利视频| 国产1区2区3区精品| 欧美老熟妇乱子伦牲交| 亚洲欧洲日产国产| 免费黄频网站在线观看国产| 精品国产乱码久久久久久小说| 狂野欧美激情性bbbbbb| 大话2 男鬼变身卡| 美国免费a级毛片| 大陆偷拍与自拍| 美女内射精品一级片tv| 中国美白少妇内射xxxbb| 少妇的逼好多水| 伊人久久国产一区二区| 欧美人与性动交α欧美精品济南到 | 亚洲婷婷狠狠爱综合网| 亚洲国产精品999| 日韩三级伦理在线观看| 女的被弄到高潮叫床怎么办| 久久99一区二区三区| 久久国产亚洲av麻豆专区| 久久久精品免费免费高清| 精品一区二区三区四区五区乱码 | 国产精品国产三级专区第一集| 亚洲成人手机| 80岁老熟妇乱子伦牲交| 久久这里只有精品19| 97精品久久久久久久久久精品| 最近中文字幕2019免费版| 丰满迷人的少妇在线观看| videos熟女内射| 一区二区三区精品91| 亚洲综合精品二区| 人人妻人人爽人人添夜夜欢视频| 久久久久久伊人网av| 韩国高清视频一区二区三区| 色网站视频免费| 欧美日韩视频高清一区二区三区二| 91成人精品电影| 日产精品乱码卡一卡2卡三| 国语对白做爰xxxⅹ性视频网站| 日韩成人伦理影院| 精品久久蜜臀av无| 十八禁高潮呻吟视频| 制服诱惑二区| 下体分泌物呈黄色| 国产熟女欧美一区二区| 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 免费在线观看黄色视频的| 久久久亚洲精品成人影院| 一二三四中文在线观看免费高清| 尾随美女入室| √禁漫天堂资源中文www| 免费黄频网站在线观看国产| 欧美bdsm另类| 秋霞在线观看毛片| 高清黄色对白视频在线免费看| 人妻人人澡人人爽人人| 免费大片黄手机在线观看| 成人影院久久| 精品一品国产午夜福利视频| 爱豆传媒免费全集在线观看| 久久精品国产亚洲av天美| 亚洲熟女精品中文字幕| 久久99热6这里只有精品| 22中文网久久字幕| 午夜福利影视在线免费观看| 伦理电影大哥的女人| 精品少妇内射三级| 捣出白浆h1v1| 久久午夜综合久久蜜桃| 免费人成在线观看视频色| 色婷婷av一区二区三区视频| 国产精品久久久久久久久免| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 欧美日韩综合久久久久久| 亚洲婷婷狠狠爱综合网| 欧美+日韩+精品| 国产视频首页在线观看| 热99久久久久精品小说推荐| 久久久国产一区二区| 国产成人精品一,二区| 亚洲一级一片aⅴ在线观看| 欧美 日韩 精品 国产| 91午夜精品亚洲一区二区三区| 久久精品久久精品一区二区三区| 熟女av电影| 美女xxoo啪啪120秒动态图| 久久精品aⅴ一区二区三区四区 | 在线观看人妻少妇| 成人国产麻豆网| 自线自在国产av| 最黄视频免费看| 久久久精品区二区三区| 日韩一本色道免费dvd| 欧美人与性动交α欧美软件 | 国产视频首页在线观看| 国产一区有黄有色的免费视频| 亚洲色图 男人天堂 中文字幕 | 女的被弄到高潮叫床怎么办| 国产欧美日韩综合在线一区二区| 久久精品国产亚洲av天美| 日韩精品免费视频一区二区三区 | 久久久久久伊人网av| 国产乱人偷精品视频| www.色视频.com| av一本久久久久| 韩国精品一区二区三区 | videossex国产| 国产在线免费精品| 亚洲 欧美一区二区三区| 十八禁网站网址无遮挡| 天堂8中文在线网| 中国美白少妇内射xxxbb| 91久久精品国产一区二区三区| 秋霞伦理黄片| 成人影院久久| 国产在线免费精品| 日韩免费高清中文字幕av| 亚洲一码二码三码区别大吗| 九色亚洲精品在线播放| 中文天堂在线官网| 少妇人妻 视频| 中文字幕免费在线视频6| 国产男人的电影天堂91| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 全区人妻精品视频| 亚洲综合精品二区| 一级,二级,三级黄色视频| 欧美 日韩 精品 国产| 国产麻豆69| 成年动漫av网址| www.熟女人妻精品国产 | 青春草亚洲视频在线观看| 亚洲伊人久久精品综合| 亚洲国产看品久久| 久久99热这里只频精品6学生| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 22中文网久久字幕| 婷婷成人精品国产| 免费播放大片免费观看视频在线观看| 夜夜爽夜夜爽视频|