程紹鉑 ，楊桂山 ，吳建楠 ，胡 艷

（1.安徽大學(xué)經(jīng)濟學(xué)院，安徽合肥 230601；2.中國科學(xué)院南京地理與湖泊研究所，江蘇南京210008；3.安徽大學(xué)長三角經(jīng)濟與社會發(fā)展研究中心，安徽合肥230601；4.中山大學(xué)廣東決策科學(xué)研究院，廣東廣州510000）

如何在現(xiàn)有的國情條件下和日益開放的國際環(huán)境下，在實現(xiàn)工業(yè)化和城市化過程中保持農(nóng)業(yè)可持續(xù)發(fā)展保護有限的耕地資源，是中國經(jīng)濟和社會可持續(xù)發(fā)展的全局性戰(zhàn)略問題[1]。中國政府在歷經(jīng)1978—1985年耕地保護制度恢復(fù)期、1986—1997年的耕地保護制度發(fā)展期和1998至今的以國土資源部成立為標志的耕地保護制度完善期[2]，先后出臺了一系列耕地保護政策，在相當大的程度上實現(xiàn)了耕地保護目標，也為社會經(jīng)濟尤其是農(nóng)業(yè)生產(chǎn)的可持續(xù)發(fā)展提供了必要的保障[3]。耕地總量動態(tài)平衡政策作為耕地保護重要政策之一，在一定程度上保證了耕地總量的平衡，有助于遏制耕地加速流失。但是，由于目前中國正處于城市化發(fā)展和經(jīng)濟發(fā)展的快速時期，對建設(shè)用地的需求不斷增長，給耕地總量動態(tài)平衡政策的實施帶來了巨大的挑戰(zhàn)。由于國情和政治體制的不同，國外對耕地總量動態(tài)平衡邊界問題的研究較少。國內(nèi)由黃賢金等率先開展了耕地總量動態(tài)平衡邊界研究[1,4]，但主要集中于區(qū)域尺度。實際上，在更小的空間尺度上研究耕地總量動態(tài)平衡邊界問題，更有助于加強對耕地總量動態(tài)平衡實現(xiàn)規(guī)律的認識，有助于地方政府對不同土地經(jīng)濟發(fā)展階段的把握，從而制定出適合于當?shù)匕l(fā)展情況的耕地保護對策。本文以昆山市為研究對象，研究在現(xiàn)有耕地保護政策實施力度保持不變趨勢下的耕地總量動態(tài)平衡邊界，有助于從實際情況出發(fā)，衡量耕地總量動態(tài)平衡政策的實施效果，增強政策的可操作性，高效地實現(xiàn)耕地保護的目標。同時，本研究對中國其他正處于經(jīng)濟高速發(fā)展面臨耕地保護問題的縣級市均有借鑒意義。

1 研究區(qū)概況

昆山市地處江蘇省東南部，土地總面積928 km2，其中水域面積占29.92%。昆山市經(jīng)濟發(fā)達，據(jù)昆山市政府網(wǎng)站公布數(shù)據(jù)顯示，2010年全市地區(qū)生產(chǎn)總值達到2100.28億元，成為國內(nèi)首個GDP總量超過2000億元的縣級市，同時按常住人口計算人均地區(qū)生產(chǎn)總值達到14.45萬元。改革開放以來，昆山市工業(yè)化和城市化進程迅速，成為以開放型經(jīng)濟為主導(dǎo)，三大產(chǎn)業(yè)協(xié)調(diào)發(fā)展的沿海新興城市之一。在國家統(tǒng)計局公布的全國百強縣（市）排名中，昆山市位居全國百強縣之首。昆山市高速的經(jīng)濟發(fā)展引起了快速的土地利用變化，成為一個典型的縣市級尺度的研究區(qū)域，為揭示經(jīng)濟高速發(fā)展和城市化過程中的耕地變化提供了一個難得的案例[4]。分析經(jīng)濟社會快速發(fā)展的昆山市耕地總量動態(tài)平衡邊界，對其他城市的耕地平衡政策的實施具有啟示性意義。

2 昆山市耕地變化態(tài)勢分析與影響因素分析

