魏玉強, 程倩雯, 單金霞, 黃秋昊

(1.南京大學 地理與海洋科學學院, 南京 210023; 2.江蘇省地理信息技術(shù)重點實驗室, 南京 210023)

?

快速城鎮(zhèn)化大都市邊緣地區(qū)耕地紅線劃定研究

魏玉強1,2, 程倩雯1,2, 單金霞1,2, 黃秋昊1,2

(1.南京大學 地理與海洋科學學院, 南京 210023; 2.江蘇省地理信息技術(shù)重點實驗室, 南京 210023)

劃定耕地紅線區(qū)域可有效減緩快速城鎮(zhèn)化大都市邊緣地區(qū)建設(shè)活動對耕地的占用，是耕地保護的重要手段。以位于快速城鎮(zhèn)化大都市邊緣地區(qū)的常州市新北區(qū)為研究對象，在分析現(xiàn)狀耕地特點的基礎(chǔ)上，根據(jù)相鄰地類之間的適宜性選取擬退耕的耕地和宜墾的農(nóng)村居民點用地。將不同土地利用類型賦予不同阻力系數(shù)，運用最小累計阻力模型，計算各地類轉(zhuǎn)化為耕地的累計最小阻力。用標準方差分類法對累計最小阻力值分類，找出耕地適宜性分區(qū)閾值，并劃定耕地空間紅線區(qū)域。結(jié)果表明：劃定的耕地紅線區(qū)域面積為299.84 km2，為研究區(qū)域總面積的66.25%。紅線范圍區(qū)域內(nèi)耕地占耕地總量的92.50%，劃定的耕地紅線區(qū)域能起到耕地保護的作用。通過最小累計阻力模型劃定出耕地紅線保護區(qū)域，對耕地紅線內(nèi)不同土地利用類型采用不同的保護標準，可以為快速城鎮(zhèn)化區(qū)域耕地的保護工作提供參考。

耕地紅線; 空間評價; 最小累積阻力模型; 常州市新北區(qū)

近年來，我國城市化水平不斷提高，土地利用/覆被變化愈發(fā)顯著，特別是在一些快速城市化地區(qū)，土地利用變化帶來的土地利用矛盾日益突出[1]。在土地利用/覆被變化過程中，大都市邊緣地區(qū)的耕地與建設(shè)用地的矛盾尤為突出。大都市邊緣區(qū)域作為受城市輻射影響顯著的過渡地帶，是耕地變化非常劇烈的地區(qū)，也是非農(nóng)建設(shè)用地擴張與耕地資源保護矛盾最為激烈的地區(qū)之一[2]。耕地被建設(shè)用地占用和空間形態(tài)破碎導致耕地規(guī)模和質(zhì)量的下降，這引起許多學者的關(guān)注和研究。

現(xiàn)有耕地研究主要包括耕地數(shù)量[3-4]、質(zhì)量[5-6]、集約利用[7-8]、驅(qū)動力[9]等多個方面?，F(xiàn)有研究多從縱向方面進行分析，很少考慮耕地與周圍地類的關(guān)系。如對耕地質(zhì)量進行研究時，分析評價單元為耕地地塊，通過構(gòu)建指標評價體系，運用地圖疊加法和邏輯規(guī)則組合法進行耕地評價。這類方法僅考慮了地塊單元的縱向疊加效應，而忽視了其水平擴展效應[10]。為更好保護耕地，本文嘗試對水平擴展的耕地適宜性進行分析，采用最小累計阻力模型進行耕地適宜性分區(qū)評價，并劃定耕地紅線區(qū)域。最小累計阻力模型依據(jù)土地利用類型構(gòu)建阻力面，從而計算“源”、“匯”之間的累計阻力值，阻力越小則越有利于耕地的集中連片管理。

因此，本文運用相鄰地類之間的相互適宜性關(guān)系分析，判斷耕地和農(nóng)村居民點與相鄰地塊單元的土地利用類型的適宜性，因地制宜開展土地整理復墾開發(fā)，有效補充耕地資源，增加耕地集中連片度。本文以位于長江三角洲發(fā)達地區(qū)的常州市新北區(qū)為例進行分析，引用最小阻力模型對耕地紅線劃定分區(qū)進行探討，以期為其他快速城鎮(zhèn)化大都市邊緣區(qū)的耕地保護提供參考和借鑒。

1　研究區(qū)概況

常州市新北區(qū)位于常州市北部，北瀕長江，與泰州市隔江相望，東與江陰市相鄰，西與揚中市、丹陽市接壤，南接常州市老城區(qū)，下轄3街道6鎮(zhèn)。新北區(qū)地勢平坦，區(qū)域內(nèi)水網(wǎng)密布，河流縱橫。四季分明，年均氣溫16.2℃。2012年新北區(qū)土地總面積為452.60 km2，耕地面積165.63 km2，占全區(qū)土地總面積的36.60%。新北區(qū)地處長江三角洲快速城市化地區(qū)，位于常州市大都市邊緣區(qū)域。經(jīng)濟快速發(fā)展對土地資源尤其是耕地資源需求量大，新北區(qū)耕地連片度差，耕地斑塊分散破碎，容易被建設(shè)用地占用，耕地保護形勢嚴峻。

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

耕地紅線劃定流程如下：首先分析現(xiàn)狀耕地特點，對耕地與周圍地類的適宜性評價選擇出適宜整理復墾的農(nóng)村居民點用地，對農(nóng)村居民點用地與周圍地類相互關(guān)系分析預測擬退耕的現(xiàn)有耕地地塊；其次，選擇“源”面并構(gòu)建阻力面生成阻力值；最后，運用最小累計阻力模型生成最小累計阻力值并根據(jù)耕地鄰近區(qū)域內(nèi)地塊轉(zhuǎn)化為耕地的難易程度，并劃定耕地保護紅線區(qū)域范圍。

2.1數(shù)據(jù)來源與處理

本文數(shù)據(jù)主要包括：2012年新北區(qū)1∶1萬土地利用變更數(shù)據(jù)，用于提取耕地等各種土地利用類型。最小累計阻力模型所需要數(shù)據(jù)均為柵格圖層，因此本文中的土地利用變更數(shù)據(jù)轉(zhuǎn)換成統(tǒng)一的10 m×10 m的柵格數(shù)據(jù)。新北區(qū)30 m空間分辨率DEM數(shù)據(jù)，來源于中國科學院計算機網(wǎng)絡(luò)信息中心國際科學數(shù)據(jù)鏡像網(wǎng)站(http: ∥www.gscloud.cn)，主要用于計算新北區(qū)地形的坡度和坡向以排除坡度>25°的耕地斑塊。經(jīng)分析，新北區(qū)無坡度>25°的耕地分布。

2.2相鄰地類適宜性關(guān)系

相鄰地類適宜性是指某用地斑塊與周圍其他用地斑塊的一致程度，如果某用地斑塊被相鄰其他用地斑塊包圍，則其相鄰用地斑塊的適宜性差。若某A類型用地斑塊周長為C，與其相鄰的B類型用地斑塊公共邊界為X，令X/C=K。則K值越大，此A類型用地斑塊被B類型用地斑塊包圍程度越大。本文用K值表示用地斑塊之間的適宜性程度。一般認為，若某塊小面積農(nóng)村居民點用地被耕地包圍或農(nóng)村居民點用地斑塊零碎分布，此類型的農(nóng)村居民點用地與周邊耕地適宜性差。這類農(nóng)村居民點用地復墾為耕地的難度較小，可劃為宜墾的農(nóng)村居民點用地。一般認為，若某塊耕地被建設(shè)用地所包圍，即使其條件再優(yōu)越，也不適宜劃入到耕地紅線區(qū)域內(nèi)，這主要是因為耕地相鄰土地利用類型在一定程度上限制了其農(nóng)耕用途的可持續(xù)性[11]，此類型的耕地可以劃為擬退耕耕地。新北區(qū)農(nóng)村居民點等建設(shè)用地散亂分布在耕地斑塊之間，通過土地復墾補充耕地可提高耕地集中區(qū)域的耕地連片度。新北區(qū)有少量現(xiàn)狀耕地布局在城鎮(zhèn)中心被建設(shè)用地包圍，形態(tài)破碎易于被建設(shè)活動占用，此類耕地不宜劃入耕地紅線區(qū)域內(nèi)。

