CN108667007B - 计及电-气耦合系统约束的电压稳定裕度计算方法 - Google Patents

计及电-气耦合系统约束的电压稳定裕度计算方法 Download PDF

Info

Publication number
CN108667007B
CN108667007B CN201810335838.4A CN201810335838A CN108667007B CN 108667007 B CN108667007 B CN 108667007B CN 201810335838 A CN201810335838 A CN 201810335838A CN 108667007 B CN108667007 B CN 108667007B
Authority
CN
China
Prior art keywords
power
gas
node
electric
power system
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.)
Active
Application number
CN201810335838.4A
Other languages
English (en)
Other versions
CN108667007A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201810335838.4A priority Critical patent/CN108667007B/zh
Publication of CN108667007A publication Critical patent/CN108667007A/zh
Priority to PCT/CN2018/113635 priority patent/WO2019200891A1/zh
Application granted granted Critical
Publication of CN108667007B publication Critical patent/CN108667007B/zh
Priority to US17/021,445 priority patent/US20200410145A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
    • H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • 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/381—Dispersed generators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • 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]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30—Reactive power compensation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Primary Health Care (AREA)
  • Marketing (AREA)
  • Human Resources & Organizations (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种计及电‑气耦合系统约束的电压稳定裕度计算方法，属于电力系统中考虑多能流耦合特性的安全分析与评估技术领域。本发明充分考虑电力系统与天然气系统的紧密耦合，求出了耦合系统中电力系统的电压稳定裕度。一方面充分考虑了天然气系统的安全与容量约束对电力系统的影响，另一方面也根据应用地区的实际情况考虑了电力负荷与天然气负荷的相关性对电压稳定裕度的影响，避免了传统的计算方法中的单纯考虑电力系统的约束而导致的计算结果过于乐观。该方法可用于电力系统的运行风险分析中，为电力系统的运行管理人员提供风险评估指标，有利于降低潜在风险，提高系统运行的安全性。

Description

计及电-气耦合系统约束的电压稳定裕度计算方法
技术领域
本发明涉及一种计及电-气耦合系统约束的电压稳定裕度计算方法，属于电力系统中考虑多能流耦合特性的安全分析与评估技术领域。
背景技术
由于燃气发电机的成本较低、对环境的破坏力较小、响应速度快，燃气站的建造周期较短等巨大优势，天然气已在世界范围内成为供电能源的重要组成部分。因此随着天然气在电力系统一次能源供应中的比例日益增加，天然气的可靠供应对于电力系统的安全具有至关重要的要作用。
然而，天然气与煤炭等可以较大规模储存的能源不同，其供应的方式大多通过管道进行远距离传输。一方面，由于压力安全约束，管道传输的天然气流量具有一定的限制。另一方面，天然气负荷在年、月、日间的波动性非常剧烈，而在很多国家规定中，其他商业、民用天然气负荷的优先级高于燃气电厂的燃气负荷。因此，电力系统的燃气供应受限于天然气系统的管道传输容量和其他天然气负荷，仅仅考虑电力系统约束的电压稳定裕度的计算方法不再适用，亟须提出一种新的计及电-气耦合系统约束的电力系统电压稳定裕度计算方法。
发明内容
本发明的目的是提出一种计及电-气耦合系统约束的电力系统电压稳定裕度计算方法，以避免因未考虑天然气系统的安全约束以及天然气负荷的影响计算出的电压稳定裕度过于乐观，存在潜在风险。
本发明提出的计及电-气耦合系统约束的电压稳定裕度计算方法，包括以下步骤：
(1)建立电-气耦合系统稳态安全运行的等式约束方程，包括：
(1-1)一个电-气耦合系统中电力系统的潮流方程如下：
其中，PGi为电力系统中第i个节点的注入有功功率，PLi为电力系统中第i个节点的注出有功功率，QGi为电力系统中第i个节点的注入无功功率，QLi为电力系统中第i个节点的注出无功功率，Vi、Vj分别为电力系统中第i个节点、第j个节点的电压幅值，θi、θj分别为电力系统中第i个节点、第j个节点的电压相角，Gij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，Bij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统调度中心获取，Ne为电力系统中所有节点的个数，NPQ为电力系统中给定了有功功率P和无功功率Q的PQ节点的个数；
(1-2)一个电-气耦合系统中天然气系统中管道的水力方程如下：
其中，fkm为天然气系统中第k个节点和第m个节点之间的管道中天然气体积流量，pk，pm分别为第k个节点和第m个节点的压强，Ckm为第k个节点和第m个节点之间管道km的阻力系数，从管道的设计报告中获取，上述天然气系统中管道的水力方程中，当时，上式中的sgnp(pk,pm)＝1，当时，sgnp(pk,pm)＝-1；
(1-3)一个电-气耦合系统中通过燃气轮机耦合的电力系统与天然气系统之间的耦合方程如下：
μG×LG×Hgas＝PG，
其中，LG为燃气轮机的燃气负荷，PG为燃气轮机的有功功率输出，Hgas为天然气的燃烧热值，取值为37.59MJ/m3，μG为燃气轮机的效率系数，由燃气轮机的出厂说明书获取；
(1-4)一个电-气耦合系统中天然气系统的节点气流平衡方程如下：
其中，Lsm为天然气系统中第m个节点的注入体积流量，LLm为天然气系统中第m个节点的注出体积流量；
(2)选择负荷裕度指标λ作为电压稳定裕度的指标，并在下述三种负荷增长方式中选择一种：①其他负荷保持不变，单一负荷的有功功率和无功功率保持原功率因数同时增加；②其他负荷保持不变，选定区域内负荷的有功功率和无功功率保持原功率因数同时增加；③所有负荷的有功功率和无功功率保持原功率因数同时增加；
(3)利用负荷裕度指标λ建立电-气耦合系统的连续潮流模型：
(3-1)建立一个电-气耦合系统中电力系统注入与注出功率变化方程如下：
其中，PLi0为初始时刻节点i的注出有功功率，PGi0为初始时刻节点i的注入有功功率，QLi0为初始时刻节点i的注入无功功率，Ne为电力系统中的节点个数，NPQ为电力系统中PQ节点的个数；
(3-2)建立一个电-气耦合系统中天然气系统中天然气负荷的变化方程如下：
LLm(λ)＝(1+rλ)LLm0，
其中，LLm0为初始时刻第m个节点的注出体积流量，通过天然气系统的运行数据获得；r为电力系统燃气负荷与天然气系统负荷的相关系数，与地域、气候、季节等相关，由当地能源统计部门的数据获得；
(3-3)将步骤(3-1)和(3-2)中的连续变化方程代入步骤(1-1)和(1-4)中的方程后得到如下方程：
(4)设定电-气耦合系统稳态安全运行的不等式约束条件，包括：
(4-1)电力系统中发电机组的输出有功功率Pgen大于或等于0，小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-2)电力系统中发电机组的输出无功功率大于或等于该发电机组出厂铭牌上给出的最小功率小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-3)电力系统第i个节点的电压幅值Ui在设定的电力系统安全运行电压的上限值和下限值U i之间运行，U i取第i个节点额定电压的0.9倍或0.95倍，为第i个节点额定电压的1.1倍或1.05倍，即：
(4-4)天然气系统中第k个节点的压强pk在设定的管道安全运行气压的上限值、下限值p k、内，即：
(4-5)天然气系统中气源供气量Ls大于或等于0，小于或等于该气源能提供天然气流的最大值Ls,max，即：
0≤Ls≤Ls,max；
(5)利用优化法(如内点法)或迭代法(如牛顿法)求解λ为0时由步骤(1)和步骤(3-3)构建的潮流方程F(X)，得到初始潮流解Xt(Vt,θt,λt)，其中下标t表示当前计算点；
(6)在上述初始解Xt处求得切向量dXt(dVt,dθt,dλt)，并设置潮流解变化的步长h，得到预测值X't+1(V't+1,θ't+1,λ't+1)，其中下标t+1表示下一个计算点：
X't+1＝Xt+h·dXt；
(7)以X't+1为初始点，重新计算由步骤(1)和步骤(3-3)构建的潮流方程，得到修正值Xt+1，判断Xt+1是否满足步骤(4)中的约束且dλt＞0，若同时满足步骤(4)的约束和dλt＞0的条件，则以Xt+1作为初始解Xt，返回执行步骤(6)，若不满足步骤(4)的约束或不满足dλt＞0的条件，则进一步判断Xt+1是否满足dλt/λt＜ε且dλt＞0，若不满足该条件，则重新调整步长h，返回执行步骤(6)，若满足该条件，输出此时的λ，作为计及电-气耦合系统约束的电压稳定裕度。
本发明涉及一种电-气耦合系统约束的电力系统电压稳定裕度计算方法，其特点和效果是：
本发明方法充分考虑电力系统与天然气系统的紧密耦合，求出了耦合系统中电力系统的电压稳定裕度。一方面充分考虑了天然气系统的安全与容量约束对电力系统的影响，另一方面也根据应用地区的实际情况考虑了电力负荷与天然气负荷的相关性对电压稳定裕度的影响，避免了传统的计算方法中的单纯考虑电力系统的约束而导致的计算结果过于乐观。该方法可用于电力系统的运行风险分析中，为电力系统的运行管理人员提供风险评估指标，有利于降低潜在风险，提高系统运行的安全性
具体实施方式
本发明提出的计及电-气耦合系统约束的电压稳定裕度计算方法，包括以下步骤：
(1)建立电-气耦合系统稳态安全运行的等式约束方程，包括：
(1-1)一个电-气耦合系统中电力系统的潮流方程如下：
其中，PGi为电力系统中第i个节点的注入有功功率，PLi为电力系统中第i个节点的注出有功功率，QGi为电力系统中第i个节点的注入无功功率，QLi为电力系统中第i个节点的注出无功功率，Vi、Vj分别为电力系统中第i个节点、第j个节点的电压幅值，θi、θj分别为电力系统中第i个节点、第j个节点的电压相角，Gij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，Bij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统调度中心获取，Ne为电力系统中所有节点的个数，NPQ为电力系统中给定了有功功率P和无功功率Q的PQ节点的个数；
(1-2)一个电-气耦合系统中天然气系统中管道的水力方程如下：
其中，fkm为天然气系统中第k个节点和第m个节点之间的管道中天然气体积流量，pk，pm分别为第k个节点和第m个节点的压强，Ckm为第k个节点和第m个节点之间管道km的阻力系数，从管道的设计报告中获取，上述天然气系统中管道的水力方程中，当时，上式中的sgnp(pk,pm)＝1，当时，sgnp(pk,pm)＝-1；
(1-3)一个电-气耦合系统中通过燃气轮机耦合的电力系统与天然气系统之间的耦合方程如下：
μG×LG×Hgas＝PG，
其中，LG为燃气轮机的燃气负荷(体积流量)，PG为燃气轮机的有功功率输出，Hgas为天然气的燃烧热值，取值为37.59MJ/m3，μG为燃气轮机的效率系数，由燃气轮机的出厂说明书获取；
(1-4)一个电-气耦合系统中天然气系统的节点气流平衡方程如下：
其中，Lsm为天然气系统中第m个节点的注入体积流量，LLm为天然气系统中第m个节点的注出体积流量；
(2)选择负荷裕度指标λ作为电压稳定裕度的指标，并在下述三种负荷增长方式中选择一种：①其他负荷保持不变，单一负荷的有功功率和无功功率保持原功率因数同时增加；②其他负荷保持不变，选定区域内负荷的有功功率和无功功率保持原功率因数同时增加；③所有负荷的有功功率和无功功率保持原功率因数同时增加；
(3)利用负荷裕度指标λ建立电-气耦合系统的连续潮流模型：
(3-1)建立一个电-气耦合系统中电力系统注入与注出功率变化方程如下：
其中，PLi0为初始时刻节点i的注出有功功率，PGi0为初始时刻节点i的注入有功功率，QLi0为初始时刻节点i的注入无功功率，Ne为电力系统中的节点个数，NPQ为电力系统中PQ节点的个数；
(3-2)建立一个电-气耦合系统中天然气系统中天然气负荷的变化方程如下：
LLm(λ)＝(1+rλ)LLm0，
其中，LLm0为初始时刻第m个节点的注出体积流量，通过天然气系统的运行数据获得；r为电力系统燃气负荷与天然气系统负荷的相关系数，与地域、气候、季节等相关，由当地能源统计部门的数据获得；
(3-3)将步骤(3-1)和(3-2)中的连续变化方程代入步骤(1-1)和(1-4)中的方程后得到如下方程：
(4)设定电-气耦合系统稳态安全运行的不等式约束条件，包括：
(4-1)电力系统中发电机组的输出有功功率Pgen大于或等于0，小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-2)电力系统中发电机组的输出无功功率大于或等于该发电机组出厂铭牌上给出的最小功率小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-3)电力系统第i个节点的电压幅值Ui在设定的电力系统安全运行电压的上限值和下限值U i之间运行，U i取第i个节点额定电压的0.9倍或0.95倍，为第i个节点额定电压的1.1倍或1.05倍，即：
(4-4)天然气系统中第k个节点的压强pk在设定的管道安全运行气压的上限值、下限值p k、内，即：
(4-5)天然气系统中气源供气量Ls大于或等于0，小于或等于该气源能提供天然气流的最大值Ls,max，即：
0≤Ls≤Ls,max；
(5)利用优化法(如内点法)或迭代法(如牛顿法)求解λ为0时由步骤(1)和步骤(3-3)构建的潮流方程F(X)，得到初始潮流解Xt(Vt,θt,λt)，其中下标t表示当前计算点；
(6)在上述初始解Xt处求得切向量dXt(dVt,dθt,dλt)，并设置潮流解变化的步长h，得到预测值X't+1(V't+1,θ't+1,λ't+1)，其中下标t+1表示下一个计算点：
X't+1＝Xt+h·dXt；
(7)以X't+1为初始点，重新计算由步骤(1)和步骤(3-3)构建的潮流方程，得到修正值Xt+1，判断Xt+1是否满足步骤(4)中的约束且dλt＞0，若同时满足步骤(4)的约束和dλt＞0的条件，则以Xt+1作为初始解Xt，返回执行步骤(6)，若不满足步骤(4)的约束或不满足dλt＞0的条件，则进一步判断Xt+1是否满足dλt/λt＜ε且dλt＞0，若不满足该条件，则重新调整步长h，如使步长h＝1/2h，返回执行步骤(6)，若满足该条件，输出此时的λ，作为计及电-气耦合系统约束的电压稳定裕度。

Claims (1)

1.一种计及电-气耦合系统约束的电压稳定裕度计算方法，其特征在于该方法包括以下步骤：
(1)建立电-气耦合系统稳态安全运行的等式约束方程，包括：
(1-1)一个电-气耦合系统中电力系统的潮流方程如下：
其中，PGi为电力系统中第i个节点的注入有功功率，PLi为电力系统中第i个节点的注出有功功率，QGi为电力系统中第i个节点的注入无功功率，QLi为电力系统中第i个节点的注出无功功率，Vi、Vj分别为电力系统中第i个节点、第j个节点的电压幅值，θij＝θi-θj，θi、θj分别为电力系统中第i个节点、第j个节点的电压相角，Gij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电导，Bij为与电力系统的节点导纳矩阵Y中第i行、第j列相对应的电纳，电力系统节点导纳矩阵Y从电力系统调度中心获取，Ne为电力系统中所有节点的个数，NPQ为电力系统中给定了有功功率P和无功功率Q的PQ节点的个数；
(1-2)一个电-气耦合系统中天然气系统中管道的水力方程如下：
其中，fkm为天然气系统中第k个节点和第m个节点之间的管道中天然气体积流量，pk，pm分别为第k个节点和第m个节点的压强，Ckm为第k个节点和第m个节点之间管道km的阻力系数，从管道的设计报告中获取，上述天然气系统中管道的水力方程中，当时，上式中的sgnp(pk,pm)＝1，当时，sgnp(pk,pm)＝-1；
(1-3)一个电-气耦合系统中通过燃气轮机耦合的电力系统与天然气系统之间的耦合方程如下：
μG×LG×Hgas＝PG，
其中，LG为燃气轮机的燃气负荷，PG为燃气轮机的有功功率输出，Hgas为天然气的燃烧热值，取值为37.59MJ/m3，μG为燃气轮机的效率系数，由燃气轮机的出厂说明书获取；
(1-4)一个电-气耦合系统中天然气系统的节点气流平衡方程如下：
其中，Lsm为天然气系统中第m个节点的注入体积流量，LLm为天然气系统中第m个节点的注出体积流量；
(2)选择负荷裕度指标λ作为电压稳定裕度的指标，并在下述三种负荷增长方式中选择一种：①其他负荷保持不变，单一负荷的有功功率和无功功率保持原功率因数同时增加；②其他负荷保持不变，选定区域内负荷的有功功率和无功功率保持原功率因数同时增加；③所有负荷的有功功率和无功功率保持原功率因数同时增加；
(3)利用负荷裕度指标λ建立电-气耦合系统的连续潮流模型：
(3-1)建立一个电-气耦合系统中电力系统注入与注出功率变化方程如下：
其中，PLi0为初始时刻节点i的注出有功功率，PGi0为初始时刻节点i的注入有功功率，QLi0为初始时刻节点i的注入无功功率，Ne为电力系统中的节点个数，NPQ为电力系统中PQ节点的个数；
(3-2)建立一个电-气耦合系统中天然气系统中天然气负荷的变化方程如下：
LLm(λ)＝(1+rλ)LLm0，
其中，LLm0为初始时刻第m个节点的注出体积流量，通过天然气系统的运行数据获得；r为电力系统燃气负荷与天然气系统负荷的相关系数，由当地能源统计部门的数据获得；
(3-3)将步骤(3-1)和(3-2)中的连续变化方程代入步骤(1-1)和(1-4)中的方程后得到如下方程：
(4)设定电-气耦合系统稳态安全运行的不等式约束条件，包括：
(4-1)电力系统中发电机组的输出有功功率Pgen大于或等于0，小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-2)电力系统中发电机组的输出无功功率大于或等于该发电机组出厂铭牌上给出的最小功率小于或等于该发电机组出厂铭牌上给出的最大功率即：
(4-3)电力系统第i个节点的电压幅值Ui在设定的电力系统安全运行电压的上限值和下限值U i之间运行，U i取第i个节点额定电压的0.9倍或0.95倍，为第i个节点额定电压的1.1倍或1.05倍，即：
(4-4)天然气系统中第k个节点的压强pk在设定的管道安全运行气压的上限值、下限值p k、内，即：
(4-5)天然气系统中气源供气量Ls大于或等于0，小于或等于该气源能提供天然气流的最大值Ls,max，即：
0≤Ls≤Ls,max；
(5)利用优化法或迭代法求解λ为0时由步骤(1)和步骤(3-3)构建的潮流方程F(X)，得到初始潮流解Xt(Vt,θt,λt)，其中下标t表示当前计算点；
(6)在上述初始解Xt处求得切向量dXt(dVt,dθt,dλt)，并设置潮流解变化的步长h，得到预测值X′t+1(V′t+1,θ′t+1,λ′t+1)，其中下标t+1表示下一个计算点：
X′t+1＝Xt+h·dXt；
(7)以X′t+1为初始点，重新计算由步骤(1)和步骤(3-3)构建的潮流方程，得到修正值Xt+1，判断Xt+1是否满足步骤(4)中的约束且dλt＞0，若同时满足步骤(4)的约束和dλt＞0的条件，则以Xt+1作为初始解Xt，返回执行步骤(6)，若不满足步骤(4)的约束或不满足dλt＞0的条件，则进一步判断Xt+1是否满足dλt/λt＜ε且dλt＞0，若不满足该条件，则重新调整步长h，返回执行步骤(6)，若满足该条件，输出此时的λ，作为计及电-气耦合系统约束的电压稳定裕度。
CN201810335838.4A 2018-04-16 2018-04-16 计及电-气耦合系统约束的电压稳定裕度计算方法 Active CN108667007B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201810335838.4A CN108667007B (zh) 2018-04-16 2018-04-16 计及电-气耦合系统约束的电压稳定裕度计算方法
PCT/CN2018/113635 WO2019200891A1 (zh) 2018-04-16 2018-11-02 计及电-气系统耦合的电力系统电压稳定裕度计算方法
US17/021,445 US20200410145A1 (en) 2018-04-16 2020-09-15 Method for calculating voltage stability margin of power system considering the coupling of electric-gas system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810335838.4A CN108667007B (zh) 2018-04-16 2018-04-16 计及电-气耦合系统约束的电压稳定裕度计算方法

