一种基于两阶段法的主配网一体化分布式状态估计方法
【专利摘要】本发明涉及一种基于两阶段法的主配网一体化分布式状态估计方法，方法包括：建立主配网一体化状态估计模型；将主配网一体化状态估计模型的配网模型解耦为若干子系统；建立解耦前主配网的状态量与解耦后子系统中状态量的线性关系；对解耦后子系统以分布式的方式独立进行状态估计；以解耦前后主配网状态量的线性关系为量测方程，建立状态估计模型进行状态估计，获取主配网的一体化状态。本发明提出的方法降低了计算量，有效且准确地解决了由于配网网络规模庞大造成的主配网一体化状态估计维数灾的问题，而且能够得到主配网一致的、精确的实时状态，为主配网一体化的分析与计算提供基础数据，进而提高了主配网一体化的运行可靠性。
【专利说明】
-种基于两阶段法的主配网一体化分布式状态估计方法
技术领域
[0001] 本发明设及配电自动化领域，具体设及一种基于两阶段法的主配网一体化分布式 状态估计方法。
【背景技术】
[0002] 分布式电源与常规电源相结合的方式将是电力系统发展的趋势，大量分布式电源 的接入，改变了配电网的特性，采用基于边界等值实现主网和配网相互独立运行分析控制 的方法已不适用，需要从主配网一体化全局考虑，对主配网进行一体化分析。
[0003] 分布式状态估计利用状态估计算法本身可分区域计算的特性，在各分区进行本区 域的状态估计，再通过与其他区域的信息交换，从而使所有分区的状态估计结果达到或接 近整体计算的效果，相对于集中式状态估计，不存在因增益矩阵的条件数过大、数值稳定性 差等问题导致的整体状态估计失败等问题，研究的核屯、在于如何协调各个子分区状态估计 结果，W期获得与集中式估计相似的计算精度。
[0004] 目前，国内外已有许多专家学者对分布式状态估计进行研究，并取得了大量成果。 按照各分区状态估计协调阶段的不同，可W分为迭代收敛后协调的分布式状态估计和迭代 期间协调的分布式状态估计两种。
[0005] (1)迭代收敛后协调的分布式状态估计最大的优点是能够将残差污染控制在本区 域内，各区域的收敛性互不影响，因此，该种分布式状态估计研究比较多，其核屯、在于如何 协调各个子分区状态估计的估计结果，W期获得与集中式估计相近的计算精度，各个分区 的协调机制对算法的计算量与计算结果影响较大，一般各区域需要反复协调，但是该种分 布式状态估计对边界节点的注入量测处理较为粗糖。
[0006] (2)迭代期间协调的分布式状态估计在各个区域每次迭代后，均要进行一次统一 协调，其优点是反复协调可W使得解更趋近最优解，缺点是协调侧与子分区或子分区之间 需要频繁地交换数据，无法使残差污染控制在本区域内，容易造成残差污染的传播与扩大， 各区域间的收敛性容易相互影响。
[0007] 主配网一体化实时分析和计算对数据断面的基本要求是全网状态量的一致性，主 配独立进行状态估计的最大问题是得到的全网状态量不一致，比如主配网互相的等值注入 的估计结果不匹配，因此有必要对主配网进行一体化状态估计，由于主配网的网络特性差 异较大、量测覆盖率相差较大，主网只需单相分析、配网需要=相分析，如果将主配网作为 一个整体进行状态估计，在数学上收敛性无法保证，况且主网的状态估计精度也会因配网 影响大大降低，所W将主配网的实时量测整合到一个系统中，进行一体化的分布式状态估 计尤为必要。
【发明内容】
[000引有鉴于此，本发明提供的一种基于两阶段法的主配网一体化分布式状态估计方 法，该方法降低了计算量，有效且准确地解决了由于配网网络规模庞大造成的主配网一体 化状态估计维数灾的问题，而且能够得到主配网一致的、精确的实时状态，为主配网一体化 的分析与计算提供基础数据，进而提高了主配网一体化的运行可靠性。
[0009] 本发明的目的是采用下述技术方案实现的：
[0010] -种基于两阶段法的主配网一体化分布式状态估计方法，其改进之处在于，包括：
[0011] 步骤1.建立主配网一体化状态估计模型；
[0012] 步骤2.在所述主配网一体化状态估计模型中添加零阻抗支路及虚拟母线，将所述 主配网一体化状态估计模型的配网模型解禪为若干子系统；
[0013] 步骤3.建立解禪前所述主配网的状态量与解禪后所述子系统中状态量的线性关 系；
[0014] 步骤4.第一阶段状态估计:对解禪后子系统W分布式的方式独立进行状态估计；
[0015] 步骤5.第二阶段状态估计：W解禪前后主配网状态量的线性关系为量测方程，建 立状态估计模型进行状态估计，获取主配网的一体化状态。
[0016] 优选的，所述步骤1中，所述主配网一体化状态估计模型包括：主网模型和配网模 型；
[0017] 其中，所述主网模型为单相模型，由主网高压母线、变压器和主网低压母线组成， 所述主网高压母线的出线依次连接所述变压器和主网低压母线；
[0018] 所述配网模型为=相模型，由馈线组成，所述馈线为所述主网低压母线的出线。
[0019] 优选的，按下式确定所述主配网一体化状态估计模型的量测方程Z:
[0020] Z = H(Xi)姑)+e (1)
[0021] 式（1)中，Z为主配网全网的量测向量，H为量测函数，e为量测误差向量，&为除边 界节点之外的节点状态向量，姑为边界节点的状态量。
[0022] 优选的，所述步骤2包括：
[0023] 将主配网一体化状态估计模型中配网模型的馈线与低压母线之间添加零阻抗支 路，并在所述馈线侧添加所述馈线对应的虚拟母线，由所述馈线及所述馈线对应的虚拟母 线构成所述子系统；
[0024] 其中，所述零阻抗支路为所述子系统之间的联络线，按下式确定解禪后所述子系 统的内部状态量与解禪前主配网的关系：
[002引 X; =Xi 12)
[0026] 式(2)中，馬为解禪后各子系统内部的状态向量;X功解禪前主配网的内部变量；
[0027] 优选的，所述步骤3包括:选择主网中非边界母线作为解禪前主配网全系统及解禪 后主网的相角参考点，选择虚拟母线作为解禪后所述子系统的相角参考点，令解禪前主配 网节点数为N，边界母线数为M，总馈线数为K，并按下式确定解禪前所述主配网的状态量Y与 解禪后所述子系统中状态量X的线性关系：
[0028]
(3>
[0029] 式(3)中，C为常数;Xb为解禪前主配网边界母线的2M维状态列向量，Yi为解禪后子 系统内部状态量2(N-M)-1维度列向量，N为解禪前主配网节点数，M为解禪前主配网边界母 线数，Yb为解禪后子系统边界状态量组成的4M维列向量，I为单位矩阵，Bi为除对应边界母线 相角元素为I外均为零的矩阵，Bb为除对应边界母线幅值的位置为I外均为零的矩阵。
[0030] 优选的，所述步骤4包括：
[0031] 确定解禪后主配网的量侧方程Zk,公式为：
[0032] zk = fk(yi,k,yb,k)+ek (4)
[0033] 式(4)中，Zk为子系统k的量测向量，fk为子系统k的量测函数，ek为子系统k的量测 误差，yi,k为主配网解禪子系统k的内部各节点变量，i为节点编号;yb,k为主配网解禪子系统 k的边界节点变量；
[0034] 其中，基于加权最小二乘状态估计方程，式(4)的求解公式为：
[0035]
(5)
[0036] 式(5)中，巧卸?，是子系统k的雅可比矩阵;j是迭代次数;巧:!=巧听F,.子系 统k的信息矩阵;卸P为子系统k第j次迭代时状态量的修正量;Wk为加权矩阵;yi,k为子系统 k的内部节点状态向量;yb,功子系统k的边界节点状态向量间为第j次迭代时的量测向量。
[0037] 优选的，所述步骤5包括：
[0038] 考虑所述子系统状态量的估计误差，按下式确定解禪前所述主配网的状态量Y与 解禪后所述子系统中状态量X的线性关系：
[0039] Y = CX 巧 （6)
[0040] 式(6)中，E为子系统状态量的误差向量。
[0041] 进一步的，当不考虑约束时，基于加权最小二乘算法，获取解禪后所述子系统中状 态量X，公式为：
[0042]
(7)
[0043] 式（7)中，；是X的估计值，Wy为Y的量测权矩阵，;，为Y的估计值，Wy为的Y量测权矩 阵；
[0044] 其中，Y量测权矩阵Wy的公式为：
[0045] Wy = COV'(V) = (8)
[0046] 式(8)中，晏解禪后子系统状态量的估计值，如为根据y计算出的信息矩阵；
[0047] 将Gf表示为内部变量与外部变量的形式，公式为：
[004引
（9)
[0049] 式（9)中，Gii、Gib、G至及Gbb均为Gf中各兀素；
[0050] 将式（7)进一步表示为：
[0051 ]
(IO)
[0052]式（10)中，Bi与Bb为C中的元素，I为单位阵，克为第二阶段状态估计对各子系统内 部状态量的估计值，；6为第二阶段状态估计对各子系统边界状态量的估计值，；,.为第一阶 段状态估计对各子系统内部状态量的估计值，为第一阶段状态估计对各子系统边界状态 量的估计值。
[0053] 进一步的，当考虑零注入约束时，求解零注入约束，公式为：
[0化4]
(11》
[0055] 式（11)中，P(X)为边界母线的零注入约束；
[0056] 添加拉格朗日乘子处理式（11 )，则有：
[0化7]
(12)
[0化引对式（12)中的Y与A分别求极值得，并对极值条件在是处展开，则有：
[0化9]
(13)
[0060] 式(13)中，X为待估计的状态量，A为拉格朗日乘子，P为等式约束的雅可比矩阵；
[0061] 消去式(13)中的Ax,求解入：
[0062]
(14)
[00创进而获取Ax:
[0064] A X = -(C^yC)"Va (15)。
[0065] 进一步的，当边界母线存在等值的馈线时，注入量测为：
[0066] Zp = p(x)+Gp (16)
[0067] 式（16)中，却为原系统边界母线的等值注入量测，P(X)为量测函数，ep为量测误差；
[0068] 构造抢格朗日函数，公式为：
[0069]
(\1)
[0070] 式(17)中，L(Y)为目标函数，X为待估计的状态量，A为拉格朗日乘子；
[0071 ]在X处进行泰勒公式展开，公式为：
[0072]
C18)
[0073] 式（18)中，A X为状态量的修正量，P为等式约束的雅可比矩阵，Wp为边界量测的量 测权矩阵为边界量测的量测函数；
[0074] 消除A X，得到A的求解公式：
[0075]
(19)
[0076] 根据求得的A计算A X:
[0077] A X = -(C^yC)"Va (20)
[0078] 从而得到满足式的解为：
[0079] X = X +At (21)。
[0080] 从上述的技术方案可W看出，本发明提供了一种基于两阶段法的主配网一体化分 布式状态估计方法，建立主配网一体化状态估计模型；在主配网一体化状态估计模型中添 加零阻抗支路及虚拟母线;根据零阻抗支路，将主配网解禪为多个子系统;建立主配网的状 态量与各子系统中状态量的线性关系;进行第一阶段估计及第二阶段估计，得到主配网的 一体化状态。本发明提出的方法降低了计算量，有效且准确地解决了由于配网网络规模庞 大造成的主配网一体化状态估计维数灾的问题，而且能够得到主配网一致的、精确的实时 状态，为主配网一体化的分析与计算提供基础数据，进而提高了主配网一体化的运行可靠 性。
[0081] 本发明的有益效果：
[0082] 1、本发明所提供的技术方案中，降低了计算量，有效且准确地解决了由于配网网 络规模庞大造成的主配网一体化状态估计维数灾的问题。
[0083] 2、本发明所提供的技术方案，主网、配网各馈线可W分别采用单相模型、=相模 型；能够得到主配网一致的、精确的实时状态，为主配网一体化的分析与计算提供基础数 据。
[0084] 3、本发明所提供的技术方案，坏数据处理在各子系统内部处理完成;提高了主配 网一体化的运行可靠性。
[0085] 4、本发明提供的技术方案，应用广泛，具有显著的社会效益和经济效益。
【附图说明】
[0086] 图1是本发明一种基于两阶段法的主配网一体化分布式状态估计方法的流程图；
[0087] 图2是本发明实施例中解禪前主配网一体化状态估计模型示意图；
[0088] 图3是本发明实施例中解禪后主配网一体化状态估计模型示意图；
[0089] 图4是本发明实施例中第二阶段状态估计流程图。
【具体实施方式】
[0090] 下面结合附图对本发明的【具体实施方式】作详细说明。
[0091] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0092] 术语及定义：
[0093] 1)状态估计:状态估计也称滤波，它是利用实时量测系统的冗余度来提高数据精 度，排除随机干扰所引起的错误信息，估计或预报系统的运行状态。电力系统状态估计分为 静态和动态状态估计。静态状态估计目前应用较为成熟，W最小二乘法等为主，静态状态估 计的局限性是没有考虑系统的动态。
[0094] 2)分布式状态估计：分布式状态估计利用状态估计算法本身可分区域计算的特 性，在各分区进行本区域的状态估计，再通过与其他区域的信息交换，从而使所有分区的状 态估计结果达到或接近整体计算的效果，相对于集中式状态估计，不存在因增益矩阵的条 件数过大、数值稳定性差等问题导致的整体状态估计失败等问题。
[0095] 3)主配网：本专利中指的是目前我国电网调度运行管辖权的规定，W220kV、 110kV、35kV变电站为界，上述变电站及电压等级高于或等于上述变电站的电网，称为主网， IOkV或6.3kV变电站与线路称为配网。
[0096] 4)零阻抗支路:零阻抗支路是指支路阻抗、对地导纳均为零的支路，其特性为支路 两端节点电压幅值与相角均相等的支路，在电力系统零阻抗支路是真实存在，比如开关中 的母联开关、同一电厂的新建机组与原有机组高压母线之间的联络线，由于电气距离只有 几十米，电气特性上也相当于零阻抗支路。
[0097] 本发明提供的一种基于两阶段法的主配网一体化分布式状态估计方法，如图1所 示，包括：
[0098] 步骤1.建立主配网一体化状态估计模型；
[0099] 步骤2.在所述主配网一体化状态估计模型中添加零阻抗支路及虚拟母线，将所述 主配网一体化状态估计模型的配网模型解禪为若干子系统；
[0100] 步骤3.建立解禪前所述主配网的状态量与解禪后所述子系统中状态量的线性关 系；
[0101] 步骤4.第一阶段状态估计:对解禪后子系统W分布式的方式独立进行状态估计；
[0102] 步骤5.第二阶段状态估计：W解禪前后主配网状态量的线性关系为量测方程，建 立状态估计模型进行状态估计，获取主配网的一体化状态。
[0103] 具体的，如图2所示，所述步骤1中，所述主配网一体化状态估计模型包括：主网模 型和配网模型；
[0104] 其中，所述主网模型为单相模型，由主网高压母线、变压器和主网低压母线组成， 所述主网高压母线的出线依次连接所述变压器和主网低压母线；
[0105] 所述配网模型为=相模型，由馈线组成，所述馈线为所述主网低压母线的出线。
[0106] 按下式确定所述主配网一体化状态估计模型的量测方程Z:
[0107] Z = H(Xi)姑)+e (1)
[010引式（1)中，Z为主配网全网的量测向量，H为量测函数，e为量测误差向量，&为除边 界节点之外的节点状态向量，姑为边界节点的状态量。
[0109] 其中，
[0110] 在图2中的配网各馈线与主网母线之间添加零阻抗支路，并为每条馈线添加虚拟 母线，如图3所示，所述步骤2包括：
[0111] 将主配网一体化状态估计模型中配网模型的馈线与低压母线之间添加零阻抗支 路，并在所述馈线侧添加所述馈线对应的虚拟母线，由所述馈线及所述馈线对应的虚拟母 线构成所述子系统，其中，所述零阻抗支路为所述子系统之间的联络线；
[0112] 令添加零阻抗后Bl、B2母线的状态量分别为Xb, 1与Xb, 2,母线B3、B4的状量分别为 ^1,3与^1,4，则与原状态量的关系如下：
[0113]
[0114] 通过增加零阻抗支路与对主网边界母线增加新状态量，实现了主网与配网各馈线 之间的解禪，则按下式确定解禪后所述子系统的内部状态量与解禪前主配网的关系：
[0115] X'；^Xj (2)
[0116] 式(2)中，iTf为解禪后各子系统内部的状态向量;Xi为解禪前主配网的内部变量；
[0117]通过添加零阻抗支路，增广了边界母线的状态量，将原系统解禪为主网、馈线I与 馈线2=个子系统，选择主网中非边界母线为解禪前主配网全系统的相角参考点，并且也作 为解禪后主网的相角参考点，即对于主网解禪前后节点相角参考点不变，解禪后选择母线 B3、B4分别为馈线1与馈线2子系统的相角参考点，则解禪后子系统状态量与原系统状态量 的关系为：
[011 引
[0119] 基于上式解禪后子系统状态量与原系统状态量的关系，选择主网中非边界母线作 为解禪前主配网全系统及解禪后主网的相角参考点，选择虚拟母线作为解禪后所述子系统 的相角参考点，令解禪前主配网节点数为N，边界母线数为M，总馈线数为K，并按下式确定解 禪前所述主配网的状态量Y与解禪后所述子系统中状态量X的线性关系：
[0120]
扮）
[0121] 式(3)中，C为常数;Xb为解禪前主配网边界母线的2M维状态列向量，Yi为解禪后子 系统内部状态量2(N-M)-1维度列向量，N为解禪前主配网节点数，M为解禪前主配网边界母 线数，Yb为解禪后子系统边界状态量组成的4M维列向量，I为单位矩阵，Bi为除对应边界母线 相角元素为1外均为零的矩阵，Bb为除对应边界母线幅值的位置为1外均为零的矩阵。
[0122] 所述步骤4包括：
[0123] 确定解禪后主配网的量侧方程Zk,公式为：
[0124] zk = fk(yi,k,yb,k)+ek (4)
[012引式(4)中，Zk为子系统k的量测向量，fk为子系统k的量测函数，ek为子系统k的量测 误差，yi,k为主配网解禪子系统k的内部各节点变量，i为节点编号;yb,k为主配网解禪子系统 k的边界节点变量；
[0126]其中，基于加权最小二乘状态估计方程，式(4)的求解公式为：
[0127]
巧)
[012引式(5)中，巧"=%/如,是子系统k的雅可比矩阵;視迭代次数;巧巧巧A子系 统k的信息矩阵；A乃W为子系统k第j次迭代时状态量的修正量;恥为加权矩阵;yi,k为子系统 k的内部节点状态向量;yb,功子系统k的边界节点状态向量间为第j次迭代时的量测向量。
[0129] 所述步骤5,如图4所示，包括：
[0130] 考虑所述子系统状态量的估计误差，按下式确定解禪前所述主配网的状态量Y与 解禪后所述子系统中状态量X的线性关系：
[0131] Y = CX 巧 （6)
[0132] 式(6)中，E为子系统状态量的误差向量。
[0133] 当不考虑约束时，基于加权最小二乘算法，获取解禪后所述子系统中状态量X，公 式为：
[0134] (7)
[013引式（7)中，；是X的估计值，町为￥的量测权矩阵，J为Y的估计值，Wy为的Y量测权矩 阵；
[0136] 其中，Y量测权矩阵Wy的公式为：
[0137]
(8)
[013引式(8)中，;，解禪后子系统状态量的估计值，如为根据；计算出的信息矩阵；
[0139] 将Gf表示为内部变量与外部变量的形式，公式为：
[0140]
(9)
[0141 ] 式（9)中，Gii、Gib、巧及G化均为如中各兀素；
[0142] 悠井(^7')讲一击亲呆责.
[01创
（脱
[0144] 式（10)中，Bi与化为C中的元素，I为单位阵，占为第二阶段状态估计对各子系统内 部状态量的估计值，；为第二阶段状态估计对各子系统边界状态量的估计值，，^，,为第一阶 段状态估计对各子系统内部状态量的估计值，：?&为第一阶段状态估计对各子系统边界状态 量的估计值。
[0145] 其中，式（10)的方程维数庞大，求解工作量很大，可W采用分步求解的方式，先求 解1? ,后求解JO。&的求解公式如下：
[0146]
[0147] 根据上式求解得到^后，将带入式(10)后，得：
[0148] Xi.=.AXi + Vi -B,Xb
[0149] 其中A Xi由下式求得：
[0150]
[0151] 由于Gii是对角分块的形式，则对于A Xi的求解可W分块进行，即：
[0152]
[0153] 下标k表示第k个子系统。
[0154] 当主网与配网的边界当没有对馈线进行等值时，考虑零注入约束，求解零注入约 束，公式为：
[0155]
ClU
[0156] 式（11)中，P(X)为边界母线的零注入约束；
[0157] 添加拉格朗日乘子处理式(11 )，则有：
[015 引
（12)
[0159] 对式(12)中的Y与A分别求极值得，并对极值条件在;^处展开，则有：
[0160]
(13)
[0161] 式(13)中，X为待估计的状态量，A为拉格朗日乘子，P为等式约束的雅可比矩阵；
[0162] 消去式(13)中的Ax,求解入：
[0163]
(14)
[0164] 进而获取Ax:
[01 化]Ax = -(CTWYC)-ipTA (15)。
[0166] 当边界母线存在等值的馈线时，考虑馈线等值注入，注入量测为：
[0167] Zp = p(x)+Gp (16)
[0168] 式（16)中，却为原系统边界母线的等值注入量测，P(X)为量测函数，ep为量测误差；
[0169] 构造拉格朗日函数，公式为：
[0170]
(17)
[0171] 式(17)中，L(Y)为目标函数，X为待估计的状态量，A为拉格朗日乘子；
[0172] 在x々h讲行泰勒公式展巧，公式为：
[017引
《戚)
[0174]式（18)中，A X为状态量的修正量，P为等式约束的雅可比矩阵，Wp为边界量测的量 测权矩阵，为边界量测的量测函数；
[01巧]消除A X，得到A的求解公式：
[0176]
(19)
[0177] 根据求得的A计算Ax:
[017 引 Ax = -(CTWYC)-ipT 入 （20)
[0179]从而得到满足式的解为：
[0180] X=;+ Ajc (21)。
[0181] 最后应当说明的是：W上实施例仅用W说明本发明的技术方案而非对其限制，尽 管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解:依然 可W对本发明的【具体实施方式】进行修改或者等同替换，而未脱离本发明精神和范围的任何 修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
【主权项】
1. 一种基于两阶段法的主配网一体化分布式状态估计方法，其特征在于，所述方法包 括： 步骤1.建立主配网一体化状态估计模型； 步骤2.在所述主配网一体化状态估计模型中添加零阻抗支路及虚拟母线，将所述主配 网一体化状态估计模型的配网模型解禪为若干子系统； 步骤3.建立解禪前所述主配网的状态量与解禪后所述子系统中状态量的线性关系； 步骤4.第一阶段状态估计:对解禪后子系统W分布式的方式独立进行状态估计； 步骤5.第二阶段状态估计：W解禪前后主配网状态量的线性关系为量测方程，建立状 态估计模型进行状态估计，获取主配网的一体化状态。2. 如权利要求1所述的方法，其特征在于，所述步骤1中，所述主配网一体化状态估计模 型包括:主网模型和配网模型； 其中，所述主网模型为单相模型，由主网高压母线、变压器和主网低压母线组成，所述 主网高压母线的出线依次连接所述变压器和主网低压母线； 所述配网模型为Ξ相模型，由馈线组成，所述馈线为所述主网低压母线的出线。3. 如权利要求2所述的方法，其特征在于，按下式确定所述主配网一体化状态估计模型 的量测方程Z: Z = H(Xi，姑)+e (1) 式（1)中，Z为主配网全网的量测向量，Η为量测函数，e为量测误差向量，X功除边界节点 之外的节点状态向量，姑为边界节点的状态量。4. 如权利要求1所述的方法，其特征在于，所述步骤2包括： 将主配网一体化状态估计模型中配网模型的馈线与低压母线之间添加零阻抗支路，并 在所述馈线侧添加所述馈线对应的虚拟母线，由所述馈线及所述馈线对应的虚拟母线构成 所述子系统； 其中，所述零阻抗支路为所述子系统之间的联络线，按下式确定解禪后所述子系统的 内部状态量与解禪前主配网的关系： X'主三文， (2 ) 式(2)中，X}'为解禪后各子系统内部的状态向量;&为解禪前主配网的内部变量。5. 如权利要求1所述的方法，其特征在于，所述步骤3包括:选择主网中非边界母线作为 解禪前主配网全系统及解禪后主网的相角参考点，选择虚拟母线作为解禪后所述子系统的 相角参考点，令解禪前主配网节点数为N，边界母线数为M，总馈线数为K，并按下式确定解禪 前所述主配网的状态量Y与解禪后所述子系统中状态量X的线性关系：(3) 式(3)中，C为常数;Xb为解禪前主配网边界母线的2M维状态列向量，Υι为解禪后子系统 内部状态量2(Ν-Μ)-1维度列向量，Ν为解禪前主配网节点数，Μ为解禪前主配网边界母线数， Υβ为解禪后子系统边界状态量组成的4Μ维列向量，I为单位矩阵，Βι为除对应边界母线相角 元素为1外均为零的矩阵，Ββ为除对应边界母线幅值的位置为1外均为零的矩阵。6. 如权利要求1所述的方法，其特征在于，所述步骤4包括： 确定解禪后主配网的量侧方程Zk，公式为： zk = fk(yi,k,yb,k)+ek (4) 式(4)中，Zk为子系统k的量测向量，fk为子系统k的量测函数，ek为子系统k的量测误差， yi,k为主配网解禪子系统k的内部各节点变量，i为节点编号;yb,k为主配网解禪子系统k的边 界节点变量； 其中，基于加权最小二乘状态估计方程，式(4)的求解公式为：(：5) 式(5)中，巧= % /却是子系统k的雅可比矩阵；j是迭代次数;巧;=巧的巧子系统k 的信息矩阵;Δ乂"为子系统k第j次迭代时状态量的修正量;Wk为加权矩阵;yi,k为子系统k的 内部节点状态向量;yb,功子系统k的边界节点状态向量间为第j次迭代时的量测向量。7. 如权利要求1所述的方法，其特征在于，所述步骤5包括： 考虑所述子系统状态量的估计误差，按下式确定解禪前所述主配网的状态量Y与解禪 后所述子系统中状态量X的线性关系： Y=CX 巧（6) 式(6)中，E为子系统状态量的误差向量。8. 如权利要求7所述的方法，其特征在于，当不考虑约束时，基于加权最小二乘算法，获 取解禪后所述子系统中状态量X，公式为：(7) 式(7)中，；是X的估计值，Wy为Y的量巧贩矩阵，；为Y的估计值，Wy为的Y量测权矩阵； 其中，Y量测权矩阵Wy的公式为：(8) 式(8)中，解禪后子系统状态量的估计值，如为根据；计算出的信息矩阵； 将Gf表不为内部变量与外部变量的形式，公式为：(9) 式（9)中，G"、G化、磅及Gbb均为如中各元素； 将式(7)进一步表示为：(10) 式（10)中，Bi与Bb为C中的元素，I为单位阵，品为第二阶段状态估计对各子系统内部状 态量的估计值，三6为第二阶段状态估计对各子系统边界状态量的估计值，为第一阶段状 态估计对各子系统内部状态量的估计值，为第一阶段状态估计对各子系统边界状态量的 估计值。9. 如权利要求8所述的方法，其特征在于，当主网与配网的边界没有对馈线进行等值时 考虑零注入约束，求解零注入约束，公式为：(11) 式（11)中，p(x)为边界母线的零注入约束； 添加拉格朗日乘子处理式(11 )，则有：(12) 对式(12)中的Y与λ分别求极值得，并对极值条件在;^处展开，则有：(13) 式(13)中，X为待估计的状态量，λ为拉格朗日乘子，Ρ为等式约束的雅可比矩阵； 消去式(13)中的Δχ，求解λ:10.如权利要求9所述的方法，其特征在于，当边界母线存在等值的馈线时，考虑馈线等 值注入，注入量测为： Zp = p(x)+ep (16) 式（16)中，zp为原系统边界母线的等值注入量测，p(x)为量测函数，ep为量测误差； 构造拉格朗日函数，公式为：(17) 式(17)中，L(Y)为目标函数，X为待估计的状态量，λ为拉格朗日乘子； 在X处进行泰勒公式展开，公式为：(18) 式（18)中，Δ X为状态量的修正量，Ρ为等式约束的雅可比矩阵，Wp为边界量测的量测权 矩阵，^(λ·)为边界量测的量测函数； 消除Αχ,得到λ的求解公式：
【文档编号】H02J3/00GK105977963SQ201610279481
【公开日】2016年9月28日
【申请日】2016年4月28日
【发明人】王少芳, 郎燕生, 赵昆, 李理, 宋旭日, 王淼, 罗雅迪
【申请人】中国电力科学研究院, 国家电网公司