耕地资源是农业生产的基础和维护粮食安全的重要保障，中央农村工作会议也明确提出了对现有耕地质量数量加强保护的基本要求^[1-2]。而占补平衡项目的实施，既缓解了现有建设用地指标不足带来的巨大压力，同时也通过新造耕地，补充了耕地数量、提高了耕地质量，保障了国家1.20×10⁸ hm²耕地红线不动摇战略决策^[3-4]。但根据现实生产实际来看，耕地数量增加已不能完全作为补充耕地的第一要务，对耕地质量的要求已变得越来越举足轻重，各地区出现了以旱地改水浇地、旱地改水田以及水浇地改水田进而提高耕地质量的举措^[5]。现阶段，耕地等别和粮食产能已成为综合反映和评价耕地质量和数量的重要因素^[6]，因此，如何提高新增耕地质量对提高粮食产能具有显著影响。而对占补平衡项目区开展耕地质量评价，有利于更好的指导占补平衡项目的开展与进行，为占补平衡项目的实施提供耕地质量参考依据。截止目前，已有较多研究从耕地数量^[7-8]、耕地质量^[9-11]以及耕地资源^[12-13]等方面对粮食产能的影响进行了报道，对科技投入、生产条件以及资源配置等^[14-16]对耕地质量的影响也进行了研究，同时也利用空间统计分析、加权平均等方法分析了区域耕地质量等别限制因素及其提升对策^[17-18]。新增耕地质量影响因素较多，基于粮食安全耕地质量影响因素对粮食产能的影响是目前耕地保护方面研究的热点和难点。因此，本文以典型的陇中黄土旱塬区甘肃省景泰县2015—2017年3个年耕地占补平衡项目为例，对新增耕地的质量等别和粮食产能进行计算，并利用计量经济学Logistic回归模型，从构成新增耕地质量的主要因素入手分析各因素对粮食产能的影响，为实现黄土旱塬区耕地保护和提高粮食安全保障提供理论依据。

1 研究区概况

研究区位于甘肃省白银市景泰县寺滩乡，地理坐标位于东经103°33′—104°43′，北纬36°43′—37°38′，该区地形复杂多样，地势西高东低，项目区为典型的侵蚀堆积黄土丘陵地貌。年平均气温10.2 ℃左右，无霜期141 d，年平均日照2 669.7 h，属温带大陆性干旱气候，海拔在1 485~1 598 m，年平均冻土深度为1 m，年平均降雨量163.7 mm，其中7—9月降雨量约占全年降雨量的70%以上。项目区内土壤以淡灰钙土为主，土层深厚，耕性良好，熟化层厚30—60 cm，为粉沙轻壤—中壤，土壤养分的平均含量：有机质1.09%，全氮0.06%，速效钾171.8 ppm，速效磷7.9 ppm，pH值7.5，土壤容重1.4 g/m³，田间持水率24%，pH值8.1左右，土壤肥沃，通透性良好，土质富含农作物生长的各种营养成分。研究区主要种植作物类型为春小麦，熟制为一年一熟。耕地类型以旱地和水浇地为主，其中旱地57.45%，水浇地42.55%，根据《景泰县2019年耕地质量分等更新项目》内容，研究区国家级耕地质量各等别集中在12—14等。新增耕地采用工程措施主要为将表层0—50 cm深度范围土壤进行表土剥离，在通过工程机械提高田块平整度，将原有条件较好的其他草地平整后进行表土回填，最后再通过旋耕施肥。项目区灌溉用水主要来源于引大入秦工程——景泰县戈壁农业调水工程，通过引大入秦工程东二干渠引水干管引入项目区，设计流量0.42 m³/s，年供水量为4.15×10⁶ m³，可以满足项目区灌溉需求。

2 数据来源与研究方法 2.1 数据来源

