CN113517697A - 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法 - Google Patents

一种基于全纯嵌入法的气电联合网络稳态潮流计算方法 Download PDF

Info

Publication number
CN113517697A
CN113517697A CN202110370536.2A CN202110370536A CN113517697A CN 113517697 A CN113517697 A CN 113517697A CN 202110370536 A CN202110370536 A CN 202110370536A CN 113517697 A CN113517697 A CN 113517697A
Authority
CN
China
Prior art keywords
gas
power
representing
natural gas
node
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.)
Pending
Application number
CN202110370536.2A
Other languages
English (en)
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202110370536.2A priority Critical patent/CN113517697A/zh
Publication of CN113517697A publication Critical patent/CN113517697A/zh
Pending legal-status Critical Current

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/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
    • 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，过程如下：输入电力系统、天然气系统的参数和运行数据；用全纯函数代表气电联合网络潮流方程中的未知变量；构建含有全纯函数及复数变量的气电联合网络潮流方程；计算全纯函数中未知的各多项式系数；对代表未知变量的全纯函数进行解析延拓，获取未知变量的数值解；将各未知变量数值解带入潮流方程，检查是否满足电力系统潮流约束，若满足约束则输出气电联合网络潮流解，若不满足约束将增加原全纯函数的级数，重复以上过程直至满足约束。本发明保证气电联合网络中的电力系统计算具有更好的收敛性，所获取解一定为高压解，避免常见方法的多解问题。

Description

