張 茹，何高波

(中煤航測(cè)遙感集團(tuán)有限公司，陜西 西安 710199)

農(nóng)村居民點(diǎn)是農(nóng)村人口居住與活動(dòng)的空間單元[1]，其空間分布狀況是一定自然、經(jīng)濟(jì)、社會(huì)發(fā)展環(huán)境下農(nóng)村居民居住活動(dòng)在其分布地區(qū)的反映[2-3]。研究農(nóng)村居民點(diǎn)的空間分布格局，對(duì)于農(nóng)村居民點(diǎn)整理規(guī)劃和統(tǒng)籌城鄉(xiāng)發(fā)展都具有重要的指導(dǎo)意義[4-5]。煤炭礦區(qū)農(nóng)村居民點(diǎn)由于受到煤炭資源的開采而遭受破壞，使原本缺乏規(guī)劃[6-8]的農(nóng)村居民點(diǎn)內(nèi)部廢棄地增加[9]，有的甚至全部廢棄。因此，研究煤礦礦區(qū)農(nóng)村居民點(diǎn)的分布規(guī)律有利于礦區(qū)農(nóng)村居民點(diǎn)的發(fā)展規(guī)劃和土地復(fù)墾，有利于促進(jìn)礦區(qū)的新農(nóng)村建設(shè)，也可為采煤沉陷區(qū)的移民搬遷提供有益的參考。

在礦區(qū)農(nóng)村居民點(diǎn)空間格局研究方面，我國的盧元清等以2009年土地利用現(xiàn)狀圖為基礎(chǔ)，對(duì)山西平朔礦區(qū)的居民點(diǎn)進(jìn)行研究，結(jié)果表明平朔礦區(qū)農(nóng)村居民點(diǎn)零散分布在中部和東南部地區(qū)，并且主要分布在露天礦區(qū)[10]；楊冰利用2013年遙感影像研究平朔礦區(qū)居民點(diǎn)，得出的結(jié)論是平朔礦區(qū)農(nóng)村居民點(diǎn)的的分布多與地形地貌有關(guān)，農(nóng)村居民點(diǎn)的優(yōu)選區(qū)位是高程1 300～1 500 m、坡度0～10°及坡向?yàn)镾、SW和SE坡[11]；王曉輝等通過解譯3期遙感影像對(duì)安徽淮南煤礦區(qū)土地利用變化進(jìn)行分析，得出的結(jié)論是淮南礦區(qū)農(nóng)村居民點(diǎn)面積呈先增加后減少趨勢(shì)，且土地綜合利用程度減少[12]。但是，這些研究的時(shí)間尺度都相對(duì)較短，不能充分說明農(nóng)村居民點(diǎn)的空間格局與時(shí)間變化之間的關(guān)系。

本文以軒崗礦區(qū)1986～2015年30年間的7期遙感影像數(shù)據(jù)為基礎(chǔ)，通過監(jiān)督分類和人工目視解譯相結(jié)合的方法，研究礦區(qū)農(nóng)村居民點(diǎn)的時(shí)空變化規(guī)律和特征，以期為采煤沉陷區(qū)移民搬遷和礦山環(huán)境恢復(fù)治理等提供理論支撐與實(shí)踐參考[13]。本研究中農(nóng)村居民點(diǎn)的概念不僅包括已有的行政村，還包括由于煤炭資源開采活動(dòng)而臨時(shí)增加或減少的居民聚集點(diǎn)。

1 研究區(qū)概況

軒崗礦區(qū)位于山西省寧武煤田的中部，北面是朔南礦區(qū)，南面是嵐縣礦區(qū)，地理位置為東經(jīng)112°9′～112°40′E，北緯38°44′～39°9′N，跨忻州市寧武縣、原平市及神池縣(圖1)。礦區(qū)內(nèi)賦存地層自上而下有第四系、新近系、侏羅系、三疊系、二疊系、石炭系及奧陶系。礦區(qū)位于寧靜向斜的北部，發(fā)育有大量的北東向正斷層，將北部切割成零星斷塊，地質(zhì)構(gòu)造較為復(fù)雜，較大的斷層有陽方口斷層。

