雄安新區(qū)耕地質(zhì)量空間分布特征與建設(shè)占用推薦分區(qū)

高 星，陳 景，劉 蕾，韓 芳，聶承靜，劉丹寧

（1. 河北經(jīng)貿(mào)大學(xué)公共管理學(xué)院，石家莊 050061；2. 河北經(jīng)貿(mào)大學(xué)社會(huì)管理德治與法治協(xié)同創(chuàng)新中心，石家莊 050061；3. 河北經(jīng)貿(mào)大學(xué)土地資源利用與評(píng)價(jià)研究所，石家莊 050061；4. 河北省第三測(cè)繪院，石家莊 050031）

0 引 言

耕地質(zhì)量是國家糧食安全的重要基礎(chǔ)。當(dāng)前中國的耕地保護(hù)已由單純的數(shù)量控制逐步向數(shù)量、質(zhì)量、生態(tài)并重轉(zhuǎn)變。盡管在新型城鎮(zhèn)化進(jìn)程中，已由傳統(tǒng)粗放擴(kuò)張型向集約內(nèi)涵轉(zhuǎn)型，但仍不可避免地要占用一些土地，特別是國家重點(diǎn)戰(zhàn)略的落地。2017年4月1日，河北雄安新區(qū)設(shè)立。雄安新區(qū)近期起步區(qū)和中長期的全面建設(shè)和發(fā)展需要占用大量土地，不可避免將占用部分耕地。同時(shí)，新區(qū)為打造優(yōu)美生態(tài)環(huán)境，實(shí)現(xiàn)藍(lán)綠交織、清新明亮、水城共融的生態(tài)城目標(biāo)，進(jìn)行新型土地綜合整治，具有迫切需求。因此，準(zhǔn)確把握新區(qū)耕地質(zhì)量的現(xiàn)狀空間分布特征和耕地等別差異，對(duì)于合理引導(dǎo)新區(qū)城市發(fā)展和經(jīng)濟(jì)布局，并在土地利用與城市建設(shè)中制定差別化政策，促進(jìn)新區(qū)土地資源合理利用和土地綜合整治具有重要意義。

目前，國內(nèi)學(xué)者基于農(nóng)用地分等成果的研究主要集中在耕地質(zhì)量評(píng)價(jià)[1-3]、耕地現(xiàn)狀分析與提升[4-5]、耕地質(zhì)量動(dòng)態(tài)監(jiān)測(cè)與預(yù)警研究[6-9]、耕地質(zhì)量變化研究[10-12]等方面，涵蓋了全國[13-14]、區(qū)域[15-18]、省級(jí)[19-20]、縣級(jí)[21-22]等尺度，并在城市建設(shè)、土地利用規(guī)劃、土地整治等方面取得了一定成果，促進(jìn)了耕地資源的質(zhì)量保護(hù)與提升[23-27]。由于雄安新區(qū)剛剛成立，尚缺乏針對(duì)全域的耕地質(zhì)量研究成果[28]，因此本文利用雄安 3縣耕地質(zhì)量等級(jí)補(bǔ)充完善數(shù)據(jù)，研究新區(qū)耕地質(zhì)量空間分布特征，為相關(guān)規(guī)劃制定提供依據(jù)。

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

1.1 研究區(qū)概況

雄安新區(qū)地處北京、天津、保定腹地，規(guī)劃范圍涉及河北省雄縣、容城、安新 3縣及周邊部分區(qū)域。新區(qū)地處冀中平原，地跨 115°38′8″～116o19′42″E，38o43′44″～39o10′14″N，北鄰定興縣、高碑店市、廊坊市固安縣，南達(dá)高陽縣、衡水市任丘市，東鄰廊坊市文安縣與霸州市、西接清苑區(qū)、徐水區(qū)。新區(qū)北部、西北部地勢(shì)稍高，南部、東南部地勢(shì)較低，地形開闊，海拔5～20 m，土層深厚。新區(qū)地處中緯度地帶，屬暖溫帶季風(fēng)型大陸性氣候，四季分明，年均氣溫11.7 ℃，年平均降雨量551.5 mm，無霜期185d左右。新區(qū)土地總面積為155 218.9 hm2，2016年土地利用現(xiàn)狀為：耕地88 894.63 hm2，園地2 356.60 hm2，林地3 245.93 hm2，其他農(nóng)用地2 479.41 hm2，城鄉(xiāng)建設(shè)用地25 076.57 hm2，交通水利用地3 725.36 hm2，其他建設(shè)用地1 041.69 hm2，其他土地28 763.31 hm2。新區(qū)共轄29個(gè)鄉(xiāng)鎮(zhèn)，人口合計(jì)112.7萬人，GDP約為215億元。

新區(qū)屬于全國標(biāo)準(zhǔn)耕作制度分區(qū)中黃淮海區(qū)一級(jí)區(qū)，雄縣和安新縣冀魯豫低洼平原區(qū)二級(jí)區(qū)，容城屬于燕山太行山山前平原區(qū)二級(jí)區(qū)，3縣標(biāo)準(zhǔn)耕作制度均為一年二熟，復(fù)種類型為小麥－玉米、小麥/棉花稻，基準(zhǔn)作物為小麥，指定作物為小麥、玉米、棉花3種。3縣共有分等單元1 169 5個(gè)，耕地面積為88 894.63 hm2。

圖1 雄安新區(qū)地理空間位置Fig.1 Regional location of Xiongan New Area

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

本文所采用的數(shù)據(jù)包括2016年雄縣、容城縣、安新縣的耕地質(zhì)量等級(jí)補(bǔ)充完善數(shù)據(jù)與土地變更調(diào)查數(shù)據(jù)庫、1∶10000行政區(qū)劃圖，河北省數(shù)字高程模型（DEM）和2016年雄縣、容城縣、安新縣統(tǒng)計(jì)年鑒。

2 研究方法

本文以雄安新區(qū)耕地為研究對(duì)象，考慮到耕地利用等具有全國可比性，體現(xiàn)耕地質(zhì)量實(shí)際狀況，選擇以耕地利用等作為區(qū)域評(píng)價(jià)基礎(chǔ)，開展耕地質(zhì)量空間分布特征分析。具體方法如下：首先，運(yùn)用耕地質(zhì)量等別面積加權(quán)平均法，計(jì)算新區(qū)耕地質(zhì)量的平均等，分析 3縣各鄉(xiāng)鎮(zhèn)的平均水平；其次，采用耕地質(zhì)量等別均值二分法，計(jì)算新區(qū)所有耕地圖斑的質(zhì)量等別與平均等的離差，繼而采用等間距分級(jí)，得出耕地質(zhì)量等別的空間分異規(guī)律；最后，采用耕地質(zhì)量局部空間自相關(guān)分析方法，得到新區(qū)耕地質(zhì)量空間聚集情況。

2.1 面積加權(quán)平均法

本文采用面積加權(quán)平均法計(jì)算雄安新區(qū)的平均利用等，用以分析新區(qū)的耕地質(zhì)量總體狀況，并為后續(xù)測(cè)算等別離差奠定基礎(chǔ)。計(jì)算公式為

2.2 等別均值二分法

耕地質(zhì)量等別均值二分法即計(jì)算各個(gè)分等單元與平均等別之間的離差，結(jié)果數(shù)值小于 0的代表該分等單元耕地質(zhì)量大于平均質(zhì)量等別，大于 0代表該分等單元耕地質(zhì)量小于平均質(zhì)量。計(jì)算公式為

式中dH 為耕地質(zhì)量等別離差；ih為第i個(gè)分等單元的耕地質(zhì)量等別。

2.3 局部空間自相關(guān)分析法

空間自相關(guān)分析是基于要素位置和要素值來度量空間的相關(guān)性最常用的地統(tǒng)計(jì)學(xué)方法之一，常用的模型有Moran’s I、Geary’s C、Join cont等[29]。以衡量耕地質(zhì)量空間分布特征為目的，為揭示耕地質(zhì)量在空間上的聚集性或離散型特征，本文利用 Geoda軟件，選取 Local Moran’s I為統(tǒng)計(jì)量進(jìn)行新區(qū)耕地質(zhì)量局部空間自相關(guān)分析，以反映相鄰分等單元之間空間變量的相關(guān)程度，便于進(jìn)一步研究耕地質(zhì)量空間分布規(guī)律。其表達(dá)式為

式中 Ii為要素 i的局部 Moran指數(shù)；hi、hj為第 i和第 j個(gè)分等單元耕地質(zhì)量等別； Wij為i和j之間的空間權(quán)重；n為3縣分等單元總個(gè)數(shù)；S2為耕地等別值方差。

Moran指數(shù)取值在-1～1之間，I越大，表明要素的空間結(jié)構(gòu)性越顯著。I等于0時(shí)，表示研究變量不存在空間自相關(guān)性，I>0表示變量在空間上呈現(xiàn)正相關(guān)，存在空間聚集；反之表示研究變量在空間上呈現(xiàn)負(fù)相關(guān)，存在空間孤立；本文利用局部空間自相關(guān)指數(shù)將耕地質(zhì)量空間分布劃分為 5種類型：HH型（高-高）、LL型（低-低）、HL型（高-低）、LH（低-高）型和NS（非顯著）型。正相關(guān)類型包括HH型和LL型，分別代表低等別耕地聚集和高等別耕地聚集（農(nóng)用地分類中，國家等的等別數(shù)值越低代表耕地質(zhì)量最高，反之代表質(zhì)量越低）；負(fù)相關(guān)類型包括HL和LH型，HL型表示低等別耕地被高等別耕地包圍，LH型表示高等別耕地地被低等別耕地包圍；NS型為非顯著型，表示耕地質(zhì)量未呈現(xiàn)出明顯聚集性規(guī)律。

空間權(quán)重是進(jìn)行空間自相關(guān)分析的前提和基礎(chǔ)[30]。基于對(duì)雄安新區(qū)的耕地分等單元鄰接性的分析發(fā)現(xiàn)區(qū)內(nèi)耕地圖斑分布較為均勻，且有相鄰圖斑的地塊占大多數(shù)，判斷采用鄰接原則構(gòu)建空間矩陣較為合理。空間鄰接原則分Rook鄰接原則（公共邊和公共點(diǎn)鄰接）和Queen鄰接原則（公共邊鄰接），對(duì)比基于Rook、Queen鄰接原則構(gòu)建的一階空間鄰接矩陣的鄰接性頻率直方圖，Queen鄰接原則的鄰接權(quán)重頻率更符合正態(tài)分布，因此選擇Queen鄰接原則構(gòu)建空間權(quán)重W矩陣。

3 結(jié)果與分析

3.1 耕地質(zhì)量平均等別分析

根據(jù)全國耕地質(zhì)量等級(jí)調(diào)查與評(píng)定成果，從高到低將全國耕地劃分為優(yōu)等(1～4等)、高等(5～8等)、中等(9～12等)和低等(13～15等)4個(gè)類型[12]。根據(jù)雄縣、容城縣、安新縣的耕地質(zhì)量等級(jí)補(bǔ)充完善數(shù)據(jù)統(tǒng)計(jì)（見表1），雄安新區(qū)耕地以中等地為主，其中高等地20 895.92 hm2，占全區(qū)耕地面積的23.51%，中等地66 638.56 hm2，占全區(qū)耕地面積的74.96%，低等地1 359.49 hm2，占全區(qū)耕地面積的1.53%。

表1 雄安新區(qū)耕地質(zhì)量等別構(gòu)成Table 1 Cultivated land quality of Xiongan New Area

依據(jù)面積加權(quán)平均法計(jì)算雄安新區(qū)耕地質(zhì)量平均等別為9.42等，較河北省耕地質(zhì)量平均利用等為10.82等約高1.4個(gè)等。其中，雄縣、安新縣、容城縣耕地質(zhì)量平均等別為分別為 9.31、9.17、10等。各鄉(xiāng)鎮(zhèn)耕地質(zhì)量平均等別在 8.04~12等之間（其中圈頭鄉(xiāng)無耕地），其中15個(gè)鄉(xiāng)鎮(zhèn)耕地大于新區(qū)平均等別。

3.2 耕地質(zhì)量等別離差分析

采用均值二分法計(jì)算新區(qū)各分等單元與平均等別之間的離差，繼而采用等間距分級(jí)，其結(jié)果見圖2。

圖2 雄安新區(qū)耕地質(zhì)量等別離差Fig.2 Deviation value of cultivated land quality gradation of Xiongan New Area

雄安新區(qū)耕地質(zhì)量等別離差值域范圍在–3.42～4.58之間。質(zhì)量離差計(jì)算結(jié)果顯示，低于平均等別 2等的耕地面積為2 462.24 hm2，占總面積的2.77%；高于平均等別2～4等的耕地面積為5 912.16，占總面積的6.65%；高于低于平均等別4等的耕地面積為545.34 hm2，占總面積的0.61%；近90%的耕地等別離差范圍在±2等之內(nèi)，面積為79974.89 hm2。

總體上而言，質(zhì)量離散程度較低耕地主要分布在新區(qū)北部和東部，離散程度較高的耕地主要分布在新區(qū)中部和南部。

3.3 局部空間自相關(guān)分析

雄安新區(qū)耕地質(zhì)量局部空間自相關(guān)指數(shù)為0.756 8，表明在空間分布上具有較強(qiáng)的正相關(guān)性，如表 3所示。屬于正相關(guān)類型HH與LL型耕地面積分別為19 458.06 hm2與22 457.11 hm2，共占耕地總面積的52.45%；屬于負(fù)相關(guān)類型的LH和HL型的耕地面積分別為1 468.23 hm2和645.13 hm2，兩者共占2.44%；NS型耕地面積為44 866.1 hm2，占耕地總面積的45.1%。因此正相關(guān)類型占主導(dǎo)地位，耕地質(zhì)量局部空間分布情況見圖3。

