一种基于磁‑热‑流耦合的变压器热稳定计算方法与流程

本发明属于电力变压器技术领域，涉及的是一种基于磁-热-流耦合的变压器热稳定计算方法。

背景技术：

变压器是电力系统能量传输和转换的关键枢纽设备。在运行工程中变压器的高损耗所引起的变压器的局部过热和温升过高，将会直接影响变压器的绝缘特性和使用寿命。

现阶段研究变压器温度场的方法主要有经验公式法、热路模型法和数值分析法。经验公式法可快速得到变压器的热点温度，在常规变压器的实际运行中具有一定的指导意义，但不适用于特殊结构的变压器，如立体卷芯变压器和新型变阻抗变压器。热路模型法以热电类比法为基础，用电路模型的方法简化变压器内部传热过程，将变压器内热交换过程简化为热路。该方法可直接得到反映变压器热交换物理过程的热路模型，并得到变压器热点温度的计算公式，计算结果的精度与热路参数的准确性密切相关。目前传统的热点温升的计算研究均假设其绕组热源的产热是均匀分布，实际上各个线饼的位置不同，每匝线圈的产热量并不相同。数值计算法是利用电磁场、传热学以及流体力学原理，通过建立变压器实体模型，研究绕组、铁心和变压器油之间的热量传递问题，并作适当简化，通过三个物理场之间的相互关系，联立求解微分方程组，得变压器温度场分布，变压器电磁损耗的计算直接影响温度场计算的准确性。

目前大多数的数值计算方法都是采用经验公式计算损耗，计算存在一定误差，其次磁-热-流三种物理场耦合分析时，将材料参数设置为常数，如铜的电导率、变压器油的比热容、热导率和粘度等，均未考虑材料属性对温度变化的依赖关系，从而温度场的计算精度无法保证。故对变压器本体在短路状况下的温度场分析时，有必要采用有限元方法计算变压器绕组、铁心和油箱等产生的电磁损耗，考虑温度变化对材料参数的影响，并且进一步采用磁-热-流耦合场的方法分别计算变压器本体的稳态温度场和暂态温度场分布。

技术实现要素：

本发明的目的在于提供一种基于磁-热-流耦合的变压器热稳定计算方法，可直接应用于实际系统中，具有较强的实用性。

本发明技术方案如下：

一种基于磁-热-流耦合的变压器热稳定计算方法，包括以下步骤：

步骤一：首先将变压器正常运行情况下的稳态电流值作为磁场的激励，进行稳态磁场计算，得出其各部件的单位体积的电磁损耗(电磁热)；

步骤二：将各部件单位体积的电磁热作为温度场分析的热源，温度场分布受变压器油流速的影响，油流速又受温度场的控制，同时计及变压器本体的物理参数，如油的热导率、动力粘度、比热容以及铜导线的电阻率依据温度变化的特性。如此循环迭代多次，直至温度和流速均不随时间变化，达到稳态，停止计算；

步骤三：将上述稳态计算的结果作为暂态初值，暂态磁场计算的瞬态激励取为流过变压器绕组的短路电流，磁场计算得到故障发生后2s内变压器本体产生的焦耳热；

步骤四：将单位体积的焦耳热作为温度场的热源，进行瞬态温度场的分析，同时考虑流速和温度的相互作用，计算过程中计及变压器本体的物理参数对温度的依赖关系，根据温度变化不断修正其材料参数；

步骤五：由电力变压器手册可知，计算变压器短路热稳定能力的持续时间是2s，故瞬态计算结束的判据是时间t>2s，最终得到变压器本体在短路故障发生2s后的温度场分布作为判别变压器本体热稳定性的依据。

本发明进一步的改进在于，在对变压器本体温度场进行模型仿真分析时，由于三绕组三相变压器绕组在结构上的对称性，且变压器的温度场沿铁心和绕组在圆周方向变化不大，所以取单相建立二维轴对称模型进行变压器本体的磁-热-流耦合场仿真计算。

本发明进一步的改进在于，将变压器本体正常工作电流的作为稳态磁场的激励，通过comsol有限元软件计算变压器本体发热和冷却过程，得到其热平衡状态的温度场。

本发明进一步的改进在于，通过采用有限元软件建立场路耦合模型。

与现有技术相比，本发明具有以下有益的技术效果：本发明采用多物理场全耦合的方式对变压器的温度场进行求解。首先将正常情况下的稳态电流值作为磁场的激励，计算得出变压器各部件单位体积的电磁热，再将其作为温度场分析的热源，考虑温度场和油流速的相互影响以及铜导线的电阻率等变压器本体的物理参数随温度变化的特性，计算得出达到稳态时温度和流速的具体数值。再以稳态计算的结果作为暂态计算的初值，暂态磁场计算的瞬态激励取为流过变压器绕组的短路电流，计算得到故障发生后2s内变压器本体产生的焦耳热。以单位体积的焦耳热作为温度场的热源，对瞬态温度场进行分析，并考虑流速和温度的相互作用；最终可得变压器本体在短路故障发生2s后的温度场分布作为判别变压器本体热稳定性的依据。

附图说明

图1为变压器本体对称模型；

图2为变压器本体2维温度场仿真模型；

图3为变压器本体稳态温度场计算流程图；

图4为变压器热交换示意图；

图5为当短路电流达到第1个峰值时的磁通密度分布；

图6为当短路电流达到第1个峰值时的电流密度分布。

具体实施方式

下面结合附图对本发明进行详细说明。

本发明为一种基于磁-热-流耦合的变压器热稳定计算方法，包括下述步骤：

因变压器本体短路温度场的初始状态是正常运行当变压器达到稳定状态的温度场，故首先对变压器正常工作状态下的稳态温度场进行计算，即以热平衡状态的温度场和速度场作为暂态耦合场的初值。再计算短路故障下变压器的暂态温升，最后得到变压器本体在短路故障发生2s后的温度场分布并以此作为判别变压器热稳定性的依据。

下面通过一个实例对本发明的步骤进行详细说明。

s01、由于三绕组三相变压器绕组在结构上的对称性，并且变压器的温度场沿铁心和绕组在圆周方向变化不大，所以取单相建立2维轴对称模型进行变压器本体的磁-热-流耦合场仿真计算，在极大的降低计算量的同时对计算精度亦不会有太大的影响，变压器本体对称模型如图1所示，仿真模型如图2所示。首先计算变压器本体正常工作下的稳态温度场，将正常情况下的稳态电流值作为磁场的激励进行稳态磁场计算，得到各部件的单位体积的电磁损耗。

s02、将各部件单位体积的电磁损耗作为温度场分析的热源，温度场分布受到变压器油流速的影响，油流速又受温度场的控制。与此同时，变压器本体的物理参数如油的热导率、动力粘度、比热容等也会随温度产生变化。经过不断循环迭代，温度和流速均不再变化，达到稳态，其过程如图3所示。

变压器油的主要物理参数如表1：

表1变压器油的主要物理参数

s03、将稳态计算的结果作为暂态初值，以流过变压器绕组的短路电流作为暂态磁场计算的瞬态激励，当短路电流达到第1个峰值t＝0.02s时的磁通密度和电流密度分布分别如图5和图6所示，磁通密度主要分布在高压绕组和中压绕组之间，和三维仿真结果一致，电流密度的分布由短路电流和每匝线圈的截面积共同决定。计算得到故障发生后2s内变压器本体产生的焦耳热。

s04、以单位体积的焦耳热作为温度场的热源，进行瞬态温度场的分析，同时考虑流速和温度的相互作用以及变压器本体的物理参数对温度的依赖关系，根据温度变化不断修正材料参数，最终可得变压器本体在短路故障发生2s后的温度场分布作为判别变压器本体热稳定性的依据。

本发明结合有限元软件，采用磁—热—流多物理场全耦合的方式对变压器的温度场进行求解，并以此为依据对变压器的热稳定性加以判断。本发明可直接应用于实际系统中，具有较强的实用性。