隋建紅 隋巖剛 李超

收稿日期：20231031；修訂日期：20240213；編輯：曹麗麗

作者簡介：隋建紅（1979—），女，山東諸城人，正高級工程師，主要從事環(huán)境地質、水文地質、遙感技術等工作；Email：sdusui@163.com? *

作者簡介：隋巖剛（1981—），男，山東乳山人，高級工程師，主要從事地質環(huán)境、水文地質等工作;Email：suiyangang@163.com

摘要：采用最大值合成法、像元二分模型和趨勢分析法，以NASA發(fā)布的2017—2022年MODIS NDVI 16d最大值合成數據為基礎數據，通過植被覆蓋度變化趨勢對山東泰山區(qū)域生態(tài)修復工程整體效果進行研究。結果表明：2017—2022年泰山區(qū)域NDVI最大值和平均值總體呈波動上升趨勢，最大值增長率為0.0055/a，平均值增長率為0.0049/a。山東省泰山區(qū)域Slope平均值為0.00478，生態(tài)修復工程絕大多數所處像元植被覆蓋度變化率大于0，說明整個區(qū)域植被覆蓋度有所改善，工程有明顯效果。

關鍵詞：山水林田湖草;生態(tài)修復工程;植被覆蓋度;像元二分法;趨勢分析法；泰山區(qū)域；山東省

中圖分類號：X37；F301??? 文獻標識碼：A??? doi：10.12128/j.issn.16726979.2024.06.007

引文格式：隋建紅，隋巖剛，李超.基于植被覆蓋度變化的山東泰山區(qū)域生態(tài)修復工程效果研究［J］.山東國土資源，2024，40（6）：5156.SUI Jianhong， SUI Yan'gang，LI Chao. Study on the Effect of Ecological Restoration Project in Taishan Area in Shandong Province Based on the Change of Vegetation Coverage［J］.Shandong Land and Resources，2024，40（6）：5156.

0? 引言

評價一個區(qū)域的生態(tài)環(huán)境變化情況，多采用區(qū)域內植被覆蓋狀態(tài)來衡量，而評價植被覆蓋狀態(tài)的重要指標是植被覆蓋度［13］（Fractional vegetation cover， FVC）。一般將植被冠層在地面上垂直投影面積之和與區(qū)域總面積的比值定義為植被覆蓋度［46］。利用遙感數據研究植被覆蓋度變化趨勢，是當前區(qū)域生態(tài)建設影響評價的重要手段。主要采用歸一化植被指數（NDVI）、像元二分模型應用和趨勢分析法等研究方法。常規(guī)工程驗收僅對單個工程生態(tài)修復效果［79］進行調查評價，不能對某區(qū)域整體生態(tài)修復效果進行評價。本文通過研究區(qū)工程建設前后植被覆蓋度變化評價整體區(qū)域生態(tài)修復效果。

1? 研究區(qū)概況與數據

1.1? 研究區(qū)概況

2017年12月，山東省泰山區(qū)域入選國家第二批山水林田湖草生態(tài)保護修復工程試點范圍［10］。山東省泰山區(qū)域生態(tài)保護修復工程涵蓋省會濟南及泰安2市的14個縣（市、區(qū)），布設工程類型有地質環(huán)境、生物多樣性、水環(huán)境、土地整治和監(jiān)管能力建設等5大類［1112］。2018年工程全面實施，截至2022年3月，總投資約159億元，修復礦山和整治土地2.8萬hm2，治理采煤塌陷地9605hm2。2022年底進行了整體工程驗收，區(qū)域生態(tài)功能系統(tǒng)得到提升。

1.2? 數據源與預處理

數據源為美國國家航空航天局（NASA，網址https：//ladsweb.modaps.eosdis.nasa.gov/）發(fā)布的2017年1月—2022年12月經過幾何精校正、輻射校正處理的16d最大值合成數據集MOD13Q1（分辨率250m）。

