泰宁新闻网

最新一种配电网数据平差方法及系统与流程

鞋

泰宁新闻网 http://www.tainingxinwen.cn 2020-10-18 01:28 出处:网络
这里介绍的最新一种配电网数据平差方法及系统与流程,为您描述接下来的东西:

这里介绍的最新一种配电网数据平差方法及系统与流程,为您描述接下来的东西:


本发明涉及电力系统领域,具体涉及一种配电网数据平差方法及系统。



背景技术:

近年来,生产生活领域对电力的需求逐年增长,配电网承担的负载越来越高,用户对供电质量的要求越来越严格。为了对电力系统的运行情况进行监控,配电网终端对电流、电压、开关位置等状态量进行了采集,因此配电网的数据质量关乎配电网的安全稳定运行,低质量的电力系统监测数据即使使用再精准的处理逻辑、再高效的并行策略都不会获得正确结果。

配电网量测系统数据主要取决于量测系统量测采样时产生的量测误差,其误差源主要包括:(1)测量仪器:一是测量表本身存在的误差,二是仪器精度不匹配产生的误差。(2)观测者或自读设备:观测者与自读设备的鉴别能力都有一定的局限性,所以在量测值读取甚至在设备安置时都会存在一定的误差,对观测结果有着重要的影响。(3)外界条件:观测时所处的外界条件,如温度、湿度、风力等都会对观测结果甚至被观测对象有一定的影响。因此只要是发生在工程实际中的观测量测都会存在这样或那样的误差,因此,配电网中测量产生的误差是不可避免的。

传统测量平差学对被测系统的属性要求较为苛刻,被测系统或对象需要具有无偏性,这意味着在多次测量之后量测量需要收敛为一个定值,换句话说也就是被测系统为一静态系统。但实际上配电网为一时变的动态系统,多次测量时其真值在将会根据实时电源出力、负荷变动情况而波动,无法统计大量的误差信息,没有足够的敏感度捕捉系统波动、暂降等信息,预测估计值不准确导致的滤波效果下降问题。



技术实现要素:

本发明提供的技术方案是:

一种配电网数据平差方法,包括:

获取历史量测数据和实时量测数据;

对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

对所述先验误差进行卡尔曼增益计算;

基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

优选的,所述历史数据包括:

节点电压幅值、相角差、支路有功、无功潮流量测、节点注入有功、无功量测、支路电流量测或节点注入电流。

优选的,所述对所述历史量测数据进行量测值预测得到先验估计值,包括:

将所述配电网历史数据通过自适应滤波处理;

将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差;

将所述滤出部分误差后的数据进行量测值预测得到先验估计值。

优选的,所述将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差后的数据如下式所示:

式中,σ为u(k)的标准差,u(k)为k时刻的配电网节点电压量测值,y(k)为自适应滤波的输出的处理后的数据信息;wt(k)为预估系数的转置;预估系数w(k)的计算式如下:

w(k+1)=w(k)+2g*e(k)*u(k)

式中,g为控制因子;e(k)为误差;

所述e(k)为误差的计算式如下:

e(k)=u(k-1)-y(k-1)。

优选的,所述先验估计的计算式如下:

x(k|k-1)=s(k)+b(k)

式中,s(k)为k时刻平滑后的值;b(k)为平滑后的趋势;

其中,所述s(k)的计算式如下:

s(k)=αhy(k)+(1-αh)x(k-1|k-2)

所述b(k)的计算式如下:

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)

式中,αh和βh为小于1的常数参数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计。

优选的,所述先验误差的计算式如下:

p(k|k-1)=bku(k)+qk

式中,p(k|k-1)为先验误差;u(k)为k时刻的配电网节点电压量测值;bk为状态转移矩阵,为先验估计计算式的泰勒展开一阶项,qk为系统误差阵。

优选的,所述卡尔曼增益的计算式如下:

式中:j(k)为卡尔曼增益;d△为量测误差的方差阵;为状态转移矩阵的转置。

优选的,所述后验估计的计算式如下:

x(k|k)=x(k|k-1)+j(k)m(k)

式中,x(k|k)为后验估计,m(k)为量测值与先验估计的差;

其中,所述量测值与先验估计的差m(k)的计算式如下:

m(k)=u(k)-bkx(k|k-1)

式中,u(k)为k时刻的配电网节点电压量测值。

优选的,所述后验方差的计算式如下:

p(k|k)=(1-j(k)bk)p(k|k-1)

其中,p(k|k)为后验方差。

一种配电网数据平差系统,包括:

获取模块:用于获取历史量测数据和实时量测数据;

预测模块:用于对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

第一计算模块:用于对所述先验误差进行卡尔曼增益计算;

第二计算模块:用于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

优选的,所述预测模块,包括:历史数据处理子模块和计算子模块;

所述历史数据处理子模块用于:

将所述配电网历史数据通过自适应滤波处理;

将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差;

将所述滤出部分误差后的数据进行量测值预测得到先验估计值;

所述计算子模块:用于计算先验估计和先验误差;

所述先验估计的计算式如下:

x(k|k-1)=s(k)+b(k)

式中,s(k)为k时刻平滑后的值;b(k)为平滑后的趋势;

其中,所述s(k)的计算式如下:

s(k)=αhy(k)+(1-αh)x(k-1|k-2)

所述b(k)的计算式如下:

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)

式中,αh和βh为小于1的常数参数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计;

所述先验误差的计算式如下:

p(k|k-1)=bku(k)+qk

式中,p(k|k-1)为先验误差;u(k)为k时刻的配电网节点电压量测值;bk为状态转移矩阵,为先验估计计算式的泰勒展开一阶项,qk为系统误差阵。

优选的,所述第一计算模块具体包括:

式中:j(k)为卡尔曼增益;d△为量测误差的方差阵;为状态转移矩阵的转置。

优选的,所述第二计算模块具体包括:

所述后验估计的计算式如下:

x(k|k)=x(k|k-1)+j(k)m(k)

式中,x(k|k)为后验估计,m(k)为量测值与先验估计的差;

其中,所述量测值与先验估计的差m(k)的计算式如下:

m(k)=u(k)-bkx(k|k-1)

式中,u(k)为k时刻的配电网节点电压量测值;

所述后验方差的计算式如下:

p(k|k)=(1-j(k)bk)p(k|k-1)

其中,p(k|k)为后验方差。

与现有技术相比,本发明的有益效果为:

本发明提供的技术方案一种配电网数据平差方法,包括:获取历史量测数据和实时量测数据;对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;对所述先验误差进行卡尔曼增益计算;基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。有效提高了预测估计值的准确性,解决了预测估计值不准确导致的滤波效果下降问题。

附图说明

图1为本发明一种配电网数据平差方法的流程示意图;

图2为本发明一种配电网数据平差方法的实施例具体流程图;

具体实施方式

为了更好地理解本发明,下面结合说明书附图和实例对本发明的内容做进一步的说明。

实施例1:

本发明提供的一种配电网数据平差方法如图1所示,包括:

s1、获取历史量测数据和实时量测数据;

s2、对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

s3、对所述先验误差进行卡尔曼增益计算;

s4、基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

如图2所示具体流程如下:

s1、获取历史量测数据和实时量测数据;

步骤1:获取配电网电压幅值量测值,并初始化状态变量。

状态变量的初始值的选择主要影响循环开始时的估计输出的准确性,然而,卡尔曼滤波较强的鲁棒性将使得整个滤波系统在初值较不精确的情况下仍能收敛到一个较为满意的结果。

假设有节点电压幅值的量测相量u1[u(k-n+1),u(k-n+2),...,u(k)],其中u(i)为第i个离散采样时间的电压量测值,我们需要根据前n-1个历史冗余量测量估计第k时间点状态量的真值。

初始化自适应滤波的控制因子g,两参数指数平滑法中的参数α和β。

s2、对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

步骤2:自适应滤波

自适应滤波器对输入的电压幅值量测数据进行第一次优化,流程图如图2所示。处理步骤遵循式(1)-(3)。由于仅使用历史数据,自适应滤波器不能对实时变化进行快速响应。因此,在式(3)中加入3σ原则判断,判断电压幅值突然跳变的异常值是误差还是系统波动。

w(k+1)=w(k)+2g*e(k)*l1(1)

e(k)=u(k-1)-y(k-1)(2)

其中w(k)为预估系数;g为控制因子;e(k)为误差,σ是u(k)的标准差,u(k)为k时刻的配电网节点电压量测值,y(k)为自适应滤波的输出。

步骤3:依据两参数指数平滑方法进行预测,得出先验估值及先验误差。

指数平滑法(也称线性外推法)是一种简单的短期负荷预测方法,具有存储变量少,计算速度快的优点,适合在线运算。预测函数主要依赖于两参数指数平滑方法,使得一方面预测系统接收相对精确的输入历史数据,另一方面仍然保持快速响应效果。

将步骤2的输出y(k)作为本步骤的输入进行处理,处理步骤遵循式(4)-(6)。

x(k|k-1)=s(k)+b(k)(4)

s(k)=αhy(k)+(1-αh)x(k-1|k-2)(5)

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)(6)

其中,参数αh和βh为小于1的常数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计;

之后计算先验误差

p(k|k-1)=bku(k)+qk(7)

