張雅杰，閆小爽，張 豐，肖展春

?

1978—2015年中國多尺度耕地壓力時空差異分析

張雅杰1，閆小爽1，張 豐2※，肖展春1

（1. 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院，武漢 430000； 2. 武漢大學(xué)遙感信息工程學(xué)院，武漢 430000）

為充分認(rèn)識中國耕地壓力的現(xiàn)狀、變化趨勢和空間分布差異，以保障中國的糧食安全，該文利用修正的耕地壓力指數(shù)模型計算了中國31個省級行政區(qū)域1978－2015年的耕地壓力指數(shù)，并在此基礎(chǔ)上從國家、四大經(jīng)濟(jì)區(qū)域、省域3個尺度分析了耕地壓力指數(shù)的時序變化和空間分布差異；該文借鑒力學(xué)原理，引入耕地壓力重心的概念，通過計算耕地壓力重心坐標(biāo)，進(jìn)一步揭示了耕地壓力指數(shù)的空間變動情況。研究發(fā)現(xiàn)，中國耕地壓力指數(shù)表現(xiàn)出明顯的地域差異性，耕地壓力重心隨時間推移逐漸向東南方向移動，具體表現(xiàn)為東部地區(qū)耕地壓力逐漸增加，西部地區(qū)耕地壓力始終較高，而中部和東北地區(qū)耕地壓力基本穩(wěn)定在安全狀態(tài)；通過對各省耕地壓力指數(shù)進(jìn)行系統(tǒng)聚類，可將中國31個省級行政區(qū)分為4種類型，分別為高壓力型、中壓力型、低壓力型和無壓力型。最后該文基于上述研究結(jié)果，對于不同耕地壓力程度和類型的地區(qū)提出針對性的建議來緩解耕地壓力，以更有效地保障各地區(qū)的糧食供給。

土地利用；糧食；模型；修正耕地壓力指數(shù)；糧食安全；耕地壓力重心；聚類分析

0 引 言

耕地資源作為影響糧食產(chǎn)量最基本的條件，其數(shù)量的多寡和質(zhì)量的優(yōu)劣必然對糧食生產(chǎn)產(chǎn)生較為深刻的影響，進(jìn)而將影響國家糧食安全、經(jīng)濟(jì)發(fā)展以及社會進(jìn)步[1]。隨著中國工業(yè)化和城鎮(zhèn)化的快速發(fā)展，耕地數(shù)量逐年減少，耕地質(zhì)量出現(xiàn)嚴(yán)重退化，這導(dǎo)致中國耕地資源正承載著前所未有的糧食生產(chǎn)壓力。此外中國各地區(qū)地形氣候差異較大，耕地后備資源分布不均衡，耕地補(bǔ)充能力在地區(qū)之間差異也較大。因此深入全面的了解中國及各省份的耕地壓力狀況及其時空差異性，將有助于政府合理規(guī)劃分配有限的土地資源，有效統(tǒng)籌耕地跨省占補(bǔ)平衡指標(biāo)，實(shí)事求是地解決耕地問題，保護(hù)中國耕地資源，進(jìn)而保障中國的糧食安全[2]。

當(dāng)前國內(nèi)外學(xué)者對于耕地壓力的研究已取得較為豐碩的成果，國外學(xué)者主要用土地承載力來反映耕地壓力狀況，提出許多有關(guān)土地承載力的估算方法[3-6]。國內(nèi)學(xué)者常采用最小人均耕地面積，人均耕地警戒值、耕地壓力指數(shù)模型等方法來研究某一特定地區(qū)的耕地壓力狀況[7-11]，趙素霞等從省域尺度運(yùn)用耕地壓力指數(shù)模型研究了河南省耕地面積的變化情況和耕地壓力狀況[7]；肖俊磊等從省域尺度研究了四川省2000－2015年間的耕地壓力指數(shù)的變化情況[8]；李倩倩等從區(qū)域尺度計算了關(guān)中地區(qū)的耕地壓力指數(shù)以了解其耕地利用狀況[11]等。通過分析發(fā)現(xiàn)已有的研究其研究尺度較單一且局限于某個區(qū)域、省份、市縣等微觀尺度，缺乏對宏觀多尺度耕地壓力狀況的全面綜合研究，且對較長時間序列上耕地壓力指數(shù)時空差異性的研究也較少。此外，國內(nèi)學(xué)者研究耕地壓力狀況時運(yùn)用較多的方法是蔡運(yùn)龍等提出的耕地壓力指數(shù)模型[12]，該模型能夠綜合農(nóng)業(yè)生產(chǎn)與經(jīng)濟(jì)發(fā)展等要素定量衡量一個地區(qū)耕地資源的緊張程度，但該模型未考慮到不同地區(qū)耕地質(zhì)量及利用狀況的差異性。

為此，本研究綜合考慮耕地質(zhì)量的地區(qū)差異性運(yùn)用修正的耕地壓力指數(shù)模型計算了1978—2015年中國多個空間尺度的耕地壓力指數(shù)，并在此基礎(chǔ)上從國家、四大經(jīng)濟(jì)區(qū)域、省域3個空間尺度來分析耕地壓力指數(shù)的時序變化和空間分布差異，充分認(rèn)識中國耕地壓力現(xiàn)狀、變化趨勢及其空間分布差異，有助于提高中國未來保障糧食安全的水平，為今后國家開展的農(nóng)業(yè)產(chǎn)業(yè)布局、土地利用規(guī)劃等決策提供參考依據(jù)。同時該研究還對國家實(shí)施跨省域補(bǔ)充耕地和城鄉(xiāng)建設(shè)用地增減掛鉤節(jié)余指標(biāo)跨省域調(diào)劑這一決策提供重要借鑒意義，通過對全國各省份耕地產(chǎn)能及耕地壓力現(xiàn)狀的全面綜合研究，有助于國家合理統(tǒng)籌確定不同地區(qū)耕地跨省域占補(bǔ)平衡的指標(biāo)和跨省域補(bǔ)充耕地資金的收取標(biāo)準(zhǔn)，切實(shí)發(fā)揮該政策的實(shí)際效用，實(shí)現(xiàn)全國范圍的耕地總量平衡以及耕地保護(hù)與經(jīng)濟(jì)發(fā)展的雙贏。

