基于农学效应的福建省水稻氮磷钾施肥限量标准

李娟; 章明清; 章赞德; 张立成; 姚宝全; 沈金泉

doi:10.19303/j.issn.1008-0384.2019.02.002

基于农学效应的福建省水稻氮磷钾施肥限量标准

1.
福建省农业科学院土壤肥料研究所, 福建福州 350013
2.
大田县农田建设与土肥技术推广站, 福建大田 366100
3.
福建省农田建设与土壤肥料技术总站, 福建福州 350003

基金项目:

国家自然科学基金项目 31572203

福建省科技计划项目——省属公益类科研院所基本科研专项 2018R1022-3

福建省农业科学院科技创新团队建设项目 STIT2017-1-9

福建省测土配方施肥项目 2010-2015

详细信息

作者简介:
李娟(1977-), 女, 硕士, 副研究员, 研究方向:作物营养与施肥(E-mail:li-95@163.com)

通讯作者:
章明清(1963-), 男, 博士, 研究员, 研究方向:作物施肥原理与技术(E-mail:zhangmq@163.com)

中图分类号: S147.2
计量
- 文章访问数: 1897
- HTML全文浏览量: 351
- PDF下载量: 55
出版历程
- 收稿日期: 2018-12-12
- 修回日期: 2019-01-28
- 刊出日期: 2019-02-27

NPK Fertilizations Recommended Based on Agronomic Benefits on Rice in Fujian

1.
Soil and Fertilizer Institute, Fujian Academy of Agricultural Science, Fuzhou, Fujian 350013, China
2.
Datian Cropland Construction and Soil and Fertilizer Station, Datian, Fujian 366100, China
3.
Fujian Cropland Construction and Soil and Fertilizer Station, Fuzhou, Fujian 350003, China

摘要

摘要:
目的为提高水稻施肥效益，实现减量、增效的目标，制定福建早稻、晚稻、单季稻、菜后稻和烟后稻的氮磷钾施肥限量标准。
方法根据近年来在不同稻作上完成的832个氮磷钾田间肥效试验结果，将福建早稻、晚稻和单季稻分为山区和沿海两个区域，菜后稻不分区域，烟后稻则归入山区区域，应用系统聚类分析方法及其类别间差异显著性检验，确定各稻作的最佳施肥类别数，然后构建各施肥类别的氮磷钾三元非结构肥效模型及其推荐施肥量，并进行田间试验验证。
结果山区区域、沿海区域和全省的早稻最佳施肥类别数分别为7、6、7个，晚稻则分别为6、4、7个；单季稻对应区域的最佳施肥类别数分别为6、4、6个；菜后稻和烟后稻的最佳施肥类别数分别为3个和4个。在此基础上，针对各稻作不同施肥类别构建的60个氮磷钾三元非结构肥效模型均属于典型式，进而建立了基于农学效应的水稻氮磷钾施肥限量标准。234个大田对比试验表明，施肥限量标准优化了氮磷钾施用量和比例，不同稻作间平均增产稻谷4.0%~12.5%，净增收达到875~2 616元·hm^-2；推荐施肥的氮磷减肥幅度平均分别为19.4%和23.9%，但增加K₂O用量27.3%。
结论氮磷钾施肥限量标准为当前水稻减肥增效技术的推广应用提供了最佳养分投入的计量依据。
- 水稻 /
- 氮磷钾 /
- 施肥类别 /
- 限量标准 /
- 增产 /
- 效益
Abstract:
Objective To improve efficiency and reduce waste on rice fertilization in the fields while protecting the environment, the upper limits on N, P and K applications for the farming of early rice, late rice, single cropping rice, rice after vegetable, and rice after tobacco in Fujian were studied.
Method A total of 832 field experiments were conducted on (1) the early rice, late rice and single cropping rice in mountainous and coastal areas, (2) the after-tobacco cropping in mountain areas only, and (3) the after-vegetable cropping in either mountainous or coastal areas. System clustering analysis and significance test of difference were employed to classify the data generated. A ternary non-structural NPK fertilizer response model was then constructed for each category prior to field testing.
Result Seven optimum fertilization categories for the early rice in the mountainous areas, 6 in the coastal areas, and 7 in the entire province were established. For the late rice, 6 in the mountainous areas, 4 in the coastal areas, and 7 in the entire province were obtained. For the single cropping rice, 6 in the mountainous areas, 4 in the coastal areas, and 6 in the entire province were resulted. And, for the after-tobacco cropping, 4 in the mountainous areas; while, for the after-vegetable cropping, 3 categories became available for analysis. Based upon these 60 models on the NPK fertilization limits set according to the agronomic effects on the rice plants, 234 field tests were carried out for final selection. It was verified that the formulations produced 4.0% to 12.5% increases on yield with an averaged net income increase from 875-2 616 Yuan·hm^-2. In comparison to the common practice by the farmers, the recommended fertilizations reduced 19.4% on N usage and 23.9% on P but increased 27.3% on K₂O.
Conclusion The fertilization recommendations would optimize the fertilizer application with significantly improved efficiency for all types of rice cultivation currently practiced in the province.
- paddy rice /
- NPK /
- fertilization category /
- application limits /
- increased yield /
- benefits

HTML全文

0. 引言

【研究意义】水稻是福建省的主要粮食作物，研究和推广应用科学施肥技术，对提高施肥效益和水稻产量、实现减肥增效的目标具有重要意义。【前人研究进展】研究和分析农业生产对养分的需求规律历来受到人们的重视^[1]。在水稻方面，宋晓梅等^[2]、戢林等^[3]应用三元二次多项式、一元二次多项式和线性加平台等模型, 建立重庆市合川区、川中丘陵区不同稻田土壤肥力等级的水稻氮磷钾最佳施肥量。毛伟等^[4]则以地力差减法预测氮肥用量、养分丰缺指标法预测磷钾肥用量，构建了江苏省江都市的县域水稻施肥指标体系等。本课题组根据十一·五期间福建省不同稻作的试验结果，建立了以肥料效应函数-土壤养分丰缺指标法为技术路线的氮磷钾施肥指标^[5-7]。【本研究切入点】上述研究在确定土壤肥力等级即施肥类别时，大都根据无肥区产量水平及其变化幅度，简单地划分为3~6个等级^[8-9]，忽略了不同土壤肥力等级间在施肥效应和产量水平上是否存在显著差异的问题；在确定各肥力等级的推荐施肥量时，采用三元二次多项式肥效模型^[10]，忽略了该模型设定偏误与多重共线性和异方差^[11-12]对推荐施肥的影响。本课题组近年来研发的非结构肥效模型克服了二次多项式肥效模型存在的上述问题，具有更高的拟合精度和更宽的适用范围^[13]。因此，有必要对福建相关研究结果进一步更新和完善，为制定适合当前福建水稻生产的氮磷钾施肥限量标准。【拟解决的关键问题】根据近8年来完成的832个水稻氮磷钾田间肥效试验结果，在施肥类别间差异显著性检验^[14]基础上，探讨水稻氮磷钾施肥类别的最佳分类数，然后应用三元非结构肥效模型建立基于农学效应的氮磷钾施肥限量标准，旨在为稻田减肥增效提供最佳养分投入量计算依据。

1. 材料与方法

1.1 水稻氮磷钾田间肥效试验

截至2017年底，笔者完成和整理了来自水稻主产区氮磷钾“3414”设计的田间肥效试验资料832个，试验设计方法同文献[5]。在此基础上，为检验水稻氮磷钾施肥限量标准在生产上的应用效果，在福建水稻产区的代表性地块进行了234个田间大区对比试验。每个对比试验设3个处理，即：(1)空白；(2)习惯施肥；(3)推荐施肥(氮磷钾施肥限量标准)。其中，农户习惯施肥是根据试验农户对应稻作的常规施肥量及其施肥方法，推荐施肥是根据课题组研发的智慧农业推荐施肥软件(网址：http://120.27.225.89:9999/index.html)的施肥模式识别结果，或者根据施肥目标产量对应的施肥限量标准，确定施肥量及其施肥方法。小区面积60 m²，不设重复，田间管理措施参考文献[5]。

上述两种试验类型的供试水稻包括早稻、晚稻、单季稻、菜后稻和烟后稻等，试验点在全省各地的设置情况见表 1。供试土壤包括灰泥田、黄泥田、灰沙田等福建稻田主要土壤类型(土属)，用常规方法^[15]测定“3414”试验点的供试土壤主要理化性状和处理(6)的施肥量见表 2，大区对比试验的基础土样主要理化性状和各处理施肥量见表 3。

表 1 福建水稻氮磷钾田间肥效试验数量及其试验点设置情况

Table 1. Number and distribution of sites for field-testing rice growth in response to varied NPK fertilizations in Fujian

稻作 Rice	试验类型 Treatments type	试验点数 Number of trials	试验点分布情况Experiment distribution/个
稻作 Rice	试验类型 Treatments type	试验点数 Number of trials	福州市 Fuzhou	莆田市 Putian	泉州市 Quanzhou	漳州市 Zhangzhou	龙岩市 Longyan	三明市 Sanming	南平市 Nanping	宁德市 Ningde
早稻Early rice	3414试验	171	27	36	0	17	36	28	17	10
早稻Early rice	大区对比试验PCE	52	0	14	0	0	17	10	2	8
晚稻Late rice	3414试验	173	18	17	0	10	39	36	23	30
晚稻Late rice	大区对比试验PCE	58	0	3	2	2	24	9	14	3
单季稻Single cropping rice	3414试验	324	25	25	0	0	20	77	90	87
单季稻Single cropping rice	大区对比试验PCE	111	16	15	12	0	8	27	18	15
菜后稻Rice after vegetable	3414试验	63	2	0	0	11	29	9	7	5
菜后稻Rice after vegetable	大区对比试验PCE	4	0	0	2	2	0	0	0	0
烟后稻Rice after tobacco	3414试验	101	0	0	0	0	52	49	0	0
烟后稻Rice after tobacco	大区对比试验PCE	9	0	0	0	0	5	1	3	0
注：PCE表示大区对比试验。 Note: PCE stands for plot check experiment.

下载: 导出CSV

| 显示表格

表 2 氮磷钾“3414”肥效试验供试土壤主要理化性状及其处理(6)施肥量

Table 2. Physiochemical properties of soil at test sites with "3414" design and fertilizer applications of treatment 6

稻作 Rice	试验点数 Number of trials	土壤主要理化性状 Soil physical and chemical properties					处理(6)施肥量 Application rate of Treatment 6/(kg·hm^-2)
稻作 Rice	试验点数 Number of trials	pH	有机质 Organic matter/(g·kg^-1)	碱解氮 Alkali-hydr N/(mg·kg^-1)	有效磷 Olsen-P/(mg·kg^-1)	速效钾 Avail K/(mg·kg^-1)	氮肥 N	磷肥 P₂O₅	钾肥 K₂O
早稻ER	171	4.7±0.8	28.4±6.2	158.8±48.7	33.4±23.4	86.3±54.4	158.0±10.3	61.2±8.8	109.2±14.2
晚稻LR	173	4.6±0.8	31.5±7.2	149.1±56.6	32.4±24.7	88.9±55.3	157.0±10.4	59.0±9.3	112.5±40.0
单季稻SR	324	4.6±0.6	28.1±13.4	147.6±60.8	28.0±21.0	87.9±48.9	165.9±16.5	67.4±12.0	114.6±12.9
菜后稻RAV	63	5.4±0.4	32.4±7.7	162.7±41.5	42.8±18.9	97.8±46.7	158.3±18.9	62.0±16.2	107.0±19.8
烟后稻RAT	101	4.9±0.2	28.4±17.2	159.7±43.2	36.1±19.2	124.0±68.1	160.1±18.6	57.9±17.5	105.0±16.5
注：表中土壤理化性状测定值和处理(6)数据以平均值±标准差表示，ER表示早稻, LR表示晚稻，SR表示单季稻, RAV表示菜后稻, RAT表示烟后稻。表 3~5、9~10同。 Note:Soil physiochemical properties and treatment 6 are expressed as mean±standard deviation, ER stands for early rice, LR for late rice, SR for single-cropping rice, RAV for after-vegetable cropping, and RAT for after-tobacco cropping.The same as table 3-5, 9-10.

下载: 导出CSV

| 显示表格

表 3 大区对比试验的土壤主要理化性状和试验处理施肥量

Table 3. Physiochemical properties of soil at comparison test sites and fertilizer applications of treatments

稻作 Rice	试验点数 Number of trials	土壤主要理化性状 Soil physical and chemical properties					推荐施肥 Recommend fertilization/(kg·hm^-2)			习惯施肥 Farmer practice/(kg·hm^-2)
稻作 Rice	试验点数 Number of trials	pH	有机质 Organic matter/(g·kg^-1)	碱解氮 Alkali-lhydr. N/(mg·kg^-1)	有效磷 Olsen-P/(mg·kg^-1)	速效钾 Avail K/(mg·kg^-1)	氮肥 N	磷肥 P₂O₅	钾肥 K₂O	氮肥 N	磷肥 P₂O₅	钾肥 K₂O
早稻ER	52	5.3±0.6	27.2±11.5	149.1±45.7	30.7±22.9	67.5±36.0	169±20	59±24	119±31	184±24	72±18	71±36
晚稻LR	58	5.2±0.5	30.9±8.9	157.9±39.6	29.5±17.2	92.6±75.2	154±25	50±15	97±30	195±21	67±24	54±30
单季稻SR	111	5.2±0.4	29.2±10.4	144.1±48.7	34.2±25.4	81.3±55.5	165±17	63±16	110±28	174±37	71±21	61±39
菜后稻RAV	4	4.9±0.5	16.7±7.1	122.2±16.8	100.3±37.8	85.6±44.8	75±86	23±26	53±60	131±103	56±43	19±38
烟后稻RAT	9	5.4±0.6	38.1±9.3	175.1±56.1	29.3±10.3	104.2±53.3	133±29	43±11	51±23	122±58	39±36	48±49

下载: 导出CSV

| 显示表格

表 4 不同稻作施用氮磷钾化肥对稻谷产量的影响

Table 4. Effect of NPK fertilizations on yield of rice under varied cultivation practices

稻作 Rice	试验点数 Number of trials	经典五处理及其平均产量 Classical five treatments and the average yield/(kg·hm^-2)					平均贡献率 Average contribution rates/%
稻作 Rice	试验点数 Number of trials	对照CK	N₀P₂K₂	N₂P₀K₂	N₂P₂K₀	N₂P₂K₂	土壤Soil	氮N	磷P	钾K
早稻ER	171	5080±1144	5706±1240	6357±1353	6199±1356	7477±1094	67.9	23.7	15.0	17.1
晚稻LR	173	5313±1220	5829±1282	6627±1413	6490±1335	7590±1368	70.0	23.2	12.7	14.5
单季稻SR	324	5869±1199	6471±1221	7067±1314	7086±1301	8346±1170	70.3	22.5	15.3	15.1
菜后稻RAV	63	6358±1299	6707±1369	7494±1290	7309±1382	8034±1460	79.1	16.5	6.7	9.0
烟后稻RAT	101	6339±813	6684±990	7601±895	7501±986	8462±916	74.9	21.0	10.2	11.4

下载: 导出CSV

| 显示表格

表 5 福建水稻氮磷钾施肥类别最佳分类数及其试验点数

Table 5. Optimal categories of NPK fertilization for rice and field test sites

稻作 Rice	区域 Area	试验点数 Number of trials	施肥类别 Number of fertilization
早稻ER	山区Mountain area	91	7
	沿海Coastal plain	80	6
	全省Province	171	7
晚稻LR	山区Mountain area	128	6
	沿海Coastal plain	45	4
	全省Province	173	7
单季稻SR	山区Mountain area	274	6
	沿海Coastal plain	50	4
	全省Province	324	6
菜后稻RAV	全省Province	63	3
烟后稻RAT	山区Mountain area	101	4

下载: 导出CSV

| 显示表格

表 9 菜后稻和烟后稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 9. Categories and optimal NPK fertilizations for after-vegetable and after-tobacco cropping rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
菜后稻RAV	1	4190±750	5329±420	137	20	117	104	0	59
	2	6179±677	7536±526	155	37	94	118	21	53
	3	7085±1106	9169±850	188	56	136	149	42	81
烟后稻RAT	1	5513±665	7355±394	159	62	100	132	10	73
	2	6308±503	8217±514	149	42	85	122	37	64
	3	6390±464	8831±536	159	65	110	137	51	77
	4	7611±744	9857±648	170	48	103	138	43	74

下载: 导出CSV

| 显示表格

表 10 氮磷钾施肥限量标准对供试水稻产量和施肥效益的影响

Table 10. Effect of optimal NPK fertilizations on rice yield and benefits

稻作 Area	试验点数 Number of trails	无肥区产量 No fertilizer area yield/(kg·hm^-2)	习惯施肥 Farmer practice		推荐施肥 Recommend fertilization		平均增收节支 Average increment and saving
稻作 Area	试验点数 Number of trails	无肥区产量 No fertilizer area yield/(kg·hm^-2)	肥料成本Cost/(Yuan·hm^-2)	产量Yield/(kg·hm^-2)	肥料成本Cost/(Yuan·hm^-2)	产量Yield/(kg·hm^-2)	增产Increase yield/%	肥料成本Cost/(Yuan·hm^-2)	净增收Net income/(Yuan·hm^-2)
早稻RR	52	5629±1162	1506	7210±919	1397	7946±905	10.2	-109	1875
晚稻LR	58	5507±1202	1596	7090±1194	1400	7373±1263	4.0	-196	875
单季稻SCR	111	5975±1173	1695	7328±989	1580	7792±1087	6.3	-115	1228
菜后稻RAV	4	6453±741	1025	7633±1067	698	8587±1322	12.5	-327	2616
烟后稻RAT	9	6814±1031	1161	7654±1129	1002	8094±1268	5.7	-159	1215

下载: 导出CSV

| 显示表格

1.2 氮磷钾施肥类别及其推荐施肥量的确定方法

为兼顾无肥区稻谷产量水平及其氮磷钾施肥对稻谷产量效应的差异，在确定稻田施肥类别时，采用欧氏距离-离差平方和系统聚类方法^[14]确定各个试验点的类别归属，类别之间差异显著性检验采用F检验^{[14, 16]}。因部分施肥类别的试验点数不能满足F值计算式的要求，本研究统一选择“3414”设计的处理(1)和处理(6)两个产量指标作为显著性检验对象，即：确保这两个产量指标具有显著差异。

确定水稻氮磷钾施肥类别后，根据各施肥类别内对应的“3414”试验点，应用笔者提出的三元非结构肥效模型构建该类别的氮磷钾肥料效应函数方程，即：

$Y = A\left( {{{\rm{N}}_0} + {\rm{N}}} \right)\left( {{{\rm{P}}_0} + {\rm{P}}} \right)\left( {{{\rm{K}}_0} + {\rm{K}}} \right){e^{ - {c_1}{\rm{N}} - {c_2}{\rm{P}} - {c_3}{\rm{K}}}}$

其中，Y为作物产量，N、P、K分别为N、P₂O₅、K₂O的施肥量；A表示施肥量X=0时试验地土壤对稻谷产量的生产能力；N₀、P₀、K₀分别表示供试土壤的氮、磷、钾供肥当量，并以N、P₂O₅、K₂O养分形态计量；c₁、c₂、c₃分别表示施用氮、磷、钾养分的增产效应系数。该模型在试验施肥量范围内存在一个产量峰值，该峰值对应的施肥量即为最高产量施肥量。因此，根据边际产量导数法可求得最高施肥量和经济施肥量的计算式。

1.3 氮磷钾施肥类别的系统聚类分析和三元肥效模型构建的计算机实现

文中涉及的欧氏距离-离差平方和法系统聚类分析、聚类谱系图绘制和施肥类别间差异显著性检验、三元非结构肥效模型参数估计和统计显著性检验等，均采用MATLAB软件完成。由于三元非结构肥效模型是非线性模型，模型参数估计采用非线性最小二乘法^[17]。计算的数学原理、计算过程和MATLAB软件相关功能函数的使用方法可参阅相关专著^[16-18]。

2. 结果与分析

2.1 水稻氮磷钾施肥的增产效应

832个水稻氮磷钾田间肥效试验结果按照稻作分别归类，经典五处理的平均稻谷产量见表 4。平均而言，稻田土壤对稻谷贡献率为71.1%；氮磷钾化肥的平均增产率分别为22.2%、13.4%和14.4%，增产效果是N>K >P。就稻作而言，无论是空白区、平衡施肥还是氮磷钾缺素区，单季稻产量均明显高于晚稻相应处理产量，而晚稻产量则高于早稻产量。

以表 4的平均产量和标准差为基础的变异系数计算结果表明，不同区域或土壤的相同稻作产量存在较大的变异系数。因土施肥、因产量水平施肥，对提高施肥效益具有重要意义。

2.2 水稻氮磷钾施肥类别的确定

为消除人为因素对氮磷钾施肥类别划分的影响，以早稻171个试验结果为例，应用欧氏距离-离差平方和法进行系统聚类分析^[14]，然后根据聚类谱系图，将早稻施肥类别划分为3~10个类别，然后进行协方差矩阵差异显著性检验^{[14, 16]}。结果表明，将早稻施肥类别分成8类或8类以上时，至少有一个施肥类别仅包含1个试验点，因试验点数太少而缺乏代表性；同时至少有5对组合的空白区产量和平衡施肥产量的协方差矩阵差异不显著，导致分类结果无效。因此，将施肥类别分成8类或以上显然不合适。将早稻施肥类别分成3~7类时，各个类别内均包含有13个以上试验点，具有较好的代表性；两两施肥类别组合间，空白区产量和平衡施肥产量的协方差矩阵均有显著水平以上的差异。因此，早稻氮磷钾施肥类别最多可划分成7个类别，此时类别1至类别7分别包括27、34、10、33、13、19、35个试验点组合。

徐志平等^[19]研究表明，就早稻、晚稻和单季稻的氮磷钾施肥效应的差异而言，福建稻田可划分为山区和沿海两个区域。因此，将南平市、三明市、龙岩市和宁德市的稻田划分为“山区稻田”，将福州市、莆田市、泉州市、漳州市的稻田划分为“沿海稻田”。福建烟区分布在龙岩市、三明市和南平市，因而烟后稻应归属山区稻田；菜后稻分布于全省各地，但面积小，试验点也少，故不分区域。按照上述施肥类别确定方法，各稻作施肥类别的最佳分类数见表 5。

2.3 基于农学效应的水稻氮磷钾施肥限量标准

根据表 5的区域划分及其施肥类别数，针对各稻作施肥类别内的相关氮磷钾田间肥效试验结果，分别建立该施肥类别的三元非结构肥效模型(结果略)。统计检验表明，60个三元非结构肥效模型的F值均达到统计显著水平，而且都属于典型肥效模型。

因此，根据三元非结构肥效模型，以N 4.3元·kg^-1、P₂O₅ 5元·kg^-1、K₂O 5元·kg^-1和稻谷2.4元·kg^-1的市场均价求得各稻作不同施肥类别的推荐施肥量，结果见表 6~9，即为基于农学效应的福建水稻氮磷钾施肥限量标准。

表 6 早稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 6. Categories and optimal NPK fertilizations for early rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	3914±676	6401±408	131	48	98	111	43	80
	2	4197±336	8382±436	134	58	115	125	55	105
	3	4656±515	6954±271	149	54	112	119	47	91
	4	5221±627	7314±228	173	61	132	144	46	85
	5	5348±787	8254±436	189	49	117	152	45	94
	6	6587±525	8813±383	172	60	143	139	50	99
	7	7544±524	9774±414	192	51	87	155	37	73
沿海Coastal plain	1	3538±589	5661±424	201	67	79	159	58	62
	2	4417±178	6382±322	166	57	101	122	49	86
	3	4985±573	6889±244	172	61	104	136	47	80
	4	5340±704	7579±360	167	65	91	138	49	76
	5	6176±698	8166±95	148	57	93	119	49	75
	6	6809±748	9411±637	147	68	91	127	59	76
全省Province	1	4025±634	6211±471	149	51	94	119	46	78
	2	4596±449	6951±262	165	58	101	134	49	84
	3	4755±716	6785±418	186	62	124	150	53	54
	4	5288±535	7348±281	158	64	106	130	42	86
	5	5597±1010	8148±446	162	60	107	138	52	72
	6	5902±745	8631±496	178	44	146	139	41	113
	7	7125±707	9490±552	163	58	89	136	48	75
注：空白区产量和平衡施肥产量均为试验产量结果的均值±标准差。表 7~10同。 Note:Rice yields at control and balanced fertilization sites are mean±standard deviation.The same as table 7-10.

下载: 导出CSV

| 显示表格

表 7 晚稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 7. Categories and optimal NPK fertilizations for late rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	3399±480	5329±626	171	74	85	139	61	57
	2	4776±5618	6793±376	155	56	118	127	46	81
	3	4934±535	7466±657	186	60	135	154	52	97
	4	5613±667	7912±411	164	70	124	132	57	97
	5	6131±823	9092±662	183	52	115	157	46	87
	6	7815±818	10261±622	179	54	97	145	44	81
沿海Coastal plain	1	4378±415	6299±324	194	51	111	138	44	91
	2	4689±577	6831±268	143	54	89	119	44	74
	3	5691±551	7504±346	148	63	92	120	43	73
	4	6830±428	8943±435	159	48	83	139	37	49
全省Province	1	3473±521	5369±617	179	63	87	143	53	61
	2	4784±652	6802±390	146	58	106	119	48	81
	3	4934±535	7466±657	186	60	135	154	53	97
	4	5253±615	7333±754	179	26	42	154	21	29
	5	5638±685	7921±416	168	66	117	136	53	92
	6	6202±839	9012±445	175	56	138	149	49	97
	7	7565±893	10243±583	186	54	94	155	42	79

下载: 导出CSV

| 显示表格

表 8 单季稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 8. Categories and optimal NPK fertilizations for single cropping rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	4068±631	5949±328	181	51	106	132	43	86
	2	4671±740	7297±421	221	64	109	176	57	87
	3	5735±711	8285±429	194	98	170	161	68	117
	4	6260±547	8515±454	166	47	92	136	42	70
	5	6540±839	9272±441	204	84	162	165	68	115
	6	7594±777	10152±664	191	63	97	153	56	74
沿海Coastal plain	1	4411±653	6447±568	171	71	108	135	58	85
	2	5681±640	7770±367	157	75	105	124	62	81
	3	6445±639	8235±712	156	60	95	117	49	74
	4	7230±107	10761±430	215	70	157	185	57	103
全省Province	1	4269±588	6614±560	225	66	113	173	57	89
	2	5039±757	7482±325	191	63	106	152	56	84
	3	5465±646	8186±494	205	98	159	167	72	114
	4	6194±641	8429±513	162	50	94	133	44	72
	5	6487±707	8925±526	185	74	155	149	61	103
	6	7454±842	10184±685	194	66	115	160	58	87

下载: 导出CSV

| 显示表格

2.4 氮磷钾施肥限量标准对水稻产量和施肥效益的影响

234个大区对比试验结果见表 10。与空白区稻谷产量相比，无论是习惯施肥还是推荐施肥，对水稻均具有显著的增产作用；基于氮磷钾施肥限量标准的推荐施肥技术，因调整了氮磷钾施肥量和养分比例，在习惯施肥基础上进一步提高了水稻产量。不同稻作间平均增产4.0%~12.5%，加权平均增产率为6.7%。

以N 4.3元·kg^-1、P₂O₅ 5元·kg^-1、K₂O 5元·kg^-1和稻谷2.4元·kg^-1的市场均价计算施肥效益。234个大区对比试验加权平均而言，与习惯施肥相比，各稻作均不同程度减少了化肥成本，平均降低化肥成本139元·hm^-2；不同稻作净增收幅度在875~2 616元·hm^-2，平均达到1 309元·hm^-2。在减肥方面，推荐施肥平均降低N 37.5kg·hm^-2和P₂O₅ 16.5 kg·hm^-2，减肥幅度分别达到19.4%和23.9%，但增加K₂O用量27.3%。结果显示，制定的水稻氮磷钾施肥限量标准具有较好的增产增收效果和降低氮磷肥用量的作用。

3. 讨论

3.1 氮磷钾施肥类别的定量分类

在测土配方施肥中，如何对多年多点田间肥效试验结果进行分类、提出不同土壤肥力等级或者施肥目标产量下的推荐施肥技术，一直是个重要研究课题。福建山区(包括南平市、三明市、龙岩市和宁德市)和沿海(包括福州市、莆田市、泉州市和漳州市)两个区域在水稻施肥效应上具有明显差异，如山区稻田土壤对早稻和晚稻产量的平均贡献率分别是71.4%和74.4%, 而闽东南沿海地区则分别是79.6%和75.2%^[19]，山区稻田碱解氮、Olsen-P、速效钾的丰缺指标值均高于沿海稻田^[5]，山区稻田氮磷钾推荐用量大都明显高于沿海稻田^[5-6]。因此，将田间肥效试验资料按照山区和沿海两个区域分类总结具有更好的针对性。因此，本研究按照山区区域、沿海区域以及将全省作为一个整体等3种条件下进行系统聚类分析，并对分类结果的协方差矩阵进行差异显著性检验。该方法的优点是，空白区和平衡施肥的稻谷产量两个指标在不同施肥类别间具有显著水平差异，确保了分类结果的有效性。

在推荐施肥中，在服务区域内的土壤肥力水平或者目标产量应该如何划分等级，即划分成多少个施肥类别，至今仍然没有定论。由于施肥量和作物产量的关系不是确定性关系，施肥类别间在统计学上应该具有显著水平的差异，才能保证分类结果的有效性。迄今不少学者研究提出了许多肥料试验资料的分类或聚类方法^{[9, 20-21]}，具有简洁、实用和易于理解的优点，但都忽略了类别间差异显著性的问题。本文提出的系统聚类分析结合类别间差异显著性检验的研究思路，从定量角度较好地解决了施肥类别划分的问题。

3.2 水稻氮磷钾施肥限量标准的应用

研究表明，在区域水稻推荐施肥建模时，三元二次多项式肥效模型出现了近半数的非典型式，影响了区域推荐施肥的可靠性^[11-12]。三元非结构肥效模型明显提高了典型肥效模型的出现几率，而且较好地克服了二次多项式肥效模型推荐施肥量普遍偏高^[22-23]的问题，推荐施肥量更好地满足减肥增效技术的客观要求。一般情况下，在应用中应首先考虑采用本试验中的经济施肥量，若种植品种增产潜力较大时，可考虑在经济施肥量和最高施肥量之间调整，但不宜超过最高施肥量。

应用系统聚类方法构建水稻氮磷钾施肥类别的另一个优点是，能根据“3414”设计的14个处理的氮、磷、钾的不同施肥量组合及其对应产量的相似程度进行归类合并^[14]，将整体施肥效应有较大区别的试验点区分开来，使同一个施肥类别内的施肥效应具有更好的相似性，克服了现有分类方法^[8-9]仅凭借空白区产量或平衡施肥产量的相近性而划分为同一个类别的缺点。

但是，本研究的结果在直接应用上也带来的另外一个问题，即：各施肥类别的空白区产量及其标准差或者施肥目标产量及其标准差都表明，在95%置信区间下，不同施肥类别之间这两个产量指标会出现一定程度的重叠，给实际应用时的施肥决策带来一定困难。幸运的是，当前统计模式识别技术的成熟应用和信息技术的飞速发展，为解决这一应用问题提供了有效技术途径。

表 1 福建水稻氮磷钾田间肥效试验数量及其试验点设置情况

Table 1 Number and distribution of sites for field-testing rice growth in response to varied NPK fertilizations in Fujian

稻作 Rice	试验类型 Treatments type	试验点数 Number of trials	试验点分布情况Experiment distribution/个
稻作 Rice	试验类型 Treatments type	试验点数 Number of trials	福州市 Fuzhou	莆田市 Putian	泉州市 Quanzhou	漳州市 Zhangzhou	龙岩市 Longyan	三明市 Sanming	南平市 Nanping	宁德市 Ningde
早稻Early rice	3414试验	171	27	36	0	17	36	28	17	10
早稻Early rice	大区对比试验PCE	52	0	14	0	0	17	10	2	8
晚稻Late rice	3414试验	173	18	17	0	10	39	36	23	30
晚稻Late rice	大区对比试验PCE	58	0	3	2	2	24	9	14	3
单季稻Single cropping rice	3414试验	324	25	25	0	0	20	77	90	87
单季稻Single cropping rice	大区对比试验PCE	111	16	15	12	0	8	27	18	15
菜后稻Rice after vegetable	3414试验	63	2	0	0	11	29	9	7	5
菜后稻Rice after vegetable	大区对比试验PCE	4	0	0	2	2	0	0	0	0
烟后稻Rice after tobacco	3414试验	101	0	0	0	0	52	49	0	0
烟后稻Rice after tobacco	大区对比试验PCE	9	0	0	0	0	5	1	3	0
注：PCE表示大区对比试验。 Note: PCE stands for plot check experiment.

下载: 导出CSV

表 2 氮磷钾“3414”肥效试验供试土壤主要理化性状及其处理(6)施肥量

Table 2 Physiochemical properties of soil at test sites with "3414" design and fertilizer applications of treatment 6

稻作 Rice	试验点数 Number of trials	土壤主要理化性状 Soil physical and chemical properties					处理(6)施肥量 Application rate of Treatment 6/(kg·hm^-2)
稻作 Rice	试验点数 Number of trials	pH	有机质 Organic matter/(g·kg^-1)	碱解氮 Alkali-hydr N/(mg·kg^-1)	有效磷 Olsen-P/(mg·kg^-1)	速效钾 Avail K/(mg·kg^-1)	氮肥 N	磷肥 P₂O₅	钾肥 K₂O
早稻ER	171	4.7±0.8	28.4±6.2	158.8±48.7	33.4±23.4	86.3±54.4	158.0±10.3	61.2±8.8	109.2±14.2
晚稻LR	173	4.6±0.8	31.5±7.2	149.1±56.6	32.4±24.7	88.9±55.3	157.0±10.4	59.0±9.3	112.5±40.0
单季稻SR	324	4.6±0.6	28.1±13.4	147.6±60.8	28.0±21.0	87.9±48.9	165.9±16.5	67.4±12.0	114.6±12.9
菜后稻RAV	63	5.4±0.4	32.4±7.7	162.7±41.5	42.8±18.9	97.8±46.7	158.3±18.9	62.0±16.2	107.0±19.8
烟后稻RAT	101	4.9±0.2	28.4±17.2	159.7±43.2	36.1±19.2	124.0±68.1	160.1±18.6	57.9±17.5	105.0±16.5
注：表中土壤理化性状测定值和处理(6)数据以平均值±标准差表示，ER表示早稻, LR表示晚稻，SR表示单季稻, RAV表示菜后稻, RAT表示烟后稻。表 3~5、9~10同。 Note:Soil physiochemical properties and treatment 6 are expressed as mean±standard deviation, ER stands for early rice, LR for late rice, SR for single-cropping rice, RAV for after-vegetable cropping, and RAT for after-tobacco cropping.The same as table 3-5, 9-10.

下载: 导出CSV

表 3 大区对比试验的土壤主要理化性状和试验处理施肥量

Table 3 Physiochemical properties of soil at comparison test sites and fertilizer applications of treatments

稻作 Rice	试验点数 Number of trials	土壤主要理化性状 Soil physical and chemical properties					推荐施肥 Recommend fertilization/(kg·hm^-2)			习惯施肥 Farmer practice/(kg·hm^-2)
稻作 Rice	试验点数 Number of trials	pH	有机质 Organic matter/(g·kg^-1)	碱解氮 Alkali-lhydr. N/(mg·kg^-1)	有效磷 Olsen-P/(mg·kg^-1)	速效钾 Avail K/(mg·kg^-1)	氮肥 N	磷肥 P₂O₅	钾肥 K₂O	氮肥 N	磷肥 P₂O₅	钾肥 K₂O
早稻ER	52	5.3±0.6	27.2±11.5	149.1±45.7	30.7±22.9	67.5±36.0	169±20	59±24	119±31	184±24	72±18	71±36
晚稻LR	58	5.2±0.5	30.9±8.9	157.9±39.6	29.5±17.2	92.6±75.2	154±25	50±15	97±30	195±21	67±24	54±30
单季稻SR	111	5.2±0.4	29.2±10.4	144.1±48.7	34.2±25.4	81.3±55.5	165±17	63±16	110±28	174±37	71±21	61±39
菜后稻RAV	4	4.9±0.5	16.7±7.1	122.2±16.8	100.3±37.8	85.6±44.8	75±86	23±26	53±60	131±103	56±43	19±38
烟后稻RAT	9	5.4±0.6	38.1±9.3	175.1±56.1	29.3±10.3	104.2±53.3	133±29	43±11	51±23	122±58	39±36	48±49

下载: 导出CSV

表 4 不同稻作施用氮磷钾化肥对稻谷产量的影响

Table 4 Effect of NPK fertilizations on yield of rice under varied cultivation practices

稻作 Rice	试验点数 Number of trials	经典五处理及其平均产量 Classical five treatments and the average yield/(kg·hm^-2)					平均贡献率 Average contribution rates/%
稻作 Rice	试验点数 Number of trials	对照CK	N₀P₂K₂	N₂P₀K₂	N₂P₂K₀	N₂P₂K₂	土壤Soil	氮N	磷P	钾K
早稻ER	171	5080±1144	5706±1240	6357±1353	6199±1356	7477±1094	67.9	23.7	15.0	17.1
晚稻LR	173	5313±1220	5829±1282	6627±1413	6490±1335	7590±1368	70.0	23.2	12.7	14.5
单季稻SR	324	5869±1199	6471±1221	7067±1314	7086±1301	8346±1170	70.3	22.5	15.3	15.1
菜后稻RAV	63	6358±1299	6707±1369	7494±1290	7309±1382	8034±1460	79.1	16.5	6.7	9.0
烟后稻RAT	101	6339±813	6684±990	7601±895	7501±986	8462±916	74.9	21.0	10.2	11.4

下载: 导出CSV

表 5 福建水稻氮磷钾施肥类别最佳分类数及其试验点数

Table 5 Optimal categories of NPK fertilization for rice and field test sites

稻作 Rice	区域 Area	试验点数 Number of trials	施肥类别 Number of fertilization
早稻ER	山区Mountain area	91	7
	沿海Coastal plain	80	6
	全省Province	171	7
晚稻LR	山区Mountain area	128	6
	沿海Coastal plain	45	4
	全省Province	173	7
单季稻SR	山区Mountain area	274	6
	沿海Coastal plain	50	4
	全省Province	324	6
菜后稻RAV	全省Province	63	3
烟后稻RAT	山区Mountain area	101	4

下载: 导出CSV

表 9 菜后稻和烟后稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 9 Categories and optimal NPK fertilizations for after-vegetable and after-tobacco cropping rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
菜后稻RAV	1	4190±750	5329±420	137	20	117	104	0	59
	2	6179±677	7536±526	155	37	94	118	21	53
	3	7085±1106	9169±850	188	56	136	149	42	81
烟后稻RAT	1	5513±665	7355±394	159	62	100	132	10	73
	2	6308±503	8217±514	149	42	85	122	37	64
	3	6390±464	8831±536	159	65	110	137	51	77
	4	7611±744	9857±648	170	48	103	138	43	74

下载: 导出CSV

表 10 氮磷钾施肥限量标准对供试水稻产量和施肥效益的影响

Table 10 Effect of optimal NPK fertilizations on rice yield and benefits

稻作 Area	试验点数 Number of trails	无肥区产量 No fertilizer area yield/(kg·hm^-2)	习惯施肥 Farmer practice		推荐施肥 Recommend fertilization		平均增收节支 Average increment and saving
稻作 Area	试验点数 Number of trails	无肥区产量 No fertilizer area yield/(kg·hm^-2)	肥料成本Cost/(Yuan·hm^-2)	产量Yield/(kg·hm^-2)	肥料成本Cost/(Yuan·hm^-2)	产量Yield/(kg·hm^-2)	增产Increase yield/%	肥料成本Cost/(Yuan·hm^-2)	净增收Net income/(Yuan·hm^-2)
早稻RR	52	5629±1162	1506	7210±919	1397	7946±905	10.2	-109	1875
晚稻LR	58	5507±1202	1596	7090±1194	1400	7373±1263	4.0	-196	875
单季稻SCR	111	5975±1173	1695	7328±989	1580	7792±1087	6.3	-115	1228
菜后稻RAV	4	6453±741	1025	7633±1067	698	8587±1322	12.5	-327	2616
烟后稻RAT	9	6814±1031	1161	7654±1129	1002	8094±1268	5.7	-159	1215

下载: 导出CSV

表 6 早稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 6 Categories and optimal NPK fertilizations for early rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	3914±676	6401±408	131	48	98	111	43	80
	2	4197±336	8382±436	134	58	115	125	55	105
	3	4656±515	6954±271	149	54	112	119	47	91
	4	5221±627	7314±228	173	61	132	144	46	85
	5	5348±787	8254±436	189	49	117	152	45	94
	6	6587±525	8813±383	172	60	143	139	50	99
	7	7544±524	9774±414	192	51	87	155	37	73
沿海Coastal plain	1	3538±589	5661±424	201	67	79	159	58	62
	2	4417±178	6382±322	166	57	101	122	49	86
	3	4985±573	6889±244	172	61	104	136	47	80
	4	5340±704	7579±360	167	65	91	138	49	76
	5	6176±698	8166±95	148	57	93	119	49	75
	6	6809±748	9411±637	147	68	91	127	59	76
全省Province	1	4025±634	6211±471	149	51	94	119	46	78
	2	4596±449	6951±262	165	58	101	134	49	84
	3	4755±716	6785±418	186	62	124	150	53	54
	4	5288±535	7348±281	158	64	106	130	42	86
	5	5597±1010	8148±446	162	60	107	138	52	72
	6	5902±745	8631±496	178	44	146	139	41	113
	7	7125±707	9490±552	163	58	89	136	48	75
注：空白区产量和平衡施肥产量均为试验产量结果的均值±标准差。表 7~10同。 Note:Rice yields at control and balanced fertilization sites are mean±standard deviation.The same as table 7-10.

下载: 导出CSV

表 7 晚稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 7 Categories and optimal NPK fertilizations for late rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	3399±480	5329±626	171	74	85	139	61	57
	2	4776±5618	6793±376	155	56	118	127	46	81
	3	4934±535	7466±657	186	60	135	154	52	97
	4	5613±667	7912±411	164	70	124	132	57	97
	5	6131±823	9092±662	183	52	115	157	46	87
	6	7815±818	10261±622	179	54	97	145	44	81
沿海Coastal plain	1	4378±415	6299±324	194	51	111	138	44	91
	2	4689±577	6831±268	143	54	89	119	44	74
	3	5691±551	7504±346	148	63	92	120	43	73
	4	6830±428	8943±435	159	48	83	139	37	49
全省Province	1	3473±521	5369±617	179	63	87	143	53	61
	2	4784±652	6802±390	146	58	106	119	48	81
	3	4934±535	7466±657	186	60	135	154	53	97
	4	5253±615	7333±754	179	26	42	154	21	29
	5	5638±685	7921±416	168	66	117	136	53	92
	6	6202±839	9012±445	175	56	138	149	49	97
	7	7565±893	10243±583	186	54	94	155	42	79

下载: 导出CSV

表 8 单季稻不同区域的氮磷钾施肥类别及其推荐施肥量

Table 8 Categories and optimal NPK fertilizations for single cropping rice at different sites

区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	最高施肥量 Max. fertilizer rate/(kg·hm^-2)			经济施肥量 Eco. fertilizer rate/(kg·hm^-2)
区域 Area	施肥类别 Number of fertilization types	空白产量 CK yield/(kg·hm^-2)	目标产量 Target yield/(kg·hm^-2)	氮肥N	磷肥P₂O₅	钾肥K₂O	氮肥N	磷肥P₂O₅	钾肥K₂O
山区Mountain area	1	4068±631	5949±328	181	51	106	132	43	86
	2	4671±740	7297±421	221	64	109	176	57	87
	3	5735±711	8285±429	194	98	170	161	68	117
	4	6260±547	8515±454	166	47	92	136	42	70
	5	6540±839	9272±441	204	84	162	165	68	115
	6	7594±777	10152±664	191	63	97	153	56	74
沿海Coastal plain	1	4411±653	6447±568	171	71	108	135	58	85
	2	5681±640	7770±367	157	75	105	124	62	81
	3	6445±639	8235±712	156	60	95	117	49	74
	4	7230±107	10761±430	215	70	157	185	57	103
全省Province	1	4269±588	6614±560	225	66	113	173	57	89
	2	5039±757	7482±325	191	63	106	152	56	84
	3	5465±646	8186±494	205	98	159	167	72	114
	4	6194±641	8429±513	162	50	94	133	44	72
	5	6487±707	8925±526	185	74	155	149	61	103
	6	7454±842	10184±685	194	66	115	160	58	87

下载: 导出CSV

参考文献(23)

[1]	李书田, 刘晓永, 何萍.当前我国农业生产中的养分需求分析[J].植物营养与肥料学报, 2017, 23(6):1416-1432. http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201706002 LI S T, LIU X Y, HE P. Analyses on nutrient requirements in current agriculture production in China[J].Plant Nutrition and Fertilizer Science, 2017, 23(6):1416-1432.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201706002
[2]	宋晓梅, 蒋太明, 刘洪斌, 等.基于"3414"田间试验的水稻施肥模型的研究——以重庆市合川区为例[J].西南大学学报(自然科学版), 2010, 32(9):94-99. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xnnydxxb201009020 SONGX M, JIANGT M, LIUH B, et al. Research on a Rice Fertilization Model Based on the"3414" Field Trial[J]. Journal of Southwest Agricultural University(Natural Science Edition), 2010, 32(9):94-99.(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xnnydxxb201009020
[3]	戢林, 张锡洲, 李廷轩.基于"3414"试验的川中丘陵区水稻测土配方施肥指标体系构建[J].中国农业科学, 2011, 44(1):84-92. DOI: 10.3864/j.issn.0578-1752.2011.01.010 JI L, ZHANG X Z, LI T X. Establishing fertilization recommendation index of paddy soil based on the "3414" field experiments in the middle of Sichuan hilly regions[J]. Scientia Agricultura Sinica, 2011, 44(1):84-92.(in Chinese) DOI: 10.3864/j.issn.0578-1752.2011.01.010
[4]	毛伟, 李文西, 唐宝国, 等.县级测土配方施肥指标体系建立研究——以江苏省江都市水稻为例[J].植物营养与肥料学报, 2014, 20(2):396-406. http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201402016 MAO W, LI W X, TANG B G, et al. Study on index system for soil testing and formulated fertilization-A case study of rice of Jiangdu City in Jiangsu Province[J].Plant Nutrition and Fertilizer Science, 2014, 20(2):396-406.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201402016
[5]	李娟, 章明清, 孔庆波, 等.福建早稻测土配方施肥指标体系研究[J].植物营养与肥料学报, 2010, 16(4):938-946. http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201004024 LI J, ZHANG M Q, KONG Q B, et al. Soil testing and formula fertilization index for early rice in Fujian Province[J]. Plant Nutrition and Fertilizer Science, 2010, 16(4):938-946.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201004024
[6]	孔庆波, 章明清, 李娟, 等.福建省晚稻测土配方施肥指标体系研究[J].热带作物学报, 2010, 31(11):1873-1880. DOI: 10.3969/j.issn.1000-2561.2010.11.001 KONG Q B, ZHANG M Q, LI J, et al. The study of soil test and formula fertilization index system on late rice in Fujian[J].Chinese Journal of Tropical Crops, 2010, 31(11):1873-1880.(in Chinese) DOI: 10.3969/j.issn.1000-2561.2010.11.001
[7]	李娟, 章明清, 姚宝全, 等.福建单季稻氮磷钾推荐施肥量研究[J].福建农业学报, 2015, 30(10):933-938. DOI: 10.3969/j.issn.1008-0384.2015.10.002 LI J, ZHANG M Q, YAO B Q, et al. Recommended N, P and K Fertilization for Single-cropping Rice in Fujian[J]. Fujian Journal of Agricultural Sciences, 2015, 30(10):933-938.(in Chinese) DOI: 10.3969/j.issn.1008-0384.2015.10.002
[8]	金耀青, 张中原.配方施肥方法及其应用[M].沈阳:辽宁科学技术出版社, 1993:1-64. JIN Y Q, ZHANG Z Y. Formula Fertilization Method and Its Application[M]. Shenyang:Liaoning Science and Technology Press, 1993:1-64.(in Chinese)
[9]	杨守春, 陈伦寿, 毛达如.黄淮海平原主要作物优化施肥和土壤培肥研究总论.黄淮海平原主要作物优化施肥和土壤培肥技术[M].北京:中国农业科技出版社, 1991:1-26. YANG S C, CHEN L S, MAO D R. Study on Main Crops Optimal Fertilization and Soil Fertilization in the Huanghuaihai Plain[M]. Beijing:China Agricultural Science and Technology Press, 1991:1-26.(in Chinese)
[10]	李仁岗.肥料效应函数[M].北京:农业出版社, 1987:1-42. LI R G. Fertilizer Response Function[M]. Beijing:Agricultural Press, 1987:1-42.(in Chinese)
[11]	章明清, 李娟, 孔庆波, 等.作物肥料效应函数模型研究进展与展望[J].土壤学报, 2016, 53(6):1143-1156. http://d.old.wanfangdata.com.cn/Periodical/trxb201606001 ZHANG M Q, LI J, KONG Q B, et al. Progress and Prospect of the Study on Crop-response-to-Fertilization Function Model[J]. Acta Pedologica Sinica, 2016, 53(6):1143-1156.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/trxb201606001
[12]	MONTGOMERY D C, PECK E A, VINING G G.线性回归分析导论[M].王辰勇, 译.第5版.北京: 机械工业出版社, 2016: 91-268. MONTGOMERY D C, PECK E A, VINING G G. Introduction of Linear Regression Analysis[M].Translation: WANG C Y. 5^th Ed. Beijing: China Machine Press, 2016: 91-268.(in Chinese)
[13]	ZHANG M Q, LI J, CHEN F, et al. Unary Non-structural Fertilizer Response Model for Rice Crops and Its Field Experimental Verification[J]. Scientific Reports, 2018, 8(11):2792. http://cn.bing.com/academic/profile?id=6a2c5ede8ed234106fb96e529c3c0f81&encoded=0&v=paper_preview&mkt=zh-cn
[14]	李娟, 章明清, 孔庆波, 等.构建县域早稻氮磷钾施肥的系统聚类方法研究[J].植物营养与肥料学报, 2017, 23(2):531-538. http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201702030 LI J, ZHANG M Q, KONG Q B, et al. Building fertilization categories of N, P and K fertilization for early rice using systematic clustering method in county territory[J]. Plant Nutrition and Fertilizer Science, 2017, 23(2):531-538.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/zwyyyflxb201702030
[15]	鲁如坤.土壤农业化学分析方法[M].北京:中国农业科技出版社, 2000:146-196. RU R K. Soil Agrochemical Analysis Methods[M]. Beijing:China Agricultural Science and Technology Press, 2000:146-196.(in Chinese)
[16]	袁志发, 宋世德.多元统计分析(第2版)[M].北京:科学出版社, 2009, 278-293. YUAN Z F, SONG S D. Multivariate Statistical Analysis (2^nd Edition)[M]. Beijing:Science Press, 2009, 278-293.(in Chinese)
[17]	马昌凤, 柯艺芬, 谢亚君.优化计算方法及其MATLAB程序实现[M].北京:国防工业出版社, 2015:31-45. MA C F, KE Y F, XIE Y J. Optimization Calculation Method and Realization of MATLAB Program[M]. Beijing:National Defense Industry Press, 2015:31-45.(in Chinese)
[18]	薛毅, 陈立萍.实用数据分析与MATLAB软件[M].北京:北京工业大学出版社, 2015:189-274, 377-394. XUE Y, CHEN L P. Practical Data Analysis and MATLAB Software[M]. Beijing:Beijing University of Technology Press, 2015:189-274, 377-394.(in Chinese)
[19]	徐志平, 姚宝全, 章明清, 等.福建主要粮油作物测土配方施肥指标体系研究, Ⅰ.土壤基础肥力对作物产量的贡献率及其施肥效应[J].福建农业学报, 2008, 23(4):396-402. DOI: 10.3969/j.issn.1008-0384.2008.04.014 XU Z P, YAO B Q, ZHANG M Q, et al. Soil testing and fertilizer formulation for major crops in Fujian-effects of soil condition and fertilization on crop yield[J]. Fujian Journal of Agricultural Sciences, 2008, 23(4):396-402.(in Chinese) DOI: 10.3969/j.issn.1008-0384.2008.04.014
[20]	王兴仁, 陈伦寿, 毛达如, 等.分类回归综合法及其在区域施肥决策中的应用[J].土壤通报, 1989, 20(1):17-21. http://www.cnki.com.cn/Article/CJFDTotal-TRTB198901005.htm WANG X R, CHEN R S, MAO D R, et al. Classification Regression Synthesis and Its Application in Regional Fertilization Decision[J].Chinese Journal of Soil Science, 1989, 20(1):17-21.(in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-TRTB198901005.htm
[21]	毛达如, 张承东.多点肥料效应函数的动态聚类方法[J].北京农业大学学报, 1991, 17(2):49-54. http://www.cnki.com.cn/Article/CJFDTOTAL-NYDX199102011.htm MAO D R, ZHANG C D. Dynamic clustering method of multi-point fertilizer effect function[J]. Journal of Beijing Agricultural University, 1991, 17(2):49-54.(in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-NYDX199102011.htm
[22]	王兴仁, 陈新平, 张福锁, 等.施肥模型在我国推荐施肥中的应用[J].植物营养与肥料学报, 1998, 4(1):67-74. DOI: 10.3321/j.issn:1008-505X.1998.01.011 WANG X R, CHEN X P, ZHANG F S, et al. Application of Fertilization Model for Fertilize Recommendation in China[J]. Plant Nutrition and Fertilizer Science, 1998, 4(1):67-74.(in Chinese) DOI: 10.3321/j.issn:1008-505X.1998.01.011
[23]	陈新平, 周金池, 王兴仁, 等.小麦-玉米轮作制中氮肥效应模型的选择——经济和环境效益分析[J].土壤学报, 2000, 37(3):346-354. DOI: 10.3321/j.issn:0564-3929.2000.03.008 CHEN X P, ZHOU J C, WANG X R, et al. Economic and environmental evaluation on models for describing crop yield response to nitrogen fertilizers at winter-wheat and summer-corn rotation system[J]. Acta Pedologica Sinica, 2000, 37(3):346-354.(in Chinese) DOI: 10.3321/j.issn:0564-3929.2000.03.008

施引文献(3)

期刊类型引用(3)

1.	章明清，李娟，张立成，姚宝全，张华. 三元肥料效应模型的整合与优化建模策略. 土壤学报. 2021(03): 755-766 . 百度学术
2.	韦玲玲. 融水县中晚稻最佳施肥量研究. 现代农业科技. 2020(15): 3-4 . 百度学术
3.	陈兴麟，刘洋，罗土炎，王昌毓，杨海宁，曾润颖. 海藻寡糖生物肥对水稻产量及土壤菌群的影响. 福建农业科技. 2020(09): 19-25 . 百度学术

其他类型引用(0)

资源附件(0)

表(10)

计量

文章访问数: 1897
HTML全文浏览量: 351
PDF下载量: 55
被引次数: 3

0. 引言
1. 材料与方法
1.1 水稻氮磷钾田间肥效试验
1.2 氮磷钾施肥类别及其推荐施肥量的确定方法
1.3 氮磷钾施肥类别的系统聚类分析和三元肥效模型构建的计算机实现
2. 结果与分析
2.1 水稻氮磷钾施肥的增产效应
2.2 水稻氮磷钾施肥类别的确定
2.3 基于农学效应的水稻氮磷钾施肥限量标准
2.4 氮磷钾施肥限量标准对水稻产量和施肥效益的影响
3. 讨论
3.1 氮磷钾施肥类别的定量分类
3.2 水稻氮磷钾施肥限量标准的应用

0. 引言
1. 材料与方法
1.1 水稻氮磷钾田间肥效试验
1.2 氮磷钾施肥类别及其推荐施肥量的确定方法
1.3 氮磷钾施肥类别的系统聚类分析和三元肥效模型构建的计算机实现
2. 结果与分析
2.1 水稻氮磷钾施肥的增产效应
2.2 水稻氮磷钾施肥类别的确定
2.3 基于农学效应的水稻氮磷钾施肥限量标准
2.4 氮磷钾施肥限量标准对水稻产量和施肥效益的影响
3. 讨论
3.1 氮磷钾施肥类别的定量分类
3.2 水稻氮磷钾施肥限量标准的应用

参考文献(23)

施引文献(3)

资源附件(0)

基于农学效应的福建省水稻氮磷钾施肥限量标准

作者简介: 李娟(1977-), 女, 硕士, 副研究员, 研究方向:作物营养与施肥(E-mail:li-95@163.com)

通讯作者: 章明清(1963-), 男, 博士, 研究员, 研究方向:作物施肥原理与技术(E-mail:zhangmq@163.com)

计量

出版历程