一种基于广域数据的风场设计方法与流程

本发明涉及风场研究，尤其涉及一种基于广域数据的风场设计方法。

背景技术：

1、空间大气风场是重要的环境因素。最大风切变在飞行器的姿态控制系统设计、制导精度分析和载荷分析中，起到尤为重要的作用。分析空间大气风切变的特性，研究其对飞行器的影响对精确飞行有重要的参考价值。风速是气象研究中一个很重要的参量，它与大气环流、大气能量、大气质量有密切的关系。目前已经有关于地表面风速场的分析，发现在近地面层年平均和季节平均风速均存在显著减弱的变化趋势，而对高空风速场还缺少一个全面系统的分析研究，本发明针对小尺度范围内风场缺乏必要的数值预报方法的现状，在研究中尺度数值预报模式wrf模式的基础上，进行高分辨率的精细化风场模拟研究。实时业务化风场预报的需求，建立基于wrf预报模式的业务化风场实时预报平台，实现初始场资料下载整理、初始场数据预处理、模式运行、预报产品制作的全自动化。开展空间大气风场预报，建立可靠的风场预报模型，有助于人类更深入地理解大气空间，为人类开展活动提供保障。

技术实现思路

1、本发明的目的是为了解决现有技术中存在的缺点，提出一种基于广域数据的风场设计方法。

2、为了实现上述目的，本发明采用了如下技术方案：一种基于广域数据的风场设计方法，具体包括以下步骤：

3、s1、风场数据的获取：

4、s1.1、数据的获取：通过气球将无线电探空仪升空，气球在空中以300～500米/分钟的速度上升，荷载质量为1～2千克，探测高度在0～40千米，在测风的同时测量空中的温度、气压和湿度，获取的数据传输至地面观测台；

5、s1.2、矢量平均法获取风速以及风向；

6、s2、风场特征的统计：

7、s2.1、收集与整理气象探空风场资料的历史数据；

8、s2.2、对这些数据进行检查鉴别，排成数列；

9、s2.3、分析时间数列，从中寻找风场随时间变化而变化的规律，得出一定的模式；

10、s3、风场模型的构建；

11、s4、风场载荷预报系统的构建；

12、s5、风场方案设计。

13、作为一种优选的实施方式，所述步骤s1.2具体为：在两个风速分量都为零时，其矢量表示为ws(0，0)，利用每秒探测得到的斜距、仰角、方位角，按照公式计算每秒风的南北分量u与东西分量v，得到平均u(v)分量后，依据以下公式来计算平均风速(ws)和风向(wd)：

14、v＞0，则

15、v＜0，则：

16、

17、v＝0，则wd＝270，u≥0，

18、wd＝90，u＜0；

19、定义wp为连续多日一天中08时或20时风速值的矢量平均值与算术平均值的比值，则：wp的取值为0到1。

20、作为一种优选的实施方式，所述步骤s2.1具体为：从互联网上下载近10年气象探空风场数据，分别统计每日探空资料早8点和晚20点的探测数据，得到早晚不同的风场强度，风向等的统计特征。

21、作为一种优选的实施方式，所述步骤s2.2具体为：对步骤s2.1中获取的气象数据进行数据重复的筛查，以及按照日期进行排列。

22、作为一种优选的实施方式，所述步骤s2.3具体为：定义11月至翌年3月为冬季风时期，6～9月为夏季风时期，4～5月和10月为夏、冬季风转换的过渡时期，基于时间序列法，以500m为高度间隔统计不同高度层的年平均月分布特征。

23、作为一种优选的实施方式，所述步骤s3具体为：

24、s3.1、wrf动力框架的构建：wrf-arw采用的是地形追随的静力气压垂直坐标，垂直坐标定义为：η＝(ph-pht)/μ，

25、μ＝phs-pht，其中，ph为静力平衡气压分量，phs为地表面气压，pht为模式顶气压，设为常数，pht取值5hpa，μ(x，y)表示的是模式区域内坐标为(x，y)处单位面积空气柱的质量，因此各变量的通量形式可表示为：其中，v＝(u，v，w)，w＝η，θ是位温，位势高度φ＝gz，气压和密度的关系为α＝1/ρ；采用以上形式的变量，通量形式的欧拉方程组为：

26、

27、

28、

29、

30、

31、其中，静力方程为：状态方程为：p＝p0(rdθ/p0α)γ，以上方程中，x，y和η表示各方向的微分：

32、其中，a表示任何一个变量，γ＝cp/cv＝1.4，rd为干空气气体常数，p0为参考气压，fu，fv，fw，和fθ分别表示模式微物理过程、湍流混合、球面投影和地球旋转引起的强迫项；

33、s3.2、预报模型的构建：预报变量表示为：

34、

35、

36、

37、由此可得到动量方程的扰动形式：

38、

39、

40、质量守恒方程和位势方程为：

41、

42、位温和水汽方程保持不变：

43、

44、扰动形式的静力平衡关系为：

45、湿空气的状态方程：p＝p0(rdθm/p0αd)γ，wrf中尺度模式的时间积分采用完全时间分裂格式，外循环使用较大时间步长的三阶或者四阶的runge-kutta算法，内循环为声波时间积分。

46、作为一种优选的实施方式，所述步骤s4具体为：

