饱和D-最优设计在高蛋白大豆施肥优化中的应用

饱和D-最优设计在高蛋白大豆施肥优化中的应用
田艺心;高凤菊;曹鹏鹏
【摘要】[目的]高蛋白大豆生产是我国农业供给侧改革发展的一个重点产业,合理
施肥是实现高蛋白大豆生产的一个必要环节.利用饱和D-最优设计对黄淮海地区高蛋白大豆氮、磷、钾肥施用量进行研究,旨在为高蛋白大豆生产中的肥料配施提供
理论依据和参考.[方法]以高蛋白大豆冀豆21为试验材料,采用氮、磷、钾三因素二次饱和D-最优设计,每个因素分别设10个处理,每个处理3次重复,共计30个小区,随机区组排列.各处理肥料均在大豆开花期分批次施用,成熟后测定大豆产量及籽粒
蛋白质含量.[结果]建立了以氮、磷、钾施用量编码值为变量因子,高蛋白大豆产量
和蛋白质含量为目标函数的三元二次多项式数学模式.通过对模型解析寻优表明,氮、磷、钾肥对高蛋白大豆和蛋白质含量均有显著影响,且氮肥＞钾肥＞磷肥;当氮、磷、钾肥用量分别为95.46、183.8和128.7 kg/hm2时,边际产量效应值为0,当氮、磷、钾肥分别为120.8、178.4和141.3 kg/hm2时,边际蛋白质含量效应值为0.
氮磷、氮钾、磷钾对高蛋白大豆产量和蛋白质含量存在明显的交互作用.本试验条
件下,大豆产量超过3104 kg/hm2,籽粒蛋白质含量46.04％以上的施肥方案为施氮量76.13～101.1 kg/hm2,施磷量131.1～168.5 kg/hm2,施钾量104.9～134.8
kg/hm2.[结论]利用饱和D建立的肥料函数模型可以很好的说明氮、磷、钾施肥与大豆产量、蛋白质含量的连续变化关系和氮、磷、钾间的互作效应,适量的氮肥、
磷肥、钾肥可有效提高高蛋白大豆产量和蛋白质含量.供试条件下,最佳肥料配比是
N76.1～101.1 kg/hm2、P 131.1～168.5 kg/hm2、K 104.9～134.8 kg/hm2.【期刊名称】《植物营养与肥料学报》
【年(卷),期】2019(025)002
【总页数】8页(P343-350)
【关键词】高蛋白大豆;饱和D-最优设计;施肥;产量;蛋白质含量
【作者】田艺心;高凤菊;曹鹏鹏
【作者单位】德州市农业科学研究院,山东德州253015;德州市农业科学研究院,山东德州253015;德州市农业科学研究院,山东德州253015
【正文语种】中文
近年来，国内外市场对大豆蛋白需求猛增[1-2]，大豆产业尤其是高蛋白大豆发展迅速，如何提高高蛋白大豆产量并获得最佳的经济效益已成为大豆研究者和生产者广泛关注的热点。

黄淮海地区是我国大豆主要产区之一，已选育出大批优良高蛋白大豆品种，由于不同农田资源配置及生态条件差异，大豆难以有效发挥优势潜力，产量层次不齐。

施用氮、磷、钾肥料是大豆的主要栽培措施之一，对产量具有提高作用。

多数生产者为提高产量，往往施用大量肥料，导致肥料过量，肥料利用率降低，不仅造成环境恶化、成本增加、甚至造成大豆减产[3-5]。

针对大豆生产盲目施肥、过量施肥的情况，许多学者进行了相关研究，但多数研究集中在单一养分氮[6-9]或磷[10-11]或钾[12-14]，氮磷钾优化配施研究相对较少[15-16]，仅在施肥量、施肥时间等方面进行研究。

由于肥料之间存在交互效应，平衡施肥才能保证大豆正常生长发育和提高产量，优化氮、磷、钾平衡施肥，对提高肥料利用率、改善生态环境、提高大豆产量具有重要意义。

饱和D-最优设计是回归方程中参数数目与试验处理组合数目相等的设计，回归自由度与总自由度相等达到饱和，这种设计试验处理数少、误差小、精确度高、获得的信息量大，对选择最佳生产措施有较强的实用性，在农业试验中已多有应用。

