电力系统暂态分析微代数方程的推导.docVIP

始资料： 3机9节点系统图： 图1 支路参数： 表1 首端 母线名末端 母线名电阻 （标幺值）电 抗 （标幺值）容纳之半 （标幺值）变压器 非标准变比450.0100.0850.088460.0170.0920.079570.0320.1610.153690.0390.1700.179780.00850.0720.0745890.01190.10080.1045140.00.05761.0270.00.06251.0390.00.05861.0 发电机参数 表2 发电机母线名TJRaXdX’dXqX’qT’d0T’q0D1147.280.00.14600.06080.09690.09698.960.02212.800.00.89580.11980.86450.19696.000.5350.0336.020.01.31250.18131.25780.25008.590.6000.0注：表中所有时间常数的单位为“S”，阻尼系数D及所有电阻、电抗均为“表幺值”。 正常运行情况下的系统潮流表： 表3 母线名电压发电机负荷幅值相角/度有功功率无功功率有功功率无功功率11.0400.00000.71640.270521.02509.28001.63000.066531.02504.66480.8500-0.108641.0258-2.216850.9956-3.98881.25000.500061.0127-3.68740.90000.300071.02583.719781.01590.72751.00000.350091.03241.9667 消去负荷节点的分析电路图： 1.消去系统的负荷节点 计算负荷的等值阻抗： 因为 所以 负荷对应的导纳为： 令用于系统潮流计算的节点导纳矩阵为。暂态稳定分析中，负荷采用恒阻抗模型，将相应节点的负荷导纳值并入系统中，即修改系统中5,6,8号节点的自导纳值为： 如此,便可得到消去系统负荷节点后的节点导纳矩阵. 2.系统的初值计算 系统暂态稳定分析过程中需要求解微分方程，因此需要计算系统初值。暂态稳定分析中，需用到系统运行状态变量和运行变量 。在本次系统暂态稳定分析中，发电机采用经典模型、负荷为恒阻抗模型、网络用节点导纳矩阵描述；因此， 据原始资料可知，已知，故需要据系统原始资料求解。 求各发电机暂态电动势： 电路图为如下： 由简单系统潮流计算公式得： ; 所以 用直角坐标可表示为： 其中 故可得 在暂态稳定分析计算过程中，认为系统发电机的调速系统可维持原动机的输入功率不变。即有 ， 3.暂态稳定分析的网络方程 暂态稳定计算过程中的网络代数方程与系统潮流计算时的网络代数方程有所区别，前者是由消去系统网络负荷节点的网络代数方程和发电机方程组合而成的。用于暂态稳定分析的代数方程可以描述为： （1） 其中为发电机的暂态电动势，为消去负荷节点和考虑发电机内节点，并消去网络中联络节点后的节点导纳矩阵，为发电机实际注入电流。 列写发电机方程： 暂态稳定计算中，考虑励磁系统具有维持发电机暂态电抗后的暂态电动势保持不变的能力，并令,即发电机采用经典模型。发电机用表示时，按固定在发电机转子上的d,q轴坐标来建立发电机电压平衡方程式有: 写成矩阵形式为： 系统潮流计算的代数方程中各节点的电压电流等相量均是以某一同步旋转轴的坐标系为基准的，今假设此坐标轴为直角坐标轴x-y。因此，为了便于发电机方程与网络方程的联立求解，需要将发电机d-q坐标系下的电压平衡方程转化为x-y坐标轴下的电压平衡方程。同时，在暂态稳定计算过程中，可以假设系统同步旋转轴和x轴重合，故在发电机经典模型下可以将发电机q轴与x轴间的夹角认为是发电机的绝对功率角。坐标间的转换关系，可参照如下相量图 。 由上图可以得知 故