基于导线径向温度和蠕变的导线弧垂计算方法.pdf

面向双碳目标的环保型电力设备DOI:10.15961/j.jsuese.202200951基于导线径向温度和蠕变的导线弧垂计算方法赵隆1，王如水1，代鹏飞1，朱超2，黄新波1(1.西安工程大学 电子信息学院，陕西 西安 710600；2.陕西电力科学研究院，陕西 西安 710100)摘要：输电线路弧垂是线路安全运行的重要指标，针对输电导线温度非均匀分布和蠕变带来的弧垂计算偏差，本文开展了径向温度和蠕变对弧垂计算模型影响程度的研究。首先，基于导线的物理构成和结构特征建立了导线的有限元模型，对不同运行状态下的输电导线进行有限元仿真，分析了不同运行状态下的导线径向温度分布规律；其次，建立了考虑径向温差的导线应力、导线径向热膨胀和蠕变的弧垂计算模型，基于力平衡原理研究了考虑径向温差的导线径向热膨胀、导线应力与伸长量的关系，分析了蠕变与伸长量的关系，并通过伸长量与弧垂的关系函数得到导线弧垂；最后，为了验证该弧垂计算模型的准确性和实用性，本文基于陕西某110 kV实际输电线路进行了验证。结果表明：对于本现场实际运行的导线来说，蠕变对弧垂计算的准确性影响最大，其次是导线径向温差，最后是径向热膨胀；在忽略导线径向温差、径向热膨胀和蠕变时，弧垂计算的误差范围为45.2%30.0%，最大误差高达45.2%，而本文的优化模型计算误差小于1.1%；当忽略导线径向温差和径向热膨胀时，导线弧垂计算存在误差，且随着径向温差增大，忽略导线径向温差和径向热膨胀带来的弧垂计算误差也逐渐增大。关键词：输电导线；径向温度；蠕变；弧垂；有限元仿真中图分类号：TM751文献标志码：A文章编号：2096-3246（2023）04-0021-09Calculation Method of Conductor Sag Based on Radial Temperature and CreepZHAO Long1，WANG Rushui1，DAI Pengfei1，ZHU Chao2，HUANG Xinbo1(1.School of Electronics Information,Xian Polytechnic Univ.,Xian 710600,China;2.State Grid Shaanxi Electric Power Research Inst.,Xian 710100,China)Abstract:The sag of the transmission line is an important index for its safe operation.Focusing on the calculation deviation of conductor sagcaused by the non-uniform temperature distribution of transmission lines and the creep,the influence degree of the radial temperature and thecreep on the sag calculation model was studied in this paper.Firstly,the finite element model of the conductor was established based on its phys-ical composition and structural characteristics.Considering different operating conditions of the transmission line,the radial temperature was ana-lyzed through the finite element simulation.Secondly,considering the radial temperature difference of stress,the radial thermal expansion and thecreep,the sag calculation model was built.Based on the force balance principle,the relationship among the radial thermal expansion consideringthe radial temperature difference,the stress and the strain of the conductor was studied and the relationship between creep and strain was ana-lyzed.Then the conductor sag was calculated through the function of the strain and the sag.Finally,in order to verify the accuracy and practicabil-ity of the sag calculation model,the model was verified by a 110 kV transmission line in Shaanxi Province.The results show that the influence onthe accuracy of the sag calculation from largest to smallest are the creep,the radial temperature difference and the radial thermal expansion basedon the 110 kV transmission line in Shaanxi Province.Without taking into account the radial temperature difference,the radial thermal expansionand creep,the error of sag was in the range from 45.2%to 30.0%.And the maximum error was 45.2%.The error range under the optimization收稿日期:2022 09 04基金项目:陕西省重点研发计划（2023-YBGY-069）；国网陕西电科院委托项目（2126KY19006K）；陕西省重点研发计划（2021GY-068）；陕西省重点研发计划（2018ZDXM-GY-040）作者简介:赵隆（1987），男，副教授，博士.研究方向：输电线路在线监测与状态评估技术；结构健康监测技术；无线传感器.E-mail：网络出版时间:2023 06 28 13:19:07 网络出版地址:https:/ http:/http:/ 第 55 卷 第 4 期工 程 科 学 与 技 术Vol.55 No.42023 年 7 月ADVANCED ENGINEERING SCIENCESJuly 2023model of sag in this paper was less than 1.1%.When the influence of the radial temperature difference and the radial thermal expansion on the sagis ignored,the error of sag is existed and the error becomes greater with the increase of the radial temperature difference.Key words:transmission line;radial temperature;creep;sag;finite element simulation 近年来，随着新能源大规模并网及社会发展对电能需求的快速增长，导线运行安全不断受到新的挑战。其中，输电线路弧垂是线路安全运行的重要指标。受温度变化，导线易出现弧垂增大现象，严重时可能产生对地放电12，威胁输电线路安全。导线弧垂计算通常考虑导线温度、档距、长度、高差、应力和比载等参数3。上述因素中，除导线温度和应力外，其余均为固定参数，而应力也受到导线温度的影响。因此，导线温度变化是运行过程中影响导线弧垂变化的主要因素。目前，针对输电线路的弧垂计算已经有许多方法。其中，Nigol4和Barrett5等提出了一种基于ACSR导体的应力应变和温度应力数据曲线的弧垂算法，但该方法忽略了导线径向温度分布不均匀的特点；Alawar等6提出了一种基于导线应力应变曲线的导线弧垂预测方法，但该方法没有考虑导线温度的影响；Albizu等7基于高温低弧垂导线分别分析了基于导线应力应变曲线的弧垂算法和考虑应变、温度、蠕变等多个因素的弧垂算法；Dong等8提出了结合悬链线函数、应力应变关系以及3种应变之间的关系和张力平衡方程的弧垂算法，但该方法仍是将导线看作统一的整体，并未考虑导线径向温度分布不均匀的特点。然而，由于导线内外层的材料不同，不同层导线的物理参数（如弹性模量、热膨胀系数等）也存在明显差异9。此外，刘刚10和李军辉11等通过有限元方法发现导线存在径向温度梯度；冯凯等12基于导线径向热路模型分析发现导线存在径向温度梯度。因此，在进行弧垂计算时不能把导线简单地当作统一整体进行分析。Liu等13提出了基于导线温度场分析的弧垂算法，即通过每一层导线的温度分析、计算导线的分层应力、弧垂和临界温度，但该方法忽略了同一层导线的不同根线股也存在较大温差，没有考虑蠕变特性对弧垂计算的影响。Albizu等14结合导线载流量、温度等运行参数，针对含空气隙的导线提出了一种应变张力的弧垂算法。Guo等15提出了基于导线径向温度分析的弧垂算法，即通过有限元仿真分析导线径向温度计算了导线弧垂。上述两种算法忽略了导线径向热膨胀的影响。由于钢芯铝绞线与其他金属一样具备蠕变特性和热膨胀性16，忽略导线蠕变和导向径向热膨胀也为导线弧垂计算带来了一定误差，因此上述算法仍需要改进。在线监测法主要为利用张力、倾角和温度传感器对导线弧垂进行实时监测17以及利用GPS定位或北斗差分定位实现对导线弧垂的精准测量18，虽然这些方法对导线弧垂进行了实时监测，但无法实现对弧垂变化的预测，不能实现故障发生前的预警功能。麻卫峰等19建立了基于航拍图像的弧垂监测方法，该方法需要通过航拍对导线弧垂进行测量和分析1921，无法满足实时回传数据、及时发现故障的需求；宰红斌等22提出了基于机器学习的导线弧垂估计方法；李嘉雨等23提出了基于PSO-BP算法的导线弧垂预测方法；刘沛轩等24提出了基于遗传算法的导线弧垂修正算法，这些方法需要大量的图片累计来完成模拟计算，才能达到合适的准确度。因此，为了更准确地计算出导线弧垂，本文提出了一种基于导线径向温度和蠕变的导线弧垂计算方法。基于输电导线的有限元仿真，考虑导线径向膨胀对计算每一根线股对输电导线弧垂的影响。同时计算导线蠕变，实现对导线弧垂算法的修正，最终结合现场监测数据验证算法的准确性。结合气象参数实现对导线弧垂的预测，更好地实现动态增容技术。1 导线径向温度分布规律由于导线内外层材料和温度均不相同，为了更准确地分析导线弧垂，本文分析了不同载流量和不同环境条件下的导线径向温度分布规律。以导线LGJ-300/40为例，建立了流固耦合3维导线的有限元模型。该导线是由7根半径为1.33 mm的钢芯和24根半径为1.98 mm的铝股绞合而成，且相邻层绞线反向绞制，导线模型长度为200 mm，图1为导线物理模型。(a)LGJ300/40 导线截面图(b)LGJ300/40 导线侧视图图 1LGJ300/40导线物理模型Fig.1Three-dimensional physical model of LGJ300/40 22工程科学与技术第 55 卷导线的钢芯钢芯之间、钢芯铝股之间以及铝股铝股之间存在空气间隙。在钢芯铝绞线的制作过程当中，相邻导体相互摩擦、挤压使其表面凹凸不平，两两导体之间存在点接触，导线相关参数见表1。在实际运行过程中，导线可视作无限场，流体域为无限大。此外，风速和风向实时发生变化；随着太阳的角度不同，导线受到日照强度的影响不同。考虑到影响导线温度的因素主要为风速、日照强度以及环境温度，对模型作出如下简化：1）由于导线外径较小，当导线热交换范围不超过流体域时，就认为导线处于无限大空间内；2）当导线的长度包含一个节距时，对导线设置symmetry面就默认导线为无限长；3）默认某时刻的风速、风向以及日照强度是固定的。导线热交换方式为接触热传递、对流散热和辐射散热。其中，接触热传递主要存在于线股之间；对流散热主要存在于线股与空气之间；辐射散热则存在于每一层线股。因此，本文将3维导线模型导入有限元仿真软件中，进行预处理和导线命名。为了更准确地模拟导线实际情况，将流体域和导线的延伸面设置为symmetry面，并设置风速出入口。其次，按照数量最少、仿真结果最准确的原则进行了网格划分。对于流体域，网格划分尺寸相对较大。而导线线股存在扭曲，且导线接触使得导线之间的空气间隙狭小并存在尖角，因此网格划分需要单独细化和加密，网格划分尺寸最小，如图2所示。其中，蓝色箭头表示风速入口，红色箭头表示风速出口。右上角方框中展示了紧密网格的局部放大结果，不同颜色的圆球代表钢芯铝绞线的线股。由于导线内外层材料不同，其电阻也不同，钢芯和铝股产生的热通量有所差异。因此，在计算导线线股热生成率时，根据导线流过的电流和电磁感应对铝股及钢芯的热生成率进行单独计算。此外，由于导线本身的物理结构特征，导线的最外层铝股会受到太阳辐射和风速的影响，故根据现场运行数据导线最外层铝股需要单独施加太阳辐射热量，根据风速出入口方向设置风速大小和角度，再设置导线所处的环境温度。本文以陕西省某110 kV输电线路的实时监测结果为例，表2给出了该线路在2020年8月2日至8月8日期间每日12：30监测所得气象、导线载流量和弧垂等数据。其中，8月2日12：30的风速为1.1 m/s，风向角为151，日照强度为850 W/m2，载流量为132.8 A，导线径向温度分布如图3所示。由图3可以看出，导线温度存在明显的径向梯度，且每一股导线的温度各不相同，最大径向温差约为2.9，主要受到风速、风向角、日照强度、导线载流量的影响。由于空气间隙、风速和风向角的存在，导线迎风侧温度最低，从迎风侧到背风侧温度逐渐升高；导线上下表面温度变化较小；钢芯散热情况最差，但由于固体传热比气体传热快，所以钢芯温度不一定最高。日照产生的能量只能辐射到导线表面，从而影响导线的径向温度分布。载流量流过导线产生的 YXZ图 23维导线模型网格Fig.2Grids of 3D conductor model 风向1.1 m/s温度/C35.6535.1034.5634.0133.4732.9232.3831.8331.2930.7430.20图 3导线径向温度分布Fig.3Distribution of radial temperature on conductor 表 1LGJ300/40导线参数Tab.1 Conductor parameters of LGJ300/40 导线数量/根半径/mm弹性模量/GPa绞入率/%最外层次外层内层铝芯241.98553.372.881.11镀锌钢芯71.33190 表 2气象、导线载流量以及弧垂数据Tab.2 Data of weather,ampacity of conductors and sag 日期载流量/A风速/(ms1)风向/（）环境温度/日照强度/(Wm2)弧垂/m8月2日132.81.115130.28502.338月3日156.51.13032.32232.098月4日136.42.015630.33642.048月5日163.51.317133.23752.158月6日151.92.43028.81371.978月7日84.11.52824.41251.888月8日137.51.716430.65022.08 第 4 期赵隆，等：基于导线径向温度和蠕变的导线弧垂计算方法23热量因导线材料而异，且空气隙传热性能不如金属接触好，故存在径向温差。因此，在计算导线弧垂时，需要考虑每一股导线的温度。2 考虑径向温度和蠕变的弧垂计算模型传统弧垂计算模型中并未考虑导线轴径向温度和蠕变的影响，因此本节将采用悬链线模型分析径向温度分布和蠕变对输电导线弧垂及应力的影响，优化导线弧垂计算模型。随着外部环境、导线架设条件和导线相关参数的变化，弧垂也随之变化。本文将导线弧垂的变化分为两个部分，包括永久性形变和非永久性形变。永久性形变主要由残余形变（绞合线相互挤压）和蠕变伸长构成；非永久性形变主要由弹性形变和导线热膨胀形变构成25。对导线考虑径向温度的弧垂计算作以下假设：1）导线的铝股和钢芯具有弹性线性应力应变行为，且导线只承受拉力，不考虑弯曲、剪切和扭转的力；2）考虑横向形变时，由于导线的径向热膨胀远大于泊松效应和相邻层挤压力产生的形变，因此忽略泊松效应和相邻层挤压力产生的形变；3）导线受到的内外力平衡，忽略残余形变。2.1 计及径向温度的应力与伸长量fL0导线应力是影响导线弧垂的重要因素之一。在导线实际运行过程中，导线弧垂 未知，导线长度 未知，故通过基础的几何关系得到导线的单位面积初始应力，图4为输电导线示意图。lhh图4中，A和B分别为导线的左右悬挂点，C为最低点。输电导线的跨距为，高差为（当右侧悬挂点B高于左侧悬挂点A时，为正，否则为负），为A、B两悬挂点连线与水平方向的夹角，为A点弧垂切线与x轴的夹角。0初始应力的计算公式13为：0lcosh(l20)l+(hLh=0)2+h2+fAC=0（1）0Lh=0fAC式中，为单位面积初始应力，为导线比载，为高差为0时的导线长度，为导线最低点的坐标值。Lh=0fACLh=0fAC在计算过程中，假定架空线比载按线向均匀分布，根据这种原则推导得到的架空线悬垂形状是悬链线型的。由式（1）可知，和均为未知量，因此根据悬链线方程分别可以求得和，其计算公式分别为：Lh=0=20sinhl20（2）fAC=0(1cosha0)（3）根据悬链线方程，未知量a的计算公式为：a=l20arcsinhhLh=0（4）0N(x)nn在一段档距内，导线单位面积下对外表现的应力是一个定值，因此导线每一处的轴向外力也是一个定值。本文将大档距的导线划分为 个小段，求解 个点的导线轴向外力，其计算公式25为：N(x)=0A01+tan(l2x)20cos 2,x=0,ln,2ln,l（5）A0式中，为导线的总截面积。在同一个点上，导线轴向外力和导线内力相等，导线内力又受到导线径向分布的影响，其计算公式13为：P(0)=mi=1piasin aia+kj=1pjssin ajs（6）mkpiapjsijaiaajs式中，和 分别代表铝股个数和钢芯个数，和分别为第 根铝股和第 根钢芯所受内力，和分别为铝股和钢芯与横截面的夹角。受到环境参数的影响，导线径向温度分布不均匀，故不同层导线以及同一层的每根导线受到的轴向内力也不相同。因此要对每根导线进行应力分析。铝股的内力计算公式为：pia=Eia(iaiaT)(Aia+(iaTRia)2)（7）EiaiaiiaTiAiaRia式中，为铝股的弹性模量，为第 根铝股的应变，为在温度作用下第 根铝股的应变，为铝股的截面积，为铝股的半径。由于钢芯和铝股的弹性模量、线股温度以及半径均不相同，故内力计算公式不同，钢芯的内力计算公式为：pjs=Ejs(jsjsT)(Ajs+(jsTRjs)2)（8）Ejsjsj式中，为钢芯的弹性模量，为第 根钢芯的应变，ABxyChfACl/2l图 4输电导线示意图Fig.4Diagrammatic sketch of the transmission line 24工程科学与技术第 55 卷jsTjAjsRjs为在温度作用下第 根钢芯的应变，为钢芯的截面积，为钢芯的半径。由于钢芯和铝股的线膨胀系数和线股温度不相同，其在温度作用下的应变也不相同，铝股在温度作用下的应变计算公式为：iaT=a(TiaTa)（9）aTiaiTa式中，为铝股线膨胀系数，为第 根铝股的温度，为环境温度。钢芯在温度作用下的应变计算公式为：jsT=s(TjsTa)（10）sTjsj式中，为钢芯线膨胀系数，为第 根钢芯的温度。0N(x)根据式（1）（4）可以得到导线的单位面积初始应力，将其代入式（5），得到导线轴向外力；将式（7）（10）代入式（6），根据导线轴向内外力相等，可以得到由每股导线伸长量。2.2 计及径向温度的径向热膨胀及蠕变与伸长量导线的横向形变的影响因素主要为3部分，分别为泊松效应、相邻层的挤压力和热膨胀，前两者占比极低，因此横向形变仅考虑导线的径向热膨胀。其中，受到径向温度的影响，每股导线的径向热膨胀量不同，其造成每股导线伸长量也不相同。以一股导线为例，将其展开得到每股导线伸长量和导线总伸长量的几何关系。一股导线的展开图如图5所示。0aiaaiaRiaRia图5中，为导线由应力和热膨胀产生的总伸长量，为铝股受到张力产生的夹角增量，为每股导线与横截面的夹角，为铝股的半径，为铝股由于导线径向热膨胀产生的半径增量。当存在径向温差时，每股导线温度不相同，导致每股导线的伸长量不相同。导线伸长量是由应力和热膨胀产生的总伸长量构成的函数，其计算公式13为：ia=0sin2(aia+aia)RiaRiacos2(aia+aia)js=0sin2(ajs+ajs)RjsRjscos2(ajs+ajs)（11）aiaajsaiaajsaiaajsRjsRjs式中：和分别为铝股和钢芯受到张力产生的夹角增量，但对于任意一股导线来说，每股导线与横截面的夹角或远大于由于张力产生的夹角增量或，因此夹角增量可忽略不计；为钢芯的半径；为钢芯由于导线径向热膨胀产生的半径增量。根据式（11）可分别得到每根铝股和钢芯产生的伸长量与导线总伸长量的关系。最后将得到的关系公式代入式（6），由于导线轴向内外力相等，可得到导线由应力和热膨胀产生的总伸长量。c蠕变是指在外部作用力下的长期形变，它是由材料的蠕变特性引起的，该部分的导线形变是塑性蠕变延伸，属于永久性形变。为了简化导线弧垂的计算，根据GB/T 365512018同心绞架空导线性能计算方法25获得10 a后钢芯铝绞线的典型蠕变值=500 mm/km。该蠕变值是由许多蠕变实验分析验证后得到的。2.3 导线弧垂计算基于径向温度和蠕变对导线伸长量的影响，得到了修正后的导线总伸长量变化量，其计算公式为：总=0+c00=0n（12）总00式中，为导线总伸长量变化量，为由引起的导线伸长量。由于初始弧长未知，单位面积下初始应力不准确，因此本文采用修正后导线弧长计算导线弧垂，导线最大弧垂fM的计算公式为：fM=(Llcos)3l8cos3（13）对于实际运行的导线，首先根据单跨架空线的几何特点得到式（1），由悬链线方程得到式（2）（4），并将式（2）（4）代入式（1），根据牛顿法求解得到单位面积下的初始应力，再将其代入式（5）得到每一段导线的轴向外力。其次，考虑了导线径向温度分布规律，即导线每股受力均不相同，因此将式（9）和（10）分别代入式（7）和（8），得到每股导线受力。由于导线径向温度分布不均匀，还考虑了径向热膨胀影响，最终将式（7）、（8）、（11）代入式（6），得到导线总轴向内力，根据导线内外力平衡得到导线由应力和径向热膨胀引起的总伸长量。最后，将典型蠕变值代入式（12），得到导线总伸长量变化量，最后再将式（12）代入式（13），得到修正后弧垂。其计算原理如图6所示。为了验证作者提出的考虑导线径向温度和蠕变的弧垂计算模型具有应用价值，本节通过改变档距 导线伸长导线展开(Ria+Ria)aiaaia+aiaiaSi0LSLSSi图 5一股导线展开示意图Fig.5Unfolding view of one wire 第 4 期赵隆，等：基于导线径向温度和蠕变的导线弧垂计算方法25分析了导线蠕变对导线弧垂的影响，改变导线径向温差分析了导线径向热膨胀和径向温度对导线弧垂的影响。其中，不改变其他条件，导线档距从50 m增加至400 m时，蠕变对弧垂的影响如图7所示。从图7中可以看出，当档距小于100 m时，蠕变对弧垂的影响较小，导线弧垂变化小于0.12 m，此时，为简化计算可以忽略蠕变。但当档距大于100 m时，蠕变对弧垂的影响较大，此时，不可以忽略蠕变。随着档距的增大，蠕变对弧垂的影响也越来越大，当档距为400 m时，弧垂的变化高达0.78 m。因此，随着输电导线电压等级的提高和输送距离的延长，不可以忽略蠕变对导线弧垂的影响。不改变其他条件，导线径向温差从0 增大至10 时，导线径向膨胀对弧垂变化的影响如图8所示。从图8中可以看出，如果忽略导线径向温差和径向热膨胀，导线弧垂计算存在误差，且随着径向温差的增大，忽略导线径向温差和径向热膨胀带来的弧垂计算误差也逐渐增大。随着电能需求的增加，导线径向温差增大，高电压等级、大档距线路数量也增加，故考虑径向温度分布和蠕变的弧垂计算模型具有实际工程应用意义。3 现场数据验证及结果分析为了验证径向温度分布规律、径向热膨胀以及导线蠕变对弧垂变化的影响，本文利用表2中2020年8月2日至8月8日的气象、导线载流量以及弧垂数据，验证该修正算法的有效性。由于导线径向温度分布受到环境温度、风速、风向以及日照强度的影响，故根据在线监测装置得到的数据模拟导线在该时刻的径向温度，再根据导线径向温度分布特征得到该时刻导线弧垂计算值，并将其与弧垂监测值做对比，其计算结果见表3。从表3中看出，改进后的弧垂算法计算误差小于1.1%。忽略径向温度和蠕变的弧垂计算值误差范围在45.2%30.0%，最大误差高达45.2%（即此时 根据图 4 计算导线单位面积下初始应力将导线划分为 n 段，计算每一段的轴向外力计算导线所受内力已知 i=0计算由应力和热膨胀引起的导线伸长量ini=i+1是否导线内外力平衡根据仿真结果得到每一根线股芯温度计算每一根铝股所受内力计算每一根钢芯所受内力计算每一根铝股径向热膨胀计算铝股和钢芯伸长量与 0 的关系 f(0)=N(x)P(0)计算由蠕变引起的导线伸长量_计算平均伸长量 0_0计算 fM总=0+c0_0=0A0EiaAia+EjsAjs 图 6弧垂计算原理Fig.6Schematic diagram of sag calculation 3.02.52.0弧垂/m1.51.00.5050100 150 200档距/mx=1.10 x=1.22y=100考虑蠕变忽略蠕变250 300 350 400 450图 7蠕变对弧垂变化的影响Fig.7Influence of creep on sag 3.0弧垂/m2.52.002468径向温差/C忽略径向膨胀考虑径向膨胀与径向温差忽略径向温差10图 8导线径向膨胀对弧垂变化的影响Fig.8Influence of radial thermal expansion on sag 26工程科学与技术第 55 卷弧垂算法得到的弧垂计算值比实际弧垂小45.2%），误差极大。针对该段输电线路，考虑了径向温度和蠕变的弧垂计算值比忽略径向温度和蠕变的弧垂计算值更准确，修正后的算法可以更准确地得到导线弧垂且符合工程应用要求。其误差对比如图9所示。从图9中可以看出，本文中提出的弧垂计算模型误差范围约为0.010.02 m，具有可实施性。此外，当忽略导线径向温度和蠕变后，误差范围约为0.650.88 m，是本文计算模型误差的4488倍。这是由于影响导线弧垂的主要因素为导线径向温度和蠕变，其中导线径向温度又包括导线径向温差和导线径向热膨胀。对输电导线实施动态增容技术时，由于忽略了径向温度和蠕变，估算的弧垂值较小，极易引发安全距离不足等故障。因此，为保证输电线路的安全运行，径向温度和蠕变条件不可忽略。本文改进的弧垂算法中分别考虑了导线蠕变、径向膨胀和径向温差对弧垂的影响，不同条件下的弧垂计算值见表4。由于影响导线弧垂计算的因素分别为：导线径向温差、导线径向热膨胀以及导线蠕变。而图9仅综合分析了径向温度和蠕变对导线弧垂计算的影响，并未逐一分析上述因素对弧垂变化的影响。故下文分别研究了导线蠕变、导线径向膨胀和导线径向温差对导线弧垂的影响，图10为不同条件下的弧垂计算值对比结果。从图10中可以看出，对于该段输电线路来说，导线蠕变、导线径向膨胀和导线径向温差对弧垂计算 2.42.22.0弧垂/m1.81.41.68 月 2 日8 月 3 日8 月 4 日8 月 5 日8 月 6 日8 月 7 日8 月 8 日1.01.2测量值忽略蠕变的计算值忽略径向膨胀的计算值忽略径向温差的计算值日期图 10不同条件下的弧垂计算值对比Fig.10Comparison of calculated sag under different con-ditions 3.02.52.0弧垂/m1.51.08 月 2 日8 月 3 日8 月 4 日8 月 5 日8 月 6 日8 月 7 日8 月 8 日测量值考虑径向温度和蠕变的计算值忽略其他条件的计算值00.5日期图 9弧垂测量值和弧垂计算值误差对比Fig.9Comparison of errors between measured and calcu-lated sag 表 3导线弧垂计算值及误差Tab.3 Calculated sag and its error 日期测量值/m考虑径向温度和蠕变计算值/m考虑径向温度和蠕变的误差/%忽略径向温度和蠕变的计算值/m忽略径向温度和蠕变的误差/%8月2日2.332.310.861.6330.08月3日2.092.100.481.3137.38月4日2.042.060.981.2936.88月5日2.152.170.931.4333.58月6日1.971.951.021.0944.78月7日1.881.870.531.0345.28月8日2.082.090.481.4331.3 表 4不同条件的弧垂计算值Tab.4 Sag calculation under different conditions 日期考虑径向温差与蠕变/m忽略径向温差/m忽略蠕变/m忽略径向膨胀/m8月2日2.311.752.272.228月3日2.101.462.082.028月4日2.061.432.032.008月5日2.171.592.132.108月6日1.951.201.911.878月7日1.871.111.861.848月8日2.091.592.072.02 第 4 期赵隆，等：基于导线径向温度和蠕变的导线弧垂计算方法27的影响不同。其中，由于气象参数和导线载流量不同，导线的径向温差也不同，随着导线径向温差的变化，导线弧垂的变化不同。其中，蠕变对弧垂计算的准确性影响最大，其次是导线径向温差，最后是导线径向热膨胀。忽略导线蠕变、导线径向膨胀和导线径向温差对导线弧垂的影响，导线弧垂计算值均偏小，可能造成事故。4 结论本文提出了一种考虑导线径向温度梯度、导线径向热膨胀和蠕变的弧垂计算模型，改进了传统导线弧垂模型，并通过现场实测数据验证了该计算模型的准确性，分析了蠕变、导线径向热膨胀对弧垂的影响，得出了以下结论：1）采用有限元仿真得到导线径向温度梯度，并基于导线径向温度梯度优化了弧垂计算模型，改进的弧垂算法误差小于1.1%，符合工程应用的要求。2）如果忽略径向温度和蠕变，导线的弧垂计算值误差范围在45.2%30.0%，与实际值偏差较大。3）忽略导线蠕变、导线径向膨胀和导线径向温差的影响，计算得到的弧垂远小于实际值，会极大地提高输电导线的运行风险。参考文献：Li Bin,Wang Guoli,Liu Lei,et al.Full scale test study on sagcharacteristics of augmented capacity conductors for trans-mission lineJ.Southern Power System Technology,2017,11(6):3541.李斌,王国利,刘磊,等.输电线路增容导线弧垂特性真型试验研究J.南方电网技术,2017,11(6):3541.1 Li Mingming,Wang Jian,Xiong Xiaofu,et al.Operating tem-perature and sag off-limit warning method for overheadtransmission lines in the conditions of high-temperatureweatherJ.Power System Protection and Control,2020,48(2):2533.李明明,王建,熊小伏,等.高温天气下架空线路运行温度与弧垂越限预警方法J.电力系统保护与控制,2020,48(2):2533.2 Liu Chen,Li Yongshuang,Tang Jian.Stress calculationmethod of overhead lines with multiple concentrated loadsand continuous gearsJ.High Voltage Engineering,2021,47(12):42884297.刘晨,李永双,唐剑.多集中荷载连续档架空线应力计算方法J.高电压技术,2021,47(12):42884297.3 Nigol O,Barrett J S.Characteristics of ACSR conductors athigh temperatures and stressesJ.IEEE Transactions onPower Apparatus and Systems,1981,PAS-100(2):485493.4 Barrett J S,Dutta S,Nigol O.A new computer model of AC-SR conductorsJ.IEEE Transactions on Power Apparatus5and Systems,1983,PAS-102(3):614621.Alawar A,Bosze E J,Nutt S R.A hybrid numerical methodto calculate the sag of composite conductorsJ.ElectricPower Systems Research,2006,76(5):389394.6 Albizu I,Mazn A J,Valverde V,et al.Aspects to take intoaccount in the application of mechanical calculation to high-temperature low-sag conductorsJ.IET Generation,Trans-mission&Distribution,2010,4(5):631640.7 Dong Xiaoyuan.Analytic method to calculate and charac-terize the sag and tension of overhead linesJ.IEEE Trans-actions on Power Delivery,2016,31(5):20642071.8 Nie Song,Cheng Yangchun,Dai Yuan.A new sag calcula-tion method based on the temperatur