一种适用于断路器的热点温度分析方法与流程

本发明涉及电力系统技术领域，尤其涉及一种适用于断路器的热点温度分析方法。

背景技术：

温升及热点温度分析在研究区域电网动态增容方面是一项关键的研究要素，断路器作为区域电网的重要变电设备，在做系统增容运行之前，必须精确地计算其热点温度，保证关键部位温度不超过国家标准所规定的数值。瓷柱式断路器是一种常见的断路器，瓷柱式短路器的结构包括灭弧室、瓷套管、法兰、上法兰接线板以及下法兰接线板，瓷柱式断路器采用sf6气体作为灭弧介质和绝缘介质。

目前，断路器的热点分析方法主要使用数值仿真计算方法和热路模型计算方法，热路模型计算方法较数值仿真分析方法，模型简单，计算量较少，能够在采集数据之后快速计算，但不足的是当前建立的热路模型当中都没有考虑环境温度、风速、sf6气压等综合因素的影响，导致计算结果不够准确。

技术实现要素：

针对现有的热路模型计算方法存在计算结果不准确的问题，本发明提供了一种提高计算结果准确率的适用于断路器的热点温度分析方法。

本发明采用如下技术方案：

一种适用于断路器的热点温度分析方法，所述断路器为瓷柱式断路器，所述热点温度分析方法包括：

步骤s1、采集所述断路器所处环境的环境温度、环境风速、环境湿度以及所述断路器的接触电阻；

步骤s2、根据所述环境温度、所述环境风速、所述环境湿度以及所述接触电阻构建一热路模型，根据所述热路模型处理得到所述断路器的热点温度。

优选的，所述瓷柱式断路器包括sf6灭弧室、瓷套管、法兰、上法兰接线板以及下法兰接线板，所述瓷柱式断路器采用sf6气体作为灭弧介质和绝缘介质。

优选的，所述步骤s2包括：

步骤s21、采用一第一处理单元处理得到所述断路器的总热损耗；

步骤s22、采用一第二处理单元处理得到所述断路器的总热阻；

步骤s23、采用一第三处理单元，所述第三处理单元根据所述环境温度、所述环境风速、所述环境湿度、所述接触电阻、所述总热损耗以及所述总热阻构建一所述热路模型；

步骤s24、根据所述热路模型处理得到所述断路器的热点温度。

优选的，所述步骤s22包括：

步骤s221、采用一第四处理单元处理得到所述瓷套管内表面的第一热阻；

步骤s222、采用一第五处理单元处理得到所述瓷套管外表面的第二热阻；

步骤s223、采用一第六处理单元处理得到所述瓷套管外表面到外部空气的第三热阻；

步骤s224、采用一第七处理单元处理得到所述上法兰接线板和所述法兰到所述外部空气的第四热阻；

步骤s225、采用一第八处理单元处理得到所述下法兰接线板和所述法兰到所述外部空气的第五热阻；

步骤s226、采用一第九处理单元，所述第九处理单元根据所述第一热阻、所述第二热阻、所述第三热阻、所述第四热阻及所述第五热阻处理得到所述总热阻。

优选的，所述步骤s223中，所述第六处理单元采用下述公式处理得到所述第三热阻：

其中，

rce-air用于表示所述第三热阻；

π用于表示圆周率；

α用于表示所述瓷套管外表面等效面积比例系数；

d用于表示所述瓷套管的外径；

l用于表示一预设的参数值；

k用于表示所述断路器所处的外部空气的对流换热系数。

优选的，所述步骤s223中，采用一第十处理单元，根据一函数模型对所述环境温度、环境风速以及所述环境湿度进行处理以得到所述对流换热系数；

获取所述函数模型的具体步骤如下：

步骤a1、获取所述环境温度、所述环境风速及所述环境湿度；

步骤a2、通过预设的仿真计算规则处理得到换热热量；

步骤a3、根据所述换热热量和预设的换热系数计算规则处理得到所述对流换热系数；

步骤a4、将所述环境温度、所述环境风速及所述环境湿度作为输入向量，将所述对流换热系数作为输出向量，通过所述输入向量和所述输出向量构建一数据集；

步骤a5、将所述数据集分为两组数据，一组数据作为训练样本，另一组数据作为预测样本，通过所述训练样本处理得到所述函数模型，通过所述预测样本对所述函数模型进行修正以得到修正后的所述函数模型。

优选的，所述步骤a3中，采用下述公式处理得到所述对流换热系数：

q′＝kα·δt；

其中，

q′用于表示所述换热热量；

k用于表示所述对流换热系数；

α用于表示所述瓷套管外表面等效面积比例系数；

δt用于表示所述瓷套管外表面与所述环境温度之间的差值。

优选的，所述步骤a4中，将所述输入向量和所述输出向量进行归一化处理后构建所述数据集。

优选的，所述步骤a5包括：

步骤a51、根据所述训练样本确定所述对流换热系数的最优区间；

步骤a52、确定所述对流换热系数的取值；

步骤a53、判断所述对流换热系数是否满足精度：

若判断结果为是，则转步骤a55；

若判断结果为否，则转步骤a54；

步骤a54、判断是否达到最大迭代次数：

若判断结果为是，则转步骤a55；

若判断结果为否，则转步骤a52；

步骤a55、根据所述对流换热系数的最优区间处理得到所述函数模型；

步骤a56、根据所述预测样本中的所述环境温度、所述环境风速及所述环境湿度和所述函数模型处理得到所述对流换热系数；

步骤a57、根据所述对流换热系数处理得到待验证换热热量；

步骤a58、判断所述待验证换热热量和所述换热热量之间的误差的绝对值是否小于换热热量的10％：

若判断结果为正，则所述函数模型为修正后的所述函数型；

若判断结果为否，则转步骤a51。

本发明的有益效果是：考虑环境温度、环境风速、环境湿度、接触电阻等因素的影响建立断路器热路模型，使得计算结果更加精确。

附图说明

图1为本发明的一种优选实施例中，适用于断路器的热点温度分析方法的流程图；

图2为本发明的一种优选实施例中，短路器的热路模型；

图3为本发明的一种优选实施例中，步骤s2的流程图；

图4为本发明的一种优选实施例中，步骤s22的流程图；

图5为本发明的一种优选实施例中，获取函数模型的流程图；

图6为本发明的一种优选实施例中，步骤a5的流程图。

具体实施方式

需要说明的是，在不冲突的情况下，下述技术方案，技术特征之间可以相互组合。

下面结合附图对本发明的具体实施方式作进一步的说明：

如图1所示，一种适用于断路器的热点温度分析方法，上述断路器为瓷柱式断路器，上述热点温度分析方法包括：

步骤s1、采集上述断路器所处环境的环境温度、环境风速、环境湿度以及上述断路器的接触电阻；

上述瓷柱式断路器包括圆柱形导体、sf6灭弧室、瓷套管、法兰、上法兰接线板以及下法兰接线板，上述瓷柱式断路器采用sf6气体作为灭弧介质和绝缘介质，圆柱形导体、sf6灭弧室、瓷套管、法兰、上法兰接线板以及下法兰接线板均为现有技术。

在本实施例中，考虑环境温度、环境风速、环境湿度、接触电阻等因素的影响建立断路器热路模型，使得计算结果更加精确。

本发明较佳的实施例中，如图2-3所示，其中，图2为断路器热路模型示意图，断路器的热量流散路径如图2所示；图3为步骤s2的流程图。

上述步骤s2包括：

步骤s21、采用一第一处理单元处理得到上述断路器的总热损耗q；

步骤s22、采用一第二处理单元处理得到上述断路器的总热阻r；

步骤s23、采用一第三处理单元，上述第三处理单元根据上述环境温度、上述环境风速、上述环境湿度、上述接触电阻、上述总热损耗q以及上述总热阻r构建一上述热路模型；

