发明内容
本发明的目的是针对上述现有技术的不足,提出一种基于氮素光谱指数法的水稻追氮调控方法,该方法可以快速、准确地估测出水稻植株氮积累量,并基于水稻目标需氮量、土壤供氮量来精确指导水稻氮肥投入。
本发明是通过以下技术方案实现的:
一种基于氮素光谱指数法的水稻追氮调控方法,通过采集样品光谱数据无损估算植株实时氮积累量,并对中后期氮肥利用率定量模拟,然后基于养分平衡原理,在综合考虑水稻目标产量需氮量、生育中后期的土壤供氮量、产量潜力以及中后期氮肥利用率和光谱估算的植株实时氮积累量基础上,按照如下公式(1)得到水稻中后期氮肥施用量:
水稻中后期氮肥施用量=[(目标产量需氮量-植株实时氮积累量)-中后期土壤供氮量]/中后期氮肥利用率 (1)。
其中,所述的水稻目标产量需氮量通过如下方法获得:先根据水稻产量目标模型确定目标产量,然后根据百公斤籽粒需氮量算法求得百公斤籽粒需氮量,进而求出水稻目标产量需氮量;计算公式如下:
目标产量需氮量=目标产量×百公斤籽粒需氮量/100 (2)。
所述的目标产量可以根据决策点前三年生态点高产田平均产量和增产系数来确定:
目标产量=前三年生态点高产田平均产量×(1+增产系数) (3)。
所述的百公斤籽粒需氮量算法计算步骤为先确定高产条件下的百公斤籽粒需氮量,然后根据目标产量影响因子和品种类型影响因子,求得百公斤籽粒需氮量,公式如下:
百公斤籽粒需氮量=高产百公斤籽粒需氮量×min(目标产量影响因子,1)-品种类型影响因子 (4),
目标产量影响因子=α×目标产量/生态点最高产量+β (5),
其中高产百公斤籽粒需氮量,通过当地高产条件下的历史数据获得;品种类型影响因子,对水稻来说籼稻取0.2,粳稻取0,杂交稻取-0.1;α、β分别为0.4773、0.50。
所述植株实时氮积累量可以根据拔节期水稻冠层获得差值植被指数无损估算求得,计算公式如式(6)所示,拔节期水稻冠层差值植被指数DVI可以通过公式(7)求得:
植株氮积累量=129.98×DVI(760,710)1.5293 (6),
DVI(760,710)=R760-R710 (7),
其中,R760和R710分别为拔节期水稻冠层Cropscan 760nm和710nm的光谱反射率,光谱测试采用光谱仪测定,测试时探头垂直于冠层1m,每小区测试5点取平均值作为该小区光谱反射率值。
所述的中后期土壤供氮量计算步骤为先确定零氮处理整个水稻生育期的土壤供氮量,然后根据中后期土壤供氮占土壤总供氮量的比例,进而求得水稻中后期土壤供氮量,中后期土壤供氮量计算公式为:
中后期土壤供氮量=零氮处理整个生育期土壤总供氮量×中后期土壤供氮量所占比例 (8)。
所述的零氮处理整个生育期的土壤供氮量和中后期土壤供氮占土壤总供氮量的比例可以根据零氮处理产量和土壤质地因子进行估算,计算公式分别如式(9)和式(10)所示:
零氮处理整个生育期土壤总供氮量=a×零氮处理产量-b (9),
中后期土壤供氮量所占比例=c×零氮处理产量-d (10),
其中,a、b、c、d为土壤质地因子,由文献资料整理确定,中间类型土壤取两者平均值。
所述的中后期氮肥利用率可以根据效应因子及其相应的相对权重来模拟求得,计算公式如式(11)所示:
所述的效应因子主要包括氮肥类型因子、氮肥施用方法因子、施肥次数因子、基蘖肥因子以及水稻品种因子。
与现有技术比较本发明的有益效果:
本发明是基于养分平衡原理,根据水稻目标产量需氮量、中后期氮肥利用率,基于反射光谱估算的植株实时氮积累量以及施氮时土壤供氮量,来快速、定量推荐水稻氮肥追施用量。本发明在实时信息获取方面,利用植株冠层群体信息,具有更快速、无损等特点,避免了“以点代面”现象的出现;在追氮调控模型方面,对中后期的土壤供氮量和中后期氮肥利用率进行了量化,较好地体现了不同土壤类型和施氮方式等因素对追氮量的影响,从而增强了追氮量估算的机理性和准确度,提高了量化水平。与现有技术比较本发明提出了一种新型的基于氮素光谱指数法的水稻追氮调控方法,并提供了中后期土壤供氮量和穗肥表观利用率的估算模型,克服了传统方法存在的耗时、费力等缺点,有助于提高水稻生产的精确化和数字化管理水平。
具体实施方式
实施例1
如图1所示流程图,于2007年、2008年在南京市农林局江宁试验站开展了不同基蘖肥的田间氮肥试验,以测试基于氮素光谱指数法的水稻追氮调控方法。供试品种为武香粳14,试验设6组不同基蘖氮肥用量的处理(N1、N2、N3、N4、N5、N6),从而人为创造出不同水稻长势状况,在此基础上使用本发明所涉及的水稻追氮调控方法进行水稻中期氮肥追施用量的定量估算和推荐,从而获得不同的氮肥施用量(N1r,N2r,N3r,N5r,N6r)。其中N4为正常高产田施氮量,作为对照,不进行调控施肥,同时设置氮空白区(全生育期不施用氮肥,N0)作为对照,按照本发明方法无损估算植株实时氮积累量,详细计算步骤如下:
水稻中后期氮肥施用量=[(目标产量需氮量-植株实时氮积累量)-中后期土壤供氮量]/中后期氮肥利用率 (1)
目标产量需氮量=目标产量×百公斤籽粒需氮量/100 (2)
目标产量=前三年生态点高产田平均产量×(1+增产系数) (3)
百公斤籽粒需氮量=高产百公斤籽粒需氮量×min(目标产量影响因子,1)-0.2 (4)
目标产量影响因子=0.4773×目标产量/生态点最高产量+0.5 (5)
根据基蘖肥用量取值不同,前三年生态点高产田平均产量和增产系数也相应变化,其中07年低基蘖肥条件下前三年生态点高产田平均产量和增产系数分别为8550kg·hm-2、15%,中基蘖肥条件下分别为9400kg·hm-2、10%,高基蘖肥条件下分别为10100kg·hm-2、5%;08年低基蘖肥条件下前三年生态点高产田平均产量和增产系数分别为8550kg·hm-2、15%,高基蘖肥条件下分别为9950kg·hm-2、5%;高产百公斤籽粒需氮量均取值2.2,生态点最高产量均取值10500kg·hm-2。
根据目标产量和百公斤籽粒需氮量可以得出N1、N2、N3、N5、N6目标产量需氮量分别为185、185、200、208、203kg·hm-2。
植株氮积累量=129.98×DVI(760,710)1.5293 (6)
DVI(760,710)=R760-R710 (7)
式中,R760和R710分别为拔节期水稻冠层760nm和710nm的光谱反射率,光谱反射率测定采用Cropscan公司生产的MSR-16型多光谱仪,测试时探头垂直于冠层1m,每小区测试5点取平均值作为该小区光谱反射率,其中N1、N2、N3、N5、N6760nm和710nm的光谱反射率分别为0.45(0.011)、0.406(0.01485)、0.479(0.0056)、0.504(0.004)、0.59125(0.0036),进而根据实时光谱反射率可求得N1、N2、N3、N5、N6植株氮积累量分别为37、30.93、41.4、45、57.65kg·hm-2。
中后期土壤供氮量=零氮处理整个生育期土壤总供氮量×中后期土壤供氮量所占比例 (8)
零氮处理整个生育期土壤总供氮量=a×零氮处理产量-b (9)
中后期土壤供氮量所占比例=c×零氮处理产量-d (10)
其中a、b、c、d为土壤质地因子,分别为0.026、-59.224、0.002、24.48,零氮处理产量为8150kg·hm-2,进而求得中后期土壤供氮量为66kg·hm-2。
效应因子分别为氮肥类型、氮肥施用方法、施肥次数、基蘖肥、水稻品种因子,效应因子相对权重(RW)可以通过式(12)求得;氮肥类型因子(FNT),控释肥、硫铵、尿素,碳铵氮肥利用率分别取值为0.95、0.9、0.85和0.8;氮肥施用方法因子(FAM),无水层撒施、水层撒施、土壤深施氮肥利用率分别取值0.85、0.9、0.95;施肥次数因子(FAT),通过式(13)求得,x为追肥次数;基蘖肥因子(FBQ),通过式(14)求得,其中FB为基蘖肥用量,ND为水稻目标需氮量;水稻品种因子(FVT),杂交稻、籼稻、粳稻分别赋值0.95、0.9、0.85。
追肥次数为2次,N1~N6基蘖肥用量分别为60、65、120、130、180、195kg·hm-2,进而求得N1、N2、N3、N5、N6中后期氮肥利用率分别为40%、44%、65%、81%、83%。
由公式(1),根据公式(2)~(14)计算结果可以求得N1、N2、N3、N5、N6中后期氮肥需用量(纯氮)分别为205、200、143、114、95kg·hm-2。
利用本发明估算的植株实时氮积累量进行追氮,将基于本发明方法指导的水稻产量和效益与对照进行比较。本发明方法指导的水稻产量与对照相比,平均增产4.9t·hm-2(图3),其中高氮低调节氮66kg·hm-2(图2),氮肥农学利用率提高3.3kg·kg-1(图4),增效3119元·hm-2(图5)(以07年为例)。
实施例2
如图6所示评价流程图,于2009年在南京市农林局江宁试验站设置2种不同基蘖肥的田间氮肥试验(低基蘖肥处理L、正常基蘖肥处理N),人为创造出不同水稻长势,在水稻穗肥追施前,利用包括本发明在内的5种氮肥调控方法(氮素光谱指数法NSI、叶面积指数法LAI【1】、实地氮肥管理SSNM【2】、氮素营养指数法NNI【3】、氮肥优化算法NFOA【4】)计算水稻实时需氮量,指导氮肥投入,从而比较评价本发明与已有技术的优缺点。供试品种为武香粳14。在此基础上设置氮空白区(全生育期不施用氮肥,N0)和常规施氮处理(SN,施氮量为270kg·hm-2)作为对照,最后将5种方法指导的水稻产量和效益进行比较。
结果如图7~图9,与常规施氮处理相比,按照5种方法指导追氮的水稻产量没有明显差异,其中低基蘖肥条件下,氮素光谱指数法的施氮量最低,产量、氮肥回收率和经济效益相对较高,实地氮肥管理施氮量略低于对照,而另外3种方法施氮量较高,其氮肥回收率及经济效益较低,其中总施氮量高低相差148kg·hm-2,氮肥农学利用率高低相差8kg·kg-1,氮肥回收率相差22%,净利润高低相差532元·hm-2,产投比相差1.3;正常基蘖肥条件下,氮肥优化算法施氮量与常规施氮基本相同,产量和经济效益略低于常规施氮处理,其它4种方法均减少了氮肥投入,提高了氮肥利用率和经济效益,其中总施氮量高低相差63kg·hm-2,氮肥农学利用率相差5kg·kg-1,氮肥回收率相差13%,净利润高低相差229元·hm-2,产投比相差0.6;综合比较表明,氮素光谱指数法无论在低基蘖肥,还是在正常基蘖肥条件下均有较好的效果。
5种调控方法计算过程如下——
氮素光谱指数法NSI:
水稻中后期氮肥施用量=[(目标产量需氮量-植株实时氮积累量)-中后期土壤供氮量]/中后期氮肥利用率 (1)
目标产量需氮量=目标产量×百公斤籽粒需氮量/100 (2)
百公斤籽粒需氮量=高产百公斤籽粒需氮量×min(目标产量影响因子,1)-0.2 (4)
目标产量影响因子=0.4773×目标产量/生态点最高产量+0.5 (5)
根据基蘖肥用量取值不同,前三年生态点高产田平均产量和增产系数也相应变化,低基蘖肥条件下前三年生态点高产田平均产量和增产系数分别为8850kg·hm-2、15%,正常基蘖肥条件下分别为9770kg·hm-2、10%;高产百公斤籽粒需氮量均取值2.3,生态点最高产量均取值11250kg·hm-2。
根据目标产量和百公斤籽粒需氮量可以得出低基蘖肥和正常基蘖肥条件下目标产量需氮量分别为205、223kg·hm-2。
植株氮积累量=129.98×DVI(760,710)1.5293 (6)
DVI(760,710)=R760-R710 (7)
式中,R760和R710分别为拔节期水稻冠层760nm和710nm的光谱反射率,光谱反射率测定采用Cropscan公司生产的MSR-16型多光谱仪,测试时探头垂直于冠层1m,每小区测试5点取平均值作为该小区光谱反射率,其中低基蘖肥和正常基蘖肥条件下760nm和710nm的光谱反射率分别为0.731(0.002)、0.86(0.0015),进而根据实时光谱反射率可求得低基蘖肥和正常基蘖肥条件下植株氮积累量分别为80.15、107.2kg·hm-2。
中后期土壤供氮量=零氮处理整个生育期土壤总供氮量×中后期土壤供氮量所占比例 (8)
零氮处理整个生育期土壤总供氮量=a×零氮处理产量-b (9)
中后期土壤供氮量所占比例=c×零氮处理产量-d (10)
其中a、b、c、d为土壤质地因子,分别为0.026、-59.224、0.002、24.48,零氮处理产量为8150kg·hm-2,进而求得中后期土壤供氮量为66kg·hm-2。
效应因子分别为氮肥类型、氮肥施用方法、施肥次数、基蘖肥、水稻品种因子,效应因子相对权重(RW)可以通过式(12)求得;氮肥类型因子(FNT),控释肥、硫铵、尿素,碳铵氮肥利用率分别取值为0.95、0.9、0.85和0.8;氮肥施用方法因子(FAM),无水层撒施、水层撒施、土壤深施氮肥利用率分别取值0.85、0.9、0.95;施肥次数因子(FAT),通过式(13)求得,x为追肥次数;基蘖肥因子(FBQ),通过式(14)求得,其中FB为基蘖肥用量,ND为水稻目标需氮量;水稻品种因子(FVT),杂交稻、籼稻、粳稻分别赋值0.95、0.9、0.85。
追肥次数为2次,低基蘖肥和正常基蘖肥条件下基蘖肥用量分别为65、135kg·hm-2,进而求得低基蘖肥和正常基蘖肥条件下中后期氮肥利用率分别为45%、74%。
由公式(1),根据公式(2)~(14)计算结果可以求得低基蘖肥、正常基蘖肥条件下中后期氮肥需用量(纯氮)分别为144.6、72kg·hm-2。
叶面积指数法LAI:利用水稻资料修订Wood等提出的基于小麦追氮算法
(1)确定目标绿叶叶面积指数(GAI),在齐穗前绿叶叶面积指数基本等同植株叶面积指数。(e.g.7.5)
(2)测定当前的GAI,GAI可以用比值植被指数RVI进行估算(e.g.1.64)。
LAI=0.08×RVI(1100,560)2-0.07×RVI(1100,560)+0.48 (15)
(3)计算当前GAI与目标GAI的差值(7.5-1.64=5.86)。
(4)设定作物每增长1个单位GAI需氮为30kg·hm-2。
(5)根据步骤(3)、(4)确定达到预期GAI的氮需求量(5.86×30=175.8kg·hm-2)。
(6)确定土壤供氮(66kg·hm-2),扣除后则为需要施入的氮(175.8-66=109.8kg·hm-2)。
(7)确定氮肥利用率(低基蘖肥条件为0.45/正常基蘖肥条件为0.74),进而求出需要氮肥施用量(低基蘖肥条件为247kg·hm-2/正常基蘖肥条件为109kg·hm-2)。
实地氮肥管理SSNM:
确定南京地区武香粳14各时期SPAD阈值及施氮量,拔节期和孕穗期阈值及施氮量相同,即SPAD阈值为45,当SPAD>45,追施40kg[N]hm-2;43<SPAD<45,施用70kg[N]hm-2;如果SPAD<43,则追施100kg[N]hm-2。
根据实时水稻叶片SPAD值,可以求得低基蘖肥条件下追氮量为200kg·hm-2,正常基蘖肥条件下追氮量为110kg·hm-2。
氮素营养指数法NNI:
NNI通过观测区差值植被指数DVI除以适宜区差值植被指数DVI值求得,当NNI<1时,表示氮肥亏缺;NNI>1时,表示氮肥充足。如缺氮,则通过公式求得:
NNI=DVIfert/DVIref (16)
ΔN=675.84×NNI-699.34 (17)
Nr=Ns+ΔN(Ns为135kg·hm-2) (18)
其中DVIfert为观测区DVI值,DVIref为适宜区DVI值,ΔN为缺少的氮肥量,Nr为需施氮肥量,Ns为适宜区中后期施氮量。
利用公式(16)~(18)可以求得低基蘖肥条件下追氮量为291kg·hm-2,正常基蘖肥条件下追氮量为102kg·hm-2。
氮肥优化算法NFOA:该方法根据作物产量潜力确定总施氮量,并根据实时植株氮积累量确定最终追氮量
将INSEY确定为DVI(760,710)与移栽后天数(DAT)的比值;水稻当季产量系数(INSEY)与产量,以及DVI(760,710)与植株氮积累量(PNA)的拟合关系见图6-3和6-4;目标产量(PGY)和目标产量需氮量(GNA)的算法如下,其中NUE取值分别为0.45、0.74,计算参照方法1。
(1)当季产量系数(INSEY)计算:INSEY=DVI(760,710)/DAT (19)
(2)根据产量系数预测方程确定作物目标产量(PGY)
PGY=13456.29×INSEY+3804.94 (20)
(3)计算氮肥总需用量(GNA):GNA=2.1×PGY/100 (21)
(4)计算当前植株氮积累量(PNA):PNA=0.0208×DVI(760,710)2.0553 (22)
(5)穗肥施用量(Nr)计算:Nr=(GNA-PNA)/NUE (23)
根据公式(19)~(23)可以求得低基蘖肥条件下追氮量为247kg·hm-2,正常基蘖肥条件下追氮量为135kg·hm-2。
参考文献
[1]Wood G A,Welsh J P,Godwin R J,et al.Real-time measures of canopy size as a basis forspatially varying nitrogen applications to winter wheat sown at different seed rates[J].BiosystemsEngineering,2003,84:513-531.
[2]Peng S,Garcia F V,Laza R C,et al.Increased N-use efficiency using a chlorophyll meter onhigh-yielding irrigated rice[J]Field Crops Research,1996,47:243-252.
[3]Xue L,Yang L. Recommendations for nitrogen fertiliser topdressing rates in rice usingcanopy reflectance spectra[J].Biosystems Engineering,2008,100(4):524-534.
[4]Lukina E V,Freeman K W,Wynn K J,et al.Nitrogen fertilization optimization algorithmbased on in-season estimates of yield and plant nitrogen uptake[J].Journal of Plant Nutrition,2001,24(6):885-898.