CN109343344A - 卷烟设备操控参数寻优方法 - Google Patents

卷烟设备操控参数寻优方法 Download PDF

Info

Publication number
CN109343344A
CN109343344A CN201811110861.XA CN201811110861A CN109343344A CN 109343344 A CN109343344 A CN 109343344A CN 201811110861 A CN201811110861 A CN 201811110861A CN 109343344 A CN109343344 A CN 109343344A
Authority
CN
China
Prior art keywords
control
data
chart
limit
subgroup
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.)
Pending
Application number
CN201811110861.XA
Other languages
English (en)
Inventor
王军
王晓东
荆素坡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Tiangong Intelligent Technology Co Ltd
Original Assignee
Beijing Tiangong Intelligent Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Tiangong Intelligent Technology Co Ltd filed Critical Beijing Tiangong Intelligent Technology Co Ltd
Priority to CN201811110861.XA priority Critical patent/CN109343344A/zh
Publication of CN109343344A publication Critical patent/CN109343344A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

卷烟设备操控参数寻优方法，采用滑动时间窗算法获取稳定度因子，当其值小于等于阈值时表明过程处于稳定工况，采用平滑滤波算法对处于稳定工况的滑动时间窗内原始数据进行算术平均，获得移动平均数据，并作为稳定工况的有效数据进行统计过程控制建模，获取控制限，绘制标准差控制图，绘制均值控制图，并分别判定过程是否稳定，获取单个参数的西格玛水平并评价过程控制能力，加权平均计算烟支质量综合过程控制能力指数，如果过程控制能力指数大于设定值，取当前批次设备参数的参数设定值作为最优工况数据，放入学习样本集；建立聚类模型；逆向标准化处理。本发明能够满足卷烟设备操控参数寻优的实际需要，为设备提供了客观、量化的参考信息。

Description

