CN107565553A - 一种基于ukf的配电网抗差动态状态估计方法 - Google Patents
一种基于ukf的配电网抗差动态状态估计方法 Download PDFInfo
- Publication number
- CN107565553A CN107565553A CN201710844451.7A CN201710844451A CN107565553A CN 107565553 A CN107565553 A CN 107565553A CN 201710844451 A CN201710844451 A CN 201710844451A CN 107565553 A CN107565553 A CN 107565553A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- mover
- mtr
- 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
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于UKF的配电网抗差动态状态估计方法,该方法将指数权函数引入到无迹卡尔曼滤波动态状态估计中,在***负荷突变引起量测误差增大的情况下,根据量测残差实现对量测权重的动态修正,进而平滑的调整卡尔曼滤波增益;同时由于引入的指数权函数是在最小二乘法中使用,为了使其在无迹卡尔曼滤波中便于实现,给出了具体的改进方法;通过采用本发明公开的方法,基于SCADA/PMU混合量测***下,对于量测存在粗差或***负荷突变时,能够有效的增强算法的抗差性能,提升UKF算法的滤波性能,提高了状态估计精度。
Description
技术领域
本发明涉及一种状态估计,具体涉及一种基于UKF的配电网抗差动态状态估计方法。
背景技术
随着欧美国家开展对智能电网的建设,中国相应地依据自己的国情,分别以特高压为骨干的输电网、构建自动化程度高的配电网作为建设方向。对配电网自动化的研究和发展,影响着智能配电网建设进度,研究高精度实时的动态状态估计对于配网管理***尤其重要。
大量分布式电源以及负荷特性变化大的电动汽车接入配电网,使配网的动态性更强。电网要求配电管理***提高监测与控制能力,为此,需要状态估计提供实时可靠完整的***运行数据。配电网三相不平衡、大量短支路的存在以及量测配置较少等特征,使在输电网应用成熟的状态估计算法不能应用在配电网中,如PQ解耦法,故研究适用于配电网状态估计的算法。
目前对配电网状态估计研究主要分为两类,一类基于最小二乘法的静态状态估计算法,以一个时间断面的量测数据进行状态估计,常使用的有基于节点电压、支路电流、支路功率作为状态变量,使***量测信息和量测估计值误差最小为收敛原则,要求量测严格服从正太分布,条件苛刻。用牛顿迭代法对状态模型求解。
另一类则是以卡尔曼滤波为基础的动态状态估计算法,考虑多个时间量测断面的量测数据进行状态估计,滤波效果好,计算速度快,而且还能预测下一时刻的***状态。在配电网量测量较少的情况下,基于两参数指数平滑法可以一步预测,补充***的伪量测。
无迹卡尔曼滤波(Unscented Kalman Filter,UKF)是以线性卡尔曼滤波理论为框架,采样无迹变换技术对非线性函数线性处理,近似程度达到四阶及以上,适用性好。而且该算法采用sigma点采样,不用计算雅克比矩阵,减少了扩展卡尔曼滤波算法的线性化误差,但也存在量测有粗差或负荷突变时,滤波性能下降,还有可能发散等缺点。
发明内容
为了解决上述UKF算法存在的不足,本发明提供了一种基于UKF的配电网抗差动态状态估计方法,以配电网三相电压幅值和相角作为状态变量,通过对称比例采样方法计算sigma点及其权值,然后利用两参数指数平滑法进行参数识别并进行一步预测,通过量测方程利用预测的状态量求取量测预测,对实时量测与量测预测值形成的新息矩阵,采用所提出的指数权函数动态修正,及时调整不良数据的权重,求取卡尔曼滤波增益,最后更新状态变量以及协方差。实验结果表明,当***正常运行时,基于指数权函数UKF算法与传统UKF算法滤波效果和时间相差不大;在***负荷突变或存在粗差时,UKF算法滤波效果极具下降且有发散的趋势,而指数权函数UKF法则不受此影响,估计精度较高,验证了改进算法的抗差性以及数值稳定性。
为了实现上述发明目的,本发明采取如下的技术方案:
提供一种基于UKF的配电网抗差动态状态估计方法,所述方法包括以下步骤:
步骤(1)输入配电网三相网络参数,包括支路参数,变压器参数,负荷参数,计算***节点导纳矩阵Y。
步骤(2)在K=0时刻,初始化状态变量节点电压的均值和协方差。
步骤(3)根据K时刻的状态估计值和协方差进行sigma点采样并计算采样点对应的权值。
步骤(4)采用Holt's两参数指数平滑法进行参数识别,得到状态变量在K时刻的预测值并利用量测方程计算得到量测预测值。
步骤(5)由K时刻量测预测值计算***的协方差及互协方差。
步骤(6)根据新息矩阵利用指数权函数对他平滑,动态修正K时刻量测值的量测权重,及时更新量测误差协方差矩阵。
步骤(7)计算卡尔曼滤波增益并更新***的状态和协方差。
步骤(8)判断***采样是否结束,若结束,执行步骤(7),若不结束,返回执行步骤(2)。
步骤(9)输出K+1时刻的状态变量估计值。
所述步骤(1)对配电网可观性进行分析,生成***的节点导纳矩阵Y。
所述步骤(2)配电网状态变量的均值和协方差作为UKF算法的初始值,其中p=abc。
所述步骤(3)在进行sigma点采样时,方法主要有单行采样、球形采样、对称比例采样。采样方式影响状态估计算法的计算效率和数值稳定性。单行采样速度最快,数值稳定性最差,球形采样的所有权系数都相同则会受状态初值影响,比例对称采样方法近似精度高。对于n维的配电***,采用对称比例采样方法,计算2n+1个sigma点和相对应的权值;选取对称采样法具有较好的数值稳定性,通过无迹变换得到2n+1个sigma采样点集及相应的权值:
参数λ=α2(n+κ)-n用来降低总的预测误差的尺度因子,参数α控制了采样点的分布状态,取值范围[0.0001,1],κ是确保矩阵(n+λ)P为半正定自由参数,β≥0可以把高阶项的影响包含在内的权系数,保证了近似精度。式(3)
中,P是状态量协方差矩阵,表示矩阵方根的第i列。
计算这些采样点相应的权值
Wm表示状态量的权值,Wc表示协方差的权值。
所述步骤(4)利用两参数指数平滑参数法识别,进行状态预测,计算速度快,精度在合理范围内。
所述步骤(5)通过加权求取状态变量的预测值和预测协方差:
其中,Q是******噪声协方差矩阵。
所述步骤(6)由sigma点集的量测预测值,通过加权求和得到***的均值和协方差、互协方差:
所述步骤(7):在配电网量测数据中,可能存在粗差,***负荷突变引起量测数据变化大,引入指数权函数对新息矩阵动态修正,改变相应量测权重,平滑了粗差或量测变化对滤波的影响,具体实现如下:
为了提高在UKF中的平滑效果,修改权函数为负指数形式:
R'k=Rkexp(-|yk-h(xk)|)
所述步骤(8):根据修正的量测权重计算量测误差方差阵,并计算卡尔曼滤波增益和更新***的状态和协方差:
其中K为卡尔曼滤波增益矩阵。
与现有技术相比,本发明的有益效果在于:
1.该方法在原有算法中实现容易,对于***负荷突变、量测存在粗差时,滤波效果好,抗差能力强,提高了状态估计精度。
2.数值稳定强,收敛效果好,能有效解决三相不平衡配电网状态估计问题。
附图说明
图1为本发明一种基于UKF的配电网抗差动态状态估计方法流程图;
图2为***正常运行时状态估计,本发明算法和传统的UKF算法比较结果;
图3为***负荷突变时的状态估计比较结果。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,但不作为对本发明的任何限制。
本发明的实施步骤如下:
步骤(1):根据配电网给定的网络接线、支路参数以及量测数据,其中网络接线如附图1,支路参数见表一,分析***的可观性并计算节点导纳矩阵。
步骤(2):在采样K=0时刻,初始化状态变量节点电压的均值和协方差。
步骤(3)根据K的状态估计值和协方差进行sigma点采样并计算采样点对应的权值。
步骤(4)采用Holt's两参数指数平滑法进行参数识别,得到状态变量在K时刻的预测值并利用量测方程计算得到量测预测值。
步骤(5)由K时刻量测预测值计算***的协方差及互协方差。
步骤(6)根据新息矩阵利用指数权函数对他平滑,动态修正K时刻量测值的量测权重,及时更新量测误差协方差矩阵。
步骤(7)计算卡尔曼滤波增益并更新***的状态和协方差。
步骤(8)判断***采样是否结束,若结束,执行步骤(7),若不结束,返回执行步骤(2)。
步骤(9)输出K+1时刻的状态变量估计值。
所述步骤(1)对配电网可观性进行分析,生成***的节点导纳矩阵Y。
所述步骤(2)配电网状态变量的均值和协方差作为UKF算法的初始值,其中p=abc。
所述步骤(3)在进行sigma点采样时,方法主要有单行采样、球形采样、对称比例采样。采样方式影响状态估计算法的计算效率和数值稳定性。单行采样速度最快,数值稳定性最差,球形采样的所有权系数都相同则会受状态初值影响,比例对称采样方法近似精度高。对于n维的配电***,采用对称比例采样方法,计算2n+1个sigma点和相对应的权值;选取对称采样法具有较好的数值稳定性,通过无迹变换得到2n+1个sigma采样点集及相应的权值:
参数设置:α=0.001,κ=0,β=2。
Wm表示状态量的权值,Wc表示协方差的权值。
所述步骤(4)利用两参数指数平滑参数法识别,进行状态预测,计算速度快,精度在合理范围内。
参数设置为:αH=0.85,βH=0.05
所述步骤(5)通过加权求取状态变量的预测值和预测协方差:
所述步骤(6)由sigma点集的量测预测值,通过加权求和得到***的均值和协方差、互协方差:
所述步骤(7):在配电网量测数据中,可能存在粗差,***负荷突变引起量测数据变化大,引入指数权函数对新息矩阵动态修正,改变相应量测权重,平滑了粗差或量测变化对滤波的影响,具体实现如下:
为了提高在UKF中的平滑效果,修改权函数为负指数形式:
R'k=Rk exp(-|yk-h(xk)|)
所述步骤(8):根据修正的量测权重计算量测误差方差阵,并计算卡尔曼滤波增益和更新***的状态和协方差:
其中K为卡尔曼滤波增益矩阵。
运用本发明的方法对IEEE33节点进行计算机仿真,定义滤波的绝对误差平均值和最大值:
其中Nbus为***节点数,k为采样时刻,上标e和t分别表示状态量的估计值和真值。
表1本发明与UKF算法在***正常运行时比较节点电压的相角和幅值的绝对误差平均值和最大值。
表2本发明与UKF算法在***负荷突变时比较节点电压的相角和幅值的绝对误差平均值和最大值。
从表1可以看出,在配电网***运行正常时,本发明的算法和UKF算法状态估计精度相近,略优于UKF算法;从表2可以看出,在***负荷突变时,本发明算法具有较高的滤波性能,而且数值稳定性好。
当然,以上只是本发明的具体应用范例,本发明还有其他的实施方式,凡采用等同替换或等效变换形成的技术方案,均落在本发明所要求的保护范围之内。
Claims (9)
1.一种基于UKF的配电网抗差动态状态估计方法,其特征在于包括如下步骤:
(1)输入配电网三相网络参数,包括支路参数,变压器参数,负荷参数,计算***节点导纳矩阵;
(2)在K=0时刻,初始化状态变量节点电压的均值和协方差;
(3)根据K时刻的状态估计值和协方差进行sigma点采样并计算sigma点相对应的权值;
(4)采用Holt's两参数指数平滑法进行参数识别,得到状态变量在K时刻的预测值并利用量测方程计算得到量测预测值;
(5)由K时刻量测预测值计算***的协方差及互协方差;
(6)利用指数权函数对残差平滑,动态修正K时刻量测值的量测权重,及时更新量测误差协方差矩阵;
(7)计算卡尔曼滤波增益并更新***的状态和协方差;
(8)判断***采样是否结束,若结束,执行步骤(9),若不结束,返回执行步骤(2);
(9)输出K+1时刻的状态变量估计值。
2.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:步骤(1)中,对配电网可观性进行分析,生成***的节点导纳矩阵。
3.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:初始化配电网状态变量的均值和协方差作为UKF算法的初始值;
计算方法如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>0</mn>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msubsup>
<mi>x</mi>
<mn>0</mn>
<mi>p</mi>
</msubsup>
<mo>&rsqb;</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>P</mi>
<mn>0</mn>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mo>(</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>0</mn>
<mi>P</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mn>0</mn>
<mi>p</mi>
</msubsup>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>0</mn>
<mi>P</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mn>0</mn>
<mi>p</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&rsqb;</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
其中p=abc表示三相。
4.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:在进行sigma点采样时,方法包括单行采样、球形采样、对称比例采样;
对于n维的配电***,采用对称比例采样方法,计算2n+1个sigma点和相对应的权值:
<mrow>
<msubsup>
<mi>&chi;</mi>
<mi>i</mi>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>p</mi>
</msup>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>p</mi>
</msup>
<mo>+</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msqrt>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
<msub>
<mi>P</mi>
<mi>x</mi>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>p</mi>
</msup>
<mo>-</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msqrt>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
<msub>
<mi>P</mi>
<mi>x</mi>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>2</mn>
<mi>n</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
参数λ=α2(n+κ)-n用来降低总的预测误差的尺度因子,参数α控制了采样点的分布状态,取值范围[0.0001,1],κ是确保矩阵(n+λ)P为半正定自由参数,β≥0可以把高阶项的影响包含在内的权系数,保证了近似精度;
表示矩阵方根的第i列;
计算这些采样点相应的权值:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>W</mi>
<mn>0</mn>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>/</mo>
<mi>n</mi>
<mo>+</mo>
<mi>&lambda;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>W</mi>
<mn>0</mn>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>/</mo>
<mi>n</mi>
<mo>+</mo>
<mi>&lambda;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>2</mn>
<mi>n</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
Wm表示状态量的权值,Wc表示协方差的权值。
5.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:利用两参数指数平滑参数法识别,进行状态预测,计算速度快,精度在合理范围内;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>H</mi>
</msub>
<msubsup>
<mi>x</mi>
<mi>k</mi>
<mi>p</mi>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mi>H</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&beta;</mi>
<mi>H</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mi>H</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
6.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:通过加权求取状态变量的预测值和预测协方差:
<mrow>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>L</mi>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>L</mi>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mi>k</mi>
<mi>p</mi>
</msubsup>
</mrow>
其中,Q是***噪声协方差矩阵。
7.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:由sigma点集的量测预测值,通过加权求和得到***的均值和协方差、互协方差:
<mrow>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>.</mo>
</mrow>
8.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:在配电网量测数据中,可能存在粗差,***负荷突变引起量测数据变化大,引入指数权函数对新息矩阵动态修正,改变相应量测权重,平滑粗差或量测变化对滤波的影响,具体实现如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>=</mo>
<mfrac>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
为了提高在UKF中的平滑效果,修改权函数为负指数形式:
R'k=Rk exp(-|yk-h(xk)|)
相应的协方差修改为:
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</munderover>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msubsup>
<msup>
<mi>R</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</msubsup>
<mo>.</mo>
</mrow>
9.根据权利要求1所述的基于UKF的配电网抗差动态状态估计方法,其特征在于:根据修正的量测权重计算量测误差方差阵,并计算卡尔曼滤波增益和更新***的状态和协方差:
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>|</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>K</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
</mrow>
其中K为卡尔曼滤波增益矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710844451.7A CN107565553A (zh) | 2017-09-19 | 2017-09-19 | 一种基于ukf的配电网抗差动态状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710844451.7A CN107565553A (zh) | 2017-09-19 | 2017-09-19 | 一种基于ukf的配电网抗差动态状态估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107565553A true CN107565553A (zh) | 2018-01-09 |
Family
ID=60980207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710844451.7A Pending CN107565553A (zh) | 2017-09-19 | 2017-09-19 | 一种基于ukf的配电网抗差动态状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107565553A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108711885A (zh) * | 2018-06-05 | 2018-10-26 | 重庆大学 | 一种针对风力场状态协同估计方法 |
CN109376910A (zh) * | 2018-09-28 | 2019-02-22 | 河海大学 | 一种基于历史数据驱动的配电网动态状态估计方法 |
CN109754013A (zh) * | 2018-12-31 | 2019-05-14 | 浙江大学 | 一种基于无迹卡尔曼滤波的电力***混合量测融合方法 |
CN110247396A (zh) * | 2019-07-17 | 2019-09-17 | 国网山东省电力公司青岛供电公司 | 基于自适应抗差无迹卡尔曼滤波的配电网状态估计方法及*** |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN111695082A (zh) * | 2020-06-30 | 2020-09-22 | 上海交通大学 | 智能配电网抗差动态状态估计方法 |
CN111723366A (zh) * | 2019-03-19 | 2020-09-29 | 中国科学院沈阳自动化研究所 | 一种基于状态估计偏差的错误数据注入攻击鲁棒检测方法 |
CN112383048A (zh) * | 2020-10-20 | 2021-02-19 | 重庆大学 | 考虑公专用配电变压器测量特点的中压馈线电网三相状态估计方法 |
CN113313339A (zh) * | 2021-03-26 | 2021-08-27 | 贵州大学 | 一种配电网动态状态估计算法 |
CN114442557A (zh) * | 2022-01-25 | 2022-05-06 | 西南交通大学 | 一种机床温度场快速辨识方法及*** |
CN115616333A (zh) * | 2022-12-20 | 2023-01-17 | 国网江西省电力有限公司电力科学研究院 | 一种配电网线损预测方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615794A (zh) * | 2009-08-05 | 2009-12-30 | 河海大学 | 基于无迹变换卡尔曼滤波的电力***动态状态估计方法 |
CN103529424A (zh) * | 2013-10-23 | 2014-01-22 | 北京工商大学 | 一种基于rfid及ukf实现室内目标快速跟踪的方法 |
CN104777426A (zh) * | 2015-04-17 | 2015-07-15 | 河海大学 | 一种基于无迹变换强跟踪的发电机动态状态估计方法 |
CN106707235A (zh) * | 2017-03-08 | 2017-05-24 | 南京信息工程大学 | 一种基于改进的无迹卡尔曼滤波的室内测距定位方法 |
-
2017
- 2017-09-19 CN CN201710844451.7A patent/CN107565553A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615794A (zh) * | 2009-08-05 | 2009-12-30 | 河海大学 | 基于无迹变换卡尔曼滤波的电力***动态状态估计方法 |
CN103529424A (zh) * | 2013-10-23 | 2014-01-22 | 北京工商大学 | 一种基于rfid及ukf实现室内目标快速跟踪的方法 |
CN104777426A (zh) * | 2015-04-17 | 2015-07-15 | 河海大学 | 一种基于无迹变换强跟踪的发电机动态状态估计方法 |
CN106707235A (zh) * | 2017-03-08 | 2017-05-24 | 南京信息工程大学 | 一种基于改进的无迹卡尔曼滤波的室内测距定位方法 |
Non-Patent Citations (3)
Title |
---|
蔡凝露等: "基于指数权函数的抗差状态估计算法", 《中国电力》 * |
赵洪山、田甜: "基于自适应无迹卡尔曼滤波的电力***动态状态估计", 《电网技术》 * |
高丽君、刘济: "基于Cauchy鲁棒函数的UKF改进算法", 《华东理工大学学报(自然科学版)》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108711885A (zh) * | 2018-06-05 | 2018-10-26 | 重庆大学 | 一种针对风力场状态协同估计方法 |
CN109376910A (zh) * | 2018-09-28 | 2019-02-22 | 河海大学 | 一种基于历史数据驱动的配电网动态状态估计方法 |
CN109376910B (zh) * | 2018-09-28 | 2021-08-31 | 河海大学 | 一种基于历史数据驱动的配电网动态状态估计方法 |
CN109754013A (zh) * | 2018-12-31 | 2019-05-14 | 浙江大学 | 一种基于无迹卡尔曼滤波的电力***混合量测融合方法 |
CN111723366A (zh) * | 2019-03-19 | 2020-09-29 | 中国科学院沈阳自动化研究所 | 一种基于状态估计偏差的错误数据注入攻击鲁棒检测方法 |
CN110247396A (zh) * | 2019-07-17 | 2019-09-17 | 国网山东省电力公司青岛供电公司 | 基于自适应抗差无迹卡尔曼滤波的配电网状态估计方法及*** |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN111695082A (zh) * | 2020-06-30 | 2020-09-22 | 上海交通大学 | 智能配电网抗差动态状态估计方法 |
CN112383048A (zh) * | 2020-10-20 | 2021-02-19 | 重庆大学 | 考虑公专用配电变压器测量特点的中压馈线电网三相状态估计方法 |
CN113313339A (zh) * | 2021-03-26 | 2021-08-27 | 贵州大学 | 一种配电网动态状态估计算法 |
CN114442557A (zh) * | 2022-01-25 | 2022-05-06 | 西南交通大学 | 一种机床温度场快速辨识方法及*** |
CN114442557B (zh) * | 2022-01-25 | 2023-05-12 | 西南交通大学 | 一种机床温度场快速辨识方法及*** |
CN115616333A (zh) * | 2022-12-20 | 2023-01-17 | 国网江西省电力有限公司电力科学研究院 | 一种配电网线损预测方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107565553A (zh) | 一种基于ukf的配电网抗差动态状态估计方法 | |
CN108155648A (zh) | 基于自适应h无穷扩展卡尔曼滤波的状态估计方法 | |
CN107590317B (zh) | 一种计及模型参数不确定性的发电机动态估计方法 | |
CN106443285B (zh) | 基于总体最小二乘法的多谐波源谐波责任量化分析方法 | |
CN109818349B (zh) | 一种基于多维状态矩阵滑动匹配的电网鲁棒状态预测方法 | |
CN107145707B (zh) | 一种计及光伏出力不确定性和全寿命周期成本的配电网变压器规划方法 | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
CN110247396A (zh) | 基于自适应抗差无迹卡尔曼滤波的配电网状态估计方法及*** | |
CN112630659A (zh) | 一种基于改进bp-ekf算法的锂电池soc估算方法 | |
CN103177188A (zh) | 一种基于特征映射的电力***负荷动特性分类方法 | |
CN111460374A (zh) | 一种考虑节点差异性的配电网d-pmu优化配置方法 | |
CN112149873A (zh) | 一种基于深度学习的低压台区线损合理区间预测方法 | |
CN113177600B (zh) | 一种电力***自适应鲁棒状态估计方法 | |
CN105184027A (zh) | 一种基于交互式多模型算法的电力负荷建模方法 | |
CN108054757A (zh) | 一种内嵌无功和电压的n-1闭环安全校核方法 | |
CN115081299A (zh) | 一种基于upf的电力***鲁棒辅助预测状态估计方法 | |
CN111709350B (zh) | 基于fcm聚类的低频振荡模态参数识别方法及*** | |
CN110783918A (zh) | 一种基于线性模型的配电三相区间状态估计求解算法 | |
CN111581768A (zh) | 一种基于混合量测的配电网分布式状态估计方法 | |
CN109218073A (zh) | 一种计及网络攻击和参数不确定性的动态状态估计方法 | |
CN106372440B (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
CN106443253A (zh) | 一种基于pmu数据的输电线路参数辨识方法 | |
CN110021931B (zh) | 一种计及模型不确定性的电力***辅助预测状态估计方法 | |
CN106599541B (zh) | 一种动态电力负荷模型的结构和参数在线辨识方法 | |
CN112670981A (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: 20180109 |