CN102222128A - 一种废塑料炼油燃烧优化的方法 - Google Patents

一种废塑料炼油燃烧优化的方法 Download PDF

Info

Publication number
CN102222128A
CN102222128A CN2011101154107A CN201110115410A CN102222128A CN 102222128 A CN102222128 A CN 102222128A CN 2011101154107 A CN2011101154107 A CN 2011101154107A CN 201110115410 A CN201110115410 A CN 201110115410A CN 102222128 A CN102222128 A CN 102222128A
Authority
CN
China
Prior art keywords
reactor
model
function
burner
generalized regression
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.)
Granted
Application number
CN2011101154107A
Other languages
English (en)
Other versions
CN102222128B (zh
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.)
QUZHOU DISTANT VIEW RESOURCE REGENERATION TECHNOLOGY Co Ltd
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN 201110115410 priority Critical patent/CN102222128B/zh
Publication of CN102222128A publication Critical patent/CN102222128A/zh
Application granted granted Critical
Publication of CN102222128B publication Critical patent/CN102222128B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Processing Of Solid Wastes (AREA)
  • Incineration Of Waste (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种废塑料炼油燃烧优化的方法。现有方法主要靠工作人员经验。本发明方法通过废塑料炼油主体工艺的数据采集,针对反应釜轴向温度分布与各燃烧器运行参数、进料量和出产品量间的相互影响关系,建立基于数据挖掘技术的模型,应用并行的优化算法寻优等手段,确立了一种废塑料炼油燃烧优化的方法。利用本发明方法可有效控制废塑料炼油过程中反应釜内温度的分布情况,提高反应效率和产品质量,既可以实施离线优化也可以进行在线实时燃烧优化。

Description

一种废塑料炼油燃烧优化的方法
技术领域
本发明属于信息与控制技术领域,涉及到自动化技术,特别是涉及一种废塑料炼油燃烧优化方法。
背景技术
废塑料炼油的燃烧优化是控制废塑料裂解反应和产品的重要技术手段,其目标是在一定生产条件和目标下,通过调整各燃烧器的运行参数而获得反应釜所需的理想的燃烧状态,使反应釜的温度分布良好,废塑料的裂解反应在适当的温度下进行,使产品符合生产需求的基础上利益最大化。加热反应釜的各燃烧器的给风和给油等运行参数的不同对反应釜内的温度分布有直接的影响,不同的燃烧器的给油和给风的配合会直接导致不同的反应釜内的温度分布的情况,尤其是在反应釜进料有扰动的情况下,温度分布更不稳定。对于一定的生产条件和产品需求下,针对反应釜所需的理想的燃烧状态特征指标,存在一种最优的各燃烧器运行参数配置方案,能够使相应燃烧状态的特征指标最优化,但是,反应釜内的温度分布与各燃烧器运行参数和进料量及出产品量间有着非常复杂的耦合关系,要找到最优的各燃烧器的运行参数的配置并不容易。废塑料炼油是一种新的生产过程,其燃烧优化问题还没有得到解决。
实际生产中废塑料炼油的燃烧主要是靠工作人员经验,其目标也只是维持生产正常进行,并没有优化的概念,因此其生产过程中燃烧状态还有很大的可提升空间。
通过数据挖掘,在大量不同的生产运行参数组合中,应用机器学习的方法,挖掘出各燃烧器的运行参数、进料量和出产品量与反应釜内轴向温度分布间的关系模型,再结合优化算法对进行燃烧优化是非常有潜力的方法。如何使该方法真正达到废塑料炼油生产实际的要求,是该技术的关键,主要难题包括,如何提高模型的预测和泛化能力,如何提高模型的自学习能力。
发明内容
本发明的目标是针对废塑料炼油燃烧优化,提出一种燃烧优化方法,优化内容兼顾了各燃烧器的运行参数、进料情况和出产品情况。
本发明的技术方案是通过废塑料炼油主体工艺的数据采集,针对反应釜轴向温度分布与各燃烧器运行参数、进料量和出产品量间的相互影响关系,建立基于数据挖掘技术的模型,应用并行的优化算法寻优等手段,确立了一种废塑料炼油燃烧优化的方法,利用该方法可有效控制废塑料炼油过程中,反应釜内温度的分布情况,提高反应效率和产品质量,既可以实施离线优化也可以进行在线实时燃烧优化。
本发明的具体步骤是:
步骤(1)采集废塑料炼油生产过程中,各燃烧器运行参数、反应釜进料量和反应釜出产品量及反应釜所需要的燃烧状态的特征指标,建立实时数据库;具体的各燃烧器运行参数、反应釜进料量、反应釜出产品量和反应釜所需要的燃烧状态的特征指标,通过废塑料炼油生产过程中实时的数据监测控制系统获取,或直接通过仪器设备测量采集。
所述的各燃烧器运行参数为各燃烧器的给风量和给油量;所述的反应釜所需要的燃烧状态的特征指标为反应釜内轴向温度分布(反应釜内壁轴向取M个温度监测点,M≥3)。
这个过程积累数据,在积累了相当的数据后,再进行后边的数据挖掘建模。
步骤(2)采用支持向量机方法和广义回归神经网络建模,建立反应釜轴向的温度分布与各燃烧器运行参数、反应釜进料量和反应釜出产品量间的模型;具体建模方法如下:
用于建模样本的输入参数及表征反应釜所需要的燃烧状态的特征指标的输出参数表示为                                               
Figure 2011101154107100002DEST_PATH_IMAGE002
,其中表示第
Figure 2011101154107100002DEST_PATH_IMAGE006
组作为输入数据的各燃烧器的运行参数 (包括各燃烧器的给油和给风量) 、反应釜进料量和反应釜出产品量向量,
Figure 2011101154107100002DEST_PATH_IMAGE008
表示第
Figure 352841DEST_PATH_IMAGE006
组作为输出参数的表征反应釜所需要的燃烧状态的特征指标的参数向量(M个反应釜轴向的温度值),为样本数量,以实际运行数据为基础建立各燃烧器运行参数、反应釜进料量和反应釜出产品量与反应釜所需要的燃烧状态的特征指标间的模型。
采用支持向量机方法建模,支持向量机核函数选为径向基函数
Figure 2011101154107100002DEST_PATH_IMAGE012
Figure 2011101154107100002DEST_PATH_IMAGE014
为映射函数,设所求的目标函数为:
Figure 2011101154107100002DEST_PATH_IMAGE016
Figure 2011101154107100002DEST_PATH_IMAGE018
为模型输出的锅炉燃烧状态的特征指标预测值,为权重系数向量,
Figure 2011101154107100002DEST_PATH_IMAGE022
为截距。引入松弛因子ξ* i≥0和ξi≥0和允许拟合误差ε,模型可以通过在约束:
Figure 2011101154107100002DEST_PATH_IMAGE024
Figure 2011101154107100002DEST_PATH_IMAGE026
,条件下,最小化:
Figure 2011101154107100002DEST_PATH_IMAGE028
获得,其中常数C>0为惩罚系数。该最小化问题为一个凸二次规划问题,引入拉格朗日函数:
Figure 2011101154107100002DEST_PATH_IMAGE030
其中:
Figure 2011101154107100002DEST_PATH_IMAGE032
≥0, 
Figure 2011101154107100002DEST_PATH_IMAGE034
≥0,为拉格朗日乘数。
在鞍点处,函数L是关于w,b,ξii *的极小点,也是
Figure 2011101154107100002DEST_PATH_IMAGE036
极大点,最小化问题转化为求其对偶问题的最大化问题。
拉格朗日函数L在鞍点处是关于w,b,ξii *极小点,得:
可得拉格朗日函数的对偶函数:
Figure 2011101154107100002DEST_PATH_IMAGE040
此时,
Figure 2011101154107100002DEST_PATH_IMAGE042
Figure 2011101154107100002DEST_PATH_IMAGE044
按照库恩-塔克(KKT)条件定理,在鞍点有下式成立:
Figure 2011101154107100002DEST_PATH_IMAGE048
由上式可见,α i ·α i * =0,α i α i * 都不会同时为非零,可得:
Figure 2011101154107100002DEST_PATH_IMAGE050
  
Figure 403404DEST_PATH_IMAGE026
从上式可求出b,获得模型。
在采用广义回归神经网络建模,对于
Figure 2011101154107100002DEST_PATH_IMAGE052
个隐节点的3层广义回归神经网络采用径向基影射函数,其输出为:
Figure 2011101154107100002DEST_PATH_IMAGE054
为权重系数, 
Figure 2011101154107100002DEST_PATH_IMAGE058
Figure 2011101154107100002DEST_PATH_IMAGE060
维输入向量,
Figure 2011101154107100002DEST_PATH_IMAGE062
Figure 2011101154107100002DEST_PATH_IMAGE064
为径向基函数参数。建立广义回归神经网络模型的关键在于确定函数的参数
Figure 619359DEST_PATH_IMAGE062
Figure 14569DEST_PATH_IMAGE064
及权重系数
Figure 828941DEST_PATH_IMAGE056
。采用遗传算法迭代训练广义回归神经网络,定义遗传算法初始群体
Figure 2011101154107100002DEST_PATH_IMAGE066
向量的各维分量,分别为隐节点个、径向基函数的参数
Figure 589086DEST_PATH_IMAGE062
和权重系数,目标函数为:
Figure 2011101154107100002DEST_PATH_IMAGE068
,其中
Figure 2011101154107100002DEST_PATH_IMAGE070
为第个样本的BP神经网络输出值,为第
Figure 31777DEST_PATH_IMAGE006
个样本的实际值。当
Figure 2011101154107100002DEST_PATH_IMAGE072
达到了最小、达到设定值或完成迭代次数时,训练完成,获得隐节点个数、径向基函数的参数和权重系数,从而获得广义回归神经网络模型。
采用新的废塑料炼油不同运行状态下的数据作为检验样本,应用支持向量机模型和广义回归神经网络模型加权平均的预测方法,对检验数据进行预测,即
Figure 2011101154107100002DEST_PATH_IMAGE074
,其中
Figure 910740DEST_PATH_IMAGE008
为第组检验样本工况的目标预测值,
Figure 2011101154107100002DEST_PATH_IMAGE076
为支持向量机模型预测值,
Figure 2011101154107100002DEST_PATH_IMAGE078
为广义回归神经网络模型预测值,
Figure 2011101154107100002DEST_PATH_IMAGE080
为支持向量机模型预测值比例系数,
Figure 2011101154107100002DEST_PATH_IMAGE082
为广义回归神经网络模型的预测比例系数,且
Figure 2011101154107100002DEST_PATH_IMAGE084
Figure 423334DEST_PATH_IMAGE080
Figure 853178DEST_PATH_IMAGE082
的确定采用粒子群算法迭代寻优确定,初始化粒子群位置向量
Figure 2011101154107100002DEST_PATH_IMAGE086
的各维分量,分别为支持向量机模型比例
Figure 32487DEST_PATH_IMAGE080
和广义回归神经网络模型权重
Figure 307611DEST_PATH_IMAGE082
,目标函数为:,其中
Figure 2011101154107100002DEST_PATH_IMAGE090
为第
Figure 75715DEST_PATH_IMAGE006
组工况实际数据与结合模型预测的反应釜轴向各测点中的最大误差,
Figure 2011101154107100002DEST_PATH_IMAGE092
为第
Figure 597832DEST_PATH_IMAGE006
组工况实际数据与结合模型预测的反应釜轴向各测点的平均误差,当
Figure 2011101154107100002DEST_PATH_IMAGE094
取得了最小、达到设定值或完成迭代次数时,寻优完成,获得最优的支持向量机模型和广义回归神经网络模型的比例系数,按最优系数将支持向量机模型和广义回归神经网络模型结合,获得废塑料炼油燃烧特性模型。
步骤(3) 利用粒子群优化算法结合所建燃烧特性模型,针对不同反应釜所需燃烧指标进行各燃烧参数配置的优化,具体步骤如下:
a.定义粒子群位置向量的各维分量分别为各燃烧器运行参数;
b.设定粒子群的搜索目标和迭代次数,搜索目标为根据实际生产状况确定的反应釜所需的最优燃烧状态指标;
c.根据实际燃烧器的设计和运行要求设定各各燃烧器运行参数的寻优范围,并初始化位置向量,然后根据上一步设定的搜索目标,用粒子群算法进行迭代计算,搜索粒子群在参数向量空间内的最优位置;
d.当粒子群算法完成迭代次数或找到设定要求的最优时,停止计算获得相应最优的位置向量,即获得最优的各燃烧器运行参数组合。
废塑料炼油是一种新的工业生产工艺,由于其主要是进行化学裂解反应的过程,所以对温度的要求比较严格,目前实际的生产过程中对燃烧的调整完全依据个人经验,在进料量的扰动下,具有不稳定性,有待提高。本发明方法具体是采集废塑料炼油的实时数据,利用数据挖掘技术和并行优化算法,针对反应釜轴向的温度分布与各燃烧器运行参数、反应釜进料量和反应釜出产品量间的关系进行建模,结合寻优算法对各燃烧器的运行进行优化,以达到最优的控制反应釜温度分布的目标。
本发明方法既可以在线优化也可以离线优化。
具体实施方式
一种废塑料炼油燃烧优化方法,具体是以下步骤:
(1) 采集废塑料炼油生产过程中,各燃烧器运行参数、反应釜进料量和反应釜出产品量及反应釜所需要的燃烧状态的特征指标,建立实时数据库;具体的各燃烧器运行参数、反应釜进料量、反应釜出产品量和反应釜所需要的燃烧状态的特征指标,通过废塑料炼油生产过程中实时的数据监测控制系统获取,或直接通过仪器设备测量采集。
所述的各燃烧器运行参数为各燃烧器的给风量和给油量;所述的反应釜所需要的燃烧状态的特征指标为反应釜内轴向温度分布(反应釜内壁轴向取M个温度监测点,M≥3),其获得的方法为成熟技术。
(2)采用支持向量机方法和广义回归神经网络建模,建立反应釜轴向的温度分布与各燃烧器运行参数、反应釜进料量和反应釜出产品量间的模型;具体建模方法如下:
用于建模样本的输入参数及表征反应釜所需要的燃烧状态的特征指标的输出参数表示为
Figure 28180DEST_PATH_IMAGE002
,其中
Figure 799827DEST_PATH_IMAGE004
表示第
Figure 750465DEST_PATH_IMAGE006
组作为输入数据的各燃烧器的运行参数 (包括各燃烧器的给油和给风量) 、反应釜进料量和反应釜出产品量向量,
Figure 570653DEST_PATH_IMAGE008
表示第组作为输出参数的表征反应釜所需要的燃烧状态的特征指标的参数向量(M个反应釜轴向的温度值),
Figure 55566DEST_PATH_IMAGE010
为样本数量,以实际运行数据为基础建立各燃烧器运行参数、反应釜进料量和反应釜出产品量与反应釜所需要的燃烧状态的特征指标间的模型。
采用支持向量机方法建模,支持向量机核函数选为径向基函数
Figure 914118DEST_PATH_IMAGE014
为映射函数,设所求的目标函数为:
Figure 324370DEST_PATH_IMAGE016
为模型输出的锅炉燃烧状态的特征指标预测值,
Figure 550001DEST_PATH_IMAGE020
为权重系数向量,
Figure 774309DEST_PATH_IMAGE022
为截距。引入松弛因子ξ* i≥0和ξi≥0和允许拟合误差ε,模型可以通过在约束:
Figure 835806DEST_PATH_IMAGE024
Figure 57840DEST_PATH_IMAGE026
,条件下,最小化:
Figure 470367DEST_PATH_IMAGE028
获得,其中常数C>0为惩罚系数。该最小化问题为一个凸二次规划问题,引入拉格朗日函数:
其中:
Figure 866899DEST_PATH_IMAGE032
≥0, ≥0,为拉格朗日乘数。
在鞍点处,函数L是关于w,b,ξii *的极小点,也是
Figure 894078DEST_PATH_IMAGE036
极大点,最小化问题转化为求其对偶问题的最大化问题。
拉格朗日函数L在鞍点处是关于w,b,ξii *极小点,得:
Figure 460188DEST_PATH_IMAGE038
可得拉格朗日函数的对偶函数:
Figure 761857DEST_PATH_IMAGE040
此时,
Figure 309382DEST_PATH_IMAGE042
Figure 696501DEST_PATH_IMAGE044
按照库恩-塔克(KKT)条件定理,在鞍点有下式成立:
Figure 894581DEST_PATH_IMAGE048
由上式可见,α i ·α i * =0,α i α i * 都不会同时为非零,可得:
Figure 691635DEST_PATH_IMAGE050
  
Figure 478633DEST_PATH_IMAGE026
从上式可求出b,获得模型。
在采用广义回归神经网络建模,对于
Figure 652125DEST_PATH_IMAGE052
个隐节点的3层广义回归神经网络采用径向基影射函数,其输出为:
Figure 662807DEST_PATH_IMAGE054
为权重系数, 
Figure 664578DEST_PATH_IMAGE058
Figure 563133DEST_PATH_IMAGE060
维输入向量,
Figure 428320DEST_PATH_IMAGE062
Figure 567178DEST_PATH_IMAGE064
为径向基函数参数。建立广义回归神经网络模型的关键在于确定函数的参数
Figure 88289DEST_PATH_IMAGE062
Figure 338005DEST_PATH_IMAGE064
及权重系数
Figure 323278DEST_PATH_IMAGE056
。采用遗传算法迭代训练广义回归神经网络,定义遗传算法初始群体向量的各维分量,分别为隐节点个、径向基函数的参数
Figure 890712DEST_PATH_IMAGE062
Figure 881801DEST_PATH_IMAGE064
和权重系数,目标函数为:
Figure 721581DEST_PATH_IMAGE068
,其中
Figure 202241DEST_PATH_IMAGE070
为第
Figure 681633DEST_PATH_IMAGE006
个样本的BP神经网络输出值,为第
Figure 170700DEST_PATH_IMAGE006
个样本的实际值。当达到了最小、达到设定值或完成迭代次数时,训练完成,获得隐节点个数、径向基函数的参数和权重系数,从而获得广义回归神经网络模型。
采用新的废塑料炼油不同运行状态下的数据作为检验样本,应用支持向量机模型和广义回归神经网络模型加权平均的预测方法,对检验数据进行预测,即
Figure 867578DEST_PATH_IMAGE074
,其中
Figure 452667DEST_PATH_IMAGE008
为第
Figure 1460DEST_PATH_IMAGE006
组检验样本工况的目标预测值,
Figure 89502DEST_PATH_IMAGE076
为支持向量机模型预测值,
Figure 294219DEST_PATH_IMAGE078
为广义回归神经网络模型预测值,
Figure 493119DEST_PATH_IMAGE080
为支持向量机模型预测值比例系数,
Figure 83369DEST_PATH_IMAGE082
为广义回归神经网络模型的预测比例系数,且
Figure 342312DEST_PATH_IMAGE084
Figure 96641DEST_PATH_IMAGE080
Figure 771336DEST_PATH_IMAGE082
的确定采用粒子群算法迭代寻优确定,初始化粒子群位置向量
Figure 294722DEST_PATH_IMAGE086
的各维分量,分别为支持向量机模型比例
Figure 911517DEST_PATH_IMAGE080
和广义回归神经网络模型权重,目标函数为:
Figure 428266DEST_PATH_IMAGE088
,其中
Figure 9420DEST_PATH_IMAGE090
为第
Figure 344586DEST_PATH_IMAGE006
组工况实际数据与结合模型预测的反应釜轴向各测点中的最大误差,
Figure 73508DEST_PATH_IMAGE092
为第
Figure 339273DEST_PATH_IMAGE006
组工况实际数据与结合模型预测的反应釜轴向各测点的平均误差,当取得了最小、达到设定值或完成迭代次数时,寻优完成,获得最优的支持向量机模型和广义回归神经网络模型的比例系数,按最优系数将支持向量机模型和广义回归神经网络模型结合,获得废塑料炼油燃烧特性模型。
步骤(3) 利用粒子群优化算法结合所建燃烧特性模型,针对不同反应釜所需燃烧指标进行各燃烧参数配置的优化,具体步骤如下:
a.定义粒子群位置向量
Figure 281001DEST_PATH_IMAGE086
的各维分量分别为各燃烧器运行参数;
b.设定粒子群的搜索目标和迭代次数,搜索目标为根据实际生产状况确定的反应釜所需的最优燃烧状态指标;
c.根据实际燃烧器的设计和运行要求设定各各燃烧器运行参数的寻优范围,并初始化位置向量
Figure 497219DEST_PATH_IMAGE086
,然后根据上一步设定的搜索目标,用粒子群算法进行迭代计算,搜索粒子群在参数向量空间内的最优位置;
d.当粒子群算法完成迭代次数或找到设定要求的最优时,停止计算获得相应最优的位置向量,即获得最优的各燃烧器运行参数组合。根据所获得的最优参数组合对各燃烧器实际燃烧进行调整,达到燃烧优化的目的。

Claims (1)

1.一种废塑料炼油燃烧优化的方法,其特征在于该方法的具体步骤包括:
步骤(1)采集废塑料炼油生产过程中,各燃烧器运行参数、反应釜进料量和反应釜出产品量及反应釜所需要的燃烧状态的特征指标,建立实时数据库;具体的各燃烧器运行参数、反应釜进料量、反应釜出产品量和反应釜所需要的燃烧状态的特征指标,通过废塑料炼油生产过程中实时的数据监测控制系统获取,或直接通过仪器设备测量采集;
所述的各燃烧器运行参数为各燃烧器的给风量和给油量;所述的反应釜所需要的燃烧状态的特征指标为反应釜内轴向M个温度监测点的温度分布,M≥3; 
步骤(2)采用支持向量机方法和广义回归神经网络建模,建立反应釜轴向的温度分布与各燃烧器运行参数、反应釜进料量和反应釜出产品量间的模型;具体建模方法如下:
用于建模样本的输入参数及表征反应釜所需要的燃烧状态的特征指标的输出参数表示为                                               
Figure 2011101154107100001DEST_PATH_IMAGE002
,其中
Figure 2011101154107100001DEST_PATH_IMAGE004
表示第
Figure 2011101154107100001DEST_PATH_IMAGE006
组作为输入数据的各燃烧器的运行参数、反应釜进料量和反应釜出产品量向量,
Figure 2011101154107100001DEST_PATH_IMAGE008
表示第组作为输出参数的表征反应釜所需要的燃烧状态的特征指标的参数向量,
Figure 2011101154107100001DEST_PATH_IMAGE010
为样本数量,以实际运行数据为基础建立各燃烧器运行参数、反应釜进料量和反应釜出产品量与反应釜所需要的燃烧状态的特征指标间的模型;
采用支持向量机方法建模,支持向量机核函数选为径向基函数
Figure 2011101154107100001DEST_PATH_IMAGE012
Figure 2011101154107100001DEST_PATH_IMAGE014
为映射函数,设所求的目标函数为:
Figure 2011101154107100001DEST_PATH_IMAGE016
Figure 2011101154107100001DEST_PATH_IMAGE018
为模型输出的锅炉燃烧状态的特征指标预测值,
Figure 2011101154107100001DEST_PATH_IMAGE020
为权重系数向量,为截距;引入松弛因子ξ* i≥0和ξi≥0和允许拟合误差ε,模型可以通过在约束:
Figure 2011101154107100001DEST_PATH_IMAGE024
Figure 2011101154107100001DEST_PATH_IMAGE026
,条件下,最小化:
Figure 2011101154107100001DEST_PATH_IMAGE028
获得,其中常数C>0为惩罚系数;引入拉格朗日函数:
Figure 2011101154107100001DEST_PATH_IMAGE030
     
其中:
Figure 2011101154107100001DEST_PATH_IMAGE032
≥0, 
Figure 2011101154107100001DEST_PATH_IMAGE034
≥0,为拉格朗日乘数;
在鞍点处,函数L是关于w,b,ξii *的极小点,也是
Figure 2011101154107100001DEST_PATH_IMAGE036
极大点,最小化问题转化为求其对偶问题的最大化问题;
拉格朗日函数L在鞍点处是关于w,b,ξii *极小点,得:
                    
可得拉格朗日函数的对偶函数:
           
此时,
Figure 2011101154107100001DEST_PATH_IMAGE042
Figure 2011101154107100001DEST_PATH_IMAGE044
按照库恩-塔克(KKT)条件定理,在鞍点有下式成立:
Figure 2011101154107100001DEST_PATH_IMAGE046
Figure 2011101154107100001DEST_PATH_IMAGE048
             
可得:
Figure 2011101154107100001DEST_PATH_IMAGE050
  
从上式可求出b,获得模型;
在采用广义回归神经网络建模,对于
Figure DEST_PATH_IMAGE052
个隐节点的3层广义回归神经网络采用径向基影射函数,其输出为:
Figure DEST_PATH_IMAGE054
Figure DEST_PATH_IMAGE056
为权重系数, 
Figure DEST_PATH_IMAGE058
维输入向量,
Figure DEST_PATH_IMAGE062
Figure DEST_PATH_IMAGE064
为径向基函数参数;建立广义回归神经网络模型的关键在于确定函数的参数
Figure 722935DEST_PATH_IMAGE064
及权重系数;采用遗传算法迭代训练广义回归神经网络,定义遗传算法初始群体
Figure DEST_PATH_IMAGE066
向量的各维分量,分别为隐节点个、径向基函数的参数
Figure 401676DEST_PATH_IMAGE064
和权重系数,目标函数为:
Figure DEST_PATH_IMAGE068
,其中
Figure DEST_PATH_IMAGE070
为第
Figure 451540DEST_PATH_IMAGE006
个样本的BP神经网络输出值,
Figure 675848DEST_PATH_IMAGE008
为第
Figure 189875DEST_PATH_IMAGE006
个样本的实际值;当
Figure DEST_PATH_IMAGE072
达到了最小、达到设定值或完成迭代次数时,训练完成,获得隐节点个数、径向基函数的参数和权重系数,从而获得广义回归神经网络模型;
采用新的废塑料炼油不同运行状态下的数据作为检验样本,应用支持向量机模型和广义回归神经网络模型加权平均的预测方法,对检验数据进行预测,即
Figure DEST_PATH_IMAGE074
,其中
Figure 411909DEST_PATH_IMAGE008
为第
Figure 11386DEST_PATH_IMAGE006
组检验样本工况的目标预测值、
Figure DEST_PATH_IMAGE076
为支持向量机模型预测值、
Figure DEST_PATH_IMAGE078
为广义回归神经网络模型预测值、
Figure DEST_PATH_IMAGE080
为支持向量机模型预测值比例系数、
Figure DEST_PATH_IMAGE082
为广义回归神经网络模型的预测比例系数,且
Figure DEST_PATH_IMAGE084
Figure 839052DEST_PATH_IMAGE080
Figure 207586DEST_PATH_IMAGE082
的确定采用粒子群算法迭代寻优确定,初始化粒子群位置向量的各维分量,分别为支持向量机模型比例
Figure 600521DEST_PATH_IMAGE080
和广义回归神经网络模型权重
Figure 500344DEST_PATH_IMAGE082
,目标函数为:
Figure DEST_PATH_IMAGE088
,其中
Figure DEST_PATH_IMAGE090
为第组工况实际数据与结合模型预测的反应釜轴向各测点中的最大误差,
Figure DEST_PATH_IMAGE092
为第组工况实际数据与结合模型预测的反应釜轴向各测点的平均误差,当
Figure DEST_PATH_IMAGE094
取得了最小、达到设定值或完成迭代次数时,寻优完成,获得最优的支持向量机模型和广义回归神经网络模型的比例系数,按最优系数将支持向量机模型和广义回归神经网络模型结合,获得废塑料炼油燃烧特性模型;
步骤(3)利用粒子群优化算法结合所建燃烧特性模型,针对不同反应釜所需燃烧指标进行各燃烧参数配置的优化,具体步骤如下:
a.定义粒子群位置向量
Figure 307127DEST_PATH_IMAGE086
的各维分量分别为各燃烧器运行参数;
b.设定粒子群的搜索目标和迭代次数,搜索目标为根据实际生产状况确定的反应釜所需的最优燃烧状态指标;
c.根据实际燃烧器的设计和运行要求设定各各燃烧器运行参数的寻优范围,并初始化位置向量
Figure 694246DEST_PATH_IMAGE086
,然后根据上一步设定的搜索目标,用粒子群算法进行迭代计算,搜索粒子群在参数向量空间内的最优位置;
d.当粒子群算法完成迭代次数或找到设定要求的最优时,停止计算获得相应最优的位置向量,即获得最优的各燃烧器运行参数组合。
CN 201110115410 2011-05-05 2011-05-05 一种废塑料炼油燃烧优化的方法 Active CN102222128B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110115410 CN102222128B (zh) 2011-05-05 2011-05-05 一种废塑料炼油燃烧优化的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110115410 CN102222128B (zh) 2011-05-05 2011-05-05 一种废塑料炼油燃烧优化的方法

Publications (2)

Publication Number Publication Date
CN102222128A true CN102222128A (zh) 2011-10-19
CN102222128B CN102222128B (zh) 2013-01-02

Family

ID=44778678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110115410 Active CN102222128B (zh) 2011-05-05 2011-05-05 一种废塑料炼油燃烧优化的方法

Country Status (1)

Country Link
CN (1) CN102222128B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102981409A (zh) * 2012-12-11 2013-03-20 富通集团有限公司 一种同轴电缆护套机优化运行方法
CN103064370A (zh) * 2012-12-17 2013-04-24 富通集团有限公司 一种多同轴电缆护套机协同生产方法
CN104318303A (zh) * 2014-09-23 2015-01-28 杭州电子科技大学 遗传算法优化的rbf神经网络的焦化炉温度预测方法
CN105808945A (zh) * 2016-03-07 2016-07-27 杭州电子科技大学 一种混合智能锅炉效率燃烧优化方法
CN109033536A (zh) * 2018-06-29 2018-12-18 杭州电子科技大学 一种燃煤分质利用清洁处理系统中热解釜的建模方法
CN109145361A (zh) * 2018-06-29 2019-01-04 杭州电子科技大学 一种燃煤分质利用与清洁预处理系统的优化运行方法
CN111157139A (zh) * 2020-02-06 2020-05-15 北京航空航天大学 一种用于单连通燃烧场温度分布的可视化测量方法
CN117170334A (zh) * 2023-11-02 2023-12-05 钥准医药科技(启东)有限公司 用于药剂快速融合的智能控制方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101498459A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的建模方法
CN101498458A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉在线燃烧优化的模型更新方法
CN101498457A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101498459A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的建模方法
CN101498458A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉在线燃烧优化的模型更新方法
CN101498457A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王春林: "《大型电站锅炉配煤及燃烧优化的支持向量机建模与实验研究》", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102981409A (zh) * 2012-12-11 2013-03-20 富通集团有限公司 一种同轴电缆护套机优化运行方法
CN103064370A (zh) * 2012-12-17 2013-04-24 富通集团有限公司 一种多同轴电缆护套机协同生产方法
CN103064370B (zh) * 2012-12-17 2015-05-06 富通集团有限公司 一种多同轴电缆护套机协同生产方法
CN104318303A (zh) * 2014-09-23 2015-01-28 杭州电子科技大学 遗传算法优化的rbf神经网络的焦化炉温度预测方法
CN105808945A (zh) * 2016-03-07 2016-07-27 杭州电子科技大学 一种混合智能锅炉效率燃烧优化方法
CN105808945B (zh) * 2016-03-07 2018-11-13 杭州电子科技大学 一种混合智能锅炉效率燃烧优化方法
CN109033536A (zh) * 2018-06-29 2018-12-18 杭州电子科技大学 一种燃煤分质利用清洁处理系统中热解釜的建模方法
CN109145361A (zh) * 2018-06-29 2019-01-04 杭州电子科技大学 一种燃煤分质利用与清洁预处理系统的优化运行方法
CN111157139A (zh) * 2020-02-06 2020-05-15 北京航空航天大学 一种用于单连通燃烧场温度分布的可视化测量方法
CN111157139B (zh) * 2020-02-06 2021-09-28 北京航空航天大学 一种用于单连通燃烧场温度分布的可视化测量方法
CN117170334A (zh) * 2023-11-02 2023-12-05 钥准医药科技(启东)有限公司 用于药剂快速融合的智能控制方法及系统
CN117170334B (zh) * 2023-11-02 2024-03-08 钥准医药科技(启东)有限公司 用于药剂快速融合的智能控制方法及系统

Also Published As

Publication number Publication date
CN102222128B (zh) 2013-01-02

Similar Documents

Publication Publication Date Title
CN102222128B (zh) 一种废塑料炼油燃烧优化的方法
CN102184287B (zh) 一种废塑料炼油燃烧优化的建模方法
CN101498457B (zh) 一种锅炉燃烧优化的方法
CN101498459B (zh) 一种锅炉燃烧优化的建模方法
CN102799778B (zh) 一种优化锅炉负荷分配方法
CN101286044B (zh) 一种燃煤锅炉系统蒸汽温度混合建模方法
CN106991539B (zh) 一种能源系统优化调度方法及装置
CN102842066A (zh) 一种生物质炉燃烧优化的建模方法
CN107016176A (zh) 一种混合智能锅炉综合燃烧优化方法
CN109886471A (zh) 基于神经网络与智能优化算法的火电机组负荷分配方法
CN102750424B (zh) 一种生物质炉燃烧优化的方法
CN102194057B (zh) 一种废塑料炼油燃烧优化的模型更新方法
CN117270483B (zh) 化工生产装置全流程动态优化控制方法、装置、电子设备
CN103631140A (zh) 基于性能评估的焦炉加热燃烧过程火道温度自动调节方法
CN105868867A (zh) 一种供热锅炉集群优化运行方法及系统
CN104102211A (zh) 一种燃煤热电厂热电联产系统的调度方法、服务器和系统
CN109507961A (zh) 一种半导体生产线动态负荷均衡投料控制方法
CN106155024A (zh) 一种无烟装煤集气管压力智能控制系统
Sun et al. Integrated control strategy of district heating system based on load forecasting and indoor temperature measurement
CN105808945B (zh) 一种混合智能锅炉效率燃烧优化方法
CN102184450A (zh) 一种多孔介质燃烧器燃烧优化的建模方法
CN208888602U (zh) 一种钢铁企业煤气管网运行优化装置
Wang et al. Soft-sensing modeling and intelligent optimal control strategy for distillation yield rate of atmospheric distillation oil refining process
CN106873377B (zh) 不确定扰动下的焦炉加热燃烧过程火道温度自动调节方法
CN109101683B (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
ASS Succession or assignment of patent right

Owner name: QUZHOU DISTANT VIEW RESOURCE REGENERATION TECHNOLO

Free format text: FORMER OWNER: HANGZHOU ELECTRONIC SCIENCE AND TECHNOLOGY UNIV

Effective date: 20150216

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 310018 HANGZHOU, ZHEJIANG PROVINCE TO: 324019 QUZHOU, ZHEJIANG PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20150216

Address after: 324019 Zhejiang province Quzhou city lotus town Qujiang District Geng Shanting

Patentee after: Quzhou Distant View Resource Regeneration Technology Co., Ltd.

Address before: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park No. 2 street

Patentee before: Hangzhou Electronic Science and Technology Univ