CN111047463B - 供热机组采用低温余热热泵供热改造后的调峰能力评估方法 - Google Patents
供热机组采用低温余热热泵供热改造后的调峰能力评估方法 Download PDFInfo
- Publication number
- CN111047463B CN111047463B CN201911212619.8A CN201911212619A CN111047463B CN 111047463 B CN111047463 B CN 111047463B CN 201911212619 A CN201911212619 A CN 201911212619A CN 111047463 B CN111047463 B CN 111047463B
- Authority
- CN
- China
- Prior art keywords
- heat supply
- unit
- steam
- flow
- heat pump
- 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.)
- Active
Links
- 239000002918 waste heat Substances 0.000 title claims abstract description 50
- 238000012986 modification Methods 0.000 title claims abstract description 18
- 230000004048 modification Effects 0.000 title claims abstract description 18
- 238000011156 evaluation Methods 0.000 title claims abstract description 17
- 230000001105 regulatory effect Effects 0.000 title claims description 7
- 238000000605 extraction Methods 0.000 claims abstract description 80
- 238000010438 heat treatment Methods 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 54
- 238000001595 flow curve Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 15
- 238000013461 design Methods 0.000 claims description 14
- 238000009833 condensation Methods 0.000 claims description 13
- 230000005494 condensation Effects 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 10
- 238000010521 absorption reaction Methods 0.000 claims description 8
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims 1
- 230000002194 synthesizing effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
本发明涉及一种供热机组采用低温余热热泵供热改造后的调峰能力评估方法,综合供热机组实际运行低压缸最小排汽量及总供热抽汽流量,通过实际运行低压缸最小排汽量下的特性曲线及各供热抽汽流量下的特性曲线确定出采用低温余热热泵供热改造后机组在供热状态下的调峰最小出力及调峰最大出力。本发明的方法简便且准确可靠,可为电网调度部门确定现有采用低温余热热泵供热改造后供热机组在采暖期不同阶段的运行方式提供重要参考数据。
Description
技术领域
本发明属于发电技术领域,具体涉及一种供热机组采用低温余热热泵供热改造后的调峰能力评估方法。
背景技术
新能源机组的迅猛发展对火电机组深度调峰提出了越来越高的要求。火电机组中,供热机组由于同时对外供应热及电两种产品,其热电耦合特性大大增加了机组参与电网调峰的难度。目前,许多供热机组供热季节调峰出力只能降低至70%左右。
为提高现有火电机组的调峰幅度,国家能源局开展了灵活性改造示范项目,并提出了灵活性改造的调峰性能提升目标。要求示范项目通过灵活性改造,热电机组增加20%额定容量的调峰能力,最小技术出力达到40%-50%额定容量,这对于现有的供热机组确定了较为苛刻的改造目标。
现有供热机组汽轮机进行低温余热热泵供热技术改造后,机组的调峰能力不仅与对外供热量相关,而且与热泵机组的工作性能紧密相关,从电网运行的角度迫切需要确定其确切的调峰能力,以保证对采用低温余热热泵供热的供热机组进行准确的调度。但目前缺乏一种实用的简化的评估技术手段来对现有供热机组采用低温余热热泵供热改造后的调峰能力进行评估。
因此迫切需要供热机组采用低温余热热泵供热改造后的调峰能力评估方法来解决目前面临的问题。
发明内容
本发明的目的是提供一种简便且准确可靠的供热机组采用低温余热热泵供热改造后的调峰能力评估方法。
本发明采用如下技术方案:
一种供热机组采用低温余热热泵供热改造后的调峰能力评估方法,综合供热机组实际运行低压缸最小排汽量及总供热抽汽流量,通过实际运行低压缸最小排汽量下的特性曲线及各供热抽汽流量下的特性曲线确定出采用低温余热热泵供热改造后机组在供热状态下的调峰最小出力及调峰最大出力。
进一步的,其具体包括如下步骤:
(A)确定供热机组热电特性曲线中的考虑适当运行安全裕度的最小等排汽流量曲线;
(B)根据采用低温余热热泵供热改造后供热机组的供热运行数据计算供热机组的实际供热抽汽流量;
(C)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最小出力Ptfmin;
(D)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最大出力Ptfmax;
(E)采用低温余热热泵供热改造后供热机组的当前调峰出力范围为:Ptfmin至Ptfmax。
进一步的,所述步骤(A)具体包括:
(a)在供热机组热电特性曲线中,将低压缸最小排汽流量限制线进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
y1=a1+b1×x (1)
式(1)中,a1、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y1表示机组功率应变量,单位MW;
(b)在供热机组汽轮机设计资料中,查找额定供热工况的以下参数:设计供热抽汽量Fcqdes、设计主蒸汽进汽量Fmsdes、设计低压缸排汽量FLPexdes、以及此工况下供热机组出力Pdes;
(c)确定机组在低压缸设计排汽量FLPexdes下的功率随主蒸汽流量的变化特性曲线;
由等式Pdes=a2+b1×Fmsdes;
可计算得到:a2= Pdes-b1×Fmsdes;
由此在低压缸设计排汽量下,机组功率随主蒸汽流量的变化特性拟合为:
y2=a2+b1×x (2)
式(2)中,a2、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y2表示机组功率应变量,单位MW;
(d)考虑适当运行安全裕度的实际运行低压缸最小等排汽流量曲线可表示为:
y=y1+(y2-y1)/(FLPexdes- FLPexmin)×(FLPex- FLPexmin)
式中,y1与y2分别根据式(1)及式(2)求得。
进一步的,考虑实际运行中低压缸排汽流量比设计低压缸最小排汽流量要偏高,并取偏高值为50t/h作为适当运行安全裕度;实际运行低压缸最小等排汽流量曲线可进一步表示为:
y=y1+(y2-y1)/(FLPexdes- FLPexmin)×50 (3)。
进一步的,所述步骤(B)具体包括:
1)收集现场机组供热运行数据:
包括采用低温余热热泵供热改造后机组热泵的热网水流量Frws、热泵热网水进水温度t1、热泵热网水出水温度t2、低温余热热泵热源水进水温度t3、低温余热热泵驱动蒸汽进口压力pds及进口温度tds、热泵凝结水箱温度tdsc以及采用低温余热热泵供热改造后机组汽轮机热网加热器抽汽流量Fcq1;流量单位均为t/h,温度单位均为℃,压力单位为MPa;
2)计算热网水在热泵中的总吸热量:
Q1=Frws×1000×4.186×(t2-t1) (4)
式(4)中,Q1为热网水在热泵中的总吸热量,单位为kJ/h;
3)计算热泵运行中的实际能效系数:
COPr=COPdes×ftrw×ftry×fpds (5)
式(5)中,COPr为实际运行热泵能效系数;COPdes为热泵设计能效系数;ftrw为热泵热网进水温度修正系数,按热泵热网水进水温度t1根据厂家提供的修正曲线查得;ftry为热泵热源水进水温度修正系数,按热泵热源水进水温度t3根据厂家提供的修正曲线查得;fpds为热泵驱动蒸汽压力修正系数,按热泵驱动蒸汽进口压力pds根据厂家提供的修正曲线查得;
4)计算单位质量热泵驱动蒸汽在热泵中的凝结热量:
q2=h1(pds,tds)- h2(tdsc) (6)
式(6)中,q2为每kg热泵驱动蒸汽在热泵中的凝结热量,单位kJ/kg;h1()为水蒸汽特性公式计算函数,得到的驱动蒸汽焓值,单位kJ/kg;h2()为水蒸气特性公式计算函数,得到的饱和水焓值,单位kJ/kg;
5)计算采用低温余热热泵在供热状态下的驱动蒸汽流量Fds:
Fds= Q1/ COPr/q2/1000 (7)
式(7)中,Fds为低温余热热泵在供热状态下的驱动蒸汽流量,单位t/h;
6)计算机组总供热抽汽流量Fcq;
Fcq= Fds+Fcq1 (8)
式(8)中,Fcq为机组总供热抽汽流量,单位t/h。
进一步的,h1()为水蒸汽特性公式计算函数,由蒸汽压力及温度计算蒸汽焓值。
进一步的,h2()为水蒸汽特性公式计算函数,由水温度计算得到的饱和水焓值。
进一步的,所述步骤(C)具体包括:
1)在供热机组热电特性曲线中,将等抽汽流量曲线进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
zi=ci+di×x (9)
式(9)中,i为由机组等抽汽流量的编号,范围由0至n,随i的增大抽汽流量Fcn(i)逐渐增大,i=0时对应供热抽汽流量为零即纯凝运行工况,i=n时对应最大供热抽汽流量的运行工况;ci、di分别为特性系数;x为机组主蒸汽流量,单位t/h;zi为对应汽轮机采暖抽汽流量Fcn(i)下的机组功率,单位MW;
2)在任意供热抽汽量Fcq时机组功率随主蒸汽流量的变化特性可表示为:
当Fcq(i)<Fcq<Fcq(i+1)时,
z=zi+(zi+1-zi)/( Fcq(i+1)- Fcq(i)) ×( Fcq- Fcq(i)) (10)
式(10)中zi及zi+1均根据公式(9)计算得到;
3)确定供热机组的实际运行低压缸最小等排汽流量曲线及汽轮机供热抽汽流量曲线的相交工况点,即为供热机组的最小调峰出力;
进一步的,所述步骤3)中,确定供热机组的最小调峰出力的具体方法包括:
a、假设主蒸汽流量初始值x0为75%额定主蒸汽流量;设定ST=100,j=1;
b、根据式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y, 根据式(10)计算机组供热抽汽量Fcq下机组功率值为z;比较y与 z;
c、如果abs(y- z)<0.01,相交工况点找到,确定的主蒸汽流量x=x0,机组功率Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S1=1,S0=1,主蒸汽流量x1= x0+SP,j=j+1;
如果abs(y- z)≥0.01且y< z时,S1=-1,S0=-1,主蒸汽流量x1= x0-SP,j=j+1;
d、根据新的主蒸汽流量按式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y,根据式(10)计算机组供热抽汽量Fcq下机组功率值为z;比较y与 z;
e、如果abs(y- z)<0.01,相交工况点找到,确定的主蒸汽流量x=x(j-1),机组功率Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S(j)=1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果abs(y- z)≥0.01且y< z时,S(j)=-1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d。
进一步的,所述步骤(D)具体包括:
1)在供热机组热电特性曲线中,将最大主蒸汽流量曲线上各点进行线性拟合,将机组功率随汽轮机抽汽流量的变化特性拟合为:
P=A0+B0×Fc (11)
式(11)中,P为机组功率,单位MW;A0,B0为线性拟合曲线系数;Fc为汽轮机抽汽流量自变量,单位t/h;
2)根据式(11)计算机组供热抽汽量Fcq下机组调峰最大出力值Ptfmax。
本发明的有益效果在于:本发明专利基于供热机组设计热电特性及额定供热设计工况参数,计算确定出汽轮机实际运行低压缸最小排汽量下的特性曲线;基于吸收式热泵系统的质量及能量平衡原理,通过采用低温余热热泵供热状态下的实际运行参数,计算确定热泵的驱动蒸汽流量及总供热抽汽流量。综合机组实际运行低压缸最小排汽量及总供热抽汽流量,通过特性曲线确定出采用低温余热热泵供热改造后机组在供热状态下的调峰最小出力及调峰最大出力。本方法解决了现有采用低温余热热泵供热改造后供热机组在任一供热状态下的调峰出力范围的评估问题,为电网实现对此类供热机组的准确可靠调度提供了重要的参考数据;采用本发明的简化评估方法,也可为电网调度部门确定现有采用低温余热热泵供热改造后供热机组在采暖期不同阶段的运行方式提供重要参考数据,提高现有采用低温余热热泵供热改造后供热机组对电网的调峰能力。
附图说明
图1为典型评估的供热机组的系统示意图。
图2为用于本发明方法的评估装置的示例图。
其中,1、1号低加;2、2号低加;3、3号低加;4、4号低加;5、除氧器;6、高压缸;7、中压缸;8、低压缸;9、1号高加;10、2号高加;11、3号高加;12、热网加热器;13、热网疏水泵;14、热网来回水;15、至热网供水;16、凝结水泵;17、给水泵;18、联通管上供热调整蝶阀;19、低温余热热泵;20、汽轮机抽汽来热泵驱动蒸汽;21、回机组的热泵驱动蒸汽疏水;22、至热泵的热网回水;23、热泵出水至热网;24、至热泵热源来水;25、热泵至热源回水;26、锅炉来主蒸汽;27、锅炉来热再热蒸汽;28、至锅炉冷再热蒸汽;29、至锅炉给水;200、获取模块;201、第一计算模块;202第二计算模块;203、第三计算模块;204、显示模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明公开了一种供热机组采用低温余热热泵供热改造后的调峰能力评估方法,具有下列基本特性:采用以下步骤来确定采用低温余热热泵供热改造后供热机组的调峰出力范围。
(1)评估的供热机组及低温余热热泵供热改造技术具有的系统配置特性
1)供热机组汽轮机为常规抽凝式汽轮机。
2)采暖供汽取自压力较低的汽轮机抽汽。抽汽压力为可调整抽汽,通过低压缸进口联通管上的调整蝶阀进行调节,压力范围0.2~1.0MPa。
3)采暖供汽至热网加热器换热后,热网加热器冷凝水通过疏水泵输送至除氧器进口凝结水管道。
4)供热机组为再热式机组时,采暖供汽均取自再热式汽轮机中压缸抽汽。
5)低温余热热泵供热改造后,热泵的驱动蒸汽取自与热网加热器对应的相同抽汽,无法满足热泵运行要求时,切换至备用的高一级抽汽汽源。
6)评估的供热机组只向热网加热器提供采暖供汽,无对外工业用汽需求,即汽轮机为单抽汽供热汽轮机。
7)采用低温余热热泵供热改造时,增加的供热能力只用于替代采暖抽汽供热。
典型评估的供热机组的系统示意如附图1所示。
(2)确定供热机组热电特性曲线中的考虑适当运行安全裕度的最小等排汽流量曲线。
1)在供热机组热电特性曲线中,将低压缸最小排汽流量限制线(对应低压缸最小排汽流量FLPexmin)进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
y1=a1+b1×x (1)
式(1)中,a1、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y1表示机组功率应变量,单位MW;
2)在供热机组汽轮机设计资料中,查找额定供热工况的以下参数:设计供热抽汽量Fcqdes、设计主蒸汽进汽量Fmsdes、设计低压缸排汽量FLPexdes、以及此工况下供热机组出力Pdes;
3)确定机组在低压缸设计排汽量FLPexdes下的功率随主蒸汽流量的变化特性曲线。方法如下:由等式Pdes=a2+b1×Fmsdes,可求得a2:
a2= Pdes-b1×Fmsdes
由此在低压缸设计排汽量下,机组功率随主蒸汽流量的变化特性拟合为:
y2=a2+b1×x (2)
式(2)中,a2、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y2表示机组功率应变量,单位MW;
4)考虑适当运行安全裕度的实际运行低压缸最小等排汽流量曲线可表示为:
y=y1+(y2-y1)/(FLPexdes- FLPexmin)×(FLPex- FLPexmin)
= y1+(y2-y1)/(FLPexdes- FLPexmin)×50 (3)
式(3)中,y1与y2分别根据公式(1)及公式(2)求得。本发明中考虑实际运行中低压缸排汽流量比设计低压缸最小排汽流量要偏高,并取偏高值为50t/h作为适当运行安全裕度。
(3)根据采用低温余热热泵供热改造后供热机组的供热运行数据计算供热机组的实际供热抽汽流量。
1)收集现场机组供热运行数据:
包括采用低温余热热泵供热改造后机组热泵的热网水流量Frws、热泵热网水进水温度t1、热泵热网水出水温度t2、低温余热热泵热源水进水温度t3、低温余热热泵驱动蒸汽进口压力pds及进口温度tds、热泵凝结水箱温度tdsc以及采用低温余热热泵供热改造后机组汽轮机热网加热器抽汽流量Fcq1;流量单位均为t/h,温度单位均为℃,压力单位为MPa;
2)计算热网水在热泵中的总吸热量:
Q1=Frws×1000×4.186×(t2-t1) (4)
式(4)中,Q1为热网水在热泵中的总吸热量,单位为kJ/h;
3)计算热泵运行中的实际能效系数:
COPr=COPdes×ftrw×ftry×fpds (5)
式(5)中,COPr为实际运行热泵能效系数;COPdes为热泵设计能效系数;ftrw为热泵热网进水温度修正系数,按热泵热网水进水温度t1根据厂家提供的修正曲线查得;ftry为热泵热源水进水温度修正系数,按热泵热源水进水温度t3根据厂家提供的修正曲线查得;fpds为热泵驱动蒸汽压力修正系数,按热泵驱动蒸汽进口压力pds根据厂家提供的修正曲线查得;
4)计算单位质量热泵驱动蒸汽在热泵中的凝结热量:
q2=h1(pds,tds)- h2(tdsc) (6)
式(6)中,q2为每kg热泵驱动蒸汽在热泵中的凝结热量,单位kJ/kg;h1()为水蒸汽特性公式计算函数,得到的驱动蒸汽焓值,单位kJ/kg;h2()为水蒸气特性公式计算函数,得到的饱和水焓值,单位kJ/kg;
5)计算采用低温余热热泵在供热状态下的驱动蒸汽流量Fds:
Fds= Q1/ COPr/q2/1000 (7)
式(7)中,Fds为低温余热热泵在供热状态下的驱动蒸汽流量,单位t/h;
6)计算机组总供热抽汽流量Fcq;
Fcq= Fds+Fcq1 (8)
式(8)中,Fcq为机组总供热抽汽流量,单位t/h;
(4)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最小出力Ptfmin。
具体步骤如下:
1)在供热机组热电特性曲线中,将等抽汽流量曲线进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
zi=ci+di×x(9)
式(9)中,i为由机组等抽汽流量的编号,范围由0至n,随i的增大抽汽流量Fcn(i)逐渐增大,i=0时对应供热抽汽流量为零即纯凝运行工况,i=n时对应最大供热抽汽流量的运行工况;ci、di分别为特性系数;x为机组主蒸汽流量,单位t/h;zi为对应汽轮机采暖抽汽流量Fcn(i)下的机组功率,单位MW。
2)在任意供热抽汽量Fcq时机组功率随主蒸汽流量的变化特性可表示为:
当Fcq(i)<Fcq<Fcq(i+1)时,
z=zi+(zi+1-zi)/( Fcq(i+1)- Fcq(i)) ×( Fcq- Fcq(i)) (10)
式(10)中zi及zi+1均根据公式(9)计算得到。
3)确定供热机组的实际运行低压缸最小等排汽流量曲线及汽轮机供热抽汽流量曲线的相交工况点,即为供热机组的最小调峰出力;步骤如下:
a、假设主蒸汽流量初始值x0为75%额定主蒸汽流量;设定ST=100,j=1;
b、根据式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y, 根据式(10)计算机组供热抽汽量Fcq下机组功率值为z;比较y与 z;
c、如果abs(y- z)<0.01,相交工况点找到,确定的主蒸汽流量x=x0,机组功率Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S1=1,S0=1,主蒸汽流量x1= x0+SP,j=j+1;
如果abs(y- z)≥0.01且y< z时,S1=-1,S0=-1,主蒸汽流量x1= x0-SP,j=j+1;
d、根据新的主蒸汽流量按式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y,根据式(10)计算机组供热抽汽量Fcq下机组功率值为z;比较y与 z;
e、如果abs(y- z)<0.01,相交工况点找到,确定的主蒸汽流量x=x(j-1),机组功率Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S(j)=1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果abs(y- z)≥0.01且y< z时,S(j)=-1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d。
(5)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最大出力。
具体步骤如下:
1)在供热机组热电特性曲线中,将最大主蒸汽流量曲线上各点进行线性拟合,将机组功率随汽轮机抽汽流量的变化特性拟合为:
P=A0+B0×Fc (11)
式(11)中,P为机组功率,单位MW;A0,B0为线性拟合曲线系数;Fc为汽轮机抽汽流量自变量,单位t/h;
2)根据式(11)计算机组供热抽汽量Fcq下机组调峰最大出力值Ptfmax。
(6)采用低温余热热泵供热改造后供热机组的当前调峰出力范围为Ptfmin至Ptfmax。
用于上述评估方法的评估装置的示例图如图2所示,该装置可以包括:获取模块200、第一计算模块201、第二计算模块202、第三计算模块203和显示模块204。
获取模块201,获取额定供热工况参数以及机组供热运行数据。
对低压缸最小排汽流量限制线进行线性拟合,利用查找得到的额定供热工况参数,确定机组在低压缸设计排汽量下的功率随主蒸汽流量的变化特性曲线,进而得到实际运行低压缸最小等排汽流量曲线。
第一计算模块201,利用机组供热运行数据,计算热网水在热泵中的总吸热量、计算热泵运行中的实际能效系数、计算单位质量热泵驱动蒸汽在热泵中的凝结热量、计算采用低温余热热泵在供热状态下的驱动蒸汽流量以及计算机组总供热抽汽流量。
第二计算模块202,根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最小出力。
第三计算模块203,根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最大出力。
显示模块204,将第二计算模块202和第三计算模块203的计算结果进行显示,即采用低温余热热泵供热改造后供热机组的当前调峰出力范围。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换或改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种供热机组采用低温余热热泵供热改造后的调峰能力评估方法,其特征在于,综合供热机组实际运行低压缸最小排汽量及总供热抽汽流量,通过实际运行低压缸最小排汽量下的特性曲线及各供热抽汽流量下的特性曲线确定出采用低温余热热泵供热改造后机组在供热状态下的调峰最小出力及调峰最大出力;
其具体包括如下步骤:
(A)确定供热机组热电特性曲线中的考虑适当运行安全裕度的最小等排汽流量曲线;
(a)在供热机组热电特性曲线中,将低压缸最小排汽流量限制线进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
y1=a1+b1×x (1)
式(1)中,a1、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y1表示低压缸最小排汽流量条件下的机组功率应变量,单位MW;
(b)在供热机组汽轮机设计资料中,查找额定供热工况的以下参数:设计主蒸汽进汽量Fmsdes、设计低压缸排汽量FLPexdes、以及此工况下供热机组出力Pdes;
(c)确定机组在低压缸设计排汽量FLPexdes下的功率随主蒸汽流量的变化特性曲线;
由等式Pdes=a2+b1×Fmsdes;
可计算得到:a2= Pdes-b1×Fmsdes;
由此在低压缸设计排汽量下,机组功率随主蒸汽流量的变化特性拟合为:
y2=a2+b1×x (2)
式(2)中,a2、b1分别为特性系数,x表示主蒸汽流量自变量,单位t/h;y2表示机组功率应变量,单位MW;
(d)考虑实际运行中低压缸排汽流量比设计低压缸最小排汽流量要偏高,并取偏高值为50t/h作为适当运行安全裕度;实际运行低压缸最小等排汽流量曲线可表示为:
y=y1+(y2-y1)/(FLPexdes- FLPexmin)×50 (3);
(B)根据采用低温余热热泵供热改造后供热机组的供热运行数据计算供热机组的实际供热抽汽流量;
1)收集现场机组供热运行数据:
包括采用低温余热热泵供热改造后机组热泵的热网水流量Frws、热泵热网水进水温度t1、热泵热网水出水温度t2、低温余热热泵热源水进水温度t3、低温余热热泵驱动蒸汽进口压力pds及进口温度tds、热泵凝结水箱温度tdsc以及采用低温余热热泵供热改造后机组汽轮机热网加热器抽汽流量Fcq1;流量单位均为t/h,温度单位均为℃,压力单位为MPa;
2)计算热网水在热泵中的总吸热量:
Q1=Frws×1000×4.186×(t2-t1)(4)
式(4)中,Q1为热网水在热泵中的总吸热量,单位为kJ/h;
3)计算热泵运行中的实际能效系数:
COPr=COPdes×ftrw×ftry×fpds(5)
式(5)中,COPr为实际运行热泵能效系数;COPdes为热泵设计能效系数;ftrw为热泵热网进水温度修正系数,按热泵热网水进水温度t1根据厂家提供的修正曲线查得;ftry为热泵热源水进水温度修正系数,按热泵热源水进水温度t3根据厂家提供的修正曲线查得;fpds为热泵驱动蒸汽压力修正系数,按热泵驱动蒸汽进口压力pds根据厂家提供的修正曲线查得;
4)计算单位质量热泵驱动蒸汽在热泵中的凝结热量:
q2=h1(pds,tds)- h2(tdsc) (6)
式(6)中,q2为每kg热泵驱动蒸汽在热泵中的凝结热量,单位kJ/kg;h1()为水蒸汽特性公式计算函数,由蒸汽压力及温度计算蒸汽焓值,单位kJ/kg;h2()为水蒸汽特性公式计算函数,由水温度计算得到饱和水焓值,单位kJ/kg;
5)计算采用低温余热热泵在供热状态下的驱动蒸汽流量Fds:
Fds= Q1/ COPr/q2/1000 (7)
式(7)中,Fds为低温余热热泵在供热状态下的驱动蒸汽流量,单位t/h;
6)计算机组总供热抽汽流量Fcq;
Fcq= Fds+Fcq1 (8)
式(8)中,Fcq为机组总供热抽汽流量,单位t/h;
(C)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最小出力Ptfmin;
1)在供热机组热电特性曲线中,将等抽汽流量曲线进行线性拟合,将功率随主蒸汽流量的变化特性拟合为:
zi=ci+di×x(9)
式(9)中,i为由机组等抽汽流量的编号,范围由0至n,随i的总供热抽汽流量Fcq(i)逐渐增大,i=0时对应供热抽汽流量为零即纯凝运行工况,i=n时对应最大供热抽汽流量的运行工况;ci、di分别为特性系数;x表示主蒸汽流量自变量,单位t/h;zi为对应汽轮机的总供热抽汽流量Fcq(i)下的机组功率,单位MW;
2)在任意总供热抽汽量Fcq时机组功率随主蒸汽流量的变化特性可表示为:
当Fcq(i)<Fcq<Fcq(i+1)时,
z=zi+(zi+1-zi)/( Fcq(i+1)- Fcq(i)) ×( Fcq- Fcq(i)) (10)
式(10)中zi及zi+1均根据公式(9)计算得到;z为总供热抽汽量Fcq时的机组功率,单位MW;
3)确定供热机组的实际运行低压缸最小等排汽流量曲线及汽轮机供热抽汽流量曲线的相交工况点,相交工况点对应的机组功率即为供热机组的最小调峰出力;
确定供热机组的最小调峰出力的具体方法包括:
a、假设主蒸汽流量初始值x0为75%额定主蒸汽流量;设定SP=100,j=1;
b、根据式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y, 根据式(10)计算机组总供热抽汽量Fcq下机组功率值为z;比较y与 z;
c、如果abs(y- z)<0.01,相交工况点找到,对应的主蒸汽流量x=x0,调峰最小出力Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S(1)=1,S(0)=1,主蒸汽流量x(1)= x(0)+SP,j=j+1;
如果abs(y- z)≥0.01且y< z时,S(1)=-1,S(0)=-1,主蒸汽流量x(1)= x(0)-SP,j=j+1;
d、根据新的主蒸汽流量按式(3)计算实际运行低压缸最小等排汽流量FLPexminp下机组功率值为y,根据式(10)计算机组总供热抽汽量Fcq下机组功率值为z;比较y与 z;
e、如果abs(y- z)<0.01,相交工况点找到,对应的主蒸汽流量x=x(j-1),调峰最小出力Ptfmin= z;结束;
如果abs(y- z)≥0.01且y> z时,S(j)=1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)+SP,j=j+1,转至步骤d;
如果abs(y- z)≥0.01且y< z时,S(j)=-1;
如果S(j)×S(j-1)>0,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d;
如果S(j)×S(j-1)<0,SP=SP/2,主蒸汽流量x(j)= x(j-1)-SP,j=j+1,转至步骤d;
(D)根据供热机组汽轮机总供热抽汽流量,在供热特性曲线上确定机组对应的调峰最大出力Ptfmax;
(E)采用低温余热热泵供热改造后供热机组的当前调峰出力范围为:Ptfmin至Ptfmax。
2.根据权利要求1所述的评估方法,其特征在于,所述步骤(D)具体包括:
1)在供热机组热电特性曲线中,将最大主蒸汽流量曲线上各点进行线性拟合,将机组功率随汽轮机抽汽流量的变化特性拟合为:
P=A0+B0×Fc (11)
式(11)中,P为机组功率,单位MW;A0,B0为线性拟合曲线系数;Fc为汽轮机抽汽流量自变量,单位t/h;
2)根据式(11)计算机组总供热抽汽量Fcq下机组调峰最大出力值Ptfmax。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911212619.8A CN111047463B (zh) | 2019-12-02 | 2019-12-02 | 供热机组采用低温余热热泵供热改造后的调峰能力评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911212619.8A CN111047463B (zh) | 2019-12-02 | 2019-12-02 | 供热机组采用低温余热热泵供热改造后的调峰能力评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111047463A CN111047463A (zh) | 2020-04-21 |
CN111047463B true CN111047463B (zh) | 2023-06-23 |
Family
ID=70234367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911212619.8A Active CN111047463B (zh) | 2019-12-02 | 2019-12-02 | 供热机组采用低温余热热泵供热改造后的调峰能力评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111047463B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112053027B (zh) * | 2020-07-17 | 2024-05-24 | 国网山东省电力公司电力科学研究院 | 一种热电厂最小运行方式的在线核定方法及系统 |
CN112231908B (zh) * | 2020-10-15 | 2023-03-21 | 国网山东省电力公司电力科学研究院 | 一种抽凝式机组抽汽流量对应负荷上下限的确定方法 |
CN113095623B (zh) * | 2021-03-12 | 2022-08-05 | 国网河北能源技术服务有限公司 | 一种双抽供热机组的调峰能力评估方法 |
CN113065745B (zh) * | 2021-03-12 | 2022-07-22 | 国网河北能源技术服务有限公司 | 双抽供热机组最小调峰出力确定方法、装置及终端 |
CN113191566B (zh) * | 2021-05-20 | 2023-03-14 | 西安热工研究院有限公司 | 热泵梯级供热机组最佳运行方式的在线确定系统及方法 |
CN113469433B (zh) * | 2021-06-29 | 2023-04-07 | 西安热工研究院有限公司 | 一种热再可调供汽方式下抽汽压力的寻优方法 |
CN113489024B (zh) * | 2021-07-13 | 2023-05-26 | 西安热工研究院有限公司 | 热电联供机组多模式抽汽辅助调峰调频控制系统及方法 |
CN113566186B (zh) * | 2021-07-28 | 2023-08-29 | 西安热工研究院有限公司 | 一种调节级压力测点损坏后的锅炉汽包水位调整方法 |
CN115660325B (zh) * | 2022-10-08 | 2023-05-30 | 国网山东省电力公司威海供电公司 | 一种电网调峰能力量化方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101798941A (zh) * | 2010-01-08 | 2010-08-11 | 华北电力大学 | 最佳冷源热网加热器及其参数的确定方法 |
CN106096310A (zh) * | 2016-06-27 | 2016-11-09 | 国网山西省电力公司电力科学研究院 | 一种抽汽供热机组供热工况热力性能计算方法 |
CN108090663A (zh) * | 2017-12-11 | 2018-05-29 | 囯网河北省电力有限公司电力科学研究院 | 供热机组深度调峰最低出力的评估方法及系统 |
CN108695898A (zh) * | 2018-06-29 | 2018-10-23 | 国网山东省电力公司电力科学研究院 | 一种热电机组调峰能力的在线监测方法 |
CN109063890A (zh) * | 2018-06-21 | 2018-12-21 | 国网山东省电力公司电力科学研究院 | 一种基于热电厂全厂调峰能力最大化的热负荷分配方法 |
CN109236391A (zh) * | 2018-09-11 | 2019-01-18 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种热电联产机组调峰范围的计算方法 |
CN109377018A (zh) * | 2018-09-27 | 2019-02-22 | 国网河北省电力有限公司电力科学研究院 | 供热机组调峰能力评估方法 |
-
2019
- 2019-12-02 CN CN201911212619.8A patent/CN111047463B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101798941A (zh) * | 2010-01-08 | 2010-08-11 | 华北电力大学 | 最佳冷源热网加热器及其参数的确定方法 |
CN106096310A (zh) * | 2016-06-27 | 2016-11-09 | 国网山西省电力公司电力科学研究院 | 一种抽汽供热机组供热工况热力性能计算方法 |
CN108090663A (zh) * | 2017-12-11 | 2018-05-29 | 囯网河北省电力有限公司电力科学研究院 | 供热机组深度调峰最低出力的评估方法及系统 |
CN109063890A (zh) * | 2018-06-21 | 2018-12-21 | 国网山东省电力公司电力科学研究院 | 一种基于热电厂全厂调峰能力最大化的热负荷分配方法 |
CN108695898A (zh) * | 2018-06-29 | 2018-10-23 | 国网山东省电力公司电力科学研究院 | 一种热电机组调峰能力的在线监测方法 |
CN109236391A (zh) * | 2018-09-11 | 2019-01-18 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种热电联产机组调峰范围的计算方法 |
CN109377018A (zh) * | 2018-09-27 | 2019-02-22 | 国网河北省电力有限公司电力科学研究院 | 供热机组调峰能力评估方法 |
Non-Patent Citations (1)
Title |
---|
蔡东旭 ; 邵志祥 ; 叶诚广 ; .300MW汽轮机更新改造中新技术的应用.发电设备.(第03期),论文正文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111047463A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111047463B (zh) | 供热机组采用低温余热热泵供热改造后的调峰能力评估方法 | |
CN111047168B (zh) | 一种供热机组高背压供热改造后的调峰能力评估方法 | |
CN110925037B (zh) | 一种考虑运行安全裕量的采暖供热机组实际调峰能力评价方法 | |
CN100354504C (zh) | 一种火电机组多级利用回热疏水余热发电装置 | |
CN100402814C (zh) | 一种天然气冷热电联供的烟气低温端热利用系统 | |
CN110930050B (zh) | 供热机组采用储热罐技术灵活性改造后的调峰能力改善评估方法 | |
CN110991877B (zh) | 供热机组采用低压缸切缸灵活性改造后供热及调峰能力改善评估方法 | |
CN104763485B (zh) | 一种补热型超高压/亚临界背压供热机组热力系统 | |
EP2423471A2 (en) | Systems and methods involving combined cycle plants | |
CN105787211A (zh) | 针对燃气透平劣化的联合循环余热锅炉压力调整方法 | |
CN204254934U (zh) | 一种利用压缩式热泵实现锅炉烟气余热深度回收的供热系统 | |
CN208124429U (zh) | 一种供热机组抽汽余热回收利用系统 | |
CN113268887B (zh) | 一种煤电机组锅炉暖风器最佳热源确定系统及方法 | |
CN106437875B (zh) | 火电机组工质分流循环调峰系统 | |
CN108167915A (zh) | 一种与调峰锅炉结合的大温差供热系统及方法 | |
CN113217974A (zh) | 一种超高压工业供热热力系统 | |
Jamil et al. | An exergoeconomic and normalized sensitivity based comprehensive investigation of a hybrid power-and-water desalination system | |
CN108843413B (zh) | 一种可再生能源调峰机组的调峰补偿基准值的计算方法 | |
CN208687705U (zh) | 一种切除低压缸进汽运行与低压省煤器耦合的系统 | |
CN108980616B (zh) | 一种针对间歇性用汽用户的长距离工业供汽系统 | |
CN108072030A (zh) | 用于高参数二次再热机组的回热系统 | |
RU62166U1 (ru) | Система охлаждения отработанного пара паровых турбин | |
CN209945054U (zh) | 一种汽轮发电机组凝汽器水冷变为空冷的系统装置 | |
CN202883026U (zh) | 发电厂回热式汽动凝结水泵系统 | |
CN208702473U (zh) | 闭式水回热加热系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |