CN111682520B - 用网络流对电-气互联系统优化分析的方法 - Google Patents
用网络流对电-气互联系统优化分析的方法 Download PDFInfo
- Publication number
- CN111682520B CN111682520B CN202010168743.5A CN202010168743A CN111682520B CN 111682520 B CN111682520 B CN 111682520B CN 202010168743 A CN202010168743 A CN 202010168743A CN 111682520 B CN111682520 B CN 111682520B
- Authority
- CN
- China
- Prior art keywords
- node
- gas
- electric
- power
- 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 62
- 238000004458 analytical method Methods 0.000 title claims abstract description 22
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 150
- 239000007789 gas Substances 0.000 claims abstract description 138
- 239000003345 natural gas Substances 0.000 claims abstract description 75
- 230000005540 biological transmission Effects 0.000 claims description 49
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 230000008878 coupling Effects 0.000 claims description 12
- 238000010168 coupling process Methods 0.000 claims description 12
- 238000005859 coupling reaction Methods 0.000 claims description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 239000003245 coal Substances 0.000 claims description 6
- 239000004020 conductor Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000000265 homogenisation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 230000001668 ameliorated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013433 optimization analysis Methods 0.000 description 1
- 239000004575 stone Substances 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
- 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)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种用网络流对电‑气互联系统优化分析的方法,所述方法的步骤包括:确定电‑气互联系统的区域参数和运行模式;建立网络流模型,所述网络流模型包括设置电力系统节点、天然气系统节点、源点和汇点的个数,基于各类节点计算包含网络流各支路的容量约束表和权值表;使用最短路径算法求解源点s到汇点t的最短路径,更新权值表和容量约束表直至所有电负荷和气负荷达到规定数值;最后根据不同的运行模型进行多能分析并求得线路功率情况和费用。本发明能对电‑气互联系统进行同质化分析,根据工程要求可将其展现为能量系统、电力系统、天然气系统,有利于对综合能源系统的求解、分析、优化和展示,具有广泛的应用场景和实用意义。
Description
技术领域
本发明涉及电力系统技术领域,尤其涉及一种用网络流对电—气互联系统优化分析的方法。
背景技术
进入21世纪以来,全球能源生产消费迅猛增长,传统化石能源大规模开发利用,资源枯竭、环境污染等问题日趋严重,研究如何实现对能源的清洁高效已迫在眉睫。能源危机的加重及能源供需多样性的发展引发了对于单能源系统弊端的思考:单能源系统束缚了需求的可变性与系统的灵活性,限制了能量的自由流动,不仅损害了能量利用效率,而且降低了系统供能经济性。可见,唯有将能源系统与单一的能源形式解绑,才能改善能源危机的窘境。相较于单能源系统,多能系统能够提高系统弹性降低特定能源的依赖性与供能风险;具备更多的调度方案,能够灵活地保持更高的能量利用效率,从而为实现区域级供能提供了一种重要解决方案。
电力系统和天然气系统都是综合能源系统中十分重要的部分。对于电—气互联系统的多能流建模及优化分析是综合能源系统研究的基石部分。目前主流方法是通过对模型进行潮流求解辅以优化算法寻找系统最优状态来寻找最优解。分析的难点主要在于寻优的速度与正确性以及模型求解的难易程度。而对于模型求解传统方法主要采用牛顿—拉夫逊法,其中又分为统一求解法和顺序求解法。统一求解法需要对电力系统、天然气系统二者整体建立雅可比矩阵和平衡方程,计算量巨大且出错率高。目前常用的是顺序求解法,对电力系统和天然气系统二者分别求解。但是通过此法,使得综合能源系统电和气之间的关系完全解列,显得泾渭分明。而对于我们一个理想中的电—气互联系统,它应该是高度耦合,配合紧密的。电—气互联系统内部的传输运行本质恰恰又是能源的传输与转化,所以可以将其通过能量这一共同核心要素串联起来,将其建立成一个统一的能量流模型来表示并求解。
发明内容
本发明针对上述问题,提供一种用网络流对电—气互联系统优化分析的方法,该方法可将电—气互联系统同质为能量系统或电力系统或天然气系统,然后通过网络流方法对其统一规划求解,实现安全性分析,并在此基础上达到优化运行效果;
为了达到上述目的,本发明提供了一种用网络流对电—气互联系统优化分析的方法,包括以下步骤:
S1:确定电—气互联系统的区域电力系统参数、区域天然气系统参数以及区域电力系统和天然气系统互联的相关参数;
S2:确定电—气互联系统的运行模式,其中,所述运行模式由两位数确定;所述运行模式决定最后计算结果同质为何种模式系统,其中能量系统为00,为模式一;电力系统为01,为模式二;天然气系统为10,为模式三;
S3:选择合适单位并将能量和功率设为标幺值,将电力系统中的电能流以及天然气系统中气流量换算成标幺值下的能流值;
S4:设所述电—气互联系统的网络流模型中电力系统节点数为ne,天然气系统节点数为ng,在此基础上增加一个源点s和一个汇点t,总节点数为n;确定电—气互联系统的网络流模型中发电机节点和气源节点,将所述发电机节点和气源节点与所述源点s相连;确定电—气互联系统的网络流模型中电力系统负荷节点和天然气系统负荷节点,将所述电力系统负荷节点和天然气系统负荷节点与所述汇点t相连;
S5:根据电力系统中发电机节点出力上下限、线路功率传输上下限、电负荷节点以及天然气系统中气源节点产气上下限、管道输气上下限、气负荷节点,形成一个包含所述网络流模型各支路n×n的容量约束表;
S6:根据发电机煤耗特性参数、气源费用参数、确定网络流模型各支路的权值表,具体求解所述权值表的方法为:
对于不考虑网损情况下的发电机节点:
式中Wsie为源点s到发电机i的权值;f为发电机煤耗特性函数;PG为发电机有功出力功率;
对于考虑网损情况下的发电机节点:
式中Ploss为网损功率;
对于气源节点:
Wsig=ci (3)
式中ci为气源节点i的单位产气费用;
对于输气管道或输电导线:
Wij=f(R,X,PLij,Cij) (4)
式中Wij为节点i到节点j的输气管道或输电导线权值;R为输电导线电阻;X为输电导线电抗;PLij为节点i到节点j的输气管道或输电导线流动功率;Cij为输气管道常数;具体公式根据所述输气管道或输电导线为输电网、配电网或输气网格而不相同,需根据实际网络给定;
对于电负荷和气负荷节点:
Wit=0 (5)
式中,Wit为电负荷或气负荷节点i到汇点t的权值;
S7:设定合适步长,采用最短路径算法求解所述电—气互联系统网络流模型中源点s到汇点t的最短路径,并根据所述最短路径中所有节点和所有支路线路功率流通情况,更新权值表和容量约束表,循环S7,直到所有电负荷节点和气负荷节点达到规定数值;
S8:根据上述电—气互联系统网络流模型以及所得到最短路径的计算结果,并根据工程要求求解并记录当前的线路功率情况和费用多少;
S9:根据S2得到的运行模式两位数,判断当前系统的运行模式,根据不同的模式进行多能分析;当系统的运行模式为模式一时,电—气互联系统表现为能量系统,线路功率即为系统能量流动情况,因此输出结果即可,结束;
S10:当系统的运行模式为模式二时,电—气互联系统表现为电力系统,线路功率可视为电功率流动;对于简单网络,用下式计算电压情况;对于复杂网络可通过牛拉法,PQ分解等方法进行求解;
式中Pij为i节点和j节点相连支路始端有功功率;Qij为i节点和j节点相连支路始端无功功率;Ui为i节点电压;Uj为j节点电压;j为虚数;
S11:若计算得到的电压不越线,则输出所述电压、线路电功率、费用等结果,结束;若计算得到的电压越线,设定电压值为临界值,通过式(6)或在牛拉法、PQ分解法等方法中设为PV节点重新计算,然后输出所述电压、线路电功率、费用等结果,结束;
S12:当系统的运行模式为模式三时,电—气互联系统表现为天然气系统,线路功率可转化为气流量流动;对于简单网络,用下式计算气压情况;对于复杂网络可通过牛拉法等方法进行求解;
式中,fij为i节点和j节点相连支路气流量;pi为i节点气压;pj为j节点气压;Cij为管道常数;
S13:若计算得到的气压不越线,则输出所述气压、线路气流量、费用等结果,结束;若计算得到的气压越线,设定气压值为临界值,通过式(7)或在牛拉法、PQ分解法等方法中设为恒压节点重新计算,然后输出所述气压、线路气流量、费用等结果,结束。
优选方式下,区域电力系统由大电网、热电联产机组、燃气轮机式发电机等供电,区域天然气系统由天然气气源、电转气设备供气,区域综合能源系统通过电、热、气能源网络传输能量以满足用户及耗能设备的能源需求。
优选方式下,所述区域电力系统参数包括但不限于:符合MATPOWER标准的区域电力系统基准功率,母线参数,发电机参数,支路参数;所述区域天然气系统参数包括但不限于:天然气节点参数矩阵,包含节点编号、节点类型、天然气负荷、节点压强;天然气管道参数,包含管道编号、管道首节点、管道末节点、管道直径、管道长度、管道绝对粗糙度、效率因子;天然气参数,包含天然气费用、天然气比重;所述区域电力系统和天然气系统互联的相关参数包括但不限于:电转换效率,燃气轮机式发电机效率,P2G设备效率。
优选方式下,对于电气耦合节点需要考虑效率或转化率等问题,因此根据效率或转化率参数对其流入流出能量进行折算,将其转化成能流,整合入电—气互联系统网络流模型的能流中,之后求解权值表和容量约束表的方法与其余节点相同;
P入=η*P出 (8)
式中P入为该电气耦合节点流入能量;P出为该电气耦合节点流出能量;η为该电气耦合节点流入流出能量的效率或转化率。
优选方式下,所述电—气互联系统网络流模型中容量约束表的具体求解方法为:
对于发电机节点:
式中为源点s到发电机i的容量约束下限;/>为发电机i有功出力下限;/>为源点s到发电机i的容量约束上限;/>发电机i有功出力上限;
对于气源节点:
式中为源点s到气源节点i的容量约束下限;/>为气源节点i有功出力下限;为源点s到气源节点i的容量约束上限;/>为气源节点i有功出力上限;
对于输气管道或输电导线:
式中为节点i到节点j的容量约束下限;/>为输气管道或输电导线传输功率下限;/>为节点i到节点j的容量约束上限;/>为输气管道或输电导线传输功率上限;
对于电负荷或气负荷节点:
式中为节点i到汇点t的容量约束下限;Pd为电负荷或气负荷值;/>为节点i到汇点t的容量约束上限。
本发明的有益效果是:本发明提出的一种对网络流对电—气互联系统优化分析的方法,能对电—气互联系统进行同质化分析,根据工程要求可将其展现为能量系统、电力系统、天然气系统,有利于对综合能源系统的求解、分析、优化和展示,具有广泛的应用场景和实用意义。
附图说明
图1为本发明提供的一种电—气互联系统示意图;
图2为本发明提供的一种电—气互联系统节点图;
图3为本发明提供的一种网络流模型示意图;
图4为本发明提供的一种电—气互联系统容量表和权值表的示意表;
图5为本发明提供的一种电—气互联系统供能节点提供功率优化结果;
图6为本发明提供的一种电—气互联系统各支路功率优化结果。
具体实施方式
本发明介绍一种用网络流对电—气互联系统优化分析的方法,该方法可将电—气互联系统同质为能量系统或电力系统或天然气系统,然后通过网络流方法对其统一规划求解,实现安全性分析,并在此基础上达到优化运行效果;以IEEE14节点电力系统模型和比利时天然气20节点耦合模型作为示例,示意图如图1所示。图1左半部分为MATPOWER上电力系统标准14节点模型,右半部分为比利时天然气20节点模型,中间连接线为电力系统与天然气系统之间的联络线。图2为图一的节点抽象图,转化后的网络流模型图如图3所示,红线代表能量从源点汇入各发电机与气源,为产能弧,蓝线为各电负荷和气负荷汇入汇点,为用能弧,黑线在电力系统中为输电管道,在天然气系统中为输气管道,黄线为电力系统与天然气系统之间的联络线。其中比利时天然气能量规模远大于IEEE14节点,为了提高耦合度,对天然气系统的负荷量折算成实际的十分之一。
本发明一种用网络流对电—气互联系统优化分析的方法,包括以下步骤:
S1:确定电—气互联系统的区域电力系统参数、区域天然气系统参数以及区域电力系统和天然气系统互联的相关参数;
S2:确定电—气互联系统的运行模式,其中,所述运行模式由两位数确定;所述运行模式决定最后计算结果同质为何种模式系统,其中能量系统为00,为模式一;电力系统为01,为模式二;天然气系统为10,为模式三;
S3:选择合适单位并将能量和功率设为标幺值,将电力系统中的电能流以及天然气系统中气流量换算成标幺值下的能流值;
S4:设所述电—气互联系统的网络流模型中电力系统节点数为ne,天然气系统节点数为ng,在此基础上增加一个源点s和一个汇点t,总节点数为n;确定电—气互联系统的网络流模型中发电机节点和气源节点,将所述发电机节点和气源节点与所述源点s相连;确定电—气互联系统的网络流模型中电力系统负荷节点和天然气系统负荷节点,将所述电力系统负荷节点和天然气系统负荷节点与所述汇点t相连;
S5:根据电力系统中发电机节点出力上下限、线路功率传输上下限、电负荷节点以及天然气系统中气源节点产气上下限、管道输气上下限、气负荷节点,形成一个包含所述网络流模型各支路n×n的容量约束表,如图4所示,电—气互联系统容量表和权值表的示意表;
S6:根据发电机煤耗特性参数、气源费用参数、确定网络流模型各支路的权值表,如图4所示,电—气互联系统容量表和权值表的示意表,具体求解所述权值表的方法为:
对于不考虑网损情况下的发电机节点:
式中Wsie为源点s到发电机i的权值;f为发电机煤耗特性函数;PG为发电机有功出力功率;
对于考虑网损情况下的发电机节点:
式中Ploss为网损功率;
对于气源节点:
Wsig=ci (3)
式中ci为气源节点i的单位产气费用;
对于输气管道或输电导线:
Wij=f(R,X,PLij,Cij) (4)
式中Wij为节点i到节点j的输气管道或输电导线权值;R为输电导线电阻;X为输电导线电抗;PLij为节点i到节点j的输气管道或输电导线流动功率;Cij为输气管道常数;具体公式根据所述输气管道或输电导线为输电网、配电网或输气网格而不相同,需根据实际网络给定;
对于电负荷和气负荷节点:
Wit=0 (5)
式中,Wit为电负荷或气负荷节点i到汇点t的权值;
S7:设定合适步长,采用Dijkstra算法,即一种最短路径算法求解所述电—气互联系统网络流模型中源点s到汇点t的最短路径,并根据所述最短路径中所有节点和所有支路线路功率流通情况,更新权值表和容量约束表,循环S7,直到所有电负荷节点和气负荷节点达到规定数值;
S8:根据上述电—气互联系统网络流模型以及所得到最短路径的计算结果,并根据工程要求求解并记录当前的线路功率情况和费用多少;
S9:根据S2得到的运行模式两位数,判断当前系统的运行模式,根据不同的模式进行多能分析;当系统的运行模式为模式一时,电—气互联系统表现为能量系统,线路功率即为系统能量流动情况,因此输出结果即可,结束;
S10:当系统的运行模式为模式二时,电—气互联系统表现为电力系统,线路功率可视为电功率流动;对于简单网络,用下式计算电压情况;对于复杂网络可通过牛拉法,PQ分解等方法进行求解;
式中Pij为i节点和j节点相连支路始端有功功率;Qij为i节点和j节点相连支路始端无功功率;Ui为i节点电压;Uj为j节点电压;j为虚数;
S11:若计算得到的电压不越线,则输出所述电压、线路电功率、费用等结果,结束;若计算得到的电压越线,设定电压值为临界值,通过式(6)或在牛拉法、PQ分解法等方法中设为PV节点重新计算,然后输出所述电压、线路电功率、费用等结果,结束;
S12:当系统的运行模式为模式三时,电—气互联系统表现为天然气系统,线路功率可转化为气流量流动;对于简单网络,用下式计算气压情况;对于复杂网络可通过牛拉法等方法进行求解;
式中,fij为i节点和j节点相连支路气流量;pi为i节点气压;pj为j节点气压;Cij为管道常数;
S13:若计算得到的气压不越线,则输出所述气压、线路气流量、费用等结果,结束;若计算得到的气压越线,设定气压值为临界值,通过式(7)或在牛拉法、PQ分解法等方法中设为恒压节点重新计算,然后输出所述气压、线路气流量、费用等结果,结束。
区域电力系统由大电网、热电联产机组(Combined Heat and Power,CHP)、燃气轮机式发电机等供电,区域天然气系统由天然气气源、电转气(Power to Gas,P2G)设备供气,区域综合能源系统通过电、热、气能源网络传输能量以满足用户及耗能设备的能源需求。
所述区域电力系统参数包括但不限于:符合MATPOWER标准的区域电力系统基准功率,母线参数,发电机参数,支路参数;所述区域天然气系统参数包括但不限于:天然气节点参数矩阵,包含节点编号、节点类型、天然气负荷、节点压强;天然气管道参数,包含管道编号、管道首节点、管道末节点、管道直径、管道长度、管道绝对粗糙度、效率因子;天然气参数,包含天然气费用、天然气比重;所述区域电力系统和天然气系统互联的相关参数包括但不限于:电转换效率,燃气轮机式发电机效率,P2G设备效率。
对于电气耦合节点需要考虑效率或转化率等问题,因此根据效率或转化率参数对其流入流出能量进行折算,将其转化成能流,整合入电—气互联系统网络流模型的能流中,之后求解权值表和容量约束表的方法与其余节点相同;
P入=η*P出 (8)
式中P入为该电气耦合节点流入能量;P出为该电气耦合节点流出能量;η为该电气耦合节点流入流出能量的效率或转化率。
所述电—气互联系统网络流模型中容量约束表的具体求解方法为:
对于发电机节点:
式中为源点s到发电机i的容量约束下限;/>为发电机i有功出力下限;/>为源点s到发电机i的容量约束上限;/>发电机i有功出力上限;
对于气源节点:
式中为源点s到气源节点i的容量约束下限;/>为气源节点i有功出力下限;/>为源点s到气源节点i的容量约束上限;/>为气源节点i有功出力上限;
对于输气管道或输电导线:
式中为节点i到节点j的容量约束下限;/>为输气管道或输电导线传输功率下限;/>为节点i到节点j的容量约束上限;/>为输气管道或输电导线传输功率上限;
对于电负荷或气负荷节点:
式中为节点i到汇点t的容量约束下限;Pd为电负荷或气负荷值;/>为节点i到汇点t的容量约束上限。
图5、图6为该示例下电—气互联系统网络流优化的计算结果,可适当参考。图5为电—气互联系统供能节点提供功率优化结果,其中2、7节点为发电机,16、20、23、28、29节点为气源,所有供能点提供能量全转换为功率形式。图六为电—气互联系统各支路功率优化结果,其中1-20支路为电力系统标准14节点模型中的所有支路,21-44为比利时天然气20节点模型中的所有支路,45-46为电力系统与天然气系统之间的联络线。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (5)
1.一种用网络流对电—气互联系统优化分析的方法,其特征在于,包括以下步骤:
S1:确定电—气互联系统的区域电力系统参数、区域天然气系统参数以及区域电力系统和天然气系统互联的相关参数;
S2:确定电—气互联系统的运行模式,其中,所述运行模式由两位数确定;所述运行模式决定最后计算结果同质为何种模式系统,其中能量系统为00,为模式一;电力系统为01,为模式二;天然气系统为10,为模式三;
S3:选择合适单位并将能量和功率设为标幺值,将电力系统中的电能流以及天然气系统中气流量换算成标幺值下的能流值;
S4:设所述电—气互联系统的网络流模型中电力系统节点数为ne,天然气系统节点数为ng,在此基础上增加一个源点s和一个汇点t,总节点数为n;确定电—气互联系统的网络流模型中发电机节点和气源节点,将所述发电机节点和气源节点与所述源点s相连;确定电—气互联系统的网络流模型中电力系统负荷节点和天然气系统负荷节点,将所述电力系统负荷节点和天然气系统负荷节点与所述汇点t相连;
S5:根据电力系统中发电机节点出力上下限、线路功率传输上下限、电负荷节点以及天然气系统中气源节点产气上下限、管道输气上下限、气负荷节点,形成一个包含所述网络流模型各支路n×n的容量约束表;
S6:根据发电机煤耗特性参数、气源费用参数、确定网络流模型各支路的权值表,具体求解所述权值表的方法为:
对于不考虑网损情况下的发电机节点:
式中Wsie为源点s到发电机i的权值;f为发电机煤耗特性函数;PG为发电机有功出力功率;
对于考虑网损情况下的发电机节点:
式中Ploss为网损功率;
对于气源节点:
Wsig=ci (3)
式中ci为气源节点i的单位产气费用;
对于输气管道或输电导线:
Wij=f(R,X,PLij,Cij) (4)
式中Wij为节点i到节点j的输气管道或输电导线权值;R为输电导线电阻;X为输电导线电抗;PLij为节点i到节点j的输气管道或输电导线流动功率;Cij为输气管道常数;具体公式根据所述输气管道或输电导线为输电网、配电网或输气网格而不相同,需根据实际网络给定;
对于电负荷和气负荷节点:
Wit=0 (5)
式中,Wit为电负荷或气负荷节点i到汇点t的权值;
S7:设定合适步长,采用最短路径算法求解所述电—气互联系统网络流模型中源点s到汇点t的最短路径,并根据所述最短路径中所有节点和所有支路线路功率流通情况,更新权值表和容量约束表,循环S7,直到所有电负荷节点和气负荷节点达到规定数值;
S8:根据上述电—气互联系统网络流模型以及所得到最短路径的计算结果,并根据工程要求求解并记录当前的线路功率情况和费用多少;
S9:根据S2得到的运行模式两位数,判断当前系统的运行模式,根据不同的模式进行多能分析;当系统的运行模式为模式一时,电—气互联系统表现为能量系统,线路功率即为系统能量流动情况,因此输出结果即可,结束;
S10:当系统的运行模式为模式二时,电—气互联系统表现为电力系统,线路功率可视为电功率流动;对于简单网络,用下式计算电压情况;对于复杂网络可通过牛拉法,PQ分解方法进行求解;
式中Pij为i节点和j节点相连支路始端有功功率;Qij为i节点和j节点相连支路始端无功功率;Ui为i节点电压;Uj为j节点电压;j为虚数;
S11:若计算得到的电压不越线,则输出所述电压、线路电功率、费用结果,结束;若计算得到的电压越线,设定电压值为临界值,通过式(6)或在牛拉法、PQ分解法方法中设为PV节点重新计算,然后输出所述电压、线路电功率、费用结果,结束;
S12:当系统的运行模式为模式三时,电—气互联系统表现为天然气系统,线路功率可转化为气流量流动;对于简单网络,用下式计算气压情况;对于复杂网络可通过牛拉法方法进行求解;
式中,fij为i节点和j节点相连支路气流量;pi为i节点气压;pj为j节点气压;Cij为管道常数;
S13:若计算得到的气压不越线,则输出所述气压、线路气流量、费用结果,结束;若计算得到的气压越线,设定气压值为临界值,通过式(7)或在牛拉法、PQ分解法等方法中设为恒压节点重新计算,然后输出所述气压、线路气流量、费用结果,结束。
2.根据权利要求1所述用网络流对电—气互联系统优化分析的方法,其特征在于,区域电力系统由大电网、热电联产机组、燃气轮机式发电机供电,区域天然气系统由天然气气源、电转气设备供气,区域综合能源系统通过电、热、气能源网络传输能量以满足用户及耗能设备的能源需求。
3.根据权利要求1所述用网络流对电—气互联系统优化分析的方法,其特征在于,所述区域电力系统参数包括:符合MATPOWER标准的区域电力系统基准功率,母线参数,发电机参数,支路参数;所述区域天然气系统参数包括:天然气节点参数矩阵,包含节点编号、节点类型、天然气负荷、节点压强;天然气管道参数,包含管道编号、管道首节点、管道末节点、管道直径、管道长度、管道绝对粗糙度、效率因子;天然气参数,包含天然气费用、天然气比重;所述区域电力系统和天然气系统互联的相关参数包括:电转换效率,燃气轮机式发电机效率,P2G设备效率。
4.根据权利要求1所述用网络流对电—气互联系统优化分析的方法,其特征在于,对于电气耦合节点需要考虑效率或转化率问题,因此根据效率或转化率参数对其流入流出能量进行折算,将其转化成能流,整合入电——气互联系统网络流模型的能流中,之后求解权值表和容量约束表的方法与其余节点相同;
P入=η*P出 (8)
式中P入为该电气耦合节点流入能量;P出为该电气耦合节点流出能量;η为该电气耦合节点流入流出能量的效率或转化率。
5.根据权利要求1所述用网络流对电—气互联系统优化分析的方法,其特征在于,所述电—气互联系统网络流模型中容量约束表的具体求解方法为:
对于发电机节点:
式中为源点s到发电机i的容量约束下限;/>为发电机i有功出力下限;/>为源点s到发电机i的容量约束上限;/>发电机i有功出力上限;
对于气源节点:
式中为源点s到气源节点i的容量约束下限;/>为气源节点i有功出力下限;为源点s到气源节点i的容量约束上限;/>为气源节点i有功出力上限;
对于输气管道或输电导线:
式中为节点i到节点j的容量约束下限;/>为输气管道或输电导线传输功率下限;/>为节点i到节点j的容量约束上限;/>为输气管道或输电导线传输功率上限;
对于电负荷或气负荷节点:
式中为节点i到汇点t的容量约束下限;Pd为电负荷或气负荷值;/>为节点i到汇点t的容量约束上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010168743.5A CN111682520B (zh) | 2020-03-11 | 2020-03-11 | 用网络流对电-气互联系统优化分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010168743.5A CN111682520B (zh) | 2020-03-11 | 2020-03-11 | 用网络流对电-气互联系统优化分析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111682520A CN111682520A (zh) | 2020-09-18 |
CN111682520B true CN111682520B (zh) | 2023-08-22 |
Family
ID=72451342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010168743.5A Active CN111682520B (zh) | 2020-03-11 | 2020-03-11 | 用网络流对电-气互联系统优化分析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111682520B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107394828A (zh) * | 2017-08-24 | 2017-11-24 | 河海大学 | 一种基于概率潮流的电气互联综合能源系统优化运行分析方法 |
CN109446545A (zh) * | 2018-09-07 | 2019-03-08 | 天津大学 | 基于通用能量母线的综合能源系统能量流分析计算方法 |
CN109740242A (zh) * | 2018-12-29 | 2019-05-10 | 重庆大学 | 考虑天然气热力过程的电-气综合能源系统统一能流计算方法 |
CN109800968A (zh) * | 2018-12-29 | 2019-05-24 | 重庆大学 | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 |
EP3582359A1 (en) * | 2018-04-20 | 2019-12-18 | Northeastern University | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
-
2020
- 2020-03-11 CN CN202010168743.5A patent/CN111682520B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107394828A (zh) * | 2017-08-24 | 2017-11-24 | 河海大学 | 一种基于概率潮流的电气互联综合能源系统优化运行分析方法 |
EP3582359A1 (en) * | 2018-04-20 | 2019-12-18 | Northeastern University | Stability evaluation and static control method for electricity-heat-gas integrated energy system |
CN109446545A (zh) * | 2018-09-07 | 2019-03-08 | 天津大学 | 基于通用能量母线的综合能源系统能量流分析计算方法 |
CN109740242A (zh) * | 2018-12-29 | 2019-05-10 | 重庆大学 | 考虑天然气热力过程的电-气综合能源系统统一能流计算方法 |
CN109800968A (zh) * | 2018-12-29 | 2019-05-24 | 重庆大学 | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 |
Non-Patent Citations (1)
Title |
---|
区域能源互联网"站-网"布局优化研究;陈娟;黄元生;鲁斌;;中国电机工程学报(第03期);116-25 * |
Also Published As
Publication number | Publication date |
---|---|
CN111682520A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104734155B (zh) | 一种获取电‑气互联能源系统可用输电能力的方法 | |
CN103490410B (zh) | 一种基于多目标优化的微电网规划和容量配置方法 | |
CN105576657B (zh) | 主动配电网与输电网集成系统随机模糊连续潮流建模方法 | |
CN104281737B (zh) | 一种海上风电场集电系统拓扑结构选型方法 | |
CN107808216B (zh) | 电-气-热互联系统弃风弃光和电气热负荷削减综合最小优化模型构建方法 | |
CN106446399B (zh) | 电缆布局方法及用于电缆布局的终端设备 | |
CN110797874B (zh) | 含电力电子变压器的交直流混合配电网状态估计方法 | |
CN114077934A (zh) | 一种综合能源微网互联系统及其调度方法 | |
CN115392035A (zh) | 电-气互联综合能源系统优化调度的改进方法 | |
CN114092277A (zh) | 柔性区域综合能源系统冷热气电联供调度方法和系统 | |
CN107404118A (zh) | 基于随机响应面法的电气互联系统概率最优潮流计算方法 | |
CN114493258A (zh) | 一种基于电气剖分的电网碳足迹追踪方法 | |
CN116388147A (zh) | 一种电-气耦合的综合能源系统碳排放监测方法 | |
CN113131513B (zh) | 考虑碳排放的电、热、气转换优化运行配置方法和存储介质 | |
CN111682520B (zh) | 用网络流对电-气互联系统优化分析的方法 | |
CN106451465B (zh) | 一种用于avc控制中基于cim的新能源模型拼接方法 | |
CN112836957A (zh) | 一种考虑供电可靠性的区域综合能源系统规划方法 | |
CN112214883A (zh) | 一种微能源网综合能源系统运行优化方法和系统 | |
CN112381267A (zh) | 一种面向工业园区的供水-供能耦合规划方法 | |
CN111724026A (zh) | 一种多能源网络与配水网耦合运行的优化方法 | |
CN106786766A (zh) | 一种基于p2g技术的提高风电最大并网容量的方法 | |
CN116231634A (zh) | 一种计及韧性提升的多能耦合配电系统故障恢复方法 | |
CN206611076U (zh) | 一种带温度测控功能的光伏交流汇流箱 | |
CN111275271B (zh) | 一种冷-热-电综合能源系统潮流优化方法 | |
CN114066056A (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 |