CN107947245B - 考虑天然气系统约束的等值最优潮流模型构建方法 - Google Patents

考虑天然气系统约束的等值最优潮流模型构建方法 Download PDF

Info

Publication number
CN107947245B
CN107947245B CN201711172326.2A CN201711172326A CN107947245B CN 107947245 B CN107947245 B CN 107947245B CN 201711172326 A CN201711172326 A CN 201711172326A CN 107947245 B CN107947245 B CN 107947245B
Authority
CN
China
Prior art keywords
gas
natural gas
flow
power
constraint
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.)
Expired - Fee Related
Application number
CN201711172326.2A
Other languages
English (en)
Other versions
CN107947245A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201711172326.2A priority Critical patent/CN107947245B/zh
Publication of CN107947245A publication Critical patent/CN107947245A/zh
Application granted granted Critical
Publication of CN107947245B publication Critical patent/CN107947245B/zh
Expired - Fee Related 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
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Control Of Eletrric Generators (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开一种考虑天然气系统约束的等值最优潮流模型构建方法。在电力系统与天然气系统数据各自独立调度运行，数据不能共享的情况下，提出天然气系统的约束等值模型，将天然气系统等值约束信息提供给电力系统而不涉密；在此基础上建立考虑天然气系统约束的等值最优潮流方法，避免了电力系统与天然气系统的交替迭代，同时还有效计及了天然气系统对电力系统运行调度的影响。最后，通过IEEE14‑NGS10电‑气互联系统的算例分析验证了所提方法的有效性。

Description

考虑天然气系统约束的等值最优潮流模型构建方法
技术领域
本发明属于电-气互联系统领域，目的是实现考虑天然气系统约束的等值最优潮流计算。具体涉及考虑气源容量约束、节点气压约束、压缩机气压变比约束的天然气系统约束的等值模型，以及基于约束等值模型提出考虑天然气系统约束的等值最优潮流方法。
背景技术
相比于燃煤/燃油发电，天然气发电具有污染低、经济性好、动作响应快、安装时间短等优点，近年来，燃气轮机装机容量在世界很多地方都持续增长。随着燃气轮机装机容量的日益增长，天然气系统与电力系统的耦合也愈加紧密。实际上，天然气系统与电力系统是独立运行调度，并没有考虑系统间的相互影响，导致优化结果不准确。例如：电力系统燃气轮机发电量较大时，大容量的天然气消耗可能会迫使天然气系统的运行状态超出安全约束(如：节点压力越限、管道流量越限)，从而影响燃气轮机的一次能源供应。目前，国内外学者对电-气互联系统的协同优化已有许多研究，可分为天然气系统与电力系统的数据完全共享和数据不能共享两类，但存在以下问题：数据完全共享方法对天然气系统与电力系统统一建模优化分析，但未考虑实际上电力系统与天然气系统分别属于不同的运营商负责管理，数据具有商业机密，获得整个电-气互联系统的全局信息目前还难以实现。数据不能共享方法在构建分解协同交互机的基础上，采用交替迭代的思想对电-气互联系统进行优化分析，然而这些方法存在着参数设定困难、收敛速度慢甚至不收敛等问题。
发明内容
本发明的目的提供一种考虑天然气系统约束的等值最优潮流模型构建方法，其特征在于，包括以下步骤:
1)输入基础数据
输入电-气互联系统基础数据，包括电力子系统设备参数(即线路的阻抗、对地电纳及其传输功率约束条件，变压器的阻抗、对地导纳、变比及其传输功率约束条件，电力系统节点的对地导纳及连接负荷功率大小，发电机出力大小及其出力约束条件)、天然气子系统设备参数(即管道的管道常数，压缩机的压变比约束、温度常数和压缩因子常数、流过流量约束条件)和拓扑结构(即包括电力系统、天然气系统中各节点连接关系，以及各耦合元件的连接关系)。
2)建立天然气系统的约束等值模型
为有效保留天然气系统约束对电力系统的影响，本发明提出天然气系统的约束等值模型。约束等值模型的建模可分为以下三个步骤：建立约束等值目标函数、建立约束等值等式约束方程、建立约束等值不等式约束方程。
2-1)建立约束等值目标函数
根据天然气系统中燃气轮机个数(N)建立多个目标函数，首先建立任意一个燃气轮机消耗天然气流量最大的目标函数，一共N个；然后建立任意两个燃气轮机消耗天然气流总流量最大的目标函数，一共C2N个；以此类推，最后建立所有燃气轮机消耗天然气总流量最大的目标函数，一共CNN个；即由Ci，Cij，…，Call组成。
Figure GDA0002603116890000021
式中，在建立任意一个燃气轮机消耗天然气流量最大的目标函数时，Ci表示第i个燃气轮机消耗天然气流量的最大值，
Figure GDA0002603116890000022
表示第i个燃气轮机消耗天然气流量。在建立任意两个燃气轮机消耗天然气总流量最大的目标函数时，Cij为第i个和第j个燃气轮机消耗天然气总流量的最大值，
Figure GDA0002603116890000023
和
Figure GDA0002603116890000024
分别表示第i个和第j个燃气轮机消耗的天然气流量。在建立所有燃气轮机消耗天然气总流量最大的目标函数时，Call为所有燃气轮机消耗天然气总流量的最大值，
Figure GDA0002603116890000025
表示第i个燃气轮机消耗天然气流量。在求解燃气轮机消耗天然气流量最大值时，其他不相关燃气轮机消耗天然气流量为0。
2-2)建立约束等值等式约束方程
天然气系统流量平衡方程为式(2)。
FG,m-FGAS,m-FL,m-Fm＝0,m＝1,2,...,Nm (2)
式中，FG,m、FGAS,m、FL,m、Fm分别为天然气节点m的气源流量、燃气轮机消耗的流量、气负荷和注入流量，Nm为天然气系统节点个数。其中节点注入流量Fm可以通过式(3)求得。
Figure GDA0002603116890000031
式中，天然气节点n表示与天然气节点m相连的节点；Fmn、Fcom,mn、τcom,mn分别为管道流量、流过压缩机的流量和压缩机损耗，sgnc(m,n) 和sgnτ(m,n)分别为压缩机流过流量和消耗流量的方向函数。这些量可以通过式(4)和(5)表示。
Figure GDA0002603116890000032
Figure GDA0002603116890000033
式中，Πm是天然气节点m的节点气压；kmn是与管道内径、长度、效率、压缩因子相关的常数；smn反应管道流量方向；
Figure GDA0002603116890000034
为压缩机消耗的功率；Bmn是压缩机系数；Zmn是压缩机压缩因子相关的常数。ατ、βτ、γτ为功率转化效率常数。
2-3)建立约束等值不等式约束方程
式(6)-(8)分别为天然气系统节点气压、气源容量、压缩机气压变比的上下限。
Figure GDA0002603116890000035
Figure GDA0002603116890000036
Figure GDA0002603116890000037
式中，Πm是天然气节点m的节点气压，FG,m为气源容量，Rr为压缩机气压变比。
Figure GDA0002603116890000038
和
Figure GDA0002603116890000039
分别为天然气节点气压的上限和下限；
Figure GDA00026031168900000310
和
Figure GDA00026031168900000311
分别为气源容量的上限和下限；
Figure GDA00026031168900000312
和
Figure GDA00026031168900000313
分别为压缩机压变比的上限和下限。Ns和Np分别为气源和压缩机的个数。
求解由公式(1)-(8)构成的优化模型可计算得到目标函数中的
Figure GDA00026031168900000314
通过燃气轮机耗量方程(9)可计算出对应的
Figure GDA00026031168900000315
再通过公式(10)计算得到任意一个燃气轮机的出力上限，一共N个；得到任意两个燃气轮机总出力上限，一共C2N个；以此类推，得到所有燃气轮机总出力上限，一共 CNN个；即由
Figure GDA0002603116890000041
组成。
Figure GDA0002603116890000042
Figure GDA0002603116890000043
式(10)中，在求解任意一个燃气轮机的出力上限时，
Figure GDA0002603116890000044
为第i个燃气轮机输出功率的最大值，
Figure GDA0002603116890000045
为第i个燃气轮机的输出功率。在求解任意两个燃气轮机总出力上限时，
Figure GDA0002603116890000046
为第i个和第j 个燃气轮机总输出功率的最大值，
Figure GDA0002603116890000047
和
Figure GDA0002603116890000048
分别表示第i个和第 j个燃气轮机输出功率。在求解所有燃气轮机总出力上限时，
Figure GDA0002603116890000049
为所有燃气轮机总输出功率的最大值，
Figure GDA00026031168900000410
为第i个燃气轮机的输出功率。GHV为天然气高热值。
通过天然气系统的约束等值模型(1)-(10)可将天然气系统气源容量约束、节点气压约束、压缩机气压变比约束等值到耦合元件燃气轮机的有功出力约束中。
3)求解天然气系统的约束等值模型
式(1)-(10)构成求解天然气约束等值参数的非线性最优化问题。通过内点法可直接求得系统中各个燃气轮机组的有功出力约束以替代天然气系统的运行约束。
4)建立考虑天然气系统约束的等值最优潮流模型
通过天然气系统的约束等值模型，天然气系统的约束等值到了燃气轮机的有功出力约束中，因此可以进一步建立考虑天然气系统约束的等值最优潮流模型。等值最优潮流模型的建模可分为以下三个步骤：建立等值最优潮流目标函数、建立等值最优潮流等式约束方程、建立等值最优潮流不等式约束方程。
4-1)等值最优潮流目标函数：
Figure GDA00026031168900000411
式中，fi(PGAS,i)为燃气轮机组的成本函数；fj(PGj)为非燃气发电机的成本函数；NG为非燃气发电机个数。
4-2)建立等值最优潮流等式约束方程：
考虑电力系统功率平衡方程(12)-(13)，构建如下等式约束：
Figure GDA0002603116890000051
Figure GDA0002603116890000052
式中，Pi和Qi分别为节点i的注入有功功率和无功功率；Ui是节点i的电压幅值；θij是节点i与j的电压相角差；Gij和Bij分别为节点导纳矩阵的第i行第j列元素的实部和虚部；NA为节点个数。
4-3)建立等值最优潮流不等式约束方程：
考虑天然气系统的等值约束
Figure GDA0002603116890000053
燃气轮机自身的容量约束、非燃气常规机组容量约束、节点电压约束和线路功率约束构建如下不等式约束。
Figure GDA0002603116890000054
Figure GDA0002603116890000055
Figure GDA0002603116890000056
Figure GDA0002603116890000057
Figure GDA00026031168900000523
Figure GDA0002603116890000058
Figure GDA0002603116890000059
其中，PGAS,i为燃气轮机i的有功功率，QGAS,i为燃气轮机i的无功功率，PGi为非燃气轮机i的有功功率，QGi为非燃气轮机i的无功功率，Ui是节点i的电压幅值，Tl是输电线路l流过的功率；NL为输电线路的总数。
Figure GDA00026031168900000510
和
Figure GDA00026031168900000511
分别为燃气轮机有功出力的上限和下限；
Figure GDA00026031168900000512
和
Figure GDA00026031168900000513
Figure GDA00026031168900000514
分别为燃气轮机无功出力的上限和下限；
Figure GDA00026031168900000515
和
Figure GDA00026031168900000516
分别为节点电压幅值的上限和下限；
Figure GDA00026031168900000517
和
Figure GDA00026031168900000518
为输电线路传输功率的上限和下限；
Figure GDA00026031168900000519
和
Figure GDA00026031168900000520
分别为非燃气常规机组有功出力的上限和下限；
Figure GDA00026031168900000521
和
Figure GDA00026031168900000522
分别为非燃气常规机组无功出力的上限和下限。
5)求解考虑天然气系统约束的等值最优潮流模型
式(11)-(20)构成求解考虑天然气系统运行约束的最优潮流的非线性最优化问题。通过内点法可直接求得电力系统各发电机的最优出力。
本发明采用上述技术方案后，主要有以下效果：
1、本发明构建的考虑天然气系统约束的等值最优潮流模型实现了在电力系统与天然气系统数据各自独立调度运行。在天然气系统与电力系统数据不能共享情况下，通过天然气系统约束等值模型将天然气系统的气源容量约束、节点气压约束、压缩机气压变比约束等保留到燃气轮机的有功出力约束中，并将该等值约束提供给电力系统具有不涉密的优势。
2、天然气系统与电力系统数据不能共享时，现有电-气互联系统采用交替迭代的优化的思想求解，存在着参数设定困难、收敛速度慢甚至不收敛等问题。而本专利提出的考虑天然气系统约束的等值最优潮流，通过约束等值模型将天然气系统约束保留到燃气轮机的有功出力约束中，避免了交替迭代存在的参数设定困难、收敛速度慢甚至不收敛等问题。
3、本发明构建的考虑天然气系统约束的等值最优潮流模型考虑了天然气系统的运行约束，在电力系统的运行调度中准确计及了天然气系统的影响，保证了优化计算结果的准确性，保证了电气互联系统的安全稳定运行。
附图说明
图1为IEEE14节点标准测试系统图
图2为NGS10节点系统结构图。
具体实施方式
下面结合实施例对本发明作进一步说明，但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的保护范围内。
本发明的技术路线可分为四步：一是天然气系统的约束等值模型的建立，二是天然气系统的约束等值模型的求解，三是考虑天然气系统约束的等值最优潮流模型的建立，四是考虑天然气系统约束的等值最优潮流模型的求解。首先，建立天然气系统的约束等式模型，以系统中所有燃气轮机消耗天然气流量的最大值作为目标函数，天然气系统节点流量平衡方程为等式约束，天然气节点气压、气源容量、压缩机气压变比的上限约束为不等式约束；其次，采用内点法对所建天然气系统的约束等式模型进行求解。再次，建立考虑天然气系统约束的等值最优潮流模型，以系统总发电成本作为目标函数，电力系统功率平衡方程为等式约束，天然气系统的等值约束、燃气轮机自身的容量约束、非燃气常规机组容量约束、节点电压约束和线路功率约束为不等式约束；最后，采用内点法对所建考虑天然气系统约束的等值最优潮流模型进行求解。结合具体实施方案，进一步说明如下：
(1)基础数据准备
图1的IEEE14节点标准系统中共有5台发电机，火电装机总容量为300MW，共有10个电负荷，功率为518MW。图2的NGS10节点天然气系统包括6条输气管道，3个燃气压缩机，2个气源站，气源容量共为 820MMCFD，4个自然气负荷L1-L4，共600MMCFD，分别是在节点6、7、9和10。天然气系统节点7和节点10分别提供电力系统节点1和节点3处燃气轮机的用气需求PGAS,1和PGAS,2。电力系统基础数据参见IEEE14节点标准系统。其中以燃气轮机和压缩机为例，燃气轮机参数取值表1所示，天然气系统中压缩机的参数取值如表2所示：
表1燃气轮机模型参数表
Figure GDA0002603116890000071
表2天然气压缩机参数表
Figure GDA0002603116890000081
(2)建立天然气系统的约束等值模型
IEEE14节点系统中共有5台发电机，其中G1和G3为燃气发电机，分别以G1消耗天然气流量、G3消耗天然气流量、G1和G3消耗天然气流量之和作为约束等值模型的目标函数；针对NGS10节点系统的所有节点可列写10个节点流量平衡方作为等式约束；考虑天然气系统节点气压、气源容量和压缩机气压变比的上下限，可列写10个天然气节点气压上下限约束、2个气源容量上下限约束、3个压缩机压变比上下限约束。其中，各变量的不等式约束的上下限如表3所示。
表3不等式约束上下限表
变量 Π<sub>m</sub>(psi) F<sub>G,m</sub>(MMCFD) R<sub>r</sub>
上限 640 700 1.6
下限 600 100 1.1
(3)求解天然气系统的约束等值模型
通过内点法求解天然气约束等值参数得到系统中各个燃气轮机组的有功出力约束以替代天然气系统的运行约束。天然气系统的约束等值计算结果：如表4所示，通过本发明提出的天然气系统的约束等值方法计算出了IEEE14-NGS10中燃气轮机的天然气耗量和有功出力最大值。可以看出PGAS,1和PGAS,2燃气轮机受天然气系统约束的有功出力上限为272MW，而自身的容量为320MW，如果不考虑天然气系统的约束可能会给电力系统的最优潮流带来误差。
表4燃气轮机天然气消耗和有功出力限值
Figure GDA0002603116890000091
(4)建立考虑天然气系统约束的等值最优潮流模型
IEEE14节点系统中共有5台发电机，其中G1和G3为燃气发电机，其他为非燃气发电机，以所有发电机的发电成本作为目标函数；针对IEEE14节点系统的所有节点可列写14个节点有功功率平衡方程和 14个无功功率平衡方程作为等式约束；考虑天然气系统等值约束(即上述表4中PGAS,1的有功出力上限约束235MW、PGAS,2的有功出力上限约束 234MW和PGAS,1、PGAS,2的有功出力之和上限约束272MW)、燃气轮机自身的容量约束(G1和G2相同)、非燃气常规机组容量约束、节点电压约束、线路功率约束，可列写3个天然气系统等值约束、2个燃气轮机自身的容量约束、3个非燃气常规机组容量约束、14个节点电压约束和20个线路功率约束。模型中变量的不等式约束上下限如表5所示：
表5不等式约束上下限表
Figure GDA0002603116890000092
(5)求解考虑天然气系统约束的等值最优潮流模型
本发明通过M0，M1和M2三种最优潮流方法计算IEEE14-NGS10系统， M0(全网统一进行最优潮流)的计算结果为准确结果，M1和M2的计算结果通过与M0的计算结果进行对比可以验证其有效性。
M0：全网统一进行最优潮流计算；
M1：本发明所提方法；
M2：电力系统独立进行最优潮流计算。
最优潮流计算结果：采用M0-M2的最优潮流方法对IEEE14-NGS10 系统进行计算，各发电机有功出力如表6所示，M2(电力系统独立进行最优潮流计算)的计算结果与M0相差很大，特别是发电机节点1处有48MW的绝对误差。M2方法中，电力系统节点1和节点3的燃气轮机的出力和为320MW，已经超出了272MW的上限，会影响电气互联系统的安全运行。M1(保留天然气系统约束的等值最优潮流方法)与M0(全网统一进行最优潮流计算)的计算结果几乎一致。因此，本发明所提出的考虑天然气系统约束的等值最优潮流能准确的计及天然气系统对电力系统的影响。
表6 M0-M2优化后的各发电机有功出力
发电机节点 M0(MW) M1(MW) M2(MW)
1 112.09 112.07 160
2 90 90 88.89
3 160 160 160
6 120 120 78.01
8 44.20 44.22 42.51
综上所述，本发明提出的一种考虑天然气系统约束的等值最优潮流计算方法。首先，考虑天然气系统的气源容量约束、节点气压约束、压缩机气压变比约束，提出天然气系统的约束等值模型；其次，采用内点法对所建天然气系统的约束等值模型进行求解，得到各燃气轮组的有功出力约束以替代天然气系统的运行约束；再次，基于约束等值模型提出考虑天然气系统约束的等值最优潮流模型；最后，采用内点法对所建考虑天然气系统约束的等值最优潮流模型进行求解，得到电力系统各发电机的最优出力。以IEEE14-NGS10电-气互联系统作为实施例实现仿真分析。通过仿真结果可知，本发明所提出的考虑天然气系统约束的等值最优潮流能准确的计及天然气系统对电力系统的影响。

