CN102879532B - 一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 - Google Patents
一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 Download PDFInfo
- Publication number
- CN102879532B CN102879532B CN201210384570.6A CN201210384570A CN102879532B CN 102879532 B CN102879532 B CN 102879532B CN 201210384570 A CN201210384570 A CN 201210384570A CN 102879532 B CN102879532 B CN 102879532B
- Authority
- CN
- China
- Prior art keywords
- node
- coefficient
- layer
- acid
- hidden layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Manufacture Of Tobacco Products (AREA)
Abstract
本发明涉及一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法。该方法是:将待测初烤烟叶通过人工撕片、剔梗、切丝三个步骤进行样品前处理;对待测样品7项化学成分(钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸)进行检测;将7项化学成分测定结果结合模型各输入层系数计算隐含层7个节点的网络值;将隐含层7个节点的网络值转换为隐含层7个节点的输出值;将隐含层7个节点输出值结合模型输出层系数计算得到烟气苯并[a]芘的释放量预测值。能够由烤烟叶原料对未来形成的卷烟成品中的苯并[a]芘累积含量进行有效预判,从而指导生产过程中的原料选择。
Description
技术领域
本发明涉及一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法,属于烟草技术领域。
背景技术
苯并[a]芘是卷烟主流烟气中的重要有害物质,被列为评价卷烟安全性的7种卷烟烟气有害成分之一。目前我国已经建立了卷烟主流烟气中苯并[a]芘测定的国家标准。但是通常测定烟气中苯并[a]芘的含量需要对样品进行卷制、吸烟机抽吸捕集及捕集物进样前处理等较为烦琐的步骤。相对而言,烟草中的一些常量化学成分的测定步骤就较为简便。而烟气中的苯并[a]芘是一种典型的稠环芳烃(简称PAH),稠环芳烃是烟草在高温缺氧条件下不完全燃烧的产物。各种有机物热解所生成的有机物碎片,经过复杂的聚合过程而形成,这些烟草化学成分变化最终会对烟气成分含量造成影响。由此可知,在烟草化学成分通过燃烧转换为烟气成分的过程中经历了一系列复杂的化学反应和物质分解、生成过程,这一切过程都属于一个复杂的非线性变化体系。但是从大量文献及相关资料来看,对于这种复杂变化体系间的关系研究还较少。
发明内容
本发明的目的在于提供一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法,对初烤烟叶的几项化学成分(由本方法确定)进行测定后,通过本方法预测其烟气中的苯并[a]芘释放量,以便于卷烟配方人员及时了解该原料的苯并[a]芘释放量水平,并在后期的选料和叶组配方设计中有效应用,从而达到选料降低苯并[a]芘释放量的一种减害措施。
现代神经网络是一种非线性统计性数据建模工具,常用来对输入和输出间复杂的关系进行建模,网络自身通常都是对自然界某种算法或者函数的逼近。通过利用大量样本所构建的烟草化学成分与烟气苯并[a]芘释放量神经网络模型,对待测样品烟气苯并[a]芘释放量水平进行预测,从而避免了进行烟气检测过程中的一些设备要求和烦琐处理步骤。
本发明提出的初烤烟叶主流烟气中苯并[a]芘释放量的预测方法具体内容如下:
一、预测方法的建立
(1)建模样品来源
2009年的186个初烤烟叶样品,覆盖全国47个产地,9个品种,上、中、下3个烟叶部位。
(2)建模样品前处理
将待测初烤烟叶进行样品前处理通过人工撕片、剔梗、切丝三个步骤。处理好的烟片(丝)马上封口、贴上标签,统一存放,不混淆。不同烟样回潮和切丝前后都对设备进行认真地清洁,以保证烟样间不错串。
所收集的初烤烟叶原料统一加工,统一编码,单一原料分别切丝并卷制成单料烟支(为尽量消除卷烟辅料对研究的干扰,选择统一透气度为60CU、克重28g/m2的同批次卷烟纸,将烟叶样品卷制成无滤嘴卷烟,并以烟支重量指标挑选烟支),不加香加料,挑选后的烟支贮藏在温度(18±1)℃;湿度(50±10)%的环境中,分析前取出按烟气国标方法平衡48小时备用。
(3)建模样品化学成分和烟气苯并[a]芘释放量测定
由于某些烟草化学成分是生成烟气成分的前体化合物(如碳水化合物、含氮化合物、有机酸、酚类等),其含量对烟气成分的含量影响较大,因此选取25种初烤烟叶化学成分(总糖、还原糖、氯、钾、总氮、烟碱、挥发碱、水分、蛋白质、挥发酸、纤维素、绿原酸、莨菪亭、芸香苷、总多酚、草酸、丙二酸、丁二酸、苹果酸、柠檬酸、棕榈酸、亚油酸、油酸、亚麻酸、硬脂酸)作为所要研究的指标。通过YC/T159-2002连续流动法测定水溶性糖;通过YC/T161-2002连续流动法测定总氮;通过YC/T217-2007连续流动法测定钾;通过YC/T162-2002连续流动法测定氯;通过GB/T 23225—2008光度法测定烟碱;YC/T 35-1996返滴定法测定挥发碱;YC/T 31-1996烘箱法测定水分;YC/T 166-2003克达尔法测定蛋白质;YC/T 202—2006高效液相色谱法测定多酚化合物(绿原酸、莨菪亭、芸香苷、总多酚);采用水蒸气蒸馏返滴定法测定测定挥发酸;烟草中粗纤维的测定的方法测定纤维素;微波辅助衍生化气相色谱法测定非挥发性有机酸(草酸、丙二酸、丁二酸、苹果酸、柠檬酸、棕榈酸、亚油酸、油酸、亚麻酸、硬脂酸),测定出的烟草化学成分单位换算为%。GB/T 21130-2007气相色谱法测定烟气苯并[a]芘释放量,测定出的烟气苯并[a]芘释放量换算为每克烟丝的烟气苯并[a]芘释放量为ng/g。
(4)建模变量和异常样品筛选
变量筛选一方面可以简化模型,另一方面是把对烟气成分影响较小的变量剔除,使模型的预测能力更强、稳健性更好。本研究变量选择的原则是通过遗传算法进行变量筛选,并根据以往的机理研究结果最大程度的保留产生烟气成分的前体化合物变量,最终确定合理的建模变量组成为7个化学成分(钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸)。
异常样品是远离模型整体的样品,对模型的回归分析影响明显,首先通过浓度残差分析剔出,然后在建模过程中重复运算,再次剔除,直到获得最佳预测效果,最终确定的进入建模的样本有97个。
(5)建模参数优化和模型内部性能评价
采用神经网络方式建模利用重复抽样方法对权衰减系数和隐含层节点数进行优化选择,最终参数确定为表1。采用所构建的预测模型的决定系数(R2)和预测标准偏差(SEC)(见式1)评价预测模型内部性能,决定系数(R2)越大,预测标准偏差(SEC)越小,则模型越好,详见表2。
表1模型参数、权重数及网络结构
表2模型内部性能指标
式(1)中,yi,actual为第i建模样品实测值,yi,perdicted为用所建模型对建模样本中第i样品预测值,n为建模样本个数。
(6)模型外部验证
采用未参与建模的28个外部验证样品对模型的预测能力进行验证。其中采用①SEP/SEC、②RPD、③配对t检验、④平均预测相对偏差,四种参数来模型的泛化能力和适用性进行验证,见表3。
①SEP/SEC
验证样本的预测标准偏差(SEP)(见式2)与建模样本的预测标准偏差(SEC)的比值小于等于1.2,即SEP不能大于1.2倍的SEC,以此来说明模型是否过拟合。
式(2)中,yi,actual为第i验证样品实测值,yi,perdicted为用验证样本预测过程中第i验证样品预测值,m为验证样本个数。
②RPD
SDv为所有验证样本实测值的标准偏差,在SDv相同的前提下,验证样本决定系数(R2)越大,模型预测准确性越高;SDv与验证样本的预测标准偏差(SEP)的比值RPD(见式3),RPD越大,模型预测准确性越高。通常认为若RPD<2,则表明预测结果是不可接受的;
式(3)中,SDv为所有验证样本实测值的标准偏差,SEP越小,RPD越大。
③配对t检验
预测验证结果通过配对t检验进行验证,当显著性水平大于0.05时,t的绝对值小于其相关的临界值(t0.05,27=2.051831),表明标准测量方法与预测方法不存在显著性差异,即两种方法的测定结果不存在系统误差。
④平均预测相对偏差
计算验证样本实测值与预测值的平均预测相对偏差(见式4),平均预测相对偏差越小,模型预测准确性越高。
式(4)中,yi,actual为第i验证样品实测值,yi,perdicted为用验证样本预测过程中第i验证样品预测值,m为验证样本个数。
表3模型外部验证情况汇总
二、预测方法的应用
(1)待测样品前处理
将待测初烤烟叶进行样品前处理通过人工撕片、剔梗、切丝三个步骤。处理好的烟片(丝)马上封口、贴上标签,统一存放,不混淆。不同烟样回潮和切丝前后都对设备进行认真地清洁,以保证烟样间不错串。
(2)待测样品7种化学成分测定
对7种初烤烟叶化学成分(钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸)通过YC/T217-2007连续流动法测定钾;通过微波辅助衍生化–气相色谱法测定非挥发性有机酸(棕榈酸、丙二酸、硬脂酸、苹果酸);通过YC/T 202—2006高效液相色谱法测定多酚化合物(绿原酸、莨菪亭),其测定单位均换算为%。
(3)待测样品通过模型系数预测苯并[a]芘释放量
通过建立的三层结构网络模型:第一层为输入层,节点数为7,对应输入变量为钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸;第二层为隐含层,节点数为7;第三层为输出层,节点数为1,对应输出变量为苯并[a]芘释放量。
表4输入层到隐含层各节点系数值
隐含层网络值按式(5)计算。
xj=ZW输入(其中j=[1,2,…,7]) (5)
式中:
Z-表示输入向量[烟丝钾,烟丝棕榈酸,烟丝丙二酸,烟丝莨菪亭烟丝硬脂酸,烟丝苹果酸,烟丝绿原酸,b],其中b为输入层偏置值,其值为1;
W输入-表示输入层系数矩阵(详见表4);
xj-为隐含层第j个节点处的网络值。
隐含层输出值按式(6)计算。
式中:
e-为自然对数的底数。
Oj-为隐含层第j个节点输出值。
注:式(6)中若xj的值大于15时,Oj按1计,若xj的值小于-15时,Oj按0计。
表5输出层系数值
输出层输出值按式(7)计算。
y苯并[a]芘=OjW输出 (7)
式中:
Oj-表示隐含层输出值向量[O1,O2,O3,O4,O5,O6,O7,b],其中b为输出层偏置值,其值为1;
W输出-为输出层系数向量[W1,W2,W3,W4,W5,W6,W7,W0](详见表5);
y苯并[a]芘-为输出层输出值,即烟气苯并[a]芘的释放量预测值。
按本发明采用初烤烟叶中的7个化学成分对其烟气中苯并[a]芘释放量进行预测,在应用过程中具备以下优点:
A.对待测烟叶样品进行7项烟草化学成分常量检测,即可预测其烟气中苯并[a]芘含量,无需进行烟支卷制及使用吸烟机设备抽吸来捕集粒相物质。
B.应用已构建的模型结构配合模型系数进行分析,只需要进行3个步骤的简单计算即可获得苯并[a]芘的释放量预测值,操作便捷,可采用计算器或是人工计算,无需使用复杂的计算软件。
C.能够依靠目前的初烤烟叶原料对未来形成的卷烟成品中可能的苯并[a]芘累积含量进行有效预判,从而指导生产过程中的原料选择,对提高卷烟产品质量安全具有重要的现实意义。
具体实施方式
实施例:
(1)待测样品前处理
将待测初烤烟叶进行样品前处理通过人工撕片、剔梗、切丝三个步骤。处理好的烟片(丝)马上封口、贴上标签,统一存放,不混淆。不同烟样回潮和切丝前后都对设备进行认真地清洁,以保证烟样间不错串。
(2)选取1份待测初烤烟叶样品进行7种化学成分测定
对7种初烤烟叶化学成分:钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸进行测定,通过YC/T217-2007连续流动法测定钾;通过微波辅助衍生化–气相色谱法测定非挥发性有机酸(棕榈酸、丙二酸、硬脂酸、苹果酸);通过YC/T 202—2006高效液相色谱法测定多酚化合物(绿原酸、莨菪亭),其结果见表6;
(3)待测样品通过模型系数预测苯并[a]芘释放量
通过建立的三层结构网络模型:第一层为输入层,节点数为7,对应输入变量为钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸;第二层为隐含层,节点数为7;第三层为输出层,节点数为1,对应输出变量为苯并[a]芘释放量。
1)将检测所得的7项化学成分结果配合表4中各输入层的系数,按式(5)进行计算得到隐含层7个节点的网络值,其结果见表7;
2)将所得的网络值按式(6)计算得到隐含层7个节点的输出值,其结果见表8;
3)将所得的隐含层输出值配合表5中的系数,按式(7)进行计算得到输出层输出值,即烟气苯并[a]芘的释放量预测值,其结果见表9。
表6待测初烤烟叶样品7项化学成分检测结果
表7待测初烤烟叶样品隐含层7个节点的网络值
表8待测初烤烟叶样品隐含层7个节点的输出值
表9待测初烤烟叶样品输出层输出值(烟气苯并[a]芘释放量预测值)
Claims (2)
1.一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法,其特征在于包含以下步骤:
1)将待测初烤烟叶通过人工撕片、剔梗、切丝三个步骤进行样品前处理;
2)对待测样品7项烟草化学成分进行检测,所述的7项烟草化学成分是钾、棕榈酸、丙二酸、莨菪亭、硬脂酸、苹果酸、绿原酸;
3)将待测样品7项烟草化学成分测定结果结合模型各输入层系数计算隐含层7个节点的网络值;
4)将计算所得隐含层7个节点的网络值转换为隐含层7个节点的输出值;
5)将计算所得隐含层7个节点输出值结合模型输出层系数计算得到烟气苯并[a]芘的释放量预测值;
所述的各输入层系数从节点1至节点7依次为:钾节点系数值-17.4676、35.2475、-1.4618、11.4246、-9.9554、-1.3241、0.0266;棕榈酸节点系数值3.1511、-29.4732、-3.1677、7.4519、-3.2488、-1.5762、0.0189;丙二酸节点系数值1.3113、-2.711、1.6466、1.6945、7.7656、2.4191、0.0175;莨菪亭节点系数值4.7053、0.1296、-5.8362、-3.9998、-0.2499、0.0677、0.0003;硬脂酸节点系数值20.5316、-1.7037、9.8412、-1.7745、1.3413、-0.3899、0.0033;苹果酸节点系数值-0.0494、0.7188、-0.26、-1.9123、0.6012、0.1265、0.2812;绿原酸节点系数值-0.9037、-13.2678、-3.8322、4.6474、6.4448、-2.5278、0.0213;输入层偏置系数节点系数值17.803、19.0837、11.7137、-12.3716、-1.9662、0.0159、0.0131,其中输入层偏置值为1;
所述的输出层系数为:隐含层节点1输出层系数8.7333;隐含层节点2输出层系数8.3806;隐含层节点3输出层系数14.3207;隐含层节点4输出层系数12.4867;隐含层节点5输出层系数13.7848;隐含层节点6输出层系数-1.3254;隐含层节点7输出层系数-0.7097;输出层偏置系数-0.6847,其中输出层偏置值为1。
2.根据权利要求1 所述的初烤烟叶主流烟气中苯并[a]芘释放量的预测方法,其特征在于所述的模型为神经网络模型,结构为7-7-1,即7个输入层节点,7个隐含层节点,1个输出层节点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210384570.6A CN102879532B (zh) | 2012-10-11 | 2012-10-11 | 一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210384570.6A CN102879532B (zh) | 2012-10-11 | 2012-10-11 | 一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102879532A CN102879532A (zh) | 2013-01-16 |
CN102879532B true CN102879532B (zh) | 2015-05-06 |
Family
ID=47480933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210384570.6A Active CN102879532B (zh) | 2012-10-11 | 2012-10-11 | 一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102879532B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316484B (zh) * | 2014-11-13 | 2017-03-29 | 中国日用化学工业研究院 | 一种中、高温煤沥青中苯并(a)芘的检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8600679D0 (en) * | 1986-01-13 | 1986-02-19 | Imp Group Plc | Chemical analysis of tobacco/smoking-related products |
CA2616631C (en) * | 2005-08-02 | 2013-05-14 | Wm. Wrigley Jr. Company | Screening method for the identification of compositions suitable for the treatment of oral cavity malodor associated with smoking a tobacco product |
CN101419207B (zh) * | 2008-10-27 | 2012-09-05 | 川渝中烟工业有限责任公司 | 烤烟烟气主要指标的预测方法 |
-
2012
- 2012-10-11 CN CN201210384570.6A patent/CN102879532B/zh active Active
Non-Patent Citations (3)
Title |
---|
BP神经网络隐含层单元数的确定;沈花玉等;《天津理工大学学报》;20081031;第24卷(第5期);第13-15页 * |
使用BP算法时应考虑的若干问题;郭海涛等;《佳木斯大学学报(自然科学版)》;20001231;第18卷(第4期);第363-365页 * |
关于BP网中隐含层层数及其节点数选取方法浅析;叶斌等;《商丘职业技术学院学报》;20040630;第3卷(第6期);第52-54页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102879532A (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879531B (zh) | 一种初烤烟叶主流烟气中氨释放量的预测方法 | |
Xu et al. | A wavelength selection method based on randomization test for near-infrared spectral analysis | |
CN101419207B (zh) | 烤烟烟气主要指标的预测方法 | |
CN104569263B (zh) | 一种快速准确评价卷烟产品质量稳定性的方法 | |
CN101498658A (zh) | 基于捕集了烟气粒相物的剑桥滤片的傅立叶变换近红外光谱的烟气化学组分预测方法 | |
CN102879530B (zh) | 一种初烤烟叶主流烟气中氢氰酸释放量的预测方法 | |
CN107677638A (zh) | 一种基于近红外光谱技术的初烤把烟含梗率快速检测方法 | |
CN101685092B (zh) | 用烟叶pH检测值判定烤烟醇化质量的方法 | |
CN102854289B (zh) | 一种初烤烟叶主流烟气中苯酚释放量的预测方法 | |
CN102854290B (zh) | 一种初烤烟叶主流烟气中巴豆醛释放量的预测方法 | |
CN110967313A (zh) | 电子烟烟油尼古丁含量的近红外光谱预测建模方法及应用 | |
CN102879532B (zh) | 一种初烤烟叶主流烟气中苯并[a]芘释放量的预测方法 | |
CN105138834A (zh) | 基于近红外光谱波数k均值聚类的烟草化学值定量方法 | |
CN111239316A (zh) | 一种基于烟气全化学成分色谱数据的卷烟感官质量分析方法 | |
CN104596979A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶纤维素的方法 | |
CN102866234B (zh) | 一种初烤烟叶主流烟气中nnk释放量的预测方法 | |
CN103868957B (zh) | 一种评价打叶复烤工序片烟感官品质稳定性的方法 | |
CN104596976A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶蛋白质的方法 | |
Liu et al. | An identification method of herbal medicines superior to traditional spectroscopy: Two-dimensional correlation spectral images combined with deep learning | |
CN108801968B (zh) | 一种基于真空红外光谱的烟丝整体质量特征分析方法 | |
CN104132904B (zh) | 一种测定烟叶燃烧性的方法 | |
CN107029038B (zh) | 一种响应面法优化龙眼叶总黄酮的提取工艺 | |
CN109060716A (zh) | 基于窗口竞争性自适应重加权采样策略的近红外特征光谱变量选择方法 | |
Liu et al. | Research on the online rapid sensing method of moisture content in famous green tea spreading | |
CN116026780A (zh) | 基于串联策略波长选择的包衣吸湿率在线检测方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |