張轉(zhuǎn)+衛(wèi)新東+王篩妮

摘要：為掌握渭北旱塬區(qū)耕地質(zhì)量等別動(dòng)態(tài)變化，加強(qiáng)耕地質(zhì)量、數(shù)量、生態(tài)三位一體管護(hù)，以陜西省黃陵縣為例，對耕地質(zhì)量等別監(jiān)測樣點(diǎn)布控進(jìn)行研究。采用多因素限制條件法、網(wǎng)格法和分層抽樣法，對比提出縣域耕地等別監(jiān)測區(qū)樣點(diǎn)布設(shè)和精度檢驗(yàn)最佳方法。結(jié)果表明，縣域耕地質(zhì)量監(jiān)測控制區(qū)劃分依據(jù)自然質(zhì)量控制區(qū)、土地利用等值區(qū)、土地經(jīng)濟(jì)等值區(qū)，通過ArcGIS空間疊加功能確定；綜合確定監(jiān)測樣點(diǎn)布控的最佳方法為多因素限制條件法；黃陵縣耕地質(zhì)量等別監(jiān)測布設(shè)45個(gè)監(jiān)測點(diǎn)，經(jīng)統(tǒng)計(jì)學(xué)檢驗(yàn)，監(jiān)測樣點(diǎn)代表性程度較好，能夠全面監(jiān)測和反映耕地質(zhì)量等別動(dòng)態(tài)變化。

關(guān)鍵詞：土地利用；監(jiān)測樣點(diǎn)；多因素限制條件法；監(jiān)測控制區(qū)

中圖分類號：F301.21 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2018）02-0048-07

DOI：10.14088/j.cnki.issn0439-8114.2018.02.013

Abstract： To master dynamic changes of cultivated land quality and enhance the quality，quantity and ecological protection of cultivated land，chosing Huangling county in Shaanxi province as an example，monitoring samples of cultivated land quality were studied. Multi-factor constraints method，grid method and stratified sampling method were used to do comparative analysis seeking for the best way to establish monitoring samples and check the acceptance at the county level. The results showed that monitoring control areas were divided according to natural quality control area，land use equivalent area and land economic equivalent area by using ArcGIS spatial overlay analysis. Multi-factor constraints method was selected as the best way for this research. The Huangling county was divided into 45 monitoring samples. After statistical examination，it shown a good representativeness of monitoring samples，which were able to monitor the dynamic changes of cultivated land comprehensively.

Key words： land use； monitoring sample； multi-factor constraints method； monitoring control area

中國人口眾多，土地資源相對缺乏，耕地資源尤為稀缺，已不足世界人均耕地的40%[1]。在人均耕地資源緊缺條件下，有限的耕地資源迫使農(nóng)戶土地利用目標(biāo)多樣，使國家政策調(diào)控的空間有限[2-4]。為有效緩解耕地?cái)?shù)量不足、質(zhì)量不高、生態(tài)環(huán)境不斷退化帶來的危機(jī)，中國的耕地保護(hù)工作已從數(shù)量管理、數(shù)量質(zhì)量并重管理，向數(shù)量、質(zhì)量和生態(tài)管護(hù)三位一體轉(zhuǎn)變[5]。同時(shí)，實(shí)施耕地質(zhì)量監(jiān)測，掌握耕地質(zhì)量動(dòng)態(tài)變化，是進(jìn)行耕地質(zhì)量宏觀精準(zhǔn)管理的基礎(chǔ)[6]。以陜西省為例，在耕地?cái)?shù)量管理方面，構(gòu)建了比較完整的土地?cái)?shù)量調(diào)查制度和技術(shù)支撐體系；但在耕地質(zhì)量管理方面，目前還處于探索階段。結(jié)合現(xiàn)有土地資源調(diào)查情況，基于地統(tǒng)計(jì)學(xué)中遙感技術(shù)布設(shè)耕地等別監(jiān)測點(diǎn)，對耕地質(zhì)量進(jìn)行動(dòng)態(tài)監(jiān)測，及時(shí)了解耕地質(zhì)量變化、掌握耕地總體質(zhì)量情況，對糧食安全、經(jīng)濟(jì)可持續(xù)發(fā)展具有積極的影響意義[7-9]。

國內(nèi)外圍繞耕地質(zhì)量及其監(jiān)測方法開展了大量研究，美國從統(tǒng)計(jì)學(xué)角度出發(fā)，建立國家資源清單計(jì)劃并進(jìn)行長期定位監(jiān)測，完成了國家資源監(jiān)測點(diǎn)的布控[10]；歐盟[11]和加拿大[12]采用“網(wǎng)格法”和“監(jiān)測區(qū)-監(jiān)測樣點(diǎn)”的方法，對全國的資源完成了監(jiān)測布點(diǎn)工作。近幾年，國內(nèi)很多研究學(xué)者在耕地質(zhì)量監(jiān)測方面進(jìn)行了深入的研究。張蚌蚌等[13]、吳克寧等[14]詳細(xì)介紹了標(biāo)準(zhǔn)樣地在耕地質(zhì)量等級監(jiān)測方面的應(yīng)用原則及方法。楊建宇等[15]提出了一種基于空間模擬退火算法的布樣及優(yōu)化方法，以Kriging預(yù)測標(biāo)準(zhǔn)差作為優(yōu)化準(zhǔn)則，初步實(shí)現(xiàn)了樣點(diǎn)布設(shè)的自動(dòng)化。殷守強(qiáng)等[16]采用定量的方法，對耕地質(zhì)量監(jiān)測中6種不同抽樣方法進(jìn)行對比，提高了樣點(diǎn)的全局代表性。張玉臻等[17]在省級標(biāo)準(zhǔn)樣地基礎(chǔ)上進(jìn)行耕地監(jiān)測樣地?cái)?shù)量和質(zhì)量兩個(gè)方面的優(yōu)化，形成省級耕地質(zhì)量監(jiān)測布點(diǎn)的技術(shù)方法。王倩等[18]和孫亞彬等[19]分別基于變異函數(shù)及潛力指數(shù)組合等方法探索縣域耕地質(zhì)量監(jiān)測樣地布控，以實(shí)現(xiàn)較小的監(jiān)測點(diǎn)反映縣域耕地質(zhì)量分布和變化趨勢。綜上所述，對于耕地質(zhì)量監(jiān)測的研究多集中于縣域尺度內(nèi)的布點(diǎn)，缺少從省和國家尺度上探討耕地質(zhì)量監(jiān)測樣地的布控；同時(shí)選取監(jiān)測點(diǎn)的方法較為單一，今后可以通過建立數(shù)據(jù)模型進(jìn)行監(jiān)測點(diǎn)的布控，使得監(jiān)測點(diǎn)的選取更加科學(xué)。

