楊建宇　張婷婷　陳　正　張　超　趙冬玲　林尚緯

(1.中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院， 北京 100083； 2.國土資源部農(nóng)用地質(zhì)量與監(jiān)控重點實驗室， 北京 100035；3.國土資源部土地整治中心， 北京 100035)

0　引言

近年來，我國經(jīng)濟飛速發(fā)展，城市化與工業(yè)化進程不斷加快[1]，我國耕地每年減少上百萬公頃，耕地資源短缺問題進一步凸顯[2]。土地整治是新增耕地的主要來源方式。土地整治是對低效利用、不合理利用和未利用土地的綜合整治，是對生產(chǎn)建設(shè)破壞和自然災(zāi)害損毀土地的恢復(fù)利用，是土地整理、開發(fā)、復(fù)墾的統(tǒng)稱[3]。我國土地整治工程技術(shù)不斷發(fā)展，土地管理政策要求逐步提高，這也促使了土地整治項目的有序開展，從而增加耕地面積，補充平衡土地供給和社會需求[4-5]。因此，對新增耕地進行時空分布規(guī)律及其來源分析，對于科學(xué)制定土地整治政策和保障京津冀區(qū)域糧食安全、可持續(xù)發(fā)展具有重要指導(dǎo)意義。

現(xiàn)階段，從研究內(nèi)容看，針對新增耕地時空分布規(guī)律的研究較少，缺乏對新增耕地來源的分析。大部分集中于土地整治方法[6-7]、土地整治項目資金效益[8-10]、土地整治質(zhì)量評價[11-12]以及潛力測算[13-15]等方面的研究。從研究對象上看主要以土地整治項目為研究目標(biāo)[16-17]，針對新增耕地時空變化規(guī)律研究較少。從研究尺度來看，宏觀上，多集中于對全國或省域的研究；微觀上，多以縣域為研究單元[18-22]；研究尺度較為單一。受研究方法和數(shù)據(jù)資料等的限制，缺乏對區(qū)域新增耕地圖斑的研究；并且，基于不同研究單元，研究會產(chǎn)生不同的結(jié)論[23]。為綜合分析研究區(qū)域的特征，以新增耕地圖斑為研究對象，采用縣域和區(qū)域兩種研究單元相結(jié)合的方法，可以避免單一尺度造成的不均衡現(xiàn)象，更充分、全面地分析新增耕地現(xiàn)狀，提出更加合理的整治建議。

京津冀是我國重要的糧食主產(chǎn)區(qū)，隨著京津冀協(xié)同發(fā)展上升為國家戰(zhàn)略，京津冀地區(qū)人口不斷增加，建設(shè)用地和生態(tài)用地面積增加，耕地面積逐漸減少，使得土地整治工程成為推動京津冀社會發(fā)展的重要組成部分，但現(xiàn)階段，缺乏對京津冀新增耕地現(xiàn)狀的研究。并且我國提出了建設(shè)雄安新區(qū)的千年大計，雄安新區(qū)的建設(shè)不僅對京津冀地區(qū)的經(jīng)濟發(fā)展有著重大意義，也會對京津冀地區(qū)的土地利用狀況和生態(tài)環(huán)境產(chǎn)生重要影響。因此本文以2013—2015年京津冀土地整治項目新增耕地數(shù)據(jù)為基礎(chǔ)，基于縣域和區(qū)域兩種尺度，采用重心模型、區(qū)位基尼系數(shù)、地理集中度及統(tǒng)計分析等方法分析現(xiàn)階段京津冀地區(qū)新增耕地時空特征及來源特征的變化規(guī)律，為京津冀耕地保護工作提供數(shù)據(jù)參考和相關(guān)建議。

1　研究區(qū)概況與數(shù)據(jù)來源

1.1　研究區(qū)概況

京津冀地區(qū)包含北京市、天津市、河北省，地處華北平原，位于東北亞中國地區(qū)環(huán)渤海心臟地帶，北接內(nèi)蒙古高原，西鄰黃土高原，東臨渤海，與遼寧省、內(nèi)蒙古自治區(qū)、山西省、河南省接壤，地理坐標(biāo)為113°04′～119°53′E、36°01′～42°37′N，總面積2.17×107hm2，見圖1。

圖1　京津冀位置圖Fig.1　Geographical position of Beijing-Tianjin-Hebei region

該區(qū)屬于典型的溫帶半濕潤半干旱大陸性季風(fēng)氣候，四季分明。地勢西北高，東南低。地貌豐富，西北部以山地、丘陵和高原為主，其間分布有盆地和谷地；東南部為廣闊的平原。整個京津冀按照地形因素可以分為后山壩上高原區(qū)、遼吉西蒙東南冀北山地區(qū)、燕山太行山山前平原區(qū)以及冀魯豫低洼平原區(qū)。雄安新區(qū)位于北京、天津、保定腹地，區(qū)位優(yōu)勢明顯、資源環(huán)境承載能力較強。根據(jù)第二次全國土地資源調(diào)查結(jié)果顯示，京津冀耕地面積共7.24×106hm2，占整個京津冀土地面積的1/3，占全國總耕地面積5.35%；其中北京市2.27×105hm2，天津市4.47×105hm2，河北省6.56×106hm2。京津冀地區(qū)耕地主要分布在燕山太行山山前平原區(qū)以及冀魯豫低洼平原區(qū)。灌溉水田集中分布在灤河下游沿海各縣，水澆地和旱地分布廣泛，在整個區(qū)域的平原和盆地地區(qū)集中連片，在山地地區(qū)沿河谷延伸。多年平均水資源總量1.61×1010m3，人均地表水量僅為全國平均水平的10%，水資源缺乏是制約河北省土地資源利用和經(jīng)濟社會發(fā)展的重要因素之一[24-27]。

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

數(shù)據(jù)包括：2013—2015年土地變更調(diào)查新增耕地管理信息核查數(shù)據(jù)，包括新增耕地圖斑面積、地類編碼、變更前地類代碼、來源結(jié)構(gòu)等信息。2013—2015年京津冀分縣行政區(qū)劃數(shù)據(jù)；2013—2015年土地整治項目區(qū)數(shù)據(jù)。

2　數(shù)據(jù)分析與研究方法

為綜合反映新增耕地的時空分布狀況，本文集合了重心模型、區(qū)位基尼系數(shù)和地理集中度3種方法。運用重心模型找到整個研究區(qū)域新增耕地的均衡點，通過重心的偏移衡量不同方位新增耕地貢獻率的差異。運用區(qū)位基尼系數(shù)，衡量各縣級行政單元新增耕地的聚集程度。為避免區(qū)位基尼系數(shù)受行政單元個數(shù)的影響，引入地理集中度的概念。地理集中度可以彌補區(qū)位基尼系數(shù)定性反映的不足，定量地反映縣域新增耕地的空間差異。

2.1　分析流程

利用京津冀行政區(qū)劃數(shù)據(jù)與新增耕地圖斑數(shù)據(jù)疊加，提取京津冀分縣新增耕地圖斑數(shù)據(jù)；再利用重心模型計算新增耕地重心分布，運用區(qū)位基尼系數(shù)和地理集中度，分析京津冀新增耕地整體聚集狀況以及二級地類聚集狀況；按新增耕地屬性信息進行統(tǒng)計，得到京津冀新增耕地地類來源；利用京津冀新增耕地圖斑數(shù)據(jù)與項目區(qū)數(shù)據(jù)疊加，得到京津冀新增耕地來源結(jié)構(gòu)。最后整理分析京津冀新增耕地一級地類特征、二級地類特征及新增耕地來源狀況，如圖2所示。

圖2　方法流程圖Fig.2　Flow chart of method

2.2　重心模型

重心的概念來源于牛頓力學(xué)，原指區(qū)域空間上存在某一點，在該點前后左右各個方向的力量對比保持相對平衡[28]。計算京津冀新增耕地重心時，將新增耕地圖斑的面積作為權(quán)重，圖斑的幾何中心為該圖斑重心，通過對各圖斑重心坐標(biāo)賦予權(quán)重，標(biāo)示出京津冀新增耕地的均值平衡點，與基準(zhǔn)點做對比，衡量京津冀不同區(qū)域新增耕地面積貢獻率的差異?；鶞?zhǔn)點即以京津冀分縣行政區(qū)劃為均質(zhì)區(qū)域，找出其幾何中心，作為衡量的基準(zhǔn)。

據(jù)此，建立計算京津冀土地整治新增耕地面積的重心模型[29]。

(1)

(2)

式中i——年份j——圖斑的編號

n——圖斑總個數(shù)

Sij——第i年、第j個圖斑的面積

Xi——第i年京津冀新增耕地的重心橫坐標(biāo)

Yi——第i年京津冀新增耕地的重心縱坐標(biāo)

Xij——第i年、第j個圖斑的重心橫坐標(biāo)

Yij——第i年、第j個圖斑的重心縱坐標(biāo)

2.3　區(qū)位基尼系數(shù)

區(qū)位基尼系數(shù)是新經(jīng)濟地理學(xué)創(chuàng)始人克魯格曼首次提出的，最初被用于產(chǎn)業(yè)聚集測度的研究[30]，通過不同年份的對比，能夠定性的反映研究目標(biāo)的時空特征和變化趨勢。通過3年間不同類別新增耕地(新增水田、新增水澆地、新增旱地)及其整體的區(qū)位基尼系數(shù)，衡量京津冀土地利用狀況的時空特征。其計算公式為

(3)

式中G——京津冀縣域區(qū)位基尼系數(shù)

k——當(dāng)年京津冀縣域個數(shù)

Sa、Sb——縣域a、b新增耕地的面積

