一种电?气耦合系统中天然气网的区间潮流计算方法
【专利摘要】本发明涉及一种电?气耦合系统中天然气网的区间潮流计算方法，属于多能流耦合系统的运行安全分析技术领域。本方法采取区间潮流的分析方法，在多种能源系统耦合运行的发展趋势下，考虑了风电不确定性对其他与电网相耦合的能源系统(天然气管网系统)稳态运行的影响，为耦合能源系统的安全稳定控制提供了依据；另一方面避免了大量的数据统计和复杂的数学模型及计算，同时该方法只需要区间信息，不需要估计隶属度信息等，避免了因人为假设的主观性带来的误差。该方法可以应用于电?气耦合多能流能量管理中，当系统存在安全问题或潜在风险时给出安全警告，有助于提高电?气耦合多能流系统运行的安全性。
【专利说明】
一种电-气耦合系统中天然气网的区间潮流计算方法
技术领域
[0001] 本发明涉及一种电-气耦合系统中天然气网的区间潮流计算方法，尤其涉及一种 基于风电不确定的电-气耦合系统中天然气网的区间潮流计算方法，属于多能流耦合系统 的运行安全分析技术领域。
【背景技术】
[0002] 由于燃气发电机的成本较低，对环境的影响较小，燃气站的建造周期较短等优势， 在世界范围内，天然气应用中用于发电的比例逐渐上升，在电力系统能源供给中也占据了 越来越重要的地位，因此，电力网络与天然气网络之间的耦合程度逐渐加深，相互之间的影 响与依赖性也愈来愈强。鉴于能源资源条件的限制，我国的能源消费结构在未来相当长的 一段时间内仍然以煤炭为主，但是环境保护及可持续发展战略都要求我国必须加快改善能 源结构的步伐。在一次能源消费结构中，提高天然气等可再生的清洁能源的比重，而降低煤 炭的比重。
[0003] 另一方面，由于能源紧缺和环境污染的问题日益突出，风能因其环境友好、技术成 熟、零燃料成本且可持续等优势，成为了最具竞争力的能源形式之一。近些年来风电在我国 得到了快速的发展，各种能源在发电中的融合成为电力系统未来发展的必然趋势。然而，风 力发电机的出力会受到实时的风力资源随机变化的限制，无法提供连续稳定的功率，因此 风电是一种具有波动性、随机性、间歇性和难以调度性的不可靠电源。在如今多能源耦合系 统统一运行模式的趋势下，各种能源之间的相互联系和交互作用更加紧密，迫切需要研究 风力发电机组并网后其不确定性对其他能源网络安全性的影响。
【发明内容】
[0004] 本发明的目的是提出一种考虑风电不确定性的电-气耦合系统天然气网区间潮流 分析方法，避免了大量的数据统计、复杂的计算以及过多人为假设造成的偏差，分析风电接 入后的电-气耦合系统中天然气网的稳态运行区间。
[0005] 本发明提出的电-气耦合系统中天然气网的区间潮流计算方法，包括以下步骤： [0006] (1)建立电-气耦合系统中耦合运行的稳态数学模型，包括：
[0007] (1-1)-个电-气親合系统中的电力系统潮流方程如下：
[0009]其中，P1为电力系统中第i个节点的注入有功功率，Q1为电力系统中第i个节点的注 入无功功率，为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，为与电力 系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统 调度中心获取；
[0010] (1-2)-个电-气耦合系统的天然气网中管道的水力方程如下：
[0012] 其中，fkm为天然气网中第k个节点和第m个节点之间的管道中天然气体积流量，pk， 别为第k个节点和第m个节点的压强，D km、Lkl别为第k个节点和第m个节点之间管道km 的管道直径和管道长度，F为管道内壁的摩擦系数，
计算得 到，Ef为管道的效率系数，Ef取值0.92,Re为雷诺数，由公式Re = ρν(1/μ计算得到，P为天然气 密度，v为天然气流速，μ为天然气黏性系数，d为天然气管道直径，γ g为天然气比重，0< γ c <l，Ta为天然气平均温度，Τη、ρη分别为标准状态下天然气的温度和压力，Τ η、ρη的取值分别 为288Κ、0 · IMpa，Zg为天然气平均可压缩系数，0 · 9<Zg< 1 · 5，天然气网水力方程中，当 (A2 - k)之0时，上式中的 sgnp (Pk，Pm) = 1，当(A - k ) < 〇时，Sgnp (Pk，Pm) = -1;
[0013] (1-3)-个电-气耦合系统中的天然气网中压缩机的能量消耗方程如下：
[0015] 其中，pk，pm分别为天然气网中第k个节点和第m个节点的压强，BHPkm为第k个节点 和第m个节点之间的压缩机的能量消耗，为该压缩机的入口体积流量，η。为该压缩机的总 效率， Ck为该压缩机的多变系数，n。和(^从压缩机的出厂说明书获取；
[0016] (1-4) 一个通过燃气轮机耦合的电力系统与天然气网之间的耦合方程如下： Ρη"·=^/1，+α/Τω.+€\,
[0017] 其中，fTur为燃气轮机的燃气的体积流量，PTur为燃气轮机的有功功率输出，&、&和 C3分别为燃气轮机的燃气系数，由燃气轮机的出厂说明书获取；
[0018] (1-5)-个电-气耦合系统中的天然气网的节点气流平衡方程如下：
[0019] Acf = L,
[0020] 其中：AC为天然气网的节点-支路矩阵，f为天然气网的支路体积流量，L为天然气 网的节点气体负荷，L根据天然气网历史运行数据得到；
[0021] (2)通过电力系统中风电场的历史纪录数据，获取风电场的历史最小风速X和最大 风速f的风速变化区间[_\:.?? ;
[0022] (3)比对电力系统中风电场的风力涡轮机的功率曲线和上述风速的变化区间 [1>]，.得到风电功率的变化区间
[0023] /?
[0024] 其中，为风电场第1个风电机组的有功功率，风力涡轮机的功率曲线由风力涡轮 机的出厂说明书获取；
[0025] (4)设定电力系统和天然气网稳态安全运行的约束条件，包括：
[0026] (4-1)电力系统中非风力发电机组的输出功率PFM、于或等于该非风力发电机组 出厂铭牌上给出的最大功率/is：:
[0027] ΡΓ < ΡΣ 5
[0028] (4-2)电力系统第i个节点的电压幅值Ui在设定的电力系统安全运行电压的上、下 限值少、巧之间运行，少为第i个节点额定电压的0.95倍，g；为第i个节点额定电压的1.05倍：
[0029] u^u^Ur，
[0030] (4-3)电力系统中第1条线路的传输容量sd、于或等于设定的电力系统安全运行 传输容量的最大值属：
[0031]
[0032] (4-4)天然气网中第k个节点的压强pk在设定的管道安全运行气压的上、下限值 ￡k、 Λ内：
[0033] Ih ^ Ρ? ^ Ρι ·
[0034] (4-5)天然气网中第b个管道的流量fb在设定的管道安全运行流量的上、下限值f b、 Λ内：
[0035] 〇 < ./；, ^ Ζ ；
[0036] (4-6)天然气网中气源供气量匕小于或等于该气源能提供天然气流的最大值 f s, max ：
[0037] fs<fs,max；
[0038] (4-7)天然气网中的压缩机安全运行时需要满足以下约束：
[0040] 其中：S是该压缩机的升压比，Smax是该压缩机的最大升压比，Smax由压缩机的出厂 铭牌获取，为该压缩机的入口的体积流量，为该压缩机的入口最大允许体积流量， 由压缩机的出厂铭牌获取，P°ut为压缩机的出口压力，P。，为压缩机的出口最大允许压 力，Pc,max由压缩机的出厂铭牌获取；
[0041] (5)利用内点法，在上述步骤(3)给定的风电功率的变化区间，根据上 述步骤（1)的电-气耦合系统耦合运行的稳态方程和上述步骤(4)的电力系统与天然气网稳 态约束方程，求解得到天然气网区间潮流，即得到电-气耦合系统中天然气网各节点的压强 P、各支路的体积流量f以及气源供气量f s的运行范围。
[0042] 本发明提出的电-气耦合系统中天然气网的区间潮流计算方法，其特点和效果是： 本发明采取区间潮流的分析方法，一方面在多种能源系统耦合运行的发展趋势下，考虑了 风电不确定性对其他能源网(天然气管网系统)稳态运行的影响，为耦合能源系统的安全稳 定控制提供了依据；另一方面避免了大量的数据统计和复杂的数学模型及计算，同时该方 法只需要区间信息，不需要估计隶属度信息等，避免了因人为假设的主观性带来的误差。该 方法可以应用于电-气耦合多能流能量管理中，当系统存在安全问题或潜在风险时给出安 全警告，有助于提高电-气耦合多能流系统运行的安全性。
【具体实施方式】
[0043] 本发明提出的电-气耦合系统中天然气网的区间潮流计算方法，包括以下步骤：
[0044] (1)建立电-气耦合系统中耦合运行的稳态数学模型，包括：
[0045] (1-1) 一个电-气耦合系统中的电力系统潮流方程如下：
[0047]其中，P1为电力系统中第i个节点的注入有功功率，Q1为电力系统中第i个节点的注 入无功功率，Gij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，Bij为与电力 系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统 调度中心获取；
[0048] (1-2)-个电-气耦合系统的天然气网中管道的水力方程如下：
[0050] 其中，fkm为天然气网中第k个节点和第m个节点之间的管道中天然气体积流量，pk， 别为第k个节点和第m个节点的压强，D km、Lkl别为第k个节点和第m个节点之间管道km 的管道直径和管道长度，F为管道内壁的摩擦系数，
计算得 到，Ef为管道的效率系数，Ef取值0.92,Re为雷诺数，由公式Re = ρν(1/μ计算得到，P为天然气 密度，v为天然气流速，μ为天然气黏性系数，d为天然气管道直径，γ g为天然气比重，0< γ c <l，Ta为天然气平均温度，Τη、ρη分别为标准状态下天然气的温度和压力，Τ η、ρη的取值分别 为288Κ、0 · IMpa，Zg为天然气平均可压缩系数，0 · 9<Zg< 1 · 5，天然气网水力方程中，当 (A2 - k)》0时，上式中的sgnP(pk，pm) = 1，当(W - ρ!) < 〇时，sgnP(pk，pm) = -1;
[0051] (1-3)-个电-气耦合系统中的天然气网中压缩机的能量消耗方程如下：
[OO53]其中，pk，pm分别为天然气网中第k个节点和第m个节点的压强，BHPkm为第k个节点 和第m个节点之间的压缩机的能量消耗，/T为该压缩机的入口体积流量，η。为该压缩机的总 效率，Ck为该压缩机的多变系数， n。和(^从压缩机的出厂说明书获取；
[0054] (1-4) 一个通过燃气轮机耦合的电力系统与天然气网之间的耦合方程如下： ^=CJl+Czfni,+C3,
[0055] 其中，fTur为燃气轮机的燃气的体积流量，PTur为燃气轮机的有功功率输出，&、&和 C3分别为燃气轮机的燃气系数，由燃气轮机的出厂说明书获取；
[0056] (1-5)-个电-气耦合系统中的天然气网的节点气流平衡方程如下：
[0057] Acf = L,
[0058] 其中：Ag为天然气网的节点-支路矩阵，f为天然气网的支路体积流量，L为天然气 网的节点气体负荷，L根据天然气网历史运行数据得到；
[0059] (2)通过电力系统中风电场的历史纪录数据，获取风电场的历史最小风速X和最大 风速f的风速变化区间[工，巧：
[0060] (3)比对电力系统中风电场的风力涡轮机的功率曲线和上述风速的变化区间 [z，f|，得到风电功率的变化区间[f:
[0061] 墨響％'，
[0062] 其中，为风电场第1个风电机组的有功功率，风力涡轮机的功率曲线由风力涡轮 机的出厂说明书获取；
[0063] (4)设定电力系统和天然气网稳态安全运行的约束条件，包括：
[0064] (4-1)电力系统中非风力发电机组的输出功率PFM、于或等于该非风力发电机组 出厂铭牌上给出的最大功率/LT:
[0065] ΙΤ?
[0066] (4-2)电力系统第i个节点的电压幅值Ui在设定的电力系统安全运行电压的上、下 限值虫、钇之间运行，虫为第i个节点额定电压的0.95倍，$为第i个节点额定电压的1.05 倍：
[0067]
[0068] (4-3)电力系统中第1条线路的传输容量Sd、于或等于设定的电力系统安全运行 传输容量的最大值属：
[0069] S<S,x
[0070] (4-4)天然气网中第k个节点的压强pk在设定的管道安全运行气压的上、下限值 ￡k、 艮内：
[0071 ] Pi ^ Pi
[0072] (4-5)天然气网中第b个管道的流量fb在设定的管道安全运行流量的上、下限值f b、 I内：
[0073] 0 < ./；, < fh ;
[0074] (4-6)天然气网中气源供气量匕小于或等于该气源能提供天然气流的最大值 f s, max ：
[0075] fs<fs,max；
[0076] (4-7)天然气网中的压缩机安全运行时需要满足以下约束：
[0078] 其中：S是该压缩机的升压比，Smax是该压缩机的最大升压比，Smax由压缩机的出厂 铭牌获取，为该压缩机的入口的体积流量，·为该压缩机的入口最大允许体积流量， /J:由压缩机的出厂铭牌获取，P° utS压缩机的出口压力，p。,^为压缩机的出口最大允许压 力，Pc,max由压缩机的出厂铭牌获取；
[0079] (5)利用内点法，在上述步骤(3)给定的风电功率的变化区间内，根据上 述步骤（1)的电-气耦合系统耦合运行的稳态方程和上述步骤(4)的电力系统与天然气网稳 态约束方程，求解得到天然气网区间潮流，即得到电-气耦合系统中天然气网各节点的压强 P、各支路的体积流量f以及气源供气量f s的运行范围。
【主权项】
1. 一种电-气禪合系统中天然气网的区间潮流计算方法，其特征在于该方法包括W下 步骤： (1)建立电-气禪合系统中禪合运行的稳态数学模型，包括： (1-1) 一个电-气禪合系统中的电力系统潮流方程如下：其中，pi为电力系统中第i个节点的注入有功功率，Q1为电力系统中第i个节点的注入无 功功率，Gu为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，Βυ为与电力系统 的节点导纳矩阵Υ中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Υ从电力系统调度 中屯、获取； (1-2)-个电-气禪合系统的天然气网中管道的水力方程如下：其中，fkm为天然气网中第k个节点和第m个节点之间的管道中天然气体积流量，Pk，Pm分 别为第k个节点和第m个节点的压强，〇km、Lkm分别为第k个节点和第m个节点之间管道km的管 道直径和管道长度，F为管道内壁的摩擦系数，F由公式计算得到，Ef 为管道的效率系数，Ef取值0.92，Re为雷诺数，由公式Re = Pvd/μ计算得到，P为天然气密度， V为天然气流速，μ为天然气黏性系数，d为天然气管道直径，丫 G为天然气比重，〇< YG<l，Ta 为天然气平均溫度，Τη、Ρη分别为标准状态下天然气的溫度和压力，Τη、Ρη的取值分别为 288Κ、0 . IMpa，Zg为天然气平均可压缩系数，0.9<Zg< 1 . 5，天然气网水力方程中，当 (拓-Λ )含 0 时，上式中的 Sgnp ( Pk，Pm) = 1，当(记-姑）< 0时，S即P ( Pk，Pm) = -1 ; (1-3)-个电-气禪合系统中的天然气网中压缩机的能量消耗方程如下：其中，Pk，Pm分别为天然气网中第k个节点和第m个节点的压强，B册km为第k个节点和第m 个节点之间的压缩机的能量消耗，记"为该压缩机的入口体积流量，η。为该压缩机的总效率， Ck为该压缩机的多变系数，ric和Ck从压缩机的出厂说明书获取； (1-4) 一个通过燃气轮机禪合的电力系统与天然气网之间的禪合方程如下：其中，fTur为燃气轮机的燃气的体积流量，PTur为燃气轮机的有功功率输出，Cl、C2和C3分 别为燃气轮机的燃气系数，由燃气轮机的出厂说明书获取； (1-5)-个电-气禪合系统中的天然气网的节点气流平衡方程如下： A 江=L， 其中:Ag为天然气网的节点-支路矩阵，f为天然气网的支路体积流量，L为天然气网的节 点气体负荷，L根据天然气网历史运行数据得到； (2) 通过电力系统中风电场的历史纪录数据，获取风电场的历史最小风速X和最大风速 V的风速变化区间[ii，巧； (3) 比对电力系统中风电场的风力满轮机的功率曲线和上述风速的变化区间[Z，巧，得 到风电功率的变化区间[C，，巧。：其中，Pi"为风电场第1个风电机组的有功功率，风力满轮机的功率曲线由风力满轮机的 出厂说明书获取； (4) 设定电力系统和天然气网稳态安全运行的约束条件，包括： (4-1)电力系统中非风力发电机组的输出功率Pigen小于或等于该非风力发电机组出厂 铭牌上给出的最大功率巧置：(4-2)电力系统第i个节点的电压幅值化在设定的电力系统安全运行电压的上、下限值 些、巧之间运行，些为第i个节点额定电压的0.95倍巧为第i个节点额定电压的1.05倍： U. < U < 0,： (4-3)电力系统中第1条线路的传输容量Si小于或等于设定的电力系统安全运行传输容 量的最大值 (4-4)天然气网中第k个节点的压强化在设定的管道安全运行气压的上、下限值￡1^、及 内： 公^公^瓦； (4-5)天然气网中第b个管道的流量fb在设定的管道安全运行流量的上、下限值旬、充 内： 当充； (4-6)天然气网中气源供气量fs小于或等于该气源能提供天然气流的最大值fs,max: f s《f s, max ; (4-7)天然气网中的压缩机安全运行时需要满足W下约束：其中：S是该压缩机的升压比，Smax是该压缩机的最大升压比，Smax由压缩机的出厂铭牌 获取，为该压缩机的入口的体积流量，怎；《X为该压缩机的入口最大允许体积流量，尤;。。1 由压缩机的出厂铭牌获取，P°ut为压缩机的出口压力，Pc,max为压缩机的出口最大允许压力， Pc,max由压缩机的出厂铭牌获取； (5) 利用内点法，在上述步骤(3)给定的风电功率的变化区间[G:,,，G：K]内，根据上述步 骤(1)的电-气禪合系统禪合运行的稳态方程和上述步骤(4)的电力系统与天然气网稳态约 束方程，求解得到天然气网区间潮流，即得到电-气禪合系统中天然气网各节点的压强P、各 支路的体积流量fW及气源供气量fs的运行范围。
【文档编号】G06F19/00GK106096269SQ201610407465
【公开日】2016年11月9日
【申请日】2016年6月12日
【发明人】孙宏斌, 郭庆来, 王彬, 乔铮, 潘昭光
【申请人】清华大学