選取黃陵縣作為渭北旱塬區(qū)的代表縣，開展耕地質(zhì)量等別監(jiān)測評價(jià)試點(diǎn)研究，即以2013年度黃陵縣耕地質(zhì)量等別年度更新評價(jià)工作為基礎(chǔ)，在耕地質(zhì)量等別漸變區(qū)域，開展耕地質(zhì)量等別監(jiān)測評價(jià)試點(diǎn)研究，探索縣域內(nèi)耕地質(zhì)量等別監(jiān)測區(qū)劃分、監(jiān)測樣點(diǎn)布設(shè)方法以及監(jiān)測樣點(diǎn)布設(shè)合理性檢驗(yàn)分析，形成縣級耕地質(zhì)量等別監(jiān)測評價(jià)試點(diǎn)成果。結(jié)合試點(diǎn)情況，對黃陵縣乃至渭北旱塬區(qū)耕地質(zhì)量等別監(jiān)測區(qū)劃分和監(jiān)測樣點(diǎn)布設(shè)方法進(jìn)行總結(jié)，為全面開展耕地質(zhì)量等別監(jiān)測提供經(jīng)驗(yàn)參考。endprint

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

1.1 研究區(qū)概況

渭北旱塬區(qū)主要包括陜西的咸陽市、銅川市部分塬區(qū)、延安南部及渭南北部等所轄區(qū)14縣3區(qū)，總面積約3萬km2，該區(qū)地形和氣候獨(dú)特，地形特點(diǎn)為黃土塬、梁、溝壑縱橫交錯(cuò)，黃土覆蓋深厚，溝谷深80～200 m，年降雨量550～600 mm。研究區(qū)選在延安市黃陵縣，該縣地貌類型復(fù)雜多樣，為典型的黃土旱塬地貌。縣境內(nèi)山、川、溝、坡、臺(tái)并存，地貌以梁峁、溝坡為主，土地利用類型分布具有明顯的地域特征，屬中溫帶大陸型季風(fēng)氣候，總體特征為干旱少雨，溫度變化大，降雨變率大，年際和季節(jié)分配不勻，無霜期長，農(nóng)作物可一年一熟或二年三熟。主要糧食作物為小麥和玉米，全年糧食播種面積10 431 hm2，其中小麥播種面積3 193 hm2，玉米播種面積4 021 hm2，其他糧食播種面積為3 217 hm2。2013年，在國家農(nóng)機(jī)具補(bǔ)貼政策優(yōu)惠下，全縣擁有各類農(nóng)業(yè)機(jī)械總動(dòng)力達(dá)256 219 kW，農(nóng)田有效灌溉面積1 447 hm2，化肥施用量（折純量）10 805 t，較上年增長10.80%?？h境內(nèi)交通發(fā)達(dá)，公路、鐵路運(yùn)輸暢通，網(wǎng)絡(luò)密布。

黃陵縣共轄6鎮(zhèn)、3社區(qū)、1鄉(xiāng)以及1個(gè)街道辦事處，共計(jì)215個(gè)行政村。根據(jù)2013年土地利用更新調(diào)查，縣耕地面積為14 668.13 hm2，占土地總面積的6.43%。其中水田62.90 hm2，占耕地面積的0.43%；水澆地276.46 hm2，占耕地面積的1.88%；旱地14 328.77 hm2，占耕地面積的97.69%。由此可知，耕地中以旱地為主，僅有少量的水田和水澆地。2011-2013年突變耕地總面積為1 546.53 hm2（新增耕地面積1 395.80 hm2，質(zhì)量建設(shè)耕地面積150.73 hm2），且全部為旱地。扣除突變耕地面積，漸變區(qū)域耕地面積為13 121.60 hm2，作為本次黃陵縣耕地質(zhì)量等別監(jiān)測對象。黃陵縣不同坡度級別的耕地面積分布情況如表1所示。

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

①黃陵縣國土、統(tǒng)計(jì)、農(nóng)業(yè)、水利、林業(yè)等部門收集糧食產(chǎn)量、糧食價(jià)格、種子、化肥、農(nóng)藥等數(shù)據(jù)；②耕地質(zhì)量年度監(jiān)測的基礎(chǔ)資料，包括黃陵縣2013年耕地質(zhì)量等別年度更新評價(jià)成果、2013年土地利用變更數(shù)據(jù)庫、2011-2013年土地利用變更調(diào)查地類面積匯總表、近3年土地開發(fā)整理項(xiàng)目規(guī)劃設(shè)計(jì)圖等資料；③2014年《黃陵縣統(tǒng)計(jì)年鑒》。