預處理采用的軟件工具為MRT批處理工具和ArcGIS軟件。首先對MOD13Q1數據進行投影和格式轉化處理，然后利用泰山區(qū)域邊界矢量數據對處理后的數據進行裁剪，再次對研究區(qū)數據進行統(tǒng)計和分析，最后采用最大化合成法（MVC）合成年NDVI數據。MOD13Q1原始數據和研究對象NDVI取值范圍不同，兩者之間轉換關系式為NDVI=DNⅹ0.0001［1314］。

2? 研究方法

2.1? 植被覆蓋度估算

像元二分模型是一種簡單實用的遙感估算模型［1516］，其計算公式如式（1）：

FVC=NDVI-NDVIsNDVIv-NDVIs（1）

式中：FVC代表計算單元的植被覆蓋度；NDVI代表計算單元的NDVI值；NDVIs代表純植物單元的NDVI值；NDVIv代表完全無植被覆蓋單元的NDVI值。

裸地表面NDVIs理論上應接近0且不易變化，但受眾多因素影響，NDVIs會隨著空間變化而變化，其變化范圍一般在0.1～0.2。同時，NDVIv值隨著植被類型不同和季節(jié)變化而變化，因此計算植被覆蓋度時，即使是對同一景影像，NDVIs和NDVIv也不能取固定值，故選取評價區(qū)域內NDVImax和NDVImin代表NDVIv和NDVIs，計算公式如式（2）：

FVC=0??????????? NDVI≤NDVIminNDVI-NDVIminNDVImax-NDVImin? NDVImin

式中：NDVImax為評價區(qū)域內累計概率為95%的NDVI值；NDVImin為評價區(qū)域內累計概率為5%的NDVI值。

2.2? 趨勢分析法

趨勢分析法以時間為變量對研究數據進行線性回歸分析，以便研究數據的變化趨勢［1718］。計算方法如式（3）：

Slope=n∑ni=1（i×NDVIi）-∑ni=1i×∑ni=1NDVIin∑ni=1i2-（∑ni=1i）2? （3）

式中：Slope為像元回歸方程的斜率;NDVIi為第1年的NDVI值;n為研究的時間長度，視自身情況而定。

當Slope>0時，表示該像元NDVI為增加趨勢;當Slope=0，表示該像元NDVI基本不變;當Slope<0時，表示該像元NDVI為減少趨勢［1920］。

3? 結果和分析

3.1? 植被覆蓋時間變化特征

3.1.1? 整個區(qū)域變化特征

2017—2022年，泰山區(qū)域NDVI最大值和平均值均總體呈波動上升趨勢，最大值增長率為0.0055/a，平均值增長率為0.0049/a，變化趨勢見圖1。2019年生態(tài)修復工程處于初期，效果尚未顯現，而區(qū)域建設仍處于發(fā)展狀態(tài)，故最大值和平均值均在該年度達到最低值。

3.1.2? 不同等級區(qū)域變化特征

參照2008年水利部頒布的《土壤侵蝕分類分級標準》［21］將劃分植被覆蓋度分為5個等級［22］，<30%為植被低覆蓋度；30%～45%為植被中低覆蓋度；45%～60%為植被中等覆蓋度；60%～75%為植被中高覆蓋度；>75%為植被高覆蓋度。

2017—2022年山東泰山區(qū)域各不同等級覆蓋度占比變化趨勢見圖2?？梢钥闯?，低覆蓋度占地呈先降后增的變化趨勢，中等覆蓋度、中低覆蓋度和中高覆蓋度占地均呈先增后降的變化趨勢，高覆蓋度占地呈逐漸升高的變化趨勢。

3.2? 植被覆蓋空間分布特征

3.2.1? 年植被覆蓋度空間分布特征

2017—2022年，泰山區(qū)域低覆蓋度區(qū)主要分布于區(qū)域北部、西部和東南部的城市建設區(qū)和東平湖、大汶河等河流、湖泊區(qū)；高覆蓋區(qū)主要分布于區(qū)域中部和南部；逐年等級植被覆蓋度分布特征變化較?。▓D3）。

