白晨陽，何菡子，賈才華，李曉華，任奕林，葉俊，汪波，蒯婕，周廣生

機(jī)械收獲方式對油菜籽粒關(guān)鍵性狀的影響

白晨陽1，何菡子1，賈才華2，李曉華1，任奕林3，葉俊4，汪波1，蒯婕1，周廣生1

1華中農(nóng)業(yè)大學(xué)植物科學(xué)技術(shù)學(xué)院/農(nóng)業(yè)農(nóng)村部長江中游作物生理生態(tài)與耕作重點(diǎn)實(shí)驗(yàn)室，武漢 430070；2華中農(nóng)業(yè)大學(xué)食品科學(xué)技術(shù)學(xué)院，武漢 430070；3華中農(nóng)業(yè)大學(xué)工學(xué)院，武漢 430070；4湖北省耕地質(zhì)量與肥料工作總站，武漢 430070

【】收獲是油菜生產(chǎn)的關(guān)鍵環(huán)節(jié)之一，影響油菜籽粒產(chǎn)量、品質(zhì)及效益。通過人工模擬聯(lián)合收獲和分段收獲方式，分析不同收獲方式對油菜籽粒關(guān)鍵性狀的影響，以期為機(jī)械收獲方式選用及配套參數(shù)確定提供依據(jù)。選用含油量不同的2個(gè)中熟油菜品種，在湖北黃岡和襄陽設(shè)置聯(lián)合收獲不同收獲時(shí)期（終花后20—44 d（黃岡）和終花后23—47 d（襄陽），每隔3 d收獲1次）以及分段收獲不同割倒時(shí)間（終花后20、26和32 d（黃岡）和終花后23、29和35 d（襄陽））和后熟時(shí)長（3、6、9和12 d）試驗(yàn)，測定千粒重、含水量、品質(zhì)等關(guān)鍵指標(biāo)。聯(lián)合收獲兩試點(diǎn)、兩品種均在終花后38 d左右達(dá)籽粒生理成熟期，所需積溫約1 100℃，千粒重和含油量最大。終花后46 d左右，籽粒含水量和葉綠素含量可快速降至最適點(diǎn)，為機(jī)械化聯(lián)合收獲適宜時(shí)期。籽粒含水量降幅與日均溫及大氣壓極顯著正相關(guān)；葉綠素含量降幅與日均溫極顯著正相關(guān)，與降雨量顯著負(fù)相關(guān)。與D1、D2相比，分段收獲在D3條件下（終花后33 d左右）割倒后熟，籽粒千粒重、含油量、出油率、整齊度及油酸含量均值最高，葉綠素含量、亞油酸含量及含水量均值最低；后熟天數(shù)對籽粒含油量、出油率、千粒重、整齊度、油酸、亞油酸含量影響不顯著，后熟6—9 d時(shí)籽粒葉綠素與含水量可快速降至最適點(diǎn)，為適宜撿拾脫粒時(shí)期。籽粒含水量降幅與太陽輻射量極顯著正相關(guān)，與降雨量極顯著負(fù)相關(guān)。終花至收獲期間，平均日均溫19—22℃、降雨量1.0—3.5 mm，大氣壓1.45—1.75 kpa條件下聯(lián)合收獲在終花后46 d左右（生理成熟后8 d左右）收獲，分段收獲在終花后33 d左右（即生理成熟前5 d左右）割倒；平均相對濕度64%—80%、降雨量0.7—3.1 mm、太陽輻射量10—13 MJ·m-2條件下分段收獲后熟6—9 d撿拾脫粒最為適宜。此條件下，籽粒含水量均可降至15%左右、葉綠素含量均可降至10 mg·kg-1以下。且2種收獲方式在籽粒千粒重、含油量、出油率及整齊度等方面均無顯著差異；但與聯(lián)合收獲相比，分段收獲籽粒油酸、亞油酸含量更高，籽粒和油脂品質(zhì)更佳。

油菜；聯(lián)合收獲；分段收獲；籽粒關(guān)鍵性狀

0 引言

【研究意義】收獲是油菜生產(chǎn)最重要的環(huán)節(jié)之一，歐美油菜生產(chǎn)大國普遍采用機(jī)械收獲方式，油菜從種到收全程機(jī)械化生產(chǎn)，生產(chǎn)效率極高。中國受傳統(tǒng)小農(nóng)經(jīng)營方式的限制，油菜機(jī)械化生產(chǎn)尚未普及，其中，最大的難點(diǎn)就是機(jī)械收獲?！厩叭搜芯窟M(jìn)展】油菜機(jī)械收獲有聯(lián)合收獲和分段收獲2種方式。聯(lián)合收獲是在油菜角果完全成熟時(shí)利用聯(lián)合收割機(jī)一次性完成切割、脫離、清選等環(huán)節(jié)的收獲方式[1-3]；分段收獲則是指在油菜角果尚未完全成熟時(shí)，利用割曬機(jī)械將油菜割倒，晾曬后熟后用撿拾脫粒機(jī)進(jìn)行撿拾、脫粒的收獲方式。加拿大和美國以分段收獲為主，德國以聯(lián)合收獲為主[4]。目前，2種收獲方式在中國油菜各產(chǎn)區(qū)均有應(yīng)用[5]。聯(lián)合收獲機(jī)械一次作業(yè)，效率高，但對角果成熟一致性要求高[6-7]，導(dǎo)致其適收期短，一般僅3—5 d[8]。收獲稍早，植株含水量高，不利于機(jī)械作業(yè)，且損失率高而影響產(chǎn)量；收獲稍晚，角果易炸裂，損失率亦高而影響產(chǎn)量[8]。與聯(lián)合收獲相比，分段收獲可提前割倒，油菜能提早騰茬而不影響下茬作物生產(chǎn)[1-2,9]，但需機(jī)械2次作業(yè)，效率不及聯(lián)合收獲。分段收獲對適宜割倒時(shí)間亦有嚴(yán)格要求。過早會(huì)導(dǎo)致產(chǎn)量和品質(zhì)下降[10]；過遲，則茬口矛盾顯現(xiàn)，影響下茬作物生產(chǎn)?！颈狙芯壳腥朦c(diǎn)】目前，聯(lián)合收獲與分段收獲相關(guān)研究主要集中在作業(yè)水平、作業(yè)效率等方面，而2種收獲方式的適宜作業(yè)時(shí)期及適宜作業(yè)時(shí)期的2種收獲方式對油菜籽粒產(chǎn)量、品質(zhì)等關(guān)鍵指標(biāo)影響的研究尚不多見，導(dǎo)致中國各油菜產(chǎn)區(qū)在選用何種收獲方式上仍有較大爭論。【擬解決的關(guān)鍵問題】本研究人工模擬油菜聯(lián)合收獲與分段收獲，分析其對籽粒含水量、千粒重、整齊度、含油量、葉綠素含量及出油率等指標(biāo)的影響，以期為油菜機(jī)械收獲方式的選擇及適宜作業(yè)時(shí)期的確定提供依據(jù)。

1 材料與方法

1.1 材料

試驗(yàn)于2018年9月至2019年5月在湖北省黃岡農(nóng)業(yè)科學(xué)研究院（114°55′5"E，30°34′24"N；編號(hào)HG）及湖北省襄陽農(nóng)業(yè)科學(xué)研究院（111°50′43"E，31°46′52"N；編號(hào)XY）進(jìn)行。采用含油量分別為41.46%和49.04%（數(shù)據(jù)來源品種審定）的華油雜62（華中農(nóng)業(yè)大學(xué)選育，編號(hào)HZ62）和中雙11（中國農(nóng)業(yè)科學(xué)研究院油料作物研究所選育，編號(hào)ZS11）為試驗(yàn)材料。

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

聯(lián)合收獲采用裂區(qū)設(shè)計(jì)，品種為主區(qū)，不同收獲時(shí)間為副區(qū)（表1）。采用離地10 cm處人工割倒后直接脫粒的方式，模擬聯(lián)合機(jī)械收獲。共計(jì)18個(gè)處理，3次重復(fù)，54個(gè)小區(qū)。

分段收獲采用裂區(qū)設(shè)計(jì)，品種為主區(qū)，不同割倒時(shí)間為副區(qū)，后熟時(shí)間為副副區(qū)（表1）。共計(jì)24個(gè)處理，3次重復(fù)，72個(gè)小區(qū)。采用離地10 cm處人工割倒、后熟一定天數(shù)后人工脫粒的方式模擬分段機(jī)械收獲。

表1 2種收獲方式收獲割倒時(shí)間及后熟時(shí)長

HG：黃岡；XY：襄陽。下同 HG: Huanggang; XY: Xiangyang. The same as below

小區(qū)面積為10 m2（2 m×5 m），以600 kg·hm-2復(fù)合肥（15%-15%-15%）及15 kg·hm-2硼砂為基肥。采用20 cm的行距條播，黃岡與襄陽試點(diǎn)播期分別為2018年9月30日和2018年10月6日。出苗后開始間苗，4—5葉期定苗，密度為37.5×104株/hm2。7—9葉期追施75 kg·hm-2尿素（46%）。其他管理同常規(guī)。

油菜終花至角果發(fā)育成熟期，兩試驗(yàn)點(diǎn)氣象數(shù)據(jù)如圖1所示。在此期間，兩試點(diǎn)日均溫基本一致，黃岡試點(diǎn)日均相對濕度、大氣壓及降雨量均高于襄陽試點(diǎn)，但日均太陽輻射量低于襄陽試點(diǎn)，因此，黃岡試點(diǎn)油菜成熟度早于襄陽點(diǎn)，收獲時(shí)間提前3 d。

1.3 測定指標(biāo)與方法

1.3.1 千粒重及整齊度測定 用油菜千粒板數(shù)取1 000粒含水量為9%且均勻一致的籽粒稱重，記為千粒重；用萬深種子考種儀SC-G（杭州萬深檢測科技有限公司）測定（1 000±100）粒籽粒的表面積，利用變異系數(shù)=標(biāo)準(zhǔn)差/均值計(jì)算籽粒整齊度。變異系數(shù)越大，則籽粒整齊度越低；反之，則籽粒整齊度越高。

1.3.2 籽粒含水量測定 稱?。?00±5）g的新鮮籽粒樣品，105℃殺青30 min，80℃干燥至恒重并稱重，計(jì)算籽粒含水量[11]。

1.3.3 關(guān)鍵品質(zhì)指標(biāo)測定 在中國農(nóng)業(yè)科學(xué)研究院油料作物研究所，采用近紅外掃描法（NYDL-3000）測定籽粒中油脂、油酸和亞油酸含量；籽粒使用磨樣機(jī)（德國Retsch MM400）磨成粉，稱取2.000 g樣品，倒入裝有3顆不銹鋼球（直徑16 mm，304材質(zhì)）的50 ml定制鋼筒內(nèi)（鋼筒直徑34 mm，高度124 mm），加入30 ml萃取液（無水乙醇﹕石油醚=1﹕3），用橡膠塞塞好瓶口，至搖床振蕩1 h（頻率240次/min）后取出直立靜置10 min，采用250 ml布氏漏斗進(jìn)行普通過濾5 min，濾液收集于15 ml離心管內(nèi)，于665、625和705 nm處測定濾液吸光值，計(jì)算葉綠素含量[12]，公式如下：

圖1 黃岡及襄陽兩試驗(yàn)點(diǎn)的氣象數(shù)據(jù)

葉綠素含量（c，mg·kg-1）=k×Acorr×V/m/l式中，k為常數(shù)13；Acorr為修正吸光度，Acorr= A665-(A705+A625)/2；V為加入鋼筒萃取液的體積，30 ml；m為試樣質(zhì)量（g）；l為比色皿光徑，10 mm。

1.3.4 籽粒出油率測定 采用德國螺旋式榨油機(jī)將去雜曬干的油菜籽樣品進(jìn)行冷榨，每份菜籽樣品大約200 g。收集冷榨出的菜籽毛油和油菜粕，稱量并記錄冷榨菜籽毛油的質(zhì)量。菜籽冷榨毛油出油率計(jì)算公式：冷榨毛油出油率=菜籽冷榨毛油質(zhì)量/冷榨菜籽總質(zhì)量。

1.3.5 數(shù)據(jù)分析與整理 各處理數(shù)據(jù)以平均值表示。采用Microsoft Excel 2016、Origin 2018進(jìn)行數(shù)據(jù)整理與作圖。SPSS 25.0進(jìn)行數(shù)據(jù)分析，LSD法進(jìn)行差異顯著性分析，皮爾遜相關(guān)系數(shù)進(jìn)行相關(guān)性分析。