Claims (1)

1.一种考虑天然气系统约束的等值最优潮流模型构建方法，其特征在于，包括以下步骤:
1)输入基础数据
输入电-气互联系统基础数据：电力子系统设备参数、天然气子系统设备参数和拓扑结构
2)建立天然气系统的约束等值模型
2-1)建立约束等值目标函数
根据天然气系统中燃气轮机个数N建立多个目标函数，首先建立任意一个燃气轮机消耗天然气流量最大的目标函数，一共N个；然后建立任意两个燃气轮机消耗天然气流总流量最大的目标函数，一共C2N个；以此类推，最后建立所有燃气轮机消耗天然气总流量最大的目标函数，一共CNN个；即由Ci，Cij，…，Call组成；
Figure FDA0002603116880000011
式中，在建立任意一个燃气轮机消耗天然气流量最大的目标函数时，Ci表示第i个燃气轮机消耗天然气流量的最大值，
Figure FDA0002603116880000012
表示第i个燃气轮机消耗天然气流量；在建立任意两个燃气轮机消耗天然气总流量最大的目标函数时，Cij为第i个和第j个燃气轮机消耗天然气总流量的最大值，
Figure FDA0002603116880000013
和
Figure FDA0002603116880000014
分别表示第i个和第j个燃气轮机消耗的天然气流量；在建立所有燃气轮机消耗天然气总流量最大的目标函数时，Call为所有燃气轮机消耗天然气总流量的最大值，
Figure FDA0002603116880000015
表示第i个燃气轮机消耗天然气流量；在求解燃气轮机消耗天然气流量最大值时，其他不相关燃气轮机消耗天然气流量为0；
2-2)建立约束等值等式约束方程
天然气系统流量平衡方程为式(2)：
FG,m-FGAS,m-FL,m-Fm＝0,m＝1,2,...,Nm (2)
式中，FG,m、FGAS,m、FL,m、Fm分别为天然气节点m的气源流量、燃气轮机消耗的流量、气负荷和注入流量，Nm为天然气系统节点个数；其中节点注入流量Fm通过式(3)求得；
Figure FDA0002603116880000021
式中，天然气节点n表示与天然气节点m相连的节点；Fmn、Fcom,mn、τcom,mn分别为管道流量、流过压缩机的流量和压缩机损耗，sgnc(m,n)和sgnτ(m,n)分别为压缩机流过流量和消耗流量的方向函数；这些量通过式(4)和(5)表示：
Figure FDA0002603116880000022
Figure FDA0002603116880000023
式中，Πm是天然气节点m的节点气压；kmn是与管道内径、长度、效率、压缩因子相关的常数；smn反应管道流量方向；Hcom,mn为压缩机消耗的功率；Bmn是压缩机系数；Zmn是压缩机压缩因子相关的常数；ατ、βτ、γτ为功率转化效率常数；
2-3)建立约束等值不等式约束方程
式(6)-(8)分别为天然气系统节点气压、气源容量、压缩机气压变比的上下限；
Figure FDA0002603116880000024
Figure FDA0002603116880000025
Figure FDA0002603116880000026
式中，Πm是天然气节点m的节点气压，FG,m为气源容量，Rr为压缩机气压变比；
Figure FDA0002603116880000027
和
Figure FDA0002603116880000028
分别为天然气节点气压的上限和下限；
Figure FDA0002603116880000029
和
Figure FDA00026031168800000210
分别为气源容量的上限和下限；
Figure FDA00026031168800000211
和
Figure FDA00026031168800000212
分别为压缩机压变比的上限和下限；Ns和Np分别为气源和压缩机的个数；
求解由公式(1)-(8)构成的优化模型，进而计算得到目标函数中的
Figure FDA00026031168800000213
通过燃气轮机耗量方程(9)计算出对应的
Figure FDA00026031168800000214
再通过公式(10)计算得到任意一个燃气轮机的出力上限，一共N个；得到任意两个燃气轮机总出力上限，一共C2N个；以此类推，得到所有燃气轮机总出力上限，一共CNN个；即由
Figure FDA0002603116880000031
组成；
Figure FDA0002603116880000032
Figure FDA0002603116880000033
式(10)中，在求解任意一个燃气轮机的出力上限时，
Figure FDA0002603116880000034
为第i个燃气轮机输出功率的最大值，
Figure FDA0002603116880000035
为第i个燃气轮机的输出功率；在求解任意两个燃气轮机总出力上限时，
Figure FDA0002603116880000036
为第i个和第j个燃气轮机总输出功率的最大值，
Figure FDA0002603116880000037
和
Figure FDA0002603116880000038
分别表示第i个和第j个燃气轮机输出功率；在求解所有燃气轮机总出力上限时，
Figure FDA0002603116880000039
为所有燃气轮机总输出功率的最大值，
Figure FDA00026031168800000310
为第i个燃气轮机的输出功率；GHV为天然气高热值；
3)求解天然气系统的约束等值模型
式(1)-(10)构成求解天然气约束等值参数的非线性最优化问题；通过内点法直接求得系统中各个燃气轮机组的有功出力约束以替代天然气系统的运行约束；
4)建立考虑天然气系统约束的等值最优潮流模型
4-1)等值最优潮流目标函数：
Figure FDA00026031168800000311
式中，fi(PGAS,i)为燃气轮机组的成本函数；fj(PGj)为非燃气发电机的成本函数；NG为非燃气发电机个数；
4-2)建立等值最优潮流等式约束方程：
考虑电力系统功率平衡方程(12)-(13)，构建如下等式约束：
Figure FDA00026031168800000312
Figure FDA00026031168800000313
式中，Pi和Qi分别为节点i的注入有功功率和无功功率；Ui是节点i的电压幅值；θij是节点i与j的电压相角差；Gij和Bij分别为节点导纳矩阵的第i行第j列元素的实部和虚部；NA为节点个数；
4-3)建立等值最优潮流不等式约束方程：
考虑天然气系统的等值约束
Figure FDA0002603116880000041
燃气轮机自身的容量约束、非燃气常规机组容量约束、节点电压约束和线路功率约束构建如下不等式约束：
Figure FDA0002603116880000042
Figure FDA0002603116880000043
Figure FDA0002603116880000044
Figure FDA0002603116880000045
Figure FDA0002603116880000046
Figure FDA0002603116880000047
Tl min≤Tl≤Tl max l＝1,2,...,NL (20)
其中，PGAS,i为燃气轮机i的有功功率，QGAS,i为燃气轮机i的无功功率，PGi为非燃气轮机i的有功功率，QGi为非燃气轮机i的无功功率，Ui是节点i的电压幅值，Tl是输电线路l流过的功率；NL为输电线路的总数；
Figure FDA0002603116880000048
和
Figure FDA0002603116880000049
分别为燃气轮机有功出力的上限和下限；
Figure FDA00026031168800000410
和
Figure FDA00026031168800000411
Figure FDA00026031168800000412
分别为燃气轮机无功出力的上限和下限；
Figure FDA00026031168800000413
和
Figure FDA00026031168800000414
分别为节点电压幅值的上限和下限；
Figure FDA00026031168800000415
和
Figure FDA00026031168800000416
为输电线路传输功率的上限和下限；
Figure FDA00026031168800000417
和
Figure FDA00026031168800000418
分别为非燃气常规机组有功出力的上限和下限；
Figure FDA00026031168800000419
和
Figure FDA00026031168800000420
分别为非燃气常规机组无功出力的上限和下限；
5)求解考虑天然气系统约束的等值最优潮流模型
式(11)-(20)构成求解考虑天然气系统运行约束的最优潮流的非线性最优化问题；通过内点法直接求得电力系统各发电机的最优出力。
CN201711172326.2A 2017-11-20 2017-11-20 考虑天然气系统约束的等值最优潮流模型构建方法 Expired - Fee Related CN107947245B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711172326.2A CN107947245B (zh) 2017-11-20 2017-11-20 考虑天然气系统约束的等值最优潮流模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711172326.2A CN107947245B (zh) 2017-11-20 2017-11-20 考虑天然气系统约束的等值最优潮流模型构建方法

Publications (2)

Publication Number Publication Date
CN107947245A CN107947245A (zh) 2018-04-20
CN107947245B true CN107947245B (zh) 2020-10-09

Family

ID=61929776

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711172326.2A Expired - Fee Related CN107947245B (zh) 2017-11-20 2017-11-20 考虑天然气系统约束的等值最优潮流模型构建方法

Country Status (1)

Country Link
CN (1) CN107947245B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108921404B (zh) * 2018-06-15 2020-05-15 浙江大学 一种基于通用生成函数的电-气耦合系统的可靠性评估方法
CN108964143B (zh) * 2018-06-28 2020-10-09 重庆大学 一种电-气综合能源系统的天然气网络静态等值模型
CN108875269B (zh) * 2018-07-09 2022-06-10 重庆大学 一种考虑电力系统多平衡机与天然气系统慢动态特性的电-气互联系统多时段能流计算方法
CN109066695A (zh) * 2018-09-29 2018-12-21 广东电网有限责任公司 一种两阶段线性化电气最优能流计算方法
CN109409705B (zh) * 2018-10-11 2021-11-23 河海大学 一种区域综合能源系统的多目标优化调度方法
CN109902854B (zh) * 2019-01-11 2020-11-27 重庆大学 电-气互联系统最优潮流全线性模型的构建方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104734155B (zh) * 2015-04-14 2017-03-01 河海大学 一种获取电‑气互联能源系统可用输电能力的方法
CN105005940A (zh) * 2015-07-09 2015-10-28 河海大学 计及相关性的电-气互联系统概率最优潮流计算方法
CN106096777B (zh) * 2016-06-12 2019-05-21 清华大学 一种电-气耦合多能流系统的优化调度方法