根據(jù)昆山市統(tǒng)計年鑒資料，昆山市1985—2009年間耕地數(shù)量數(shù)量變化如圖1所示。昆山市耕地總量1985—2009年間呈減少趨勢，25年間總共減少了68.2%，2009年耕地面積僅相當于1985年耕地面積的32.8%。除1999年耕地面積略有上升外，其他年份耕地面積均不斷下降，耕地面積的減少是一個長期積累的過程。

圖1 昆山市耕地數(shù)量變化趨勢圖Fig.1 The trend of cultivated land change

耕地面積變化受多種因素影響，昆山市自然因素對耕地面積的影響較小，耕地面積變化主要受政策因素和社會經(jīng)濟因素影響。根據(jù)昆山市實際情況以及資料的可獲得性和可靠性，并且參考已有學(xué)者的研究成果[5-8]，選用總戶籍人口X1、非農(nóng)業(yè)人口X2、地區(qū)生產(chǎn)總值X3、人均地區(qū)生產(chǎn)總值X4、第二產(chǎn)業(yè)生產(chǎn)總值X5、糧食總產(chǎn)量X6、社會固定資產(chǎn)投資額X7、城市化率X8等8種指標來反映社會經(jīng)濟變化對耕地面積的影響，數(shù)據(jù)均來自于昆山市統(tǒng)計年鑒；選用耕地保護政策實施強度X9反映耕地面積對保護政策實施力度的反饋。為了量化政策因素對耕地面積的影響，按照已有專家的評價方法[9]并根據(jù)昆山市的實際情況，將不同時期耕地保護政策實施的強度分別定義：將1985年前耕地保護政策強度賦值為“-1”，此時為弱耕地保護狀態(tài)；將1986—1994年賦值為“-3”，這一時期開始實施《土地管理法》，可以定義為較強耕地保護狀態(tài)；將1995—2003年賦值為“-5”，開始實施基本農(nóng)田保護區(qū)制度，定義為強耕地資源保護狀態(tài)；將2003年后定義為“-8”，雖然自1997年開始實施土地用途管制制度，但結(jié)合昆山市具體情況以及考慮制度實施的時間滯后問題，將2003年后定義為最嚴格的耕地資源保護狀態(tài)。

從昆山市社會經(jīng)濟因素和政策因素驅(qū)動變量與耕地數(shù)量的相關(guān)分析可以看出，各驅(qū)動變量與耕地數(shù)量均呈極顯著相關(guān)，說明所選指標的合理性，這些指標能夠反映耕地變化的驅(qū)動。其中，總戶籍人口、非農(nóng)業(yè)人口、地區(qū)生產(chǎn)總值、人均地區(qū)生產(chǎn)總值、第二產(chǎn)業(yè)生產(chǎn)總值、社會固定資產(chǎn)投資額和城市化率與耕地數(shù)量呈負相關(guān)，即這些變量增加會導(dǎo)致耕地數(shù)量的減少。糧食總產(chǎn)量和耕地保護政策實施強度與耕地數(shù)量呈正相關(guān)，即這些變量值越大，耕地面積越大。從耕地保護政策實施強度與耕地數(shù)量呈現(xiàn)極顯著相關(guān)，且相關(guān)系數(shù)達到0.937，可以看出，耕地保護政策實施力度的增強能夠有效遏制耕地數(shù)量的減少。

表1 昆山市耕地數(shù)量與驅(qū)動變量的相關(guān)分析結(jié)果Tab.1 Correlation result of cultivated land and driving factors of Kunshan

以耕地數(shù)量為因變量，表1中所列的驅(qū)動變量為自變量，對耕地數(shù)量進行線性回歸擬合。得到昆山市耕地數(shù)量變化模型為：

模型的復(fù)相關(guān)系數(shù)為0.997，調(diào)整后的復(fù)相關(guān)系數(shù)的平方為0.991，表示自變量總體對因變量的解釋程度為99.1%，也就是說，所選的自變量幾乎完全解釋了耕地數(shù)量的變化。模型的顯著性檢驗中F值為302.874，通過檢驗。實際上，由于所有的變量都是時間序列變量，本身隨時間的變化具有一定的規(guī)律，使得變量間相關(guān)程度增加[10]。因此，抽取其中少數(shù)幾個自變量的組合對耕地數(shù)量進行多元線性回歸，仍然能得到較好的結(jié)果：

