CN112540748B - 一种中尺度风能资源分析自动化运行系统 - Google Patents

一种中尺度风能资源分析自动化运行系统 Download PDF

Info

Publication number
CN112540748B
CN112540748B CN202011237572.3A CN202011237572A CN112540748B CN 112540748 B CN112540748 B CN 112540748B CN 202011237572 A CN202011237572 A CN 202011237572A CN 112540748 B CN112540748 B CN 112540748B
Authority
CN
China
Prior art keywords
mode
script file
script
file
wrf
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.)
Active
Application number
CN202011237572.3A
Other languages
English (en)
Other versions
CN112540748A (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 Langrun Zhitian Technology Co ltd
Huaneng Group Technology Innovation Center Co Ltd
Huaneng Renewables Corp Ltd
Original Assignee
Beijing Langrun Zhitian Technology Co ltd
Huaneng Group Technology Innovation Center Co Ltd
Huaneng Renewables Corp 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 Langrun Zhitian Technology Co ltd, Huaneng Group Technology Innovation Center Co Ltd, Huaneng Renewables Corp Ltd filed Critical Beijing Langrun Zhitian Technology Co ltd
Priority to CN202011237572.3A priority Critical patent/CN112540748B/zh
Publication of CN112540748A publication Critical patent/CN112540748A/zh
Application granted granted Critical
Publication of CN112540748B publication Critical patent/CN112540748B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F8/00Arrangements for software engineering
    • G06F8/20Software design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F8/00Arrangements for software engineering
    • G06F8/30Creation or generation of source code
    • G06F8/35Creation or generation of source code model driven

Landscapes

  • Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开一种基于Linux系统bash脚本控制的中尺度风能资源分析自动化运行系统,包括参数配置脚本文件、中尺度模式WRF前处理WPS运行脚本文件、中尺度模式WRF运行脚本、资料同化脚本文件、模式后处理资料同化脚本文件、模式产品生成LRPP脚本文件、运行控制脚本文件;避免了气象业务中常用的rocoto、ECflow、Airflow等复杂的作业流程控制软件,采用简单的bash脚本实现作业的提交、运行、监控、补算等功能;所有的参数配置都设置在参数配置脚本文件中,其余的脚本都采取模块化方法,实现任意方式运行;实现了将大格点数的区域范围分隔为若干个相对较小的格点区域进行中尺度模拟,在产品阶段进行区域拼接。

Description

一种中尺度风能资源分析自动化运行系统
技术领域
本发明涉及中尺度风能资源分析技术领域,具体涉及一种基于Linux系统bash脚本控制的中尺度风能资源分析自动化运行系统。
背景技术
气象模拟系统运行的流程控制与作业管理方式一般分为预报业务与科研两大类。其中预报业务应用软件主要包括rocoto(美国国家环境预报中心业务流程管理软件)和ECFlow(欧洲中期天气预报中心业务流程管理软件);科研的作业管理方式比较随意,有采用上述两种软件的,也有采用Apache AirFlow的,但大部分的科研人员和科研项目采用的是手动管理或者python脚本、bash等脚本来进行控制。对于大部分气象业务和科研人员来说,不采用成熟的作业流管理软件的原因是:1)这些软件需要有一定的计算机知识才能正确安装,也就是说安装比较麻烦;2)是这些软件的安装都需要管理员账户,在一些超级计算机上,气象人员只是小部分用户,很难要求该系统使用此类管理软件;3)是这些作业流管理软件需要花时间去学习掌握。
发明内容
本发明要解决的技术问题是解决上述现有技术的不足,提供一种中尺度风能资源分析自动化运行系统。
为了解决上述技术问题,本发明采用的技术方案为:一种基于Linux系统bash脚本控制的中尺度风能资源分析自动化运行系统,包括参数配置脚本文件、中尺度模式WRF前处理WPS运行脚本文件、中尺度模式WRF运行脚本、资料同化脚本文件、模式后处理资料同化脚本文件、模式产品生成LRPP脚本文件、运行控制脚本文件;
参数配置脚本文件参数包括试验运行目录、工作路径、源数据、WPS参数、WRF模式参数、GSI同化模式参数、UPP后处理参数;中尺度模式WRF前处理WPS运行脚本文件能自动获取时间、区域参数,根据参数配置脚本文件生成WRF模式的初始场、边界条件、格点同化数据中尺度模式WRF运行脚本,该脚本自动获取时间、区域参数,根据参数配置脚本文件运行WRF模式,生成模拟结果;资料同化脚本文件包括对已经处理好的观测资料进行数值资料同化质量控制、生成GSI同化所需要的preBUFR格式文件、进行GSI资料同化;模式后处理资料同化脚本文件运行UPP模块,从WRF模拟结果计算等地面高度层上的风能资源要素;模式产品生成LRPP脚本文件读取模式后处理资料同化脚本文件中的结果,生成风能资源产品;运行控制脚本文件根据计算机状态提交作业、监控作业状态、补交失败的作业。
进一步的,作业运行管理包括作业是否需要运行、是否正在运行、计算资源是否允许现在运行。
进一步的,在模式结束后,我们需要对数值模式输出进行处理,并输出等地面高度层上的风资源数据,运行步骤包括:
1)运行UPP模块,将数值模式数据转换到等地面高度层上的风资源要素grib2文件;
2)将每一个区域的grib2结果文件都处理到同样的经纬度网格点上:
3)将每个经纬度文件依次进行合并。
从上述技术方案可以看出本发明具有以下优点:避免了气象业务中常用的rocoto、ECflow、Airflow等复杂的作业流程控制软件,采用简单的bash脚本实现作业的提交、运行、监控、补算等功能;所有的参数配置都设置在参数配置脚本文件中,其余的脚本都采取模块化方法,实现任意方式运行;实现了将大格点数的区域范围分隔为若干个相对较小的格点区域进行中尺度模拟,在产品阶段进行区域拼接。
附图说明
图1为本发明的功能框图。
具体实施方式
以下结合附图对本发明的具体实施方式做具体说明。
本发明在Linux系统下,以bash语言为基础,实现风能资源中尺度再分析任务的作业运行、控制管理以及出错后的作业重新提交。如图1所示,其包括7个核心文件:(1)参数配置脚本文件,该参数文件包含中尺度模式网格、区域、投影、文件系统、计算机信息等内容;(2)中尺度模式WRF前处理WPS运行脚本文件,该文件能自动获取时间、区域参数,根据参数配置脚本文件生成WRF模式的初始场、边界条件、格点同化数据;3)中尺度模式WRF运行脚本,该脚本自动获取时间、区域参数,根据参数配置脚本文件运行WRF模式,生成模拟结果;4)资料同化脚本文件,该文件包括对已经处理好的观测资料进行数值资料同化质量控制、生成GSI同化所需要的preBUFR格式文件、进行GSI资料同化;5)模式后处理资料同化脚本文件,该脚本文件运行UPP模块,从WRF模拟结果计算等地面高度层上的风能资源要素;6)模式产品生成LRPP脚本文件,该脚本文件读取5)中的结果,生成风能资源产品;7)运行控制脚本文件,该文件根据计算机状态提交作业、监控作业状态、补交失败的作业。
本发明的自动化运行系统避免了气象业务中常用的rocoto、ECflow、Airflow等复杂的作业流程控制软件,采用简单的bash脚本实现作业的提交、运行、监控、补算等功能;所有的参数配置都设置在参数配置脚本文件中,其余的脚本都采取模块化方法,实现任意方式运行;实现了将大格点数的区域范围分隔为若干个相对较小的格点区域进行中尺度模拟,在产品阶段进行区域拼接。
该系统的核心流程包括如下几个方面。
1.参数配置:
所有的参数配置都放在一个参数配置脚本文件func_configuration.sh中,该文件中的参数包括试验运行目录、工作路径、源数据、WPS参数、WRF模式参数、GSI同化模式参数、UPP后处理参数等。示例如下:
2.分区运行:
由于全国高分辨率的需求以及计算资源的限制,很难满足我们用一个区域来进行全国的计算,主要的原因包括:1)运行全国2km格点模式运行时,格点数将达到3000x4000的量级,我们使用的超级计算机很难满足这样的资源需求;2)采用这样大的区域,需要很多节点来并行。一般来说,节点越多,并行效率越低下;3)我国地型复杂,西部有地形高度会超过7km,有些地方地形特别陡峭,两个格点之间地形高度差可以超过1000米。在这些复杂地形区域,为了保持系统运行的稳定性,需要采用短时间步长。如果全国一个区域的话,短时间步长将需要多一到五倍的计算资源。
分区的设置同样只需要在参数配置脚本文件中设置,示例如下:
max_dom=2
time_step=30
wrf_dx=(12000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 20002000)
grids_we=(457 487 601 403 601 601 781 487 601 601 781 721 361)
grids_sn=(412 415 511 631 451 661 601 661 511 571 601 721 361)
i_parent_start=(1 153 223 313 223 293 293 223 131 131 17 17 202)
j_parent_start=(1 57 21 66 95 165 262 165 117 191 233 127 165)
center=(37.5 105.0)
parent_id=(1 1 1 1 1 1 1 1 1 1 1 1 1 1 1)
parent_grid_ratio=(1 6 6 6 6 6 6 6 6 6 6 6 6 6 6)
map_proj="lambert"#map projection
truelat1=30.0
truelat2=60.0
stand_lon=110.0
当运行模式系统时,将通过区域编码来配置各模块参数,实现分区运行。如WPS模块:parent_id=${parent_id[0]},${parent_id[${ND}]},parent_grid_ratio
=${parent_grid_ratio[0]},${parent_grid_ratio[${ND}]},
i_parent_start=1,${i_parent_start[${ND}]},${i_parent_start[2]},${i_parent_start[3]},
j_parent_start=1,${j_parent_start[${ND}]},${j_parent_start[2]},${j_parent_start[3]},
s_we =1,1,1,1,
e_we =${grids_we[0]},${grids_we[${ND}]},${grids_we[2]},${grids_we[3]},
s_sn =1,1,1,1,
e_sn =${grids_sn[0]},${grids_sn[${ND}]},${grids_sn[2]},${grids_sn[3]},
3.区域合并:
在模式结束后,我们需要对数值模式输出进行处理,并输出等地面高度层上的风资源数据。运行步骤包括:1)运行UPP模块,将数值模式数据转换到等地面高度层上的风资源要素grib2文件;2)将每一个区域的grib2结果文件都处理到同样的经纬度网格点上:
3)将每个经纬度文件依次进行合并:
4.作业运行管理:
作业运行管理包括作业是否需要运行、是否正在运行、计算资源是否允许现在运行等。分别示例如下:
1)作业是否需要运行
2)是否正在运行、是否有僵尸进程等
3)计算资源是否允许现在运行

Claims (2)

1.一种基于Linux系统bash脚本控制的中尺度风能资源分析自动化运行系统,其特征在于,包括:
参数配置脚本文件,其参数包括试验运行目录、工作路径、源数据、WPS参数、WRF模式参数、GSI同化模式参数、UPP后处理参数;
中尺度模式WRF前处理WPS运行脚本文件,该文件能自动获取时间、区域参数,根据参数配置脚本文件生成WRF模式的初始场、边界条件、格点同化数据;
中尺度模式WRF运行脚本,该脚本自动获取时间、区域参数,根据参数配置脚本文件运行WRF模式,生成模拟结果;
资料同化脚本文件,该文件包括对已经处理好的观测资料进行数值资料同化质量控制、生成GSI同化所需要的preBUFR格式文件、进行GSI资料同化;
模式后处理资料同化脚本文件,该脚本文件运行UPP模块,从WRF模拟结果计算等地面高度层上的风能资源要素;
模式产品生成LRPP脚本文件,该脚本文件读取模式后处理资料同化脚本文件中的结果,生成风能资源产品;
运行控制脚本文件,该文件根据计算机状态提交作业、监控作业状态、补交失败的作业;
所述根据参数配置脚本文件运行WRF模式,生成模拟结果,以及模式后处理资料同化脚本文件运行UPP模块,从WRF模拟结果计算等地面高度层上的风能资源要素,包括:
所有的参数配置都设置在参数配置脚本文件中,其余的脚本都采取模块化方法,分区的设置在参数配置脚本文件中设置,当运行模式系统时,将通过区域编码来配置各模块参数,实现分区运行;在模式结束后运行UPP模块,将数值模式数据转换到等地面高度层上的风资源要素grib2文件并将每一个区域的grib2结果文件都处理到同样的经纬度网格点上,将每个经纬度文件依次进行合并。
2.根据权利要求1所述的基于Linux系统bash脚本控制的中尺度风能资源分析自动化运行系统,其特征在于,作业运行管理包括作业是否需要运行、是否正在运行、计算资源是否允许现在运行。
CN202011237572.3A 2020-11-09 2020-11-09 一种中尺度风能资源分析自动化运行系统 Active CN112540748B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011237572.3A CN112540748B (zh) 2020-11-09 2020-11-09 一种中尺度风能资源分析自动化运行系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011237572.3A CN112540748B (zh) 2020-11-09 2020-11-09 一种中尺度风能资源分析自动化运行系统