一种基于全纯嵌入法的气电联合网络稳态潮流计算方法
技术领域
本发明涉及气电联合网络稳态运行与分析技术领域，具体涉及一种基于全纯嵌入法的气电联合网络稳态潮流计算方法。
背景技术
由于可再生能源份额增加导致的对于电力系统的更高灵活度的要求，以及减少发电环节污染排放的要求，越来越多的天然气发电机组接入了电力系统。这加深了电力系统与天然气系统的耦合关系，也使得电力系统和天然气系统常作为一个联合网络(即“气电联合网络”)进行控制、优化、分析。对气电联合网络进行统一的潮流计算及分析，可以有效地获取该系统的运行状态，为系统运行人员提供了基本数据。因此，如何快速、准确地计算气电联合网络的潮流是一个具有重要意义的研究课题。
对于气电联合网络的统一潮流计算，已有方法为牛顿拉夫逊法。然而，牛顿拉夫逊法的收敛能力比较有限：该方法对于计算的初始值比较敏感，初始值的选取不佳会直接增加计算过程的循环次数和计算耗时，不同的初始值会产生不同的潮流解；潮流计算过程中出现不收敛的情况时，无法判断不收敛的原因是系统无解还是牛顿拉夫逊法本身能力有限，令运行人员无法判断系统实际情况。以上问题导致现有的气电联合网络的潮流计算方法存在一定的失败风险，难以保证迅速、准确的潮流计算效果。
发明内容
本发明的目的是为了解决现有气电联合网络潮流计算技术的缺陷，提供一种基于全纯嵌入法的气电联合网络的潮流计算方法，通过重构气电联合网络的潮流方程、采用全纯函数代替变量、求解全纯函数的未知系数以及对全纯函数进行解析延拓，实现准确、快速的气电联合网络潮流计算。
本发明的目的可以通过采取如下技术方案达到：
一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，所述潮流计算方法包括以下步骤：
S1、输入电力系统、天然气系统的参数和运行数据；
S2、用全纯函数代表气电联合网络潮流方程中的未知变量；
S3、构建含有全纯函数及复数变量的气电联合网络潮流方程；
S4、计算全纯函数中的未知系数；
S5、对代表未知变量的全纯函数进行解析延拓，获取未知变量的数值解；
S6、将各未知变量的数值解带入气电联合网络潮流方程，检查是否满足气电联合网络的运行约束，若满足约束则输出气电联合网络潮流解，若不满足约束将增加原全纯函数的阶数，重复步骤S2到S6直至满足气电联合网络潮流约束。
进一步地，所述电力系统的参数包括电力系统网络拓扑结构、线路导纳参数，所述电力系统的运行数据包括各节点的电力负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角。
进一步地，所述天然气系统的参数包括天然气系统网络拓扑结构、管道压降系数、电驱动压缩机耗电量和天然气驱动压缩机耗气量计算公式中各项参数、燃气轮机气电转换效率，所述天然气系统的运行数据包括各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比。
进一步地，所述气电联合网络的未知变量包括：电力系统中各节点的电压向量、电力系统中各燃气轮机所消耗的天然气气量、电力系统各燃气轮机发电量、电力系统为天然气系统中各电驱动的压缩机所提供电量、各天然气驱动的压缩机耗气量、天然气系统各管道中的天然气流量；
所述步骤S2过程如下：
在电力系统中，用含有复数变量s的全纯函数V(s)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(s)写为以下形式：
Figure BDA0003009107640000031
在天然气系统中，用含有复数变量s的全纯函数代表各未知数和部分运行参数，所使用的全纯函数包括：代表天然气系统各管道中的天然气流量的FP(s)、代表电力系统燃气轮机所消耗天然气气量的FGT(s)、代表电力系统各燃气轮机发电量的PGT(s)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的PEC(s)、代表各天然气驱动的压缩机耗气量的FGC(s)，代表以上全纯函数写为幂级数形式，总共N阶，n代表阶数，fn、
Figure BDA0003009107640000032
Figure BDA0003009107640000033
分别为幂级数FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)第n阶项的系数，各全纯函数的幂级数表达式写为以下形式：
Figure BDA0003009107640000034
Figure BDA0003009107640000035
Figure BDA0003009107640000036
Figure BDA0003009107640000037
Figure BDA0003009107640000038
进一步地，，所述步骤S3过程如下：
利用式(1)至(6)中的全纯函数及复数变量s表达气电联合网络的潮流方程组中的各项方程如下：
Figure BDA0003009107640000041
Figure BDA0003009107640000042
Figure BDA0003009107640000043
Figure BDA0003009107640000044
Figure BDA0003009107640000045
Figure BDA0003009107640000046
Figure BDA0003009107640000047
Figure BDA0003009107640000048
式(7)至(9)中，
Figure BDA0003009107640000049
为代表电力系统PQ节点电压相量的全纯函数的共轭，计算公式为
Figure BDA00030091076400000410
Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(s)为代表电力系统所有节点电压相量的全纯函数；
Figure BDA00030091076400000411
代表电力系统PQ节点电压相量初始值的共轭值，ck代表电力系统所有节点电压相量初始值，
Figure BDA00030091076400000412
代表电力系统PQ节点功率注入量的共轭值，
Figure BDA00030091076400000413
代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(s)和
Figure BDA00030091076400000414
分别代表电力系统发电机节点电压相量的全纯函数及其共轭，
Figure BDA0003009107640000051
计算公式为
Figure BDA0003009107640000052
cigen和cigen*分别代表电力系统发电机节点电压相量初始值及其共轭，Vigen,sp代表电力系统发电机节点电压相量预定值；Vipv(s)和
Figure BDA0003009107640000053
分别代表电力系统PV节点电压相量的全纯函数及其共轭，
Figure BDA0003009107640000054
计算公式为
Figure BDA0003009107640000055
Yipv,k和
Figure BDA0003009107640000056
分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，cipv和
Figure BDA0003009107640000057
分别代表电力系统PV节点电压相量初始值及其共轭，代表电力系统PV节点的有功注入功率，
Figure BDA0003009107640000058
代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
式(10)至式(14)中，Aing,bng代表天然气系统中节点-管道关联矩阵对应节点ing、管道bng的元素，
Figure BDA0003009107640000059
代表天然气系统管道bng内天然气流量的全纯函数，mbng代表天然气系统中管道bng内天然气流量初始值，
Figure BDA00030091076400000510
和
Figure BDA00030091076400000511
分别代表天然气系统中节点ing处用气量和供气量，
Figure BDA00030091076400000512
代表电力系统第igt台燃气轮机所消耗的天然气气量，
Figure BDA00030091076400000513
为天然气系统节点-燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素，
Figure BDA00030091076400000514
代表天然气系统中节点ing处向天然气驱动的压缩机耗气量的全纯函数；Bl,bng代表天然气系统中环路-管道关联矩阵对应环路l、节点bng的元素，Cbng代表天然气系统管道bng内的管道压降系数；aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数，
Figure BDA00030091076400000515
代表电力系统中第igt台燃气轮机的发电量；
Figure BDA00030091076400000516
代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数，
Figure BDA00030091076400000517
代表电驱动压缩机耗电量计算参数，
Figure BDA00030091076400000518
代表天然气系统中管道iecg内天然气流量；
Figure BDA00030091076400000519
代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数，
Figure BDA00030091076400000520
代表天然气系统中管道igcb内天然气流量。
进一步地，所述步骤S4过程如下：
S401、将式(1)到(6)代入式(7)到(14)，获得含有全纯函数的非线性方程组；
S402、提取出步骤S401得到的非线性方程组中所有的s项，获得仅包括全纯函数一阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)的一阶项系数；
S403、将已获得的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)一阶项系数代入步骤S401得到的非线性方程组，提取出非线性方程组中所有的s2项，获得仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)二阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的二阶项系数；
S404、将已经计算出的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的系数代入非线性方程组，提取更高阶的s项，获取仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)中该阶数系数的线性方程组进行求解，直至完成第N阶项系数的计算。
进一步地，所述步骤S5过程如下：
S501、取正整数M和L，其中，
Figure BDA0003009107640000061
用
Figure BDA0003009107640000062
表示气电联合网络的各未知变量对应的全纯函数V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)，i为气电联合网络未知变量的编号，j代表阶数，xi,j为Xi(s)中第j阶项系数；
S502、为每个气电联合网络未知变量设置未知幂级数函数
Figure BDA0003009107640000063
其中j1代表阶数，yi,j1为Yi(s)中第j1阶项系数，yi,0＝1，利用Yi(s)和Xi(s)构造以下方程组：
Figure BDA0003009107640000071
通过求解公式(15)可获得Yi(s)的各阶项系数yi,1、…、yi,j1、…、yi,M；
S503、为每个气电联合网络未知变量i设置未知幂级数函数
Figure BDA0003009107640000072
其中j2代表阶数，zi,j2为Zi(s)中第j2阶系数，利用
Figure BDA0003009107640000073
和
Figure BDA0003009107640000074
计算Zi(s)的各阶项系数zi,0、…、zi,j2、…、zi,L：
Figure BDA0003009107640000075
S504、令s取值为1，将公式(15)、(16)获得的Yi(s)和Zi(s)的各阶项系数代入公式(17)，得到第i个变量的数值解
Figure BDA0003009107640000076
Figure BDA0003009107640000077
进一步地，所述步骤S6过程如下：
S601、将式(17)获得的各未知变量数值解
Figure BDA0003009107640000078
带入潮流方程公式(7)到(14)：
S602、计算潮流方程公式(7)到(14)左右两侧的差值，若差值低于预设标准，则输出潮流解，计算结束；若高于预设标准，则增加各变量对应的全纯函数的项数，重复步骤S2到S6，直至方程(7)到(14)左右两侧的差值低于预设标准，或步骤重复次数超过预设的最大循环次数。
本发明相对于现有技术具有如下的优点及效果：
1、本发明提出的气电联合网络稳态潮流计算方法，使用全纯函数与复数变量重新构造气电联合网络的潮流方程组，可保证气电联合网络中的电力系统计算所获取解一定为高压解，避免了常用的牛顿拉夫逊方法的多解问题；本方法收敛性能更好，能够更好地区分当前气电联合网络潮流方程是否存在可行解。
2、相对于已有方法，本发明需要的迭代计算次数少，和已有算法相比，该方法的在接近潮流实际解时仍然能保持较快的收敛速度，计算耗时更低；计算过程中该方法需要求逆的矩阵个数明显降低，计算量低，计算所需的存储空间小。
3、本发明提出的气电联合网络稳态潮流计算方法在气电联合网络的规划、优化、控制及分析等各方面具有广泛的使用空间，适用于各种实际运行的电力系统模型，通用性强。
附图说明
图1是本发明公开的一种基于全纯嵌入法的气电联合网络的潮流计算方法的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
如图1所示，本实施例提供一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，该方法用全纯函数代替原系统潮流方程中的变量，通过构造新的潮流方程求解全纯函数的未知系数并进行解析延拓，求得系统潮流方程中各未知变量的解，最终获取系统的潮流状态数据。具体过程如下：
S1、输入电力系统、天然气系统的参数和运行数据；
其中，电力系统的参数包括电力系统网络拓扑结构、线路导纳参数，电力系统的运行数据包括各节点的电力负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角等。天然气系统的参数包括天然气系统网络拓扑结构、管道压降系数、压缩机耗电量和耗气量计算公式中各项参数、燃气机组气电转换效率，电力系统的运行数据包括各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比等。
S2、用全纯函数代表气电联合网络潮流方程中的未知变量及部分已知量，具体如下：
在电力系统中，用含有复数变量s的全纯函数V(s)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(s)写为以下形式：
Figure BDA0003009107640000091
在天然气系统中，用含有复数变量s的全纯函数代表各未知数和部分运行参数，所使用的全纯函数包括：代表天然气系统各管道中的天然气流量的FP(s)、代表电力系统燃气轮机所消耗天然气气量的FGT(s)、代表电力系统各节点处燃气轮机发电量的PGT(s)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的PEC(s)、代表各天然气驱动的压缩机耗气量的FGC(s)，代表以上全纯函数写为幂级数形式，总共N阶，n代表阶数，fn、
Figure BDA0003009107640000092
分别为幂级数FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)第n阶项的系数，各全纯函数的幂级数表达式写为以下形式：
Figure BDA0003009107640000093
Figure BDA0003009107640000094
Figure BDA0003009107640000101
Figure BDA0003009107640000102
Figure BDA0003009107640000103
S3、构建含有全纯函数及复数变量的气电联合网络潮流方程，具体过程如下：
S301、利用式(1)至(7)中的全纯函数及复数变量s表达气电联合网络的潮流方程组中的各项方程
Figure BDA0003009107640000104
Figure BDA0003009107640000105
Figure BDA0003009107640000106
Figure BDA0003009107640000107
Figure BDA0003009107640000108
Figure BDA0003009107640000109
Figure BDA00030091076400001010
Figure BDA00030091076400001011
式(7)至(9)中，
Figure BDA00030091076400001012
为代表电力系统PQ节点电压的全纯函数的共轭，计算公式为
Figure BDA0003009107640000111
Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(s)为代表电力系统所有节点电压的全纯函数；
Figure BDA0003009107640000112
代表电力系统PQ节点电压初始值的共轭值，ck代表电力系统所有节点电压初始值，
Figure BDA0003009107640000113
代表电力系统PQ节点功率注入量的共轭值，
Figure BDA0003009107640000114
代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(s)和
Figure BDA0003009107640000115
分别代表电力系统发电机节点电压的全纯函数及其共轭，
Figure BDA0003009107640000116
计算公式为
Figure BDA0003009107640000117
cigen和cigen*分别代表电力系统发电机节点电压初始值及其共轭，Vigen,sp代表电力系统发电机节点电压预定值；Vipv(s)和
Figure BDA0003009107640000118
分别代表电力系统PV节点电压的全纯函数及其共轭，
Figure BDA0003009107640000119
计算公式为
Figure BDA00030091076400001110
Yipv,k和
Figure BDA00030091076400001111
分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，cipv和
Figure BDA00030091076400001112
分别代表电力系统PV节点电压初始值及其共轭，代表电力系统PV节点的有功注入功率，
Figure BDA00030091076400001113
代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
式(10)至式(14)中，Aing,bng代表天然气系统中节点-管道关联矩阵对应节点ing、管道bng的元素，
Figure BDA00030091076400001114
代表天然气系统管道bng内天然气流量的全纯函数，mbng代表天然气系统中管道bng内天然气流量初始值，
Figure BDA00030091076400001115
和
Figure BDA00030091076400001116
分别代表天然气系统中节点ing处用气量和供气量，
Figure BDA00030091076400001117
代表电力系统第igt台燃气轮机所消耗的天然气气量，
Figure BDA00030091076400001118
为天然气系统节点-燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素，
Figure BDA00030091076400001119
代表天然气系统中节点ing处向天然气驱动的压缩机耗气量的全纯函数；Bl,bng代表天然气系统中环路-管道关联矩阵对应环路l、节点bng的元素，Cbng代表天然气系统管道bng内的管道压降系数；aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数，
Figure BDA0003009107640000121
代表电力系统中第igt台燃气轮机的发电量；
Figure BDA0003009107640000122
代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数，
Figure BDA0003009107640000123
代表电驱动压缩机耗电量计算参数，
Figure BDA0003009107640000124
代表天然气系统中管道iecg内天然气流量；
Figure BDA0003009107640000125
代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数，
Figure BDA0003009107640000126
代表天然气系统中管道igcb内天然气流量。
S4、计算全纯函数中未知的各多项式系数，具体过程如下：
S401、将式(1)到(6)代入式(7)到(14)，获得含有幂级数的非线性方程组；
S402、提取出S401得到的非线性方程组中所有的s项，获得仅包括全纯函数一阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)一阶项系数；
S403、将已获得的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)一阶项系数代入步骤S401得到的非线性方程组，提取出非线性方程组中所有的s2项，获得仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)二阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的二阶项系数；
S404、将已经计算出的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的系数代入非线性方程组，提取更高阶的s项，获取仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)中该阶数系数的线性方程组进行求解。该步骤持续进行，直至完成第N阶项系数的计算。
S5、对代表未知变量的全纯函数进行解析延拓，获取未知变量的数值解，具体过程如下：
S501、取正整数M和L，其中，
Figure BDA0003009107640000131
用
Figure BDA0003009107640000132
表示气电联合网络的各未知变量对应的全纯函数V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)，i为气电联合网络未知变量的编号，j代表阶数，xi,j为Xi(s)中第j阶项系数；
S502、为每个气电联合网络未知变量设置未知幂级数函数
Figure BDA0003009107640000133
其中j1代表阶数，yi,j1为Yi(s)中第j1阶项系数，yi,0＝1，利用Yi(s)和Xi(s)构造以下方程组：
Figure BDA0003009107640000134
通过求解公式(15)可获得Yi(s)的各阶项系数yi,1、…、yi,j1、…、yi,M；
S503、为每个气电联合网络未知变量i设置未知幂级数函数
Figure BDA0003009107640000135
其中j2代表阶数，zi,j2为Zi(s)中第j2阶系数，利用
Figure BDA0003009107640000136
和
Figure BDA0003009107640000137
计算Zi(s)的各阶项系数zi,0、…、zi,j2、…、zi,L：
Figure BDA0003009107640000138
S504、令s取值为1，将公式(15)、(16)获得的Yi(s)和Zi(s)的各阶项系数代入公式(17)，得到第i个变量的数值解
Figure BDA0003009107640000139
Figure BDA00030091076400001310
S6将各未知变量数值解带入潮流方程，检查是否满足系统潮流约束，具体步骤如下：
S601、将式(17)获得的各未知变量数值解
Figure BDA00030091076400001311
带入潮流方程公式(8)到(15)：
S602、计算潮流方程公式(7)到(14)左右两侧的差值，若差值低于预设标准，则输出潮流解，计算结束；若高于预设标准，则增加各变量对应的全纯函数的项数，重复步骤S2到S6，直至潮流方程公式(7)到(14)左右两侧的差值低于预设标准，或步骤重复次数超过预设的最大循环次数。
综上所述，本发明通过重构气电联合网络潮流方程组、采用全纯函数代替变量、求解全纯函数的未知系数以及对全纯函数进行解析延拓，能够快速、准确地完成气电联合网络潮流计算，并可用于不同系统，具有较强的实用性。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (8)

