亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于遙感監(jiān)測的2000—2020年三屯河流域植被覆蓋時空演化分析

        2023-07-03 01:45:12張艷嬌宿彥鵬陶洪飛姜有為馬合木江艾合買提魏建群楊文新
        長江科學院院報 2023年6期
        關鍵詞:屯河覆蓋度植被

        張艷嬌,李 巧,宿彥鵬,陶洪飛,姜有為,馬合木江·艾合買提,魏建群,楊文新

        (1.新疆農業(yè)大學 水利與土木工程學院,烏魯木齊 830052; 2.新疆農業(yè)大學 新疆水利工程安全與水災害防治重點實驗室,烏魯木齊 830052; 3.甘肅省臨夏回族自治州廣河縣水務局,甘肅 臨夏 731300 )

        1 研究背景

        植被是人類賴以生存發(fā)展的物質基礎[1-2],是陸地生態(tài)系統(tǒng)的重要組成部分[3],對生態(tài)系統(tǒng)中的能量流動、物質循環(huán)和信息傳遞起著先導的作用。植被與周圍環(huán)境密切相關,對生態(tài)環(huán)境的變化反應十分敏感[4],人們常根據(jù)植被生長狀況推測區(qū)域生態(tài)環(huán)境的變化情況,因而從時間及空間尺度上監(jiān)測植被的變化情況對于揭示生態(tài)環(huán)境變化、評價區(qū)域生態(tài)質量具有重要的意義[5-6]。隨著科技的進步,覆蓋范圍廣、時間有效性強且易獲取[7]的遙感監(jiān)測技術成為大范圍內植被監(jiān)測的主要手段。

        目前基于遙感反演手段的植被覆蓋度估算方法有很多,如植被指數(shù)法、像元分解法、經驗模型法等[8],其中通過提取歸一化植被指數(shù)計算植被覆蓋度最為常用[9-10]。 Neinavaz等[11]通過收集和評述利用熱紅外遙感數(shù)據(jù)獲取的葉和冠層水平陸地植被最相關的信息,討論了在熱紅外區(qū)域植被研究中利用遙感的效果,結果表明植被熱紅外遙感是其他遙感數(shù)據(jù)源的補充,融合來自光譜不同部分的數(shù)據(jù)具有很高的潛力。 Zhang等[12]利用歸一化植被指數(shù)(NDVI )數(shù)據(jù)分析2000—2019年洞庭湖濕地植被變化的時空格局和驅動因素。結果表明人類活動有利于濕地植被恢復,氣候變化威脅濕地植被。田智慧等[13]以NDVI作為植被覆蓋及生長狀況指標,研究了2000—2020年黃河流域植被時空演化驅動機制。得出研究區(qū)降水對植被的影響高于氣溫,人類活動和氣候變化對黃河流域植被變化貢獻度分別為72%和28%。

        新疆屬于典型的干旱地區(qū),對生態(tài)環(huán)境變化反應十分敏感,具有氣候干燥、降水量少且分布不均、蒸發(fā)性強、水資源匱乏、水土流失、植被退化等自然屬性及生態(tài)環(huán)境問題[14-15],因此新疆的植被覆蓋變化情況一直以來備受學者的關注[16-17]。同時,新疆作為中國面積最大的省級行政區(qū),具有很大的發(fā)展空間,三屯河流域位于新疆天山北坡經濟帶的中部,是以農業(yè)為傳統(tǒng)經濟基礎的昌吉市的經濟保證[18]。因此研究三屯河流域長時間序列的植被覆蓋變化對掌握三屯河流域環(huán)境的發(fā)展變化規(guī)律、生態(tài)環(huán)境保護及促進經濟發(fā)展具有重要意義。目前三屯河流域長時間尺度植被覆蓋變化方面的研究還處于空白階段,本研究以三屯河流域為研究區(qū),以Landsat TM/ETM+/OLI衛(wèi)星遙感數(shù)據(jù)、像元二分模型為基礎,通過圖像差值法、變異系數(shù)、轉移矩陣等手段,探討自然因素與人類因素對三屯河流域植被覆蓋的影響,分析三屯河流域植被覆蓋時空變化特征,為促進三屯河流域生態(tài)恢復、經濟發(fā)展提供科學依據(jù)。

        2 材料與方法

        2.1 研究區(qū)概況

        三屯流域發(fā)源于天山北坡中段的天格爾峰,由大氣降水和冰雪融水進行補給,冰川補給是三屯河徑流較穩(wěn)定的補給來源之一[19]。三屯河流域面積為4 466 km2,流域整體呈條帶狀,其地理位置在86°20′E—86°58′E、43°08′N—44°44′N[20]之間,是典型的內陸半干旱地區(qū)[21]。三屯河流域位于新疆天山北坡經濟帶的中部,承擔著供應昌吉市生產生活用水任務[22],是昌吉市的經濟保證。三屯河流域的東西南北4個方向上分別對應以頭屯河、呼圖壁河、阿斯克達、古爾班通古特沙漠作為邊界[23]。三屯河自發(fā)源地起始流經南部山區(qū)、中部平原,消失于北部沙漠[24],呈現(xiàn)出南高北低,由西南向東北傾斜的地勢。三屯河流域生態(tài)環(huán)境較為脆弱,整個流域自上游至下游以天然林、干草原、半荒漠、平原、荒漠的生態(tài)環(huán)境進行過渡,三屯河流域年平均降水量為193.75 mm,年平均氣溫為6.87 ℃,屬中溫帶大陸性半荒漠干旱氣候,具有降雨量小、蒸發(fā)量大等特點[25]。

        2.2 數(shù)據(jù)來源及數(shù)據(jù)處理

        綜合考慮三屯河流域的農作物[26]、自然植被生長期、沙塵暴、云量對植被覆蓋度的影響,本次研究選用2000年9月17日、2004年9月4日、2008年8月6日、2012年8月17日、2015年8月18日及2020年8月31日6期、12景分辨率為30 m的遙感影像圖,遙感數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http:∥www.gscloud.cn)及美國地質勘探局USGS(http:∥glovis.usgs.gov)。研究所用到的用于繪制三屯河流域邊界的30 m分辨率數(shù)字高程模型(Digital Elevation Model,DEM)高程數(shù)據(jù)也來源于地理空間數(shù)據(jù)云。氣象數(shù)據(jù)來自國家氣象信息中心(http:∥data.cma.cn/),采用2000年、2004年、2008年、2012年及2020年的年均氣溫及降水量數(shù)據(jù)。土地利用數(shù)據(jù)來自資源環(huán)境科學與數(shù)據(jù)中心(https:∥www.resdc.cn/),采用2000年及2020年的土地利用類型數(shù)據(jù),分辨率為1 km。

        首先在ArcGIS 10.6軟件中通過提取河網、分水嶺等操作將三屯河流域邊界圖繪制出來,然后運用ENVI 5.3軟件對本研究用到的遙感影像依次去條帶、輻射定標、大氣校正、鑲嵌及裁剪和掩膜等處理,得到三屯河流域遙感影像圖。運用ArcGIS 10.6軟件處理已下載的土地利用衛(wèi)星數(shù)據(jù)得到三屯河流域土地利用變化情況。

        2.3 研究方法

        2.3.1 NDVI計算

        NDVI是所有植被指數(shù)中應用最為廣泛的一種,其能夠涵蓋90%以上的植被信息[27-28],對植被生長、覆蓋變化等反應十分靈敏,因此常用于植被監(jiān)測。NDVI的計算公式為

        (1)

        式中:NIR表示近紅外波段的反射值;RED表示紅外波段的反射值。Landsat TM/ETM+的NIR對應的是Band4,RED對應的是Band3;Landsat OLI的NIR對應的是Band5,RED對應的是Band4。

        不同研究區(qū)域的置信區(qū)間的取值范圍也不盡相同,常見的置信區(qū)間取值范圍有1%~99%、2%~98%、5%~95%[29-32]。研究處理得到的NDVI像元數(shù)據(jù)在2%左右時有遞增的變化趨勢,故本次研究取2%~98%為置信區(qū)間,即取頻次為98%的NDVI作為NDVIv的值,取頻次為2%的NDVI作為NDVIs的值。為保證NDVI最大值與最小值的準確性,在ENVI 5.3軟件中利用波段運算掩膜去除NDVI>1與NDVI<-1的異常值。

        根據(jù)半干旱地區(qū)植被特征并結合《土地利用現(xiàn)狀調查技術規(guī)程》《草場資源調查技術規(guī)程》《土壤侵蝕分類分級標準》等標準規(guī)范,將研究區(qū)植被覆蓋度VFC劃分為4個等級:裸地或極低(0≤VFC<0.2)、低(0.2≤VFC<0.3)、中(0.3≤VFC<0.6)、高(0.6≤VFC<1)。

        2.3.2 植被覆蓋度估算

        利用像元二分模型提取植被覆蓋度是常用的植被覆蓋度評估方法,其在解譯過程中能夠減少大氣、水體、土壤背景對遙感影像的影響[33]。像元二分模型是將傳感器監(jiān)測到的光譜信息一分為二,即植被覆蓋的像元信息和沒有植被覆蓋的像元信息。像元二分模型的計算公式為

        (2)

        式中:S表示混合像元信息;Ss表示無植被覆蓋像元信息;Sv表示植被覆蓋像元信息。

        由于NDVI對植被生長狀態(tài)和植被覆蓋等信息非常靈敏,且能夠減少部分因輻射引起的誤差且其與植被覆蓋度具有強相關性[34],故在像元二分模型解析光譜波段和植被指數(shù)信息進行遙感提取時,常選取NDVI進行估測。用NDVI替換像元信息,將NDVI也分為兩部分即有植被覆蓋的NDVIv和無植被覆蓋的NDVIs,得到如下植被覆蓋度的計算公式為

        (3)

        基于研究區(qū)2000—2020年20 a間6期12景遙感影像,通過ENVI 5.3軟件對數(shù)據(jù)進行處理,分別得到2000年、2004年、2008年、2012年、2015年及2020年6組年均植被覆蓋度數(shù)據(jù)用于分析20 a間三屯河流域植被覆蓋的變化趨勢。

        2.3.3 圖像差值比較法

        圖像差值法是基于遙感影像灰度值的差值原理而產生的一種差值比較法,可以反映研究區(qū)植被的空間動態(tài)變化,差值為正時表示研究區(qū)植被增加,差值為負時表示研究區(qū)植被減少,差值為0時表示植被覆蓋狀態(tài)未發(fā)生改變。圖像差值法計算公式為

        ΔVFC=VFCyear2-VFCyear1。

        (4)

        式中:ΔVFC表示植被覆蓋的變化量;VFCyear1表示前一個時相的植被覆蓋數(shù)據(jù);VFCyear2表示后一個時相的植被覆蓋數(shù)據(jù)。

        圖像差值法能夠很好地反映研究區(qū)植被覆蓋的變化情況,根據(jù)圖像差值法將三屯河流域2000年與2020年的植被覆蓋度進行比較,分析三屯河流域植被覆蓋度絕對空間變化。

        根據(jù)傳統(tǒng)標準差值,將三屯河流域研究區(qū)植被覆蓋空間動態(tài)變化分為7個等級[35]:極度改善(0.67≤ΔVFC<1)、中度改善低(0.33≤ΔVFC<0.67)、輕微改善(0≤ΔVFC<0.33)、未改變(ΔVFC=0)、輕微衰退(-0.33≤ΔVFC<0)、中度衰退(-0.67≤ΔVFC<-0.33)、極度衰退(-1≤ΔVFC<-0.67)。

        2.3.4 植被覆蓋度穩(wěn)定性

        變異系數(shù)(Coefficient of Variation,CV)是一種能夠客觀反映數(shù)據(jù)評價結果的方法,計算結果均來自數(shù)據(jù)本身所包含的各種信息,不受主觀因素影響,與德爾菲法等相比,能更好、更客觀地反映數(shù)據(jù)的評價結果[36]。變異系數(shù)可以反映數(shù)據(jù)的波動程度,可以用來表達植被變化的穩(wěn)定性[37]。變異系數(shù)的計算公式為

        (5)

        (6)

        植被覆蓋的波動性與變異系數(shù)呈正比,變異系數(shù)值越大,植被覆蓋越敏感[38]。在ArcGIS軟件中選擇地圖代數(shù)模塊,通過柵格計算器計算2000—2020年植被覆蓋度的變異系數(shù),得到三屯河流域植被覆蓋的穩(wěn)定性狀況。

        3 結果與分析

        3.1 植被覆蓋度時空變化特征

        3.1.1 植被覆蓋時間分布特征

        三屯河流域20 a間植被覆蓋變化趨勢如圖1所示,2000—2020年,三屯河流域年均植被覆蓋度總體上呈現(xiàn)逐年增長的趨勢,在0.30~0.50之間變化,年均增長速度為4.93%/(10 a),植被覆蓋度最高的年份出現(xiàn)在2020年,植被覆蓋度為0.470 2;植被覆蓋度最低年份出現(xiàn)在2000年,植被覆蓋度為0.371 7。2012年的植被覆蓋度出現(xiàn)了輕微的下降趨勢,除2012年外,其余年份植被覆蓋較前一期相比均為上升趨勢,且2012—2015年植被覆蓋度的增長速度最快,年均增長速度為18.00%/(10 a)。

        圖1 植被覆蓋度年際變化Fig.1 Interannual change of vegetation coverage

        通過查閱相關資料發(fā)現(xiàn)2012年前后,三屯河流域內并未發(fā)生大規(guī)模自然災害,結合研究區(qū)氣溫降水資料,2012年氣溫為6.38 ℃,降水量為179.80 mm,2000年、2004年、2008年、2015年及2020年的氣溫分別為6.05、6.79、7.21、7.48、7.14 ℃,相應的降水量分別為198.30、259.20、174.50、218.00、77.70 mm。

        發(fā)現(xiàn)6個時相中2012年的氣溫與降水量與其他年份相比處于較為低下水平,這可能是造成2012年植被覆蓋度出現(xiàn)下降的原因。

        3.1.2 植被覆蓋空間動態(tài)變化特征

        通過對三屯河流域2000年、2004年、2008年、2012年、2015年及2020年6期植被覆蓋度分級圖(圖2)分析,發(fā)現(xiàn)研究區(qū)各年份植被覆蓋度分布表現(xiàn)出一定的規(guī)律性,即三屯河流域植被覆蓋總體表現(xiàn)為流域內兩端高、中間低的大面積片狀分布,上下游地區(qū)植被覆蓋度高,中游地區(qū)植被覆蓋度低,中、低植覆蓋區(qū)域被包裹著高植被覆蓋區(qū)域生長。以高植被覆蓋區(qū)域為中心向外發(fā)散,植被覆蓋等級逐步降低。

        圖2 三屯河流域6期植被覆蓋度分級Fig.2 Rating of vegetation coverage in Santun River Basin in six periods

        3.2 植被覆蓋度變化趨勢及穩(wěn)定性分析

        3.2.1 植被覆蓋度變化趨勢分析

        2000—2020年三屯河流域植被覆蓋度絕對空間變化分布如圖3所示,三屯河流域植被覆蓋情況整體呈現(xiàn)出改善的趨勢,植被覆蓋表現(xiàn)為改善的區(qū)域占比39.58%,其中極度改善的區(qū)域占比3.49%,主要以斑塊的形式分布于下游多處,中游少量區(qū)域及上游周邊區(qū)域;中度改善的區(qū)域占比7.35%,主要分布在極度改善周圍,上、中、下游均有分布;輕微改善的區(qū)域占比28.74%,遍布整個三屯河流域,其中以上游與中游地區(qū)分布最為廣泛。植被覆蓋表現(xiàn)為衰退的區(qū)域占比11.80%,其中極度衰退的區(qū)域占比1.16%,主要分布在下游區(qū)域,呈大面積斑塊狀;中度衰退的區(qū)域占比2.52%,仍主要分布的下游地區(qū),上游地區(qū)也存在零星分布;輕微衰退的區(qū)域占比8.12%,分布在上游流域邊界周圍,呈半包圍狀包圍上游,于中下游交接處大面積片狀分布;植被覆蓋未改變區(qū)域站總區(qū)域的48.62%,可見改善大于衰退,三屯河流域植被覆蓋呈改善變化趨勢。

        圖3 三屯河流域植被覆蓋度絕對空間變化分布Fig.3 Distribution of absolute spatial change of vegetation coverage in Santun River Basin

        3.2.2 植被覆蓋度穩(wěn)定性分析

        為分析研究區(qū)2000—2020年植被覆蓋的的穩(wěn)定性,繪制三屯河流域植被覆蓋穩(wěn)定性圖(圖4),經統(tǒng)計分析發(fā)現(xiàn)在三屯河流域內大部分地區(qū)的植被覆蓋度處于弱變異與中等變異,只有零星少部分地區(qū)處于強變異。弱變異地區(qū)主要分布在上游地區(qū),中等變異近乎涵蓋整個流域,強變異主要散亂分布在上游高海拔流域邊界處與下游居民居住地。由此可見,三屯河流域植被覆蓋整體相對穩(wěn)定,植被覆蓋表現(xiàn)為強變異,除了受到自然因素的影響,還受到人為因素的影響。

        圖4 三屯河流域植被覆蓋穩(wěn)定性分析Fig.4 Stability of vegetation coverage in Santun River Basin

        4 植被覆蓋變化影響因素的討論

        三屯河流域是典型的半干旱地區(qū),其植被覆蓋變化受到自然因素與人類活動因素的共同的影響。

        4.1 自然因素影響

        三屯河流域屬中溫帶大陸性半荒漠干旱氣候,是典型的內陸半干旱地區(qū)。如圖5所示,2000—2020年20 a間降水與氣溫均呈波動變化,降水量整體趨勢表現(xiàn)為輕微下降,氣溫整體趨勢呈上升狀態(tài);20 a間三屯河流域植被覆蓋度呈上升趨勢。據(jù)前人研究結果可知,氣溫的增高和降水的增加對新疆干旱區(qū)植被覆蓋具有正向的促進作用[39-40],三屯河流域的主要補給方式是大氣降水和冰川融水,氣溫升高導致的冰川消融對三屯河流域水資源的補給貢獻了一定力量,同時,2011年中央一號文件提出的水資源政策有效控制了地下水的過度開采,也為區(qū)域植被的增加提供了水資源保障。

        圖5 三屯河流域氣溫、降水及植被覆蓋變化趨勢Fig.5 Changes in temperature,precipitation and vegetation coverage in the Santun River Basin

        4.2 人類活動因素影響

        人類因素對植被覆蓋的影響主要表現(xiàn)在土地利用類型的改變上。為分析研究區(qū)20 a間土地利用的變化情況,對2000年與2020年土地利用數(shù)據(jù)進行處理,得到2000—2020年土地利用類型面積轉移矩陣(表1),用以分析三屯河流域2000—2020年20 a間不同地類之間的面積轉換情況。如表1所示,2020年較2000年相比,經過轉入轉出換算后可知,草地面積增加141.88 km2,城鄉(xiāng)工礦面積增加34.97 km2,耕地面積增加36.98 km2,林地面積減少166.87 km2,水域面積減少118.90 km2,未利用土地面積增加71.93 km2。根據(jù)轉移矩陣的計算結果可知,轉入為草地的各類土地中,林地的貢獻率最大,高達59.95%,通過查閱資料可知20 a間三屯河流域未發(fā)生大面積森林火災,故林地林地向草地的轉入可認定為人為因素對植被覆蓋的逆向影響。隨著城市化的發(fā)展,由草地和耕地的轉入導致城鄉(xiāng)工礦面積的增加也可視為人為因素對植被覆蓋的逆向影響。20 a間林地轉入面積全部來自草地轉入,轉入為耕地的各類型土地中,草地的貢獻率高達64.06%,而研究區(qū)主要的種植作物包括玉米、籽用葫蘆及葡萄,單位面積植被覆蓋度明顯高于草地,這表明在人為因素對植被覆蓋情況具有正向影響。未利用土地向耕地的轉移可以視為人為因素對植被覆蓋的正向影響??傮w上看,各類型的土地類型的轉移是交替進行的,但這種轉移對植被覆蓋的正向影響大于逆向影響,故研究區(qū)植被覆蓋呈增長趨勢。因此除了自然因素影響下植被短期的變化外,人類活動也是影響植被覆蓋長期演化趨勢的一個重要因素。

        表1 2000—2020年土地利用類型面積轉移矩陣Table 1 Transfer matrix of land use type area in 2000-2020

        5 結 論

        本文定量分析了三屯河流域2000—2020年20 a間植被覆蓋的時空變化特征,并對其變化趨勢與穩(wěn)定性進行了分析,得到如下結論:

        (1)時間上,三屯河流域植被覆蓋度呈現(xiàn)出逐年增長趨勢,特別是2012—2015年的年均增長速度為4.93%/(10 a),除2012年外,其余年份植被覆蓋較前一期相比均為上升趨勢,2012年植被覆蓋出現(xiàn)下降趨勢可能是受到降水與氣溫的影響;空間上,三屯河流域各年份植被覆蓋度分布表現(xiàn)出一定的規(guī)律性,流域內植被覆蓋兩端高、中間低,呈大面積片狀分布,低植被覆蓋包裹著高植被覆蓋生長。

        (2)三屯河流域植被覆蓋整體的變化趨勢為改善的趨勢,改善區(qū)域占比39.58%,衰退區(qū)域占比11.80%,未改變區(qū)域占比48.62%;三屯河流域內植被覆蓋整體相對穩(wěn)定,弱變異與中等變異分布于流域內大部分地區(qū),強變異零星分布在上游邊界處與下游居住地。研究區(qū)植被覆蓋表現(xiàn)為衰退的區(qū)域面積雖然較小,但應引起研究者及相關部門的關注。

        (3)氣溫和降水對三屯河流域植被覆蓋變化在短期內有一定的影響,但人類活動因素對研究區(qū)植被覆蓋長期演化及趨勢的影響更為直接,主要體現(xiàn)在土地利用類型的改變上。城市化的擴大及經濟的發(fā)展造成部分林地植被退化是人類活動對三屯河流域植被覆蓋的逆向影響;耕地的開墾擴張、林地種植是人類活動對三屯河流域植被覆蓋的正向影響??傮w上看,人類活動對三屯河流域植被覆蓋的正向影響大于逆向影響,研究區(qū)植被覆蓋呈增長趨勢。研究區(qū)植被覆蓋的增加有一部分貢獻于由開墾草地導致的耕地的增加,但開墾會降低區(qū)域生態(tài)環(huán)境平衡,長期來說不利于區(qū)域植被穩(wěn)定性及生態(tài)環(huán)境的長期可持續(xù)發(fā)展。

        猜你喜歡
        屯河覆蓋度植被
        呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
        基于植被復綠技術的孔植試驗及應用
        河北地質(2022年2期)2022-08-22 06:24:04
        基于NDVI的晉州市植被覆蓋信息提取
        低覆蓋度CO分子在Ni(110)面的吸附研究
        綠色植被在溯溪旅游中的應用
        新疆頭屯河洪水特征分析
        頭屯河水庫溢洪道除險加固初探
        新疆頭屯河防洪工程水土流失預測
        新疆頭屯河流域水資源配置
        基于原生植被的長山群島植被退化分析
        久久久国产精品免费a片3d| 精品人妻一区二区三区狼人| 国内自拍愉拍免费观看| 精品免费久久久久久久| japanese无码中文字幕| 国产成人AV乱码免费观看| 伊人久久亚洲精品中文字幕| 亚洲日韩中文字幕在线播放| 免费观看又污又黄的网站| 欧美综合图区亚洲综合图区| 久久久黄色大片免费看| 在线播放五十路熟妇| 国产精品卡一卡二卡三| 欧洲亚洲色一区二区色99| 国产一级黄色片在线播放| 成人精品天堂一区二区三区 | 精品少妇人妻久久免费| 精品人妻中文av一区二区三区| 久久久受www免费人成| 日中文字幕在线| 久久蜜桃一区二区三区| 国产精品无码制服丝袜| 男人扒开女人下面狂躁小视频| 国产成人精品三上悠亚久久| 国产亚洲3p一区二区| 久久午夜羞羞影院免费观看| 成人无码视频| 91国产自拍视频在线| 日本高清乱码中文字幕| 国产精品爽黄69天堂a| 久久成人永久免费播放| 东京热加勒比国产精品| 精品久久久久久久久午夜福利| 亚洲国产中文在线二区三区免| 国产激情视频免费观看| 国产av一区二区三区性入口 | 中文字幕亚洲乱码熟女在线| 久久国产精品亚洲婷婷片| 午夜亚洲www湿好爽| 国产极品视觉盛宴在线观看| 成熟妇女毛茸茸性视频|