黃艷利，郭亞超，齊文躍，李俊孟，王佳奇，歐陽神央，吳來偉

(1.中國礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116；2.中國礦業(yè)大學(xué) 煤炭資源與安全開采國家重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221116；3.新疆工程學(xué)院 礦業(yè)工程與地質(zhì)學(xué)院，新疆 烏魯木齊 830023；4. 燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島 066004)

西部是我國內(nèi)陸向西開放和建設(shè)“絲綢之路經(jīng)濟(jì)帶”的關(guān)鍵節(jié)點(diǎn)，西部5省(陜蒙甘寧新)已成為我國煤炭資源的主產(chǎn)區(qū)，產(chǎn)量約占全國總產(chǎn)量的70%[1-2]。西部礦區(qū)煤炭資源多，具有煤層埋藏淺、厚度大、地質(zhì)構(gòu)造簡單等特點(diǎn)，礦區(qū)多屬水資源短缺、地表生態(tài)脆弱的干旱-半干旱區(qū)域[3-5]，礦區(qū)高強(qiáng)度開采極易引起地表塌陷[6]、地下水位下降[7]、植被退化[8]、土地荒漠化等生態(tài)環(huán)境損傷問題[9-10]。植被作為表征礦區(qū)生態(tài)環(huán)境變化的綜合指示器[11-12]，可以較好地反映采動(dòng)損傷對(duì)生態(tài)環(huán)境變化的影響。

近年來，關(guān)于采礦活動(dòng)對(duì)植被發(fā)育狀態(tài)的影響國內(nèi)外學(xué)者做了大量的研究工作，并取得了許多有價(jià)值的研究成果。ERENER[13]分析了不同植被指數(shù)在土耳其礦區(qū)植被生長狀況監(jiān)測(cè)的差異；雷少剛等[14-17]以神東礦區(qū)為研究對(duì)象，研究了采區(qū)與非采區(qū)植被的差異，揭示了采煤塌陷對(duì)植物葉綠素、植物高光譜特征、植物光合作用速率以及礦區(qū)植被覆蓋度的影響規(guī)律；李長春等[18-19]通過遙感技術(shù)，研究了典型露天礦區(qū)植被覆蓋情況，對(duì)礦區(qū)環(huán)境進(jìn)行了評(píng)價(jià)；李茂林等[20-21]研究了內(nèi)蒙古典型干旱半干旱礦區(qū)植被覆蓋度變化趨勢(shì)，分析了地形地貌、采礦活動(dòng)、土地利用類型與植被覆蓋度變化的關(guān)系；吳立新等[22]通過線性回歸的方法研究了神東礦區(qū)的土壤沙化現(xiàn)象和植被覆蓋變化特征；YANG等[23]通過野外調(diào)查與遙感相結(jié)合的方法，研究了補(bǔ)連河流域植被群落結(jié)構(gòu)及其對(duì)地下開采的響應(yīng)，發(fā)現(xiàn)開采影響下植被群落指數(shù)下降了0～21.5%，但整個(gè)流域尺度的NDVI增加了15%；黨曉宏等[24]發(fā)現(xiàn)相比于非沉陷區(qū)，沉陷區(qū)土壤累積蒸發(fā)量、沙柳和油蒿累計(jì)蒸騰量分別下降28.8%，18.6%，16.1%。針對(duì)礦區(qū)植被遙感監(jiān)測(cè)的有關(guān)報(bào)道，主要集中在礦區(qū)植被演變以及植被生長特征的監(jiān)測(cè)，而關(guān)于生態(tài)脆弱區(qū)采煤沉陷區(qū)地表植被演化規(guī)律及其退化機(jī)制的研究仍較少。

鑒于此，筆者以棗泉煤礦為研究對(duì)象，基于SentineL-1A雷達(dá)數(shù)據(jù)和光學(xué)數(shù)據(jù)，利用SBAS-InSAR技術(shù)研究采礦擾動(dòng)下地表沉陷規(guī)律，分析工作面沉陷區(qū)與非沉陷區(qū)植被NDVI時(shí)序演化特征，評(píng)估采礦擾動(dòng)下礦區(qū)植被覆蓋度時(shí)序演化趨勢(shì)，并對(duì)沉陷區(qū)植被發(fā)育影響因素進(jìn)行剖析。研究結(jié)果對(duì)深入認(rèn)識(shí)西部礦區(qū)煤炭資源開采對(duì)生態(tài)環(huán)境的影響具有重要意義。

1 研究區(qū)概況

