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

        ?

        基于遙感生態(tài)指數(shù)的宜昌市區(qū)生態(tài)環(huán)境變化研究

        2022-06-06 12:57:42羅瑪詩藝張城芳陳興芳
        地理空間信息 2022年5期
        關(guān)鍵詞:綠度干度宜昌市

        羅瑪詩藝,陶 榮*,張城芳,陳興芳

        (1. 湖北省空間規(guī)劃研究院,湖北 武漢 430071;2. 武漢華夏理工學(xué)院 土木建筑工程學(xué)院,湖北 武漢 430223;3. 武漢理工大學(xué)資源與環(huán)境工程學(xué)院,湖北 武漢 430070;4. 華中師范大學(xué) 城市與環(huán)境科學(xué)學(xué)院,湖北 武漢 430079)

        自然生態(tài)系統(tǒng)是人類賴以生存和發(fā)展的基礎(chǔ)。近年來,隨著大量以城鎮(zhèn)化為主的人類活動的增加,土地利用情況發(fā)生了很大的變化,對各種尺度上的生態(tài)系統(tǒng)均造成了嚴(yán)重擾動。因此,有效評價某地區(qū)的生態(tài)環(huán)境質(zhì)量和動態(tài)監(jiān)測生態(tài)環(huán)境狀況的時空變化,一方面符合生態(tài)可持續(xù)發(fā)展的實際需求,另一方面將為區(qū)域規(guī)劃人員進(jìn)行國土空間規(guī)劃提供參考。

        目前,遙感技術(shù)以其方便、快捷和客觀等優(yōu)勢,成為生態(tài)環(huán)境研究中不可或缺的一部分。遙感技術(shù)運(yùn)用于生態(tài)環(huán)境質(zhì)量評價方面的研究可分為3 類:①基于遙感獲取的單個指標(biāo)生態(tài)環(huán)境評價方法,如根據(jù)遙感解譯獲取的土地利用/覆被的分布特征和變化趨勢分析生態(tài)環(huán)境變化[1],土地利用/覆被變化雖可作為生態(tài)環(huán)境變化的一個方面,但生態(tài)環(huán)境是一個復(fù)雜的整體,單一的生態(tài)環(huán)境指標(biāo)不能反映生態(tài)環(huán)境全面且真實的變化情況;②運(yùn)用不同指數(shù)所構(gòu)建的指標(biāo)體系綜合反映區(qū)域生態(tài)環(huán)境的變化情況,如國家頒布的《生態(tài)環(huán)境狀況評價技術(shù)規(guī)范》則利用遙感反演獲取的生物豐度指標(biāo)、植被覆蓋指標(biāo)、水網(wǎng)密度指標(biāo)、土地脅迫指標(biāo)和污染負(fù)荷指標(biāo)等來綜合反映區(qū)域生態(tài)環(huán)境狀況[2],但各指標(biāo)在體系中的權(quán)重是人為賦予的,因此其結(jié)果存在不合理性;③基于遙感信息的空間化,從遙感影像中直接獲取與生態(tài)環(huán)境息息相關(guān)的指數(shù),從而減少對指標(biāo)權(quán)重的人為干預(yù),實現(xiàn)生態(tài)環(huán)境變化的可視化,如對于長株潭地區(qū)生態(tài)環(huán)境的研究,肖瑛[3]選取了植被覆蓋度、土壤含水量、地形3 個生態(tài)環(huán)境要素作為評價指標(biāo),徐涵秋[4]采用RSEI評價了一個地區(qū)的生態(tài)環(huán)境質(zhì)量,該體系中主要包括綠度指標(biāo)、濕度指標(biāo)、干度指標(biāo)和熱度指標(biāo)。

        宜昌市作為長江流域的重要節(jié)點(diǎn)城市,其生態(tài)環(huán)境關(guān)乎長江整體流域的發(fā)展。在全球氣候變暖的背景下,當(dāng)?shù)剡^度挖掘資源和毀林開荒等活動使得當(dāng)?shù)氐乃亮魇Ш退|(zhì)污染等環(huán)境問題越來越嚴(yán)重。本文采用RSEI 對宜昌市西陵區(qū)、點(diǎn)軍區(qū)、猇亭區(qū)、夷陵區(qū)和伍家崗區(qū)(以下簡稱宜昌市區(qū))的生態(tài)環(huán)境進(jìn)行了評價,并分析了該區(qū)域RSEI 的時間變化情況和空間分布情況。

        1 研究區(qū)概況與數(shù)據(jù)處理

        1.1 研究區(qū)概況

        本文選取的研究對象是宜昌市區(qū),位于湖北省西南部。該區(qū)域?qū)賮啛釒Ъ撅L(fēng)氣候,年均氣溫為16.9℃,年均降水量為992.1~1 404.1 mm,植被類型以常綠闊葉落葉林為主;地形西北高、東南低,西北部為高山植被覆蓋區(qū),東南部丘陵平原相間分布。研究區(qū)全域林業(yè)和耕地面積占比較大,中南部城市建設(shè)用地占比較高,建成區(qū)已達(dá)170 km2。社會經(jīng)濟(jì)方面以電力、化工、食品醫(yī)藥、裝配制造為四大支柱產(chǎn)業(yè),2018年5個市轄區(qū)所創(chuàng)造的生產(chǎn)總值占宜昌全市生產(chǎn)總值的36.8%[5],經(jīng)濟(jì)實力較強(qiáng)。

        1.2 數(shù)據(jù)來源與預(yù)處理

        本文以Landsat8 遙感影像為主要數(shù)據(jù)源,為保證影像能最大限度地觀測到綠度信息,影像獲取時間盡量處于該區(qū)域主要植被的生長期;同時為避免時間/季節(jié)變化效應(yīng)引起不同光譜的響應(yīng)問題,不同年份的多期影像盡量保證在一年中相近的時間段內(nèi),使植被具有相近的物候和生長狀態(tài)。宜昌市夏季多云多雨,因此大量5—10月的Landsat8影像數(shù)據(jù)云量較大,不能滿足應(yīng)用條件,最終只有2013年8月、2015年8月和2019年10月3期影像能滿足云量較少的條件,影像詳細(xì)信息如表1 所示。Landsat8 影像數(shù)據(jù)均從美國地質(zhì)勘探局(USGS)官網(wǎng)(https://earthexplorer.usgs.gov)下載。

        表1 影像數(shù)據(jù)詳細(xì)參數(shù)

        基于ENVI 軟件,結(jié)合LaSRC 2.0,本文對Landsat OLI數(shù)據(jù)進(jìn)行輻射定標(biāo)和大氣校正,采用云量最小的策略,無縫鑲嵌遙感影像;并通過影像裁剪生產(chǎn)特定區(qū)域的可用于地表參數(shù)反演的影像數(shù)據(jù),詳細(xì)如圖1所示。

        圖1 數(shù)據(jù)預(yù)處理路線圖

        2 研究方法

        2.1 指標(biāo)選取與反演

        本文借助徐涵秋等學(xué)者提出的RSEI 去評價宜昌市區(qū)的生態(tài)環(huán)境。RSEI由綠度、濕度、干度和熱度4個與人類生產(chǎn)生活息息相關(guān)的指標(biāo)耦合而成。4 個指標(biāo)可直接從遙感影像中獲取,可綜合反映某一地區(qū)的生態(tài)環(huán)境狀況[6]。

        綠度指標(biāo)是表示研究區(qū)內(nèi)植被狀況的一個綜合感知指標(biāo)。本文選取歸一化植被指數(shù)(NDVI)來表示綠度。

        式中 ,ρNIR和 ρRed為Landsat8 OLI 影像近紅外波段和紅光波段經(jīng)大氣校正后的反射率數(shù)據(jù)。

        濕度指標(biāo)反映的是水體、土壤和植被的含水量情況[7]。本文采用Landsat8 OLI傳感器所獲影像的纓帽變換的濕度分量——地表濕度(LSM)表示濕度指標(biāo)。詳細(xì)推導(dǎo)過程見參考文獻(xiàn)[8]。

        式中, ρBlue、 ρGreen、 ρRed、 ρNIR、 ρSWIR1、 ρSWIR2分別為Landsat8 OLI 影像藍(lán)、綠、紅、近紅外、短波紅外1 和短波紅外2 六個波段經(jīng)大氣校正后的反射率數(shù)據(jù)。

        熱度指標(biāo)反映地表熱量的分布情況,可采用從熱紅外影像中獲取的地表溫度(LST)來表示。本文利用Jiménez-Munoz J C[9]等提出的單通道算法來反演LST,即

        式中, Lsen為衛(wèi)星傳感器測得的輻射強(qiáng)度,單位為W/m2/sr/μm;Tsen為亮溫;γ、δ 為基于普朗克函數(shù)的兩個參數(shù),其中對于TIRS 10波段而言bγ為1 324;ε 為地表比輻射率,具體數(shù)字見參考文獻(xiàn)[9];ψ1、ψ1和ψ1是3個大氣水含量的函數(shù),具體函數(shù)表達(dá)和推導(dǎo)過程見參考文獻(xiàn)[10]。

        Lsen可通過輻射校正將影像的亮度值(DN)轉(zhuǎn)換為大氣頂部的光譜輻射值,即

        式中, ML為波段 λ 的絕對輻射定標(biāo)系數(shù);AL為波段λ的調(diào)整偏移量,這兩個值均能從影像的元數(shù)據(jù)文件(MTL文件)中獲取。

        Tsen可根據(jù)黑體輻射能量的分布定律將輻射能量轉(zhuǎn)換為亮溫,即

        式中,對于TIRS 10波段而言,K1和K2是熱紅外波段的定標(biāo)常數(shù),分別為774.89 W/m2/sr/μm和1 321.08 K。

        干度指標(biāo)反映的是裸土、人工建設(shè)的道路、建筑等不透水面造成的地表“干化”情況,能在一定程度上反映人類活動對生態(tài)環(huán)境狀況的影響,一般由反映裸土(SI)和不透水面(IBI)[11]的指數(shù)綜合而成。本文參照參考文獻(xiàn)[12],采用一個綜合植被指數(shù)NDBSI來估算。針對Landsat8影像,其計算方法為:

        式中 ,ρGreen、 ρRed、 ρNIR、 ρSWIR1、 ρSWIR2分別為Landsat8 OLI影像綠、紅、近紅外、短波紅外1和短波紅外2五個波段經(jīng)大氣校正后的反射率數(shù)據(jù)。

        2.2 RSEI的計算

        由于宜昌市區(qū)內(nèi)存在長江、西北口水庫等大面積水體,為避免大片水域影響主成分變換的荷載分布,使?jié)穸戎笜?biāo)能更準(zhǔn)確地反映土壤和植被的含水量信息,在計算RSEI 前,本文采用MNDWI 指數(shù)[13]來提取水體面積,用于掩膜水體信息。其計算方式為:

        式中, ρGreen、 ρMIR分別為影像數(shù)據(jù)中綠、中紅外波段經(jīng)大氣校正后的反射率數(shù)據(jù)。

        RSEI 由綠度、濕度、干度、熱度4 個指標(biāo)構(gòu)成。為了消除各分量指標(biāo)量綱和數(shù)值范圍的不同,在執(zhí)行主成分分析前,本文對4個指標(biāo)進(jìn)行歸一化處理,即

        式中,NI 為歸一化后的指標(biāo)值; I 為指標(biāo)的數(shù)值;Imax和Imin分別為指標(biāo)的最大值和最小值。

        將4個歸一化后的指標(biāo)借助ENVI軟件融合為一個圖像,并進(jìn)行主成分分析;再根據(jù)各主成分對其不同的貢獻(xiàn)率,保留貢獻(xiàn)率最大的主成分。由于第一個成分(PC1)可解釋超過80%的RSEI,集中了4個指標(biāo)的大部分特征,可選擇PCA的PC1來表示RSEI,即

        其中用1 減去PC1是為了使大的數(shù)值可以代表好的生態(tài)狀況。

        為了比較不同時期生態(tài)環(huán)境質(zhì)量的變化情況,將實際的RSEI值進(jìn)一步標(biāo)準(zhǔn)化為0~1的范圍,即

        RSEI值越高代表生態(tài)環(huán)境越好。

        3 研究結(jié)果與分析

        3.1 RSEI評價模型檢驗

        本文通過4個分量指標(biāo)和RSEI指數(shù)的相關(guān)性以及分量指標(biāo)相互之間的相關(guān)性來檢驗RSEI 的可行性與綜合代表程度。由表2 可知,在4 個分量指標(biāo)中,干度指標(biāo)的平均相關(guān)度最高,其均值為0.775 5,綠度、熱度和濕度指標(biāo)的3期平均相關(guān)度均值分別為0.679 5、0.575 0和0.650 1;3期RSEI對各指標(biāo)的平均相關(guān)度均值為0.731,而4 個分量指標(biāo)的平均相關(guān)度僅為0.67,表明RSEI 指數(shù)相對于單一指標(biāo)來說更具代表性,能從整體上評價區(qū)域生態(tài)質(zhì)量;從各指標(biāo)分量與RSEI的相關(guān)性來看,干度與熱度指標(biāo)對RSEI 的貢獻(xiàn)為負(fù)值,綠度與濕度指標(biāo)對RSEI 的貢獻(xiàn)為正值,表明干度與熱度指標(biāo)對城市生態(tài)環(huán)境起消極作用,而綠度與濕度指標(biāo)的提高可以更好地促進(jìn)城市生態(tài)環(huán)境的優(yōu)化。

        表2 RSEI與各指標(biāo)的相關(guān)性分析

        3.2 生態(tài)環(huán)境質(zhì)量狀況評價

        宜昌市區(qū)2013—2019 年RSEI 和4 個分量指標(biāo)統(tǒng)計如表3所示,可以發(fā)現(xiàn),2013年、2015年和2019年的RSEI 均值分別為0.633、0.620 和0.529,生態(tài)環(huán)境質(zhì)量呈逐漸下降趨勢,整體降幅達(dá)到16.4%。從空間分布特征來看(圖2),宜昌市區(qū)西部山地地區(qū)的生態(tài)質(zhì)量明顯高于東南部的平原地區(qū),西部的夷陵區(qū)與點(diǎn)軍區(qū)由于植被覆蓋度較好,人口壓力和經(jīng)濟(jì)開發(fā)力度較小,大部分地區(qū)的生態(tài)環(huán)境質(zhì)量略高于其他三區(qū)。生態(tài)質(zhì)量較差的地區(qū)主要集中在東南沿江地區(qū)以及黃泊河流域附近。結(jié)合2011—2030年宜昌市中心城區(qū)空間結(jié)構(gòu)規(guī)劃圖(圖3)發(fā)現(xiàn),宜昌市規(guī)劃發(fā)展呈現(xiàn)沿江沿河向南部發(fā)展的特點(diǎn),規(guī)劃了伍家崗組團(tuán)和白洋組團(tuán)等多個建設(shè)集聚區(qū)。隨著多個工業(yè)加工園和產(chǎn)業(yè)園相繼建成投入使用,一方面對耕地和林地等生態(tài)用地有較大地破壞,另一方面工業(yè)園與產(chǎn)業(yè)園所產(chǎn)生的工業(yè)廢物難以進(jìn)入生態(tài)循環(huán),對當(dāng)?shù)氐纳鷳B(tài)環(huán)境造成了較大影響。

        圖3 2011—2030年宜昌市中心城區(qū)空間結(jié)構(gòu)圖注:來源于中國城市規(guī)劃設(shè)計研究院。

        表3 各指標(biāo)均值分析

        圖2 宜昌市區(qū)2013—2019年RSEI區(qū)域分布圖

        從4個分量指標(biāo)來看,與RSEI呈正相關(guān)的濕度與綠度指標(biāo)呈先升后降的趨勢,整體降幅分別為2.8%和1.5%;而對RSEI 的貢獻(xiàn)為負(fù)的干度與熱度指標(biāo)均呈上升趨勢,增幅分別為12.5%和2.7%,表明宜昌市區(qū)生態(tài)環(huán)境的下降與植被、土壤和植被含水量的減少、地表熱量和不透水面的增加有關(guān)系。植被的減少和城市不透水面的增加在一定程度上導(dǎo)致了宜昌市區(qū)生態(tài)環(huán)境質(zhì)量的下降。

        3.3 生態(tài)環(huán)境動態(tài)變化分析

        根據(jù)徐涵秋[6]提出的分級研究方法,本文以0.2為分級間隔,將RSEI分為[0 0.2]、[0.2 0.4]、[0.4 0.6]、[0.6 0.8]、[0.8 1]五個等級,分別借鑒等級差值運(yùn)算和土地利用變化研究的轉(zhuǎn)移矩陣法,探究了生態(tài)環(huán)境質(zhì)量的詳細(xì)變化情況。2013—2015 年RSEI 等級轉(zhuǎn)移矩陣如表4所示,可以看出,第四等級和第三等級存在明顯的轉(zhuǎn)出,第一等級和第五等級轉(zhuǎn)入與轉(zhuǎn)出不明顯;經(jīng)統(tǒng)計分析可知,2013—2015年生態(tài)環(huán)境等級為第三等級地區(qū)的面積通過轉(zhuǎn)化增加了120.212 km2,其原因是部分第四等級地區(qū)轉(zhuǎn)化為第三等級,第四等級地區(qū)的面積減少了190.282 km2;對角線上的值表示生態(tài)環(huán)境質(zhì)量不變,對角線上方區(qū)域為高等級指標(biāo)轉(zhuǎn)化為低等級指標(biāo),表示生態(tài)環(huán)境質(zhì)量下降,而下方區(qū)域為低等級指標(biāo)轉(zhuǎn)化為高等級指標(biāo),表示生態(tài)環(huán)境質(zhì)量有所提高,對角線上方的面積為377.583 km2,指標(biāo)下降面積占變化值的71.3%,因此2013—2015年宜昌市生態(tài)質(zhì)量變差區(qū)域大于變好區(qū)域。

        表4 2013—2015年RSEI等級轉(zhuǎn)移矩陣/km2

        本文將 2015 年的 RSEI 值減去 2013 年的 RSEI 值,得差值為正、零、負(fù)分別代表生態(tài)環(huán)境變好、不變和變差。由圖4 可知,2013—2015 年大部分地區(qū)生態(tài)環(huán)境質(zhì)量沒有大的變化;生態(tài)環(huán)境變好的區(qū)域主要在西北口水庫地區(qū),查閱2015 年宜昌市的氣象資料可知,宜昌市區(qū)經(jīng)歷了大面積降雨,西北口水庫庫容增加,為水庫周圍林地和草地的恢復(fù)提供了充足的水源,因此該區(qū)生態(tài)環(huán)境質(zhì)量有所轉(zhuǎn)變;生態(tài)環(huán)境變差的區(qū)域主要集中在點(diǎn)軍區(qū)的中部、夷陵區(qū)的西南部以及黃泊河流域附近,且呈點(diǎn)狀分布,這是由于宜昌市區(qū)在進(jìn)行經(jīng)濟(jì)建設(shè)時,對局部地區(qū)過度開發(fā),使得當(dāng)?shù)氐脑忌鷳B(tài)環(huán)境發(fā)生破壞所致。

        圖4 2013—2015年等級變化圖

        2015—2019年RSEI等級轉(zhuǎn)移矩陣如表5所示,可以看出,第二等級與第三等級地區(qū)的面積有所增加,其增加面積分別為387.175 km2和1 392.179 km2,第二等級地區(qū)面積增加的原因主要是第三等級和第四等級不斷轉(zhuǎn)化為第二等級;而第四等級地區(qū)的面積減少了1 758.7 km2,其原因是大部分第四等級轉(zhuǎn)化為第三等級和第二等級。

        表5 2015—2019年RSEI等級轉(zhuǎn)移矩陣/km2

        由2015—2019 年的等級變化圖(圖5)來看,區(qū)域生態(tài)環(huán)境變差的區(qū)域大面積增長,變差區(qū)域主要集中在北部與中部地區(qū),且呈局部地區(qū)集中分布的特點(diǎn);北部地區(qū)土地覆蓋原為森林草地,但由于人們的大肆開發(fā),土地覆蓋發(fā)生明顯變化,生態(tài)環(huán)境變化較為明顯,同時由于2019 年的影像獲取時期為10 月,一方面該地區(qū)的植被處于凋落期,相對的綠度指標(biāo)有所降低,另一方面該時期亞熱帶地區(qū)降水相對減少,使得干度指標(biāo)有所增加,使得RSEI值有所下降。

        圖5 2015—2019年等級變化圖

        4 結(jié) 語

        本文基于 2013 年、2015 年和 2019 年 3 期遙感數(shù)據(jù),反演獲取了干度、熱度、濕度和綠度4 個指標(biāo),通過主成分分析法構(gòu)建了RSEI;并結(jié)合差值運(yùn)算和等級轉(zhuǎn)移矩陣研究方法,全面分析了宜昌市區(qū)生態(tài)環(huán)境的時空變化情況,為宜昌市城市發(fā)展與規(guī)劃提供了有效的依據(jù)。

        1)濕度、綠度、干度、熱度4個指標(biāo)所構(gòu)成的生態(tài)環(huán)境指標(biāo)體系比單一指標(biāo)更能反映生態(tài)環(huán)境的變化。從單一指標(biāo)與RSEI 的相關(guān)性來看,綠度與濕度指標(biāo)對生態(tài)環(huán)境的貢獻(xiàn)為正,而干度與熱度指標(biāo)對生態(tài)環(huán)境的貢獻(xiàn)為負(fù)。

        2)由3期遙感指數(shù)的均值來看,宜昌市區(qū)生態(tài)環(huán)境質(zhì)量呈下降態(tài)勢。從時間跨度來看,2013年、2015年和2019 年宜昌市區(qū)RSEI 均值分別為0.633、0.620 和0.529,表明2013—2019 年宜昌市區(qū)生態(tài)環(huán)境質(zhì)量在逐漸變差。從空間跨度來看,東南方向整體的生態(tài)環(huán)境質(zhì)量較差,而西北方向的RSEI 指數(shù)相對其他地區(qū)較高。

        3)在RSEI 等級變化研究中,以中等等級和較好等級變化最為明顯,且在轉(zhuǎn)化中,等級下降比例明顯高于等級上升比例,這也表明宜昌市區(qū)生態(tài)環(huán)境質(zhì)量呈下降態(tài)勢。

        猜你喜歡
        綠度干度宜昌市
        寧夏草地綠度時空變化圖譜及驅(qū)動因素分析
        全球森林綠度異常對干旱事件的響應(yīng)
        2003~2018年干旱對中國森林綠度異常的影響
        宜昌市舉辦中小學(xué)實驗室安全管理培訓(xùn)班
        淺談唐卡中綠度母的造型特點(diǎn)
        湖北省宜昌市第一中學(xué)2014—2015學(xué)年第一學(xué)期高二年級期末考試第21題
        注汽鍋爐應(yīng)用干度自控系統(tǒng)提高稠油采收率
        務(wù)實創(chuàng)新 開拓奉獻(xiàn)——湖北宜昌市第三批援藏干部援助西藏加查兩周年掃描
        民族大家庭(2015年5期)2015-08-21 08:01:36
        我國青少年五人制足球推廣研究——以宜昌市為例
        稠油注汽鍋爐蒸汽干度測控裝置的應(yīng)用
        国产精品久久久免费精品| 久久久久成人精品免费播放| 午夜国产小视频在线观看黄| 亚洲国产综合在线亚洲区亚洲av | 国产极品美女高潮无套在线观看| 亚洲无码中文字幕日韩无码| 成人国产高清av一区二区三区| 国产精品久免费的黄网站 | 久久精品中文闷骚内射| 小12箩利洗澡无码视频网站| 国产精彩刺激对白视频| 亚洲久无码中文字幕热| 在线观看黄片在线播放视频| 人妻少妇69久久中文字幕| 国产色在线 | 亚洲| 精品人妻无码视频中文字幕一区二区三区 | 99久久综合九九亚洲| 亚洲一区二区三区最新视频 | 亚洲精品中文字幕一二| 精品乱人伦一区二区三区| 99精品国产兔费观看久久99| www.av在线.com| 最新国产激情视频在线观看| 国产午夜福利精品一区二区三区 | 香蕉人妻av久久久久天天| 亚洲AV无码日韩一区二区乱| 男男亚洲av无一区二区三区久久| 少妇被粗大的猛烈进出69影院一| 在线一区不卡网址观看| 亚洲视频精品一区二区三区| 精品人妻一区二区三区浪人在线| 丰满岳妇乱一区二区三区| 无码人妻精品一区二区三区下载| 99视频一区二区日本| 国产精品爽爽久久久久久竹菊| 精品少妇一区二区三区视频| 扒下语文老师的丝袜美腿| 久久精品色福利熟妇丰满人妻91| 99精品欧美一区二区三区| 免费黄色福利| 亚洲精品中字在线观看|