目的:研究北方半夏在人工栽培条件下产量与施肥的数学模型,以指导大田科学施肥。方法:试验采用3因素5水平2次通用旋转组合设计。结果与结论:只有氮(N),磷(P),钾(K)配施,方可实现本试验条件下栽培半夏的最大产量:29 9214 kg·hm-2;此时施肥量N,五氧化二磷(P2O5),氧化钾(K2O)分别为41379,22479,16401 kg·hm-2,施肥纯收益为475 80921 元/hm2。肥料的主效应P>N>K。
Objective: To study the mathematical model of optimum fertilizer application to Pinellia by field experiments. Method: Three factors and 5 levels quadratic rotational combination design was applied to study the effects of N, P and K. Result and Conclusion: N, P and K combination application can achieve the maximum yield of 29 9214 kg·hm-2, the amount of N, P2O5, K2O were 41379, 22479, 16401 kg·hm-2, respectively. The effects of the 3 factors were in following orders: P>N>K.
全 文 :半夏高产高效栽培最佳施肥数学模型研究
王 鹏1,裴建文2,孙万仓1,杨少平3,刘文臖4,孙新荣1
(1.甘肃农业大学 农学院,甘肃 兰州 730070;
2.天水师范学院 现代农业与陇东南动植物资源研究保护开发中心,甘肃 天水 741001;
3.天水市麦积区农牧局 农业技术推广中心,甘肃 天水 741020;
4.天水市农业学校,甘肃 天水 741400)
[摘要] 目的:研究北方半夏在人工栽培条件下产量与施肥的数学模型,以指导大田科学施肥。方法:试验采用3因素5
水平2次通用旋转组合设计。结果与结论:只有氮(N),磷(P),钾(K)配施,方可实现本试验条件下栽培半夏的最大产量:
299214kg·hm-2;此时施肥量N,五氧化二磷(P2O5),氧化钾(K2O)分别为41379,22479,16401kg·hm
-2,施肥纯收益为
47580921元/hm2。肥料的主效应P>N>K。
[关键词] 半夏;产量;施肥
[收稿日期] 20080408
[通讯作者] 裴建文,Tel:(0941)3623266,13909415316,Email:
ibh203@126.com
随着半夏野生资源的不断减少和保护力度的加
大,人工栽培半夏逐渐已成为重要的栽培药材之
一[1]。目前,半夏栽培仅停留在群众经验的基础
上,技术模糊,资料奇缺,生产效益差异悬殊。为了
改变目前的生产现状,在当地主栽区作者研究氮
(N)、磷(P)、钾(K)施肥配方[25],旨在找到施肥与
产量的相关数学模型,为半夏施肥提供科学依据。
1 材料与方法
1.1 试验地基本情况
2007年试验设在甘肃省天水市清水县城郊川
水地上,海拔1378m;土质中壤,肥力中等。土壤养
分全氮063g·kg-1,全P2O5129g·kg
-1,全K2O
1685g·kg-1;碱解氮 680mg·kg-1,速效 P2O5
227mg·kg-1,速效K2O1735mg·kg
-1。前茬以
小麦做匀地试验,试验区均匀一致。
1.2 材料
1.2.1 试验用化肥 氮肥用尿素(含N46%,兰州
石化公司产),磷肥用过磷酸钙(含 P2O512%,白银
绿源磷复合肥有限责任公司产),钾肥用氯化钾(含
K2O60%,青海利源化肥厂产)。
1.2.2 供试半夏品种 为采挖的清水县百家乡野
生半夏。
1.3 方法
1.3.1 试验设计 采用3因子2次通用旋转组合
设计。
1.3.2 播种、田间管理及收获 为便于田间管理和
降低因化肥水平渗透作用对试验的影响,每小区做
一畦,四周打埂,埂宽10cm,高15cm。小区净面积
1225m2(35m×35m),共20处理,分2个区组
随机排列(“0”点处理重复6次)。播种时先铲出畦
内表土,整平畦底,将全部化肥混合后一次性施入畦
底,搂入土中混匀后撒播种茎,最后均匀覆土 15
cm。整个生育期内全试验未施有机肥,也不追施其
他肥料。2006年11月13日播种,播量(块茎)4200
kg·hm-2,密度510万粒/hm2。出苗前4月1日浇
透水1次,并地表覆盖麦草保持土壤松软湿润。5
月3日开始出苗,5月25日植株叶片已封住地面。
全生育期视苗情、土壤和大气湿度随时浅灌,并做到
田间无杂草。待全田倒苗,叶片干枯后于2007年9
月15日一次性采挖。每小区四周除去30cm宽的
面积不计产(以进一步消除肥料不同水平在小区间
的横向渗透影响),剩余部分严格计产。生育期间
病虫采用统一方案防治,效果理想。
1.3.3 土壤化验方法 全氮用半微量开氏法、全磷
用氢氧化钠熔融钼锑抗比色法、全钾用氢氧化钠熔
融火焰光度法、碱解氮用碱解扩散法、速效磷用碳
酸氢钠熔融钼锑抗比色法、速效钾用乙酸铵提取
火焰光度法[6]。
·966·
第34卷第6期
2009年3月
Vol.34,Issue 6
March,2009
1.4 统计方法
结果数据用软件DPS进行处理。
2 结果与分析
2.1 回归方程
具体的试验方案及结果见表1。
由表1资料得方程:
表1 试验方案及结果
处理
编码值
X1(N) X2(P2O5) X3(K2O)
施肥量/kg·hm-2
N P2O5 K2O
鲜块茎产量/kg·hm-2
Y Y′
产值/
元/hm2
1 -1 -1 -1 12165 6075 6075 203602 204951 32576320
2 -1 -1 1 12165 6075 23925 227614 226680 36418240
3 -1 1 -1 12165 23925 6075 254627 257106 40740320
4 -1 1 1 12165 23925 23925 260130 260576 41620800
5 1 -1 -1 47835 6075 6075 255628 259961 40900480
6 1 -1 1 47835 6075 23925 283642 281690 45382720
7 1 1 -1 47835 23925 6075 281641 287353 45062560
8 1 1 1 47835 23925 23925 291646 290823 46663360
9 -1682 0 0 0 150 150 214607 213900 34337120
10 1682 0 0 600 150 150 288644 285601 46183040
11 0 -1682 0 300 0 150 242121 241738 38739360
12 0 1682 0 300 300 150 296648 293281 47463680
13 0 0 -1682 300 150 0 259630 252661 41540800
14 0 0 1682 300 150 300 270635 273854 43301600
15 0 0 0 300 150 150 286143 285500 45782880
16 0 0 0 300 150 150 289645 285500 46343200
17 0 0 0 300 150 150 287644 285500 46023040
18 0 0 0 300 150 150 282141 285500 45142560
19 0 0 0 300 150 150 283642 285500 45382720
20 0 0 0 300 150 150 283142 285500 45302720
注:Y为实际值;Y′为依方程计算出的理论值;产值由各处理产量、半夏干块茎48元/kg及鲜块茎出干率3333%折算而来。
Y′=28549964+2131417X1 +1532205
X2+629981X3 -1263619X1
2 -635898X22 -
786192X3
2 -619050X1X2 + 106300X1X3 -
456475X2X3 (1)
方程经失拟性检验,拟和良好[F1=321,F005
(5,5)=505,即F1<F005],并达极显著水平(F2=
8508>F001=494)。
各回归系数经显著性检验,除 t13外亦均达显著
或极显著水平(t0=16569,t1=1864,t2=
134,t3=551,t12=414,t13=071,t23=
306,t11=1135,t22=571,t33=706)。
说明该方程可以用作施肥对产量的预报。
2.2 施肥对鲜块茎产量的影响
2.2.1 N,P,K配施 原方程⑴剔除不显著项后为:
Y′=28549964+2131417X1+1532205X2+
629981X3 - 1 263619X1
2 - 635898X2
2 -
786192X3
2-619050X1X2-456475X2X3 (2)
对方程⑵求X1,X2,X3的一阶偏导数,并令其为
0得方程组:
2131417-2527238X1-61905X2=0
1 532205 -61905X1 -1 271796X2 -
456475X3=0
629981-456475X2-1572384X3=0
求解知,当编码值 X1=0638,X2=0838,X3=
0157时,方程实现鲜块茎的最大产量299214kg
·hm-2。此时施肥量为 N41379kg·hm-2,P2O5
22479kg·hm-2,K2O16401kg·hm
-2;施肥成本
293328元/hm2,投产比为 16321,施肥纯收益
47580912元/hm2(现市场价 N380元/kg,P2O5
396元/kg,K2O287元/kg;半夏干块茎48元/kg,
栽培半夏的出干率3333%)。
求解经济最佳施肥量知,当N,P,K配施时栽培半
夏的最大产量施肥量与经济最佳施肥量间相差无几。
2.2.2 N,P配施 将 X3=-1682(不施 K)代入
产量方程⑵得:
Y′(N,P)=25266101+2131417X1+2299996
X2-619050X1X2-1263619X1
2-635898X2
2(3)
并计算出N,P配施的交互效应值见表2。
由表2可看出,不施 K肥的情况下,施 N量在
-1682≤X1≤0时,产量随 P2O5施量的增大而增
·076·
第34卷第6期
2009年3月
Vol.34,Issue 6
March,2009
表2 N,P肥对半夏产量的交互效应 kg·hm-2
X1(N)编码值
X2(P2O5)
-1682 -1 0 1 1682
-1682 1068713 1412899 1810612 20811462192705
-1 1516220 1831612 2187107 2415421 2498187
0 1959847 2233021 2526610 2693020 2733566
1 2150751 2381706 2613391 2717895 2716222
1682 2135995 2338156 2527621 2589907 2559440
加;施N量在1≤X1≤1682时,产量先随P2O5施量
的增大而增加,当P2O5施量越过1水平时,则下降。
P2O5施量在-1682≤X2≤1时,产量随施N量的增
大而增加,但N施量越过1水平时,却下降;当P2O5
量取1682水平时,产量也先随 N施量的增大而增
加,一越过0水平时,就下降。
N肥(或 P肥)的增产效果因 P肥(或 N肥)
的施量水平不同而不同。在本试验施 P范围内,
当 P2O5取 -1682,-1,0,1,1682水平时,平均
每 1kgN增产分别为 1779,1542,1195,
848,611kg,可见随P2O5施量增加,平均每1kg
N肥的增产效应逐渐减小;当 N取 -1682,-1,
0,1,1682水平时,平均每 1kgP2O5分别增产
3747,3273,2579,1885,1412kg,随 N施量
增加,平均每 kgP2O5肥增产效应也逐渐减小。可
见 N,P配施时,若 P2O5(或 N)施量较小,则单位
N肥(或 P2O5)增产效应较大,反之,则单位 N肥
(或 P2O5)增产效应较小。在不施任何化肥(X1,
X2,X3均取 -1682)时,产量最低,仅 1068713
kg·hm-2。
当编码值 X1=0455,X2=1587时,方程⑶则
实现N,P两素配施的最大产量施肥量和经济最佳
施肥量。
2.2.3 N,K配施 将 X2=-1682(不施 P)代入
产量方程(2)得:
Y′(N,K) = 24 173761 + 3 172659X1 +
1397772X3-1263619X1
2-786192X3
2 (4)
并计算出N,K配施的交互效应值见表3。
表3 N,K肥对半夏产量的交互效应 kg·hm-2
X1(N)编码值
X3(K2O)
-1682 -1 0 1 1682
-1682 1068713 1307845 1526241 1587399 1538923
-1 1516220 1755352 1973748 2034906 1986430
0 1959847 2198980 2417376 2478534 2430058
1 2150751 2389884 2608280 2669438 2620962
1682 2135995 2375128 2593524 2654682 2606206
由表3可见,不施P肥时,施N(或K2O)量在本
试验设计范围内,产量先随 K2O(或 N)施量的增大
而增加,但K2O(或 N)越过1水平时,则下降,二者
趋势相同;不管K2O(或N)的施量水平如何变化,平
均每1kgN都增产1779kg,而平均每1kgK2O都
增产1567kg。可见N,K配施时,增产效应N>K。
当编码值 X1=1255,X3=0889时,方程⑷则
达N,K两素配施的最大产量施肥量和经济最佳施
肥量。
2.2.4 P,K配施 将 X1=-1682(不施 N)代入
产量方程(2)得:
Y′(P,K) = 21 389986 + 2 573447X2 +
629981X3 - 456475X2 X3 - 635898X2
2 -
786192X3
2 (5)
并计算出P,K配施的交互效应值见表4。
从表4看出,不施N肥,当P2O5量在-1682≤
X2≤-1时,产量先随K2O施量的增大而增加,K2O
施量越过1水平时则下降;而在0≤X2≤1682时,
产量也先随K2O施量的增大而增加,K2O施量一越
过0水平,就下降。但在本试验设计范围内,K2O施
量不管取何水平产量均随P2O5施量的增大而增加。
这说明如果只是P,K配施,P2O5施量如高于本试验
设计范围,产量还可增加。
·176·
第34卷第6期
2009年3月
Vol.34,Issue 6
March,2009
表4 P,K肥对半夏产量的交互效应 kg·hm-2
X2(P2O5)编码值
X3(K2O)
-1682 -1 0 1 1682
-1682 1068713 1307845 1526241 1587399 1538923
-1 1412899 1630799 1818064 1848091 1778383
0 1810612 1997381 2138999 2123378 2022538
1 2081146 2236784 2332754 2271485 2139514
1682 2192705 2327111 2391949 2299549 2146346
K肥(或P肥)的增产效果因 P肥(或 K肥)的
施量水平不同而不同:当P2O5取-1682,-1,0,1
时,平均每 1kgK2O可分别增产 1567,1218,
706,195kg,当P2O5取1682时,平均每1kgK2O
却减产155kg(比不施K肥时X3=-1682),可见
随P2O5施量的增加,平均每1kgK2O的增产效应在
逐渐减小直至减产;与此相似,当K2O取-1682,-
1,0,1,1682时,平均每 1kgP2O5增产量分别为
3747,3398,2886,2374,2025kg,可见随 K2O施
量的增加,每1kgP2O5的平均增产效应也在递减。
当编码值 X1=2098,X3=-0209时,方程⑸
则达 P,K两素配施的最大产量施肥量和经济最佳
施肥量。
2.2.5 N,P,K单施 对方程⑵中 3因子分别取
-1682(不施),使方程降维得单因子产量方程:
Y′(N) = 19 598474 + 3 172659X1 -
1263619X1
2 (6)
Y′(P)=18106124+3341238X2-635898X2
2
(7)
Y′(K)=15262414+1397772X3-786192X3
2
(8)
可解出当编码值 X1=1255,X2=2627(已超
出本试验设计上限),X3=0889时,方程(6),(7),
(8)则分别达 N,P,K单施的最大产量施肥量和经
济最佳施肥量,见图1。依方(6),(7),(8)再算出
N,P,K单施的效应值见表5。
图1 N,P,K单施效应曲线
表5 N,P,K单施对半夏产量的效应 kg·hm-2
施肥
编码值
-1682 -1 0 1 1682
N 10687131516220195984721507512135995
P2O5 10687131412899181061220811462192705
K2O 10687131307845152624115873991538923
由表5看出,本试验设计范围内,单施时平均每
kgN,P2O5,K2O分别可增产鲜块茎 1779,3747,
1567kg,增产效应P>N>K。
3 结论
综上分析,计算出不同配施方式的经济效益见
表6。
表6 栽培半夏施肥配合方式的效益比较
方案 施肥方式
最大产量施肥量/kg·hm-2
(N∶P2O5∶K2O)
最大产量
/kg·hm-2
施肥成本
/元/hm2
施肥纯收益
/元/hm2
施肥投产比
最大产量施肥量 N,P,K配施 41379∶22479∶16401 299214 293328 47580912 16321
N,P配施 38115∶29164 275758 260326 43860954 16948
N,K配施 52383∶22934 267865 264876 42593524 16181
P,K配施 33725∶13135 240243 171248 38267632 22447
单施N 52383 215899 199055 34344785 17354
单施P 38446 224951 152246 35839914 23641
单施K 22934 158837 65821 25348099 38611
经济最佳施肥量 N,P,K配施 41379∶22470∶16401 299214 293292 47581017 16323
N,P配施 38115∶29164 275758 260326 43860954 16948
N,K配施 52383∶22934 267865 264876 42593524 16181
P,K配施 33725∶13135 240243 171248 38267632 22447
单施N 52383 215899 199055 34344785 17354
单施P 38446 224951 152246 35839914 23641
单施K 22934 158837 65821 25348099 38611
·276·
第34卷第6期
2009年3月
Vol.34,Issue 6
March,2009
从表6可见,最大产量施肥量和经济最佳施肥
量2种方案,只有 N,P,K配施时二者才略有差异,
其余完全相同。不同方式的施肥成本 N,P,K>N,
K>N,P>N>P,K>P>K;最大产量及施肥纯收益
N,P,K>N,P>N,K>P,K>P>N>K;投产比K>
P>P,K>N>N,P>N,P,K>N,K。还可看出,单施
K虽施肥成本最小,投产比最大,但每66667m2施
肥纯收益却最小;N,P,K配施虽施肥成本最高,投
产比偏小,但每 66667m2施肥纯收益却最高。可
见要实现半夏的高产高效栽培,单位面积上创造最
大经济效益,以N,P,K配施最好。
[参考文献]
[1] 冉懋雄.名贵中药材绿色栽培技术———半夏,水半夏,附子
[M].北京:科学技术文献出版社,2002:74.
[2] 陈伦寿,李仁岗.农田施肥原理与实践[M].北京:农业出版
社,984:62.
[3] 裴建文,王成家,蒲建刚.高山冷凉区马铃薯高产高效栽培施
肥数学模型研究[J].西北农业学报,1998,7(6):12.
[4] 裴建文,吕 汰,柴晓芹,等.天水市蔬菜主产区耕层土壤养分
评价及施肥建议[J].甘肃农业科技,2002(7):38.
[5] 张 明,钟国跃,王德立,等.半夏种植密度及底肥施用水平的
研究[J].现代中药研究与实践,2003,17(4):23.
[6] 中国科学院南京土壤研究所.土壤理化分析[M].上海:上海
科学技术出版社,1978:62.
Studyonmathematicalmodelofoptimumfertilizerapplicationto
Pineliaternatecultivation
WANGPeng1,PEIJianwen2,SUNWanchang1,YANGShaoping3,LIUWenjun4,SUNXinrong1
(1.ColegeofAgriculture,GansuAgriculturalUniversity,Lanzhou730070,China;
2.ProtectionandResearchInstituteofModernAgricultureandAnimalandPlantResourceofEastandSouthofGansu,
TianshuiNormalUniversity,Tianshui741001,China;
3.AgriculturalTechnologyPopularizationCenter,AgriculturalandAnimalHusbandryBureauof
MaijiCounty,Tianshui741020,China;
4.TianshuiAgriculturalSchool,Tianshui741400,China)
[Abstract] Objective:TostudythemathematicalmodelofoptimumfertilizerapplicationtoPineliabyfieldexperiments.
Method:Threefactorsand5levelsquadraticrotationalcombinationdesignwasappliedtostudytheefectsofN,PandK.Resultand
Conclusion:N,PandKcombinationapplicationcanachievethemaximumyieldof299214kg·hm-2,theamountofN,P2O5,
K2Owere41379,22479,16401kg·hm
-2,respectively.Theefectsofthe3factorswereinfolowingorders:P>N>K.
[Keywords] Pineliaternate;yield;applicationfertilizer
[责任编辑 吕冬梅]
·376·
第34卷第6期
2009年3月
Vol.34,Issue 6
March,2009