做电力系统二次这一行最吊诡的一个场景，我愿称之为“薛定谔的电表”：调度中心屏幕上那个节点电压，你不去看它的时候，总觉得它应该是稳的；你真的去查它，发现要么量测点根本不在那儿，要么那个测点的通道已经坏了一个星期，要么它的遥测精度让人怀疑人生。系统越大，这种“既存在又不存在”的量测就越扎眼。
这套题有个正经解法叫MMSE，最小均方误差估计。我不是说它是万能的，但在“量测数量不够、质量还拉胯”的有限测量困局里，它确实比传统加权最小二乘更抗造。这篇博文就是聊聊我实际做状态估计、做CDT规约解析、被RTU和调度主站来回折腾时，怎么理解MMSE，怎么把它落地到电力系统自动化工程里。适合刚接触状态估计的新人，也适合正在跟有限量测死磕的现场工程师。
1. 先聊聊“薛定谔的电表”：有限测量困局到底被困在哪
1.1 量测永远不完整：一个真实的二次系统视角
先说一个反常识的事实：调度中心看到的电网，从来不是完整电网。你以为SCADA系统里有全网的实时数据，实际上变电站侧受投资、通道带宽、设备老旧程度影响，很多节点压根没有装高精度测量装置。就算装了，RTU的上送周期、死区设置、通信规约的帧格式，都会让部分数据要么过时、要么缺失、要么干脆是上个周期的“死值”。
我见过很多地区级的调度主站，状态估计覆盖的节点可能就占全网的六七成。剩下那些节点，既没有遥测也没有遥信，基本处于“盲区”。你让调度员盲调？肯定不行。所以有限测量困局的第一层意思，就是物理层面上测量设备的稀缺和分布不均。电网那么大，不可能给每个母线和线路都配上高精度的同步测量单元，那成本谁也扛不住。
第二层意思，是数据质量维度的“有效测量”比名义测量更少。一个RTU通道里出来的数据，经过CDT或者101、104规约解析以后，能不能直接进状态估计？能不能直接进MMSE？真不一定。坏数据、死数据、时标错乱，这些光是做检测就是一堆活。所以很多时候，真正的有效量测数量比SCADA里看到的遥测总数低一截。
第三层，也是最容易被忽略的，是量测方程本身的非线性。电力系统的量测模型是电压幅值、相角和支路潮流的非线性函数。节点多了、网络复杂了，量测配置稍微稀疏一点，传统求解方法就可能遇到雅可比矩阵病态甚至秩亏的问题。量测一少，系统就算不过去；就算算过去，结果也可能离谱得没法看。
1.2 经典方法为什么会在稀疏量测下翻车
传统电力系统状态估计最常用的主角是加权最小二乘，也就是WLS。WLS的思路很简单：找一个状态向量，让量测残差的加权平方和最小。这在量测冗余度足够、数据质量尚可的时候，表现得相当不错。
但WLS有个硬前提：量测方程的数量要大于等于待估计的状态量个数，而且量测矩阵要满秩。一旦量测配置不满足这个条件，雅可比矩阵会奇异或者接近奇异，方程要么无解，要么数值上极不稳定。你拿一套稀疏量测数据去跑WLS，经常会出现估计值发散、迭代不收敛，或者收敛到一堆火星值的情况。
更麻烦的是，WLS没有把“我们对系统已有的了解”用起来。比如说，某条母线的负荷虽然没量测，但调度运行人员心里清楚它大概是多少，历史上这个时段它就在这个区间附近波动。WLS对这种信息直接无视，全靠量测硬推；量测一缺，它就抓瞎。这就像让一个考生闭卷考试但不给草稿纸，对知识储备和题目覆盖度要求极高。
所以稀疏量测场景下的核心矛盾就出来了：系统真实状态只有一个，但量测远远不足以把它唯一确定下来。这时候你要么接受“无解”，要么换一种思路——把“不确定性”本身变成数学上的资源，用统计方法从概率意义上把状态“估计”出来。这就是MMSE登场的原因。
2. 从最小二乘到MMSE：为什么偏偏是它
2.1 先看传统状态估计到底怎么干的
为了说清MMSE的价值，简单回顾一下WLS的套路是值得的。电力系统状态估计选择的状态变量通常是各节点的电压幅值V和相角θ，量测向量z包括节点注入功率、支路潮流、电压幅值等。量测方程写成：
z = h(x) + v
其中v是量测噪声，一般假设零均值高斯噪声，协方差矩阵为R。WLS的目标函数是min [z - h(x)]ᵀ R⁻¹ [z - h(x)]，通过对h(x)线性化，用高斯-牛顿法迭代求解：
x_{k+1} = x_k + (Hᵀ R⁻¹ H)⁻¹ Hᵀ R⁻¹ [z - h(x_k)]
这里H是量测雅可比矩阵。看到那个(Hᵀ R⁻¹ H)⁻¹没有？它要求Hᵀ R⁻¹ H可逆。量测不足或者系统弱可观时，这个矩阵就是奇异的。我在实际运行里见过最典型的现象，就是某几个边远节点的电压幅值估计结果在连续几个调度周期里来回蹦，幅度大到一度以为系统真的在振荡，查了半天发现是量测矩阵条件数已经上千了，纯粹数值问题。
2.2 MMSE的核心：不把未知当敌人，把未知变成先验
MMSE的全称是Minimum Mean Square Error，最小均方误差估计。它的目标不是找一个让残差平方和最小的点，而是找一个让估计误差的期望平方和最小的估计量，也就是最小化：
E[ ||x - x̂||² ]
这个目标函数里出现了一个WLS里没有的东西——数学期望。这意味着MMSE要求我们把状态量x本身看成随机变量，并且给它一个分布。听起来抽象，但工程上这个分布恰恰来自我们手里那些“没人要”的常识和历史数据。
在线性高斯模型下，MMSE有一个非常优美的闭合解：估计值等于在量测条件下的条件期望。如果已知先验均值x̄、先验协方差P₀，量测噪声协方差为R，量测矩阵为H，那么：
x̂ = x̄ + P₀Hᵀ (HP₀Hᵀ + R)⁻¹ (z - Hx̄)
看一眼这个公式的形态，本质就是“先验估计 + 增益 × (量测残差)”。当量测不够时，先验项x̄还在起作用；当量测充足高质量时，增益会把估计推向量测。这种“两边下注”的结构，天生就是给稀疏量测准备的。
用个生活化的类比：你在一间几乎全黑的屋子里找一把椅子。WLS的思路是只靠手往前摸，摸到桌子角也算一个约束，摸不到就抓瞎。MMSE的思路是你记得这个房间大概布局，椅子七成在靠窗那个位置，然后你再用手去摸，摸到一点线索就和记忆融合，越摸越准。没有记忆的人摸黑找椅子，和带着模糊记忆的人摸黑找椅子，后者显然更稳，而MMSE干的就是这件事。
2.3 MMSE和WLS到底是什么关系
有个关键点很多人容易忽略：MMSE不是要彻底推翻WLS。在没有任何先验信息的极端情况下，把先验协方差P₀推向无穷大，MMSE的解就会退化成WLS的形式。反过来说，WLS只是MMSE家族里一个“无先验”的特例。
这个关系的工程意义很大。我在做方案的时候不会跟老工程师说“换一套方法”，而是说“在原有状态估计框架上引入先验融合”。这样项目评审好过、被试人接受度也高。而且卡尔曼滤波本身就是递归形式的MMSE，你平时可能没意识到，但扩展卡尔曼滤波在电力系统动态状态估计里的应用，本质上就是MMSE思想的工程化延伸。
| 对比维度 | 传统WLS状态估计 | MMSE状态估计 |
|---|---|---|
| 目标函数 | 量测残差加权平方和最小 | 估计误差期望平方和最小 |
| 对量测数量要求 | 量测冗余，否则矩阵奇异 | 稀疏量测也可用，先验兜底 |
| 对先验信息利用 | 完全不用 | 融合历史/预测/拓扑信息 |
| 与卡尔曼滤波关系 | 无直接关联 | 卡尔曼滤波是递推MMSE |
| 工程落地难度 | 实现简单但易发散 | 需要先验建模，实现稍复杂 |
| 抗稀疏量测能力 | 弱 | 强 |
3. 工程落地：MMSE状态估计的建模与实现细节
3.1 量测方程：先把你手里有什么牌列清楚
落地MMSE，第一步不是调公式，是理清楚你手里到底有哪些量测。按电力系统自动化模板的说法，量测体系一般分三类：模拟量遥测、状态量遥信、电能量。状态估计主要吃遥测，包括节点注入有功无功、支路有功无功潮流、节点电压幅值；少数高级应用还会吃相角量测。
每个量测都要对应一个量测方程。比如支路i-j的有功潮流，方程是：
P_ij = Vi² g_ij - Vi Vj (g_ij cosθ_ij + b_ij sinθ_ij)
节点注入有功就是相连支路潮流的求和。这些方程是物理世界到量测空间的映射，MMSE也好、WLS也好，都绕不开它。我在实际项目中习惯先把量测清单列成表，字段包括：量测点编号、量测类型、数值、量测方差、时标、通道状态。这张表既是数据预处理的对象，也是状态估计的输入档案。
量测方差的设定也大有讲究。很多工程人员直接拍脑袋给个0.01或者0.001，这在量测足够时看不出毛病，但稀疏场景下方差大小直接决定量测和先验之间谁说了算。我常说，方差是状态估计里的“信任权重”，给错了就等于裁判瞎吹黑哨。通常的做法是结合互感器精度等级和A/D转换量化误差估算，比如0.2级的电压互感器，量测标准差可以取额定值的0.2%到0.5%；没把握时宁可取大一点，让估计器别太自信。
3.2 先验信息从哪来：三大实用来源
MMSE说要有先验，可先验不是天上掉下来的。我在工程里总结了三类最靠谱的来源。
第一类是负荷预测信息。调度中心一般都有短期负荷预测系统，能给出未来时段各片区或各主变的负荷期望值。把负荷预测折算成功率注入的期望值，就成了节点注入功率的先验均值。误差范围可以从预测系统的历史精度统计出来，通常10%到20%的相对误差算是常见水平，这就对应先验协方差P₀的对角元素。
第二类是拓扑约束和物理限值。一个变压器节点的功率不会超出额定容量太多，一台发电机的有功出力受运行上下限约束。这些物理边界虽然不是精确的先验均值，但可以作为强约束或软约束揉进先验协方差。我在做220kV变电站估计时，就经常用主变容量来限制高压侧注入功率的取值范围，效果立竿见影，估计值再也没出现过“主变容量翻三倍”的闹剧。
第三类是历史相似日数据。把过去几个月同类型日、同时段的历史遥测求均值和方差，作为当前时段的先验。这个方法在负荷规律性强的地区非常有效。比如工业园区周一到周五的负荷曲线重复性很高，用前四周同工作日同一时刻的数据做先验，节点注入的估计精度能有明显提升。
P₀的建模也要谨慎。先验协方差矩阵的维度等于状态量维度，在实际电网里可能是几千乘几千的大矩阵。但好在拓扑上相距很远的节点之间先验相关性往往很弱，可以近似按区域分块处理。这一招能把矩阵运算量降一个量级。
3.3 非线性扩展与实现步骤
前面给的MMSE公式是线性的，电力系统量测方程是非线性的，所以不能拿公式直接套。工程上常见做法是照着扩展卡尔曼滤波的思路走：在当前估计点附近对量测方程做线性化，得到雅可比矩阵H，然后用线性MMSE的框架迭代更新。
具体迭代流程大概是：
- 初始化先验均值x̄和先验协方差P₀，读取当前量测向量z及其噪声协方差R；
- 计算在当前状态x̂下的量测函数h(x̂)，以及它的雅可比矩阵H；
- 计算新息矩阵S = HP₀Hᵀ + R，求增益矩阵K = P₀Hᵀ S⁻¹；
- 更新状态估计：x̂_new = x̂ + K [z - h(x̂)]；
- 更新估计协方差：P_new = (I - KH)P₀；
- 检查收敛条件，不满足就回到第2步，用x̂_new替代x̂。
用伪代码表示，核心逻辑不到二十行就能说清：
读入先验均值 x_prior, 先验协方差 P, 量测 z, 噪声协方差 R x = x_prior 循环直到收敛： 计算 h = H_function(x), H = Jacobian(x) 新息：S = H * P * H^T + R 增益：K = P * H^T * inv(S) 更新：x_new = x + K * (z - h) P = (I - K * H) * P x = x_new 输出 x 和 P我在项目里直接用Python的NumPy就能跑这套迭代，算一个几百节点的区域网络，单次迭代性能完全够在几秒内出结果。真正要花心思的反而是数据接入和矩阵整形，这些过程性代码的量反而是算法本身的好几倍。
3.4 一个三节点小系统的数字演练
为了说明MMSE到底怎么化腐朽为神奇，我搭一个极简三节点系统手算一遍流程。假设系统里节点1是平衡节点，电压幅值固定1.0，相角0；节点2和节点3是负荷节点。拓扑上1-2和2-3各有支路，线路参数是标幺值，电抗都是0.1，电阻忽略。
量测配置很寒酸：只有支路1-2的有功潮流P12 = 2.0，以及节点2的电压幅值V2 = 0.98，其他全部缺失。这样两个量测要估计的节点2幅值、节点2相角、节点3幅值、节点3相角，明显方程个数不够。传统WLS此刻已经宣告阵亡，矩阵秩亏跑不起来。
MMSE的做法是用先验补：假设根据历史负荷预测，节点2和节点3的注入有功期望分别是-1.5和-1.0，标准差20%。把这些转换成功率注入方程的先验约束，加上量测方程，形成增广估计。迭代几步后，节点3的相角被合理地估计出来，而不是因为没有量测就卡死。关键在于，这种估计结果不会出现节点3功率反向这种物理上荒谬的结论，因为先验协方差把它拽在合理范围内。
这个三节点例子的意义不是精确数字，而是展示一种范式：当量测不够时，MMSE用概率约束替代硬性量测，把不可观系统变成了弱可观系统，再把弱可观变成可估计。这是我实际调试中觉得它最值钱的地方。
4. 喂给MMSE的数据从哪来：CDT规约与数据质量实战
4.1 CDT规约解析是绕不开的基本功
聊MMSE的上游，必然要聊数据来源。国内老一代变电站自动化设备大量采用CDT规约，也就是循环式远动规约。它不像IEC 60870-5-104那种问答式，而是RTU按固定周期主动循环上送数据。主站被动接收，所以叫“循环式”。
CDT规约的帧结构一般包括同步字、控制字、信息字。同步字用EB90H开头，控制字里包含帧类别、信息字数等，信息字再根据功能码区分遥测、遥信、电度量和事件顺序记录。解析CDT规约，基本的功夫是先找同步字，确认帧同步，再校验长度和校验码，最后按信息字地址把数据还原成物理量。
我在现场调试中最常用的一张表就是功能码映射表：哪个功能码对应哪个遥测点，哪个地址对应哪路遥信。很多数据“丢了”根本不是规约问题，是主站侧转发表和RTU侧的信息字地址没对齐。CDT规约解析看着是个低技术含量的活，但恰恰是它决定了MMSE吃到嘴里的是干净饭还是馊饭。
4.2 数据质量的三种坑：死数据、坏数据、时标错位
第一坑是死数据。RTU的采集模块如果卡死，或者通道中断但主站没发现，CDT帧照样按照周期发过来，但里面的数值是某个历史时刻的快照。主站如果直接把死数据送进状态估计，MMSE会以为量测和先验冲突剧烈，结果像抽风一样来回摆动。判断死数据的方法很简单：连续多个周期数值纹丝不动，超过正常波动范围的时间窗，就要打上“疑似死数据”的标签。
第二坑是坏数据，比如A/D采样芯片受温度漂移影响、互感器饱和导致大电流工况下数据失真，或者通信链路干扰导致某位翻转。坏数据在CDT解析环节往往能通过校验和发现一部分，但校验通过的数据也可能被干扰引入合理范围内的偏差。这类数据混进MMSE后，因为MMSE有先验兜底，单点坏数据的影响会被摊薄，但多个坏数据协同作用时还是会把估计值带跑。
第三坑是时标错位。CDT规约循环上送的数据本身就是“准实时”，它没有高精度的时标字段。如果解析时不做缓存对齐，有的数据是这秒的、有的是上一周期的，混在一起做状态估计，相当于把不同时刻的系统状态拼在一起算，结果自然失真。我的处理原则是：给每帧数据打上接收时间戳，再按通道延时估算进行粗补偿，宁愿数据新鲜度差一点，也不能让学生吃夹生饭。
4.3 实用的数据预处理技巧
有了这三类坑，数据预处理就有了方向。我的做法是在MMSE前面加一个三兄弟模块：
第一兄弟是合理性检查。每个遥测必须落在上下限区间内，比如电压幅值在0.9到1.1倍额定电压之间，功率不超过线路热稳极限。超限直接标记可疑，不参与估计或降权处理。
第二兄弟是突变检测。对比当前值与前后几个周期的值，如果变化幅度超过物理允许的最大爬坡速率，比如一秒钟内主变负荷翻倍，基本可以断定是坏数据。检测出来以后不删除，而是把协方差R对应元素调大，把它从“可信量测”降级为“模糊信息”。这个思路跟MMSE特别搭，直接把数据质量问题映射成了噪声参数问题。
第三兄弟是零值处理。电网设备停运或者通道故障时遥测会变0，但这种0不代表真实潮流为0，处理不当会让估计器把系统算成“全网停电”。实际做法是结合遥信位置判断设备状态，遥控信为合位但遥测为0时，按死数据处理。
5. 现场工程中的常见问题与排查速查表
5.1 先验设错会怎样
先验是MMSE的双刃剑。先验均值设偏了，比如负荷预测系统因为天气突变给出严重偏离实际的注入功率期望，MMSE会把估计结果往错误方向拽。量测充足时它还能靠量测纠正，量测一稀疏，误差就会一直存在。我在一个项目里就吃过这个亏，节点3的估计功率连续一周都比SCADA后来核对的实际值低15%，查到最后是负荷预测模型没有更新节假日模板。
解决办法是先验自适应。把估计残差的统计量实时反馈给先验模块：如果残差均值持续不为零，说明先验有偏，要动态修正先验均值；如果残差方差比理论方差大很多，说明先验协方差给小了，需要放大P₀。说白了，MMSE的参数不能设完就撒手不管，它需要一套自动调参的闭环。
5.2 数值稳定性与计算性能的坑
稀疏量测下，增益矩阵K的计算涉及对S = HP₀Hᵀ + R求逆。当P₀过于自信，或者H的某些行与其他行高度相关时，S可能接近奇异，求逆结果就会出现巨大的数值波动。我也翻过车，表现为估计值在某个迭代步突然跳到几十亿，那纯粹是S矩阵条件数太差。
应对手段首先是加正则化，在S对角线上加一个小量εI，工程里ε取1e-6到1e-8量级基本管用。其次是换更稳的分解方式，比如Cholesky分解或奇异值分解替代直接求逆。最后是把模型分块，每个区域独立做MMSE，然后边界协调。这样既能提高数值稳定性，也能减少单次矩阵运算量，实时性更好。
5.3 常见问题速查表
| 故障现象 | 可能原因 | 处理办法 |
|---|---|---|
| 估计值持续跳变 | 死数据未处理、量测方差过小 | 加死数据检测，重新标定量测方差 |
| 估计值整体偏移 | 先验均值有偏、负荷预测未更新 | 启动残差反馈，修正先验均值 |
| 迭代不收敛 | 雅可比矩阵病态、S接近奇异 | 加正则化，改用SVD或分块求解 |
| 某节点估计结果荒诞 | 先验协方差给太小、量测通道反向 | 检查该节点拓扑和先验建模 |
| 数据时标混乱 | CDT帧未打时间戳、通道延时未补偿 | 统一接收时间戳，按通道补偿 |
| 估计结果对量测噪声过敏感 | 噪声协方差R设置过小 | 按互感器精度重新估算R |
5.4 一个普适的电力系统自动化落地模板
最后给一个可以抄作业的实现模板，这是我做多个项目后固定下来的流程框架，从原始数据到MMSE估计输出共七步：
第一步数据接入，通过CDT解析库把RTU数据变成内存中的量测表。第二步数据清洗，执行合理性检查、突变检测、零值处理，生成带质量标签的干净量测集。第三步量测聚合，把分散的遥测按节点和支路归类，形成量测向量和量测矩阵。第四步先验构建，从负荷预测、拓扑库、历史库生成当前时段的先验均值和协方差矩阵。第五步参数配置，设定量测噪声协方差R、收敛门限。第六步迭代估计，按扩展MMSE流程循环更新状态。第七步结果校验与输出，对比估计值与已知开关状态、限值，输出到调度自动化画面和后续高级应用。
这套模板最大的好处是每一层都能单独调试。哪个阶段出了问题，日志里能看得很清楚。比如“量测聚合”阶段发现信息字地址对不上，就不用去怀疑迭代算法，排查范围直接从整个系统缩小到一个模块。
做状态估计这几年，我最深的体会是：电网调度不缺理论最优的算法，缺的是能在有限量测、有限数据质量、有限算力条件下稳定输出结果的那一套东西。MMSE的妙处在于它把“不知道”和“知道一点”都摆到桌面上来算，而不是假装问题不存在。先验不是胡猜，是对系统运行规律的事后总结；量测不是全知，是对真实状态的有噪观测。两者一融合，薛定谔的电表虽然还戴着面纱，但至少你不再是对着全黑房间干瞪眼的那个人了。
