直流牵引供电谐波治理:单调谐滤波器Simulink仿真与工程实践
2026/9/8 8:18:48 网站建设 项目流程

做电气铁路供电设计的人,基本都经历过这种焦虑:牵引变电所投运后,后台频报谐波超标,变压器响声异常,滤波电容器鼓包甚至炸裂,排查到最后发现根子出在方案设计阶段。直流牵引供电系统看起来比交流系统简单,没有牵引网无功损耗、没有电分相,但整流机组产生的谐波污染一点不比交流系统省心。我前段时间帮某地铁线路的牵引所做谐波治理,动手翻设备铭牌之前,先打开Simulink把单调谐无源滤波器方案完整仿真了一遍。这篇博文就是整个过程的记录:直流电气铁路TPSS下谐波从哪里来,单调谐滤波器怎么选参数、怎么建模、怎么分析,以及最后落地时踩过的那些坑。

无论你是正在做轨道交通供电方向课程设计的学生,还是刚接触牵引变电所设计、需要用仿真验证滤波器方案的工程师,这篇文章都能直接拿来当参考。文里所有参数、步骤、计算结果都是按工程习惯来的,你照着改数据就能用到自己的项目里。

1. 直流牵引供电系统的谐波源头:整流机组与列车变流器到底贡献了什么

1.1 交流侧特征谐波:6脉波和12脉波整流的那笔账

直流电气铁路的牵引变电所,核心功能就是把电网的三相交流电降压、整流成直流,送到接触网或第三轨上。城市轨道交通常用的直流等级是750V、1500V,老一些的矿山铁路还有3000V。整流机组通常由整流变压器加二极管或晶闸管整流桥组成。

以最基础的三相桥式6脉波整流器为例,交流侧电流波形是典型的阶梯波。对这种波形做傅里叶分解,谐波次数满足这个规律:

h = 6k ± 1,k = 1, 2, 3...

也就是说,5次、7次、11次、13次、17次、19次……这些“特征谐波”会源源不断地从整流器注入到牵引变电所的交流母线,然后顺着馈线反灌到上级电网。谐波电流的有效值大致等于基波电流除以谐波次数,I_h ≈ I_1 / h,所以5次谐波是幅值最大的那个,设计滤波器时一般优先照顾它。

如果用12脉波整流,也就是用两台整流变压器,一组Yy接法、一组Yd接法,加上两组整流桥串联,则5次和7次谐波会在变压器网侧互相抵消,剩下的特征谐波变成:

h = 12k ± 1,k = 1, 2, 3...

也就是11次、13次、23次、25次……谐波次数越高,电流幅值越小,处理起来更容易。这也是现代牵引所普遍采用12脉波甚至24脉波整流的原因之一。

但注意一个坑:12脉波整流对变压器参数的一致性非常敏感,两组整流桥的负荷分配不完全均衡时,5次、7次谐波并不会完全消失,只是幅值大大降低。所以工程上做谐波评估时,不能完全按理想12脉波去估算,最好留出余量。

1.2 直流侧纹波与列车变流器的传导干扰

谐波不只是交流侧有。整流输出的直流电压本身含有纹波,6脉波整流的直流侧特征谐波次数是6k次,也就是300Hz、600Hz及其倍数;12脉波整流则是12k次,最低是600Hz。这些纹波电流会沿着接触网传输,通过列车取流口进入车载变流器。

现代列车基本都采用PWM牵引逆变器驱动交流异步电机或永磁电机,逆变器开关频率通常在几百赫兹到几千赫兹,产生的谐波频带非常宽。还有再生制动回馈装置、辅助电源系统(蓄电池充电机、空调、照明),每一个都是宽带谐波源。这些高频分量虽然单次幅值不大,但频带宽、衰减难,很容易耦合到轨道电路和信号系统,造成轨道区段占用检测误动作。

我在现场遇到过一种典型情况:站台门联动、信号系统误报,查来查去,最后定位到一台老式充电机的高次谐波干扰。所以谐波治理不能只盯着整流机的特征谐波,还要看整个直流系统的传导骚扰。当然,针对特征谐波做单调谐滤波仍然是最基础、最必要的一步。

