CN113361828A - 多机组、多供热模式的热电厂供热负荷分配优化方法 - Google Patents

多机组、多供热模式的热电厂供热负荷分配优化方法 Download PDF

Info

Publication number
CN113361828A
CN113361828A CN202110831724.0A CN202110831724A CN113361828A CN 113361828 A CN113361828 A CN 113361828A CN 202110831724 A CN202110831724 A CN 202110831724A CN 113361828 A CN113361828 A CN 113361828A
Authority
CN
China
Prior art keywords
unit
entering
back pressure
heat supply
load
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.)
Granted
Application number
CN202110831724.0A
Other languages
English (en)
Other versions
CN113361828B (zh
Inventor
马汀山
石慧
吕凯
居文平
王妍
许朋江
薛朝囡
邓佳
张建元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Thermal Power Research Institute Co Ltd
Xian Xire Energy Saving Technology Co Ltd
Original Assignee
Xian Thermal Power Research Institute Co Ltd
Xian Xire Energy Saving Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Thermal Power Research Institute Co Ltd, Xian Xire Energy Saving Technology Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202110831724.0A priority Critical patent/CN113361828B/zh
Publication of CN113361828A publication Critical patent/CN113361828A/zh
Priority to PCT/CN2022/101319 priority patent/WO2023000927A1/zh
Application granted granted Critical
Publication of CN113361828B publication Critical patent/CN113361828B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06311—Scheduling, planning or task assignment for a person or group
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E20/00—Combustion technologies with mitigation potential
    • Y02E20/14—Combined heat and power generation [CHP]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)

Abstract

本发明公开了多机组、多供热模式的热电厂供热负荷分配优化方法，本发明提出了一种适用于多机组、多供热模式的热电联产电厂在分机组电负荷和总供热负荷实时调度以及供热回水温度、供热循环水流量等多变量约束条件下的厂级运行优化方法，以在供热和供电负荷条件下厂级标煤消耗总量最低值为寻优目标函数，得出高背压供热机组运行背压、其他机组承担的热负荷以及相应的供热模式。本发明寻优目标直观且符合生产实际，大幅简化寻优过程，扩大了热电联产机组运行台数以及供热模式，具有广阔的应用前景。

Description

