高铁牵引供电系统运行数据分析及等值建模（2020 电工杯 A 题）

历年真题解法库 · 本题为公开赛题要点整理，完整原题请到官方渠道获取（链接见下）；
下方「练习数据」为本站生成的合成数据集，仅用于跑通方法，非官方原题数据。

一、赛题要点

• 赛事：2020 年"电工杯"全国大学生电工数学建模竞赛 · A 题（真实题名"高铁牵引供电系统运行数据分析及等值建模"）
• 问题实质：高铁牵引变电所把 220 kV 三相电压经牵引变压器变换为 27.5 kV 单相电压，经馈线、接触网向动车组供电，再由钢轨与回流线回流。这类负荷单相、大功率、剧烈波动且非线性，会同时引入负序不平衡、谐波注入与冲击功率三类问题；同时列车再生制动的能量若无处消纳就会被制动电阻烧掉。赛题给出该变电所 20 kHz 采样的三相电压电流录波（空载、牵引、制动三个 4000 点窗口 + 全程约 200 s）与全天 86400 点功率数据，要求依次回答四问：①评估三种典型工况及全程数据的电能质量；②提出节能降耗措施并量化收益；③对动态牵引负荷做精确预测；④建立牵引供电系统等值数学模型并校验其与实测的偏差。

二、需要产出什么

1. 电能质量评估：对称分量法分解正序/负序/零序，算电流与电压不平衡度；整周期 DFT 求各次谐波含有率与总谐波畸变率 THD；给出瞬时有功、无功与功率因数曲线。
2. 节能降耗：全天曲线正负分离积分求牵引电量与再生制动电量；设计车载储能与变电所地面储能两套回收方案，比投资、年节约电费、静态回收期与节电率。
3. 动态负荷预测：从全天曲线中提取牵引脉冲，按编组（8/16 节）、方向（上/下行）、车型（CRH380B/CRH3C/CRH2A）分类辨识，回归峰值功率模型，再按时刻表外推次日负荷。
4. 等值建模：工频层次做戴维南等值（等值电势 + 内阻抗），谐波层次做集中 RLC 等效并算谐波阻抗与串联谐振点，最后把模型与实测的偏差分解归因。

三、数据说明（官方来源 + 练习数据集）

• 官方数据/赛题：中国电机工程学会 / 电工杯组委会当年官网发布的题面与附件（附件为多个 20 kHz 录波文件与全天功率文件），请以当年官方发布为准。真实建模应直接读入官方录波。
• 官方获取渠道：http://www.mcm.edu.cn（电工杯以组委会当年通知为准）
• 本站练习数据集：data/dgcup2020a.csv —— 由真源模型 tools/gen_dgcup2020a.py 产出的 三工况电能质量 + 全程录波 + 全天电量 + 两方案经济性 + 负荷辨识 + 等值参数 面板（84 行"参数—数值—单位—说明"）。
• 合成口径：脚本以 random.seed(2020) 固定种子、纯标准库实现，严格沿用赛题给定的采样率 20 kHz、窗口长度 4000 点（= 0.2 s = 10 个工频周期）与窗口位置（[1,4000] / [1600000,1604000] / [4000000,4004000]），全程 200 s、全天 86400 点。数值量级贴近实际牵引变电所，但是确定性合成、非官方实测。
• 字段说明：CSV 数值全部由真源脚本一次算出，与三篇范文的正文、配图、附录四路一致，可直接复现与练手。

四、可用解法 → 对应手册

本题横跨信号与频谱分析 + 预测 + 回归/分类 + 优化决策 + 评价。下面按题型给出对应手册跳转。

题型：预测类（动态牵引负荷）

先看选型 👉 🧭 预测类选型指南 · 🧭 时序类选型指南

• 📘 多元线性回归
• 📘 一元线性回归
• 📘 时间序列 ARIMA
• 📘 灰色预测 GM(1,1)
• 📘 BP 神经网络

题型：分类 / 判别类（编组、方向、车型辨识）

先看选型 👉 🧭 分类类选型指南

• 📘 K-Means 聚类
• 📘 K 近邻 KNN
• 📘 主成分分析 PCA
• 📘 支持向量机 SVM

题型：优化类（储能容量与功率配置）