1 數(shù)據(jù)來源與研究方法

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

本文所用耕地、人口、經(jīng)濟(jì)數(shù)據(jù)來源于國土資源公報和各級統(tǒng)計年鑒，31個省級行政區(qū)省會地理坐標(biāo)數(shù)據(jù)來源于高德地圖。

由于統(tǒng)計技術(shù)落后、統(tǒng)計標(biāo)準(zhǔn)不統(tǒng)一等原因，導(dǎo)致中國1978—1995年所統(tǒng)計的耕地面積數(shù)據(jù)準(zhǔn)確性較低，因此為提高數(shù)據(jù)的精度，本文以1996年全國土地利用詳查數(shù)據(jù)作為基準(zhǔn)點(diǎn)，通過耕地資源的增減變化對1978—1995年的耕地面積數(shù)據(jù)進(jìn)行重建，得到一套完整連續(xù)且準(zhǔn)確性較高的耕地面積數(shù)據(jù)。此外由于本研究涉及數(shù)據(jù)較多且時間跨度大，導(dǎo)致安徽省和福建省1979年、重慶市1981—1984年的耕地面積數(shù)據(jù)缺失。為了保證時間序列的延續(xù)性和研究結(jié)果的全面性，本文運(yùn)用移動平均的方法將缺失數(shù)據(jù)補(bǔ)充完整[13]。

1.2 研究方法

1.2.1 修正耕地壓力指數(shù)

耕地壓力指數(shù)能夠定量反映一個地區(qū)耕地資源所承受的壓力和能承載的最大人口數(shù)量[12]。而常用的耕地壓力指數(shù)模型僅考慮了耕地數(shù)量與人口的關(guān)系問題，忽略了耕地質(zhì)量和利用效率的異質(zhì)性，因此，本文在耕地壓力指數(shù)模型的基礎(chǔ)上，結(jié)合耕地的區(qū)域差異性，引入質(zhì)量系數(shù)，將耕地壓力指數(shù)拓展為集耕地數(shù)量和質(zhì)量為一體的指標(biāo)[14]，修正后的耕地壓力指數(shù)模型如（1）所示

式中K為原始耕地壓力指數(shù)；min為最小人均耕地面積，hm2；S為實(shí)際人均耕地面積，hm2；為糧食自給率；G為人均糧食需求量，kg；為糧食單產(chǎn)，kg/hm2；為糧食作物種植面積占農(nóng)作物種植面積的比例，以下簡稱糧作比；為復(fù)種指數(shù)；為修正后的耕地壓力指數(shù)；為質(zhì)量系數(shù)；P為各省份糧食單產(chǎn)，kg/hm2；k為各省份復(fù)種指數(shù)；P為全國糧食單產(chǎn)，kg/hm2；k為全國復(fù)種指數(shù)。

當(dāng)=1時，說明區(qū)域內(nèi)現(xiàn)有耕地面積能夠滿足糧食生產(chǎn)的需求，糧食需求與供應(yīng)基本處于平衡狀態(tài)。當(dāng)<1時，實(shí)際人均耕地面積占有量大于滿足生產(chǎn)所需的最小人均耕地面積，說明現(xiàn)有耕地資源的生產(chǎn)力能夠?yàn)楫?dāng)?shù)厝丝谔峁┏渥愕募Z食，耕地基本無壓力。當(dāng)>1時，區(qū)域所擁有的耕地面積已不能滿足保障糧食安全所需要的最少耕地數(shù)量，說明區(qū)域耕地資源較為稀缺且壓力很大，與1相比，值越大，耕地壓力就越大，相應(yīng)的糧食安全保障程度就越低。此外在計算值時，為了科學(xué)地反映區(qū)域耕地壓力情況并能在全國進(jìn)行比較，結(jié)合相關(guān)研究成果，將糧食自給率均設(shè)定為95%[15]，人均糧食需求量按照時間序列依次設(shè)定為1978－1987年為350 kg，1988－1999年為380 kg，2000－2015年為420 kg[16-17]。

1.2.2 耕地壓力重心

為了進(jìn)一步揭示耕地壓力指數(shù)的空間變動情況，本文借鑒力學(xué)重心坐標(biāo)原理[18]，引入耕地壓力重心的概念，計算方法如式（2）所示。

式中M表示各省份的耕地壓力指數(shù)，（X，Y）表示各省會的地理坐標(biāo)，（X，Y）表示第年中國耕地壓力指數(shù)的重心坐標(biāo)。

2 結(jié)果與分析

2.1 中國耕地壓力時序變化分析

在對耕地壓力指數(shù)相關(guān)要素進(jìn)行分析的基礎(chǔ)上，利用修正的耕地壓力指數(shù)模型，計算1978－2015年中國的耕地壓力指數(shù)，結(jié)果如圖1所示。

圖1 中國耕地壓力指數(shù)變化趨勢

從圖1可以看出，1978－2015年中國耕地壓力指數(shù)波動較大，根據(jù)變動趨勢大致可以分為4個階段：1978－1984年為波動降低期；1985－1997年為相對穩(wěn)定期；1998－2003年為快速增長期；2004－2015年為迅速降低期。

