PQ分解法计算潮流
更新时间:2024-05-24 13:25:01 阅读量: 综合文库 文档下载
一、PQ分解法的原理
P-Q分解法是牛顿-拉夫逊法潮流计算的一种简化方法。 P-Q分解法利用了电力系统的一些特有的运行特性,对牛顿-拉夫逊法做了简化,以改进和提高计算速度。
的基本思想是根据电力系统实际运行特点:通常网络上的电抗远大于电阻,则系统母线电压幅值的微小变化对用功功率的改变影响很小。同样,母线电压相角的的改变对无功功率的影响较小。因此,节点功率方程在用极坐标形式表示时。它的修正方程式可简化为:
??P??H0???????Q???0L???UU? ??????将P、Q分开来迭代计算,因此大大地减少了计算工作量。但是H、L在迭代过程中仍
将不断变化,而且又都是不对称矩阵。对牛顿法的进一步简化。为把上式中的系数矩阵简化成迭代过程中不变的对称矩阵。
在一般情况下线路两端的电压相角?ij是不大的,因此可以认为:
cos?ij?1Gijsin?ij=Bij
Qi=Ui2Bii
考虑到上述关系,可以得到:
Hij?UiBijUjLij?UiBijUj节点的功率增量为:
n
?Pi?Pis?Ui?Uj(Gijcos?ij?Bijsin?ij)j?1?Qi?Qis?Ui?Uj(Gijsin?ij?Bijcos?ij)j?1n
P-Q分解法的特点:以一个n-1阶和一个n-m-1阶线性方程组代替原有的2n-m-1阶线性方程组;修正方程的系数矩阵B’和B”为对称常数矩阵,且在迭代过程中保持不变;P-Q分解法具有线性收敛特性,与牛顿-拉夫逊法相比,当收敛到同样的精度时需要的迭代次数较多。
二、程序说明
1.数据说明
Branch1.txt:支路参数矩阵
第1列为支路的首端编号;第2列为支路的末端编号(首端编号小于末端编号);第3列为之路的阻抗;第4为支路的对地容抗;第5列为支路的变比;第6列为折算到那一侧的标志 Branch2.txt:节点参数矩阵
第1列为节点所接发电机的功率;第2列为节点负荷的功率;第3列为节点电压的初始值;第4列为PV节点的电压V给定值;第5列为节点所接的无功补偿设备的容量;第6列为节点分类标号igl,其中igl=1为平衡节点,igl=2为PQ节点,igl=3为PV节点。 n为节点数、sb为平很母线节点号(固定为1)、pr为误差精度。
2.运行
将上述三个文档及主程序文件放在同一个路径下,程序运行时,按照相应的提示输入如下数据:n为节点数、sb为平很母线节点号、pr为误差精度,运行,得出结果。
三、流程图
开 始 输入原始数据 形成矩阵B’和B” 设置PQ节点电压初值,各节点电压相角的初值 置迭代计数n=0 Kp=1,Kq=1 计算不平衡有功功率 是 max{?Pik}??p?是 解修正方程,并修正后电压相角 置Kp=0 Kq=0? 是 否 置Kq=1 计算不平衡无功功率 无功收敛? 否 解修正方程,求得修正后电压 是 置Kq=0 是 Kp=0? 否 置Kp=1,k=k+1 计算平衡节点功率及全部线路功率 输出
四、源程序及运行结果 x = size(B2); n = x(1);
m = length(find(B2(:,1)==1));
y = size(B1); G = zeros(n); B = zeros(n); for i=1:y(1)
QS = B1(i,1); %起始节点 ZZ = B1(i,2); %终止节点 R = B1(i,3); X =B1(i,4); BB = B1(i,5);
K = 1/B1(i,6); %变比 RX2 = R^2+X^2;
G(QS,ZZ) = -K*R/RX2; G(ZZ,QS) = -K*R/RX2; B(QS,ZZ) = K*X/RX2; B(ZZ,QS) = K*X/RX2;
G(QS,QS) = G(QS,QS) + R/RX2; B(QS,QS) = B(QS,QS) - X/RX2; B(QS,QS) = B(QS,QS) + BB/2;
G(ZZ,ZZ) = G(ZZ,ZZ) + K^2*R/RX2; B(ZZ,ZZ) = B(ZZ,ZZ) - K^2*X/RX2; B(ZZ,ZZ) = B(ZZ,ZZ) + BB/2; end;
for i=1:n
B(i,i) = B(i,i) +B2(i,8); end;
Delta_P = zeros(n,1); Delta_Q = zeros(n,1); Delta_V_Ample = zeros(n,1); Delta_V_Angle = zeros(n,1);
PP = find(B2(:,1)<3); % PQ/PV 节点 QQ = find(B2(:,1)==1); % PQ 节点
%V_Ample = B2(:,2);
%V_Angle = B2(:,3)*pi/180; V_Ample = B2(:,2); V_Ample(QQ) = 1.0;
V_Angle = zeros(n,1); Pis = B2(:,4)/100; Qis = B2(:,5)/100; Pld = B2(:,6)/100; Qld = B2(:,7)/100;
Pis(PP) = Pis(PP) - Pld(PP); Qis(QQ) = Qis(QQ) - Qld(QQ);
kp = 0; kq = 0;
ep = 1e-7; %误差 k = 0;
while(~(kp*kq)) k =k + 1;
%----------------------------------------------------------------------
Delta_V_Angle_ij = V_Angle*ones(1,n) - ( V_Angle*ones(1,n) )'; % 计算相角差矩阵
S = G.*cos(Delta_V_Angle_ij) + B.*sin(Delta_V_Angle_ij); % Sij = Gij *cos(delta_ij) + Bij *sin(delta_ij)
P = V_Ample.*(S*V_Ample); % P为计算出来的母线功率
Delta_P(PP) = Pis(PP) - P(PP); % 有功不平衡量
%[Pis(PP), P(PP),Pis(PP) - P(PP),(Pis(PP) - P(PP))./Pis(PP)*100]
if max(abs(Delta_P)) if kp==0 % 有功不平衡量大于允许范围就作修正 Delta_V_Angle(PP) = -(B(PP,PP)\\( Delta_P(PP)./V_Ample(PP) ))./V_Ample(PP); % 解修正方程 V_Angle(PP) = V_Angle(PP) + Delta_V_Angle(PP); % 修正相角 kq = 0; % 这是必须的,认为无功不平衡量还是不满足 end;
正在阅读:
PQ分解法计算潮流05-24
石首市文峰中学2018年质检二考试复习训练(03)09-24
双语生活与梦想03-07
乡镇开展一村一策年度工作总结报告08-04
基于单片机的身高体重测量仪 - 图文04-14
醉在诗歌之上的现代诗03-21
qq空间自我介绍模板04-16
诚信和企业文化之间有着怎样的关系05-18
难忘的中秋节作文300字精选06-14
财务会计实习生的报告范文02-22
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 解法
- 潮流
- 计算
- 临时污水处理站基坑支护
- 用于银行贷款2013年衣领净洗涤剂项目可行性研究报告(甲级资质+
- 2013版 普通试卷答案及评分标准建设法规教程A
- 人教版初中体育与健康—跨越式跳高(助跑起跳)教学设计 - 图文
- 系统科学与系统思维
- 基于单片机可编程作息时间控制器设计与制作
- 政府和社会资本合作(PPP)-城市道路及隧道大修工程项目可行性研究
- 生产实习总结表 - 图文
- 2011年普通高中信息技术会考笔试考点解析
- 排水管网课设任务书(2010) - 图文
- 婚姻财产关系的思考与论述
- 《信息管理》教学大纲
- 注会讲义《财管》第二章财务报表分析02
- ASP程序设计课程设计报告 毕业设计论文
- 安全生产重大事故应急救援预案
- 中国电信 FTTH 终端设备技术要求( 1FE+1POTS SFU )
- 《教育心理学》汇总
- 新所得税准则资产负债表债务法探析(1)
- 山师英语阅读二试题D及答案
- 国家林业局关于印发《中央财政造林补贴试点检查验收管理办法(试