多机组、多供热模式的热电厂供热负荷分配优化方法
技术领域
本发明属于燃煤发电机组热电联产供热领域，具体涉及一种多机组、多供热模式的热电厂供热负荷分配优化方法。
背景技术
对于热电联产企业，发电和供热的合理分配对降低热电联产的运营成本意义重大，而拥有多机组、多供热模式的热电联产企业面临着厂级的运行优化问题，其中供热负荷在多机组、多供热模式间的分配对热电厂的节能降耗指标影响很大。
目前，已有相关专业人员对热电联产机组的运行优化进行了一些研究，但并无关于2台及以上、2种及以上供热模式供热机组间的热负荷优化分配问题的深入研究。文献《配置储热罐后热电联产机组运行优化》介绍了热电联产机组配置储热罐后的运行优化问题，研究对象为热电厂的两台机组，机组型式完全相同，供热方式仅有中低压连通管打孔抽汽，研究对象过于单一；硕士论文《多机组联合抽汽供热时热负荷分配优化研究》研究了案例电厂在不同发电负荷下热负荷的分配原则，但各计算工况供热负荷一定，并没有考虑热电厂实际供热时在严寒期、平均期和初末期供热负荷的变化；文献《双机联调抽汽-高背压联合供热分析与优化》研究了抽汽供热、抽汽-高背压联合供热、双机联调抽汽-高背压联合供热技术三者的供热能力，但研究侧重点在三种供热技术供热能力的对比，并没有对供热负荷的优化分配进行深入研究；硕士论文《高背压供热机组节能分析与运行优化》研究了发电负荷、供热负荷、供回水温度和背压对高背压供热机组性能的影响，研究对象为某300MW热电联产机组，供热方式仅有高背压供热，研究范围相对狭小。
因此，为指导多机组(2台及以上)、多供热模式(2种及以上)的热电联产企业供热负荷的分配优化，亟需提出一种适用于多机组、多供热模式的热电厂供热负荷分配优化方法。
发明内容
本发明的目的在于克服上述不足，提供一种多机组、多供热模式的热电厂供热负荷分配优化方法，以对外供热负荷和分机组电负荷给定条件下，热电联产机组厂级标煤消耗总量最低值对应的工况作为最优运行工况，列出高背压供热机组运行背压、其他机组承担的热负荷以及相应的供热模式。
为了达到上述目的，包括以下步骤：
S1，确定多机组、多供热模式的热电厂供热负荷分配优化方法的寻优目标函数，寻优目标函数包括全厂标煤消耗总量；
S2，确定各热电联产机组不同供热方式下的电负荷-供热负荷、电负荷-供热负荷-标煤消耗的关联特性；
S3，根据全厂标煤消耗总量，进行供热负荷和不同分机组电负荷给定条件下的厂级运行优化。
S1中，全厂总盈利值M的计算方法如下：
M＝H+E-C＝Qt×h+(Nge1+Nge2+Nge3)×e-(B1+B2+B3)×c
其中，H为售热收入；Qt为总对外供热负荷，h为出厂热价，E为售电收入，Nge1、Nge2、Nge3分别为三个机组的实时发电负荷，e为上网电价，C为标煤消耗成本，B1、B2、B3分别为三个机组的实时标煤消耗量，c为标煤单价。
S2中，当机组采用高背压梯级供热，供热负荷Q计算方法如下：
Figure BDA0003175691790000021
ts＝f(Pc)＝-0.029×Pc 2+2.28×Pc+26.13
其中，mcw为进出高背压热网凝汽器的热网循环水流量，Cp为热网循环水定压比热容，t0、t1分别为进出高背压热网凝汽器的热网循环水温度，δt为高背压热网凝汽器的换热端差，ts为乏汽在高背压热网凝汽器放热后的凝结水温度，Pc为高背压机组的运行背压；
供热负荷Q，是电负荷Nge1、进出高背压热网凝汽器的热网循环水温度t0及t1、热网循环水流量mcw和运行背压Pc的多元函数，通过下式表征：
Q＝f1(Nge1,Pc,mcw,t0)
高背压供热机组电、热双供应条件下的标煤消耗量B是电负荷Nge1和运行背压Pc的二元函数，通过下式表征：
B＝F1(Nge1,Pc)。
S2中，当机组采用中排抽汽模式时，电负荷Nge确定，供热负荷介于Q在0和最大值Qmax之间：
0≤Q≤Qmax＝f2(Nge)
标煤消耗量B是电负荷Nge和供热负荷Q的二元函数，通过下式表征：
B＝F2(Nge,Q)。
S2中，当机组采用低压缸零出力供热模式时，机组以热定电运行，供热负荷Q是电负荷Nge的一元线性函数，通过下式表征：
Q＝f3(Nge)
标煤消耗量B是电负荷Nge和供热负荷Q的二元函数，通过下式表征：
B＝F3(Nge)
锅炉设计出力Dms0条件下，机组采用低压缸零出力供热模式，电负荷运行范围为：
0≤Nge≤Nge,max＝f4(Dms0)。
S3的具体优化方法如下：
S301，输入边界参数，边界参数包括供热回水温度t0、热网循环水流量mcw、总供热负荷Qt以及一号机组发电负荷Nge1、二号机组发电负荷Nge2和三号机组发电负荷Nge3；
S302，制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S303，计算Q1＝f1(Nge,Pc0,t0,mcw)，其中，Q1为一号机组供热负荷，Nge为实时发电负荷，Pc0为运行背压，t0为进高背压热网凝汽器的热网循环水温度，mcw为进出高背压热网凝汽器的热网循环水流量；
S304，针对二号机组，判断是否具备投运中排抽汽和低压缸零出力两种模式的电负荷条件，若是，则进入S305，若否，则执行S315；
S305，针对三号机组，判断是否具备投运中排抽汽或低压缸零出力两种模式的电负荷条件，若是，则进入S306，若否，则执行S311和S315；
S306，对二号机组进行热负荷判定，判断二号机组是否具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，若是，则进入S307，若否，则执行S315；
S307，对三号机组进行热负荷判定，判断三号机组是否具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，若是，则进入S308，若否，则执行S315；
S308，进行热负荷判定，判断f2-2(Nge2)+f3-3(Nge3)是否大等于Qt-Q1，若是，则进入S309，若否，则终止寻优迭代流程；
S309，进行热负荷判定，判断f2-3(Nge3)+f3-2(Nge2)是否大于等于Qt-Q1，若是，则进入S310，若否，则终止寻优迭代流程；
S310，二号机组和三号机组均具备中排抽汽或低压缸零出力两种供热模式，进入第一迭代寻优模式，输出最优结果；
S311，判断二号机组是否具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，同时，三号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，则进入S312，若否，则终止寻优迭代流程；
S312，进入第二迭代寻优模式，输出最优结果；
S313，判断二号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，同时，三号机组是否具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，若是，则进入S314，若否，则终止寻优迭代流程；
S314，进行第三寻优模式，输出最优结果；
S315，判断二号机组和三号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，若是，则进入S316，若否，则终止寻优迭代流程；
S316，进行第四迭代寻优模式，输出最优结果。
第一迭代寻优模式包括以下步骤：
S10-1-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升5℃对应背压Pc0，作为迭代寻优的背压基准。
S10-1-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-1-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S10-1-4：判断Q30是否大于0，若是，则进入S10-1-5；若否，则转入S10-1-11；
S10-1-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))，进入S10-1-7；
S10-1-7：判断Q21是否小于0，或Q31是否大于等于f3-3(Nge3)，若是，则寻优迭代流程终止；若否，则转入S10-1-8；
S10-1-8：判断Q31是否小于等于f2-3(Nge3)，若是，则进入S10-1-9；若否，则转入S10-1-15；
S10-1-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-1-6；
S10-1-11：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S10-1-12；若否，则转入S10-1-14；
S10-1-12：令Q20＝Qt-Q10，则Q30＝0，进入S10-1-13；
S10-1-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-14：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S10-1-13；
S10-1-15：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F3-3(Nge3)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组的热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-1-17；
S10-1-17：令Pc1＝Pc0+1kPa，进入S10-1-18；
S10-1-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-1-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-1-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b；
S10-2-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S10-2-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-2-3：令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S10-2-4：判断Q20是否大于等于0，若是，则进入S10-2-5；若否，则转入S10-2-11；
S10-2-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))，进入S10-2-7；
S10-2-7：判断Q31是否小于等于0，或Q21大于等于f3-2(Nge2)，若是，则终止寻优迭代流程；若否，则转入S10-2-8；
S10-2-8：判断Q21是否小于等于f2-2(Nge2)，若是，则进入S10-2-9；若否，则转入S10-2-15；
S10-2-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-2-6；
S10-2-11：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S10-2-12；若否，则转入S10-2-14；
S10-2-12：令Q30＝Qt-Q10，则Q20＝0，进入S10-2-13；
S10-2-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-14：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f2-3(Nge3)，进入S10-2-13；
S10-2-15：计算：B10＝F1(Nge,Pc0)、B21＝F3-2(Nge2、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-2-17；
S10-2-17：令Pc1＝Pc0+1kPa，进入S10-2-18；
S10-2-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-2-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-2-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b；
第二寻优模式的具体方法如下：
S12-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S12-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S12-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S12-4：判断Q30是否大于等于0，若是，则进入S12-5；若否，则转入S12-10；
S12-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S12-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))；
S12-7：判断Q21是否小于等于0，或Q31大于等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S12-8；
S12-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S12-9；
S12-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S12-6；
S12-10：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S12-11；若否，则转入S12-13；
S12-11：令Q20＝Qt-Q10，则Q30＝0，进入S12-12；
S12-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S12-6；
S12-13：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S12-12；
S12-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S12-15；
S12-15：令Pc1＝Pc0+1kPa，进入S12-16；
S12-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S12-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S12-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
第三寻优模式的具体方法如下：
S14-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S14-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S14-3令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S14-4判断Q20是否大于等于0，若是，则进入S14-5；若否，则转入S14-10；
S14-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30；
S14-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))；
S14-7：判断Q31是否小于等于0，或Q21是否大于等于f2-2(Nge2)，若是，则终止寻优迭代流程；若否，则进入S14-8；
S14-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S14-9；
S14-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S14-6；
S14-10：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S14-11；若否，则转入S14-13；
S14-11：令Q30＝Qt-Q10，则Q20＝0，进入S14-12；
S14-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S14-6；
S14-13：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)，进入S14-12；
S14-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S14-15；
S14-15：令Pc1＝Pc0+1kPa，进入S14-16；
S14-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S14-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S14-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
第四寻优模式的具体方法如下：
S16-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S16-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S16-3令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)；
S16-4判断Q30是否大于等于0，若是，则进入S16-5；若否，则转入S16-10；
S16-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S16-6：令Q21＝Q20-0.1×Q20则Q31＝Q30+0.1×Q20；
S16-7：判断Q21是否等于0，或Q31是否等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S16-8；
S16-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S16-9；
S16-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S16-6；
S16-10：令Q20＝Qt-Q10、Q30＝0，进入S16-11；
S16-11：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S16-6；
S12-12：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S16-13；
S16-13：令Pc1＝Pc0+1kPa，进入S16-14；
S16-14：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S16-2；否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S16-15：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
与现有技术相比，本发明提出了一种适用于多机组(2台以上)、多供热模式(2种以上)的热电联产电厂在分机组电负荷和总供热负荷实时调度以及供热回水温度、供热循环水流量等多变量约束条件下的厂级运行优化方法，以在供热和供电负荷条件下厂级标煤消耗总量最低值为寻优目标函数，得出高背压供热机组运行背压、其他机组承担的热负荷以及相应的供热模式。现有的采用热耗率、煤耗等为目标的传统寻优方法，仅适用于机组台数不超过2台、供热模式单一的场景，计算过程繁琐，且难以为相关领域技术人员理解和接受。与传统方法相比，本发明寻优目标直观且符合生产实际，大幅简化寻优过程，扩大了热电联产机组运行台数以及供热模式，具有广阔的应用前景。
附图说明
图1为全厂供热系统示意图；
图2为本发明的流程图；
图3为第一寻优模式的流程图；
图4为第一寻优模式的流程图；
图5为第二寻优模式的流程图；
图6为第三寻优模式的流程图；
图7为第四寻优模式的流程图；
其中，1、高中压缸，2、低压缸，3、高背压热网凝汽器，4、热网循环水泵，5、热网加热器，6、供热蝶阀。
具体实施方式
下面结合附图对本发明做进一步说明。
实施例：
该供热电厂由四台机组组成，编号分别为一号机组、二号机组、三号机组和四号机组，其中一号机组和二号机组参数及容量相同，称为一期机组；三号机组和四号机组参数及容量参数相同，称为二期机组；一期和二期的机组参数及容量可以相同，也可以不同。为使得本发明提出的方法更具实用性，约定案例中一期和二期的机组参数及容量不同。
一号机组高背压供热，二号机组有中排抽汽和低压缸零出力供热两种模式，三号机组有中排抽汽和低压缸零出力供热两种模式，四号机组纯凝运行，不承担对外供热。全厂供热系统流程见附图1。供热回水经热网循环水泵4升压后先进入高背压热网凝汽器3，吸收一号机组汽轮机低压缸2排汽余热，再进入热网加热器5由二号机组和三号机组高中压缸1排汽尖峰提温后，对外供出。
本实施例提出的多机组、多供热模式的热电厂供热负荷分配优化方法，不涉及四号机组。
步骤1：确定多机组、多供热模式的热电厂供热负荷分配优化方法的寻优目标函数。
燃煤供热电厂生产盈利值M按式(1)计算。
M＝H+E-C＝Qt×h+(Nge1+Nge2+Nge3)×e-(B1+B2+B3)×c (1)
式中，H为售热收入；Qt为总对外供热负荷，MW；h为出厂热价，元/MW；
E为售电收入；Nge1、Nge2、Nge3分别为一号机组、二号机组和三号机组的实时发电负荷，MW；e为上网电价，元/MW；
C为标煤消耗成本；B1、B2、B3分别为一号机组、二号机组和三号机组的实时标煤消耗量，t；c为标煤单价，元/t；B＝B1+B2+B3。
燃煤供热电厂对外供热和供电，总供热负荷Qt由热用户根据需求实时调度，分机组供电负荷Nge1、Nge2和Nge3由电网根据供需关系实时调度，燃煤供热电厂没有独立的热、电负荷调度权限。
在出厂热价h、上网电价e、标煤单价c等经济性关联数据受市场供需关系影响，短期内一般不发生变化。
综上，多机组、多供热模式的热电厂供热负荷分配优化方法的寻优目标函数为全厂标煤消耗总量B，在总供热负荷Qt和分机组电负荷Nge1、Nge2、Nge3约束条件下，通过改变总供热负荷Qt在各机组的分配，使得全厂标煤消耗总量B最低，可达到盈利能力最高的目的。
步骤2：得出各热电联产机组不同供热方式下的电负荷-供热负荷、电负荷-供热负荷-标煤消耗关联特性。
通过性能试验和理论计算相结合的方法，分别得出一号机组高背压梯级供热、二号机组中排抽汽和低压缸零出力供热、三号机组中排抽汽和低压缸零出力供热的的电负荷-供热负荷和电负荷-供热负荷-标煤消耗关联特性
1)一号机组高背压梯级供热，供热负荷Q按式(2)计算。
Figure BDA0003175691790000141
式中，mcw为进出高背压热网凝汽器的热网循环水流量，t/h；Cp为热网循环水定压比热容，kJ/kg·K；t0、t1分别为进出高背压热网凝汽器的热网循环水温度，℃；δt为高背压热网凝汽器的换热端差，℃；ts为乏汽在高背压热网凝汽器放热后的凝结水温度，℃；Pc为高背压机组的运行背压，kPa。
供热负荷Q，是电负荷Nge1、进出高背压热网凝汽器的热网循环水温度t0及t1、热网循环水流量mcw和运行背压Pc的多元函数，以式(3)表征。
Q＝f1(Nge1,Pc,mcw,t0) (3)
高背压供热机组电、热双供应条件下的标煤消耗量B是电负荷Nge1和运行背压Pc的二元函数，以式(4)表征。
B＝F1(Nge1,Pc) (4)
2)2、三号机组同时具备中排抽汽和低压缸零出力供热两种模式。
中排抽汽模式，电负荷Nge给定，供热负荷Q在0和最大值Qmax之间灵活可调：0≤Q≤Qmax＝f2(Nge)，此时标煤消耗量B是电负荷Nge和供热负荷Q的二元函数：B＝F2(Nge,Q)。
切换为低压缸零出力供热模式，机组以热定电运行，供热负荷Q是电负荷Nge的一元线性函数：Q＝f3(Nge)，此时标煤消耗量B是电负荷Nge和供热负荷Q的二元函数：B＝F3(Nge)。锅炉设计出力Dms0条件下，机组采用低压缸零出力供热模式，电负荷运行范围：0≤Nge≤Nge,max＝f4(Dms0)。
为区分二号机组、三号机组，在函数f2、F2、f3、F3增设下标。二号机组为：f2-2、F2-2、f3-2、F3-2。4号机组为：f2-3、F2-3、f3-3、F3-3。
步骤3：以全厂标煤消耗总量B最低，进行总供热负荷Qt和分机组电负荷Nge1、Nge2、Nge3给定条件下的厂级运行优化。
附图2-7给出了本方法的实施过程。
附图2给出了本方法的实施总图。
S1：输入边界参数，边界参数包括供热回水温度t0、热网循环水流量mcw、总供热负荷Qt以及一号机组发电负荷Nge1、二号机组发电负荷Nge2和三号机组发电负荷Nge3；
S2：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升5℃对应背压Pc0，作为迭代寻优的背压基准。
S3：计算Q1＝f1(Nge,Pc0,t0,mcw)
S4：针对二号机组，进行判断Nge2是否小于等于f4-2(Dms0-2)，若是，则认定二号机组具备投运中排抽汽和低压缸零出力两种模式的电负荷条件，并进入S5；否转入S15。
S5：针对三号机组，进行判断Nge3是否小于等于f4-3(Dms0-3)，若是，则认定三号机组具备投运中排抽汽或低压缸零出力两种模式的电负荷条件，并进入S6；否转入S11和S15。
S6：热负荷判定：f3-2(Nge2)是否小于等于Qt-Q1，若是，则认定二号机组具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，并进入S7；否转入S15
S7：热负荷判定：f3-3(Nge3)是否小于等于Qt-Q1，若是，则认定三号机组具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，并进入S8；否转入S11和S15
S8：进行热负荷判定，判断f2-2(Nge2)+f3-3(Nge3)是否大等于Qt-Q1，若是，则进入S309，若否，则终止寻优迭代流程；
S9：进行热负荷判定，判断f2-3(Nge3)+f3-2(Nge2)是否大于等于Qt-Q1，若是，则进入S310，若否，则终止寻优迭代流程；
S10：二号机组和三号机组均具备中排抽汽或低压缸零出力两种供热模式，进入迭代寻优模式1-1和1-2，输出最优结果
S10-1-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升5℃对应背压Pc0，作为迭代寻优的背压基准。
S10-1-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-1-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S10-1-4：判断Q30是否大于0，若是，则进入S10-1-5；若否，则转入S10-1-11；
S10-1-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))，进入S10-1-7；
S10-1-7：判断Q21是否小于0，或Q31是否大于等于f3-3(Nge3)，若是，则寻优迭代流程终止；若否，则转入S10-1-8；
S10-1-8：判断Q31是否小于等于f2-3(Nge3)，若是，则进入S10-1-9；若否，则转入S10-1-15；
S10-1-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-1-6；
S10-1-11：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S10-1-12；若否，则转入S10-1-14；
S10-1-12：令Q20＝Qt-Q10，则Q30＝0，进入S10-1-13；
S10-1-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-14：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S10-1-13；
S10-1-15：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F3-3(Nge3)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组的热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-1-17；
S10-1-17：令Pc1＝Pc0+1kPa，进入S10-1-18；
S10-1-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-1-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-1-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
S10-2寻优模式1-2的实施流程见附图4
S10-2-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S10-2-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-2-3：令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S10-2-4：判断Q20是否大于等于0，若是，则进入S10-2-5；若否，则转入S10-2-11；
S10-2-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))，进入S10-2-7；
S10-2-7：判断Q31是否小于等于0，或Q21大于等于f3-2(Nge2)，若是，则终止寻优迭代流程；若否，则转入S10-2-8；
S10-2-8：判断Q21是否小于等于f2-2(Nge2)，若是，则进入S10-2-9；若否，则转入S10-2-15；
S10-2-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-2-6；
S10-2-11：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S10-2-12；若否，则转入S10-2-14；
S10-2-12：令Q30＝Qt-Q10，则Q20＝0，进入S10-2-13；
S10-2-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-14：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f2-3(Nge3)，进入S10-2-13；
S10-2-15：计算：B10＝F1(Nge,Pc0)、B21＝F3-2(Nge2、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-2-17；
S10-2-17：令Pc1＝Pc0+1kPa，进入S10-2-18；
S10-2-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-2-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-2-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
S10-3：比较S10-1和S10-2迭代寻优过程的全厂总耗煤量最低值，取小值对应工况作为最优运行工况：一号机组运行背压Pcb、2和三号机组热负荷分配Q2b、Q3b。
S11：S4的是判断和S6的是判断，共同转入S11，进行热负荷判断：f2-2(Nge2)+f2-3(Nge3)是否大于等于Qt-Q1，若是，则认定二号机组具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，三号机组具备投运中排抽汽模式的电负荷和热负荷条件，进入S12；见附图2的圆点线所示；若否，则寻优迭代流程终止；
S12：进入迭代寻优模式2，输出最优结果。
寻优模式2的实施流程见附图5。
S12-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S12-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S12-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S12-4：判断Q30是否大于等于0，若是，则进入S12-5；若否，则转入S12-10；
S12-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S12-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))；
S12-7：判断Q21是否小于等于0，或Q31大于等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S12-8；
S12-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S12-9；
S12-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S12-6；
S12-10：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S12-11；若否，则转入S12-13；
S12-11：令Q20＝Qt-Q10，则Q30＝0，进入S12-12；
S12-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S12-6；
S12-13：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S12-12；
S12-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S12-15；
S12-15：令Pc1＝Pc0+1kPa，进入S12-16；
S12-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S12-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S12-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
S13：S5的是判断和S7的是判断，共同转入S13，进行热负荷判断：f2-2(Nge2)+f2-3(Nge3)是否大于等于Qt-Q1，若是，则认定二号机组具备投运中排抽汽模式的电负荷和热负荷条件，三号机组具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，进入S14；见附图2的虚线所示；若否，则寻优迭代流程终止
S14：进入迭代寻优模式3，输出最优结果
寻优模式3的实施流程见附图6
S14-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S14-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S14-3令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S14-4判断Q20是否大于等于0，若是，则进入S14-5；若否，则转入S14-10；
S14-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30；
S14-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))；
S14-7：判断Q31是否小于等于0，或Q21是否大于等于f2-2(Nge2)，若是，则终止寻优迭代流程；若否，则进入S14-8；
S14-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S14-9；
S14-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S14-6；
S14-10：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S14-11；若否，则转入S14-13；
S14-11：令Q30＝Qt-Q10，则Q20＝0，进入S14-12；
S14-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S14-6；
S14-13：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)，进入S14-12；
S14-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S14-15；
S14-15：令Pc1＝Pc0+1kPa，进入S14-16；
S14-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S14-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S14-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
S15：进行热负荷判断：f2-2(Nge2)+f2-3(Nge3)是否大于等于Qt-Q1，若是，则认定二号机组和三号机组具备投运中排抽汽模式的电负荷和热负荷条件，进入S16；否，寻优迭代流程终止
S16：进入迭代寻优模式4，输出最优结果
寻优模式4的实施流程见附图7
S16-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S16-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S16-3令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)；
S16-4判断Q30是否大于等于0，若是，则进入S16-5；若否，则转入S16-10；
S16-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S16-6：令Q21＝Q20-0.1×Q20则Q31＝Q30+0.1×Q20；
S16-7：判断Q21是否等于0，或Q31是否等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S16-8；
S16-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S16-9；
S16-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S16-6；
S16-10：令Q20＝Qt-Q10、Q30＝0，进入S16-11；
S16-11：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S16-6；
S12-12：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S16-13；
S16-13：令Pc1＝Pc0+1kPa，进入S16-14；
S16-14：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S16-2；否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S16-15：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
S10、S12、S14、S16的输出结果为最优运行方式：高背压供热机组运行背压、其他机组承担的热负荷以及相应的供热模式。

Claims (10)

1.多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，包括以下步骤：
S1，确定多机组、多供热模式的热电厂供热负荷分配优化方法的寻优目标函数，寻优目标函数包括全厂标煤消耗总量；
S2，确定各热电联产机组不同供热方式下的电负荷-供热负荷、电负荷-供热负荷-标煤消耗的关联特性；
S3，根据全厂标煤消耗总量，进行供热负荷和不同分机组电负荷给定条件下的厂级运行优化。
2.根据权利要求1所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，S1中，全厂总盈利值M的计算方法如下：
M＝H+E-C＝Qt×h+(Nge1+Nge2+Nge3)×e-(B1+B2+B3)×c
其中，H为售热收入；Qt为总对外供热负荷，h为出厂热价，E为售电收入，Nge1、Nge2、Nge3分别为三个机组的实时发电负荷，e为上网电价，C为标煤消耗成本，B1、B2、B3分别为三个机组的实时标煤消耗量，c为标煤单价。
3.根据权利要求1所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，S2中，当机组采用高背压梯级供热，供热负荷Q计算方法如下：
Figure FDA0003175691780000011
ts＝f(Pc)＝-0.029×Pc 2+2.28×Pc+26.13
其中，mcw为进出高背压热网凝汽器的热网循环水流量，Cp为热网循环水定压比热容，t0、T1分别为进出高背压热网凝汽器的热网循环水温度，δT为高背压热网凝汽器的换热端差，ts为乏汽在高背压热网凝汽器放热后的凝结水温度，Pc为高背压机组的运行背压；
供热负荷Q，是电负荷Nge1、进出高背压热网凝汽器的热网循环水温度t0及T1、热网循环水流量mcw和运行背压Pc的多元函数，通过下式表征：
Q＝f1(Nge1，Pc，mcw，t0)
高背压供热机组电、热双供应条件下的标煤消耗量B是电负荷Nge1和运行背压Pc的二元函数，通过下式表征：
B＝F1(Nge1，Pc)。
4.根据权利要求1所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，S2中，当机组采用中排抽汽模式时，电负荷Nge确定，供热负荷介于Q在0和最大值Qmax之间：
0≤Q≤Qmax＝f2(Nge)
标煤消耗量B是电负荷Nge和供热负荷Q的二元函数，通过下式表征：
B＝F2(Nge，Q)。
5.根据权利要求1所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，S2中，当机组采用低压缸零出力供热模式时，机组以热定电运行，供热负荷Q是电负荷Nge的一元线性函数，通过下式表征：
Q＝f3(Nge)
标煤消耗量B是电负荷Nge和供热负荷Q的二元函数，通过下式表征：
B＝F3(Nge)
锅炉设计出力Dms0条件下，机组采用低压缸零出力供热模式，电负荷运行范围为：
0≤Nge≤Nge，mas＝f4(Dms0)。
6.根据权利要求1所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，S3的具体优化方法如下：
S301，输入边界参数，边界参数包括供热回水温度t0、热网循环水流量mcw、总供热负荷Qt以及一号机组发电负荷Nge1、二号机组发电负荷Nge2和三号机组发电负荷Nge3；
S302，制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S303，计算Q1＝f1(Nge,Pc0,t0,mcw)，其中，Q1为一号机组供热负荷，Nge为实时发电负荷，Pc0为运行背压，t0为进高背压热网凝汽器的热网循环水温度，mcw为进出高背压热网凝汽器的热网循环水流量；
S304，针对二号机组，判断是否具备投运中排抽汽和低压缸零出力两种模式的电负荷条件，若是，则进入S305，若否，则执行S315；
S305，针对三号机组，判断是否具备投运中排抽汽或低压缸零出力两种模式的电负荷条件，若是，则进入S306，若否，则执行S311和S315；
S306，对二号机组进行热负荷判定，判断二号机组是否具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，若是，则进入S307，若否，则执行S315；
S307，对三号机组进行热负荷判定，判断三号机组是否具备投运中排抽汽和低压缸零出力两种模式的热负荷条件，若是，则进入S308，若否，则执行S315；
S308，进行热负荷判定，判断f2-2(Nge2)+f3-3(Nge3)是否大等于Qt-Q1，若是，则进入S309，若否，则终止寻优迭代流程；
S309，进行热负荷判定，判断f2-3(Nge3)+f3-2(Nge2)是否大于等于Qt-Q1，若是，则进入S310，若否，则终止寻优迭代流程；
S310，二号机组和三号机组均具备中排抽汽或低压缸零出力两种供热模式，进入第一迭代寻优模式，输出最优结果；
S311，判断二号机组是否具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，同时，三号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，则进入S312，若否，则终止寻优迭代流程；
S312，进入第二迭代寻优模式，输出最优结果；
S313，判断二号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，同时，三号机组是否具备投运中排抽汽或低压缸零出力两种模式的电负荷和热负荷条件，若是，则进入S314，若否，则终止寻优迭代流程；
S314，进行第三寻优模式，输出最优结果；
S315，判断二号机组和三号机组是否具备投运中排抽汽模式的电负荷和热负荷条件，若是，则进入S316，若否，则终止寻优迭代流程；
S316，进行第四迭代寻优模式，输出最优结果。
7.根据权利要求6所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，第一迭代寻优模式包括以下步骤：
S10-1-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升5℃对应背压Pc0，作为迭代寻优的背压基准；
S10-1-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-1-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S10-1-4：判断Q30是否大于0，若是，则进入S10-1-5；若否，则转入S10-1-11；
S10-1-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))，进入S10-1-7；
S10-1-7：判断Q21是否小于0，或Q31是否大于等于f3-3(Nge3)，若是，则寻优迭代流程终止；若否，则转入S10-1-8；
S10-1-8：判断Q31是否小于等于f2-3(Nge3)，若是，则进入S10-1-9；若否，则转入S10-1-15；
S10-1-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-1-6；
S10-1-11：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S10-1-12；若否，则转入S10-1-14；
S10-1-12：令Q20＝Qt-Q10，则Q30＝0，进入S10-1-13；
S10-1-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-1-6；
S10-1-14：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S10-1-13；
S10-1-15：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F3-3(Nge3)、B1＝B10+B21+B31，进入S10-1-10；
S10-1-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组的热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-1-17；
S10-1-17：令Pc1＝Pc0+1kPa，进入S10-1-18；
S10-1-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-1-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-1-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b；
S10-2-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S10-2-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S10-2-3：令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S10-2-4：判断Q20是否大于等于0，若是，则进入S10-2-5；若否，则转入S10-2-11；
S10-2-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))，进入S10-2-7；
S10-2-7：判断Q31是否小于等于0，或Q21大于等于f3-2(Nge2)，若是，则终止寻优迭代流程；若否，则转入S10-2-8；
S10-2-8：判断Q21是否小于等于f2-2(Nge2)，若是，则进入S10-2-9；若否，则转入S10-2-15；
S10-2-9：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-10：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S10-2-6；
S10-2-11：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S10-2-12；若否，则转入S10-2-14；
S10-2-12：令Q30＝Qt-Q10，则Q20＝0，进入S10-2-13；
S10-2-13：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S10-2-6；
S10-2-14：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f2-3(Nge3)，进入S10-2-13；
S10-2-15：计算：B10＝F1(Nge,Pc0)、B21＝F3-2(Nge2、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S10-2-10；
S10-2-16：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S10-2-17；
S10-2-17：令Pc1＝Pc0+1kPa，进入S10-2-18；
S10-2-18：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S10-2-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S10-2-19：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
8.根据权利要求6所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，第二寻优模式的具体方法如下：
S12-1：制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S12-2：计算Q10＝f1(Nge,Pc0,t0,mcw)；
S12-3：令Q20＝f3-2(Nge2)，则Q30＝Qt-Q10-f3-2(Nge2)；
S12-4：判断Q30是否大于等于0，若是，则进入S12-5；若否，则转入S12-10；
S12-5：计算：B10＝F1(Nge,Pc0)、B20＝F3-2(Nge2)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S12-6：令Q21＝Q20-(f3-2(Nge2)-f2-2(Nge2))，则Q31＝Q30+(f3-2(Nge2)-f2-2(Nge2))；
S12-7：判断Q21是否小于等于0，或Q31大于等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S12-8；
S12-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S12-9；
S12-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S12-6；
S12-10：判定f2-2(Nge2)是否大于等于Qt-Q10，若是，则进入S12-11；若否，则转入S12-13；
S12-11：令Q20＝Qt-Q10，则Q30＝0，进入S12-12；
S12-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S12-6；
S12-13：令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)，进入S12-12；
S12-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S12-15；
S12-15：令Pc1＝Pc0+1kPa，进入S12-16；
S12-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S12-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S12-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
9.根据权利要求6所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，第三寻优模式的具体方法如下：
S14-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S14-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S14-3令Q30＝f3-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)；
S14-4判断Q20是否大于等于0，若是，则进入S14-5；若否，则转入S14-10；
S14-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F3-3(Nge3)、B0＝B10+B20+B30；
S14-6：令Q31＝Q30-(f3-3(Nge3)-f2-3(Nge3))，则Q21＝Q20+(f3-3(Nge3)-f2-3(Nge3))；
S14-7：判断Q31是否小于等于0，或Q21是否大于等于f2-2(Nge2)，若是，则终止寻优迭代流程；若否，则进入S14-8；
S14-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S14-9；
S14-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S14-6；
S14-10：判定f2-3(Nge3)是否大于等于Qt-Q10，若是，则进入S14-11；若否，则转入S14-13；
S14-11：令Q30＝Qt-Q10，则Q20＝0，进入S14-12；
S14-12：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S14-6；
S14-13：令Q30＝f2-3(Nge3)，则Q20＝Qt-Q10-f3-3(Nge3)，进入S14-12；
S14-14：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S14-15；
S14-15：令Pc1＝Pc0+1kPa，进入S14-16；
S14-16：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S14-2；若否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S14-17：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
10.根据权利要求6所述的一种多机组、多供热模式的热电厂供热负荷分配优化方法，其特征在于，第四寻优模式的具体方法如下：
S16-1制定高背压供热机组运行背压Pc0初始值，以高背压热网凝汽器温升预设阈值对应背压Pc0，作为迭代寻优的背压基准；
S16-2计算Q10＝f1(Nge,Pc0,t0,mcw)；
S16-3令Q20＝f2-2(Nge2)，则Q30＝Qt-Q10-f2-2(Nge2)；
S16-4判断Q30是否大于等于0，若是，则进入S16-5；若否，则转入S16-10；
S16-5：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30；
S16-6：令Q21＝Q20-0.1×Q20则Q31＝Q30+0.1×Q20；
S16-7：判断Q21是否等于0，或Q31是否等于f2-3(Nge3)，若是，则终止寻优迭代流程；若否，则进入S16-8；
S16-8：计算：B10＝F1(Nge,Pc0)、B21＝F2-2(Nge2,Q21)、B31＝F2-3(Nge3,Q31)、B1＝B10+B21+B31，进入S16-9；
S16-9：判断B1是否小于等于B0，若否，则认定原基准工况仍为基准工况；若是，则认定新工况为最优工况，令B1＝Bb0，转入S16-6；
S16-10：令Q20＝Qt-Q10、Q30＝0，进入S16-11；
S16-11：计算：B10＝F1(Nge,Pc0)、B20＝F2-2(Nge2,Q20)、B30＝F2-3(Nge3,Q30)、B0＝B10+B20+B30，进入S16-6；
S12-12：以Pc0背压下全厂总耗煤量最低值Bb0，得出二号机组和三号机组热负荷分配Q2、Q3，作为下一次背压迭代寻优的对比基准，进入S16-13；
S16-13：令Pc1＝Pc0+1kPa，进入S16-14；
S16-14：判断Pc1是否小于等于min(Pc,s、Pc,max)，若是，则以1kPa的幅度提升一号机组背压，并进入S16-2；否，则终止寻优迭代流程，其中Pc,s和Pc,max分别为高背压机组安全运行最高背压值和背压调节手段能达到的实际最高值；
S16-15：以全厂总耗煤量最低值Bb为寻优目标，得出一号机组运行背压Pcb、二号机组和三号机组热负荷分配Q2b、Q3b。
CN202110831724.0A 2021-07-22 2021-07-22 多机组、多供热模式的热电厂供热负荷分配优化方法 Active CN113361828B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110831724.0A CN113361828B (zh) 2021-07-22 2021-07-22 多机组、多供热模式的热电厂供热负荷分配优化方法
PCT/CN2022/101319 WO2023000927A1 (zh) 2021-07-22 2022-06-24 多机组、多供热模式的热电厂供热负荷分配优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110831724.0A CN113361828B (zh) 2021-07-22 2021-07-22 多机组、多供热模式的热电厂供热负荷分配优化方法

