[发明专利]一种内部热耦合精馏塔的自适应非线性控制系统及方法

权利要求书

1.一种内部热耦合精馏塔的自适应非线性控制系统，包括与内部热耦合精馏塔直接连接的现场智能仪表和DCS系统，所述DCS系统包括存储装置、控制站和上位机，所述现场智能仪表与存储装置、控制站和上位机相连，其特征在于：所述的上位机包括用以计算输出内部热耦合精馏塔控制变量值的自适应非线性控制器，所述自适应非线性控制器包括：

组分推断模块，用以从现场智能仪表获取温度，压强数据，计算内部热耦合精馏塔各块塔板的组分浓度，并将组分浓度计算结果存储在历史数据库当中，采用式(1)(2)得到：

Xi(k)=Pr(k)×α×10Ti(k)+cb-a-1α-1,i=1,2,......,f-1---(1)]]>

Xi(k)=Ps×α×10Ti(k)+cb-a-1α-1,i=f,f+1,......,n---(2)]]>

其中，k为当前采样时刻，下脚标i为塔板编号，1为塔顶编号，f为进料板编号，n为塔底编号，X_i(k)为k采样时刻液相轻组分浓度，P_r(k)为k采样时刻精馏段压强、P_s提馏段压强，T_i(k)为k采样时刻塔内各块塔板的温度，α为相对挥发度，a、b、c为安东尼常数；

模型参数自适应校正模块，用以实现模型参数的在线更新，在线拟合模型函数如下(3)、(4)：

X^i=Xmin,r+Xmax,r-Xmin,r1+e-kr(i-Sr),i=1,2,......,f-1---(3)]]>

X^i=Xmin,s+Xmax,s-Xmin,s1+e-ks(i-Ss),i=f,f+1,......,n---(4)]]>

其中，为第i块塔板处液相组分浓度预估浓度，X_min，r、X_max，r、k_r，X_min，s、X_max，s、k_s为拟合参数，S_r、S_s分别为内部热耦合精馏塔精馏段、提馏段液相组分浓度分布的位置；

自适应非线性控制律求解模块，用以根据当前组分浓度数据，参考轨迹和当前时刻操作变量值求取当前的控制变量的理想改变值，采用式(5)-(14)得到：

Yi(k)=αXi(k)(α-1)Xi(k)+1,i=1,2,......,n---(5)]]>

Qi(k)=UA×b(1a-ln{(Pr(k)+ΔPr(k)×t)/[Xi(k)+(1-Xi(k))/α]}]]>i＝1，2，……，f-1    (6)

-1a-ln{pS/[Xi+f-1(k)+(1-Xi+f-1(k))/α]})]]>

V₁(k)＝F(1-q(k)-Δq(k)×t)    (7)

L_n(k)＝F(q(k)+Δq(k)×t)      (8)

Lf-1(k)=Σi=1f-1Qi(k)λ---(9)]]>

V_f(k)＝V₁(k)+L_f-1(k)          (10)

Sr*=1+1krln(Xmax,r-Y1*/[α-(α-1)Y1*]Y1*/[α-(α-1)Y1*]-Xmin,r),i=1,2,......,f-1---(11)]]>

Ss*=n+1ksln(Xmax,s-Xn*Xn*-Xmin,s),i=f,f+1,......,n---(12)]]>

-V1(k)Y1(k)+Vf(k)Yf(k)-Lf-1(k)Xf-1(k)HΣj=1f-1kr(Xmin,r-Xj(k))(Xj(k)-Xmax,r)Xmax,r-Xmin,r---(13)]]>

=K1(Sr*-Sr(k))+K2Σi=1k(Sr*-Sr(i))t]]>

-Vf(k)Yf(k)+Lf-1(k)Xf-1(k)+FZf-Ln(k)Xn(k)HΣj=fnks(Xmin,s-Xj(k))(Xj(k)-Xmax,s)Xmax,s-Xmin,s---(14)]]>

=K3(Ss*-Ss(k))+K4Σi=1k(Ss*-Ss(i))t]]>

其中，k为当前采样时刻，t为采样周期，X_i(k)、Y_i(k)分别为k采样时刻第i块塔板轻组分液相轻组分浓度和汽相轻组分浓度，Q_i(k)为第i块塔板之间的热耦合量，UA为传热速率，X_i+f-1(k)为k采样时刻第i+f-1块塔板液相轻组分浓度，q(k)为k采样时刻进料热状况，P_r(k)为当k采样时刻精馏段压强，F为进料流率，Z_f为进料组分浓度，V₁(k)、V_f(k)分别为k采样时刻塔顶和进料板的汽相流率，L₁(k)、L_f-1(k)、L_n(k)分别为k采样时刻塔顶、第f-1块塔板和塔底的液相流率，H为持液量，λ为汽化潜热，Y₁^*、X_n^*分别为塔顶的汽相轻组分浓度Y₁、塔底的液相轻组分浓度X_n的设定值，X_j(k)、X_n(k)、X_f-1(k)分别为k采样时刻第j块塔板、塔底、第f-1块塔板的液相轻组分浓度，Y₁(k)、Y_f(k)分别为k采样时刻塔顶和进料板的汽相轻组分浓度，K₁，K₂，K₃，K₄为控制律参数，S_r^*，S_s^*分别为精馏段提馏段拐点参考轨迹，S_r(k)、S_r(i)分别为k和i采样时刻内部热耦合精馏塔精馏段液相组分浓度分布的位置，S_s(k)、S_s(i)分别为k和i采样时刻内部热耦合精馏塔提馏段液相组分浓度分布的位置，Δq(k)，ΔP_r(k)，分别为当前时刻内部热耦合精馏塔控制变量即进料热状况和精馏段压强的当前理想改变值。