全 文 ：收稿日期：!""#$%%$"! 接受日期：!""&$"’$%’
基金项目：国家科技支撑计划（!""()*+%")",，!""()*+%")"&）；“’-&”项目（!""($.("）；测土配方施肥项目资助。
作者简介：孙义祥（%’#%—），男，安徽省霍邱县人，副研究员，博士研究生，研究方向为养分资源综合管理。/01234：567839327:;%!(< =>1
! 通讯作者 /01234：=63?4; =26@ AB6@ =7；=CA79D;=26@ AB6@ =7
应用“!"#"”试验建立冬小麦测土配方施肥
指标体系
孙义祥%，!，郭跃升,，于舜销章,，蒋庆功,，程琳琳%，崔振岭%!，
陈新平%!，江荣风%，张福锁%
（%中国农业大学资源与环境学院，北京 %"""’-；! 安徽省农业科学院土壤肥料研究所，安徽合肥 !,"",%；
, 山东省土壤肥料总站，山东济南 !E"%""）
摘要：以山东省冬小麦为例，对应用“,-%-”完全或部分实施试验结果建立磷、钾施肥技术指标体系的技术环节进行
探讨，并就当前建立指标体系中遇到的问题进行解析。本研究以相对产量 #EF、’"F和 ’EF为标准将山东省土壤
有效磷、速效钾划分为低、中、高和极高 -级；并分别用三元二次、一元二次和线性加平台模型对不同土壤养分分级
范围内施肥量与产量关系进行模拟，计算最佳肥料用量。结果表明，当山东省土壤有效磷、速效钾含量处于低等级
（G45A70H I %" 1: J K:、LM-G*=0N I E" 1: J K:）时，磷、钾肥（H!GE和 N!G）用量分别应为 %""!%,"和 %!"!%E" K: J C1!；
中等级（G45A70H %"!," 1: J K:、LM-G*=0N E"!%"" 1: J K:）时，磷、钾肥用量为 &"!%""和 %""!%!" K: J C1!；高等
级（G45A70H ,"!E" 1: J K:、LM-G*=0N %""!%-" 1: J K:）时，磷、钾肥用量为 ("!&"和 ("!&" K: J C1!；当土壤有效
磷高于 E" 1: J K:，速效钾高于 %-" 1: J K: 时，无需磷钾肥施用。
关键词：肥力分级指标；相对产量；推荐施肥
中图分类号：O%-#!"#$%&’"(’)* +(,"+(,-." $)/ +,#$""’.0 12-#’&’3$#’,) -24,002)/$#’,)
’)/25 %$"2/ ,) #(2“,-%-”1’2&/ 25+2-’02)#"
OQL R309327:%，!，.QG R6A05CA7:,，RQ OC670?C27:,，ST*L. U37:0:>7:,，VM/L. W370437%，
VQT XCA70437:%!，VM/L P370D37:%!，ST*L. Y>7:0ZA7:%，XM*L. [6056>%
（! "#$$%&% #’ (%)#*+,% -./ 0.12+#.3%.4-$ 5,2%.,%)，"62.- 7&+2,*$4*+% 8.21%+)249，:%2;2.& !<<<=>，"62.-；? 5#2$ -./ @%+42$2A%+ B.)424*4%，
7.6*2 7,-/%39 #’ 7&+2,*$4*+% 5,2%.,%)，C%’%2 ?D<6%"#-$4#：YA=A7\48，761A]>65 Z3A4B A9DA]31A7\5 2]A =>7B6=\AB 37 72\3>724 5>34 \A5\37: 27B ZA]\343?2\3>7 ]A=>11A7B2\3>7
D]>^A=\，C>_A‘A]，C>_ \> 65A \CA5A A9DA]31A7\24 B2\2 \> A5\2a435C ZA]\343?2\3>7 ]A=>11A7B2\3>7 37BA9 35 6]:A7\48 7AABAB@
b2K37: \CA _37\A] _CA2\“,-%-”A9DA]31A7\5 37 OC27B>7: H]>‘37=A 25 27 A921D4A，\CA D]>=A55 >Z A5\2a435C37: ZA]\343?2\3>7
]A=>11A7B2\3>7 37BA9 _25 37‘A5\3:2\AB 27B B35=655AB@ O>34 G45A70H 27B LM-G*=0N _A]A =42553Z3AB 25 4>_（ I #EF ]A420
\3‘A 83A4B），1AB361（#EF $ ’"F ]A42\3‘A 83A4B），C3:C（’"F $ ’EF ]A42\3‘A 83A4B）27B ‘A]8 C3:C（ c ’EF ]A42\3‘A
83A4B）a25AB >7 \CA ]A42\3‘A 83A4B5 @ [A]\343?A] ]A=>11A7BAB ]2\A5 _A]A 531642\AB a8 \CA 1>BA45 >Z \C]AA0Z2=\>] >] >7A0Z2=0
\>] Z>] A2=C“,-%-”Z3A4B A9DA]31A7\ @ bCA ]A564\5 37B3=2\A \C2\ DC>5DC>]65 ZA]\343?A] ]2\A5（H!GE）2]A %""$%,"，&"$%""，
("$&"，27B " K: J C2 _CA7 5>34 ZA]\343\8 2]A =42553Z3AB 25 4>_，1AB361，C3:C 27B ‘A]8 C3:C（G45A70H I %"，%"$,"，,"$
E" 27B c E" 1: J K:），]A5DA=\3‘A48，_C34A D>\255361 ZA]\343?A] ]2\A5（N!G）2]A %!"$%E"，%""$%!"，("$&"，27B "
K: J C2（LM-G*=0N I E"，E"$%""，%""$%-" 27B c %-" 1: J K:）@ bCA]AZ>]A，A5\2a435C37: ZA]\343?2\3>7 ]A=>11A7B2\3>7
37BA9 a25AB >7“,-%-”Z3A4B A9DA]31A7\24 B2\2 35 D]2=\3=2a4A @
植物营养与肥料学报 !""’，%E（%）：%’#$!",
""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
H427\ L6\]3\3>7 27B [A]\343?A] O=3A7=A
!"# $%&’(：!"#$%&!# ’()**%+%’)#%,! %!-&.；$&()#%/& 0%&(-；+&$#%(%1)#%,! $&’,22&!-)#%,!
土壤肥力指标法是国际上最通用的测土配方施
肥方法［345］，也是当前全国测土配方施肥工作中主
推的施肥方法之一。在 67世纪 87年代我国第二次
土壤普查工作中，土壤养分丰缺指标协作组曾系统
地建立了我国主要土壤类型、主要大田作物的土壤
肥力分级指标［948］，这为测土结果的应用、对分析结
果的校验与解释和确定施肥量等方面奠定了基
础［:437］，在很大程度上推动了当时配方施肥技术的
发展。
自第二次土壤普查以来，由于缺少全国性田间
试验结果的积累，我国测土配方施肥工作一直沿用
以往的土壤养分丰缺指标。然而，我国大田作物品
种特性、产量水平、栽培方式、农民施肥方式和土壤
肥力等要素都发生了很大变化［33436］，原有的技术指
标体系已不能适应当前生产的需求［35］；而且，当前
测土配方施肥目标已由以往单一的追求高产向作物
高产、高效、土壤培肥、环境保护等多目标过渡。因
此，建立我国当前不同农业生态区、不同作物的土壤
肥力分级指标及相应的推荐施肥量是做好我国测土
配方施肥工作的关键。
随着全国测土配方施肥项目实施规模的不断增
加，全国不同农业生态区、不同作物体系积累了大量
的“5939”田间试验数据，而且数据量呈逐年递增趋
势，这为建立我国当前土壤肥力分级指标体系提供
了很好的机遇。但受技术力量的限制，很多科技人
员，特别是基层技术人员对利用“5939”田间试验数
据建立土壤肥力分级指标的技术环节尚不清楚，实
际操作过程遇到种种问题，限制了该项工作的快速
推进。为此，以山东省冬小麦为例，对应用“5939”完
全或部分实施方案试验结果建立磷、钾推荐施肥技
术指标体系的技术环节进行了系统探讨，并就当前
建立指标体系建立过程中遇到的问题进行解析，旨
在评价建立指标体系的必要性和能否利用“5939”试
验建立磷、钾推荐施肥技术指标体系，为全国测土配
方施肥工作中指标体系的建立提供技术支持。
) 材料与方法
)*) 试验设计
试验数据来源于山东省测土配方项目实施县
677;年冬小麦试验，分为“5939”完全实施和部分实
施两种试验方案。本研究共涉及田间试验 3<8 个，
其中完全实施试验 ==个，部分实施试验 336 个，分
布在成武、聊城、广饶、乐陵、寿光、菏泽、无棣、文登、
招远、新泰、枣庄、临淄、章丘、郯城等地。
“5939”完全实施试验，即：氮、磷、钾 5因素，每
个因素 9水平，共 39 个处理的肥料试验设计方案。
9个水平的含义：7水平指不施肥，6水平指当地最
佳施肥量的近似值，3水平 > 6水平 ? 7@;，5水平 > 6
水平 ? 3@;（该水平为过量施肥水平）。39 个处
理［3943=］分 别 为： 3） A7B7C7； 6） A7B6C6； 5）
A3B6C6；9）A6B7C6；;）A6B3C6；=）A6B6C6；<）
A6B5C6；8）A6B6C7；:）A6B6C3；37）A6B6C5；33）
A5B6C6；36）A3B3C6；35）A3B6C3；39）A6B3C3。
“5939”部分实施试验设计有：AB、AC、BC和参
数试验几种方案。其中 AB试验包括处理 6、5、9、;、
=、<、33、36；AC试验包括处理 6、5、=、8、:、37、33、35；
BC试验包括处理 9、;、=、<、8、:、37、39；参数试验包
括处理 3、6、9、=和 8。
)*+ 数据统计方法
3@6@3 土壤磷钾肥力分级指标的建立方法 冬小
麦缺素区相对产量具体计算方法如下：
设施用磷肥的 9 个水平分别为 B7、B3、B6、B5
（氮、钾均为 6个水平）；对应产量分别为 DB7、DB3、
DB6、DB5。如果 DB7最大，缺磷相对产量为 377E；
如果 DB3、DB6或 DB5中之一最大，缺磷相对产量为
DB7与最高产量比值。
用同样的方法计算冬小麦缺钾区相对产量。
以“对数”类型获得相对产量与对应土壤养分测
试值之间的数学关系式，分别以相对产量为 ;7E、
<;E、:7E和 :;E计算对应的土壤养分含量，根据
这些值划分土壤养分丰缺指标。
本文定义冬小麦缺素区相对产量低于 ;7E的
养分值为极低，介于 ;7E至 <;E范围内的为低，介
于 <;E至 :7E范围内的为中，介于 :7E!:;E之间
的为高，大于 :;E的为极高。
3@6@6 推荐施肥量的计算方法 对山东省所有
“5939”试验数据分别用三元二次、一元二次和线性
加平台（一元）的模型模拟，根据散点图趋势和不同
方程拟合的决定系数选择最适模型，三元二次、一元
二次模型可通过边际效应分析确定每个试验点的最
佳磷、钾肥施用量，线性加平台模型可以直接计算最
佳施肥量。
采用三元二次肥料效应模型进行拟合时，所采
用的方程为：
8:3 植 物 营 养 与 肥 料 学 报 3;卷
! " #$ % #&’& % #(’&( % #)’( % #*’(( % #+’) %
#,’)( % #-’&’( % #.’&’) % #/’(’)
式中，’&、’(、’)分别为 0、1(2+、3(2用量。如果上述
方程拟合成功（二次项前系数为负值，一次项前系数
为正值，4检验显著），根据边际收益等于边际成本，
即 5!·6! " 5’·6’ 的原则计算经济最佳产量施肥
量［&*］，分别以 ’&、’(、’) 为变量，对方程两边求导，得
到方程组：
#& % (#(’& % #-’( % #.’) " 1’& 7 1!
#) % (#*’( % #-’& % #/’) " 1’( 7 1!
#+ % (#,’) % #.’& % #/’( " 1’)
{
7 1!
式中，1’&、1’(、1’) 和 1!分别代表 0、1(2+、3(2和冬
小麦的价格，三元二次方程系数 #&、#(、#)、#*、#+、#,、
#-、#.、#/ 的值和肥料 0、1(2+、3(2的价格代入上述
方程组，解方程组即可获取最佳氮、磷、钾肥用量。
采用一元二次模型进行拟合时，方程为：
! " 8 % #’ % 9’(
式中，! 为子粒产量（ :; 7 <=(）； ’ 为肥料用量
（:; 7 <=(）；8为截距；#为一次回归系数；9 为二次
回归系数。选用处理 *、+、,、-的产量结果模拟磷肥
的推荐用量，选用处理 ,、.、/、&$的产量结果模拟钾
肥的推荐用量，根据边际收益等于边际成本，即 5!·
6! " 5’·6’计算经济最佳施肥量。
采用直线加平台模型进行拟合时，方程为：
! " 8 % #’（’!>）
! " 1（’ ? >）
式中，! 为子粒产量（ :; 7 <=(）； ’ 为肥料用量
（:; 7 <=(）；8为截距；#为回归系数；>为直线与平
台的交点，即最佳施肥量；1为平台产量（:; 7 <=(），
即最高子粒产量。
对山东省所有不完全试验数据用一元二次和线
性加平台（一元）的模型模拟。
计算时 1和 3肥以 1(2+和 3(2计，价格分别为
*@+和 )@$元 7 :;，冬小麦为 &@* 元 7 :;。当试验结果
表明施肥增产效果不显著时，则推荐施肥量为 $，如
果推荐施肥量高于试验的最高施肥量，则以试验的
最高施肥量为推荐施肥量。
将每个试验点的最佳施肥量与土壤养分含量对
应进行统计（磷肥用量对应土壤有效磷含量，钾肥用
量对应土壤交换性钾含量），按上述建立的土壤肥力
分级指标进行分类（低、中、高、极高）汇总统计，确定
在不同养分水平下适宜的推荐用量。
! 结果与分析
!"# 山东省土壤有效磷钾分级指标
(@&@& 山东省土壤有效磷分级指标 根据山东省
冬小麦试验缺磷处理的相对产量与土壤有效磷含量
作散点图（图 &），并作对数回归方程：! " &&@*()AB
（’）% +$@-.+（C( " $@(.$)""）。土壤有效磷含量
与缺磷处理冬小麦相对产量的相关性达到极显著性
水平，将相对产量 +$D、-+D、/$D和 /+D代入对数
方程，求出对应的土壤有效磷含量数值分别为 &、.、
)&和 *. =; 7 :;，即为土壤有效磷的丰缺指标值。
所有试验点冬小麦相对产量的平均值为 .,D，
最大值为 &$$D，最小值为 +.D。也就是说，在本研
究条件下缺磷处理冬小麦相对产量没有低于 +$D
的。在本研究中冬小麦相对产量为 +$D时对应的
土壤有效磷含量（& =; 7 :;）为模型外推值，没有实际
意义。因此，本研究将土壤有效磷的 +级指标简化
为 * 级，即冬小麦相对产量低于 -+D时，土壤有效
磷为低等级（. =; 7 :;）；相对产量在 -+D至 /$D之
间，土壤有效磷为中等级（.!)& =; 7 :;）；相对产量
在 /$D至 /+D之间，土壤有效磷为高等级（)&!*.
=; 7 :;）；相对产量大于 /+D时，土壤有效磷为极高
等级（ ? *. =; 7 :;）。考虑到实际操作的方便性，可
将土壤有效磷丰缺指标值进一步简化，即土壤有效
磷含量小于 &$ =; 7 :;的土壤为低等级磷肥力土壤，
介于 &$!)$ =;:; 的土壤为中等级磷肥力土壤，介
于 )$!+$ =; 7 :;的土壤为高等级磷肥力土壤，大于
+$ =; 7 :;的土壤为极高等级磷肥力土壤。
图 # 土壤有效磷含量与缺磷处理冬小麦相对产量的关系
$%&’# ()* +*,-.%/01)%2 3*.4**0 1/%, 5,1*067 -08
+*,-.%9* &+-%0 :%*,8 /; 4%0.*+ 4)*-.
(@&@( 土壤速效钾丰缺指标的建立 山东省冬小麦
试验缺钾处理的相对产量与土壤速效钾含量的对数
//&&期 孙义祥，等：应用“)*&*”试验建立冬小麦测土配方施肥指标体系
回归 方 程 为： ! " #$%&’()*（ +）, &$%-./（0& "
.%1-$1!!）（图 &）。缺钾处理冬小麦相对产量与土壤
速效钾含量的相关性达到极显著性水平，根据方程计
算，获得缺钾区冬小麦相对产量 ’.2、3’2、/.2和
/’2时对应的土壤速效钾含量分别为 (、1$、/3和 #13
45 6 75，这组数值即为土壤速效钾的丰缺指标值。
与土壤有效磷丰缺指标建立遇到相似的问题，
山东省小麦田间试验中低土壤速效钾的试验点很少，
图 ! 土壤速效钾含量与缺钾处理冬小麦相对产量的关系
"#$%! &’( )(*+,#-./’#0 1(,2((. /-#* 3456789: +.;
)(*+,#<( $)+#. =#(*; -> 2#.,() 2’(+,
全部试验中没有一个缺钾处理的小麦相对产量低于
’.2。因此，’.2小麦相对产量对应的土壤速效钾
含量为方程外推获得，具有很大不确定性。考虑到
实际操作的方便性和土壤速效钾含量现状，可以将
山东省的土壤钾素肥力指标分为 $级，速效钾含量
小于 ’. 45 6 75的土壤为低等级钾肥力土壤，速效钾
含量介于 ’.!#.. 45 6 75的土壤为中等级钾肥力土
壤，#..!#$. 45 6 75的土壤为高等级钾肥力土壤，速
效钾含量大于 #$. 45 6 75的土壤为极高等级钾肥力
土壤。
!?! 磷钾肥推荐施肥指标体系的建立
对山东“1$#$”完全实施试验进行三元二次回归
分析发现，模拟成功率仅有 /2，磷钾分别进行一元
二次拟合的成功率为 132和 112（表 #）。
对于三元二次回归方程拟合不成功的试验点分
别采用一元二次方程和线性加平台方程进行拟合，
方程的选择根据散点图趋势和不同方程拟合的决定
系数选择。对所有方程拟合不成功且增产效果不明
显的点，本研究设定最佳施肥量为 .；对于推荐施肥
量高于试验的最高施肥量的点，设定试验的最高施
肥量为推荐施肥量。
表 @ 山东省 AA个“B5@5”试验方程拟合的成功率比较
&+1*( @ &’( )+,#- -> +<+#*+1*( (C0()#D(.,/ 8+*8E*+,(; 1= ;#>>()(., >(),#*#F() D-;(*/
项目
89:4
三元二次方程
;<=::>?@A9B= 4BC:)
一元二次方程
D*:>?@A9B= 4BC:)
E F
模型拟合成功的百分数（2）
;<: G:=A:*9@5: B? HIAA:HH?I) 4BC:)H /%. 13%1 1&%-
模型拟合不成功的百分数（2）
;<: G:=A:*9@5: B? I*>HIAA:HH?I) 4BC:)H /#%. (&%3 (3%&
&%&%# 磷肥推荐指标体系的建立 在山东省 ((个
“1$#$”完全实施试验和 #(个部分实施试验中，试验
点土壤有效磷含量集中于中等级范围（#.! 1.
45 6 75）的试验点 ’(个，肥料效应方程计算的平均磷
肥推荐用量（E&D’）为 -/ J 3 75 6 <4&；土壤有效磷含
量为低等级（ K #. 45 6 75）的试验点有 ##个，肥料效
应方程计算的平均磷肥推荐用量（E&D’）为 ##’ J #$
75 6 <4&；土壤有效磷含量处于高等级（1.! ’.
45 6 75）的试验点 #1个，肥料效应方程计算的平均磷
肥推荐用量（E&D’）为 3# J #1 75 6 <4&。土壤有效磷
在 ’. 45 6 75以上的有 &个，肥料效应方程计算的平
均磷肥推荐用量（E&D’）为 #. J #. 75 6 <4&。以磷肥
的当季效应为主，再综合考虑养分收支平衡对土壤
磷肥力的长期影响和磷肥恒量监控的思路［#3］，给出
不同土壤有效磷等级下的磷肥推荐量（表 &）。
&%&%& 钾肥推荐指标体系的建立 山东省 (( 个
“1$#$”完全实施试验和 #3个部分实施试验中，土壤
速效钾含量为中等级（’.!#.. 45 6 75）的试验点有
&-个，肥料效应方程计算的平均钾肥用量（F&D）为
##. J #1 75 6 <4&。1$个试验点的土壤速效钾含量处
于高等级（#..!#$. 45 6 75），肥料效应方程计算的
平均钾肥用量（F&D）为 (- J / 75 6 <4&。&.个试验点
土壤速效钾高于 #$. 45 6 75，肥料效应方程计算的平
均钾肥用量（F&D）为 ’$ J #. 75 6 <4&；但考虑到实际
生产中小麦对钾素不敏感，并参考以往相关研究结
果［#-］，建议在土壤速效钾含量高于 #$. 45 6 75的极
..& 植 物 营 养 与 肥 料 学 报 #’卷
表 ! 山东省土壤有效磷丰缺指标及冬小麦磷肥推荐用量表
"#$%& ! ’()% #*#)%#$%& + ,-#.&/ #0. + 1&-2)%)3&-
-&4(55&0.#2)(0 1(- 6)02&- 67
土壤有效磷等级
!"#$% &’
(&)* +*(%,-.
相对产量
/%*#0)1% 2)%*$
（3）
土壤有效磷
4&)* +*(%,-.
（56 7 86）
磷肥推荐用量
/%9&55%,$ . "#0%
（.:+; 86 7 <5:）
极高 =%"2 <)6< > ?; > ;@ @!
高 A)6< ?@!?; B@!;@ C@!D@
中 E%$)F5 G;!?@ H@!B@ D@!H@@
低 I&J K G; K H@ H@@!HB@
!或根据当地实际情况施种肥（.:+;）H;!B@ 86 7 <# &"（.:+;）L#(#*
’%"0)*)M%" #0 N*#,0),6 H;OB@ 86 7 <# #NN*)%$ #99&"$),6 0& *&9#* ()0F#0)&,P
高等级土壤上，不施用钾肥。此外，土壤速效钾含量
为低等级（ K ;@ 56 7 86）的试验点仅有 H 个，肥料效
应方程计算的钾肥推荐用量（Q:+）为 ::; 86 7 <5:，由
于点数太少，难以代表该水平土壤上的钾肥推荐用
量。根据“磷钾肥恒量监控技术”［HG］进行钾肥用量
推荐，当土壤速效钾含量为低钾时，建议钾肥用量为
作物收获带走钾素的 :倍，在保证施钾增产的前提
下，土壤速效钾含量能稳步提高，在当地秸秆还田的
生产条件下，推荐低钾土壤的钾肥用量（Q:+）为 H:@
!H;@ 86 7 <5:。不同土壤速效钾等级下的钾肥推荐
量见表 B。
表 8 山东省土壤速效钾丰缺指标及冬小麦钾肥推荐用量表
"#$%& 8 ’()% 9:;<=4>? ,-#.&/ #0. ? 1&-2)%)3&-
-&4(55&0.#2)(0 1(- 6)02&- 67
土壤速效钾等级
!"#$% &’
(&)* RAS+T9-Q
相对产量
/%*#0)1% 2)%*$
（3）
土壤速效钾
4&)* RAS+T9-Q
（56 7 86）
钾肥推荐用量
/%9&55%,$
Q "#0%
（Q:+ 86 7 <5:）
极高 =%"2 <)6< > ?; > HS@ @
高 A)6< ?@!?; H@@!HS@ C@!D@
中 E%$)F5 G;!?@ ;@!H@@ H@@!H:@
低 I&J K G; K ;@ H:@!H;@
8 讨论
8@A 土壤肥力分级指标
建立合理的土壤肥力分级指标是有效利用土壤
养分测试值进行推荐施肥的基础［B］，缺素区相对产
量是制定土壤养分肥力分级指标的重要参数［H?］。
国内外科研工作者曾用不同的相对产量值进行土壤
分级指标的划分：U&N% 等［HD］用相对产量值 ;@3、
G;3和 H@@3为标准进行土壤肥力等级的划分；在
国内，黄德明等［S］以缺素小麦相对产量 G;3、D;3、
?;3和 H@@3为标准，将土壤有效养分划分为 ;个级
别；林守宗等［;］用缺磷区小麦相对产量 ;@3、G@3、
?@3和 ?;3对土壤有效磷肥力分级进行划分。
V#6%")#等［:@］以相对产量值 G@3、?;3和 H@@3为标
准将水稻田的土壤有效磷划分为极低、中低、中和高
S个级别；陈新平建议我国农田土壤大量元素养分
的划分级别可以采用 > ?;3为高，G;3!?;3为中，
;@3!G;3为低，K ;@3为极低［H;OHC］。
自上世纪 D@年代以来，华北平原土壤有效磷由
于大量施用磷肥而普遍提高。对当前山东省小麦试
验点的土壤有效养分含量进行分析，大部分试验点
土壤有效磷含量集中在 H@!B@ 56 7 86之间，占 GD3；
超过 B@ 56 7 86的试验点 H; 个，仅占 HD3；低于 H@
56 7 86 的试验点只有 HH个。同时，在小麦钾肥的相
关试验点中，土壤速效钾含量低于 ;@ 56 7 86的试验
点只有 H个，不足 :3；在 ;@!H@@ 56 7 86的试验点
有 :D个，占 BS3；大于 H@@ 56 7 86的试验点 ;S个，占
CC3。因此，考虑到指标划分的准确性和实用性，应
弱化低和极低肥力水平 :个级别的划分，同时，考虑
到过高的土壤磷肥力可能带来的环境问题，细分高
肥力水平土壤为高和极高 :个水平，最终将山东省
土壤磷钾肥力指标划分为 S级：相对产量 > ?;3为
极高、?@3!?;3为高、G;3!?@3为中、K G;3为
低。
随着作物品种、产量水平和土壤肥力的变化，不
同年代的土壤肥力指标也会发生变化［HB］。在美国
W&J#州，目前把土壤速效钾含量低于 ?@ 56 7 86的土
壤分级为极低水平，但该养分含量在以往的分级指
标中是中肥力水平［:H］。在第二次土壤普查中，林守
宗等［;，D］将山东棕壤区土壤有效磷划分为 S级，K B
56 7 86为极低，B!D 56 7 86为低，D!:H 56 7 86为中，:H
!:G 56 7 86为高；将潮土区土壤有效磷划分为 ;
级，K BX; 56 7 86为极低，BX;!D 56 7 86为低，D!H;
56 7 86为中，H;!::X; 56 7 86为高，> ::X; 56 7 86为极
高。本试验建立的新指标与之相比，不仅极低级别
已经取消，其他等级范围内土壤有效磷含量也大幅
度地提高。造成这一现状的主要原因可能是，与第
二次土壤普查相比，当前作物产量和土壤养分含量
大幅度提高，导致作物磷素需求量增加［HH］。
相比而言，山东省土壤速效钾的指标变化不大。
例如，在第二次土壤普查中，林守宗等［;，D］将山东棕
壤区的速效钾分级指标分为 S 级，K :G 56 7 86为极
H@:H期 孙义祥，等：应用“BSHS”试验建立冬小麦测土配方施肥指标体系
低，!"!#" $% & ’%为低，#"!() $% & ’%为中，()!**+
$% & ’%为高；将潮土区速效钾在 *,+!*"+ $% & ’%的
土壤定义为丰钾区，这与本试验中山东当前产量水
平下土壤速效钾含量达到 - )+ $% & ’%为低，)+!*++
为中，*++!*#+ $% & ’%为高，高于 *#+ $% & ’% 为极高
的标准差异不大。
!"# 模型拟合与推荐施肥指标
“,#*#”试验方案采用回归最优设计，处理少、效
率高，试验可建立三元二次肥料效应回归方程，一次
模拟可以给出氮、磷和钾肥的最佳施用量，如果三元
二次回归方程模拟成功率太低，也就失去了该试验
设计的优势。在本研究条件，三元二次模型拟合的
成功率只有 ./（表 *），比王圣瑞等［*#］报道的“,#*#”
试验的模型拟合成功率还要低。其中一个非常重要
的原因是，相当多的试验田块由于土壤有效养分较
高，氮、磷、钾三要素中一个或多个元素施肥不增产
或增产不显著，因此使总体的肥料效应不能按三元
二次肥料效应回归方程拟合。
在此情况下，“,#*#”试验方案的优点是可以对
单因素进行一元肥料效应分析。本研究发现，通过
单因素分析，磷、钾的一元二次模型拟合成功率可以
达到 ,+/以上（表 *）。而在一元肥料效应函数拟合
不成功的试验中，有相当一部分是由于施肥不增产，
也有一部分可能与试验设计或管理有关，如所有试
验没有设置重复等。考虑到部分田块土壤磷、钾肥
力较高的实际情况，对出现非典型式或拟合最佳产
量施肥量不理想的试验数据可根据作物产量与施肥
量的关系，结合专业判断确定最佳施肥量，如对于施
肥不增产或增产幅度过小的田块，可以考虑将最佳
施肥量定为 +或仅施少量“启动肥”。
总之，本试验结果表明，与 !+ 世纪 (+ 年代相
比，由于土壤肥力和作物产量水平发生了变化，土壤
有效养分的丰缺指标、特别是土壤有效磷的丰缺指
标发生了很大变化，而应用“,#*#”试验建立大田作
物磷、钾的测土配方施肥指标体系是可行的；同时，
在最佳施肥量的确定中，需要注意施肥模型的选择
和专业知识的运用。
参 考 文 献：
［*］ 01%2341 5 6，7184%13 9 :; <2=>?@=2 ?A B?$$?@ C21@，D>81@E 34B2，
B?3@，FG21H，1@E =?IC21@ H? =?48 A23H484HI ?A 1@ ?J4=?8［K］; K ; L81@H 5DM
H3 ;，*.."，!+：*!".N*!(.O
［!］ P3?2G 0 <，PG?$>=?@ Q R; S?48 1@E =?48 A23H484HI（S4JHG 2E4H4?@）［R］;
T?F1，U$34B1：781B’F288 LDC84=G4@%，!++)O !)*N!V(O
［,］ 01%2341 5 6，7184%13 9 :，K?@2= : U; W41%@?=H4B H2BG@4XD2= A?3 @DH34M
H4?@18 E4=?3E23=; Y3?FHG 1@E $4@2318 @DH34H4?@ ?A A428E B3?>=［R］; 52F
Z?3’，U$234B1：R13B28 W2’’23，T@B;，*.."O (,N*,#
［#］ 黄德明 ;我国农田土壤养分肥力状况及丰缺指标［K］;华北农学
报，*.((，,（!）：#VN),O
[D1@% W R; S?48 A23H484HI 1@E A23H484HI 4@E2J A?3 B3?>81@E =?48= 4@ :G4@1
［K］; UBH1 U%34B; 7?32184MS4@;，*.((，,（!）：#VN),O
［)］ 林守宗，阎华，赵树慧 ;山东棕壤土养分丰缺指标研究初报［K］;
山东农业科学，*.()，（)）：*)N*.O
Q4@ S \，Z1@ [，\G1? S [; L3284$4@13I 32>?3H ?@ =?48 A23H484HI 4@E2J
A?3 C3?F@ =?48= 4@ SG1@E?@% >3?]4@B2［ K］; SG1@E?@% U%34B; SB4 ;，
*.()，（)）：*)N*.O
［V］ 周广业，丁宁平，黄世伟 ;用肥料田间试验与测土指标作为指导
合理施肥的探讨［K］;土壤通报，*.()，*V（*）：,(N#*O
\G?D Y Z，W4@% 5 L，[D1@% S ^; W4=BD==4?@ ?@ 321=?@1C82 A23H484_23
32B?$$2@E1H4?@ C1=2E ?@ HG2 =?48 H2=H4@% 4@E2J 1@E A428E 2J>234$2@H=
［K］; :G4@; K; S?48 SB4 ;，*.()，*V（*）：,(N#*O
［"］ 黄德明，徐建铭，武书敏，等 ;北京郊区小麦测土施肥技术研究
—土壤养分丰缺指标法在测土施肥中的应用［K］;华北农学通
报，*.(V，*（*）：#*N#"O
[D1@% W R，‘D K R，^ D S R !" #$ % SHDE42= ?A =?48 H2=H= 1@E A23H484_23
32B?$$2@E1H4?@= A?3 FG21H B3?> 4@ 724a4@% 1321［K］; UBH1 U%34B; 7?32M
184MS4@;，*.(V，*（*）：#*N#"O
［(］ 林守宗，赵树慧，阎华 ;潮土土壤养分丰缺指标与施肥的研究
［K］;山东农业科学，*.("，（)）：!)N!"，#+O
Q4@ S \，\G1? S [，Z1@ [; SHDE42E ?@ HG2 =?48 A23H484HI 4@E2J 1@E A23M
H484_23 32B?$$2@E1H4?@= ?@ 08D]?M1XD4B =?48 1321=［K］; SG1@E?@% U%34B;
SB4 ;，*.("，（)）：!)N!"，#+O
［.］ 陆允甫，吕晓男 ;中国测土施肥的工作进展和展望［K］; 土壤学
报，*..)，,!（,）：!#*N!)!O
QD Z 0，Qb ‘ 5; L3?%32== 1@E >3?a2BH 4@ A23H484_23 32B?$$2@E1H4?@
C1=2E ?@ =?48 H2=H4@% 4@ BG4@1［K］; UBH1 L2E?8 ; S4@;，*..)，,!（,）：
!#*N!)!O
［*+］ 黄德明 ;十年来我国测土施肥的进展［K］; 植物营养与肥料学
报，!++,，.（#）：#.)N#..O
[D1@% W R; S?48 H2=H4@% 1@E A23H484_23 32B?$$2@E1H4?@= 4@ :G4@1 ED3M
4@% HG2 >1=H E2B1E2［K］; L81@H 5DH3 ; 023H ; SB4 ;，!++,，.（#）：#.)N
#..O
［**］ 甄兰，崔振岭，陈新平，等 ; !)年来种植业结构调整驱动的县
域养分平衡状况的变化—以山东惠民县为例［K］;植物营养与
肥料学报，!++"，*,（!）：!*,N!!!O
\G2@ Q，:D4 \ Q，:G2@ ‘ L !" #$ % PG2 BG1@%2= ?A @DH342@H C181@B2=
E34]2@ CI 1EaD=H4@% B3?> =H3DBHD32 1H B?D@HI 82]28 ED34@% !) I213=—U
B1=2 =HDEI 4@ [D4$4@ :?D@HI，SG1@E?@% L3?]4@B2［K］; L81@H 5DH3 ;
023H ; SB4 ;，!++"，*,（!）：!*,N!!!O
［*!］ 张世熔，黄元仿，李保国，等 ;黄淮海冲积平原区土壤速效磷、
钾的时空变异特征［K］;植物营养与肥料学报，!++,，.（*）：,N
(O
\G1@% S <，[D1@% Z 0，Q4 7 Y !" #$ % P2$>?318M=>1H418 ]1341C484HI ?A
=?48 1]1481C82 >G?=>G?3D= 1@E >?H1==4D$ 4@ HG2 188D]418 32%4?@ ?A HG2
[D1@%M[D14M[14 L814@［K］; L81@H 5DH3 ; 023H ; SB4 ;，!++,，.（*）：,
N(O
［*,］ :?J 0 <; <1@%2 4@ =?48 >G?=>G?3D= B34H4B18 82]28= F4HG H4$2［K］; S?48
!+! 植 物 营 养 与 肥 料 学 报 *)卷
!"#$ %&$ ’$，())*，+,（+）：(+-./(+-)0
［(.］ 王圣瑞，陈新平，高祥照，等 $“1.(.”肥料试验模型拟合的探
讨［’］$植物营养与肥料学报，*--*，2（.）：.-)/.(10
3456 ! 7，89:5 ; <，=4" ; > !" #$ % !?@AB "5 CD&@E4?D"5 "F“1.(.”
F:G?DEDH:G :IJ:GD&:5?C［’］$ /.(10
［(+］ 陈新平，张福锁 $通过“1.(.”试验建立测土配方施肥技术指标
体系［’］$中国农技推广，*--,，**（.）：1,/1)0
89:5 ; <，>9456 L !$ MC?4NEDC9D56 F:G?DEDH4?D"5 G:#"&&:5A4?D"5
D5A:I N4C:A "5“1.(.”FD:EA :IJ:GD&:5?C［’］$ 89D54 %6GD#$ O:#9$
MI?:5C$，*--,，**（.）：1,/1)0
［(,］ 张福锁 $测土配方施肥技术要览［P］$北京：中国农业大学出
版社，*--,0 )1/((-0
>9456 L !$ P4D5 J"D5?C "5 C"DE ?:C?D56 45A F:G?DEDH:G G:#"&&:5A4?D"5
?:#95"E"6B［P］$ Q:DRD56：89D54 %6GD#@E?@G4E S5DT:GCD?B )1/((-0
［(U］ 陈新平，张福锁，崔振岭，等 $小麦—玉米轮作体系养分资源综
合管理理论与实践［P］$北京：中国农业大学出版社，*--,0 U+
/2.0
89:5 ; <，>9456 L !，8@D > V !" #$ % O9:"GB 45A JG4#?D#: "F 5@?GD:5?
G:C"@G#:C &4546:&:5? D5 W9:4?X#"G5 #"GJ"G4?D"5 CBC?:&［P］$ Q:DRD56：
89D54 %6GD#@E?@G4E S5DT:GCD?B ［(2］ 李秋梅 $高肥力土壤上冬小麦 Y夏玉米轮作中磷钾肥合理施用
的研究［Z］$北京：中国农业大学硕士学位论文，*--(0
VD [ P$ !?@AB "5 J9"CJ9"G@C 45A J"?4CCD@& 4JJED#4?D"5 F"G WD5?:G
W9:4? 45A C@&&:G &4DH: G"?4?D"5 D5 9D69 F:G?DED?B C"DE［Z］$ Q:DRD56：
PC ?9:CDC "F 89D54 %6GD#@E?@G4E S5DT:GCD?B，*--(0
［()］ 8"J: ’ O，MT45C 8 M$ !"DE ?:C?D56［’］$ %AT45$ !"DE !#D $，()2+（(）：
*-(/**20
［*-］ L46:GD4 K \，!45?"C % Q，Q4ED64G ] 8$ <9"CJ9"G@C C"DE ?:C? #4EDNG4X
?D"5 F"G E"WE45A GD#: "5 45 D5#:J?DC"E［’］$ %6G"5$ ’$，())U，2)（+）：
U1U/U.*0
［*(］ P4EE4GD5" % <，3D??GB Z ’，Q4GN46:E4?4 < %$ K:W C"DE ?:C? D5?:GJG:X
?4?D"5 #E4CC:C F"G J"?4CCD@&［’］$ Q:? $ 8G"JC （.）：(*/(.0
1-*(期 孙义祥，等：应用“1.(.”试验建立冬小麦测土配方施肥指标体系