2 結(jié)果

2.1 機(jī)械收獲關(guān)鍵性狀

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，兩試點(diǎn)、兩品種千粒重均呈先增后略降趨勢（圖2-A），黃岡與襄陽試點(diǎn)千粒重分別于終花后35和41 d達(dá)峰值，為籽粒生理成熟期。與生理成熟期相比，黃岡試點(diǎn)HZ62、ZS11千粒重在終花后44 d分別下降4.7%和8.3%，襄陽試點(diǎn)兩品種千粒重在終花后47 d分別下降8.2%和3.7%。兩試點(diǎn)均為ZS11千粒重高于HZ62。

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，含水量快速下降（圖2-A）。生理成熟期千粒重最大，但含水量仍在40%—50%。黃岡試點(diǎn)，ZS11和HZ62分別于終花后41和44 d含水量降至14.10%和10.58%；襄陽試點(diǎn)，ZS11和HZ62于終花后47 d含水量分別降至15.50%和15.97%。籽粒含水量降幅與日均溫和大氣壓極顯著正相關(guān)（表2）。

不同字母表示p＜0.05水平上差異性顯著。A：聯(lián)合收獲方式；B：分段收獲方式。下同

籽粒生理成熟前，干物質(zhì)持續(xù)累積、千粒重增加。分段收獲方式下，隨割倒時(shí)間推遲千粒重顯著增加（圖2-B），兩試點(diǎn)、兩品種千粒重均值均為D3處理高于D2、D1處理，分別高13.48%和35.98%。割倒后熟過程中，隨后熟時(shí)間延長，千粒重呈先增后略降趨勢。D1、D2處理均于后熟3—6 d時(shí)達(dá)最大值；D3處理?xiàng)l件下除襄陽試點(diǎn)ZS11外，其余處理不同后熟天數(shù)間均無顯著性差異。

分段收獲方式下，隨割倒時(shí)間推遲及后熟時(shí)間延長，籽粒含水量均快速下降（圖2-B），兩試點(diǎn)、兩品種含水量均值均為D3處理低于D2、D1處理，后熟0 d時(shí)分別低15.92%和22.75%，后熟12 d時(shí)分別低73.17%和69.13%。且后熟過程中，D3處理含水量降幅大于D2、D1處理。D3處理?xiàng)l件下，黃岡及襄陽試點(diǎn)兩品種籽粒含水量分別于后熟9—12 d、6—9 d時(shí)降至15%左右。含水量降幅與相對濕度及降雨量極顯著負(fù)相關(guān)、與太陽輻射量極顯著正相關(guān)（表2）。

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，兩試點(diǎn)、兩品種含油量均呈先增后降趨勢（圖2-A）。黃岡與襄陽試點(diǎn)分別于終花后29和35 d達(dá)峰值。分段收獲方式下，隨割倒時(shí)間推遲，籽粒含油量逐漸增加（圖2-B）。兩試點(diǎn)、兩品種含油量均值均為D3處理高于D2、D1處理，分別高5.69%和17.13%；后熟時(shí)間對籽粒含油量影響在不同割倒處理間存在差異。D1、D2處理?xiàng)l件下，隨后熟時(shí)間延長，籽粒含油量呈先升后降趨勢，后熟3—6 d可達(dá)最大值；D3處理下，隨后熟時(shí)間延長，含油量呈下降趨勢，但不同后熟天數(shù)間差異均未達(dá)顯著水平。

2.2 籽粒商品性

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，兩試點(diǎn)、兩品種葉綠素含量均呈先下降后穩(wěn)定趨勢（圖3-A），黃岡及襄陽試點(diǎn)分別于終花后35 d、41 d達(dá)到穩(wěn)定。葉綠素含量降幅與日均溫極顯著正相關(guān)，與降雨量顯著負(fù)相關(guān)（表2）。聯(lián)合收獲方式下，黃岡試點(diǎn)終花后44 d、襄陽試點(diǎn)終花后47 d時(shí)葉綠素含量趨于一致且均為最低點(diǎn)。分段收獲方式下，隨割倒時(shí)間推遲，葉綠素含量顯著下降（圖3-B）。兩試點(diǎn)、兩品種葉綠素含量均值均為D3處理低于D2、D1處理，后熟0 d時(shí)分別低68.74%和81.11%，后熟12 d時(shí)分別低61.45%和75.12%。隨后熟時(shí)間的延長，籽粒葉綠素含量顯著降低。兩試點(diǎn)、兩品種葉綠素降幅均表現(xiàn)為D3

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，襄陽試點(diǎn)ZS11籽粒整齊度呈增加趨勢，終花后47 d達(dá)最低值；黃岡試點(diǎn)兩品種及襄陽試點(diǎn)HZ62均呈先增后降趨勢，分別于終花后38和44 d時(shí)達(dá)峰值（圖3-A）。兩試點(diǎn)ZS11整齊度均高于HZ62；且黃岡試點(diǎn)品種間差異更大。分段收獲方式下，隨割倒時(shí)間推遲，籽粒整齊度增加（圖3-B）。兩試點(diǎn)、兩品種均為D3處理高于D2、D1處理，其中D3處理籽粒面積變異系數(shù)均值較D2、D1處理分別低5.70%、14.83%。后熟過程中，黃岡試點(diǎn)兩品種D3處理籽粒整齊度均呈下降趨勢，后熟0 d時(shí)達(dá)最大值；黃岡試點(diǎn)HZ62 D2處理呈先增后降趨勢，后熟3 d時(shí)達(dá)最大值；其余處理均呈增加趨勢，后熟3—6 d時(shí)達(dá)最大值。

表2 籽粒含水量及葉綠素含量下降幅度與氣象數(shù)據(jù)相關(guān)性

*表示＜0.05，**表示＜0.01。由于表中含水量及葉綠素含量下降幅度指每3天含水量的變化，所有氣象數(shù)據(jù)為3 d內(nèi)各氣象因子的總和

* Significant difference at＜0.05. ** Significant difference at＜0.01. Because the moisture content and chlorophyll content declines in the table refers to the change of moisture content every 3 days, all the meteorological data are the sum of the three days