2 耕地質(zhì)量監(jiān)測樣點(diǎn)布控方法

在進(jìn)行監(jiān)測樣點(diǎn)布控前，首先劃定監(jiān)測控制區(qū)。采用ArcGIS空間疊置方法，對影響耕地等級變化的自然狀況、利用水平和收益水平進(jìn)行疊加分析確定監(jiān)測區(qū)。通過對比分析常用的多因素限制條件法、方格網(wǎng)法和分層抽樣法等監(jiān)測點(diǎn)布控方法，在監(jiān)測控制區(qū)內(nèi)布設(shè)監(jiān)測樣點(diǎn)，進(jìn)行監(jiān)測樣點(diǎn)數(shù)量確定。

2.1 監(jiān)測控制區(qū)劃分方法

監(jiān)測控制區(qū)是根據(jù)影響耕地質(zhì)量的各分等因素及其組合，依綜合主導(dǎo)因素原則和地域分異原則劃分區(qū)域，區(qū)域內(nèi)耕地質(zhì)量相對均勻，區(qū)域間有明顯差別。黃陵縣耕地質(zhì)量監(jiān)測控制區(qū)劃分方法如下。

1）影響耕地質(zhì)量的自然因素有地貌、氣候、水文、土壤等，它們相互聯(lián)系、相互作用、相互制約，以不同方式，從不同側(cè)面，按不同程度，獨(dú)立或綜合地影響著土地資源的綜合特征。因此根據(jù)地貌、氣候、土壤等自然因素，考慮環(huán)境污染狀況，劃分耕地自然質(zhì)量分區(qū)。

2）影響耕地利用水平的因素有種植業(yè)結(jié)構(gòu)、主要作物產(chǎn)量水平、農(nóng)田基礎(chǔ)設(shè)施及土地整治等工程，在耕地分等中，耕地利用水平由種植作物的產(chǎn)量直接反映，土地利用系數(shù)是耕地分等單元作物實(shí)際產(chǎn)量與耕作制度區(qū)域內(nèi)最大糧食產(chǎn)量的比值，是耕地利用水平的直接量化指標(biāo)。因此根據(jù)土地利用系數(shù)劃分利用水平分區(qū)。

3）影響耕地收益水平的因素有種植業(yè)投入-產(chǎn)出和效益等，土地經(jīng)濟(jì)系數(shù)是耕地分等單元的“產(chǎn)量-成本”指數(shù)與耕作制度區(qū)域內(nèi)最大“產(chǎn)量-成本”指數(shù)的比值，是耕地收益水平的直接量化指標(biāo)。因此根據(jù)土地經(jīng)濟(jì)系數(shù)劃分經(jīng)濟(jì)水平分區(qū)。

4）綜合耕地自然質(zhì)量控制區(qū)、土地利用等值區(qū)和土地經(jīng)濟(jì)等值區(qū)形成縣域耕地質(zhì)量監(jiān)測控制區(qū)。

2.2 監(jiān)測樣點(diǎn)數(shù)量確定方法

綜合國內(nèi)耕地監(jiān)測的研究進(jìn)展，胡曉濤等[20]采用變異函數(shù)模型進(jìn)行監(jiān)測樣點(diǎn)布局，結(jié)合耕地等別類型進(jìn)行監(jiān)測點(diǎn)調(diào)整和優(yōu)化；蔡鷺斌等[21]結(jié)合中國耕地質(zhì)量整體狀況和借鑒外國經(jīng)驗(yàn)的基礎(chǔ)上，提出利用“分區(qū)組合法”對耕地質(zhì)量監(jiān)測樣點(diǎn)進(jìn)行布設(shè)；楊建宇等[22]提出基于空間平衡法的縣域耕地質(zhì)量監(jiān)測布樣方法。結(jié)合生產(chǎn)工作實(shí)際情況，黃陵縣耕地質(zhì)量等別監(jiān)測點(diǎn)布設(shè)主要采用3種方法，分別是多因素限制條件法、方格網(wǎng)法和分層抽樣法。

2.2.1 多因素限制條件法 黃陵縣地勢西北高、東南低，加之境內(nèi)地貌類型復(fù)雜，所以在監(jiān)測樣點(diǎn)的布設(shè)中主要考慮耕地等別、耕地地類、監(jiān)測控制區(qū)、土地利用規(guī)劃狀況、監(jiān)測便利程度（對于山區(qū)而言主要包括坡度和交通便利度情況）等多種因素[23]。在實(shí)際情況中，由于高等別耕地多為水田，而低等別耕地多為旱地，因此對不同等別耕地監(jiān)測樣點(diǎn)個(gè)數(shù)進(jìn)行合理選取，實(shí)現(xiàn)對耕地類型選取的合理分配。

在監(jiān)測點(diǎn)的選取過程中，首先將監(jiān)測樣點(diǎn)按耕地等別類型落位在集中連片的耕地上，然后根據(jù)土地利用規(guī)劃、坡度、監(jiān)測便利度等因素對監(jiān)測樣點(diǎn)數(shù)量和位置進(jìn)行微調(diào)，保證每個(gè)監(jiān)測樣區(qū)里每一個(gè)等別類型上均至少布設(shè)一個(gè)監(jiān)測樣點(diǎn)，監(jiān)測樣點(diǎn)都應(yīng)處于土地利用總體規(guī)劃確定的基本農(nóng)田上。

