含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法技术

本发明专利技术公开了一种小阻抗支路PQ端点变雅可比矩阵的牛顿法潮流计算方法，包括以下步骤：根据支路电阻和电抗的大小确定两端节点所连支路类型T；小阻抗支路的PQ端点与其他节点采用不同的方法计算雅可比矩阵；解修正方程。本发明专利技术通过在迭代过程中小阻抗支路的PQ端点采用与其他节点不同的雅可比矩阵计算方法，解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路电力系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时，本发明专利技术能够可靠收敛，且比现有专利技术迭代次数少。由于本发明专利技术不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题，同时也能对正常电力系统进行潮流计算，没有不良影响。

Newton Raphson method for power flow calculation of power network with small impedance branches

The invention discloses a small impedance matrix of the PQ endpoint Jacobi Newton power flow calculation method, which comprises the following steps: according to the branch resistance and reactance is determined by the size of the end nodes connected branch type T; PQ endpoint small impedance branches with other nodes using different methods to calculate the Jacobi matrix solution of the correction equation. The PQ endpoint iterative process of small impedance branches using Jacobi matrix and other node calculation method to solve the rectangular coordinate Newton power flow calculation convergence problem with small impedance branches in the analysis of power system. When the power flow calculation by the conventional Cartesian coordinate method does not converge, the invention can converge reliably and has fewer iterations than the existing patent technology. The invention not only can effectively solve the conventional rectangular coordinate Newton power flow calculation and analysis of the convergence problem with small impedance branches of power system, but also can compute the power flow of power system normal, no adverse effect.

【技术实现步骤摘要】
含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法
本专利技术涉及含小阻抗支路电力系统的潮流计算方法，特别是一种电力系统的极坐标牛顿法潮流计算方法。
技术介绍
电力系统潮流计算是研究电力系统稳态运行的一项基本计算，它根据电力系统给定的运行条件和网络结构确定整个电力系统的运行状态。潮流计算也是电力系统其他分析的基础，如安全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度较快及内存需求适中的优点，牛顿法成为当前潮流计算的主流方法。牛顿法分为极坐标和直角坐标两种形式，两种形式的牛顿法潮流计算都在电力系统中得到了广泛的应用。在极坐标牛顿法潮流计算中，节点i的电压采用极坐标表示为：对正常电力网络，极坐标牛顿法潮流计算具有良好的收敛性，但遇到含有小阻抗支路的病态网络时，极坐标牛顿法潮流计算就可能发散。电力系统小阻抗支路可分为小阻抗线路和小阻抗变压器支路，在数学模型上线路可以看作变比为1:1的变压器，因此下面分析时仅以小阻抗变压器支路为例分析。小阻抗变压器模型见图1，变压器的非标准变比k位于节点i侧，阻抗位于标准变比侧。变压器阻抗zij＝rij+jxij很小，导纳为式中，yij、gij、bij分别为节点i和节点j之间小阻抗支路的导纳、电导和电纳；rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗。由于小阻抗支路li-j的阻抗很小，支路的电压降也很小，因此变压器两端节点的电压应满足：式中，Ui、θi分别为节点i的电压幅值和相角；Uj、θj分别为节点j的电压幅值和相角。如图2所示，现有极坐标牛顿法潮流计算方法，主要包括以下步骤：A、输入原始数据和初始化电压根据电力系统节点的特点，潮流计算把电力系统节点分成3类：节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点；节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点；节点电压幅值和电压相角已知，节点有功功率和无功功率未知的节点称为平衡节点。电压初始化采用平启动，即PV节点和平衡节点的电压幅值取给定值，PQ节点的电压幅值取1.0；所有电压相角都取0.0。这里相角单位为弧度，其他量采用标幺值。B、形成节点导纳矩阵设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0，在它们之间增加一条小阻抗支路后的自导纳和互导纳分别为：式中，Yii、Yjj分别为节点i和节点j的自导纳；Yij为节点i和节点j之间的互导纳；rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗；k为节点i和节点j之间小阻抗支路的变比(如果是输电线支路，变比为1)；C、计算功率偏差节点功率公式为：式中，Pi、Qi分别为节点i的节点有功功率和无功功率；Ui、Uk分别为节点i和节点k的节点电压幅值；θik＝θi-θk，θi和θk分别为节点i和节点k的节点电压相角；Gik、Bik分别为节点导纳矩阵元素Yik的实部和虚部；n为电力系统的节点数。节点功率偏差计算公式为：式中，ΔPi、ΔQi分别为节点i的节点有功功率偏差和无功功率偏差；Pis、Qis分别为节点i给定的节点注入有功功率和注入无功功率；m为PQ节点数。PV节点的无功功率不是给定值，不需要计算无功功率偏差；平衡节点不参与迭代计算，不需要计算有功功率偏差和无功功率偏差。求各节点功率偏差中绝对值最大的值，称为最大不平衡量，如果最大不平衡量的绝对值小于给定的收敛精度，转步骤F，否则执行步骤D。D、形成雅可比矩阵J雅可比矩阵J的元素(i≠j时)计算公式如下：雅可比矩阵J的元素(i＝j时)计算公式如下：E、解修正方程及修正电压幅值U、相角θ潮流计算的基本方程(7)是非线性方程组，通常采用逐次线性化方法迭代求解。线性化得到的方程称为修正方程，用来求电压幅值和相角的修正量。修正方程为：式中，J为雅可比矩阵；H、N、M、L为雅可比矩阵的分块子矩阵；ΔP和ΔQ分别为有功功率和无功功率偏差列向量；Δθ为电压相角修正量列向量；ΔU/U为电压幅值修正量与电压幅值之比的列向量。电压修正公式为：式中，上标t表示第t次迭代。F、输出节点及支路数据。对正常电力网络，牛顿法潮流计算具有良好的收敛性，但遇到含有小阻抗支路的病态网络时，牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在，收敛性是电力系统潮流计算这类非线性问题的最重要指标，计算不收敛就无法得到方程的解。因此改善极坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。专利技术者在早期研究中提出了一种变雅可比牛顿法潮流计算方法，该方法采用给定值Pis和Qis代替计算值Pi和Qi计算雅可比矩阵的Nii和Lii元素，而雅可比矩阵的Hii和Mii元素仍由Pi和Qi计算。变雅可比牛顿法潮流计算方法有效解决了极坐标牛顿法潮流计算分析包含小阻抗支路的电力系统时的发散问题。现在仍然在使用的一些极坐标牛顿法潮流计算程序不具备处理小阻抗支路的能力，但对这些老的潮流计算程序进行改造往往比较麻烦。有些程序可能采用比较老的编程语言编写，不易修改，有的潮流计算版本甚至没有源文件，无法修改。
技术实现思路
为解决现有技术存在的上述问题，本专利技术要提出一种含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法，以改善极坐标牛顿法潮流计算方法分析含有小阻抗支路电力系统的潮流计算的收敛性。为了实现上述目的，本专利技术的技术方案如下：含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法，在进行潮流计算之前，先对小阻抗支路进行处理，把此小阻抗支路变成两个支路串联的形式，其中一个支路的阻抗为z1＝rij+j(xij+xc)，另一个支路的阻抗为z2＝-jxc，这里xc为电力系统中电抗绝对值大于小电抗阈值的支路的电抗绝对值的平均值，称为电力系统正常电抗均值xav。经过串联补偿把一个小阻抗变成两个阻抗较大的支路，可以提高潮流计算的收敛性。方案包括以下步骤：A、输入原始数据；B、对小阻抗支路进行串联补偿；用电力系统正常电抗均值对小阻抗支路进行补偿，得到两个阻抗较大的支路，其中一个支路的阻抗为z1＝rij+j(xij+xc)，另一个支路的阻抗为z2＝-jxc，这里xc为电力系统正常电抗均值xav。小阻抗支路串联补偿的方法，包括以下步骤：B1、读入电力系统所有线路和变压器支路数据，设置小电阻阈值rmin和小电抗阈值xmin。B2、计算电力系统正常电抗均值xav。B3、设置支路计数初值m＝1。B4、设置新增加支路和节点计数初值p＝0。B5、取支路m的首末节点号i和j、电阻rij、电抗xij、变比k。B6、判断是否满足rij≤rmin且xij≤xmin的条件，如果不满足转至步骤B12。B7、令xc＝xav。B8、令p＝p+1。B9、增加节点号为n+p的节点，节点类型设置为PQ节点，节点电源有功功率和无功功率及负荷有功功率和无功功率都设为0。B10、设支路m的末节点号为n+p、电抗为xij+xc，其它不变。B11、增加支路n+p，令其首末节点号分别为n+p和j、电阻为0、电抗为-xc、变比为1.0。B12、令m＝m+1。B13、判断m是否大于支路数l，如果m不大于l，则返回到步骤B5；否则，转至步骤C。C、初始化电压；D、形成节点导纳矩阵；E、设置迭代计数t＝0；F、计算功率偏差，求最大不平衡量ΔW本文档来自技高网...

