CN109472397B - 基于粘度变化的聚合工艺参数调节方法 - Google Patents
基于粘度变化的聚合工艺参数调节方法 Download PDFInfo
- Publication number
- CN109472397B CN109472397B CN201811221084.6A CN201811221084A CN109472397B CN 109472397 B CN109472397 B CN 109472397B CN 201811221084 A CN201811221084 A CN 201811221084A CN 109472397 B CN109472397 B CN 109472397B
- Authority
- CN
- China
- Prior art keywords
- elm
- particle
- value
- polymerization process
- follows
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000006116 polymerization reaction Methods 0.000 title claims abstract description 48
- 230000008859 change Effects 0.000 title claims abstract description 26
- 239000002245 particle Substances 0.000 claims abstract description 81
- 239000000155 melt Substances 0.000 claims abstract description 32
- 238000012549 training Methods 0.000 claims abstract description 30
- 238000005469 granulation Methods 0.000 claims abstract description 27
- 230000003179 granulation Effects 0.000 claims abstract description 27
- 230000006870 function Effects 0.000 claims abstract description 23
- 230000013016 learning Effects 0.000 claims abstract description 20
- 239000013598 vector Substances 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 41
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 230000004913 activation Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 230000006872 improvement Effects 0.000 claims description 2
- 229920000642 polymer Polymers 0.000 claims description 2
- 230000002776 aggregation Effects 0.000 abstract 1
- 238000004220 aggregation Methods 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000007380 fibre production Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 238000006068 polycondensation reaction Methods 0.000 description 2
- 229920000728 polyester Polymers 0.000 description 2
- 230000009326 social learning Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013485 heteroscedasticity test Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 239000004753 textile Substances 0.000 description 1
- 230000037303 wrinkles Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于粘度变化的聚合工艺参数调节方法，将过去连续时间段Ta内的熔体粘度数据转化为时间序列I后对其进行模糊信息粒化处理得到三个特征向量Low、R和Up，再将三个特征向量输入到改进的ELM中由其输出与未来连续时间段Tb内的熔体粘度数据的波动范围相关的数值，将波动范围与参考范围进行比较后根据比较结果调节聚合工艺参数，改进的ELM是采用粒子群算法优化输入层权值和隐含层阈值后的ELM，粒子群算法的适应度函数为训练样本的均方误差。本发明将模糊信息粒化、极限学习机及粒子群算法有机结合，样本维数少，实时性好，预测精度高；同时预测结果可对工艺参数进行调节。
Description
技术领域
本发明属于工艺优化领域，涉及一种基于粘度变化的聚合工艺参数调节方法，特别是涉及一种基于模糊信息粒化和改进的ELM的预测熔体特性粘度进而调节聚合工艺参数的方法。
背景技术
聚酯纤维不仅弹性高、强度大，还具有良好的化学稳定性、耐磨性、拉伸性和抗折皱性，因此在民用纺织品和工业纤维品领域都有非常广泛的应用。聚合过程是聚酯纤维生产的主体和龙头，聚合过程的水平会直接影响产品的数量和质量。其中，熔体的特性粘度是直接体现聚合水平的指标。因此，对特性粘度的预测能为后续的生产环节提供一定的指导。
目前，国内外相关的理论研究主要是利用反应釜的温度、压强及流量等容易测得的变量作为输入，将特性粘度作为输出，利用支持向量机或人工神经网络等方法进行静态预测。然而实际生产数据是动态的，其具有一定的时序性，因而静态预测的准确率不高。所以针对特性粘度构成的时间序列进行分析是有必要的。粘度变化预测侧重于对粘度的走向进行预测，而并不局限于某个时间点的粘度值的精确预测。如生产人员能了解聚合过程中熔体的性能指标的变化，对指导生产是十分有帮助的，可以基于粘度变化相应调整工艺参数以优化其产品品质。因此，进行粘度变化的预测就显得更有意义。而现阶段关于这方面的研究在国内外还是空白。
信息粒化这一概念是由Fuzzy集创始人L.A.Zadeh在1979年首次提出的。L.A.Zadeh认为人类的认识和推理是由三个基本概念构成：粒化、组织和因果。粒化是将整体分解成部分，组织则是综合部分成为整体，因果是指因果关系。目前，国内外主要有三种信息粒化模型：基于模糊集理论的模型、基于粗糙集理论的模型和基于商空间理论的模型。用模糊集对时间序列进行模糊粒化，主要可以分为两个步骤：窗口划分和模糊化。窗口划分就是将所给的时间序列分割成一个个的小子列，作为一个个的操作窗口；而模糊化则是将第一步产生的每一个窗口进行模糊化，生成一个个模糊集，也就是模糊粒子。这两种广义模式结合在一起就是模糊信息粒化。模糊信息粒化可以将时间序列缩小为窗口序列，应用于对变化范围的预测。
传统的时序预测的方法是根据历史趋势预测未来时间序列的趋势，即通过建立相应的数学模型来拟合历史时间趋势曲线，并根据模型趋势曲线预测未来时间序列，常用的模型有自回归和移动平均模式、自回归综合移动平均模型、向量自回归模型、阈值自回归模型和自回归条件异方差模型等，但是它们只适用于处理线性或简单的非线性数据，无法应用于复杂非线性的工业生产数据；近年来随着深度学习的兴起，深层神经网络被广泛用于时序预测，但是它们的时间代价太高，硬件要求比较严格，不符合工业要求的实时性。
极限学习机是一种针对单隐含层前馈神经网络的新算法。极限学习机克服了传统的前馈神经网络训练速度慢、容易陷入局部极小值点和学习率参数的选择敏感的缺点，极限学习机随机产生输入层与隐含层的连接权值及隐含层神经元的阈值，且在训练过程中无需调整，只需要设置隐含层神经元的个数，其具有易于实现、速度快和泛化能力强的优点。然而，采用随机生成初始化参数的方式促使传统极限学习机算法模型不可避免地带来了隐含层节点冗余和预测结果波动性小的缺点，阻碍了它在高精度工程和高指标实验上的应用。因此，需要选择一种智能算法对输入层权重和隐含层阈值进行优化，提高模型的预测精度。
粒子群算法是一种基于群体智能的优化算法。该算法源于对鸟类捕食行为的研究，当鸟类捕食时，每只鸟找到食物最简单有效地方法就是搜寻当前距离食物最近的鸟的周围区域。其优点在于：(1)粒子群算法没有交叉和变异运算，依靠粒子速度完成搜索，并且在迭代进化中只有最优的粒子把信息传递给其它粒子，搜索速度快；(2)粒子群算法具有记忆性，粒子群体的历史最好位置可以记忆并传递给其它粒子；(3)需要调整的参数较少，结构简单，易于工程实现；(4)采用实数编码，直接由问题的解决定，问题解的变量数目直接作为粒子的维数。
因此，开发一种基于模糊信息粒化和改进的ELM的预测熔体特性粘度进而调节聚合工艺参数的方法极具现实意义。
发明内容
本发明的目的是克服现有技术的不足，提供一种基于模糊信息粒化和改进的ELM的预测熔体特性粘度进而调节聚合工艺参数的方法。本发明的方法将模糊信息粒化方法、粒子群算法及极限学习机有机结合起来，首先通过对采集到的特性粘度数据进行模糊信息粒化构成窗口序列，实现了区间预测；然后通过采用粒子群算法对极限学习机的输入层权值和隐含层阈值进行优化，提高了预测精度和速度；最后基于预测得到的熔体粘度调节聚合工艺参数。
为了实现上述目的，本发明采用如下技术方案：
基于粘度变化的聚合工艺参数调节方法，将由过去连续时间段Ta内的熔体粘度数据预测得到的未来连续时间段Tb内的熔体粘度数据的波动范围与参考范围进行比较后根据比较结果调节聚合工艺参数；
所述预测的方法为：将过去连续时间段Ta内的熔体粘度数据转化为时间序列I后对其进行模糊信息粒化处理得到三个特征向量，再将三个特征向量输入到改进ELM中由其输出三个代表未来连续时间段Tb内的熔体粘度数据最大值、平均值和最小值的数值，即得到未来连续时间段Tb内的熔体粘度数据的波动范围；所述未来连续时间段Tb的长度等于对时间序列I进行模糊信息粒化处理时窗口的大小；
所述改进的ELM是采用粒子群算法优化输入层权值和隐含层阈值后的ELM，所述粒子群算法的适应度函数为训练样本的均方误差；
所述训练样本的构建过程为：首先采集过去连续时间段Tc内的熔体粘度数据构成时间序列II，然后对时间序列II进行模糊信息粒化处理得到三个特征向量，最后分别根据三个特征向量构建输入矩阵P和输出矩阵T得到样本集(P,T)，再由样本集得到训练样本。输入矩阵P和输出矩阵T的表达式如下：
T＝[En+1 En+2…En+1 …E[M/W]]；
式中，El代表特征向量中第l个元素，l＝1,2,...,[M/w]-n，En+l代表特征向量中第n+l个元素，[M/w]代表特征向量中元素的总个数，输入矩阵P中的第l列元素与输出矩阵T中的第l列元素构成训练样本；
所述参考范围是已经生产出来的、性能满足使用要求的聚合物在连续时间段Tb内的熔体粘度数据的波动范围；
所述根据比较结果调节聚合工艺参数是指根据生产经验获得聚合工艺参数与熔体粘度的相关性后调节聚合工艺参数提高或降低熔体粘度进而降低波动范围与参考范围的偏差。
作为优选的技术方案：
如上所述的基于粘度变化的聚合工艺参数调节方法，所述熔体粘度为特性粘度，单位为dl/g；所述聚合工艺参数为终缩聚反应釜内的温度、压强和传质系数，单位分别为℃、MPa和m/s。
如上所述的基于粘度变化的聚合工艺参数调节方法，所述过去连续时间段Ta和过去连续时间段Tc的长度为M，M大于等于15min。过去连续时间段Ta和过去连续时间段Tc的熔体粘度数据来源于某纤维生产工艺的聚合过程，按照每分钟的采样频率，从分散控制系统(Distributed Control System,DCS)中收集得到的。
如上所述的基于粘度变化的聚合工艺参数调节方法，所有的时间序列记为矩阵X，表达式为X＝[x1,x2,...,xt,...,xM]T，其中，xt为第t分钟的粘度，t＝1,2,...,M。
如上所述的基于粘度变化的聚合工艺参数调节方法，所有的模糊信息粒化处理的步骤如下：
(1)划分窗口；
确定粒化窗口的大小为w，将X以每w个数据为单位长度划分为[M/w]个子列窗口，记为ΔXi，i＝1,2,...,[M/w]，[M/w]为M/w向前取整数，即如M/w的结果为整数，则[M/w]为M/w，反之，则为(M/w)+1，(M/w)代表M/w的整数部分；
(2)建立模糊粒子；
分别在每个子列窗口ΔXi上建立模糊粒子Ai；
(3)确定模糊粒子的参数；
分别确定每个模糊粒子Ai的支撑下限ai、支撑上限bi和核mi得到三个特征向量Low、R和Up；
ai和bi的计算公式如下：
mi和mi'的取值方式为：将ΔXi＝(xi1,xi2,...,xij,...xid)按从小到大重新排序，当d为偶数时，mi＝xi(d/2)，mi'＝xi[(d+2)/2]，当d为奇数时，mi＝xi[(d+1)/2]，mi'＝mi，记mi＝x(e)，mi'＝x(e')，其中d表示第i个子序列中元素的个数，xij表示第i个子序列中第j个元素的值，e和e'分别表示m和m'取值时对应的下标；
Low、R和Up的表达式如下：
Low＝[a1,a2,...,a[M/w]]T
R＝[m1,m2,...,m[M/w]]T；
Up＝[b1,b2,...,b[M/w]]T
其中，Low、R和Up参数描述的分别是X经过模糊信息粒化后在每个子列窗口上变化的最小水平、平均水平和最大水平。
如上所述的基于粘度变化的聚合工艺参数调节方法，所述模糊粒子的形式为三角型、高斯型或抛物线型，本发明中模糊粒子A的形式包括但不限于三角型、高斯型和抛物线型，优选为三角型，因为采用三角型作为模糊粒子A的形式能够得知窗口的最小值、平均值和最大值，有利于实现对特性粘度波动趋势和范围的预测。
如上所述的基于粘度变化的聚合工艺参数调节方法，所述模糊粒子的形式为三角型，其隶属函数如下：
式中，xi是特性粘度的值。
如上所述的基于粘度变化的聚合工艺参数调节方法，所述训练样本还进行归一化处理，归一化公式如下：
xrs'＝(k max-k min)*(xrs-xr min)/(xr max-xr min)+k min；
yrs'＝(k max-k min)*(yrs-yr min)/(yr max-yr min)+k min；
其中，k max和k min为归一化区间[-1,1]的最大值和最小值，xrs为输入矩阵P的第r行第s列的元素，xr min和xr max分别是输入矩阵P的第r行元素的最小值和最大值，xrs'为归一化处理后的xrs，yrs为输出矩阵T的第r行第s列的元素，yr min和yr max分别是输出矩阵T的第r行元素的最小值和最大值，yrs'为归一化处理后的yrs。
如上所述的基于粘度变化的聚合工艺参数调节方法，所述ELM的改进步骤如下：
(1)粒子群算法初始化；
设定种群规模、最大迭代次数、最大惯性权值、最小惯性权值、学习因子、粒子的最大速度、粒子的最小速度、最大位置、最小位置和粒子的维度D，D＝(I+1)×K+(K+1)×O，I为ELM输入层节点的个数，K为ELM隐含层节点的个数，O为ELM输出层节点的个数；
(2)选择适应度函数；
种群中的每一个粒子都是由ELM的输入层权值矩阵和隐含层阈值矩阵两部分组成，可以根据种群的信息，采用ELM对训练样本进行计算，以训练样本预测结果的均方误差作为适应度函数f，适应度函数f的表达形式如下：
其中，N为训练样本的个数，Tu为训练样本输入到ELM后ELM的实际输出值，
为训练样本输入到ELM后ELM的理论输出值，适应度值为函数值，在训练过程中应使适应度值越小越好；
(3)初始化迭代次数z＝1；
(4)寻找个体极值Pub和群体极值Pg；
每个粒子根据预设的适应度函数计算各自的适应度值，并和当前的个体极值和群体极值进行比较，在Pub和Pg的选择中，具有较小的均方误差的粒子被选中，其计算公式如下：
其中，
和
分别是第u个粒子的均方误差、第u个粒子最优的均方误差和所有粒子中最优的均方误差；
(5)更新种群的惯性权值ω；
惯性权值可以调节对解空间的搜索能力，采用更新公式使得惯性权值随着迭代次数的增加而减小，更新公式如下：
其中，ωmax和ωmin分别为初始设置的最大惯性权值和最小惯性权值，z和T分别为当前的迭代次数和总的迭代次数；
(6)更新粒子的位置和速度向量；
进行更新时，组成粒子的所有元素的取值都要限制在区间[-1,1]内，粒子的速度和位置的更新规则如下：
其中，c1和c2是学习因子，分别代表了粒子的自身学习能力和社会学习能力，rand1和rand2是区间[0,1]内的随机数，Vu(z)和Xu(z)分别代表当前时刻粒子的速度和位置，Vu(z+1)和Xu(z+1)分别代表下一时刻粒子的速度和位置；
(7)判断当前迭代次数z是否达到最大迭代次数，如果是则进入步骤(9)，反之，进入步骤(8)；
(8)令z＝z+1，返回步骤(4)；
(9)得到最优的ELM的输入层权值和隐含层偏置存在于群体极值pg中，并通过Moore-Penrose广义逆公式β＝H+T求解输出层权值，其中，β为ELM的输出层权值，T为ELM的输出值，H+为隐含层输出矩阵H的广义逆计算，H＝g(AX+B)，A为输入层权值，B为隐含层阈值，X为样本的输入矩阵，g(·)为隐含层的激活函数。
有益效果：
(1)本发明的基于粘度变化的聚合工艺参数调节方法，采用模糊信息粒化能够降低样本维数，有效地将原始数据转换成窗口序列，为预测提供了样本，减小了数据处理量；
(2)本发明的基于粘度变化的聚合工艺参数调节方法，采用极限学习机进行预测计算，符合工业数据的复杂非线性和实时性要求，且耗时短；
(3)本发明的基于粘度变化的聚合工艺参数调节方法，采用粒子群算法寻优取代极限学习机权值和阈值随机赋值，提高了预测的稳定性、精确度和泛化性能；
(4)本发明的基于粘度变化的聚合工艺参数调节方法，填补了特性粘度预测在变化趋势和空间范围这一方面的研究空白，根据特性粘度的预测结果可对聚合工艺参数进行调节以保证产品质量，极具应用前景。
附图说明
图1为本发明的基于粘度变化的聚合工艺参数调节方法中的预测流程图；
图2为本发明采用粒子群算法优化极限学习机的流程图；
图3、4和5分别为本发明Low序列、R序列和Up序列的预测值与真实值的对比示意图。
具体实施方式
下面结合具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
一种基于粘度变化的聚合工艺参数调节方法，其步骤如下：
(1)获取过去连续时间段Ta内和过去连续时间段Tc内的熔体粘度数据，其中过去连续时间段Ta和过去连续时间段Tc的长度为M，M为5000min，熔体粘度为特性粘度，单位为dl/g；
(2)分别将过去连续时间段Ta内和过去连续时间段Tc内的熔体粘度数据构成时间序列，其中时间序列构成方法如下：
将时间序列记为矩阵X，表达式为X＝[x1,x2,...,xt,...,xM]T，其中，xt为第t分钟的粘度，t＝1,2,...,M；
(3)分别对Ta和Tc的时间序列进行模糊信息粒化处理得到三个特征向量Low、R和Up，步骤如下：
(3.1)划分窗口；
确定粒化窗口的大小为w，w＝5，将X以每5个数据为单位长度划分为1000个子列窗口，记为ΔXi，i＝1,2,...,1000；
(3.2)建立模糊粒子；
分别在每个子列窗口ΔXi上建立模糊粒子Ai，模糊粒子Ai的形式为三角型，其隶属函数如下：
式中，xi是特性粘度的值，ai为Ai的支撑下限，bi为Ai的支撑上限，mi为Ai的核；
(3.3)确定模糊粒子的参数；
分别确定每个模糊粒子Ai的支撑下限ai、支撑上限bi和核mi得到三个特征向量Low、R和Up；
ai、bi和mi的计算公式如下：
mflag和mflag′的取值方式为：将ΔXi＝(xi1，xi2，...，xij，...xid)按从大到小重新排序，
其中
代表向上取整函数，
代表向下取整函数。d表示第i个子序列中元素的个数，xij表示第i个子序列中第j个元素的值。
Low、R和Up的表达式如下：
Low＝[a1，a2，...，a1000]T
R＝[m1，m2，...，m1000]T
Up＝[b1，b2，...，b1000]T；
其中，Low、R和Up参数分别描述的是X经过模糊信息粒化后在每个子列窗口上变化的最小水平、平均水平和最大水平；
(4)构建训练样本及测试样本：
(4.1)采取单步滚动预测的方式，针对Tc对应的三个特征向量Low、R和Up，分别用前3个窗口的数据预测后面1个窗口的数据，以最小值序列Low＝[a1，a2，...，a1000]T为例，利用其构建得到的输入矩阵P和输出矩阵T表达形式如下：
T＝[a4 a5 a6 … a1000]；
式中，矩阵P的第i列与矩阵T的第i列构成第i个样本，i＝1,2,...,997，P与T构成样本集；
(4.2)从样本集中选取75％的样本作为训练样本，余下25％作为测试样本；
(4.3)分别利用R和Up序列构建训练样本和测试样本，步骤同(4.1)和(4.2)；
(4.4)对所有的样本进行归一化处理，其中归一化公式如下：
xrs'＝(k max-k min)*(xrs-xr min)/(xr max-xr min)+k min；
yrs'＝(k max-k min)*(yrs-yr min)/(yr max-yr min)+k min；
其中，k max和k min为归一化区间[-1,1]的最大值和最小值，xrs为输入矩阵P的第r行第s列的元素，xr min和xr max分别是输入矩阵P的第r行元素的最小值和最大值，xrs'为归一化处理后的xrs，yrs为输出矩阵T的第r行第s列的元素，yr min和yr max分别是输出矩阵T的第r行元素的最小值和最大值，yrs'为归一化处理后的yrs；
(5)改进ELM，改进ELM的步骤如图2所示：
(5.1)粒子群算法初始化；
设定种群规模为100，最大迭代次数为200，最大惯性权值和最小惯性权值分别为0.8和0.4，两个学习因子都为2，限定粒子的最大速度和最小速度分别为-1和1，最大位置和最小位置分别为-1和1，粒子的维度为D＝(I+1)×K+(K+1)×O，其中，I为极限学习机输入层节点的个数，K为隐含层节点的个数，O为输出层节点的个数；
(5.2)选择适应度函数；
种群中的每一个粒子都是由ELM的输入层权值矩阵和隐含层阈值矩阵两部分组成，可以根据种群的信息，采用ELM对训练样本进行计算，以训练样本预测结果的均方误差作为适应度函数f，适应度函数f的表达形式如下：
其中，N为训练样本的个数，Tu为训练样本输入到ELM后ELM的实际输出值，
为训练样本输入到ELM后ELM的理论输出值，适应度值为函数值，在训练过程中应使适应度值越小越好；
(5.3)初始化迭代次数z＝1；
(5.4)寻找个体极值pub和群体极值Pg；
每个粒子根据预设的适应度函数计算各自的适应度值，并和当前的个体极值和群体极值进行比较，在pub和pg的选择中，具有较小的均方误差的粒子被选中，其计算公式如下：
其中，
和
分别是第u个粒子的均方误差、第u个粒子最优的均方误差和所有粒子中最优的均方误差；
(5.5)更新种群的惯性权值；
惯性权值可以调节对解空间的搜索能力，采用更新公式使得惯性权值随着迭代次数的增加而减小，更新公式如下：
其中，ωmax和ωmin分别为初始设置的最大惯性权值和最小惯性权值，z和T分别为当前的迭代次数和总的迭代次数；
(5.6)更新粒子的位置和速度向量；
进行更新时，组成粒子的所有元素的取值都要限制在区间[-1,1]内，粒子的速度和位置的更新规则如下：
其中，ω是惯性权值，c1和c2是学习因子，分别代表了粒子的自身学习能力和社会学习能力，rand1和rand2是区间[0,1]内的随机数，Vu(z)和Xu(z)分别代表当前时刻粒子的速度和位置，Vu(z+1)和Xu(z+1)分别代表下一时刻粒子的速度和位置；
(5.7)判断当前迭代次数z是否达到最大迭代次数，如果是则进入步骤(5.9)，反之，进入步骤(5.8)；
(5.8)令z＝z+1，返回步骤(5.4)；
(5.9)得到最优的ELM的输入层权值和隐含层偏置存在于群体极值pg中，并通过广义逆公式β＝H+T求解输出层权值，其中，β为ELM的输出层权值，T为ELM的输出值，H+为隐含层输出矩阵H的广义逆计算，H＝g(AX+B)，A为输入层权值，B为隐含层阈值，X为样本的输入矩阵，g(·)为隐含层的激活函数；
(6)如图1所示，将步骤(3)得到的Ta对应的Low、R和Up输入到步骤(5)得到改进ELM中由其输出三个代表未来连续时间段Tb内的熔体粘度数据最大值、平均值和最小值的数值，即得到未来连续时间段Tb内的熔体粘度数据的波动范围，未来连续时间段Tb的长度等于Ta的长度；
(7)经未来连续时间段Tb内的熔体粘度数据的波动范围与参考范围进行比较后根据比较结果调节聚合工艺参数，参考范围是已经生产出来的、性能满足使用要求的聚合物在连续时间段Tb内的熔体粘度数据的波动范围，根据比较结果调节聚合工艺参数是指根据生产经验获得聚合工艺参数(终缩聚反应釜内的温度、压强和传质系数，单位分别为℃、MPa和m/s)与熔体粘度的相关性后调节聚合工艺参数提高或降低熔体粘度进而降低波动范围与参考范围的偏差，具体为：根据生产经验可知终缩聚反应釜内的温度和传质系数与熔体粘度正相关，压强与熔体粘度负相关，假设参考范围为[Qa,Za]，波动范围为[Qb,Zb]，依次比较Qa、Qb和Za、Zb的大小。如果Qa≤Qb且Za≥Zb，则说明波动范围在参考范围内，粘度值属于正常情况，无需调整工艺参数；如果Qa≥Qb且Za≥Zb，则说明波动范围超出了参考范围，粘度值低于正常值，则可以采取增大温度、增大传质系数或降低压强等措施；如果Qa≤Qb且Za≤Zb，则说明波动范围超出了参考范围，粘度值高于正常值，则可以采取降低温度、降低传质系数或增大压强等措施；如果Qa≥Qb且Za≤Zb，则说明波动范围超出了参考范围，该段时间内的粘度值既出现了低于正常值的情况，也出现了高于正常值的情况，则可以先采取增大温度、增大传质系数或降低压强等措施消除粘度值低于正常值的情况，再采取增大压强、降低传质系数或降低温度等措施。根据经验，将0.01dl/g设为粘度波动范围与参考范围的偏差单位值，则相应的工艺参数的调整单位值为15℃、0.075MPa和0.04m/s。
为了验证改进ELM的准确性，本发明还将由过去连续时间段Tc对应的三个特征向量Low、R和Up得到的测试样本输入到改进ELM中，由其输出预测值，并结合真实值计算均方误差MSE和相关系数R1，计算公式如下：
其中，N为测试样本的个数，Yi为预测值，
为真实值，
为Y、
的协方差，D(Y)、
分别为Y、
的方差，计算结果如表1所示，其中Low、R和Up序列预测值与真实值的比较示意图如图3、4和5所示。
表1
| 评价指标 | 测试时间 | MSE(均方误差) | R1(相关系数) |
| Low序列 | 0.009s | 0.00690 | 0.88899 |
| R序列 | 0.009s | 0.00522 | 0.91805 |
| Up序列 | 0.009s | 0.00584 | 0.91023 |
| 平均结果 | 0.009s | 0.00597 | 0.90576 |
由图3～图5可以看出，本发明的预测方法所得的预测结果符合特性粘度实际的变化状态(真实值所对应的曲线)，其均方误差小，相关系数较高，结合表1可以发现，本发明的基于粘度变化的聚合工艺参数调节方法对未来时刻粘度的变化预测精确，稳定性好，反应迅速，根据特性粘度的预测结果可对聚合工艺参数进行调节以保证产品质量，极具应用前景。
Claims (5)
1.基于粘度变化的聚合工艺参数调节方法，其特征是：将由过去连续时间段Ta内的熔体粘度数据预测得到的未来连续时间段Tb内的熔体粘度数据的波动范围与参考范围进行比较后，根据比较结果调节聚合工艺参数；
所述预测的方法为：将过去连续时间段Ta内的熔体粘度数据转化为时间序列I后对其进行模糊信息粒化处理得到三个特征向量，再将三个特征向量输入到改进ELM中由其输出三个代表未来连续时间段Tb内的熔体粘度数据最大值、平均值和最小值的数值，即得到未来连续时间段Tb内的熔体粘度数据的波动范围；所述未来连续时间段Tb的长度等于对时间序列I进行模糊信息粒化处理时窗口的大小；
所述改进ELM是采用粒子群算法优化输入层权值和隐含层阈值后的ELM，所述粒子群算法的适应度函数为训练样本的均方误差；
所述ELM的改进步骤如下：
(1)粒子群算法初始化；
设定种群规模、最大迭代次数、最大惯性权值、最小惯性权值、学习因子、粒子的最大速度、粒子的最小速度、最大位置、最小位置和粒子的维度D，D＝(S+1)×K+(K+1)×O，S为ELM输入层节点的个数，K为ELM隐含层节点的个数，O为ELM输出层节点的个数；
(2)选择适应度函数；
适应度函数f的表达形式如下：
其中，N为训练样本的个数，Tu为训练样本输入到ELM后ELM的实际输出值，
为训练样本输入到ELM后ELM的理论输出值；
(3)初始化迭代次数z＝1；
(4)寻找个体极值Pub和群体极值Pg；
计算公式如下：
其中，
和
分别是第u个粒子的均方误差、第u个粒子最优的均方误差和所有粒子中最优的均方误差；
(5)更新种群的惯性权值ω；
更新公式如下：
其中，ωmax和ωmin分别为初始设置的最大惯性权值和最小惯性权值，z和T分别为当前的迭代次数和总的迭代次数；
(6)更新粒子的位置和速度向量；
粒子的速度和位置的更新规则如下：
其中，c1和c2是学习因子，rand1和rand2是区间[0,1]内的随机数，Vu(z)和Xu(z)分别代表当前时刻粒子的速度和位置，Vu(z+1)和Xu(z+1)分别代表下一时刻粒子的速度和位置；
(7)判断当前迭代次数z是否达到最大迭代次数，如果是则进入步骤(9)，反之，进入步骤(8)；
(8)令z＝z+1，返回步骤(4)；
(9)得到最优的ELM的输入层权值和隐含层偏置存在于群体极值pg中，并通过广义逆公式β＝H+T求解输出层权值，其中，β为ELM的输出层权值，T为ELM的输出值，H+为隐含层输出矩阵H的广义逆计算，H＝g(AX+B)，A为输入层权值，B为隐含层阈值，X为样本的输入矩阵，g(·)为隐含层的激活函数；
所述训练样本的构建过程为：首先采集过去连续时间段Tc内的熔体粘度数据构成时间序列II，然后对时间序列II进行模糊信息粒化处理得到三个特征向量，最后分别根据三个特征向量构建输入矩阵P和输出矩阵T得到样本集P,T再由样本集得到训练样本，输入矩阵P和输出矩阵T的表达式如下：
T＝[En+1 En+2...En+l...E[M/W]]
式中，El代表特征向量中第l个元素，l＝1,2,...,[M/w]-n，En+l代表特征向量中第n+l个元素，[M/w]代表特征向量中元素的总个数，输入矩阵P中的第l列元素与输出矩阵T中的第l列元素构成训练样本；
所述的模糊信息粒化处理的步骤如下：
(1)划分窗口；
确定粒化窗口的大小为w，将X以每w个数据为单位长度划分为[M/w]个子列窗口，记为ΔXi，i＝1,2,...,[M/w]，[M/w]为M/w向前取整数，即如M/w的结果为整数，则[M/w]为M/w，反之，则为(M/w)+1，(M/w)代表M/w的整数部分；
(2)建立模糊粒子；
分别在每个子列窗口ΔXi上建立模糊粒子Ai；
(3)确定模糊粒子的参数；
分别确定每个模糊粒子Ai的支撑下限ai、支撑上限bi和核mi得到三个特征向量Low、R和Up；
ai和bi的计算公式如下：
mi和mi'的取值方式为：将ΔXi＝(xi1,xi2,...,xij,...xid)按从小到大重新排序，当d为偶数时，mi＝xi(d/2)，mi'＝xi[(d+2)/2]，当d为奇数时，mi＝xi[(d+1)/2]，mi'＝mi，记mi＝x(e)，mi'＝x(e')，其中d表示第i个子序列中元素的个数，xij表示第i个子序列中第j个元素的值，e和e'分别表示m和m'取值时对应的下标；
Low、R和Up的表达式如下：
所述模糊粒子的形式为三角型，其隶属函数如下：
式中，xi是特性粘度的值；
所述参考范围是已经生产出来的、性能满足使用要求的聚合物在连续时间段Tb内的熔体粘度数据的波动范围；
所述根据比较结果调节聚合工艺参数是指根据生产经验获得聚合工艺参数与熔体粘度的相关性后调节聚合工艺参数提高或降低熔体粘度进而降低波动范围与参考范围的偏差。
2.根据权利要求1所述的基于粘度变化的聚合工艺参数调节方法，其特征在于，所述熔体粘度为特性粘度，单位为dl/g；所述聚合工艺参数为聚合反应釜内的温度、压强和传质系数，单位分别为℃、MPa和m/s。
3.根据权利要求2所述的基于粘度变化的聚合工艺参数调节方法，其特征在于，所述过去连续时间段Ta和过去连续时间段Tc的长度为M，M大于等于15min。
4.根据权利要求3所述的基于粘度变化的聚合工艺参数调节方法，其特征在于，所述的时间序列记为矩阵X，表达式为X＝[x1，x2，...，xt，...，xM]T，其中，xt为第t分钟的粘度，t＝1，2，...，M。
5.根据权利要求1所述的基于粘度变化的聚合工艺参数调节方法，其特征在于，所述训练样本还进行归一化处理，归一化公式如下：
xrs′＝(k max-k min)*(xrs-xrmin)/(xrmax-xrmin)+k min；
yrs′＝(k max-k min)*(yrs-yrmin)/(yrmax-yrmin)+k min；
其中，k max和k min为归一化区间[-1，1]的最大值和最小值，xrs为输入矩阵P的第r行第s列的元素，xrmin和xrmax分别是输入矩阵P的第r行元素的最小值和最大值，xrs′为归一化处理后的xrs，yrs为输出矩阵T的第r行第s列的元素，yrmin和yrmax分别是输出矩阵T的第r行元素的最小值和最大值，yrs'为归一化处理后的yrs。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811221084.6A CN109472397B (zh) | 2018-10-19 | 2018-10-19 | 基于粘度变化的聚合工艺参数调节方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811221084.6A CN109472397B (zh) | 2018-10-19 | 2018-10-19 | 基于粘度变化的聚合工艺参数调节方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109472397A CN109472397A (zh) | 2019-03-15 |
| CN109472397B true CN109472397B (zh) | 2021-12-14 |
Family
ID=65665019
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811221084.6A Active CN109472397B (zh) | 2018-10-19 | 2018-10-19 | 基于粘度变化的聚合工艺参数调节方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109472397B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111861050B (zh) * | 2019-04-25 | 2024-02-20 | 富联精密电子(天津)有限公司 | 生产制程管控装置、方法及计算机可读存储介质 |
| CN110955147A (zh) * | 2019-12-16 | 2020-04-03 | 泉州华中科技大学智能制造研究院 | 一种施胶工艺参数优化方法 |
| CN113239654B (zh) * | 2021-05-21 | 2022-03-25 | 电子科技大学 | 基于fig和ipso算法的igbt剩余寿命预测方法 |
| CN113240096B (zh) * | 2021-06-07 | 2022-03-11 | 北京理工大学 | 基于粗糙集和神经网络的铸造气缸盖微组织预测方法 |
| CN116027756B (zh) * | 2023-02-24 | 2023-07-04 | 季华实验室 | 生产参数在线最优化方法、装置、电子设备及存储介质 |
| CN116501003B (zh) * | 2023-06-29 | 2023-10-24 | 张家港广大特材股份有限公司 | 一种用于特种钢材冶炼的加工管控方法及系统 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102279928B (zh) * | 2011-07-20 | 2013-04-03 | 北京航空航天大学 | 基于支持向量机和模糊信息粒化的产品性能退化区间预测方法 |
| CN103761423B (zh) * | 2013-12-31 | 2016-06-29 | 中南大学 | 一种基于pso-elm的热轧板材组织-性能预测方法 |
| CN104200268B (zh) * | 2014-09-03 | 2017-02-15 | 辽宁大学 | 一种基于粒子群优化极限学习机的带钢出口厚度预测方法 |
| CN105354646B (zh) * | 2015-12-04 | 2017-04-12 | 福州大学 | 一种混合粒子群优化与极限学习机的电力负荷预测方法 |
| CN106651012A (zh) * | 2016-12-02 | 2017-05-10 | 东华大学 | 一种基于改进极限学习机的作物蒸腾量预测方法 |
| CN107400935B (zh) * | 2017-08-23 | 2019-10-18 | 东华大学 | 基于改进elm的熔融纺丝工艺的调节方法 |
-
2018
- 2018-10-19 CN CN201811221084.6A patent/CN109472397B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109472397A (zh) | 2019-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109472397B (zh) | 基于粘度变化的聚合工艺参数调节方法 | |
| CN112116147A (zh) | 一种基于lstm深度学习的河流水温预测方法 | |
| CN106600059A (zh) | 基于改进rbf神经网络的智能电网短期负荷预测方法 | |
| Deng et al. | Incremental learning of dynamic fuzzy neural networks for accurate system modeling | |
| Luo et al. | A laguerre neural network-based ADP learning scheme with its application to tracking control in the Internet of Things | |
| Han et al. | Hierarchical extreme learning machine for feedforward neural network | |
| Liu et al. | Melt index prediction based on fuzzy neural networks and PSO algorithm with online correction strategy | |
| CN113191092A (zh) | 一种基于正交增量随机配置网络的工业过程产品质量软测量方法 | |
| Tomasiello et al. | On fractional Tikhonov regularization: application to the adaptive network-based fuzzy inference system for regression problems | |
| CN108594793A (zh) | 一种改进的rbf飞控系统故障诊断网络训练方法 | |
| CN116484747A (zh) | 一种基于自适应优化算法与深度学习的污水智能监控方法 | |
| Assaleh et al. | Modeling nonlinear behavior of Buckling-Restrained Braces via different artificial intelligence methods | |
| CN115689070B (zh) | 基于帝王蝶算法优化bp神经网络模型的能源预测方法 | |
| CN113947182A (zh) | 基于双阶段堆叠图卷积网络的交通流预测模型构建方法 | |
| CN110032069B (zh) | 一种基于误差补偿的聚酯纤维纺丝工序分段参数配置方法 | |
| CN110263380B (zh) | 一种纺丝工序级联建模分段区间参数配置方法 | |
| Dong et al. | Quality prediction of fused deposition molding parts based on improved deep belief network | |
| Chen et al. | Temperature distribution prediction in control cooling process with recurrent neural network for variable-velocity hot rolling strips | |
| Alessandri et al. | A recursive algorithm for nonlinear least-squares problems | |
| Zheng et al. | Dynamic rolling force prediction of reversible cold rolling mill based on BP neural network with improved PSO | |
| CN113365298B (zh) | 基于高阶模糊认知图的mr网络信号强度预测方法及系统 | |
| Zhu et al. | From numeric to granular models: A quest for error and performance analysis | |
| Liu et al. | Prediction of Biochemical Oxygen Demand Based on VIP-PSO-Elman Model in Wastewater Treatment | |
| CN111583990A (zh) | 一种结合稀疏回归和淘汰规则的基因调控网络推断方法 | |
| CN114970343B (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 |