4.1 能量管理系统（EMS）与配网管理系统（DMS）

4.1 能量管理系统（EMS）与配网管理系统（DMS）

「Energy Management System」，能量管理系统——调度专业把这个词写在每间机房门口，但它的内核其实是一串算法：从原始量测里估出真实状态，再用状态去解潮流、查越限、发调节指令。本节作为决策中枢的开篇，负责把这套算法栈拆开讲透；上一章的数据由此完成从"原料"到"决策"的惊险一跃。

算法栈：四级台阶撑起调度台

EMS 的功能可以排成一条从底向上的台阶，每一级都在回答一个更难的问题。

第一级是数据采集与监视（SCADA）：遥测遥信上墙，越限告警，回答"系统现在报上来什么"。

第二级是状态估计：量测永远有误差、有丢失、有坏点，且量测数量少于全状态维数，需要用加权最小二乘法把成千上万条带噪量测"熔炼"成一组自洽的电压幅值与相角。这一步还兼职侦探——某条量测与全网估计明显冲突时，它会被标记为坏数据剔除。第 2 章说"感知端的质量决定闭环上限"，兑现之处就在这里：估计器的输出是后续一切决策的底图。

第三级是潮流计算与安全分析：拿到底图后，解出全网电流电压分布，逐条校核线路与变压器是否越限；再按 N-1 原则预演"任意一台设备退出后会怎样"，把预案做在事故前面。

第四级是闭环控制：自动发电控制（AGC）管频率与联络线功率，自动电压控制（AVC）管无功与电压，指令周期性下发到电厂与变电站执行。

配电网侧的对应物是配网管理系统（DMS）。配网的三相不平衡、辐射状拓扑、量测稀疏使它无法照搬输电算法，但问题同构：状态估计换成了三相与伪量测版本，安全分析换成了故障定位、隔离与负荷转供（俗称 FLISR），外加一个输电网没有的特色业务——馈线重构，通过改变分段开关的分合状态均衡负荷、降低网损。第 5 章讲配网自动化时，DMS 就是那只"看得见的手"。

AGC：一个每四秒转一圈的闭环

AGC 值得单独讲，因为它是调度里最像"反射"的部分。电网频率反映发电与负荷的瞬时平衡：负荷突然增加，频率立刻下跌，机组调速器在数秒内按下垂特性自动增出力——这叫一次调频，快而粗糙；AGC 在此之上每几秒钟计算一次区域控制偏差（ACE，综合频率偏差与联络线计划偏差），向调频机组下发功率修正指令，把频率拉回标准值并守住与其他区域的交换计划——这叫二次调频，慢而精准。再往上是小时级的经济调度与机组组合，负责"让谁发电最省钱"，属于三次调优。

三层结构的时间尺度差了三个数量级，这正是第 1 章"闭环分层"论断的现场证据。

图 4-1 频率调节的三层时间尺度

演练：亲手解一次三节点潮流

潮流计算是 EMS 的心脏。它的数学面貌是：已知网络参数与边界条件，求各节点电压，进而得到全部支路潮流。方程非线性，工程上用牛顿-拉夫逊法迭代求解——线性化、解修正方程、更新、重复，直到失配量足够小。下面是一个三节点系统（1 号平衡节点、2 号带负荷、3 号带电源）的完整求解：

import cmath, math

Yc = [[-9j, 5j, 4j], [5j, -13j, 8j], [4j, 8j, -12j]]   # 节点导纳矩阵（纯电抗线路）
Bm = [[-9, 5, 4], [5, -13, 8], [4, 8, -12]]            # 其虚部，用于极坐标功率式
P_spec = {2: -0.8, 3: 1.0}                              # 2 号带负荷、3 号带电源（标幺值）
Q_spec = {2: -0.4, 3: 0.5}

def pq(delta, vm):                                      # 计算节点 2、3 的注入功率
    out = {}
    for i in (2, 3):
        p = q = 0.0
        for j in (1, 2, 3):
            ang = delta[i] - delta[j]
            p += vm[i] * vm[j] * Bm[i-1][j-1] * math.sin(ang)
            q -= vm[i] * vm[j] * Bm[i-1][j-1] * math.cos(ang)
        out[i] = (p, q)
    return out