Also Published As

Publication number Publication date
CN107947245A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN107947245B (zh) 考虑天然气系统约束的等值最优潮流模型构建方法
CN107291990B (zh) 基于电-气互联综合能源系统暂态模型的能量流仿真方法
CN105046369B (zh) 一种基于能源中心的电气混联系统建模和优化调度方法
CN109830955B (zh) 考虑柔性约束与全周期成本的电-气配网柔性规划方法
CN109492325B (zh) 一种基于扩展能源集线器的多能耦合系统潮流分析方法
CN105576657A (zh) 能源互联下主动配网与输电网集成系统的随机模糊连续潮流建模及算法
CN111950122A (zh) 园区综合能源系统运行优化方法
CN110728032B (zh) 一种计及环网的电-热互联综合能源系统快速潮流计算方法
Qu et al. Distributed optimal control of energy hubs for micro-integrated energy systems
CN112035984A (zh) 一种电-气-储区域综合能源系统协同规划方法
Huang et al. A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale
CN110647040B (zh) 一种综合能源系统的安全控制方法及装置
CN115392035A (zh) 电-气互联综合能源系统优化调度的改进方法
Dong et al. Hierarchical multi-objective planning for integrated energy systems in smart parks considering operational characteristics
Zhang et al. Steady-state security region of integrated energy system considering thermal dynamics
Wang et al. Decoupled optimization of integrated energy system considering CHP plant based on energy hub model
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN114221346A (zh) 一种综合能源系统潮流计算方法
Varganova et al. Energy-efficient distribution of heat between the boiler units of industrial power palnts with the use of computers
Yang et al. Probabilistic power and gas flow analysis for electricity-gas coupled networks considering uncertainties in pipeline parameters
CN113011655A (zh) 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法
CN111783309A (zh) 基于内部守恒的蒸汽供热网络动态仿真方法
CN111310310A (zh) 一种用于量调节的热力系统静态潮流快速解耦计算方法
CN114549232A (zh) 电-热综合能源系统混合能流计算方法
CN113901653A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201009

Termination date: 20211120