1.3 谐波造成的实际破坏与治理思路取舍

谐波的真实危害往往比教科书里写的更直接。谐波电流流过变压器绕组,产生附加铜耗;涡流损耗随频率平方递增,所以高频谐波会让变压器局部过热,绝缘加速老化。电容器对谐波更敏感,容抗随频率降低,高频谐波电流大部分涌进电容,导致电容器过流、鼓包、甚至炸裂。继电保护也可能被谐波影响——有些微机保护的采样计算采用傅里叶算法,谐波含量高时测量阻抗发生偏移,保护范围变化,整定好的动作边界不可靠了。

治理谐波的主流思路有两条线。一是主动治理,用APF有源滤波器实时产生反向补偿电流,效果好、动态快,但单台容量有限、成本高,在牵引所这种功率等级下往往只做局部补充。二是被动治理,也就是无源滤波器,利用LC回路的谐振特性为特征谐波提供低阻抗通路,把它“锁”在滤波支路里,不让它流入电网。单调谐无源滤波器结构最简单、损耗小、容量可以做得很大,是牵引变电所谐波治理的经典方案。

当然,无源滤波器的缺点也很明显:只能针对固定次数谐波、可能和系统阻抗产生并联谐振、滤波效果依赖元件参数稳定性。正因为这些缺点,才需要在设计阶段做充分的仿真验证,把风险提前暴露出来。

2. 单调谐滤波器参数计算:从谐振原理到一组可复现的数值

2.1 谐振极点的物理意义:为什么LC串联就能“吸”谐波

单调谐滤波器的电路结构特别简单,一个电抗器L和一个电容器C串联,再并联到母线上。它的阻抗表达式写出来是这样的:

Z = R + j(ωL - 1/ωC)

当频率等于谐振频率f_n时,感抗和容抗相等,j项里ωL = 1/ωC,阻抗只剩下电阻R这一项。也就是说,在这个频率点上,滤波器支路呈纯电阻特性,且阻抗值很小。母线上该次谐波电流遇到这条低阻抗支路,自然就往里钻,而不是流向电网。

这个原理可以用一个生活场景来理解:你推一个秋千,每次都在秋千荡到最低点时顺着它的运动方向加一点力,秋千越荡越高。LC串联回路也是这个道理,在谐振频率附近,感抗和容抗的能量交换互相抵消,外部只需要很小的电压就能建立很大的电流,表现出来就是低阻抗。

但要记住,低阻抗是“窄带”的。偏离谐振频率,滤波器阻抗迅速增大,滤波效果急剧下降。这既是单调谐滤波器的特点,也是它的软肋——后面讲仿真和失谐分析时,这个特性会反复出现。

2.2 参数计算的完整流程:先定无功,再算电容电感

设计单调谐滤波器,核心任务就是确定C、L、R三个参数。具体步骤我按工程习惯整理如下。

第一步,确定目标谐波次数。根据整流类型来定,6脉波整流选5次,12脉波整流选11次或13次,优先选幅值大的那一次。

第二步,确定滤波器的基波无功补偿容量Q_f。单调谐滤波器在基波下呈现容性,会向系统注入无功。容量不能随便定,需要考虑母线无功需求、电压支撑要求、谐波电流大小三个因素。一般按牵引变压器容量的15%到30%来估算,具体值要通过潮流计算或者仿真校核。

第三步,计算电容值。三相星形接法下:

C = Q_f / (3 × ω_1 × U_ph²)

其中U_ph是母线相电压有效值,ω_1是基波角频率。这里用相电压是因为每相电容承受的是相电压。

第四步,计算电感值。根据谐振条件nω_1L = 1/(nω_1C),得到:

L = 1 / (n² × ω_1² × C)

第五步,确定电阻值。R决定滤波器通带宽度和品质因数,工程上定义品质因数Q值:

Q = nω_1L / R

调谐滤波器的Q值一般取30到80之间,Q值越高,滤波越尖锐,但对元件参数的漂移越敏感。实际电抗器的电阻不可能做到很小,通常会把电抗器的品质因数直接当成滤波器支路的品质因数来用。

第六步,校核电容的过电压和过电流能力。电容承受的电压是基波分量和各次谐波分量的叠加,谐波电流大时电容端电压可能显著升高,选型时额定电压至少留出10%到20%的裕度。

2.3 一个完整的算例:10kV母线2Mvar的5次单调谐滤波器

我拿一个典型牵引所的参数来做具体计算。假设交流母线线电压U_L = 10kV,整流机组容量12.5MVA,需要补偿无功Q_f = 2Mvar,目标治理5次谐波,系统频率50Hz。

先算相电压:

U_ph = 10000 / √3 ≈ 5774V

再算电容:

C = 2000000 / (3 × 314 × 5774²) = 2000000 / (942 × 33.34 × 10⁶) ≈ 2000000 / 31.4 × 10⁹ ≈ 63.7 × 10⁻⁶ F ≈ 63.7μF

这个值落在标准电容器容量的常见范围内。接下来算电感:

L = 1 / (25 × 314² × 63.7 × 10⁻⁶) = 1 / (25 × 98596 × 63.7 × 10⁻⁶) ≈ 1 / 157.05 ≈ 6.37mH

品质因数取Q = 50,则电阻:

R = 5 × 314 × 6.37 × 10⁻³ / 50 ≈ 0.2Ω

好,这样一组参数就出来了。为了核对电容器是否在安全范围,我再算一下基波电容电流I_C1。基波下每相容抗X_C1 = 1/(ω₁C) = 1/(314 × 63.7 × 10⁻⁶) ≈ 50Ω,于是I_C1 = U_ph / X_C1 ≈ 5774/50 ≈ 115.5A。三相总的无功为√3 × 10kV × 115.5A ≈ 2Mvar,与预设补偿容量一致。5次谐波电流如果按100A估算,电容上5次谐波电压为100/(5×314×63.7×10⁻⁶) ≈ 1000V,叠加在基波电压上,电容峰值电压约√(5774² + 1000²) ≈ 5860V,比额定相电压高出约1.5%,属于可接受范围。但如果谐波电流达到200A,电容电压会进一步升高,选型时就该提高电容额定电压等级。这个校核步骤在工程上不能省。

2.4 品质因数和失谐:决定滤波器现场表现的关键点

品质因数这个参数,设计时常被忽略,但它对现场表现影响非常大。Q值高,谐振点尖锐,滤波效率高,但是电网频率稍有波动、电感温度漂移、电容老化容值下降,谐振频率就会偏离实际谐波频率,滤波效果断崖式下跌。Q值低,通带变宽,对参数漂移容忍度高,但同样的C和L下滤波效果变差,且电阻损耗增大。

工程上有个经验做法:如果电网频率波动较大或者环境温度变化明显,Q值取30到50;如果谐波频谱相对稳定,可以取60到80。另外电容器通常有负温度系数,温度升高容值下降;电抗器电感量也会随电流和温度变化。以5次滤波器为例,谐振频率f_n = 1/(2π√(LC)),如果L和C各漂移2%,谐振频率会偏移约2%,也就是10Hz左右。5次谐波是250Hz,偏移10Hz后,滤波器阻抗会明显变大,滤波效果肉眼可见地变差。这个现象在仿真里很好复现,后面我会专门做一次失谐仿真。

3. Simulink仿真建模:先把系统简化对,再谈滤波效果

3.1 模型结构与简化边界

搭建Simulink仿真模型时,最容易犯的错误是“什么都要仿真,结果什么都仿不准”。牵引供电系统是一个庞大的网络,完整建模几乎不可能。我的做法是先明确仿真要回答什么问题,再决定简化的尺度。