Publications (2)

Publication Number Publication Date
CN113361828A true CN113361828A (zh) 2021-09-07
CN113361828B CN113361828B (zh) 2023-04-07

Family

ID=77540143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110831724.0A Active CN113361828B (zh) 2021-07-22 2021-07-22 多机组、多供热模式的热电厂供热负荷分配优化方法

Country Status (2)

Country Link
CN (1) CN113361828B (zh)
WO (1) WO2023000927A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114548659A (zh) * 2022-01-07 2022-05-27 华能国际电力股份有限公司大连电厂 一种厂级多类型供热机组电、热负荷分配方法和设备
WO2023000927A1 (zh) * 2021-07-22 2023-01-26 西安热工研究院有限公司 多机组、多供热模式的热电厂供热负荷分配优化方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116011789B (zh) * 2023-02-07 2023-10-31 华能济南黄台发电有限公司 一种城市电厂供热机组热水蓄能系统
CN116341852B (zh) * 2023-03-27 2024-04-26 湖北清江水电开发有限责任公司 一种用于水电厂的多机组负荷分配方法
CN116663870B (zh) * 2023-08-02 2023-10-03 北京世纪黄龙技术有限公司 基于云计算的供热系统调度方法及系统
CN117078047B (zh) * 2023-10-16 2024-02-23 华能济南黄台发电有限公司 一种基于lstm的热负荷预测与分配优化方法及系统
CN117787653B (zh) * 2024-01-22 2024-06-04 上海金联热电有限公司 母管制热电联产机组负荷分配方法、系统、介质及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2865176A1 (en) * 2012-03-07 2013-09-12 Siemens Corporation Online heurisitc algorithm for combined cooling heating and power plant optimization
CN111401771A (zh) * 2020-03-25 2020-07-10 西安热工研究院有限公司 一种热电联产机组多供热方式的经济调度方法
CN112633560A (zh) * 2020-12-11 2021-04-09 西安交通大学 一种含燃煤热电联产机组的电站优化调度方法
CN112964492A (zh) * 2021-03-22 2021-06-15 西安热工研究院有限公司 一种适用于高背压梯级供热机组的供热耗煤在线测定方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009150121A1 (de) * 2008-06-10 2009-12-17 Siemens Aktiengesellschaft Verfahren zur optimierten kraftwerkseinsatzplanung
CN111287811B (zh) * 2020-03-31 2021-06-15 西安热工研究院有限公司 一种高背压梯级供热机组最佳运行真空在线寻优方法
CN112906292B (zh) * 2021-01-26 2024-02-23 西安热工研究院有限公司 热电联产机组厂级热电负荷在线优化分配的方法、系统、设备及存储介质
CN113361828B (zh) * 2021-07-22 2023-04-07 西安热工研究院有限公司 多机组、多供热模式的热电厂供热负荷分配优化方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2865176A1 (en) * 2012-03-07 2013-09-12 Siemens Corporation Online heurisitc algorithm for combined cooling heating and power plant optimization
CN111401771A (zh) * 2020-03-25 2020-07-10 西安热工研究院有限公司 一种热电联产机组多供热方式的经济调度方法
CN112633560A (zh) * 2020-12-11 2021-04-09 西安交通大学 一种含燃煤热电联产机组的电站优化调度方法
CN112964492A (zh) * 2021-03-22 2021-06-15 西安热工研究院有限公司 一种适用于高背压梯级供热机组的供热耗煤在线测定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DONG ANYOU ETC.: "Research on Heat/Power Load Distribution and Coordinated Control of Cogeneration Units", 《2018 INTERNATIONAL CONFERENCE ON ENGINEERING SIMULATION AND INTELLIGENT CONTROL (ESAIC)》 *
吕凯等: "330MW高背压供热机组热力特性研究", 《汽轮机技术》 *
邹淼: "区域热电厂负荷优化分配与热网优化研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023000927A1 (zh) * 2021-07-22 2023-01-26 西安热工研究院有限公司 多机组、多供热模式的热电厂供热负荷分配优化方法
CN114548659A (zh) * 2022-01-07 2022-05-27 华能国际电力股份有限公司大连电厂 一种厂级多类型供热机组电、热负荷分配方法和设备