從空間分布上來看，LL型耕地大多分布在新區(qū)東北部和西南部，局部區(qū)域之間空間正相關(guān)性強(qiáng)烈，整體空間連通性和農(nóng)業(yè)基礎(chǔ)設(shè)施條件較好、耕作限制因素少、耕地質(zhì)量普遍較高，大體上分布在新蓋房溢洪區(qū)東北部及唐河新道兩岸，集中了區(qū)內(nèi)大部分高于平均等的耕地。HH型耕地分布較為分散，整體空間連通性和灌溉條件較差、土壤鹽漬化程度較高、耕地質(zhì)量普遍較低，空間覆蓋范圍相對(duì)來說不如LL型集中，主要分布在新區(qū)的西北部、東部及南部的部分區(qū)域。LH型及HL型耕地分布無集中趨勢(shì)，通常與HH型和HL型相鄰，分別主要分布在新區(qū)中部和西南部。綜上，通過采用局部空間自相關(guān)分析法對(duì)新區(qū)耕地質(zhì)量空間分布進(jìn)行分析，得出的結(jié)果與實(shí)際情況相符，表明該方法適用。

表3 雄安新區(qū)耕地質(zhì)量等別局部空間自相關(guān)類型Table 3 Local indices of spatial autocorrelation types of cultivated land qualitygradation in Xiongan New Area

圖3 雄安新區(qū)耕地質(zhì)量局部空間自相關(guān)聚集圖Fig.3 Local indices of spatial autocorrelation of cultivated land quality gradation in Xiongan New Area

4 建設(shè)占用耕地區(qū)域劃定與合理性評(píng)價(jià)

4.1 建設(shè)占用耕地區(qū)域劃定規(guī)則

在耕地質(zhì)量單要素約束下，根據(jù)耕地質(zhì)量空間分布聚集特點(diǎn)制定新區(qū)建設(shè)占用耕地劃定規(guī)則。HH型耕地是低質(zhì)量等別地塊高度聚集區(qū)域，質(zhì)量離差較小，且內(nèi)部耕地等別存在較強(qiáng)的空間正相關(guān)性，耕地質(zhì)量障礙因素較多，在新區(qū)建設(shè)不可避免占用耕地的前提下，占用地等別的耕地對(duì)糧食產(chǎn)能的影響最小，因此作為建設(shè)占用優(yōu)先區(qū)域；LL型耕地是高等別耕地在空間上聚集分布，質(zhì)量離差同樣較小，內(nèi)部存在較強(qiáng)的空間正相關(guān)性，耕地質(zhì)量障礙因素較少，且耕地產(chǎn)能要明顯高于其他地塊，未來在新區(qū)進(jìn)行土地綜合整治時(shí)的資源投入量會(huì)明顯小于其他耕地，應(yīng)長期保持耕地利用的現(xiàn)狀；NS型耕地局部區(qū)域內(nèi)耕地質(zhì)量高低不等，質(zhì)量離差較大，且耕地質(zhì)量影響因素差別較大，內(nèi)部之間空間無相關(guān)性，未來進(jìn)行土地綜合整治投入資源較大，不利于集中連片土地綜合整治，因此建議作為建設(shè)占用后備區(qū)；HL與LH型耕地局部區(qū)域內(nèi)耕地質(zhì)量高等別地塊被低等別地塊包圍或低等別被高等別耕地，內(nèi)部離差較大，且存在一定的空間負(fù)相關(guān)性，但面積較小，因此就近并入相近區(qū)域；白洋淀水域在新區(qū)建設(shè)中的地位至關(guān)重要，作為新區(qū)的生態(tài)景觀核心，未來擴(kuò)展的可能性較大，將水域周圍海拔低于湖面的耕地以及現(xiàn)狀類型為內(nèi)陸灘涂、沼澤地、坑塘水面作為水系調(diào)整區(qū)。

4.2 建設(shè)占用耕地推薦分區(qū)結(jié)果

根據(jù)建設(shè)占用耕地區(qū)域劃定規(guī)則，在考慮村域行政界限盡量完整的情況下，將新區(qū)建設(shè)用地占用耕地分為建設(shè)優(yōu)先占用區(qū)、建設(shè)后備占用區(qū)和耕地保持區(qū)和水系調(diào)整區(qū)（圖4）。

圖4 雄安新區(qū)建設(shè)用地占用耕地推薦分區(qū)Fig.4 Recommending zoning of cultivated land occupied by construction land in Xiongan New Area