圖3 黃岡襄陽試點(diǎn)2種收獲方式下籽粒葉綠素含量、籽粒整齊度變化

2.3 籽粒加工品質(zhì)

2.3.1 油酸及亞油酸含量 聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，兩試點(diǎn)、兩品種油酸含量均呈先增后降趨勢（圖4-A）。襄陽試點(diǎn)HZ62于終花后29 d達(dá)峰值；黃岡試點(diǎn)兩品種及襄陽試點(diǎn)ZS11均于終花后38 d達(dá)峰值。兩試點(diǎn)，兩品種間的籽粒油酸含量差異均不顯著。分段收獲方式下，隨割倒時(shí)間延遲，籽粒油酸含量呈增加趨勢（圖4-B）。兩試點(diǎn)、兩品種油酸含量均值均為D3處理高于D2、D1處理，分別高2.39%和4.89%；后熟時(shí)長亦影響籽粒油酸含量，襄陽試點(diǎn)ZS11在D1、D2、D3割倒后，隨后熟時(shí)間延長，油酸含量增加，其余處理則呈下降趨勢，其中襄陽試點(diǎn)HZ62 D1處理于后熟3 d時(shí)達(dá)最低值，其余處理不同后熟天數(shù)間均無顯著性差異。兩試點(diǎn)、兩品種規(guī)律一致。

圖4 不同收獲方式下油菜籽粒油酸與亞油酸含量變化

聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，黃岡試點(diǎn)ZS11籽粒亞油酸含量在終花后44 d內(nèi)均持續(xù)下降；其余處理均呈先降后增趨勢（圖4-A），黃岡HZ62及ZS11分別于終花后38和44 d達(dá)最低值；襄陽試點(diǎn)兩品種分別于35和41 d達(dá)最低值。兩試點(diǎn)均為HZ62高于ZS11。分段收獲方式下，隨割倒時(shí)間推遲，亞油酸含量呈下降趨勢（圖4-B）。兩試點(diǎn)、兩品種亞油酸含量均值均為D3處理低于D2、D1處理，分別低5.29%和8.34%；后熟時(shí)長對籽粒亞油酸含量亦有影響。D1、D2條件下割倒，隨后熟時(shí)間延長，籽粒亞油酸含量均呈增加趨勢，后熟9—12 d均可達(dá)到最大值；D3條件下割倒，均呈先增加后降低趨勢，后熟6—9 d可達(dá)到最大值。

2.3.2 出油率 聯(lián)合收獲方式下，隨著收獲時(shí)間的推遲，籽粒出油率呈先升后降趨勢（圖5-A）。黃岡及襄陽試點(diǎn)出油率分別于終花后35和38 d達(dá)到最高點(diǎn)，略晚于含油率峰值。兩試點(diǎn)出油率均為ZS11高于HZ62，且黃岡試點(diǎn)差異更大。分段收獲方式下，隨割倒時(shí)間推遲，籽粒出油率增加（圖5-B）。兩試點(diǎn)、兩品種出油率均值均為D3處理高于D2、D1處理，分別高11.02%和41.48%。不同后熟時(shí)長對籽粒出油率影響顯著。D1、D2條件下，隨后熟時(shí)間延長，籽粒出油率量總體呈先升后降趨勢，后熟3—6 d達(dá)最大值；D3條件下，隨后熟時(shí)間延長，出油率呈下降趨勢，但處理間差異均未達(dá)顯著水平。

圖5 不同收獲方式下油菜籽粒出油率含量變化

2.4 適宜作業(yè)時(shí)期2種機(jī)械作業(yè)方式對籽粒關(guān)鍵指標(biāo)的影響

綜合機(jī)械收獲關(guān)鍵性狀（籽粒千粒重及含水量），籽粒商品性（籽粒整齊度及葉綠素含量）和加工品質(zhì)（籽粒出油率）在不同收獲方式間的差異，確定了2種收獲方式作業(yè)的適宜時(shí)期。聯(lián)合收獲在生理成熟后8 d左右開始機(jī)械作業(yè)；分段收獲在生理成熟前5 d左右割倒、后熟6—9 d開始撿拾脫粒。在最適作業(yè)時(shí)間條件下，2種收獲方式的籽粒千粒重、含水量、含油量、出油率、籽粒整齊度以及葉綠素含量均無顯著差異；分段收獲方式下籽粒油酸和亞油酸含量均略高于聯(lián)合收獲方式，其中，兩試點(diǎn)ZS11處理間差異性顯著（表3）。

表3 聯(lián)合收獲機(jī)械化收獲最適點(diǎn)與分段收獲最適點(diǎn)比較

不同字母表示＜0.05水平上差異性顯著。表中黃岡試點(diǎn)的CK指終花后44 d所取樣品；襄陽試點(diǎn)的CK指終花后47 d所取樣品

Different letters indicate significant differences at the＜0.05 level. The CK refers to the sample taken 44 days after the final flower in Huanggang; The CK refers to the sample taken 47 days after the flowers in Xiangyang

3 討論

歐美等國家（地區(qū)）已實(shí)現(xiàn)油菜生產(chǎn)機(jī)械化，甚至自動(dòng)化、智能化[13]，而中國油菜生產(chǎn)機(jī)械化水平較低且效益不高[14]。推廣油菜機(jī)械化收獲技術(shù)是提高生產(chǎn)效益的必要途徑。產(chǎn)量是油菜生產(chǎn)效益的基礎(chǔ)，千粒重、含水量以及含油量是機(jī)械化收獲關(guān)鍵性狀；籽粒整齊度和葉綠素含量是重要商品指標(biāo)；出油率是加工企業(yè)的關(guān)注點(diǎn)；油酸、亞油酸是消費(fèi)者關(guān)注的重要營養(yǎng)指標(biāo)。機(jī)械收獲方式及作業(yè)時(shí)間對瓜類及辣椒種子產(chǎn)量、質(zhì)量等關(guān)鍵指標(biāo)影響顯著[15-16]。本研究針對聯(lián)合收獲和分段收獲2種方式，綜合考慮收獲方式及作業(yè)時(shí)間對籽粒關(guān)鍵指標(biāo)的影響，以期為收獲方式的選擇及適宜參數(shù)的確定提供依據(jù)。本試驗(yàn)選在湖北省黃岡和襄陽兩地開展，兩地油菜生育進(jìn)程差異較大。黃岡地區(qū)油菜生育進(jìn)程快，與長江上游接近；襄陽地區(qū)油菜生育進(jìn)程較慢，與長江下游接近。且本試驗(yàn)選用含油量差異較大的華油雜62與中雙11為試驗(yàn)材料。因此，本試驗(yàn)點(diǎn)及試驗(yàn)材料的選擇具有較強(qiáng)代表性，基于此得到的聯(lián)合收獲及分段收獲適宜作業(yè)參數(shù)亦有應(yīng)用價(jià)值，可為油菜機(jī)械收獲提供技術(shù)參數(shù)。

籽粒發(fā)育是油菜籽產(chǎn)量、品質(zhì)形成的重要環(huán)節(jié)。油菜籽粒發(fā)育分3個(gè)階段，第一階段籽粒重量較輕，主要成分為淀粉及乙醇等化合物；第二階段，籽粒重量增加，油脂和蛋白質(zhì)累積；第三階段籽粒重量繼續(xù)增加，但成分不再變化[17-18]，籽粒干重最大時(shí)達(dá)生理成熟[19]。本研究中，油菜終花后20 d時(shí)油脂大量合成使得籽?？焖偬畛洌蚜８芍丶罢R度顯著增加，含水量及葉綠素含量快速下降；終花后35 d左右油脂累積結(jié)束；終花后38 d左右干物質(zhì)累積結(jié)束，籽粒千粒重最大，達(dá)生理成熟期，所需有效積溫1 100℃。后熟對提高籽粒產(chǎn)量和品質(zhì)亦具有重要意義[20-21]。D1、D2處理在后熟過程中含油量、千粒重及整齊度略增，有后熟效應(yīng)；但D3處理后熟效應(yīng)不明顯，可能由于終花后33 d左右（D3），籽粒成熟度高，且植株已處衰老狀態(tài)，后熟過程中，籽粒含水量及葉綠素含量快速下降，光合能力減弱，導(dǎo)致儲(chǔ)存物質(zhì)合成速率下降[22]。

籽粒含水量是判斷適宜機(jī)械作業(yè)時(shí)期的重要指標(biāo)。油菜籽粒小，收獲時(shí)含水量過低會(huì)導(dǎo)致?lián)p失率較高而影響產(chǎn)量[23]；但含水量偏高，即使籽粒已達(dá)生理成熟期，也會(huì)因收獲時(shí)籽粒破碎或脫粒困難導(dǎo)致產(chǎn)量下降[24]。含水量降至15%左右，方可收獲[8,25]。本試驗(yàn)中，聯(lián)合收獲條件下，黃岡與襄陽試點(diǎn)分別在終花后41—44 d和47 d，均為生理成熟后6—9 d含水量降至適宜。分段收獲方式下，D3處理?xiàng)l件下，黃岡試點(diǎn)HZ62與ZS11分別于后熟9和12 d，襄陽試點(diǎn)分別于后熟6和9 d含水量降至15%左右。籽粒含水量降幅與氣象因素密切相關(guān)[26-27]，但影響不同收獲方式下籽粒含水量降幅的氣象因子略有差異。溫度及大氣壓高，有利于聯(lián)合收獲方式下籽粒含水量的下降；降雨量小、太陽輻射量高、空氣相對濕度低，有利于分段收獲后熟過程中籽粒含水量的下降；因此，聯(lián)合收獲可據(jù)當(dāng)?shù)販囟?、大氣壓高低預(yù)測機(jī)械作業(yè)時(shí)間，分段收獲據(jù)降雨量、太陽輻射量及空氣相對濕度調(diào)整撿拾脫粒時(shí)間。

葉綠素含量是判斷籽粒成熟度的重要指標(biāo)。籽粒發(fā)育過程中較高葉綠素含量有利油脂積累[28]，但成熟收獲時(shí)葉綠素含量較高則降低籽粒質(zhì)量、加速油脂氧化劣敗[29]、縮短籽粒壽命[30]。加拿大將籽粒進(jìn)行分級：一級、二級及三級籽粒的葉綠素含量分別低于10、43和101 mg·kg-1[31]。籽粒成熟及后熟過程中葉綠素含量均顯著下降，2種收獲方式最佳作業(yè)時(shí)間下籽粒葉綠素含量無顯著性差異，均低于10 mg·kg-1，達(dá)到加拿大的一級籽粒標(biāo)準(zhǔn)。葉綠素分解受基因型和環(huán)境的影響，本試驗(yàn)中，聯(lián)合收獲條件下葉綠素含量降幅與日均溫極顯著正相關(guān)，與降雨量顯著負(fù)相關(guān)，因此，在日均溫較高、降雨量較小的地區(qū)采用聯(lián)合收獲可獲得低葉綠素含量的優(yōu)質(zhì)油菜籽。此外，葉綠素分解受多種內(nèi)源物質(zhì)的調(diào)控，其中，（reactive oxygen species，ROS）在葉綠素降解方面發(fā)揮著重要作用[32]，葉綠素含量的降低與葉綠素結(jié)合蛋白有關(guān)[33]，過量的ROS攻擊生物大分子，導(dǎo)致葉綠素結(jié)合蛋白喪失其功能，葉綠素含量顯著降低。分段收獲割倒時(shí)，植物遭受機(jī)械損傷，ROS大量積累[34-35]，葉綠素的降解速率加快。因此，籽粒成熟前期聯(lián)合收獲處理葉綠素降幅顯著低于分段收獲，聯(lián)合收獲過早收獲可能導(dǎo)致綠籽率增加。而籽粒發(fā)育后期，由于植株衰老也會(huì)導(dǎo)致ROS的積累[34]，進(jìn)而促進(jìn)的葉綠素降解。

此外，籽粒發(fā)育過程中隨著含油量的增加，出油率顯著增加，不飽和脂肪酸含量增加，飽和脂肪酸含量下降[36]，油脂產(chǎn)量和品質(zhì)均顯著提高。本試驗(yàn)結(jié)果表明分段收獲方式下籽粒油酸及亞油酸含量均優(yōu)于聯(lián)合收獲，且油酸及亞油酸含量累積速率均大于聯(lián)合收獲。油酸由硬脂酸經(jīng)脫不飽和產(chǎn)生，亞油酸由油酸脫不飽和產(chǎn)生，SAD、FAD2、FAD3分別是控制油酸和亞油酸合成及降解的關(guān)鍵酶[37]。前人研究表明，非生物脅迫條件下，如缺氧、低溫以及脫水等可誘導(dǎo)表達(dá)量上調(diào)，FAD表達(dá)量下調(diào)，油酸含量顯著增加[38-39]。在水稻中亞油酸含量也會(huì)發(fā)生一定的積累[40]。本試驗(yàn)中，分段收獲在割倒后，由于植株受到外部機(jī)械損傷，可能誘導(dǎo)了的表達(dá)并下調(diào)了FAD和FAD的表達(dá)，從而促進(jìn)了油酸和亞油酸的累積。