耕地质量等别评价是对各等别耕地的空间分布、面积以及各分等单元的自然、利用和经济属性等的综合评价，对明确土地整治项目区耕地质量的限制性因素以及项目实施后所要达到的耕地质量级别具有重要作用，因此耕地质量等别评价在土地整治项目中已得到广泛的应用，可作为土地利用总体规划的基础依据和资料。本研究所使用的主要数据来源包括：①景泰县耕地质量等别评定成果。景泰县2015—2017年度完成的耕地质量等别补充完善成果，包括图件、数据库、表册和报告等，由此获得土地利用系数、土地经济系数、作物产量比系数、光温生产力潜力指数、气候生产潜力指数等数据。②土地变更调查数据库。景泰县2017年更新评价的土地变更调查数据库。③土地利用现状数据。景泰县自然资源局提供的基于全国第二次土地调查土地利用调查结果的现状图、地形图、等高图、土壤类型图等。④土地利用总体规划数据来源于景泰县土地利用总体规划(2006—2020年)。⑤其他资料：包含项目立项、测量、规划设计、竣工、复核报告、土壤检测化验报告、水质监测报告、环境监测报告等相关资料。

2.2 制定作物分等因素及相关参数获取

依据《甘肃省农用地分等技术规定》，景泰县属中黄土高原区陇中青东丘陵农牧区，其标准耕作制度为春小麦，一年一熟。春小麦的光温生产潜力指数为1 581，气候生产潜力为262。景泰县以春小麦为基准作物，其指定作物产量比系数为1.00。有灌溉设施的耕地采用光温生产潜力，缺乏灌排保障设施的耕地采用气候生产潜力。依据《农用地质量分等推荐因素及其分级、分值和权重》，汇总出景泰县农用地分等补充成果中指定作物的限制因子指标值、权重、级别和自然质量分。选取地形坡度(权重为0.27)、地表岩石露头度(权重为0.05)、有效土层厚度(权重为0.16)、表层土壤质地(权重为0.15)、土壤有机质含量(权重为0.13)、土壤pH值(权重为0.06)以及灌溉保证率(权重为0.18)作为分等因素。

2.3 样品采集及分析方法

在白银市景泰县2015—2017年选取不同耕种年限的3个占补平衡项目〔分别记为Z₁(1 a)，Z₂(2 a)，Z₃(3 a)〕范围内进行土样选择，选择10个样地面积为20 m 20 m的研究小区，按对角线法设调查样方3个(即3次重复)，样方面积为2 m×2 m，在10块标准地的30个样地内，每个样地采用直径为10 cm的土钻在每个样方中心位置采样，取土样深度为0—30 cm，取1 000 g土样，用四分法将多余的土壤弃去。全氮(TN)采用半微量凯氏法测定；土壤有效磷采用紫外可见分光光度计测定；速效钾用乙酸铵提取—火焰光度计法；有机质用重铬酸钾氧化外加热法(检测仪器：油浴锅DU-30G)。同时对表层0—30 cm范围内的原状土样，并保存其原有结构状态进行土壤物理性质测定^[19-20]。土壤粒度及微团聚体采用粒度分析仪测定。

2.4 研究方法

① 获取2015—2017年度内新增耕地图斑，将最近一期耕地质量等别图与最新土地利用现状图的耕地图层进行叠加，提取3个典型新增耕地图斑，形成年度更新评价数据。②评价新增耕地的质量等别信息通过项目资料、确认新增耕地的质量等别情况及相关属性信息，再通过外业补充调查进行检验，按照《农用地质量分等规程(GB/T28407-2012)》和耕地质量等别补充完善工作确定的技术方法和参数，对新增耕地进行耕地质量等别评定。③对占补平衡项目新增耕地土壤理化特性进行分析，揭示不同占补平衡项目年限对耕地土壤物理性质和化学性质的影响。④对新增耕地粮食产能影响因素进行分析。基于新增耕地粮食产能角度来看，新增耕地面积、新增耕地等别、新增耕地基础设施(含道路条件、灌排设施等)、新增耕地单位面积投资、新增耕地年限等对新增耕地粮食产能具有较大意义。本文以有序多分类Logistic回归模型对新增耕地粮食产能进行回归分析。

2.4.1 新增耕地产能计算方法

新增耕地产能的计算公式为：

$ \begin{array}{l} 粮食产能=新增耕地面积\times \\ \left( {D - 新增耕地平均质量等别} \right) \times 15 \times 100\;{\rm{kg}} \end{array} $

(1)

式中：D为指产能计算常数，D≤16(当产能为0时，D=16)。

2.4.2 新增耕地粮食产能影响因素分析

基于新增耕地粮食产能角度来看，新增耕地面积、新增耕地等别、新增耕地基础设施(含道路条件、灌排设施等)、新增耕地单位面积投资、新增耕地年限等对新增耕地粮食产能具有较大意义。有序多分类Logistic回归模型通常适用于因变量有序且是多分类的，本文将粮食产能提升情况划分为3个等级，则该Logistic回归模型为3分类Logistic回归模型(见表 1)。

$ \ln = \frac{{p(y - 1)}}{{1 - p(y = 1)}} + \alpha + \sum\limits_{i = 1}^k {{\beta _i}} {X_i} $

(2)

式中：因变量Y有3个等级，低产能为“1”，中产能为“2”，高产能为“3”。自变量即影响因素有：新增耕地面积X₁，新增耕地等别X₂，新增耕地基础设施X₃，新增耕地单位面积投资X₄，新增耕地年限X₅。本研究回归测算通过SPSS统计软件完成。

2.5 数据统计分析

采用Excel分析软件处理试验数据，同时采用SPSS 11.5软件进行统计学分析，并对相关指标进行显著性分析，显著性水平为(p＜0.05)。

3 结果与分析 3.1 景泰县2015-2017年新增耕地等别计算结果

本文选取景泰县2015—2017年3个耕地占补平衡项目，整治总面积115.56 hm²，新增耕地105.98 hm²，其中水浇地97.33 hm²，旱地8.65 hm²。本文采取“因素法”对项目新增耕地等别进行评价，采取的评价方式严格按照《农用地质量分等规程》进行评价因素指标体系和评价系数沿用耕地质量等级补充完善成果工作中采用的参数。表 2为2015—2017年选取的3个景泰县占补平衡项目国家级自然等指数、利用等指数、经济等指数计算结果。整体来看，2015—2017年在景泰县实施占补平衡项目新增耕地平均等别范围在11—12等之间。根据景泰县2015—2017年度耕地质量等别年度更新评价成果，景泰县项目区周边耕地质量国家级自然等别11—12等，国家级利用等别11—12等，国家级经济等别12—13等，平均等别为11—12等，因此可以看出，通过占补平衡项目的实施可以将新增耕地提升1—2个等别，此外，根据白银市耕地占补平衡动态监管系统查询可知，本文选取的3个补充耕地项目被占用的耕地平均等别为12等，因此占补平衡项目的实施得到新增水浇地等别高于被占用耕地1个等别，旱地等别与被占用耕地等别相同。这主要是在相关工程措施的投入下，增加了有效耕地面积，且提高了田间水利配套设施和改善了道路条件，通过土壤旋耕施肥的进行，增加了土壤养分和改善了土壤理化特性，有利于作物产量的提高，由于耕地等别受耕地面积、有效土层厚度、土壤养分含量、灌溉保证率以及交通条件等因素影响，因此，通过占补平衡项目的实施耕地等别较未实施前显著提高。此外，可以看到的是，新增耕地地类对耕地等别具有较大的影响，水浇地面积占比与新增耕地等别呈显著的正相关关系。

3.2 景泰县2015-2017年新增耕地粮食产能计算结果

表 3为2015—2017年景泰县内3个耕地占补平衡项目粮食产能变化情况。由表 3可以看出，粮食产能受耕地面积和耕地等别影响较大。2015—2017年耕地占卜平衡项目实现的耕地面积分别为38.51，27.24以及40.23 hm²，粮食产能分别为284.55，202.25以及295.08 t，粮食产能随新增耕地面积的增大而增大，粮食产能与新增耕地平均等别呈显著的正相关关系。此外，水浇地项目的粮食产能显著高于旱地项目的粮食产能，2015—2017年3个项目，占补平衡项目实现的水浇地面积为97.33 hm²，对应实现的产能为729.98 t，而旱地面积为8.65 hm²，对应产能为51.90 t。可以看出，水浇地以占总新增耕地91.84%的面积，创造了占比为93.36%的总粮食产能，旱地以占总新增耕地8.16%的面积，创造了占比为6.64%的总粮食产能。

