直流牵引供电系统数学模型与短路计算研究

直流牵引供电系统数学模型与短路计算研究摘要直流牵引供电系统在城市轨道交通等领域应用广泛，其数学模型的建立与短路计算对于系统设计、设备选型及安全运行意义重大。本文深入剖析直流牵引供电系统各组件，如整流机组、接触轨、走行轨等，构建精确数学模型。同时，对短路计算方法进行研究，涵盖稳态与暂态短路电流计算，并借助实例与仿真验证模型及计算方法的准确性，为直流牵引供电系统的优化设计与可靠运行提供有力支撑。关键词直流牵引供电系统；数学模型；短路计算；城市轨道交通一、引言随着城市规模的扩张与交通需求的增长，城市轨道交通凭借其高效、环保、大运量等优势，成为缓解城市交通拥堵的关键手段。直流牵引供电系统作为城市轨道交通的核心供电方式，为列车运行提供稳定电能。准确掌握其数学模型与短路计算方法，是保障系统安全、可靠、经济运行的基础。短路故障是直流牵引供电系统运行中可能面临的严重问题，短路电流的大小和特性直接影响设备选型、保护配置以及系统的整体稳定性。通过建立精确的数学模型并开展短路计算研究，能够为系统设计提供关键依据，提升系统运行的可靠性与安全性。二、直流牵引供电系统组成及工作原理2.1系统组成直流牵引供电系统主要由牵引变电所、接触网（接触轨）、走行轨和回流线等部分构成。牵引变电所将中压交流电转换为直流电，通过接触网（接触轨）为列车供电，列车电流经走行轨和回流线返回牵引变电所。牵引变电所内包含整流机组、变压器、开关设备、保护装置等设备。整流机组是核心设备，将三相交流电转换为直流电，常见的有6脉波、12脉波、24脉波等整流形式，脉波数越高，输出直流电压的谐波含量越低。变压器实现电压变换，满足整流机组的输入电压要求。开关设备用于控制电路的通断，保护装置则对系统进行短路、过载等故障保护。接触网（接触轨）是向列车输送电能的导电装置，分为架空式接触网和接触轨式两种。架空式接触网通过绝缘子悬挂在轨道上方，接触轨则安装在轨道一侧或中间。走行轨不仅作为列车的运行轨道，还承担着回流的功能，将列车电流返回牵引变电所。回流线则进一步降低走行轨上的电位差，减少杂散电流对周边金属结构的腐蚀。2.2工作原理在交流侧，外部电网的三相交流电经变压器降压后送入整流机组。以三相桥式6脉波整流机组为例，其工作原理基于二极管的单向导电性。在一个周期内，通过不同二极管的导通组合，将三相交流电转换为脉动直流电。其输出电压波形在一个周期内有6个脉波，表达式为：U_d=1.35U_2\cos\alpha其中，U_d为整流输出直流电压，U_2为变压器二次侧线电压有效值，\alpha为触发延迟角。对于12脉波整流机组，通常由两个6脉波整流器通过移相变压器组合而成，输出电压波形在一个周期内有12个脉波，可有效降低谐波含量。24脉波整流机组则是在12脉波整流机组基础上进一步组合，谐波抑制效果更佳。在直流侧，整流机组输出的直流电通过接触网（接触轨）供给列车。列车根据运行需求，通过受电弓（集电靴）从接触网（接触轨）获取电能，驱动电动机运行。同时，列车运行产生的电流经走行轨和回流线返回牵引变电所，形成完整的供电回路。三、直流牵引供电系统数学模型建立3.1整流机组模型3.1.124脉波整流机组工作原理及等效电路24脉波整流机组由四个6脉波整流桥组成，通过移相变压器实现不同整流桥之间的相位差。移相变压器的一次侧绕组采用三角形接法，二次侧绕组采用延边三角形接法，通过合理设计绕组匝数比和移相角度，可使四个6脉波整流桥输出电压的相位依次相差15°。其等效电路可视为由多个理想电压源、电阻和电感组成。其中，理想电压源模拟整流桥输出的直流电压，电阻和电感分别考虑整流桥内部电阻以及变压器漏感等因素的影响。在实际建模中，为简化计算，可将移相变压器和整流桥组合视为一个整体，用戴维南等效电路表示，即一个理想电压源E和一个内阻R串联。3.1.2多段线性化外特性模型整流机组的输出特性并非完全线性，为更准确描述其工作特性，采用多段线性化外特性模型。通过对不同工况下整流机组输出电压和电流的测量与分析，将其外特性曲线划分为多个线性段。例如，在低电流段，输出电压基本保持恒定，可表示为U=U_{d0}（U_{d0}为空载输出电压）；随着电流增大，输出电压因内阻压降等因素逐渐下降，可采用线性方程U=U_{d0}-R_{eq}I描述，其中R_{eq}为等效内阻，I为输出电流。通过确定各线性段的分界点及相应的线性方程系数，建立起多段线性化外特性模型，能更精确地模拟整流机组在不同工况下的输出特性。3.2接触轨与走行轨模型3.2.1直流电阻计算接触轨和走行轨的直流电阻是影响供电系统性能的重要参数。接触轨的直流电阻R_{c}可根据其材质、截面积S_{c}和长度L，利用公式R_{c}=\rho_{c}\frac{L}{S_{c}}计算，其中\rho_{c}为接触轨材料的电阻率。走行轨由于存在集肤效应，其直流电阻计算较为复杂。在低频情况下，可近似采用直流电阻公式R_{r}=\rho_{r}\frac{L}{S_{r}}，\rho_{r}为走行轨材料电阻率，S_{r}为走行轨截面积。但在高频或短路暂态过程中，集肤效应显著，电阻增大，此时可采用经验公式或通过有限元软件（如Ansoft）仿真计算修正电阻值。3.2.2电感计算钢轨内电感L_{in}主要与钢轨的几何尺寸、材质以及电流分布有关。对于单根钢轨，其内部电感可通过公式L_{in}=\frac{\mu_0}{8\pi}(1+\frac{4}{\pi}\ln\frac{D}{r})估算，其中\mu_0为真空磁导率，D为钢轨等效直径，r为钢轨半径。供电回路外电感L_{out}则与接触轨和走行轨之间的距离、相对位置以及线路长度等因素相关。可利用电磁感应原理，通过建立磁场模型计算外电感，或采用经验公式L_{out}=\frac{\mu_0}{2\pi}\ln\frac{d}{r_{eq}}，d为接触轨与走行轨中心线距离，r_{eq}为等效半径。综合考虑钢轨内电感和供电回路外电感，得到走行轨和接触轨的电感参数L=L_{in}+L_{out}。将计算得到的电感参数与实测数据对比验证。通过在实际线路上布置测量点，采用高精度电感测量仪器测量电感值，对比计算值与测量值，若偏差在合理范围内，则验证了计算方法的准确性；若偏差较大，则分析原因，如是否忽略了某些影响因素，对计算模型进行修正。3.3系统整体数学模型将整流机组、接触轨、走行轨等各部分模型整合，建立直流牵引供电系统的整体数学模型。采用电路分析方法，以基尔霍夫电压定律（KVL）和基尔霍夫电流定律（KCL）为基础，建立系统的电路方程。设系统中各支路电流为I_1,I_2,\cdots,I_n，各元件电压降为U_1,U_2,\cdots,U_n，根据KVL，对于任意闭合回路有\sum_{k=1}^{m}U_{k}=0（m为回路中元件个数）；根据KCL，对于任意节点有\sum_{i=1}^{n}I_{i}=0。结合各元件的数学模型，如整流机组的外特性模型、接触轨和走行轨的电阻电感模型等，构建联立方程组，求解方程组可得到系统在不同工况下的电流、电压分布。例如，在正常运行工况下，考虑列车负载电流，通过求解方程组可得到接触轨电压、走行轨电位差等参数；在短路故障工况下，改变短路点位置和故障类型，重新求解方程组，得到短路电流大小及分布情况。四、直流牵引供电系统短路计算4.1短路类型分析直流牵引供电系统常见短路类型包括接触轨对地短路、接触轨相间短路（双轨制供电时）、接触轨与走行轨短路等。接触轨对地短路通常由于接触轨绝缘损坏、外物搭接等原因引起。短路发生后，电流经接触轨通过短路点流向大地，再通过走行轨和回流线返回牵引变电所，可能导致接触轨电压骤降，影响列车正常运行，同时可能对周边人员和设备造成安全隐患。接触轨相间短路在双轨制供电系统中可能发生，如因异物跨越两根接触轨、电气设备故障等导致。相间短路时，短路电流较大，会对供电设备和线路造成严重冲击。接触轨与走行轨短路较为常见，可能由于列车受电弓故障、接触轨与走行轨间距过小等原因引发。短路电流通过走行轨回流，会使走行轨电流增大，可能导致走行轨发热、电位升高，影响轨道结构和周边金属设施。4.2短路电流稳态值计算4.2.1计算方法采用基于系统数学模型的解析计算方法。在短路发生稳态后，将系统数学模型中的暂态项忽略，简化为直流电阻电路模型。根据不同短路类型，确定短路回路，利用欧姆定律和基尔霍夫定律计算短路电流稳态值。例如，对于接触轨对地短路，设短路点电阻为R_{f}，根据系统数学模型，可得到短路回路的总电阻R_{total}（包括整流机组内阻、接触轨电阻、走行轨电阻以及短路点电阻等），短路电流稳态值I_{sc}=\frac{E}{R_{total}}，其中E为整流机组等效电源电动势。4.2.2计算程序与界面开发编写短路电流稳态值计算程序，采用高级编程语言（如Python）结合数值计算库（如NumPy）实现。程序流程如下：首先输入系统参数，包括整流机组参数、接触轨和走行轨参数、短路点位置等；然后根据短路类型确定短路回路，调用相应的计算模型计算总电阻；最后根据公式计算短路电流稳态值。为方便用户操作，开发图形用户界面（GUI）。利用Python的GUI框架（如Tkinter）设计界面，界面上设置输入框用于输入系统参数，按钮用于触发计算，文本框用于显示计算结果。用户在界面上输入相关参数后，点击计算按钮，程序即可快速计算并显示短路电流稳态值，提高计算效率和易用性。4.3短路电流暂态过程分析4.3.1暂态过程特点短路发生瞬间，系统进入暂态过程，电流电压迅速变化。由于系统中存在电感，电流不能突变，短路电流从初始值逐渐上升至稳态值。在暂态过程中，整流机组的输出特性、接触轨和走行轨的电感效应等因素相互作用，使得短路电流波形复杂。例如，在短路初始阶段，由于电感的阻碍作用，短路电流上升较慢，随着时间推移，电感储能逐渐释放，短路电流快速上升。同时，整流机组的外特性也会随着短路电流变化而改变，进一步影响暂态电流的波形和幅值。4.3.2仿真模型建立利用SIMULINK中电力系统仿真模块建立包含24脉波整流机组在内的直流牵引供电系统的仿真模型。在模型中，根据实际系统参数设置整流机组、接触轨、走行轨、列车负载等模块参数。整流机组模块采用前面建立的多段线性化外特性模型，通过Simulink中的模块搭建实现。接触轨和走行轨模块则根据其电阻电感模型，利用RLC串联电路模块模拟。列车负载模块根据列车运行特性，设置不同工况下的电流需求。在模型中设置不同的短路点位置和短路类型，如在接触轨不同位置设置对地短路、相间短路等故障，通过仿真运行，观察短路发生后电流的暂态变化过程，获取短路电流的波形、幅值、上升时间等参数，深入分析短路暂态特性。五、案例分析与验证5.1实际地铁线路案例选取某实际地铁线路作为案例，该线路采用直流1500V供电系统，配备24脉波整流机组。收集该线路的系统参数，包括整流机组的额定容量、变比、内阻，接触轨和走行轨的材质、规格、长度，以及列车的运行参数等。根据所建立的数学模型和短路计算方法，计算该线路在不同短路情况下的短路电流稳态值和暂态过程。如计算接触轨在距离牵引变电所不同位置发生对地短路时的短路电流稳态值，以及短路发生后0-1s内的暂态电流变化曲线。5.2与试验数据对比在该地铁线路进行短路试验，在不同位置设置短路点，测量短路电流的实际值。将计算得到的短路电流稳态值和暂态过程数据与试验测量数据进行对比。对于短路电流稳态值，对比计算值与测量值的误差。若误差在5%以内，认为计算方法准确可靠；若误差较大，分析原因，如是否参数取值不准确、计算模型简化过度等，对模型和计算方法进行修正。对于短路电流暂态过程，对比计算得到的电流波形与实测波形的相似度。通过计算波形相关系数等指标，评估仿真模型对暂态过程模拟的准确性。若相似度较高，验证了仿真模型的有效性；若差异较大，进一步分析模型中各元件参数设置是否合理，对仿真模型进行优化。六、结论本文深入研究了直流牵引供电系统的数学模型与短路计算方法。通过对系统各组成部分的详细分析，建立了精确的数学模型，包括整流机组的多段线性化外特性模型、接触轨和走行轨的电阻电感模型等，并整合形