卷烟设备操控参数寻优方法
技术领域
本发明涉及数据处理技术领域，具体涉及一种卷烟设备操控参数寻 优方法。
背景技术
烟草工业企业顺应要求，正在大力落实“中国智造”，加速新一代技 术的集成实施，积极推进智能制造单元、智能车间、智能工厂和智能制 造系统升级，着力提升智能研发水平，打造柔性化、智能化、高端化、 绿色化的现代生产体系和发展模式。
卷烟机组是烟草工业企业最为重要的生产装备之一，也是卷烟生产 工艺流程的重点工序和瓶颈工序。卷烟机组的产能直接关系到卷烟生产 计划的完成情况，卷烟机组的加工水平直接影响卷烟产品的质量水平， 对产品市场营销和消费者满意度有很大的影响。卷烟机组的运行效率是 烟草行业重点关注的最为重要指标之一。
随着近年来烟草行业全球化、新的法律法规以及外部环境的重大变 化，中国烟草面临愈来愈严峻的竞争压力和社会压力。在保证产品质量 的前提下，有效提升设备的智能化水平和高效运行能力成为卷烟工厂关 注的焦点。智能化，即生产设备具有感知、分析、推理和决策等功能； 高效运行，即在规定时间内消耗更少的资源生产更多、更好的产品。卷烟工厂设备的智能化水平和高效运行能力与其自动化程度密不可分，设 备操控参数寻优技术是提升智能化水平、保证高效运行的重要手段。
发明内容
本发明实施例的目的在于提供一种卷烟设备操控参数寻优方法，面 向生产过程的工况特性，基于稳定工况识别、稳定工况SPC统计建模、 K-Means聚类等关键技术，能够满足卷烟设备操控参数寻优的实际需要。
为实现上述目的，本发明实施例提供一种卷烟设备操控参数寻优方 法，所述寻优方法包括：
1)历史运行数据的获取：选择监测变量，收集正常工况下各监测变 量的过程数据集Xn×m，其中，n为采样点个数，m为监测变量个数，建立 统计模型；
2)稳定工况的识别：基于设备运行速度，采用滑动时间窗算法获取 稳定度因子Fs，采用平滑滤波算法对处于稳定工况的滑动时间窗内原始 数据进行算术平均，获得一组移动平均数据，将移动平均数据作为稳定 工况的有效数据X′(I×m)，其中，I为有效数据的个数，m为监测变量个 数；
3)对稳定工况数据进行SPC建模：对稳定工况数据进行统计过程控 制建模，获取控制限，绘制标准差控制图，并根据标准差控制图判定工 况数据是否稳定，绘制均值控制图，并根据均值控制图判定工况数据是 否稳定，根据过程满足质量要求能力的单个度量参数即西格玛水平值评 价过程控制能力，加权平均计算烟支质量综合过程控制能力指数CPK， 如果过程控制能力指数大于设定值，取当前批次设备参数的参数设定值 作为最优工况数据，放入学习样本集Xtrain；
4)建立K-Means聚类模型：对学习样本集Xtrain进行中心化和标准 化处理，获得处理后的建模数据Xs，对建模数据Xs建立K-Means聚类模 型；
5)逆向标准化处理：设备操控参数寻优，通过获取各个聚类结果集 的最大、最小和中心样本，然后对所得样本进行逆向标准化处理，得到 最终操控参数寻优结果。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤2)中，当设 备运行速度的值小于等于阈值ζ时表明过程处于稳定工况，反之表明过 程处于过渡工况；
基于设备运行速度，采用滑动时间窗算法获取稳定度因子Fs的公式 为：
其中，VHmax表示滑动时间窗内卷烟机车速的最大值，VHmin表示滑 动时间窗内卷烟机车速的最小值，Vmax表示卷烟设备允许的车速最大值， H表示滑动时间窗口的长度。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤3)中，对稳 定工况数据进行统计过程控制建模的步骤为：
301)确定控制对象：将反映产品质量的关键质量特性作为控制对象， 控制对象包括计量型和计数型；
302)收集预备数据：以特定时间作为1个子组进行数据的汇总，按 照时间顺序，从历史数据中选择若干子组进行正态性检验。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤3)中，计量 型控制对象的数据处理方法为：
获取控制限；
获取每个子组均值标准差s和所有观测值的总均值平均标准差
获取标准差控制图的中心线和上下控制限：
式中：CLs为标准差控制图中心线；为子组标准差的平均值；UCLs为 标准差控制图上控制限；LCLs为标准差控制图下控制限；B3、B4为控制 限系数；
获取均值控制图的中心线和上下控制限：
式中：为均值控制图中心线；为子组平均值的平均值；为子组 标准差的平均值；为均值控制图上控制限；为均值控制图下控 制限；A3为控制限系数；
绘制标准差控制图，在标准差控制图上标记中心线和上下控制限， 点绘每一个子组的标准差；
观察标准差控制图所有的点是否落在控制界限内，判定过程是否处 于稳定状态，如标准差控制图显示过程处于稳定状态，则用来获取均值 图的控制限；如过程处于非稳定状态，剔除异常点，重新计算并修正控 制限直至过程处于稳定状态；
绘制均值控制图，在均值控制图上标记中心线和上下控制限，点绘 每个子组的平均值；
观察均值控制图所有的点是否落在控制界限内，并且在控制界限内 的点是否随机排列，如过程处于稳定状态，则将分析用控制图的控制限 确定为控制用控制图的控制限；如过程处于非稳定状态，剔除异常点， 重新计算并修正控制限直至过程处于稳定状态。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤3)中，计数 型控制对象的数据处理方法为：
获取控制限：
获取每个子组缺陷数和平均缺陷数，获取缺陷数控制图的中心线和 上下控制限：
式中：CLc为缺陷数控制图中心线；c为子组缺陷数的平均值；UCLc为 缺陷数控制图上控制限；UCLc为缺陷数控制图下控制限；
获取每个子组不合格品率和平均不合格品率，获取不合格品率控制 图的中心线和上下控制限：
式中：CLp为不合格品率控制图中心线；为子组不合格品率的平均 值；UCLp为不合格品率控制图上控制限；LCLp为不合格品率控制图下控 制限；
获取每个子组单位缺陷数和平均单位缺陷数，获取单位缺陷数控制 图的中心线和上下控制限：
式中：CLu为单位缺陷数控制图中心线；为子组单位缺陷数的平均 值；UCLu为单位缺陷数控制图上控制限；LCLu为单位缺陷数控制图下控 制限；
绘制缺陷数控制图，在缺陷数控制图上标记中心线和上下控制限， 点绘每一个子组的缺陷数，制成缺陷数控制图；
绘制不合格品率控制图，在不合格品率控制图上标记中心线和上下 控制限，点绘每一个子组的不合格品率，制成不合格品率控制图；
绘制单位缺陷数控制图，在单位缺陷数控制图上标记中心线和上下 控制限，点绘每一个子组的单位缺陷数，制成单位缺陷数控制图；
观察控制图所有的点是否落在控制界限内，判定过程是否处于稳定 状态，如控制图显示过程处于稳定状态，则将分析用控制图的控制限确 定为控制用控制图的控制限；如过程处于非稳定状态，剔除异常点，重 新计算并修正控制限直至过程处于稳定状态。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤3)中，当控 制图均处于稳定状态时，收集给定数量以上的子组数据获取过程西格玛 水平，评价过程控制能力指数CPK的公式为；
其中，USL为控制对象特性的规格上限，LSL为控制对象特性的规 格下限，X为所有取样数据的平均值，σ取样数据的标准差；
获取所取样数据的标准差σ，再获取规格公差T及规格中心值U；
T＝USL-LSL；U＝(USL+LSL)/2；
获取制程准确度Ca：Ca＝(X-U)/(T/2)；
获取制程精密度Cp：Cp＝T/6σ；
获取制程能力指数：Cpk＝Cp(1-|Ca|)；
加权平均获取烟支质量综合过程控制能力CPK。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤4)中，对学 习样本集Xtrain进行中心化和标准化处理，获得处理后的建模数据Xs的公 式为：
其中，xi，j′相当于Xs，i代表样本、j代表变量，表示第j个变量的 均值，sj表示第j个变量的标准差。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤4)中，对Xs建 立K-Means聚类模型的步骤为：
401)从数据集Xs随机选择k个样本作为初始的质心向量： {μ1，μ2，…，μk}；
402)对于n＝1，2，…，N，其中N为最大迭代次数
a)将簇划分C初始化为Ct＝Φt＝1，2，…，k；
b)对于i＝1，2，…，m，获取样本xi和各个质心向量μj(j＝ 1，2，…，k)的距离：将xi标记最小的为dij所对应的类 别λi，此时更新
c)对于j＝1，2，…，k，对Cj中所有的样本点重新计算新的质心
d)如果所有的k个质心向量都没有发生变化，则转到步骤3)；
403)输出簇划分C＝{C1，C2，…，Ck}。
作为卷烟设备操控参数寻优方法的优选方案，所述步骤5)中，获取 各个聚类结果集的最大、最小和中心样本，然后对所得样本进行逆向标 准化处理，逆向标准化处理的方式为：
ci×std(Xtrain)+mean(Xtrain)
其中，ci代表簇划分；
std(Xtrain)代表标准学习样本集；
mean(Xtrain)代表平均学习样本集。
本发明实施例具有如下优点：基于稳定工况识别、稳定工况SPC(统 计过程控制，Statistical Process Control)统计建模、K-Means聚类等技术， 对海量历史数据通过滑动时间窗离线计算稳定度因子，实现了卷烟机稳 定工况和过渡工况的有效识别，采用统计过程控制实现稳定工况的统计 建模以及控制限的求取，通过聚类算法获取各个聚类结果集的最大、最 小和中心样本，满足卷烟设备操控参数寻优的实际需要，能有效提高卷 烟机的可靠性和准确性，进而为设备预防维修提供了客观、量化的参考 信息。
附图说明
图1为本发明实施例提供的卷烟设备操控参数寻优方法流程图；
图2为本发明实施例提供的工况识别前的卷烟机车速图；
图3为本发明实施例提供的稳定度因子示意图；
图4为本发明实施例提供的经过平滑滤波后稳定工况的卷烟机车速 图；
图5为本发明实施例提供的稳定工况有效数据的聚类图。
具体实施方式
以下实施例用于说明本发明，但不用来限制本发明的范围。
参见图1，提供一种卷烟设备操控参数寻优方法，所述寻优方法包括：
S1：历史运行数据的获取；
S2：稳定工况的识别；
S3：对稳定工况数据进行SPC建模；
S4：建立K-Means聚类模型；
S5：逆向标准化处理。
本发明实施例具体的提供一种卷烟设备操控参数寻优方法，
S1：历史运行数据的获取：选择监测变量，监测变量包括烟支质量 指标和设备操控参数，收集正常工况下各监测变量的过程数据集Xn×m， 其中，n为采样点个数，m为监测变量个数，建立统计模型；
S2：稳定工况的识别：基于设备运行速度，采用滑动时间窗算法获 取稳定度因子Fs，当稳定度因子Fs值小于等于阈值ζ时表明过程处于稳定 工况，反之表明过程处于过渡工况，采用平滑滤波算法对处于稳定工况 的滑动时间窗内原始数据进行算术平均，获得一组移动平均数据，将移 动平均数据作为稳定工况的有效数据X′(I×m)，其中，I为有效数据的个 数，m为监测变量个数；
S3：对稳定工况数据进行SPC建模：对稳定工况数据进行统计过程 控制建模，获取控制限，绘制标准差控制图，并根据标准差控制图判定 工况数据是否稳定，绘制均值控制图，并根据均值控制图判定工况数据 是否稳定，根据过程满足质量要求能力的单个度量参数西格玛水平评价 过程控制能力，加权平均计算烟支质量综合过程控制能力指数CPK，如果烟支质量综合综合过程控制能力指数CPK>1.6，则认为当前批次为优 秀批次，取当前批次设备参数的参数设定值作为最优工况数据，放入学 习样本集Xtrain；
S4：建立K-Means聚类模型：对学习样本集Xtrain进行中心化和标准 化处理，获得处理后的建模数据Xs，对建模数据Xs建立K-Means聚类模 型；
S5：逆向标准化处理：设备操控参数寻优，通过获取各个聚类结果 集的最大、最小和中心样本，然后对所得样本进行逆向标准化处理，得 到最终操控参数寻优结果。
实施例中，步骤S3判定工况数据是否稳定的标准如下表：
步骤S3中，西格玛水平是将过程输出的平均值、标准差与质量要求 的目标值、规格限联系起来进行比较，是对过程满足质量要求能力的一 种度量。西格玛水平越高，过程满足质量要求的能力越高；反之，西格 玛水平越低，过程满足质量要求的能力越低。
6个西格玛＝3.4失误/百万机会，意味着卓越的管理，强大的竞争力 和忠诚的客户。
5个西格玛＝230失误/百万机会，意味着优秀的管理、很强的竞争力 和比较忠诚的客户。
4个西格玛＝6,210失误/百万机会，意味着较好的管理和运营能力， 满意的客户。
3个西格玛＝66,800失误/百万机会，意味着平平常常的管理，缺乏 竞争力。
2个西格玛＝308,000失误/百万机会，意味着企业资源每天都有三分 之一的浪费。
1个西格玛＝690,000失误/百万机会，意味着每天有三分之二的事情 做错，企业无法生存。
卷烟设备操控参数寻优方法的一个实施例中，所述S2中，当设备运 行速度的值小于等于阈值ζ时表明过程处于稳定工况，反之表明过程处 于过渡工况；
基于设备运行速度，采用滑动时间窗算法获取稳定度因子Fs的公式 为：
其中，VHmax表示滑动时间窗内卷烟机车速的最大值，VHmin表示滑 动时间窗内卷烟机车速的最小值，Vmax表示卷烟设备允许的车速最大值， H表示滑动时间窗口的长度。
所述S3中，对稳定工况数据进行统计过程控制建模的步骤为：
S301：确定控制对象：将反映产品质量的关键质量特性作为控制对 象，控制对象包括计量型和计数型。其中，质量、圆周、吸阻、通风率、 空头、漏气，其中质量、圆周、吸阻、通风率为计量型，空头、漏气为 计数型。
S302：收集预备数据：以特定时间作为1个子组进行数据的汇总， 按照时间顺序，从历史数据中选择若干子组进行正态性检验。本实施例 中每5分钟作为1个子组进行数据的汇总(随机取样10支)，按照时间 顺序，从近期历史数据(一个月为宜)中选择30个子组，进行正态性检验， 数据应符合正态分布。收集预备数据前应跟踪每一组数据，有异常数据时应及时标记，以备后续数据异常剔除分析用。
所述S3中，计量型控制对象的数据处理方法为：
获取控制限；
获取每个子组均值标准差s和所有观测值的总均值平均标准差
获取标准差控制图的中心线和上下控制限：
式中：CLs为标准差控制图中心线；为子组标准差的平均值；UCLs为 标准差控制图上控制限；LCLs为标准差控制图下控制限；B3、B4为控制 限系数；
获取均值控制图的中心线和上下控制限：
式中：为均值控制图中心线；为子组平均值的平均值；为子组 标准差的平均值；为均值控制图上控制限；为均值控制图下控 制限；A3为控制限系数；
绘制标准差控制图，用左侧纵坐标表示各子组标准差，用横坐标表 示子组号；在标准差控制图上标记中心线和上下控制限，点绘每一个子 组的标准差；
观察标准差控制图所有的点是否落在控制界限内，并依据A类判稳 准则判定过程是否处于稳定状态。如标准差控制图显示过程处于稳定状 态，则可用来计算均值图的控制限；如过程处于非稳定状态，应剔除异 常点，重新计算并修正控制限直至过程处于稳定状态，同时应查明异常 点产生的原因，并采取适当的补救措施防止其再次出现。如异常点的数 量超过子组数量的20％，应重新收集数据。
绘制均值控制图，用左侧纵坐标表示各子组均值，用横坐标表示子 组号；在均值控制图上标记中心线和上下控制限，点绘每个子组的平均 值；
观察均值控制图所有的点是否落在控制界限内，并且在控制界限内 的点是否随机排列，依据A类和B类判稳准则判定过程是否处于稳定状 态。如过程处于稳定状态，则将分析用控制图的控制限确定为控制用控 制图的控制限；如过程处于非稳定状态，应剔除异常点，重新计算并修 正控制限直至过程处于稳定状态，同时应查明异常点产生的原因，并采 取适当的补救措施防止其再次出现。如异常点的数量超过子组数量的 20％，应重新收集数据。
所述S3中，计数型控制对象的数据处理方法为：
获取控制限：
获取每个子组缺陷数和平均缺陷数，获取缺陷数控制图的中心线和 上下控制限：
式中：CLc为缺陷数控制图中心线；c为子组缺陷数的平均值；UCLc为 缺陷数控制图上控制限；UCLc为缺陷数控制图下控制限；
获取每个子组不合格品率和平均不合格品率，获取不合格品率控制 图的中心线和上下控制限：
式中：CLp为不合格品率控制图中心线；为子组不合格品率的平均 值；UCLp为不合格品率控制图上控制限；LCLp为不合格品率控制图下控 制限；
获取每个子组单位缺陷数和平均单位缺陷数，获取单位缺陷数控制 图的中心线和上下控制限：
式中：CLu为单位缺陷数控制图中心线；为子组单位缺陷数的平均 值；UCLu为单位缺陷数控制图上控制限；LCLu为单位缺陷数控制图下控 制限；由于控制限不可能为负，如计数下控制限计算为负时不应标出。
绘制缺陷数控制图，用左侧纵坐标表示各子组缺陷数，用横坐标表 示子组号；在缺陷数控制图上标记中心线和上下控制限，点绘每一个子 组的缺陷数，制成缺陷数控制图；
绘制不合格品率控制图，用左侧纵坐标表示各子组不合格品率，用 横坐标表示子组号；在不合格品率控制图上标记中心线和上下控制限， 点绘每一个子组的不合格品率，制成不合格品率控制图；
绘制单位缺陷数控制图，用左侧纵坐标表示各子组单位缺陷数，用 横坐标表示子组号；在单位缺陷数控制图上标记中心线和上下控制限， 点绘每一个子组的单位缺陷数，制成单位缺陷数控制图；
观察控制图所有的点是否落在控制界限内，并依据A类判稳准则判 定过程是否处于稳定状态。如控制图显示过程处于稳定状态，则将分析 用控制图的控制限确定为控制用控制图的控制限；如过程处于非稳定状 态，应剔除异常点，重新计算并修正控制限直至过程处于稳定状态，同 时应查明异常点产生的原因，并采取适当的补救措施防止其再次出现。 如异常点的数量超过子组数量的20％，应重新收集数据。
卷烟设备操控参数寻优方法的一个实施例中，所述S3中，当控制图 均处于稳定状态时，收集给定数量以上的子组数据获取过程西格玛水平， 评价过程控制能力指数CPK的公式为；
其中，USL为控制对象特性的规格上限，LSL为控制对象特性的规 格下限，X为所有取样数据的平均值，σ取样数据的标准差；
获取所取样数据的标准差σ，再获取规格公差T及规格中心值U；
T＝USL-LSL；U＝(USL+LSL)/2；
获取制程准确度Ca：Ca＝(X-U)/(T/2)；
获取制程精密度Cp：Cp＝T/6σ；
获取制程能力指数：Cpk＝Cp(1-|Ca|)；
加权平均获取烟支质量综合过程控制能力CPK。
卷烟设备操控参数寻优方法的一个实施例中，所述S4中，对学习样 本集Xtrain进行中心化和标准化处理，获得处理后的建模数据Xs的公式为：
其中，xi，j′相当于Xs，i代表样本、j代表变量，表示第j个变量的 均值，sj表示第j个变量的标准差。
所述S4中，对Xs建立K-Means聚类模型的步骤为：
S401：从数据集Xs随机选择k个样本作为初始的质心向量： {μ1，μ2，…，μk}；
S402：对于n＝1，2，…，N，其中N为最大迭代次数
a)将簇划分C初始化为Ct＝Φt＝1，2，…，k；
b)对于i＝1，2，…，m，获取样本xi和各个质心向量μj(j＝ 1，2，…，k)的距离：将xi标记最小的为dij所对应的类 别λi，此时更新
c)对于j＝1，2，…，k，对Cj中所有的样本点重新计算新的质心
d)如果所有的k个质心向量都没有发生变化，则转到步骤3)；
S403：输出簇划分C＝{C1，C2，…，Ck}。
具体的，S5中，获取各个聚类结果集的最大、最小和中心样本，然 后对所得样本进行逆向标准化处理，逆向标准化处理的方式为：
ci×std(Xtrain)+mean(Xtrain)
其中，ci代表簇划分；
std(Xtrain)代表标准学习样本集；
mean(Xtrain)代表平均学习样本集。
具体的，选取ZJ17E机台的卷包机，采集底层PLC中的卷烟机运行 数据，包含以下变量：针辊供丝补偿量、针辊停止、针辊供丝量系数、 风室负压、圆周修正、圆周控制电机、启动位置、剔出极限、废品极限、 重量校正、内部目标重量、起始位置、过轻烟支重量极限、过重烟支重 量极限、软点重量极限、硬点重量极限、松头重量极限、烟丝端部大小、 滤嘴端大小、压实端位置、修整器位置、烙铁1温度、烙铁2温度、胶 枪温度、搓板温度、水松纸温度、梗签分离压力、小风机正压、前轨密 度、后轨密度、管道风速、紧密度、卷烟机车速、梗签剔除量、弹丝辘 转速、回丝量、重量、吸丝带位置、吸阻、通风度、漏气、松头、光学 外观、圆周值等43个监测变量。监测变量的采样频率为1秒/次，获得的 原始数据为D(490542x44)。
采用本实施例的技术方案，滑动时间窗的长度H取为10，允许的车 速最大值Vmax取为8000，阈值ζ取为0.02，获得的有效数据为D(65314 X 44)。以前1000个原始采样数据为例，工况识别前卷烟机车速如图2 所示，相应的稳定度因子如图3所示，经过平滑滤波后稳定工况的卷烟 机车速如图4所示，最优工况数据的聚类图如图5所示。实现了卷烟机 稳定工况和过渡工况的有效识别，采用统计过程控制实现稳定工况的统 计建模以及控制限的求取，通过聚类算法获取各个聚类结果集的最大、 最小和中心样本,满足卷烟设备操控参数寻优的实际需要，能有效提高卷 烟机的可靠性和准确性，进而为设备预防维修提供了客观、量化的参考 信息。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的 描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技 术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的 这些修改或改进，均属于本发明要求保护的范围。

