CN103267539B - 一种卧式三段式给水加热器上端差及下端差的测算方法 - Google Patents
一种卧式三段式给水加热器上端差及下端差的测算方法 Download PDFInfo
- Publication number
- CN103267539B CN103267539B CN201310142718.XA CN201310142718A CN103267539B CN 103267539 B CN103267539 B CN 103267539B CN 201310142718 A CN201310142718 A CN 201310142718A CN 103267539 B CN103267539 B CN 103267539B
- Authority
- CN
- China
- Prior art keywords
- ntu
- temperature
- cooling section
- well heater
- feedwater
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001816 cooling Methods 0.000 claims abstract description 102
- 238000012546 transfer Methods 0.000 claims abstract description 54
- 238000009833 condensation Methods 0.000 claims abstract description 47
- 230000005494 condensation Effects 0.000 claims abstract description 47
- 230000002209 hydrophobic effect Effects 0.000 claims description 84
- 238000000605 extraction Methods 0.000 claims description 46
- 238000005259 measurement Methods 0.000 claims description 14
- 210000003296 saliva Anatomy 0.000 claims description 12
- 229920006395 saturated elastomer Polymers 0.000 claims description 12
- 238000012886 linear function Methods 0.000 abstract 3
- 238000012360 testing method Methods 0.000 description 26
- 238000004364 calculation method Methods 0.000 description 5
- 238000009530 blood pressure measurement Methods 0.000 description 3
- 238000010977 unit operation Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004781 supercooling Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种卧式三段式给水加热器上端差及下端差的测算方法,分别对三段式加热器的疏水冷却段、凝结段及过热蒸汽冷却段进行数据拟合得到各段传热单元数与某一参数的简化线性函数关系,据此计算得到上、下端差;本发明基于量纲分析原理和运行可测参数,分别对加热器疏水冷却段、凝结段及过热蒸汽冷却段进行分析,利用三段式加热器各段传热单元数与某一参数的简化线性函数关系,最终可求得任意给定运行工况下加热器进出口各热力参数,得到上、下端差,为确定加热器的变工况性能提供参考,测算方法简捷、应用范围广、精度高。
Description
技术领域
本发明涉及一种卧式三段式给水加热器,具体涉及一种卧式三段式给水加热器上端差及下端差的测算方法。
背景技术
现代大型火电厂均采用回热系统减少冷源损失,提高机组循环热效率,一般可使热经济性提高10%~20%。其中,给水加热器是回热系统的核心,其运行性能的好坏对整个机组产生极大的影响,现场通常采用端差的大小来评价加热器运行特性。
目前大型火电机组高压加热器大都采用卧式三段式U形管加热器,给水从给水进口管进入水室,通过U形管下端进入传热区段,首先进入疏水冷却段,再进入蒸汽凝结段,最后是过热蒸汽冷却段,进入水室上部后完成该级的加热工作。蒸汽首先进入过热蒸汽冷却段,将一部分过热度传给给水,蒸汽流出过热蒸汽冷却段后进入加热器的主要传热区蒸汽凝结段,将蒸汽的潜热传给给水,蒸汽凝结后成为凝结水进入疏水冷却段,通过疏水冷却段使疏水具有一定的过冷度。
在机组实际运行过程中,虽然设有加热器的出水和疏水温度测点,可以获得其变工况特性,但因运行条件恶劣、检修维护薄弱等原因,普遍存在测量可靠性和经济性较差的问题:热工测量系统中常采用电阻式传感器,与之相应的数据采集系统测量成本较高;工况变化较大时,水温响应表现出较大的热惯性,影响测量精度;现场安装位置复杂,一旦传感器出现故障会导致测量数据错误或缺失,因此需要一种适用于现场快速测算端差的方法。
传统的加热器端差计算方法中,有的在机组负荷、加热器入口水温等条件变化时,假设加热器传热系数不变,虽然计算过程简化,但与实际情况不符,有一定误差;有的用传热学原理及热平衡理论,考虑到加热器入口水温、汽轮机负荷等对传热系数的影响,提出低压上、下端差应达值的计算方法,但没有考虑疏水冷却段的过热蒸汽冷却段,只适用于低压加热器;还有的方法在计算端差时将加热器各段分界点处汽侧参数进行假定,影响到计算的精度;其他的计算方法亦或局限于两段式低压加热器,或局限于上端差的求解。
发明内容
发明目的:本发明的目的在于为了克服现有技术的不足,提供一种测量方法简单、测量值可靠精确且成本低的卧式三段式给水加热器上端差及下端差的测算方法。
技术方案:本发明所述的一种卧式三段式给水加热器上端差及下端差的测算方法,包括以下步骤:
(1)计算加热器疏水冷却段的传热单元数NTU1与抽汽流量Gs、给水流量Gw的函数关系;
(2)计算加热器凝结段的传热单元数NTU2与给水流量Gw、饱和压力ph的函数关系;
(3)求任意运行工况下加热器的上端差及下端差:
(31)在某一给定运行工况下,测量加热器的饱和压力ph和给水进口温度tw1,由此计算得到加热器的饱和温度th和饱和汽焓hh;同时测量加热器的给水流量Gw、给水压力pw、疏水压力pod、抽汽压力ps和抽汽温度ts,根据给水压力pw和给水进口温度tw1计算得到疏水冷却段的给水定压比热Cpl1,根据抽汽压力ps和抽汽温度ts计算得到加热器的抽汽焓值hs;
步骤(32)~(37)为迭代Ⅱ:
(32)设定本次迭代Ⅱ的抽汽流量初始值Gs k2,其中,k2为迭代Ⅱ次数,未开始迭代Ⅱ时k2=0,Gs 0为步骤(31)测量得到的给水流量Gw,k2大于0时,Gs k2为上一次迭代Ⅱ求得的加热器抽汽流量Gs k2′,由Gs k2和步骤(31)中测得的给水流量Gw根据步骤(1)得到的函数关系,求得本次迭代Ⅱ中加热器的疏水冷却段的传热单元数NTU1 k2;
(33)将步骤(31)中得到的饱和温度th、给水进口温度tw1和步骤(32)求得的疏水冷却段的传热单元数NTU1 k2带入公式tod=th-[1-exp(-NTU1)](th-tw1),得到本次迭代Ⅱ中加热器的疏水温度tod k2,根据本次迭代Ⅱ中加热器的疏水温度tod k2和步骤(31)测得的疏水压力pod计算本次迭代Ⅱ中疏水冷却段的疏水定压比热C'pl1 k2;
步骤(34)~(35)为迭代Ⅲ:
(34)设定疏水冷却段给水出口温度t1t k3,其中,k3为迭代Ⅲ次数,未开始迭代Ⅲ时k3=0,t1t 0为步骤(31)测算得到的饱和温度th,k3大于0时,t1t k3=th-0.001k3,将步骤(31)中测得的饱和温度th、给水进口温度tw1、步骤(33)中计算得到的本次迭代Ⅱ中加热器的疏水温度tod k2以及疏水冷却段给水出口温度t1t k3代入公式计算得到本次迭代Ⅲ中的疏水冷却段传热单元数NTU1 k3;
(35)比较步骤(32)中得到的NTU1 k2和步骤(34)中得到的NTU1 k3,判断是否满足|NTU1 k3-NTU1 k2|≤r3,其中,r3为工程应用的精度要求,0≤r3≤10-3,满足则进入步骤(36),否则令k3=k3+1并回到步骤(34);
(36)将步骤(31)得到的给水流量Gw、饱和温度th、给水进口温度tw1、给水定压比热Cpl1、步骤(33)得到的疏水温度tod k2、疏水定压比热C'pl1 k2和步骤(34)得到的疏水冷却段给水出口温度t1t k3带入疏水冷却段能量平衡方程求得加热器抽汽流量Gs k2′;
(37)比较步骤(36)中求得的加热器抽汽流量Gs k2′与步骤(32)中的本次迭代Ⅱ的抽汽流量初始值Gs k2,判断是否满足|Gs k2’-Gs k2|≤r2,其中,r2为工程应用的精度要求,0≤r2≤10-3,满足则进入步骤(38),否则令k2=k2+1并回到步骤(32);
最终根据上述迭代Ⅱ和迭代Ⅲ步骤,计算得到某一给定运行工况下的加热器抽汽流量Gs、疏水冷却段传热单元数NTU1、加热器疏水温度tod以及疏水冷却段给水出口温度t1t;
(38)根据步骤(31)测得的给水流量Gw和饱和压力ph,通过步骤(2)得到的函数关系,计算得到该给定运行工况下的加热器凝结段的传热单元数NTU2;
(39)将步骤(31)得到的饱和温度th、步骤(34)中的疏水冷却段给水出口温度t1t k3和步骤(38)求得的凝结段的传热单元数NTU2代入公式t2t=th-(th-t1t)exp(-NTU2),求得该给定运行工况下的加热器凝结段给水出口温度t2t;
(310)根据步骤(31)得到的给水压力pw和步骤(39)得到的加热器凝结段给水出口温度t2t,求得过热蒸汽冷却段的定压比热Cpl3;
(311)根据步骤(31)得到的饱和汽焓hh、抽汽焓值hs、给水流量Gw、步骤(36)求得的加热器抽汽流量Gs k2′、步骤(39)求得的加热器凝结段给水出口温度t2t和步骤(310)求得的过热蒸汽冷却段的定压比热Cpl3,通过能量平衡方程求得该给定运行工况下的过热蒸汽冷却段给水出口温度tw2;
(312)将步骤(31)中得到的饱和温度th和步骤(311)中得到的过热蒸汽冷却段给水出口温度tw2代入公式δ=th-tw2,计算得到该给定运行工况下加热器的上端差δ;将步骤(31)中得到的给水进口温度tw1和步骤(33)中得到的疏水温度tod k2代入公式θ=tod-tw1,计算得到该给定运行工况下加热器的下端差θ。
优选的,步骤(1)包括如下步骤:
(11)分别在至少五种不同功率的运行工况下测量得到加热器给水压力pw、饱和压力ph、疏水温度tod、给水进口温度tw1、加热器抽汽流量Gs和给水流量Gw,计算加热器的饱和温度th;
(12)将每一个运行工况下的饱和温度th、疏水温度tod和给水进口温度tw1代入公式计算得到每一个运行工况下加热器疏水冷却段的传热单元数NTU1;
(13)根据步骤(11)和步骤(12)得到的不同运行工况下的加热器抽汽流量Gs、给水流量Gw和疏水冷却段的传热单元数NTU1,利用最小二乘方法对NTU1值进行线性拟合,得到第一常数C1和第二常数C2,从而得到NTU1值与Gs/Gw的函数关系为
更进一步,步骤(2)包括以下步骤:
步骤(21)~(22)为迭代Ⅰ:
(21)设定疏水冷却段出口水温t1t k1,其中,k1为迭代Ⅰ的次数,未开始迭代Ⅰ时k1=0,t1t 0为测量得到的饱和温度th,k1大于0时,t1t k1=th-0.001k1,将步骤(11)中测量得到加热器的饱和温度th、疏水温度tod、给水进口温度tw1,连同疏水冷却段出口水温t1t k1代入公式计算得到本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1;
(22)比较步骤(21)中求得的本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1与步骤(12)中求得的每一个运行工况下加热器疏水冷却段的传热单元数NTU1进行比较,判断是否满足|NTU1 k1-NTU1|≤r1,其中,r1为工程应用的精度要求,0≤r1≤10-3,如满足则进入步骤(23),否则令k1=k1+1并回到步骤(21);
(23)根据步骤(21)求得的满足精度要求的疏水冷却段出口水温t1t k1和步骤(11)测得的给水压力pw,计算得到疏水冷却段出口水焓值h1t;
(24)根据步骤(11)测得的饱和压力ph,计算得到加热器饱和汽焓hh和饱和水焓值h′h,将加热器饱和汽焓hh、饱和水焓值h′h、步骤(11)中测得的加热器抽汽流量Gs、给水流量Gw以及步骤(23)中得到的疏水冷却段出口水焓值h1t带入公式计算得到凝结段出口水焓值h2t;
根据计算得到的凝结段出口水焓值h2t和步骤(11)测量得到的给水压力pw,计算得到凝结段出口水温度t2t;
(25)将步骤(11)中测得的饱和温度th、步骤(21)中的疏水冷却段出口水温t1t k1和步骤(24)中求得的凝结段出口水温度t2t代入公式计算得到加热器凝结段的传热单元数NTU2;
(26)根据步骤(11)和步骤(25)得到的不同运行工况下的给水流量Gw、饱和压力ph和加热器凝结段的传热单元数NTU2值,利用最小二乘方法对NTU2值进行线性拟合,得到第三常数C3和第四常数C4,从而得到NTU2值与Gw/ph的函数关系为
有益效果:1、本发明基于量纲分析原理和运行可测参数,分别对加热器疏水冷却段、凝结段及过热蒸汽冷却段进行分析,利用三段式加热器各段传热单元数与某一参数的简化线性函数关系,通过迭代方法求得任意给定运行工况下加热器进出口各热力参数,准确且误差小,最后得到上、下端差,为确定加热器的变工况性能提供参考,测算方法简捷、应用范围广、精度高;2、只需采用常见简单的测量仪器得到准确的测量数据,简化测量现场,保证测量可靠性,节约成本,而后对加热器拟合测量得到的数据即可得到疏水冷却段及凝结段传热单元数与某一参数的线性关系,据此测算出任意给定运行工况下加热器的上、下端差,而不需要了解结构参数,模型简捷;3、提出的方法按照疏水冷却段、凝结段、过热段的顺序分段依次计算,因此既适用于三段式高压加热器,又适用于不含过热段的两段式低压加热器,整个计算过程的推导公式具有一定的通用性,适应于变工况特性,不受热惯性的影响,方法应用范围广;4、提出的方法既能计算上端差,又能计算下端差,弥补了一些传统计算方法中局限于上端差求解的不足。
附图说明
图1为本发明实施例中的23种试验工况下的汽轮机发电机组电功率变化图;
图2为本发明疏水冷却段传热单元数NTU1与抽汽流量Gs、给水流量Gw的比值Gs/Gw的线性关系图;
图3为本发明凝结段传热单元数NTU2与给水流量Gw、饱和压力ph的比值Gw/ph的线性关系图。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
实施例:本实施例中测算哈汽某330MW汽轮机组#1高压加热器上端差及下端差的方法,包括以下步骤:
(1)计算加热器疏水冷却段的传热单元数NTU1与抽汽流量Gs、给水流量Gw的比值Gs/Gw的线性关系:
(11)设定23种不同功率的试验工况,23种试验工况对应汽轮机发电机组电功率变化如图1所示,分别在每个试验工况下,利用温度压力测量仪表测量得到加热器给水压力pw、饱和压力ph、疏水温度tod和给水进口温度tw1,利用流量测量仪表测量得到加热器抽汽流量Gs和给水流量Gw,根据饱和压力ph通过IFC-67工业用水和水蒸汽热力性质模型计算得到加热器的饱和温度th,各参数见表1;
表1 23种不同功率试验工况的测算值
(12)将每一个试验工况下的饱和温度th、疏水温度tod和给水进口温度tw1代入公式计算得到每一个试验工况下加热器疏水冷却段的传热单元数NTU1,见表2;
表2 23种不同功率试验工况的疏水冷却段NTU1值
(13)根据步骤(11)和步骤(12)得到的不同试验工况下的加热器抽汽流量Gs、给水流量Gw和疏水冷却段的传热单元数NTU1,利用最小二乘方法对NTU1值进行线性拟合,得到第一常数C1=-30.216、第二常数C2=3.172,从而得到NTU1值与Gs/Gw的函数关系为拟合得到的线性关系如图2所示;
(2)计算加热器凝结段的传热单元数NTU2与给水流量Gw、饱和压力ph的比值Gw/ph的线性关系:
步骤(21)~(22)为迭代Ⅰ:
(21)设定疏水冷却段出口水温t1t k1,其中,k1为迭代Ⅰ的次数,未开始迭代Ⅰ时k1=0,t1t 0为测量得到的饱和温度th,k1大于0时,t1t k1=th-0.001k1,将步骤(11)中测量得到加热器的饱和温度th、疏水温度tod、给水进口温度tw1,连同疏水冷却段出口水温t1t k1代入公式计算得到本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1;
(22)比较步骤(21)中求得的本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1与步骤(12)中求得的每一个试验工况下加热器疏水冷却段的传热单元数NTU1进行比较,判断是否满足|NTU1 k1-NTU1|≤r1,其中,r1为工程应用的精度要求,0≤r1≤10-3,如满足则进入步骤(23),否则令k1=k1+1并回到步骤(21),通过迭代Ⅰ计算得到的满足精度要求的疏水冷却段出口水温t1t,见表3;
表3 23种不同功率试验工况的疏水冷却段出口水温t1t
(23)根据步骤(21)求得的满足精度要求的疏水冷却段出口水温t1t k1和步骤(11)测得的给水压力pw,通过IFC-67工业用水和水蒸汽热力性质模型计算得到疏水冷却段出口水焓值h1t;
(24)根据步骤(11)测得的饱和压力ph,通过IFC-67工业用水和水蒸汽热力性质模型计算得到加热器饱和汽焓值hh和饱和水焓值h′h,将加热器饱和汽焓值hh、饱和水焓值h′h、步骤(11)中测得的加热器抽汽流量Gs、给水流量Gw以及步骤(23)中得到的疏水冷却段出口水焓值h1t带入公式计算得到凝结段出口水焓值h2t;
根据计算得到的凝结段出口水焓值h2t和步骤(11)测量得到的给水压力pw,通过IFC-67工业用水和水蒸汽热力性质模型计算得到凝结段出口水温度t2t,23种不同试验工况下的凝结段计算结果见表4;
表4 23种不同功率试验工况的凝结段计算结果
(25)将步骤(11)中测得的饱和温度th、步骤(21)中的疏水冷却段出口水温t1t k1和步骤(24)中求得的凝结段出口水温度t2t代入公式计算得到加热器凝结段的传热单元数NTU2,结果见表5;
表5 23种不同功率试验工况的凝结段段NTU2值
(26)根据步骤(11)和步骤(25)得到的不同试验工况下的给水流量Gw、饱和压力ph和加热器凝结段的传热单元数NTU2值,利用最小二乘方法对NTU2值进行线性拟合,得到第三常数C3=0.000660和第四常数C4=0.778,从而得到NTU2值与Gw/ph的函数关系为,拟合得到的线性关系如图3所示;
(3)求取任意5个试验工况下加热器的上端差及下端差:
(31)在每个给定试验工况下,利用温度压力测量仪表测量加热器的饱和压力ph和给水进口温度tw1,根据饱和压力ph,通过IFC-67工业用水和水蒸汽热力性质模型计算得到加热器的饱和温度th和饱和汽焓hh;
同时,利用流量测量仪表测量加热器的给水流量Gw,利用压力测量仪表测量给水压力pw和疏水压力pod,根据给水压力pw和给水进口温度tw1,通过IFC-67工业用水和水蒸汽热力性质模型计算疏水冷却段给水定压比热Cpl1;
利用温度压力测量仪表测量抽汽压力ps和抽汽温度ts,根据抽汽压力ps和抽汽温度ts,通过IFC-67工业用水和水蒸汽热力性质模型计算得到加热器的抽汽焓值hs;
5种不同试验工况下测算的部分热力参数如表1所示:
表1 5种不同试验工况下测算得到的各热力参数
步骤(32)~(37)为迭代Ⅱ:
(32)设定本次迭代Ⅱ的抽汽流量初始值Gs k2,其中,k2为迭代Ⅱ次数,未开始迭代Ⅱ时k2=0,Gs 0为步骤(31)测量得到的给水流量Gw,k2大于0时,Gs k2为上一次迭代Ⅱ求得的加热器抽汽流量Gs k2′,由Gs k2和步骤(31)中测得的给水流量Gw根据步骤(1)得到的函数关系,求得本次迭代Ⅱ中加热器的疏水冷却段的传热单元数NTU1 k2;
(33)将步骤(31)中得到的饱和温度th、给水进口温度tw1和步骤(32)求得的疏水冷却段的传热单元数NTU1 k2带入公式tod=th-[1-exp(-NTU1)](th-tw1),得到本次迭代Ⅱ中加热器的疏水温度tod k2,根据本次迭代Ⅱ中加热器的疏水温度tod k2和步骤(31)测得的疏水压力pod,通过IFC-67工业用水和水蒸汽热力性质模型计算本次迭代Ⅱ中疏水冷却段的疏水定压比热C'pl1 k2;
步骤(34)~(35)为迭代Ⅲ:
(34)设定疏水冷却段给水出口温度t1t k3,其中,k3为迭代Ⅲ次数,未开始迭代Ⅲ时k3=0,t1t 0为步骤(31)测算得到的饱和温度th,k3大于0时,t1t k3=th-0.001k3,将步骤(31)中测得的饱和温度th、给水进口温度tw1、步骤(33)中计算得到的本次迭代Ⅱ中加热器的疏水温度tod k2以及疏水冷却段给水出口温度t1t k3代入公式计算得到本次迭代Ⅲ中的疏水冷却段传热单元数NTU1 k3;
(35)比较步骤(32)中得到的NTU1 k2和步骤(34)中得到的NTU1 k3,判断是否满足|NTU1 k3-NTU1 k2|≤r3,其中,r3为工程应用的精度要求,0≤r3≤10-3,满足则进入步骤(36),否则令k3=k3+1并回到步骤(34);
(36)将步骤(31)得到的给水流量Gw、饱和温度th、给水进口温度tw1、给水定压比热Cpl1、步骤(33)得到的疏水温度tod k2、疏水定压比热C'pl1 k2和步骤(34)得到的疏水冷却段给水出口温度t1t k3带入疏水冷却段能量平衡方程求得加热器抽汽流量Gs k2′;
(37)比较步骤(36)中求得的加热器抽汽流量Gs k2′与步骤(32)中的本次迭代Ⅱ的抽汽流量初始值Gs k2,判断是否满足|Gs k2’-Gs k2|≤r2,其中,r2为工程应用的精度要求,0≤r2≤10-3,满足则进入步骤(38),否则令k2=k2+1并回到步骤(32);
根据上述迭代Ⅱ和迭代Ⅲ步骤,计算5个试验工况下满足精度要求的加热器抽汽流量Gs、疏水冷却段传热单元数NTU1、加热器疏水温度tod以及疏水冷却段给水出口温度t1t,如表7所示:
表7 5种不同试验工况下通过迭代Ⅱ和迭代Ⅲ计算得到的各热力参数
(38)根据步骤(31)测得的给水流量Gw和饱和压力ph,通过步骤(26)得到的线性关系式计算得到5个试验工况下的加热器凝结段的传热单元数NTU2,分别为2.071、2.154、2.572、3.056和3.861;
(39)将步骤(31)得到的饱和温度th、步骤(34)求得的疏水冷却段给水出口温度t1t k3和步骤(38)求得的NTU2代入公式t2t=th-(th-t1t)exp(-NTU2),求得5个给定试验工况下的加热器凝结段给水出口温度t2t,分别为270.086℃、265.593℃、258.944℃、251.789℃和243.226℃;
(310)根据步骤(31)得到的给水压力pw和步骤(39)得到的加热器凝结段给水出口温度t2t,通过IFC-67工业用水和水蒸汽热力性质模型求得过热蒸汽冷却段的定压比热Cpl3;
(311)根据步骤(31)得到的饱和汽焓hh、抽汽焓值hs、给水流量Gw、步骤(36)求得的加热器抽汽流量Gs k2′、步骤(39)求得的加热器凝结段给水出口温度t2t和步骤(310)求得的过热蒸汽冷却段的定压比热Cpl3,通过能量平衡方程求得5个给定试验工况下的过热蒸汽冷却段给水出口温度tw2,分别为274.550℃、270.182℃、263.070℃、255.481℃、247.500℃;
(312)将步骤(31)中得到的饱和温度th和步骤(311)中得到的过热蒸汽冷却段给水出口温度tw2代入公式δ=th-tw2,计算得到该给定运行工况下加热器的上端差δ;将步骤(31)中得到的给水进口温度tw1和步骤(33)中得到的疏水温度tod k2代入公式θ=tod-tw1,计算得到该给定运行工况下加热器的下端差θ;
将计算得到的加热器机组的上端差及下端差与机组运行前的试验测量数据进行对比,如表8所示:
表8 计算上、下端差与机组试验数据的对比
可以看出,本发明提供的测算方法能够根据机组运行数据检测加热器在运行过程中任意工况下的上端差及下端差,且测算过程简单清晰,所需测点较少,测算所需的原始机组运行数据易于得到且准确,因此应用在不同负荷运行工况下计算得到的上、下端差值都有较高的精度。
如上,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (3)
1.一种卧式三段式给水加热器上端差及下端差的测算方法,其特征在于:包括以下步骤:
(1)计算加热器疏水冷却段的传热单元数NTU1与抽汽流量Gs、给水流量Gw的函数关系
(2)计算加热器凝结段的传热单元数NTU2与给水流量Gw、饱和压力ph的函数关系
(3)求任意运行工况下加热器的上端差及下端差:
(31)在某一给定运行工况下,测量加热器的饱和压力ph和给水进口温度tw1,由此计算得到加热器的饱和温度th和饱和汽焓hh;同时测量加热器的给水流量Gw、给水压力pw、疏水压力pod、抽汽压力ps和抽汽温度ts,根据给水压力pw和给水进口温度tw1计算得到疏水冷却段的给水定压比热Cpl1,根据抽汽压力ps和抽汽温度ts计算得到加热器的抽汽焓值hs;
步骤(32)~(37)为迭代Ⅱ:
(32)设定本次迭代Ⅱ的抽汽流量初始值Gs k2,其中,k2为迭代Ⅱ次数,未开始迭代Ⅱ时k2=0,Gs 0为步骤(31)测量得到的给水流量Gw,k2大于0时,Gs k2为上一次迭代Ⅱ求得的加热器抽汽流量Gs k2′,由抽汽流量初始值Gs k2和步骤(31)中测得的给水流量Gw根据步骤(1)得到的函数关系,求得本次迭代Ⅱ中加热器的疏水冷却段的传热单元数NTU1 k2;
(33)将步骤(31)中得到的饱和温度th、给水进口温度tw1和步骤(32)求得的疏水冷却段的传热单元数NTU1 k2带入公式tod=th-[1-exp(-NTU1)](th-tw1),得到本次迭代Ⅱ中加热器的疏水温度tod k2,根据本次迭代Ⅱ中加热器的疏水温度tod k2和步骤(31)测得的疏水压力pod计算本次迭代Ⅱ中疏水冷却段的疏水定压比热C'pl1 k2;
步骤(34)~(35)为迭代Ⅲ:
(34)设定疏水冷却段给水出口温度t1t k3,其中,k3为迭代Ⅲ次数,未开始迭代Ⅲ时k3=0,t1t 0为步骤(31)测算得到的饱和温度th,k3大于0时,t1t k3=th-0.001k3,将步骤(31)中测得的饱和温度th、给水进口温度tw1、步骤(33)中计算得到的本次迭代Ⅱ中加热器的疏水温度tod k2以及疏水冷却段给水出口温度t1t k3代入公式计算得到本次迭代Ⅲ中的疏水冷却段传热单元数NTU1 k3;
(35)比较步骤(32)中得到的NTU1 k2和步骤(34)中得到的NTU1 k3,判断是否满足|NTU1 k3-NTU1 k2|≤r3,其中,r3为工程应用的精度要求,0≤r3≤10-3,满足则进入步骤(36),否则令k3=k3+1并回到步骤(34);
(36)将步骤(31)得到的给水流量Gw、饱和温度th、给水进口温度tw1、给水定压比热Cpl1、步骤(33)得到的疏水温度tod k2、疏水定压比热C'pl1 k2和步骤(34)得到的疏水冷却段给水出口温度t1t k3带入疏水冷却段能量平衡方程求得加热器抽汽流量Gs k2′;
(37)比较步骤(36)中求得的加热器抽汽流量Gs k2′与步骤(32)中的本次迭代Ⅱ的抽汽流量初始值Gs k2,判断是否满足|Gs k2’-Gs k2|≤r2,其中,r2为工程应用的精度要求,0≤r2≤10-3,满足则进入步骤(38),否则令k2=k2+1并回到步骤(32);
(38)根据步骤(31)测得的给水流量Gw和饱和压力ph,通过步骤(2)得到的函数关系,计算得到该给定运行工况下的加热器凝结段的传热单元数NTU2;
(39)将步骤(31)得到的饱和温度th、步骤(34)中的疏水冷却段给水出口温度t1t k3和步骤(38)求得的凝结段的传热单元数NTU2代入公式t2t=th-(th-t1t)exp(-NTU2),求得该给定运行工况下的加热器凝结段给水出口温度t2t;
(310)根据步骤(31)得到的给水压力pw和步骤(39)得到的加热器凝结段给水出口温度t2t,求得过热蒸汽冷却段的定压比热Cpl3;
(311)根据步骤(31)得到的饱和汽焓hh、抽汽焓值hs、给水流量Gw、步骤(36)求得的加热器抽汽流量Gs k2′、步骤(39)求得的加热器凝结段给水出口温度t2t和步骤(310)求得的过热蒸汽冷却段的定压比热Cpl3,通过能量平衡方程求得该给定运行工况下的过热蒸汽冷却段给水出口温度tw2;
(312)将步骤(31)中得到的饱和温度th和步骤(311)中得到的过热蒸汽冷却段给水出口温度tw2代入公式δ=th-tw2,计算得到该给定运行工况下加热器的上端差δ;将步骤(31)中得到的给水进口温度tw1和步骤(33)中得到的疏水温度tod k2代入公式θ=tod-tw1,计算得到该给定运行工况下加热器的下端差θ。
2.根据权利要求1所述的卧式三段式给水加热器上端差及下端差的测算方法,其特征在于:步骤(1)包括如下步骤:
(11)分别在至少五种不同功率的运行工况下测量得到加热器给水压力pw、饱和压力ph、疏水温度tod、给水进口温度tw1、加热器抽汽流量Gs和给水流量Gw,计算加热器的饱和温度th;
(12)将每一个运行工况下的饱和温度th、疏水温度tod和给水进口温度tw1代入公式计算得到每一个运行工况下加热器疏水冷却段的传热单元数NTU1;
(13)根据步骤(11)和步骤(12)得到的不同运行工况下的加热器抽汽流量Gs、给水流量Gw和疏水冷却段的传热单元数NTU1,利用最小二乘方法对NTU1值进行线性拟合,得到第一常数C1和第二常数C2,从而得到NTU1值与Gs/Gw的函数关系为
3.根据权利要求2所述的卧式三段式给水加热器上端差及下端差的测算方法,其特征在于:步骤(2)包括以下步骤:
步骤(21)~(22)为迭代Ⅰ:
(21)设定疏水冷却段出口水温t1t k1,其中,k1为迭代Ⅰ的次数,未开始迭代Ⅰ时k1=0,t1t 0为测量得到的饱和温度th,k1大于0时,t1t k1=th-0.001k1,将步骤(11)中测量得到加热器的饱和温度th、疏水温度tod、给水进口温度tw1,连同疏水冷却段出口水温t1t k1代入公式计算得到本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1;
(22)比较步骤(21)中求得的本次迭代Ⅰ中的疏水冷却段传热单元数NTU1 k1与步骤(12)中求得的每一个运行工况下加热器疏水冷却段的传热单元数NTU1进行比较,判断是否满足|NTU1 k1-NTU1|≤r1,其中,r1为工程应用的精度要求,0≤r1≤10-3,如满足则进入步骤(23),否则令k1=k1+1并回到步骤(21);
(23)根据步骤(21)求得的满足精度要求的疏水冷却段出口水温t1t k1和步骤(11)测得的给水压力pw,计算得到疏水冷却段出口水焓值h1t;
(24)根据步骤(11)测得的饱和压力ph,计算得到加热器饱和汽焓hh和饱和水焓值h′h,将加热器饱和汽焓hh、饱和水焓值h′h、步骤(11)中测得的加热器抽汽流量Gs、给水流量Gw以及步骤(23)中得到的疏水冷却段出口水焓值h1t带入公式计算得到凝结段出口水焓值h2t;
根据计算得到的凝结段出口水焓值h2t和步骤(11)测量得到的给水压力pw,计算得到凝结段出口水温度t2t;
(25)将步骤(11)中测得的饱和温度th、步骤(21)中的疏水冷却段出口水温t1t k1和步骤(24)中求得的凝结段出口水温度t2t代入公式计算得到加热器凝结段的传热单元数NTU2;
(26)根据步骤(11)和步骤(25)得到的不同运行工况下的给水流量Gw、饱和压力ph和加热器凝结段的传热单元数NTU2值,利用最小二乘方法对NTU2值进行线性拟合,得到第三常数C3和第四常数C4,从而得到NTU2值与Gw/ph的函数关系为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310142718.XA CN103267539B (zh) | 2013-04-23 | 2013-04-23 | 一种卧式三段式给水加热器上端差及下端差的测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310142718.XA CN103267539B (zh) | 2013-04-23 | 2013-04-23 | 一种卧式三段式给水加热器上端差及下端差的测算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103267539A CN103267539A (zh) | 2013-08-28 |
CN103267539B true CN103267539B (zh) | 2015-06-17 |
Family
ID=49011180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310142718.XA Expired - Fee Related CN103267539B (zh) | 2013-04-23 | 2013-04-23 | 一种卧式三段式给水加热器上端差及下端差的测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103267539B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020189A (zh) * | 2014-05-23 | 2014-09-03 | 国家电网公司 | 一种确定电厂汽轮机回热系统加热器端差的监测方法 |
CN106248133B (zh) * | 2016-08-15 | 2019-05-24 | 上海交通大学 | 一种加热器全工况上端差和下端差应达值的在线估计方法 |
CN106295203B (zh) * | 2016-08-15 | 2020-01-07 | 上海交通大学 | 基于上端差应达值实时计算的机组热经济性在线评估方法 |
CN112799447B (zh) * | 2021-01-06 | 2022-03-04 | 国家电投集团河南电力有限公司 | 一种高压加热器的控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827429A (en) * | 1987-06-16 | 1989-05-02 | Westinghouse Electric Corp. | Turbine impulse chamber temperature determination method and apparatus |
CN101825502A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 |
CN101832545A (zh) * | 2010-04-16 | 2010-09-15 | 东南大学 | 汽机带蒸汽冷却器的加热器出水及疏水温度测算方法 |
-
2013
- 2013-04-23 CN CN201310142718.XA patent/CN103267539B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827429A (en) * | 1987-06-16 | 1989-05-02 | Westinghouse Electric Corp. | Turbine impulse chamber temperature determination method and apparatus |
CN101825502A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 |
CN101832545A (zh) * | 2010-04-16 | 2010-09-15 | 东南大学 | 汽机带蒸汽冷却器的加热器出水及疏水温度测算方法 |
Non-Patent Citations (3)
Title |
---|
杨涛,胥建群,周克毅等.电站给水加热器水位的优化运行.《东南大学学报(自然科学版)》.2012,第42卷第314-319页. * |
火电机组高压给水加热器动态过程的数值分析;张欣刚,徐治皋,李勇等;《中国动力工程学报》;20050415;第25卷(第2期);第262-268、288页 * |
电站给水加热器水位优化运行研究;杨涛,胥建群,周克毅等;《2012电站自动化信息化学术和技术交流会议论文集 》;20121113;第231-237页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103267539A (zh) | 2013-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101825502B (zh) | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 | |
CN104048842A (zh) | 一种基于软测量技术的汽轮机热耗率在线监测方法 | |
CN102967464B (zh) | 凝汽式汽轮机高背压改造后的性能评价方法 | |
CN107201921B (zh) | 一种汽轮机热耗率在线监测系统及测量方法 | |
CN103267539B (zh) | 一种卧式三段式给水加热器上端差及下端差的测算方法 | |
CN103048020B (zh) | 基于性能试验数据的发电厂主蒸汽流量在线计算方法 | |
CN105956370B (zh) | 一种用于燃煤电厂实时计算煤耗等经济指标的方法 | |
CN103487272A (zh) | 直接空冷机组空冷凝汽器进汽焓的计算方法 | |
CN102338568A (zh) | 基于清洁系数指标的电厂凝汽器性能在线监测系统及方法 | |
CN102749156B (zh) | 一种汽轮机排汽焓的检测方法 | |
CN107543733A (zh) | 一种凝汽器热负荷在线监测系统及测量方法 | |
CN104615857A (zh) | 一种凝汽式汽轮机凝汽器热负荷的确定方法 | |
CN103134834A (zh) | 一种湿蒸汽干度测量装置及方法 | |
CN108446465B (zh) | 通过工质分解在线测算火电机组厂用蒸汽量的方法 | |
CN106018730A (zh) | 基于磨煤机入口一次风修正的煤质水分测量装置和方法 | |
CN103063436A (zh) | 汽轮机热耗率指标计算结果判断方法 | |
CN106097151A (zh) | 一种基于数据协调的降低电厂数据不确定度的方法 | |
CN103759769A (zh) | 燃煤电站再热蒸汽流量软测量方法和装置 | |
CN104655673A (zh) | 一种三分仓回转式空气预热器换热性能在线诊断方法 | |
CN102589824B (zh) | 一种凝汽器漏空气检测方法 | |
CN102261936B (zh) | 高加危急疏水泄漏流量确定方法 | |
CN101825503B (zh) | 汽机表面式加热器的出水及疏水温度测算方法 | |
CN213980887U (zh) | 一种确定核电汽轮机抽汽参数的装置 | |
CN105445028A (zh) | 汽轮机热耗率的测定及计算方法 | |
Li et al. | Study on the effect of regenerative system on power type relative internal efficiency of nuclear steam turbine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20171110 Address after: Jing'an District of Nanjing city streets MA Qixia 210059 Jiangsu Province Du Village Co-patentee after: SOUTHEAST University Patentee after: Datang Nanjing power plant Address before: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210018 Patentee before: Southeast University |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 |
|
CF01 | Termination of patent right due to non-payment of annual fee |