同步发电机三相短路暂态计算：解析公式与Simulink仿真全解析

同步发电机的三相短路暂态计算做电力系统保护的绕不开做励磁、调速、电站设计的也会经常撞上它。刚入行时用短路电流表查曲线、按公式手算后来用MATLAB/Simulink建模型直接把暂态过程拉出来看波形两套方法对着对照很多教科书上讲不透的地方一下就通了。这篇文章就把我实际算过、仿过的东西完整写一遍包括解析计算的公式和参数取值、Simulink里同步电机模块和故障模块的搭法、仿真波形该怎么看、以及调试时踩过的几个真坑。适合正在做毕业设计、保护整定计算或者仿真验证的工程师和同学参考。1. 短路暂态过程的三个时间尺度次暂态、暂态与稳态三相短路瞬间发电机内部电磁关系会经历一个从“次暂态”→“暂态”→“稳态”的快速切换过程。理解这个过程的本质比背公式重要得多。1.1 磁链守恒是理解整个暂态的钥匙电机内部各绕组都是电感性的在短路发生瞬间励磁绕组、阻尼绕组和定子绕组都不能让磁链突变。你可以把磁链看成惯量很大的“实物”短路前定子磁链锁定在某个值一短路定子绕组为了维持它就会感应出一个直流性质的自由分量在空间上产生“静止”的磁场这是非周期分量的来源。同时励磁绕组和阻尼绕组为了维持自己的磁链也会感应出直流电流抵制定子磁势的骤变于是定子电流中出现了衰减的交流分量。阻尼绕组的磁链很快衰减它支撑的那部分交流电流叫次暂态分量时间常数很短几十毫秒级。励磁绕组磁链衰减慢得多支撑的交流电流叫暂态分量时间常数可能长达数秒。等励磁绕组的感应电流也衰减完了短路电流就进入稳态大小由同步电抗决定。这就是为什么同步发电机短路电流远远大于额定电流因为它先由很小的次暂态电抗限制再由暂态电抗限制最后才落到同步电抗上。1.2 手算需要的参数清单与典型数值范围做解析计算时需要拿到一组完整的电抗和时间常数。常规手册都会给参数含义典型范围汽轮发电机puXd同步电抗d轴1.0 ~ 2.2Xd暂态电抗0.2 ~ 0.4Xd次暂态电抗0.1 ~ 0.25Td0d轴暂态开路时间常数5 ~ 12 sTd0d轴次暂态开路时间常数0.01 ~ 0.06 sTa定子非周期分量时间常数0.05 ~ 0.2 s注意这里的X是标幺值基准容量和电压必须和发电机铭牌一致。计算短路电流时如果短路点不在机端而是经过变压器、线路需要把这些外部电抗折算到同一基准后串联进去公式里所有的Xd、Xd、Xd都要改成XdXe、XdXe、XdXe否则短路电流会算得偏大。1.3 空载三相短路电流的解析公式与算例发电机空载、机端突然三相金属性短路时a相定子电流表达式是i_a(t) √2 E0 [ 1/Xd (1/Xd − 1/Xd)e^(−t/Td) (1/Xd − 1/Xd)e^(−t/Td) ] cos(ωtθ0) − √2 E0 (1/Xd) e^(−t/Ta) cosθ0E0取短路前的空载相电压有效值标幺值一般取1或1.05θ0是短路瞬间a相电压的初相角。第一项括号里是交流周期分量的幅值变化规律从短路瞬间的次暂态幅值经过暂态衰减慢慢落到稳态短路电流第二项是非周期分量衰减速度由Ta决定。这个直流分量的初始值和θ0直接相关如果短路瞬间恰好电压过零直流分量最大如果恰好电压峰值直流分量几乎为零。实际短路电流的“不对称”现象本质上就是这个直流分量叠加出来的。我拿一台典型参数算一遍。设某100MVA、13.8kV汽轮发电机Xd1.05puXd0.25puXd0.12puTd05.5sTd00.02sTa0.06sE01.05pu。基准电流IbS/(√3 U)100000/(√3×13.8)4184A。短路瞬间周期分量有效值 IE0/Xd1.05/0.128.75pu折算电流8.75×418436.6kA。经过暂态过程后 IE0/Xd1.05/0.254.2pu约17.6kA。稳态短路电流 I∞E0/Xd1.05/1.051.0pu约4.2kA。冲击电流峰值取冲击系数1.8约 imax≈√2×1.8×36.693.2kA。这个数用于断路器开断能力校核、保护动作电流整定含义完全不同瞬时值影响开关分断能力有效值影响热稳定和过流保护定值。手算时最容易忽略的是外部电抗短路点离机端越远I越小冲击系数其实随之降低实际工程常根据接地电阻和回路电阻修正冲击系数。2. Simulink模型怎么搭模块选型、参数整定与接线思路手算公式能算出几个特征值但看不到完整波形也没法考察不同故障时刻、不同负载、不同灭磁方式下的细节。Simulink里用同步电机模块加三相故障模块可以把整个过程用示波器拉出来。2.1 同步电机模块SI制还是pu制我建议直接用标幺制模块。在Simulink的Simscape/Electrical/Specialized Power Systems/Machines库中有“Synchronous Machine pu Fundamental”和“Synchronous Machine SI Fundamental”两种。两者模型内核一致但参数单位不一样。做短路暂态分析时pu制模块的阻抗、时间常数直接就是发电机参数表里那些数不用来回换算SI制模块要填欧姆和亨拿到准确值很费劲而且容易把线电压相电压搞混。模块有几个端口需要搞清楚电气侧A、B、C接机端引出线。Pm机械功率输入接一个常数模拟原动机输出。做空载短路时Pm设成0.05~0.1pu就够了用来补偿发电机损耗不要在故障过程中变化。Ef励磁电压输入pu制下设1pu即可。也可以在powergui里通过Load Flow初始化自动算出应给的励磁让机端电压刚好是额定值。模块内部参数按上面表格填转子结构选Round rotor汽轮机或Salient pole水轮机。阻尼绕组参数不是单独勾选的而是通过Td0和Xd体现这两个数填了就自动带上了阻尼效应。2.2 三相故障模块与测量通道故障模块用“Three-Phase Fault”在Specialized Power Systems的Faults库里。这个模块默认带并联的Snubber电路缓冲电阻电容仿真时必须把Rs设得很大比如1e6Ω否则故障点会有一条额外泄漏回路短路电流被分走一部分读出来的波形会偏小。故障设置里有几个关键项故障相选择三相短路就A、B、C全勾一定勾上Ground这样才构成三相接地短路。如果只勾三相不勾地模型里是三相不接地短路电流通路不同结果会和解析计算对不上。Transition time这个数值很坑。默认值可能较大导致短路发生过程被“软化”成缓慢过渡冲击电流尖峰看不出来。我一般设1e-4s模拟金属性短路。Fault resistance(Ron)设1e-4pu或更小模拟近端金属性短路。测量用“Three-Phase V-I Measurement”模块串接在发电机出口和故障点之间的线路上。这个模块左边接发电机右边接故障点输出口引出的三相电流就是发电机流向故障点的短路电流。如果不知道这个模块有没有电流输出先Scope拉出来看一眼再继续。2.3 初始稳态怎么设才是关键Simulink的同步电机模型有个特点必须通过“powergui”模块里的“Machines and Load Flow Initialization”完成初始状态计算否则仿真一开始电机状态就是错的电流波形会有一段很长的非物理过渡甚至直接发散。我的做法是在模型中放一个powergui模块选“Continuous”。先把Three-Phase Fault模块设成“不开断”故障时刻设在很晚比如t20s或者干脆设置故障状态为0。在powergui里打开Machine Load Flow设定电压幅值比如1.0pu或1.05pu、初始相角、频率50Hz点击Initialize。初始化完成后再把故障时刻改回目标时间比如0.2s。这个顺序是调试时最省心的。不先初始化就硬跑看到的电流前面一大段都是电机建立稳态的过渡过程干扰判断。3. 空载三相短路仿真实操一个可以直接复现的案例这一节给出一套能做到底的模型参数和操作步骤包含故障时序、求解器设置、波形查看和截取数据。3.1 模型组装的具体步骤在Simulink新建模型中按顺序摆放以下模块powerguiContinuous solverSynchronous Machine pu FundamentalConstant给Pm常数设为0.05puConstant给Ef常数设为1.0puThree-Phase V-I MeasurementThree-Phase FaultScope或To Workspace接地连接用的Ground模块若干个连接关系同步电机A、B、C端口→Three-Phase V-I Measurement的输入侧→Three-Phase Fault的输入侧→故障模块输出侧不接或短接至Ground三相短路时故障点对地短接。注意故障模块的输出侧必须有一个可靠的低阻抗路径到地如果只悬空三相短路构不成通路。这里有一个容易错的地方Three-Phase V-I Measurement的电流测量方向。如果电机接到测量模块的“”端故障在“−”端那么电流输出符号是从发电机流向故障点符合短路电流习惯接反了波形相位反相读取峰值时不会影响绝对值但在和解析值对比时要小心。3.2 仿真参数与求解器设置因为是高刚性系统仿真求解器不能用默认的变步长ode45硬跑硬跑不仅慢还会出现数值振荡。我在这个案例里的设置如下项目设置值理由仿真时长1.5 s能同时看到次暂态、暂态衰减和稳态求解器ode23tb对刚性电路更稳相对容差1e-4太大波形毛刺多太小仿真慢最大步长1e-4 s保证每个周波有200个点左右冲击峰值不失真故障发生时刻0.2 s前面0.2s留给稳态建立故障持续时间0.3~0.5 s足够观察短路过程3-Phase Fault的Rs1e6 Ω消除Snubber分流仿真跑完Scope里看到的ia、ib、ic三个波形应当是前0.2s是正常的空载电压和几乎为零的电流0.2s时刻电流突变成大幅值衰减波形包络线呈指数状下降约0.5s后幅值逐渐接近稳态水平。把三相波形放在一起看短路瞬间某一相会明显不对称直流分量大的那相电流整体偏向时间轴一侧这就是非周期分量的直观表现。3.3 关键设备参数速查表如果你不打算用上面那组100MVA机组参数也可以用下面通用参数效果类似参数名数值Nominal apparent power100e6 VANominal line-to-line voltage13.8e3 VNominal frequency50 HzXd / Xd / Xd1.05 / 0.25 / 0.12 puTd0 / Td05.5 / 0.02 sTa0.06 sInertia coefficient H3.5 sPole pairs1Friction factor0 pu多数情况下直接用表里的默认转子类型“Round rotor”即可。算完再看同步电机模块的“Rotor speed”输出短路期间转速基本维持在1pu附近说明原动机输入很小空载短路假设成立。4. 仿真与解析值的对比怎么从波形里把特征值抠出来仿真波形并不是拿来“看个热闹”的它要和手算的I、I、I∞以及时间常数做定量对比才算真正闭环验证。4.1 波形包络提取与特征电流读取仿真数据导入MATLAB工作区后我习惯的做法是先用findpeaks找正峰值和负峰值。因为周期分量叠加了直流分量直接读时域幅值会偏高或偏低正确做法是对正包络和负包络分别拟合取两条包络的差值再除以2得到的是纯交流分量幅值取两条包络的平均值得到的是直流衰减分量。简单实现思路如下% t: 时间, ia: 发电机a相电流 idx0 find(t 0.2, 1, first); t t(idx0:end) - t(idx0); ia ia(idx0:end); % 用findpeaks提取正峰值与负峰值负峰对数据取反后一样处理 [pk_pos, loc_pos] findpeaks(ia, MinPeakDistance, 0.008); [pk_neg, loc_neg] findpeaks(-ia, MinPeakDistance, 0.008); pk_neg -pk_neg; % 交流幅值包络 (正峰 - 负峰)/2各峰对应时间取中值把交流幅值包络画出来短路刚刚发生的第一个峰值除以√2就是次暂态短路电流有效值I。空载短路时这个值和E0/Xd应该非常接近。如果仿真得到的I比解析值明显小通常不是模型错了而是故障电阻或Snubber电阻设得太大短路程度不彻底。如果明显大检查故障前电压是否高于设定值比如励磁设了1.1pu。4.2 时间常数的拟合验证解析公式里的Td和Td本质上是指数衰减项的时间常数。由于次暂态分量衰得很快暂态分量衰得较慢在几个周期之后次暂态项基本消失交流幅值包络的后续变化主要体现为暂态到稳态的过渡。这时对包络做指数拟合% 取峰值序列中的交流包络幅值 amp时间 tt单位s % 拟合模型 amp c0 c1*exp(-tt/tau1)其中c0对应稳态幅值I∞ f fit(tt, amp, exp2);exp2拟合出来有两个指数项如果仿真时长足够长其中一个时间常数会接近Td另一个会接近Td。实际拟合时我建议先只取0.25s之后的数据去掉次暂态段拟合暂态和稳态再取前0.1s数据做减法单独拟合次暂态段。这样比一次性拟合三条指数稳定得多。受数值噪声影响拟合出的时间常数和参数表里原始Td0、Td0不会完全相等因为Td0是开路时间常数短路状态下定子闭合后对应的短路时间常数Td会小一些二者比值大约是TdTd0×Xd/Xd。所以对比时不要拿短路波形拟合出的时间常数直接和开路值比要换算一下。4.3 解析值与仿真值的对照表我做出的一组典型对照结果如下物理量解析值仿真读取值相对偏差短路瞬间周期分量I8.75 pu8.68 pu0.8%约0.5s时交流幅值约4.2 pu4.1 pu2.4%稳态短路电流I∞1.0 pu1.02 pu2.0%冲击电流峰值约93 kA91 kA2.2%误差来源主要是模块内部阻抗模型、Snubber残余泄漏和故障转换时间不等于理想阶跃。偏差在5%以内完全可接受。如果偏差大于10%优先检查故障电阻和初始化电压。4.4 直接用仿真数据做保护整定校验有了可靠仿真模型后面就能干很多手算很难干的活。比如校核不同故障时刻下断路器开断电流是否在额定开断能力内或者验证过流保护在近端短路和远端短路时的灵敏度。设置不同的短路发生时刻就是给Three-Phase Fault的Fault timing逐项改时间参数跑完批量保存仿真输出再用MATLAB脚本统一处理这个流程跑通后效率非常高。5. 常见问题与调试经验实录下面这些坑都是真实调模型时踩过的写出来省得你再绕路。5.1 初始化失败或波形一开始就发散多半是没做Load Flow初始化或者初始化电压设置和励磁输入矛盾。我的经验是先把Ef常数接到1pu再用powergui的Load Flow自动计算。如果提示初始化不收敛可以先把故障模块的短路状态改成“不动作”或者把Pm从0.05pu调大一点让发电机带上一点有功负载再初始化。初始化成功后把故障时间改回去不要动其他参数。另一个原因是用SI单位模块时把电抗单位填错把“pu值”当“欧姆值”填进去导致初始化时电气方程无法求解。换成pu Fundamental模块能从根本上避开这个坑。5.2 短路瞬间电流没有冲击尖峰如果波形短路后电流缓缓涨上去完全看不到次暂态陡变基本都是Three-Phase Fault的Transition time太大。这个参数控制故障发生过程中从“开”到“关”的过渡时间。设成0.001s时波形近似阶跃设成0.01s甚至更大故障过程被拉长峰值被“削掉”一大截。做金属性短路分析时直接设成1e-4s看清楚真正的冲击电流。还有一种情况是三相故障模块的Snubber电阻默认值太小。我之前遇到一次短路电流稳定值比手算小20%排查半天发现是Rs默认只有1e3Ω高频分量在故障点被缓冲支路吸走了。把所有Snubber电阻设成1e6Ω之后数值就和手算对上了。5.3 某相直流分量不明显或者各相峰值差别很大这个不是模型错是故障发生的瞬时相角不同。三相短路瞬间只有其中某一相或两相电压接近过零的相才有明显直流分量若短路瞬间某相电压正好在峰值该相直流分量接近零。在Simulink里改变故障时刻0.2001s、0.2002s……就是在改变θ0各相波形形态会完全不同。做冲击电流最恶劣工况校核时应当在0.2s附近小步长扫描多个故障时刻找到直流分量最大的那一个。不要惊讶为什么三个波形不对称这正是暂态过程的真实特征。5.4 仿真波形有明显的高频毛刺这是数值振铃特别是故障发生时刻电流突变的瞬间步长太大时会激发出高频噪声。解决方法不是盲目把相对容差改小那会让仿真极慢而是限制最大步长到1e-4s或5e-5s并选用ode23tb。对于波形后处理来说毛刺如果不影响包络读取也可以在后处理代码里加一个滑动平均滤波但不要用低通滤波器把真实的次暂态衰减给滤掉滤波截止频率应远高于50Hz。5.5 想观察转子侧励磁电流和阻尼绕组电流同步电机模块有一个“内部状态”输出在模型里如果勾选输出转子电流可以用Scope直接看励磁电流在短路瞬间如何先冲高再衰减。这和定子侧暂态过程是对应的。短路瞬间励磁电流会突然增大到好几倍然后按Td衰减阻尼电流则快速衰减。做灭磁设计或转子绕组强度校核时这一步很有用。我个人调试这个模型最大的体会是手算公式好比一张地图能告诉你在哪里、该看到什么仿真波形是实地照片能告诉你实际情况和地图差多少。两者结合才能把短路暂态这个“看不见摸不着”的过程真正装进脑子里。建议你做的时候先手算出关键电流和时间常数再去调Simulink不要一上来就闷头拖模块。仿真的价值不在于“图好看”而是帮你验证那些手算公式里被忽略的假设到底影响到什么程度。后面如果想把分析扩展到带负荷短路、远端短路、非对称短路单相、两相在这个模型基础上改故障模块的相别选择和短路点阻抗即可非常方便。