袁 媛,羅志軍,*,齊 松,張 田,劉夢婷,許 巾,劉洋洋
1 江西農業大學 國土資源與環境學院, 南昌 330045 2 江西省鄱陽湖流域農業資源與生態重點實驗室, 南昌 330045 3 贛州市自然資源局, 贛州 341000
土地資源是人類賴以生存并為人類生產生活提供所需基礎物質的不可再生資源[1- 2],其可持續利用是社會可持續發展的關鍵[3- 4]。土地適宜性是土地利用情況和土地利用規劃的基本依據[5],目前對土地適宜性的研究多集中于損毀礦山的土地復墾適宜性評價[6- 7],而對建設用地開發適宜性研究少,且研究區多為山地等地勢較高的地區,如趙筱青等[8]對云南喀斯特山區進行了建設用地開發適宜性評價及分區研究。因此,對于經濟社會快速發展、城鎮化進程加快和城鄉土地利用結構失衡問題突出的城市地區,亟需開展以城鄉建設用地開發適宜性評價為內容的相關研究。
土地承載力是指在一定的物質生活水準下,土地資源所能承載的最大人口和社會活動規模[9-10]。20世紀80年代末,我國開始展開圍繞人、地、糧關系的土地資源承載力研究[11-13],隨著土地開發程度的加深[14-15],土地利用空間不均衡等問題突出,相關研究逐漸拓展到土地綜合承載力[16-18]。近年來,隨著土地承載力評價受到愈來愈多的關注,其評價對象、評價方法、評價體系等都愈加豐富。評價對象方面,主要為城市土地承載力[19]、區域土地承載力[20-22]以及特殊地區土地承載力[23-24]。研究方法上,系統動力學[25]、生態足跡模型[26]、多指標綜合評價[27- 28]等方法被廣泛應用到不同地區的土地資源承載力評價。從評價體系來看,土地承載力評價主要從環境容量[29]、生態承載力[30]、自然基礎要素[31]等方面選取評價指標,較少考慮到區位條件及人口密度等社會經濟因素,無法反映在人類活動對土地承載力的影響。
不同學者分別從不同方面開展了土地適宜性和土地承載力評價研究,但以土地適宜性評價為基礎進行城鄉建設用地承載力研究還較少,更缺少城鄉建設用地承載力和現狀開發狀態之間的空間匹配關系研究。因此,理清城鄉建設用地承載力和現狀開發狀態之間的空間匹配關系,對于合理利用土地和優化國土空間布局,以及加快生態文明建設和實現區域可持續利用具有重要的意義[32-34]。以江西省南昌市為例,在城鄉建設用地開發適宜性評價基礎上,測算城鄉建設用地承載力和可承載城鄉建設用地豐度,并探討城鄉建設用地承載力和現狀開發狀態之間的空間匹配關系,以期為南昌市城鄉建設用地合理開發利用和國土空間優化提供參考。
南昌市是江西省的省會,地處江西中部偏北,位于115°27′—116°35′E、28°09′—29°11′N之間,總面積達7402.36 km2,城鄉建設用地面積為1237.57 km2,占總面積的16.72%,下轄五區(東湖區、西湖區、青云譜區、灣里區、青山湖區)和四縣(南昌縣、新建縣、安義縣、進賢縣),地理位置如圖1所示。境內地形東南平坦,西北丘陵起伏,氣候屬于亞熱帶濕潤季風氣候,光照充足,雨水充沛,自然資源豐富。近年來,南昌市經濟社會發展迅猛,2018年全市生產總值達到5274.67億元,常住人口達554.55萬人,城鎮化率74.2%。

