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

    基于遙感生態(tài)指數(shù)(RSEI)的水電開發(fā)區(qū)生態(tài)環(huán)境變化分析:以清江中下游地區(qū)為例

    2019-11-25 05:24:12劉索玄袁艷斌
    關(guān)鍵詞:綠度清江梯級

    劉索玄,袁艷斌,趙 皞,李 倩

    (武漢理工大學(xué)資源與環(huán)境工程學(xué)院,湖北 武漢 430070)

    隨著對能源和水資源的需求增加,水力發(fā)電作為一種可再生、無污染的能源被廣泛開發(fā)利用[1]。大型水利水電工程往往建于生態(tài)脆弱區(qū),工程建設(shè)在創(chuàng)造巨大經(jīng)濟(jì)效益的同時,也會對生態(tài)環(huán)境系統(tǒng)產(chǎn)生脅迫[2]。如何準(zhǔn)確分析評價(jià)水利工程對生態(tài)環(huán)境的影響,將對區(qū)域生態(tài)環(huán)境保護(hù)和可持續(xù)發(fā)展產(chǎn)生重要意義。

    區(qū)域生態(tài)環(huán)境評價(jià)方法主要有模糊評判法、層次分析法、網(wǎng)絡(luò)分析法和綜合指數(shù)法[3]。以往生態(tài)評價(jià)方法所選指標(biāo)獲取難度大,指標(biāo)權(quán)重設(shè)定具有主觀性,不能準(zhǔn)確反映整個研究區(qū)的評價(jià)結(jié)果。遙感技術(shù)的發(fā)展推動了區(qū)域生態(tài)環(huán)境評價(jià)研究的發(fā)展:付剛等[4]利用遙感數(shù)據(jù)和GIS對北京市生態(tài)脆弱性進(jìn)行評價(jià);周雪欣等[5]選取植被覆蓋度、土壤侵蝕程度和石漠化敏感性程度指標(biāo),采用綜合指數(shù)法定量評價(jià)北盤江流域生態(tài)環(huán)境質(zhì)量。利用遙感技術(shù)獲取評價(jià)指標(biāo),便于快速、全面評價(jià)生態(tài)環(huán)境時空變化,已成為水電開發(fā)區(qū)生態(tài)評價(jià)的創(chuàng)新點(diǎn)。

    作為全國水電梯級開發(fā)的首個試點(diǎn),清江中下游成為我國生態(tài)環(huán)境建設(shè)的重點(diǎn)區(qū)域。關(guān)于清江水電開發(fā)的研究,主要針對水庫汛限水位[6]、洪水概率[7]或水電梯級聯(lián)合調(diào)度[8]以及運(yùn)營模式[9]等方面,而面向整個水電開發(fā)區(qū)長期的綜合性生態(tài)環(huán)境評價(jià)研究還較少。因此,根據(jù)水電梯級開發(fā)的時間順序,選擇清江中下游地區(qū)1987、2004和2015年3期遙感數(shù)據(jù)反演各期生態(tài)環(huán)境指標(biāo),運(yùn)用主成分分析法構(gòu)建新型遙感生態(tài)指數(shù)(remote sensing ecology index,RSEI,IRSE),對該地區(qū)1987—2015年生態(tài)環(huán)境時空變化進(jìn)行分析,以期為區(qū)域生態(tài)環(huán)境保護(hù)和綜合治理提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    清江是湖北省境內(nèi)第2大長江支流,發(fā)源于利川市齊岳山龍洞溝。地勢西高東低,流經(jīng)區(qū)域以山地為主(山地面積占比>80%),總海拔高差為1 330 m。流域內(nèi)年平均無霜期為285.4 d,年均降水量為1 335.5 mm,水能資源較為豐富,其中,水資源主要集中在恩施以下的干流河段(占比>80%)。清江梯級開發(fā)工程始于1987年,規(guī)劃在清江中下游地區(qū)建設(shè)“水布埡—隔河巖—高壩洲”3個梯級工程。截至2014年,清江13座水電站全面竣工。其中,下游的隔河巖和高壩洲水電站分別建成于1994和2003年;位于恩施州巴東縣境內(nèi)的水布埡水電站于2009年6月底建成。此次研究區(qū)域?yàn)榍褰邢掠渭捌渲苓叺貐^(qū),包括宜都市、五峰縣、長陽縣和巴東縣部分地區(qū),面積共計(jì)8 783.12 km2(圖1)。

    圖1 研究區(qū)概覽Fig.1 Overview of the study area

    2 數(shù)據(jù)來源與研究方法

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

    以清江中下游為界,以水電站梯級建設(shè)起始時間、下游梯級建成時間、水電梯級全面建成時間為依據(jù),統(tǒng)一選用美國地質(zhì)調(diào)查局(USGS)網(wǎng)站(http:∥glovis.usgs.gov/)提供的1987年4月17日、2004年5月17日的Landsat 5 TM影像及2015年4月14日的Landsat 8 OLI影像。所選影像質(zhì)量完好,云量較少,時相一致,能保證研究結(jié)果具有可比性。采用ENVI 5.3軟件對遙感影像數(shù)據(jù)進(jìn)行預(yù)處理,主要包括輻射定標(biāo)、大氣校正、波段組合以及圖像配準(zhǔn)和裁剪等。非遙感數(shù)據(jù)包括地貌分布圖、土壤植被覆蓋圖、湖北省省市縣行政區(qū)劃圖和清江統(tǒng)計(jì)文本資料。

    2.2 RSEI的構(gòu)建

    綠度、濕度、熱度和干度能直觀反映生態(tài)條件優(yōu)劣,常用于評價(jià)生態(tài)環(huán)境。采用改進(jìn)型土壤植被指數(shù)(modified adjusted soil vegetation index,MASVI)、濕度分量(WET)、地表溫度(land surface temperature,LST)和裸土指數(shù)(normalized difference soil index,NDSI)分別代表綠度、濕度、熱度和干度來綜合評價(jià)清江水電開發(fā)區(qū)生態(tài)環(huán)境?;跀?shù)據(jù)標(biāo)準(zhǔn)化和主成分變換方法將IRSE轉(zhuǎn)換為各項(xiàng)指標(biāo)的函數(shù)[10-11]為

    IRSE=f(IMASV,TWE,TLS,INDS)=

    (1)

    式(1)中,IMASV為綠度;TWE為濕度;TLS為熱度;INDS為干度;Mu、Wu、Lu和Nu分別為柵格圖像中第u個柵格MASVI、WET、LST和NDSI值;ωM、ωW、ωL和ωN分別為MASVI、WET、LST和NDSI權(quán)重;n為柵格數(shù)。

    2.2.1分量指標(biāo)的選取

    (1)綠度指標(biāo)。植被對區(qū)域生態(tài)環(huán)境狀況指示極為敏感。改進(jìn)型土壤植被指數(shù)通過降低土壤背景影響來提高植被指數(shù)的敏感度。由于研究區(qū)主要地貌分布類型為山地,故選用IMASV可以更好地反映此類型用地的綠度情況,其計(jì)算公式[12]為

    (2)

    式(2)中,ρNIR和ρred分別為TM和OLI近紅外和紅光波段的反射率。

    (2)濕度指標(biāo)。土壤濕度是用于土壤退化等生態(tài)環(huán)境研究的重要指標(biāo)。遙感纓帽變換提取的濕度分量反映了地表土壤和植被的濕度狀況?;赥M和OLI反射率數(shù)據(jù)計(jì)算濕度分量(WET,TWE)的公式[13-14]分別為

    TWE,TM=0.031 5ρblue+0.202 1ρgreen+0.310 2ρred+

    0.159 4ρNIR-0.680 6ρSWIR1-0.610 9ρSWIR2,

    (3)

    TWE,OLI=0.151 1ρblue+0.197 2ρgreen+0.328 3ρred+

    0.340 7ρNIR-0.711 7ρSWIR1-0.455 9ρSWIR2。

    (4)

    式(3)~(4)中,ρblue、ρgreen、ρred、ρNIR、ρSWIR1和ρSWIR2分別為TM和OLI藍(lán)、綠、紅、近紅外、短波紅外1和短波紅外2波段的反射率。

    (3)熱度指標(biāo)。地表溫度與植被的生長與分布、地表水資源蒸發(fā)循環(huán)密切相關(guān),是反映地表環(huán)境的重要參數(shù)之一。由于NASA的陸地過程分布式數(shù)據(jù)檔案中心MODIS現(xiàn)有產(chǎn)品地表溫度分辨率為1 000 m,為確保用GIS對分量指標(biāo)進(jìn)行疊加分析的精度一致,可通過地表熱輻射傳輸方程反演獲取?;赥M6提取地表真實(shí)溫度Ts的計(jì)算公式[15]為

    Ts=Tb/[1+(λTb/α)lnε6]。

    (5)

    式(5)中,Tb為傳感器處溫度,K;中心波長λ為11.435 μm;α取1.438×10-2mK;ε6為基于TM6波段的地表比輻射率。

    同時考察Landsat 8熱紅外波段10和11,由于波段11對大氣剖面的誤差敏感[16],故基于波段10來反演地表溫度,通過普朗克定律反函數(shù)求得Ts:

    Ts=K2/ln [K1/B10(Ts)+1]。

    (6)

    式(6)中,B10(Ts)為溫度Ts的黑體在熱紅外波段10的輻射亮度值;K1和K2取值分別為774.89 W·m-2·sr-1·μm-1和1 321.08 K。地表比輻射率ε6計(jì)算方法見文獻(xiàn)[17]。

    (4)干度指標(biāo)。自然裸土和人為施工建筑用地均會造成地表“干化”,干度指標(biāo)由裸土指數(shù)(soil index,SI,IS)[18]和建筑指數(shù)(index-based built-up index,IBI,IIB)[19]組成,記為裸土指數(shù)(NDSI,INDS)[11],其計(jì)算公式為

    INDS=(IS+IIB)/2。

    (7)

    為了消除量綱差異對RSEI的影響,對各指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理:

    In=(I-Imin)/(Imax-Imin)。

    (8)

    式(8)中,In為標(biāo)準(zhǔn)化后指標(biāo);I為原始指標(biāo);Imax和Imin為指標(biāo)最大值和最小值。

    2.2.2主成分分析

    生態(tài)環(huán)境評價(jià)的關(guān)鍵是將改進(jìn)型土壤植被指數(shù)、濕度分量、地表溫度和裸土指數(shù)轉(zhuǎn)化為綜合評價(jià)指標(biāo)。由于研究區(qū)存在清江、長江大面積水體,為了避免大片水域影響主成分分析(principal component analysis,PCA)的荷載分布,使?jié)穸戎笜?biāo)能夠刻畫地表真實(shí)濕度,采用MNDWI(modified normalized difference water index)指數(shù)掩膜掉水體信息[20]。此次研究將各標(biāo)準(zhǔn)化分量指標(biāo)進(jìn)行波段合成后,基于水體掩膜信息進(jìn)行主成分分析,對清江IRSE進(jìn)行轉(zhuǎn)換,其表達(dá)式為

    (9)

    (10)

    (11)

    式(11)中,IMASV和TWE變量系數(shù)為正,說明其對生態(tài)環(huán)境起積極作用;而TLS和INDS變量系數(shù)為負(fù),說明其對生態(tài)環(huán)境起消極作用。利用式(8)對3個時期IRSE進(jìn)行標(biāo)準(zhǔn)化,采用ArcGIS 10.5軟件中相等間隔(0.2)的重分類函數(shù),將IRSE劃分為差(0~0.2)、較差(>0.2~0.4)、中等(>0.4~0.6)、良(>0.6~0.8)和優(yōu)(>0.8~1.0)5個生態(tài)等級[21],提取研究區(qū)生態(tài)環(huán)境變化分級圖。

    3 結(jié)果與分析

    3.1 模型檢驗(yàn)

    為了判定RSEI的綜合代表程度,通過統(tǒng)計(jì)各年份4個分量指標(biāo)和RSEI的相關(guān)系數(shù)以及各指標(biāo)自身之間的相關(guān)性來檢驗(yàn)?zāi)P瓦m宜性(表1)。干度分量平均相關(guān)度最高,3 a平均相關(guān)度達(dá)0.685,綠度、濕度和熱度的3 a平均相關(guān)度則分別為0.626、0.642和0.505。而RSEI與各指標(biāo)的平均相關(guān)度均大于0.750,3 a平均相關(guān)度達(dá)0.799,比4個指標(biāo)的平均值(0.614)高30.1%。因此,RSEI與各指標(biāo)間的相關(guān)度更高,可反映清江水電開發(fā)區(qū)生態(tài)環(huán)境狀況。

    表1 各年份清江水電開發(fā)區(qū)生態(tài)環(huán)境評價(jià)分量指標(biāo)與遙感生態(tài)指數(shù)(RSEI)相關(guān)性Table 1 Correlation between the indicators of ecological environment assessment and remote sensing ecological index (RSEI)of Qingjiang hydropower development zone in each year

    MASVI為綠度,WET為濕度,LST為熱度,NDSI為干度。

    3.2 水電開發(fā)區(qū)域生態(tài)環(huán)境評價(jià)分析

    統(tǒng)計(jì)分析1987—2015年4個分量指標(biāo)和RSEI(表2),清江水電開發(fā)區(qū)1987、2004和2015年RSEI均值分別為0.443、0.664和0.515,呈逐階段先上升后下降趨勢,變化幅度分別為49.9%和22.4%。其中,代表生態(tài)變好的綠度和濕度均值呈先上升后下降趨勢;代表生態(tài)變差的熱度和干度均值呈先降低后上升趨勢,單一指標(biāo)的綜合作用與RSEI基本吻合。

    表2 清江水電開發(fā)區(qū)生態(tài)環(huán)境評價(jià)各項(xiàng)分量指標(biāo)和遙感生態(tài)指數(shù)(RSEI)統(tǒng)計(jì)值Table 2 Four statistical indicators of ecological environment assessment and the statistical value of remote sensing ecological index (RSEI)in Qingjiang hydropower development zone

    MASVI為綠度,WET為濕度,LST為熱度,NDSI為干度。

    通過分析各分量指標(biāo)對PCA載荷值的貢獻(xiàn)發(fā)現(xiàn),在1987—2004年間,綠度和濕度分量對PCA貢獻(xiàn)度較大,由于森林覆蓋面積與水域面積擴(kuò)大,MASVI與WET分別上升41.19%和16.27%,清江中下游地區(qū)環(huán)境質(zhì)量明顯提高;而2004—2015年間,綠度和干度分量對PCA貢獻(xiàn)度較大,植被覆蓋度減少及裸土干化使MASVI下降15.51%,NDSI上升9.40%,研究區(qū)生態(tài)質(zhì)量開始惡化。在水電梯級開發(fā)的兩個時間段,綠度貢獻(xiàn)度最大,說明植被覆蓋度對生態(tài)環(huán)境的影響最大,成為影響清江中下游地區(qū)生態(tài)環(huán)境的決定因素。

    為了進(jìn)一步分析RSEI體現(xiàn)的生態(tài)效應(yīng),統(tǒng)計(jì)1987、2004和2015年各生態(tài)等級面積和占比(表3)。結(jié)果表明,1987—2015年間,清江中下游地區(qū)生態(tài)等級變化明顯,1987年生態(tài)等級為中等以下面積所占比例為86.64%,2004年為37.89%,2015年為72.11%,說明1987—2004年間生態(tài)環(huán)境呈良好的發(fā)展態(tài)勢,但2004—2015年間生態(tài)環(huán)境呈下降趨勢。因此,在水電梯級項(xiàng)目全面建成之后,對清江中下游地區(qū)環(huán)境進(jìn)行整治尤為必要。

    為了探究生態(tài)環(huán)境的空間變化特征,對水電梯級開發(fā)各階段生態(tài)環(huán)境進(jìn)行空間分級(圖2)。在水電梯級開發(fā)起始年,清江中下游地區(qū)生態(tài)環(huán)境整體處于中等水平,環(huán)境惡劣地區(qū)(差、較差等級)主要集中在西部裸土嚴(yán)重的喀斯特丘陵、山地以及中部、西南部大起伏山地區(qū),環(huán)境良好(良、優(yōu)等級)地區(qū)主要分布在東部林地、農(nóng)業(yè)發(fā)達(dá)的耕地地區(qū)及清江流域沿岸地帶。2004年清江中下游地區(qū)生態(tài)環(huán)境整體處于良水平,環(huán)境惡劣地區(qū)主要集中在隔河巖、高壩洲兩個水電梯級建設(shè)項(xiàng)目的下游地區(qū)及長江沿岸丘陵區(qū)。2015年清江中下游地區(qū)生態(tài)環(huán)境整體處于中等水平,環(huán)境惡劣地區(qū)以水布埡水電站為起始點(diǎn),圍繞清江干、支流呈放射狀分布。

    表3 1987、2004和2015年清江水電開發(fā)區(qū)各生態(tài)等級面積及占比Table 3 Area and proportion of different ecological grades in Qingjiang hydropower development zone in 1987,2004 and 2015

    3.3 清江水電開發(fā)區(qū)域生態(tài)環(huán)境動態(tài)變化分析

    為了對生態(tài)狀況時空變化特征進(jìn)行分析,在5個RSEI生態(tài)等級的基礎(chǔ)上,利用ArcGIS 10.5軟件對清江中下游地區(qū)生態(tài)狀況差值變化進(jìn)行檢測(表4)。1987—2004年間,該區(qū)域生態(tài)條件變化顯著,生態(tài)條件變差區(qū)域面積為989.34 km2,約占總面積的11.50%,主要分布在清江沿岸溝壑以及新建隔河巖、高壩洲下游區(qū)域〔圖3(a)中紅色圖斑〕;生態(tài)條件改善區(qū)域面積達(dá)6 015.50 km2,占69.92%,主要分布在低山丘陵、喀斯特山地,土地利用以常綠闊葉林、常綠針葉林和針闊葉混交林為主。

    圖2 1987、2004和2015年清江中下游地區(qū)遙感生態(tài)指數(shù)(RSEI)等級分布Fig.2 Distribution of remote sensing ecological index (RSEI)in the middle and lower reaches of Qingjiang River in 1987,2004 and 2015

    表4 1987、2004和2015年清江中下游地區(qū)生態(tài)狀況差值變化Table 4 Changes of ecological conditions in the middle and lower reaches of the Qingjiang River in 1987,2004 and 2015

    隨著下游梯級大壩的建成,日調(diào)節(jié)型隔河巖水庫正常蓄水位為200 m,庫容30.18億m3,年調(diào)節(jié)型高壩洲水庫正常蓄水位為80 m,庫容4.03億m3。更多的水資源被保留在峽谷型庫容內(nèi),為農(nóng)業(yè)、漁業(yè)和林業(yè)等的發(fā)展以及林地、草地等的恢復(fù)提供更多水源。2003年植被混凝土護(hù)坡綠化技術(shù)的成功應(yīng)用[22],使得2004年綠度指數(shù)(0.593)和土壤濕度(0.636)為3 a最高,促進(jìn)了該地區(qū)農(nóng)業(yè)升級改造的發(fā)展及植被的恢復(fù)。

    圖3 1987、2004和2015年清江中下游地區(qū)遙感生態(tài)指數(shù)(RSEI)變化Fig.3 Changes of remote sensing ecological index (RSEI)in the middle and lower reaches of the Qingjiang River in 1987,2004 and 2015

    2015年清江中下游地區(qū)生態(tài)環(huán)境整體呈退化趨勢,生態(tài)環(huán)境惡化區(qū)域面積達(dá)4 991.12 km2,占總面積的58.01%,而生態(tài)環(huán)境改善區(qū)域面積為618.48 km2,僅占7.19%。生態(tài)條件改善區(qū)域主要分布在東部低山丘陵區(qū)和平坦平原區(qū)〔圖3(b)中綠色圖斑〕。生態(tài)環(huán)境退化區(qū)域廣泛分布在以新建水布埡樞紐為中心的上下游區(qū)域,梯級水電站建設(shè)對生態(tài)環(huán)境的影響作用是疊加的,包括時間和空間的累積。工程區(qū)層間剪切帶發(fā)育,斷裂構(gòu)造較發(fā)育,有巖溶洞穴及巖溶通道,地質(zhì)條件較為復(fù)雜。巖性為灰?guī)r、泥灰?guī)r和頁巖等軟硬互層,土層淺薄,易受侵蝕。

    水電梯級開發(fā)活動引起水土流失,表現(xiàn)為擾動面積大、強(qiáng)度高,建設(shè)周期長。水土保持生態(tài)效應(yīng)評價(jià)主要包括水、土和植被3個方面,即保水、保土和植被恢復(fù)[23]。水電梯級項(xiàng)目全面建成后,土壤濕度(0.477)為3 a最低,比2004年降低25.00%。土壤含水量是陸表生態(tài)系統(tǒng)水循環(huán)的重要組成部分,影響土壤性質(zhì)和植被生長。土壤抗沖蝕能力差等自然因素和主體工程施工活動等人為因素,壓占土地、開挖表土和擾動土地范圍較大,均不利于植被恢復(fù)[24]。2015年干度指數(shù)(0.512)比2004年升高9.40%,綠度指數(shù)(0.501)比2014年降低15.51%,過度開發(fā)使清江中下游地區(qū)進(jìn)入快速調(diào)整和不穩(wěn)定階段。因此,在水電梯級建設(shè)的同時,有效防止水土流失造成的地表干化和提高植被覆蓋度成為生態(tài)環(huán)境恢復(fù)的關(guān)鍵。

    1987—2015年清江中下游地區(qū)各縣市生態(tài)環(huán)境空間特征穩(wěn)定比例均超過30%(表5),從空間〔圖3(c)〕上來看,生態(tài)環(huán)境惡化區(qū)域面積比例相對較高(26.89%)的為宜都市,生態(tài)環(huán)境改善區(qū)域面積比例較高的為長陽自治縣和巴東縣,比例分別為47.77%和51.46%。在清江下游水電開發(fā)階段(1987—2004年),宜都市生態(tài)環(huán)境變化劇烈,環(huán)境退化區(qū)域面積比例達(dá)26.05%,長陽自治縣、五峰自治縣和巴東縣環(huán)境改善區(qū)域面積比例均>70.00%。隔河巖和高壩洲水利樞紐工程興建引發(fā)的水庫蓄水及部分河段的擴(kuò)寬工程使得水域面積明顯增加,植被環(huán)境恢復(fù)與建設(shè)的關(guān)鍵就是解決水源問題?;谶b感估算植被覆蓋度[25]的結(jié)果也表明,1987年清江中下游地區(qū)林地覆蓋率僅為21.50%,2004年上升為57.66%,植被覆蓋面積呈現(xiàn)增長趨勢。

    水電梯級項(xiàng)目全面建成后,長陽自治縣、五峰自治縣和巴東縣呈現(xiàn)不同程度的環(huán)境惡化趨勢,宜都市惡化區(qū)域面積比例雖有所減緩,但也超過環(huán)境恢復(fù)速率。其中,巴東縣中部及東南部、五峰自治縣北部以及長陽自治縣西北部和東南部差等級增加幅度最大,以隔河巖為中心的長陽自治縣東南部和宜都市東南部環(huán)境有所改善,說明清江中下游地區(qū)受水電梯級建設(shè)人類活動干擾的影響,生態(tài)環(huán)境等級變化明顯。大壩及水庫群的建設(shè)占用大量土地資源,植被覆蓋度降低至34.29%,比2004年降低23.37百分點(diǎn)。由于裸土造成的地表干化程度(祼土指數(shù))提高9.40%,流域梯級的過度開發(fā)會破壞天然的地表地貌,導(dǎo)致不同程度的水土流失。2015年清江中下游地區(qū)地表溫度比2004年升高16.95%,水庫蓄水后,庫區(qū)水體面積及覆被變化,使其與空氣之間的能量交換方式和強(qiáng)度發(fā)生改變,導(dǎo)致氣溫變化,影響原來穩(wěn)定的群落結(jié)構(gòu)完整性[26]。因此,在發(fā)揮水電梯級項(xiàng)目建設(shè)巨大經(jīng)濟(jì)效益和生態(tài)價(jià)值的同時,提高植被覆蓋度及防范因過度開發(fā)引起的水土流失、地表干化已成為流域治理和生態(tài)修復(fù)的重點(diǎn)。

    表5 1987、2004和2015年清江中下游地區(qū)各縣市生態(tài)環(huán)境變化面積比例Table 5 Proportion of ecological environment changes in counties and cities alongside the middle and lower reaches of the Qingjiang River in 1987,2004 and 2015

    4 討論與結(jié)論

    傳統(tǒng)的生態(tài)環(huán)境評價(jià)方法中所選河流水質(zhì)、水文流量和底棲動物群落等數(shù)據(jù),需人工采集,獲取耗時、耗力且具有滯后性。完全基于遙感的RSEI,采用主成分分析方法確定權(quán)重,克服了以往生態(tài)評價(jià)指標(biāo)難以空間化和權(quán)重確定具有主觀性的缺點(diǎn),通過遙感空間變量的解譯和耦合,能定量揭示區(qū)域生態(tài)環(huán)境的空間分布與變化,為水電開發(fā)區(qū)生態(tài)環(huán)境提供實(shí)時、客觀的監(jiān)測技術(shù)。

    依據(jù)水電梯級開發(fā)的時間順序,選取Landsat 5和Landsat 8遙感數(shù)據(jù),通過建立新型RSEI,對清江水電開發(fā)區(qū)1987—2015年間生態(tài)環(huán)境時空變化進(jìn)行分析,以期為區(qū)域生態(tài)環(huán)境保護(hù)與治理提供科學(xué)依據(jù),得出以下結(jié)論。

    (1)清江中下游地區(qū)1987、2004和2015年RSEI均值分別為0.443、0.664和0.515,呈先上升后下降趨勢。這是因?yàn)樵谙掠翁菁壦娬窘ǔ珊?水庫蓄水、部分河段的擴(kuò)寬工程和防護(hù)坡綠化技術(shù),促進(jìn)了植被的恢復(fù),生態(tài)環(huán)境得到一定改善;水電梯級項(xiàng)目全面建成后,過度開發(fā)擾動大量地表,原生植被遭到破壞,土壤抗侵蝕能力變差,水土流失加劇,導(dǎo)致生態(tài)環(huán)境惡化。

    (2)清江中下游地區(qū)受水電梯級建設(shè)活動干擾的影響,生態(tài)環(huán)境等級變化明顯。從空間上來看,1987—2015年間,生態(tài)環(huán)境改善區(qū)域主要分布于低山丘陵區(qū)和平坦平原區(qū),生態(tài)環(huán)境退化區(qū)域主要分布在水利樞紐中心及其附屬庫區(qū)、沿岸喀斯特山地以及清江匯入長江口的丘陵區(qū)。在清江下游水電開發(fā)階段,宜都市生態(tài)環(huán)境變化劇烈,環(huán)境退化比例最大;梯級水電站全面建成后,長陽自治縣、五峰自治縣和巴東縣生態(tài)環(huán)境均呈現(xiàn)惡化趨勢,宜都市生態(tài)環(huán)境惡化趨勢雖有所減緩,但超過環(huán)境恢復(fù)速率。

    (3)為了提高生態(tài)承載力,加強(qiáng)水利樞紐沿岸地帶水土流失及地表干化的治理力度,各縣市應(yīng)該積極發(fā)展水土保持型生態(tài)農(nóng)業(yè),保護(hù)好天然林和退耕林,合理調(diào)度水資源,提高植被覆蓋度,以改善清江中下游地區(qū)生態(tài)環(huán)境。

    猜你喜歡
    綠度清江梯級
    寧夏草地綠度時空變化圖譜及驅(qū)動因素分析
    全球森林綠度異常對干旱事件的響應(yīng)
    2003~2018年干旱對中國森林綠度異常的影響
    梯級水電站多目標(biāo)聯(lián)合經(jīng)濟(jì)運(yùn)行初探
    清江引
    影劇新作(2018年1期)2018-05-26 09:00:52
    淺談唐卡中綠度母的造型特點(diǎn)
    魚躍清江 廣場舞
    文化交流(2017年8期)2017-09-14 22:03:27
    同飲清江水 共護(hù)母親河——首個“清江保護(hù)日”在長陽舉行
    跨流域梯級電站群遠(yuǎn)程集控管理模式的探索與實(shí)踐
    跨流域梯級電站群集中調(diào)控的探索與思考
    国内精品久久久久精免费| 婷婷精品国产亚洲av| 精品国产亚洲在线| 国产成人系列免费观看| 久久久精品欧美日韩精品| 成年人黄色毛片网站| 精品久久久久久成人av| av电影中文网址| 国产精品99久久99久久久不卡| 亚洲中文字幕日韩| 成人av一区二区三区在线看| 午夜福利在线在线| 午夜福利18| 一区二区日韩欧美中文字幕| 亚洲成av片中文字幕在线观看| 成在线人永久免费视频| 一个人观看的视频www高清免费观看 | 国产成人系列免费观看| 一进一出好大好爽视频| 变态另类丝袜制服| 久久久久久国产a免费观看| 亚洲欧美一区二区三区黑人| 久久九九热精品免费| 中文字幕精品亚洲无线码一区 | 亚洲人成网站在线播放欧美日韩| 97碰自拍视频| 午夜免费观看网址| 欧美黑人精品巨大| 又紧又爽又黄一区二区| 少妇粗大呻吟视频| 国产在线观看jvid| 久久国产精品男人的天堂亚洲| 88av欧美| 国产精品久久视频播放| 亚洲av成人不卡在线观看播放网| 婷婷六月久久综合丁香| or卡值多少钱| 国产高清有码在线观看视频 | 男女床上黄色一级片免费看| 亚洲狠狠婷婷综合久久图片| av超薄肉色丝袜交足视频| 久久香蕉激情| 一卡2卡三卡四卡精品乱码亚洲| 最新美女视频免费是黄的| 精品国产国语对白av| 国产黄片美女视频| 老司机靠b影院| 一a级毛片在线观看| 美女高潮到喷水免费观看| 日韩高清综合在线| 国产精品美女特级片免费视频播放器 | 老熟妇仑乱视频hdxx| 丰满的人妻完整版| 精品国产乱码久久久久久男人| 岛国视频午夜一区免费看| 亚洲精品中文字幕一二三四区| 满18在线观看网站| 给我免费播放毛片高清在线观看| 免费在线观看黄色视频的| 成熟少妇高潮喷水视频| 一区二区三区国产精品乱码| 白带黄色成豆腐渣| 国产黄片美女视频| 国产精品久久久久久亚洲av鲁大| 久久99热这里只有精品18| 琪琪午夜伦伦电影理论片6080| 99久久无色码亚洲精品果冻| 久久久精品国产亚洲av高清涩受| 中文字幕最新亚洲高清| 久久精品亚洲精品国产色婷小说| 日本黄色视频三级网站网址| 搡老岳熟女国产| 欧美激情极品国产一区二区三区| 色播亚洲综合网| 亚洲成人久久爱视频| 日本在线视频免费播放| 韩国av一区二区三区四区| 麻豆久久精品国产亚洲av| 亚洲第一av免费看| 久久久久久国产a免费观看| 夜夜爽天天搞| 久久性视频一级片| 狂野欧美激情性xxxx| 免费观看精品视频网站| 女性被躁到高潮视频| 身体一侧抽搐| 少妇被粗大的猛进出69影院| 黄色视频,在线免费观看| 桃红色精品国产亚洲av| 久久久久久亚洲精品国产蜜桃av| 久久伊人香网站| 欧美一级a爱片免费观看看 | 窝窝影院91人妻| 国产成+人综合+亚洲专区| 免费在线观看视频国产中文字幕亚洲| 久久久久久人人人人人| 国内揄拍国产精品人妻在线 | 十分钟在线观看高清视频www| 人妻丰满熟妇av一区二区三区| 女性被躁到高潮视频| 波多野结衣高清无吗| 岛国视频午夜一区免费看| 99国产精品一区二区蜜桃av| 精品一区二区三区四区五区乱码| 天堂影院成人在线观看| av电影中文网址| 国内少妇人妻偷人精品xxx网站 | 又黄又粗又硬又大视频| 免费在线观看影片大全网站| 久久99热这里只有精品18| 一夜夜www| 黄色片一级片一级黄色片| 亚洲欧美日韩高清在线视频| 午夜福利18| 人成视频在线观看免费观看| 中文在线观看免费www的网站 | 黄色视频不卡| 精品欧美一区二区三区在线| 999精品在线视频| 成人亚洲精品一区在线观看| 99久久综合精品五月天人人| 欧美日韩黄片免| 国产一区二区在线av高清观看| 非洲黑人性xxxx精品又粗又长| 国产伦在线观看视频一区| 午夜福利欧美成人| 亚洲五月婷婷丁香| 中文资源天堂在线| 国产91精品成人一区二区三区| 欧美成狂野欧美在线观看| 一区二区三区国产精品乱码| 少妇被粗大的猛进出69影院| 久久久精品国产亚洲av高清涩受| 久久国产精品影院| 国产一区二区激情短视频| 两人在一起打扑克的视频| 女性被躁到高潮视频| 香蕉av资源在线| 国产欧美日韩一区二区精品| 日本熟妇午夜| 亚洲精品美女久久久久99蜜臀| 久久这里只有精品19| 老司机深夜福利视频在线观看| avwww免费| 一区二区三区精品91| 脱女人内裤的视频| 国内揄拍国产精品人妻在线 | 一边摸一边抽搐一进一小说| 啦啦啦观看免费观看视频高清| 日韩欧美国产一区二区入口| 中亚洲国语对白在线视频| 黄片播放在线免费| 丁香欧美五月| 给我免费播放毛片高清在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲精品美女久久久久99蜜臀| 99久久精品国产亚洲精品| 美女国产高潮福利片在线看| 在线天堂中文资源库| av片东京热男人的天堂| 国产精品亚洲一级av第二区| 国产午夜精品久久久久久| 亚洲一区高清亚洲精品| 搡老熟女国产l中国老女人| 欧美日韩瑟瑟在线播放| 99久久综合精品五月天人人| 在线av久久热| 我的亚洲天堂| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| av在线天堂中文字幕| 午夜激情福利司机影院| 国产激情欧美一区二区| 老汉色∧v一级毛片| 亚洲男人天堂网一区| 一本一本综合久久| 最近最新中文字幕大全电影3 | 长腿黑丝高跟| av福利片在线| 中文字幕久久专区| 12—13女人毛片做爰片一| 久久久国产精品麻豆| 97超级碰碰碰精品色视频在线观看| 欧美日韩瑟瑟在线播放| 国产不卡一卡二| 琪琪午夜伦伦电影理论片6080| a在线观看视频网站| 久久久久国产一级毛片高清牌| 窝窝影院91人妻| 妹子高潮喷水视频| 老司机靠b影院| 国产精品九九99| 高清毛片免费观看视频网站| 久久久久久亚洲精品国产蜜桃av| 黑人操中国人逼视频| 欧美+亚洲+日韩+国产| 变态另类丝袜制服| 最近在线观看免费完整版| 国产精品一区二区免费欧美| 免费无遮挡裸体视频| 国产又色又爽无遮挡免费看| 18禁美女被吸乳视频| www国产在线视频色| 99国产精品一区二区三区| 人人妻人人澡欧美一区二区| 日韩视频一区二区在线观看| 免费在线观看成人毛片| 在线国产一区二区在线| 亚洲av成人不卡在线观看播放网| 久久中文字幕一级| 精品日产1卡2卡| 亚洲自偷自拍图片 自拍| 国产又色又爽无遮挡免费看| 亚洲精华国产精华精| 亚洲天堂国产精品一区在线| 日本黄色视频三级网站网址| 国产精品一区二区三区四区久久 | 国产三级在线视频| 视频区欧美日本亚洲| 99精品久久久久人妻精品| 色尼玛亚洲综合影院| 欧美不卡视频在线免费观看 | 国产一区在线观看成人免费| 丝袜人妻中文字幕| 夜夜躁狠狠躁天天躁| 男人舔奶头视频| 亚洲精品美女久久久久99蜜臀| 非洲黑人性xxxx精品又粗又长| 国产精品乱码一区二三区的特点| 亚洲成人精品中文字幕电影| 国产精品久久视频播放| 亚洲国产精品久久男人天堂| 亚洲 欧美 日韩 在线 免费| 国产一级毛片七仙女欲春2 | 狂野欧美激情性xxxx| 一本大道久久a久久精品| 美女大奶头视频| 亚洲av第一区精品v没综合| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看 | 伦理电影免费视频| 欧美三级亚洲精品| 美女高潮到喷水免费观看| 日韩一卡2卡3卡4卡2021年| 精品国产国语对白av| 日韩欧美 国产精品| 天堂√8在线中文| 99久久国产精品久久久| 熟妇人妻久久中文字幕3abv| av在线播放免费不卡| 免费看十八禁软件| 亚洲欧洲精品一区二区精品久久久| 亚洲三区欧美一区| АⅤ资源中文在线天堂| 亚洲成a人片在线一区二区| 欧美乱色亚洲激情| 久久精品国产亚洲av香蕉五月| 给我免费播放毛片高清在线观看| 久久国产亚洲av麻豆专区| 精品人妻1区二区| 国产精品久久久久久人妻精品电影| 久久热在线av| 日本黄色视频三级网站网址| 久久久久久大精品| 精品午夜福利视频在线观看一区| 色综合婷婷激情| 中文亚洲av片在线观看爽| 久久天堂一区二区三区四区| 在线观看日韩欧美| 国产精品一区二区免费欧美| 黄色成人免费大全| 亚洲激情在线av| 欧美日韩一级在线毛片| 欧美绝顶高潮抽搐喷水| 精品久久久久久久末码| 国内久久婷婷六月综合欲色啪| 欧美性长视频在线观看| www日本在线高清视频| 国内揄拍国产精品人妻在线 | 久久久精品国产亚洲av高清涩受| 看免费av毛片| 久久久久久九九精品二区国产 | 男女之事视频高清在线观看| 成人三级做爰电影| 少妇 在线观看| 色综合亚洲欧美另类图片| 丝袜人妻中文字幕| 夜夜躁狠狠躁天天躁| 老熟妇仑乱视频hdxx| 天天添夜夜摸| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品一卡2卡三卡4卡5卡| 久久天堂一区二区三区四区| 男人操女人黄网站| 久久精品国产综合久久久| 日韩精品免费视频一区二区三区| a在线观看视频网站| 丁香六月欧美| 男人舔女人下体高潮全视频| 91国产中文字幕| 中文资源天堂在线| 丰满的人妻完整版| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 岛国视频午夜一区免费看| 国产精品久久视频播放| 国产av又大| 一级片免费观看大全| 国产av一区在线观看免费| 一区二区三区高清视频在线| videosex国产| 免费在线观看亚洲国产| 成人永久免费在线观看视频| 国产不卡一卡二| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 男女做爰动态图高潮gif福利片| 夜夜夜夜夜久久久久| 日本成人三级电影网站| 日韩欧美免费精品| www.熟女人妻精品国产| 国内少妇人妻偷人精品xxx网站 | 黑人欧美特级aaaaaa片| 国内少妇人妻偷人精品xxx网站 | 欧美亚洲日本最大视频资源| 亚洲九九香蕉| 亚洲第一电影网av| 久久婷婷成人综合色麻豆| 国产精品精品国产色婷婷| 欧美另类亚洲清纯唯美| 久久久久久免费高清国产稀缺| 亚洲性夜色夜夜综合| 在线播放国产精品三级| 国产高清视频在线播放一区| 国产黄色小视频在线观看| 久久热在线av| 午夜福利18| 嫁个100分男人电影在线观看| 亚洲av成人一区二区三| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 国产伦在线观看视频一区| 美女高潮喷水抽搐中文字幕| xxxwww97欧美| 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| 成人三级黄色视频| 精品欧美一区二区三区在线| 久久狼人影院| 成人18禁在线播放| 久久精品国产清高在天天线| 观看免费一级毛片| 国产av一区二区精品久久| 夜夜夜夜夜久久久久| 老司机在亚洲福利影院| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| av天堂在线播放| 禁无遮挡网站| 国产成人精品久久二区二区免费| 久久人人精品亚洲av| 午夜免费激情av| 成人永久免费在线观看视频| 88av欧美| 亚洲精品国产一区二区精华液| 亚洲全国av大片| 一边摸一边抽搐一进一小说| 91九色精品人成在线观看| 欧美激情 高清一区二区三区| 十分钟在线观看高清视频www| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 两个人看的免费小视频| 波多野结衣av一区二区av| 久久九九热精品免费| 色老头精品视频在线观看| 国产一区二区在线av高清观看| 日韩欧美在线二视频| 最好的美女福利视频网| 成人手机av| 亚洲午夜精品一区,二区,三区| 90打野战视频偷拍视频| 久久热在线av| 精华霜和精华液先用哪个| 日韩有码中文字幕| 黑人巨大精品欧美一区二区mp4| 一夜夜www| 精品卡一卡二卡四卡免费| 成在线人永久免费视频| 757午夜福利合集在线观看| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区mp4| 男女之事视频高清在线观看| 一区二区三区国产精品乱码| 国产精品久久电影中文字幕| 亚洲精品国产区一区二| 亚洲精品国产精品久久久不卡| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 国产精品久久视频播放| 中文字幕高清在线视频| 久久天堂一区二区三区四区| 免费看日本二区| 97碰自拍视频| 香蕉国产在线看| 岛国在线观看网站| 色综合亚洲欧美另类图片| 美女 人体艺术 gogo| 久热爱精品视频在线9| 久久中文字幕人妻熟女| 亚洲欧美日韩高清在线视频| 他把我摸到了高潮在线观看| 国产单亲对白刺激| 久久精品人妻少妇| 国产色视频综合| 老熟妇乱子伦视频在线观看| 黄色视频,在线免费观看| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 熟女少妇亚洲综合色aaa.| 日韩欧美国产一区二区入口| 成年人黄色毛片网站| 日本成人三级电影网站| 国产一区二区激情短视频| 最近最新中文字幕大全电影3 | www.999成人在线观看| www.熟女人妻精品国产| 午夜福利欧美成人| e午夜精品久久久久久久| 欧美三级亚洲精品| 一级片免费观看大全| 不卡av一区二区三区| 免费高清在线观看日韩| 成人永久免费在线观看视频| a级毛片在线看网站| 精品国产超薄肉色丝袜足j| 精品熟女少妇八av免费久了| 国产精品1区2区在线观看.| 99精品在免费线老司机午夜| 午夜福利一区二区在线看| 少妇的丰满在线观看| av天堂在线播放| 黄色片一级片一级黄色片| 国产一区二区三区视频了| 这个男人来自地球电影免费观看| 白带黄色成豆腐渣| 欧美日韩乱码在线| 久久精品国产综合久久久| 国产亚洲欧美在线一区二区| 黄色a级毛片大全视频| 18禁国产床啪视频网站| 法律面前人人平等表现在哪些方面| 国产私拍福利视频在线观看| 亚洲av中文字字幕乱码综合 | 精品卡一卡二卡四卡免费| 日本a在线网址| 又紧又爽又黄一区二区| 可以在线观看的亚洲视频| 色在线成人网| 中文字幕最新亚洲高清| videosex国产| 亚洲av五月六月丁香网| 欧美色视频一区免费| 波多野结衣高清无吗| 婷婷丁香在线五月| 亚洲精品色激情综合| 色播在线永久视频| 一个人观看的视频www高清免费观看 | 国产麻豆成人av免费视频| 精品久久久久久,| 黄色丝袜av网址大全| 亚洲精品粉嫩美女一区| 熟妇人妻久久中文字幕3abv| 桃红色精品国产亚洲av| 国产精品亚洲一级av第二区| 亚洲五月天丁香| 国产精品免费视频内射| 美女午夜性视频免费| 国内毛片毛片毛片毛片毛片| 夜夜躁狠狠躁天天躁| 最好的美女福利视频网| 国产极品粉嫩免费观看在线| 欧美亚洲日本最大视频资源| 真人一进一出gif抽搐免费| 亚洲精品在线美女| 大香蕉久久成人网| 亚洲熟女毛片儿| 色精品久久人妻99蜜桃| 白带黄色成豆腐渣| cao死你这个sao货| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频| 少妇裸体淫交视频免费看高清 | 在线看三级毛片| 最好的美女福利视频网| 成年免费大片在线观看| 91成人精品电影| 少妇的丰满在线观看| 亚洲国产精品sss在线观看| 精品欧美一区二区三区在线| www日本黄色视频网| 午夜福利成人在线免费观看| avwww免费| www.精华液| 久久精品国产清高在天天线| 99久久99久久久精品蜜桃| 老鸭窝网址在线观看| 久久久久久九九精品二区国产 | 午夜激情av网站| 欧美日韩一级在线毛片| 国产精品香港三级国产av潘金莲| 亚洲 欧美一区二区三区| 亚洲欧美一区二区三区黑人| 老汉色∧v一级毛片| 免费看日本二区| 欧美黑人欧美精品刺激| 精品国产国语对白av| 欧美成人一区二区免费高清观看 | 亚洲精品国产区一区二| 久久精品国产99精品国产亚洲性色| 18禁国产床啪视频网站| 国产伦在线观看视频一区| 嫁个100分男人电影在线观看| 午夜免费成人在线视频| 国产aⅴ精品一区二区三区波| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 在线看三级毛片| 亚洲激情在线av| 中文资源天堂在线| 大香蕉久久成人网| 欧美av亚洲av综合av国产av| 国内精品久久久久精免费| 久久精品国产亚洲av高清一级| 欧美亚洲日本最大视频资源| 国产亚洲欧美精品永久| 日韩三级视频一区二区三区| 白带黄色成豆腐渣| 三级毛片av免费| 日韩欧美一区二区三区在线观看| 精品国产国语对白av| 十八禁人妻一区二区| 国产野战对白在线观看| 久久久久久久久免费视频了| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品在线美女| 精品午夜福利视频在线观看一区| 变态另类丝袜制服| 国产精品一区二区精品视频观看| x7x7x7水蜜桃| 欧美黄色片欧美黄色片| 亚洲最大成人中文| 免费在线观看黄色视频的| 久久精品国产亚洲av香蕉五月| 午夜日韩欧美国产| 久久久久国内视频| 精品久久久久久,| 人妻丰满熟妇av一区二区三区| 国产成人精品无人区| 999久久久国产精品视频| 国产精品永久免费网站| 亚洲av电影在线进入| 亚洲五月色婷婷综合| 成人国语在线视频| АⅤ资源中文在线天堂| 国产精品99久久99久久久不卡| 国内精品久久久久久久电影| 亚洲,欧美精品.| 无限看片的www在线观看| 麻豆国产av国片精品| 亚洲精品久久成人aⅴ小说| 午夜日韩欧美国产| 91国产中文字幕| 欧美在线黄色| 波多野结衣巨乳人妻| 国产成人影院久久av| 日韩成人在线观看一区二区三区| 91老司机精品| 亚洲国产精品成人综合色| 精品国产超薄肉色丝袜足j| 美女 人体艺术 gogo| 精品国产一区二区三区四区第35| 久久精品国产清高在天天线| 国产又黄又爽又无遮挡在线| 亚洲精品色激情综合| 久久久久久久久中文| 免费女性裸体啪啪无遮挡网站| 婷婷精品国产亚洲av| 好男人在线观看高清免费视频 | cao死你这个sao货| 精品福利观看| 一区二区三区精品91| 精品电影一区二区在线| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费| 国产国语露脸激情在线看| 精品久久久久久,| svipshipincom国产片| 国产欧美日韩精品亚洲av| 一本综合久久免费| 色精品久久人妻99蜜桃| 免费在线观看影片大全网站| 国产97色在线日韩免费| 久久精品国产亚洲av香蕉五月| 精品国产乱子伦一区二区三区| 日韩视频一区二区在线观看| 色老头精品视频在线观看| 在线观看免费日韩欧美大片| 亚洲国产欧洲综合997久久, | 亚洲精品美女久久av网站| 午夜日韩欧美国产| 久久亚洲真实| 99久久久亚洲精品蜜臀av| 成人国语在线视频|