如
应用二次饱和D-最优设计，刘伟明等建立了甘薯产量回归模型，对甘薯扦插密度、氮钾肥用量进行了优化[17]，董英等对青椒干燥因子进行了优化[18]，石云平等对龙鳞百合小鳞茎诱导培养基因子进行了优化[19]，腾云等对大豆水分和氮肥用量进行了优化[20]等。

肥料函数模型可以更好的说明氮、磷、钾施肥与大豆产量的连续变化关系，也可以研究氮、磷、钾间的互作效应。

因此，本研究利用饱和D-最优
设计对黄淮海地区高蛋白大豆氮、磷、钾肥施用量进行研究，以期获得最高大豆产量的适宜施肥量，为高蛋白大豆生产中的肥料配施提供理论依据和参考。

1 材料与方法
1.1 试验材料
供试大豆品种为‘冀豆21’，由ms1雄性不育轮回群体选择培育而成，蛋白质含量45.38% (超过国家高蛋白45%标准)，脂肪含量为17.87%，为高蛋白大豆品种，适宜黄淮海及华北地区种植，由河北省农林科学研究院粮油作物研究所提供。

1.2 试验设计
试验于2016—2017在山东省德州市农业科技园内进行。

供试大豆品种于2016
年6月11日播种育苗，行距为40 cm、株距为20 cm，小区面积15 m2。

供试
土壤为壤土，pH 7.53，有机质含量11.3 g/kg、全氮0.380 g/kg、碱解氮48.2 mg/kg、速效磷20.0 mg/kg、速效钾87.0 mg/kg。

前茬为冬小麦，麦收后秸秆还田。

试验采用氮、磷、钾三因素二次饱和D-最优设计，每个因素分别设10个
处理，每个处理重复3次，共计30个试验小区，随机区组排列。

各肥料处理均在大豆开花期分批次施用，除施肥量不同外，其他均按照当地生产管理水平进行田间管理，各小区田间操作保证一致。

成熟后对试验小区分别实收测产，脱粒后自然晒干，称量小区籽粒产量折合成公顷产量。

籽粒蛋白质含量采用近红外9100谷物分析仪测定。

具体试验因素水平编码及肥料用量如表1所示。

1.3 数据分析
利用SPSS15.0和Excel2007对所得数据进行统计分析，用MATLAB 7.0.1绘制三维图。

2 结果与分析
2.1 模型建立
为了揭示不同数量氮肥、磷肥、钾肥相互配合对大豆产量的影响，以表1中施氮量 (x1)、施磷量 (x2)、施钾量 (x3) 编码值为自变量，表2中籽粒产量 (y1)、籽粒蛋白质含量 (y2) 分别为目标函数，得到大豆产量、籽粒蛋白质含量与各因素之间的回归方程分别为：
表1 试验设计水平编码及肥料用量Table 1 Horizontal encoding of experimental design and fertilization amount处理Treatment肥料用量(kg/hm2)Fertilizer amount N P K N P K 1 -1 -1 -1 30 50 40 2 1 -1 -1 120 50 40 3 -1 1 -1 30 200 40 4 -1 -1 1 30 50 160 5 -1 0.1925 0.1925 30 139.44 111.55 6 0.1925-1 0.1925 83.66 50 111.55 7 0.1925 0.1925-1 83.66 139.44 40 8 -0.29121 1 61.90 200 160 9 1 -0.2912 1 120 103.16 160 10 1 1 -0.2912 120 200 82.53水平编码Level code
表2 各处理大豆的产量和蛋白质含量Table 2 Yield and protein content of soybean with different treatments处理Treatment蛋白质含量 (%)Protein content 1 2515 45.53 2 2870 45.86 3 2807 45.51 4 2833 45.82 5 3025 45.92 6 3180 46.05 7 3090 45.97 8 3250 46.21 9 3242 46.22 10 3260 46.32大豆产量 (kg/hm2)Yield
对两个方程分别进行F检验，F产量 = 2075＞ F0.01(9，18) = 3.60，F蛋白 = 2322＞ F0.01 (9，18)，说明模型的回归关系达到极显著水平。

进一步对各回归系
数进行T检验，产量方程回归系数除T (1，3) 差异达到5%显著，其余均达到1%极显著水平。

蛋白方程各回归系数T检验值则均达到1%极显著水平。

