劉艷芳　雷莎莎　方然

摘要：以大冶市為例，基于大冶市2005和2010年的土地利用現(xiàn)狀圖及土地變更登記數(shù)據(jù)，分析了大冶市5年間的耕地動態(tài)變化情況，構(gòu)建并定量分析了耕地侵占因子指標體系，并結(jié)合 GIS空間分析技術(shù)對耕地易被侵占區(qū)進行了風險等級劃分。結(jié)果表明，建設用地侵占是大冶市耕地減少的主要原因；大冶市有48.56%的耕地位于被侵占風險較高的區(qū)域，耕地保護工作形勢嚴峻；耕地被侵占的風險存在地區(qū)差異，陳貴鎮(zhèn)、東岳街道、金湖街道、羅家橋街道、茗山鄉(xiāng)等鄉(xiāng)（鎮(zhèn)、辦事處）的耕地所面臨的被侵占風險要高于其他鄉(xiāng)（鎮(zhèn)、辦事處）。

關鍵詞：耕地動態(tài)變化；侵占因子；易被侵占區(qū)；風險等級；大冶市

中圖分類號：F323.211 文獻標識碼：A 文章編號：0439-8114（2015）05-1053-06

DOI：10.14088/j.cnki.issn0439-8114.2015.05.007

Abstract： Taking Daye City as an example， combining with the current land use map and the land alteration registration data in the year of 2005 and 2010， the dynamic changes of cultivated land were analyzed. The index system of cultivated land occupation factors were established and quantified. GIS spatial analysis technology was used for the division of vulnerable cultivated land area. The results showed that the construction land occupation was the main factor for the decline of cultivated land quantity in Daye City . 48.56% of cultivated land located in the area had a high risk of being invaded. Cultivated land protection faced serious situation. There were regional differences in the risk of being invaded that cultivated land was facing. The chance of being invaded， cultivated land in Chengui town， Dongyue street， Luojiaqiao street and Mingshan town was higher than that of others.

Key words： dynamic change of cultivated land； occupation factors； vulnerable area； risk class； Daye City

耕地是農(nóng)業(yè)生產(chǎn)最重要的生產(chǎn)資料，是輕工業(yè)生產(chǎn)原料的主要來源，豐富的耕地資源是一國國民經(jīng)濟發(fā)展和糧食安全的重要保障。中國是人口大國，人均耕地數(shù)量低于世界平均水平，且耕地后備資源不足。隨著社會經(jīng)濟的發(fā)展和城市化進程的推進，非農(nóng)用地大量占用耕地，耕地資源緊缺狀況愈發(fā)凸顯。為此，國內(nèi)外專家和學者從不同角度，就耕地利用動態(tài)變化趨勢及驅(qū)動力展開了深入研究和探討，研究領域主要集中于耕地動態(tài)變化的人文驅(qū)動力的定性分析，一些學者認為人口因素在耕地利用動態(tài)變化中起到了突出作用[1-4]，還有一些學者認為人口變化、國家制度、經(jīng)濟發(fā)展、政治經(jīng)濟政策、價值取向等人文因素是影響耕地利用動態(tài)變化的直接和間接因素[5，6]；建立模型對耕地利用動態(tài)變化的驅(qū)動力進行定量分析，如CLUE模型，Adams模型，CTM模型等；不同尺度下研究耕地動態(tài)變化驅(qū)動因子，其中以全國尺度和省域尺度居多，如封志明等[7]通過研究新中國成立至2003年這55年里的中國耕地數(shù)量變化后認為國家政策對耕地數(shù)量動態(tài)變化影響較大，也有學者認為社會經(jīng)濟的發(fā)展、農(nóng)業(yè)科技的提高是省域耕地動態(tài)變化的人文驅(qū)動力[8-10]。這些研究多基于對宏觀尺度的耕地利用動態(tài)變化的驅(qū)動力研究[11-13]，而對微觀尺度下根據(jù)某區(qū)域的耕地變化狀況來預判該區(qū)域耕地保護狀況的應用研究較少。前者能為國家制定相關政策提供科學的指導方向和決策依據(jù)，但對各土地管理基層單位的耕地保護工作的指導作用不明顯。本研究通過對大冶市耕地利用動態(tài)變化進行分析，探究該市耕地侵占因子指標體系，并結(jié)合層次分析法的研究原理對侵占因子的權(quán)重進行測算，運用GIS空間分析功能進行耕地易被侵占區(qū)域的風險等級劃分，得出了大冶市耕地所面臨的被侵占風險等級圖，旨在為大冶市的耕地保護工作提供科學依據(jù)。

1 研究區(qū)概況

大冶市位于湖北省東南部，隸屬黃石市，地跨東經(jīng)114°31′—115°20′，北緯29°40′—30°15′。大冶市土地總面積156 629.21 hm2，行政區(qū)域范圍包括羅橋街道辦事處、金湖街道辦事處、金山街道辦事處、茗山鄉(xiāng)、保安鎮(zhèn)、金山店鎮(zhèn)、殷祖鎮(zhèn)、大箕鋪鎮(zhèn)、劉仁八鎮(zhèn)、金牛鎮(zhèn)、靈鄉(xiāng)鎮(zhèn)、還地橋鎮(zhèn)、陳貴鎮(zhèn)、汪仁鎮(zhèn)、東風農(nóng)場等11個鄉(xiāng)（鎮(zhèn)）、3個街道辦事處和1個農(nóng)場。大冶市屬典型的亞熱帶大陸性季風氣候，冬冷夏熱，四季分明，光照充足，雨量充沛，氣候條件適宜耕作，截至2010年，大冶市的耕地保有量為47 780 hm2。

2 數(shù)據(jù)來源和研究方法

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

本研究所引用的數(shù)據(jù)包括：①大冶市2005和2010年土地利用現(xiàn)狀圖；②大冶市2005-2010年土地利用類型變更登記資料；③大冶市2005年和2010年統(tǒng)計年鑒。

2.2 研究方法

2.2.1 耕地侵占因子測算 根據(jù)2005年和2010年大冶市土地利用現(xiàn)狀圖，通過GIS空間分析技術(shù)，得到大冶市2005～2010年這6年的土地利用結(jié)構(gòu)的空間變化和土地利用類型變更數(shù)據(jù)，在此基礎上分析大冶市耕地動態(tài)變化的原因，并進行耕地侵占因子識別，構(gòu)建侵占因子指標體系。結(jié)合專家打分法和層次分析法（AHP）的原理[14，15]，測算耕地侵占因子的權(quán)重。層次方法原理簡單，有數(shù)學依據(jù)，已有廣泛應用，其方法大體可分為四個步驟：第一步，建立問題的遞階層次結(jié)構(gòu)；第二步，構(gòu)造兩兩判斷矩陣；第三步，由判斷矩陣計算被比較評價因子的相對權(quán)值；第四步，計算各層次因子的組合權(quán)重；第五步，進行一致性檢驗。層次分析法的計算過程如下：

2.2.2 耕地易被侵占的區(qū)域風險等級劃分 根據(jù)層次分析法計算的各個侵占因子的權(quán)重，結(jié)合GIS空間分析技術(shù)進行緩沖區(qū)分析。具體方法是參考有關研究的成果和大冶市的具體情況，圍繞各個侵占耕地的地類確定合理的緩沖距離分類標準，依次生成多個緩沖距離圖層，用ArcGIS9.2軟件依次進行柵格運算，計算各緩沖區(qū)域的得分，再將多個圖層加權(quán)疊加匯總，確定耕地易侵占區(qū)的風險等級并劃分。公式為：

式中，Sij為第ij個侵占因子的易侵占分值（等級）；gij（m）表示第m個侵占因子在第ij個格網(wǎng)的易侵占分值（等級）；Wm為表示第m個評價因子的權(quán)重。

3 結(jié)果與分析

3.1 耕地動態(tài)變化分析

運用ArcGIS9.2軟件對大冶市2005年和2010年土地利用現(xiàn)狀圖的對比分析，得到了大冶市土地利用結(jié)構(gòu)變化情況（表3）及耕地轉(zhuǎn)變?yōu)槠渌仡惖南嚓P數(shù)據(jù)（圖1、圖2）。5年間大冶市的耕地保有量由48 093.52 hm2下降為47 780.00 hm2，耕地面積減少了313.52 hm2。其他主要用地類型的面積都有不同程度的增加，其中園地面積增加了772.84 hm2，林地增加了682.26 hm2，城鎮(zhèn)建設用地增加了207.56 hm2，工礦用地增加了178.67 hm2，農(nóng)村居民點面積增加了89.59 hm2，交通水利用地增加了658.20 hm2。