2.3耕地紅線劃定方法

2.3.1最小累計阻力模型最小累計阻力模型是指從源地經(jīng)過不同阻力的景觀類型所累積耗費的最小費用或者克服阻力做的最小的功[12]。它最早由Knaapen等[12]于1992年提出，經(jīng)國內(nèi)俞孔堅等[13]的修改，用下式表示：

式中：MCR——最小累計阻力值；f——一個最小累計阻力與土地類型變化的正相關(guān)關(guān)系，是一個單調(diào)增長函數(shù)；min——被評價的斑塊單元對不同的“源”取累計阻力最小值；∑——單元i與源j之間穿越所有單元距離和阻力的累計；Dij——斑塊單元從源j到斑塊單元i的空間距離；Ri——斑塊單元i轉(zhuǎn)化為“源”的阻力系數(shù)[14]。

2.3.2源地的識別“源”是事物或者事件向外擴散的起點和基地，具有內(nèi)部同質(zhì)性和向四周擴張或吸引的能力[15]。在耕地適宜性評價過程中，源地是集中連片面積較大的耕地。本文選取2012年現(xiàn)狀耕地斑塊合并后面積超過10 hm2且空間上集聚的地塊作為源地。

2.3.3阻力面構(gòu)建本文將不同的土地利用類型賦予不同的阻力系數(shù)來構(gòu)建阻力面。土地的空間異質(zhì)性使得耕地在不同方向吸引其他地類轉(zhuǎn)化為耕地的過程中受到的阻力不同。一般認為園地容易轉(zhuǎn)化為耕地，而公路用地和建制鎮(zhèn)建設(shè)用地等很難轉(zhuǎn)化為耕地。在耕地保護方面，阻力系數(shù)就是其他土地利用類型轉(zhuǎn)化為耕地時的難易程度或者適宜程度，它是一個相對概念并不是絕對值。其他土地利用類型轉(zhuǎn)化為耕地的難度越大，此種土地利用類型越不適宜劃入耕地紅線，其阻力系數(shù)越大；反之，則阻力系數(shù)越小。

參考劉孝富[14]、陳燕飛[15]等的阻力系數(shù)設(shè)定方法，設(shè)置阻力系數(shù)。在模型計算過程中，將“源”所在的柵格賦阻力系數(shù)為1，表示其轉(zhuǎn)化為耕地的阻力最小。不同連片度的耕地，其阻力系數(shù)不同。將源地之外的耕地根據(jù)其規(guī)模采用自然斷點法分3級并局部微調(diào)后賦給不同的阻力系數(shù)(表1)。城市和建制鎮(zhèn)等用地轉(zhuǎn)化為耕地的阻力最大，本文將其阻力系數(shù)設(shè)置為最高值100。園地的耕作層基本沒有遭到破壞，易于轉(zhuǎn)化為耕地；通過地類適宜性評價可以選擇出宜墾的農(nóng)村居民點用地，這類土地可以通過復墾轉(zhuǎn)化為耕地；設(shè)施農(nóng)用地、農(nóng)村道路、小型坑塘水面等其他農(nóng)用地，轉(zhuǎn)化為耕地阻力較小；面積較大的坑塘水面，轉(zhuǎn)化為耕地的阻力相比小型坑塘水面大；通過分析新北區(qū)坑塘水面現(xiàn)狀，一般認為，當坑塘水面面積小于0.3 hm2時，坑塘水面不會影響耕地的集中連片度，轉(zhuǎn)化為耕地的阻力相對較小。將林地、牧草地、河流水面、湖泊水面視為生態(tài)用地類。國家推進生態(tài)文明建設(shè)生態(tài)用地被嚴格保護，轉(zhuǎn)化為耕地的阻力很大。將新北區(qū)土地利用類型重分類，分類后組成標準及所賦值阻力系數(shù)如表1所示。

