本发明公开了一种基于全纯函数的气电联合网络动态仿真计算方法，包括：输入气电联合网络的参数和运行数据、动态仿真的时段与步长；用时间变量的全纯函数代表气电联合网络动态方程中的未知变量，并计算气电联合网络动态方程中各全纯函数的未知系数；将当前时刻的时间变量代入全纯函数获取未知变量的取值，若当前时刻全纯函数满足计算精度要求，则继续用原全纯函数计算下一时刻的气电联合网络运行状态，直至整个仿真时段结束，反之则在下一时刻起重新计算全纯函数。本发明能够实现气电联合网络在运行中任意时刻的运行状态的准确计算，计算量少，在气电联合网络的运行分析与控制中具有较好的应用前景。
1.一种基于全纯函数的气电联合网络动态仿真计算方法，其特征在于，包括以下步骤：S1、输入气电联合网络的参数和运行数据、动态仿真的时段与步长；其中，设置动态仿end *真的时段为[0，t ]，步长用变量Δt代表，设置当前时刻t＝0；
S2、基于气电联合网络的参数和运行数据、动态仿真的时段与步长，用时间变量的全纯函数代表气电联合网络动态方程中的未知变量，构建含有全纯函数的气电联合网络动态方程，并计算在当前时刻各全纯函数的未知系数，包括以下步骤：S201、在电力系统中，用含有时间变量t的全纯函数V(t)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(t)写为以下形式：在天然气系统中，用含有时间变量t的全纯函数代表管道内天然气流量和节点处天然Pin气压强，具体包括四类全纯函数：代表天然气系统各管道入口处的天然气流量的F (t)、代Pout表天然气系统各管道出口处的天然气流量的F (t)、代表天然气系统各管道入口处的天Pin Pout然气压强的Π (t)、代表天然气系统各管道出口处的天然气压强的Π (t)，这四类全纯Pin函数写为幂级数形式，总共N阶，n代表阶数， 分别为幂级数F
(t)、F (t)、Π (t)、Π (t)第n阶项的系数，fn、Πn的幂级数表达式写为以下形式：气电联合网络其它的未知变量同样用含有时间变量t的全纯函数代表，具体包括四类GT全纯函数：电力系统燃气轮机所消耗天然气气量的F (t)、代表电力系统各燃气轮机发电量GT EC的P (t)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的P (t)、代表各天GC然气驱动的压缩机耗气量的F (t)，这四类全纯函数写为幂级数形式，总共N阶，n代表阶数，GT GT EC GC分别为幂级数F (t)、P (t)、P (t)、F (t)第n阶项的系数，各全纯函数的幂级数表达式写为以下形式：
式(1)至(9)中的全纯函数的零阶项系数全部已知，零阶项系数的值来自于气电联合网络在初始时刻各未知变量的值；
S202、利用式(1)至(9)中的全纯函数及时间变量t表达气电联合网络的动态方程组中的各项方程如下：式(10)至(12)中， 为代表电力系统PQ节点电压相量的全纯函数的共轭，计算公式为 Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(t)为代表电力系统所有节点电压相量的全纯函数； 代表电力系统PQ节点功率注入量的共轭值，代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(t)和 分别代表电力系统发电机节点电压相量的全纯函数及其共轭， 计算公式为Vigen,sp代表电力系统发电机节点电压相量预定值；Vipv(t)和 分别代表电力系统PV节点电压相量的全纯函数及其共轭， 计算公式为 Yipv,k和 分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，Pipv代表电力系统PV节点的有功注入功率， 代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
式(13)至式(17)中， 代表天然气系统中节点‑从该节点流出管道关联矩阵中对应节点ing、管道bng的元素， 代表天然气系统中节点‑流入该节点的管道关联矩阵中对应节点ing、管道bng的元素， 代表天然气系统管道bng入口处的天然气流量全纯函数、 代表天然气系统管道bng出口处的天然气流量全纯函数， 和 分别代表天然气系统中节点ing处用气量和供气量， 代表电力系统第igt台燃气轮机所消耗的天然气气量， 为天然气系统节点‑天然气驱动压缩机关联矩阵中对应天然气节点ing、天然气驱动压缩机igc的元素， 为天然气系统节点‑燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素； 分别代表天然气系统中配有压缩机的管道入口、出口处压强， 代表管道上压缩机出口、入口处压强比例； 代表以气源节点为出发点的管道的入口处压强， 为气源节点处预先设定的固定不变的压强值；
代表天然气系统未配有压缩机的管道bng1入口处、出口处压强，
代表管道bng1入口处、出口处流量，Abng1代表管道bng1的横截面面积，lbng1代表管道bng1的长度，c代表声波在天然气中的传播速度；λbng1代表管道bng1的摩擦系数，Dbng1代表管道bng1的直径；
式(18)至(20)中，aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数， 代表电力系统中第igt台燃气轮机的发电量； 代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数， 代表电驱动压缩机耗电量计算参数，代表天然气系统中管道iecg入口处天然气流量； 代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数， 代表天然气系统中管道igcb入口处天然气流量；
S203、将式(1)到(9)代入式(10)到(20)，获得含有全纯函数的非线性方程组；
S204、提取出步骤S203得到的非线性方程组中所有的t项，获得仅包括全纯函数一阶项Pin Pout Pin Pout系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、Π (t)、GT GT EC GCF (t)、P (t)、P (t)、F (t)的一阶项系数；
Pin Pout Pin Pout GT GT EC GC
S205、将已获得的V(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F2
(t)的一阶项系数代入步骤S203得到的非线性方程组，提取出非线性方程组中所有的t 项，Pin Pout Pin Pout GT GT EC GC获得仅包括V(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)二阶项Pin Pout Pin Pout系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、Π (t)、GT GT EC GCF (t)、P (t)、P (t)、F (t)的二阶项系数；
Pin Pout Pin Pout GT GT EC GC
S206、将V(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)中的已Pin计算出的系数代入步骤S203得到的非线性方程组，提取更高阶的t项，获取仅包括V(t)、FPout Pin Pout GT GT EC GC(t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)中该阶数系数的线性方程组进行求解，直至完成第N阶项系数的计算；
S3、将当前时刻的时间变量代入气电联合网络动态方程中的全纯函数，获取未知变量在当前时刻的取值，将未知变量的取值代入气电联合网络动态方程，判断全纯函数是否满足计算的精度要求；
S4、若当前时刻全纯函数满足计算精度要求，则继续用原全纯函数计算下一时刻的气电联合网络运行状态，直至整个仿真时段结束，若当前时刻全纯函数不满足计算精度要求，则在下一时刻起重新计算全纯函数，即重复步骤S2至步骤S4。
2.根据权利要求1所述的一种基于全纯函数的气电联合网络动态仿真计算方法，其特征在于：在步骤S1中，所述气电联合网络的参数包括：电力系统的拓扑结构、线路导纳参数；
天然气系统网络拓扑结构、管道压降系数；耦合元件的各项参数，即电驱动压缩机耗电量和天然气驱动压缩机耗气量计算公式中各项参数、燃气轮机气电转换效率；
所述气电联合网络的运行数据包括：电力系统中各节点的负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角；天然气系统中各节点的负荷、各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比；
气电联合网络的未知变量包括：电力系统中各节点的电压向量、电力系统中各燃气轮机所消耗的天然气气量、电力系统各燃气轮机发电量、电力系统为天然气系统中各电驱动的压缩机所提供电量、各天然气驱动的压缩机耗气量、天然气系统各管道出入口处的天然气流量、天然气系统各节点的天然气压强。
3.根据权利要求1所述的一种基于全纯函数的气电联合网络动态仿真计算方法，其特征在于：所述步骤S3包括以下步骤：* *
S301、采用变量t代表当前时刻，将t代入步骤S2中的式(1)到(9)，获取当前时刻各未知变量的值，并将未知变量取值代入步骤S3中的式(10)到(20)，计算式(10)到(20)左右两侧的差值；
S302、若式(10)到(20)左右两侧的差值全部低于预设标准，则认为全纯函数满足计算的精度要求，反之则认为全纯函数不满足计算的精度要求；
S401、如果t≥t 时，动态仿真完成，输出计算数据；反之则进入步骤S402；
S402、全纯函数满足计算的精度要求，则直接进入下一仿真时刻气电联合网络运行状*态的计算，将下一仿真时刻对应的时间变量t+Δt代入式(10)到(20)，重复步骤S3；如果全* *纯函数未满足计算的精度要求，则设置t ＝t+Δt，返回步骤S2重新构造全纯函数、含有全纯函数的气电联合网络动态方程，重复步骤S2至S4。
一种基于全纯函数的气电联合网络动态仿真计算方法
技术领域
[0001] 本发明涉及气电联合网络动态仿真的技术领域，尤其是指一种基于全纯函数的气电联合网络动态仿真计算方法。
背景技术
[0002] 近年来，传统的发电机组逐渐由可再生能源、清洁能源驱动的发电机组代替，为了保证系统运行的可靠性，系统接入越来越多的天然气发电机组以实现对于突变的负荷、发电量的快速响应。在实际运行中，供能系统需要同时考虑电力、天然气的供应以及两者之间的关联影响，通常电力系统、天然气系统被视作为一个整体(即气电联合网络)进行建模、计算、分析等。天然气系统中，当管道一侧的压强或流量改变时，管道内部以及另一侧的压强、流量也会发生波动；在变化结束后，系统需要一定时间达到新的稳态。为了更加准确地描述气电联合网络在系统运行条件发生变化后的波动过程，常采用动态仿真计算计算气电联合网络的运行状态，通过计算气电联合网络的状态变量的值刻画出系统的波动，更加准确地描述系统从受到扰动到进入新稳态的动态过程。
[0003] 气电联合网络的动态仿真计算中，通常采用偏微分方程组对气电联合网络进行建模，使用有限差分法、有限元法等数值解法实现求解。这些数值解法存在以下问题：求解过程中计算量大小取决于时间、空间的差分大小，通常计算耗时较长；动态仿真计算的结果为一段仿真时段内多个特定的时刻的数据，各时刻之间的间隔在仿真计算前已确定，若需要计算其他时刻的数据则需要使用插值法或重新计算，导致整体计算效率较低。需要引入新的方法实现准确、快速的动态仿真。
发明内容
[0004] 本发明的目的在于克服现有技术的缺点与不足，提出了一种基于全纯函数的气电联合网络动态仿真计算方法，使用时间变量的全纯函数代替气电联合网络的状态变量并进行气电联合网络动态方程的求解，从而获取气电联合网络运行过程中各时刻的运行状态。
[0005] 为实现上述目的，本发明所提供的技术方案为：一种基于全纯函数的气电联合网络动态仿真计算方法，包括以下步骤：
[0006] S1、输入气电联合网络的参数和运行数据、动态仿真的时段与步长；其中，设置动end *态仿真的时段为[0，t ]，步长用变量Δt代表，设置当前时刻t＝0
[0007] S2、基于气电联合网络的参数和运行数据、动态仿真的时段与步长，用时间变量的全纯函数代表气电联合网络动态方程中的未知变量，构建含有全纯函数的气电联合网络动态方程，并计算在当前时刻各全纯函数的未知系数；
[0008] S3、将当前时刻的时间变量代入气电联合网络动态方程中的全纯函数，获取未知变量在当前时刻的取值，将未知变量的取值代入气电联合网络动态方程，判断全纯函数是否满足计算的精度要求；
[0009] S4、若当前时刻全纯函数满足计算精度要求，则继续用原全纯函数计算下一时刻的气电联合网络运行状态，直至整个仿真时段结束，若当前时刻全纯函数不满足计算精度要求，则在下一时刻起重新计算全纯函数，即重复步骤S2至步骤S4。
[0010] 进一步，在步骤S1中，所述气电联合网络的参数包括：电力系统的拓扑结构、线路导纳参数；天然气系统网络拓扑结构、管道压降系数；耦合元件的各项参数，即电驱动压缩机耗电量和天然气驱动压缩机耗气量计算公式中各项参数、燃气轮机气电转换效率；
[0011] 所述气电联合网络的运行数据包括：电力系统中各节点的负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角；天然气系统中各节点的负荷、各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比；
[0012] 气电联合网络的未知变量包括：电力系统中各节点的电压向量、电力系统中各燃气轮机所消耗的天然气气量、电力系统各燃气轮机发电量、电力系统为天然气系统中各电驱动的压缩机所提供电量、各天然气驱动的压缩机耗气量、天然气系统各管道出入口处的天然气流量、天然气系统各节点的天然气压强。
[0013] 进一步，所述步骤S2包括以下步骤：
[0014] S201、在电力系统中，用含有时间变量t的全纯函数V(t)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(t)写为以下形式：
[0015]
[0016] 在天然气系统中，用含有时间变量t的全纯函数代表管道内天然气流量和节点处Pin天然气压强，具体包括四类全纯函数：代表天然气系统各管道入口处的天然气流量的FPout
(t)、代表天然气系统各管道出口处的天然气流量的F (t)、代表天然气系统各管道入口Pin Pout
处的天然气压强的Π (t)、代表天然气系统各管道出口处的天然气压强的Π (t)，这四类全纯函数写为幂级数形式，总共N阶，n代表阶数， 分别为幂级
Pin Pout Pin Pout
数F (t)、F (t)、Π (t)、Π (t)第n阶项的系数，fn、Πn的幂级数表达式写为以下形式：
[0017]
[0018]
[0019]
[0020]
[0021] 气电联合网络其它的未知变量同样用含有时间变量t的全纯函数代表，具体包括GT四类全纯函数：电力系统燃气轮机所消耗天然气气量的F (t)、代表电力系统各燃气轮机发GT EC
电量的P (t)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的P (t)、代表GC
各天然气驱动的压缩机耗气量的F (t)，这四类全纯函数写为幂级数形式，总共N阶，n代表GT GT EC GC
阶数， 分别为幂级数F (t)、P (t)、P (t)、F (t)第n阶项的系数，
各全纯函数的幂级数表达式写为以下形式：
[0022]
[0023]
[0024]
[0025]
[0026] 设置整个仿真时段的初始时刻为t0，t0＝0；
[0027] 式(1)至(9)中的全纯函数的零阶项系数全部已知，零阶项系数的值来自于气电联合网络在初始时刻各未知变量的值；
[0028] S202、利用式(1)至(9)中的全纯函数及时间变量t表达气电联合网络的动态方程组中的各项方程如下：
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
[0040] 式(10)至(12)中， 为代表电力系统PQ节点电压相量的全纯函数的共轭，计算公式为 Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(t)为代表电力系统所有节点电压相量的全纯函数； 代表电力系统PQ节点功率注入量的共轭值， 代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(t)和分别代表电力系统发电机节点电压相量的全纯函数及其共轭， 计算公式为
Vigen,sp代表电力系统发电机节点电压相量预定值；Vipv(t)和 分别
代表电力系统PV节点电压相量的全纯函数及其共轭， 计算公式为
Yipv,k和 分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，Pipv代表电力系统PV节点的有功注入功率， 代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
[0041] 式(13)至式(17)中， 代表天然气系统中节点‑从该节点流出管道关联矩阵中对应节点ing、管道bng的元素， 代表天然气系统中节点‑流入该节点的管道关联矩阵Pin中对应节点ing、管道bng的元素，F (t)代表天然气系统管道bng入口处的天然气流量全纯Pout
函数、F (t)代表天然气系统管道bng出口处的天然气流量全纯函数， 和 分别代表天然气系统中节点ing处用气量和供气量， 代表电力系统第igt台燃气轮机所消耗的天然气气量， 代表天然气系统中节点ing处向天然气驱动的压缩机耗气量的全纯函数， 为天然气系统节点‑天然气驱动压缩机关联矩阵中对应天然气节点ing、天然气驱动压缩机igc的元素， 为天然气系统节点‑燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素； 分别代表天然气系统中配有压缩机的管道入口、
出口处压强， 代表管道上压缩机出口、入口处压强比例； 代表以气源节点为出
发点的管道的入口处压强， 为气源节点处预先设定的固定不变的压强值；
代表天然气系统未配有压缩机的管道bng1入口处、出口处压强，
代表管道bng1入口处、出口处流量，Abng1代表管道bng1的横截面面积，
lbng1代表管道bng1的长度，c代表声波在天然气中的传播速度；λbng1代表管道bng1的摩擦系数，Dbng1代表管道bng1的直径；
[0042] 式(18)至(20)中，aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数，代表电力系统中第igt台燃气轮机的发电量； 代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数， 代表电驱动压缩机耗电量计算参数， 代表天然气系统中管道iecg入口处天然气流量； 代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数， 代表天然气系统中管道igcb入口处天然气流量；
[0043] S203、将式(1)到(9)代入式(10)到(20)，获得含有全纯函数的非线性方程组；
[0044] S204、提取出步骤S203得到的非线性方程组中所有的t项，获得仅包括全纯函数一Pin Pout Pin Pout阶项系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、ΠGT GT EC GC
(t)、F (t)、P (t)、P (t)、F (t)的一阶项系数；
[0045] S205、将已获得的V(t)、FPin(t)、FPout(t)、ΠPin(t)、ΠPout(t)、FGT(t)、PGT(t)、PECGC(t)、F (t)的一阶项系数代入步骤S203得到的非线性方程组，提取出非线性方程组中所有
2 Pin Pout Pin Pout GT GT EC GC
的t项，获得仅包括V(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)Pin Pout Pin
二阶项系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、Pout GT GT EC GC
Π (t)、F (t)、P (t)、P (t)、F (t)的二阶项系数；
[0046] S206、将V(t)、FPin(t)、FPout(t)、ΠPin(t)、ΠPout(t)、FGT(t)、PGT(t)、PEC(t)、FGC(t)中的已计算出的系数代入步骤S203得到的非线性方程组，提取更高阶的t项，获取仅包括VPin Pout Pin Pout GT GT EC GC(t)、F (t)、F (t)、П (t)、Π (t)、F (t)、P (t)、P (t)、F (t)中该阶数系数的线性方程组进行求解，直至完成第N阶项系数的计算；
[0047] 所述步骤S3包括以下步骤：
[0048] S301、采用变量t*代表当前时刻，将t*代入步骤S3中的式(1)到(9)，获取当前时刻各未知变量的值，并将未知变量取值代入步骤S3中的式(10)到(20)，计算式(10)到(20)左右两侧的差值；
[0049] S302、若式(10)到(20)左右两侧差值全部低于预设标准，则认为全纯函数满足计算的精度要求，反之则认为全纯函数不满足计算的精度要求；
[0050] 所述步骤S4包括以下步骤：
[0051] S401、如果t*≥tend时，动态仿真完成，输出计算数据；反之则进入步骤S402；
[0052] S402、全纯函数满足计算的精度要求，则直接进入下一仿真时刻气电联合网络运*行状态的计算，将下一仿真时刻对应的时间变量t+Δt代入式(10)到(20)，重复步骤S3；如* *
果全纯函数未满足计算的精度要求，则设置t＝t+Δt，返回步骤S2重新构造全纯函数、含有全纯函数的气电联合网络动态方程，重复步骤S2至S4。
[0053] 本发明与现有技术相比，具有如下优点与有益效果：
[0054] 1、本发明方法使用时间变量的全纯函数代替气电联合网络的未知变量，可以实现气电联合网络的运行状态求解，在气电联合网络运行状态随时可能发生变化的条件下准确获取电压、天然气压强、天然气流量的值，为系统的正常、稳定、安全、经济运行提供了详细全面的运行数据。
[0055] 2、本发明方法最重要特性即为气电联合网络的各未知变量由时间的全纯函数表示，通过代入时间变量便能直接获取各未知变量的值。和常用的有限差分法相比，本发明方法无需在每个仿真时刻进行全纯函数的计算，本发明中每间隔一个仿真步长边进行一次气电联合网络动态方程组的不平衡量计算，当全纯函数不满足计算的精度要求时才重新进行全纯函数的计算，显著地减少了动态仿真的计算量，提高了计算速度。
[0056] 3、和常用的有限差分法(仿真计算仅在离散时刻进行)相比，本发明可以实现面向气电联合网络的运行状态的连续动态仿真。对于未获取仿真数据的时刻，常用方法需要采用插值法进行估算，或重新设置仿真步长再次进行仿真计算，本发明提出的方法则可以代入时间直接计算该时刻的未知变量值，明显地提高了计算效率，保证了计算的准确性。
[0057] 4、本发明方法能够实现气电联合网络在日常运行中的动态仿真，准确地反映气电联合网络运行条件发生变化时各变量的波动过程，为气电联合网络的运行优化、计算分析提供了详细全面的数据，同时该方法通用性强，可以用于含有其他电力系统、天然气系统及耦合部分的元件的气电联合网络的动态仿真。
附图说明
[0058] 图1为本发明方法逻辑流程示意图。
[0059] 图2为气电联合网络的架构图。
具体实施方式
[0060] 下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
[0061] 如图1所示，本实施例所提供的基于全纯函数的气电联合网络动态仿真计算方法，该方法使用时间变量的全纯函数代替气电联合网络的状态变量并进行气电联合网络动态方程的求解，从而获取气电联合网络运行过程中各时刻的运行状态，其具体情况如下：
[0062] S1、输入气电联合网络的参数和运行数据、动态仿真的时段与步长。
[0063] 所述气电联合网络的参数包括：电力系统的拓扑结构、线路导纳参数；天然气系统网络拓扑结构、管道压降系数；耦合元件的各项参数，即电驱动压缩机耗电量和天然气驱动压缩机耗气量计算公式中各项参数、燃气轮机气电转换效率。
[0064] 所述气电联合网络的运行数据包括：电力系统中各节点的负荷、除松弛节点以外的各发电机组出力与电压幅值、松弛节点的电压幅值与相角；天然气系统中各节点的负荷、各节点处的天然气负荷、除松弛节点以外各气源的供气量、松弛节点的压强、压缩机进出口压强比。
[0065] 气电联合网络的未知变量包括：电力系统中各节点的电压向量、电力系统中各燃气轮机所消耗的天然气气量、电力系统各燃气轮机发电量、电力系统为天然气系统中各电驱动的压缩机所提供电量、各天然气驱动的压缩机耗气量、天然气系统各管道出入口处的天然气流量、天然气系统各节点的天然气压强。
[0066] 动态仿真的时段为[0，tend](tend＝100s)，步长用变量Δt代表(Δt＝1s)，设置当*前时刻t＝0。
[0067] 如图2所示，天然气系统含有7个节点、6个管道；电力系统含有14个节点、其中5个节点为PV节点；天然气系统和电力系统通过燃气轮机(连接天然气系统的节点3和电力系统的节点2)、天然气驱动的压缩机(连接天然气系统的压缩机和电力系统的节点5)互相耦合。
[0068] S2、基于气电联合网络的参数和运行数据、动态仿真的时段与步长，用时间变量的全纯函数代表气电联合网络动态方程中的未知变量，构建含有全纯函数的气电联合网络动态方程，并计算在当前时刻各全纯函数的未知系数，具体步骤如下：
[0069] S201、在电力系统中，用含有时间变量t的全纯函数V(t)代表各节点的电压向量，全纯函数写为总共N阶的幂级数形式，n代表阶数，vn为幂级数第n阶项的系数，V(t)写为以下形式：
[0070]
[0071] 在天然气系统中，用含有时间变量t的全纯函数代表管道内天然气流量和节点处Pin天然气压强，具体包括四类全纯函数：代表天然气系统各管道入口处的天然气流量的FPout
(t)、代表天然气系统各管道出口处的天然气流量的F (t)、代表天然气系统各管道入口Pin Pout
处的天然气压强的Π (t)、代表天然气系统各管道出口处的天然气压强的Π (t)，这四类全纯函数写为幂级数形式，总共N阶，n代表阶数， 分别为幂级
Pin Pout Pin Pout
数F (t)、F (t)、Π (t)、Π (t)第n阶项的系数，fn、Πn的幂级数表达式写为以下形式：
[0072]
[0073]
[0074]
[0075]
[0076] 气电联合网络其它的未知变量同样用含有时间变量t的全纯函数代表，具体包括GT四类全纯函数：电力系统燃气轮机所消耗天然气气量的F (t)、代表电力系统各燃气轮机发GT EC
电量的P (t)、代表电力系统为天然气系统中各电驱动的压缩机所提供电量的P (t)、代表GC
各天然气驱动的压缩机耗气量的F (t)，这四类全纯函数写为幂级数形式，总共N阶，n代表GT GT EC GC
阶数， 分别为幂级数F (t)、P (t)、P (t)、F (t)第n阶项的系数，
各全纯函数的幂级数表达式写为以下形式：
[0077]
[0078]
[0079]
[0080]
[0081] 设置整个仿真时段的初始时刻为t0，t0＝0。
[0082] 式(1)至(9)中的全纯函数的零阶项系数全部已知，零阶项系数的值来自于气电联合网络在t0时各未知变量的值。
[0083] S202、利用式(1)至(9)中的全纯函数及时间变量t表达气电联合网络的动态方程组中的各项方程如下：
[0084]
[0085]
[0086]
[0087]
[0088]
[0089]
[0090]
[0091]
[0092]
[0093]
[0094]
[0095] 式(10)至(12)中， 为代表电力系统PQ节点电压相量的全纯函数的共轭，计算公式为 Yipq,k代表电力系统的导纳矩阵中PQ节点对应所有节点的元素，Vk(t)为代表电力系统所有节点电压相量的全纯函数； 代表电力系统PQ节点功率注入量的共轭值， 代表电力系统PQ节点处电驱动压缩机耗电量的全纯函数；Vigen(t)和
分别代表电力系统发电机节点电压相量的全纯函数及其共轭， 计算公式为
Vigen,sp代表电力系统发电机节点电压相量预定值；Vipv(t)和 分别
代表电力系统PV节点电压相量的全纯函数及其共轭， 计算公式为
Yipv,k和 分别代表电力系统的导纳矩阵元素中PV节点对应所有节点的元素及其共轭，Pipv代表电力系统PV节点的有功注入功率， 代表电力系统PV节点处为天然气系统中电驱动压缩机所提供电量的全纯函数；
[0096] 式(13)至式(17)中， 代表天然气系统中节点‑从该节点流出管道关联矩阵中对应节点ing、管道bng的元素， 代表天然气系统中节点‑流入该节点的管道关联矩阵Pin中对应节点ing、管道bng的元素，F (t)代表天然气系统管道bng入口处的天然气流量全纯Pout
函数、F (t)代表天然气系统管道bng出口处的天然气流量全纯函数， 和 分别代表天然气系统中节点ing处用气量和供气量， 代表电力系统第igt台燃气轮机所消耗的天然气气量， 代表天然气系统中节点ing处向天然气驱动的压缩机耗气量的全纯函数，为天然气系统节点‑天然气驱动压缩机关联矩阵中对应天然气节点ing、天然气驱动压缩机igc的元素， 为天然气系统节点‑燃气轮机关联矩阵中对应天然气节点ing、燃气轮机igt的元素； 分别代表天然气系统中配有压缩机的管道入口、出
口处压强， 代表管道上压缩机出口、入口处压强比例； 代表以气源节点为出发
点的管道的入口处压强， 为气源节点处预先设定的固定不变的压强值；
代表天然气系统未配有压缩机的管道bng1入口处、出口处压强，
代表管道bng1入口处、出口处流量，Abng1代表管道bng1的横截面面积，
lbng1代表管道bng1的长度，c代表声波在天然气中的传播速度；λbng1代表管道bng1的摩擦系数，Dbng1代表管道bng1的直径。
[0097] 式(18)至(20)中，aigt、bigt、cigt代表燃气轮机耗气量计算公式中各项参数，代表电力系统中第igt台燃气轮机的发电量； 代表电力系统节点iec处向天然气系统管道iecg上的电驱动压缩机所提供电量的全纯函数， 代表电驱动压缩机耗电量计算参数， 代表天然气系统中管道iecg入口处天然气流量； 代表天然气系统节点igc处向天然气系统管道igcb上的天然气驱动压缩机耗气量的全纯函数，αigc、βigc、γigc代表天然气驱动压缩机耗气量计算公式中的各项参数， 代表天然气系统中管道igcb入口处天然气流量。
[0098] S203、将式(1)到(9)代入式(10)到(20)，获得含有全纯函数的非线性方程组；
[0099] S204、提取出步骤S301得到的非线性方程组中所有的t项，获得仅包括全纯函数一Pin Pout Pin Pout阶项系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、ΠGT GT EC GC
(t)、F (t)、P (t)、P (t)、F (t)的一阶项系数；
[0100] S205、将已获得的V(t)、FPin(t)、FPout(t)、ΠPin(t)、ΠPout(t)、FGT(t)、PGT(t)、PECGC(t)、F (t)的一阶项系数代入步骤S203得到的非线性方程组，提取出非线性方程组中所有
2 Pin Pout Pin Pout GT GT EC GC
的t项，获得仅包括V(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)Pin Pout Pin
二阶项系数的线性方程组，通过求解该线性方程组获得V(t)、F (t)、F (t)、Π (t)、Pout GT GT EC GC
Π (t)、F (t)、P (t)、P (t)、F (t)的二阶项系数；
[0101] S206、将V(t)、FPin(t)、FPout(t)、ΠPin(t)、ΠPout(t)、FGT(t)、PGT(t)、PEC(t)、FGC(t)中的已经计算出的系数代入步骤S203得到的非线性方程组，提取更高阶的t项，获取仅包括Pin Pout Pin Pout GT GT EC GCV(t)、F (t)、F (t)、Π (t)、Π (t)、F (t)、P (t)、P (t)、F (t)中该阶数系数的线性方程组进行求解，直至完成第N阶项系数的计算。
[0102] S3、将当前时刻的时间变量代入气电联合网络动态方程中的全纯函数，获取未知变量在当前时刻的取值，将未知变量的取值代入气电联合网络动态方程，判断全纯函数是否满足计算的精度要求，具体步骤如下：
[0103] S301、采用变量t*代表当前时刻，将t*代入式(1)到(9)，获取当前时刻各未知变量的值，并将未知变量取值代入式(10)到(20)，计算式(10)到(20)左右两侧的差值。
[0104] S302、若式(10)到(20)左右两侧不平衡量全部低于预设标准，则认为全纯函数满足计算的精度要求，反之则认为全纯函数不满足计算的精度要求。
[0105] S4、若当前时刻全纯函数满足计算精度要求，则继续用原全纯函数计算下一时刻的气电联合网络运行状态，直至整个仿真时段结束，若当前时刻全纯函数不满足计算精度要求，则在下一时刻起重新计算全纯函数，即重复步骤S2至步骤S4。具体步骤如下：
[0106] S401、如果t*≥tend时，动态仿真完成，输出计算数据；反之则进入步骤S402；
[0107] S402、全纯函数满足计算的精度要求，则直接进入下一仿真时刻气电联合网络运*行状态的计算，将下一仿真时刻对应的时间变量t+Δt代入式(10)到(20)，重复步骤S3；如* *
果全纯函数未满足计算的精度要求，则设置t＝t+Δt，返回步骤S2重新构造全纯函数、含有全纯函数的气电联合网络动态方程，重复步骤S2至S4。
[0108] 上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。
