一种电缆导体温度确定方法
【技术领域】
[0001] 本发明涉及电力系统运行控制技术，具体涉及一种电缆导体温度确定方法。
【背景技术】
[0002] 随着城市规模的扩大和经济建设的快速发展，城市电缆化率明显提升，电缆线路 的电压等级不断提高，输送能量持续增加，电力电缆的运行管理和监测维护工作变得愈加 重要。
[0003] 电缆导体的温度是电缆的一个重要参数，当电缆在额定负荷下运行时，导体温度 在允许值范围内，一旦过负荷，导体温度将急剧上升，加速绝缘老化，使得泄露电流增加，温 度持续升高，有可能导致绝缘击穿，甚至发生火灾。温度升高是造成电缆线路故障的一个重 要原因，因此可以将其作为表征电缆运行情况的一个重要参数。无论出于电力电缆自身运 行安全考虑，还是从电力系统全局角度出发，都需要对电力电缆导体温度进行实时监测。
[0004] 目前采用的电缆温度监测方式主要有光信号传感系统和电信号传感系统两种。国 内部分城市在电缆线路上应用了分布式光纤测温系统，能够实时监测电缆护套表面的温 度，但要由其推算导体温度就需要建立护套表面温度与导体温度之间的对应关系，这又涉 及到电力电缆的热路和热场问题。
[0005] 电缆导体温度的计算主要有以下几种方法：建立热传学模型，根据红外热像仪拍 摄电缆的表面热图像，反演计算出电缆导体温度；根据电缆等效热路与电路在数学形式上 相同的特点，建立电缆热路模型，通过实测电缆外护套表面温度推算出实时电缆导体温度； 利用数值分析法建立导体温度场，如采用有限元分析方法及有限差分法建立导体的短时暂 态温度场。
[0006] 值得注意的是，电缆外护套测温结果容易受到外界环境的影响，如：有风和无风的 环境下测得的温度会有所不同，这直接影响到导体温度计算的准确性，电缆结构和敷设参 数也会影响计算结果。
[0007] 电缆导体的温度对导体电阻的影响很大，已知导体电阻，可推算出导体温度。
【发明内容】
[0008] 针对现有技术的不足，本发明的目的是提供一种电缆导体温度确定方法，该方法 根据实时采集的电压、电流、有功/无功功率计算导体电阻，进而根据导体电阻推算导体温 度，可为电缆线路过热或过负荷风险辨识及预警提供参考依据。
[0009] 本发明的目的是采用下述技术方案实现的：
[0010] 本发明提供一种电缆导体温度确定方法，其改进之处在于，所述方法包括下述步 骤：
[0011] (1)建立电缆线路首末端电压、电流、有功功率、无功功率与电阻的四元线性回归 模型；
[0012] (2)根据四元线性回归模型确定电缆线路电阻；
[0013] (3)将电缆线路电阻换算为单位长度的导体电阻；
[0014] (4)确定电缆导体温度。
[0015] 进一步地，所述步骤（1)包括：忽略电缆线路首末端电压降的横向分量，电缆线路 首末端的电压损耗等于电压降的纵向分量，表达式如下：
[0016]
[0017] 式中：P。、Q。分别为电缆线路首端有功功率、无功功率；Pn、Qn分别为电缆线路末端 有功功率、无功功率；RS、XS分别为电缆线路的总电阻和总电抗；Rq、Xq分别为引入的虚拟 电阻和电抗；￡为服从正态分布的常数；U。、^分别为电缆线路首端、末端的电压；
[0018] 式1)为以AU为因变量，
为自变量的四元线性回 归模型，RS、XS、Rq、Xq为回归系数；回归系数Rs即为电缆线路电阻。
[0019] 进一步地，所述步骤（2)包括：采用最小二乘法确定四元线性回归模型的回归系 数，进而确定电缆线路电阻：
[0020] 设在时间段[0,T]内，采集间隔为T/m，采集时刻为kt2,. ..，k，共采集m组数据： 线路首端电压[UQ1，UQ2，…，UQJ;线路末端电压[Unl，Un2，…，UJ;线路首端有功功率[P01， PQ2，…，PJ;线路首端无功功率[QM，9。2，…，Q。」；线路末端有功功率[Pnl，Pn2，…，Pnn]; 线路末端无功功率[Qnl，Qn2，…，QJ;
[0021] 其中：1^,1^，…，1]。"分别表示第一个采集时刻、第二个采集时刻以及第m个采集 时刻的线路首端电压；Unl，Un2，…，Unni分别表示第一个采集时刻、第二个采集时刻以及第m 个采集时刻的线路末端电压；PM，Pc2，…，Pcm分别表示第一个采集时刻、第二个采集时刻以 及第m个采集时刻的线路首端有功功率；QM，〇。2，…，Q&分别表示第一个采集时刻、第二个 采集时刻以及第m个采集时刻的线路首端无功功率；Pnl，Pn2，…，Pnni分别表示第一个采集 时刻、第二个采集时刻以及第m个采集时刻的线路末端有功功率；Qnl，Qn2，…，Qnni分别表示 第一个采集时刻、第二个采集时刻以及第m个采集时刻的线路末端有功功率；
[0022] 为表述方便，引入Ylk、Y2k、Y3k、Y4k，不代表任何实际电气量；
[0023] 令：
[0024]
[0025] 式中：k=l，2,. . .m;PQk、QQk分别为第k个采集时刻电缆线路首端有功功率、无功功 率;Pnk、Qnk分别为第k个采集时刻电缆线路末端有功功率、无功功率；
[0026] 贝 1J:AUk =YlkRs+Y2kXs+Y3kReq+Y4kXeq+e;
[0027] 式中：AUk为第k个采集时刻的首末端电压降；RS、Xs分别为电缆线路的总电阻 和总电抗况q、L分别为引入的虚拟电阻和电抗；e为服从正态分布的常数；
[0028] 根据最小二乘法原理，令：
[0032] 有：Ab= B ;
[0033] A表7K-阶偏导矩阵、b表7K回归系数向量、B表7K值向量；
[0034] 则：b = A 屯；
[0035] 式中：A 1表示矩阵A的逆阵；RS即为电缆线路的总电阻。
[0036] 进一步地，所述步骤（3)中，将电缆线路电阻换算为单位长度的导体电阻的表达式 如下：
[0037]
5)；
[0038] 式中：RT表示单位长度的导体电阻；L表示导体长度；RS表示电缆线路电阻。
[0039] 进一步地，所述步骤（4)中，根据单位长度的导体电阻确定电缆导体温度，表达式 如下：
[0040]
6);
[0041] 式中：R2。表示20°C时电缆导体单位长度的电阻，单位：Q/km，通过心计算； P2。表示电缆导体材料在2〇°C时的电阻率，单位：Qgmm2/km;S表示导体标称截面积，单位：mm2;a表示电缆导体材料的电阻温度系数，1/°C;T表示电缆导体温度，°C，RT表示工作温 度T时单位长度的导体电阻。
[0042] 进一步地，所得到的电缆导体温度为电缆导体的等效温度。
[0043] 与现有技术比，本发明达到的有益效果是：
[0044] (1)本发明提供的基于电缆导体电阻的温度确定方法，建立电缆线路首末端电压、 电流、有功功率、无功功率与电阻的四元线性回归模型，采用最小二乘法计算回归系数，进 而获得电缆线路电阻，再根据电阻推算电缆导体实时等效温度。
[0045] (2)本发明避免了复杂的热传学模型或热路模型的建模过程；
[0046] (3)本发明无需测量电缆外护套温度，无需考虑外界环境对温度计算结果的影 响；
[0047] (4)本发明原理简单、计算快速，可为过热或过负荷风险辨识及预警提供参考依 据。
【附图说明】
[0048] 图1是本发明提供的基于电缆电阻估算电缆导体温度的总体流程图。
【具体实施方式】
[0049] 下面结合附图对本发明的【具体实施方式】作进一步的详细说明。
[0050] 本发明提供的基于电缆电阻估算电缆导体温度的总体流程图如图1所示，包括下 述步骤：
[0051] (1)建立电缆线路首末端电压、电流、有功功率、无功功率与电阻的四元线性回归 模型：
[0052] -般情况下，电缆线路首末两端电压相角差不大，可忽略电压降的横向分量，近似 认为电压损耗就等于电压降的纵向分量，表达式如下：
[0053]
[0054] 式中：P。、Q。分别为电缆线路首端有功功率、无功功率;Pn、Qn分别为电缆线路末端 有功功率、无功功率；RS、XS分别为电缆线路的总电阻和总电抗；Rq、Xq分别为引入的虚拟 电阻和电抗；e为服从正态分布的常数；U。、^分别为电缆线路首端、末端的电压；
[0055] 式1)为以AU为因变量，
为自变量的四元线性回 归模型，Rs、Xs、Rq、X%为回归系数；回归系数1即为电缆线路电阻，利用多组量测数据， 采用最小二乘法即可求得电阻Rs。
[0056] (2)根据四元线性回归模型确定电缆线路电阻：采用最小二乘法确定四元线性回 归模型的回归系数，进而确定电缆线路电阻。
[0057] 设在时间段[0,T]内，采集