本发明属于电力系统和天然气系统协同运行领域，涉及一种电力系统和天然气系统协同运行的建模方法。
背景技术：
为实现经济社会的可持续发展，大力推进清洁能源和可再生能源的使用成为必然趋势。天然气作为一种清洁能源，有很大的使用前景，“以气代煤”成为能源结构调整的大趋势。2015年，我国的天然气用量大约1900亿立方米，天然气消费比重仅占5％左右，相较之下，全球水平为21％，美国则达到28％，大力发展天然气已经成为我国新能源发展的大势所趋。国务院办公厅在能源发展战略中指出，到2020年，天然气在一次能源消耗中占比要提升至10％，根据iea在2012年的预测结果，我国天然气年产量将于2025年突破2000×108m3，达到2170×108m3；2035年将突破3000×108m3，达到3180×108m3，未来天然气有望成为煤炭、石油后的第三大主体能源。
在传统的社会供能系统中，电力供能系统和天然气供能系统一般是独立运行的，各个能源系统按照自身的运行状况做出调度策略、调整峰谷的规划，这使得能源和设备的利用效率下降。随着各种能源之间联系越来越密切，电力系统和天然气系统之间的耦合也越来越强。电力系统和天然气系统的实时协同运行状态分析可以作为两个系统在优化调度、市场决策和安全保护等方面的研究基础，是能源使用情况的重要参考。
在电力系统分析的角度，研究电力系统和天然气系统协同运行时，需要对天然气系统进行合理有效的建模。已有的天然气建模存在以下问题：
(1)对实际的管道特性考虑不足，一般是以一个人为给定的数值代替管道性质对流过天然气流量的影响，这缺失实际的物理含义而且和实际运行情况具有较大的误差，是一种理想化的建模。
(2)对天然气系统中的压缩机建模过于复杂，且不利于同电力系统统一求解，其非线性方程组在协同运行求解时由于线性化，导致误差过大，且计算量大，计算用时长。
发明目的
本发明的目的在于针对电力系统和天然气系统协同运行的建模方法偏离实际的问题，同时针对现有的压缩机建模方法的不足，提出一种电力系统和天然气系统协同运行的建模方法，对天然气管道特性的建模更加贴近实际，从实际管道的物理特性抽象得到；对天然气系统中的压缩机建模提出了实用性处理方法，使得天然气系统和电力系统的建模更加统一，求解过程中计算量更小，计算速度更快。
为实现上述目的，本发明采用的技术方案是：
一种电力系统和天然气系统协同运行的建模方法，包括以下步骤：
步骤i：等值天然气管网参数：
天然气管网的拓扑结构以及管道中流过天然气的流动方向用一个矩阵表征，称之为标识函数矩阵；天然气管道对天然气流量的影响特性等值为另一个矩阵，称之为管道导纳矩阵；
步骤ii：等值电网输电线路参数：
将电力系统的输电线路和变压器等值为节点导纳矩阵yb；
步骤iii：天然气系统中处理压缩机的实用化方法：
根据压缩机的耗能类型将压缩机分为燃气压缩机和电动压缩机，引入压缩机耗能的度量，并将压缩机从天然气管网中去除，改变管网结构，相应的站点等效天然气气负荷和天然气气源；
步骤iv：建立天然气管网和电网运行的数学模型：
建立天然气管道中流过天然气流量和电网线路中流过电功率的数学模型表达。
本发明进一步的改进在于，步骤i中，在天然气系统中，为了表征管道内天然气的流动方向，引入一个标识函数amn，其函数值与天然气站点m和n的站点压力相对大小有关；若两个站点没有管道相连，则amn＝0；若两个站点有管道相连，相对于站点m而言，如果站点m的站点压力pm比与其相连站点n的站点压力pn大，那么amn＝+1，反之amn＝-1；
每条天然气管道对应的一个标识函数amn，能够生成标识函数矩阵amn，amn是一个方阵，其阶数等于天然气系统的管道数量。
本发明进一步的改进在于，步骤i中，天然气输配管道存在管道阻力，用rmn来表征连接天然气站点m和站点n的管道对流过天然气的阻力大小，称之为管道阻抗；管道阻抗rmn随天然气管道压力等级的不同，存在数量级上的差异，低压管道阻抗表达式如式(2)所示，中高压管道阻抗表达式如(3)所示：
在式(2)和式(3)中，d为天然气管道内径；l为天然气管道长度；ρ0为天然气密度；t为天然气温度；t0为标准状态绝对温度，λ为天然气管道内壁摩擦阻力系数。
本发明进一步的改进在于，步骤i中，天然气管道为地埋式天然气管道。
本发明进一步的改进在于，步骤i中，相连的两个天然气站点之间的管道阻抗能够用rmn来表征，不相连的两个天然气站点不存在管道阻抗，为简化天然气管网的描述，用管道导纳cmn表示，定义管道导纳cmn与管道阻抗rmn的大小关系如式(4)所示：
每条天然气管道对应的一个管道导纳cmn，能够生成管道导纳矩阵cmn，cmn是一个方阵，其阶数等于天然气系统的管道数量。
本发明进一步的改进在于，步骤iii中，在天然气系统中，随天然气输配距离的增大，输气管道内的压力下降；因此，在天然气系统中要增设压缩机，用来增加天然气的压力，来补偿燃气输配过程中的压力损耗，压缩机属于天然气系统的非管道支路；
将含压缩机的支路去掉，改变天然气管网的拓扑结构，等效得到新的天然气管网拓扑结构；在去掉压缩机后，将压缩机的出口站点处，等效注入压缩机的输出天然气流量，作为该站点的天然气源注入；将压缩机的入口站点处，等效施加压缩机的流入天然气流量，作为该站点的天然气负荷。
本发明进一步的改进在于，将天然气系统中的压缩机分为两种：一种是燃气压缩机；另一种是电动压缩机；
燃气压缩机升压消耗的能量egas与流过燃气压缩机的天然气流量和升压大小的关系如式(5)所示：
式(5)中，egas为燃气压缩机在单位时间内消耗能量；ηgas为燃气压缩机的效率；λg为天然气的比热容；t0为标准状态绝对温度；tg为天然气温度；p0为天然气的标准压力；zg为天然气的压缩系数；fmn为流出燃气压缩机的流量；pin和pout分别为燃气压缩机入口压力和出口压力；
所需的燃气压缩机在单位时间内消耗能量egas由拖动燃气压缩机的燃气轮机消耗天然气供给，燃气轮机的耗量特性方程为式(6)：
式(6)中，α、β和γ为燃气轮机的耗量系数；fcom为燃气轮机用于带动燃气压缩机消耗的天然气流量；
将燃气压缩机去掉后，压缩机入口站点处添加等效天然气负荷fload＝fin，其中，fin是流入压缩机的天然气流量，压缩机出口站点处添加等效的天然气气源fsource＝fout，其中，fout是流出压缩机的天然气流量，fload＞fsource，差值fcom＝fload-fsource用于供升压耗能；
电动压缩机单位时间消耗能量eee与流过电动压缩机的流量和升压大小的关系，如式(7)所示：
式(7)中，eee为电动压缩机单位时间消耗能量；ηee为电动压缩机的效率；
所需的电动压缩机单位时间消耗能量eee由电动机供给，电动机消耗的电功率pcom和电动压缩机在单位时间内消耗能量eee的关系如式(8)所示：
由于电动压缩机不消耗天然气流量，所以流入电动压缩机的流量和流出压缩机的流量大小是相等的，即fin＝fout；在电动压缩机入口站点等效的气负荷fload和出口站点等效的气源fsource也是相等的，即fload＝fsource。
本发明进一步的改进在于，对于天然气系统：
天然气管网中流过的天然气流量需要通过水力分析，由水力方程推导得出天然气管网的运行方程，流体动力学方程、连续性方程和流体状态方程分别如式(9)～(11)所示：
流体动力学方程：
连续性方程：
流体状态方程：
p＝zρrt(11)
其中：ω为天然气流速；ρ为天然气密度；p为天然气压力；t为天然气温度；t为时间；x为管道轴向坐标；g为重力加速度；α为管道与水平面夹角；λ为管道内壁摩擦系数；d为管道内径；z为压缩系数；r为气体常数。
稳定流量基本方程如式(12)所示。
式(12)中，f表示天然气流量，下标“0”表示基准状态：基准压力p0＝101325pa，基准温度t0＝293.15k，基准压缩系数z0是p0和t0的函数；
由于不考虑管道敷设高程变化，所以认为管道与水平面夹角α＝0°；
通过两边取积分，并用管道导纳表示实际管道情况对天然气的影响，得到天然气系统中的管道支路的流量方程，如式(13)所示。
其中，fmn为天然气系统中的管道支路的流量，hmn为管道压力降，它表征两个站点之间绝对压力的差值，管道压力降的表达式如式(14)所示。
其中，各个站点的压力值和各个管道中的流量值均有上下限的要求，如式(15)和式(16)所示。
pimin≤pi≤pimax(15)
fmnmin≤fmn≤fmnmax(16)
式中，pi为任意一个站点的压力值，i取1,2,…ngas，ngas为天然气站点的总数，pimin为第i个站点压力值的下限值，pimax为第i个站点压力值的上限值；
fmn为站点m和站点n之间连接管道中流过的天然气流量，fmnmin为站点m和站点n之间连接管道中流过的天然气流量的下限值，fmnmax为站点m和站点n之间连接管道中流过的天然气流量上限值；
对于电力系统：
建立电网潮流方程的理论基础是电路理论中的节点分析法，节点电压用极坐标形式表示，如式(17)和式(18)所示。
pij＝uiuj(gijcosδij+bijsinδij)(17)
qij＝uiuj(gijsinδij-bijcosδij)(18)
式中，pij为电力系统中母线i和母线j之间连接上线路流过的有功功率，qij为电力系统中母线i和母线j之间连接上线路流过的无功功率ui为电力系统母线i的电压幅值，uj为电力系统母线j的电压幅值，gij为电力系统中母线i和母线j之间连接上线路的电导，bij为电力系统中母线i和母线j之间连接上线路的电纳，δij为电力系统母线i和母线j之间的电压相角的差值，gij和bij都是节点导纳矩阵yb中的元素；
其中，各个母线的电压值和各个线路上的功率均有上下限的要求，如式(19)～(21)所示。
uimin≤ui≤uimax(19)
pijmin≤pij≤pijmax(20)
qijmin≤qij≤qijmax(21)
其中，ui为电力系统中任意母线的电压幅值，i为1,2,…nee，nee为电力系统母线的总数，uimin为电力系统中母线i的电压幅值的下限值，uimax为电力系统中母线i的电压幅值的上限值，pijmin为电力系统中母线i和母线j之间连接上线路流过的有功功率的下限值，pijmax为电力系统中母线i和母线j之间连接上线路流过的有功功率的上限值，qijmin为电力系统中母线i和母线j之间连接上线路流过的无功功率的下限值，qijmax为电力系统中母线i和母线j之间连接上线路流过的无功功率的上限值。
与现有技术相比，本发明具有的有益效果：本发明采用的天然气管道特性建模方法，考虑了实际的天然气管道物理特性，使得管道对其中流过的天然气流量的影响更加切合实际情况，天然气系统运行状态的分析更加准确有效；此外，本发明提出的处理天然气系统中压缩机的实用化方法，通过改变天然气管网的拓扑结构，有效地化简了天然气系统中站点之间的连接类型，使得天然气系统中的站点均通过管道支路连接，非管道支路从管网的拓扑结构中去除，等效方式合理准确，减小了压缩机复杂的非线性方程线性化后，对天然气系统运行求解的工程误差，避免了求解时，计算量大和计算时间长的问题。
附图说明
图1为协同运行时电网用电量高低时刻对应天然气站点压力对比图。
具体实施方式
下面结合附图和具体算例对本发明作进一步详细说明。
采取分步建模，统一求解的思想，首先对天然气管网和电网输电线路进行等值参数的建模，重点在于对天然气管网建模内容的改进；在此基础上，对天然气系统和电力系统的耦合单元进行建模，重点在于对压缩机处理方法的改进；然后建立天然气管网和电网运行的数学模型，并进行协同运行的计算，验证建模方法的有效性和实用性。
具体包括以下步骤：
步骤i：等值天然气管网参数：
天然气管网的拓扑结构以及管道中流过天然气的流动方向可以用一个矩阵表征，称之为标识函数矩阵；天然气管道对天然气流量的影响特性等值为另一个矩阵，称之为管道导纳矩阵；
具体过程如下：
标识函数矩阵：
在天然气系统中，为了表征管道内天然气的流动方向，本发明引入一个标识函数，其主要作用是解决天然气流量方程中含有开方问题，为保证开方数为正，必须采用一个标识函数来予以修正。
标识函数amn，其函数值与天然气站点m和n的站点压力相对大小有关。若两个站点没有管道相连，则amn＝0；若两个站点有管道相连，相对于站点m而言，如果站点m的站点压力pm比与其相连站点n的站点压力pn大，那么amn＝+1，反之amn＝-1，具体定义式如式(1)所示。
每条天然气管道对应的一个标识函数amn，可以生成标识函数矩阵amn，amn是一个方阵，其阶数等于天然气系统的管道数量。
管道导纳矩阵：
天然气输配管道存在管道阻力，可以用rmn来表征连接天然气站点m和站点n的管道对流过天然气的阻力大小，称之为管道阻抗。管道阻抗rmn的大小与天然气管道的长度、内径、倾斜度和内壁摩擦系数有关，同样也受管道内天然气的温度、密度和压缩系数的影响；除了以上管道和天然气物理特性的影响以外，管道阻抗rmn还随天然气管道压力等级的不同，存在数量级上的差异，低压管道(p≤0.01mpa)阻抗表达式如式(2)所示，中高压管道(p≥0.01mpa)阻抗表达式如(3)所示。
在式(2)和式(3)中，d为天然气管道内径，m；l为天然气管道长度，km；ρ0为天然气密度，kg/m3；t为天然气温度，k；t0为标准状态绝对温度，即273.15k，λ为天然气管道内壁摩擦阻力系数。值得强调的是，由于城市天然气管道大多属于地埋式管道，其管道首尾两端的高度差较管道长度很小，因此可以近似认为天然气管道是水平，而管道倾斜度主要是针对给城市高层建筑供气时的管道阻抗的求解，由于本发明是针对电力系统和天然气系统协同运行的问题，所以只考虑城市地埋式天然气管道。
相连的两个天然气站点之间的管道阻抗可以用rmn来表征，不相连的两个天然气站点不存在管道阻抗，但对其进行数学描述时会产生困难，仅能用无穷大来刻画，这不利于天然气流量方程的建模。为了解决这个问题，本发明提出了管道导纳的概念，来简化天然气管网的描述，用cmn表示，定义管道导纳cmn与管道阻抗rmn的大小关系如式(4)所示。
每条天然气管道对应的一个管道导纳cmn，可以生成管道导纳矩阵cmn，cmn是一个方阵，其阶数等于天然气系统的管道数量。
通过标识函数矩阵amn和管道导纳矩阵cmn，可以表征整个天然气系统的管网结构和实际运行时天然气管道对天然气流量的影响特性。
步骤ii：等值电网输电线路参数：
将电力系统的输电线路和变压器等值为节点导纳矩阵，具体过程如下：
电力系统在正常的运行状态下三相对称，所以在等值时认为三相线路参数完全相同，单相等值电路代表三相。短电力线路一般指长度不大于100km、额定电压小于60kv的架空电力线路，由于电压等级低，所以可以只考虑阻抗，忽略输电线路上的电导和电纳(g＝0，b＝0)，即z＝r+jx；中等长度电力线路是指长度在100～300km之间、线路电压为110～220kv的架空电力线路，或者长度不大于100km的电缆电力线路，该线路进行等值时需要考虑电阻r、电抗x、电导g和电纳b四个参数，即z＝r+jx，y＝g+jb，由这四个参数可以构成π形两种等值电路；长电力线路一般指长度大于300km的架空电力线路或者长度大于100km的电缆电力线路，对这种线路需要考虑其分布参数特性。本发明主要针对电力系统和天然气系统协同运行的问题，涉及的电力系统为中等长度线路。电力系统中的变压器通过二端口网络特性可以得到变压器的π型等效电路。
对电力线路和变压器的等效模型都是π型等效电路，使得电力系统中不同元件的等效网络参数表达形式相同，所以可以用节点导纳矩阵对电力系统网络参数进行数学描述，它将电力系统输电线路的电阻r、电抗x、电导g和电纳b和变压器的阻抗zt和变比k等参数归纳为一个矩阵yb。yb是一个n×n阶的节点导纳矩阵，阶数n等于网络中除参考节点之外的节点总数。
通过节点导纳矩阵yb，可以表征电力系统的输电线路参数。
步骤iii：天然气系统中处理内压缩机的实用化方法：
根据压缩机的耗能类型可以将压缩机分为燃气压缩机和电动压缩机，处理压缩机的实用化方法引入压缩机耗能的度量，并将压缩机从天然气管网中去除，改变管网结构，相应的站点等效天然气气负荷和天然气气源；具体过程如下：
在天然气系统中，随天然气输配距离的增大，输气管道内的压力也会下降。因此，在天然气系统中要增设压缩机，用来增加天然气的压力，来补偿燃气输配过程中的压力损耗。压缩机可以类比变压器，前者为提升天然气系统中站点的压力，后者为改变电力系统中母线的电压，但压缩机的处理方式不同于变压器，在电力系统建模时，将变压器等值为π型电路，作为网络参数的一部分，而压缩机的特性复杂，不能作为管道网络参数融入管道导纳矩阵cmn中，因此要对压缩机进行单独建模。
压缩机属于天然气系统的非管道支路，由于升压大小和流过压缩机的流量大小对压缩机耗能都有影响，而燃气轮机的耗能又对天然气流量有反影响，所以将压缩机支路的流量方程求解混在天然气系统中，同管道支路混合计算会造成交叉求解，使得计算复杂，运算量大，不利于处理规模较大的天然气系统。因此，本发明提出的实用化方法是将含压缩机的支路去掉，从而改变天然气管网的拓扑结构，等效得到新的天然气管网拓扑结构。在去掉压缩机后，将压缩机的出口站点处，等效注入压缩机的输出天然气流量，作为该站点的天然气源注入；将压缩机的入口站点处，等效施加压缩机的流入天然气流量，作为该站点的天然气负荷。
现有天然气系统中的压缩机一般分为两种：一种是由燃气轮机供能，消耗天然气流量，提升管道内天然气的压力；另一种是由电动机供能，消耗电功率，提升管道内天然气的压力。本发明提出的实用化方法，既不直接研究燃气压缩机的升压大小和消耗天然气流量关系，也不直接研究电动压缩机的升压大小和耗电量的关系，而是引入压缩机耗能的度量，寻找两种压缩机的升压大小与所需能量e之间的关系，使得两种压缩机具有同等的处理结果，在此基础上，构建为压缩机供能的燃气轮机和电动机输出的能量e与消耗天然气流量和电功率的关系。
燃气压缩机升压消耗的能量egas与流过燃气压缩机的天然气流量和升压大小的关系如式(5)所示。
式(5)中，egas为燃气压缩机在单位时间内消耗能量，mw；ηgas为燃气压缩机的效率；λg为天然气的比热容，λg＝1.309[j/(kg·k)]；t0为标准状态绝对温度，即273.15k；tg为天然气温度，k；p0为天然气的标准压力，即1.01325bar；zg为天然气的压缩系数；fmn为流出压缩机的流量，m3/h；pin和pout分别为压缩机入口压力和出口压力，bar。
所需的燃气压缩机在单位时间内消耗能量egas由拖动压缩机的燃气轮机消耗天然气供给，燃气轮机的耗量特性方程为式(6)：
式(6)中，α、β和γ为燃气轮机的耗量系数；fcom为燃气轮机用于带动压缩机消耗的天然气流量，m3/h。
利用实用化方法将燃气压缩机去掉后，压缩机入口站点处应添加等效天然气负荷fload＝fin(fin是流入压缩机的天然气流量)，压缩机出口站点处应添加等效的天然气气源fsource＝fout(fout是流出压缩机的天然气流量)，值得强调的是，fload＞fsource，差值fcom＝fload-fsource用于供升压耗能。经过该处理使得去掉压缩机等效变换后的系统流量情况与原拓扑结构的系统流量情况相同。
电动压缩机单位时间消耗能量eee与流过压缩机的流量和升压大小的关系，如式(7)所示。
式(7)中，eee为电动压缩机单位时间消耗能量，mw；ηee为电动压缩机的效率。
所需的电动压缩机单位时间消耗能量eee由电动机供给，电动机消耗的电功率pcom和电动压缩机在单位时间内消耗能量eee的关系如式(8)所示。
由于电动压缩机不消耗天然气流量，所以流入电动压缩机的流量和流出电动压缩机的流量大小应该是相等的，即fin＝fout。进而，在电动压缩机入口站点等效的气负荷fload和出口站点等效的气源fsource也是相等的，即fload＝fsource。
步骤iv：建立天然气管网和电网运行的数学模型：
建立天然气管道中流过天然气流量和电网线路中流过电功率的数学模型表达；
天然气系统：
天然气管网中流过的天然气流量需要通过水力分析，由水力方程可以推导得出天然气管网的运行方程，流体动力学方程、连续性方程和流体状态方程分别如式(9)～(11)所示。
流体动力学方程：
连续性方程：
流体状态方程：
p＝zρrt(11)
其中：ω为天然气流速，m/s；ρ为天然气密度，kg/m3；p为天然气压力，pa；t为天然气温度，k；t为时间，s；x为管道轴向坐标，m；g为重力加速度，m/s2；α为管道与水平面夹角，°；λ为管道内壁摩擦系数；d为管道内径，m；z为压缩系数；r为气体常数，j/(kg·k)。
稳定流量基本方程如式(12)所示。
式(12)中，f表示天然气流量，m3/h，下标“0”表示基准状态：基准压力p0＝101325pa，基准温度t0＝293.15k，基准压缩系数z0是p0和t0的函数。由于本发明提出的方法不考虑管道敷设高程变化，所以认为管道与水平面夹角α＝0°。
通过两边取积分，并用管道导纳表示实际管道情况对天然气的影响，可以得到天然气系统中的管道支路的流量方程，如式(13)所示。
其中，fmn为天然气系统中的管道支路的流量，hmn为管道压力降，它表征两个站点之间绝对压力的差值，管道压力降的表达式如式(14)所示。
其中，各个站点的压力值和各个管道中的流量值均有上下限的要求，如式(15)和式(16)所示。
pimin≤pi≤pimax(15)
fmnmin≤fmn≤fmnmax(16)
式中，pi为任意一个站点的压力值，i取1,2,…ngas，ngas为天然气站点的总数，pimin为第i个站点压力值的下限值，pimax为第i个站点压力值的上限值；
fmn为站点m和站点n之间连接管道中流过的天然气流量，fmnmin为站点m和站点n之间连接管道中流过的天然气流量的下限值，fmnmax为站点m和站点n之间连接管道中流过的天然气流量上限值；
电力系统：
建立电网潮流方程的理论基础是电路理论中的节点分析法，节点电压用极坐标形式表示，如式(17)和式(18)所示。
pij＝uiuj(gijcosδij+bijsinδij)(17)
qij＝uiuj(gijsinδij-bijcosδij)(18)
式中，pij为电力系统中母线i和母线j之间连接上线路流过的有功功率，qij为电力系统中母线i和母线j之间连接上线路流过的无功功率ui为电力系统母线i的电压幅值，uj为电力系统母线j的电压幅值，gij为电力系统中母线i和母线j之间连接上线路的电导，bij为电力系统中母线i和母线j之间连接上线路的电纳，δij为电力系统母线i和母线j之间的电压相角的差值，gij和bij都是节点导纳矩阵yb中的元素；
其中，各个母线的电压值和各个线路上的功率均有上下限的要求，如式(19)～(21)所示。
uimin≤ui≤uimax(19)
pijmin≤pij≤pijmax(20)
qijmin≤qij≤qijmax(21)
其中，ui为电力系统中任意母线的电压幅值，i为1,2,…nee，nee为电力系统母线的总数，uimin为电力系统中母线i的电压幅值的下限值，uimax为电力系统中母线i的电压幅值的上限值，pijmin为电力系统中母线i和母线j之间连接上线路流过的有功功率的下限值，pijmax为电力系统中母线i和母线j之间连接上线路流过的有功功率的上限值，qijmin为电力系统中母线i和母线j之间连接上线路流过的无功功率的下限值，qijmax为电力系统中母线i和母线j之间连接上线路流过的无功功率的上限值。
步骤v：天然气管网和电网协同运行计算：
对天然气管网和电网的实时运行状态进行统一计算，分析每个时刻天然气管网和电网协同运行的状态，验证该方法的有效性和实用性。具体过程如下：
天然气管网和电网协同运行的计算来验证本发明提出方法的有效性和实用性。通过牛顿-拉夫逊法，对天然气管网和电网协同运行进行求解。天然气管网和电网的不平衡方程如式(22)～(24)所示。
其中，fi为天然气站点注入的天然气流量，fee是燃气轮机组发电消耗的天然气流量，只有接燃气轮机组的天然气站点才存在该项；pgas和pcom分别是燃气轮机组发电向电力系统注入的有功功率和电动压缩机给天然气流量升压从电力系统中消耗的电功率，只有连接相应耦合单元的电力系统节点的不平衡方程才存在相应项。
下面以一个实施例说明本发明的具体实施步骤。
采用12节点天然气系统和ieee14节点电力系统的协同运行计算进行算例验证。
表112节点天然气系统节点原始数据
表212节点天然气系统管道原始数据
表3压缩机原始数据
通过协同运行分析，从图1中可以分析得出，在高用电量时，天然气系统中流量节点的压力值比低用电量时小很多，这是由于燃气轮机组给电网的发电量增大导致其耗气量的增大，需要天然气系统供应更多的流量，导致输配流量增大，在管道首端压力恒定时，管道末端的压力值越小，流过的流量越大，所以天然气系统整体的节点压力值都有下降。
管道中的天然气流量计算结果如表4所示。
表4天然气流量计算结果
管道流量单位：105m3/h
最后，着重强调压缩机处理实用化方法的有效性。管道2和管道4的流量存在一定的差值，这是因为管道2的末端接燃气压缩机的入口，流量大小即为压缩机入口等效的气负荷流量，管道4的首端接燃气压缩机的出口，流量大小即为压缩机出口等效的气源流量，该差值就是燃气压缩机用于升压而消耗掉的流量fcom。同样连接电动压缩机的管道8和管道9流过的流量则相等，这是因为电动压缩机升压消耗的能量来自于电网的电功率，不消耗天然气流量。通过比较分析，本发明提出的压缩机处理实用化方法具有适用性和有效性。
本发明有两大改进：第一、对天然气管道特性的建模更加贴近实际，从实际管道的物理特性抽象得到；第二、对天然气系统中的压缩机建模提出了实用性处理方法，使得天然气系统和电力系统的建模更加统一，求解过程中计算量更小，计算速度更快。