氮磷硒配施对大豆产量及籽粒硒含量的影响*

陈 鹏， 孙亚波， 张世文， 周 毅， 周 成,， 谢 越， 汪建飞,， 邢素芝

(1.安徽科技学院农业农村部生物有机肥创制重点实验室 安徽蚌埠 233400；2.安徽理工大学地球与环境学院 安徽淮南 232001)

大豆[Glycinemax(Linn.) Merr.]是我国乃至世界上重要的豆科经济作物，由于其营养价值颇高，被冠以“豆中之王”“田中之肉”“绿色的牛乳”等美誉。大豆不仅含有丰富的蛋白质、碳水化合物、矿物质以及维生素，而且有促进肠道蠕动、降血糖、降血脂等功效，因此成为人们餐桌上不可或缺的一种食物[1]。近年来，一些具有保健功能的大豆食品被开发与生产出来，满足了人们日益增长的需求。

硒是人体必需的微量元素，被称为元素中的“抗癌之王”。当人体缺硒时，可能会引发大骨节病[2]、克山病[3]等多种疾病。有研究表明，中国是世界上缺硒最严重的国家之一，全国有72%的土地面积处于缺硒地带，其中30%为严重缺硒地区[4]。研究发现，人体主要摄入植物中的硒，这是动物和人类摄入硒最为安全有效的方式[5]。

相对于禾谷类粮食作物，大豆富硒能力较强，有研究表明，富硒大豆中生物可给性高达90%，建议将其作为富硒农产品供居民日常补硒[6]。无论是土壤施用硒肥，还是叶面喷施亚硒酸钠，均可显著提高大豆籽粒中硒的含量。同时，硒还能增强植物受到胁迫后的抗性和抗氧化能力，促进植物的生长发育[7]。此前研究关注的重点是大豆对外源硒的吸收与富集情况，而对大豆生长的影响及籽粒产量相关性研究则较少[8-9]。因此，本文采用二次回归通用旋转组合设计，研究氮、磷、硒肥配施对大豆产量以及籽粒硒含量的影响，以获得最优施肥方案，为富硒大豆的实际生产提供理论参考。

1 材料与方法

1.1 试验材料

1.1.1 供试土壤

盆栽试验的土壤取自安徽科技学院种植园试验田表土，土壤为下蜀黄土发育而成的黄褐土。将土壤风干、粉碎、混匀，反复过5 mm筛以确保土壤肥力一致。该土壤基本理化性状：pH为6.9，有机质质量分数为19.6 g/kg，碱解氮质量分数为112.4 mg/kg，速效磷质量分数为25.7 mg/kg，速效钾质量分数为196.7 mg/kg。

1.1.2 供试材料

供试大豆：品种为“中黄13”。

供试用盆：试验选用普通塑料盆钵，规格为22 cm(上口直径)×14 cm(下口直径)×17 cm(高)，每盆装6 kg土壤。

试验用肥料：尿素，N质量分数为46%，安庆市氮肥厂；过磷酸钙，P2O5质量分数为12%，安徽六国化工股份有限公司；亚硒酸钠，Na2SeO3质量分数为45.6%，广东先导稀贵金属材料有限公司。

1.2 试验方法

试验以氮肥施用量(X1)、磷肥施用量(X2)、硒肥施用量(X3)为试验编码因素，采用三因素五水平通用旋转组合设计，共计20个处理，每个处理设2次重复，施肥水平见表1。

表1 试验因素水平编码及施肥水平

根据设计方案，将各处理的氮、磷肥称好备用，亚硒酸钠配成对应用量的100 mL溶液(即按照表1中各水平编码所对应的1 kg土壤中亚硒酸钠的用量乘以每盆土壤的质量，然后溶解并在100 mL容量瓶中定容至100 mL)。将预先处理的土壤按每份6 kg分装，共40份，并分为2组，每组20份。按照每个处理3种肥料用量依次加入，每加入1种肥料都充分混匀土壤，装盆。每盆喷淋500 mL蒸馏水，放置5 d。

