一种高炉系统与 TRT 装置分布式协调控制方法 【技术领域】
本发明涉及能源与动力工程技术,特别是涉及一种高炉系统和 TRT 装置协调优 化控制的方法。背景技术
钢铁工业是关系国民经济命脉的重要产业,但是高炉冶炼能耗高和污染严重的 特点影响钢铁工业的持续发展。 高炉冶炼是炼铁过程中最耗能的部分,而且高炉冶炼中 经常会出现干扰,产生波动,影响高炉的稳定。 高炉冶炼主要是把烧结矿、焦炭放到高 炉中,通过进风口向高炉中鼓吹热风,反应中产生的还原性气体将矿石中的铁元素还原 出来。 矿石的成分变化,烧结矿转鼓强度的变化,以及热风压力变化都会对高炉炉内的 温度分布产生较大影响,主要通过进风操作来调节炉内温度。 高炉煤气余压透平发电装 置 ( 以下简称 TRT 装置 ) 是一种工业二次能源回收装置,利用原来浪费在减压阀组上的 高炉煤气余压进行发电,产生电能,消除减压阀组上的噪声,还能调节高炉顶压稳定, 保障高炉冶炼的正常运行。 TRT 装置对钢铁企业进行节能减排技术具有重要意义。
TRT 装置与高炉系统紧密相连,相互影响。 高炉系统对 TRT 装置的影响很大, 高炉产生的煤气量越大,压力越大,煤气在透平机中做功越多,发电量越大。 如果高炉 系统产生的煤气量波动较大,就需要对透平机进行较多的调节,包括对快开旁通阀的调 整,影响 TRT 装置的正常运行。
反过来, TRT 装置对高炉系统也有很大的关联影响作用。 尤其是在 TRT 装置 启动、升速和紧急停止过程中高炉顶压会发生不稳定,主要通过 TRT 透平机静叶开度来 调节煤气流量,达到调节高炉顶压的目的。 顶压是高炉冶炼的重要参数,如果顶压不能 稳定,会引起炉内反应剧烈波动,温度分布发生变化,高炉出现异常。 当炉内压力高于 额定值时,会使煤气气流分布不均,引起崩料,严重时会损坏设备。 当炉内压力低于额 定值时,会引起炉内煤气体积增大,气流压力损失增大,煤气流速上升,使炉喉磨损严 重,更有甚者,会引起炉顶设备事故。 就 TRT 装置的发电量与高炉生产相比较而言, TRT 装置处于从属地位,当高炉炉顶压力波动而影响其正常生产时, TRT 装置会被迫停 机,保证高炉冶炼的正常运行。
高炉系统和 TRT 装置紧密相连,相互影响。 实际炼铁厂中采取的优化控制技 术都是单独针对高炉系统或者 TRT 装置本身来研究的,各自为政,缺乏协调,将来自于 另一系统的输入看作已知量,相当于切断了相互之间的关联耦合作用,孤立的对本地系 统进行优化控制,没有考虑与相关工序的协调优化控制,得到的控制效果只是局部最优 ( 即对本工序而言是最优的 ),而非整体最优 ( 对整个高炉系统和 TRT 装置而言最优 ), 就会导致控制作用过量,影响正常生产,造成资源浪费和效率不高。 发明内容
本发明的目的在于提供一种高炉系统与 TRT 装置分布式协调控制方法,即从全局的、整体的角度对高炉系统和 TRT 装置进行协调优化控制,通过数据服务器进行信息 交换,高炉系统和 TRT 装置分别在本地进行优化控制,能够有效缓解原料成分波动和其 他扰动作用对生产操作和产品质量的影响,减少不必要能源浪费,保障高炉系统和 TRT 装置平稳、高效的运行。
本发明采用的技术方案的步骤如下 :
1 高炉系统的描述和动态建模 :
首先对高炉系统进行描述,再根据关联模型总体架构,在对高炉系统描述和分 析的基础上,建立高炉系统的动态关联模型 ;
2TRT 装置的描述和动态建模 :
首先对 TRT 装置进行描述,再根据关联模型总体架构,在对 TRT 装置的描述与 分析基础上,建立 TRT 装置的动态关联模型 ;
3 高炉系统和 TRT 装置分布式协调控制 :
采用分布式协调控制方法,首先确定高炉系统和 TRT 装置的采样周期,然后 对两个模型进行优化求解,采用可变窗口长度的分布式预测控制方法,建立协调控制机 制,实现高炉系统和 TRT 装置协调控制。
所述的对高炉系统进行描述,即考虑影响高炉温度分布和正常运行的主要因 素,按照协调优化的关联模型结构,以高炉顶压作为关联输入,以热风流量、热风温 度、富氧量和喷煤量为控制量,以转鼓强度、热风压力、矿石成分和焦炭质量为干扰, 以高炉煤气流量和煤气温度为关联输出,以高炉炉内温度为内部状态,建立高炉系统的 动态关联模型。
所述的对 TRT 装置进行描述,通过机理分析,按照协调优化的关联模型结构, 以高炉煤气流量和温度为关联输入,以静叶开度和旁通阀开度为控制量,以雷达料线和 叶片积灰为干扰,以高炉顶压为关联输出,以炉顶压力、透平机进口压力、透平机进口 压力容器入口温度、透平机进口压力容器内部温度、透
平机进口压力容器内部压力、透平机出口压力容器温度和透平机出口压力容器 内部压力为内部状态,建立 TRT 装置的动态关联模型。
根据高炉系统 S1 和 TRT 装置 S2 的动态响应特性确定它们的采样周期为 T1 和 T2,其中 T1 > T2,利用前向差分法对高炉系统和 TRT 装置的连续模型进行近似离散 化 ;其中高炉系统 S1 和 TRT 装置 S2 的自身变量以各自的采样周期为准,关联变量的 采样方法如下 :高炉系统的关联输入高炉顶压为 :z1(t) = h1(x2(t)),离散化后变为 k = k0, k0+1, ..., kf, k0 和 kf 为自然数,且 k0 < kf, 表示以 T1 为采样周期对 TRT 装置 S2 的状态变量 x2(t) 进行虚拟采样得到的值,得到 S1 的离散模型 M′ 1,因为 T2 < T1,则不存在采样值缺失的问题 ;同理,TRT 装置的关联输入为煤气流 量和煤气温度 :z2(t) = h2(x1(t)),离散化后变为 kf, k0 和 kf 为自然数,且 k0 < kf, k = k0,k0+1,..., 表示以 T2 为采样周期对 S1 的状态变量 x1(t) 进行虚拟采样得到的值,得到 S2 的离散模型 M′ 2,由于 T2 < T1,则在采样不重合的时间 点,将出现采样值缺失的情况,即无法获取煤气流量和煤气温度的数据,缺失的采样值 用最近一次采样值代替。
采用二次型优化算法对两个动态关联模型进行优化求解 ;高炉系统以高炉模型和生产工艺参数允许范围为约束条件,以高炉内部温度分布与设定值之间的偏差最小以 及控制作用最小为优化目标 ;TRT 装置以 TRT 模型和生产工艺参数允许范围为约束条 件,以高炉顶压稳定、发电量最大和控制作用最小为优化目标,分别对高炉系统和 TRT 装置进行优化求解,解算出控制量发送给现场控制器。
可变窗口长度的分布式预测控制方法,是采用预测控制的滚动优化方法,在 t = kTi 时刻,对高炉系统 S1 和 TRT 装置 S2 的离散模型在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时 间段的系统行为进行预测和优化,其中 Pi 为 Si 在 t = kTi 时刻进行滚动优化的预测时域, Ls = Pi ·Ti 为 t = kTi 时刻滚动优化的时间窗口长度,其中 Si 为 S1 或 S2,Pi 为 P1 或 P2, Ti 为 T1 或 T2 ;当高炉系统与 TRT 装置协调优化时,预测时间窗口长度由高炉系统 S1 决 定。
可变窗口长度的分布式预测控制方法,是采用预测控制的模型预测和反馈校正 方法,在 t = kTi 时刻,利用 Si 的离散模型对 Si 在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时间段 的系统行为进行预测,利用实测值对预测结果进行反馈校正,减小偏差。 使用二次型优 化算法求解 Si 在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时间段的最优控制问题 Ji, k,其中 Si 为 S1 或 S2, P1 为 P1 或 P2, Ti 为 T1 或 T2 ;实现可变窗口长度的分布式预测控制。 所述的协调控制机制中协调优化过程是建在对系统未来行为预测的基础上,包 含一系列变长度预测窗口内的优化计算 ;每一个预测窗口内的优化计算需迭代多次,直 至收敛,协调优化机制包括 :
启动协调优化 :本地要求启动协调优化 ;关联部分要求启动协调优化 ;
停止协调优化 :每一个预测窗口内,高炉系统和 TRT 装置的控制量都保持稳 定,停止本窗口迭代优化,移动预测窗口 ;当高炉系统与 TRT 装置本地状态都恢复正 常,停止预测窗口的移动,整个协调优化过程结束。
本发明具有的有益效果是 :
本发明使高炉炉内温度分布偏差小,高炉炉顶压力稳定, TRT 装置发电量最 大,控制作用最小,实现两个系统整体性能最优,而且具有收敛速度快、计算效率高、 优化性能好等显著优点。
附图说明
图 1 是高炉系统与 TRT 装置关联模型总体结构图。
图 2 是 TRT 装置综合模型示意图。
图 3 是多输入单输出管接头模型。
图 4 是单输入多输出管接头模型。
图 5 是透平机模型示意图。 具体实施方式
下面将结合附图和具体实例对本实用新型作进一步的说明。
本发明具体步骤如下 :
一、高炉系统的描述与动态建模 :
高炉生产是在高温、高压、密闭条件下的物理、化学、动力学过程。 由于影响高炉安全稳定生产、提高铁水产量和质量的因素很多,炉内的气体温度、压力是相互影 响,都和炉况息息相关。 作为衡量高炉炉况和铁水质量的重要指标之一的温度信息,钢 铁厂主要从物理温度 ( 铁水温度 ) 和化学温度 ( 铁水硅含量 ) 中获得。 尽管高炉的炉壁 及炉缸位置装有为数不少的热电偶,其测得的温度可以在中控室计算机上显示,但因热 电偶无法测得炉体内部反应区的温度,平时通常只将其作为判断操作炉型的依据,而不 作为炉体温度参考。 将铁水温度、铁水硅含量作为温度指标的一个很大的缺点就是有滞 后性,如能通过某种方式实时获取炉内温度变化情况,并通过改变风温、喷煤量、风量 等手段调节及时进行调节,在原燃料波动较大的现状下,将对减少异常炉况的发生、提 高高炉冶炼质量带来很大帮助。
通过研究影响高炉温度分布和正常运行的主要因素,建立高炉的动态关联模 型。 按照图 1 所示的协调优化模型结构,以高炉顶压作为关联输入,以热风温度、富 氧量、喷煤量和热风流量为控制量,以转鼓强度、热风压力、矿石成分和焦炭质量为干 扰,以高炉煤气流量和煤气温度为关联输出,以高炉炉内温度为内部状态,建立高炉系 统的动态关联模型。
高炉偏差状态空间方程模型为 :
变量说明见下表 :二、 TRT 装置的描述和动态建模 :
1、建模基本假设
TRT 装置主要包括快开旁通阀组,减压阀组和透平发电机等几个主要部件,且 相互之间通过管道并联,综合考虑 TRT 装置的气体流程,可对其进行如下的简化 :
(1)TRT 装置中各个部件之间的距离较近,因此可以忽略相连管路的动态特性。
(2) 管道接头处近似为定容容器,进出口的气体质量流量基本一致。此处与外界 基本不发生能量交换,可以忽略进出气体之间的热容变化。
(3)TRT 装置中的旁通阀、减压阀和透平静叶的作用都是通过开度来调节流量, 这里将其近似为线性调节阀。
根据上述假设,在建模过程中,采用静态和动态相结合的方法将其分成两部 分。 静态部分主要是部件的特性方程,主要包括透平机性能曲线以及控制阀流量压差关 系等。 动态部分主要是由部件前后的管接头近似所得的压力容器建立。 TRT 装置的综合 模型如图 2 所示,其中,进口压力容器与高炉的炉顶压力容器相连,出口压力容器则直 接与煤气管网相连。
2、 TRT 装置机理模型 (1) 高炉压力容器模型
考虑 TRT 装置建模时,可以对高炉模型进行适当的简化。 这里将高炉简化为反 应容器和压力容器的串联模型。
根据压力容器入口和出口的质量流量的变化与容器中的密度变化的关系可以得 到:
式中, G2 为压力容器入口煤气的质量流量, Go_BF 为压力容器出口即高炉炉顶出 口的煤气质量流量, V 为压力容器的体积, ρ 为压力容器中气体的平均密度。
根据理想气体状态方程和忽略温度变化的影响,上式简化为 :
V = V0+S0H
式中,P2 为高炉顶压,M 为压力容器中气体平均分子量,R 为气体常数,T2,0 为 压力容器气体温度稳态值, H 为雷达料线高度,其中压力容器的体积受到料柱变化的影 响,将料柱达到料线时的体积视为常数 V0, S 为料柱对炉顶压力容器体积的影响系数。
(2) 除尘设备和管路简化模型
G1 为透平机进口煤气流量, T2, T3 为管道进出口温度, P3 为透平机进口压力, K1, K2, K3, K4 为常数, τ1, τ2, τ3 为时间常数,与实际工艺参数有关。(3)TRT 数学模型 多输入单输出型管接头模型,如图 3 所示 :Gout = Gin1+Gin2
式中, T 和 P 分别为压力容器的温度和压力, V 为压力容器体积, Gin1, Gin2, Gout 分别为进出口管道的气体质量流量, Pin1, Pin2 分别为进口管道的气体压力, Tin1, Tin2 分别为进口管道的气体温度,γ = Cp/Cv 为绝热指数,也叫比热比,Rg = R/M = 8.314/ M(J/Kg · K) 为气体常数。
单输入多输出型管接头模型,如图 4 所示 :
Gin1 = Gout1+Gout2
透平机模型
根据上述分析,可以认为透平机前后均与压力容器相连,如图 5 所示。 分析透 平膨胀机工作原理,可以得到透平机特性曲线方程 :
式中, Tout* 为透平机出口温度, Ψ 为透平机速度系数,一般取 0.96。
透平机的流量主要由透平机静叶开度来调节,其作用可以等效为控制阀,因此 根据控制阀模型可以得到透平机流量方程 :
式中, Gout* 为透平机出口流量, Kt 为透平机静叶开度的函数。 根据 TRT 煤气余压发电的透平出力计算公式,可以得到 :
式中, N 为透平机功率,也称透平出力, Gt, out* 为高炉煤气质量流量 (Kg/s) ; Cp- 高炉煤气定压比热 (kJ/m3) ;Tt, in 为透平入口煤气绝对温度 (K) ;Pt, in 为透平入 口煤气压力 (kPa,) ;Pt,out 为透平出口煤气压力 (kPa) ;γ 为绝热指数,高炉煤气 γ = 1.4 ;fd 为汽化潜热的热量修正系数 ( 干法除尘为 1.0) ;ηT 为透平机效率 ( 轴流反动式取 0.85) ;ηG 为发电机效率 (0.95 ~ 0.97)。
综上可得 TRT 综合模型 :
根据煤气流量物料守恒可得 :Gi_trt = Go_trt = Gt+Gp+Gj。
3、离散化状态空间方程模型
模型化简处理假设 :由于管接头入口压力和管接头内部压力相差不大,可近似 相等,即 Pi,trt = Pt,in,透平机进口压力容器和出口压力容器的体积近似相等,即 Vt,in = Vt, out。
化为离散化状态空间方程模型,首先选取高炉顶压和 TRT 功率为输出,高炉顶 压、透平机进口压力、透平机进口压力容器入口温度、透平机进口压力容器内部温度、 透平机进口压力容器内部压力、透平机出口压力容器温度、透平机出口压力容器内部压 力为状态,静叶开度和旁通阀开度为控制量,高炉压力容器进口煤气流量 ( 高炉煤气流 量 ) 和煤气温度为关联输入,雷达料线和透平机叶片积灰为本系统内部干扰,其中透平 机叶片积灰作为随机干扰,变量说明如下表 :
变量名 状态 x1 状态 x2 状态 x3物理意义 炉顶压力 透平机进口压力符号 P2 P3单位 Kpa Kpa ℃稳态值 160 155 149变量类型 内部状态 内部状态 内部状态透平机进口压力容器入口温度 Ti_trt10CN 102021254 A CN 102021268 A说明书℃ 148 内部状态8/11 页状态 x4透平机进口压力容器内部温度 Ti_in状态 x5透平机进口压力容器内部压力 Pi_inKpa152内部状态状态 x6透平机出口压力容器温度Ti_out℃141内部状态状态 x7透平机出口压力容器内部压力 Pt_outKpa10内部状态输入 z1 输入 z2 控制 u1 控制 u2 干扰 v1 干扰 v2 输出 y1 输出 y2
高炉压力容器入口煤气流量 高炉压力容器入口煤气温度 静叶开度 旁通阀开度 雷达料线 透平机静叶叶片积灰 炉顶压力 透平机功率G2 T2 Lt Lp H v2 P2 Nm3/h ℃12500 150 45 30关联输入 关联输入 控制输入 控制输入 干扰输入 干扰输入m0.8658Kpa kW160 2885关联输出 本地输出在 TRT 装置生产稳态操作点,由于各种扰动和关联输入的偏差变化,引起高炉 顶压和发电量的波动。 通过偏差变化来判断 TRT 装置是否正常工作。 针对上述模型, 通过泰勒级数展开并取一次项,实现模型线性化,可以得到 TRT 装置的偏差模型,写成 状态空间方程形式为 :
偏差状态空间方程如下 :
其中,偏差量为 ΔX = X-X0, ΔU = U-U0, ΔZ = Z-Z0, Δv = v-v0,初始 时刻,所有变量和状态的偏差都为零,ΔX(0) = 0,ΔU(0) = 0,ΔZ(0) = 0,Δv(0) = 0。
输出方程 :
Y = C′ ΔX+D′ ΔU+E′
通过一阶前向差分11(Ts 为采样周期 ) 将微分方程化为代数CN 102021254 A CN 102021268 A说明书9/11 页方程,实现模型离散化。
离散化偏差状态空间方程如下 :
ΔX(k+1) = AΔX(k)+BΔU(k)+HΔZ(k)+JΔv(k)
输出方程 :
Y(k) = CΔX(k)+DΔU(k)+E
三、高炉和 TRT 装置分布式协调控制
1、确定采样周期
根据高炉系统 S1 和 TRT 装置 S2 的动态响应特性确定它们的采样周期为 T1 和 T2, 其中 T1 > T2,利用前向差分法对高炉系统和 TRT 装置的连续模型进行近似离散化 ;其中 高炉系统 S1 和 TRT 装置 S2 的自身变量以各自的采样周期为准,关联变量的采样方法如 下 :高炉系统的关联输入高炉顶压为 :z1(t) = h1(x2(t)),离散
化后变为k = k0, k0+1, ..., kf, k0 和 kf 为 自 然 数, 且k0 < kf,表示以 T1 为采样周期对 TRT 装置 S2 的状态变量 x2(t) 进行虚拟采样得到的值,得到 S1 的离散模型 M ′ 1,因为 T2 < T1,则不存在采样值缺失的问题 ;同 理, TRT 装置的关联输入为煤气流量和煤气温度 :z2(t) = h2(x1(t)),离散化后变为 k = k0, k0+1, ..., kf, k0 和 kf 为自然数,且 k0 < kf, 表示以 T2 为采样周期对 S1 的状态变量 x1(t) 进行虚拟采样得到的值,得到 S2 的离散模型 M′ 2, 由于 T2 < T1,则在采样不重合的时间点,将出现采样值缺失的情况,即无法获取煤气流 量和煤气温度的数据,缺失的采样值用最近一次采样值代替。
2、模型优化求解方法
采用二次型优化算法对两个模型进行优化求解 ;高炉系统以高炉模型和生产工 艺参数允许范围为约束条件,以高炉内部温度分布与设定值之间的偏差最小以及控制作 用最小为优化目标 ;TRT 装置以 TRT 模型和生产工艺参数允许范围为约束条件,以高炉 顶压稳定、发电量最大和控制作用最小为优化目标,分别对高炉系统和 TRT 装置进行优 化求解,解算出控制量发送给现场控制器。
优化目标函数 :
高炉系统 :内部温度分布稳定,控制作用最小
TRT 装置 :顶压稳定,发电量大,控制作用最小3、可变窗口长度的分布式预测控制方法
采用预测控制的滚动优化方法,在 t = kTi 时刻,对高炉系统 S1 和 TRT 装置 S2 的离散模型在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时间段的系统行为进行预测和优化,其中 Pi 为 Si 在 t = kTi 时刻进行滚动优化的预测时域, Ls = Pi · Ti 为 t = kTi 时刻滚动优化的 时间窗口长度,其中 Si 为 S1 或 S2, Pi 为 P1 或 P2, Ti 为 T1 或 T2 ;当高炉系统与 TRT 装 置协调优化时,预测时间窗口长度由高炉系统 S1 决定。
采用预测控制的模型预测和反馈校正方法,在 t = kTi 时刻,利用 Si 的离散模型
对 Si 在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时间段的系统行为进行预测,利用实测值对预测 结果进行反馈校正,减小偏差。 使用二次型优化算法求解 Si 在未来 t = (k+1)Ti 至 t = (k+Pi)Ti 时间段的最优控制问题 Ji, k,其中 Si 为 S1 或 S2, Pi 为 P1 或 P2, Ti 为 T1 或 T2, 实现可变窗口长度的分布式预测控制。
4、协调控制机制
协调控制机制主要指高炉系统与 TRT 装置优化的协调控制,即协调优化。 协 调优化过程建在对系统未来行为预测的基础上,包含一系列变长度预测窗口内的优化计 算。 每一个预测窗口内的优化计算需迭代多次,直至收敛。 协调优化机制包括 :
启动协调优化 :本地要求启动协调优化 ;关联部分要求启动协调优化 ;
停止协调优化 :每一个预测窗口内,高炉系统和 TRT 装置的控制量都保持稳 定,停止本窗口迭代优化,移动预测窗口 ;当高炉系统与 TRT 装置本地状态都恢复正 常,停止预测窗口的移动,整个协调优化过程结束。
高炉系统和 TRT 装置协调优化步骤 :
1) 到服务器共享缓冲区查看是否有外界协调优化请求 ;
2) 检查本地状态是否需要启动协调优化 ;
3) 如果有外界要求,或本身需要 :
①将响应优化应答或启动优化请求写入服务器共享缓冲区 ;
②到服务器共享缓冲区查看是否其他子系统均已响应 ;
4) 如已具备协调优化条件,则启动优化 :
①到服务器共享缓冲区查看所需的预测窗口长度 ;
②计算本地的优化采样周期 ;
③到服务器共享缓冲区查看给定窗口内所需的关联输入预测值,根据需要进行 插值处理 ;
④计算给定窗口内的本地控制量预测值 ( 优化 )、本地关联输出预测值 ( 模 型);
⑤将给定窗口内的本地关联输出预测值写到服务器共享缓冲区 ;
⑥计算给定窗口内的本地控制量预测值是否满足本地停止窗口内迭代条件 ;
⑦如满足本地停止窗口内迭代条件,将本地停止窗口内迭代要求写到服务器共 享缓冲区 ;
⑧查看其它子系统是否要求停止窗口内迭代 ;如不满足全体停止窗口内迭代条 件,返回到第③步。 否则输出第一个计算时刻控制量到本地系统 ;
⑨计算本地状态是否满足停止窗口移动条件 ;如不满足本地停止窗口移动条 件,移动预测窗口,返回到第①步。 否则将本地停止窗口移动要求写到服务器共享缓冲 区;
⑩查看其它子系统是否要求停止窗口移动 ;如不满足全体停止窗口移动迭代条 件,返回第①步。 否则返回第 1 步。
采用分布式协调控制方法,高炉系统和 TRT 装置通过数据服务器进行数据交 互,读取关联系统的数据信息,分别进行优化求解,大大减小协调优化计算量,从全局 的、整体的角度对高炉系统和 TRT 装置进行协调优化控制,能够有效缓解原料成分波动和其他扰动作用对生产操作和产品质量的影响,减少不必要能源浪费,保证高炉炉内温 度分布正常,炉顶压力稳定, TRT 装置发电量最大,两个系统平稳、高效的运行。