2. 同步发电机数学模型：转子运动方程、励磁系统模型、原动机及调速器模型
同步发电机，说白了就是电力系统的“心脏”。它的数学模型搞不清楚，后面的稳定性分析全是空中楼阁。我这些年做电网仿真，踩过最多的坑，就是模型参数没给对，结果仿真结果跟现场实测差了十万八千里。
今天咱们就掰开揉碎，把同步发电机的三个核心模型讲透。你想想看，一台发电机并网运行，它要解决三个问题：转起来、励上磁、带上负荷。对应的就是转子运动、励磁系统、原动机及调速器这三块。
2.1 转子运动方程
转子运动方程，描述的是发电机转子“转不转得动”的问题。说白了就是牛顿第二定律在旋转机械上的应用。
我习惯用标幺值形式来写这个方程，因为实际工程中参数范围差异太大，用标幺值一眼就能看出问题。
dδ/dt = ω - ω₀
dω/dt = (Pm - Pe - D·Δω) / (2H)这里：
- δ — 转子角，单位是弧度（电角度）
- ω — 转子角速度，单位是rad/s
- ω₀ — 同步角速度，50Hz系统就是314.159 rad/s
- Pm — 机械功率，来自原动机
- Pe — 电磁功率，发电机输出到电网
- D — 阻尼系数，别小看它，有时候能救命
- H — 惯性时间常数，单位是秒
关键点：H值越大，转子越“沉”，系统越稳定。大型汽轮发电机H一般在4-8秒，水轮发电机在2-4秒。我曾经遇到一个项目，业主把H值填错了，结果暂态稳定算出来怎么都不对，查了三天才发现是参数单位搞混了。
为什么会用两个一阶方程而不是一个二阶方程？我个人习惯用这种形式，因为便于数值积分。你想想看，在时域仿真里，每一步都要更新δ和ω，拆开写程序逻辑更清晰。
工程小技巧：实际做仿真时，阻尼系数D往往很难准确获取。我建议先取D=0做保守估计，如果系统能稳定，加上阻尼只会更好。但如果是做精确校核，一定要查厂家提供的实测数据。
2.2 励磁系统模型
励磁系统，说白了就是给发电机转子“充磁”的。它控制机端电压，也直接影响系统的无功功率和暂态稳定性。
励磁系统的模型种类很多，IEEE标准里列了一大堆。但实际工程中，最常用的就三种：
| 模型类型 | 典型代表 | 适用场景 |
|---|---|---|
| 直流励磁机 | IEEE DC1A | 老机组改造项目 |
| 交流励磁机 | IEEE AC1A, AC4A | 大中型汽轮发电机 |
| 静态励磁 | IEEE ST1A | 新建机组，响应速度快 |
我重点讲一下最常用的静态励磁系统ST1A模型。它的核心就是一个一阶惯性环节加一个限幅器：
dEfd/dt = (Ka·(Vref - Vt + Vs) - Efd) / Ta其中：
- Efd — 励磁电压，直接决定无功输出
- Ka — 励磁增益，典型值200-400
- Ta — 励磁时间常数，典型值0.01-0.05秒
- Vref — 电压参考值
- Vt — 机端电压实测值
- Vs — 附加稳定信号（PSS输出）
注意：励磁系统最怕的就是“过调”。我曾经在西南某水电站调试，励磁增益Ka设得太高，结果系统出现低频振荡，频率大概0.8Hz。后来把Ka从400降到250，问题就解决了。所以参数整定一定要留有余量。
实际工程中，励磁系统还包含限幅环节。励磁电压上限一般是额定值的4-5倍，下限是0或者负值（取决于能不能强励）。这些限幅在暂态稳定计算中绝对不能忽略，否则仿真结果会过于乐观。
2.3 原动机及调速器模型
原动机提供机械功率，调速器控制转速。这两者通常放在一起建模，因为它们是“一个萝卜一个坑”的关系。
对于汽轮机，我习惯用如下简化模型：
dPm/dt = (Pg - Pm) / Tch其中：
- Pg — 调速器输出的功率指令
- Pm — 实际机械功率
- Tch — 蒸汽容积时间常数，典型值0.2-0.5秒
调速器模型稍微复杂一点，它包含测量、比较、放大、执行几个环节。我常用的简化模型是：
dPg/dt = (Kp·(ωref - ω) + Ki·∫(ωref - ω)dt - Pg) / Ts这里：
- Kp — 比例系数，决定响应速度
- Ki — 积分系数，消除稳态误差
- Ts — 伺服电机时间常数
实际经验：水轮机的调速器模型和汽轮机差别很大。水轮机有“水锤效应”，就是导叶动作后，水流惯性会导致功率先反向变化再正向变化。这个效应在长引水隧洞的电站特别明显。我做过一个项目，水锤时间常数Tw达到2秒，如果不考虑这个效应，调速器参数怎么调都稳不住。
调速器还有一个重要参数——调差系数（也叫速度变动率）。它决定了机组在孤网运行时的频率-功率特性。典型值在4%-5%之间。调差系数越小，机组对频率变化越敏感，并联运行时承担的负荷变化就越大。
2.4 三个模型的耦合关系
这三个模型不是孤立的，它们通过功率和电压耦合在一起。我画了一张图，帮你理清关系：
从这张图你能看到：
- 转子运动方程输出转速ω和功角δ，这两个量同时影响励磁系统和调速器
- 励磁系统通过控制励磁电压影响电磁功率Pe，Pe又反过来影响转子运动
- 调速器根据转速偏差调整机械功率Pm，Pm是转子运动的驱动力
嗯，这里要注意：在做小扰动稳定分析时，我们通常把这三个模型线性化，然后组合成状态空间方程。但在做暂态稳定分析时，必须保留非线性，因为大扰动下这些方程都是强非线性的。
我的建议：刚开始学的时候，先分别把每个模型的参数灵敏度跑一遍。比如固定其他模型，只改变H值，看看系统阻尼比怎么变。这样你就能直观感受到每个参数的作用。我当年就是这么干的，效果比死记硬背公式好得多。
最后说一句，这三个模型的参数获取，最靠谱的来源是厂家提供的型式试验报告。如果拿不到，才用典型值估算。但估算的结果一定要留有余地，最好做敏感性分析，看看参数偏差对结果的影响有多大。