圖1 軒崗礦區(qū)地理位置

礦區(qū)總面積約978 km2，煤炭資源總量約102億t。在煤炭開采過程中，長(zhǎng)期存在“多、小、散、亂”的開采格局，資源利用率較低，環(huán)境污染嚴(yán)重。2006年山西省政府發(fā)布《山西省煤炭資源整合和有償使用辦法》，煤炭資源整合在山西省全面實(shí)行。煤炭資源整合有利于優(yōu)化煤炭產(chǎn)業(yè)布局、保護(hù)生態(tài)環(huán)境和促進(jìn)行業(yè)安全生產(chǎn)[14]。但也存在一系列遺留問題，鄉(xiāng)鎮(zhèn)小煤礦是山西省煤炭資源整合的主體，它們的被整合對(duì)農(nóng)村居民點(diǎn)變化產(chǎn)生很大影響。

軒崗礦區(qū)劃分為3個(gè)井田，1個(gè)勘探區(qū)和3個(gè)資源整合區(qū)。礦區(qū)屬于中溫帶大陸性季風(fēng)氣候，年均溫為7.2℃，年均降水量為485 mm。礦區(qū)內(nèi)鐵路有北同蒲線及寧苛支線，公路主要有寧武至婁煩的寧白線(S215)、大同至忻州的大忻線(S206)、崞陽至五寨的崞五線(S305)三條干線，交通較為便利。礦區(qū)內(nèi)城鎮(zhèn)有寧武縣、陽方口鎮(zhèn)、余莊鄉(xiāng)、薛家洼鄉(xiāng)、段家堡鄉(xiāng)、軒崗鎮(zhèn)和長(zhǎng)梁溝鎮(zhèn)。

2 研究基礎(chǔ)和技術(shù)方法

2.1 基礎(chǔ)數(shù)據(jù)和資料

本研究收集采用了7期存檔的Landsat TM/OLI遙感影像，軌道號(hào)為126-33，成像時(shí)間分別為1986年6月8日、1990年8月22日、1994年5月29日、2000年7月1日、2006年6月15日、2010年7月12日和2015年6月8日，影像云量覆蓋小，成像質(zhì)量好，數(shù)據(jù)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(tái)(http://www.gscloud.cn)，圖幅號(hào)分別為L(zhǎng)T51260331986159BJC00、LT51260331990234BJC00、LT51260331994149BJC00、LT51260332000182BJC00、LT51260332006166BJC00、LT51260332010193IKR00和LC81260332015159LGN00。TM/OLI遙感影像的空間分辨率為30 m，地圖投影為UTM-WGS84，地球橢球?yàn)閃GS_1984。還收集到了一些其他輔助數(shù)據(jù)，包括礦區(qū)矢量圖、礦區(qū)煤炭開發(fā)利用現(xiàn)狀圖、2000年土地利用圖等。

2.2 技術(shù)路徑與方法

首先，對(duì)收集的遙感影像進(jìn)行預(yù)處理。因?yàn)門M/OLI遙感影像是L1T數(shù)據(jù)產(chǎn)品，已經(jīng)過系統(tǒng)輻射校正和地面控制點(diǎn)幾何校正，且已利用DEM進(jìn)行了地形校正[15]。所以，本研究對(duì)遙感影像的預(yù)處理只有大氣校正和研究區(qū)域提取。采用ENVI5.1軟件中的FLAASH大氣校正模塊進(jìn)行大氣校正，該模塊廣泛適用于多光譜和高光譜數(shù)據(jù)的大氣校正中。再用軒崗礦區(qū)矢量圖分別裁剪經(jīng)過大氣校正的7期遙感影像，提取研究區(qū)影像。

