傅建春 李鋼

doi：10.15889/j.issn.1002-1302.2017.13.075[HT9.]

摘要：基于江蘇省沛縣各行政村2001、2007、2013年的居民點用地面積，利用位序-規(guī)模法則和空間集聚等分形理論，從不同角度對沛縣村級體系的規(guī)模和空間結(jié)構(gòu)進行研究。結(jié)果表明，沛縣的村級規(guī)模體系遵循位序-規(guī)模法則，聚集分形和形態(tài)分形特征顯著。在全縣尺度上，縣城南部農(nóng)業(yè)發(fā)展區(qū)農(nóng)村居民點趨于均勻分布，沛北礦區(qū)呈現(xiàn)集聚性特征，且空間形狀特征日益復(fù)雜，穩(wěn)定性降低。該結(jié)論為沛縣農(nóng)村居民點用地整理、布局優(yōu)化提供參考依據(jù)。

關(guān)鍵詞：農(nóng)村居民點用地；位序-規(guī)模法則；空間集聚；分形理論；江蘇沛縣；村級體系；空間結(jié)構(gòu)

中圖分類號： F291文獻標(biāo)志碼： A[HK]

文章編號：1002-1302（2017）13-0291-04[HS）][HT9.SS]

[HJ1.4mm]

收稿日期：2016-12-26

基金項目：國家自然科學(xué)基金（編號：41401610）；國土資源部公益性行業(yè)科研專項（編號：201211050）。

作者簡介：傅建春（1976—），男，江西新余人，博士研究生，講師，從事城鎮(zhèn)化、國土資源規(guī)劃與評價研究。E-mail：fjchpu@163.com 。

通信作者：李鋼，教授，博士生導(dǎo)師，從事城鎮(zhèn)化、礦區(qū)土地復(fù)墾、數(shù)字化測繪及地理信息工程研究。E-mail：cumtlig@263.net。[HJ]

[ZK）]

城鎮(zhèn)分布的復(fù)雜性可以通過分形進行研究，因受人類活動影響較大，分形理論在解釋城鎮(zhèn)體系規(guī)模-位序分布和空間特征等方面可以收到顯著的效果[1-2]，同樣受人類活動影響較大的農(nóng)村居民點的空間復(fù)雜性也可以通過分形來研究。李傳武以重慶市云陽縣南溪鎮(zhèn)24個村莊為研究對象，得出各村莊人口位序-規(guī)模服從Zipf定律；李玉華等研究了重慶不同生態(tài)類型區(qū)農(nóng)村居民點的分維數(shù)在一定時間跨度內(nèi)是逐漸上升的[3]；宋志軍等以北京城郊區(qū)總體、近郊平原區(qū)以及平原區(qū)為研究對象，發(fā)現(xiàn)城郊區(qū)不具有而近郊平原區(qū)和平原區(qū)所具有的明顯的分形特征[4]；車明亮等以山東省平邑山區(qū)農(nóng)村居民點為例，發(fā)現(xiàn)居民點的面積、周長、人口和分維具有一定的聯(lián)系[5]。農(nóng)村居民點空間特征在宏觀和微觀上分別體現(xiàn)在空間結(jié)構(gòu)和個體形態(tài)特征上[6]，研究空間特征的宏觀和微觀2種存在形式，有助于總結(jié)農(nóng)村居民點的用地特點。目前，利用分形理論研究農(nóng)村居民點的大多針對山區(qū)、丘陵和平原地區(qū)，而以礦區(qū)農(nóng)村居民點為研究對象的卻鮮有報道，因此本研究明確采煤塌陷區(qū)江蘇省沛縣農(nóng)村居民點空間分布所呈現(xiàn)的分形特征，以深化沛北一體化過程中農(nóng)村居民點的空間格局優(yōu)化進程。

1研究區(qū)概況與數(shù)據(jù)來源

1.1研究區(qū)概況

沛縣位于116°41′～117°09′E、34°28′～34°59′N，坐落在黃淮海區(qū)域江蘇省徐州市西北部，蘇、魯、豫、皖4省交界之地，是華東地區(qū)最大的煤炭工業(yè)基地的重要組成部分，總面積1 805.77 km2，其中徽山湖水域面積達320.21 km2，本研究范圍不包括徽山湖水面的陸地范圍（圖1）。全縣轄龍固、楊屯、大屯、安國等15個鎮(zhèn)，分布有8對礦井，多年開采形成的采煤塌陷地主要分布在沛縣北部的龍固、楊屯、大屯、安國、鹿樓、朱寨、沛城等鎮(zhèn)，累計塌陷2.26萬hm2，涉及村莊190多個，煤炭資源的開發(fā)和利用在促進沛縣經(jīng)濟發(fā)展的同時，造成土地的嚴(yán)重破壞，導(dǎo)致村莊的搬遷甚至多次搬遷，居民的生活、生產(chǎn)以及社會經(jīng)濟受到嚴(yán)重影響，進一步加深了人地之間的矛盾。同時，由于農(nóng)村宅基地使用管理制度及農(nóng)村自然發(fā)展的歷史原因，沛縣農(nóng)村居民點用地普遍存在“散、亂、大、空”現(xiàn)象，土地資源浪費嚴(yán)重。截至2013年年底，沛縣城鎮(zhèn)建設(shè)用地43 99.09 hm2，農(nóng)村居民點用地為18 549.49 hm2，占全縣城鄉(xiāng)居民點用地面積的80.83%，戶均占地面積 688.16 m2，人均達到181.11 m2，遠遠超過國家和江蘇省的用地標(biāo)準(zhǔn)，可見農(nóng)村居民點用地在全縣建設(shè)用地中所占比重極大，在強調(diào)加快城市化進程中城鎮(zhèn)建設(shè)用地擴大的同時，必須加強對農(nóng)村居民點用地的格局優(yōu)化和管理[7-8]。

1.2數(shù)據(jù)來源

從中國科學(xué)院地理空間數(shù)據(jù)云獲取了行列號為122/036的2001、2007、2013年空間分辨率為30 m的TM影像數(shù)據(jù)，衛(wèi)星過境時間分別為4月18日、4月3日、5月21日，由于此時間段的植被生長通常達到其最大程度上在該地區(qū)的生長季節(jié)，對識別居民點用地更為有利[9]。

對獲取的影像進行坐標(biāo)轉(zhuǎn)換、輻射校正和幾何校正處理后，利用沛縣的行政區(qū)矢量圖從中切割出沛縣的遙感影像范圍，借助eCognition 8.7軟件，根據(jù)各地物的波譜特征、形狀、紋理等特征，采用決策樹分類法，由易到難先后自動識別出水體、植被、道路和居民點用地類型[10]。由于本研究只須提取居民點用地，在用eCognition 8.7軟件自動提取居民點用地時，會誤提取部分裸地及旱地，后期利用ArcGIS結(jié)合野外調(diào)查，刪除誤判的圖斑，最終得到各年居民點用地圖斑（圖2）。

2研究方法

由法國數(shù)學(xué)家曼得爾布羅特提出的分形理論已經(jīng)廣泛應(yīng)用于各個領(lǐng)域，毫無疑問，分形理論也為土地利用空間結(jié)構(gòu)和空間形態(tài)的研究提供了新的理論支撐，本試驗對農(nóng)村居民點的分形研究基于以下幾個方面：用位序-規(guī)模法則研究居民點體系的規(guī)模分布特征，用集聚維數(shù)、形態(tài)維數(shù)分別反映農(nóng)村居民點的空間結(jié)構(gòu)、空間形態(tài)。

2.1位序-規(guī)模計算

村鎮(zhèn)居民點體系的規(guī)模與等級結(jié)構(gòu)能較好地表達規(guī)模的分布規(guī)律，通常可以用村鎮(zhèn)居民點用地規(guī)模來表征，而從村鎮(zhèn)規(guī)模和規(guī)模位序的相互關(guān)系表達村鎮(zhèn)體系的規(guī)模分布，就是村鎮(zhèn)位序-規(guī)模法則（別稱Zipf模型法則）。具體公式為endprint

[JZ（]lnP（k）=lnP-qlnk。[JZ）][JY]（1）

式中：k表示村莊的位序；P（k） 表示居民點用地規(guī)模排序后第k個居民點的用地面積；P表示常數(shù)；q體現(xiàn)出分維的性質(zhì)，常被稱作Zipf維數(shù)，和分位數(shù)D互為倒數(shù)，即D=1/q。當(dāng)系數(shù)q即分維數(shù)D為1時，規(guī)模分布呈等級規(guī)模分布，結(jié)構(gòu)達到最優(yōu)狀態(tài)；當(dāng)q﹥1、D﹤1時，居民點體系等級結(jié)構(gòu)較松散，相互之間等級規(guī)模差異較大；當(dāng)q﹤1、D﹥1，表示規(guī)模結(jié)構(gòu)較集中，處于中間位序的居民點數(shù)較多，其分布比較均勻[11]。