步骤s24、根据上述热路模型处理得到上述断路器的热点温度。

在本实施例中，可以根据热路模型处理得到断路器的最高热点温度，第一处理单元、第二处理单元、第三处理单元可集成于同一硬件设备(如同一处理器)中。

本发明较佳的实施例中，采用下述公式(1)表示上述热路模型：

ts＝qr-te；(1)

其中，

ts用于表示上述热点温度；

q用于表示上述断路器的总热损耗；

r用于表示上述断路器的总热阻；

te用于表示上述环境温度。

如图4所示，本发明较佳的实施例中，上述步骤s22包括：

步骤s221、采用一第四处理单元处理得到上述瓷套管内表面的第一热阻rsf6；

步骤s222、采用一第五处理单元处理得到上述瓷套管外表面的第二热阻rce；

步骤s223、采用一第六处理单元处理得到上述瓷套管外表面到外部空气的第三热阻rce-air；

步骤s224、采用一第七处理单元处理得到上述上法兰接线板和上述法兰到上述外部空气的第四热阻rt-air；

步骤s225、采用一第八处理单元处理得到上述下法兰接线板和上述法兰到上述外部空气的第五热阻rb-air；

步骤s226、采用一第九处理单元，上述第九处理单元根据上述第一热阻rsf6、上述第二热阻rce、上述第三热阻rce-air、上述第四热阻rt-air及上述第五热阻rb-air处理得到上述总热阻r；

第一到第九处理单元均可集成在同一硬件设备(如同一处理器)中其中，采用下述公式(2)表示r：

r＝(rsf6+rce+rce-air)||rt-air||rb-air；(2)

其中，

rsf6用于表示上述sf6气体到上述瓷套管内表面的热阻；

rce用于表示上述瓷套管内表面到上述瓷套管外表面的热阻；

rce-air用于表示上述瓷套管外表面到外部空气的热阻；

rt-air用于表示上述上法兰接线板和上述法兰到上述外部空气的热阻；

rb-air用于表示上述下法兰接线板和上述法兰到上述外部空气的热阻。

进一步的，采用下述公式(3)表示rsf6：

rsf6＝rsf6-c||rsf6-r；(3)

其中，

rsf6-c用于表示上述sf6气体对流换热热阻；

rsf6-r用于表示上述sf6气体到上述瓷套管内表面的辐射热阻。

进一步的，采用下述公式(4)表示rsf6-c：

其中，

q1用于表示对流换热过程中所交换的热量；

t2用于表示圆形导电体温度，单位为k；

t1用于表示上述瓷套管内表面的温度，单位为k。

进一步的，采用下述公式(5)表示q1：

其中，

g＝9.8m/s2，为万有引力常数；

β用于表示上述sf6气体在特定温度下的热膨胀系数，单位为k^－1；

p用于表示上述sf6气体的气体压力；

μ用于表示上述sf6气体在特定温度下的动力粘度，单位为kg/(m·s)；

cp用于表示上述sf6气体在特定温度下的比热容，单位为j/(kg·k)；

v用于表示上述sf6气体在特定温度下的运动粘度，单位为m²/s；

d1用于表示上述瓷套管的内径，单位为m；

d2用于表示表示上述圆形导电体的外径，单位为m；

l用于表示上述灭弧室的长度，单位为m；

λ用于表示上述sf6气体在特定温度下的导热系数，单位为w/(m·k)；

t2用于表示圆形导电体温度，单位为k；

t1用于表示上述瓷套管内表面的温度，单位为k。

进一步的，采用下述公式(6)表示rsf6-r：

其中，

t1用于表示上述瓷套管内表面的温度，单位为k；

t2用于表示上述圆形导电体的温度，单位为k；

q2用于表示辐射换热过程中所交换的热量。

进一步的，采用下述公式(7)表示上述q2：

其中：

π用于表示圆周率；

t1用于表示上述瓷套管内表面的温度，单位为k；

t2用于表示上述圆形导电体的温度，单位为k；

d1用于表示上述瓷套管的内径，单位为m；

d2用于表示表示上述圆形导电体的外径，单位为m；

ε1用于表示上述瓷套管内表面的黑度；

ε2用于表示上述圆形导电体表面的黑度。

进一步的，采用下述公式(8)表示rce：

其中，

π用于表示圆周率；

r1用于表示上述瓷套管内表面的半径，单位为m；

r2用于表示上述瓷套管外表面的半径，单位为m；

λt用于表示上述瓷套管的陶瓷热导率，单位为w/(m·k)；

l用于表示上述瓷套管的套管长度，单位为m。

本发明较佳的实施例中，上述步骤s223中，上述第六处理单元采用下述公式(9)处理得到上述第三热阻rce-air：

其中，

π用于表示圆周率；

α用于表示上述瓷套管外表面等效面积比例系数；

d用于表示上述瓷套管的外径，单位为m；

l用于表示一预设的参数值；

k用于表示上述外部空气的对流换热系数；

α、d、l为固定值，rce-air的取值主要取决于k，k随环境温度、环境风速及环境湿度的变化而变化。

如图5所示，本发明较佳的实施例中，上述步骤s223中，采用一第十处理单元，根据一函数模型对上述环境温度、环境风速以及上述环境湿度进行处理以得到上述对流换热系数；

获取上述函数模型的具体步骤如下：

步骤a1、获取上述环境温度te、上述环境风速v及上述环境湿度s；

步骤a2、通过预设的仿真计算规则处理得到换热热量q’；

步骤a3、根据上述换热热量q’和预设的换热系数计算规则处理得到上述对流换热系数k；

步骤a4、将上述环境温度te、上述环境风速v及上述环境湿度s作为输入向量，将上述对流换热系数k作为输出向量，通过上述输入向量和上述输出向量构建一数据集；

步骤a5、将上述数据集分为两组数据，一组数据作为训练样本，另一组数据作为预测样本，通过上述训练样本处理得到上述函数模型，通过上述预测样本对上述函数模型进行修正以得到修正后的上述函数模型；利用数值方法对函数模型进行修正，进而对热路模型进行修正，使得计算结果更加精确。

上述步骤a4中，将上述输入向量和上述输出向量进行归一化处理后构建上述数据集。

在本实施例中，将环境温度te、环境风速v、环境湿度s作为输入向量，将对流换热系数k作为输出向量，通过支持向量机泛得出k与te、v、s的函数关系。其中，环境温度te、环境风速v、环境湿度s可通过微气象站测试得出，q’可通过仿真计算得出。

进一步的，上述步骤a3中，采用下述公式(10)处理得到上述对流换热系数：

q′＝kα·δt；(10)

其中，

q′用于表示上述换热热量；

k用于表示上述对流换热系数；

α用于表示上述瓷套管外表面等效面积比例系数；

δt用于表示上述瓷套管外表面与上述环境温度之间的差值。

进一步的，采用下述公式(10)表示k：

其中，

用于表示上述瓷套管外表面与上述外部空气对流的平均努塞尔数；

λair用于表示上述外部空气的空气导热率；

te用于表示上述环境温度；

v用于表示上述sf6气体在特定温度下的运动粘度；

s用于表示上述环境湿度；

b用于表示一预设的参数值。

进一步的，采用下述公式(12)表示上述

其中，

c用于表示一预设的参数值；

red用于表示雷诺数；

pr用于表示普朗特数。

进一步的，采用下述公式(13)表示上述red：

其中，

ρair用于表示上述外部空气的空气密度，单位为kg/m³；

v用于表示上述环境风速，单位为m/s；

d用于表示一预设的参数值；

μair用于表示一预设的参数值。

在不同温度下，空气的动力粘度不同：

当t<283时，μair＝14.16×10^(-6)，pr＝0.705，λair＝2.51×10^(-2)，ρair＝1.247；

