CN111049135B - 一种多区域电气耦合系统的分布式两阶段合作运行方法 - Google Patents

一种多区域电气耦合系统的分布式两阶段合作运行方法 Download PDF

Info

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
Application number
CN201911390588.5A
Other languages
English (en)
Other versions
CN111049135A (zh
Inventor
王彬
孙勇
郭庆来
李振元
孙宏斌
曲绍杰
潘昭光
张磊
尹冠雄
陈彬彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
State Grid Jilin Electric Power Corp
Original Assignee
Tsinghua University
State Grid Jilin Electric Power Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University, State Grid Jilin Electric Power Corp filed Critical Tsinghua University
Priority to CN201911390588.5A priority Critical patent/CN111049135B/zh
Publication of CN111049135A publication Critical patent/CN111049135A/zh
Application granted granted Critical
Publication of CN111049135B publication Critical patent/CN111049135B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)
  • Entrepreneurship & Innovation (AREA)
  • Power Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种多区域电气耦合系统的分布式两阶段合作运行方法，属于电‑气耦合系统运行和控制技术领域。本发明方法通过在优化过程中引入协调层，在优化过程中考虑各区域电气耦合系统合作优化运行后所降低成本的再分配，进一步得到合作优化运行下各区域内部的运行参数，为多区域电气耦合系统的分布式合作运行提供支持。本方法保证个体理性与集体理性一致，采用分布式方法进行求解，保护了区域隐私信息，对大型电气耦合系统的优化运行具有非常重要的意义。

Description

一种多区域电气耦合系统的分布式两阶段合作运行方法
技术领域
本发明涉及一种多区域电气耦合系统的分布式两阶段合作运行方法，属于电-气耦合系统运行和控制技术领域。
背景技术
为提高能源总体使用效率，扩大可再生能源消纳能力，对多类能源互联集成和互补融合的需求增长迅速，综合能源系统成为能源系统发展的重要趋势之一。其中，电气耦合系统是综合能源系统的一个重要组成部分。
由于行政、地理等因素，一个大型的电气耦合系统往往由多个子区域电气耦合系统组成，区域之间存在有联络管道和联络线进行能量传输。为获得更大的社会效益，不同区域之间会考虑合作优化运行，而如何对新增社会福利进行再分配，保证个体理性与集体理性一致，成为亟待解决的重要问题。与此同时，由于各区域运营商相互独立，为保护区域隐私信息，集中优化的方难以适应大型电气耦合系统的运行需求，需要采用分布式方法进行求解。
因此，研究一种可以保护各子区域隐私信息，同时兼顾个体理性和集体理性的多区域电气耦合系统合作优化运行方法对大型电气耦合系统的优化运行具有非常重要的意义。
发明内容
本发明的目的是提出一种多区域电气耦合系统的分布式两阶段合作运行方法，在优化过程中自动考虑新增社会效益的再分配，兼顾了个体理性、社会理性以及用户隐私。
本发明提出的多区域电气耦合系统的分布式两阶段合作运行方法，包括以下步骤：
(1)根据电气耦合系统的联络线和联络管道，将一个大型的电气耦合系统分为S个区域，区域编号为s＝1,2,…,S，不同的区域之间通过联络管道传输天然气，通过联络线传输电能；
(2)从S个区域中任取两个相连的区域，假设两个相连区域的编号分别为s1和s2，定义
Figure BDA00023448336500000210
为两个相连区域的联络线与s1相连节点的相角在区域s1的取值，
Figure BDA00023448336500000211
为两个相连区域的联络线与s1相连节点的相角在区域s1取值，
Figure BDA00023448336500000212
为两个相连区域的联络管道中天然气流量在区域s1的取值，
Figure BDA00023448336500000213
为两个相连区域的联络线与s2相连节点的相角在区域s2取值，
Figure BDA00023448336500000214
为两个相连区域的联络线与s1相连节点的相角在区域s2的取值，
Figure BDA00023448336500000215
为两个相连区域的联络管道中的天然气流量在区域s2的取值，设定相连两个区域上同一条联络线上的耦合变量满足：
Figure BDA0002344833650000021
Figure BDA0002344833650000022
同一条联络管道上的耦合变量满足：
Figure BDA0002344833650000023
(3)建立每个区域s的内部优化模型，内部优化模型的约束条件包括：
(3-1)区域s的电网约束包括：
(3-1-1)电网节点电量平衡约束：
Figure BDA0002344833650000024
其中，k为电网中的节点编号，Ne,s为区域s内的电网节点总数；上标cfu表示燃煤电站，Pi cfu表示燃煤电站i的发电有功功率，为待求变量，
Figure BDA0002344833650000025
表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和；上标gfu表示燃气电站，
Figure BDA00023448336500000216
表示燃气电站j的发电有功功率，为待求变量，
Figure BDA0002344833650000026
表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和；上标ld表示电负荷，Pl ld表示负荷l的有功功率，Lsl表示负荷l的弃电有功功率，为待求变量，
Figure BDA0002344833650000027
表示区域s中所有与电网节点k相连的负荷实际有功功率之和；pfmk表示区域s中节点m和节点k所连接支路的有功功率，为待求变量，规定从节点m流向节点k为正，反向为负；
Figure BDA0002344833650000028
表示所有流入节点k的联络线支路的有功功率；
(3-1-2)电网直流潮流约束：
Figure BDA0002344833650000029
其中，θm、θk分别表示区域s中节点m和节点k的电压相角，xmk表示节点m和节点k间所连支路的电抗；
(3-1-3)电网中参考节点的电压相角约束：
θh＝0,h∈refp
其中，θh表示电网节点h的相角，refp表示电网选取的参考节点集合；
(3-1-4)电网支路的有功功率上限和下限约束：
Figure BDA0002344833650000031
其中，
Figure BDA0002344833650000034
表示区域s中节点m和节点k间所连支路的有功功率上限，由电网调度给定；
(3-1-5)电网燃煤电站的有功功率上限和下限约束：
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中，Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限，由电网调度给定；
(3-1-6)电网燃气电站有功功率上限和下限约束：
Figure BDA0002344833650000032
其中，
Figure BDA0002344833650000035
分别表示燃气电站j的发电有功功率下限和上限，由电网调度给定；
(3-1-7)电网弃负荷功率上限和下限约束：
Figure BDA0002344833650000033
其中，
Figure BDA0002344833650000036
分别表示区域s中电负荷l的允许弃负荷功率上限，由电网调度给定；
(3-2)区域s的天然气网约束，包括：
(3-2-1)天然气网压力参考节点约束：
ωq＝const,q＝refg
其中，ωq表示天然气网节点q的压力，refg表示天然气网选取的参考节点集合；
(3-2-2)天然气网节点的天然气流量平衡约束：
Figure BDA0002344833650000041
其中，w为天然气网中的节点编号，Ng,s为区域s中的天然气网节点总数；上标well表示天然气井，
Figure BDA0002344833650000049
表示天然气井d的出气流量，为待求变量，
Figure BDA0002344833650000042
表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和；上标res表示居民负荷，
Figure BDA00023448336500000410
表示居民负荷e消耗的天然气流量，由天然气网调度给定，
Figure BDA0002344833650000043
表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和；gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量，为待求变量，规定从天然气网节点n流向节点w为正，反向为负；
Figure BDA0002344833650000044
表示所有流入节点w的联络管道中的流量；
(3-2-3)天然气网节点压力的上限和下限约束：
Figure BDA0002344833650000045
其中，
Figure BDA00023448336500000411
分别表示天然气网节点q的压力下限和上限，由天然气网调度给定；
(3-2-4)天然气网中天然气井出气流量的上限和下限约束：
Figure BDA0002344833650000046
其中，
Figure BDA00023448336500000412
分别表示区域s中天然气井d的出气流量的下限和上限，由天然气网调度给定；
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束：
Figure BDA0002344833650000047
其中，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)的取值为一个二元变量，引入整数变量
Figure BDA00023448336500000413
将该约束转化为以下表达式：
Figure BDA0002344833650000048
Figure BDA0002344833650000051
Figure BDA0002344833650000052
Figure BDA0002344833650000053
Figure BDA0002344833650000054
Figure BDA0002344833650000055
其中，将上述约束
Figure BDA0002344833650000056
进一步松弛为：
Figure BDA0002344833650000057
(3-2-6)天然气网中天然气管道流量的上限约束：
Figure BDA0002344833650000058
其中，
Figure BDA00023448336500000512
表示从节点n和节点w的连接管道的流量上限，由天然气网调度给定；
(3-3)区域s的电网和天然气网之间的耦合约束如下：
Figure BDA0002344833650000059
其中，
Figure BDA00023448336500000513
和
Figure BDA00023448336500000514
为区域s中的燃气电站s的天然气发电量和耗气量二次关系的二次项系数、一次项系数、常数项，由燃气电站给出；
(3-4)区域s独立优化运行时，区域s的联络线和联络管道上的约束为：
Figure BDA00023448336500000510
gflink,s＝0
(3-5)区域s的内部优化模型的目标函数为：
每个区域进行独立优化运行时，目标函数表示为：
Figure BDA00023448336500000511
其中，OCs表示区域s的运行成本，
Figure BDA00023448336500000517
表示区域s的电负荷l的弃负荷单位成本，
Figure BDA00023448336500000515
表示天然气网的天然气井d的产气单位成本，η为罚函数因子，取值范围为1-100，
Figure BDA00023448336500000516
表示电网中燃煤电站的i的发电单位成本，发电单位成本为发电量的二次函数，二次函数中的各项系数
Figure BDA0002344833650000066
由燃煤电站给定：
Figure BDA0002344833650000061
(3-5)利用分支定界法和内点法，求解上述约束条件和目标函数组成的区域s的内部优化模型，得到各区域s独立优化时的优化成本OCs；
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型，包括：
(4-1)基于讨价还价的博弈模型，分布式两阶段合作运行优化模型的目标函数表示为：
Figure BDA0002344833650000062
其中，TPs为电气耦合系统中区域s的转移支付，即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配，OCs为各区域s单独运行的最小成本，为步骤(3-6)中所求解得到的OCs，TOCs表示区域电气耦合系统s的运行成本，表示为：
Figure BDA0002344833650000063
(4-2)分布式两阶段合作运行优化模型的约束条件为：
(4-2-1)多区域电气耦合系统各区域s，满足步骤(3-1)、(3-2)、(3-3)中的条件约束，以及与区域s相连的联络线的约束：
Figure BDA0002344833650000064
其中，θf,s为联络线的一个节点与s区域相连节点的相角在区域s的取值，
Figure BDA0002344833650000067
为该联络线另一节点的相角在区域s的取值；
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量，满足步骤(2)中的条件约束；
(4-3)多区域电气耦合系统中的所有区域的转移支付满足：
Figure BDA0002344833650000065
(5)采用带协调层的两阶段交替方向乘子法(ADMM)，求解上述步骤(4)的优化模型，包括以下步骤：
(5-1)建立一个信号处理协调层，使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯；
(5-2)多区域电气耦合系统中的各区域s分别同时向信号处理协调层发送一个协调向量
Figure BDA0002344833650000075
协调向量
Figure BDA0002344833650000076
中包括所有相应区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量；
(5-3)对协调层进行初始化，设初始化迭代次数k＝0，初始化协调层发送给各区域s一个协调变量
Figure BDA0002344833650000077
向量
Figure BDA0002344833650000078
的维数与步骤(5-2)的
Figure BDA0002344833650000079
相同，初始化各区域s的拉格朗日乘子
Figure BDA00023448336500000710
向量
Figure BDA00023448336500000711
的维数与
Figure BDA00023448336500000712
相同，惩罚因子
Figure BDA00023448336500000713
向量
Figure BDA00023448336500000714
的维数与
Figure BDA00023448336500000715
相同，向量中元素取值范围为0.1-10，设定收敛阈值
Figure BDA00023448336500000716
取值范围为0.001-0.1，并将初始化的
Figure BDA00023448336500000717
传递给各区域s；
(5-4)各区域s收到协调层初始化的
Figure BDA00023448336500000718
后，建立本区域的第一阶段优化模型，该第一阶段优化模型的约束为步骤(4-2-1)的条件约束，第一阶段优化模型目标函数为：
Figure BDA0002344833650000071
其中，
Figure BDA0002344833650000072
分别表示
Figure BDA0002344833650000073
的第v个元素；
(5-5)采用分支定界法和内点法，各区域s求解步骤(5-4)建立的第一阶段优化模型，得到各区域s合作运行的运行成本TOCs、协调变量
Figure BDA00023448336500000719
以及各约束中的待求变量，并将
Figure BDA00023448336500000720
发送回协调层；
(5-6)协调层收到所有区域s的协调变量
Figure BDA00023448336500000721
后，建立协调层第一阶段优化模型，协调层第一阶段优化模型的目标函数为：
Figure BDA0002344833650000074
不同区域s的协调变量
Figure BDA00023448336500000722
包括所有区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量，相邻区域联络管道和联络线上的相角和流量满足步骤(2)中的条件约束；
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型，>得到协调变量
Figure BDA00023448336500000723
将该协调变量
Figure BDA0002344833650000086
判断与收敛阈值
Figure BDA0002344833650000087
进行比较，若
Figure BDA0002344833650000088
则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本，相应各约束中的待求变量的取值即为设备运行情况，进行步骤(5-8)，若
Figure BDA0002344833650000089
则令
Figure BDA0002344833650000081
协调层将
Figure BDA00023448336500000810
发送给区域s，令k＝k+1，返回步骤(5-4)；
(5-8)对协调层进行第二次初始化，初始化迭代次数k＝0，初始化所有区域s的拉格朗日乘子
Figure BDA00023448336500000811
ρs＝0，取值范围为0.1-10，收敛阈值εs，取值范围为0.001-0.1，初始化区域s的转移支付在协调层的取值
Figure BDA00023448336500000812
并将耦合变量
Figure BDA00023448336500000813
拉格朗日乘子
Figure BDA00023448336500000814
发送给各区域s；
(5-9)各区域s收到协调层的耦合变量
Figure BDA00023448336500000815
后，建立本区域的第二阶段优化模型，该第二优化模型无条件约束，第二优化模型的目标函数为：
Figure BDA0002344833650000082
其中，TOCs为各区域s的运行成本，由步骤(5-5)得到，
Figure BDA00023448336500000816
表示区域s的转移支付在区域s的取值，为待求变量。各区域s求解本区域的第二阶段优化模型，得到转移支付
Figure BDA00023448336500000817
并将
Figure BDA00023448336500000818
发送回协调层；
(5-10)协调层收到
Figure BDA00023448336500000819
后，建立协调层的第二阶段优化模型，第二阶段目标函数为：
Figure BDA0002344833650000083
约束为：
Figure BDA0002344833650000084
(5-11)协调层采用内点法，求解步骤(5-10)建立的优化模型，得到
Figure BDA00023448336500000820
判断下面不等式是否成立：
Figure BDA0002344833650000085
若不等式成立，则迭代终止，
Figure BDA0002344833650000092
为最后各系统的转移支付，记为TPs，
若不等式不成立，则令
Figure BDA0002344833650000091
协调层将λi k+1发送给各区域s，令k＝k+1，返回步骤(5-9)；
(6)根据步骤(5-5)得到的各区域s约束中的待求变量取值，即区域s中各燃煤电站i的发电有功功率Pi cfu，各燃气电站j的发电有功功率
Figure BDA0002344833650000093
各电网负荷l的弃电有功功率Lsl，各天然气井d的出气流量
Figure BDA0002344833650000094
各天然气网节点q的压力ωq，作为各区域s中分布式合作优化运行的参数，步骤(5-11)的得到的TPs为各区域s的转移支付，实现多区域电气耦合系统的分布式两阶段合作运行。
本发明提出的一种多区域电气耦合系统的分布式两阶段合作运行方法，其优点是：
本发明方法通过在优化过程中引入协调层，在优化过程中考虑各区域电气耦合系统合作优化运行后所降低成本的再分配，进一步得到合作优化运行下各区域内部的运行参数，为多区域电气耦合系统的分布式合作运行提供支持。本方法保证个体理性与集体理性一致，采用分布式方法进行求解，保护了区域隐私信息，对大型电气耦合系统的优化运行具有非常重要的意义。
附图说明
图1为本发明方法涉及的多区域电气耦合系统与协调层的关系示意图。
具体实施方式
本发明提出的多区域电气耦合系统的分布式两阶段合作运行方法，包括以下步骤：
(1)根据电气耦合系统的联络线和联络管道，将一个大型的电气耦合系统分为S个区域，区域编号为s＝1,2,…,S，不同的区域之间通过联络管道传输天然气，通过联络线传输电能，如图1所示；
(2)从S个区域中任取两个相连的区域，假设两个相连区域的编号分别为s1和s2，定义
Figure BDA0002344833650000095
为两个相连区域的联络线与s1相连节点的相角在区域s1的取值，
Figure BDA0002344833650000097
为两个相连区域的联络线与s1相连节点的相角在区域s1取值，
Figure BDA0002344833650000096
为两个相连区域的联络管道中天然气流量在区域s1的取值，
Figure BDA00023448336500001010
为两个相连区域的联络线与s2相连节点的相角在区域s2取值，
Figure BDA00023448336500001011
为两个相连区域的联络线与s1相连节点的相角在区域s2的取值，
Figure BDA00023448336500001012
为两个相连区域的联络管道中的天然气流量在区域s2的取值，设定相连两个区域上同一条联络线上的耦合变量满足：
Figure BDA0002344833650000101
Figure BDA0002344833650000102
同一条联络管道上的耦合变量满足：
Figure BDA0002344833650000103
(3)建立每个区域s的内部优化模型，内部优化模型的约束条件包括：
(3-1)区域s的电网约束包括：
(3-1-1)电网节点电量平衡约束：
Figure BDA0002344833650000104
其中，k为电网中的节点编号，Ne,s为区域s内的电网节点总数；上标cfu表示燃煤电站，Pi cfu表示燃煤电站i的发电有功功率，为待求变量，
Figure BDA0002344833650000105
表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和；上标gfu表示燃气电站，
Figure BDA00023448336500001013
表示燃气电站j的发电有功功率，为待求变量，
Figure BDA0002344833650000106
表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和；上标ld表示电负荷，Pl ld表示负荷l的有功功率，Lsl表示负荷l的弃电有功功率，为待求变量，
Figure BDA0002344833650000107
表示区域s中所有与电网节点k相连的负荷实际有功功率之和；pfmk表示区域s中节点m和节点k所连接支路的有功功率，为待求变量，规定从节点m流向节点k为正，反向为负；
Figure BDA0002344833650000108
表示所有流入节点k的联络线支路的有功功率；
(3-1-2)电网直流潮流约束：
Figure BDA0002344833650000109
其中，θm、θk分别表示区域s中节点m和节点k的电压相角，xmk表示节点m和节点k间所连支路的电抗；
(3-1-3)电网中参考节点的电压相角约束：
θh＝0,h∈refp
其中，θh表示电网节点h的相角，refp表示电网选取的参考节点集合；
(3-1-4)电网支路的有功功率上限和下限约束：
Figure BDA0002344833650000111
其中，
Figure BDA0002344833650000115
表示区域s中节点m和节点k间所连支路的有功功率上限，由电网调度给定；
(3-1-5)电网燃煤电站的有功功率上限和下限约束：
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中，Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限，由电网调度给定；
(3-1-6)电网燃气电站有功功率上限和下限约束：
Figure BDA0002344833650000112
其中，
Figure BDA0002344833650000116
分别表示燃气电站j的发电有功功率下限和上限，由电网调度给定；
(3-1-7)电网弃负荷功率上限和下限约束：
Figure BDA0002344833650000113
其中，
Figure BDA0002344833650000117
分别表示区域s中电负荷l的允许弃负荷功率上限，由电网调度给定；
(3-2)区域s的天然气网约束，包括：
(3-2-1)天然气网压力参考节点约束：
ωq＝const,q＝refg
其中，ωq表示天然气网节点q的压力，refg表示天然气网选取的参考节点集合；
(3-2-2)天然气网节点的天然气流量平衡约束：
Figure BDA0002344833650000114
其中，w为天然气网中的节点编号，Ng,s为区域s中的天然气网节点总数；上标well表示天然气井，
Figure BDA0002344833650000129
表示天然气井d的出气流量，为待求变量，
Figure BDA0002344833650000121
表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和；上标res表示居民负荷，
Figure BDA00023448336500001210
表示居民负荷e消耗的天然气流量，由天然气网调度给定，
Figure BDA0002344833650000122
表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和；gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量，为待求变量，规定从天然气网节点n流向节点w为正，反向为负；
Figure BDA0002344833650000123
表示所有流入节点w的联络管道中的流量；
(3-2-3)天然气网节点压力的上限和下限约束：
Figure BDA0002344833650000124
其中，
Figure BDA00023448336500001211
分别表示天然气网节点q的压力下限和上限，由天然气网调度给定；
(3-2-4)天然气网中天然气井出气流量的上限和下限约束：
Figure BDA0002344833650000125
其中，
Figure BDA00023448336500001212
分别表示区域s中天然气井d的出气流量的下限和上限，由天然气网调度给定；
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束：
Figure BDA0002344833650000126
其中，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)的取值为一个二元变量，引入整数变量
Figure BDA00023448336500001213
将该约束转化为以下表达式：
Figure BDA0002344833650000127
Figure BDA0002344833650000128
Figure BDA0002344833650000131
Figure BDA0002344833650000132
Figure BDA0002344833650000133
Figure BDA0002344833650000134
其中，将上述约束
Figure BDA0002344833650000135
进一步松弛为：
Figure BDA0002344833650000136
(3-2-6)天然气网中天然气管道流量的上限约束：
Figure BDA0002344833650000137
其中，
Figure BDA00023448336500001311
表示从节点n和节点w的连接管道的流量上限，由天然气网调度给定；
(3-3)区域s的电网和天然气网之间的耦合约束如下：
Figure BDA0002344833650000138
其中，
Figure BDA00023448336500001312
和
Figure BDA00023448336500001313
为区域s中的燃气电站s的天然气发电量和耗气量二次关系的二次项系数、一次项系数、常数项，由燃气电站给出；
(3-4)区域s独立优化运行时，区域s的联络线和联络管道上的约束为：
Figure BDA0002344833650000139
gflink,s＝0
(3-5)区域s的内部优化模型的目标函数为：
每个区域进行独立优化运行时，目标函数表示为：
Figure BDA00023448336500001310
其中，OCs表示区域s的运行成本，
Figure BDA00023448336500001314
表示区域s的电负荷l的弃负荷单位成本，
Figure BDA00023448336500001315
表示天然气网的天然气井d的产气单位成本，η为罚函数因子，取值范围为1-100，本发明一个实施例中取值10，
Figure BDA00023448336500001316
表示电网中燃煤电站的i的发电单位成本，发电单位成本为发电量的二次函数，二次函数中的各项系数
Figure BDA00023448336500001317
由燃煤电站给定：
Figure BDA0002344833650000141
(3-5)利用分支定界法和内点法，求解上述约束条件和目标函数组成的区域s的内部优化模型，得到各区域s独立优化时的优化成本OCs；
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型，包括：
(4-1)基于讨价还价的博弈模型，分布式两阶段合作运行优化模型的目标函数表示为：
Figure BDA0002344833650000142
其中，TPs为电气耦合系统中区域s的转移支付，即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配，OCs为各区域s单独运行的最小成本，为步骤(3-6)中所求解得到的OCs，TOCs表示区域电气耦合系统s的运行成本，表示为：
Figure BDA0002344833650000143
(4-2)分布式两阶段合作运行优化模型的约束条件为：
(4-2-1)多区域电气耦合系统各区域s，满足步骤(3-1)、(3-2)、(3-3)中的条件约束，以及与区域s相连的联络线的约束：
Figure BDA0002344833650000144
其中，θf,s为联络线的一个节点与s区域相连节点的相角在区域s的取值，
Figure BDA0002344833650000146
为该联络线另一节点的相角在区域s的取值；
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量，满足步骤(2)中的条件约束；
(4-3)多区域电气耦合系统中的所有区域的转移支付满足：
Figure BDA0002344833650000145
(5)采用带协调层的两阶段交替方向乘子法(ADMM)，求解上述步骤(4)的优化模型，包括以下步骤：
(5-1)建立一个信号处理协调层，使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯；
(5-2)多区域电气耦合系统中的各区域s分别同时向信号处理协调层发送一个协调向量
Figure BDA0002344833650000153
协调向量
Figure BDA0002344833650000154
中包括所有相应区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量；
(5-3)对协调层进行初始化，设初始化迭代次数k＝0，初始化协调层发送给各区域s一个协调变量
Figure BDA0002344833650000155
向量
Figure BDA0002344833650000156
的维数与步骤(5-2)的
Figure BDA0002344833650000157
相同，初始化各区域s的拉格朗日乘子
Figure BDA0002344833650000158
向量
Figure BDA0002344833650000159
的维数与
Figure BDA00023448336500001510
相同，惩罚因子
Figure BDA00023448336500001511
向量
Figure BDA00023448336500001512
的维数与
Figure BDA00023448336500001513
相同，向量中元素取值范围为0.1-10，本发明一个实施例中所有元素取值1，设定收敛阈值
Figure BDA00023448336500001514
取值范围为0.001-0.1，本发明一个实施例中取值0.01，并将初始化的
Figure BDA00023448336500001515
传递给各区域s；
(5-4)各区域s收到协调层初始化的
Figure BDA00023448336500001516
后，建立本区域的第一阶段优化模型，该第一阶段优化模型的约束为步骤(4-2-1)的条件约束，第一阶段优化模型目标函数为：
Figure BDA0002344833650000151
其中，
Figure BDA00023448336500001517
分别表示
Figure BDA00023448336500001518
的第v个元素；
(5-5)采用分支定界法和内点法，各区域s求解步骤(5-4)建立的第一阶段优化模型，得到各区域s合作运行的运行成本TOCs、协调变量
Figure BDA00023448336500001519
以及各约束中的待求变量，并将
Figure BDA00023448336500001520
发送回协调层；
(5-6)协调层收到所有区域s的协调变量
Figure BDA00023448336500001521
后，建立协调层第一阶段优化模型，协调层第一阶段优化模型的目标函数为：
Figure BDA0002344833650000152
不同区域s的协调变量
Figure BDA00023448336500001522
包括所有区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量，相邻区域联络管道和联络线上的相角和流量满足步骤(2)中的条件约束；
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型，>得到协调变量
Figure BDA00023448336500001523
将该协调变量
Figure BDA0002344833650000166
判断与收敛阈值
Figure BDA0002344833650000167
进行比较，若
Figure BDA0002344833650000168
则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本，相应各约束中的待求变量的取值即为设备运行情况，进行步骤(5-8)，若
Figure BDA0002344833650000169
则令
Figure BDA0002344833650000161
协调层将
Figure BDA00023448336500001610
发送给区域s，令k＝k+1，返回步骤(5-4)；
(5-8)对协调层进行第二次初始化，初始化迭代次数k＝0，初始化所有区域s的拉格朗日乘子
Figure BDA00023448336500001611
ρs＝0，取值范围为0.1-10，本发明一个实施例中取值1，收敛阈值εs，取值范围为0.001-0.1，本发明一个实施例中取值0.01，初始化区域s的转移支付在协调层的取值
Figure BDA00023448336500001612
并将耦合变量
Figure BDA00023448336500001613
拉格朗日乘子
Figure BDA00023448336500001614
发送给各区域s；
(5-9)各区域s收到协调层的耦合变量
Figure BDA00023448336500001615
后，建立本区域的第二阶段优化模型，该第二优化模型无条件约束，第二优化模型的目标函数为：
Figure BDA0002344833650000162
其中，TOCs为各区域s的运行成本，由步骤(5-5)得到，
Figure BDA00023448336500001616
表示区域s的转移支付在区域s的取值，为待求变量。各区域s求解本区域的第二阶段优化模型，得到转移支付
Figure BDA00023448336500001617
并将
Figure BDA00023448336500001618
发送回协调层；
(5-10)协调层收到
Figure BDA00023448336500001619
后，建立协调层的第二阶段优化模型，第二阶段目标函数为：
Figure BDA0002344833650000163
约束为：
Figure BDA0002344833650000164
(5-11)协调层采用内点法，求解步骤(5-10)建立的优化模型，得到
Figure BDA00023448336500001620
判断下面不等式是否成立：
Figure BDA0002344833650000165
若不等式成立，则迭代终止，
Figure BDA0002344833650000172
为最后各系统的转移支付，记为TPs，
若不等式不成立，则令
Figure BDA0002344833650000171
协调层将λi k+1发送给各区域s，令k＝k+1，返回步骤(5-9)；
(6)根据步骤(5-5)得到的各区域s约束中的待求变量取值，即区域s中各燃煤电站i的发电有功功率Pi cfu，各燃气电站j的发电有功功率
Figure BDA0002344833650000174
各电网负荷l的弃电有功功率Lsl，各天然气井d的出气流量
Figure BDA0002344833650000173
各天然气网节点q的压力ωq，作为各区域s中分布式合作优化运行的参数，步骤(5-11)的得到的TPs为各区域s的转移支付，实现多区域电气耦合系统的分布式两阶段合作运行。

Claims (1)

1.一种多区域电气耦合系统的分布式两阶段合作运行方法，其特征在于该方法包括以下步骤：
(1)根据电气耦合系统的联络线和联络管道，将一个大型的电气耦合系统分为S个区域，区域编号为s＝1,2,…,S，不同的区域之间通过联络管道传输天然气，通过联络线传输电能；
(2)从S个区域中任取两个相连的区域，假设两个相连区域的编号分别为s1和s2，定义
Figure FDA0003905959200000011
为两个相连区域的联络线与s1相连节点的相角在区域s1的取值，
Figure FDA0003905959200000012
为两个相连区域的联络线与s1相连节点的相角在区域s1取值，
Figure FDA0003905959200000013
为两个相连区域的联络管道中天然气流量在区域s1的取值，
Figure FDA0003905959200000014
为两个相连区域的联络线与s2相连节点的相角在区域s2取值，
Figure FDA0003905959200000015
为两个相连区域的联络线与s1相连节点的相角在区域s2的取值，
Figure FDA0003905959200000016
为两个相连区域的联络管道中的天然气流量在区域s2的取值，设定相连两个区域上同一条联络线上的耦合变量满足：
Figure FDA0003905959200000017
Figure FDA0003905959200000018
同一条联络管道上的耦合变量满足：
Figure FDA0003905959200000019
(3)建立每个区域s的内部优化模型，内部优化模型的约束条件包括：
(3-1)区域s的电网约束包括：
(3-1-1)电网节点电量平衡约束：
Figure FDA00039059592000000110
其中，k为电网中的节点编号，Ne,s为区域s内的电网节点总数；上标cfu表示燃煤电站，Pi cfu表示燃煤电站i的发电有功功率，为待求变量，
Figure FDA00039059592000000111
表示区域s中所有与节点k相连的燃煤电站的发电有功功率之和；上标gfu表示燃气电站，
Figure FDA00039059592000000112
表示燃气电站j的发电有功功率，为待求变量，
Figure FDA00039059592000000113
表示区域s中所有与电网节点k相连的燃气电站的发电有功功率之和；上标ld表示电负荷，Pl ld表示负荷l的有功功率，Lsl表示负荷l的弃电有功功率，为待求变量，
Figure FDA00039059592000000114
表示区域s中所有与电网节点k相连的负荷实际有功功率之和；pfmk表示区域s中节点m和节点k所连接支路的有功功率，为待求变量，规定从节点m流向节点k为正，反向为负；
Figure FDA0003905959200000021
表示所有流入节点k的联络线支路的有功功率；
(3-1-2)电网直流潮流约束：
Figure FDA0003905959200000022
其中，θm、θk分别表示区域s中节点m和节点k的电压相角，xmk表示节点m和节点k间所连支路的电抗；
(3-1-3)电网中参考节点的电压相角约束：
θh＝0,h∈refp
其中，θh表示电网节点h的相角，refp表示电网选取的参考节点集合；
(3-1-4)电网支路的有功功率上限和下限约束：
Figure FDA0003905959200000023
其中，
Figure FDA0003905959200000024
表示区域s中节点m和节点k间所连支路的有功功率上限，由电网调度给定；
(3-1-5)电网燃煤电站的有功功率上限和下限约束：
Pi cfu,min≤Pi cfu≤Pi cfu,max
其中，Pi cfu,min、Pi cfu,max分别表示燃煤电站i的发电有功功率下限和上限，由电网调度给定；
(3-1-6)电网燃气电站有功功率上限和下限约束：
Figure FDA0003905959200000025
其中，
Figure FDA0003905959200000026
分别表示燃气电站j的发电有功功率下限和上限，由电网调度给定；
(3-1-7)电网弃负荷功率上限和下限约束：
Figure FDA0003905959200000027
其中，
Figure FDA0003905959200000028
分别表示区域s中电负荷l的允许弃负荷功率上限，由电网调度给定；
(3-2)区域s的天然气网约束，包括：
(3-2-1)天然气网压力参考节点约束：
ωq＝const,q＝refg
其中，ωq表示天然气网节点q的压力，refg表示天然气网选取的参考节点集合；
(3-2-2)天然气网节点的天然气流量平衡约束：
Figure FDA0003905959200000031
其中，w为天然气网中的节点编号，Ng,s为区域s中的天然气网节点总数；上标well表示天然气井，
Figure FDA0003905959200000032
表示天然气井d的出气流量，为待求变量，
Figure FDA0003905959200000033
表示区域s中的所有与天然气网中节点w相连的天然气井的出气流量之和；上标res表示居民负荷，
Figure FDA0003905959200000034
表示居民负荷e消耗的天然气流量，由天然气网调度给定，
Figure FDA0003905959200000035
表示区域s中所有与节点w相连的居民负荷消耗的天然气流量之和；gfnw表示从天然气网节点n与天然气网节点w相连管道的天然气流量，为待求变量，规定从天然气网节点n流向节点w为正，反向为负；
Figure FDA0003905959200000036
表示所有流入节点w的联络管道中的流量；
(3-2-3)天然气网节点压力的上限和下限约束：
Figure FDA0003905959200000037
其中，
Figure FDA0003905959200000038
分别表示天然气网节点q的压力下限和上限，由天然气网调度给定；
(3-2-4)天然气网中天然气井出气流量的上限和下限约束：
Figure FDA0003905959200000039
其中，
Figure FDA00039059592000000310
分别表示区域s中天然气井d的出气流量的下限和上限，由天然气网调度给定；
(3-2-5)天然气网中天然气管道的流量和压力之间关系约束：
Figure FDA00039059592000000311
其中，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)的取值为一个二元变量，引入整数变量
Figure FDA00039059592000000312
将该约束转化为以下表达式：
Figure FDA00039059592000000313
Figure FDA00039059592000000314
Figure FDA0003905959200000041
Figure FDA0003905959200000042
Figure FDA0003905959200000043
Figure FDA0003905959200000044
其中，将上述约束
Figure FDA0003905959200000045
进一步松弛为：
Figure FDA0003905959200000046
(3-2-6)天然气网中天然气管道流量的上限约束：
Figure FDA0003905959200000047
其中，
Figure FDA0003905959200000048
表示从节点n和节点w的连接管道的流量上限，由天然气网调度给定；
(3-3)区域s的电网和天然气网之间的耦合约束如下：
Figure FDA0003905959200000049
其中，
Figure FDA00039059592000000410
和
Figure FDA00039059592000000411
为区域s中的燃气电站s的天然气发电量和耗气量二次关系的二次项系数、一次项系数、常数项，由燃气电站给出；
(3-4)区域s独立优化运行时，区域s的联络线和联络管道上的约束为：
Figure FDA00039059592000000412
gflink,s＝0
(3-5)区域s的内部优化模型的目标函数为：
每个区域进行独立优化运行时，目标函数表示为：
Figure FDA00039059592000000413
其中，OCs表示区域s的运行成本，
Figure FDA00039059592000000414
表示区域s的电负荷l的弃负荷单位成本，
Figure FDA00039059592000000415
表示天然气网的天然气井d的产气单位成本，η为罚函数因子，取值范围为1-100，
Figure FDA00039059592000000416
表示电网中燃煤电站i的发电单位成本，发电单位成本为发电量的二次函数，二次函数中的各项系数
Figure FDA00039059592000000417
由燃煤电站给定：
Figure FDA00039059592000000418
(3-6)利用分支定界法和内点法，求解上述约束条件和目标函数组成的区域s的内部优化模型，得到各区域s独立优化时的优化成本OCs；
(4)建立基于讨价还价博弈模型的多区域电气耦合系统的分布式两阶段合作运行优化模型，包括：
(4-1)基于讨价还价的博弈模型，分布式两阶段合作运行优化模型的目标函数表示为：
Figure FDA0003905959200000051
其中，TPs为电气耦合系统中区域s的转移支付，即进行合作运行后多区域电气耦合系统总降低的运行成本对各区域的再分配，OCs为各区域s单独运行的最小成本，为步骤(3-6)中所求解得到的OCs，TOCs表示区域电气耦合系统s的运行成本，表示为：
Figure FDA0003905959200000052
(4-2)分布式两阶段合作运行优化模型的约束条件为：
(4-2-1)多区域电气耦合系统各区域s，满足步骤(3-1)、(3-2)、(3-3)中的条件约束，以及与区域s相连的联络线的约束：
Figure FDA0003905959200000053
其中，θf,s为联络线的一个节点与s区域相连节点的相角在区域s的取值，θt,s为该联络线另一节点的相角在区域s的取值；
(4-2-2)多区域电气耦合系统中任意两个相连区域s1、s2之间联络线两端节点的相角和联络管道上的流量，满足步骤(2)中的条件约束；
(4-3)多区域电气耦合系统中的所有区域的转移支付满足：
Figure FDA0003905959200000054
(5)采用带协调层的两阶段交替方向乘子法(ADMM)，求解上述步骤(4)的优化模型，包括以下步骤：
(5-1)建立一个信号处理协调层，使该信号处理协调层与多区域电气耦合系统中的各区域之间进行独立通讯；
(5-2)多区域电气耦合系统中的各区域s分别同时向信号处理协调层发送一个协调向量
Figure FDA0003905959200000055
协调向量
Figure FDA0003905959200000056
中包括所有相应区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量；
(5-3)对协调层进行初始化，设初始化迭代次数k＝0，初始化协调层发送给各区域s一个协调变量
Figure FDA0003905959200000057
向量
Figure FDA0003905959200000058
的维数与步骤(5-2)的
Figure FDA0003905959200000059
相同，初始化各区域s的拉格朗日乘子
Figure FDA00039059592000000510
向量
Figure FDA00039059592000000511
的维数与
Figure FDA00039059592000000512
相同，惩罚因子
Figure FDA00039059592000000513
向量
Figure FDA00039059592000000514
的维数与
Figure FDA00039059592000000515
相同，向量中元素取值范围为0.1-10，设定收敛阈值
Figure FDA0003905959200000061
取值范围为0.001-0.1，并将初始化的
Figure FDA0003905959200000062
Figure FDA0003905959200000063
传递给各区域s；
(5-4)各区域s收到协调层初始化的
Figure FDA0003905959200000064
后，建立本区域的第一阶段优化模型，该第一阶段优化模型的约束为步骤(4-2-1)的条件约束，第一阶段优化模型目标函数为：
Figure FDA0003905959200000065
其中，
Figure FDA0003905959200000066
分别表示
Figure FDA0003905959200000067
的第v个元素；
(5-5)采用分支定界法和内点法，各区域s求解步骤(5-4)建立的第一阶段优化模型，得到各区域s合作运行的运行成本TOCs、协调变量
Figure FDA0003905959200000068
以及各约束中的待求变量，并将
Figure FDA0003905959200000069
发送回协调层；
(5-6)协调层收到所有区域s的协调变量
Figure FDA00039059592000000610
后，建立协调层第一阶段优化模型，协调层第一阶段优化模型的目标函数为：
Figure FDA00039059592000000611
不同区域s的协调变量
Figure FDA00039059592000000612
包括所有区域s相连的多区域电气耦合系统中联络管道和联络线上的相角和流量，相邻区域联络管道和联络线上的相角和流量满足步骤(2)中的条件约束；
(5-7)协调层求解步骤(5-6)建立的协调层第一阶段优化模型，>得到协调变量
Figure FDA00039059592000000613
将该协调变量
Figure FDA00039059592000000614
判断与收敛阈值
Figure FDA00039059592000000615
进行比较，若
Figure FDA00039059592000000616
则将步骤(5-5)的各区域s合作运行的运行成本TOCs作为区域s的运行成本，相应各约束中的待求变量的取值即为设备运行情况，进行步骤(5-8)，若
Figure FDA00039059592000000617
则令
Figure FDA00039059592000000618
协调层将
Figure FDA00039059592000000619
发送给区域s，令k＝k+1，返回步骤(5-4)；
(5-8)对协调层进行第二次初始化，初始化迭代次数k＝0，初始化所有区域s的拉格朗日乘子
Figure FDA00039059592000000620
ρs＝0，取值范围为0.1-10，收敛阈值εs，取值范围为0.001-0.1，初始化区域s的转移支付在协调层的取值
Figure FDA00039059592000000621
并将耦合变量
Figure FDA00039059592000000622
拉格朗日乘子
Figure FDA00039059592000000623
发送给各区域s；
(5-9)各区域s收到协调层的耦合变量
Figure FDA0003905959200000071
后，建立本区域的第二阶段优化模型，该第二阶段优化模型无条件约束，第二阶段优化模型的目标函数为：
Figure FDA0003905959200000072
其中，TOCs为各区域s的运行成本，由步骤(5-5)得到，
Figure FDA0003905959200000073
表示区域s的转移支付在区域s的取值，为待求变量； 各区域s求解本区域的第二阶段优化模型，得到转移支付
Figure FDA0003905959200000074
并将
Figure FDA0003905959200000075
发送回协调层；
(5-10)协调层收到
Figure FDA0003905959200000076
后，建立协调层的第二阶段优化模型，第二阶段目标函数为：
Figure FDA0003905959200000077
约束为：
Figure FDA0003905959200000078
(5-11)协调层采用内点法，求解步骤(5-10)建立的优化模型，得到
Figure FDA0003905959200000079
判断下面不等式是否成立：
Figure FDA00039059592000000710
若不等式成立，则迭代终止，
Figure FDA00039059592000000711
为最后各系统的转移支付，记为TPs，
若不等式不成立，则令
Figure FDA00039059592000000712
协调层将
Figure FDA00039059592000000713
发送给各区域s，令k＝k+1，返回步骤(5-9)；
(6)根据步骤(5-5)得到的各区域s约束中的待求变量取值，即区域s中各燃煤电站i的发电有功功率Pi cfu，各燃气电站j的发电有功功率
Figure FDA00039059592000000714
各电网负荷l的弃电有功功率Lsl，各天然气井d的出气流量
Figure FDA00039059592000000715
各天然气网节点q的压力ωq，作为各区域s中分布式合作优化运行的参数，步骤(5-11)的得到的TPs为各区域s的转移支付，实现多区域电气耦合系统的分布式两阶段合作运行。
CN201911390588.5A 2019-12-30 2019-12-30 一种多区域电气耦合系统的分布式两阶段合作运行方法 Active CN111049135B (zh)

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 (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113642165B (zh) * 2021-07-29 2024-04-12 西安理工大学 考虑需求侧不确定性的电-气耦合配网的分布式优化方法

Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN108173282B (zh) 一种考虑电转气运行成本综合能源系统优化调度方法
CN110348602B (zh) 计及天然气管网和热力管网特性的综合能源系统优化方法
CN111582658B (zh) 考虑多区域互联协同的能源站设备配置及管线规划方法
CN112668188B (zh) 多园区综合能源系统的分布式鲁棒协同优化调度方法
CN106992555B (zh) 一种辅助服务市场环境下的热电联合系统分布式调度方法
Voropai et al. Modeling of integrated energy supply systems: main principles, model, and applications
CN109389248B (zh) 基于区域热网的综合能源系统热电功率协调调度方法
CN112861357B (zh) 一种多微能网耦合系统分布式协同调度方法及系统
CN111049135B (zh) 一种多区域电气耦合系统的分布式两阶段合作运行方法
CN112035984A (zh) 一种电-气-储区域综合能源系统协同规划方法
US20230318300A1 (en) Feedforward Dynamic Multifunctional Structure with Distributed Energy Storage System
CN112200695A (zh) 一种城市级综合能源系统优化调度方法
CN111931977A (zh) 一种考虑电—热能传输价值投入的虚拟电厂扩展规划模型构建方法
CN112531788A (zh) 考虑多重不确定性和自趋优运行的透明微网群规划方法
CN111724026A (zh) 一种多能源网络与配水网耦合运行的优化方法
CN115082235B (zh) 一种虚拟能源站内部多能共享的交易方法及系统
CN116542447A (zh) 一种电热系统优化调度方法
Yu et al. Energy hub model for energy management in energy internet
CN114118803A (zh) 一种多区域综合能源系统集成与协同优化方法及装置
CN110728441B (zh) 基于序列线性规划的电-气联合市场集中式出清方法
CN111062598B (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