47、s4.1、大涡模拟流程：根据步骤s3中，结合wrf-les大涡模拟，设置在一定的区域范围内，模拟时长24–48小时情况下，风速和风向的演变特性，通过调整wrf-les大涡模式的初始场参数，对比模拟分析不同初始场情况下的风场特性，具体如下：在wrf模式和les模拟区域之间采用双向嵌套，运行双向嵌套是在不同的网格分辨率下同时运行多个域并且各个区域之间存在相互作用，粗网格提供了边界值，细网格则会将一部分计算值返回给粗网格，从而形成反馈，其中，采用单向嵌套的ndown，wps预处理系统个中geogrid.exe为地形处理程序，在水平方向上进行网格划分并对地面静态数据进行插值，生成地形文件，地形文件与初始气象数据met文件经过wps预处理系统生成的垂直方向分层met文件，met文件作为运行wrf模式real.exe的输入进行初始化插值处理，把运行real.exe的输出作为ndown.exe程序的输入，最终将ndown程序的输出作为wrf模式wrf.exe的输入进行模拟运算；

48、s4.2、风场载荷预报系统的构建：基于mpi并行运算方式，通过过编写定时脚本程序并调用shell脚本实现预报平台的自动化运行，shell是命令解释性语言，是unix/linux系统的用户界面，提供了用户与内核交互操作的接口，接受用户输入的命令并把它送入内核去执行，为公用程序，在用户登录的时候启动，采用bash shell版本，并且融入了c shell和k shell的功能。

49、作为一种优选的实施方式，步骤s4.2中，风场载荷预报系统的具体工作流程如下：

50、s4.2.1、模式参数的传递与模式的定时启动：通过设定参数件路径、起始时间和预报长度，实现对模式运行过程中的namelist自动修改；

51、s4.2.2、模式的定时启动：将crontab作为unix/linux系统下的定时任务触发器，通过预先创建crontab文件，crontab文件有六个域，每个域之间用空格或者tab键分开，之后通过cron服务使用户能够在指定的时间执行某些程序；

52、s4.2.3、初值资料的自动获取：连接美国ncep数据服务器，通过shell的wget命令每日定时下载gfs模式数据，wget是unix/linx环境下用于从互联网提取文件的工具，支持代理服务器和断点续传功能，能够自动地递归远程主机的目录，找到符合条件的文件并将其下载到本地服务器；

53、s4.2.4、预报产品制作的自动化：模式自动运行结束后，对模式预报结果进行后处理，之后生成预报的产品。

54、作为一种优选的实施方式，根据步骤s4建立的预报模型，得到风场在某段时间存在的大致的确定的方向∝0，分别计算最多风向、99％概率最大风速、平均条件风、最小条件风和最大风切变，进而给出综合矢量风，基于一定概率的风场，进行大涡模拟，补充说明高层大气风场特性，具体为：

55、s5.1、计算最大概率风速分布：

56、此时，风速w的分布条件为：是ωs标准正态分布函数；设定pω＝99％，则：根据ωs＝aω*-b/a得到最多风向时，最大风速ω*；

57、hi为参考高度，ui和vi分别为hi高度的纬向风和经向风，h为其他高度，u和v分别为h高度的纬向风和经向风，当h＞hi时：u′＝u-ui，v′＝v-vi，当h＜hi时：u′＝ui-u，v′＝vi-v，条件风切变的数学期望为：标准偏差为：u*和v*分别为最大风速下的纬向分量和经向分量；u′、v′、σu′、σv′分别为纬向风切变、经向风切变及其相对应的数学期望和标准偏差，ρuu′、ρvv′为风速与切变间的相关系数，由计算可知，σ(u′|u*)≈σ(v′|v*)，而且ρvv′很小，所以切边(实为条件切变)为二维圆分布，条件风标准偏差为：

58、条件风的统计特征量：在参考高度hi上出现最大风条件下，其他高度h上的风都是由风切边而得来的，定义为条件风，条件风的数学期望：

59、在h＜hi，e(u|u*)＝u*-e(u′|u*)，

60、e(v|v*)＝v*-e(v′|v*)，

61、在h＞hi，e(u|u*)＝u*+e(u′|u*)，

62、e(v|v*)＝v*+e(v′|v*)，

63、在h＝hi，e(u′|u*)＝0，

64、e(u|u*)＝u*，

65、e(v′|v*)＝0，

66、e(v|v*)＝v*；

67、s5.2、m*为经向风分量和纬向风分量比值，计算最小条件风，公式为：m*＝tanα*，α*＝270°-θ*，最多风向平面与对顶锥的交点坐标为：v*＝m*u*，其中，u*和v*分别为最小条件风的纬向分量和经向分量，设定p(λc)为出现概率，取值99％，则，最小条件风w*为：根据计算后的风场数据给出飞行建议。

68、与现有技术相比，本发明的优点和积极效果在于：

69、本发明利用新的风场计算方法及新的统计方法得到风场特性，并且建立了一个具有高时空分辨率的历史风场数据集，为检验风场数值模拟结果及实时预报产品提供历史统计风场数据集。

70、本发明根据不同地点是否受季风的影响，划分不同的四季，非季风区按照一年四季，春夏秋冬来划分，而即使都是季风区，但每个地点受到季风影响时间不同根据实际情况划分季节，得到更合理的风场特性研究分析。

71、基于历史气象风场数据，结合高分辨率wrf-les大涡模型模拟分析风场特性，搭建所需的风场预报系统，实现业务化风场实时预报，且实现自动化的实时风场预报。