圖1 南昌市地理位置示意圖Fig.1 Geographical location map of Nanchang City
本研究數據主要包括:南昌市2017年土地利用變更調查數據庫、永久基本農田劃定成果、土地利用總體規劃(2006—2020年)數據庫、生態紅線劃定成果等,還有江西省地面塌陷分布圖、地質災害調查評價圖、南昌市統計年鑒、道路矢量數據等,以及通過地理空間數據云(http://www.gscloud.cn/)獲取的數字高程模型(DEM)(空間分辨率為30 m)。通過ArcGIS 10.2軟件,將數據投影到統一坐標系下,并劃分為30 m30 m的格網作評價單元,以便對數據進行空間分析。
隨著經濟快速發展和城鎮化不斷加快,大量耕地開發為建設用地,國土空間開發利用格局發生明顯重構,自然資源的過度開發以及土地利用不合理等使區域空間開發失衡問題日益凸顯。為解決區域土地利用空間失衡問題,協調人口、資源和社會發展的關系,需要開展城鄉建設用地承載力和現狀開發狀態之間的空間匹配關系研究。本文研究思路如圖2所示,具體如下:
(1)指標選取及量化處理:適宜性評價從生態安全、自然條件、社會經濟等方面選取11個評價指標構建評價體系,生態紅線、永久基本農田區、地面塌陷為強限制性指標,高程、蓄滯洪區、土地利用類型、區位條件為較強限制性指標,坡度大于25°為強限制性指標,小于25°為較強限制性指標,并對指標賦值,進行量化處理。
(2)劃分適宜性等級:以ArcGIS 10.2軟件為處理平臺,進行建設用地開發適宜性評價,并通過SPSS 22.0軟件對適宜性分值進行聚類分析,劃分適宜性等級。
(3)計算城鄉建設用地承載力:利用ArcGIS 10.2軟件進行空間疊加分析,計算城鄉建設用地可承載臨界值和可承載城鄉建設用地豐度。
(4)空間匹配度分析和相關分析:計算城鄉建設用地空間匹配度,并將城鄉建設用地可承載臨界值、可承載城鄉建設用地豐度和空間匹配度進行相關分析,研究現狀利用狀態與其合理開發狀態之間的關系。

圖2 南昌市城鄉建設用地承載能力評價及空間匹配研究技術路線圖Fig.2 Technical roadmap of carrying capacity evaluation and spatial matching of urban-rural construction land in Nanchang
城鄉建設用地開發適宜性評價是指評定一定條件下土地開發為建設用地的適宜性程度[35]。根據數據的可獲取性和南昌市實際,從生態安全、自然條件、社會經濟等3方面選取11個評價指標,構建城鄉建設用地開發適宜性評價指標體系,并將指標分為強限制性指標和較強限制性指標兩類,如表1所示。
在評價時,考慮到經濟發展與生態保護的矛盾,將生態紅線作為強限制性指標處理。生態紅線是指為維持生態安全和生態系統完整性,通過科學合理的方法劃定的必須實施特殊保護的國土空間界線[36]。同時為貫徹最嚴格的耕地保護政策,選取永久基本農田區為強限制性指標。地面塌陷區開發難度大且多為自然災害易發區,所以將地面塌陷區也列為強限制性指標。此外,坡度直接影響建設用地開發成本和開發難度,坡度過大會影響建筑物的穩定性及增加潛在危險性。相關研究表明,在進行建設用地開發時,坡度大于25°的陡坡地一般被視為“禁止建設區”[37-38]。地質災害直接威脅到人類的生命財產安全,區位條件和人口密度影響社會經濟發展以及信息的可達性和居民生活的便利程度。因此,在適宜性評價時,都將其作為評價指標考慮在內。
對于強限制性指標,根據“一票否決”原則進行0和1賦值。對于較強限制性指標,在參考相關研究[39-41]的基礎上,根據適宜性程度進行0—100賦值,采用層次分析法和專家打分法確定各較強限制性指標的權重。由于坡度對建設用地開發影響特殊,對于坡度大于25°的陡坡地進行0和1賦值,而坡度小于25°的地區則按限制性等級進行0—100賦值(見表1)。
采用限制系數評價模型,計算各評價單元的城鄉建設用地開發適宜性分值,計算公式如下:
(1)
式中E為綜合適宜性分值;i為強限制性指標序號;j為較強限制性指標序號;Fi為第i個強限制性指標適宜性分值;fj為第j個較強限制性指標適宜性分值;wj為第j個較強限制性指標權重;m為強限制性指標個數;n為較強限制性指標個數。

表1 南昌市城鄉建設用地開發適宜性評價指標及其劃分標準
城鄉建設用地可承載臨界值是指區域土地資源所能承載的人類活動的最大規模和強度。基于城鄉建設用地開發適宜性評價結果,將最適宜開發建設區視為最適宜承載區,較適宜開發建設區視為較適宜承載區,結合現狀建設用地情況,通過ArcGIS 10.2軟件,測算城鄉建設用地可承載臨界值。計算公式如下:
(2)
式中,LDI為城鄉建設用地可承載臨界值,E1為最適宜承載區面積,E2為較適宜承載區面積,C為現狀城鄉建設用地面積,U為研究區總面積。
城鄉建設用地已開發強度反映區域現狀建設用地利用狀況,以各縣(區)現狀城鄉建設用地總量(C)占各縣(區)區域總面積(UT)的比例表示。
DI=C/UT
(3)
以城鄉建設用地可承載臨界值為城鄉建設用地極限開發強度,將其與已開發強度相減,差值為剩余開發強度,即城鄉建設用地的開發潛力。
可承載城鄉建設用地豐度(A)反映的是研究區內可承載建設用地的豐富程度[42],主要從數量和空間分布情況考慮。數量方面(A1),以各縣(區)可承載城鄉建設用地面積(E1+E2)占縣域總面積(UT)的比例表示;空間方面(A2),以各縣(區)可承載城鄉建設用地面積(E1+E2)占研究區可承載城鄉建設用地面積(US)的比例表示,數據標準化處理后求和得到各縣(區)可承載城鄉建設用地豐度,計算公式如下:
(4)