沿海地區耕地占補時空演變研究
——以寧波市為例

張志軍，趙志剛，邱俊武

(1.天津市測繪院，天津 300381； 2.武漢大學資源與環境科學學院，湖北武漢 430079)

沿海地區耕地占補時空演變研究
——以寧波市為例

張志軍1?，趙志剛2，邱俊武1

(1.天津市測繪院，天津 300381； 2.武漢大學資源與環境科學學院，湖北武漢 430079)

目前，建設占用耕地，后備資源補充模式是被占用耕地得到補償的一種有效方式，然而這種單一模式是否長期有效是值得深思的一個問題。本文使用Logistic回歸模型和二分類因變量的多元化分析方法對此模式在沿海地區的時空演變進行研究。以寧波市為例，模擬了建設用地、耕地和后備資源三者在時空演變中的關系，使用現實數據對模型進行了驗證，證明其可行性和有效性。最后預測了2011年～2021年建設占用耕地和后備資源補充耕地的趨勢。結果表明建設用地擴展將向南部遷移，占用耕地數量相對減少，后備資源將會枯竭等后果。

Logistic回歸模型；多元量化分析；時空模擬；耕地占補平衡

1 引 言

隨著1996年耕地占補平衡政策提出后，中國大部分省市特別是沿海地區在解決這個問題上遇到重大的挑戰。目前，建設占用耕地，后備資源補充的模式被認為是應對這種挑戰的一種有效方式之一[1，2]，該模式在沿海地區的耕地占補平衡中初見成效。然而該模式是否長期有效，能否滿足耕地占補平衡的要求和建設用地擴張的趨勢是一個值得探討的問題。本文將以寧波市為例對此模式下的建設用地、耕地和后備資源的時空演變進行研究，通過對2011年和2021年狀況進行預測，對該模式的可持續性進行探討。

目前，建設用地擴張和耕地被占用的研究方法基本都是基于LUCC的研究框架下采用不同的模型開展的，大致分4大類:統計模型、動態模型、綜合模型和馬爾柯夫模型[3～7]。其中統計模型的優勢在于:簡化問題，容易抓住復雜系統的主要矛盾[8]，更適合本文的研究，因此本文采用Logistic回歸模型作為時空演變模型[9，10]。與此同時，耕地占補平衡分為3個層次:一是數量平衡，二是質量平衡，三是生產力平衡[11]。本文將“與淡水域距離”等衡量耕地質量和生產力因子加入到演變模型中，使用二分類因變量的多元量化分析方法，分別對建設占用耕地和后備資源補充耕地進行模擬，通過對模擬結果的分析，驗證了模型的可行性和有效性。

2 研究區域

都很快的城市，其城市化水平從1977年的13.6%提高到1998年的41%，到2005年已達到56%，表現出了強勁的增長態勢。這與它的區位、自然條件和社會經濟條件是分不開的。寧波在地勢上西南高，東北低，自西南向東北傾斜，三面環海、一面依山，地貌以平原和丘陵為主。耕地和建成區基本分布在平原地區，且呈現出耕地包圍建設用地的圈層結構。

3 模型和數據

其中:P是因變量，表示事件發生的概率，X1，X2……Xm為自變量，參數α為常數，與參數β1，β2……βm是回歸系數，βi決定Xi的影響大小和敏感性[12]。

數據是寧波市域土地利用分類圖，分為因變量和自變量(驅動因子)。因變量是指建設用地、耕地和后備資源，自變量包括DEM、坡度、道路、水系、人口密度、GDP、城市規劃圖等。建設用地、耕地和自變量為自1991年～2005年歷史數據，后備資源為自1998年整理入庫的數據。數據格式為柵格，格網大小為100 m×100 m。

Logistic回歸模型普遍應用于針對二分類因變量(y=1或y=0)的多元量化分析方法。利用該模型可以研究影響事件發生的概率、二分性變量的主要因子，模型表達式如式(1):

寧波是我國沿海地區一個經濟和城市化發展速度

4 研究方法

本文的研究路線如下:首先建立建設用地擴展的時空演變模型，預測2011年和2021年建設用地分布情況，從而計算出2011年和2021年建設用地占用耕地的情況；然后建立后備資源補充占用耕地的時空演變模型，從而得到2011年和2021年的補充情況。同時建立相應的評價模型對時空演變過程的有效性進行驗證。

4.1 建設用地占用耕地時空演變

(1)指標選取

本文采用二分類因變量多元分析方法，使用SPSS軟件的Hosmer-Lemeshow指標評價篩選結果[13，14]。其中自變量選取了7個驅動因子:與道路的距離(x1)、與組團中心(縣級以上)的距離(x2)、人口密度(x3)、人均GDP(x4)、坡度圖(x5)、DEM(x6)和城市規劃因子(x7)；因變量選取1991年～2005年占用建設用地作為二分類因變量。各變量的空間分布如圖1所示。

圖1 時空演變模型中各變量空間分布圖

(2)禁建區域

在時空演變中需考慮禁建區，即不可建設擴展的區域。本文中的演變模型將各類飲用水源水庫(包括水庫、湖泊、大型河流)、坡度大于25°或高程大于200 m的區域、濕地保護區、風景名勝區等作為禁建區。

(3)建設用地占用耕地時空演變

首先進行數據采樣，因變量將y=1和y=0各取一半隨機采樣；然后利用Logistic模型進行演變模擬，得到預測分類結果，如表1所示；為規避隨機采樣帶來的影響，本文采集6套數據，并對結果進行比較，詳情見表2；最后選取預測率最高的反演模擬公式，如式(2)。

預測分類與原始數據分類表 表1

采樣擬合匯總表 表2

(4)驗證模型

以2005年～2007年數據進行預測，并與真實數據進行比較，結果如圖2所示。結果顯示預測建設占用面積為137.566 7 km2，真實面積為136.107 km2，因此該演變模型視為有效。

(5)建設用地占用耕地預測

基于1991年～2005年的歷史數據，采用時序擬合的方法對2011年和2021年占用建設用地情況進行預測，結果表明2005年～2011年需新增建設用地309.80 km2，占用耕地241.65 km2；2011年～2021年需新增建設用地377.89 km2，占用耕地294.75 km2，其空間分布如圖3所示。

圖2 2005年～2007年建設用地擴展預測圖和真實圖

圖3 預測2005年～2011和2005年～2021年建設用地擴展和占用耕地空間分布圖

(6)結果分析

對建設用地占用耕地情況進行分析，得到如下結論:

(1)2011年～2021年，寧波市城市擴張主要分布在中心城區和北部余慈地區，南部的寧海和象山縣的新增建設用地有增加的趨勢。

(2)預測期間新增建設用地75%左右來自于耕地，尤其集中于公路和居民點周邊等區位條件較優的肥沃的耕地。

4.2 后備資源補充占用耕地時空演變

(1)指標選取

經過驅動力分析，后備資源補充占用耕地時空演變模型選取6個指標作為自變量驅動因子:與道路的距離(x1)、與組團中心的距離(x2)、人口密度(x3)、坡度圖(x5)、DEM(x6)，與淡水域的距離(x8)，如圖4所示。因變量為1998年～2005年后備資源已開發為耕地和未開發為耕地的分布圖，如圖5所示。

圖4 與淡水域距離 圖5 1998年～2005年后備補充耕地圖

(2)模型擬合和驗證

后備資源補充占用耕地時空演變模型公式為:

驗證模型與建設用地占用耕地的時空演變驗證模型類似的，在此不做過多贅述。

(3)后備資源補充占用耕地預測

根據前面預測得到的2011年建設用地占用耕地的面積(249.95 km2)，確定2011年后備資源庫補充耕地的概率閾值(0.68)，對2021年后備資源補充占用耕地進行預測。發現后備資源已不能全額補充建設占用耕地，如采取占一補半的策略(即概率閾值取0.5)，后備資源分布情況如圖6所示。

圖6 2011年后備資源補充耕地和2021年補半分布圖

(4)結果分析

對后備資源補充占用耕地情況進行分析，得到如下結論:

①預測閾值比建設占用耕地時空演變模型偏低。

②2021年以后，寧波市后備資源潛力基本消失。

③在后備資源中，沿海海涂面積較大，其他資源地塊面積較小，且分布零散。

5 結論與對策

5.1 結論

在后備資源補充建設占用耕地模式下，基于logistic模型分別對寧波市建設用地占用耕地和后備資源補充占用耕地進行時空演變模擬，對模擬結果進行分析后得到如下結論:

(1)寧波市建設用地擴展速度很快，建設用地在區位和土質上都偏重于利用率高、效益好的耕地，使得耕地逐年呈現負增長的態勢。

(2)從模擬建設用地占用耕地的時空圖可直觀地發現，到2021年寧波市“中心”耕地擴展空間所剩無幾，南部建設占用耕地開始增長。

(3)從模擬后備資源補充耕地的結果發現，到2021年寧波市將無地可補，后備資源補充建設占用耕地單一模式無法持續維持耕地占用平衡，亟須尋求新的解決方案。

5.2 對策

對于如何維持耕地占用平衡問題，筆者提出了如下參考對策:

(1)從政府的角度，完善法規制度，制定更加科學合理的政策措施；

(2)從經濟的角度，在政策允許的范圍內通過耕地占用權交易，可緩解部分問題；

(3)從效益的角度，提高耕地的利用效益來彌補占補不平衡問題；

(4)從技術的角度，可嘗試摸索表土剝離利用實現耕地質量占補平衡。

6 總結與展望

