变工况亚临界火电机组的锅炉‑汽机控制模型构建方法技术

本发明专利技术提供了一种变工况亚临界火电机组的锅炉‑汽机控制模型构建方法，所述建模步骤包括：1、获得亚临界火电机组的结构参数，建立适用于亚临界机组的工质物性参数数据库；2、根据燃煤机组内部能量转换与传递过程，将锅炉‑汽机系统全流程简化为给煤‑热量转化过程、热量‑蒸汽转化过程和蒸汽‑电功率转化过程，基于工质的质量、能量动态平衡，分别建立这三个过程的集总参数动态模型；3、从DCS实时数据库获取机组的历史运行数据；4、利用获取的机组DCS历史运行数据，对模型的中间变量进行回归分析，通过参数辨识得到模型参数。本发明专利技术可以为变工况运行的机组动态特性仿真分析、协调控制系统设计等提供技术支撑。

【技术实现步骤摘要】

本专利技术涉及火力发电控制领域的一种变工况亚临界火电机组的锅炉-汽机控制模型构建方法，具体地，涉及一种基于锅炉-汽轮机系统机理建模和动态参数辨识的亚临界火电机组协调控制模型构建方法。
技术介绍
随着我国电网的日益复杂和电力需求的多样化，对负荷和频率调节提出了越来越高的要求。对于火力电站而言，协调控制系统是影响机组安全经济运行的重要因素。但由于锅炉-汽机系统的复杂性，大多基于简单线性模型设计的协调控制系统无法达到满意的控制效果。改善机组协调控制系统设计的一个重要手段是建立能反映机组全工况动态特性的非线性动态过程模型。建立一个复杂度和精度适中的锅炉-汽机控制模型，对于进一步改善协调控制系统设计和分析机组动态特性具有重要意义。目前常见的协调控制模型是线性传递函数模型，但是线性模型忽略了锅炉-汽机系统本身的非线性和耦合特性，因此只适用于机组运行工况变化不大的情况。对现有技术的检索发现，韩忠旭(韩忠旭等.机炉协调多变量受控系统线性增量形数学模型.中国电机工程学报.25(7),24-29,2005)提出了一种火电厂单元机组线性增量式集总参数模型，推导出了系统输入与输出之间的传递函数矩阵。然而针对需要参与调峰的大型火电机组，其通常需要在30％～100％额定负荷的变工况条件下运行，故线性模型无法表征锅炉-汽机系统的全局动态特性；同时文献模型没有考虑到模型参数的时变性，部分模型参数物理意义不明确，因此该模型缺乏可移植性，尤其不适用于工况变化大的亚临界机组。
技术实现思路
针对现有技术的缺陷，本专利技术的目的是提供一种变工况亚临界火电机组的锅炉-汽机控制模型构建方法，充分考虑了变工况运行的锅炉-汽机系统的非线性特性，对锅炉机组的各个阶段能量转换和工质参数变化进行了机理分析，通过适当的模型简化和参数辨识，推导出复杂度适中的机组协调控制模型。为实现以上目的，本专利技术提供一种变工况亚临界火电机组的锅炉-汽机控制模型构建方法，包括以下步骤：步骤一、根据锅炉、汽轮机运行设计规程，获得亚临界火电机组(包括水冷壁、汽包、各级换热器、)以下结构参数：汽包、水冷壁和过热器的管径、管厚、管长及管数量。建立可在线调用的适用于亚临界机组的工质物性参数(包括工质的密度、比焓)数据库，用于步骤二中的工质参数计算。优选地，所述工质物性参数库，是指根据水和水蒸气热力性质工业公式(IAPWS-IF97)开发的用于在线计算的工质物性参数数据库，可参考文献：王旭辉、于彤、惠兆宇、袁景淇，用于火电全范围仿真的工质物性参数数据库，控制工程，2011；18:131-133。步骤二、为了便于分析，根据燃煤机组内部能量转换与传递过程，将锅炉-汽机系统全流程(见图1)简化为三个子过程，即给煤-热量转化过程、热量-蒸汽转化过程和蒸汽-电功率转化过程。基于工质的质量、能量动态平衡。分别建立这三个子过程的集总参数动态模型，建模流程见图2，建模过程中有关工质参数的计算采用步骤一中的工质物性参数数据库：①给煤-热量转化过程：是指原煤经过磨煤机后磨成煤粉，并随一次风送入炉膛燃烧释放热量的过程。其中，进入炉膛的煤粉质量动态平衡方程为TfdDfdt=-Df+μB(t-τ)---(1)]]>炉膛内燃烧过程的动态能量平衡方程为Qs＝DfQnet,ar(1-q4)-Dfg(hfg-hin)＝DfQnet,ar(1-q4)-Qfg(2)其中，烟气能量与发电功率之间的关系用线性多项式来逼近：Qfg＝λ10+λ11N(3)式(1)～(3)中，μB、Df分别表示给煤质量流量和入炉煤质量流量，Tf和τ表示时间常数；Qs表示水冷壁内工质的吸热速率，q4表示机械不完全燃烧损失，Qnet,ar表示燃煤的低位发热量，Qfg、Dfg、hfg分别表示炉膛出口烟气携带的能量、炉膛出口烟气的质量流量和比焓，hin表示烟气在环境温度下的比焓；λ10、λ11为多项式函数的常系数；N是发电机的发电功率。②热量-蒸汽转化过程：是指工质在水冷壁和各级过热器吸收煤粉燃烧所释放出的热量，沿着汽水流程工质参数逐步升高最终形成主蒸汽的过程。其中，水冷壁和汽包中工质的质量、能量平衡方程为Dgs-Ds=d(VDρD)dt=VDdρDdpDdpDdt---(4)]]>Qs＝Dgs(hs-hfw)(5)其中，汽包内饱和蒸汽和给水的热力学参数可用二次多项式表示ρD=λ20+λ21pD+λ22pD2---(6)]]>hs=λ30+λ31pD+λ32pD2---(7)]]>hfw=λ40+λ41pD+λ42pD2---(8)]]>过热器内工质的质量动态平衡方程为Ds-Dms=d(VTρT)dt=VTdρTdt=VT∂ρT∂pTdpTdt---(9)]]>当主蒸汽温度变化范围很小时，主蒸汽密度与主蒸汽压力的关系可用线性函数表示：ρT＝λ50+λ51pT(10)根据Darcy-Weisbach方程，可计算离开汽包的饱和蒸汽质量流量Ds=C0ρDpD-pT---(11)]]>由式(9)和(11)可推导出主蒸汽压力的动态特性方程为VT∂ρT∂pTdpTdt=C0ρDpD-pT-Dms---(12)]]>式(4)～(12)中，Dgs、Ds分别表示水冷壁和汽包中产生的饱和蒸汽质量流量和离开汽包的饱和蒸汽质量流量，VD、ρD、hs、pD分别表示汽包和水冷壁中饱和蒸汽的体积、饱和蒸汽密度、饱和蒸汽比焓和汽包压力，hfw表示给水比焓；Dms、ρT、pT表示主蒸汽质量流量、主蒸汽密度、主蒸汽压力，VT是过热器中过热蒸汽的总体积；C0是过热器管道惯性阻尼系数；Qs是指水冷壁中工质的吸热速率；λ20、λ21、λ22、λ30、λ31、λ32、λ40、λ41、λ42、λ50、λ51均为多项式函数的常系数，用于表达中间量与系统状态变量之间的关系。③蒸汽-电功率转化过程：是指主蒸汽进入汽轮机高压缸做功，然后高压缸排气流经再热器吸收热量形成再热蒸汽并进入汽轮机中、低压缸做功的过程。其中，汽轮机调节阀门的开度定义为μT=p1pT---(13)]]>汽机阀门的动态特性简化为一个一阶惯性环节Tedp1dt+p1=μTpT---(14)]]>根据进出汽轮机系统的总能量平衡，发电机的发电功率可计算如下N=ηmηgEin=ηmηg(Dmshms+Dr(hr-hHE)-Σi=18Dihi-DLEhLE)=ηmηg[Dms(hms-Σi=19αihi)+Dr(hr-hHE)]---(15)]]>αi＝Di/Dms,i＝1,2...9(16)其中，再热蒸汽在再热通道中吸收的总热量与发电功率之间存在线性相关，可表示如下：Dr(hr-hHE)＝k1N(17)因而式(15)可改写为N=(hms-Σi=19αihi)Dms1ηmηg-k1=k2Dms---(18)]]>其中，主蒸汽流量可按下式计算Dms本文档来自技高网...

【技术保护点】
一种变工况亚临界火电机组的锅炉‑汽机控制模型构建方法，其特征在于，所述模型的建模步骤包括：步骤一、根据锅炉、汽轮机运行设计规程，获得亚临界火电机组以下结构参数：汽包、水冷壁和过热器各自的内体积，建立可在线调用的适用于亚临界机组的工质物性参数数据库，所述工质物性参数包括工质的密度、比焓，用于步骤二中的工质参数计算；步骤二、根据燃煤机组内部能量转换与传递过程，将锅炉‑汽机系统全流程简化为三个子过程，即给煤‑热量转化过程、热量‑蒸汽转化过程和蒸汽‑电功率转化过程，基于工质的质量、能量动态平衡，分别建立这三个子过程的集总参数动态模型；①给煤‑热量转化过程：是指原煤经过磨煤机后形成煤粉，并随一次风被送入炉膛燃烧释放热量的过程；②热量‑蒸汽转化过程：是指工质在水冷壁和各级过热器中吸收燃烧过程所释放出的热量，工质参数沿着汽水流程逐步升高最终形成主蒸汽的过程；③蒸汽‑电功率转化过程：是指主蒸汽进入汽轮机高压缸做功，然后高压缸排气流经再热器吸收热量形成再热蒸汽并再次进入汽轮机中、低压缸做功的过程；以汽机调门开度和给煤质量流量作为系统输入，发电功率和主蒸汽压力作为系统输出，以发电功率、主蒸汽压力、汽包压力和入炉煤的质量流量作为系统状态变量，建立适用于变工况亚临界火电机组的锅炉‑汽机系统控制模型，建模过程中有关工质参数的计算采用步骤一中的工质物性参数数据库；步骤三、通过OPC协议从DCS实时数据库获取机组的历史运行数据，包括：给煤质量流量、主蒸汽压力、主蒸汽温度、主蒸汽质量流量、再热蒸汽压力、再热蒸汽温度、再热蒸汽质量流量、汽轮机首级压力、汽包压力、发电功率、给水温度、给水压力、烟气温度、烟气质量流量；步骤四、在步骤三的基础上，利用一段时间的DCS历史运行数据，对步骤二中锅炉‑汽机系统控制模型的中间变量进行回归分析，确定中间变量与系统状态变量之间的多项式函数关系式，通过参数辨识得到模型参数。...