棗泉煤礦位于寧夏回族自治區(qū)靈武市的東南區(qū)域，地處毛烏素沙漠附近(106.50°E—106.58°E，37.86°N—38.03°N)，井田南北方向長13 km，東西方向長4 km，總面積52 km2。目前主要開發(fā)的煤層為2煤和6煤，煤層平均埋深約200 m，煤礦地表多為黃土，風(fēng)積沙等松散介質(zhì)。在井田部分區(qū)域，采動(dòng)影響下地表裂縫較發(fā)育，裂縫周邊植被退化現(xiàn)象較明顯[25-26]。礦區(qū)位置及礦井工作面規(guī)劃布置如圖1所示。

圖1 棗泉煤礦地理位置圖及地表狀況Fig.1 Geographical location map and surface condition of Zaoquan Coal Mine

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

2.1 研究數(shù)據(jù)選取與預(yù)處理

采用2016-01-20—2016-09-16，覆蓋研究區(qū)且連續(xù)的10景SentineL-1A影像數(shù)據(jù)?；诓罘指缮鏈y(cè)量短基線集時(shí)序分析技術(shù)(SBAS-InSAR)，結(jié)合精密軌道數(shù)據(jù)及外部DEM數(shù)據(jù)，對(duì)原始影像進(jìn)行預(yù)處理，生成沉陷影像圖。

根據(jù)沉陷監(jiān)測(cè)時(shí)間，選取覆蓋研究區(qū)的2014-01—2021-12(采前-采中-采后)的93景月度LANDSAT8 OLI系列數(shù)據(jù)(2016年4月、2018年12月、2020年6月數(shù)據(jù)覆蓋較差，未選取)，數(shù)據(jù)選取時(shí)控制影像云量小于20%，并對(duì)影像進(jìn)行輻射定標(biāo)、大氣校正和幾何校正等預(yù)處理。

結(jié)合礦井工作面采掘接替資料，統(tǒng)計(jì)監(jiān)測(cè)期間棗泉煤礦工作面具體信息見表1。

表1 監(jiān)測(cè)期間棗泉煤礦工作面統(tǒng)計(jì)Table 1 Statistical table of Zaoquan Coal Mine working face during monitoring period

2.2 礦區(qū)植被覆蓋度演化特征分析

基于歸一化植被指數(shù)(NDVI)結(jié)果，根據(jù)像元二分模型對(duì)礦區(qū)植被覆蓋度進(jìn)行估算，并通過線性回歸分析及相關(guān)系數(shù)對(duì)礦區(qū)地表植被發(fā)育時(shí)序演化趨勢(shì)進(jìn)行表征，主要流程如圖2所示。其中，θslope為一元線性回歸方程斜率；ti為年份序號(hào)；xi為第i年生長季植被蓋度均值；R為相關(guān)系數(shù)；n為總年份數(shù)。

圖2 植被覆蓋度演化趨勢(shì)分析流程Fig.2 Analysis process of FVC evolution trend

3 采礦擾動(dòng)區(qū)地表植被發(fā)育特征

3.1 工作面開采沉陷特征

監(jiān)測(cè)期間工作面地表3個(gè)典型沉降區(qū)域，如圖3所示。提取3個(gè)典型沉陷區(qū)下沉-10 mm邊界范圍，并在3個(gè)沉降邊界范圍內(nèi)，沿走向中心線布置觀測(cè)線(L1，L2，L3)，分別繪制130202，110204以及220602工作面不同開采階段下沉曲線，3個(gè)工作面參數(shù)以及沉降特征見表2，下沉曲線如圖4所示。

圖3 工作面沉降結(jié)果及測(cè)線布置Fig.3 Face settlement results and survey line arrangement

由圖4可知，隨著工作面推進(jìn)距離不斷增大，地表下沉量最大值逐漸增加，沉陷速度呈現(xiàn)出先增大后減小的趨勢(shì)。結(jié)合表2可知，監(jiān)測(cè)期間，130202工作面總開采面積較110204工作面增加114%，最大下沉陷量增加46.5%，沉降盆地移動(dòng)距離增加300 m，說明較大的工作面開采面積，在地表產(chǎn)生的沉陷影響也較大；130202工作面與220602工作面沉降盆地移動(dòng)距離大致相等，工作面130202深厚比(58)小于工作面220602(214)，深厚比增加260%情況下，工作面地表最大沉陷量減少62%，說明深厚比較大的工作面，開采在地表產(chǎn)生的影響較小，增大工作面深厚比可有效降低采動(dòng)對(duì)地表的影響。

圖4 3個(gè)典型工作面不同開采階段下沉曲線 Fig.4 Subsidence curves of three typical working faces at different mining stages

表2 3個(gè)典型沉降區(qū)工作面開采參數(shù)和下沉特征統(tǒng)計(jì)Table 2 Mining parameters and subsidence characteristics of working faces in three typical subsidence areas

3.2 工作面沉陷區(qū)植被NDVI時(shí)序變化

由于礦區(qū)3個(gè)典型沉降區(qū)的形成時(shí)間均相同，為分析沉陷區(qū)植被NDVI時(shí)序變化特征，選取井田范圍植被長勢(shì)均勻、不受人為干擾的區(qū)域?yàn)閷?duì)照區(qū)，繪制2014—2021年3個(gè)沉陷區(qū)觀測(cè)線和對(duì)照區(qū)植被NDVI時(shí)序特征及其差值，如圖5所示，其中紅色箭頭表示開采時(shí)間。

根據(jù)圖5可以看出，對(duì)照區(qū)及沉陷區(qū)的NDVI均呈現(xiàn)出周期性變化趨勢(shì)，在每年的7、8月左右NDVI達(dá)到最大值，而在12月和1月NDVI達(dá)到最小值，這與區(qū)域夏季植被生長旺盛、冬季植被覆蓋較低的自然規(guī)律相吻合。綜合分析可以看出，對(duì)照區(qū)與沉陷區(qū)的NDVI差值在開采前期相對(duì)變化較小，在開采影響下，兩者的差值持續(xù)增大，說明采礦活動(dòng)對(duì)沉陷區(qū)地表的植被生長有較大影響。此外，在2018年后沉陷區(qū)與對(duì)照區(qū)的NDVI有了較大的提升，究其原因可能是棗泉煤礦于2018年對(duì)煤礦塌陷區(qū)進(jìn)行覆土和植被重建[23]，主要是在130202，110204以及220602工作面沉陷區(qū)進(jìn)行適生植被種植等修復(fù)措施，促使礦區(qū)的NDVI有了較大的提高，說明植被修復(fù)重建可以改善礦區(qū)的植被覆蓋度。

圖5 對(duì)照區(qū)及沉陷區(qū)植被NDVI及其差值時(shí)序變化特征Fig.5 Temporal variation characteristics of NDVI difference between control area and subsidence areas

以130202工作面為例，開采前期，沉陷區(qū)的NDVI與對(duì)照區(qū)NDVI的差值較小，幾乎在0附近波動(dòng)，說明在采礦活動(dòng)前，工作面沉陷區(qū)與對(duì)照區(qū)的NDVI變化趨勢(shì)基本一致；在采中階段，沉陷區(qū)的NDVI開始低于對(duì)照區(qū)，NDVI的差值隨著時(shí)間增大，說明開采的擾動(dòng)已經(jīng)波及地表植被，對(duì)植被覆蓋度產(chǎn)生影響；在采后階段，NDVI差值持續(xù)變大，并達(dá)到最大值0.113 6，說明采礦活動(dòng)的影響明顯增大；而在2018年礦區(qū)進(jìn)行了大范圍植被修復(fù)，130202工作面沉陷區(qū)與對(duì)照區(qū)的NDVI均實(shí)現(xiàn)提高，兩者之間的差值逐步趨向于0，說明植被修復(fù)不僅可以改善礦區(qū)的植被覆蓋度，還可以補(bǔ)償采煤對(duì)塌陷區(qū)以及地表植被的擾動(dòng)效應(yīng)，降低采煤活動(dòng)對(duì)生態(tài)環(huán)境的影響。

通過分析發(fā)現(xiàn)，井下煤炭開采過程會(huì)破壞上覆巖層，進(jìn)而導(dǎo)致地表沉陷，破壞地表植被，降低沉陷工作面的NDVI，并且隨著采煤活動(dòng)的持續(xù)進(jìn)行，對(duì)地表的擾動(dòng)持續(xù)增加，地表植被一般在采后2～3 a受采動(dòng)影響最大，植被退化最為明顯，而隨著覆巖沉陷穩(wěn)定，地表植被逐漸開始恢復(fù)，開采擾動(dòng)區(qū)植被經(jīng)歷輕微退化期、顯著退化期以及逐漸恢復(fù)期等主要變化階段。結(jié)合工作面沉陷特征，發(fā)現(xiàn)不同沉陷特征區(qū)域植被受損程度不同，如130202工作面總開采面積大于110204工作面，工作面地表最大下沉陷也較大，地表植被NDVI下降幅度越明顯，也說明受擾動(dòng)程度就越大；220602工作面的深厚比大于130202工作面，采動(dòng)對(duì)地表的影響相對(duì)較小，植被NDVI受擾動(dòng)程度也較??；而在礦區(qū)植被修復(fù)方面，塌陷區(qū)的植被修復(fù)難度要大于非塌陷區(qū)，而且塌陷區(qū)的植被自維持能力也低于非塌陷區(qū)，由此可知，井下開采產(chǎn)生的地表沉陷會(huì)加速地表植被的退化。