3.3 占补平衡项目对新增耕地土壤理化性质的影响

土壤物理性质是土壤最基本的性状和指标，也是土壤质量评价的重要参数。表 4反映了不同年限占补平衡项目新增耕地土壤物理性状统计情况。由表 4可以看出，新增耕地总孔隙度、大团聚体、水稳性团聚体和土壤团聚度随新增耕地年限的增长而增加，到新增耕地3 a时，涨幅分别为3.15%，1.51%，35.23%和23.96%，而＜1 μm微团聚体随新增耕地年限增长表现为显著的降低趋势，经过3 a的变化，降幅达8.65%。旱塬区地形复杂，土壤结构输松，土壤结构破坏率最高，随着新增耕地年限的增加，提高了保水保土的能力，部分物理特性得到有效改善，土壤抗侵蚀能力也得到一定提升。我国目前主要采用土壤结构系数或分散系数评价土壤结构性指标，结构系数是指参与团聚化的黏粒与黏粒总量的比值，结构系数越大，微团聚体的水稳性越强，形成团聚体的潜在能力越大。土壤结构系数在新增耕地年限时间不断增大在3 a时达到18.47%，土壤分散系数随增耕地年限时间逐渐减小，在3 a时为81.53%。从土壤结构破坏率来看，新增耕地年限越长，其破坏率则越低，新增耕地年限3 a时的结构破坏率较1 a时降低了2.05%。主要由于新增耕地年限的增加，在经过不断的土壤翻耕和培肥等基础上，土壤水肥逐渐得到积累，促进了土壤颗粒的团聚作用，土壤中水稳性团聚体也越来越稳定，土壤水肥气热条件也得到有效改善，土壤抗侵蚀能力得到显著提高。

3.4 占补平衡项目对新增耕地土壤养分的影响

土壤化学性质是土壤质量的重要表征，对土壤营养状况和作物生长潜力有重要影响。表 5反映了不同年限新增耕地土壤化学性状。由表 5可以看出，随新增耕地年限的增加，全氮、有机质含量、速效磷含量等均表现为在新增耕地年限3 a时达到最大，全氮、有机质及速效钾分别较新增耕地年限1 a时提高了10.67%，15.39%以及13.62%。速效钾含量表现为先增大后减小的趋势，在新增耕地年限2 a时达到最大为199.21 mg/kg。土壤化学性质的改变受自然环境因素和人为因素影响较大，全氮、有效磷和速效钾的变化范围分别在0.75~0.83 g/kg，13.58~15.67 mg/kg以及178.62~199.21 mg/kg之间。占补平衡项目的实施，降低了水土流失和土壤侵蚀现象，水肥可以在垂直方向入渗及运移，再加上耕作过程中进行的土地翻耕、有机肥和无机肥的施用、合理灌溉和田间管理，促进了土壤中营养元素的积累，氮磷钾元素和有机物含量显著增加，土壤化学性质也得到了有效改善。

3.5 实证检验

土地生产能力是自然条件与土地开发利用程度共同作用的结果，粮食产能高低也是新增耕地质量各因素综合作用的结果，本文采用Logistic回归模型对新增耕地面积、新增耕地单位面积投资、新增耕地年限、新增耕地等别以及新增耕地基础设施等对新增耕地粮食产能的影响进行模型检验，结果详见表 6。由表 6可以看出，新增耕地面积越大，创造的粮食产能也越高，二者间显著性水平为0.021，虽然新增耕地面积是影响粮食产能提高的重要因素，但需要看到的是单位耕地面积上粮食产能的增量是有限的，因此，要依据当地土地质量的自然、经济因素优化作物种植结构，切实做好耕地保护，同时，结合当地气候条件等，适宜的提作物高复种比例，也有利于粮食产能的提高。另外，由于农业从业人员数量的不断减少，出现了耕地撂荒和弃荒等消极管理等现象，也会造成新增耕地粮食产能的降低。因此要做好项目后期管护工作，落实主体责任并提高农业种植积极性，避免出现耕地撂荒和弃荒现象。新增耕地单位面积投资的高低，一方面体现了项目的实施难度，另一方面也反映了项目技术手段和工程措施投入的多少。一般而言，项目单位面积投资越大，新增耕地粮食产能也越高，从分析结果来看，单位面积投资对新增耕地粮食产能的影响系数为0.193，低于新增耕地面积，但其显著性达到0.042，高于新增耕地面积。新增耕地年限对粮食产能的影响程度为0.038，是5项影响因素中数值最低，这是由于经过必要的工程施工，虽然对新增耕地进行了土壤翻耕、土壤培肥等措施以熟化土壤，但由于熟化过程存在一定的滞后性，部分营养物质和有机质含量未能在第一年就进行较好的分布，因此表现为新增耕地年限对粮食产能的影响程度较低。