因此，产量和蛋白方程能够客观的反映施肥与大豆产量、施肥与籽粒蛋白质含量之间的关系，可以作为大豆产量和籽粒蛋白质含量预测的依据。

2.2 模型解析
2.2.1 主因素效应分析试验设计因素已经过量纲化处理，偏回归系数已经标准化，回归系数绝对值大小可直接反映出相应因素对因变量的影响程度。

综合考虑方程中各因素的偏回归系数与T检验结果，从方程 (1) 一次项的回归系数可以看出，本试验中三个因素氮、磷、钾肥偏回归系数绝对值分别为147.3、99.76、121.9，方程 (2) 一次项偏回归系数绝对值分别为0.21、0.11、0.16，说明氮、磷、钾肥对大豆产量和籽粒蛋白质含量作用大小的顺序为：氮肥主效应＞钾肥主效应＞磷肥主效应，高蛋白大豆花期施肥需氮肥、钾肥量要相对多于需磷量。

2.2.2 单因素效应分析将大豆产量和籽粒蛋白质含量回归模型中三个自变量中的任意两个固定在0码值，可得到剩余自变量与目标函数的关系，即氮、磷、钾肥对大豆产量 [方程(3)、(4)、(5)]和籽粒蛋白质含量 [方程(6)、(7)、(8)] 的单因素效应模型。

1) 大豆产量：
2) 大豆蛋白质含量：
将各单因素效应方程绘图，由图1可以看出，氮、磷、钾的产量单因子效应曲线均开口向下，尤其是氮肥、钾肥施用量较少的情况下，大豆产量均随着氮、磷、钾施用量的增加而增加，达到最高产量后，大豆产量增幅逐渐减少，特别是氮、钾编码值超过0.5以后，已经没有明显的增产效果。

籽粒蛋白质含量在较低施肥量条件
下，随氮、磷、钾施用量的增加而提高，但过量投入，尤其是钾编码值超过0以后，磷编码值超过0.5之后，籽粒蛋白质含量降低。

图1 各单因素对大豆产量和籽粒蛋白质含量的效应Fig. 1 Effects of single-factors on soybean yield and grain protein content
进一步对单因素效应方程求一阶偏导数进行产量和蛋白质含量边际效应分析。

边际效应分析的目的是对单位投入增量的收益进行分析。

由此可求得：
将各一阶偏导数方程绘图，由图2可以看出，氮、磷、钾三因素对大豆产量和蛋
白质含量的边际效益均呈递减趋势。

从大豆产量偏导数方程 (9)、(10)、(11) 斜率
绝对值可以看出，氮的边际效益递减率＞钾的边际效益递减率＞磷的边际效益
递减率；
从大豆蛋白质含量偏导数方程 (12)、(13)、(14) 斜率绝对值可以看出，钾的边际
效益递减率＞氮的边际效益递减率＞磷的边际效益递减率。

求解偏导数方程，
可得出最适氮肥投入量为0.4546 (编码值)，实际用量95.46 kg/hm2，此时大豆
产量为3327 kg/hm2；最适磷肥投入量为0.7838 (编码值)，实际用量183.8
kg/hm2，此时大豆产量为3334 kg/hm2；最适钾肥投入量为0.4782 (编码值)，实际用量128.7 kg/hm2，此时大豆产量为3323 kg/hm2。

同理，最适氮肥投入
量为1.0172 (编码值)，实际用量120.8 kg/hm2，此时籽粒蛋白质含量为46.29%；最适磷肥投入量为0.7120(编码值)，实际用量178.4 kg/hm2，此时籽粒蛋白质含量为46.22%，最适钾肥投入量为0.6887 (编码值)，实际用量141.3 kg/hm2，此时籽粒蛋白质含量为46.24%。

图2 各单因素对大豆产量和籽粒蛋白质含量的边际效应Fig. 2 Marginal effects
of single factors on soybean yield and grain protein content
2.2.3 因素耦合效应分析大豆产量和籽粒蛋白质含量受多因素的影响，由方程(1)、
(2) 可知，氮磷、氮钾、磷钾交互项其偏回归系数均达到显著水平，说明氮、磷、钾两两之间的交互效应对大豆产量和籽粒蛋白质含量产生显著影响。

采用降维法将回归模型 (1)、(2) 中的氮、磷、钾素固定为0码值，研究其余两个因素的交互作用，并且作曲面图。

