CN113176420A - 一种针对电网杆塔点的风速预报订正系统 - Google Patents

一种针对电网杆塔点的风速预报订正系统 Download PDF

Info

Publication number
CN113176420A
CN113176420A CN202110172682.4A CN202110172682A CN113176420A CN 113176420 A CN113176420 A CN 113176420A CN 202110172682 A CN202110172682 A CN 202110172682A CN 113176420 A CN113176420 A CN 113176420A
Authority
CN
China
Prior art keywords
point
wind speed
module
forecast
tower
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
Application number
CN202110172682.4A
Other languages
English (en)
Other versions
CN113176420B (zh
Inventor
韩乐琼
徐驰
张永山
官丽
王丙兰
赵留学
王香云
郭鹏
赵晓栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jiutian Jiutian Meteorological Technology Co ltd
State Grid Beijing Electric Power Co Ltd
Original Assignee
Beijing Jiutian Jiutian Meteorological Technology Co ltd
State Grid Beijing Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Jiutian Jiutian Meteorological Technology Co ltd, State Grid Beijing Electric Power Co Ltd filed Critical Beijing Jiutian Jiutian Meteorological Technology Co ltd
Priority to CN202110172682.4A priority Critical patent/CN113176420B/zh
Publication of CN113176420A publication Critical patent/CN113176420A/zh
Application granted granted Critical
Publication of CN113176420B publication Critical patent/CN113176420B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Environmental & Geological Engineering (AREA)
  • Algebra (AREA)
  • Atmospheric Sciences (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请公开一种针对电网杆塔点的风速预报订正系统。所述系统包括地形复杂度最优格点筛选装置，计算最优格点序号并发送至数据重构装置和综合预报装置；数据重构装置，根据最优格点序号提取对应格点的观测和预报值，并发送至订正模型装置；订正模型装置，接收格点观测和预报值，计算模式预报样本阶段的风速频率，并计算累积风频，根据累积风频设置风速阈值并构建大风订正模型和普态订正模型并发送至综合预报装置；综合预报装置，采集模式预报值，根据接收的格点序号，提取对应预报值，根据风速阈值和订正模型计算最终的格点预报订正值。本申请能够获取到最接近真是的杆塔点风速观测信息，并建立误差订正模型有效提高风速预报的准确度。

Description

一种针对电网杆塔点的风速预报订正系统
技术领域
本申请涉及电力生产和气象预测的交叉领域，尤其涉及一种针对电网杆塔点的风速预报订正系统。
背景技术
电网系统中最关键最易受损坏的部分就是输电线路，其架设后的风险除了少部分人为破坏外，更多来自于外部气象因素的影响，特别是受大风天气影响导致的异物侵扰、风偏等灾害形成的跳线放电，线路跳闸、短路等，以及因为该类天气导致后续的检查、施工、施救等工作的不通畅。基于数值预报模式输出的风力要素是分析是否会产生异物侵扰、风偏现象的关键，也是预判输电线路是否会出现气象灾害的主要参数。
输电线路架设于杆塔间，基于现阶段的数值模式空间精度，难以描绘出每条线路上每一米所受的风力强度大小及其变化，因此针对不同杆塔点所在位置的精准风力预报及相应的预警判断，是判断输电线路上是否产生异物侵扰、风偏现象的关键。
近年来，数值预报模式虽能较好地预测出边界层、近地层等不同垂直层、不同尺度的大气环流特征，并输出格点化预报产品，但产品的准确率和空间分辨率有限，难以针对更小区域或某一指定点位的风力要素，根据其地形特征做出更准确的预测。而针对大风这种相对极端的天气现象，数值模式更是难以预报出其峰值特征，因此实现从网格化预测产品到服务对象点位的产品后释用方法即为订正技术，该技术是获取预报与实况偏差的重要途径。目前主要有两类订正技术：(1) 通过计算一段时间内，实况和预报的平均值误差，将预报加上该误差值，作为订正预报结果；(2)通过构建实况和预报的数学统计经验关系的线性回归技术。
前者的先决条件是模式预报偏差较稳定，长期保持恒定的相对实况偏大或者偏小的状态，否则如果长期偏大的预报，误差值为负的情况下，在大风天气时，预报本身偏小，再减去误差后，会导致预报风速更小，对后续的异物、风偏等的预警有较大影响，无法做到准确的预报、预警。而第二种方法构建了二者的相关关系，但属于纯数理统计的范畴，同样在大风事件这种极端情况时，较难预报出峰值。
因此，需要一种能够通过订正技术对杆塔点的风力要素做出精确预测，且最重要的是能在大风事件中较好的预报出峰值，从而能为后释用的异物预警和风偏预警提供更好的基础数据的方法。。
发明内容
本申请提供了一种针对电网杆塔点的风速预报订正系统，包括地形复杂度最优格点筛选装置、数据重构装置、订正模型装置和综合预报装置；
地形复杂度最优格点筛选装置，用于计算代替单个杆塔点的最优格点序号，并将最优格点序号发送至数据重构装置和综合预报装置；
数据重构装置，用于采集历史观测和模式预报值，并根据接收的最优格点序号，提取对应格点的观测和预报值，并发送至订正模型装置；
订正模型装置，用于接收数据重构装置给出的格点观测和预报值，计算模式预报样本阶段的风速频率，并计算累积风频，根据累积风频设置风速阈值并构建大风订正模型和普态订正模型，将大风订正模型、普态订正模型和风速阈值发送至综合预报装置；
综合预报装置，用于采集模式预报值，根据接收的地形复杂度最优格点筛选装置给的格点序号，提取对应格点的预报值，根据接收的订正模型构建模块的风速阈值和订正模型，计算最终的格点预报订正值。
如上所述的针对电网杆塔点的风速预报订正系统，其中，地形复杂度最优格点筛选装置具体包括距离权重系数模块、下垫面权重系数模块、高度权重系数模块和综合复杂度模块；
距离权重系数模块，具体用于根据杆塔点的经纬度与周围最近的n个格点，计算反距离权重系数，得到其距离权重系数，并发送至综合地形复杂度模块；
下垫面权重系数模块，具体用于根据高分辨率的下垫面数据，得到杆塔点的下垫面类型，计算杆塔点附近的n个格点的下垫面类型，得到各个格点的下垫面权重系数，并发送至综合地形复杂度模块；
高度权重系数模块，具体用于根据杆塔点的海拔高度，以及各个格点的海拔高度，计算杆塔点附近的n个格点的高度权重系数，并发送至综合地形复杂度模块；
综合复杂度模块，具体用于根据距离权重系数模块发送的距离权重系数、下垫面权重系数模块发送的下垫面权重系数、以及高度权重系数模块发送的高度权重系数计算综合相似度。
如上所述的针对电网杆塔点的风速预报订正系统，其中，在距离权重系数模块中，以每个杆塔点的经纬度和周围的格点计算距离，计算距离的公式如下：
Figure RE-GDA0003027264540000031
其中，loni为格点i的经度，lati为格点i的纬度，glon和glat分别是杆塔点的经度和纬度，di为杆塔点和格点的距离；
然后选择和杆塔点距离最小的n个格点计算的各个格点的距离权重系数计算公式如下：
Figure RE-GDA0003027264540000041
其中，n为选取的距离杆塔点最近的格点数目，Zi为第i点的距离权重系数，di为第i个格点到杆塔点的距离。
如上所述的针对电网杆塔点的风速预报订正系统，其中，在下垫面权重系数模块中，下垫面权重系数的计算公式如下：
Figure RE-GDA0003027264540000042
其中，n为选取的距离杆塔点最近的格点数目，Di为第i点的下垫面权重系数，dg为杆塔点的下垫面类型，di为第i个格点的下垫面类型。
如上所述的针对电网杆塔点的风速预报订正系统，其中，在高度权重模块中，每个杆塔点周围n个最近格点的高度权重系数计算公式如下：
Figure RE-GDA0003027264540000043
其中，n为选取的距离杆塔点最近的格点数目，Hi为第i点的高度权重系数，hg为杆塔点的海拔高度，hi为第i个格点的海拔高度。
如上所述的针对电网杆塔点的风速预报订正系统，其中，数据重构装置具体包括数据采集模块和数据整合模块；
数据采集模块，用于采集历史观测数据和历史模式预报数据，并发送至数据整合模块；
数据整合模块，用于根据接收地形复杂度最优格点筛选装置给出的最优格点序号，提取该序号的观测值和预报值，作为杆塔点历史观测和预报值发送到订正模型装置。
如上所述的针对电网杆塔点的风速预报订正系统，其中，订正模型装置，具体用于计算累积风频达到0.9时所对应的风速值作为风速阈值，根据风速阈值分类，大于等于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建大风订正模型，小于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建普态订正模型，并将上述两个模型及风速阈值发送至综合预报装置。
如上所述的针对电网杆塔点的风速预报订正系统，其中，订正模型装置具体包括风频计算模块和订正模型构建模块；
风频计算模块，用于接收所述数据重构装置的杆塔点历史预报数据，并分级计算历史模式预报数据样本每个风速段的风速频率，并计算累积频率，计算累积频率达到阈值所对应的风速阈值，并发送风频数据到订正模型构建模块；
订正模型构建模块，用于接收所述数据重构装置的杆塔点历史观测数据和历史预报数据，并接收风频计算模块的风频数据，计算累积频率达到0.9 所对应的风速阈值，根据风速阈值对预报进行分类，构建大风订正模型和普态订正模型，并将风速阈值和两个订正模型发送到综合预报装置。
如上所述的针对电网杆塔点的风速预报订正系统，其中，综合预报装置具体包括预报值提取模块和预报值订正模块；
预报值提取模块，用于采集模式预报风速值，并接收地形复杂度最优格点装置的格点序号，提取对应杆塔点的预报值，并将所述杆塔点的预报值发送至预报值订正模块；
预报值订正模块，用于接收所述订正模型装置的模型风速阈值文件和订正模型，并接收预报值提取模块的杆塔点预报值，根据风速阈值和预报值判断调用的订正模型后，计算订正后的预报值。
如上所述的针对电网杆塔点的风速预报订正系统，其中，预报值订正模块，具体用于根据提取的杆塔点预报值和接收的风速阈值判断选取大风订正模型还是普态订正模型，将杆塔点预报值代入对应的订正模型公式计算预报订正值。
本申请实现的有益效果如下：采用本申请技术方案，能够获取到最接近真是的杆塔点风速观测信息，并建立误差订正模型有效提高风速预报的准确度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。
图1是本申请实施例一提供的一种针对电网杆塔点的风速预报订正系统示意图；
图2是2018年12月6日至2019日1月5日的逐时均值示意图；
图3是2018年12月6日至2019日1月5日的逐时绝对误差示意图；
图4是2018年12月6日20时-7日20时的逐时均值示意图；
图5是2018年12月6日20时-7日20时的逐时绝对误差示意图。
具体实施方式
下面结和本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
实施例一
本申请实施例一提供一种针对电网杆塔点的风速预报订正系统，应用于单一杆塔点和预报的单一时次，如图1所示，所述系统具体包括地形复杂度最优格点筛选装置、数据重构装置、订正模型装置和综合预报装置；以下对各装置进行详细描述：
(一)地形复杂度最优格点筛选装置，用于计算代替单个杆塔点的最优格点序号，并将最优格点序号发送至数据重构装置和综合预报装置。
具体地，地形复杂度最优格点筛选装置具体包括距离权重系数模块、下垫面权重系数模块、高度权重系数模块和综合复杂度模块；
①距离权重系数模块，具体用于根据杆塔点的经纬度与周围最近的n个格点，计算反距离权重系数，得到其距离权重系数，并发送至综合地形复杂度模块；
其中，在距离权重系数模块中，以每个杆塔点的经纬度和周围的格点计算距离，计算距离的公式如下：
Figure RE-GDA0003027264540000071
式(1)中，loni为格点i的经度，lati为格点i的纬度，glon和glat分别是杆塔点的经度和纬度，di为杆塔点和格点的距离；
然后选择和杆塔点距离最小的n个格点计算的各个格点的距离权重系数计算公式如下：
Figure RE-GDA0003027264540000072
式(2)中，n为选取的距离杆塔点最近的格点数目，Zi为第i点的距离权重系数，di为第i个格点到杆塔点的距离。距离权重系数的取值范围从0到1；当某点与杆塔点经纬度完全重合时，该点权重系数为1。
②下垫面权重系数模块，具体用于根据高分辨率的下垫面数据，得到杆塔点的下垫面类型，计算杆塔点附近的n个格点的下垫面类型，得到各个格点的下垫面权重系数，并发送至综合地形复杂度模块；
其中，在下垫面权重系数模块中，下垫面权重系数的计算公式如下：
Figure RE-GDA0003027264540000081
式(3)中，n为选取的距离杆塔点最近的格点数目，Di为第i点的下垫面权重系数，dg为杆塔点的下垫面类型，di为第i个格点的下垫面类型；下垫面权重系数的取值范围从0到1，当某点与杆塔点下垫面类型一致时，该点的下垫面权重系数为1。
③高度权重系数模块，具体用于根据杆塔点的海拔高度，以及各个格点的海拔高度，计算杆塔点附近的n个格点的高度权重系数，并发送至综合地形复杂度模块；
其中，在高度权重模块中，每个杆塔点周围n个最近格点的高度权重系数计算公式如下：
Figure RE-GDA0003027264540000082
式(4)中，n为选取的距离杆塔点最近的格点数目，Hi为第i点的高度权重系数，hg为杆塔点的海拔高度，hi为第i个格点的海拔高度；高度权重系数的取值范围从0到1，当某点与杆塔点下海拔高度相同时，该点的高度权重系数为 1。
④综合复杂度模块，具体用于根据距离权重系数模块发送的距离权重系数、下垫面权重系数模块发送的下垫面权重系数、以及高度权重系数模块发送的高度权重系数计算综合相似度；
在综合地形复杂度模块中，区域内的格点从西往东，从南往北做循环递增排序，即第一行第一列的格点序号为1，第一行第二列的格点序号为2，以此类推至最后一行最后一列的格点序号为M；
其中，综合相似度公式如下：
Ei＝zi+Di+Hi (5)
式(5)中，Ei表示格点综合相似度，Zi为距离权重系数，Di为下垫面权重系数，Hi为高度权重系数；
本申请实施例中，优选通过地形复杂度最接近法计算杆塔点最优格点序号，地形复杂度最接近算法为计算n个格点的综合相似度，综合相似度越大，该点越与杆塔点的地形复杂度越接近，因此选取综合相似度最大的格点来代表杆塔点，该格点序号即为杆塔点的最优格点序号。
(二)数据重构装置，用于采集历史观测和模式预报值，并根据接收的地形复杂度最优格点筛选装置的最优格点序号，提取对应格点的观测和预报值，并发送至订正模型装置。
具体地，数据重构装置具体包括数据采集模块和数据整合模块；
①数据采集模块，用于采集历史观测数据和历史模式预报数据，并发送至数据整合模块；
②数据整合模块，用于根据接收地形复杂度最优格点筛选装置给出的最优格点序号，提取该序号的观测值和预报值，作为杆塔点历史观测和预报值发送到订正模型装置。
例如，根据模式的预报时效T，会出现T个样本区间，即假设历史数据存在 N日，则第一个样本区间，对应模式预报的第一时次N日的值，如第一个时次对应08时，则对应N日的观测08时的值。
(三)订正模型装置，用于接收数据重构装置给出的格点观测和预报值，计算模式预报样本阶段的风速频率，并计算累积风频，根据累积风频设置风速阈值并构建大风订正模型和普态订正模型，将大风订正模型、普态订正模型和风速阈值发送至综合预报装置。
本申请实施例中，计算累积风频达到0.9时所对应的风速值作为风速阈值，根据风速阈值分类，大于等于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建大风订正模型，小于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建普态订正模型，并将上述两个模型及风速阈值发送至综合预报装置。
具体地，订正模型装置具体包括风频计算模块和订正模型构建模块；
①风频计算模块，用于接收所述数据重构装置的杆塔点历史预报数据，并分级计算历史模式预报数据样本每个风速段的风速频率，并计算累积频率，计算累积频率达到阈值所对应的风速阈值，并发送风频数据到订正模型构建模块；
其中，在风频计算模块中，对历史模式预报样本时段中风速的分级如下：
Figure RE-GDA0003027264540000101
式(6)中，i表示级别，i＝1，2，3…50，ws为模式预报风速，单位m/s，即风速小于0.5m/s时为1级，风速大于等于0.5m/s且小于1m/s时为2级，以此类推，50级时风速大于等于24.5m/s；
计算每个风速分级的风速频率公式如下：
Figure RE-GDA0003027264540000102
式(7)中，csik表示第k个样本区间内其风速为i级的次数，fik表示第k个样本区间内i级的风速频率，n为风速的分级数，n＝1，2，3…50。
②订正模型构建模块，用于接收所述数据重构装置的杆塔点历史观测数据和历史预报数据，并接收风频计算模块的风频数据，计算累积频率达到0.9所对应的风速阈值，根据风速阈值对预报进行分类，构建大风订正模型和普态订正模型，并将风速阈值和两个订正模型发送到综合预报装置；
其中，在订正模型构建模块中，根据接收的风频数据，计算累积风频数据公式如下：
Figure RE-GDA0003027264540000111
式(8)中，m为分级，Pmk为第k个样本区间的m级累积风频，fik为第k 个样本区间的i级的风速频率；根据Pmk第一次大于等于0.9时的m值，计算得模型风速阈值为m*0.5，后文中记为Ymk。
然后根据风速阈值对预报进行分类，构建的两类订正模型公式如下：
Figure RE-GDA0003027264540000112
式(9)中，x为模式预报值，a1k和b1k为在第k个历史样本时段内，当模式预报值大于等于风速阈值Ymk时，预报和观测构建的线性回归方程中的系数，x ≥Ymk时的模型公式为大风订正模型；a2k和b2k为在历史样本时段内，当模式预报值小于风速阈值Ymk时，预报和观测构建的线性回归方程中的系数，x<Ymk时的模型公式为普态订正模型；Wk(x)为第k个时次的订正预报值。
(四)综合预报装置，用于采集模式预报值，根据接收的地形复杂度最优格点筛选装置给的格点序号，提取对应格点的预报值，根据接收的订正模型构建模块的风速阈值和订正模型，计算最终的格点预报订正值。
具体地，综合预报装置具体包括预报值提取模块和预报值订正模块；
①预报值提取模块，用于采集模式预报风速值，并接收地形复杂度最优格点装置的格点序号，提取对应杆塔点的预报值，并将所述杆塔点的预报值发送至预报值订正模块；
②预报值订正模块，用于接收所述订正模型装置的模型风速阈值文件和订正模型，并接收预报值提取模块的杆塔点预报值，根据风速阈值和预报值判断调用的订正模型后，计算订正后的预报值。
其中，在预报值订正模块中，根据提取的杆塔点预报值，和接收的风速阈值 Ymk判断选取大风订正模型还是普态订正模型，并将杆塔点预报值代入对应的订正模型公式计算预报订正值。
以下结合图1所示的针对电网杆塔点的风速预报订正系统对北京区域内的电网杆塔点的风速预报进行订正的具体案例为例来说明本申请技术方案的实际应用：
具体地，先在地形复杂度最优格点筛选装置中对北京区域内的数值模式预报格点进行了排序，模式分辨率为2.5公里，北京范围内格点数为47492个；针对单个杆塔点，根据其经纬度计算出距离其最近的16个格点，经过距离权重模块、下垫面权重模块、高度权重模块的计算，最终得到综合地形复杂度最接近的最优格点，各个杆塔点的最优格点序号如下表1：
表1杆塔点的最优格点序号
杆塔点序号 格点序号
1 2111
2 4304
3 26584
4 7942
5 9111
6 11902
7 35029
8 38495
然后在数据重构装置中，采集了2018年11月1日至12月05日的总共35 天的历史观测和预报数据，选取模式预报的每日前24个预报时次，根据最优格点序号提取了对应格点的每个时次的历史观测和预报数据作为样本时段，即每个小时的35天数据均为一个样本区间。针对单个格点，总共有24个样本区间，每个样本区间包括35个风速观测数据和预报数据；
在订正模型装置中，对单个格点每个样本区间内的35个预报数据按照风速分级(如下表2所示)计算风速频率，并计算累积风速频率第一次大于等于0.9 时的风速阈值。根据风速阈值对模式该预报样本区间内的风速进行划分，得到大于等于风速阈值和小于风速阈值的两个分样本，并对两个分样本内的观测和预报分别构建线性回归方程，形成大风订正模型和普态订正模型。即单个格点有24 个风速阈值，每个风速阈值作为2个不同的订正模型的选取判别条件，并将对应杆塔点的阈值和对应的订正模型发送到了综合预报装置；
表2风速分级
级别 风速范围
1级 [0,0.5)
2级 [0.5,1)
3级 [1,1.5)
4级 [1.5,2)
5级 [2,2.5)
6级 [2.5,3)
7级 [3,3.5)
8级 [3.5,4)
9级 [4,4.5)
… …
47级 [23,23.5)
48级 [23.5,24)
49级 [24,24.5)
50级 [24.5,∞)
在综合预报装置中采集了2018年12月6日至2019日1月5日，31天的模式预报做检验时段，并根据地形复杂度筛选最优格点装置给出的序号，采集了对应格点的预报值，根据预报时次对应的格点风速阈值，调用对应的订正模型，计算了对应格点的订正值。
如14时，第一个杆塔点的风速阈值为4.5m/s，该点的订正模型见下方：
Figure RE-GDA0003027264540000141
上式中，
Figure RE-GDA0003027264540000142
为订正值，ws为模式预报值，第一个公式为普态订正模型，第二个公式为大风订正模型。
此时该点模式预报值为3.2m/s，满足第一个公式要求，则调用该点该时次的普态订正模型计算订正值为2.51m/s。其他格点和时次按照样本时段计算的对应阈值，根据模式预报值判断所采用订正模型，并计算得到最终订正值。
根据计算，检验时段内杆塔点的观测均值为3.31m/s，模式预报均值为 4.06m/s，订正均值为3.24m/s，三者的逐时变化见图2。模式预报与观测的绝对误差为1.3m/s，订正与观测的绝对误差为0.9m/s，模式预报的绝对误差和订正的绝对误差逐时变化见图3。从均值和绝对误差上表现出，订正后的值比模式本身的预报值更接近实况。
在检验时段中，提取出单日预报做对比检验结果如下：
2018年12月6日20时-12月7日20时，针对单日预报的24小时预报订正检验，观测平均值为2.98m/s，模式预报平均值为3.96m/s，订正平均值为3.17m/s，逐时变化见图4，模式预报的绝对误差为1.7m/s，订正的绝对误差为1.02m/s，逐时变化见图5。在单日预报的检验中，订正的风速值相对模式预报来说更贴近观测。
采用本申请的技术方案能够达到如下技术效果：
(1)本申请针对电网杆塔点的风速预报订正系统从网格化数值预报产品风速预报特征入手，通过考虑三种常用的替代杆塔点最优格点的筛选方法，综合距离最近、高度最相似、下垫面最相似三个方面，筛选出代替杆塔点的最优格点，这是获取最接近真实的杆塔点的风速观测信息；
(2)本申请针对电网杆塔点的风速预报订正系统对每个预报时次分别构建了订正模型，没有采用常用的一个回归模型对所有预报时次做订正的方法，可以有效抓住不同时次的预报和观测的误差和特性。
(3)本申请针对电网杆塔点的风速预报订正系统生成了大风订正模型和普态订正模型，两种模型结合使用，弥补了常规订正方法里极端事件无法体现的特性。
以上所述实施例，仅为本申请的具体实施方式，用以说明本申请的技术方案，而非对其限制，本申请的保护范围并不局限于此，尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围。都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种针对电网杆塔点的风速预报订正系统，其特征在于，包括地形复杂度最优格点筛选装置、数据重构装置、订正模型装置和综合预报装置；
地形复杂度最优格点筛选装置，用于计算代替单个杆塔点的最优格点序号，并将最优格点序号发送至数据重构装置和综合预报装置；
数据重构装置，用于采集历史观测和模式预报值，并根据接收的最优格点序号，提取对应格点的观测和预报值，并发送至订正模型装置；
订正模型装置，用于接收数据重构装置给出的格点观测和预报值，计算模式预报样本阶段的风速频率，并计算累积风频，根据累积风频设置风速阈值并构建大风订正模型和普态订正模型，将大风订正模型、普态订正模型和风速阈值发送至综合预报装置；
综合预报装置，用于采集模式预报值，根据接收的地形复杂度最优格点筛选装置给的格点序号，提取对应格点的预报值，根据接收的订正模型构建模块的风速阈值和订正模型，计算最终的格点预报订正值。
2.如权利要求1所述的针对电网杆塔点的风速预报订正系统，其特征在于，地形复杂度最优格点筛选装置具体包括距离权重系数模块、下垫面权重系数模块、高度权重系数模块和综合复杂度模块；
距离权重系数模块，具体用于根据杆塔点的经纬度与周围最近的n个格点，计算反距离权重系数，得到其距离权重系数，并发送至综合地形复杂度模块；
下垫面权重系数模块，具体用于根据高分辨率的下垫面数据，得到杆塔点的下垫面类型，计算杆塔点附近的n个格点的下垫面类型，得到各个格点的下垫面权重系数，并发送至综合地形复杂度模块；
高度权重系数模块，具体用于根据杆塔点的海拔高度，以及各个格点的海拔高度，计算杆塔点附近的n个格点的高度权重系数，并发送至综合地形复杂度模块；
综合复杂度模块，具体用于根据距离权重系数模块发送的距离权重系数、下垫面权重系数模块发送的下垫面权重系数、以及高度权重系数模块发送的高度权重系数计算综合相似度。
3.如权利要求2所述的针对电网杆塔点的风速预报订正系统，其特征在于，在距离权重系数模块中，以每个杆塔点的经纬度和周围的格点计算距离，计算距离的公式如下：
Figure 287542DEST_PATH_IMAGE001
其中，loni为格点i的经度，lati为格点i的纬度，glon和glat分别是杆塔点的经度和纬度，di为杆塔点和格点的距离；
然后选择和杆塔点距离最小的n个格点计算的各个格点的距离权重系数计算公式如下：
Figure 417172DEST_PATH_IMAGE002
其中，n为选取的距离杆塔点最近的格点数目，Zi为第i 点的距离权重系数，di为第i 个格点到杆塔点的距离。
4.如权利要求2所述的针对电网杆塔点的风速预报订正系统，其特征在于，在下垫面权重系数模块中，下垫面权重系数的计算公式如下：
Figure 231544DEST_PATH_IMAGE003
其中，n为选取的距离杆塔点最近的格点数目，Di为第i 点的下垫面权重系数，dg为杆塔点的下垫面类型，di为第i 个格点的下垫面类型。
5.如权利要求2所述的针对电网杆塔点的风速预报订正系统，其特征在于，在高度权重模块中，每个杆塔点周围n个最近格点的高度权重系数计算公式如下：
Figure 788427DEST_PATH_IMAGE004
其中，n为选取的距离杆塔点最近的格点数目，Hi为第i 点的高度权重系数，hg为杆塔点的海拔高度，hi为第i 个格点的海拔高度。
6.如权利要求1所述的针对电网杆塔点的风速预报订正系统，其特征在于，数据重构装置具体包括数据采集模块和数据整合模块；
数据采集模块，用于采集历史观测数据和历史模式预报数据，并发送至数据整合模块；
数据整合模块，用于根据接收地形复杂度最优格点筛选装置给出的最优格点序号，提取该序号的观测值和预报值，作为杆塔点历史观测和预报值发送到订正模型装置。
7.如权利要求1所述的针对电网杆塔点的风速预报订正系统，其特征在于，订正模型装置，具体用于计算累积风频达到0.9时所对应的风速值作为风速阈值，根据风速阈值分类，大于等于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建大风订正模型，小于风速阈值的预报风速样本段和对应时次的观测建立线性回归方程，从而构建普态订正模型，并将上述两个模型及风速阈值发送至综合预报装置。
8.如权利要求7所述的针对电网杆塔点的风速预报订正系统，其特征在于，订正模型装置具体包括风频计算模块和订正模型构建模块；
风频计算模块，用于接收所述数据重构装置的杆塔点历史预报数据，并分级计算历史模式预报数据样本每个风速段的风速频率，并计算累积频率，计算累积频率达到阈值所对应的风速阈值，并发送风频数据到订正模型构建模块；
订正模型构建模块，用于接收所述数据重构装置的杆塔点历史观测数据和历史预报数据，并接收风频计算模块的风频数据，计算累积频率达到0.9所对应的风速阈值，根据风速阈值对预报进行分类，构建大风订正模型和普态订正模型，并将风速阈值和两个订正模型发送到综合预报装置。
9.如权利要求1所述的针对电网杆塔点的风速预报订正系统，其特征在于，综合预报装置具体包括预报值提取模块和预报值订正模块；
预报值提取模块，用于采集模式预报风速值，并接收地形复杂度最优格点装置的格点序号，提取对应杆塔点的预报值，并将所述杆塔点的预报值发送至预报值订正模块；
预报值订正模块，用于接收所述订正模型装置的模型风速阈值文件和订正模型，并接收预报值提取模块的杆塔点预报值，根据风速阈值和预报值判断调用的订正模型后，计算订正后的预报值。
10.如权利要求9所述的针对电网杆塔点的风速预报订正系统，其特征在于，预报值订正模块，具体用于根据提取的杆塔点预报值和接收的风速阈值判断选取大风订正模型还是普态订正模型，将杆塔点预报值代入对应的订正模型公式计算预报订正值。
CN202110172682.4A 2021-02-08 2021-02-08 一种针对电网杆塔点的风速预报订正系统 Active CN113176420B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110172682.4A CN113176420B (zh) 2021-02-08 2021-02-08 一种针对电网杆塔点的风速预报订正系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110172682.4A CN113176420B (zh) 2021-02-08 2021-02-08 一种针对电网杆塔点的风速预报订正系统

Publications (2)

Publication Number Publication Date
CN113176420A true CN113176420A (zh) 2021-07-27
CN113176420B CN113176420B (zh) 2024-03-12

Family

ID=76921759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110172682.4A Active CN113176420B (zh) 2021-02-08 2021-02-08 一种针对电网杆塔点的风速预报订正系统

Country Status (1)

Country Link
CN (1) CN113176420B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330641A (zh) * 2021-11-09 2022-04-12 国网山东省电力公司应急管理中心 基于复杂地形的深度学习的短期风速订正模型的建立方法
CN115204712A (zh) * 2022-07-26 2022-10-18 中国气象局上海台风研究所(上海市气象科学研究所) 一种海上和沿海风电场选址评估方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013050417A (ja) * 2011-08-31 2013-03-14 Seiko Epson Corp 風向風速情報提供システム及び風向風速情報提供方法
CN103389388A (zh) * 2012-05-08 2013-11-13 华锐风电科技(集团)股份有限公司 风电场风速预测方法及其装置和功率预测方法及其系统
CN103514341A (zh) * 2012-06-14 2014-01-15 华锐风电科技(集团)股份有限公司 基于数值天气预报和计算流体动力学的风资源评估方法
KR101551739B1 (ko) * 2015-05-12 2015-09-10 (주)이지스 3차원 지형상에 구조물을 위치시키는 방법
CN107300853A (zh) * 2017-06-26 2017-10-27 中船第九设计研究院工程有限公司 一种分区计算空间相关性的脉动风速时程输入方法
CN107526852A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 一种核设施事故场外后果实时在线评价方法及系统
CN107688906A (zh) * 2017-09-04 2018-02-13 北京玖天气象科技有限公司 多方法融合的输电线路气象要素降尺度分析系统及方法
CN107885968A (zh) * 2017-11-06 2018-04-06 国电联合动力技术有限公司 基于周期叠加外推的风电场年平均风速预测方法及系统
CN108564263A (zh) * 2018-04-02 2018-09-21 国网安徽省电力有限公司电力科学研究院 一种用于飑线环境下输电线路杆塔受灾预测方法
CN110503064A (zh) * 2019-08-28 2019-11-26 国网湖南省电力有限公司 一种电网覆冰微地形自动识别方法及系统
CN111220986A (zh) * 2020-02-28 2020-06-02 中国民航大学 回波功率筛选与dlcd辅助的低空风切变风速估计方法
CN111880242A (zh) * 2020-07-22 2020-11-03 中国气象局气象探测中心 一种高铁沿线大风监测点布设方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013050417A (ja) * 2011-08-31 2013-03-14 Seiko Epson Corp 風向風速情報提供システム及び風向風速情報提供方法
CN103389388A (zh) * 2012-05-08 2013-11-13 华锐风电科技(集团)股份有限公司 风电场风速预测方法及其装置和功率预测方法及其系统
CN103514341A (zh) * 2012-06-14 2014-01-15 华锐风电科技(集团)股份有限公司 基于数值天气预报和计算流体动力学的风资源评估方法
KR101551739B1 (ko) * 2015-05-12 2015-09-10 (주)이지스 3차원 지형상에 구조물을 위치시키는 방법
CN107526852A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 一种核设施事故场外后果实时在线评价方法及系统
CN107300853A (zh) * 2017-06-26 2017-10-27 中船第九设计研究院工程有限公司 一种分区计算空间相关性的脉动风速时程输入方法
CN107688906A (zh) * 2017-09-04 2018-02-13 北京玖天气象科技有限公司 多方法融合的输电线路气象要素降尺度分析系统及方法
CN107885968A (zh) * 2017-11-06 2018-04-06 国电联合动力技术有限公司 基于周期叠加外推的风电场年平均风速预测方法及系统
CN108564263A (zh) * 2018-04-02 2018-09-21 国网安徽省电力有限公司电力科学研究院 一种用于飑线环境下输电线路杆塔受灾预测方法
CN110503064A (zh) * 2019-08-28 2019-11-26 国网湖南省电力有限公司 一种电网覆冰微地形自动识别方法及系统
CN111220986A (zh) * 2020-02-28 2020-06-02 中国民航大学 回波功率筛选与dlcd辅助的低空风切变风速估计方法
CN111880242A (zh) * 2020-07-22 2020-11-03 中国气象局气象探测中心 一种高铁沿线大风监测点布设方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
D.Q DANG 等: "A multi-objective optimal nonlinear control of variable speed wind turbine", 《2009 IEEE INTERNATIONAL CONFERENCE ON CONTROL AND AUTOMATION》, pages 17 - 22 *
侯慧 等: "基于空间多源异构数据的台风下输电杆塔风险评估", 《电力系统自动化》, vol. 44, no. 10, pages 127 - 134 *
王丙兰 等: "河南电网输电线路舞动的气象要素指标研究", 《气象》, vol. 43, no. 1, pages 108 - 114 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330641A (zh) * 2021-11-09 2022-04-12 国网山东省电力公司应急管理中心 基于复杂地形的深度学习的短期风速订正模型的建立方法
CN115204712A (zh) * 2022-07-26 2022-10-18 中国气象局上海台风研究所(上海市气象科学研究所) 一种海上和沿海风电场选址评估方法
CN115204712B (zh) * 2022-07-26 2023-02-03 中国气象局上海台风研究所(上海市气象科学研究所) 一种海上和沿海风电场选址评估方法

Also Published As

Publication number Publication date
CN113176420B (zh) 2024-03-12

Similar Documents

Publication Publication Date Title
CN111045117B (zh) 一种气候监测预测平台
Allen et al. A severe thunderstorm climatology for Australia and associated thunderstorm environments
CN113919231B (zh) 基于时空图神经网络的pm2.5浓度时空变化预测方法及系统
CN111897030A (zh) 一种雷暴预警系统及方法
CN106229972B (zh) 一种基于多气象源集成并分段建模的风电功率预测方法
Mandement et al. Contribution of personal weather stations to the observation of deep-convection features near the ground
CN112114384B (zh) 一种输电线覆冰发生概率预报方法
CN104182594A (zh) 一种电力系统风区图的绘制方法
CN109472283B (zh) 一种基于多重增量回归树模型的危险天气预测方法和装置
CN113176420A (zh) 一种针对电网杆塔点的风速预报订正系统
CN112886923B (zh) 一种雷电天气的光伏电站运维方法及设备
CN109740195A (zh) 一种基于气象站观测数据的极值台风风速概率分布模型及设计台风风速的评估方法
CN105974495A (zh) 利用分类拟合法预判目标区域未来平均云量的方法
CN107403004B (zh) 一种基于地形数据的遥测雨量站点可疑数值检验方法
Supari Spatiotemporal characteristics of extreme rainfall events over Java Island, Indonesia
CN114048944A (zh) 一种暴雨诱发地质灾害下应撤离人口及损毁房屋的预估方法
CN113836808A (zh) 一种基于重污染特征约束的pm2.5深度学习预测方法
Alerskans et al. Optimizing spatial quality control for a dense network of meteorological stations
CN114444020A (zh) 一种针对电网杆塔点的气温预报订正系统
Bajamgnigni Gbambie et al. Added value of alternative information in interpolated precipitation datasets for hydrology
Jia et al. An application of the LTP_DSEF model to heavy precipitation forecasts of landfalling tropical cyclones over China in 2018
Miller et al. A preliminary assessment of using spatiotemporal lightning patterns for a binary classification of thunderstorm mode
CN112948352A (zh) 一种构建大气光学湍流时空特征及概率化数据库的方法
CN115936238B (zh) 一种全域风电中长期出力预测方法、系统、设备及介质
CN112632799B (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