《电力系统仿真》讲义

1、KMDt X ITCFF<第一章系统仿真的基本概念1- 1系统仿真的定义系统仿真定义建立系统模型（数学模型、物理效应模型或数学-物理效应混合模型），并在模型上进行试验。系统仿真范例水利方面在葛洲坝和三峡大坝建设前，许多单位都建立起缩小了的长江水道模型和大坝模型，并进行冲水试验，获得流体力学和结构力学等相关参数，为工程设计提供依据。电力方面在动模实验室中用电动机-发电机组对实际发电厂进行模拟，并且用电感-电容n形结构模拟输电线路；在高压实验室中用冲击发生器模拟雷电；在计算机上用EMTP软件对电力系统电磁暂态过程进行数值分析，用 PSASP软件对电力系统机电动态过程进行数值分析。国防方面用巨

2、型计算机按数学特性相同的原则模拟核裂变（聚变）反应过程，这样可以减少真正核试验的次数，甚至不进行核爆炸试验也能发展出先进的核武器。系统仿真依据最基本的依据是相似定理。从应用角度分析几种相似等效方法：几何比例相似。军事指挥员的沙盘演习。特性比例相似。两个系统运动的物理本质完全不同，但具有相似的微分方程，且参数一一对应，我们称这两个系统的特性比例相似。电气系统机械系统微分方程d2Xdt2DdXKX F(t)Ldq Rdqdt2dt1Cq E(t)参数对应距离Xh电荷q速度dX/dt/电流dq/dt外力F（t）/电源E(t)质量M/电感L阻尼系数D 电阻R弹簧系数K1/电容1/C图1-1两个系统特性

3、比例相似注：动模试验也是根据特性比例相似的原则，这个原则可理解为真实系统与模拟系统具有相同的无量 纲（标幺值）方程。感觉相似。主要是视觉、听觉、触觉和运动感觉相似，是人在模拟环境中的仿真，特别是用 各类模拟器件对操作人员进行训练的依据。如航天员在宇航中心培训，宇航中心就是一个虚 拟太空环境。逻辑思维方法相似。对获取的信息进行分析、归纳、综合、判断直至操作控制的方法相似。 例如机器人。微分方程的数值解法、离散相似法。如数值方法中的龙格-库塔法、差分法、非线性曲线的分段线性化等。这是数字仿真的基础。电力系统仿真就是以电力系统为仿真对象的系统仿真。1- 2系统仿真的作用人类认识或研究、开发一个系统可

4、以通过理论推演或实物试验的方法进行，但对于一些大的、复杂系 统，如一个电网，无法得到其数学模型的解析解，有的子系统甚至无法得到可信的数学模型，而由于多种 条件限制，实物试验不可能做或做起来困难很大，这样就只能借助于模型试验（即仿真）来达到认识或研 发一个系统的目的。通常在下列情况时一般考虑用模型试验（仿真）而不用实物试验：1）系统还处于设计阶段，并没有真正建立起来，因此不可能在真实系统上进行实验；2）在真实系统上做试验会破坏系统运行。例如在正常运行的电网中做一个短路故障试验，有可能造 成系统剧烈振荡乃至崩溃；3）如果人是系统的一部分时，由于知道自己是试验的一部分，行为往往会和平常不同，因此会影

5、响 试验效果，这时最好将人也建立模型。例如一些基于专家系统的软件或装置；4）在实际系统上做多次试验时，很难保证每次的操作条件都相同，因而无法对试验结果做出正确判 断。例如一个电器产品为什么要通过型式试验和出厂试验，就是为了人为地设置一些确切的操作 条件，以获得对未来使用性能的正确判断；5）试验时间太长或费用太高或有危险。例如“神州 3号、4号”上的模拟人；6）系统无法复原。例如，想要知道变电站运行中的变压器绝缘到底能够耐受多高冲击电压，不可能 真的对正在运行的变压器施加很高的冲击电压。有些学者以较为学术化的术语归纳了仿真方法的适用情况：1）不存在完整的数学公式，或者还没有一套解答数学模型公式的