模型的復(fù)相關(guān)系數(shù)為0.987，調(diào)整后的復(fù)相關(guān)系數(shù)的平方為0.970，這表示自變量總體對因變量的解釋程度為97.0%，表明地區(qū)生產(chǎn)總值X3、第二產(chǎn)業(yè)生產(chǎn)總值X5和耕地保護政策實施強度X9這3個自變量基本上能夠完全解釋耕地數(shù)量的變化。模型的顯著性檢驗中F值為227.229，通過檢驗?？傮w來看，模型的擬合效果較好。也就是說，這3個自變量能夠反映與耕地數(shù)量變化的伴隨關(guān)系，用地區(qū)生產(chǎn)總值、第二產(chǎn)業(yè)生產(chǎn)總值和政策保護實施強度來表示多因素綜合作用下驅(qū)動耕地數(shù)量變化的機制，是比較符合實際的，用此式替代式1是可行的。

這種驅(qū)動力模型雖然解釋程度較高，能夠通過所選的自變量基本上完全解釋耕地面積的變化，但這是在不考慮耕地變化的邊際收益傾向隨時間變化的前提下得到的，實際上隨著經(jīng)濟的不斷發(fā)展與技術(shù)的進步，耕地變化的邊際收益隨時間不斷減少，呈現(xiàn)出逐漸遞減的趨勢。

以式2中的經(jīng)濟發(fā)展驅(qū)動變量地區(qū)生產(chǎn)總值X3和第二產(chǎn)業(yè)生產(chǎn)總值X5來考察耕地面積變化和經(jīng)濟發(fā)展之間的數(shù)量轉(zhuǎn)換關(guān)系。以“Δ耕地/Δ地區(qū)生產(chǎn)總值”和“Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值”作為經(jīng)濟發(fā)展對耕地資源占用的邊際傾向分析指標，考察這兩個變量隨時間的變化趨勢，通過這一分析可以看出昆山市地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值每增加一單位，耕地數(shù)量的變化情況。

由昆山市耕地變化的邊際占用傾向趨勢分析圖2和圖3可以看出，雖然由于昆山市耕地變化的邊際收益傾向出現(xiàn)階段性的波動，但是總體上遞減的趨勢還是主要的，增加一個單位地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值所需占用的耕地資源總體上呈現(xiàn)出減少的趨勢，表明昆山市經(jīng)濟發(fā)展正逐漸向節(jié)約和集約利用土地方向發(fā)展。

3 昆山市耕地資源邊際占用傾向的經(jīng)濟邊界分析

昆山市地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值的增加對耕地資源的邊際占用傾向是遞減的，在維持目前相對穩(wěn)定高速的經(jīng)濟發(fā)展趨勢與最為嚴格的耕地保護政策實施強度不變的前提下，可以認為當邊際數(shù)值接近于零的時候，耕地資源將不再被占用。此時的耕地面積為耕地資源占用的經(jīng)濟邊界。

3.1 耕地資源的地區(qū)生產(chǎn)總值邊際占用傾向的邊界分析

以耕地數(shù)量為因變量Y，以Δ耕地/Δ地區(qū)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式3中，決定系數(shù)R=0.843，說明方程的擬合效果不錯。當X=0時的得到昆山市耕地面積的地區(qū)生產(chǎn)總值發(fā)展邊界，此時耕地面積為250843.7 hm2。也就是說，在不考慮其他因素的影響下，當耕地面積減少到250843.7 hm2之后，耕地面積就不再隨地區(qū)生產(chǎn)總值的增加而減少。這一面積相當于2006年昆山市耕地面積的83.2%，相當于在2006年的基礎(chǔ)上耕地再減少50727.27 hm2。

以時間為因變量T，Δ耕地/Δ地區(qū)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式4中，決定系數(shù)R=0.826。當X=0時可以等到昆山市耕地面積的地區(qū)生產(chǎn)總值邊際占用停滯時間為2011.8，也就是會說在只考慮地區(qū)生產(chǎn)總值對耕地占用的情況下，在2011—2012年耕地資源將停止被占用。

