基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法技术方案

本发明专利技术公开了基于迭代多输出‑马尔科夫链的空分系统氩馏分变量多步预测方法，包括以下步骤：S1：采集数据，构建数据集，将数据集进行归一化处理，并划分为训练集、验证集和测试集；S2：构建迭代多输出预测模型；S3：采用马尔科夫链进行误差校正；本发明专利技术将迭代多输出法与马尔科夫链结合，一方面通过参数寻优的方式确定预测模型的最佳输入输出维度，有效减少模型的迭代误差累积，增强模型的适用性与灵活性；另一方面，马尔科夫链能够通过计算预测值的误差状态确定其校正范围和方向，使预测值与真实值的偏离程度更小，预测精度更高，本发明专利技术对于空分系统的变量时间序列多步预测具有较好的适用性。

【技术实现步骤摘要】
基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法
本专利技术涉及自动化
，尤其涉及一种基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法。
技术介绍
空气分离是指通过膨胀、精馏等物理或化学反应，在低温条件下将空气分离为氮气、氧气和氩气的过程。空气分离过程是一种流程工业，生产过程连续进行且机理复杂，其中的一些关键指标变量数据是设备和工艺工况的重要体现。工业现场人员通常依据关键指标的变化趋势调节相关设备的运行，从而保障生产过程安全平稳运行。然而在实际的生产现场，一些重要指标难以实时检测或全面检测，使整体的工况调节和计划调度受到负面影响。如果能够通过数学建模方法对关键指标变量进行实时预测并预知其未来的变化趋势，将会为工业生产的正常运作与实时调度提供十分重要的信息，因此建立精确稳定的生产过程变量短期预测模型是必要且有意义的。传统的生产过程主要依据能量守恒定律、热力学定律及动力学方程等机理分析建立指标预测模型，但是空分生产过程涉及到多种物理化学反应，过程复杂、时变性强，预测变量的机理建模难度较大、可靠性低，并且建模方法的泛化能力一般。现代化工业生产现场的控制系统能够采集大量包含生产信息的数据，基于数据驱动的预测建模方法通过建立输入与输出变量的关系模型达到预测的目的，而无需对系统的机理信息过多关注，因此成为了实现生产过程指标变量预测的有效途径。变量的短期预测通常根据预测步长分为单步预测和多步预测。相较于单步预测，多步预测能够实现更长时间的指标趋势预测，提供更多的有效信息，但同时也伴随着更多的不确定性和随机性。例如，中国专利CN95116365.5公开了空气分离系统的模型预测控制方法。该申请对于空气分离采用中间产物预测，但是其预测结果依然存在较大的不确定性和随机性。
技术实现思路
本专利技术主要解决现有的技术中空气分离对于预测结果存在较大的不确定性和随机性的问题；提供基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，降低预测结果的随机性，提高预测精度，具有较好的适用性。本专利技术的上述技术问题主要是通过下述技术方案得以解决的：基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，包括以下步骤：S1：采集数据，构建数据集，将数据集进行归一化处理，并划分为训练集、验证集和测试集；S2：构建迭代多输出预测模型；S3：采用马尔科夫链进行误差校正。从空分工业现场的监控与数据采集系统中获取氩馏分指标变量一段连续的时间序列数据，将数据集进行归一化，并划分为训练集、验证集和测试集，将训练集用于预测模型的训练，将测试集用于验证预测模型的效果，验证集用于评估训练模型，经过多输出预测模型和误差校正后提高空气分离结果的预测精度。作为优选，步骤S2中构建迭代多输出预测模型具体包括以下步骤：S21：建立迭代多输出神经网络预测模型；S22：确定迭代多输出预测模型的最优参数。作为优选，步骤S21中建立迭代多输出神经网络预测模型的具体方法为：采用具有较强的非线性拟合能力的人工神经网络作为指标变量的迭代多输出预测模型，设定模型的输入维度为m，隐含层节点数为p，输出维度为s，训练模型表示为：[yt+s,…,yt+1]＝F(yt,…,yt-m+1)+w,t∈{m,…,N-s}式中，s>1时，F:Rr→Rs表示向量函数，w∈Rs表示噪声向量，设训练好的预测模型为则多步预测值的获取方法为：其中q∈{1,…,n}，ceil(X)表示X向上取整，[y1,…,yN]表示多个历史实际值的变量时间序列，[yN+1,…,yN+H]表示预测数值，N代表时间序列的历史数值的数量，H(H>1)代表预测总步长。预测过程中，预测模型在已有的真实输入数据基础上，进行一次输出步长为s的多步预测，模型的多个输出值将成为下次预测的输入变量。从第二次预测开始，模型的输入变量包含了预测数据，以此展开迭代直到完成H步的数据预测。作为优选，步骤S22中确定迭代多输出预测模型的最优参数的具体方法为：确定隐含层节点数：在不同输入维度m、隐含层节点数p、输出维度s的参数组合下，根据平均绝对误差(MeanAbsoluteError,MAE)、均方根误差(RootMeanSquareError,RMSE)和平均绝对百分比误差(MeanAbsolutePercentageError,MAPE)进行误差评价，得到最小误差，获得最优网络结构参数，其具体计算式为：其中，n是样本量，yi和分别代表真实值和预测值。采用多个评价指标方法进行误差评价，提高准确度。作为优选，步骤S3中采用马尔科夫链进行误差校正的具体方法包括以下步骤：S31：建立转态转移矩阵；S32：预测值校正。作为优选，步骤S31中将预测值与真实值的相对误差根据其分布密度划分为n个状态区间，记为S1,S2,…,Sn，S1,S2,…,Sn通常按照区间分割点的数值从小到大排列，表示预测值被划分到不同区间时，处于从相对低估到相对高估的不同状态，每个相对误差数据都有其对应的误差状态，统计相对误差状态转移的概率时，用表示状态Si经过m时步转移到状态Sj的概率：式中，Mij(m)表示间隔时步为m时Si转移到Sj的次数，Mi表示误差状态为Si的数据个数，i,j∈{1,2,…,n}，Mi的统计过程中末位时序的变化状态是未知的，在统计时需要将最后的m个状态移除，由此得出第m时步的状态转移矩阵P(m)：作为优选，步骤S32中从预测数据点的相邻实测点开始，连续取N个实际值，N的取值与状态区间的个数相同(N＝n)，根据实测点距离预测点的时步间隔，将转移步数设定为i＝1,2,…n，转移步数i对应的实测点的起始状态为Sk，k＝1,2,…n；每个转移步数都有各自的起始状态和状态转移矩阵，按转移步数从小到大递增的规则，依次取步数i对应转移矩阵P(i)中的第k行向量组成新的矩阵，并对该矩阵的列向量求和，其中和最大的列数即是预测值与所对应真实值的相对误差所处的状态区间数，相对误差状态区间能够描述预测值的校正范围和校正方向，取预测值对应状态区间的中点作为校正系数定义预测校正值F(x)的表达式为：式中，f(x)表示预测模型的预测值；ΔU和ΔD分别表示预测值对应的误差状态区间的上限和下限。本专利技术的有益效果是：本专利技术将迭代多输出法与马尔科夫链结合，一方面通过参数寻优的方式确定预测模型的最佳输入输出维度，有效减少模型的迭代误差累积，增强模型的适用性与灵活性；另一方面，马尔科夫链能够通过计算预测值的误差状态确定其校正范围和方向，使预测值与真实值的偏离程度更小，预测精度更高，本专利技术对于空分系统的变量时间序列多步预测具有较好的适用性。附图说明图1是本专利技术实施例的多步预测方法的流程图。图2是本专利技术实施例的输出维度s＝3时不同网络结构的均方根误差示意图。具体实施方式下面通过实本文档来自技高网...