1）1978－1984年耕地壓力指數(shù)波動降低期

這一時期，最小人均耕地面積表現(xiàn)出波動降低的態(tài)勢，最小人均耕地面積由1978年的0.146 9 hm2下降到1984年的0.127 9 hm2。耕地壓力指數(shù)除1980和1981年有上升，1978－1984年總體上呈現(xiàn)出減少的趨勢，與最小人均耕地面積變動趨勢基本一致。

2）1985－1997年耕地壓力指數(shù)相對穩(wěn)定期

該時期耕地壓力指數(shù)較為波動，但整體表現(xiàn)為相對穩(wěn)定的狀態(tài)，而最小人均耕地面積表現(xiàn)出緩慢降低的態(tài)勢。耕地壓力指數(shù)除1988年外其他年份均小于1，且均在0.85～1.0范圍內(nèi)。

3）1998－2003耕地壓力指數(shù)快速增長期

這一時期，中國耕地壓力指數(shù)增長迅速，由1998年的0.879 2迅速上升到2003年的1.197 2，較1998年增加0.326 0，耕地狀態(tài)由無壓力轉(zhuǎn)變?yōu)閴毫^明顯。最小人均耕地面積也呈現(xiàn)出較大的增長幅度，由1998年的0.091 1 hm2上升至2003年的0.114 3 hm2，累計增加0.023 2 hm2。

4）2004－2015耕地壓力指數(shù)快速降低期

2004－2015年耕地壓力指數(shù)呈直線下降，耕地壓力指數(shù)由2004年的1.1048下降至2015年的0.8826，下降20.11個百分點(diǎn)。最小人均耕地面積也呈現(xiàn)出下降趨勢，由2004年的0.104 1 hm2縮減至2015年的0.086 7 hm2，下降幅度為16.71%，與耕地壓力指數(shù)變化態(tài)勢大體呈現(xiàn)相似的特點(diǎn)。

2.2 四大經(jīng)濟(jì)區(qū)域耕地壓力時空差異分析

依據(jù)《中國區(qū)域經(jīng)濟(jì)統(tǒng)計年鑒（2010）》，中國可以分為東部、中部、西部和東北4個經(jīng)濟(jì)區(qū)域，這種劃分方式是分析中國人口和社會經(jīng)濟(jì)發(fā)展的常用分區(qū)。由于耕地資源在各經(jīng)濟(jì)區(qū)域內(nèi)的分布極不均衡，且區(qū)域間自然條件與社會經(jīng)濟(jì)發(fā)展的差異也較大，因此本文以經(jīng)濟(jì)區(qū)域?yàn)檠芯繂卧獊矸治龈鱾€區(qū)域耕地壓力的時空變化情況。

本文在計算了31個省級行政區(qū)域耕地壓力指數(shù)的基礎(chǔ)上，發(fā)現(xiàn)各經(jīng)濟(jì)區(qū)域內(nèi)所包含省份的耕地壓力指數(shù)值相差不大，因此取其平均值作為各經(jīng)濟(jì)區(qū)域的耕地壓力指數(shù)，則1978－2015年四大經(jīng)濟(jì)區(qū)域的耕地壓力指數(shù)如圖2所示。

圖2 四大經(jīng)濟(jì)區(qū)域耕地壓力指數(shù)變化趨勢

從圖2中可以看出，耕地壓力表現(xiàn)出明顯的地域差異性，東部地區(qū)耕地壓力呈現(xiàn)逐漸增加的趨勢，而西部地區(qū)耕地壓力始終較大，中部和東北地區(qū)耕地壓力基本穩(wěn)定在安全狀態(tài)。

其中東部地區(qū)耕地壓力逐漸增大的原因在于中國東部地區(qū)經(jīng)濟(jì)發(fā)展水平逐漸提高，工業(yè)化和城鎮(zhèn)化發(fā)展速度逐漸加快，導(dǎo)致大量耕地資源被占用。另外東部地區(qū)因其良好的社會經(jīng)濟(jì)發(fā)展基礎(chǔ)不斷吸引大量人口流入，使得該地區(qū)耕地壓力面臨著嚴(yán)峻的形勢。

西部地區(qū)的耕地壓力始終穩(wěn)定在較高水平，這主要是受到西部地區(qū)自然和社會經(jīng)濟(jì)發(fā)展以及技術(shù)條件的限制，西部地區(qū)有著超過中國耕地資源總量1/3的耕地，幾乎相當(dāng)于東北地區(qū)和東部地區(qū)耕地總量的加和，但這些耕地多以中低產(chǎn)田為主，作物產(chǎn)量很大程度上取決于“天意”[19]。此外西部地區(qū)生態(tài)環(huán)境脆弱，水土流失、土地沙化十分嚴(yán)重，因此耕地生產(chǎn)能力極其低下，耕地壓力較大。

而中部和東北地區(qū)的耕地壓力除個別年份外基本穩(wěn)定在安全狀態(tài)，原因在于該區(qū)域主要位于中國的平原地區(qū)，地勢平坦，耕地資源相對豐富，而且土壤肥沃，光照、溫度、降水等農(nóng)業(yè)生產(chǎn)的條件適宜，使得耕地糧食產(chǎn)量較高，另外其經(jīng)濟(jì)發(fā)展水平不及東部地區(qū)，建設(shè)占用耕地的風(fēng)險相對較低，因此耕地壓力程度相對較小。

2.3 省域尺度耕地壓力時空差異分析

中國31個省級行政區(qū)自然地理環(huán)境差異較大、區(qū)域社會經(jīng)濟(jì)發(fā)展程度極不平衡，致使人口與土地、糧食矛盾地域差異十分明顯，因此研究中國31個省級行政區(qū)的耕地壓力狀況及其時空差異性對于有效保護(hù)耕地和保障糧食安全具有重要的理論和實(shí)踐意義。

2.3.1 省域尺度耕地壓力時序變化分析

由于本文研究時間跨度較大，數(shù)據(jù)較復(fù)雜，因此依據(jù)1978－2015年國家尺度耕地壓力指數(shù)變化趨勢所劃分的4個階段（圖1），這里也從1978－1984、1985－1997、1998－2003、2004－2015年這4個時段來分析31個省級行政區(qū)耕地壓力指數(shù)的變化情況，以每個時段各省份耕地壓力指數(shù)的平均值來反映某一時段耕地所承擔(dān)的壓力，并結(jié)合耕地壓力水平標(biāo)準(zhǔn)的評判法則，劃定耕地壓力指數(shù)平均值在（0，1）為無壓力，[1，1.5）為輕度壓力，[1.5，2）為中度壓力，[2，+∞）為重度壓力。此外，為了進(jìn)一步揭示耕地壓力在4個時段內(nèi)的空間變動情況，本文還借助重心模型分別計算了4個時段的耕地壓力重心坐標(biāo)，并將其分別展會在相應(yīng)時段的耕地壓力程度分布圖中，結(jié)果如表1、圖3所示。

表1 耕地壓力指數(shù)和耕地壓力重心分時段統(tǒng)計Table 1 Sub-period statistics of cultivated land pressure index and cultivated land pressure barycentric coordinates

圖3 31個省級行政區(qū)耕地壓力程度分時段動態(tài)變化

從圖3中可以看出，中國耕地壓力重心隨時間推移逐漸向東南方向移動。從4個時段來看，浙江、福建、廣東、海南等東南沿海地區(qū)的耕地壓力逐步增加，原因在于這些地區(qū)耕地后備資源匱乏，且社會經(jīng)濟(jì)發(fā)展迅速，隨著城市化的快速推進(jìn)，建設(shè)用地占用了大量耕地資源，耕地補(bǔ)充能力不足，同時這些地區(qū)匯聚了大量的人口，導(dǎo)致耕地壓力逐漸增大；而內(nèi)蒙古、新疆、西藏、甘肅、寧夏等位于中國西北地區(qū)的省份的耕地壓力在4個時段內(nèi)均有所下降，原因是這些地區(qū)社會經(jīng)濟(jì)發(fā)展速度遠(yuǎn)低于東南沿海地區(qū)，耕地面積保持相對穩(wěn)定，且隨著科技不斷進(jìn)步、農(nóng)業(yè)投入不斷增加以及開展大規(guī)模的土地開發(fā)整理復(fù)墾項(xiàng)目[20]，使得這些地區(qū)糧食單產(chǎn)和總產(chǎn)不斷提高，耕地壓力有所緩解；此外，山西、陜西、青海等耕地壓力在4個時段內(nèi)均較大，原因是這些地區(qū)山地高原面積多，地形坡度大，不易于開展大規(guī)模的農(nóng)業(yè)生產(chǎn)，且技術(shù)水平落后，復(fù)種指數(shù)較低，糧食產(chǎn)量不高，因此耕地壓力一直較大。而河北、山東、江蘇、河南、湖北、湖南、安徽、江西、四川、黑龍江、吉林、遼寧等位于中部和東北糧食主產(chǎn)區(qū)的省份的耕地壓力指數(shù)在4個時段內(nèi)始終為無壓力或輕度壓力，原因是它們主要位于中國的平原地區(qū)，光照、溫度、降水等農(nóng)業(yè)生產(chǎn)條件適宜，耕地資源相對豐富，生產(chǎn)率水平較高，且糧食產(chǎn)量隨著農(nóng)業(yè)技術(shù)的進(jìn)步而持續(xù)顯著增加，使得這些地區(qū)耕地始終處于輕度壓力甚至無壓力狀態(tài)，并呈現(xiàn)出逐漸減小的趨勢。

2.3.2 省域尺度耕地壓力空間聚類分析

為了在省域尺度上揭示中國耕地壓力的空間差異性，本文采用Ward聚類法[21]對31個省級行政區(qū)1978－2015年的耕地壓力指數(shù)進(jìn)行空間系統(tǒng)聚類分析，根據(jù)聚類結(jié)果，結(jié)合各省份耕地壓力指數(shù)情況，可將31個省級行政區(qū)劃分為4種類型，分別為高壓力型、中壓力型、低壓力型和無壓力型。本文利用ArcGIS對聚類結(jié)果進(jìn)行空間直觀表達(dá)，以清楚反映出各省所屬的耕地壓力類型。

圖4 耕地壓力分類

從圖4中可以看出，京、津、青屬于高壓力型，其耕地壓力指數(shù)遠(yuǎn)大于1，耕地長期承受著巨大壓力；甘、陜、蒙、晉、藏、云屬于中壓力型，這類型耕地壓力指數(shù)相對大些，多分布在1.5～2.5之間，耕地壓力較明顯；粵、貴、瓊、冀、豫、黑、遼、寧、滬、閩、新、浙屬于低壓力型，該類型耕地壓力指數(shù)相對較小，大部分位于1.0～1.5之間，且耕地壓力指數(shù)波動較大；皖、桂、鄂、湘、吉、蘇、贛、魯、川、渝屬于無壓力型，這些地方的耕地壓力指數(shù)較小，多數(shù)地方耕地壓力指數(shù)在1.0以下，耕地基本無壓力。

3 結(jié) 論

1）從國家尺度來看：1978－2015年中國耕地壓力指數(shù)的變化趨勢大致可以劃分為4個階段：1978－1984年的波動降低期；1985－1997年的相對穩(wěn)定期；1998－2003年的快速增長期；2004－2015年的迅速降低期。

2）從經(jīng)濟(jì)區(qū)域尺度來看：耕地壓力指數(shù)表現(xiàn)出明顯的地域差異性，具體表現(xiàn)為東部地區(qū)耕地壓力逐漸增加，西部地區(qū)耕地壓力始終較高，中部和東北地區(qū)耕地壓力基本穩(wěn)定在安全狀態(tài)。

3）從省域尺度來看：對31個省級行政區(qū)1978－2015年的耕地壓力指數(shù)進(jìn)行系統(tǒng)聚類可將其分為4種類型，分別為高壓力型、中壓力型、低壓力型和無壓力型；1978－2015年，中國耕地壓力重心隨時間推移逐漸向東南方向移動。

4 討 論

從研究結(jié)論出發(fā)，本文提出以下幾點(diǎn)建議：

由于各地區(qū)耕地壓力程度不同，且造成耕地壓力的原因不一，因此需因地制宜，針對不同耕地壓力程度和類型的地區(qū)分別采取相應(yīng)的措施來緩解耕地壓力，更有效地保障各地區(qū)的糧食供給。

對于位于東部地區(qū)的北京、上海、浙江、福建、廣東等經(jīng)濟(jì)較發(fā)達(dá)的省份而言，應(yīng)嚴(yán)格實(shí)施國家制定的一系列保護(hù)耕地的政策規(guī)定，如土地用途管制制度，禁止占用基本農(nóng)田，耕地占補(bǔ)平衡制度等[22-24]，其中由于這些省份經(jīng)濟(jì)較發(fā)達(dá)，人口密度較大，對建設(shè)用地的需求較高，導(dǎo)致耕地后備資源已接近枯竭，難以在省域內(nèi)完成耕地占補(bǔ)平衡的要求，因此可考慮向國務(wù)院提出補(bǔ)充耕地國家統(tǒng)籌的申請，經(jīng)國務(wù)院批準(zhǔn)后，可通過購買外省耕地占補(bǔ)平衡指標(biāo)，來保障經(jīng)濟(jì)的快速持續(xù)發(fā)展。此外還應(yīng)按照各省城鎮(zhèn)化發(fā)展的階段適當(dāng)控制人口流入，以免人均耕地面積逐漸減少[25]。

對于位于西部地區(qū)的內(nèi)蒙古、新疆、西藏、青海、陜西、山西等耕地壓力較大的省份而言，應(yīng)積極開展土地整治，加強(qiáng)農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)，加大農(nóng)業(yè)技術(shù)投入，不斷提高耕地質(zhì)量，增加糧食產(chǎn)量，以緩解該地區(qū)的耕地壓力狀況[26]。

對于位于中部和東北糧食主產(chǎn)區(qū)的河北、山東、江蘇、河南、湖北、湖南、安徽、江西、四川、黑龍江、吉林、遼寧等耕地壓力較小甚至無壓力的省份而言，應(yīng)當(dāng)繼續(xù)實(shí)施強(qiáng)農(nóng)利農(nóng)惠農(nóng)政策，確保農(nóng)民進(jìn)行農(nóng)業(yè)生產(chǎn)的熱情不減退，保證糧食播種面積基礎(chǔ)[27]，同時也要繼續(xù)實(shí)施土地整治，完善農(nóng)田基礎(chǔ)設(shè)施，以提高糧食產(chǎn)量，對于有條件的省份，還應(yīng)當(dāng)通過土地流轉(zhuǎn)實(shí)施規(guī)?；?jīng)營，在滿足自身糧食供給的條件下，能夠?yàn)槿奔Z地區(qū)輸送糧食，為國家的糧食安全提供保障[28]。此外，由于這些省份內(nèi)集中連片的耕地后備資源較充足，可以為建設(shè)用地緊張的大城市和缺少耕地的發(fā)達(dá)地區(qū)提供跨省域補(bǔ)充耕地的指標(biāo)，保證全國范圍內(nèi)耕地總量的動態(tài)平衡，并獲得經(jīng)濟(jì)發(fā)展的資金，發(fā)揮經(jīng)濟(jì)發(fā)達(dá)地區(qū)和資源豐富地區(qū)資金資源互補(bǔ)的優(yōu)勢。

