CN113239579A - 一种电网风区分布图绘制方法 - Google Patents

一种电网风区分布图绘制方法 Download PDF

Info

Publication number
CN113239579A
CN113239579A CN202110741596.0A CN202110741596A CN113239579A CN 113239579 A CN113239579 A CN 113239579A CN 202110741596 A CN202110741596 A CN 202110741596A CN 113239579 A CN113239579 A CN 113239579A
Authority
CN
China
Prior art keywords
wind speed
simulation
area
maximum
minute average
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
CN202110741596.0A
Other languages
English (en)
Other versions
CN113239579B (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.)
Zhanjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Research Institute of Southern Power Grid Co Ltd
Original Assignee
Zhanjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Research Institute of Southern Power Grid 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 Zhanjiang Power Supply Bureau of Guangdong Power Grid Co Ltd, Research Institute of Southern Power Grid Co Ltd filed Critical Zhanjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority to CN202110741596.0A priority Critical patent/CN113239579B/zh
Publication of CN113239579A publication Critical patent/CN113239579A/zh
Application granted granted Critical
Publication of CN113239579B publication Critical patent/CN113239579B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/06—Wind turbines or wind farms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明提供了一种电网风区分布图绘制方法，属于电网防风技术领域。该方法包括：获取绘制区域内各气象台站的气象资料；对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速；进行台风风场模拟，得到各模拟区域的最大10分钟平均风速和不同重现期的基本风速；进行基本风速校验，判断是否满足预设校验要求，若是，则执行下一步骤，若否，则调整模拟参数后重新模拟；根据各气象台站的风速数据和各模拟区域的风速数据进行差值计算并绘制风区分布图。本发明的方法能够增加空间插值点的数目，解决了使用现有的绘制方法绘制风区分布图导致有效空间分辨率不高的问题。

Description

