第2章数学模型基础

数学模型基础：同步发电机数学模型

说实话，搞电力系统仿真这么多年，我最大的体会就是——数学模型是根。根扎不牢，后面什么优化策略、智能控制都是空中楼阁。今天咱们就聊聊同步发电机的数学模型，这块内容我当年啃了整整三遍才真正吃透。

你想想看，三相交流电在abc坐标系下是时变的，正弦波来回振荡。直接分析？太复杂了。Park变换说白了就是「换个角度看问题」——把定子侧的abc三相绕组，等效成跟转子一起旋转的dq0坐标系。

为什么要这么干？因为变换之后，那些时变的电感参数就变成了常数。嗯，这里要注意：Park变换不是近似，是严格的数学等价。

核心公式：

[fd, fq, f0]^T = P * [fa, fb, fc]^T

其中变换矩阵P为：

P = (2/3) * 
[cosθ    cos(θ-120°)   cos(θ+120°)
 sinθ    sin(θ-120°)   sin(θ+120°)
 1/2     1/2           1/2       ]

我在项目中遇到过一件事：有同事直接用MATLAB自带的Park变换函数，结果仿真结果怎么都对不上。后来一查，原来是功率守恒型和幅值守恒型两种变换没搞清楚。我个人习惯用功率守恒型，因为后续算功率时不用再乘系数。

避坑指南：我曾经在调试一个水电站仿真模型时，发现励磁电压波形异常。查了两天，最后发现是Park变换的初始角度设错了。记住：θ = ωt + θ0，θ0取决于转子初始位置。

派克方程是在Park变换基础上建立的。它描述了发电机在dq0坐标系下的电压、磁链关系。说白了，就是告诉你「给转子加多少励磁，定子能出多少电」。

电压方程（标幺值形式）：

ud = -Rs*id - ψq*ω + dψd/dt
uq = -Rs*iq + ψd*ω + dψq/dt
uf = Rf*if + dψf/dt

磁链方程：

ψd = -Xd*id + Xad*if
ψq = -Xq*iq
ψf = -Xad*id + Xf*if

看着有点晕？我刚开始学的时候也晕。后来我师傅教我一招：把派克方程想象成「直流电机」的方程。因为经过Park变换后，交流量变成了直流量，分析起来就跟直流电机一样直观。

注意：派克方程中的电抗参数（Xd、Xq、Xad等）都是饱和值。我在做大型汽轮发电机仿真时，曾经忽略了饱和效应，结果算出的短路电流比实测值大了15%。后来加了饱和修正，精度才上去。

励磁系统的作用，说白了就是控制发电机端电压。我见过最经典的模型是IEEE推荐的ST1A型（静止励磁系统）。

典型传递函数：

Efd(s) = Ka/(1 + Ta*s) * (Vref - Vt + Vs)

其中：

我个人习惯把励磁系统分成三类：

经验之谈：我曾经调试过一个励磁系统，PSS投入后反而振荡加剧。后来发现是励磁系统的限幅环节没建模——实际中励磁电压不可能无限大，上限一般是4~5倍额定值。这个限幅不加上，仿真结果就是错的。

原动机和调速器，管的是发电机的转速和功率。水轮机、汽轮机、燃气轮机，模型各有不同。我重点说说最常见的汽轮机-调速器模型。

调速器模型（简化版）：

Pm = P0 - (1/R) * (ω - ω0)

R是调差系数，一般取4%~5%。什么意思？就是转速每变化1%，功率变化20%~25%。

汽轮机模型（串联多级）：

G(s) = 1/(1 + Tch*s) * (Fhp + Fip/(1 + Trh*s) + Flp/(1 + Tco*s))

这里Tch是蒸汽容积时间常数，Fhp、Fip、Flp分别是高、中、低压缸的功率占比。我在做一次调频仿真时，发现如果忽略再热器时间常数Trh，动态响应会偏快很多。

关键参数表（典型值）：

说了这么多，我画个图帮你理清思路。同步发电机数学模型，核心就是三大块：发电机本体（派克方程）、励磁系统、原动机+调速器。三者通过电气量和机械量耦合在一起。

这张图我画了好几次才满意。你看，三大模型通过转子运动方程耦合在一起——原动机给机械功率Pm，发电机输出电磁功率Pe，差值决定了转速变化。励磁系统则通过调节励磁电压来影响机端电压和无功功率。

最后提醒一句：模型再漂亮，参数不准也是白搭。我曾经接手过一个项目，仿真结果跟实测总是差一截。后来发现是发电机Xd参数用了出厂值，但机组运行十年后饱和特性已经变了。所以，模型参数一定要跟现场实测数据对标，这是仿真可信度的底线。

公众号:蓝海资料掘金营，微信deep3321