一种基于变量时段分解的间歇过程无模型摄动优化方法
【技术领域】
[0001]本发明属于化工流程制造业领域，涉及一种针对间歇过程的操作轨线变量时段分解的无模型摄动优化方法，适用于包括间歇反应器、间歇精馏塔、间歇干燥、间歇发酵，间歇结晶以及其它采用间歇方式操作的过程和系统的操作轨线优化。
【背景技术】
[0002]随着激烈的市场竞争和对多品种、多规格、高附加值产品的需求，间歇生产方式逐渐成为企业生存的关键。一方面，尽管人们在过程控制领域的研究已经取得了众多成果，但大多集中在连续工业过程领域中，间歇过程的自动化水平仍然处于比较低级的阶段;另一方面，间歇过程的特殊操作方式也使得大部分现有方法难以直接应用。因此，在间歇生产过程中推行有效的优化与控制方法，以获得更满意的产品质量已成为间歇生产过程的普遍需求。寻找最优操作状态或轨迹，并控制过程在最优操作条件下运行，是产品质量控制研究领域的重要研究和亟待解决的问题。
[0003]然而，间歇生产过程的开/停车、产品和原料的切换、进料的变化以及批量生产过程，是一个一直在“进行中”的动态过程。另外，间歇过程的优化面对的一般是非线性、多目标、机理复杂的优化问题。由于间歇过程的机理模型很难得到，因而基于数据驱动的优化方法凸显出其必要性和重要性。尽管在间歇过程数据驱动的优化方面，已经取得了一系列研究成果，但目前仍缺乏适合工业应用的有效方法。因此提出一个能解决实际工业问题的间歇过程优化策略和框架，可以提高过程效率，降低能耗，提高经济效益，直接地促进间歇过程优化技术的工业化应用。
【发明内容】
[0004]本发明提出的方法，将待优化轨线在时间轴上进行划分，通过方差分析建立摄动优化策略。完全基于生产过程的操作数据，不需要过程机理的先验知识和机理模型，对一般间歇过程优化问题具有通用性和有效性。
[0005]本发明为实现上述目的，采用如下技术方案:
[0006]本发明数据收集和计算步骤如下:
[0007]步骤一:针对操作完整的间歇过程，按批次收集待优化变量和最终质量或产率指标。数据的收集时间间隔可以是等时间间隔或非等时间间隔，要保证在一个时间间隔内，过程的待优化变量没有显著变化，或其变化不会对最终质量或产率指标有显著影响。一般要求30-50组有效数据。
[0008]步骤二:对所采集的数据，按批次操作轨线为变量进行主元分析并在主元模式图中剔除奇异点，使得所有数据点在一个可信度之内。
[0009]步骤三:将剔除奇异点后的剩余批次操作曲线在时间轴上进行等间隔划分或不等间隔划分。
[0010]步骤四:将每一个间隔所包含的各批次数据表达为一个连续变量，这些变量被称为分解后的时段变量。时段变量的值由待优化变量在一个特定时间区间的各个批次数据所组成。
[0011]步骤五:将步骤四中所对应的每个批次质量或收率指标，称为指标变量。指标变量的值是由各个批次最终质量或收率形成的连续变量。
[0012]步骤六:将步骤四和步骤五中形成的时段变量和指标变量进行合并，形成时段变量和指标变量的联合变量数据矩阵。
[0013]步骤七:对上述联合变量矩阵做主元分析，形成主元载荷图。
[0014]步骤八:对步骤七中的主元载荷图按时段变量对指标变量的作用方向和大小进行分类。一般可以分为正作用，反作用及无(微)作用三类。
[0015]步骤九:按上述划分的不同类别，分别计算出各个时段变量的均值和标准差。按以下摄动法建立各时段变量的优化策略:
[0016]J(i)=M(i)+sign(i) Χ3σ(?)
[0017]此处的1(0，1(0和0(0分别是第1个时段变量的优化目标值，均值和标准差；sign(i)是第i个时段变量和指标变量所形成的夹角余弦符号。如图5分类图上，夹角小于90度时为+1，夹角大于90度时为-1，夹角为90度时为0。
[0018]步骤十:将步骤九中所得各时段的优化目标值，按时段顺序i= l，2，…N排列，组合成一条针对整个批次过程的新优化变量曲线。
[0019]步骤十一:一般地需要将上述优化曲线进行数字滤波，使得新的优化曲线比较光滑，易于跟踪控制。
[0020]本发明方法首先采集实际运行数据，在全局过程主元分析的基础上，依照各个时间段内操作条件对最终质量指标的灵敏度，对整个过程做时间轴上的元素变量划分并构造成一系列数据子集。分别对这些数据子集做方差分析，形成在各个子集上的摄动优化作用量，最后综合成为全局优化策略。本方法完全基于间歇过程的正常操作数据形成优化策略，不需要过程本身机理的先验知识和模型，具有对一般间歇过程优化问题的通用性和有效性，对实际生产过程具有指导意义。
【附图说明】
[0021]图1一个间歇过程的温度操作曲线示例
[0022]图2—个间歇过程温度为优化变量的主元模式图
[0023]图3时段变量的构成
[0024]图4时段变量和指标变量主元载荷图
[0025]图5时段变量对指标变量作用分类图
[0026]图6—个结晶过程的优化温度曲线和原始温度曲线比较
[0027]图7实施步骤框图
[0028]图8滑动平均滤波处理后的优化曲线与滤波前优化曲线的比较
[0029]图9指标变量优化策略生成计算
[0030]图10—个间歇结晶过程的优化结果
【具体实施方式】
[0031]以下以一个间歇结晶过程为例，说明【具体实施方式】。这个示例不构成对本发明方法的范围限制。
[0032]本发明分为三个部分。第一部分为数据收集和预处理。第二部分是联合变量数据矩阵的构造。第三部分是摄动优化策略的计算。
[0033]本发明实施步骤框图如图7所示，具体实施步骤和算法如下:
[0034]步骤1:针对操作完整的间歇结晶过程，温度是与产品收率密切相关的待优化变量，按批次收集50组温度变量以及最终收率指标数据。数据的收集时间间隔为1分钟。图1是一个间歇结晶过程的温度曲线数据收集示例，为清晰起见，图中仅画出2个批次的温度曲线。
[0035]步骤2:对所采集的50批次的温度数据，按批次操作温度为变量进行主元分析并在主