一种电网风区分布图绘制方法
技术领域
本发明属于电网防风技术领域，具体涉及一种电网风区分布图绘制方法。
背景技术
台风造成电网导地线断线、倒塔等，容易引起大面积停电，严重威胁电网安全稳定运行。
台风造成电网导地线断线、倒塔等，容易引起大面积停电，严重威胁电网安全稳定运行。
为应对台风灾害，电网公司绘制了电网风区分布图，电网风区分布图规定了不同区域的不同重现期风速，是输变电设备防风改造和防风设计重要依据。
电网风区分布图关系到电网抗风能力的水平以及电网改造建设投资，这就要求电网风区分布图能够准确反映区域实际风速分布，一方面要保证电网的抗风能力达到抵御台风的要求，另一方面要避免过大的提升防风水平，会增大投资成本，造成浪费。电网风区分布图的绘制方法决定了电网风区分布图是否准确，因此，电网风区分布图的绘制方法就显得尤为重要。
目前，电网风区分布图绘制的一般流程包括：1）收集气象台站风速等基础资料；2）气象台站基本风速推算；3）空间插值与成图。现有绘制方法中，空间插值点大多依赖气象台站，由于气象台站空间分布和数量的局限，风区分布图的有效空间分辨率难以进一步提高。
发明内容
有鉴于此，本发明旨在解决使用现有的绘制方法绘制风区分布图导致有效空间分辨率不高的问题。
为了解决上述技术问题，本发明第一方面提供了一种电网风区分布图绘制方法，包括如下步骤：
获取绘制区域内各气象台站的气象资料；
基于各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速；
在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速；
基于各模拟区域的最大10分钟平均风速进行模拟区域的基本风速推算，得到各模拟区域不同重现期的基本风速；
基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速进行基本风速校验，判断是否满足预设校验要求，若是，则执行下一步骤，若否，则调整模拟参数，返回执行在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速及其后续步骤，直至满足所述预设校验要求；
根据各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速进行差值计算，并根据差值计算的结果绘制风区分布图。
进一步的，基于各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速具体包括：
从各气象台站的气象资料中获取各气象台站的年最大10分钟平均风速；
基于各气象台站的年最大10分钟平均风速利用极值I型分布函数推算各气象台站不同重现期的基本风速。
进一步的，极值I型分布函数具体为：
Figure 729767DEST_PATH_IMAGE001
式中：
Figure 326839DEST_PATH_IMAGE002
——各气象台站的年最大10分钟平均风速，
Figure 532693DEST_PATH_IMAGE003
；
Figure 728182DEST_PATH_IMAGE004
——分布的位置函数，即分布众值；
Figure 328665DEST_PATH_IMAGE005
——分布尺度函数；
分布尺度函数采用下式计算，
Figure 570291DEST_PATH_IMAGE006
分布的位置函数采用下式计算，
Figure 48676DEST_PATH_IMAGE007
式中：
Figure 364251DEST_PATH_IMAGE008
为样本均值，
Figure 964997DEST_PATH_IMAGE009
为样本标准差。
进一步的，各气象台站不同重现期的基本风速采用下式计算：
Figure 631602DEST_PATH_IMAGE010
式中：记各气象台站在重现期T的基本风速位
Figure 146634DEST_PATH_IMAGE011
，T为重现期，下标
Figure 644612DEST_PATH_IMAGE012
，
Figure 88363DEST_PATH_IMAGE013
为各气象台站编号。
进一步的，在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速具体包括：
将绘制区域分成多个模拟区域；
利用YanMeng风场模型对经过每个模拟区域的历史全部台风开展模拟，得到每个模拟区域的各历史台风模拟结果；
在每个模拟区域均设置模拟参数，模拟参数包括模拟边界设置参数、地表粗糙度设置参数以及分辨率设置参数；
基于每个模拟区域的分辨率设置参数得到每个模拟区域内的模拟子区域；
基于每个模拟区域的各历史台风模拟结果得到每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速；
将每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速以台风登陆顺序进行筛选排列，得到每个模拟区域内各模拟子区域的最大10分钟平均风速。
进一步的，基于各模拟区域的平均风速进行模拟区域的基本风速推算，得到各模拟区域的基本风速具体包括：
基于每个模拟区域内各模拟子区域的最大10分钟平均风速利用极值I型分布函数推算各模拟区域不同重现期的基本风速。
进一步的，基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速进行基本风速校验具体包括：
基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速计算模拟子区域的最大10分钟平均风速与距离模拟子区域最近的气象台站的年最大10分钟平均风速的相关性
Figure 507843DEST_PATH_IMAGE014
，计算公式如下，
Figure 390348DEST_PATH_IMAGE015
式中，
Figure 414936DEST_PATH_IMAGE016
为各模拟子区域的最大10分钟平均风速，
Figure 91905DEST_PATH_IMAGE002
为距离每个模拟子区域最近的气象台站的年最大10分钟平均风速，
Figure 966058DEST_PATH_IMAGE017
为
Figure 589937DEST_PATH_IMAGE002
与
Figure 734610DEST_PATH_IMAGE016
的协方差，
Figure 848060DEST_PATH_IMAGE018
为
Figure 976553DEST_PATH_IMAGE002
与
Figure 200861DEST_PATH_IMAGE016
的方差；
基于各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速计算二者的差值的绝对值
Figure 698576DEST_PATH_IMAGE019
，计算公式如下：
Figure 920610DEST_PATH_IMAGE020
式中，
Figure 67558DEST_PATH_IMAGE021
为各模拟区域不同重现期的基本风速，
Figure 298819DEST_PATH_IMAGE022
为距离每个模拟区域最近的气象台站的年最大10分钟平均风速。
进一步的，预设校验要求具体为：
相关性
Figure 152505DEST_PATH_IMAGE014
和差值的绝对值
Figure 607757DEST_PATH_IMAGE019
同时满足
Figure 943799DEST_PATH_IMAGE014
≥0.6和
Figure 509909DEST_PATH_IMAGE019
≤4。
进一步的，调整模拟参数具体包括：
调整模拟边界设置参数和地表粗糙度设置参数。
本发明第二方面提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如第一方面的一种电网风区分布图绘制方法。
综上，本发明提供了一种电网风区分布图绘制方法，通过在电网风区绘制区域内以气象台站为准获得其年最大10分钟平均风速并推算得到气象台站的不同重现期基本风速，再对绘制区域划分台风模拟区域，利用风场模型推算得到模拟区域的最大10分钟平均风速和不同重现期基本风速，然后结合气象站的风速推算结果以及模拟区域的风速推算结果进行校验，在满足校验要求的基础上，从而可以利用气象台站和模拟区域的风速数值综合进行差值计算，从而使得绘制分布图时在原来利用气象台站数据插值的基础上还可以使用各个模拟区域的数据进行差值计算，解决了使用现有的绘制方法绘制风区分布图导致有效空间分辨率不高的问题。采用本发明提供的绘制方法，可以增加空间插值点数目，有效提升分布图的有效空间分辨率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的一种电网风区分布图绘制方法的流程示意图；
图2为本发明实施例提供的采用一种电网风区分布图绘制方法绘制的风区分布图。
具体实施方式
为使得本发明的目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
请参阅图1，本实施例提供了一种电网风区分布图绘制方法，包括如下步骤：
S101：获取绘制区域内各气象台站的气象资料。
需要说明的是，绘制区域内各气象台站的观测资料既包括各种基础气象资料（只收集具有30年以上观测资料的气象台站的观测资料），同时也包括绘制区域内的历史台风资料，包括历史台风等级以及历史台风路径资料等。
S102：基于各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速。
可以理解的是，各气象台站的年最大10分钟平均风速即可从上一步骤中的基础气象资料中获取。
S103：在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速。
可以理解的是，进行台风风场模拟即是基于历史台风资料，利用台风风场模拟模型进行台风风场模拟。
S104：基于各模拟区域的最大10分钟平均风速进行模拟区域的基本风速推算，得到各模拟区域不同重现期的基本风速。
需要说明的是，在实际实施中，重现期具体取30年一遇、50年一遇、100年一遇。
S105：基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速进行基本风速校验，判断是否满足预设校验要求，若基本风速校验的结果满足预设校验要求，则执行下一步骤，若基本风速校验的结果不满足预设校验要求，则调整模拟参数，返回执行在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速及其后续步骤，直至满足所述预设校验要求。
S106：根据各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速进行差值计算，并根据差值计算的结果绘制风区分布图。
本实施例提供了一种电网风区分布图绘制方法，通过在电网风区绘制区域内以气象台站为准获得其年最大10分钟平均风速并推算得到气象台站的不同重现期基本风速，再对绘制区域划分台风模拟区域，利用风场模型推算得到模拟区域的最大10分钟平均风速和不同重现期基本风速，然后结合气象站的风速推算结果以及模拟区域的风速推算结果进行校验，在满足校验要求的基础上，从而可以利用气象台站和模拟区域的风速数值综合进行差值计算，从而使得绘制分布图时在原来利用气象台站数据插值的基础上还可以使用各个模拟区域的数据进行差值，解决了使用现有的绘制方法绘制风区分布图导致有效空间分辨率不高的问题。
以上是对本发明提供的一种电网风区分布图绘制方法的一个实施例进行详细的描述，以下将对本发明提供的一种电网风区分布图绘制方法的另一个实施例进行详细的描述。
S201：获取绘制区域内各气象台站的气象资料。
需要说明的是，该步骤与前一实施例的步骤S101相同，在此不再赘述。
S202：基于各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速。
需要说明的是，本步骤具体包括从各气象台站的气象资料中获取各气象台站的年最大10分钟平均风速；然后基于各气象台站的年最大10分钟平均风速利用极值I型分布函数推算各气象台站不同重现期的基本风速。
极值I型分布函数具体为：
Figure 483681DEST_PATH_IMAGE023
式中：
Figure 47518DEST_PATH_IMAGE002
——各气象台站的年最大10分钟平均风速，
Figure 169058DEST_PATH_IMAGE003
；
Figure 742121DEST_PATH_IMAGE004
——分布的位置函数，即分布众值；
Figure 898296DEST_PATH_IMAGE005
——分布尺度函数；
分布尺度函数采用下式计算，
Figure 865990DEST_PATH_IMAGE024
分布的位置函数采用下式计算，
Figure 412509DEST_PATH_IMAGE025
式中：
Figure 320422DEST_PATH_IMAGE008
为样本均值，
Figure 534366DEST_PATH_IMAGE009
为样本标准差。
由此，各气象台站重现期T的基本风速采用下式计算：
Figure 174426DEST_PATH_IMAGE010
式中：记各气象台站在重现期T的基本风速位
Figure 270558DEST_PATH_IMAGE011
，T为重现期，下标
Figure 418380DEST_PATH_IMAGE012
，
Figure 486830DEST_PATH_IMAGE013
为各气象台站编号。
S203：在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速。
需要说明的是，本步骤具体包括：
S2031：将绘制区域分成多个模拟区域，记为
Figure 625687DEST_PATH_IMAGE026
，
Figure 881219DEST_PATH_IMAGE027
，每个模拟区域的大小控制在200km*200km~600km*600km之间。
S2032：利用YanMeng风场模型对经过每个模拟区域的历史全部台风开展模拟，得到每个模拟区域的各历史台风模拟结果。
S2033：在每个模拟区域均设置模拟参数S，模拟参数S包括模拟边界设置参数S1、地表粗糙度设置参数S2以及分辨率设置参数S3，其中分辨率设置参数S3控制在10km*10km~50km*50km之间。
S2034：基于每个模拟区域的分辨率设置参数S3得到每个模拟区域内的模拟子区域，即模拟区域内模拟子区域数目为
Figure 599777DEST_PATH_IMAGE028
。
S2035：基于每个模拟区域的各历史台风模拟结果得到每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速。
可以理解的是，由于对每个模拟区域均划分为了若干个模拟子区域，因此根据每个模拟区域的历史台风模拟结果，就可以获得每个历史台风在每个模拟子区域内的最大10分钟平均风速。
S2036：将每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速以台风登陆顺序进行筛选排列，得到每个模拟区域内各模拟子区域的最大10分钟平均风速。
需要说明的是，每个模拟区域内历史台风的资料各不相同，有些区域同一年可能发生多起台风，而有些区域某一年份并未发生台风。故以台风登陆顺序进行筛选排列的过程中，当年份发生多次台风的取当年各台风模拟计算后得到的多个最大10分钟平均风速中的最大值，当年份没有台风的取相近3年的最大10分钟平均风速的平均值，从而得到每个模拟子区域最大10min平均风速样本，记为
Figure 585050DEST_PATH_IMAGE016
，
Figure 65448DEST_PATH_IMAGE029
。
S204：基于各模拟区域的最大10分钟平均风速进行模拟区域的基本风速推算，得到各模拟区域不同重现期的基本风速。
需要说明的是，热带气旋影响下的风速服从极值I型分布。同样的基于每个模拟区域内各模拟子区域的最大10分钟平均风速利用极值I型分布函数推算各模拟区域不同重现期的基本风速，即按步骤202得到
Figure 870593DEST_PATH_IMAGE016
的样本均值
Figure 861683DEST_PATH_IMAGE008
，样本标准差
Figure 904725DEST_PATH_IMAGE009
，计算得到模拟区域重现期T的基本风速
Figure 385385DEST_PATH_IMAGE030
，T为重现期（取30年一遇、50年一遇、100年一遇），下标
Figure 615509DEST_PATH_IMAGE026
，
Figure 643246DEST_PATH_IMAGE027
为模拟区域。
S205：基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速进行基本风速校验，判断是否满足预设校验要求，若基本风速校验的结果满足预设校验要求，则执行下一步骤，若基本风速校验的结果不满足预设校验要求，则调整模拟参数，返回执行在绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速及其后续步骤，直至满足所述预设校验要求。
需要说明的是，基本风速校验包括
1）基于各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速计算模拟子区域的最大10分钟平均风速与距离模拟子区域最近的气象台站的年最大10分钟平均风速的相关性
Figure 603111DEST_PATH_IMAGE014
，计算公式如下，
Figure 192356DEST_PATH_IMAGE031
式中，
Figure 972093DEST_PATH_IMAGE016
为各模拟子区域的最大10分钟平均风速，
Figure 304985DEST_PATH_IMAGE002
为距离每个模拟子区域最近的气象台站的年最大10分钟平均风速，
Figure 322620DEST_PATH_IMAGE032
为
Figure 145082DEST_PATH_IMAGE002
与
Figure 582755DEST_PATH_IMAGE016
的协方差，
Figure 781655DEST_PATH_IMAGE033
为
Figure 388217DEST_PATH_IMAGE002
与
Figure 584843DEST_PATH_IMAGE016
的方差；
2）基于各气象台站不同重现期的基本风速和各模拟区域不同重现期的基本风速计算二者的差值的绝对值
Figure 339172DEST_PATH_IMAGE019
，计算公式如下：
Figure 450085DEST_PATH_IMAGE034
式中，
Figure 973471DEST_PATH_IMAGE030
为各模拟区域不同重现期的基本风速，
Figure 75419DEST_PATH_IMAGE011
为距离每个模拟区域最近的气象台站的年最大10分钟平均风速。
当相关性
Figure 254727DEST_PATH_IMAGE014
和差值的绝对值
Figure 795430DEST_PATH_IMAGE019
同时满足
Figure 376584DEST_PATH_IMAGE014
≥0.6和
Figure 147969DEST_PATH_IMAGE019
≤4时，则认为满足预设校验要求，否则就需要调整模拟边界设置参数S1和地表粗糙度设置参数S2后，重新进行台风模拟。
S206：根据各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速进行差值计算，并根据差值计算的结果绘制风区分布图。
需要说明的是，需要采用风区分布图绘制区域的全部
Figure 611311DEST_PATH_IMAGE002
与
Figure 893388DEST_PATH_IMAGE016
进行差值计算，绘制风区分布图，即在基本风速校验时要求所有的模拟子区域和各气象站的数据均要满足校验要求。
请参阅图2，图中线条为原风区分布图等值线（由气象台站数据插值而成），分别为35m/s 、37m/s、 39m/s。
底图为校验合格的模拟结果，通过对比可以发现，从下至上（从海面到内陆）二者风速均逐渐减小，在大部分区域，二者风速值基本一致，说明模拟结果和气象台站观测数据大趋势是一致的，但也可以发现，通过模拟，进一步提升了风区图绘制原始数据的空间分辨率，采用气象台站数据加模拟数据的方法可以进一步风区分布图的准确性。
本实施例提供的一种电网风区分布图绘制方法，通过对模拟区域根据分辨率设置进一步划分为多个模拟子区域，然后采用极值I型分布函数作为概率模型，推断各气象台站和模拟子区域的不同重现期基本风速，然后分别对模拟子区域的最大平均风速数据和距离该模拟子区域最近的气象台站的最大平均风速数据以及模拟区域的不同重现期的基本风速数据和距离该模拟区域最近的气象台站的不同重现期的基本风速数据进行风速校验，从而增加了空间插值时模拟数据的数量，并且通过校验使得最终空间插值的数据有效，即实现了通过增加模拟仿真结果，增加空间插值点数目，有效提升分布图的有效空间分辨率。
以上是对本发明提供的一种电网风区分布图绘制方法的实施例进行详细的描述，以下将对本发明提供的一种计算机可读存储介质的实施例进行详细的描述。
一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如前述实施例的一种电网风区分布图绘制方法。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种电网风区分布图绘制方法，其特征在于，包括如下步骤：
获取绘制区域内各气象台站的气象资料；
基于所述各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速；
在所述绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速；
基于所述各模拟区域的最大10分钟平均风速进行模拟区域的基本风速推算，得到各模拟区域不同重现期的基本风速；
基于所述各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及所述各气象台站不同重现期的基本风速和所述各模拟区域不同重现期的基本风速进行基本风速校验，判断是否满足预设校验要求，若是，则执行下一步骤，若否，则调整模拟参数，返回执行在所述绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速及其后续步骤，直至满足所述预设校验要求；
根据所述各气象台站的年最大10分钟平均风速和所述各模拟区域的最大10分钟平均风速进行差值计算，并根据所述差值计算的结果绘制风区分布图。
2.根据权利要求1所述的一种电网风区分布图绘制方法，其特征在于，所述基于所述各气象台站的气象资料对各气象台站的年最大10分钟平均风速进行基本风速推算，得到各气象台站不同重现期的基本风速具体包括：
从所述各气象台站的气象资料中获取各气象台站的年最大10分钟平均风速；
基于所述各气象台站的年最大10分钟平均风速利用极值I型分布函数推算各气象台站不同重现期的基本风速。
3.根据权利要求2所述的一种电网风区分布图绘制方法，其特征在于，所述极值I型分布函数具体为：
Figure 62609DEST_PATH_IMAGE001
式中：
Figure 763586DEST_PATH_IMAGE002
——各气象台站的年最大10分钟平均风速，
Figure 47937DEST_PATH_IMAGE003
；
Figure 398147DEST_PATH_IMAGE004
——分布的位置函数，即分布众值；
Figure 426146DEST_PATH_IMAGE005
——分布尺度函数；
所述分布尺度函数采用下式计算，
Figure 342149DEST_PATH_IMAGE006
所述分布的位置函数采用下式计算，
Figure 672768DEST_PATH_IMAGE007
式中：
Figure 8809DEST_PATH_IMAGE008
为样本均值，
Figure 512602DEST_PATH_IMAGE009
为样本标准差。
4.根据权利要求3所述的一种电网风区分布图绘制方法，其特征在于，所述各气象台站不同重现期的基本风速采用下式计算：
Figure 751954DEST_PATH_IMAGE010
式中：记各气象台站在重现期T的基本风速位
Figure 846949DEST_PATH_IMAGE011
，T为重现期，下标
Figure 873548DEST_PATH_IMAGE012
，
Figure 243350DEST_PATH_IMAGE013
为各气象台站编号。
5.根据权利要求4所述的一种电网风区分布图绘制方法，其特征在于，所述在所述绘制区域内进行台风风场模拟，得到各模拟区域的最大10分钟平均风速具体包括：
将所述绘制区域分成多个模拟区域；
利用YanMeng风场模型对经过每个所述模拟区域的历史全部台风开展模拟，得到每个模拟区域的各历史台风模拟结果；
在每个所述模拟区域均设置所述模拟参数，所述模拟参数包括模拟边界设置参数、地表粗糙度设置参数以及分辨率设置参数；
基于每个所述模拟区域的分辨率设置参数得到每个所述模拟区域内的模拟子区域；
基于所述每个模拟区域的各历史台风模拟结果得到每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速；
将所述每个模拟区域内各历史台风在各模拟子区域的最大10分钟平均风速以台风登陆顺序进行筛选排列，得到所述每个模拟区域内各模拟子区域的最大10分钟平均风速。
6.根据权利要求5所述的一种电网风区分布图绘制方法，其特征在于，所述基于所述各模拟区域的平均风速进行模拟区域的基本风速推算，得到各模拟区域的基本风速具体包括：
基于所述每个模拟区域内各模拟子区域的最大10分钟平均风速利用极值I型分布函数推算各模拟区域不同重现期的基本风速。
7.根据权利要求1所述的一种电网风区分布图绘制方法，其特征在于，所述基于所述各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速以及所述各气象台站不同重现期的基本风速和所述各模拟区域不同重现期的基本风速进行基本风速校验具体包括：
基于所述各气象台站的年最大10分钟平均风速和各模拟区域的最大10分钟平均风速计算模拟子区域的最大10分钟平均风速与距离所述模拟子区域最近的气象台站的年最大10分钟平均风速的相关性
Figure 71629DEST_PATH_IMAGE014
，计算公式如下，
Figure 603104DEST_PATH_IMAGE015
式中，
Figure 211940DEST_PATH_IMAGE016
为各模拟子区域的最大10分钟平均风速，
Figure 526378DEST_PATH_IMAGE002
为距离每个模拟子区域最近的气象台站的年最大10分钟平均风速，
Figure 537059DEST_PATH_IMAGE017
为
Figure 941233DEST_PATH_IMAGE002
与
Figure 975048DEST_PATH_IMAGE016
的协方差，
Figure 624336DEST_PATH_IMAGE018
为
Figure 489523DEST_PATH_IMAGE002
与
Figure 362801DEST_PATH_IMAGE016
的方差；
基于所述各气象台站不同重现期的基本风速和所述各模拟区域不同重现期的基本风速计算二者的差值的绝对值
Figure 585710DEST_PATH_IMAGE019
，计算公式如下：
Figure 835426DEST_PATH_IMAGE020
式中，
Figure 758382DEST_PATH_IMAGE021
为各模拟区域不同重现期的基本风速，
Figure 802562DEST_PATH_IMAGE011
为距离每个模拟区域最近的气象台站的年最大10分钟平均风速。
8.根据权利要求7所述的一种电网风区分布图绘制方法，其特征在于，所述预设校验要求具体为：
所述相关性
Figure 873286DEST_PATH_IMAGE014
和所述差值的绝对值
Figure 67638DEST_PATH_IMAGE019
同时满足
Figure 343636DEST_PATH_IMAGE014
≥0.6和
Figure 824296DEST_PATH_IMAGE019
≤4。
9.根据权利要求5所述的一种电网风区分布图绘制方法，其特征在于，所述调整模拟参数具体包括：
调整所述模拟边界设置参数和所述地表粗糙度设置参数。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-9中任一项所述的一种电网风区分布图绘制方法。
CN202110741596.0A 2021-07-01 2021-07-01 一种电网风区分布图绘制方法 Active CN113239579B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110741596.0A CN113239579B (zh) 2021-07-01 2021-07-01 一种电网风区分布图绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110741596.0A CN113239579B (zh) 2021-07-01 2021-07-01 一种电网风区分布图绘制方法

