直流牵引供电系统数学模型构建与短路计算方法的深度剖析

直流牵引供电系统数学模型构建与短路计算方法的深度剖析一、引言1.1研究背景与意义随着全球城市化进程的加速，城市人口不断增长，交通拥堵问题日益严重。城市轨道交通作为一种高效、便捷、环保的公共交通方式，在各大城市得到了广泛的应用和快速的发展。例如，截至2023年底，中国内地累计有53个城市开通城市轨道交通运营线路292条，运营里程9584公里，车站5760座。直流牵引供电系统作为城市轨道交通的关键组成部分，承担着为列车提供动力的重要任务。直流牵引供电系统具有功率大、效率高、结构简单、运行稳定等优点，能够满足城市轨道交通列车频繁启动、加速、制动等运行工况的需求。在城市轨道交通中，直流牵引供电系统将交流高压通过整流机组转换为直流电源，再通过接触网或接触轨将直流电输送给列车，为列车的牵引电机提供动力。其典型的供电电压等级有DC750V和DC1500V，不同的电压等级适用于不同的线路和车型。然而，直流牵引供电系统在运行过程中，不可避免地会面临各种电气故障问题，其中短路故障是较为常见且危害严重的一种。短路故障的发生会导致电流急剧增大，可能引发设备损坏、火灾等严重后果，对列车的运行安全和乘客的生命财产安全构成巨大威胁。例如，2018年某城市地铁线路曾发生一起直流牵引供电系统短路故障，导致列车紧急停车，部分设备烧毁，造成了数小时的运营中断，给城市交通和居民生活带来了极大的不便。短路计算在直流牵引供电系统中具有至关重要的作用，它是保障系统安全稳定运行的关键环节。精确的短路计算能够为系统的设计、设备选型、保护装置整定以及故障分析等提供重要依据。在系统设计阶段，通过短路计算可以合理确定牵引变电所的容量、布局以及供电线路的参数，确保系统能够满足列车运行的需求，同时提高系统的经济性和可靠性。在设备选型方面，短路计算结果是选择直流快速开关、熔断器等设备的重要依据，只有选择合适的设备，才能在短路故障发生时迅速切断故障电流，保护设备和系统的安全。对于直流保护装置的准确整定，短路计算更是不可或缺。通过计算不同运行工况下的短路电流，能够为保护装置设定合理的动作值和动作时间，使其在故障发生时能够快速、准确地动作，切除故障，避免事故的扩大。短路计算还可以帮助分析故障原因，为故障的排查和修复提供指导，提高系统的维护效率和运行可靠性。1.2国内外研究现状国外对直流牵引供电系统的研究起步较早，德国、日本等轨道交通发达国家在早期主要通过现场试验获取电气参数，如德国在16.7Hz单相交流供电系统建设与运营中，就对线路电感、电阻等参数进行了细致测量。在数学模型和短路计算方面，早期的研究主要集中在对系统基本原理和简单模型的构建上。随着技术的发展，数值计算方法逐渐兴起，有限元仿真等技术被广泛应用于电气参数计算，为直流牵引供电系统的研究提供了更精确的手段。国内对直流牵引供电系统的研究也取得了丰硕的成果。在系统数学模型建立方面，学者们深入分析了直流牵引变电所中整流机组的工作原理，提出了多段线性化外特性模型，并采用理想电压源、内阻的戴维南等效电路代替直流牵引变电所，从而建立起直流牵引供电系统的数学模型。在短路计算方法上，除了传统的电路图法，还结合计算机仿真技术，利用MATLAB等专业软件对系统进行建模和仿真分析，提高了短路计算的准确性和效率。然而，现有研究仍存在一些不足之处。一方面，在数学模型的精确性方面，虽然已经取得了一定进展，但由于直流牵引供电系统的复杂性，部分模型在考虑实际运行中的多种因素时还存在一定的局限性。例如，在考虑系统中电力电子器件的非线性特性、线路参数的分布性以及不同运行工况下的动态特性等方面，模型的精确性还有待进一步提高。另一方面，在短路计算方面，虽然各种计算方法和仿真技术不断涌现，但在计算的快速性和准确性之间还难以达到完美的平衡。部分计算方法在处理复杂网络和多故障情况时，计算速度较慢，无法满足实时性要求；而一些追求快速计算的方法，其计算结果的准确性又可能受到一定影响。现有研究在将理论成果应用于实际工程时，还存在一定的障碍，需要进一步加强理论与实践的结合，提高研究成果的实用性。1.3研究内容与方法本文围绕直流牵引供电系统数学模型与短路计算展开研究，主要研究内容如下：直流牵引供电系统数学模型构建：深入剖析直流牵引供电系统的电路结构与工作原理，涵盖直流牵引变电所、接触网、走行轨以及列车等关键组成部分。例如，详细分析直流牵引变电所中整流机组的工作过程，研究其将交流高压转换为直流电源的原理，并采用多段线性化外特性模型来描述整流机组的特性。通过对各组成部分的精确分析，建立起全面、准确的直流牵引供电系统数学模型，明确系统中电压、电流、功率等相关变量的算式和方程式，为后续的短路计算奠定坚实基础。直流牵引供电系统短路计算方法研究：全面统计和深入分析直流牵引供电系统的短路情况，研究短路发生时的特性和规律。在传统电路图法的基础上，结合现代计算机技术和数值计算方法，提出一种创新的、适合于直流牵引供电系统的短路计算方法。该方法充分考虑系统中各种元件的特性以及不同运行工况的影响，能够更快速、准确地计算出短路电流的大小和分布，有效解决现有计算方法在计算速度和准确性方面存在的不足。数据采集与仿真实验：通过现场实验测量，采集直流牵引供电系统在实际运行中的各种电气参数，如电压、电流、电阻、电感等。利用这些实测数据，建立系统的仿真模型，并运用MATLAB、SIMULINK等专业仿真软件进行电气仿真实验分析。将仿真结果与实际测量数据进行对比验证，评估新提出的短路计算方法的有效性和精度，确保研究成果的可靠性和实用性。计算机模拟与优化分析：基于MATLAB等专业计算机软件工具，构建直流牵引供电系统的数学模型，进行计算机模拟和优化分析。通过模拟不同运行工况和故障条件下的系统响应，深入研究系统的性能和稳定性。根据模拟结果，提出相应的改进建议和优化方案，为直流牵引供电系统的设计、运行和维护提供科学依据，提高系统的安全性和可靠性。在研究方法上，本文综合运用理论分析、案例研究和仿真验证相结合的方式：理论分析：从直流牵引供电系统的基本原理出发，运用电路理论、电磁学等相关知识，对系统的数学模型和短路计算方法进行深入的理论推导和分析，揭示系统的内在规律和本质特性。案例研究：选取实际的城市轨道交通线路中的直流牵引供电系统作为案例，对其进行详细的分析和研究。通过对实际案例的研究，获取真实的数据和运行经验，验证理论研究成果的可行性和实用性，并为进一步的研究提供实践基础。仿真验证：利用计算机仿真技术，对直流牵引供电系统进行建模和仿真分析。通过仿真实验，可以在虚拟环境中模拟各种复杂的运行工况和故障情况，快速、准确地获取系统的响应数据。将仿真结果与理论分析和实际案例进行对比验证，进一步完善和优化研究成果，提高研究的科学性和可靠性。二、直流牵引供电系统的构成与工作原理2.1系统主要组成部分直流牵引供电系统主要由牵引变电所、接触网、走行轨、回流线以及电力监控系统等部分组成，各部分相互协作，共同为列车提供稳定可靠的电能。牵引变电所是直流牵引供电系统的核心部分，其主要功能是将城市电网的三相交流电降压并整流为适合列车使用的直流电。以常见的DC1500V直流牵引供电系统为例，牵引变电所首先通过降压变压器将110kV或35kV的高压交流电降至合适的电压等级，如33kV或10kV。然后，利用整流机组将降压后的交流电转换为直流电，整流机组通常采用多脉波整流技术，如12脉波或24脉波整流，以减少谐波对电网的影响。在实际运行中，为了提高供电的可靠性和灵活性，牵引变电所一般配备两台或多台整流机组，当其中一台发生故障时，其他机组仍能继续运行，保证对列车的供电。接触网是直流牵引供电系统向列车传输电能的重要通道，它通过绝缘子悬挂在轨道上方或敷设在轨道侧面，列车通过受电弓与接触网接触，获取电能。接触网的结构形式多样，主要包括刚性接触网和柔性接触网。刚性接触网具有结构紧凑、占用空间小、维护方便等优点，常用于城市轨道交通的地下线路；柔性接触网则具有弹性好、受流性能优良等特点，适用于地面和高架线路。在接触网的设计和建设中，需要考虑诸多因素，如导线的材质、张力、悬挂高度以及接触网的分段和绝缘等，以确保其能够满足列车高速、频繁启停等运行工况的需求，为列车提供稳定、可靠的电能。走行轨不仅是列车运行的轨道，还承担着回流的重要作用，它构成了牵引供电回路的一部分。在列车运行过程中，从接触网获取的电流通过列车的牵引电机后，经走行轨返回牵引变电所。走行轨通常采用特殊的钢材制造，具有良好的导电性和耐磨性，能够满足大电流回流的要求。为了降低走行轨的电阻，减少电能损耗，在实际工程中，会采取一些措施，如增加走行轨的截面积、采用并联回流线等。同时，为了防止走行轨上的电流对地下金属设施造成腐蚀，还需要采取有效的迷流防护措施，如安装排流装置、设置绝缘层等。回流线是连接走行轨和牵引变电所的导线，其作用是将走行轨上的回流电流顺利地引回牵引变电所，形成完整的供电回路。回流线一般采用截面积较大的导线，以降低电阻，减少电能损耗。在一些情况下，为了提高回流的可靠性，还会设置多条回流线。回流线与走行轨之间通过接地装置连接，确保回流电流能够安全、稳定地返回牵引变电所。电力监控系统是直流牵引供电系统的“大脑”，它对整个系统的运行状态进行实时监测和控制。电力监控系统通过各种传感器和监测设备，采集牵引变电所、接触网、走行轨等部分的电压、电流、温度等运行参数，并将这些数据传输到监控中心。监控中心的工作人员可以根据这些数据，及时了解系统的运行情况，发现并处理故障。电力监控系统还可以实现对牵引变电所内设备的远程控制，如开关的分合闸、整流机组的调节等，提高系统的运行效率和可靠性。在现代化的城市轨道交通中，电力监控系统还与列车运行控制系统、综合监控系统等进行集成，实现信息共享和协同工作，为城市轨道交通的安全、高效运行提供有力保障。2.2工作原理阐述直流牵引供电系统的工作原理涉及多个环节，从交流高压到直流供电的转换过程以及电能传输和分配的原理如下：交流高压转换为直流电源：在牵引变电所内，首先通过降压变压器将城市电网的三相高压交流电（如110kV或35kV）降至适合整流的电压等级，如33kV或10kV。以常见的12脉波整流机组为例，其由两个6脉波整流器组成，通过变压器的特殊接线方式，使两个6脉波整流器的输出电压在相位上相差30°，从而合成12脉波的直流输出，有效减少了谐波含量。整流机组采用的整流技术主要有二极管整流和晶闸管整流。二极管整流具有结构简单、成本低、可靠性高的优点，但其输出电压不可调节；晶闸管整流则可以通过控制晶闸管的触发角来调节输出电压，灵活性较高，但电路相对复杂，成本也较高。在实际应用中，会根据系统的具体需求选择合适的整流技术。电能传输：经过整流机组转换得到的直流电，通过正馈电线传输至接触网。接触网作为向列车传输电能的通道，根据线路的不同特点和需求，分为刚性接触网和柔性接触网。刚性接触网通常由汇流排和接触线组成，汇流排采用铝合金材质，具有良好的导电性和机械强度，能够支撑和固定接触线；接触线则采用高导电率的铜合金材料，以确保电能的高效传输。柔性接触网由承力索、接触线、吊弦等部件组成，承力索主要承受接触网的机械负荷，接触线负责向列车供电，吊弦用于调节接触线的高度和张力，使接触线与受电弓保持良好的接触。列车通过受电弓与接触网接触，获取电能，受电弓的滑板采用碳滑板或铜基粉末冶金滑板，具有良好的耐磨性和导电性，能够在高速运行和频繁启停的工况下，稳定地从接触网获取电能。电能分配与回流：电能进入列车后，根据列车的运行需求，分配到牵引电机、辅助设备等各个部分。牵引电机将电能转换为机械能，驱动列车运行；辅助设备则包括空调、照明、控制系统等，为列车的正常运行和乘客的舒适提供保障。在列车运行过程中，电流经走行轨和回流线返回牵引变电所，形成完整的供电回路。走行轨不仅是列车运行的轨道，还承担着回流的重要作用，为了降低走行轨的电阻，减少电能损耗，通常会采取增加走行轨截面积、采用并联回流线等措施。回流线则将走行轨上的回流电流引回牵引变电所，确保电流能够安全、稳定地返回。2.3系统在不同场景下的应用实例直流牵引供电系统凭借其独特的优势，在多种交通场景中得到了广泛应用，下面将以城市轨道交通和铁路货运等场景为例，详细说明其特点和适用性。在城市轨道交通领域，以上海地铁为例，其广泛采用了直流牵引供电系统，电压等级主要为DC1500V。上海地铁的线路复杂，客流量巨大，高峰时段部分线路每小时的客流量可达数万人次。直流牵引供电系统能够满足列车频繁启动、加速、制动等运行工况的需求，确保列车在不同的运行条件下都能稳定运行。在高峰时段，列车需要频繁地启动和加速，以满足大量乘客的出行需求。直流牵引供电系统能够快速响应列车的功率需求，提供足够的电能，使列车能够迅速达到运行速度，减少乘客的等待时间。直流牵引供电系统的可靠性高，能够有效降低列车运行故障的发生率，保障地铁的安全、高效运营。上海地铁通过对直流牵引供电系统的精心维护和管理，确保了系统的稳定运行，使得地铁的正点率始终保持在较高水平，为市民的出行提供了可靠的保障。在铁路货运场景中，大秦铁路是一个典型的应用案例。大秦铁路是我国重要的煤炭运输专线，其运输任务繁重，年运量可达数亿吨。该铁路采用了直流牵引供电系统，电压等级为DC3000V。这种高电压等级的直流牵引供电系统能够满足铁路货运列车大功率、长距离运输的需求。大秦铁路上的货运列车通常载重量较大，需要强大的动力来牵引。直流牵引供电系统能够提供稳定的大功率电能，确保列车在重载情况下能够顺利启动和运行，实现长距离的煤炭运输。直流牵引供电系统在应对恶劣环境条件方面具有优势。大秦铁路部分路段地处山区，气候条件复杂，冬季寒冷，夏季炎热，且多风沙。直流牵引供电系统的设备具有良好的适应性，能够在恶劣的环境条件下正常运行，保证了铁路货运的连续性和稳定性，为我国的能源运输做出了重要贡献。三、直流牵引供电系统数学模型构建3.1电气参数计算在直流牵引供电系统中，电气参数的准确计算对于系统数学模型的构建和短路计算至关重要。这些参数直接影响着系统的性能和运行特性，下面将详细介绍直流电阻和电感参数的计算方法。3.1.1直流电阻计算以第三轨供电为例，走行轨和接触轨的直流电阻计算是构建系统数学模型的基础。走行轨通常采用特殊的钢材制造，其直流电阻的计算需要考虑材料的电阻率、截面积以及长度等因素。根据电阻的基本计算公式R=\rho\frac{L}{S}（其中R为电阻，\rho为电阻率，L为长度，S为截面积），对于走行轨，其电阻率\rho与钢材的材质有关，一般可通过材料手册获取。截面积S则根据走行轨的具体形状和尺寸进行计算，如常见的工字形走行轨，其截面积需要考虑腹板和翼缘的面积。接触轨的直流电阻计算方法与走行轨类似，但由于接触轨的材料和结构可能与走行轨不同，其电阻率和截面积的计算也有所差异。例如，接触轨常用的材料有钢铝复合轨，这种材料结合了钢的高强度和铝的良好导电性。对于钢铝复合轨，其直流电阻的计算需要考虑钢层和铝层的电阻率以及它们在整个截面中所占的比例。假设钢铝复合轨中钢层的电阻率为\rho_1，截面积为S_1；铝层的电阻率为\rho_2，截面积为S_2，则钢铝复合轨的等效电阻率\rho_{eq}可通过公式\rho_{eq}=\frac{\rho_1S_1+\rho_2S_2}{S_1+S_2}计算得到，再根据上述电阻计算公式即可求得接触轨的直流电阻。在实际工程中，为了降低接触轨的电阻，还会采取一些措施，如增加接触轨的截面积、优化接触轨的材料配比等。3.1.2电感参数计算在直流牵引供电系统中，电感参数的计算对于准确分析系统的电磁特性和短路电流有着关键作用，下面将探讨钢轨内电感和供电回路外电感的计算方法，并对比计算值与实测数据。钢轨内电感的计算较为复杂，它与钢轨的几何形状、材质以及电流的分布等因素密切相关。对于圆形截面的钢轨，其单位长度内电感L_{in}可通过公式L_{in}=\frac{\mu_0}{8\pi}计算，其中\mu_0为真空磁导率，其值为4\pi\times10^{-7}H/m。然而，实际的钢轨通常不是圆形截面，如常见的工字形钢轨，其内部磁场分布较为复杂，需要采用更为精确的方法进行计算。一种常用的方法是基于有限元分析，通过建立钢轨的三维模型，利用有限元软件求解磁场分布，进而计算出内电感。以某实际地铁线路的工字形钢轨为例，采用有限元软件计算得到的单位长度内电感为0.68\times10^{-6}H/m。供电回路外电感的计算主要考虑接触轨与走行轨之间的互感以及它们与周围环境的影响。假设接触轨与走行轨之间的距离为d，单位长度供电回路外电感L_{out}可近似通过公式L_{out}=\frac{\mu_0}{2\pi}\ln\frac{d}{r}计算，其中r为接触轨和走行轨的等效半径。在实际计算中，还需要考虑轨道的铺设方式、周围金属结构的影响等因素对电感的修正。为了验证计算值的准确性，将计算得到的电感参数与实测数据进行对比。通过在实际的直流牵引供电系统中布置传感器，测量不同位置的电感值。以某地铁线路为例，在一段长度为1000m的线路上进行实测，得到的钢轨内电感平均值为0.70\times10^{-6}H/m，供电回路外电感平均值为1.25\times10^{-6}H/m。与计算值相比，钢轨内电感的计算值与实测值相对误差约为2.86\%，供电回路外电感的计算值与实测值相对误差约为3.2\%。误差产生的原因主要包括计算模型的简化、实际线路中存在的一些难以精确量化的因素，如轨道连接处的接触电阻、周围环境中的杂散磁场等。通过对比可以发现，虽然计算值与实测值存在一定误差，但在工程允许的范围内，上述计算方法能够较为准确地估算电感参数，为直流牵引供电系统的数学模型构建和短路计算提供了可靠的依据。3.2整流机组模型3.2.1工作原理分析整流机组作为直流牵引供电系统的关键设备，其核心组成部分包括整流变压器和硅整流器，它们协同工作，实现了将交流电转换为直流电的重要功能。整流变压器的工作基于电磁感应原理，与普通变压器类似，它通过初、次级绕组在交变磁通的作用下，实现电压的变换。其主要功能是为整流系统提供合适的电压，并降低因整流系统导致的波形畸变对电网的污染。以常见的12脉波整流机组中的整流变压器为例，它通常采用特殊的接线方式，如Dd0/d11或Yd11/d11接线。在这种接线方式下，变压器的阀侧（二次侧）会输出两组相位相差30°的电压，为后续的整流电路提供输入。例如，对于Dd0/d11接线的整流变压器，其一次侧绕组采用三角形连接，二次侧有两个绕组，一个为三角形连接，另一个为星形连接，且这两个绕组的线电压相位相差30°。这种特殊的接线方式使得整流变压器能够输出多相交流电，为实现多脉波整流创造了条件，有效减少了整流后直流电压中的谐波含量。硅整流器则利用了半导体的单向导电性，将整流变压器输出的交流电转换为直流电。常见的硅整流器由多个硅整流二极管组成，这些二极管按照特定的电路结构连接，如三相桥式整流电路、双反星形带平衡电抗器整流电路等。以三相桥式整流电路为例，它由六个硅整流二极管组成，其中三个二极管的阳极连接在一起，另外三个二极管的阴极连接在一起。在交流电压的正半周，阳极电位较高的三个二极管导通，电流从变压器的二次侧绕组流出，经过导通的二极管流向负载；在交流电压的负半周，阴极电位较低的三个二极管导通，电流从负载流回变压器的二次侧绕组。通过这种方式，交流电被转换为直流电，实现了电能的整流。硅整流器具有结构简单、可靠性高、成本较低等优点，在直流牵引供电系统中得到了广泛的应用。3.2.2多段线性化外特性模型为了准确描述整流机组的工作特性，建立多段线性化外特性模型是一种有效的方法。整流机组的输出特性受到多种因素的影响，如整流变压器的参数、硅整流器的特性以及负载的变化等，呈现出非线性的特点。通过对整流机组在不同工作条件下的输出特性进行分析，可以将其外特性曲线划分为多个线性段，每个线性段可以用一个简单的数学模型来描述。在建立多段线性化外特性模型时，通常会根据整流机组的实际工作情况，选择合适的分段点。例如，根据整流机组的额定电压、额定电流以及不同负载下的输出特性，将外特性曲线划分为三段：轻载段、额定负载段和过载段。在轻载段，整流机组的输出电压较为稳定，电流较小，此时可以用一个线性函数来描述其外特性，如U=U_0-I\timesr_1，其中U为输出电压，U_0为空载电压，I为输出电流，r_1为轻载段的等效内阻。在额定负载段，整流机组工作在最佳状态，输出电压和电流相对稳定，其外特性可以表示为U=U_{rated}-I\timesr_2，其中U_{rated}为额定电压，r_2为额定负载段的等效内阻。当整流机组处于过载段时，由于变压器的饱和以及硅整流器的非线性特性，输出电压会随着电流的增加而迅速下降，此时外特性可以用U=U_{max}-I\timesr_3来描述，其中U_{max}为最大输出电压，r_3为过载段的等效内阻。为了便于分析和计算，通常会采用戴维南等效电路来替代整流机组。戴维南等效电路由一个理想电压源和一个内阻串联组成，其等效电压源的电压等于整流机组的开路电压，等效内阻等于整流机组在特定条件下的输出电阻。通过将整流机组用戴维南等效电路代替，可以将复杂的整流机组模型简化为一个简单的电路模型，方便后续的电路分析和计算。例如，在计算直流牵引供电系统的短路电流时，可以将整流机组用戴维南等效电路代替，然后利用电路分析的方法求解短路电流。这种方法不仅能够提高计算效率，还能够更直观地反映整流机组在不同工作条件下的输出特性，为直流牵引供电系统的设计、运行和维护提供了有力的支持。3.3系统整体数学模型建立在构建直流牵引供电系统的整体数学模型时，需要将前面分析得到的各部分模型进行有机整合。首先，将直流牵引变电所的多段线性化外特性模型，通过戴维南等效电路，以理想电压源和内阻串联的形式融入到整个系统模型中。例如，对于某实际的直流牵引供电系统，假设其整流机组在额定负载段的戴维南等效电路中，理想电压源的电压为U_{rated}，内阻为r_2。接触网和走行轨的模型则基于之前计算得到的直流电阻和电感参数进行构建。根据基尔霍夫定律，在系统的任意节点，流入节点的电流之和等于流出节点的电流；对于任意闭合回路，各元件上的电压降之和为零。以一个简单的双边供电的直流牵引供电系统为例，假设在某一时刻，有列车在接触网的某位置取流，根据基尔霍夫电流定律，在牵引变电所与接触网的连接节点处，从牵引变电所流出的电流等于列车取流电流与流向其他部分的电流之和；根据基尔霍夫电压定律，在包含接触轨、走行轨和列车的闭合回路中，接触轨和走行轨上的电压降与列车两端的电压之和为零。将列车的负载模型也纳入系统整体模型中，列车的负载特性会随着列车的运行状态而变化，如列车的启动、加速、匀速行驶和制动等不同工况下，其功率需求和电流特性都有所不同。在建立系统整体数学模型时，需要考虑这些因素，通过相应的数学表达式来描述列车负载与系统其他部分的相互关系。通过以上步骤，将直流牵引供电系统的各个组成部分的模型进行整合，建立起系统整体数学模型。该模型能够全面、准确地描述直流牵引供电系统在不同运行工况下的电气特性，为后续的短路计算提供了坚实的基础。在实际应用中，可以根据具体的系统参数和运行条件，对该模型进行进一步的优化和调整，以提高其计算精度和实用性。四、直流牵引供电系统短路计算方法研究4.1短路计算的意义与内容短路计算在直流牵引供电系统中具有举足轻重的地位，它是确保系统安全稳定运行、实现设备合理选型以及继电保护准确整定的关键环节。在设备选型方面，短路计算结果是选择直流快速开关、熔断器等关键设备的重要依据。以直流快速开关为例，其额定电流、分断能力等参数的选择必须基于准确的短路电流计算结果。如果所选直流快速开关的分断能力不足，在短路故障发生时，将无法迅速切断故障电流，可能导致设备烧毁、火灾等严重后果，威胁整个供电系统的安全。若分断能力过大，则会造成设备成本增加，资源浪费。对于继电保护整定，短路计算更是不可或缺。继电保护装置的作用是在系统发生故障时，迅速、准确地动作，切除故障部分，保护系统的其他部分正常运行。通过短路计算，可以确定不同运行工况下的短路电流大小和分布情况，从而为继电保护装置设定合理的动作值和动作时间。例如，对于过电流保护，需要根据短路计算结果确定保护装置的动作电流，使其在正常运行时不动作，而在发生短路故障时能够迅速动作，切除故障线路。如果动作电流设定过小，可能导致保护装置误动作，影响系统的正常运行；若设定过大，则在发生短路故障时，保护装置可能无法及时动作，使故障范围扩大。直流牵引供电系统短路计算的内容丰富多样，主要包括以下几个方面：在正常情况下双边供电时，需计算各供电区间任一点的直流短路电流。双边供电是直流牵引供电系统常见的供电方式，通过计算各供电区间任一点的短路电流，可以全面了解系统在正常运行状态下的短路特性，为设备选型和继电保护整定提供基础数据。当任一中间牵引变电所解列时，由相邻牵引变电所构成大双边供电，此时需要计算这种特殊供电方式下区间任一点的直流短路电流。中间牵引变电所解列可能是由于设备故障、检修等原因引起的，计算大双边供电时的短路电流，有助于分析系统在这种异常情况下的运行稳定性，确保继电保护装置能够正确动作，保障系统的安全。端头牵引变电所解列时，由次端头牵引变电所单边供电，同样需要计算该区间任一点的直流短路电流。端头牵引变电所解列会改变系统的供电结构和电流分布，计算单边供电时的短路电流，对于评估系统在极端情况下的可靠性，制定相应的应急预案具有重要意义。4.2现有短路计算方法分析4.2.1电路图法电路图法是针对城市轨道交通直流牵引供电系统电源多、供电回路多、供电方式多、回路参数多的特点而发展起来的一种短路计算方法。该方法首先按照实际供电网络画出等效电路图，由于在直流系统中，电感对稳态短路电流的影响较小，因此在等效电路图中主要考虑电阻元件，将网络变换为只包含电阻的形式。例如，对于一个双边供电的直流牵引供电系统，其等效电路图会包含牵引变电所的内阻、接触网电阻、走行轨电阻以及短路点的位置等信息。在完成等效电路图的绘制和网络变换后，利用基本的电路定律，如欧姆定律和基尔霍夫定律进行计算。欧姆定律描述了电流、电压和电阻之间的关系，即I=\frac{U}{R}，其中I为电流，U为电压，R为电阻。基尔霍夫定律包括电流定律和电压定律，电流定律指出在电路的任一节点，流入节点的电流之和等于流出节点的电流；电压定律则表明在电路的任一闭合回路中，各段电压的代数和等于零。通过应用这些定律，可以列出关于电路中各节点电流和各回路电压的方程组，进而求解出短路电流。然而，电路图法存在一定的局限性。它只能计算稳态短路电流，无法计算供电回路的时间常数和短路电流上升率di/dt。在实际的直流牵引供电系统中，短路故障发生后，电流并不是瞬间达到稳态值，而是有一个上升的过程，短路电流上升率对于评估系统的暂态特性和保护装置的动作时间具有重要意义。时间常数反映了电路中电感和电阻的综合作用，影响着短路电流的上升速度和暂态过程。由于电路图法忽略了这些因素，在分析系统的暂态性能和设计快速响应的保护装置时，其计算结果的参考价值有限。例如，在设计直流快速开关时，需要准确知道短路电流的上升率，以确保开关能够在规定的时间内切断故障电流，而电路图法无法提供这方面的信息，可能导致开关选型不合理，影响系统的安全运行。4.2.2示波图法示波图法是一种基于工程实践的短路计算方法，其计算原理是通过对现场短路试验所拍摄的示波图进行数理分析，从而计算出相关参数。在进行现场短路试验时，会在直流牵引供电系统的特定位置设置测量设备，如电流互感器、电压互感器等，以获取短路发生时的电流和电压信号。这些信号被记录在示波图上，示波图直观地展示了短路电流和电压随时间的变化曲线。以某地铁线路的现场短路试验为例，在试验过程中，当短路故障发生时，测量设备迅速捕捉到电流和电压的变化，并将其记录在示波图上。通过对示波图的分析，可以获取短路电流的初始值、上升过程以及稳态值等信息。根据电路理论，短路电流的上升过程符合一定的数学规律，通过对示波图上电流曲线的分析，可以利用相关的数学公式计算出短路电流的上升率di/dt。例如，根据短路电流上升的指数函数关系i=I_k(1-e^{-\frac{t}{\tau}})（其中i为短路电流随时间t的变化值，I_k为稳态短路电流，\tau为时间常数），通过在示波图上选取合适的时间点和对应的电流值，代入公式中进行计算，就可以得到短路电流的上升率和时间常数。示波图法的优点在于它能够计算直流回路的多个重要参数，包括稳态短路电流、回路时间常数和短路电流上升率，这是电路图法所无法比拟的。通过准确获取这些参数，可以更全面地了解直流牵引供电系统在短路故障发生时的暂态特性，为系统的设计、保护装置的整定以及故障分析提供更丰富、准确的数据支持。然而，示波图法也存在一些不足之处。现场短路试验的实施难度较大，需要投入大量的人力、物力和时间，而且试验过程中存在一定的安全风险。试验条件的控制也较为困难，不同的试验条件可能会导致测量结果的差异，从而影响计算结果的准确性。由于示波图法依赖于现场试验数据，对于尚未建设或改造的系统，无法提前进行短路计算分析，具有一定的局限性。4.3提出改进的短路计算方法针对现有短路计算方法存在的局限性，本文提出一种改进的基于节点电压法和动态时间步长的短路计算方法，该方法综合考虑系统的动态特性和元件参数的变化，能够更准确地计算短路电流。改进思路主要体现在以下几个方面：充分考虑系统中电感和电容等动态元件对短路电流的影响，摒弃传统方法中对电感和电容的简化处理。在直流牵引供电系统中，虽然电感对稳态短路电流的影响相对较小，但在短路暂态过程中，电感和电容的作用不可忽视。通过建立包含电感和电容的电路模型，能够更真实地反映系统在短路瞬间的电气特性。考虑系统运行工况的动态变化，如列车的启动、加速、制动等过程，这些过程会导致系统的负载特性发生变化，进而影响短路电流的大小和变化规律。传统的短路计算方法往往假设系统处于稳态运行，无法准确反映这些动态变化对短路电流的影响。新方法通过实时监测系统的运行工况，动态调整计算参数，以适应不同工况下的短路计算需求。引入动态时间步长的概念，根据短路电流的变化率自动调整计算时间步长。在短路发生初期，电流变化迅速，采用较小的时间步长能够更精确地捕捉电流的变化；随着短路过程的发展，电流逐渐趋于稳定，适当增大时间步长可以提高计算效率，减少计算时间。具体计算方法如下：基于节点电压法，对直流牵引供电系统的电路进行分析。将系统中的各个元件，如牵引变电所、接触网、走行轨、列车等，用相应的电路模型表示，并将其连接成完整的电路网络。在电路网络中，定义各个节点的电压和支路的电流，根据基尔霍夫定律列出节点电压方程和支路电流方程。对于一个包含n个节点的直流牵引供电系统电路网络，其节点电压方程可以表示为YV=I，其中Y为节点导纳矩阵，V为节点电压向量，I为节点注入电流向量。节点导纳矩阵Y中的元素Y_{ij}表示节点i和节点j之间的导纳，当i=j时，Y_{ii}为节点i的自导纳，等于连接到该节点的所有支路导纳之和；当i\neqj时，Y_{ij}为节点i和节点j之间的互导纳，若节点i和节点j之间有支路直接相连，则Y_{ij}等于该支路导纳的相反数，否则Y_{ij}=0。在计算过程中，采用动态时间步长技术。根据短路电流的变化率\frac{di}{dt}来调整时间步长\Deltat。当\frac{di}{dt}较大时，说明短路电流变化迅速，此时减小时间步长\Deltat，例如取\Deltat=\frac{\Deltat_{0}}{k}（其中\Deltat_{0}为初始时间步长，k>1为调整系数），以提高计算精度；当\frac{di}{dt}较小时，说明短路电流趋于稳定，适当增大时间步长\Deltat，如取\Deltat=k\Deltat_{0}，以加快计算速度。在每个时间步长内，根据系统的当前状态和参数，更新节点导纳矩阵Y和节点注入电流向量I，然后求解节点电压方程YV=I，得到各节点的电压值。根据节点电压和支路导纳，计算各支路的电流，从而得到短路电流的大小和分布。在实际计算中，利用计算机编程实现上述算法。采用MATLAB等编程语言，编写短路计算程序。程序首先读取直流牵引供电系统的参数，如牵引变电所的内阻、接触网和走行轨的电阻和电感、列车的负载特性等，并根据这些参数初始化节点导纳矩阵和节点注入电流向量。在计算过程中，通过循环迭代的方式，按照动态时间步长逐步计算短路电流的变化。每计算一步，判断是否满足计算结束条件，如短路电流已趋于稳定或达到设定的计算时间。若满足条件，则结束计算，输出短路电流的计算结果；否则，继续进行下一步计算。通过这种改进的短路计算方法，能够更全面、准确地计算直流牵引供电系统在不同工况下的短路电流，为系统的设计、保护装置的整定以及故障分析提供更可靠的依据。五、案例分析与仿真验证5.1实际工程案例选取本文选取了某城市地铁线路的直流牵引供电系统作为实际工程案例，该线路采用DC1500V电压等级，双边供电方式，具有典型性和代表性。线路全长30km，共设20座车站，其中包含10座地下站和10座高架站。牵引变电所分布在沿线，相邻牵引变电所之间的平均距离约为3km。该线路的牵引变电所采用12脉波整流机组，由两台整流变压器和六台整流器组成。整流变压器的容量为3150kVA，一次侧电压为35kV，二次侧电压为10.5kV。整流器采用晶闸管整流，通过控制晶闸管的触发角来调节输出电压。牵引变电所的直流输出额定电压为1500V，额定电流为3000A。接触网采用柔性接触网，导线材质为铜合金，截面积为120mm²。接触网的悬挂高度为4.5m，跨距为5m。走行轨采用60kg/m的钢轨，其直流电阻在20℃时为0.15Ω/km，电感为0.7mH/km。回流线采用截面积为150mm²的铜导线，其直流电阻为0.12Ω/km，电感为0.6mH/km。该地铁线路的列车采用4动2拖的编组形式，每节动车的牵引功率为200kW，每节拖车的牵引功率为100kW。列车的最大运行速度为80km/h，平均运行速度为35km/h。在正常运行情况下，列车的启动电流为额定电流的1.5倍，加速电流为额定电流的1.2倍，匀速运行电流为额定电流，制动电流为额定电流的0.8倍。通过对该实际工程案例的系统参数进行详细分析，为后续的短路计算和仿真验证提供了真实可靠的数据基础，有助于更准确地评估本文所提出的数学模型和短路计算方法在实际工程中的应用效果。5.2基于数学模型的短路计算运用前文建立的直流牵引供电系统数学模型和改进的短路计算方法，对不同短路情况下的电流进行计算。以所选实际工程案例的地铁线路为例，分别计算正常双边供电、中间牵引变电所解列大双边供电以及端头牵引变电所解列单边供电时的短路电流。在正常双边供电情况下，假设短路点位于距离某一牵引变电所1km处，利用改进的短路计算方法，考虑系统中电感、电容等动态元件以及列车运行工况的动态变化，通过节点电压法和动态时间步长技术进行计算。根据系统参数，初始化节点导纳矩阵和节点注入电流向量，在计算过程中，按照动态时间步长逐步更新计算参数。经过计算，得到该短路点的短路电流初始值为I_{01}=4500A，短路电流上升率\frac{di}{dt}=1500A/ms，在短路发生0.1s后，短路电流达到稳态值I_{s1}=5000A。当中间牵引变电所解列，形成大双边供电时，假设短路点位于两个相邻牵引变电所中间位置。此时，由于供电结构的变化，系统的参数和电流分布发生改变。重新计算节点导纳矩阵和节点注入电流向量，按照改进的短路计算方法进行计算。计算结果表明，短路电流初始值为I_{02}=4800A，短路电流上升率\frac{di}{dt}=1600A/ms，在短路发生0.1s后，短路电流稳态值为I_{s2}=5200A。在端头牵引变电所解列，由次端头牵引变电所单边供电时，设短路点距离次端头牵引变电所2km处。同样，根据系统的实际情况，调整计算参数，运用改进的短路计算方法进行计算。得到短路电流初始值为I_{03}=4200A，短路电流上升率\frac{di}{dt}=1400A/ms，在短路发生0.1s后，短路电流稳态值为I_{s3}=4800A。通过对不同短路情况下电流的计算，可以清晰地了解系统在各种故障条件下的电气特性，为后续的仿真验证和实际工程应用提供了具体的数据支持，有助于评估系统的安全性和可靠性，为直流快速开关的选型和直流保护装置的整定提供依据。5.3仿真模型建立与结果分析利用MATLAB的SIMULINK模块建立直流牵引供电系统的仿真模型，以进一步验证前面计算结果的准确性和所提方法的有效性。在SIMULINK中，直流牵引供电系统的各组成部分都有对应的模块，通过合理设置参数和连接方式，构建出完整的系统模型。首先，根据实际工程案例的参数，对牵引变电所中的整流机组进行建模。以12脉波整流机组为例，在SIMULINK中使用相应的电力电子模块搭建整流电路，设置整流变压器的参数，如一次侧和二次侧的电压、匝数比等，以及晶闸管的触发角等参数，使其符合实际系统的工作特性。利用电压源和内阻模块来模拟整流机组的戴维南等效电路，其中电压源的电压设置为整流机组的开路电压，内阻设置为根据多段线性化外特性模型确定的等效内阻。对于接触网和走行轨，根据前面计算得到的直流电阻和电感参数，在SIMULINK中使用电阻和电感模块进行建模。将接触网和走行轨的电阻和电感按照实际线路的分布情况进行连接，考虑到线路的长度和分段情况，采用多个电阻和电感模块串联或并联的方式来模拟实际的电气特性。例如，对于一条较长的接触网线路，可以将其分为多个小段，每段用一个电阻和电感模块来表示，然后将这些模块依次串联起来，以更准确地模拟接触网的电阻和电感分布。列车的负载模型根据其运行工况进行设置。在SIMULINK中，通过设置电流源或功率源模块来模拟列车在不同运行工况下的负载特性。例如，在列车启动时，设置电流源的电流值为启动电流，根据实际情况，启动电流通常为额定电流的1.5倍左右；在列车加速时，电流源的电流值按照加速电流的规律变化，如加速电流为额定电流的1.2倍；在列车匀速运行时，电流源的电流值保持为额定电流；在列车制动时，电流源的电流值设置为制动电流，如制动电流为额定电流的0.8倍。通过这种方式，能够真实地模拟列车在不同运行工况下对直流牵引供电系统的影响。在建立好仿真模型后，设置不同的短路故障场景，如正常双边供电时不同位置的短路、中间牵引变电所解列大双边供电时的短路以及端头牵引变电所解列单边供电时的短路等，进行仿真分析。以正常双边供电时距离某一牵引变电所1km处短路为例，运行仿真模型，得到短路电流随时间的变化曲线。将仿真结果与前面基于数学模型的短路计算结果进行对比，如图1所示。从图中可以看出，仿真结果与计算结果在趋势上基本一致，短路电流的初始值、上升率以及稳态值都较为接近。仿真得到的短路电流初始值为I_{01sim}=4480A，与计算值I_{01}=4500A相比，相对误差约为0.44\%；短路电流上升率\frac{di}{dt}_{sim}=1490A/ms，与计算值\frac{di}{dt}=1500A/ms相比，相对误差约为0.67\%；在短路发生0.1s后，仿真得到的短路电流稳态值为I_{s1sim}=4980A，与计算值I_{s1}=5000A相比，相对误差约为0.4\%。对于中间牵引变电所解列大双边供电和端头牵引变电所解列单边供电时的短路情况，同样进行仿真分析并与计算结果对比。在中间牵引变电所解列大双边供电时，仿真得到的短路电流初始值为I_{02sim}=4780A，与计算值I_{02}=4800A相比，相对误差约为0.42\%；短路电流上升率\frac{di}{dt}_{sim}=1590A/ms，与计算值\frac{di}{dt}=1600A/ms相比，相对误差约为0.63\%；短路电流稳态值为I_{s2sim}=5180A，与计算值I_{s2}=5200A相比，相对误差约为0.38\%。在端头牵引变电所解列单边供电时，仿真得到的短路电流初始值为I_{03sim}=4180A，与计算值I_{03}=4200A相比，相对误差约为0.48\%；短路电流上升率\frac{di}{dt}_{sim}=1390A/ms，与计算值\frac{di}{dt}=1400A/ms相比，相对误差约为0.71\%；短路电流稳态值为I_{s3sim}=4780A，与计算值I_{s3}=4800A相比，相对误差约为0.42\%。通过以上对比分析可以看出，仿真结果与基于数学模型的短路计算结果具有较高的一致性，验证了所建立的数学模型和改进