甘蔗(Soccharum officinarum)是重要的制糖原料，原产于热带和亚热带地区。广西是我国重要的产糖基地，甘蔗种植面积和产量均位居全国之首^[1]。近年来，由于甘蔗生产成本不断攀升、机械化程度低^[2]、优良品种缺乏等因素的影响，广西蔗糖产业发展受到一定程度的冲击。选育和推广优良甘蔗品种是保障蔗糖产业健康持续发展的有效途径^[3]，而科学合理的甘蔗综合性状评价方法和结果可为甘蔗品种的选育和推广提供重要参考^[4]。

灰色关联度分析可衡量系统内各因素之间的关联程度，对研究对象的贡献度作出评价^[5]。TOPSIS法作为一种解决多属性决策问题的方法，已应用于多种作物的综合评价^[6]。目前已有部分研究采用灰色关联度对甘蔗性状作评价，而TOPSIS法则较少运用在甘蔗综合性状的评价中，这可能会导致一些具备优良综合性状的甘蔗品种难以得到应用。因此，综合运用灰色关联度及TOPSIS法对甘蔗综合性状进行客观评价，对选育和推广高产、高糖、抗逆性强的甘蔗优良新品种具有重要意义。桃联安等^[7]采用灰色关联度对云南82-114割手密F1代的9个性状进行分析后发现，蔗茎产量对产糖量的影响最大，提高甘蔗产量能更有效地提高单位面积产糖量。赵俊等^[8]对19个国外引进的甘蔗品种的9个主要性状进行灰色关联度分析后也发现，公顷蔗茎产量与公顷含糖量之间的关联度最大。黄婷等^[9]对玉米(Zea mays)自交系各性状进行灰色关联度分析，结果发现秃尖长性状与产量构成性状关联度最密切。李小玉等^[10]对12份油菜(Brassica napus)品种在不同处理条件下的多个指标进行灰色关联度分析，得出干旱胁迫下与单株产量关联度最密切的是产量因子的结论。姜敏等^[11]对18份贵州地区花生(Arachis hypogaea)种质资源的10个主要农艺性状与花生产量的相关性进行灰色关联度分析发现，与产量关联度最密切的指标为成苗率，其次为单株生产力，最后为结果数，同时指出为提高产量，在花生栽培过程中要注重对成苗率、单株生产力等指标的选择。曹燕燕等^[12]以4个小麦(Triticum aestivum)品种为材料，测定其植株在受到低温胁迫影响下，叶片中超氧化物歧化酶(SOD)等生理指标的活性变化并测定其单株产量等指标，经灰色关联度分析发现，脯氨酸和可溶性蛋白含量与小麦产量降幅关联度最大。单晓敏等^[13]运用TOPSIS法与其他方法进行联合分析，建立并验证了甘蔗联合收割机液压系统设计方案的有效性及合理性。梁秋萍等^[14]通过对10个不同品种甜樱桃(Prunus avium)理化性质及总酚物质进行测定并构建TOPSIS模型，筛选出抗氧化能力最强且综合营养价值最高的樱桃优良品种。包淑英等^[15]将TOPSIS法用于13个绿豆(Vigna radiata)品种(系)的综合评价，获得4个可示范推广的优良品种。雷涌涛等^[16]将11个大麦(Hordeum vulgare)新品系的产量、基本苗等10个性状作为正向指标，以白粉病等3个性状作为负向指标，进行TOPSIS评价，筛选出保大麦10-J5等4个适合在云南推广种植的大麦新品系。安东等^[17]将33个原产于南部地区的大豆(Glycine max)品种引种到黄土高原地区，对其粗蛋白等多个指标进行TOPSIS分析，发现华夏3号等2个大豆品种产量和品质表现优良，适宜作为优质饲草大豆品种推广应用。目前，对甘蔗农艺性状的评价大多采用方差分析方法，运用灰色关联度和TOPSIS方法对甘蔗农艺性状进行评价的研究鲜见报道。本研究以桂热1号(GR1)和桂热2号(GR2)甘蔗品种(系)为试验材料，以新台糖22号(ROC22)为对照，进行6次新植4次宿根(6新4宿)种植试验，对其产糖量等7个主要农艺性状进行灰色关联度分析，筛选出与产糖量关联度最大的性状因子，同时将黑穗病等3个病虫害发病率作为负向指标、产糖量等8个农艺性状作为正向指标，运用TOPSIS法进行综合分析，以筛选出综合性状最优的甘蔗品种，为甘蔗新品种的选育和推广提供参考依据。