表1不同土地利用類型阻力系數(shù)

重分類土地利用類型阻力系數(shù)具體地類組成標準源地1規(guī)模較大布局集中耕地,本文確定的源地耕地1級5特點為具有一定規(guī)模且布局集中的耕地2級10特點為規(guī)模相對二級減少布局分散3級15遠離耕地集中區(qū)或者面積較小的耕地斑塊園地20園地宜墾地類30通過計算確定的宜墾的農(nóng)村居民點用地其他農(nóng)用地類40設(shè)施農(nóng)用地、農(nóng)村道路、農(nóng)田水利用地、田坎、小型坑塘水面大型坑塘水面60面積超過0.3hm2的坑塘水面規(guī)模較大村莊60排除宜墾地外的農(nóng)村居民點用地生態(tài)用地類70林地、牧草地、其他土地擬退耕耕地90通過計算和分析確定的擬退耕耕地城鎮(zhèn)用地類100城鎮(zhèn)用地、采礦用地、交通水利用地等

3　結(jié)果與分析

3.1耕地現(xiàn)狀評價

2012年新北區(qū)耕地總面積為167.01 km2，主要分布在春江鎮(zhèn)、羅溪鎮(zhèn)、孟河鎮(zhèn)和西夏墅鎮(zhèn)，占新北區(qū)耕地總面積的86.32%；薛家鎮(zhèn)和新橋鎮(zhèn)也有少量耕地和外圍鎮(zhèn)耕地集中毗連。新北區(qū)仍有少量耕地零散分布在城市和建制鎮(zhèn)的核心區(qū)域。

新北區(qū)耕地形態(tài)較破碎。耕地共有1.04萬個斑塊，其中：面積小于1 hm2的斑塊數(shù)量占總斑塊數(shù)量的66.44%，其面積僅占耕地總面積的11.58%。耕地斑塊面積差異懸殊，最大單個斑塊面積為59.16 hm2，平均斑塊面積為1.60 hm2。新北區(qū)平均耕地斑塊形狀指數(shù)為1.48，斑塊形狀不規(guī)則。耕地斑塊破碎化區(qū)域主要布局在城鎮(zhèn)建設(shè)用地周邊和道路兩側(cè)。耕地斑塊形狀不規(guī)則的原因主要是農(nóng)村居民點用地、坑塘水面等塊狀要素和農(nóng)村道路、河流等條狀要素分布在成片耕地中，呈現(xiàn)不規(guī)則的形態(tài)導致耕地形狀指數(shù)高。

3.2耕地變動區(qū)域分析

3.2.1宜墾的農(nóng)村居民點用地分析結(jié)合新北區(qū)實際情況，2012年新北區(qū)農(nóng)村戶均人口為3.1人[16]，農(nóng)村人均居民點面積為154 m2。一般認為，小于4戶或者是超過4戶小于15戶但被耕地包圍的農(nóng)村居民點用地適宜復墾為耕地。

