CN107658881A - 基于戴维南等值方法的电压稳定临界点判断方法 - Google Patents
基于戴维南等值方法的电压稳定临界点判断方法 Download PDFInfo
- Publication number
- CN107658881A CN107658881A CN201710889796.4A CN201710889796A CN107658881A CN 107658881 A CN107658881 A CN 107658881A CN 201710889796 A CN201710889796 A CN 201710889796A CN 107658881 A CN107658881 A CN 107658881A
- Authority
- CN
- China
- Prior art keywords
- thevenin
- voltage
- equivalent
- load
- impedance
- 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
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
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
-
- 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]
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于戴维南等值方法的电压稳定临界点判断方法，其步骤为：S1：将***等值为一个两节点戴维南等值模型；S2：根据戴维南等值模型，得出负荷节点的量测数据和戴维南等值参数，建立量测矩阵方程；S3：对戴维南等值参数进行求解：根据量测矩阵方程，状态方程，利用递归最小二乘法，求解戴维南等值参数的状态量的递归表达式，更新状态量；S4：根据状态方程得出残差方程，利用线性估计状态量的不确定性，对量测数据的虚假数据进行检测处理；S5：根据等阻抗模判据，建立评价指标，并对电压稳定临界点进行判断。本发明具有良好的精度，且适用于大扰动下暂态过程中戴维南等值参数的计算；快速有效地在线实时跟踪戴维南等值参数。
Description
技术领域
本发明涉及一种电力***电压稳定分析与控制领域，特别涉及一种基于戴维南等值方法的电压稳定临界点判断方法。
背景技术
随着电力需求的快速增长，电网之间互联进程的加快，以及特高压交直流工程的相继投运，电力***的运行方式更加多变复杂，并且逐渐向电压稳定运行边界靠近。
在电力***电压稳定分析与控制中，电压崩溃点的计算具有极其重要的作用。给定一个基态的电力***，并给出一种负荷与发电机同步增长的方式，利用连续潮流法可估算出静态电压崩溃的临界点，从而确定稳定裕度，指导***调度与运行。然而该方法存在临界点附近不收敛的问题，难以得到精确的崩溃点，并且随着***规模的扩大，计算量也相当可观。因此，找到一种能快速有效地判断电压稳定临界点的方法十分重要。
电力***的电压崩溃具有隐蔽性和突发性，但通过研究发现整个***的稳定性往往与某一节点或某一区域的稳定有着密切的联系。采用面向负荷节点的戴维南等值方法，可以有效地对***网络进行简化，快速分析某一负荷节点的电压稳定性，其机理简单，物理意义明确，因而受到越来越多的关注。此外，戴维南等值方法还可用于寻找复杂***中的电压薄弱节点，从而为***的运行和调控提供重要依据。
20世纪90年代初，基于全球定位***(Global Position System，GPS)的相量测量单元(Phasor Measurement Unit，PMU)的成功研制，标志着同步相量技术的诞生。在高速通信网络的支持下，各PMU采集的带时标的数据能以较小的延时传至数据中心站，完成同步处理和分析，从而构成广域测量***(Wide Area Measurement System，WAMS)。因此，利用广域测量***的同步数据对***电压稳定进行在线分析和监测具有重要的现实意义和研究前景。
发明内容
本发明的目的是提供一种基于戴维南等值方法的电压稳定临界点判断方法，其是一种时域仿真的戴维南等值参数跟踪算法，具有良好的精度，并且适用于大扰动下暂态过程中戴维南等值参数的计算。
为了达到上述目的，本发明的基于戴维南等值方法的电压稳定临界点判断方法，其步骤为：
S1：将***等值为一个两节点戴维南等值模型；
S2：根据所述两节点戴维南等值模型，基于相量测量单元的数据采集，得出负荷节点的量测数据和戴维南等值参数，并建立负荷节点的量测数据的量测矩阵方程；
S3：对戴维南等值参数进行求解：根据所述量测矩阵方程，得出戴维南等值模型的状态方程，利用递归最小二乘法，求解戴维南等值参数的状态量的递归表达式，并更新所述状态方程的状态量；
S4：根据戴维南等值模型的状态方程得出相应的残差方程，并根据状态方程的状态量，利用线性估计状态量的不确定性，对负荷节点的量测数据的虚假数据进行检测处理；
S5：根据等阻抗模判据及负荷等值阻抗模和戴维南等值阻抗模，建立评价指标，对电压稳定临界点进行判断。
优选地，在所述步骤S1中，所述两节点戴维南等值模型方法为：将***中节点之外的电网其余部分用戴维南等效，用一个等效电压源和一个等效阻抗表示；对于任何一个时刻k，从某一负荷节点向***内看进去，***等值为一个两节点戴维南等值模型。
优选地，在所述步骤S2中，所述两节点戴维南等值模型中，负荷节点的量测数据包含k时刻负荷节点的测量电压相量和k时刻负荷节点的电流相量戴维南等值参数包含k时刻戴维南等值电势和戴维南等值阻抗Zk；
由电路理论可得：
将公式(1)表示为复数形式：
式中，VRk为测量电压相量的实部，jVIk为测量电压相量的虚部；ERk为等值电势的实部，jEIk为等值电势的实虚，RTk为等值阻抗Zk的实部，jXTk为等值阻抗Zk的虚部；IRk为电流相量的实部，jIIk为电流相量的虚部；
拆开负荷节点的量测数据和戴维南等值参数的实部和虚部，将写成矩阵形式，得负荷节点的量测数据的量测矩阵方程；
对于任意数据窗内的第k组量测数据和假定所述戴维南等值参数值不变；通过选取合适的数据窗长度，联立量测方程。
优选地，在所述步骤S3中，对戴维南等值参数进行求解的具体过程为：将所述量测矩阵方程写成戴维南等值模型的状态方程的形式：
yk＝Hkxk+εk (4)
其中，yk＝[VRk VIk]T，xk＝[ERk EIk RTk XTk]T，yk表示负荷节点量测数据中电压的实部和虚部组成的矩阵，xk表示k时刻解戴维南等值参数的的状态量；Hk是与xk对应的迭代矩阵；εk是由于负荷节点量测数据的不确定性及不良数据造成的误差；
利用递归最小二乘法，k时刻的状态量xk的递归表达式为：
其中，Gk,Hk和Pk均为迭代矩阵，用于更新状态量xk；xk-1为k-1时刻的状态量，Pk-1为k-1时刻的对应的迭代矩阵，I表示电流；
选取合适的量测长度K，具体算法如下：
a、初始化数据：x0＝[1 0 0 0]T，P0＝aI4，a>0；其中，x0为初始状态、p0为初始状态时的迭代矩阵Pk的数据，a为选定的一常数值；
b、对k＝1,2,…,K，计算中间量并更新迭代矩阵Gk，使Gk＝Pk-1HkSk；
c、更新迭代矩阵Pk，使当前k时刻的迭代状态量xk满足
对于负荷节点阻抗满足
优选地，在步骤S4中，负荷节点量测数据的不确定性会对戴维南等值参数有影响，具体分析方法为：
戴维南等值模型的状态方程为yk＝Hkxk+εk；其中，εk表示因负荷节点量测数据的不确定性及不良数据造成的误差；
可得残差方程：
r(x)＝y-Hx (6)
其中，r(x)为由于量测数据的不确定性及不良数据造成的误差；
记分别表示量测数据VRk、VIk、IRk、IIk中各自对应的实部或虚部的不确定量，表示k时刻数据窗中状态估计量xk的不确定性；记MK＝{VRk,VIk,IRk,IIk|k＝1,2,…,K}为从开始直至截止到K时刻的全部量测数据的集合，利用线性估计表示状态量的不确定性为：
其中，m是属于集合Mk的任意一量测数据，σm表示对应于量测数据m的不确定量；{R,I}表示一种集合，其包含量测数据的实部或虚部；β是集合{R,I}的元素；Iβj表示j时刻时，β所代表的量测数据的虚部；为Iβj代表的量测数据的虚部的不确定量；Vβj表示j时刻时，β所代表的量测数据的实部；δVβj为Vβj代表的量测数据的实部的不确定量；且j介于1到K之间。
优选地，在所述步骤S5中，所述等阻抗模判据为：当负荷等值阻抗模与戴维南等值阻抗模相等时，负荷功率达到极限时对应于静态电压稳定临界点；
根据负荷等值阻抗模与戴维南等值阻抗模，建立阻抗模比和阻抗模裕度来监测***是否达到电压稳定临界点；
当评价指标为阻抗模比时：
阻抗模比α＝|Zt,Thev|/|Zt,Li|，其中，|Zt,Thev|为戴维南等值阻抗模，|Zt,Li|为负荷等值阻抗模；当阻抗模比α趋向于1时，***越接近电压稳定临界点；当阻抗模比α＝1时，***运行状态处在电压稳定临界点；当α越小，负荷节点的电压越稳定；
当评价指标为阻抗模裕度时：
阻抗模裕度|ΔZ|＝(|Zt,Li|-|Zt,Thev|)/|Zt,Li|，当阻抗模裕度|ΔZ|趋向于0时，***越接近电压稳定临界点；当阻抗模裕度|ΔZ|＝0时，***运行状态处在电压稳定临界点；阻抗模裕度|ΔZ|越大，节点电压越稳定。
优选地，在所述步骤S5中，根据阻抗模比或阻抗模裕度来判断***中的薄弱的节点，具体为：
***进行戴维南等值之后，通过计算单个负荷节点的奇异值进行电压稳定性的分析和电压薄弱节点的识别；根据两节点戴维南等值模型，计算潮流雅可比矩阵其中：
其中，A、O、C和D分别为雅可比矩阵中的有功功率与相角的一阶偏导子矩阵、有功功率与电压的一阶偏导的子矩阵、无功功率与相角的一阶偏导子矩阵和无功功率与电压的一阶偏导子矩阵；PL、QL分别为负荷的有功功率和无功功率；P1、Q1分别为等值的有功功率和无功功率；Es为等值电压的幅值；δ1为等值电压的相角；P2、Q2分别为负荷节点的有功功率和无功功率；；UL为负荷节点电压的幅值，δ2为负荷节点电压的相角；G为等效电导；B为等效电纳；
再求出所述雅可比矩阵的奇异值，对奇异值进行比较，找到最小奇异值，对应的节点是最薄弱的节点。
与现有技术相比，本发明的有益效果为：
(1)本发明主要通过戴维南等值，将复杂庞大的***简化为简单的两节点***，运用等阻抗模判据分析负荷节点的电压稳定性及临界运行情况；本发明具有良好的精度，且适用于大扰动下暂态过程中戴维南等值参数的计算。
(2)本发明根据电网节点电压电流相量特性，基于PMU测量数据利用递归最小二乘法计算出电网侧等值阻抗，能够快速有效地在线实时跟踪戴维南等值参数，还进一步提出了电压稳定临界点判断指标。
附图说明
图1本发明的两节点戴维南等值模型图；
图2本发明的存在不良数据的戴维南等值结果示意图；
图3本发明的两节点***示意图；
图4a和图4b为本发明的IEEE-39***的节点15戴维南等值结果示意图。
具体实施方式
本发明提供了一种基于戴维南等值方法的电压稳定临界点判断方法，为使本发明更加明显易懂，以下结合附图和具体实施方式做进一步说明。
本发明是一种基于戴维南等值的广域测量***的电压稳定分析与监测方法，该方法主要为：
(一)从某一负荷节点向***内看进去，将整个***等值为一个简单两节点***；
在实际电网中，可以把配置有PMU(Phasor Measurement Unit，相量测量单元)的指定负荷节点之外的电网看作一个黑箱，按照电路原理，把节点之外的电网其余部分用戴维南等效，即用一个等效电压源和一个等效阻抗表示。对于任何一个时刻k，从某一负荷节点向***内看进去，可以将整个***等值为一个简单两节点***。
(二)再通过各种电压稳定指标来评估电压稳定性。具体如下所述：
如图1所示是一种两节点戴维南等值模型，其中，为k时刻负荷节点的测量电压相量，为k时刻负荷节点的电流相量，为k时刻戴维南等值电势，Zk为k时刻戴维南等值阻抗。
由电路理论可得：
将公式(1)表示为复数形式：
式中，VRk为测量电压相量的实部，jVIk为测量电压相量的虚部；ERk为等值电势的实部，jEIk为等值电势的实虚，RTk为等值阻抗Zk的实部，jXTk为等值阻抗Zk的虚部；IRk为电流相量的实部，jIIk为电流相量的虚部。
拆开上述量测数据的实部与虚部，将写成矩阵形式，可得以下量测矩阵方程：
对于任意数据窗内的第k组量测数据和假定戴维南等值参数值不变。通过选取合适的数据窗长度，将量测方程联立，可通过数学方法估算出戴维南等值参数。
以下通过最小二乘法进行戴维南等值参数求解，具体过程如下：
将量测矩阵方程写成以下状态方程的形式：
yk＝Hkxk+εk (4)
其中，yk＝[VRk VIk]T，xk＝[ERk EIk RTk XTk]T，yk表示负荷节点量测数据中的电压的实部和虚部组成的矩阵，xk表示k时刻的戴维南等值参数的状态量；Hk是与xk对应的迭代矩阵；εk是由于PMU量测数据的不确定性以及不良数据造成的误差。
对于公式(4)，可利用递归最小二乘法(The recursive least square RLS)，则k时刻的状态量xk的递归表达式为：
其中，Gk,Hk和Pk均为迭代矩阵，用于更新状态量xk；xk-1为k-1时刻的状态量，Pk-1为k-1时刻的对应的迭代矩阵，I表示电流。
选取合适的量测长度K，具体算法如下：
(a)初始化数据：x0＝[1 0 0 0]T，P0＝aI4，a>0；其中，x0为初始状态、p0为初始状态时的迭代矩阵Pk的数据，a为选定的一常数值。
(b)对k＝1,2,…,K，计算中间量并更新迭代矩阵Gk，使Gk＝Pk-1HkSk；
(c)更新迭代矩阵Pk，使当前k时刻的迭代状态量xk满足
对于负荷节点阻抗，满足再根据等阻抗模判据，进而进一步对电压稳定与崩溃点分析。
其中，PMU数据不确定性会对等值参数的有影响，具体分析如下：
戴维南等值模型的状态方程为yk＝Hkxk+εk，其中εk表示因PMU量测数据的不确定性以及不良数据造成的误差；可得残差方程：
r(x)＝y-Hx (6)
其中，r(x)为由于PMU量测数据的不确定性以及不良数据造成的误差。
在现有电压稳定分析研究中，多将PMU数据处理为方差较小的高斯白噪声，但当出现重大误差时(如脉冲噪声、数据异步性、网络攻击等)，会导致戴维南等值参数的重大误差从而造成电压稳定与崩溃点判别的错误。
记分别表示PMU量测数据VRk、VIk、IRk、IIk中各自对应的实部或虚部的不确定量，表示当前的k时刻数据窗中状态估计量xk的不确定性；记MK＝{VRk,VIk,IRk,IIk|k＝1,2,…,K}为从开始直至截止到K时刻的全部量测数据的集合，则利用线性估计表示状态量的不确定性为：
其中，m是属于集合Mk的任意一量测数据，σm表示对应于量测数据m的不确定量；{R,I}表示一种集合，其包含PMU量测数据的实部或虚部；β是集合{R,I}的元素；Iβj表示j时刻时，β所代表的量测数据的虚部；为Iβj代表的量测数据的虚部的不确定量；Vβj表示j时刻时，β所代表的量测数据的实部；δVβj为Vβj代表的量测数据的实部的不确定量；且j介于1到K之间。
当数据中存在不良数据(即虚假数据)时，若攻击者在已知***结构与状态方程情况下发起数据攻击，会导致等值结果出现极端的情况。
如图2所示，横坐标为负荷增长率(％)，纵坐标为阻抗值。其中，曲线L1是负荷增长与通过鲁棒递归最小二乘法((Robust RLS))求解的戴维南等值阻抗的关系曲线，曲线L2是负荷增长与通过递归最小二乘法(RLS)求解的戴维南等值阻抗的关系曲线，曲线L2是与曲线L1重合的一部分。曲线L3是负荷增长和负荷节点的等值阻抗的关系曲线。
当负荷增长率为183％时：曲线L2的纵坐标等于曲线L3的纵坐标，而曲线L1的纵坐标小于曲线L3的纵坐标。此时会出现当由递归最小二乘法求解的戴维南等值判别出***电压失稳，而实际***稳定，或者；由递归最小二乘法求解的戴维南等值判别出***稳定，而实际已达电压崩溃点。
所以，对于电力***的量测数据中不良虚假数据的检测与处理仍是电力***状态估计与分析的重要方向。
本发明中，基于戴维南等值的电压稳定临界点判断方法如下：
(1)关于电压稳定临界点的判断
根据等阻抗模判据可知，当负荷等值阻抗模与***戴维南等值阻抗模相等时，负荷功率达到极限，此时对应于静态电压稳定临界点。因此可以直接通过比较负荷等值阻抗模与***戴维南等值阻抗模进行电压稳定临界点的判断。
根据负荷等值阻抗模与***戴维南等值阻抗模，还可以建立多种指标，如阻抗模比、阻抗模裕度等，用来监测***是否达到电压稳定临界点。
(a)阻抗模比
定义阻抗模比α＝|Zt,Thev|/|Zt,Li|，其中，|Zt,Thev|为戴维南等值阻抗模，|Zt,Li|为负荷等值阻抗模。当阻抗模比α趋向于1时，则***越来越接近电压稳定临界点；当阻抗模比α＝1时，***运行状态正好处在电压稳定临界点。所以当α越小，负荷节点的电压越稳定。
(b)阻抗模裕度
定义阻抗模裕度|ΔZ|＝(|Zt,Li|-|Zt,Thev|)/|Zt,Li|，当阻抗模裕度|ΔZ|趋向于0时，***越来越接近电压稳定临界点；当阻抗模裕度|ΔZ|＝0时，***运行状态正好处在电压稳定临界点。与阻抗模α相比，阻抗模裕度|ΔZ|能够体现节点正常工作状态“距离”其临界点有多“远”，且|ΔZ|阻抗模裕度越大，则表示节点电压越稳定。
另，还可以建立其他类似的指标，但本质都是采用等阻抗模判据来判断***的电压稳定临界点。
(2)关于电压薄弱点的识别
理论上，根据阻抗模比或阻抗模裕度两种指标，即可判断***中哪些节点的电压比较薄弱(即阻抗模比大或者阻抗模裕度小为负荷节点)。
从另一个角度看，由于把***进行戴维南等值之后，则网络结构得到了极大地简化，计算维数大大降低，从而可通过计算单个负荷节点的奇异值来进行电压稳定性的分析和电压薄弱节点的识别。
如图3所示为两节点***示意图，可计算潮流雅可比矩阵，其为 其中：
其中，A、O、C和D分别为雅可比矩阵中的有功功率与相角的一阶偏导子矩阵、有功功率与电压的一阶偏导的子矩阵、无功功率与相角的一阶偏导子矩阵和无功功率与电压的一阶偏导子矩阵；PL、QL分别为负荷的有功功率和无功功率；P1、Q1分别为等值的有功功率和无功功率；Es为等值电压的幅值；δ1为等值电压的相角；P2、Q2分别为负荷节点的有功功率和无功功率；；UL为负荷节点电压的幅值，δ2为负荷节点电压的相角；G为等效电导；B为等效电纳。
所以，在任意运行状态，对***的所有负荷节点进行戴维南等值，然后计算等效的雅可比矩阵，再求出雅可比矩阵的奇异值，最后经过比较找到最小奇异值，对应的节点就是最薄弱的节点。
示例地，为了验证所提算法的等值精度及有效性，可对IEEE-39节点***进行了计算分析。如图4a和图4b是IEEE-39中的15号节点的戴维南等值结果示意图。在IEEE-39节点***中，设定其增长方式为：***中所有负荷的功率(功率因数保持与初始状态相同)和发电机的有功出力均以同一比例增长，直到潮流不再收敛。
如图4a所示，横坐标表示负荷增长指数λ，纵坐标表示阻抗值。其中，曲线S1表示戴维南等值阻抗的变化趋势，曲线S2表示负荷阻抗的变化趋势。曲线S1中，当负荷增长指数λ逐渐增大时，戴维南等值阻抗的值逐渐减小；曲线S2中，当负荷增长指数λ逐渐增大时，负荷阻抗的值逐渐增大；且戴维南等值阻抗小于负荷阻抗的值。
如图4b所示为阻抗模比与负荷增长指数λ的变化曲线示意图，横坐标为负荷增长指数λ，纵坐标为阻抗模比的值。当负荷增长指数λ逐渐增大时，阻抗模比的值逐渐由0增大接近1，且阻抗模比的值始终介于0和1之间。即当***的潮流无法收敛(重负荷电压崩溃)时，戴维南等值阻抗模与负荷阻抗模在误差范围内近似相等，即阻抗模比约为1。
若当对该***中所有的负荷节点进行戴维南等值，并计算电压临界指标，指标变化明显的对应的节点就是在该***接近电压稳定临界点时最薄弱的节点。***还要考虑到数据传输过程中的噪声误差，与如脉冲噪声、数据异步性、网络攻击等对等值结果造成的较大误差的影响。
本发明中的戴维南等值参数的电压稳定临界点判断方法在大***中仍然是适用的，在判断***电压稳定临界点的同时，还能识别出电压比较薄弱的节点，可以为电力***的运行规划、离线计算、事故分析等提供参考依据。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
Claims (7)
1.一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，该方法的步骤为：
S1：将***等值为一个两节点戴维南等值模型；
S2：根据所述两节点戴维南等值模型，基于相量测量单元的数据采集，得出负荷节点的量测数据和戴维南等值参数，并建立负荷节点的量测数据的量测矩阵方程；
S3：对戴维南等值参数进行求解：根据所述量测矩阵方程，得出戴维南等值模型的状态方程，利用递归最小二乘法，求解戴维南等值参数的状态量的递归表达式，并更新所述状态方程的状态量；
S4：根据戴维南等值模型的状态方程得出相应的残差方程，并根据状态方程的状态量，利用线性估计状态量的不确定性，对负荷节点的量测数据的虚假数据进行检测处理；
S5：根据等阻抗模判据及负荷等值阻抗模和戴维南等值阻抗模，建立评价指标，对电压稳定临界点进行判断。
2.如权利要求1所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在所述步骤S1中，所述两节点戴维南等值模型方法为：将***中节点之外的电网其余部分用戴维南等效，用一个等效电压源和一个等效阻抗表示；对于任何一个时刻k，从某一负荷节点向***内看进去，***等值为一个两节点戴维南等值模型。
3.如权利要求2所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在所述步骤S2中，所述两节点戴维南等值模型中，负荷节点的量测数据包含k时刻负荷节点的测量电压相量和k时刻负荷节点的电流相量戴维南等值参数包含k时刻戴维南等值电势和戴维南等值阻抗Zk；
由电路理论可得：
<mrow>
<mover>
<msub>
<mi>V</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mover>
<msub>
<mi>E</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>k</mi>
</msub>
<mover>
<msub>
<mi>I</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
将公式(1)表示为复数形式：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>V</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>jV</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>E</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>jE</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>T</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mrow>
<mi>T</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>I</mi>
<mi>k</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>jI</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，VRk为测量电压相量的实部，jVIk为测量电压相量的虚部；ERk为等值电势的实部，jEIk为等值电势的实虚，RTk为等值阻抗Zk的实部，jXTk为等值阻抗Zk的虚部；IRk为电流相量的实部，jIIk为电流相量的虚部；
拆开负荷节点的量测数据和戴维南等值参数的实部和虚部，将写成矩阵形式，得负荷节点的量测数据的量测矩阵方程；
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>I</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>I</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mrow>
<mi>R</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>I</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
对于任意数据窗内的第k组量测数据和假定所述戴维南等值参数值不变；通过选取合适的数据窗长度，联立量测方程。
4.如权利要求3所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在所述步骤S3中，对戴维南等值参数进行求解的具体过程为：
将所述量测矩阵方程写成戴维南等值模型的状态方程的形式：
yk＝Hkxk+εk (4)
其中，yk＝[VRk VIk]T，xk＝[ERk EIk RTk XTk]T，yk表示负荷节点量测数据中电压的实部和虚部组成的矩阵，xk表示k时刻解戴维南等值参数的的状态量；Hk是与xk对应的迭代矩阵；εk是由于负荷节点量测数据的不确定性及不良数据造成的误差；
利用递归最小二乘法，k时刻的状态量xk的递归表达式为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>+</mo>
<msubsup>
<mi>H</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>H</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>k</mi>
</msub>
<msubsup>
<mi>H</mi>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Gk,Hk和Pk均为迭代矩阵，用于更新状态量xk；xk-1为k-1时刻的状态量，Pk-1为k-1时刻的对应的迭代矩阵，I表示电流；
选取合适的量测长度K，具体算法如下：
a、初始化数据：x0＝[1 0 0 0]T，P0＝aI4，a>0；其中，x0为初始状态、p0为初始状态时的迭代矩阵Pk的数据，a为选定的一常数值；
b、对k＝1,2,…,K，计算中间量并更新迭代矩阵Gk，使Gk＝Pk- 1HkSk；
c、更新迭代矩阵Pk，使当前k时刻的迭代状态量xk满足
对于负荷节点阻抗满足
5.如权利要求4所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在步骤S4中，负荷节点量测数据的不确定性会对戴维南等值参数有影响，具体分析方法为：
戴维南等值模型的状态方程为yk＝Hkxk+εk；其中，εk表示因负荷节点量测数据的不确定性及不良数据造成的误差；
可得残差方程：
r(x)＝y-Hx (6)
其中，r(x)为由于量测数据的不确定性及不良数据造成的误差；
记分别表示量测数据VRk、VIk、IRk、IIk中各自对应的实部或虚部的不确定量，表示k时刻数据窗中状态估计量xk的不确定性；记MK＝{VRk,VIk,IRk,IIk|k＝1,2,…,K}为从开始直至截止到K时刻的全部量测数据的集合，利用线性估计表示状态量的不确定性为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>x</mi>
<mi>K</mi>
</msub>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mi>K</mi>
</msub>
</mrow>
</munder>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>m</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mi>m</mi>
</msub>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msqrt>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&beta;</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mi>R</mi>
<mo>,</mo>
<mi>I</mi>
<mo>}</mo>
</mrow>
</munder>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>&beta;</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>I</mi>
<mrow>
<mi>&beta;</mi>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&beta;</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mi>R</mi>
<mo>,</mo>
<mi>I</mi>
<mo>}</mo>
</mrow>
</munder>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>&beta;</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>V</mi>
<mrow>
<mi>&beta;</mi>
<mi>j</mi>
</mrow>
</msub>
</msub>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，m是属于集合Mk的任意一量测数据，σm表示对应于量测数据m的不确定量；{R,I}表示一种集合，其包含量测数据的实部或虚部；β是集合{R,I}的元素；Iβj表示j时刻时，β所代表的量测数据的虚部；为Iβj代表的量测数据的虚部的不确定量；Vβj表示j时刻时，β所代表的量测数据的实部；为Vβj代表的量测数据的实部的不确定量；且j介于1到K之间。
6.如权利要求5所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在所述步骤S5中，所述等阻抗模判据为：当负荷等值阻抗模与戴维南等值阻抗模相等时，负荷功率达到极限时对应于静态电压稳定临界点；
根据负荷等值阻抗模与戴维南等值阻抗模，建立阻抗模比和阻抗模裕度来监测***是否达到电压稳定临界点；
当评价指标为阻抗模比时：
阻抗模比α＝|Zt,Thev|/|Zt,Li|，其中，|Zt,Thev|为戴维南等值阻抗模，/|Zt,Li|为负荷等值阻抗模；当阻抗模比α趋向于1时，***越接近电压稳定临界点；当阻抗模比α＝1时，***运行状态处在电压稳定临界点；当α越小，负荷节点的电压越稳定；
当评价指标为阻抗模裕度时：
阻抗模裕度|ΔZ|＝(|Zt,Li|-|Zt,Thev|)/|Zt,Li|，当阻抗模裕度|ΔZ|趋向于0时，***越接近电压稳定临界点；当阻抗模裕度|ΔZ|＝0时，***运行状态处在电压稳定临界点；阻抗模裕度|ΔZ|越大，节点电压越稳定。
7.如权利要求6所述的一种基于戴维南等值方法的电压稳定临界点判断方法，其特征在于，
在所述步骤S5中，根据阻抗模比或阻抗模裕度来判断***中的薄弱的节点，具体为：
***进行戴维南等值之后，通过计算单个负荷节点的奇异值进行电压稳定性的分析和电压薄弱节点的识别；根据两节点戴维南等值模型，计算潮流雅可比矩阵其中：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mo>=</mo>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>/</mo>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>E</mi>
<mi>s</mi>
</msub>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>G</mi>
<mi> </mi>
<msub>
<mi>sin&delta;</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<mi>B</mi>
<mi> </mi>
<msub>
<mi>cos&delta;</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>O</mi>
<mo>=</mo>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>/</mo>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mi>G</mi>
<mo>+</mo>
<msub>
<mi>E</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>G</mi>
<mi> </mi>
<msub>
<mi>cos&delta;</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<mi>B</mi>
<mi> </mi>
<msub>
<mi>sin&delta;</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>C</mi>
<mo>=</mo>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>/</mo>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>E</mi>
<mi>s</mi>
</msub>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>G</mi>
<mi> </mi>
<msub>
<mi>cos&delta;</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<mi>B</mi>
<mi> </mi>
<msub>
<mi>sin&delta;</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mo>=</mo>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>/</mo>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>U</mi>
<mi>L</mi>
</msub>
<mi>G</mi>
<mo>+</mo>
<msub>
<mi>E</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>G</mi>
<mi> </mi>
<msub>
<mi>sin&delta;</mi>
<mn>21</mn>
</msub>
<mo>+</mo>
<mi>B</mi>
<mi> </mi>
<msub>
<mi>cos&delta;</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，A、O、C和D分别为雅可比矩阵中的有功功率与相角的一阶偏导子矩阵、有功功率与电压的一阶偏导的子矩阵、无功功率与相角的一阶偏导子矩阵和无功功率与电压的一阶偏导子矩阵；PL、QL分别为负荷的有功功率和无功功率；P1、Q1分别为等值的有功功率和无功功率；Es为等值电压的幅值；δ1为等值电压的相角；P2、Q2分别为负荷节点的有功功率和无功功率；；UL为负荷节点电压的幅值，δ2为负荷节点电压的相角；G为等效电导；B为等效电纳；
再求出所述雅可比矩阵的奇异值，对奇异值进行比较，找到最小奇异值，对应的节点是最薄弱的节点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710889796.4A CN107658881A (zh) | 2017-09-27 | 2017-09-27 | 基于戴维南等值方法的电压稳定临界点判断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710889796.4A CN107658881A (zh) | 2017-09-27 | 2017-09-27 | 基于戴维南等值方法的电压稳定临界点判断方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107658881A true CN107658881A (zh) | 2018-02-02 |
Family
ID=61117166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710889796.4A Pending CN107658881A (zh) | 2017-09-27 | 2017-09-27 | 基于戴维南等值方法的电压稳定临界点判断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107658881A (zh) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108599184A (zh) * | 2018-05-16 | 2018-09-28 | 福州大学 | 一种电力***的有效无功补偿地点的离线选择方法 |
| CN108899897A (zh) * | 2018-07-12 | 2018-11-27 | 哈尔滨工业大学 | 数据驱动的电力***电压薄弱节点在线辨识方法 |
| CN108923439A (zh) * | 2018-07-10 | 2018-11-30 | 福州大学 | 一种电力***有效无功补偿地点的在线选择方法 |
| CN109524960A (zh) * | 2018-11-23 | 2019-03-26 | 国网江苏省电力有限公司 | 自治-协调的电力***两级分布式静态电压稳定判断方法 |
| CN110034558A (zh) * | 2019-04-17 | 2019-07-19 | 国家电网有限公司 | 基于同步相量量测的配电网戴维南等值参数鲁棒估计方法 |
| CN110138241A (zh) * | 2019-06-06 | 2019-08-16 | 西南交通大学 | 一种基于雅可比理论的单相电压源的稳定性分析方法 |
| CN110571790A (zh) * | 2019-03-21 | 2019-12-13 | 中国电力科学研究院有限公司 | 基于戴维南等值的在线电压稳定预警的方法、*** |
| CN110781638A (zh) * | 2019-11-07 | 2020-02-11 | 山东大学 | 基于节点端口伏安特性的最优戴维南等值参数计算方法 |
| CN111797436A (zh) * | 2020-09-10 | 2020-10-20 | 深圳华工能源技术有限公司 | 一种配用电***节能设备节能数据伪造识别方法 |
| CN112260273A (zh) * | 2020-10-13 | 2021-01-22 | 武汉大学 | 一种基于全纯嵌入法的电力***薄弱节点可视化识别方法 |
| CN112421619A (zh) * | 2020-11-09 | 2021-02-26 | 广东电网有限责任公司江门供电局 | 一种智能低压台区配电网电压稳定性实时监测方法 |
| CN112859638A (zh) * | 2019-11-28 | 2021-05-28 | Abb瑞士股份有限公司 | 确定转换器***的戴维南等效模型 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1134867A1 (en) * | 2000-03-10 | 2001-09-19 | ABB Power Automation AG | Method and device for assessing the stability of an electric power transmission network |
| CN103310390A (zh) * | 2013-07-05 | 2013-09-18 | 武汉大学 | 一种电网安全综合评价方法 |
| CN106682407A (zh) * | 2016-12-19 | 2017-05-17 | 国网上海市电力公司 | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 |
-
2017
- 2017-09-27 CN CN201710889796.4A patent/CN107658881A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1134867A1 (en) * | 2000-03-10 | 2001-09-19 | ABB Power Automation AG | Method and device for assessing the stability of an electric power transmission network |
| CN103310390A (zh) * | 2013-07-05 | 2013-09-18 | 武汉大学 | 一种电网安全综合评价方法 |
| CN106682407A (zh) * | 2016-12-19 | 2017-05-17 | 国网上海市电力公司 | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 |
Non-Patent Citations (2)
| Title |
|---|
| CHENCHEN等: ""PMU Uncertainty Quantification in Voltage Stability Analysis"", 《IEEE》 * |
| 李日波: "一种电网戴维南等值及静态电压稳定性研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108599184B (zh) * | 2018-05-16 | 2021-04-27 | 福州大学 | 一种电力***的有效无功补偿地点的离线选择方法 |
| CN108599184A (zh) * | 2018-05-16 | 2018-09-28 | 福州大学 | 一种电力***的有效无功补偿地点的离线选择方法 |
| CN108923439A (zh) * | 2018-07-10 | 2018-11-30 | 福州大学 | 一种电力***有效无功补偿地点的在线选择方法 |
| CN108899897A (zh) * | 2018-07-12 | 2018-11-27 | 哈尔滨工业大学 | 数据驱动的电力***电压薄弱节点在线辨识方法 |
| CN109524960A (zh) * | 2018-11-23 | 2019-03-26 | 国网江苏省电力有限公司 | 自治-协调的电力***两级分布式静态电压稳定判断方法 |
| CN109524960B (zh) * | 2018-11-23 | 2021-10-15 | 国网江苏省电力有限公司 | 自治-协调的电力***两级分布式静态电压稳定判断方法 |
| CN110571790B (zh) * | 2019-03-21 | 2021-11-16 | 中国电力科学研究院有限公司 | 基于戴维南等值的在线电压稳定预警的方法、*** |
| CN110571790A (zh) * | 2019-03-21 | 2019-12-13 | 中国电力科学研究院有限公司 | 基于戴维南等值的在线电压稳定预警的方法、*** |
| CN110034558A (zh) * | 2019-04-17 | 2019-07-19 | 国家电网有限公司 | 基于同步相量量测的配电网戴维南等值参数鲁棒估计方法 |
| CN110138241B (zh) * | 2019-06-06 | 2020-09-25 | 西南交通大学 | 一种基于雅可比理论的单相电压源的稳定性分析方法 |
| CN110138241A (zh) * | 2019-06-06 | 2019-08-16 | 西南交通大学 | 一种基于雅可比理论的单相电压源的稳定性分析方法 |
| CN110781638A (zh) * | 2019-11-07 | 2020-02-11 | 山东大学 | 基于节点端口伏安特性的最优戴维南等值参数计算方法 |
| CN112859638A (zh) * | 2019-11-28 | 2021-05-28 | Abb瑞士股份有限公司 | 确定转换器***的戴维南等效模型 |
| CN111797436B (zh) * | 2020-09-10 | 2020-12-25 | 深圳华工能源技术有限公司 | 一种配用电***节能设备节能数据伪造识别方法 |
| CN111797436A (zh) * | 2020-09-10 | 2020-10-20 | 深圳华工能源技术有限公司 | 一种配用电***节能设备节能数据伪造识别方法 |
| CN112260273A (zh) * | 2020-10-13 | 2021-01-22 | 武汉大学 | 一种基于全纯嵌入法的电力***薄弱节点可视化识别方法 |
| CN112421619A (zh) * | 2020-11-09 | 2021-02-26 | 广东电网有限责任公司江门供电局 | 一种智能低压台区配电网电压稳定性实时监测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107658881A (zh) | 基于戴维南等值方法的电压稳定临界点判断方法 | |
| Liu et al. | Comparisons on Kalman-filter-based dynamic state estimation algorithms of power systems | |
| CN107577870B (zh) | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 | |
| CN110676940B (zh) | 一种提高参数辨识精度的配电网pmu配置方法及*** | |
| CN108155648A (zh) | 基于自适应h无穷扩展卡尔曼滤波的状态估计方法 | |
| CN109459705A (zh) | 一种抗野值鲁棒无迹卡尔曼滤波的动力电池soc估计方法 | |
| CN107453484B (zh) | 一种基于wams信息的scada数据校准方法 | |
| CN103886193B (zh) | 一种电力***模糊自适应抗差估计方法 | |
| CN109754013B (zh) | 一种基于无迹卡尔曼滤波的电力***混合量测融合方法 | |
| CN110190600A (zh) | 一种基于ami量测近邻回归的三相配电网拓扑辨识方法 | |
| CN115000947A (zh) | 基于智能电表量测的配电网拓扑结构与线路参数辨识方法 | |
| Sheng et al. | Online measurement-based estimation of dynamic system state matrix in ambient conditions | |
| CN106372440B (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
| CN107025609A (zh) | 基于奇异值分解cdkf的电力***动态状态估计方法 | |
| CN109146336B (zh) | 一种基于t分布的电力***鲁棒状态估计方法 | |
| CN107958120B (zh) | 一种基于幂级数展开的***戴维南等值参数计算方法 | |
| CN109638811B (zh) | 基于模型等值的配电网电压功率灵敏度鲁棒估计方法 | |
| CN105932670A (zh) | 一种基于半定规划的配电网状态估计方法 | |
| CN112865096A (zh) | 一种考虑pmu量测相角偏差的配电网状态估计方法及*** | |
| CN109193639B (zh) | 一种电力***抗差估计方法 | |
| CN115293090A (zh) | 基于重构数据处理的多谐波源责任量化方法 | |
| CN113489072B (zh) | 一种考虑支路功率伪量测的配电网状态估计方法 | |
| CN115795222A (zh) | 基于同步相量测量的电网不良参数辨识和校正方法 | |
| CN106355029B (zh) | 一种配电***关键量测解耦快速辨识的方法 | |
| CN113553538B (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180202 |