4. 接口变量与交互机制

电磁暂态（EMT）与机电暂态（TS）协同仿真的核心难点，并不在于两个仿真器各自能否稳定运行，而在于两者在边界处“说什么、什么时候说、怎么说、说得多准”。接口变量、交互时序、数据交换格式与接口误差，共同决定了协同仿真的精度、稳定性和计算效率。本节围绕这四个要素展开，重点说明工程实现中如何选择接口变量、如何安排交互时序，以及误差从何而来、如何抑制。

一、接口变量的选择

接口变量是两侧仿真器在边界处互相传递的物理量。常见候选包括电压、电流、功率和角度（相角）。选择哪一组，取决于两侧的建模习惯、边界元件的物理特性以及数值稳定性要求。

电磁暂态侧通常以三相瞬时值描述网络，天然适合以电压、电流作为接口量；机电暂态侧以基波相量、有功/无功和转子角描述，天然适合以功率、角度作为接口量。若直接以瞬时值对接相量，需要额外的相量提取或瞬时值重构环节，会引入滤波延迟和幅相误差。因此，工程上常见的做法是：在EMT侧边界处提取基波正序电压相量与电流相量，再转换为功率和角度送入TS侧；反向则由TS侧给出边界电压幅值与相角，EMT侧将其重构为三相瞬时电压源。

从数值稳定性看，电压型接口（一侧给电压、另一侧回电流）通常比功率型接口更稳定，因为电压是“强加”量，电流是“响应”量，符合电路因果。功率接口虽然物理意义直观，但功率对电压和电流都是非线性乘积，容易在迭代中放大误差。角度接口则主要用于同步发电机等需要转子角信息的场景，单独使用角度作为接口量并不常见，通常与电压幅值配合构成“幅值—相角”接口。

二、交互时序

交互时序决定了两侧在什么时刻交换数据、以什么步长推进。EMT仿真步长通常在微秒级（如10～50 μs），TS仿真步长在毫秒级（如1～10 ms），两者相差两到三个数量级。若每个EMT步都与TS交互，通信开销巨大；若交互过疏，则边界动态被严重平滑。

常见的时序方案有三种。第一种是“TS主导、EMT跟随”：TS以自身步长推进，在每个TS步的起点向EMT提供边界电压，EMT在该步长内以多个小步积分，步末将边界电流的平均值或基波相量返回TS。这种方式实现简单，但EMT侧在步内看不到TS侧的变化，属于单向耦合。第二种是“EMT主导、TS插值”：EMT每步都推进，TS在较粗的时间网格上计算，中间时刻通过插值获得边界量。这种方式对TS侧动态的捕捉较弱。第三种是“并行迭代交互”：在一个交互窗口内，两侧分别求解，然后交换边界量并迭代，直到边界残差小于阈值。这种方式精度最高，但计算量和通信次数显著增加。

实际工程中，常采用“固定交互步长+步末平均”的折中方案：交互步长取TS步长的整数倍，例如1 ms；EMT在1 ms内以20 μs步长积分50步，步末将边界电流的基波相量送给TS。这样既控制了通信频率，又保留了EMT侧的主要动态。

三、数据交换格式

数据交换格式决定了接口量如何被编码、传输和解析。在离线协同仿真中，两侧可能运行在同一进程或不同进程，格式可以是内存数组、共享文件或标准化的数据流。在实时仿真中，则常通过共享内存、反射内存或以太网协议交换。

无论采用何种载体，接口数据至少应包含：时间戳、边界标识、变量类型、数值、单位、参考方向和质量标志。时间戳用于对齐两侧时间轴；参考方向必须统一，例如规定“电流从EMT侧流入TS侧为正”，否则功率符号会出错；质量标志用于标识该帧数据是否有效、是否经过插值或外推。

一个常见的交换帧结构如下（以文本格式示意，实际可用二进制）：

{
  "t": 0.001,
  "boundary": "BUS1",
  "vars": {
    "V_mag": 1.02,
    "V_ang": 0.05,
    "I_mag": 0.85,
    "I_ang": -0.12
  },
  "unit": {"V": "pu", "I": "pu", "ang": "rad"},
  "direction": "EMT_to_TS",
  "quality": "valid"
}

若两侧步长不同，还需在帧中注明该数据对应的时间窗口，例如“本帧电流为0～1 ms内的基波相量平均值”。缺少这一信息，TS侧就无法判断该数据的有效时间范围，容易造成时间错位。

四、接口误差来源与抑制

接口误差主要来自四个方面。第一是相量提取误差：EMT侧的瞬时值含有谐波和直流分量，若提取算法带宽不足或窗长不合适，得到的基波相量幅值和相角会有偏差。第二是时间离散误差：两侧步长不同，边界量在交互窗口内被平均或插值，高频动态被平滑。第三是参考方向与坐标变换误差：两侧坐标系、正方向定义不一致，或派克变换角度未对齐，会导致功率和角度出现系统性偏差。第四是迭代收敛误差：在并行迭代交互中，若迭代次数不足或松弛因子不当，边界残差未降到阈值以下就推进，误差会累积。

抑制措施包括：在EMT侧采用与TS侧基波频率匹配的滑动窗DFT或二阶广义积分器提取相量；在交互窗口内采用步末平均而非瞬时采样，减少混叠；在接口初始化阶段显式校验参考方向和坐标变换矩阵；对迭代交互设置残差阈值和最大迭代次数，并在超限时给出警告或降级为单向耦合。

五、一个可操作的接口配置示例

假设某测试系统在BUS1处划分EMT与TS子系统，EMT步长20 μs，TS步长1 ms，交互步长1 ms。接口变量选择为：TS→EMT传递电压幅值（pu）和相角（rad），EMT→TS传递基波电流幅值（pu）和相角（rad）。参考方向规定为“电流从EMT侧流入TS侧为正”。

初始化时，先由TS侧潮流计算得到BUS1的电压幅值和相角，送入EMT侧作为初始边界电压；EMT侧在第一个1 ms窗口内积分50步，步末提取基波电流相量，送入TS侧。TS侧用该电流更新边界注入功率，推进下一个1 ms。如此循环。判断接口是否正常，可观察：稳态下TS侧收到的电流相量与潮流计算值的偏差应小于0.5%；动态过程中边界功率的波形应与两侧内部功率波形趋势一致，无持续发散或高频振荡。

若发现边界功率出现周期性尖峰，优先检查相量提取窗长是否与基波周期匹配；若出现缓慢漂移，检查参考方向和坐标变换；若出现发散，检查交互步长是否过大或迭代是否收敛。通过这样一套“变量—时序—格式—误差”的闭环检查，接口问题通常可以定位到具体环节，而不是笼统地归因于“协同仿真不稳定”。