6、方法；2）虽然可以有解析法方法，但数学过程太复杂，仿真可以提供比较简单的求解方法；3）解析解存在而且是可能的，但超出了个人的数学能力，因而应该估计一下，建立仿真模型、检查 并且运行模型的费用比起向外求助以获得解析解，何者合算；4）希望在一段较短的时间内能观测到过程的全部历史，以及估计某些参数对系统行为的影响；5）难于在实际的环境中进行实验观测；6）需要对系统或过程进行长期运行的比较，而仿真则可以随意控制时间，使它加快或减慢。1-3系统仿真的分类仿真可以有多种分类方法。按系统模型的类型，可分为：1）连续系统仿真一一系统模型以微分方程描述；2）间断（事件）系统仿真系统模型以面向事件、面向进程、面向

7、活动的方式描述；3）连续/间断（事件）混合系统仿真；4）定性系统仿真一一系统模型以模糊理论等描述。按仿真的实现方法和手段，可分为：1）物理仿真；2）计算机仿真，又称数学仿真；3）实物在回路中的仿真，一般称为半实物仿真；4）人在回路中的仿真。物理仿真要求模型与原型有相同的物理属性，其优点是模型能最真实全面地体现原系统特性，缺点是 模型制作复杂、成本高、周期长、灵活性差；计算机仿真的优缺点正好与物理仿真相反。1-4计算机仿真的三个要素与三项基本活动计算机仿真（即数学仿真）采用数学模型，它是用数学语言描述系统行为的特性。计算机仿真的三个要素是：系统、模型、计算机。联系它们的三项基本活动是：模型建立、

8、仿真模型建立（又称二次建模）、仿真试验。三个要素和三项基本活动的相互关系如图1-2。图1-2计算机仿真的三要素和三项基本活动从下一章起，主要内容是针对连续系统的计算机仿真模型的建立。第二章集中参数网络的瞬态分析概念理解一：连续系统仿真模型的建立，就是将系统的微分方程模型按某种数值计算方法处理成时间 离散的代数表达式，易于计算机逐个步长地递推求解。概念理解二：微分方程初值问题数值求解，一般指的是，由变量在初始时刻to的值(已知)，按确定时间步长At(通常设为固定步长)并根据系统计算机模型依次求出变量在ti、t2、.、tk、.的值。2- 1电阻元件R电阻元件的模型(支路电压-电流关系)为u(t)=

9、Ri(t)(2-1)其时间离散的表达式为uk=Ri k(2-2)仿真模型等值电路的建立如图2-1所示I*u(t)离散近似.FoikUkKt)R图2-1电阻的计算机仿真模型建立例2-1图2-2中，Us(t)为u-t函数已知的电压源，求系统的计算机仿真模型。G4图2-2由电阻元件与已知电源组成的系统Us(t)解由电阻的仿真模型等值电路列系统节点方程U1.kG1 G 2G 2G2G2 G3 G4 u2.k这是一个很容易由计算机编程求数值解的离散形式的线性代数方程组，此模型即为所求。 推论由电阻和已知电源构成的线性电路，可直接作为求数值解的计算机模型电路。2- 2电感元件LL的电压-电流关系为u(t)

10、或 i(t)Ldi(t)dt1 ti(t t) t tu(t)dt(2 3)为方便起见，用 表示法与电流值一样，2-3)写成:ik表示t (或tk)时刻的电流值、ik-1表示t- At或tk-1)时刻的电流值，电压值的下标 将1 ttk udtL tk 1(2 4)tkt udttk 12 (Uk 1Uk)将（2-4）写成：t /k 1U k )i ki k 1-(U2Lt(ikt)Uk1U k 1 )2L2L上式又可写成：ik gLUkIsl（2 5）其中gLt2L和 Isl ik 1gLUk1公式（2-5）计算机数值递推模型也称为电感元件如图2-4所示：用梯形近似积分法（图 2-3）有:图

