基于wams量测的多区域电力系统状态估计方法及系统的制作方法
基于wams量测的多区域电力系统状态估计方法及系统的制作方法
【技术领域】
[0001] 本发明设及电力系统分析控制领域，尤其是一种基于WAMS量测的多区域电力系统 状态估计方法及系统。
【背景技术】
[0002] 电力系统是覆盖迂阔地域的大互联系统，精确估计其运行状态，为其它高级应用 软件提供实时可靠的系统状态信息，一直是电力工业界长久不衰的研究课题。随着基于相 量量测单元（Phasor Measurement 加 it,PMU)的广域量测系统（Wide-area Monitor System, WAMS)在技术上的成熟与推广应用，利用WAMS量测技术已可实现电力系统的线性状 态估计。但电力系统状态量的维数很高，量测的数据量过于庞大，会导致集中处理变得非常 烦杂。
[0003] 为克服大规模系统集中式状态估计难的问题，目前大多数研究会采用基于星型通 信拓扑的分层式状态估计算法，即需要一个协调中屯、与区域电网间进行通，信W协调整合 各区域局部估计值。随着电力系统自动化水平的提高，上述星型通信拓扑也会随之被取代， 逐渐发展成一种更一般化的通信拓扑结构一一网状或链状结构。目前，基于网状或链状通 信拓扑的分布式状态估计方法已有部分研究，但大部分算法存在收敛速度慢、算法参数不 易确定或对通信拓扑有特定要求等缺点，且都没有考虑到WAMS存在量测数据缺失的情况， 容易导致状态估计值存在较大的偏差，不够准确。
【发明内容】
[0004] 为解决上述技术问题，本发明的目的在于:提供一种收敛速度快、适用性广且准确 的，基于WAMS量测的多区域电力系统状态估计方法。
[000引本发明的另一目的在于:提供一种收敛速度快、适用性广且准确的，基于WAMS量测 的多区域电力系统状态估计系统。
[0006 ]本发明解决其技术问题所采取的技术方案是：
[0007] 基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括W下步骤：
[0008] S1、建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；
[0009] S2、各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得 到本地状态估计值；
[0010] S3、各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行 修正，所述分布式一致性协议W各子区域的电网调度中屯、作为智能体，各个智能体之间通 过链状或网状通信拓扑联通；
[0011] S4、判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情 况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计 值，反之，则直接将各子区域修正后的本地状态估计值进行输出。
[0012] 进一步，所述步骤S1，其包括：
[0013] 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0014]
[0015] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0016] 依次通过第一方程组、第二方程组和第Ξ方程组求出多区域状态估计模型的最优 解，所述第一方程组的表达式为：
，所述第二方程组的表达式为：
其中，式为区域i未经 修正的本地状态估计值，λ为拉格朗日乘子，為为区域i的最优状态估计值，di为区域i的估 计协调增量。
[0017] 进一步，所述步骤S2包括：
[0018] S21、各子区域获取当前时刻本地上传的WAMS量测；
[0019] S22、各子区域根据获取的WAMS量测W及多区域状态估计模型采用分布式最小二 乘估计算法估计本地状态，得到本地状态估计值。
[0020] 进一步，所述步骤S22包括：
[0021] S221、进行初始化，为各区域i置一致性协议初始值天
，然后各区域i 执行有限时间平均一致性协议得到^馬；
[0022] S222、各区域i根据和第一方程组独立计算禾和公，式，并将马克，的值赋予曰1,〇;
[0023] S223、各区域i根据赋值后的新曰1,日执行有限时间平均一致性协议，得質
[0024] S224、各区域i根巧
和第二方程组独立计算拉格朗日乘子λ，然后由λ和第 Ξ方程组计算本地状态估计值衣 〇
[00巧]进一步，所述步骤S3包括：
[0026] 各子区域执行分布式一致性协议，得到各子区域的估计协调增量；
[0027] 各子区域根据得到的估计协调增量对本地状态估计值进行修正。
[0028] 进一步，所述步骤S4包括：
[0029] S41、各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的值来判断本 地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定本地上传的 WAMS量测有量测缺失，此时执行步骤S42;若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测无量测缺失，此时直接跳至步骤S43;
[0030] S42、迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计 值收敛至剔除缺失量测后的集中式估计值，然后执行步骤S43;
[0031] S43、将各子区域修正后的本地状态估计值进行输出。
[0032] 进一步，所述步骤S42包括：
[0033] S421、给定收敛阔值eend，并为各区域i置状态估计值初值= -?·, W及时刻初值t =0；
[0034] S422、各区域i根据时刻t的状态估计值\,独立计算
和 巧'i,,，并将巧韦,,的值赋予曰1,〇;
[003引S423、各区域i根据赋值后的新ai,日执行有限时间平均一致性协议，得到
[0036] S424、各区域i根据
和第二方程组独立计算时刻t的拉格朗日乘子λ*，并 计算各区域i在时刻t+1的状态估计值，计算公式为：
[0037]
[003引S425、各区域i计算巧二max(|i,.灿，然后置ei,0 = ei，执行最大一致性协议， 所述最大一致性协议的表达式为：
[0039]
[0040] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0041] S426、根据最大一致性协议的执行结果判断吗,<馬,是否成立，若是，则结束迭 代，将时刻t的状态估计值馬^作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤S422。
[0042 ]本发明解决其技术问题所采取的另一技术方案是：
[0043]基于WAMS量测的多区域电力系统状态估计系统，包括W下模块：
[0044] 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状 态估计模型的最优解求解公式；
[0045] 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模 型估计本地状态，得到本地状态估计值；
[0046] 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行 的结果对本地状态估计值进行修正，所述分布式一致性协议W各子区域的电网调度中屯、作 为智能体，各个智能体之间通过链状或网状通信拓扑联通；
[0047] 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS 量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的 集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出；
[0048] 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布 式一致性协议执行模块进而与判断模块的输入端连接。
[0049] 进一步，所述多区域状态估计模型构建模块包括：
[0050] 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0051]
[0052] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0053] 最优解求解单元，用于依次通过第一方程组、第二方程组和第Ξ方程组求出多区 域状态估计模型的最优解，所述第一方程组的表达式为：
，所述第二方程组的 表达式为：
所述第Ξ方程组的表达式为：
其中，睾为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，兩为区域i的最优状态估计值，di为 区域i的估计协调增量；
[0054] 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端 连接。
[00巧]进一步，所述判断模块包括：
[0056] 判 断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的 值来判断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定 本地上传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵Si对角元素的值 为1，则判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元；
[0057] 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后 的本地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元；
[0058] 输出单元，用于将各子区域修正后的本地状态估计值进行输出；
[0059] 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单 元的输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的 输出端还与输出单元的输入端连接。
[0060] 本发明的方法的有益效果是:各子区域先根据本地上传的WAMS量测估计其辖区内 系统状态，然后执行一次分布式一致性协议来求解修正本地状态估计值，并判断地上传的 WAMS量测是否有量测缺失，若当前时刻WAMS有部分量测缺失，则通过WAMS量测缺失情况下 的迭代修正算法来进行修正；引入了分布式一致性协议，在WAMS量测无缺失的情况下，能快 速收敛至集中式估计值，收敛速度快，且算法参数易确定，对通信拓扑无任何特定要求，适 用性广;在出现WAMS量测缺失时，无需改变原有的信息矩阵结构，只需执行若干次迭代修正 算法，就能使各子区域本地估计值收敛至剔除缺失量测后的集中式估计值，更加准确。
[0061] 本发明的系统的有益效果是:各子区域先根据本地上传的WAMS量测估计其辖区内 系统状态，然后执行一次分布式一致性协议来求解修正本地状态估计值，并判断地上传的 WAMS量测是否有量测缺失，若当前时刻WAMS有部分量测缺失，则通过WAMS量测缺失情况下 的迭代修正算法来进行修正；引入了分布式一致性协议，在WAMS量测无缺失的情况下，能快 速收敛至集中式估计值，收敛速度快，且算法参数易确定，对通信拓扑无任何特定要求，适 用性广;在出现WAMS量测缺失时，无需改变原有的信息矩阵结构，只需执行若干次迭代修正 算法，就能使各子区域本地估计值收敛至剔除缺失量测后的集中式估计值，更加准确。
【附图说明】
[0062] 图1为本发明基于WAMS量测的多区域电力系统状态估计方法的整体流程图；
[0063] 图2为本发明基于WAMS量测的多区域电力系统状态估计系统的整体框图；
[0064] 图3为二分区互联电力系统的结构示意图；
[0065] 图4为实施例一互联电力系统的状态估计方法流程图；
[0066] 图5为IE邸118节点标准测试系统的7分区示意图；
[0067] 图6为本发明图5系统对应的通信拓扑。
【具体实施方式】
[0068] 参照图1，基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括W下 步骤：
[0069] S1、建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；
[0070] S2、各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得 到本地状态估计值；
[0071] S3、各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行 修正，所述分布式一致性协议W各子区域的电网调度中屯、作为智能体，各个智能体之间通 过链状或网状通信拓扑联通；
[007引S4、判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情 况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计 值，反之，则直接将各子区域修正后的本地状态估计值进行输出。
[0073] 进一步作为优选的实施方式，所述步骤S1，其包括：
[0074] 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0075]
[0076] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[oow]依次通过第一方程组、第二方程组和第Ξ方程组求出多区域状态估计模型的最优 解，所述第一方程组的表达式为：
，所述第二方程组的表达式为：
所述第Ξ方程组的表达式为： ，其中，^,为区域i未经 ， ， 修正的本地状态估计值，λ为拉格朗日乘子，X为区域i的最优状态估计值，di为区域i的估 计协调增量。
[0078] 进一步作为优选的实施方式，所述步骤S2包括：
[0079] S21、各子区域获取当前时亥陆地上传的WAMS量测；
[0080] S22、各子区域根据获取的WAMS量测W及多区域状态估计模型采用分布式最小二 乘估计算法估计本地状态，得到本地状态估计值。
[0081 ]进一步作为优选的实施方式，所述步骤S22包括：
[0082] S221、进行初始化，为各区域i置一致性协议初始值为巧,0=马巧1《，然后各区域i 执行有限时间平均一致性协议得到；
[0083] S222、各区域i根据和第一方程组独立计算兩和巧兩，并将公式的值赋予曰1,〇;
[0084] S223、各区域i根据赋值后的新曰1,〇执行有限时间平均一致性协议，得到
[008引 S224、各区域i根巧
和第二方程组独立计算拉格朗日乘子λ，然后由λ和第Ξ 方程组计算本地状态估计值為： .0.
[0086] 进一步作为优选的实施方式，所述步骤S3包括：
[0087] 各子区域执行分布式一致性协议，得到各子区域的估计协调增量；
[0088] 各子区域根据得到的估计协调增量对本地状态估计值进行修正。
[0089] 进一步作为优选的实施方式，所述步骤S4包括：
[0090] S41、各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的值来判断本 地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定本地上传的 WAMS量测有量测缺失，此时执行步骤S42;若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测无量测缺失，此时直接跳至步骤S43;
[0091] S42、迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计 值收敛至剔除缺失量测后的集中式估计值，然后执行步骤S43;
[0092] S43、将各子区域修正后的本地状态估计值进行输出。
[0093] 进一步作为优选的实施方式，所述步骤S42包括：
[0094] S421、给定收敛阔值eend，并为各区域i置状态估计值初值=韦.W及时刻初值t =0；
[009引 S422、各区域i根据时刻t的状态估计值卷,独立计算
和 ，并将巧韦,，的值赋予ai,0;
[0096] S423、各区域i根据赋值后的新ai,日执行有限时间平均一致性协议，得到
[0097] S424、各区域i根据
巧第二方程组独立计算时刻t的拉格朗日乘子λ*，并 计算各区域i在时刻t+1的状态估计值·%κ;ι，计算公式为：
[009引
[0099] S425、各区域i计算
然后置ei,〇 = ei，执行最大一致性协议， 所述最大一致性协议的表达式为：
[0100]
[0101] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0102] S426、根据最大一致性协议的执行结果判断< Ceud是否成立，若是，则结束迭 代，将时刻t的状态估计值作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤S422。
[0103] 参照图2,基于WAMS量测的多区域电力系统状态估计系统，包括W下模块：
[0104] 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状 态估计模型的最优解求解公式；
[0105] 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模 型估计本地状态，得到本地状态估计值；
[0106] 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行 的结果对本地状态估计值进行修正，所述分布式一致性协议W各子区域的电网调度中屯、作 为智能体，各个智能体之间通过链状或网状通信拓扑联通；
[0107] 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS 量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的 集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出；
[0108] 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布 式一致性协议执行模块进而与判断模块的输入端连接。
[0109] 参照图2,进一步作为优选的实施方式，所述多区域状态估计模型 构建模块包括：
[0110] 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0111]
[0112] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，xi为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0113] 最优解求解单元，用于依次通过第一方程组、第二方程组和第Ξ方程组求出多区 域状态估计模型的最优解，所述第一方程组的表达式为
所述第二方程组的 表达式为：
，其中，疋为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，为区域i的最优状态估计值，dl为 区域i的估计协调增量；
[0114] 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端 连接。
[0115] 参照图2,进一步作为优选的实施方式，所述判断模块包括：
[0116] 判断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的 值来判断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定 本地上传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵Si对角元素的值 为1，则判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元；
[0117] 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后 的本地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元；
[0118] 输出单元，用于将各子区域修正后的本地状态估计值进行输出；
[0119] 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单 元的输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的 输出端还与输出单元的输入端连接。
[0120] 下面结合说明书附图和具体实施例对本发明作进一步详细说明。
[0121] 实施例一
[0122] 本实施例对本发明所设及的相关理论及本发明的具体实现过程进行详细介绍。
[0123] 针对传统WAMS具有庞大数据量、存在量测数据缺失的问题和新型通信拓扑（链状 或网状)的发展，本发明提出了一种针对WAMS量测缺失的多区域电力系统状态估计方法，如 图4所示，无需协调中屯、，能适应任意通信拓扑;在WAMS正常情况(WAMS量测无缺失)下，能快 速收敛至集中式估计值;在WAMS出现部分PMU量测缺失时，无需改变原有信息矩阵结构，只 需执行若干次分布式一致性协议能使各子区域本地估计值收敛至剔除缺失量测后的集中 式估计值。
[0124] (一)多区域电力系统状态估计模型及相关定义
[0。5]如图3所示的二分区互联电力系统，区域i的状态量矢i
定义 为：
[0126] a)区域i内部状态XI,int:区域i内不与其它区域相连接的节点状态量；
[0127] b)区域i内部边界状态XI,bint:区域i内与其它区域相连接的节点状态量；
[0128] C)区域i外部边界状态XI,bext:其它区域内与区域i相连接的节点状态量。
[0129 ]类似地，可定义出区域j的状态量矢量。
[0130] 本发明建立的多区域电力系统状态估计模型为：
[0131]
[013引其中，Z功区域i的量测;Wi为区域i的权重矩阵;B为4bXn阶矩阵，B每行中有两个 非零元素为1或-1，分别对应两个相邻区域的边界状态量；b为联络线支路数。令B = [Bi，...，Bn]，则有：
[0133]
[0134] 引入拉格朗日乘子λ，上述多区域电力系统状态估计模型的最优求解过程可通过 W下Ξ个方程组来实现：
[0138] 其中，為为区域i的最优状态估计值，di为区域i的估计协调增量。
[0139] 由W上Ξ个方程组可知，λ需要各区域联合求解。
[0140] 故最优求解过程的第二方程组可等效改写为：
[0141]
[0142] 各区域i仅已知Bi，即仅已知矩阵Β中关于区域i边界状态的部分列向量，而区域i 的初始状态设为：
[0145]分布式一致性协议可采用现有的有限时间平均一致性协议来实现。各区域执行分 布式一致性协议后可得到一致平均值
，然后再通过分布式加权最小二乘估 计算法即可得到本地状态估计值，因为为常数矩阵，初始化只需一次。
[0146 ](二)分布式加权最小二乘估计算法
[0147]本发明的分布式加权最小二乘估计算法的具体过程为：
[014引（1)进行初始化，为各区域i置一致性协议初始值为巧G; I巧'，然后各区域i执 行有限时间平均一致性协议得到^辟；
[0149] (2)各区域i根据!：&和(一）的第一方程组独立计算韦和度式，并将公,?,的值赋予 过i，〇;
[0150] (3)各区域i根据赋值后的新ai，〇执行有限时间平均一致性协议，得到
；
[0151] (4)各区域i根据
巧(一)的第二方程组独立计算拉格朗日乘子λ，然后由λ和 (一)的第Ξ方程组计算本地状态估计值4。
[0152] (S)WAMS量测缺失情况下的迭代修正算法
[0153] 定理设缺失量测剔除后的系统能保持可观性，给定初值：
[0154]
[015引由
爱代可收敛至：
[0156]
[0157] 其中，S为对角线矩阵，对应异常量测（即有量测缺失）的对角元素为0,对应正常量 巧。（即无量测缺失）的对角元素为l;Zm、Am和Rm分别为不包含异常量测的量测向量、量测方程 矩阵和量测协方差矩阵。
[0158] 根据定理，在(二）的算法步骤中引入对角线矩阵Si，则（二）中最优求解过程的第 一方程组第一式子修正为：
[0159]
[0160] 并置。,。=昕巧，巧Δ?,·.,)叩。Tapi为区域i的WAMS量测缺失标志，当区域i有WAMS量测 缺失时，Tapi = l;否则，Tapi = 0。执行完(二）的算法后，各区域i检查Tapi是否为0,如不为0， 则继续执行W下的迭代修正算法：
[0161] (1)给定收敛阔值eend，并为各区域i置状态估计值初值韦,0 =韦W及时刻初值t = 0；
[0162] (2)各区域i根据时刻t的状态估计值為独立计算 和马％， 并将巧..iy的值赋予ai,0;
[0163] (3)各区域i根据赋值后的新ai,〇执行有限时间平均一致性协议，得到
[0164] (4)各区域i根据
和(一）的第二方程组独立计算时刻t的拉格朗日乘子 λ*，并计算各区域i在时刻t+1的状态估计值，计算公式为：
[0165]
[0166] (5)各区域i计算
，然后置ei,〇 = ei，执行最大一致性协议，所 述最大一致性协议的表达式为：
[0167]
[0168] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0169] (6)根据最大一致性协议的执行结果判断6,為是否成立，若是，则结束迭 代，将时刻t的状态估计值焉>,作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤(2)。
[0170] 与现有技术相比，本实施例具有W下优点：
[0171] (1)根据互联电网已有的分区方式，将子区域电网调度中屯、与分布式一致性协议 中的智能体概念相联系，并通过链状或网状通信拓扑联通，W实现无协调中屯、的多区域状 态估计，效率更高。
[0172] (2)在估计时刻，各子区域间通过与相邻区域交换边界信息，得到估计结果得到全 局协调，保证各子区域本地估计值与集中式估计相同，更加科学。
[0173] (3)建立各子区域电网间状态估计模型的禪合关系，然后基于该关系合理引入分 布式一致性协议来实现全局协调求解，更加全面。
[0174] (4)扩展推导出WAMS量测存在部分缺失情况下的修正算法，对WAMS量测缺失具有 较强的适应性，更加全面和准确。
[017引实施例二
[0176] 1)W图5所示7分区的IEEE 118节点标准测试系统为例，所有节点均配置PMU，相应 节点电压相量和与该节点相联的支路电流相量是可量测的。节点电压实部与虚部量测值为 潮流真值叠加标准差σ = 0.008的高斯白噪声，支路电流实部与虚部量测值为潮流真值叠加 标准差〇 = 0.0005的高斯白噪声，权重矩阵W取为。图6为对应图5系统的通信拓扑。 巧
[01W]根据先验知识与实际测试，经过有限时间平均一致性迭代，各区域节点电压相量 的实虚部均收敛至集中式估计值。
[0178] 假设各区域均出现一个电压相量实部的量测缺失，即全系统共有7个量测缺失，并 将收敛阔值eend设为0.0001。若采用本发明的迭代修正算法进行测试，则收敛时所需的迭代 步数心=3。表1为其它缺失数据情况下采用本发明的迭代修正算法的收敛 结果。从表1我们 可W看到，t。随着系统缺失量测总数的增多而增大，但增幅并不明显。一般地，WAMS可靠性 较高，一次采样中的缺失量测数相对较少，故采用本发明的算法能通过较少的迭代步数迭 代至收敛。
[0179] 表1WAMS量测缺失情况下采用迭代修正算法的收敛结果
[0180]
~w上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施胃 例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替 换，运些等同的变形或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括以下步骤： 51、 建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式； 52、 各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得到本 地状态估计值； 53、 各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行修正， 所述分布式一致性协议以各子区域的电网调度中心作为智能体，各个智能体之间通过链状 或网状通信拓扑联通； 54、 判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情况下 的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计值， 反之，则直接将各子区域修正后的本地状态估计值进行输出。2. 根据权利要求1所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤Sl，其包括： 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：其中，J为多区域状态估计模型的目标函数，N为多区域的总数，Z1为区域i的WAMS量测， Ai为区域i的量测方程矩阵，Xi为区域i的状态矢量，Wi为区域i的权重矩阵;x=[X1，. . .，XN] 为η维全系统状态矢量，B为4b X η阶矩阵，b为联络线支路数，B = [Bi，. . .，Bn]，满足依次通过第一方程组、第二方程组和第三方程组求出多区域状态估计模型的最优解， 所述第一方程组的表达式为：所述第二方程组的表达式为：所述第三方程组的表达式为：其中，g为区域i 未经修正的本地状态估计值，λ为拉格朗日乘子，为区域i的最优状态估计值，Cl1为区域i 的估计协调增量。3. 根据权利要求2所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S2包括： 521、 各子区域获取当前时刻本地上传的WAMS量测； 522、 各子区域根据获取的WAMS量测以及多区域状态估计模型采用分布式最小二乘估 计算法估计本地状态，得到本地状态估计值。4. 根据权利要求3所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S22包括： 5221、 进行初始化，为各区域i置一致性协议初始值为然后各区域i执行 有限时间平均一致性协议得到; 5222、 各区域i根据告Gs和第一方程组独立计算笔和尽足，并将尽矣的值赋予ai,〇; 5223、 各区域i根据赋值后的新〇执行有限时间平均一致性协议，得到5224、 各区域i根据和第二方程组独立计算拉格朗日乘子λ，然后由λ和第三方 程组计算本地状态估计值· 〇5. 根据权利要求2所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S3包括： 各子区域执行分布式一致性协议，得到各子区域的估计协调增量； 各子区域根据得到的估计协调增量对本地状态估计值进行修正。6. 根据权利要求4所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S4包括： 541、 各子区域i根据本地上传的WAMS量测中对角线矩阵S1对角元素的值来判断本地上 传的WAMS量测是否有量测缺失，若对角线矩阵S i对角元素的值为1，则判定本地上传的WAMS 量测有量测缺失，此时执行步骤S42;若对角线矩阵Sdt角元素的值为1，则判定本地上传的 WAMS量测无量测缺失，此时直接跳至步骤S43; 542、 迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收 敛至剔除缺失量测后的集中式估计值，然后执行步骤S43; 543、 将各子区域修正后的本地状态估计值进行输出。7. 根据权利要求6所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S42包括： 5421、 给定收敛阈值eend，并为各区域i置状态估计值初值毛，。=天以及时刻初值t = 0; 5422、 各区域i根据时刻t的状态估计值&.^虫立计算并将尽足7的值赋予ai,o; 5423、 各区域i根据赋值后的新〇执行有限时间平均一致性协议，得到5424、 各区域i根据和第二方程组独立计算时刻t的拉格朗日乘子&，并计 算各区域i在时刻t+1的状态估计值7.+|计算公式为：5425、 各区域i计算_然后置ei,〇 = ei，执行最大一致性协 议，所述最大一致性协议的表达式为：其中，Q1表示与区域i相邻的区域编号集合，dc为简单无向图G(v，〇的直径，V为图G(v， O的顶点集合，ε为图G(v，〇的边集合； 5426、 根据最大一致性协议的执行结果判断是否成立，若是，则结束迭 代，将时刻t的状态估计值作为各子区域修正后的本地状态估计值进行输出，反之，则 令t自加1然后返回步骤S422。8. 基于WAMS量测的多区域电力系统状态估计系统，其特征在于:包括以下模块： 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状态估 计模型的最优解求解公式； 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模型估 计本地状态，得到本地状态估计值； 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行的结 果对本地状态估计值进行修正，所述分布式一致性协议以各子区域的电网调度中心作为智 能体，各个智能体之间通过链状或网状通信拓扑联通； 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测 缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中 式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出； 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布式一 致性协议执行模块进而与判断模块的输入端连接。9. 根据权利要求8所述的多区域电力系统状态估计系统，其特征在于:所述多区域状态 估计模型构建模块包括： 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：其中，J为多区域状态估计模型的目标函数，N为多区域的总数，Z1为区域i的WAMS量测， Ai为区域i的量测方程矩阵，Xi为区域i的状态矢量，Wi为区域i的权重矩阵;X=[X1，. . .，XN] 为η维全系统状态矢量，B为4b X η阶矩阵，b为联络线支路数，B = [Bi，. . .，Bn]，满足最优解求解单元，用于依次通过第一方程组、第二方程组和第三方程组求出多区域状 态估计模型的最优解，所述第一方程组的表达式为：，所述第二方程组的表达 式为：所述第三方程组的表达式为：其中，天为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，$为区域i的最优状态估计值，Cl1为 区域i的估计协调增量； 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端连接。10.根据权利要求8所述的多区域电力系统状态估计系统，其特征在于:所述判断模块 包括： 判断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵S1对角元素的值来判 断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵SM角元素的值为1，则 判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元； 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本 地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元； 输出单元，用于将各子区域修正后的本地状态估计值进行输出； 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单元的 输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的输出 端还与输出单元的输入端连接。
【专利摘要】本发明公开了一种基于WAMS量测的多区域电力系统状态估计方法及系统，方法包括：建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得到本地状态估计值；各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行修正；判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出。本发明具有收敛速度快、适用性广且准确的优点，可广泛应用于电力系统分析控制领域。
【IPC分类】G06F19/00, H02J3/00
【公开号】CN105490269
【申请号】CN201511028254
【发明人】陈亦平, 蔡永智, 侯君, 陈皓勇, 杨林, 郭曼兰, 王巍, 万楚林, 刘起兴
【申请人】中国南方电网有限责任公司, 华南理工大学
【公开日】2016年4月13日
【申请日】2015年12月30日
【技术领域】
[0001] 本发明设及电力系统分析控制领域，尤其是一种基于WAMS量测的多区域电力系统 状态估计方法及系统。
【背景技术】
[0002] 电力系统是覆盖迂阔地域的大互联系统，精确估计其运行状态，为其它高级应用 软件提供实时可靠的系统状态信息，一直是电力工业界长久不衰的研究课题。随着基于相 量量测单元（Phasor Measurement 加 it,PMU)的广域量测系统（Wide-area Monitor System, WAMS)在技术上的成熟与推广应用，利用WAMS量测技术已可实现电力系统的线性状 态估计。但电力系统状态量的维数很高，量测的数据量过于庞大，会导致集中处理变得非常 烦杂。
[0003] 为克服大规模系统集中式状态估计难的问题，目前大多数研究会采用基于星型通 信拓扑的分层式状态估计算法，即需要一个协调中屯、与区域电网间进行通，信W协调整合 各区域局部估计值。随着电力系统自动化水平的提高，上述星型通信拓扑也会随之被取代， 逐渐发展成一种更一般化的通信拓扑结构一一网状或链状结构。目前，基于网状或链状通 信拓扑的分布式状态估计方法已有部分研究，但大部分算法存在收敛速度慢、算法参数不 易确定或对通信拓扑有特定要求等缺点，且都没有考虑到WAMS存在量测数据缺失的情况， 容易导致状态估计值存在较大的偏差，不够准确。
【发明内容】
[0004] 为解决上述技术问题，本发明的目的在于:提供一种收敛速度快、适用性广且准确 的，基于WAMS量测的多区域电力系统状态估计方法。
[000引本发明的另一目的在于:提供一种收敛速度快、适用性广且准确的，基于WAMS量测 的多区域电力系统状态估计系统。
[0006 ]本发明解决其技术问题所采取的技术方案是：
[0007] 基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括W下步骤：
[0008] S1、建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；
[0009] S2、各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得 到本地状态估计值；
[0010] S3、各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行 修正，所述分布式一致性协议W各子区域的电网调度中屯、作为智能体，各个智能体之间通 过链状或网状通信拓扑联通；
[0011] S4、判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情 况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计 值，反之，则直接将各子区域修正后的本地状态估计值进行输出。
[0012] 进一步，所述步骤S1，其包括：
[0013] 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0014]
[0015] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0016] 依次通过第一方程组、第二方程组和第Ξ方程组求出多区域状态估计模型的最优 解，所述第一方程组的表达式为：
，所述第二方程组的表达式为：
其中，式为区域i未经 修正的本地状态估计值，λ为拉格朗日乘子，為为区域i的最优状态估计值，di为区域i的估 计协调增量。
[0017] 进一步，所述步骤S2包括：
[0018] S21、各子区域获取当前时刻本地上传的WAMS量测；
[0019] S22、各子区域根据获取的WAMS量测W及多区域状态估计模型采用分布式最小二 乘估计算法估计本地状态，得到本地状态估计值。
[0020] 进一步，所述步骤S22包括：
[0021] S221、进行初始化，为各区域i置一致性协议初始值天
，然后各区域i 执行有限时间平均一致性协议得到^馬；
[0022] S222、各区域i根据和第一方程组独立计算禾和公，式，并将马克，的值赋予曰1,〇;
[0023] S223、各区域i根据赋值后的新曰1,日执行有限时间平均一致性协议，得質
[0024] S224、各区域i根巧
和第二方程组独立计算拉格朗日乘子λ，然后由λ和第 Ξ方程组计算本地状态估计值衣 〇
[00巧]进一步，所述步骤S3包括：
[0026] 各子区域执行分布式一致性协议，得到各子区域的估计协调增量；
[0027] 各子区域根据得到的估计协调增量对本地状态估计值进行修正。
[0028] 进一步，所述步骤S4包括：
[0029] S41、各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的值来判断本 地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定本地上传的 WAMS量测有量测缺失，此时执行步骤S42;若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测无量测缺失，此时直接跳至步骤S43;
[0030] S42、迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计 值收敛至剔除缺失量测后的集中式估计值，然后执行步骤S43;
[0031] S43、将各子区域修正后的本地状态估计值进行输出。
[0032] 进一步，所述步骤S42包括：
[0033] S421、给定收敛阔值eend，并为各区域i置状态估计值初值= -?·, W及时刻初值t =0；
[0034] S422、各区域i根据时刻t的状态估计值\,独立计算
和 巧'i,,，并将巧韦,,的值赋予曰1,〇;
[003引S423、各区域i根据赋值后的新ai,日执行有限时间平均一致性协议，得到
[0036] S424、各区域i根据
和第二方程组独立计算时刻t的拉格朗日乘子λ*，并 计算各区域i在时刻t+1的状态估计值，计算公式为：
[0037]
[003引S425、各区域i计算巧二max(|i,.灿，然后置ei,0 = ei，执行最大一致性协议， 所述最大一致性协议的表达式为：
[0039]
[0040] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0041] S426、根据最大一致性协议的执行结果判断吗,<馬,是否成立，若是，则结束迭 代，将时刻t的状态估计值馬^作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤S422。
[0042 ]本发明解决其技术问题所采取的另一技术方案是：
[0043]基于WAMS量测的多区域电力系统状态估计系统，包括W下模块：
[0044] 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状 态估计模型的最优解求解公式；
[0045] 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模 型估计本地状态，得到本地状态估计值；
[0046] 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行 的结果对本地状态估计值进行修正，所述分布式一致性协议W各子区域的电网调度中屯、作 为智能体，各个智能体之间通过链状或网状通信拓扑联通；
[0047] 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS 量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的 集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出；
[0048] 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布 式一致性协议执行模块进而与判断模块的输入端连接。
[0049] 进一步，所述多区域状态估计模型构建模块包括：
[0050] 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0051]
[0052] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0053] 最优解求解单元，用于依次通过第一方程组、第二方程组和第Ξ方程组求出多区 域状态估计模型的最优解，所述第一方程组的表达式为：
，所述第二方程组的 表达式为：
所述第Ξ方程组的表达式为：
其中，睾为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，兩为区域i的最优状态估计值，di为 区域i的估计协调增量；
[0054] 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端 连接。
[00巧]进一步，所述判断模块包括：
[0056] 判 断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的 值来判断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定 本地上传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵Si对角元素的值 为1，则判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元；
[0057] 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后 的本地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元；
[0058] 输出单元，用于将各子区域修正后的本地状态估计值进行输出；
[0059] 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单 元的输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的 输出端还与输出单元的输入端连接。
[0060] 本发明的方法的有益效果是:各子区域先根据本地上传的WAMS量测估计其辖区内 系统状态，然后执行一次分布式一致性协议来求解修正本地状态估计值，并判断地上传的 WAMS量测是否有量测缺失，若当前时刻WAMS有部分量测缺失，则通过WAMS量测缺失情况下 的迭代修正算法来进行修正；引入了分布式一致性协议，在WAMS量测无缺失的情况下，能快 速收敛至集中式估计值，收敛速度快，且算法参数易确定，对通信拓扑无任何特定要求，适 用性广;在出现WAMS量测缺失时，无需改变原有的信息矩阵结构，只需执行若干次迭代修正 算法，就能使各子区域本地估计值收敛至剔除缺失量测后的集中式估计值，更加准确。
[0061] 本发明的系统的有益效果是:各子区域先根据本地上传的WAMS量测估计其辖区内 系统状态，然后执行一次分布式一致性协议来求解修正本地状态估计值，并判断地上传的 WAMS量测是否有量测缺失，若当前时刻WAMS有部分量测缺失，则通过WAMS量测缺失情况下 的迭代修正算法来进行修正；引入了分布式一致性协议，在WAMS量测无缺失的情况下，能快 速收敛至集中式估计值，收敛速度快，且算法参数易确定，对通信拓扑无任何特定要求，适 用性广;在出现WAMS量测缺失时，无需改变原有的信息矩阵结构，只需执行若干次迭代修正 算法，就能使各子区域本地估计值收敛至剔除缺失量测后的集中式估计值，更加准确。
【附图说明】
[0062] 图1为本发明基于WAMS量测的多区域电力系统状态估计方法的整体流程图；
[0063] 图2为本发明基于WAMS量测的多区域电力系统状态估计系统的整体框图；
[0064] 图3为二分区互联电力系统的结构示意图；
[0065] 图4为实施例一互联电力系统的状态估计方法流程图；
[0066] 图5为IE邸118节点标准测试系统的7分区示意图；
[0067] 图6为本发明图5系统对应的通信拓扑。
【具体实施方式】
[0068] 参照图1，基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括W下 步骤：
[0069] S1、建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；
[0070] S2、各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得 到本地状态估计值；
[0071] S3、各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行 修正，所述分布式一致性协议W各子区域的电网调度中屯、作为智能体，各个智能体之间通 过链状或网状通信拓扑联通；
[007引S4、判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情 况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计 值，反之，则直接将各子区域修正后的本地状态估计值进行输出。
[0073] 进一步作为优选的实施方式，所述步骤S1，其包括：
[0074] 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0075]
[0076] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，XI为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[oow]依次通过第一方程组、第二方程组和第Ξ方程组求出多区域状态估计模型的最优 解，所述第一方程组的表达式为：
，所述第二方程组的表达式为：
所述第Ξ方程组的表达式为： ，其中，^,为区域i未经 ， ， 修正的本地状态估计值，λ为拉格朗日乘子，X为区域i的最优状态估计值，di为区域i的估 计协调增量。
[0078] 进一步作为优选的实施方式，所述步骤S2包括：
[0079] S21、各子区域获取当前时亥陆地上传的WAMS量测；
[0080] S22、各子区域根据获取的WAMS量测W及多区域状态估计模型采用分布式最小二 乘估计算法估计本地状态，得到本地状态估计值。
[0081 ]进一步作为优选的实施方式，所述步骤S22包括：
[0082] S221、进行初始化，为各区域i置一致性协议初始值为巧,0=马巧1《，然后各区域i 执行有限时间平均一致性协议得到；
[0083] S222、各区域i根据和第一方程组独立计算兩和巧兩，并将公式的值赋予曰1,〇;
[0084] S223、各区域i根据赋值后的新曰1,〇执行有限时间平均一致性协议，得到
[008引 S224、各区域i根巧
和第二方程组独立计算拉格朗日乘子λ，然后由λ和第Ξ 方程组计算本地状态估计值為： .0.
[0086] 进一步作为优选的实施方式，所述步骤S3包括：
[0087] 各子区域执行分布式一致性协议，得到各子区域的估计协调增量；
[0088] 各子区域根据得到的估计协调增量对本地状态估计值进行修正。
[0089] 进一步作为优选的实施方式，所述步骤S4包括：
[0090] S41、各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的值来判断本 地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定本地上传的 WAMS量测有量测缺失，此时执行步骤S42;若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测无量测缺失，此时直接跳至步骤S43;
[0091] S42、迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计 值收敛至剔除缺失量测后的集中式估计值，然后执行步骤S43;
[0092] S43、将各子区域修正后的本地状态估计值进行输出。
[0093] 进一步作为优选的实施方式，所述步骤S42包括：
[0094] S421、给定收敛阔值eend，并为各区域i置状态估计值初值=韦.W及时刻初值t =0；
[009引 S422、各区域i根据时刻t的状态估计值卷,独立计算
和 ，并将巧韦,，的值赋予ai,0;
[0096] S423、各区域i根据赋值后的新ai,日执行有限时间平均一致性协议，得到
[0097] S424、各区域i根据
巧第二方程组独立计算时刻t的拉格朗日乘子λ*，并 计算各区域i在时刻t+1的状态估计值·%κ;ι，计算公式为：
[009引
[0099] S425、各区域i计算
然后置ei,〇 = ei，执行最大一致性协议， 所述最大一致性协议的表达式为：
[0100]
[0101] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0102] S426、根据最大一致性协议的执行结果判断< Ceud是否成立，若是，则结束迭 代，将时刻t的状态估计值作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤S422。
[0103] 参照图2,基于WAMS量测的多区域电力系统状态估计系统，包括W下模块：
[0104] 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状 态估计模型的最优解求解公式；
[0105] 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模 型估计本地状态，得到本地状态估计值；
[0106] 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行 的结果对本地状态估计值进行修正，所述分布式一致性协议W各子区域的电网调度中屯、作 为智能体，各个智能体之间通过链状或网状通信拓扑联通；
[0107] 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS 量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的 集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出；
[0108] 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布 式一致性协议执行模块进而与判断模块的输入端连接。
[0109] 参照图2,进一步作为优选的实施方式，所述多区域状态估计模型 构建模块包括：
[0110] 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：
[0111]
[0112] 其中，J为多区域状态估计模型的目标函数，N为多区域的总数，zi为区域i的WAMS 量测，Ai为区域i的量测方程矩阵，xi为区域i的状态矢量，Wi为区域i的权重矩阵；x = [XI，. . .，XN]为η维全系统状态矢量，B为4bXn阶矩阵，b为联络线支路数，B=[Bi，. . .，Bn]， 满足
[0113] 最优解求解单元，用于依次通过第一方程组、第二方程组和第Ξ方程组求出多区 域状态估计模型的最优解，所述第一方程组的表达式为
所述第二方程组的 表达式为：
，其中，疋为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，为区域i的最优状态估计值，dl为 区域i的估计协调增量；
[0114] 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端 连接。
[0115] 参照图2,进一步作为优选的实施方式，所述判断模块包括：
[0116] 判断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵Si对角元素的 值来判断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为0,则判定 本地上传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵Si对角元素的值 为1，则判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元；
[0117] 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后 的本地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元；
[0118] 输出单元，用于将各子区域修正后的本地状态估计值进行输出；
[0119] 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单 元的输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的 输出端还与输出单元的输入端连接。
[0120] 下面结合说明书附图和具体实施例对本发明作进一步详细说明。
[0121] 实施例一
[0122] 本实施例对本发明所设及的相关理论及本发明的具体实现过程进行详细介绍。
[0123] 针对传统WAMS具有庞大数据量、存在量测数据缺失的问题和新型通信拓扑（链状 或网状)的发展，本发明提出了一种针对WAMS量测缺失的多区域电力系统状态估计方法，如 图4所示，无需协调中屯、，能适应任意通信拓扑;在WAMS正常情况(WAMS量测无缺失)下，能快 速收敛至集中式估计值;在WAMS出现部分PMU量测缺失时，无需改变原有信息矩阵结构，只 需执行若干次分布式一致性协议能使各子区域本地估计值收敛至剔除缺失量测后的集中 式估计值。
[0124] (一)多区域电力系统状态估计模型及相关定义
[0。5]如图3所示的二分区互联电力系统，区域i的状态量矢i
定义 为：
[0126] a)区域i内部状态XI,int:区域i内不与其它区域相连接的节点状态量；
[0127] b)区域i内部边界状态XI,bint:区域i内与其它区域相连接的节点状态量；
[0128] C)区域i外部边界状态XI,bext:其它区域内与区域i相连接的节点状态量。
[0129 ]类似地，可定义出区域j的状态量矢量。
[0130] 本发明建立的多区域电力系统状态估计模型为：
[0131]
[013引其中，Z功区域i的量测;Wi为区域i的权重矩阵;B为4bXn阶矩阵，B每行中有两个 非零元素为1或-1，分别对应两个相邻区域的边界状态量；b为联络线支路数。令B = [Bi，...，Bn]，则有：
[0133]
[0134] 引入拉格朗日乘子λ，上述多区域电力系统状态估计模型的最优求解过程可通过 W下Ξ个方程组来实现：
[0138] 其中，為为区域i的最优状态估计值，di为区域i的估计协调增量。
[0139] 由W上Ξ个方程组可知，λ需要各区域联合求解。
[0140] 故最优求解过程的第二方程组可等效改写为：
[0141]
[0142] 各区域i仅已知Bi，即仅已知矩阵Β中关于区域i边界状态的部分列向量，而区域i 的初始状态设为：
[0145]分布式一致性协议可采用现有的有限时间平均一致性协议来实现。各区域执行分 布式一致性协议后可得到一致平均值
，然后再通过分布式加权最小二乘估 计算法即可得到本地状态估计值，因为为常数矩阵，初始化只需一次。
[0146 ](二)分布式加权最小二乘估计算法
[0147]本发明的分布式加权最小二乘估计算法的具体过程为：
[014引（1)进行初始化，为各区域i置一致性协议初始值为巧G; I巧'，然后各区域i执 行有限时间平均一致性协议得到^辟；
[0149] (2)各区域i根据!：&和(一）的第一方程组独立计算韦和度式，并将公,?,的值赋予 过i，〇;
[0150] (3)各区域i根据赋值后的新ai，〇执行有限时间平均一致性协议，得到
；
[0151] (4)各区域i根据
巧(一)的第二方程组独立计算拉格朗日乘子λ，然后由λ和 (一)的第Ξ方程组计算本地状态估计值4。
[0152] (S)WAMS量测缺失情况下的迭代修正算法
[0153] 定理设缺失量测剔除后的系统能保持可观性，给定初值：
[0154]
[015引由
爱代可收敛至：
[0156]
[0157] 其中，S为对角线矩阵，对应异常量测（即有量测缺失）的对角元素为0,对应正常量 巧。（即无量测缺失）的对角元素为l;Zm、Am和Rm分别为不包含异常量测的量测向量、量测方程 矩阵和量测协方差矩阵。
[0158] 根据定理，在(二）的算法步骤中引入对角线矩阵Si，则（二）中最优求解过程的第 一方程组第一式子修正为：
[0159]
[0160] 并置。,。=昕巧，巧Δ?,·.,)叩。Tapi为区域i的WAMS量测缺失标志，当区域i有WAMS量测 缺失时，Tapi = l;否则，Tapi = 0。执行完(二）的算法后，各区域i检查Tapi是否为0,如不为0， 则继续执行W下的迭代修正算法：
[0161] (1)给定收敛阔值eend，并为各区域i置状态估计值初值韦,0 =韦W及时刻初值t = 0；
[0162] (2)各区域i根据时刻t的状态估计值為独立计算 和马％， 并将巧..iy的值赋予ai,0;
[0163] (3)各区域i根据赋值后的新ai,〇执行有限时间平均一致性协议，得到
[0164] (4)各区域i根据
和(一）的第二方程组独立计算时刻t的拉格朗日乘子 λ*，并计算各区域i在时刻t+1的状态估计值，计算公式为：
[0165]
[0166] (5)各区域i计算
，然后置ei,〇 = ei，执行最大一致性协议，所 述最大一致性协议的表达式为：
[0167]
[0168] 其中，Ωι表示与区域i相邻的区域编号集合，山为简单无向图G(v，e)的直径，V为图 G(v，e)的顶点集合，ε为图G(v，e)的边集合；
[0169] (6)根据最大一致性协议的执行结果判断6,為是否成立，若是，则结束迭 代，将时刻t的状态估计值焉>,作为各子区域修正后的本地状态估计值进行输出，反之，则令 t自加1然后返回步骤(2)。
[0170] 与现有技术相比，本实施例具有W下优点：
[0171] (1)根据互联电网已有的分区方式，将子区域电网调度中屯、与分布式一致性协议 中的智能体概念相联系，并通过链状或网状通信拓扑联通，W实现无协调中屯、的多区域状 态估计，效率更高。
[0172] (2)在估计时刻，各子区域间通过与相邻区域交换边界信息，得到估计结果得到全 局协调，保证各子区域本地估计值与集中式估计相同，更加科学。
[0173] (3)建立各子区域电网间状态估计模型的禪合关系，然后基于该关系合理引入分 布式一致性协议来实现全局协调求解，更加全面。
[0174] (4)扩展推导出WAMS量测存在部分缺失情况下的修正算法，对WAMS量测缺失具有 较强的适应性，更加全面和准确。
[017引实施例二
[0176] 1)W图5所示7分区的IEEE 118节点标准测试系统为例，所有节点均配置PMU，相应 节点电压相量和与该节点相联的支路电流相量是可量测的。节点电压实部与虚部量测值为 潮流真值叠加标准差σ = 0.008的高斯白噪声，支路电流实部与虚部量测值为潮流真值叠加 标准差〇 = 0.0005的高斯白噪声，权重矩阵W取为。图6为对应图5系统的通信拓扑。 巧
[01W]根据先验知识与实际测试，经过有限时间平均一致性迭代，各区域节点电压相量 的实虚部均收敛至集中式估计值。
[0178] 假设各区域均出现一个电压相量实部的量测缺失，即全系统共有7个量测缺失，并 将收敛阔值eend设为0.0001。若采用本发明的迭代修正算法进行测试，则收敛时所需的迭代 步数心=3。表1为其它缺失数据情况下采用本发明的迭代修正算法的收敛 结果。从表1我们 可W看到，t。随着系统缺失量测总数的增多而增大，但增幅并不明显。一般地，WAMS可靠性 较高，一次采样中的缺失量测数相对较少，故采用本发明的算法能通过较少的迭代步数迭 代至收敛。
[0179] 表1WAMS量测缺失情况下采用迭代修正算法的收敛结果
[0180]
~w上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施胃 例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替 换，运些等同的变形或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 基于WAMS量测的多区域电力系统状态估计方法，其特征在于:包括以下步骤： 51、 建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式； 52、 各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得到本 地状态估计值； 53、 各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行修正， 所述分布式一致性协议以各子区域的电网调度中心作为智能体，各个智能体之间通过链状 或网状通信拓扑联通； 54、 判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情况下 的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计值， 反之，则直接将各子区域修正后的本地状态估计值进行输出。2. 根据权利要求1所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤Sl，其包括： 建立多区域状态估计模型，所述多区域状态估计模型的表达式为：其中，J为多区域状态估计模型的目标函数，N为多区域的总数，Z1为区域i的WAMS量测， Ai为区域i的量测方程矩阵，Xi为区域i的状态矢量，Wi为区域i的权重矩阵;x=[X1，. . .，XN] 为η维全系统状态矢量，B为4b X η阶矩阵，b为联络线支路数，B = [Bi，. . .，Bn]，满足依次通过第一方程组、第二方程组和第三方程组求出多区域状态估计模型的最优解， 所述第一方程组的表达式为：所述第二方程组的表达式为：所述第三方程组的表达式为：其中，g为区域i 未经修正的本地状态估计值，λ为拉格朗日乘子，为区域i的最优状态估计值，Cl1为区域i 的估计协调增量。3. 根据权利要求2所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S2包括： 521、 各子区域获取当前时刻本地上传的WAMS量测； 522、 各子区域根据获取的WAMS量测以及多区域状态估计模型采用分布式最小二乘估 计算法估计本地状态，得到本地状态估计值。4. 根据权利要求3所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S22包括： 5221、 进行初始化，为各区域i置一致性协议初始值为然后各区域i执行 有限时间平均一致性协议得到; 5222、 各区域i根据告Gs和第一方程组独立计算笔和尽足，并将尽矣的值赋予ai,〇; 5223、 各区域i根据赋值后的新〇执行有限时间平均一致性协议，得到5224、 各区域i根据和第二方程组独立计算拉格朗日乘子λ，然后由λ和第三方 程组计算本地状态估计值· 〇5. 根据权利要求2所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S3包括： 各子区域执行分布式一致性协议，得到各子区域的估计协调增量； 各子区域根据得到的估计协调增量对本地状态估计值进行修正。6. 根据权利要求4所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S4包括： 541、 各子区域i根据本地上传的WAMS量测中对角线矩阵S1对角元素的值来判断本地上 传的WAMS量测是否有量测缺失，若对角线矩阵S i对角元素的值为1，则判定本地上传的WAMS 量测有量测缺失，此时执行步骤S42;若对角线矩阵Sdt角元素的值为1，则判定本地上传的 WAMS量测无量测缺失，此时直接跳至步骤S43; 542、 迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收 敛至剔除缺失量测后的集中式估计值，然后执行步骤S43; 543、 将各子区域修正后的本地状态估计值进行输出。7. 根据权利要求6所述的基于WAMS量测的多区域电力系统状态估计方法，其特征在于： 所述步骤S42包括： 5421、 给定收敛阈值eend，并为各区域i置状态估计值初值毛，。=天以及时刻初值t = 0; 5422、 各区域i根据时刻t的状态估计值&.^虫立计算并将尽足7的值赋予ai,o; 5423、 各区域i根据赋值后的新〇执行有限时间平均一致性协议，得到5424、 各区域i根据和第二方程组独立计算时刻t的拉格朗日乘子&，并计 算各区域i在时刻t+1的状态估计值7.+|计算公式为：5425、 各区域i计算_然后置ei,〇 = ei，执行最大一致性协 议，所述最大一致性协议的表达式为：其中，Q1表示与区域i相邻的区域编号集合，dc为简单无向图G(v，〇的直径，V为图G(v， O的顶点集合，ε为图G(v，〇的边集合； 5426、 根据最大一致性协议的执行结果判断是否成立，若是，则结束迭 代，将时刻t的状态估计值作为各子区域修正后的本地状态估计值进行输出，反之，则 令t自加1然后返回步骤S422。8. 基于WAMS量测的多区域电力系统状态估计系统，其特征在于:包括以下模块： 多区域状态估计模型构建模块，用于建立多区域状态估计模型，并获得多区域状态估 计模型的最优解求解公式； 本地状态估计模块，用于各子区域根据本地上传的WAMS量测和多区域状态估计模型估 计本地状态，得到本地状态估计值； 分布式一致性协议执行模块，用于各子区域执行分布式一致性协议，并根据执行的结 果对本地状态估计值进行修正，所述分布式一致性协议以各子区域的电网调度中心作为智 能体，各个智能体之间通过链状或网状通信拓扑联通； 判断模块，用于判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测 缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中 式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出； 所述多区域状态估计模型构建模块的输出端依次通过本地状态估计模块和分布式一 致性协议执行模块进而与判断模块的输入端连接。9. 根据权利要求8所述的多区域电力系统状态估计系统，其特征在于:所述多区域状态 估计模型构建模块包括： 建立单元，用于建立多区域状态估计模型，所述多区域状态估计模型的表达式为：其中，J为多区域状态估计模型的目标函数，N为多区域的总数，Z1为区域i的WAMS量测， Ai为区域i的量测方程矩阵，Xi为区域i的状态矢量，Wi为区域i的权重矩阵;X=[X1，. . .，XN] 为η维全系统状态矢量，B为4b X η阶矩阵，b为联络线支路数，B = [Bi，. . .，Bn]，满足最优解求解单元，用于依次通过第一方程组、第二方程组和第三方程组求出多区域状 态估计模型的最优解，所述第一方程组的表达式为：，所述第二方程组的表达 式为：所述第三方程组的表达式为：其中，天为 区域i未经修正的本地状态估计值，λ为拉格朗日乘子，$为区域i的最优状态估计值，Cl1为 区域i的估计协调增量； 所述建立单元的输出端通过最优解求解单元进而与本地状态估计模块的输入端连接。10.根据权利要求8所述的多区域电力系统状态估计系统，其特征在于:所述判断模块 包括： 判断单元，用于各子区域i根据本地上传的WAMS量测中对角线矩阵S1对角元素的值来判 断本地上传的WAMS量测是否有量测缺失，若对角线矩阵Si对角元素的值为1，则判定本地上 传的WAMS量测有量测缺失，此时执行迭代修正单元;若对角线矩阵SM角元素的值为1，则 判定本地上传的WAMS量测无量测缺失，此时直接跳至输出单元； 迭代修正单元，用于迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本 地状态估计值收敛至剔除缺失量测后的集中式估计值，然后执行输出单元； 输出单元，用于将各子区域修正后的本地状态估计值进行输出； 所述分布式一致性协议执行模块的输出端与判断单元的输入端连接，所述判断单元的 输出端分别与迭代修正单元的输入端和输出单元的输入端连接，所述迭代修正单元的输出 端还与输出单元的输入端连接。
【专利摘要】本发明公开了一种基于WAMS量测的多区域电力系统状态估计方法及系统，方法包括：建立多区域状态估计模型，并获得多区域状态估计模型的最优解求解公式；各子区域根据本地上传的WAMS量测和多区域状态估计模型估计本地状态，得到本地状态估计值；各子区域执行分布式一致性协议，并根据执行的结果对本地状态估计值进行修正；判断本地上传的WAMS量测是否有量测缺失，若是，则迭代执行WAMS量测缺失情况下的迭代修正算法，直到修正后的本地状态估计值收敛至剔除缺失量测后的集中式估计值，反之，则直接将各子区域修正后的本地状态估计值进行输出。本发明具有收敛速度快、适用性广且准确的优点，可广泛应用于电力系统分析控制领域。
【IPC分类】G06F19/00, H02J3/00
【公开号】CN105490269
【申请号】CN201511028254
【发明人】陈亦平, 蔡永智, 侯君, 陈皓勇, 杨林, 郭曼兰, 王巍, 万楚林, 刘起兴
【申请人】中国南方电网有限责任公司, 华南理工大学
【公开日】2016年4月13日
【申请日】2015年12月30日
最新回复(0)