1 材料与方法 1.1 材料

参试甘蔗品种(系)包括GR1和GR2，以ROC22为对照(CK)。其中，甘蔗品系GR1原编号为LZ03-75，其母本为粤糖91-976^[18]，父本为ROC10^[19]；甘蔗品种GR2原编号为LZ02-6，其母本同为粤糖91-976，而父本为ROC20^[20]。GR1和GR2是由其父、母本分别通过传统杂交育种方式选育，并经过单株筛选、株系选择、初级品比试验及在攀西蔗区进行区域性试验筛选出的甘蔗优良新品种(系)，由广西壮族自治区亚热带作物研究所(以下简称“广西热作所”)与四川省凉山州亚热带作物研究所合作，于2012年引入广西蔗区种植^[21]。

1.2 方法 1.2.1 试验设计

在2012-2017年，共进行6新4宿的小区品比试验。试验地分别设在广西热作所第二试验区、金光农场、新兴农场、良圻农场和滨海农场。各试验品种(系)种植小区设3个重复，随机排列，每小区行长7.0 m，共设5行，行距为1.2 m，小区面积约为42.0 m²，边际设保护行。田间栽培管理遵循当地蔗田栽培管理模式，同时按时对参试品种(系)进行数据调查和收集。

1.2.2 调查项目及方法

下种后30 d开始调查统计各品种(系)的出苗率，且每隔15 d再调查1次，共调查3次；自蔗苗分蘖开始，每隔15 d调查各小区各品种(系)分蘖率，共调查3次；于5月中旬调查各品种(系)甘蔗的枯心苗数、黑穗病苗数和梢腐病苗数；新植蔗收获后的次年调查各小区各品种(系)宿根蔗的分蘖率；第二年开始于每年1月测定各小区各品种(系)的株高、茎径、有效茎数、蔗茎产量、蔗糖分和产糖量^[3]。

病虫害率(%)=病虫害植株数/植株总数×100%，

有效茎数(株/hm²)=小区有效茎数/小区面积×10⁴，

蔗茎产量(t/hm²)=单茎重×公顷有效茎数，

单茎重(kg)=(茎径)²×茎长×0.7854×1/1000(其中1表示茎比重)，

蔗糖分(%)=(平均田间锤度×1.0825-7.703)×100%，

产糖量(t/hm²)=蔗糖分×蔗茎产量，

式中，平均田间锤度为随机测量各小区各品种(系)中间行15-20株甘蔗田间锤度的平均值。

1.2.3 分析方法

根据邓聚龙^[22]的灰色系统理论，将3个甘蔗品种(系)6新4宿的6个主要农艺性状调查结果作为子序列，以产糖量作为母序列，进行灰色关联度分析，再对3个病虫害发病率(负向指标)及产糖量等8个农艺性状(正项指标)进行TOPSIS分析。灰色关联度有关指标的计算方法如下。

$ \text { 数据的标准化处理: } x_i(k)^{\prime}=\left[x_i(k)-\bar{x}_i\right] / s \text {, } $

$ \text {绝对差值: } \Delta_{0 i}(k)=\left|x_0(k)-x_i(k)\right|, 1 \leqslant$ $i \leqslant m, $

$ \text { 关联系数: } L_{0 i}(k)=\frac{\Delta_{\min }+\Delta_{\max }}{\Delta_{0 i}(k)+\rho \Delta_{\max }} \text {, } $

式中，x_i表示比较数列，x₀表示参考序列，x_i为该性状的平均值，s为该性状的标准差；Δ_max和Δ_min分别表示比较序列绝对差值中的最大值和最小值，ρ为分辨系数，取值为0.1。

TOPSIS法分析步骤为

① 数据同趋势化处理：将调查得到的黑穗病病害率、枯心苗率和梢腐病率进行逆向化(取倒数)处理。

② 数据的平方和归一化处理：

$ A_{i j}=X_{i j} / \sqrt{\sum_{i=1}^m x_{i j}^2} \text { 。} $

③ 找出最优矩阵向量A⁺和最劣矩阵向量A^－：

$ \begin{aligned} & A^{+}=\left(A_1^{+}, A_2^{+}, A_3^{+}, \cdots, A_n^{+}\right), A^{-}=\left(A_1^{-}, A_2^{-}, A_3^{-}, \cdots, \right. \\ & \left.A_n^{-}\right) 。\end{aligned} $