1.一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述潮流计算方法包括以下步骤：
S1、输入电力系统、天然气系统的参数和运行数据；
S2、用全纯函数代表气电联合网络潮流方程中的未知变量；
S3、构建含有全纯函数及复数变量的气电联合网络潮流方程；
S4、计算全纯函数中的未知系数；
S5、对代表未知变量的全纯函数进行解析延拓，获取未知变量的数值解；
S6、将各未知变量的数值解带入气电联合网络潮流方程，检查是否满足气电联合网络的运行约束，若满足约束则输出气电联合网络潮流解，若不满足约束将增加原全纯函数的阶数，重复步骤S2到S6直至满足气电联合网络潮流约束。
2.根据权利要求1所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述电力系统的参数包括电力系统网络拓扑结构、线路导纳参数，所述电力系统的运行数据包括各节点的电力负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角。
3.根据权利要求1所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述天然气系统的参数包括天然气系统网络拓扑结构、管道压降系数、电驱动压缩机耗电量和天然气驱动压缩机耗气量计算公式中各项参数、燃气轮机气电转换效率，所述天然气系统的运行数据包括各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比。
4.根据权利要求1所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述气电联合网络的未知变量包括：电力系统中各节点的电压向量、电力系统中各燃气轮机所消耗的天然气气量、电力系统各燃气轮机发电量、电力系统为天然气系统中各电驱动的压缩机所提供电量、各天然气驱动的压缩机耗气量、天然气系统各管道中的天然气流量；
所述步骤S2过程如下：
在电力系统中，用含有复数变量s的全纯函数V(s)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(s)写为以下形式：
Figure FDA0003009107630000021
在天然气系统中，用含有复数变量s的全纯函数代表各未知数和部分运行参数，所使用的全纯函数包括：代表天然气系统各管道中的天然气流量的FP(s)、代表电力系统燃气轮机所消耗天然气气量的FGT(s)、代表电力系统各燃气轮机发电量的PGT(s)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的PEC(s)、代表各天然气驱动的压缩机耗气量的FGC(s)，代表以上全纯函数写为幂级数形式，总共N阶，n代表阶数，fn、
Figure FDA0003009107630000022
Figure FDA0003009107630000023
分别为幂级数FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)第n阶项的系数，各全纯函数的幂级数表达式写为以下形式：
Figure FDA0003009107630000031
Figure FDA0003009107630000032
Figure FDA0003009107630000033
Figure FDA0003009107630000034
Figure FDA0003009107630000035
5.根据权利要求1所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述步骤S3过程如下：
利用式(1)至(6)中的全纯函数及复数变量s表达气电联合网络的潮流方程组中的各项方程如下：
Figure FDA0003009107630000036
Figure FDA0003009107630000037
Figure FDA0003009107630000038
Figure FDA0003009107630000039
Figure FDA00030091076300000310
Figure FDA0003009107630000041
Figure FDA0003009107630000042
Figure FDA0003009107630000043
式(7)至(9)中，
Figure FDA0003009107630000044
为代表电力系统PQ节点电压相量的全纯函数的共轭，计算公式为
Figure FDA0003009107630000045
Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(s)为代表电力系统所有节点电压相量的全纯函数；
Figure FDA0003009107630000046
代表电力系统PQ节点电压相量初始值的共轭值，ck代表电力系统所有节点电压相量初始值，
Figure FDA0003009107630000047
代表电力系统PQ节点功率注入量的共轭值，
Figure FDA0003009107630000048
代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(s)和
Figure FDA0003009107630000049
分别代表电力系统发电机节点电压相量的全纯函数及其共轭，
Figure FDA00030091076300000410
计算公式为
Figure FDA00030091076300000411
cigen和cigen*分别代表电力系统发电机节点电压相量初始值及其共轭，Vigen,sp代表电力系统发电机节点电压相量预定值；Vipv(s)和
Figure FDA00030091076300000412
分别代表电力系统PV节点电压相量的全纯函数及其共轭，
Figure FDA00030091076300000413
计算公式为
Figure FDA00030091076300000414
Yipv,k和
Figure FDA00030091076300000415
分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，cipv和
Figure FDA00030091076300000416
分别代表电力系统PV节点电压相量初始值及其共轭，代表电力系统PV节点的有功注入功率，
Figure FDA00030091076300000417
代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
式(10)至式(14)中，Aing,bng代表天然气系统中节点-管道关联矩阵对应节点ing、管道bng的元素，
Figure FDA00030091076300000418
代表天然气系统管道bng内天然气流量的全纯函数，mbng代表天然气系统中管道bng内天然气流量初始值，
Figure FDA00030091076300000419
和
Figure FDA0003009107630000051
分别代表天然气系统中节点ing处用气量和供气量，
Figure FDA0003009107630000052
代表电力系统第igt台燃气轮机所消耗的天然气气量，
Figure FDA0003009107630000053
为天然气系统节点-燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素，
Figure FDA0003009107630000054
代表天然气系统中节点ing处向天然气驱动的压缩机耗气量的全纯函数；Bl,bng代表天然气系统中环路-管道关联矩阵对应环路l、节点bng的元素，Cbng代表天然气系统管道bng内的管道压降系数；aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数，
Figure FDA0003009107630000055
代表电力系统中第igt台燃气轮机的发电量；
Figure FDA0003009107630000056
代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数，
Figure FDA0003009107630000057
代表电驱动压缩机耗电量计算参数，
Figure FDA0003009107630000058
代表天然气系统中管道iecg内天然气流量；
Figure FDA0003009107630000059
代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数，
Figure FDA00030091076300000510
代表天然气系统中管道igcb内天然气流量。
6.根据权利要求5所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述步骤S4过程如下：
S401、将式(1)到(6)代入式(7)到(14)，获得含有全纯函数的非线性方程组；
S402、提取出步骤S401得到的非线性方程组中所有的s项，获得仅包括全纯函数一阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)的一阶项系数；
S403、将已获得的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)一阶项系数代入步骤S401得到的非线性方程组，提取出非线性方程组中所有的s2项，获得仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)二阶项系数的线性方程组，通过求解该线性方程组获得V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的二阶项系数；
S404、将已经计算出的V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)的系数代入非线性方程组，提取更高阶的s项，获取仅包括V(s)、FP(s)、FGT(s)、PGT(s)、PEC(s)、FGC(s)中该阶数系数的线性方程组进行求解，直至完成第N阶项系数的计算。
7.根据权利要求6所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述步骤S5过程如下：
S501、取正整数M和L，其中，
Figure FDA0003009107630000061
用
Figure FDA0003009107630000062
表示气电联合网络的各未知变量对应的全纯函数V(s)、FP(s)、FGT(s)、FEC(s)、FGC(s)，i为气电联合网络未知变量的编号，j代表阶数，xi,j为Xi(s)中第j阶项系数；
S502、为每个气电联合网络未知变量设置未知幂级数函数
Figure FDA0003009107630000063
其中j1代表阶数，yi,j1为Yi(s)中第j1阶项系数，yi,0＝1，利用Yi(s)和Xi(s)构造以下方程组：
Figure FDA0003009107630000064
通过求解公式(15)可获得Yi(s)的各阶项系数yi,1、…、yi,j1、…、yi,M；
S503、为每个气电联合网络未知变量i设置未知幂级数函数
Figure FDA0003009107630000065
其中j2代表阶数，zi,j2为Zi(s)中第j2阶系数，利用
Figure FDA0003009107630000071
和
Figure FDA0003009107630000072
计算Zi(s)的各阶项系数zi,0、…、zi,j2、…、zi,L：
Figure FDA0003009107630000073
S504、令s取值为1，将公式(15)、(16)获得的Yi(s)和Zi(s)的各阶项系数代入公式(17)，得到第i个变量的数值解
Figure FDA0003009107630000074
Figure FDA0003009107630000075
8.根据权利要求7所述的一种基于全纯嵌入法的气电联合网络稳态潮流计算方法，其特征在于，所述步骤S6过程如下：
S601、将式(17)获得的各未知变量数值解
Figure FDA0003009107630000076
带入潮流方程公式(7)到(14)：
S602、计算潮流方程公式(7)到(14)左右两侧的差值，若差值低于预设标准，则输出潮流解，计算结束；若高于预设标准，则增加各变量对应的全纯函数的项数，重复步骤S2到S6，直至公式(7)到(14)左右两侧的差值低于预设标准，或步骤重复次数超过预设的最大循环次数。
CN202110370536.2A 2021-04-07 2021-04-07 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法 Pending CN113517697A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110370536.2A CN113517697A (zh) 2021-04-07 2021-04-07 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110370536.2A CN113517697A (zh) 2021-04-07 2021-04-07 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法

