发明名称 基于转子温度场模拟计算的电站锅炉空预器热点检测方法
摘要 本发明公开了一种基于转子温度场模拟计算的电站锅炉空预器热点检测方法,该方法按以下步骤实施:确定所要监控的空预器的参数,并设置温度参考值;实时采集各种参数并将所得参数值输入到计算模型中;边界温度值初始化,并确定边界条件;用一个差分方程,得到该段、该分仓内部及其出口边界上流体和受热面的温度;若迭代过程前后两次得到的温度误差符合设定的温度参考值的要求,则计算结束,否则更新边界温度值重新迭代。当现有热点检测系统的红外热点传感器或执行机构故障,系统发出切换命令,用本发明方法模拟计算得到的温度场参数作为控制的参考值,提高了空预器热点检测装置的可靠性,保证该检测系统安全稳定运行。
申请公布号 CN101382458B 申请公布日期 2010.11.10
申请号 CN200810231772.0 申请日期 2008.10.16
申请人 西安理工大学 发明人 刘丁;刘涵;杨延西
分类号 G01K13/00(2006.01)I 主分类号 G01K13/00(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 罗笛
主权项 1.一种基于转子温度场模拟计算的电站锅炉空预器热点检测方法,其特征在于,该方法按以下步骤实施:步骤1:利用现有的空预器热点检测装置,确定所要监控的空预器的结构参数:半径R(m)、转子高度h(m)、旋转角速度ω(rad/s),设置该空预器热点检测装置的温度参考值;步骤2:实时采集以下参数:流速、密度、比热、热传导率、换热系数、蓄热元件单位体积上的换热面积、温度;并将所得参数值输入到方程式(1)和方程式(2)中:<img file="FA20177972200810231772001C00011.GIF" wi="1094" he="162" /><img file="FA20177972200810231772001C00012.GIF" wi="1589" he="180" />(1)<img file="FA20177972200810231772001C00013.GIF" wi="1156" he="157" />(2)其中,v为空预器内部气体的流速,单位为m/s;ρ为密度,单位为kg/m<sup>3</sup>;c为比热,单位为J/kg·K,密度、比热指空预器内部气体的密度和比热;<img file="FA20177972200810231772001C00014.GIF" wi="36" he="43" />为蓄热元件横截面积比,单位为m<sup>2</sup>/m<sup>2</sup>;T为温度,单位为K;ω为旋转角速度,单位为rad/s;σ为蓄热元件体积比,单位为m<sup>3</sup>/m<sup>3</sup>;λ为空预器内部金属蓄热元件的热传导率,单位为W/m·K;r,z,θ分别为径向、轴向和周向坐标变量,单位分别为m、m、rad,α为流体与金属蓄热元件换热面的换热系数,单位为W/m<sup>2</sup>·K;β为蓄热元件单位体积上换热面积,单位为m<sup>2</sup>/m<sup>3</sup>;下标g和m分别表示烟气和蓄热元件;步骤3:对边界温度值进行初始化,设置T<sub>m</sub>、T<sub>g</sub>为任意值,并确定边界条件如下:径向:i=0时,令T<sub>m</sub>(-1,j,k)=T<sub>m</sub>(0,j,k);i=N<sub>r</sub>时,考虑热损失系数α′,令T<sub>m</sub>(N<sub>r</sub>+1,j,k)=T<sub>m</sub>(N<sub>r</sub>,j,k)-Δrα′/λ<sub>g</sub>[T<sub>m</sub>(N<sub>r</sub>,j,k)-T′]周向:j=0时,令<img file="FA20177972200810231772001C00021.GIF" wi="744" he="71" />+号定义为上一个测算边的结果;j=N<sub>θ</sub>时,令<img file="FA20177972200810231772001C00022.GIF" wi="812" he="71" />+号定义为上一个测算边的结果;轴向:k=0时,令T<sub>m</sub>(i,j,-1)=T<sub>m</sub>(i,j,0);k=N<sub>z</sub>时,令T<sub>m</sub>(i,j,N<sub>z</sub>+1)=T<sub>m</sub>(i,j,N<sub>z</sub>);步骤4:将步骤2所得的方程式(1)和方程式(2)合并,并写成差分方程式(3)的形式:αβT<sub>m</sub>(i,j,k)-αβT<sub>g</sub>(i,j,k)+ωσ<sub>m</sub>ρ<sub>m</sub>c<sub>m</sub>×Δ<sub>1</sub>T<sub>m</sub>+Δ<sub>2</sub>T<sub>m</sub>+Δ<sub>3</sub>T<sub>m</sub>+Δ<sub>4</sub>T<sub>m</sub>+Δ<sub>5</sub>T<sub>m</sub>=0(3)其中,<img file="FA20177972200810231772001C00023.GIF" wi="943" he="146" /><img file="FA20177972200810231772001C00031.GIF" wi="1266" he="134" /><img file="FA20177972200810231772001C00032.GIF" wi="1326" he="132" /><img file="FA20177972200810231772001C00033.GIF" wi="971" he="137" /><img file="FA20177972200810231772001C00034.GIF" wi="1226" he="128" />方程式(3)的抛物线型偏差分方程用标准的数值方法进行测算,柱面坐标r,θ,z用变量i,j,k表示,分别代表其径向、周向和轴向的位置,i,j,k的取值范围分别为:i∈[0,N<sub>r</sub>],j∈[0,N<sub>θ</sub>],k∈[0,N<sub>z</sub>]    (4)其中N表示在相应坐标分割点的数目,通过迭代计算,可得到该段(i,j,k)流体和受热面的温度;步骤5:对方程式(3)进行迭代,若迭代过程前后两次得到的温度误差符合设定的温度参考值的要求,则计算结束,将所得温度场数值输入到空预器系统主控制器中,用计算的温度值取代实际检测温度值,完成空预器热点的检测与判别;否则更新步骤3中的边界温度值,回到步骤4继续。
地址 710048 陕西省西安市金花南路5号