1）建設(shè)優(yōu)先占用區(qū)

該區(qū)耕地質(zhì)量等別在10～14等之間，局部空間自相關(guān)類型為HH型，區(qū)域面積為40 981.03 hm2，占全區(qū)面積的26.40%。從空間分布上來看，該區(qū)域主要分布新區(qū)北部，區(qū)位條件較好，距離容城鎮(zhèn)和雄州鎮(zhèn)較近。另外，新區(qū)東部和西部也有HH型區(qū)域。受灌溉條件、土壤鹽漬化等因素影響，該地區(qū)耕地質(zhì)量較低，可作為優(yōu)先占用區(qū)。

2）建設(shè)后備占用區(qū)

后備占用區(qū)耕地質(zhì)量等別在 9～11等之間，局部空間自相關(guān)類型為NS型，區(qū)域面積為57 304.07 hm2，占全區(qū)面積的36.92%。這類耕地平均質(zhì)量等別較低，但高于建設(shè)優(yōu)先占用區(qū)耕地。本區(qū)的NS型耕地大多分布在HH型周圍，且與之形成了良好的連片性，因此可以作為建設(shè)后備占用的耕地，與優(yōu)先占用區(qū)形成集中連片式開發(fā)建設(shè)，可為建設(shè)用地選址提供更加廣闊的選擇余地。

3）耕地保持區(qū)

耕地保持區(qū)耕地質(zhì)量等別在6～9等之間，局部空間自相關(guān)類型為LL型，區(qū)域面積為39 008.29 hm2，占全區(qū)面積的25.13%，主要分布在新區(qū)東北部和南部。該區(qū)域內(nèi)的耕地綜合質(zhì)量最高，且該區(qū)域整體耕地連片性較高，基礎(chǔ)設(shè)施完善，耕地產(chǎn)能高，土地整治難度較小，未來可以投入較小的人力、物力和財(cái)力，即可建設(shè)成為高標(biāo)準(zhǔn)農(nóng)田、城市生態(tài)農(nóng)業(yè)和旅游休閑地，因此應(yīng)劃入耕地保持區(qū)，嚴(yán)格禁止非農(nóng)建設(shè)。

4）水系調(diào)控區(qū)

該類型耕地面積比例較小，主要分布于在白洋淀水域周邊，區(qū)域面積為17 925.5 hm2，占全區(qū)面積的11.55%。

4.3 建設(shè)占用耕地推薦分區(qū)合理性評(píng)價(jià)

從保護(hù)優(yōu)質(zhì)耕地視角分析，優(yōu)先占用區(qū)主要位于雄安新區(qū)中北部，集中連片性較好，面積能夠滿足100 km2起步區(qū)建設(shè)的需要，且南接白洋淀水域，有利于打造水城交融的新城。其余的優(yōu)先占用區(qū)分散在新區(qū)東部和中西部，可滿足城市多中心的規(guī)劃發(fā)展需求。后備占用區(qū)主要分布在優(yōu)先占用區(qū)周圍，在未來拓展中，利于通過利用后備占用區(qū)向東西 2個(gè)方向擴(kuò)展，選擇余地較大。耕地保持區(qū)主要位于新區(qū)南部和東北部，區(qū)域位置遠(yuǎn)離新區(qū)中心，其中南部耕地保持區(qū)地處白洋淀上游，開發(fā)建設(shè)對(duì)白洋淀產(chǎn)生生態(tài)風(fēng)險(xiǎn)較大。因此，利用耕地質(zhì)量等別空間分布特征劃分建設(shè)占用耕地推薦分區(qū)符合新區(qū)建設(shè)需要，具有較強(qiáng)的合理性。

5 結(jié)論

本文主要基于雄縣、安新縣和容城縣耕地質(zhì)量等級(jí)補(bǔ)充完善數(shù)據(jù)，采用面積加權(quán)平均法、耕地質(zhì)量二分法、局部空間自相關(guān)分析法對(duì)新區(qū)耕地質(zhì)量進(jìn)行了測(cè)算，分析了耕地質(zhì)量空間分布特征，劃定了建設(shè)占用推薦分區(qū)。結(jié)論如下：

1）采用耕地質(zhì)量等別面積加權(quán)平均法計(jì)算了雄安新區(qū)的國家平均利用等，經(jīng)計(jì)算得到新區(qū)平均耕地質(zhì)量等別為9.42等，雄縣、安新縣、容城縣耕地質(zhì)量平均等別為分別為9.31、9.17、10等，均高于河北省10.82等的平均等別。

