李俊穎,吳克寧※,宋 文
(1. 中國地質(zhì)大學(北京)土地科學技術學院,北京 100083;2. 國土資源部土地整治重點實驗室,北京 100035;3. 山東農(nóng)業(yè)大學資源與環(huán)境學院,泰安 271018)
農(nóng)用地分等定級是通過定量指標科學量化農(nóng)用地數(shù)量、質(zhì)量及其分布[1],對農(nóng)用地的質(zhì)量優(yōu)劣進行綜合評定并劃分等級[2],是實現(xiàn)耕地占補平衡、保證耕地總量動態(tài)平衡的有效措施。1989年原國家土地管理局制定《農(nóng)用地分等定級規(guī)程(征求意見稿)》,并以此為依據(jù)在全國進行試點研究;2003年國土資源部在試點經(jīng)驗總結(jié)的基礎上經(jīng)過多次研討和修正,正式頒布《農(nóng)用地分等規(guī)程》(TD/T1004-2003)[3],并于2012年修改頒布《農(nóng)用地質(zhì)量分等規(guī)程》(GB/T28407-2012)[4],奠定了全國農(nóng)用地分等工作的基礎。在此過程中,相關專家學者針對農(nóng)用地分等工作中出現(xiàn)的問題進行了大量研究。陳安慶、王建國等[5-6]探討了農(nóng)用土地分等定級的基本理論和方法;張鳳榮等[7-8]分析制定了農(nóng)用地分等定級野外診斷指標體系,并且構(gòu)建了農(nóng)用地分等土壤質(zhì)量指標體系;安萍莉等[9]提出了標準耕作制度的定義,并探討了確定標準耕作制度的原則和方法。楊建鋒等[10]分析了自然質(zhì)量分的計算方法以及自然質(zhì)量分的可比對分等成果可比的作用;孔祥斌、伍育鵬等[11-12]對產(chǎn)量比系數(shù)的求算方法進行了對比分析。
農(nóng)用地分等的核心是全國可比的分等指數(shù)計算[13]和等級體系建立[14]。在分等指數(shù)計算過程中產(chǎn)量比系數(shù)是直接影響農(nóng)用地分等結(jié)果及全國可比性的重要參數(shù),其作用是把糧食折算為標準糧,以實現(xiàn)全國可比[15]。在《農(nóng)用地質(zhì)量分等規(guī)程》(GB/T28407-2012)(下稱《規(guī)程》)中列舉了利用作物最高單產(chǎn)求取產(chǎn)量比系數(shù)的計算方法,最高產(chǎn)量比法因其簡單易行而被廣泛利用。產(chǎn)量比系數(shù)的計算與各地區(qū)的標準耕作制度以及耕作環(huán)境等條件有很大的關系[16],由于不同的工作者對《規(guī)程》中所規(guī)定的產(chǎn)量比系數(shù)的理解不同,以及不同地區(qū)的分等因素和標準耕作制度不同,其基準作物和指定作物有所不同,導致了農(nóng)用地分等成果出現(xiàn)了分值相同但內(nèi)涵不同的現(xiàn)象[17],無法完全達到全國可比,因而改進產(chǎn)量比系數(shù)的計算方法是有必要的。王秋香[18]提出為保證分等結(jié)果的可比性,建議以糧食作物的卡路里熱量和干蛋白含量等生物單位作為折算依據(jù)。以此為出發(fā)點,考慮到不同地區(qū)基準作物和指定作物雖然不同,但在作物供給人體所需能量和營養(yǎng)的營養(yǎng)素層面是可以統(tǒng)一可比的。以此為依據(jù),本文選取蛋白質(zhì)、脂肪、碳水化合物、礦物質(zhì)、維生素等20種營養(yǎng)素指標和9種標準耕作制度規(guī)定作物,構(gòu)建營養(yǎng)當量修正法,通過計算不同作物“營養(yǎng)當量”來修正產(chǎn)量比系數(shù),并選取北方以小麥為基準作物的河北省曲周縣土地整治項目區(qū)和南方以水稻為基準作物的江蘇省宜興市丁蜀鎮(zhèn)作為案例進行實證,對比通過營養(yǎng)當量修正法和最高產(chǎn)量比法計算產(chǎn)量比系數(shù)所得出的農(nóng)用地自然等的差異。以期為完善農(nóng)用地質(zhì)量分等提供借鑒,為相關研究提供參考。
研究區(qū)選擇河北省邯鄲市曲周縣東北部土地整治項目區(qū)和江蘇省宜興市丁蜀鎮(zhèn)。曲周縣位于河北省南部,邯鄲市東北部,36°35′~36°57′N,114°50′~115°13′E 之間,總面積676.695 km2。在全國標準耕作制度分區(qū)中屬于黃淮海區(qū)下的冀魯豫低洼平原區(qū)。項目區(qū)為河北沖積平原,受河流擺動淤積等影響,地表形態(tài)除多為低平地外,還有局部的二坡地、緩崗、決口沖積扇和洼地,海拔為32.7 m至45.4 m。年平均氣溫13.1 ℃,≥10 ℃積溫4 472 ℃,無霜期201 d,年均降水量556.2 mm,土壤以潮土為主,屬河北省農(nóng)用地分等因素三級指標區(qū)的“冀南低平原冬小麥夏玉米一年兩熟區(qū)”。丁蜀鎮(zhèn)位于江蘇省宜興市,31°07′~31°37′N,119°31′~120°03′E 之間,地處長江三角洲經(jīng)濟開發(fā)區(qū),東瀕太湖,西部為天目山余脈,總面積205 km2。地處北亞熱帶季風氣候,全年溫暖濕潤,氣候溫和、干濕冷暖、四季分明、雨量充沛、日照充足、無霜期長,年平均氣溫 15.7 ℃,無霜期 240多天,生長期可達 250 d左右,年平均降水量為1 199.5 mm。在全國標準耕作制度分區(qū)中屬于江蘇省農(nóng)用地分等因素三級指標區(qū)的“寧鎮(zhèn)揚丘陵小麥水稻一年兩熟區(qū)”。
1.2.1 農(nóng)用地分等指標及權(quán)重
遵循農(nóng)用地分等評價指標選擇的“主導因素原則”“穩(wěn)定性原則”“差異性原則”和“實際性原則”[19],依據(jù)兩地的不同的自然條件、標準耕作制度,綜合考慮評價農(nóng)用地等級的影響因子確定相應的分等指標,并分別采用德爾菲法[20]確定指標權(quán)重。曲周縣項目區(qū)選取表層土壤質(zhì)地、排水條件、鹽漬化程度、灌溉保證率、土壤有機質(zhì)含量、剖面構(gòu)型6個分等指標,數(shù)據(jù)來源于曲周縣土壤圖、土壤有機質(zhì)含量分布圖、土地利用現(xiàn)狀圖等,指標分值及權(quán)重見表1;丁蜀鎮(zhèn)選取表層土壤質(zhì)地、耕作層厚度、障礙層距地表深度、土壤有機質(zhì)含量、土壤pH值、灌溉保證率、排水條件、土壤侵蝕程度等8個分等指標,數(shù)據(jù)來源于宜興市 2011年度土地利用現(xiàn)狀變更調(diào)查數(shù)據(jù)庫、農(nóng)田設施基本狀況、農(nóng)業(yè)經(jīng)濟資料及土壤自然質(zhì)量資料等,指標分值及權(quán)重見表2。
表1 曲周縣項目區(qū)分等指標分值及權(quán)重Table 1 Score and weight of classification indexes in Quzhou county project area
表2 丁蜀鎮(zhèn)分等指標分值及權(quán)重Table 2 Score and weight of classification indexes in Dingshu town
1.2.2 營養(yǎng)素指標和分等作物的選取
人體正常運轉(zhuǎn)需要很多的營養(yǎng)素,本文選取營養(yǎng)素的標準一是為目前已明確研究并納入中國食物成分表檢測范圍且已給出參考攝入量標準的營養(yǎng)素[21];二是在維持人類健康所必需的蛋白質(zhì)、脂類、糖類、維生素、礦物質(zhì)、膳食纖維等營養(yǎng)素之中。以18歲及以上輕體力活動成年男子膳食營養(yǎng)素平均需要量為參考攝入量標準,以100 g為計算單位,依據(jù)《中國食物成分表(第二版)》[22]和《中國居民膳食營養(yǎng)素參考攝入量(2013版)》[23]選取脂肪、碳水化合物、蛋白質(zhì)、膳食纖維、維生素 A、維生素B1、維生素B2、維生素C、維生素E、鈉、鈣、鐵、煙酸、磷、鉀、鎂、鋅、硒、錳、銅共20種營養(yǎng)素指標,同時根據(jù)《規(guī)程》選取 9種標準耕作制度規(guī)定作物,對作物的營養(yǎng)價值進行評價。
1.2.3 作物營養(yǎng)素含量及參考攝入量
表 3為作物中營養(yǎng)素含量及其日均參考攝入量,營養(yǎng)素參考攝入量標準參考《中國居民膳食營養(yǎng)素參考攝入量(20l3版)》中l(wèi)8歲及以上輕體力活動成年男子膳食營養(yǎng)素日均平均需要量。如無日均平均需要量,則可依次選取日均推薦攝入量、日均適宜攝入量或日均特定建議值。
表3 100 g作物中營養(yǎng)素含量及各營養(yǎng)素日均參考攝入量Table 3 Nutrient content of crop (100 g) and nutrients reference daily value
營養(yǎng)當量,即營養(yǎng)相當?shù)牧?,可以反映食物中所含營養(yǎng)素的總體情況,以食物中營養(yǎng)素密度的算數(shù)均數(shù)計算營養(yǎng)當量值,是根據(jù)食物鮮質(zhì)量中的各營養(yǎng)素的種類及含量作為衡量食物中不同營養(yǎng)素總體含量的評價指標,評價食物對機體所需營養(yǎng)素的綜合供給能力??赏ㄟ^比較不同食物的營養(yǎng)當量大小,以對不同食物營養(yǎng)價值做出綜合評價和判斷[24],作物營養(yǎng)素密度及營養(yǎng)當量計算公式[21]如下:
營養(yǎng)素密度=100 g食物鮮質(zhì)量中某營養(yǎng)素含量/該營養(yǎng)素參考攝入量標準
營養(yǎng)當量=∑營養(yǎng)素密度/n,其中n=20。
表4 作物的營養(yǎng)素密度及營養(yǎng)當量(100 g)Table 4 Nutrient density and nutrition equivalent unit of crops (100 g)
營養(yǎng)當量系數(shù)(γ)=指標區(qū)內(nèi)基準作物的營養(yǎng)當量/指定作物的營養(yǎng)當量
根據(jù)表4作物營養(yǎng)素密度及營養(yǎng)當量,以小麥的營養(yǎng)當量為基準計算其他作物的營養(yǎng)當量系數(shù),結(jié)果如表5。
表5 營養(yǎng)當量系數(shù)Table 5 Nutrition equivalent unit coefficient
2.2.1 確定光溫生產(chǎn)潛力指數(shù)
有灌溉條件的水田、水澆地用光溫生產(chǎn)潛力評價,地形、地貌及水源等限制因素導致無灌溉條件的旱地用氣候生產(chǎn)潛力評價[25],曲周項目區(qū)與丁蜀鎮(zhèn)兩地農(nóng)用地均具備灌溉條件,因而采用光溫生產(chǎn)潛力進行評價,通過查詢《規(guī)程》確定兩個研究區(qū)所處二級指標區(qū)指定作物的光溫生產(chǎn)潛力。
2.2.2 確定產(chǎn)量比系數(shù)
1)最高產(chǎn)量比法產(chǎn)量比系數(shù)
指定作物產(chǎn)量比系數(shù)=基準作物最高單產(chǎn)/指定作物最高單產(chǎn)
產(chǎn)量比系數(shù)是表示在相同光、溫、水、土資源條件下,不同作物產(chǎn)量之間的當量關系。因此,如基準作物與指定作物是在同一生長季節(jié)內(nèi)播種生長,則作物單產(chǎn)可直接用來計算產(chǎn)量比系數(shù);如基準作物與指定作物是在不同季節(jié)播種生長,則將指定作物或基準作物的最高單產(chǎn)調(diào)整為基準作物或指定作物生長發(fā)育所對應光、溫、水條件下的最高單產(chǎn),再通過上述公式計算作物的產(chǎn)量比系數(shù)[26]。
2)營養(yǎng)當量修正法產(chǎn)量比系數(shù)
營養(yǎng)當量修正法是在原有的最高產(chǎn)量比系數(shù)基礎上乘以營養(yǎng)當量系數(shù),將產(chǎn)量比系數(shù)由單位面積產(chǎn)量比修正為單位面積營養(yǎng)比,營養(yǎng)當量修正法產(chǎn)量比系數(shù)的計算公式為:式中β新ij為修正后第j種作物的產(chǎn)量比系數(shù);βij為最高產(chǎn)量比法確定的第j種作物的產(chǎn)量比系數(shù):γij為第i個分等單元內(nèi)第j種指定作物的營養(yǎng)當量系數(shù);j為指定作物種數(shù)。
2.2.3 計算分等單元自然質(zhì)量分
采用加權(quán)平均法計算分等單元指定作物的農(nóng)用地自然質(zhì)量分,計算公式為
式中CLij為分等單元指定作物的農(nóng)用地自然質(zhì)量分;i為分等單元編號;j為指定作物編號;k為分等因素編號:m為分等指標的數(shù)目,ωk為分等指標的權(quán)重;fijk為第 i個分等單元內(nèi)第j種指定作物第k個分等指標的分值。
2.2.4 計算分等單元自然質(zhì)量等指數(shù)
分別采用最高產(chǎn)量比法、營養(yǎng)當量修正法計算分等單元的自然質(zhì)量等指數(shù)
最高產(chǎn)量比法自然等指數(shù)的計算公式為
式中Rij為第i單元第j種指定作物的自然質(zhì)量等指數(shù);αij為第j種作物的光溫生產(chǎn)潛力指數(shù):CLij為分等單元指定作物的農(nóng)用地自然質(zhì)量分;i為分等單元編號,j為指定作物編號,k為分等因素編號;βij為第j種作物的產(chǎn)量比系數(shù):Ri為第i個分等單元的農(nóng)用地自然質(zhì)量等指數(shù);j為指定作物種數(shù)。
營養(yǎng)當量修正法的計算公式為
式中β新ij為修正后第j種作物的產(chǎn)量比系數(shù)。
按照《規(guī)程》確定曲周項目區(qū)農(nóng)用地分等工作中的基準作物為小麥,指定作物為玉米,標準耕作制度為一年兩熟,主要糧食作物是小麥、玉米。項目區(qū)基準作物小麥光溫生產(chǎn)潛力指數(shù)為1 096,玉米光溫生產(chǎn)潛力指數(shù)為2 700。按照最高產(chǎn)量比法求取小麥、玉米的產(chǎn)量比系數(shù)分別為1.0、0.7。按前文方法求取小麥、玉米的營養(yǎng)當量系數(shù)為 1.0、1.41。通過最高產(chǎn)量比法計算得出曲周項目區(qū)農(nóng)用地自然等指數(shù)分布在2 060~2 896之間,通過等間距法劃分等別序列為11~15等,共5個等別;通過營養(yǎng)當量修正法計算得出項目區(qū)的自然等指數(shù)分布在2 516~3580之間,通過等間距法劃分等別序列為13~18等,共6個等別,與最高產(chǎn)量比法相比最低等提高2等,最高等提高3等,等指數(shù)和等別整體呈上升趨勢(圖1、表6),等別差異擴大。
按照《規(guī)程》確定丁蜀鎮(zhèn)農(nóng)用地分等工作中的基準作物為水稻,指定作物為小麥,標準耕作制度為一年兩熟,主要糧食作物是水稻、小麥。水稻光溫生產(chǎn)潛力指數(shù)為1 977,小麥為930。按照最高產(chǎn)量比法求取水稻、小麥的產(chǎn)量比系數(shù)分別為1.0、1.3。按前文方法求取水稻、小麥的營養(yǎng)當量系數(shù)分別為 1.0、0.63。采用最高產(chǎn)量比法計算可得,丁蜀鎮(zhèn)農(nóng)用地自然等指數(shù)分布在 2 460~2 918之間,通過等間距法劃分等別序列為13~15等,共 3個等;通過營養(yǎng)當量修正法計算得出丁蜀鎮(zhèn)的自然等指數(shù)分布在2 013~2 393,通過等間距法劃分等別序列為11~12等,共2個等,與最高產(chǎn)量比法相比最低等降低2等,最高等降低3等(表6),等別整體呈下降趨勢(圖2、表6),等別差異縮小。
圖1 曲周項目區(qū)自然質(zhì)量等指數(shù)散點圖Fig.1 Scatter plot of natural quality index in Quzhou county project area
圖2 丁蜀鎮(zhèn)自然質(zhì)量等指數(shù)散點圖Fig.2 Scatter plot of natural quality index in Dingshu town
表6 不同產(chǎn)量比系數(shù)計算農(nóng)用地分等指數(shù)Table 6 Statistics of farmland classification index with different methods
此外,丁蜀鎮(zhèn)與曲周縣項目區(qū)相比,采用最高產(chǎn)量比法求得農(nóng)用地自然等,最高等均為15等,最低等高于項目區(qū) 2等;采用營養(yǎng)當量修正法求得農(nóng)用自然等,最高等低于項目區(qū)6等,最低等低于項目區(qū)2等,兩個案例區(qū)間的等別差異增大。
農(nóng)用地分等是依據(jù)構(gòu)成土地質(zhì)量穩(wěn)定的自然條件和經(jīng)濟條件,在全國范圍內(nèi)進行的農(nóng)用地質(zhì)量綜合評定[11]。由于各省自然、經(jīng)濟、技術條件不同,農(nóng)用地的地域性比較強,農(nóng)用地分等所采取的方法和所選技術參數(shù)不一致,同一評價方法往往很難適用于不同的區(qū)域,分等成果不具備直觀可比性[27-28]。產(chǎn)量比系數(shù)是農(nóng)用地分等中的一個基礎參數(shù),產(chǎn)量比系數(shù)的正確與否對農(nóng)用地分等起著關鍵作用,鑒于農(nóng)用地分等結(jié)果需要進行全國匯總,因此全國范圍內(nèi)應該采用統(tǒng)一的產(chǎn)量比系數(shù)計算方法[29]。如何選擇合理的產(chǎn)量比系數(shù),將各種不同的作物換算到同一基準上(可比的標準糧)進行比較或計算,保證農(nóng)用地分等成果的全國可比是一項重要內(nèi)容[30]。目前計算產(chǎn)量比系數(shù)的方法中最高產(chǎn)量比法以簡單易行而得到廣泛使用,但其確定的最高單產(chǎn)具有一定的區(qū)域性,所獲得的數(shù)據(jù)也受到時間和社會經(jīng)濟條件的限制,對于計算的精確度產(chǎn)生一定的影響;生產(chǎn)潛力比法通過作物的光溫(氣候)生產(chǎn)潛力計算產(chǎn)量比系數(shù),較適用于單一指標區(qū)內(nèi)的耕地生產(chǎn)能力計算,如同一個區(qū)域內(nèi)跨兩個或以上的指標區(qū)就存在選取何種生產(chǎn)潛力參加計算的問題;市場價格比法受市場影響強烈,不能保證時間上的持續(xù)性,穩(wěn)定性較差,精確度較低;熱量比法采用作物熱量這一的生物能量指標進行計算,不受因素指標區(qū)的限制,但選用這種方法進行產(chǎn)量比系數(shù)的計算,對外業(yè)調(diào)查的精度要求較高,數(shù)據(jù)獲取的難度較大,數(shù)據(jù)精度相對較低;生產(chǎn)潛力—熱量綜合比法,就是在耕地本底生產(chǎn)能力的基礎上,以作物的生物能量指標加以控制,剔除了人為因素的影響,提高了計算結(jié)果的精度,但是由于兩個指標對外業(yè)調(diào)查數(shù)據(jù)要求精度較高,適用范圍也有一定的限制,可操作性較差[16-17]。營養(yǎng)當量修正法中所選的蛋白質(zhì)、脂肪和碳水化合物不僅具有給人體供給熱量的功能,其分解產(chǎn)生的氨基酸、脂肪酸等物質(zhì)還具有構(gòu)成及修補人體組織、調(diào)節(jié)生理功能的作用,例如脂肪除為機體提供高深度的熱量和必要的熱量儲備外,還可以提供人體所必需脂肪酸,促進維生素的吸收。蛋白質(zhì)、脂肪和碳水化合物與其他所選營養(yǎng)素共同促進人體生長,保持代謝平衡,維持體內(nèi)環(huán)境穩(wěn)態(tài)。種植作物的主要目的是為了獲取其中的營養(yǎng)以滿足人體機能的需求,營養(yǎng)當量修正法不僅涵蓋了熱量,還包含了維持人類人體生理機能所必須營養(yǎng)素,與熱量法相比,其數(shù)據(jù)獲取的難度小、精度高,且更為完整、全面?,F(xiàn)有農(nóng)用地分等工作主要是基于耕地的自然屬性、利用水平等相關因素,從產(chǎn)出的角度看,盡管考慮到產(chǎn)量差異的問題,但是對于不同作物對人體所需營養(yǎng)素的貢獻程度,也就是土地產(chǎn)出作物的營養(yǎng)價值沒有被科學考慮,即使不同地域具有不同的標準耕作制度,基準作物、指定作物也有所不同,但不同作物供給人體所需能量和營養(yǎng)是統(tǒng)一可比的,不存在區(qū)域性的影響,那么用營養(yǎng)當量來修正產(chǎn)量比系數(shù)較現(xiàn)有的產(chǎn)量比系數(shù)計算方法更具科學性、準確性和實用性。因此在《規(guī)程》的框架下,可通過構(gòu)建不同作物的“營養(yǎng)當量”來修正現(xiàn)有的產(chǎn)量比系數(shù),建立更為規(guī)范的等級體系,服務于農(nóng)用土地開發(fā)整理設計、耕地占補平衡[31]。
與最高產(chǎn)量法相比,曲周縣項目區(qū)營養(yǎng)當量修正法計算得出的農(nóng)用地自然等等別提高,而丁蜀鎮(zhèn)通過營養(yǎng)當量修正法計算得出的自然等等別降低,由于營養(yǎng)當量修正法是將產(chǎn)量比系數(shù)由單位面積的最高產(chǎn)量比修正為單位面積生產(chǎn)作物所產(chǎn)出的營養(yǎng)比,曲周項目區(qū)為小麥玉米一年兩熟制,其基準作物為小麥,丁蜀鎮(zhèn)為小麥水稻一年兩熟制,其基準作物為水稻,而水稻的脂肪、膳食纖維、礦物質(zhì)、維生素的含量均低于小麥,營養(yǎng)素含量的差異是造成 2個案例區(qū)等別差異變化的主要原因。但根據(jù)“營養(yǎng)當量”計算產(chǎn)量比系數(shù),僅適用于可食用的作物,對于標準耕作制度為棉花等不可食用作物的區(qū)域,該方法則不適用。因此在后續(xù)的研究中,可將“營養(yǎng)當量”這一方法與其他計算方法相結(jié)合,形成合理有效的產(chǎn)量比系數(shù)綜合計算方法,并選擇研究區(qū)進行實證研究,為查清資源質(zhì)量,實現(xiàn)國土資源從數(shù)量管護到數(shù)量質(zhì)量管護并重,耕地的異地占補平衡提供科學依據(jù)。
本文通過構(gòu)建不同作物的“營養(yǎng)當量”來計算產(chǎn)量比系數(shù),選擇以小麥為基準作物的河北省曲周縣土地整治項目區(qū)和以水稻為基準作物的江蘇省宜興市丁蜀鎮(zhèn)南北方兩個項目區(qū)檢驗不同方法對農(nóng)用地自然質(zhì)量等等別的影響,研究結(jié)果表明:
1)營養(yǎng)當量修正法不存在區(qū)域性影響,較現(xiàn)有產(chǎn)量比系數(shù)計算方法更具科學性;
2)與最高產(chǎn)量比法相比,營養(yǎng)當量修正法得出的自然質(zhì)量等在曲周項目區(qū)呈現(xiàn)整體上升趨勢,由11~15等提升為13~18等;在丁蜀鎮(zhèn)呈現(xiàn)整體下降趨勢,由13~15等降為11~12等;
3)營養(yǎng)當量修正法增加了位于北方的曲周項目區(qū)和位于南方的丁蜀鎮(zhèn)兩個研究區(qū)間自然質(zhì)量等的等別差異。