耕地因地類變更而減少的313.52 hm2中，有22.43%被城鎮(zhèn)建設用地侵占，18.80%被工礦用地侵占，16.04%被農(nóng)村居民點侵占，13.88%被道路用地侵占，11.92被園地侵占，7.73%被水利設施用地侵占，6.39%被林地侵占，還有2.81%轉(zhuǎn)變?yōu)榱似渌仡?。根?jù)實地走訪調(diào)查，原因如下：城鎮(zhèn)建設用地對耕地的侵占，是由于大冶市的發(fā)展導致城鎮(zhèn)建設用地的增加，而城鎮(zhèn)周邊的土地因交通便利、市場需求大的緣故大多種植糧油蔬菜等作物，耕地比例大，建設用地的擴張逐步蠶食了城鎮(zhèn)周邊的耕地；大冶市是著名的礦區(qū)，2005年大冶市的工礦用地規(guī)模是3 389.60 hm2，到2010年工礦用地面積是3 568.27 hm2，增長了178.67 hm2，而工礦企業(yè)多分布在郊區(qū)，周邊多為農(nóng)用地，工礦用地的擴張勢必侵占部分耕地；農(nóng)村居民點用地侵占耕地是因為在大冶市農(nóng)民在自己的承包地上修建樓房的現(xiàn)象十分普遍，加上大冶市近年來的新農(nóng)村建設工作積極推進，這兩大因素導致農(nóng)村居民點的擴張和耕地資源相應減少；道路用地對耕地的侵占是因為在2005～2010年間，大冶市為了提高對外交通通達性修建了大量道路，如大廣高速湖北南段、武漢—大冶—陽新一級公路、武漢—黃石城際鐵路大冶部分等等，而耕地大多沿道路分布，侵占耕地等現(xiàn)象就不可避免了；耕地變更為園地是因為為了發(fā)展農(nóng)村經(jīng)濟，許多鄉(xiāng)鎮(zhèn)按照“一村一產(chǎn)業(yè)”的思路，積極發(fā)展農(nóng)業(yè)產(chǎn)業(yè)化，許多農(nóng)民放棄種糧，通過轉(zhuǎn)包、承包土地的方式大面積種植果樹，到2010年，園地面積由2005年的2 248.89 hm2上升到了3 021.73 hm2，增加了772.84 hm2；水利設施對耕地的侵占是因為大冶市修建引水灌溉工程，溝渠、泵房等設施占用了少量耕地；耕地變更為林地是因為大冶市按照土地整治的要求將一部分坡度較大而不適宜耕作的耕地進行了退耕還林。

3.2 侵占因子指標選取與量化

根據(jù)大冶市2005～2010年土地利用動態(tài)類型的分析結(jié)果，確定城鎮(zhèn)建設用地侵占、工礦用地侵占、農(nóng)村居民點侵占、道路用地侵占、園地侵占、水利設施用地侵占、林地侵占這7個侵占因子作為耕地動態(tài)變化的影響因子。將這7個侵占因子分為兩類，構(gòu)建侵占因子指標體系（表4）：一類為生態(tài)侵占，包括林地侵占和園地侵占；另一類為建設侵占，包括城鎮(zhèn)建設用地侵占、工礦用地侵占、農(nóng)村居民點侵占、道路用地侵占和水利設施用地侵占。

結(jié)合專家打分法和層次分析法計算上述侵占因子的權(quán)重，利用層次分析軟件YAAHP 7.5運算，并進行一致性檢驗，結(jié)果如表5。

3.3 易被侵占區(qū)域風險等級劃分

根據(jù)侵占因子的權(quán)重，運用GIS空間分析技術(shù)進行緩沖區(qū)分析，并將各緩沖區(qū)圖層進行疊加，得出大冶市耕地易被侵占區(qū)域風險等級劃分結(jié)果（圖3）。其中各侵占因子的緩沖距離分類標準結(jié)合了大冶市的實際情況，參考了相關研究及《水土保持綜合治理技術(shù)規(guī)范》（GB/T16453.1-1996—GB/T16453.6-1996），得到耕地動態(tài)變化侵占因子緩沖距離分類標準（表6）。

3.3.1 各風險等級區(qū)域的比例 對大冶市不同風險等級的耕地易被侵占區(qū)進行面積統(tǒng)計，各風險等級區(qū)的面積比例如下：高風險區(qū)的耕地面積占耕地總面積的21.32%，較高風險區(qū)的耕地面積占耕地總面積的27.24%，中風險區(qū)耕地面積占耕地總面積的36.03%，低風險區(qū)的耕地面積占耕地總面積的15.40%。其中，高風險區(qū)與較高風險區(qū)主要分布在建制鎮(zhèn)周邊，共有48.56%的耕地位于被侵占風險較高的區(qū)域，大冶市的耕地保護工作形勢較為嚴峻。

3.3.2 各鄉(xiāng)鎮(zhèn)耕地易被侵占區(qū)域風險等級的分布

圖4中Ⅰ級表示高風險區(qū)，Ⅱ級表示較高風險區(qū)，Ⅲ級表示中風險區(qū)，IV級表示低風險區(qū)。從圖4的數(shù)據(jù)中可以看出，陳貴鎮(zhèn)、東岳路街道、金湖街道、羅家橋街道、茗山鄉(xiāng)的耕地所面臨的被侵占的風險較其他鄉(xiāng)鎮(zhèn)略高；金牛鎮(zhèn)、靈鄉(xiāng)鎮(zhèn)、劉仁八鎮(zhèn)、汪仁鎮(zhèn)、殷祖鎮(zhèn)的耕地所面臨的被侵占的風險較其他鄉(xiāng)鎮(zhèn)略低。

陳貴鎮(zhèn)等鄉(xiāng)鎮(zhèn)和街道轄區(qū)是人口分布較為集中的行政區(qū)，經(jīng)濟比較活躍，且都位于大冶市中心及周邊，大冶市區(qū)的擴張和周邊鄉(xiāng)鎮(zhèn)的發(fā)展逐漸蠶食這些地區(qū)周邊的耕地資源；金牛鎮(zhèn)等地耕地易被侵占風險低是因為這些鄉(xiāng)鎮(zhèn)分布于大冶市的行政轄區(qū)邊緣，人口分布較少且分散，經(jīng)濟活躍度不高，城鎮(zhèn)化速度較為緩慢，且多丘陵地形，道路建設也相對滯后，耕地所面臨的被侵占的風險低。綜上，大冶市的耕地動態(tài)監(jiān)測的重點區(qū)域應為陳貴鎮(zhèn)、東岳路街道、金湖街道、羅家橋街道、茗山鄉(xiāng)等地，需密切關注這些地區(qū)的耕地動態(tài)變化，進行土地用途管制，從源頭上遏制耕地資源流失。

4 結(jié)論與討論

①2005-2010年大冶市的耕地動態(tài)變化的大體趨勢是耕地資源數(shù)量下降，建設用地對耕地的侵占是主要因素，其中城鎮(zhèn)建設用地、工礦用地、農(nóng)村居民點、道路用地等用地類型的擴張對耕地的侵占較為嚴重。為了遏制耕地資源非農(nóng)化轉(zhuǎn)變，大冶市應嚴格按照土地利用總體規(guī)劃的控制指標進行土地用途管制，摒棄城市發(fā)展“攤大餅”的思路，集約節(jié)約利用土地，積極推進城中村改造工程、農(nóng)村居民點拆并項目和工礦用地復墾，提高土地利用效率，嚴控耕地“非農(nóng)化”現(xiàn)象。

