一种电缆导体温度间接测量方法
【专利摘要】一种电缆导体温度间接测量方法，在电缆外表皮选取测量点E，电缆通流后用温度传感器采集E点的温度值，得到测量点温度测量值集Tm(t)；建立电缆温度场有限元模型，对电缆温度场进行分析计算，得到E点温度计算值TE；通过对E点温度实测值Tm(t)与计算值TE之间的误差进行处理，对温度场模型中的初始热参数进行修正；用修正后的热参数代替原先的初始热参数载入到有限元热模型中，从而获取表面温度和导体温度的对应关系表，测量到表面温度后利用查表方式得到导体温度。本发明通过对电缆内部材料的热参数进行辨识从而校正电缆的温度场模型，在校正后的温度场模型基础上得到正确的表面温度与导体温度的对应关系，从而通过对表面温度的测量间接得出导体温度。
【专利说明】
-种电缆导体溫度间接测量方法
技术领域
[0001] 本发明设及传输电缆技术领域，尤其是一种电缆导体溫度间接测量方法。
【背景技术】
[0002] 电力电缆是城市输配电网络的核屯、载体，近些年来使用数量不断增加。电缆运行 中由于电流流过导体线忍，引起线忍焦耳损耗、主绝缘介质损耗W及金属护套损耗等，造成 电缆导体溫度升高。导体溫度过高将加速材料老化、缩短电缆使用寿命，甚至是发生击穿事 故。因此，从电力电缆自身的安全运行角度出发，需要对电力电缆的导体溫度进行实时监 测。
[0003] 通过传感器直接测量电缆导体溫度在电缆运行中难W实现，目前主要采用一些间 接测量电缆导体溫度的方法，主要是解析计算方法和数值计算方法。
[0004] 解析法WIEC标准为代表，根据电缆的结构详细分析每个介质层的结构和热学特 征，逆向分析热流传导过程，推算导体溫度。运种方法往往基于W下假设：（1)大地表面为等 溫面；（2)电缆表面为等溫面；（3)叠加原理适用。但实际中电缆敷设环境复杂，往往无法满 足W上的计算条件。
[0005] 数值计算方法W有限元法为代表，通过对电缆溫度场进行有限元计算，获取不同 载荷、边界条件下的表面溫度和导体溫度的对应关系表，测量到表面溫度后利用查表方式 得到导体溫度。该方法的不足之处是电缆溫度场建模中需要标定电缆内部不同材料的热参 数值，而厂家往往不能提供运些参数值。即使可W提供，运些值在电缆的长期运行中也会由 于老化等原因发生很大的变化。
【发明内容】
[0006] 本发明所要解决的技术问题是提供一种电缆导体溫度间接测量方法，通过对电缆 内部材料的热参数进行辨识从而校正电缆的溫度场模型，在校正后的溫度场模型基础上得 到正确的表面溫度与导体溫度的对应关系，从而通过对表面溫度的测量间接得出导体溫 度。
[0007] 本发明所采用的技术方案是：
[000引一种电缆导体溫度间接测量方法，包括W下步骤：
[0009] 步骤1、在电缆外表皮选取测量点E，电缆通流后用溫度传感器采集E点的溫度值， 得到测量点溫度测量值集Tm(t);
[0010] 步骤2、根据电缆的结构参数、初始热参数、实测时电缆的载荷值和所处环境建立 电缆溫度场有限元模型，对电缆溫度场进行分析计算，得到E点溫度计算值TE;
[0011] 步骤3、通过对E点溫度实测值Tm(t)与计算值TE之间的误差进行处理，对溫度场模 型中的初始热参数进行修正；
[0012] 步骤4、用修正后的热参数代替原先的初始热参数载入到有限元热模型中，在此模 型的基础上进行溫度场计算，从而获取不同载荷、边界条件下的表面溫度和导体溫度的对 应关系表，测量到表面溫度后利用查表方式得到导体溫度。
[0013] 进一步的，本发明方法中，所述步骤1中的溫度测量是在电缆表面选定点敷设热电 偶传感器，传感器通过连接导线接入测控仪，测控仪可控制传感器溫度采集的周期，为保证 实测数据的数量可W有效用于后续步骤的开展，需保证一定的溫度测量时长。
[0014] 进一步的，本发明方法中，步骤2中电缆溫度场计算包括W下步骤：
[0015] (2.1)、根据结构参数G对电缆进行几何建模；
[0016] (2.2)、指定电缆内部不同材料的初始热参数X0,包括材料热传导率、密度、比热 容；
[0017] (2.3)、对几何模型进行网格划分；
[0018] (2.4)、指定热载荷（负荷电流I)、初始溫度To和边界条件B(边界对流换热系数和 环境溫度）；
[0019] (2.5)、求解计算电缆的溫度场并提取E点的溫度计算值。
[0020] 进一步的，本发明方法中，步骤3中对溫度场模型中初始热参数进行修正的过程 为：
[0021 ] (3.1 )、基于计算值与实测值的目标函数的构建；
[0022] 电缆的瞬态热特性取决于W下5个参数：1)负荷电流I; 2)初始溫度To; 3)电缆表面 的边界条件B;4)电缆结构参数G;5)电缆材料的热参数X。所W，E的溫度计算值TE的函数表 达式可表达为TE = f(I，T〇，B，G，x)，式子f表明TE是输入变量组I，To，B，G，X非线性映射的输 出结果。
[0023] 上述5个参数中，热参数X在电缆长期运行中受材料老化等因素影响变化较大，无 法准确获取;负荷电流I、初始溫度To可由相应设备进行测量;边界条件B可由经验公式计算 得到;结构参数G厂家会有手册说明。结合分析，TE函数表达式的输入变量组中除热参数X之 外另外4个参数T〇，B，G，I都可视为常数，于是可将E的溫度计算值TE视为热参数X的映射，记 为了6 = 了￡(1，义）。式中义即为电缆内部材料的热参数，1为时间。
[0024] 由W上分析可看出，热参数X是影响溫度场计算结果的重要因素，而往往溫度场有 限元计算中选用的初始热参数/只是经验数值，可能与实际值偏差较大，必须进行修正。本 发明结合E点的计算溫度值与实测溫度值Tm(t)来构造目标函数F，通过对F进行最小化处理 来实现初始热参数/中相关元素的修正，其中目标函数
巧中，X为所 有热参数集合，X = (XI，X2，…，xn);N是实测值总数;测量时间间隔，= TD/N，其中TD为 测量时间。
[0025] (3.2)、W所述目标函数取得全局最小值为停止条件，通过最速下降法进行迭代求 解，得到热参数的修正值X*，具体步骤如下：
[0026] (3.2.1) W初始热参数/作为初点，设置允许误差ε = 0.1，置k = 1;
[0027] (3.2.2)计算捜索方向(1(^，其中：；
式中，目 标函数F相对X中的第j个参数^的偏微分可由式(1)求得：
[0031 ]式中：xe= (XI ,X2,... ,Xj+e ,,Χη) ;xd = (XI ,X2,... ,χ广ε , ,χη) ;Te( j A t ,xe)和Te( j A t，xd)的值通过将xe与xd代入有限元热模型中进行溫度场计算得到；
[0032] (3.2.3)若||(1^||《6，停止计算，输出义《作为修正值义*;否则，从义〇〇出发，沿4? 进行一维捜索，求λιο使得：
[0033]
(3)
[0034] 式中:Ak是从xW出发沿方向dW进行一维捜索的步长.
[0035] (3.2.4):令x(k+i) = x(k)+Akd(k)，置k = k+l，转步骤(3.2.2)。
[0036] 进一步的，本发明方法中，步骤4中构建表面溫度与导体溫度的映射表格是指将修 正值代替原先的初始热参数载入到有限元热模型中，在此模型的基础上进行溫度场计算并 提取出每个时刻的E点溫度值与导体溫度值，由其构建表面溫度与导体溫度的映射表格。
[0037] 本发明一种电缆导体溫度间接测量方法，有益效果是：
[0038] (1)针对利用解析法间接获取导体溫度时假设条件要求难W满足的缺陷，本方法 基于有限元法实现导体溫度的测量；（2)针对利用有限元法对导体溫度进行间接测量时由 于不能正确获取热参数从而影响测量结果准确性运一缺陷，本发明设计了一种间接测量方 法。此方法首先对电缆表面选定点的溫度进行实测，然后通过对E点溫度实测值Tm(t)与计 算值TE之间的误差进行处理，对有限元溫度场模型中的热参数进行了修正从而使计算值最 接近于实测值。将修正后的热参数载入溫度场有限元模型中，在此模型的基础上进行溫度 场计算，由计算结果构建表面溫度与导体溫度的映射表格，测量表面溫度后即可通过查询 表格的方式知道此时的导体溫度。
【附图说明】
[0039] 下面结合附图和实施例对本发明作进一步说明：
[0040] 图1为本发明实例中的电缆截面图。
[0041 ]图2为本发明的测量方法流程图。
[0042] 图3(a)为本发明实施例的电缆几何模型图。
[0043] 图3(b)为本发明实施例的有限元模型图。
[0044] 图4为本发明实施例的电缆材料组成图。
[0045] 图5(a)为本发明实施例的电缆溫度化时的溫度分布云图。
[0046] 图5(b)为本发明实施例的电缆溫度化时的溫度分布云图。
[0047] 图6为本发明实施例的E点溫度实测值与计算值的曲线图。
[0048] 图7为本发明实施例的热参数修正流程图。
[0049 ]图8为本发明实施例的800A工况下表面溫度实测值与计算值的比较。
[0050]图9为本发明实施例的表面溫度与导体溫度的对应曲线图。
【具体实施方式】
[0051] -种电缆导体溫度间接测量方法，包括W下步骤：
[0052] 步骤1、在电缆外表皮选取测量点E，电缆通流后用溫度传感器采集E点的溫度值， 得到测量点溫度测量值集Tm(t);
[0053] 步骤2、根据电缆的结构参数、初始热参数、实测时电缆的载荷值和所处环境建立 电缆溫度场有限元模型，对电缆溫度场进行分析计算，得到E点溫度计算值TE;
[0054] 步骤3、通过对E点溫度实测值Tm(t)与计算值TE之间的误差进行处理，对溫度场模 型中的初始热参数进行修正；
[0055] 步骤4、用修正后的热参数代替原先的初始热参数载入到有限元热模型中，在此模 型的基础上进行溫度场计算，从而获取不同载荷、边界条件下的表面溫度和导体溫度的对 应关系表，测量到表面溫度后利用查表方式得到导体溫度。
[0056] 进一步的，本发明方法中，所述步骤1中的溫度测量是在电缆表面选定点敷设热电 偶传感器，传感器通过连接导线接入测控仪，测控仪可控制传感器溫度采集的周期，为保证 实测数据的数量可W有效用于后续步骤的开展，需保证一定的溫度测量时长。
[0057] 进一步的，本发明方法中，步骤2中电缆溫度场计算包括W下步骤：
[0058] (2.1)、根据结构参数G对电缆进行几何建模；
[0059] (2.2)、指定电缆内部不同材料的初始热参数X0,包括材料热传导率、密度、比热 容；
[0060] (2.3)、对几何模型进行网格划分；
[0061] (2.4)、指定热载荷（负荷电流I)、初始溫度To和边界条件B(边界对流换热系数和 环境溫度）；
[0062] (2.5)、求解计算电缆的溫度场并提取E点的溫度计算值。
[0063] 进一步的，本发明方法中，步骤3中对溫度场模型中初始热参数进行修正的过程 为：
[0064] (3.1 )、基于计算值与实测值的目标函数的构建；
[0065] 电缆的瞬态热特性取决于W下5个参数：1)负荷电流I; 2)初始溫度To; 3)电缆表面 的边界条件B;4)电缆结构参数G;5)电缆材料的热参数X。所W，E的溫度计算值TE的函数表 达式可表达为TE = f(I，T〇，B，G，x)，式子f表明TE是输入变量组I，To，B，G，X非线性映射的输 出结果。
[0066] 上述5个参数中，热参数X在电缆长期运行中受材料老化等因素影响变化较大，无 法准确获取;负荷电流I、初始溫度To可由相应设备进行测量;边界条件B可由经验公式计算 得到;结构参数G厂家会有手册说明。结合分析，TE函数表达式的输入变量组中除热参数X之 外另外4个参数T〇，B，G，I都可视为常数，于是可将E的溫度计算值TE视为热参数X的映射，记 为了6 = 了￡(1，义）。式中义即为电缆内部材料的热参数，1为时间。
[0067] 由W上分析可看出，热参数X是影响溫度场计算结果的重要因素，而往往溫度场有 限元计算中选用的初始热参数/只是经验数值，可能与实际值偏差较大，必须进行修正。本 发明结合E点的计算溫度值与实测溫度值Tm(t)来构造目标函数F，通过对F进行最小化处理 来实现初始热参数^*^中相关元素的修正，其中目标函邀
式中，X为所 有热参数集合，x=(xl，x2，…，xn);N是实测值总数;^t现慢时间间隔，^t = TD/N，其中TD为 测量时间。
[0068] (3.2 )、W所述目标函数取得全局最小值为停止条件，通过最速下降法进行迭代求 解，得到热参数的修正值具体步骤如下：
[0069] (3.2.1) W初始热参数/作为初点，设置允许误差ε = 0.1，置k= 1;
[0070] (3.2.2)计算捜索方向dW，其中：<
式中，目 标函数F相对X中的第j个参数^的偏微分可由式(1)求得：
[0074] 式中：xe= (xi,X2,... ,Xj+e ,,Χη) ;xd = (xi,X2,... ,Xj-e , ,Χη) ;Te( j A t ,xe)和 Te( j A t，xd)的值通过将xe与xd代入有限元热模型中进行溫度场计算得到；
[0075] (3.2.3)若||(1?||《6，停止计算，输出义(^作为修正值义*;否则，从义〇〇出发，沿4? 进行一维捜索，求λιο使得：
[0076]
(3)
[0077] 式中:Ak是从xW出发沿方向dW进行一维捜索的步长.
[007引（3.2.4):令x(k+i) = x(k)+Akd(k)，置k = k+l，转步骤(3.2.2)。
[0079] 进一步的，本发明方法中，步骤4中构建表面溫度与导体溫度的映射表格是指将修 正值代替原先的初始热参数载入到有限元热模型中，在此模型的基础上进行溫度场计算并 提取出每个时刻的E点溫度值与导体溫度值，由其构建表面溫度与导体溫度的映射表格。
[0080] 下面结合附图给出本发明较佳实施例，W详细说明本发明的技术方案。
[0081] 图1为本发明中需要测量导体溫度的一种llOkV交联聚乙締电缆的结构示意图。如 图所示:需要测量导体溫度的一种1 lOkV交联聚乙締电缆包括铜导体1，导体1表面有导体屏 蔽2，导体屏蔽2的外层为化PE绝缘3，化PE绝缘3的外层设有绝缘屏蔽4，绝缘屏蔽4的外层依 次设有皱纹侣护套5和外护套6。
[0082] 图2为本发明的方法流程图。实施本发明的步骤如下：
[0083] 步骤1、在电缆外表皮选取测量点E，电缆通流后用溫度传感器采集E点的溫度值， 得到测量点溫度测量值集Tm(t)。本实施例中电缆的负荷电流为400A、环境溫度为23°C。测 溫模块使用K型热电偶作为溫度传感器，可实现精度±2°C、范围0-200°C的电缆外表面溫度 测量。将热电偶传感器贴敷在测量点E处，启动测溫模块，便可进行溫度值的实时连续采集， 采集到的溫度数据通过USB接口的方式传输到后台的测溫软件中予W显示保存。另外，本实 施例中设置采集周期Imin，采集时长化。
[0084] 步骤2、建立电缆溫度场有限元模型，对电缆溫度场进行分析计算，得到E点溫度计 算值TE。具体过程如下列子步骤：
[0085] (2.1)、根据结构参数G对电缆进行几何建模。本实施例中电缆的结构参数如表1所 示，建立的几何模型如图3a)所示；
[0086] (2.2)、指定电缆内部不同材料的初始热参数X0,包括材料热传导率、密度、比热 容。本实施例中电缆的材料类型如图4所示，可W看出其中有5种材料:铜(导体）、侣(皱纹侣 护套）、交联聚乙締(绝缘层）、半导体材料(导体屏蔽层与绝缘屏蔽层）、橡胶与聚乙締的混 合材料(外护套）。对应运5种材料，共需要15种热参数:λ。。。(导体热导率）、狂。《(导体密度）、 Ccnn(导体比热）；λΑ1(侣护套热导率）、:勉（侣护套密度）、CA1(侣护套比热）；、ns(绝缘层热导 率）、扔?(绝缘层密度）、Cins(绝缘层比热）；λ^3Ε(外护套热导率(外护套密度）、Cja。(外 护套比热）；Ascr(屏蔽层热导率）、批α (屏蔽层密度）、Csu(屏蔽层比热）。综上， Jf_ 心又"Μ、恥〇)，、Cc。"·'华姑、C ΑΓ，文'.ins、巧ins、C i"s. Aji化'、φjuc、Cβκ.' /iscr、.巧sc,'、C …〇
[0087] 在对电缆进行溫度场计算中，导体、侣护套、屏蔽层的热参数在电缆运行中变化较 小，可通过规格书得到，如表1所示。至于其它的六个热参数(、ns、fW、Cins ; Ajac、妈。。、Cjac)， 其在运行中受老化等因素影响较大因此不能准确获得，可先根据经验数据设定它们的初始 值，在通过后续步骤对它们进行修正，设定初始值如表2所示。表1与表2的数据构成了热参 数初始值
[0088] 表1已知电缆热参数
[0092] (2.3)、对几何模型进行网格划分得到有限元模型。本实施例中有限元模型如图 3b)所示；
[0093] (2.4)、指定热载荷（负荷电流I)、初始溫度To和边界条件B(边界对流换热系数和 环境溫度）。本实施例中按照实测溫度时的条件给有限元模型加载，即负荷电流400A、环境 溫度23°C、初始溫度To与环境溫度相同也为23°C、边界对流换热系数由式(4)计算为7.371。 其中，模型中的热源实际上是电缆中各部分的损耗，由负荷电流结合经验公式可计算得到 导体损耗化= 48.351W/m;绝缘损耗化= 0.14抓/m;金属护套损耗化= 0.614W/m，将它们载入 到模型中再进行下一步计算。
[0094] (2.5)、求解计算电缆的溫度场并提取E点的溫度计算值。为与实测溫度进行对照， 本实例中的溫度场计算时长也为化，溫度场分布云图如图5所示。其中图5a)是计算进行到 化时电缆的溫度分布云图，图5b)是计算进行到化时电缆的溫度分布云图。图中标注的E点 即溫度监测点，提取计算时间内E点的溫度值，并W时间（单位:h)为横轴、溫度值(单位：°C) 为纵轴绘制E点溫度值随时间的变化曲线，即图6中的表面溫度初始计算值曲线所示。另外， 图6中其它的曲线在下续步骤中描述。
[0095] 步骤3、通过对E点溫度实测值Tm(t)与计算值TE之间的误差进行处理，对溫度场模 型中的初始热参数进行修正。为与步骤2中得到的计算值进行比较，本实例中将步骤1中得 到的E点溫度实测值也绘制曲线示于图6中，由图可看出在初始热参数/的条件下计算值与 实测值有一定的偏差，需要对初始热参数进行修正。
[0096] 对热参数的修正流程如图7所示。本发明WE点溫度实测值Tm(t)与计算值TE为基 础构造目标函数，通过选择最优化方法对目标函数F进行最小化来对热参数进行修正，将溫 度场有限元计算程序作为最优化方法的子程序，从而将有限元程序与主程序优化算法相衔 接，最后输出修正值本发明中采用最速下降法作为优化方法，具体步骤如下：
[0097] (3.1 )、基于计算值与实测值的目标函数的构建；
[0098] 电缆的瞬态热特性取决于W下5个参数：1)负荷电流I; 2)初始溫度To; 3)电缆表面 的边界条件B;4)电缆结构参数G;5)电缆材料的热参数X。所W，E的溫度计算值TE的函数表 达式可表达为TE = f(I，T〇，B，G，x)，式子f表明TE是输入变量组I，To，B，G，X非线性映射的输 出结果。
[0099] 上述5个参数中，热参数X在电缆长期运行中受材料老化等因素影响变化较大，无 法准确获取;负荷电流I、初始溫度To可由相应设备进行测量;边界条件B可由经验公式计算 得到;结构参数G厂家会有手册说明。结合分析，TE函数表达式的输入变量组中除热参数X之 外另外4个参数T〇，B，G，I都可视为常数，于是可将E的溫度计算值TE视为热参数X的映射，记 为TE = TE(t，x)。式中X即为热参数，t为时间。
[0100] 由W上分析可看出，热参数X是影响溫度场计算结果的重要因素，而往往溫度场有 限元计算中选用的初始热参数/只是经验数值，可能与实际值偏差较大，必须进行修正。本 发明结合E点的计算溫度值与实测溫度值Tm(t)来构造目标函数FW拟合出热参数集合X中 所有元素的真实值，目标函鑽
式中，X为所有热参数集合，x=(xi， X2，…，Xn); N是实测值总数;现慢时间间隔，= TD/N，其中TD为测量时间。
[0101] (3.2)、W所述目标函数取得全局最小值为停止条件，通过最速下降法进行迭代求 解，得到热参数的修正值具体步骤如下：
[0102] (3.2.1)^初始热参数/作为初点，设置允许误差6 = 〇.1，置4=1;
[0103] (3.2.2)计算捜索方向dW，其中
式中，目 标函数F相对X中的第j个参数^的偏微分可由式(1)求得：
[0104]
U)
[0105] 其中TE(jAt，x)相对xj的偏微分由式(2)求得：
[0106]
(2)
[0107]式中：xe= (xi,X2,... ,Xj+e ,,Χη) ;xd = (xi,X2,... ,χ广 ε , ,χη) ;Te( j A t ,xe)和 Te( j A t，xd)的值通过将xe与xd代入有限元热模型中进行溫度场计算得到；
[010引（3.2.3)若||(1?||《6，停止计算，输出义(^作为修正值义*;否则，从义《出发，沿(1? 进行一维捜索，求λιο使得：
[0109]
C3J
[0110] 式中:Ak是从出发沿方向dW进行一维捜索的步长；
[0111] (3.2.4):令义化")=义化)+人1^(1化)，置1^ = 4+1，转步骤(3.2.2)。
[0112] 为对修正后的热参数勺有效性进行考证，本实例另选取负载电流800A、环境溫 度29°C的工况，将优化参数代入有限元程序进行计算并仍提取E点的溫度值绘制曲线，如图 8所示。实测值与计算值的最大溫差为1. rC，误差在可接受范围内，证明了优化后参数的有 效性。
[0113] 步骤4、用修正后的热参数代替原先的初始热参数载入到有限元热模型中，在此模 型的基础上进行溫度场计算，从而获取不同载荷、边界条件下的表面溫度和导体溫度的对 应曲线，测量到表面溫度后在曲线上进行捜寻即可得到导体溫度。本实例中针对环境溫度 23°C、电缆负载电流400A工况下的电缆，在对其热参数进行修正后得到正确的有限元热模 型，在此模型基础上进行溫度场计算得到不同时间点的溫度场分布，然后提取各时间点电 缆表面E点的溫度计算值与内部导体上某点溫度的计算值，在此基础上W表面E点的溫度计 算值为横坐标、导体上某点的计算值为纵坐标绘制对应曲线，如图9所示。通过实测得到表 面E点溫度后利用其在曲线上捜寻对应点，此点对应的纵轴值即导体溫度。
[0114] W上所述仅是本发明的优选实施方式，但运些说明不能被理解为限制本发明的范 围，本发明的保护范围由随附的权利要求书限定。任何在本发明权利要求基础上的改动也 应视在本发明的保护范围内。
【主权项】
1. 一种电缆导体温度间接测量方法，其特征在于包括以下步骤： 步骤1、在电缆外表皮选取测量点E，电缆通流后用温度传感器采集E点的温度值，得到 测量点温度测量值集Tm(t); 步骤2、根据电缆的结构参数、初始热参数、实测时电缆的载荷值和所处环境建立电缆 温度场有限元模型，对电缆温度场进行分析计算，得到E点温度计算值TE; 步骤3、通过对E点温度实测值Tm(t)与计算值TE之间的误差进行处理，对温度场模型中 的初始热参数进行修正； 步骤4、用修正后的热参数代替原先的初始热参数载入到有限元热模型中，在此模型的 基础上进行温度场计算，从而获取不同载荷、边界条件下的表面温度和导体温度的对应关 系表，测量到表面温度后利用查表方式得到导体温度。2. 根据权利要求1所述的一种电缆导体温度间接测量方法，其特征在于，步骤1中的温 度测量是在电缆表面选定点敷设热电偶传感器，传感器通过连接导线接入测控仪，测控仪 可控制传感器温度采集的周期，为保证实测数据的数量可以有效用于后续步骤的开展，需 保证一定的温度测量时长。3. 根据权利要求1所述的一种电缆导体温度间接测量方法，其特征在于，步骤2中电缆 温度场计算包括以下步骤： (2.1 )、根据结构参数G对电缆进行几何建模； (2.2)、指定电缆内部不同材料的初始热参数/，包括材料热传导率、密度、比热容； (2.3 )、对几何模型进行网格划分； (2.4) 、指定热载荷(负荷电流I)、初始温度To和边界条件B(边界对流换热系数和环境温 度)； (2.5) 、求解计算电缆的温度场并提取E点的温度计算值。4. 根据权利要求1所述的一种电缆导体温度间接测量方法，其特征在于，步骤3中对温 度场模型中的初始热参数进行修正的过程为： (3.1 )、基于计算值与实测值的目标函数的构建； 电缆的瞬态热特性取决于以下5个参数：1)负荷电流I; 2)初始温度To; 3)电缆表面的边 界条件B;4)电缆结构参数G;5)电缆材料的热参数X。所以，E的温度计算值TE的函数表达式 可表达为TE = f (I，To，B，G，X)，式子f表明TE是输入变量组I，To，B，G，X非线性映射的输出结 果； 上述5个参数中，热参数X在电缆长期运行中受材料老化等因素影响变化较大，无法准 确获取；负荷电流I、初始温度To可由相应设备进行测量;边界条件B可由经验公式计算得 到;结构参数G厂家会有手册说明；结合分析，TE函数表达式的输入变量组中除热参数X之外 另外4个参数T()，B，G，I都可视为常数，则E的温度计算值TE可视为热参数 X的映射，记为TE = 1^(七0)。式中1即为热参数4为时间； 由以上分析可看出，热参数X是影响温度场计算结果的重要因素，而往往温度场有限元 计算中选用的初始热参数/只是经验数值，可能与实际值偏差较大，必须进行修正；结合E 点的实测温度值!^(〇来构造目标函数F以拟合出热参数集合X中所有元素的真实值，目标 函数：。式中，X为所有热参数集合，1=(11，1 2，一，^1)4是实测值总数； 测量时间间隔，= TD/N，其中TD为测量时间； (3.2)、以所述目标函数取得全局最小值为停止条件，通过最速下降法进行迭代求解， 得到热参数的修正值，，具体步骤如下： (3.2.1) 以初始热参数X*3作为初点，设置允许误差ε =0.1，置k= 1; (3.2.2) 计算搜索方向d(k\其中：.式中，目标函 数F相对X中的第j个参数^的偏微分可由式(1)求得：(1) 其中Te (j Δ t，X)相对Xj的偏微分由式(2)求得： (2) j式中：xe = (XI，X2，…，χ」+ε，…，xn) ;xd = (XI，X2，…，Xj-ε，，xn);Te( j Δ t，xe)和Te( j Δ t， xd)的值由将xe与xd代入有限元热模型中进行温度场计算得到； (3.2.3) 若||(1(1{)||〇，停止计算，输出，)作为修正值，否则，从，)出发，沿(1 (1{)进行 一维搜索，求Ak，使得： F (X ⑷ + /L?i ?/")) = /(X ⑷ + /Li/⑷） （.4.). 式中Ak是从x(k)出发沿方向d(k)进行一维搜索的步长； (3.2.4) :令叉(1^1)=叉(1〇+人1^(1〇，置1^ = 1^+1，转步骤(3.2.2)。5.根据权利要求1所述的一种电缆导体温度间接测量方法，其特征在于，步骤4中构建 表面温度与导体温度的映射表格是指将修正后的热参数代替原先的初始热参数载入到有 限元热模型中，通过不同载荷、边界条件下的计算可获得完整的表面温度和导体温度的对 应关系表，测量到表面温度后利用查表方式即可间接得到导体温度。
【文档编号】G01K7/02GK106066212SQ201610362296
【公开日】2016年11月2日
【申请日】2016年5月27日 公开号201610362296.0, CN 106066212 A, CN 106066212A, CN 201610362296, CN-A-106066212, CN106066212 A, CN106066212A, CN201610362296, CN201610362296.0
【发明人】张宇娇, 黄雄峰, 吴馨, 徐彬昭
【申请人】三峡大学