本文尚且存在不足之處,對未來的展望如下：1）中國糧食主產(chǎn)區(qū)與非主產(chǎn)區(qū)之間自然條件與社會經(jīng)濟(jì)發(fā)展差異較大，兩大區(qū)域分別具有不同的發(fā)展優(yōu)勢[29-30]，其中糧食主產(chǎn)區(qū)耕地資源豐富，表現(xiàn)出良好的農(nóng)作物種植優(yōu)勢，糧食產(chǎn)量較高，其糧食生產(chǎn)狀況直接影響著全國十三多億人的吃飯問題，而非主產(chǎn)區(qū)二三產(chǎn)業(yè)較發(fā)達(dá)，經(jīng)濟(jì)發(fā)展水平較高，耕地資源相對稀缺，相比糧食主產(chǎn)區(qū)，糧食安全的風(fēng)險更易發(fā)生在非糧食主產(chǎn)區(qū)，由此可見，糧食主產(chǎn)區(qū)與非主產(chǎn)區(qū)之間耕地壓力也表現(xiàn)出明顯的差異性，未來可就此展開更加深入細(xì)致的研究。2）本文所引耕地質(zhì)量系數(shù)是通過地區(qū)耕地與全國耕地的實(shí)際平均產(chǎn)能計算而得，該指標(biāo)反映了耕地的實(shí)際糧食產(chǎn)量，較能體現(xiàn)耕地的質(zhì)量狀況。但近來國家倡導(dǎo)實(shí)施“藏糧于地”戰(zhàn)略[31]，這表明中國不再一味強(qiáng)調(diào)耕地實(shí)際產(chǎn)量的提升，而是提倡通過土地休耕以提高耕地地力，實(shí)現(xiàn)耕地的可持續(xù)利用，促進(jìn)農(nóng)業(yè)的可持續(xù)發(fā)展，保障國家糧食的穩(wěn)定供給和安全，因此未來可對耕地質(zhì)量系數(shù)做進(jìn)一步完善，可以考慮運(yùn)用耕地地力指數(shù)來更好地反映耕地的質(zhì)量狀況。

