[发明专利]一种含小阻抗支路电网串联补偿快速分解法潮流计算方法

摘要

本发明公开了一种含小阻抗支路电网串联补偿快速分解法潮流计算方法，在进行潮流计算之前，先对电阻值大于电抗绝对值的小阻抗支路进行处理，把此小阻抗支路变成两个支路串联的形式，其中一个支路的阻抗为z₁＝r+j(x+x_av)，另一个支路的阻抗为z₂＝‑jx_av，这里x_av为电力系统正常电抗均值。经过串联补偿把一个小阻抗变成两个阻抗较大的支路，可以提高潮流计算的收敛性。能够计算任意阻抗的小阻抗支路。本发明仅对r>|x|的小阻抗支路进行补偿，有效减少了新增加的节点和支路个数，内存占用量比现有专利方法少，计算速度更快。本发明对数据进行预处理，不需要改变潮流计算程序。

主权项

1.一种含小阻抗支路电网串联补偿快速分解法潮流计算方法，其特征在于：采用串联补偿方式对电阻大于电抗绝对值的小阻抗支路进行补偿，把此小阻抗支路变成两个正常阻抗的支路，然后进行快速分解法潮流计算；具体包括以下步骤：A、原始数据输入；B、对小阻抗支路进行串联补偿；用电力系统正常电抗均值对电阻值大于电抗绝对值的小阻抗支路进行补偿，得到两个阻抗较大的支路，其中一个支路的阻抗为z1＝r+j(x+xav)，另一个支路的阻抗为z2＝‑jxav，这里xav为电力系统正常电抗均值；小阻抗支路串联补偿的方法，包括以下步骤：B1、读入电力系统所有线路和变压器支路数据，设置小阻抗阈值电阻rmin和阈值电抗xmin；B2、计算电力系统正常电抗均值xav；B3、设置支路计数初值q＝1；B4、设置新增加支路和节点计数初值p＝0；B5、取支路q的首末节点号i和j、电阻r、电抗x、变比k；B6、判断是否满足r≤rmin且|x|≤xmin且r>|x|的条件，如果不满足转至步骤B11；B7、令p＝p+1；B8、增加节点号为n+p的节点，节点类型设置为PQ节点，节点电源有功功率和无功功率及负荷有功功率和无功功率都设为0；B9、设支路q的末节点号为n+p、电抗为x+xav，其它不变；B10、增加支路n+p，令其首末节点号分别为n+p和j、电阻为0、电抗为‑xav、变比为1.0；B11、令q＝q+1；B12、判断q是否大于支路数l，如果q不大于l，则返回到步骤B5；否则，转至步骤C；C、初始化电压；根据电力系统节点的特点，潮流计算把电力系统节点分成3类：节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点；节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点；节点电压幅值和电压相角已知，节点有功功率和无功功率未知的节点称为平衡节点；电压初始化采用平启动，即PV节点和平衡节点的电压幅值取给定值，PQ节点的电压幅值取1.0；所有电压的相角都取0.0；这里相角单位为弧度，其他量单位采用标幺值；D、形成节点导纳矩阵；E、形成快速分解法修正方程的系数矩阵B′和B″并进行因子表分解；快速分解法修正方程为：B′Δθ＝ΔP/V  (1)B″ΔV＝ΔQ/V  (2)式中，ΔP/V和ΔQ/V分别为有功功率和无功功率不平衡量除以电压幅值后的列向量；ΔV和Δθ分别为电压幅值和电压相角修正量列向量；B′为导纳矩阵的虚部，但计算时不计及支路电阻、对地导纳和非标准变比，导纳矩阵中包含PQ节点和PV节点相关的行和列；B″为导纳矩阵的虚部，仅包括与PQ节点有关的行和列；F、设置迭代计数t＝0，收敛标志KP＝0，KQ＝0；G、计算PQ节点的无功功率不平衡量ΔQ，并求无功功率最大不平衡量ΔQmax；PQ节点的无功功率不平衡量为：式中，Qis为节点i的给定无功功率；Vi为节点i的电压幅值；Vj为节点j的电压幅值，θij＝θi‑θj，θi、θj分别为节点i和节点j的电压相角；Gij和Bij分别为导纳矩阵元素的电导部分和电纳部分；n为节点数；m为PQ节点数；求各PQ节点中无功功率不平衡量绝对值最大的值，称为无功功率最大不平衡量，记为ΔQmax；H、判断无功功率最大不平衡量绝对值|ΔQmax|是否小于收敛精度ε；如果小于收敛精度ε，令KQ＝1，转到步骤I；否则，解快速分解法修正方程B"ΔV＝ΔQ/V，修正电压幅值，令KQ＝0，转到步骤J；求解快速分解法修正方程B″ΔV＝ΔQ/V，得到ΔV，按下式修正电压幅值：V(t+1)＝V(t)‑ΔV(t)  (4)式中，上标(t)表示第t次迭代；I、判断KP是否等于1；如果KP＝1，转到步骤N；J、计算PQ节点和PV节点的有功功率不平衡量ΔP，并求有功功率最大不平衡量ΔPmax；PQ节点和PV节点的有功功率不平衡量为：式中，Pis为节点i的给定有功功率；n为节点数；求各PQ节点和PV节点中有功功率不平衡量绝对值最大的值，称为有功功率最大不平衡量，记为ΔPmax；K、判断有功功率最大不平衡量绝对值|ΔPmax|是否小于收敛精度ε；如果小于收敛精度ε，令KP＝1，转到步骤L；否则，解快速分解法修正方程B'Δθ＝ΔP/V，修正电压相角，令KP＝0，转到步骤M；求解快速分解法修正方程B′Δθ＝ΔP/V，得到Δθ，按下式修正电压相角：θ(t+1)＝θ(t)‑Δθ(t)  (6)L、判断KQ是否等于1；如果KQ＝1，转到步骤N；M、令t＝t+1，返回步骤G进行下一次迭代；N、计算平衡节点功率及PV节点的无功功率，计算支路功率，结束。