CN111342452B - 一种多区域电气综合能源系统能量与备用分散式调度方法 - Google Patents
一种多区域电气综合能源系统能量与备用分散式调度方法 Download PDFInfo
- Publication number
- CN111342452B CN111342452B CN202010182634.9A CN202010182634A CN111342452B CN 111342452 B CN111342452 B CN 111342452B CN 202010182634 A CN202010182634 A CN 202010182634A CN 111342452 B CN111342452 B CN 111342452B
- Authority
- CN
- China
- Prior art keywords
- standby
- gas
- constraint
- region
- energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000003993 interaction Effects 0.000 claims abstract description 6
- 239000007789 gas Substances 0.000 claims description 113
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 92
- 238000005457 optimization Methods 0.000 claims description 63
- 239000003345 natural gas Substances 0.000 claims description 46
- 230000008878 coupling Effects 0.000 claims description 36
- 238000010168 coupling process Methods 0.000 claims description 36
- 238000005859 coupling reaction Methods 0.000 claims description 36
- 230000009194 climbing Effects 0.000 claims description 15
- 238000010248 power generation Methods 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 10
- 230000009977 dual effect Effects 0.000 claims description 10
- 230000009191 jumping Effects 0.000 claims description 5
- 230000010485 coping Effects 0.000 claims description 4
- 230000002040 relaxant effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000000295 complement effect Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 241001334134 Rugopharynx epsilon Species 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
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/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
- 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
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种多区域电气综合能源系统能量与备用分散式调度方法,首先对多区域电气综合能源系统能量‑备用协同优化调度问题进行数学建模,然后基于联络线与联络管道解耦的思想,将多区域集中式优化问题分解为区域优化调度子问题,最后根据交替方向乘子法,设计了一种分布式的能量‑备用优化调度模型求解策略。所提方法仅需区域间交互有限的联络线与联络管道信息,不需要交换每个区域内部的能源系统信息,从而可以较好地保护不同区域的信息隐私,解决了传统的多区域电气综合能源系统能量‑备用决策难以保护不同区域的信息隐私的问题。
Description
技术领域
本发明属于能源系统运行优化技术领域,具体为一种多区域电气综合能源系统能量与备用分散式调度方法。
背景技术
电力系统与天然气系统相互耦合构成的电气综合能源系统由于其能源利用效率高、有助于波动性清洁能源消纳、减少碳排放量等优点受到了政府、企业和科研机构等单位广泛地关注,而制定合理的电气综合能源系统优化调度策略是充分发挥能源系统协同互补优势的基石。近年来,以光伏和风电为代表的强波动性、间歇性能源快速发展,综合能源系统中清洁能源的渗透比例不断增加,同时,负荷侧以电动汽车等为代表的不确定性负荷也处于快速增长的时期。电源和负荷的双重不确定性使得电气综合能源系统运行面对的风险不断增加。因此,如何在管控不确定性因素带来的风险的前提下,科学合理地制定优化调度策略是当今综合能源系统领域关注的重点问题之一。
备用作为一种有效地应对不确定性风险的主要方式,已在能源系统中得到了广泛地实践和应用。由于供能装置如发电机、气井的备用可调节能力大小与其实时运行状态有关,因此,备用优化调度通常与能量优化调度进行协同优化。由于区域互联的能源系统可以实现不同区域的能量-备用协同互补,并降低系统运行成本和运行风险,多区域互联电气综合能源系统的协调调度问题开始受到了关注。传统多区域互联电气综合能源系统协同优化调度均采用一种集中式地优化决策方式,即假设一个独立的系统运营商负责整个多区域能源系统优化调度策略的制定,系统运营商需掌握所有互联区域内部和相关连接通道的信息,然而,不同区域的能源系统往往隶属于不同的运营或管理主体,运营或管理主体间出于行业竞争或其它方面的考虑,不会将内部的能源系统信息共享给其它主体。因此,迫切需要一种去中心化、分散式地优化决策方法来解决多区域互联电气综合能源系统能量-备用协同优化问题。
发明内容
本发明提供了一种多区域电气综合能源系统能量与备用分散式调度方法,用于解决现有的不同区域的能源系统往往隶属于不同的运营或管理主体,运营或管理主体间出于行业竞争或其它方面的考虑,不会将内部的能源系统信息共享给其它主体,从而影响多区域能源之间协同互补的问题。
为实现上述目的,本发明所采用的技术方案是:
一种多区域电气综合能源系统能量与备用分散式调度方法,包括以下步骤:
步骤1:建立多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型;
步骤2:基于联络线与联络管道解耦的思想,将步骤1中的多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型进行分解,得到区域能量-备用协同优化调度子问题分解模型;进而将多区域集中式优化问题分解为区域优化调度子问题;
步骤3:根据步骤2得到的区域能量-备用协同优化调度子问题分解模型,基于交替方向乘子法,设计了一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略。
优选的,所述多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型的建立步骤如下:
步骤1.1,建立目标函数,以调度周期内所有区域运行成本之和最小化为目标,包括常规机组发电成本、气井供气成本、弃风成本、向上备用容量成本与向下备用容量成本:
优化模型以调度周期内所有区域运行成本之和最小化为目标,包括常规机组发电成本、气井供气成本、弃风成本、向上备用容量成本与向下备用容量成本。
式中:a为区域标识;φ为区域集合;i与NE,c为发电机组标识与常规发电机组集合;gw与MG为气井标识与集合;w与NW为风电场标识与集合;r为备用源标识,包括发电机与气井;与/>分别表示常规机组单位发电成本、气井单位供气成本与单位弃风惩罚成本;Pa,i,t、Fa,gw,t与/>分别表示在时段t常规发电机组的输出功率、气井的供气量与弃风功率;/>和/>分别表示单位向上和向下备用容量成本;/>与/>分别表示备用供应源r在t时段提供的上下备用容量,基于书写方便,在此节后续表征系统运行约束时,均忽略区域标识a。
步骤1.2,在目标函数中增加电力系统运行约束、天然气系统运行约束、备用约束与耦合约束进行优化。
优选的,所述电力系统运行约束包括:
常规机组与燃气轮机输出功率约束:
风电场输出功率约束:
机组爬坡约束:
功率平衡约束:
线路潮流约束:
式中:Pi,t与Pj,t分别表示常规机组i与燃气轮机j在时段t的输出功率;Pw,t与分别表示风电场w在时段t的输出功率与预测功率;ri u与ri d为常规机组i以时段t为间隔的上下爬坡速率;rj u与rj d为燃气轮机j以时段t为间隔的上下爬坡速率;/>为电负荷d在时段t的负荷需求;Ne为电负荷集合;Kl,i、Kl,j、Kl,w与Kl,d分别为常规机组、燃气轮机、风电场与负荷对第l条线路的功率转移分布系数,其可由直流潮流模型计算得到;efl,t为线路l在时段t流过的功率;(·)min/max表示最小/最大限值。
优选的,所述天然气系统运行约束包括:
气井输出气量约束;
气井爬坡约束:
节点气压约束:
管道气流约束:
节点气流平衡约束:
式中:Fgw,t表示气井在时段t的输出气量;与/>分别表示气井上下爬坡速率;πm,t为节点m的气压;gfmr为管道mr流过的气流;Cmr为管道综合参数,其由管道长度、直径、摩擦力等综合决定;/>为天然气负荷d在时段t的负荷需求值;/>为燃气轮机在时段t的耗气量;r∈m表示与节点m相连的节点r。
优选的,所述备用约束在电力系统侧时,备用用于应对电负荷需求与风电输出功率不确定性带来的风险。基于机会约束理论,设置一定的置信水平来表征备用约束成立的概率,从而避免备用决策结果过于保守,约束包括:
上、下备用约束:
电力系统上备用与下备用机会约束的物理意义可表征为:在所有可能的风电输出功率和电负荷需求波动情况下,系统所提供的备用容量满足向上和向下功率调整需求的概率不低于置信水平β1和β2,即上、下备用约束成立的概率不低于β1和β2。
常规机组和燃气轮机可提供的备用容量约束:
所述备用约束在天然气系统侧时,备用用于应对天然气负荷需求不确定性和燃气轮机的天然气备用量需求。基于机会约束理论,设置一定的置信水平来表征备用约束成立的概率,从而避免备用决策结果过于保守,约束包括:
上、下备用约束:
天然气系统上备用与下备用机会约束的物理意义可表征为:在所有可能的燃气轮机备用需求和气负荷需求波动情况下,系统所提供的备用容量满足向上和向下供气调整需求的概率不低于置信水平β3和β4,即上、下备用约束成立的概率不低于置信水平β3和β4。
气井可提供的备用容量约束:
式中:NE,c为常规发电机组集合;NE,g为燃气轮机机组集合;Ne为电负荷集合;NW为风电场集合;Mg为气井集合;与/>分别表示常规机组与燃气轮机在时段t提供的向上调整备用容量;/>与/>分别表示常规机组与燃气轮机在时段t提供的向下调整备用容量;与ΔPw,t分别表示负荷需求与风电场预测功率偏离预测值的波动量;β1与β2分别表示电力系统上下备用约束成立的置信水平;/>与/>分别表示气井在时段t提供的向上与向下调整备用容量;/>为天然气负荷偏离预测值的波动量;/>与/>分别表示燃气轮机在电力系统侧提供向上与向下备用容量所需的天然气量;β3与β4分别表示天然气系统上下备用约束成立的置信水平。
优选的,所述耦合约束包括:
燃气轮机发电耗气量约束:
燃气轮机提供备用所耗天然气量约束:
式中:φj燃气轮机单位发电所耗天然气量。
优选的,所述步骤2还包括以下子步骤:
步骤2.1,基于联络线与联络管道等效解耦的思想,将多区域耦合互联电气综合能源系统分解为独立的区域子系统,并在此基础上,构建区域间耦合一致性约束。
步骤2.2,将一致性约束进行松弛并将其作为惩罚项加入到每个区域的优化目标函数中。
步骤2.3,分解后相应的区域内功率平衡和备用约束需要进行重构以计及区域间的能量和备用交换,从而可将原有多区域集中优化模型分解为多区域优化子模型。
优选的,所述一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略的过程如下:
步骤3.1,初始化;
步骤3.2,设置迭代次数r=0,设置原始和对偶残差收敛阈值,设置区域间交互变量信息初始值,拉格朗日乘子初始值与设置惩罚参数初始值;
步骤3.3,交换区域间耦合变量信息并其平均值;
步骤3.4,各个区域求解区域能量-备用优化调度子问题,并更新区域间耦合变量信息;
步骤3.5,每个区域计算原始残差与对偶残差,
步骤3.6,若最大原始残差与对偶残差满足收敛阈值要求,则终止迭代,输出优化调度结果;反之,跳转到步骤3.7;
步骤3.7,更新惩罚参数与拉格朗日乘子,设置r=r+1,跳转到步骤3.3。
本发明的有益技术效果是:
与传统集中式多区域互联电气综合能源系统能量-备用优化方法相比,本专利所提出的多区域互联电气综合能源系统能量-备用分散式协同优化调度方法,可以有效地保护不同区域能源子系统的信息隐私,求解方法仅需互联区域间交互有限的联络线与联络管道耦合信息。
附图说明
图1显示为本发明实施例的区域分解示意图;
图2显示为本发明的求解策略过程示意图;
图3显示为本发明的仿真系统示意图;
图4显示为本发明的分散式求解策略迭代收敛过程;
图5显示为本发明的不同区域优化目标函数收敛过程示意图。
具体实施方式
下面结合本发明的附图1-5,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种多区域电气综合能源系统能量与备用分散式调度方法,包括以下步骤:
步骤1:建立多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型;
步骤2:基于联络线与联络管道解耦的思想,将步骤1中的多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型进行分解,得到区域能量-备用协同优化调度子问题分解模型;进而将多区域集中式优化问题分解为区域优化调度子问题;
步骤3:根据步骤2得到的区域能量-备用协同优化调度子问题分解模型,基于交替方向乘子法,得到一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略。
进一步的,所述多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型的建立步骤如下:
步骤1.1,建立目标函数:
式中:a为区域标识;φ为区域集合;i与NE,c为发电机组标识与常规发电机组集合;gw与MG为气井标识与集合;w与NW为风电场标识与集合;r为备用源标识,包括发电机与气井;与/>分别表示常规机组单位发电成本、气井单位供气成本与单位弃风惩罚成本;Pa,i,t、Fa,gw,t与/>分别表示在时段t常规发电机组的输出功率、气井的供气量与弃风功率;/>和/>分别表示单位向上和向下备用容量成本;/>与/>分别表示备用供应源r在t时段提供的上下备用容量。
步骤1.2,在目标函数中增加电力系统运行约束、天然气系统运行约束、备用约束与耦合约束。
进一步的,所述电力协同运行约束包括:
常规机组与燃气轮机输出功率约束:
风电场输出功率约束:
机组爬坡约束:
功率平衡约束:
线路潮流约束:
式中:Pi,t与Pj,t分别表示常规机组i与燃气轮机j在时段t的输出功率;Pw,t与分别表示风电场w在时段t的输出功率与预测功率;ri u与ri d为常规机组i以时段t为间隔的上下爬坡速率;rj u与rj d为燃气轮机j以时段t为间隔的上下爬坡速率;/>为电负荷d在时段t的负荷需求;Ne为电负荷集合;Kl,i、Kl,j、Kl,w与Kl,d分别为常规机组、燃气轮机、风电场与负荷对第l条线路的功率转移分布系数,其可由直流潮流模型计算得到;efl,t为线路l在时段t流过的功率;(·)min/max表示最小/最大限值。
进一步的,所述天然气系统运行约束包括:
井输出气量约束;
气井爬坡约束:
节点气压约束:
管道气流约束:
节点气流平衡约束:
式中:Fgw,t表示气井在时段t的输出气量;与/>分别表示气井上下爬坡速率;πm,t为节点m的气压;gfmr为管道mr流过的气流;Cmr为管道综合参数,其由管道长度、直径、摩擦力等综合决定;/>为天然气负荷d在时段t的负荷需求值;/>为燃气轮机在时段t的耗气量;r∈m表示与节点m相连的节点r。
进一步的,在电力系统侧,备用用于应对电负荷需求与风电输出功率不确定性带来的风险。本专利基于机会约束理论,设置一定的置信水平来表征备用约束成立的概率,从而避免备用决策结果过于保守,约束包括上、下备用约束(13)以及常规机组和燃气轮机可提供的备用容量约束(14)。在天然气系统侧,备用用于应对天然气负荷需求不确定性和燃气轮机的天然气备用量需求,包括上、下备用约束(15)与气井可提供的备用容量约束(16):
上、下备用约束:
电力系统上备用与下备用机会约束的物理意义可表征为:在所有可能的风电输出功率和电负荷需求波动情况下,系统所提供的备用容量满足向上和向下功率调整需求的概率不低于置信水平β1和β2,即上、下备用约束成立的概率不低于β1和β2。
常规机组和燃气轮机可提供的备用容量约束:
所述备用约束在天然气系统侧包括:
上、下备用约束:
天然气系统上备用与下备用机会约束的物理意义可表征为:在所有可能的燃气轮机备用需求和气负荷需求波动情况下,系统所提供的备用容量满足向上和向下供气调整需求的概率不低于置信水平β3和β4,即上、下备用约束成立的概率不低于置信水平β3和β4。
气井可提供的备用容量约束:
式中:NE,c为常规发电机组集合;NE,g为燃气轮机机组集合;Ne为电负荷集合;NW为风电场集合;Mg为气井集合与/>分别表示常规机组与燃气轮机在时段t提供的向上调整备用容量;/>与/>分别表示常规机组与燃气轮机在时段t提供的向下调整备用容量;与ΔPw,t分别表示负荷需求与风电场预测功率偏离预测值的波动量;β1与β2分别表示电力系统上下备用约束成立的置信水平;/>与/>分别表示气井在时段t提供的向上与向下调整备用容量;/>为天然气负荷偏离预测值的波动量;/>与/>分别表示燃气轮机在电力系统侧提供向上与向下备用容量所需的天然气量;β3与β4分别表示天然气系统上下备用约束成立的置信水平。
进一步的,电力系统与天然气系统以燃气轮机作为关键耦合元件,耦合约束包括燃气轮机发电耗气量约束(17)与提供上下调整备用容量所需的天然气量约束(18)。
燃气轮机发电耗气量约束:
燃气轮机提供备用所耗天然气量约束:
式中:φj燃气轮机单位发电所耗天然气量,所建立的能量-备用协同优化调度模型中的机会约束(13)(15)可通过采样平均估计的方法转化为确定性的表达式,非线性管道气流约束(10)(11)可通过增量线性化方法转化为线性约束。
进一步的,所述步骤2还包括以下子步骤:
步骤2.1,基于联络线与联络管道等效解耦的思想,对多区域耦合互联电气综合能源系统分解为独立的区域子系统,并在此基础上,构建区域间耦合一致性约束;
多区域电气综合能源系统物理层面通过联络线与联络管道相连接,各个区域由独立的系统运行商负责,不同系统运行商之间可以交互联络线与联络管道耦合相关变量信息,包括能量流、备用流、节点气压等信息。通过基于联络线与联络管道等效解耦的思想,多区域协同优化调度问可划分为子问题,分解示意图如图1所示。这样分解过后,每个区域可以单独决策自己区域内的能量-备用协同优化调度子问题,并将相关决策后的联络线与联络管道上的能量流动、备用流动、节点气压等变量信息发布给相邻区域。显然,当互联区域交换的联络线与联络管道耦合的变量信息相等(一致)时(如下耦合一致性式(19)(20)所示),则原有多区域协同优化调度问题达到收敛。
式中:efa,kl,t与efb,kl,t分别为由区域A与区域B决策的在时段t流过联络线kl的功率;与/>分别为由区域A与区域B决策的在时段t交换的电力系统向上调整备用容量;与/>分别为由区域A与区域B决策的在时段t交换的电力系统向下调整备用容量;gfa,mn,t与gfb,mn,t分别为由区域A与区域B决策的在时段t流过联络管道mn的天然气量;/>与/>分别为由区域A与区域B决策的在时段t交换的天然气系统向上调整备用容量;与/>分别为由区域A与区域B决策的在时段t交换的天然气系统向下调整备用容量;πa,m,t与πb,m,t分别为由区域A与区域B决策的联络管道节点m的气压。
步骤2.2,将区域间耦合一致性约束式(19)-(20)进行松弛,并将其加入到各个区域自身优化目标函数中。以区域A为例,其增广拉格朗日优化目标函数为:
其中:
式中:与/>为区域间交换变量的拉格朗日乘子;/>与/>为惩罚系数;/>与/>为交换变量信息的平均值;/>为与区域A相连接的联络线集合;/>表示与区域A相连接的联络管道集合;/>为与区域A相连接的联络管道节点集合。
步骤2.3,在上述分解后,相应的区域内功率平衡和备用等约束需要进行重构以计及区域间的能量和备用交换,具体约束如下式所示:
式中:efa,kl,t与gfa,mn,t分别为区域A在时段t流出的功率和天然气;与/>分别为区域A在时段t流出的电力系统向上和向下调整备用容量;/>与/>分别为区域A在时段t流出的天然气系统向上和向下调整备用容量。
进一步的,所述一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略的过程如下:
步骤3.1,初始化;
步骤3.2,设置迭代次数r=0,设置原始和对偶残差收敛阈值εpri、εdual,设置区域间交互变量信息初始值,包括与/>设置拉格朗日乘子初始值,包括/>与/>设置惩罚参数初始值,包括与/>
步骤3.3,交换区域间耦合变量信息并其平均值。每个区域向相连区域发布其最新的区域间耦合变量信息与/>并通过式(22)-(23)计算信息交互后的变量平均值;
步骤3.4,各个区域求解区域能量-备用优化调度子问题,并更新区域间耦合变量信息为与/>计算区域间耦合变量信息的平均值/>与/>
步骤3.5,收敛校核。每个区域通过式(28)-(29)计算原始残差和通过式(30)-(31)计算对偶残差。
步骤3.6,若最大原始与对偶残差满足式(32)-(33),则终止迭代,输出优化调度结果,输出优化调度结果;反之,跳转到步骤3.7;
/>
步骤3.7,更新惩罚参数与拉格朗日乘子。以惩罚参数与拉格朗日乘子/>为例,惩罚参数通过式(34)更新,拉格朗日乘子通过式(35)更新,其余惩罚参数与拉格朗日乘子均可通过同样的计算方式更新,设置r=r+1,跳转到步骤3.3。
实施例2:
在实施例1的基础之上,本专利采用3区域互联的电气综合能源系统作为仿真系统来验证所提方法的有效性,每个区域均由6节点电力系统与6节点天然气系统组成,电力系统与天然气系统间通过燃气轮机进行耦合连接,仿真系统示意图如图3所示,具体算例数据见https://github.com/scugw/Case-Study-Parameter/blob/master/Multi-area% 20IEGS.pdf。算例调度周期T为一天。设置惩罚参数初始值为0.1,τin与τde均为2,μ设置为10,分别设置残差收敛阈值为0.1、0.05,对比不同收敛阈值下的优化结果。上述模型中的机会约束采用采样平均估计的方法可转化为线性约束,天然气潮流约束可通过增量线性化方法转化为线性约束,从而可将原有混合整数非线性优化问题转化为混合整数线性规划问题。在MATLAB R2016a中调用商业求解器Gurobi 8.1对上述优化问题进行求解。
为了表明本专利所提分散式求解方法的有效性,本专利首先对比了集中式优化方法与分散式优化方法求解多区域电气综合能源系统能量与备用协同优化调度模型的结果,如表1所示。集中式优化方法即通过商业求解器Gurobi求解优化问题(1)-(18)(机会约束采用采样平均估计的方法可转化为线性约束,天然气潮流约束可通过增量线性化方法转化为线性约束)。
表1
/>
由表1可知,在不同的收敛阈值,分散式优化得到的总运行成本与集中式优化得到的总运行成本相近,差异均小于1%;并且收敛阈值设置的越小,分散式优化获得的成本与集中式优化结果更加相近。由此可见,基于交替方向乘子法的多区域分散式优化求解策略可以实现调度优化问题可靠求解的同时保护不同区域的信息隐私。分散式求解策略迭代收敛过程如图4所示,不同区域优化目标函数收敛过程如图5所示。
进一步,本专利对比分析多区域电气综合能源系统分别在能量与备用优化调度上的协同优势,分别设置下4个场景,如下所示,其中场景2-4均采用上述分散式求解策略进行求解,收敛阈值设为0.5。
场景1:各个区域间不交互能量与备用,仅单独优化各自区域的能量与备用
场景2:各个区域间仅交互备用容量,即考虑多区域电气综合能源系统在备用优化方面的协同互补
场景3:各个区域仅交能量,即考虑多区域电气综合能源系统在能量优化方面的协同互补
场景4:各个区域交互能量与备用,即考虑多区域电气综合能源系统能量与备用互补协同优化。
不同场景下的运行总成本与各区域间的能量、备用总交互情况如表2所示,由表2可知,本专利所提多区域电气综合能源能量与备用一体化协同优化调度可以充分利用不同区域间的协调互补优势,从而降低系统运行成本。
表2
在本发明的描述中,需要理解的是,术语“逆时针”、“顺时针”“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
Claims (4)
1.一种多区域电气综合能源系统能量与备用分散式调度方法,其特征在于,包括以下步骤:
步骤1:建立多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型;
步骤2:基于联络线与联络管道解耦的思想,将步骤1中的多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型进行分解,得到区域能量-备用协同优化调度子问题分解模型;进而将多区域集中式优化问题分解为区域优化调度子问题;
步骤3:根据步骤2得到的区域能量-备用协同优化调度子问题分解模型,基于交替方向乘子法,设计了一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略;
所述多区域耦合互联电气综合能源系统能量与备用一体化优化调度模型的建立步骤如下:
步骤1.1,建立目标函数,以调度周期内所有区域运行成本之和最小化为目标,包括常规机组发电成本、气井供气成本、弃风成本、向上备用容量成本与向下备用容量成本:
式中:a为区域标识;φ为区域集合;i与NE,c为常规机组标识与常规机组集合;gw与MG为气井标识与集合;w与NW为风电场标识与集合;r为备用源标识,包括发电机与气井; 与/>分别表示常规机组单位发电成本、气井单位供气成本与单位弃风惩罚成本;Pa,i,t、Fa,gw,t与/>分别表示在时段t常规机组的输出功率、气井的供气量与弃风功率;/>和/>分别表示单位向上和向下备用容量成本;/>与/>分别表示备用供应源r在t时段提供的上下备用容量;
步骤1.2,在模型中增加电力系统运行约束、天然气系统运行约束、备用约束与耦合约束进行优化;
所述电力系统运行约束包括:
常规机组与燃气轮机输出功率约束:
风电场输出功率约束:
常规机组爬坡约束:
功率平衡约束:
线路潮流约束:
式中:Pi,t与Pj,t分别表示常规机组i与燃气轮机j在时段t的输出功率;Pw,t与分别表示风电场w在时段t的输出功率与预测功率;ri u与ri d为常规机组i以时段t为间隔的上、下爬坡速率;/>与/>为燃气轮机j以时段t为间隔的上下爬坡速率;/>为电负荷d在时段t的负荷需求;Ne为电负荷集合;Kl,i、Kl,j、Kl,w与Kl,d分别为常规机组、燃气轮机、风电场与负荷对第l条线路的功率转移分布系数,其可由直流潮流模型计算得到;efl,t为线路l在时段t流过的功率;
所述天然气系统运行约束包括:
气井输出气量约束;
气井爬坡约束:
节点气压约束:
管道气流约束:
节点气流平衡约束:
式中:Fgw,t表示气井在时段t的输出气量;与/>分别表示气井上下爬坡速率;πm,t为时段t内节点m的气压;gfmr为管道mr流过的气流;Cmr为管道综合参数,其由管道长度、直径、摩擦力综合决定;/>为天然气负荷d在时段t的负荷需求值;/>为燃气轮机在时段t的耗气量;r∈m表示与节点m相连的节点r;
所述备用约束在电力系统侧时,备用用于应对电负荷需求与风电输出功率不确定性带来的风险,基于机会约束理论,设置一定的置信水平来表征备用约束成立的概率,从而避免备用决策结果过于保守,约束包括:
上、下备用约束:
电力系统上备用与下备用机会约束的物理意义可表征为:在所有可能的风电输出功率和电负荷需求波动情况下,系统所提供的备用容量满足向上和向下功率调整需求的概率不低于置信水平β1和β2,即上、下备用约束成立的概率不低于β1和β2;
常规机组和燃气轮机可提供的上、下备用容量约束:
所述备用约束在天然气系统侧时,备用用于应对天然气负荷需求不确定性和燃气轮机的天然气备用量需求,基于机会约束理论,设置一定的置信水平来表征备用约束成立的概率,从而避免备用决策结果过于保守,约束包括:
上、下备用约束:
天然气系统上备用与下备用机会约束的物理意义可表征为:在所有可能的燃气轮机备用需求和气负荷需求波动情况下,系统所提供的备用容量满足向上和向下供气调整需求的概率不低于置信水平β3和β4,即上、下备用约束成立的概率不低于置信水平β3和β4;
气井可提供的备用容量约束:
式中:NE,c为常规机组集合;NE,g为燃气轮机机组集合;Ne为电负荷集合;NW为风电场集合;Mg为气井集合;与/>分别表示常规机组与燃气轮机在时段t提供的向上调整备用容量;/>与/>分别表示常规机组与燃气轮机在时段t提供的向下调整备用容量;/>与ΔPw,t分别表示负荷需求与风电场预测功率偏离预测值的波动量;β1与β2分别表示电力系统上下备用约束成立的置信水平;/>与/>分别表示气井在时段t提供的向上与向下调整备用容量;/>为天然气负荷偏离预测值的波动量;/>与/>分别表示燃气轮机在电力系统侧提供向上与向下备用容量所需的天然气量;β3与β4分别表示天然气系统上下备用约束成立的置信水平。
2.根据权利要求1所述的一种多区域电气综合能源系统能量与备用分散式调度方法,其特征在于,所述耦合约束包括:
燃气轮机发电耗气量约束:
燃气轮机提供备用所耗天然气量约束:
式中:φj燃气轮机单位发电所耗天然气量。
3.根据权利要求1所述的一种多区域电气综合能源系统能量与备用分散式调度方法,其特征在于,所述步骤2还包括以下子步骤:
步骤2.1,基于联络线与联络管道等效解耦的思想,将多区域耦合互联电气能源系统分解为独立的区域子系统,并在此基础上,构建区域间耦合一致性约束;
步骤2.2,将区域间耦合一致性约束进行松弛并将其作为惩罚项加入到每个区域的优化目标函数中;
步骤2.3,分解后相应的区域内功率平衡和备用约束需要进行重构以计及区域间的能量和备用交换,从而可将原有多区域集中优化模型分解为多区域优化子模型。
4.根据权利要求1所述的一种多区域电气综合能源系统能量与备用分散式调度方法,其特征在于,所述一种多区域互联电气综合能源系统能量-备用优化调度分散式求解策略的过程如下:
步骤3.1,初始化;
步骤3.2,设置迭代次数r=0,设置原始和对偶残差收敛阈值,设置区域间交互变量信息初始值,拉格朗日乘子初始值与设置惩罚参数初始值;
步骤3.3,交换区域间耦合变量信息并其平均值;
步骤3.4,各个区域求解区域能量-备用优化调度子问题,并更新区域间耦合变量信息;
步骤3.5,每个区域计算原始残差与对偶残差,
步骤3.6,若最大原始残差与对偶残差满足收敛阈值要求,则终止迭代,输出优化调度结果;反之,跳转到步骤3.7;
步骤3.7,更新惩罚参数与拉格朗日乘子,设置r=r+1,跳转到步骤3.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010182634.9A CN111342452B (zh) | 2020-03-16 | 2020-03-16 | 一种多区域电气综合能源系统能量与备用分散式调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010182634.9A CN111342452B (zh) | 2020-03-16 | 2020-03-16 | 一种多区域电气综合能源系统能量与备用分散式调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111342452A CN111342452A (zh) | 2020-06-26 |
CN111342452B true CN111342452B (zh) | 2023-09-12 |
Family
ID=71186095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010182634.9A Active CN111342452B (zh) | 2020-03-16 | 2020-03-16 | 一种多区域电气综合能源系统能量与备用分散式调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111342452B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112258077A (zh) * | 2020-09-28 | 2021-01-22 | 国网浙江省电力有限公司经济技术研究院 | 一种多区互联综合能源系统的参数确定方法及相关装置 |
CN112380681A (zh) * | 2020-09-28 | 2021-02-19 | 国网浙江省电力有限公司经济技术研究院 | 基于分布式优化的综合能源系统参数确定方法及相关装置 |
CN112803463B (zh) * | 2021-03-08 | 2022-07-08 | 浙江浙能电力股份有限公司萧山发电厂 | 计及电网运行约束的电化学储能应用场景评估方法 |
CN113111578B (zh) * | 2021-04-01 | 2023-04-07 | 上海晨翘智能科技有限公司 | 电力负荷预测方法、装置、计算机设备及存储介质 |
CN113762632B (zh) * | 2021-09-10 | 2024-07-02 | 国网四川省电力公司经济技术研究院 | 一种电气综合能源系统的协同优化运行方法及系统 |
CN114066151A (zh) * | 2021-10-12 | 2022-02-18 | 国网福建省电力有限公司经济技术研究院 | 基于共识变量的城市综合能源系统优化调度方法及终端 |
CN114374219B (zh) * | 2021-11-29 | 2023-09-15 | 山东大学 | 基于合作博弈的园区综合能源系统分布式优化方法及系统 |
CN115271267B (zh) * | 2022-09-27 | 2023-04-25 | 中国电力科学研究院有限公司 | 多区域电能与辅助服务联合出清方法、装置、设备及介质 |
CN115860241B (zh) * | 2022-12-12 | 2024-02-02 | 河海大学 | 基于salr的互联多区域综合能源系统分布式调度方法及装置 |
CN116805800B (zh) * | 2023-02-28 | 2024-05-07 | 南方电网科学研究院有限责任公司 | 区域互联电网调度方法、装置和非易失性存储设备 |
CN118246710A (zh) * | 2024-05-28 | 2024-06-25 | 山东大学 | 一种计及多区域互动策略的电-热系统调度方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103248043A (zh) * | 2013-04-28 | 2013-08-14 | 西南交通大学 | 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法 |
CN109217291A (zh) * | 2018-08-28 | 2019-01-15 | 南京理工大学 | 考虑削峰填谷模型的电气互联系统多目标优化方法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
CN110429663A (zh) * | 2019-07-18 | 2019-11-08 | 中国电力科学研究院有限公司 | 一种利用储能系统辅助电力调峰的调度方法及系统 |
CN110620402A (zh) * | 2019-10-21 | 2019-12-27 | 山东大学 | 基于分布式多场景的电-气混联系统规划运行联合优化方法及系统 |
CN110707755A (zh) * | 2019-10-11 | 2020-01-17 | 南京师范大学 | 一种考虑非理想通信条件下基于能量枢纽的综合能源系统超短期调度方法 |
CN110729766A (zh) * | 2019-09-27 | 2020-01-24 | 国网河北省电力有限公司 | 一种多区域电-气综合能源系统分布式鲁棒优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9568901B2 (en) * | 2012-08-27 | 2017-02-14 | Nec Corporation | Multi-objective energy management methods for micro-grids |
-
2020
- 2020-03-16 CN CN202010182634.9A patent/CN111342452B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103248043A (zh) * | 2013-04-28 | 2013-08-14 | 西南交通大学 | 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法 |
CN109217291A (zh) * | 2018-08-28 | 2019-01-15 | 南京理工大学 | 考虑削峰填谷模型的电气互联系统多目标优化方法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
CN110429663A (zh) * | 2019-07-18 | 2019-11-08 | 中国电力科学研究院有限公司 | 一种利用储能系统辅助电力调峰的调度方法及系统 |
CN110729766A (zh) * | 2019-09-27 | 2020-01-24 | 国网河北省电力有限公司 | 一种多区域电-气综合能源系统分布式鲁棒优化方法 |
CN110707755A (zh) * | 2019-10-11 | 2020-01-17 | 南京师范大学 | 一种考虑非理想通信条件下基于能量枢纽的综合能源系统超短期调度方法 |
CN110620402A (zh) * | 2019-10-21 | 2019-12-27 | 山东大学 | 基于分布式多场景的电-气混联系统规划运行联合优化方法及系统 |
Non-Patent Citations (1)
Title |
---|
"计及中长期合同电量分解和风电不确定性的电-气综合能源系统日前优化调度";吴刚,刘俊勇,向月,沈晓东,马宇航;《电力自动化设备》;第39卷(第08期);第246-253页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111342452A (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111342452B (zh) | 一种多区域电气综合能源系统能量与备用分散式调度方法 | |
CN111900732B (zh) | 考虑风电不确定性及无功潮流与电压约束的风电集群接入方式规划方法 | |
CN106786799B (zh) | 一种直流联络线功率阶梯化发电计划优化方法 | |
CN112701687B (zh) | 考虑价格型联合需求响应的气电配网系统鲁棒优化运行方法 | |
CN108599268B (zh) | 一种考虑风电场时空关联约束的日前优化调度方法 | |
CN111082466B (zh) | 考虑风电不确定性的新能源接入与网架扩建优化方法 | |
CN116341881B (zh) | 计及热网灵活性的电-热系统鲁棒超前调度方法及系统 | |
CN106786753A (zh) | 多用户的区域能源互联网的系统及其调节方法 | |
CN112990523A (zh) | 基于多目标模型预测控制的区域综合能源系统分层优化运行方法 | |
CN113792969A (zh) | 考虑气网动态特性和电转气综合能源系统优化调度方法 | |
CN114865631A (zh) | 源荷协同降碳的综合能源系统最优分布鲁棒经济调度方法 | |
CN114077934A (zh) | 一种综合能源微网互联系统及其调度方法 | |
CN115860413A (zh) | 基于负荷需求响应与双层可调鲁棒优化的并网冷热电联产微网系统经济调度方法 | |
CN112380681A (zh) | 基于分布式优化的综合能源系统参数确定方法及相关装置 | |
CN113762632B (zh) | 一种电气综合能源系统的协同优化运行方法及系统 | |
CN107330551A (zh) | 一种优化安排能源基地外送通道的送电方法 | |
CN111342453B (zh) | 计及多类型备用资源的电气综合能源系统备用决策方法 | |
CN104239960A (zh) | 考虑抽水蓄能机组的发电计划优化方法 | |
CN117477542A (zh) | 一种考虑灵活性约束的源储协同规划方法 | |
CN117455034A (zh) | 一种用于海上风电输电方案优选的编码优化方法及终端 | |
CN115021327B (zh) | 基于虚拟储能的互联微电网系统的协同日前优化调度方法 | |
CN115425697B (zh) | 基于交替方向乘子法的分布式跨区跨省调度方法及系统 | |
CN115719132A (zh) | 含多座抽蓄电站的省级电网日前计划分层递进优化方法 | |
CN115051406A (zh) | 一种计及源网不确定性的多阶段协调控制风险调度方法 | |
CN107196348A (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 |