MATLAB在排汽焓线性拟合中的应用研究.pdfVIP

MATLAB在排汽焓线性拟合中的应用 全炳丈 (中国国电天津第一热电厂天津300171) 【摘要1本文利用汽轮机变工况热力特性数据，通过利用MATLAB软件的曲线拟台功能对 汽轮机低压缸排汽蛤进行线性回归，得到排汽压力一排汽培线性曲线及其对应函数关系．最 后，利用某热电厂200MW机组的汽轮机热力特性数据给出排汽焓线性回归的实例并进行误差 分析．通过样本误差分析得到排汽压力一排汽焙的合理线性关系，用该方法得到的排汽焓可 满足在线计算． 【关键词】排汽培MATLAB拟合误差分析 l引言 目前，大型火电机组在线能损分析计算中，低压缸的效率计算及在应用等效热降计算机 组变工况一直都是很麻烦。当最末几级抽汽处于湿蒸汽区时，其压力和温度不再是独立参数， 其焓值不能根据压力、温度来确定，原因是尚不具备在线测量湿蒸汽干度的手段。因此，机 组在线热平衡计算或等效热降计算中，其焓值不能由焓熵图或水蒸气表查得，从而低压缸排 汽焓的计算不能采用常规方法。另外，机组在线经济性分析计算要求较高的实时性，这就要 求整个计算周期要小于数据采集周期，以避免出现计算的滞后。 虽然汽轮机排汽点在湿蒸汽区，其焓值不能根据压力、温度来确定。但可以根据汽轮机 变工况运行中有关数据，找出排汽焓与某一个或几个运行参数之间的相互关系，通过这种方 法将排汽点的焓值计算出来。本文中提出了排汽焓线性拟合计算方法。其基本思路是根据汽 轮机变工况数据利用机组排汽焓与排汽压力和湿蒸汽干度之间的内在关系，进而通过 MATLAB软件的曲线拟合功能求解该机组排汽焓的一般表达式。 2 目前理论上计算汽轮机排汽焓的主要方法 (1)能量平衡法“1由ASMEPTC6A1982年提出，将汽轮机及回热系统看作一个封闭的热 力单元，利用热平衡方程、物质平衡方程和汽轮机功率方程，根据单元能量守恒原则，计算 排汽焓。 (2)曲线外推方法 曲线外推法“1也由ASMEPTC6A--1982年提出，根据汽轮机在过热蒸汽区入口蒸汽状态 点和抽汽状态点做热力过程线，并平滑外推到湿蒸汽区，由此确定排汽焓。 (3)余速损失法 余速损失法。“”1，认为机组背压变化，将主要影响末级相对内效率，而这一末级相对 内效率的变化主要是由于末级余速损失的变化引起的。 (4)曲线迭代法 曲线迭代法”…‘”类似于曲线外推法，以机组热平衡方程为收敛条件，迭代计算。 文献[5]中借助常规计算法的基本原理，在加热器单元的基础上，以反应热平衡方程式 的拓扑矩阵为基本工具，并与文献[6]中采用的膨胀线外推法相结合，进行机组排汽焓的在 线计算。但该方法一方面涉及计算参数较多，另一方面要进行多次迭代且对初值选择的要求 较高。 上述方法的优点是：理论上可以精确计算汽轮机排汽焓；缺点是：很难全面考虑单元的 能量进出，使用极为不便。因此，目前很少使用。方法2的优点是：简单、易于计算机在线 计算；缺点是：精度不够，尤其是在汽轮机低负荷时。方法3需要事先知道汽轮机在额定参 数时末级级前压力，拟合“末级前后压比—排汽余速损失”曲线。方法4精度与所选取迭代 初值密切相关，同时也会影响热经济性在线计算的时间周期。 3排汽焓在线计算方法的提出 由烩熵图(如图l所示)可知对于湿蒸汽，其焓值可由压力和干度唯一确定。即排汽焓 (h)是排汽压力(P)和湿蒸汽干度(x)的函数，用数学方法可表示为： h=f(e，∞ (1) h S 图1 由于目前尚不具备湿蒸汽干度在线测量的手段，因此要在线通过排汽压力和湿蒸汽于 度确定排汽焓很困难。从焓熵图可看出在等熵条件下湿蒸汽的干度随着压力而变化，如果能 找出压力与湿蒸汽干度之间的关系则可把排汽焓看作是排汽压力的复合函数，用数学方法可 表示为： 103 (2) h=厂(P，j(P))