一种计及多种网络结构的配电网量测配置方法

文档序号:10596436阅读:598来源:国知局
一种计及多种网络结构的配电网量测配置方法
【专利摘要】本发明公开了一种计及多种网络结构的配电网量测配置方法,包括预估配电网的运行状态,确定量测量与状态量的关系;根据配电网实际的运行状态,在量测配置时构建多种不同状态下的网络结构;分析配电网系统在不同网络结构下的饱和量测特性,基于启发式算法确定饱和量测数量;在饱和量测数量为量测数量的安装数量上限的约束条件下,以多种网络结构估计误差的加权和最小为目标,构建量测配置模型;利用层次分析法求解量测配置模型的权重系数,遍历所有可选方案,求解最优解,输出最优量测配置方案。本发明以权重系数的方式分析了不同网络结构对量测配置的影响差异,量测配置结果更切实际。
【专利说明】
一种计及多种网络结构的配电网量测配置方法
技术领域
[0001] 本发明涉及一种计及多种网络结构的配电网量测配置方法。
【背景技术】
[0002] 智能配电网高级应用软件技术及其自愈控制技术的发展,大大提高了配电网的自 动化水平,在线网络重构成为可能,网络重构将作为主动配电网中降低网络损耗、提高电压 水平、以及对故障后负荷恢复供电的常见措施。因此,网络重构对量测配置的影响将不容忽 视,基于一种网络结构的量测配置,已不能达到配电网状态估计精度的要求。因而设计一种 兼顾多种网络结构的鲁棒性量测系统是有必要的。
[0003] 合理的量测配置是保证系统可观测的重要条件,系统可观测是状态估计的前提。 目前可观性分析的方法主要是针对输电网,输电网因其量测冗余度高,一般能满足可观性 要求。而配电网规模较大、网络结构变化较快,量测冗余度不足,从经济角度考虑,大范围地 进行量测配置达到传统意义下的可观性是不切实际的。为实现系统的可观性,配电网中一 般需添加误差较大的伪量测数据,这影响了状态估计的精度,而状态估计结果的准确性对 其他应用有很大的影响。因此,在保证系统可观性及经济性基础上,量测配置应尽可能提高 状态估计的精度及鲁棒性。
[0004] 近年来,配电网的量测配置研究倍受关注,但大多是针对单一网络结构。虽然国外 已有少数文献在量测配置过程中兼顾多种网络结构的估计精度,但各网络结构的影响权重 相同,并未考虑不同网络结构对量测配置的影响差异。事实上,不同网络结构在配电网运行 中所占的比重不同,比如优化重构后的网络结构运行的时间可能相对于故障重构后的网络 结构要长一些。因此,在量测配置中,各网络结构应赋予不同的影响权重,以表征这种影响 差异。除此之外,目前以量测数量为约束的量测配置方法并没有给出量测数量上限的依据, 大都只是根据已有经验确定。
[0005] 目前的量测配置方法存在以下问题:
[0006] (1)针对单一网络结构,最终的量测配置方案仅仅保证单一网络结构的状态估计 精度满足要求,若网络结构发生变化,其估计精度难以保证,量测配置的鲁棒性不足;
[0007] (2)少数针对多种网络结构的量测配置方法,各网络结构的影响权重相同,未考虑 不同网络结构对量测配置的影响差异;
[0008] (3)没有给出量测数量上限的依据。

【发明内容】

[0009] 本发明为了解决上述问题,提出了一种计及多种网络结构的配电网量测配置方 法,本方法首先提出了量测系统的饱和特性,并确定了系统的饱和量测数量,以此作为量测 数量上限约束;然后基于启发式算法确定最终量测配置方案,能够兼顾多种网络结构的估 计精度,鲁棒性提高。
[0010] 为了实现上述目的,本发明采用如下技术方案:
[0011] -种计及多种网络结构的配电网量测配置方法,包括以下步骤:
[0012] (1)预估配电网的运行状态,确定量测量与状态量的关系;
[0013] (2)根据配电网实际的运行状态,在量测配置时构建多种不同状态下的网络结构;
[0014] (3)分析配电网系统在不同网络结构下的饱和量测特性,基于启发式算法确定饱 和量测数量;
[0015] (4)在饱和量测数量为量测数量的安装数量上限的约束条件下,以多种网络结构 估计误差的加权和最小为目标,构建量测配置模型;
[0016] (5)利用层次分析法求解量测配置模型的权重系数,遍历所有可选方案,求解最优 解,输出最优量测配置方案。
[0017]所述步骤(1)中,计及的状态矢量具体包括所有节点的电压幅值及除平衡节点外 的电压相角,量测矢量类型包括实时电压幅值量测、实时功率量测、虚拟量测和伪量测。
[0018] 所述步骤(2)中,具体包括正常运行、弱环网、优化重构后的网络结构以及故障重 构后的网络结构。
[0019] 所述步骤(2)中,弱环网网络结构为闭合某一联络开关后形成的网络。
[0020] 所述步骤(2)中,以网损最小、电压偏移及负荷均衡为目标,通过遗传算法求得优 化重构后的网络结构。
[0021] 所述步骤(2)中,故障重构后的网络结构由选取某一支路故障,通过遗传算法求 得。
[0022]所述步骤(3)中,具体步骤包括:
[0023] (3-1)初始化系统数据,计算在缺省的量测配置时系统状态估计总方差;
[0024] (3-2)每次增加一个量测,遍历所有可选方案,使状态估计总方差最小的方案即为 新增量测的位置和类型;
[0025] (3-3)继续增加量测,直到量测数量使系统状态估计总方差变化量与系统状态估 计总方差的比值小于设定值。
[0026] 所述步骤(4)中,各网络结构的估计误差的加权和为每个网络结构对量测配置的 影响权重与相应网络结构的系统状态估计总误差乘积之和。
[0027] 所述步骤(4)中,采用M。次蒙特卡洛法计算的均值作为各状态变量的估计值以表 征量测系统的不确定性。
[0028] 所述步骤(5)中,确定各网络结构对量测配置的影响权重的具体步骤包括:
[0029] (5-1)分析各因素之间的关系,建立递阶层次结构;
[0030] (5-2)将各因素对同一准则或目标的重要性进行两两比较,依据判断矩阵标度表, 对各准则进行评分,形成判断矩阵;
[0031] (5-3)由判断矩阵计算被比较元素对上一层准则的相对权重,并进行一致性检验。 [0032]进一步的,已知判断矩阵,用特征值法求得各权重的值,对判断矩阵进行一致性检 验,检测判断矩阵是否可接受,否则需对判断该矩阵做适当调整。
[0033]本发明的有益效果为:
[0034] (1)本发明分析了配电网中量测系统的饱和特性,确定了饱和量测数量;
[0035] (2)本发明以饱和量测数量作为量测设备的安装数量上限,对配电系统量测配置 的建设投资给出指导性方案,避免了盲目安装大量量测装置而估计精度又没有明显提高所 造成的经济浪费,具有实际意义;
[0036] (3)本发明的量测配置模型考虑了多种网络结构,最终量测配置能够兼顾多种网 络结构的估计精度,鲁棒性提高;
[0037] (4)本发明以权重系数的方式分析了不同网络结构对量测配置的影响差异,量测 配置结果更切实际。
【附图说明】
[0038] 图1本发明提供的设计方案流程图;
[0039] 图2(a)和图2(b)分别为本发明提供的层次分析法的结构图及确定权重流程图;
[0040] 图3本发明提供的基于启发式算法确定最终量测配置方案;
[0041 ]图4本发明提供的IEEE33节点系统接线图;
[0042]图5本发明提供的不同网络结构下,基于启发式算法确定的IEEE33节点系统的饱 和量测特性。
【具体实施方式】:
[0043]下面结合附图与实施例对本发明作进一步说明。
[0044] 如图1所示,计及多种网络结构的配电网量测配置方法,包括以下步骤:
[0045] (1)状态估计方法选用加权最小二乘法;
[0046] (2)根据配电网实际的运行状态,量测配置考虑的网络结构包括正常运行网络结 构、弱环网、优化重构后的网络结构以及故障重构后的网络结构四类;
[0047] (3)分析系统在不同网络结构下的饱和量测特性,基于启发式算法确定系统的饱 和量测数量,并以此作为量测配置模型中量测数量的安装数量上限;
[0048] (4)建立量测配置的数学模型,以多种网络结构估计误差的加权和作为衡量系统 估计精度的标准,并以此为目标。约束条件为量测设备的安装数量上限约束。
[0049] (5)前述步骤(4)的权重系数由层次分析法(Analytic Hierarchy Process,AHP) 求得;
[0050] (6)基于启发式算法确定最终的量测配置方案,即遍历所有可选方案,新增量测位 置及类型的选取使所述步骤(4)的目标函数最小。重复该过程,直到达到量测设备的安装数 量上限;
[0051] (7)输出最优量测配置方案。
[0052] 前述步骤(1)状态估计中量测量与状态量之间的关系为:
[0053] z = f (x)+e
[0054] 式中:x为状态矢量,jsf [ n,巧,v-,』f,包含所有节点的电压幅值 及除平衡节点外的电压相角;z为量测矢量,量测类型包括实时电压幅值量测,实时功率量 测、虚拟量测(零注入节点)和伪量测(节点负荷);f(x)为量测方程;e为量测误差矢量,e~N (0,R),貧[CT A CT ?,…,GT.J],为量测误差的协方差矩阵。
[0055] 另有:
[0056] H(x) = <^-~)- dx
[0057] C=(htR-咕)一1
[0058]式中:H为系统的雅克比矩阵,C为系统状态估计误差的协方差矩阵,是对称非负定 矩阵;Ckj或Cjk是状态变量Xk和Xj的协方差;Cii为C中第i个对角元素,对应状态变量Xi的方 差。
[0059]随着Cll的增大,状态变量^的误差随之增大,反之亦然。定义系统状态估计总方差 为:
[0061 ]式中:N为系统中的节点总数,(2N-1)为系统中的状态变量总数。
[0062] 前述步骤(2)弱环网、优化重构后的网络结构以及故障重构后的网络结构:
[0063] 1)闭合某一联络开关,形成弱环网;
[0064] 2)以网损最小、电压偏移及负荷均衡为目标,通过遗传算法求得优化重构后的网 络结构;
[0065]优化重构的目标函数为:
[0067] 式中:Nbr是系统的总支路数;Ii是第i条支路的电流是所有支路电流的平均 值;VjPVNl*别为第i个节点电压的真值和额定值; ri是第i条支路的电阻;PlcissQ是正常网络 结构下的有功网损;m〇ad,u?it age和Uloss是三个子目标的权重。所有的值都是标么值的形式。
[0068] 3)选取某一支路故障,通过遗传算法求得故障重构后的网络结构。
[0069] 前述步骤(3)基于系统状态估计总方差与量测数量的关系,推导了量测系统的饱 和特性,并确定了饱和量测数量,具体推导过程如下:
[0070] 假定当前有m个量测量,对应的系统雅可比矩阵、量测误差的协方差矩阵及状态估 计的协方差矩阵分别为化、心和(^。新增一个量测量,对应的系统雅可比矩阵、量测误差的协 方差矩阵及状态估计的协方差矩阵分别为H m+1、UPCm+1。其中,Hm+1可表示为:
L0072」式中:hm+1为雅可比矩阵中新增量测对应的行向量。
[0073] 根据Sherman-Morrison方程,Cm和Cm+i间的关系可表示为:
[0076]假定新增量测为功率量测,则其与四个状态变量相关,hm+冲有四个非零元素,分 别表示为:
[0081] 式中:fm+i(x)为第(m+1)个量测量对应的量测方程,€4、(1」、(11及(^分别为]1 111+1中第:[、 」、1及1"个元素。
[0082] 矩阵A中的对角元素即为新增一个量测量,(^中对角元素的减小量,SP:
[0084] 式中:akk是A的第k个对角元素,cki、Ckj、c kl和ckr分别为Cm中第i、j、1及r列中的第k 个元素 ,kG 1,2,3,…,2N-1 .Wm+i = l/(〇m+i2)是第(m+1)个量测量的权重.d是计算过程中的 中间变量,可表示为:
[0085] d = ^ +tXjCJI, +ajcln +arcm) nav
[0086] 式中:v={i,j,l,r} ?
[0087] 显然,akk的正负由d的符号决定。为d证明的正负,可将Cij表示为:
[0088]
[0089]式中:-1彡Pij彡1,是第i和第j个状态变量的相关系数。
[0090] 取|P |=1,则d可等价为:
[0091 ] " ={(U: ypl、± (U.. ) ±[a, j,',, ) ±(a,.:sp^ 0
[0092]由此可得,d》0。因此,akk>0,贝撕增一个量测量,系统状态估计总方差变化量总大 于0,即:
[0094] 式中:Csys,jPCsys,m+1分别为量测数量为m和m+1时系统状态估计总方差。
[0095] -般的,若新增量测为电压量测,hm+1中只有一个非零元素,ACsys, m>0仍然成立。 [0096] 另外,A Csys具有子模块饱和特性(Submodularity),即当已有量测很少时,新增一 个量测A Csys的下降量要大于已有较多量测情况下新增一个量测时A Csys的下降量。也即:
[0097] A Csys(Sl U S2)+ A Csys(Sl n S2)^ A Csys(Sl)+ A Csys(S2)
[0098] 式中:&、S2 e 0,0是所有可能的量测位置及类型的备选方案。
[0099] 对于及£爲£ ft s e 0/B2,A Csys满足:
[0100] A Csys(Bl U s)- A Csys(Bl)^ A Csys(B2 u s)- A Csys(B2)
[0101] 因此,随着量测数量的增加 ,a csys逐渐减小。定义当满足下式时,对应的量测数量 为饱和量测数量Ns。
[0103] 式中:Csys,o是基准值,为在缺省的量测配置时(变电站处的一个功率量测和一个电 压量测)系统状态估计总方差;A Csys为新增一个量测时,Csys的下降量;e为与状态估计精度 有关的变量。
[0104] 在此,对于某一网络结构下,可基于启发式算法确定饱和量测数量,具体步骤如 下:
[0105] 1)初始化系统数据;
[0106] 2)计算 Csys,o;
[0107] 3)每次增加一个量测,遍历所有可选方案,使Csys最小方案即为新增量测的位置和 类型;
[0108] 4)增加量测,直到量测数量使
满足。
[0109] 不同网络结构的饱和量测数量均可按上述方案求得。
[0110] 前述步骤(4)的目标函数F及约束条件如下: min ^^ = > p: K f Com] com sJ. M < N&
[0112] 式中:E_为系统综合估计误差別为第j种网络结构对量测配置的影响权重,0<Pj 彡1,且,=1 N s t r为所考虑的网络结构总数,E s, j为第j个网络结构的系统状态估计总误 , 差。
[0113] 其中,各网络结构的系统状态估计精度由各状态变量的相对估计误差之和表示:
[0115] 式中:ES为系统状态估计总误差;2N-1为状态变量的个数;Xi,true3和xipst分别为第i 个状态变量的真值和估计值,其中各状态变量真值由潮流计算求得。
[0116] 为表征量测系统的不确定性,采用M。次蒙特卡洛法计算的均值作为各状态变量的 估计值,即:
[0118] 式中:Mc = 400,为蒙特卡洛的运行次数;Xi,j为第j次蒙特卡洛第i个状态变量的估 计值。
[0119] 前述步骤(5)由AHP确定各网络结构对量测配置的影响权重pi的步骤如下:
[0120] 1)分析系统中各因素之间的关系,建立递阶层次结构,如图2(a)所示;
[0121] 2)将各因素对同一准则或目标的重要性进行两两比较,依据判断矩阵标度表,对 各准则进行评分,形成判断矩阵;
[0122] 3)由判断矩阵计算被比较元素对上一层准则的相对权重,并进行一致性检验;
[0123] 已知判断矩阵,可用特征值法求得各权重的值。对判断矩阵进行一致性检验,检测 判断矩阵是否可接受,否则需对判断该矩阵做适当调整。
[0124] 前述步骤(6)基于启发式算法求解量测配置的过程如下:
[0125] 1)计算量测在所有可选方案的E。?并保存;
[0126] 2)比较上述所有E_,使E_最小的情况即为新增量测的安装类型和位置。
[0127] 3)重复步骤1) 一 3),直到达到量测设备的安装数量上限。
[0128] 由图5可以看出,对于同一系统,不同网络结构下的量测饱和特性基本一致。随着 量测数量的增加,系统的状态估计总误差随之呈非线性下降,但其下降量越来越小。当量测 数量增加到某一值,量测系统达到饱和。虽然增加量测,系统状态估计总方差会略有减小, 但下降幅度很小,可以忽略。因此,以饱和量测数量作为量测配置的安装数量上限可以保证 配电网的状态估计精度。
[0129] 当e = 0.01时,饱和量测数量为9。以饱和量测数量为量测设备的安装数量上限,本 发明确定的最终量测配置方案包含5个功率量测和4电压幅值量测。功率量测分别安装在支 路2、5、10、18和21,电压幅值量测分别安装在节点9、19、26和29。基于本发明确定的量测配 置方案,系统综合估计误差为E_=0.004348,估计精度较高。因此,本发明提出的计及多种 网络结构的量测配置方法,能够兼顾多种网络结构运行情况下的状态估计精度,同时又能 保证正常运行网络结构的状态估计精度,具有鲁棒性。
[0130] 上述虽然结合附图对本发明的【具体实施方式】进行了描述,但并非对本发明保护范 围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不 需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
【主权项】
1. 一种计及多种网络结构的配电网量测配置方法,其特征是:包括以下步骤: (1) 预估配电网的运行状态,确定量测量与状态量的关系; (2) 根据配电网实际的运行状态,在量测配置时构建多种不同状态下的网络结构; (3) 分析配电网系统在不同网络结构下的饱和量测特性,基于启发式算法确定饱和量 测数量; (4) 在饱和量测数量为量测数量的安装数量上限的约束条件下,以多种网络结构估计 误差的加权和最小为目标,构建量测配置模型; (5) 利用层次分析法求解量测配置模型的权重系数,遍历所有可选方案,求解最优解, 输出最优量测配置方案。2. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(1)中,计及的状态矢量具体包括所有节点的电压幅值及除平衡节点外的电压相角,量 测矢量类型包括实时电压幅值量测、实时功率量测、虚拟量测和伪量测。3. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(2)中,具体包括正常运行、弱环网、优化重构后的网络结构以及故障重构后的网络结 构。4. 如权利要求3所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(2)中,弱环网网络结构为闭合某一联络开关后形成的网络。5. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(2)中,以网损最小、电压偏移及负荷均衡为目标,通过遗传算法求得优化重构后的网 络结构; 所述步骤(2)中,故障重构后的网络结构由选取某一支路故障,通过遗传算法求得。6. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(3)中,具体步骤包括: (3-1)初始化系统数据,计算在缺省的量测配置时系统状态估计总方差; (3-2)每次增加一个量测,遍历所有可选方案,使状态估计总方差最小的方案即为新增 量测的位置和类型; (3-3)继续增加量测,直到量测数量使系统状态估计总方差变化量与系统状态估计总 方差的比值小于设定值。7. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(4)中,各网络结构的估计误差的加权和为每个网络结构对量测配置的影响权重与相 应网络结构的系统状态估计总误差乘积之和。8. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(4)中,采用M。次蒙特卡洛法计算的均值作为各状态变量的估计值以表征量测系统的 不确定性。9. 如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所述 步骤(5)中,确定各网络结构对量测配置的影响权重的具体步骤包括: (5-1)分析各因素之间的关系,建立递阶层次结构; (5-2)将各因素对同一准则或目标的重要性进行两两比较,依据判断矩阵标度表,对各 准则进行评分,形成判断矩阵; (5-3)由判断矩阵计算被比较元素对上一层准则的相对权重,并进行一致性检验。10.如权利要求1所述的一种计及多种网络结构的配电网量测配置方法,其特征是:所 述步骤(5)中,已知判断矩阵,用特征值法求得各权重的值,对判断矩阵进行一致性检验,检 测判断矩阵是否可接受,否则需对判断该矩阵做适当调整。
【文档编号】H02J3/00GK105958487SQ201610431772
【公开日】2016年9月21日
【申请日】2016年6月15日
【发明人】张文, 王红
【申请人】山东大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1