2.2集聚維數(shù)計算

空間集聚維數(shù)可以表達空間分布的自相似性，各居民點之間是按照某種自相似性圍繞某一中心成集聚狀態(tài)分布，且各個方向上的分形體變化均勻，由此有如下公式[12]：

[JZ（]Rs=（[SX（]1S[SX）]∑[DD（]ni=1[DD）]ri2）1/2；[JZ）][JY]（2）

遵循分維關(guān)系

[JZ（]Rs∝S1/D。[JZ）][JY]（3）

式中：Rs表示平均半徑；ri表示第i個居民點到中心居民點的歐式距離；S表示居民點個數(shù)；D表示分維，反映的是居民點圍繞中心居民點隨機集聚的特征，因此，也可將D表述為集聚維數(shù)。

公式（3）可以變換為如下公式：

[JZ（]lnRs=（1/D）lnS+C。[JZ）][JY]（4）

聚集維數(shù)D表示農(nóng)村居民點分布從中心向周圍的衰減特征。當(dāng)D值小于2 時，各農(nóng)村居民點圍繞中心居民點呈聚集狀態(tài)分布，且D值越小，其集聚程度越高，即農(nóng)村居民點分布由里向外逐漸衰減，此時在理論上居民點空間結(jié)構(gòu)處于優(yōu)化狀態(tài)；當(dāng)D值為2 時，各農(nóng)村居民點分布較均勻，空間結(jié)構(gòu)有進一步優(yōu)化的可能；當(dāng)D大于2 時，其分布由里向外遞增，一般來說這種情況不會出現(xiàn)。

2.3形態(tài)維數(shù)的計算

形態(tài)維數(shù)能夠通過定量分析圖斑邊界的曲折性和面積的大小關(guān)系表征圖斑的自相似性，具體公式如下：

[JZ（]lnA（r）=（2/D）lnP（r）+C。[JZ）][JY]（5）

式中：A（r） 表示某一居民點的面積；P（r） 表示同一居民點的周長；C表示截距；D表示分形維數(shù)，介于1～2之間。D的值與居民點形態(tài)的復(fù)雜性和穩(wěn)定性有關(guān)，D值越大，圖斑形狀越復(fù)雜，達到2時最復(fù)雜；D值越小，圖斑形狀越簡單，當(dāng)D值為1時，居民點圖斑為正方形；當(dāng)D值為1.5時，居民點圖斑處于最不穩(wěn)定狀態(tài)，也是其穩(wěn)定性的臨界點，即D值越接近15，居民點圖斑就越不穩(wěn)定，由此定義不穩(wěn)定指數(shù)S：

[JZ（]S=1.5-D。 [JZ）][JY]（6）

同理，S的值介于0～0.5，值越大表示居民點圖斑越穩(wěn)定，反之則越不穩(wěn)定。

3結(jié)果與分析

3.1農(nóng)村居民點位序-規(guī)模研究

根據(jù)位序-規(guī)模法則，利用遙感影像解譯的沛縣2001、2007、2013年的居民點圖斑，結(jié)合行政區(qū)劃界線，分別提取沛北礦區(qū)和沛南農(nóng)業(yè)發(fā)展區(qū)的居民點面積數(shù)據(jù)，用Excel軟件排序，以位序為橫軸、居民點面積規(guī)模為縱軸作雙對數(shù)坐標(biāo)線性擬合（圖3），各年份的擬合度良好，具體分維數(shù)q、分維度D、測定系數(shù)r2見表1。

由圖3可知，沛北礦區(qū)和沛南農(nóng)業(yè)發(fā)展區(qū)的居民點規(guī)模分布都服從Zipf定律，沛北首位居民點即城區(qū)的坐標(biāo)點遠位于回歸直線的上方，表明首位居民點的規(guī)模較大，在沛北甚至全縣范圍內(nèi)有較強的輻射帶動作用，只有2001年的雙對數(shù)線拖下一條長尾巴，說明存在很多規(guī)模較小的居民點，隨著農(nóng)村居民點的整治及新農(nóng)村建設(shè)的開展，部分規(guī)模太小的居民點逐漸被整合，這種現(xiàn)象可以從2007、2013年的雙對數(shù)曲線看出。

