李慧璇,黃林峰,胡 鋒,宋善海,劉 蕓

(貴州省生態(tài)氣象和衛(wèi)星遙感中心,貴州 貴陽(yáng) 550002)

0 引言

隨著我國(guó)對(duì)生態(tài)文明建設(shè)的高度重視,全社會(huì)對(duì)生態(tài)環(huán)境保護(hù)的認(rèn)識(shí)不斷加深,諸多學(xué)者也在該研究領(lǐng)域進(jìn)行了積極的探索。進(jìn)一步加大區(qū)域生態(tài)環(huán)境全方位監(jiān)測(cè)和綜合治理力度,利用衛(wèi)星遙感技術(shù)對(duì)區(qū)域生態(tài)環(huán)境進(jìn)行全天候、精準(zhǔn)、快速和廣視域的環(huán)境監(jiān)測(cè)是各監(jiān)管部門(mén)區(qū)域生態(tài)保護(hù)和生態(tài)文明建設(shè)的主要措施和手段,也是當(dāng)今社會(huì)的重點(diǎn)研究課題[1-3]。

區(qū)域生態(tài)環(huán)境狀況遙感監(jiān)測(cè)是生態(tài)環(huán)境治理的一項(xiàng)重要工作。目前,各界學(xué)者對(duì)生態(tài)環(huán)境的評(píng)價(jià)技術(shù)方法已從20世紀(jì)60年代初對(duì)生態(tài)環(huán)境的單要素調(diào)查發(fā)展為綜合多源數(shù)據(jù)、多環(huán)境要素指標(biāo)評(píng)價(jià)[4-9]。此外,生態(tài)環(huán)境遙感監(jiān)測(cè)已成為區(qū)域生態(tài)環(huán)境評(píng)價(jià)的重要手段之一,被廣泛應(yīng)用于城市、森林、沙漠、草地和濕地等區(qū)域,取得顯著成效[10-14]。赤水河流域作為貴州省重要的生態(tài)水源涵養(yǎng)地、生態(tài)屏障和經(jīng)濟(jì)地帶,由于流域內(nèi)喀斯特地貌廣泛分布,生態(tài)敏感性和脆弱性并存,其生態(tài)環(huán)境保護(hù)問(wèn)題長(zhǎng)久以來(lái)備受關(guān)注。隨著國(guó)產(chǎn)高分辨率衛(wèi)星遙感數(shù)據(jù)自給率大幅提升,為貴州赤水河流域土地利用分類信息提取提供了可靠的數(shù)據(jù)支撐。且利用衛(wèi)星影像對(duì)生態(tài)環(huán)境進(jìn)行監(jiān)測(cè)的技術(shù)方法在我國(guó)環(huán)境治理中已經(jīng)越來(lái)越成熟,能更好地對(duì)流域生態(tài)環(huán)境狀況進(jìn)行高效判讀[15-19]。

以高分一號(hào)和高分二號(hào)衛(wèi)星影像為主要數(shù)據(jù)源,通過(guò)對(duì)影像進(jìn)行輻射定標(biāo)、大氣校正、正射校正、融合、鑲嵌和波段運(yùn)算等預(yù)處理,運(yùn)用面向?qū)ο髨D像分類法中基于規(guī)則的面向?qū)ο蠓诸惙ㄅc基于樣本的面向?qū)ο蠓诸惙ㄟM(jìn)行流域土地利用分類的提取,并依據(jù)國(guó)家生態(tài)環(huán)境部發(fā)布的《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范(HJ192-2015)》[20](以下簡(jiǎn)稱《規(guī)范》),利用生物豐度、植被覆蓋、水網(wǎng)密度、土地脅迫和污染負(fù)荷等多個(gè)指數(shù)對(duì)研究區(qū)的生態(tài)環(huán)境狀況進(jìn)行綜合評(píng)估,以期為赤水河流域生態(tài)環(huán)境保護(hù)和綜合治理提供科學(xué)依據(jù)和數(shù)據(jù)支撐。

1 研究區(qū)概況

赤水河屬中國(guó)長(zhǎng)江上游的一級(jí)支流,全流域覆蓋云南、貴州和四川3省的4個(gè)地區(qū)級(jí)和13個(gè)縣級(jí)行政單位。貴州境內(nèi)赤水河流域地理位置為26°49′～28°54′N(xiāo)、104°09′～107°10′E,主要分布于播州、七星關(guān)、大方、金沙、仁懷、習(xí)水、赤水、桐梓等多個(gè)縣區(qū)。

本次監(jiān)測(cè)流域面積13 926.09 km2,流域海拔在211～2213 m之間,高度落差>2000 m。研究范圍、地理位置以及區(qū)域DEM如圖1所示。

圖1 赤水河流域監(jiān)測(cè)范圍及地理位置圖

流域上游所在的畢節(jié)市七星關(guān)區(qū)、大方縣和金沙縣平均海拔最高,低海拔的區(qū)域基本都是河谷地帶。流域下游赤水市平均海拔最低,呈現(xiàn)出較大面積的平緩河谷。流域整體水能豐富,年均降水1090 mm,降雨多集中在6—9月,降水量約占全年的75%。降水量在上游七星關(guān)區(qū)偏低,在下游赤水市偏高,分別為974.7 mm和1289.8 mm。

2 資料方法

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

監(jiān)測(cè)區(qū)資源和環(huán)境數(shù)據(jù)主要來(lái)源于中國(guó)環(huán)境監(jiān)測(cè)總站的全國(guó)城市空氣質(zhì)量實(shí)時(shí)發(fā)布平臺(tái)、貴州省宏觀經(jīng)濟(jì)數(shù)據(jù)庫(kù)和《貴州統(tǒng)計(jì)年鑒2021》。流域監(jiān)測(cè)矢量范圍來(lái)源于貴州省環(huán)境科學(xué)研究設(shè)計(jì)院。遙感影像數(shù)據(jù)為國(guó)家資源衛(wèi)星中心提供的貴州省區(qū)域高分一號(hào)和高分二號(hào)影像資料,融合后影像分辨率分別為1 m、2 m。此外,影像處理過(guò)程中采用的DEM數(shù)據(jù)空間分辨率為12.5 m。

通過(guò)中國(guó)資源衛(wèi)星與應(yīng)用中心獲得原始影像后,利用Titan批量處理軟件進(jìn)行數(shù)據(jù)預(yù)處理,處理過(guò)程包括影像解壓、選片、建立金字塔、區(qū)域網(wǎng)平差像點(diǎn)量測(cè)、輻射定標(biāo)、正射校正、融合、波段運(yùn)算、勻色、鑲嵌、裁剪等。

2.2 研究區(qū)生態(tài)環(huán)境狀況評(píng)價(jià)方法

2.2.1 指標(biāo)體系構(gòu)建及計(jì)算方法 根據(jù)《規(guī)范》和流域特點(diǎn),構(gòu)建研究區(qū)生態(tài)環(huán)境狀況的綜合評(píng)價(jià)指標(biāo)體系,最終利用生態(tài)環(huán)境EI指數(shù)來(lái)反映研究區(qū)生態(tài)環(huán)境現(xiàn)狀,數(shù)值范圍取值0～100。文章所構(gòu)建的5個(gè)指標(biāo)體系共有和具體計(jì)算方法如下。

生態(tài)環(huán)境狀況指數(shù)(EI)=0.35×生物豐度指數(shù)+0.25×植被覆蓋指數(shù)+0.15×水網(wǎng)密度指數(shù)+0.15×(100-土地脅迫指數(shù))+0.1×(100-污染負(fù)荷指數(shù))+環(huán)境限制指數(shù)。各項(xiàng)指標(biāo)的權(quán)重以及歸一化系數(shù)如表1所示。

表1 流域生態(tài)環(huán)境狀況指標(biāo)體系及指標(biāo)權(quán)重

2.2.2 生態(tài)環(huán)境狀況分級(jí) 在分析過(guò)程中,結(jié)合生態(tài)環(huán)境狀況指數(shù)計(jì)算結(jié)果,參考現(xiàn)有研究[22],將生態(tài)環(huán)境分為優(yōu)、良、一般、差4個(gè)等級(jí)。4個(gè)等級(jí)的分級(jí)標(biāo)準(zhǔn)分別為EI≥75、55≤EI<75、35≤EI<55、EI<35。

3 結(jié)果及分析

3.1 流域地類提取及結(jié)果分析

得到研究區(qū)全域預(yù)處理結(jié)果影像后,進(jìn)一步借助ENVI影像處理平臺(tái),采用面向?qū)ο蠓?根據(jù)高分辨率影像的光譜統(tǒng)計(jì)特征、紋理特征和形狀大小特征等和實(shí)地采集的地類樣本,對(duì)流域的林地、草地、水域濕地、耕地、建設(shè)用地和未利用地等6大地類進(jìn)行地類分割和提取。

人工修正和檢驗(yàn)過(guò)程中,結(jié)合前期無(wú)人機(jī)外業(yè)采集的多個(gè)植被、建筑和水體等樣本,使用人機(jī)判讀法進(jìn)一步完成流域空間分布信息提取,得到的地類分布結(jié)果如圖4所示。

圖4 2021年赤水河流域土地利用信息

根據(jù)提取得到的圖斑類型,利用ArcGIS軟件進(jìn)行統(tǒng)計(jì),得到流域監(jiān)測(cè)范圍內(nèi)各地類的土地面積及土地利用類型占比(圖略)。

赤水河流域貴州段的土地利用類型主要以林地和耕地為主,占總面積的93.7%。其中,林地占63.28%,包括灌木林、有林地、疏林地和其他林地,耕地占30.26%。建設(shè)用地、草地、未利用地次之,分別占總面積的3.94%、1.45%、0.007%。

3.2 評(píng)價(jià)指數(shù)計(jì)算及結(jié)果分析

3.2.1 生物豐度指數(shù) 本次赤水河流域監(jiān)測(cè)區(qū)內(nèi)生物豐度指數(shù)主要通過(guò)計(jì)算林地、草地、水域濕地、耕地、建設(shè)用地和未利用地等不同生態(tài)因子所占區(qū)域的面積比重,并作歸一化處理得到。根據(jù)地類面積提取結(jié)果,計(jì)算得到流域監(jiān)測(cè)區(qū)內(nèi)生物豐度指數(shù)為134.14。

根據(jù)劉華等[23,24]研究,把生境質(zhì)量IRHQ最終的評(píng)價(jià)結(jié)果劃分為優(yōu)、良、一般、差、極差5個(gè)等級(jí)。5個(gè)等級(jí)的分級(jí)標(biāo)準(zhǔn)分別為IRHQ≥75、60≤IRHQ﹤75、45≤IRHQ﹤60、30≤IRHQ﹤45、IRHQ<30。

結(jié)果可見(jiàn),赤水河流域監(jiān)測(cè)區(qū)內(nèi)生境質(zhì)量等級(jí)為優(yōu),生物多樣性豐富,以流域西北部的赤水市為最佳。

3.2.2 植被覆蓋指數(shù) 結(jié)合高分遙感全年影像,計(jì)算監(jiān)測(cè)區(qū)NDVI值,計(jì)算結(jié)果值中,巖石或裸土用0值表示,有植被覆蓋的,數(shù)值越大,表示該地的植被覆蓋情況越好。計(jì)算監(jiān)測(cè)區(qū)NDVI結(jié)果如圖5所示。其中,最大值為0.67,平均值為0.42。

圖5 赤水河流域(貴州段)2021年NDVI計(jì)算結(jié)果

結(jié)果可見(jiàn),流域植被覆蓋度高,固土護(hù)坡、涵養(yǎng)水源能力較強(qiáng),流域西北及西南部植被覆蓋能力最強(qiáng),森林多保持原始狀態(tài),人類活動(dòng)干擾少。

3.2.3 水網(wǎng)密度指數(shù) 結(jié)合高分遙感影像和統(tǒng)計(jì)資料可知,赤水河流域貴州段多年平均降雨量約1100 mm,多年平均蒸發(fā)量為456 mm,多年平均徑流量56.23億m3,多年平均地下水資源量15.72億m3,多年平均地表水資源量約為70.76億m3。此外,赤水河流域河網(wǎng)密布,三省范圍內(nèi)面積超過(guò)1000 km2的河流便有5條之多,其中二道河、桐梓河和習(xí)水河3條支流分布于貴州段的上中下游。

提取水域信息后,利用ArcGIS計(jì)算得到監(jiān)測(cè)區(qū)的大小河流的總長(zhǎng)度達(dá)2403.58 km,其中赤水河長(zhǎng)度約300 km,水域面積127.06 km2。從研究區(qū)總面積基數(shù)大的角度看,河網(wǎng)整體面積占比少,計(jì)算得到水網(wǎng)密度指數(shù)結(jié)果為6.81,但從圖6可見(jiàn),流域河網(wǎng)密布,水系相對(duì)發(fā)達(dá)。

圖6 赤水河流域(貴州段)河流分布

3.2.4 土地脅迫指數(shù) 根據(jù)《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》(SL190-2007),將土壤侵蝕強(qiáng)度分為輕度、中度、強(qiáng)烈、極強(qiáng)烈和劇烈5級(jí)(表2),把強(qiáng)烈、極強(qiáng)烈和劇烈侵蝕認(rèn)為是重度侵蝕,利用土地類型提取信息、監(jiān)測(cè)區(qū)坡度信息和植被覆蓋指數(shù),完成土壤侵蝕面積的計(jì)算。

表2 土地侵蝕強(qiáng)度分級(jí)

計(jì)算流域坡度,結(jié)果如圖7所示。

圖7 赤水河流域(貴州段)坡度信息提取結(jié)果

赤水河流域的地帶性土壤類型以石灰土、紫色土和粗骨土為主,土質(zhì)整體透水性雖強(qiáng),但由于土層較薄,加之地形起伏大,流域整體相對(duì)容易被侵蝕。利用AcrGIS的空間分析模塊得到監(jiān)測(cè)區(qū)中度侵蝕面積4025.39 km2,占總面積28.9%,重度侵蝕面積2915.7 km2,占總面積20.9%。通過(guò)計(jì)算,土地脅迫指數(shù)計(jì)算結(jié)果為35.34。

