3. 同步发电机暂态建模与参数获取
暂态稳定分析的核心是“能量—转矩—功角”的动态平衡,而同步发电机正是这一平衡的物理载体。工程上并不需要把发电机内部所有电磁细节都搬进仿真,而是要在“足够准确”和“可计算”之间取一个平衡点。目前电力系统机电暂态仿真普遍采用d-q轴派克模型:把定子三相量通过派克变换投影到随转子旋转的d轴(直轴)和q轴(交轴)上,从而把时变电感问题转化为常系数问题。理解这一模型,是理解后续励磁、调速、稳定判据的基础。
一、d-q轴模型:从物理结构到状态方程
派克变换的参照系约定必须首先明确:d轴取转子主磁极中心线(直轴),q轴沿转子旋转方向超前d轴90°电角度,转子角速度ω对应同步转速。定子电压方程在d-q坐标系下写为:
u_d = -R_a·i_d - ω·ψ_q + dψ_d/dt
u_q = -R_a·i_q + ω·ψ_d + dψ_q/dt
其中ψ_d、ψ_q为d、q轴磁链,R_a为定子电阻。磁链由各绕组自感与互感耦合决定,工程建模的关键在于保留哪些绕组。常用的三阶模型保留励磁绕组f和d轴阻尼绕组D、q轴阻尼绕组Q,忽略定子暂态(即令dψ_d/dt、dψ_q/dt为0),得到“三阶实用模型”;若再忽略阻尼绕组,则为“二阶模型”,仅保留励磁绕组动态。阶数越高,对次同步振荡、励磁快速响应等细节刻画越准,但参数需求也越多。
需要特别注意的是,模型中所有电感、磁链均以标幺值表示,且d、q轴采用各自基准。工程参数通常以“不饱和同步电抗X_d、X_q”“暂态电抗X'_d、X'_q”“次暂态电抗X''_d、X''_q”以及对应时间常数T'_d0、T''_d0、T'_q0、T''_q0给出。这些量并非独立,而是由绕组漏抗与互感组合而成,辨识时需保持物理一致性。
二、励磁系统与调速器模型
发电机模型给出“本体”,励磁和调速则给出“控制”。励磁系统决定机端电压E_q'的动态,工程上常用IEEE AC1A/AC4A/ST1A等标准模型,核心环节包括:测量滤波、电压调节器(含放大倍数K_A和时间常数T_A)、励磁机(K_E、T_E)、励磁限幅与软反馈(K_F、T_F)。以AC1A为例,其传递函数可简写为:
V_err = V_ref - V_t + V_stab
V_R = K_A · V_err / (1 + T_A·s)
E_fd = f(V_R, K_E, T_E, S_E)
其中V_stab为励磁附加稳定信号,S_E为励磁机饱和函数。参数获取时,K_A、T_A通常由励磁系统厂家试验报告给出;K_E、T_E需通过励磁机空载特性曲线拟合;限幅值则直接取自调节器整定单。调速器模型常用IEEE GGOV1或IEEEG1,核心是“转速偏差→调门开度→机械功率”的传递链,关键参数包括调差系数R、伺服时间常数T_servo、再热时间常数T_rh、高压缸功率占比F_hp等。这些参数直接影响暂态过程中机械功率的恢复速度,进而影响加速面积。
三、从铭牌、试验与辨识报告获取工程参数
工程参数获取遵循“先铭牌、再试验、后辨识”的递进逻辑。铭牌提供额定容量S_N、额定电压U_N、额定功率因数cosφ_N、额定转速n_N、直轴同步电抗X_d(有时给出饱和值与非饱和值)、惯性时间常数H(部分铭牌给出飞轮转矩GD²,需换算)。换算关系为:
H = (GD² · n_N²) / (4 · S_N) [单位:秒,GD²单位t·m²,S_N单位MVA]
若铭牌只给X_d的饱和值,需根据空载特性曲线反推非饱和值。短路试验是获取暂态参数的主要手段:三相突然短路电流的包络线可分解为次暂态、暂态和稳态三个分量,分别对应X''_d、X'_d、X_d,时间常数T''_d、T'_d由包络线衰减斜率提取。频率响应试验则用于辨识阻尼绕组参数和q轴参数,通过在不同频率下注入小信号,拟合得到X_q、X''_q、T''_q0等。参数辨识报告通常由电科院或厂家提供,包含辨识方法、拟合误差和推荐值,工程使用时应核对辨识工况(饱和程度、温度)是否与仿真工况一致。
四、参数一致性检查与典型示例
拿到参数后不能直接填入仿真,必须做一致性检查。常用判据包括:X''_d < X'_d < X_d,X''_q < X'_q < X_q,T''_d0 < T'_d0,且X_d/X_q比值与机组类型相符(凸极机约1.0~1.2,隐极机约1.0~1.1)。若出现X'_d > X_d或时间常数倒置,说明参数来源或单位换算有误。以下为一台300MW汽轮发电机的示例参数(合理假设,非实测):
| 参数 | 数值 | 单位 | 来源 |
|---|---|---|---|
| X_d | 1.85 | pu | 铭牌(非饱和) |
| X'_d | 0.28 | pu | 三相短路试验 |
| X''_d | 0.19 | pu | 三相短路试验 |
| T'_d0 | 6.5 | s | 短路试验包络线 |
| T''_d0 | 0.045 | s | 短路试验包络线 |
| H | 4.2 | s | 铭牌GD²换算 |
| K_A | 200 | pu | 励磁厂家报告 |
| T_A | 0.02 | s | 励磁厂家报告 |
将上述参数代入三阶模型,在0.1s机端三相短路、0.2s切除的工况下,可观察到功角δ先快速增大后回摆,励磁电压E_fd在故障期间被顶到限幅值,故障切除后按T'_d0衰减恢复。若把T'_d0误填为0.65s(少一个数量级),仿真中功角首摆会明显偏小,导致稳定裕度判断偏乐观。因此,参数获取的最后一步不是“填表”,而是用典型扰动工况做一次合理性校核:检查功角首摆幅值、励磁顶值倍数、机械功率恢复时间是否与同类型机组经验范围一致。只有通过校核的参数,才具备进入暂态稳定计算的资格。