然后，利用ArcMap10.2軟件首先對(duì)2015年遙感影像進(jìn)行監(jiān)督分類，然后去除過大或過小斑塊以提高分類精度[16]，人工檢查并修改漏分錯(cuò)分斑塊，完成2015年土地利用分類圖，再從中提取2015年農(nóng)村居民點(diǎn)，生成礦區(qū)農(nóng)村居民點(diǎn)矢量圖[17-18]。以此為基礎(chǔ)，人工解譯其余影像農(nóng)村居民點(diǎn)。統(tǒng)計(jì)每個(gè)居民點(diǎn)在7期遙感影像中的面積，利用SigmaPlot軟件下Dynamic Fit Wizard命令中的線性函數(shù)進(jìn)行顯著性分析。若P<0.05，說明該居民點(diǎn)面積變化顯著，反之，不顯著。若擬合函數(shù)自變量系數(shù)a>0，說明居民點(diǎn)面積呈增加趨勢(shì)，反之，呈減少趨勢(shì)。

最后，利用2000年土地利用圖和實(shí)地調(diào)查的方法，對(duì)提取出的2000年農(nóng)村居民點(diǎn)進(jìn)行精度驗(yàn)證。精度驗(yàn)證包括點(diǎn)位精度驗(yàn)證和面積精度驗(yàn)證兩方面[19-20]。點(diǎn)位精度驗(yàn)證時(shí)，從2000年土地利用圖上均勻選取50個(gè)農(nóng)村居民點(diǎn)進(jìn)行檢驗(yàn)，其中有8個(gè)居民點(diǎn)沒有提取出來。進(jìn)一步從提取出的2000年農(nóng)村居民點(diǎn)中選取50個(gè)居民點(diǎn)進(jìn)行實(shí)地調(diào)查，經(jīng)調(diào)查發(fā)現(xiàn)有3個(gè)居民點(diǎn)為誤提居民點(diǎn)，可以得到點(diǎn)位精度為78%。面積精度驗(yàn)證時(shí)，對(duì)上述從土地利用圖上選取的50個(gè)居民點(diǎn)面積進(jìn)行統(tǒng)計(jì)，與提取出的農(nóng)村居民點(diǎn)相應(yīng)面積進(jìn)行對(duì)比，得到面積精度為83%。

3 數(shù)據(jù)分析與特征研究

3.1 農(nóng)村居民點(diǎn)數(shù)量變化情況分析

由圖2可以看出，1986～2000年軒崗礦區(qū)農(nóng)村居民點(diǎn)的數(shù)量呈緩慢增加趨勢(shì)。2000～2010年居民點(diǎn)數(shù)量沒有增減，保持不變。而2010～2015年居民點(diǎn)數(shù)量快速減少。

圖2 1986～2015年軒崗礦區(qū)農(nóng)村居民點(diǎn)數(shù)量變化

在整個(gè)研究期間都存在的農(nóng)村居民點(diǎn)共135個(gè)。不同研究時(shí)段農(nóng)村居民點(diǎn)數(shù)量的變化情況如表1、圖3所示。1986～1900年增加3個(gè)居民點(diǎn)，結(jié)合圖3可知，其中兩個(gè)居民點(diǎn)分別鄰近劉家梁煤礦和長(zhǎng)梁溝煤礦，另一個(gè)居民點(diǎn)在馬家灣村東邊工礦建筑物周圍。分析其原因，主要是由于煤礦建成初期，需要大量勞動(dòng)人員，煤礦為此修建礦工宿舍或礦工小區(qū)，在一定程度上促進(jìn)居民點(diǎn)數(shù)量的增加。1990～1994年沒有發(fā)生居民點(diǎn)增加或減少。1994～2000年增加3個(gè)居民點(diǎn)，分別鄰近寧白線、崞五線和大忻線。其原因主要是交通道路的改善帶動(dòng)沿線地區(qū)的發(fā)展，使一些居民點(diǎn)自發(fā)聚集在此。2000～2006年沒有發(fā)生居民點(diǎn)增加或減少。2006～2010年增加了1個(gè)居民點(diǎn)，位置鄰近同軒金海煤礦。此間，減少1個(gè)居民點(diǎn)，鄰近焦家寨煤礦。焦家寨煤礦開采較早，煤炭的大量開采使地表出現(xiàn)裂縫、塌陷等地質(zhì)災(zāi)害，對(duì)周邊居民點(diǎn)造成破壞。2010～2015年增加1個(gè)居民點(diǎn)，位于神達(dá)花溝煤礦西700 m處，減少10個(gè)居民點(diǎn)。減少的居民點(diǎn)主要為丁家梁村、前石湖村、染峪村和黃草坡村，分別鄰近程坑煤礦、老窯溝煤礦、榆樹坡煤礦和長(zhǎng)梁溝煤礦，西地村和西紅河村位于煤礦密集區(qū)。這些地區(qū)隨著煤炭資源的大量開采，采空區(qū)面積逐年擴(kuò)大，不利于礦區(qū)農(nóng)村居民點(diǎn)的發(fā)展。

圖3 1986～2015年不同時(shí)間段軒崗礦區(qū)

表1 1986～2015年不同時(shí)間段軒崗礦區(qū)農(nóng)村居民點(diǎn)數(shù)量變化

3.2 農(nóng)村居民點(diǎn)面積變化分析

如圖4所示，1986～2015年軒崗礦區(qū)農(nóng)村居民點(diǎn)總面積呈現(xiàn)先增加后減少的變化特征。1986～2000年居民點(diǎn)總面積呈增加態(tài)勢(shì)，年均增加0.122 km2。2000～2015年居民點(diǎn)總面積顯著減小，年均減少0.089 km2。

圖4 1986～2015年軒崗礦區(qū)農(nóng)村居民點(diǎn)總面積變化

對(duì)整個(gè)研究期間都存在的農(nóng)村居民點(diǎn)面積變化進(jìn)行分析，其中發(fā)生顯著變化的居民點(diǎn)如圖5所示。其中，面積顯著增加的居民點(diǎn)主要分布在礦區(qū)西部北辛窯井田和勘查區(qū)。馬圈村、瓦窯村和馮家墕村鄰近東南-西北走向的崞五線，苗莊村、暖水灣村和陽方村鄰近南-北走向的寧白線，這說明農(nóng)村地區(qū)交通基礎(chǔ)設(shè)施建設(shè)的不斷加快，對(duì)沿線農(nóng)村居民點(diǎn)的發(fā)展有巨大的推動(dòng)作用。莊子上村、半溝村和劉家園村緊臨寧武縣，袁家窯村位于陽方口鎮(zhèn)東200 m處，這表明城鎮(zhèn)的快速發(fā)展對(duì)周邊居民點(diǎn)有較好的帶動(dòng)作用。城鎮(zhèn)化的發(fā)展使大量勞動(dòng)力聚集到該地區(qū)，從而不斷加大對(duì)城鎮(zhèn)周圍住房的需求，使城鎮(zhèn)周圍農(nóng)村居民點(diǎn)用地規(guī)模也逐漸擴(kuò)大。

圖5 1986～2015年軒崗礦區(qū)顯著增加農(nóng)村居民點(diǎn)或顯著減少農(nóng)村居民點(diǎn)空間分布