从图3可以看出，在编码值范围内，氮、磷、钾肥对大豆的产量和蛋白质含量效应都呈抛物线性。

氮磷肥互作时，施磷量一定时，大豆产量随氮肥量的增加而增加；施氮量一定时，随磷肥施用量的增加，大豆产量在磷肥编码-1到0.5区间增加，在0.5～1区间则降低。

氮钾肥互作时，施钾量一定时，随着氮肥用量的增加，大豆产量在氮肥编码-1到0.5区间增加，在0.5～1区间降低；反之钾肥亦然。

磷钾互作时，施钾量一定时，大豆产量在磷肥编码-1到0区间随磷肥增加而增加，在0～1区间则降低；施磷量一定时，大豆产量随钾肥量增加而增加。

氮磷、氮钾、磷钾交互项对蛋白质含量的影响相似，两两之间均具有明显的交互促进作用。

由图3分析得出，大豆产量≥ 3104 kg/hm2，氮磷理想互作编码区间为氮肥取0.027～0.551，磷肥取-0.277～0.382，即施氮量76.24～99.82 kg/hm2，施磷量104.2～153.7 kg/hm2；氮钾理想互作编码区间为氮肥和钾肥均取0.002～0.598，即施氮量75.09～101.9 kg/hm2，施钾量100.1～135.9 kg/hm2；磷钾理想互作区间为磷肥取-0.227～0.427，钾肥取-0.056～0.506，即施磷量108.0～157.0 kg/hm2，施钾量96.66～130.3 kg/hm2。

同理，由图3分析得出，籽粒蛋白质含量≥46.04%，氮磷理想互作编码区间为氮肥取0.079～0.627，磷肥取-0.190～0.484，即施氮量78.54～103.2 kg/hm2，施磷量110.7～161.3 kg/hm2；氮钾理想互作编码区间为氮肥取0.052～0.636，钾肥取-0.105～0.542，即施氮量77.33～103.6 kg/hm2，施钾量93.68～132.6 kg/hm2；磷钾理想互作区间为磷肥和钾肥均取0.002～0.598，即施磷量125.2～169.9 kg/hm2，施钾量 100.1～135.9 kg/hm2。

2.3 模型寻优
对具有同等重要性的试验指标，通过频率出现的高低可优先排序，频率高的因素水平可被优先选择，从而获得最佳置信区间和配比方案。

采用频率分析法对模型进行寻优，由表3可知，高蛋白大豆冀豆21要获得超过3104 kg/hm2的产量时，施用氮、磷、钾肥最优组合取值范围为：氮肥施用量75.09～101.9 kg/hm2，磷肥施用量125.2～169.9 kg/hm2，钾肥施用量100.1～135.9 kg/hm2。

籽粒蛋白质含量超过46.04%时，施用氮、磷、钾肥最优组合取值范围为：氮肥施用量
75.16～103.8 kg/hm2，磷肥施用量129.6～173.9 kg/hm2，钾肥施用量103.7～139.2 kg/hm2。

为了提高大豆产量的同时保证籽粒蛋白质含量，采用交集法可以得到大豆产量超过3104 kg/hm2，籽粒蛋白质含量在46.04%以上的施肥量为施氮量75.16～101.9 kg/hm2，施磷量129.6～169.9 kg/hm2，施钾量103.7～135.9 kg/hm2。

图3 氮、磷、钾与大豆产量和蛋白质含量的两因素互作效应Fig. 3 Two-factor interactive effect between N, P, K with yield and protein content of soybean
表3 大豆产量大于3104 kg/hm2的因素取值频率分布Table 3 Frequency distribution of the factors for soybean yield over 3104 kg/hm2水平编码Horizontal encoding X1 X2 X3次数Count频数Frequency 次数Count频数Frequency 次数Count频数Frequency-1.0 0 0 0 0 0 0-0.5 3 20 3 20 3 20 0 4 26.7 4 26.7 4 26.7 0.5 4 26.7 4 26.7 4 26.7 1.0 4 26.7 4 26.7 4 26.7加权均数Weighted average 0.3000 0.3000 0.3000 SD 0.1398 0.1398 0.1398 95% 置信区间 95% confidence interval 0.0020～0.5980 0.0020～0.5980 0.0020～
0.5980最优配比方案 Optimal matching scheme 75.09～101.91 125.15～169.85 100.12～135.88
3 讨论与结论
大豆是需肥量较大的作物，氮、磷、钾对大豆产量及品质形成具有重要影响。