圖2 昆山市地區(qū)生產(chǎn)總值的耕地資源邊際占用傾向趨勢圖Fig.2 The trend of the marginal effect on cultivated land occupation due to GDP growth

圖3 昆山市第二產(chǎn)業(yè)生產(chǎn)總值的耕地資源邊際占用傾向趨勢圖Fig.3 The trend of the marginal effect on cultivated occupation due to Gop growth in the secondary sector

3.2 耕地資源的第二產(chǎn)業(yè)生產(chǎn)總值邊際占用傾向的邊界分析

以耕地面積為因變量Y，以Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值為自變量X建立邊界分析函數(shù)，得到：

式5中，決定系數(shù)R=0.786。當X=0時，得到昆山市耕地面積的第二產(chǎn)業(yè)生產(chǎn)總值發(fā)展邊界，此時的耕地面積為276933.2 hm2。即在不考慮其他因素的影響下，當耕地面積減少到276933.2 hm2之后，耕地面積就不再隨第二產(chǎn)業(yè)生產(chǎn)總值的增加而減少。這一面積相當于昆山市2006年耕地面積的91.8%，比昆山市地區(qū)生產(chǎn)總值邊界的結(jié)果要小，說明相比于地區(qū)生產(chǎn)總值昆山市第二產(chǎn)業(yè)生產(chǎn)總值更早達到發(fā)展邊界，工業(yè)用地更早的實現(xiàn)土地資本的存量調(diào)劑，零增量占用耕地資源的發(fā)展模式。

以時間為自變量T，Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式6中，決定系數(shù)R=0.843，擬合效果較好。當X=0時可以得到昆山市耕地面積的第二產(chǎn)業(yè)生產(chǎn)總值邊際占用停滯時間為2010.9。也就是會說在只考慮第二產(chǎn)業(yè)生產(chǎn)總值對耕地占用的情況下，在2010—2011年間耕地資源可能會停止被占用。這一停滯時間同樣比由地區(qū)生產(chǎn)總值停滯邊界得到的停滯時間要早，進一步證實了相對于地區(qū)生產(chǎn)總值來說，昆山市第二產(chǎn)業(yè)生產(chǎn)總值更早達到發(fā)展邊界。

3.3 耕地資源占用態(tài)勢分析

通過建立耕地數(shù)量隨時間變化的回歸模型來驗證前面對耕地占用停滯邊界的分析。以年份為因變量T，耕地面積為自變量X建立函數(shù)得到：

式7中，決定系數(shù)R=0.961，擬合效果很好。令耕地面積為地區(qū)生產(chǎn)總值占用邊界值即X=250843.7時得到T=2012.93，這個值與由式4得到的只考慮地區(qū)生產(chǎn)總值對耕地占用的情況下耕地會停止被占用的時間接近。可以認為式3—4與本式得到了互為印證，可以認為昆山市耕地數(shù)量的地區(qū)生產(chǎn)總值邊界值為250843.7 hm2，當耕地數(shù)量減少到此邊界值之后，耕地數(shù)量不再隨地區(qū)生產(chǎn)總值的增加而減少，其達到時間大致為2012年左右。

令耕地面積為第二產(chǎn)業(yè)生產(chǎn)總值占用邊界值即X=276933.2時得到T=2011.2，這個值與由式6得到的只考慮第二產(chǎn)業(yè)生產(chǎn)總值對耕地占用的情況下耕地會停止被占用的時間非常接近?？梢哉J為式5—6與式7得到了相互印證，因此基于本分析可以認為昆山市耕地數(shù)量的第二產(chǎn)業(yè)生產(chǎn)總值邊界值為276933.2 hm2，當耕地數(shù)量減少到此邊界值之后，耕地數(shù)量不再隨地區(qū)生產(chǎn)總值的增加而減少，其達到時間大致為2011年。

4 結(jié)論