Publications (2)

Publication Number Publication Date
CN112540748A CN112540748A (zh) 2021-03-23
CN112540748B true CN112540748B (zh) 2024-02-27

Family

ID=75013946

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011237572.3A Active CN112540748B (zh) 2020-11-09 2020-11-09 一种中尺度风能资源分析自动化运行系统

Country Status (1)

Country Link
CN (1) CN112540748B (zh)

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310002A (zh) * 2013-06-27 2013-09-18 曙光信息产业(北京)有限公司 用于气象预报计算系统MM5的Web服务封装方法以及系统
CN103309675A (zh) * 2013-06-27 2013-09-18 曙光信息产业(北京)有限公司 用于气象预报计算系统WRF的Web服务封装方法以及系统
WO2016036806A1 (en) * 2014-09-02 2016-03-10 Vega Rolando E Method and system for hours-ahead and day-ahead forecasting of solar energy
KR101642226B1 (ko) * 2015-06-23 2016-07-22 경북대학교 산학협력단 레이더의 반사도 및 시선속도를 모의하기 위한 전자 장치 및 방법, 이를 수행하기 위한 기록 매체
CN107247837A (zh) * 2017-06-06 2017-10-13 云南电网有限责任公司 一种输电线路的杆塔信息建模方法
CN107609702A (zh) * 2017-09-15 2018-01-19 郑州云海信息技术有限公司 一种气象数据处理方法及装置
CN107784667A (zh) * 2016-08-24 2018-03-09 中国海洋大学 基于并行的全球海洋中尺度涡快速识别算法
KR20180079101A (ko) * 2016-12-30 2018-07-10 한국환경정책평가연구원 기후변화시나리오의 역학적 상세화 방법을 이용한 기후변화 고해상도 자료생상 및 표출시스템
CN108802856A (zh) * 2018-03-30 2018-11-13 南京大学 一种基于ai的源数据动态修正预报系统及其工作方法
JP2018194497A (ja) * 2017-05-19 2018-12-06 国立研究開発法人海洋研究開発機構 気象予測システム、気象予測方法、気象予測プログラム
CN109492756A (zh) * 2018-11-19 2019-03-19 中国气象局公共气象服务中心 基于深度学习的多要素导线舞动预警方法及相关装置
CN109871986A (zh) * 2019-01-25 2019-06-11 同济大学 一种用于中尺度模式cnop求解识别观测敏感区的架构及方法
CN110275224A (zh) * 2019-05-24 2019-09-24 兰州大学 精细化近地层气象要素预报系统及其预报方法
CN110427594A (zh) * 2019-07-18 2019-11-08 北京中科资环信息技术研究院 适合小型实验室的气象要素数据获取系统
CN110456356A (zh) * 2019-08-27 2019-11-15 江苏省气象探测中心 一种面向资料同化的天气雷达径向速度质量控制方法
CN110717252A (zh) * 2019-09-20 2020-01-21 哈尔滨工业大学 一种基于wrf和sedris的虚拟大气环境数据生成方法及系统
CN111414723A (zh) * 2020-04-29 2020-07-14 中国水利水电科学研究院 一种利用wrf模式分析面雨量误差的方法
CN111639054A (zh) * 2020-05-29 2020-09-08 中国人民解放军国防科技大学 一种海洋模式与资料同化的数据耦合方法、系统及介质
KR20200112645A (ko) * 2019-03-22 2020-10-05 (주)해안해양기술 실시간 검증이 가능한 고정밀 파랑 예측 시스템

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9104451B2 (en) * 2012-11-26 2015-08-11 Oracle International Corporation Dynamic communication between script and execution layers

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103309675A (zh) * 2013-06-27 2013-09-18 曙光信息产业(北京)有限公司 用于气象预报计算系统WRF的Web服务封装方法以及系统
CN103310002A (zh) * 2013-06-27 2013-09-18 曙光信息产业(北京)有限公司 用于气象预报计算系统MM5的Web服务封装方法以及系统
WO2016036806A1 (en) * 2014-09-02 2016-03-10 Vega Rolando E Method and system for hours-ahead and day-ahead forecasting of solar energy
KR101642226B1 (ko) * 2015-06-23 2016-07-22 경북대학교 산학협력단 레이더의 반사도 및 시선속도를 모의하기 위한 전자 장치 및 방법, 이를 수행하기 위한 기록 매체
CN107784667A (zh) * 2016-08-24 2018-03-09 中国海洋大学 基于并行的全球海洋中尺度涡快速识别算法
KR20180079101A (ko) * 2016-12-30 2018-07-10 한국환경정책평가연구원 기후변화시나리오의 역학적 상세화 방법을 이용한 기후변화 고해상도 자료생상 및 표출시스템
JP2018194497A (ja) * 2017-05-19 2018-12-06 国立研究開発法人海洋研究開発機構 気象予測システム、気象予測方法、気象予測プログラム
CN107247837A (zh) * 2017-06-06 2017-10-13 云南电网有限责任公司 一种输电线路的杆塔信息建模方法
CN107609702A (zh) * 2017-09-15 2018-01-19 郑州云海信息技术有限公司 一种气象数据处理方法及装置
CN108802856A (zh) * 2018-03-30 2018-11-13 南京大学 一种基于ai的源数据动态修正预报系统及其工作方法
CN109492756A (zh) * 2018-11-19 2019-03-19 中国气象局公共气象服务中心 基于深度学习的多要素导线舞动预警方法及相关装置
CN109871986A (zh) * 2019-01-25 2019-06-11 同济大学 一种用于中尺度模式cnop求解识别观测敏感区的架构及方法
KR20200112645A (ko) * 2019-03-22 2020-10-05 (주)해안해양기술 실시간 검증이 가능한 고정밀 파랑 예측 시스템
CN110275224A (zh) * 2019-05-24 2019-09-24 兰州大学 精细化近地层气象要素预报系统及其预报方法
CN110427594A (zh) * 2019-07-18 2019-11-08 北京中科资环信息技术研究院 适合小型实验室的气象要素数据获取系统
CN110456356A (zh) * 2019-08-27 2019-11-15 江苏省气象探测中心 一种面向资料同化的天气雷达径向速度质量控制方法
CN110717252A (zh) * 2019-09-20 2020-01-21 哈尔滨工业大学 一种基于wrf和sedris的虚拟大气环境数据生成方法及系统
CN111414723A (zh) * 2020-04-29 2020-07-14 中国水利水电科学研究院 一种利用wrf模式分析面雨量误差的方法
CN111639054A (zh) * 2020-05-29 2020-09-08 中国人民解放军国防科技大学 一种海洋模式与资料同化的数据耦合方法、系统及介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于SMS的华中区域中尺度数值天气预报系统设计和实现;赖安伟;王明欢;陈晓霞;李武阶;;气象科技(01);全文第52-56页,附图1-2 *
广西GRAPES中尺度模式系统的开发和建立;刘泽军;基础科学(第3期);全文 *
气象计算网格模式预报系统的建立与优化;王彬;周斌;魏敏;;计算机应用研究(11);全文 *

