一种用于输电线路动态增容的线路温度辨识方法
一种用于输电线路动态增容的线路温度辨识方法
【技术领域】
[0001] 本发明涉及一种辨识方法，具体涉及一种用于输电线路动态增容的线路温度辨识 方法。
【背景技术】
[0002] 电力负荷的迅猛增长，使得一些输电线路接近或达到热稳定限值，限制了输电线 路的输送容量，影响了供电能力，采用新建或改造线路是增强供电能力的主要办法，但新建 线路面临投资巨大、建设周期长、输电走廊资源紧张、短路电流超标等问题和困难。输电线 路增容技术可以在不进行线路改造的情况下，合理提升输电断面热稳定限值，而热稳定限 值是目前制约经济发达地区线路输送能力的主要约束因素，因而输电线路动态增容技术近 期得到了学者、电网公司的重点关注，并逐渐有其投入实际运行的实例。
[0003] 传统的调度能源管理系统通过SCADA(数据采集与监控系统）对现场运行设备进 行监视和控制。在采集数据的种类上，主要包括电网潮流（如支路功率、母线电压等）以及 电网运行频率等电气量信息。随着智能电网的快速发展，输变电设备在线监测技术的发展 为调度中心提供了设备的各种非电气量信息，如设备状态数据和（微）气象数据等，设备状 态数据包括导线温度等。（微）气象数据如输电线路在线温度、线路所在区域环境温度、风 速、日照以及线路载流量等。这些输电线路在线监测信息接入为输电线路动态增容的实际 应用提供了数据基础和技术支撑。
[0004] 导线温度信息是输电线路增容过程安全性的重要评估指标，其准确可信是动态增 容执行的前提和安全保障。目前，由于产品缺乏统一标准、外部环境多变、通讯干扰等原因， 目前的输电线路在线监测数据的实时型、准确性和稳定性方面都相对存在不足，以目前的 监测信息作为动态增容的基础存在较大的安全隐患。
【发明内容】
[0005] 为了弥补现有技术的不足，本发明提供一种用于输电线路动态增容的线路温度辨 识方法，实现了对导线温度的预测和辨识，在保证系统稳定运行和设备安全前提下，提高输 电线路的输送容量和输电线路监测信息的可信度。
[0006] 本发明的目的是采用下述技术方案实现的：
[0007] -种用于输电线路动态增容的线路温度辨识方法，所述方法包括：
[0008] (1)获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CM/E格 式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；
[0009] (2)预处理CM/E格式监测数据；
[0010] (3)构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行辨识。
[0011] 优选的，所述步骤（1)的控制参数，包括连续采样次数和采样间隔；所述监控信 息，包括系统信息、线路信息、杆塔信息和采集信息。
[0012] 优选的，所述步骤（2)预处理的具体步骤包括：
[0013] 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测 数据是否可用；
[0014] 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范 围的不良数据。
[0015] 进一步地，所述判断某一时间窗口的CIM/E格式监测数据是否可用包括，以该时 间窗口中连续测点的n个CIM/E格式监测数据作为输电线路温度辨识的采样数据，n多2 ; 并以最大采样数据作为各测点的温度预测值；同一输电线路上各温度预测值满足输电线路 热平衡方程，且同一时刻的温度预测值偏差小于所述预设阀值TDiff_。
[0016] 进一步地，所述输电线路热平衡方程包括，当导线电流或者外界环境发生变化时， 导线温度随之改变；其关系构成输电线路热平衡方程：
[0018] 式⑴中，M为单位长度导线的质量；Cp为导线热容系数，X/(kg*°C);
为温度 对时间的导数；Q。为导线的对流散热功率，W/m 导线的辐射散热功率，W/m;〇3为导线的 日照吸热功率，W/m;1为导线电流，A;R(T。）为温度T。时导线的交流电阻，D/m ;
[0019] 当导线温度在某一时刻达到稳定时
?为0,其表达式为：
[0020] 0 =Qs+I2R(TC)-(QC+Qr) (2)；
[0021] 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Qs，获得架空 导线载流量i2r(t。）；
[0022] Qc+Qr=Qs+I2R(Tc) (3) 〇
[0023] 进一步地，根据导线表面的辐射系数e、环境温度!；和导线温度T。，计算所述导线 的辐射散热功率A，其表达式为：
[0024] Qr= 31eD〇 [(Tc+273)4-(Ta+273)4] (4)
[0025] 式⑷中，D为导线直径，m;若所述导线为光亮新线，e的取值范围为0. 23~ 0.43;若为旧线或者涂黑色防腐剂的导线：则e的取值范围为0.90~0.95; 〇为斯蒂 芬-包尔兹曼常数5. 67X10_8WAm2 ?K4)。
[0026] 进一步地，采用简化摩尔根公式计算导线的空气对流散热功率Q。，其表达式如下：
[0027] Qc=AEuJi(Tc-Ta) (5)
[0028] 式（5)中，人为与导体接触的空气膜导热系数，取值范围为?rSEu 为欧拉数，其表达式如下：
[0030] 式（6)中，&为雷诺数，其值可表示为：
[0031] Re= 1. 644X10 9VD[Ta+0. 5(Tc_Ta)r178 (7)
[0032] 式（7)中，V为风速，m/s;D为导线直径。
[0033] 进一步地，根据导线表面的吸热系数a、导线的直径D和日照辐射强度S，计算所 述导线的日照吸热功率Qs;
[0034]Qs=aDS(8)
[0035] 进一步地，当所述架空导线载流量I2R(T。）发生集肤效应时，其交流电阻R(T。）随 导线温度变化，其表达式为：
[0036]R(TC) = (l+k)Rd= (l+k)R20[l+a20(Tc-20)] (9)
[0037] 式（9)中，Rd为温度为T。时导线的直流电阻，D/m;a2(|为20°C的导线材料温度 系数，对于铝的取值为〇. 004031/°C;k为集肤效应系数，若导线截面小于等于400_2取 0. 0025,若大于 400mm2则取 0. 01。
[0038] 进一步地，设采样间隔为At，输电线路的采集装置个数为m，m多1 ;则连续采样时 长为（n-1)XAt，温度实测数据量为mXn个，温度实测数据为；
[0040]与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；
[0042]当采样间隔At较小时，式（1)可简化为：
[0044] 优选的，根据温度预测值与实测值的差值加权平方和作为加权最小二乘的目标函 数，构建导线温度辨识模型，其表达式为：
[0046] 式（11)中，Tm/，为第j个温度采集装置第i时刻的温度量测，I为相应的温度预 测值，为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置；
[0047] 约束条件为：
[0049]式（12)中，1彡i彡m，1彡j彡n-1 ; 为第j+1个温度采集装置第i时刻的温 度预测值；l、&和义分别为第i个采集装置上导线的日照吸热功率、对流散热功率和 辐射散热功率；I2R(T。）为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。
[0050] 进一步地，一次辨识结束后，输出辨识结果；如果所述终止辨识提示为"是"，则终 止辨识程序；若为"否"，则等待采样间隔At后进行下一次辨识。
[0051] 与最接近的现有技术比，本发明提供的技术方案的有益效果是：
[0052] 该方法采用遵循标准的文件接口，增加了软件的通用性；采用数据预处理和数据 修正功能，提升了辨识的辨识效率和可信度；
[0053] 采用动态温度变化模型，并利用连续多采样周期量测数据作为数据源，并且以加 权最小二乘的目标函数，通过权重的合理调整实现不同量测装置精度、不同可信度的量测 数据对目标函数的合理影响，提高辨识结果的准确性和实用性。
[0054] 利用导线热平衡特性建立采集信息的估计和辨识算法，实现对导线温度采集信息 合理辨识和估计，为线路增容的过程控制提供技术基础，保障线路动态增容过程的安全性。
【附图说明】
[0055] 图1为本发明提供的输电线路导线温度在线辨识方法流程图；
[0056]图2为本发明提供的输电线路动态增容系统硬件配置示意图；
[0057] 图3为本发明提供的单位时间内导线温度实测值与预测值的示意图。
【具体实施方式】
[0058] 下面结合附图对本发明的【具体实施方式】作进一步的详细说明。
[0059] 如图1所示，一种用于输电线路动态增容的线路温度辨识方法，所述方法包括：
[0060] (1)获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格 式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；所述步骤 (1)的控制参数，包括连续采样次数和采样间隔；所述监控信息，包括系统信息、线路信息、 杆塔彳目息和米集彳目息。
[0061]
[0063] (2)预处理CM/E格式监测数据；
[0064] 所述步骤（2)预处理的具体步骤包括：
[0065] 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测 数据是否可用；
[0066] 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范 围的不良数据。所述可信值范围，包括历史最大值、历史最小值、前一个月极大值、前一个月 极小值、最大升变化率、最大降变化率等；根据以上特征值检测所述CIM/E格式监测数据是 否在合理范围；
[0067] 当最大变化率、最大变化幅值，遥测变化幅值和变化率超过历史最高值或其某一 倍数，则判断为不良数据。此外还包括额定值检查，当数据超过额定值的k倍；k一般可取 1.1~1.3,则判断为不良数据。
[0068] 并对数据中的空点、毛刺、突变、阶跃等进行分析优化与修正，以提升后续辨识的 可信度。其中，空点一般采用前后一帧数据插值的办法完成填充。
[0069] 毛刺指明显超过正常值范围，包括超出可信值范围或是超过额定值的k倍等情 况，类似于空点的插值处理；
[0070] 突变用于通过平滑的方式执行插值修改将突变的前后几帧数据。
[0071] 所述判断某一时间窗口的CIM/E格式监测数据是否可用包括，以该时间窗口中连 续测点的n个CIM/E格式监测数据作为输电线路温度辨识的采样数据，n多2;并以最大 采样数据作为各测点的温度预测值；同一输电线路上各温度预测值满足输电线路热平衡方 程，且同一时刻的温度预测值偏差小于所述预设阀值TDiff_。
[0072] 所述输电线路热平衡方 程包括，当导线电流或者外界环境发生变化时，导线温度 随之改变；其关系构成输电线路热平衡方程：
[0074]式⑴中，M为单位长度导线的质量，kg;Cp为导线热容系数，Jパkg?°C);
为 温度对时间的导数；Q。为导线的对流散热功率，W/m 导线的辐射散热功率，W/m;QS为导 线的日照吸热功率，W/m;1为导线电流，A;R(T。）为温度T。时导线的交流电阻，D/m ;
[0075] 当导线温度在某一时刻达到稳定时
为0,其表达式为：
[0076] 0 =Qs+I2R(TC)-(QC+Qr) (2)；
[0077] 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Qs，获得架空 导线载流量i2r(t。）；
[0078]Qc+Qr=Qs+I2R(Tc) (3)〇
[0079] 根据导线表面的辐射系数e、环境温度Ta和导线温度T。，计算所述导线的辐射散 热功率义，其表达式为：
[0080]Qr= 31eD〇 [(Tc+273)4-(Ta+273)4] (4)
[0081] 式⑷中，D为导线直径，m;若所述导线为光亮新线，e的取值范围为0. 23~ 0.43;若为旧线或者涂黑色防腐剂的导线：则e的取值范围为0.90~0.95; 〇为斯蒂 芬-包尔兹曼常数5. 67X10_8WAm2 ?K4)。
[0082] 采用简化摩尔根公式计算导线的空气对流散热功率Q。，其表达式如下：
[0083]Qc=AEuJi(Tc-Ta) (5)
[0084] 式（5)中，人为与导体接触的空气膜导热系数，取值范围为?rSEu 为欧拉数，其表达式如下：
[0085]En =0.65/^°-2+0.23/1;i0-61 (6)
[0086] 式（6)中，&为雷诺数，其值可表示为：
[0087] Re= 1. 644X10 9VD[Ta+0. 5(TC-Ta)r178 (7)
[0088] 式（7)中，V为风速，m/s;D为导线直径。
[0089] 根据导线表面的吸热系数a、导线的直径D和日照辐射强度S，计算所述导线的日 照吸热功率Qs;
[0090] Qs= a DS (8)
[0091] 当所述架空导线载流量I2R(T。）发生集肤效应时，其交流电阻R(T。）随导线温度变 化，其表达式为：
[0092]R(TC) = (l+k)Rd= (l+k)R20[l+a20(Tc-20)] (9)
[0093] 式（9)中，Rd为温度为T。时导线的直流电阻，D/m;a2(|为20°C的导线材料温度 系数，对于铝的取值为〇. 004031/°C;k为集肤效应系数，若导线截面小于等于400_2取 0. 0025,若大于 400mm2则取 0. 01。
[0094] 如图3所示，设采样间隔为At，输电线路的采集装置个数为m，m多1;则连续采样 时长为（n-1)XAt，温度实测数据量为mXn个，温度实测数据为；
[0096]与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；
[0098] 当采样间隔At较小时，式（1)可简化为：
[0100] 如图2所示，（3)构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行 辨识。
[0101] 根据温度预测值与实测值的差值加权平方和作为加权最小二乘的目标函数，构建 导线温度辨识模型，其表达式为：
[0103] 式（11)中，Tm,,为第j个温度采集装置第i时刻的温度量测，Tp_、，为相应的温度预 测值，为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置；
[0104] 约束条件为：
[0106] 式（12)中，1彡i彡m，1彡j彡n-1 ; \ ,+1为第j+1个温度采集装置第i时刻的温 度预测值；I、义和义分别为第i个采集装置上导线的日照吸热功率、对流散热功率和 辐射散热功率；I2R(T。）为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。
[0107] -次辨识结束后，输出辨识结果；如果所述终止辨识提示为"是"，则终止辨识程 序；若为"否"，则等待采样间隔At后进行下一次辨识。
[0108] 最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽 管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然 可以对本发明的【具体实施方式】进行修改或者等同替换，而未脱离本发明精神和范围的任何 修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。
【主权项】
1. 一种用于输电线路动态增容的线路温度辨识方法，其特征在于，所述方法包括： (1) 获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格式监 测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息； (2) 预处理CM/E格式监测数据； (3) 构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行辨识。2. 根据权利要求1所述的方法，其特征在于，所述步骤（1)的控制参数，包括连续采样 次数和采样间隔；所述监控信息，包括系统信息、线路信息、杆塔信息和采集信息。3. 根据权利要求1所述的方法，其特征在于，所述步骤（2)预处理的具体步骤包括： 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测数据 是否可用； 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范围的 不良数据。4. 根据权利要求3所述的方法，其特征在于，所述判断某一时间窗口的CIM/E格式监测 数据是否可用包括，以该时间窗口中连续测点的η个CIM/E格式监测数据作为输电线路温 度辨识的采样数据，η多2 ;并以最大采样数据作为各测点的温度预测值；同一输电线路上 各温度预测值满足输电线路热平衡方程，且同一时刻的温度预测值偏差小于所述预设阀值 TDifT-。5. 根据权利要求4所述的方法，其特征在于，所述输电线路热平衡方程包括，当导线电 流或者外界环境发生变化时，导线温度随之改变；其关系构成输电线路热平衡方程：(1) 式（1)中，M为单位长度导线的质量；Cp为导线热容系数，JAkg · °C ) ; f为温度对 时间的导数；Q。为导线的对流散热功率，W/m ;Qr为导线的辐射散热功率，W/m ;QS为导线的日 照吸热功率，W/m ;1为导线电流，A ;R(T。）为温度T。时导线的交流电阻，Ω/m ; 当导线温度在某一时刻达到稳定时f为〇,其表达式为： at 0. Qs+I2R(TC)-(QC+Qr) (2); 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Q s，获得架空导线 载流量I2R(Tc); Qc+Qr= Qs+I2R(Tc) ⑶。6. 根据权利要求4所述的方法，其特征在于，根据导线表面的辐射系数ε、环境温度Ta和导线温度Τ。，计算所述导线的辐射散热功率A，其表达式为： Qr= π ε?σ [(Tc+273)4-(Ta+273)4] (4) 式（4)中，D为导线直径，m ;若所述导线为光亮新线，ε的取值范围为0. 23~0. 43 ; 若为旧线或者涂黑色防腐剂的导线：则ε的取值范围为0.90~0.95; 〇为斯蒂芬-包尔 兹曼常数 5. 67 X KT8WAm2 · K4)。7. 根据权利要求4所述的方法，其特征在于，采用简化摩尔根公式计算导线的空气对 流散热功率Q。，其表达式如下： 0=入 E J (Tc-Ta) (5) 式（5)中，λ为与导体接触的空气膜导热系数，取值范围为〇. 02585W · πΓ1 · IT15Eu为 欧拉数，其表达式如下： L1 = 0. 65/^0-2 + 0. ^ (H) 式（6)中，&为雷诺数，其值可表示为： Re= 1.644X 10 9VD [Ta+0. 5 (Tc-Tjr1.78 (7) 式（7)中，V为风速，m/s ;D为导线直径。8. 根据权利要求4所述的方法，其特征在于，根据导线表面的吸热系数α、导线的直径 D和日照辐射强度S，计算所述导线的日照吸热功率Qs; Qs= a DS (8) 〇9. 根据权利要求4所述的方法，其特征在于，当所述架空导线载流量12R (Τ。）发生集肤 效应时，其交流电阻R(T。）随导线温度变化，其表达式为： R(Tc) = (1+k) Rd= (1+k) R20 [1+a20 (Tc-20)] (9) 式（9)中，Rd为温度为T。时导线的直流电阻，Ω/πι;α 2(|为20°C的导线材料温度系数， 对于铝的取值为0. 004031/°C ;k为集肤效应系数，若导线截面小于等于400mm2取0. 0025， 若大于400mm2则取0. 01。10. 根据权利要求4所述的方法，其特征在于，设采样间隔为Λ t，输电线路的采集装置 个数为m，m多1 ;则连续采样时长为（n-1) X At，温度实测数据量为mXn个，温度实测数 据为；与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；当采样间隔△ t较小时，式（1)可简化为：(10)〇11. 根据权利要求1所述的方法，其特征在于，根据温度预测值与实测值的差值加权平 方和作为加权最小二乘的目标函数，构建导线温度辨识模型，其表达式为：(11) 式（11)中，Tm,,为第j个温度采集装置第i时刻的温度量测，Tew为相应的温度预测 值，Wu为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置； 约束条件为：(12； 式（12)中，1彡i彡m，1彡j彡n-1 ; \(i+u为第j+Ι个温度采集装置第i时刻的温度预 测值；I、Mp &分别为第i个采集装置上导线的日照吸热功率、对流散热功率和辐射 散热功率；I2R(Tc)为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。12.根据权利要求4所述的方法，其特征在于，一次辨识结束后，输出辨识结果；如果所 述终止辨识提示为"是"，则终止辨识程序；若为"否"，则等待采样间隔At后进行下一次辨 识。
【专利摘要】本发明涉及一种用于输电线路动态增容的线路温度辨识方法，所述方法包括：获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；预处理CIM/E格式监测数据；构建导线温度辨识模型，对预处理后的CIM/E格式监测数据进行辨识；从而提升输电线路监测信息的可信度，在保证系统稳定运行和设备安全前提下，提高了输电线路的输送容量。
【IPC分类】G01K13/00
【公开号】CN104897304
【申请号】CN201510342596
【发明人】冯树海, 王勇, 姚建国, 朱斌, 苏大威, 杨胜春, 杨争林, 姜宁, 於益军, 李峰, 石飞, 刘俊, 王礼文, 汤必强
【申请人】中国电力科学研究院, 国家电网公司, 江苏省电力公司, 江苏省电力公司南京供电公司
【公开日】2015年9月9日
【申请日】2015年6月18日
【技术领域】
[0001] 本发明涉及一种辨识方法，具体涉及一种用于输电线路动态增容的线路温度辨识 方法。
【背景技术】
[0002] 电力负荷的迅猛增长，使得一些输电线路接近或达到热稳定限值，限制了输电线 路的输送容量，影响了供电能力，采用新建或改造线路是增强供电能力的主要办法，但新建 线路面临投资巨大、建设周期长、输电走廊资源紧张、短路电流超标等问题和困难。输电线 路增容技术可以在不进行线路改造的情况下，合理提升输电断面热稳定限值，而热稳定限 值是目前制约经济发达地区线路输送能力的主要约束因素，因而输电线路动态增容技术近 期得到了学者、电网公司的重点关注，并逐渐有其投入实际运行的实例。
[0003] 传统的调度能源管理系统通过SCADA(数据采集与监控系统）对现场运行设备进 行监视和控制。在采集数据的种类上，主要包括电网潮流（如支路功率、母线电压等）以及 电网运行频率等电气量信息。随着智能电网的快速发展，输变电设备在线监测技术的发展 为调度中心提供了设备的各种非电气量信息，如设备状态数据和（微）气象数据等，设备状 态数据包括导线温度等。（微）气象数据如输电线路在线温度、线路所在区域环境温度、风 速、日照以及线路载流量等。这些输电线路在线监测信息接入为输电线路动态增容的实际 应用提供了数据基础和技术支撑。
[0004] 导线温度信息是输电线路增容过程安全性的重要评估指标，其准确可信是动态增 容执行的前提和安全保障。目前，由于产品缺乏统一标准、外部环境多变、通讯干扰等原因， 目前的输电线路在线监测数据的实时型、准确性和稳定性方面都相对存在不足，以目前的 监测信息作为动态增容的基础存在较大的安全隐患。
【发明内容】
[0005] 为了弥补现有技术的不足，本发明提供一种用于输电线路动态增容的线路温度辨 识方法，实现了对导线温度的预测和辨识，在保证系统稳定运行和设备安全前提下，提高输 电线路的输送容量和输电线路监测信息的可信度。
[0006] 本发明的目的是采用下述技术方案实现的：
[0007] -种用于输电线路动态增容的线路温度辨识方法，所述方法包括：
[0008] (1)获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CM/E格 式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；
[0009] (2)预处理CM/E格式监测数据；
[0010] (3)构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行辨识。
[0011] 优选的，所述步骤（1)的控制参数，包括连续采样次数和采样间隔；所述监控信 息，包括系统信息、线路信息、杆塔信息和采集信息。
[0012] 优选的，所述步骤（2)预处理的具体步骤包括：
[0013] 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测 数据是否可用；
[0014] 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范 围的不良数据。
[0015] 进一步地，所述判断某一时间窗口的CIM/E格式监测数据是否可用包括，以该时 间窗口中连续测点的n个CIM/E格式监测数据作为输电线路温度辨识的采样数据，n多2 ; 并以最大采样数据作为各测点的温度预测值；同一输电线路上各温度预测值满足输电线路 热平衡方程，且同一时刻的温度预测值偏差小于所述预设阀值TDiff_。
[0016] 进一步地，所述输电线路热平衡方程包括，当导线电流或者外界环境发生变化时， 导线温度随之改变；其关系构成输电线路热平衡方程：
[0018] 式⑴中，M为单位长度导线的质量；Cp为导线热容系数，X/(kg*°C);
为温度 对时间的导数；Q。为导线的对流散热功率，W/m 导线的辐射散热功率，W/m;〇3为导线的 日照吸热功率，W/m;1为导线电流，A;R(T。）为温度T。时导线的交流电阻，D/m ;
[0019] 当导线温度在某一时刻达到稳定时
?为0,其表达式为：
[0020] 0 =Qs+I2R(TC)-(QC+Qr) (2)；
[0021] 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Qs，获得架空 导线载流量i2r(t。）；
[0022] Qc+Qr=Qs+I2R(Tc) (3) 〇
[0023] 进一步地，根据导线表面的辐射系数e、环境温度!；和导线温度T。，计算所述导线 的辐射散热功率A，其表达式为：
[0024] Qr= 31eD〇 [(Tc+273)4-(Ta+273)4] (4)
[0025] 式⑷中，D为导线直径，m;若所述导线为光亮新线，e的取值范围为0. 23~ 0.43;若为旧线或者涂黑色防腐剂的导线：则e的取值范围为0.90~0.95; 〇为斯蒂 芬-包尔兹曼常数5. 67X10_8WAm2 ?K4)。
[0026] 进一步地，采用简化摩尔根公式计算导线的空气对流散热功率Q。，其表达式如下：
[0027] Qc=AEuJi(Tc-Ta) (5)
[0028] 式（5)中，人为与导体接触的空气膜导热系数，取值范围为?rSEu 为欧拉数，其表达式如下：
[0030] 式（6)中，&为雷诺数，其值可表示为：
[0031] Re= 1. 644X10 9VD[Ta+0. 5(Tc_Ta)r178 (7)
[0032] 式（7)中，V为风速，m/s;D为导线直径。
[0033] 进一步地，根据导线表面的吸热系数a、导线的直径D和日照辐射强度S，计算所 述导线的日照吸热功率Qs;
[0034]Qs=aDS(8)
[0035] 进一步地，当所述架空导线载流量I2R(T。）发生集肤效应时，其交流电阻R(T。）随 导线温度变化，其表达式为：
[0036]R(TC) = (l+k)Rd= (l+k)R20[l+a20(Tc-20)] (9)
[0037] 式（9)中，Rd为温度为T。时导线的直流电阻，D/m;a2(|为20°C的导线材料温度 系数，对于铝的取值为〇. 004031/°C;k为集肤效应系数，若导线截面小于等于400_2取 0. 0025,若大于 400mm2则取 0. 01。
[0038] 进一步地，设采样间隔为At，输电线路的采集装置个数为m，m多1 ;则连续采样时 长为（n-1)XAt，温度实测数据量为mXn个，温度实测数据为；
[0040]与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；
[0042]当采样间隔At较小时，式（1)可简化为：
[0044] 优选的，根据温度预测值与实测值的差值加权平方和作为加权最小二乘的目标函 数，构建导线温度辨识模型，其表达式为：
[0046] 式（11)中，Tm/，为第j个温度采集装置第i时刻的温度量测，I为相应的温度预 测值，为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置；
[0047] 约束条件为：
[0049]式（12)中，1彡i彡m，1彡j彡n-1 ; 为第j+1个温度采集装置第i时刻的温 度预测值；l、&和义分别为第i个采集装置上导线的日照吸热功率、对流散热功率和 辐射散热功率；I2R(T。）为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。
[0050] 进一步地，一次辨识结束后，输出辨识结果；如果所述终止辨识提示为"是"，则终 止辨识程序；若为"否"，则等待采样间隔At后进行下一次辨识。
[0051] 与最接近的现有技术比，本发明提供的技术方案的有益效果是：
[0052] 该方法采用遵循标准的文件接口，增加了软件的通用性；采用数据预处理和数据 修正功能，提升了辨识的辨识效率和可信度；
[0053] 采用动态温度变化模型，并利用连续多采样周期量测数据作为数据源，并且以加 权最小二乘的目标函数，通过权重的合理调整实现不同量测装置精度、不同可信度的量测 数据对目标函数的合理影响，提高辨识结果的准确性和实用性。
[0054] 利用导线热平衡特性建立采集信息的估计和辨识算法，实现对导线温度采集信息 合理辨识和估计，为线路增容的过程控制提供技术基础，保障线路动态增容过程的安全性。
【附图说明】
[0055] 图1为本发明提供的输电线路导线温度在线辨识方法流程图；
[0056]图2为本发明提供的输电线路动态增容系统硬件配置示意图；
[0057] 图3为本发明提供的单位时间内导线温度实测值与预测值的示意图。
【具体实施方式】
[0058] 下面结合附图对本发明的【具体实施方式】作进一步的详细说明。
[0059] 如图1所示，一种用于输电线路动态增容的线路温度辨识方法，所述方法包括：
[0060] (1)获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格 式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；所述步骤 (1)的控制参数，包括连续采样次数和采样间隔；所述监控信息，包括系统信息、线路信息、 杆塔彳目息和米集彳目息。
[0061]
[0063] (2)预处理CM/E格式监测数据；
[0064] 所述步骤（2)预处理的具体步骤包括：
[0065] 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测 数据是否可用；
[0066] 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范 围的不良数据。所述可信值范围，包括历史最大值、历史最小值、前一个月极大值、前一个月 极小值、最大升变化率、最大降变化率等；根据以上特征值检测所述CIM/E格式监测数据是 否在合理范围；
[0067] 当最大变化率、最大变化幅值，遥测变化幅值和变化率超过历史最高值或其某一 倍数，则判断为不良数据。此外还包括额定值检查，当数据超过额定值的k倍；k一般可取 1.1~1.3,则判断为不良数据。
[0068] 并对数据中的空点、毛刺、突变、阶跃等进行分析优化与修正，以提升后续辨识的 可信度。其中，空点一般采用前后一帧数据插值的办法完成填充。
[0069] 毛刺指明显超过正常值范围，包括超出可信值范围或是超过额定值的k倍等情 况，类似于空点的插值处理；
[0070] 突变用于通过平滑的方式执行插值修改将突变的前后几帧数据。
[0071] 所述判断某一时间窗口的CIM/E格式监测数据是否可用包括，以该时间窗口中连 续测点的n个CIM/E格式监测数据作为输电线路温度辨识的采样数据，n多2;并以最大 采样数据作为各测点的温度预测值；同一输电线路上各温度预测值满足输电线路热平衡方 程，且同一时刻的温度预测值偏差小于所述预设阀值TDiff_。
[0072] 所述输电线路热平衡方 程包括，当导线电流或者外界环境发生变化时，导线温度 随之改变；其关系构成输电线路热平衡方程：
[0074]式⑴中，M为单位长度导线的质量，kg;Cp为导线热容系数，Jパkg?°C);
为 温度对时间的导数；Q。为导线的对流散热功率，W/m 导线的辐射散热功率，W/m;QS为导 线的日照吸热功率，W/m;1为导线电流，A;R(T。）为温度T。时导线的交流电阻，D/m ;
[0075] 当导线温度在某一时刻达到稳定时
为0,其表达式为：
[0076] 0 =Qs+I2R(TC)-(QC+Qr) (2)；
[0077] 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Qs，获得架空 导线载流量i2r(t。）；
[0078]Qc+Qr=Qs+I2R(Tc) (3)〇
[0079] 根据导线表面的辐射系数e、环境温度Ta和导线温度T。，计算所述导线的辐射散 热功率义，其表达式为：
[0080]Qr= 31eD〇 [(Tc+273)4-(Ta+273)4] (4)
[0081] 式⑷中，D为导线直径，m;若所述导线为光亮新线，e的取值范围为0. 23~ 0.43;若为旧线或者涂黑色防腐剂的导线：则e的取值范围为0.90~0.95; 〇为斯蒂 芬-包尔兹曼常数5. 67X10_8WAm2 ?K4)。
[0082] 采用简化摩尔根公式计算导线的空气对流散热功率Q。，其表达式如下：
[0083]Qc=AEuJi(Tc-Ta) (5)
[0084] 式（5)中，人为与导体接触的空气膜导热系数，取值范围为?rSEu 为欧拉数，其表达式如下：
[0085]En =0.65/^°-2+0.23/1;i0-61 (6)
[0086] 式（6)中，&为雷诺数，其值可表示为：
[0087] Re= 1. 644X10 9VD[Ta+0. 5(TC-Ta)r178 (7)
[0088] 式（7)中，V为风速，m/s;D为导线直径。
[0089] 根据导线表面的吸热系数a、导线的直径D和日照辐射强度S，计算所述导线的日 照吸热功率Qs;
[0090] Qs= a DS (8)
[0091] 当所述架空导线载流量I2R(T。）发生集肤效应时，其交流电阻R(T。）随导线温度变 化，其表达式为：
[0092]R(TC) = (l+k)Rd= (l+k)R20[l+a20(Tc-20)] (9)
[0093] 式（9)中，Rd为温度为T。时导线的直流电阻，D/m;a2(|为20°C的导线材料温度 系数，对于铝的取值为〇. 004031/°C;k为集肤效应系数，若导线截面小于等于400_2取 0. 0025,若大于 400mm2则取 0. 01。
[0094] 如图3所示，设采样间隔为At，输电线路的采集装置个数为m，m多1;则连续采样 时长为（n-1)XAt，温度实测数据量为mXn个，温度实测数据为；
[0096]与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；
[0098] 当采样间隔At较小时，式（1)可简化为：
[0100] 如图2所示，（3)构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行 辨识。
[0101] 根据温度预测值与实测值的差值加权平方和作为加权最小二乘的目标函数，构建 导线温度辨识模型，其表达式为：
[0103] 式（11)中，Tm,,为第j个温度采集装置第i时刻的温度量测，Tp_、，为相应的温度预 测值，为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置；
[0104] 约束条件为：
[0106] 式（12)中，1彡i彡m，1彡j彡n-1 ; \ ,+1为第j+1个温度采集装置第i时刻的温 度预测值；I、义和义分别为第i个采集装置上导线的日照吸热功率、对流散热功率和 辐射散热功率；I2R(T。）为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。
[0107] -次辨识结束后，输出辨识结果；如果所述终止辨识提示为"是"，则终止辨识程 序；若为"否"，则等待采样间隔At后进行下一次辨识。
[0108] 最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽 管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然 可以对本发明的【具体实施方式】进行修改或者等同替换，而未脱离本发明精神和范围的任何 修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。
【主权项】
1. 一种用于输电线路动态增容的线路温度辨识方法，其特征在于，所述方法包括： (1) 获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格式监 测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息； (2) 预处理CM/E格式监测数据； (3) 构建导线温度辨识模型，对预处理后的CM/E格式监测数据进行辨识。2. 根据权利要求1所述的方法，其特征在于，所述步骤（1)的控制参数，包括连续采样 次数和采样间隔；所述监控信息，包括系统信息、线路信息、杆塔信息和采集信息。3. 根据权利要求1所述的方法，其特征在于，所述步骤（2)预处理的具体步骤包括： 首先，根据测点、预设阈值和终止辨识提示判断某一时间窗口的CIM/E格式监测数据 是否可用； 其次，根据历史数据设定可信值范围，删除所述CIM/E格式监测数据中超出该范围的 不良数据。4. 根据权利要求3所述的方法，其特征在于，所述判断某一时间窗口的CIM/E格式监测 数据是否可用包括，以该时间窗口中连续测点的η个CIM/E格式监测数据作为输电线路温 度辨识的采样数据，η多2 ;并以最大采样数据作为各测点的温度预测值；同一输电线路上 各温度预测值满足输电线路热平衡方程，且同一时刻的温度预测值偏差小于所述预设阀值 TDifT-。5. 根据权利要求4所述的方法，其特征在于，所述输电线路热平衡方程包括，当导线电 流或者外界环境发生变化时，导线温度随之改变；其关系构成输电线路热平衡方程：(1) 式（1)中，M为单位长度导线的质量；Cp为导线热容系数，JAkg · °C ) ; f为温度对 时间的导数；Q。为导线的对流散热功率，W/m ;Qr为导线的辐射散热功率，W/m ;QS为导线的日 照吸热功率，W/m ;1为导线电流，A ;R(T。）为温度T。时导线的交流电阻，Ω/m ; 当导线温度在某一时刻达到稳定时f为〇,其表达式为： at 0. Qs+I2R(TC)-(QC+Qr) (2); 通过计算导线的对流散热功率Q。、辐射散热功率A和日照吸热功率Q s，获得架空导线 载流量I2R(Tc); Qc+Qr= Qs+I2R(Tc) ⑶。6. 根据权利要求4所述的方法，其特征在于，根据导线表面的辐射系数ε、环境温度Ta和导线温度Τ。，计算所述导线的辐射散热功率A，其表达式为： Qr= π ε?σ [(Tc+273)4-(Ta+273)4] (4) 式（4)中，D为导线直径，m ;若所述导线为光亮新线，ε的取值范围为0. 23~0. 43 ; 若为旧线或者涂黑色防腐剂的导线：则ε的取值范围为0.90~0.95; 〇为斯蒂芬-包尔 兹曼常数 5. 67 X KT8WAm2 · K4)。7. 根据权利要求4所述的方法，其特征在于，采用简化摩尔根公式计算导线的空气对 流散热功率Q。，其表达式如下： 0=入 E J (Tc-Ta) (5) 式（5)中，λ为与导体接触的空气膜导热系数，取值范围为〇. 02585W · πΓ1 · IT15Eu为 欧拉数，其表达式如下： L1 = 0. 65/^0-2 + 0. ^ (H) 式（6)中，&为雷诺数，其值可表示为： Re= 1.644X 10 9VD [Ta+0. 5 (Tc-Tjr1.78 (7) 式（7)中，V为风速，m/s ;D为导线直径。8. 根据权利要求4所述的方法，其特征在于，根据导线表面的吸热系数α、导线的直径 D和日照辐射强度S，计算所述导线的日照吸热功率Qs; Qs= a DS (8) 〇9. 根据权利要求4所述的方法，其特征在于，当所述架空导线载流量12R (Τ。）发生集肤 效应时，其交流电阻R(T。）随导线温度变化，其表达式为： R(Tc) = (1+k) Rd= (1+k) R20 [1+a20 (Tc-20)] (9) 式（9)中，Rd为温度为T。时导线的直流电阻，Ω/πι;α 2(|为20°C的导线材料温度系数， 对于铝的取值为0. 004031/°C ;k为集肤效应系数，若导线截面小于等于400mm2取0. 0025， 若大于400mm2则取0. 01。10. 根据权利要求4所述的方法，其特征在于，设采样间隔为Λ t，输电线路的采集装置 个数为m，m多1 ;则连续采样时长为（n-1) X At，温度实测数据量为mXn个，温度实测数 据为；与Tm相对应，各时刻点的温度预测数据量为mXn个，温度预测数据为；当采样间隔△ t较小时，式（1)可简化为：(10)〇11. 根据权利要求1所述的方法，其特征在于，根据温度预测值与实测值的差值加权平 方和作为加权最小二乘的目标函数，构建导线温度辨识模型，其表达式为：(11) 式（11)中，Tm,,为第j个温度采集装置第i时刻的温度量测，Tew为相应的温度预测 值，Wu为相应测点的权重；该权重根据量测装置精度和量测历史可信状况进行设置； 约束条件为：(12； 式（12)中，1彡i彡m，1彡j彡n-1 ; \(i+u为第j+Ι个温度采集装置第i时刻的温度预 测值；I、Mp &分别为第i个采集装置上导线的日照吸热功率、对流散热功率和辐射 散热功率；I2R(Tc)为架空导线载流量;M为单位长度导线的质量；Cp为导线热容系数。12.根据权利要求4所述的方法，其特征在于，一次辨识结束后，输出辨识结果；如果所 述终止辨识提示为"是"，则终止辨识程序；若为"否"，则等待采样间隔At后进行下一次辨 识。
【专利摘要】本发明涉及一种用于输电线路动态增容的线路温度辨识方法，所述方法包括：获取输电线路的在线监测数据，并将其转化为遵循IEC61970标准的CIM/E格式监测数据，保存至内存；所述CIM/E格式监测数据包括控制参数和监控信息；预处理CIM/E格式监测数据；构建导线温度辨识模型，对预处理后的CIM/E格式监测数据进行辨识；从而提升输电线路监测信息的可信度，在保证系统稳定运行和设备安全前提下，提高了输电线路的输送容量。
【IPC分类】G01K13/00
【公开号】CN104897304
【申请号】CN201510342596
【发明人】冯树海, 王勇, 姚建国, 朱斌, 苏大威, 杨胜春, 杨争林, 姜宁, 於益军, 李峰, 石飞, 刘俊, 王礼文, 汤必强
【申请人】中国电力科学研究院, 国家电网公司, 江苏省电力公司, 江苏省电力公司南京供电公司
【公开日】2015年9月9日
【申请日】2015年6月18日
最新回复(0)