一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法技术

本发明专利技术公开一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法，旨在解决如何避免信息丢失的情况下利用SAE逐层提取丁烷塔过程测量数据的非线性特征，并充分有效利用这些非线性特征对塔底部产品的丁烷浓度进行在线软测量问题。首先，本发明专利技术方法在训练模型时，各层自编码器的输入都使用了原过程数据，因此在特征提取的过程中避免了输入数据的信息丢失问题。其次，本发明专利技术方法将各层神经网络的隐层输出向量做为单层神经网络的输入，从而较全面的利用到了SAE提取的非线性特征。与传统方法相比，本发明专利技术方法在线监测脱丁烷塔底部产品丁烷浓度时，其实时软测量精度是能得到保障的。

【技术实现步骤摘要】
一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法
本专利技术涉及一种软测量技术，特别涉及一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法。
技术介绍
脱丁烷塔是炼油工业过程中非常重要的设备，其在炼油过程中起到的作用是分离石脑油和脱硫处理。脱丁烷塔的的工艺流程一般包含6个化工设备：塔顶部冷凝器，热交换器，塔顶部回流泵，塔底部再沸器，回流储液器，和进料泵。在石脑油的流经通道内，需要去除石脑油中的丙烷和丁烷。因此，对于丁烷塔质量控制要求而言，最小化塔底部产品的丁烷含量是非常重要的。为保证质量控制的实时性，在线实时测量丁烷的浓度是最直接，也是最有效的技术手段。然而，丁烷含量的测量跟化工过程测量温度、压力、流量是不同的，需要使用成分分析仪。此外，在线成分分析仪价格高昂，且后期维护成本较高，测量丁烷浓度一般采用离线分析仪器。但是，离线分析仪器无法保证对丁烷浓度的实时测量，不能满足丁烷浓度测量的实时性要求。因此，需要一种软测量技术来实时在线监测丁烷的浓度。近年来，栈式自编码器(StackedAuto-Encoder，缩写：SAE)在深度学习领域得到了广泛的关注，其优势在于逐层提取输入数据中的特征，已被广泛用于模式识别领域。由于脱丁烷塔运行过程中动态变化性强，运行环境中还存在随机干扰，因此对脱丁烷塔过程测量变量进行深度特征提取是非常有必要的。值得强调的是，SAE虽然可以分层递进的提取输入数据中的特征，但是在逐层提取特征的过程中会出现信息丢失问题。从这个角度上讲，无法做到对特征的充分有效挖掘。此外，SAE是一种无监督型的建模技术，而软测量模型归属于有监督模型。因此，如何有效的利用这些提取的特征用于实时软测量脱丁烷塔底部产品的丁烷浓度还是一个需要商榷的技术难题。
技术实现思路
本专利技术所要解决的主要技术问题是：如何避免信息丢失的情况下利用SAE逐层提取丁烷塔过程测量数据的非线性特征，并充分有效利用这些非线性特征对塔底部产品的丁烷浓度进行在线软测量。具体来讲，本专利技术方法通过将脱丁烷塔的过程测量数据皆做为各层自编码器的输入，在逐层提取非线性特征时，可以避免信息的丢失问题。此外，利用各层自编码器提取的非线性特征，与丁烷浓度的输出数据之间建立单层神经网络，从而实现对丁烷浓度的实时在线监测。本专利技术方法解决上述问题所采用的技术方案为：一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法，包括以下步骤：步骤(1)：利用丁烷塔设备中安装的7个测量仪表实时测量过程变量数据，直至采集得到N个样本数据，分别记录成7个列向量x1，x2，…，x7，并利用离线分析仪获取相同测量时刻塔底部产品丁烷含量的N个数据，对应记录成列向量y，其中xi由第i个过程变量的N个样本数据组成，i∈{1，2，…，7}分别对应于塔顶部温度，塔顶部压力，回流流量，底部产品出口流量，第6层塔板温度，塔底部温度A和塔底部温度B。由于涉及商业机密，塔底部温度测量点有两个，具体测量位置使用A和B来表示。步骤(2)：根据如下所示公式对x1，x2，…，x7和y分别实施归一化处理，对应得到7个数据向量和1个输出向量并将数据向量组建成一个数据矩阵其中，RN×5表示N×5维的实数矩阵，表示归一化处理后的第i个数据向量，xi(min)表示列向量xi中的最小值，xi(max)表示列向量xi中的最大值，y(min)与y(max)分别表示列向量y中的最小值和最大值，下标号i∈{1，2，…，7}。步骤(3)：根据如下所示公式组建输入矩阵Z中的第k列输入向量zk∈R28×1：zk＝[X(k+3)，X(k+2)，X(k+1)，X(k)]T②上式中，X(k+3)，X(k+2)，X(k+1)，X(k)分别表示数据矩阵X中的第k+3行，第k+2行，第k+1行和第k行的行向量，k∈{1，2，…，n}，R28×1表示28×1维的实数向量，n＝N-3，上标号T表示矩阵或向量的转置。步骤(4)：搭建一个由M层自编码器串联组成的SAE，并确定隐层神经元激活函数f(u)，输出层神经元激活函数ζ(u)，和各层自编码器的隐层神经元个数h1，h2，…，hM；其中u表示函数自变量。步骤(5)：利用反向传播(Back-Propagation，缩写为BP)算法依次训练SAE中第1层自编码器，第2层自编码器，直至第M层自编码器的隐层和输出层的权重系数W1，W2，…，WM和以及隐层和输出层的阈值b1，b2，…，bM和具体的实施过程如下所示。步骤(5.1)：第1层自编码器的输入层有28个神经元，隐层有h1个神经元，输出层有28个神经元，分别初始化隐层和输出层的权重系数和阈值。步骤(5.2)：以z1，z2，…，zn做为第1层自编码器的输入，同时以z1，z2，…，zn做为第1层自编码器的输出，利用BP算法训练得到第1层自编码器的隐层和输出层的权重系数和以及隐层和输出层的阈值和后，初始化m＝1，其中分别表示输出向量中的第4个，第5个，至第N个元素。步骤(5.3)：第m+1层自编码器的输入层有hm+28个神经元，隐层有hm+1个神经元，输出层有28个神经元，并分别初始化隐层和输出层的权重系数和阈值。步骤(5.4)：将第m层自编码器隐层的输出向量g1(m)，g2(m)，…，gn(m)组建成矩阵Gm＝[g1(m)，g2(m)，…，gn(m)]T，并将矩阵Gm与输入矩阵Z合并成一个矩阵Hm＝[Gm，ZT]T，其中，g1(m)，g2(m)，…，gn(m)的计算方式如下所示。步骤(5.5)：以矩阵Hm的n个列向量做为第m+1层自编码器的输入，同时以z1，z2，…，zn做为第m+1层自编码器的输出，再次利用BP算法训练得到第m+1层自编码器的隐层和输出层的权重系数和以及隐层和输出层的阈值和其中表示(hm+28)×hm+1维的实数矩阵，表示hm+1×1维的实数向量。步骤(5.6)：判断是否满足m+1＜M；若是，则设置m＝m+1后返回步骤(5.3)；若否，则训练结束，保留栈式自编码器所有的权重系数W1，W2，…，WM和以及阈值b1，b2，…，bM和步骤(6)：搭建一个只有输入层与输出层的单层神经网络，其中输入层包含H＝h1+h2+…+hM个神经元，输出层包含1个神经元，输出层的激活函数为φ(u)。步骤(7)：将矩阵G1，G2，…，GM合并成一个矩阵G＝[G1，G2，…，GM]T后，将矩阵G∈RH×n中的n个列向量做为单层神经网络的输入，将做为单层神经网络的输出，利用BP算法训练得到输出层的权重系数向量w0∈RH×1和阈值b0∈R1×1。上述实施步骤完成了对丁烷塔底部产品丁烷浓度的软测量建模，接下来就是利用在线实时测量的数据，实现对脱丁烷塔底部产品丁烷浓度的在线监测，具体实施步骤如下所示。步骤(8)：在最新采样时刻t，对7个采样数据v1(t)，v2(t)，…，v7(t)进行归一化处理，得到归一化后的数据具体的归一化方式如下所示：上式中，i∈{1，2，…，7}。步骤(9)：根据at本文档来自技高网...