x = {2: [0.0, 1.0], 3: [0.0, 1.0]}   # 每个待求节点存 [相角, 幅值]，平启动

def states():
    delta = {1: 0.0, 2: x[2][0], 3: x[3][0]}
    vm    = {1: 1.0, 2: x[2][1], 3: x[3][1]}
    return delta, vm

def solve4(J, b):                    # 4x4 高斯消元（含列主元）
    A = [row[:] + [b[i]] for i, row in enumerate(J)]
    for c in range(4):
        p = max(range(c, 4), key=lambda r: abs(A[r][c]))
        A[c], A[p] = A[p], A[c]
        for r in range(4):
            if r != c and A[c][c] != 0:
                f = A[r][c] / A[c][c]
                A[r] = [A[r][k] - f * A[c][k] for k in range(5)]
    return [A[i][4] / A[i][i] for i in range(4)]

eps = 1e-6
for it in range(1, 12):
    delta, vm = states()
    res = pq(delta, vm)
    F = [P_spec[2]-res[2][0], Q_spec[2]-res[2][1], P_spec[3]-res[3][0], Q_spec[3]-res[3][1]]
    if max(abs(f) for f in F) < 1e-10:
        print(f"第 {it-1} 次迭代收敛"); break
    vec = [x[2][0], x[2][1], x[3][0], x[3][1]]
    J = []                            # 数值雅可比：逐列扰动求差分
    for c in range(4):
        vp = vec[:]; vp[c] += eps
        dp, mp = {1:0.0, 2:vp[0], 3:vp[2]}, {1:1.0, 2:vp[1], 3:vp[3]}
        rp = pq(dp, mp)
        Fp = [P_spec[2]-rp[2][0], Q_spec[2]-rp[2][1], P_spec[3]-rp[3][0], Q_spec[3]-rp[3][1]]
        J.append([(Fp[k]-F[k])/eps for k in range(4)])
    JT = [[J[c][k] for c in range(4)] for k in range(4)]
    dx = solve4(JT, [-f for f in F])  # 牛顿修正量沿失配的负梯度方向
    for k, key in enumerate((2, 2, 3, 3)):
        x[key][k % 2] += dx[k]

print(f"V2 = {x[2][1]:.4f} ∠ {math.degrees(x[2][0]):.2f}°")
print(f"V3 = {x[3][1]:.4f} ∠ {math.degrees(x[3][0]):.2f}°")
Vc = [1+0j, x[2][1]*cmath.exp(1j*x[2][0]), x[3][1]*cmath.exp(1j*x[3][0])]
S1 = Vc[0] * sum(Yc[0][j]*Vc[j] for j in range(3)).conjugate()
print(f"平衡节点注入 S1 = {S1.real:.3f} + j{S1.imag:.3f}")
​

程序第三次迭代即收敛：2 号节点电压约 0.981∠负一度（负荷侧电压略低、相角滞后），3 号节点约 1.025∠正四度（电源侧相角超前），平衡节点补入约负零点二的有功——正好补齐负荷与电源之间的差额，无功平衡与网损也完全自洽。读者可以自己验证两个性质：把 2 号负荷加大一倍再看电压跌落；把 3 号电源出力归零，观察潮流如何反灌。调度的日常，就是把这个三节点演算放大到上万节点，并且每五分钟算一遍。

运行经验：算法栈的三个现实

其一，坏数据比缺数据更危险：状态估计能容忍少量量测缺失，但一条被污染的关键量测若未剔除，会让整张底图失真，现场要求对估计残差做常态化监视。其二，模型与实际必须对表：线路参数漂移、拓扑图实不符都会让潮流"算得出来但不对"，参数辨识与图模校核是调度自动化的隐形战场。其三，配网算法要接受稀疏：DMS 的状态估计大量依赖伪量测与历史数据补位，精度目标应设为"够用且可解释"，追求输电级精度只会得到虚假的自信。

本节要点

• EMS 算法栈四级台阶：采集、估计、分析、控制，每一级都在把"原始数据"提炼成更高层级的"判断"。
• AGC 的三层时间尺度互不越位，分层解耦是电网控制区别于一般工业控制的关键特征。
• 牛顿-拉夫逊的"线性化-修正-迭代"套路值得亲手跑一遍：它是调度台每五分钟上演一次的世界建模。