Publications (2)

Publication Number Publication Date
CN113239579A true CN113239579A (zh) 2021-08-10
CN113239579B CN113239579B (zh) 2023-03-17

Family

ID=77141264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110741596.0A Active CN113239579B (zh) 2021-07-01 2021-07-01 一种电网风区分布图绘制方法

Country Status (1)

Country Link
CN (1) CN113239579B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003095010A2 (en) * 2002-05-10 2003-11-20 Oriel Therapeutics, Inc. Dry powder inhalers for use with piezoelectric polymer-driven delivery means, and associated blister package comprising a piezoelectric polymer material
US20090132176A1 (en) * 2002-12-09 2009-05-21 Verisae, Inc. Method and system for tracking and managing destruction, reconstitution, or reclamation of regulated substances
CN102128725A (zh) * 2010-12-02 2011-07-20 李惠 用于大跨度空间结构健康监测与安全预警的方法
CN102221448A (zh) * 2011-05-20 2011-10-19 大连海事大学 一种确定海上搜救区域的方法、***及搜救模拟器
US20120136496A1 (en) * 2010-11-30 2012-05-31 General Electric Company System and method for estimating demand response in electric power systems
CN102593828A (zh) * 2012-02-28 2012-07-18 重庆大学 一种含wtg的配电网可靠性评估方法
EP2924598A1 (en) * 2014-03-28 2015-09-30 Airbus Operations GmbH A method for determining a structural response of a flow body to an atmospheric disturbance
CN107420269A (zh) * 2016-05-23 2017-12-01 远景能源（江苏）有限公司 识别转子平面上的风力分布模式的方法以及实现该方法的风力涡轮机

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003095010A2 (en) * 2002-05-10 2003-11-20 Oriel Therapeutics, Inc. Dry powder inhalers for use with piezoelectric polymer-driven delivery means, and associated blister package comprising a piezoelectric polymer material
US20090132176A1 (en) * 2002-12-09 2009-05-21 Verisae, Inc. Method and system for tracking and managing destruction, reconstitution, or reclamation of regulated substances
US20120136496A1 (en) * 2010-11-30 2012-05-31 General Electric Company System and method for estimating demand response in electric power systems
CN102128725A (zh) * 2010-12-02 2011-07-20 李惠 用于大跨度空间结构健康监测与安全预警的方法
CN102221448A (zh) * 2011-05-20 2011-10-19 大连海事大学 一种确定海上搜救区域的方法、***及搜救模拟器
CN102593828A (zh) * 2012-02-28 2012-07-18 重庆大学 一种含wtg的配电网可靠性评估方法
EP2924598A1 (en) * 2014-03-28 2015-09-30 Airbus Operations GmbH A method for determining a structural response of a flow body to an atmospheric disturbance
CN107420269A (zh) * 2016-05-23 2017-12-01 远景能源（江苏）有限公司 识别转子平面上的风力分布模式的方法以及实现该方法的风力涡轮机

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MING YING 等: "An Overview of the China Meteorological Administration Tropical Cyclone Database", 《JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY》 *
吴新桥 等: "基于台风风场模型的广东沿海电网风区图优化方法", 《南方电网技术》 *

Also Published As

Publication number Publication date
CN113239579B (zh) 2023-03-17

Similar Documents

Publication Publication Date Title
AU2017352549B2 (en) Method and device for calculating power generation of wind farm
CN110598290B (zh) 考虑气候变化的流域未来水电发电能力预测方法和***
CN105512799B (zh) 一种基于海量在线历史数据的电力***暂态稳定评估方法
CN112699623B (zh) 基于非结构网格规则化重构技术的高精度热流计算方法
US20170338802A1 (en) Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation
CN103870845A (zh) 点云聚类去噪过程中新的k值优化方法
CN103871102A (zh) 一种基于高程点和道路轮廓面的道路三维精细建模方法
CN110717619A (zh) 一种自下而上配电网多尺度时空负荷预测方法及***
CN113435630B (zh) 一种产流模式自适应的流域水文预报方法及***
KR20110016523A (ko) Ｇｉｓ 자료를 활용한 풍속할증 계수 산정 시스템, 방법, 및 기록매체
CN113792430A (zh) 复杂地形风电场风切变指数计算方法及***
CN116401327A (zh) 无资料地区中小流域设计暴雨洪水计算辅助***
CN112700349A (zh) 测风塔选址方法和装置
CN113239579B (zh) 一种电网风区分布图绘制方法
CN107944466A (zh) 一种基于分段思想的降雨偏差纠正方法
CN110688805A (zh) 一种汽车外流场流动数值的分析方法
KR101286909B1 (ko) 수치모델에 접합 가능한 유전알고리즘 인터페이스 시스템
CN113516303A (zh) 一种提取多模式最优信息进行降水集合预测的方法、装置和介质
CN103793339B (zh) 基于内存访问堆栈距离的数据Cache性能探索方法
CN113434495A (zh) 一种基于ArcGIS的中尺度风速数据订正方法及***
CN110472260B (zh) 一种基于涡旋参数的流域地形截留区域确定方法及***
CN112365469B (zh) 一种基于Otsu阈值的二维波浪谱分割方法
CN117595260B (zh) 一种基于5g短切片的电力设备远程监控方法及***
CN111625877B (zh) 一种基于cad等高线提高大气污染预测准确性的方法
CN116306026B (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