先看选型 👉 🧭 优化类选型指南

• 📘 非线性规划
• 📘 多目标规划
• 📘 遗传算法

题型：评价 / 综合指数类（方案比选）

先看选型 👉 🧭 评价类选型指南

• 📘 熵权法
• 📘 TOPSIS（两方案排序）

题型：不确定性 / 稳健类（电价与补偿容量灵敏度）

先看选型 👉 🧭 不确定性选型指南

• 📘 蒙特卡洛仿真

五、建模思路提示

• 分窗必须取整周期：20 kHz 下 4000 点正好是 0.2 s、10 个工频周期，整周期 DFT 无频谱泄漏；若随意截断 3000 点，谐波幅值会出现数个百分点的虚假误差。
• 对称分量法：正序 = 三相相量按 1、a、a² 加权平均，负序把 a 与 a² 互换，零序取三相之和的三分之一（a 为 120° 旋转因子）。V/v 接线三相侧无零序通路，算出的零序应接近 0，可当作程序自检。
• 区分两个不平衡度：电流不平衡度在空载时会因正序基数极小而被放大到 40% 以上，不能据此判越限；真正考核并网影响的是电压不平衡度（限值 2%）。
• THD 与负荷电流反向：牵引时基波大、谐波被"稀释"，THD 最低；空载时基波仅几安培，励磁谐波占比被放大，THD 最高但绝对危害小。
• 负荷辨识用物理判据而非黑箱：16 编组功率约为 8 编组的 1.92 倍，可用峰值阈值一刀切；上行（上坡）牵引附加功率大而可回馈动能少，故"制动峰值/牵引峰值"之比能干净地分出方向。回归模型务必带编组 × 车型交互项，否则车型影响在长编组上会被系统性低估。
• 等值建模要分层：工频层次拟合 U = E0 − Zeq·I 的下垂直线；谐波层次用系统电感与补偿电容的串联回路算 Z(h)，谐振次数与补偿电容平方根成反比。
• 让误差可归因：把模型误差拆成"谐波截断"与"工频压降"两项，就能直接判断该增加保留谐波次数还是该把集中阻抗改成分段阻抗。

六、常见坑与规避

• 单位陷阱：逐秒功率（MW）积分成电量要除以 3.6 才是 kWh；不要把 86400 点直接当小时数用。
• 窗口位置别搞错：制动窗口起点 4000000 对应 200 s，恰是全程录波的末尾，因此制动功率是全程最深的负值点；若误按 4000 s 定位会与全程时长矛盾。
• 功率因数的符号：制动工况有功为负、功率因数为负，是能量回馈的直接证据，不能取绝对值。
• 再生能量不能全算收益：必须乘"可利用份额"与"储能综合效率"两道折扣，否则收益会被高估三成以上。
• 回收期别只算一个电价：做 0.45~0.85 元/kWh 的区间扫描，确认方案排序不反转再下结论。
• 谐振点必须校核：补偿电容容量变化会让谐振次数穿越 11 次与 13 次两个危险点，选容量时要与滤波器调谐协同。
• 合成练习数据只验证逻辑，正式参赛务必用官方录波并做敏感性与误差分析。

七、优秀范文（本站手写，非生成器）

三篇均为 ≥2500 字、含可独立运行复现附录的优秀范文，正文 / 配图 / 附录 / 真源四路数字一致。

视角	链接	字数 · 配图 · 附录
一、电能质量评估（序分量 / 谐波 / 冲击功率）	dgcup2020a-1.md	手写 · 8 图 · 1 附录
二、节能降耗与动态牵引负荷预测	dgcup2020a-2.md	手写 · 8 图 · 1 附录
三、等值建模、谐波阻抗与误差分解	dgcup2020a-3.md	手写 · 8 图 · 1 附录

权威数字速览（四路一致，已验证）

• 问题一 · 三工况电能质量（窗口均为 4000 点 = 0.2 s = 10 个工频周期）：
  • 空载（[1,4000]，t=0 s）：正序 0.390 A、负序 0.162 A、零序 0.0000 A，电流不平衡度 41.50%、电压不平衡度 0.003%；馈线电流 5.4 A、THD 18.40%（3 次 13.34% / 5 次 9.01% / 7 次 6.98%）；有功 0.013 MW、无功 0.148 Mvar、功率因数 0.086。
  • 牵引（[1600000,1604000]，t=80 s）：正序 21.795 A、负序 6.190 A，电流不平衡度 28.40%、电压不平衡度 0.740%（三工况最大，限值 2%）；馈线电流 302.0 A、THD 5.73%；视在 8.305 MVA、有功 7.989 MW、无功 2.268 Mvar、功率因数 0.962。
  • 制动（[4000000,4004000]，t=200 s）：正序 7.794 A、负序 2.541 A，电流不平衡度 32.60%、电压不平衡度 0.257%；馈线电流 108.0 A、THD 8.36%；有功 −2.795 MW、功率因数 −0.941（能量回馈）。
  • 全程 200 s：最大牵引功率 8.971 MW（第 62 s）、最大制动功率 −2.818 MW、最大馈线电流 339.8 A、最低母线电压 27.016 kV；窗口内牵引 172.40 kWh、再生 26.75 kWh。
• 问题二 · 节能降耗：全天牵引电量 35889.4 kWh、再生制动电量 6515.3 kWh（占比 18.15%），年电费基数 825.3 万元；现状邻车吸收率 38.6%，每日弃用 4000.4 kWh。
  • 方案 A 车载超级电容（24 列 × 22.5 kWh）：效率 0.912、份额 0.742，日回收 4408.9 kWh，投资 636.0 万元，年节约 101.4 万元，回收期 6.27 年，节电率 12.28%。
  • 方案 B 变电所地面储能（640 kWh + 4 MW 变流器）：效率 0.884、份额 0.836，日回收 4815.0 kWh，投资 620.0 万元，年节约 110.7 万元，回收期 5.60 年，节电率 13.42% → 推荐方案 B。
  • 电价灵敏度：0.45→0.85 元/kWh 时 B 由 7.84 年缩至 4.15 年、A 由 8.78 年缩至 4.65 年，排序不反转。
• 问题三 · 动态牵引负荷：全天 160 列（16 编组 65 / 8 编组 95；上行 86 / 下行 74；CRH380B 48 / CRH3C 50 / CRH2A 62）；阈值 1.5 MW 提取 160 个脉冲、成功率 100.00%；编组辨识（阈值 6.30 MW）100.00%、方向辨识（制动比阈值 0.380）100.00%、六类最近质心 99.38%。
  • 峰值回归：P = 4.2123 + 3.8609·G16 + 0.8560·UP + 4.2130·(Kt−1) + 3.9834·G16(Kt−1)，R² = 0.9969、RMSE 0.1179 MW、MAPE 1.69%；次日预测 RMSE 0.1682 MW、MAPE 2.46%、逐小时电量 MAPE 0.75%。
  • 全天最大功率 10.336 MW、平均 1.224 MW、负荷率 11.84%。
• 问题四 · 等值建模：戴维南等值电势 27.4965 kV、内阻抗 1.3132 Ω（177 样本，R² = 0.9833，残差 RMSE 20.03 V）；24.6 km 供电臂等效 R = 4.5067 Ω、L = 33.308 mH、C = 2.1218 μF（接触网 0.2718 + 补偿 1.85）；串联谐振 598.7 Hz = 11.97 次，紧邻 11/13 次特征谐波。
  • 误差分解（保留 1/3/5/7/11/13 次谐波源）：空载 4.11%（谐波截断 4.11% + 压降 0.03%）、牵引 2.13%（1.57% + 1.44%）、制动 2.29%（2.23% + 0.52%）；平均 2.84%、最大 4.11%，低于 5% 工程阈值。
  • 补偿电容灵敏度：1.0→3.5 μF 时谐振次数由 15.46 次降至 8.98 次，途中穿越 13 次与 11 次危险点。

八、配套资源

• 回到 📚 历年真题解法库首页
• 总入口：🗺️ 专题手册库
• 论文怎么写：📝 论文模板使用说明专区
• 避坑总览：⚠️ 建模避坑指南专区
• 真源模型：🔧 tools/gen_dgcup2020a.py · 配图脚本：🖼️ tools/fig_dgcup2020a.py