由于黃陵縣處于渭北旱塬區(qū)，地貌類型復(fù)雜，絕大多數(shù)耕地分布在道路、溝壑附近，坡度和交通通達(dá)度對于監(jiān)測點(diǎn)的選取有著很大的制約性，為了保證監(jiān)測點(diǎn)的實(shí)際可行性，在監(jiān)測點(diǎn)選取過程中設(shè)置了坡度小于15°，距離道路溝壑1 km范圍以內(nèi)的限制性條件。選點(diǎn)過程借助ArcGIS中的按屬性選擇功能，首先提取坡度小于15°的耕地圖斑，然后對道路溝壑圖層進(jìn)行緩沖區(qū)分析，緩沖距為1 km，再將上述兩個(gè)圖層與黃陵縣耕地圖斑層進(jìn)行疊加分析，選取符合條件的圖斑，設(shè)立監(jiān)測樣點(diǎn)。endprint

2.2.2 方格網(wǎng)法 方格網(wǎng)法選取監(jiān)測樣點(diǎn)是根據(jù)研究區(qū)的實(shí)際情況，按照一定的區(qū)域面積，在研究區(qū)的土地利用現(xiàn)狀圖上繪制方格網(wǎng)，然后將落在耕地圖斑上的格網(wǎng)交點(diǎn)或中心點(diǎn)作為監(jiān)測樣點(diǎn)[24]。采用方格網(wǎng)法劃分耕地質(zhì)量等別類型主要是根據(jù)黃陵縣轄區(qū)面積、地形地貌、耕地分布等情況。具體方法為基于ArcGIS軟件平臺(tái)，通過與千米格網(wǎng)的疊加分析，采用4 km×4 km網(wǎng)格大小，劃分190個(gè)網(wǎng)格單元，最后將格網(wǎng)中心點(diǎn)與耕地相重疊的耕地圖斑作為監(jiān)測樣點(diǎn)。

2.2.3 分層抽樣法 針對縣域耕地包含的各等級，樣點(diǎn)容量y與抽樣精度x的函數(shù)關(guān)系類似分層抽樣模型[25]，如下：

式中，h為耕地等級；nh為縣域內(nèi)h等級耕地的圖斑數(shù)；n為全縣耕地圖斑總數(shù)；Wh為h等級耕地權(quán)重；Sh為h等級耕地等指數(shù)的真實(shí)標(biāo)準(zhǔn)差；Ch為調(diào)查h等級耕地的單樣本費(fèi)用。

將等級成果中的數(shù)據(jù)帶入式（1），可得因變量樣本容量與自變量抽樣精度的表達(dá)式y(tǒng)=f（x），對其求導(dǎo)[26]。

式中，導(dǎo)數(shù)表征樣本容量y隨抽樣精度x的變化率，設(shè)定誤差變化率閾值x0，則當(dāng)樣本容量導(dǎo)數(shù)的絕對值小于等于x0時(shí)，表明誤差滿足精度變化要求后，不需要增加樣點(diǎn)數(shù)，解出該不等式即可確定合適的最少樣點(diǎn)數(shù)量。

2.3 監(jiān)測樣點(diǎn)代表性檢驗(yàn)及分析方法

數(shù)據(jù)離散程度采用的測度值主要有方差、標(biāo)準(zhǔn)差、偏度和峰度等。集中趨勢的度量包括均值、中位數(shù)、重?cái)?shù)等，集中趨勢是指一組數(shù)據(jù)向中心值靠攏的程度，取決于該組數(shù)據(jù)的離散水平。本研究通過分析監(jiān)測樣點(diǎn)與耕地總體等別的集中、離散趨勢是否一致，反映監(jiān)測樣點(diǎn)對耕地總體的代表性程度高低。

對比分析監(jiān)測樣點(diǎn)在各主要屬性值上的數(shù)量比例與耕地面積比例的分布趨勢，對監(jiān)測樣點(diǎn)的代表性程度進(jìn)行檢驗(yàn)。主要的參考屬性值包括耕地的等別、地類、土壤類型、土地利用系數(shù)、土地經(jīng)濟(jì)系數(shù)等。

3 結(jié)果與分析

3.1 監(jiān)測控制區(qū)劃分

1）自然質(zhì)量分區(qū)。決定黃陵縣耕地質(zhì)量的主要自然因素是地貌和土壤。黃陵縣西部為梁峁山區(qū)，包括雙龍、腰坪、店頭3個(gè)鄉(xiāng)（鎮(zhèn)），海拔1 400 m；中部主要為河流沖積而成的河谷川道區(qū)，平均海拔800 m；東部為黃土塬區(qū)，由于沮河切割，分為南北二塬，平均海拔1 100 m。北部塬區(qū)包括隆坊、阿黨、太賢和倉村4個(gè)鄉(xiāng)（鎮(zhèn)）；南部塬區(qū)包括田莊、侯莊兩個(gè)鄉(xiāng)（鎮(zhèn)），地勢較平坦，坡度一般為0°～2°。按照地貌類型將黃陵縣劃分為3個(gè)地貌控制區(qū)。黃陵縣土壤類型分區(qū)參考《陜西省黃陵縣農(nóng)業(yè)資源調(diào)查和農(nóng)業(yè)區(qū)劃報(bào)告集》劃分為a.以黑壚土為主的塬區(qū)土壤，主要分布在田莊、侯莊、隆坊等鄉(xiāng)（鎮(zhèn)）；b.以淤土為主的河谷土壤，主要分布在洛河、沮河、河谷和主要溝道；c.以黃膳土為主的塬坡地土壤，主要分布在太賢社區(qū)以東和田莊鎮(zhèn)以東；d.以褐土為主的林灌土壤，主要分布在雙龍、腰坪和店頭等鄉(xiāng)（鎮(zhèn)）。綜合地貌控制區(qū)和土壤類型控制區(qū)，利用ArcGIS空間疊加功能，形成自然質(zhì)量控制區(qū)。

