CN111049135B - 一种多区域电气耦合系统的分布式两阶段合作运行方法 - Google Patents
一种多区域电气耦合系统的分布式两阶段合作运行方法 Download PDFInfo
- Publication number
- CN111049135B CN111049135B CN201911390588.5A CN201911390588A CN111049135B CN 111049135 B CN111049135 B CN 111049135B CN 201911390588 A CN201911390588 A CN 201911390588A CN 111049135 B CN111049135 B CN 111049135B
- Authority
- CN
- China
- Prior art keywords
- region
- node
- natural gas
- power
- area
- 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
- 230000008878 coupling Effects 0.000 title claims abstract description 76
- 238000010168 coupling process Methods 0.000 title claims abstract description 76
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000005457 optimization Methods 0.000 claims abstract description 86
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 210
- 239000003345 natural gas Substances 0.000 claims description 105
- 239000007789 gas Substances 0.000 claims description 48
- 238000012546 transfer Methods 0.000 claims description 21
- 239000002343 natural gas well Substances 0.000 claims description 14
- 238000010248 power generation Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 7
- 238000012887 quadratic function Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 5
- 239000003245 coal Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000033772 system development 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
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种多区域电气耦合系统的分布式两阶段合作运行方法,属于电‑气耦合系统运行和控制技术领域。本发明方法通过在优化过程中引入协调层,在优化过程中考虑各区域电气耦合系统合作优化运行后所降低成本的再分配,进一步得到合作优化运行下各区域内部的运行参数,为多区域电气耦合系统的分布式合作运行提供支持。本方法保证个体理性与集体理性一致,采用分布式方法进行求解,保护了区域隐私信息,对大型电气耦合系统的优化运行具有非常重要的意义。
Description
技术领域
本发明涉及一种多区域电气耦合系统的分布式两阶段合作运行方法,属于电-气耦合系统运行和控制技术领域。
背景技术
为提高能源总体使用效率,扩大可再生能源消纳能力,对多类能源互联集成和互补融合的需求增长迅速,综合能源系统成为能源系统发展的重要趋势之一。其中,电气耦合系统是综合能源系统的一个重要组成部分。
由于行政、地理等因素,一个大型的电气耦合系统往往由多个子区域电气耦合系统组成,区域之间存在有联络管道和联络线进行能量传输。为获得更大的社会效益,不同区域之间会考虑合作优化运行,而如何对新增社会福利进行再分配,保证个体理性与集体理性一致,成为亟待解决的重要问题。与此同时,由于各区域运营商相互独立,为保护区域隐私信息,集中优化的方难以适应大型电气耦合系统的运行需求,需要采用分布式方法进行求解。
因此,研究一种可以保护各子区域隐私信息,同时兼顾个体理性和集体理性的多区域电气耦合系统合作优化运行方法对大型电气耦合系统的优化运行具有非常重要的意义。
发明内容
本发明的目的是提出一种多区域电气耦合系统的分布式两阶段合作运行方法,在优化过程中自动考虑新增社会效益的再分配,兼顾了个体理性、社会理性以及用户隐私。
本发明提出的多区域电气耦合系统的分布式两阶段合作运行方法,包括以下步骤:
(1)根据电气耦合系统的联络线和联络管道,将一个大型的电气耦合系统分为S个区域,区域编号为s=1,2,…,S,不同的区域之间通过联络管道传输天然气,通过联络线传输电能;
(2)从S个区域中任取两个相连的区域,假设两个相连区域的编号分别为s1和s2,定义为两个相连区域的联络线与s1相连节点的相角在区域s1的取值,为两个相连区域的联络线与s1相连节点的相角在区域s1取值,为两个相连区域的联络管道中天然气流量在区域s1的取值,为两个相连区域的联络线与s2相连节点的相角在区域s2取值,为两个相连区域的联络线与s1相连节点的相角在区域s2的取值,为两个相连区域的联络管道中的天然气流量在区域s2的取值,设定相连两个区域上同一条联络线上的耦合变量满足:
同一条联络管道上的耦合变量满足:
(3)建立每个区域s的内部优化模型,内部优化模型的约束条件包括:
(3-1)区域s的电网约束包括:
(3-1-1)电网节点电量平衡约束:
其中,k为电网中的节点编号,Ne,s为区域s内的电网节点总数;上标cfu表示燃煤电站,Pi cfu表示燃煤电站i的发电有功功率,为待求变量,表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和;上标gfu表示燃气电站,表示燃气电站j的发电有功功率,为待求变量,表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和;上标ld表示电负荷,Pl ld表示负荷l的有功功率,Lsl表示负荷l的弃电有功功率,为待求变量,表示区域s中所有与电网节点k相连的负荷实际有功功率之和;pfmk表示区域s中节点m和节点k所连接支路的有功功率,为待求变量,规定从节点m流向节点k为正,反向为负;表示所有流入节点k的联络线支路的有功功率;
(3-1-2)电网直流潮流约束:
其中,θm、θk分别表示区域s中节点m和节点k的电压相角,xmk表示节点m和节点k间所连支路的电抗;
(3-1-3)电网中参考节点的电压相角约束:
θh=0,h∈refp
其中,θh表示电网节点h的相角,refp表示电网选取的参考节点集合;
(3-1-4)电网支路的有功功率上限和下限约束:
(3-1-5)电网燃煤电站的有功功率上限和下限约束:
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中,Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限,由电网调度给定;
(3-1-6)电网燃气电站有功功率上限和下限约束:
(3-1-7)电网弃负荷功率上限和下限约束:
(3-2)区域s的天然气网约束,包括:
(3-2-1)天然气网压力参考节点约束:
ωq=const,q=refg
其中,ωq表示天然气网节点q的压力,refg表示天然气网选取的参考节点集合;
(3-2-2)天然气网节点的天然气流量平衡约束:
其中,w为天然气网中的节点编号,Ng,s为区域s中的天然气网节点总数;上标well表示天然气井,表示天然气井d的出气流量,为待求变量,表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和;上标res表示居民负荷,表示居民负荷e消耗的天然气流量,由天然气网调度给定,表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和;gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量,为待求变量,规定从天然气网节点n流向节点w为正,反向为负;表示所有流入节点w的联络管道中的流量;
(3-2-3)天然气网节点压力的上限和下限约束:
(3-2-4)天然气网中天然气井出气流量的上限和下限约束:
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束:
其中,gfnw表示从天然气网节点n与节点w相连管道的天然气流量,规定从节点n流向节点w为正,反向为负;ωn、ωw分别表示节点n、节点w的压力;sgn(ωn,ωw)为关于ωn、ωw的函数,当ωn>ωw时,sgn(ωn,ωw)取1,当ωn≤ωw,sgn(ωn,ωw)取值为0;Cnw为从节点n、节点w所连管道的韦茅斯常数,由天然气网提供,由于流量和压力之间关系约束中的sgn(ωn,ωw)的取值为一个二元变量,引入整数变量将该约束转化为以下表达式:
(3-2-6)天然气网中天然气管道流量的上限约束:
(3-3)区域s的电网和天然气网之间的耦合约束如下:
(3-4)区域s独立优化运行时,区域s的联络线和联络管道上的约束为:
gflink,s=0
(3-5)区域s的内部优化模型的目标函数为:
每个区域进行独立优化运行时,目标函数表示为:
其中,OCs表示区域s的运行成本,表示区域s的电负荷l的弃负荷单位成本,表示天然气网的天然气井d的产气单位成本,η为罚函数因子,取值范围为1-100,表示电网中燃煤电站的i的发电单位成本,发电单位成本为发电量的二次函数,二次函数中的各项系数由燃煤电站给定:
(3-5)利用分支定界法和内点法,求解上述约束条件和目标函数组成的区域s的内部优化模型,得到各区域s独立优化时的优化成本OCs;
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型,包括:
(4-1)基于讨价还价的博弈模型,分布式两阶段合作运行优化模型的目标函数表示为:
其中,TPs为电气耦合系统中区域s的转移支付,即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配,OCs为各区域s单独运行的最小成本,为步骤(3-6)中所求解得到的OCs,TOCs表示区域电气耦合系统s的运行成本,表示为:
(4-2)分布式两阶段合作运行优化模型的约束条件为:
(4-2-1)多区域电气耦合系统各区域s,满足步骤(3-1)、(3-2)、(3-3)中的条件约束,以及与区域s相连的联络线的约束:
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量,满足步骤(2)中的条件约束;
(4-3)多区域电气耦合系统中的所有区域的转移支付满足:
(5)采用带协调层的两阶段交替方向乘子法(ADMM),求解上述步骤(4)的优化模型,包括以下步骤:
(5-1)建立一个信号处理协调层,使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯;
(5-3)对协调层进行初始化,设初始化迭代次数k=0,初始化协调层发送给各区域s一个协调变量向量的维数与步骤(5-2)的相同,初始化各区域s的拉格朗日乘子向量的维数与相同,惩罚因子向量的维数与相同,向量中元素取值范围为0.1-10,设定收敛阈值取值范围为0.001-0.1,并将初始化的传递给各区域s;
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型,>得到协调变量将该协调变量判断与收敛阈值进行比较,若则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本,相应各约束中的待求变量的取值即为设备运行情况,进行步骤(5-8),若则令
(5-8)对协调层进行第二次初始化,初始化迭代次数k=0,初始化所有区域s的拉格朗日乘子ρs=0,取值范围为0.1-10,收敛阈值εs,取值范围为0.001-0.1,初始化区域s的转移支付在协调层的取值并将耦合变量拉格朗日乘子发送给各区域s;
约束为:
若不等式不成立,则令
协调层将λi k+1发送给各区域s,令k=k+1,返回步骤(5-9);
(6)根据步骤(5-5)得到的各区域s约束中的待求变量取值,即区域s中各燃煤电站i的发电有功功率Pi cfu,各燃气电站j的发电有功功率各电网负荷l的弃电有功功率Lsl,各天然气井d的出气流量各天然气网节点q的压力ωq,作为各区域s中分布式合作优化运行的参数,步骤(5-11)的得到的TPs为各区域s的转移支付,实现多区域电气耦合系统的分布式两阶段合作运行。
本发明提出的一种多区域电气耦合系统的分布式两阶段合作运行方法,其优点是:
本发明方法通过在优化过程中引入协调层,在优化过程中考虑各区域电气耦合系统合作优化运行后所降低成本的再分配,进一步得到合作优化运行下各区域内部的运行参数,为多区域电气耦合系统的分布式合作运行提供支持。本方法保证个体理性与集体理性一致,采用分布式方法进行求解,保护了区域隐私信息,对大型电气耦合系统的优化运行具有非常重要的意义。
附图说明
图1为本发明方法涉及的多区域电气耦合系统与协调层的关系示意图。
具体实施方式
本发明提出的多区域电气耦合系统的分布式两阶段合作运行方法,包括以下步骤:
(1)根据电气耦合系统的联络线和联络管道,将一个大型的电气耦合系统分为S个区域,区域编号为s=1,2,…,S,不同的区域之间通过联络管道传输天然气,通过联络线传输电能,如图1所示;
(2)从S个区域中任取两个相连的区域,假设两个相连区域的编号分别为s1和s2,定义为两个相连区域的联络线与s1相连节点的相角在区域s1的取值,为两个相连区域的联络线与s1相连节点的相角在区域s1取值,为两个相连区域的联络管道中天然气流量在区域s1的取值,为两个相连区域的联络线与s2相连节点的相角在区域s2取值,为两个相连区域的联络线与s1相连节点的相角在区域s2的取值,为两个相连区域的联络管道中的天然气流量在区域s2的取值,设定相连两个区域上同一条联络线上的耦合变量满足:
同一条联络管道上的耦合变量满足:
(3)建立每个区域s的内部优化模型,内部优化模型的约束条件包括:
(3-1)区域s的电网约束包括:
(3-1-1)电网节点电量平衡约束:
其中,k为电网中的节点编号,Ne,s为区域s内的电网节点总数;上标cfu表示燃煤电站,Pi cfu表示燃煤电站i的发电有功功率,为待求变量,表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和;上标gfu表示燃气电站,表示燃气电站j的发电有功功率,为待求变量,表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和;上标ld表示电负荷,Pl ld表示负荷l的有功功率,Lsl表示负荷l的弃电有功功率,为待求变量,表示区域s中所有与电网节点k相连的负荷实际有功功率之和;pfmk表示区域s中节点m和节点k所连接支路的有功功率,为待求变量,规定从节点m流向节点k为正,反向为负;表示所有流入节点k的联络线支路的有功功率;
(3-1-2)电网直流潮流约束:
其中,θm、θk分别表示区域s中节点m和节点k的电压相角,xmk表示节点m和节点k间所连支路的电抗;
(3-1-3)电网中参考节点的电压相角约束:
θh=0,h∈refp
其中,θh表示电网节点h的相角,refp表示电网选取的参考节点集合;
(3-1-4)电网支路的有功功率上限和下限约束:
(3-1-5)电网燃煤电站的有功功率上限和下限约束:
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中,Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限,由电网调度给定;
(3-1-6)电网燃气电站有功功率上限和下限约束:
(3-1-7)电网弃负荷功率上限和下限约束:
(3-2)区域s的天然气网约束,包括:
(3-2-1)天然气网压力参考节点约束:
ωq=const,q=refg
其中,ωq表示天然气网节点q的压力,refg表示天然气网选取的参考节点集合;
(3-2-2)天然气网节点的天然气流量平衡约束:
其中,w为天然气网中的节点编号,Ng,s为区域s中的天然气网节点总数;上标well表示天然气井,表示天然气井d的出气流量,为待求变量,表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和;上标res表示居民负荷,表示居民负荷e消耗的天然气流量,由天然气网调度给定,表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和;gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量,为待求变量,规定从天然气网节点n流向节点w为正,反向为负;表示所有流入节点w的联络管道中的流量;
(3-2-3)天然气网节点压力的上限和下限约束:
(3-2-4)天然气网中天然气井出气流量的上限和下限约束:
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束:
其中,gfnw表示从天然气网节点n与节点w相连管道的天然气流量,规定从节点n流向节点w为正,反向为负;ωn、ωw分别表示节点n、节点w的压力;sgn(ωn,ωw)为关于ωn、ωw的函数,当ωn>ωw时,sgn(ωn,ωw)取1,当ωn≤ωw,sgn(ωn,ωw)取值为0;Cnw为从节点n、节点w所连管道的韦茅斯常数,由天然气网提供,由于流量和压力之间关系约束中的sgn(ωn,ωw)的取值为一个二元变量,引入整数变量将该约束转化为以下表达式:
(3-2-6)天然气网中天然气管道流量的上限约束:
(3-3)区域s的电网和天然气网之间的耦合约束如下:
(3-4)区域s独立优化运行时,区域s的联络线和联络管道上的约束为:
gflink,s=0
(3-5)区域s的内部优化模型的目标函数为:
每个区域进行独立优化运行时,目标函数表示为:
其中,OCs表示区域s的运行成本,表示区域s的电负荷l的弃负荷单位成本,表示天然气网的天然气井d的产气单位成本,η为罚函数因子,取值范围为1-100,本发明一个实施例中取值10,表示电网中燃煤电站的i的发电单位成本,发电单位成本为发电量的二次函数,二次函数中的各项系数由燃煤电站给定:
(3-5)利用分支定界法和内点法,求解上述约束条件和目标函数组成的区域s的内部优化模型,得到各区域s独立优化时的优化成本OCs;
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型,包括:
(4-1)基于讨价还价的博弈模型,分布式两阶段合作运行优化模型的目标函数表示为:
其中,TPs为电气耦合系统中区域s的转移支付,即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配,OCs为各区域s单独运行的最小成本,为步骤(3-6)中所求解得到的OCs,TOCs表示区域电气耦合系统s的运行成本,表示为:
(4-2)分布式两阶段合作运行优化模型的约束条件为:
(4-2-1)多区域电气耦合系统各区域s,满足步骤(3-1)、(3-2)、(3-3)中的条件约束,以及与区域s相连的联络线的约束:
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量,满足步骤(2)中的条件约束;
(4-3)多区域电气耦合系统中的所有区域的转移支付满足:
(5)采用带协调层的两阶段交替方向乘子法(ADMM),求解上述步骤(4)的优化模型,包括以下步骤:
(5-1)建立一个信号处理协调层,使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯;
(5-3)对协调层进行初始化,设初始化迭代次数k=0,初始化协调层发送给各区域s一个协调变量向量的维数与步骤(5-2)的相同,初始化各区域s的拉格朗日乘子向量的维数与相同,惩罚因子向量的维数与相同,向量中元素取值范围为0.1-10,本发明一个实施例中所有元素取值1,设定收敛阈值取值范围为0.001-0.1,本发明一个实施例中取值0.01,并将初始化的传递给各区域s;
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型,>得到协调变量将该协调变量判断与收敛阈值进行比较,若则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本,相应各约束中的待求变量的取值即为设备运行情况,进行步骤(5-8),若则令
(5-8)对协调层进行第二次初始化,初始化迭代次数k=0,初始化所有区域s的拉格朗日乘子ρs=0,取值范围为0.1-10,本发明一个实施例中取值1,收敛阈值εs,取值范围为0.001-0.1,本发明一个实施例中取值0.01,初始化区域s的转移支付在协调层的取值并将耦合变量拉格朗日乘子发送给各区域s;
约束为:
若不等式不成立,则令
协调层将λi k+1发送给各区域s,令k=k+1,返回步骤(5-9);
Claims (1)
1.一种多区域电气耦合系统的分布式两阶段合作运行方法,其特征在于该方法包括以下步骤:
(1)根据电气耦合系统的联络线和联络管道,将一个大型的电气耦合系统分为S个区域,区域编号为s=1,2,…,S,不同的区域之间通过联络管道传输天然气,通过联络线传输电能;
(2)从S个区域中任取两个相连的区域,假设两个相连区域的编号分别为s1和s2,定义为两个相连区域的联络线与s1相连节点的相角在区域s1的取值,为两个相连区域的联络线与s1相连节点的相角在区域s1取值,为两个相连区域的联络管道中天然气流量在区域s1的取值,为两个相连区域的联络线与s2相连节点的相角在区域s2取值,为两个相连区域的联络线与s1相连节点的相角在区域s2的取值,为两个相连区域的联络管道中的天然气流量在区域s2的取值,设定相连两个区域上同一条联络线上的耦合变量满足:
同一条联络管道上的耦合变量满足:
(3)建立每个区域s的内部优化模型,内部优化模型的约束条件包括:
(3-1)区域s的电网约束包括:
(3-1-1)电网节点电量平衡约束:
其中,k为电网中的节点编号,Ne,s为区域s内的电网节点总数;上标cfu表示燃煤电站,Pi cfu表示燃煤电站i的发电有功功率,为待求变量,表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和;上标gfu表示燃气电站,表示燃气电站j的发电有功功率,为待求变量,表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和;上标ld表示电负荷,Pl ld表示负荷l的有功功率,Lsl表示负荷l的弃电有功功率,为待求变量,表示区域s中所有与电网节点k相连的负荷实际有功功率之和;pfmk表示区域s中节点m和节点k所连接支路的有功功率,为待求变量,规定从节点m流向节点k为正,反向为负;表示所有流入节点k的联络线支路的有功功率;
(3-1-2)电网直流潮流约束:
其中,θm、θk分别表示区域s中节点m和节点k的电压相角,xmk表示节点m和节点k间所连支路的电抗;
(3-1-3)电网中参考节点的电压相角约束:
θh=0,h∈refp
其中,θh表示电网节点h的相角,refp表示电网选取的参考节点集合;
(3-1-4)电网支路的有功功率上限和下限约束:
(3-1-5)电网燃煤电站的有功功率上限和下限约束:
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中,Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限,由电网调度给定;
(3-1-6)电网燃气电站有功功率上限和下限约束:
(3-1-7)电网弃负荷功率上限和下限约束:
(3-2)区域s的天然气网约束,包括:
(3-2-1)天然气网压力参考节点约束:
ωq=const,q=refg
其中,ωq表示天然气网节点q的压力,refg表示天然气网选取的参考节点集合;
(3-2-2)天然气网节点的天然气流量平衡约束:
其中,w为天然气网中的节点编号,Ng,s为区域s中的天然气网节点总数;上标well表示天然气井,表示天然气井d的出气流量,为待求变量,表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和;上标res表示居民负荷,表示居民负荷e消耗的天然气流量,由天然气网调度给定,表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和;gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量,为待求变量,规定从天然气网节点n流向节点w为正,反向为负;表示所有流入节点w的联络管道中的流量;
(3-2-3)天然气网节点压力的上限和下限约束:
(3-2-4)天然气网中天然气井出气流量的上限和下限约束:
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束:
其中,gfnw表示从天然气网节点n与节点w相连管道的天然气流量,规定从节点n流向节点w为正,反向为负;ωn、ωw分别表示节点n、节点w的压力;sgn(ωn,ωw)为关于ωn、ωw的函数,当ωn>ωw时,sgn(ωn,ωw)取1,当ωn≤ωw,sgn(ωn,ωw)取值为0;Cnw为从节点n、节点w所连管道的韦茅斯常数,由天然气网提供,由于流量和压力之间关系约束中的sgn(ωn,ωw)的取值为一个二元变量,引入整数变量将该约束转化为以下表达式:
(3-2-6)天然气网中天然气管道流量的上限约束:
(3-3)区域s的电网和天然气网之间的耦合约束如下:
(3-4)区域s独立优化运行时,区域s的联络线和联络管道上的约束为:
gflink,s=0
(3-5)区域s的内部优化模型的目标函数为:
每个区域进行独立优化运行时,目标函数表示为:
其中,OCs表示区域s的运行成本,表示区域s的电负荷l的弃负荷单位成本,表示天然气网的天然气井d的产气单位成本,η为罚函数因子,取值范围为1-100,表示电网中燃煤电站i的发电单位成本,发电单位成本为发电量的二次函数,二次函数中的各项系数由燃煤电站给定:
(3-6)利用分支定界法和内点法,求解上述约束条件和目标函数组成的区域s的内部优化模型,得到各区域s独立优化时的优化成本OCs;
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型,包括:
(4-1)基于讨价还价的博弈模型,分布式两阶段合作运行优化模型的目标函数表示为:
其中,TPs为电气耦合系统中区域s的转移支付,即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配,OCs为各区域s单独运行的最小成本,为步骤(3-6)中所求解得到的OCs,TOCs表示区域电气耦合系统s的运行成本,表示为:
(4-2)分布式两阶段合作运行优化模型的约束条件为:
(4-2-1)多区域电气耦合系统各区域s,满足步骤(3-1)、(3-2)、(3-3)中的条件约束,以及与区域s相连的联络线的约束:
其中,θf,s为联络线的一个节点与s区域相连节点的相角在区域s的取值,θt,s为该联络线另一节点的相角在区域s的取值;
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量,满足步骤(2)中的条件约束;
(4-3)多区域电气耦合系统中的所有区域的转移支付满足:
(5)采用带协调层的两阶段交替方向乘子法(ADMM),求解上述步骤(4)的优化模型,包括以下步骤:
(5-1)建立一个信号处理协调层,使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯;
(5-3)对协调层进行初始化,设初始化迭代次数k=0,初始化协调层发送给各区域s一个协调变量向量的维数与步骤(5-2)的相同,初始化各区域s的拉格朗日乘子向量的维数与相同,惩罚因子向量的维数与相同,向量中元素取值范围为0.1-10,设定收敛阈值取值范围为0.001-0.1,并将初始化的 传递给各区域s;
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型,>得到协调变量将该协调变量判断与收敛阈值进行比较,若则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本,相应各约束中的待求变量的取值即为设备运行情况,进行步骤(5-8),若则令
(5-8)对协调层进行第二次初始化,初始化迭代次数k=0,初始化所有区域s的拉格朗日乘子ρs=0,取值范围为0.1-10,收敛阈值εs,取值范围为0.001-0.1,初始化区域s的转移支付在协调层的取值并将耦合变量拉格朗日乘子发送给各区域s;
约束为:
若不等式不成立,则令
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911390588.5A CN111049135B (zh) | 2019-12-30 | 2019-12-30 | 一种多区域电气耦合系统的分布式两阶段合作运行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911390588.5A CN111049135B (zh) | 2019-12-30 | 2019-12-30 | 一种多区域电气耦合系统的分布式两阶段合作运行方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111049135A CN111049135A (zh) | 2020-04-21 |
CN111049135B true CN111049135B (zh) | 2023-02-03 |
Family
ID=70241445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911390588.5A Active CN111049135B (zh) | 2019-12-30 | 2019-12-30 | 一种多区域电气耦合系统的分布式两阶段合作运行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111049135B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113642165B (zh) * | 2021-07-29 | 2024-04-12 | 西安理工大学 | 考虑需求侧不确定性的电-气耦合配网的分布式优化方法 |
CN118281961B (zh) * | 2024-06-04 | 2024-08-09 | 太原理工大学 | 一种电-热耦合系统的分布式协同机组组合方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2869377A1 (en) * | 2013-10-29 | 2015-05-06 | Total Marketing Services | Power generation system |
CN106096269A (zh) * | 2016-06-12 | 2016-11-09 | 清华大学 | 一种电‑气耦合系统中天然气网的区间潮流计算方法 |
CN106960282A (zh) * | 2017-03-27 | 2017-07-18 | 武汉大学 | 一种气电综合能源配网系统的协调运行方法 |
CN107800158A (zh) * | 2017-10-17 | 2018-03-13 | 清华大学 | 一种兼顾经济与能效的电‑热耦合多能流系统优化调度方法 |
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法 |
CN110046750A (zh) * | 2019-03-25 | 2019-07-23 | 国网江苏省电力有限公司经济技术研究院 | 一种能源互联网协同优化运行方法 |
CN110503250A (zh) * | 2019-08-08 | 2019-11-26 | 燕山大学 | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150263589A1 (en) * | 2014-03-11 | 2015-09-17 | Varnell M. Castor | Rail Barrel Direct Energy Transferor Piezoelectricity (RBDETP) |
-
2019
- 2019-12-30 CN CN201911390588.5A patent/CN111049135B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2869377A1 (en) * | 2013-10-29 | 2015-05-06 | Total Marketing Services | Power generation system |
CN106096269A (zh) * | 2016-06-12 | 2016-11-09 | 清华大学 | 一种电‑气耦合系统中天然气网的区间潮流计算方法 |
CN106960282A (zh) * | 2017-03-27 | 2017-07-18 | 武汉大学 | 一种气电综合能源配网系统的协调运行方法 |
CN107800158A (zh) * | 2017-10-17 | 2018-03-13 | 清华大学 | 一种兼顾经济与能效的电‑热耦合多能流系统优化调度方法 |
CN110046750A (zh) * | 2019-03-25 | 2019-07-23 | 国网江苏省电力有限公司经济技术研究院 | 一种能源互联网协同优化运行方法 |
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法 |
CN110503250A (zh) * | 2019-08-08 | 2019-11-26 | 燕山大学 | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 |
Non-Patent Citations (2)
Title |
---|
考虑不确定性及电/热储能的综合能源系统两阶段规划-运行联合优化方法;雷金勇等;《电力自动化设备》;20190815;第39卷(第8期);第169-175页 * |
考虑网络约束和源荷不确定性的区域综合能源系统两阶段鲁棒优化调度;郭尊等;《电网技术》;20190915;第43卷(第9期);第3090-3100页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111049135A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers | |
AU2020100429A4 (en) | A dynamic optimal energy flow computing method for the combined heat and power system | |
Lv et al. | Model predictive control based robust scheduling of community integrated energy system with operational flexibility | |
Singh et al. | Optimal scheduling of water distribution systems | |
CN110348602B (zh) | 计及天然气管网和热力管网特性的综合能源系统优化方法 | |
CN111049135B (zh) | 一种多区域电气耦合系统的分布式两阶段合作运行方法 | |
CN107968439B (zh) | 基于混合整数线性规划的主动配电网联合优化算法 | |
CN111582658A (zh) | 考虑多区域互联协同的能源站设备配置及管线规划方法 | |
CN112668188B (zh) | 多园区综合能源系统的分布式鲁棒协同优化调度方法 | |
CN112861357B (zh) | 一种多微能网耦合系统分布式协同调度方法及系统 | |
CN106055773A (zh) | 结合热网模型的多区域综合能源系统配置模型的建立方法 | |
CN106992555B (zh) | 一种辅助服务市场环境下的热电联合系统分布式调度方法 | |
US11011909B1 (en) | Feedforward dynamic and distributed energy storage system | |
CN109389248B (zh) | 基于区域热网的综合能源系统热电功率协调调度方法 | |
CN108053104A (zh) | 区域电力-天然气-热力综合能源系统能量流优化方法 | |
CN110084394A (zh) | 一种考虑温控设备和电动汽车的工业用户智能用电优化方法 | |
US20230318300A1 (en) | Feedforward Dynamic Multifunctional Structure with Distributed Energy Storage System | |
CN118052420B (zh) | 计及多供热网络互动策略的电-热系统调度方法及系统 | |
CN116542447A (zh) | 一种电热系统优化调度方法 | |
CN105956693A (zh) | 一种基于分布式梯度算法的虚拟电厂经济调度方法 | |
CN111931977B (zh) | 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法 | |
CN113837589A (zh) | 一种综合能源系统n-1调度方法及装置 | |
CN112531788A (zh) | 考虑多重不确定性和自趋优运行的透明微网群规划方法 | |
CN112200695A (zh) | 一种城市级综合能源系统优化调度方法 | |
CN115222102B (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 |