基于3S的平朔礦區土地利用及景觀格局演變研究

韓武波2，賈 薇，孫泰森

（1.山西農業大學資源環境學院，山西太谷030801；2.山西省國土資源廳，山西太原030024；3.山西大學環境與資源學院，山西太原030006）

煤炭資源作為社會經濟發展重要的基礎材料與能源工業，為中國經濟建設做出了巨大的貢獻，但也造成了嚴重的生態環境問題與社會問題。主要表現為景觀生態功能喪失或降低，景觀生態環境惡化，生態平衡受到破環。現今礦區生產活動引起的土地利用和景觀格局蛻變的現象已受到重視，利用3S技術對監測土地利用變化，及時掌握土地演變信息并做出相應預測、優化、調控體系顯得更為必要。為了使礦區土地利用景觀格局演變這一動態問題處于可持續、平穩狀態，依據中國露天礦土地復墾與生態重建方法技術，應充分、相對完整地研究礦區的土地利用和景觀格局演變規律，對開展礦區人工復合生態系統建設具有重要的意義。

1 研究區域概況及區域景觀類型

1.1 研究區域概況

平朔礦區地處黃土高原晉陜蒙接壤的黑三角地帶，山西省北部的朔州市境內,地理坐標為東徑112°10′58"—113°30′，北緯39°23′—39°37′，包括安太堡露天煤礦、安家嶺露天煤礦和東露天煤礦三個大型露天礦田。本文中所指的平朔礦區是包含平朔礦區所在的7個鄉鎮，總面積約為826 km2。

1.2 研究區域景觀類型

本研究主要參考國土資源部制定的《全國土地分類（試行）》，同時綜合考慮遙感圖像分辨率和各景觀類型在礦產開采中的功能特征，將研究區景觀生態劃分為10種類型：草地、林地、耕地、建設用地、水域、未利用地、礦坑采掘區、排土場（已復墾）、礦坑邊坡、礦坑剝離區。

2 研究方法

首先從ENVI軟件中的Change Detection功能對3個時相的景觀組分圖進行柵格疊加，得到不同年份間景觀類型的轉移矩陣，從而定量、定位分析景觀組分在各斑塊類型間地相互轉變方向、強度和位置。在此基礎上，采用空間格局指數法，利用3個時相的景觀組分圖和景觀格局指數分析景觀變化，從斑塊類型水平和景觀水平兩個各層面上進行研究。利用景觀格局分析軟件FRAGSTATS，得出以下具有代表性的指標：斑塊類型面積、斑塊密度、最大斑塊指數、邊緣密度、斑塊平均大小、平均分維數、散布與并列指數、景觀形狀指數、香農多樣性指數、香農均勻度指數、蔓延度指數。

3 研究結果分析

3.1 平朔礦區的土地利用動態分析

通過對2001年、2004年、2009年3個時相的TM遙感影像進行處理，利用ArcGIS的統計工具可獲得礦區的土地利用轉移矩陣（表1、表2）。這兩個時段的土地利用動態變化過程分析如下：

（1）2001—2004年間，整體趨勢是耕地、水域、礦坑剝離區持續減少，耕地減少的面積最大，而其余地類都在增加，林地、草地增加的面積較大。從各類轉化關系來看，耕地33.83%轉化為草地，12.87%轉化為林地，主要原因是礦區積極開展退耕還林還草。林地、礦坑邊坡、礦坑采掘區的動態變化最為顯著，礦坑邊坡和礦坑采掘區的增長趨勢是由于煤炭的不斷生產所引起的。剝離區雖然減少的面積不大，但是每年的動態變化中下降的最快，是由于礦坑剝離區有35.78%轉化為已復墾的排土場，主要原因是礦區采取及時跟進的土地復墾措施。

表1 2001—2004年土地利用轉換矩陣Tab.1 Land use transition matrix from 2001 to 2004

（2）2004—2009年間，耕地、林地、水域、未利用地面積持續減少，草地、礦坑采掘區、礦坑邊坡、礦坑剝離區、排土場的面積都在增加。主要原因是東露天礦的開采建設，企業擴張兼并其周圍部分村莊；耕地的轉出率為42.78%，主要轉化為草地，轉出率為33.1%，主要原因除了退耕還林還草等保護措施以外，主要是企業在開采東露天礦區時注重邊開采邊復墾的保護措施，將土地復墾和生態重建的意識貫穿于整個建設過程。通過以上分析可以看出，2001年到2009年平朔礦區各土地類型的總轉移變化情況為：建設用地一直保持較高的轉入量和極低轉出量，說明人類的生產生活活動較為頻繁，擴大礦區規模諸如東露天礦的開采建設以及城鎮化的步伐加快等。而耕地則相反，草地和林地的轉入量較高。這主要是因為退耕還林還草政策的推行。因此，對于東露天礦的開采建設，今后應該更加積極開展土地復墾措施，以增加耕地的有效面積，將糧食保護放在安全要位。

3.2 平朔礦區的土地利用景觀格局演變分析

平朔礦區景觀格局演變的主要特征景觀破碎度和異質性增加，以及隨著東露天礦的建設發展，煤炭資源開采量的增加，工礦區規模不斷擴大，部分村鎮轉化為工礦建設用地。

3.2.1 景觀斑塊類型指數變化分析 采用FRAGSTATS軟件和EXCEL計算各時相景觀水平景觀指數（圖1—7）。

（1）斑塊數量和面積變化分析。耕地是最大的斑塊類型，面積呈遞減趨勢，2001—2004年期間減少了30%，同時耕地的斑塊密度增加，最大斑塊指數減少，從而表明耕地的破碎化隨著礦產資源的開發變得加劇；以及城鎮和農村居民點鑲嵌在廣大的耕地區域中，城鎮和農村居民點的擴張，人口增長，這些都會直接或間接的破壞耕地大斑塊的面積。建設用地斑塊、最大斑塊指數、斑塊密度均有所增加，這主要是由于東露天礦的生產建設，煤炭資源產量增加，引起企業規模擴大和礦區工業廣場建設面積增加，與此同時，礦區的擴張對周邊地區的征地政策，促使當地失地居民加速對其他土地的使用。這反映了人為開采活動的強干擾性，以及礦區的擴張建設所引發的土地社會職能問題。

表2 2004—2009年土地利用轉換矩陣Tab.2 Land use transition matrix from 2004 to 2009

（2）斑塊類型空間構型特征分析。在邊緣密度變化方面，林地、草地、未利用地的邊緣密度2001—2004年增加到2009年卻減少，說明林地、草地、未利用地的布局先呈分散化趨勢后呈集聚化趨勢，而建設用地、剝離區則相反。其中，林地的邊緣密度增加的幅度最為顯著，表明林地的破碎化速度和程度都最大。耕地、水域的邊緣密度呈遞減趨勢，表明耕地的布局逐漸分散，而礦坑采掘區、礦坑邊坡、排土場的邊緣密度呈增長趨勢，說明礦區作業區集聚化趨勢更加明顯。從各時相的散布與并列指數來看，在2001—2009年，排土場的IJI值為最大，說明現今排土場與其他所有類型近乎等量相鄰，亦說明礦區作業區的復墾狀況良好。

（3）斑塊形狀特征變化分析。從2001年到2009年，林地、草地、耕地、未利用地的平均分維數均有所增加，水域、礦坑采掘區是從2001年到2009年先減少后增加，說明斑塊形狀變的越來越復雜，景觀受人類的生產和生活的影響，穩定性變差。水域的平均分維數最大，說明其形狀最不規則。礦坑邊坡、礦坑剝離區是從2001—2009年先增加后減少，建設用地、排土場是呈遞減趨勢，說明隨著礦區采取的土地復墾措施，其形狀變得較為規則，景觀的穩定性較好。

圖1 2001—2009年斑塊類型面積指數Fig.1 The indices of CA in different years from 2001 to 2009

圖2 2001—2009年斑塊密度指數Fig.2 The indices of PD in different years from 2001 to 2009

圖3 2001—2009年最大斑塊指數Fig.3 The indices of LPI in different years from 2001 to 2009

圖4 2001—2009年邊緣密度指數Fig.4 The indices of ED in different years from 2001 to 2009

圖5 2001—2009年斑塊平均大小指數Fig.5 The indices of AREA_MN in different years from 2001 to 2009

圖6 2001—2009年平均分維數指數Fig.6 The indices of FRAC_MN in different years from 2001 to 2009

圖7 2001—2009年散布與并列指數Fig.7 The indices of IJI in different years from 2001 to 2009

3.2.2 景觀水平景觀指數變化分析 采用FRAGSTATS軟件和數據庫軟件計算各時相景觀水平景觀指數（表3）。

表3 2001—2009年各時相景觀水平指數Tab.3 The indices of landscape types in different years from 2001 to 2009

從表3得出：在研究期，景觀形狀指數先增后減，表明斑塊在自然與人為的干擾下，形狀先變得復雜化而后變得相對規整；香農多樣性指數增大，說明斑塊類型增加或各斑塊類型在景觀格局中的空間分布更加均衡化；香農均勻度指數增大，表明各景觀斑塊類型在空間上分布的均勻程度有所增加；散布與并列指數增大，說明景觀斑塊各類型間相鄰的類型有增大的趨勢，接壤的景觀類型豐富度有所增加；蔓延度指數減少，表明研究區小斑塊個數較多或優勢斑塊連接度低，景觀的破碎化較大，呈破碎化和異質化趨勢。

產生這些變化的原因主要為：礦區在脆弱的黃土高原生態環境氣候下，采取土地復墾和生態重建措施，土地利用多樣化和均勻化發展，而且近年來的退耕還林還草措施政策以及平朔煤礦的工業旅游規劃的建設等。

4 結論

平朔礦區整個景觀向著復雜化、多樣性、破碎化、邊界不規則化的方向發展，其主要原因為人類開采活動的強干擾性，但研究得出在此不斷的干擾演變中，土地復墾措施為礦區的生態重建奠定了重要基礎，有效地指引礦區的可持續規劃并且防止由于礦區生產生活活動引起的土地利用和景觀的過度蛻變。

本研究將工礦區的情況置身于當地社會、經濟、生態環境等各個方面，以追求一種長遠的、穩定的、可持續的狀態。生態建設目標的多元性和動態性，決定了礦區土地復墾內容和標準應隨其不斷調整優化。因此，隨著煤礦生產活動的推進，必須要綜合考慮礦區的開采活動對土地和植被等生態環境的影響，也要注重礦山開采活動對當地人口遷徙、村鎮搬遷以及生產生活方式的影響變化，加強以人口遷徙布局為核心內容對動態調控理論、技術和整體方案的研究，實施有效的恢復治理措施和政策，促進礦區生產生活環境的和諧發展。

（

）：

［1］ 白中科，趙景逵，段永紅，等.工礦區土地復墾與生態重建［M］.中國農業科技出版社，2000：7-14.

［2］ Forman R T T Godron M.Landscape ecology［M］.NewYork: John Willey&Sons，1986.

［3］ 傅伯杰，陳利璜，馬克明，等.景觀生態學原理及應用［M］.北京：科學出版牡，2001：59-65.

［4］ 鄔建國.景觀生態掌——格局、過程、尺度與等級［M］.北京：高等教育出版社，2000：16-34.

［5］ 肖篤寧，李秀珍，高峻，等.景觀生態學［M］.北京：科學出版社，2003：1-10.

［6］ 卞正富，張燕平.徐州礦區土地利用格局演變分析［J］.地理學報，2006，61（4）：349-358.

［7］ 田維淵.基于RS與GIS的雅安市土地利用景觀格局變化分析研究［D］.成都：成都理工大學，2010：3-8.

［8］ 王月榮.洞庭湖退田還湖區錢糧湖垸景觀格局、土壤質量與土地承載力研究［D］.華中農業大學，2010：4-30.

［9］ 田維淵，曾世斌，楊武年，等.雅安市生態環境遙感動態監測及景觀格局變化分析［J］.現代測繪，2010，33（2）：23-26.

［10］王娟，崔保山，姚華榮.云南瀾滄江流域景觀格局時空動態研究［J］.水土保持學報，2007，21（4）：85-97.