劉美茹，張仙美，張權(quán)，鄭磊，吳鶴敏，戴晉，李天奇，劉立峰

（漯河市農(nóng)業(yè)科學(xué)院，河南 漯河 462000）

芝麻是我國重要的油料作物之一，種植面積逐年增加。我國芝麻種植廣泛，但分布極不平衡，以黃淮、江漢和長(zhǎng)江中下游為主產(chǎn)區(qū)，其中河南、湖北和安徽種植面積最廣［1，2］，芝麻產(chǎn)量占全國的75%左右［3］。2003年，河南、安徽、湖北等芝麻產(chǎn)區(qū)受長(zhǎng)期降雨影響而減產(chǎn)，其中河南駐馬店、漯河及安徽阜陽等大部分地區(qū)甚至絕收，導(dǎo)致當(dāng)年我國芝麻產(chǎn)量大幅下降［1］。因此，選育綜合性狀良好的芝麻品種對(duì)我國芝麻產(chǎn)業(yè)發(fā)展具有重要意義。

灰色關(guān)聯(lián)分析是基于灰色系統(tǒng)理論的一種重要的統(tǒng)計(jì)分析方法，此方法的最大特點(diǎn)就是將雜亂無章數(shù)據(jù)的隨機(jī)過程看作是灰色過程，將整理的標(biāo)準(zhǔn)化數(shù)據(jù)看作是由眾多因素共同影響和作用的結(jié)果，最后用較少的數(shù)據(jù)較直觀、準(zhǔn)確地反映問題［4，5］。目前該方法已廣泛應(yīng)用于多種農(nóng)作物產(chǎn)量性狀的關(guān)聯(lián)分析和品種的綜合評(píng)價(jià)以及區(qū)試多點(diǎn)、多品種結(jié)果分析等方面。楊志剛等［6］運(yùn)用灰色關(guān)聯(lián)分析法對(duì)小麥區(qū)試結(jié)果進(jìn)行分析及綜合評(píng)價(jià)，篩選出適宜本地種植的最佳小麥品種。蔡立夫等［7］用灰色關(guān)聯(lián)分析法對(duì)東北地區(qū)13個(gè)試點(diǎn)的14個(gè)品種進(jìn)行綜合評(píng)價(jià)，得出單株產(chǎn)量也影響花生產(chǎn)量的結(jié)果。柳青山等［8］對(duì)全國區(qū)試的4個(gè)糯高粱品種進(jìn)行綜合評(píng)價(jià)，得出灰色關(guān)聯(lián)度分析能夠充分表示各性狀對(duì)產(chǎn)量影響的大小。李愛國等［9］利用灰色系統(tǒng)理論對(duì)漯麥26的10個(gè)主要農(nóng)藝性狀與產(chǎn)量進(jìn)行灰色關(guān)聯(lián)分析，得出了小麥10個(gè)農(nóng)藝性狀與產(chǎn)量之間的關(guān)聯(lián)度。郭建秋等［10］采用灰色關(guān)聯(lián)分析方法，就15個(gè)大豆品種的主要性狀對(duì)產(chǎn)量的影響程度進(jìn)行綜合評(píng)價(jià)得出，灰色關(guān)聯(lián)分析能真實(shí)地反映大豆新品種的優(yōu)劣表現(xiàn)，并能夠充分利用產(chǎn)量以及產(chǎn)量相關(guān)的其它多個(gè)性狀的全部信息，使參試品種排序更加客觀合理。鄭磊等［11］利用相關(guān)性和灰色關(guān)聯(lián)分析方法，對(duì)芝麻單株產(chǎn)量及其構(gòu)成因素與相關(guān)性狀的關(guān)系進(jìn)行研究認(rèn)為，育種中要注重選育開花早、生育期長(zhǎng)、果軸長(zhǎng)、果節(jié)數(shù)多、單株蒴果多的品種。

綜上看出，灰色關(guān)聯(lián)分析可綜合多個(gè)性狀對(duì)品種進(jìn)行定量評(píng)價(jià)，而非單靠產(chǎn)量性狀評(píng)價(jià)品種優(yōu)劣，因此能夠更加全面、準(zhǔn)確地評(píng)價(jià)品種的優(yōu)劣［6］。目前用灰色關(guān)聯(lián)分析篩選適合某一地區(qū)最佳芝麻品種的研究較少，本研究采用該方法對(duì)多個(gè)芝麻品種進(jìn)行綜合評(píng)價(jià)，以期篩選出適合河南地區(qū)種植的品種。

1 材料與方法

1.1 試驗(yàn)材料

供試芝麻品種共20個(gè)（表1），其中皖芝11號(hào)、皖芝12號(hào)和皖芝5號(hào)來自安徽省農(nóng)業(yè)科學(xué)院作物研究所；項(xiàng)芝秸葉黃來自周口市農(nóng)業(yè)科學(xué)院；鄭太芝1號(hào)和鄭雜芝3號(hào)來自河南省農(nóng)業(yè)科學(xué)院；中芝14號(hào)、中芝23號(hào)和中芝30號(hào)來自中國農(nóng)業(yè)科學(xué)院油料作物研究所；周芝11號(hào)和周芝12號(hào)來自周口市農(nóng)業(yè)科學(xué)院；駐芝21號(hào)來自駐馬店農(nóng)業(yè)科學(xué)院，豫芝4號(hào)來自駐馬店地區(qū)農(nóng)科所；漯12、漯芝15號(hào)、漯芝16號(hào)、漯芝18號(hào)、漯芝19號(hào)、漯芝21號(hào)和漯芝23號(hào)來自漯河市農(nóng)業(yè)科學(xué)院。

1.2 試驗(yàn)方法

試驗(yàn)于2019年在漯河市農(nóng)業(yè)科學(xué)院試驗(yàn)基地進(jìn)行。前茬作物為小麥。土壤類型為黏土，肥力中上等。隨機(jī)區(qū)組排列，重復(fù)3次。小區(qū)凈面積20 m2。條播，10行區(qū)，行長(zhǎng)5 m，行距0.4 m，666.7m2種植密度為1.0萬株。每個(gè)重復(fù)間留有走道，試驗(yàn)地周圍設(shè)保護(hù)區(qū)。

成熟后每個(gè)小區(qū)隨機(jī)選取10株進(jìn)行田間考種，測(cè)定株高、始蒴高度、空梢尖長(zhǎng)度和單株蒴果數(shù)。收獲后室內(nèi)考種測(cè)定每蒴粒數(shù)、千粒重和單株產(chǎn)量。

1.3 單株產(chǎn)量及相關(guān)性狀的灰色關(guān)聯(lián)分析

按照灰色系統(tǒng)理論［12］，將20個(gè)芝麻品種作為一個(gè)灰色系統(tǒng)，設(shè)單株產(chǎn)量為參考數(shù)列X0，株高、始蒴高度、空梢長(zhǎng)度、果軸長(zhǎng)、單株蒴果數(shù)、每蒴粒數(shù)、千粒重分別為比較數(shù)列Xi，即（Xi＝X1，X2，…Xi，i＝7）。首先對(duì)原始數(shù)據(jù)進(jìn)行無量綱化處理，即采用均值化生成法對(duì)原始序列進(jìn)行無量綱化處理得到一個(gè)新數(shù)列，然后進(jìn)行標(biāo)準(zhǔn)化處理，根據(jù)標(biāo)準(zhǔn)化處理結(jié)果計(jì)算出X0與Xi差的絕對(duì)值，計(jì)算公式如下：

式中，△i（k）表示絕對(duì)差值，X0（k）表示參考數(shù)列值，Xi（k）表示比較數(shù)列標(biāo)準(zhǔn)化結(jié)果。

進(jìn)一步計(jì)算出單株產(chǎn)量與其它性狀的關(guān)聯(lián)系數(shù)，計(jì)算公式如下：

式中，ξi（k）表示關(guān)聯(lián)系數(shù)，△i（k）min和△i（k）max分別表示最小和最大差值，ρ表示分辨系數(shù)，取值0.5。

1.4 品種綜合評(píng)價(jià)的灰色關(guān)聯(lián)分析

研究品種優(yōu)劣時(shí)，需根據(jù)供試品種的各項(xiàng)性狀指標(biāo)，構(gòu)建一個(gè)各性狀表現(xiàn)優(yōu)良的理想品種作為參考品種［13，14］，參考品種標(biāo)記為Y0。參考1.3中的方法，首先對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，計(jì)算供試品種與參考品種各性狀的△j（k），并得到△j（k）min和△j（k）max，參考公式（2）計(jì)算供試品種與參考品種的關(guān)聯(lián)系數(shù)，并計(jì)算同一品種各性狀關(guān)聯(lián)系數(shù)的平均值得到等權(quán)關(guān)聯(lián)度，再根據(jù)各性狀的權(quán)重系數(shù)計(jì)算加權(quán)關(guān)聯(lián)度。

1.5 數(shù)據(jù)分析

用Microsoft Excel 2007和SPSS 19.0軟件進(jìn)行試驗(yàn)數(shù)據(jù)處理和灰色關(guān)聯(lián)分析。

2 結(jié)果與分析

2.1 不同芝麻品種主要性狀表現(xiàn)

由表1可知，20個(gè)芝麻品種株高為163.6～176.5 cm，平均為170.38 cm。始蒴高度為40.2～56.9 cm，平均為50.81 cm。空梢長(zhǎng)度為3.9～6.6 cm，平均為5.06 cm。果軸長(zhǎng)為106.6～131.3 cm，平均為114.5 cm。單株蒴果數(shù)為80.8～105.6個(gè)，平均為94.01個(gè)。每蒴粒數(shù)為61.8～70.9個(gè)，平均為66.77個(gè)。千粒重為2.69～3.20 g，平均為2.97 g。單株產(chǎn)量為10.3～17.2 g，平均為14.27 g。

表1 不同芝麻品種主要性狀表現(xiàn)

2.2 不同芝麻品種主要性狀的相關(guān)分析

由表2可知，果軸長(zhǎng)、單株蒴果數(shù)、千粒重、單株產(chǎn)量與株高、單株蒴果數(shù)與果軸長(zhǎng)呈極顯著正相關(guān)（P＜0.01）；每蒴粒數(shù)與單株蒴果數(shù)顯著正相關(guān)；千粒重與單株蒴果數(shù)極顯著正相關(guān)，與果軸長(zhǎng)、每蒴粒數(shù)顯著正相關(guān)（P＜0.05）；單株產(chǎn)量與果軸長(zhǎng)、單株蒴果數(shù)、每蒴粒數(shù)、千粒重極顯著正相關(guān)。由此看出，單株蒴果數(shù)、每蒴粒數(shù)、千粒重與芝麻單株產(chǎn)量有直接關(guān)系，株高、果軸長(zhǎng)對(duì)單株蒴果數(shù)有直接影響，從而間接影響單株產(chǎn)量。

表2 不同芝麻品種主要性狀間的相關(guān)系數(shù)

2.3 不同芝麻品種單株產(chǎn)量及其相關(guān)性狀的灰色關(guān)聯(lián)分析

對(duì)不同芝麻品種各性狀測(cè)定值進(jìn)行無量綱化處理，結(jié)果見表3。對(duì)表3數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化，結(jié)果見表4。根據(jù)表4數(shù)據(jù)計(jì)算每個(gè)性狀的絕對(duì)差值，結(jié)果見表5。根據(jù)表5結(jié)果可知，△i（k）min為0.0005，△i（k）max為0.9153，進(jìn)一步計(jì)算得出各性狀與單株產(chǎn)量的關(guān)聯(lián)系數(shù)，見表6。根據(jù)灰色關(guān)聯(lián)度分析原理，關(guān)聯(lián)度可代表系統(tǒng)中各性狀的重要性，其中關(guān)聯(lián)度越大，則表示該性狀的重要性越高，同時(shí)可表示與參考數(shù)列的關(guān)系越密切，反之關(guān)系則越遠(yuǎn)［15］。由表6可知，各性狀與單株產(chǎn)量的關(guān)聯(lián)系數(shù)排序?yàn)椋簡(jiǎn)沃贻艄麛?shù)、果軸長(zhǎng)、千粒重、每蒴粒數(shù)、株高、始蒴高度和空梢長(zhǎng)度。由此可知，單株蒴果數(shù)與芝麻單株產(chǎn)量關(guān)聯(lián)度最大，其次是果軸長(zhǎng)和千粒重。

表3 不同芝麻品種產(chǎn)量及其構(gòu)成因素?cái)?shù)據(jù)無量綱化處理

表4 不同芝麻品種單株產(chǎn)量及其構(gòu)成因素灰色關(guān)聯(lián)分析的標(biāo)準(zhǔn)化結(jié)果

表5 不同芝麻品種單株產(chǎn)量及其構(gòu)成因素灰色關(guān)聯(lián)分析的△i（k）

表6 不同芝麻品種單株產(chǎn)量與其相關(guān)性狀的關(guān)聯(lián)系數(shù)

2.4 不同芝麻品種綜合評(píng)價(jià)

參照楊志剛［6］和闞躍峰［14］等的方法選取參考品種即理想品種各性狀的平均表型值，然后對(duì)20個(gè)芝麻品種各性狀測(cè)定值進(jìn)行標(biāo)準(zhǔn)化處理并計(jì)算△j（k），見表7，可知，△j（k）min為0，△j（k）max為0.6923，進(jìn)而得到不同芝麻品種與參考品種的關(guān)聯(lián)系數(shù)，見表8。

表7 基于灰色關(guān)聯(lián)分析的不同芝麻品種綜合評(píng)價(jià)的△j（k）

表8 不同芝麻品種與參考品種的關(guān)聯(lián)系數(shù)

據(jù)表8中關(guān)聯(lián)系數(shù)進(jìn)一步計(jì)算可得到不同芝麻品種與理想品種的等權(quán)關(guān)聯(lián)度，等權(quán)關(guān)聯(lián)度越高，越接近理想品種，其綜合評(píng)價(jià)越高。由表9可知，鄭雜芝3號(hào)綜合評(píng)價(jià)最接近理想品種，等權(quán)關(guān)聯(lián)度為0.8789；其次是皖芝5號(hào)，等權(quán)關(guān)聯(lián)度為0.8334；第3位是漯芝16號(hào)，等權(quán)關(guān)聯(lián)度為0.8091。

實(shí)際生產(chǎn)中，各性狀指標(biāo)的重要性不同，因此要結(jié)合實(shí)際對(duì)不同性狀指標(biāo)賦予不同的權(quán)重系數(shù)，進(jìn)而計(jì)算加權(quán)關(guān)聯(lián)度。由表9可知，皖芝5號(hào)綜合評(píng)價(jià)最接近理想品種，加權(quán)關(guān)聯(lián)度最大，為0.9118；其次是鄭雜芝3號(hào)，加權(quán)關(guān)聯(lián)度為0.8711；第3位是項(xiàng)芝秸葉黃，加權(quán)關(guān)聯(lián)度為0.8415。

表9 不同芝麻品種與參考品種的關(guān)聯(lián)度及排序

對(duì)等權(quán)關(guān)聯(lián)度和加權(quán)關(guān)聯(lián)度進(jìn)行比較看出，二者最終結(jié)果相差不大，等權(quán)關(guān)聯(lián)度前5位分別為鄭雜芝3號(hào)、皖芝5號(hào)、漯芝16號(hào)、項(xiàng)芝秸葉黃和皖芝11號(hào)，加權(quán)關(guān)聯(lián)度前5位為皖芝5號(hào)、鄭雜芝3號(hào)、項(xiàng)芝秸葉黃、漯芝16號(hào)和駐芝21號(hào)。加權(quán)處理后，由于產(chǎn)量方面優(yōu)勢(shì)更加突出，因此產(chǎn)量較高的皖芝5號(hào)、項(xiàng)芝秸葉黃和駐芝21號(hào)位次提前，這一結(jié)果與孫建等［16］的研究結(jié)果相一致。

3 討論與結(jié)論

本研究就20個(gè)芝麻品種相關(guān)農(nóng)藝性狀對(duì)單株產(chǎn)量的影響進(jìn)行灰色關(guān)聯(lián)分析發(fā)現(xiàn)，單株蒴果數(shù)、果軸長(zhǎng)、千粒重、每蒴粒數(shù)、株高對(duì)產(chǎn)量影響較大，并且與單株產(chǎn)量呈極顯著正相關(guān)，單株產(chǎn)量與其相關(guān)性狀的關(guān)聯(lián)系數(shù)依次是單株蒴果數(shù)＞果軸長(zhǎng)＞千粒重＞每蒴粒數(shù)＞株高＞始蒴高度＞空梢長(zhǎng)度。因此，實(shí)際育種中應(yīng)把單株蒴果數(shù)、果軸長(zhǎng)、千粒重、每蒴粒數(shù)等性狀作為重要性狀進(jìn)行選擇。本研究采用灰色關(guān)聯(lián)分析法對(duì)20個(gè)芝麻品種進(jìn)行綜合評(píng)價(jià)得出，皖芝5號(hào)、鄭雜芝3號(hào)、漯芝16號(hào)和項(xiàng)芝秸葉黃綜合評(píng)價(jià)較好，更適宜河南地區(qū)芝麻大規(guī)模種植。

灰色關(guān)聯(lián)分析法已經(jīng)應(yīng)用在農(nóng)業(yè)科研的多個(gè)領(lǐng)域［17－19］，是一個(gè)以發(fā)展變化系統(tǒng)進(jìn)行發(fā)展動(dòng)態(tài)量化比較的分析方法，不同時(shí)間、地點(diǎn)、品種及氣候條件均能影響產(chǎn)量性狀，因此育種中應(yīng)用該方法時(shí)應(yīng)對(duì)不同環(huán)境、不同品種做具體分析，才能選育出適宜當(dāng)?shù)胤N植的優(yōu)良品種。