高天毅

摘 要：本文通過對基本農田劃定與評價相關文獻進行參考歸納，選取耕地地類、自然等、利用等、經濟等、坡度、面積、與公路距離、與城鎮(zhèn)距離、與農村居民點距離、耕地圖斑形狀指數(shù)、耕地圖斑規(guī)整度11個指標，在神木市永久基本農田成果基礎上，運用灰色關聯(lián)度模型和主成分分析方法分析各指標的影響程度。結果表明，圖斑面積、坡度、耕地類型、與農村居民點距離對永久基本農田劃定相關性最大；使用主成分分析，得到3個綜合因子，能夠為未來基本農田保護和建設提供一定的指導意義。

關鍵詞：永久基本農田；土地評價；灰色關聯(lián)度分析；主成分分析；神木市

中圖分類號：S-3 文獻標識碼：A

DOI：10.19754/j.nyyjs.20190415007

引言

基本農田是按照一定時期人口和社會經濟發(fā)展對農產品的需求，依據(jù)土地利用總體規(guī)劃確定的不得占用的耕地[1]，而永久基本農田是依據(jù)法律實行永久和特殊保護的基本農田。劃定永久基本農田是實施糧食安全戰(zhàn)略、推進農業(yè)結構調整、保障農民利益的重要手段。劃定工作要依據(jù)當?shù)氐慕洕鐣l(fā)展對農產品的需求預測，合理確定基本農田數(shù)量，同時以農用地分等定級成果為參考，確保劃定結果在數(shù)量和質量上具有保證。永久基本農田劃定要求耕地集中連片、生產設施配套、作物高產穩(wěn)產、農業(yè)生態(tài)良好。但在實際基本農田劃定工作中在多地會出現(xiàn)劃定結果以次充好、以劣補優(yōu)、以遠代近等問題，因此，分析永久基本農田劃定成果中不同因素產生的影響狀況是十分有意義的。

國內學者在基本農田劃定與評價等方面做了大量的研究，如董秀茹等人為實現(xiàn)在基本農田劃定中數(shù)量和質量并重的目標，從耕地的自然條件和立地條件2方面出發(fā)，選取耕地等別、坡地、道路通達度、水利設施狀況以及與相鄰土地利用類型的適宜性等7個因素作為評價指標，建立一種土地評價為基礎的劃定方法，借助GIS技術劃定東港市基本農田[2]；錢鳳魁等人將LESA方法與農用地分等成果相結合，建立評價體系，進行耕地質量與立地條件評價分析，在此基礎上對凌源市進行基本農田劃定[3]；章詩捷在引入CLUE-S模型對上饒市廣豐區(qū)未來土地利用變化預測的基礎上，參考農用地分等成果將優(yōu)質耕地優(yōu)先劃入基本農田，并對其連片性進行優(yōu)化調整，與規(guī)劃基本農田進行對比分析[4]；劉慧芳等人以物-場模型為基礎，將基本農田劃定產生的影響理解為基本農田保護區(qū)劃定、基本農田功能要素以及基本農田保護區(qū)影響場3要素有機組合所形成的系統(tǒng)效應，并借鑒空間相互作用的引力模型，構建包含基本農田劃定適宜性、基本農田功能要素協(xié)調性以及距離系數(shù)的基本農田劃定合理性評價模型[5]。但是現(xiàn)有研究大多都是在對耕地或者基本農田建立評價指標體系，對其進行評價分析，得出劃定或優(yōu)化布局后的結果。選取指標方式過于主觀，僅僅依據(jù)基本農田劃定的目標、原則以及國外評價相關經驗的借鑒等確定指標及計算方法。評價的目標和側重不同，評價的結果也會有偏差。

本文通過對基本農田劃定與評價相關文獻進行參考歸納，在以神木市永久基本農田劃定實踐結果的基礎上，選取耕地地類、自然等、利用等、經濟等、坡度、面積、與公路距離、與城鎮(zhèn)距離、與農村居民點距離、耕地圖斑形狀指數(shù)、耕地圖斑規(guī)整度11個普遍使用的指標，利用灰色關聯(lián)度和主成分分析方法對評價指標進行驗證，分析各指標在實際劃定工作中的影響程度，以期為基本農田劃定結果和評價指標體系的優(yōu)化提供理論基礎和方法支撐。