3.2.2? 年植被覆蓋度變化率空間分布特征

采用趨勢分析法對2017—2022年山東泰山區(qū)域植被覆蓋度以時間作為變量因子進行趨勢分析（圖4，表1），分析結果表明：6年來，部分區(qū)域植被覆蓋度處于增加趨勢，植被覆蓋度顯著改善區(qū)域面積為810.17km2，占比達5.2%；90.5%面積占比區(qū)域植被覆蓋度基本保持不變。

3.3? 生態(tài)修復工程與覆蓋度變化率關系

將實施的山東泰山區(qū)域生態(tài)修復工程（水環(huán)境治理類除外）與2017—2022年植被覆蓋度變化率疊合，獲取到其所處像元的變化率，進行分析（圖5）。

由圖5可以看出，2017年以來實施的生態(tài)修復工程絕大多數所處像元植被覆蓋度變化率大于0，且均不小于0.05（顯著退化），平均值為0.014，進一步說明了實施工程有效性，顯著改善區(qū)域與生態(tài)修復工程中植樹、綠化區(qū)也基本對應。

4? 結語

采用最大值合成法、像元二分模型和趨勢分析法等研究方法，以NASA發(fā)布的2017—2022年MODIS NDVI 16d最大值合成數據為基礎數據，對山東泰山區(qū)域近6年來（生態(tài)保護修復工程實施前后）植被覆蓋度時間和空間變化趨勢進行分析，以便驗證山東省泰山區(qū)域生態(tài)修復工程整體效果。

（1）2017—2022年，泰山區(qū)域NDVI最大值和平均值總體均呈上升趨勢，最大值增長率為0.0055/a，平均值增長率為0.0049/a。

（2）2017—2022年，山東泰山區(qū)域低覆蓋度占地呈先降后增的變化趨勢，中等覆蓋度、中低覆蓋度和中高覆蓋度占地均呈先增后降的變化趨勢，高覆蓋度占地呈逐漸升高的變化趨勢。

（3）部分區(qū)域植被覆蓋度處于增加趨勢，植被覆蓋度明顯改善區(qū)域面積占比達5.2%；90.5%面積占比區(qū)域植被覆蓋度基本保持不變。

（4）山東泰山區(qū)域Slope平均值為0.00478，生態(tài)修復工程絕大多數所處像元植被覆蓋度變化率大于0，且均不小于0.05（顯著退化），說明整個區(qū)域植被覆蓋度呈現改善趨勢，工程有明顯效果。

參考文獻：

［1］? 王柳，段英.利用遙感影像進行植被分布分析［J］.測繪與空間地理信息，2012，35（3）：140142.

［2］? 李世偉，尼加提·穆合塔爾，周迪，等.濟南市植被覆蓋時空變化特征研究［J］.綠色科技，2020，14：813.

［3］? 金凱，王飛，韓劍橋，等.1982—2015年中國氣候變化和人類活動對植被NDVI變化的影響［J］.地理學報，2020，75（5）：961974.

［4］? 白欣，賀炳彥，高婷，等.基于GEE云平臺的榆林市植被覆蓋度時空變化研究［J］.湖北農業(yè)科學，2022，61（16）：7075.

［5］? 李玉虎，滕佳華，馬文濤.京津冀地區(qū)植被覆蓋度時空演變特征分析［J］.北京測繪，2019，33（11）：13021307.

［6］? 張志強，劉歡，左其亭，等.2000—2019年黃河流域植被覆蓋度時空變化［J］.資源科學，2021，43（4）：849858.

［7］? 郝啟勇，王忠，胡玉祿，等.沂蒙山區(qū)域山水林田湖草一體化保護和修復工程規(guī)劃實踐［J］.山東國土資源，2022，38（5）：4146.

