寧夏耕地質(zhì)量分等成果在耕地占補(bǔ)平衡中的應(yīng)用

摘 要：耕地是人類賴以生存和發(fā)展的基礎(chǔ)，如何充分利用耕地質(zhì)量成果更好地服務(wù)社會(huì)和保護(hù)耕地，成為當(dāng)下的重要研究課題。本文通過(guò)建立利用等指數(shù)—標(biāo)準(zhǔn)糧食產(chǎn)量線性回歸模型，以寧夏回族自治區(qū)平原、山地丘陵及山地丘陵與平原耕地3種類型為例，研究耕地利用等成果在耕地占補(bǔ)平衡中的應(yīng)用，得到不同耕地類型區(qū)自治區(qū)級(jí)耕地占補(bǔ)平衡利用等別折算系數(shù)，以期為土地管理和耕地占補(bǔ)平衡政策的貫徹落實(shí)提供理論依據(jù)。

關(guān)鍵詞：耕地;質(zhì)量;占補(bǔ)平衡;線性回歸模型;寧夏

改革開放以來(lái)，我國(guó)經(jīng)濟(jì)出現(xiàn)了長(zhǎng)期的高速增長(zhǎng)態(tài)勢(shì)，極大地刺激了社會(huì)經(jīng)濟(jì)發(fā)展對(duì)土地資源的需求，導(dǎo)致耕地?cái)?shù)量不斷減少。與此同時(shí)，耕地保護(hù)制度不斷被建立和完善。1999—2009年，我國(guó)第一輪全國(guó)農(nóng)用地分等完成。2011—2013年，我國(guó)開展了耕地質(zhì)量等別補(bǔ)充完善工作，形成了基于二次調(diào)查及變更調(diào)查數(shù)據(jù)的最新耕地質(zhì)量等別成果。自2014年起，對(duì)耕地質(zhì)量實(shí)行年度更新評(píng)價(jià)，建立了能反映耕地自然質(zhì)量、平均利用和經(jīng)濟(jì)效益水平于一體的耕地質(zhì)量等別評(píng)價(jià)體系[1]。根據(jù)國(guó)家最新相關(guān)政策部署，將適時(shí)更新和完善耕地質(zhì)量等別成果，其成果應(yīng)用的重要性不言而喻[2]。

隨著社會(huì)經(jīng)濟(jì)的發(fā)展，寧夏回族自治區(qū)耕地后備資源不足和儲(chǔ)備數(shù)量不斷下降，中南部落實(shí)“先補(bǔ)后占”“占優(yōu)補(bǔ)優(yōu)”的困難不斷加大。在此背景下，寧夏回族自治區(qū)出臺(tái)了相關(guān)文件，政策的出臺(tái)對(duì)于滿足全區(qū)社會(huì)經(jīng)濟(jì)建設(shè)用地需求和耕地保護(hù)找到了新的突破口，但政策的落地需要一定的制度保障，目前跨縣（市、區(qū)）、跨地級(jí)市補(bǔ)充耕地和耕地改造后的耕地質(zhì)量平衡缺乏統(tǒng)一的折算標(biāo)準(zhǔn)和參考依據(jù)。本文根據(jù)寧夏回族自治區(qū)實(shí)際，提出基于新形勢(shì)下耕地質(zhì)量等別評(píng)價(jià)成果應(yīng)用的對(duì)策和方向，以期對(duì)寧夏回族自治區(qū)耕地保護(hù)相關(guān)政策制定起到一定啟發(fā)。

一、研究區(qū)概況

寧夏回族自治區(qū)是我國(guó)西北地區(qū)重要的產(chǎn)糧基地，隨著社會(huì)經(jīng)濟(jì)的發(fā)展，建設(shè)占用耕地和對(duì)耕地后備資源的持續(xù)開發(fā)，“十二五”期間寧夏回族自治區(qū)建設(shè)占用耕地14 470 hm2、補(bǔ)充耕地27 556 hm2及土地整治288 558 hm2。根據(jù)2014年寧夏回族自治區(qū)土地利用現(xiàn)狀變更調(diào)查成果，耕地總面積為1 290 121.88 hm2，其中水田188 002.77 hm2、水澆地320 329.22 hm2、旱地781 789.89 hm2。根據(jù)《農(nóng)用地分等規(guī)程》（TD-T 1004—2016），寧夏回族自治區(qū)劃分為3個(gè)三級(jí)分等因素指標(biāo)區(qū)，基準(zhǔn)作物均為小麥，指定作物為玉米、水稻和馬鈴薯，農(nóng)作物耕作制度為一年一熟。

二、研究方法

（一）樣點(diǎn)選取方法

根據(jù)樣點(diǎn)選取要遵循的代表性、適量性原則，本研究采用分層隨機(jī)抽樣調(diào)查的方法從2014年寧夏回族自治區(qū)耕地質(zhì)量等別年度評(píng)價(jià)成果分等單元中提取樣本。首先分地類、分等別統(tǒng)計(jì)面積及其圖斑個(gè)數(shù)，其次以圖斑個(gè)數(shù)10%作為總量，按地類、等別面積加權(quán)分配總量圖斑作為各等別樣本個(gè)數(shù)，最后在Excel表格中根據(jù)寧夏回族自治區(qū)匯總的耕地質(zhì)量等別評(píng)價(jià)成果，按地類、等別隨機(jī)篩選樣本。

（二）標(biāo)準(zhǔn)糧產(chǎn)量計(jì)算方法

標(biāo)準(zhǔn)糧產(chǎn)量是基準(zhǔn)作物的產(chǎn)量，根據(jù)調(diào)查樣地指定作物產(chǎn)量，運(yùn)用產(chǎn)量比系數(shù)計(jì)算分等單元標(biāo)準(zhǔn)糧的實(shí)際產(chǎn)量[3]，如式（1）所示：

（1）

式（1）中，P表示樣點(diǎn)分等單元標(biāo)準(zhǔn)糧產(chǎn)量;n指樣點(diǎn)分等單元個(gè)數(shù);Wi指第i個(gè)樣點(diǎn)分等單元指定作物實(shí)際產(chǎn)量;fi指第i個(gè)樣點(diǎn)分等單元指定作物的產(chǎn)量比系數(shù)。

（三）平均利用等指數(shù)計(jì)算方法

為了準(zhǔn)確計(jì)算樣點(diǎn)分等單元等別平均利用等指數(shù)，利用面積加權(quán)平均算法計(jì)算得到樣點(diǎn)在每個(gè)利用等別條件下的平均利用等指數(shù)[4]，如式（2）所示：

（2）

式（2）中，Yj表示第j個(gè)利用等別的平均利用等指數(shù);n代表樣點(diǎn)分等單元的個(gè)數(shù);Yij是第j個(gè)利用等別中第i個(gè)樣點(diǎn)的利用等指數(shù);Sij是第j個(gè)利用等別中第i個(gè)樣點(diǎn)的單元面積;Sj是第j個(gè)利用等別中樣點(diǎn)的總面積。

（四）等級(jí)折算系數(shù)計(jì)算方法

根據(jù)標(biāo)準(zhǔn)糧產(chǎn)量與利用等指數(shù)之間的函數(shù)關(guān)系編制等別折算系數(shù)表[5]，公式為：

式（3）中，Mi表示新補(bǔ)充耕地及被占用耕地的等別折算系數(shù);V占i表示被占用的第i等耕地的標(biāo)準(zhǔn)糧產(chǎn)量;V占j則為補(bǔ)充的第j等耕地的標(biāo)準(zhǔn)糧產(chǎn)量。

三、結(jié)果與分析

（一）樣本分析

根據(jù)耕地自治區(qū)級(jí)利用等指數(shù)的計(jì)算方法，寧夏回族自治區(qū)耕地自治區(qū)級(jí)利用等指數(shù)分布在2～1 290，按照等間距法，每100分一個(gè)等，等指數(shù)越大，等別越高，耕地質(zhì)量越好，寧夏回族自治區(qū)耕地自治區(qū)級(jí)利用等共劃分為13個(gè)等別，等別分布為1～13等。

按照樣本的選取方法，從寧夏回族自治區(qū)181 301個(gè)有效分等單元中選擇了18 131個(gè)自治區(qū)區(qū)級(jí)利用等分等單元，建立分析樣本數(shù)據(jù)庫(kù)，樣本數(shù)占全部有效分等單元比例10%，符合統(tǒng)計(jì)檢驗(yàn)樣本數(shù)的要求。

（二）不同區(qū)域的耕地利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量之間的關(guān)系分析

以耕地利用等指數(shù)作為自變量X，標(biāo)準(zhǔn)糧產(chǎn)量作為因變量Y，根據(jù)調(diào)查數(shù)據(jù)，每一個(gè)樣點(diǎn)對(duì)應(yīng)一個(gè)利用等指數(shù)和標(biāo)準(zhǔn)糧產(chǎn)量，得出利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量之間的數(shù)據(jù)散點(diǎn)圖。

1.平原耕地利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量的關(guān)系分析

通過(guò)線性回歸分析方法計(jì)算得出，平原耕地總體樣點(diǎn)利用等指數(shù)與標(biāo)準(zhǔn)糧之間的函數(shù)關(guān)系模型Ⅰ為：Y=5.168X+

1 142。其中，判定系數(shù)R2=0.810，即變量Y的變化中有81.00%是由變量X引起的，二者具有一定的相關(guān)性，線性擬合度較好，相關(guān)性顯著（見圖1）。

2.山地丘陵坡耕地利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量的關(guān)系分析

通過(guò)線性回歸分析方法計(jì)算得出，山地丘陵坡總體樣點(diǎn)利用等指數(shù)與標(biāo)準(zhǔn)糧之間的函數(shù)關(guān)系模型Ⅱ?yàn)椋篩=6.762X+841.5。其中，判定系數(shù)R2=0.800，即變量Y的變化中有80.00%是由變量X引起的，二者具有一定的相關(guān)性，線性擬合度較好，相關(guān)性顯著（見圖2）。

