1. 引言
自然资源资产核算是推动自然资源保护和生态文明建设的基础性工作 [1] 。2019年4月,中办、国办印发《关于统筹推进自然资源资产产权制度改革的指导意见》,提出“研究建立自然资源资产核算评价制度”,并进一步指出要“完善自然资源资产分等定级价格评估制度和资产审核制度” [2] 。耕地资源作为自然资源资产的重要组成部分,其数量和质量的长期稳定,对于维护国家粮食安全、保护自然生态环境和统筹城乡发展有着其他资源无法替代的作用 [3] [4] 。同时,耕地是粮源之基,在确保耕地面积稳定情况下,提高耕地产能是确保粮食安全的根本 [5] 。因此,研究耕地产能对于摸清粮食安全家底、制定针对性耕地保护政策,确保国家粮食安全具有重要意义 [6] 。国内外对于耕地产能核算的研究总体上可以分为以下3个阶段:土地自然生产力研究阶段、光温(气候)生产力研究阶段、耕地健康产能研究阶段。常用的耕地产能核算方法主要有趋势外推法、机制法、作物生长模拟法、经验法、Storie指数法、The Muencheberg Soil Quality Rating (M-SQR)法等 [7] 。
本文通过建立农用地分等单元标准粮理论单产模型与可实现单产模型,从而获得分等单元的理论单产及可实现单产,以此为基础核算区域理论产能与可实现产能 [8] 。以陕北长城沿线风沙区为例,基于农用地的自然质量等指数与利用等指数,建立产能核算模型 [9] ,结合实际调查资料,通过GIS手段来分析研究指标区的理论、可实现和实际产能空间分布状况,并为该地区提高农用地产能提出相关建议,同时,本文构建的数学模型对其它地区的耕地产能核算也有一定的参考价值。
2. 研究方法
2.1. 理论产能
理论产能是在农业生产条件得到充分保证,光、热、水、土等环境因素均处于最优状态,技术因素所决定的农作物所能达到的最高产量 [10] 。主要研究思路:以指标区为研究对象,基于分等单元建立农用地标准粮理论单产与相应的农用地自然质量等指数的函数关系,再利用函数方程获取分等单元的农用地理论单产与理论产能。
2.1.1. 理论产能模型
以二级指标区为单位,建立指定作物理论单产样本值与样本地块相应的自然质量等指数的函数关系模型:
(1)
式(1)中Yi为第i个单元标准粮理论单产样本值;Xi为第i个分等单元自然质量等指数;a、b为回归系数值。
2.1.2. 理论单产与理论产能
依据a和b,将该指标区内所有农用地分等单元的自然质量等指数代入公式(1)中获得各单元的年均标准粮理论单产(YF)。
分等单元理论单产乘以分等单元耕地面积获得农用地分等单元理论产能。公式如下:
(2)
式(2)中WFi为第i个分等单元理论产能,YFi为第i个分等单元理论单产,Si为第i个分等单元耕地面积。
2.2. 可实现产能
可实现产能核算指在农业生产条件得到基本保证,光、热、水、土等环境因素均处于正常状态、技术条件可以满足,由政策、投入等因素决定的正常年景下农作物能够获得的最高产量 [11] 。主要研究思路:以指标区为研究对象,建立抽样单元的可实现单产与相应的农用地利用等指数的函数关系,再利用函数获取各分等单元的农用地可实现单产与可实现产能。
2.2.1. 可实现产能模型
以二级指标区为单位,建立指定作物可实现单产样本值与样本地块相应的自然利用等指数的函数关系模型:
(3)
式(3)中Yi为第i个单元标准粮可实现单产样本值;Xi为第i个分等单元利用量等指数;c、d为回归系数值。
2.2.2. 可实现单产与可实现产能
依据c和d,将该指标区内所有农用地分等单元的利用等指数代入公式(3)中获得各单元的年均标准粮可实现单产(YF)。
分等单元理论单产乘以分等单元耕地面积获得农用地分等单元可实现产能。公式如下:
(4)
式(4)中WFi为第i个分等单元可实现产能,YFi为第i个分等单元可实现单产,Si为第i个分等单元耕地面积。
2.3. 实际产能
根据2006~2008年各县(区)统计年鉴和农业统计年报等资料,将各县(区)乡(镇)级的指定作物单位面积产量视作乡(镇)指定作物的实际单产;然后通过标准粮换算系数获得标准粮实际单产;最后基于实际耕作制度核算实际产能。
2.4. 利用强度与利用潜力
理论利用强度等于可实现单产与理论单产的比值,可实现利用强度等于实际单产与可实现单产的比值;理论利用潜力等于理论单产与可实现单产的差值,可实现利用潜力等于可实现单产与实际单产的差值。
3. 研究区模型的建立
3.1. 研究区概况
陕北长城沿线风沙区主要包括府谷县、神木县、榆阳区、横山县、靖边县和定边县,总土地面积33716.41平方公里,占全省土地总面积的16.38%。本区域地貌类型主要是由沙丘、河谷、湖盆滩地及盖沙黄土丘陵组成的,属温带半干旱大陆性季风气候,具有寒冷与半干旱的特征。陕北长城沿线风沙区耕作制度为一年一熟,指定作物为春玉米和马铃薯,共确定了72个理论样点,670个可实现样点。
3.2. 理论产能模型的建立
根据本指标区的耕作制度,分别建立春玉米和马铃薯的自然质量等指数和标准粮产量之间的关系模型,模型如下:
3.2.1. 春玉米理论产能模型
以春玉米调查样点的理论单产为纵轴,其各自对应的自然质量等指数为横轴,建立线性回归模型。如图1所示。
(5)
式中:y为理论单产;x为自然质量等指数。
3.2.2. 马铃薯理论产能模型
以马铃薯换调查样点理论单产为纵轴,其各自对应的自然质量等指数为横轴,建立线性回归模型。
(6)
式中:y为理论单产;x为自然质量等指数。
![](//html.hanspub.org/file/8-2420646x14_hanspub.png?20240304082455251)
Figure 1. Specify the relationship between the theoretical yield of the standard grain of the crop and the index of natural quality
图1. 指定作物标准粮理论单产与自然质量等指数的关系
3.3. 可实现产能模型的建立
根据本指标区的耕作制度,分别建立春玉米和马铃薯的利用等指数和标准粮产量之间的关系模型,模型如下:
3.3.1. 春玉米可实现产能模型
以春玉米调查样点可实现单产为纵轴,其各自对应的利用等指数为横轴,建立线性回归模型。
(7)
式中:y——可实现单产;x——利用等指数。
3.3.2. 马铃薯可实现产能模型
以马铃薯调查样点可实现单产为纵轴,其各自对应的利用等指数为横轴,建立线性回归模型。如图2所示。
(8)
式中:y——可实现单产;x——利用等指数。
![](//html.hanspub.org/file/8-2420646x17_hanspub.png?20240304082455251)
Figure 2. The relationship between yield and utilization index can be realized by specifying crop standard grains
图2. 指定作物标准粮可实现单产与利用等指数的关系
3.4. 标准粮单产核算
通过产能核算模型计算获得春玉米与马铃薯各分等单元的单产;再分别将两种作物的单产乘以标准粮转换系数,得到分等单元分作物标准粮理论单产;考虑到该县属陕北长城沿线风沙区,春玉米与马铃薯的种植制度为一年一熟,再将玉米与马铃薯的标准粮理论单产取平均值作为单元标准粮理论单产。
(9)
其中:
为第i个分等单元的标准粮理论单产;yij为第i个单元第j种作物标准粮理论单产。
3.5. 模型检验
对本指标区理论单产模型和可实现单产模型进行显著性检验,第一自由度f1 = 1,第二自由度f2 = 样本数-2,F检验结果见表1。
df为自由度,SS为平方和,MS为均方,查F表得
,由表1可知,以上模型检验结果中的F值
,由此可得该标准粮产量与分等指数之间有极显著的线性相关,其置信度均在95%以上。
![](Images/Table_Tmp.jpg)
Table 1. Linear regression analysis F-test
表1. 线性回归分析F检验
4. 结果分析
4.1. 农用地单产分析
4.1.1. 理论单产分析
陕北长城沿线风沙区的理论单产最大值为6687千克/公顷,最小值为4360千克/公顷,理论单产分布如图3。陕北长城沿线风沙区理论单产分布图如图3所示。从空间分布来看,陕北长城沿线风沙区地貌特征是理论单产呈现北高南低的趋势。从县域来看,各个县理论单产高的区域主要分布于地势平坦的川道以及沿河两岸的河谷阶地上。
![](//html.hanspub.org/file/8-2420646x22_hanspub.png?20240304082455251)
Figure 3. Theoretical yield distribution map of aeolian sand areas along the Great Wall of northern Shaanxi
图3. 陕北长城沿线风沙区理论单产分布图
4.1.2. 可实现单产分析
陕北长城沿线风沙区的可实现单产最大值为4265千克/公顷,最小值为3080千克/公顷,可实现单产分布如图4所示。从空间分布来看,可实现单产空间分布与理论单产分布一致,呈现北高南低的趋势。从县域来看,各个县的可实现单产高的区域主要分布于交通便利,农业设施较为齐全的地区,极大的调动了农民的生产积极性,因此土地的利用情况比较充分,产量稳定。
![](//html.hanspub.org/file/8-2420646x23_hanspub.png?20240304082455251)
Figure 4. Distribution map of achievable yields in aeolian sand areas along the Great Wall of northern Shaanxi
图4. 陕北长城沿线风沙区可实现单产分布图
4.1.3. 实际单产分析
陕北长城沿线风沙区的实际单产最大值为2100千克/公顷,最小值为1201千克/公顷。从空间分布来看,实际单产与理论单产、可实现单产空间分布基本一致,但空间分布趋于零碎。从县域分布来看,与县内各乡镇的微地形地貌、气候、土壤、灌溉、农业投入、社会经济发展状况等有关,从图5可以看出,指标区内各乡镇实际单产分布不均,整体上看北部比南部偏高。
![](//html.hanspub.org/file/8-2420646x24_hanspub.png?20240304082455251)
Figure 5. Actual yield distribution map of aeolian sand areas along the Great Wall of northern Shaanxi
图5. 陕北长城沿线风沙区实际单产分布图
4.2. 农用地产能分析
通过产能核算模型计算,从而获得农用地的单产和产能,结果如表2所示。
![](Images/Table_Tmp.jpg)
Table 2. Yield and production capacity results of wind-blown sand areas along the Great Wall of northern Shaanxi (Unit: hectares, ten thousand tons, kg/hectare)
表2. 陕北长城沿线风沙区的单产与产能结果表(单位:公顷、万吨、千克/公顷)
陕北长城沿线风沙区的理论产能为295.33万吨,可实现产能为215.61万吨,实际产能为94.26万吨。六个县区的理论产能在29.39~83.41万吨之间,可实现产能在22.26~60.28万吨之间,实际产能在9.73~27.09万吨之间。实际产能与可实现产能之间差值大于理论产能与可实现产能之间的差值,表明实际人为因素和基础设施完善情况是对农用地产能影响很大。
4.3. 利用强度与利用潜力分析
4.3.1. 利用强度与利用潜力分析
从表3可知,陕北长城沿线风沙区的理论利用强度0.73,可实现利用强度0.44。本区的理论强度在0.71~0.76之间,可实现强度在0.42~0.45之间。从六个县区的利用强度表明以下趋势:耕地紧张、交通顺畅、区位优越、土地肥沃、水资源丰富等自然、社会经济条件优越的区域,农用地理论产能利用强度偏高,反之偏低。因此,各县区应重视农业生产,加大农业基础设施建设,提高粮食综合生产能力。
![](Images/Table_Tmp.jpg)
Table 3. Utilization intensity and utilization potential of aeolian sand area along the Great Wall of northern Shaanxi (Unit: hectares, kg/hectare)
表3. 陕北长城沿线风沙区利用强度与利用潜力表(单位:公顷、千克/公顷)
陕北长城沿线风沙区的理论利用潜力1239.26千克/公顷,可实现利用潜力为1886.46千克/公顷;六个区县的理论利用潜力在1051.53~1378.11千克/公顷之间,可实现利用潜力在1751.44~2121.12千克/公顷之间。农用地理论产能潜力纵向地反映可实现产能的提升程度,农用地可实现产能潜力表明指标区的农业发展不均衡。
4.3.2. 利用强度与利用潜力相关分析
根据利用强度和利用潜力的定义,可知利用强度和利用潜力呈现负相关,说明土地利用强度越高,其利用潜力越小,土地利用强度越小,其利用潜力越大。为了进一步分析农用地利用强度和利用潜力相互影响关系,计算公式为如下:
(10)
式中:y为评价区域理论产能利用潜力;x为评价区域理论产能利用强度。
(11)
式中:y为评价区域可实现产能利用潜力;x为评价区域可实现产能利用强度。
![](//html.hanspub.org/file/8-2420646x27_hanspub.png?20240304082455251)
Figure 6. Correlation between theoretical capacity utilization intensity and utilization potential
图6. 理论产能利用强度与利用潜力相关关系
![](//html.hanspub.org/file/8-2420646x28_hanspub.png?20240304082455251)
Figure 7. The correlation between capacity utilization intensity and utilization potential can be realized
图7. 可实现产能利用强度和利用潜力相关关系
如图6、图7,从计算模型可以看出,理论产能利用强度每提高1个百分点,理论产能利用潜力可以挖掘73.87千克/公顷;可实现产能利用强度每提高1个百分点,可实现产能利用潜力可以挖掘42.42千克/公顷。
5. 结论与建议
5.1. 结论
(1) 从空间分布来看,陕北长城沿线风沙区理论单产呈现北高南低的趋势。可实现单产最大值为4265千克/公顷,最小值为3080千克/公顷。可实现单产空间分布与理论单产分布一致,呈现北高南低的趋势;实际单产最大值为2100千克/公顷,最小值为1201千克/公顷。实际单产与理论单产、可实现单产空间分布基本一致,但空间分布趋于零碎。
(2) 陕北长城沿线风沙区的理论利用强度0.73,可实现利用强度0.44。本区的理论强度在0.71~0.76之间,可实现强度在0.42~0.45之间。
(3) 陕北长城沿线风沙区的理论利用潜力为1239.26千克/公顷,可实现利用潜力为1886.46千克/公顷。且利用强度和利用潜力呈现负相关,说明土地利用强度越高,其利用潜力越小,土地利用强度越小,其利用潜力越大,从计算模型可以看出,理论产能利用强度每提高1个百分点,理论产能利用潜力可以挖掘73.87千克/公顷。
5.2. 建议
5.2.1. 严格控制耕地减少趋势,保护基本农田不受侵占
坚决贯彻执行《中华人民共和国土地管理法》,严禁一切形式的非法占用耕地特别是基本农田,合理规划乡、村建设,遏制土地撂荒现象,制止耕地非农化、耕地非粮化,认真落实国家对农业的各项扶持政策,加大对农业基础设施的投入,加强农田水利基本建设,不断提高耕地质量和粮食产能。
5.2.2. 有效控制农资价格
参照“三补贴”的成功做法,实行“第四补贴”即农资补贴,让农民得到更多的实惠,让国家的“惠农”政策不折不扣得到充分贯彻。
5.2.3. 加强对农民的技术培训
按照科学方法耕种管理。农业技术部门应经常送技术下乡,积极引导农民采用良种良法,科学种田,合理配量农资,避免盲目使用化肥、种子、农药;大力推行精耕细作,努力实现农业生产经营由粗放化向精细化的转变,以达到最低的投入成本取得最大的经济效益。
5.2.4. 力争耕地占补平衡,确保粮食生产安全
由于该区土地利用较差,净耕地系数小,耕地质量差,但未利用土地面积较大有待今后开发整理。耕地增加主要靠未利用地开发、农用地和建设用地整理、工矿废弃地复垦等来实现,今后应加大土地整治力度,以此来实现本区的耕地占补平衡。
基金项目
第三次全国国土调查陕西省耕地资源质量分类调查与评价(211535230406);
陕西省园地林地定级与基准地价研究(211535230407)。
NOTES
*通讯作者。