耕地等别受当地土壤、气候、水分条件以及社会经济条件等综合因素影响程度较大，从分析结果可以看出，新增耕地等别由14等到11等的回归系数不断增大，而Wald值逐渐降低，以耕地等别11等为参照，14级=1的优势比OR值为0.089，13级=2的优势比OR值为0.134，较14级高50.56%，12级=3的优势比OR值为0.266，较13级高98.51%，优势比OR值与耕地等别呈显著的正相关关系，也就是说，耕地等别越高，对应新增耕地的粮食产能也越大。另外可以看到的是，新增耕地道路和水利配套设施通达、完善与否对新增耕地粮食产能的提高具有十分显著的影响。

新增耕地基础设施与粮食产能产生负向影响，随着基础设施完善程度的提高，道路和排灌设施完善程度对新增耕地粮食产能的影响系数为呈正相关关系，优势比OR值也不断增大，说明道路和排灌设施越完善，新增耕地粮食产能越大。这主要是由于道路的通达和配套设施的完善，提供了农业生产的基础保障，有利于机械化耕种，增强了群众耕作的意愿，且极大的节约了生产实践成本，因此对粮食产能的提高具有显著的促进作用。项目实施前种植玉米，耕地单位面积产能为7.65 t/hm²，经过3~5 a种植后，在土壤熟化和肥料施用以及集约管理的基础上，耕种5 a，4 a和3 a的项目单位面积产能分别达到8.85，8.40，7.95 t/hm²，较未实施时粮食产能分别提高15.69%，9.80%和3.92%。根据分析结果来看，在5项对新增耕地粮食产能影响较为突出的因素排序为：新增耕地等别＞新增耕地基础设施＞新增耕地面积＞新增耕地单位面积投资＞新增耕地年限。

4 结论

(1) 占补平衡项目的实施将新增耕地等别提升了1~2个等别，且水浇地面积占比与新增耕地等别呈显著的正相关关系。粮食产能随新增耕地面积的增大而增大，粮食产能与新增耕地平均等别呈显著的正相关关系。水浇地以占总新增耕地91.84%的面积，创造了占比为93.36%的总粮食产能。

(2) 占补平衡项目对新增耕地土壤理化性质和养分积累具有显著的改善作用，新增耕地总孔隙度、大团聚体、水稳性团聚体和土壤团聚度随新增耕地年限的增长而增加，到新增耕地3 a时，涨幅分别为3.15%，1.51%，35.23%和23.96%。新增耕地年限越长，其破坏率则越低，新增耕地年限3 a时的结构破坏率较1 a时降低了2.05%，土壤抗侵蚀能力得到显著提高。新增耕地年限3 a时，土壤全氮、有机质及速效钾分别较新增耕地年限1 a时提高了10.67%，15.39%以及13.62%。

(3) 经过3~5 a种植后，耕种5 a，4 a和3 a的项目单位面积产能较未实施时粮食产能分别提高15.69%，9.80%和3.92%，对新增耕地粮食产能影响较为突出的因素排序为：新增耕地等别＞新增耕地基础设施＞新增耕地面积＞新增耕地单位面积投资＞新增耕地年限，因此要将新增耕地粮食产能提高的着力点放在耕地等别的提高。