基于地形起伏度的耕地與聚落耦合關系演變研究
——以都江堰市為例

王會豪， 任 平*, 張智波

(1. 四川師范大學 西南土地資源評價與監測教育部重點實驗室, 四川 成都 610066;2. 四川師范大學 地理與資源科學學院, 四川 成都 610066)

基于地形起伏度的耕地與聚落耦合關系演變研究
——以都江堰市為例

王會豪1,2， 任 平1,2*, 張智波1,2

(1. 四川師范大學 西南土地資源評價與監測教育部重點實驗室, 四川 成都 610066;2. 四川師范大學 地理與資源科學學院, 四川 成都 610066)

基于DEM數據提取都江堰市地形起伏度，利用都江堰市2005年與2011年耕地、聚落數據，借助ArcGIS軟件平臺，定量分析了都江堰市耕地、聚落及耕地-聚落耦合類型的時空分布特征及其分異規律.結果表明：1) 都江堰市的地形起伏度以平原為主、臺地和丘陵次之，山地面積較小；2) 耕地、聚落在平原地區分布面積最大，并且隨地形起伏度的增加不斷減少；3) 耕地偏少型主要分布在城區及建制鎮，耕地偏多型主要分布在臺地地區，耕地—聚落均衡類型主要分布在都江堰市城區周邊；4) 在不同地形起伏度上，以小起伏山地為界，耕地—聚落不同耦合類型在平原、臺地、丘陵地區變化較明顯，耕地偏少型都處于增加趨勢，而其它類型則都處于減少趨勢，在小起伏山地至極大起伏山地地區，各類型的變化幅度小且變化的趨勢比較復雜.

地形起伏度； 耕地； 聚落； 空間耦合； 都江堰市

耕地是人類最基本的生產資料，是人類賴以生存與發展的物質基礎，在我國人多地少的實際國情下耕地資源顯得尤為重要[1].對耕地的研究一直是地理學研究的熱點之一，特別是隨著我國經濟社會的快速發展、工業化的加速推進，耕地資源不可避免地被大量占用，對耕地資源的研究成果不斷增多，研究內容主要集中在耕地的時空分布特征[2-3]、景觀格局變化特征[4]、耕地變化軌跡及驅動力分析[5-8]、耕地非農化及保護對策[9-11]等方面.聚落是人類生產、生活的主要場所，按其人口規模與用地性質可以分為城市聚落與鄉村聚落[12].目前大部分研究主要集中在耕地或聚落本身，而探討耕地與聚落間耦合關系的研究則較為少見.耕地和聚落受人類活動的影響均比較深刻且二者間關系比較密切，特別是近年來，隨著我國城鎮化的快速推進，聚落面積迅速向外擴張，大量耕地被聚落用地侵占，耕地資源在城鎮化過程中呈現銳減趨勢，雖然我國實行了世界上最嚴格的耕地保護制度，但并沒有從根本上遏制住耕地不斷減少的態勢，耕地與聚落之間的矛盾日益尖銳.

都江堰市位于成都市西北部，地跨龍門山和成都平原岷江沖積扇扇頂部位，境內地形地貌復雜多樣，耕地、聚落空間分布差異明顯，對其進行研究具有重要的現實指導意義.本文以都江堰市為例，基于地形起伏度，運用GIS空間分析方法，探究了都江堰市耕地、聚落的時空演變特征，并揭示了耕地與聚落耦合關系的時空分異規律，試圖為都江堰市耕地與聚落的協調發展提供科學參考.

1 研究區概況

都江堰市位于北緯31°44′54″～31°02′9″，東經103°25′42″～103°47′0″之間，轄區總面積1 208 km2，下轄19個鄉鎮與一個街道辦事處.地勢西北高，東南低，山地、丘陵、臺地和平原呈階梯分布，受地貌類型影響，耕地與聚落主要集中分布在東南部地區的丘陵、臺地、平原區，山地地區分布較少.近年來都江堰市經濟社會發展迅速，2014年全市實現地區生產總值251.58億元，同比增長8%，城鎮化率超50%，但人均耕地面積僅0.042 hm2，低于聯合國糧農組織(FAO)劃定的0.053 3 hm2人均耕地警戒線，人地矛盾比較突出.

2 數據來源與研究方法

2.1 數據來源及處理 本文使用的數據包括都江堰市DEM數據(30 m×30 m)、都江堰市2005年、2011年土地利用現狀數據以及都江堰市鄉鎮級行政區劃矢量圖.首先利用ArcGIS 10.1參照國土資源部頒布的《第二次全國土地利用現狀調查技術規程》(TD/T1014—2007)，將都江堰市土地利用類型分為耕地、園地、林地、草地、聚落用地、采礦風景名勝及特殊用地、交通運輸用地、水域及水利設施用地、其它用地，并從中提取都江堰市耕地、聚落圖斑，最終得到本研究所需基礎數據.

2.2 研究方法

2.2.1 地形起伏度的提取 地形起伏度是描述一個區域地形特征的指標之一，是指在這個區域內海拔高差與地表切割程度的綜合特征[13]，其計算公式為

Hij=hij,max-hij,min,

i=1，2，3…n, j=1，2，3…n,

其中，Hij為該鄰域內的高差，hij為鄰域內每一個像元的高程值，hij,max為鄰域內最高像元的高程值，hij,min為鄰域內最低像元的高程值，由此得到網格單元n×n(2×2,3×3，…，15×15)與地形起伏度間的關系.

2.2.2 最佳統計單元的計算 關于最佳統計單元的確定，部分學者將像元網格面積與對應分析窗口中的平均起伏度數值生成曲線，采用目視判斷法[14]確定最佳統計單元，但其缺陷在于受主觀影響大;本文采用均值變點法[15-16]確定都江堰市的最佳統計單元，利用ArcGIS10.1中的焦點統計工具對DEM數據進行矩形窗口分析，以n×n像元作為模板算子，從2×2開始，首先計算像元內的最大值并記為專題層Max，接著計算最小值并記為專題層Min，然后利用柵格計算器計算最大值與最小值的差值，得到2×2像元地勢起伏度專題層，在計算過程中記錄網格面積的大小與平均起伏度數值，如此循環一直計算到15×15，經過計算最終確定都江堰市最統計窗口為4×4，由此得到都江堰市的地形起伏度.

3 研究結果

3.1 都江堰市地形起伏度規律分析 在ArcGIS10.1支持下，首先對4×4窗口下提取的都江堰市地形起伏度進行校準，然后依據我國1∶10×105地貌圖制圖規范[17]將都江堰市地形起伏度劃分為7個等級：平原(≤30m)、臺地(31～70m)、丘陵(71～200 m)、小起伏山地(201～500 m)、中起伏山地(501～1 000 m)、大起伏山地(1 000～2 500 m)、極大起伏山地(≥2 500 m)，如圖1和表1所示.

圖 1 都江堰市高程和地形起伏度分級

由圖1和表1可知,都江堰市地形起伏度空間分布差異明顯：極大、大、中、小起伏度山地主要分布在北部及西部地區，丘陵、臺地集中分布在中部地區，平原主要分布在東南部以及地勢較低的山麓與溝谷地區.整體而言，都江堰市地形起伏度主要以平原為主，臺地和丘陵次之.其中，平原所占比例最大為43.27%，臺地占30.85%，丘陵占22.04%，小、中、大起伏山地分別占1.89%、1.17%、0.76%，極大起伏山地所占比例最小為0.02%.說明都江堰市地形起伏度較大，地形地貌復雜多樣，山地、丘陵、平原等均有分布，其原因在于都江堰市北部及西部處于龍門山區，而東南部地區處于成都平原沖積扇.

3.2 都江堰市耕地時空分異特征 在ArcGIS技術支持下，將都江堰市地形起伏度圖層與2005年、2011年耕地斑塊圖層疊加(圖2)，并統計出各地形起伏度上耕地的分布情況，如表2所示.

圖 2 都江堰市2005年和2011年不同地形起伏度下耕地分布

從圖2和表2可知，以小起伏山地為界，隨著海拔不斷升高，都江堰市的耕地面積在不同地形起伏度上的變化大致呈先減少后增加的趨勢;超過99%的耕地分布在平原、臺地和丘陵3大地形起伏度單元，其中，全市近90%的耕地都分布在平原地區，分布在臺地、丘陵的耕地面積依次減少；小起伏山地及以上地形起伏度單元內盡管也有耕地分布，但面積很小，不足全市耕地面積的0.2%.2011年與2005年相比，耕地的變化也主要發生在平原、臺地和丘陵地區，其中，平原地區的耕地凈減少最多，達2 474.50 hm2；臺地和丘陵的耕地變化則表現為凈增加，分別為24.85 hm2和33.20 hm2，其耕地變化的絕對值遠遠小于平原地區；其它地形起伏度單元下的耕地則變化不大.

3.3 都江堰市聚落時空分異特征 在ArcGIS技術支持下，將都江堰市地形起伏度圖層與2005年、2011年聚落斑塊圖層疊加(圖3)，并統計出各地形起伏度上聚落的分布情況，如表3所示.

圖 3 都江堰市2005年和2011年不同地形起伏度下聚落分布

由圖3和表3可知，隨著海拔不斷升高，都江堰市聚落在各地形起伏度上的分布面積不斷減少，超過99%的聚落分布在平原、臺地與丘陵地區，其中全市近85%的聚落都分布在平原地區，10%左右的聚落分布在臺地地區，2%左右的聚落分布在丘陵地區.在小、中、大、極大起伏山地上雖有聚落分布，但面積很小，總和不足全市聚落總面積的0.2%.2011年與2005年相比，都江堰市的聚落面積在各地形起伏度上都有所增加，增加量主要發生在平原地區，達2 986.29 hm2;臺地、丘陵地區的增加量次之，分別為554.70、174.51 hm2;小、中、大起伏山地地區的增加量較小,分別為6.70、4.93、12.35 hm2;極大起伏度山地增加量最小,僅為0.12 hm2.

3.4 耕地聚落耦合類型時空分異特征 為定量表征都江堰市耕地與聚落的空間組合關系，本文運用耕地與聚落面積比指數(K)[18]分析二者間的耦合關系.計算方法為研究單元內耕地面積與聚落面積之比.該指數可以反映研究區域內耕地資源的豐缺程度與人地矛盾問題:比值越大說明研究區內的耕地資源比較豐富，人地矛盾比較緩和;反之,則說明耕地資源比較緊缺，人地矛盾比較尖銳.由于耕地與聚落空間分布的差異形成了不同耕地聚落耦合類型，為合理確定耕地—聚落耦合類型的劃定標準.本文將世界糧農組織劃定的人均耕地警戒線533 m2作為耕地緊缺閥值,現有世界人均耕地面積1 920 m2作為耕地盈余閥值,將《國家新型城鎮化規劃(2014—2020)》規定人均建設用地范圍[19](城市人均100 m2、農村人均150 m2，城市化率為60%，折算后人均約130 m2)作為緊缺閥值,現有世界人均建設用地面積83 m2作為盈余閥值[12]，根據耕地與聚落面積比指數計算得到K的盈余值為23.1、緊缺值為4.1.據此將都江堰市耕地—聚落耦合類型劃分為耕地偏多型、耕地偏少型，介于二者之間則為耕地—聚落均衡型.

由于都江堰市的行政區劃面積差異較大，以此研究都江堰市耕地聚落耦合類型的空間差異，無法準確揭示其分布規律[20]，本文依據都江堰市面積的大小，采用1 km×1 km的格網做為研究單元[21]，意義在于其空間粒度大大小于行政單元，所揭示的地理事物的詳盡程度比基于行政單元分析的結果詳盡的多[22].在ArcGIS技術支持下，得到覆蓋都江堰市的格網，最后將格網圖層與耕地、聚落圖層疊加，分別計算每個格網內的耕地面積與聚落面積的比值，并依據上述標準劃分為不同的類型，如圖4所示，最后將都江堰市耕地—聚落耦合類型圖層與都江堰市起伏度圖層疊加得到不同地形起伏度上耕地—聚落耦合類型分布狀況，如表4所示.