【技术保护点】
1.基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，其特征在于，包括以下步骤：/nS1：采集数据，构建数据集，将数据集进行归一化处理，并划分为训练集、验证集和测试集；/nS2：构建迭代多输出预测模型；/nS3：采用马尔科夫链进行误差校正。/n

【技术特征摘要】
1.基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，其特征在于，包括以下步骤：
S1：采集数据，构建数据集，将数据集进行归一化处理，并划分为训练集、验证集和测试集；
S2：构建迭代多输出预测模型；
S3：采用马尔科夫链进行误差校正。


2.根据权利要求1所述的基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，其特征在于，
步骤S2中构建迭代多输出预测模型具体包括以下步骤：
S21：建立迭代多输出神经网络预测模型；
S22：确定迭代多输出预测模型的最优参数。


3.根据权利要求2所述的基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，其特征在于，
步骤S21中建立迭代多输出神经网络预测模型的具体方法为：采用具有较强的非线性拟合能力的人工神经网络作为指标变量的迭代多输出预测模型，设定模型的输入维度为m，隐含层节点数为p，输出维度为s，训练模型表示为：
[yt+s,…,yt+1]＝F(yt,…,yt-m+1)+w,t∈{m,…,N-s}
式中，s>1时，F:Rr→Rs表示向量函数，w∈Rs表示噪声向量，设训练好的预测模型为则多步预测值的获取方法为：



其中q∈{1,…,n}，ceil(X)表示X向上取整，[y1,…,yN]表示多个历史实际值的变量时间序列，[yN+1,…,yN+H]表示预测数值，N代表时间序列的历史数值的数量，H(H>1)代表预测总步长。


4.根据权利要求2所述的基于迭代多输出-马尔科夫链的空分系统氩馏分变量多步预测方法，其特征在于，
步骤S22中确定迭代多输出预测模型的最优参数的具体方法为：确定隐含层节点数：



在不同输入维度m、隐含层节点数p、输出维度s的参数组合下，根据平均绝对误差(MeanAbsoluteError,MAE)、均方根误差(RootMeanSquareError,RMSE)和平均绝对百分比误差(MeanAbsolutePercentageError,MAPE)进行误差评价，得到最小误差，获得最优网络结构参数，其具体计算式为：



<...

【专利技术属性】
技术研发人员：郑松，史佳霖，葛铭，郑小青，魏江，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33