［8］? 徐風，焦玉國，唐麗偉，等.泰安市山水林田湖草生態(tài)修復區(qū)廢棄露天礦山治理模式與技術體系研究［J］.山東國土資源，2022，38（6）：6371.

［9］? 邊微.泰山區(qū)域山水林田湖草保護修復試點：泰安市土地綜合整治模式探索［J］.山東國土資源，2022，38（3）：7075.

［10］? 山東省人民政府辦公廳.關于推進泰山區(qū)域山水林田湖草生態(tài)保護修復工程的實施意見［N］.山東省人民政府公報，2018（16）：1420.

［11］? 徐風，唐麗偉，魏凱，等.山東省泰安市山水林田湖草生態(tài)保護修復工程平臺設計與實現［J］.山東國土資源，2023，39（12）：3948.

［12］? 于楚榕.泰山區(qū)域山水林田湖草生態(tài)保護修復工程調查研究［J］.低碳世界，2020（2）：216217.

［13］? 徐勇，奚硯濤，劉欣婷.2000—2013年江蘇省植被覆蓋度時空變化特征［J］.西南農業(yè)學報，2016，29（2）：437444.

［14］? 徐勇.植被覆蓋動態(tài)變化及其與氣候因子的時空響應特征研究［D］.徐州：中國礦業(yè)大學，2015：120.

［15］? 王愛蕓，趙志芳.基于遙感的植被覆蓋度估算方法研究進展［J］.綠色科技，2015（3）：1015.

［16］? 鐘琪，胡晉山，康建榮.基于像元二分法的大寧礦區(qū)植被覆蓋度研究［J］.金屬礦山，2021（11）：197203.

［17］? 王曉蕾，石守海，陳江朝霞.黃河流域植被覆蓋度變化及驅動因素［J］.中國環(huán)境科學，2022，42（11）：53585368.

［18］? 馮莉莉，何貞銘，劉學鋒，等.基于MODISNDVI數據的吉林省植被覆蓋度及其時空動態(tài)變化［J］.中國科學院大學學報，2014，31（4）：492499.

［19］? 徐勇，奚硯濤.江蘇省植被覆蓋動態(tài)變化及其與氣候因子的關系［J］.水土保持通報，2015（4）：195201.

［20］? 郭帥，裴艷茜，胡勝，等.黃河流域植被指數對氣候變化的響應及其與水沙變化的關系［J］.水土保持通報，2020，40（3）：17.

［21］? 中華人民共和國水利部.土壤侵蝕分類分級標準：SL1902007［S］.北京：中國水利水電出版社，2008.

［22］? 李鈺溦，賈坤，魏香琴，等.中國北方地區(qū)植被覆蓋度遙感估算及其變化分析［J］.國土資源遙感，2015（2）：112117.

Study on the Effect of Ecological Restoration Project in Taishan Area in Shandong Province Based on the Change of Vegetation Coverage

SUI Jianhong， SUI Yan'gang，LI Chao

（Shandong Institute of Coalfield Geological Engineering and Exploration， Shandong Ji'nan 250100， China）

Abstract：According to the change trend of vegetation coverage， the overall effect of ecological restoration project in Taishan area in Shandong province has been studied by maximum synthesis method， trend analysis method and pixel dichotomy model. The study data comes from MODIS NDVI data products released by NASA from 2017 to 2022. As showed by the results， from 2017 to 2022， the maximum value and average value of NDVI in Taishan region show a fluctuating upward trend， with the maximum growth rate of 0.0055/a and the average growth rate of 0.0049/a. The average slope value is 0.00478. The vegetation coverage change rate of most of the restoration projects is greater than 0. It is indicated that the vegetation coverage of the whole region shows an improvement trend， and the project has obvious effect.

Key words：Mountain，water，forest，field，lake and grass; ecological environment restoration project; vegetation coverage; pixel dichotomy; trend analysis; Taishan region; Shandong province