11、2-3梯形积分L的“瞬态伴随模型”，由数学模型转化为数值模型1blUk!+ 吵:o iAS1«1/ikgL U(t) MS)图2-4电感的仿真模型建立Isl图中右边的数值模型等效电路（电感 L的“瞬态伴随模型”等效电路），由一个等效电阻（电导）gL和一个等效电流源Isl并联组成，在进行第 k步Uk、ik计算时，Isl已由前一步的Uk-i、ik-1得到，是已知的。2-3 电容元件Cc的电压-电流关系为i(t)C du(t) dt或 u(t)1u(t t) Cti(t)dt(26)同样，用梯形积分近似法，可导出C的数值递推公式（电容元件 C的“瞬态伴随模型”）：ik gcUk1 SC(2

12、 7)其中gc 2Ct其等效电路为和IscgcUk i+ i(t)u(t)Uk2-4互感元件电压-电流关系为U山伯山 其中u Ui, UnT,i ii, inT dt不难导出：t 1 tgjL1r (2 8)2 2数值模型的递推公式与电感L类似，只是矩阵形式而已。2-5非线性电阻电力系统中主要的非线性电阻元件是避雷器，其电压-电流关系为u CiKu1/解非线性代数方程，一般用牛顿迭代法 设非线性方程为F ( x) =0又假设x(0)为方程解的初始“猜测值”，则对真值的第一次逼近 x(1)= x(0)-Ax(0)则 F (x(0) - AX) =0按泰勒级数展开并略去高阶项，得(0)F(x)F(

13、x(0)(1)(0)F(x(0)F(x(0)经过(k-1)次迭代后(k)(k 1)xF(x(k 1)F(2 9)实际计算中预先规定精度刍当达到(k)(k 1)(k 1)停止迭代，x(k)就是非线性方程的解。例 2-2图2-6 含非线性电阻的回路非线性回路方程F(u) U0 R1i u U0 R1Ku1/u 0套牛顿迭代公式（2-9)得u(k) u(k1)F(u(k 1)dFduu(k 1)(2 10)u（k-1）点处的考察di/du具有电导量纲，将其第（k-1）次迭代时的值以Gd（k-1）表示，它相当于非线性电阻在 动态电导。因此dFdU u(k 1)RiGdk1)1 代入(2-10)得(k)

14、(ku1)Uoi(k1)R1 u(k1)1 RQdk 1(2 11)进一步有(丄 Gdk 1)u(k)& (i(k1)Gdk1)u(k 1)R1R1这是一个标准的节点方程形式，不难看出非线性电阻通过牛顿法被“线性化”了（如图2-7 ）。ls=i(k-1)-Gd(k-1)u(k-1)图2-7 非线性电阻的“线性化”计算模型这里特别要注意，牛顿迭代法中的u（k）和前面梯形积分法中的uk区别：在由uk-1求uk的一个步长中，由于非线性元件的存在使得差分方程为非线性，因此必须用牛顿迭代法 来解方程。一个可行的处理方法是：设 u(0)=uk-1.牛顿迭代.令 uk=u(k)IIII1 11 1u

15、k-12-6非线性电感元件特性f(i)(2 12)磁链是电压的积分(t)(ttt) t tu(t)dt由梯形积分法得(t)-tu(t)2(t t)其中(tt)t2u(t t) (t t)将上式代入(2-12 )得(tt)f(i(t)经过整理，可写成非线性电阻和一个已知电源的形式2 2u(t)屮)(t t) (2 13)例2-3已知非线性电感在to时刻的初值io、uo以及特性f(i)，问如何算出此后时刻的值?0 f(io)2uo在tl时刻的方程表达式为Ui2tf(ii)由牛顿迭代解非线性电阻方法，与电路的其它方程式联立得到比、i12-7三相常规n形电路图2-8为实际计算中常用的三相耦合n形电路，

16、经常作为结构对称的交流输电线路的集中参数模型,线路耦合参数R、L、C分别称为串联电阻、串联电感、并联电容。RLRRsRmR MR MRsR MR MRmRs自电阻Rs1(2RposRzero)互电阻 Rm 3(Rzero Rpos)LLsL ML ML MLsL ML ML MLs自电感Ls-(2L3posiL zero )互电感 L M - (L zero L pos )CCsCMCMCMCsCMCMCMCs对角元Cs3(2CposCzero）非对角元Cm1 (Czero C pos )实际系统的参数均以正序、 以C矩阵中的Cm为负数，其数值为相间电容。至于n形电路的数值模型（瞬态伴随模型）

17、，其公式可由单相 n形电路导出，只是 R、L、C参数换成零序给出，串联阻抗参数有Zpos<Z zero, 而并联电容参数有Cpos> Czero, 所矩阵，而i、u变量换成列向量。第三章 均匀单导线中暂态过程的计算在这一章里，主要对单导线-地系统的波过程建立数值模型。3- 1长线过渡过程的解一、长线的微分方程式Ro、Lo、Co、Go,该线路以长度设有一条单导线线路，它的单位长度电阻、电感、电容、电导分别为 微元表示的等效电路如图3-1所示。Rodx Lodx-amrvGodx Codxii dx u 1 xu dxxx%图dx3-1单导线线路取电流的正方向为 x增加的方向,KVL、

18、KCL可写出如下方程:Roix-GoUxL。Co(3 1)这就是长线方程。严格来说，导线存在集肤效应和电晕效应，大地也非理想导体，因此R。、L0、C。、G0应该是电压、电流波形的函数。但在这里，为掌握一般规律起见，Ro、Lo、Co、Go均视为常数，这也能够符合绝大多数研究项目的精度要求。二、零初始条件下长线方程的解(象函数形式) 设S为拉普拉斯算子，dU(x,S)dxdI(x,S)由拉氏变换将(3-1 )变成常微分方程(考虑零初始条件)乙(S)l(x,S)(3 2)再令dx其中：Zo(S)Y°(S)U(x,S)Ro SLo Yo(S) GoSCo.Zo(S)Y 0(S)Y(S)Z(S

19、).乙(S)/ Yo(S)则(3-2)式的象函数形式的通解为F2(S)eY(S)xU(x,S) R(S)e Y(S)x I(x,S) F1(S)e Y(S)xZ(S)式中F1(S)、F2(S)为待定象函数，由边界条件确定。Y ( S)、Z ( S)为拉氏运算形式的传播常数和特征Y(S)xF2(S)e(3 3)阻抗，其含义为:传播常数丫（S）特性阻抗Z(S),Zo(S)Yo(S)( Ro SLo)(Go SCo)1 1 (S )v. Zo(S)/Yo(S).( Ro SLo) (GoSCo)其中 VZc为传播速度RoLoRoLoLo£G oCo为衰减系数为畸变系数为无损线的特性阻抗由于

20、Ro=o， Go=o故Y(S)Z(S)S. L oCoS/vZCL o / C o代入（3-3）式得xS_U(x,S)E(S)evF2(S)ev1xxSS-I(x,S)vF2(S)ev利用拉氏变换基本定理f(t)变换eSF(S)写出时域中的通解形式三、无损长线微分方程的通解形式(34)U(X,t) fi(tX)f2(tvx1xxi(x,t)fi(t ) f2(t )ZcvV(35)其中，f1、f2由边界条件确定。从这个解的形式，我们可归纳如下几点物理概念：x ）表示前行电压波，沿 x的正方向传播;V1) f1(tf2（t ）表示反行电压波，沿 x的负方向传播。2)前行波、反行波在无损导线中传播

21、时不发生畸变和衰减。3)前行电压波U1伴随一个前行电流波i1,它们之间关系U1i1Zc(3 6)反行电压波U2伴随一个前行电流波i2,它们之间关系U2Zc(3 7)注意：负号是由于反行电流波方向与电流参考方向相反i24)无论对于前行波或反行波，电压波与电流波的比值在数量上等于|Zc|，但对于线路上某点的电压和电流由前、反行波合成时却不存在这个关系，即U(x,t)i(x,t)Zc(3 8)3- 2无损单长线的贝杰龙模型这里针对无损单长线，介绍如何将分布参数的输电线路化为集中参数计算模型，因为该方法出自于贝 杰龙数学模型，所以又称为贝杰龙法。ikm(t)imk(t)Uk(t)Um(t)图3-2长度

22、、波速、波阻抗分别为I、v、Z的单长线如图所示的无损长线，两端节点分别为(3-5)式得u(x,t)Zi(x,t) 2fi(tu(x,t)Zi(x,t) 2f2(tV)V)(3 9)上式表明，当tx禾口t x各为固定值时,Vxxf 1 (t)和f 2 (t )也各为固定值，进而VVu(x,t)+Zi(x,t) u(x,t)-Zi(x,t)也各为固定值。其含义是：若观察者沿x以速度V移动，则在他所在位置观测到的u(x,t) ±(x,t)对该观察者而言分别始终不变。设Tl/V，假如观察者在t- T时刻从k点出发，在t时刻到达m点，根据上述推理并结合考虑(3-9) 式，有Uk(t ) Zkm

23、(t ) Um(t) Z imk(t)将上式改写为imk(t)1 -u Z同理可得其中Im(tm(t) Im(t|uk(tikm(t)(3 10)其中ikm(t)ZUk(t)Ik(t(3 11)Ik(timk (t)1) 計住I m图3-3贝杰龙模型的等效电路由贝杰龙模型的数学表达式（3-10 ）、（ 3-11 ）及其等效电路图3-3可总结其特点:1）k和m是一个线路元件的两端节点，在计算时可按两个分离节点计算；2）无损单长线的贝杰龙模型由两个分开的诺顿等效电路构成，其中诺顿等效电阻值等于线路波阻抗，诺顿等效电流源由另一端前T时刻的电压、电流值确定，从数值角度看，对当前计算步长而言，电流源是已

24、知的。第四章多导线系统中的暂态过程计算本章主要针对三相架空输电线路。4- 1三相输电线路微分方程的矩阵形式图4-1为对称三相输电线路一个dx长度元的等值电路。Lt4dx图4-1三相线路的一个长度元注意：图中的相间电容K在以下均表示为 C。省略与单相长线方程一样的推导过程，我们直接得出方程的矩阵形式UxixRiGu（4 1）其中，变量的列向量u UA（X,t） UB（x,t）Uc（x,t）Ti iA（x,t） iB（x,t） ic（x,t）T参数矩阵R AAR ABR ACRR BAR BBR BCR CAR CBR CCL AAL ABL ACLL baL BBL bcL CAL CBL cc

25、C A0CabC ACRCABC ACG A0G abG ACGG ABG ACCabC acC B0C abC bcC bcC bcC c 0 C acC bcG ABG ACG B0 G ABG BCG bcGbcG C 0G ACG BC零初始条件下（4-1）的象函数矩阵形式为dU(x,S)dxdl(x,S)dx进一步对x求导得Z(S)l(x,S)Y( S)U(x,S)(42)d2U 今 ZYU PU dx2 嘎 YZIPTIdx2一般情况下，P护T，但如果Z、Y方阵为平衡阵，则(43)P=PT。注：所谓“平衡”方阵就是所有对角元素相等、所有非对角元素也相等的方阵。例如三相全换位线路 的

26、Z、Y就可视为平衡阵。4- 2多导线系统的相-模变换一、问题的提出（4-2）、（4-3）与单导线微分方程形式类似，但方程中的变量为向量、系数是矩阵。由于系数矩阵中 有非对角元素，解某一相要受到其它相的牵连，所以直接解方程是十分困难的，如果我们借助于相似变换 将矩阵中的非对角元素化为零，则矩阵方程组中的每个方程只有一个电压、电流，方程形式与单导线微分 方程完全相同，解答方法与前述的一样。设Um、Im表示变换后的电压、电流向量（称为模量） ，即U=TuUmI=TiIm其中Tu、Ti分别为电压、电流的变换矩阵。将（4-3）式用模量表示d Um Tu1PTuUmmdx2d2I(4 4)mdx2对于三相

27、全换位系统，P为平衡阵，P=PT,所以Tu=Ti=T，变量的变换矩阵是相同的。Ti1pTTiIm、变换矩阵T根据矩阵理论中的相似变换基本定理, 非正交变换矩阵略去烦琐的推导过程，直接得变换后各模量功率之和不等于变换前各相量功率之和正交变换矩阵1T J3 1T121_212032:32变换后各模量功率之和等于变换前各相量功率之和无论是非正交变换矩阵 概念的理解，T用于实际的数值计算。Q或正交变换矩阵T,均可用于相模变换的计算。但实际上，Q 一般用于物理三、各模的物理意义三相系统的相-模变换又称 用变换矩阵Q来写电流的相-模变换,ABC? a30) °1 m1Q 1II m321 260

28、IaIbIc进一步写成单个方程形式11 m13(IaIb Ic)1 m21AIb Ic)1 m32( ibIc)2第一个模分量是以大地为回路的“地中模量”，这一分量的波的传播与大地有关，是a（0系统中的0模，可简称为“地模”；第二、三个模分量是以线间为回路的“空间模量”，这两分量的波的传播与大地无关，是a 0系统中的a B模，可简称为“线模”；模分量的构成如图 4-2所示。0ABBCC0模分量Ba模分量四、相-模变换与对称分量参数 由于4-2(X、B、0模分量1TZT 1- 131201ZsZMZ M11132Z MZsZ M2121.3Z MZMZS03_22_2Zs 2Zm000ZS Z

29、M000ZS Z M那么对（4-2）中的第一式（第二式方法相同，后不重复叙述）进行相 -模变换，有dUdxTZT 1Im写成三个独立的模量方程dUm0dxdUdxdUdx(Zs2ZM ) 1 m00(ZsZM )l m0(ZsZ M )l m0mm(4 5)既然Zs、Zm分别为三相对称线路的自阻抗和互阻抗，且有零序阻抗Zzero=Zs+2ZM=Zo正序阻抗Zpos=Zs-ZM=Z 1所以，地模分量方程是由零序参数确定的微分方程，线模分量方程是由正序参数确定的微分方程。这与我们所熟悉的对称分量法十分相似，只是对称分量法用于单频复域的相量值稳态计算，而相-模变换用于时域的瞬时值暂态计算。五、三相输

30、电线路方程的通解Z。丫。00由于TPT 10Zi Yi000ZiYiPm 所以由（4-4 ）得dU mo(x S)戶ZoYoUm0(X,S) 0dx2dU2m (x,S)占乙YiUm (X,S) 0 dxdU2m (x,S)V 乙YiUm (x,S)0dx上述模量方程的通解为(4 6)Um0(x,S)F eZoYox厂m0ieF e ZoYox厂 m02eUm (x,S)ZiYi XFm ieF e Zi丫1x m 2CUm (x,S)F e ZiYix mFm 2e ZiYix将上式写成矩阵形式U mFmi exFm2 ex这样，三相线路在 ABC坐标中的电压通解为(4 7)U TUmTF

31、miTF m2ex(48)同理，可求得电流的通解为I TI mT Z 1 Fmi人 T Z iFlm2ex(4 9)其中，模波阻抗为Z0 00Z 丫000Z 0 Zi00Zi Yi00 0Z00乙Yi当线路为无损线路时乙 L Jc0由单位长度零序电感、电容确定Z L1 C1由单位长度正序电感、电容确定六、n阶正交变换矩阵的一般形式相-模变换可推广到任意 N相多导线系统，其变换矩阵的通用式为11126J(J 1)11126.J(J 1)016T0J 1J(J 1)0N(N 1)1N(N 1)_1_NN(N 1)当N=3时，131 31 3T121一20161616这似乎与前面得到的变换矩阵不同，

32、但如果将ABC相序换作C B A、将OaB模量顺序换作0( - B) ( -a),则变换矩阵就一样了。这再次说明相-模变换矩阵不是唯一的，其原因在于：在将实对称矩阵P转化为与其相似的对角矩阵入时，特征方程detP-入=0有重根。另外，当N=2时，1_ 122T /2适用于长距离直流线路的暂态计算。112 24-3三相系统计算中的贝杰龙法在这一节里，针对无损三相线路，推导出其贝杰龙数值模型。设无损三相线路的两端分别为s、r,由单相贝杰龙模型表达式(3-10)、(3-11)直接得到三相系统中的模量方程lSr(t)1Z 000Us(t)Is(t0)lsr(t)01 Zi0Us(t)Is(t1)lsr

33、(t)0 01 ZUs(t)Is(t1)iS(t 0)1乙000Ur(t0)其中Is(t1)01 Z0Ur(t1)Is(t1)001 ZUr(t1)l0s(t)1Z0000Ur(t)I0(t0)lrs(t)01 z 0Ur(t)Ir(t1)lrs(t)001 ZUr(t)Ir(t1)i0(t 0)1 Z0000Us(t0)其中Ir(t1)01 Z0Us(t1)Ir(t1)、r» 、 . r、 r - *001 ZUs(t1)将上两式用间洁方式表示l sr.m (t)ZmUs.m(t)I s.m (tm)l rs.m (t)ZXm(t)I r.m (tm)变换到ABC系统得lsr(t)

34、TZmT 1Us(t)TI s.m(tm)Y *Us(t) TIs.m (tm)(48)lrs(t)Y *Ur(t) TIr.m (tm)1 21111Z0乙乙乙乙ZTZmT1 131 11211其中，Y*乙乙乙乙乙Z1 11112乙乙乙乙乙Z0 tlrs(tJs(tlrs(tlsr(tlsr(t-0 z .lsr(tl)l)l)l)显然，三相全换位无损线的Y*是一个平衡矩阵，其对角元素相等，非对角元素也相等，即1122y330Co：G122113y3123y32CoC101BC在等效电路中，各导纳参数Ya Yb Yc y*1 y*2y；、C 等于0模波阻抗的倒数V LoY12C1Co3 L1

35、Lo而等效电流源IAsI：(to)IBTls(t1)ICsIs(t1)IArio(to)IBTIr(t1)ICrIr(t1)y AB y BC y CA第五章 电力系统物理模拟的理论和方法电力系统的物理模拟是用相应的电气设备作为模型来对实际系统原型进行仿真。在物理 模拟中，模型与原型的区别只是大小比例的不同，而所进行过程的物理本质是完全相同的。 如实际系统中的电感必须用电感实物元件来模拟、电容用电容来模拟等。物理模拟技术是建 立在物理现象相似理论基础上的。§ 5.1 物理模拟的建模依据5.1.1 相似第一定理定理 ：相似系统有同样的相似判据，即相似指标等于 1。什么是相似判据和相似指

36、标呢？由于相似第一定理是由牛顿创立的，我们不妨以两个相 似的最简单的动力学系统为例进行研究，其中一个系统为原型，另一个系统为仿真模型。根据牛顿第二定律，两个系统的物理本质为第1个系统:d2ldt(5 1)第2个系统:F2M2d2l2dt；(5 2)如果要是这两个系统相似，则两个系统的对应变量和参数在整个动态过程中应分别保持 一个固定的比例值（又称模拟比），即FF2mFm2mM1112mit2mt代入（5-1）式得F2mFM 2mMd2l2dt：mi2mtmMmi2mt,d2l2将上式与（5-2）式比较,mMmi2 mt或呼！ 1mjmFmF(53)式（5-3）等号左边的表达式称为动力学系统的相

37、似指标。很明显相似指标表示系统各物理量 模拟比的关系。将各物理量模拟比表达式代入（5-3）式，得MM2I2F1t2F>t2这个无量纲等式就是这个简单动力学系统的相似判据。等号两边由原型系统或模拟系统的物理量构成的关系表达式为相似判据表达式，一般用n表示，式（5-1 ）或（5-2）表征的简单动力学系统的判据表达式为MlFt2idem(5 4)式中，idem为“相同”的意思，表示两个或多个系统的物理过程（现象）相同。相似判据表达式的特点是无量纲，但又说明了系统中各参与的物理量间应保持的某种固 定关系。很显然，对于复杂的动力学系统，数学描述方程为多项式，相应的相似定理也可推 论如下：有多个相似

38、判据n 1，n 2，冗k，n n。且相似系统的所有n都分别满足n =idem存在相似变换的运算法则kk j idemkidemk j1idemkKk idem相似第二定理n定理定理：假设任意物理系统是由n个量纲不同的物理量所组成，物理过程的关系由如下方 程式决定F(x 1,X2,Xn)=0若上式n个物理量中有K (K< n)个是互相独立的，如选出相互独立的 K个物理量作为基本 量，则另外n-K个物理量与选定的基本量所组成的 n-K个无量纲的比例数n i,n 2，n n-K 可以用算式完全地表达出来，而这些无量纲的比例数n就是相似判据。n定理指出如何利用量纲分析法找到描述一个物理现象的相似

39、判据的个数，并确定这些判据的表达式。对于力学系统而言，所有的力学量都是由长度I、质量M、时间t这三个基本量组成，如果用表示物理量的量纲，则任何物理量 Q的量纲都可写成Q=l “ M 卩 tY 式中a、B、丫为量纲指数。举两个简单例子来帮助理解n定理是很有必要的 例1:物理过程的关系为牛顿第二定律dt2参与的物理量为F、I、M、t四个，其中选取I、M、t三个相互独立的物理量为基本量, 则F=MIt -2进而 理，即相似判据个数为1。Ft2假如选取长度I、速度V、密度p为基本量，则有M= p I3t=Iv -1进而F= p l2v2,I2v2F可见，基本量选取的个数是唯一的，但物理性质不唯一，只要

40、相互独立即可。这一点能 够增强复杂系统物理模拟时的建模灵活性。例2:物理过程的关系为运动方程IVot 1at22参与的物理量为I、v、a、t四个，其中独立数为2，选取I、t为基本量，不难得出有vtat2两个相似判据即使不知道所研究过程的数学方程式,由上分析可推断，应用基于n定理的量纲分析法, 我们也可导出相似判据。相似第三定理定理：如果两个现象的单值条件相似，并且从单值条件引出的相似判据数值相等，那么 这两个现象相似。所谓单值条件是指一个现象从一群现象中区别出来时所需要的条件。模拟电力系统时主 要包括下列几个因素：(1) 几何相似。如果物理过程是在一个有限的、具有确定形状和大小的空间之内进行的

41、, 则在几何相似的系统中任何相应点的坐标应满足模拟比mI1112对于模拟具有集中参数的系统，则可以不要求几何相似。(2) 物理参数相似。系统原型的物理参数(主要为电阻 R、电感L、电容C)与模型中相应的物理参数应该分别满足R1R2mRLiL2mLCiC2me(3) 状态变量及其起始条件相似。物理过程一方面取决于过程的性质，另一方面也取决于起始条件。两个相似系统的相应起始条件之间的比例应该等于该物理量的模拟比。如电压、电流应满足土也mV20 v汕L mj12012（4）边界条件相似。所研究的任何一个现象都有其活动范围，在边界上可能有别的现象 存在，这些现象本身并不依赖于被研究的现象，但却对其有影

42、响，这样就不能只依靠被研究 现象的特性来判断各变量间的相互作用，而必须同时考虑边界上所进行的所有现象，即相应 时刻、相应边界点上的干扰因数应保持一定的比例系数。（5）时间相似。在随时间变化的过程里（包括暂态过渡过程），每一时刻都对应着各物理 量的一系列确定的数值。为了使物理量在模型中 t2时刻的值与在原型中tl时刻的值相似，需 保持时间模拟比固定不变。即mtt1 t2当mt=1时，模型对原型的模拟为实时模拟，否则为时谐模拟。即(5 5)mt 1实时性mt 1时谐性§ 5.2基尔霍夫系统相似判据与物理模拟比的确定方法模拟必须以相似判据为依据，求得相似判据就等于得出相似指标，也就决定了物

43、理模拟 比的选择条件。确定相似判据的方法有两种：一为基于相似第二定理（n定理）的量纲分析法，它适用 于包括特性方程无法得知的所有系统；二是基于相似第三定理的分析方程式法（主要是标幺值相等法），它适用于特性方程已知系统。以下综合两种方法推导相似判据: 任意满足基尔霍夫定律的系统，其特性均可表示为(5 6)F(R,L,C,v,i,t) 0式中R，L，C，v，i，t分别表示电阻、电感、电容、电压、电流和时间。从更广泛的动力学系统角度考虑，由于长度 I、质量M、时间t为三个基本量且R、L、C 的物质属性与I、M有关，所以上述基尔霍夫表达式可写成(5 7)F(R,L,C,v,i,t,I,M) 0不妨暂且

44、称式（5-7）为广义基尔霍夫系统。首先的问题是，上式有几个独立物理量？为 对此求解，我们可利用系统或其中各部分所遵循的物理定律和物理参数定义表达式结合量纲 分析法来推断。在式（5-7）的8个物理量中，除了 M、丨、t外，暂时假设电气量i也为基本量，则由库仑定律F q1q2 r 0F 12 2214并考虑到力F的量纲式F=Mlt -2和i=dq/dt得出量纲平衡式C=M -1l-2t4i2(5-8)进一步由电容定义C=q/v得112 丄-3. -11v=M l t i (5-9)同样，由安培力定律F12O 0l2 l 1i1dl1 (i2dl2 r0)2r得 L=M 1l2t-2i-2(5-10

45、)又根据线电荷分布导体的焦耳定律F?dl2dldtlt l式中，p为电阻率，s为线电流密度。即得出电阻R量纲式R=M 1l2t-3i-2(5-11)由(5-8)、(5-9)、(5-10)、(5-11)可列出C、v、L、R的量纲指数矩阵A1242123112221232用行列式行变换的方法容易得出行列式的秩=3由此得出结论：所有满足基尔霍夫定律的电力系统，其独立的物理量数为3个。这一点与前述牛顿力学系统相同，证实了电力系统属于动力学系统一个子集的基本概念。接下来继续推导相似判据的个数和形式。系统的电气特性是由式( 5-6)描述的，如果在 该式的6个物理量中选取3个作为基本量(这意味着广义系统中的 M、丨不作为基本量)，根 据n定理，相似判据的个数等于 3。按照电路基础理论，不论多么复杂的基尔霍夫系统都是由 R、L、C构成的，并且式(5-6) 可由R、L、C元件的欧姆定律表达式复合而成。因此基尔霍夫系统的基本表达式为diL - dt C巴 dt这样，系统的3个相似判据为RivL