CN103487554B - 一种卷烟主流烟气中焦油量及一氧化碳量测定方法 - Google Patents

Abstract

本发明公开了一种卷烟主流烟气中焦油量及一氧化碳量测定方法，该测定方法在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油量回归方程，对卷烟样品烟气中焦油量进行校正；在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳量回归方程，对卷烟样品烟气中一氧化碳量进行校正。本发明建立大气压力为90kPa和100kPa的实验室烟气分析结果之间的修正关系，通过对结果的修正，提高了烟草检测实验室分析数据的准确性和可比性，对烟草工业企业的产品设计和卷烟盒标标注具有重要指导的意义。

Description

一种卷烟主流烟气中焦油量及一氧化碳量测定方法

技术领域

本发明属于卷烟质量检测技术领域，尤其涉及一种卷烟主流烟气中焦油量及一氧化碳量测定方法。

背景技术

卷烟主流烟气是对卷烟质量进行评价的主要内容之一，随着我国开始履行《烟草控制框架公约》，卷烟烟气分析的准确性对烟气释放物的披露具有非常大的影响。在我国卷烟盒标中要标注：焦油量、烟气一氧化碳量和烟气烟碱量三个重要化学成份，准确检测这三个化学成份是烟草检测实验室的一项重要工作。

由于我国地理环境是东低西高，全国有近三分之一烟草检测实验室地处高海拔地区、海拔波动幅度大，大气压力一般为84.8～91.7KPa，大气压力较低。

而资料显示，大气压力每降低1KPa，空气中每立方米氧含量也同时降低3.219克。卷烟主流烟气分析是一燃吸过程，氧是助燃剂，大气压力对卷烟主流烟气分析结果具有一定的影响，这从每年的国际亚太卷烟烟气共同实验和国内卷烟烟气共同实验的数据中，大气压力偏低的烟草检测实验室分析数据偏低可以得到验证。

探寻在90kPa左右的大气压力条件下，烟气指标受影响的程度和规律，建立大气压力为90kPa和100kPa的实验室烟气分析结果之间的修正关系，通过结果的修正，提高烟草检测实验室分析数据的准确性和可比性，对烟草工业企业的产品设计和卷烟盒标标注具有重要指导的意义。

发明内容

本发明提供了一种卷烟主流烟气中焦油量及烟气一氧化碳量测定方法，旨在解决现有技术在实验室检测分析烟草焦油及烟气一氧化碳数据的准确性和可比性较差，同时检测分析过程中受大气压力条件影响较大的问题。

本发明的目的在于提供一种卷烟主流烟气中焦油量及烟气一氧化碳量测定方法，该测定方法包括以下步骤：

步骤一，在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油含量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正；

步骤二，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正。

进一步，在步骤一中，通过卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量；

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量。

进一步，在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标焦油量小于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数，取Δδ=2.39。

进一步，在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量大于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量大于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标焦油量大于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数，取Δδ=12。

进一步，在步骤二中，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量；

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量。

进一步，在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量小于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ=2.28。

进一步，在非标准大气压下，采集盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量大于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量大于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量大于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ=47。

本发明提供的卷烟主流烟气中焦油量及一氧化碳量测定方法，该测定方法在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正；在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正；本发明建立大气压力为90kPa和100kPa的实验室烟气分析结果之间的修正关系，通过对结果的修正，提高了烟草检测实验室分析数据的准确性和可比性，对烟草工业企业的产品设计和卷烟盒标标注具有重要指导的意义。

附图说明

图1是本发明实施例提供的卷烟主流烟气中焦油量及一氧化碳量测定方法的实现流程图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步的详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定发明。

图1示出了本发明实施例提供的卷烟主流烟气中焦油量及一氧化碳量测定方法的实现流程。

该测定方法包括以下步骤：

步骤S101，在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油含量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正；

步骤S102，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正。

在本发明实施例中，在步骤S101中，通过卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量；

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量。

在本发明实施例中，在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标焦油量小于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ为焦油量校正常数，取Δδ=2.39。

在本发明实施例中，在非标准大气压下，采集盒标焦油量大于等于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量大于等于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量大于等于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量大于等于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标焦油量大于等于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数，取Δδ=12。

在本发明实施例中，在步骤S102中，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量；

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量。

在本发明实施例中，在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量小于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ=2.28。

在本发明实施例中，在非标准大气压下，采集盒标一氧化碳量大于等于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量大于等于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量大于等于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量大于等于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量大于等于11mg/支的卷烟样品运用

LINEST函数进行回归分析的数据代入多重回归公式：

y=m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y=m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量大于等于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ=47。

下面结合附图及具体实施例对本发明的应用原理作进一步描述。

实测焦油量小于8mg/支的样品的烟气焦油量校正

分段对盒标焦油量小于8mg/支的样品运用LINEST函数进行回归的分析数据见表1.1所示。

表1.1焦油量小于8mg/支样品的烟气焦油含量多元分析数据表

将表1.1所得的结果保留到小数点后四位，代入多重回归公式y=m1*x1+m2*x2+b可得到：

Y＝0.4241X₁+2.7210X₂-0.9582X₃-2.3899     （1.1）

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式（1.1）系数保留小数点后两位，得该段焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(2.72-0.96\sqrt{{\mathrm{TAR}}_0})+0.42\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ}---\left(1.2\right)$

式中：P₀—标准大气压或指定大气压（指定中心站大气压为99.9kPa）,kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值（修正值），mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数（推荐Δδ=2.39）。

将各站测定值带入公式1.2，修正烟气中焦油量结果见下表1.2所示。

表1.2焦油量小于8mg/支样品的烟气焦油含量校正值

由表1.2的修正结果可以看出，利用公式（1.2）对烟气的焦油含量进行修正其最小偏差为-0.05%，最大偏差为2.44%。运用公式（1.2）所得的焦油量校正值与2实验室的测定值十分吻合。

实测焦油量大于等于8mg/支的样品（含8mg/支）的烟气焦油量校正

分段对盒标焦油量大于等于8mg/支的样品运用LINEST函数进行回归的分析数据见表1.3所示。

表1.3焦油量大于等于8mg/支样品的烟气焦油量多元分析数据表

将表1.3所得的结果保留到小数点后四位，代入多重回归公式y=m1*x1+m2*x2+b可得到：

Y＝-0.1491X₁-2.4078X₂+0.8214X₃+12.4557    （1.3）

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{TAR}}_0*\frac{P_0}{P};X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式（1.1）系数保留小数点后两位，得该段焦油回归方程：

$\mathrm{TAR}={\mathrm{TRA}}_0[\frac{P_0}{P}(0.82\sqrt{{\mathrm{TAR}}_0}-2.41)-0.15\sqrt{{\mathrm{TAR}}_0}]+\mathrm{\Δ\δ}---\left(1.4\right)$

式中：P₀—标准大气压或指定大气压（指定2实验室大气压为99.9kPa）,kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值（修正值），mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油校正常数（推荐Δδ=12）。

将各站测定值带入公式1.4，修正烟气中焦油量结果见下表1.4所示。

表1.4焦油量大于等于8mg/支样品（含8mg/支）的烟气焦油量校正值

由表1.4的修正结果可以看出，利用公式（1.4）对烟气的焦油量进行修正其最小偏差为0.02%，最大偏差为-3.02%。运用公式（1.4）所得的焦油量校正值与2实验室的测定值十分吻合。

非标准大气压下烟气CO测定结果修正公式

非标准大气压实验室的数据按照CO量11mg/支分两段进行修正（按盒标含量分段）。结合经验公式：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}-(\frac{P_0}{P}-1)\sqrt{{\mathrm{CO}}_0/C_{\mathrm{CO}}}],$

设定自变量 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ 因变量Y＝CO，利用表1.4中的LINEST函数寻找自变量X1,X2,X3与Y的回归关系。

烟气CO量小于11mg/支的样品的烟气CO量校正

分段对盒标CO量小于11mg/支的样品运用LINEST函数进行回归的分析数据见表2.1所示。

表2.1烟气CO量小于11mg/支样品的烟气CO量多元分析数据表

将表2.1所得的结果保留到小数点后四位，代入多重回归公式y=m1*x1+m2*x2+b可得到：

Y＝0.2459X₁+2.1802X₂-0.5624X₃-2.2812   （2.1）

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式（2.1）系数保留小数点后两位，得该段焦油回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(2.18-0.56\sqrt{{\mathrm{CO}}_0})+0.25\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ}---\left(2.2\right)$

式中：P₀为标准大气压或指定大气压（指定2实验室大气压为99.9kPa）,kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值（修正值），mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳校正常数（推荐Δδ=2.28）。

将各站测定值带入公式2.2，修正烟气中CO含量结果见下表2.2所示。

表2.2CO量小于11mg/支样品的烟气CO量校正值

由表2.2的修正结果可以看出，利用公式（2.2）对烟气的CO含量进行修正其最小偏差为0.13%，最大偏差为-2.60%。运用公式（2.2）所得的CO量校正值与2实验室的测定值十分吻合。

烟气CO量大于等于11mg/支的样品的烟气CO量校正

分段对盒标CO量大于等于11mg/支的样品运用LINEST函数进行回归的分析数据见表2.3所示。

表2.3烟气CO量大于等于11mg/支样品的烟气CO量多元分析数据表

将表2.3所得的结果保留到小数点后四位，代入多重回归公式y=m1*x1+m2*x2+b可得到：

Y＝-0.6515X₁-7.0927X₂+1.9115X₃+46.8898    （2.3）

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};X_2={\mathrm{CO}}_0*\frac{P_0}{P};X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式（2.1）系数保留小数点后两位，得该段焦油回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(1.91\sqrt{{\mathrm{CO}}_0}-7.09)-0.65\sqrt{{\mathrm{CO}}_0}]+\mathrm{\Δ\δ}---\left(2.4\right)$

式中：P0为标准大气压或指定大气压（指定2实验室大气压为99.9kPa）,kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值（修正值），mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳校正常数（推荐Δδ=47）。

