CN108038268B - 一种廊管内gil运行温度评估方法 - Google Patents
一种廊管内gil运行温度评估方法 Download PDFInfo
- Publication number
- CN108038268B CN108038268B CN201711149403.2A CN201711149403A CN108038268B CN 108038268 B CN108038268 B CN 108038268B CN 201711149403 A CN201711149403 A CN 201711149403A CN 108038268 B CN108038268 B CN 108038268B
- Authority
- CN
- China
- Prior art keywords
- conductor
- gil
- shell
- heat
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000004020 conductor Substances 0.000 claims abstract description 199
- 230000005855 radiation Effects 0.000 claims abstract description 57
- 238000004364 calculation method Methods 0.000 claims abstract description 38
- 239000003570 air Substances 0.000 claims description 52
- 238000012546 transfer Methods 0.000 claims description 51
- 230000017525 heat dissipation Effects 0.000 claims description 40
- 230000020169 heat generation Effects 0.000 claims description 35
- 230000002500 effect on skin Effects 0.000 claims description 24
- 238000004519 manufacturing process Methods 0.000 claims description 16
- 238000010438 heat treatment Methods 0.000 claims description 13
- 239000012080 ambient air Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 239000011343 solid material Substances 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 9
- 238000011156 evaluation Methods 0.000 abstract description 7
- 230000005540 biological transmission Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 238000004134 energy conservation Methods 0.000 description 2
- 239000012212 insulator Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 238000009421 internal insulation Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
一种廊管内GIL运行温度评估方法，建立GIL物理模型，计算外壳初始功率损耗，计算外壳与空气的热对流和热辐射量，计算导体与外壳间的热对流和热辐射量，本发明通过在Windows操作系统的MATLAB环境中设计的廊管内特高压GIL热特性MATLAB计算程序，利用初始工况下的GIL相关特性参数来对其稳定运行时的温度进行评估计算；本发明可评估计算不同的初始工况条件下廊管内特高压GIL的稳态运行温度；本发明评估计算速度快，精度高，可嵌入到GIL温度监测系统中，对提高GIL热可靠性具有重要意义。
Description
技术领域
本发明涉及特高压电气输电、气体绝缘输电线路(GIL)热特性、长距离特高压GIL设备热可靠性等领域，特别是涉及一种廊管内GIL运行温度评估方法。
背景技术
气体绝缘输电线路(GIL)是一种采用气体作为绝缘介质的长距离输电装置，与传统电缆输电相比，具有容量大，损耗低，使用寿命长等优点。常见的GIL故障如放电、绝缘子击穿等发生时都伴随有导体或外壳温度异常等现象。不仅如此，GIL母线温度过高还会发生母线拱顶，甚至盆式绝缘子破裂，盆子防水胶开裂等故障。GIL热特性的研究对提高GIL的热可靠性和保障GIL安全稳定运行具有重要意义。现阶段多采用有限元计算方法对GIL热特性进行研究，但是有限元仿真方法计算量大，一般不能作为子程序嵌入到GIL温度监测系统中。
因此，如何创新的提出一种快速有效的廊管内GIL运行温度评估方法以满足实际应用的需求是重中之重。
发明内容
本发明所要解决的技术问题在于解决现有有限元数值模拟计算方法不能快速得到廊管内GIL运行温度，并且不能作为子程序嵌入到GIL温度实时监测系统中的问题，从而提出一种廊管内GIL运行温度评估方法,采用解析方法对GIL的热特性进行求解，不仅能够快速的对GIL热特性问题进行求解，还可作为子程序嵌入到GIL温度监测系统中。对GIL运行温度的快速掌握更有助于提高GIL的热可靠性和更好保障GIL安全稳定运行，该方法还能作为子程序嵌入到GIL温度实时监测系统中。
为解决上述问题，本发明的技术方案为：
一种廊管内GIL运行温度评估方法，包括以下步骤：
步骤S1：建立GIL物理模型，模型涉及电阻焦耳发热效应、外界空气的自然对流层流流动或强迫湍流对流、内部绝缘气体(如SF6气体)自然对流的层流流动、辐射传热等电、热、流、固多物理场，是一个多物理场耦合模型。因此本模型评估廊管内GIL运行温度参数指标包括导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量。
步骤S2：计算导体初始功率损耗，根据上述的物理模型，涉及计算的流体和固体材质属性均考虑随局部温度的变化，即外界空气、内部绝缘气体、导体和外壳材质热物性参数均为温度的函数公式，所述的外壳材质热物性参数包括粘性系数，传导系数，定常比压，比热容、热膨胀系数，由此使得计算时热物性参量的数值更贴近运行实际情况；在实际求解中根据导体的热物性参数和集肤效应系数，确定导体的初始温度，假定初始温度与环境温度相等，根据给定的初始工况参数计算GIL管道启动运行工况下迭代求解出相应公况下导体的生热量；
步骤S3：计算外壳初始功率损耗，计算外壳电阻焦耳发热量时考虑外壳电阻的集肤效应，即外壳截面上，中心电流密度小，靠近表面的电流密度大，根据步骤S1的物理模型，在确定了外壳的热物性参数和集肤效应系数后，确定导体的初始温度，根据给定的初始工况参数计算GIL管道启动工况下迭代求解出相应公况下导体的生热量。
步骤S4：计算外壳与空气的热对流和热辐射量，判断GIL导体整体的产热量与散热量误差是否小于程序设定误差，误差不超过5％，GIL导体整体的产热量由导体的产热量和外壳的产热量组成；GIL导体的散热通过与外界空气直接接触的外壳进行，故GIL导体整体的散热量由外壳与外界空气间的的对流传热量和辐射传热量组成；计算涉及的外壳和外界空气的热物性参数均为温度的函数，通过判断GIL导体整体的产热量和散热量误差是否小于设定误差来判断GIL导体整体的产热和散热是否达到平衡，即GIL导体运行是否已经达到稳态。
步骤S5：若步骤S4误差小于设定误差，则进行下一步骤，否则外壳初始温度增大0.1℃返回步骤S3；
步骤S6：计算导体与外壳间的热对流和热辐射量，判断导体的产热量与散热量误差是否小于设定误差；导体的功率损耗的直接表示形式即为导体的产热量，导体的散热主要通过绝缘气体进行，故导体的散热量由导体与外壳间的自然对流传热量和辐射传热量组成，导体和外壳间的热对流和热辐射量的计算，主要涉及导体和外壳的运行温度，以及绝缘气体的动力粘度、导热系数、定压比热容、容积膨胀系数等热物性参数，其中所有的热物性参数均是温度的函数。通过判断导体的产热量与散热量误差是否小于设定误差来判断导体产热和散热是否达到平衡，即导体运行是否已经达到稳态。
步骤S7：若步骤S6误差小于设定误差，则计算结束输出结果，否则导体温度增大0.1℃返回步骤S2。当导体的产热量与散热量误差小于设定误差时，说明导体的运行达到稳态，可以得到稳定运行时导体的运行温度。若导体的产热量与散热量误差大于设定误差，即导体运行未达到稳态，此时需要增加导体温度再进行迭代计算，直至误差满足条件即可跳出循环，输出GIL运行稳定时接近实际的导体和外壳的温度。
所述步骤S1中GIL物理模型涉及的未知量包括空气与绝缘气体的对流换热系数、热扩散率、热导率、动力粘度、比热容等为温度函数的热物性参数，为温度函数的导体和外壳的导热系数和集肤效应系数，通过采用与温度联系的参数，可以使结果更加贴近实际。
所述步骤S1中初始工况参数包括导体负载电流、廊管内环境温度、GIL管道内压、廊管内空气流速、GIL管道结构参数。
所述步骤S2中，导体存在集肤效应，计算导体电阻焦耳发热量时考虑了导体的集肤效应；确定导体初始温度后，可根据初始工况条件确定导体初始功率损耗，其中相关计算公式为：
式中，Rc为导体的电阻，Kfc为导体的集肤效应系数，ρc20为20℃时导体的直流电阻率，Tc为导体的运行温度，αc20为导体的电阻温度系数，Sc为导体的横截面积，Dc为导体外径，Cc为导体壁厚，Pc为导体的每米损耗，I为导体额定电流。
所述步骤S3中，外壳存在集肤效应，计算外壳电阻焦耳发热量时考虑了外壳的集肤效应；确定外壳初始温度后，可根据初始工况条件确定外壳初始功率损耗。其中相关计算公式为：
式中，Rt为外壳的电阻，Kft为外壳的集肤效应系数，ρta为20℃时外壳的直流电阻率，Tt为外壳的运行温度，αta为外壳的电阻温度系数，St为外壳的横截面积，Dt为外壳外径，Ct为外壳壁厚，Pt为外壳的每米损耗，I为外壳额定电流。
所述步骤S4中，GIL导体整体的产热量由导体的产热量和外壳的产热量组成，散热量由外壳与外界空气间的的对流传热量和辐射传热量组成，相关计算公式为：
式中，Qtc为外壳每米对流散热量，Qtr为外壳每米辐射散热量，αa为空气对流换热系数，Ta为环境空气温度，εt为外壳内表面黑度，Nu为外界空气的努塞尔数，λa为空气的导热系数，Gra为空气的普朗特数，g为重力加速度，βa为空气的容积膨胀系数，μa为空气的动力粘度，Tda为空气定型温度，Δ1为GIL导体的产热量与散热量误差。
所述步骤S6中，导体的功率损耗即为导体的产热量，导体的散热量由导体与外壳间的自然对流传热量和辐射传热量组成，相关计算公式为：
式中，Qcc为导体每米对流散热量，Qcr为导体每米辐射散热量，εc为导体外表面黑度，εe为全组物体的相当黑度，λe为当量导热系数，λ为绝缘气体导热系数，Grs为绝缘气体格拉晓夫数，Prs为绝缘气体普朗特数，Cp为绝缘气体定压比热容，μ为绝缘气体粘度，βs绝缘气体的容积膨胀系数，Rti为外壳内半径，Rco为导体外半径，Tds为绝缘气体的定型温度，Δ2为导体的产热量与散热量误差。
步骤S5中，增大的温度单位值可更具实际情况中GIL测温装置的进度进行设置，计算精度高。
步骤S7中，增大的温度单位值可更具实际情况中GIL测温装置的进度进行设置；输出的结果为GIL运行到稳态时导体、外壳的温度以及需要的绝缘气体的定性温度。
本发明的廊管内GIL运行温度评估方法，基于能量守恒定律建立的以为节点热量传递网络模型，通过在Windows操作系统的MATLAB环境中设计的廊管内特高压GIL热特性MATLAB计算程序，可以利用不同初始工况下的GIL相关特性参数来对其稳定运行时的温度进行评估计算。
本发明通过设计MATLAB程序进行评估计算，为将预测GIL热特性的程序嵌入到GIL温度监测系统中提供了可行的手段。本方法可以方便快捷的求得不同工况下GIL的传热特性以及导体和外壳的温度，进而对GIL从开始运行至稳定时的热特性有一个很好的研究，其实现的计算程序，操作方便，实用性强，计算速度快，可依据不同的工况设置初始参数，为对GIL热特性的研究提供了很好的理论基础，而且可以作为子程序嵌入到GIL的温度监测系统中，对提高实现GIL的可靠性有重要作用。
附图说明
图1是本发明的廊管内GIL运行温度评估方法的流程示意图。
具体实施方式
下面结合幅图对本发明做详细叙述。
本发明主要是想对特高压廊管GIL热特性产热和传热过程通过解析计算进行模拟，首先对进行计算的GIL结构进行简化处理：一是忽略触头部分的功率损耗与传热，二是只对一米长的单相GIL管道进行热特性的研究。GIL母线通电运行后，不仅导体上存在电流，外壳上也有与负载电流等大反向的感应电流，因此导体和外壳均会产生功率损耗。导体和外壳间的热量传递是通过两者间的绝缘气体的对流和辐射进行的，外壳与外界空气的热量传递是通过外界空气的对流和辐射进行的。从能量守恒定律出发，处于热稳定状态时，导体的损耗以辐射方式和自然对流方式传递给外壳，GIL总体损耗以辐射和自然对流方式传递给周围空气，二者处于平衡状态。本模型评估廊管内GIL运行温度的主要参数指标包括导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量。
GIL母线能量守恒：
Pc+Pt＝Qtc+Qtr (5)
Pc＝Qcc+Qcr (6)
其中Pc和Pt分别为导体和外壳的每米功率损耗，W/m；Qtc为外壳与空气的自然对流换热量，W/m；Qtr为外壳外表面的辐射散热量，W/m；Qcc为导体和外壳的对流换热量，W/m；Qcr为导体外表面的辐射散热量，W/m。
主要评估参数的计算：
Pc＝I2Rc (7)
Pt＝I2Rt (8)
Qtc＝αaπDt(Tt-Ta) (11)
其中，Tc为导体的温度，Tt为外壳温度，Ta为初始环境温度。可以从式中得到，导体温度、外壳温度、环境温度与GIL的损耗与传热量息息相关。
参见图1给出的廊管内GIL运行温度评估方法示意图，提供一种廊管内GIL运行温度评估方法，用于评估计算不同初始工况条件下廊管内特高压GIL的稳态运行温度，其评估计算步骤如下：
步骤S1：根据空气、绝缘气体、导体和外壳的相关热物性参数建立GIL物理模型；母线导体和外壳为发热元件,外壳通过自然对流和辐射向外界散发热量。为了反映真实传热过程，导体与外壳之间考虑自然对流和辐射的热量传递方式，即内部绝缘气流的自然对流换热、导体外表面及外壳内表面的热量辐射与吸收，因为假设导体及外壳轴向和周向温度均匀，只考虑导体和外壳在固体域内的径向热传导，由此基于能量守恒定律建立一维节点热量传递网络模型。模型涉及电阻焦耳发热效应、外界空气的自然对流层流流动或强迫湍流对流、内部绝缘气体(如SF6气体)自然对流的层流流动、辐射传热等电、热、流、固多物理场，是一个高度复杂的多物理场耦合模型。因此本模型评估廊管内GIL运行温度参数指标包括导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量。
步骤S2：确定导体初始温度与初始工况条件，计算导体初始功率损耗；与大多数常用电缆计算电阻取为常数不同，计算导体电阻焦耳发热量时考虑了导体电阻的集肤效应(即导体截面上，中心电流密度小，靠近表面的电流密度大)。根据上述的物理模型，涉及计算的流体和固体材质属性均考虑随局部温度的变化，即外界空气、内部绝缘气体、导体和外壳材质热物性参数均为温度的函数公式，所述的外壳材质热物性参数包括粘性系数，传导系数，定常比压，比热容、热膨胀系数，由此使得计算时热物性参量的数值更贴近运行实际情况；在实际求解中根据导体的热物性参数和集肤效应系数，确定导体的初始温度，假定初始温度与环境温度相等，根据给定的初始工况参数计算GIL管道启动运行工况下迭代求解出相应公况下导体的生热量；
步骤S3：确定外壳初始温度与初始工况条件，计算外壳初始功率损耗；计算外壳电阻焦耳发热量时考虑外壳电阻的集肤效应，即外壳截面上，中心电流密度小，靠近表面的电流密度大，根据步骤S1的物理模型，在确定了外壳的热物性参数和集肤效应系数后，确定导体的初始温度(初始温度与环境温度相等)，根据给定的初始工况参数计算GIL管道启动工况下迭代求解出相应公况下导体的生热量。
步骤S4：计算外壳与空气的热对流和热辐射量，判断GIL导体整体的产热量与散热量误差是否小于程序设定误差(一般不超过5％)；GIL导体整体的产热量由导体的产热量和外壳的产热量组成；GIL导体的散热通过与外界空气直接接触的外壳进行，故GIL导体整体的散热量由外壳与外界空气间的的对流传热量和辐射传热量组成；计算涉及的外壳和外界空气的热物性参数均为温度的函数，通过判断GIL导体整体的产热量和散热量误差是否小于设定误差来判断GIL导体整体的产热和散热是否达到平衡，即GIL导体运行是否已经达到稳态。
步骤S5：若步骤S4误差小于设定误差，则进行下一步骤，否则外壳初始温度增大0.1℃返回步骤S3；当GIL导体整体的产热量和散热量误差小于设定误差时，说明GIL导体运行已经达到稳态，可以得到稳定运行时外壳的运行温度；若GIL导体整体的产热量和散热量误差大于设定误差，即GIL导体运行未达到稳态，则需要增加外壳温度再进行迭代计算，直至误差满足条件跳出循环进行下一步。
步骤S6：计算导体与外壳间的热对流和热辐射量，判断导体的产热量与散热量误差是否小于设定误差；步骤S7：若步骤S6误差小于设定误差，则计算结束输出结果，否则导体温度增大0.1℃返回步骤S2。
本发明的廊管内GIL运行温度评估方法，通过在Windows操作系统的MATLAB环境中设计的廊管内特高压GIL热特性MATLAB计算程序，利用初始工况下的GIL相关特性参数来对其稳定运行时的温度进行评估计算。
上述工况初始工况参数包括导体负载电流、廊管内环境温度、GIL管道母线内压、廊管内空气流速、GIL管道结构参数。
给定不同工况的初始参数值，经过计算可得GIL导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量、导体和外壳的温度。
验证示例：
为验证上述方法的正确性，以导体内径为160mm、导体壁厚20mm、外壳内径860mm、外壳壁厚16mm的GIL作为实验对象，针对导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量等数据进行计算。实验GIL的外界环境温度为25℃，实验工况为额定电流8000A。
GIL外界环境温度为25℃，假定GIL初始运行时刻导体的初始温度为24.9℃，外壳的初始温度为25℃，由式7～12便可计算得到初始时刻的导体与外壳的焦耳热损分别为187.96W/m、80.47W/m，导体与外壳间的对流传热量和热辐射量分别为2.19×10-3W/m、-2.48×10-2W/m，外壳与空气间的对流传热量和热辐射量均为0W/m。由计算结果可明显观察到，该时刻GIL并未运行稳定。此时需要增大外壳或者导体的温度进行循环迭代计算，直到设定误差得到满足：
通过计算得到GIL稳定运行时导体与外壳的焦耳热损分别为177.88W/m、41.15W/m，导体与外壳间的对流传热量和热辐射量分别为93.77W/m、75.39W/m，外壳与空气间的对流传热量和热辐射量分别为167.13W/m、15.56W/m。由GIL稳定运行时得到的导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量便可得到稳定运行时导体和外壳的温度分别为66.7℃、39.3℃。而在实际实验中，GIL导体和外壳的温度分别为68.8℃、40.6℃，与本评估方法得到的导体和外壳温度的误差均为3.2％。可见本计算评估方法具有较高的精度，是实际可行的。
本发明的廊管内GIL运行温度评估方法评估计算得到的导体和外壳温度可以表征GIL母线实际正常运行时的导体和外壳温度，在GIL温度监测系统中嵌入本评估计算模型，可以快速发现GIL温度异常现象，并对GIL进行检修，可以有效提高GIL的运行可靠性。
本发明提供的廊管内GIL运行温度评估方法，通过设计MATLAB程序进行评估计算，为将预测GIL热特性的程序嵌入到GIL温度监测系统中提供了可行的手段。本方法可以方便快捷的求得不同工况下GIL的传热特性以及导体和外壳的温度，进而对GIL从开始运行至稳定时的热特性有一个很好的研究，其实现的计算程序，操作方便，实用性强，计算速度快，可依据不同的工况设置初始参数，为对GIL热特性的研究提供了很好的理论基础，而且可以作为子程序嵌入到GIL的温度监测系统中，对提高实现GIL的可靠性有重要作用。
以上对本发明所提供的廊管内GIL运行温度评估方法进行了详细介绍，并应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例子的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，再具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应该理解为对本发明的限制。
Claims (7)
1.一种廊管内GIL运行温度评估方法，其特征在于，包括以下步骤：
步骤S1：建立GIL物理模型，模型涉及电阻焦耳发热效应、外界空气的自然对流层流流动或强迫湍流对流、内部绝缘气体自然对流的层流流动、辐射传热、电、热、流、固多物理场，是一个多物理场耦合模型；因此本模型评估廊管内GIL运行温度参数指标包括导体与外壳的焦耳热损、导体与外壳间的对流传热量和热辐射量、外壳与空气间的对流传热量和热辐射量；
步骤S2：计算导体初始功率损耗，根据上述的物理模型，涉及计算的流体和固体材质属性均考虑随局部温度的变化，即外界空气、内部绝缘气体、导体和外壳材质热物性参数均为温度的函数公式，所述的外壳材质热物性参数包括粘性系数，传导系数，定常比压，比热容、热膨胀系数，由此使得计算时热物性参量的数值更贴近运行实际情况；在实际求解中根据导体的热物性参数和集肤效应系数，确定导体的初始温度，假定初始温度与环境温度相等，根据给定的初始工况参数计算GIL管道启动运行工况下迭代求解出相应工况下导体的生热量；
步骤S3：计算外壳初始功率损耗，计算外壳电阻焦耳发热量时考虑外壳电阻的集肤效应，即外壳截面上，中心电流密度小，靠近表面的电流密度大，根据步骤S1的物理模型，在确定了外壳的热物性参数和集肤效应系数后，确定导体的初始温度，根据给定的初始工况参数计算GIL管道启动工况下迭代求解出相应工况下导体的生热量；
步骤S4：计算外壳与空气的热对流和热辐射量，判断GIL导体整体的产热量与散热量误差是否小于程序设定误差，误差不超过5％，GIL导体整体的产热量由导体的产热量和外壳的产热量组成；GIL导体的散热通过与外界空气直接接触的外壳进行，故GIL导体整体的散热量由外壳与外界空气间的的对流传热量和辐射传热量组成；计算涉及的外壳和外界空气的热物性参数均为温度的函数，通过判断GIL导体整体的产热量和散热量误差是否小于设定误差来判断GIL导体整体的产热和散热是否达到平衡，即GIL导体运行是否已经达到稳态；
步骤S5：若步骤S4误差小于设定误差，则进行下一步骤，否则外壳初始温度增大0.1℃返回步骤S3；
步骤S6：计算导体与外壳间的热对流和热辐射量，判断导体的产热量与散热量误差是否小于设定误差；导体的功率损耗的直接表示形式即为导体的产热量，导体的散热主要通过绝缘气体进行，故导体的散热量由导体与外壳间的自然对流传热量和辐射传热量组成，导体和外壳间的热对流和热辐射量的计算，涉及导体、外壳的运行温度和热物性参数，热物性参数包括绝缘气体的动力粘度、导热系数、定压比热容、容积膨胀系数，其中所有的热物性参数均是温度的函数；通过判断导体的产热量与散热量误差是否小于设定误差来判断导体产热和散热是否达到平衡，即导体运行是否已经达到稳态；
步骤S7：若步骤S6误差小于设定误差，则计算结束输出结果，否则导体温度增大0.1℃返回步骤S2；当导体的产热量与散热量误差小于设定误差时，说明导体的运行达到稳态，可以得到稳定运行时导体的运行温度；若导体的产热量与散热量误差大于设定误差，即导体运行未达到稳态，此时需要增加导体温度再进行迭代计算，直至误差满足条件即可跳出循环，输出GIL运行稳定时接近实际的导体和外壳的温度。
2.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S1中GIL物理模型涉及的未知量为温度函数的热物性参数，包括空气与绝缘气体的对流换热系数、热扩散率、热导率、动力粘度、比热容，为温度函数的导体和外壳的导热系数和集肤效应系数，通过采用与温度联系的参数，可以使结果更加贴近实际。
3.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S1中初始工况参数包括导体负载电流、廊管内环境温度、GIL管道内压、廊管内空气流速、GIL管道结构参数。
4.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S2中，导体存在集肤效应，计算导体电阻焦耳发热量时考虑了导体的集肤效应；确定导体初始温度后，可根据初始工况条件确定导体初始功率损耗，其中相关计算公式为：
式中，Rc为导体的电阻，Kfc为导体的集肤效应系数，ρc20为20℃时导体的直流电阻率，Tc为导体的运行温度，αc20为导体的电阻温度系数，Sc为导体的横截面积，Dc为导体外径，Cc为导体壁厚，Pc为导体的每米损耗，I为导体额定电流。
5.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S3中，外壳存在集肤效应，计算外壳电阻焦耳发热量时考虑了外壳的集肤效应；确定外壳初始温度后，可根据初始工况条件确定外壳初始功率损耗；其中相关计算公式为：
式中，Rt为外壳的电阻，Kft为外壳的集肤效应系数，ρt20为20℃时外壳的直流电阻率，Tt为外壳的运行温度，αt20为20℃时外壳的电阻温度系数，St为外壳的横截面积，Dt为外壳外径，Ct为外壳壁厚，Pt为外壳的每米损耗，I为外壳额定电流。
6.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S4中，GIL导体整体的产热量由导体的产热量和外壳的产热量组成，散热量由外壳与外界空气间的的对流传热量和辐射传热量组成，相关计算公式为：
式中，Qtc为外壳每米对流散热量，Qtr为外壳每米辐射散热量，αa为空气对流换热系数，Ta为环境空气温度，εt为外壳内表面黑度，Nu为外界空气的努塞尔数，λa为空气的导热系数，Gra为空气的普朗特数，g为重力加速度，βa为空气的容积膨胀系数，μa为空气的动力粘度，Tda为空气定型温度，Δ1为GIL导体的产热量与散热量误差。
7.根据权利要求1所述的一种廊管内GIL运行温度评估方法，其特征在于，所述步骤S6中，导体的功率损耗即为导体的产热量，导体的散热量由导体与外壳间的自然对流传热量和辐射传热量组成，相关计算公式为：
式中，Qcc为导体每米对流散热量，Qcr为导体每米辐射散热量，εc为导体外表面黑度，εe为全组物体的相当黑度，λe为当量导热系数，λ为绝缘气体导热系数，Grs为绝缘气体格拉晓夫数，Prs为绝缘气体普朗特数，Cp为绝缘气体定压比热容，μ为绝缘气体粘度，βs绝缘气体的容积膨胀系数，Rti为外壳内半径，Rco为导体外半径，Tds为绝缘气体的定型温度，Δ2为导体的产热量与散热量误差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149403.2A CN108038268B (zh) | 2017-11-18 | 2017-11-18 | 一种廊管内gil运行温度评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149403.2A CN108038268B (zh) | 2017-11-18 | 2017-11-18 | 一种廊管内gil运行温度评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108038268A CN108038268A (zh) | 2018-05-15 |
| CN108038268B true CN108038268B (zh) | 2020-06-02 |
Family
ID=62092757
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711149403.2A Expired - Fee Related CN108038268B (zh) | 2017-11-18 | 2017-11-18 | 一种廊管内gil运行温度评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108038268B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597748B (zh) * | 2020-05-22 | 2024-05-17 | 中国矿业大学 | 一种基于gil热特性实现故障判定的方法 |
| CN112103923B (zh) * | 2020-05-22 | 2021-07-06 | 中国矿业大学 | 一种基于热特性异常实现gil保护的方法 |
| CN113835948B (zh) * | 2020-06-23 | 2024-10-18 | 华为技术有限公司 | 温度检测方法、温度检测装置及电子设备 |
| CN111950797B (zh) * | 2020-08-21 | 2023-03-10 | 中国科学院合肥物质科学研究院 | 一种带连接头的大功率水冷母线局部温度预测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792258A (zh) * | 2015-04-23 | 2015-07-22 | 华北电力大学 | 气体绝缘金属封闭输电线路监测方法及装置 |
| CN104965946A (zh) * | 2015-06-26 | 2015-10-07 | 合肥工业大学 | 一种高铁轴承稳态温度场的有限元计算方法 |
| CN106124928A (zh) * | 2016-06-21 | 2016-11-16 | 中国电力科学研究院 | 特高压交流混合气体gil线段温度监测方法 |
-
2017
- 2017-11-18 CN CN201711149403.2A patent/CN108038268B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792258A (zh) * | 2015-04-23 | 2015-07-22 | 华北电力大学 | 气体绝缘金属封闭输电线路监测方法及装置 |
| CN104965946A (zh) * | 2015-06-26 | 2015-10-07 | 合肥工业大学 | 一种高铁轴承稳态温度场的有限元计算方法 |
| CN106124928A (zh) * | 2016-06-21 | 2016-11-16 | 中国电力科学研究院 | 特高压交流混合气体gil线段温度监测方法 |
Non-Patent Citations (3)
| Title |
|---|
| "Thermal Problems in Gas-Insulated Transmission lines solved by the Finite-Element Method";Xiaowen Wu.etc;《2012 Sixth International Conference on Electromagnetic Field Problems and Applications》;20120924;Page1-4 * |
| "基于热力耦合分析的GIL热致伸缩特性及其影响因素";王健等;《高电压技术》;20170228;第43卷(第2期);第429-437页 * |
| "户内气体绝缘输电线路温升数值计算模型的研究";张改杰等;《陕西电力》;20140331;第42卷(第3期);第45-48页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108038268A (zh) | 2018-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108038268B (zh) | 一种廊管内gil运行温度评估方法 | |
| CN108052697B (zh) | 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 | |
| CN104217061A (zh) | 低压配电柜的温度场仿真设计方法 | |
| Sedaghat et al. | Enhanced thermal model of power cables installed in ducts for ampacity calculations | |
| Zhang et al. | Experiments and modeling of heat transfer in oil transformer winding with zigzag cooling ducts | |
| Rebzani et al. | Numerical multiphysics modeling of temperature rises in gas insulated busbars | |
| Liu et al. | Numerical and experimental investigation of temperature distribution for oil-immersed transformer winding based on dimensionless least-squares and upwind finite element method | |
| CN102880747A (zh) | 一种光电复合海缆温度场建模计算分析方法 | |
| CN110333443B (zh) | 感应电机定子绕组温升测试方法 | |
| Zhao et al. | Real-time reconstruction of temperature field for cable joints based on inverse analysis | |
| Lin et al. | Theoretical research on temperature field of power cable joint with FEM | |
| Dao et al. | A study of hot-spot localization in distribution transformers | |
| Boukrouche et al. | Mock-up study of the effect of wall distance on the thermal rating of power cables in ventilated tunnels | |
| Fu et al. | Parameter identification for current–temperature relationship of contact wire under natural convection | |
| CN115061039B (zh) | 一种高压断路器内接线端子特性评估方法及相关设备 | |
| Hua et al. | Lifetime Prediction of Aviation Cables Based on an Improved Electromagnetic-Thermal Coupling Model | |
| Ryzhov et al. | Three-dimensional thermal stator model of a fully air-cooled turbogenerator | |
| Ai et al. | Coupled Electromagnetic and Thermal Analysis of a 15kW Cryogenic Induction Motor for Submerged Liquefied Natural Gas Pumps | |
| Ediriweera et al. | Localizing hotspot in an oil immersed distribution transformer using finite element analysis | |
| Luo et al. | Investigation of contact resistance influence on power cable joint temperature based on 3-D coupling model | |
| Lewandowska et al. | Experimental study of steady-state transverse heat transfer in a single channel CICC | |
| Xu et al. | Calculation of the equivalent thermal parameters of external environment for cable based on the real-time measurement data | |
| Chen et al. | A Fast Calculation Model of Buried Cable Ratings based on Distributed Temperature Sensors | |
| Park et al. | A study for the estimation of temperature and thermal life of traction motor for commercial HEV | |
| Liang et al. | An optimal heat line simulation method to calculate the steady-stage temperature and ampacity of buried cables |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190212 Address after: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant after: XI'AN JIAOTONG University Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. MAINTENANCE BRANCH Address before: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant before: Xi'an Jiaotong University |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200602 |