CN101383511B - 基于数据采集系统量测数据的电力系统状态估计方法 - Google Patents
基于数据采集系统量测数据的电力系统状态估计方法 Download PDFInfo
- Publication number
- CN101383511B CN101383511B CN2008102237660A CN200810223766A CN101383511B CN 101383511 B CN101383511 B CN 101383511B CN 2008102237660 A CN2008102237660 A CN 2008102237660A CN 200810223766 A CN200810223766 A CN 200810223766A CN 101383511 B CN101383511 B CN 101383511B
- Authority
- CN
- China
- Prior art keywords
- node
- power
- theta
- value
- generator
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000005259 measurement Methods 0.000 title claims description 52
- 238000013480 data collection Methods 0.000 title 1
- 238000012544 monitoring process Methods 0.000 claims abstract description 3
- 230000011664 signaling Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000010248 power generation Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000013499 data model Methods 0.000 description 2
- SDNXQWUJWNTDCC-UHFFFAOYSA-N 2-methylsulfonylethanamine Chemical compound CS(=O)(=O)CCN SDNXQWUJWNTDCC-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
Abstract
基于数据采集系统量测数据的电力系统状态估计方法属于电力系统分析技术领域,其特征在于,基于电力系统网络参数,以及数据采集和监控系统实时测量的量测数据,把所述量测数据分配到正在运行的一个电力系统中,形成一个计算模型,把联络节点注入功率为零以及发电机最大有功出力和最大无功出力作为约束,建立估计用的数学模型,求解各节点电压幅值和电压相角的最优解,并用状态变量最优解计算电压,有功和无功的估计值,若量测值和估计值之差小于设定的不同电压等级下的参数α,则估计值合格。本发明提出的估计准则有很强的抗差性,不易受不良数据的影响,提高了估计的正确性。
Description
技术领域
本技术发明属于电力系统分析计算领域。
背景技术
状态估计也被称为滤波,它是利用实时量测系统的冗余度来提高数据精度,自动排除随机干扰所引起的错误信息,估计或预报系统的运行状态。状态估计的功能是将低精度、不完整、有不良数据的生数据转为完整、合理、误差小的熟数据。其主要任务是利用实时量测系统的遥测、遥信信息,按开关状态建立网络模型,并对实时量测数据进行粗检测,剔除其中的不良数据,对不足的量测点补充伪量测,以保证冗余度,从而估计出系统的状态,计算出系统的潮流分布,同时利用计算结果分析系统的运行状态,加强全网的可观测性。
传统的状态估计方法于1970年由MIT的Fred Charles Schweppe创立,迄今已有30多年的历史。根据求解状态变量所用的目标函数的不同,状态估计器有不同的估计准则,目前电力系统状态估计的估计准则有:加权最小二乘准则(WLS)、非二次(non-quadratic)准则、加权最小绝对值(WLAV)、LMS(least median of squares)和LTS(least trimmed squares)。当前应用最广泛的是WLS估计准则,它的优点是模型简单,计算量小,对理想正态分布的量测量,估计具有最优一致且无偏等优良统计特性,缺点是抗差能力差,即估计结果易受不良数据的影响而偏离真值较远。非二次准则、WLAV、LMS和LTS估计器都属于抗差估计器,它们的显著优点就是具有抗差能力,但有计算量大、需要主观确定加权因子等缺点,这也是妨碍抗差估计器发展的重要原因。
本发明提出了一种电力系统状态估计方法,该方法提出了一种估计准则,该估计准则具有很强的抗差性,可反映量测的真实分布情况,建立了在此估计准则基础上的优化问题数学模型,将联络节点的注入功率作为等式约束放在约束条件中,提高了估计的正确性,给出了目标函数中参数的设定原则,以及求解此优化问题时初值的设定方法。
发明内容
本发明的目的在于提供一种基于实时量测的电力系统状态估计方法。
1.基于数据采集系统量测数据的电力系统状态估计方法其特征在于所述方法是在计算机中依次按以下步骤实现的:
步骤(1),从能量管理系统中获得电力系统的网络参数,包括:输电线路的支路号、首节点和末节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器的变比和阻抗;
步骤(2),从数据采集和监视系统SCADA的网络数据服务器获得遥测量测值和遥信量测值,其中
遥测量测值包括:电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率、线路末端无功功率、线路电流幅值,
遥信量测值包括:开关、刀闸开合状态以及变压器分接头位置,开关或刀闸闭合的量测值为1,断开的量测值为0,
所述SCADA系统采用轮询方式采集数据,轮询周期为5~10s,每一个测点至少包括:测点类型、测点号以及量测值;
步骤(3),构造计算模型
步骤(3.1),把开关、输电线、变压器、发电机、电容器以及负荷的端口定义为节点,遍历所述这些节点,根据得到的开关量把电网中通过闭合的开关或刀闸连接在一起的节点合并成一个计算用节点,下面所述节点均指计算用节点,
步骤(3.2),把步骤(2)得到的电压幅值量测值、发电机功率量测值以及负荷功率量测值按以下方法分配到对应的节点上:
若所述节点上有多个电压幅值量测值,则取平均值作为一个电压幅值量测值,
若只有所述节点上所有的发电机和负荷都有功率量测时,才把发电机功率减去负荷功率作为该节点上的注入功率;
步骤(4),以节点电压幅值与电压相角作为状态变量x,并取某一节点作为电压相角参考节点,便按下式求所述状态变量最优解x*:
步骤(4.1),把没有发电机也没有负荷的节点定义为联络节点,该联络节点的集合为A,而发电机节点集合为E,用Pkmax、Qkmax分别表示节点号为k的发电机节点最大有功出力和最大无功出力,
步骤(4.2),采用下述优化模型求解状态变量最优解x*:
其中,A为联络节点集合,E为发电机节点集合,
vt,vj,vk,vl分别为节点t、节点j、节点k、节点l的电压幅值,
j∈t表示节点j与节点t有线路相连,
l∈k表示节点l与节点k有线路相连,
节点t为联络节点,节点k为发电机节点,
Gtj,Btj,Gkl,Bkl为导纳矩阵元素,
θtj为节点t和节点j的电压相角差,
θkl为节点k和节点l的电压相角差,
zi为第i个量测的量测值,
hi(x)为第i个量测的量测函数,
αi是和第i个量测有关的参数:
对于电压测点:0.02×|zibase|,
对于有功测点:0.02×|zibase|,
对于无功测点:0.03×|zibase|,
zibase的大小和第i个量测有关,对于不同的电压等级,zibase是设定的,
βi是一个与所述αi相关的一个参数,
βi=c/αi,c是一个常数,c=5~10;
步骤(5),按照步骤(4)得到的最优解x*,计算出量测值的估计值z*=h(x*),节点i电压幅值,
节点i有功注入功率,
节点i无功注入功率,
2.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法其特征在于:
对于线路功率量测,所述zibase按以下方式取值:500kv等级取1083MVA,330kv等级取686MVA,220kv等级取305MVA,110kv等级取114MVA,66kv电压等级取69.7MVA。
3.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法其特征在于:
对于电压幅值量测,所述zibase按以下方式取值:500kv等级取600kv,330kv等级取396kv,220kv等级取264kv,110kv等级取132kv,66kv等级取79.2kv。
4.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法其特征在于:
对于发电机功率量测,所述zibase取发电机视在功率。
本发明提出的电力系统状态估计计算方法可以应用到我国各大区域和省级电网能量管理系统(EMS)中,也可以应用于基于实时量测量的电力系统监视、分析和控制系统。
附图说明
图1是元件等值电路图,其中图(a)是线路π形等值电路图,图(b)是变压器π形等值电路图。
图2是本发明所述方法的流程框图。
图3是硬件架构示意图。
图4是四节点系统图。
图5是四节点系统计算模型图。
具体实施方式
SCADA(Supervisory Control And Data Acquisition)系统,即数据采集与监视控制系统是电力系统自动化的实时数据源,提供电力系统实时量测数据。SCADA量测可分为遥测量和遥信量。遥测量是指模拟量的量测,遥信量是指离散量的量测。电力系统的遥测量方程可表示为
z=h(x)+v (1)
(1)式中z为量测值矢量;h(x)为量测函数矢量;v为量测误差矢量;设量测量共m个,则上述矢量均为m维;x为状态量,设系统节点数为n,以节点电压幅值与电压相角作为状态变量,并取某一个节点作为电压相角参考节点,则x为2n-1维,(2)~(12)式分别为几种常见量测的量测方程:
节点i注入功率Pi,Qi:
线路i-j上节点i侧潮流:
Pij=vi 2g-vivjgcosθij-vivjbsinθij(4)
Qij=-vi 2(b+yc)-vivjgsinθij+vivjbcosθij(5)
线路i-j上节点j侧潮流:
Pji=vj 2g-vivjgcosθij+vivjbsinθij(6)
Qji=-vj 2(b+yc)+vivjgsinθij+vivjbcosθij(7)
变压器支路i-j上节点i侧潮流:
变压器支路i-j上节点j侧潮流:
节点i电压幅值量测:
Ui=vi(12)
上面量测方程中,vi,vj分别为节点i和节点j的电压幅值,θij为节点i和节点j电压相角差,Gij,Bij为导纳矩阵元素,g,b,yc为线路π形模型中的参数,K为变压器的非标准变比,bT为变压器标准测的电纳,线路和变压器的π形等值电路如图1所示。
既没有发电机也没有负荷的节点称为联络节点,设联络节点集合为A,发电机节点集合为E,用Pkmax、Qkmax分别表示节点号为k的发电机节点最大有功出力和最大无功出力,采用如下优化问题数据模型:
(13)式中αi,βi是和第i个量测有关的参数,参数的取值原则如下:
1.αi的取值
●电压测点:0.02×|zibase|
●有功测点:0.02×|zibase|
●无功测点:0.03×|zibase|
zibase的大小和第i个量测有关,取值原则如下:
a)对于线路功率量测,500kv等级取1083MVA,330kv等级取686MVA,220kv等级取305MVA,110kv等级取114MVA,66kv电压等级取69.7MVA
b)对于电压幅值量测,500kv等级取600kv,330kv等级取396kv,220kv等级取264kv,110kv等级取132kv,66kv等级取79.2kv
c)对于发电机功率量测取发电机视在功率
2.βi的取值
βi=c/αi,c是一个常数,本发明中c的取值范围为5~10。
根据目前国家电网公司对测点估计值是否合格的判定标准,当第i个量测的估计值和量测值之间的偏差满足|zi-hi(x)|≤αi时,测点估计值是合格的。
设求解得到的最优解为x*,用x*计算出量测值的估计值z*=h(x*),整个计算过程结束。
本发明应用于某实际电力系统实时控制系统——AEMS,程序流程框图如图2所示,该系统硬件架构图如图3所示。该系统的服务器包括两台数据服务器和两台应用服务器,服务器操作系统均采用Solaris8,两台数据服务器及两台应用服务器均互为热备用。该系统需要同时进入现有能量管理系统(EMS)网络和SCADA网络,并保持一定的独立性。该系统启用了两个交换机,通过它们分别连接到EMS网络和SCADA网络。数据服务器都采用双网卡设计,分别与上述两个交换机相联接。这样就可以实现AEMS分区网络与EMS网络、SCADA网络之间的相互访问了。
以图4所示的4节点系统为例,具体步骤说明如下:
步骤1:初始化
从现有的能量管理系统中获得电力系统的网络参数。包括输电线路的串联电阻、串联电抗、并联电导和并联电纳;变压器的变比和阻抗。
4节点系统线路参数如表1所示:
表1 线路参数表
支路号 | 首节点 | 末节点 | 电阻R(Ω) | 电抗R(Ω) | 对地电纳y<sub>c</sub> | 变比K |
1 | 1 | 2 | 0.52 | 2.66 | 0.00014 | — |
2 | 1 | 3 | 3.05 | 8.37 | 0.00056 | — |
3 | 2 | 3 | 0.41 | 2.15 | 0.00044 | — |
4 | 3 | 4 | 0.00 | 7.50 | — | 1.05 |
步骤2:获得SCADA量测值
量测量包括遥测量测和遥信量测。遥测量测主要包括电压幅值量测,发电机有功功率量测,发电机无功功率量测,负荷有功功率量测,负荷无功功率量测,线路首端有功功率量测,线路首端无功功率量测,线路末端有功功率量测,线路末端无功功率量测,线路电流幅值量测。遥信量测主要是开关、刀闸开合状态的量测以及变压器分接头位置量测。每一个测点至少要包括三个内容:测点类型,测点号以及量测值。离线状态下可以通过文件给出上述量测量,在线应用时可以通过SCADA系统提供的API或使用TCP/IP通信协议从SCADA系统中获得。SCADA系统采用轮询的方式采集数据,轮询的周期一般为5~10s。
设某一个时刻4节点系统的量测数据如表2所示,其中遥信量的量测值为1表示开关闭合,0表示开关断开,遥测量的量测值使用有名值:有功功率的单位是MW,无功功率的单位是MVAR,电压幅值的单位是kv,角度的单位是度。
表2 量测数据
步骤3:形成计算模型
开关、输电线、变压器、发电机、电容器以及负荷的端口被称为节点(node),遍历所有这些节点,根据得到的开关量将电网中通过闭合的开关或刀闸连接在一起的节点合并成一个计算用节点,以下所述节点均指计算用节点。
将第二步中得到的电压幅值量测值、发电机功率量测以及负荷功率量测重新分配到节点上,具体分配办法如下:1)如果节点上有多个电压幅值量测,则取平均值作为一个电压量测2)只有节点上所有的发电机和负荷都有功率量测时才将发电机功率减去负荷功率作为节点上注入功率量测进入下一步。
4节点系统中母线段1-1和1-2通过开关K2连接,K2闭合情况下,1-1和1-2合并成一个bus,该bus的电压幅值量测是1-1和1-2上电压幅值量测的平均值。计算模型如图5所示。量测重新分配后如表3所示。
表3 量测值
量测项目 | P<sub>1</sub> | Q<sub>1</sub> | P<sub>12</sub> | Q<sub>12</sub> | P<sub>13</sub> | Q<sub>13</sub> | P<sub>31</sub> | Q<sub>31</sub> |
量测值 | 18.71 | 34.05 | -41.79 | 37.96 | -19.1 | -2.24 | 17.55 | -10.74 |
量测项目 | P<sub>32</sub> | Q<sub>32</sub> | P<sub>3</sub> | Q<sub>3</sub> | P<sub>34</sub> | Q<sub>34</sub> | U<sub>1</sub> | U<sub>3</sub> |
量测值 | 132.57 | 42.61 | -49.49 | -41.93 | -193.2 | -71.08 | 112.15 | 110.99 |
步骤4:求解状态变量
设联络节点集合为A,发电机节点集合E,用Pkmax、Qkmax分别表示发电机节点k最大有功出力和最大无功出力,采用如下优化问题数据模型:
离线应用时可以用WLS估计器估计后得到的状态量作为初值,在线实时应用时可以使用上一次估计得到的状态量作为初值。使用非线性优化问题的数值求解方法求解,设得到的最优解为x*。
4节点系统取节点1作为电压相角参考节点,量测值矢量z以及参数αi,βi见表4。状态量计算结果如表5所示。
表4 量测值矢量及参数
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
z<sub>i</sub> | 18.71 | 34.05 | -41.79 | 37.96 | -19.1 | -2.24 | 17.55 | -10.74 |
α<sub>i</sub> | 0.777 | 1.1656 | 2.28 | 3.42 | 2.28 | 3.42 | 2.28 | 3.42 |
β<sub>i</sub> | 12.869 | 8.5796 | 4.38596 | 2.924 | 4.386 | 2.92398 | 4.386 | 2.92398 |
i | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
z<sub>i</sub> | 132.57 | 42.61 | -49.49 | -41.93 | -193.2 | -71.08 | 112.15 | 110.99 |
α<sub>i</sub> | 2.28 | 3.42 | 2.28 | 3.42 | 2.28 | 3.42 | 2.64 | 2.64 |
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
β<sub>i</sub> | 4.386 | 2.924 | 4.38596 | 2.924 | 4.386 | 2.92398 | 3.7879 | 3.78788 |
表5 估计结果
节点号 | 幅值 | 相角 |
1 | 111.31 | 0 |
2 | 110.17 | -0.398 |
3 | 111.52 | 0.82 |
4 | 112.1 | 7.92 |
步骤5:计算所有量测量的估计值
用上一步求得的x*计算出量测值的估计值z*=h(x*),并保存计算结果,求得的结果如表6所示。如果是在线实时计算,等待数秒,返回步骤2;如果是离线计算,整个计算过程结束。
表6 量测估计值
量测项目 | 量测值 | 估计值 |
P<sub>1</sub> | 18.71 | 20.33 |
Q<sub>1</sub> | 34.05 | 35.9 |
P<sub>12</sub> | -41.79 | 39.87 |
Q<sub>12</sub> | 37.96 | 38.3 |
P<sub>13</sub> | -19.1 | -19.5 |
Q<sub>13</sub> | -2.24 | -2.39 |
P<sub>31</sub> | 17.55 | 19.6 |
Q<sub>31</sub> | -10.74 | -11.23 |
P<sub>32</sub> | 132.57 | 130.2 |
Q<sub>32</sub> | 42.61 | 40.77 |
P<sub>3</sub> | -49.49 | -46.43 |
Q<sub>3</sub> | -41.93 | -41.74 |
P<sub>34</sub> | -193.22 | -196.27 |
Q<sub>34</sub> | -71.08 | -71.27 |
U<sub>1</sub> | 112.15 | 111.31 |
U<sub>3</sub> | 110.99 | 111.52 |
Claims (4)
1.基于数据采集系统量测数据的电力系统状态估计方法,其特征在于,所述方法是在计算机中依次按以下步骤实现的:
步骤(1),从能量管理系统中获得电力系统的网络参数,包括:输电线路的支路号、首节点和末节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器的变比和阻抗;
步骤(2),从数据采集和监视系统SCADA的网络数据服务器获得遥测量测值和遥信量测值,其中
遥测量测值包括:电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率、线路末端无功功率、线路电流幅值,
遥信量测值包括:开关、刀闸开合状态以及变压器分接头位置,开关或刀闸闭合的量测值为1,断开的量测值为0,
所述SCADA系统采用轮询方式采集数据,轮询周期为5~10s,每一个测点至少包括:测点类型、测点号以及量测值;
步骤(3),构造计算模型
步骤(3.1),把开关、输电线、变压器、发电机、电容器以及负荷的端口定义为节点,遍历所述这些节点,根据得到的开关量把电网中通过闭合的开关或刀闸连接在一起的节点合并成一个计算用节点,下面所述节点均指计算用节点,
步骤(3.2),把步骤(2)得到的电压幅值量测值、发电机功率量测值以及负荷功率量测值按以下方法分配到对应的节点上:
若节点上有多个电压幅值量测值,则取平均值作为一个电压幅值量测值,
若只有节点上所有的发电机和负荷都有功率量测时,才把发电机功率减去负荷功率作为节点上的注入功率;
步骤(4),以节点电压幅值与电压相角作为状态变量x,并取某一节点作为电压相角参考节点,便按下式求所述状态变量最优解x*:
步骤(4.1),把没有发电机也没有负荷的节点定义为联络节点,该联络节点的集合为A,而发电机节点集合为E,用Pkmax、Qkmax分别表示节点号为k的发电机节点最大有功出力和最大无功出力,
步骤(4.2),采用下述优化模型求解状态变量最优解x*:
其中,A为联络节点集合,E为发电机节点集合,
vt,vj,vk,vl分别为节点t、节点j、节点k、节点l的电压幅值,
j∈t表示节点j与节点t有线路相连,
l∈k表示节点l与节点k有线路相连,
节点t为联络节点,节点k为发电机节点,
Gtj,Btj,Gkl,Bkl为导纳矩阵元素,
θtj为节点t和节点j的电压相角差,
θkl为节点k和节点l的电压相角差,
zi为第i个量测的量测值,
hi(x)为第i个量测的量测函数,
αi是和第i个量测有关的参数:
对于电压测点:0.02×|zibase|,
对于有功测点:0.02×|zibase|,
对于无功测点:0.03×|zibase|,
zibase的大小和第i个量测有关,对于不同的电压等级,zibase是设定的,
βi是一个与所述αi相关的一个参数,
βi=c/αi,c是一个常数,c=5~10;
步骤(5),按照步骤(4)得到的最优解x*,计算出量测值的估计值z*=h(x*),
节点i电压幅值,h(x*)=vi *,
节点i有功注入功率,
节点i无功注入功率,
vi *,vj *,θi *,θj *均为状态变量最优解,用x*表示;
步骤(6),按下式判断每一个量测的估计值zi *是否合格:
|zi-zi *|≤αi,则估计值合格。
2.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法,其特征在于:对于线路功率量测,所述zibase按以下方式取值:500kv等级取1083MVA,330kv等级取686MVA,220kv等级取305MVA,110kv等级取114MVA,66kv电压等级取69.7MVA。
3.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法,其特征在于:对于电压幅值量测,所述zibase按以下方式取值:500kv等级取600kv,330kv等级取396kv,220kv等级取264kv,110kv等级取132kv,66kv等级取79.2kv。
4.根据权利要求1所述的基于数据采集系统量测数据的电力系统状态估计方法,其特征在于:对于发电机功率量测,所述zibase取发电机视在功率。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102237660A CN101383511B (zh) | 2008-10-10 | 2008-10-10 | 基于数据采集系统量测数据的电力系统状态估计方法 |
US12/475,223 US8150640B2 (en) | 2008-10-10 | 2009-05-29 | Method for state estimation based on the measured data of the data-acquisition system in power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102237660A CN101383511B (zh) | 2008-10-10 | 2008-10-10 | 基于数据采集系统量测数据的电力系统状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101383511A CN101383511A (zh) | 2009-03-11 |
CN101383511B true CN101383511B (zh) | 2010-08-04 |
Family
ID=40463198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102237660A Expired - Fee Related CN101383511B (zh) | 2008-10-10 | 2008-10-10 | 基于数据采集系统量测数据的电力系统状态估计方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8150640B2 (zh) |
CN (1) | CN101383511B (zh) |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599643B (zh) * | 2009-04-23 | 2011-04-27 | 清华大学 | 一种基于指数型目标函数的电力系统抗差状态估计方法 |
CN101635456B (zh) * | 2009-06-23 | 2013-07-10 | 上海申瑞继保电气有限公司 | 输配电网联合状态估计的方法 |
US8744638B2 (en) | 2009-09-11 | 2014-06-03 | General Electric Company | Method and system for demand response in a distribution network |
CN101661069B (zh) * | 2009-09-25 | 2011-07-20 | 北京四方继保自动化股份有限公司 | 不依赖状态矩阵的弱可观非pmu测点动态过程实时估计方法 |
CN101964525B (zh) * | 2010-06-25 | 2013-04-03 | 清华大学 | 一种支持大规模电流量测的配电网状态估计方法 |
CN101924393B (zh) * | 2010-08-11 | 2013-04-03 | 徐文远 | 一种电力系统负荷模型参数在线估计方法及其系统 |
CN102005758B (zh) * | 2010-11-16 | 2012-11-07 | 重庆大学 | 一种基于内网实测信息的两端口外网静态等值方法 |
CN102255310A (zh) * | 2011-07-15 | 2011-11-23 | 广东电网公司电力科学研究院 | Avc评估方法及系统 |
CN102495936B (zh) * | 2011-12-20 | 2013-04-03 | 东方电子股份有限公司 | 一种avc系统自动建模免维护方法 |
WO2013187975A1 (en) | 2012-06-15 | 2013-12-19 | Abb Research Ltd. | Parallel computation of dynamic state estimation for power system |
EP2874273B1 (en) * | 2012-07-10 | 2017-09-06 | Hitachi, Ltd. | Power system control system, and power system control method |
CN102801381B (zh) * | 2012-08-15 | 2014-10-29 | 上海交通大学 | 电压幅值与相角独立控制的可控变压器装置及其控制方法 |
CN102982379B (zh) * | 2012-11-16 | 2016-08-03 | 清华大学 | 基于集合论估计模型的电力系统状态估计方法及系统 |
CN102928655B (zh) * | 2012-11-26 | 2015-08-05 | 慈溪市供电局 | 一种电力异常数据检测方法 |
CN103107536B (zh) * | 2013-01-29 | 2014-11-05 | 中国海洋石油总公司 | 一种海上油田群电网的状态估计方法 |
CN103345570B (zh) * | 2013-06-05 | 2016-05-25 | 中国南方电网有限责任公司 | 基于抗差估计和拉格朗日乘子法的电力系统参数辨识方法 |
CN103605856B (zh) * | 2013-11-28 | 2016-08-31 | 东北石油大学 | 基于分数阶线路模型的输电线路参数估计方法 |
US20150160670A1 (en) * | 2013-12-09 | 2015-06-11 | Georgia Tech Research Corporation | Methods and systems for using distributed energy resources in an electric network |
CN104184144B (zh) * | 2014-09-05 | 2016-08-24 | 国家电网公司 | 一种用于多电压等级电网模型的抗差状态估计方法 |
CN104267298B (zh) * | 2014-10-23 | 2017-11-24 | 上海华群实业股份有限公司 | 一种双绕组电力变压器的参数计算方法 |
CN104330678B (zh) * | 2014-11-20 | 2017-08-11 | 国家电网公司 | 一种基于低压配电网三相状态估计的台区用电监测方法 |
CN105046038A (zh) * | 2015-04-03 | 2015-11-11 | 国家电网公司 | 一种线路电纳参数的估计方法 |
CN105046039A (zh) * | 2015-04-03 | 2015-11-11 | 国家电网公司 | 一种线路电阻参数的估计方法 |
CN104766248B (zh) * | 2015-04-16 | 2019-05-03 | 东南大学 | 一种有功坏数据辨识与修正方法 |
CN104809340A (zh) * | 2015-04-23 | 2015-07-29 | 国家电网公司 | 一种多端柔性直流输电系统运行点的快速计算方法 |
CN104899435A (zh) * | 2015-05-25 | 2015-09-09 | 清华大学 | 计及零注入约束的电力系统动态状态估计方法 |
JP6412841B2 (ja) * | 2015-09-10 | 2018-10-24 | 株式会社日立製作所 | 分散電源発電量推定装置および方法 |
FR3055481B1 (fr) * | 2016-08-24 | 2018-08-17 | Ifp Energies Now | Procede et systeme de gestion d'energie d'une supercapacite au moyen d'un modele de vieillissement et d'une prediction de la houle |
CN106443347B (zh) * | 2016-10-09 | 2018-05-15 | 中国矿业大学 | 一种基于故障暂态行波衰减分量的电网故障区段定位方法 |
CN106570774A (zh) * | 2016-10-11 | 2017-04-19 | 国网山东省电力公司烟台供电公司 | 一种基于学习理论的渐进学习电网调控方法 |
CN106789296B (zh) * | 2016-12-29 | 2019-07-05 | 山东鲁能智能技术有限公司 | 一种站所终端遥信与线路开关信号配置的处理系统及方法 |
CN106887844B (zh) * | 2017-04-20 | 2019-05-21 | 清华大学 | 一种基于稀疏恢复的快速分解配电网三相状态估计方法 |
CN107451333B (zh) * | 2017-06-30 | 2020-06-05 | 上海建工集团股份有限公司 | 一种预应力空间结构张力体系的内力测量方法 |
CN109873415B (zh) * | 2017-12-01 | 2020-12-25 | 南方电网科学研究院有限责任公司 | 一种电力系统等值方法和装置 |
CN108695862B (zh) * | 2018-04-25 | 2021-06-04 | 湖南大学 | 一种基于pmu实测数据的电网惯量特征在线评估方法 |
CN109167568B (zh) * | 2018-10-26 | 2019-09-27 | 江苏云上电力科技有限公司 | 一种基于冗余量测的光伏电站系统状态估计方法 |
CN109412857A (zh) * | 2018-11-12 | 2019-03-01 | 山东鲁能智能技术有限公司 | 一种泛型通用数据结构、网关机 |
CN109557390B (zh) * | 2018-11-16 | 2021-02-09 | 陕西千山航空电子有限责任公司 | 一种离散量输出接口的自检测电路 |
CN109639355B (zh) * | 2018-12-04 | 2022-07-01 | 马鞍山钢铁股份有限公司 | 一种工业光纤环网无扰接入ems在线系统的方法 |
CN109638823B (zh) * | 2018-12-19 | 2022-08-02 | 国网河南省电力公司电力科学研究院 | 用于辐射型电网状态估计的自动分区方法及装置 |
TWI693768B (zh) * | 2019-02-18 | 2020-05-11 | 英捷力電機股份有限公司 | 電力系統的電容器之電容值監控方法 |
CN110380409B (zh) * | 2019-07-16 | 2020-11-13 | 山东大学 | 考虑通信失败的主动配电网分布式鲁棒状态估计方法及系统 |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN110428185B (zh) * | 2019-08-08 | 2020-11-03 | 河海大学 | 基于伪量测模型的电-热互联综合能源系统抗差状态估计方法 |
CN110912267A (zh) * | 2019-11-22 | 2020-03-24 | 苏州浪潮智能科技有限公司 | 远程控制服务器机柜接入和离开电网的方法及装置 |
CN111314006B (zh) * | 2019-12-05 | 2022-01-28 | 国网辽宁省电力有限公司锦州供电公司 | 基于无时标越限数据传送的时间延时曲线拟合方法 |
CN112636325B (zh) * | 2020-11-11 | 2023-04-18 | 广东电网有限责任公司广州供电局 | 一种配电网抗差状态计算方法、装置、终端及存储介质 |
CN112653182B (zh) * | 2020-12-05 | 2023-07-28 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于极大似然估计优化清洁能源并网模型参数方法 |
CN112506909B (zh) * | 2020-12-11 | 2022-06-21 | 清华大学 | 配电系统测量数据补全方法及系统、电子设备及存储介质 |
CN112993989B (zh) * | 2021-03-05 | 2022-12-16 | 广东电网有限责任公司广州供电局 | 一种主动配电系统抗差状态估计数据处理方法 |
CN113139288A (zh) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | 基于gpu二次加速的配电网分布式抗差状态估计方法及装置 |
CN113505468B (zh) * | 2021-06-08 | 2022-12-13 | 西安交通大学 | 区间线性交流最大传输容量模型的构建和解析方法 |
CN113489072B (zh) * | 2021-07-05 | 2022-12-06 | 合肥工业大学 | 一种考虑支路功率伪量测的配电网状态估计方法 |
CN113746097B (zh) * | 2021-09-27 | 2023-09-29 | 杭州电力设备制造有限公司 | 一种电力系统分析方法、系统、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3903402A (en) * | 1971-08-26 | 1975-09-02 | Westinghouse Electric Corp | Digital computer program system employed in a hybrid loadflow computer arrangement for monitoring the security of an electric power system |
US7519506B2 (en) * | 2002-11-06 | 2009-04-14 | Antonio Trias | System and method for monitoring and managing electrical power transmission and distribution networks |
-
2008
- 2008-10-10 CN CN2008102237660A patent/CN101383511B/zh not_active Expired - Fee Related
-
2009
- 2009-05-29 US US12/475,223 patent/US8150640B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US20100094574A1 (en) | 2010-04-15 |
CN101383511A (zh) | 2009-03-11 |
US8150640B2 (en) | 2012-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101383511B (zh) | 基于数据采集系统量测数据的电力系统状态估计方法 | |
CN111026927B (zh) | 一种低压台区运行状态智能监测系统 | |
CN101924364B (zh) | 变电站-调度中心两级分布式电网的非线性状态估计方法 | |
CN101599643B (zh) | 一种基于指数型目标函数的电力系统抗差状态估计方法 | |
CN107453357A (zh) | 一种基于分层求解的配电网状态估计方法 | |
CN106026086B (zh) | 一种电网中运行状态的动态估计方法 | |
CN103020726A (zh) | 面向全pmu量测的抗差状态估计方法 | |
CN101615213B (zh) | 基于扩展不确定度的电力系统状态估计结果评价方法 | |
CN100554976C (zh) | 基于同步相量测量的区域电压稳定性监视方法 | |
CN106487010A (zh) | 一种配电网可靠性评估方法及系统 | |
CN113078630A (zh) | 一种基于实时量测数据的低压配电网拓扑辨识方法 | |
CN102590652B (zh) | 基于电气信息的设备性能评价系统及方法 | |
CN106229970A (zh) | 基于换流器控制特性的微电网状态估计方法 | |
CN103997041B (zh) | 一种交直流混联系统的在线电压稳定判别系统 | |
CN106159941B (zh) | 一种考虑实际量测误差传递特性的电力系统状态估计方法 | |
CN105139286A (zh) | 基于无功不平衡量的scada数据有效性评估方法 | |
CN114113885B (zh) | 一种基于冗余校验的异常低压分相线损精准定位方法 | |
CN104022501B (zh) | 基于模糊理论的配电网状态估计方法 | |
CN109901022A (zh) | 基于同步量测数据的配电网区域定位方法 | |
Liu et al. | PMU based voltage stability analysis for transmission corridors | |
CN112103956A (zh) | 基于智能电表动态测量点的配网状态估计方法 | |
Qin et al. | A Method for Distribution Line Impedance Calculation Based on Metering and Distribution Data Integration | |
CN107546742B (zh) | 一种日前计划潮流有功及电压计算误差分析方法 | |
Zhou et al. | Automatic identification method of topology structure of low-voltage distribution station area based on multi-point measurement data | |
CN106374462B (zh) | 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100804 |
|
CF01 | Termination of patent right due to non-payment of annual fee |