2）利用系數(shù)等值區(qū)。依據(jù)農(nóng)用地分等中確定的指定作物冬小麥、春玉米土地利用系數(shù)，劃分各指定作物土地利用系數(shù)等值區(qū)，采用空間疊加的方法，計(jì)算綜合土地利用系數(shù)，劃分利用系數(shù)等值區(qū)，根據(jù)等值區(qū)劃分土地利用水平分區(qū)。其中土地利用系數(shù)范圍為0.600 6～0.880 6，將其分為3個(gè)區(qū)間，分別為低（0.600 6～0.693 7）、中（0.693 8～0.776 5）、高（0.776 6～0.880 6）。

3）經(jīng)濟(jì)系數(shù)等值區(qū)。依據(jù)農(nóng)用地分等中確定的指定作物冬小麥、春玉米土地經(jīng)濟(jì)系數(shù)，劃分各指定作物土地經(jīng)濟(jì)系數(shù)等值區(qū)，采用與利用水平分區(qū)相同的辦法，計(jì)算綜合土地經(jīng)濟(jì)系數(shù)，劃分經(jīng)濟(jì)系數(shù)等值區(qū)，根據(jù)等值區(qū)劃分土地經(jīng)濟(jì)水平分區(qū)。其中土地經(jīng)濟(jì)系數(shù)范圍為0.665 3～0.955 8，將其分為3個(gè)區(qū)間，分別為低（0.665 3～0.771 1）、中（0.771 2～0.847 5）、高（0.847 6～0.955 8）。

4）綜合分區(qū)。采用ArcGIS空間疊置法，對“自然質(zhì)量分區(qū)”、“利用系數(shù)等值區(qū)”、“經(jīng)濟(jì)系數(shù)等值區(qū)”進(jìn)行空間疊加分析，初步劃定“監(jiān)測分區(qū)”；根據(jù)主導(dǎo)性原則，對初步劃分的“監(jiān)測區(qū)”邊界進(jìn)行適當(dāng)調(diào)整，確保各監(jiān)測區(qū)邊界不打破村界，并將面積較小的單元進(jìn)行臨近區(qū)的歸并，從而確定最終的監(jiān)測評價(jià)綜合分區(qū)。

考慮研究區(qū)行政界線，將其劃分為10個(gè)耕地質(zhì)量等別監(jiān)測控制區(qū)，并以“JC-01”的形式依次編號。每個(gè)監(jiān)測控制區(qū)具有一個(gè)因素組合類型，形式為“地貌—土壤類型—利用等值區(qū)—經(jīng)濟(jì)等值區(qū)”，如監(jiān)測控制區(qū)JC-01的為“西部梁峁山區(qū)-褐土-高-高”，主要屬性見表2。

3.2 監(jiān)測樣點(diǎn)數(shù)量確定

3.2.1 多因素限制條件法確定樣點(diǎn)數(shù)量 根據(jù)監(jiān)測樣點(diǎn)布設(shè)原則及方法，黃陵縣共布設(shè)了45個(gè)監(jiān)測樣點(diǎn)，其中11個(gè)監(jiān)測樣點(diǎn)是已設(shè)定的標(biāo)準(zhǔn)樣地，其余34個(gè)監(jiān)測樣點(diǎn)是通過多因素限制條件法新設(shè)立的監(jiān)測樣點(diǎn)。布設(shè)的45個(gè)監(jiān)測樣點(diǎn)均屬于縣域內(nèi)永久基本農(nóng)田，且均處于土地利用總體規(guī)劃劃定的禁止建設(shè)區(qū)，監(jiān)測樣點(diǎn)可長期固定監(jiān)測，且空間分布均勻，具有代表性，一定程度上能反映黃陵縣內(nèi)耕地質(zhì)量等別的動(dòng)態(tài)變化趨勢。結(jié)合年度更新評價(jià)結(jié)果，黃陵縣耕地國家自然等為10～13等，以12等地為主，占黃陵縣耕地總面積的93.05%，因此在布設(shè)監(jiān)測樣點(diǎn)時(shí)，主要布設(shè)在12等地，其他等別監(jiān)測樣點(diǎn)分布相對較少。監(jiān)測樣點(diǎn)面積最小為0.62 hm2，面積最大為20.08 hm2。具體監(jiān)測等別對應(yīng)的樣點(diǎn)情況見表3。

3.2.2 方格網(wǎng)法確定樣點(diǎn)數(shù)量 采用方格網(wǎng)法（交點(diǎn)），共布設(shè)了32個(gè)監(jiān)測樣點(diǎn)，其中11個(gè)監(jiān)測樣點(diǎn)是已設(shè)定的標(biāo)準(zhǔn)樣地作為固定樣點(diǎn)，其余21個(gè)監(jiān)測樣點(diǎn)是根據(jù)方格網(wǎng)法獲取的耕地圖斑作為動(dòng)態(tài)樣點(diǎn)。采用方格網(wǎng)法（中心點(diǎn)），共布設(shè)了37個(gè)監(jiān)測樣點(diǎn)，其中11個(gè)監(jiān)測樣點(diǎn)是以耕地質(zhì)量等級補(bǔ)充完善工作的標(biāo)準(zhǔn)樣地作為固定樣點(diǎn)，其余26個(gè)監(jiān)測樣點(diǎn)是根據(jù)方格網(wǎng)法（交點(diǎn)）獲取。endprint

