1. 项目背景与核心需求
直流电气铁路牵引供电系统（TPSS）作为现代轨道交通的核心基础设施，其电能质量直接影响列车运行的稳定性和电网安全性。在实际运行中，由于电力电子器件的开关动作和负载突变，系统会产生大量特征谐波（如5次、7次、11次等）。这些谐波不仅会导致变压器过热、继电保护误动作，还会对通信系统造成电磁干扰。
单调谐无源滤波器（Passive Harmonic Filter）因其结构简单、成本低廉、可靠性高等特点，成为治理牵引供电系统谐波的主流方案。其核心原理是利用LC电路的谐振特性，在特定频率点形成低阻抗通路，从而分流该次谐波电流。本项目通过Simulink仿真平台，完整演示从参数计算到性能验证的滤波器设计全流程。
关键提示：牵引供电系统谐波治理需同时考虑滤波效果与系统兼容性，滤波器参数设计必须满足GB/T 14549-93《电能质量 公用电网谐波》的限值要求。
2. 滤波器设计原理与参数计算
2.1 单调谐电路工作原理
单调谐滤波器的核心是由电感L、电容C和电阻R组成的串联谐振电路。其阻抗特性可用公式表示：
$$
Z = R + j(\omega L - \frac{1}{\omega C})
$$
当工作频率等于谐振频率（$f_n = \frac{1}{2\pi\sqrt{LC}}$）时，电路阻抗达到最小值（约等于电阻R），此时该频率的谐波电流将主要通过滤波器回路分流。对于直流牵引系统，典型需要滤除的谐波次数为6k±1次（如5、7、11、13次等），这是由于12脉波整流器的工作特性决定的。
2.2 关键参数设计步骤
以某实际铁路牵引站为例，设计针对5次谐波（250Hz）的单调谐滤波器：
-
基础参数确定：
- 系统标称电压：27.5kV
- 基准容量：100MVA
- 需滤除谐波次数：5次（250Hz）
- 目标滤波容量：8Mvar
-
电容值计算：
根据无功补偿需求计算电容基波容抗：
$$
Q_c = \frac{V^2}{X_c} \Rightarrow X_c = \frac{27.5^2}{8} = 94.53\Omega
$$
电容值：
$$
C = \frac{1}{2\pi f X_c} = \frac{1}{314 \times 94.53} = 33.7\mu F
$$ -
电感值计算：
谐振条件：
$$
L = \frac{1}{(2\pi f_n)^2 C} = \frac{1}{(2\pi \times 250)^2 \times 33.7 \times 10^{-6}} = 12.06mH
$$ -
电阻值选取：
品质因数Q一般取30-100，此处选Q=50：
$$
R = \frac{X_L}{Q} = \frac{2\pi \times 250 \times 12.06 \times 10^{-3}}{50} = 0.379\Omega
$$
设计要点：实际电容需考虑谐波电压耐受能力，通常选用额定电压为1.5倍系统电压的电容组；电感建议采用空心电抗器以避免磁饱和。
3. Simulink建模与仿真实现
3.1 系统整体架构搭建
在Simulink中构建包含以下模块的完整系统：
- 电源模块：采用Three-Phase Programmable Voltage Source模拟27.5kV电网
- 负载模块：用Controlled Current Source模拟12脉波整流器负载
- 滤波器模块：通过Series RLC Branch实现单调谐电路
- 测量模块：包含电压/电流传感器、FFT分析仪等
关键建模技巧：
- 整流器负载采用基于开关函数的等效模型，通过注入6k±1次谐波电流模拟实际工况
- 设置Solver为ode23tb（适用于电力电子系统仿真）
- 采样时间设为1e-6s以保证谐波分析精度
3.2 参数化建模实现
通过MATLAB脚本实现参数自动配置：
matlab复制% 滤波器参数
fn = 250; % 谐振频率(Hz)
Q = 50; % 品质因数
Vn = 27.5e3; % 额定电压(V)
Qc = 8e6; % 滤波容量(var)
% 计算元件参数
Xc = Vn^2/Qc;
C = 1/(2*pi*50*Xc); % 基波频率50Hz
L = 1/((2*pi*fn)^2*C);
R = (2*pi*fn*L)/Q;
% 设置Simulink模块参数
set_param('TPSS_Filter/RLC', 'Resistance', num2str(R));
set_param('TPSS_Filter/RLC', 'Inductance', num2str(L));
set_param('TPSS_Filter/RLC', 'Capacitance', num2str(C));
3.3 典型仿真结果分析
对比接入滤波器前后的系统谐波畸变率（THD）：
| 谐波次数 | 未滤波电流(A) | 滤波后电流(A) | 衰减率(%) |
|---|---|---|---|
| 5次 | 85.6 | 9.2 | 89.3 |
| 7次 | 63.4 | 41.7 | 34.2 |
| 11次 | 27.8 | 25.1 | 9.7 |
系统总THD从12.7%降至4.3%，满足国标≤5%的要求。值得注意的是，单调谐滤波器对非谐振频率谐波（如7次、11次）也有一定衰减作用，这是通过改变系统阻抗分布特性实现的。
4. 工程实施要点与故障排查
4.1 现场安装注意事项
-
电容组配置：
- 采用H型结构多电容串联/并联组合
- 配置熔断器保护（每相至少2组冗余）
- 安装放电电阻（5分钟内降至50V以下）
-
电抗器布置：
- 空心电抗器相间距离≥1.5倍直径
- 垂直安装时需考虑地震载荷
- 避免与钢结构形成闭合磁路
-
系统调试流程：
mermaid复制graph TD A[绝缘测试] --> B[空载投切试验] B --> C[谐波频谱测量] C --> D[参数微调] D --> E[72小时连续运行测试]
4.2 常见故障处理指南
| 故障现象 | 可能原因 | 排查方法 |
|---|---|---|
| 滤波器过电流 | 电容击穿 | 测量各电容单元电压分布 |
| 谐振频率偏移 | 电感值漂移（温度影响） | 采用温度补偿型电抗器 |
| 系统THD不降反升 | 滤波器与系统阻抗谐振 | 修改滤波器分组方案 |
| 保护频繁动作 | 谐波放大效应 | 加装阻尼电阻或改为C型滤波器 |
5. 设计优化与进阶方案
5.1 参数灵敏度分析
通过蒙特卡洛仿真评估元件容差对滤波效果的影响：
matlab复制% 元件参数±5%随机变化
C_var = C*(0.95 + 0.1*rand(100,1));
L_var = L*(0.95 + 0.1*rand(100,1));
for i = 1:100
set_param('model/RLC','Capacitance',num2str(C_var(i)));
set_param('model/RLC','Inductance',num2str(L_var(i)));
simout = sim('model');
THD(i) = simout.THD.Data(end);
end
histogram(THD) % 显示THD分布统计
分析表明，电容变化对谐振频率影响更显著，建议选用±1%精度电容。
5.2 混合滤波方案
对于谐波频谱复杂的情况，可采用"单调谐+高通"的混合滤波器结构：
- 单调谐支路：针对主导谐波（如5次）
- 高通支路：滤除13次以上谐波
- 典型参数比：70%容量分配单调谐支路
这种方案在沪宁城际铁路改造项目中使THD进一步降至3.1%，同时降低成本约15%。