CN114970154A - 一种暂态电压支撑能力量化评估指标构建方法 - Google Patents

一种暂态电压支撑能力量化评估指标构建方法 Download PDF

Info

Publication number
CN114970154A
CN114970154A CN202210582096.1A CN202210582096A CN114970154A CN 114970154 A CN114970154 A CN 114970154A CN 202210582096 A CN202210582096 A CN 202210582096A CN 114970154 A CN114970154 A CN 114970154A
Authority
CN
China
Prior art keywords
subsystem
input
power system
evaluation index
quantitative evaluation
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
CN202210582096.1A
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.)
Xian Jiaotong University
Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd
Original Assignee
Xian Jiaotong University
Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd
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 Xian Jiaotong University, Electric Power Research Institute of State Grid Shaanxi Electric Power Co Ltd filed Critical Xian Jiaotong University
Priority to CN202210582096.1A priority Critical patent/CN114970154A/zh
Publication of CN114970154A publication Critical patent/CN114970154A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • 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
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及电力系统技术领域，公开了一种暂态电压支撑能力量化评估指标构建方法，通过电力系统的各个子系统的数学模型获取状态变量x和输出量y，并根据各个子系统的状态变量x和输出量y计算得到各个子系统的属性；根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当满足稳定性时结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程的安全状态，本发明将互联系统中子系统的稳定分析解耦，通过两个代数不等式分析系统的稳定性。该方法计算量小，计算效率高。

Description

一种暂态电压支撑能力量化评估指标构建方法
技术领域
本发明涉及电力系统技术领域，具体为一种暂态电压支撑能力量化评估指标构建方法。
背景技术
在当前能源短缺、环境污染和气候变化等问题日趋严重的背景下，大规模开发利用风能、太阳能等清洁可再生的新能源成为世界各国发展的共识，新能源的投入使用对于我国能源结构的调整具有非常重要的推动作用，也成为我国电力系统未来的重要发展方向。
然而，由于新能源发电设备采用电力电子接口并网，其固有的弱阻尼和低惯性特点会恶化系统的动态特性，引发电力系统的电压崩溃，该问题严重制约了可再生能源在电力系统中的大规模应用。因此，亟需发展一种灵活高效的暂态电压稳定性分析方法，准确评估高比例新能源接入场景下电网的暂态电压支撑能力，为电力的系统安全稳定运行提供指导。
现有电力系统暂态电压稳定性分析方法主要包括时域仿真法与直接法。
由于时域仿真法可以考虑复杂的系统模型，能够获得准确的稳定性分析结果，该方法得到了广泛应用。时域仿真法存在计算量大、难以提供系统稳定裕度的定量信息等不足，不适用于新能源接入场景下电网的安全稳定分析。
直接法在新能源接入场景下电网安全稳定分析中已有应用，通过构造暂态能量函数，定量地分析系统的暂态电压稳定性。由于新能源接入场景下电网具有强非线性和不确定性等特征，构造Lyapunov函数或能量函数并精确刻画电网安全稳定边界面临挑战，且直接法难以量化外部扰动对电网安全稳定运行的影响。
发明内容
针对现有技术中存在的问题，本发明提供一种暂态电压支撑能力量化评估指标构建方法。
本发明是通过以下技术方案来实现：
一种暂态电压支撑能力量化评估指标构建方法，包括如下步骤：
步骤1，根据电力系统的各个子系统的数学模型获取状态变量x和输出量y；
步骤2，通过所获取的状态变量和输出量计算得到各个子系统的属性；
步骤3，根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当电力系统受扰后无法回到稳定状态，则电力系统失稳；当电力系统受扰后回到稳定状态，则执行步骤4；
步骤4，结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程是否满足安全稳定约束条件；当电力系统暂态过程不满足安全稳定约束条件，则电力系统运行在不安全状态；反之，则电力系统运行在安全状态。
优选的，步骤1中，各个子系统包括发电机、分布式电源、电动机、恒阻抗负载和恒功率负载；其中发电机和分布式电源选取机端电压的标幺值作为输出；电动机、恒阻抗负载和恒功率负载选取端电流的标幺值作为输出。
优选的，步骤2中，所获取的状态变量和输出量的计算公式如下：
Figure BDA0003664350770000021
y＝h(x，u)；
其中，
Figure BDA0003664350770000022
为，输入u∈U∈Rm，函数f：D→Rn，g：D→Rn×m；f和g关于x连续且满足局部李普希茨条件；f(0，0)＝0，h(0，0)＝0，D和U分别表示状态变量和外部输入的局部域。
进一步的，所获取的状态变量和输出量的计算公式计算得到各个子系统的属性，包括如下步骤：
对含有外部输入的非线性系统，存在
Figure BDA0003664350770000031
则对于任何初始状态和外部输入，若以下不等式成立：
Figure BDA0003664350770000032
其中，x(t)表示系统状态变量随时间的变化情况，x0表示系统状态变量初始值，t表示时间，γ函数表示扰动对系统的影响，|·|表示Euclidean范数，||·||表示矩阵的范数，β函数则可以量化系统状态变量随时间变化的动态过程，||u||∞是使得|u(t)|≤a对于所有的时间t均成立的最小的a；
Figure BDA0003664350770000033
为比较函数；
根据积分-积分估计，定义公式如下：
Figure BDA0003664350770000034
Figure BDA0003664350770000035
其中，α为增益函数，α0为增益函数初始值，s为拉普拉斯算子；
计算相应的局部输入范围、初值范围以及增益函数α、α0、γ时，各个状态变量、输入、输出需要以相应的平衡点作为参考，定义公式如下：
Figure BDA0003664350770000036
Figure BDA0003664350770000037
Figure BDA0003664350770000038
其中，ue、xe、ye分别为子系统输入、状态变量以及输出的平衡点；
给定输入信号u，固定系统的初始状态x0，得到输出y；
固定|x0|＝0，此时α0(0)＝0，则近似估计输入-输出增益的计算公如下：
Figure BDA0003664350770000041
其中，
Figure BDA0003664350770000042
为子系统输出积分能量，
Figure BDA0003664350770000043
为子系统输入积分能量；
改变输入的大小或采用不同形式的输入，计算输入-输出增益可以得到一系列的γ，从中取最大的γmax近似子系统的输入-输出增益；若系统状态变量和输出的稳态运行区间型积分能量值超过上限值，则电力系统失稳；逐步减小输入直至电力系统稳定，估计子系统输入信号的范围。
更进一步的，确定γmax后，在规定的范围内改变系统的初值以及输入，计算公司如下：
Figure BDA0003664350770000044
其中，γmax为输入-输出增益最大值。
优选的，步骤3，基于得到的子系统LISS/LIOS属性，考虑由n个子系统组成的动态系统，第i个子系统的数学模型表达式如下：
Figure BDA0003664350770000045
yi＝hi(xi，ui，ωi)
根据各个子系统的属性基于电气网络连接关系所组成的数学模型如下：
Figure BDA0003664350770000046
Figure BDA0003664350770000047
0＝g(y，u)；
其中，x＝[x1 … xn]T∈RN，N＝n1+…+nn，u＝[u1 … un]T∈Rm，m＝m1+…+mn，y＝[y1…yn]T∈Rp，p＝l1+…+ln；
Figure BDA0003664350770000048
是第i个子系统的状态变量，
Figure BDA0003664350770000049
和
Figure BDA00036643507700000410
分别是子系统的输入与输出，ωi是子系统受到的外部扰动输入。
优选的，步骤3中，判断电力系统受扰后能否回到稳定状态所满足的条件如下：
A，当|xoi|≤vi，||ui||∞≤τi，||ωi||∞≤εi时，每个子系统是LISS和LIOS的，并且具有线性的渐近增益；
B，函数g(y，u)满足隐函数定理，在电力系统的应用中，该条件表示系统潮流方程有解；存在zij≥0，di≥0，使得下式满足
|ui(t)|≤∑zij(|yj(t)|)+di
C，小增益条件满足，即
ρ(GIOS)＜1
其中，GIOS＝ΓIOSZ，ΓIOS是形如
Figure BDA0003664350770000051
的输入/输出增益矩阵；ρ表示矩阵GIOS的谱半径。
进一步的，提炼出反映互联系统稳定运行能力的量化评估指标如下：
a＝1-ρ(GLIOS)；
其中，a为定义的稳定性量化评估指标。
优选的，步骤4中，判断电力系统暂态过程是否满足安全稳定约束条件如下：
(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d≤τ；
其中，τ＝[τ1 … τn]T，d＝[d1 … dn]T，
Figure BDA0003664350770000052
其中，τi为第i个子系统输入上限。
优选的，计算安全性量化评估指标公式如下：
Figure BDA0003664350770000053
其中，b＝Z(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d，bi和τi分别为b和τ中的第i个元素；当安全稳定约束条件不满足，即Z·(Id-GIOS)-1·βIOS((|xo|)c,0)>τ时，系统运行在不安全状态。
与现有技术相比，本发明具有以下有益的技术效果：
本发明提供了一种暂态电压支撑能力量化评估指标构建方法，通过电力系统的各个子系统的数学模型获取状态变量x和输出量y，并根据各个子系统的状态变量x和输出量y计算得到各个子系统的属性；根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当满足稳定性时结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程的安全状态，本发明将互联系统中子系统的稳定分析解耦，通过两个代数不等式分析系统的稳定性。该方法计算量小，计算效率高。
进一步的，本发明所提暂态电压稳定性分析方法能量化评估高比例新能源接入场景下电网的暂态电压稳定性，为衡量新能源的接入规模提供有效途径；所提暂态电压支撑能力量化评估指标统筹考虑系统稳定性和安全性，能够全面准确地反映高比例新能源接入场景下电网的暂态电压稳定水平。
附图说明
图1为本发明中暂态电压支撑能力量化评估指标构建方法的流程图；
图2为输入状态稳定理论示意图；
图3为实施例中电力系统地理接线图；
图4为系统故障时不同新能源渗透率下C区域750kV母线电压动态变化示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述：
参见图1，本发明一个实施例中，提供了一种暂态电压支撑能力量化评估指标构建方法，包括如下步骤：
步骤1，根据电力系统的各个子系统的数学模型获取状态变量x和输出量y；
具体的，各个子系统包括发电机、分布式电源、电动机、恒阻抗负载和恒功率负载；其中发电机和分布式电源选取机端电压的标幺值作为输出；电动机、恒阻抗负载和恒功率负载选取端电流的标幺值作为输出。
步骤2，通过所获取的状态变量和输出量计算得到各个子系统的属性；
具体的，输入、输出以及系统状态积分形式的衡量与能量具有直接关系，因此，在介绍子系统积分能量的积分形式的LISS/LIOS仿真分析方法之前给出积分能量定义，其目的主要在以下两个方面：1)便于估计子系统的LISS稳定域、渐近增益等属性。2)便于判断系统状态、输出变量是否在安全运行范围内。
首先介绍输入状态稳定理论(Input-to-state Stability，ISS)理论。不失一般性，考虑含有外部输入的非线性系统：
Figure BDA0003664350770000071
y＝h(x，u)；
其中，
Figure BDA0003664350770000072
为，输入u∈U∈Rm，函数f：D→Rn，g：D→Rn×m；f和g关于x连续且满足局部李普希茨条件；f(0，0)＝0，h(0，0)＝0，D和U分别表示状态变量和外部输入的局部域。
所获取的状态变量和输出量的计算公式计算得到各个子系统的属性，包括如下步骤：
对含有外部输入的非线性系统，存在
Figure BDA0003664350770000081
则对于任何初始状态和外部输入，若以下不等式成立：
Figure BDA0003664350770000082
则该含有外部输入的非线性系统是ISS的；
其中，x(t)表示系统状态变量随时间的变化情况…，x0表示系统状态变量初始值，t表示时间，γ函数表示扰动对系统的影响，|·|表示Euclidean范数，||·||表示矩阵的范数，β函数则可以量化系统状态变量随时间变化的动态过程，||u||∞是使得|u(t)|≤a对于所有的时间t均成立的最小的a；
Figure BDA0003664350770000083
为比较函数；
以下给出比较函数的定义：一个函数γ：R≥0→R≥0，若满足连续、严格单增以及γ(0)＝0，则称之为
Figure BDA00036643507700000811
函数；进一步，若γ函数满足当s→∞时，γ(s)→∞，则称之为
Figure BDA0003664350770000084
函数；一个函数β：R≥0×R≥0→R≥0，若满足连续条件，且对于任一固定的t≥0，函数β(·，t)是
Figure BDA0003664350770000085
函数，β(s，t)关于t是递减的，并且当t→∞时，β(s，t)→0，则称β(s，t)为
Figure BDA0003664350770000086
函数。
ISS范数描述表达式的右边由两项组成，一项是与初值相关的
Figure BDA0003664350770000087
函数，该项随着时间趋于无穷衰减至0；另一项是描述终值的
Figure BDA0003664350770000088
函数，式(2)中的γ称为输入-状态增益。
Figure BDA0003664350770000089
函数描述了系统状态到其吸引子的吸引速度，
Figure BDA00036643507700000810
函数则描述了此吸引子的半径。
ISS几何意义如图2所示，初值为x0的系统，其状态轨迹以一定的收敛速度，最终收敛至一个球的内部，该球的半径与输入||u||∞的无穷范数成“正比”。实际工程应用中绝大多数系统仅是局部稳定的，以下介绍局部输入状态稳定理论(Local Input-to-stateStability，LISS)的定义。
对于传统技术而言，对于任意x0∈Ω∈Rn，u∈U∈Rm，如果存在比较函数
Figure BDA0003664350770000091
则对于任何初始状态和外部输入，若以下不等式成立：
Figure BDA0003664350770000092
则上述非线性系统是LISS的。
其中，Ω和U分别表示初始状态和外部输入的LISS域。
而本发明根据积分-积分估计，定义公式如下：
Figure BDA0003664350770000093
Figure BDA0003664350770000094
其中，α为增益函数，α0为增益函数初始值，s为拉普拉斯算子。
计算相应的局部输入范围、初值范围以及增益函数α、α0、γ时，各个状态变量、输入、输出需要以相应的平衡点作为参考，定义公式如下：
Figure BDA0003664350770000095
Figure BDA0003664350770000096
Figure BDA0003664350770000097
其中，ue、xe、ye分别为子系统输入、状态变量以及输出的平衡点；
给定输入信号u，固定系统的初始状态x0，得到输出y；
固定|x0|＝0，此时α0(0)＝0，则近似估计输入-输出增益的计算公如下：
Figure BDA0003664350770000101
其中，
Figure BDA0003664350770000102
为子系统输出积分能量，
Figure BDA0003664350770000103
为子系统输入积分能量。
改变输入的大小或采用不同形式的输入，计算输入-输出增益可以得到一系列的γ，从中取最大的γmax近似子系统的输入-输出增益；若系统状态变量和输出的稳态运行区间型积分能量值超过上限值，则电力系统失稳；逐步减小输入直至电力系统稳定，估计子系统输入信号的范围。
确定γmax后，在规定的范围内改变系统的初值以及输入，计算公式如下：
Figure BDA0003664350770000104
若改变初值，则可以得到不同的α0函数。
其中，γmax为输入-输出增益最大值。
步骤3，根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当电力系统受扰后无法回到稳定状态，则电力系统失稳；当电力系统受扰后回到稳定状态，则执行步骤4；
具体的，考虑由n个子系统组成的动态系统，第i个子系统的数学模型表达式如下：
Figure BDA0003664350770000105
yi＝hi(xi，ui，ωi)
其中，
Figure BDA0003664350770000106
是第i个子系统的状态变量，
Figure BDA0003664350770000107
和
Figure BDA0003664350770000108
分别是子系统的输入与输出，ωi是子系统受到的外部扰动输入。
由这些子系统互联而成的系统数学模型如下：
Figure BDA0003664350770000111
Figure BDA0003664350770000112
0＝g(y，u)
其中，x＝[x1 … xn]T∈RN，N＝n1+…+nn，u＝[u1 … un]T∈Rm，m＝m1+…+mn，y＝[y1… yn]T∈Rp，p＝l1+…+ln。
具体的，判断电力系统受扰后能否回到稳定状态所满足的条件如下：
A，当|xoi|≤vi，||ui||∞≤τi，||ωi||∞≤εi时，每个子系统是LISS和LIOS的，并且具有线性的渐近增益；
B，函数g(y，u)满足隐函数定理，在电力系统的应用中，该条件表示系统潮流方程有解；存在zij≥0，di≥0，使得下式满足
|ui(t)|≤∑zij(|yj(t)|)+di
C，小增益条件满足，即
ρ(GIOS)＜1
其中，GIOS＝ΓIOSZ，ΓIOS是形如
Figure BDA0003664350770000113
的输入/输出增益矩阵；ρ表示矩阵GIOS的谱半径。
具体的，提炼出反映互联系统稳定运行能力的量化评估指标如下：
a＝1-ρ(GLIOS)；
其中，a为定义的稳定性量化评估指标。
步骤4，结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程是否满足安全稳定约束条件；当电力系统暂态过程不满足安全稳定约束条件，则电力系统运行在不安全状态；反之，则电力系统运行在安全状态。
具体的，判断电力系统暂态过程是否满足安全稳定约束条件如下：
(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d≤τ；
其中，τ＝[τ1 … τn]T，d＝[d1 … dn]T，
Figure BDA0003664350770000121
其中，τi为第i个子系统输入上限。
具体的，计算安全性量化评估指标公式如下：
Figure BDA0003664350770000122
其中，b＝Z(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d，bi和τi分别为b和τ中的第i个元素；当安全稳定约束条件不满足，即Z·(Id-GIOS)-1·βIOS((|xo|)c，0)＞τ时，系统运行在不安全状态。
上述稳定性和安全性量化评估指标能够充分反映电力系统的安全稳定运行水平，量化评估高比例新能源接入场景下系统的暂态支撑能力，为分析新能源接入对电力系统暂态电压稳定性的影响提供理论依据和评价标准。
下面通过一个实例对本发明实施例进行进一步说明，以下仅为本发明实施例的一个实例，本发明实施例并不以此为限。
某地区典型的750kV电网结构一如图3所示，传统电源总装机为1175万千瓦，风电总装机985万千瓦，光伏总装机1292万千瓦，负荷水平为3130万千瓦。其中，归算到750kV下各母线的新能源装机分别为：D区域750kV供电区风电装机为636万千瓦，光伏装机为259万千瓦；A区域750kV供电区风电装机为125万千瓦，光伏装机为127万千瓦；C区域750kV供电区风电装机为27万千瓦，光伏装机为75万千瓦；F区域750kV供电区风电装机为115万千瓦，光伏装机为121万千瓦；J区域750kV供电区风电装机为21万千瓦，光伏装机为132万千瓦I区域750kV供电区光伏装机为71万千瓦；H区域750kV供电区风电装机为21万千瓦，光伏装机为108万千瓦；G区域750kV供电区风电装机为33万千瓦，光伏装机为15万千瓦。此时，某地区电网新能源渗透率为35％。
研究系统发生短路故障时不同新能源渗透率对电压动态过程的影响。三相短路故障发生在陕榆横750kV母线上。基于提出的子系统LIOS属性估计方法和互联系统稳定性分析思路，对简化的某地区电力系统暂态电压稳定性进行分析。具体地，设置新能源渗透率分别为30％、50％和70％，判断系统在发生三相短路故障后是否满足小增益条件：ρ(GLIOS)＜1和安全稳定约束条件：Z·(Id-GIOS)-1·βIOS((|xo|)c，0)＜τ。若小增益条件不满足，系统失稳；若小增益条件满足，计算稳定性量化评估指标：1-ρ(GLIOS)，该指标越大，系统的稳定裕度越大，稳定性越好。若安全稳定约束条件不满足，系统运行在不安全状态；若安全稳定约束条件满足，计算安全性量化评估指标
Figure BDA0003664350770000131
该指标越大，系统的安全运行裕度越大，安全性越好。计算结果如表1所示。
Figure BDA0003664350770000132
表1不同新能源渗透率下系统量化评估指标
该结果表明：当系统发生短路故障时，随着新能源渗透率增加，小增益条件和安全稳定约束条件始终满足，但系统稳定性量化评估指标和安全性量化评估指标逐渐减小，系统的稳定裕度和安全裕度不断降低。时域仿真验证了分析结果，C区域750kV母线电压动态变化如图4所示。
综上所述，本发明提供了一种暂态电压支撑能力量化评估指标构建方法，通过电力系统的各个子系统的数学模型获取状态变量x和输出量y，并根据各个子系统的状态变量x和输出量y计算得到各个子系统的属性；根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当满足稳定性时结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程的安全状态，本发明将互联系统中子系统的稳定分析解耦，通过两个代数不等式分析系统的稳定性。该方法计算量小，计算效率高。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (10)

1.一种暂态电压支撑能力量化评估指标构建方法，其特征在于，包括如下步骤：
步骤1，根据电力系统的各个子系统的数学模型获取状态变量x和输出量y；
步骤2，通过所获取的状态变量和输出量计算得到各个子系统的属性；
步骤3，根据各个子系统的属性基于电气网络连接关系，计算电力系统稳定性量化评估指标，判断电力系统受扰后能否回到稳定状态，当电力系统受扰后无法回到稳定状态，则电力系统失稳；当电力系统受扰后回到稳定状态，则执行步骤4；
步骤4，结合电力系统初始运行状态，计算安全性量化评估指标，判断电力系统暂态过程是否满足安全稳定约束条件；当电力系统暂态过程不满足安全稳定约束条件，则电力系统运行在不安全状态；反之，则电力系统运行在安全状态。
2.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，步骤1中，各个子系统包括发电机、分布式电源、电动机、恒阻抗负载和恒功率负载；其中发电机和分布式电源选取机端电压的标幺值作为输出；电动机、恒阻抗负载和恒功率负载选取端电流的标幺值作为输出。
3.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，步骤2中，所获取的状态变量和输出量的计算公式如下：
Figure FDA0003664350760000012
y＝h(x,u)；
其中，
Figure FDA0003664350760000011
为,输入u∈U∈Rm,函数f:D→Rn,g:D→Rn×m；f和g关于x连续且满足局部李普希茨条件；f(0,0)＝0，h(0,0)＝0，D和U分别表示状态变量和外部输入的局部域。
4.根据权利要求3所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，所获取的状态变量和输出量的计算公式计算得到各个子系统的属性，包括如下步骤：
对含有外部输入的非线性系统，存在
Figure FDA0003664350760000021
则对于任何初始状态和外部输入，若以下不等式成立：
Figure FDA0003664350760000022
其中，x(t)表示系统状态变量随时间的变化情况，x0表示系统状态变量初始值，t表示时间，γ函数表示扰动对系统的影响，|·|表示Euclidean范数，||·||表示矩阵的范数，β函数则可以量化系统状态变量随时间变化的动态过程，||u||∞是使得|u(t)|≤a对于所有的时间t均成立的最小的a；
Figure FDA0003664350760000023
为比较函数；
根据积分-积分估计，定义公式如下：
∫0 tα(|x(s)|)ds≤α0(|x0|)+∫0 tγ(|u(s)|)ds
Figure FDA0003664350760000024
其中，α为增益函数，α0为增益函数初始值，s为拉普拉斯算子；
计算相应的局部输入范围、初值范围以及增益函数α、α0、γ时，各个状态变量、输入、输出需要以相应的平衡点作为参考，定义公式如下：
Figure FDA0003664350760000025
Figure FDA0003664350760000026
Figure FDA0003664350760000027
其中，ue、xe、ye分别为子系统输入、状态变量以及输出的平衡点；
给定输入信号u，固定系统的初始状态x0，得到输出y；
固定|x0|＝0，此时α0(0)＝0，则近似估计输入-输出增益的计算公如下：
Figure FDA0003664350760000028
其中，
Figure FDA0003664350760000031
为子系统输出积分能量，
Figure FDA0003664350760000032
为子系统输入积分能量；
改变输入的大小或采用不同形式的输入，计算输入-输出增益可以得到一系列的γ，从中取最大的γmax近似子系统的输入-输出增益；若系统状态变量和输出的稳态运行区间型积分能量值超过上限值，则电力系统失稳；逐步减小输入直至电力系统稳定，估计子系统输入信号的范围。
5.根据权利要求4所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，确定γmax后，在规定的范围内改变系统的初值以及输入，计算公司如下：
Figure FDA0003664350760000033
其中，γmax为输入-输出增益最大值。
6.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，步骤3，基于得到的子系统LISS/LIOS属性，考虑由n个子系统组成的动态系统，第i个子系统的数学模型表达式如下：
Figure FDA0003664350760000034
yi＝hi(xi，ui，ωi)
根据各个子系统的属性基于电气网络连接关系所组成的数学模型如下：
Figure FDA0003664350760000035
Figure FDA0003664350760000036
0＝g(y，u)；
其中，x＝[x1 … xn]T∈RN，N＝n1+…+nn，u＝[u1 … un]T∈Rm，m＝m1+…+mn，y＝[y1 …yn]T∈Rp，p＝l1+…+ln；
Figure FDA0003664350760000037
是第i个子系统的状态变量，
Figure FDA0003664350760000038
和
Figure FDA0003664350760000039
分别是子系统的输入与输出，ωi是子系统受到的外部扰动输入。
7.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，步骤3中，判断电力系统受扰后能否回到稳定状态所满足的条件如下：
A，当|xoi|≤vi，||ui||∞≤τi，||ωi||∞≤εi时，每个子系统是LISS和LIOS的，并且具有线性的渐近增益；
B，函数g(y，u)满足隐函数定理，在电力系统的应用中，该条件表示系统潮流方程有解；存在zij≥0，di≥0，使得下式满足
|ui(t)|≤∑zij(|yj(t)|)+di
C，小增益条件满足，即
ρ(GIOS)＜1
其中，GIOS＝ΓIOSZ，ΓIOS是形如
Figure FDA0003664350760000041
的输入/输出增益矩阵；ρ表示矩阵GIOS的谱半径。
8.根据权利要求7所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，提炼出反映互联系统稳定运行能力的量化评估指标如下：
a＝1-ρ(GLIos)；
其中，a为定义的稳定性量化评估指标。
9.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，步骤4中，判断电力系统暂态过程是否满足安全稳定约束条件如下：
(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d≤τ；
其中，τ＝[τ1 … τn]T，d＝[d1 … dn]T，
Figure FDA0003664350760000042
其中，τi为第i个子系统输入上限。
10.根据权利要求1所述的一种暂态电压支撑能力量化评估指标构建方法，其特征在于，计算安全性量化评估指标公式如下：
Figure FDA0003664350760000051
其中，b＝Z(Id-GIOS)-1(βIOS((|xo|)c，0)+ΓIOS·d)+d，bi和τi分别为b和τ中的第i个元素；当安全稳定约束条件不满足，即Z·(Id-GIOS)-1·βIOS((|xo|)c，0)＞τ时，系统运行在不安全状态。
CN202210582096.1A 2022-05-26 2022-05-26 一种暂态电压支撑能力量化评估指标构建方法 Pending CN114970154A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210582096.1A CN114970154A (zh) 2022-05-26 2022-05-26 一种暂态电压支撑能力量化评估指标构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210582096.1A CN114970154A (zh) 2022-05-26 2022-05-26 一种暂态电压支撑能力量化评估指标构建方法

Publications (1)

Publication Number Publication Date
CN114970154A true CN114970154A (zh) 2022-08-30

Family

ID=82955691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210582096.1A Pending CN114970154A (zh) 2022-05-26 2022-05-26 一种暂态电压支撑能力量化评估指标构建方法

Country Status (1)

Country Link
CN (1) CN114970154A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115588983A (zh) * 2022-10-25 2023-01-10 广东工业大学 一种暂态电压稳定评估方法
CN117215205A (zh) * 2023-11-09 2023-12-12 国网经济技术研究院有限公司 基于决策树和iss理论的直流系统控制参数分析方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115588983A (zh) * 2022-10-25 2023-01-10 广东工业大学 一种暂态电压稳定评估方法
CN115588983B (zh) * 2022-10-25 2023-04-18 广东工业大学 一种暂态电压稳定评估方法
CN117215205A (zh) * 2023-11-09 2023-12-12 国网经济技术研究院有限公司 基于决策树和iss理论的直流系统控制参数分析方法
CN117215205B (zh) * 2023-11-09 2024-02-06 国网经济技术研究院有限公司 基于决策树和iss理论的直流系统控制参数分析方法

Similar Documents

Publication Publication Date Title
Liu et al. Online voltage stability assessment for load areas based on the holomorphic embedding method
Ju et al. Analytical assessment for transient stability under stochastic continuous disturbances
CN114970154A (zh) 一种暂态电压支撑能力量化评估指标构建方法
Mohn et al. Tracing PV and QV curves with the help of a CRIC continuation method
CN108667005B (zh) 一种计及新能源影响的电网静动态结合脆弱性评估方法
Hou et al. Cluster computing-based trajectory sensitivity analysis application to the WECC system
CN109474014B (zh) 一种对双馈风场接入电网友好性的量化评估方法
CN108565852B (zh) 一种三阶段递进的大电网电压稳定评估的故障筛选与排序方法
CN104766142B (zh) 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法
Bidram et al. Capacitor design considering first swing stability of distributed generations
CN112907075A (zh) 一种电力系统综合负荷模型参数辨识方法
Ejebe et al. Online dynamic security assessment in an EMS
Djari et al. Voltage stability assessment using the VQ sensitivity and modal analyses methods
Ma et al. Impact analysis of wind generation on voltage stability and system load margin
CN116742632A (zh) 一种配电网最大承载力估算方法及系统
CN112510685B (zh) 一种城市电网风险规避预测控制方法及装置
CN111555311B (zh) 一种电力系统即插即用的稳定性分析与控制方法
Zhou et al. A novel approach to trace time-domain trajectories of power systems in multiple time scales
CN108448569B (zh) 一种大电网机电暂态仿真数值收敛性的识别和提高方法
Cao et al. Polynomial Approximation of Transient Voltage Stability Region Boundary in the Parameter Space Based on Collocation Method
CN113852123B (zh) 一种电力系统自动电压控制方法、装置、电子设备和存储介质
Li et al. Systematic Evaluation of Deep Neural Network Based Dynamic Modeling Method for AC Power Electronic System
Zhang et al. An uncertain harmonic power-flow algorithm based on cloud model for energy internet
CN117215205B (zh) 基于决策树和iss理论的直流系统控制参数分析方法
CN118297473A (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