Claims (9)

1.卷烟设备操控参数寻优方法，其特征在于：所述寻优方法包括：
1)历史运行数据的获取：选择监测变量，收集正常工况下各监测变量的过程数据集Xn×m，其中，n为采样点个数，m为监测变量个数，建立统计模型；
2)稳定工况的识别：基于设备运行速度，采用滑动时间窗算法获取稳定度因子Fs，当其值小于等于阈值ζ时表明过程处于稳定工况，采用平滑滤波算法对处于稳定工况的滑动时间窗内原始数据进行算术平均，获得一组移动平均数据，将移动平均数据作为稳定工况的有效数据X′(I×m)，其中，I为有效数据的个数，m为监测变量个数；
3)对稳定工况数据进行SPC建模：对稳定工况数据进行统计过程控制建模，获取控制限，绘制标准差控制图，并根据标准差控制图判定工况数据是否稳定，绘制均值控制图，并根据均值控制图判定工况数据是否稳定，根据过程满足质量要求能力的单个度量参数即西格玛水平值评价过程控制能力，加权平均计算烟支质量综合过程控制能力指数CPK，如果过程控制能力指数大于设定值，取当前批次设备参数的参数设定值作为最优工况数据，放入学习样本集Xtrain；
4)建立K-Means聚类模型：对学习样本集Xtrain进行中心化和标准化处理，获得处理后的建模数据Xs，对建模数据Xs建立K-Means聚类模型；
5)逆向标准化处理：设备操控参数寻优，通过获取各个聚类结果集的最大、最小和中心样本，然后对所得样本进行逆向标准化处理，得到最终操控参数寻优结果。
2.根据权利要求1所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤2)中，当设备运行速度的值小于等于阈值ζ时表明过程处于稳定工况，反之表明过程处于过渡工况；
基于设备运行速度，采用滑动时间窗算法获取稳定度因子Fs的公式为：
其中，VHmax表示滑动时间窗内卷烟机车速的最大值，VHmin表示滑动时间窗内卷烟机车速的最小值，Vmax表示卷烟设备允许的车速最大值，H表示滑动时间窗口的长度。
3.根据权利要求1所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤3)中，对稳定工况数据进行统计过程控制建模的步骤为：
301)确定控制对象：将反映产品质量的关键质量特性作为控制对象，控制对象包括计量型和计数型；
302)收集预备数据：以特定时间作为1个子组进行数据的汇总，按照时间顺序，从历史数据中选择若干子组进行正态性检验。
4.根据权利要求3所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤3)中，计量型控制对象的数据处理方法为：
获取控制限；
获取每个子组均值标准差s和所有观测值的总均值平均标准差
获取标准差控制图的中心线和上下控制限：
式中：CLs为标准差控制图中心线；为子组标准差的平均值；UCLs为标准差控制图上控制限；LCLs为标准差控制图下控制限；B3、B4为控制限系数；
获取均值控制图的中心线和上下控制限：
式中：为均值控制图中心线；为子组平均值的平均值；为子组标准差的平均值；为均值控制图上控制限；为均值控制图下控制限；A3为控制限系数；
绘制标准差控制图，在标准差控制图上标记中心线和上下控制限，点绘每一个子组的标准差；
观察标准差控制图所有的点是否落在控制界限内，判定过程是否处于稳定状态，如标准差控制图显示过程处于稳定状态，则用来获取均值图的控制限；如过程处于非稳定状态，剔除异常点，重新计算并修正控制限直至过程处于稳定状态；
绘制均值控制图，在均值控制图上标记中心线和上下控制限，点绘每个子组的平均值；
观察均值控制图所有的点是否落在控制界限内，并且在控制界限内的点是否随机排列，如过程处于稳定状态，则将分析用控制图的控制限确定为控制用控制图的控制限；如过程处于非稳定状态，剔除异常点，重新计算并修正控制限直至过程处于稳定状态。
5.根据权利要求3所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤3)中，计数型控制对象的数据处理方法为：
获取控制限：
获取每个子组缺陷数和平均缺陷数，获取缺陷数控制图的中心线和上下控制限：
式中：CLc为缺陷数控制图中心线；c为子组缺陷数的平均值；UCLc为缺陷数控制图上控制限；UCLc为缺陷数控制图下控制限；
获取每个子组不合格品率和平均不合格品率，获取不合格品率控制图的中心线和上下控制限：
式中：CLp为不合格品率控制图中心线；为子组不合格品率的平均值；UCLp为不合格品率控制图上控制限；LCLp为不合格品率控制图下控制限；
获取每个子组单位缺陷数和平均单位缺陷数，获取单位缺陷数控制图的中心线和上下控制限：
式中：CLu为单位缺陷数控制图中心线；为子组单位缺陷数的平均值；UCLu为单位缺陷数控制图上控制限；LCLu为单位缺陷数控制图下控制限；
绘制缺陷数控制图，在缺陷数控制图上标记中心线和上下控制限，点绘每一个子组的缺陷数，制成缺陷数控制图；
绘制不合格品率控制图，在不合格品率控制图上标记中心线和上下控制限，点绘每一个子组的不合格品率，制成不合格品率控制图；
绘制单位缺陷数控制图，在单位缺陷数控制图上标记中心线和上下控制限，点绘每一个子组的单位缺陷数，制成单位缺陷数控制图；
观察控制图所有的点是否落在控制界限内，判定过程是否处于稳定状态，如控制图显示过程处于稳定状态，则将分析用控制图的控制限确定为控制用控制图的控制限；如过程处于非稳定状态，剔除异常点，重新计算并修正控制限直至过程处于稳定状态。
6.根据权利要求1或3所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤3)中，当控制图均处于稳定状态时，收集给定数量以上的子组数据获取过程西格玛水平，评价过程控制能力指数CPK的公式为；
其中，USL为控制对象特性的规格上限，LSL为控制对象特性的规格下限，X为所有取样数据的平均值，σ取样数据的标准差；
获取所取样数据的标准差σ，再获取规格公差T及规格中心值U；
T＝USL-LSL；U＝(USL+LSL)/2；
获取制程准确度Ca：Ca＝(X-U)/(T/2)；
获取制程精密度Cp：Cp＝T/6σ；
获取制程能力指数：Cpk＝Cp(1-|Ca|)；
加权平均获取烟支质量综合过程控制能力CPK。
7.根据权利要求1所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤4)中，对学习样本集Xtrain进行中心化和标准化处理，获得处理后的建模数据Xs的公式为：
其中，xi，j′相当于Xs，i代表样本、j代表变量，表示第j个变量的均值，sj表示第j个变量的标准差。
8.根据权利要求1所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤4)中，对Xs建立K-Means聚类模型的步骤为：
401)从数据集Xs随机选择k个样本作为初始的质心向量：{μ1，μ2，…，μk}；
402)对于n＝1，2，…，N，其中N为最大迭代次数
a)将簇划分C初始化为Ct＝Φt＝1，2，…，k；
b)对于i＝1，2，…，m，获取样本xi和各个质心向量μj(j＝1，2，…，k)的距离：将xi标记最小的为dij所对应的类别λi，此时更新
c)对于j＝1，2，…，k，对Cj中所有的样本点重新计算新的质心
d)如果所有的k个质心向量都没有发生变化，则转到步骤3)；
403)输出簇划分C＝{C1，C2，…，Ck}。
9.根据权利要求1所述的卷烟设备操控参数寻优方法，其特征在于：所述步骤5)中，获取各个聚类结果集的最大、最小和中心样本，然后对所得样本进行逆向标准化处理，逆向标准化处理的方式为：
ci×std(Xtrain)+mean(Xtrain)
其中，ci代表簇划分；
std(Xtrain)代表标准学习样本集；
mean(Xtrain)代表平均学习样本集。
CN201811110861.XA 2018-09-21 2018-09-21 卷烟设备操控参数寻优方法 Pending CN109343344A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811110861.XA CN109343344A (zh) 2018-09-21 2018-09-21 卷烟设备操控参数寻优方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811110861.XA CN109343344A (zh) 2018-09-21 2018-09-21 卷烟设备操控参数寻优方法