3.平原與山地丘陵耕地利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量的關(guān)系分析

通過(guò)線性回歸分析方法計(jì)算得出，平原與山地丘陵耕地總體樣點(diǎn)利用等指數(shù)與標(biāo)準(zhǔn)糧之間的函數(shù)關(guān)系模型Ⅲ為：Y=5.234X+1 125。其中，判定系數(shù)R2=0.825，即變量Y的變化中有82.50%是由變量X引起的，二者具有一定的相關(guān)性，線性擬合度較好，相關(guān)性顯著（見圖3）。

通過(guò)式（2）計(jì)算出各等別平均利用等指數(shù)之后，將其代入函數(shù)關(guān)系模型Ⅰ、Ⅱ、Ⅲ，即可得到平原耕地、山地丘陵坡耕地以及平原與山地丘陵耕地3種耕地類型區(qū)各利用等別所對(duì)應(yīng)的理論標(biāo)準(zhǔn)糧產(chǎn)量，見表1。

根據(jù)等級(jí)折算公式（3），分別測(cè)算出平原耕地、山地丘陵耕地及平原與山地丘陵耕地的等別折算系數(shù)，最終得到寧夏回族自治區(qū)山地丘陵耕地?cái)?shù)量質(zhì)量等別折算系數(shù)平均級(jí)差約為0.20，4等地與8等地之間的折算系數(shù)為1.78，即占用8等地需要補(bǔ)充1.78倍的4等耕地;寧夏回族自治區(qū)平原耕地質(zhì)量4等地與8等地之間的折算系數(shù)為1.70;寧夏回族自治區(qū)平原與山地丘陵耕地質(zhì)量4等地與8等地之間的折算系數(shù)為1.71，其他等別間的折算以此類推，可以看出等級(jí)折算系數(shù)隨著土地利用水平的提高逐漸增大。

四、結(jié)論與討論

（1）該文利用統(tǒng)計(jì)、分析的方法建立標(biāo)準(zhǔn)糧產(chǎn)量與利用等指數(shù)線性回歸模型，充分驗(yàn)證了寧夏回族自治區(qū)利用等指數(shù)與標(biāo)準(zhǔn)糧產(chǎn)量顯著相關(guān)。通過(guò)測(cè)算得出了區(qū)內(nèi)和區(qū)間不同耕地類型利用等級(jí)耕地折算系數(shù)，與寧夏回族自治區(qū)的實(shí)際情況相符。從變化規(guī)律來(lái)看，隨著土地利用水平的提高，其等級(jí)折算系數(shù)呈現(xiàn)遞增趨勢(shì)，這一變化規(guī)律與土地的利用規(guī)律相吻合。利用此等級(jí)折算系數(shù)執(zhí)行占用耕地補(bǔ)償?shù)南嚓P(guān)制度，能有效防止占優(yōu)補(bǔ)劣、占多補(bǔ)少，保證了耕地占補(bǔ)過(guò)程中的數(shù)量和質(zhì)量平衡。

（2）該文對(duì)耕地占補(bǔ)平衡折算系數(shù)計(jì)算是理論模型，實(shí)際占補(bǔ)平衡過(guò)程中可能存在占用一塊耕地補(bǔ)充幾塊耕地或者占用多塊耕地補(bǔ)充一塊耕地及占多塊補(bǔ)多塊的情況。這就需要在確定被占耕地等別與面積的情況下，按照補(bǔ)充耕地實(shí)際，結(jié)合各耕地類型區(qū)域內(nèi)和耕地類型區(qū)域間的等別折算系數(shù)表，綜合確定補(bǔ)充耕地等別和面積。這無(wú)疑增大了操作難度，折算系數(shù)如何科學(xué)合理應(yīng)用有待于進(jìn)一步研究。

參考文獻(xiàn)：

[1]魏洪斌.基于土地整治的耕地質(zhì)量評(píng)價(jià)與提升研究[D].北京：中國(guó)地質(zhì)大學(xué)，2015.

[2]程鋒，王洪波，鄖文聚.中國(guó)耕地質(zhì)量等級(jí)調(diào)查與評(píng)價(jià)[J].中國(guó)土地科學(xué)，2014（2）：83-89.

[3]匡麗花，葉英聰，趙小敏，等.基于農(nóng)用地分等修正的土地整治項(xiàng)目對(duì)耕地質(zhì)量的影響評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào)，2016（17）：198-205.

[4]楊偉，謝德體，廖和平，等.基于高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)模式的農(nóng)用地整治潛力分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2013（7）：219-229.

[5]路婕，吳克寧，鄖文聚，等.全國(guó)統(tǒng)一的耕地占補(bǔ)平衡按等級(jí)折算系數(shù)研究[J].資源與產(chǎn)業(yè)，2012（4）：128-134.