經(jīng)分析選取后，新北區(qū)宜墾的農(nóng)村居民點用地面積為466.28 hm2，占新北區(qū)耕地總面積的11.60%。宜墾的農(nóng)村居民點用地平均斑塊面積為0.11 hm2，復墾難度小，可操作性較高。結(jié)合新北區(qū)農(nóng)村居民點用地復墾增加耕地經(jīng)驗，預計可新增耕地373.03 hm2。宜墾的農(nóng)村居民點用地在空間布局上多分布在較大農(nóng)村居民點用地邊緣，經(jīng)過復墾后可以使農(nóng)村居民點用地集約度提高，同時耕地的整體連片度增加。利用核密度分析法分析，新北區(qū)宜墾的農(nóng)村居民點主要分布在春江鎮(zhèn)、羅溪鎮(zhèn)和西夏墅鎮(zhèn)的交界處(圖1)。

圖1新北區(qū)宜墾的農(nóng)村居民點用地核密度

3.2.2擬退耕地塊預測新北區(qū)部分現(xiàn)狀耕地布局在建設(shè)用地內(nèi)部，此類型耕地空間上分布比較分散。結(jié)合新北區(qū)耕地現(xiàn)狀，本文將被建設(shè)用地包圍并且面積小于1 hm2的耕地斑塊和單塊面積不足0.05 hm2的耕地斑塊作為擬退耕耕地。分布于中心城區(qū)的耕地難以發(fā)揮其耕地的功能，本文將河海街道、三井街道、龍虎塘街道等分布在建設(shè)用地核心區(qū)內(nèi)的耕地亦作為擬退耕耕地。新北區(qū)擬退耕耕地面積為256.34 hm2，占現(xiàn)狀耕地總面積的1.53%。除位于中心城區(qū)的擬退耕斑塊外，其他擬退耕耕地斑塊均較小且分散。擬退耕耕地布局在建成區(qū)和道路兩側(cè)。位于建成區(qū)內(nèi)的擬退耕地形狀多不規(guī)則，而布局在道路兩側(cè)的多成條狀，其形成原因是建設(shè)用地的擴展對耕地的切割所致。

3.3耕地紅線邊界確定

3.3.1最小累計阻力表面計算參照表1將新北區(qū)土地利用斑塊重新分類并賦予不同阻力系數(shù)值的矢量數(shù)據(jù)轉(zhuǎn)化為柵格數(shù)據(jù)即可以得到圖2A所示的阻力系數(shù)空間分布圖。研究區(qū)內(nèi)阻力低值占主導類型，低值區(qū)域集中連片；阻力高值區(qū)域主要分布在新北區(qū)中心城區(qū)和各鎮(zhèn)的鎮(zhèn)中心區(qū)域；部分高值區(qū)域成條帶狀嵌入在低值區(qū)域內(nèi)。將源地矢量圖(圖2B)和阻力面柵格圖輸入到ArcGIS下的cost-distance模塊，計算新北區(qū)每個柵格到源地的最小累計阻力值，得到如圖2C所示的最小累計阻力分布圖，最小累計阻力值基于標準方差法分成8類。

圖2新北區(qū)最小累計阻力值分析

3.3.2最小累計阻力閾值確定最小累計阻力模型在空間上是一種連續(xù)的概念，每一個像元或者圖斑都有其對應的阻力值，而耕地保護分區(qū)則是一種整體上分割而在局部上連續(xù)的一種分布格局[17]。將計算出的最小累計阻力利用標準方差法[18]進行分類，每類之間間隔為1/2方差，并建立最小累計阻力值方差分類與類別的柵格數(shù)目之間的對應關(guān)系(圖3)。

從圖3中可以看出，從C1到C2的過程中有一個突變的過程，柵格數(shù)量差別特別巨大，說明這里的最小累計阻力值產(chǎn)生了突變，C1類型柵格數(shù)目占到總柵格數(shù)目的66.25%，超過總數(shù)目的一半以上；C2和C3柵格數(shù)量差距也較大，只是相比C1與C2的差距?。粡腃3到C7，柵格數(shù)目變化不明顯；在C7到C8過程中，柵格數(shù)量略有增長。因此本文將C1對應的最小累計阻力值作為分割耕地紅線區(qū)域的閾值，其閾值為23 580。

圖3柵格數(shù)量與像元值對應關(guān)系

3.3.3耕地紅線區(qū)域保護根據(jù)確定的耕地紅線區(qū)域劃定閾值，可以將耕地紅線區(qū)域在空間上確定(圖4)。劃定耕地紅線是為了確定耕地在空間上的保護區(qū)域，劃定的耕地紅線區(qū)域面積為299.84 km2，為區(qū)域總面積的66.25%，主要分布在孟河鎮(zhèn)、西夏墅鎮(zhèn)、春江鎮(zhèn)、羅溪鎮(zhèn)等。耕地紅線區(qū)域主要是耕地的集中區(qū)域，耕地形態(tài)較好、連片度高、單塊耕地斑塊較大。紅線區(qū)域內(nèi)部現(xiàn)狀耕地面積為153.20 km2，占耕地總面積的92.50%。劃定的耕地紅線能起到對耕地的保護作用。耕地紅線區(qū)域范圍內(nèi)，包括耕地、農(nóng)村居民點用地、坑塘水面、建制鎮(zhèn)和園地等。其中區(qū)域內(nèi)耕地占耕地紅線區(qū)域總面積的51.10%，農(nóng)村居民點用地和坑塘水面分別占11.83%和11.20%。

注：C1為耕地紅線區(qū)域。

圖4新北區(qū)耕地紅線區(qū)域分布

耕地紅線區(qū)域內(nèi)部地類包括耕地和較容易轉(zhuǎn)化為耕地的園地、其他農(nóng)用地、農(nóng)村居民點和坑塘水面等地類，這些地類劃入到耕地紅線內(nèi)有利于區(qū)域景觀格局的穩(wěn)定。耕地紅線區(qū)域劃定后，紅線內(nèi)部的不同地類要根據(jù)實際采取不同管理措施。對耕地紅線區(qū)域內(nèi)部的現(xiàn)狀耕地，要執(zhí)行嚴格的耕地保護政策；對于區(qū)域內(nèi)部的農(nóng)村居民點用地，根據(jù)其所處位置和規(guī)模判斷保留還是復墾，并通過復墾提高耕地連片度和建設(shè)用地集約度；坑塘可以為耕地提供灌溉水源，因此紅線區(qū)域內(nèi)部的坑塘保持現(xiàn)狀有利于提高耕地產(chǎn)量。大面積連片坑塘水面可以作為生態(tài)用地美化環(huán)境，提高人們的生活質(zhì)量。

最小累計阻力值較高的區(qū)域分布在中心城區(qū)和春江鎮(zhèn)北部的長江生態(tài)區(qū)。中心城區(qū)距離耕地集中連片區(qū)域較遠，長江沿岸區(qū)域是重要的生態(tài)用地區(qū)并且遠離源地區(qū)域，不易變?yōu)楦貜亩纬勺枇Ω咧祬^(qū)域。位于高值區(qū)域內(nèi)部的零碎耕地，將其變更開發(fā)為生態(tài)用地或者是公園等城市綠地，有利于城市生態(tài)環(huán)境質(zhì)量的提高。

4　結(jié)論與討論

(1) 耕地紅線的劃定是為了更好地保護耕地資源。通過最小累計阻力模型可選出耕地集中連片區(qū)域來進行保護。相比于其他方法，此方法劃定的區(qū)域是一個連片范圍，而不是將耕地紅線細分到具體的地塊。在劃定過程中考慮了耕地水平擴展效應，分析其他土地利用類型轉(zhuǎn)化為耕地的適宜性程度，方法上可行適用。