Publications (1)

Publication Number Publication Date
CN113517697A true CN113517697A (zh) 2021-10-19

Family

ID=78062342

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110370536.2A Pending CN113517697A (zh) 2021-04-07 2021-04-07 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法

Country Status (1)

Country Link
CN (1) CN113517697A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114336635A (zh) * 2022-01-05 2022-04-12 国网内蒙古东部电力有限公司通辽供电公司 基于常项值和先验节点的全纯嵌入潮流计算方法、装置
WO2024103221A1 (zh) * 2022-11-14 2024-05-23 中山大学 一种快速灵活全纯嵌入式电网越限预测及稳定性评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10594136B1 (en) * 2015-07-02 2020-03-17 Antonio Trias System and method to simulate controls in the holomorphic embedding loadflow calculation
CN110957728A (zh) * 2018-09-27 2020-04-03 国网甘肃省电力公司电力科学研究院 一种三相有源配电网的全纯嵌入潮流方法及装置
CN111082427A (zh) * 2020-01-07 2020-04-28 湘潭大学 一种基于全纯函数的微电网潮流计算方法
CN111950135A (zh) * 2020-07-28 2020-11-17 重庆大学 一种基于网络解耦的电-气互联系统概率能流计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10594136B1 (en) * 2015-07-02 2020-03-17 Antonio Trias System and method to simulate controls in the holomorphic embedding loadflow calculation
CN110957728A (zh) * 2018-09-27 2020-04-03 国网甘肃省电力公司电力科学研究院 一种三相有源配电网的全纯嵌入潮流方法及装置
CN111082427A (zh) * 2020-01-07 2020-04-28 湘潭大学 一种基于全纯函数的微电网潮流计算方法
CN111950135A (zh) * 2020-07-28 2020-11-17 重庆大学 一种基于网络解耦的电-气互联系统概率能流计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HSIAO-DONG CHIANG等: "A Novel Fast and Flexible Holomorphic Embedding Power Flow Method", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
TONG ZHANG等: "Power Flow Analysis of Integrated Gas and Electricity Systems Using the Fast and Flexible Holomorphic Embedding Method", 《2020 IEEE POWER & ENERGY SOCIETY GENERAL MEETING (PESGM)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114336635A (zh) * 2022-01-05 2022-04-12 国网内蒙古东部电力有限公司通辽供电公司 基于常项值和先验节点的全纯嵌入潮流计算方法、装置
CN114336635B (zh) * 2022-01-05 2023-07-25 国网内蒙古东部电力有限公司通辽供电公司 基于常项值和先验节点的全纯嵌入潮流计算方法、装置
WO2024103221A1 (zh) * 2022-11-14 2024-05-23 中山大学 一种快速灵活全纯嵌入式电网越限预测及稳定性评估方法

