我国是饲料生产、消费大国，随着畜牧业发展，大豆、玉米等粮食的进口量逐年增加[1]。青稞在西藏粮食作物种植中具有重要地位，2020年西藏粮食作物的播种面积为183 300公顷，青稞的种植比例达78.2%，种植面积为143 400公顷[2]。大麦是一种优良的饲料，蛋白质、氨基酸含量高于玉米，在西藏的需求量较高[3]。在青稞灌浆期至乳熟期刈割作为青贮饲料或晒干作饲草，具有营养价值高、适口性好的特点[4]。研究表明，不同青稞秸秆的留茬高度和刈割时期其营养品质存在差异[5-6]。关于青稞饲用品质研究多集中在粗蛋白质（CP）、粗纤维（CF）、干物质、钙（Ca）、磷（P）等方面[7-9]，原因是不同品种（系）的青稞秸秆营养成分不同，造成青稞品种（系）的秸秆营养饲用价值存在差异[10]。营养成分是饲用价值的关键因素[11]，有必要对青稞秸秆的营养成分进行系统研究，明确不同青稞品种秸秆的营养特性，提高秸秆的饲用比例有利于缓解西藏牲畜饲草料短缺问题，促进青稞秸秆的资源化利用，延长产业链，增加附加值。西藏天然草场过载、退化严重，发展农区畜牧业能够缓解草场压力，饲草又是限制农区畜牧业发展的瓶颈，因此从青稞秸秆饲料化角度进行研究，研制秸秆饲料对畜牧业发展有利[12]。本研究选择7个不同青稞品种（系），对试验材料的CF、Ca及P等营养成分进行检测，利用数据分析相关软件进行主成分及相关性分析，以期找出营养特性较好的青稞品种，为青稞的饲料化种植提供参考。1材料与方法1.1试验地概况试验在西藏自治区农牧科学院农业研究所5号试验地进行，前茬作物为小麦，土质为砂壤土，播种前精耕细作，每667 m2人工撒施基肥磷酸二铵7.5 kg、尿素5 kg，于2021年4月20日人工小锄头开沟条播。1.2试验材料试验材料分别为区内主推新品种（系）14-3492、QTB13、QTB25、喜拉22、13-6927、藏青3000、藏青2000，为西藏自治区农牧科学院农业研究所自行培育的高生物量、强分蘖的新品种（系）。1.3试验设计采用随机区组试验设计，重复3次，小区面积为20.01 m2，行长为6.67 m，行宽为3 m，每小区12行，行距为0.25 m，四周设有保护行，全生育期进行田间统一管理。在6月28日至7月19日，共取样5次，每隔7 d随机采集每个品种（系）的秸秆，采用四分法取样100 g。1.4测定指标及方法秸秆取样后，在60 ℃烘箱中处理24 h，称重，烘2 h称重，两次重量差无明显差异，此时样品已非常干燥，易于粉碎，使用研钵粉碎，过40目筛，制成样品。CP采用GB 5009.5—2016中凯氏定氮法测定；EE采用GB 5009.6—2016索氏抽提法测定；Ash的测定方法参考GB/T 6438—1992；CF的测定方法参考 GB/T 6434—1994；NDF按GB/T 20806—2006测定；ADF按NY/T 1459—2007测定；Ca按照国标GB 5009.92—2016测定；P含量测定采用钼锑抗比色法按国标NY/T 2421—2013测定。1.5数据统计与分析数据采用Excel 2019及SPSS 19软件进行相关性分析、方差分析和主成分分析。根据公式（1）对数据进行转换，采用SPSS软件进行主成分分析[13]，获取每个主成分分值（Fi），样品主成分的分值Fi乘以其对应的贡献率（Ei），相加可得综合得分（D），见公式（2）。Sin=(Xin-Ximax)/(Ximax-Ximin) (1)D=∑Fi×Eii=1,2,3,4……  (2)式中：Sin为隶属函数值，指第n个样品的i项指标经过对原始数据转化后所得；Xin为样品原始测定数据；Ximax、Ximin为样品中第i个指标的最大值和最小值。2结果与分析2.1青稞秸秆的营养品质分析（见表1）由表1可知，QTB13的CP含量显著高于13-6927、喜拉22、藏青2000（P0.05），秸秆的CP平均含量为10.39%，变异系数0.17%；QTB13的EE含量显著高于藏青3000、藏青2000、14-3492、喜拉22（P0.05），EE的平均含量为1.92%，变异系数0.28%；14-3492的Ash含量显著高于藏青3000、QTB25、QTB13（P0.05），Ash的平均含量为10.22%，变异系数0.24%；喜拉22的CF含量显著高于藏青3000、QTB25、QTB13（P0.05），秸秆的CF平均含量为34.02%，变异系数0.09%；藏青2000的NDF含量显著高于除13-6927外的秸秆（P0.05），NDF的平均含量59.32%，变异系数0.08%；藏青2000的ADF含量显著高于藏青3000、QTB25、QTB13、喜拉22（P0.05），秸秆的ADF平均含量36.78%，变异系数0.11%；藏青3000的Ca显著高于藏青2000、QTB25、13-6927、QTB13、喜拉22（P0.05），秸秆的Ca平均含量2.10 g/kg，变异系数0.76%；喜拉22的P含量显著高于QTB13、藏青3000、14-3492、藏青2000（P0.05），秸秆的P平均含量2.09 g/kg，变异系数0.23%。10.13557/j.cnki.issn1002-2813.2023.06.017.T001表1青稞秸秆的营养品质分析项目CP/%EE/%Ash/%CF/%NDF/%ADF/%Ca/(g/kg)P/(g/kg)QTB1311.69±2.37a2.23±0.33a8.88±0.69c28.88±0.69c55.60±2.18c32.75±1.54c1.32±0.39d2.05±0.36bQTB2511.68±3.10a2.15±0.10a8.90±1.19c32.48±1.51b59.05±1.09b35.80±1.25b2.07±0.43bc2.10±0.26ab喜拉229.39±5.18bc1.62±0.16c10.86±3.46ab37.30±4.05a56.26±6.11c31.62±3.84c1.07±0.32d2.34±0.19a13-692710.47±1.09b2.00±0.30a11.38±3.48a34.52±3.05ab62.78±2.80a38.10±1.92ab1.70±0.60c2.33±0.76a藏青300010.68±2.68ab1.90±0.35b9.56±2.84bc33.80±1.52b56.54±6.40c36.78±1.92b3.10±3.59a2.04±0.53b14-349211.65±1.19a1.75±0.33c11.50±0.78a34.45±1.75ab59.88±2.68b38.95±1.79ab2.94±0.76ab2.04±0.33b藏青20008.35±2.76c1.85±0.30b10.20±2.16ab35.25±2.05a63.58±3.79a41.85±3.21a2.45±1.27b1.81±0.57c平均值10.391.9210.2234.0259.3236.782.102.09变异系数/%0.170.280.240.090.080.110.760.23注：同列数据肩标字母不同表示差异显著（P0.05），相同字母或无字母表示差异不显著（P0.05）。2.2青稞秸秆的营养品质指标相关性分析结果（见表2）由表2可知，秸秆中CP含量与EE、Ash、Ca含量呈正相关，与秸秆中P含量呈显著正相关（P0.05），而与秸秆中NDF、CF、ADF含量呈负相关；秸秆中EE含量与Ca、P含量呈正相关，与秸秆中Ash含量呈显著正相关（P0.05），与CF含量呈极显著负相关（P0.01）；秸秆中Ash含量与Ca含量呈正相关，与P含量呈显著正相关（P0.05），与秸秆中CF、NDF、ADF含量呈负相关；秸秆中CF含量与NDF、ADF含量呈正相关，与秸秆中Ca、P含量呈负相关；秸秆中NDF与ADF含量呈极显著正相关（P0.01），与Ca含量呈显著负相关（P0.05），与P含量呈负相关；秸秆中ADF含量与Ca及P的含量呈显著负相关（P0.05），Ca与P含量呈显著正相关（P0.05）。10.13557/j.cnki.issn1002-2813.2023.06.017.T002表2秸秆相关性分析项目CPEEAshCFNDFADFCaPCP1.000EE0.2081.000Ash0.2380.348*1.000CF-0.320-0.459**-0.1611.000NDF-0.105-0.075-0.0950.1001.000ADF-0.218-0.063-0.0670.0500.722**1.000Ca0.0960.0220.134-0.089-0.420*-0.355*1.000P0.345*0.2100.406*-0.062-0.242-0.408*0.296*1.000注：“*”表示相关性显著（P0.05）；“**”表示相关性极显著（P0.05）。2.3青稞秸秆的营养品质性状主成分分析结果（见表3~表5）由表3可知，对于存在相关性的一组指标，在进行主成分分析时需要将各指标的相关性消除，删去重叠信息，降低信息维度，减少信息干扰。因此为发现青稞秸秆品质的主要影响指标，需要进行以下分析：KMO值为0.8360.500，表示变量之间存在相关性。10.13557/j.cnki.issn1002-2813.2023.06.017.T003表3KMO和巴特利特检验KMO取样适切性量数巴特利特球形度检验近似卡方自由度显著性0.836102.881360.000由表4可知，每个载荷量表示主成分与对应变量间的相关系数，用主成分载荷矩阵中的数据除以主成分相对应的特征值再开平方根，计算出主成分中每一个指标所对应的系数及特征向量A1、A2、A3。10.13557/j.cnki.issn1002-2813.2023.06.017.T004表4因子载荷矩阵结果项目主成分特征向量123A1A2A3CP-0.418-0.3310.6820.2520.2080.057EE0.201-0.7900.1550.074-0.3210.429Ash-0.2950.5670.6000.068-0.410-0.088CF0.0060.763-0.3210.1180.0830.496NDF0.7450.2480.438-0.0580.436-0.056ADF0.8060.2180.3040.352-0.0850.111Ca0.633-0.0210.1820.360-0.0310.035P-0.7100.2320.3570.242-0.106-0.044注：提取方法为主成分分析法。由表5可知，根据特征值1的原则提取了3个主成分，特征值分别为2.411、1.802、1.398，其累计方差贡献率达到70.127%，其中主成分1、主成分2、主成分3的方差贡献率为30.134%、22.523%、17.470%。10.13557/j.cnki.issn1002-2813.2023.06.017.T005表5主成分提取分析主成分初始特征值特征值方差百分比/%累计贡献率/%12.41130.13430.13421.80222.52352.65731.39817.47070.12740.7639.54079.66650.5586.97686.64260.4725.90492.54670.3664.57097.11680.2312.884100.0002.4青稞秸秆营养品质综合评价结果（见表6）根据公式F1=-0.173X1+0.083X2-0.123X3+0.003X4+0.309X5+0.334X6+0.262X7-0.294X8，F2=-0.184X1-0.439X2+0.315X3+0.424X4+0.138X5+0.121X6-0.012X7+0.129X8，F3=0.488X1+0.111X2+0.429X3-0.23X4+0.314X5+0.218X6+0.13X7+0.256X8建立品质评价的数学模型F总=（30.134F1+22.523F2+17.470F3）/100。由表6可知，各秸秆的营养成分排序为14-3492QTB2513-6927藏青3000藏青2000QTB13喜拉22。10.13557/j.cnki.issn1002-2813.2023.06.017.T006表6秸秆营养品质综合评价结果项目得分F总排序F1F2F3喜拉2225.0223.7928.4117.867藏青300027.2926.2129.4319.274QTB1324.8328.6726.1418.516QTB2528.9528.9731.1120.682藏青200027.2127.0228.4419.25513-692728.7628.6931.2420.58314-349231.6029.9530.1321.5312.5青稞秸秆营养品质聚类分析（见图1）由图1可知，青稞秸秆的营养品质可分为三大类，14-3492、QTB25、13-6927聚为一类，藏青3000、藏青2000、QTB13聚为一类，喜拉22单独为一类。根据试验材料秸秆的综合评价及聚类分析可知，秸秆14-3492、QTB25、13-6927的营养品质相似且饲用品质相对较高。10.13557/j.cnki.issn1002-2813.2023.06.017.F001图1秸秆营养品质的聚类分析3讨论3.1青稞秸秆的饲用品质分析西藏年平均日照时数为2 900~3 500 h，平均降雨量220~650 mm，无霜期在125 d以上，全年平均气温在3~17 ℃，天然草场产草量较低，草地退化严重。对于西藏农区畜牧业发展，一些专家较提倡进行饲草种植[14]。青稞在西藏具有良好的适应性，青稞秸秆可作为饲草使用，但需要明确秸秆的营养成分以选择饲用的青稞品种。秸秆营养成分的种类及含量可评价饲草营养价值[15]，饲用品质经常由ADF、NDF指标进行评价[16]。Ash过高表明饲料原料中杂质含量过多，CF是产生热能的物质，能够提高秸秆饲料对牲畜的适口性[17]，CP是评价饲料营养的重要指标[18]，是牲畜生长发育过程中蛋白质的主要来源。动物对EE需求量不大，但也不能缺少。CP和CF含量对饲草品质具有决定性作用，一般认为CP含量较高、CF含量较低是优质牧草，饲草转化率也会相对较高，具有较好的饲用价值[19]。3.2青稞秸秆的饲用品质综合评价分析对饲用品质系统进行描述时，一般会尽可能选取多的指标，从不同角度对系统进行描述，以获得尽可能多的信息。但是这些信息存在重叠，变量彼此间会有联系，很难掌握主要信息，使用SPSS软件进行主成分分析，能够从许多信息中提取信息特征，进而筛选信息，减少其他信息对有用信息的影响。主成分分析在对牧草营养品质评价使用较多，能够较好地发现营养成分综合评价较好的饲草品种（系）。本试验经过对样品秸秆进行主成分分析发现，各秸秆的营养成分排序为14-3492QTB2513-6927藏青3000藏青2000QTB13喜拉22，在种植饲用青稞时可优先选择这些品种，对饲草的转化率及畜牧业的发展具有促进作用。青稞品种（系）秸秆CP含量变幅8.35%~11.69%，CF含量变幅28.88%~37.30%，玉米CP含量是8.15%[20]，青稞秸秆CP含量高于玉米全株，表明青稞秸秆具有良好的饲用价值，是比较优质的饲用牧草。西藏青稞秸秆的数量约186.49万t。青藏高原整体高寒、干旱、生态脆弱、不利于饲草生长，天然草场面临超载、退化的风险，若能够充分利用青稞秸秆，可有效补充牲畜饲料蛋白来源，缓解高原饲草料不足的状况。4结论本试验发现，14-3492、QTB25、13-6927青稞秸秆的饲用品质相对较高，综合考虑秸秆营养品质的情况下区内种植饲用青稞时可优先选择上述品种（系）。