[1] 魏洪斌，吳克寧，趙華甫，等. 未來10年我國耕地數(shù)量質(zhì)量變化對區(qū)域糧食產(chǎn)能影響研究[J]. 廣東農(nóng)業(yè)科學(xué)，2014，41(19)：213－219，237.

Wei Hongbin, Wu Kening, Zhao Huafu, et al. Research on the impact of the change of quantity of cultivated land on regional grain productivity in the next 10 years in China[J]. Guangdong Agricultural Sciences, 2014, 41(19): 213－219, 237. (in Chinese with English abstract)

[2] 國務(wù)院. 國務(wù)院辦公廳印發(fā)跨省域補(bǔ)充耕地國家統(tǒng)籌管理辦法和城鄉(xiāng)建設(shè)用地增減掛鉤節(jié)余指標(biāo)跨省域調(diào)劑管理辦法的通知[EB/OL]. [2018-03-26]. http://www.gov.cn/ zhengce/content/2018-03/26/content_5277477.htm.

[3] William Rees. Ecological footprints and appropriated carrying capacity: What urban economics leaves out[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1992, 6(2): 121－130.

[4] William Rees.Urban ecological footprints: Why cities cannot be sustainable: And why they are a key to sustainability[J]. Environmental Impact Assessment Review, 1996, 16(4): 223－248.

[5] Vogt W. Road to survival[M]. New York: William Sloan, 1948.

[6] Alan W. The African Husbandman[M]. Edinburgh: Oliver and Boyd, 1965.

[7] 趙素霞，牛海鵬. 基于灰色馬爾科夫模型的河南省耕地壓力狀況研究[J]. 干旱區(qū)資源與環(huán)境，2015，29(8)：46－51.

Zhao Suxia, Niu Haipeng. Study on pressure of cultivated land in Henan Province based on gray markov model[J]. Journal of Arid Land Resources and Environment, 2015, 29(8): 46－51. (in Chinese with English abstract)

[8] 肖俊磊. 基于耕地壓力指數(shù)模型的四川省糧食安全分析與趨勢預(yù)測[D]. 成都：西南交通大學(xué)，2017.

Xiao Junlei. Food Security Analysis and Trend Forecast of Sichuan Province Based on Cultivated Land Pressure Index Model[D]. Chengdu: Southwest Jiaotong University, 2017. (in Chinese with English abstract)

[9] 冉清紅，岳云華，謝德體，等. 中國耕地警戒值的測算與討論[J]. 資源科學(xué)，2007，29(3)：158－164.

Ran Qinghong, Yue Yunhua, Xie Deti, et al. Estimation and discussion on the warning value of cultivated land in China[J]. Resource Science, 2007, 29(3): 158－164. (in Chinese with English abstract)

[10] 郭巍，宋戈. 基于糧食安全的黑龍江省耕地壓力動態(tài)變化定量分析[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報，2009，14(2)：47－51.

Guo Wei, Song Ge. Quantitative analysis of dynamic change of cultivated land pressure in Heilongjiang province based on food security[J]. Journal of China Agricultural University, 2009, 14(2): 47－51. (in Chinese with English abstract)

[11] 李倩倩，陳印軍. 關(guān)中地區(qū)耕地壓力指數(shù)分析及預(yù)測[J].中國農(nóng)學(xué)通報，2011，27(29)：241－245.

Li Qianqian, Chen Yinjun. Analysis and forecast of cultivated land pressure index in Guanzhong region[J]. Chinese Agricultural Science Bulletin, 2011, 27(29): 241－245. (in Chinese with English abstract)

[12] 蔡運(yùn)龍，傅澤強(qiáng)，戴爾阜. 區(qū)域最小人均耕地面積與耕地資源調(diào)控[J]. 地理學(xué)報，2002，12(2)：127－134.

Cai Yunlong, Fu Zeqiang, Dai Erfu. The minimum area per capita of cultivated land and its implication for the optimization of land resource allocation[J]. Acta Geographica Sinica, 2002, 12(2): 127－134. (in Chinese with English abstract)

[13] 李建平，上官周平. 陜西省糧食生產(chǎn)與糧食安全趨勢預(yù)測[J].干旱地區(qū)農(nóng)業(yè)研究，2011，29(4)：245－251.

Li Jianping, Shangguan Zhouping. Forecast of grain production and food security in Shaanxi province[J]. Arid Area Agriculture Research, 2011, 29(4): 245－251. (in Chinese with English abstract)

[14] 劉成武，黃利民. 農(nóng)地邊際化過程中農(nóng)戶土地利用行為變化及其對糧食生產(chǎn)的影響[J]. 地理研究，2015，34(12)：2268－2282.

Liu Chengwu, Huang Limin. The changes of farmers’ behavior in land use and its impacts on the food production during the process of marginalization of arable land[J]. Geographical Research, 2015, 34(12): 2268－2282. (in Chinese with English abstract)

[15] 張利婭. 中國糧食自給率研究[D]. 鄭州：河南農(nóng)業(yè)大學(xué)，2006.

Zhang Liya. Research on the Self-sufficiency Rate of Grain in China[D]. Zhengzhou: Henan Agricultural University, 2006. (in Chinese with English abstract)

[16] 張利國，陳蘇. 中國人均糧食占有量時空演變及驅(qū)動因素[J]. 經(jīng)濟(jì)地理，2015，35(3)：171－177.

Zhang Liguo, Chen Su. The temporal and spatial evolution of the per capita grain possession in China and its driving factors[J]. Economic Geography, 2015, 35(3): 171－177. (in Chinese with English abstract)

[17] 唐華俊，李哲敏. 基于中國居民平衡膳食模式的人均糧食需求量研究[J]. 中國農(nóng)業(yè)科學(xué)，2012，45(11)：2315－2327.

Tang Huajun, Li Zhemin. Per capita food demand based on the balanced dietary patterns of Chinese residents[J]. Chinese Agricultural Sciences, 2012, 45(11): 2315－2327. (in Chinese with English abstract)

[18] 賈科利，馬欣，張俊華. 寧夏耕地壓力與社會經(jīng)濟(jì)發(fā)展耦合時空特征分析[J]. 干旱區(qū)地理，2014，37(4)：812－819.

Jia Keli, Ma Xin, Zhang Junhua. Analysis of spatial and temporal characteristics of cultivated land pressure and social economic development in Ningxia[J]. Arid area geography, 2014, 37(4): 812－819. (in Chinese with English abstract)

[19] 張光宏，崔許鋒. 耕地資源非農(nóng)化驅(qū)動機(jī)制及其區(qū)域差異性[J]. 中國農(nóng)業(yè)科學(xué)，2015，48(8)：1632－1640.

Zhang Guanghong, Cui Xufeng. Non-agricultural driving forces of cultivated land resources and its regional differences[J]. Chinese Agricultural Sciences, 2015, 48(8): 1632－1640. (in Chinese with English abstract)

[20] 陸汝成，黃賢金. 中國省域土地整理復(fù)墾開發(fā)補(bǔ)充耕地的空間效應(yīng)[J]. 水土保持研究，2011，18(6)：172－176.

Lu Rucheng, Huang Xianjin. Spatial effects of supplementary cultivated land reclamation and exploitation in provincial land consolidation in China[J]. Soil and Water Conservation Research, 2011, 18(6): 172－176. (in Chinese with English abstract)

[21] 楊志恒. 基于Ward法的區(qū)域空間聚類分析[J]. 中國人口·資源與環(huán)境，2010，20(增刊1)：382－386.

Yang Zhiheng. Regional spatial clustering analysis based on ward method[J]. China Population·Resources and Environment, 2010, 20(Supp.1): 382－386. (in Chinese with English abstract)

[22] 劉彥隨，喬陸印. 中國新型城鎮(zhèn)化背景下耕地保護(hù)制度與政策創(chuàng)新[J]. 經(jīng)濟(jì)地理，2014，34(4)：1－6.

Liu Yansui, Qiao Luyin. Farmland protection system and policy innovation under the background of new urbanization in China[J]. Economic Geography, 2014, 34(4): 1－6. (in Chinese with English abstract)

[23] 劉彥隨. 中國土地資源研究與學(xué)術(shù)交流新進(jìn)展[J]. 自然資源學(xué)報，2013，28(9)：1479－1487.

Liu Yansui. New development of land resources research and academic exchange in China[J]. Journal of Natural Resources, 2013, 28(9): 1479－1487. (in Chinese with English abstract)

[24] 劉彥隨. 中國土地資源研究進(jìn)展與發(fā)展趨勢[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報，2013，21(1)：127－133.

Liu Yansui. Research progress and development trend of land resources in China[J]. Chinese Journal of Eco-Agriculture, 2013, 21(1): 127－133. (in Chinese with English abstract)

[25] 王琪. 基于糧食安全、經(jīng)濟(jì)發(fā)展和生態(tài)保護(hù)目標(biāo)的成縣耕地保護(hù)研究[D]. 蘭州：甘肅農(nóng)業(yè)大學(xué)，2012.

Wang Qi. Study on Cultivated Land Protection in Cheng County Based on Food Security, Economic Development and Ecological Protection[D]. Lanzhou: Gansu Agricultural University, 2012. (in Chinese with English abstract)

[26] 張帆. 科技投資方向?qū)μ岣呒Z食產(chǎn)量的影響研究[D]. 南京：南京農(nóng)業(yè)大學(xué)，2010.

Zhang Fan. Research on the Impact of Investment in Science and Technology on Improving Grain Output[D]. Nanjing: Nanjing Agricultural College, 2010. (in Chinese with English abstract)

[27] 李本軍.“三農(nóng)”問題與政府行為[D]. 北京：中國農(nóng)業(yè)科學(xué)院，2005.

Li Benjun. Problems Concerning Agriculture, Countryside and Farmers and Governmental Behavior[D]. Beijing: Chinese Academy of Agricultural Sciences, 2005. (in Chinese with English abstract)

[28] 梁宇. 土地流轉(zhuǎn)對糧食安全的影響研究[D]. 荊州：長江大學(xué)，2015.

Liang Yu. Impact of Land Transfer on Food Security[D]. Jingzhou: Yangtze University, 2015. (in Chinese with English abstract)

[29] 崔奇峰，周寧，蔣和平. 糧食主產(chǎn)區(qū)利益補(bǔ)償必要性分析：基于主產(chǎn)區(qū)與非主產(chǎn)區(qū)糧食生產(chǎn)及經(jīng)濟(jì)發(fā)展水平差距的視角[J]. 中國農(nóng)學(xué)通報，2013，29(32)：118－124.

Cui Qifeng, Zhou Ning, Jiang Heping. Necessary analysis of benefits compensation in main grain producing areas:From the perspective of gap between grain production and economic development in main producing areas and non-main producing areas[J]. Chinese Agricultural Science Bulletin, 2013, 29(32): 118－124. (in Chinese with English abstract)

[30] 趙波. 中國糧食主產(chǎn)區(qū)利益補(bǔ)償機(jī)制的構(gòu)建與完善[J]. 中國人口·資源與環(huán)境，2011，21(1)：85－90.

Zhao Bo. Construction and improvement of interest compensation mechanism in main grain producing areas in China[J]. China Population·Resources and Environment, 2011, 21(1): 85－90. (in Chinese with English abstract)

[31] 國務(wù)院. 國務(wù)院關(guān)于建立糧食生產(chǎn)功能區(qū)和重要農(nóng)產(chǎn)品生產(chǎn)保護(hù)區(qū)的指導(dǎo)意見[EB/OL]. [2017-04-10]. http://www.gov.cn/zhengce/content/2017-04/10/content_5184613.htm.

Analysis on temporal-spatial difference of cultivated land pressure at multiple scales in China from 1978 to 2015

Zhang Yajie1, Yan Xiaoshuang1, Zhang Feng2※, Xiao Zhanchun1

(1.,,430000,; 2.,,430000,)

This paper aims to investigate the spatial distribution, temporal change and policy proposal of the cultivated land pressure in China. Considering the relationship between the quantity of cultivated land and population, the model of cultivated land pressure index put forward by Professor Cai Yunlong, can quantitatively measure the tension of cultivated land resources in a region. But the difference of cultivated land quality and utilization in different regions was not considered in this model. Therefore, this paper adopted the quality coefficient to expand the cultivated land pressure index, which can reflect both the quantity and quality of cultivated land. Moreover, most of the existing relating studies focused on shorter time series or smaller research fields, few studies have been conducted with long time series data like the years of 1978-2015 in this paper, and little research about spatial-temporal change of cultivated land pressure has been done at different scales across the whole China. So, this paper used the modified model to calculate the cultivated land pressure index at different levels i.e. the whole country, 4 major economic regions and 31 provincial administrative regions from 1978 to 2015. And based on the above-mentioned research data, by using the methods of comparative analysis, clustering analysis and barycenter analysis, this paper analyzed the temporal changes and spatial distribution differences of cultivated land pressure index at the 3 spatial scales and put forward some proposals for cultivated land protection. The study showed that: 1) The changing trend of the cultivated land pressure index in China can be approximately divided into 4 stages: fluctuant declining period from 1978 to 1984, relatively stable period from 1985 to 1997, rapidly rising period from 1998 to 2003 and rapidly declining period from 2004 to 2015. 2) The regional difference of cultivated land pressure is very obvious. In more detail, the cultivated land pressure in the eastern region of China increases gradually, and in the western region of China, the cultivated land pressure is always high, while it is basically stable at a safe state in the middle and northeastern regions of China. 3) The barycenter of cultivated land pressure in China gradually moves to the southeast over time. 4) Thirty-one provincial administrative regions can be divided into 4 types by using the hierarchical clustering method, namely high-pressure type, medium-pressure type, low-pressure type and non-pressure type. 5) Finally, aiming at different types or degrees of cultivated land pressure, suitable countermeasures should be taken to alleviate the cultivated land pressure effectively so as to guarantee the grain supply in all regions of China, such as purchasing the indicators of cultivated land compensation in other provinces or controlling the population inflow moderately for economically developed provinces of high pressure, carrying out land consolidation and rehabilitation actively or increasing agricultural investment for ecologically fragile provinces of high or medium pressure, improving farmland infrastructure or transporting grain to the food shortage areas or offering the indicators of cultivated land compensation to other provinces with insufficient cultivated land reserve resources for low-pressure or non-pressure provinces located in major grain producing areas, and so on.