G的取值在0～1之間，取值越小聚集程度越低；反之，聚集程度越高。區(qū)位基尼系數(shù)大于0.5，即視為高度聚集，系數(shù)越大說明聚集程度越高。

2.4　地理集中度

地理集中度是衡量新增耕地空間差異的有效指標(biāo)。選取新增耕地面積最大的前p個縣域計算地理集中度，公式為[31]

(4)

式中Cp——研究區(qū)新增耕地地理集中度

Sq——第q個縣域的新增耕地面積

S——整個研究區(qū)新增耕地面積之和

Cp的取值是0～1,值越大說明地理集中度越高，根據(jù)不同的地類選擇適合的p值，探究其新增耕地的集中區(qū)。

3　結(jié)果與分析

3.1　京津冀新增耕地時空特征

為對京津冀新增耕地進行深入、多元化的分析，并且占補平衡嚴(yán)禁占優(yōu)補劣、占水田補旱地，說明一級地類和二級地類的分布對于土地整治新增耕地等工作的開展有著同等重要的作用，因此分別從區(qū)域和縣域兩個角度分別對2013—2015年京津冀新增耕地一級地類和二級地類分布情況進行分析，以便提出服務(wù)區(qū)域和縣域的不同管理建議。

3.1.1京津冀新增耕地一級地類特征

3.1.1.1京津冀新增耕地一級地類時序特征

匯總2013—2015年3年京津冀新增耕地總面積，新增耕地主要集中在中東部平原區(qū)、西部山區(qū)和北部高原區(qū)3大區(qū)域。東部山前平原區(qū)是高產(chǎn)農(nóng)業(yè)區(qū)，逐漸成為京津冀的重要糧食產(chǎn)區(qū)；北部和西部的高原區(qū)和山區(qū)是京津冀生態(tài)保護的重點區(qū)域，后備資源薄弱，在進行土地整治的同時需注意生態(tài)建設(shè)和保護。

從不同年份分析，2013—2015年新增耕地數(shù)量呈波動態(tài)勢。2013—2015年新增耕地面積分別為1.76×104、1.08×104、1.20×104hm2。借助ArcGIS自然段點法，將京津冀縣域新增耕地面積分為5級，如圖3。2013年面積較大的縣域多集中在中東部偏北的遼吉西蒙東南冀北山地區(qū)和燕山太行山山前平原區(qū)；至2014年，面積較大的新增耕地縣域數(shù)量下降，京津冀南部的冀魯豫低洼平原區(qū)新增耕地面積明顯減少；2015年南部平原地區(qū)縣域新增耕地面積均較高，北部地區(qū)縣域新增耕地面積都維持在較低水平。

究其原因，首先，近幾年耕地的重心開始向質(zhì)量方向轉(zhuǎn)移，對新增耕地數(shù)量的要求有所降低。其次，北部地區(qū)多為高原、山地，受區(qū)位因素影響，可開發(fā)、復(fù)墾的面積有限；從占補平衡政策推出以來，北部地區(qū)可開發(fā)面積逐年遞減，因此，補充新增耕地的任務(wù)逐漸落到南部地區(qū)。

圖3　2013—2015年京津冀新增耕地分布圖Fig.3　Distribution maps of newly increased cultivated land in Beijing-Tianjin-Hebei region from 2013 to 2015

3.1.1.2京津冀新增耕地一級地類空間特征

2013—2015年新增耕地重心不斷向西南方向遷移，詳見圖4，西部和南部新增耕地貢獻率上升。2013年新增耕地重心位于京津冀中東部偏北的地區(qū)，說明京津冀北部地區(qū)新增耕地貢獻率大于南部地區(qū)，東部地區(qū)大于西部地區(qū)；2014年、2015年西部和南部地區(qū)新增耕地貢獻率比重不斷增加，新增耕地重心南移。這與新增耕地面積分布變化趨勢一致。

圖4　2013—2015年新增耕地重心遷移圖Fig.4　Changing map of gravity center of newly increased cultivated land from 2013 to 2015

區(qū)位基尼系數(shù)和地理集中度能夠綜合分析新增耕地的時空變化趨勢。選取p=30，計算京津冀新增耕地地理集中度，結(jié)果如圖5所示。2013—2015年京津冀土地整治新增耕地的區(qū)位基尼系數(shù)和地理集中度平穩(wěn)上升，并沒有太大程度的變化。表明京津冀新增耕地高度聚集，并且聚集程度逐年增加，但整體狀況穩(wěn)定。結(jié)合新增耕地分布圖分析，2013年有新增耕地面積較大的縣域集中在中東部地區(qū)；2014年無新增耕地縣域逐漸集中，有新增耕地縣域逐漸向京津冀周圍擴散；至2015年，無新增耕地縣域主要集中在東北西南軸線上，有新增耕地的縣域呈環(huán)形，集中分布在京津冀周圍；使得3年新增耕地區(qū)位基尼系數(shù)逐漸增加。通過計算地理集中度發(fā)現(xiàn)，整個京津冀地區(qū)新增耕地大量集中于行唐縣、靜?？h、張北縣、威縣、圍場滿族蒙古族自治縣、盧龍縣、平山縣、寶坻區(qū)、密云縣、順平縣等縣域。