2018年4月27日播种，每盆播种5粒，5月19日间苗，每盆保苗2株，8月26日收获。

1.3 数据分析

试验数据采用Excel 2016和DPS数据处理系统进行分析。

2 结果与分析

2.1 数学模型的建立

二次通用旋转组合试验方案及结果见表2。采用DPS数据处理系统对单株籽粒产量和籽粒硒含量进行多元回归分析，得到产量(Y1)、籽粒硒含量(Y2)与氮肥施用量(X1)、磷肥施用量(X2)、硒肥施用量(X3)三因素的三元二次回归数学模型见式(1)和式(2)：

表2 各处理试验方案及结果

(1)

(2)

对式(1)进行拟合度检验，计算出F1=1.605F0.05(9，10)=3.02，表明回归系数达显著水平。同样对式(2)进行拟合度检验和显著性检验，得出F1=4.360F0.05(9，10)=3.02，回归系数也达到显著水平。由此可以看出，式(1)和式(2)均失拟不显著，且回归系数都达到显著水平，所以2个二次回归模型具有代表性，说明试验所选的3个因素对大豆单株籽粒产量和籽粒硒含量均有显著影响，可以作为预测模型[10]。

在α=0.10显著水平剔除不显著项后，简化后的回归方程见式(3)和式(4)：

(3)

(4)

2.2 单因素效应分析

2.2.1 单株籽粒产量的单因素效应

将式(1)中的任意2个因素固定在零水平，就可以得到单因素一元偏回归子模型，见式(5)～(7)；然后根据方程绘制N、P、Se三因素对大豆产量的影响图，见图1。

图1 N、P、Se三因素对大豆产量的影响

(5)

(6)

(7)

由图1可以看出，大豆单株籽粒产量受氮肥施用量的影响较磷肥和硒肥施用量的大。在设计范围内，氮肥施用量在-1.682～1.682水平上，产量随氮肥施用量的增加而升高；磷肥施用量在-1.682～0水平上，产量随之升高；磷肥施用量在0～1.682水平上，产量下降；硒肥施用量则在-1.682～1.682水平上，产量随硒肥施用量的增加而缓慢下降。这3个因素对大豆单株籽粒产量的影响顺序为N>P>Se。

2.2.2 籽粒硒含量的单因素效应

将式(2)中的任意2个因素固定在零水平，得到单因素一元偏回归子模型，见式(8)～(10)；然后绘制N、P、Se三因素对大豆籽粒硒含量的影响图，见图2。

图2 N、P、Se三因素对大豆籽粒硒含量的影响

(8)

(9)

(10)

由图2可以看出，大豆籽粒硒含量主要受硒肥施用量影响。在设计范围内，氮肥施用量在-1.682～1水平上，籽粒硒含量先随之升高然后下降；磷肥施用量在-1.682～0水平上，籽粒硒含量随磷肥施用量的增加缓慢升高，然后磷肥施用量在0～1.682水平上籽粒硒含量缓慢下降；在-1.682～1.682水平上，籽粒硒含量则随着硒肥施用量的增加快速升高。这3个因素对大豆籽粒硒含量的影响顺序为Se>N>P。

2.3 两因素交互效应

经过对式(1)的交互项检验后发现，氮肥与硒肥(X1X3)、磷肥与硒肥(X2X3)的交互效应不显著，故分析氮肥与磷肥(X1X2)的交互作用。将式(1)中硒肥施用量固定在零水平，得到大豆单株籽粒产量与氮肥和磷肥的二因素子模型见式(11)：

(11)

通过式(11)可以计算出氮肥与磷肥的交互效应，见图3。

图3 氮肥与磷肥的交互效应