当t<293时，μair＝15.06×10^(-6)，pr＝0.703，λair＝2.59×10^(-2)，ρair＝1.205；

当t<303时，μair＝16.00×10^(-6)，pr＝0.701，λair＝2.67×10^(-2)，ρair＝1.165；

当t<313时，μair＝16.96×10^(-6)，pr＝0.699，λair＝2.76×10^(-2)，ρair＝1.128；

除此之外，μair＝16.96×10^(-6)，pr＝0.698，λair＝2.83×10^(-2)，ρair＝1.1；

其中，c、n可以查表得到，red在数值在不同范围时，c、n的取值不同。

如图6所示，本发明较佳的实施例中，上述步骤a5包括：

步骤a51、根据上述训练样本确定上述对流换热系数的最优区间；

步骤a52、确定上述对流换热系数的取值；

步骤a53、判断上述对流换热系数是否满足精度(此处是否满足精度的判断条件为对流换热系数是否接近无穷小，若是则，满足精度，若否，则不满足精度)：

若判断结果为是，则转步骤a55；

若判断结果为否，则转步骤a54；

步骤a54、判断是否达到最大迭代次数：

若判断结果为是，则转步骤a55；

若判断结果为否，则转步骤a52；

步骤a55、根据上述对流换热系数的最优区间处理得到上述函数模型；

步骤a56、根据上述预测样本中的上述环境温度、上述环境风速及上述环境湿度和上述函数模型处理得到上述对流换热系数；

步骤a57、根据上述对流换热系数处理得到待验证换热热量；

步骤a58、判断上述待验证换热热量和上述换热热量之间的误差的绝对值是否小于换热热量的10％：

若判断结果为正，则上述函数模型为修正后的上述函数型；

若判断结果为否，则转步骤a51。

本发明较佳的实施例中，rt-air为上述上法兰接线板及上述法兰表面对流换热热阻与辐射换热热阻并联，分别采用下述公式(14)和(15)表示上述法兰表面对流换热热阻和上述辐射换热热阻：

其中，

r11用于表示上述对流换热热阻；

π用于表示圆周率；

d用于表示一预设的参数值；

k1用于表示散热系数；

r12用于表示上述辐射换热热阻；

ts用于表示上述热点温度；

te用于表示上述环境温度；

qt-air用于表示上法兰向空气辐射换热过程中交换的热量。

进一步的，采用下述公式(16)表示上述k1：

其中，

用于表示上述法兰表面对空气对路的平均努塞尔数；

qt-air用于表示上法兰向空气辐射换热过程中交换的热量；

ε用于表示接线板与上述法兰外表面发射率；

a1用于表示等效散热面积。

本发明较佳的实施例中，rb-air为上述上法兰接线板及上述法兰表面对流换热热阻r21与辐射换热热阻r22并联，分别采用下述公式(19)和(20)表示r21和r22；

其中，

qb-air为下法兰向空气辐射换热过程中交换的热量；

a2为等效散热面积。

本发明较佳的实施例中，采用下述公式(22)表示上述q：

q＝i²r；(22)

其中，

r≈rc+rb；(23)

rb＝r20[1+4.33(te-20)]；(24)

rc为上述接触电阻；

rb为导体电阻；

r20为导体在20℃的直流电阻；

i为当前运行电流。

本申请中的所有处理单元均可集成在同一硬件设备(如同一处理器)中，上述所有参数的计算均可通过该硬件设备处理得到。

通过说明和附图，给出了具体实施方式的特定结构的典型实施例，基于本发明精神，还可作其他的转换。尽管上述发明提出了现有的较佳实施例，然而，这些内容并不作为局限。

对于本领域的技术人员而言，阅读上述说明后，各种变化和修正无疑将显而易见。因此，所附的权利要求书应看作是涵盖本发明的真实意图和范围的全部变化和修正。在权利要求书范围内任何和所有等价的范围与内容，都应认为仍属本发明的意图和范围内。