綜上，油菜終花后46 d左右，全田角果基本干枯變黃、籽粒含水量及葉綠素含量均快速下降，分別降至15%左右及10 mg·kg-1以下，可聯(lián)合收獲；終花后33 d左右，主莖角果基本褪綠呈枇杷黃，籽粒干物質(zhì)累積緩慢，可割倒后熟，后熟過程中籽粒含水量及葉綠素含量快速下降，6—9 d后分別降至15%左右及10 mg·kg-1以下，可撿拾脫粒。但如遇降雨，籽粒含水量降幅受到明顯影響，2種收獲方式的適宜作業(yè)時(shí)期均應(yīng)適度順延。在最適條件下，2種收獲方式均可獲得高產(chǎn)優(yōu)質(zhì)油菜籽。2種收獲方式下籽粒千粒重及含油量無顯著性差異，但分段收獲籽粒及油脂品質(zhì)更佳。本試驗(yàn)為油菜機(jī)械收獲提供了技術(shù)參數(shù)及理論依據(jù)，但未考慮不同收獲方式對機(jī)收損失率的影響，后續(xù)可開展相關(guān)試驗(yàn)，測定不同收獲方式對產(chǎn)量的影響。此外，成熟后的干籽粒內(nèi)含有大量mRNA，后熟過程中多數(shù)mRNA豐度降低，蛋白質(zhì)含量增加[41]，這些變化是否影響油菜籽粒萌發(fā)與抗逆性？后熟過程中葉綠素降解較快受哪些物質(zhì)調(diào)控？油脂色澤的差異是否與種皮內(nèi)花青素等色素含量有關(guān)？2種收獲方式小分子代謝又有哪些差異等諸多問題還有待進(jìn)一步探索。

4 結(jié)論

終花至收獲期間，平均日均溫19—22℃、降雨量1.0—3.5 mm，大氣壓1.45—1.75 kpa條件下聯(lián)合收獲在生理成熟后8 d左右收獲，分段收獲在生理成熟前5 d左右割倒；平均相對濕度64%—80%、降雨量0.7—3.1 mm、太陽輻射量10—13 MJ·m-2條件下分段收獲后熟6—9 d撿拾脫粒最為適宜。此條件下，籽粒含水量均可降至15%左右、葉綠素含量低于10 mg·kg-1。且2種收獲方式在籽粒千粒重、含油量、油酸、亞油酸、整齊度及出油率方面均無顯著差異。但與聯(lián)合收獲相比，分段收獲籽粒油酸、亞油酸含量更高，籽粒和油脂品質(zhì)更佳。

[1] 涂安富, 吳崇友, 盧晏. 加拿大油菜的收獲管理(分段與聯(lián)合收割)(一). 農(nóng)業(yè)裝備技術(shù), 2010, 36(1): 30-33.

TU A F, WU C Y, LU Y. Harvesting management of rapeseed in Canada (section and combined harvesting) (1). Agricultural Equipment Technology, 2010, 36(1): 30-33. (in Chinese)

[2] 涂安富, 吳崇友, 盧晏. 加拿大油菜的收獲管理(分段與聯(lián)合收割)(二). 農(nóng)業(yè)裝備技術(shù), 2010, 36(2): 49-51.

TU A F, WU C Y, LU Y. Harvest management of rapeseed in Canada (section and combined harvesting) (2). Agricultural Equipment Technolog, 2010, 36(2): 49-51. (in Chinese)

[3] 羅海峰, 湯楚宙, 官春云, 吳明亮. 油菜機(jī)械化收獲研究進(jìn)展. 農(nóng)機(jī)化研究, 2015, 37(1): 1-8.

LUO H F, TANG C Z, GUAN C Y, WU M L. Research progress on mechanized harvesting of rapeseed. Journal of Agricultural Mechanization Research, 2015, 37(1): 1-8. (in Chinese)

[4] 秦順創(chuàng). 油菜機(jī)械化技術(shù)發(fā)展影響因素分析及推廣策略研究[D]. 杭州: 浙江大學(xué), 2015.

QIN S C. Analysis of factors affecting the development of rape mechanization technology and research on its promotion strategy[D]. Hangzhou: Zhejiang University, 2015. (in Chinese)

[5] SHI Z X, WU M L, YANG W M, SHEN Y F, WANG X S, HE J, WANG J Y. Research status and development measures of rape segment harvester in China. Agricultural Engineering, 2015, 5(5): 1-4.

[6] 吳崇友, 肖圣元, 金梅. 油菜聯(lián)合收獲與分段收獲效果比較. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(17): 10-16.

WU C Y, XIAO S Y, JIN M. Comparison of combined harvesting and two-staged harvesting of rapeseed. Journal of Agricultural Engineering, 2014, 30(17): 10-16. (in Chinese)

[7] 周可金, 官春云, 肖文娜, 譚太龍. 催熟劑對油菜角果光合特性、品質(zhì)及產(chǎn)量的影響. 作物學(xué)報(bào),2009, 35(7): 1369-1373.

ZHOU K J, GUAN C Y, XIAO W N, TAN T L. Effects of ripening agent on photosynthetic characteristics, quality and yield of rapeseed carambola.Acta Agronomica Sinica, 2009, 35(7): 1369-1373. (in Chinese)

[8] 左青松, 黃海東, 曹石, 楊士芬, 廖慶喜, 冷鎖虎, 吳江生, 周廣生. 不同收獲時(shí)期對油菜機(jī)械收獲損失率及籽粒品質(zhì)的影響. 作物學(xué)報(bào), 2014, 40(4): 650-656.

ZUO Q S, HUANG H D, CAO S, YANG S F, LIAO Q X, LENG S H, WU J S, ZHOU G S. Effects of different harvesting periods on mechanical harvesting loss rate and grain quality of rapeseed. Acta Agronomica Sinica, 2014, 40(4): 650-656. (in Chinese)

[9] WU C Y, XIAO S Y, JIN M. Comparation on rape combine harvesting and two-stage harvesting. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(17): 10-16.

[10] 伊淑麗, 梁穎, 代柳亭, 諶利, 柴友榮, 李加納. 高溫對甘藍(lán)型油菜籽粒后熟相關(guān)特性的影響. 西南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 30(2): 48-50.

YI S L, LIANG Y, DAI L T, CHEN L, CHAI Y R, LI J N. Effects of high temperature on related characteristics of rape seed rape post- ripening. Journal of Southwest University (Natural Science Edition), 2008, 30(2): 48-50. (in Chinese)

[11] ISTA.Chapter 9: Determination of moisture content. International Rules for Seed Testing. Baaserdorf, Switzerland. The International Seed Testing Association (ISTA), 2017.

[12] WARD K A, SCARTH R, DAUN J, THORSTEINSON C T. A comparison of high-performance liquid chromatography and spectrophotometry to measure chlorophyll in canola seed and oil. Journal of the American Oil Chemists’ Society, 1994, 71(9): 931-934.

[13] XU Q C, ZHOU Y. Consideration on accelerating the development of rape harvesting mechanization. Chinese Agricultural Mechanization, 2008, 18(2): 54-56.

[14] WU C Y, YI Z Y. Selection of technology route of oilseed rape mechanization in entire production proceeding in China. Chinese Agricultural Mechanization, 2009, 3(2): 3-6.

[15] ALAN ?, ESER B. The effect of fruit maturity and post-harvest ripening on seed quality in hot and conic pepper cultivars. Seed Science and Technolog, 2008, 36(2): 467-474.

[16] ALOHO K P, JOHNSON O A. The effects of leaf colour at fruit harvest and fruit after-ripening duration on (Naudin.) seed quality. Journal of Biology Agriculture & Healthcare, 2013, 3(2): 190-200.

[17] NORTON G, HARRIS J F. Compositional changes in developing rape seed (L.). Planta, 1975, 123(2): 163-174.

[18] NI F, LIU J H, ZHANG J, KHAM M N, LUO T, XU Z H, HU L Y. Effect of soluble sugar content in silique wall on seed oil accumulation during the seed-filling stage in. Crop and Pasture Science,2018, 69(12): 1251-1263.

[19] TEKRONY D M, EGLI D B, BALLES J, PFEIFFER T, FELLOWS R J. Physiological maturity in soybean1. Agronomy Journal, 1979, 71(5): 771.

[20] SUN S R, CHEN Y L, ZHAN X M, ZENG K, LI Y. The effect of different fruit age and post-ripening period on the quality index of squash seed.Acta Agriculturae Boreali-sinica, 2004, 19(4): 29-32.

[21] BLAY E T, DANQUAH E Y, ABABU A.Effect of time of harvest, stage of fruit ripening, and post-harvest ripening on seed yield and germinability of local garden egg (Radii). Ghana Journal of Agricultural Science, 1999, 32(2): 159-167.

[22] AARON F, RUTHIE A, HADAR L, ITZHAK O, EWA U W, ALISDAIR R F, GAD G.seed development and germination is associated with temporally distinct metabolic switches. Pant Physiology, 2006, 142: 839-854.

[23] 盧晏, 吳崇友, 金誠謙, 涂安富. 油菜機(jī)械化收獲方式的選擇. 農(nóng)機(jī)化研究, 2008(11):240-242, 245.

LU Y, WU C Y, JIN C Q, TU A F. Choice of mechanized harvesting methods for rape. Agricultural Mechanization Research, 2008(11): 240-242, 245. (in Chinese)

[24] HU Q, HUA W, YIN Y, ZHANG X K, LIU L J, SHI J Q, ZHAO Y G, QIN L, CHEN C. Rapeseed research and production in China. The Crop Journal, 2017, 5(2): 35-43.

[25] ZHANG M, WU C Y, JIN M, MU S L, LIANG S N, TANG Q. Effects of harvesting method and date on yield loss and seed quality of rapeseed. Chinese Journal of Oil Crop Sciences,2019, 4(3): 166-174.

[26] SINGH B, LINVILL D E. Field drying of navy beans in the harvesting period. Transactions of the Asae, 1977, 20(2): 228-231.

[27] SAMARAH N H, MULLEN R E, ALQUDAH A M. An index to quantify the relationship of seed moisture loss rate to seed desiccation tolerance in common vetch. Seed Technology, 2008, 30(1): 94.

[28] HUA S J, CHEN Z H, ZHANG Y F, YU H S, LIN B G, ZHANG D Q. Chlorophyll and carbohydrate metabolism in developing silique and seed are prerequisite to seed oil content ofL.. Botanical Studies, 2014, 55(1): 34.

[29] APARICIO-RUIZ R, MINGUZE-MOSQUERAO M I, GANDUL- ROJAS B. Thermal degradation kinetics of chlorophyll pigments in virgin olive oils:1. Compounds of series a. Journal of Agricultural and Food Chemistry, 2010, 58(10): 6200-6208.

[30] NAKAJIMA S, ITO H, TANAKA R, TANAKA A. Chlorophyll b reductase plays an essential role in maturation and storability ofseeds.Plant Physiology, 2012, 160(1): 261-273.

[31] Canadian Grain Commission. Quality of western Canadian Canola. 2018. http://www.grainscanada.gc.ca/oggg-gocg/ggg-gcg-eng. htm.

[32] SONG Y, YANG C G, GAO S, ZHANG W, LI L, KUAI B K. Age-triggered and dark-induced leaf senescence require the bHLH transcription factors PIF3, 4, and 5. Molecular Plant, 2014, 7(12): 1776-1787.

[33] YANG S, ZENG X, LI T, LIU M, ZHANG S C, GAO S J, WANG Y Q, PENG C L, LI L, YANG C W. AtACDO1, an ABC1-like kinase gene, is involved in chlorophyll degradation and the response to photooxidative stress in. Journal of Experimental Botany, 2012, 63(10): 3959-3973.

[34] ROSENVASSER S, MAYAK S, FRIEDMAN H. Increase in reactive oxygen species (ROS) and in senescence-associated gene transcript (SAG) levels during dark-induced senescence ofcuttings, and the effect of gibberellicacid. Plant Science, 2006, 170(4): 873-879.

[35] SEOL Y B, KIM J, PARK S H, CHANG H Y. Atmospheric pressure pulsed plasma induces cell death in photosynthetic organs via intracellularly generated ROS. Scientific Reports, 2017, 7(1): 589.

[36] BAUD S, BOUTIN J P, MIQUEL M, LEPINIEC L, ROCHAT C. An integrated overview of seed development inecotype WS. Plant Physiology and Biochemistry (Paris), 2002, 40(2): 151-160.

[37] FOFANA B, CLOUTIER S, DUGUID S, CHING J, RAMPITSCH C. Gene expression of stearoyl-ACP desaturase and Δ12 fatty acid desaturase 2 is modulated during seed development of flax (). Lipids, 2006, 41(7): 705-712.

[38] KLINKENBERG J, FAIST H, SAUPE S, LAMBERTS S, KRISCHKE M, STINGL N, FEKETE A, MUELLER M J, FEUSSNER I, HEDRICH R, DEEKEN R. Two fatty acid desaturases, SAD6and FAD3, are involved in drought and hypoxia stress signaling incrown galls. Plant Physiology, 2013, 164(2): 570-583.

[39] YUAN S W, WU X L, LIU Z H, LUO H B, HUANG R Z. Abiotic stresses and phytohormones regulate expression ofFADgene in. Journal of Integrative Agriculture, 2012(1): 62-72.

[40] VARTAPETIAN B B, MAZLIAK P, LANCE C. Lipid biosynthesis in rice coleoptiles grown in the presence or in the absence of oxygen. Plant Science Letters,1978, 13(4): 321-328.

[41] HOLDSWORTH M J, BENTSINK L, SOPPE W J J. Molecular networks regulatingseed maturation, after-ripening, dormancy and germination. New Phytologist, 2008, 179(1): 33-54.

effect of the mechanical harvesting methods on the key traits of rapeseed

Bai Chenyang1, He Hanzi1, Jia Caihua2, Li Xiaohua1, Ren Yilin3, Ye Jun4, Wang Bo1, Kuai jie1, Zhou Guangsheng1

1College of Plant Science and Technology, Huazhong Agricultural University/Key Laboratory of Crop Ecophysiology and Farming System in the Middle Reaches of the Yangtze River, Ministry of Agriculture, Wuhan 430070;2College of food Science and Technology, Huazhong Agricultural University, Wuhan 430070;3College of engineering, Huazhong Agricultural University, Wuhan 430070;4Hubei Provincial Cultivated Land Quality and Fertilizer Station, Wuhan 430070

【】Harvesting is one of the key links in rapeseed production, which affects the yield and quality. In this experiment, artificially simulated combined harvesting and two-stage harvesting were used and analyzed the impacts on key traits of rapeseeds. This study provides a basis for the selection of mechanical harvesting methods and the determination of supporting parameters.【】Two middle-ripening rapeseed varieties with different oil content were selected, and different harvesting periods were used. For combined harvesting the period is 20-44 d in Huanggang and 23-47 d in Xiangyang, and the two-stage harvesting is practiced on different time 20, 26, and 32 d in Huanggang while, 23, 29, 35 d in Xiangyang, after final flowering, whereas 3, 6, 9, and 12 d post-ripening, in Huanggang and Xiangyang, Hubei. Furthermore, key indicators such as1000-grain weight, water content, quality, etc were measured and analyzed. 【】The two pilots and two varieties of the combined harvesting reached the physiological maturity about 38 days after final flowering, the required accumulated temperature was about 1 100℃, and the thousand-grain weight and oil content were the largest. Approximately 46 days after the final flowering, the seed moisture and chlorophyll content rapidly decreased to the optimum point, which is a suitable period for mechanized combined harvesting. Moreover, reduction in seed moisture content is highly significantly positively correlated with daily average temperature and atmospheric pressure, while decreased chlorophyll content is highly significantly positively correlated with daily average temperature, but significantly negatively correlated with rainfall. Under D3 (about 33 d after final flowering), the highest average content of 1000-grain weight, oil contents, rate of oil yield, uniformity and oleic acid, and the lowest average content of chlorophyll, linoleic acid and moisture under the two-stage harvest while comparing D3 with D1 and D2, Our study showed that the number of days of post-ripening has no significant effect on the oil content, oil yield, thousand-grain weight, uniformity, oleic acid, and linoleic acid content, but the moisture content rapidly decreased to the optimum point after 6-9 days, which is a suitable time for picking and threshing. The seed moisture content decreased and significantly positively correlated with solar radiation but highly significantly negatively correlated with rainfall. 【】During the period of final flowering to the harvesting, the average daily temperature is 19-22℃, the rainfall is 1.0-3.5 mm, and the atmospheric-pressure is 1.45-1.75 kpa. However, under prevailing climatic conditions, the combined harvesting is done on approximately 46th day after final flowering (which is about 8 days after physiological maturity) and the two-staged harvesting is used to cut down on about 33rd day after the final flowering (which is about 5 days before the physiological maturity), then post-ripening practice is carried out (average relative humidity is 64%-80%, the rainfall is 0.7-3.1 mm, and the solar radiation is 10-13 MJ·m-2), to picked the harvested crop after 6-9 days of two-stage harvesting. Under this condition, the moisture content of seed can be reduced to about 15%, and the chlorophyll content can be reduced to below 10 mg·kg-1. There was no significant difference between the two harvesting methods in terms of 1000-grain weight, oil content, rate of oil yield and seed uniformity, however, the oleic acid and linoleic acid contents of seed were higher, seed and oil quality were better in two-stage harvesting as compared with the combined harvesting.

rapeseed; combined harvesting; two-stage harvesting; key traits of seed

10.3864/j.issn.0578-1752.2021.14.006

2020-09-08；

2020-09-21

國家重點(diǎn)研發(fā)計(jì)劃（2018YFD1000900）、湖北省技術(shù)創(chuàng)新專項(xiàng)重大項(xiàng)目（2017ABA064）

白晨陽，E-mail：baicy@webmail.hzau.edu.cn。通信作者蒯婕，E-mail：kuaijie@mail.hzau.edu.cn

（責(zé)任編輯 李莉）