[发明专利]一种利用WRF模式分析面雨量误差的方法

权利要求书

1.一种利用WRF模式分析面雨量误差的方法，其特征在于，包括以下步骤：

S1，针对待分析区域构建WRF模式，确定WRF模式模拟的时间范围；

S2，获取WRF模式驱动数据，启动运行WRF模式，运行完成后获取降水场信息并存储在文本文件中；

S3，利用特定插值方法，包括雨量站算数平均插值法、泰森多边形插值法、等值线插值法中的一种，生成WRF模式各模拟时段的面雨量P_wrf，P_wrf＝{P_{wrf_1,}P_{wrf_2,…,}P_{wrf_m}}，其中m为时段个数；

S4，获取区域内雨量站位置信息，选定雨量站算数平均插值法、泰森多边形插值法、等值线插值法中的其中一种，采用选定的插值方法获取各个模拟时段的面雨量P_inter，P_inter＝{P_{inter_1,}P_{inter_2,…,}P_{inter_m}}，其中m为时段个数；

S5，对比计算P_wrf和P_inter两者之差：分别将采用相同差值法获得的P_wrf和P_inter做差，分析每一种插值方法下，面雨量误差的特征特性，从而确定最佳的面雨量插值方法。

2.根据权利要求1所述的利用WRF模式分析面雨量误差的方法，其特征在于，步骤S1具体包括：

S11，选定模拟区域：根据研究的实际需要，选定模拟区域，需要提供区域边界的经纬度序列boundary＝{{lat₁,lon₁},{lat₂,lon₂},…,{lat_n,lon_n}},其中lat为纬度，lon为经度，n为边界上的点的总个数；

S12，针对选定区域构建WRF模式：针对选定的经纬度序列构建WRF模式，使WRF模式覆盖的区域能够包含经纬度序列boundary；其空间分辨率不应低于3km×3km；

S13，确定模拟时间范围：根据研究的实际需要，确定WRF模式模拟的时间范围。

3.根据权利要求1所述的利用WRF模式分析面雨量误差的方法，其特征在于，步骤S2中获取的WRF驱动数据为FNL驱动数据；

获取降水场数据是通过NCL工具，从wrfout文件中获取得到，并以栅格形式存储于文本文件中。

4.根据权利要求1所述的利用WRF模式分析面雨量误差的方法，其特征在于，步骤S4具体包括：

S41，获取区域内雨量站位置信息：获取boundary内的雨量站经纬度位置，一般为已建雨量站经纬度位置，若区域内没有已建雨量站，也可以设置虚拟雨量站；

S42，构建雨量站与WRF输出降水场格点的对照关系：基于空间位置最近原则，构建雨量站与WRF输出降水场格点的对照关系，即获取每一个雨量站分别与WRF模式输出降水场的哪个格点对应；

S43，根据步骤S42中雨量站与输出降水场的对照关系，选定插值方法，计算各个模拟时段的面雨量P_inter，P_inter＝{P_{inter_1},P_{inter_2},…,P_{inter_m}}，其中m为时段个数。

5.根据权利要求4所述的利用WRF模式分析面雨量误差的方法，其特征在于，步骤S43中的插值方法包括雨量站算数平均插值法、泰森多边形插值法、等值线插值法中的至少一种。