Also Published As

Publication number Publication date
CN112540748A (zh) 2021-03-23

Similar Documents

Publication Publication Date Title
CN109934402A (zh) 一种风电场集控中心集中风功率预测系统及其设计方法
CN113642069B (zh) 基于bim和异构系统的建筑风荷载快速迭代设计方法
CN108777637A (zh) 一种支持服务器异构的数据中心综合管理系统和方法
DE112018000125T5 (de) Quantifizieren einer kombinierten Auswirkung untereinander abhängiger unsicherer Ressourcen in einem elektrischen Stromnetz
Zhang et al. Mlmodelci: An automatic cloud platform for efficient mlaas
CN113193881A (zh) 一种基于hplc深化应用的智能化功能验证检测方法
Morshed et al. Recommending environmental knowledge as linked open data cloud using semantic machine learning
CN106933572B (zh) 一种基于llvm中间表示程序切片的度量模型
CN112540748B (zh) 一种中尺度风能资源分析自动化运行系统
CN112948353B (zh) 一种应用于DAstudio的数据分析方法、系统及存储介质
CN110032570A (zh) 一种基于b/s架构的空间数据动态更新系统
Gómez et al. Design and development of a Java-based graphical user interface to monitor/control a meteorological real-time forecasting system
CN111666368B (zh) 一种数据处理方法、装置、存储介质及终端
CN109917776B (zh) 风力发电机组的故障智能分析方法及装置
CN110717623B (zh) 融合多种天气状况的光伏发电功率预测方法及装置、设备
CN111639381A (zh) 城市建筑太阳能资源评估信息系统及其工作方法
Zhang et al. The E3SM Diagnostics Package (E3SM Diags v2. 6): A Python-based Diagnostics Package for Earth System Models Evaluation
CN111079982A (zh) 风电场的电缆路径的规划方法、系统、介质及电子设备
Souto et al. Processing mesoscale climatology in a grid environment
Lovas et al. Application of P-GRADE development environment in meteorology
Huang et al. Research on intelligent cloud test platform
CN102541564A (zh) 一种ui复用方法及装置
Ku et al. Self-learning mechanism for prediction of energy consumption and generation
CN114755746A (zh) 一种自动生成适用于闪电等多种资料同化预报的装置
Davidović et al. Implementation of the WRF-ARW prognostic model on the Grid

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