由表1可知，隨著時間的推移，全縣范圍、沛北及沛南的居民點用地規(guī)模擬合的相關(guān)系數(shù)r2均在0.7以上，證明擬合程度均較好；q值均小于1，表明沛縣居民點規(guī)模結(jié)構(gòu)符合位序-規(guī)模分布特征；分維數(shù)D均大于1且分別從2001年的1.181 3、1.141 0、1.396 8下降到2013年的1.173 0、1.054 5、1.300 9，其分維數(shù)特別是沛北礦區(qū)逐漸減小并靠近理想值1，說明各居民點之間規(guī)模較集中，中間位序的居民點較多。其中，沛北礦區(qū)居民點規(guī)模分布為壟斷性較強的首位型分布，而沛南農(nóng)業(yè)發(fā)展區(qū)首位型分布較弱。

3.2農(nóng)村居民點空間聚集的分形分析

借助ArcGIS軟件，將遙感解譯的居民點圖斑重心轉(zhuǎn)換成中心點圖層后，測算沛縣城區(qū)居民點與各居民點中心點之間的直線距離，可以求出其平均半徑Rs，并以lnS、lnRs分別為橫軸、縱軸作雙對數(shù)圖（圖4）?？梢钥闯鰺o標(biāo)度區(qū)間明顯，2001、2013年分別近似為[2.3-5.4]、[3.0-5.5]，測定系數(shù)r2分別為0.967 6、0.980 5，兩者線性相關(guān)顯著，表明沛縣農(nóng)村居民點符合分形特征，其空間結(jié)構(gòu)能用集聚維數(shù)描述。分形維數(shù)D分別為1.962 3、1.797 9，居民點體系集聚維數(shù)均小于2，表明沛縣農(nóng)村居民點已形成由縣城區(qū)向四周密度遞減的農(nóng)村居民點體系，同時驗證了縣城居民點的首位效應(yīng)。2013年的分形維數(shù)較2004年有下降趨勢，說明2001—2013年在沛縣城鎮(zhèn)化的帶動下，沛縣農(nóng)村居民點體系分形特征正在逐漸加強。但分維值還是較大，并接近于2，說明中心居民點周圍的集聚度還不夠，各個方向都均勻地分布著居民點，居民點體系分形結(jié)構(gòu)還不夠成熟。由解譯的居民點分布圖可以看出，有部分農(nóng)民將房屋建在自己的承包地上，特別是沿道路分布明顯，導(dǎo)致同一個村莊由多處住宅地組成，居民點破碎程度較高，這在一定程度上阻礙了農(nóng)業(yè)的規(guī)模經(jīng)濟，造成巨大的土地浪費。

3.3農(nóng)村居民點空間形態(tài)的分析

利用ArcGIS提取農(nóng)村居民點圖斑的周長和面積測算2001、2013年的形態(tài)維數(shù)（圖5），可知無標(biāo)度區(qū)間明顯存在，決定系數(shù)r2分別為0.946 2、0.957 5，表明沛縣的農(nóng)村居民點空間形態(tài)中周長與面積的關(guān)系存在分形特征。圖5可以計算出2001、2013年沛縣農(nóng)村居民點的形態(tài)維數(shù)D分別為 1191 2、1.232 2，分形維數(shù)值較小，說明屬于人工景觀農(nóng)村居民點受人類活動的干擾較強，圖斑自我相似性較強。同比之下，2013年較2001年的居民點空間形態(tài)更復(fù)雜且呈擴張之勢，空間形態(tài)越復(fù)雜說明空間穩(wěn)定性越小，表明農(nóng)村居民點擴建的隨意性越大。endprint

上述研究只能反映沛縣農(nóng)村居民點的總體形態(tài)結(jié)構(gòu)，無法識別沛北礦區(qū)與沛南農(nóng)業(yè)發(fā)展區(qū)的區(qū)域差異，所以進一步研究2個區(qū)域的居民點空間形態(tài)特征。由表2可知，沛北部礦區(qū)農(nóng)村居民點空間形態(tài)較沛南要復(fù)雜，表明該區(qū)域的農(nóng)村居民點圖斑邊界曲折程度大，形狀不規(guī)則，并且隨時間變化其復(fù)雜性逐漸增大，而沛南相對較穩(wěn)定，變化不大。

4結(jié)論與討論

