[发明专利]一种基于极大似然估计的配电网状态估计方法

权利要求书

1.一种基于极大似然估计的配电网状态估计方法，其特征在于，该方法包括以下步骤：

(1)从配电网的配电管理系统数据库中获取系统状态变量的实时量测值从历史数据库中获取同一时刻的配电网负荷，利用短期负荷预测方法生成负荷伪量测值其中j为配电网的节点编号，通过下式计算量测误差：

或

其中，v_j为节点j的量测误差，为节点j的量测真值，由已知配电网节点的注入功率和部分电压值，通过潮流计算得到潮流分布，作为量测真值；

(2)对于步骤(1)中得到的伪量测误差v_j，计算伪量测误差v_j概率密度函数为：

$pdf\left(v_j\right)=\frac{1}{{\mathrm{nd}}_j}\underset{i=1}{\overset{n}{\Σ}}K\left(\frac{v_j-v_j^i}{d_j}\right)$

其中，是上述v_j的第i个采样值，i＝1,2,...,n，d_j为节点j核密度估计函数的宽带，σ_j是节点j核密度估计函数的标准差，K(.)是标准高斯核函数，定义为：

其中，y为标准高斯核函数的自变量取值，R为实数集合；

(3)或设步骤(1)中的实时量测满足正态分布，则实时量测误差的概率密度函数为：

$pdf\left(v_j\right)=\frac{1}{{\mathrm{\σ}}_j\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2}$

其中，x为状态变量，h_j(x)是节点j真值的概率分布；

(4)建立用于配电网状态估计极大似然估计目标函数如下：

$L\left(x\right)=\underset{x}{\max }\left(\ln \right(\underset{j=1}{\overset{m_1}{\Π}}\frac{1}{{\mathrm{nd}}_j}\underset{i=1}{\overset{n}{\Σ}}K\left(-\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j}\right))+\ln \underset{j=1}{\overset{m_2}{\Π}}\frac{1}{{\mathrm{\σ}}_j\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2})$

s.t.

c(x)＝0，为零注入母线的约束方程，

其中，x为状态变量矢量，h_ij(x)为步骤(1)中的量测真值对状态变量的函数，该函数的具体表达式由配电网络拓扑和潮流方程决定，是负荷伪量测值的第i个采样值；

(5)建立配电网状态估计模型，步骤如下：

(5-1)模型简化:

$\begin{array}{c}L\left(x\right)=\underset{x}{\max }\left(\ln \right(\underset{j=1}{\overset{m_1}{\Π}}\frac{1}{{\mathrm{nd}}_j}\underset{i=1}{\overset{n}{\Σ}}K\left(-\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j}\right))+\ln \underset{j=1}{\overset{m_2}{\Π}}\frac{1}{{\mathrm{\σ}}_j\sqrt{2\mathrm{\π}}}{e}^{-\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2})\\ =\underset{x}{\max }(\underset{j=1}{\overset{m_1}{\Σ}}\[\ln \frac{1}{{\mathrm{nd}}_j}+\ln \underset{i=1}{\overset{n}{\Σ}}K(-\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j})\]+\underset{j=1}{\overset{m_2}{\Σ}}-\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2-\frac{m_2}{2}\ln 2\mathrm{\π}-\underset{j=1}{\overset{m_2}{\Σ}}{\mathrm{ln\σ}}_j)\\ =\underset{x}{\max }\left(\underset{j=1}{\overset{m_1}{\Σ}}\ln \underset{i=1}{\overset{n}{\Σ}}K\right(-\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j})-\underset{j=1}{\overset{m_2}{\Σ}}\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2)\end{array}$

(5-2)根据琴生不等式:

$\begin{array}{c}L\left(x\right)=\underset{x}{\max }\left(\underset{j=1}{\overset{m_1}{\Σ}}\ln \underset{i=1}{\overset{n}{\Σ}}K\right(-\frac{h_{ij}-Z_{ij}^{pseudo}}{d_j})-\underset{j=1}{\overset{m_2}{\Σ}}\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2)\\ \≥\underset{x}{\max }\left(\underset{j=1}{\overset{m_1}{\Σ}}\underset{i=1}{\overset{n}{\Σ}}{W}_j\right(i\left)\ln \right(\frac{K(-\frac{h_{ij}-Z_{ij}^{pseudo}}{d_j})}{W_j\left(i\right)})-\underset{j=1}{\overset{m_2}{\Σ}}\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2)\end{array}$

而:

不等式取等号条件：

c为常数

则：

$\underset{x}{\min }\left(\underset{j=1}{\overset{m_1}{\Σ}}\underset{i=1}{\overset{n}{\Σ}}\frac{1}{2}{W}_j\right(i){\left(\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j}\right)}^2+\underset{j=1}{\overset{m_2}{\Σ}}\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2)$

(5-3)得到配电网状态估计模型如下：

$\underset{x}{\min }\left(\underset{j=1}{\overset{m_1}{\Σ}}\underset{i=1}{\overset{n}{\Σ}}\frac{1}{2}{W}_j\right(i){\left(\frac{h_{ij}\left(x\right)-Z_{ij}^{pseudo}}{d_j}\right)}^2+\underset{j=1}{\overset{m_2}{\Σ}}\frac{1}{2}{\left(\frac{h_j\left(x\right)-Z_j^{rt}}{{\mathrm{\σ}}_j}\right)}^2)$

s.t.c(x)＝0

其中，c(x)＝0是零注入母线的约束方程，其中包括配电网的有功功率零注入和无功功率零注入方程，x是状态变量矢量，包括支路首端功率P_lj，Q_lj和支路电流幅值I_lj；

其中：

$\begin{array}{c}{W}_j\left(i\right)=\frac{P\left(Z_{ij}^{pseudo}\right|i)}{\underset{i=1}{\overset{n}{\Σ}}P\left(Z_{ij}^{pseudo}\right|i)}=\frac{\frac{1}{{\mathrm{nd}}_j}K\left(\frac{Z_j^{true}-Z_{ij}^{pseudo}}{d_j}\right)}{\frac{1}{{\mathrm{nd}}_j}\underset{i=1}{\overset{n}{\Σ}}K\left(\frac{Z_j^{true}-Z_{ij}^{pseudo}}{d_j}\right)}\\ =\frac{K\left(\frac{Z_j^{true}-Z_{ij}^{pseudo}}{d_j}\right)}{\underset{i=1}{\overset{n}{\Σ}}K\left(\frac{Z_j^{true}-Z_{ij}^{pseudo}}{d_j}\right)}\end{array}$

支路电流幅值量测方程：

${\left(I_{lj}^m\right)}^2=I_{lj}^2+v_{I_{lj}^2}$

支路首端功率量测方程：

$P_{lj}^m=P_{lj}+v_{P_{lj}}$

$Q_{lj}^m=Q_{lj}+v_{Q_{lj}}$

支路末端功率量测方程:

$P_{jl}^m=P_{lj}+I_{lj}^2{R}_{lj}+v_{P_{lj}}$

$Q_{jl}^m=Q_{lj}+I_{lj}^2{X}_{lj}+v_{Q_{lj}}$

节点注入功率量测方程:

$P_j^m=\underset{l\∈j}{\Σ}{P}_{lj}-\underset{l\∈j}{\Σ}{I}_{lj}^2{R}_{lj}-\underset{k\∈j}{\Σ}{P}_{jk}+v_{P_j}$

$Q_j^m=\underset{l\∈j}{\Σ}{Q}_{lj}-\underset{l\∈j}{\Σ}{I}_{lj}^2{R}_{lj}-\underset{k\∈j}{\Σ}{Q}_{jk}-\frac{P_{jk}^2+Q_{jk}^2}{I_{jk}^2{X}_{cj}}+v_{Q_j}$

节点电压平方的量测方程：

${\left(U_j^m\right)}^2=\frac{{(P_{lj}-I_{lj}^2{R}_{lj})}^2+{(Q_{lj}-I_{lj}^2{X}_{lj})}^2}{I_{lj}^2}+v_{{\left(U_j^m\right)}^2}$

节点j的电压约束方程：

$\frac{{(P_{lj}-I_{lj}^2{R}_{lj})}^2+{(Q_{lj}-I_{lj}^2{X}_{lj})}^2}{I_{lj}^2}-\frac{P_{jk}^2+Q_{jk}^2}{I_{jk}^2}=0$

其中，l,k∈j表示配电网中与节点j相连的节点，X_cj为节点j的充电电容的电抗值，U_j为节点电压幅值，I_lj和I_jk分别为支路lj和支路jk电流幅值，P_lj、P_jk、Q_lj和Q_jk分别为支路lj与支路jk的首端有功功率和首端无功功率，P_j，Q_j分别表示节点注入有功和无功功率，是节点j的电压量测误差，是支路lj的电流量测的误差，和分别是支路lj的有功功率量测误差和有功功率量测误差，和分别是节点j的有功功率量测误差和无功功率量测误差，R_lj，X_lj分别为支路电阻和电抗，是出现的概率。