水稻水肥耦合效應(yīng)研究

張 超，林彥宇，聶堂哲

(東北農(nóng)業(yè)大學(xué)農(nóng)業(yè)部農(nóng)業(yè)水資源高效利用重點(diǎn)實(shí)驗(yàn)室，哈爾濱 150030)

水和肥是水稻生長(zhǎng)過(guò)程中最重要的兩大影響因子[1,2]，但在生產(chǎn)中如果不合理利用這兩大資源，不僅會(huì)造成資源的浪費(fèi)，而且也會(huì)對(duì)環(huán)境造成嚴(yán)重的威脅，因此，水肥耦合效應(yīng)的研究引起了國(guó)內(nèi)外許多學(xué)者的關(guān)注與重視[3-5]。朱慶森等人對(duì)不同土壤水分狀況下氮素對(duì)水稻產(chǎn)量的影響及其機(jī)理進(jìn)行了研究，認(rèn)為在輕度干旱下(土壤水勢(shì)在-30 kPa)，適當(dāng)?shù)脑黾邮┑坑欣谒井a(chǎn)量的提高[6,7]。但以往的研究都集中在水肥單因子或水氮耦合方面，而對(duì)水肥多因子的互作效應(yīng)并未做深入系統(tǒng)的研究。

目前，人工神經(jīng)網(wǎng)絡(luò)模型有很多種，BP(Back Propagation)[8]神經(jīng)網(wǎng)絡(luò)是人工神經(jīng)網(wǎng)絡(luò)的重要模型之一，應(yīng)用極為廣泛，其具有很好的分布存儲(chǔ)性和容錯(cuò)性，非常適合解決非線性問(wèn)題，而且在實(shí)測(cè)資料不夠完備的情況下仍可以進(jìn)行計(jì)算和預(yù)測(cè)；遺傳算法(Genetic Algorithm,GA)[9]是一種新型的優(yōu)化算法，能有效的利用歷史信息來(lái)推斷下一代期望性能有所提高的期望點(diǎn)集。這樣一代代的不斷進(jìn)化最后收斂到一個(gè)最適應(yīng)環(huán)境的個(gè)體上，求得問(wèn)題的最優(yōu)解，而基于實(shí)數(shù)編碼的遺傳算法(Real Coded Accelerating Genetic Algorithm,RAGA)[10]又克服了二進(jìn)制編碼的缺點(diǎn)，使算法的尋優(yōu)能力大大加強(qiáng)；這兩種方法都很適合求得非線性問(wèn)題，但又有各自的缺點(diǎn)，BP網(wǎng)絡(luò)的誤差函數(shù)為平方型，存在局部極小值問(wèn)題且收斂速度較慢；遺傳算法在建立系統(tǒng)層次結(jié)構(gòu)等方面有著獨(dú)自的特點(diǎn)，但在學(xué)習(xí)、訓(xùn)練能力上遠(yuǎn)不如BP模型，因此將這兩種模型有機(jī)結(jié)合起來(lái)，取長(zhǎng)補(bǔ)短，可更準(zhǔn)確、高效的求出最優(yōu)值。

基于RAGA的BP模型在水文地質(zhì)、食品醫(yī)藥等領(lǐng)域已廣泛應(yīng)用[11,12]，但在作物水肥耦合研究方面，國(guó)內(nèi)外研究還是很罕見(jiàn)。為了探討水稻產(chǎn)量與水肥多因子之間的復(fù)雜關(guān)系，找出最優(yōu)的水肥配施方案，本文選取施氮量、施鉀量、施磷量、分蘗末期土壤含水率占飽和含水率的百分比4個(gè)設(shè)計(jì)因子，以產(chǎn)量為目標(biāo)因子，在2011年和2012年進(jìn)行了連續(xù)兩年的盆栽試驗(yàn)，在2011年試驗(yàn)中，經(jīng)RAGA-BP神經(jīng)網(wǎng)絡(luò)模型計(jì)算出的最佳水肥施入量，于2012年在相同條件下進(jìn)行了試驗(yàn)驗(yàn)證，其預(yù)測(cè)結(jié)果在合理的試驗(yàn)誤差范圍內(nèi)，表明此模型對(duì)生產(chǎn)實(shí)踐具有指導(dǎo)意義，并對(duì)制定水肥配施方案提供相應(yīng)的參考。

