电力系统稳态实验报告..docx

电力系统稳态潮流计算上机实验报告一、问题如下图所示的电力系统网络，分别用牛顿拉夫逊法、PQ解耦法、高斯赛德尔法、保留非线性法计算该电力系统的潮流。该电力系统的电路参数如下，名称电阻（pu）电抗(pu)-B/2(pu)线路10.040.25-0.25线路20.10.350线路30.080.30-0.25变压器参数如下，名称电阻（pu）电抗(pu)变比变压器100.031.05:1变压器200.0151:1.05发电机的参数如下，名称电压（pu）相角（rad）有功(pu)发电机11.050*发电机21.05*5*表示任意值负荷参数如下，名称有功(pu)无功(pu)负荷11.60.8负荷22.01.0负荷33.71.3二、问题分析如上图所示的电力系统，可以看出，节点1、2、3是PQ节点，节点4是PV节点，而将节点5作为平衡节点。根据问题所需，采用牛顿拉夫逊法、PQ解耦法、高斯赛德尔法、保留非线性法，通过对每次修正量的收敛判据的判断，得出整个电力系统的潮流，并分析这四种方法的收敛速度等等。算法分析1.牛顿拉夫逊法节点5为平衡节点，不参加整个的迭代过程，节点1、2、3为PQ节点，节点4为PV节点，计算修正方程中各量，进而得到修正量，判断修正量是否收敛，如果不收敛，迭代继续，如果收敛，算出PQ节点的电压幅值以及电压相角，得出PV节点的无功量以及电压相角，得出平衡节点的输出功率。潮流方程的直角坐标形式，直角坐标形式的修正方程式，修正方程式中的各量值的计算，Jacobi矩阵的元素计算，牛顿拉夫逊法潮流计算的流程图如下，2.PQ解耦法如同牛顿拉夫逊法，快速解耦法的前提是，输电线路的阻抗要比电阻大得多，并且输电线路两端的电压相角相差不大，此时可利用PQ快速解耦法，来计算整个电力系统网络的潮流。快速解耦法的迭代方程组，?P＝－H?θ?Q＝－L(?U/U)快速解耦法潮流计算的流程图如下，3.高斯赛德尔法高斯赛德尔法原理较前两种方法简单，程序设计十分容易，占内存小，是所有的潮流计算方法中迭代计算量最小的。高斯赛德尔法的迭代格式为，高斯赛德尔法的收敛判据如下，在高斯赛德尔法中，不应对PV节点的幅值进行修正，只对其电压相角进行一定的修正。高斯赛德尔法潮流计算的流程图为，4.保留非线性法保留非线性法主要是在牛顿拉夫逊法的基础上，通过泰勒展开，保留到二阶项，由于三阶导数值等于零，所以泰勒展开式是准确的，无截断误差，与牛顿拉夫逊法不同的是，保留非线性法只需算一次jacobi矩阵，每次迭代得到的修正量都是在初始值上的修正量，因此，保留非线性法的计算量小于牛顿拉夫逊法的计算量，大大节约计算机的内存空间，提高计算机的计算速度。保留非线性的迭代格式为，式中，k表示迭代次数；J为按x＝x(0)估计而得。收敛判据为，也可采用相继二次迭代的二阶项之差作为收敛判据(更合理)，相应的收敛判据如下，保留非线性法的流程图如下，三、MATLAB仿真结果1.牛顿拉夫逊法迭代次数k=6；各节点的电压值、有功功率以及无功功率见下表。名称电压幅值电压相角电压向量有功功率无功功率节点10.8683-0.08290.8653 - 0.0719i-1.6-0.8节点21.07830.31081.0267 + 0.3298i-2-1.0节点31.0370-0.07461.0341 - 0.0773i-3.7-1.3节点41.05000.38040.9749 + 0.3899i51.7857节点51.050001.05002.57602.28032.PQ解耦法迭代次数k=13；各节点的电压值、有功功率以及无功功率见下表。名称电压幅值电压相角电压向量有功功率无功功率节点10.8683-0.08290.8653 - 0.0719i-1.6-0.8节点21.07830.31081.0267 + 0.3298i-2-1.0节点31.0370-0.07461.0341 - 0.0773i-3.7-1.3节点41.05000.38040.9749 + 0.3899i51.7857节点51.050001.05002.57602.28033.高斯赛德尔法迭代次数k=137;各节点的电压值、有功功率以及无功功率见下表。名称电压幅值电压相角电压向量有功功率无功功率节点10.8688-0.08470.8657 - 0.0735i-1.6-0.8节点21.07840.30771.0277 + 0.3266i-2-1.0节点31.0372-0.07491.0343 - 0.0776i-3.7-1.3节点41.05000.37720.9762 + 0.3868i51.7857节点51.050001.05002.57602.28034.保留非线性法 迭代次数k=11。各节点的电压值、有功功率以及无功功率见下表。名称电压幅值电压相角电