圖5　2013—2015年京津冀新增耕地地理集中度與區(qū)位基尼系數(shù)Fig.5　Gini-coefficient and geographic concentration of newly increased cultivated land in Beijing-Tianjin-Hebei region from 2013 to 2015

3.1.2京津冀新增耕地二級地類特征

新增耕地二級地類包括新增水田、新增水澆地、新增旱地3種類型。京津冀是水資源利用限制區(qū)域，水田面積不大，水澆地和旱地占很高的比例。從2013—2015年新增耕地二級地類整體狀況看，除新增水澆地面積在2015年有所增長以外，新增水田和新增旱地均減少。

新增水田、新增水澆地和新增旱地高度聚集且具有明顯的區(qū)域差異。新增水田僅分布在寶坻區(qū)、唐?？h、樂亭縣、灤南縣、寧河縣等幾個縣域，其地理集中度接近于1，3年中并無太大變化，因此不在圖6中進行標(biāo)示。根據(jù)新增水澆地和新增旱地的區(qū)位基尼系數(shù)，分別選取p=30、p=20，計算地理集中度。3年中，新增水澆地的面積最大，也呈現(xiàn)高度聚集狀態(tài)，但是分布比新增水田和新增旱地均勻，主要分布于行唐縣、順平縣、晉州市、圍場滿族蒙古族自治縣、張北縣、密云縣、雄縣、寶坻區(qū)等縣域。新增旱地的面積和聚集程度介于新增水田和新增水澆地之間，主要集中在靜?？h、盧龍縣、撫寧縣、大城縣、豐寧滿族自治縣、威縣、隆化縣等縣域。

圖6　2013—2015年京津冀新增耕地二級地類區(qū)位基尼系數(shù)和地理集中度Fig.6　Gini-coefficient and geographic concentration of three types of newly increased cultivated land in Beijing-Tianjin-Hebei region from 2013 to 2015

平原區(qū)是各類新增耕地的主要分布區(qū)，但不同類別新增耕地分布有較大差異，如圖7所示。新增水田主要集中在燕山太行山山前平原區(qū)，2013年和2015年在冀魯豫低洼平原區(qū)也有分布，但數(shù)量較少。因地形限制，新增水田在遼吉西蒙東南冀北山區(qū)和后山壩上高原區(qū)分布極少。新增水澆地是京津冀地區(qū)新增耕地的主要地類，因平原區(qū)的地勢、交通、水資源等的優(yōu)勢，新增水澆地主要分布在冀魯豫低洼平原區(qū)和燕山太行山山前平原區(qū)；2013年遼吉西蒙東南冀北山區(qū)有一定數(shù)量的新增水澆地分布；但隨著中北部山區(qū)的開發(fā)、整理，新增耕地潛力逐漸降低，2014年和2015年遼吉西蒙東南冀北山區(qū)新增水澆地面積逐漸減少，新增水澆地主要來源于西南部平原區(qū)。除平原區(qū)外，新增旱地在遼吉西蒙東南冀北山地區(qū)有一定比例的分布，但在2014年和2015年所占比例逐漸下降。山區(qū)和高原區(qū)耕地雖然集中，但是自然條件和基礎(chǔ)設(shè)施薄弱，分布零散[32]，不利于大型機械的普及，生產(chǎn)能力較差，因此不作為土地整治的重點發(fā)展區(qū)域。

圖7　2013—2015年京津冀新增耕地二級地類分布情況圖Fig.7　Distribution maps of three types’ newly increased cultivated land in Beijing-Tianjin-Hebei region in 2013—2015

3.2　京津冀新增耕地來源分析

為了解2013—2015年京津冀土地整治新增耕地的成效及其對京津冀生態(tài)環(huán)境的影響，分別從新增耕地來源結(jié)構(gòu)和變更來源兩方面進行分析，提出針對未來土地整治的建議。

3.2.1京津冀新增耕地來源結(jié)構(gòu)

新增耕地共分為6類，各類別及含義如表1所示。

表1　新增耕地分類及含義Tab.1　Different kinds and implication of newly increased cultivated land

表中ZZB為在農(nóng)村土地整治監(jiān)測監(jiān)管系統(tǒng)備案的本年驗收土地整治項目本年變更新增耕地；ZZWB為驗收本年變更新增耕地；ZZBJ為在建未驗收項目本年變更新增耕地，其他均為類別首字母的縮寫。為進一步統(tǒng)一土地整治項目管理與土地變更，2014年起不再允許項目未驗收即開展地類變更，取消了ZZBJ。