本次仿真要回答三个问题:第一,不加滤波器时,牵引所交流母线电流和电压的谐波畸变率是多少;第二,加上5次单调谐滤波器后,谐波降到什么水平;第三,在失谐和过负荷工况下,滤波器表现如何。

围绕这三个目标,模型做如下简化。

  • 电网侧用理想三相电压源串联短路阻抗表示。短路容量取500MVA,折算到10kV母线,阻抗约0.2Ω,X/R比取7。
  • 整流变压器用Simulink自带的三相双绕组变压器模块,容量12.5MVA,变比10kV/1.2kV,短路阻抗8%,连接组别Dyn11。
  • 整流器用Universal Bridge模块,选择Diode六脉波整流。虽然后面我们要提12脉波,但模型用6脉波更容易让读者看清5次谐波的主导地位。
  • 直流侧负载用一个电阻串联小电感表示,模拟列车牵引工况。电阻0.35Ω、电感1mH,对应直流电流约1600V/0.35Ω≈4570A,直流功率约7.3MW,对于试验性牵引所来说合理。

滤波器用三个串联RLC支路模块,分别接到三相母线和地之间,星形连接。参数就是上一节算出来的C=63.7μF、L=6.37mH、R=0.2Ω。

3.2 逐步搭建流程与关键设置

搭建顺序我建议从电源往负载方向走,避免遗漏连接点。

第一步,新建模型,从Simscape Electrical Specialized Power Systems库中拉入Powergui模块,仿真类型选择离散,采样时间1e-5秒。离散仿真在变流器系统中比连续仿真稳定得多,跑起来也快。

第二步,搭建电网和变压器部分。三相电压源模块设线电压10kV,频率50Hz,内阻抗按前面的短路容量折算。变压器模块设好容量、变比和短路阻抗后,记得把连接组别设为Dyn11,否则波形相位关系可能不对。

第三步,搭整流器和直流负载。Universal Bridge模块选择Diode,其他参数保持默认即可。直流侧接RL串联负载,加一个直流电压表和电流表,方便记录平均电压和电流。

第四步,加入滤波器支路。滤波器使用Series RLC Branch模块,连接方式选串联,填入C、L、R三个数值。这个模块默认的是并联形式,一定要仔细检查连接方式,我见过不少同事在这一步选错导致仿真结果完全不对。

第五步,加入测量点。在滤波器的接入点电压、输入电流、滤波器支路电流位置分别放置电压电流测量模块。我给输入电流命名为I_source,滤波器电流命名为I_filter,母线电压命名为U_bus,方便后续在FFT分析时选择信号。

第六步,设置仿真时长。我通常跑到0.4秒,前0.2秒让系统进入稳态,取后0.2秒的波形做FFT分析,避开启动暂态。

3.3 测量与FFT分析配置:别让频谱泄漏骗了你

Simulink自带的Scope只能看波形,算不了THD。正确的做法是Powergui模块里的FFT Analysis工具,或者用To Workspace把波形数据导到MATLAB工作区,自己写代码分析。

用Powergui的FFT工具时注意三个坑。

第一,信号必须是稳态后的一个完整周期倍数。取0.02秒的整数倍,比如取0.1到0.3秒这0.2秒,正好10个周期。窗口长度不是工频周期的整数倍时,频谱会泄漏,谐波幅值看起来比实际小,分布却更散。

第二,采样步长必须能覆盖要分析的谐波。如果要分析到31次谐波(1550Hz),采样率至少是它的10倍,也就是16kHz以上,对应采样时间6.25e-5秒。我用的1e-5秒采样时间,完全没有问题。

第三,看三相电流时,建议只分析A相或者先做三相平衡检验。整流系统在三相完全对称时三相波形是120度互差关系,任何一相的数据都代表整个系统。但实际模型中由于数值计算的微小不对称,三相波形可能有细微区别,取三相平均值更稳妥。

3.4 仿真参数设置与常见报错处理

离散仿真的采样时间不要拍脑袋设。我测试过1e-4秒和1e-5秒对THD结果的影响,1e-4秒时整流器换相过程失真明显,THD偏高;1e-5秒时结果收敛。再加密到1e-6秒,结果几乎不变,但仿真时间成倍增加。所以1e-5秒是精度和速度的平衡点。

仿真跑到一半不收敛怎么办?最常见的原因是变压器励磁涌流导致初始阶段数值震荡。处理办法有两种:一种是把仿真返回到初值设置为“Start from steady state”,让Powergui先自动计算稳态初值;另一种是在变压器参数里把励磁电阻和励磁电感加到合理值,别用0。还有一种情况是整流器二极管默认带Snubber电路,在某些条件下引起数值刚性,可以把Snubber电阻改成1e5Ω以上或者直接设成inf,Snubber电容改成0,变流器模块照样正常工作。

4. 仿真结果对比与失谐工况分析:数据会告诉你该怎么调

4.1 滤波前后频谱对比:从波形看不出东西,FFT才说实话

仿真跑完先看整体波形,再进FFT分析。滤波前的模拟结果显示,交流母线输入电流波形明显不是正弦波,顶部比较平坦,底部有阶梯状突变,这是典型的六脉波整流特征。

用Powergui FFT工具看频谱,基波电流约460A,5次谐波电流约92A,相对基波约20%;7次谐波约65A,约14%;11次、13次也有但幅值小了很多。总谐波畸变率THD大约27%左右,远超国标对公用电网5%的要求。即便牵引负荷接入的是专用供电系统,这么大的畸变率也过不了验收。

投入5次单调谐滤波器后,220Hz到270Hz这个频段的电流分布发生显著变化。5次谐波电流从92A降到了约18A,降幅约80%。7次谐波几乎不受影响,因为滤波器在350Hz处阻抗远大于5次处,谐波不会往滤波器里跑。总THD从27%降到约10%。如果现场要求把THD降到5%以下,就需要再增加一组7次滤波器,或者把滤波支路设计成5次+7次组合。这个结论也验证了单调谐滤波器“一对一治理”的特点。

4.2 无功补偿与电压支撑效果:滤波器的第二贡献

滤波器投入后,母线电压有效值也发生了变化。滤波前,因为无功不足,母线电压大约是9.6kV左右,低于10kV额定值。滤波器投入后,它向系统注入2Mvar容性无功,母线电压升高到9.85kV,改善了约2.5%。

功率因数从0.88左右提升到了0.97以上。这正好说明单调谐滤波器在工程上常被当作“谐波治理+无功补偿”一机两用的装置。设计时如果只看THD不看功率因数,可能浪费了这部分无功能力。

但这里要提一个设计时要警惕的现象:如果滤波器容量过大,无功过剩会造成母线电压抬升过高。尤其在牵引所轻载工况,负载消耗的无功小,滤波器的无功注入占主导,电压可能超过允许上限。所以滤波器容量不能单看谐波治理需求,还要校核轻载时的电压水平。这也是为什么不少牵引所给滤波器配了自动投切开关,负荷小了就切除一组。

4.3 失谐和过载模拟:把最坏的工况在仿真里先走一遍

这部分是我觉得整个仿真过程最有价值的地方。我故意把电感值从6.37mH调到6.69mH,模拟电感因温升发生5%的正偏差。这个偏移量不算大,但谐振频率从250Hz下移到约243Hz,偏离了5次谐波频率7Hz。

滤波效果随即恶化:5次谐波电流从18A回升到约40A,THD从10%涨到18%。这就是典型的失谐问题。工程现场经常出现这种情况:滤波器投运初期效果很好,运行一两年后效果越来越差,排查发现是电抗器温升导致电感量漂移,或者电容器老化容值下降。

模拟过载工况时,我把直流侧负载电阻从0.35Ω降到0.25Ω,模拟列车加速时的重载工况。直流电流从4570A升到约6400A,交流侧基波电流同步增大。此时5次谐波电流的绝对值会增加,但由于基波也变大了,谐波含有率变化不大。真正需要注意的是滤波器支路本身的电流是否超限。在重载工况下,滤波器支路的基波电流和谐波电流叠加后的有效值达到约140A,接近了电容器的额定电流。如果谐波电流再大一点,电容器就处于过载边缘了。这个仿真结果提醒我,选电容器额定电流时不能只看基波,一定要把谐波分量算进去。

5. 从仿真到牵引变电所的工程落地经验

5.1 仿真与实际的差距:为什么现场数据和模型对不上

仿真一切正常、现场问题一堆,这种情况我见得太多了。原因主要集中在以下几点。

一是电网背景谐波。仿真里假设电网电压是纯正弦,但实际牵引所接入的电网本身就有一定谐波背景,特别是同一母线上还挂着其他非线性负荷时,背景谐波会与滤波器发生交互。这个现象在仿真中很难复现,除非你能拿到接入点实际的背景谐波频谱,我建议有条件的话做电能质量实测后再校核滤波器设计。

二是元件实际参数偏差。电容器出厂容差通常是-5%到+10%,电抗器电感量也有制造偏差。设计时的谐振点本来就是个理论值,实际投运后可能需要现场微调电抗器的抽头来对准谐波频率。这就需要在设计阶段预留调节空间。

三是温度与老化效应。电容器容值随温度和工作时间缓慢下降,电抗器电感量随温度上升略有增加,两者叠加会让谐振频率在多年运行周期里逐渐偏移。所以工程上优选品质因数不要过高的设计,用适度的通带宽度给失谐留出余地。

5.2 多组滤波器并联时的谐振风险:5次和7次会相互干扰吗

当系统需要同时治理5次和7次谐波时,两组单调谐滤波器并联运行,问题就来了。在5次谐振频率附近,7次滤波器呈现容性;在7次谐振频率附近,5次滤波器也呈现容性。两者组合后可能在某个频率点形成一个并联谐振峰,如果这个峰值恰好落在某个非特征谐波频率上,反而会放大该次谐波。

这个风险在Simulink里很好验证。我把5次和7次两组滤波器都搭进去,扫描1到20次的阻抗频率特性,发现在4.6次附近出现了一个明显的阻抗峰值。也就是说,如果系统里恰巧存在4.6次左右的间谐波(比如某些变频设备产生的),这个间谐波会被放大。处理办法有两种:一是每组滤波器串联一定比例的电抗,抬高并联谐振频率;二是采用C型滤波器或高通阻尼滤波器替代其中一组,增加阻尼消耗掉谐振能量。具体选哪种,要结合系统背景谐波频谱来定。

5.3 滤波方案选型:单调谐、C型和高压APF怎么取舍

从事无源滤波多年的经验告诉我,方案选型最忌讳“一把尺子量到底”。

单调谐滤波器的优势是结构简单、损耗小、单组容量大、运维门槛低。缺点是只能治理单次谐波、失谐敏感、有并联谐振风险。它最适合特征谐波明显、谐波频谱相对稳定的场景,比如传统牵引变电所的5次、7次谐波治理。

C型滤波器在单调谐的基础上增加了一个电容和电阻串联的阻尼支路,对失谐的耐受能力更强,高次谐波抑制效果也更好,但元件数量多、占地大、造价高。通常用在谐波频谱复杂或者电网阻抗波动较大的场景。

高压APF技术这几年进步很快,但当容量做到Mvar级别时成本依然远超无源方案,且电力电子设备自身的可靠性、电磁兼容问题在牵引所里又会产生新的课题。我的实践经验是:优先考虑无源方案做大容量治理,无源搞不定的特殊频次或动态波动大的部分,再用APF补充。

最后再分享一个落地时的小技巧:滤波器投运后不要急着验收波形。先让系统稳定运行24小时,用便携式电能质量分析仪在滤波器的进线端和母线PT处各测一轮数据,对比投运前后的THD、各次谐波电流、功率因数变化。如果发现某次谐波不降反升,八成是并联谐振或参数失谐,及时调整电抗器抽头或者重新整定品质因数,现场大部分问题都能解决。仿真做得再漂亮,都不如实测数据有信心。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询