CN103838143B - 多模态全局最优丙烯聚合生产过程最优软测量系统和方法 - Google Patents
多模态全局最优丙烯聚合生产过程最优软测量系统和方法 Download PDFInfo
- Publication number
- CN103838143B CN103838143B CN201310659204.1A CN201310659204A CN103838143B CN 103838143 B CN103838143 B CN 103838143B CN 201310659204 A CN201310659204 A CN 201310659204A CN 103838143 B CN103838143 B CN 103838143B
- Authority
- CN
- China
- Prior art keywords
- network
- sub
- ant
- rbf
- solution
- 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.)
- Expired - Fee Related
Links
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 title claims abstract description 42
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 title claims abstract description 42
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 37
- 238000006116 polymerization reaction Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims description 25
- 238000005259 measurement Methods 0.000 claims abstract description 55
- 238000000513 principal component analysis Methods 0.000 claims abstract description 47
- 238000005457 optimization Methods 0.000 claims abstract description 38
- 238000013528 artificial neural network Methods 0.000 claims abstract description 27
- 238000000691 measurement method Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000012986 modification Methods 0.000 claims abstract description 5
- 230000004048 modification Effects 0.000 claims abstract description 5
- 238000012549 training Methods 0.000 claims description 37
- 239000011159 matrix material Substances 0.000 claims description 29
- 230000001537 neural effect Effects 0.000 claims description 29
- 238000003062 neural network model Methods 0.000 claims description 19
- 238000002844 melting Methods 0.000 claims description 15
- 230000008018 melting Effects 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 13
- 230000002087 whitening effect Effects 0.000 claims description 13
- 239000004743 Polypropylene Substances 0.000 claims description 12
- 229920001155 polypropylene Polymers 0.000 claims description 12
- -1 polypropylene Polymers 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 7
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 7
- 241000257303 Hymenoptera Species 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 239000012141 concentrate Substances 0.000 claims description 7
- 230000002068 genetic effect Effects 0.000 claims description 7
- 230000007774 longterm Effects 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 7
- 230000001105 regulatory effect Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 230000017105 transposition Effects 0.000 claims description 7
- 239000003054 catalyst Substances 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims 1
- 239000010410 layer Substances 0.000 description 7
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000007334 copolymerization reaction Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011127 biaxially oriented polypropylene Substances 0.000 description 1
- 229920006378 biaxially oriented polypropylene Polymers 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- HGAZMNJKRQFZKS-UHFFFAOYSA-N chloroethene;ethenyl acetate Chemical compound ClC=C.CC(=O)OC=C HGAZMNJKRQFZKS-UHFFFAOYSA-N 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000000088 plastic resin Substances 0.000 description 1
- 239000004800 polyvinyl chloride Substances 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 229920005992 thermoplastic resin Polymers 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种多模态全局最优丙烯聚合生产过程最优软测量系统,包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、多模态全局最优软测量系统以及熔融指数软测量值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连,与DCS数据库相连;最优软测量系统与DCS数据库及软测量值显示仪相连。所述的多模态全局最优软测量系统包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络多模优化模块。以及提供了一种用软测量系统实现的软测量方法。本发明实现在线测量、在线参数优化、软测量速度快、模型自动更新、抗干扰能力强、精度高。
Description
技术领域
本发明涉及一种最优软测量系统及方法,具体是一种多模态全局最优丙烯聚合生产过程最优软测量系统及方法。
背景技术
聚丙烯是由丙烯聚合而制得的一种热塑性树脂,丙烯最重要的下游产品,世界丙烯的50%,我国丙烯的65%都是用来制聚丙烯,是五大通用塑料之一,与我们的日常生活密切相关。聚丙烯是世界上增长最快的通用热塑性树脂,总量仅仅次于聚乙烯和聚氯乙烯。为使我国聚丙烯产品具有市场竞争力,开发刚性、韧性、流动性平衡好的抗冲共聚产品、无规共聚产品、BOPP和CPP薄膜料、纤维、无纺布料,及开发聚丙烯在汽车和家电领域的应用,都是今后重要的研究课题。
熔融指数是聚丙烯产品确定产品牌号的重要质量指标之一,它决定了产品的不同用途,对熔融指数的测量是聚丙烯生产中产品质量控制的一个重要环节,对生产和科研,都有非常重要的作用和指导意义。
然而,熔融指数的在线分析测量目前很难做到,一方面是在线熔融指数分析仪的缺乏,另一方面是现有的在线分析仪由于经常会堵塞而测量不准甚至无法正常使用所导致的使用上的困难。因此,目前工业生产中MI的测量,主要是通过人工取样、离线化验分析获得,而且一般每2-4小时只能分析一次,时间滞后大,给丙烯聚合生产的质量控制带来了困难,成为生产中急需解决的一个瓶颈问题。聚丙烯熔融指数的在线软测量系统及方法研究,从而成为学术界和工业界的一个前沿和热点。
发明内容
为了克服目前已有的丙烯聚合生产过程的测量精度不高、易受人为因素的影响的不足,本发明的目的在于提供一种在线测量、在线参数优化、软测量速度快、模型自动更新、抗干扰能力强、精度高的多模态全局最优丙烯聚合生产过程熔融指数最优软测量系统及方法。
本发明解决其技术问题所采用的技术方案是:
一种多模态全局最优丙烯聚合生产过程最优软测量系统,包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库、多模态全局最优软测量系统以及熔融指数软测量显示仪,所述现场智能仪表、控制站与丙烯聚合生产过程连接,所述现场智能仪表、控制站与DCS数据库连接,所述DCS数据库与多模态全局最优软测量系统的输入端连接,所述多模态全局最优软测量系统的输出端与熔融指数软测量显示仪连接,其特征在于:所述多模态全局最优软测量系统包括:
(1)、数据预处理模块,用于将从DCS数据库输入的模型输入变量进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(2)、PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵。若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(3)、神经网络模型模块,用于采用RBF神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;需要建立若干子神经网络,第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子网络预报结果,d(·)为实际结果。
从第二个子网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子网络的训练目标,Fi(·)为第i个网络的预报结果;d(·)为实际结果;F(·)为前i-1个子网络的综合结果;λ为调节参数,N为样本数目。
训练的终止条件为将得到的新的子网络加入多模神经网络后,网络群的预报误差不再减小。
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(i=1,2,…,n)
n为初始解的个数,sn为第n个初始解,k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据P来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回。
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数。
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能。
(4)、神经网络多模优化模块,用于对步骤(5.4)中的每个子网络进行构建
O(·)为模型输出,i为总的子网络数目,x为输入变量,Fj(·)为第j个子网络的输出;即最终多模RBF神经网络的预报结果为各个子网络预报结果的平均值。
作为优选的一种方案,所述多模态全局最优软测量模型还包括:模型更新模块,用于模型的在线更新,将定期将离线化验数据输入到训练集中,更新神经网络模型。
作为优选的再一种方案:在所述的智能连续空间蚁群算法训练多模RBF神经网络模型中,训练子RBF神经网络,然后将其构建起来形成神经网络群;由于子网络的选取标准是预报误差小、与其他的子网络差异大,所以这些预报效果好、又各不相同的子神经网络的综合预报效果能够具有更好的预报精度和稳定性。
作为优选的再一种方案:在PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
一种多模态全局最优聚丙烯生产过程最优软测量系统实现的软测量方法,所述软测量方法具体实现步骤如下:
(1)对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
(2)对样本数据进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(3)PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵。若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(4)基于模型输入、输出数据建立若干个初始子神经网络模型,采用RBF神经网络,通过误差最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子网络预报结果,d(·)为实际结果。
从第二个子网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子网络的训练目标,Fi(·)为第i个网络的预报结果;d(·)为实际结果;F(·)为前i-1个子网络的综合结果;λ为调节参数,N为样本数目。
训练的终止条件为将得到的新的子网络加入多模神经网络后,网络群的预报误差不再减小。
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(i=1,2,…,n)
n为初始解的个数,sn为第n个初始解,k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据P来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回。
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数。
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能。
(5)、神经网络多模优化模块,用于对步骤(4)中的每个子网络进行构建
O(·)为模型输出,i为总的子网络数目,x为输入变量,Fj(·)为第j个子网络的输出;即最终多模RBF神经网络的预报结果为各个子网络预报结果的平均值。
本发明的技术构思为:对丙烯聚合生产过程的重要质量指标熔融指数进行在线最优软测量,克服已有的聚丙烯熔融指数测量仪表测量精度不高、易受人为因素的影响的不足,通过智能连续空间蚁群算法训练多模RBF神经网络的方法来建立预报精度高、稳定性好的预报模型来得到最优的软测量结果。
本发明的有益效果主要表现在:1、在线测量;2、在线参数自动优化;3、软测量速度快;4、模型自动更新;5、抗干扰能力强;6、精度高。
附图说明
图1是多模态全局最优丙烯聚合生产过程最优软测量系统及方法的基本结构示意图;
图2是多模态全局最优软测量系统结构示意图;
图3是丙烯聚合生产过程Hypol工艺生产流程图。
具体实施方式
下面结合附图对本发明做进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
1.参照图1、图2和图3,一种多模态全局最优丙烯聚合生产过程最优软测量系统,包括丙烯聚合生产过程1、用于测量易测变量的现场智能仪表2、用于测量操作变量的控制站3、存放数据的DCS数据库4、多模态全局最优软测量系统5以及熔融指数软测量值显示仪6,所述现场智能仪表2、控制站3与丙烯聚合生产过程1连接,所述现场智能仪表2、控制站3与DCS数据库4连接,所述DCS数据库4与多模态全局最优软测量系统5的输入端连接,所述多模态全局最优软测量系统5的输出端与熔融指数软测量值显示仪6连接,所述多模态全局最优软测量系统包括:
(1)、数据预处理模块,用于将从DCS数据库输入的模型输入变量进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(2)、PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵。若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(3)、神经网络模型模块,用于采用RBF神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;需要建立若干子神经网络,第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子网络预报结果,d(·)为实际结果。
从第二个子网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子网络的训练目标,Fi(·)为第i个网络的预报结果;d(·)为实际结果;F(·)为前i-1个子网络的综合结果;λ为调节参数,N为样本数目。
训练的终止条件为将得到的新的子网络加入多模神经网络后,网络群的预报误差不再减小。
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(i=1,2,…,n)
n为初始解的个数,sn为第n个初始解,k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据P来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回。
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数。
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能。
(4)、神经网络多模优化模块,用于对步骤(3)中的每个子网络进行构建
O(·)为模型输出,i为总的子网络数目,x为输入变量,Fj(·)为第j个子网络的输出;即最终多模RBF神经网络的预报结果为各个子网络预报结果的平均值。
在PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
2.丙烯聚合生产过程流程图如图3所示,根据反应机理以及流程工艺分析,考虑到聚丙烯生产过程中对熔融指数产生影响的各个因素,取实际生产过程中常用的九个操作变量和易测变量作为模型输入变量,有:三股丙烯进料流率,主催化剂流率,辅催化剂流率,釜内温度、压强、液位,釜内氢气体积浓度。
表1列出了作为多模态全局最优软测量系统5输入的9个模型输入变量,分别为釜内温度(T)、釜内压强(p)、釜内液位(L)、釜内氢气体积浓度(Xv)、3股丙烯进料流率(第一股丙烯进料流率f1,第二股丙烯进料流率f2,第三股丙烯进料流率f3)、2股催化剂进料流率(主催化剂流率f4,辅催化剂流率f5)。反应釜中的聚合反应是反应物料反复混合后参与反应的,因此模型输入变量涉及物料的过程变量采用前若干时刻的平均值。此例中数据采用前一小时的平均值。熔融指数离线化验值作为多模态全局最优软测量系统5的输出变量。通过人工取样、离线化验分析获得,每4小时分析采集一次。
表1多模态全局最优软测量系统所需模型输入变量
现场智能仪表2及控制站3与丙烯聚合生产过程1相连,与DCS数据库4相连;最优软测量系统5与DCS数据库4及软测量值显示仪6相连。现场智能仪表2测量丙烯聚合生产对象的易测变量,将易测变量传输到DCS数据库4;控制站3控制丙烯聚合生产对象的操作变量,将操作变量传输到DCS数据库4。DCS数据库4中记录的变量数据作为多模态全局最优软测量系统5的输入,软测量值显示仪6用于显示多模态全局最优软测量系统5的输出,即软测量值。
多模态全局最优软测量系统5,包括:
(1)数据预处理模块7,用于对模型输入进行预处理,即中心化和归一化。对输入变量中心化,就是减去变量的平均值,使变量为零均值的变量,从而简化算法;对输入变量归一化,就是除以输入变量值的变化区间,是变量的值落到-0.5~0.5之内,进一步简化。
(2)PCA主成分分析模块8,用于对输入变量预白化处理即变量去相关,对输入变量施加一个线性变换,使得变换后的变量各个分量间互不相关,同时其协方差矩阵为单位阵,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵。若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵。
(3)神经网络模型模块9,采用RBF神经网络,多层前馈神经网络在网络结构上通常由输入层、隐含层和输出层组成。在网络特征上主要表现为既无层内神经元的互联,也无层间的反联络。这种网络实质上是一种静态网络,其输出只是现行输入的函数,而于过去和将来的输入或输出无关。RBF神经网络模型有一个输入层、一个输出层和一个隐藏层。理论上可以证明,RBF神经网络可以任意逼近非线性系统。RBF神经网络训练算法通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;需要建立若干子神经网络,第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子网络预报结果,d(·)为实际结果。
从第二个子网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子网络的训练目标,Fi(·)为第i个网络的预报结果;d(·)为实际结果;F(·)为前i-1个子网络的综合结果;λ为调节参数,N为样本数目。
训练的终止条件为将得到的新的子网络加入多模神经网络后,网络群的预报误差不再减小。
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(i=1,2,…,n)
n为初始解的个数,sn为第n个初始解,k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据P来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回。
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数。
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能。
(4)、神经网络多模优化模块10,用于对步骤(3)中的每个子网络进行构建
O(·)为模型输出,i为总的子网络数目,x为输入变量,Fj(·)为第j个子网络的输出;即最终多模RBF神经网络的预报结果为各个子网络预报结果的平均值。
在PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
(5)模型更新模块11,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型。
实施例2
1.参照图1、图2和图3,一种多模态全局最优丙烯聚合生产过程最优软测量方法包括以下步骤:
(1)对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量由DCS数据库获得;
(2)对样本数据进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(3)PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵。若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置。当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(4)基于模型输入、输出数据建立初始神经网络模型,采用RBF神经网络,通过误差最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子网络预报结果,d(·)为实际结果。
从第二个子网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子网络的训练目标,Fi(·)为第i个网络的预报结果;d(·)为实际结果;F(·)为前i-1个子网络的综合结果;λ为调节参数,N为样本数目。
训练的终止条件为将得到的新的子网络加入多模神经网络后,网络群的预报误差不再减小。
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(i=1,2,…,n)
n为初始解的个数,sn为第n个初始解,k为迭代次数。初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据P来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回。
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数。
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能。
(5)、多模所有的子神经网络,用于对步骤(4)中的每个子网络进行构建
O(·)为模型输出,i为总的子网络数目,x为输入变量,Fj(·)为第j个子网络的输出;即最终多模RBF神经网络的预报结果为各个子网络预报结果的平均值。
进一步,在所述的步骤(3)中采用PCA主成分分析方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
2.本实施例的方法具体实施步骤如下:
步骤1:对丙烯聚合生产过程对象1,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入。
步骤2:对样本数据进行预处理,由数据预处理模块7完成。
步骤3:对经过预处理的数据进行主成分分析,由PCA主成分分析模块8完成。
步骤4:模块9基于模型输入、输出结合步骤(4)建立若干初始神经网络模型。输入数据如步骤1所述获得,输出数据由离线化验获得。
步骤5:模块10结合步骤(5)根据子网络预报误差将所有的子神经网络的构建起来;
步骤6:模型更新模块11定期将离线化验数据输入到训练集中,更新神经网络模型,多模态全局最优软测量系统5建立完成。
步骤7:建立好的多模态全局最优软测量系统5基于DCS数据库4传来的实时模型输入变量数据对丙烯聚合生产过程1的熔融指数进行多模态全局最优软测量。
步骤8:熔融指数软测量显示仪6显示多模态全局最优软测量系统5的输出,完成对丙烯聚合生产过程熔融指数的最优软测量的显示。
Claims (2)
1.一种多模态全局最优丙烯聚合生产过程最优软测量系统,包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库、多模态全局最优软测量系统以及熔融指数软测量显示仪,所述现场智能仪表、控制站与丙烯聚合生产过程连接,所述现场智能仪表、控制站与DCS数据库连接,所述DCS数据库与多模态全局最优软测量系统的输入端连接,所述多模态全局最优软测量系统的输出端与熔融指数软测量显示仪连接,其特征在于:所述多模态全局最优软测量系统包括:
(1)、数据预处理模块,用于将从DCS数据库输入的模型输入变量进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(2)、PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵;若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置;当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(3)、神经网络模型模块,用于采用RBF神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;需要建立若干子神经网络,第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子RBF网络预报结果,d(·)为实际结果;
从第二个子RBF网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子RBF网络的训练目标,Fi(·)为第i个子RBF网络的预报结果;d(·)为实际结果;F(·)为前i-1个子RBF网络的综合结果;λ为调节参数,N为样本数目;
训练的终止条件为将得到的新的子RBF网络加入多模神经网络后,网络群的预报误差不再减小;
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(k),i=1,2,…,n,
n为初始解的个数,k为迭代次数;初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据Pi(k)来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回;
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优次数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数;
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能;
(4)、神经网络多模优化模块,用于对步骤(3)中的每个子RBF网络进行构建:
O(x)为模型输出,i为总的子RBF网络数目,x为输入变量,Fk(x)为第k个子RBF网络的输出,wk为第k个子RBF网络的输出Fk(x)的权值;即最终多模RBF神经网络的预报结果为各个子RBF网络预报结果的平均值;
所述多模态全局最优软测量系统还包括:模型更新模块,用于模型的在线更新,定期将离线化验数据输入到训练集中,更新神经网络模型;
所述多模态全局最优软测量系统训练子RBF神经网络,然后将其的构建起来形成神经网络群;由于子RBF网络的选取标准是预报误差小、与其他的子RBF网络差异大,所以这些预报效果好、又各不相同的子神经网络的综合预报效果能够具有更好的预报精度和稳定性;PCA主成分分析模块中,PCA方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
2.一种用如权利要求1所述的多模态全局最优聚丙烯生产过程最优软测量系统实现的软测量方法,其特征在于:所述软测量方法具体实现步骤如下:
(5.1)对丙烯聚合生产过程对象,根据工艺分析和操作分析,选择操作变量和易测变量作为模型的输入,操作变量和易测变量取温度、压力、液位、氢气气相百分数、3股丙烯进料流速和2股催化剂进料流速这些变量,由DCS数据库获得;
(5.2)对样本数据进行预处理,对输入变量中心化,即减去变量的平均值;再进行归一化处理,即除以变量值的变化区间;
(5.3)PCA主成分分析模块,用于将输入变量预白化处理及变量去相关,通过对输入变量施加一个线性变换实现,即主成分由C=xU得到,其中x为输入变量,C为主成分得分矩阵,U为载荷矩阵;若对原始数据进行重构,可由x=CUT计算,其中上标T表示矩阵的转置;当选取的主成分数目小于输入变量的变量个数时,x=CUT+E,其中E为残差矩阵;
(5.4)基于模型输入、输出数据建立若干个初始子神经网络模型,采用RBF神经网络,通过误差最小化来完成输入到输出的一种高度非线性映射,映射中保持拓扑不变性;第一个子RBF网络的训练目标是预报结果与实际结果差距J1最小;
N为样本数目,x为输入变量,l为样本点序号,F1(·)为子RBF网络预报结果,d(·)为实际结果;
从第二个子RBF网络开始,训练目标变为使得网络的预报误差尽可能小,同时网络的预报结果与之前的网络预报结果又尽可能大的差异,目标函数如下:
Ji为前i个子RBF网络的训练目标,Fi(·)为第i个子RBF网络的预报结果;d(·)为实际结果;F(·)为前i-1个子RBF网络的综合结果;λ为调节参数,N为样本数目;
训练的终止条件为将得到的新的子RBF网络加入多模神经网络后,网络群的预报误差不再减小;
采用一种智能连续空间蚁群算法对每个RBF网络进行训练和最优化,具体步骤为:
(a)算法初始化,根据待优化的RBF神经网络结构构造出初始的解集S=(s1,s2,…,sn),n为初始解的个数,sn为第n个初始解,确定蚁群的大小M,设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen=0;
(b)计算出解集S对应的适应度值Gi(i=1,2,…,n),适应度值越大代表解越好;再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率Pi(k),i=1,2,…,n,
n为初始解的个数,k为迭代次数;初始化执行寻优算法的蚂蚁编号a=0;
(c)蚂蚁a选取S中的一个解作为寻优的初始解,选取规则是根据Pi(k)来做轮盘选;
(d)蚂蚁a在选取的初始解的基础上进行寻优,找到更好的解sa';
(e)如果a<M,则a=a+1,返回步骤c;否则继续向下执行步骤f;
(f)如果gen<MaxGen,则gen=gen+1,使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解,返回步骤b;否则向下执行步骤g;
(g)计算出解集S对应的适应度值Ga(a=1,2,…,n),选取适应度值最大的解作为算法的最优解,结束算法并返回;
每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数,如果本次循环得到了更好的解,则在下次循环中会基于该解并保持搜索方向不变;否则在下次循环中仍基于原来的解但会调整搜索方向;
同时随着整个蚁群寻优代数的增加,蚂蚁搜索的步长会智能地减小,以适合整个蚁群寻优的收敛:
delk=Random·kr(4)
式中,delk为蚂蚁第k代迭代的初始步长,k为迭代代数,Random为随机向量,r为负常实数;
对于解集S中长期不被蚂蚁选作寻优初始解的那些解,会采用遗传算法中的变异和交叉策略进行处理,从而提高算法的全局寻优性能;
(5.5)神经网络多模优化模块,用于对步骤(5.4)中的每个子RBF网络进行构建:
O(x)为模型输出,i为总的子RBF网络数目,x为输入变量,Fk(x)为第k个子RBF网络的输出,wk为第k个子RBF网络的输出Fk(x)的权值;即最终多模RBF神经网络的预报结果为各个子RBF网络预报结果的平均值;
所述软测量方法还包括:定期将离线化验数据输入到训练集中,更新神经网络模型;所述软测量方法训练子RBF神经网络,然后将其构建起来形成神经网络群;由于子RBF网络的选取标准是预报误差小、与其他的子RBF网络差异大,所以这些预报效果好、又各不相同的子神经网络的综合预报效果能够具有更好的预报精度和稳定性;在所述的步骤(5.3)中采用PCA主成分分析方法实现输入变量的预白化处理,能够简化神经网络模型的输入变量,进而提高模型的性能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310659204.1A CN103838143B (zh) | 2013-12-09 | 2013-12-09 | 多模态全局最优丙烯聚合生产过程最优软测量系统和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310659204.1A CN103838143B (zh) | 2013-12-09 | 2013-12-09 | 多模态全局最优丙烯聚合生产过程最优软测量系统和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103838143A CN103838143A (zh) | 2014-06-04 |
CN103838143B true CN103838143B (zh) | 2016-03-30 |
Family
ID=50801774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310659204.1A Expired - Fee Related CN103838143B (zh) | 2013-12-09 | 2013-12-09 | 多模态全局最优丙烯聚合生产过程最优软测量系统和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103838143B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108873838A (zh) * | 2018-06-28 | 2018-11-23 | 浙江大学 | 一种群智寻优的丙烯聚合生产过程最优软测量仪表 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59180348A (ja) * | 1983-03-30 | 1984-10-13 | Mitsubishi Petrochem Co Ltd | 自動メルトインデクス測定装置 |
US5110214A (en) * | 1990-07-09 | 1992-05-05 | Phillips Petroleum Company | Apparatus and method for evaluating the propensity of polymers for smoking during processing |
CN1916791A (zh) * | 2006-09-12 | 2007-02-21 | 浙江大学 | 一种丙烯聚合工业生产熔融指数软测量方法 |
CN101315556A (zh) * | 2008-06-25 | 2008-12-03 | 浙江大学 | 基于Chaos混沌寻优的丙烯聚合生产过程最优软测量仪表及方法 |
CN101901006A (zh) * | 2010-06-30 | 2010-12-01 | 浙江大学 | 一种双层优化的工业过程最优控制系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6544193B2 (en) * | 1996-09-04 | 2003-04-08 | Marcio Marc Abreu | Noninvasive measurement of chemical substances |
-
2013
- 2013-12-09 CN CN201310659204.1A patent/CN103838143B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59180348A (ja) * | 1983-03-30 | 1984-10-13 | Mitsubishi Petrochem Co Ltd | 自動メルトインデクス測定装置 |
US5110214A (en) * | 1990-07-09 | 1992-05-05 | Phillips Petroleum Company | Apparatus and method for evaluating the propensity of polymers for smoking during processing |
CN1916791A (zh) * | 2006-09-12 | 2007-02-21 | 浙江大学 | 一种丙烯聚合工业生产熔融指数软测量方法 |
CN101315556A (zh) * | 2008-06-25 | 2008-12-03 | 浙江大学 | 基于Chaos混沌寻优的丙烯聚合生产过程最优软测量仪表及方法 |
CN101901006A (zh) * | 2010-06-30 | 2010-12-01 | 浙江大学 | 一种双层优化的工业过程最优控制系统及方法 |
Non-Patent Citations (2)
Title |
---|
一种基于新型蚁群算法的聚丙烯熔融指数预报模型;张志猛等;《化工学报》;20110831;第62卷(第8期);第2270-2274页 * |
群智能优化LSSVM最优聚丙烯熔融指数预报;蒋华琴等;《化工学报》;20120930;第63卷(第9期);第2794-2798页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103838143A (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101315557B (zh) | 基于遗传算法优化bp神经网络的丙烯聚合生产过程最优软测量仪表及方法 | |
CN101382801B (zh) | 基于ega优化的丙烯聚合生产过程最优软测量仪表及方法 | |
CN101315556B (zh) | 基于Chaos混沌寻优的丙烯聚合生产过程最优软测量仪表及方法 | |
CN103823964A (zh) | 一种全局最优丙烯聚合生产过程最优软测量仪表和方法 | |
CN103823430B (zh) | 智能加权丙烯聚合生产过程最优软测量系统和方法 | |
CN103839103B (zh) | 丙烯聚合生产过程bp最优预报系统和方法 | |
CN109886464B (zh) | 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法 | |
CN103942461A (zh) | 基于在线贯序极限学习机的水质参数预测方法 | |
CN110457369A (zh) | 一种模型的训练方法及相关设备 | |
CN103838142B (zh) | 基于混合寻优的丙烯聚合生产过程最优软测量系统和方法 | |
Hernández-Lobato et al. | Designing neural network hardware accelerators with decoupled objective evaluations | |
CN103838206A (zh) | 最优bp多模丙烯聚合生产过程最优软测量仪表和方法 | |
CN103824121A (zh) | 基于多模群智的丙烯聚合生产过程最优预报系统和方法 | |
CN103955170A (zh) | 基于群智寻优的丙烯聚合生产过程在线预报系统和方法 | |
CN109858785A (zh) | 一种评价智能电能表运行状态的方法和系统 | |
CN103675006A (zh) | 最小二乘的工业熔融指数软测量仪表和方法 | |
CN103838955A (zh) | 最优混合丙烯聚合生产过程最优软测量仪表和方法 | |
CN103823369B (zh) | 基于bp多模网络的丙烯聚合生产过程预报系统和方法 | |
CN103838205A (zh) | Bp全局最优丙烯聚合生产过程最优软测量仪表和方法 | |
CN103838143B (zh) | 多模态全局最优丙烯聚合生产过程最优软测量系统和方法 | |
CN103838209A (zh) | 丙烯聚合生产过程自适应最优预报系统和方法 | |
CN103838957A (zh) | 丙烯聚合生产过程径向基最优软测量仪表和方法 | |
CN103838958A (zh) | 模糊智能最优丙烯聚合生产过程最优软测量仪表和方法 | |
CN109829197A (zh) | 改进的布谷鸟寻优的丙烯聚合生产过程最优软测量仪表 | |
CN109856971A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160330 Termination date: 20211209 |