CN110661254A - 一种区域综合能源系统冷热电互补效益量化方法 - Google Patents
一种区域综合能源系统冷热电互补效益量化方法 Download PDFInfo
- Publication number
- CN110661254A CN110661254A CN201910915611.1A CN201910915611A CN110661254A CN 110661254 A CN110661254 A CN 110661254A CN 201910915611 A CN201910915611 A CN 201910915611A CN 110661254 A CN110661254 A CN 110661254A
- Authority
- CN
- China
- Prior art keywords
- power
- heat
- cost
- gas
- mathematical model
- 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.)
- Pending
Links
- 230000008901 benefit Effects 0.000 title claims abstract description 52
- 230000000295 complement effect Effects 0.000 title claims abstract description 44
- 238000010438 heat treatment Methods 0.000 title claims abstract description 39
- 238000001816 cooling Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000013178 mathematical model Methods 0.000 claims abstract description 48
- 230000005611 electricity Effects 0.000 claims abstract description 46
- 230000008878 coupling Effects 0.000 claims abstract description 9
- 238000010168 coupling process Methods 0.000 claims abstract description 9
- 238000005859 coupling reaction Methods 0.000 claims abstract description 9
- 238000011002 quantification Methods 0.000 claims abstract description 7
- 239000007789 gas Substances 0.000 claims description 79
- 238000012423 maintenance Methods 0.000 claims description 32
- 239000002918 waste heat Substances 0.000 claims description 26
- 238000010521 absorption reaction Methods 0.000 claims description 23
- 238000003860 storage Methods 0.000 claims description 23
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 18
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000010248 power generation Methods 0.000 claims description 10
- 239000000446 fuel Substances 0.000 claims description 9
- 239000003345 natural gas Substances 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 5
- 238000007599 discharging Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 3
- 101150067055 minC gene Proteins 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000002330 electrospray ionisation mass spectrometry Methods 0.000 claims description 2
- 239000002737 fuel gas Substances 0.000 claims description 2
- 238000005286 illumination Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Power Engineering (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Air Conditioning Control Device (AREA)
Abstract
本发明涉及一种区域综合能源系统冷热电互补效益量化方法,包括:S1:收集所研究区域基础数据;S2:建立CCHP型综合能源系统供能结构的数学模型;S3:确定约束条件及目标函数的确定;S4:对CCHP型综合能源系统供能结构的数学模型进行求解;S5:冷热电互补效益量化指标计算。本发明所述系统运行能完全满足区域内电、冷、热负荷的需求,基于各设备的处理情况计算冷热电互补效益量化指标,通过与冷热电分供系统经济性对比,可以量化出能源间互补耦合的价值;通过对比分析各个区域中综合能源系统冷热电互补的效益,可以有针对性地解决地区中哪个区域适合建设综合能源系统,对综合能源系统规划问题具有指导意义。
Description
技术领域
本发明涉及城市区域综合能源系统规划与优化领域,更具体地,涉及一种区域综合能源系统冷热电互补效益量化方法。
背景技术
综合能源系统作为一种包含可再生能源的分布式电源综合集成技术,是分布式发电的有效管理单元和重要组织形式,对于推进配电系统的环保性和经济性具有重要的意义。CCHP系统凭借其能效高、环境效益好等优势已成为综合能源系统的主要研究方向之一。
目前国内外对CCHP型综合能源系统优化调度模型和控制方法均有一定的研究,但都尚未明确量化综合能源互补的效益。研究的CCHP型综合能源系统模型结构包括典型CCHP系统母线式结构、集成太阳能和联供系统的太阳能冷热电三联供系统结构以及风光气储互补发电的CCHP模型结构;在进行优化调度时主要考虑了燃气费用和购电费用;控制方式包含电跟随和热跟随两种策略。一方面,现有研究考虑综合能源系统中的设备不全面;再有,优化过程中考虑的目标因素不全面,从而优化结果误差较大;最后,都未明确给出综合能源系统中能源耦合的价值。
发明内容
本发明为克服上述现有技术所述的对CCHP型综合能源系统优化调度模型和控制方法尚未明确量化综合能源互补的效益的缺陷,提供一种区域综合能源系统冷热电互补效益量化方法。
为了达到上述目的,本发明建立了包括燃气轮机、余热锅炉、燃气锅炉、电制冷机、吸收式制冷机、光伏和蓄电池等主要设备的CCHP型综合能源系统模型,采用以冷热电负荷平衡和总体经济性成本最优来确定设备出力的运行策略,针对模型特性使用基于0-1混合整数线性规划软件对模型进行求解,最后量化分析综合能源系统冷热电互补效益。通过分析各CCHP型综合能源系统内各设备出力、冷热电负荷平衡情况、经济性运行成本来验证所提出的模型和量化方法的有效性,同时对各区域是否适合建设综合能源系统进行深入探究。
所述方法包括以下步骤:
S1:收集所研究区域的冷、热、电负荷功率数据,区域的分时电价和天然气价格,以及系统内各种设备的运行参数;
S2:系统模型构建:建立CCHP型综合能源系统供能结构的数学模型;
S3:约束条件及目标函数的确定:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束;以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用;
S4:模型求解:根据S2和S3所建立的数学模型,对CCHP型综合能源系统供能结构的数学模型进行求解;
S5:冷热电互补效益量化指标计算:从经济性对综合能源系统冷热电互补耦合效益进行量化分析,包含日购电费用、日燃气费用、设备日运行维护费用和设备初始日投资费用四个指标将综合能源系统与冷热电分供系统进行对比分析,得出冷热电互补效益。
本发明提供一种区域综合能源系统冷热电互补效益量化方法。其包括收集所研究区域夏季、冬季和过渡期典型日24小时的冷、热、电负荷数据;综合能源系统模型构建;确定约束条件及目标函数;模型求解;冷热电互补效益量化指标(包括日购电费用、日燃气费用、设备日运行维护费用和设备初始日投资费用)计算等步骤。针对不同的能源供需应用场景,提出适应价值最优调控策略,深度挖掘多能源综合应用价值,发挥多能源的综合运营优势,是量化综合能源系统冷热电互补效益的目的所在。以中国北方某生态城中的居民区、商业区、办公区和工业区为算例,本发明提供的基于运行模拟的区域综合能源系统冷热电互补效益量化方法能够通过结合地区综合能源系统内各设备出力、冷热电负荷平衡情况,对比量化冷热电分供系统与冷热电联供系统的效益,突出综合能源系统冷热电互补的经济性优势,可为各区域是否适合建设综合能源系统提供参考依据。
优选地,S2中CCHP型综合能源系统包括:电网、光伏机组、蓄电池、电制冷机、冷负荷、燃气轮机、燃气锅炉、余热锅炉、转换装置、吸收式制冷机、热负荷、电负荷;
电网、光伏机组、燃气轮机、蓄电池、电负荷、电制冷机构成电力系统;光伏机组、燃气轮机为电力系统提供电能;电网和蓄电池既能为电力系统提供电能也能从电力系统吸收电能;电负荷和电制冷机从电力系统吸收电能;
燃气轮机为余热锅炉提供热能,余热锅炉为热转换装置、吸收式制冷机提供热能;燃气锅炉和热转换装置为热负荷提供热能;
吸收式制冷机、电制冷机为冷负荷提供冷能。
优选地,S2中CCHP型综合能源系统供能结构的数学模型包括燃气轮机数学模型、燃气锅炉数学模型、余热锅炉数学模型、热转换装置数学模型、电制冷机数学模型、吸收式制冷机数学模型、光伏机组数学模型、蓄电池数学模型;
(1)燃气轮机数学模型为:
ηGT=aβ3-bβ2+cβ+d
其中,ηGT为燃气轮机发电效率;PGT(t)为燃气轮机发电功率;QGT(t)为燃气轮机排出烟气中的余热量;a、b、c、d为正的常数,由燃气轮机发电效率的拟合曲线决定;β为机组电负荷率,即平均电负荷与最大电负荷的比值;
(2)燃气锅炉数学模型为:
QGB(t)=FGB(t)LNGηGB
其中,QGB(t)为燃气锅炉输出热量;FGB(t)为燃气锅炉的消耗燃气量;ηGB为燃气锅炉效率;LNG为燃气热值;
(3)余热锅炉数学模型为:
QWH(t)=QWH,in(t)ηWH
其中,QWH(t)为余热锅炉输出热量;QWH,in(t)为输入热量;ηWH为余热锅炉效率;
(4)换热装置数学模型为:
QHX(t)=QWH,heat(t)ηHX
其中,QHX(t)为换热装置输出热量;ηHX为换热装置效率;QWH,heat(t)为余热锅炉输出热量中用于供给系统热负荷的热量;
(5)电制冷机数学模型为:
QEC(t)=PEC(t)λEC
其中,QEC(t)为电制冷机输出制冷量;PEC(t)为电制冷机输入电功率;λEC为能效比;
(6)吸收式制冷机数学模型为:
QAC(t)=QWH,cool(t)ηAC
其中,QAC(t)为吸收式制冷机输出制冷量;ηAC为吸收式制冷机效率;QWH,cool(t)为余热锅炉输出热量中用于供给系统冷负荷的热量;
(7)光伏机组数学模型为:
(8)蓄电池数学模型为:
优选地,S3中约束条件及目标函数的确定方法包括下列步骤:
S3.1:确定目标函数:以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用:
minC=Ffuel+FRm+Fgrid
式中:C为CCHP型综合能源系统的总运行成本;Ffuel为日燃气费用;FRm为设备日运行维护费用;Fgrid为日购电费用;
S3.2:确定约束条件:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束,其中:
电功率平衡约束如公式为:
PGT(t)+PPV(t)+Pgrid(t)+PESD(t)-PESC(t)-PEC(t)=LE(t)
其中,PPV(t)为光伏发电输出的电功率,Pgrid(t)为综合能源系统与配电网交互的电功率,PESD(t)为蓄电池的放电功率,PESC(t)为蓄电池的充电功率,LE(t)为电负荷功率;
冷功率平衡约束如公式为:
λECPEC(t)+QAC(t)=LCA(t)
其中,LCA(t)为冷负荷功率;
热功率平衡约束如公式为:
QHX(t)+QGB(t)=LTW(t)
式中:LTW(t)分别为热负荷功率。
优选地,S5中冷热电互补效益量化指标计算方法为:
S5.1:计算日购电费用,计算公式为:
S5.2:计算日燃气费用,计算公式为:
S5.3:计算设备日运行维护费用,计算公式为:
kACQAC(t)+kECPEC(t)+kPVPPV(t)+kBTPBT(t)
式中:kGT,kGB和kWH分别为燃气轮机、燃气锅炉和余热锅炉单位时段单位功率运行维护费用;kHX,kAC和kEC分别为换热装置、吸收式制冷机和电制冷机单位时段单位功率运行维护费用;kPV和kBT分别为光伏发电和蓄电池单位时段单位功率运行维护费用,QWH(t)为余热锅炉输出的热功率,QHX(t)为换热装置输出的热功率,QAC(t)为吸收式制冷机输出的冷功率,PBT(t)为蓄电池输出(或输入)的电功率。
S5.4:计算设备初始日投资费用:
式中:Cinv,m为设备m的单位容量投资成本;λm为各设备安装容量;μcrf(r,Y)为资金回收系数,表示一次性投资成本折算到每年的费用支出;Y为设备的使用寿命;r为基准折现率;
S5.5:计算效益量化总费用:
F=Fgrid+Ffuel+FRm+Finv
其中,Finv为系统内所有设备的总初始日投资费用。
与现有技术相比,本发明技术方案的有益效果是:本发明提供的基于运行模拟的区域综合能源系统冷热电互补效益量化方法,系统运行能完全满足区域内电、冷、热负荷的需求,基于各设备的处理情况计算冷热电互补效益量化指标,通过与冷热电分供系统经济性对比,可以量化出能源间互补耦合的价值;通过对比分析各个区域中综合能源系统冷热电互补的效益,可以有针对性地解决地区中哪个区域适合建设综合能源系统,对综合能源系统规划问题具有指导意义。
附图说明
图1为实施例1所述区域综合能源系统冷热电互补效益量化方法流程图。
图2为CCHP型综合能源系统供能结构示意图。
图3a、图3b、图3c分别为居民区综合能源系统电、冷、热负荷平衡及设备出力曲线图。
图4a、图4b分别为商业区综合能源系统电、冷负荷平衡及设备出力曲线图。
图5a、图5b分别为办公区综合能源系统电、冷负荷平衡及设备出力曲线图。
图6a、图6b、图6c分别为工业区综合能源系统电、冷、热负荷平衡及设备出力曲线图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
本实施例提供一种区域综合能源系统冷热电互补效益量化方法,
如图1所示,所述方法包括以下步骤:
S1:收集包括所研究区域的冷、热、电负荷功率数据,区域的分时电价和天然气价格,以及系统内各种设备的运行参数等基础数据;
S2:系统模型构建:在Matlab建模界面中建立CCHP型综合能源系统供能结构的数学模型;
S3:约束条件及目标函数的确定:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束;以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用;
S4:模型求解:根据S2和S3所建立的数学模型,编写了基于Matlab线性规划模型的0-1混合整数规划程序对其求解,由于约束条件中包含耦合变量,如蓄电池的充放电功率,因此在模型中引入0-1变量;
S5:冷热电互补效益量化指标计算:从经济性对综合能源系统冷热电互补耦合效益进行量化分析,包含日购电费用、日燃气费用、设备日运行维护费用和设备初始日投资费用四个指标将综合能源系统与冷热电分供系统进行对比分析,得出冷热电互补效益。
在S1中,所述的收集整理数据的方法为:
收集所研究区域的夏季典型日冷、热、电负荷24小时功率数据,所在区域的分时电价和天然气价格,以及系统内各种设备的运行参数;
以我国北方某一生态城为例,生态城中包含居民区、商业区、办公区和工业区。收集到的居民区、商业区、办公区和工业区中夏季典型日冷、热、电负荷24小时功率数据如表1所示;各区域的分时电价和天然气价格数据如表2所示;综合能源系统内各种设备的运行参数数据如表3所示;综合能源系统内各种设备的安装投资成本及使用寿命数据如表4所示。
表1我国北方某一生态城内各区域的夏季典型日冷、热、电负荷24小时功率数据表
表2我国北方某一生态城内各区域的夏季典型日分时电价和天然气价格数据表
表3综合能源系统内各种设备的运行参数数据表
参数 | 数值及单位 |
燃气轮机额定效率 | 0.3 |
余热锅炉额定效率 | 0.8 |
燃气锅炉额定效率 | 0.9 |
换热装置额定效率 | 0.9 |
吸收式制冷机额定效率 | 1.2 |
电制冷机能效比 | 4.0 |
蓄电池充电效率 | 0.95 |
蓄电池放电效率 | 0.95 |
蓄电池自放电率 | 0.04 |
蓄电池容量 | 2000kWh |
蓄电池最大充电功率 | 400kW |
蓄电池最大放电功率 | 800kW |
光伏使用维护成本 | 0.08元/kWh |
蓄电池使用维护成本 | 0.02元/kWh |
燃气轮机使用维护成本 | 0.03元/kWh |
余热锅炉使用维护成本 | 0.02元/kWh |
燃气锅炉使用维护成本 | 0.02元/kWh |
换热装置使用维护成本 | 0.025元/kWh |
电制冷机使用维护成本 | 0.01元/kWh |
吸收式制冷机使用维护成本 | 0.025元/kWh |
表4综合能源系统内各种设备的安装投资成本及使用寿命数据表
在S2中,根据图2所示的CCHP型综合能源系统供能结构示意图,在Matlab建模界面中建立其数学模型,其数学模型表达式如表5所示;
表5综合能源系统中主要设备及其数学模型
在S3中,所述的约束条件及目标函数的确定方法包括下列步骤:
S3.1:确定目标函数:以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用:
minC=Ffuel+FRm+Fgrid (1)
式中:C为CCHP型综合能源系统的总运行成本;Ffuel为日燃气费用;FRm为设备日运行维护费用;Fgrid为日购电费用。
S3.2:确定约束条件:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束,其中:
电功率平衡约束如公式(2)所示:
PGT(t)+PPV(t)+Pgrid(t)+PESD(t)-PESC(t)-PEC(t)=LE(t) (2)
其中,PPV(t)为光伏发电输出的电功率,Pgrid(t)为综合能源系统与配电网交互的电功率,PESD(t)为蓄电池放出的电功率,PESC(t)为蓄电池吸收的电功率,LE(t)为电负荷功率;
冷功率平衡约束如公式(3)所示:
λECPEC(t)+QAC(t)=LCA(t) (3)
热功率平衡约束如公式(4)所示:
QHX(t)+QGB(t)=LTW(t) (4)
式中:LE(t),LCA(t),LTW(t)分别为电、冷、热负荷功率。
在S4中,所述的模型求解的方法为:根据S2和S3所建立的数学模型,编写了基于Matlab线性规划模型的0-1混合整数规划程序对其求解,由于约束条件中包含耦合变量,如蓄电池的充放电功率,因此在模型中引入0-1变量;运行程序,耗时0.9918秒,得到的居民区中各设备输出结果如图3a、图3b、图3c所示,得到的商业区中各设备输出结果如图4a、图4b所示,得到的办公区中各设备输出结果如图5a、图5b所示,得到的工业区中各设备输出结果如图6a、图6b、图6c所示。
在S5中,所述的冷热电互补效益量化指标计算方法为:
S5.1:计算日购电费用,如公式(5)所示:
S5.2:计算日燃气费用,如公式(6)所示:
S5.3:计算设备日运行维护费用,如公式(7)所示:
kACQAC(t)+kECPEC(t)+kPVPPV(t)+kBTPBT(t)
式中:kGT,kGB和kWH分别为燃气轮机、燃气锅炉和余热锅炉单位时段单位功率运行维护费用;kHX,kAC和kEC分别为换热装置、吸收式制冷机和电制冷机单位时段单位功率运行维护费用;kPV和kBT分别为光伏发电和蓄电池单位时段单位功率运行维护费用。QWH(t)为余热锅炉输出的热功率,QHX(t)为换热装置输出的热功率,QAC(t)为吸收式制冷机输出的冷功率,PBT(t)为蓄电池输出(或输入)的电功率。
S5.4:计算设备初始日投资费用:
式中:Cinv,m为设备m的单位容量投资成本;λm为各设备安装容量;μcrf(r,Y)为资金回收系数,表示一次性投资成本折算到每年的费用支出;Y为设备的使用寿命;r为基准折现率,本实施例取6.7%。
S5.5:计算效益量化总费用:
F=Fgrid+Ffuel+FRm+Finv (10)
其中,Finv为系统内所有设备的总初始日投资费用。
根据S4的输出结果以及S5的计算公式,得出我国北方某一生态城中居民区、商业区、办公区和工业区中冷热电互补效益量化结果,结果如表6所示;为了便于比较分析,将传统的冷热电分供系统(即电负荷由电网提供,冷负荷由电制冷机提供,热负荷由燃气锅炉提供)作为案例1,将冷热电互补耦合系统作为案例2。
表6不同区域不同案例费用对比
用包含日购电费用、日燃气费用、设备日运行维护费用和设备初始日投资费用的总费用来量化综合能源系统冷热电互补的效益。由表6可知,案例2相对于案例1节省了购电费用,增加了初始设备投资费用、燃气费用和设备运行维护费用。居民区案例2比案例1整体节约了382.1元,节省率约为1.68%;商业区案例2比案例1整体节约了6597元,节省率约为9.07%;办公区案例2比案例1整体节约了2371元,节省率约为8.09%;工业区案例2比案例1整体节约了3152元,节省率约为2.92%。由于案例2下燃气轮机出力增大,燃气费用略有上升,但同时其购电费用显著降低,且由于商业区和办公区均属于电价较高的地区,所以其在冷热电互补的效益上更加明显,更适合建设综合能源系统。
附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (5)
1.一种区域综合能源系统冷热电互补效益量化方法,其特征在于,所述方法包括以下步骤:
S1:收集所研究区域的冷、热、电负荷功率数据,区域的分时电价和天然气价格,以及系统内各种设备的运行参数;
S2:系统模型构建:建立CCHP型综合能源系统供能结构的数学模型;
S3:约束条件及目标函数的确定:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束;以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用;
S4:模型求解:根据S2和S3所建立的数学模型,对CCHP型综合能源系统供能结构的数学模型进行求解;
S5:冷热电互补效益量化指标计算:从经济性对综合能源系统冷热电互补耦合效益进行量化分析,包含日购电费用、日燃气费用、设备日运行维护费用和设备初始日投资费用四个指标将综合能源系统与冷热电分供系统进行对比分析,得出冷热电互补效益。
2.根据权利要求1所述的区域综合能源系统冷热电互补效益量化方法,其特征在于,S2中CCHP型综合能源系统包括:电网、光伏机组、蓄电池、电制冷机、冷负荷、燃气轮机、燃气锅炉、余热锅炉、转换装置、吸收式制冷机、热负荷、电负荷;
电网、光伏机组、燃气轮机、蓄电池、电负荷、电制冷机构成电力系统;光伏机组、燃气轮机为电力系统提供电能;电网和蓄电池既能为电力系统提供电能也能从电力系统吸收电能;电负荷和电制冷机从电力系统吸收电能;
燃气轮机为余热锅炉提供热能,余热锅炉为热转换装置、吸收式制冷机提供热能;燃气锅炉和热转换装置为热负荷提供热能;
吸收式制冷机、电制冷机为冷负荷提供冷能。
3.根据权利要求2所述的区域综合能源系统冷热电互补效益量化方法,其特征在于,S2中CCHP型综合能源系统供能结构的数学模型包括燃气轮机数学模型、燃气锅炉数学模型、余热锅炉数学模型、热转换装置数学模型、电制冷机数学模型、吸收式制冷机数学模型、光伏机组数学模型、蓄电池数学模型;
(1)燃气轮机数学模型为:
ηGT=aβ3-bβ2+cβ+d
其中,ηGT为燃气轮机发电效率;PGT(t)为燃气轮机发出的电功率;QGT(t)为燃气轮机排出烟气中的余热量;a、b、c、d为正的常数,由燃气轮机发电效率的拟合曲线决定;β为机组电负荷率,即平均电负荷与最大电负荷的比值;
(2)燃气锅炉数学模型为:
QGB(t)=FGB(t)LNGηGB
其中,QGB(t)为燃气锅炉输出热量;FGB(t)为燃气锅炉的消耗燃气量;ηGB为燃气锅炉效率;LNG为燃气热值;
(3)余热锅炉数学模型为:
QWH(t)=QWH,in(t)ηWH
其中,QWH(t)为余热锅炉输出热量;QWH,in(t)为输入热量;ηWH为余热锅炉效率;
(4)换热装置数学模型为:
QHX(t)=QWH,heat(t)ηHX
其中,QHX(t)为换热装置输出热量;ηHX为换热装置效率;QWH,heat(t)为余热锅炉输出热量中用于供给系统热负荷的热量;
(5)电制冷机数学模型为:
QEC(t)=PEC(t)λEC
其中,QEC(t)为电制冷机输出制冷量;PEC(t)为电制冷机输入电功率;λEC为能效比;
(6)吸收式制冷机数学模型为:
QAC(t)=QWH,cool(t)ηAC
其中,QAC(t)为吸收式制冷机输出制冷量;ηAC为吸收式制冷机效率;QWH,cool(t)为余热锅炉输出热量中用于供给系统冷负荷的热量;
(7)光伏机组数学模型为:
其中,为光伏机组预测出力;ηPV为太阳能电池板效率;S为电池板面积;为单位面积光照强度;
(8)蓄电池数学模型为:
4.根据权利要求3所述的区域综合能源系统冷热电互补效益量化方法,其特征在于,S3中
约束条件及目标函数的确定方法包括下列步骤:
S3.1:确定目标函数:以调度周期内总的运行成本最小为目标函数,包括日燃气费用、日购电费用和设备日运行维护费用:
minC=Ffuel+FRm+Fgrid
式中:C为CCHP型综合能源系统的总运行成本;Ffuel为日燃气费用;FRm为设备日运行维护费用;Fgrid为日购电费用;
S3.2:确定约束条件:CCHP型综合能源系统运行约束条件包括电、冷、热平衡约束,设备出力上下限约束,其中:
电功率平衡约束如公式为:
PGT(t)+PPV(t)+Pgrid(t)+PESD(t)-PESC(t)-PEC(t)=LE(t)
其中,PPV(t)为光伏发电输出的电功率,Pgrid(t)为综合能源系统与配电网交互的电功率,PESD(t)为蓄电池放出的电功率,PESC(t)为蓄电池吸收的电功率,LE(t)为电负荷功率;
冷功率平衡约束如公式为:
λECPEC(t)+QAC(t)=LCA(t)
其中,LCA(t)为冷负荷功率;
热功率平衡约束如公式为:
QHX(t)+QGB(t)=LTW(t)
式中:LTW(t)分别为热负荷功率。
5.根据权利要求4所述的区域综合能源系统冷热电互补效益量化方法,其特征在于,S5中冷热电互补效益量化指标计算方法为:
S5.1:计算日购电费用,计算公式为:
S5.2:计算日燃气费用,计算公式为:
S5.3:计算设备日运行维护费用,计算公式为:
式中:kGT,kGB和kWH分别为燃气轮机、燃气锅炉和余热锅炉单位时段单位功率运行维护费用;kHX,kAC和kEC分别为换热装置、吸收式制冷机和电制冷机单位时段单位功率运行维护费用;kPV和kBT分别为光伏发电和蓄电池单位时段单位功率运行维护费用,QWH(t)为余热锅炉输出的热功率,QHX(t)为换热装置输出的热功率,QAC(t)为吸收式制冷机输出的冷功率,PBT(t)为蓄电池输出(或输入)的电功率;
S5.4:计算设备初始日投资费用:
式中:Cinv,m为设备m的单位容量投资成本;λm为各设备安装容量;μcrf(r,Y)为资金回收系数,表示一次性投资成本折算到每年的费用支出;Y为设备的使用寿命;r为基准折现率;
S5.5:计算效益量化总费用:
F=Fgrid+Ffuel+FRm+Finv
其中,Finv为系统内所有设备的总初始日投资费用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910915611.1A CN110661254A (zh) | 2019-09-26 | 2019-09-26 | 一种区域综合能源系统冷热电互补效益量化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910915611.1A CN110661254A (zh) | 2019-09-26 | 2019-09-26 | 一种区域综合能源系统冷热电互补效益量化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110661254A true CN110661254A (zh) | 2020-01-07 |
Family
ID=69039208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910915611.1A Pending CN110661254A (zh) | 2019-09-26 | 2019-09-26 | 一种区域综合能源系统冷热电互补效益量化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110661254A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111382902A (zh) * | 2020-03-02 | 2020-07-07 | 天津大学 | 基于运行效益增量的区域综合能源系统储能优化配置方法 |
CN111524028A (zh) * | 2020-04-03 | 2020-08-11 | 中煤科工(天津)清洁能源研究院有限公司 | 一种多能互补耦合系统配置方法 |
CN111860966A (zh) * | 2020-06-24 | 2020-10-30 | 东南大学 | 基于模糊相关机会规划的含储能综合能源系统调度方法 |
CN111953018A (zh) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | 一种分布式多能互补供能系统及功率分配方法 |
CN111967659A (zh) * | 2020-08-03 | 2020-11-20 | 华北电力大学 | 一种基于光伏消纳的区域综合能源系统配置优化方法 |
CN112329263A (zh) * | 2020-11-24 | 2021-02-05 | 国网河北省电力有限公司经济技术研究院 | 一种综合能源系统优化方法及终端设备 |
CN112507507A (zh) * | 2020-10-12 | 2021-03-16 | 上海电力大学 | 一种基于经济性与可靠性的综合能源设备优化配置方法 |
CN113128868A (zh) * | 2021-04-20 | 2021-07-16 | 云南电网有限责任公司电力科学研究院 | 一种区域综合能源系统调度优化方法及装置 |
CN113469430A (zh) * | 2021-06-29 | 2021-10-01 | 华电电力科学研究院有限公司 | 综合用能园区的多能互补容量配置方法 |
CN114971000A (zh) * | 2022-05-18 | 2022-08-30 | 东南大学 | 一种计及供需不平衡率的综合能源系统配置优化方法 |
WO2024109327A1 (zh) * | 2022-11-21 | 2024-05-30 | 云南电网有限责任公司德宏供电局 | 一种基于多能互补的综合能源运行控制方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108717594A (zh) * | 2018-04-16 | 2018-10-30 | 东南大学 | 一种冷热电联供型多微网系统经济优化调度方法 |
CN109193772A (zh) * | 2018-09-29 | 2019-01-11 | 东北大学 | 一种基于风光微网的储能优化配置系统及方法 |
CN109492815A (zh) * | 2018-11-15 | 2019-03-19 | 郑州大学 | 一种市场机制下面向电网的储能电站选址定容优化方法 |
CN110163450A (zh) * | 2019-05-31 | 2019-08-23 | 国网山东省电力公司经济技术研究院 | 一种基于运行受限的配电网规划双层优化模型构建方法 |
-
2019
- 2019-09-26 CN CN201910915611.1A patent/CN110661254A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108717594A (zh) * | 2018-04-16 | 2018-10-30 | 东南大学 | 一种冷热电联供型多微网系统经济优化调度方法 |
CN109193772A (zh) * | 2018-09-29 | 2019-01-11 | 东北大学 | 一种基于风光微网的储能优化配置系统及方法 |
CN109492815A (zh) * | 2018-11-15 | 2019-03-19 | 郑州大学 | 一种市场机制下面向电网的储能电站选址定容优化方法 |
CN110163450A (zh) * | 2019-05-31 | 2019-08-23 | 国网山东省电力公司经济技术研究院 | 一种基于运行受限的配电网规划双层优化模型构建方法 |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111382902A (zh) * | 2020-03-02 | 2020-07-07 | 天津大学 | 基于运行效益增量的区域综合能源系统储能优化配置方法 |
CN111524028A (zh) * | 2020-04-03 | 2020-08-11 | 中煤科工(天津)清洁能源研究院有限公司 | 一种多能互补耦合系统配置方法 |
CN111860966A (zh) * | 2020-06-24 | 2020-10-30 | 东南大学 | 基于模糊相关机会规划的含储能综合能源系统调度方法 |
CN111860966B (zh) * | 2020-06-24 | 2024-05-10 | 东南大学 | 基于模糊相关机会规划的含储能综合能源系统调度方法 |
CN111967659B (zh) * | 2020-08-03 | 2023-10-31 | 华北电力大学 | 一种基于光伏消纳的区域综合能源系统配置优化方法 |
CN111967659A (zh) * | 2020-08-03 | 2020-11-20 | 华北电力大学 | 一种基于光伏消纳的区域综合能源系统配置优化方法 |
CN111953018A (zh) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | 一种分布式多能互补供能系统及功率分配方法 |
CN112507507A (zh) * | 2020-10-12 | 2021-03-16 | 上海电力大学 | 一种基于经济性与可靠性的综合能源设备优化配置方法 |
CN112329263A (zh) * | 2020-11-24 | 2021-02-05 | 国网河北省电力有限公司经济技术研究院 | 一种综合能源系统优化方法及终端设备 |
CN112329263B (zh) * | 2020-11-24 | 2023-02-03 | 国网河北省电力有限公司经济技术研究院 | 一种综合能源系统优化方法及终端设备 |
CN113128868A (zh) * | 2021-04-20 | 2021-07-16 | 云南电网有限责任公司电力科学研究院 | 一种区域综合能源系统调度优化方法及装置 |
CN113128868B (zh) * | 2021-04-20 | 2023-04-07 | 云南电网有限责任公司电力科学研究院 | 一种区域综合能源系统调度优化方法及装置 |
CN113469430A (zh) * | 2021-06-29 | 2021-10-01 | 华电电力科学研究院有限公司 | 综合用能园区的多能互补容量配置方法 |
CN113469430B (zh) * | 2021-06-29 | 2024-02-09 | 华电电力科学研究院有限公司 | 综合用能园区的多能互补容量配置方法 |
CN114971000A (zh) * | 2022-05-18 | 2022-08-30 | 东南大学 | 一种计及供需不平衡率的综合能源系统配置优化方法 |
CN114971000B (zh) * | 2022-05-18 | 2024-04-30 | 东南大学 | 一种计及供需不平衡率的综合能源系统配置优化方法 |
WO2024109327A1 (zh) * | 2022-11-21 | 2024-05-30 | 云南电网有限责任公司德宏供电局 | 一种基于多能互补的综合能源运行控制方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110661254A (zh) | 一种区域综合能源系统冷热电互补效益量化方法 | |
CN108537409B (zh) | 一种考虑多能源耦合特性的工业园区配电网协同规划方法 | |
Pan et al. | Multi-objective and two-stage optimization study of integrated energy systems considering P2G and integrated demand responses | |
CN109659927B (zh) | 一种考虑储能参与度的综合能源微网储能容量配置方法 | |
CN110458353B (zh) | 电热联合微网能量梯级优化方法及系统 | |
CN110288152B (zh) | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 | |
Guo et al. | Multi-objective optimization design and multi-attribute decision-making method of a distributed energy system based on nearly zero-energy community load forecasting | |
Ren et al. | A novel planning method for design and dispatch of hybrid energy systems | |
CN114091913A (zh) | 考虑热网及p2g多园区综合能源系统低碳经济调度方法 | |
Fan et al. | Two-layer collaborative optimization for a renewable energy system combining electricity storage, hydrogen storage, and heat storage | |
CN111668878A (zh) | 一种可再生微能源网的优化配置方法和系统 | |
Liu | Research on optimal placement of low‐carbon equipment capacity in integrated energy system considering carbon emission and carbon trading | |
Liu et al. | Capacity allocation for regional integrated energy system considering typical day economic operation | |
CN111382902A (zh) | 基于运行效益增量的区域综合能源系统储能优化配置方法 | |
CN115099007B (zh) | 基于综合成本-能耗曲线的综合能源系统优化运行方法 | |
CN115864376A (zh) | 一种多时间尺度电力保供分析方法及系统 | |
Ma et al. | Collaborative optimization scheduling of integrated energy system considering user dissatisfaction | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
Wakui et al. | Impact analysis of sampling time interval and battery installation on optimal operational planning of residential cogeneration systems without electric power export | |
CN114580746A (zh) | 一种基于低碳经济效益量化的综合能源站复合储能配置优化方法 | |
Guo et al. | Operation optimization of integrated energy system from the perspective of sustainable development | |
Zhao et al. | Multi-objective bi-level planning of the integrated energy system considering uncertain user loads and carbon emission during the equipment manufacturing process | |
CN113128799A (zh) | 能源管控方法、装置、电子设备及计算机存储介质 | |
CN111126675A (zh) | 多能互补微网系统优化方法 | |
Qiao et al. | Energy hub‐based configuration optimization method of integrated energy system |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200107 |
|
RJ01 | Rejection of invention patent application after publication |