本文摘要：摘要：采煤塌陷積水是東部高潛水位礦區(qū)最主要的生態(tài)環(huán)境問題，但塌陷湖的形成給生態(tài)環(huán)境提供了難得的水資源，因此，研究塌陷積水后對周邊生態(tài)環(huán)境的影響特征，對于指導礦區(qū)塌陷地生態(tài)恢復治理與濕地資源保護利用具有現實意義。本文以任樓煤礦塌陷積水區(qū)為研究區(qū)，以建
摘要：采煤塌陷積水是東部高潛水位礦區(qū)最主要的生態(tài)環(huán)境問題，但塌陷湖的形成給生態(tài)環(huán)境提供了難得的水資源，因此，研究塌陷積水后對周邊生態(tài)環(huán)境的影響特征，對于指導礦區(qū)塌陷地生態(tài)恢復治理與濕地資源保護利用具有現實意義。本文以任樓煤礦塌陷積水區(qū)為研究區(qū)，以建礦至今的5期遙感影像為數據源，通過多光譜計算和主成分分析獲取RSEI指數，并提取北、東、南、西四個方向數據，通過數據擬合與秩和檢驗。研究結果表明：①自2000年以來，研究區(qū)生態(tài)環(huán)境有較大改善，礦區(qū)土地復墾與生態(tài)恢復成效顯著;②塌陷積水對周邊生態(tài)環(huán)境總體上產生正面影響，影響具有顯著的方向差異性和年度變化特征;③塌陷積水對生態(tài)環(huán)境的影響隨距離呈“凸型”拋物線變化，平均影響范圍約3.7km，南北向約3.5km，東西向不超過7.8km，塌陷積水區(qū)形成后需要經土地整治工程才能最大程度發(fā)揮水域生態(tài)環(huán)境功能。
關鍵詞：高潛水位礦區(qū);遙感生態(tài)指數;塌陷積水;生態(tài)環(huán)境響應;任樓煤礦
由于煤炭開采造成的生態(tài)環(huán)境問題長期以來受到廣泛關注，礦區(qū)生態(tài)環(huán)境的監(jiān)測、評價、恢復治理等也是研究的熱點問題[1]。
我國東部煤田分布在高潛水位地區(qū)，井下開采造成地表塌陷和積水，尤其是近年來疊加深部開采，造成地面大范圍沉陷，逐漸形成塌陷湖泊濕地[2-3]。盡管地表沉陷積水損失了大量耕地，但是濕地具有調節(jié)氣候、涵養(yǎng)水源等生態(tài)功能，合理開發(fā)利用能產生顯著生態(tài)效益和經濟社會效益[4-6]，是東部礦區(qū)生態(tài)恢復的有利條件。皖北地區(qū)很多濕地資源沒有得到相應的保護，大多數煤礦企業(yè)對于濕地資源開發(fā)和保護重視程度不夠[7]，因此研究塌陷水域的生態(tài)環(huán)境效應是對礦區(qū)生態(tài)環(huán)境問題研究的進一步深化。衛(wèi)星遙感具有長時序、多信息的優(yōu)點，是地表過程研究的重要手段[8]。
近年來，遙感生態(tài)指數(RSEI)被應用于礦區(qū)生態(tài)環(huán)境變化研究中，成功地反映了礦區(qū)和局部尺度生態(tài)環(huán)境和優(yōu)劣等級的多年變化[9]。但是時序RSEI只能反映礦區(qū)生態(tài)環(huán)境的時空變化，無法進行生態(tài)環(huán)境變化歸因。本文嘗試將遙感提取數據與統(tǒng)計檢驗相結合，在研究煤礦塌陷積水區(qū)域一定范圍生態(tài)環(huán)境變化的基礎上，進一步探究塌陷積水區(qū)的形成對周邊生態(tài)環(huán)境的影響規(guī)律，可以為礦區(qū)生態(tài)環(huán)境恢復治理及濕地資源保護利用提供指導。
1研究區(qū)概況
研究區(qū)以皖北煤電集團任樓煤礦塌陷積水區(qū)為中心，東西長6.5km，南北長10km。任樓煤礦位于安徽省北部，隸屬于濉溪縣南坪鎮(zhèn)和蒙城縣許疃鎮(zhèn)，地理坐標為東經116°42′19″～116°48′13″，北緯33°25′36″～33°32′46″。井田東西寬4～7km，南北長8～11km，面積42.07km2，研究區(qū)土地覆被情況。
任樓煤礦自1997年12月31日建成投產以來，已經形成了塌陷積水區(qū)面積約340hm2，水深0～9.0m。研究區(qū)地處溫帶季風氣候區(qū)，夏季晴熱濕潤，冬季寒冷干燥，常年主導風向為ENE，次主導風向為E，夏季主要為ESE，溫度和降雨由東南向西北遞減。研究區(qū)內多年平均氣溫14.5℃，多年平均降水量847.2mm。研究區(qū)內地勢平坦，地表水體主要為采空塌陷積水，形狀不規(guī)則，積水主要來源為大氣降水和地表潛水，兩者為互補關系，以水面蒸發(fā)為主要排泄途徑。
2數據與方法
2.1數據來源
采用Landsat系列遙感影像數據，時點分別為2000年6月、2005年8月、2010年9月、2015年7月、2018年6月，其中2000年、2005年、2010年的遙感影像為Landsat5TM數據，2015年和2018年遙感影像為Landsat8OLI數據。空間分辨率為30m，云層覆蓋面積小于10%。
2.2遙感生態(tài)指數
RSEI模型見式(1)[10]。RSEI=(NDVI，WET，NDSI，LST)(1)式中：NDVI為綠度指數;WET為濕度指數;NDSI為干度指數;LST為熱度指數。其中，NDVI采用歸一化植被指數表示;WET采用纓帽變化中的濕度分量表示;NDSI由裸土指數和建筑指數構成;LST根據研究區(qū)實際情況采用地表溫度反演算法[11-14]。遙感影像應用ENVI5.3BandMath工具去除背景值干擾，經標準化處理和LayerStacking工具波段合成，通過主成分分析法和指標耦合方法，獲取各分指數的權重[15]。
遙感生態(tài)指數計算公式見式(2)。RSEI=0.504×NDVI+0.102×WET-0.307×NDSI-0.090×LST(2)
2.3統(tǒng)計檢驗
以塌陷積水區(qū)最大范圍邊界為起點，在北、東、南、西4個方向上按像元提取RSEI值，其中，北方向73個像元，距離2190m，南方向49個像元，距離1470m，東西方向各115個像元，距離各3450m。以2000年數據為基期，假定之后生態(tài)環(huán)境的變化僅由塌陷積水引起，采用Wilcoxon檢驗比較遙感生態(tài)指數在方向上、年度間的差異，采用二次函數擬合水域影響與距離的關系。
3結果與分析
3.1研究區(qū)全域RSEI時空變化分析
3.1.1RSEI及分指數的時間變化
經計算，研究區(qū)各期遙感生態(tài)指數、綠度指數、濕度指數、干度指數和熱度指數標準化后的年度變化�？傮w上，研究區(qū)4個分指數表現為2個 上升、2個下降;2個波動大、2個波動小;2個均值大，在0.5以上，2個均值小，在0.5以下。綠度指數呈上升趨勢，波動較小，2000年最低，均值為0.732，2018年最大，均值為0.844。
濕度指數呈總體上升趨勢，但波動較大，均值最低為0.46，也出現在2000年，最高均值為2015年的0.907，2018年略有下降。干度指數最低均值和最高均值分別為0.117和0.436，雖然2000—2015年、2015—2018年均表現為上升趨勢，但由于2010—2015年劇烈下降，總體上呈現下降趨勢，多項式擬合的二次項系數為-0.038。熱度指數呈總體下降趨勢，但波動不大，最小均值為2018年的0.288，最大均值出現在2010年，為0.358。
遙感生態(tài)指數與綠度指數的變化相似，2000年來一直呈增長趨勢，平均年度增長0.007，2010—2015年增長較快，年平均增加0.013，其他年度均在均值以下。由結果可知，2000—2018年，研究區(qū)生態(tài)環(huán)境質量總體上逐漸向好，積水區(qū)周邊綠度情況整體良好并逐年改善，濕度有較大幅度提升，干度和熱度總體下降。濕度與干度受年度降水量變化影響有較大幅度的波動，塌陷積水對改善區(qū)域濕度有較大貢獻。監(jiān)測資料顯示，研究區(qū)2000年以來，年平均氣溫逐漸上升，熱度變化平穩(wěn)下降應與塌陷積水有關。研究區(qū)除水域外，土地利用變化不大，建設用地、裸地 占比較低，耕地面積大，對于穩(wěn)定和提升區(qū)域生態(tài)環(huán)境均發(fā)揮了重要作用。
3.1.2RSEI與土地利用的空間變化
2000—2005年研究區(qū)處于塌陷積水初期，總體而言，生態(tài)質量變化不大，略微有上升。其北部建設用地面減少以及未利用地面積減少，生態(tài)環(huán)境質量好轉，中部塌陷積水區(qū)面積增加，生態(tài)環(huán)境質量顯著下降。
到2010年，煤礦塌陷積水區(qū)域進一步擴大，積水區(qū)周邊零星建設用地減少，耕地結構有所改善，東北部生態(tài)環(huán)境質量有所提高，塌陷積水區(qū)局部下降，周邊區(qū)域生態(tài)環(huán)境有所上升。到2015年，塌陷積水區(qū)進一步擴大，建設用地逐漸集中并略有增加，建設用地、耕地、未利用地減少，研究區(qū)生態(tài)環(huán)境仍不斷改善，積水區(qū)周圍耕地生態(tài)質量提升，中心積水區(qū)生態(tài)質量惡化。到2018年，塌陷積水區(qū)面積略有增加，建設用地、耕地面積基本無變化，積水區(qū)生態(tài)環(huán)境質量惡化，但西北區(qū)域以及南部區(qū)域耕地生態(tài)環(huán)境質量提升，研究區(qū)整體生態(tài)質量改善。
3.2研究區(qū)方向RSEI變化分析
3.2.1距離積水區(qū)遠近變化特征以塌陷積水區(qū)最大范圍邊界為起點，取遠離水域方向RSEI指數的北、東、南、西4個方向的均值，遙感生態(tài)指數呈波動變化，年度間的變化呈現逐年向好趨勢。一次擬合遙感生態(tài)指數隨距離增加而升高，截距為0.7668，斜率為1.955e-5，塌陷積水區(qū)外圍3200m范圍內遠離水域生態(tài)環(huán)境質量更好。在距離積水區(qū)距離分別為700m、1700m、3200m出現3組低值區(qū)，應與某些方向上該位置出現建設用地有關。
3.2.2方向上的變化特征
以南向距離長度為標準，各方向取塌陷水域外緣1470m(49像元)距離數據，將各年度RSEI值平均后數據特征統(tǒng)計，遙感生態(tài)質量在方向上差異大，南北比較相似。從中值和平均值、極值來看，其優(yōu)劣順序為南>西>北>東，從數據的分布來看，南、西、北向數據集中在高值區(qū)，東方向數據分布跨度大，高值區(qū)數據相對較少。遙感生態(tài)指數這種分布特征與研究區(qū)土地利用特征相符，該塌陷積水區(qū)南北方向為農田，東部為煤礦建成區(qū)，西部為村莊。
3.2.3時間上的變化特征
取各年度北、東、南、西4個方向數據的均值。數據的中值和平均值優(yōu)劣順序一致，均為2018年>2015年>2005年>2010年>2000年，整體上隨著時間推移逐漸提高。2010年前后是煤炭產能快速提升的階段，加劇了生態(tài)環(huán)境損害，2000年后礦區(qū)大力推進土地復墾工作，也使2005年遙感生態(tài)指數高于2010年。從數據分布看，2005年與2015年、2018年多位于高值區(qū)，2000年相對最差。遙感生態(tài)環(huán)境質量取決于煤炭開采強度和土地復墾與生態(tài)恢復力度之間的動態(tài)平衡。
3.3塌陷積水形成后對生態(tài)環(huán)境的影響
3.3.1距離對塌陷積水影響的差異分析
以2000年數據為基準，其他年度數據與之做差。塌陷積水對RSEI指數總體上產生正面影響，且隨距離由近及遠影響增強;4個方向上的影響距離和特點不同，具有明顯的差異性，即南北向影響隨距離變化大、東西向影響隨距離變化小。
由于南向距離較短，與北向在1500m范圍內同樣為增加的趨勢，因此選用北向影響距離作為南北向變化趨勢，二次擬合結果顯示，其影響距離3568m，在1784m即一半距離處的影響達到最大。東西向影響比較平穩(wěn)，二次擬合的影響距離為3480m，在1740m距離處影響最大，西向影響距離7844m，在3922m處影響最大。取四個方向的平均值，塌陷積水對周邊遙感生態(tài)指數影響范圍為3720m，在1860m處產生的影響最大。
由于西向1500～1700m處于2005年之后沿路建設了村莊，因此扣除2000年RSEI中沒有包含該變化，造成了數據在此處大幅度下降，其均值擬合曲線呈凹型，與其他方向上的趨勢相反。東向在500～1400m范圍經過煤礦廠區(qū)，也造成之后的生態(tài)指數較2000年下降，整體趨勢受到一定影響，造成影響范圍擬合結果偏大。而南向由于距離短數據少，上升趨勢受一定影響，但從數據分布看，其數據應與北向相似。由此推斷塌陷積水的影響范圍在南北向上約為3568m，在東西向上不超過7844m。
4結論與討論
4.1結論
本文采用遙感和GIS技術，獲取礦區(qū)塌陷地自采煤以來近20年的RSEI作為評估生態(tài)環(huán)境優(yōu)劣的指標，對研究區(qū)生態(tài)環(huán)境變化進行了總體時空變化分析，對提取的4個方向上RSEI數據進行了其距離、方向、時間上變化分析，經扣除基期干擾后分析了塌陷積水對周邊生態(tài)環(huán)境的影響特征。得出主要結論如下所述。
1)皖北煤電任樓礦域近20年來，生態(tài)環(huán)境逐漸改善，土地復墾與生態(tài)恢復工作成效顯著，生態(tài)環(huán)境的變化決定于煤炭開采活動強度和土地復墾與生態(tài)恢復力度的平衡。2)采煤塌陷積水總體上對周邊生態(tài)環(huán)境產生正面影響，影響強度隨距離增加呈“凸型”拋物線變化，平均影響范圍約3.7km，南北向和東西向影響范圍不同，南北向影響范圍約3.5km，東西向影響范圍不超過7.8km。3)采煤塌陷積水對周邊生態(tài)環(huán)境的影響具有顯著的方向性差異和年度變化特征。
4.2討論
水是生態(tài)系統(tǒng)不可或缺的重要環(huán)境要素，水分充足的地方往往容易形成良好的生態(tài)環(huán)境。我國東部高潛水位礦區(qū)，井下采煤形成大面積地表塌陷積水區(qū)，在損毀耕地的同時也為礦區(qū)生態(tài)環(huán)境建設提供了潛在的水資源條件。但是地表塌陷形成的積水湖與天然湖泊或人工湖有顯著不同，積水區(qū)周邊拉坡地土壤漬水嚴重，造成農作物減產或無法耕種。從影響強度隨距離的“凸型”曲線也可以看出，在靠近塌陷積水區(qū)生態(tài)環(huán)境受到的正面影響較小，甚至出現下降的情況，最主要的原因可能是拉坡地土壤漬水。因此，塌陷積水區(qū)形成后，必須通過土地復墾與生態(tài)恢復工程措施，抬高積水區(qū)周邊一定范圍的地面標高，才能更好地改善周邊生態(tài)環(huán)境，發(fā)揮地表水資源的巨大作用。
積水區(qū)東西兩邊由于村莊和煤礦工業(yè)區(qū)的存在和建設，造成水域對生態(tài)環(huán)境影響趨于平緩。這種變化可能是建設區(qū)位于塌陷盆地短軸，拉坡地不明顯，同時建設區(qū)保持了地勢較高、地形平整的原因。積水區(qū)西向對生態(tài)環(huán)境的影響大于東向，應與該地夏季主導風向為東風和東南風有很大關系，這種影響在南北向也有體現。這與林浩等[16]、王詠薇等[17]研究結果一致。因此可以推斷，塌陷積水對周邊生態(tài)環(huán)境的影響主要是通過改變地形提高土壤水分含量和主導風向輸送水汽產生影響;地形改變產生的影響變化大，風向的影響變化平穩(wěn)，塌陷積水區(qū)下風向更有利于生態(tài)環(huán)境改善。
參考文獻
[1]王雪芹，盛舞.我國礦區(qū)生態(tài)環(huán)境研究知識圖譜分析[J].河南理工大學學報(社會科學版)，2019，20(2)：47-53.WANGXueqin，SHENGWu.AknowledgemappinganalysisoftheresearchonecologicalenvironmentinminingareasinChina[J].JournalofHenanPolytechnicUniversity(SocialSciences)，2019，20(2)：47-53.
[2]周書東，王小霞，李廷芥.煤田開采誘發(fā)環(huán)境地質問題及防治對策[J].水土保持研究，2007，14(3)：351-354.ZHOUShudong，WANGXiaoxia，LITingjie.Geologicalenvi-ronmentalproblemsinducedbycoalfieldexploitionandcon-trolcountermeasures[J].ResearchofSoilandWaterConser-vation，2007，14(3)：351-354.
[3]卞正富，張國良.高潛水位礦區(qū)土地復墾模式及其決策方法[J].煤礦環(huán)境保護，1993，7(5)：2-5.BIANZhengfu，ZHANGGuoliang.Landreclamationmodelanddecision-makingmethodinhighgroundwaterlevelminingarea[J].CoalmineEnvironmentalProtection，1993，7(5)：2-5.
[4]XUJiaxing，YINPengcheng，HUWenmin，etal.Assessingtheecologicalregimeandspatialspillovereffectsofare-claimedminingsubsidedlake：acasestudyofthePan'anlakewetlandinXuzhou[J].PlosOne，2020，15(8)：e0238243.ht-tps：∥doi.org/10.1371/journal.pone.0238243.
作者：畢衛(wèi)華1，錢倬珺2，王輝2，姜念念2，黃文亭2，花逢春3
轉載請注明來自發(fā)表學術論文網：http:///jjlw/29117.html
