一种二元精馏塔集中式控制方法技术

本发明专利技术公开一种二元精馏塔集中式控制方法及系统。通过辨识各个被控变量与控制变量之间的通道传递函数，建立二元精馏塔的模型。对获得的二元精馏塔模型求取并分析其逆模型的频域特性以确定集中式PID控制器中的滤波器参数，再根据期望的动态性能和鲁棒性确定PID控制器的调节参数。利用计算出的集中式PID控制器调节控制信号，使得二元精馏塔塔顶与塔底轻组分与重组分的纯度达到要求。本方法思路简单，易于实现，可对存在耦合作用的二元精馏塔实施有效的控制。

【技术实现步骤摘要】
一种二元精馏塔集中式控制方法及系统
本专利技术属于石油化工精馏塔的自动控制领域，特别是二元精馏塔集中式控制方法及系统。
技术介绍
精馏是石化工业中应用最为广泛的传质过程，约90％的产品的提纯与回收是通过精馏完成。精馏分离的机理主要是利用混合物中不同组分的挥发度不同进行分离。混合物经过精馏塔后，从塔顶或塔底采出的产品满足一定的纯度要求。二元精馏塔主要用来分离含两种组分的混合物。二元精馏塔是一个多变量过程，内部机理复杂，动态响应迟缓且控制变量与被控变量之间的耦合作用较为严重，因此对控制方案提出了较高的要求。在精馏塔的控制方案中，较为简单的是通过变量配对的方法将多变量的精馏过程拆分成多个单回路进行控制，这种方法在单端组分控制时较为有效，但是其并未对变量之间的耦合作用进行有效的补偿，而且配对的准则往往是依照静态的相对增益矩阵(RGA)，因此当对塔进行两端组分控制时，就要对过程进行综合分析。解耦控制是一个消除变量间耦合的好办法，但现有的解耦控制中，静态解耦虽容易实现，但同样面临忽视过程动态信息的问题，而理想解耦或简单解耦等动态解耦方法其解耦器或解耦后过程又难于求取，这都对设计控制器带来了难度。现在工业中约90％的控制器仍为PID控制器，如何设计有效的PID控制器，使其良好地控制二元精馏塔是一个值得研究的问题。
技术实现思路
有鉴于此，本专利技术的主要目的是提供一种二元精馏塔集中式控制方法。包括以下步骤：S1、测量二元精馏塔的各个输入和各个输出的数值，并将测量得到的值变送为电信号；S2、将S1步骤得到的电信号存储起来，经过存储时间t后，将电信号传送至系统辨识模块；S3、采用最小二乘辨识算法辨识各个控制通道的传递函数，建立二元精馏塔的传递函数矩阵模型；S4、根据S3步骤得到的传递函数模型，计算得到归一化逆模型的频率特性，然后通过分析逆模型频率特性设计滤波器。S5、根据S3步骤得到的传递函数模型和S4步骤得到的滤波器，依照给定的鲁棒性和动态性能要求，计算出PID控制器参数，并确定最终的集中式PID控制器；S6、利用由S5步骤得到的集中式PID控制器调节输入。所述存储时间t为从加入用于辨识各通道的阶跃信号时刻起，到输出达到新的稳态值为止的时间段。S3的具体过程为：根据所述阶跃信号的幅值h、采样周期Ts以及采集对象的输出Y(k)，按照下式计算得到各通道的传递函数中的参数θ＝(PT·P)-1PTZ其中P和Z分别是由h、Ts、Y(k)构成的系数矩阵，θ为用于标识各支路的传递函数的参数。S4的具体步骤为：S401、求取过程的稳态增益逆矩阵K。S402、根据过程传递函数矩阵分别求出与输出1对应的所有通道中的最小时滞τ1和与输出2对应的所有通道中的最小时滞τ2.。S403、将与输出1对应的所有通道的时滞部分减小τ1，与输出2对应的所有通道的时滞部分减小τ2后，得到用于设计滤波器的模型G0。S404、选定最大频率ωmax,在频率范围[0,ωmax]上画出归一化的逆模型的每一个元素的频率特性曲线(Nyquist图)。最大频率ωmax要满足归一化逆模型的所有元素的频率特性曲线都在复平面的右半平面。归一化的逆模型元素的频率特性按下式求出其中，j为虚数单位，ω为频率，K为过程的稳态增益逆矩阵，为逆模型频率特性矩阵，为归一化逆模型频率特性矩阵，下标i和k分别代表矩阵的第i行和k列的元素。S405、对归一化逆模型的元素的频率特性曲线进行逐列分析以确定控制器参数γk的值，下标k表示参数γk对应归一化逆模型的第k列元素，参数γk的选择标准是选择一个较小值使得经过补偿后的归一化逆模型的第k列所有元素的频率特性曲线都在稳定区域，归一化逆模型的第k列元素按照下式补偿其中，j为虚数单位，ω为频率，γk为对应的控制器参数，为归一化逆模型频率特性矩阵，为补偿后的归一化逆模型频率特性矩阵，下标i和k分别代表矩阵的第i行和第k列的元素，S406、利用得到的补偿后的归一化逆模型，在其每个元素的频率特性曲线上取N个点，这些点所对应的频率分别为ω1，ω2，…ωN，其中，ω1＝0。再利用复曲线拟合的方法得到对应的滤波器的传递函数其中，Qik(s)为由补偿后的归一化逆模型的第i行第k列元素拟合获得的滤波器，Aik和Bik为滤波器对应的参数值，Kik为过程的稳态增益逆矩阵的第i行第k列元素，s为拉普拉斯算子，参数按照如下的复曲线拟合方法求取：其中，Qik(jωm)为第i行第k列元素，频率特性曲线上在频率取ωm时对应的复数，|Qik(jωm)|是求取该复数的模值，Re(Qik(jωm))是求取该复数的实部，Im(Qik(jωm))是求取该复数的虚部，N为选取的频率点的个数，a，b，c，d为中间变量。S5的具体过程为S501、给定鲁棒性指标Ms的值，控制器调节参数λ1和λ2均按照下式求取S502、计算子PID控制器参数其中，Kp为比例增益，KI为积分增益，Kd为微分增益，Tf为滤波器时间常数，这些参数的具体求取方法如下：可得到最终的带滤波器的集中式PID控制器为：其中，C为总的集中式PID控制器，C1和C2为上述求取的子PID控制器，Qik为之前求取的滤波器。第二方面，本专利技术提供一种二元精馏塔集中式控制系统，包括依次连接的测量变送器和数据存储与输出单元，数据存储与输出单元之后还顺序连接有系统辨识单元、参数分析单元和控制器单元，其中测量变送器，用于采集输入信号和输出信号；数据存储与输出单元，将所述测量变送器传送来的数据进行存储，并输出至系统辨识单元；系统辨识单元，采用辨识算法识别出各通道的传递函数，构成传递函数矩阵；参数分析单元，对过程的传递函数矩阵进行分析，设计出集中式滤波器，并利用定鲁棒性指标计算出PID控制器参数；控制器单元，利用得到的滤波器参数和PID控制器参数构成集中式PID控制器，从而输出控制信号。依照本方法设计的控制器，可以有效的降低耦合对于系统输出的影响，从而使得有效的提升塔顶和塔底采出的产品质量。且控制器结构简单，易于实现，控制器参数可以直接获取，避免了根据经验整定的过程。附图说明图1为本专利技术实施例二元精馏塔集中式控制方法的流程图。图2为本专利技术实施例二元精馏塔集中式控制系统的结构图。图3为本专利技术实施例归一化逆模型的每个元素在规定频段内的频率特性曲线。图4为本专利技术实施例归一化逆模型的每个元素的稳定性判断原理图。图5为本专利技术实施例补偿后的归一化逆模型的每个元素在规定频段内的频率特性曲线。具体实施方式下面结合附图和具体实施例对本专利技术进行详细说明。以组分为甲醇和水的二元精馏塔为例，该过程选取的输入(控制变量)为塔顶的回流量和塔底的蒸汽流量，输出(被控变量)选择为塔顶甲醇组分的摩尔分数和塔底甲醇的摩尔分数。过程的输入与输出之间存在较强的耦合，控制目标为塔顶与塔底的产品中的甲醇的摩尔分数均要达到设定的要求。二元精馏塔集中式控制方法实施例如图1所示，包括以下步骤：步骤1、从过程的输入端测量塔顶回流量值，塔底蒸汽流量值，从过程输出端测量塔顶与塔底的甲醇摩尔分数，将流量值和摩尔分数变送为电信号；步骤2、存储步骤1中的电信号，直到满足辨识各通道所需要的数量后，将电信号传送至系统辨识模块；步骤3、依照设定的最小二乘辨识算法识别各通道的传递函数，并把辨识出的传递函数组成传递函数矩阵；经本文档来自技高网...