将各站测定值带入公式2.4，修正烟气中CO量结果见下表2.4所示。

表2.4烟气CO量大于等于11mg/支样品的烟气CO量校正值

由表2.4的修正结果可以看出，利用公式（2.4）对烟气的一氧化碳含量进行修正其最小偏差为0.07%，最大偏差为3.07%。运用公式（2.4）所得的CO量校正值与2实验室的测定值均比较吻合。

对上述线性回归、多元分段回归校正后相对偏差分析。采用多元分段回归分析、分析数据高于线性回归分析，相对偏差也小于线性回归分析，主流烟气焦油量、CO量通过修正方程修正值与2实验室的测定值均达到满意结果。

本发明实施例提供的卷烟主流烟气中焦油量及一氧化碳量测定方法，该测定方法在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油含量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正；在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳含量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正；本发明建立大气压力为90kPa和100kPa的实验室烟气分析结果之间的修正关系，通过对结果的修正，提高了烟草检测实验室分析数据的准确性和可比性，对烟草企业的产品设计和卷烟盒标标注具有重要指导的意义。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种卷烟主流烟气中焦油量及一氧化碳量测定方法，其特征在于，该测定方法包括以下步骤：

步骤一，在非标准大气压下，采集卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油量回归方程，对卷烟样品烟气中焦油量进行校正；

步骤二，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳量回归方程，对卷烟样品烟气中一氧化碳量进行校正；

在步骤一中，通过卷烟样品运用LINEST函数进行焦油量回归分析的数据，获得卷烟样品的焦油回归方程，对卷烟样品烟气中焦油量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量；

在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量大于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量大于8mg/支的卷烟样品烟气中的焦油量；

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量小于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量小于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量小于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y＝m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};$ $X_2={\mathrm{TAR}}_0*\frac{P_0}{P};$ $X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y＝m1*X1+m2*X2+m3*X3+b，获得盒标焦油量小于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TAR}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数，取Δδ＝2.39；

在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标焦油量大于8mg/支的卷烟样品的焦油回归方程，校正盒标焦油量大于8mg/支的卷烟样品烟气中的焦油量的具体实现方法为：

在非标准大气压下，采集盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标焦油量大于8mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y＝m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}};$ $X_2={\mathrm{TAR}}_0*\frac{P_0}{P};$ $X_3={\left({\mathrm{TAR}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝TAR代入公式Y＝m1*X1+m2*X2+m3*X3+b，获得盒标焦油量大于8mg/支的卷烟样品焦油回归方程：

$\mathrm{TAR}={\mathrm{TAR}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{TAR}}_0})+C\sqrt{{\mathrm{TAR}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

TAR—焦油量校正指标值，mg/支；

TAR₀—焦油量实验室实测值，mg/支；

Δδ—焦油量校正常数，取Δδ＝12；

在步骤二中，在非标准大气压下，采集卷烟样品运用LINEST函数进行一氧化碳量回归分析的数据，获得卷烟样品的一氧化碳回归方程，对卷烟样品烟气中一氧化碳量进行校正时，具体包括以下两个步骤：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量；

在非标准大气压下，采集盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量大于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量大于11mg/支的卷烟样品烟气中的一氧化碳量；

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量小于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量小于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量小于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y＝m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};$ $X_2={\mathrm{CO}}_0*\frac{P_0}{P};$ $X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y＝m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量小于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ＝2.28；

在非标准大气压下，采集盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据，获得盒标一氧化碳量大于11mg/支的卷烟样品的一氧化碳回归方程，校正盒标一氧化碳量大于11mg/支的卷烟样品烟气中的一氧化碳量的具体实现方法为：

在非标准大气压下，采集盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据；

将采集的盒标一氧化碳量大于11mg/支的卷烟样品运用LINEST函数进行回归分析的数据代入多重回归公式：y＝m1*x1+m2*x2+m3*x3+b，获得m1、m2、m3、b，式中，因变量y是自变量x的函数值，m1、m2、m3是与每组x1、x2、x3值相对应的系数，b为常量；

将 $X_1={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}};$ $X_2={\mathrm{CO}}_0*\frac{P_0}{P};$ $X_3={\left({\mathrm{CO}}_0\right)}^{\frac{3}{2}}*\frac{P_0}{P};$ Y＝CO代入公式Y＝m1*X1+m2*X2+m3*X3+b，获得盒标一氧化碳量大于11mg/支的卷烟样品一氧化碳回归方程：

$\mathrm{CO}={\mathrm{CO}}_0[\frac{P_0}{P}(A-B\sqrt{{\mathrm{CO}}_0})+C\sqrt{{\mathrm{CO}}_0}]-\mathrm{\Δ\δ},$ 其中A、B、C由m1、m2、m3计算获得；

式中：P₀—标准大气压或指定大气压，kPa；

P—实验室实测大气压，kPa；

CO—一氧化碳量校正指标值，mg/支；

CO₀—一氧化碳量实验室实测值，mg/支；

Δδ—一氧化碳量校正常数，取Δδ＝47。