发明名称 一种适用于态势分析的配电网均匀接线图的生成方法
摘要 本发明涉及一种适用于态势分析的配电网均匀接线图的生成方法。本发明首先获取包含变电站和负荷的位置及线路的原始信息;再将变电站和负荷看为节点,并对每个节点遍历与其相连节点,作引力合力计算;再根据点边和边边规避得到斥力矢量和;得到合力后将其代入有牵引力的运动方程以移动节点;迭代至节点的最终位移量小于定值后终止本次拓扑中的节点移动;最后导出节点最终位置,写入SVG,完成该配电网均匀接线图的生成。本发明可满足系统对SVG语法格式的要求;成图效率高;对其进行电压等高线渲染后,图形的渲染效果层次分明,能够大大提升配电网运行人员对电网态势的把控。
申请公布号 CN103605838B 申请公布日期 2016.09.28
申请号 CN201310537677.4 申请日期 2013.11.01
申请人 杭州电子科技大学 发明人 章坚民;陈昊
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 杭州君度专利代理事务所(特殊普通合伙) 33240 代理人 杜军
主权项 一种适用于态势分析的配电网均匀接线图的生成方法,其特征在于,该方法包括以下步骤:步骤1.原始信息获取,包含变电站和负荷的地理位置信息及线路构成信息,获取源:电网EMS平台系统;步骤2.将变电站和负荷当做均匀接线图中的节点,并对其中的节点p<sub>i</sub>根据其线路连接情况遍历与其相连节点,从而对其做引力的合力<img file="FDA0000974053260000011.GIF" wi="70" he="87" />计算,得到<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mrow><mi>a</mi><mi>i</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msub><msup><mover><mi>d</mi><mo>&RightArrow;</mo></mover><mn>2</mn></msup><mrow><msub><mi>p</mi><mi>i</mi></msub><msub><mi>p</mi><mi>j</mi></msub></mrow></msub></mrow><mi>k</mi></mfrac></mrow>]]></math><img file="FDA0000974053260000012.GIF" wi="348" he="175" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><mi>C</mi><mo>&times;</mo><msqrt><mrow><mi>W</mi><mo>&times;</mo><mi>H</mi><mo>/</mo><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msqrt></mrow>]]></math><img file="FDA0000974053260000013.GIF" wi="605" he="91" /></maths>式中f<sub>ai</sub>为该点所受引力的矢量和,n和m分别为对点p<sub>i</sub>有引力和斥力的点的数量和,d为点p<sub>j</sub>与点p<sub>i</sub>间的距离向量,C为常量,W和H分别表示布局区域的宽度和长度;步骤3.根据点边规避对节点p<sub>i</sub>进行点边斥力<img file="FDA0000974053260000014.GIF" wi="194" he="77" />计算如下式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>R</mi><mo>&times;</mo><mo>|</mo><mi>d</mi><mo>-</mo><mi>d</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>&lt;</mo><mi>d</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>d</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000974053260000015.GIF" wi="1007" he="170" /></maths>式中,R为点边斥力常量,d(p<sub>i</sub>,e<sub>j</sub>)为节点和边之间的垂直距离;d为点和边间的最小距离;步骤4.根据边边规避对节点p<sub>i</sub>进行点边斥力<img file="FDA0000974053260000016.GIF" wi="198" he="78" />的矢量和计算如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mrow><mi>r</mi><mi>i</mi></mrow></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>m</mi></munderover><mfrac><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>&times;</mo><msub><mi>r</mi><mn>1</mn></msub></mrow><msub><mover><mi>d</mi><mo>&RightArrow;</mo></mover><mrow><msub><mi>p</mi><mi>i</mi></msub><msub><mi>p</mi><mi>j</mi></msub></mrow></msub></mfrac><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>&times;</mo><msub><mi>r</mi><mn>2</mn></msub></mrow><msub><mover><mi>d</mi><mo>&RightArrow;</mo></mover><mrow><msub><mi>p</mi><mi>i</mi></msub><msub><mi>p</mi><mi>j</mi></msub></mrow></msub></mfrac><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>q</mi></munderover><mfrac><msup><mi>k</mi><mn>2</mn></msup><msub><mover><mi>d</mi><mo>&RightArrow;</mo></mover><mrow><msub><mi>p</mi><mi>i</mi></msub><msub><mi>p</mi><mi>j</mi></msub></mrow></msub></mfrac><mo>+</mo><msub><mover><mi>f</mi><mo>&RightArrow;</mo></mover><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><msub><mi>e</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000974053260000017.GIF" wi="1349" he="183" /></maths>式中,m为与当前节点不同馈线上的节点数;n为当前节点上一级别支路上的节点数;q为剩余的节点数;r<sub>1</sub>为与当前节点不同馈线上的节点对其的斥力权重;r<sub>2</sub>为当前节点上一级别支路上的节点对其的斥力权重;一般情况下,1&lt;r<sub>1</sub>&lt;r<sub>2</sub>&lt;1.3,<img file="FDA0000974053260000018.GIF" wi="199" he="76" />为该点所受的点边斥力;步骤5.节点的引力和斥力相合成合力后将其带入引入牵引力参数的节点运动方程,从而移动节点;其带牵引力参数的方程为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>x</mi><mrow><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>x</mi><mi>l</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>f</mi><mi>x</mi><mi>l</mi></msubsup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>g</mi><mo>&times;</mo><msup><mi>x</mi><mi>l</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>T</mi><mi>l</mi></msub></mrow>]]></math><img file="FDA0000974053260000021.GIF" wi="867" he="71" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>y</mi><mrow><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>y</mi><mi>l</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>f</mi><mi>y</mi><mi>l</mi></msubsup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>g</mi><mo>&times;</mo><msup><mi>y</mi><mi>l</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>T</mi><mi>l</mi></msub></mrow>]]></math><img file="FDA0000974053260000022.GIF" wi="832" he="71" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>T</mi><mi>l</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>j</mi><mo>&lt;</mo><mo>=</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>t</mi><mo>&times;</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mi>n</mi><mo>)</mo></mrow><mo>/</mo><mi>m</mi></mrow></mtd><mtd><mrow><mi>j</mi><mo>&gt;</mo><mi>n</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000974053260000023.GIF" wi="470" he="134" /></maths>式中l为节点pi的迭代次数,g为重力常量,f为该节点所受合力,T<sub>l</sub>为牵引力位移,t为牵引力常量,j为p<sub>i</sub>点所在该条线路中的序号,m和n分别为p<sub>i</sub>点在该条线路的负荷数和变电站线路平均负荷数;步骤6.对均匀接线图中的剩余的节点重复步骤2~步骤5,若节点遍历完则从起始节点重新开始节点位置的移动,直到节点的最终位移量小于当前拓扑分布总长度的1/20,说明各个节点所受合力趋向于零,则认为系统达到平衡,即在可以接受的范围之内,从而终止本次拓扑中节点的移动;步骤7.导出节点最终位置,写入SVG;在循环迭代结束之后,将得到各个节点和线路的最终位置,在将其写入SVG后,最终完成该配电网均匀接线图的生成。
地址 310027 浙江省杭州市江干区2号大街1号