④ 分别计算各个评价对象与正、负理想解的距离D⁺和D^－：

$ D^{+}=\sqrt{\sum_{j=1}^m W_j \times\left(A_j^{+}-A_{i j}\right)^2}, $

$ D^{-}=\sqrt{\sum_{j=1}^m W_j \times\left(A_j^{-}-A_{i j}\right)^2} 。$

⑤ 结合距离值计算得出评价对象与最优方案的接近程度C值并进行排序：

$ C=D^{-} /\left(D^{+}+D^{-}\right) 。$

以上式子中，假定为i个评价对象、m个评价指标，A_ij表示规范化决策矩阵中的元素，X_ij表示矩阵中的元素，W_j表示权重值。

1.3 统计分析

试验数据采用Excel 2003进行整理，采用SPSS 18.0进行方差分析和TOPSIS分析，采用DPS V9.01进行灰色关联度分析。

2 结果与分析 2.1 各甘蔗品种(系)抗病虫害能力比较

由表 1可知，各甘蔗品种(系) 6新4宿的病虫害发生率中，CK的平均枯心苗率、平均黑穗病率和平均梢腐病率均最高，而GR2均为最低。值得注意的是，2012-2013年广西热作所第二试验区3个品种(系)新植及宿根蔗的枯心苗、黑穗病和梢腐病发病率均为0.00%，可能与试验区的土壤环境及当时的温度气候等有一定的关系。

2.2 各甘蔗品种(系)的主要农艺性状表现

由表 2可知，3个品种(系)的新植蔗平均出苗率、平均分蘖率、平均株高及平均茎径之间没有显著差异，平均出苗率最高的为GR2，比CK高6.14%；平均分蘖率最高的为GR1，比CK高11.66%，GR2的平均分蘖率比CK高8.23%；平均株高最高的为CK，比GR2高0.41 cm；平均茎径最大的为GR1，其次为GR2；GR2的平均有效茎数显著高于GR1和CK，较CK的增幅达22.00%。2016-2017年，GR2在新兴农场和滨海农场的宿根蔗的有效茎数均显著高于CK；GR1在新兴农场的新植蔗的有效茎数显著高于CK。综上所述，在多年多点的跟踪种植试验中，GR2在有效茎数方面的表现优势明显。

2.3 各甘蔗品种(系)的蔗茎产量和产糖量比较

表 3显示，GR2的平均蔗茎产量最高，较CK增产25.04 t/hm²，增幅达29.14%；GR1的平均蔗茎产量次之，较CK增产8.04 t/hm²，增幅为9.36%。GR2的平均蔗茎产量显著高于GR1和CK。平均蔗糖分最高的为CK，其次为GR2；GR2和CK的平均蔗糖分均超过13%。平均产糖量最高的为GR2，较CK增产3.35 t/hm²，增幅达28.03%；GR1的平均产糖量较CK增产0.67 t/hm²，增幅为5.61%。GR2新植蔗在2012-2013年广西热作所第二试验区、2013年及2016-2017年新兴农场蔗茎产量显著高于CK，在2013年及2016-2017年新兴农场产糖量显著高于CK；相较于新植蔗，GR2宿根蔗的表现则更突出，其蔗茎产量均显著高于CK，蔗糖分平均值与CK差异不显著，而在宿根蔗产糖量方面，除2012-2013年广西热作所第二试验区表现与CK差异不显著外，其余年份均显著高于CK，总体上GR2宿根蔗的蔗茎产量和产糖量均表现出显著优势。

2.4 产糖量与主要农艺性状的关联度分析

根据灰色关联度分析法，以产糖量为母序列，蔗茎产量、有效茎数、株高、分蘖率、蔗糖分和茎径等主要农艺性状为子序列，计算序列之间的关联系数。由表 4可知，与产糖量关系最密切的是蔗茎产量，其次为有效茎数，再次为株高；而茎径与产糖量的关联系数最小。这说明蔗茎产量对产糖量的影响最大，其次为有效茎数、株高和分蘖率，茎径对产糖量影响最小。

2.5 不同品种(系)甘蔗农艺性状的综合评价

根据TOPSIS法首先计算3个品种(系)甘蔗6新4宿11个指标的正、负理想解(A⁺和A^－)，结果如表 5所示，正理想解A⁺表示各评价指标的最大值，而负理想解A^－表示各评价指标的最小值；然后计算出各评价对象分别与正、负理想解的距离值D⁺和D^－。与正理想解的几何距离最近、负理想解的几何距离最远的为最佳方案(品种)，结果如表 6所示。经TOPSIS法计算得出，3个品种(系)甘蔗农艺性状的综合评价得分排序为GR2＞GR1＞CK。