2）采用耕地質(zhì)量等別均值二分法測(cè)算了雄安新區(qū)耕地與平均等別之間的離差，分析可知新區(qū)耕地質(zhì)量離差值域范圍在–3.42～4.58之間，近90%的耕地等別離差范圍在±2等之間。

3）局部空間自相關(guān)分析結(jié)果表明，雄安新區(qū)耕地質(zhì)量空間聚集程度，分析可知新區(qū)LL型高等地明顯聚集分布于新區(qū)中部和南部，HH型耕地面積較小且分布較為分散，負(fù)相關(guān)類型（HL和LH型）耕地則無明顯集中區(qū)域，零星分布于縣域范圍內(nèi)。

4）并針對(duì)雄安新區(qū)建設(shè)可能占用耕地的事實(shí)，根據(jù)耕地局部空間自相關(guān)類型，制定了建設(shè)占用耕地區(qū)域劃定規(guī)則，劃定了優(yōu)先占用區(qū)、后備占用區(qū)和耕地保持區(qū)和水系調(diào)整區(qū)等新區(qū)建設(shè)占用耕地推薦分區(qū)，并對(duì)推薦分區(qū)合理性進(jìn)行了評(píng)價(jià)。

雄安新區(qū)未來將是生態(tài)宜居城市的典范，一方面新區(qū)的建設(shè)不可避免的占用耕地，耕地的生產(chǎn)功能將逐漸弱化，另一方面耕地的景觀文化功能、生態(tài)服務(wù)功能等將更加凸顯，建議未來從耕地多功能角度對(duì)雄安新區(qū)的耕地空間分布特征及功能增強(qiáng)進(jìn)行深入研究。

[1] 付國珍，擺萬奇. 耕地質(zhì)量評(píng)價(jià)研究進(jìn)展及發(fā)展趨勢(shì)[J].資源科學(xué)， 2015， 37（2）: 226－236.Fu Guozhen, Bai Wanqi. Advances and prospects of evaluating cultivated land quality[J]. Resources Science,2015, 37（2）: 226－236.(in Chinese with English abstract)

[2] 秦元偉，趙庚星，董超，等. 鄉(xiāng)鎮(zhèn)級(jí)耕地質(zhì)量綜合評(píng)價(jià)及其時(shí)空演變分析[J]. 自然資源學(xué)報(bào)，2010，25(3)：454－464.Qin Yuanwei, Zhao Gengxing, Dong Chao, et al.Comprehensive evaluation and spatio-temporal variation analysis of cultivated land quality at town scale[J]. Journal of Natural Resources, 2010, 25(3): 454－464. (in Chinese with English abstract)

[3] 方琳娜，宋金平. 基于SPOT多光譜影像的耕地質(zhì)量評(píng)價(jià)—以山東省即墨市為例[J]. 地理科學(xué)進(jìn)展，2008，27(5)：71－78.Fang Linna, Song Jinping. Cultivated land quality assessmentbased on SPOT multispectral remote sensing Image: A case study in Jimo city of Shandong province[J]. Progress In Geography,2008, 27(5): 71－78. (in Chinese with English abstract)

[4] 陳印軍，肖碧林，方琳娜，等. 中國耕地質(zhì)量狀況分析[J].中國農(nóng)業(yè)科學(xué)，2011，44(17)：3557－3564.Chen Yinjun, Xiao Bilin, Fang Linna, et al. The quality analysis of cultivated land in China[J]. ScientiaAgriculturaSinica, 2011,44(17): 3557－3564. (in Chinese with English abstract)

[5] 徐明崗，盧昌艾，張文菊，等. 我國耕地質(zhì)量狀況與提升對(duì)策[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2016，7(37)：8－14.Xu Minggang, Lu Changai, Zhang Wenju, et al. Situation of The quality of arable land in China and improvement strategy[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2016, 7(37), 8－14. (in Chinese with English abstract)

[6] 王倩，尚月敏，馮銳，等. 基于變異函數(shù)的耕地質(zhì)量等別監(jiān)測(cè)點(diǎn)布設(shè)分析：以四川省中江縣和北京市大興區(qū)為例[J].中國土地科學(xué)，2012，26(8)：80－86.Wang Qian, Shang Yuemin, Feng Rui, et al. Study on location sampling for monitoring the quality of arable land based on variation functions: Case studies in Zhongjiang County and Daxing District[J]. China Land Sciences, 2012, 26(8): 80－86. (in Chinese with English abstract)

[7] 張蚌蚌，孔祥斌，鄖文聚，等. 我國耕地質(zhì)量與監(jiān)控研究綜述[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2015，20(2)：216－222.Zhang Bangbang, Kong Xiangbin, Yun Wenju, et al. A Review on quality and monitoring of arable land in China[J].Journal of China Agricultural University, 2015, 20(2): 216－222. (in Chinese with English abstract)

[8] 任艷，陳蘭康，尹秋月，等. 耕地質(zhì)量監(jiān)測(cè)漸變分布范圍確定方法研究[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2017，38(1)：38-44.Ren Yan, Chen Lankang, Yin Qiuyue, et al. Gradient Distribution of cultivated land quality monitoring[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2017, 38(1):38－44. (in Chinese with English abstract)

[9] 馬建輝. 耕地質(zhì)量等別監(jiān)測(cè)成果應(yīng)用方法研究[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2016，37(10)：55－66.Ma Jianhui. Methods of monitoring cultivated land quality based on variation function and cultivated land gradation type[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2016, 37(10): 55－66. (in Chinese with English abstract)

[10] 陳朝，呂昌河. 基于綜合指數(shù)的湖北省耕地質(zhì)量變化分析[J]. 自然資源學(xué)報(bào)，2010，25(12)：2018－2029.Chen Zhao, Lü Changhe. An analysis on the variation of farmland quality using an integrated index approach in Hubei Province[J]. Journal of Natural Resources, 2010, 25(12):2018－2029. (in Chinese with English abstract)

[11] 姜廣輝，趙婷婷，段增強(qiáng)，等. 北京山區(qū)耕地質(zhì)量變化及未來趨勢(shì)模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(10)：304－311.Jiang Guanghui, Zhao Tingting, Duan Zengqiang, et al.Cultivated land quality change and its future trend modeling in Beijing mountainous area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(10): 304－311. (in Chinese with English abstract)

[12] 肖麗群，陳偉，吳群，等. 未來10a長江三角洲地區(qū)耕地?cái)?shù)量變化對(duì)區(qū)域糧食產(chǎn)能的影響：基于耕地質(zhì)量等別的視角[J]. 自然資源學(xué)報(bào)，2012，27(4)：565－576.Xiao Liqun, Chen Wei, Wu Qun, et al. Impact of cultivated land quantity change on grain productive capacity of Yangtze River delta in the next 10 years: From the perspective of cultivated land quality grade[J]. Journal of Natural Resources,2012, 27(4): 565－576. (in Chinese with English abstract)

[13] 王洪波，程鋒，張中帆，等. 中國耕地等別分異特性及其對(duì)耕地保護(hù)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2011，27(11)：1－8.Wang Hongbo, Cheng Feng, Zhang Zhongfan, et al. Differential characteristics of cultivated land grade and its effect on cultivated land protection in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(11): 1－8. (in Chinese with English abstract)

[14] 程鋒，王洪波，鄖文聚. 中國耕地質(zhì)量等級(jí)調(diào)查與評(píng)定[J].中國土地科學(xué)，2014，28(2)：75－84.Cheng Feng, Wang Hongbo, Yun Wenju. Study on investigation and assessment of cultivated land quality grade in China[J].China Land Sciences, 2014, 28(2): 75－84. (in Chinese with English abstract)

[15] 孔祥斌，張青璞. 中國西部區(qū)耕地等別空間分布特征[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28(22)：1－7.Kong Xiangbin, Zhang Qingpu. Spatial distribution characteristics of arable lang grade in Western China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(22): 1－7. (in Chinese with English abstract)

[16] 楊瑞珍，陳印軍. 東北地區(qū)耕地質(zhì)量狀況及變化態(tài)勢(shì)分析[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2014，6(35)：19－24.Yang Ruizhen, Chen Yinjun. Analysis of quality of cultivated and land change trend in northeast area[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2014, 6(35):19－24. (in Chinese with English abstract)

[17] 李鵬山，杜振博，張超. 京津冀地區(qū)耕地質(zhì)量等別空間差異分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2017，48(2)：150－157.Li Pengshan, Du Zhenbo, Zhang Chao. Space difference analysis of cultivated land quality grade in Beijing-Tianjin-Hebei region[J]. Agricultural Engineering, 2017, 48(2): 150－157.(in Chinese with English abstract)

[18] 魏洪斌，吳克寧，趙華甫，等. 中國中部糧食主產(chǎn)區(qū)耕地等別空間分布特征[J]. 資源科學(xué)，2015，37(8)：1552－1560.Wei Hongbin, Wu Kening, Zhao Huafu, et al. Spatial distribution characteristics of cultivated land quality gradation in the main grain production area of central China[J]. Resources Science, 2015, 37(8): 1552－1560. (in Chinese with English abstract)

[19] 張紅富，周生路，吳紹華，等. 省域尺度耕地產(chǎn)能空間分異規(guī)律及其影響因子[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(8)：308－314.Zhang Hongfu, Zhou Shenglu, Wu Shaohua, et al. Provincial scale spatial variation of cultivated land production capacity and its impact factors[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2010, 26(8): 308－314. (in Chinese with English abstract)

[20] 李武艷，朱從謀，王華，等. 浙江省耕地質(zhì)量多尺度空間自相關(guān)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(23)：239－245.Li Wuyan, Zhu Congmou, Wang Hua, et al. Multi-scale spatial autocorrelation analysis of cultivated land quality in Zhejiang province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(23): 239－245. (in Chinese with English abstract)

[21] 王立為，安萍莉，潘志華，等. 半干旱區(qū)氣候變化背景下近 20年內(nèi)蒙古武川縣耕地質(zhì)量變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(11)：224－231.Wang Liwei, An Pingli, Pan Zhihua, et al. Cultivated land quality change of Wuchuan county in Inner Mongolia under background of climate change in semi-arid regions during recent 20 years[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013,29(11): 224－231. (in Chinese with English abstract)

[22] 奉婷，張鳳榮，李燦，等. 基于耕地質(zhì)量綜合評(píng)價(jià)的縣域基本農(nóng)田空間布局[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30(1)：200－210.Feng Ting, Zhang Fengrong, Li Can, et al. Spatial distribution of prime farmland based on cultivated land quality comprehensive evaluation at county scale[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(1): 200－210. (in Chinese with English abstract)

[23] 龔媛，李飛雪，王麗妍，等. 耕地空間優(yōu)化配置研究：以常州市新北區(qū)為例[J]. 水土保持研究，2016，23(4)：199－205.Gong Yuan, Li Feixue, Wang Liyan, et al. Studies on spatial optimization allocation of cultivated land: Taking Xinbei district,Changzhou City as the example[J]. Research of Soil and Water Conservation, 2016, 23(4): 199－205. (in Chinese with English abstract)

[24] 高星，吳克寧，鄖文聚，等. 縣域耕地后備資源與規(guī)劃期內(nèi)數(shù)量質(zhì)量并重的占補(bǔ)平衡分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(12)：213－219.Gao Xing, Wu Kening, Yun Wenju, et al. Analysis on county based reserved resource for cultivated land and quality-quantity requisition-compensation balance in planning period[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2015, 31(12): 213－219. (in Chinese with English abstract)

[25] 高星，吳克寧，陳學(xué)砧，等. 土地整治項(xiàng)目提升耕地質(zhì)量可實(shí)現(xiàn)潛力測(cè)算[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(16)：233－240.Gao Xing, Wu Kening, Chen Xuezhen, et al. Feasible potential of cultivated land quality promoted by land consolidation project[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(16): 233－240. (in Chinese with English abstract)

[26] 楊建宇，杜貞容，杜振博，等. 基于耕地質(zhì)量評(píng)價(jià)和局部空間自相關(guān)的高標(biāo)準(zhǔn)農(nóng)田劃定[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2017，48(6)：109－115.Yang Jianyu, Du Zhenrong, Du Zhenbo, et al.Well-facilitied capital farmland assignment based onland quality evaluation and LISA[J].Agricultural Engineering, 2017, 48(6): 109－115. (in Chinese with English abstract)

[27] 宋文，吳克寧，張敏，等. 基于村域耕地質(zhì)量均勻度的高標(biāo)準(zhǔn)農(nóng)田建設(shè)時(shí)序分區(qū)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(9)：250－259.Song Wen, Wu Kening, Zhang Min, et al. High standard farmland construction time sequence division based on cultivated land quality uniformity in administrative village scale[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(9): 250－259. (in Chinese with English abstract)

[28] 姜魯光，呂佩憶，封志明，等. 雄安新區(qū)土地利用空間特征及起步區(qū)方案比選研究[J]. 資源科學(xué)，2017，39(6)：991－998.Jiang Luguang, Lü Peiyi, Feng Zhiming, et al. Land use patterns of the Xiongan New Area and comparison among potential choices of start zone[J]. Resources Science, 2017,39(6): 991－998. (in Chinese with English abstract)

[29] 陳彥光. 基于 Moran 統(tǒng)計(jì)量的空間自相關(guān)理論發(fā)展和方法改進(jìn)[J]. 地理研究，2009，28(6)：1449–1463.Chen Yanguang. Reconstructing the mathematical process of spatial autocorrelation based on Moran’s statistics[J].Geographical Research, 2009, 28(6): 1449–1463. (in Chinese with English abstract)

[30] Anselin L. Spatial Econometrics, Methods and Models[M].Boston: Kluwer Academic, 1988: 16－21.