[發(fā)明專(zhuān)利]一種電力系統(tǒng)快速抗差狀態(tài)估計(jì)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410333797.7
• 申請(qǐng)日: 2014-07-14
• 公開(kāi)（公告）號(hào): CN104102836B
• 公開(kāi)（公告）日: 2018-07-17
• 發(fā)明（設(shè)計(jì)）人: 賈育培;郭子明;趙昆;王東升;陳利杰;張昊;鄒昱;宋磊;羅雅迪;閻博;張印;李靜;竇成龍;徐杰
• 申請(qǐng)（專(zhuān)利權(quán)）人: 國(guó)家電網(wǎng)公司;中國(guó)電力科學(xué)研究院;國(guó)網(wǎng)冀北電力有限公司
• 主分類(lèi)號(hào): G06F19/00
• 分類(lèi)號(hào): G06F19/00
• 代理公司: 北京安博達(dá)知識(shí)產(chǎn)權(quán)代理有限公司 11271
• 代理人: 徐國(guó)文

1.一種電力系統(tǒng)快速抗差狀態(tài)估計(jì)方法，其特征在于，所述方法包括下述步驟：

步驟1：依據(jù)電網(wǎng)模型和輸電斷面的數(shù)據(jù)獲取電力系統(tǒng)量測(cè)量方程Z的雅可比矩陣H；

步驟2：對(duì)所述雅可比矩陣H進(jìn)行正交變換得到有功迭代上三角矩陣L_P和無(wú)功迭代上三角矩陣L_Q；

步驟3：對(duì)所述矩陣L_P和L_Q分別進(jìn)行有功迭代計(jì)算和無(wú)功迭代計(jì)算獲取狀態(tài)量估計(jì)值和量測(cè)量殘差v；

步驟4：依據(jù)所述量測(cè)量殘差v對(duì)量測(cè)量權(quán)重進(jìn)行判斷，若需要進(jìn)行權(quán)重調(diào)整，則采用因子表增量修正法修正所述矩陣L_P和L_Q并返回步驟3；

所述步驟4中采用IGG抗差估計(jì)法依據(jù)所述量測(cè)量殘差v分別獲取處于保全區(qū)、降權(quán)區(qū)和淘汰區(qū)的量測(cè)量集合；

所述保全區(qū)的量測(cè)量值為正常量測(cè)值，不需進(jìn)行權(quán)重調(diào)整；所述降權(quán)區(qū)的量測(cè)量值為可

疑量測(cè)值，依據(jù)所述量測(cè)量值可信度進(jìn)行降權(quán)處理，并計(jì)算第i個(gè)量測(cè)量值的權(quán)重調(diào)整系數(shù)a_i；

所述淘汰區(qū)的量測(cè)量值為不良數(shù)據(jù)；

對(duì)不良數(shù)據(jù)進(jìn)行權(quán)重調(diào)整，將不良數(shù)據(jù)的權(quán)重置零后，對(duì)所述矩陣L_P和L_Q重新進(jìn)行迭代計(jì)算，得到新的狀態(tài)量估計(jì)值和量測(cè)量殘差；

所述步驟4中對(duì)可疑量測(cè)值進(jìn)行權(quán)重調(diào)整包括：

在所述矩陣L_P和L_Q的最后一行后增加j行修正元素；所述矩陣L_P和L_Q中新增加的對(duì)角元素值為(a_i-1)·r_i^-1；其中，1≤i≤m,1≤j≤m,r_i^-1為第i行的對(duì)角元素值，j為可疑量測(cè)值的個(gè)數(shù)，m為量測(cè)量個(gè)數(shù)；

分別對(duì)所述矩陣L_P和L_Q中增加的修正元素進(jìn)行有功迭代計(jì)算和無(wú)功迭代計(jì)算，得到新的狀態(tài)量估計(jì)值和量測(cè)量殘差。

2.如權(quán)利要求1所述的一種電力系統(tǒng)快速抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟2中采用吉文斯變換法對(duì)所述雅可比矩陣H進(jìn)行正交變換。

3.如權(quán)利要求1所述的一種電力系統(tǒng)快速抗差狀態(tài)估計(jì)方法，其特征在于，所述IGG抗差估計(jì)法的量測(cè)量殘差閾值包括k₁和k₂；其中，0＜k₁＜k₂,σ為單位權(quán)方差；

當(dāng)量測(cè)量殘差|ν|≤k₁σ時(shí)，所述量測(cè)量值為正常量測(cè)值，權(quán)重調(diào)整系數(shù)a_i＝1；

當(dāng)k₁σ≤|ν|≤k₂σ時(shí)，所述量測(cè)量值為可疑量測(cè)值，權(quán)重調(diào)整系數(shù)a_i＝k₁σ/|ν|；

當(dāng)|ν|≥k₂σ時(shí)，所述量測(cè)值為不良數(shù)據(jù)，權(quán)重調(diào)整系數(shù)a_i＝0。