1 研究區(qū)概況

神木市隸屬陜西省，位于陜西省北部，秦晉蒙3?。▍^(qū)）接壤地帶，國土面積達7635km2，是陜西省面積最大的縣?？h域西北部為風沙草灘區(qū)，地勢較為平坦，中部屬黃土丘陵溝壑區(qū)，地形較為破碎，東南部為土石山區(qū)，地面斜度較大，山大溝深，石多土薄，水土流失嚴重。神木市屬中溫帶半干旱大陸性季風氣候，春季干旱多風，夏季短促炎熱，秋季多雨集中，冬季寒冷漫長，多西北風，風沙頻繁，無霜期短，日照豐富。神木市轄19個鎮(zhèn)（辦事處），326個行政村，縣域經濟綜合競爭力居全國百強縣第21位、西北第1位。2014年農業(yè)總產值達21億元，糧食產量22.2萬t，農民人均純收入達到1.4萬元以上，農業(yè)科技貢獻率達到56%。

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

本文研究數(shù)據(jù)分為空間數(shù)據(jù)和統(tǒng)計數(shù)據(jù)。空間數(shù)據(jù)來源包括：行政區(qū)劃數(shù)據(jù)、神木市永久基本農田數(shù)據(jù)、2014年神木市土地利用現(xiàn)狀數(shù)據(jù)、土地利用總體規(guī)劃數(shù)據(jù)、坡度數(shù)據(jù)、農用地分等成果等。統(tǒng)計數(shù)據(jù)主要來源于《2014年神木縣統(tǒng)計年鑒》和《2014年榆林市統(tǒng)計年鑒》。

3 永久基本農田劃定影響因素分析方法

3.1 研究方法

3.1.1 灰色關聯(lián)度模型

灰色關聯(lián)度模型是用于定量分析因子之間相關程度的有效統(tǒng)計方法[6]?；舅枷雽儆趲缀翁幚淼姆懂?，其實質是對反映各因子變化特性的數(shù)據(jù)序列所進行的幾何比較，根據(jù)序列幾何曲線的擬合相似程度來判斷各因子間的相關程度。從變化大小、方向、速度等的相對性方面來描述一個序列中因子間相對變化的情況，如果兩者相對一致，則兩者之間的相關性大，反之亦然。通過對關聯(lián)度的分析，可以看出各個因子對神木市永久基本農田劃定的影響程度，從中發(fā)現(xiàn)主要的影響因素和指標?；疑P聯(lián)度分析具體步驟如下。

3.1.1.1 確定比較序列和參考序列

設評價對象數(shù)量為m個，評價指標數(shù)量為n個。

式中，a為原始變量X的標準化變量（即每個原始變量減去樣本均數(shù)再除以樣本標準差）；cij，i，j=1，…，m為線性組合系數(shù)，被稱為因子載荷量，其大小及前面的正負號直接反映了主成分與相應變量之間關系的密切程度和方向。主成分所反映的是所有樣本的總信息，信息量由Z1到Zm逐漸減少。第i個主成分的貢獻率為λi/m×100；λi為與第i個主成分對應的特征值，可以通過特征方程R-λI=0進行求解，其中R為標準化變量的協(xié)方差矩陣（即相關矩陣），I為與相關矩陣同階的單位矩陣。由此可得，前P個主成分的累積貢獻率是（∑Pi=1λi/m）×100。在應用時，一般取累計貢獻率為70%～85%或以上所對應的前P個主成分即可。

在研究復雜問題時，使用主成分分析方法，往往只需考慮少數(shù)幾個主成分就行，并且不會損失太多信息。這樣做更容易抓住主要矛盾，揭示事物內部變量之間的規(guī)律，同時簡化問題，提高分析效率。

3.2 確定研究單元

考慮到土地利用規(guī)劃管理工作和耕地保護工作共同落實的基礎是永久基本農田保護成果，以及神木市基礎數(shù)據(jù)獲取難易和成果的應用度。本文選取在2014年神木市土地利用現(xiàn)狀數(shù)據(jù)中的耕地圖斑為基礎數(shù)據(jù)，將土地利用規(guī)劃數(shù)據(jù)中已規(guī)劃確定為建設用地的圖斑剔除，并疊加農用地分等成果以及永久基本農田數(shù)據(jù)成果，形成矢量耕地單元79456個。