圖 4 都江堰市2005年和2011年耕地—聚落耦合類型分布

Fig. 4 Coupling type of settlements and cultivated lands distribution in 2005 and 2011 in Dujiangyan City

由圖4可知，都江堰市耕地—聚落耦合類型空間分布差異明顯，且隨著時間的推移發生著顯著的變化.耕地偏少類型區集中分布在都江堰市主城區及各建制鄉鎮且面積顯著增加，主要以中心城區向周邊擴張，尤以向東南部及南部地區擴張最為明顯；耕地—聚落均衡類型區主要分布在都江堰城區周邊且分布面積明顯減少；耕地偏多類型區主要分布在都江堰市主城區的東北、西北及西南部地區，即臺地地區較多，且數量有所減少；無耕地或聚落分布類型區主要分布在都江堰市北部地區，西部地區分布較少，整體處于減少趨勢.

表 4 都江堰市2005年和2011年不同地形起伏度下耕地—聚落耦合類型

由表4可知，都江堰市耕地—聚落耦合類型在不同地形起伏度上分布差異明顯，并隨著時間的推移也發生著較顯著的變化：

1) 2005年，從整體上來看，耕地偏少型、耕地—聚落均衡型、耕地偏多型以及無耕地或聚落分布的面積分別占全區總面積的30.37%、33.35%、5.17%和31.09%，耕地偏多型明顯處于劣勢地位，說明都江堰市的人地矛盾較為突出，其他3種類型的面積比例基本相當.從不同地形起伏度下的耦合類型來看，平原地區耕地—聚落均衡型占主導地位，面積比例占到了全區總面積的近1/4；其次是平原地區的耕地偏少型和丘陵地區無耕地或聚落的分布情況，其面積占比均在15%左右；臺地也存在相當數量的耕地偏少型與無耕地或聚落分布的情況，其面積占比也大致相當，分別為10.15%和11.02%.另外，隨著海拔的不斷升高，耕地偏少型與耕地—聚落均衡型在平原至丘陵這一地形起伏度變化范圍內，其面積占比出現了明顯的下降趨勢，進入山地以后，在小、中、大起伏山地范圍內，保持低位穩定，而至極大起伏山地則幾乎沒有分布.耕地偏多型以及無耕地或聚落分布的情形則表現出與前兩者不同的變化趨勢：(a) 以臺地為界，耕地偏多型的面積占比呈現出先上升后下降，然后趨于穩定，并最終消失的緩波峰變化形態；(b) 而無耕地或聚落分布的情形，則是以丘陵為界，先后經歷了急劇上升、平緩上升、急劇下降、平緩下降直至消失的陡波峰變化趨勢.

2) 2011年，耕地偏少型、耕地—聚落均衡型、耕地偏多型以及無耕地或聚落分布的面積分別占全區總面積的41.05%、24.95%、3.72%、30.28%，其中，耕地偏少型所占比例最大，耕地偏多型比例最小.從不同地形起伏度下的耦合類型來看，平原地區耕地偏少型所占比例最大為23.08%，耕地—聚落均衡型、無耕地或聚落分布分別為17.72%、1.17%，耕地偏多型所占比例最小為0.99%.臺地、丘陵地區無耕地或聚落分布、耕地偏少型所占比例較大，其它類型所占比例則相對較小.在山地地區4種類型均有分布，但所占比例較小且相差不大.

3) 2011年與2005年相比，總體而言，耕地—聚落不同耦合類型在不同地形起伏度上均有變化，但變化幅度又不盡相同.耕地偏少型與耕地—聚落均衡型變化幅度比較大，分別由2005年的30.37%、33.35%變為2011年的41.05%、24.95%，而耕地偏多型與無耕地或聚落分布型的變化則相對比較小.從不同的地形起伏度下的耦合類型看，以小起伏山地為界，平原地區耕地偏少型增加幅度最大為7.33%，耕地—聚落均衡型減小幅度最大為6.63%，其它兩種類型變化幅度較小.臺地與丘陵地區耕地偏少型都處于增加趨勢，而其它3種類型都處于減少趨勢.進入山地后，隨著地形起伏度的增大，耕地—聚落耦合各類型的變化趨勢雖不盡相同但變化的幅度都不大，耕地偏少型在不同起伏度山地地區都處于增加趨勢，而無耕地或聚落分布類型則都處于減少趨勢，耕地偏多型與耕地—聚落均衡型的變化趨勢則不具規律性.

