热电机组灵活性提升——储热罐建模与特性分析.docVIP

热电机组灵活性提升——储热罐建模与特性分析 杨大锚 袁福祥 华能国际电力开发公司铜川照金电厂 X 关注成功！ 加关注后您将方便地在 我的关注中得到本文献的被引频次变化的通知！ 新浪微博 腾讯微博 人人网 开心网 豆瓣网 网易微博 摘????要： 热电联产机组“以热定电”运行所导致的调峰能力不足时造成我国三北地区大量弃风的一个重要原因。基于储热罐的蓄能技术是实现热电联产机组热电解耦的一种有效途径。分析了温度分层型储热罐内的换热机理, 建立了储热罐内温度场的一维隐式差分数学模型。模型包含纯导热和纯对流数学模型, 在求解纯导热数学模型的过程中, 引入混合因子来表征水槽内冷温水体之间的混合因素, 通过仿真研究储热罐中温度场的动态特性, 为储热罐的实际运行提供指导依据。 关键词： 储热罐; 热电机组; 非线性模型; 收稿日期：2017.7.4 Received： 2017.7.4 目前, 国外关于配置储热消纳风电的研究主要集中在小型分布式背压燃气机组方面, 包括:峰谷电价下的运营策略, 机组容量与储热罐容量配置计算等。对于储热罐内部温度环境的计算则相对较少;在国内储热技术已经应用于空调系统等领域, 如文献[5]中提及的水蓄冷槽, 其蓄能原理与储热罐类似, 所使用的计算方法可以为储热罐温度环境计算提供指导。本文分析了温度分层型水蓄冷槽内换热机理, 在不同的流量和蓄冷温差下, 对储热罐进行温度场进行激励分析的基础上建立了一维隐式差分数学模型。 1 储热罐温度场建模 储热罐是有温度分层型、迷宫型、空实槽多槽切换型、隔膜型等类型。 1.1 储热罐传热机理建模 以直径D=26m、高h=41.5m的储热罐为例, 数值求解储热罐内温度场布置情况。储热开始前, 罐内充满38℃, 储热开始后, 91℃热水从顶部流入, 38℃冷水从底部流出。将储热罐在高度方向上离散成1000等份, 形成1000个单元, 每个单元厚度△X=0.0415m。每个单元内的水具有相同的温度 (忽略半径方向的温度变化) 。罐内中间温度为64.5℃, 按64℃水的物理特性:密度ρ=984kg/m, 比热C=4.192k J/kg℃, 导热系数0.65w/m℃。 以罐内初始环境全部为冷水为例, 斜温层构造过程应尽量构造充入热水的水平移动, 避免高度方向上的扰动和混合, 这要求布水器出口流型近似重力流, 充入热水的流速小而均匀。设计罐内经历24小时斜温层到达罐中部, 则水流注满一个单元需时172.8s, 将此作为计算时间步长, 在迭代计算中, 每经历一个步长便用i号单元温度取代i+1号单元温度。忽略储热罐直径方向上的温度变化, 忽略储热罐对外散热, 罐内温度场为高度方向的一维温度场, 单元序号自下向上编制。分析i号单元能量平衡关系如图1所示。 图1 储热罐单元能量平衡关系 ??下载原图 即单元中共质的内能变化等于上下端传热功率之差, 如式 (1) 所示: 式 (1) 可进一步表示为: 式中:ti和ti分别为单元i在n时刻和n+1时刻的温度, ti+1和ti-1分别为单元i+1和单元i-1在n+1时刻的温度, S为储热罐横截面积, V为每个单元的体积, V=S*△X, ρ为水密度, C为比热容。整理得: 上述大型稀疏矩阵方程组可利用高斯-赛德尔迭代法进行求解, 对构造斜温层过程进行迭代计算结果如图2所示。 图2 储热罐温度场迭代计算结果 ??下载原图 1.2 充/放热混合因子计算 以上是将储热罐看做纯导热模型进行离散后的代数求解, 求解中没有考虑混合对斜温层的影响, 而在实际过程中, 混合是不可避免的。文献[8]提出混合因子EDDY的概念, 其大小表征了混合的强度。引入混合因子EDDY= (a+3) /a, 其中3取决于混合的强度, 其物理意义是放大了导温系数, 即把流体混合引起的热交换等效于导热。计算储热罐构造斜温层过程相当于将流动看做层流无混合情况, 即3=0, EDDY=1。文献[9]对混合因子分布规律进行了大量的实验验证, 结论是混合因子与入口水流速有关, 且在入口处最大, 从入口到出口呈指数递减规律, 如式 (4) 所示: 式中EDDYi为第i层的混合因子, EDDYinlet为入口处混合因子, 储热过程为储热罐顶端, 放热过程为储热罐底端。 流量的大小可以用雷诺数Re表征;不同蓄热温差表示储热罐内冷热水密度差异, 其可用理查森数Ri表征, 其计算公式分别为: 式中D为特征长度, 取储热罐直径, V为水的运动粘度系数, g为重力加速度, 9.8m/s, β为水的膨胀体积系数, △T为储热温差。分析式 (4) 、 (5) 发现, EDDY受流量变化影响更大, 即EDDY对Re数变化更敏感。储热罐被设计为储热时间18小时, 放热时间6小时, 冷热端温差为定值