3.3 指標選取

通過近年來對耕地綜合評價體系進行的大量研究：奉婷等人從耕地自然質量、利用條件、空間形態(tài)和生態(tài)安全4個方面選取分析指標，構建耕地質量綜合評價指標體系[8]；任艷敏等人選擇地形坡度、土壤侵蝕模數(shù)、生態(tài)系統(tǒng)服務價值、地質災害率、連片度等與耕地生態(tài)質量密切相關的指標來構建耕地入選基本農田評價的指標體系，其強調耕地生態(tài)質量的重要性[9]。在前人研究的基礎之上，并結合神木市的實際情況，本文選取11個指標：耕地類型、自然等、利用等、經濟等、坡度、面積、與公路距離、與城鎮(zhèn)距離、與農村居民點距離、耕地斑塊形狀指數(shù)和田塊規(guī)整度，作為研究因子（詳見表1）。

4 結果與分析

4.1 灰色關聯(lián)計算結果分析

通過對研究單元各指標計算，并對各數(shù)據(jù)進行無量綱處理，再通過Matlab軟件編程求解獲得關聯(lián)度結果，如表2所示。

由灰色關聯(lián)度值可知，指標排序為：X6>X5>X1>X9>X8>X3>X7>X2>X4>X11>X10。圖斑面積與永久基本農田劃定的關聯(lián)度最大，關聯(lián)值達到0.706，說明當?shù)卦谶M行永久基本農田劃定工作時，優(yōu)先將集中連片度高，面積較大的耕地劃入，保證農業(yè)集約化生產的需求；坡度對永久基本農田劃定也具有很大的影響程度，關聯(lián)度僅次于圖斑面積，關聯(lián)度值達到0.687，坡度越低越方便進行種植耕作；耕地類型關聯(lián)值為0.608，神木市耕地地類分為水澆地和旱地，水澆地較于旱地耕作條件優(yōu)越，更加有力于農業(yè)耕種作業(yè)，保障糧食生產；與農村居民點距離和城鎮(zhèn)距離的關聯(lián)度分別為0.594和0.582，兩者的相對影響程度基本一致，與農村居民點距離較高，符合耕地耕作實際條件；利用等關聯(lián)度值為0.572，在劃定永久基本農田時，對利用等較為重視，保證劃定的永久基本農田平均利用等高于全市耕地平均利用等，符合工作實際；與公路距離的關聯(lián)度值為0.571，與城鎮(zhèn)距離和居民點距離相比最低，分析其原因，由于神木市中南部處于黃土高原丘陵溝壑區(qū)，在交通的通達度方面不如地形平坦區(qū)域；自然等和經濟等關聯(lián)度值分別為0.546和0.518，由于神木市耕地自然條件一般，對永久基本農田劃定的影響程度稍?。桓貓D斑規(guī)整度和形狀指數(shù)對劃定的影響程度最小，原因在于黃土高原丘陵溝壑區(qū)地形破碎，對圖斑影響主要表現(xiàn)在圖斑邊界更復雜，耕地圖斑規(guī)整度和形狀指數(shù)較低。

4.2 主成分分析

由于所選取指標的計量單位不同，原始數(shù)據(jù)各因子間存在量綱和數(shù)量級不一致的現(xiàn)象，首先需要對原始數(shù)據(jù)進行正向化處理和標準化處理。使用SPSS19.0軟件中的標準差標準化方法，對數(shù)據(jù)進行標準化處理，以此消除量綱差異和數(shù)量級的影響。

計算指標間相關系數(shù)，得出兩兩指標間的相關系數(shù)矩陣，即計算原始指標之間的相關性，是否存在多重共線性現(xiàn)象，以此來判斷數(shù)據(jù)是否適合于因子分析。

對相關系數(shù)矩陣進行KMO（Kaiser-Meyer-Olkin）檢驗統(tǒng)計，結果為0.741，說明指標體系的整體相關程度很高，Bartlett球形檢驗結果為P<0.0001，說明樣本數(shù)量足夠大。根據(jù)上述檢驗結果，認為適合進行主成分分析。

通過分析直接計算得出主成分矩陣特征值和貢獻率、累計貢獻率以及主成分載荷矩陣。利用軟件對數(shù)據(jù)進行主成分分析，最終提取出3個主因子，各主因子Zi的特征值及貢獻率見表3，前3個成分的累積方差貢獻率為85.272%，說明用前3個成分作為新的評價因素可以用來反映神木市永久基本農田劃定影響效果的信息數(shù)據(jù)可靠性在85%以上，由此提取這3個成分作為主成分進行分析，其因子載荷矩陣見表4。

由表4可知，Z1在X1耕地類型、X2自然等、X3利用等、X4經濟等、X5坡度5個因素的載荷量很大，這些指標主要表現(xiàn)農業(yè)耕作生產條件方面對永久基本農田劃定的影響；Z2在X10耕地圖斑形狀指數(shù)、X11耕地圖斑規(guī)整度載荷量很大，主要表現(xiàn)耕地景觀生態(tài)學方面對永久基本農田劃定的影響；Z3在X6圖斑面積、X8城鎮(zhèn)距離、X7與公路距離載荷量很大，說明這3個因素能夠比較準確的反映主成分Z3的信息。提取的該3個主成分基本反映原全部指標的信息，因此可用這3個新變量來代替原來的11個影響因素。

5 結論

通過已有文獻的評價指標體系建立永久基本農田劃定影響因素體系，采用灰色關聯(lián)度模型和主成分分析模型對神木市永久基本農田劃定實際結果進行影響因素分析。主要的成果和結論如下：

通過對11個指標與永久基本農田劃定結果進行灰色關聯(lián)度分析和相關性排序表明，圖斑面積、坡度、耕地類型、與農村居民點距離對永久基本農田劃定相關性最大。

使用主成分分析，得到3個綜合因子，3個主成分的累計方差貢獻率達到85%以上，由于每個主成分都是原始變量的線性組合，主成分之間沒有相關性，能夠通過這3個主成分分析得出更加準確更加清晰的結果。

結果表明，在神木市未來的農業(yè)生產和永久基本農田保護方面，應當因地制宜，注重優(yōu)質耕地的合理利用，減少耕地破碎化，加強集約利用強度和規(guī)模。加大對永久基本農田投入和建設，進一步提高農業(yè)生產效率以及農業(yè)機械化現(xiàn)代化生產水平，不斷推動農業(yè)生產技術進步和創(chuàng)新，逐步優(yōu)化農業(yè)生產方式和保障農民生產生活，提高農民和政府對永久基本農田利用和保護的積極性。

參考文獻

[1] 國務院.基本農田保護條例[Z].1998-12-24.

[2] 董秀茹，尤明英，王秋兵.基于土地評價的基本農田劃定方法[J].農業(yè)工程學報，2011，27（04）：336-339.

[3] 錢鳳魁，王秋兵.基于農用地分等和LESA方法的基本農田劃定[J].水土保持研究，2011，18（02）：251-255.

[4] 章詩捷. 廣豐區(qū)土地利用變化空間模擬及其對基本農田劃定影響分析[D].江西農業(yè)大學，2017.

[5] 劉慧芳，畢如田.基于物-場模型的基本農田劃定合理性評價[J].中國農業(yè)資源與區(qū)劃，2017，38（03）：107-114.

[6] 鄧聚龍.灰色系統(tǒng)基本方法[M].武漢：華中理工大學出版社，1987：74-106.

[7] 王鵬，況福民，鄧育武，等.基于主成分分析的衡陽市土地生態(tài)安全評價[J].經濟地理，2015（1）：124-128.

[8] 奉婷，張鳳榮，李燦，曲衍波，朱鳳凱.基于耕地質量綜合評價的縣域基本農田空間布局[J].農業(yè)工程學報，2014，30（01）：200-210，293.

[9] 任艷敏，孫九林，劉玉，潘瑜春.縣域永久基本農田劃定方法研究[J].農業(yè)機械學報，2017，48（04）：135-141.