从图3可以看出：在低氮水平时，大豆产量随磷肥施用量的增加呈先升高后下降的趋势，说明氮肥施用量少时不宜多施磷肥；在中、高氮水平时，大豆产量则随着磷肥施用量的增加逐渐升高，表明氮肥施用量较多时，应适当增加磷肥的施用量；当X1=1.682、X2=1时，X1X2的交互效应达到最大值，即产量达到最大值。

对式(2)的交互项检验后发现，氮肥与磷肥(X1X2)、磷肥与硒肥(X2X3)的交互效应不显著，所以只分析氮肥与硒肥(X1X3)的交互作用。将式(2)中的磷肥施用量固定在零水平，得到大豆籽粒硒含量与氮肥和硒肥的二因素子模型见式(12)：

(12)

通过式(12)可以计算出氮肥与硒肥的交互效应，见图4。

图4 氮肥与硒肥的交互效应

由图4可知：硒肥施用量在-1.682～0水平时，籽粒硒含量随着氮肥施用量的增加呈先升高后下降的趋势，表明硒肥施用量少时不宜多施氮肥；硒肥施用量在0～1.682水平时，籽粒硒含量随着氮肥施用量的增加而逐渐升高，表明硒肥施用量较多时，应同时增加氮肥的投入量；当X1=1.682、X3=1.682时，X1X3的交互效应达到最大值，即籽粒硒含量达到最大值。

2.4 最优组合分析

2.4.1 产量最优组合

经过计算机模拟，在试验范围内，大豆的最高单株籽粒产量为12.86 g，其处理组合为X1=1.682、X2=1、X3=-1.682，即氮肥施用量为0.80 g/kg，磷肥施用量为7.97 g/kg，硒肥施用量为1.20 mg/kg。采用频率分析法对产量回归模型进行再解析，结果见表3。

表3 单株籽粒产量高于11.13 g的数值与频率

从表3可以看出，大豆单株籽粒产量大于11.13 g的组合方案有50个，X1、X2和X3变量的95%置信区间分别为0.902～1.244、0.177～0.696和-0.343～0.343，相当于氮肥施用量为0.616～0.699 g/kg，磷肥施用量为5.526～7.067 g/kg，硒肥施用量为2.282～2.838 mg/kg。

2.4.2 硒含量最优组合

经模拟计算，在试验范围内，大豆籽粒最高硒含量为9.84 μg/g，其处理组合为X1=1.682、X2=-1.682、X3=1.682，即氮肥施用量为0.80 g/kg，磷肥施用量为0，硒肥施用量为3.92 mg/kg。同样用频率分析法对硒含量回归模型进行解析，结果见表4。

表4 大豆籽粒硒含量高于5.09 μg/g的数值与频率

由表4可知，大豆籽粒硒含量大于5.09 μg/g的组合方案有70个，X1、X2和X3变量的95%置信区间分别为-0.159～0.400、-0.290～0.290和0.801～1.115，相当于氮肥施用量为0.362～0.496 g/kg，磷肥施用量为4.139～5.861 g/kg，硒肥施用量为3.209～3.463 mg/kg。

3 结语

近年来，二次通用旋转组合设计作为一种高效的数学模型被广泛应用于各种领域[11-13]，在计算黄秋葵[14]、玉米[15]、豇豆[16]等作物的最佳施肥量方面也很常见。有研究表明，硒不仅能提高大豆中蛋白质和氨基酸含量[17]，而且可以提高油脂的氧化稳定性和脂质含量[18]。由于硒元素在人体内不可合成，而大豆具有极其丰富的营养价值，将二者有机结合，既提高了大豆的营养价值，又保障了人体健康，应用前景非常广阔。

本文通过计算机模拟得到的数学模型反映了大豆生产中氮、磷、硒等3种肥料的配比与产量和硒含量的关系，为后续大豆富硒试验以及生产提供了一定的科学依据。但试验结果是在特定的试验环境下得到的，具有一定的局限性。另外，大豆籽粒中硒的形态及其与人体吸收之间的关系还不很明确，若应用于实际生产，有待进一步深入研究。