3 讨论

选育高产、高糖和抗逆性强的甘蔗品种是提高甘蔗产量和蔗糖分的有效途径^[3]。目前，我国经过长期多单位的大量选育种已获得较多甘蔗新品种，但对这些新品种大多采用方差分析方法进行某一性状的比较^{[23, 24]}。甘蔗品种(系)GR1和GR2由四川省攀西蔗区引入广西，目前已有近10年的种植时间，期间朱鹏锦等^[25]通过低温胁迫处理耐寒性不同的甘蔗品种发现，GR2等品种具有较强的低温耐受能力；欧克纬等^[21]通过对10个甘蔗品种(系)在柳州地区种植后的各性状进行比较后发现，GR2在产量和产糖量方面均表现出优势；单彬等^[26]使用5种除草剂对GR2做安全性测试后，得出GR2对敌草隆、莠灭净等常见除草剂具有较好的耐受性的结论。上述试验多采用方差分析对甘蔗的某单一性状作评价，这可能使得甘蔗品种(系)的综合性状难以得到客观全面的体现。为评估GR2等甘蔗品种(系)在广西蔗区的适应性及产量潜力，本研究特对其进行跨地点、跨年份的跟踪种植调查并采用灰色关联度和TOPSIS法对其综合性状进行评价。

灰色关联度在分析作物多种性状方面的评价效果较好，其评价结果有利于在开展作物品种选育时把握好方向，已应用于包括甘蔗在内的多种作物产量与相关农艺性状的关系研究^[27-30]。王尊欣^[31]对31个甘蔗品种的公顷蔗产量与其他10个性状进行灰色关联度分析后发现，公顷含糖量与公顷蔗茎产量之间的关联度最大。刘晓燕等^[32]将10个引进的甘蔗新品种的公顷含糖量与其他7个性状进行灰色关联度分析后发现，与公顷含糖量关联度最大的是公顷实产。TOPSIS法能对作物品种多个评价指标与理想化目标的接近程度作出分析，进而以此作为依据评估作物品种的优劣，是一种较全面且客观的评价方法^[33]。宁东贤等^[34]利用TOPSIS法对12个花生品种(系)作出综合评价，获得冀花13号等4个可推广种植的品种。上述研究结果均可作为选育作物品种的侧重点参考，以及评估作物品种综合性状和推广价值的参考依据。

本研究中，3个甘蔗品种(系) 6新4宿试验结果表明，GR2的平均产糖量达15.30 t/hm²，GR1达12.62 t/hm²，分别较CK增产3.35 t/hm²和0.67 t/hm²，增幅分别达28.03%和5.61%。GR2在以各试验地为代表的桂中地区、桂北地区以及桂南地区的种植中均表现出产量及产糖量优势，这与其他同类研究的结果一致，如在来宾蔗区GR2的产量和产糖量均为最高^[35]，在柳州蔗区GR2宿根蔗的产量和产糖量方面也表现出较大优势^[21]；而GR1则在桂南地区的种植中表现出优势。灰色关联度分析与TOPSIS法评价的结果较一致。TOPSIS法评价结果表明，3个品种综合性状的评价排序为GR2＞GR1＞ROC22。灰色关联度分析结果表明，与产糖量关联度最高的是蔗茎产量，这与桃联安等^[7]及赵俊等^[8]的研究结果一致。综合评价排名第一的GR2的蔗茎产量平均值最高，达110.96 t/hm²，其有效茎数平均值也显著高于GR1和CK，达77 958株/hm²。此外，从11月至次年1月的甘蔗蔗糖分变化趋势可以看出，GR2属于晚熟品种^[36]。1月之后，GR2在一些地区的平均蔗糖分在14%以上，如2013年新兴农场GR2新植蔗1月的蔗糖分在15%以上。广西蔗区的榨季为11月中下旬一直持续到次年3月，十分有利于GR2等晚熟甘蔗品种的推广。

4 结论

甘蔗新品种GR2具有蔗茎产量高、有效茎数多、宿根性强等特点，经灰色关联度和TOPSIS法综合评价，其综合性状优于GR1和ROC22，是优良的晚熟甘蔗品种，适宜在广西蔗区进一步推广种植。