單從經(jīng)濟效益的角度來看，耕地的經(jīng)濟產(chǎn)出效益遠小于建設(shè)用地經(jīng)濟效益，所以耕地被占用轉(zhuǎn)變?yōu)榻ㄔO(shè)用地是符合經(jīng)濟理性的，屬于自發(fā)、主動的經(jīng)濟行為。但隨著耕地占用成本的不斷增加、科學(xué)技術(shù)的進步、區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整以及政府耕地保護政策實施力度的加強，耕地占用的邊際收益傾向正在不斷減小，昆山市經(jīng)濟發(fā)展正逐漸向節(jié)約和集約利用土地方向發(fā)展。從地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值來看，昆山市經(jīng)濟發(fā)展將由大量占用土地資本向相對高效占用土地資本轉(zhuǎn)變，且最終將實現(xiàn)土地資本存量調(diào)劑，經(jīng)濟發(fā)展零占用耕地資源的發(fā)展模式。單從地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值發(fā)展角度來說，在維持目前的最為嚴格的耕地保護政策實施強度不變的情況下，昆山市耕地數(shù)量將在24637.8—250843.7 hm2之間達到邊界值。此后地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值的發(fā)展將不再依賴于耕地占用，達到邊界值的時間大致為2011—2012年之間，且第二產(chǎn)業(yè)生產(chǎn)總值的耕地占用比地區(qū)生產(chǎn)總值的耕地占用更早達到邊界值。屆時，昆山市耕地總量動態(tài)平衡有望得到更為高效的實現(xiàn)。

當然，這只是在不考慮耕地省內(nèi)異地占補平衡以及假設(shè)耕地保護政策實施強度維持不變的情況下得出的，事實上耕地占用的實際情況更為復(fù)雜，并且當經(jīng)濟社會進入一個新的發(fā)展時期后，這一預(yù)測結(jié)果將發(fā)生相應(yīng)的變化。從昆山市近25年來耕地面積變化與經(jīng)濟高速增長的相互關(guān)系中可以看出，與經(jīng)濟增長條件下環(huán)境的庫茲涅茲倒U形曲線相似，耕地面積的減少也存在著同樣的動態(tài)平衡點。該平衡點即為曲線的拐點，動態(tài)平衡點達到的時間除受到社會經(jīng)濟發(fā)展變化的影響外，與耕地保護政策的實施力度也有很大關(guān)系。同時耕地保護政策實施力度也會對社會經(jīng)濟發(fā)展產(chǎn)生影響，如何建立它們之間的反饋關(guān)系，合理調(diào)整耕地保護政策，使其能盡快達到耕地總量的動態(tài)平衡點，是今后需要進一步研究的問題。

（

）：

［1］ 翟文俠，黃賢金.我國耕地保護政策運行效果分析［J］.中國土地科學(xué)，2003，17（2）：8-13.

［2］ 張效軍，歐名豪，李景剛.我國耕地保護制度變遷及其績效分析［J］.社會科學(xué)，2007，（8）：13-20.

［3］ 徐夢潔，陳慧中，趙紫玉，等.耕地總量動態(tài)平衡評價模型初探［J］.資源科學(xué)，2009，（12）：2137-2142.

［4］ 黃賢金，濮勵杰，彭補拙.城市土地利用變化及其響應(yīng)：模型構(gòu)建與實證研究［M］.北京：科學(xué)出版社，2008.

［5］ 楊桂山.長江三角洲近50年耕地數(shù)量變化的過程及驅(qū)動機制研究［J］.自然資源學(xué)報，2001，16（2）.

［6］ Pu Li-jie,Yang Gui-shan,Xu Meng-jie.Characteristics and Process of Land Use Changes in the Yangtze River Delta,China［J］.Pedosphere,2001,11（3）:193-198.

［7］ 濮勵杰，黃賢金，彭補拙.耕地總量動態(tài)平衡政策體系研究［J］.南京大學(xué)學(xué)報（自然科學(xué)），2001，37（6）：757-765.

［8］ 丁洪建，吳次芳，梁留科.耕地保護理念的創(chuàng)新研究［J］.中國土地科學(xué)，2002，16（4）：14-19.

［9］ 黃賢金，濮勵杰，周峰，等.長江三角洲地區(qū)耕地總量動態(tài)平衡政策目標實現(xiàn)的可能性分析［J］.自然資源學(xué)報，2002，17（6）：670-676.

［10］姚鑫，楊桂山，萬榮榮.昆山市 22年來耕地流失研究［J］.長江流域資源與環(huán)境，2010，19（9）：1063-1068.