【技术保护点】
1.一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法，其特征在于，包括以下步骤：首先，离线软测量建模阶段包括如下所示步骤(1)至步骤(7)；/n步骤(1)：利用丁烷塔设备中安装的7个测量仪表实时测量过程变量数据，直至采集得到N个样本数据，分别记录成7个列向量x

【技术特征摘要】
1.一种基于SAE的脱丁烷塔底部产品丁烷浓度实时监测方法，其特征在于，包括以下步骤：首先，离线软测量建模阶段包括如下所示步骤(1)至步骤(7)；
步骤(1)：利用丁烷塔设备中安装的7个测量仪表实时测量过程变量数据，直至采集得到N个样本数据，分别记录成7个列向量x1，x2，…，x7，并利用离线分析仪获取相同测量时刻塔底部产品丁烷含量的N个数据，对应记录成列向量y，其中xi由第i个过程变量的N个样本数据组成，i∈{1，2，…，7}分别对应于塔顶部温度，塔顶部压力，回流流量，底部产品出口流量，第6层塔板温度，塔底部温度A和塔底部温度B，其中A和B分别表示温度的测量位置；
步骤(2)：根据如下所示公式对x1，x2，…，x7和y分别实施归一化处理，对应得到7个数据向量和1个输出向量并将数据向量组建成一个数据矩阵



