CN115146562A - 一种基于中尺度气象要素边界速度场赋值方法 - Google Patents
一种基于中尺度气象要素边界速度场赋值方法 Download PDFInfo
- Publication number
- CN115146562A CN115146562A CN202210824461.5A CN202210824461A CN115146562A CN 115146562 A CN115146562 A CN 115146562A CN 202210824461 A CN202210824461 A CN 202210824461A CN 115146562 A CN115146562 A CN 115146562A
- Authority
- CN
- China
- Prior art keywords
- mesoscale
- boundary
- height
- wind speed
- grid
- 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.)
- Pending
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
本发明名称为:一种基于中尺度气象要素边界速度场赋值方法。技术领域涉及风资源评估领域,具体而言,涉及一种基于中尺度气象要素边界速度场赋值方法。包括:获取中尺度模式输出数据计算得到的中尺度气象要素数据;根据中尺度气象要素建立边界速度参考场uc、摩擦风速参考场us;使用边界速度参考场对近地层以上边界速度场进行赋值;使用中尺度气象要素数据对近地层以下边界速度场进行赋值。发明解决了在风资源评估过程中为模拟区域在设定边界条件时没有实测数据基础,完全用主观设定的环境风况构造模拟区域的边界条件的问题,采用与实际发生过的历史风况相近的气象要素进行边界条件赋值,提高了模拟计算结果的准确性。
Description
技术领域
本发明涉及风资源评估领域,具体而言,涉及一种基于中尺度气象要素边界速度场赋值方法。
背景技术
利用数值模拟技术进行风资源评估是近年来兴起的新的技术方向,在基于计算流体力学(CFD)数值模拟的风资源评估过程中,需要设置合适的入口边界条件,尤其是速度场的边界初始值,合理的初始边界速度场能够提高模拟计算收敛的速度,提高模拟计算结果的准确性。但是在目前主流的CFD风资源评估计算软件中,都是简单地使用风速廓线进行边界速度场赋值,或者是用主观设定的环境风况构造模拟区域的边界条件,即边界条件的设定没有实测数据基础,这就使得计算结果不能完全反映真实的实际风资源情况。
本发明实现了一种基于中尺度气象要素边界速度场赋值方法,使用经过再分析处理的中尺度输出数据计算得到的中尺度要素数据,采用大气边界层理论以及插值技术生成边界条件,以此对边界速度场进行赋值。用基于中尺度气象要素数据对模拟区域的边界速度场进行赋值,解决了在设定边界条件时没有实测数据基础,完全用主观设定的环境风况构造模拟区域的边界条件的问题。特别地,当所述赋值方法用于风力发电风资源评估时,在模拟区域范围的边界处,采用与实际发生过的历史风况相近的气象要素进行边界条件赋值,提高了模拟计算结果的准确性。
发明内容
本发明实施例提供了一种基于中尺度气象要素边界速度场赋值方法。特别地,还包括如何将所述方法用于对模拟区域边界速度场进行赋值,至少解决了在设定边界条件时没有实测数据基础,完全用主观设定的环境风况构造模拟区域的边界条件的问题。特别地,当所述赋值方法用于风力发电风资源评估时,在模拟区域范围的边界处,采用与实际发生过的历史风况相近的气象要素进行边界条件赋值,提高了模拟计算结果的准确性。
根据发明实施例的一个方面,提供一种基于中尺度气象要素边界速度场赋值方法,该方法包括:获取中尺度模式输出数据计算得到的中尺度气象要素数据;根据中尺度气象要素建立边界速度参考场uc、摩擦风速参考场us;使用边界速度参考场对近地层以上边界速度场进行赋值;使用中尺度气象要素数据对近地层以下边界速度场进行赋值。
进一步地,获取中尺度模式输出数据计算得到的中尺度气象要素数据,包括:当所述方法用于风力发电风资源评估时,中尺度气象要素数据包括以下内容:各个中尺度格点的风速分量ux,uy,uz、地面温度T、摩擦风速u*、铅直热通量HFX。
中尺度气象要素数据是按照中尺度位置格点进行组织的,每个位置格点至少包含15个高层以上的数据。
加载所有边界面网格,获取各个边界面网格的坐标、距地面高度h、地面的海拔高度zg;
找到与各个边界面网格水平距离小于R的所有中尺度格点,为每个边界面网格生成一个临时对象,每个临时对象单元包括:水平距离小于R的所有中尺度格点的下标、中尺度格点到网格的距离di、中尺度格点风速分量ux,uy,uz、风速分量在垂直方向的梯度
进一步地,根据中尺度气象要素建立边界速度参考场,包括:当所述方法用于风资源评估时,包括:
建立边界速度参考场时会同时考虑地形的影响和海拔高度的影响;
在计算底层高度的速度参考场时,充分考虑地形及地表高度h对风速的影响;在计算高层的速度参考场时将减小地表高度的影响,主要考虑海拔高度z对风速的影响;这两个高度层的分界线距地面高度和海拔高度分别为图2中hs和zs;
基于变换后的高度h′建立参考场uc;对于每个边界面网格,由上至下遍历临时对象,找到第一个低于边界面网格高度h′的中尺度层级,获取该高度层中尺度格点对象单元的高度h0、风速分量ux0,uy0,uz0、风速分量在垂直方向的梯度通过加权计算得到uc的三个分量,保存到uc参考场;
本发明中使用的加权方法为反距离平方加权算法:
其中N为与本网格及边界面水平距离小于R的所有中尺度格点个数;di为第i个中尺度格点到本网格的距离;uxi,uyi,uzi为第i个中尺度格点在h′高度的插值:
使用中尺度气象要素摩擦风速u*建立摩擦风速参考场us:
其中u*i为第i个中尺度格点的摩擦风速。。
进一步地,使用边界速度参考场对近地层以上边界速度场进行赋值,包括:
设定图2中的近地层高度bs,在近地层高度以上及以下采用不同的方法对边界速度场进行赋值;
使用参考场uc中对应的边界面网格数据对近地层以上边界速度场进行赋值。
进一步地,使用中尺度气象要素数据对近地层以下边界速度场进行赋值,包括:
使用中尺度气象要素数据和摩擦风速参考场us建立近地层非中性风廓线:
其中:
其中:
其中z0为地表粗糙度值,设为固定值,或者通过地形文件设置;u*为摩擦风速参考场us中对应边界面网格的数据值;L为莫宁-奥布霍夫长度,计算如下:
其中HFX为中尺度气象要素铅直热通量;T为中尺度气象要素地面温度;cp为空气的等压比热容,量纲为J/(kg*K);g为重力加速度,量纲为m/s2;κ为冯·卡门系数;ρ为空气密度值。。
进一步地,在对近地层以下边界速度场进行赋值时,会通过修正近地层风廓线的方式对接近地层边界,使得近地层上下平滑过渡,并使用修正后的风速值对近地层以下边界速度场进行赋值,包括:
图2中修正后的近地层风廓线实现方法为:
其中us(bs)为边界面网格位置在近地层边界处通过中尺度格点数据反距离平方加权得到的风速;un(bs)为非中性风廓线相同高度的风速;un(h)为近地层非中性风廓线在h高度的风速值;u(h)为h高度修正后的风速值。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是基于中尺度气象要素边界速度场赋值方法的流程图;
图2是根据本发明实施例的参数定义示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及说明书附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
术语解释
模拟区域,是表示在风资源评估分析过程中,采用CFD技术进行风资源计算的区域。
地表粗糙度,是表示地球表面粗糙程度并具有长度量纲的特征参数。
下面结合图1基于中尺度气象要素边界速度场赋值方法的流程图对本发明的具体实施方式作进一步的详细说明。
步骤101,获取中尺度模式输出数据计算得到的中尺度气象要素数据。中尺度气象要素数据包括以下内容:各个中尺度格点的风速分量ux,uy,uz、地面温度T、摩擦风速u*、铅直热通量HFX。可选地,每个位置格点至少包含15个高层以上的数据。
加载所有边界面网格,获取各个边界面网格的坐标、距地面高度h、地面的海拔高度zg。
找到与各个边界面网格水平距离小于R的所有中尺度格点,为每个边界面网格生成一个临时对象,每个临时对象单元包括:水平距离小于R的所有中尺度格点的下标、中尺度格点到网格的距离di、中尺度格点风速分量ux,uy,uz、风速分量在垂直方向的梯度可选地,R设置为6000米。
步骤102,根据中尺度气象要素建立边界速度参考场。在计算底层高度的速度参考场时,充分考虑地形及地表高度h对风速的影响;在计算高层的速度参考场时将减小地表高度的影响,主要考虑海拔高度z对风速的影响;这两个高度层的分界线距地面高度和海拔高度分别为图2中的hs和zs。可选地,该分界线可以设置为中尺度倒数第二层,获取高度层的平均距地面高度和海拔高度设置为hs和zs。
基于变换后的高度h′建立参考场uc;对于每个边界面网格,由上至下遍历临时对象,找到第一个低于边界面网格高度h′的中尺度层级,获取该高度层中尺度格点对象单元的高度h0、风速分量ux0,uy0,uz0、风速分量在垂直方向的梯度通过加权计算得到uc的三个分量,保存到uc参考场;
本发明中使用的加权方法为反距离平方加权算法:
其中N为与本网格及边界面水平距离小于R的所有中尺度格点个数;di为第i个中尺度格点到本网格的距离;uxi,uyi,uzi为第i个中尺度格点在h′高度的插值:
使用中尺度气象要素摩擦风速u*建立摩擦风速参考场us:
其中u*i为第i个中尺度格点的摩擦风速。
步骤103,使用边界速度参考场对近地层以上边界速度场进行赋值。首先设定近地层高度bs,在近地层高度以上及以下采用不同的方法对边界速度场进行赋值。可选地,近地层高度bs可以设置为150米、200米或300米。使用参考场uc中对应的边界面网格数据对近地层以上的边界速度场进行赋值。
步骤104,使用中尺度气象要素数据对近地层以下边界速度场进行赋值。首先使用中尺度气象要素数据和摩擦风速参考场us建立近地层非中性风廓线:
其中:
其中:
其中z0为地表粗糙度值;u*为摩擦风速参考场us中对应边界面网格的数据值;L为莫宁-奥布霍夫长度,计算如下:
其中HFX为中尺度气象要素铅直热通量;T为中尺度气象要素地面温度;cp为空气的等压比热容,量纲为J/(kg*K);g为重力加速度,量纲为m/s2;κ为冯·卡门系数;ρ为空气密度值。可选地,空气的等压比热容cp可以设为经验值1000。可选地,重力加速度g可以设为9.8。可选地,冯·卡门系数κ可设为0.4。可选地,地表粗糙度值z0可以设为固定值,如0.1、0.01,或者通过读取地理粗糙度文件设置。可选地,空气密度ρ值可以通过公式计算得到,其中P为大气压力(单位为Pa),R为气体常数(单位为J/(kg*K)),取值为287.05,T为温度(单位为K)。可选地,空气密度ρ值可以设置为标准空气密度1.293kg/m3。
步骤105,修正近地层风廓线以对接近地层边界,使得近地层上下平滑过渡,并使用修正后的风速值对近地层以下边界速度场进行赋值。图2中修正后的近地层风廓线实现方法为:
其中us(bs)为边界面网格位置在近地层边界处通过中尺度格点数据反距离平方加权得到的风速;un(bs)为非中性风廓线相同高度的风速;un(h)为近地层非中性风廓线在h高度的风速值;u(h)为h高度修正后的风速值。
使用修正后的风速值对近地层以下边界速度场进行赋值,完成边界速度场赋值操作。
Claims (6)
1.一种基于中尺度气象要素边界速度场赋值方法,其特征在于,包括:
获取中尺度模式输出数据计算得到的中尺度气象要素数据;
根据中尺度气象要素建立边界速度参考场uc、摩擦风速参考场us;
使用边界速度参考场对近地层以上边界速度场进行赋值;
使用中尺度气象要素数据对近地层以下边界速度场进行赋值。
2.按照权利要求1所述的方法,其特征在于:获取中尺度模式输出数据计算得到的中尺度气象要素数据,包括:
中尺度气象要素数据包括各个中尺度格点的风速分量ux,uy,uz、地面温度T、摩擦风速u*、铅直热通量HFX;
中尺度气象要素数据包括各个中尺度水平位置格点进行组织的,每个水平位置格点至少包含15个高层以上的数据;
加载所有边界面网格,获取各个边界面网格的坐标、距地面高度h、地面的海拔高度zg;
3.按照权利要求1所述的方法,其特征在于:根据中尺度气象要素建立边界速度参考场,包括:
建立边界速度参考场时会同时考虑地形的影响和海拔高度的影响;
在计算底层高度的速度参考场时,充分考虑地形及地表高度h对风速的影响;在计算高层的速度参考场时将减小地表高度的影响,主要考虑海拔高度z对风速的影响;具体而言,这两个高度层的分界线距地面高度和海拔高度分别为hs和zs;
基于变换后的高度h′建立参考场uc;对于每个边界面网格,由上至下遍历临时对象,找到第一个低于边界面网格高度h′的中尺度层级,获取该高度层中尺度格点对象单元的高度h0、风速分量ux0,uy0,uz0、风速分量在垂直方向的梯度通过加权计算得到uc的三个分量,保存到uc参考场;
本发明中使用的加权方法为反距离平方加权算法:
其中N为与本网格及边界面水平距离小于R的所有中尺度格点个数;di为第i个中尺度格点到本网格的距离;uxi,uyi,uzi为第i个中尺度格点在h′高度的插值:
使用中尺度气象要素摩擦风速u*建立摩擦风速参考场us:
其中u*i为第i个中尺度格点的摩擦风速。
4.按照权利要求1所述的方法,其特征在于:使用边界速度参考场对近地层以上边界速度场进行赋值,包括:
设定近地层高度bs;
使用参考场uc中对应的边界面网格数据对近地层以上边界速度场进行赋值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210824461.5A CN115146562A (zh) | 2022-07-14 | 2022-07-14 | 一种基于中尺度气象要素边界速度场赋值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210824461.5A CN115146562A (zh) | 2022-07-14 | 2022-07-14 | 一种基于中尺度气象要素边界速度场赋值方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115146562A true CN115146562A (zh) | 2022-10-04 |
Family
ID=83412509
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210824461.5A Pending CN115146562A (zh) | 2022-07-14 | 2022-07-14 | 一种基于中尺度气象要素边界速度场赋值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115146562A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117057164A (zh) * | 2023-10-11 | 2023-11-14 | 中国船舶集团风电发展有限公司 | 风资源评估方法、装置、电子设备及存储介质 |
-
2022
- 2022-07-14 CN CN202210824461.5A patent/CN115146562A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117057164A (zh) * | 2023-10-11 | 2023-11-14 | 中国船舶集团风电发展有限公司 | 风资源评估方法、装置、电子设备及存储介质 |
CN117057164B (zh) * | 2023-10-11 | 2024-01-12 | 中国船舶集团风电发展有限公司 | 风资源评估方法、装置、电子设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109685152B (zh) | 一种基于dc-spp-yolo的图像目标检测方法 | |
CN106779154B (zh) | 基于卫星云图的区域分布式光伏出力预测方法 | |
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN105576650B (zh) | 基于同化技术与同风带理论的风电功率预测方法及系统 | |
CN103136270B (zh) | 一种获得数据插值的方法及系统 | |
CN110851790B (zh) | 一种基于深度学习算法的海流动力优化预报方法 | |
CN115146562A (zh) | 一种基于中尺度气象要素边界速度场赋值方法 | |
CN109444955B (zh) | 三维地震射线追踪的双线性走时扰动插值方法 | |
WO2013121340A1 (en) | Digital elevation model | |
JP7153330B2 (ja) | 土砂災害予測装置、コンピュータプログラム、土砂災害予測方法及び地図情報 | |
CN113095552B (zh) | 一种海洋平台泄漏可燃气云团体积预测方法及系统 | |
CN107544098B (zh) | 地表粗糙度的生成方法、装置、存储介质和处理器 | |
Navarro-Pons et al. | Characterization of wind-blown sediment transport with height in a highly mobile dune (SW Spain) | |
CN112347657B (zh) | 一种风能理论储量评估方法 | |
CN109033181A (zh) | 一种复杂地形地区风场地理数值模拟方法 | |
CN117351165A (zh) | 一种基于ArcGIS的三维地层模型的构建及测算方法 | |
CN109557588B (zh) | 一种煤矿井下二维矿震波速反演降维方法 | |
CN107526910A (zh) | 一种核设施事故场外后果评价中的风场诊断方法 | |
JP2021012078A (ja) | 温度分布推定装置及び温度分布推定方法 | |
CN112945198B (zh) | 一种基于激光lidar点云的输电线路铁塔倾斜自动检测方法 | |
CN113985490B (zh) | 利用地形和地壳密度数据进行地表重力仿真的方法及装置 | |
CN115512076A (zh) | 一种网格重建方法、装置、设备和存储介质 | |
US20160138984A1 (en) | Apparatus and method for calculating wind load | |
KR102218179B1 (ko) | 초단기 예보 모델의 복사 물리 모수화를 에뮬레이션 하는 시스템 및 방법 | |
CN114740540A (zh) | 一种基于方向约束的洋中脊区磁异常图构建方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication |