发明名称 一种环状多源配电网的电能质量监测器优化配置方法
摘要 适用于环状多源电网结构拓扑分析的新概念定义,包括正方向、双向等效线路和端点;环状多源电网结构表征矩阵的构建,包括修正邻接矩阵、系统结构矩阵、潮流方向矩阵、监测可观矩阵;权重因子的定义,包括电能质量监测器布置一般性准则归纳、不同端点属性分类及对应权重系数定义;模糊度指标定义,包括系统整体模糊度、部分最大模糊度和监测系统模糊度;电能质量监测器优化配置模型的建立;优化模型的求解,确定系统电能质量监测器最优配置方案。
申请公布号 CN105977968A 申请公布日期 2016.09.28
申请号 CN201610393772.5 申请日期 2016.06.06
申请人 浙江工业大学 发明人 翁国庆;黄飞腾;朱双双;林雨丰;方敏佳;王妍彦
分类号 H02J3/00(2006.01)I;G06F17/50(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种环状多源配电网的电能质量监测器优化配置方法,包括如下步骤:步骤1、定义适用于环状多源电网结构拓扑分析的新概念:与辐射状电网结构特点不同,环状多源电网运行方式和状态具有可变性,部分线路电流流向存在多种可能,从而影响各节点配置电能质量监测器的监测范围,要求将各种可能情况均纳入考虑;电能质量监测装置简称PQM;步骤101,正方向的规定:为匹配电能质量扰动源方向判定方法,需要建立一个固定坐标系,设定由母线指向线路方向为正方向;电能质量扰动源简称PQDS;步骤102,双向等效线路的定义:将除电源、负载引出线外的每条线路等效成两条并行双向的具有固定电流流向的线路;步骤103,端点的定义:将线路电流流入端设置为可配置PQM的端点来表示该线路,即每条线路可等效表示为两个端点;步骤2、环状多源系统结构表征矩阵的构建:基于步骤1环网拓扑分析,分别构建修正邻接矩阵、系统结构矩阵、潮流方向矩阵、监测可观矩阵描述网络结构,为PQM优化配置提供数据支撑;步骤201,修正邻接矩阵的构建:邻接矩阵A<sub>org</sub>常用于电力网络拓扑分析,其形式为<img file="FDA0001010187620000011.GIF" wi="694" he="71" />为母线节点编号,N<sub>B</sub>为系统母线节点数量;矩阵元素<img file="FDA0001010187620000012.GIF" wi="74" he="71" />的赋值原则如式(1)所示;因为网络化电能质量监测系统中负载和电源都是重要的监测对象,需要对邻接矩阵进行适当修正:用母线与等效为母线0的大地间的连接线表示负载,用母线与其自身的连接线表示电源;修正邻接矩阵<img file="FDA0001010187620000013.GIF" wi="77" he="71" />的形式为<img file="FDA0001010187620000014.GIF" wi="477" he="95" /><img file="FDA0001010187620000015.GIF" wi="38" he="71" />为增加了表示大地的母线0后的系统节点编号;元素<img file="FDA0001010187620000016.GIF" wi="70" he="78" />表示环状多源网络中每个端点的结构编号,其取值表示该端点是否存在,赋值原则如式(2)所示;<img file="FDA0001010187620000017.GIF" wi="1190" he="148" /><img file="FDA0001010187620000018.GIF" wi="1190" he="149" />步骤202,系统结构矩阵构建:将<img file="FDA0001010187620000019.GIF" wi="93" he="82" />矩阵中所有有效端点,即<img file="FDA00010101876200000110.GIF" wi="175" he="79" />的端点,按行顺序提取出,并将其依次按序标记为新的系统端点符号t<sub>j</sub>,其中j=1,2,...T为依次按序标记后t<sub>j</sub>的编号,T为有效端点总数量;利用其与母线可共同构成描述各端点与母线间固定位置关系的系统结构矩阵S<sub>SM</sub>;S<sub>SM</sub>是一个T×N<sub>B</sub>的矩阵,其元素赋值原则如式(3)所示;<img file="FDA00010101876200000111.GIF" wi="1286" he="227" />步骤203,潮流方向矩阵构建:尽管环状多源系统运行方式和状态存在多种可能,但其绝大多数时间处于正常运行状态,PQM配置需优先保证系统正常运行状态下拓扑结构的全网可观;通过比较电流方向和设定正方向,可对S<sub>SM</sub>矩阵进行修正,得到包含系统潮流信息的潮流方向矩阵C<sub>DM</sub>;C<sub>DM</sub>矩阵元素c<sub>ji</sub>的赋值原则如式(4)所示,其中j为t<sub>j</sub>的编号,i为系统中母线编号;<img file="FDA0001010187620000021.GIF" wi="1478" he="283" />步骤204,系统单距区块划分:考虑系统拓扑结构中,线路密集度对监测点配置的影响;首先基于系统正常运行状态下反映潮流方向的拓扑结构,对系统线路进行等效间距处理,选取拓扑结构中相邻母线间最短的线路作为单元距离的标尺,各母线间线路的等效间距为整数倍单元距离,并以单元间距划分区块,命名为单距区块;单距区块内单元线路数量指不包含电源和负荷引出线的线路数量;步骤205,监测可观矩阵构建:目标电网完成监测器配置,根据正常运行状态下的潮流拓扑结构,可构建监测可观矩阵M<sub>CM</sub>,用于描述各PQM与各线路间的位置关系,表征各监测器的可观范围;M<sub>CM</sub>矩阵是一个T×M的矩阵,M是监测器数量,其元素m<sub>jr</sub>赋值原则如式(5)所示,其中r为系统中PQM编号;<img file="FDA0001010187620000022.GIF" wi="1342" he="159" />步骤3、权重因子的定义:基于PQM布置一般性准则,系统中不同端点具有不同性质,存在不同的重要程度;合理定义权重因子量化描述各个端点的重要程度,使各端点的性质能在优化配置算法中得到体现;步骤301,PQM优化布置原则归纳:根据环状多源网络结构特点并结合考虑PQDS定位需求,归纳出6条PQM布置一般性准则;A1)准则1:(N‑1)准则,N为同一母线上连接的线路总数;A2)准则2:负载支路上具有更高的权重;A3)准则3:电源引出线具有最高的权重;A4)准则4:同一条线路上优先只装一个PQM;A5)准则5:除负载引出线外与该母线相连的端点小于等于两个,此类端点权重较低;相应地,与该类端点属于同一线路的端点权重较高;A6)准则6:优先满足正常运行状态下全网可观;步骤302,有效端点进行运算编号和属性分类:为提高算法效率,对所有筛选出的有效端点进行运算编号和属性分类;端点的属性分为电源端点、负荷端点、对向端点和普通端点;对向端点指除去负荷外与母线相连的端点数量小于等于2的端点,负荷端点包括负荷引出端点和与对向端点属于同一线路的端点,普通端点即不归属于其余3种类型的端点;运算编号即端点按考虑母线编号和S<sub>SM</sub>矩阵中有效端点排序进行编号;若某一端点同时包含两种属性,则对其属性值进行均一化处理;步骤303,定义权重因子:基于步骤301的一般性准则的说明,系统中不同端点具有不同性质,存在不同的重要程度;合理定义权重因子量化描述各个端点的重要程度,使各端点性质能在优化配置算法中得到体现;结合端点属性和编号,定义α<sub>ij</sub>,β<sub>ij</sub>,γ<sub>ij</sub>三个权重系数量化表征PQM配置一般性准则:β<sub>ij</sub>表征准则1;α<sub>ij</sub>表征准则2,3,5;γ<sub>ij</sub>表征准则4,6;α<sub>ij</sub>,γ<sub>ij</sub>,β<sub>ij</sub>的数学表达式分别如式(6)~(12)所示;<maths num="0001"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>T</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>C</mi><mi>i</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000031.GIF" wi="1028" he="149" /></maths><img file="FDA0001010187620000032.GIF" wi="1406" he="203" /><img file="FDA0001010187620000033.GIF" wi="1414" he="370" />其中,<maths num="0002"><math><![CDATA[<mrow><msub><mi>A</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><msub><mi>N</mi><mi>i</mi></msub></mrow><msub><mi>N</mi><mi>B</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000034.GIF" wi="1078" he="215" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>i</mi></msub><mo>=</mo><msub><mi>&omega;</mi><msub><mi>t</mi><mn>1</mn></msub></msub><mo>&CenterDot;</mo><msubsup><mi>N</mi><mn>1</mn><mi>i</mi></msubsup><mo>+</mo><msub><mi>&omega;</mi><msub><mi>t</mi><mn>2</mn></msub></msub><mo>&CenterDot;</mo><msubsup><mi>N</mi><mn>2</mn><mi>i</mi></msubsup><mo>+</mo><msub><mi>&omega;</mi><msub><mi>t</mi><mn>3</mn></msub></msub><mo>&CenterDot;</mo><msubsup><mi>N</mi><mn>3</mn><mi>i</mi></msubsup><mo>+</mo><msub><mi>&omega;</mi><msub><mi>t</mi><mn>4</mn></msub></msub><mo>&CenterDot;</mo><msubsup><mi>N</mi><mn>4</mn><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000035.GIF" wi="1485" he="79" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>D</mi><mi>k</mi></msub><mo>=</mo><mfrac><msubsup><mi>N</mi><mi>D</mi><mi>k</mi></msubsup><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>N</mi><mi>D</mi><mi>k</mi></msubsup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000036.GIF" wi="1094" he="218" /></maths><img file="FDA0001010187620000037.GIF" wi="1254" he="327" />式中:i表示母线编号;j表示有效端点序号;a<sub>ij</sub>为端点运算编号,表示i母线上连接的j序号有效端点;T<sub>ij</sub>表示运算编号为a<sub>ij</sub>的端点为不同性质时对应的属性值,如式(12)所示;<img file="FDA0001010187620000038.GIF" wi="365" he="63" />分别为电源端点、负荷端点、普通端点和对向端点对应的属性值,且其关系为<img file="FDA00010101876200000312.GIF" wi="574" he="70" /><img file="FDA00010101876200000313.GIF" wi="293" he="62" />和<img file="FDA00010101876200000310.GIF" wi="59" he="62" />分别表示母线i上电源端点、负荷端点、普通端点和对向端点数量;C<sub>i</sub>表示母线i上所有不同性质端点数量与属性值乘积的累加和,如式(10)所示;<img file="FDA00010101876200000311.GIF" wi="66" he="70" />表示单距区块k内单元线路数量;k表示单距区块编号;L表示单距区块数量;D<sub>k</sub>表示单距区块k内单元线路数量与系统所有线路数量的比值,如式(11)所示;N<sub>Bij</sub>表示端点a<sub>ij</sub>下区域包含的母线数量;N<sub>i</sub>表示母线i上端点总数量;A<sub>ij</sub>表示系统所有端点数量与系统母线数量的比值,如式(9)所示;B<sub>ij</sub>表示端点a<sub>ij</sub>下区域包含的端点数量;步骤4、模糊度指标定义:系统中布置的监测器数量不够,将导致部分端点不能被有效监测,端点识别时将出现拓扑模糊;定义不同模糊度指数形式进行系统端点模糊度描述;步骤401,规定PQM监测模糊的范畴;B1)规定1:安装有监测器的端点不模糊;B2)规定2:与同一条母线相连接的N<sub>i</sub>个端点,安装有监测器的端点数量大于或等于N<sub>i</sub>‑1,则该N<sub>i</sub>个端点不模糊;B3)规定3:除去已安装监测器的端点后仍存在两个与母线相连的端点时,若与其相连的端点不模糊,则这两个端点也不模糊;步骤402,环状多源系统监测模糊指标定义:定义两种模糊度指数,系统整体模糊度和部分最大模糊度;系统整体模糊度简称IEA,部分最大模糊度简称IMA;分别如式(13),(14)所示;<maths num="0005"><math><![CDATA[<mrow><mi>I</mi><mi>E</mi><mi>A</mi><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msubsup><mi>N</mi><mrow><mi>N</mi><mi>D</mi></mrow><mi>k</mi></msubsup></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><msub><mi>N</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000041.GIF" wi="1078" he="275" /></maths><maths num="0006"><math><![CDATA[<mrow><mi>I</mi><mi>M</mi><mi>A</mi><mo>=</mo><mfrac><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>{</mo><msub><mi>b</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>&Element;</mo><mo>&lsqb;</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>K</mi><mo>&rsqb;</mo><mo>}</mo></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>b</mi><mi>k</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000042.GIF" wi="1174" he="184" /></maths>其中:<img file="FDA0001010187620000043.GIF" wi="85" he="64" />表示k区块不能被识别的线路端点数量;b<sub>k</sub>表示区块k的整体模糊程度,<img file="FDA0001010187620000044.GIF" wi="222" he="134" />M<sub>k</sub>表示对区块k进行监测的监测器数量;K为系统区块总数;对IMA进行[0,1]归整处理,得到如式(15)所示归一化的部分最大模糊度IMA';<maths num="0007"><math><![CDATA[<mrow><msup><mi>IMA</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msubsup><mi>N</mi><mrow><mi>N</mi><mi>D</mi></mrow><mi>k</mi></msubsup><mo>&times;</mo><mi>I</mi><mi>M</mi><mi>A</mi><mo>-</mo><mn>1</mn></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msubsup><mi>N</mi><mrow><mi>N</mi><mi>D</mi></mrow><mi>k</mi></msubsup><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000045.GIF" wi="1150" he="258" /></maths>IEA和IMA'的取值范围都已归整到[0,1]区间;IEA说明当电网任意点发生电能质量扰动时,监测系统判定扰动源位置的能力;IMA'说明监测系统对各区域监测能力的均等程度;根据这两个指标,可计算配置N<sub>M</sub>个监测点的监测系统模糊度,监测系统模糊度简称IA,如式(16)所示;<maths num="0008"><math><![CDATA[<mrow><mi>I</mi><mi>A</mi><mo>=</mo><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><msub><mi>N</mi><mi>i</mi></msub><mo>-</mo><msub><mi>N</mi><mi>M</mi></msub><mo>)</mo></mrow><msqrt><mrow><mi>I</mi><mi>E</mi><mi>A</mi><mrow><mo>(</mo><mi>I</mi><mi>E</mi><mi>A</mi><mo>+</mo><msup><mi>IMA</mi><mo>&prime;</mo></msup><mo>)</mo></mrow></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000046.GIF" wi="1198" he="143" /></maths>步骤403,使用监测可观矩阵检查模糊度:监测系统的各模糊度指标可通过监测可观矩阵进行获取和校验,其主要流程为:1)构建系统监测可观矩阵M<sub>CM</sub>;2)对M<sub>CM</sub>矩阵进行分区域处理;3)根据模糊规定1,2,3从各区域中排除各PQM;4)计算各区块的b<sub>k</sub>;5)计算IEA,IMA',IA;步骤5、PQM优化配置模型的建立:基于网络特性分析,以及特征矩阵、权重系数和模糊度指标的定义,建立PQM最优布置的优化模型,目标函数和约束条件分别如式(17),(18)所示;<maths num="0009"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></munder><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>i</mi></msub></munderover><msub><mi>&alpha;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>&beta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>&gamma;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000051.GIF" wi="1270" he="168" /></maths><maths num="0010"><math><![CDATA[<mrow><msub><mi>NM</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&le;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><msub><mi>N</mi><mi>i</mi></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>B</mi></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>i</mi></msub></munderover><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&lt;</mo><msub><mi>NM</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010187620000052.GIF" wi="1285" he="151" /></maths>其中:F(x)为目标函数;NM<sub>min</sub>和NM<sub>max</sub>为监测点的最小数量和最大数量;决策变量x<sub>ij</sub>的赋值原则如式(19)所示;<img file="FDA0001010187620000053.GIF" wi="1270" he="166" />步骤6、优化模型的求解:PQM优化配置模型求解的基本思路是,首先计算权重因子,基于其建立目标函数并在给定监测器数量下进行优化,确定监测器位置;然后,基于已知的监测器位置,构建监测可观矩阵,进而基于模糊度检测流程计算IA指数;经过多次迭代,最后确定监测器的最优数量和位置,使IA值在给定约束条件下最小化。
地址 310014 浙江省杭州市下城区潮王路18号浙江工业大学科技处