通過方格網(wǎng)法（交點(diǎn)）獲取監(jiān)測樣點(diǎn)均為旱地，耕地自然質(zhì)量等別均為12等，動(dòng)態(tài)監(jiān)測樣點(diǎn)在城區(qū)街道辦和阿黨鎮(zhèn)沒有分布，動(dòng)態(tài)監(jiān)測樣點(diǎn)中面積最小為0.34 hm2，面積最大為38.41 hm2。同樣通過方格網(wǎng)法（中心點(diǎn)）獲取監(jiān)測樣點(diǎn)均為旱地，耕地自然質(zhì)量等別也均為12等，動(dòng)態(tài)監(jiān)測樣點(diǎn)在各鄉(xiāng)鎮(zhèn)均有分布，動(dòng)態(tài)監(jiān)測樣點(diǎn)中面積最小為0.09 hm2，面積最大為26.61 hm2。

3.2.3 分層抽樣法確定樣點(diǎn)數(shù)量 將黃陵縣耕地各等級成果數(shù)據(jù)代入式（2）可得，

由式（3）可以看出，y隨著x的增大而逐漸減少，也就是說抽樣精度隨著樣點(diǎn)數(shù)量的增加而逐步提高。本研究將變化率閾值設(shè)為0.5，在黃陵縣內(nèi)除去異常監(jiān)測點(diǎn)外，最少應(yīng)該布設(shè)42個(gè)監(jiān)測樣點(diǎn)。

3.2.4 監(jiān)測樣點(diǎn)最終確定 對上述3種方法獲取的監(jiān)測樣點(diǎn)結(jié)果進(jìn)行對比分析，不同之處在于通過多因素限制條件法所確定的監(jiān)測樣點(diǎn)在每個(gè)監(jiān)測控制區(qū)中數(shù)量不少于兩個(gè)，國家自然等各等別中均有分布，監(jiān)測樣點(diǎn)中耕地地類有水田、水澆地和旱地，并且交通通達(dá)度較好，便于實(shí)地監(jiān)測數(shù)據(jù)獲取。方格網(wǎng)法所確定的監(jiān)測樣點(diǎn)，在部分監(jiān)測控制區(qū)中無監(jiān)測樣點(diǎn)，國家自然等只有12等地，監(jiān)測樣點(diǎn)地類均為旱地，部分監(jiān)測樣點(diǎn)交通通達(dá)度不好，不便監(jiān)測數(shù)據(jù)獲取，樣點(diǎn)隨機(jī)性太大，這種方法適合于耕地分布均勻的地區(qū)，山區(qū)、丘陵地區(qū)一般耕地較少，不適合通過此方法進(jìn)行監(jiān)測樣點(diǎn)布設(shè)。層次分析法是對定性問題進(jìn)行定量分析的一種決策方法，權(quán)重向量計(jì)算存在偏差，得出的結(jié)果空間關(guān)聯(lián)度不高，在實(shí)際應(yīng)用中建議結(jié)合其他檢測樣點(diǎn)確定方法。

綜上所述，結(jié)合黃陵縣實(shí)際，本次黃陵縣耕地質(zhì)量等別監(jiān)測評價(jià)的監(jiān)測樣點(diǎn)選取應(yīng)用多因素限制條件法，確定監(jiān)測樣點(diǎn)數(shù)量。黃陵縣耕地等別監(jiān)測樣點(diǎn)分布情況如表4所示，監(jiān)測樣點(diǎn)分布如圖1所示。

3.3 監(jiān)測樣點(diǎn)代表性檢驗(yàn)及分析

3.3.1 監(jiān)測樣點(diǎn)代表性檢驗(yàn) 對黃陵縣耕地監(jiān)測樣點(diǎn)的代表性進(jìn)行統(tǒng)計(jì)檢驗(yàn)，借助SPSS統(tǒng)計(jì)分析軟件對黃陵縣全部耕地分等單元和監(jiān)測樣點(diǎn)的耕地等別進(jìn)行統(tǒng)計(jì)檢驗(yàn)，結(jié)果見表5。根據(jù)表5，對比分析黃陵縣全部耕地分等單元和監(jiān)測樣點(diǎn)等別的集中、離散趨勢測度值可知，兩者中值、全距、極小值、極大值均相同，分別為6、3、5、8；監(jiān)測樣點(diǎn)等別的均值為6.27，標(biāo)準(zhǔn)誤差為0.133，黃陵縣全部耕地分等單元等別的均值為6.08，標(biāo)準(zhǔn)誤差為0.004，因此總體與樣本的均值接近，差別不大，說明兩者集中趨勢一致；監(jiān)測樣點(diǎn)等別的標(biāo)準(zhǔn)差、方差分別為0.889和0.791，黃陵縣全部耕地分等單元耕地總體的標(biāo)準(zhǔn)差、方差分別為0.712和0.670，兩者差別不大；監(jiān)測樣點(diǎn)等別的偏度、峰度分別為1.059和0.346，黃陵縣監(jiān)測樣點(diǎn)和耕地分等單元總體的偏度分別為1.059和1.131，兩者偏差不太明顯，說明兩者的離散趨勢較一致?？傮w看來，黃陵縣監(jiān)測樣點(diǎn)對耕地總體的代表性較好。

3.3.2 監(jiān)測樣點(diǎn)代表性分析 對監(jiān)測樣點(diǎn)的代表性進(jìn)行分析，選取影響和反映耕地質(zhì)量的土地利用系數(shù)、土地經(jīng)濟(jì)系數(shù)、地類和等別因素，將監(jiān)測樣點(diǎn)的各因素與黃陵縣總體對比，分析其數(shù)量比例和面積比例趨勢是否一致。其中，選取冬小麥和春玉米的土地利用系數(shù)和土地經(jīng)濟(jì)系數(shù)[27]，分別進(jìn)行相關(guān)性分析及耕地的等別和地類因素的代表性分析，結(jié)果見圖2。由土地利用系數(shù)（圖2A）、土地經(jīng)濟(jì)系數(shù)（圖2B）、等別和地類的相關(guān)性（圖2C）可知，監(jiān)測樣點(diǎn)的數(shù)量分布趨勢與全縣耕地總體的面積分布趨勢基本一致，因此，監(jiān)測樣點(diǎn)對黃陵縣耕地總體的代表性程度較高。

4 小結(jié)與討論

4.1 小結(jié)

1）采用ArcGIS空間疊置方法，對影響耕地等級變化的自然狀況、利用水平和收益水平進(jìn)行疊加分析，確定監(jiān)測區(qū)。對初步劃分的“監(jiān)測區(qū)”邊界進(jìn)行適當(dāng)調(diào)整，確保各監(jiān)測區(qū)邊界不打破村界，并將面積較小的單元進(jìn)行臨近區(qū)的歸并，從而確定最終的監(jiān)測評價(jià)綜合分區(qū)，符合統(tǒng)計(jì)學(xué)空間處理方式。

2）在監(jiān)測控制區(qū)內(nèi)布設(shè)監(jiān)測樣點(diǎn)，采用多因素限制條件法、網(wǎng)格法和分層抽樣法進(jìn)行監(jiān)測樣點(diǎn)的布設(shè)，研究結(jié)果顯示，多因素限制條件法配合分層抽樣法確定監(jiān)測樣點(diǎn)，可更好地滿足監(jiān)測樣點(diǎn)布設(shè)代表性的精度要求。

3）黃陵縣共劃定45個(gè)監(jiān)測樣點(diǎn)，其中11個(gè)監(jiān)測樣點(diǎn)是以耕地質(zhì)量等級補(bǔ)充完善工作的標(biāo)準(zhǔn)樣地作為固定樣點(diǎn)，其余34個(gè)監(jiān)測樣點(diǎn)是通過多因素限制條件法獲取。對研究區(qū)監(jiān)測樣點(diǎn)等別的集中、離散程度趨勢與研究區(qū)耕地總體等別的分布趨勢一致，在影響耕地質(zhì)量主要因素上，監(jiān)測樣點(diǎn)數(shù)量的分布趨勢與耕地總體的面積分布趨勢較一致。

4.2 討論

1）方格網(wǎng)法所確定的監(jiān)測樣點(diǎn)適合于耕地分布均勻的地區(qū)，山區(qū)、丘陵地區(qū)一般耕地較少；由于國家自然等只有12等地，監(jiān)測樣點(diǎn)地類全為旱地，導(dǎo)致本研究中在部分監(jiān)測控制區(qū)無監(jiān)測樣地，且存在部分監(jiān)測樣點(diǎn)交通通達(dá)度不好、不便監(jiān)測數(shù)據(jù)獲取、樣點(diǎn)隨機(jī)性太大等問題，因此，這種方法不適合在黃陵縣進(jìn)行監(jiān)測樣點(diǎn)布設(shè)。分層抽樣法得出的結(jié)果空間關(guān)聯(lián)度不高，權(quán)重向量計(jì)算存在偏差，所以并未采用。本研究提出多因素限制條件法，監(jiān)測樣點(diǎn)中耕地地類有水田、水澆地和旱地，并且交通通達(dá)度較好，便于實(shí)地監(jiān)測數(shù)據(jù)獲取，結(jié)合分層抽樣對樣點(diǎn)數(shù)量下限的確定，更能滿足黃陵縣實(shí)際需要。

2）考慮到黃陵縣地貌復(fù)雜，不同位置的耕地對于農(nóng)民耕作也有一定影響。因此，建議增加排水條件、灌溉水源、交通便捷度和耕作便利度4個(gè)分等因素參與耕地自然質(zhì)量分等的計(jì)算，權(quán)重根據(jù)各年度的監(jiān)測數(shù)據(jù)具體分析確定，其準(zhǔn)確性更符合實(shí)際。

3）根據(jù)《農(nóng)用地質(zhì)量分等規(guī)程》（GB/T 28407-2012）要求，每2～3年應(yīng)調(diào)整土地經(jīng)濟(jì)系數(shù)，但為了保障耕地質(zhì)量等別年度更新評價(jià)成果的延續(xù)性和可比性，建議對土地經(jīng)濟(jì)系數(shù)和土地利用系數(shù)在每6年全面更新耕地質(zhì)量分等成果時(shí)再重新調(diào)查樣點(diǎn)，根據(jù)樣點(diǎn)數(shù)據(jù)分析修正土地經(jīng)濟(jì)系數(shù)和土地利用系數(shù)等值區(qū)。endprint

4）建立縣域耕地質(zhì)量等別監(jiān)測樣點(diǎn)體系，對監(jiān)測樣點(diǎn)進(jìn)行定位監(jiān)測和動(dòng)態(tài)監(jiān)測，可反映耕地質(zhì)量變化情況。同時(shí)，當(dāng)年新增耕地、毀損土地等的監(jiān)測需另外布設(shè)動(dòng)態(tài)監(jiān)測樣點(diǎn)，可為實(shí)現(xiàn)實(shí)時(shí)動(dòng)態(tài)監(jiān)測提供可能。

參考文獻(xiàn)：

[1] 孫江鋒，張合兵，劉明遠(yuǎn)，等.耕地動(dòng)態(tài)監(jiān)測應(yīng)用管理系統(tǒng)框架設(shè)計(jì)及實(shí)現(xiàn)[J].中國農(nóng)業(yè)資源與區(qū)劃，2015，36（2）：54-60.

