CN107862419A - 一种多能源网络的低碳优化调度方法及系统 - Google Patents
一种多能源网络的低碳优化调度方法及系统 Download PDFInfo
- Publication number
- CN107862419A CN107862419A CN201711292668.8A CN201711292668A CN107862419A CN 107862419 A CN107862419 A CN 107862419A CN 201711292668 A CN201711292668 A CN 201711292668A CN 107862419 A CN107862419 A CN 107862419A
- Authority
- CN
- China
- Prior art keywords
- network
- power
- carbon
- mrow
- object function
- 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
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 92
- 238000005457 optimization Methods 0.000 title claims abstract description 9
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 58
- 238000010438 heat treatment Methods 0.000 claims abstract description 44
- 238000009826 distribution Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000008878 coupling Effects 0.000 claims abstract description 14
- 238000010168 coupling process Methods 0.000 claims abstract description 14
- 238000005859 coupling reaction Methods 0.000 claims abstract description 14
- 230000005611 electricity Effects 0.000 claims description 28
- 230000005540 biological transmission Effects 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 8
- 238000013213 extrapolation Methods 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 239000008400 supply water Substances 0.000 claims description 4
- 239000004744 fabric Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 3
- 230000009467 reduction Effects 0.000 description 5
- 230000005619 thermoelectricity Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 239000002028 Biomass Substances 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000003763 carbonization Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
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/06313—Resource planning in a project environment
-
- 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
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种多能源网络的低碳优化调度方法及系统。该方法包括:采用支路潮流模型描述多能源网络中的电网的潮流分布;采用水力热力耦合模型描述多能源网络中的区域供热网络的潮流分布;确定多能源网络的低碳调度模型;所述低碳调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布;确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合;根据所述非劣解集合,利用纳什协商问题确定最佳权衡解;根据所述最佳权衡解对所述多能源网络进行调度。本发明的方法及系统供多能源网络运行调度人员使用,优化调度过程考虑成本与碳排双目标,客观快捷。
Description
技术领域
本发明涉及多能源网络的运行调控领域,特别是涉及一种多能源网络的低碳优化调度方法及系统。
背景技术
在当前环境气候挑战不断加剧、节能减排形势日益严峻的背景下,实现能源需求供应的低碳化迫在眉睫。从终端需求来看,热能和电能是最主要的能源需求形式,占终端能源需求的60%以上。一方面,当前终端电能需求主要由区域配电网(Power DistributionNetwork,PDN)供给,PDN可有效融合主动配电网、微电网等技术,高效接纳分布式可再生能源机组,实现电力供应的清洁化。另一方面,结合热泵、生物质、光热等技术,区域集中式供热系统(District Heating Network,DHN)日渐成为实现热能供应清洁化的主要手段之一。特别地,在能源互联网背景下,随着电气化供热技术的普及与能源网络互联进程的深入,PDN、DHN等能源网络将紧密耦合,形成更具运行灵活性的多能源网络。实现该类多能源网路的节能减排优化调度是发挥其灵活性的必要手段。
尽管目前已有涉及区域电网节能减排调度或区域供热管网节能减排优化调度的研究,但对于多能源网路的节能减排优化调度研究还很鲜见。节能减排优化调度是发挥多能源网络低碳环保优势的重要措施,目前常用的加权平均法等手段往往具有较强主观性。特别地,多能源网络运行灵活性更强,网络特性更为复杂,用于区域电网节能减排调度或区域供热管网节能减排优化调度方式不适用于多能源网络节能减排的调度,因此多能源网络的节能减排优化调度更具挑战。
发明内容
本发明的目的是提供一种多能源网络的低碳优化调度方法及系统,以解决传统调度方法主观性强且不适用于多能源网络的问题。
为实现上述目的,本发明提供了如下方案:
一种多能源网络的低碳优化调度方法,所述方法包括:
采用支路潮流模型描述所述多能源网络中的电网的潮流分布;
采用水力热力耦合模型描述所述多能源网络中的区域供热网络的潮流分布;
确定所述多能源网络的低碳调度模型;所述低碳调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布;
确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合;
根据所述非劣解集合,利用纳什协商问题确定最佳权衡解;
根据所述最佳权衡解对所述多能源网络进行优化调度。
可选的,所述采用支路潮流模型描述多能源网络中的电网的潮流分布,具体包括:
利用模型描述多能源网络中的电网的潮流分布;
其中,t表示调度时刻,Pij,t为线路传输有功功率,Qij,t为线路传输无功功率,为节点j的有功负荷需求;为节点j的无功负荷需求;为电源有功功率出力;为电源的无功功率出力;π(j)为节点j的子节点集合;xij为线路电抗,rij为线路电阻,zij为线路阻抗,Uj,t为节点电压幅值的平方;Iij,t为线路电流的平方,上标l表示下界,上标u表示上界。
可选的,所述采用水力热力耦合模型描述多能源网络中的区域供热网络的潮流分布,具体包括:
利用模型描述多能源网络中的区域供热网络的潮流分布
其中,c为载热流体的比热容,分别为供回水管网出口及入口温度,b为管道索引,t为调度时刻索引; 分别为供回水管道节点温度,为环境温度;分别为供回水管网载热流体质量流速;λb为管道温度损耗系数,Lb为管道长度;F(i)表示以节点i为首节点的管道集合,T(i)表示以节点i为末节点的管道集合;表示热源供热量,表示热负荷热需求。
可选的,所述确定所述多能源网络的低碳调度模型,具体包括:
运行成本:
其中,fC为所述多能源网络的运行成本,ai与bi为电网机组i的运行成本系数,为电网机组i的电源有功功率出力;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足 表示热网机组j的热源供热量;
碳排量:
其中,fE为所述多能源网络的碳排量,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数;
将运行成本最低确定为第一目标函数;
将碳排量最低确定为第二目标函数;
将所述第一目标函数、所述第二目标函数及所述约束条件确定为所述多能源网络的低碳调度模型。
可选的,所述确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合,具体包括:
利用正交边界交叉法获取均匀分布的帕累托前沿解,获得所述非劣解集合,所述非劣解集合的元素为所述帕累托前沿解。
可选的,所述根据所述非劣解集合,利用纳什协商问题确定最佳权衡解,具体包括:
确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
确定所述纳什协商问题的等效二阶锥松弛形式为
根据所述纳什协商问题的等效二阶锥松弛形式获得纳什协商解;
将所述纳什协商解确定为所述最佳权衡解。
一种多能源网络的低碳优化调度系统,所述系统包括:
电网潮流确定模块,用于采用支路潮流模型描述所述多能源网络中的电网的潮流分布;
区域供热网络潮流确定模块,用于采用水力热力耦合模型描述所述多能源网络中的区域供热网络的潮流分布;
调度模型确定模块,用于确定所述多能源网络的低碳调度模型,所述低碳调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布;
非劣解集合确定模块,用于确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合;
最佳权衡解确定模块,用于根据所述非劣解集合,利用纳什协商问题确定最佳权衡解;
优化调度模块,用于根据所述最佳权衡解对所述多能源网络进行优化调度。
可选的,所述调度模型确定模块,具体包括:
运行成本确定单元,用于利用模型确定所述多能源网络的运行成本;其中,fC为所述多能源网络的运行成本,ai与bi为电网机组i的运行成本系数,为电网机组i的电源有功功率出力;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足 表示热网机组j的热源供热量;
碳排量确定单元,用于利用模型确定所述多能源网络的碳排量;其中,fE为所述多能源网络的碳排量,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数;
第一目标函数确定单元,用于将运行成本最低确定为第一目标函数;
第二目标函数确定单元,用于将碳排量最低确定为第二目标函数;
调度模型确定单元,用于将所述第一目标函数、所述第二目标函数及所述约束条件确定为所述多能源网络的调度模型。
可选的,所述非劣解集合确定模块利用正交边界交叉法获取均匀分布的帕累托前沿解,获得所述非劣解集合,所述非劣解集合的元素为所述帕累托前沿解。
可选的,所述最佳权衡解确定模块,具体包括:
纳什协商问题确定单元,用于确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
等效二阶锥形式确定单元,用于确定所述纳什协商问题的等效二阶锥松弛形式为
纳什协商解获得单元,用于根据所述纳什协商问题的等效二阶锥松弛形式获得所述纳什协商解;
最佳权衡解确定单元,用于将所述纳什协商解确定为所述最佳权衡解。。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
计及多能源网络多能流损耗,采用支路潮流模型描述电网潮流分布,采用水力热力耦合模型描述区域供热网络潮流分布,具有运行成本最小与碳排最小两个目标。采用正交边界交叉法获得均匀分布的帕累托前沿,为进一步获取最佳权衡解,基于合作博弈思想,将双目标转化为纳什协商问题,并采用二阶锥松弛规划方法求解纳什协商解,从而实现基于合作博弈的成本与碳排双目标间的最佳权衡,供多能源网络运行调度人员使用,优化调度过程客观快捷。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明多能源网络的低碳优化调度方法的流程示意图;
图2为本发明多能源网络的低碳优化调度系统的结构示意图;
图3为本发明具体实施方式的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明多能源网络的低碳优化调度方法的流程示意图。如图1所示,所述方法包括:
步骤100:确定多能源网络中电网的潮流分布。采用支路潮流模型确定多能源网络中的电网的潮流分布,具体包括:
其中,t表示调度时刻,Pij,t为线路传输有功功率,Qij,t为线路传输无功功率,为节点j的有功负荷需求;为节点j的无功负荷需求;为电源有功功率出力;为电源的无功功率出力;π(j)为节点j的子节点集合;xij为线路电抗,rij为线路电阻,zij为线路阻抗,Uj,t为节点电压幅值的平方;Iij,t为线路电流的平方,上标l表示下界,上标u表示上界。
步骤200:确定多能源网络中的区域供热网络的潮流分布。采用水力热力耦合模型确定所述多能源网络中的区域供热网络的潮流分布,具体包括:
其中,c为载热流体的比热容,分别为供回水管网出口及入口温度,b为管道索引,t为调度时刻索引;分别为供回水管道节点温度,为环境温度;分别为供回水管网载热流体质量流速;λb为管道温度损耗系数,Lb为管道长度;F(i)表示以节点i为首节点的管道集合,T(i)表示以节点i为末节点的管道集合;表示热源供热量,表示热负荷热需求。
步骤300:确定多能源网络的调度模型。多能源网络节能减排优化调度具有最小化运行成本与最小化碳排两个目标,因此,所述调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布。确定多能源网络的调度模型具体包括:
利用模型确定所述多能源网络的运行成本;其中,ai与bi为电网机组i的运行成本系数;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足
利用模型确定所述多能源网络的碳排量;其中,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数,表示热网机组j的热源供热量;
将运行成本最低确定为第一目标函数,即第一目标函数为minfC;
将碳排量最低确定为第二目标函数,即第二目标函数为minfE;
将所述第一目标函数、所述第二目标函数及所述的区域电网与供热网络的潮流分布为所述多能源网络的调度模型,约束条件为步骤100和步骤200确定的潮流分布。
步骤400:确定多能源网络调度模型中第一目标函数和第二目标函数的非劣解集合。采用正交边界交叉法获取均匀分布的帕累托前沿解。非劣解集合的元素为所述帕累托前沿解。
步骤500:利用纳什协商问题确定最佳权衡解。基于合作博弈思想,将节能减排优化调度模型转为纳什协商问题,并采用二阶锥松弛规划方法进行求解。具体如下:
确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述热源包括热电联产机组,所述电源包括燃气轮机,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
确定所述纳什协商问题的等效二阶锥松弛形式为
根据所述纳什协商问题的等效二阶锥松弛形式获得纳什协商解;
将所述纳什协商解确定为所述最佳权衡解。
步骤600:对多能源网络的机组出力进行优化调度。步骤500中获得最佳权衡解,最佳权衡解包括运行成本和碳排量的具体值,即最佳的成本与碳排权衡解。每个数值对应调度模型中的决策变量的参数。决策变量包括:各热源出力、电源出力、配电网从电网购电量、电网运行状态变量(节点电压、线路传输有功无功等)、热网运行状态变量(节点供回水温度等)。从而根据最佳权衡解对调度模型中的决策变量的参数进行优化调度。
图2为本发明多能源网络的低碳优化调度系统的结构示意图。如图2所示,所述系统包括:
电网潮流确定模块201,用于采用支路潮流模型确定多能源网络中的电网的潮流分布。
区域供热网络潮流确定模块202,用于采用水力热力耦合模型确定所述多能源网络中的区域供热网络的潮流分布。
调度模型确定模块203,用于确定所述多能源网络的调度模型,所述调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布。
调度模型确定模块203,具体包括:
运行成本确定单元,用于利用模型确定所述多能源网络的运行成本;其中,ai与bi为电网机组i的运行成本系数;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足
碳排量确定单元,用于利用模型确定所述多能源网络的碳排量;其中,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数;
第一目标函数确定单元,用于将运行成本最低确定为第一目标函数;
第二目标函数确定单元,用于将碳排量最低确定为第二目标函数;
调度模型确定单元,用于将所述第一目标函数、所述第二目标函数及所述电网潮流分布与所述热网潮流分布确定为所述多能源网络的调度模型。
非劣解集合确定模块204,用于确定所述调度模型中所述第一目标函数和所述第二目标函数的非劣解集合。
最佳权衡解确定模块205,用于根据所述非劣解集合,利用纳什协商问题确定最佳权衡解。
最佳权衡解确定模块205,具体包括:
纳什协商问题确定单元,用于确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
等效二阶锥形式确定单元,用于确定所述纳什协商问题的等效二阶锥松弛形式为
纳什协商解获得单元,用于根据所述纳什协商问题的等效二阶锥松弛形式获得所述纳什协商解;
最佳权衡解确定单元,用于将所述纳什协商解确定为所述最佳权衡解。
优化调度模块206,用于根据所述最佳权衡解对所述多能源网络的机组出力进行优化调度。
图3为本发明具体实施方式的流程示意图。如图3所示,具体实施时,针对热电多能流网络,采用支路潮流法描述配电网潮流,采用水力热力耦合模型描述区域供热网络潮流,计及多能源网络多能流损耗,构建基于运行成本最小与碳排最小的多能源网络节能减排优化调度方法。采用正交边界交叉法获得均匀分布的帕雷托前沿,为进一步获取成本与碳排最佳权衡解,基于合作博弈思想,将双目标转化为纳什协商问题,并采用二阶锥松弛规划方法获取纳什协商最优解,从而实现成本与碳排双目标间的最佳权衡,供多能源网络运行调度人员使用。包括以下内容:
(1)计及多能源网络损耗,采用支路潮流描述电网潮流,采用水力热力耦合模型描述区域供热管网潮流。同时也是后续调度模型的约束条件。
采用支路潮流模型描述热电多能源网络中配电网潮流分布,其优势在于可监测区域电网的节点电压、线路无功等,同时该模型可被线性化或采用二阶锥松弛进行高效求解;
支路潮流模型如下所示:
其中,t表示调度时刻,Pij,t为线路传输有功功率,Qij,t为线路传输无功功率,为节点j的有功负荷需求;为节点j的无功负荷需求;为电源有功功率出力;为电源的无功功率出力;π(j)为节点j的子节点集合;xij为线路电抗,rij为线路电阻,zij为线路阻抗,Uj,t为节点电压幅值的平方;Iij,t为线路电流的平方,上标l表示下界,上标u表示上界。
采用水力热力耦合模型描述热电多能源网络中区域供热网络的潮流分布,如下所示:
其中,c为载热流体的比热容,分别为供回水管网出口及入口温度,b为管道索引,t为调度时刻索引;分别为供回水管道节点温度,为环境温度;分别为供回水管网载热流体质量流速;λb为管道温度损耗系数,Lb为管道长度;F(i)表示以节点i为首节点的管道集合,T(i)表示以节点i为末节点的管道集合;表示热源供热量,表示热负荷热需求。
(2)多能源网络节能减排优化调度具有最小化运行成本与最小化碳排两个目标。决策变量主要包括各热源出力、电源出力、配电网从电网购电量、电网运行状态变量(节点电压、线路传输有功无功等)、热网运行状态变量(节点供回水温度等)。
(3)多能源网络节能减排优化调度模型求解方法:
运行成本与碳排量两个目标一般无法同时达到最小值,实际热电多能源网络中运行调度人员需要两种调度策略,一是一系列非劣解集合(帕累托前沿);二是在帕累托前沿上获取一个最佳权衡解。
针对第一种情况,采用正交边界交叉法获取均匀分布的帕累托前沿前沿解。此时,调度人员根据工作需求可以从非劣解集合中选取一个帕累托前沿前沿解进行调度。
对于第二种情况,在得到非劣解集合的基础上,综合考虑各种因素,确定一个最佳权衡解,座位调度的策略。基于合作博弈思想,将节能减排优化调度模型转为Nash协商问题并采用二阶锥松弛规划方法进行求解。
本具体实施方式的优点是:计及多能源网络多能流损耗,采用支路潮流模型描述配电网潮流分布,采用水力热力耦合模型描述区域集中供热网络潮流分布,具有运行成本最小与碳排最小两个目标。采用正交边界交叉法获得均匀分布的帕累托前沿,为进一步获取最佳权衡解,基于合作博弈思想,将双目标转化为纳什协商问题,并采用二阶锥松弛规划方法求解纳什协商解,从而实现基于合作博弈的成本与碳排双目标间的最佳权衡,供多能源网络运行调度人员使用。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种多能源网络的低碳优化调度方法,其特征在于,所述方法包括:
采用支路潮流模型描述所述多能源网络中的电网的潮流分布;
采用水力热力耦合模型描述所述多能源网络中的区域供热网络的潮流分布;
确定所述多能源网络的低碳调度模型;所述低碳调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布;
确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合;
根据所述非劣解集合,利用纳什协商问题确定最佳权衡解;
根据所述最佳权衡解对所述多能源网络进行优化调度。
2.根据权利要求1所述的方法,其特征在于,所述采用支路潮流模型描述多能源网络中的电网的潮流分布,具体包括:
利用模型描述多能源网络中的电网的潮流分布;
其中,t表示调度时刻,Pij,t为线路传输有功功率,Qij,t为线路传输无功功率,为节点j的有功负荷需求;为节点j的无功负荷需求;为电源有功功率出力;为电源的无功功率出力;π(j)为节点j的子节点集合;xij为线路电抗,rij为线路电阻,zij为线路阻抗,Uj,t为节点电压幅值的平方;Iij,t为线路电流的平方,上标l表示下界,上标u表示上界。
3.根据权利要求1所述的方法,其特征在于,所述采用水力热力耦合模型描述多能源网络中的区域供热网络的潮流分布,具体包括:
利用模型描述多能源网络中的区域供热网络的潮流分布
其中,c为载热流体的比热容,分别为供回水管网出口及入口温度,b为管道索引,t为调度时刻索引; 分别为供回水管道节点温度,为环境温度;分别为供回水管网载热流体质量流速;λb为管道温度损耗系数,Lb为管道长度;F(i)表示以节点i为首节点的管道集合,T(i)表示以节点i为末节点的管道集合;表示热源供热量,表示热负荷热需求。
4.根据权利要求1所述的方法,其特征在于,所述确定所述多能源网络的低碳调度模型,具体包括:
运行成本:
其中,fC为所述多能源网络的运行成本,ai与bi为电网机组i的运行成本系数,为电网机组i的电源有功功率出力;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足 表示热网机组j的热源供热量;
碳排量:
其中,fE为所述多能源网络的碳排量,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数;
将运行成本最低确定为第一目标函数;
将碳排量最低确定为第二目标函数;
将所述第一目标函数、所述第二目标函数及所述约束条件确定为所述多能源网络的低碳调度模型。
5.根据权利要求1所述的方法,其特征在于,所述确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合,具体包括:
利用正交边界交叉法获取均匀分布的帕累托前沿解,获得所述非劣解集合,所述非劣解集合的元素为所述帕累托前沿解。
6.根据权利要求1所述的方法,其特征在于,所述根据所述非劣解集合,利用纳什协商问题确定最佳权衡解,具体包括:
确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
确定所述纳什协商问题的等效二阶锥松弛形式为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mi>max</mi>
<mrow>
<mi>x</mi>
<mo>&Element;</mo>
<msub>
<mi>X</mi>
<mi>P</mi>
</msub>
</mrow>
</munder>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>&GreaterEqual;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>f</mi>
<mi>C</mi>
<mi>d</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>f</mi>
<mi>E</mi>
<mi>d</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>E</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
根据所述纳什协商问题的等效二阶锥松弛形式获得纳什协商解;
将所述纳什协商解确定为所述最佳权衡解。
7.一种多能源网络的低碳优化调度系统,其特征在于,所述系统包括:
电网潮流确定模块,用于采用支路潮流模型描述所述多能源网络中的电网的潮流分布;
区域供热网络潮流确定模块,用于采用水力热力耦合模型描述所述多能源网络中的区域供热网络的潮流分布;
调度模型确定模块,用于确定所述多能源网络的低碳调度模型,所述低碳调度模型中包括运行成本最低的第一目标函数和碳排量最低的第二目标函数,约束条件为所述电网的潮流分布和所述区域供热网络的潮流分布;
非劣解集合确定模块,用于确定所述低碳调度模型中所述第一目标函数和所述第二目标函数的非劣解集合;
最佳权衡解确定模块,用于根据所述非劣解集合,利用纳什协商问题确定最佳权衡解;
优化调度模块,用于根据所述最佳权衡解对所述多能源网络进行优化调度。
8.根据权利要求7所述的系统,其特征在于,所述调度模型确定模块,具体包括:
运行成本确定单元,用于利用模型确定所述多能源网络的运行成本;其中,fC为所述多能源网络的运行成本,ai与bi为电网机组i的运行成本系数,为电网机组i的电源有功功率出力;θt为上级电网售电电价;ξj与ωj为热网机组j的运行成本系数;为从上级电网的购电量,满足 表示热网机组j的热源供热量;
碳排量确定单元,用于利用模型确定所述多能源网络的碳排量;其中,fE为所述多能源网络的碳排量,αi,βi分别为电网机组i的碳排系数;ζj与为热网机组j的碳排系数;
第一目标函数确定单元,用于将运行成本最低确定为第一目标函数;
第二目标函数确定单元,用于将碳排量最低确定为第二目标函数;
调度模型确定单元,用于将所述第一目标函数、所述第二目标函数及所述约束条件确定为所述多能源网络的调度模型。
9.根据权利要求7所述的系统,其特征在于,所述非劣解集合确定模块利用正交边界交叉法获取均匀分布的帕累托前沿解,获得所述非劣解集合,所述非劣解集合的元素为所述帕累托前沿解。
10.根据权利要求7所述的系统,其特征在于,所述最佳权衡解确定模块,具体包括:
纳什协商问题确定单元,用于确定纳什协商问题为第三目标函数其中为运行成本最大值,为碳排量最大值,Xp为所述非劣解集合中帕累托前沿解对应的决策变量集合,所述决策变量集合包括各热源出力、电源出力、从电网购电量、电网运行状态变量、热网运行状态变量,所述电网运行状态变量包括节点电压、线路传输有功功率和无功功率,所述热网运行状态包括节点供回水温度;
等效二阶锥形式确定单元,用于确定所述纳什协商问题的等效二阶锥松弛形式为
纳什协商解获得单元,用于根据所述纳什协商问题的等效二阶锥松弛形式获得所述纳什协商解;
最佳权衡解确定单元,用于将所述纳什协商解确定为所述最佳权衡解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711292668.8A CN107862419A (zh) | 2017-12-08 | 2017-12-08 | 一种多能源网络的低碳优化调度方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711292668.8A CN107862419A (zh) | 2017-12-08 | 2017-12-08 | 一种多能源网络的低碳优化调度方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107862419A true CN107862419A (zh) | 2018-03-30 |
Family
ID=61705184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711292668.8A Pending CN107862419A (zh) | 2017-12-08 | 2017-12-08 | 一种多能源网络的低碳优化调度方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107862419A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102104A (zh) * | 2018-06-29 | 2018-12-28 | 清华大学 | 一种基于碳排放流的多能源系统协同低碳规划方法 |
CN110210104A (zh) * | 2019-05-28 | 2019-09-06 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
CN115375183A (zh) * | 2022-09-19 | 2022-11-22 | 浙江大学 | 考虑能流-碳流耦合约束的电力系统配网优化调度方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236694A (zh) * | 2013-04-28 | 2013-08-07 | 广东电网公司电力科学研究院 | 一种适用于主动配电网最优潮流求解的方法 |
US20140222665A1 (en) * | 2011-11-28 | 2014-08-07 | Expanergy, Llc | Energy search engine methods and systems |
CN105046353A (zh) * | 2015-07-06 | 2015-11-11 | 国家电网公司 | 一种电力系统低碳化水平的评价方法 |
CN106159974A (zh) * | 2016-08-02 | 2016-11-23 | 清华大学 | 一种输配协调的分布式无功电压优化方法 |
CN106253350A (zh) * | 2016-08-11 | 2016-12-21 | 清华大学 | 基于供热管网储热效益的热‑电联合机组组合方法 |
CN106339772A (zh) * | 2016-08-11 | 2017-01-18 | 清华大学 | 基于供热管网储热效益的热‑电联合优化调度方法 |
CN106779361A (zh) * | 2016-12-02 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 并网混合可再生能源系统的多目标优化设计方法 |
CN106815661A (zh) * | 2017-02-22 | 2017-06-09 | 清华大学 | 一种热电联合系统的分解协调调度方法 |
CN107368961A (zh) * | 2017-07-12 | 2017-11-21 | 东南大学 | 一种适用于新能源接入背景下的区域电网碳排放管理方法 |
-
2017
- 2017-12-08 CN CN201711292668.8A patent/CN107862419A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140222665A1 (en) * | 2011-11-28 | 2014-08-07 | Expanergy, Llc | Energy search engine methods and systems |
CN103236694A (zh) * | 2013-04-28 | 2013-08-07 | 广东电网公司电力科学研究院 | 一种适用于主动配电网最优潮流求解的方法 |
CN105046353A (zh) * | 2015-07-06 | 2015-11-11 | 国家电网公司 | 一种电力系统低碳化水平的评价方法 |
CN106159974A (zh) * | 2016-08-02 | 2016-11-23 | 清华大学 | 一种输配协调的分布式无功电压优化方法 |
CN106253350A (zh) * | 2016-08-11 | 2016-12-21 | 清华大学 | 基于供热管网储热效益的热‑电联合机组组合方法 |
CN106339772A (zh) * | 2016-08-11 | 2017-01-18 | 清华大学 | 基于供热管网储热效益的热‑电联合优化调度方法 |
CN106779361A (zh) * | 2016-12-02 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 并网混合可再生能源系统的多目标优化设计方法 |
CN106815661A (zh) * | 2017-02-22 | 2017-06-09 | 清华大学 | 一种热电联合系统的分解协调调度方法 |
CN107368961A (zh) * | 2017-07-12 | 2017-11-21 | 东南大学 | 一种适用于新能源接入背景下的区域电网碳排放管理方法 |
Non-Patent Citations (4)
Title |
---|
WEI WEI 等: "Nash Bargain and Complementarity Approach Based Environmental/Economic Dispatch", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
WEI WEI 等: "Network Equilibrium of Coupled Transportation and Power Distribution Systems", 《IEEE TRANSACTIONS ON SMART GRID》 * |
WEI WEI 等: "Optimal Power Flow of Radial Networks and Its Variations: A Sequential Convex Optimization Approach", 《IEEE TRANSACTIONS ON SMART GRID》 * |
谢鑫鹏 等: "低碳供应链生产及交易决策机制", 《控制与决策》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102104A (zh) * | 2018-06-29 | 2018-12-28 | 清华大学 | 一种基于碳排放流的多能源系统协同低碳规划方法 |
CN109102104B (zh) * | 2018-06-29 | 2022-02-18 | 清华大学 | 一种基于碳排放流的多能源系统协同低碳规划方法 |
CN110210104A (zh) * | 2019-05-28 | 2019-09-06 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
CN115375183A (zh) * | 2022-09-19 | 2022-11-22 | 浙江大学 | 考虑能流-碳流耦合约束的电力系统配网优化调度方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106339772B (zh) | 基于供热管网储热效益的热-电联合优化调度方法 | |
CN104734168B (zh) | 一种基于电热联合调度的微电网运行优化系统及方法 | |
CN105914758B (zh) | 考虑分布式电源无功补偿成本的电网无功优化方法 | |
CN102520675B (zh) | 燃气联合循环与太阳能发电联合制热系统及其调度方法 | |
CN106253350B (zh) | 基于供热管网储热效益的热-电联合机组组合方法 | |
CN102506519B (zh) | 热电联产机组与风力发电联合供热系统及调度方法 | |
CN110175311B (zh) | 一种基于多能耦合模型的优化潮流计算方法 | |
CN107862419A (zh) | 一种多能源网络的低碳优化调度方法及系统 | |
CN106786753A (zh) | 多用户的区域能源互联网的系统及其调节方法 | |
CN106055773B (zh) | 结合热网模型的多区域综合能源系统配置模型的建立方法 | |
CN106447529A (zh) | 考虑热水管网的分布式能源系统建模和运行优化方法 | |
CN108960503A (zh) | 基于内点法的综合能源系统多场景优化分析方法 | |
CN105958537B (zh) | 面向能源互联网的能源转换系统及其优化控制方法 | |
CN111582658A (zh) | 考虑多区域互联协同的能源站设备配置及管线规划方法 | |
CN104362650B (zh) | 一种考虑成本因素的电力系统无功优化方法 | |
CN109740824A (zh) | 一种考虑供热网和热负荷的多能源系统优化规划方法 | |
CN112182887A (zh) | 一种综合能源系统规划优化仿真方法 | |
CN107990398A (zh) | 一种基于热力网、电力网和物联网的电蓄热供热系统 | |
CN112035984A (zh) | 一种电-气-储区域综合能源系统协同规划方法 | |
CN109449934A (zh) | 一种基于高压配电网分区重构的输电网阻塞管控方法 | |
CN109256766A (zh) | 一种能源互联网中电热耦合系统 | |
CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
CN113837589A (zh) | 一种综合能源系统n-1调度方法及装置 | |
CN108491973A (zh) | 一种考虑需求侧响应的电-热联合系统优化调度方法 | |
CN104319822B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180330 |
|
RJ01 | Rejection of invention patent application after publication |