基于近红外光谱净信号的烟叶加料均匀性检测.pdf

2023 年 9 月第 56 卷第 9 期烟草科技Tobacco Science&TechnologySep.2023Vol.56 No.9摘要：利用近红外光谱技术与净信号分析方法检测烟叶加料的均匀性。将加料后样品的近红外光谱正交投影至加料前样品光谱构建的子空间，得到被测成分料液的净信号，通过计算净信号的范数表征加料后样品间的相似程度。结果表明，净信号范数的大小与料液的施加浓度（质量分数）正相关。通过测量生产线上加料后样品的近红外光谱，计算出各样品的净信号范数进而计算得到的加料均匀系数与行业标准方法所测结果具有相同趋势。基于近红外光谱的净信号分析方法可快速评价样品的加料均匀性，操作简便，检测效率高。运用净信号分析方法分析不同间套和阀芯位置下样品的加料均匀性，能够实现加料机喷嘴的间套和阀芯的位置优化。关键词：烟叶；加料；均匀性；近红外光谱；净信号中图分类号：TS452文献标志码：A文章编号：1002-0861（2023）09-0078-06收稿日期：2022-11-02录用日期：2023-07-20基金项目：贵州省科技计划支持项目“复杂的生物源料Primatone和Yeastolate的分析方法研究”（黔科合基础 2018 1130）；国家烟草专卖局科技重大专项“烟草近红外大数据构建及应用研究”110201901023（SJ-02）；贵州中烟工业有限责任公司科技项目“卷烟加料均匀性表征方法的研究与应用”（GZZYKJ/JZ2019GSY0110）。第一作者：胡芸（1978），女，博士，工程师，从事近红外光谱数据分析方法研究。E-mail：*通信作者：李博岩（1979），男，博士，副教授，从事化学计量学算法研究。E-mail：boyan_引用本文：胡芸，李博岩，庞红蕊，等.基于近红外光谱净信号的烟叶加料均匀性检测 J.烟草科技，2023，56（9）：78-83，90.（HUYun，LI Boyan，PANG Hongrui，et al.Detection of tobacco casing uniformity based on net analyte signal of near-infrared spectra J.Tobacco Science&Technology，2023，56（9）：78-83，90.DOI：10.16135/j.issn1002-0861.2022.0805）基于近红外光谱净信号的烟叶加料均匀性检测胡芸1，李博岩*2，庞红蕊1，张辞海1，刘剑11.贵州中烟工业有限责任公司技术中心，贵阳市开发大道96号5500092.贵州医科大学公共卫生与健康学院，贵阳市花溪大学城550025Detection of tobacco casing uniformity based on net analyte signalof near-infrared spectraHU Yun1,LI Boyan*2,PANG Hongrui1,ZHANG Cihai1,LIU Jian11.Technology Center,China Tobacco Guizhou Industrial Co.,Ltd.,Guiyang 550009,China2.School of Public Health,Guizhou Medical University,Guiyang 550025,ChinaAbstract：Near-infrared(NIR)spectroscopy combined with net analyte signal(NAS)processing wereused to evaluate the uniformity of tobacco casing.The NIR spectra from processed tobacco sampleswith casings were orthogonally projected into the spectral space spanned by the samples withoutcasing to obtain the NAS of the casing.The similarity between the samples with and without casingswas then evaluated by calculating the norm of the obtained NAS.The results showed that:1)TheNAS norm was positively correlated with the mass fraction of the casing applied.2)By on-linemeasurement of NIR spectra of tobacco samples with casings in the production line,the NAS normof each processed tobacco sample and its corresponding coefficient of variation were calculated.Thecasing uniformity coefficient obtained had the same trend as that measured by the standard industrialmethod.The NAS analysis method based on NIR spectra was easy to operate and efficient,and itcould be used in the rapid evaluation of tobacco casing uniformity.3)The use of NAS analysismethod to evaluate the casing uniformity of tobacco samples under different spacer sleeve and valvecore positions could achieve the optimization of the spacer sleeve and valve core position of thespraying nozzle in the casing machine.第 56 卷第 9 期加料是卷烟加工过程的重要环节之一，不仅可以改进烟叶的感官性能，还能平衡烟叶中的内在化学成分，调整烟香并赋予品牌卷烟的风格特征1-3。因此，加料工序过程控制对卷烟产品的质量稳定性有着极为重要的作用。关于烟叶加料均匀性评价的方法很多，例如传统的检测料液流量的稳定性方法、标记物法4-5、近红外光谱法6-7、图像法8、过程能力指数法9等。其中，标记物法是评价加料均匀性应用较广的方法之一。烟草行业标准10中常采用1，2-丙二醇作为标记物，通过检测其含量，计算加料均匀系数，以评价加料的均匀性。该方法检测耗时长、过程繁琐，且要求料液中必须含有一定量的1，2-丙二醇，而实际料液成分不是均能满足上述条件，很难在实际生产中真正发挥作用。近红外光谱技术具有简单、快速、非接触测量以及对产品无损等优点，不同物质在近红外线区域多有丰富的吸收光谱，每种质量特征都对应特定的吸收特性，与 C-H、N-H 及O-H等分子基团相关联，因此被广泛用于卷烟制丝工艺过程中的检测分析。目前，近红外光谱技术在加料过程均匀性评价方面的应用主要有两种方式。一种是运用偏最小二乘法建立料液与烟叶比例的定量关系模型，再通过测量实际加料后的各取样点烟叶近红外光谱，计算出各样品中加料量，进而评价烟叶加料的均匀性6。另一种是选择烟叶中的特定标记物，构建其含量的预测模型，再检测加料过程中施加到烟叶表面的该标记物含量，计算含量的变异系数，以此评价加料过程的均匀性7。这两种方法都需要定量分析料液成分，继而才能评价加料均匀性，均涉及标记物选择和建模等问题，前期工作量大，模型亦要求及时维护与更新，因此，方法的检测效率受到制约。鉴于此，基于净信号方法11-12，对加料后样品的近红外光谱进行投影分析，得到不含烟叶、噪声背景等无关信息的料液光谱的净信号，再利用净信号的范数作为样品相似性判断依据，表征加料的均匀性，旨在建立快速高效的卷烟加料均匀性评价方法，为加料过程参数的调整和研究提供技术支撑。1材料与方法1.1材料和仪器烟叶样品和料液样品来自贵州中烟工业有限责任公司生产线上的“贵烟（硬黄精品）”卷烟。Thermo Antaris傅立叶近红外光谱仪（美国Thermo Scientific公司）；QS-5A烟叶切丝机（开封捷利美嘉机械设备有限公司）。1.2方法1.2.1取样方法根据烟叶在加料机上的运行速度，确定加料前后取样的时间间隔为 3 min 21 s，即经过该时间间隔，烟叶从加料工序的入口处运行至加料工序的出口处。因此，加料设备稳定运行后，按照取样间隔4 min，采集加料前后的烟叶样品各 11 份，每份样品的取样量约为100 g，然后将每份样品切丝后分为6份，待测。1.2.2光谱采集在光谱采集过程中，环境温度控制在（222），相对湿度控制在40%10%。近红外光谱仪开机预热时间不低于1 h，然后用仪器自带的自检程序校验，合格后采集各烟丝样品的近红外光谱。将装有烟丝（压平后，样品厚度 20 mm）的样品杯放置于旋转台，采用积分球漫反射的方式测量光谱，扫描范围为4 00010 000 cm-1，分辨率为8 cm-1，扫描次数为64次。每次采集样品光谱前均需测量背景光谱。1.2.3净信号净信号定义为样品信息中待分析成分的有用组分信息，即光谱中与其他干扰信息无关而只与样品分析组分相关的光谱信号11-13。对于样品中的某一待分析成分而言，其净信号是唯一的，不包含其他组分的干扰信息，图1为净信号分析的原理示意图。物质的净信号可以准确地描述光谱与被测成分浓度之间的定量关系。假设样品的光谱数据矩阵Xmn（其中，m表示样品数量，n表示波长或者波数）由待分析指标数据矩阵和与分析指标无关的数据矩阵两部分构成，那么加料后烟丝的光谱X可表示为加料前烟丝的光谱X-k与料液的光谱Xk之和，即：X=Xk+X-k（1）Keywords：Tobacco leaf;Casing;Uniformity;Near-infrared spectroscopy;Net analyte signal胡芸，等：基于近红外光谱净信号的烟叶加料均匀性检测图1净信号分析的原理示意图Fig.1The principle of net analyte signal analyzing 792023 年烟草科技净信号算法就是找到一个包含所有干扰信息的子空间X-k，将加料后烟丝样品的光谱信号X正交投影于该空间，以消除干扰信息的影响，从而得到净信号Xk，即表征待分析指标料液的真正信息。净信号是正交于加料前烟丝光谱子空间信息的，因此，料液的净信号（NAS）可以表示为：NASxk，i=（I-X-k+X-k）xi（i=1，m）（2）式中：I为单位矩阵；X-k为加料前烟丝样品光谱组成的矩阵；X-k+为X-k的广义逆矩阵；xi为第i个加料后烟丝样品的光谱。若净信号不受其他非目标因素的影响，则采用其2范数来衡量料液的净信号大小。这样，将原来的全谱矢量信息xi简化为确定的范数标量ri，大大减少了模型的复杂度，也可以用于表征加料后样品间的相似度。ri=|NASxk，i|（i=1，m）（3）1.2.4加料均匀系数加料均匀系数由式（4）计算：i=1m（ri-r）2（m-1）CU=1-r（4）式中：CU表示加料均匀系数；m表示加料后烟丝样品的数量；ri为第i个加料后烟丝样品光谱净信号的2范数；r为m个加料后烟丝样品光谱净信号范数的平均值。CU值越大，表示加料均匀性越好。2结果与讨论2.1加料前后烟丝样品的近红外光谱烟用料液中含有羟基化合物，其官能团在近红外光谱区具有特征吸收，导致施加料液后烟丝样品的近红外光谱发生显著变化，如图2。随着加料比例的增加，谱图的吸收强度逐渐增大，谱图的基线也向上漂移。为减少或消除烟丝样品固体颗粒大小、表面散射效应引起的基线漂移对光谱数据的影响，利用标准正态变量变换法（SNV）对原始光谱数据进行预处理。图3为料液施加比例为3%（质量分数，下图2施加不同比例料液后烟丝样品的近红外图谱Fig.2NIR spectra of cut tobacco samples with casings added at different ratios图3加料前后烟丝样品（料液施加比例3%）和料液样品的预处理后的图谱Fig.3Pre-processed NIR spectra of cut tobacco samples with/without casing at the addition ratio of 3%and the casing sample 80第 56 卷第 9 期胡芸，等：基于近红外光谱净信号的烟叶加料均匀性检测图4加料前后烟丝样品近红外光谱数据的主成分得分图（a）和载荷图（b）Fig.4Principal component scores（a）and loadings（b）of the pre-processed NIR spectraof cut tobacco samples with/without casing图5不同施加比例料液的烟丝样品的净信号Fig.5NAS of cut tobacco samples with casings added at different ratios同）的加料前后烟丝样品的近红外光谱。从图3中可以看出，在 6 5007 100 cm-1与 4 4005 400 cm-1区域，其谱带强度变化明显。料液样品在该区域有较强的光谱吸收，O-H键的伸缩振动在5 155 cm-1的吸收处于饱和状态。对预处理后的加料前后的烟丝样品（各66个）光谱数据进行主成分分析，结果显示加料前后烟丝样品的主成分得分区别明显（图4a），对应的第1载荷矢量轮廓反映出加料后样品光谱的变化主要在6 5007 100 cm-1与4 4005 400 cm-1区域（图4b）。这说明烟丝中料液含量不同，将直接影响近红外光谱吸收谱带的变化，为加料均匀性的检测与判断提供了可能性。2.2净信号分析方法的建立与验证净信号分析方法的本质是寻找一个包含除料液组分外的各因素构成的子空间X-k。这里采用主成分分析的方法来构建X-k。首先，对加料前的样品光谱数据进行主成分分析，确定主成分数并计算各主成分的累计方差贡献率，使得前q个主成分的累计方差贡献率达到99.0%以上。选取前q个主成分对应的得分矩阵Tq和载荷矩阵Pq构建投影矩阵X-k。矩阵X-k有助于消除测量误差与噪声的影响。对空白烟丝样品分别施加不同比例（2%、3%、5%、10%、20%、30%）的料液，并检测加料后烟丝样品的近红外光谱。按照公式（1），计算加料后烟丝样品的净信号（图5）。从图5中可以看出，施加不同比例料液烟丝样品的净信号矢量变化明显，尤其在6 5007 100 cm-1、4 9005 400 cm-1波数区域与料液的施加比例呈正相关变化；在5 4406 500 cm-1、4 0504 900 cm-1波数区 812023 年烟草科技域内呈负相关变化。计算净信号的2范数，结果如图6。由图6可以看出，净信号范数的大小随着烟丝中料液施加比例的增加而增大。空白烟丝样品净信号范数的平均值标准偏差为0.080.02。当料液施加比例为5%时，其净信号范数值增加到1.000.07。这表明，净信号的2范数能表征料液施加比例的大小。因此，可以通过净信号的2范数来衡量烟丝样品中料液的含量。在构建投影矩阵X-k中，主成分数的选择极为重要。如果选择过少的主成分数，投影矩阵不能合理包含加料前烟丝的光谱信息；若主成分数目过大，会带入不必要的噪声和误差信息，得到的净信号也不准确。针对生产线上两个批次烟丝样品的加料均匀性数据，采用留一法交叉验证（leave-one-out crossvalidation）主成分分析方法，前9个主成分的数据方差累计贡献率达99.10%（图7）。同时，对应的交互验证均方根误差（RMSECV）数值也较小，其数值由此趋于稳定。因此，选取9个主成分构建净信号投影矩阵。由此计算其加料均匀系数分别为 0.895 与0.938。通过行业标准方法测定加料后烟丝样品中1，2-丙二醇含量，得到的加料均匀系数为0.919和0.939。两种方法的检测结果具有相同的趋势，验证了该方法能有效表征烟丝加料的均匀性，但更加简单、快捷，且不需要构建特定标记物的定量模型。图6烟丝样品净信号的2范数与料液施加比例的关系Fig.6Relationship between the 2-norm of NAS of cuttobacco sample and casing addition ratio图7烟丝样品光谱主成分的累计方差贡献率和RMSECV图Fig.7Cumulative variance contribution rate and RMSECV of principal components of spectra of cut tobacco samples2.3净信号分析方法的应用在加料过程中，料液施加的均匀度直接取决于双介质混合喷嘴的喷射角度和雾化效果。加料喷嘴的调校是通过调节阀芯与间套的位置来实现的。调节喷嘴的阀芯能改变喷射区域锥面的大小和雾化效果，可使到达的料液不被作为废气抽走并均匀地施加到烟叶上。实验通过调节加料机喷嘴间套和阀芯的位置，收集不同位置下加料前后的烟叶样品。图8为加料机喷嘴间套与阀芯的位置示意图，喷射出口方向“+”为正值，反方向“-”为负值。设置3个阀芯位置（0、1.0、1.5 mm处），2个间套位置（-1.0和0 mm处），考察烟叶加料效果。分别获取不同组合位置 间套位置（mm），阀芯位置（mm）加料前样品的光谱数据，用10个主成分构建投影矩阵X-k，计算不同组合位置样品的净信号范数，分别得到 6 个组合位置（-1.0，0）、（-1.0，1.0）、（-1.0，1.5）、（0，0）、（0，1.0）、（0，1.5）样品的加料均匀系数分别为0.922、0.927、0.879、0.935、0.928和0.903。该 82第 56 卷第 9 期结果表明喷嘴在组合位置（-1.0，0）、（-1.0，1.0）、（0，0）及（0，1.0）时，烟叶的加料效果都比较好，故均可作为加料喷嘴的调校位置。3结论 基于加料前后烟丝样品近红外光谱的变化，建立了光谱净信号分析方法，计算料液光谱净信号范数，该指标可以表征加料过程中料液施加的均匀性。基于净信号范数计算的加料均匀系数与行业标准方法所测结果具有相同趋势。该方法既不需要寻找标记物也无需建立近红外定量分析模型，具有简单、快速、高效的特点。通过调校加料机喷嘴间套和阀芯的位置，收集不同位置下样品的近红外光谱，运用净信号方法分析加料均匀性，能够实现加料机喷嘴的间套和阀芯的位置优化，有效保障烟叶的均匀加料。参考文献1 庞红蕊，杨佳玫，张志航，等.卷烟加工过程加料均匀性研究进展 J.轻工科技，2019，35（5）：23-25.PANG Hongrui，YANG Jiamei，ZHANG Zhihang，et al.Research progress on casing uniformity during cigarettemanufacturing J.LightIndustryScienceandTechnology，2019，35（5）：23-25.2孙志涛，魏明杰，王东，等.叶丝加料和片烟加料对一、二、三类卷烟感官质量的影响 J.烟草科技，2017，50（1）：74-79.SUN Zhitao，WEI Mingjie，WANG Dong，et al.Effectsof cut strip casing and strip casing on sensory qualityof cigarettes of Category ，and J.TobaccoScience&Technology，2017，50（1）：74-79.3廖惠云，张映，郝喜良.卷烟制丝加料工序均匀性的评价与控制 J.烟草科技，2008（8）：5-7，21.LIAO Huiyun，ZHANG Ying，HAO Xiliang.Evaluationand control of casing uniformity in primary processingJ.Tobacco Science&Technology，2008（8）：5-7，21.4 温若愚，赵维一，曾建，等.卷烟加料工序均匀性的测定方法 J.江西农业学报，2010，22（4）：36-38，42.WEN Ruoyu，ZHAO Weiyi，ZENG Jian，et al.Methodfor measuring uniformity of cigarette casing processJ.Acta Agriculturae Jiangxi，2010，22（4）：36-38，42.5刘泽，朱勇，邱昌桂，等.储叶过程中片烟吸收料液效果的评价方法 J.中国烟草学报，2015，21（1）：22-27.LIU Ze，ZHU Yong，QIU Changgui，et al.Method ofevaluating casing absorption by tobacco strips duringstoring J.Acta Tabacaria Sinica，2015，21（1）：22-27.6 于海洋，丁超，王宇，等.基于近红外光谱的卷烟加料均匀性检测 J.农业工程，2014，4（1）：50-52，58.YU Haiyang，DING Chao，WANG Yu，et al.Test ofcigarette charging uniformity based on near-infraredspectrum J.Agriculture Engineering，2014，4（1）：50-52，58.7 韩明，王建民，朱建新，等.烟片含水率对丙二醇近红外预测值及加料均匀性系数的影响 J.中国农学通报，2018，34（32）：150-156.HAN Ming，WANG Jianmin，ZHU Jianxin，et al.EffectofmoisturecontentonNIRpredictedvalueofpropylene glycol in tobacco leaf and casing uniformitycoefficient J.Chinese Agricultural Science Bulletin，2018，34（32）：150-156.8 郑飞，李媛，刘德强，等.基于图像处理的烟片加料均匀性评价方法 J.烟草科技，2015，48（11）：65-68.ZHENG Fei，LI Yuan，LIU Deqiang，et al.Method forevaluating casing uniformity of tobacco strips based onimage processing J.Tobacco Science&Technology，2015，48（11）：65-68.9 刘泽春，黄华发，洪伟龄，等.过程能力分析在卷烟加料均匀性评价中的应用 J.烟草科技，2010（10）：5-7，21.LIU Zechun，HUANG Huafa，HONG Weiling，et al.Application of process capability analysis in evaluationof casing uniformityJ.Tobacco Science&Technology，2010（10）：5-7，21.10 YC/T 3532010卷烟加料均匀性的测定 S.YC/T 3532010CigarettesDetermination of casinguniformity S.11 韩广，刘蓉，徐可欣.近红外无创血糖检测中基于差动式浮动基准测量的有效信号提取 J.光谱学与光谱分析，2018，38（5）：1599-1604.HAN Guang，LIU Rong，XU Kexin.Extraction of胡芸，等：基于近红外光谱净信号的烟叶加料均匀性检测图8加料机喷嘴间套和阀芯的位置图Fig.8Position of spacer sleeve and valve core of thespraying nozzle in a casing machine（下转第90页）832023 年烟草科技nanofibrillated cellulose to increase the precipitatedcalcium carbonate retention and physical propertiesduringreconstitutedtobaccosheetpreparation J.Industrial Crops and Products，2019，130：592-597.17 GB/T 292852012纸浆实验室湿解离机械浆解离 S.GB/T 292852012PulpsLaboratory wet disintegrationDisintegration of mechanical pulps S.18 GB/T 53992004纸浆浆料浓度的测定 S.GB/T53992004PulpsDetermination of stockconcentration S.19 GB/T 292872012纸浆实验室打浆PFI 磨法S.GB/T 292872012PulpsLaboratory beatingPFImill method S.20 GB/T 33322004纸浆打浆度的测定（肖伯尔-瑞格勒法）S.GB/T 33322004PulpsDetermination of beatingdegree（Schopper-Riegler method）S.21 GB/T 164472004烟草及烟草制品调节和测试的大气环境 S.GB/T 164472004Tobacco and tobacco productsAtmosphere for conditioning and testing S.22 GB/T 451.22002纸和纸板定量的测定 S.GB/T 451.22002Paper and boardDeterminationof grammage S.23 GB/T 451.32002纸和纸板厚度的测定 S.GB/T 451.32002Paper and boardDeterminationof thickness S.24 金岚峰，蒋宇凡，薛冬，等.造纸法再造烟叶原料的打浆性能 J.烟草科技，2014（8）：10-13.JINLanfeng，JIANG Yufan，XUEDong，etal.Performance of raw material beating in paper-makingprocess reconstituted tobacco productionJ.TobaccoScience&Technology，2014（8）：10-13.25 俞京，袁广翔，蒋宇凡，等.一种多元助留体系在再造烟叶生产中的应用 J.烟草科技，2021，54（4）：81-88.YU Jing，YUAN Guangxiang，JIANG Yufan，et al.Application of multiple retention system to reconstitutedtobacco production J.Tobacco Science&Technology，2021，54（4）：81-88.责任编辑周雅宁（上接第83页）effective signal in non-invasive blood glucose sensingwith near-infrared spectroscopy J.Spectroscopy andSpectral Analysis，2018，38（5）：1599-1604.12 耿颖，相秉仁，何兰.净分析信号算法用于近红外模型优化的研究 J.光谱学与光谱分析，2015，35（10）：2730-2733.GENG Ying，XIANG Bingren，HE Lan.Study on theapplication of NAS-based algorithm in the NIR modeloptimization J.Spectroscopy and Spectral Analysis，2015，35（10）：2730-2733.13 Lorber A，Faber K，Kowalski B R.Net analyte signalcalculation in multivariate calibrationJ.AnalyticalChemistry，1997，69（8）：1620-1626.责任编辑周雅宁 90