Also Published As

Publication number Publication date
WO2023000927A1 (zh) 2023-01-26
CN113361828B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN113361828B (zh) 多机组、多供热模式的热电厂供热负荷分配优化方法
CN108316980B (zh) 一种火电机组熔盐蓄热放热调峰系统
CN109118017B (zh) 热负荷优化分配方法、电子设备和存储介质
CN113255160B (zh) 一种直接空冷机组低真空供热运行背压寻优系统及方法
CN113343490B (zh) 一种耦合熔盐储热的工业供汽电站运行优化方法及系统
CN111706411A (zh) 一种背压机组改造为抽凝机组的热力系统及工作方法
CN113565591B (zh) 热电联产机组耦合空气储能膨胀发电系统及优化运行方法
CN110005488B (zh) 一种高背压供热系统节能优化方法
CN110288135A (zh) 一种高压加热系统疏水水位节能优化方法
CN111520204A (zh) 一种提高再热式汽轮机供汽能力及机组效率的系统和方法
CN110991877A (zh) 供热机组采用低压缸切缸灵活性改造后供热及调峰能力改善评估方法
CN113283121A (zh) 一种熔盐储热工业供汽系统的流程及容量设计方法及系统
CN113898428B (zh) 分析双抽热电联产机组热、电、汽负荷调整空间的方法
CN112070358A (zh) 低真空供热机组电负荷调整区间确定方法及系统
CN112964492B (zh) 一种适用于高背压梯级供热机组的供热耗煤在线测定方法
CN204476487U (zh) 中间抽汽采暖加热用后置式汽轮机系统
CN113325713B (zh) 采用匹配抽汽外供技术的供热机组最佳运行方式确定方法
CN111463773A (zh) 一种区域型综合能源系统能量管理优化方法和装置
CN114462163B (zh) 一种考虑能源品质的㶲集线器模型
CN116253333A (zh) 一种氨合成系统及氨合成控制方法
CN212406830U (zh) 一种背压机组改造为抽凝机组的热力系统
CN112700350A (zh) 一种热电联产机组供热影响供电煤耗计算方法及装置
CN114459069B (zh) 一种基于网源协调的一次管网最佳供水温度的控制方法
Van Riet et al. Hydronic Optimization Of Hybrid Heating Systems: A Methodology Based On Base Circuits
CN115095898A (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