本课题组分别在2002年、2012年和2017年1-3月在重庆市内12个植烟区县进行试验，将重庆现代烟草农业规划定义的基本种植单元设置为一个采样点，每个样点面积约为6.67 hm²，依据海拔、地形、母质类型、土壤类型等因素，采用“S”型采样方法，采集20个点作为一个混合样，采用四分法获得2 kg土样作为最终样品，带回实验室.在试验区域内采用GPS定位技术进行采样，其中2002，2012和2017年分别采集土壤样品2 100，2 777和3 800个左右，采集的土壤样点基本是按照现代烟草农业规划一个种植单元一个采样点，所以采集的土壤样品大部分(约90%)是重合的.为了反映烟田的真实肥力状态，土壤样品采集时间为烟草底肥施用前，样品土壤经过自然风干、研磨、过筛之后混匀保存，用于测定各项指标.

土壤样品均送由西南大学样品测试中心进行分析检测，测定项目为：pH值、土壤有机质、碱解氮、速效磷、速效钾、有效硼、交换性镁、有效锌和水溶性氯，测定方法参考《土壤农化分析》^[4].

为全面评价土壤养分肥力状况，我们引入烟田土壤综合肥力指标值(Integrated Fertility Index，IFI)的概念，按照该值给土壤综合肥力划分等级，计算公式如(1)：

式中：W_i和N_i分别为第i种养分指标的权重系数和隶属度值，求得的IFI值的范围在0~1之间，值越大，说明土壤的养分肥力状况越好^[5].

隶属度函数采用模糊统计得出，养分指标隶属度值的估算采用下述步骤.首先标准化处理各项肥力指标值，以避免或消除各指标数量级和单位差异产生的误差.两种类型的隶属度函数被本研究采用，分别为S型和抛物线型(图 1)，公式(2)为S型隶属度函数表达式，公式(3)为抛物线型隶属度函数表达式.

由于不同参评因子对土壤肥力状况贡献度不同，因此必需确定评价指标的权重，同时各个因子的标准化值需要使用权重进行修正.本研究应用偏相关分析方法在SPSS 20.0中计算出参评指标的相关系数矩阵，建立相关系数矩阵.从而明确不同肥力指标对土壤肥力的实际贡献率，以尽量避免其他因素带来的影响.

由相关系数矩阵求出其逆矩阵R^-1：

由逆矩阵中的相关元素计算偏相关系数r_ij：

计算各因子与其他因子偏相关系数的平均值(R_i)，该单项肥力指标表征的土壤肥力贡献率是指该平均值所占所有肥力指标偏相关系数平均值之和(∑R_i)的百分率，即土壤肥力因子的权重值^[5].

采用SPSS 20.0软件对试验数据进行统计分析，用Excel 2010软件制图.

以2017年采样数据进行分析，重庆市烟区土壤肥力描述性指标如表 1所示，一般认为烤烟最适宜pH值范围在5.5~7.0，重庆市烟田pH值的变化范围为3.77~8.62，有部分烟田的土壤pH值已经低于4，酸化严重，pH＜5.0的酸性土壤占29.0%，烟田土壤有机质平均质量分数为27.5 g/kg.重庆植烟土壤的碱解氮质量分数为156.8 mg/kg，质量分数范围为52.6~386.9 mg/kg，其变异程度很大，质量分数过高和过低同时存在.土壤中速效磷平均质量分数为41.8 mg/kg，质量分数范围为2.5~95.6 mg/kg，变异系数45.2%；土壤中速效钾平均质量分数为315.7 mg/kg，质量分数范围为19.0~589.3 mg/kg，可以看出烟田土壤速效磷和速效钾的变异程度较大.烟田土壤中有效锌质量分数为2.26 mg/kg，重庆植烟土壤有效硼的质量分数为0.25 mg/kg，有效硼的分布区间广泛，变异系数为61.7%，有效硼分布不均匀.从表 2可以看出，土壤pH值在15年后下降了0.43个单位，水溶性氯下降了3.0 mg/kg，其他指标均有不同程度的增加.

根据烤烟研究资料和重庆烟田实际情况，并结合专家建议，综合确定各隶属函数的下限、下优、上优和上限，将这些所测得的指标带入对应的表达式(2)或(3)，计算出每个采样点的隶属度值，该值表示该指标在土壤中的质量分数状态，隶属度小的指标可能是烟田的障碍因子，需要针对性改良，隶属度值越高表明越适宜烤烟生长(表 3).

首先计算某肥力因子i与其他因子偏相关系数的平均值A_i，再计算出A_i占所有偏相关系数平均值之和的百分比，以表示土壤各肥力指标对土壤肥力的贡献率，即得出各项肥力指标的权重值(W_i).结果如表 4所示，各权重系数总和为1.由表 4可见，有机质的权重值达到了0.180 0，表明有机质是提高土壤肥力的关键，其次是pH值，权重为0.161 8，碱解氮、速效磷、交换性镁、有效硼依次递减，有效锌对土壤肥力的贡献率最小.

将土壤肥力综合指标值分为5个等级，即Ⅰ级的IFI≥0.8为高，Ⅱ级的0.6≤IFI＜0.8为较高，Ⅲ级0.4≤IFI＜0.6时为中等，Ⅳ级0.2≤IFI＜0.4时较低，Ⅴ级IFI＜0.2时肥力低.由表 5可以看出，重庆市植烟土壤综合肥力指数平均值是0.69，说明重庆市土壤肥力为中等较高水平，其变异系数17.71%，说明烟田土壤养分肥力水平分布不均匀.土壤肥力综合指数较高占63.53%，极高的占16.63%，可见重庆市的植烟土壤养分总体上处于丰富水平. 2012年和2017年综合肥力指数均为高(表 6)，而2002年综合肥力指数为中等.

养分肥力综合指数如表 7所示，不同植烟区县植烟土壤养分肥力水平差异有统计学意义，12个植烟县中以彭水县、武隆区、丰都县、涪陵区和巫溪县的植烟土壤养分肥力最高，IFI分别为0.77，0.76，0.76，0.70和0.69，土壤肥力适宜(高和极高)的土壤占80%以上.万州区土壤养分综合指数IFI最低，为0.53，仅56.4%土壤属较高肥力水平.

烟区土壤肥力的综合评价可为烤烟精准施肥和烟草品种改良提供依据.土壤肥力综合评价指标体系，主要考虑对烟草生长影响较大的土壤因素.本课题组以往研究证实，土壤综合肥力指数IFI与烟叶产量、产值、均价、上等烟产量的正相关关系均极有统计学意义，模糊综合评价法得出的重庆烟田土壤肥力指标具有重要的现实意义，对于植烟土壤肥力状况的评价是可行的^[5].土壤肥力综合指数评价结果可为重庆市烟草可持续发展提供可靠的理论及实践依据.比如在土壤综合肥力较高的地区有针对性采用降低肥料用量措施和引种耐肥性较高的品种；对土壤综合肥力较差的地区采用有机肥增施培土，适当增加肥料用量改善土壤养分状况，引进推广肥料吸收能力较好的烟草品种.

重庆市土壤肥力总体为较高水平，但烟田还存在一系列问题，比如变异系数17.71%，说明烟田土壤养分肥力水平分布不均匀.植烟土壤常年连作15年pH值下降了0.43个单位，水溶性氯下降了3.0 mg/kg，其他指标均有不同程度的增加.张福锁团队成员在《Science》上报道了我国主要农田出现显著酸化现象，农田耕作20年，pH值平均下降了约0.5个单位，重庆市植烟土壤酸化面积大，pH＜5.0的强酸性土壤占29.0%，在秦巴山区如黔江、石柱等区县强酸性烟田所占的比例较大.烟田土壤酸化的原因已基本明确，烤烟每年吸收带走的盐基阳离子(K⁺，Ca²⁺和Mg²⁺等)没有得到有效的归还，氮肥滥用也促使了土壤酸化的加剧，这为防治烟田进一步酸化提供了坚实的理论基础.烟杆和烟叶每年带走了土壤中大量的钙镁，却没有及时得到归还，致使烟田土壤中交换性钙镁离子缺乏. 表 2可见，2017年重庆烟区土壤重交换性钙离子质量分数比10年前减少了4.32 cmol/kg，且整个重庆市植烟区交换性钙质量分数缺乏^[6].防治烟田酸化本项目组采用分类治理的策略，即对于pH值在5.0~5.5的烟田，推荐施用钙镁磷肥、硅钙肥和生物炭等弱碱性物料，同时可以归还烟田钙和镁元素来改善烟田土壤酸性.对于pH＜5.0的烟田，推荐采用生石灰进行改良，具体是在整地前撒施750 kg/hm²，起垄前30天撒施750 kg/hm²，烟草移栽期60天翻压.