1 材料與方法

1.1 試驗(yàn)區(qū)概況

本試驗(yàn)在東北農(nóng)業(yè)大學(xué)園藝試驗(yàn)站溫室內(nèi)進(jìn)行，供試水稻品種為東農(nóng)427。供試土壤為黑土，其基本理化性質(zhì)如下：有機(jī)質(zhì)28.56 g/kg、全氮1.45 g/kg、全磷0.78 g/kg、全鉀20.76 g/kg、速效氮115.82 mg/kg、速效磷54.3 mg/kg、速效鉀182.6 mg/kg、pH值為6.52，土壤體積飽和含水率為51.47%。

1.2 試驗(yàn)材料

供試肥料為尿素(含N46%)、 鉀肥(含K2O 40% )、 磷酸二銨(含N18%，含P2O546% )。試驗(yàn)盆缽為可再生性密閉圓桶，上口直徑30 cm，圓桶高30 cm，每桶裝土15 kg，定苗4株。

1.3 試驗(yàn)方法

1.3.1二次飽和D-416最優(yōu)設(shè)計(jì)

本研究采用二次飽和D-416最優(yōu)設(shè)計(jì)[13]制定試驗(yàn)方案，選取施氮量、施鉀量、施磷量、分蘗末期土壤含水率占飽和含水率的百分比4個(gè)因素為設(shè)計(jì)因子，編制因子水平編碼表，見(jiàn)表1。試驗(yàn)共16個(gè)處理(見(jiàn)表3)，每個(gè)處理重復(fù)5次，共計(jì)80次。

表1 因素水平編碼表Tab.1 Factor level codes

施氮量為0～1.06 g/盆(折合公頃量為0～150 kg/hm2)，比例為基肥：分蘗肥∶拔節(jié)肥∶穗肥=3∶3∶3∶1，基肥在插秧前施入；分蘗肥在水稻移栽后7～10 d施入；拔節(jié)肥在水稻的幼穗長(zhǎng)1～2 mm施入；穗肥在抽穗后5～10 d施入。

施鉀量為0～0.71 g/盆(折合公頃量為0～100 kg/hm2)，施入比例為基肥∶穗肥=1∶1。

施磷量為0～0.53 g/盆(折合公頃量為0～75 kg/hm2)，作為基肥一次性施入。

計(jì)算施入土壤中的化肥量，其計(jì)算公式為施肥量=推薦施肥量/化肥的有效含量。

以分蘗末期土壤含水率占飽和含水率的下限值作為基本設(shè)計(jì)參數(shù)，其他各生育時(shí)期與其的比例關(guān)系為分蘗(前∶中∶后) ∶拔節(jié)孕穗期(前∶后)∶抽穗開(kāi)花期∶乳熟期=(1. 3∶1.15∶1)∶(1.15∶1.3)∶1.3∶1.15(見(jiàn)表2)[14]。

表2 各生育時(shí)期水分管理表Tab.2 Water management at different growth stages

注：分蘗后期曬田5～7 d；表2中單位為%的數(shù)字表示田間無(wú)水層情況下，土壤水分占飽和含水率的百分比。

1.3.2產(chǎn)量的測(cè)定

在成熟期應(yīng)考種并測(cè)定其有效穗數(shù)、穗粒數(shù)、千粒重、結(jié)實(shí)率，最終計(jì)算理論產(chǎn)量。

(1)有效穗數(shù)。應(yīng)查每一處理的結(jié)實(shí)的穗數(shù)，取其平均值。

(2)穗粒數(shù)。以每盆水稻計(jì)數(shù)，自然落干后統(tǒng)計(jì)其每穗粒數(shù)，然后計(jì)算出每一處理平均穗粒數(shù)。

(3)千粒重。將其晾曬干后使籽粒充分混合，任意取出3組，每組300粒，分別稱重，當(dāng)各組的質(zhì)量相差不到3%時(shí)，平均重即為千粒重，如差值超過(guò)3%，再取300粒稱重，用最為接近的3組數(shù)值平均值作為千粒重。

(4)結(jié)實(shí)率。將其晾曬干后使籽粒充分混合，任意取出3組，每組300粒，數(shù)其空癟數(shù)，當(dāng)各組的質(zhì)量相差不到3%時(shí)，實(shí)粒數(shù)與總粒數(shù)的比值即為結(jié)實(shí)率，如差值超過(guò)3%，再取300粒數(shù)其空癟數(shù)，用比值為最接近的3組數(shù)值平均值作為結(jié)實(shí)率。

(5)理論產(chǎn)量。產(chǎn)量=有效穗數(shù)×穗粒數(shù)×千粒重×結(jié)實(shí)率。

1.3.3RAGA-BP神經(jīng)網(wǎng)絡(luò)模型

RAGA-BP神經(jīng)網(wǎng)絡(luò)模型[15]尋優(yōu)主要分為BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練擬合和遺傳算法極值尋優(yōu)兩步，算法流程如圖1所示。

(1)輸入層節(jié)點(diǎn)的確定。根據(jù)影響產(chǎn)量的主要因素施氮量、施鉀量、施磷量、分蘗末期土壤含水率占飽和含水率的百分比等因素建立BP模型，這4個(gè)因素作為模型的輸入節(jié)點(diǎn)。

(2)隱含層節(jié)點(diǎn)數(shù)的確定。隱含層節(jié)點(diǎn)數(shù)的確定由于沒(méi)有明確的方法，其理論計(jì)算較為復(fù)雜，一般有如下公式來(lái)確定隱含層神經(jīng)元個(gè)數(shù)的范圍，即：

(1)

式中：l為隱含層節(jié)點(diǎn)數(shù)；n為輸入節(jié)點(diǎn)數(shù)；m為輸出節(jié)點(diǎn)數(shù)；a為1～10之間的調(diào)節(jié)常數(shù)。

本文通過(guò)計(jì)算隱含層節(jié)點(diǎn)數(shù)在3～12之間，然后通過(guò)不同神經(jīng)元的網(wǎng)絡(luò)進(jìn)行訓(xùn)練對(duì)比，發(fā)現(xiàn)當(dāng)網(wǎng)絡(luò)的隱含層節(jié)點(diǎn)數(shù)為10時(shí)，網(wǎng)絡(luò)具有足夠的泛化能力和輸出精度，且網(wǎng)絡(luò)的訓(xùn)練步數(shù)較少，因此確定網(wǎng)絡(luò)的隱含層節(jié)點(diǎn)數(shù)為10。

(3)輸出層節(jié)點(diǎn)的確定。輸出層節(jié)點(diǎn)的設(shè)置根據(jù)需要預(yù)測(cè)的性能參數(shù)來(lái)確定，在此網(wǎng)絡(luò)中以產(chǎn)量作為輸出節(jié)點(diǎn)。

(4)傳遞函數(shù)及算法。隱含層傳遞函數(shù)為T(mén)ansig和函數(shù)Logsig，輸出層傳遞函數(shù)為Purelin。誤差算法采用Levenberg-Marquardt，該算法較傳統(tǒng)BP算法而言，其梯度下降要快得多，從而在整個(gè)網(wǎng)絡(luò)的收斂上能以很少的迭代次數(shù)達(dá)到誤差要求。

(5)遺傳算法極值尋優(yōu)。將訓(xùn)練后的BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)結(jié)果作為個(gè)體適應(yīng)度值，通過(guò)選擇、交叉和變異操作尋找全局最優(yōu)值時(shí)對(duì)應(yīng)的輸入值。

圖1 算法流程圖Fig.1 Flow diagram of algorithm

2 結(jié)果與分析

2.1 飽和設(shè)計(jì)結(jié)果與分析

依據(jù)表3試驗(yàn)結(jié)果可以看出，當(dāng)施氮量為0.84 g/盆，施鉀量為0.56 g/盆，施磷量為0.42 g/盆，分蘗末期土壤含水率占飽和含水率的75%時(shí)，水稻的產(chǎn)量最高，達(dá)72.34 g/盆。但由于二次飽和設(shè)計(jì)僅考慮這幾個(gè)試驗(yàn)點(diǎn)間的部分試驗(yàn)結(jié)果，并未從整個(gè)條件區(qū)間進(jìn)行考慮，存在偶然現(xiàn)象或容易忽略最優(yōu)條件，所以本文嘗試?yán)肦AGA-BP神經(jīng)網(wǎng)絡(luò)來(lái)找出它們之間的變化關(guān)系，并最終求出最優(yōu)值。

表3 二次飽和D-416最優(yōu)設(shè)計(jì)表Tab.3 Quadratic saturation D-416 design Table

2.2 BP神經(jīng)網(wǎng)絡(luò)模型

用函數(shù)輸入輸出數(shù)據(jù)訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)，使訓(xùn)練后的網(wǎng)絡(luò)能夠擬合非線性函數(shù)輸出，保存訓(xùn)練好的網(wǎng)絡(luò)用于計(jì)算個(gè)體適應(yīng)度值。利用飽和試驗(yàn)設(shè)計(jì)中的16組數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的輸入信號(hào)，對(duì)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，選擇最大訓(xùn)練次數(shù)500次，訓(xùn)練時(shí)間0.1 s，期望誤差10-5，再通過(guò)訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)對(duì)隨機(jī)產(chǎn)生的16組數(shù)據(jù)進(jìn)行預(yù)測(cè)，檢查網(wǎng)絡(luò)的外推性能，試驗(yàn)測(cè)得值與網(wǎng)絡(luò)預(yù)測(cè)值結(jié)果如圖2所示，預(yù)測(cè)誤差如圖3所示。從圖2和圖3可以看出，BP神經(jīng)網(wǎng)絡(luò)可以準(zhǔn)確的預(yù)測(cè)非線性函數(shù)的輸出，可以將網(wǎng)絡(luò)預(yù)測(cè)輸出近似看成函數(shù)的實(shí)際輸出。

圖2 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)Fig.2 BP neural network forecast

圖3 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)誤差Fig.3 Error of BP neural network forecast

2.3 遺傳算法主函數(shù)

遺傳算法采用實(shí)數(shù)編碼的加速遺傳算法(RAGA)，由于尋優(yōu)函數(shù)只有4個(gè)輸入?yún)?shù)，所以個(gè)體長(zhǎng)度為4，個(gè)體適應(yīng)度值為BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值，適應(yīng)度值越小，個(gè)體越優(yōu)。本文RAGA中種群規(guī)模選為20，交叉概率0.4，變異概率0.2，進(jìn)化次數(shù)100次。優(yōu)化過(guò)程中最優(yōu)個(gè)體適應(yīng)度值變化曲線如圖4所示。RAGA得到的最優(yōu)個(gè)體適應(yīng)度值為73.55 g/盆，對(duì)應(yīng)的最優(yōu)個(gè)體為，即施氮量為1.01 g/盆、施鉀量為0.63 g/盆、施磷量為0.46 g/盆、分蘗末期土壤含水率占飽和含水率的75.2%；

圖4 適應(yīng)度變化曲線Fig.4 Changing curves of fitness

2.4 模型的驗(yàn)證

將RAGA-BP模型中得出的水肥參數(shù)在2012年度相同條件下進(jìn)行試驗(yàn)與計(jì)算，得到的水稻產(chǎn)量實(shí)測(cè)值與預(yù)測(cè)值見(jiàn)表4，與模型的相對(duì)誤差僅1.68%，說(shuō)明基于RAGA-BP神經(jīng)網(wǎng)絡(luò)模型用于利用水肥施入量去預(yù)測(cè)水稻產(chǎn)量方面是可行的。

3 結(jié) 語(yǔ)

(1)通過(guò)二次飽和D-416最優(yōu)設(shè)計(jì)結(jié)合RAGA-BP神經(jīng)網(wǎng)絡(luò)，得出水稻產(chǎn)量最高時(shí)最優(yōu)水肥方案，即RAGA得到的最優(yōu)個(gè)體適應(yīng)度值為73.55 g/盆，對(duì)應(yīng)的最優(yōu)個(gè)體為，即施氮量為1.01 g/盆、施鉀量為0.63 g/盆、施磷量為0.46 g/盆、分蘗末期土壤含水率占飽和含水率的75.2%；此參數(shù)在2012年度試驗(yàn)中對(duì)應(yīng)實(shí)際產(chǎn)量達(dá)74.78 g/盆，與網(wǎng)絡(luò)預(yù)測(cè)值相差為1.68%。