沛縣農(nóng)村居民點空間結(jié)構(gòu)分形特征明顯。以農(nóng)村居民點用地面積作為衡量農(nóng)村規(guī)模的指標(biāo)，其規(guī)模分布符合位序-規(guī)模法則，沛北礦區(qū)和沛南農(nóng)業(yè)發(fā)展區(qū)擬合曲線的決定系數(shù)都在0.9以上，各農(nóng)村居民點規(guī)模較集中。從穩(wěn)定性來看，研究區(qū)農(nóng)村居民點用地形態(tài)相對穩(wěn)定，從穩(wěn)定性變化來看，沛北礦區(qū)由于受煤炭開采、經(jīng)濟和社會因素的影響，其農(nóng)村居民點具有較復(fù)雜的空間分布，其形態(tài)維數(shù)高于沛南農(nóng)業(yè)區(qū)，考慮生態(tài)安全和發(fā)展的限制條件，應(yīng)規(guī)劃對受塌陷及潛在塌陷的村莊進行遷建，遷建至鎮(zhèn)區(qū)或周圍新型農(nóng)村社區(qū)。以大屯、龍固鎮(zhèn)為代表的經(jīng)濟較發(fā)達的鎮(zhèn)和以楊屯、安國、朱寨鎮(zhèn)為代表的采煤塌陷區(qū)，村莊規(guī)模較大且聚集程度較高；以農(nóng)業(yè)為主導(dǎo)的沛縣南部各鎮(zhèn)村莊規(guī)模偏小，聚集程度較低，南北差異明顯。因此，在推進沛北一體化過程中，應(yīng)嚴(yán)把規(guī)劃關(guān)，推進“多規(guī)融合”、搭建礦地統(tǒng)籌平臺，兼顧沛南農(nóng)業(yè)發(fā)展區(qū)，綜合開發(fā)黃河故道，優(yōu)化沛縣布局，節(jié)約集約用地。

[HS2]參考文獻：

[1][ZK（#]車前進，曹有揮，馬曉冬，等. 基于分形理論的徐州城市空間結(jié)構(gòu)演變研究[J]. 長江流域資源與環(huán)境，2010（8）：859-866.

[2][JP2]傅建春，李鋼，趙華. 基于分形理論的徐州市城鎮(zhèn)體系規(guī)模結(jié)構(gòu)及空間分布研究[J]. 干旱區(qū)資源與環(huán)境，2015（10）：25-30.[JP]

[3]李玉華，高明，呂煊，等. 重慶市農(nóng)村居民點分形特征及影響因素分析[J]. 農(nóng)業(yè)工程學(xué)報，2014，30（12）：225-232.

[4]宋志軍，關(guān)小克，朱戰(zhàn)強. 北京農(nóng)村居民點的空間分形特征及復(fù)雜性[J]. 地理科學(xué)，2013（1）：52-60.

[5]車明亮，聶宜民，姜曙千，等. 平邑山區(qū)農(nóng)村居民點分形特征及影響因素[J]. 農(nóng)業(yè)工程學(xué)報，2010，26（增刊2）：360-365.

[6]李傳武，楊本俊，葉宋鈴. 基于分形理論的山區(qū)村鎮(zhèn)體系研究——以大別山區(qū)南溪鎮(zhèn)為例[J]. 皖西學(xué)院學(xué)報，2013，29（2）：12-15.

[7]李裕瑞，劉彥隨，龍花樓. 中國農(nóng)村人口與農(nóng)村居民點用地的時空變化[J]. 自然資源學(xué)報，2010（10）：1629-1638.

[8]楊存建，周成虎. TM影像的居民地信息提取方法研究[J]. 遙感學(xué)報，2000，4（2）：146-150，166.

[9]郭鵬，徐麗萍，王玲，等. 面向?qū)ο蟮默敿{斯河流域TM影像分類技術(shù)研究[J]. 水土保持研究，2013，20（3）：120-124.[ZK）]

[10][ZK（#]黃婷. 近15年來荊州市城鄉(xiāng)居民點用地規(guī)模時空變化特征研究[D]. 武漢：華中師范大學(xué)，2012.

[11]劉繼生，陳彥光. 城鎮(zhèn)體系空間結(jié)構(gòu)的分形維數(shù)及其測算方法[J]. 地理研究，1999，18（2）：60-67.

[12]莊至鳳，姜廣輝，何新，等. 基于分形理論的農(nóng)村居民點空間特征研究——以北京市平谷區(qū)為例[J]. 自然資源學(xué)報，2015，30（9）：1534-1546.endprint