2. 电力系统仿真基础

在PSCAD中建立任何模型之前，必须先回答三个基础问题：要研究的是哪一类暂态过程？求解器如何把连续时间的微分方程变成可计算的离散序列？仿真步长和初始化参数应当怎样设置才能既保证精度又不至于让仿真慢到无法接受？这三个问题分别对应电磁暂态与机电暂态的概念区分、PSCAD求解器原理，以及步长与初始化设置。理解它们，是后续所有建模工作的前提。

一、电磁暂态与机电暂态：两类时间尺度的过程

电力系统中的“暂态”并不是单一现象。按所关注物理量的变化速度和所保留的元件模型精度，通常划分为电磁暂态和机电暂态两类。

电磁暂态关注的是电压、电流瞬时值随时间的快速变化，时间尺度通常在微秒到毫秒级。典型过程包括：雷击或操作过电压引起的行波、断路器开断时的电弧与重击穿、换流阀换相过程、变压器励磁涌流、电缆与架空线上的波过程。这类过程必须保留线路的分布参数特性（用贝杰龙模型或频率相关模型表示）、电感和电容的储能效应，以及三相之间的电磁耦合。PSCAD正是为电磁暂态仿真而设计的工具，其核心求解对象是完整的微分方程组。

机电暂态关注的是发电机转子角、转速、有功/无功功率等“机械—电气”耦合量的变化，时间尺度通常在几十毫秒到几十秒。典型过程包括：发电机功角摇摆、低频振荡、电压稳定、频率稳定。这类过程通常忽略线路的波过程，用准稳态的相量模型表示网络，重点保留发电机转子运动方程和励磁、调速系统的动态。机电暂态一般由BPA、PSASP或PSS/E等机电暂态程序处理，PSCAD也可以做，但需要简化网络模型并采用较大的步长。

两者的关键区别可以归纳为：电磁暂态保留瞬时值和分布参数，步长小、频带宽；机电暂态使用相量和集中参数，步长大、频带窄。选择哪一类，取决于研究目标——若关心过电压、谐波、换相失败，必须用电磁暂态；若关心功角稳定和系统频率，机电暂态更合适。

二、PSCAD求解器原理：从微分方程到离散序列

PSCAD的求解器本质上是一个“数值积分器”。电力系统元件（电阻、电感、电容、变压器、电机、线路）在PSCAD中都被描述为微分方程或代数方程，整个系统构成一组形如 dx/dt = f(x, t) 的状态方程。求解器的工作，就是把这组连续方程在时间轴上离散化，逐步求出每个时间点的状态量。

PSCAD默认使用隐式梯形积分法。以电感为例，其伏安关系为 v = L·di/dt，用梯形法离散后可得：

i(t) = i(t-Δt) + (Δt / 2L) · [v(t) + v(t-Δt)]

整理后，电感在每一步都可以等效为一个电阻与一个历史电流源的并联：

i(t) = v(t) / R_eq + I_hist
R_eq = 2L / Δt
I_hist = i(t-Δt) + v(t-Δt) / R_eq

电容和电阻同理，都可以转化为诺顿等效支路。这样，每一步求解就变成了解一个线性方程组 [Y]·[V] = [I]，其中 [Y] 是节点导纳矩阵，[I] 是包含历史项和电源的注入电流向量。PSCAD在每一步先更新历史项，再求解节点电压，然后回代求出各支路电流，最后推进到下一步。

隐式梯形法的优点是数值稳定、无衰减，适合电力系统这种刚性方程；缺点是每一步都要解方程组，计算量比显式方法大。PSCAD还提供后退欧拉法作为可选方案，它更稳定但会引入数值阻尼，适合处理强非线性或开关频繁动作的场景。

三、仿真步长选择：精度、稳定与效率的平衡

步长 Δt 是PSCAD仿真中最重要的参数之一。它决定了仿真的时间分辨率，也直接影响计算量和数值稳定性。步长选择需要同时满足三个条件：

• 捕捉最快动态：步长应小于所关心最高频率分量周期的1/10~1/20。例如要观察2 kHz的谐波，周期为0.5 ms，步长应取25~50 μs。
• 满足数值稳定：对于隐式梯形法，步长过大可能导致开关事件被“跳过”，产生虚假的数值振荡。一般要求步长小于最小时间常数的1/5。
• 控制计算量：步长减半，计算量约翻倍。在保证精度的前提下，应尽量取较大步长。

工程中常用的经验值如下：

需要特别注意的是，PSCAD中还有绘图步长（Plot Step）和通信步长（Communication Step）的概念。绘图步长只影响输出曲线的采样密度，不影响求解精度；通信步长用于与外部程序交互。求解步长（Time Step）才是真正决定数值积分精度的参数，三者不要混淆。

四、初始化设置：让仿真从稳态出发

电磁暂态仿真通常从稳态开始，否则会因初始条件不匹配而产生虚假的暂态过程。PSCAD提供了两种初始化方式：

1. 自动初始化（Auto Initialization）：PSCAD先运行一个“初始化阶段”，把系统在工频稳态下求解一遍，得到各节点电压、支路电流和状态变量的初值，再以此为起点开始暂态仿真。这是最常用的方式，适用于大多数线性或弱非线性系统。

2. 手动初始化：对于含有强非线性元件（如避雷器、饱和变压器、电力电子开关）的系统，自动初始化可能不收敛。此时可以手动指定关键元件的初值，例如发电机转子角、励磁电压、变压器分接头位置，或者先运行一个较长的“预热”阶段让系统自然进入稳态。

在PSCAD的Project Settings中，与初始化相关的关键设置包括：

• Start-up Time：仿真开始时间，通常设为0。
• Initialization Time：初始化阶段时长，默认0.1~0.5 s，复杂系统可延长。
• Snapshot File：保存稳态快照，后续仿真可直接从快照启动，节省初始化时间。
• Time Step：求解步长，需与初始化阶段一致。

五、一个可复现的参数设置示例

假设要仿真一条10 kV配电线路在断路器合闸时的励磁涌流，线路长度5 km，变压器容量2 MVA，空载合闸。合理的设置如下：

Time Step (Δt)        = 5 μs
Solution Time         = 0.5 s
Initialization Time   = 0.2 s
Integration Method    = Implicit Trapezoidal
Plot Step             = 50 μs
Snapshot              = 启用，保存于0.2 s
线路模型              = Bergeron（贝杰龙）
变压器模型            = 带饱和特性的UMEC模型
断路器合闸时刻        = 0.2 s

这样设置的理由是：5 μs步长可以分辨涌流中约2 kHz的高频分量；0.2 s初始化让变压器铁芯进入稳态；0.2 s合闸后观察0.3 s的涌流衰减过程；绘图步长50 μs使曲线文件不至于过大。仿真结束后，应检查合闸瞬间的电流峰值是否在合理范围（通常为额定电流的6~8倍），并观察涌流是否在几个周波内衰减。若峰值异常或衰减过慢，应检查变压器饱和参数和线路模型。

电磁暂态与机电暂态的区分决定了模型精度和步长量级，PSCAD求解器通过隐式积分把微分方程转化为逐步求解的代数方程组，步长选择需要在精度、稳定和效率之间取得平衡，而初始化设置则保证仿真从正确的稳态出发。这四点共同构成了PSCAD仿真的基础，后续的建模、调试和结果分析都建立在此之上。