发明名称 基于道路矩阵的配电网三相解耦潮流计算方法
摘要 本发明公开了一种基于道路矩阵的配电网三相解耦潮流计算方法,首先,采用对称分量法对三相不平衡配电网进行序分量解耦,得到零序、正序和负序配电序网络,接着,采用基于道路矩阵的回路分析法进行一相序分量网络潮流计算,求出三序网的潮流,其次,通过对称分量法的逆变换原理把序网潮流变换到相分量模型中,求出三相潮流,通过本发明把三相不平衡的配电网系统解耦为零序、正序和负序网络,避免了三相潮流计算过程中大矩阵运算,减少了计算量,提高了计算效率。整个发明的计算过程清晰,编程简单,计算速度快。最后,通过6母线测试算例验证了本发明的正确性和良好的收敛性,同时,具有很好的通用性和实用性。
申请公布号 CN102842907A 申请公布日期 2012.12.26
申请号 CN201210334037.9 申请日期 2012.09.11
申请人 河海大学 发明人 杨雄;卫志农;孙国强;孙永辉;韦延方;袁阳;陆子刚
分类号 H02J3/00(2006.01)I;G06F19/00(2011.01)I 主分类号 H02J3/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 1.一种基于道路矩阵的配电网三相解耦潮流计算方法,其特征在于,包括以下步骤:1)设首节点是电源且作为参考节点,电源节点三相电压相量矩阵为<img file="FDA00002125894600011.GIF" wi="131" he="59" />各节点三相电压相量矩阵为<img file="FDA00002125894600012.GIF" wi="131" he="59" />在配电系统序网络中,得出电源节点的三序电压矩阵为<img file="FDA00002125894600013.GIF" wi="327" he="59" />各节点三序电压矩阵为<img file="FDA00002125894600014.GIF" wi="126" he="59" />其中,令a=e<sup>j2π/3</sup>,<maths num="0001"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>a</mi></mtd><mtd><msup><mi>a</mi><mn>2</mn></msup></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><msup><mi>a</mi><mn>2</mn></msup></mtd><mtd><mi>a</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msup><mi>A</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><msup><mi>a</mi><mn>2</mn></msup></mtd><mtd><mi>a</mi></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>a</mi></mtd><mtd><msup><mi>a</mi><mn>2</mn></msup></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>n为独立节点的个数,则独立支路条数为b=n;即针对具有N个节点的三相辐射状配电网,假设首节点是电源且作为参考节点,则独立节点个数为n=N-1,独立支路条数b=n;2)计算各序网络参数<img file="FDA00002125894600017.GIF" wi="154" he="52" />为基于支路i的序阻抗<img file="FDA00002125894600018.GIF" wi="56" he="52" />形成的对角阵,其中,上标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;支路i的三相阻抗为<img file="FDA00002125894600019.GIF" wi="134" he="59" />则有<img file="FDA000021258946000110.GIF" wi="388" he="59" />其中,<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mi>bi</mi><mn>0,1,2</mn></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>Z</mi><mi>bi</mi><mn>0</mn></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msubsup><mi>Z</mi><mi>bi</mi><mn>1</mn></msubsup></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>Z</mi><mi>bi</mi><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mi>bi</mi><mrow><mi>a</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Z</mi><mi>iaa</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>iab</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>iac</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>iba</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>ibb</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>ibc</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>ica</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>icb</mi></msub></mtd><mtd><msub><mi>Z</mi><mi>icc</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>3)在解耦的各序网模型电路中,计算各序网络的道路矩阵T<sub>s</sub>;另对于注入序电流为零的节点,在各序网的道路矩阵T<sub>s</sub>中把该节点所对应行删去后形成新矩阵为T<sub>sg</sub>。其中,下标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络;4)计算各序网中阻抗灵敏性矩阵<img file="FDA000021258946000113.GIF" wi="300" he="63" />5)给配电网各节点三相电压赋初始值<img file="FDA000021258946000114.GIF" wi="360" he="60" />其中E<sub>n</sub>=[E,E,…,E]<sup>T</sup>,共n个E,E为3×3单位矩阵;6)计算d次迭代时节点i注入的各相电流<img file="FDA000021258946000115.GIF" wi="662" he="77" />其中<img file="FDA000021258946000116.GIF" wi="52" he="60" />是节点i各相注入功率,<img file="FDA000021258946000117.GIF" wi="53" he="52" />是节点i各相并联导纳之和,p=a,b,c,i=1,2,…,m;m为节点注入序电流不为零的节点个数,d为迭代次数变量;7)计算d次迭代时节点i注入的各序电流<img file="FDA00002125894600021.GIF" wi="314" he="72" />i=1,2,…,m;8)计算d次迭代时的<img file="FDA00002125894600022.GIF" wi="378" he="78" />其中,<img file="FDA00002125894600023.GIF" wi="57" he="66" />为d次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵,m为节点注入序电流不为零的节点个数,上标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;9)计算d次迭代时的<img file="FDA00002125894600024.GIF" wi="454" he="76" />其中,1<sub>n</sub>=[1,1,…,1]<sup>T</sup>,共n个1;s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;10)基于逆变换计算d次迭代时节点i三相电压相量<img file="FDA00002125894600025.GIF" wi="377" he="72" />i=1,2,…,n;11)判断<img file="FDA00002125894600026.GIF" wi="106" he="66" />和<img file="FDA00002125894600027.GIF" wi="106" he="71" />幅值之差是否满足收敛精度要求,满足结束迭代;不满足转步骤6)。
地址 210098 江苏省南京市鼓楼区西康路1号