一种基于数据驱动的单元机组锅炉汽轮机控制方法及装置制造方法及图纸

本发明专利技术公开了一种基于数据驱动的单元机组锅炉汽轮机控制方法，包括，建立单元机组非线性模型，采用动态线性化技术得到线性数据模型，估计数据模型中的伪雅各比矩阵，设计基于数据的自适应控制算法及选择控制器参数并验证算法稳定性等五步。本发明专利技术不依赖单元机组的模型信息，避免了复杂工况和运行环境下单元机组的建模问题，同时对于单元机组未建模动态等因素具有较强的鲁棒性。算法设计简单，计算量小，易于实际运行，本发明专利技术所述的技术方案可以被推广应用于其他难以建立系统复杂模型的工业过程被控对象。

【技术实现步骤摘要】

本专利技术涉及一种基于数据驱动的单元机组锅炉汽轮机控制方法及装置，属电力系统和工业过程控制

技术介绍
电力工业是国民经济发展中最重要的基础能源产业，“十三五”期间我国电力需求将逐年增加。为满足国民经济发展需求，实现电力工业可持续发展战略和节能减排目标，采用先进的控制技术和方法提高现有火电厂单元机组的运行效率具有重要的意义。火电厂发电机组普遍采用一台锅炉和一台汽轮机的单元制机组，整个生产流程中能量的传递和转换由锅炉系统和汽轮机系统合作承担。单元机组锅炉汽轮机系统是一类具有典型非线性特征的多变量被控对象，它具有大迟延、强耦合、多不确定干扰等特点。实际运行过程中，机组动态特性除了受其内在物理结构属性影响，还与机组的运行方式、运行工况、外部环境等因素有关。由于上述因素的影响，锅炉汽轮机系统几乎不可能测得运行的所有参数，即使相同的设备、相同的运行工况也会因为某一个操作参数发生变化而影响系统运行过程，因此从物质、能量等角度建立的过程数学模型是不准确的，基于近似模型设计的控制器在实际应用中很难达到期望的控制效果。目前火电厂单元机组协调控制多采用网络化控制系统，系统中每天都产生并存储着大量的过程数据，这些数据隐含着单元机组的运行信息。因此，提出一种仅利用这些大量的离在线数据，在不需要建立锅炉汽轮机系统数学模型的情况下实现对单元机组发电过程的控制和优化，具有重要的实际意义。
技术实现思路
本专利技术的目的是要提供一种基于数据驱动的单元机组锅炉汽轮机控制方法及装置。为达到上述目的，本专利技术是按照以下技术方案实施的：一种基于数据驱动的单元机组锅炉汽轮机控制方法，包括如下步骤：第一步，建立单元机组非线性模型，针对单元机组运行工况，建立以锅炉汽轮机燃油调节阀开度、燃气调节阀开度和给水调节阀开度为控制输入量，汽包压力、功率输出和液位高度变化量为系统输出量的三输入三输出非线性模型；其表达式为：y(k+1)＝f(y(k),u(k)),其中：y＝[y1,y2,y3]T,u＝[u1,u2,u3]T；y1表示汽包压力；y2表示功率输出；y3表示水位高度变化量；f(y(k),u(k))为未知非线性函数；第二步，采用动态线性化技术得到线性模型，针对第一步中建立的三输入三输出非线性模型，采用动态线性化技术将其转化成如下线性数据模型：△y(k+1)＝Φ(k)△u(k),式中△y(k+1)＝y(k+1)-y(k),△u(k)＝u(k)-u(k-1)，Φ(k)为伪雅各比矩阵：且满足||Φ(k)||≤b，b为一个正的常数；第三步，估计数据模型中的伪雅各比矩阵，定义估计准则函数，根据最优化原理得到伪雅各比矩阵的递推估计算法。设定雅各比矩阵估值的初值，根据单元机组实际运行输入输出数据依次更新伪雅各比矩阵的估计值估计准则函数为：式中μ>0为权重因子；为Φ(k)的估计值；雅各比矩阵估计算法：η∈(0,1]为步长因子；为预先设定初始值；第四步，设计基于数据的自适应控制算法，定义控制输入准则函数，根据最优化原理得到控制输入量的更新算法。设定控制输入量的初始值，根据单元机组实际运行输入输出数据以及第三步中得到的雅各比矩阵估计值依次更新控制输入；其中控制输入准则函数为：J(u(k))＝||y*(k+1)-y(k+1)||2+λ||u(k)-u(k-1)||2,；λ>0为权重因子；为期望的汽包压力；为期望的功率输出；为期望的液位高度变化量；控制输入量更新算法为：ρ∈(0,1]为步长因子；为Φ(k)估计值；u(0)为预先设定的控制输入初始值。第五步，选择控制器参数并验证算法稳定性，首先定义跟踪误差，然后将跟踪误差带入到第一步建立的输出动力学数据模型和第三步的自适应控制算法函数中，并得到跟踪误差的递推关系，其中：跟踪误差为：e(k)＝[e1(k)e2(k)e3(k)]；为汽包压力跟踪误差；为功率输出跟踪误差；为液位高度变化跟踪误差；跟踪误差的递推关系为：进一步的，所选择的控制器参数ρ,λ满足如下关系本专利技术不依赖单元机组的模型信息，避免了复杂工况和运行环境下单元机组的建模问题，同时对于单元机组未建模动态等因素具有较强的鲁棒性。算法设计简单，计算量小，易于实际运行，本专利技术所述的技术方案可以被推广应用于其他难以建立系统复杂模型的工业过程被控对象。附图说明下面结合附图和具体实施方式来详细说明本专利技术图1：单元机组数据驱动控制系统原理结构图；图2：单元机组的典型工作点表；图3：单元机组汽包压力跟踪输出曲线；图4：单元机组功率输出跟踪输出曲线；图5：单元机组液位高度变化曲线；图6：本专利技术操作流程图。具体实施方式为使本专利技术实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本专利技术。结合某160MW三输入三输出单元机组为例对本专利技术进行详细说明，如图1-6所示，一种基于数据驱动的单元机组锅炉汽轮机控制方法，包括如下步骤：第一步，建立单元机组非线性模型，针对单元机组运行工况，建立以锅炉汽轮机燃油调节阀开度、燃气调节阀开度和给水调节阀开度为控制输入量，汽包压力、功率输出和液位高度变化量为系统输出量的三输入三输出非线性模型；其表达式为：y(k+1)＝f(y(k),u(k)),其中：y＝[y1,y2,y3]T,u＝[u1,u2,u3]T；y1表示汽包压力；y2表示功率输出；y3表示水位高度变化量；f(y(k),u(k))为未知非线性函数；第二步，采用动态线性化技术得到线性模型，针对第一步中建立的三输入三输出非线性模型，采用动态线性化技术将其转化成如下线性数据模型：△y(k+1)＝Φ(k)△u(k),式中△y(k+1)＝y(k+1)-y(k),△u(k)＝u(k)-u(k-1)，Φ(k)为如下伪雅各比矩阵：且满足||Φ(k)||≤b，b为一个正的常数；第三步，估计数据模型中的伪雅各比矩阵，定义估计准则函数，根据最优化原理得到伪雅各比矩阵的递推估计算法。设定雅各比矩阵估值的初值，根据单元机组实际运行输入输出数据依次更新伪雅各比矩阵的估计值估计准则函数为：式中μ>0为权重因子；为Φ(k)的估计值；雅各比矩阵估计算法：η∈(0,1]为步长因子；为预先设定初始值；第四步，设计基于数据的自适应控制算法，定义控制输入准则函数，根据最优化原理得到控制输入量的更新算法。设定控制输入量的初始值，根据单元机组实际运行输入输出数据以及第三步中得到的雅各比矩阵估计值依次更新控制输入；其中控制输入准则函数为：J(u(k))＝||y*(k+1)-y(k+1)||2+λ||u(k)-u(k-1)||2,；λ>0为权重因子；为期望的汽包压力；为期望的功率输出；为期望的液位高度变化量；控制输入量更新算法为：ρ∈(0,1]为步长因子；为Φ(k)估计值；u(0)为预先设定的控制输入初始值。第五步，选择控制器参数并验证算法稳定性，首先定义跟踪误差，然后将跟踪误差带入到第一步建立的输出动力学数据模型和第三步的自适应控制算法函数中，并得到跟踪误差的递推关系，其中：跟踪误差为：e(k)＝[e1(k)e2(k)e3(k)]；为汽包压力跟踪误差；为功率输出跟踪误差；为液位高度变化跟踪误差；跟踪误差的递推关系为：本实施例中，所选择的控制器参数ρ,λ满足如下关系本专利技术在实际实施时，按照以本文档来自技高网...