Similar Documents

Publication Publication Date Title
Ziogou et al. Optimal production of renewable hydrogen based on an efficient energy management strategy
Chen et al. Fast decoupled multi-energy flow calculation for integrated energy system
CN108416507B (zh) 面向电-气耦合综合能源系统的静态灵敏度分析方法
CN101976837A (zh) 一种电力系统潮流无解调整方法
CN113517697A (zh) 一种基于全纯嵌入法的气电联合网络稳态潮流计算方法
CN112736923A (zh) 一种考虑电网故障影响的天然气网连锁故障评估控制方法
Suyono et al. Analysis of power losses due to distributed generation increase on distribution system
CN110429591B (zh) 一种基于电力系统时序耦合性的输电网利用率评估方法
Yun et al. Multi-time collaborative restoration for integrated electrical-gas distribution system based on rolling optimization
Webster et al. Implementation of a non-discretized multiphysics PEM electrolyzer model in modelica
Yang et al. Multi-objective optimal control approach for static voltage stability of power system considering interval uncertainty of the wind farm output
Ni et al. The study of energy filtering management process for microgrid based on the dynamic response model of vanadium redox flow battery
Behera et al. Voltage sag mitigation using distributed generation for an industrial distribution system
Aly et al. Voltage stability assessment for radial distribution power system with wind power penetration
CN110048428B (zh) 基于概率守恒原理的电力系统概率潮流计算方法
CN106410788A (zh) 基于传输介数的电网关键线路辨识方法
Ayan et al. An examination of the effects of distributed generation on distribution systems by load flow analysis
Al-Sakkaf et al. Power Flow Analysis of Weakly Meshed Distribution Network Including DG.
Saini et al. Digital Twin of a commercial building microgrid: Economic & environmental sustainability analysis
Zhao et al. Static sensitivity analysis method of integrated energy system based on power flow model
Oktaviani et al. Continuous power flow and time domain analysis for assessing voltage stability
CN112018787B (zh) 分布式电源的三相不平衡时序仿真方法
Zhang et al. Power flow analysis of integrated gas and electricity systems using the fast and flexible holomorphic embedding method
Lei et al. A Novel Methodology for Electric-Thermal Mixed Power Flow Simulation and Transmission Loss Analysis in Multi-Energy Micro-Grids
Jiang et al. Combined heat and power dispatch using simplified district heat flow model

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20211019