这里的u(k)为量测相量u1中的u(k)。bk为状态转移矩阵,为状态方程(式(4)的泰勒展开一阶项,qk为系统误差阵,可认为是常数矩阵。

步骤4:计算后验估值及后验误差并更新

s3、对所述先验误差进行卡尔曼增益计算;

首先将步骤3计算得出的先验误差p(k|k-1)带入式中,计算卡尔曼增益,

其中,j(k)为卡尔曼增益,d△为量测误差的方差阵。

再进行新息计算

m(k)=u(k)-bkx(k|k-1)(9)

m(k)为新息,即量测值与预测值(先验估计)的差。

s4、基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

最后进行后验估计计算

x(k|k)=x(k|k-1)+j(k)m(k)(10)

其中,x(k|k)为后验估值,即对k时间点真值的估计值

以及后验方差计算

p(k|k)=(1-j(k)bk)p(k|k-1)(11)

其中p(k|k)为后验方差。

最终,配电网节点电压幅值平差的输出为式(10)及式(11)。分别表示对真值的估计值及估计值的误差。

同理,支路有功和无功潮流量测pij、qij;节点注入有功和无功量测pi、qi;

支路电流量测iij和节点注入电流ii等数据也可应用此方法进行平差与预测。

实施例2:

基于同一种发明构思,本发明还提供了一种配电网数据平差系统,包括:

获取模块:用于获取历史量测数据和实时量测数据;

预测模块:用于对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

第一计算模块:用于对所述先验误差进行卡尔曼增益计算;

第二计算模块:用于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

所述预测模块,包括:历史数据处理子模块和计算子模块;

所述历史数据处理子模块用于:

将所述配电网历史数据通过自适应滤波处理;

将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差;

将所述滤出部分误差后的数据进行量测值预测得到先验估计值;

所述计算子模块:用于计算先验估计和先验误差;

所述先验估计的计算式如下:

x(k|k-1)=s(k)+b(k)

式中,s(k)为k时刻平滑后的值;b(k)为平滑后的趋势;

其中,所述s(k)的计算式如下:

s(k)=αhy(k)+(1-αh)x(k-1|k-2)

所述b(k)的计算式如下:

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)

式中,αh和βh为小于1的常数参数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计;

所述先验误差的计算式如下:

p(k|k-1)=bku(k)+qk

式中,p(k|k-1)为先验误差;u(k)为k时刻的配电网节点电压量测值;bk为状态转移矩阵,为先验估计计算式的泰勒展开一阶项,qk为系统误差阵。

所述第一计算模块具体包括:

式中:j(k)为卡尔曼增益;d△为量测误差的方差阵;为状态转移矩阵的转置。

所述第二计算模块具体包括:

所述后验估计的计算式如下:

x(k|k)=x(k|k-1)+j(k)m(k)

式中,x(k|k)为后验估计,m(k)为量测值与先验估计的差;

其中,所述量测值与先验估计的差m(k)的计算式如下:

m(k)=u(k)-bkx(k|k-1)

式中,u(k)为k时刻的配电网节点电压量测值;

所述后验方差的计算式如下:

p(k|k)=(1-j(k)bk)p(k|k-1)

其中,p(k|k)为后验方差。

显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。

本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。


技术特征:

1.一种配电网数据平差方法,其特征在于,包括:

获取历史量测数据和实时量测数据;

对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

对所述先验误差进行卡尔曼增益计算;

基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

2.如权利要求1所述的一种配电网数据平差方法,其特征在于,所述历史数据包括:

节点电压幅值、相角差、支路有功、无功潮流量测、节点注入有功、无功量测、支路电流量测或节点注入电流。

3.如权利要求1所述的一种配电网数据平差方法,其特征在于,所述对所述历史量测数据进行量测值预测得到先验估计值,包括:

将所述配电网历史数据通过自适应滤波处理;

将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差;

将所述滤出部分误差后的数据进行量测值预测得到先验估计值。

4.如权利要求3所述的一种配电网数据平差方法,其特征在于,所述将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差后的数据如下式所示:

式中,σ为u(k)的标准差,u(k)为k时刻的配电网节点电压量测值,y(k)为自适应滤波的输出的处理后的数据信息;wt(k)为预估系数的转置;预估系数w(k)的计算式如下:

w(k+1)=w(k)+2g*e(k)*u(k)

式中,g为控制因子;e(k)为误差;

所述e(k)为误差的计算式如下:

e(k)=u(k-1)-y(k-1)。

5.如权利要求4所述的一种配电网数据平差方法,其特征在于,所述先验估计的计算式如下:

x(k|k-1)=s(k)+b(k)

式中,s(k)为k时刻平滑后的值;b(k)为平滑后的趋势;

其中,所述s(k)的计算式如下:

s(k)=αhy(k)+(1-αh)x(k-1|k-2)

所述b(k)的计算式如下:

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)

式中,αh和βh为小于1的常数参数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计。

6.如权利要求5所述的一种配电网数据平差方法,其特征在于,所述先验误差的计算式如下:

p(k|k-1)=bku(k)+qk

式中,p(k|k-1)为先验误差;u(k)为k时刻的配电网节点电压量测值;bk为状态转移矩阵,为先验估计计算式的泰勒展开一阶项,qk为系统误差阵。

7.如权利要求6所述的一种配电网数据平差方法,其特征在于,所述卡尔曼增益的计算式如下:

式中:j(k)为卡尔曼增益;d△为量测误差的方差阵;为状态转移矩阵的转置。

8.如权利要求7所述的一种配电网数据平差方法,其特征在于,所述后验估计的计算式如下:

x(k|k)=x(k|k-1)+j(k)m(k)

式中,x(k|k)为后验估计,m(k)为量测值与先验估计的差;

其中,所述量测值与先验估计的差m(k)的计算式如下:

m(k)=u(k)-bkx(k|k-1)

式中,u(k)为k时刻的配电网节点电压量测值。

9.如权利要求7所述的一种配电网数据平差方法,其特征在于,所述后验方差的计算式如下:

p(k|k)=(1-j(k)bk)p(k|k-1)

其中,p(k|k)为后验方差。

10.一种配电网数据平差系统,其特征在于,包括:

获取模块:用于获取历史量测数据和实时量测数据;

预测模块:用于对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;

第一计算模块:用于对所述先验误差进行卡尔曼增益计算;

第二计算模块:用于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。

11.如权利要求10所述的一种配电网数据平差系统,其特征在于,所述预测模块,包括:历史数据处理子模块和计算子模块;

所述历史数据处理子模块用于:

将所述配电网历史数据通过自适应滤波处理;

将通自适应滤波处理的配电网数据信息加入对统计误差的判断滤出部分误差;

将所述滤出部分误差后的数据进行量测值预测得到先验估计值;

所述计算子模块:用于计算先验估计和先验误差;

所述先验估计的计算式如下:

x(k|k-1)=s(k)+b(k)

式中,s(k)为k时刻平滑后的值;b(k)为平滑后的趋势;

其中,所述s(k)的计算式如下:

s(k)=αhy(k)+(1-αh)x(k-1|k-2)

所述b(k)的计算式如下:

b(k)=βh(s(k)-s(k-1))+(1-βh)b(k-1)

式中,αh和βh为小于1的常数参数;x(k|k-1)为利用k-1时刻及k-1时刻之前的信息对k时刻的先验估计;

所述先验误差的计算式如下:

p(k|k-1)=bku(k)+qk

式中,p(k|k-1)为先验误差;u(k)为k时刻的配电网节点电压量测值;bk为状态转移矩阵,为先验估计计算式的泰勒展开一阶项,qk为系统误差阵。

12.如权利要求10所述的一种配电网数据平差系统,其特征在于,所述第一计算模块具体包括:

式中:j(k)为卡尔曼增益;d△为量测误差的方差阵;为状态转移矩阵的转置。

13.如权利要求10所述的一种配电网数据平差系统,其特征在于,所述第二计算模块具体包括:

所述后验估计的计算式如下:

x(k|k)=x(k|k-1)+j(k)m(k)

式中,x(k|k)为后验估计,m(k)为量测值与先验估计的差;

其中,所述量测值与先验估计的差m(k)的计算式如下:

m(k)=u(k)-bkx(k|k-1)

式中,u(k)为k时刻的配电网节点电压量测值;

所述后验方差的计算式如下:

p(k|k)=(1-j(k)bk)p(k|k-1)

其中,p(k|k)为后验方差。

技术总结
本发明提供的技术方案一种配电网数据平差方法,包括:获取历史量测数据和实时量测数据;对所述历史量测数据进行量测值预测得到先验估计值,并基于实时量测数据得到先验误差;对所述先验误差进行卡尔曼增益计算;基于卡尔曼增益和先验估计、先验误差计算后验估计和后验方差。有效提高了预测估计值的准确性,解决了预测估计值不准确导致的滤波效果下降问题。

技术研发人员:盛万兴;胡丽娟;王思博;刘科研;林挚;董伟杰;贾东梨;刁赢龙;何开元;叶学顺;白牧可;张淼;吕琛;张稳;刘杨涛;张怀天
受保护的技术使用者:中国电力科学研究院有限公司;国家电网有限公司;国网北京市电力公司
技术研发日:2019.11.21
技术公布日:2020.06.26

最新一种配电网数据平差方法及系统与流程的相关内容如下:

本文标题:最新一种配电网数据平差方法及系统与流程
http://www.tainingxinwen.cn/qitaxinxi/520675.html

0

精彩评论

暂无评论...
验证码 换一张
取 消