其中，RN×7表示N×7维的实数矩阵，表示归一化处理后的第i个数据向量，xi(min)表示列向量xi中的最小值，xi(max)表示列向量xi中的最大值，y(min)与y(max)分别表示列向量y中的最小值和最大值，下标号i∈{1，2，…，7}；
步骤(3)：根据如下所示公式组建输入矩阵Z中的第k列输入向量zk∈R28×1：
zk＝[X(k+3)，X(k+2)，X(k+1)，X(k)]T②
上式中，X(k+3)，X(k+2)，X(k+1)，X(k)分别表示数据矩阵X中的第k+3行，第k+2行，第k+1行和第k行的行向量，k∈{1，2，…，n}，R28×1表示28×1维的实数向量，n＝N-3，上标号T表示矩阵或向量的转置；
步骤(4)：搭建一个由M层自编码器串联组成的SAE，并确定隐层神经元激活函数f(u)，输出层神经元激活函数ζ(u)，和各层自编码器的隐层神经元个数h1，h2，…，hM；其中u表示函数自变量；
步骤(5)：利用BP算法依次训练SAE第1层自编码器，第2层自编码器，直至第M层自编码器的隐层和输出层的权重系数W1，W2，…，WM和以及隐层和输出层的阈值b1，b2，…，bM和具体的实施过程如下所示；
步骤(5.1)：第1层自编码器的输入层有28个神经元，隐层有h1个神经元，输出层有28个神经元，分别初始化隐层和输出层的权重系数和阈值；
步骤(5.2)：以z1，z2，…，zn做为第1层自编码器的输入，同时以z1，z2，…，zn做为第1层自编码器的输出，利用BP算法训练得到第1层自编码器的隐层和输出层的权重系数和以及隐层和输出层的阈值和后，初始化m＝1，其中分别表示输出向量中的第4个，第5个，至第N个元素；
步骤(5.3)：第m+1层自编码器的输入层有hm+28个神经元，隐层有hm+1个神经元，输出层有28个神经元，并分别初始化隐层和输出层的权重系数和阈值；
步骤(5.4)：将第m层自编码器隐层的输出向量g1(m)，g2(m)，…，gn(m)组建成一个矩阵Gm＝[g1(m)，g2(m)，…，gn(m)]T，并将矩阵Gm与输入矩阵Z合并成一个矩阵Hm＝[G...

【专利技术属性】
技术研发人员：葛英辉，朱莹，其他发明人请求不公开姓名，
申请(专利权)人：宁波大学，
类型：发明
国别省市：浙江;33