②2005-2010年大冶市的耕地動態(tài)變化的大體趨勢是耕地資源數(shù)量下降，建設用地對耕地的侵占是主要因素，其中城鎮(zhèn)建設用地、工礦用地、農(nóng)村居民點、道路用地等用地類型的擴張對耕地的侵占較為嚴重。為了遏制耕地資源非農(nóng)化轉(zhuǎn)變，大冶市應嚴格按照土地利用總體規(guī)劃的控制指標進行土地用途管制，摒棄城市發(fā)展“攤大餅”的思路，集約節(jié)約利用土地，積極推進城中村改造工程、農(nóng)村居民點拆并項目和工礦用地復墾，提高土地利用效率，嚴控耕地“非農(nóng)化”現(xiàn)象。

③大冶市有48.56%的耕地位于被侵占風險較高的區(qū)域，耕地保護工作形勢嚴峻，且耕地被侵占風險存在地區(qū)差異。其中，陳貴鎮(zhèn)、東岳路街道、金湖街道、羅家橋街道、茗山鄉(xiāng)等鄉(xiāng)（鎮(zhèn)、辦事處）的耕地所面臨的被侵占的風險要高于其他鄉(xiāng)（鎮(zhèn)、辦事處）。土地管理部門在工作中應有針對性地對高風險區(qū)域進行耕地動態(tài)監(jiān)測，及時監(jiān)測耕地的變更狀況，開展非法侵占耕地清理工作，堅決懲治侵占耕地尤其是基本農(nóng)田的違規(guī)違法行為。

④微觀尺度下的耕地利用動態(tài)變化研究是一個非常具有實用性的應用研究，如何科學地從耕地變化空間數(shù)據(jù)中提取耕地動態(tài)變化的驅(qū)動因子并構(gòu)建指標體系；探究合適的模型對指標體系進行量化測算；結(jié)合驅(qū)動因子指標體系進行耕地預警研究，提高耕地保護工作的高效性和科學性，是今后深入探究的主要方向。

參考文獻：

[1] 王秀蘭，蘇 忠.基于GIS的內(nèi)蒙古耕地時空變化研究[J].干旱區(qū)地理，1999，22（2）：71-76.

[2] 王秀蘭.土地利用/覆蓋變化中的人口因素分析[J].資源科學，2000，22（3）：39-42.

[3] 王 波，唐志剛，濮勵杰，等.區(qū)域土地利用動態(tài)變化及人文驅(qū)動力初步研究——以無錫馬山區(qū)為例[J].土壤，2001，33（2）：86-91.

[4] 劉旭華，王勁鋒，劉紀遠，等.國家尺度耕地變化驅(qū)動力的定量分析方法[J].農(nóng)業(yè)工程學報，2005，21（4）：56-60.

[5] STERN P C， YOUNG O R， DRUCKMAN D. Global Environmental Change： Understanding the Human Dimensions[M].Washington， DC： National Academies Press，1992.

[7] 封志明，劉寶勤，楊艷昭.中國耕地資源數(shù)量變化的趨勢分析與數(shù)據(jù)重建：1949～2003[J].自然資源學報，2005，20（1）：35-42.

[8] 張 建，程久苗，王 棟.安徽省耕地資源變化及驅(qū)動力研究[J].國土資源科技管理，2009，26（3）：74-78.

[9] 王海力.朱麗東.浙江省耕地動態(tài)變化及驅(qū)動力分析[J].廣東農(nóng)業(yè)科學，2011（17）：229-232.

[10] 許風娟.濟南耕地變化分析及驅(qū)動因子研究[J].海南師范大學學報（自然科學版），2007，20（4）：369-374.

[11] 楊瑞珍，陳印軍，郭淑敏.中國耕地資源流失深層原因及對策[J].中國農(nóng)業(yè)資源與區(qū)劃，2005，26（6）：37-41.

[12] 李 平，李秀彬，劉學軍.我國現(xiàn)階段土地利用變化驅(qū)動力的宏觀分析[J].地理研究，2001，20（2）：129-138.

[13] 楊 君，周子英，劉文潔.湖南省耕地資源動態(tài)變化驅(qū)動力分析[J].農(nóng)機化研究，2011（6）：37-40.

[14] 鄧 雪，李家銘，曾浩健，等.層次分析法權(quán)重計算方法分析及其應用研究[J].數(shù)學的實踐與認識，2012，42（7）：93-99.

[15] 常建娥，蔣太立.層次分析法確定權(quán)重的研究[J].武漢理工大學學報（信息與管理工程版），2007，29（1）：153-155.