2008年以后，取消了土地整理、土地復(fù)墾、土地開發(fā)3類土地整治分類，在開展土地變更調(diào)查新增耕地時，將新增耕地細分為6類。如圖8所示，分析不同來源地類的變化情況，為后續(xù)發(fā)展提供參考。土地整治項目開展以來取得了豐碩的成果，其新增耕地一直是全國補充耕地的主要來源，為保障國家糧食安全發(fā)揮了極為重要的作用。農(nóng)業(yè)結(jié)構(gòu)調(diào)整新增耕地占比逐年上升，說明通過政策引導(dǎo)，京津冀農(nóng)用地結(jié)構(gòu)整理取得顯著成效，農(nóng)用地結(jié)構(gòu)逐漸合理；城鄉(xiāng)建設(shè)用地增減掛鉤新增耕地占比較少，且主要分布在平原區(qū)，說明京津冀平原區(qū)地區(qū)仍具備增減掛鉤補充耕地整治潛力，但在進行補充時不僅要注重數(shù)量的平衡，也要注意耕地質(zhì)量的補充平衡，嚴(yán)防占優(yōu)補劣。工礦廢棄地復(fù)墾新增耕地面積不大，但是在山區(qū)有一定的比例，可以通過整治改善山區(qū)生態(tài)環(huán)境。

圖8　2013—2015年新增耕地來源結(jié)構(gòu)分布狀況Fig.8　Distributions of different areas’ newly increased cultivated land source structure in 2013—2015

圖9　2013—2015年京津冀新增耕地來源Fig.9　Source of newly increased cultivated land in Beijing-Tianjin-Hebei region in 2013—2015

3.2.2京津冀新增耕地變更來源

為了進一步分析新增耕地的變更情況，清楚地掌握近年來新增耕地的來源，利用2013—2015年度京津冀新增耕地數(shù)據(jù)進行地理分析以及統(tǒng)計分析，如圖9所示，得出以下結(jié)論：

(1)3年間新增耕地主要來源于草地、水域及水利設(shè)施用地、其他土地、園地、林地、城鎮(zhèn)村及工礦用地、交通運輸用地及混合地類。其中，草地、水域及水利設(shè)施用地、其他土地是新增耕地最重要的3大來源。特別是草地中的其他草地、水域及水利設(shè)施用地中的內(nèi)陸灘涂、其他土地中的鹽堿地以及園地中的果園，在每年新增耕地來源中都占很大比例。

(2)新增水田、新增水澆地、新增旱地的來源不同。新增旱地和新增水澆地主要來源于草地中的其他草地；新增水田是二級地類中質(zhì)量較高的類型，由于水田對水資源的要求較高，因此新增水田的來源和新增旱地、新增水澆地差異較大，從3年間新增水田的來源情況分析，新增水田主要來源于水域及水利設(shè)施用地中的溝渠。

(3)不同年份，新增耕地來源略有差異。2014年、2015年新增耕地中溝渠占比較低，而2013年溝渠占比達到8.23%，這也是2013年新增水田面積較大的原因之一。從二級地類的角度分析，新增旱地在2013—2015年主要來源地類分別是草地、水域及水利設(shè)施用地和林地。2013年及2015年新增水田主要來源為水域及水利設(shè)施用地，2014年新增水田的主要來源為其他土地。

4　討論與建議

保證數(shù)量是落實藏糧于地、藏糧于技戰(zhàn)略，優(yōu)化用地結(jié)構(gòu)布局、落實節(jié)約優(yōu)先戰(zhàn)略，推動山水林田湖綜合整治、筑牢生態(tài)安全屏障等一系列工作的基礎(chǔ)。綜合上述現(xiàn)狀分析，提出推進京津冀地區(qū)土地整治工作的建議：

(1)從新增耕地分布來看，京津冀平原區(qū)整治開發(fā)潛力較大。根據(jù)第二次全國土地調(diào)查結(jié)果顯示，平原區(qū)自然潛力好，耕地質(zhì)量較高，并且有利于配套基礎(chǔ)設(shè)施與附屬設(shè)施，發(fā)展大型機械化農(nóng)業(yè)。山區(qū)和后山壩上高原區(qū)生態(tài)環(huán)境脆弱、整治難度大，投入大量資金但是難以獲得較高的收益。應(yīng)著力推進平原區(qū)農(nóng)用地整理、農(nóng)村建設(shè)用地及工礦廢棄地復(fù)墾。

(2)大力發(fā)展設(shè)施農(nóng)業(yè)等現(xiàn)代化農(nóng)業(yè)，在水資源利用受限的條件下，必須考慮高效農(nóng)業(yè)，節(jié)約水資源、提高農(nóng)業(yè)附加值。

