发明名称 含多变压器支路的配电网三相解耦潮流计算方法
摘要 本发明公开了一种含多变压器支路的配电网三相解耦潮流计算方法,根据配电网的三相线路参数基本对称,三相负荷不平衡和配电网树状结构等特点,首先,运用对称分量法对配电网三相不平衡网络进行序分量解耦。接着,在配电三序网络模型中,运用相位变换技术将变压器支路变换成普通支路。最后,采用基于道路矩阵的回路分析法计算配电系统序网潮流,并运用逆变换原理将序分量潮流变换成相分量潮流,从而实现含多变压器支路的配电网三相不平衡系统的三相解耦潮流计算,减少了计算量。整个发明的计算过程清晰,编程简单,计算速度快。最后,通过34母线测试算例验证了本发明的正确性和良好的收敛性,同时,具有很好的通用性和实用性。
申请公布号 CN102842908B 申请公布日期 2014.11.12
申请号 CN201210335586.8 申请日期 2012.09.11
申请人 河海大学 发明人 杨雄;卫志农;孙国强;孙永辉;韦延方;袁阳;陆子刚
分类号 H02J3/00(2006.01)I;G06F19/00(2011.01)I 主分类号 H02J3/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种含多变压器支路的配电网三相解耦潮流计算方法,其特征在于,包括以下步骤:1)设首节点是电源且作为参考节点,电源节点三相电压相量矩阵为<img file="FDA0000556783130000011.GIF" wi="175" he="80" />各节点三相电压相量矩阵为<img file="FDA0000556783130000012.GIF" wi="143" he="80" />在配电系统序网络中,得出电源节点的三序电压矩阵为<img file="FDA0000556783130000013.GIF" wi="349" he="80" />各节点三序电压矩阵为<img file="FDA0000556783130000014.GIF" wi="145" he="80" />其中,令a=e<sup>j2π/3</sup>,<maths num="0001" id="cmaths0001"><math><![CDATA[<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><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><img file="FDA0000556783130000015.GIF" wi="846" he="234" /></maths>n为独立节点的个数,则独立支路条数为b=n;即针对具有N个节点的三相辐射状配电网,假设首节点是电源且作为参考节点,则独立节点个数为n=N‑1,独立支路条数b=n;2)根据配电网中变压器的个数K,相应的把配电网划分为K块区域,依据参考节点,各变压器的接线方式依次计算出各块区域的相位变换矩阵Θ<sub>k</sub>,并计算出各块区域的解耦变换矩阵A<sub>k</sub>=Θ<sub>k</sub>A;其中,k表示配电网中第k块区域,k∈{1,2,…,K};Θ为相位变换矩阵,是一个3×3对角阵,且<img file="FDA0000556783130000016.GIF" wi="502" he="77" />θ<sub>0</sub>、θ<sub>1</sub>和θ<sub>2</sub>分别为三序网络系统中的零序、正序和负序相移量,下标“0”、“1”和“2”分别代表三序网中的零序、正序和负序;同时求出<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>A</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msup><mi>A</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mi>&Theta;</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000556783130000017.GIF" wi="292" he="78" /></maths>3)计算各序网络参数<img file="FDA0000556783130000018.GIF" wi="95" he="71" /><img file="FDA0000556783130000019.GIF" wi="69" he="71" />为基于支路i的序阻抗<img file="FDA00005567831300000110.GIF" wi="70" he="71" />形成的对角阵,其中,上标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;支路i的三相阻抗为<img file="FDA00005567831300000111.GIF" wi="141" he="78" />假设其属于第k块区域,则有<img file="FDA00005567831300000112.GIF" wi="410" he="78" />其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<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><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><img file="FDA00005567831300000113.GIF" wi="1117" he="249" /></maths>其中A<sub>k</sub>为支路i所在区域的解耦变换矩阵;4)在解耦的各序网模型电路中,计算各序网络的道路矩阵T<sub>s</sub>;另对于注入序电流为零的节点,在各序网的道路矩阵T<sub>s</sub>中把该节点所对应行删去后形成新矩阵为T<sub>sg</sub>;其中,下标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络;5)计算各序网中阻抗灵敏性矩阵<img file="FDA0000556783130000021.GIF" wi="316" he="82" />6)给配电网各节点三相电压赋初始值<img file="FDA0000556783130000022.GIF" wi="372" he="80" />其中E<sub>n</sub>=[E,E,…,E]<sup>T</sup>,共n个E,E为3×3单位矩阵;7)计算d次迭代时节点i注入的各相电流<img file="FDA0000556783130000023.GIF" wi="677" he="97" />其中<img file="FDA0000556783130000024.GIF" wi="68" he="81" />是节点i各相注入功率,<img file="FDA0000556783130000025.GIF" wi="139" he="95" />是第d‑1次迭代时求出的节点i各相电压相量,Y<sub>i</sub><sup>p</sup>是节点i各相并联导纳之和,p=a,b,c,i=1,2,…,m;m为节点注入序电流不为零的节点个数,d为迭代次数变量;8)计算d次迭代时节点i注入的各序电流<img file="FDA0000556783130000026.GIF" wi="349" he="94" />i=1,2,…,m;其中,A<sub>k</sub>为节点i所在区域的解耦变换矩阵;9)计算d次迭代时的<img file="FDA0000556783130000027.GIF" wi="391" he="97" />其中,<img file="FDA0000556783130000028.GIF" wi="84" he="100" />为d次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵,m为节点注入序电流不为零的节点个数,上标s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;10)计算d次迭代时的<img file="FDA0000556783130000029.GIF" wi="492" he="97" />其中,1<sub>n</sub>=[1,1,…,1]<sup>T</sup>,共n个1;s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型;11)基于逆变换计算d次迭代时节点i三相电压相量<img file="FDA00005567831300000210.GIF" wi="402" he="93" />i=1,2,…,n;其中,A<sub>k</sub>为节点i所在区域的解耦变换矩阵,<img file="FDA00005567831300000211.GIF" wi="128" he="93" />为第d次迭代时节点i三序电压相量矩阵;12)判断<img file="FDA00005567831300000212.GIF" wi="121" he="87" />和<img file="FDA00005567831300000213.GIF" wi="125" he="91" />幅值之差是否满足收敛精度要求;满足结束迭代;不满足转步骤7)。
地址 210098 江苏省南京市鼓楼区西康路1号