Publications (1)

Publication Number Publication Date
CN109343344A true CN109343344A (zh) 2019-02-15

Family

ID=65306125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811110861.XA Pending CN109343344A (zh) 2018-09-21 2018-09-21 卷烟设备操控参数寻优方法

Country Status (1)

Country Link
CN (1) CN109343344A (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110442092A (zh) * 2019-07-05 2019-11-12 柳州钢铁股份有限公司 一种基于热连轧生产的自动统计和分析方法
CN110794787A (zh) * 2019-11-18 2020-02-14 马鞍山钢铁股份有限公司 一种适合多品种混合生产的spc方法
CN111184251A (zh) * 2019-12-20 2020-05-22 深圳市鸿捷源自动化系统有限公司 一种诊断吸丝带故障的处理方法
CN111667139A (zh) * 2020-04-22 2020-09-15 广东中烟工业有限责任公司 一种烟丝制造质量检测预警模块及方法
CN112132722A (zh) * 2020-08-20 2020-12-25 彭涛 政务热线数量趋势异常确定方法、装置、电子设备及介质
CN112420132A (zh) * 2020-10-29 2021-02-26 重庆大学 一种汽油催化裂化过程产品质量优化控制方法
CN112434867A (zh) * 2020-11-30 2021-03-02 红云红河烟草(集团)有限责任公司 一种叶片段水份智能预测模型及应用
CN113139078A (zh) * 2021-03-31 2021-07-20 青岛奥利普自动化控制系统有限公司 控制图生成方法及电子设备
CN113325796A (zh) * 2021-06-09 2021-08-31 内蒙古昆明卷烟有限责任公司 一种包装机监测方法、设备及介质
CN113592314A (zh) * 2021-08-02 2021-11-02 红云红河烟草(集团)有限责任公司 基于西格玛水平的制丝工艺质量评价方法
CN114324764A (zh) * 2021-12-30 2022-04-12 红云红河烟草(集团)有限责任公司 基于物料特性值的卷烟制丝过程制品质量多维评测方法
CN115060351A (zh) * 2022-06-20 2022-09-16 河南中烟工业有限责任公司 一种电子秤在线监测方法、系统、电子设备及计算机可读存储介质
CN115659535A (zh) * 2022-10-25 2023-01-31 苏州玖钧智能装备有限公司 一种单轴龙门机械手灌胶机的设备评估方法及系统
CN115857455A (zh) * 2022-12-19 2023-03-28 河南中烟工业有限责任公司 一种卷烟生产过程物测指标过程质量控制方法
CN115877805A (zh) * 2022-12-07 2023-03-31 河南中烟工业有限责任公司 一种基于烟支圆周数据的质量控制方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005173911A (ja) * 2003-12-10 2005-06-30 Trecenti Technologies Inc 工程管理システムおよび工程管理方法
CN101414183A (zh) * 2008-12-03 2009-04-22 山东中烟工业公司 基于灰色关联分析的卷烟工序质量综合评价系统及其方法
CN104317285A (zh) * 2014-11-10 2015-01-28 广东中烟工业有限责任公司 一种卷烟生产过程中异常原因确定方法及装置
CN104360677A (zh) * 2014-12-18 2015-02-18 厦门烟草工业有限责任公司 一种卷烟加工过程质量评价与诊断方法
CN104503402A (zh) * 2014-12-13 2015-04-08 中国烟草总公司郑州烟草研究院 一种卷烟加工烟支卷制质量稳定性的检验方法
CN105118332A (zh) * 2015-09-30 2015-12-02 中国民用航空总局第二研究所 一种基于聚类分析法的空中交通管制模拟仿真异常检测方法和装置
CN105573290A (zh) * 2015-12-16 2016-05-11 浙江中烟工业有限责任公司 卷烟工厂超高速条盒包装机多工况过程在线监测与故障诊断方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005173911A (ja) * 2003-12-10 2005-06-30 Trecenti Technologies Inc 工程管理システムおよび工程管理方法
CN101414183A (zh) * 2008-12-03 2009-04-22 山东中烟工业公司 基于灰色关联分析的卷烟工序质量综合评价系统及其方法
CN104317285A (zh) * 2014-11-10 2015-01-28 广东中烟工业有限责任公司 一种卷烟生产过程中异常原因确定方法及装置
CN104503402A (zh) * 2014-12-13 2015-04-08 中国烟草总公司郑州烟草研究院 一种卷烟加工烟支卷制质量稳定性的检验方法
CN104360677A (zh) * 2014-12-18 2015-02-18 厦门烟草工业有限责任公司 一种卷烟加工过程质量评价与诊断方法
CN105118332A (zh) * 2015-09-30 2015-12-02 中国民用航空总局第二研究所 一种基于聚类分析法的空中交通管制模拟仿真异常检测方法和装置
CN105573290A (zh) * 2015-12-16 2016-05-11 浙江中烟工业有限责任公司 卷烟工厂超高速条盒包装机多工况过程在线监测与故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
夏国华: "《汽车制造业实施ISO/TS16949:2002实践》", 31 July 2004 *
蔡慧丽: "零陵卷烟厂基于统计过程控制SPC系统的质量管理体系研究", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110442092A (zh) * 2019-07-05 2019-11-12 柳州钢铁股份有限公司 一种基于热连轧生产的自动统计和分析方法
CN110794787A (zh) * 2019-11-18 2020-02-14 马鞍山钢铁股份有限公司 一种适合多品种混合生产的spc方法
CN111184251A (zh) * 2019-12-20 2020-05-22 深圳市鸿捷源自动化系统有限公司 一种诊断吸丝带故障的处理方法
CN111667139A (zh) * 2020-04-22 2020-09-15 广东中烟工业有限责任公司 一种烟丝制造质量检测预警模块及方法
CN112132722A (zh) * 2020-08-20 2020-12-25 彭涛 政务热线数量趋势异常确定方法、装置、电子设备及介质
CN112132722B (zh) * 2020-08-20 2023-12-26 彭涛 政务热线数量趋势异常确定方法、装置、电子设备及介质
CN112420132A (zh) * 2020-10-29 2021-02-26 重庆大学 一种汽油催化裂化过程产品质量优化控制方法
CN112434867A (zh) * 2020-11-30 2021-03-02 红云红河烟草(集团)有限责任公司 一种叶片段水份智能预测模型及应用
CN113139078A (zh) * 2021-03-31 2021-07-20 青岛奥利普自动化控制系统有限公司 控制图生成方法及电子设备
CN113139078B (zh) * 2021-03-31 2022-11-29 青岛奥利普奇智智能工业技术有限公司 控制图生成方法及电子设备
CN113325796A (zh) * 2021-06-09 2021-08-31 内蒙古昆明卷烟有限责任公司 一种包装机监测方法、设备及介质
CN113592314A (zh) * 2021-08-02 2021-11-02 红云红河烟草(集团)有限责任公司 基于西格玛水平的制丝工艺质量评价方法
CN113592314B (zh) * 2021-08-02 2024-03-26 红云红河烟草(集团)有限责任公司 基于西格玛水平的制丝工艺质量评价方法
CN114324764A (zh) * 2021-12-30 2022-04-12 红云红河烟草(集团)有限责任公司 基于物料特性值的卷烟制丝过程制品质量多维评测方法
CN115060351A (zh) * 2022-06-20 2022-09-16 河南中烟工业有限责任公司 一种电子秤在线监测方法、系统、电子设备及计算机可读存储介质
CN115659535B (zh) * 2022-10-25 2023-10-13 苏州玖钧智能装备有限公司 一种单轴龙门机械手灌胶机的设备评估方法及系统
CN115659535A (zh) * 2022-10-25 2023-01-31 苏州玖钧智能装备有限公司 一种单轴龙门机械手灌胶机的设备评估方法及系统
CN115877805A (zh) * 2022-12-07 2023-03-31 河南中烟工业有限责任公司 一种基于烟支圆周数据的质量控制方法
CN115857455A (zh) * 2022-12-19 2023-03-28 河南中烟工业有限责任公司 一种卷烟生产过程物测指标过程质量控制方法
CN115857455B (zh) * 2022-12-19 2024-10-22 河南中烟工业有限责任公司 一种卷烟生产过程物测指标过程质量控制方法

Similar Documents

Publication Publication Date Title
CN109343344A (zh) 卷烟设备操控参数寻优方法
CN109670676A (zh) 基于支持向量数据描述的配网台区风险预警方法及系统
CN101516099B (zh) 一种传感器网络异常检测方法
CN105467975A (zh) 一种设备故障诊断方法
CN108304567B (zh) 高压变压器工况模式识别与数据分类方法及系统
CN110555058A (zh) 基于改进决策树的电力通信设备状态的预测方法
CN116148753A (zh) 一种智能电能表运行误差监测系统
CN116485020B (zh) 一种基于大数据的供应链风险识别预警方法、系统及介质
CN105974356A (zh) 一种电力计量自动化检定流水线故障诊断方法
CN111882188A (zh) 基于Birch聚类算法的工艺质量同质化水平评价方法及系统
CN110555619A (zh) 一种基于智能配电网的供电能力评估方法
CN115510302A (zh) 基于大数据统计的智能工厂数据分类方法
CN117056688A (zh) 一种基于数据分析的新材料生产数据管理系统及方法
CN117171157A (zh) 基于数据分析的清算数据采集清洗方法
CN117933822A (zh) 一种基于边缘计算的物联网数据管理系统及方法
CN109960780A (zh) 一种机泵运行经济性能和稳定状态的评判系统及评判方法
CN113726558A (zh) 基于随机森林算法的网络设备流量预测系统
CN107274025B (zh) 一种实现用电模式智能识别与管理的系统和方法
CN116975670A (zh) 一种短期调荷场景下线损特征分类处理方法
CN115392663A (zh) 一种基于大数据的数据采集处理方法
CN115658772A (zh) 一种无人机光伏巡检数据类资产管理方法及系统
CN115392710A (zh) 一种基于数据过滤的风电机组运行决策方法及系统
CN110175705B (zh) 一种负荷预测方法及包含该方法的存储器、系统
CN108335231A (zh) 一种自动匹配的配电网数据诊断方法
CN115879680A (zh) 钢铁表面缺陷判定规则管理系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190215

RJ01 Rejection of invention patent application after publication