面積顯著減小的居民點(diǎn)主要集中在礦區(qū)北部二號(hào)資源整合區(qū)和中南部劉家梁井田、焦家寨井田。礦區(qū)北部二號(hào)資源整合區(qū)是煤礦密集區(qū)，除煤炭開采沉陷對(duì)農(nóng)村居民點(diǎn)的不利影響外，煤炭資源整合也是一個(gè)影響因素。鄉(xiāng)鎮(zhèn)小煤礦的被整合一方面使原來煤礦為農(nóng)村公益事業(yè)和公共基礎(chǔ)設(shè)施建設(shè)提供的資金出現(xiàn)缺口，另一方面使當(dāng)?shù)匾慌迕袷ピ诿旱V工作的機(jī)會(huì)，加劇村民外出務(wù)工數(shù)量，不利于礦區(qū)農(nóng)村居民點(diǎn)的發(fā)展。

4 研究結(jié)論

本研究基于軒崗礦區(qū)1986～2015年30年間7期遙感影像數(shù)據(jù)，通過監(jiān)督分類和人工目視解譯相結(jié)合的方法提取農(nóng)村居民點(diǎn)，并對(duì)農(nóng)村居民點(diǎn)的數(shù)量和面積的時(shí)序和空間發(fā)展變化進(jìn)行顯著性分析，研究礦區(qū)農(nóng)村居民點(diǎn)的時(shí)序空間分布規(guī)律和影響變化的主要因素。通過研究分析得出以下幾點(diǎn)結(jié)論：

(1)1986～2000年軒崗礦區(qū)農(nóng)村居民點(diǎn)數(shù)量呈緩慢增加趨勢(shì)，2000～2010年居民點(diǎn)數(shù)量保持不變，2010～2015年居民點(diǎn)數(shù)量顯著減少。其中，1986～1900年增加3個(gè)居民點(diǎn)；1990～1994年沒有發(fā)生居民點(diǎn)增加或減少；1994-2000年增加3個(gè)居民點(diǎn)；2000～2006 年沒有發(fā)生居民點(diǎn)增加或減少；2006～2010年增加1個(gè)居民點(diǎn)，減少1個(gè)居民點(diǎn)；2010～2015年增加1個(gè)居民點(diǎn)，減少10個(gè)居民點(diǎn)。

(2)1986～2015年軒崗礦區(qū)農(nóng)村居民點(diǎn)總面積呈先增加后減少的趨勢(shì)。1986～2000年呈增加態(tài)勢(shì)，2000～2015年呈減小態(tài)勢(shì)。面積顯著增加12個(gè)居民點(diǎn)，主要分布在礦區(qū)西部北辛窯井田和勘查區(qū)；面積顯著減小25個(gè)居民點(diǎn)，主要集中在礦區(qū)北部二號(hào)資源整合區(qū)和中南部劉家梁井田、焦家寨井田。

(3)煤炭開采的不利影響造成的自然減少、區(qū)域的交通發(fā)展帶來的自然增加和政府的移民搬遷政策和行動(dòng)是影響礦區(qū)居民點(diǎn)時(shí)空格局變化的主要因素。城鎮(zhèn)和交通的快速發(fā)展對(duì)周邊居民點(diǎn)有較好的帶動(dòng)作用，而煤炭資源的大力開采在促進(jìn)經(jīng)濟(jì)發(fā)展的同時(shí)，對(duì)礦區(qū)居民點(diǎn)產(chǎn)生一些不利影響。煤炭采空區(qū)面積的逐年增大，造成大量居民住宅受損；煤炭資源整合加劇村民外出務(wù)工，減少農(nóng)村基礎(chǔ)建設(shè)資金來源，這些都不利于當(dāng)?shù)鼐用顸c(diǎn)的發(fā)展。

本研究尚存在一些不足之處。一方面是所采用的遙感影像空間分辨率為30 m，對(duì)一些較小農(nóng)村居民點(diǎn)的解譯存在困難，在測(cè)算居民點(diǎn)數(shù)量和面積時(shí)可能存在較大誤差；另一方面是，研究中對(duì)造成居民點(diǎn)增加或減少的因素可能存在片面性，因?yàn)檗r(nóng)村人口遷徙和居民點(diǎn)的變化影響因素很多。