4. 实时仿真基本原理

RTDS 的核心任务，是在真实时钟推进下求解电力系统的电磁暂态过程。它与离线仿真最大的区别不是“算得更准”，而是“必须在下一个仿真步到来之前算完”。理解这一点，才能理解步长、实时性和数值稳定性之间的相互制约关系。

一、电磁暂态仿真的本质

电磁暂态仿真关注的是系统在微秒到毫秒尺度上的电压、电流变化，典型对象包括线路分布参数、变压器饱和、开关动作、换流阀换相、故障暂态等。其数学本质是求解一组微分—代数方程：

dx/dt = f(x, u, t)      （状态方程，描述电感电流、电容电压等）
0     = g(x, u, t)      （代数方程，描述节点电压约束）

RTDS 采用 Dommel 提出的伴随模型法，把电感、电容等储能元件在每一个步长内等效为电阻与历史电流源的并联支路，从而把微分方程离散成节点导纳方程 Y·V = I。以电感为例，其梯形积分伴随模型为：

i(t) = (Δt / 2L)·[v(t) - v(t-Δt)] + i(t-Δt)
等效为：G = Δt/(2L)，I_hist = i(t-Δt) + G·v(t-Δt)

电容同理，G = 2C/Δt。这意味着：步长 Δt 直接进入元件等效导纳。步长变了，整个网络的 Y 矩阵和历史电流项都要重算。这是 RTDS 中步长不能随意切换的根本原因。

二、步长选择：精度、稳定与实时的三角约束

步长选择不是“越小越好”，而是要在三个约束之间取平衡：

• 精度约束：步长应远小于被仿真动态的最小时间常数。工程经验是 Δt ≤ τ_min/10。例如 50 Hz 系统基波周期 20 ms，若只关心工频，Δt = 50 μs 已足够；但若含 2 kHz 电力电子开关，开关周期 500 μs，则 Δt 通常取 1～5 μs。
• 稳定性约束：显式积分法存在稳定域限制，梯形法虽为 A 稳定，但对刚性系统仍可能产生数值振荡。线路行波传播时间 τ = l/v 是硬约束，Δt 必须小于最短线路的传播时间，否则行波模型失真。
• 实时性约束：RTDS 每个步长内必须完成采样、求解、输出全过程。若单步计算耗时 T_calc，则必须满足 T_calc < Δt，否则实时性被破坏，仿真进入“过载”状态。

三、实时性约束的工程含义

实时性不是“算得快”，而是“算得准时”。RTDS 采用固定步长、等间隔推进，每个步长的时间预算被严格切分：

单步时间预算 Δt = T_采样 + T_求解 + T_输出 + T_余量
若 T_求解 + T_输出 > Δt - T_余量 → 过载（Overrun）

过载的典型表现是：仿真时钟落后于真实时钟，波形出现周期性“卡顿”或相位漂移，与外部物理设备接口时产生数据错位。排查时应先看 RTDS 的过载计数器，再定位是网络规模过大、开关器件过多，还是接口通信占用了过多时间。

四、数值稳定性：从“能算”到“算得可信”

数值稳定性指误差在迭代中是否被放大。RTDS 中常见的失稳来源有三类：

1. 刚性系统：系统中同时存在快、慢时间常数（如大电网含小电感支路），显式法会因快模态而失稳。RTDS 采用隐式梯形法，具备 A 稳定性，可抑制此类失稳。
2. 开关与非线性元件：开关动作造成网络拓扑突变，若步长过大，会在开关时刻产生虚假尖峰。解决方法是减小步长或在开关附近采用插值处理。
3. 数值振荡：梯形法对电感电流在突变后可能产生 2Δt 周期的振荡。工程上通过串联小阻尼电阻、或改用后向欧拉起步来抑制。

五、步长、实时性与稳定性的联动关系

三者构成一个闭环：步长减小 → 精度提高、稳定性改善，但单步计算量增大 → 实时性压力上升；步长增大 → 实时性宽松，但精度下降、可能失稳。RTDS 工程实践中的典型决策路径是：先按被仿真对象的最小时间常数确定精度上限，再按网络规模和硬件资源确定实时性下限，最后在两者之间取满足稳定性的折中值。

六、工程实践中的参数示例

以一个含双馈风机的 110 kV 并网系统为例，说明步长与稳定性的实际取舍。系统含 2 kHz 变流器开关，最小时间常数为变流器 LCL 滤波器的谐振周期约 200 μs。按 Δt ≤ τ_min/10 原则，步长应不大于 20 μs；但变流器开关周期 500 μs，为准确捕捉开关纹波，实际取 Δt = 5 μs。此时单步计算量约为 20 μs 步长下的 4 倍，需评估 RTDS 机架的计算余量。若过载，可将变流器部分放入小步长接口，主网仍用 20 μs，通过接口变压器解耦。

判断仿真是否可信，不能只看波形“像不像”，而应检查三点：过载计数是否为零、开关时刻电流是否出现非物理尖峰、能量是否随时间单调衰减或守恒。若三者均满足，说明步长、实时性与稳定性达到了工程可接受的平衡。