本研究力求以圖形表達的方式來簡化建設用地、耕地和后備用地三者之間復雜的演變系統，提供直觀的決策依據。同時，本文提出的時空演變模型可以通過修改邊界條件或閾值等參數，快速模擬出新的結果，具備輔助決策功能。但該模型只考慮了表現耕地質量和生產力的自然因子，在后面的研究工作中需考慮社會因子等。

[1] 孟媛.“堅守18億畝耕地紅線”研討會摘要.中國農業大學學報[J].2007，12(3)：25～26.

[2] 張琳，張鳳榮，薛永森等.中國各省耕地數量占補平衡趨勢預測.資源科學[J].2007，29(6)：114～119.

[3] 史培軍，宮鵬，李曉兵等.土地利用/覆蓋變化研究的方法與實踐[M].北京：科學出版社，2000.

[4] 擺萬奇，趙士洞.土地利用變化驅動力系統分析[J].資源科學，2001，21(3)，39～41.

[5] 韓超峰，陳仲新.LUCC驅動力模型研究綜述[J].中國農業通報，2008，24(4)：365～368.

[6] 郭斌，陳佑啟，姚艷敏等.土地利用與土地覆被變化驅動力研究綜述[J].中國農學通報，2008，(4)：81～86.

[7] 蒙吉軍，嚴汾.大城市邊緣區的LUCC驅動力時空分異研究——以北京冒平區為例[J].北京大學學報(自然科學版)，2008，3：137～144.

[8] 錢銘杰.區域建設用地集約利用模式與評價研究[D].北京：中國地質大學，2011.

[9] 李全，李霖，姜文亮等.基于logistic模型城市擴展時空模擬——以深圳龍崗區為例[J].華中師范大學學報·自然科學版，2006，40(3)：442～446.

[10] James D.A.Millington，George L.W.Perry，and Rau′Romero-Calcerrada.Regression Techniques for Examining Land Use/Cover Change：A Case Study of a Mediterranean Landscape[J].Ecosystems，2007，10：562～578.

[11] 萬江波.基于多源數據的土地利用信息圖譜研究[D].武漢：武漢大學，2008.

[12] 王濟川，郭志剛.Logistic回歸模型方法與應用[M].北京：高等教育出版社，2001.

[13] Bo Huang，L Zhang，and BWu.Spatio-temporal Analysis of Rural-urban Land Conversion[J].International Journal of Geographical Information Science，2009，23(10)：1～27.

[14] 焦鋒，秦伯強.GIS支持下的小尺度土地驅動力研究——以宜興市湖滏小流域為例[J].長江流域資源與環境，2003，12(3)：205～210.

上海市測繪院為上海市政府更新工作用圖

近日，上海市測繪院圓滿完成了為上海市政府領導更新工作用圖的保障任務。

本次工作用圖更新任務周期緊、任務重，從接到通知到完成任務只有短短幾天時間。上海市政府有關部門不僅要求更新地圖信息，還對裝幀提出了嚴格要求。接到任務后，上海市測繪院立即啟動了領導用圖保障機制，各相關部門積極行動，迅速制定更新方案。任務實施過程中，計劃業務處和第四分院按照上海市測繪院“服務政府和服務社會公眾”的宗旨，多次赴實地踏勘測量，精在點上、細到實處，力爭呈現最好的效果，并提前完成了本次工作用圖更新任務，得到了上級的肯定。

多年來，上海市測繪院堅持為上海市委、市政府、市人大、市政協等部門制作和更新工作用圖，及時反映本市最新的城市建設和發展狀況，為領導決策提供測繪地理信息支撐。

(來源：http：//www.shsmi.cn)

Exam ining Temporal and Spatial Evolution for Arable Land Occupied and Compensated at Coastal Areas——A Case Study for Ningbo

Zhang Zhijun1，Zhao Zhigang2，Qiu Junwu1
(1.Tianjin Institute of Surveying and Mapping，Tianjin 300381，China；2.School of Resource and Environment Science，Wuhan University，Wuhan 430079，China)

Themode thatusing back-up resources to supply farmland thatused for construction is an effectiveway of compensation for occupied farm land at presently.But is it continuable?This paper studied the temporal and spatial evolution at coastal areaswith the logistic regression model and the two categories dependent variable analysismethod.Take Ningbo as an example，we simulated space-time evolution of the relationship among the construction land，cultivated land and reserve resources.Using real data validated model to prove its feasibility and effectiveness.And lastly，predicted the trend that construction occupying arable land and land reserving resource supplementing arable land in 2011and 2021. Then，the result showed the construction land would expanse more quickly in the south，the number of occupied arable land would reduce relatively，and the reserving resource would exhaust in 2021.

Logistic regressionmodel；Two categories dependent variable analysis；Spatiotemporal simulation；Requisition-compensation quality balance of arable land

1672-8262(2013)02-30-05

P208.2

A

2012—05—18

張志軍(1983—)，男，工程師，主要從事智慧城市及城市發展理論與技術相關研究。