【技术保护点】
用于含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法，其特征在于：采用串联补偿方式对小阻抗支路进行补偿，把每个小阻抗支路变成两个正常阻抗的支路，然后进行极坐标牛顿法潮流计算；具体包括以下步骤：A、原始数据输入；B、对小阻抗支路进行串联补偿；用电力系统正常电抗均值对小阻抗支路进行补偿，得到两个阻抗较大的支路，其中一个支路的阻抗为z

【技术特征摘要】
1.用于含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法，其特征在于：采用串联补偿方式对小阻抗支路进行补偿，把每个小阻抗支路变成两个正常阻抗的支路，然后进行极坐标牛顿法潮流计算；具体包括以下步骤：A、原始数据输入；B、对小阻抗支路进行串联补偿；用电力系统正常电抗均值对小阻抗支路进行补偿，得到两个阻抗较大的支路，其中一个支路的阻抗为z1＝r+j(x+xc)，另一个支路的阻抗为z2＝-jxc，这里xc为电力系统正常电抗均值xav；小阻抗支路串联补偿的方法，包括以下步骤：B1、读入电力系统所有线路和变压器支路数据，设置小阻抗阈值电阻rmin和阈值电抗xmin；B2、计算电力系统正常电抗均值xav；B3、设置支路计数初值m＝1；B4、设置新增加支路和节点计数初值p＝0；B5、取支路m的首末节点号i和j、电阻r、电抗x、变比k；B6、判断是否满足r≤rmin且x≤xmin的条件，如果不满足转至步骤B12；B7、令xc＝xav；B8、令p＝p+1；B9、增加节点号为n+p的节点，节点类型设置为PQ节点，节点电源有功功率和无功功率及负荷有功功率和无功功率都设为0；B10、设支路m的末节点号为n+p、电抗为x+xc，其它不变；B11、增加支路n+p，令其首末节点号分别为n+p和j、电阻为0、电抗为-xc、变比为1.0；B12、令m＝m+1；B13、判断m是否大于支路数l，如果m不大于l，则返回到步骤B5；否则，转至步骤C；C、初始化电压；根据电力系统节点的特点，潮流计算把电力系统节点分成3类：节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点；节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点；节点电压幅值和电压相角已知，节点有功功率和无功功率未知的节点称为平衡节点；电压初始化采用平启动，即PV节点和平衡节点的电压幅值取给定值，PQ节点的电压幅值取1.0；所有电压的相角都取0.0；这里相角单位为弧度，其他量单位采用标幺值；D、形成节点导纳矩阵；E、设置迭代计数t＝0；F、计算功率偏差，求最大不平衡量ΔWmax；节点功率公式为：式中，Pi、Qi分别为节点i的节点有功功率和无功功率；Ui、Uk分别为节点i和节点k的节点电压幅值；θik＝θi-θk，θi和θk分别为节点i和节点k的节点电压相角；Gik、Bik分别为节点导纳矩阵元素Yik的实部和虚部；n为电力系统的节点数；节点功率偏差计算公式为：

【专利技术属性】
技术研发人员：姚玉斌，吴志良，王丹，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：辽宁,21