一种基于指数函数的电力系统鲁棒状态估计方法技术方案

本发明专利技术提供一种基于指数函数的电力系统鲁棒状态估计方法，首先根据电力系统拓扑结构和支路阻抗参数，计算得到节点导纳矩阵；根据传感器布局策略，在电力系统相应节点安装PMU，形成量测矩阵H；读取当前时刻节点量测值z(k)，节点量测值包括节点电压相量，支路电流相量，其中k表示采样时刻；构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计；本发明专利技术提供的方法，基于指数函数的目标函数，求偏导得到的权重值，能够根据量测残差的大小，自动调整，有助于抑制坏数据和非高斯噪声等影响，提高算法的鲁棒性和状态估计精度抗干扰能力强，鲁棒性好。鲁棒性好。鲁棒性好。

【技术实现步骤摘要】
一种基于指数函数的电力系统鲁棒状态估计方法


[0001]本专利技术涉及电力系统分析
，特别是指一种基于指数函数的电力系统鲁棒状态估计方法。

技术介绍

[0002]由于监测控制和数据采集(Supervisory Control And Data Acquisition，SCADA)系统或同步相量测量单元(Phasor Measurement Unit，PMU)采集到的原始量测值存在误差和不良数据，因此这些量测值需要通过电力系统状态估计算法剔除不良数据之后，才能为电网控制系统提供可靠准确的电网状态变量。因此，电力系统状态估计，为实现电网安全可靠运行提供的重要的支撑作用。
[0003]状态估计过程中通常将量测噪声模型假定成高斯分布，但此假设有时不符合实际情况。根据太平洋西北国家实验室的研究报告，PMU的电压相量和电流相量的量测误差甚至遵循非高斯分布。

技术实现思路

[0004]本专利技术的主要目的在于克服现有技术中的上述缺陷，提出一种基于指数函数的电力系统鲁棒状态估计方法，能够提高估计精度，提供更加准确的电网状态，且抗干扰能力强，鲁棒性好。
[0005]本专利技术采用如下技术方案：
[0006]一种基于指数函数的电力系统鲁棒状态估计方法，包括如下：
[0007]根据电力系统拓扑结构和支路阻抗参数，计算得到节点导纳矩阵；
[0008]根据传感器布局策略，在电力系统相应节点安装PMU，形成量测矩阵H；
[0009]读取当前时刻节点量测值z(k)，节点量测值包括节点电压相量，支路电流相量，其中k表示采样时刻；
[0010]构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计。
[0011]具体地，所述构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计，具体为：
[0012]已知量测z(k)，电力系统在第k时刻的状态x(k)与量测值存在如下关系式：
[0013]z(k)＝Hx(k)+v(k)
[0014]其中，H是量测矩阵，v(k)是量测噪声；
[0015]基于指数函数的鲁棒状态估计通过最小化如下所述目标函数求得：
[0016][0017]其中，i＝1,2,3....m，m是第k时刻所采集到的量测值个数；ρ
MCC
(e
i
(k))是：
[0018][0019]其中a
i
是门限值，e
i
＝z
i
‑
H
i
x是量测残差，σ
i
是噪声的标准差；
[0020]最小化上述目标函数，通过对ρ
MCC
求导数，使得导数为0，进而可推导得出状态估计值为：
[0021][0022]其中和分别为在时刻k的第u+1和u次迭代的状态估计值，w
u
是第u次计算得到的对角权重矩阵，e
u
(k)为在时刻k的第u次计算得到的量测残差。
[0023]若算法收敛，输出当前时刻节点状态估计值，跳转到下一时刻读取新量测值，否则继续迭代。
[0024]由上述对本专利技术的描述可知，与现有技术相比，本专利技术具有如下有益效果：
[0025]本专利技术提供一种基于指数函数的电力系统鲁棒状态估计方法，首先根据电力系统拓扑结构和支路阻抗参数，计算得到节点导纳矩阵；根据传感器布局策略，在电力系统相应节点安装PMU，形成量测矩阵H；读取当前时刻节点量测值z(k)，节点量测值包括节点电压相量，支路电流相量，其中k表示采样时刻；构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计；本专利技术提供的方法，基于指数函数的目标函数，求偏导得到的权重值，能够根据量测残差的大小，自动调整，有助于抑制坏数据和非高斯噪声等影响，提高算法的鲁棒性和状态估计精度抗干扰能力强，鲁棒性好。
附图说明
[0026]图1为本专利技术实施例提供的一种基于指数函数的电力系统鲁棒状态估计方法流程图；
[0027]图2是本专利技术实施例提供的利用本专利技术方法一节点的电压相角实部估计结果并与其它方法的对比结果图；
[0028]以下结合附图和具体实施例对本专利技术作进一步详述。
具体实施方式
[0029]本专利技术提供一种基于指数函数的电力系统鲁棒状态估计方法，基于指数函数的目标函数，求偏导得到的权重值，能够根据量测残差的大小，自动调整，有助于抑制坏数据和
非高斯噪声等影响，提高算法的鲁棒性和状态估计精度抗干扰能力强，鲁棒性好。
[0030]如图1，示出了一种基于指数函数的电力系统鲁棒状态估计方法流程图，包括如下：
[0031]S101：根据电力系统拓扑结构和支路阻抗参数，计算得到节点导纳矩阵；
[0032]对于各种微波网络，在选定的网络参考面上，定义出每个端口的电压和电流后，由于线性网络的电压和电流之间是线性关系，故选定不同的自变量和因变量，可以得到不同的线性组合。类似于低频双端口网络理论，这些不同变量的线性组合可以用不同的网络参数来表征，节点导纳矩阵就是微波网络中用端口电压(自变量)表示端口电流(因变量)的参量矩阵。因此，节点导纳矩阵的计算需要基于电力系统拓扑结构和支路阻抗参数
[0033]S102：根据传感器布局策略，在电力系统相应节点安装PMU，形成量测矩阵H；
[0034]PMU的安装节点位置会影响量测信息的准确性，现有技术中存在不少关于PMU的安装节点位置的优化算法，包括枚举法、启发式算法、搜索算法，来实现组合优化，这部分内容是不是本专利技术的创新点，因此不再赘述。
[0035]S103：读取当前时刻节点量测值z(k)，节点量测值包括节点电压相量，支路电流相量，其中k表示采样时刻；
[0036]S104：构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计。
[0037]具体地，所述构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计，具体为：
[0038]已知量测z(k)，电力系统在第k时刻的状态x(k)与量测值存在如下关系式：
[0039]z(k)＝Hx(k)+v(k)
[0040]其中，H是量测矩阵，v(k)是量测噪声；
[0041]基于指数函数的鲁棒状态估计通过最小化如下所述目标函数求得：
[0042][0043]其中，i＝1,2,3....m，m是第k时刻所采集到的量测值个数；ρ
MCC
(e
i
(k))是：
[0044][0045]其中a
i
是门限值，e
i
＝z
i
‑
H
i
x是量测残差，σ
i
是噪声的标准差；...

【技术保护点】

【技术特征摘要】
1.一种基于指数函数的电力系统鲁棒状态估计方法，其特征在于，包括如下：根据电力系统拓扑结构和支路阻抗参数，计算得到节点导纳矩阵；根据传感器布局策略，在电力系统相应节点安装PMU，形成量测矩阵H；读取当前时刻节点量测值z(k)，节点量测值包括节点电压相量，支路电流相量，其中k表示采样时刻；构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计。2.根据权利要求1所述的一种基于指数函数的电力系统鲁棒状态估计方法，其特征在于，所述构建基于指数函数的电力系统鲁棒状态估计算法，判断收敛条件，得到基于指数函数电力系统鲁棒状态的估计，具体为：已知量测z(k)，电力系统在第k时刻的状态x(k)与量测值存在如下关系式：z(k)＝Hx(k)+v(k)其中，H是量测矩阵，v(k)是量测噪声；基于指数函数的鲁棒状态估计通过最小化如下所述目标函数求得：其...

【专利技术属性】
技术研发人员：陈腾鹏，任和，李璐，朱真慧，高卓浩，王津聿，
申请(专利权)人：厦门大学，
类型：发明
国别省市：