重庆烟区土壤90%以上有机质质量分数超过15 g/kg，与国外一些优质烟叶生产国相比，其质量分数偏高，如美国北卡州烟区有机质低于15 g/kg的土壤占80%以上^[7].土壤有机质与土壤供氮能力呈正相关，重庆烟区土壤有机质矿化后能为烟草提供大量氮素，特别是烟草生长后期高温高湿的环境条件下，大量供应氮素使烟草贪青晚熟，上中部烟叶烟碱质量分数提高，大幅度降低烟叶吸味品质；相反，有机质质量分数过低则使烤烟香气量不足，因此适宜重庆烟田的土壤有机质质量分数范围还有待进一步研究.重庆植烟区县缺磷和磷丰富的土壤并存，磷主要参与植物新陈代谢过程，促进部分重要化合物的合成，增强烟株抗逆能力，提高烟叶品质^[8].重庆市万州和酉阳土壤磷质量分数较低，对于该地区土壤应注重增施磷肥，以满足烟草正常新陈代谢和品质形成的需要.重庆山地烟区多雨，钙镁淋失也非常严重，而镁的补充往往被农户所忽视，导致重庆植烟土壤镁质量分数严重缺乏.严重缺钙镁的土壤主要为酸性土壤，也是土传病害易发区，建议施用土壤调理剂(富含活性钙镁如硅钙肥、钙镁磷肥等)或增加生物有机肥施用量，在补充钙镁元素的同时，对土壤生物肥力进行修复和重建.在烟田大面积生产中，考虑在提苗肥或追肥时施用硝酸钙、硫酸钾镁肥等，可以协调平衡烟草氮、钾、钙和镁营养；烟草是需硼较多的作物，随着氮、磷和钾大量元素的大量施用，土壤有效硼不足已成为限制烤烟产质量的重要因素之一，硼元素主要生理功能是参与糖的转运与代谢，缺硼对烟株主根和侧根生长影响最大，由于根系受阻，烟叶中烟碱质量分数和糖分合成均受阻碍^[9].从土壤有效硼来看，重庆市烟田土壤有效硼质量分数较10年前显著降低，所以重庆市烟区更应该注重硼肥的施用，通过土壤施用硼砂等补充硼(注意不能连年用)，应根据不同的土壤特性及养分情况选择相应的措施.硼肥施用在酸性砂质土壤中易淋失而降低肥效，在黏土中肥效可持续数年.推荐3.75~7.5 kg/hm²作基肥，或0.2%在团棵期进行叶面喷施，或在烤烟专用复合肥中加0.5%硼砂.重庆市烟田土壤肥力总体状况较好，根据烟田实际情况调整烟草专用肥配方，严格控制氮肥用量，因地制宜调整磷钾配比，微量元素因缺补缺，构建基于基地单元土壤肥力的“大配方”施肥方案.根据不同区县或种植单元土壤养分质量分数调整氮磷钾比例，实现种植单元“小调整”.在肥料配方上多样化，设置高磷、中磷、低磷以及高钾、中钾、低钾配方，改变目前生产上1个配方的局面，增加配方肥种类，做到因地制宜施用肥料，实现按区域平衡施肥^[5].

氯对于烟草品质具有重要作用，烟草需要适当施氯，氯能够增加烟叶的弹性和油润性^[10].植烟土壤水溶性氯一直呈现下降趋势，与15年前相比，土壤水溶性氯下降了3.0 mg/kg单位，其变异系数为43.6%，属于强变异性.结合重庆地区土壤现状，烟草生长所需水溶性氯最佳质量分数在10~30 mg/kg之间^[11].土壤中水溶性氯易随水分流失，重庆市水溶性氯小于10 mg/kg的土壤所占比例为63.5%，水溶性氯偏低的植烟区县有万州、巫山、巫溪和奉节等区县，这些地区急需要补充氯.由于烟草对氯特别敏感，为防止氯化钾的施用不当，建议在烟草配方肥中每50 kg配入3~5 kg氯化钾为宜，这样可以满足烟草对氯的需要，又可防止氯化钾施用不当对烟草造成不良影响^[12].

烟田土壤肥力主要是指烟田土壤中各种养分的供应能力，是土壤物理、化学和生物因子作用于土壤的综合体现^[13-14].本研究仅采用单一的统计方法对重庆市植烟土壤土壤肥力进行了综合评价，尚存在不完善的地方，比如是否需要考虑生物指标，是土壤肥力工作者一直思考的问题，由于生物指标受环境的影响较大，到目前为止，还未见在综合评价时考虑土壤生物指标的报道，因此今后应加强综合、全面评价土壤肥力指标方面的研究.土壤肥力评价还应该将当地气候、栽培水平等因素加以综合考量，以便更好地为烤烟综合养分管理提供科学参考依据.

1) 重庆市植烟土壤肥力因子变异强度较大，多处于中等至强的变异强度，土壤速效磷、速效钾和有效锌变异强度较大.

2) 对土壤综合肥力指数贡献最大的为有机质，其次是碱解氮、pH值、有效硼、速效磷，贡献最低的是有效锌.

3) 重庆市植烟区县整体土壤肥力较为适宜，不同区县对于土壤肥力水平较低的烟田需加强测土配方施肥，采用综合养分管理技术，进一步提升烟田地力水平.