(2)為提高寒地黑土區(qū)水稻的產(chǎn)量，快速、準(zhǔn)確的選擇最優(yōu)水肥組合，本文在飽和設(shè)計(jì)的基礎(chǔ)上，提出了基于RAGA-BP神經(jīng)網(wǎng)絡(luò)模型對(duì)產(chǎn)量進(jìn)行預(yù)測(cè)的方法，在2011年試驗(yàn)的基礎(chǔ)上于2012年重新進(jìn)行了試驗(yàn)，并對(duì)試驗(yàn)結(jié)果在同等條件下進(jìn)行了驗(yàn)證，誤差僅為1.68%。結(jié)果表明，運(yùn)用飽和設(shè)計(jì)與RAGA-BP神經(jīng)網(wǎng)絡(luò)模型能較好的反映出產(chǎn)量與水肥施入量之間的復(fù)雜非線性關(guān)系，對(duì)指導(dǎo)農(nóng)業(yè)生產(chǎn)試驗(yàn)具有一定的參考意義。

[1] Gan Y T,Lafond G P,May W E. Grain yield and water use relative performance of winter v s. spring cereals in east-central Saskatchewan[J].Canadian Journal of Plant Science,2000,80：533-541.

[2] 王小彬，高緒科，蔡典雄.旱地農(nóng)田水肥相互作用的研究[J].干旱地區(qū)農(nóng)業(yè)研究，1993，11(3)：6-11.

[3] Karasov V G,Irrigation efficiency in water delivery[J].Technology,1982,2：62-74.

[4] Terry A Howell.Enhancing water use efficiency in irrigated agriculture[J].Agronomy Journal,2000,2：690-697.

[5] 文宏達(dá)，劉玉柱，李曉麗，等.水肥耦合與旱地農(nóng)業(yè)持續(xù)發(fā)展[J].土壤與環(huán)境，2002，11(3)：315-318.

[6] 朱慶森，邱澤森，姜長(zhǎng)鑒. 水稻各生育期不同土壤水勢(shì)對(duì)產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)科學(xué)，1994，27(6)：15-22.

[7] 楊建昌，朱慶森，王志琴，等.不同土壤水分狀況下氮素營(yíng)養(yǎng)對(duì)水稻產(chǎn)量的影響及其生理機(jī)制的研究[J].中國(guó)農(nóng)業(yè)科學(xué)，1995，29(4)：58-66.

[8] 高 雋.人工神經(jīng)網(wǎng)絡(luò)原理及仿真實(shí)例[M].北京：機(jī)械工業(yè)出版社，2003.

[9] Holland J H.Genetic algorithms[M].Science American,1992.

[10] 金菊良，楊曉華，丁 晶.標(biāo)準(zhǔn)遺傳算法的改進(jìn)方案----加速遺傳算法[J].系統(tǒng)工程理論與實(shí)踐，2001，(4)：8-13.

[11] 葉 咸，許 模，廖曉超，等.遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)在求解水文地質(zhì)參數(shù)中的應(yīng)用[J].水電能源科學(xué)，2013，31(12)：55-58.

[12] 基于BP人工神經(jīng)網(wǎng)絡(luò)和改進(jìn)遺傳算法的鉤藤堿的提取工藝優(yōu)化研究[J].計(jì)算機(jī)與現(xiàn)代化，2012，(8)：17-20.

[13] 徐仲儒.農(nóng)業(yè)試驗(yàn)最優(yōu)回歸設(shè)計(jì)[M].哈爾濱：黑龍江科學(xué)技術(shù)出版社，1988.

[14] 魏永霞，何雙紅.控制灌溉條件下水肥耦合對(duì)水稻產(chǎn)量及其構(gòu)成因子的影響[J].灌溉排水學(xué)報(bào)，2010，29(5)：98-102.

[15] Matlab中文論壇.Matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析[M].北京：北京航空航天大學(xué)出版社，2010.