(3)加大力度實施城鄉(xiāng)建設(shè)用地增減掛鉤，調(diào)整區(qū)域利用形態(tài)，提高土地利用效率，促進城鄉(xiāng)一體化發(fā)展，通過改變城鄉(xiāng)建設(shè)用地布局，進一步拓展空間減少建設(shè)占用新增耕地，從根本上解決耕地保護的難題。

(4)在補充耕地的同時也要注重京津冀地區(qū)生態(tài)環(huán)境的保護。從京津冀一體化發(fā)展的區(qū)域整體功能定位尤其是雄安新區(qū)建設(shè)角度綜合考慮，新增耕地不是本區(qū)域主要需求。在京津冀未來發(fā)展和雄安新區(qū)規(guī)劃中，要進一步注重生態(tài)與環(huán)境安全，合理地實施土地整治項目來補充耕地，宜耕則耕、宜林則林、宜草則草，通過土地整治提高土地利用效率，保證整個京津冀地區(qū)的生態(tài)安全。

不斷推進土地整治工作進程、嚴(yán)守耕地保護紅線才能保證耕地數(shù)量、質(zhì)量、生態(tài)“三位一體”建設(shè)有序開展。

5　結(jié)論

(1)從整體看，2013—2015年京津冀新增耕地主要位于中東部平原區(qū)、西部山區(qū)和北部高原區(qū)；從不同年份看，西部和南部縣域新增耕地貢獻率逐漸增加，新增耕地的聚集程度也逐漸增加。

(2)2013—2015年京津冀地區(qū)新增耕地中占比最大的是新增水澆地，其次為新增旱地。水澆地對水資源要求較高，而京津冀是水資源限制區(qū)，近幾年區(qū)域內(nèi)地下水下降嚴(yán)重；旱地質(zhì)量不高，不能作為發(fā)展重點。因此，要控制農(nóng)業(yè)生產(chǎn)用水區(qū)域，加大節(jié)水力度，提高水資源利用率。

(3)2013—2015年京津冀地區(qū)新增耕地90%來源于土地整治新增耕地，尤其是農(nóng)村土地整治新增耕地。2013年農(nóng)村土地整治新增耕地占新增耕地的比例為87.73%，到2014年，下降為68.58%，2015年為67.20%，近3年來占比雖有所下降，但仍是新增耕地的主要組成部分。

(4)3年中，草地是新增耕地變更來源的重要組成部分，新增旱地和新增水澆地都主要來源于草地，不利于京津冀地區(qū)生態(tài)平衡。從目標(biāo)達成度分析，京津冀地區(qū)環(huán)境承載能力已達上限，要更加注重生態(tài)環(huán)境和人類活動的和諧穩(wěn)定。

1王文剛，龐笑笑，宋玉祥，等.中國建設(shè)用地變化的空間分異特征[J].地域研究與開發(fā)，2012,31(1):110-115.

WANG Wengang, PANG Xiaoxiao, SONG Yuxiang, et al. The spatial differential features of construction land changes in China [J]. Areal Research and Development, 2012, 31(1): 110-115. (in Chinese)

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

ZHAO Xiaoli, ZHANG Zengxiang, WANG Xiao, et al.Analysis of Chinese cultivated land’s spatial-temporal changes and causes in recent 30 years [J]. Transactions of the CASE,2014,30(3):1-11.(in Chinese)

3鹿心社.論中國土地整理的總體方略[J].農(nóng)業(yè)工程學(xué)報,2002,18(1):1-5.

LU Xinshe. General strategy of land consolidation in China [J]. Transactions of the CSAE, 2002, 18(1): 1-5. (in Chinese)

4鄖文聚.土地整治規(guī)劃概論[M].北京：地質(zhì)出版社，2011.

5龍華樓.論土地整治與鄉(xiāng)村空間重構(gòu)[J].地理學(xué)報，2013,68(8):1019-1028.

LONG Hualou. Land consolidation and rural spatial restructuring [J]. Acta Geographica Sinica, 2013, 68(8):1019-1028. (in Chinese)

6GALIULIN R V, BASHKIN V N, GALIULINA R A. Application of peat to land reclamation in a gas production area [J]. Solid Fuel Chemistry, 2013, 47(4): 249-251.

7CAY T, ISCAN F. Fuzzy expert system for land reallocation in land consolidation [J]. Expert System with Applications, 2011, 38(9): 11055-11071.

8胡業(yè)翠，鄭新奇，徐勁原，等.中國土地整治資金效益及其空間分異[J].中國土地科學(xué)，2012,26(2):34-43.

HU Yecui, ZHENG Xinqi, XU Jinyuan, et al. Analysis on the monetary benefit and spatial disparity of land consolidation in China [J]. China Land Sciences, 2012, 26(2): 34-43. (in Chinese)

9谷曉坤，代兵，陳百明.土地整理投資區(qū)域差異、原因及建議[J].中國土地科學(xué)，2007,21(5): 49-53.

GU Xiaokun, DAI Bing, CHEN Baiming. The regional differences, reasons and suggestions of investment in land consolidation [J]. China Land Science, 2007, 21(5): 49-53. (in Chinese)

