电力系统计算报告

本文格式为Word版，下载可任意编辑——电力系统计算报告

电力系统计算报告

院（系）电气工程及自动化授课教师胡林献

姓名张远实学号13S106052

P-Q分解潮流法

简述P-Q分解潮流法的基本原理、计算过程、计算框图。用C语言编制P-Q分解法潮流程序，并用电科院6机22节点系统加以验证。要求采用稀疏技术、因子表技术和节点优化技术，并考虑负荷静态特性。

P-Q分解潮流法的基本原理：

P-Q分解法即是基于采用极坐标形式表示的牛顿法，其根据电力系统实际运行状态的线路参数R/X寻常很小的状况，对求解修正量的修正方程系数矩阵加以简化，使其变为常数阵(即所谓的等斜率)，且P、Q迭代解耦。这样可减少每次迭代的计算时间，提高计算速度，又不影响最终结果，因此是寻常选用的一种方法。但在低电压配电网中，当线路R/X比值很大时，可能出现不收敛状况，此时应考虑更换其它方法。

计算过程：

1、形成有功迭代和无功迭代的简化雅克比矩阵B和B2、给定PQ节点的初值和各节点电压相角初值

(k)3、作有功迭代?P/Vi(k)，解修i?Pis?Pi?Pis?Vi?Vj(Gijcos?ij?Bijsin?ij)，计算?Pij?1n'''正方程?i(k?1)??i(k)???i(k)，得各节点电压相角的修正值。

4、作无功迭代?Qi?Qis?Qi?Qis?Vi?Vj(Gijsin?ij?Bijcos?ij)，计算?Qi(k)/Vi(k)，解

j?1n修正方程Vi(k?1)?Vi(k)??Vi(k)，得各节点电压幅值的修正量。

5、返回第三步，继续迭代到满足要求为止。

计算框图流程：

输入原始数据形成矩阵B及B并进行三角分解'''设PQ节点电压初值，各节点电压相角初值制迭代计数k=0KP?1,KQ?1用公式计算不平衡功率?Pi(k)，计算?Pi(k)/Vi(k)Ymax{?Pi(k)}??P?置Kp?0解修正方程求?i?(k)KQ?0?Y?i(k?1)??i(k)???i(k)置KQ?1(k)用公式计算不平衡功率?Qi，计算?Qi(k)/Vi(k)max{?Qi(k)}??Q?Y置KQ?0解修正方程求?Vi(k)KP?0?YVi(k?1)?Vi(k)??Vi(k)置KP?1k?1?k计算平衡节点功率及全部线路输出算例描述

用电科院6机22节点算例进行验证，详细参数和网络拓扑图如下所示

支路号首末端节点号17-827-938-949-22512-13614-19716-19816-20916-211016-181119-211220-221321-22148-221511-12支路号首末端节点号17-129-2322-3419-4518-5617-6710-9811-10915-121017-131115-141216-17节点号12161112

表1线路参数

支路电阻支路电抗0.01060.07400.01470.10400.00340.01310.05590.21800.002450.02550.00340.02000.05780.21800.01630.06620.03740.17800.00330.03330.01140.03700.02140.08590.01500.06070.05370.19000.0033

0.0343

表2变压器支路数据

电阻电抗0.00.01500.00.02170.00.01240.00.06400.00.03750.00.03370.0-0.0020.00.01800.00.01800.00.01000.0-0.0020.0

0.0010

表3并联电容器数据

电纳-1.36650.50176-1.3665-1.3665

对地电纳/20.00.00.00.19541.3950.00.18870.23530.1640.00.00.30080.21980.16531.08797

变比1.0501.0751.1001.0251.0501.0001.0001.0001.0001.0001.0001.027

intisgn(inta,intb){

//****本函数功能返回值为a的绝对值b的符号****////参数1提供值，参数2提供符号//if(b0)a=-a;returna;}

intfind(intk,inta[],int*z){

//****本函数查找a[]中是否有fabs(k)有则返回0，无则返回1****////参数1为待查找量，参数2待探寻数组，参数3返回k在a[]中的次序号//inti;for(i=1;in