4. 同步发电机实用模型

同步发电机的完整数学模型（Park 模型）包含定子电压方程、转子励磁绕组方程、直轴阻尼绕组方程、交轴阻尼绕组方程以及转子运动方程，共七阶。在电力系统动态仿真中，直接使用完整模型会带来较大的计算负担，而且很多动态过程对仿真结果影响很小。因此，工程上常根据研究目的对模型进行降阶简化，形成六阶、五阶、四阶、三阶和经典二阶等实用模型。理解每种模型的简化条件，是正确选择模型并解释仿真结果的前提。

一、从完整模型到实用模型的简化逻辑

Park 模型在 dq0 坐标系下描述同步发电机，其核心状态量包括：转子角速度 ω、转子功角 δ、励磁绕组磁链 ψ_f、直轴阻尼绕组磁链 ψ_D、交轴阻尼绕组磁链 ψ_Q，以及定子磁链 ψ_d、ψ_q。定子绕组的暂态过程对应电磁暂态，其时间常数远小于转子机械和励磁暂态，因此在机电暂态仿真中通常忽略定子磁链的暂态项，即令 dψ_d/dt = dψ_q/dt = 0。这一步简化把定子方程从微分方程变为代数方程，是几乎所有实用模型的共同基础。

在此基础上，是否保留阻尼绕组、是否保留励磁绕组暂态、是否考虑次暂态过程，决定了模型的阶数。六阶模型保留全部转子绕组动态；五阶模型忽略 q 轴阻尼绕组；四阶模型进一步忽略 d 轴阻尼绕组；三阶模型忽略全部阻尼绕组，只保留励磁绕组动态；经典二阶模型则进一步忽略励磁绕组暂态，把励磁磁链视为恒定，只保留转子运动方程。

二、各模型的简化条件与适用场景

六阶模型保留励磁绕组、d 轴阻尼绕组和 q 轴阻尼绕组的暂态，能够描述定子电流的次暂态变化，适用于研究短路电流初始阶段、次同步振荡以及需要精确描述阻尼绕组作用的场合。其状态量为 ω、δ、ψ_f、ψ_D、ψ_Q，加上定子代数方程，共六阶微分方程。

五阶模型忽略 q 轴阻尼绕组，即令 ψ_Q = 0。它适用于 q 轴阻尼效应不明显、或研究重点在 d 轴动态的场合，例如某些汽轮发电机的励磁调节研究。五阶模型仍保留 d 轴阻尼绕组，因此 d 轴次暂态过程仍可体现。

四阶模型同时忽略 d 轴和 q 轴阻尼绕组，只保留励磁绕组暂态和转子运动方程。它适用于阻尼绕组影响较小、研究时间尺度较长的机电暂态问题，如电力系统低频振荡分析。四阶模型的状态量为 ω、δ、ψ_f。

三阶模型在四阶模型基础上进一步忽略励磁绕组暂态，即认为 ψ_f 恒定或由励磁系统快速调节。实际上三阶模型通常指保留励磁绕组动态但忽略阻尼绕组，而经典二阶模型才忽略励磁暂态。这里需要区分：三阶模型保留 ω、δ、ψ_f 三个状态量，适用于励磁系统动态不可忽略的场合。

经典二阶模型只保留转子运动方程，即 ω 和 δ 两个状态量，认为暂态电抗后的电动势 E′ 恒定。它适用于远离故障点、励磁调节作用较慢或研究重点为功角稳定性的初步分析。经典二阶模型计算量最小，但无法反映励磁和阻尼绕组的快速动态。

三、MATLAB 代码实现示例

下面以经典二阶模型为例，给出 MATLAB 代码实现。假设单机无穷大系统，发电机经变压器和线路接入无穷大母线，参数为：H = 5 s，D = 0.5，X′_d = 0.3，E′ = 1.2，无穷大母线电压 U = 1.0，线路电抗 X_L = 0.4。状态量为 δ 和 ω，输入为机械功率 P_m。

function dydt = gen2nd(t, y, Pm, H, D, Eq, U, X)
% 经典二阶模型：y(1)=delta, y(2)=omega
delta = y(1);
omega = y(2);
Pe = Eq * U / X * sin(delta);   % 电磁功率
ddelta = omega - 1;              % 转子角速度偏差
domega = (Pm - Pe - D*(omega-1)) / (2*H);
dydt = [ddelta; domega];
end

% 仿真调用
H = 5; D = 0.5; Eq = 1.2; U = 1.0;
X = 0.3 + 0.4;   % 总电抗
Pm = 0.8;
y0 = [30*pi/180; 1.0];   % 初始功角30度，初始转速1.0
[t, y] = ode45(@(t,y) gen2nd(t,y,Pm,H,D,Eq,U,X), [0 10], y0);
plot(t, y(:,1)*180/pi);
xlabel('时间 t / s'); ylabel('功角 delta / deg');
title('经典二阶模型功角响应');

运行上述代码后，可以观察到功角在扰动后经历衰减振荡并趋于新的稳态值。若将模型替换为三阶或四阶，只需在状态方程中增加励磁磁链或阻尼磁链的微分方程，并相应修改电磁功率和励磁电压的计算。例如三阶模型需要增加励磁绕组方程：dψ_f/dt = ω₀(E_fd - ψ_f/T′_d0)，其中 E_fd 为励磁电动势，T′_d0 为 d 轴开路暂态时间常数。

四、模型选择与参数准备要点

在实际仿真中，模型阶数的选择还需要考虑参数的可获得性。六阶模型需要 d 轴次暂态电抗 X″_d、q 轴次暂态电抗 X″_q、d 轴开路次暂态时间常数 T″_d0、q 轴开路次暂态时间常数 T″_q0 等参数；而经典二阶模型只需要 H、D、X′_d 和 E′。如果参数手册中缺少次暂态参数，强行使用六阶模型反而会引入不准确的结果。

此外，不同模型对初始条件的处理也不同。经典二阶模型需要给定初始功角 δ₀ 和初始转速 ω₀；三阶及以上模型还需要给定初始励磁磁链 ψ_f0，通常由稳态运行条件反推。在 MATLAB 中，可以通过潮流计算得到发电机端电压和电流，再计算各磁链初值，确保仿真从稳态出发，避免初始瞬间的数值振荡。

总之，六阶、五阶、四阶、三阶和经典二阶模型构成了同步发电机从精细到简化的模型谱系。理解每一步简化的物理含义和适用边界，才能在仿真中做出合理选择，并用 MATLAB 正确实现。