(2) 通過對現(xiàn)狀耕地分析可知：新北區(qū)耕地主要分布在新北區(qū)春江鎮(zhèn)等鎮(zhèn)，耕地形態(tài)破碎、平均斑塊形狀指數(shù)為1.48。耕地斑塊面積差異大且細小斑塊多。耕地形態(tài)破碎原因主要是道路、農(nóng)村居民點等建設(shè)用地和河流對耕地的切割造成。新北區(qū)宜墾的農(nóng)村居民點用地面積為466.28 hm2，宜墾的農(nóng)村居民點空間集聚高密度區(qū)位于春江鎮(zhèn)、西夏墅鎮(zhèn)和羅溪鎮(zhèn)交界處。擬退耕耕地僅占現(xiàn)狀耕地面積的1.53%，主要布局在道路兩側(cè)和建設(shè)用地集中區(qū)域。

(3) 將現(xiàn)狀土地利用變更數(shù)據(jù)重分類，選擇出源地。根據(jù)不同土地利用類型轉(zhuǎn)為耕地的難易程度，設(shè)置不同的相對阻力系數(shù)。通過分類確定耕地紅線區(qū)域范圍，耕地紅線區(qū)域面積為區(qū)域總面積的66.15%，高于現(xiàn)狀耕地在研究區(qū)內(nèi)的比例。劃定耕地紅線區(qū)域后對紅線內(nèi)部的各現(xiàn)狀地類采用不同管理措施。耕地紅線區(qū)域劃定結(jié)果可為土地利用總體規(guī)劃時建設(shè)用地與耕地的合理布局提供參考依據(jù)。

(4) 確定耕地紅線區(qū)域需要綜合考慮規(guī)劃、耕地利用現(xiàn)狀、經(jīng)濟社會發(fā)展等因素。本文運用最小累計阻力模型考慮了水平方向其他地類轉(zhuǎn)化為耕地的適宜性，在阻力系數(shù)設(shè)置方面，人為主觀因素較大，進一步研究要結(jié)合專家打分綜合取得分。對于社會經(jīng)濟因素對耕地的布局影響和規(guī)劃對空間安排考慮不足，這有待后期研究中進一步深入探討。

[1]羅鼎,月卿,邵曉梅,等.土地利用空間優(yōu)化配置研究進展與展望[J].地理科學進展,2009,28(5):791-797.

[2]唐常春,陳烈,魏成.大都市邊緣區(qū)域耕地數(shù)量變化的時空特征及動力機制:以廣州市花都區(qū)為例[J].資源科學,2006,28(5):43-49.

[3]楊建云.基于EKC模型的河南省工業(yè)化、城鎮(zhèn)化水平與耕地面積關(guān)系研究[J].水土保持研究,2013,20(2):213-216.

[4]譚術(shù)魁,張紅霞.基于數(shù)量視角的耕地保護政策績效評價[J].中國人口·資源與環(huán)境,2010,20(4):153-158.

[5]韋仕川,熊昌盛,欒喬林,等.基于耕地質(zhì)量指數(shù)局部空間自相關(guān)的耕地保護分區(qū)[J].農(nóng)業(yè)工程學報,2014,30(18):249-256.

[6]郭洪峰,韓蕾,許月卿,等.云南省陸良縣耕地等級變化及其影響因素研究[J].資源科學,2014,36(10):2075-2083.

[7]龍冬冬,趙宏波,宋戈,等.河南省耕地集約利用時空變化特征研究[J].水土保持研究,2012,19(4):223-228.

[8]楊佳惠,雷國平,張健.煤炭資源枯竭型城市雞西市耕地集約利用評價[J].水土保持研究,2013,20(2):250-253.

[9]王瑞發(fā),夏非,張永戰(zhàn).青島市近10年來耕地變化及其驅(qū)動力分析[J].水土保持研究,2013,20(2):108-114.

[10]楊緒紅,金曉斌,郭貝貝,等.基于最小費用距離模型的高標準基本農(nóng)田建設(shè)區(qū)劃定方法[J].南京大學學報:自然科學,2014,50(2):202-210.

[11]董秀茹,尤明英,王秋兵.基于土地評價的基本農(nóng)田劃定方法[J].農(nóng)業(yè)工程學報,2011,27(4):336-339.

[12]Knaapen J P, Scheffer M, Harms B. Estimating habitat isolation in landscape planning [J]. Landscape and Urban Planning,1992,1(23):1-16.

[13]Yu K J. Land scape ecological security patterns in biological conservation [J]. Acta Ecologica Sinica,1999,1(19):8-15.

[14]劉孝富,舒儉民,張林波.最小累積阻力模型在城市土地生態(tài)適宜性評價中的應用:以廈門為例[J].生態(tài)學報,2010,30(2):421-428.

[15]陳燕飛,杜鵬飛.基于最小累積阻力模型的城市用地擴展分析[C]∥中國城市規(guī)劃協(xié)會和諧城市規(guī)劃:2007中國城市規(guī)劃年會論文集.哈爾濱:黑龍江科學技術(shù)出版社,2007.

[16]常州市統(tǒng)計局,國家統(tǒng)計局常州調(diào)查隊.常州統(tǒng)計年鑒[Z].北京:中國統(tǒng)計出版社,2013.

[17]李紀宏,劉雪華.基于最小費用距離模型的自然保護區(qū)功能分區(qū)[J].自然資源學報,2006,21(2):217-224.

[18]俞孔堅,黃剛,李迪華,等.景觀網(wǎng)絡(luò)的構(gòu)建與組織:石花洞風景名勝區(qū)景觀生態(tài)規(guī)劃探討[J].城市規(guī)劃學刊,2005(3):76-81.

Study on the Demarcating Method of Farmland Redline in Metropolitan Fringe Area with Rapid Urbanization

WEI Yuqiang1,2, CHENG Qianwen1,2, SHAN Jinxia1,2, HUANG Qiuhao1,2

(1.SchoolofGeographicandOceanographicSciences,NanjingUniversity,Nanjing210023,China; 2.JiangsuProvincialKeyLaboratoryofGeographicInformationScienceandTechnology,NanjingUniversity,Nanjing210023,China)

Demarcating farmland redline is an important means for farmland protection, which can slow down farmland occupation in metropolis fringe areas where massive construction activities frequently take place during urbanization proceess. Xinbei District in Changzhou City which is located in metropolitan fringe area with rapid urbanization was selected as study area in this research. Based on the analysis of the present farmland status and its characteristics, the returning farmland and reclaim-suitable rural residential land were then selected according to suitability of adjacent land patch. A specific resistance coefficient was then assigned to each land use type, and accumulated minimum resistance for other land use type convertedto farmland was calculated using MCR model. Accumulated minimum resistance was classified using standard deviation classification method in order to find a farmland property threshold, which demonstrates the location and shape of farmland redline. The result shows that the demarcated farmland redline region in Xinbei District has an area of 299.84 km2, which occupies 66.25% of the whole study area and contains 92.50% of current farmland, shall help to protect most of the farmland. This research shows that MCR model could be used to demarcate farmland redline effectively. Adopting different protection standards within different regions, which split by redline, can provide an important reference for farmland protection.

farmland redline; spatial evaluation; minimum cumulative resistance model; Xinbei District of Changzhou City

2015-03-05

2015-03-20

國土資源部公益性行業(yè)專項經(jīng)費“基于市域尺度的國土空間利用質(zhì)量評價提升技術(shù)研究與示范”(201411014-3)

魏玉強(1989—),男,山東日照人,碩士研究生,研究方向為土地利用規(guī)劃與GIS研究。E-mail:weiyq2013@163.com

黃秋昊(1977—),男,江蘇無錫人,副教授,主要從事土地利用規(guī)劃研究。E-mail:qiuhao.huang@gmail.com

F301.21

A

1005-3409(2016)01-0080-06