Publications (2)

Publication Number Publication Date
CN108667007A CN108667007A (zh) 2018-10-16
CN108667007B true CN108667007B (zh) 2019-12-13

Family

ID=63783514

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810335838.4A Active CN108667007B (zh) 2018-04-16 2018-04-16 计及电-气耦合系统约束的电压稳定裕度计算方法

Country Status (3)

Country Link
US (1) US20200410145A1 (zh)
CN (1) CN108667007B (zh)
WO (1) WO2019200891A1 (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108667007B (zh) * 2018-04-16 2019-12-13 清华大学 计及电-气耦合系统约束的电压稳定裕度计算方法
CN109378834A (zh) * 2018-11-01 2019-02-22 三峡大学 基于信息最大相关的大规模电网电压稳定裕度评估系统
CN109242366B (zh) * 2018-11-06 2020-08-07 国电南瑞科技股份有限公司 一种电-气互联综合能源系统的多时段潮流优化方法
CN110070213B (zh) * 2019-03-28 2022-05-20 广东工业大学 一种电-气综合能源系统的日前调度方法
CN110277785B (zh) * 2019-06-26 2021-08-03 国网浙江省电力有限公司电力科学研究院 基于连续多能流的电气耦合系统负荷裕度计算方法及系统
CN111125880B (zh) * 2019-11-25 2022-07-22 国网四川省电力公司电力科学研究院 一种暂态稳定视角下电力系统仿真数据生成方法
CN111241479B (zh) * 2020-01-10 2024-03-01 河海大学 基于交叉熵及客观熵权法的电-热互联综合能源系统风险评估方法
CN111769603B (zh) * 2020-07-13 2022-04-08 国网天津市电力公司 一种基于电-气耦合系统安全裕度的机组优化调度方法
CN112101726B (zh) * 2020-08-14 2022-10-21 国网宁夏电力有限公司中卫供电公司 天然气网络受电压暂降影响的评估方法、介质及系统
CN112528480B (zh) * 2020-12-02 2022-08-12 清华大学 一种提升负荷裕度的热电耦合系统灾前预防方法
CN112861292B (zh) * 2021-01-12 2022-08-05 浙江大学 一种电-气综合能源系统恢复改善方法
CN112928753B (zh) * 2021-02-03 2023-01-06 东北电力大学 一种多能协同的配电网主动解列控制方法
CN112950031A (zh) * 2021-03-03 2021-06-11 重庆大学 一种天然气管网稳态能流计算的改进迭代法
CN113326605B (zh) * 2021-05-08 2022-07-26 华南理工大学 考虑灵活性冷负荷调控的多模式集中供冷系统优化方法
CN113221358B (zh) * 2021-05-13 2022-09-27 浙江大学 基于可靠性参数的电-气耦合系统的备用出力优化方法
CN114060825B (zh) * 2021-11-15 2024-01-26 南方电网电力科技股份有限公司 一种垃圾焚烧炉中多污染物协同处置方法及相关装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103336882A (zh) * 2013-05-23 2013-10-02 国家电网公司 一种基于时域仿真的全过程动态电压稳定裕度评估方法
CN106056478A (zh) * 2016-06-12 2016-10-26 清华大学 一种电‑热耦合系统中热网的区间潮流计算方法
CN106096269A (zh) * 2016-06-12 2016-11-09 清华大学 一种电‑气耦合系统中天然气网的区间潮流计算方法
US9921602B2 (en) * 2013-05-14 2018-03-20 Rensselaer Polytechnic Institute Methods of computing steady-state voltage stability margins of power systems

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE60144367D1 (de) * 2001-05-21 2011-05-19 Abb Research Ltd Stabilität-Vorhersage für elektrisches Energieversorgungsnetz
US7603203B2 (en) * 2006-10-09 2009-10-13 Electric Power Research Institute, Inc. Method for voltage instability load shedding using local measurements
JP6412822B2 (ja) * 2015-04-22 2018-10-24 株式会社日立製作所 電力系統電圧無効電力監視制御装置及び方法
CN108667007B (zh) * 2018-04-16 2019-12-13 清华大学 计及电-气耦合系统约束的电压稳定裕度计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9921602B2 (en) * 2013-05-14 2018-03-20 Rensselaer Polytechnic Institute Methods of computing steady-state voltage stability margins of power systems
CN103336882A (zh) * 2013-05-23 2013-10-02 国家电网公司 一种基于时域仿真的全过程动态电压稳定裕度评估方法
CN106056478A (zh) * 2016-06-12 2016-10-26 清华大学 一种电‑热耦合系统中热网的区间潮流计算方法
CN106096269A (zh) * 2016-06-12 2016-11-09 清华大学 一种电‑气耦合系统中天然气网的区间潮流计算方法

Also Published As

Publication number Publication date
US20200410145A1 (en) 2020-12-31
WO2019200891A1 (zh) 2019-10-24
CN108667007A (zh) 2018-10-16

Similar Documents

Publication Publication Date Title
CN108667007B (zh) 计及电-气耦合系统约束的电压稳定裕度计算方法
CN104734147B (zh) 一种综合能源系统概率能量流分析方法
CN108846507B (zh) 基于混合整数二阶锥规划的电-气耦合系统日前经济调度方法
CN107665384B (zh) 一种含多区域能源站的电力-热力综合能源系统调度方法
CN106096269B (zh) 一种电-气耦合系统中天然气网的区间潮流计算方法
CN109031952B (zh) 一种电-气互联综合能源系统混合控制方法
CN113129164B (zh) 一种天然气管网天然气流量压力调度决策指标的计算方法
CN104935017B (zh) 基于改进轻鲁棒优化模型的风电与火电机组组合方法
CN105958531B (zh) 一种电-气耦合多能流网络状态估计方法
CN110232217A (zh) 一种综合能源配电系统运行域建模方法
CN104242314A (zh) 一种风电场汇集区域电网电压安全域值的确定方法
WO2023015923A1 (zh) 一种基于序列凸规划的电气互联系统最优能流计算方法
CN104899659A (zh) 一种智能园区多级能耗传递模型
CN103500997B (zh) 基于混合多目标λ迭代法和牛顿法的电力系统调度方法
CN115099567A (zh) 一种电气互联综合能源系统优化调度方法及终端
CN111241479A (zh) 基于交叉熵及客观熵权法的电-热互联综合能源系统风险评估方法
Mu et al. Research on sensitivity analysis of wind power consumption capability of integrated energy system based on unified optimal power flow model
CN112862157A (zh) 一种考虑耦合度的综合能源系统优化调度方法及系统
Dai et al. Modeling and analysis of electrical heating system based on entransy dissipation-based thermal resistance theory
Aly et al. Voltage stability assessment for radial distribution power system with wind power penetration
CN103544654A (zh) 电网经济调度局部最小解确定和全局最小解搜索方法
CN115545340A (zh) 一种考虑双重不确定性的综合能源站-网协同规划方法
Liu et al. A method to calculate energy station’s output in community energy planning considering the attenuation and delay of pipe network
Chen et al. Optimization of two-stage closure law of turbine wicket gates and its application
CN105160063B (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