大豆高产优质，主要取决于氮磷钾肥的平衡吸收，优化和平衡大豆施肥可促进大豆正常生长发育、干物质积累和共生同氮作用，进而提高大豆产量和品质。

大豆开花至鼓粒期处于养分吸收和物质积累的高峰期，营养生长与生殖生长旺盛，60%～70%的氮磷钾养分在这个时期被吸收，故应特别注重花荚期养分的充分供给。

本试验对主因子效应研究表明，氮、磷、钾肥对大豆产量和籽粒蛋白质含量均有显著影响，均表现为氮肥最大，钾肥次之，磷肥最小。

氮是大豆细胞质、细胞核和各种酶的组成成分，与籽粒蛋白质、脂肪和碳水化合物形成有密切关系，花期施用氮肥有利于减少花荚脱落，构建大豆良好个体及群体形态，增加光能利用率，促进有机物质积累和养分供给[21]，因而施氮有利于大豆产量和籽粒蛋白质含量的提高。

钾能促进光合作用，并促进碳水化合物向籽粒的运输，最终促进大豆产量和籽粒蛋白质的形成。

磷可以激活并释放以ATP形式贮存的能量，还可促进光合产物的合成转化和运输，使有机物更多地向根系分配[22]，对提高大豆产量和品质也具有重大作用。

但也有学者研究表明氮、磷、钾肥对大豆产量的效应排序为氮肥＞磷肥＞钾肥[23-24]，对蛋白质含量的效应排序为磷肥＞氮肥＞钾肥[25-26]，本试验中磷肥对大豆产量的影响相对较小，原因可能在于土壤本身含有的磷肥一定程度上掩盖了施用磷肥的效果，可能与不同地区土壤肥力、气候环境不同有极大关系，另一方面可能与大豆品种本身营养特性、大豆不同生育期有关。

本试验中，对大豆产量和籽粒蛋白质含量的单因子效应分析结果表明，大豆产量和籽粒蛋白质含量均先随着氮、磷、钾施用量的增加而增加，达到最高产量后增幅逐渐减少，这符合肥料的报酬递减定律，过低或过量施肥均不利于产量和品质提高。

另外，三种肥料的边际产量效应和边际蛋白质效应规律类似，均随着肥料投入量的
增加而不断减小，当氮、磷、钾肥用量分别达95.46、183.8和128.7 kg/hm2时，边际产量效应值降至0；当氮、磷、钾肥用量分别达120.8、178.4和141.3
kg/hm2时，边际蛋白质含量值降至0。

氮、磷、钾对植物的影响并不是单独效应的线性累加，存在着交互作用。

因此，在施肥过程中，需要综合考虑交互效应，寻求施肥的平衡点，才能得到合理的施肥量。

本试验中，氮磷、氮钾、磷钾互作对大豆产量和籽粒蛋白质含量均有显著影响，且两两存在一定范围的相互促进作用。

籽粒蛋白质含量是高蛋白大豆品种中最关键的有效成分，结合产量及蛋白质含量进行肥料优化，本试验得出，大豆产量超过3104 kg/hm2，籽粒蛋白质含量在46.04%以上的施肥量为施氮量76.13～101.1 kg/hm2，施磷量 131.1～168.5 kg/hm2，施钾量 104.9～134.8 kg/hm2。

用码值表示的饱和D-最优设计等现代回归设计方案作为一种数学工具，其信息矩
阵行列式最大，在设置重复的条件下试验效率最高，可普遍应用于包括农业在内的各个领域[27-28]。

利用其所产生的肥料效应方程，能够客观地反映施肥量与产量
及品质的关系，从而得到适宜施肥量。

但由于实际农业中，供试土壤肥力、栽种品种、植株受到气候、水分、栽培措施等多种因素差异很大，因此在应用中还需根据实际情况具体分析。

