专利名称：基于环境因子监测和有限元的地下电力电缆载流量在线预测系统及其方法
技术领域：
本发明涉及一种基于环境因子监测和有限元的地下电力电缆载流量在线预测系 统及其方法。
背景技术：
地下高压电力电缆是城市输配电线路的主要形式，由于其安装成本较高、不易维 修的特点，在保证其长期安全可靠运行的基础上，合理确定高压电力电缆的载流量，最大限 度地发挥其输送电能的能力对于电力部门具有重要的意义。目前地下电力电缆载流量确定的方法有两种，一种是离线计算，即根据当地的土 壤性质、最热月份的最高气温，根据IEC60287或有限元进行计算；另一种是在线计算，通过 在线监测电缆表面的温度，通过热路或有限元的方法推算电缆线芯温度，并在此基础上估 算电缆的载流量。电缆实时载流量是与电缆结构、周围土壤的导热系数、地表空气等密切相关的物 理量，上述两种方法都不能准确地在线预测电缆的载流量。第一种方法是根据最高气温和 给定土壤导热系数的情况计算而得的，当气温或土壤导热系数发生变化时，电缆可承受的 最大负荷电流，即载流量也应进行相应调整；对于第二种方法，通过电缆表面的温度可以 较准确的计算出当前负荷下的电缆线芯温度，而要预测最大可承受电流即载流量必须考虑 气温、土壤导热系数的影响。目前，地下电力电缆载流量的计算方法主要以有限元计算为主，首先假定一个电 流，根据边界条件计算电力电缆绝缘层的最高温度，如果不超过绝缘层长期工作耐受温度， 则增加电流，如果超过耐受温度，则减少电流，循环迭代直到电缆绝缘层最高温度等于绝缘 层长期工作耐受温度为止，此时的电流即为电缆的载流量。
发明内容
本发明所要解决的技术问题是提供一种基于环境因子监测和有限元的地下电力 电缆载流量在线预测系统及其方法。本发明解决其技术问题所采用的技术方案一、技术方案1 本发明由用于测量空气温度的温度传感器、风速传感器、湿度传感器、单片机系 统、电力部门的调度主机组成；所述温度传感器、风速传感器、湿度传感器的输出端分别接 所述单片机系统的相应输入端，所述单片机系统通过GSM网络与所述调度主机相连接；所 述温度传感器和风速传感器安装在所述电力电缆的正上方，所述湿度传感器安装在所述电 力电缆沿线的土壤中。二、技术方案2:本发明方法的具体步骤如下
(1)建模a. 土壤直埋式电缆电力电缆距地面lm，电力电缆两侧边界距电缆10m，计算中，所述两侧边界法向方 向热流变化率为0，即土壤温度不再变化；电力电缆下面至土壤深层边界距电缆10m，取土 壤深层边界温度不变，且温度值为土壤深层温度即383K ；b.排管式电缆排管顶部距地面lm，电力电缆敷设于排管内，排管两侧边界距排管10m，计算中， 所述两侧边界法向方向热流变化率为0，即土壤温度不再变化；排管下面至土壤深层边界 距排管10m，取土壤深层边界温度不变，且温度值为土壤深层温度即383K ；C.沟槽式电缆沟槽深距地面lm，电力电缆敷设于沟槽内，室内沟槽顶部常敷设铁板，室外沟槽顶 部常敷设水泥板，模型两侧边界距沟槽10m，计算中，所述两侧边界法向方向热流变化率为 0，即土壤温度不再变化；沟槽下面至土壤深层边界距沟槽10m，取土壤深层边界温度不变， 且温度值为土壤深层温度即383K ；(2)剖分将电力电缆及其周围土壤剖分为小的单元；在直埋式、排管式和沟槽式三种方式下，采用三角形或四边形单元将整个模型剖 分为小的单元；电缆区域温度变化较大，剖分密度较高，排管内电缆外表面与排管内壁间空 气域和沟槽内电缆外表面与沟槽内壁间空气域，温度变化较大，且需计算流速，剖分密度较 高，土壤区域温度变化较小，剖分密度较小；(3)读取实时测量的下述数据空气温度T&、空气风速v&、土壤湿度rh ；并计算 地表空气对流换热系数α和土壤的导热系数λ ；a.利用下述公式(1)计算地表空气对流换热系数α α = 7. 371+6. 43vair0·75 (1)式中，α为地表空气对流换热系数，单位W/(πι2·Κ) ；为风速传感器实时测量的 地表风速，单位m/s ；b.利用公式(2)计算土壤导热系数入A = (l 1.33 IogrA-3.4)10° 64Vrf"2 (2)式中，λ为土壤导热系数，单位W/(m2· K) ;Pd为土壤密实度，单位kg/m3;rh为湿 度传感器实时测量的土壤相对湿度。(4)首先，选取初始电流11，单位A ；(5)调用有限元温度场计算程序，计算电流Il下的电力电缆绝缘层最高温度Tl， 单位K;(6)当Tl小于电力电缆绝缘层长期工作绝缘层耐受温度Tin时，且Tin-Tl > 0. 3K，则取 12 = 211，转至第(7)步；当 Tl 大于 Tin 时，且 Tl-Tin >0. 3K，则取 12 = 0. 511， 转至第(7)步；当ITl-TinI彡0.5K时，Il即为所求电力电缆载流量，退出；(7)调用有限元温度场计算程序，计算电流12下的电力电缆绝缘层最高温度T2 ；(8)当T2-Tin >0. 3K时，转至第(9)步；当T2_Tin彡0· 3K时，12即为所求 电力电缆载流量，退出；(9)取I = 12+(12-11) · (Tin_T2) / (T2-T1)，调用有限元温度场计算程序，计算电流I下的电力电缆绝缘层最高温度T ；(10)当 T-Tin > 0. 3K 时，取 Il = 12，12 = I,Tl = T2，T2 = Τ，转至上述第(9) 步；当IT-TinI彡0.3Κ时，I即为所求电力电缆载流量，退出。本发明的有益效果是能够较准确的在线预测地下电力电缆的载流量，电力部门据 此能够最大限度发挥其输送电能的能力，并且能够在保证安全可靠运行的前提下，不会出 现大马拉小车的资源浪费现象。
图1为本发明的硬件配置图。图2-1为土壤直埋式地下电力电缆敷设方式的示意图。图2-2为排管式地下电力电缆敷设方式的示意图。图2-3为沟槽式地下电力电缆敷设方式的示意图。图3-1为直埋式电缆及周围土壤的有限元整场剖分结果示意图。图3-2为直埋式电缆及周围土壤的有限元电缆部分的剖分结果示意图。在图1、2、3中，1温度传感器(型号DS18B20)，2风速传感器(型号:CQ2-FC_1)，3 湿度传感器(型号BQ8-CHR-01)，4单片机系统(单片机的型号89C52)，5GSM网络，6电力 部门的调度机(台式计算机)，7 土壤，8电力电缆线，9排管，10盖板，11支撑板，12沟槽。
具体实施例方式实施例1 (本系统的实施例，参见图1、2)本实施例由用于测量空气温度的温度传感器1、风速传感器2、湿度传感器3、单片 机系统4、电力部门的调度主机6组成；所述温度传感器1、风速传感器2、湿度传感器3的 输出端分别接所述单片机系统的相应输入端，所述单片机系统通过GSM网络5与所述调度 主机6相连接；所述温度传感器1和风速传感器2安装在所述电力电缆的正上方，所述湿度 传感器3安装在所述电力电缆沿线的土壤中。实施例2 (本方法的实施例，参见图1、2、3)本实施例的具体步骤如下(1)建模a. 土壤直埋式电缆电力电缆距地面lm，电力电缆两侧边界距电缆10m，计算中，所述两侧边界法向方 向热流变化率为0，即土壤温度不再变化；电力电缆下面至土壤深层边界距电缆10m，取土 壤深层边界温度不变，且温度值为土壤深层温度即383K ；b.排管式电缆排管顶部距地面lm，电力电缆敷设于排管内，排管两侧边界距排管10m，计算中， 所述两侧边界法向方向热流变化率为0，即土壤温度不再变化；排管下面至土壤深层边界 距排管10m，取土壤深层边界温度不变，且温度值为土壤深层温度即383K ；C.沟槽式电缆沟槽深距地面lm，电力电缆敷设于沟槽内，室内沟槽顶部常敷设铁板，室外沟槽顶 部常敷设水泥板，模型两侧边界距沟槽10m，计算中，所述两侧边界法向方向热流变化率为
60，即土壤温度不再变化；沟槽下面至土壤深层边界距沟槽10m，取土壤深层边界温度不变， 且温度值为土壤深层温度即383K ；(2)剖分将电力电缆及其周围土壤剖分为小的单元；在直埋式、排管式和沟槽式三种方式下，采用三角形或四边形单元将整个模型剖 分为小的单元；电缆区域温度变化较大，剖分密度较高，排管内电缆外表面与排管内壁间空 气域和沟槽内电缆外表面与沟槽内壁间空气域，温度变化较大，且需计算流速，剖分密度较 高，土壤区域温度变化较小，剖分密度较小；(3)读取实时测量的下述数据空气温度Tai,、空气风速Vai,、土壤湿度rh ；并计算 地表空气对流换热系数α和土壤的导热系数λ ；a.利用下述公式(1)计算地表空气对流换热系数α α = 7. 371+6. 43vair0·75 (1)式中，α为地表空气对流换热系数，单位W/(m2 ·Κ) ；Vair为风速传感器实时测量的 地表风速，单位m/s ；b.利用公式(2)计算土壤导热系数入/I = (11.33 log^-3.4)100·649""-2 (2)式中，λ为土壤导热系数，单位W/(m2· K) ;Pd为土壤密实度，单位kg/m3;rh为湿 度传感器实时测量的土壤相对湿度。(4)首先，选取初始电流11，单位A ；(5)调用有限元温度场计算程序，计算电流Il下的电力电缆绝缘层最高温度Tl， 单位K;(6)当Tl小于电力电缆绝缘层长期工作绝缘层耐受温度Tin时，且Tin-Tl > 0. 3K，则取 12 = 211，转至第(7)步；当 Tl 大于 Tin 时，且 Tl-Tin >0. 3K，则取 12 = 0. 511， 转至第(7)步；当ITl-TinI彡0.5K时，Il即为所求电力电缆载流量，退出；(7)调用有限元温度场计算程序，计算电流12下的电力电缆绝缘层最高温度T2 ；(8)当T2-Tin >0. 3K时，转至第(9)步；当T2_Tin彡0· 3K时，12即为所求 电力电缆载流量，退出；(9)取I = 12+(12-11) · (Tin_T2) / (T2-T1)，调用有限元温度场计算程序，计算电 流I下的电力电缆绝缘层最高温度T ；(10)当 T-Tin > 0. 3K 时，取 Il = 12，12 = I,Tl = T2，T2 = Τ，转至上述第(9) 步；当IT-TinI彡0.3Κ时，I即为所求电力电缆载流量，退出。在本实施例中，收敛条件0. 3Κ由工程实际和温度测量误差决定。温度测量精度为 0. 1-0. 2Κ，结合工程实际要求，在本实施例中收敛条件取0. 3Κ。以石家庄夏季为例，实测地表空气温度为35度，风速为1. 5m/s, 土壤湿度为 18. 5%，由此计算得土壤导热系数为0. 8ff/(m2 ·K)，地表空气对流换热系数为16W/(m2 ·Κ)， 图2-1所示土壤直埋单回路一字型排列800mm2YJLW02XLPE电力电缆，单端接地，埋地深度 1000mm,电缆间距200mm，计算得电缆载流量为892A。图2-2所示排管敷设800mm2YJLW02XLPE电力电缆的计算载流量为858A。图2-3所示沟槽敷设800mm2YJLW02XLPE电力电缆的计算载流量为1085A。本发明涉及的温度场有限元计算法为现有技术(见上述背景技术部分)，下面对温度场有限元计算法介绍如下由于地下电力电缆长度相对于横截面近似为无穷大，因而地下电力电缆群的稳态 温度场属于二维稳态导热问题。有热源区域(如电缆导体、绝缘介质、金属屏蔽层和铠装层)的温度控制方程如下 式(3)所示。+(3)
8x2 Sy2式中，qv为体积发热率，W/m3。无热源固体介质区域(如电力电缆其他层、土壤和排管等)的温度控制方程如式 ⑷所示。 对于排管敷设和沟槽敷设两种方式下，电力电缆和周围固体介质间存在空气，空
气存在受热自然对流的特性，这部分区域的温度场计算需耦合求解连续性方程、动量方程
和能量方程。分别如下式(5)、(6)、(7)所示。 du Sv
2+2 = 0(5)
dx dy 式中，U、ν为流场速度向量在χ和y轴的分量，单位为m/St
.du du、 dp ,d2u 32m、^ N 八
dx dy dx dx2 Sy2
,dv 3v、 dp ,δ2ν 32v、彻 ^ 、 . λ‘
dx dy dy 8χ2 dy2
式中， ；为流体参考温度；α为体积膨胀系数，单位为IT15P为流体密度，单位为
kg/m3 ；ρ为流场的压力，单位为Pa ； η为流体动力粘度，单位为Pa · S。
Γ δτ θτ .rd2T a2r、Λ( Λw —+ ν —+ —) = 0(7)
dx dy dx2 dy2求解温度场需要确定整个温度场的边界，土壤直埋敷设情况下，假定深层土壤 度恒定(即图2-1中的下边界)，如下式(8)所示
d2T B1T &+ = °
■ (8)
τ (χ, y) \r =Tsoil式中，Tstjil为深层土壤温度，取383K。图2-1中两侧边界的法向梯度为0，即远离电力电缆区域土壤的温度场不受电力 电缆发热的影响，如下式(9)所示
d2T | d2T Γ π dx2 dy2^
Sn 1
9
/V 图2-1中的上边界为地表边界，土壤中的热量通过对流换热的形式扩散到空气 中，如下式(10)所示
8 在上式(8)、(9)和(10)中，Γ为积分边界。对于土壤直埋电力电缆，整个场域进行剖分后，利用伽略金公式，可得有限元计算方程，如下式(11)所示KT = P(11)式中，K是与导热系数和单元形状有关的全场域单元矩阵，T是全场域单元温度矩 阵，P是与发热有关的全场域单元矩阵。通常整个场域采用小三角形进行剖分，剖分结果如图3所示。每个三角形单元上 的温度T是坐标χ和y的线性函数，表示为T = ai+a2X+a3y， 、a2和a3由三角形三个顶点上的温度值决定。利用该线性函数 对公式(3)、(8)、(9)、(10)的变分式延三角形单元进行积分，可得每个单元的k值和ρ值， 最后集成得上式(11)，其展开形式为 用高斯消去法对上式(12)进行求解，即可解得每个节点的温度值，从而得到整个 场域的温度场分布和电力电缆的最高温度。对于排管敷设(见图2-2)和沟槽敷设电力电缆(见图2-3)，整个场域进行剖分 后，需要用迭代方法进行计算。列出固体介质区域和空气流体区域的接触耦合边界条件，其中耦合边界上的条件 需满足以下两个条件耦合边界w上温度连续
TJ1 = T,
w I 2
(13)
耦合边界w上的热流密度连续
Qnll = Qnl2(14)假定耦合边界上的温度分布，对其中一个区域1(如固体介质)进行求解，得出耦 合边界上的局部热流密度和温度梯度，然后应用式(12)和式(13)求解另一个区域2 (如空 气介质)，以得出耦合边界上新的温度分布。再以此分布作为区域1的输入，重复上述计算 直到收敛，最后给出整个场域的温度场分布，并得到电力电缆绝缘的最高温度。对于排管敷设和沟槽敷设两种情况，都在电力电缆和外部土壤之间存在空气层， 其温度场计算过程包含固体和流体的耦合计算，可以采用分区计算，边界耦合方式进行整 个场域的温度场计算。分区计算、边界耦合方法的实施步骤是1.分别对电力电缆本体、空气、土壤中的物理问题建立控制方程；2.列出每个区域的边界条件，其中固体区域和空气区域耦合边界上的条件为耦合 边界上温度连续和热流密度连续，即满足上式(13)和式(14)。
9
3.对电力、空气和土壤均按传导传热进行有限元求解(即整场按与土壤直埋的方 式进行求解)，得到初始温度场分布。4.根据上述第3步的计算结果，求解空气域内的对流扩散方程和辐射方程，以得 出耦合边界上新的温度分布。5.以上述第4步计算得固体和流体的边界温度作为电力电缆和土壤的输入，求解 传导控制方程，计算电力电缆和土壤区域的温度场。6.重复上述第3-5步，直到满足收敛条件为止。即当温度和流体速度两者的变化 量均小于10,时，即可中止迭代。给出最终的温度场分布和电力电缆的最高温度。由上述可知，只要给定空气温度、空气风速、太阳光强、土壤湿度，就可计算出地表 对流换热系数，土壤的导热系数，从而计算出整个场域的温度场分布，最终计算出电缆绝缘 层的最高温度。然后通过迭代方法，就可预测出当前环境条件下的电力电缆的载流量。
权利要求
基于环境因子监测和有限元的地下电力电缆载流量在线预测系统，其特征在于它由用于测量空气温度的温度传感器(1)、风速传感器(2)、湿度传感器(3)、单片机系统(4)、电力部门的调度主机(6)组成；所述温度传感器(1)、风速传感器(2)、湿度传感器(3)的输出端分别接所述单片机系统的相应输入端，所述单片机系统通过GSM网络(5)与所述调度主机(6)相连接；所述温度传感器(1)和风速传感器(2)安装在所述电力电缆的正上方，所述湿度传感器(3)安装在所述电力电缆沿线的土壤中。
2.根据权利要求1所述的基于环境因子监测和有限元的地下电力电缆载流量在线预 测方法，其特征在于具体步骤如下(1)建模a.土壤直埋式电缆电力电缆距地面lm，电力电缆两侧边界距电缆10m，计算中，所述两侧边界法向方向热 流变化率为0，即土壤温度不再变化；电力电缆下面至土壤深层边界距电缆10m，取土壤深 层边界温度不变，且温度值为土壤深层温度即383K ；b.排管式电缆排管顶部距地面lm，电力电缆敷设于排管内，排管两侧边界距排管10m，计算中，所述 两侧边界法向方向热流变化率为0，即土壤温度不再变化；排管下面至土壤深层边界距排 管10m，取土壤深层边界温度不变，且温度值为土壤深层温度即383K ；c.沟槽式电缆沟槽深距地面lm，电力电缆敷设于沟槽内，室内沟槽顶部常敷设铁板，室外沟槽顶部常 敷设水泥板，模型两侧边界距沟槽10m，计算中，所述两侧边界法向方向热流变化率为0，即 土壤温度不再变化；沟槽下面至土壤深层边界距沟槽10m，取土壤深层边界温度不变，且温 度值为土壤深层温度即383K ；(2)剖分将电力电缆及其周围土壤剖分为小的单元；在直埋式、排管式和沟槽式三种方式下，采用三角形或四边形单元将整个模型剖分为 小的单元；(3)读取实时测量的下述数据空气温度T&、空气风速v&、土壤湿度rh；并计算地表 空气对流换热系数α和土壤的导热系数入；a.利用下述公式(1)计算地表空气对流换热系数αα = 7. 371+6. 43vair0·75(1)式中，α为地表空气对流换热系数，单位W/(πι2·Κ) ；为风速传感器实时测量的地表 风速，单位m/s ；b.利用公式(2)计算土壤导热系数入Λ = (11.33 log rh - 3.4)100.649内 (2)式中，λ为土壤导热系数，单位W/(m2· K) ;Pd为土壤密实度，单位kg/m3 ；rh为湿度传感器实时测量的土壤相对湿度。(4)首先，选取初始电流II，单位A；(5)调用有限元温度场计算程序，计算电流Il下的电力电缆绝缘层最高温度Tl，单位K ；(6)当Tl小于电力电缆绝缘层长期工作绝缘层耐受温度Tin时，且Tin-Tl> 0. 3K，则取12 = 211，转至第(7)步；当Tl大于Tin时，且Tl-Tin > 0. 3K，则取12 = 0. 511，转至 第⑵步；当Tl-Tin彡0.5K时，Il即为所求电力电缆载流量，退出；(7)调用有限元温度场计算程序，计算电流12下的电力电缆绝缘层最高温度T2；(8)当T2-Tin>0. 3K时，转至第(9)步；当T2-Tin彡0. 3K时，12即为所求电力 电缆载流量，退出；(9)取I= I2+(I2-I1) · (Tin-T2)/(T2-Tl)，调用有限元温度场计算程序，计算电流I 下的电力电缆绝缘层最高温度T ；(10)当T-Tin > 0. 时，取 Il = 12，12 = I,Tl = T2，T2 = Τ，转至上述第(9)步； 当IT-TinI彡0.3Κ时，I即为所求电力电缆载流量，退出。
全文摘要
本发明涉及一种基于环境因子监测和有限元的地下电力电缆载流量在线预测系统及其方法。本发明所述系统由用于测量空气温度的温度传感器、风速传感器、湿度传感器、单片机系统、电力部门的调度主机组成；本发明所述方法是将实时测量的温度、风速、湿度等数据与有限元温度场计算相结合来预测地下电力电缆的载流量。本发明的有益效果是能够较准确的在线预测地下电力电缆的载流量，电力部门据此能够最大限度发挥其输送电能的能力，并且能够在保证安全可靠运行的前提下，不会出现大马拉小车的资源浪费现象。
文档编号G01R31/00GK101900773SQ20101020947
公开日2010年12月1日 申请日期2010年6月25日 优先权日2010年6月25日
发明者梁永春 申请人:河北科技大学