[2] 孔祥斌，張青璞.中國西部區(qū)耕地等別空間分布特征[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（22）：1-7.

[3] 孔祥斌，張青璞，鄖文聚，等.基于平衡轉(zhuǎn)換的國家級農(nóng)用地分等成果匯總方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29（16）：238-246.

[4] 鄖文聚.我國耕地資源開發(fā)利用的問題與整治對策[J].中國科學(xué)院院刊，2015，30（4）：484-491.

[5] 張 超，喬 敏，鄖文聚，等.耕地?cái)?shù)量、質(zhì)量、生態(tài)三位一體綜合監(jiān)管體系研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2017，48（1）：1-6.

[6] 余述瓊，張蚌蚌，相 慧，等.基于因素組合的耕地質(zhì)量等級監(jiān)測樣點(diǎn)布控方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（24）：288-297.

[7] 吳克寧，韓春建，馮新偉，等.基于3S技術(shù)的土地生態(tài)敏感性分區(qū)研究[J].土壤，2008，40（2）：293-298.

[8] 吳克寧，史原軻，馮新偉，等.基于農(nóng)用地分等的城區(qū)擴(kuò)展用地空間布局優(yōu)化研究[J].中國土地科學(xué)，2007，21（6）：17-22，16.

[9] 程 鋒，王洪波，鄖文聚.中國耕地質(zhì)量等級調(diào)查與評定[J].中國土地科學(xué)，2014，28（2）：75-82，97.

[10] NUSSER S M，BREIDT F J，F(xiàn)ULLER W A. Design and estimation for investigating the dynamics of natural resources[J].Journal of Ecological Applications，1998，8（2）：234-245.

[11] JONES R J A，SPOOR G，THOMASSON A J. Assessing the vulnerability of subsoils in Europe to compaction：A preliminary analysis[J].Soil Tillage Research，2003，73（1）：131-143.

[12] HUFFMAN E，EILERS R G，PADBURY G，et al. Canadian agri-environmental indicators related to land quality： Integrating census and biophysical data to estimate soil cover， wind erosion and soil salinity[J].Agriculture， Ecosystems and Environment，2000，81（2）：113-123.

[13] 張蚌蚌，孔祥斌，相 慧，等.中國西部區(qū)國家級標(biāo)準(zhǔn)樣地分布特征及代表性[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（7）：227-235.

[14] 吳克寧，焦雪瑾，梁思源，等.基于標(biāo)準(zhǔn)樣地國家級匯總的耕地質(zhì)量動(dòng)態(tài)監(jiān)測點(diǎn)構(gòu)架研究[J].農(nóng)業(yè)工程學(xué)報(bào)，2008，24（10）：74-79.

[15] 楊建宇，岳彥利，宋海榮，等.基于空間模擬退火算法的耕地質(zhì)量布樣及優(yōu)化方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2015，31（20）：253-261.

[16] 殷守強(qiáng)，王 鑫，賀文龍，等.耕地質(zhì)量監(jiān)測中不同抽樣方法的精度比較——以河北省黃驊市為例[J].資源科學(xué)，2016，38（11）：2049-2057.

[17] 張玉臻，孔祥斌，劉 焱，等.基于標(biāo)準(zhǔn)樣地的省級耕地質(zhì)量監(jiān)測樣地布設(shè)方法——以內(nèi)蒙古自治區(qū)為例[J].資源科學(xué)，2016，38（11）：2037-2048.

[18] 王 倩，尚月敏，馮 銳，等.基于變異函數(shù)的耕地質(zhì)量等別監(jiān)測點(diǎn)布設(shè)分析——以四川省中江縣和北京市大區(qū)為例[J].中國土地科學(xué)，2012，26（8）：80-86.

[19] 孫亞彬，吳克寧，胡曉濤，等.基于潛力指數(shù)組合的耕地質(zhì)量等級監(jiān)測布點(diǎn)方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29（4）：245-254.

[20] 胡曉濤，吳克寧，馬建輝，等.北京市大興區(qū)耕地質(zhì)量等級監(jiān)測控制點(diǎn)布設(shè)[J].資源科學(xué)，2012，34（10）：1891-1896.

[21] 蔡鷺斌，孔祥斌，段建南，等.國外經(jīng)驗(yàn)對中國耕地質(zhì)量監(jiān)測布點(diǎn)的啟示[J].中國農(nóng)學(xué)通報(bào)，2014，30（14）：192-197.

[22] 楊建宇，岳彥利，宋海榮，等.基于空間平衡法的縣域耕地質(zhì)量監(jiān)測布樣方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2015，31（24）：274-280.

[23] 余述瓊，張蚌蚌，相 慧，等.基于因素組合的耕地質(zhì)量等級監(jiān)測樣點(diǎn)布控方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（24）：288-297.

[24] 黨福江，戈素芬.溝壑密度調(diào)查應(yīng)采取水文網(wǎng)絡(luò)法[J].水土保持通報(bào)，1998，18（6）：31-33.

[25] 王勁峰.空間抽樣與統(tǒng)計(jì)推斷[M].北京：科學(xué)出版社，2009.

[26] 孫樹東.初等函數(shù)的幾點(diǎn)探討[J].西昌學(xué)院學(xué)報(bào)（自然科學(xué)版），2010，24（4）：32-34.

[27] 馮蓉曄，談志浩，黃勁松，等.農(nóng)用地分等中土地利用系數(shù)與經(jīng)濟(jì)系數(shù)計(jì)算方法改進(jìn)探討——以江蘇省無錫市為例[J].經(jīng)濟(jì)地理，2004，24（2）：246-249.endprint