一种修正系数矩阵的快速分解法潮流计算方法
一种修正系数矩阵的快速分解法潮流计算方法
【技术领域】
[0001] 本发明涉及一种电力系统的快速分解法潮流计算方法，特别是一种适合含小阻抗 支路系统的快速分解法潮流计算方法。
【背景技术】
[0002] 电力系统潮流计算是研宄电力系统稳态运行的一项基本计算，它根据给定的运行 条件和网络结构确定整个网络的运行状态。潮流计算也是其他电力系统分析的基础，如安 全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度快及内存需求少 的优点，快速分解法成为当前潮流计算的主流算法之一。
[0003] 当电力网络不存在小阻抗支路或电力网络中的小阻抗支路的电阻很小时，快速分 解法潮流计算具有良好的收敛性，但电网中存在电阻相对较大的小阻抗支路时，快速分解 法潮流计算就可能发散。电力系统小阻抗支路可以分为小阻抗线路和小阻抗变压器支路， 在数学模型上线路可以看作变比为1:1的变压器，因此下面分析时仅以小阻抗变压器支路 为例分析。变压器支路、模型见图4,变压器的非标准变比k位于节点i侧，阻抗位于标 准变比侧。当变压器的电阻ru和电抗Xij都很小时，此变压器支路即为小阻抗变压器支路。 变压器阻抗Zij=rij+jXij，导纳为：
[0005] 由于小阻抗变压器支路的阻抗很小，其阻抗上的电压降也很小，因此小阻抗变压 器支路两端的电压相量应满足：
[0007] 如果电压相量用极坐标表示，贝1J有
[0009] 式中，K^分别为节点i和节点j的电压相量；Vi、\分别为节点i和节点j的 电压幅值；0i、0 0>别为节点i和节点j的电压相角。
[0010] 如图1所示，现有快速分解法潮流计算方法，主要包括以下步骤：
[0011] A、原始数据输入和电压初始化；
[0012] 电压初始化一般采用平启动，即PV节点和平衡节点的电压幅值取给定值，PQ节点 的电压幅值取1. 〇 ;所有电压的相角都取〇. 0。这里单位采用标幺值。
[0013]B、形成节点导纳矩阵；
[0014] 设节点i和节点j原来的自电导与自电纳分别为61(|』^|、6」(|』」 (|，在它们之间增加 一条小阻抗支路后的自导纳和Y"、互导纳分别为：
[0018] C、形成修正方程的系数矩阵B'和B"并进行因子表分解；
[0019] 快速分解法修正方程为：
[0020] B，A0 =AP/V (7)
[0021] B"AV=AQ/V(8)
[0022] 式中，AP/V和AQ/V分别为有功功率和无功功率不平衡量除以电压幅值后的列 向量；AV和△ 0分别为电压幅值和电压相角修正量列向量；B'为导纳矩阵的虚部，但计 算时不计及支路电阻、对地导纳和非标准变比，矩阵中包含PQ节点和PV节点相关的行和 列；B"为导纳矩阵的虚部，仅包括与PQ节点有关的行和列。
[0023] 与小阻抗支路L相关的系数矩阵元素为：
[0030] 式中，B'n、B'jj、B'u是快速分解法系数矩阵B'的元素；B' j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素；B"wB" "、B"u是快速分解法系数矩 阵B"的元素；B" #是快速分解法系数矩阵B"中不含小阻抗支路时的元素。
[0031] D、设置迭代计数t= 0,收敛标志KP= 0,KQ= 0 ;
[0032] E、计算有功功率不平衡量AP;
[0033] 有功功率不平衡量（不包含平衡节点）为：
[0035] 式中，PiS为节点i的给定有功功率；VJP0汾别为节点i的电压幅值和电压相 角；Gu和Bu分别为导纳矩阵的电导部分和电纳部分；n为节点数。
[0036] F、判断有功功率最大不平衡量|APmax|是否小于收敛精度e;如果小于收敛精度 e，令KP= 1，转到步骤G;否则，解方程B'A0 =AP/V，修正电压相角，令KP= 0,转到步 骤H;
[0037]求解方程B'A0 =AP/V，得到A0，按下式修正电压相角：
[0038] 0 (t+D= 0W_A0 (t) (IQ)
[0039] G、判断心是否等于1 ;如果KQ= 1，转到步骤K;
[0040] H、计算无功功率不平衡量AQ;
[0041] 无功功率不平衡量（仅包含PQ节点）为：
[0043] 式中，QiS为节点i的给定无功功率；m为PQ节点数。
[0044] I、判断无功功率最大不平衡量|AQ_|是否小于收敛精度e;如果小于收敛精度 e，令&= 1 ;否则，解方程8〃八￥ =AQ/V，修正电压幅值，令KQ= 0 ;
[0045] 求解方程B〃AV=AQ/V，得到AV，按下式修正电压幅值：
[0046] y(t+1)=V(t)-AV(t) (18)
[0047] J、判断是否同时满足KP= 1和KQ= 1，如果不满足，令t=t+1，返回步骤E进行 下一次迭代；
[0048] K、计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
[0049] 步骤E和步骤F为P~0迭代，即通过AP求A0进而修正0 ;步骤H和步骤I 为Q~V迭代，即通过AQ求AV进而修正V。主流快速分解法都是按上述步骤设计算法， 即先进行P~9迭代，后进行Q~V迭代。也有文献采用先进行Q~V迭代，后进行P~ 9迭代的算法。
[0050] 对正常电力网络或含有电阻非常小的小阻抗支路的网络，快速分解法潮流计算具 有良好的收敛性，但遇到含有电阻较大的小阻抗的病态网络时，快速分解法潮流计算就可 能发散。电力系统中小阻抗支路普遍存在，潮流计算的收敛性是电力系统潮流计算这类非 线性问题的最重要指标，计算不收敛就无法得到问题的解。因此改善快速分解法潮流计算 针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。中国专利ZL201410314990. 6 提出的一种适合含小阻抗支路系统的快速分解法潮流计算方法如下：
[0051] (1)采用先进行Q~V迭代，后进行P~0迭代的方法；
[0052](2)对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：修改小阻抗变压 器支路非标准变比侧节点i对应的系数矩阵B'元素B'n、B'y标准变比侧节点j对应的 元素B'ji、B'不变。
[0057] 式中，8'"、8、、8\、8'^是快速分解法系数矩阵"的元素出'i(l、B' j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素。
[0058] 中国专利ZL201410314990. 6所提出方法有效提高了含有小阻抗支路电力系统潮 流计算的收敛性，但当小阻抗支路的电阻较大（rijAXij)时，该方法仍可能发散。
【发明内容】
[0059] 为解决现有技术存在的上述问题，本发明要提出一种修正系数矩阵的快速分解法 潮流计算方法，以改善其分析含有电阻较大的小阻抗支路电力系统的潮流计算收敛性。
[0060] 为了实现上述目的，本发明提出了适合含电阻较大的小阻抗支路系统的快速分解 法潮流计算方法来改善潮流计算收敛性。本发明的技术方案如下：一种修正系数矩阵的快 速分解法潮流计算方法，采用先进行Q~V迭代，后进行P~0迭代的方法，并对P~0 迭代修正方程的系数矩阵B'的部分元素进行修改。方案包括以下步骤：
[0061] A、原始数据输入和电压初始化；
[0062] B、形成节点导纳矩阵；
[0063] C、形成修正方程的系数矩阵B'和B"并进行因子表分解；
[0064] 对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：形成系数矩阵B'时， 给小阻抗变压器支路的电抗乘一个修正因子a，正常支路的电抗不乘这个修正因子。小阻 抗支路相关的B'元素如下：
[0069] 式中，8'"、8、、8\、8'^是快速分解法系数矩阵"的元素出'j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素；k为变压器非标准变比；a为系数矩阵 B'的修正因子。
[0070] 形成系数矩阵B'的方法，包括以下步骤：
[0071]C1、读入支路数据，设置小阻抗阈值电阻1^和电抗xmin、设置修正因子a。
[0072] C2、设置支路计数初值m= 1。
[0073] C3、取支路m的首末节点号i和j、电阻rij、电抗Xij、变比k。
[0074] C4、判断支路的电阻和电抗是否同时满足和x ，如果不满足转至步骤 C7〇
[0075]C5、令叉。=ax。。
[0076] C6、令B'B' ^-1/(10^)、B'l/(kx。），转至步骤C8。
[0077] C7、令B，B，io-l/Xj、B，。.= 1/xu。
[0078] C8、令B'B'jo-l/Xij、B'1/xij〇
[0079] C9、令m=m+l。
[0080] C10、判断m是否大于支路数1，如果m不大于1，则返回到步骤C3 ;否则，形成系数 矩阵B"并进行因子表分解。
[0081]D、设置迭代计数t= 0,收敛标志KP= 0,KQ= 0;
[0082]E、计算无功功率不平衡量AQ;
[ 0083] F、判断无功功率最大不平衡量|AQ_|是否小于收敛精度e;如果小于收敛精度 e，令&= 1，转到步骤G;否则，解方程^八乂 =AQ/V，修正电压幅值，令KQ= 0,转到步骤 H；
[0084] G、判断KP是否等于1 ;如果KP= 1，转到步骤K;
[0085]H、计算有功功率不平衡量AP;
[0086] I、判断有功功率最大不平衡量|APmax|是否小于收敛精度e;如果小于收敛精度 e，令KP= 1 ;否则，解方程B'A0 =AP/V，修正电压相角，令KP= 0 ;
[0087] J、判断是否同时满足KP= 1和KQ= 1，如果不满足，贝令t=t+1，返回步骤E进 行下一次迭代；
[0088] K、计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
[0089] 与现有技术相比，本发明具有以下有益效果：
[0090] 1、本发明通过修改快速分解法修正方程系数矩阵B'中与小阻抗相关的部分元 素，改善了快速分解法潮流计算在分析含有电阻较大的小阻抗支路系统时的收敛性问题。 采用常规快速分解法潮流计算或现有专利算法不收敛时，本发明能够可靠收敛。
[0091] 2、由于本发明不仅能有效解决常规快速分解法潮流算法分析含有电阻较大的小 阻抗支路系统的收敛性问题，同时也能对正常系统进行潮流计算，因此没有不良影响。
【附图说明】
[0092] 本发明共有附图5张。其中：
[0093] 图1是现有快速分解法潮流计算的流程图。
[0094] 图2是本发明快速分解法潮流计算的流程图。
[0095]图3是本发明形成系数矩阵B'的流程图。
[0096] 图4是电力系统变压器模型示意图。
[0097] 图5是IEEE14节点电力系统算例的接线图。
【具体实施方式】
[0098] 下面结合附图对本发明进行进一步地说明，按照图2所示流程对IEEE14节点系统 (电气和电子工程师协会14节点系统测试数据）和一个445节点实际系统两个算例进行了 计算，作为对比同时采用常规快速分解法和已有专利算法对该算例进行了计算。潮流计算 的收敛精度为0.00001。
[0099] 两种常规快速分解法分别为：
[0100] 常规算法一：先进行P~0迭代，后进行Q~V迭代；
[0101] 常规算法二：先进行Q~V迭代，后进行P~0迭代；
[0102] 已有专利算法：专利ZL201410314990. 6所提出方法。
[0103] 图5是IEEE14节点系统，为了验证小阻抗支路对算法收敛性的影响，把算例中节 点4与节点7之间的支路1 47改为小阻抗支路，支路1 47的变比k= 0. 978,位于节点4侧。
[0104] -、IEEE14节点算例的计算结果
[0105] 支路147的阻抗取不同值时，四种算法潮流计算的迭代结果见表1。
[0106] 表1IEEE14节点算例不同支路阻抗时四种算法的迭代结果
[0108] 由表1可见，对于IEEE14节点系统算例，当小阻抗支路147的r/x= 1(T2时，常规 算法一就已经不收敛了，而本发明算法能够收敛；当小阻抗支路1 47的1~八=10 3时，常规算 法二已经不收敛了；当小阻抗支路147的r/x= 5X10 4时，已有专利算法还收敛，但迭代次 数很多，本发明算法通过适当的修正因子可以快速地收敛。可见本发明算法能有效处理电 阻较大的小阻抗支路。
[0109] IEEE14节点系统算例结果表明，当小阻抗支路的电阻较大时，两种常规算法不收 敛，但本发明算法能够收敛，迭代次数小于已有专利算法。
[0110] (2)445节点实际算例的计算结果
[0111] 445节点实际大型电网有445个节点，含有大量的小阻抗支路。其中，0.01的 小阻抗支路有118条，x< 0. 001的小阻抗支路有49条，x< 0. 0001的小阻抗支路有41 条，x< 0. 00001的小阻抗支路有22条。其中阻抗值最小的是节点118和节点125之间的 小阻抗支路为x= 0. 00000001，变比k= 0. 9565,k位于节点118侧。为了验证小阻抗支 路对算法收敛性的影响，使小阻抗支路1118_ 125的阻抗取不同值，四种算法潮流计算的迭代 结果见表2。
[0112] 表2 455节点系统不同支路阻抗时四种算法的迭代结果
[0113]
[0114] 由表2可见，对于455节点系统算例，当小阻抗支路1118_125的r/x= 1(T2时，常规 算法一就已经不收敛了，而本发明算法能够收敛；当小阻抗支路1 118%5的r/x= 10 3时，常 规算法二已经不收敛了；当小阻抗支路1118_125的1~八=10 4时，已有专利算法还收敛，但迭 代次数很多，本发明算法通过适当的修正因子可以快速地收敛。当小阻抗支路1118^ 125的r/ x =5X104时，已有专利算法不收敛，本发明算法通过适当的修正因子仍可快速地收敛。可 见本发明算法能有效处理电阻较大的小阻抗支路。
[0115] 455节点系统算例结果表明，当小阻抗支路的电阻较大时，两种常规算法不收敛， 但本发明算法能够收敛，迭代次数小于已有专利算法。
[0116] 本算法可以米用任何一种编程语目和编程环境实现，如C语目、C++、FORTRAN、 Delphi等。开发环境可以米用VisualC++、BorlandC++Builder、VisualFORTRAN等。
【主权项】
1. 一种修正系数矩阵的快速分解法潮流计算方法，采用先进行Q~V迭代，后进行P~ Θ迭代的方法，具体方法包括以下步骤： A、 原始数据输入和电压初始化； B、 形成节点导纳矩阵； C、 形成修正方程的系数矩阵B'和B"并进行因子表分解； 对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：形成系数矩阵B'时，给小 阻抗变压器支路的电抗乘一个修正因子α，正常支路的电抗不乘这个修正因子；小阻抗支 路相关的B'元素如下：式中，8'"、8\_、8、、8、是快速分解法系数矩阵8'的元素出，j(l是快速分解 法系数矩阵B'中不含小阻抗支路时的元素；k为变压器非标准变比；α为系数矩阵B'的 修正因子；\」是小阻抗变压器支路的电抗； 形成系数矩阵B'的方法，包括以下步骤： C1、读入支路数据，设置小阻抗阈值电阻rmin和电抗Xmin、设置修正因子α ; C2、设置支路计数初值m = 1 ; C3、取支路m的首末节点号i和j、电阻Α」、电抗xiP变比k ; C4、判断支路的电阻和电抗是否同时满足i^〈rmin和X yXmin，如果不满足转至步骤C7 ; C5、令 Xij= α X ij; C6、令 B' U= B' id-l/GiXij)、B' i」=l/(kx υ)，转至步骤 C8 ; C7、令 B' U= B' io-l/Xij'B' u= 1/x ij; C8、令 B'jj=B'』0-1/χυ、Β'Μ= l/xij; C9、令 m = m+1 ; CIO、判断m是否大于支路数1，如果m不大于1，则返回到步骤C3 ;否则，形成系数矩阵 B"并进行因子表分解； D、 设置迭代计数t = 0,收敛标志Kp= 0, K Q= 0 ; E、 计算无功功率不平衡量AQ; F、 判断无功功率最大不平衡量I AQmaxI是否小于收敛精度ε ;如果小于收敛精度ε， 令Kq= 1，转到步骤G ;否则，解方程Β〃 Λ V = Λ Q/V，修正电压幅值，令K Q= 0,转到步骤H ; G、 判断Kp是否等于1 ;如果Kp= 1，转到步骤K ; H、 计算有功功率不平衡量ΔΡ ; I、 判断有功功率最大不平衡量I ΔΡ_|是否小于收敛精度ε ;如果小于收敛精度ε， 令Kp= 1 ;否则，解方程Β' Λ θ = ΛΡ/ν，修正电压相角，令Kp= O ; J、 判断是否同时满足Kp= 1和K Q= 1，如果不满足，则令t = t+Ι，返回步骤E进行下 一次迭代； K、 计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
【专利摘要】本发明公开了一种修正系数矩阵的快速分解法潮流计算方法，采用先进行Q～V迭代，后进行P～θ迭代的方法，并对P～θ迭代修正方程的系数矩阵B′的部分元素进行修改。形成系数矩阵B′时，给小阻抗变压器支路的电抗乘一个修正因子α，正常支路的电抗不乘这个修正因子。本发明通过修改修正方程系数矩阵B′中与小阻抗相关的部分元素，改善了快速分解法潮流计算在分析含有电阻较大的小阻抗支路系统时的收敛性问题。采用常规快速分解法潮流计算或现有专利算法不收敛时，本发明能够可靠收敛。由于本发明不仅能有效解决常规快速分解法潮流算法分析含有电阻较大的小阻抗支路系统的收敛性问题，同时也能对正常系统进行潮流计算，因此没有不良影响。
【IPC分类】G06F17/50
【公开号】CN104899396
【申请号】CN201510346974
【发明人】姚玉斌, 阳义青, 吴志良, 王丹
【申请人】大连海事大学
【公开日】2015年9月9日
【申请日】2015年6月19日
【技术领域】
[0001] 本发明涉及一种电力系统的快速分解法潮流计算方法，特别是一种适合含小阻抗 支路系统的快速分解法潮流计算方法。
【背景技术】
[0002] 电力系统潮流计算是研宄电力系统稳态运行的一项基本计算，它根据给定的运行 条件和网络结构确定整个网络的运行状态。潮流计算也是其他电力系统分析的基础，如安 全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度快及内存需求少 的优点，快速分解法成为当前潮流计算的主流算法之一。
[0003] 当电力网络不存在小阻抗支路或电力网络中的小阻抗支路的电阻很小时，快速分 解法潮流计算具有良好的收敛性，但电网中存在电阻相对较大的小阻抗支路时，快速分解 法潮流计算就可能发散。电力系统小阻抗支路可以分为小阻抗线路和小阻抗变压器支路， 在数学模型上线路可以看作变比为1:1的变压器，因此下面分析时仅以小阻抗变压器支路 为例分析。变压器支路、模型见图4,变压器的非标准变比k位于节点i侧，阻抗位于标 准变比侧。当变压器的电阻ru和电抗Xij都很小时，此变压器支路即为小阻抗变压器支路。 变压器阻抗Zij=rij+jXij，导纳为：
[0005] 由于小阻抗变压器支路的阻抗很小，其阻抗上的电压降也很小，因此小阻抗变压 器支路两端的电压相量应满足：
[0007] 如果电压相量用极坐标表示，贝1J有
[0009] 式中，K^分别为节点i和节点j的电压相量；Vi、\分别为节点i和节点j的 电压幅值；0i、0 0>别为节点i和节点j的电压相角。
[0010] 如图1所示，现有快速分解法潮流计算方法，主要包括以下步骤：
[0011] A、原始数据输入和电压初始化；
[0012] 电压初始化一般采用平启动，即PV节点和平衡节点的电压幅值取给定值，PQ节点 的电压幅值取1. 〇 ;所有电压的相角都取〇. 0。这里单位采用标幺值。
[0013]B、形成节点导纳矩阵；
[0014] 设节点i和节点j原来的自电导与自电纳分别为61(|』^|、6」(|』」 (|，在它们之间增加 一条小阻抗支路后的自导纳和Y"、互导纳分别为：
[0018] C、形成修正方程的系数矩阵B'和B"并进行因子表分解；
[0019] 快速分解法修正方程为：
[0020] B，A0 =AP/V (7)
[0021] B"AV=AQ/V(8)
[0022] 式中，AP/V和AQ/V分别为有功功率和无功功率不平衡量除以电压幅值后的列 向量；AV和△ 0分别为电压幅值和电压相角修正量列向量；B'为导纳矩阵的虚部，但计 算时不计及支路电阻、对地导纳和非标准变比，矩阵中包含PQ节点和PV节点相关的行和 列；B"为导纳矩阵的虚部，仅包括与PQ节点有关的行和列。
[0023] 与小阻抗支路L相关的系数矩阵元素为：
[0030] 式中，B'n、B'jj、B'u是快速分解法系数矩阵B'的元素；B' j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素；B"wB" "、B"u是快速分解法系数矩 阵B"的元素；B" #是快速分解法系数矩阵B"中不含小阻抗支路时的元素。
[0031] D、设置迭代计数t= 0,收敛标志KP= 0,KQ= 0 ;
[0032] E、计算有功功率不平衡量AP;
[0033] 有功功率不平衡量（不包含平衡节点）为：
[0035] 式中，PiS为节点i的给定有功功率；VJP0汾别为节点i的电压幅值和电压相 角；Gu和Bu分别为导纳矩阵的电导部分和电纳部分；n为节点数。
[0036] F、判断有功功率最大不平衡量|APmax|是否小于收敛精度e;如果小于收敛精度 e，令KP= 1，转到步骤G;否则，解方程B'A0 =AP/V，修正电压相角，令KP= 0,转到步 骤H;
[0037]求解方程B'A0 =AP/V，得到A0，按下式修正电压相角：
[0038] 0 (t+D= 0W_A0 (t) (IQ)
[0039] G、判断心是否等于1 ;如果KQ= 1，转到步骤K;
[0040] H、计算无功功率不平衡量AQ;
[0041] 无功功率不平衡量（仅包含PQ节点）为：
[0043] 式中，QiS为节点i的给定无功功率；m为PQ节点数。
[0044] I、判断无功功率最大不平衡量|AQ_|是否小于收敛精度e;如果小于收敛精度 e，令&= 1 ;否则，解方程8〃八￥ =AQ/V，修正电压幅值，令KQ= 0 ;
[0045] 求解方程B〃AV=AQ/V，得到AV，按下式修正电压幅值：
[0046] y(t+1)=V(t)-AV(t) (18)
[0047] J、判断是否同时满足KP= 1和KQ= 1，如果不满足，令t=t+1，返回步骤E进行 下一次迭代；
[0048] K、计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
[0049] 步骤E和步骤F为P~0迭代，即通过AP求A0进而修正0 ;步骤H和步骤I 为Q~V迭代，即通过AQ求AV进而修正V。主流快速分解法都是按上述步骤设计算法， 即先进行P~9迭代，后进行Q~V迭代。也有文献采用先进行Q~V迭代，后进行P~ 9迭代的算法。
[0050] 对正常电力网络或含有电阻非常小的小阻抗支路的网络，快速分解法潮流计算具 有良好的收敛性，但遇到含有电阻较大的小阻抗的病态网络时，快速分解法潮流计算就可 能发散。电力系统中小阻抗支路普遍存在，潮流计算的收敛性是电力系统潮流计算这类非 线性问题的最重要指标，计算不收敛就无法得到问题的解。因此改善快速分解法潮流计算 针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。中国专利ZL201410314990. 6 提出的一种适合含小阻抗支路系统的快速分解法潮流计算方法如下：
[0051] (1)采用先进行Q~V迭代，后进行P~0迭代的方法；
[0052](2)对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：修改小阻抗变压 器支路非标准变比侧节点i对应的系数矩阵B'元素B'n、B'y标准变比侧节点j对应的 元素B'ji、B'不变。
[0057] 式中，8'"、8、、8\、8'^是快速分解法系数矩阵"的元素出'i(l、B' j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素。
[0058] 中国专利ZL201410314990. 6所提出方法有效提高了含有小阻抗支路电力系统潮 流计算的收敛性，但当小阻抗支路的电阻较大（rijAXij)时，该方法仍可能发散。
【发明内容】
[0059] 为解决现有技术存在的上述问题，本发明要提出一种修正系数矩阵的快速分解法 潮流计算方法，以改善其分析含有电阻较大的小阻抗支路电力系统的潮流计算收敛性。
[0060] 为了实现上述目的，本发明提出了适合含电阻较大的小阻抗支路系统的快速分解 法潮流计算方法来改善潮流计算收敛性。本发明的技术方案如下：一种修正系数矩阵的快 速分解法潮流计算方法，采用先进行Q~V迭代，后进行P~0迭代的方法，并对P~0 迭代修正方程的系数矩阵B'的部分元素进行修改。方案包括以下步骤：
[0061] A、原始数据输入和电压初始化；
[0062] B、形成节点导纳矩阵；
[0063] C、形成修正方程的系数矩阵B'和B"并进行因子表分解；
[0064] 对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：形成系数矩阵B'时， 给小阻抗变压器支路的电抗乘一个修正因子a，正常支路的电抗不乘这个修正因子。小阻 抗支路相关的B'元素如下：
[0069] 式中，8'"、8、、8\、8'^是快速分解法系数矩阵"的元素出'j(l是快速 分解法系数矩阵B'中不含小阻抗支路时的元素；k为变压器非标准变比；a为系数矩阵 B'的修正因子。
[0070] 形成系数矩阵B'的方法，包括以下步骤：
[0071]C1、读入支路数据，设置小阻抗阈值电阻1^和电抗xmin、设置修正因子a。
[0072] C2、设置支路计数初值m= 1。
[0073] C3、取支路m的首末节点号i和j、电阻rij、电抗Xij、变比k。
[0074] C4、判断支路的电阻和电抗是否同时满足和x ，如果不满足转至步骤 C7〇
[0075]C5、令叉。=ax。。
[0076] C6、令B'B' ^-1/(10^)、B'l/(kx。），转至步骤C8。
[0077] C7、令B，B，io-l/Xj、B，。.= 1/xu。
[0078] C8、令B'B'jo-l/Xij、B'1/xij〇
[0079] C9、令m=m+l。
[0080] C10、判断m是否大于支路数1，如果m不大于1，则返回到步骤C3 ;否则，形成系数 矩阵B"并进行因子表分解。
[0081]D、设置迭代计数t= 0,收敛标志KP= 0,KQ= 0;
[0082]E、计算无功功率不平衡量AQ;
[ 0083] F、判断无功功率最大不平衡量|AQ_|是否小于收敛精度e;如果小于收敛精度 e，令&= 1，转到步骤G;否则，解方程^八乂 =AQ/V，修正电压幅值，令KQ= 0,转到步骤 H；
[0084] G、判断KP是否等于1 ;如果KP= 1，转到步骤K;
[0085]H、计算有功功率不平衡量AP;
[0086] I、判断有功功率最大不平衡量|APmax|是否小于收敛精度e;如果小于收敛精度 e，令KP= 1 ;否则，解方程B'A0 =AP/V，修正电压相角，令KP= 0 ;
[0087] J、判断是否同时满足KP= 1和KQ= 1，如果不满足，贝令t=t+1，返回步骤E进 行下一次迭代；
[0088] K、计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
[0089] 与现有技术相比，本发明具有以下有益效果：
[0090] 1、本发明通过修改快速分解法修正方程系数矩阵B'中与小阻抗相关的部分元 素，改善了快速分解法潮流计算在分析含有电阻较大的小阻抗支路系统时的收敛性问题。 采用常规快速分解法潮流计算或现有专利算法不收敛时，本发明能够可靠收敛。
[0091] 2、由于本发明不仅能有效解决常规快速分解法潮流算法分析含有电阻较大的小 阻抗支路系统的收敛性问题，同时也能对正常系统进行潮流计算，因此没有不良影响。
【附图说明】
[0092] 本发明共有附图5张。其中：
[0093] 图1是现有快速分解法潮流计算的流程图。
[0094] 图2是本发明快速分解法潮流计算的流程图。
[0095]图3是本发明形成系数矩阵B'的流程图。
[0096] 图4是电力系统变压器模型示意图。
[0097] 图5是IEEE14节点电力系统算例的接线图。
【具体实施方式】
[0098] 下面结合附图对本发明进行进一步地说明，按照图2所示流程对IEEE14节点系统 (电气和电子工程师协会14节点系统测试数据）和一个445节点实际系统两个算例进行了 计算，作为对比同时采用常规快速分解法和已有专利算法对该算例进行了计算。潮流计算 的收敛精度为0.00001。
[0099] 两种常规快速分解法分别为：
[0100] 常规算法一：先进行P~0迭代，后进行Q~V迭代；
[0101] 常规算法二：先进行Q~V迭代，后进行P~0迭代；
[0102] 已有专利算法：专利ZL201410314990. 6所提出方法。
[0103] 图5是IEEE14节点系统，为了验证小阻抗支路对算法收敛性的影响，把算例中节 点4与节点7之间的支路1 47改为小阻抗支路，支路1 47的变比k= 0. 978,位于节点4侧。
[0104] -、IEEE14节点算例的计算结果
[0105] 支路147的阻抗取不同值时，四种算法潮流计算的迭代结果见表1。
[0106] 表1IEEE14节点算例不同支路阻抗时四种算法的迭代结果
[0108] 由表1可见，对于IEEE14节点系统算例，当小阻抗支路147的r/x= 1(T2时，常规 算法一就已经不收敛了，而本发明算法能够收敛；当小阻抗支路1 47的1~八=10 3时，常规算 法二已经不收敛了；当小阻抗支路147的r/x= 5X10 4时，已有专利算法还收敛，但迭代次 数很多，本发明算法通过适当的修正因子可以快速地收敛。可见本发明算法能有效处理电 阻较大的小阻抗支路。
[0109] IEEE14节点系统算例结果表明，当小阻抗支路的电阻较大时，两种常规算法不收 敛，但本发明算法能够收敛，迭代次数小于已有专利算法。
[0110] (2)445节点实际算例的计算结果
[0111] 445节点实际大型电网有445个节点，含有大量的小阻抗支路。其中，0.01的 小阻抗支路有118条，x< 0. 001的小阻抗支路有49条，x< 0. 0001的小阻抗支路有41 条，x< 0. 00001的小阻抗支路有22条。其中阻抗值最小的是节点118和节点125之间的 小阻抗支路为x= 0. 00000001，变比k= 0. 9565,k位于节点118侧。为了验证小阻抗支 路对算法收敛性的影响，使小阻抗支路1118_ 125的阻抗取不同值，四种算法潮流计算的迭代 结果见表2。
[0112] 表2 455节点系统不同支路阻抗时四种算法的迭代结果
[0113]
[0114] 由表2可见，对于455节点系统算例，当小阻抗支路1118_125的r/x= 1(T2时，常规 算法一就已经不收敛了，而本发明算法能够收敛；当小阻抗支路1 118%5的r/x= 10 3时，常 规算法二已经不收敛了；当小阻抗支路1118_125的1~八=10 4时，已有专利算法还收敛，但迭 代次数很多，本发明算法通过适当的修正因子可以快速地收敛。当小阻抗支路1118^ 125的r/ x =5X104时，已有专利算法不收敛，本发明算法通过适当的修正因子仍可快速地收敛。可 见本发明算法能有效处理电阻较大的小阻抗支路。
[0115] 455节点系统算例结果表明，当小阻抗支路的电阻较大时，两种常规算法不收敛， 但本发明算法能够收敛，迭代次数小于已有专利算法。
[0116] 本算法可以米用任何一种编程语目和编程环境实现，如C语目、C++、FORTRAN、 Delphi等。开发环境可以米用VisualC++、BorlandC++Builder、VisualFORTRAN等。
【主权项】
1. 一种修正系数矩阵的快速分解法潮流计算方法，采用先进行Q~V迭代，后进行P~ Θ迭代的方法，具体方法包括以下步骤： A、 原始数据输入和电压初始化； B、 形成节点导纳矩阵； C、 形成修正方程的系数矩阵B'和B"并进行因子表分解； 对系数矩阵B'中与小阻抗支路有关的部分元素进行修改：形成系数矩阵B'时，给小 阻抗变压器支路的电抗乘一个修正因子α，正常支路的电抗不乘这个修正因子；小阻抗支 路相关的B'元素如下：式中，8'"、8\_、8、、8、是快速分解法系数矩阵8'的元素出，j(l是快速分解 法系数矩阵B'中不含小阻抗支路时的元素；k为变压器非标准变比；α为系数矩阵B'的 修正因子；\」是小阻抗变压器支路的电抗； 形成系数矩阵B'的方法，包括以下步骤： C1、读入支路数据，设置小阻抗阈值电阻rmin和电抗Xmin、设置修正因子α ; C2、设置支路计数初值m = 1 ; C3、取支路m的首末节点号i和j、电阻Α」、电抗xiP变比k ; C4、判断支路的电阻和电抗是否同时满足i^〈rmin和X yXmin，如果不满足转至步骤C7 ; C5、令 Xij= α X ij; C6、令 B' U= B' id-l/GiXij)、B' i」=l/(kx υ)，转至步骤 C8 ; C7、令 B' U= B' io-l/Xij'B' u= 1/x ij; C8、令 B'jj=B'』0-1/χυ、Β'Μ= l/xij; C9、令 m = m+1 ; CIO、判断m是否大于支路数1，如果m不大于1，则返回到步骤C3 ;否则，形成系数矩阵 B"并进行因子表分解； D、 设置迭代计数t = 0,收敛标志Kp= 0, K Q= 0 ; E、 计算无功功率不平衡量AQ; F、 判断无功功率最大不平衡量I AQmaxI是否小于收敛精度ε ;如果小于收敛精度ε， 令Kq= 1，转到步骤G ;否则，解方程Β〃 Λ V = Λ Q/V，修正电压幅值，令K Q= 0,转到步骤H ; G、 判断Kp是否等于1 ;如果Kp= 1，转到步骤K ; H、 计算有功功率不平衡量ΔΡ ; I、 判断有功功率最大不平衡量I ΔΡ_|是否小于收敛精度ε ;如果小于收敛精度ε， 令Kp= 1 ;否则，解方程Β' Λ θ = ΛΡ/ν，修正电压相角，令Kp= O ; J、 判断是否同时满足Kp= 1和K Q= 1，如果不满足，则令t = t+Ι，返回步骤E进行下 一次迭代； K、 计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。
【专利摘要】本发明公开了一种修正系数矩阵的快速分解法潮流计算方法，采用先进行Q～V迭代，后进行P～θ迭代的方法，并对P～θ迭代修正方程的系数矩阵B′的部分元素进行修改。形成系数矩阵B′时，给小阻抗变压器支路的电抗乘一个修正因子α，正常支路的电抗不乘这个修正因子。本发明通过修改修正方程系数矩阵B′中与小阻抗相关的部分元素，改善了快速分解法潮流计算在分析含有电阻较大的小阻抗支路系统时的收敛性问题。采用常规快速分解法潮流计算或现有专利算法不收敛时，本发明能够可靠收敛。由于本发明不仅能有效解决常规快速分解法潮流算法分析含有电阻较大的小阻抗支路系统的收敛性问题，同时也能对正常系统进行潮流计算，因此没有不良影响。
【IPC分类】G06F17/50
【公开号】CN104899396
【申请号】CN201510346974
【发明人】姚玉斌, 阳义青, 吴志良, 王丹
【申请人】大连海事大学
【公开日】2015年9月9日
【申请日】2015年6月19日
最新回复(0)