一种过热汽温快速节能预测控制方法技术

本发明专利技术公开了一种过热汽温快速节能预测控制方法，该方法通过过热汽温对象的阶跃试验，获得被控对象预测控制模型，过热汽温控制系统在线实时记录减温喷水阀门开度控制量序列和过热汽温测量值序列，控制系统根据本发明专利技术所提出的新型性能指标，通过简单计算得出当前减温喷水阀门开度。此控制方法解决了传统受限预测控制中的优化问题复杂，计算量大，且缺乏有效的控制参数对系统控制性能进行调节的问题，能使过热汽温控制稳定、快速，过程动态偏差小，且无稳态偏差。该方法在保证控制品质的同时还能够减少控制过程中的减温喷水量，提高机组运行的经济性。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了，该方法通过过热汽温对象的阶跃试验，获得被控对象预测控制模型，过热汽温控制系统在线实时记录减温喷水阀门开度控制量序列和过热汽温测量值序列，控制系统根据本专利技术所提出的新型性能指标，通过简单计算得出当前减温喷水阀门开度。此控制方法解决了传统受限预测控制中的优化问题复杂，计算量大，且缺乏有效的控制参数对系统控制性能进行调节的问题，能使过热汽温控制稳定、快速，过程动态偏差小，且无稳态偏差。该方法在保证控制品质的同时还能够减少控制过程中的减温喷水量，提高机组运行的经济性。【专利说明】
本专利技术属于热工自动控制
，具体涉及。
技术介绍
过热汽温是火电厂热力系统中的重要参数，过热汽温控制品质的优劣直接影响到整个机组的安全和经济运行，过热汽温控制系统是机组的重要控制系统之一。由于过热汽温对象是大滞后大延迟热工对象，并且随着机组容量的增加，这一特点会更加显著。传统PID控制算法往往难以满足即快又稳的控制要求，影响机组安全经济运行。因此研究设计先进节能的过热汽温控制算法具有特别重要的意义。 预测控制算法是使用过程模型来控制对象未来行为的一类计算机算法，是工业过程控制中最具价值的先进控制策略之一。但是传统的预测控制，系统算法复杂，计算量大，难以在实际工程中实现在线实时控制。同时传统预测控制，缺乏有效的控制参数对系统控制性能进行调节，不利于工程应用。 为了提高过热汽温控制品质，本专利技术在深入研究传统预测控制算法的基础上，针对过热汽温系统的特性，提出了一种新的预测控制性能指标，并根据该性能指标设计新的快速节能预测控制算法。该算法计算量小，控制性能调整方便有效，在保证控制效果的同时能够减少控制过程中的减温喷水量。将该控制算法应用于过热汽温控制系统中，得到了很好的控制效果。
技术实现思路
专利技术目的:针对上述问题，提出，解决传统预测控制方法算法复杂，计算量大和缺乏有效的控制参数对系统控制性能进行有效调节，难于工程实现的缺点。将该控制方法应用于过热汽温控制系统，能够实现稳定快速无偏差的过程控制，并能通过控制参数对控制系统的控制性能进行有效调节，同时在控制过程中能够减少减温喷水量，提高机组的经济性。 技术方案:本专利技术依据过热汽温系统的特性，提出了一种新的预测控制性能指标，该指标主要考虑预测时域最后一步的预测偏差和预测偏差变化率的情况，同时包含减温喷水量的优化项。基于该指标设计预测控制器，能够实现快速稳定控制，并能通过调节控制参数对控制系统的控制性能进行有效调节；同时能够减少控制过程中的减温喷水量，提高机组经济性。 ，包括如下步骤: 步骤1:获取过热汽温对象模型，在稳态的工况下，作减温喷水阀门开度的开环阶跃响应试验，获取系统响应数据，使用阶跃响应模型辨识方法，辨识得到减温喷水阀门开度一过热汽温的传递函数模型； 步骤2:选择采样周期T，T首先满足香农定理，然后使得T95/T取值在10到25之间，T95为过热汽温过渡过程上升到95%的调节时间。对步骤I中的过热汽温传递函数模型进行离散化，得到其CARIMA模型为“(z-OyGO = B (z—1) u (k_l) + ε (k)/A，式中 【权利要求】1.，其步骤如下: 步骤1:获取过热汽温对象模型，在稳态的工况下，作减温喷水阀门开度的开环阶跃响应试验，获取系统响应数据，使用阶跃响应模型辨识方法，辨识得到减温喷水阀门开度一过热汽温的传递函数模型； 步骤2:选择采样周期T，T首先满足香农定理，然后使得T95/T取值在10到25之间，T95为过热汽温过渡过程上升到95%的调节时间。对步骤I中的过热汽温传递函数模型进行离散化，得到其CARIMA模型为:A (z—1) y(k) =B (z—1) u (k_l) + ε (k) / Λ，式中.傘s) =，z-1为后移算子，Δ = l-z_1,y(k)为k时刻的过热汽温测量值，u(k-l)为k-Ι时刻减温喷水阀门开度控制量，ε (k)为白噪声，Bi,匕为多项式A(P) H中^的系数，na，nb分别为多项式A(^BM)的阶次； 步骤3:根据丟番图方程I = Ej (z-1) A (ζ-1) Δ +Z-jFj (ζ-1)，解得控制变量多项式:,式中 \ 分别为多项式 EJz-1) ,FJz-1)中 ζ-1 的系数，j = N-ι, N, N为预测时域长度，其选择需大于过热汽温系统被控对象的纯延迟时间； 步骤4:根据算SGj (ζ—1) =Ej (P)B(P),求得多项式且令多项式为多项式Gj(z’中ζ"的系数，hJ；i为多项式HJz-1)中的系数，j = N-1, N ； 步骤5:控制系统在线实时记录该控制系统从开始运行至采样k时刻的控制信号{u (k-1)，u (k-2)，…u (I)}，过热汽温的测量数据{y (k)，y (k-1)，…y (I)}，同时记录控制量的增量{Au(k-l)，Au(k-2)，…AU(1)}，并接收预测控制系统的过热汽温期望值yjk);计算得到中间变量步骤6:按下式对过热汽温系统未来第k+Ν-Ι和k+N时刻的输出进行预测:y (k+j I k)=Guj.Δ U+fj (j = N-l.N)，式中y (k+j | k)为未来k+j时刻过热汽温预测值，Guj为一维行向量,Guj= ，Δυ = ,其中 Nu 为控制时域长度，取之间的数值，一般取I或2，Au(k+i), i = O…Nu-1表示在k时刻对未来k+i时刻减温喷水阀门开度增量的估计值； 步骤7:根据新型性能指标来求取减温喷水阀门开度控制量增量Au(k): 式中e (k+N) = y (k+N I k) -yr (k)为第k+N时刻的预测偏差，对应e (k+N-1)为k+N-1时刻的预测偏差；ce(k+N) =e(k+N)-e (k+N-1)为k+N时刻的预测偏差变化率；u(k+i)为未来k+i时刻减温喷水阀门开度的估计值；将步骤6中的未来预测值带入上述性能指标中，并通过极值必要条件求得: AU= (Ρτ.P+QT.Rt.R.Q) ―1.Pt. 式中，P = Gun- λ ? (Gun-Gu^1)，入 ι 和 λ 2 为控制器参数，根据控制系统对系统稳定性和快速性的综合要求选择，入1取区间内数值，A2取区间内数值； 步骤8:计算当前k时刻控制增量Au(k) = Θ.AU, Θ = 1XNu ;由当前k时刻的减温喷水阀门开度控制增量的约束上下限，计算该时刻实际控制增量的约束:取Λ Umin和Umin-U (k-Ι)之间的大者记为Ua，取Δ Umax和Umax-U (k-Ι)之间的小者记为Ub ;如果Au(k)小于Ua则令Au(k) =Ua,如果Au(k)大于ub,则令Δ u (k) = Ub ； 步骤9:令u(k) =u(k-l) + Au(k);其中，u(k)即为当前预测控制器的输出控制量。输出u (k)，在之后的每个采样周期内重复执行步骤5到步骤9。2.根据权利要求1所述的过热汽温快速节能预测控制方法，其特征在于:该预测控制系统预测时域长度N选择上只需大于纯延迟时间；同时控制时域长度Nu取之间内的数值，一般取I或2; 当Nu取I时，所述向量D，P, 本文档来自技高网...

【技术保护点】
一种过热汽温快速节能预测控制方法，其步骤如下：步骤1：获取过热汽温对象模型，在稳态的工况下，作减温喷水阀门开度的开环阶跃响应试验，获取系统响应数据，使用阶跃响应模型辨识方法，辨识得到减温喷水阀门开度—过热汽温的传递函数模型；步骤2：选择采样周期T，T首先满足香农定理，然后使得T95/T取值在10到25之间，T95为过热汽温过渡过程上升到95％的调节时间。对步骤1中的过热汽温传递函数模型进行离散化，得到其CARIMA模型为：A(z‑1)y(k)＝B(z‑1)u(k‑1)+ε(k)/Δ，式中z‑1为后移算子，Δ＝1‑z‑1，y(k)为k时刻的过热汽温测量值，u(k‑1)为k‑1时刻减温喷水阀门开度控制量，ε(k)为白噪声，ai、bi为多项式A(z‑1)，B(z‑1)中z‑i的系数，na，nb分别为多项式A(z‑1)，B(z‑1)的阶次；步骤3：根据丟番图方程1＝Ej(z‑1)A(z‑1)Δ+z‑jFj(z‑1)，解得控制变量多项式：式中ej,i、fj,i分别为多项式Ej(z‑1)，Fj(z‑1)中z‑i的系数，j＝N‑1,N，N为预测时域长度，其选择需大于过热汽温系统被控对象的纯延迟时间；步骤4：根据算式Gj(z‑1)＝Ej(z‑1)B(z‑1)，求得多项式且令多项式Hj(z-1)=zj-1(Gj(z-1)-gj,j-1z-(j-1)-...-gj,0)=Σi=1nbhj,iz-i;]]>式中gj,i为多项式Gj(z‑1)中z‑i的系数，hj,i为多项式Hj(z‑1)中z‑i的系数，j＝N‑1,N；步骤5：控制系统在线实时记录该控制系统从开始运行至采样k时刻的控制信号{u(k‑1),u(k‑2),…u(1)}，过热汽温的测量数据{y(k),y(k‑1),…y(1)}，同时记录控制量的增量{Δu(k‑1),Δu(k‑2),…Δu(1)}，并接收预测控制系统的过热汽温期望值yr(k)；计算得到中间变量fj=Σi=0nafjiy(k-i)+Σi=1nbhjiΔu(k-i),]]>其中j＝N‑1,N；步骤6：按下式对过热汽温系统未来第k+N‑1和k+N时刻的输出进行预测：y(k+j|k)＝Guj·ΔU+fj(j＝N‑1.N)，式中y(k+j|k)为未来k+j时刻过热汽温预测值，Guj为一维行向量，Guj＝[gj,j‑1,gj,j‑2,…,gj,j‑Nu]，ΔU＝[Δu(k),…,Δu(k+Nu‑1)]，其中Nu为控制时域长度，取[1，5]之间的数值，一般取1或2，Δu(k+i),i＝0…Nu‑1表示在k时刻对未来k+i时刻减温喷水阀门开度增量的估计值；步骤7：根据新型性能指标J=[e(k+N)+λ1ce(k+N)]2+λ2Σi=1Nu[u(k+i)-u(k-1)]2]]>来求取减温喷水阀门开度控制量增量Δu(k)：式中e(k+N)＝y(k+N|k)‑yr(k)为第k+N时刻的预测偏差，对应e(k+N‑1)为k+N‑1时刻的预测偏差；ce(k+N)＝e(k+N)‑e(k+N‑1)为k+N时刻的预测偏差变化率；u(k+i)为未来k+i时刻减温喷水阀门开度的估计值；将步骤6中的未来预测值带入上述性能指标中，并通过极值必要条件求得：ΔU＝(PT·P+QT·RT·R·Q)‑1·PT·[yr(k)‑fN‑λ1·(fN‑fN‑1)]式中，P＝GuN‑λ1(GuN‑GuN‑1)，λ1和λ2为控制器参数，根据控制系统对系统稳定性和快速性的综合要求选择，λ1取[0，5]区间内数值，λ2取[0，1]区间内数值；步骤8：计算当前k时刻控制增量Δu(k)＝θ·ΔU，θ＝[1 0…0]1×Nu；由当前k时刻的减温喷水阀门开度控制增量的约束上下限[Δumin,Δumax]和控制量的约束上下限[umin,umax]，计算该时刻实际控制增量的约束：取Δumin和umin‑u(k‑1)之间的大者记为ua，取Δumax和umax‑u(k‑1)之间的小者记为ub；如果Δu(k)小于ua则令Δu(k)＝ua，如果Δu(k)大于ub，则令Δu(k)＝ub；步骤9：令u(k)＝u(k‑1)+Δu(k)；其中，u(k)即为当前预测控制器的输出控制量。输出u(k)，在之后的每个采样周期内重复执行步骤5到步骤9。...

【技术特征摘要】

【专利技术属性】
技术研发人员：雎刚，陈欢乐，陈雪峰，
申请(专利权)人：南京归图科技发展有限公司，
类型：发明
国别省市：江苏;32