CN112288219A - 一种基于㶲损率分析的园区热力系统能效提升方法 - Google Patents

一种基于㶲损率分析的园区热力系统能效提升方法 Download PDF

Info

Publication number
CN112288219A
CN112288219A CN202011000200.9A CN202011000200A CN112288219A CN 112288219 A CN112288219 A CN 112288219A CN 202011000200 A CN202011000200 A CN 202011000200A CN 112288219 A CN112288219 A CN 112288219A
Authority
CN
China
Prior art keywords
loss
load
rate
total
subsystem
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
CN202011000200.9A
Other languages
English (en)
Other versions
CN112288219B (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.)
Tianjin University
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd
Original Assignee
Tianjin University
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power 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 Tianjin University, State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd, Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd filed Critical Tianjin University
Priority to CN202011000200.9A priority Critical patent/CN112288219B/zh
Publication of CN112288219A publication Critical patent/CN112288219A/zh
Application granted granted Critical
Publication of CN112288219B publication Critical patent/CN112288219B/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/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/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Computational Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Mathematical Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Game Theory and Decision Science (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Primary Health Care (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种基于
Figure DDA0002694025230000011
损率分析的园区热力系统能效提升方法，包括如下步骤：步骤1、对热力系统进行热负荷分配；步骤2、建立各子系统
Figure DDA0002694025230000012
损失及
Figure DDA0002694025230000013
损率模型，计算各子系统的
Figure DDA0002694025230000014
损率、
Figure DDA0002694025230000015
损敏感性，最后计算系统总
Figure DDA0002694025230000016
损率；步骤3、根据
Figure DDA0002694025230000017
损敏感性调节负荷分配，根据前一日的运行数据可得改变系统
Figure DDA0002694025230000018
损率对应的负荷调节方向，根据实际设备情况设定调节步长ΔD，每次调节时优先对系统承载能力与总系统负荷需求进行判断，判断为合理后对系统是否需要电加热补偿进行判断，之后进行系统总
Figure DDA0002694025230000019
损率变化的计算；步骤4、重复步骤2和步骤3，直至系统总
Figure DDA00026940252300000110
损率达到最小。本发明通过分配调节各系统的供热负荷，实现整个系统在满足负荷需求下的最小
Figure DDA00026940252300000111
损率运行。

Description

一种基于㶲损率分析的园区热力系统能效提升方法
技术领域
本发明属于多种能源综合利用的供能技术领域，涉及园区热力系统能效提升方法，尤其是一种基于
Figure BDA0002694025210000012
损率分析的园区热力系统能效提升方法。
背景技术
我国能源生产量和消费量均已居世界前列，但在能源供给和利用方式上存在着一系列突出问题，如能源结构不合理、能源利用效率不高、可再生能源开发利用比例低、能源安全利用水平有待进一步提高等。针对传统能源布局中在智慧城市建设增速和综合能源交易市场发展的背景下，综合能源服务成为各个能源服务公司争夺终端用户的主要竞争手段。在此背景下，综合能源系统应运而生，其理念突破了传统能源系统的技术、市场和管理壁垒，是电、气、热、冷等各类能源统一规划、统一调度的综合性能源系统，对于推动能源结构转型，推进我国能源革命具有重要意义。综合能源系统势必将成为未来能源系统的主要形态。
目前传统的能源系统的冷、热、电等不同供能子系统一般是相互独立的，其功能结构设计与运行控制往往不相关联，使得整个能源系统很难实现最大化的能源利用以及最小化的能源损失。采用多能互补的综合统筹能效分析能够将冷、热、电等各供能子系统进行配合、协调以及优化，从各系统之间的联系入手考虑整个能源总系统的综合性能。
目前传统的能源系统，包括多能互补的综合系统多是以能源利用率的角度进行评估，从热力系统角度来看，即多是从热效率的角度进行评估。
Figure BDA0002694025210000013
正是稳定流动的工质，遵循卡诺循环所能做出的最大机械功。或者说，在实质上是表示工质在某一状态下可作的最大机械功。虽然迄今为止完全按卡诺循环进行工作的热机和工质还不存在，但是
Figure BDA0002694025210000014
却有着理论和实践上的重要意义。
Figure BDA0002694025210000021
是理论上工质从某一状态变化到另一状态所能做出的最大机械功。因此
Figure BDA0002694025210000022
可以作为衡量工质质量或者品位的量度，从而使我们根据工质
Figure BDA0002694025210000023
值的大小和
Figure BDA0002694025210000024
损的大小，来确定工质品位的高低和鉴定热力设备(过程)热力学的完善程度，从而指导节能工作。
发明内容
本发明的目的在于克服现有技术的不足，提供一种基于
Figure BDA00026940252100000216
损率分析的园区热力系统能效提升方法，能够实现整个综合供能系统
Figure BDA0002694025210000025
损率的最小化。
本发明解决其技术问题是采取以下技术方案实现的：
一种基于
Figure BDA0002694025210000026
损率分析的园区热力系统能效提升方法，包括如下步骤：
步骤1、结合约束条件，根据前一日的热负荷分配比例对热力系统进行热负荷分配；
步骤2、建立各子系统
Figure BDA0002694025210000027
损失及
Figure BDA0002694025210000028
损率模型，根据步骤1的热负荷分配结果，采集地源热泵、燃气锅炉、余热锅炉和电加热的各子系统的运行数据，并根据采集到的数据计算各子系统的
Figure BDA0002694025210000029
损率、
Figure BDA00026940252100000210
损敏感性，最后计算系统总
Figure BDA00026940252100000211
损率；
步骤3、在保证总负荷满足的前提下，首先根据
Figure BDA00026940252100000212
损敏感性调节负荷分配，根据前一日的运行数据可得改变系统
Figure BDA00026940252100000213
损率对应的负荷调节方向，根据实际设备情况设定调节步长ΔD，每次调节时优先对系统承载能力与总系统负荷需求进行判断，判断为合理后对系统是否需要电加热补偿进行判断，之后进行系统总
Figure BDA00026940252100000214
损率变化的计算；
步骤4、重复步骤2和步骤3，直至系统总
Figure BDA00026940252100000215
损率达到最小。
而且，所述步骤1的约束条件包括：
(1)子系统设备负荷承载能力
客观存在的约束，每个子系统设备具有负荷承载上限，安排设备在承载能力之外运行是不现实的。
(2)用户热负荷总需求
首先需要满足的最基本约束，根据客观的子系统负荷承载能力合理设置调配使总系统负荷承载能力满足用户热负荷需求。
(3)辅助系统
优先使用地源热泵、燃气锅炉、燃气轮机余热锅炉，以下情况时使用电加热设备进行补偿：各子系统满负荷后系统总负荷仍不能满足需求符合时；为满足总负荷需求而增大某子系统负荷时产生的
Figure BDA0002694025210000032
损率增加大于直接用电加热补偿的
Figure BDA0002694025210000033
损率增加时。
而且，所述步骤2的具体步骤包括：
(1)建立各子系统
Figure BDA0002694025210000034
损失及
Figure BDA0002694025210000035
损率模型：
①燃气锅炉
Figure BDA0002694025210000036
损失
Lg＝Lgout+Lgin，
其中外部
Figure BDA0002694025210000037
损失Lgout包括排烟
Figure BDA0002694025210000038
损失Lga，气体未完全燃烧
Figure BDA00026940252100000316
损失Lgb，散热
Figure BDA00026940252100000310
损失Lgc，灰渣物理
Figure BDA0002694025210000039
损失Lgd，即
Lgout＝Lga+Lgb+Lgc+Lgd，
内部
Figure BDA00026940252100000311
损失Lgin包括燃烧不可逆
Figure BDA00026940252100000312
损失Lge，传热
Figure BDA00026940252100000313
损失Lgf，即
Lgin＝Lge+Lgf，
燃气锅炉总
Figure BDA00026940252100000314
损失
Lg＝Lga+Lgb+Lgc+Lgd+Lge+Lgf，
燃气锅炉总
Figure BDA00026940252100000315
损率
Figure BDA0002694025210000031
其中，Loadg表示燃气锅炉承载的热负荷。
②余热锅炉
Figure BDA0002694025210000044
损失
Lw＝Lwout+Lwin，
其中外部
Figure BDA0002694025210000045
损失Lwout包括排烟
Figure BDA0002694025210000046
损失Lwa，散热
Figure BDA0002694025210000047
损失Lwb，灰渣物理
Figure BDA0002694025210000048
损失Lwc，即
Lwout＝Lwa+Lwb+Lwc，
内部
Figure BDA0002694025210000049
损失Lwin为传热
Figure BDA00026940252100000410
损失Lwd，即
Lwin＝Lwd，
余热锅炉总
Figure BDA00026940252100000411
损失
Lw＝Lwa+Lwb+Lwc+Lwd，
余热锅炉总
Figure BDA00026940252100000412
损率
Figure BDA0002694025210000041
其中，Loadw表示余热锅炉承载的热负荷。
③地源热泵
Figure BDA00026940252100000413
损失
地源热泵系统
Figure BDA00026940252100000414
损失包括地埋管换热器
Figure BDA00026940252100000416
损失Lpa、空气处理单元
Figure BDA00026940252100000417
损失Lpb、热泵机组
Figure BDA00026940252100000415
损失Lpc、系统中各水泵的总能耗∑Wpump、系统中末端风机的总能耗∑Wfan，即
Lp＝Lpa+Lpb+Lpc+∑Wpump+∑Wfan，
地源热泵总
Figure BDA0002694025210000043
损率
Figure BDA0002694025210000042
其中，Loadp表示地源热泵承载的热负荷。
④直接式电加热设备
Figure BDA0002694025210000057
损失
Figure BDA0002694025210000051
其中，T0、Ts分别为环境温度、设备供暖时的绝对温度，在
Figure BDA0002694025210000058
损失计算时认定其为定值，因而认为电加热设备
Figure BDA0002694025210000059
损失只与电加热设备提供热负荷Loaded有关系，因此，电加热设备的总
Figure BDA00026940252100000510
损率为
Figure BDA0002694025210000052
(2)计算
Figure BDA00026940252100000511
损敏感性
定义系统
Figure BDA00026940252100000512
损敏感性，为
Figure BDA00026940252100000513
损率的变化量与负荷变化量的比值，即
Figure BDA0002694025210000053
针对各子系统：优先调节运行于
Figure BDA00026940252100000514
损率较高工况的子系统，当多个子系统运行于
Figure BDA00026940252100000515
损率较高工况时，应优先调节
Figure BDA00026940252100000516
损敏感性强的子系统；当多个子系统运行于
Figure BDA00026940252100000517
损率较低工况时，应优先调节
Figure BDA00026940252100000518
损敏感性弱的子系统。
(3)计算系统总
Figure BDA00026940252100000519
损失
Figure BDA0002694025210000054
系统总
Figure BDA00026940252100000520
损率
Figure BDA0002694025210000055
其中，Loadtotal表示系统承载的总热负荷；
而且，所述步骤3的计算系统总
Figure BDA00026940252100000521
损率变化的计算公式为：
Figure BDA0002694025210000056
其中，下标1表示当次优化，下标0表示前一次优化，每日首次优化时，前一次优化的系统总
Figure BDA0002694025210000061
损率定义为0，根据计算结果判断是否需要进行下一次优化。
而且，所述步骤3的根据
Figure BDA0002694025210000062
损敏感性调节负荷分配的原则为：优先调节运行于
Figure BDA0002694025210000063
损率较高工况的子系统，当多个子系统运行于
Figure BDA0002694025210000064
损率较高工况时，应优先调节
Figure BDA0002694025210000065
损敏感性强的子系统；当多个子系统运行于
Figure BDA0002694025210000066
损率较低工况时，应优先调节
Figure BDA0002694025210000067
损敏感性弱的子系统。
而且，所述步骤3的针对系统承载能力与总系统负荷需求的判断原则为：存在某一个或某几个子系统分配负荷超过其承载能力，判断为不合理；总系统负荷不满足用户热负荷需求，判断为不合理；同时保证各子系统分配负荷不超额，且总系统负荷满足用户热负荷需求，判断为合理。
而且，所述步骤3的是否需要电加热补偿的判断原则：各子系统满负荷后系统总负荷仍不能满足需求符合时、为满足总负荷需求而增大某子系统负荷时产生的
Figure BDA0002694025210000068
损率增加大于直接用电加热补偿的
Figure BDA0002694025210000069
损率增加时判断为需要电加热补偿。
而且，所述步骤4的判断系统总
Figure BDA00026940252100000610
损率最小的方法为：ΔELtotal等于0或者趋近0，这表示在该方法下系统总
Figure BDA00026940252100000611
损率已经达到了最小。
本发明的优点和积极效果是：
1、本发明以子系统为地源热泵、燃气锅炉、燃气轮机余热锅炉、电加热的多能互补生态园区供热系统为研究对象，以热负荷分配为手段，以
Figure BDA00026940252100000612
为评价标准，提出了一种系统能效提升方法。该方法从
Figure BDA00026940252100000613
的角度进行分析，从热力学第二定律的高度，不仅考虑能量的数量，也考虑能量的品质，基于
Figure BDA00026940252100000614
损率最小的原则对园区热力系统热负荷分配进行优化，进一步统筹优化各种供能子系统工作的协调配合，通过分配调节各系统的供热负荷，实现整个系统在满足负荷需求下的最小
Figure BDA0002694025210000071
损率运行。
2、本发明能够从
Figure BDA0002694025210000072
损率分析的角度出发，通过不断重新分配园区热力系统中地源热泵、燃气锅炉、余热锅炉、电加热等各子系统的工作中的热负荷，实现整个综合供能系统
Figure BDA0002694025210000073
损率的最小化。
附图说明
图1为本发明的处理流程图。
具体实施方式
以下结合附图对本发明实施例作进一步详述：
一种基于
Figure BDA0002694025210000074
损率分析的园区热力系统能效提升方法，如图1所示，包括如下步骤：
步骤1、结合约束条件，根据前一日的热负荷分配比例对热力系统进行热负荷分配。
在本实施例中，所述步骤1的约束条件包括：
(1)子系统设备负荷承载能力
客观存在的约束，每个子系统设备具有负荷承载上限，安排设备在承载能力之外运行是不现实的。
(2)用户热负荷总需求
首先需要满足的最基本约束，根据客观的子系统负荷承载能力合理设置调配使总系统负荷承载能力满足用户热负荷需求。
(3)辅助系统
优先使用地源热泵、燃气锅炉、燃气轮机余热锅炉，以下情况时使用电加热设备进行补偿：各子系统满负荷后系统总负荷仍不能满足需求符合时；为满足总负荷需求而增大某子系统负荷时产生的
Figure BDA0002694025210000075
损率增加大于直接用电加热补偿的
Figure BDA0002694025210000082
损率增加时。
步骤2、建立各子系统
Figure BDA0002694025210000083
损失及
Figure BDA0002694025210000084
损率模型，根据步骤1的热负荷分配结果，采集地源热泵、燃气锅炉、余热锅炉和电加热的各子系统的运行数据，并根据采集到的数据计算各子系统的
Figure BDA0002694025210000085
损率、
Figure BDA0002694025210000086
损敏感性，最后计算系统总
Figure BDA0002694025210000087
损率ELtotal1。
在本实施例中，所述步骤2的具体步骤包括：
(1)建立各子系统
Figure BDA0002694025210000088
损失及
Figure BDA0002694025210000089
损率模型：
①燃气锅炉
Figure BDA00026940252100000810
损失
Lg＝Lgout+Lgin，
其中外部
Figure BDA00026940252100000811
损失Lgout包括排烟
Figure BDA00026940252100000814
损失Lga，气体未完全燃烧
Figure BDA00026940252100000815
损失Lgb，散热
Figure BDA00026940252100000812
损失Lgc，灰渣物理
Figure BDA00026940252100000813
损失Lgd，即
Lgout＝Lga+Lgb+Lgc+Lgd，
内部
Figure BDA00026940252100000818
损失Lgin包括燃烧不可逆
Figure BDA00026940252100000816
损失Lge，传热
Figure BDA00026940252100000817
损失Lgf，即
Lgin＝Lge+Lgf，
燃气锅炉总
Figure BDA00026940252100000819
损失
Lg＝Lga+Lgb+Lgc+Lgd+Lge+Lgf，
燃气锅炉总
Figure BDA00026940252100000820
损率
Figure BDA0002694025210000081
其中，Loadg表示燃气锅炉承载的热负荷。
②余热锅炉
Figure BDA00026940252100000821
损失
Lw＝Lwout+Lwin，
其中外部
Figure BDA00026940252100000822
损失Lwout包括排烟
Figure BDA00026940252100000823
损失Lwa，散热
Figure BDA00026940252100000824
损失Lwb，灰渣物理
Figure BDA00026940252100000825
损失Lwc，即
Lwout＝Lwa+Lwb+Lwc，
内部
Figure BDA0002694025210000094
损失Lwin为传热
Figure BDA0002694025210000095
损失Lwd，即
Lwin＝Lwd，
余热锅炉总
Figure BDA0002694025210000096
损失
Lw＝Lwa+Lwb+Lwc+Lwd，
余热锅炉总
Figure BDA0002694025210000097
损率
Figure BDA0002694025210000091
其中，Loadw表示余热锅炉承载的热负荷。
③地源热泵
Figure BDA0002694025210000098
损失
地源热泵系统
Figure BDA0002694025210000099
损失包括地埋管换热器
Figure BDA00026940252100000910
损失Lpa、空气处理单元
Figure BDA00026940252100000911
损失Lpb、热泵机组
Figure BDA00026940252100000912
损失Lpc、系统中各水泵的总能耗∑Wpump、系统中末端风机的总能耗∑Wfan，即
Lp＝Lpa+Lpb+Lpc+∑Wpump+∑Wfan，
地源热泵总
Figure BDA00026940252100000913
损率
Figure BDA0002694025210000092
其中，Loadp表示地源热泵承载的热负荷。
④直接式电加热设备
Figure BDA00026940252100000914
损失
Figure BDA0002694025210000093
其中，T0、Ts分别为环境温度、设备供暖时的绝对温度，在
Figure BDA00026940252100000915
损失计算时认定其为定值，因而认为电加热设备
Figure BDA00026940252100000916
损失只与电加热设备提供热负荷Loaded有关系，因此，电加热设备的总
Figure BDA0002694025210000106
损率为
Figure BDA0002694025210000101
(2)计算
Figure BDA0002694025210000107
损敏感性
定义系统
Figure BDA0002694025210000108
损敏感性，为
Figure BDA0002694025210000109
损率的变化量与负荷变化量的比值，即
Figure BDA0002694025210000102
针对各子系统：优先调节运行于
Figure BDA00026940252100001010
损率较高工况的子系统，当多个子系统运行于
Figure BDA00026940252100001011
损率较高工况时，应优先调节
Figure BDA00026940252100001012
损敏感性强的子系统；当多个子系统运行于
Figure BDA00026940252100001013
损率较低工况时，应优先调节
Figure BDA00026940252100001014
损敏感性弱的子系统。
(3)计算系统总
Figure BDA00026940252100001015
损失
Figure BDA0002694025210000103
系统总
Figure BDA00026940252100001016
损率
Figure BDA0002694025210000104
其中，Loadtotal表示系统承载的总热负荷；
步骤3、在保证总负荷满足的前提下，首先根据
Figure BDA00026940252100001017
损敏感性调节负荷分配，根据前一日的运行数据可得改变系统
Figure BDA00026940252100001018
损率对应的负荷调节方向，根据实际设备情况设定调节步长ΔD，每次调节时优先对系统承载能力与总系统负荷需求进行判断，判断为合理后对系统是否需要电加热补偿进行判断，之后进行系统总
Figure BDA00026940252100001019
损率变化的计算
Figure BDA0002694025210000105
其中，下标1表示当次优化，下标0表示前一次优化，每日首次优化时，前一次优化的系统总
Figure BDA0002694025210000113
损率定义为0，根据计算结果判断是否需要进行下一次优化。
在本实施例中，所述步骤3的根据
Figure BDA0002694025210000114
损敏感性调节负荷分配的原则为：优先调节运行于
Figure BDA0002694025210000116
损率较高工况的子系统，当多个子系统运行于
Figure BDA0002694025210000115
损率较高工况时，应优先调节
Figure BDA0002694025210000117
损敏感性强的子系统；当多个子系统运行于
Figure BDA0002694025210000119
损率较低工况时，应优先调节
Figure BDA0002694025210000118
损敏感性弱的子系统。
在本实施例中，所述步骤3的针对系统承载能力与总系统负荷需求的判断原则为：存在某一个或某几个子系统分配负荷超过其承载能力，判断为不合理；总系统负荷不满足用户热负荷需求，判断为不合理；同时保证各子系统分配负荷不超额，且总系统负荷满足用户热负荷需求，判断为合理。
在本实施例中，所述步骤3的是否需要电加热补偿的判断原则：各子系统满负荷后系统总负荷仍不能满足需求符合时、为满足总负荷需求而增大某子系统负荷时产生的
Figure BDA00026940252100001110
损率增加大于直接用电加热补偿的
Figure BDA00026940252100001111
损率增加时判断为需要电加热补偿。
步骤4、重复步骤2和步骤3，直至系统总
Figure BDA00026940252100001113
损率
Figure BDA0002694025210000111
达到最小，判断系统总
Figure BDA00026940252100001112
损率最小的方法为：
Figure BDA0002694025210000112
等于0或者趋近0，这表示在该方法下系统总
Figure BDA00026940252100001114
损率已经达到了最小。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。

Claims (7)

1.一种基于
Figure FDA0002694025200000011
损率分析的园区热力系统能效提升方法，其特征在于：包括如下步骤：
步骤1、结合约束条件，根据前一日的热负荷分配比例对热力系统进行热负荷分配；
步骤2、建立各子系统
Figure FDA0002694025200000012
损失及
Figure FDA0002694025200000013
损率模型，根据步骤1的热负荷分配结果，采集地源热泵、燃气锅炉、余热锅炉和电加热的各子系统的运行数据，并根据采集到的数据计算各子系统的
Figure FDA0002694025200000014
损率、
Figure FDA0002694025200000015
损敏感性，最后计算系统总
Figure FDA00026940252000000111
损率；
步骤3、在保证总负荷满足的前提下，首先根据
Figure FDA0002694025200000016
损敏感性调节负荷分配，根据前一日的运行数据可得改变系统
Figure FDA0002694025200000017
损率对应的负荷调节方向，根据实际设备情况设定调节步长ΔD，每次调节时优先对系统承载能力与总系统负荷需求进行判断，判断为合理后对系统是否需要电加热补偿进行判断，之后进行系统总
Figure FDA00026940252000000110
损率变化的计算；
步骤4、重复步骤2和步骤3，直至系统总
Figure FDA0002694025200000018
损率达到最小。
2.根据权利要求1所述的一种基于
Figure FDA0002694025200000019
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤1的约束条件包括：
(1)子系统设备负荷承载能力
客观存在的约束，每个子系统设备具有负荷承载上限，安排设备在承载能力之外运行是不现实的。
(2)用户热负荷总需求
首先需要满足的最基本约束，根据客观的子系统负荷承载能力合理设置调配使总系统负荷承载能力满足用户热负荷需求。
(3)辅助系统
优先使用地源热泵、燃气锅炉、燃气轮机余热锅炉，以下情况时使用电加热设备进行补偿：各子系统满负荷后系统总负荷仍不能满足需求符合时；为满足总负荷需求而增大某子系统负荷时产生的
Figure FDA0002694025200000022
损率增加大于直接用电加热补偿的
Figure FDA0002694025200000027
损率增加时。
3.根据权利要求1所述的一种基于
Figure FDA0002694025200000023
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤2的具体步骤包括：
(1)建立各子系统
Figure FDA0002694025200000026
损失及
Figure FDA0002694025200000024
损率模型：
①燃气锅炉
Figure FDA0002694025200000025
损失
Lg＝Lgout+Lgin，
其中外部
Figure FDA0002694025200000028
损失Lgout包括排烟
Figure FDA00026940252000000211
损失Lga，气体未完全燃烧
Figure FDA00026940252000000216
损失Lgb，散热
Figure FDA0002694025200000029
损失Lgc，灰渣物理
Figure FDA00026940252000000210
损失Lgd，即
Lgout＝Lga+Lgb+Lgc+Lgd，
内部
Figure FDA00026940252000000217
损失Lgin包括燃烧不可逆
Figure FDA00026940252000000212
损失Lge，传热
Figure FDA00026940252000000213
损失Lgf，即
Lgin＝Lge+Lgf，
燃气锅炉总
Figure FDA00026940252000000214
损失
Lg＝Lga+Lgb+Lgc+Lgd+Lge+Lgi，
燃气锅炉总
Figure FDA00026940252000000215
损率
Figure FDA0002694025200000021
其中，Loadg表示燃气锅炉承载的热负荷。
②余热锅炉
Figure FDA0002694025200000033
损失
Lw＝Lwout+Lwin，
其中外部
Figure FDA0002694025200000034
损失Lwout包括排烟
Figure FDA0002694025200000036
损失Lwa，散热
Figure FDA0002694025200000037
损失Lwb，灰渣物理
Figure FDA00026940252000000317
损失Lwc，即
Lwout＝Lwa+Lwb+Lwc，
内部
Figure FDA0002694025200000035
损失Lwin为传热
Figure FDA0002694025200000038
损失Lwd，即
Lwin＝Lwd，
余热锅炉总
Figure FDA0002694025200000039
损失
Lw＝Lwa+Lwb+Lwc+Lwd，
余热锅炉总
Figure FDA00026940252000000310
损率
Figure FDA0002694025200000031
其中，Loadw表示余热锅炉承载的热负荷。
③地源热泵
Figure FDA00026940252000000311
损失
地源热泵系统
Figure FDA00026940252000000312
损失包括地埋管换热器
Figure FDA00026940252000000316
损失Lpa、空气处理单元
Figure FDA00026940252000000318
损失Lpb、热泵机组
Figure FDA00026940252000000313
损失Lpc、系统中各水泵的总能耗∑Wpump、系统中末端风机的总能耗∑Wfan，即
Lp＝Lpa+Lpb+Lpc+∑Wpump+∑Wfan，
地源热泵总
Figure FDA00026940252000000314
损率
Figure FDA0002694025200000032
其中，Loadp表示地源热泵承载的热负荷。
④直接式电加热设备
Figure FDA00026940252000000315
损失
Figure FDA0002694025200000041
其中，T0、Ts分别为环境温度、设备供暖时的绝对温度，在
Figure FDA00026940252000000420
损失计算时认定其为定值，因而认为电加热设备
Figure FDA00026940252000000419
损失只与电加热设备提供热负荷Loaded有关系，因此，电加热设备的总
Figure FDA00026940252000000418
损率为
Figure FDA0002694025200000042
(2)计算
Figure FDA00026940252000000416
损敏感性
定义系统
Figure FDA00026940252000000417
损敏感性，为
Figure FDA00026940252000000415
损率的变化量与负荷变化量的比值，即
Figure FDA0002694025200000043
针对各子系统：优先调节运行于
Figure FDA00026940252000000413
损率较高工况的子系统，当多个子系统运行于
Figure FDA00026940252000000411
损率较高工况时，应优先调节
Figure FDA00026940252000000412
损敏感性强的子系统；当多个子系统运行于
Figure FDA00026940252000000410
损率较低工况时，应优先调节
Figure FDA00026940252000000414
损敏感性弱的子系统。
(3)计算系统总
Figure FDA0002694025200000049
损失
Figure FDA0002694025200000044
系统总
Figure FDA0002694025200000048
损率
Figure FDA0002694025200000046
其中，Loadtotal表示系统承载的总热负荷；
而且，所述步骤3的计算系统总
Figure FDA0002694025200000047
损率变化的计算公式为：
Figure FDA0002694025200000045
其中，下标1表示当次优化，下标0表示前一次优化，每日首次优化时，前一次优化的系统总
Figure FDA0002694025200000051
损率定义为0，根据计算结果判断是否需要进行下一次优化。
4.根据权利要求1所述的一种基于
Figure FDA0002694025200000052
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤3的根据
Figure FDA0002694025200000053
损敏感性调节负荷分配的原则为：优先调节运行于
Figure FDA0002694025200000054
损率较高工况的子系统，当多个子系统运行于
Figure FDA00026940252000000514
损率较高工况时，应优先调节
Figure FDA0002694025200000055
损敏感性强的子系统；当多个子系统运行于
Figure FDA00026940252000000515
损率较低工况时，应优先调节
Figure FDA0002694025200000056
损敏感性弱的子系统。
5.根据权利要求1所述的一种基于
Figure FDA0002694025200000057
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤3的针对系统承载能力与总系统负荷需求的判断原则为：存在某一个或某几个子系统分配负荷超过其承载能力，判断为不合理；总系统负荷不满足用户热负荷需求，判断为不合理；同时保证各子系统分配负荷不超额，且总系统负荷满足用户热负荷需求，判断为合理。
6.根据权利要求1所述的一种基于
Figure FDA0002694025200000058
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤3的是否需要电加热补偿的判断原则：各子系统满负荷后系统总负荷仍不能满足需求符合时、为满足总负荷需求而增大某子系统负荷时产生的
Figure FDA0002694025200000059
损率增加大于直接用电加热补偿的
Figure FDA00026940252000000513
损率增加时判断为需要电加热补偿。
7.根据权利要求1所述的一种基于
Figure FDA00026940252000000510
损率分析的园区热力系统能效提升方法，其特征在于：所述步骤4的判断系统总
Figure FDA00026940252000000511
损率最小的方法为：
Figure FDA00026940252000000516
等于0或者趋近0，这表示在该方法下系统总
Figure FDA00026940252000000512
损率已经达到了最小。
CN202011000200.9A 2020-09-22 2020-09-22 一种园区热力系统能效提升方法 Active CN112288219B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011000200.9A CN112288219B (zh) 2020-09-22 2020-09-22 一种园区热力系统能效提升方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011000200.9A CN112288219B (zh) 2020-09-22 2020-09-22 一种园区热力系统能效提升方法

Publications (2)

Publication Number Publication Date
CN112288219A true CN112288219A (zh) 2021-01-29
CN112288219B CN112288219B (zh) 2023-04-07

Family

ID=74422235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011000200.9A Active CN112288219B (zh) 2020-09-22 2020-09-22 一种园区热力系统能效提升方法

Country Status (1)

Country Link
CN (1) CN112288219B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109118017A (zh) * 2018-09-03 2019-01-01 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 热负荷优化分配方法、电子设备和存储介质
CN109858759A (zh) * 2018-12-29 2019-06-07 陕西鼓风机（集团）有限公司 一种工业园区综合能源平衡调度方法
CN110070460A (zh) * 2019-04-16 2019-07-30 东南大学 多套燃气-蒸汽联合循环机组热电负荷优化调度系统
CN110276489A (zh) * 2019-06-21 2019-09-24 华北电力大学 基于*损最低的多能量枢纽分布式优化调度方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109118017A (zh) * 2018-09-03 2019-01-01 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 热负荷优化分配方法、电子设备和存储介质
CN109858759A (zh) * 2018-12-29 2019-06-07 陕西鼓风机（集团）有限公司 一种工业园区综合能源平衡调度方法
CN110070460A (zh) * 2019-04-16 2019-07-30 东南大学 多套燃气-蒸汽联合循环机组热电负荷优化调度系统
CN110276489A (zh) * 2019-06-21 2019-09-24 华北电力大学 基于*损最低的多能量枢纽分布式优化调度方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘升男等: "基于㶲分析的联合运行热泵机组负荷优化分配研究", 《建筑科学》 *
赵兴杰: "基于热力学第二定律的火电厂热力系统节能研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技II辑》 *

Also Published As

Publication number Publication date
CN112288219B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN109858759B (zh) 一种工业园区综合能源平衡调度方法
CN101667013B (zh) 微型燃气轮机冷电联供分布式供能系统优化运行控制方法
Wang et al. Analysis of inlet air throttling operation method for gas turbine in performance of CCHP system under different operation strategies
CN110070460A (zh) 多套燃气-蒸汽联合循环机组热电负荷优化调度系统
CN109146182A (zh) 计及多种储能的分布式冷热电联供系统的经济调度方法
CN106593771B (zh) 一种热电联产系统中的电热解耦系统及其控制方法
CN107358345B (zh) 计及需求侧管理的分布式冷热电联供系统优化运行方法
Chang et al. Multi-objective optimization of a novel combined cooling, dehumidification and power system using improved M-PSO algorithm
CN112464477A (zh) 计及需求响应的多能耦合综合能源运行仿真方法
CN105117557B (zh) 一种楼宇型冷热电联供系统动力装置选型方法
CN109523065A (zh) 一种基于改进量子粒子群算法的微能源网优化调度方法
CN110361969B (zh) 一种冷热电综合能源系统优化运行方法
CN112531696B (zh) 一种考虑用户需求响应的冷热电联供系统优化调度方法
CN113435095A (zh) 一种用于综合能源系统的优化调度的方法及系统
CN112364556A (zh) 基于多能互补的智慧能源优化配置方法及终端设备
CN111799772A (zh) 一种考虑机组深度调峰的电热系统优化调度方法
CN111382902A (zh) 基于运行效益增量的区域综合能源系统储能优化配置方法
CN113642802A (zh) 一种基于碳计量模型的综合能源站能源优化调度方法和系统
CN107730032B (zh) 一种泛能网电冷热三联供优化调度系统及方法
CN113190941A (zh) 一种适用于电热联合调度的集中供热管网建模方法
Chen et al. A novel and optimized operation strategy map for CCHP systems considering optimal thermal energy utilization
CN115864376A (zh) 一种多时间尺度电力保供分析方法及系统
CN113128868B (zh) 一种区域综合能源系统调度优化方法及装置
CN106096248A (zh) 不同环境下的最佳供暖方式的甄选方法
CN112288219B (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