【技术特征摘要】
1.一种变工况亚临界火电机组的锅炉-汽机控制模型构建方法，其特征在于，所述模型的建模步骤包括：步骤一、根据锅炉、汽轮机运行设计规程，获得亚临界火电机组以下结构参数：汽包、水冷壁和过热器各自的内体积，建立可在线调用的适用于亚临界机组的工质物性参数数据库，所述工质物性参数包括工质的密度、比焓，用于步骤二中的工质参数计算；步骤二、根据燃煤机组内部能量转换与传递过程，将锅炉-汽机系统全流程简化为三个子过程，即给煤-热量转化过程、热量-蒸汽转化过程和蒸汽-电功率转化过程，基于工质的质量、能量动态平衡，分别建立这三个子过程的集总参数动态模型；①给煤-热量转化过程：是指原煤经过磨煤机后形成煤粉，并随一次风被送入炉膛燃烧释放热量的过程；②热量-蒸汽转化过程：是指工质在水冷壁和各级过热器中吸收燃烧过程所释放出的热量，工质参数沿着汽水流程逐步升高最终形成主蒸汽的过程；③蒸汽-电功率转化过程：是指主蒸汽进入汽轮机高压缸做功，然后高压缸排气流经再热器吸收热量形成再热蒸汽并再次进入汽轮机中、低压缸做功的过程；以汽机调门开度和给煤质量流量作为系统输入，发电功率和主蒸汽压力作为系统输出，以发电功率、主蒸汽压力、汽包压力和入炉煤的质量流量作为系统状态变量，建立适用于变工况亚临界火电机组的锅炉-汽机系统控制模型，建模过程中有关工质参数的计算采用步骤一中的工质物性参数数据库；步骤三、通过OPC协议从DCS实时数据库获取机组的历史运行数据，包括：给煤质量流量、主蒸汽压力、主蒸汽温度、主蒸汽质量流量、再热蒸汽压力、再热蒸汽温度、再热蒸汽质量流量、汽轮机首级压力、汽包压力、发电功率、给水温度、给水压力、烟气温度、烟气质量流量；步骤四、在步骤三的基础上，利用一段时间的DCS历史运行数据，对步骤二中锅炉-汽机系统控制模型的中间变量进行回归分析，确定中间变量与系统状态变量之间的多项式函数关系式，通过参数辨识得到模型参数。2.根据权利要求1所述的一种变工况亚临界火电机组的锅炉-汽机控制模型构建方法，其特征在于，步骤二中所述的锅炉-汽机系统控制模型包括：①给煤-热量转化过程，其中：进入炉膛的煤粉质量动态平衡方程为：TfdDfdt=-Df+μB(t-τ)---(1)]]>炉膛内燃烧过程的动态能量平衡方程为：Qs＝DfQnet,ar(1-q4)-Dfg(hfg-hin)＝DfQnet,ar(1-q4)-Qfg(2)烟气能量与发电功率之间的关系用线性多项式来逼近：Qfg＝λ10+λ11N(3)式(1)～(3)中，μB、Df分别表示给煤质量流量和入炉煤质量流量，Tf和τ表示时间常数；Qs表示水冷壁内工质的吸热速率，q4表示机械不完全燃烧损失，Qnet,ar表示燃煤的低位发热量，Qfg、Dfg、hfg分别表示炉膛出口烟气携带的能量、炉膛出口烟气的质量流量和比焓，hin表示烟气在环境温度下的比焓；λ10、λ11为多项式函数的常系数；N是发电机的发电功率；②热量-蒸汽转化过程，其中：水冷壁和汽包中工质的质量、能量平衡方程为：Dgs-Ds=d(VDρD)dt=VDdρDdpDdpDdt---(4)]]>Qs＝Dgs(hs-hfw)(5)其中，汽包内饱和蒸汽和给水的热力学参数可用二次多项式表示：ρD=λ20+λ21pD+λ22pD2---(6)]]>hs=λ30+λ31pD+λ32pD2---(7)]]>hfw=λ40+λ41pD+λ42pD2---(8)]]>过热器内工质的质量动态平衡方程为：Ds-Dms=d(VTρT)dt=VTdρTdt=VT∂ρT∂pTdpTdt---(9)]]>当主蒸汽温度变化范围很小时，主蒸汽密度与主蒸汽压力的关系用线性函数表达：ρT＝λ50+λ51pT(10)根据Darcy-Weisbach方程，计算离开汽包的饱和蒸汽质量流量：Ds=C0ρDpD-pT---(11)]]>由式(9)和(11)推导出主蒸汽压力的动态特性方程为：VT∂ρT∂pTdpTdt=C0ρDpD-pT-Dms---(12)]]>式(4)～(12)中，Dgs、Ds分别表示水冷壁和汽包中产生的饱和蒸汽质量流量和离开汽包的饱和蒸汽质量流量，VD、ρD、hs、pD分别表示汽包和水冷壁中饱和蒸汽的体积、饱和蒸汽密度、饱和蒸汽比焓和汽包压力，hfw表示给水比焓；Dms、ρT、pT表示主蒸汽质量流量、主蒸汽密度、主蒸汽压力，VT是过热器中过热蒸汽的总体积；C0是过热器管道惯性阻尼系数；Qs是指水冷壁中工质的吸热速率；λ20、λ21、λ22、λ30、λ31、λ32、λ40、λ41、λ42、λ50、λ51均为多项式函数的常系数，用于表达中间量与系统状态变量之间的关系；③蒸汽-电功率转化过程，其中：汽轮机调节阀门的开度定义为：μT=p1pT---(13)]]>汽机阀门的动态特性简化为一个一阶惯性环...

【专利技术属性】
技术研发人员：田震，袁景淇，徐亮，张翔，孔磊，王润杉，王景成，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海;31