CN110826229B - 基于长期预测的烘丝干头过程建模和优化设定控制方法 - Google Patents
基于长期预测的烘丝干头过程建模和优化设定控制方法 Download PDFInfo
- Publication number
- CN110826229B CN110826229B CN201911083453.4A CN201911083453A CN110826229B CN 110826229 B CN110826229 B CN 110826229B CN 201911083453 A CN201911083453 A CN 201911083453A CN 110826229 B CN110826229 B CN 110826229B
- Authority
- CN
- China
- Prior art keywords
- drying
- tps
- tobacco
- follows
- model
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A24—TOBACCO; CIGARS; CIGARETTES; SIMULATED SMOKING DEVICES; SMOKERS' REQUISITES
- A24B—MANUFACTURE OR PREPARATION OF TOBACCO FOR SMOKING OR CHEWING; TOBACCO; SNUFF
- A24B9/00—Control of the moisture content of tobacco products, e.g. cigars, cigarettes, pipe tobacco
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Manufacture Of Tobacco Products (AREA)
- Drying Of Solid Materials (AREA)
Abstract
本发明公开了一种基于长期预测的烘丝干头过程建模和优化设定控制方法。该方法,首先利用不同工况下采集的多段烘丝干头过程历史数据进行TPS‑ARX建模。在TPS‑ARX模型的参数优化过程中,通过合理定义目标函数,得到了可使TPS‑ARX模型长期预测误差最小的模型参数。并设计了基于多段烘丝干头历史数据的、可使模型预测输出与设定值之间累积误差最小的、各输入变量最优设定曲线参数优化方法。该方法优化出的烘丝干头过程最优设定曲线,可使出口烟丝含水率尽快上升至设定值附近。
Description
技术领域
本发明涉及烘丝机出口水分自动控制方法,尤其涉及一种基于长期预测的烘丝机烘丝干头过程的建模和优化设定控制方法。
背景技术
按照滚筒式烘丝机烟草烘丝生产的工艺要求,整个烟草烘丝过程分为3个阶段:干头、中间和干尾。在烘丝机启动后,当检测到入口烟丝流量时,烘丝干头过程开始;当检测到的出口烟丝含水率基本稳定在设定值时,烘丝干头过程结束。然后,进入一段相对较长的中间过程。当入口烟丝流量由正常值变为零时,标志着烘丝干尾过程开始。在烘丝干头过程中,入口烟丝流量不断增加,但并无出口烟丝含水率的检测值,因此无法进行反馈控制器设计,极易造成烟草烘丝干头过程的出口烟丝含水率控制品质差、干料多,从而造成烟丝原料的浪费,严重影响企业生产的效率。目前,“干头”问题也是烘丝过程出口烟丝含水率控制的难点所在。
考虑到滚筒式烘丝机干头过程的特殊性:缺少出口烟丝含水率的检测值,采用机理建模无法获取该过程的精确数学模型。目前,已有技术主要采用基于数据驱动的辨识模型来建立该过程的数学模型。专利“201210376264.8”提出了一种“基于模型预测的烘丝机出口水分控制方法”,该方法采用一种基于高斯基函数的状态相依模型来描述烘丝过程的干头过程。专利“201310659839.1”提出了“一种烘丝机头尾段工艺变量优化控制方法”,该方法采用一种基于三次函数的状态相依模型来描述烘丝过程的干头过程。但上述方法均为基于某一次烟草烘丝过程的运行历史数据建立的干头过程一步向前预测模型,其均是基于辨识出的一步向前预测模型设计的干头过程设定控制方法。而在实际的烟草烘丝过程中,出口烟丝含水率极易受到外界气温、湿度和工况条件等的影响。考虑到实际烘丝干头过程的特殊性,如何充分利用不同工况下采集的多段输入输出数据来建立干头过程的长期预测模型,并基于该长期预测模型进行控制器设计仍是亟待解决的关键技术难题。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种基于长期预测的烘丝干头过程建模和优化设定控制方法。该方法,首先基于不同工况下的多段烘丝干头过程采样数据建立基于薄板样条函数的状态相依自回归(TPS-ARX)模型的干头过程长期预测模型,然后基于该长期预测模型设计干头过程的优化设定控制方法。
为解决上述技术问题,本发明所采用的技术方案是:基于长期预测的烘丝干头过程建模和优化设定控制方法,包括以下步骤:
1)建立用于对烟草烘丝干头过程进行预测的TPS-ARX模型结构为:
其中:yt为烘丝干头过程t时刻的出口烟丝含水率,为烘丝干头过程t时刻的排潮风门开度,/>为t时刻的热风风温,/>为t时刻的入口烟丝流量,/>为t时刻的入口烟丝含水率,ξt为t时刻的建模误差;xt-1为TPS-ARX模型的状态量,且/>φ0(xt-1)和均为关于xt-1的薄板样条(ThinPlate Spline,TPS)函数型系数,且具体结构为:/> 上式中,/>为TPS-ARX模型的非线性参数集;为TPS-ARX模型的线性参数集。
3)优化步骤2)中建立的烘丝干头过程长期预测TPS-ARX模型的参数集θL和θN。
4)基于步骤3)中建立烘丝干头过程TPS-ARX模型,通过选择待优化的控制变量函数曲线,设计烘丝干头过程的最优设定控制方法。
步骤3)中,优化步骤2)中建立的TPS-ARX模型的参数集θL和θN,具体方法如下:
(1)定义模型参数优化的目标函数如下:
上式中,为第β组烘丝干头过程采样数据集对应的TPS-ARX模型n步向前预测的建模误差向量;/>是采集的第β组烘丝干头过程数据集,在第τ+p时刻的出口烟丝含水率实际值,且τ=6,n为该数据集的长度,且β=1,2,…,α;/>是采集的第β组烘丝干头过程数据集,在第τ时刻的基于步骤2)中TPS-ARX模型的p步向前预测的输出值。
(2)采用高斯牛顿法求解如下非线性优化问题:
步骤4)中,烘丝干头过程最优设定控制方法的设计步骤如下:
其中,{x1,x2,x3,x4},{x5,x6,x7,x8}和{z1,z2,z3,z4}分别为函数f1(t),f2(t)和f3(t)中的待优化参数集,并定义x={x1,x2,x3,x4,x5,x6,x7,x8},z={z1,z2,z3,z4}。则参数集x和z的具体优化方法如下:
Step 1,将上述控制变量的设定函数曲线代入已建立的干头过程TPS-ARX长期预测模型,可得在该输入曲线情况下的模型预测输出为:则基于采集的第β组烘丝干头过程数据序列,得到的相应TPS-ARX模型的预测输出为:
Step 2,定义优化参数集x和z的目标函数如下:
Step 3,采用高斯牛顿法求解如下非线性优化问题:
与现有技术相比,本发明所具有的有益效果为:考虑到在实际的烟草烘丝干头过程中,出口烟丝含水率极易受到外界气温、湿度和工况条件等的影响。本发明充分利用不同工况下采集的多段历史数据对烘丝干头过程进行TPS-ARX建模。并通过合理定义目标函数,得到了可使TPS-ARX模型长期预测误差最小的模型参数。该方法有效解决了已有技术中仅仅采用某次工况下采集历史建立的干头过程一步向前预测模型,存在的模型泛化性差,且未考虑模型长期预测性能的问题。并设计了基于多段烘丝干头历史数据的,可使模型预测输出与设定值之间累积误差最小的各输入变量最优设定曲线参数优化方法。该方法可使出口烟丝含水率尽快上升至设定值附近,可有效解决烘丝干头过程中经常出现的干料过多的问题。
附图说明
图1为本发明涉及的烟草烘丝干头过程工艺示意图。
具体实施方式
本发明所针对的滚筒式烘丝机烘丝干头过程工艺示意图如图1所示。首先,待烘丝的烟丝经传送带进入滚筒,经烘丝加工后传出滚筒。图1中,1为排潮风门开度检测;2为热风风温检测点;3为入口烟丝流量检测点;4入口烟丝含水率检测点;5为出口烟丝含水率检测点。当3处检测到入口烟丝流量时,本发明所针对的烘丝干头过程开始,当5处检测到出口烟丝含水率高于13%时,本发明所针对的烘丝干头过程结束。针对该烟草烘丝干头过程,本发明所述基于长期预测的烘丝干头过程建模和优化设定控制方法,具体包括以下步骤:
步骤1:用于对烟草烘丝干头过程进行长期预测的TPS-ARX模型结构为:
其中:yt为烘丝干头过程t时刻的出口烟丝含水率,为烘丝干头过程t时刻的排潮风门开度,/>为t时刻的热风风温,/>为t时刻的入口烟丝流量,/>为t时刻的入口烟丝含水率,ξt为t时刻的建模误差;xt-1为TPS-ARX模型的状态量,且/>φ0(xt-1)和/>均为关于xt-1的薄板样条(ThinPlate Spline,TPS)函数型系数,且具体结构如下:
步骤:2:建立基于步骤1中烘丝干头过程TPS-ARX模型的多步向前预测模型:
步骤3:优化步骤2中建立的烘丝干头过程长期预测TPS-ARX模型(4)的参数集θL和θN:
首先,定义模型参数优化的目标函数如下:
上式中,为第β组烘丝干头过程采样数据集对应的TPS-ARX模型(3)的n步向前预测的建模误差向量;/>是采集的第β组烘丝干头过程数据集,在第τ+p时刻的出口烟丝含水率实际值,且τ=6;n为该数据集的长度,本具体实施例中n=700;且β=1,2,…,α,本具体实施例中α=18;是采集的第β组烘丝干头过程数据集,在第τ时刻的基于步骤2中TPS-ARX模型(3)的p步向前预测的输出值。
其次,采用高斯牛顿法求解如下非线性优化问题:
上述高斯牛顿法优化过程的终止条件如下:
步骤4:基于TPS-ARX模型长期预测的烘丝干头过程优化设定控制方法:
上式中,{x1,x2,x3,x4},{x5,x6,x7,x8}和{z1,z2,z3,z4}分别为函数f1(t),f2(t)和f3(t)中的待优化参数集。定义x={x1,x2,x3,x4,x5,x6,x7,x8},z={z1,z2,z3,z4},则参数集x和z的具体优化方法如下:
Step 1,将上述控制变量的设定函数曲线代入已建立的干头过程TPS-ARX长期预测模型,可得在该输入曲线情况下的模型预测输出为:则基于采集的第β组烘丝干头过程数据序列,得到的相应TPS-ARX模型的预测输出为:
Step 2,定义优化参数集x和z的目标函数如下:
Step 3,采用高斯牛顿法求解如下非线性优化问题:
上述高斯牛顿法优化过程的终止条件为:
Claims (3)
1.一种基于长期预测的烘丝干头过程建模和优化设定控制方法,其特征在于,包括以下步骤:
1)建立用于对烟草烘丝干头过程进行预测的TPS-ARX模型结构如下:
上式中,yt为烘丝干头过程t时刻的出口烟丝含水率,为烘丝干头过程t时刻的排潮风门开度,/>为t时刻的热风风温,/>为t时刻的入口烟丝流量,/>为t时刻的入口烟丝含水率,ξt为t时刻的建模误差;xt-1为TPS-ARX模型的状态量,且/>φ0(xt-1)和均为关于xt-1的薄板样条(ThinPlate Spline,TPS)函数型系数,且具体结构为:
3)优化步骤2)中建立的烘丝干头过程长期预测TPS-ARX模型的参数集θL和θN;
4)基于步骤2)中建立的烘丝干头过程长期预测TPS-ARX模型,通过选择待优化的控制变量函数曲线,设计烘丝干头过程的最优设定控制方法。
2.根据权利要求1所述的基于长期预测的烘丝干头过程建模和优化设定控制方法,其特征在于,步骤3)中,用于优化步骤2)中建立的TPS-ARX模型的参数集θL和θN,具体方法如下:
(1)定义模型参数优化的目标函数如下:
(2)采用高斯牛顿法求解非线性优化问题:
3.根据权利要求1所述的基于长期预测的烘丝干头过程建模和优化设定控制方法,其特征在于,步骤4)中,烘丝干头过程最优设定控制方法的具体设计步骤如下:
上式中,{x1,x2,x3,x4},{x5,x6,x7,x8}和{z1,z2,z3,z4}分别为函数f1(t),f2(t)和f3(t)中的待优化参数集,并定义x={x1,x2,x3,x4,x5,x6,x7,x8},z={z1,z2,z3,z4},且参数集x和z的具体优化方法如下:
Step 1,将上述控制变量的设定函数曲线代入已建立的干头过程TPS-ARX长期预测模型,可得在该设定函数曲线代入情况下的模型预测输出为:则基于采集的第β组烘丝干头过程数据序列,得到的相应TPS-ARX模型的预测输出为:
Step 2,定义优化参数集x和z的目标函数如下:
Step 3,采用高斯牛顿法求解如下非线性优化问题:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911083453.4A CN110826229B (zh) | 2019-11-07 | 2019-11-07 | 基于长期预测的烘丝干头过程建模和优化设定控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911083453.4A CN110826229B (zh) | 2019-11-07 | 2019-11-07 | 基于长期预测的烘丝干头过程建模和优化设定控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110826229A CN110826229A (zh) | 2020-02-21 |
CN110826229B true CN110826229B (zh) | 2023-06-09 |
Family
ID=69553405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911083453.4A Active CN110826229B (zh) | 2019-11-07 | 2019-11-07 | 基于长期预测的烘丝干头过程建模和优化设定控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110826229B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113017132A (zh) * | 2021-04-09 | 2021-06-25 | 红云红河烟草(集团)有限责任公司 | 一种基于烘丝机工艺参数预测的烟丝质量优化的方法 |
CN115167582B (zh) * | 2022-07-11 | 2023-07-14 | 浪潮工业互联网股份有限公司 | 一种基于数字孪生的玫瑰烘干处理监控方法及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2003282218A1 (en) * | 2002-11-08 | 2004-06-07 | British American Tobacco (Investments) Limited | Method and system for predicting constituent yields in tobacco smoke using a multivariate regression model |
CN102871214A (zh) * | 2012-10-08 | 2013-01-16 | 秦皇岛烟草机械有限责任公司 | 基于模型预测的烘丝机出口水分控制方法 |
CN203369376U (zh) * | 2013-06-27 | 2014-01-01 | 红云红河烟草(集团)有限责任公司 | 烘叶丝干头干尾流动回潮装置 |
CN103610227A (zh) * | 2013-12-09 | 2014-03-05 | 中南大学 | 一种烘丝机头尾段工艺变量优化控制方法 |
CN109034388A (zh) * | 2018-07-27 | 2018-12-18 | 湖北中烟工业有限责任公司 | 一种基于遗传算法优化神经网络的卷烟材料与主流烟气成分的预测模型 |
-
2019
- 2019-11-07 CN CN201911083453.4A patent/CN110826229B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2003282218A1 (en) * | 2002-11-08 | 2004-06-07 | British American Tobacco (Investments) Limited | Method and system for predicting constituent yields in tobacco smoke using a multivariate regression model |
CN102871214A (zh) * | 2012-10-08 | 2013-01-16 | 秦皇岛烟草机械有限责任公司 | 基于模型预测的烘丝机出口水分控制方法 |
CN203369376U (zh) * | 2013-06-27 | 2014-01-01 | 红云红河烟草(集团)有限责任公司 | 烘叶丝干头干尾流动回潮装置 |
CN103610227A (zh) * | 2013-12-09 | 2014-03-05 | 中南大学 | 一种烘丝机头尾段工艺变量优化控制方法 |
CN109034388A (zh) * | 2018-07-27 | 2018-12-18 | 湖北中烟工业有限责任公司 | 一种基于遗传算法优化神经网络的卷烟材料与主流烟气成分的预测模型 |
Also Published As
Publication number | Publication date |
---|---|
CN110826229A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826229B (zh) | 基于长期预测的烘丝干头过程建模和优化设定控制方法 | |
CN103610227B (zh) | 一种烘丝机头尾段工艺变量优化控制方法 | |
CN100584224C (zh) | 一种烘丝过程控制方法 | |
CN110946314B (zh) | 一种气流烘丝的烟丝含水量的控制方法及系统 | |
CN103202528B (zh) | 基于ht水分补偿的前馈和反馈的叶丝干燥控制方法及系统 | |
CN112914139B (zh) | 一种松散回潮工序的加水量的控制方法及系统 | |
CN106773657B (zh) | 一种滚筒干燥非稳态加工过程头料和尾料阶段的调控优化方法 | |
CN111109643A (zh) | 一种薄板烘丝的烟丝含水量的控制方法及系统 | |
CN102488308A (zh) | 烘丝机水分先进协调控制系统 | |
CN110360838B (zh) | 用于气流干燥系统的稳定控制产品水分含量方法 | |
CN114115393A (zh) | 一种制丝线薄板烘丝机出口水分和温度的控制方法 | |
CN112790418A (zh) | 一种叶丝超级回潮的加水控制方法及系统 | |
CN110209217B (zh) | 一种基于不定周期的pid的蒸箱温度的控制方法 | |
CN101480271A (zh) | 一种对烟丝进行滚筒式干燥精确性加工的方法及其装置 | |
CN112890260B (zh) | 一种基于滑窗预测的薄板烘丝出口含水率的控制方法 | |
CN105404143B (zh) | 燃气炉窑炉膛压力计算机智能模糊控制节能方法 | |
CN110673490A (zh) | 烘丝干尾过程的长期预测建模及优化设定控制方法 | |
CN112790421B (zh) | 一种基于滑窗预测的梗丝加料出口含水率的控制方法 | |
CN203243922U (zh) | 基于ht水分补偿的前馈和反馈的叶丝干燥控制系统 | |
CN113040411B (zh) | 一种参数自优化的烟机蒸汽加热控制模型及其构建方法 | |
CN110580326B (zh) | 烘丝机烟草烘丝中间过程的建模方法 | |
CN112826119B (zh) | 一种水浸式烟梗回潮的水分平衡控制方法及系统 | |
CN112826120B (zh) | 一种气流式叶丝干燥的前馈控制方法及系统 | |
CN113925190B (zh) | 一种烘丝机水分调节的控制方法 | |
CN113017137B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |