110KV电力网最大负荷时的潮流计算说明书

1、110KV电力网最大负荷时的潮流计算 专业名称；电气工程及其自动化 学号： 学生姓名: 指导教师: 班级： 电力系统的潮流计算在电力系统稳态分析和电力系统设讣中有很重要的作用，潮 流计算也是电力系统暂态分析的基础，潮流计算是根据给定的系统运行条件，来计算 系统各个部分的运行状况，主要包括电压和功率的计算。利用电子计算机进行电力系 统的潮流计算的算法，已经出现了很多，其中应用最为广泛的是基于牛顿拉夫逊法的 潮流计算方法。 本设计采用牛顿拉夫逊法的潮流计算方法，首先对网络节点进行划分和编号，建 立电力网络的数学模型即电力系统的网络方程式。将电力网络进行等值电路的变换 （主要包括输电线路和变压器的等

2、值电路的变换），利用计算，计算出各个节点的导 纳，形成节点导纳矩阵，建立雅克比矩阵。最后对雅克比矩阵进行多次迭代，求出迭 代后的PV、PQ以及平衡节点的功率和电压的计算值，并对值进行精确修正，与实际 的值进行比较。开发工具釆用Mat lab编程语言进行软件的编写，本文采用一个6节 点的网络进行分析。 关键词：潮流计算：节点导纳矩阵：牛顿拉夫逊：Matlab IV Abstract Power flow calculation has a very important role in power system steady-state analysis and power system desi

3、gn, and it is also the basis of transient analysis in power system. Flow calculation is based on given conditions of the power system and calculates the operational status of every part of the system including voltage and power So far there are kinds of algorithm which use the electronic coinpufer i

4、n power flow calculation, the most widely used algorithm is the Newton-Raphson power flow calculation method. This design uses the Newton Raphson power flow calculation method, the first network node and divide numbers networking equation mathematical model of the electricity network that is the pow

5、er system. Will transform the electricity network equivalent circuit (including conversion equivalent circuit of transmission lines and transformers), using the calculated calculated admittance of each node, the node admittance matrix formation, the establishment of Jacobian matrix. Finally, multipl

6、e iterations Jacobian matrix obtained calculated value PV.PQ and the balance of power and voltage node iteration, the value of accurate and corrected, with the actual value, Deve!opinent tools using Matlab programming language to write software, use a 6-node network for analysis. Key Words： flow cal

7、culation; node admittance matrix; Newton-Rap son Gauss seidel; Matiab 59 外文原文 60 外文译文 66 1绪论 1.1概述 电力是衡量一个国家经济发展的主要指标也是反映人民生活水平的重要标志，它 己成为现代工农业生产交通运输以及城乡生活等许多方面不可或缺的能源和动力电 力系统，是山发电输电变电配电和用电等环节组成的电能生产与消费系统，它的功能 是将自然界的一次能源通过发电动力装置转化成电能再经输电变电和配电，将电能供 应到各用户为实现这一功能电力系统在各个环节和不同层次，还具有相应的信息与控 制系统对电能的生产过程进行测

8、量调节控制保护通信和调度，以保证用户获得安全经 济优质的电能。 电力系统的出现使电能得到广泛应用推动了社会生产各个领域的变化，开创了电 力时代出现了近代史上的第二次技术革。20世纪以来电力系统的发展使动力资源得 到更充分的开发，工业布局也更为合理使电能的应用不仅深刻地影响着社会物质生产 的各个侧面，也越来越广地渗透到人类日常生活的各个层面，电力系统的发展程度和 技术水准已成为各国经济发展水平的标志之一。 MATLAB自1980年问世以来它的强大的矩阵处理功能，给电力系统的分析计算带 来许多方便，在处理潮流计算时其计算机软件的速度已无法满足大电网模拟和实时控 制的仿真要求。而高效的潮流问题相关软

9、件的研究己成为大规模电力系统仿真计算的 关键，随着计算机技术的不断发展和成熟对MATLAB潮流计算的研究，为快速详细地 解决大电网的计算问题开辟了新思路。 1.2选题背景 电力系统已经与我们的生活息息相关不可分割，进行电力系统潮流计算是保证电 力系统正常运行的必要，计算具体来讲电力系统潮流计算具有以下意义。在电网规划 阶段通过潮流计算合理规划电源容量及接入点合理规划网架选择，无功补偿方案满足 规划水平的大小方式下潮流交换控制调峰调相调压的要求，在编制年运行方式时，在 预计负荷增氏及新设备投运基础上选择典型方式进行潮流计算，发现电网中薄弱环节 供调度员日常调度控制参考，并对规划基建部门提出改进网

10、架结构加快基建进度的建 议，正常检修及下的潮流计算，用于日运行方式的编制指导发电厂开机方式有功无功 调整方案及负荷调整方案，满足线路变压器热稳定要求及电压质量要求。 预想事故设备退出运行对静态安全的影响分析，及作出预想的运行方式调整方案 总结，为在和规划方案的研究中都需要进行潮流计算，以比较运行方式或规划供电方 案的可行性，可靠性和经济性同时为了电力系统的运行状态也需要进行大量而快速的 潮流计算。因此，潮流计算是电力系统中应用最广泛最基本和最重要的一种电气运算, 利用电子计算机进行潮流计算从20世纪50年代中期就己经开始此后潮流计算，曾采 用了各种不同的方法这些方法的发展主要是H绕着对潮流汁算

11、的一些基本要求进行 的电力系统潮流计算，属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。 因此，其数学模型不包含微分方程是一组高阶非线性方程非线性代数方程组的解法， 离不开迭代因此潮流计算方法。首先要求它是能可靠的收敛并给出正确答案，随着电 力系统规模的不断扩大，潮流问题的方程式阶数越来越高，U前已达到儿千阶，祺至 上万阶对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的，这种 悄况，促使电力系统的研究人员不断寻求新的更可靠的计算方法潮流算法的研究，仍 然非常活跃但是大多数研究，都是H绕改进牛顿法和P-Q分解法进行的。此外，随着 人工智能理论的发展遗传算法模糊算法也逐渐被引入

12、潮流计算，但是到U前为止这些 新的模型和算法还不能取代牛顿法和P-Q分解法的地位，山于电力系统规模的不断扩 大对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的 应用，成为重要的研究领域。 1.3国内外发展现状 电力系统潮流计算是研究电力系统稳定运行悄况的一种计算，它根据给定的运行 条件及系统的接线悄况确定整个电力系统各部分的运行状态，即各母线的电压，各元 件中通过的功率，系统的功率损耗等等。在电力系统运行方式和规划方案的研究中， 都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济 性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。通

13、过潮流计 算，不仅能够确定整个系统的运行状态，而且为其他电力系统分析软件提供初值条件, 诸如配电重构、状态佔计、故障处理、系统无功优化、设计规划等。 电力系统潮流计算经历了儿个阶段：开始是手算，只能求解简单系统的潮流而且 还需要采用一些简化假设。其后曾广泛采用计算台（包括直流计算台和交流计算台） 解算潮流，这类实体仿真方法精度差，规模小。上世纪50年代电子计算机发明以来, 借助计算机为解决潮流计算问题提供了新的技术手段。早期利用计算机求解潮流计算 方法是以导纳矩阵为基础的简单迭代法，称为高斯迭代法。这种方法原理简单，内存 需求较少，但算法收敛性极差。后来发展了以阻抗矩阵为基础的算法。这种方法收

14、敛 性好，但内存占用量大大增加，限制了解题规模。牛顿-拉夫逊(Xewton-Rsphson即 X-R)方法是解非线性代数方程组的一种基本方法，在潮流计算中也得到了十分广泛 的应用。60年代中后期，牛顿-拉夫逊潮流算法采用了 Tinney提出的稀疏矩阵技术 和节点优化编号技术，使牛顿-拉夫逊法成为电力系统潮流计算中广泛采用的优秀算 法，而且至今它仍是潮流汁算中的一种广泛使用的基本算法。70年代中期,Stott在 大量进行潮流计算实践的基础上，提出了快速分解法2的潮流计算模式，使得潮流 汁算的速度大大提高，并可以应用于在线。80年代末期，对快速分解算法潮流计算 的收敛机理给出了比较满意的解释。至今

15、，电力科技工作者对潮流计算投入了大量的 研究出现了大量的研究成果。同时，新的潮流模型，如约束潮流、连续潮流、动态潮 流、随机潮流、最优潮流、开断潮流等，也进行了大量的研究。 利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后，潮流计 算曾采用了各种不同的方法，这些方法的发展主要是H绕着对潮流计算的一些基本要 求进行的。对潮流计算的要求可以归纳为下面儿点： (1) 算法的可靠性或收敛性 A升圧站T-1) 变电騎2 (T2) 变电所3 (T3) a号 无ttSMS变压 63000 有 SMS 交 S8 31S00 WnM压变压ff 20000 定0* (MVA) 90 3LS 20

16、台ft 2 2 2 空耗(kW) S5.8 3 但在国外，有名制的使用也很普遍。至于变圧器模型的使用范ffl,则泾消分明。手算 时，都是用形或T型等值电路模型；计算机计算时，都是用等值变压器或n型等值 电路模型。 此外，在制定电力网络等值电路模型时，有时还同时作某些简化，常见的有: 线路的电导通常都被略去；变压器的电导有时以具有定值的有功功率损耗的形式 出现在电路中；lOOkm以下架空线路的电纳被略去；100300km架空线路或变压 器的电纳有时以具有定值的容性或感性无功功率损耗的形式出现在电路中。有时，整 个元件，其至部分系统都可能不包括在等值电路中。例如，将某些发电厂的高压母线 看作为可维

17、持给定电压、输出给定功率的等值电源时，这些发电厂内部的元件就不再 包括在等值电路中。 2. 4节点导纳矩阵 在电路原理课程中，已导出了运用节点导纳矩阵的节点电压方程 B B B(2.11) 上式中，厶是节点注入电流的列向量，可理解为某个节点的电源电流与负荷电 流之和，并规定电源流向网络的注入电流为正。因此，仅有负荷的负荷节点注入电流 就具有负值。是节点电压的列向量。因通常以大地作参考节点，网络中有接地支 路时，节点电压通常就指该节点的对地电压；网络中没有接地支路时，各节点电压可 指各该节点与某一个被选定参考节点之间的电压差。坨是一个节点导纳矩阵，它的 阶数n等于网络中除参考节点外的节点数。 它

18、可展开为 A 旷 A 厶 * 5” * 上1 Y nn * (2.12) 2.4.1节点导纳矩阵的形成 根据定义直接求取节点导纳矩阵时，注意以下儿点： (1) 节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考 节点一般取大地，编号为零。 (2) 节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应 节点所连接的不接地支路数。 (3) 节点导纳矩阵的对角元素就等于该节点所连接导纳的总和。因此，与没 有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。 (4) 节点导纳矩阵的非对角元素等于连接节点i, j支路导纳的负值。因此， 在一般悄况下，节点导纳矩阵的对

19、角元素往往大于非对角元素的负值。 (5) 节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而， 一般只要求取这个矩阵的上三角或下三角部分。 (6) 网络中的变压器。 242节点导纳矩阵的修改 (1) 从原有网络引出一支路，同时增加一节点。 设i为原有网络中的节点，j为新增加的节点，新增加支路导纳为必严 则因新增 一节点，节点导纳矩阵将增加一阶。 新增的对角元岭，由于在节点j上只有一个支路儿，将为二;新增的非 对角元Y, = Yji = -yij ；原有矩阵中的对角元将增加儿，儿二为。 (2) 在原有网络的节点i、j之间增加一支路。 这时山于仅增加支路不增加节点，节点导纳矩阵阶数不变，

20、但与节点i、j有关的 元素应作一下修改，其增量为： 儿二儿二y沪七二二一为 (3) 在原有网络的节点i, j之间切除一支路。 切除一导纳为的支路，相当于增加一导纳为一必7的支路，从而与节点i、j 有关的元素应作如下修改： Xi = -yijAYjj =-儿,岭=岭=Xj (4) 原有网络的节点i、j之间的导纳由打改变为y- O 这种悄况相当于切除一导纳为的支路，并增加一导纳为x/的新支路。从而 与节点i、j有关的元素应作如下修改： =叫=为-为 (5) 原有网络节点i、j之间变压器的变比山化改变为 这种悄况相当于在i、j节点之间并联一个变比为-k*的变压器，再并联一个变比 为的变压器，即相当于

21、修改变压器。修改前，i、j节点之间的自导纳和互导纳为: 修改后，引用理想变压器”的n型等值电路，变压器变比山改变为时, 原网中与节点i、j有关的元素应作如下修改： 乂 =o，M = 2.5潮流计算节点的类型 用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内 的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力 系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般 是给出发电机母线上发电机的有功功率P和母线电压的幅值U,给出负荷母线上负荷 消耗的有功功率P和无功功率Q。主要U的是由这些已知量去求电力系统内的各种电 气量。所以，根

22、据电力系统中各节点性质的不同，很自然地把节点分成三类： (l) PQ节点 对这类节点，等值负荷功率/in、Qu和等值电源功率耳i、Qi是给定的，从而注 入功率R、Qj是给定的，待求的则是节点电压的大小E和相位角q。属于这一类节点的有按给定有功无功功率发电的发电厂母线和没有其他电源的变电所母线。 (2) PV节点 对这类节点，等值负荷和等值电源的有功功率3。再将求得的 牢代入，乂可以求得纣、丿的新值，从而解得心，以及）=甲+(q -B)+q(G也 +B再) /-I (3. 13a) (3. 13b) (3. 13c) 占=(；+/?) 当7知时，山于对特定的j,只有该特定节点的fj和勺是变量，山

23、式(3.12)、 式(3. 13)可得: 血=-举+G/； Mj= 二 q+B-Z- (3. 14a) cQ- dQ ,沁uf A c 沁U： A Ri= =0:S= =0 当j = i时，为使这些偏导数的表示武更为简洁，先引入节点诸如电流的表示式 如下: B也)曲泸川 ?! 0 73 然后山式(3-12)、式(3-13)和上式可得 Gijfi - Bq 筋 iqj - Bij九+G代 B 2. nl=6（支路数）; 3. 输入修正值：ip=O. 00001： 电源侧的功率 负荷侧的功率 该点的电压值 节点的类型 0 0 1 2 表4-3 由各支路参数形成的矩阵 0 0 1 2 0 0 1 2

24、 0 3. 7+jl. 3 1 2 0 5+jO 1.05 3 0 0 1.05 1 四、修正方程 计算节点1的不平衡量和AQ， 甲）凯-皆= SAV； J = 以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对PQ节点来说, 耳和2是给定的，W而可以写出 、P产P工IG启-Bjf ）f Xgj眉厂 。广 2-丿+勺工（Gj+B“）=o 对PV节点来说，给定量是只才【乂J因此可以列出 当时，雅可比矩阵中非对角元素为 云一# 一+吋） 0； 迢 八 泌=泌=0 当八：时,雅可比矩阵中对角元素为: - （G内-BJ）- G启一 BJ. % 1 PAP V2（久人 + Be） - G/ + B尼

25、Ofi Jt-I =若（GiJk + Be - G-J. + B启 6边=-9凶-色门+ % + % M 陆 ， 代入数值后的修正方程为: -1.25 -5,5 0.5 3 e：W -2 -5.5 1.25 3 -0.5 -1 0.5 3 -1.3 -7 0.5 _ 0 0 -2 0 _ 0 _ 求解修正方程得: -0.2547 -03611 0 Ml -0J015 5 6节点潮流计算程序 x =xlsread(* sj. xls*, A3 ：A3) y=xlsread(* sj. xls, B3 ：B3) z=xlsread(* sj. xlsC3：C3) pr=xlsread(* sj.

26、xls, B4 ：B4) zhilu =xlsread(* sj. xls * E4:K100) point =xlsread(* sj. xls, M4:V100) duidi=xlsread(* sj. xls*, * X4: YlOO) n=x; nl=y; ph=z; Bl=zhilu; B2=point; A=duidi; Tl=zeros(n, 2); T2=zeros(nb 3); i=sqrt(-l); for j=l:n Tl(j. l)=Bl(j.3)+Bl(j.4)*i; end for j=l:nl T2(j, l)=B2(j,2)+B2(j,3)*i; T2(j,2)

27、=B2(j,4)+B2(j,5)*i; T2(j,3)=B2(j,6)+B2(j,7)*i; end Bl=zeros(n, 6); B2=zeros(nb 6); for j=l:n Bl(j. l)=zhilu(j. 1); end for B2(j. l)=T2(j. 1)； B2(j,2)=T2(j,2); B2(j,3)=T2(j,3); B2(j,4)=point(j,8); B2(j.5)=point(j.9); B2(j.6)=point(j. 10); end %以上为从excel中导入初值的程序 Y=zeros(n) ; e=zeros(l n) ;f=zeros(l n)

28、 ;V=zeros(l n); 0=zeros(l,n): Sl=zeros(nl); for i=l:n if A(i,2)=0; P=A(ia)； Y(p,p)=l/A(i,2); end end if Bl(i,6)=0 p=Bl(i. l)；q=Bl(i.2)： else p=Bl(i,2) ;q=Bl(i, 1); end Y(p, q)=Y(p, q)-l/(Bl(i,3)*Bl(i,5); Y(q, p)=Y(p. q); Y(q, q)=Y(q, q)+1./(Bl (i, 3)*B1 (i, 5厂2)+B1 (i, 4)/2; Y(p, p)=Y(p. p)+l/Bl(i,3

29、)+Bl(i,4)/2; end %求导纳矩阵 dispC导纳矩阵Yi); disp(Y); G=real(Y);B=imag(Y); for i=l:n e(i)=real(B2(i.3); end for i=l:n S(i)=B2(i, 1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5) ;%对角线上节点对地导纳(包含各节点对地参数) end %求各节点功率偏移量 P=real(S);Q=imag(S); k=0;IT=l;a=0;%k一一迭代次数；IT 一一没有达到精度要求个数 while IT=0 IT=0;a=a+l;num2=0:numl=0; for i=l:n

30、 if i=ph C(i)=0; D(i)=O; for j=l:n C(i)= C(i)+G(i, j)*e(j)-B(i, j)*f (j); D(i)= D(i)+G(i, j)*f (j)+B(i, j)*e(j); end Pl=C(i)*e(i)+f(i)*D(i); Ql=f(i)*C(i)-D(i)*e(i); end %针对PQ节点求雅克比矩阵各参数H, N, J,L if B2(i,6)=l %1平衡节点,2PQ节点,3PV节点 DP=P(i)-Pl; DY(numb 1)=DP; %各节点功率偏移量形成的矩阵: 10*1.先从 PQ DQ=Q(i)-Ql; %的有功功率偏

31、移量开始 numl=numl+l; DY (numb 1)=DQ; dispC各节点功率不平衡量:DY=); disp(DY); for j=l:n if j=i X2=-G(i, j)*e(i)-B(i, j)*f(i); %X2=N(i, j) X3=-X2; %X3=J(i. j)=-X(i, j)=-X2 X4=X1; %X4=L(i, j)=H(i, j)=Xl J(P, q)=Xl; m=q+l; p=p+l; J(P q)=X3; q=q+l； J(P q)=X4; elseif j=i%Xl=H(i i) X2=-C(i)-G(i, i)*e(i)-B(i,;%X2=N(i,

32、i) X3=-C(i)+G(i, i)*e(i)+B(i, i)*f(i) ;%X3=J(i. i) X4=D(i)+B(i, i)*e(i)-G(i,; %X4=L(i, i) J(P q)=Xl; m=q+l; P=P+1； J(P q)=X3; q=q+l； J(P, q)=X4; end end elseif B2(i,6)=l numl=2*i-l; DY (numb 1)=DP; DV=V(i 厂 2-Vl; numl=numl+l; DY (numb 1)=DV; for jHlm xihg? j)羞氏(i; X5HP X6HP J (9 q)HX5- nFP+l _ J(n q

33、) qHq+u J(9q)HX9J(n q)HXW elseif jHisfph sihh? j) 娄2hn? j) 费 x5HR(i9j) 决 x6HS(i 二) X5v2*e(i)- %x5HR(i9 i) X6V2沃f(i)- J(pq)HX= nFP+1 _ J(n q)HXU qHq 十一 - J(p q)HXWJ(目 q)HX0 end end end end 73 disp(J)： %求解修正方程，即各节点电压修正量DX=Df, De,DY=J*DX -DX=inv (J) *DY DX=JDY; dispC各节点电压修正量:DX=) disp(DX) %求各节点电压新值 for

34、 i=l:(n-1) 73 p=2*i-l; Df=DX(p); q=2*i; De=DX(q); e(i)=e(i)-De; end for i=2*(n-l); DETl(i)=abs(DX(i); 辻 DETKiXpr IT=0; else %判断是否满足精度误差，如满足则跳出，否则返回继续迭代 IT=IT+1; end end k=k+l; %计算迭代次数 if k100 break end end dispC迭代次数：k=) disp(k) dispC没有达到精度要求的个数：IT=） disp(IT) %求各节点实际电压和功率 for i=l:n V(i)=sqrt(e(i 厂 2+

35、f(i 厂 2); shita(i)=atan(f(i). /e(i)*180/pi; j=sqrt(-l); end dispC各节点的实际电压标么值I；为（节点号从小到大排列）：） disp(U) dispC各节点的电压大小V为（节点号从小到大排列）：） disp(V) dispC各节点的电压相角shita为（节点号从小到大排列）：） disp(shita) for p=l:n C(p)=0; for q=l：n C (p) =C (p) +conj (Y (p, q) *conj (U (q); end S (p) =U (p) *C (p); end dispC各节点的功率S为（节点号

36、从小到大排列）：） disp(S) dispC各条支路的首端功率Si为（顺序同您输入Bl时一样）：） for i=l:nl if Bl(i,6)=0 p=Bl(i. l)；q=Bl(i.2); else p=Bl(i,2) ;q=Bl(i, 1); end Si (p, q) =U(p)*(conj (U(p)*conj (Bl (i, 4). /2) + (conj (U(p) *B1 (i, 5) -conj (U(q) )*conj(l./(Bl(i,3)*Bl(i,5); disp(Si (p, q); end disp （各条支路的末端功率Sj为（顺序同您输入B1时一样）：） if

37、Bl(i,6)=0 p=Bl(i. l)；q=Bl(i.2); else p=Bl(i2) ;q=Bl(i 1); end Sj (q, p) =U(q) * (conj (U(q) *conj (Bl (i, 4) /2) + (conj (U(q). /Bl (i, 5) -conj (U (p )*conj(l./(Bl(i,3)*Bl(i,5); disp(Sj(q, p); end dispC各条支路的功率损耗DS为（顺序同您输入Bl时一样）：） for i=l:nl if Bl(i,6)=0 p=Bl(i. l)：q=Bl(i.2); else p=Bl(i,2) ;q=Bl(i,

38、 1); end DS (i) =Si (p, q) +Sj (q, p); disp(DS (i) end fid=fopen(* jisuanjieguo. txt, wt); fprintf(fid, * 121408131*n*); fprintf (fid, * 潮流计算输出结果 *n); fprintf(fid/ 迭代次数 k 为：%d n,k); fprintf (fid, fprintf （fid,各节点的实际电压标么值U为（节点号从小到大排列）:n） for m=l:n fprintf(fid, 第 %d 个 实 际 电 压 标 么 值 U: %f+i*(%f) n , m,

39、 real (U(m), imag(U(m); end fprintf(fid, fprintf （fid,各节点的电压大小V为（节点号从小到大排列）:n ） for m=l:n fprintf （fid,第d个节点的电压大小V为（节点号从小到大排 列)：%fn , m, V(in); end fprintf (fid, fprintf （fid,各节点的电压相角时shita为（节点号从小到大排列）:n ） for m=l:n fprintf （fid,第d个节点的电压相角时shita为（节点号从小到大排 列):%fn*, m, shita(in); end fprintf (fid, fpr

40、intf （fid,各节点的功率S为（节点号从小到大排列）：n ） for m=l:n fprintf (fid,第d 个节点复功率：n*, m, real (S (tn), imag(S (m); end n ); fprintf (fid, fprintf （fid,节点有功功率P为：n）; for m=l:n fprintf (fid,第%ug FoTdei 2012-5-17 d：丄1 2012-5-15 3： 二jg 翁 Fold 柳 UGUI FOTdo*KOf iioiipajd JO tw: aip na:iartltnh miicKjo iwwinnsdai prti!；(i

41、|d dtp ot pd农抑 81 ftsigdq duydp 01 Sa s ”卿 stqi o|tq 0us)d 即 oi J04i indino imtpnixud m|1 m oSruinp uMxy pox pop isoa 押();笊 v9HH!)n 9tp Xq p却!gid 8!3$ 刘伽 3q)0)dpojd ue) wiiddns 災八郴 x|i imp 5npj, pnppn弓 U *定辆ui pppop 何)Muiod JO XiMidds Mdu oqi ui luojd 屛 lueuoduu XyetnxuBd s| sxxuoisrd pui w h03 m1 心

42、!pnO 环0占 9塀 jo luiqoMl 沁 iKOqv xj为/帔 7 doy Ip QtpTK肝 MJ9 npru,! y 4; zp/V 07必 jp 8询叙gd HpUXUffZ Q 冲gq J 5血与 Q suawoisno nviaisfiaNi ki isoo AxnvnO h3aoi no AaAuns v f 第 I - fewer Ottiny(nrcr of short* inlcmiptions and voltage 呻.the actual nei value of the dama was csiimaied in I 8 Million of USD 1,2

43、 can also become in the ftiUirc 1/3 with Lbv expected increase of sensitive loach. In tile presence of industrial process susceptible to voltage sags (as paper-mill, glass works, mkro-clUvnic products, etc J the pay Mck lime of a DVR is estimated in a period varying from 3 10 4.5 years 121- If the P

44、Q problems are due to harmonic or intcrhannoflic components in the industrial system currents and voltages, the compensation devices that can be inserted in the iKlwork to reduce the disnirhance. are passive or active fillers (14 J5(. To 薊 ve an example of the technical-economical cvaluaiimi of this

45、 gludoft 4 some customenk selected on the besis of the renmied qucsuonnarc and of indicaUon given by (he AssociiMjon, have been further!) contacted in order lo deepen the invarigation by means of personal interviews. The questionnaire has bom structured m sovcral sodiwi, each aimed at inquirwg infbn

46、naliun regarding the (bllowmg data rekvarn lo the individual customer: - data for identifyu ihe type of activity, the number of enployces and ibe work daily time-schedule; . data on the electric endw (voltage Jevel, maximum contractual demand, annual electric energy consumption) m order to character

47、ue the user as ckclric load; data on availability in the plant of means fix miEigMiDg Power Quality probtenu such as motor-generaror (M-G) sets, unintOTUptible power supply UPSX harmonic fikas; data relevant to provision of PQ monttoring equipment or recordii% devices allowing perfbnn statislicE of

48、PQ events inacimg production process; - data on the averse annual occutrence of PQ ewnu and lheir typical duiaiton; data on economical evatuatKm of the disturtanocs impact, wilt tbu possibilily of ghloK a smmaf* qualitative cost, ehetw the exact costs were not available; - data on the crirical outag

49、e duration that prevents the regular (kvelopfnent of the activity process* Prcsewly, the PQ sunty is being lo be exiended lo several further Italian Districts, such as Vcneio (Northern Disvkl and Campaaia (Southern Wsirict) 1 lowcw, for the survey already completed in the above mentioned Districts.

50、93 bdu$tnal customers have returned questionnaires. About 50% of (he sample takes part of the Lombardia Dbtrkf. The rcmaini: 50% includes indusinal manuuXurers of Lazio and Abruzzo. The sample inctudes manutacturers involved m lhe activities listed in the Tabic 1 TatfclA npe rf iW.ftCERN whnj in 就 P

51、Q idaKiwiM fof oil jndw$tiv (vfeUr ecodKtioa dedne huuuhold appJMXci pMec rmner wood tt；sdov VBCS clecujdicabki median Lcal look (Udrio) devices food pvo4i*cti uxatd pnKinn k Mructurci Kxtilc fibres chfoucali notirboMt ctomic products iron iMUlurev pkirnucMlich IV. SURVEY RESULTS The returned questi

52、onruire forms have been processed m order to obtain some main relt iHuarated below. In Table H the br 曲 4 S : SOO-IOX) 2CX60 flO-tOO 3I0M 2CX60 I fl叭 Ur 1000, must ereipy gisii m (WAt Fig 2 Annual (ggy ccmumiKM gImsci found in (he ample 82%Mihc tanci- 35 20100 206300 OWN. onioctL3i demnd Mnt tXW) Fi

53、g. 3. Maximum antnUttd JanaikJ Umit cUmcs found (99% vC (he tnplc) Availability in Ihe pUni of motor-generator (M-G) seis and UPS has beat also invcsligated* with the results reported in Tables V and VL UPS iutabiliiv 1%) rci 39 Ko 61 TZe VI 26 TT TabteV N F0 5, UPS lypical rations (46 VPS found). 7

54、5-150 150-300 Cl E Harmonic filters haw bgn found available in in 1971 In 1972 be joined the OipeitmcnU) di Ekiumccn心 of the MHccnico di Milaiw where be is Fu(! Profcffior at preseni His area of research inclttdc power ckctfockics and power syacm harmonics, lie i$ a mnnber of (he hilhia National Re:

55、wch Cuundl (CNIt) group of Eleotkal Power Systons and of the luliaa Ekcvicdl Association (AEI). E-mail Qftdreu: tirenvbettaniwtpoiimi it Dario Zoolidli (SM%) rweivcd Ihe PhD degree in Electrical Enginccrins in 1989 from Ihe Holitccnico dl Milano Since 1998 be has been Asodatc Proftswr in Qie UtpvUnK

56、fito di EkRMeciUa of PoJilccnko di Milano. Hb areas of research include povtr syiicm harmcnics. and power system Ally血 DrZanhwlli h 8 senior member of IEEE PCS and IAS, member of the Italian National Research Council (CH.lt) jremp of Rkctrical Powtr Systems and of the Italian BIcctried Assodalion (A

57、Elh EaftgUrM、sanig”晦b0t5.c忆ept或mzt O-78O3-6672-7. or $10.00 (O 2001 IEEE943 73 外文译文 工业用户电能质量成本的调查 摘要：在当今竞争激烈的电力市场的新风景下电能质量(PQ)成本问题在工业客 户中显得尤为重要。事实上，它代表了一种附加的价值，即电力供应商可以提供给用 户，以提高他们的生产率和可靠性的工业过程。行业必须比较与来自损坏的产量和成 本的公用事业提供的电能质量服务成本。而后者是容易界定，因为它是一些设备的物 理替代关系，生产成本不是简单的赚钱。本文是在意大利的工业化地区的中等大型工 厂进行的调査，用我们的U标

58、去确立是什么类型的干扰导致了生产力的损失，并估计 其相关的损失。 一、引言 自从过去十年，人们越来越关注在电能质量在最终使用的各个环节存在问题的调 査。 笫一次相关的努力在这个领域应该重新去在美国和加拿大进行调査。随后的儿次 调査U前正在进行。随着电力行业的放松管制和舆论不断增长的全球，电力公司电能 质量可以确定新的机遇并且发挥新的竞争作用。 不断增加的工业用途大大提高生产自动化渗透生产过程电能质量扰动的敏感性。 因此，对客户来说干扰通常涉及高成本，主要是III于可能损坏设备，产品损失、恢复 时间的负担。 正确的识别成本归因于电能质量问题的确是一个至关重要的一点，因为它在客户 合作关系在这个领

59、域中代表了非常初步的步骤。 本文的主要U的是说明位于三个意大利区的一项调査在某些工业客户(意大利北 部和中部)，为了开始判断真正的低廉的电能质量成本。 这项调査是基于问卷调査的各个方面的问题。特别是，第二部分一些技术解决方 案,通常釆用工业设施最小化或消除客户接触电能质量问题，简要总结了成本的关注。 第三部分，详细说明了组织调查和问卷的主要结构。最后。第四部分,对处理活动的 结果应用于样品山返回的问卷调査进行了说明和讨论。 应该强调的是，从机构协会的生产商和个人企业本身的能够积极性很好的接受。 这种接受，征求延长调查，进一步了解北部和南部在全国各地以讣划短期电能质量的 可用性状态的方式在意大利

60、工业部分电动最终用途。 一、电能质量的改善行动 供电公司在供应的质量中起着关键的作用，即使后者还受其他企业的约束，为用 户的工作条件下，设备、丄厂和当局设定的限制标准。图1中总结了电能质量电动分 销公司的管理过程。为代表的主要任务是监测活动(a) 8,收集的数据来自网络3) 和用户(f)。在这个领域，作者提出的问卷调査可以代表一个有趣的经济贡献,干扰对 客户的操作的经济影响的佔计。 假设信息监控给定起源的状态变量(b),因为它们描述了完整体系的“网络用户”。 质量不满童时，应可能的行动计划(e)。 Fig 1. Power Quality rnanncment Migran 图一:电能质量管理