【技术保护点】
一种二元精馏塔集中式控制方法，其特征在于，包括以下步骤：S1、测量二元精馏塔的各个输入和各个输出的数值，并将测量得到的值变送为电信号；S2、将S1步骤得到的电信号存储起来，经过存储时间t后，将电信号传送至系统辨识模块；S3、采用最小二乘辨识算法辨识各个控制通道的传递函数，建立二元精馏塔的传递函数矩阵模型；S4、根据S3步骤得到的传递函数模型，计算得到归一化逆模型的频率特性，然后通过分析逆模型频率特性设计滤波器，所述归一化逆模型的频率特性按下式求出

【技术特征摘要】
1.一种二元精馏塔集中式控制方法，其特征在于，包括以下步骤：S1、测量二元精馏塔的各个输入和各个输出的数值，并将测量得到的值变送为电信号；S2、将S1步骤得到的电信号存储起来，经过存储时间t后，将电信号传送至系统辨识模块；S3、采用最小二乘辨识算法辨识各个控制通道的传递函数，建立二元精馏塔的传递函数矩阵模型；S4、根据S3步骤得到的传递函数模型，计算得到归一化逆模型的频率特性，然后通过分析逆模型频率特性设计滤波器，所述归一化逆模型的频率特性按下式求出其中，j为虚数单位，ω为频率，K为过程的稳态增益逆矩阵，为逆模型频率特性矩阵，为归一化逆模型频率特性矩阵，下标i和k分别代表矩阵的第i行和第k列的元素。S5、根据S3步骤得到的传递函数模型和S4步骤得到的滤波器，依照给定的鲁棒性和动态性能要求，计算出PID控制器参数，并确定最终的集中式PID控制器；S6、利用由S5步骤得到的集中式PID控制器调节输入。2.根据权利要求1所述的二元精馏塔集中式控制方法，其特征在于，所述存储时间t为从加入用于辨识各通道的阶跃信号时刻起，到输出达到新的稳态值为止的时间段。3.根据权利要求1所述的二元精馏塔集中式控制方法，其特征在于，所述S3的具体步骤为：根据所述阶跃信号的幅值h、采样周期Ts以及采集对象的输出Y(k)，按照下式计算得到各通道的传递函数中的参数θ＝(PT·P)-1PTZ其中P和Z分别是由h、Ts、Y(k)构成的系数矩阵，θ为用于标识各支路的传递函数的参数。4.根据权利要求1所述的二元精馏塔集中式控制方法，其特征在于，S4的具体过程为：S401、求取过程的稳态增益逆矩阵K，S402、根据过程传递函数矩阵分别求出与输出1对应的所有通道中的最小时滞τ1和与输出2对应的所有通道中的最小时滞τ2.，S403、将与输出1对应的所有通道的时滞部分减小τ1，与输出2对应的所有通道的时滞部分减小τ2后，得到用于设计滤波器的模型G0，S404、选定最大频率ωmax,在频率范围[0,ωmax]上画出归一化的逆模型的每一个元素的频率特性曲线，即Nyquist图，最大频率ωmax要满足归一化逆模型的所有元素的频率特性曲线都在复平面的右半平面，归一化的逆模型元素的频率特性按下式求出其中，j为虚数单位，ω为频率，K为过程的稳态增益逆矩阵，为逆模型频率特性矩阵，为归一化逆模型频率特性矩阵，下标i和k分别代表矩阵的第i行和第k列的元素，S405、对归一化逆模型的元素的频率特性曲线进行逐列分析以确定控制器参数γk的值，下标k表示参数γk对应归一化逆模型的第k列元素，参数γk的选择标准是选择一个较小值，使得经过补偿后的归一化逆模型的第k列所有元素的频率特性曲线都在稳定区域，归一化逆模型的第k列元素按照下式补偿其中，...

【专利技术属性】
技术研发人员：靳其兵，杜星瀚，蒋北艳，周星，
申请(专利权)人：北京化工大学，
类型：发明
国别省市：北京,11