【技术保护点】
一种基于数据驱动的单元机组锅炉汽轮机控制方法，其特征在于，所述的基于数据驱动的单元机组控制方法包括如下步骤：第一步，建立单元机组非线性模型，针对单元机组运行工况，建立以锅炉汽轮机燃油调节阀开度、燃气调节阀开度和给水调节阀开度为控制输入量，汽包压力、功率输出和液位高度变化量为系统输出量的三输入三输出非线性模型；其表达式为：y(k+1)＝f(y(k),u(k)),其中：y＝[y1,y2,y3]T,u＝[u1,u2,u3]T；y1表示汽包压力；y2表示功率输出；y3表示水位高度变化量；f(y(k),u(k))为未知非线性函数；第二步，采用动态线性化技术得到线性模型，针对第一步中建立的三输入三输出非线性模型，采用动态线性化技术将其转化成如下线性数据模型：△y(k+1)＝Φ(k)△u(k),式中△y(k+1)＝y(k+1)‑y(k),△u(k)＝u(k)‑u(k‑1)，Φ(k)为伪雅各比矩阵：且满足||Φ(k)||≤b，b为一个正的常数；第三步，估计数据模型中的伪雅各比矩阵，定义估计准则函数，根据最优化原理得到伪雅各比矩阵的递推估计算法。设定雅各比矩阵估值的初值，根据单元机组实际运行输入输出数据依次更新伪雅各比矩阵的估计值估计准则函数为：式中μ>0为权重因子；为Φ(k)的估计值；雅各比矩阵估计算法：η∈(0,1]为步长因子；为预先设定初始值；第四步，设计基于数据的自适应控制算法，定义控制输入准则函数，根据最优化原理得到控制输入量的更新算法。设定控制输入量的初始值，根据单元机组实际运行输入输出数据以及第三步中得到的雅各比矩阵估计值依次更新控制输入；其中控制输入准则函数为：J(u(k))＝||y*(k+1)‑y(k+1)||2+λ||u(k)‑u(k‑1)||2,；λ>0为权重因子；为期望的汽包压力；为期望的功率输出；为期望的液位高度变化量；控制输入量更新算法为：ρ∈(0,1]为步长因子；为Φ(k)估计值；u(0)为预先设定的控制输入初始值。第五步，选择控制器参数并验证算法稳定性，首先定义跟踪误差，然后将跟踪误差带入到第一步建立的输出动力学数据模型和第三步的自适应控制算法函数中，并得到跟踪误差的递推关系，其中：跟踪误差为：e(k)＝[e1(k) e2(k) e3(k)]；为汽包压力跟踪误差；为功率输出跟踪误差；为液位高度变化跟踪误差；跟踪误差的递推关系为：...

【技术特征摘要】
1.一种基于数据驱动的单元机组锅炉汽轮机控制方法，其特征在于，所述的基于数据驱动的单元机组控制方法包括如下步骤：第一步，建立单元机组非线性模型，针对单元机组运行工况，建立以锅炉汽轮机燃油调节阀开度、燃气调节阀开度和给水调节阀开度为控制输入量，汽包压力、功率输出和液位高度变化量为系统输出量的三输入三输出非线性模型；其表达式为：y(k+1)＝f(y(k),u(k)),其中：y＝[y1,y2,y3]T,u＝[u1,u2,u3]T；y1表示汽包压力；y2表示功率输出；y3表示水位高度变化量；f(y(k),u(k))为未知非线性函数；第二步，采用动态线性化技术得到线性模型，针对第一步中建立的三输入三输出非线性模型，采用动态线性化技术将其转化成如下线性数据模型：△y(k+1)＝Φ(k)△u(k),式中△y(k+1)＝y(k+1)-y(k),△u(k)＝u(k)-u(k-1)，Φ(k)为伪雅各比矩阵：且满足||Φ(k)||≤b，b为一个正的常数；第三步，估计数据模型中的伪雅各比矩阵，定义估计准则函数，根据最优化原理得到伪雅各比矩阵的递推估计算法。设定雅各比矩阵估值的初值，根据单元机组实际运行输入输出数据依次更新伪雅各比矩阵的估计值估计准则函数为：式...

【专利技术属性】
技术研发人员：卜旭辉，崔立志，余琼霞，尹艳玲，王森，
申请(专利权)人：河南理工大学，
类型：发明
国别省市：河南;41