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

    基于Landsat影像的千將坪滑坡植被覆蓋度變化監(jiān)測

    2023-06-09 01:41:00付珺琳徐堅趙登忠李國忠趙保成徐健
    水利水電快報 2023年3期
    關鍵詞:覆蓋度滑坡植被

    付珺琳 徐堅 趙登忠 李國忠 趙保成 徐健

    摘要:

    監(jiān)測滑坡區(qū)植被生態(tài)恢復過程,可為滑坡區(qū)生態(tài)恢復提供科學依據(jù)。以湖北省秭歸縣千將坪滑坡為研究區(qū),以多期空間分辨率為30 m的Landsat TM/OLI影像為數(shù)據(jù)源,計算歸一化植被指數(shù)(NDVI),采用像元二分模型估算2002~2020年間9個時相的植被覆蓋度數(shù)值,對滑坡區(qū)植被覆蓋度變化趨勢進行分析。結(jié)果表明:滑坡區(qū)平均植被覆蓋度2002~2020年呈現(xiàn)先下降后上升的趨勢,在2003年發(fā)生滑坡后平均植被覆蓋度明顯下降,之后逐漸上升,至2014年平均植被覆蓋度已超過滑坡前水平?;掳l(fā)生后,低植被覆蓋區(qū)面積明顯增加,高植被覆蓋區(qū)毀損嚴重,經(jīng)歷近20 a生態(tài)恢復,至2020年研究區(qū)低植被覆蓋區(qū)面積比例低于滑坡前,高植被覆蓋區(qū)面積比例高于滑坡前,千將坪滑坡區(qū)植被生態(tài)總體恢復良好。

    關鍵詞:

    植被覆蓋度; Landsat; 動態(tài)監(jiān)測; 像元二分模型; 千將坪滑坡; 三峽庫區(qū)

    中圖法分類號:X171.4

    文獻標志碼:A

    DOI:10.15974/j.cnki.slsdkb.2023.03.018

    文章編號:1006-0081(2023)03-0101-06

    0 引 言

    植被是生態(tài)環(huán)境的指示因子,通過植被信息監(jiān)測,可以直觀地展示環(huán)境變化情況[1]。植被覆蓋度是植被垂直投影面積占地面總面積的比例,在生態(tài)環(huán)境監(jiān)測、氣候變化、土地利用、大氣污染等多個領域有廣泛應用。滑坡的發(fā)生必然會改變植被覆蓋度,遙感技術能做到大范圍、長時間且高效的監(jiān)測,可規(guī)避傳統(tǒng)方法精確性低以及耗費大量人力物力等缺點,目前已經(jīng)成為植被覆蓋度監(jiān)測的重要方法。

    本文以Landsat TM/OLI多時相影像作為數(shù)據(jù)源,結(jié)合像元二分模型,對湖北省秭歸縣千將坪滑坡2002~2020年滑坡發(fā)生前后陸地植被覆蓋度進行長時序遙感監(jiān)測,分析植被覆蓋度變化規(guī)律,探究滑坡區(qū)域植被生態(tài)恢復狀況,為滑坡區(qū)的生態(tài)恢復提供科學依據(jù)和數(shù)據(jù)支撐。

    1 研究區(qū)概況

    千將坪滑坡發(fā)生于2003年7月13日00∶20,地處湖北省秭歸縣沙鎮(zhèn)溪鎮(zhèn)千將坪村,E110°35′58″~110°36′32″,N30°58′03″~30°58′31″之間[2]?;戮嗲喔珊雍涌诩s3 km,距三峽大壩約41 km,所處斜坡向東南傾斜,坡度自上而下為 35°~15°。千將坪滑坡導致14人遇難10人失蹤,青干河河道淤積斷航,334省道宜巴公路長2.2 km的路段被毀,房屋倒塌量數(shù)以百計,22艘停泊在碼頭的漁船被毀,山林植被被局部破壞,經(jīng)濟損失高達5 000余萬元。

    千將坪滑坡地處秭歸向斜盆地西南邊緣,巖性為黃綠色泥巖砂巖互層[3-4]。千將坪滑坡平面似圈椅形,前緣高程約100 m,后緣高程約405 m,寬度410~480 m,總面積近68萬m2[5]。該地區(qū)屬亞熱帶大陸性季風氣候,氣候溫暖濕潤,四季分明,日照充足,雨水充沛,雨季集中在5~10月,全年總降水天數(shù)為120~140 d,年平均降水量為1 025 mm,年平均氣溫為18 ℃[6]。

    千將坪滑坡是由三峽水庫一期蓄水和強降雨共同作用導致的[7]。該斜坡地區(qū)滑體整體下滑后,勢能降低,部分應力釋放,因此滑坡現(xiàn)今整體穩(wěn)定性良好,即使三峽水庫水位每年175 m→145 m→175 m周期性變化,10多年來千將坪滑坡體依然保持穩(wěn)定,再次發(fā)生失穩(wěn)破壞的可能性不大。

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)選擇及預處理

    2.1.1 數(shù)據(jù)選擇

    Landsat遙感影像因獲取簡單、分辨率較高和在軌運行時間長等特點,在陸地植被覆蓋度提取方面被廣泛應用。本文選用滑坡前(2002年)、滑坡年(2003年)、恢復初期(2004,2005年)、恢復中后期(2008,2011,2014,2017,2020年)植被生長旺盛時期且云覆蓋量少的9景遙感影像作為研究數(shù)據(jù)源,如表1所示。數(shù)據(jù)來源自地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/)。

    2.1.2 數(shù)據(jù)預處理

    用Arcgis10.7和ENVI 5.3軟件對影像進行預處理。先進行輻射定標,通過輻射定標處理把遙感影像的亮度灰度值轉(zhuǎn)換為絕對輻射亮度值,以減少不同傳感器、不同時間上獲取影像的誤差。再進行Flaash大氣校正,消除由大氣條件、太陽方位和傾角而引起的輻射誤差,獲取地物真實反射率數(shù)據(jù)。最后以千將坪滑坡區(qū)矢量邊界為掩膜,裁剪經(jīng)過預處理后的遙感影像,進而獲得千將坪滑坡區(qū)植被覆蓋影像。

    2.2 歸一化植被指數(shù)計算

    歸一化植被指數(shù)(NDVI)是植被生長狀態(tài)及植被覆蓋度最佳指示因子,能消除部分與太陽高度角、衛(wèi)星觀測角、地形、云影等與大氣條件有關的輻射變化的影響[8]。定義為近紅外波段(NIR)和可見光紅光波段(RED)的反射率數(shù)值之差與反射率數(shù)值之和的比值[9]:

    NDVI=NIR-RED/NIR+RED(1)

    式中:NIR,RED依次對應近紅外波段反射率和紅光波段反射率。在Landsat TM影像中,紅光波段和近紅外波段是第3和第4波段;在Landsat OLI影像中,紅光波段和近紅外波段是第4和第5波段。

    2.3 像元二分模型估算植被覆蓋度

    像元二分模型是一種常用的線性模型,因其計算簡單結(jié)果可靠,廣泛應用于估算植被覆蓋度[10-11]。其基本原理為假定每一個像元光譜信息都由植被信息和非植被(裸土)信息組成,假定衛(wèi)星傳感器獲得的像元信息由S表示,S由植被信息Sv和裸土信息Ss組成,即:

    S=Sv+Ss(2)

    假設一個像元中有植被覆蓋的面積所占比例為FVC,則可認為此像元的植被覆蓋度是FVC,裸土覆蓋度則是1-FVC;設純植被覆蓋像元的遙感信息是Sveg,純裸土覆蓋像元的遙感信息是Ssoil[12-13],可得:

    Sv=FVC×Sveg(3)

    Ss=1-FVC×Ssoil(4)

    把公式(3)和(4)與(2)合并,轉(zhuǎn)換后可得出:

    FVC=S-Ssoil/Sveg-Ssoil(5)

    在實際應用中,把NDVI值直接代入式(5),得到用NDVI估算植被覆蓋度的公式[14]:

    FVC=NDVI-NDVIsoil/NDVIveg-NDVIsoil(6)

    式中:FVC取值范圍通常為 [0,1];NDVIsoil為無植被覆蓋區(qū)域的NDVI值;NDVIveg為純植被覆蓋區(qū)域的NDVI值。

    NDVIveg理論上應該為1,但由于空間分布不均和植被類型差異等影響因素的存在,該值往往隨時間空間變化;NDVIsoil理論上應接近于0,但由于受到大氣、土壤類型、土壤顏色、地表濕度及粗糙度等影響,其值也會發(fā)生相應的變化,一般為-0.1~0.2[15-17]。因此,為減小誤差,去掉極端值,給NDVIsoil和NDVIveg設立一個置信區(qū)間,本文設定置信度區(qū)間為5%~95%,區(qū)間最小值為NDVIsoil,區(qū)間最大值為NDVIveg。以此計算可得出千將坪滑坡區(qū)不同時期的植被覆蓋度。

    3 結(jié)果與分析

    3.1 平均植被覆蓋度變化分析

    通過對滑坡區(qū)植被覆蓋度的提取,從定量分析角度分別計算出2002~2020年9個時相的年平均植被覆蓋度(圖1)。總體上滑坡區(qū)平均植被覆蓋度從2002~2020年近20 a間呈現(xiàn)先下降后上升的趨勢。2002年平均植被覆蓋度為0.827 8,說明滑坡前研究區(qū)植被覆蓋情況良好。2003年滑坡后平均植被覆蓋度下降至0.715 1,降幅為0.112 7,從而可推斷滑坡地質(zhì)災害使該地區(qū)植被破壞嚴重。2003 ~2005年,平均植被覆蓋度快速升高,由0.715 1增至0.802 4,逐漸接近滑坡前的水平。2005~2020年,平均植被覆蓋度總體呈穩(wěn)中有升趨勢,至2014年平均植被覆蓋度達到0.836 5,已超過滑坡前水平,2020年平均植被覆蓋度達到最高值0.853 8。

    3.2 不同等級植被覆蓋度變化分析

    為更好分析植被覆蓋度變化特征,對照SL 190-2007《土壤侵蝕分類分級標準》 ,根據(jù)千將坪滑坡區(qū)實際情況將植被覆蓋度分為5級:低植被覆蓋度(0~0.30)、中低植被覆蓋度(0.30~0.45)、中植被覆蓋度(0.45~0.60)、中高植被覆蓋度(0.60~0.75)、高植被覆蓋度(0.75~1.00)。ENVI5.3和Arcgis10.7軟件生成的分級結(jié)果如圖2所示。

    統(tǒng)計2002~2020年千將坪滑坡區(qū)各等級植被覆蓋度面積和所占比例,如表2所示,并繪制了各等級植被覆蓋度面積所占比例的變化趨勢折線圖,如圖3所示。

    由上述圖表可知,滑坡區(qū)2002年中高、高植被覆蓋區(qū)域共占91.54%,低、中低植被覆蓋區(qū)域僅占3.36%。與2002年滑坡前相比,2003年滑坡后低、中低、中植被覆蓋區(qū)域占比明顯升高,面積比例依次從0.94%升至4.43%,2.42%升至15.57%,5.10%升至11.95%;中高植被覆蓋區(qū)域基本保持不變;高植被覆蓋區(qū)域占比顯著下降,由77.18%降至53.42%。由此可見,千將坪滑坡發(fā)生后,低植被覆蓋度面積大幅增加,高植被覆蓋區(qū)毀損嚴重,滑坡對植被覆蓋造成顯著破壞。

    (1) 2003~2005年,作為恢復初期,植被恢復速度較快,逐漸向滑坡前水平發(fā)展。表現(xiàn)為低、中低植被覆蓋區(qū)面積比例明顯下降,低植被覆蓋區(qū)面積占比從最大值4.43%降至0.27%,中低植被覆蓋區(qū)面積占比從最大值15.57%降至5.50%,中高、高植被覆蓋區(qū)面積比例明顯上升,中高植被覆蓋區(qū)面積占比從14.63%上升到17.18%,高植被覆蓋區(qū)面積占比從最小值53.42%上升到66.98%。

    (2) 2005~2011年,各等級植被覆蓋區(qū)面積變化幅度均低于2003~2005年滑坡恢復初期,恢復速度降低。中低植被覆蓋區(qū)域占比由5.50%降低到2.02%,中植被覆蓋度區(qū)域占比由10.07%降低到7.92%,中高植被覆蓋度區(qū)域占比由17.18%上升到22.01%,高植被覆蓋度區(qū)域面積比例波動上升,由66.98%上升到67.38%,但高植被覆蓋度區(qū)域面積始終小于滑坡前水平。

    (3) 2011~2020年,滑坡區(qū)植被覆蓋恢復較為平穩(wěn),穩(wěn)中有升,表現(xiàn)為高植被覆蓋區(qū)占比大幅上升,由67.38%升至82.01%,而其他等級植被覆蓋區(qū)占比逐步降低,至2014年,各等級植被覆蓋區(qū)面積占比與滑坡前基本相同。至2020年,低、中低、中、中高植被覆蓋區(qū)面積均小于滑坡前,高植被覆蓋區(qū)面積超過滑坡前,植被生態(tài)恢復良好。

    3.3 原因分析

    從現(xiàn)有的研究結(jié)果可知,在恢復過程中植被由早期的快速增長到后期逐步穩(wěn)定。千將坪滑坡發(fā)生后,5 a來整體植被覆蓋度快速接近滑坡前水平,但滑坡體上部依然為中低植被覆蓋度,說明植被大多都能自然恢復,但不同區(qū)域的恢復速度和恢復程度會存在一定差異。2008~2020年,整體植被覆蓋度穩(wěn)定增長,逐漸超過滑坡前,滑坡體上部植被損壞最嚴重的地方也基本恢復至滑坡前水平,說明此地區(qū)未受到明顯人為干擾,無大型工程建設或開發(fā)等破壞環(huán)境的現(xiàn)象。隨著植被的自然恢復,滑坡體生態(tài)環(huán)境也得到改善。

    4 結(jié) 論

    本文采用像元二分模型獲取了千將坪滑坡區(qū)2002~2020年間9個時相的植被覆蓋度,并將植被覆蓋度劃分成5個等級。據(jù)此分析了滑坡區(qū)年平均植被覆蓋度變化情況、不同等級植被覆蓋區(qū)面積變化情況,得出結(jié)論如下。

    (1) 通過分析千將坪滑坡區(qū)平均植被覆蓋度的變化可知,總體上滑坡區(qū)平均植被覆蓋度從2002~2020年呈現(xiàn)先下降后上升的趨勢。2002年植被平均覆蓋度為0.827 8,2003年滑坡后下降至0.715 1,說明滑坡對植被覆蓋破壞顯著,多年生態(tài)恢復后,至2014年平均植被覆蓋度達到0.836 5,超過滑坡前。2020年平均植被覆蓋度為0.853 8,是歷史最高。

    (2) 統(tǒng)計分析研究區(qū)不同等級植被覆蓋度的變化可知,與2002年滑坡前相比,2003年滑坡發(fā)生后高植被覆蓋區(qū)面積比例由77.18%下降到53.42%,低、中低植被覆蓋區(qū)面積合計占比由3.36%上升到20%,說明滑坡導致低植被覆蓋區(qū)面積顯著增加,高植被覆蓋區(qū)面積大幅減少。經(jīng)過近20 a漫長的生態(tài)重建,至2020年滑坡區(qū)低、中低、中、中高等級植被覆蓋區(qū)面積比例均小于滑坡前,高植被覆蓋區(qū)面積比例大于滑坡前值。低、中低、中植被覆蓋區(qū)面積共降低了11 784.48 m2,高植被覆蓋區(qū)面積增長了32 525.17 m2,植被恢復良好。

    (3) 千將坪滑坡植被穩(wěn)定增長的主要原因是植被的自然恢復,滑坡后無人為干擾,無大型工程建設或開發(fā)等破壞環(huán)境的現(xiàn)象。

    參考文獻:

    [1] 張嵩,馬保東,陳玉騰,等.超貧磁鐵礦區(qū)植被覆蓋度變化的遙感監(jiān)測[J].測繪與空間地理信息,2018,41(6):126-128,131.

    [2] 楊日紅,楊金中,王治華.千將坪滑坡特征信息遙感提取及滑坡成因機理分析[J].地球與環(huán)境,2007(1):85-90.

    [3] 侯珍珠,胡志強,肖詩榮,等.三峽庫區(qū)秭歸縣千將坪滑坡軟化模型實驗研究[J].三峽大學學報(自然科學版),2019,41(3):36-40.

    [4] 肖詩榮,劉德富,姜福興,等.三峽庫區(qū)千將坪滑坡地質(zhì)力學模型試驗研究[J].巖石力學與工程學報,2010,29(5):1023-1030.

    [5] 李會中,潘玉珍,王團樂,等.三峽庫區(qū)千將坪滑坡成因與機制分析[J].人民長江,2006(7):12-14,20.

    [6] 李會中,王團樂,孫立華,等.三峽庫區(qū)千將坪滑坡地質(zhì)特征與成因機制分析[J].巖土力學,2006,27(增2):1239-1244.

    [7] 張振華,錢明明,位偉.基于改進破壞接近度的千將坪岸坡失穩(wěn)機制分析[J].巖石力學與工程學報,2018,37(6):1371-1384.

    [8] 曾強國,高子讓,丁鋒.基于Landsat8 OLI數(shù)據(jù)的植被覆蓋度時空變化分析[J].國土資源導刊,2021,18(3):30-34,91.

    [9] 苗正紅,劉志明,王宗明,等.基于MODIS NDVI的吉林省植被覆蓋度動態(tài)遙感監(jiān)測[J].遙感技術與應用,2010,25(3):387-393.

    [10] 張珍珍,熊康寧,黃登紅.近30年來梵凈山植被覆蓋時空變化及影響因素分析[J].水土保持研究,2018,25(2):183-189,389.

    [11] 李苗苗.植被覆蓋度的遙感估算方法研究[D].北京:中國科學院(遙感應用研究所),2003.

    [12] 彭文甫,周介銘,徐新良,等.成都平原及其周邊區(qū)域植被覆蓋動態(tài)監(jiān)測[J].地球與環(huán)境,2017,45(2):193-202.

    [13] 溫小樂,姜興強,徐涵秋.基于TM影像的長樂市植被覆蓋度變化研究[J].遙感信息,2013,28(6):56-61.

    [14] 彭文甫,王廣杰,周介銘,等.基于多時相Landsat5/8影像的岷江汶川-都江堰段植被覆蓋動態(tài)監(jiān)測[J].生態(tài)學報,2016,36(7):1975-1988.

    [15] 陳峰.基于Landsat TM/OIL的濟寧市植被覆蓋度動態(tài)變化分析[J].河南科技,2021,40(22):129-133.

    [16] 吳昌廣,周志翔,肖文發(fā),等.基于MODIS NDVI的三峽庫區(qū)植被覆蓋度動態(tài)監(jiān)測[J].林業(yè)科學,2012,48(1):22-28.

    [17] 楊磊,張梅,羅明良,等.基于MODIS NDVI的川中丘陵區(qū)植被覆蓋度景觀格局變化[J].生態(tài)學雜志,2013,32(1):171-177.

    (編輯:李 晗)

    Monitoring vegetation cover changes on Qianjiangping landslide based on Landsat images

    FU Junlin1,XU Jian1,ZHAO Dengzhong2,LI Guozhong1,ZHAO Baocheng1,XU Jian1

    (1.Institute of Spatial Information Technology Application,Changjiang River Scientific Research Institute,Wuhan 430010,China; 2.Science and Technology Exchange and International Cooperation Division,Changjiang River Scientific Research Institute,Wuhan 430010,China)Abstract:

    In order to monitor the ecological recovery process of vegetation in landslide areas,and provide a scientific basis for the ecological recovery of landslide areas,Qianjiangping landslide in Zigui county,Hubei province was taken as the research area for study.The Landsat TM/OLI imagery with a spatial resolution of 30m was used as the data source to calculate the normalized vegetation index (NDVI),and an image dichotomous model was used to estimate the vegetation cover values for 9 time phases from 2002 to 2020,and to analyze the change of vegetation cover of this landslide area.The results showed that the average vegetation cover in the landslide area from 2002 to 2020 showed a trend of first decreasing and then increasing,with the average vegetation cover decreasing significantly after the landslide occurred in 2003,and then gradually increasing,until the average vegetation cover exceeded the level before the landslide in 2014.After the occurrence of the landslide,the area of low vegetation cover increased significantly,while the area of high vegetation cover was severely damaged.After nearly 20 years of ecological recovery,by 2020,the proportion of the area of low vegetation cover in the study area was lower than that before the landslide,while the proportion of the area of high vegetation cover was higher than that before the landslide.A generally recovery was shown in Qinajiangping landslide area.

    Key words:

    vegetation cover; Landsat; dynamic monitoring; image element dichotomous model; Qianjiangping landslide; Three Gorges Reservoir area

    收稿日期:

    2022-05-10

    基金項目:

    湖北省自然資源廳科研計劃項目“三峽庫區(qū)典型水土流失與滑坡區(qū)生態(tài)修復關鍵技術與應用示范”(ZRZY2021KJ10)

    作者簡介:

    付珺琳,女,碩士,主要從事遙感、地理信息系統(tǒng)研究。E-mail:fujunlin12@163.com

    猜你喜歡
    覆蓋度滑坡植被
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
    基于植被復綠技術的孔植試驗及應用
    基于NDVI的晉州市植被覆蓋信息提取
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    低覆蓋度CO分子在Ni(110)面的吸附研究
    綠色植被在溯溪旅游中的應用
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    基于Fluent的滑坡入水過程數(shù)值模擬
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    基于原生植被的長山群島植被退化分析
    午夜福利视频在线观看免费| 精品一品国产午夜福利视频| 97在线人人人人妻| 精品亚洲成国产av| 97超碰精品成人国产| 精品国产国语对白av| 欧美精品一区二区大全| 亚洲精品aⅴ在线观看| 一级片免费观看大全| 母亲3免费完整高清在线观看 | 日韩欧美一区视频在线观看| 国产精品久久久av美女十八| 国产av一区二区精品久久| 国产成人一区二区在线| 精品酒店卫生间| 国产国拍精品亚洲av在线观看| 久久热在线av| 黄色 视频免费看| 尾随美女入室| av在线观看视频网站免费| 夫妻午夜视频| a级毛色黄片| 热99久久久久精品小说推荐| av播播在线观看一区| 中文字幕最新亚洲高清| 超碰97精品在线观看| kizo精华| 黄片无遮挡物在线观看| 一区二区三区四区激情视频| 亚洲美女搞黄在线观看| 99视频精品全部免费 在线| 少妇人妻 视频| 高清不卡的av网站| av在线app专区| 国产精品一区二区在线观看99| 美女主播在线视频| 王馨瑶露胸无遮挡在线观看| 欧美成人午夜免费资源| 99国产综合亚洲精品| 久久这里只有精品19| 啦啦啦视频在线资源免费观看| 黄色 视频免费看| 99国产精品免费福利视频| 免费高清在线观看视频在线观看| 午夜久久久在线观看| 国产亚洲最大av| 在现免费观看毛片| 欧美成人精品欧美一级黄| 18禁国产床啪视频网站| 国产极品天堂在线| 一级毛片 在线播放| 黄片播放在线免费| 美女国产高潮福利片在线看| 内地一区二区视频在线| 在线观看一区二区三区激情| 18+在线观看网站| 成人亚洲精品一区在线观看| 女人久久www免费人成看片| 国产欧美亚洲国产| 免费观看在线日韩| 日本与韩国留学比较| 久久久久久伊人网av| 成人影院久久| 丝袜美足系列| 国产黄色视频一区二区在线观看| 中文字幕免费在线视频6| 久久狼人影院| 日韩熟女老妇一区二区性免费视频| 精品福利永久在线观看| 国产成人免费观看mmmm| 曰老女人黄片| 午夜视频国产福利| 欧美亚洲 丝袜 人妻 在线| 青春草亚洲视频在线观看| 国产日韩欧美视频二区| a级片在线免费高清观看视频| 国产精品国产三级专区第一集| 国产视频首页在线观看| 青春草国产在线视频| 高清毛片免费看| 欧美97在线视频| 熟女av电影| 老司机影院毛片| 国产精品无大码| 久久这里只有精品19| 国产一区二区在线观看av| 99久久中文字幕三级久久日本| 伦理电影免费视频| 国语对白做爰xxxⅹ性视频网站| 天堂俺去俺来也www色官网| 久久久国产欧美日韩av| 尾随美女入室| 成人毛片60女人毛片免费| 国产欧美日韩综合在线一区二区| 男女边吃奶边做爰视频| 性色avwww在线观看| 成人二区视频| 成年动漫av网址| 插逼视频在线观看| 少妇被粗大猛烈的视频| 国产成人精品婷婷| 三级国产精品片| 亚洲久久久国产精品| 99香蕉大伊视频| 精品少妇内射三级| 国产高清三级在线| 汤姆久久久久久久影院中文字幕| 日日爽夜夜爽网站| 国产亚洲av片在线观看秒播厂| 欧美日韩一区二区视频在线观看视频在线| 自拍欧美九色日韩亚洲蝌蚪91| av女优亚洲男人天堂| 亚洲人成网站在线观看播放| 精品国产乱码久久久久久小说| 美女内射精品一级片tv| 成人午夜精彩视频在线观看| av福利片在线| 日本欧美国产在线视频| 亚洲,欧美,日韩| av福利片在线| 一本久久精品| 久久婷婷青草| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区| 在线观看免费视频网站a站| 男女免费视频国产| 欧美日韩视频高清一区二区三区二| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 日韩三级伦理在线观看| 在线免费观看不下载黄p国产| 色5月婷婷丁香| 日韩伦理黄色片| 精品久久久精品久久久| 97超碰精品成人国产| 99香蕉大伊视频| 日韩欧美精品免费久久| 赤兔流量卡办理| 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频| av在线老鸭窝| 欧美日韩成人在线一区二区| 午夜福利影视在线免费观看| av视频免费观看在线观看| 最新的欧美精品一区二区| 日本av手机在线免费观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲四区av| 亚洲精品久久成人aⅴ小说| 日本av免费视频播放| 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| 国产精品秋霞免费鲁丝片| 99久久精品国产国产毛片| 国产成人午夜福利电影在线观看| 精品一区二区免费观看| 秋霞在线观看毛片| 精品国产一区二区三区久久久樱花| 日本猛色少妇xxxxx猛交久久| 欧美3d第一页| 亚洲av日韩在线播放| 国产片内射在线| 999精品在线视频| 岛国毛片在线播放| 国产精品三级大全| 97在线人人人人妻| 国产一区二区三区av在线| 丰满少妇做爰视频| 尾随美女入室| 777米奇影视久久| 天堂俺去俺来也www色官网| 国产成人精品在线电影| 欧美性感艳星| 晚上一个人看的免费电影| 中文字幕人妻丝袜制服| 天美传媒精品一区二区| 免费人成在线观看视频色| 久久99蜜桃精品久久| 精品国产国语对白av| 成人午夜精彩视频在线观看| 亚洲美女视频黄频| 天堂俺去俺来也www色官网| 国产男女超爽视频在线观看| 色哟哟·www| 国产成人欧美| 男的添女的下面高潮视频| 免费观看性生交大片5| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| 咕卡用的链子| 亚洲精华国产精华液的使用体验| 国产一区亚洲一区在线观看| 久久精品久久久久久久性| 熟女电影av网| 最近中文字幕高清免费大全6| 国产色爽女视频免费观看| 亚洲精品,欧美精品| 国产免费一区二区三区四区乱码| 少妇被粗大的猛进出69影院 | 一边摸一边做爽爽视频免费| 国产av码专区亚洲av| 国产亚洲最大av| av线在线观看网站| 久久99蜜桃精品久久| 最近手机中文字幕大全| 大香蕉久久网| 久久久久人妻精品一区果冻| 国产日韩欧美在线精品| 亚洲美女视频黄频| 亚洲人与动物交配视频| 日本-黄色视频高清免费观看| 9色porny在线观看| 男女午夜视频在线观看 | 黑人巨大精品欧美一区二区蜜桃 | 久久久久久久国产电影| 在线观看免费日韩欧美大片| 男人爽女人下面视频在线观看| 十八禁网站网址无遮挡| 黑人猛操日本美女一级片| 男人舔女人的私密视频| 久久久国产精品麻豆| 精品国产乱码久久久久久小说| 国产激情久久老熟女| 一级毛片我不卡| 国产不卡av网站在线观看| 精品一区二区三区四区五区乱码 | 午夜福利在线观看免费完整高清在| 国产成人av激情在线播放| 天天影视国产精品| 久久久久久伊人网av| 99热这里只有是精品在线观看| 国产精品久久久久久久电影| 日日爽夜夜爽网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 99热网站在线观看| 三上悠亚av全集在线观看| 中文天堂在线官网| 卡戴珊不雅视频在线播放| 热99久久久久精品小说推荐| 欧美精品国产亚洲| 97在线人人人人妻| 免费播放大片免费观看视频在线观看| 国产乱来视频区| 在线天堂最新版资源| 青青草视频在线视频观看| 天天影视国产精品| 十分钟在线观看高清视频www| 国产成人aa在线观看| 国产不卡av网站在线观看| 中文字幕亚洲精品专区| 久久精品国产亚洲av天美| 考比视频在线观看| 国产熟女欧美一区二区| 成人国产麻豆网| 看免费av毛片| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 国产成人av激情在线播放| 成年动漫av网址| 少妇人妻久久综合中文| 亚洲精品456在线播放app| 男的添女的下面高潮视频| 久久精品国产综合久久久 | 大香蕉久久成人网| 色视频在线一区二区三区| 内地一区二区视频在线| 美女脱内裤让男人舔精品视频| 久久人人爽av亚洲精品天堂| 久久久欧美国产精品| 亚洲一级一片aⅴ在线观看| 国产成人欧美| 国产精品无大码| a 毛片基地| 99国产综合亚洲精品| 日本爱情动作片www.在线观看| 日韩精品免费视频一区二区三区 | 激情视频va一区二区三区| 最近最新中文字幕大全免费视频 | 欧美精品人与动牲交sv欧美| 大香蕉97超碰在线| 九九在线视频观看精品| 嫩草影院入口| 看免费成人av毛片| 人人妻人人添人人爽欧美一区卜| 中国美白少妇内射xxxbb| 免费久久久久久久精品成人欧美视频 | 久久久精品区二区三区| 免费观看av网站的网址| 国产精品三级大全| a级片在线免费高清观看视频| 亚洲经典国产精华液单| 国产精品国产三级国产av玫瑰| 毛片一级片免费看久久久久| 丁香六月天网| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲日产国产| 9热在线视频观看99| 国产av一区二区精品久久| 欧美日韩成人在线一区二区| av天堂久久9| 晚上一个人看的免费电影| 色网站视频免费| 大码成人一级视频| a级毛片黄视频| 国产av一区二区精品久久| 成人黄色视频免费在线看| 日韩欧美精品免费久久| 国产成人精品一,二区| av播播在线观看一区| 久久韩国三级中文字幕| 夜夜爽夜夜爽视频| 国产精品无大码| 日本欧美国产在线视频| 国产成人精品久久久久久| 日韩欧美精品免费久久| 九九在线视频观看精品| 婷婷色综合www| 免费av不卡在线播放| 亚洲精品av麻豆狂野| 亚洲久久久国产精品| 久久影院123| 亚洲丝袜综合中文字幕| 夫妻午夜视频| 久久毛片免费看一区二区三区| 国产男女内射视频| 午夜福利乱码中文字幕| 久久久精品94久久精品| 国产午夜精品一二区理论片| 国产高清国产精品国产三级| 宅男免费午夜| 美女xxoo啪啪120秒动态图| 久久毛片免费看一区二区三区| 午夜久久久在线观看| 免费少妇av软件| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 亚洲伊人久久精品综合| 日本91视频免费播放| 亚洲美女视频黄频| 日韩制服丝袜自拍偷拍| 亚洲综合色惰| 国产高清三级在线| 啦啦啦在线观看免费高清www| 久久久久久久亚洲中文字幕| 亚洲av福利一区| 天天影视国产精品| 国产麻豆69| 久久人人爽人人片av| 你懂的网址亚洲精品在线观看| 制服人妻中文乱码| 久久精品国产a三级三级三级| 91午夜精品亚洲一区二区三区| 亚洲,一卡二卡三卡| 午夜精品国产一区二区电影| 国产福利在线免费观看视频| 日本av手机在线免费观看| 国产成人精品一,二区| 国产 一区精品| 大香蕉97超碰在线| 国产成人一区二区在线| 国产精品人妻久久久影院| 婷婷色av中文字幕| 女人久久www免费人成看片| 黑人高潮一二区| 香蕉精品网在线| 一级毛片黄色毛片免费观看视频| 九色成人免费人妻av| 韩国av在线不卡| 99国产精品免费福利视频| 亚洲av.av天堂| 熟女人妻精品中文字幕| 另类亚洲欧美激情| 国产欧美日韩综合在线一区二区| 欧美人与性动交α欧美软件 | 超色免费av| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 久久国内精品自在自线图片| 亚洲国产成人一精品久久久| 女性生殖器流出的白浆| 七月丁香在线播放| a 毛片基地| 国产一区二区三区综合在线观看 | 一二三四中文在线观看免费高清| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| 丰满乱子伦码专区| 亚洲精品中文字幕在线视频| 国产成人精品在线电影| 多毛熟女@视频| 香蕉精品网在线| 亚洲欧美精品自产自拍| 中文字幕最新亚洲高清| 伊人久久国产一区二区| 亚洲精品久久久久久婷婷小说| 亚洲成人手机| 久久女婷五月综合色啪小说| videossex国产| 国产av精品麻豆| 91精品伊人久久大香线蕉| 久久精品人人爽人人爽视色| videossex国产| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 不卡视频在线观看欧美| 91精品伊人久久大香线蕉| 国产在线免费精品| 精品少妇内射三级| 国产男女超爽视频在线观看| 美国免费a级毛片| 成人亚洲欧美一区二区av| 99久久综合免费| 免费大片黄手机在线观看| 国产综合精华液| 97超碰精品成人国产| 伦理电影免费视频| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看| 精品少妇黑人巨大在线播放| 亚洲 欧美一区二区三区| 欧美激情国产日韩精品一区| 天堂8中文在线网| 亚洲国产日韩一区二区| av在线播放精品| 久久精品人人爽人人爽视色| 亚洲欧美清纯卡通| 国产永久视频网站| 婷婷成人精品国产| 久久这里只有精品19| 国产高清三级在线| www.熟女人妻精品国产 | 午夜老司机福利剧场| 肉色欧美久久久久久久蜜桃| 欧美xxxx性猛交bbbb| 黄片无遮挡物在线观看| 乱人伦中国视频| 亚洲一级一片aⅴ在线观看| 日韩精品免费视频一区二区三区 | 成人综合一区亚洲| 国产视频首页在线观看| 岛国毛片在线播放| 2021少妇久久久久久久久久久| 伊人亚洲综合成人网| 天堂俺去俺来也www色官网| av福利片在线| 国产精品麻豆人妻色哟哟久久| 少妇的逼好多水| 尾随美女入室| 18禁在线无遮挡免费观看视频| 永久网站在线| 久久99热6这里只有精品| 26uuu在线亚洲综合色| 十八禁高潮呻吟视频| 高清av免费在线| 尾随美女入室| 国产免费一区二区三区四区乱码| 永久免费av网站大全| 三级国产精品片| 久久久久精品久久久久真实原创| 国产又爽黄色视频| 9色porny在线观看| 亚洲成人av在线免费| 丝袜脚勾引网站| 18禁动态无遮挡网站| 桃花免费在线播放| 久久这里有精品视频免费| 国产综合精华液| 老司机影院毛片| 曰老女人黄片| 日本午夜av视频| 精品亚洲成a人片在线观看| 高清av免费在线| 久久久久久久久久成人| 一区二区三区乱码不卡18| 纵有疾风起免费观看全集完整版| 九色成人免费人妻av| 欧美变态另类bdsm刘玥| 久久 成人 亚洲| 国产成人精品久久久久久| 三级国产精品片| 一区二区av电影网| 亚洲三级黄色毛片| 日本爱情动作片www.在线观看| 美女脱内裤让男人舔精品视频| 九色亚洲精品在线播放| 日韩伦理黄色片| 免费少妇av软件| 欧美性感艳星| 国产熟女欧美一区二区| 秋霞在线观看毛片| 极品人妻少妇av视频| 久久这里有精品视频免费| 最近中文字幕高清免费大全6| 全区人妻精品视频| 久久精品夜色国产| 视频区图区小说| 青春草视频在线免费观看| 青青草视频在线视频观看| 亚洲性久久影院| 免费大片18禁| 免费av不卡在线播放| 国产 精品1| 亚洲一区二区三区欧美精品| 在线免费观看不下载黄p国产| 一区二区三区精品91| 韩国精品一区二区三区 | 十八禁网站网址无遮挡| 国产亚洲精品久久久com| 18在线观看网站| 乱人伦中国视频| 一级毛片电影观看| 五月开心婷婷网| 国产片特级美女逼逼视频| 成人黄色视频免费在线看| 国产高清三级在线| 天堂8中文在线网| 欧美日韩一区二区视频在线观看视频在线| 建设人人有责人人尽责人人享有的| www.色视频.com| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 国产男女超爽视频在线观看| 精品久久久精品久久久| 免费观看在线日韩| 一本大道久久a久久精品| 午夜福利在线观看免费完整高清在| 亚洲国产色片| 国产精品 国内视频| 久久人人爽人人片av| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 制服诱惑二区| 欧美成人午夜精品| 91国产中文字幕| 免费黄色在线免费观看| 高清黄色对白视频在线免费看| 免费黄色在线免费观看| 国产精品成人在线| 久久人人爽人人片av| 免费高清在线观看视频在线观看| 韩国高清视频一区二区三区| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 免费久久久久久久精品成人欧美视频 | 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 男女高潮啪啪啪动态图| 黑丝袜美女国产一区| 欧美成人午夜精品| 亚洲成国产人片在线观看| 在线精品无人区一区二区三| 国产精品嫩草影院av在线观看| 欧美日韩视频精品一区| 免费少妇av软件| 久久97久久精品| 国产精品久久久久久久电影| 亚洲精品日本国产第一区| 日韩人妻精品一区2区三区| 欧美xxxx性猛交bbbb| 亚洲精品视频女| 国产 精品1| 午夜精品国产一区二区电影| 亚洲国产精品一区二区三区在线| 最近中文字幕高清免费大全6| 久久鲁丝午夜福利片| 嫩草影院入口| 熟女电影av网| 精品一区二区三区视频在线| 亚洲精品视频女| 免费看av在线观看网站| 免费在线观看黄色视频的| 涩涩av久久男人的天堂| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 在线观看美女被高潮喷水网站| 少妇高潮的动态图| 午夜老司机福利剧场| 国产福利在线免费观看视频| 亚洲欧美日韩另类电影网站| 香蕉国产在线看| 欧美 日韩 精品 国产| 国产一区二区在线观看av| 天天影视国产精品| 丰满乱子伦码专区| 国产探花极品一区二区| 久久精品国产综合久久久 | 99久久精品国产国产毛片| 男女啪啪激烈高潮av片| 日本91视频免费播放| 久久精品久久精品一区二区三区| 边亲边吃奶的免费视频| 久久精品国产鲁丝片午夜精品| 黄片播放在线免费| 亚洲精品aⅴ在线观看| 看非洲黑人一级黄片| 在线观看人妻少妇| 久久97久久精品| 亚洲av国产av综合av卡| 精品亚洲成国产av| 一区二区三区四区激情视频| 亚洲精品日韩在线中文字幕| 久久久久久久亚洲中文字幕| 久久国产亚洲av麻豆专区| 亚洲一区二区三区欧美精品| 久久久久久久亚洲中文字幕| 卡戴珊不雅视频在线播放| av在线观看视频网站免费| 精品国产一区二区三区四区第35| 91精品国产国语对白视频| 在线观看免费日韩欧美大片| 美女中出高潮动态图| 成人无遮挡网站| 欧美人与性动交α欧美精品济南到 | 午夜精品国产一区二区电影| 中文字幕最新亚洲高清|