10金曉斌，丁寧，張志宏，等.中國土地整治資金在省際間分配及土地整治效果[J].農(nóng)業(yè)工程學(xué)報,2012,28(16): 1-9.

JIN Xiaobin, DING Ning, ZHANG Zhihong, et al. Inter-provincial allocation of land consolidation fund and effects of land consolidation in China[J]. Transactions of the CSAE, 2012, 28(16): 1-9. (in Chinese)

11曲衍波，張鳳榮，郭麗娜，等. 農(nóng)村居民點整理后耕地質(zhì)量評價與應(yīng)用[J].農(nóng)業(yè)工程學(xué)報，2012,28(2):226-233.

QU Yanbo，ZHANG Fengrong, GUO Li’na, et al. Estimation of farmland quality after rural residential land consolidation and its application [J]. Transactions of the CSAE, 2012, 28(2): 226-233. (in Chinese)

12高明秀，李占軍，趙庚星.面向土地整理的項目尺度耕地質(zhì)量評價[J].農(nóng)業(yè)工程學(xué)報，2008,24(增刊1)：128-132.

GAO Mingxiu, LI Zhanjun, ZHAO Gengxing. Evaluation of land consolidation oriented cultivated land quality at project scale [J]. Transactions of the CSAE, 2008, 24(Supp.1): 128-132. (in Chinese)

13劉明香,關(guān)欣,徐鄒華,等.土地綜合整治背景下的農(nóng)村居民點整理潛力分析與評價——以花垣縣為例[J].中國農(nóng)學(xué)通報,2013,29(29):103-106.

LIU Mingxiang, GUAN Xin, XU Zouhua, et al. Analysis and evaluation of rural residential reclamation potential under the background of rural land comprehensive reclamation：take Huayuan County as an example [J]. Chinese Agricultural Science Bulletin, 2013, 29(29): 103-106. (in Chinese)

14唐秀美,潘瑜春,劉玉,等.基于耕地系數(shù)和預(yù)評價法的耕地整治潛力測算方法[J].農(nóng)業(yè)工程學(xué)報,2014,30(1)：211-218.

TANG Xiumei, PAN Yuchun, LIU Yu, et al. Calculation method of cultivated land consolidation potential based on cultivated land coefficient and pre-evaluation of farmland classification [J]. Transactions of the CSAE, 2014, 30(1): 211-218. (in Chinese)

15張仕超,魏朝富,李萍. 區(qū)域土地開發(fā)整理新增耕地潛力及其貢獻分析[J]. 農(nóng)業(yè)工程學(xué)報, 2010,26(增刊2): 312-319.

ZHANG Shichao, WEI Chaofu, LI Ping. Analysis of potentialities and contribution rate of newly-increased cultivated land from land development and consolidation [J]. Transactions of the CSAE, 2010,26(Supp.2): 312-319. (in Chinese)

16楊慶媛，馮應(yīng)斌，楊華均，等.土地整理項目空間分異及其與經(jīng)濟發(fā)展的耦合關(guān)系[J].農(nóng)業(yè)工程學(xué)報，2010,26(3): 323-331.

YANG Qingyuan, FENG Yingbin, YANG Huajun, et al. Spatial distribution of land consolidation projects and correlation between them and economic development in Chongqing city[J]. Transactions of the CSAE, 2010, 26(3):323-331. (in Chinese)

17楊緒紅，金曉斌，郭貝貝，等, 2006—2012年中國土地整治項目投資時空分析[J]. 農(nóng)業(yè)工程學(xué)報, 2014, 30(8): 227-235.

YANG Xuhong, JIN Xiaobin, GUO Beibei, et al. Spatio-temporal differentiation of land consolidation investment in China from 2006 to 2012[J]. Transactions of the CSAE, 2014, 30(8): 227-235. (in Chinese)

18楊緒紅,金曉斌，管栩，等. 2006—2012年中國土地整治項目空間特征分析[J]. 資源科學(xué)學(xué)報, 2013, 35(8): 1535-1541.

YANG Xuhong, JIN Xiaobin, GUAN Xu, et al. Spatial characteristics of land consolidation projects in China from 2006 to 2012[J]. Resource Science, 2013, 35(8): 1535-1541.(in Chinese)

19周建,張鳳榮,王秀麗,等. 中國土地整治新增耕地時空變化及其分析[J].農(nóng)業(yè)工程學(xué)報，2014，30(19)：282-289.

ZHOU Jian, ZHANG Fengrong, WANG Xiuli, et al. Spatial-temporal change and analysis of land consolidation’s newly increased cultivated land in China [J]. Transactions of the CSAE, 2014, 30(19): 282-289. (in Chinese)

20魏鳳娟，李江風(fēng)，劉艷中. 湖北縣域土地整治新增耕地的時空特征及其影響因素分析[J].農(nóng)業(yè)工程學(xué)報,2014,30(14):267-276.

WEI Fengjuan, LI Jiangfeng, LIU Yanzhong. Spatial-temporal characteristics and impact factors of newly increased farmland by land consolidation in Hubei province at county level [J]. Transactions of the CSAE, 2014, 30(14): 267-276. (in Chinese)

21杜國明,張繼心,于鳳榮,等.黑龍江省土地整治項目及新增耕地時空格局分析[J].農(nóng)業(yè)現(xiàn)代化研究, 2016, 37(4): 794-801.

DU Guoming, ZHANG Jixin, YU Fengrong, et al. Spatial-temporal pattern of land consolidation programs and newly increased cultivated land in Heilongjiang province [J]. Research of Agricultural Modernization, 2016, 37(4): 794-801. (in Chinese)

22危小建,劉耀林,王娜.湖北省土地整治項目空間分異格局[J].農(nóng)業(yè)工程學(xué)報,2014,30(4)：195-203.

WEI Xiaojian, LIU Yaolin, WANG Na. Spatial disparity pattern of land consolidation projects in Hubei province [J]. Transactions of the CSAE, 2014, 30(4): 195-203. (in Chinese)

23胡業(yè)翠，鄭新奇，徐勁原，等．中國土地整治新增耕地面積的區(qū)域差異[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(2):1-6.

HU Yecui, ZHENG Xinqi, XU Jingyuan, et al. Regional difference for newly increased cultivated land area through land consolidation in China [J]. Transactions of the CSAE, 2012, 28(2): 1-6. (in Chinese)

24魏成林，劉輝．中國耕地質(zhì)量等級調(diào)查與評定(北京卷)[M]．北京:中國大地出版社,2010．

25吳延龍．中國耕地質(zhì)量等級調(diào)查與評定(天津卷)[M]．北京:中國大地出版社，2010．

26蔡博廣，許暤．中國耕地質(zhì)量等級調(diào)查與評定(河北卷)[M]．北京:中國大地出版社，2010．

27李鵬山，呂雅慧，張超，等.基于核密度估計的京津冀地區(qū)耕地破碎化分析[J/OL]．農(nóng)業(yè)機械學(xué)報,2016,47(5): 281-287.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160538&flag=1.DOI:10.6041/j.issn.1000-1298.2016.05.038.

LI Pengshan,Lü Yahui, ZHANG Chao，et al．Analysis of cultivated land fragmentation in Beijing-Tianjin-Hebei region based on kernel density estimation[J/OL]．Transactions of the Chinese Society for Agricultural Machinery,2016,47(5):281-287.(in Chinese)

28陳志剛，王青，黃賢金，等. 長三角城市群重心移動及其驅(qū)動因素研究[J]. 地理科學(xué)，2007,27(4):457-462.

CHEN Zhigang, WANG Qing, HUANG Xianjin, et al. Movement of urban agglomentration gravity and its driving forces in the Changjiang (Yangtze) Delta of China [J]. Scientia Geographica Sinica, 2007, 27(4):457-462. (in Chinese)

29徐建華，岳文澤.近20年來中國人口重心與經(jīng)濟重心的演變及其對比分析[J]. 地理科學(xué)，2001,21(5):385-389.

XU Jianhua, YUE Wenze.Evolvement and comparative analysis of the population center gravity and the economy gravity center in recent twenty years in China [J]. Scientia Geographica Sinica,2001，21(5):385-389. (in Chinese)

30戴平生.區(qū)位基尼系數(shù)的計算、性質(zhì)及其應(yīng)用[J]. 數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究, 2015,32(7):149-161.

DAI Pingsheng. Calculations and characters of the locational gini index and its application[J]. The Journal of Quantitative & Technical Economics,2015,32(7):149-161.(in Chinese)

31張維陽，段學(xué)軍，高金龍，等.中國高新技術(shù)產(chǎn)業(yè)的地理格局與地理集中[J].長江流域資源與環(huán)境，2011,20(7)：830-836.

ZHANG Weiyang, DUAN Xuejun, GAO Jinlong, et al.Geographical pattern and geographical concentration of Chinese new and high tech industries[J]. Resources and Environment in the Yangtze Basin, 2011, 20(7): 830-836.(in Chinese)

32李鵬山，杜振博，張超，等.京津冀地區(qū)耕地質(zhì)量等別空間差異分析[J/OL].農(nóng)業(yè)機械學(xué)報,2017,48(2): 150-157．http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20170220&flag=1&journal_id=jcsam. DOI:10.6041/j.issn.1000-1298.2017.02.020.

LI Pengshan, DU Zhenbo, ZHANG Chao, et al.Space different analysis of cultivated land quality grade in Beijing-Tianjin-Hebei region[J/OL]. Transactions of the Chinese Society for Agricultural Machinery,2017,48(2): 150-157.(in Chinese)