3.3 礦區(qū)植被覆蓋度時(shí)序波動(dòng)程度

通過采用一元線性回歸法得到不同格點(diǎn)植被覆蓋度一元線性回歸斜率θslope，利用ArcGIS疊置分析功能對(duì)θslope和對(duì)應(yīng)相關(guān)系數(shù)R進(jìn)行疊加計(jì)算得到礦區(qū)植被覆蓋度空間變化趨勢(shì)[17]，并采用F分布檢驗(yàn)2個(gè)變量是否相關(guān)以及相關(guān)的顯著性水平。根據(jù)θslope和F檢驗(yàn)結(jié)果，基于生長季(4—8月)植被覆蓋度均值數(shù)據(jù)，將礦區(qū)植被覆蓋度變化情況分為：①θslope<0，P≤0.01，R∈[-0.87,-0.40),極顯著退化區(qū)；②θslope<0，0.010.05，R∈[-0.01，0)或θslope>0，P>0.05，R∈[0，0.01)，基本穩(wěn)定區(qū)；④θslope>0，0.010，P≤0.01，R∈[0.40,0.97]，極顯著改善區(qū)，如圖6所示。統(tǒng)計(jì)礦區(qū)植被覆蓋度變化不同程度面積及比例見表3。

由圖6和表3可知，2014—2021年研究區(qū)地表植被覆蓋度均值整體呈現(xiàn)小幅增加趨勢(shì)，植被覆蓋度呈升高與降低的區(qū)域面積分別占礦區(qū)總面積的53.91%，35.08%，其中極顯著退化區(qū)域植被受損最為嚴(yán)重，面積為7.33 km2，占礦區(qū)面積的12.9%，顯著退化區(qū)面積為12.59 km2，占礦區(qū)面積的22.18%。退化的區(qū)域主要分布在棗泉煤礦工作面上覆地表小區(qū)域、工業(yè)廣場(chǎng)周邊、矸石山等區(qū)域，值得注意的是，大約占整個(gè)礦區(qū)面積11.01%的區(qū)域地表植被覆蓋度變化呈現(xiàn)基本穩(wěn)定狀態(tài)，分布在礦區(qū)中南部。大約占礦區(qū)總面積的53.91%、面積為30.6 km2的區(qū)域地表植被呈現(xiàn)改善趨勢(shì)，分布在礦區(qū)北部和南部未受采礦影響的區(qū)域。分析發(fā)現(xiàn)，礦區(qū)植被覆蓋度的改善區(qū)面積顯著大于退化區(qū)面積，說明實(shí)施植被修復(fù)工程會(huì)提升礦區(qū)整體的植被指數(shù)，改善礦區(qū)生態(tài)環(huán)境。

表3 礦區(qū)植被覆蓋度變化不同程度面積及比例Table 3 Area and proportion of FVC changes in different degrees in mining area

圖6 礦區(qū)植被覆蓋度時(shí)序變化趨勢(shì)Fig.6 Time-series variation trend of vegetation coverage in mining area

4 礦區(qū)采損地表植被退化機(jī)制

4.1 植被NDVI與土壤理化性質(zhì)的關(guān)系

土壤理化性質(zhì)對(duì)地表植物的發(fā)育有著明顯的影響。豐菲等[27-28]發(fā)現(xiàn)土壤容重對(duì)灌木物種多樣性貢獻(xiàn)較大，土壤水分含量、有機(jī)質(zhì)與草本植物多樣性相關(guān)性較大。張玉秀等[29-30]認(rèn)為西部風(fēng)沙區(qū)開采引起的地表裂縫，使得土壤水分降低、養(yǎng)分流失，造成植被枯死，地表植被受影響。

在棗泉煤礦礦區(qū)范圍布置5個(gè)10 m×10 m的樣方，現(xiàn)場(chǎng)調(diào)研測(cè)試不同樣方區(qū)域的土壤理化性質(zhì)與植被生物數(shù)量，得到植被生長量與土壤理化性質(zhì)的關(guān)系見表4。從表4可知，植被的生長狀況與土壤有機(jī)質(zhì)含量和土壤含水率顯著相關(guān)，與土壤全氮含量和土壤全磷含量也呈正相關(guān)的關(guān)系。

表4 礦區(qū)植被生存狀況與土壤理化性質(zhì)相關(guān)性分析Table 4 Correlation analysis between vegetation survival and soil physicochemical properties in mining areas

在棗泉煤礦220602工作面沉陷穩(wěn)定的區(qū)域，選取位于工作面推進(jìn)位置的3條寬7，14，25 cm動(dòng)態(tài)裂縫(LF-1，LF-2，LF-3)及位于回采巷道兩側(cè)的3條寬25，38，131 cm的邊緣裂縫(LF-4，LF-5，LF-6)作為分析對(duì)象，裂縫區(qū)土壤樣品采集點(diǎn)布置，如圖7所示，分析不同埋深處6條裂縫的土壤含水率、有機(jī)質(zhì)、全氮、全磷含量，如圖8所示。

圖7 棗泉煤礦裂縫區(qū)土壤樣品采集點(diǎn)布置Fig.7 Soil sample collection site arrangement in fractured area of Zaoquan Coal Mine

由圖8(a)，(b)可知，6條裂縫周圍土壤平均含水率在垂直方向上整體呈現(xiàn)出隨著埋深的增加而增大的趨勢(shì)，并且邊緣裂縫對(duì)土壤含水率的影響范圍和影響程度都遠(yuǎn)大于相同寬度的動(dòng)態(tài)裂縫，此外，針對(duì)同一條裂縫而言，埋深越淺，受采礦活動(dòng)影響越大，含水率越低。土壤有機(jī)質(zhì)含量在埋深0～20 cm和20～40 cm，隨著裂縫寬度的增大，有機(jī)質(zhì)含量呈下降趨勢(shì)，且邊緣裂縫周圍土壤有機(jī)質(zhì)含量小于動(dòng)態(tài)裂縫，在埋深40～60 cm土壤有機(jī)質(zhì)含量較為接近，表明裂縫對(duì)土壤有機(jī)質(zhì)含量的影響范圍大致在埋深0～40 cm的淺層土壤。

由圖8(c)，(d)可知，土壤全氮含量與全磷含量整體隨裂縫寬度的增大而下降，6條裂縫在不同埋深時(shí)土壤的氮、磷含量均出現(xiàn)明顯下降，且埋深越大、含量越低。對(duì)于全氮含量分布來講，在垂直方向上，裂縫發(fā)育對(duì)表層土壤的影響程度高于埋藏較深土壤，且埋深0～20 cm的土壤平均全氮含量遠(yuǎn)大于20～60 cm；對(duì)于全磷含量分布，隨著裂縫寬度的增大，其對(duì)土壤全磷含量的影響程度也逐漸增大；隨著埋深的增大，裂縫對(duì)土壤全磷含量的影響逐漸減弱。

圖8 裂縫不同深度處土壤理化指標(biāo)Fig.8 Distribution of soil physicochemical indexes at different crack depths

綜上可知，采煤活動(dòng)產(chǎn)生的地裂縫會(huì)明顯降低土壤的含水率，改變土壤中有機(jī)質(zhì)、全氮、全磷的含量，且地裂縫的寬度越大，對(duì)土壤含水率、有機(jī)質(zhì)、全氮、全磷含量影響越大。而植被的生長狀況與土壤有機(jī)質(zhì)，土壤含水率顯著正相關(guān)。

因此可以推斷，采煤產(chǎn)生的地裂縫會(huì)影響土壤的理化性質(zhì)，進(jìn)而對(duì)礦區(qū)地表植被生長產(chǎn)生影響。

4.2 植被NDVI與地下水位埋深的關(guān)系

對(duì)干旱半干旱煤礦區(qū)而言，地下水是控制植被發(fā)育的決定性因素。馬雄德等[31-32]認(rèn)為隨著地下水位埋深的不斷增加，地下水對(duì)植被發(fā)育貢獻(xiàn)值減少，榆神礦區(qū)最適水位埋深上限為0.5 m，下限為4.0 m。王強(qiáng)民等[33]認(rèn)為只有當(dāng)干旱煤礦區(qū)采前地下水位埋深為1.0～2.5 m時(shí)，地下水位下降才會(huì)引發(fā)植被生態(tài)退化。金曉媚等[34]以海流兔河流域?yàn)檠芯繉?duì)象，認(rèn)為地下水位埋深在1.0～5.0 m時(shí)才對(duì)植被產(chǎn)生影響，當(dāng)?shù)叵滤怀^5 m時(shí)，植被變化主要取決于氣候因素。由此可見，地下水位埋深對(duì)植被的影響范圍基本處在0.5～5.0 m。

井下煤炭開采后，一旦導(dǎo)水裂隙帶直接導(dǎo)通含水層，將會(huì)引起地下水位發(fā)生明顯下降，范立民等[35]研究認(rèn)為，榆神礦區(qū)部分區(qū)域地下水位下降超過8 m，侯恩科等[36]發(fā)現(xiàn)開采影響下地下含水層水位出現(xiàn)“先下降后回升”并恢復(fù)至初始水位的動(dòng)態(tài)變化過程。

通過調(diào)研棗泉煤礦潛水水文觀測(cè)井、周邊民井水位數(shù)據(jù)，結(jié)合礦區(qū)植被覆蓋情況，得到了礦區(qū)植被蓋度與地下水位埋深關(guān)系，如圖9所示。從圖9可知，棗泉礦區(qū)的植被適生地下水埋深在3.5 m左右。

圖9 植被蓋度與地下水位埋深關(guān)系Fig.9 Relationship between vegetation cover and water table depth of burial

根據(jù)棗泉煤礦地表Z1潛水長觀孔(地表高程1 310.5 m)數(shù)據(jù)資料，得到110207工作面開采前后水位變化曲線，如圖10所示。

圖10 110207工作面Z1潛水長觀孔水位變化曲線Fig.10 Variation curve of water level of long viewing hole on 110207 working face

從圖10可以看出，采前地下水埋深在3 m左右，適宜植被的生長，當(dāng)工作面開采位置過鉆孔時(shí)出現(xiàn)約3 m的下降，水位埋深約6 m，隨后水位逐漸恢復(fù)，恢復(fù)后的地下水埋深波動(dòng)在4 m左右，影響地表植被的生長。因此可以推斷，采礦活動(dòng)會(huì)降低水位標(biāo)高，導(dǎo)致地下水位下降，水位埋深增加，從而引起地表植被蓋度的降低。

4.3 植被NDVI與氣候因子的關(guān)系

通常來講，溫度和降雨量是影響區(qū)域植被生長最為重要的氣候因素[16]。通過鄰近的梅花井煤礦氣象觀測(cè)站收集2014—2021年礦區(qū)的月平均降雨量和平均氣溫，并基于礦區(qū)220602工作面的NDVI數(shù)據(jù)，分析礦區(qū)降雨量和溫度與工作面植被NDVI的關(guān)系，如圖11所示。從圖11可以看出氣溫和降雨量都與220602工作面的NDVI有相似的周期性變化，其中降雨量和NDVI的相關(guān)性系數(shù)為0.543 06，氣溫和NDVI的相關(guān)性系數(shù)為0.618 57，氣候因子中的溫度和降雨量與工作面的NDVI相關(guān)性較高，說明氣候因子也是影響地表植被發(fā)育的又一重要因素。

圖11 220602工作面NDVI與月平均氣溫和降雨量關(guān)系Fig.11 Relation between NDVI of 220602 working face and monthly mean temperature and rainfall

5 結(jié) 論

(1)2014—2021年3個(gè)開采沉陷區(qū)地表植被NDVI時(shí)序上呈物候性周期變化趨勢(shì)，工作面開采當(dāng)年地表植被開始出現(xiàn)退化，在采后2～3 a退化程度最為明顯，植被NDVI呈現(xiàn)先降低再緩慢升高的趨勢(shì)，總體經(jīng)歷輕微退化期、顯著退化期以及逐漸恢復(fù)期的變化過程。

(2)礦區(qū)植被覆蓋度呈升高與降低的區(qū)域面積分別占礦區(qū)總面積的53.91%和35.08%，植被覆蓋度總體以極顯著改善為主。

(3)土壤有機(jī)質(zhì)、含水率是影響地表植被發(fā)育的最重要土壤理化因子。采礦擾動(dòng)區(qū)地表下沉、裂縫發(fā)育導(dǎo)致的地表土壤理化性質(zhì)的改變以及地下水位下降是造成工作面植被退化的重要原因，溫度和降雨量也是影響地表植被發(fā)育的重要因素。