land use; grain; models; modified cultivated land pressure index; grain security; barycenter of cultivated land pressure; cluster analysis

張雅杰，閆小爽，張 豐，肖展春. 1978—2015年中國多尺度耕地壓力時空差異分析[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(13)：1－7.doi：10.11975/j.issn.1002-6819.2018.13.001 http://www.tcsae.org

Zhang Yajie, Yan Xiaoshuang, Zhang Feng, Xiao Zhanchun. Analysis on temporal-spatial difference of cultivated land pressure at multiple scales in China from 1978 to 2015[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 1－7. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.13.001 http://www.tcsae.org

2018-01-26

2018-05-08

國家自然科學(xué)基金項(xiàng)目—面向城市擴(kuò)展過程的動態(tài)空間格局指標(biāo)與自組織空間結(jié)構(gòu)分析模型（41571385）；湖北省自然科學(xué)基金項(xiàng)目—綠地生態(tài)網(wǎng)絡(luò)與城鎮(zhèn)擴(kuò)張空間協(xié)同耦合與優(yōu)化模擬（2017CFC847）

張雅杰，吉林人，博士，副教授，研究方向?yàn)橥恋乩门c經(jīng)濟(jì)管理。Email：yjzhang@whu.edu.cn

張 豐，遼寧人，高級實(shí)驗(yàn)師，研究方向?yàn)镚IS與攝影測量。Email：7570056@qq.com

10.11975/j.issn.1002-6819.2018.13.001

F323.21

A

1002-6819(2018)-13-0001-07