土地脅迫指數(shù)與生態(tài)環(huán)境狀況評(píng)價(jià)指數(shù)呈負(fù)相關(guān),計(jì)算得到的土地脅迫指數(shù)值越大,代表流域內(nèi)的土地質(zhì)量受脅迫程度越大,對(duì)應(yīng)的流域生態(tài)環(huán)境狀況越差。

3.2.5 污染負(fù)荷指數(shù) 近年來(lái),政府監(jiān)管部門(mén)對(duì)赤水河流域水污染、土壤污染、固體廢物污染和空氣質(zhì)量等的監(jiān)管力度不斷加大,且數(shù)據(jù)開(kāi)放渠道越來(lái)越規(guī)范化。據(jù)地方相關(guān)資料,赤水河流域釀酒、煤礦等主要工業(yè)企業(yè)和現(xiàn)有規(guī)?；男笄蒺B(yǎng)殖業(yè)相對(duì)較多,存在一定數(shù)量的點(diǎn)污染源。以區(qū)縣為單位,監(jiān)測(cè)區(qū)總磷、SO2排放量、氨氮排放量和固體廢物排放量均可從地方公開(kāi)數(shù)據(jù)獲取。但在計(jì)算過(guò)程中,因檢測(cè)區(qū)域范圍較大,受專業(yè)資料限制,監(jiān)測(cè)區(qū)污染數(shù)據(jù)無(wú)法全面獲取,根據(jù)可獲取的數(shù)據(jù),得到的流域污染負(fù)荷指數(shù)計(jì)算結(jié)果值為0.53。

4 結(jié)論與討論

綜合各指標(biāo)計(jì)算結(jié)果,得到赤水河流域生態(tài)環(huán)境狀況EI指數(shù)為68.15,生態(tài)環(huán)境狀況評(píng)價(jià)為良。可見(jiàn),流域整體森林保護(hù)能力和水源涵養(yǎng)能力都比較強(qiáng),主要表現(xiàn)為以下4點(diǎn):生境質(zhì)量?jī)?yōu),生物多樣性豐富;植被覆蓋率高,固土護(hù)坡、水源涵養(yǎng)能力強(qiáng);水資源豐富,生態(tài)水景好;生態(tài)區(qū)監(jiān)測(cè)區(qū)污染防治成效顯著,流域生態(tài)環(huán)境良好。

文章依托現(xiàn)有遙感在生態(tài)保護(hù)中的監(jiān)測(cè)能力,收集GF系列衛(wèi)星遙感數(shù)據(jù),結(jié)合赤水河的區(qū)域特點(diǎn)及環(huán)境保護(hù)要點(diǎn)確定了基于高分衛(wèi)星數(shù)據(jù)的流域地類分類提取方法。采用面向?qū)ο蠓诸惙椒?通過(guò)構(gòu)建研究區(qū)光譜特征、紋理特征及形狀特征,完成研究區(qū)地類分類提取和地類面積計(jì)算。進(jìn)一步利用地面數(shù)據(jù)和生態(tài)環(huán)境狀況指數(shù)EI客觀、真實(shí)地反映反映了貴州赤水河流域生態(tài)環(huán)境的整體狀況,并制作了相應(yīng)的專題圖。該方法可有效提升流域生態(tài)環(huán)境感知和生態(tài)環(huán)境數(shù)據(jù)標(biāo)準(zhǔn)化水平,也能客觀衡量生態(tài)環(huán)境保護(hù)工作的狀況和水平。但是,貴州赤水河流域的生態(tài)環(huán)境是復(fù)雜的系統(tǒng)性問(wèn)題,流域生態(tài)環(huán)境遙感監(jiān)測(cè)也面臨著諸多難題,比如EI指數(shù)綜合評(píng)價(jià)中污染負(fù)荷指數(shù)的計(jì)算容易受監(jiān)測(cè)范圍和專業(yè)資料限制,無(wú)法全面獲取數(shù)據(jù)。此外,由于數(shù)據(jù)時(shí)間域較窄,本文只能對(duì)流域生態(tài)環(huán)境現(xiàn)狀進(jìn)行評(píng)價(jià),后續(xù)研究將進(jìn)行多方位、多時(shí)相、長(zhǎng)時(shí)序、多數(shù)據(jù)源的綜合監(jiān)測(cè),關(guān)注流域生態(tài)環(huán)境的變化規(guī)律和區(qū)域經(jīng)濟(jì)的和諧可持續(xù)發(fā)展。