摘要：为了探明水肥配施对枸杞生长、产量和品质等方面的综合影响，以获取最适水肥配施方案，设置3个不同灌溉定额(194，254，313mm)和3个不同N-P2O5-K2O配比施肥量(135-45-90，180-60-120，225-75-90kg/hm2)9个组合，采用完全随机设计，利用AHP和熵权法相结合的综合评判方法对枸杞的4类一级指标和12类二级指标进行综合评判，以构建枸杞水肥配施综合评判模型.通过最小二乘法对AHP法和熵权法进行优化组合，计算综合评判权重，以自然间断法将综合评判权重进行分类，确定水肥配施方案的决策依据.以综合评判指数为因变量，以灌溉定额和总施肥量为自变量进行逐步二项式回归分析，建立回归方程，分析水肥耦合下对综合评价指数的影响，并确定在综合评价指数最大时对应的水肥配施方案.结果表明，综合评判指数随水肥施量增加呈先升高后降低的趋势，与生产实际相符;各类因素对综合评判指数的影响中，粒度、百粒重、干果产量和总糖这4个因素对枸杞生长的综合评判影响较大，占总评判指数73.3%以上.基于AHP和熵权法组合优化后得出的枸杞最优水肥方案：灌溉定额为252～278mm(不含春冬灌)，配比施肥量(N-P2O5-K2O)为197-66-131kg/hm2～203-68-135kg/hm2，枸杞综合效益最佳.
关键词：水肥方案;枸杞;综合评判;层次分析法;熵权法
宁夏是枸杞的原产地，已有500多a的栽培历史，果实、根皮及嫩叶均可入中药.截至2019年底，宁夏的枸杞种植面积超过2.33万hm2，占全国枸杞种植面积一半以上，是当地助力脱贫的重要经济作物.作物产量和品质不仅需要水分和肥料这2个关键因素，气候环境、温度和降水等因素的影响也尤为重要[1].因此，如何制定出科学、合理的水肥方案是水资源高效利用前提下，作物高产需要面临的重要问题.目前，对枸杞水肥方案制定的研究工作多集中在利用产量、品质或水肥利用效率等单方面的效应为依据进行[2-3].而利用多指标综合评判的方法制定灌溉施肥方案，目前已广泛应用在温室番茄[4]、小麦[5]、水稻[6]等作物上
.常用的综合评判方法有主成分分析法[7]、层次分析法(analytichierarchyprocess，AHP)[8]、模糊综合决策[9]、熵权法[10]等.AHP法权重的确定是由专家根据1～9的标度法进行打分，主要依赖于专家经验;熵权法则是根据各评判指标的实测数据确定其权重的客观评判方法，明显的缺点是过度依赖实测数据，没有考虑专家经验.AHP和熵权法相结合的评判方法已成功应用于水环境安全、水资源[11]、富水性预测[12]等领域.
综上所述，文中采用基于AHP和熵权法相结合的赋权方法对枸杞从气候、生长、产量和品质这4个方面进行综合评判，同时分析各项指标对不同水肥方案的响应，旨在为枸杞优质高产水肥方案优化提供借鉴.
1材料与方法
1.1试验基地概况
试验于2018年和2019年4—10月在宁夏中部干旱区同心县润德枸杞庄园(105°54′E，36°58′N，海拔1240m)试验基地开展.该地区属典型大陆季风型气候，多年平均降水量为270mm，且多集中在7—9月，干旱指数为8.4，地下水埋深大于30m.距离试验区域5m空旷处安装Decagon微型气象监测站，监测研究区域气象情况，2018年和2019年生育期内降水总量分别为172.9和175.5mm，≥5mm的有效降雨量分别为151.65，154.90mm.试验基地2018—2019年枸杞全生育期内气候变化如图1所示，图中物理量为日平均气温tda、降雨量P、参考作物蒸发蒸腾量ET0;土壤理化性质见表1，表中ρ为体积质量.
试验基地供试作物品种为4a生“宁杞7号”枸杞，长势均匀，株行距为75cm×300cm，种植密度为4440株/hm2.采用滴灌方式，铺设壁厚0.15mm、内径16mm的内嵌式滴灌管，滴头流量为3.0L/h，滴头间距30cm，每棵树左右各有滴头1个;滴灌管距枸杞树主干5cm，灌水量由水表控制.春、冬水分别为300和450m3/hm2.
1.2试验设计
田间试验设计采用完全随机试验设计方法，设定灌溉定额(W)和施肥量(F)2个因素，各3个水平，1个灌水不施肥梯度作为对照(CK)，共10个处理.每个处理重复3次，每个小区10株枸杞树.各小区灌水次数相同，灌水7次，萌芽期、春梢生长期和开花期各灌水1次，果熟期4次.每次灌水定额相同，均为灌溉定额的1/7.灌溉定额的3个水平W1，W2，W3分别为194，254，313mm;配比施肥量(N-P2O5-K2O)的3个水平F1，F2，F3分别为135-45-90，180-60-120，225-75-150kg/hm2.施肥方式采用水肥一体化，肥料选用尿素(含氮量46%)、过磷酸钙(含磷量12%)、硫酸钾(含钾量50%);全生育期施肥共7次，占总施肥比例为春梢生长期20%(1次)、开花初期20%(1次)、果熟期50%(分4次等量施用)、落叶期10%(1次).
1.3试验测定项目及评判说明
1.3.1气候指标
枸杞全生育期气候指标主要考虑界限温度为≥10℃的积温和灌溉条件下的有效降水量[13].采用Decagon微型气象站监测试验区域气温、降水量、蒸发量和风速等气象参数，数据步长为1h.根据前人[13]研究发现，枸杞全生育期≥10℃积温为3200～3600℃，则枸杞能够获得正常生长所需的热量，此时热量将不受枸杞产量和品质的限制.在灌溉条件下，全生育期有效降水量为100～170mm，则枸杞产量不受降水量的影响;降水量小于100mm时，对枸杞的产量产生不利影响;但当降水量达到240mm以上时，果实则因水分太多导致采摘时果实裂口、黑果病等，丰产不丰收[14].因此，在选取气象评判指标时，只考虑温度和降水量这2个因素.
1.3.2生长指标
监测全生育期枸杞株高和叶面积生长指标，每个小区选取3棵样株，每隔10d测量1次.株高选用米尺测量从枸杞根部到顶部的距离，选用整个生育期的平均株高.叶面积测量选用便携式叶面积仪(UPA-210APB)进行测量，选用最大平均叶面积值.
1.3.3产量指标
产量指标包括枸杞的鲜果产量(简称：鲜产)、干果产量(简称：干产)、百粒重和粒度.
鲜产和干产：分茬采摘每个处理小区10株全部枸杞夏果，用电子秤(精度为0.01g)称重测定枸杞鲜果产量和烘干后的干果产量.
百粒重及粒度：每个处理测定干果产量后，随机选择100粒用电子秤(精度0.01g)称重，每处理重复3次取最大平均值;用电子秤(精度0.1g)称重50g，并清查枸杞干果的粒数，重复3次.
1.3.4品质指标
每个处理选取0.5kg枸杞样品送交专业检测机构，检测不同处理枸杞果中总糖、枸杞多糖、蛋白质和甜菜碱等品质指标.枸杞多糖和总糖含量根据标准《枸杞》(GB/T18672—2014)测定并评判，蛋白质含量根据标准《食品安全国家标准：食品中蛋白质的测定》(GB/T5009.5—2016)测定，甜菜碱含量根据标准《甜菜中甜菜碱的测定比色法》(NY/T1746—2009)测定并评判.
1.4基于AHP和熵权法的枸杞综合评判体系
1.4.1评判指标体系的构建
1)枸杞综合评判层次模型.以科学合理性和全面性相结合的原则，查阅相关文献和《枸杞》(GB/T18672—2014)选定评判指标.通过分析，建立了以气候、作物生长、产量和品质4个一级指标作为准则层、12个二级指标作为评判因素层、3个决策指标为备选方案层，共12项具有代表性的评判指标，构建枸杞综合评判层次模型.具体模型层次结构如图2所示.
2结果与分析
2.1基于AHP和熵权法的枸杞综合评判
2.1.1基于AHP和熵权法评判因素权重的确定
1)基于AHP的因素权重确定.通过邀请10位相关专家以调查问卷的方式，对不同层级中的多个因素两两比较，给出判断值.将所有专家的意见平均后，构建各个层次的判断矩阵.利用Yaahp10.2计算一级指标和二级指标的权重，获得AHP逐层指标权重，见表2，ek和ei分别为一级指标和二级指标的局部权重，所有指标的一致性比率CR的值均小于0.1则认为判断矩阵都通过了一致性检验，层级分析法赋权合理.
2)基于熵权法的因素权重确定.将2018—2019年枸杞大田试验的实测值标准化后再进行权重计算.熵权法权重的具体计算过程可按公式(3)—(5)进行，计算结果见表3.
3)基于AHP和熵权法的综合权重评判值的确定.在Matlab中利用公式(6)计算综合权重wi，结果见表3.由于AHP法与熵权法是按照不同的评判标准进行计算的，因此12个评判指标在不同的方法下，计算的权值对总目标的影响不尽相同.选取2种评判方法中得分最高的4个指标为例，分析其对枸杞综合评判结果的影响.
AHP法结果得出对评判结果影响最大的4个指标：按影响由大至小为粒度u8，蛋白质u11，总糖u9，干果产量u7，总计占总评判值的54.6%.由熵权法得出的2018年和2019年对总目标影响程度最高的4个指标：2018年按影响由大至小为干果产量u6，百粒重u7，粒度u8，总糖u9，总计占总评判值的66%;2019年按影响由大至小为干果产量u6，粒度u8，鲜果产量u5，总糖u9，总计占总评判值69.9%.综合法得出的2018年和2019年影响程度最高的4个指标：2018年按影响由大至小为粒度u8，百粒重u7，干果产量u6，总糖u9，总计占总评判值78%;2019年按影响由大至小为粒度u8，干果产量u6，总糖u9，百粒重u7，总计占总评判值73.3%.可以发现粒度、干果产量、百粒重是对枸杞评判影响最大的指标，所占比例高于73.3%.另一方面，从评判指标数值的差异也可以看出，综合法能更好地反映各指标对水肥配施的响应，更符合枸杞的生产实际.
2.1.2基于AHP和熵权法的枸杞水肥方案决策结果
综合评判的结果决定着枸杞种植相对较优的水肥配施方案.因此在进行决策时，采用综合评判的方式，通过自然间断法对枸杞综合评判指数进行分类，使各个决策类别差异明显.综合评判指数的间断分类见表4.
各水肥处理方案最终的综合评判指数按公式(7)进行计算，并根据表4自然间断法进行决策分类，基于AHP和熵权法的枸杞水肥配施最终综合评判指数Sj和决策结果见表5.总体上，综合评判的权值均随水肥用量增加呈先升高后降低的趋势.2018与2019年可以推广的水肥方案分别是处理W2F2和W2F3以及W2F3和W3F1，存在评判指数最优区间，表明此区间内的水肥方案最有利于枸杞生长.
2.2基于AHP和熵权法评判的枸杞最优水肥方案构建
2.2.1综合评判下枸杞水肥对评判指标的影响
以综合评判指数为依据，灌溉定额、总施肥量与综合评判指数关系如图3所示，图中I为灌溉定额，F为总施肥量.由图可见，2018年指标高的区域主要在中水中肥区域，灌溉定额和总施肥量分别为250～280mm和320～400kg/hm2;2019年指标高的区域则出现在中水高肥和高水低肥区域，灌水定额和总施肥量分别为240～280mm和420～450kg/hm2，300～313mm和270～290kg/hm2.
3讨论
以枸杞为研究对象，通过前人[2-3]对枸杞产量和品质的研究可知，气候、枸杞植株的生长、灌溉定额及肥料用量都是影响枸杞产量和品质的重要因素.目前，多数研究人员均利用一种方法对作物进行评判，鲜有利用多个方法从多个方面对作物水肥方案进行综合评判.
为了综合分析各类指标，通过主客观评判相结合的方法对枸杞4类一级指标12类二级指标进行建模运算.根据评判指数，利用综合评判法得出的枸杞水肥方案，能够准确地反应枸杞的需水需肥变化，更为符合生产实际.
由综合评判法得出，2018年的综合评判指数随肥料量增加而先增大后降低，这一规律与张智等[4]在樱桃番茄上的研究结果一致.这反映了不同水肥方案对综合评判指数的影响，同时也说明了在中水中肥区间内水肥对枸杞的生长具有协同效应，在高水低肥和低水高肥区间则存在拮抗效应.其中2018年产量最高的为处理W2F2，其增产率为16.41%;2019年产量最高的为处理W2F3，其增产率则高达34.38%.其增产率与李永梅等[2]在水肥一体化下研究枸杞水肥配施的增产率相近.但2019年的综合评判指数则随肥料用量增加而增大，这与前人[6，9]的研究结果不相一致，是由于2019年全年有效积温偏低，呼吸和蒸腾强度相对较小，植株自身养分的有效积累量小，而有效降雨量却较2018年偏高3.0mm，水肥的协同效应使肥料中更多的有机质以水为介质被植株吸收，导致了肥料用量增加，使得相同灌溉定额条件下综合评判指标的最大值均出现在高肥F3处理中.
前人[16-17]研究的结果显示，滴灌条件下枸杞适宜的灌溉定额和总施肥量分别为237～435mm和357～552kg/hm2时，枸杞的产量最高、生长和品质能够达到最优.
通过综合评判得出的灌溉定额和总施肥量区间分别为252～278mm和359～405kg/hm2，范围更精确.以此为依据制定枸杞的水肥方案也更为可靠，便于准确实施.
4结论
1)通过AHP和熵权法相结合的方法构建了枸杞综合评判模型.该方法从影响枸杞综合效益的气候、生长、产量和品质的14类12项指标，主客观并行对不同水肥方案进行评判、决策，最终得出的可推广水肥方案符合实际，且更为可靠.
2)在以枸杞产量和品质为主导的枸杞综合评判中，发现粒度、百粒重、干果产量和总糖这4个指标的评判指数对最终决策起着决定性的作用，影响占总评判指数高达73.3%以上.因此，当评判指标数据缺失或不全面时，推荐可使用这4个指标对枸杞进行综合评判，同样能得出较为客观的结论.
3)该模型通过综合评判法对各水肥处理进行评判决策，通过回归分析，得出枸杞的最优水肥方案，灌溉定额为252～278mm(不含春、冬灌)，配比施肥量(N-P2O5-K2O)为197-67-131kg/hm2～203-68-135kg/hm2.
4)本模型仅利用试验期的数据分析了枸杞水肥方案，且评判指标仍有不足.下一步将更加深入研究不同水文年型下的综合评判在枸杞水肥方案中的应用.——论文作者：邓箴1，2，尹娟1，3，4，5?，吴军斌1，张海军1
相关期刊推荐：《排灌机械工程学报》创刊于1982年，主要报道党和国家有关排灌行业的方针政策、最新标准;报道排灌机械行业的新产品、新工艺及工程设计的新论点、新发明及具有研究意义的前沿课题;报道排灌机械行业的新产品、新工艺，报道机电排灌站新建、改造、管理等的新思路、新办法;并为各类企事业单位承办单位形象及产品宣传。
