本发明公开了一种针对产品浓度曲线的内部热耦合空分塔控制装置，包括与内部热耦合空分塔直接连接的智能仪表、控制器和DCS系统。所述DCS系统包括上位机、控制站、存储装置、现场总线、数据接口、存储装置。智能仪表通过温度检测元件、压力检测元件、流量检测元件测量相关参数，上位机用以实现控制参数的求解，控制站根据得到的控制参数对控制器进行调整，控制器实现对内部热耦合空分塔的直接控制调整。本发明提供的控制装置能够很好地处理内部热耦合空分塔的强非线性特征，具有高效的在线运算速度，同时针对产品浓度曲线特性，减少中间塔板累积误差，并具有非常好的伺服跟踪控制效果、干扰抑制效果。
1.一种针对产品浓度曲线的内部热耦合空分塔控制装置，其特征在于，包括与内部热耦合空分塔直接连接的智能仪表、控制器和DCS系统；所述DCS系统包括上位机、控制站、存储装置、现场总线和数据接口，存储装置、控制站及上位机通过现场总线与数据接口相连接；所述智能仪表通过温度检测元件、压力检测元件、流量检测元件测量相关参数，并与数据接口连接；所述上位机用以实现控制参数的求解，包括浓度曲线描述模块、设定值转换模块、控制参数求解模块，并将求解的控制参数通过现场总线传递给控制站；所述控制站根据得到的控制参数通过与现场总线相连的数据接口，对控制器进行调整；所述控制器实现对内部热耦合空分塔的直接控制调整；所述上位机包括：
1)浓度曲线描述模块，通过智能仪表中的温度检测元件采集温度信息以及压力检测元件采集压强数据，采集到的数据通过数据接口进行电信号转换，并通过现场总线将检测信号输送到所述浓度曲线描述模块，用以推断出各塔板浓度分布曲线的静态描述函数如下：其中， 为t采样时刻第j块塔板的液相i组分预测浓度，Si,h(t)、Si,l(t)分别为t采样时刻内部热耦合空分塔高压塔、低压塔浓度曲线的表征位置，Xi,h_min表示高压塔i组分浓度曲线的最小浓度值，Xi,h_max表示高压塔i组分浓度曲线的最大浓度值，γi,h表示高压塔i组分浓度曲线表征位置处的斜率，Xi,l_min表示低压塔i组分浓度曲线的最小浓度值，Xi,l_max表示低压塔i组分浓度曲线的最大浓度值，γi,l表示低压塔i组分浓度曲线表征位置处的斜率，N为总塔板数；
2)设定值转换模块，根据浓度曲线描述模块所获得的浓度曲线参数，将浓度设定值转换为表征位置设定值，转换公式如下：
其中， 分别为塔顶的汽相氮组分浓度的设定值和塔底的液相氧组分浓度的设定值， 分别为高压塔和低压塔浓度曲线表征位置的设定值，ki,j为第j块塔板i组分的气液平衡系数，可以由Peng-Robinson状态方程计算而来，最终的计算公式如下：其中每层塔板上的气液两相的逸度系数 均可以由下式计算：
其中P是压强，T是温度，v是摩尔体积，R是气体常数，取8.3145，xi为混合物中i组分浓度， 为i1组分浓度， 为i2组分浓度，ai是i组分的引力参数， 是i1和i2两种组分间的引力参数，a是所有组分分子间引力参数的加权和，bi是i组分的范德华体积，b为所有组分范德华体积的加权和，A为由(9)式定义的系数，B为由(10)式定义的系数，Z为压缩因子；
3)控制参数求解模块，用以通过浓度曲线描述模块和设定值转换模块得到的当前时刻信息求取下一时刻的控制参数，采用如下代数方程组可求得控制参数：
yi,j(t)＝ki,jxi,j(t) (14)Qj(t)＝UovAΔTj(t) (15)其中，yi,j(t)为t采样时刻第j块塔板的气相i组分浓度，xi,j(t)为t采样时刻第j块塔板的液相i组分浓度，zi,j(t)为t采样时刻第j块塔板进料的i组分浓度，Qj(t)为t采样时刻第j块塔板的传热量，UovA为传热系数，ΔTj(t)为t采样时刻第j组塔板间温差，λ为汽化潜热，Lj(t)为t采样时刻第j块塔板的液相流量，Fj(t)为t采样时刻的进料流量，Vj(t)为t采样时刻第j块塔板的气相流量，Uj(t)为t采样时刻第j块塔板的液相采出流量，Gj(t)为t采样时刻第j块塔板的气相采出流量，qj(t)为t采样时刻的进料热状况；压强P的作用包含在气液平衡系数ki,j中，推导关系如2)中所述，Si,h、Si,l分别为内部热耦合空分塔高压塔和低压塔浓度曲线的表征位置， 分别为高压塔和低压塔浓度曲线表征位置的设定值，K11、K12、K21、K22为系统参数，可根据实际的控制需求调整，通常K11和K21在10-100之间取值，K12和K22在
100-1000之间取值，qF(t+1)、Ph(t+1)分别为第t+1采样时刻的进料热状况和高压塔压强，同时也是控制装置下一时刻的控制参数；
4)所述的上位机还用于设定系统参数K11、K12、K21、K22的值，以及设定高压塔塔顶的汽相氮组分浓度和低压塔塔底的液相氧组分浓度的设定值 显示当前时刻浓度测量值和控制参数求解模块求解出的下一时刻的控制参数，并将控制参数通过现场总线传递给控制站，控制站再通过数据接口对控制器进行调整，从而完成控制装置的控制作用；同时所述的上位机还将以上信息通过现场总线传递给存储装置，方便操作人员查阅历史记录，提高生产控制品质。
一种针对产品浓度曲线的内部热耦合空分塔控制装置
技术领域
[0001] 本发明涉及工业节能控制的非线性控制领域，特别地，涉及内部热耦合空分过程的控制装置设计。
背景技术
[0002] 空分装置是对空气进行分离，并得到氧、氮、氩等高纯工业气体的装置，它被广泛应用于石油、化工、冶金、电子、能源、航空航天、食品饮料、医疗保健等各种工业领域。所得到的氧、氮和氩产品在一个国家国民经济中的应用十分广泛。自从上世纪70年代的两次“石油危机”以来，能源危机加深，强烈地要求许多领域能源的有效利用。在能耗很大的空气分离工业中，能源成本占了空气产品价格的75％。于是出现这样的情况，一方面，由于现代工业的发展，一些大型工业项目如钢铁工业、化学工业、石油开采等都需要由大型空分装置提供空气制品，需求量也越来越大。另一方面，能耗成本也随着能源危机，变得越来越大。因此在这样的形势下，提高空气分离技术的能量效率显得刻不容缓。
[0003] 内部热耦合空分技术比常规空分技术节能40％以上，节能效果显著。然而，由于内部热耦合空分过程具有强耦合、强病态、强不对称性、强反向响应等复杂的非线性动态特性，该塔的控制策略设计显得尤为困难。传统的PID、内模控制方案等已经不能满足要求，在内部热耦合空分塔的过程控制当中，这些方案已经很难使空分过程稳定。而基于线性辨识模型的控制方案只能工作在稳态工作点附近，稍微增大干扰幅度，或者设定值阶跃变化，系统控制质量则出现明显下降。准确把握内部热耦合空分塔的非线性特性，并在此基础上实现内部热耦合空分塔高效节能过程有效的的非线性控制方案，是提高该过程的生产控制品质的前提，已经成为一项关键的空分节能技术，具有十分重要的意义。
发明内容
[0004] 为了克服目前内部热耦合空分塔的控制装置抑制干扰能力差、控制效果差、难以实现精确设定值跟踪的不足，本发明的目的在于提供一种抑制干扰能力良好、控制效果好、可以实现精确迅速的设定值跟踪的针对产品浓度曲线特性的内部热耦合空分塔的非线性控制装置。
[0005] 本发明解决其技术问题所采用的技术方案是：
[0006] 一种针对产品浓度曲线的内部热耦合空分塔控制装置，包括与内部热耦合空分塔直接连接的智能仪表、控制器和DCS系统。所述DCS系统包括上位机、控制站、存储装置、现场总线和数据接口，存储装置、控制站及上位机通过现场总线与数据接口相连接。所述智能仪表通过温度检测元件、压力检测元件、流量检测元件测量相关参数，并与数据接口连接。所述上位机用以实现控制参数的求解，包括浓度曲线描述模块、设定值转换模块、控制参数求解模块，并将求解的控制参数通过现场总线传递给控制站。所述控制站根据得到的控制参数通过与现场总线相连的数据接口，对控制器进行调整。所述控制器实现对内部热耦合空分塔的直接控制调整。所述上位机包括：
[0007] 1)浓度曲线描述模块，通过智能仪表中的温度检测元件采集温度信息以及压力检测元件采集压强数据，采集到的数据通过数据接口进行电信号转换，并通过现场总线将检测信号输送到所述模块，用以推断出各塔板浓度分布曲线的静态描述函数如下：
[0008]
[0009]
[0010] 其中， 为t采样时刻第j块塔板的液相i组分(氧气、氮气或者氩气)预测浓度，Si,h(t)、Si,l(t) 分别为t采样时刻内部热耦合空分塔高压塔、低压塔浓度曲线的表征位置，Xi,h_min表示高压塔i组分浓度曲线的最小浓度值，Xi,h_max表示高压塔i组分浓度曲线的最大浓度值，γi,h表示高压塔i组分浓度曲线表征位置处的斜率，Xi,l_min表示低压塔i组分浓度曲线的最小浓度值， Xi,l_max表示低压塔i组分浓度曲线的最大浓度值，γi,l表示低压塔i组分浓度曲线表征位置处的斜率，N为总塔板数。
[0011] 2)设定值转换模块，根据浓度曲线描述模块所获得的浓度曲线参数，将浓度设定值转换为表征位置设定值，转换公式如下：
[0012]
[0013]
[0014] 其中， 分别为塔顶的汽相氮组分浓度的设定值和塔底的液相氧组分浓度的设定值， 分别为高压塔和低压塔浓度曲线表征位置的设定值，ki,j为第j块塔板i组分的气液平衡系数，可以由Peng-Robinson状态方程计算而来，最终的计算公式如下：
[0015]
[0016] 其中每层塔板上的气液两相的逸度系数 均可以由下式计算：
[0017]
[0018] 混合物a和b的混合规则为：
[0019]
[0020]
[0021]
[0022]
[0023]
[0024] 其中P是压强，T是温度，v是摩尔体积，R是气体常数，取8.3145，xi为混合物中i组分(氧气、氮气或者氩气)的浓度， 为i1组分的浓度， 为i2组分浓度，ai是i组分的引力参数， 是i1和i2两种组分间的引力参数，a是所有组分分子间引力参数的加权和，bi是i组分的范德华体积，b为所有组分范德华体积的加权和，A为由(9)式定义的系数，B为由(10)式定义的系数，Z为压缩因子。
[0025] 3)控制参数求解模块，用以通过浓度曲线描述模块和设定值转换模块得到的当前时刻信息求取下一时刻的控制参数，采用如下代数方程组可求得控制参数：
[0026]
[0027]
[0028] yi,j(t)＝ki,jxi,j(t) (14)
[0029] Qj(t)＝UovAΔTj(t) (15)
[0030]
[0031]
[0032] 其中，yi,j(t)为t采样时刻第j块塔板的气相i组分浓度，xi,j(t)为t采样时刻第j块塔板的液相i组分浓度，zi,j(t)为t采样时刻第j块塔板进料的i组分浓度，Qj(t)为t采样时刻第j 块塔板的传热量，UovA为传热系数，ΔTj(t)为t采样时刻第j组塔板间温差，λ为汽化潜热， Lj(t)为t采样时刻第j块塔板的液相流量，Fj(t)为t采样时刻的进料流量，Vj(t)为t采样时刻第j块塔板的气相流量，Uj(t)为t采样时刻第j块塔板的液相采出流量，Gj(t)为t采样时刻第j块塔板的气相采出流量，qj(t)为t采样时刻的进料热状况；压强P的作用包含在气液平衡系数ki,j中，推导关系如2)中所述，Si,h、Si,l分别为内部热耦合空分塔高压塔和低压塔浓度曲线的表征位置， 分别为高压塔和低压塔浓度曲线表征位置的设定值，K11、K12、 K21、K22为系统参数，可根据实际的控制需求调整，通常K11和K21在10-100之间取值，K12和K22在100-1000之间取值，qF(t+1)、Ph(t+1)分别为第t+1采样时刻的进料热状况和高压塔压强，同时也是控制装置下一时刻的控制参数。
[0033] 作为优选的一种方案：所述的上位机还用于设定系统参数K11、K12、K21、K22的值，以及设定高压塔塔顶的汽相氮组分浓度和低压塔塔底的液相氧组分浓度的设定值显示当前时刻浓度测量值和控制参数求解模块求解出的下一时刻的控制参
数，并将控制参数通过现场总线传递给控制站，控制站再通过数据接口对控制器进行调整，从而完成控制装置的控制作用。同时所述的上位机还将以上信息通过现场总线传递给存储装置，方便操作人员查阅历史记录，提高生产控制品质。
[0034] 本发明的技术构思为：对内部热耦合空分过程中的浓度曲线特性进行精确描述，成功准确地把握内部热耦合空分塔的非线性动态特性，克服已有的控制装置抑制干扰能力差、控制效果差、难以实现精确设定值跟踪的不足，针对产品浓度曲线特性，通过对浓度曲线静态和动态特征的准确把握，减小中间塔板累积误差，从而设计出内部热耦合空分过程的抑制干扰能力良好、控制效果好、可以实现精确迅速的设定值跟踪的非线性控制装置。
[0035] 本发明的有益效果主要表现在：1.非线性控制方案建立在高精度非线性模型基础上，能够及时抑制干扰作用；2.控制方案较好地处理了耦合问题，能够快速准确地跟踪设定值变化。
附图说明
[0036] 图1是本发明所提出的内部热耦合空分塔的非线性控制装置的结构图。
[0037] 图2是上位机实现方法的原理图。
具体实施方式
[0038] 下面根据附图具体说明本发明。
[0039] 参照图1和图2，一种针对产品浓度曲线的内部热耦合空分塔控制装置，包括与内部热耦合空分塔1直接连接的智能仪表2、控制器8和DCS系统。所述DCS系统包括上位机6、控制站5、存储装置4、现场总线7和数据接口3，存储装置4、控制站5及上位机6通过现场总线7与数据接口3相连接。所述智能仪表2通过温度检测元件、压力检测元件、流量检测元件测量相关参数，并与数据接口3连接。所述上位机6用以实现控制参数的求解，包括浓度曲线描述模块9、设定值转换模块10、控制参数求解模块11，并将求解的控制参数通过现场总线7传递给控制站5。所述控制站5根据得到的控制参数通过与现场总线7相连的数据接口3，对控制器8进行调整。所述控制器8实现对内部热耦合空分塔1的直接控制调整。所述上位机6包括：
[0040] 1)浓度曲线描述模块9，通过智能仪表2中的温度检测元件采集温度信息以及压力检测元件采集压强数据，采集到的数据通过数据接口3进行电信号转换，并通过现场总线7将检测信号输送到所述浓度曲线描述模块9，用以推断出各塔板浓度分布曲线的静态描述函数如下：
[0041]
[0042]
[0043] 其中， 为t采样时刻第j块塔板的液相i组分(氧气、氮气或者氩气)预测浓度，Si,h(t)、Si,l(t) 分别为t采样时刻内部热耦合空分塔高压塔、低压塔浓度曲线的表征位置，Xi,h_min表示高压塔i组分浓度曲线的最小浓度值，Xi,h_max表示高压塔i组分浓度曲线的最大浓度值，γi,h表示高压塔i组分浓度曲线表征位置处的斜率，Xi,l_min表示低压塔i组分浓度曲线的最小浓度值， Xi,l_max表示低压塔i组分浓度曲线的最大浓度值，γi,l表示低压塔i组分浓度曲线表征位置处的斜率，N为总塔板数。
[0044] 2)设定值转换模块10，根据浓度曲线描述模块9所获得的浓度曲线参数，将浓度设定值转换为表征位置设定值，转换公式如下：
[0045]
[0046]
[0047] 其中， 分别为塔顶的汽相氮组分浓度的设定值和塔底的液相氧组分浓度的设定值， 分别为高压塔和低压塔浓度曲线表征位置的设定值，ki,j为第j块塔板i组分的气液平衡系数，可以由Peng-Robinson状态方程计算而来，最终的计算公式如下：
[0048]
[0049] 其中每层塔板上的气液两相的逸度系数 均可以由下式计算：
[0050]
[0051] 混合物a和b的混合规则为：
[0052]
[0053]
[0054]
[0055]
[0056]
[0057] 其中P是压强，T是温度，v是摩尔体积，R是气体常数，取8.3145，xi为混合物中i组分(氧气、氮气或者氩气)的浓度， 为i1组分的浓度， 为i2组分浓度，ai是i组分的引力参数， 是i1和i2两种组分间的引力参数，a是所有组分分子间引力参数的加权和，bi是i组分的范德华体积，b为所有组分范德华体积的加权和，A为由(9)式定义的系数，B为由(10)式定义的系数，Z为压缩因子。
[0058] 3)控制参数求解模块11，用以通过浓度曲线描述模块9和设定值转换模块10得到的当前时刻信息求取下一时刻的控制参数，采用如下代数方程组可求得控制参数：
[0059]
[0060]
[0061] yi,j(t)＝ki,jxi,j(t) (14)
[0062] Qj(t)＝UovAΔTj(t) (15)
[0063]
[0064]
[0065] 其中，yi,j(t)为t采样时刻第j块塔板的气相i组分浓度，xi,j(t)为t采样时刻第j块塔板的液相i组分浓度，zi,j(t)为t采样时刻第j块塔板进料的i组分浓度，Qj(t)为t采样时刻第j 块塔板的传热量，UovA为传热系数，ΔTj(t)为t采样时刻第j组塔板间温差，λ为汽化潜热， Lj(t)为t采样时刻第j块塔板的液相流量，Fj(t)为t采样时刻的进料流量，Vj(t)为t采样时刻第j块塔板的气相流量，Uj(t)为t采样时刻第j块塔板的液相采出流量，Gj(t)为t采样时刻第j块塔板的气相采出流量，qj(t)为t采样时刻的进料热状况；压强P的作用包含在气液平衡系数ki，j中，推导关系如2)中所述，Si，h、Si，l分别为内部热耦合空分塔高压塔和低压塔浓度曲线的表征位置， 分别为高压塔和低压塔浓度曲线表征位置的设定值，K11、K12、 K21、K22为系统参数，可根据实际的控制需求调整，通常K11和K21在10-100之间取值，K12和K22在100-1000之间取值，qF(t+1)、Ph(t+1)分别为第t+1采样时刻的进料热状况和高压塔压强，同时也是控制装置下一时刻的控制参数。
[0066] 所述的上位机6还用于设定系统参数K11、K12、K21、K22的值，以及设定高压塔塔顶的汽相氮组分浓度和低压塔塔底的液相氧组分浓度的设定值 显示当前时刻浓度测量值和控制参数求解模块11求解出的下一时刻的控制参数，并将控制参数通过现场总线
7传递给控制站5，控制站5再通过数据接口3对控制器8进行调整，从而完成控制装置的控制作用。同时所述的上位机还将以上信息通过现场总线7传递给存储装置4，方便操作人员查阅历史记录，提高生产控制品质。
[0067] 上述实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。