4 結論

本文利用都江堰市2005年與2011年的耕地、聚落數據，基于地形起伏度，分析了都江堰市耕地、聚落的時空分異特征，并在此基礎上探究了耕地—聚落耦合類型的時空演變特征，結論如下：

1) 都江堰市的地形起伏度空間分布差異明顯，平原、臺地和丘陵面積占全市面積的96%以上且主要分布在都江堰市的東南部地區，而山地面積不足全市面積的4%，主要分布在都江堰市的北部及西部地區.

2) 都江堰市的耕地主要分布在東南部的平原地區，占全市耕地總面積的89%以上，臺地、丘陵地區次之，隨著海拔的不斷升高，耕地的分布面積不斷減少，2011年與2005年相比，耕地的變化也主要發生在平原、臺地和丘陵地區，其它地形起伏度單元下的耕地則變化不大，表明地形對耕地資源的分布具有顯著影響.

3) 都江堰市超過99%的聚落分布在平原、臺地與丘陵地區，在山地地區雖有聚落分布，但面積很小，且隨著時間的推移，在各地形起伏度地形區，聚落的面積都處于增加趨勢.

4) 從空間分布看，2005年與2011年都江堰市耕地—聚落耦合類型發生了顯著的變化，耕地偏少類型區主要分布在城區及建制鎮且不斷向外擴張，耕地—聚落均衡類型區主要分布在都江堰城區周邊且分布面積明顯減少，耕地偏多類型區主要分布在都江堰市城區的東北、西北及西南部地區，無耕地或聚落分布類型區主要分布在都江堰市北部地區，且面積都有所減少；從耕地—聚落不同耦合類型在不同地形區上的變化結果及趨勢看，耕地—聚落不同耦合類型在平原、臺地、丘陵地區變化較明顯，耕地偏少型都處于增加趨勢，而其它類型則都處于減少趨勢，在小起伏山地至極大起伏山地地區，各類型的變化幅度較小且增加或減少的趨勢比較復雜.

[1] 王華,甘昉. 基于洛倫茨曲線及景觀格局指數的江漢平原地區耕地利用分區研究[J]. 國土與自然資源研究,2013(3):39-42.

[2] 田光進,莊大方,劉明亮. 近10年中國耕地資源時空變化特征[J]. 地球科學進展,2003,18(1):30-36.

[3] 梁海鷗. 巴彥縣耕地資源利用變化時空特征研究[D]. 哈爾濱:東北農業大學,2012.

[4] 劉淼,胡遠滿,常禹,等. 崛江上游汶川縣耕地景觀變化及空間驅動力[J]. 應用生態學報,2007,18(3):569-574.

[5] 關興良,方創琳,魯莎莎. 中國耕地變化的空間格局與重心曲線動態分析[J]. 自然資源學報,2010,25(12):1997-2006.

[6] 張孝宇,謝新朋,張安錄. 武漢市耕地非農化的空間非均衡發展與空間擴散路徑分析[J]. 自然資源學報,2014,29(10):1649-1659.

[7] 劉桂林,張落成,張倩. 1985—2010年南京市耕地變化軌跡及驅動力分析[J]. 生態與農村環境學報,2013,29(6):688-694.

[8] 趙曉麗,張增祥,汪瀟. 中國近30 a耕地變化時空特征及其主要原因分析[J]. 農業工程學報,2014,30(3):1-11.

[9] 任平,吳濤,周介銘. 基于耕地保護價值空間特征的非農化區域補償方法[J]. 農業工程學報,2014,30(20):277-287.

[10] 王雨濛. 土地用途管制與耕地保護及補償機制研究[D]. 武漢:華中農業大學,2010.

