CN108734342A - 区域天气预报用于地区电网架空线路热定值概率预测方法 - Google Patents
区域天气预报用于地区电网架空线路热定值概率预测方法 Download PDFInfo
- Publication number
- CN108734342A CN108734342A CN201810399366.9A CN201810399366A CN108734342A CN 108734342 A CN108734342 A CN 108734342A CN 201810399366 A CN201810399366 A CN 201810399366A CN 108734342 A CN108734342 A CN 108734342A
- Authority
- CN
- China
- Prior art keywords
- meteorological
- transmission line
- deviation
- definite value
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 88
- 238000013277 forecasting method Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012549 training Methods 0.000 claims abstract description 13
- 239000000470 constituent Substances 0.000 claims abstract description 7
- 239000004020 conductor Substances 0.000 claims description 20
- 238000005286 illumination Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 5
- 230000017525 heat dissipation Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000009423 ventilation Methods 0.000 claims description 5
- 230000007774 longterm Effects 0.000 claims description 3
- 230000002950 deficient Effects 0.000 claims 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims 1
- 230000005611 electricity Effects 0.000 abstract description 14
- 238000005259 measurement Methods 0.000 description 8
- 238000012544 monitoring process Methods 0.000 description 6
- 238000005096 rolling process Methods 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 102000000584 Calmodulin Human genes 0.000 description 2
- 108010041952 Calmodulin Proteins 0.000 description 2
- 241001269238 Data Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000739 chaotic effect Effects 0.000 description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241001123248 Arma Species 0.000 description 1
- 241000127225 Enceliopsis nudicaulis Species 0.000 description 1
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 210000004744 fore-foot Anatomy 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了区域天气预报用于地区电网架空线路热定值概率预测方法，包括：对区域天气预报历史数据作为训练集进行预处理；在训练集中提取相同季度下气象要素的区域天气预报值以及对应小时内的架空线路的微气象数据；计算在各气象区间内的各气象要素预报值与微气象之间的偏差，建立气象偏差的样本集合；基于气象偏差样本集合，估计各季度下各个气象要素在各预报区间内的气象偏差的概率密度分布；估计预测时段内架空线路热定值的概率密度分布。本申请方法可依据区域天气预报信息给出预测时段内热定值可能的波动范围及概率分布，预测结果对帮助调度高效利用输电线路载荷能力，提高电力系统安全经济运行水平具有实际意义。
Description
技术领域
本发明涉及电网技术领域，特别是涉及区域天气预报用于地区电网架空线路热定值概率预测方法。
背景技术
当前，随着发电、负荷的不断增长、新能源发电的大规模接入使电网载荷能力面临挑战，成为制约电力系统经济运行及新能源消纳的重要因素之一，如何充分利用现有输电元件载荷能力，提高其使用效率已成为电力工作者关注的焦点问题。
架空输电线路的最大允许载流量(热定值)与气象环境密切相关，目前工程中广泛使用的静态热定值(Static thermal rating,STR)是在假设较为不利的气象条件(高气温、低风速，强日照)同时发生的情况下计算得到的，具有明显的保守性。对此，上世纪70年代美国学者Davis提出了动态热定值(Dynamic thermal rating,DTR)的技术框架，旨在根据输电线路微气象量测实时计算其热定值。DTR体现了输电元件载荷能力的时变性，研究表明在多数情况下架空线路DTR均显著高于STR。
在此基础上，文献“王孔森，盛戈皞.基于径向基神经网络的输电线路动态容量在线预测[J].电网技术，2013，37(6)：1719-1725.”基于架空线路微气象量测的历史数据，利用径向基神经网络进行风速、光照强度及环境温度变化规律的在线学习和预测，进而实现了输电线路动态热定值的点预测。
文献“任丽佳，江秀臣，盛戈皞，等.输电线路允许输送容量的混沌预测[J].中国电机工程学报，2009，29(25)：86-91”则基于动态热定值的历史数据，运用混沌预测理论实现架空线路热定值的短期预测。
文献“Junpeng Zhan.Time series modelling for dynamic thermal rating ofoverhead lines.IEEE Transactions on Power Systems,2017,32(3):2172-2182”进一步考虑到对热定值预测结果不确定性的描述，基于微气象历史数据，结合累积式自回归积分滑动平均(ARIMA)时间序列模型以及多项式分位点回归的方法实现了架空线路热定值的概率密度预测。
此外，文献“Fulin Fan,Keith Bell,David Infield.Probabilistic real-timethermal rating forecasting for overhead lines by conditionallyheteroscedastic auto-regressive models.IEEE Transactions on Power Delivery，2017，32(4)：1881-1890”采用自回归条件异方差模型滚动预测架空线路半小时内每十分钟的各气象要素，并实现架空线路动态热定值的概率密度预测。
以上热定值预测的研究为DTR技术与电网运行调控的结合奠定了基础，对实现输电线路载荷能力的充分利用具有重要意义。然而，现有架空线路热定值预测的研究主要依据架空线路微气象的在线监测数据，虽然架空线在线监测的相关技术已较为成熟(如DTR系统)但其实施也必然增加相关量测、通信设备的投资以及日常的运维成本，这在一定程度上限制了架空线路动态热定值计算及预测技术的广泛应用。
目前，随着全国范围内气象监测数据的日趋完善及开放平台的建设，区域天气预报(如各城市天气预报或网格化区域的数值天气预报)已能够提供多种气象要素的预报信息且易于获取，而区域内或附近的微气象变化与区域整体天气变化之间存在密切关联。
发明内容
为了解决现有技术的不足，本发明提供了区域天气预报用于地区电网架空线路热定值概率预测方法，本发明的方法可摆脱对在线监测设备的依赖，便于工程应用。
区域天气预报用于地区电网架空线路热定值概率预测方法，包括：
对区域天气预报历史数据作为训练集进行预处理剔除天气预报不良数据；
在训练集中提取相同季度下气象要素环境温度、风速、风向以及光照强度的区域天气预报值以及对应小时内的架空线路的微气象数据；
计算在各气象区间内的各气象要素预报值与微气象之间的偏差，建立气象偏差的样本集合；
基于气象偏差样本集合，估计各季度下各个气象要素在各预报区间内的气象偏差的概率密度分布；
在获得气象偏差的概率密度分布曲线后，拟合得到不同季度下各个气象要素在各预报区间内的气象偏差的概率密度函数；
根据预测时段内各个气象要素的区域天气预报值选择其所在区间的气象偏差的概率密度函数，通过随机抽样获得各个气象要素偏差，与预测时段天气预报叠加，获得预测时段内架空线路关键线档的微气象数据，计算预测时段内的架空线路热定值；
估计预测时段内架空线路热定值的概率密度分布。
进一步优选的技术方案，对区域天气预报历史数据作为训练集进行预处理时，采用模糊等价矩阵的聚类分析方法，以区域天气预报误差为特征对区域天气预报历史数据进行聚类分析，确定不良数据检测门槛值以辨识并剔除天气预报不良数据。
进一步优选的技术方案，计算在各气象区间内的各气象要素预报值与微气象之间的偏差时，具体的计算方式为：
式中，向量yi为某气象要素区域天气预报在第i个区间内时的气象偏差统计向量；为某气象要素第i个预报区间中第n个小时内第m个微气象量测数据，设某气象要素在该季度中共有n个小时的区域天气预报值落在第i区间内；I为单位向量，λi,n为某气象要素第i个预报区间中第n小时的区域天气预报值，b为区间的个数。
进一步优选的技术方案，基于气象偏差样本集合，估计各季度下各个气象要素在各预报区间内的气象偏差的概率密度分布时，使用高斯函数作为Parzen窗估计的窗口函数，概率密度估计式的表示如下：
其中，pi(x)为某气象要素第i区间中气象偏差x的概率密度；k(·)为高斯核函数；yi.u为某气象要素第i预报区间中第u个偏差样本；h为带宽系数。
进一步优选的技术方案，计算预测时段内的架空线路热定值，计算公式为：
其中，Tmax为架空导线长期连续运行的最大允许温度，qs为日照吸热量(w/m)，主要与日射强度有关，qc为单位长度导体的空气对流散热量(w/m)，与风速、风向以及导体与环境温差等因素有关；qr为热辐射散热量(w/m)，主要与导体和环境的温差有关，R(Tma)为温度Tma下单位导线长度电阻。
进一步优选的技术方案，在获得气象偏差的概率密度分布曲线后，利用高斯函数拟合得到不同季度下4个气象要素在各预报区间内的气象偏差的概率密度函数。
进一步优选的技术方案，计算预测时段内的架空线路热定值时，反复对预测时段内的微气象数据进行抽样并计算热定值，直至热定值的概率密度分布基本不变时计算结束。
进一步优选的技术方案，各气象要素按区域天气预报数值分区间统计气象偏差的分布。
进一步优选的技术方案，区域天气预报数据与周边架空线微气象数据偏差的分布呈近似正态分布。
进一步优选的技术方案，利用Parzen窗法估计预测时段内架空线路热定值的概率密度分布。
与现有技术相比，本发明的有益效果是：
本发明首先结合实测微气象及区域天气预报数据分析了两者之间的相关性及偏差的分布规律，在此基础上提出了将区域天气预报用于架空线路热定值概率预测的思路及实施步骤，最后结合某架空线实测微气象数据及小时级区域天气预报验证了所提出热定值概率预测方法的有效性。在根据历史微气象及区域天气预报数据统计两者偏差分布规律的基础上，本申请方法可依据区域天气预报信息给出预测时段内热定值可能的波动范围及概率分布，预测结果对帮助调度高效利用输电线路载荷能力，提高电力系统安全经济运行水平具有实际意义。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1为3个线档的区域天气预报与微气象的相关系数；
图2(a)-图2(d)为线档1各气象要素偏差的频率分布直方图；
图3(a)-图3(d)为线档2各气象要素偏差的频率分布直方图；
图4(a)-图4(d)为线档3各气象要素偏差的频率分布直方图；
图5区域天气预报用于热定值概率预测的思路；
图6(a)-图6(d)各气象要素偏差的概率密度分布；
图7架空线路热定值概率密度分布；
图8预测时段内不同分位数预测值、静态热定值与实测值；
图9 24小时架空线路动态热定值不同分位数滚动预测；
图10 24小时架空线实际热定值与预测值分位数滚动预测；
图11(a)-图11(d)各气象要素区域天气预报与实际气象数据对比；
图12 2015年DTR预测结果采样与实际热定值对比；
图13 2016年DTR预测结果采样与实际热定值对比；
图14 2017年DTR预测结果采样与实际热定值对比。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本申请使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
区域内或附近的微气象变化与区域整体天气变化之间存在密切关联，本申请利用这一关联关系，在统计分析区域天气预报相对架空线微气象偏差的分布规律的基础上，提出区域天气预报用于输电线路动态热定值概率预测的方法，为在一定程度上摆脱架空线路热定值预测对在线监测系统的依赖提供支持，并为运行人员合理设定架空线路热定值，改善STR的保守性提供依据。
架空输电线路热定值计算模型：根据IEEE标准，架空线路稳态热平衡模型可由式(1)描述：
qi(T(t))+qs(t)＝qc(T(t))+qr(T(t)) (1)
式(1)中，t表示时间(s)，T表示导线温度(℃)；等式左侧qi为考虑导体电阻温度效应下的单位长度导线的电阻发热量(w/m)，与导线载流和温度有关；qs为日照吸热量(w/m)，主要与日射强度有关。等式右侧qc为单位长度导体的空气对流散热量(w/m)，与风速、风向以及导体与环境温差等因素有关；qr为热辐射散热量(w/m)，主要与导体和环境的温差有关。式(1)描述了导体热稳态下吸热与散热之间的平衡关系，式中各吸、散热项的具体计算式如下：
qi(T)＝I2R(T) (2)
qs＝αQsesin(θ)A′ (4)
式(2)～式(6)中，I为架空线路载流(A)，R(T)为温度T下单位导线长度电阻(Ω)，[Tlow,Thigh]为电阻-温度线性关系的存续区间，本申请设定Tlow为25℃,Thigh为75℃；D为导线直径(mm),ε为导体的散热系数,α为导体光照吸收率，A’为单位长度的导体投影面积(mm/m2)，由导线材料及型号决定；kf为空气导热系数(W/(m·℃))；qc1为低风速时计算的单位长度导体的空气对流散热量，qc2适用于高风速。在任何风速下，应使用两个计算式中对流散热量较大者。在零风速下，发生自然对流，应采用qc3计算空气对流散热量；Qse为光照强度(W/m2)，ρf为空气密度(kg/m3)，μf为空气的动态粘度(Pa·s)，θ为有效的太阳光线入射角(°)，与导线所处地理位置及海拔有关；Ta为导线周围环境温度(℃)；Vw为风速(m/s)；kangle为风向因子，和架空线路与风向的夹角有关。
由式(2)～式(6)可见，一旦架空线路架设地点及导线型号确定，输电线路热定值主要取决于环境温度(Ta)、风速(Vw)、风向(kangle)以及光照强度(Qse)4个气象要素。考虑到架空导线沿线气象分布的差异性，通常在分析架空线沿线弧垂及气象分布规律的基础上选择少量关键线档安装在线监测装置(如DTR设备)以节省成本，并以关键线档中的最小热定值代表整条线路的热定值。在已知某关键线档4要素微气象数据的情况下，架空线路热定值可由式(7)计算：
式(7)中Tmax为架空导线长期连续运行的最大允许温度，我国通常取为70℃。
架空线路微气象与区域天气预报的关联关系分析：如上所述，架空线热定值的计算需要关键线档的微气象信息，本节结合实际数据分析架空线微气象与附近区域天气预报的关联关系，以说明将区域天气预报用于架空线路热定值概率预测的可行性。本申请采用的区域天气预报为城市气象站小时级的天气预报，本申请方法也适用于更为精细的网格化、分层化的区域数值天气预报。
城市天气预报主要反映距地高度2～10m范围内的城市及周边整体气象变化情况[18]，而架空线横担高度多在15～50m范围内，受高差及微地形的影响，架空线沿线微气象与区域气象之间必然存在差异。然而，由于架空线关键线档与城市的相对位置、所处的微地形状况及高度一定，且环境温度、光照强度本身具有较强的区域一致性，因此，架空导线微气象与附近区域天气预报之间也必然存在联系，以下结合实测架空线路微气象及附近城市区域天气预报数据分析两者之间的相关性，并发掘两者相对偏差的统计规律。
本申请收集了某220kV架空线路沿线3个关键线档(分别记为1、2、3号线档，横担高度分别为18m、24m和33m，距离附近某城市中心分别约为20km、30km和50km)2013～2017年分钟级微气象数据以及附近城市气象站的小时级历史天气预报数据。其中，以2013～2014年气象数据作为训练集，2015～2017年气象数据作为测试集验证本申请方法。
基于本申请方法对测试集内不同年份(2015年-2017年)热定值概率密度分布进行预测，架空线热定值不同分位数预测结果与实际热定值均值(黑色实线)对比如下所示(在不同季度内共随机抽取3000个点)：
由图12、图13和图14可见，由于夏季光照强度及温度偏高而风速偏低，因此夏季热定值总体看来较低。本申请方法的预测结果较好的反映了架空线动态热定值的变化趋势和波动范围，可为电力系统调度决策提供帮助。
考虑到在特殊情况下，如当气象环境波动剧烈时，天气预报可能出现较大误差(不良数据)。因此，本申请结合模糊等价矩阵的聚类分析方法，对历史区域天气预报数据进行预处理，剔除不良数据。在此基础上，经统计，训练集中3个线档不同气象要素区域天气预报与微气象的相关系数变化如图1所示。可见，各气象要素区域天气预报与微气象的相关性较强(相关系数高于0.5)，其中环境温度和光照强度的区域天气预报与微气象的相关性最强，且随着线档与城市之间距离的增大相关系数下降并不明显。相比之下，风速、风向的区域天气预报值与微气象之间的相关系数则有较为明显的下降，但从数值上看相关性仍然较强。
为进一步分析区域天气预报与周边架空线档微气象之间偏差(下文中简称气象偏差)的分布规律，本申请对各气象要素按区域天气预报数值分区间统计气象偏差的分布。例如，根据历史区域天气预报中夏季气温数据，将气温预报分为[10℃,20℃]、(20℃,30℃]及(30℃,45℃]3个区间，其中第3区间代表夏季高气温区间，而后分别统计区域天气预报在此3个区间内的气象偏差分布情况。同理，其它气象要素，如风速分为[0m/s,3m/s]、(3m/s,6m/s]、(6m/s,9m/s]及9m/s以上4个区间；风向按与导线夹角分为[0,90°]、(90°,180°]、(180°,270°]及(270°，360°]4个区间；光照强度分为[0w/m2,200w/m2]、(200w/m2,400w/m2]、(400w/m2,600w/m2]、(600w/m2,800w/m2]、800w/m2以上5个区间。
图2(a)-图2(d)、图3(a)-图3(d)和图4(a)-图4(d)分别给出了架空线沿线3个关键线档在区域天气预报处于夏季高气温(第4区间)、低风速(第1区间)、强光照(第4区间)以及小风向夹角(第1区间)的情况下，气象偏差的分布情况。
可见，距离城市较近的各气象要素偏差呈近似正态分布，以距离市中心20km的线档1为例，该极端气象条件下环境温度的区域天气预报与微气象偏差主要集中在[-1℃,2℃](偏差落在该区间的概率大于0.95)，风速的偏差主要集中在[-2m/s,2m/s]，光照强度的偏差主要集中在[-40w/m2，40w/m2]；当关键线档距离市区较远时，风速、风向的气象偏差则出现偏离正态分布的趋势。以上分析表明，区域天气预报与周边微气象之间存在较强相关性，且在区域周边一定范围内的气象偏差分布具有规律性，可为热定值的概率预测提供参考。
区域天气预报用于架空线路动态热定值概率密度预测的方法：由式(7)可知，在设定的最大允许运行温度下，架空导线的热定值随4个气象要素的变化而变化。从架空线在线监测的角度，目前对架空线路微气象要素的测量周期可达分钟级。然而，从调度决策的角度则需要在一定决策周期内(如实时调度决策周期为15～60min)具有代表性的、唯一的热定值预测值以构建输电线路载荷能力约束。因此，实现未来决策周期内架空线路热定值可能的变化范围及其概率分布预测，合理确定热定值，对高效利用输电设备进而提高电力系统的安全经济运行水平具有重要意义。
区域天气预报与周边架空线微气象偏差的分布呈近似正态分布。因此，在通过历史数据学习该分布规律的基础上，结合区域天气预报即可获得未来架空线各气象要素微气象的分布情况，进而实现热定值的概率预测。预测思路如下图5所示。
按此思路，本申请结合分钟级架空线路实测微气象历史数据及附近城市气象站的小时级天气预报历史数据，提出基于区域天气预报的架空线热定值概率密度预测方法，具体计算步骤如下：
(1)首先对区域天气预报历史数据进行预处理。本申请采用模糊等价矩阵的聚类分析方法[19]，以区域天气预报误差为特征对区域天气预报历史数据进行聚类分析，确定不良数据检测门槛值以辨识并剔除天气预报不良数据；
本申请利用模糊聚类分析法辨识并剔除区域天气预报中的不良数据。首先根据区域历史天气预报及实测气象的历史数据，统计训练集中区域天气预报误差并标准化，并以此作为特征指标结合模糊聚类分析法建立模糊关系矩阵，对良数据和不良数据进行分类，获得不良数据检测门槛值，超过检测门槛值即判为不良数据。
图11(a)-图11(d)给出了在数据预处理前区域天气预报与气象站实测气象数据的对比(在全年不同季节内随机采样2000次)。其中，蓝色曲线为区域气象实际值，红色曲线为区域气象站预报数据。表1、表2分别给出了预处理前后的区域天气预报准确率。
表1.数据预处理前各气象要素未来1小时城市气象站预报准确率(％)
表2.数据预处理后各气象要素预报数据准确率％
由以上数据预处理前的准确率情况可知，本申请建立气象偏差概率模型所采用的未来1小时城市气象站天气预报在大多数情况下是可信的。在此基础上，本申请对少数不良数据进行了辨识和剔除，进一步提高了数据源的准确性。
(2)在训练集中提取相同季度下4个气象要素的区域天气预报值以及对应小时内的架空线路的微气象数据，而后利用式(8)分别计算在各气象区间i(i＝1…b)内的各气象要素预报值与微气象之间的偏差，建立气象偏差的样本集合。
式中，向量yi为某气象要素区域天气预报在第i个区间内时的气象偏差统计向量；为某气象要素第i个预报区间中第n个小时内第m个微气象量测数据(设某气象要素在该季度中共有n个小时的区域天气预报值落在第i区间内)；I为单位向量，λi,n为某气象要素第i个预报区间中第n小时的区域天气预报值；。
(3)基于气象偏差样本集合，采用Parzen窗法分别估计各季度下4个气象要素在各预报区间内的气象偏差的概率密度分布。本申请使用高斯函数作为Parzen窗估计的窗口函数，概率密度估计式的表示如下：
其中，pi(x)为某气象要素第i区间中气象偏差x的概率密度；k(·)为高斯核函数；yi.u为某气象要素第i预报区间中第u个偏差样本；h为带宽系数；
(4)在获得气象偏差的概率密度分布曲线后，利用高斯函数拟合得到不同季度下4个气象要素在各预报区间内的气象偏差的概率密度函数；
(5)根据预测时段内4个气象要素的区域天气预报值选择其所在区间的气象偏差的概率密度函数，通过随机抽样获得4个气象要素偏差，与预测时段天气预报叠加，获得预测时段内架空线路关键线档的微气象数据，并利用式(7)计算预测时段内的架空线路热定值。反复对预测时段内的微气象数据进行抽样并计算热定值，直至热定值的概率密度分布基本不变时计算结束，经大量计算发现当模拟次数达到5000次后基本能够获得稳定的热定值概率密度函数。因此，本申请设定图5中N＝5000次时计算结束。
(6)利用Parzen窗法估计预测时段内架空线路热定值的概率密度分布。
为了使得本领域技术人员能够更加清楚地了解本申请的技术方案，以下将结合具体的算例分析详细说明本申请的技术方案。
本申请以预测某220kV输电线路关键线档(距市区约30km)2017年夏季7月15日10点至11点时段内架空线路热定值为例，根据前文所述的预测流程，首先统计训练集(2013～2014年夏季6月1日-8月31日)中4个气象要素预报值在各气象区间内的气象偏差分布，而后根据测试集中预测时段内区域天气预报的数值选择其所在区间的气象偏差概率密度函数。在本算例中，预测时段(2017年7月15日10点到11点)内4要素区域气象预报信息分别为气温26℃，风速2m/s，风向东东南，光照强度730w/m2，据此分别选择对应区间的气象偏差概率密度曲线如图6(a)-图6(d)所示(4气象要素区域天气预报的区间划分如前文所述)。
对以上4个气象要素偏差的概率密度分布曲线进行高斯函数拟合，得到各气象要素偏差的概率密度函数表达式如下：
环境温度：
风速：
风向：
光照强度：
基于概率密度函数(10)～(13)随机抽样预测时段内各气象要素偏差并与区域天气预报数据相加，而后利用式(7)计算架空线路热定值(共抽样5000次)，可得到预测时段内架空线路热定值概率密度分布以及不同分位数如图7所示，整体计算耗时约50s。图7中0、0.1、0.3、0.5为概率值。
预测时段内的真实热定值(分钟级)，不同分位数的预测值及其与STR(STR计算对应的气象条件为：气温35℃，风速0.5m/s，风向垂直于导线，光照强度800w/㎡)的对比如图8所示。
由图8可见，预测时段内的实际热定值高于STR(650A)且波动性较强，若使用期望预测值(0.5分位数)构建输电线路载荷安全约束用于调度决策会带来较高的安全风险。此外，本申请方法的预测结果显示由0、1分位数构成的热定值区间可完全覆盖真实热定值的变化范围。不同分位数的预测结果表明该时段内架空线热定值高于755.3A(0分位数)的概率为1(相比STR提高16.2％)，高于985.5A(0.1分位数)的概率为0.9(相比STR提高61.6％)，高于1138.0A(0.3分位数)的概率为0.7(相比STR提高75.1％)，预测的期望值为1205.6A。以上述预测结果为例，保守起见运行人员可选取0～0.1分位数作为整个预测时段的热定值，即便如此在多数情况下仍可显著改善STR的保守性。
基于每小时更新的区域天气预报数据，对2017年7月15日24小时架空线不同分位数热定值进行滚动预测的结果如图9所示，与实际热定值对比如图10所示。此外，其它年份不同季节的预测结果见附录B。由图9、图10可见，本申请方法预测结果能够较好的反映架空线路动态热定值的变化趋势和波动范围，可帮助调度人员根据风险偏好选择合适的分位数热定值用于调度决策，提高电网利用效率。
本申请针对架空线路热定值的时变性及难以被准确预测的特点，在把握区域天气预报与微气象之间的关联关系的基础上提出区域天气预报用于架空线路热定值概率密度预测方法，有结论如下：(1)不同位置、横担高度的架空线路关键线档的微气象与附近区域天气预报之间存在较强的相关性；(2)本申请方法在充分考虑预测时段内架空线路微气象变化的基础上，结合区域天气预报，能够给出预测时段内热定值概率密度分布和不同分位数下架空线路热定值预测信息，为调度充分利用输电线路载荷能力提供支持；(3)经过对微气象与地区天气预报偏差统计分析建模后，本申请方法可摆脱对在线监测设备的依赖，便于工程应用。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (10)
1.区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，包括：
对区域天气预报历史数据作为训练集进行预处理剔除天气预报不良数据；
在训练集中提取相同季度下气象要素环境温度、风速、风向以及光照强度的区域天气预报值以及对应小时内的架空线路的微气象数据；
计算在各气象区间内的各气象要素预报值与微气象之间的偏差，建立气象偏差的样本集合；
基于气象偏差样本集合，估计各季度下各个气象要素在各预报区间内的气象偏差的概率密度分布；
在获得气象偏差的概率密度分布曲线后，拟合得到不同季度下各个气象要素在各预报区间内的气象偏差的概率密度函数；
根据预测时段内各个气象要素的区域天气预报值选择其所在区间的气象偏差的概率密度函数，通过随机抽样获得各个气象要素偏差，与预测时段天气预报叠加，获得预测时段内架空线路关键线档的微气象数据，计算预测时段内的架空线路热定值；
估计预测时段内架空线路热定值的概率密度分布。
2.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，对区域天气预报历史数据作为训练集进行预处理时，采用模糊等价矩阵的聚类分析方法，以区域天气预报误差为特征对区域天气预报历史数据进行聚类分析，确定不良数据检测门槛值以辨识并剔除天气预报不良数据。
3.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，计算在各气象区间内的各气象要素预报值与微气象之间的偏差时，具体的计算方式为：
式中，向量yi为某气象要素区域天气预报在第i个区间内时的气象偏差统计向量；为某气象要素第i个预报区间中第n个小时内第m个微气象量测数据，设某气象要素在该季度中共有n个小时的区域天气预报值落在第i区间内；I为单位向量，λi,n为某气象要素第i个预报区间中第n小时的区域天气预报值，b为区间的个数。
4.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，基于气象偏差样本集合，估计各季度下各个气象要素在各预报区间内的气象偏差的概率密度分布时，使用高斯函数作为Parzen窗估计的窗口函数，概率密度估计式的表示如下：
其中，pi(x)为某气象要素第i区间中气象偏差x的概率密度；k(·)为高斯核函数；yi.u为某气象要素第i预报区间中第u个偏差样本；h为带宽系数。
5.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，计算预测时段内的架空线路热定值，计算公式为：
其中，Tmax为架空导线长期连续运行的最大允许温度，qs为日照吸热量(w/m)，主要与日射强度有关，qc为单位长度导体的空气对流散热量(w/m)，与风速、风向以及导体与环境温差等因素有关；qr为热辐射散热量(w/m)，主要与导体和环境的温差有关，R(Tma)为温度Tma下单位导线长度电阻。
6.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，在获得气象偏差的概率密度分布曲线后，利用高斯函数拟合得到不同季度下4个气象要素在各预报区间内的气象偏差的概率密度函数。
7.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，计算预测时段内的架空线路热定值时，反复对预测时段内的微气象数据进行抽样并计算热定值，直至热定值的概率密度分布基本不变时计算结束。
8.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，各气象要素按区域天气预报数值分区间统计气象偏差的分布。
9.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，区域天气预报数据与周边架空线微气象数据偏差的分布呈近似正态分布。
10.如权利要求1所述的区域天气预报用于地区电网架空线路热定值概率预测方法，其特征是，利用Parzen窗法估计预测时段内架空线路热定值的概率密度分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810399366.9A CN108734342B (zh) | 2018-04-28 | 2018-04-28 | 区域天气预报用于地区电网架空线路热定值概率预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810399366.9A CN108734342B (zh) | 2018-04-28 | 2018-04-28 | 区域天气预报用于地区电网架空线路热定值概率预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108734342A true CN108734342A (zh) | 2018-11-02 |
| CN108734342B CN108734342B (zh) | 2020-12-25 |
Family
ID=63940086
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810399366.9A Active CN108734342B (zh) | 2018-04-28 | 2018-04-28 | 区域天气预报用于地区电网架空线路热定值概率预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108734342B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109934379A (zh) * | 2019-01-21 | 2019-06-25 | 平安科技（深圳）有限公司 | 空气污染物浓度的预测方法、装置及电子设备 |
| CN110084403A (zh) * | 2019-03-27 | 2019-08-02 | 山东信通电子股份有限公司 | 架空导线载流量的多时段联合概率密度预测方法及系统 |
| CN111340302A (zh) * | 2020-03-05 | 2020-06-26 | 国惠环保新能源有限公司 | 一种基于局部气象预测的精准供热负荷预测方法及系统 |
| CN112580845A (zh) * | 2019-09-30 | 2021-03-30 | 北京金风慧能技术有限公司 | 气象数据处理方法、装置、设备及计算机可读存储介质 |
| CN112733354A (zh) * | 2020-11-06 | 2021-04-30 | 国网山东省电力公司聊城供电公司 | 一种气象要素时间序列模拟方法、系统、介质及电子设备 |
| CN113222276A (zh) * | 2021-05-27 | 2021-08-06 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 覆冰情况预测方法、装置、计算机设备和存储介质 |
| CN114493171A (zh) * | 2021-12-31 | 2022-05-13 | 国网山东省电力公司临沂供电公司 | 一种动态增容设备安装选址方案生成方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008017581A1 (en) * | 2006-08-11 | 2008-02-14 | Abb Research Ltd | Parameter estimation for a thermal model of a power line |
| CN104849776A (zh) * | 2014-12-08 | 2015-08-19 | 国家电网公司 | 一种结合动态修正的电网高低温精细化预警方法 |
| CN105574616A (zh) * | 2015-12-16 | 2016-05-11 | 国网山东省电力公司电力科学研究院 | 一种架空输电线路动态容量预测方法 |
| US20170075035A1 (en) * | 2015-09-11 | 2017-03-16 | Kabushiki Kaisha Toshiba | Probabilistic weather forecasting device, probabilistic weather forecasting method, and non-transitory computer readable medium |
| CN106845663A (zh) * | 2015-12-07 | 2017-06-13 | 中国电力科学研究院 | 一种基于数值天气预报进行热负荷预测的方法 |
| CN107292427A (zh) * | 2017-06-07 | 2017-10-24 | 北京国网富达科技发展有限责任公司 | 输电线路最大允许载流量概率性预测方法及装置 |
| CN107330243A (zh) * | 2017-05-25 | 2017-11-07 | 国网浙江省电力公司电力科学研究院 | 一种基于典型气象年的架空输电线路载流量计算方法 |
-
2018
- 2018-04-28 CN CN201810399366.9A patent/CN108734342B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008017581A1 (en) * | 2006-08-11 | 2008-02-14 | Abb Research Ltd | Parameter estimation for a thermal model of a power line |
| CN104849776A (zh) * | 2014-12-08 | 2015-08-19 | 国家电网公司 | 一种结合动态修正的电网高低温精细化预警方法 |
| US20170075035A1 (en) * | 2015-09-11 | 2017-03-16 | Kabushiki Kaisha Toshiba | Probabilistic weather forecasting device, probabilistic weather forecasting method, and non-transitory computer readable medium |
| CN106845663A (zh) * | 2015-12-07 | 2017-06-13 | 中国电力科学研究院 | 一种基于数值天气预报进行热负荷预测的方法 |
| CN105574616A (zh) * | 2015-12-16 | 2016-05-11 | 国网山东省电力公司电力科学研究院 | 一种架空输电线路动态容量预测方法 |
| CN107330243A (zh) * | 2017-05-25 | 2017-11-07 | 国网浙江省电力公司电力科学研究院 | 一种基于典型气象年的架空输电线路载流量计算方法 |
| CN107292427A (zh) * | 2017-06-07 | 2017-10-24 | 北京国网富达科技发展有限责任公司 | 输电线路最大允许载流量概率性预测方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 张斌等: "随环境条件变化的输电线路输送容量概率建模研究", 《电工电能新技术》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109934379A (zh) * | 2019-01-21 | 2019-06-25 | 平安科技（深圳）有限公司 | 空气污染物浓度的预测方法、装置及电子设备 |
| CN110084403A (zh) * | 2019-03-27 | 2019-08-02 | 山东信通电子股份有限公司 | 架空导线载流量的多时段联合概率密度预测方法及系统 |
| CN110084403B (zh) * | 2019-03-27 | 2022-03-01 | 山东信通电子股份有限公司 | 架空导线载流量的多时段联合概率密度预测方法及系统 |
| CN112580845A (zh) * | 2019-09-30 | 2021-03-30 | 北京金风慧能技术有限公司 | 气象数据处理方法、装置、设备及计算机可读存储介质 |
| CN111340302A (zh) * | 2020-03-05 | 2020-06-26 | 国惠环保新能源有限公司 | 一种基于局部气象预测的精准供热负荷预测方法及系统 |
| CN112733354A (zh) * | 2020-11-06 | 2021-04-30 | 国网山东省电力公司聊城供电公司 | 一种气象要素时间序列模拟方法、系统、介质及电子设备 |
| CN113222276A (zh) * | 2021-05-27 | 2021-08-06 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 覆冰情况预测方法、装置、计算机设备和存储介质 |
| CN113222276B (zh) * | 2021-05-27 | 2023-09-26 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 覆冰情况预测方法、装置、计算机设备和存储介质 |
| CN114493171A (zh) * | 2021-12-31 | 2022-05-13 | 国网山东省电力公司临沂供电公司 | 一种动态增容设备安装选址方案生成方法及系统 |
| CN114493171B (zh) * | 2021-12-31 | 2023-11-28 | 国网山东省电力公司临沂供电公司 | 一种动态增容设备安装选址方案生成方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108734342B (zh) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108734342A (zh) | 区域天气预报用于地区电网架空线路热定值概率预测方法 | |
| Han et al. | A PV power interval forecasting based on seasonal model and nonparametric estimation algorithm | |
| WO2021238505A1 (zh) | 基于联邦学习的区域光伏功率概率预测方法及协同调控系统 | |
| Ohunakin et al. | Generation of a typical meteorological year for north–east, Nigeria | |
| US20140180637A1 (en) | Irradiance mapping leveraging a distributed network of solar photovoltaic systems | |
| CN106779129A (zh) | 一种考虑气象因素的短期电力负荷预测方法 | |
| Lawal et al. | Dynamic line rating forecasting algorithm for a secure power system network | |
| CN109740195B (zh) | 一种基于气象站观测数据的极值台风风速概率分布模型及设计台风风速的评估方法 | |
| Chen et al. | A new method utilizing smart meter data for identifying the existence of air conditioning in residential homes | |
| CN110135649A (zh) | 短期建筑能耗区间预测方法、系统、介质及设备 | |
| CN116205425A (zh) | 一种基于典型数据库的低碳园区冷热电气负荷预测方法 | |
| Tavarov et al. | Model and Algorithm of Electricity Consumption Management for Household Consumers in the Republic of Tajikistan. | |
| CN108615087A (zh) | 一种光伏组件温度短期预测方法及其预测系统 | |
| CN117033927A (zh) | 地铁站碳排放实时监测预测分析方法及预测分析方法 | |
| Sun et al. | Spatio-temporal weather model-based probabilistic forecasting of dynamic thermal rating for overhead transmission lines | |
| Xu et al. | Influence of urban morphological factors on building energy consumption combined with photovoltaic potential: A case study of residential blocks in central China | |
| Wang et al. | Uncertainty modeling method of weather elements based on deep learning for robust solar energy generation of building | |
| Guo et al. | A new approach for interval forecasting of photovoltaic power based on generalized weather classification | |
| CN110084403A (zh) | 架空导线载流量的多时段联合概率密度预测方法及系统 | |
| Cabezon et al. | Comparison of methods for power curve modelling | |
| CN116633016B (zh) | 一种基于高频数据采集技术的电力管理系统和方法 | |
| Zhang et al. | A study on the electric power load of Beijing and its relationships with meteorological factors during summer and winter | |
| Mulero et al. | Data-driven energy resource planning for Smart Cities | |
| CN111144628A (zh) | 一种分布式供能的冷热电负荷预测模型系统及方法 | |
| CN109409597A (zh) | 一种基于成因分析的中长期大洪水预报方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |