徐壽軍 王 磊 德木其格 王金波 齊海祥 朱國(guó)立 ，5

（1內(nèi)蒙古民族大學(xué)農(nóng)學(xué)院，028043，內(nèi)蒙古通遼；2內(nèi)蒙古自治區(qū)高校蓖麻產(chǎn)業(yè)工程技術(shù)研究中心，028043，內(nèi)蒙古通遼；3內(nèi)蒙古自治區(qū)蓖麻育種重點(diǎn)實(shí)驗(yàn)室，028043，內(nèi)蒙古通遼；4內(nèi)蒙古自治區(qū)蓖麻產(chǎn)業(yè)協(xié)同創(chuàng)新培育中心，028043，內(nèi)蒙古通遼；5內(nèi)蒙古通遼市農(nóng)業(yè)科學(xué)研究院，028015，內(nèi)蒙古通遼）

虛擬作物研究是借助計(jì)算機(jī)模擬作物生長(zhǎng)發(fā)育和空間形態(tài)變化過(guò)程，為作物研究提供了一種新的方法，在農(nóng)業(yè)、林業(yè)、生態(tài)和遙感等領(lǐng)域得到廣闊的應(yīng)用，已成為農(nóng)業(yè)信息技術(shù)研究的熱點(diǎn)之一[1-4]。而作物器官形態(tài)模擬模型是虛擬作物研究的關(guān)鍵技術(shù)，它是基于三維空間實(shí)現(xiàn)作物可視化的基礎(chǔ)。目前，眾多學(xué)者以多種作物為對(duì)象，構(gòu)建了大量的作物地上器官形態(tài)建成預(yù)測(cè)模型[5-9]。其中，關(guān)于水稻[10]、小麥[11]和玉米[12]等禾本科作物器官形態(tài)建成的動(dòng)態(tài)模型的研究較為系統(tǒng)深入，不僅模擬了莖稈、葉片和果穗等重要器官的生長(zhǎng)動(dòng)態(tài)，而且構(gòu)建了葉鞘、莖鞘夾角、根系和葉片曲線(xiàn)乃至葉色的生長(zhǎng)變化動(dòng)態(tài)模型，為這些作物的可視化研究奠定了良好基礎(chǔ)。而棉花[13-14]、大豆[15]和蓖麻[16]等雙子葉作物因具有無(wú)限生長(zhǎng)習(xí)性，器官種類(lèi)多且結(jié)構(gòu)復(fù)雜，對(duì)其器官形態(tài)建成的模擬研究還處于初步探索階段。尤其是關(guān)于蓖麻主莖形態(tài)建成的預(yù)測(cè)模型研究鮮見(jiàn)報(bào)道。本研究選擇地上器官差異較大的蓖麻品種，以生理發(fā)育時(shí)間為模擬尺度，在系統(tǒng)觀(guān)測(cè)的基礎(chǔ)上，構(gòu)建了蓖麻地上部器官形態(tài)建成的預(yù)測(cè)模型，為開(kāi)展蓖麻的可視化研究提供理論基礎(chǔ)。

1 材料與方法

1.1 供試材料

參試品種為通蓖9號(hào)、通蓖10號(hào)、通蓖11號(hào)和通蓖13號(hào)，均由內(nèi)蒙古通遼市農(nóng)業(yè)科學(xué)研究院提供。其中通蓖11號(hào)和通蓖13號(hào)植株較矮，莖稈較細(xì)；通蓖9號(hào)和通蓖10號(hào)植株較高，莖稈較粗。

1.2 試驗(yàn)設(shè)計(jì)

1.2.1 播期試驗(yàn) 于2016年在內(nèi)蒙古通遼市農(nóng)業(yè)科學(xué)研究院進(jìn)行播期試驗(yàn)，設(shè)4月17日、5月7日和5月27日3個(gè)播期。分2次施純氮120kg/hm2，其中基肥84kg/hm2，追肥36kg/hm2，于現(xiàn)蕾前施用。磷肥和鉀肥作為基肥一次性施用，用量分別為30和60kg/hm2。行距65cm，株距70cm，小區(qū)面積90m2。隨機(jī)區(qū)組設(shè)計(jì)，3次重復(fù)。栽培管理與當(dāng)?shù)卮筇锷a(chǎn)相同。

1.2.2 氮肥試驗(yàn) 于2017年在通遼市科爾沁區(qū)農(nóng)牧業(yè)高新科技示范園區(qū)進(jìn)行氮肥試驗(yàn)，設(shè)0、60、120和180kg/hm2共4個(gè)氮肥水平（純氮），基追比7∶3，追肥在現(xiàn)蕾期前施用。磷肥和鉀肥作為基肥一次性施用，用量分別為30和60kg/hm2。5月3日播種，行距65cm，株距60cm，小區(qū)面積60m2。隨機(jī)區(qū)組設(shè)計(jì)，3次重復(fù)。栽培管理與當(dāng)?shù)卮筇锷a(chǎn)相同。

1.3 測(cè)定項(xiàng)目和方法

1.3.1 生育期 在蓖麻主莖第1片復(fù)葉開(kāi)始展開(kāi)時(shí)，每小區(qū)選擇株高、莖粗和葉片擴(kuò)展等長(zhǎng)勢(shì)接近的3株蓖麻進(jìn)行標(biāo)記。記錄標(biāo)記植株第3復(fù)葉、主莖各節(jié)間、主莖果穗以及蒴果生長(zhǎng)開(kāi)始和結(jié)束的日期。

1.3.2 農(nóng)藝性狀 在每個(gè)節(jié)間生長(zhǎng)開(kāi)始至生長(zhǎng)結(jié)束的時(shí)間段內(nèi)，每3d測(cè)定1次標(biāo)記植株主莖的節(jié)間長(zhǎng)度和直徑，使用直尺測(cè)量節(jié)間長(zhǎng)度，使用游標(biāo)卡尺測(cè)量各節(jié)直徑（節(jié)間最粗處）。蓖麻同一品種各個(gè)節(jié)位復(fù)葉的形態(tài)基本相似。選取主莖第3復(fù)葉作為模擬對(duì)象，在第3片復(fù)葉開(kāi)始出現(xiàn)至停止生長(zhǎng)的時(shí)間段內(nèi)，每3d測(cè)量1次葉長(zhǎng)（葉片中心的最長(zhǎng)主脈的長(zhǎng)度）和葉寬（與最長(zhǎng)主脈垂直的兩葉裂間的最大距離）。自主莖果穗出現(xiàn)開(kāi)始，每3d用直尺測(cè)量1次果穗長(zhǎng)。選取主莖果穗的某一蒴果標(biāo)記，每3d用游標(biāo)卡尺測(cè)量1次蒴果長(zhǎng)和蒴果直徑（蒴果最粗處）。

1.3.3 植株含氮量 各處理分別在出苗期、現(xiàn)蕾期、開(kāi)花期、結(jié)果期和成熟期取樣，105℃下殺青0.5h，80℃下烘48h至恒重，粉碎后采用凱氏定氮法[17]測(cè)定氮含量。

1.4 數(shù)據(jù)處理

用播期試驗(yàn)中5月7日和5月27日播種的數(shù)據(jù)及氮肥試驗(yàn)數(shù)據(jù)構(gòu)建模型，并調(diào)試參數(shù)，用播期試驗(yàn)中4月17日播種的數(shù)據(jù)檢驗(yàn)?zāi)Ｐ汀?/p>

1.5 模型構(gòu)建方法與步驟

以生理發(fā)育時(shí)間作為模擬尺度。生理發(fā)育時(shí)間（physiological development time，PDT）是指在最適溫光條件下，作物完成某一發(fā)育階段所需的最短時(shí)間。對(duì)于某一特定基因型品種來(lái)說(shuō)，完成某一特定的生育階段所需的生理發(fā)育時(shí)間基本恒定[18]。蓖麻各生育階段所對(duì)應(yīng)的生理發(fā)育時(shí)間見(jiàn)文獻(xiàn)[19]。

在分析節(jié)間長(zhǎng)度、節(jié)直徑、葉片長(zhǎng)度、葉片寬度、果穗長(zhǎng)度、蒴果長(zhǎng)度和蒴果直徑隨PDT動(dòng)態(tài)變化的基礎(chǔ)上，運(yùn)用Richards方程W=A/（1+be-kt）m構(gòu)建各器官形態(tài)建成預(yù)測(cè)模型，采用顧世梁等[20]提出的縮張算法對(duì)b、k和m等參數(shù)進(jìn)行確定。

2 結(jié)果與分析

2.1 主莖生長(zhǎng)預(yù)測(cè)模型

2.1.1 節(jié)間伸長(zhǎng)模型 根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，蓖麻第一節(jié)間開(kāi)始伸長(zhǎng)的生理發(fā)育時(shí)間為10.9d，主莖各節(jié)間從開(kāi)始伸長(zhǎng)到定長(zhǎng)需要22.1個(gè)生理日，蓖麻主莖各節(jié)間的生長(zhǎng)具有一定的重疊性，存在疊次生長(zhǎng)的現(xiàn)象，相鄰節(jié)間疊次伸長(zhǎng)間隔1.7個(gè)生理日，見(jiàn)式（1）。

式（1）為節(jié)間伸長(zhǎng)預(yù)測(cè)模型，PDT代表生理發(fā)育時(shí)間，CINL(PDT)n代表蓖麻主莖某一生理發(fā)育時(shí)間第n節(jié)節(jié)間長(zhǎng)度（m），CINLmaxn代表蓖麻主莖第n節(jié)節(jié)間長(zhǎng)度極值（m），為品種遺傳參數(shù)；bnl、knl和mnl為模型參數(shù)；NDF為氮肥影響因子，計(jì)算方法見(jiàn)文獻(xiàn)[19]。式（1）為3段函數(shù)：當(dāng)PDT小于某一節(jié)間開(kāi)始伸長(zhǎng)的時(shí)間時(shí)，其節(jié)間長(zhǎng)度為0；當(dāng)PDT大于某一節(jié)間停止伸長(zhǎng)的時(shí)間時(shí)，其節(jié)間長(zhǎng)度為節(jié)間長(zhǎng)度極值CINLmaxn；當(dāng)PDT處于某一節(jié)間開(kāi)始伸長(zhǎng)和停止伸長(zhǎng)時(shí)間范圍內(nèi)時(shí)，其節(jié)間長(zhǎng)度可用Richards方程計(jì)算，并受氮肥因子影響。

本試驗(yàn)觀(guān)測(cè)結(jié)果表明，隨著蓖麻主莖節(jié)間數(shù)由基部到頂部的遞增，CINLmaxn呈逐漸增加趨勢(shì)，第1節(jié)間長(zhǎng)度極值最小，頂部節(jié)間長(zhǎng)度極值最大，第1節(jié)間至第5節(jié)間節(jié)間長(zhǎng)度極值增加緩慢，第5節(jié)間至第11節(jié)間的節(jié)間長(zhǎng)度極值增加較快，第12節(jié)間之后，節(jié)間長(zhǎng)度極值增加又趨于緩慢（圖1）。因此，蓖麻主莖各節(jié)間長(zhǎng)度極值可用式（2）計(jì)算。

圖1 節(jié)間長(zhǎng)度極值變化Fig.1 Internode length variation

模型參數(shù)bnl與CINLmaxn呈指數(shù)函數(shù)關(guān)系，可通過(guò)式（3）計(jì)算。

模型參數(shù)knl隨CINLmaxn呈“S”曲線(xiàn)變化，可用Logistic方程表達(dá)，見(jiàn)式（4）。

mnl與CINLmaxn亦呈指數(shù)函數(shù)關(guān)系，可通過(guò)式（5）計(jì)算。

2.1.2 各節(jié)直徑預(yù)測(cè)模型 根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，蓖麻主莖各節(jié)直徑開(kāi)始變粗的時(shí)間間隔為1.7個(gè)生理日，從開(kāi)始變粗到粗度固定時(shí)間跨度為22.1個(gè)生理日。各節(jié)直徑變化動(dòng)態(tài)由式（6）計(jì)算：

式（6）中，CINT(PDT)n代表蓖麻主莖某一生理發(fā)育時(shí)間第n節(jié)直徑（m），CINTmaxn代表蓖麻主莖第n節(jié)直徑極值，是品種遺傳參數(shù)。bnt、knt和mnt為模型參數(shù)。

隨著蓖麻主莖節(jié)間數(shù)由基部到頂部遞增，CINTmaxn呈逐漸降低的趨勢(shì)，第1節(jié)直徑極值最大，最頂部節(jié)直徑極值最?。▓D2）。各節(jié)直徑極值用式（7）計(jì)算。

圖2 各節(jié)直徑極值變化Fig.2 Change in diameter maximum of each section

式（7）中，CINTmax為頂部節(jié)間極值，是品種遺傳參數(shù)。式（7）表明，某一節(jié)間極值可以根據(jù)頂部節(jié)間極值，運(yùn)用Richards方程來(lái)計(jì)算。

bnt、knt和mnt等模型參數(shù)分別是CINTmaxn的指數(shù)函數(shù)、線(xiàn)性函數(shù)和對(duì)數(shù)函數(shù)關(guān)系，見(jiàn)式（8）、（9）和（10）。

2.2 葉片擴(kuò)展的預(yù)測(cè)模型

2.2.1 葉片伸長(zhǎng)的動(dòng)態(tài)模型 葉片伸長(zhǎng)動(dòng)態(tài)見(jiàn)圖3，其增長(zhǎng)動(dòng)態(tài)可用式（11）表達(dá)。

圖3 葉長(zhǎng)極值的動(dòng)態(tài)Fig.3 Dynamic of the maximum of leaf length

式中，CLLEN(PDT)代表某一生理發(fā)育時(shí)間蓖麻主莖葉長(zhǎng)（m），CLLENmax代表蓖麻主莖第3復(fù)葉的長(zhǎng)度極值（m），定為品種遺傳參數(shù)。根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，第3復(fù)葉出現(xiàn)的生理發(fā)育時(shí)間是18.55個(gè)生理日，停止生長(zhǎng)的生理發(fā)育時(shí)間是37.56個(gè)生理日。式（11）表明，在18.55～37.56個(gè)生理日的時(shí)間范圍內(nèi)，蓖麻主莖第3復(fù)葉的長(zhǎng)度可用Richards方程來(lái)計(jì)算，并受氮肥效應(yīng)的影響。

2.2.2 葉片增寬的動(dòng)態(tài)模型 葉片增寬動(dòng)態(tài)見(jiàn)圖4，葉片寬度極值在不同PDT階段變化不同，當(dāng)PDT＜5或者PDT＞15時(shí)，其增加較為緩慢，當(dāng)5＜PDT＜15時(shí)，其增加的幅度較大。

圖4 葉寬極值的動(dòng)態(tài)Fig.4 Dynamic of the maximum of leaf width

這種變化趨勢(shì)可用式（12）表達(dá)。

其中，CLWID(PDT)為某一生理發(fā)育時(shí)間蓖麻主莖第3復(fù)葉的寬度（m），CLWIDmax表示蓖麻主莖第3復(fù)葉的寬度極值（m），可定為品種遺傳參數(shù)。式（12）表明，在18.55～37.56個(gè)生理日的時(shí)間范圍內(nèi)，蓖麻主莖第3復(fù)葉的寬度可用Richards方程來(lái)計(jì)算，并受氮肥效應(yīng)的影響。

2.3 果穗生長(zhǎng)的預(yù)測(cè)模型

蓖麻主莖果穗長(zhǎng)變化的動(dòng)態(tài)見(jiàn)圖5，主莖果穗長(zhǎng)在前期（PDT≤15）增加幅度不大，后期增加幅度較大，可用式（13）表達(dá)。

圖5 穗長(zhǎng)的生長(zhǎng)動(dòng)態(tài)Fig.5 Dynamic of ear length

式中，CBLEN(PDT)為某一生理發(fā)育時(shí)間蓖麻主莖果穗長(zhǎng)（m），CBLENmax為蓖麻主莖果穗長(zhǎng)極值（m），可定義為品種遺傳參數(shù)。根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，主莖果穗開(kāi)始生長(zhǎng)的生理發(fā)育時(shí)間是33.48個(gè)生理日，停止生長(zhǎng)的生理發(fā)育時(shí)間是57.41個(gè)生理日。式（13）表明，在33.48～57.41個(gè)生理日的時(shí)間范圍內(nèi)，蓖麻主莖第3復(fù)葉的寬度可用Richards方程來(lái)計(jì)算，并受氮肥效應(yīng)的影響。

2.4 蒴果生長(zhǎng)的預(yù)測(cè)模型

2.4.1 蒴果長(zhǎng)變化的動(dòng)態(tài)模型 蓖麻蒴果長(zhǎng)動(dòng)態(tài)變化（圖6）顯示，其增加的時(shí)期主要在生育前期（PDT＜15），生育后期幾乎沒(méi)有變化，可用方程（14）表示：

圖6 蒴果長(zhǎng)極值的動(dòng)態(tài)Fig.6 Dynamic of the maximum of capsule length

式中，CBWID(PDT)為某一生理發(fā)育時(shí)間蓖麻主莖蒴果長(zhǎng)（m），CBWIDmax為蒴果長(zhǎng)極值（m），是品種遺傳參數(shù)。根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，蒴果開(kāi)始增長(zhǎng)的生理發(fā)育時(shí)間是33.48個(gè)生理日，停止增長(zhǎng)的生理發(fā)育時(shí)間是57.41個(gè)生理日。式（14）表明，在33.48～57.41個(gè)生理日的時(shí)間范圍內(nèi)，蓖麻蒴果長(zhǎng)變化呈不對(duì)稱(chēng)“S”曲線(xiàn)，分漸增、激增和緩增3個(gè)階段，可用Richards方程來(lái)計(jì)算，并受氮肥效應(yīng)的影響。

2.4.2 蒴果直徑變化的動(dòng)態(tài)模型 蒴果直徑動(dòng)態(tài)變化見(jiàn)圖7，可用式（15）表達(dá)。

圖7 蒴果直徑極值的動(dòng)態(tài)Fig.7 Dynamic of the maximum of capsule diameter

式中，CBTHI(PDT)為某一生理發(fā)育時(shí)間蓖麻蒴果直徑（m），CBTHImax為蓖麻蒴果直徑極值（m），是品種遺傳參數(shù)。根據(jù)文獻(xiàn)[19]的計(jì)算方法，運(yùn)用本試驗(yàn)的觀(guān)測(cè)數(shù)據(jù)計(jì)算可知，蒴果開(kāi)始增粗的生理發(fā)育時(shí)間是33.48個(gè)生理日，停止增粗的生理發(fā)育時(shí)間是57.41個(gè)生理日。

2.5 品種遺傳參數(shù)的確定及模型的檢驗(yàn)

2.5.1 品種遺傳參數(shù)的確定 表1顯示，4個(gè)供試蓖麻品種遺傳參數(shù)中，頂部節(jié)間長(zhǎng)在18.68～26.52cm，第1節(jié)間直徑在2.96～3.38cm，葉長(zhǎng)在13.72～17.00cm，葉寬在12.32～14.04cm，果穗長(zhǎng)在44.00～56.84cm，蒴果長(zhǎng)在2.30～2.54cm，蒴果直徑在1.90～2.04cm。通蓖10號(hào)各項(xiàng)參數(shù)值最高，通蓖13號(hào)最低。

表1 4個(gè)蓖麻品種遺傳參數(shù)Table 1 Genetic parameters of four castor varieties cm

2.5.2 模型的檢驗(yàn) 選擇通遼市農(nóng)業(yè)科學(xué)研究院試驗(yàn)點(diǎn)現(xiàn)蕾期的觀(guān)測(cè)數(shù)據(jù)，采用根均方差RMSE方法[18]對(duì)蓖麻主莖器官形態(tài)建成的預(yù)測(cè)模型進(jìn)行檢驗(yàn)。

式（16）中，OBSi為觀(guān)測(cè)值，SIMi為模擬值，n為樣本容量。

選擇各品種蓖麻主莖第3、5、7和9節(jié)間，對(duì)其長(zhǎng)度和直徑預(yù)測(cè)模型進(jìn)行檢驗(yàn)，結(jié)果（表2）顯示，所檢驗(yàn)節(jié)間長(zhǎng)度預(yù)測(cè)絕對(duì)誤差在0.10～0.72cm，RMSE在0.14～0.60cm；各節(jié)直徑的模擬誤差在0.03～0.16cm，RMSE在0.06～0.11cm。

表2 節(jié)間長(zhǎng)和主莖直徑預(yù)測(cè)誤差（n=16）Table 2 Prediction errors of internode length and main stem diameter on castor (n=16) cm

選擇苗期和灌漿期的觀(guān)測(cè)數(shù)據(jù)對(duì)蓖麻主莖第3復(fù)葉的生長(zhǎng)動(dòng)態(tài)模型和果穗長(zhǎng)、蒴果長(zhǎng)和蒴果直徑生長(zhǎng)動(dòng)態(tài)模型進(jìn)行檢驗(yàn)，結(jié)果（表3）顯示，葉長(zhǎng)的預(yù)測(cè)誤差在0.19～0.73cm，RMSE為0.47cm；葉寬的預(yù)測(cè)誤差在0.30～0.60cm，RMSE為0.38cm；果穗長(zhǎng)的模擬誤差在0.89～1.85cm，RMSE為1.42cm；蒴果長(zhǎng)的模擬誤差在0.11～0.21cm，RMSE為0.16cm；蒴果直徑的模擬誤差在0.05～0.12cm，RMSE為0.07cm。

表3 蓖麻葉長(zhǎng)、葉寬、果穗長(zhǎng)度、蒴果長(zhǎng)和蒴果直徑預(yù)測(cè)誤差（n=20）Table 3 Prediction errors of leaf length, leaf width, ear length, capsule length and capsule diameter of castor (n=20) cm

運(yùn)用各品種蓖麻主莖莖稈長(zhǎng)度、粗度、葉長(zhǎng)、葉寬、果穗長(zhǎng)、蒴果長(zhǎng)和蒴果直徑的觀(guān)測(cè)值和模擬值數(shù)據(jù)，繪制兩者之間的1∶1關(guān)系圖（圖8），各觀(guān)測(cè)點(diǎn)均集中在1∶1直線(xiàn)附近，表明模擬值與觀(guān)測(cè)值間誤差較小，顯示模型具有一定的可靠性。

圖8 蓖麻主莖莖稈、葉片及果實(shí)觀(guān)測(cè)值與模擬值比較Fig.8 Comparison of simulated values with observed values on stem, leaf and fruit of castor

3 討論

3.1 關(guān)于模型的建模思路

目前，作物器官形態(tài)建成模擬模型的構(gòu)建有以下4種思路：一是以積溫（生長(zhǎng)度日，GDD）為驅(qū)動(dòng)變量，構(gòu)建積溫與器官生長(zhǎng)的動(dòng)態(tài)關(guān)系。如周娟等[21]構(gòu)建了棉花主莖葉長(zhǎng)和葉寬、主莖葉柄長(zhǎng)、主莖節(jié)間長(zhǎng)和粗、果枝葉長(zhǎng)和寬、果枝葉柄長(zhǎng)、果節(jié)長(zhǎng)和粗以及棉鈴高度和直徑等隨GDD變化的動(dòng)態(tài)模型；孫愛(ài)珍等[22]構(gòu)建了水稻葉片幾何形態(tài)參數(shù)與有效積溫之間的動(dòng)態(tài)關(guān)系模型。二是以輻熱積為驅(qū)動(dòng)變量，構(gòu)建其與器官生長(zhǎng)的動(dòng)態(tài)關(guān)系，這類(lèi)模型多見(jiàn)于蔬菜和鮮花等溫室作物。如常毅博等[23]建立了番茄株高、單株展開(kāi)葉數(shù)及單葉面積隨輻熱積的變化模型；倪紀(jì)恒等[24]構(gòu)建了基于輻熱積的黃瓜單果果長(zhǎng)、果徑、開(kāi)花節(jié)位的動(dòng)態(tài)變化模型；李剛等[25]構(gòu)建了溫室切花百合出葉數(shù)、株高、第1花蕾長(zhǎng)度和第1花蕾直徑隨幅熱積的動(dòng)態(tài)變化模型。三是以生物量為驅(qū)動(dòng)變量，構(gòu)建其與作物器官生長(zhǎng)的動(dòng)態(tài)關(guān)系，這是目前大田作物器官形態(tài)建成模型中應(yīng)用最多的一類(lèi)，通過(guò)系統(tǒng)觀(guān)測(cè)作物器官生物量、體積、密度和直徑等指標(biāo)數(shù)據(jù)與器官生長(zhǎng)的動(dòng)態(tài)趨勢(shì)來(lái)建立模型。如宋有洪等[26]構(gòu)建了基于器官生物量的植株形態(tài)的玉米虛擬模型；陳超等[27]在系統(tǒng)分析生物量增加與器官增長(zhǎng)關(guān)系的基礎(chǔ)上，構(gòu)建了棉花葉片和棉鈴生長(zhǎng)動(dòng)態(tài)模型；展志崗等[28]描述了冬小麥葉片、莖稈和穗等器官隨其生物量變化的動(dòng)態(tài)方程；張偉欣等[29]構(gòu)建了基于生物量的油菜越冬前植株葉片空間形態(tài)結(jié)構(gòu)模型。四是以生理發(fā)育時(shí)間為驅(qū)動(dòng)變量，構(gòu)建其與作物器官生長(zhǎng)的動(dòng)態(tài)關(guān)系。如徐壽軍等[30-31]以生理發(fā)育時(shí)間為模型尺度，構(gòu)建了大麥和大豆地上部主要器官形態(tài)建成預(yù)測(cè)模型；楊恒山等[32]構(gòu)建了基于生理發(fā)育時(shí)間的玉米器官形態(tài)建成的模擬模型。

3.2 關(guān)于模型的機(jī)理性

本研究預(yù)測(cè)了常規(guī)栽培條件下蓖麻地上部葉片、莖稈和蒴果等器官的生長(zhǎng)過(guò)程，檢驗(yàn)結(jié)果表明所建模型具有較好的預(yù)測(cè)性和可靠性。與既有模型比較，本研究在模擬尺度確定和模型參數(shù)確定方面有所差異。以生理發(fā)育時(shí)間而非作物實(shí)際生長(zhǎng)天數(shù)作為模擬尺度，消除了以實(shí)際天數(shù)作為模擬尺度造成的溫光差異對(duì)蓖麻生育進(jìn)程的影響，統(tǒng)一了各器官開(kāi)始生長(zhǎng)和結(jié)束生長(zhǎng)的時(shí)間節(jié)點(diǎn)以及器官間疊次生長(zhǎng)的時(shí)間間隔。在品種遺傳參數(shù)確定方面，對(duì)不同器官采取了不同的處理方法。對(duì)于節(jié)間生長(zhǎng)的模型參數(shù)處理較為系統(tǒng)，探明了各節(jié)間同伸規(guī)律，量化了各節(jié)間同步生長(zhǎng)及其與生理發(fā)育時(shí)間的關(guān)系，找出了各節(jié)間長(zhǎng)度和粗度的變化規(guī)律，確定了最優(yōu)條件下各品種節(jié)間長(zhǎng)度和粗度極值，在此基礎(chǔ)上，以節(jié)間長(zhǎng)極值和主莖直徑極值作為品種遺傳參數(shù)來(lái)計(jì)算其他模型參數(shù)，較好地體現(xiàn)了節(jié)間生長(zhǎng)過(guò)程中的同步關(guān)系。對(duì)于葉片、果穗和蒴果的模型參數(shù)處理則相對(duì)簡(jiǎn)化，將葉長(zhǎng)極值、葉寬極值、果穗長(zhǎng)極值和蒴果直徑極值直接確定為品種遺傳參數(shù)，經(jīng)驗(yàn)性和機(jī)理性有效地結(jié)合，使模型更加簡(jiǎn)便實(shí)用。

蓖麻分枝多，形態(tài)復(fù)雜，本研究構(gòu)建了地上部主莖器官形態(tài)建成模型。如何描述分枝器官的生長(zhǎng)動(dòng)態(tài)，是下一步要進(jìn)行的工作，模型也需更多品種和更大地域范圍的檢驗(yàn)。

4 結(jié)論

以生理發(fā)育時(shí)間為模擬尺度，構(gòu)建了蓖麻地上部主莖器官形態(tài)建成的預(yù)測(cè)模型，檢驗(yàn)結(jié)果表明，模型表現(xiàn)出較好的預(yù)測(cè)性、可靠性和一定的機(jī)理性。