[11] 柯新利,楊柏寒,丁璐,等. 基于目標責任區際優化的耕地保護補償[J]. 中國人口·資源與環境,2015,25(1):142-151.

[12] 陳永林,謝炳庚,李曉青,等. 長株潭地區聚落空間演化及其與耕地的空間關系研究[J]. 人文地理,2015(6):106-112.

[13] 楊雪婷,蘇維詞. 云南省縣域地形起伏度與人口分布的耦合關系研究[J]. 重慶師范大學學報(自然科學版),2016,33(1):123-129.

[14] 唐飛,陳曦,程維明,等. 基于DEM的準噶爾盆地及其西北山區地勢起伏度研究[J]. 干旱區地理,2006,29(3):388-392.

[15] 陳學兄,常慶瑞,郭碧云,等. 基于SRTM DEM數據的中國地形起伏度分析研究[J]. 應用基礎與工程科學學報,2013,21(4):670-678.

[16] 張軍,李曉東,陳春艷,等. 新疆地勢起伏度的分析研究[J]. 蘭州大學學報(自然科學版),2008,44(S1):10-13,19.

[17] 中國科學院地理研究所. 中國1:1 000 000地貌圖制圖規范(試行)[M]. 北京:科學出版社,1987:33-34.

[18] 甘彩紅. 基于拼地與泉落女問拓合的本區度地人地關系研究[D]. 重慶:重慶師范大學,2015.

[19] 新華網. 國家新型城鎮化規劃[EB/OL]. http://politics.people.com.cn/n/2014/0317/c1001-24649809.html.

[20] 黃瑩,包安明,陳曦,等. 基于綠洲土地利用的區域GDP公里格網化研究[J]. 冰川凍土,2009,31(1):58-65.

[21] 熊俊楠,韋方強,蘇鵬程,等. 基于多源數據的四川省GDP公里格網化研究[J]. 應用基礎與工程科學學報,2013,21(2):317-327.

[22] 王黎明,張少輝,張大泉. 基于公里格網的中國人地關系緊張度定量模擬[J]. 地理研究,2007,26(3):425-432.

(編輯 周 俊)

Research on the Development of the Coupling Relationship between Cultivated Lands and Settlements Based on Relief Amplitude: Take the Case of Dujiangyan City

WANG Huihao1,2, REN Ping1,2, ZHANG Zhibo1,2

(1.KeyLaboratoryofLandResourcesEvaluationandMonitoringinSouthwestofMinistryofEducation,SichuanNormalUniversity,Chengdu610066,Sichuan; 2.GeographyandResourcesScienceCollege,SichuanNormalUniversity,Chengdu610066,Sichuan)

Based on the extracted relief amplitude from DEM data and the data of cultivated lands and settlements in 2005 and 2011 of Dujiangyan City, and the use of ArcGIS software, this paper gave a quantitative analysis of the spatiotemporal distribution and its difference law of cultivated lands, settlements and the coupling type of settlements and cultivated lands. The results are as follow: 1) among the relief amplitude zone in Dujiangyan City, the plain covers the largest area, and mesa and hills second, and the mountains is the smallest; 2) in the plain, the area of cultivated lands and settlements is the largest, and decreases with the increase of relief amplitude; 3) the urban and the organic town have the less cultivated lands which mainly distribute in the mesa area and the balanced type of settlements and cultivated lands are mostly located in the surroundings of Dujiangyan City; 4) in the interior of each kind of landscape in every relief zone, take the low relief amplitude mountain as a boundary, the different coupling patterns of cultivated lands and settlements in plain, mesa and hills have a significant change. The type of more cultivated lands tends to increase and other types are decreasing. From the low relief amplitude mountain to the higher, each type varies within a small range and has a complex changing trend.

relief amplitude; cultivated lands; settlement; spatial coupling; Dujiangyan City

2016-06-18

國家自然科學基金(40301196)

K901

A

1001-8395(2017)04-0536-08

10.3969/j.issn.1001-8395.2017.04.018

*通信作者簡介：任 平(1978—),男,教授,主要從事國土資源管理與生態環境效應的研究,E-mail:pren121680@126.com