参考文献：
【相关文献】
[1]殷瑞锋, 徐雪高, 冯学静. 2016年以来大豆市场分析与后市展望[J]. 市场预测, 2017, 9(6): 11-16.Yin R F, Xu X G, Feng X J. Soybean market since 2016 and its future prospect[J]. Market Forecasting, 2017, 9(6): 11-16.
[2]Ainsworth E A, Yendrek C R, Skoneczka J A, et al. Accelerating yield potential in soybean: potential targets for biotechnological improvement[J]. Plant Cell&Environment,
2012, 35(1): 38-52.
[3]白由路. 植物营养与肥料研究的回顾与展望[J]. 中国农业科学,2015, 48(17): 3477-3492.Bai Y L. Review on research in plant nutrition and fertilizers[J].Scientia Acricultura Sinica, 2015,
48(17): 3477-3492.
[4]Sebilo M, Mayer B, Nicolardot B, et al. Long-term fate of nitrate fertilizer in agricultural soils[J]. Proceedings of the National Academy of the Sciences of the United States of America, 2013,110(45): 18185-18189.
[5]王火焰, 周健民. 肥料养分真实利用率计算与施肥策略[J]. 土壤学报, 2014, 51(2): 216-
225.Wang H Y, Zhou J M. Calculation of real fertilizer use efficiency and discussion on fertilization strategies[J]. Acta Pedologica Sinica, 2014,51(2): 216-225.
[6]谌俊旭, 黄山, 范元芳, 等. 单作套作大豆叶片氮素积累与光谱特征[J]. 作物学报, 2017, 43(12): 1835-1844.Shen J X, Huang S, Fan Y F, et al. Remote detection of canopy leaf nitrogen status in soybean by hyper spectral data under monoculture and intercropping systems[J]. Acta Agronomica Sinica, 2017, 43(12):1835-1844.
[7]宋柏权. 不同基因型大豆产质量形成及其氮素调控[D]. 哈尔滨: 东北农业大学博士学位论文, 2009.Song B Q. Differences between yield and quality of different genotypes soybean and its adjustment of nitrogen[D]. Harbin: PhD Dissertation of Northeast Agricultural University, 2009.
[8]夏玄, 龚振平. 氮素与豆科作物固氮关系研究进展[J]. 东北农业大学学报, 2017, 48(1): 79-88.Xia X, Gong Z P. Research advance on the relationship between nitrogen and leguminous nitrogen fixation[J]. Journal of Northeast Agricultural University, 2017, 48(1): 79-88.
[9]姜妍, 王浩, 王绍东, 等. 不同氮素水平对7S亚基缺失大豆品质及相关性状的影响[J]. 核农学报, 2017, 31(3): 574-579.Jiang Y, Wang H, Wang S D, et al. Effect of different nitrogen level to quality characteristics and other related traits in the 7S subunit deficient soybean[J]. Journal of Nuclear Agricultural Sciences, 2017,31(3): 574-579.
[10]尹元萍, 董文汉, 王明君, 等. 2种磷水平下大豆种质资源品质性状的遗传差异分析[J]. 西南农业学报, 2017, 30(10): 2185-2190.Yin Y P, Dong W H, Wang M J, et al. Analysis of genetic variation of soybean germplasm resources under two phosphorus levels[J].Southwest China Journal of Agricultural Sciences, 2017, 30(10):2185-2190.
[11]钟鹏, 刘杰, 王建丽. 不同磷效率基因型大豆产量的构成因素及籽粒产量通径[J]. 农业工程, 2016, 6(2): 110-115.Zhong P, Liu J, Wang J L. Yield component factors and grain yield path of genotype soybean with different phosphorus efficiency[J].Agricultural Engineering, 2016, 6(2): 110-115.
[12]朱晓玲, 陈海峰, 王程, 等. 大豆钾转运体基因GmKT12的克隆和信息学分析[J]. 作物学报, 2013, 39(9): 1701-1709.Zhu X L, Chen H F, Wang C, et al. Molecular cloning and bioinformatics analysis of K+ transporter gene(GmKT12) from soybean (Glycine max L. Merri)[J]. Acta Agronomica Sinica, 2013,39(9): 1701-1709.
[13]Parvej M R, Williams A S, Holshouser D L, et al. Double-crop soybean response to potassium on mid-atlantic coastal plain and piedmont soils[J]. Agronomy Journal, 2018, 110(1): 399-410.
[14]Carmona M A, Sautua F J, Grijalba P E, et al. Effect of potassium and manganese phosphates in the control of Pythium damping-off in soybean: a feasible alternative to fungicide seed treatments[J]. Pest Management Science, 2018, 74(2): 366-374.
[15]王海艳. 优化施肥对大豆氮磷钾吸收与分配的影响[D]. 哈尔滨：东北农业大学硕士学位论文, 2011.Wang H Y. Effects of optimized fertilization on NPK uptake and distribution of soybean[D]. Harbin: MS Thesis of Northeast Agricultural University, 2011.
[16]Wei Y, Yuting G, Xiaochun W, et al. Temporal variations of soil microbial community under compost addition in black soil of Northeast China[J]. Applied Soil Ecology, 2017, 121(12): 214-222.
[17]刘伟明. 二次饱和D-最优设计在甘薯栽培试验中的应用[J]. 安徽农业科学, 2011, 39(35): 21625-21626, 21633.Liu W M. Application of the secondary saturation and D-optimal design in the experiment in sweet potato cultivation[J]. Journal of Anhui Agricultural Science, 2011, 39(35): 21625-21626, 21633.
[18]董英, 陈立, 付西光. 饱和D最优设计在青椒干燥试验中的应用[J].农业机械学报, 2002, 33(3): 70-73.Dong Y, Chen L, Fu X G. Study on D-saturating optimum design for sweet pepper drying[J]. Transactions of the Chinese Society for Agricultural Machinery, 2002, 33(3): 70-73.
[19]石云平, 杨美纯, 李锋. 饱和D-最优设计在龙牙百合小鳞茎诱导中的应用[J]. 广西农业科学, 2009, 40(12): 1519-1523.Shi Y P, Yang M C, Li F. Application of saturation D-optimum design for induction of Lilium brownii var viridulum bulblet derived from squamule of bulb[J]. Journal of Guangxi Agricultural, 2009,40(12): 1519-1523.
[20]滕云, 郭亚芬, 张忠学, 等. 东北半干旱区大豆水肥耦合模式试验研究[J]. 东北农业大学学报, 2005, 36(4): 639-644.Teng Y, Guo Y F, Zhang Z X, et al. Study on couple effect of water and fertilizer of soybean in northeast semiarid area[J]. Journal of Northeast Agricultural University, 2005, 36(4): 639-644.
[21]Zhu X M, Han B J. Root water uptake and root nitrogen mass distributions of soybean under different nitrogen and water supply[J].Interaction Journal of Agriculture and Biology, 2018, 20(1): 52-56.
[22]Barbosa N C, Pereira H S, Arruda E M, et al. Spatial distribution of phosphorus in the soil and soybean yield as function of fertilization methods[J]. Bioscience Journal, 2018, 34(1): 88-94.
[23]孙振宁, 李晶, 段兴武, 等. 氮磷钾施肥水平对大豆产量及性状的影响[J]. 作物杂志, 2012, 5(35): 135-143.Sun Z N, Li J, Duan X W, et al. Effects of combined fertilization of N, P and K on soybean yield and yield traits[J]. Crops, 2012, 5(35):135-143.
[24]吴书平, 马永学, 纪永民. 氮磷钾优化施肥对夏大豆产量的影响[J].大豆科技, 2008, 2(11): 36-38.Wu S P, Ma Y X, Ji Y M. Effect of optimal fertilization of N, P and K on summer soybean yield[J]. Soybean Science & Technology,2008, 2(11): 36-38.
[25]纪永民, 张存岭, 曹杰. 氮磷钾追施对夏大豆产量的影响[J]. 磷肥与复肥, 2012, 27(2): 74-77.Ji Y M, Zhang C L, Cao J. The effect of dressing NPK on the yield of summer soybean[J]. Phosphate & Compound Fertilizer, 2012, 27(2):74-77.
[26]陈峰, 张开文, 刘波, 等. 氮磷钾配施对大豆产量的影响[J]. 大豆科技, 2011, 3(11): 18-22.Chen F, Zhang K W, Liu B, et al. Effect of different proportion of N,P, K on soybean yield[J]. Soybean Science & Technology, 2011,3(11): 18-22.
[27]Xiao Z, Jiang W, Yu H, et al. Substrate water content and nitrogen interactions in growing media: yield, fruit quality, water consumption and water use efficiency on tomato[J]. International Symposium on Soilless Culture and Hydroponics, 2009, 843(6): 57-64.
[28]任露泉. 试验优化设计与分析[M]. 北京: 高等教育出版社, 2002.344-347.Ren L Q. Experimental optimization design and analysis[M]. Beijing:High Education Press, 2002. 344-347.。