一种抑制短路电流的联络线电抗快速计算方法与流程

本发明属于电力系统调度自动化
技术领域：
，具体涉及一种抑制短路电流的联络线电抗快速计算方法。
背景技术：
：城市电网在500kv网架不断发展和完善的推动下，逐步形成了主网为500kv环网，区域电网为220kv电网的网络构造。在此基础上，推进实行以500kv变电站为核心的电网分区运行。城市电网采用分区运行方式后，能够大幅提高系统输电能力，有效控制并降低网损，便于系统调整运行方式以适应电网发展，降低事故预案编制的难度。可见，实施分区方案是城市电网发展的必经阶段，可以显著促进电网的安全稳定运行。但是，对于具有非独立供电可靠性的分区，当分区内发生故障时如主变n-2等，分区内的主变可能会发生过载情况，需要考虑投入分区联络线以获取其他分区的功率支援，此时，由于联络线的投入会导致分区互联，分区内的厂站短路电流会上升，有可能超出该电压等级下的遮断电流。为了满足系统安全稳定运行的要求，需要采取一定的短路电流限制措施。技术实现要素：本发明的目的是提出一种抑制短路电流的联络线电抗快速计算方法，其特征在于，包括以下步骤：步骤a.基于系统网络参数形成节点导纳矩阵，对节点导纳矩阵求逆得到描述网络的节点阻抗矩阵，求的系统各节点的短路电流，获取短路水平最高的节点；步骤b.基于以上计算结果，由该节点引入一条支路，支路电抗为xk,利用支路追加法推导新网络的节点电抗矩阵，计算此时该节点的短路电流，在能够选支路中获得最优支路。所述步骤a包括：步骤a1：任一网络用节点阻抗矩阵表示的节点电压方程为u＝zi，i和u分别为：用矩阵表示如下：式(1)中，电压相量u为网络各节点对地电压；电流相量i为网络外部向各节点的注入电流；系数矩阵即节点阻抗矩阵z；上述矩阵的矩阵对角元素zii称为自阻抗，表达为：即当除节点i外所有其他节点注入电流均为零时，节点i电压与注入电流之比，也就是从节点i看进网络的等效阻抗；非对角元素zij称为互阻抗，表示为：即当除节点i外所有其他节点注入电流均为零时，节点j电压与节点i电流之比，或者说节点i单位电流在节点j引起的电压值；显然zij是不会为零的，即节点阻抗矩阵总是满阵；步骤a2：基于步骤a1形成的节点阻抗矩阵，假设某一点f发生短路故障时，短路点f的接地阻抗为zf＝0；该网络中只有短路点f有注入电流-if，即if由f点流向地；故节点电压方程为短路点电压故障分量为：由此求得短路电流：当zf＝0时，考虑到在输电网规划中，关心的是电力系统中有功潮流的分布，而不需要计算各节点的电压幅值，且对于超高压电力网，线路电阻比线路电抗小得多，此时可忽略所有的电阻，则式(7)简化为：所述步骤b包括如下步骤：步骤b1：原输电系统网络的节点阻抗矩阵为x＝b0-1,支路k两端的节点为l,m；这里的支路是指两节点间各线路的并联，线路是支路的一个元件；当支路k闭合时，相当于在l,m间增加一条电抗为xk的线路，形成新的网络；新网络的节点电抗矩阵为x′，并且设网络变化前后节点注入功率不变；其节点注入电流相量i和电压相量u，由于注入功率不变，故满足：u＝x′i(9)原网络的注入电流表示为：式中：ek为支路k的节点-支路关联矢量，仅在l,m位置对应元素为+1和-1，其余元素为0。由原网络知：节点l,m间的电压为：将式(11)代入式(12)，得：由此得到：将式(14)代入到式(11)中，得：将式(15)与(9)比较，得：从而节点电抗矩阵的变化量为：将xk换为xlm，ek换为alm，则式(16)为:实际追加的支路其xk为正值时追加的是电抗支路，xk为负值时追加的是电容支路。由公式(8)、(18)可知，若要减小节点短路电流，则x′要增大，应当追加电容支路；步骤b2：由新的节点阻抗矩阵可得短路电流最高节点此时的自阻抗为x′ff，代入公式(8)可得：选择不同的支路与xlm，即可得到取其最小值i′fmin即为所求最优解。本发明的有益效果是实际电网经常会由于网架结构不合理、负荷增长而出现短路电流超过最大值的情况，从而影响系统和设备的安全稳定运行。此时电网公司需要在合适的位置新增联络线以降低短路电流。本发明给出了一种可以快速计算得到联络线电阻的方法，可有效降低短路电流。附图说明图1为制短路电流的联络线电抗快速计算方法流程图。图2为单支路等值网络示意图。图3为14节点系统结构示意图。具体实施方式本发明的目的是提出一种抑制短路电流的联络线电抗快速计算方法，下面结合附图和实施例予以说明。图1所示为制短路电流的联络线电抗快速计算方法流程图。图中所示步骤为：步骤a.基于系统网络参数形成节点导纳矩阵，对节点导纳矩阵求逆得到描述网络的节点阻抗矩阵，求的系统各节点的短路电流，获取短路水平最高的节点；步骤b.基于以上计算结果，由该节点引入一条支路，支路电抗为xk,利用支路追加法推导新网络的节点电抗矩阵，计算此时该节点的短路电流，在能够选支路中获得最优支路(如图2所示)。以下面实施例详细描述本发明，图3所示为14节点系统结构示意图。本次算例基于matlab的14节点系统进行的；为了简化计算过程，在计算过程中有以下计算条件：1)假设系统有无限大的容量，短路后，系统母线电压维持不变；2)在计算短路电流时，只考虑发电机、线路等的电抗而忽略其电阻，且不考虑变压器支路的影响即令变压器变比为1，计算过程中的参数均为标幺值；3)短路电流计算过程中采取的短路形式为三相短路，因为单项短路与两相短路的短路电流均小于三相短路电流。步骤a1：基于系统的网络参数与拓扑结构形成节点导纳矩阵y＝g+jb，y是14×14阶的方阵，其中g为电导，b为电纳，令g＝0得y＝jb，对y求导得节点电抗矩阵x＝y-1，取x的对角元即可得各节点的自阻抗xii＝diag(x)。在对节点导纳矩阵求逆的中，由于导纳阵为奇异阵无法求逆，故去掉系统的平衡节点1，得到13×13阶的节点导纳阵y，求得的节点自阻抗为：xii＝[0.04960.14950.09540.08690.24850.22420.40030.23360.27930.31400.37880.31160.3656]t；步骤a2：代入公式(8)可得节点短路电流：if＝[20.176.6910.4911.514.024.462.504.283.583.182.673.212.73]t，ifmax＝if2＝20.17，因此节点2为短路水平最高的节点，针对该节点进行支路追加。步骤b1：原14节点系统共有20条支路，现以节点2作为支路一端节点新增一条支路，另一端可选节点为节点6、节点7、节点8、节点9、节点10、节点11、节点12、节点13、节点14九种情况。以节点6为支路末端追加支路，支路电抗xlm取为0.2。代入公式(18)得即可得新的节点电抗矩阵，取其对角元得以下节点自阻抗：xii＝[0.04920.14880.09120.08180.14310.20590.38210.20390.23960.24590.27610.21970.3129]t，步骤b2：代入公式(8)可得节点短路电流：if＝[20.316.7210.9612.236.984.862.624.904.174.073.624.553.19]t；i′f2＝20.31，短路电流增大，应投入电容支路，取xlm＝-0.2。分别以节点6、节点7、节点8、节点9、节点10、节点11、节点12、节点13、节点14为末端节点，重复步骤b1,b2得到节点2短路电流如下表所示：表1节点2路电流对比末端节点67891011121314短路电流if223.1812.8120.4433.9421.2820.8320.5820.8620.56由上表可知，if2的最小值if2min＝12.81，在末端节点为节点7时取得。由以上结果可知，通过支路追加法可以降低节点2的短路电流且最优支路为节点2至节点7支路。当前第1页12