气象要素缺失时典型气象年的生成方法及装置.pdf

《气象要素缺失时典型气象年的生成方法及装置.pdf》由会员分享，可在线阅读，更多相关《气象要素缺失时典型气象年的生成方法及装置.pdf（15页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910872192.8 (22)申请日 2019.09.16 (71)申请人 北京工业大学 地址 100124 北京市朝阳区平乐园100号 (72)发明人 张晓静盖世博谢静超刘加平 (74)专利代理机构 北京华专卓海知识产权代理 事务所(普通合伙) 11664 代理人 彭锐 (51)Int.Cl. G06F 30/20(2020.01) G06F 111/10(2020.01) (54)发明名称 一种气象要素缺失时典型气象年的生成方 法及装置 (57)摘要 本申请提供了一种。

2、气象要素缺失时典型气 象年的生成方法。 包括以下步骤： 对气象要素进 行相关性分析， 得到与所述缺失气象要素相关性 不小于阈值的其它气象要素； 使用所述其它气象 要素替代所述缺失气象要素， 通过桑迪亚国家实 验室统计方法获取表示长期累积频率分布函数 接近程度的FS值； 通过所述每个月的累积频率分 布函数和FS值， 计算得到逐年各月份的平均加权 总值WS； 根据所述WS值生成典型气象年。 以此方 式， 可以解决部分气象要素缺失， 难以生成典型 气象年的问题。 权利要求书2页 说明书10页 附图2页 CN 110781577 A 2020.02.11 CN 110781577 A 1.一种气象要素。

3、缺失时典型气象年的生成方法， 其特征在于， 包括： 对气象要素进行相关性分析， 得到与所述缺失气象要素相关性不小于阈值的其它气象 要素； 使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国家实验室统计方法获取 表示长期累积频率分布函数接近程度的FS值； 通过所述表示长期累积频率分布函数接近程度的FS值， 计算得到逐年各月份的平均加 权总值WS； 根据所述WS值生成典型气象年。 2.根据权利要求1所述的方法， 其特征在于， 所述对气象要素进行相关性分析， 得到与 所述缺失气象要素相关性不小于阈值的其它气象要素包括： 通过相关性分析公式计算所述缺失气象要素与其它气象要素之间的确定系数R2； 。

4、得到R2大于等于阈值的其它气象要素。 3.根据权利要求2所述的方法， 其特征在于， 所述相关性分析公式为： 其中， R2为确定系数； Yi为某气象要素第i个观测值； Xi为另一气象要素第i个观测值； 为某气象要素观测平均值。 4.根据权利要求3所述的方法， 其特征在于， 根据下式计算所述FS值： 其中， n为参数值选取个数； FSx(y,m)为第x个气象要素在xi范围的FS(y,m)统计值， y为研究对象年， m为研究对象 月； CDFy,m(xi)为第x个气象要素在xi范围内的累积频率分布函数； CDFm(xi)为对于月份m， 第x个气象要素长期统计值域在xi范围的累积频率分布函数。 5.根。

5、据权利要求4所述的方法， 其特征在于， 根据下式计算所述WS： 其中， WSx(y,m)为y年m月的平均加权和； WFx为第x个气象要素的加权因子， 替代的气象要素的加权因子按照被替代的缺失的气 象要素的加权因子选取； m为气象要素选取个数。 6.根据权利要求5所述的方法， 其特征在于， 所述根据所述WS值生成典型气象年包括： 选取最小WS值的月份作为典型月； 权利要求书 1/2 页 2 CN 110781577 A 2 将12个所述典型月连接得到第一典型气象年。 7.根据权利要求6所述的方法， 其特征在于， 还包括： 按照气象要素所对应的气象参数将所述气象要素进行分组； 使用同组的一个或多个。

6、其它气象要素替代所述缺失的气象要素， 生成第二典型气象 年； 根据第一典型气象年和第二典型气象年生成各要素的典型月CDF曲线； 分别将所述典型月CDF曲线与长期实测数据CDF曲线进行对比， 选取吻合度更好的典型 气象年作为最终的典型气象年。 8.一种气象要素缺失时典型气象年的生成装置， 其特征在于， 包括： 分析模块， 用于对气象要素进行相关性分析， 得到与所述缺失气象要素相关性不小于 阈值的其它气象要素； 替换模块， 用于使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国家实验 室统计方法获取表示长期累积频率分布函数接近程度的FS值； 计算模块， 用于通过所述表示长期累积频率分布函数接。

7、近程度的FS值， 计算得到逐年 各月份的平均加权总值WS； 第一生成模块， 用于根据所述WS值生成典型气象年。 9.根据权利要求8所述的装置， 其特征在于， 还包括： 第二生成模块， 用于按照气象要素所对应的气象参数将所述气象要素进行分组； 使用同组的一个或多个其它气象要素替代所述缺失的气象要素， 生成第二典型气象 年； 对比模块， 用于将根据第一典型气象年和第二典型气象年生成各要素的典型月CDF曲 线， 分别与长期实测数据CDF曲线进行对比； 选取模块， 用于选取与长期实测数据CDF曲线吻合度更好的典型气象年作为最终的典 型气象年。 权利要求书 2/2 页 3 CN 110781577 A 。

8、3 一种气象要素缺失时典型气象年的生成方法及装置 技术领域 0001 本发明涉及建筑能耗领域， 特别地涉及一种气象要素缺失时典型气象年的生成方 法及装置。 背景技术 0002 能源对于中国的经济和社会发展至关重要。 能源保护已经成为人们关注的热门话 题。 建筑物的空调能耗是主要的能源消耗者之一， 其目的是移除室内负荷。 而建筑物的负荷 除了建筑物本身的设计以外， 还与外界气候条件密不可分。 在暖通空调行业， 进行科学研究 或工程设计， 都需依靠气象数据对建筑物的冷热能耗及建筑热环境进行计算分析。 气象参 数对指导地区的建筑设计标准及能耗模拟评估起着基础、 直接的作用， 所以在研究建筑物 能耗模。

9、拟前尚需掌握大量的各种气象要素， 才能真实地反映建筑物的热负荷状况。 0003 所有能耗模拟程序都需要输入气象数据来驱动模拟程序中的热模型。 由于每年的 气象条件可能会发生显著变化， 因此需要得出一个可以很好地代表长期平均气象状况的定 制气象数据集。 典型气象年TMY(Typical Meteorological Year)， 就是一个地区具有典型 气候特征的一年， 它是由12个均具有气候代表性的典型月(TMM)组成一个 “假想” 气象年， 是 一年的代表性数据库。 0004 典型气象年的生成方法有Sandia国家实验室法、 Danish法、 Festa-Ratto法、 CSWD 法、 CTY。

10、W法、 主成分法等。 上述方法通过对比所选月份各气象要素与长期气象数据的接近程 度来挑选典型气象月， 从而构成典型气象年。 不同的方法利用的气象要素类型略有差异， 一 般利用温度、 湿度、 风速、 太阳辐射等多种气象要素。 0005 然而， 由于各地区气象观测站点建站时间不一， 且部分地区气候要素资料数据记 录存在缺失问题， 尤其是气象辐射资料数据还需进一步完善， 导致缺少部分气象要素， 难以 进行气象年的生成， 从而无法真实地模拟和反映建筑物的热负荷。 发明内容 0006 本发明的目的在于提供一种气象要素缺失时典型气象年的生成方法及装置， 用以 解决部分气象要素缺失， 难以生成典型气象年的问。

11、题。 0007 为此， 根据本发明的第一方面提供了一种气象要素缺失时典型气象年的生成方 法， 包括以下步骤： 0008 对气象要素进行相关性分析， 得到与所述缺失气象要素相关性不小于阈值的其它 气象要素； 0009 使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国家实验室统计方法 获取表示长期累积频率分布函数接近程度的FS值； 0010 通过所述表示长期累积频率分布函数接近程度的FS值， 计算得到逐年各月份的平 均加权总值WS； 0011 根据所述WS值生成典型气象年。 说明书 1/10 页 4 CN 110781577 A 4 0012 进一步地， 所述对气象要素进行相关性分析， 得。

12、到与所述缺失气象要素相关性不 小于阈值的其它气象要素包括： 0013 通过相关性分析公式计算所述缺失气象要素与其它气象要素之间的确定系数R2； 0014 得到R2不小于阈值的其它气象要素。 0015 进一步地， 所述相关性分析公式为： 0016 0017 其中， R2为确定系数； 0018 Yi为某气象要素第i个观测值； 0019 Xi为另一气象要素第i个观测值； 0020为某气象要素观测平均值。 0021 进一步地， 根据下式计算所述FS值： 0022 0023 n为参数值选取个数； 0024 FSx(y,m)为第x个气象要素在xi范围的FS(y,m)统计值， y为研究对象年， m为研究 对。

13、象月； 0025 CDFy,m(xi)为第x个气象要素在xi范围内的累积频率分布函数； 0026 CDFm(xi)为对于月份m， 第x个气象要素长期统计值域在xi范围的累积频率分布函 数。 0027 进一步地， 根据下式计算所述WS： 0028 0029 0030 其中， WSx(y,m)为y年m月的平均加权和； 0031 WFx为第x个气象要素的加权因子， 替代的气象要素的加权因子按照被替代的缺失 的气象要素的加权因子选取； 0032 m为气象要素选取个数。 0033 进一步地， 所述根据所述WS值生成典型气象年包括： 0034 选取最小WS值的月份作为典型月； 0035 将12个典型月连接。

14、得到第一典型气象年。 0036 进一步地， 还包括： 0037 按照气象要素所对应的气象参数将所述气象要素进行分组； 0038 使用同组的一个或多个其它气象要素替代所述缺失的气象要素， 生成第二典型气 象年； 0039 根据第一典型气象年和第二典型气象年生成各要素的典型月CDF曲线； 0040 分别将所述典型月CDF曲线与长期实测数据CDF曲线进行对比， 选取吻合度更好的 说明书 2/10 页 5 CN 110781577 A 5 典型气象年作为最终的典型气象年。 0041 根据本发明的第二方面， 提供了一种气象要素缺失时典型气象年的生成装置， 包 括： 0042 分析模块， 用于对气象要素进。

15、行相关性分析， 得到与所述缺失气象要素相关性不 小于阈值的其它气象要素； 0043 替换模块， 用于使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国家 实验室统计方法获取表示长期累积频率分布函数接近程度的FS值； 0044 计算模块， 用于通过所述表示长期累积频率分布函数接近程度的FS值， 计算得到 逐年各月份的平均加权总值WS； 0045 第一生成模块， 用于根据所述WS值生成典型气象年。 0046 进一步地， 还包括： 0047 第二生成模块， 用于按照气象要素所对应的气象参数将所述气象要素进行分组； 0048 使用同组的一个或多个其它气象要素替代所述缺失的气象要素， 生成第二典型。

16、气 象年。 0049 对比模块， 用于将根据第一典型气象年和第二典型气象年生成各要素的典型月 CDF曲线， 分别与长期实测数据CDF曲线进行对比； 0050 选取模块， 用于选取与长期实测数据CDF曲线吻合度更好的典型气象年作为最终 的典型气象年。 0051 本发明通过对缺失气象要素进行分析， 得到补充所述缺失气象要素的其它气象要 素， 然后通过桑迪亚国家实验室统计方法， 计算得到典型气象年， 从而在缺失部分气象要素 时， 也能够生成典型气象年。 同时还能较为准确、 真实地模拟和反映建筑物的热负荷。 0052 本发明的附加方面和优点将在下面的描述部分中给出， 部分将从下面的描述中变 得明显， 。

17、或通过本发明的实践了解到。 附图说明 0053 本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得 明显和容易理解， 其中： 0054 图1是根据本发明的气象要素缺失时典型气象年的生成方法的流程示意图。 0055 图2是根据本发明的气象要素缺失时典型气象年的生成装置的框架示意图。 0056 图3是根据本发明一实施例的西沙地区在2003-2014年间实测日照时数与日总辐 射数据的散点图。 0057 图4是根据本发明一实施例的CDF曲线对比结果分类示意图。 具体实施方式 0058 为了能够更清楚地理解本发明的上述目的、 特征和优点， 下面结合附图和具体实 施方式对本发明进行进一。

18、步的详细描述。 需要说明的是， 在不冲突的情况下， 本申请的实施 例及实施例中的特征可以相互组合。 0059 在下面的描述中阐述了很多具体细节以便于充分理解本发明， 但是， 本发明还可 以采用其他不同于在此描述的其他方式来实施， 因此， 本发明的保护范围并不受下面公开 说明书 3/10 页 6 CN 110781577 A 6 的具体实施例的限制。 0060 图1示出了根据本发明的一种气象要素缺失时典型气象年的生成方法， 包括以下 步骤： 0061 S11、 对气象要素进行相关性分析， 得到与所述缺失气象要素相关性不小于阈值的 其它气象要素； 0062 其中， 通过观测站得到的气象要素信息， 。

19、由于观测站建站时间、 数量、 分布等问题， 会出现缺失典型气象年计算所需要的部分气象要素信息的情况， 此时则需要用其它气象要 素对所述缺失气象要素进行补充。 0063 根据相关性分析公式得到与所述缺失气象要素相关性较高的其它气象要素。 具体 地， 将确定系数R2与预设阈值相比较， 若所述确定系数R2大于等于预设阈值(例如为0.8)， 即 可判定所述其它气象要素与所述缺失气象要素相关性较高,由此从观测得到的气象要素信 息中选出可替代缺失气象要素的其它气象要素。 对于一个缺失气象要素， 可能存在多个可 替代的其它气象要素， 此时可以从中任意选择一个， 也可以选择相关系数最高的一个。 具体 地， 相。

20、关性分析公式如下所示： 0064 0065 R2为确定系数； 0066 Yi为某气象要素第i个观测值； 0067 Xi为另一气象要素第i个观测值； 0068为某气象要素观测平均值。 0069 S12、 使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国家实验室统计 方法获取表示长期累积频率分布函数接近程度的FS值； 0070 其中， 利用桑迪亚(Sandia)国家实验室统计方法(又称Filkenstein-Schafer(FS) 统计方法)， 计算表征干球温度、 露点温度、 风速和太阳辐射等各气象参数的各月份累积频 率分布函数(Cumulative Distribution Functi。

21、on， 简称CDF)与长期累积频率分布函数的 接近程度的FS值， 计算中缺失的气象要素用步骤S11中得到相关性较高的其它气象要素替 代， 具体计算公式如下： 0071 0072 n为参数值选取个数； 0073 FSx(y,m)为第x个气象要素在xi范围的FS(y,m)统计值， y为研究对象年， m为研究 对象月； 0074 CDFy,m(xi)为第x个气象要素在xi范围内的累积频率分布函数； 0075 CDFm(xi)为对于月份m， 第x个气象要素长期统计值域在xi范围的累积频率分布函 数。 0076 S13、 通过所述表示长期累积频率分布函数接近程度的FS值， 计算得到逐年各月份 的平均加权。

22、总值WS； 0077 其中， 根据下述公式计算得到逐年各月份的平均加权总值WS： 说明书 4/10 页 7 CN 110781577 A 7 0078 0079 0080 其中， WSx(y,m)为y年m月的平均加权和； 0081 WFx为第x个气象要素的加权因子， 替代的气象要素的加权因子按照被替代的缺失 的气象要素的加权因子选取； 0082 m为气象要素选取个数。 0083 S14、 根据所述WS值生成典型气象年； 0084 其中， 在某些情况下会导致没有候选月份或筛选掉一些有用的月份， 因而选择具 有最小WS值的月份作为典型月， 将12个典型月连接得到第一典型气象年。 0085 进一步地。

23、， 所述方法还包括： 0086 按照气象要素所对应的气象参数将所述气象要素进行分组； 0087 使用同组的一个或多个其它气象要素替代所述缺失的气象要素， 生成第二典型气 象年； 0088 根据第一典型气象年和第二典型气象年生成各要素的典型月CDF曲线； 0089 分别将所述典型月CDF曲线与长期实测数据CDF曲线进行对比， 选取吻合度更好的 典型气象年作为最终的典型气象年。 0090 其中， 当缺失的气象要素为某气象参数的子集时， 可将该缺失气象要素的权重累 加到同组气象要素中， 进而生成第二典型气象年。 例如， 温度为气象参数， 干球温度日最高 值与最低值和干球温度日平均值为所述温度气象参数。

24、的子集， 当干球温度日最高值与最低 值缺失时， 可将其权重累加到干球温度日平均值的权重中。 0091 进一步地， 将第一典型气象年和第二典型气象年中各要素的各典型月CDF曲线分 别与长期实测数据CDF曲线比对， 选取吻合度更好地典型气象年作为最终的典型气象年。 0092 本发明实施例还提供了一种气象要素缺失时典型气象年的生成装置， 如图2所示， 包括： 0093 分析模块21， 用于对气象要素进行相关性分析， 得到与所述缺失气象要素相关性 不小于阈值的其它气象要素； 0094 其中， 通过观测站得到的气象要素信息， 由于观测站建站时间、 数量、 分布等问题， 会出现缺失典型气象年计算所需要的部。

25、分气象要素信息的情况， 此时则需要用其它气象要 素对所述缺失气象要素进行补充。 0095 根据相关性分析公式得到与所述缺失气象要素相关性较高的其它气象要素。 具体 地， 将确定系数R2与预设阈值相比较， 若所述确定系数R2大于等于预设阈值(例如为0.8)， 即 可判定所述其它气象要素与所述缺失气象要素相关性较高,由此从观测得到的气象要素信 息中选出可替代缺失气象要素的其它气象要素。 对于一个缺失气象要素， 可能存在多个可 替代的其它气象要素， 此时可以从中任意选择一个， 也可以选择相关系数最高的一个。 具体 地， 相关性分析公式如下所示： 0096 说明书 5/10 页 8 CN 110781。

26、577 A 8 0097 R2为确定系数； 0098 Yi为某气象要素第i个观测值； 0099 Xi为另一气象要素第i个观测值； 0100为某气象要素观测平均值。 0101 替换模块22， 用于使用所述其它气象要素替代所述缺失气象要素， 通过桑迪亚国 家实验室统计方法获取表示长期累积频率分布函数接近程度的FS值； 0102 其中， 利用桑迪亚(Sandia)国家实验室统计方法(又称Filkenstein-Schafer(FS) 统计方法)， 计算表征干球温度、 露点温度、 风速和太阳辐射等各气象参数的各月份累积频 率分布函数(Cumulative Distribution Function， 。

27、简称CDF)与长期累积频率分布函数的 接近程度的FS值， 计算中缺失的气象要素用分析模块21中得到相关性较高的其它气象要素 替代， 具体计算公式如下： 0103 0104 n为参数值选取个数； 0105 FSx(y,m)为第x个气象要素在xi范围的FS(y,m)统计值， y为研究对象年， m为研究 对象月； 0106 CDFy,m(xi)为第x个气象要素在xi范围内的累积频率分布函数； 0107 CDFm(xi)为对于月份m， 第x个气象要素长期统计值域在xi范围的累积频率分布函 数。 0108 计算模块23， 用于通过所述表示长期累积频率分布函数接近程度的FS值， 计算得 到逐年各月份的平均。

28、加权总值WS； 0109 其中， 根据下述公式计算得到逐年各月份的平均加权总值WS： 0110 0111 0112 其中， WSx(y,m)为y年m月的平均加权和； 0113 WFx为第x个气象要素的加权因子， 替代的气象要素的加权因子按照被替代的缺失 的气象要素的加权因子选取； 0114 m为气象要素选取个数。 0115 第一生成模块24， 用于根据所述WS值生成典型气象年。 0116 其中， 在某些情况下会导致没有候选月份或筛选掉一些有用的月份， 因而选择具 有最小WS值的月份作为典型月， 将12个典型月连接得到第一典型气象年。 0117 进一步地， 所述装置还包括： 0118 第二生成模。

29、块， 用于按照气象要素所对应的气象参数将所述气象要素进行分组； 0119 使用同组的一个或多个其它气象要素替代所述缺失的气象要素， 生成第二典型气 象年。 0120 其中， 当缺失的气象要素为某气象参数的子集时， 可将该缺失气象要素的权重累 加到同组气象要素中， 进而生成第二典型气象年。 例如， 温度为气象参数， 干球温度日最高 说明书 6/10 页 9 CN 110781577 A 9 值与最低值和干球温度日平均值为所述温度气象参数的子集， 当干球温度日最高值与最低 值缺失时， 可将其权重累加到干球温度日平均值的权重中。 0121 对比模块， 用于将根据第一典型气象年和第二典型气象年生成各要。

30、素的典型月 CDF曲线， 分别与长期实测数据CDF曲线进行对比； 0122 将第一典型气象年和第二典型气象年中各要素的各典型月CDF曲线分别与长期实 测数据CDF曲线比对。 0123 选取模块， 用于选取与长期实测数据CDF曲线吻合度更好的典型气象年作为最终 的典型气象年。 0124 通过对比模块的对比结果， 选取与所述长期实测数据CDF曲线吻合度更好的典型 气象年作为最终的典型气象年。 0125 下面给出根据本发明一个具体实施例的生成典型气象年流程： 0126 气象数据来源于中国气象局提供的西沙台站(112.33 E， 16.83 N)实测气象数据 集。 0127 需要用到公式如下所示： 0。

31、128(相关性分析公式) 0129 0130 0131 0132 原始气象数据的类型、 要素、 精度及记录长度列于表1。 其中， 四次定时观测时间为 北京时间02、 08、 14和20时， 各气象要素的日平均值统一用四次定时的平均值代表。 0133 本实施例实测数据取自西沙太阳辐射数据记录质量较好的年段2000-2014年。 0134 表1西沙历史观测气象数据质量概况 说明书 7/10 页 10 CN 110781577 A 10 0135 0136 Sandia国家实验室法可保证挑选的典型气象年气象数据与长期气象数据具有结 构上的一致性， 使之具有长期气象数据的代表性。 基于Sandia国家。

32、实验室法挑选的典型气 象年与长期气象资料在结构上具有较好的接近度并且去除了极端的气候状况， 使得典型气 象年具有较好的长期气候代表性。 0137 本实施例提出的用于挑选典型年的两种气象要素组合及各权重因子具体介绍见 表2。 0138 表2挑选典型年所需要的气象要素及其权重因子 0139 0140 通过相关性分析公式对实测日照时数与水平面日总辐射数据进行相关性分析， 结 果如图3所示。 可见西沙2000-2014年实测日照时数与日总辐射之间呈高度正相关(r 0.912， p0.001； 即R20.832， 所述实测日照时数与水平面日总辐射相关性高)， 即当日照 时数越长时， 水平面日总辐射量也越。

33、高。 0141 进一步地， 研究表明， 当气象要素部分参数缺失时， 也可将同属一类的气象要素权 重进行合并， 所生成的典型气象年依然具有代表性。 说明书 8/10 页 11 CN 110781577 A 11 0142 由表1和表2可知， 实测数据中缺少Sandia法所需的直射辐射数据， 本实施例尝试 用日照时数或水平面日总辐射替代Sandia法中的水平面总辐射和直接辐射要素， 并将两种 缺失的气象要素的权重因子进行合并后赋值给日照时数或水平面日总辐射， 即日照时数与 水平面日总辐射的权重因子均设置为10/20。 0143 通过公式1和公式2计算得到表征各气象要素的观测数据在各月份的累积频率分。

34、 布函数(Cumulative Distribution Function， 简称CDF)与长期累积频率分布函数的接近 程度的FS值， 再结合表2所列的各气象要素的权重因子， 通过公式3计算FS值的加权总值WS； 0144 进一步地， 根据WS值为每个月份挑选出5个候选月， 通过评估候选月气温和太阳辐 射的连续性， 然后按照得分最高的月份挑选与长期观测数据最为 “接近” 的典型气象月TMM (Typical Meteorological Month)， 将12个典型气象月连接得到典型气象年。 进一步地， 根 据所述表2组合中数据计算得到的典型气象年为第一典型气象年， 根据所述表2组合中 数据计。

35、算得到的典型气象年为第二典型气象年。 0145 基于西沙2000-2014年实测数据， 利用上述两组气象要素分别挑选西沙的典型气 象年。 典型气象年的挑选结果于表3， 可见有7个典型月的挑选结果是一致的， 分别为3、 6、 7、 8、 10、 11和12月。 0146 表3两组气象要素挑选典型气象年结果对比 0147 0148 为了比较典型气象年的准确性与代表性， 计算了两组典型年的RMSE与MBE值及相 应的百分比， 见表4。 可以看出两组TMY的干球温度、 露点温度、 水平面日总辐射、 日照时数、 相对湿度都具有较小的RMSE值， 这表明与长期观测值有较好的接近度； 而两组TMY的风速 R。

36、MSE值都偏大， 这可能是由于风速的随机性造成的。 两组TMY的日照时数的MBE值为负， 说明 典型年的日照时数平均值略小于长期观测平均值， 其他气象要素MBE值为正， 说明其略大于 长期观测平均值。 0149 综上来看， 两组TMY除了风速以外均具有较小的RMSE值和MBE值， 说明两组典型年 数据均与长期观测值有着较好的接近程度， 具有长期观测数据的代表性。 0150 表4两组气象要素挑选典型气象年与对应WS值 0151 0152 如图4所示， 为直观对比两种典型气象年的挑选结果， 统计了各气象要素典型月数 据的CFD曲线以及相应FS值。 同时， 以西沙台站长期观测气象数据(2000-20。

37、14)的CDF作为参 照曲线评价典型气象年挑选的准确性。 以风速的某4个典型月挑选结果为例， 对比结果大致 分为四类： 以日照时数挑选的典型月CDF曲线与长期实测数据的CDF曲线吻合度好， 且明显 优于以水平面日总辐射挑选的典型月结果(FSFS ， 图示)。 0153 按照上述分类标准， 各气象要素12个典型月的CDF曲线对比结果， 见表5。 可见， 以 日照时数等要素挑选典型气象年时， 得到的各要素的各典型月CDF曲线与长期实测数据CDF 曲线吻合度均很好， 干球温度的优选月份数量为12个， 风速、 相对湿度的优选月份数量均为 9个， 水平面日总辐射的优选月份有8个。 可见， 以日照时数替代。

38、水平面日总辐射挑选的典型 气象年， 准确度较高， 该方法可行。 0154 表5各气象要素的两组典型年CDF曲线与实测数据CDF曲线的对比结果分类统计 0155 0156 本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可 以通过程序来指令相关的硬件来完成， 该程序可以存储于一计算机可读存储介质中， 存储 介质包括只读存储器(Read-Only Memory， ROM)、 随机存储器(Random Access Memory， RAM)、 可编程只读存储器(Programmable Read-only Memory， PROM)、 可擦除可编程只读存 储器(Erasable P。

39、rogrammable Read Only Memory， EPROM)、 一次可编程只读存储器(One- time Programmable Read-Only Memory， OTPROM)、 电子抹除式可复写只读存储器 (Electrically-Erasable Programmable Read-Only Memory， EEPROM)、 只读光盘(Compact Disc Read-Only Memory， CD-ROM)或其他光盘存储器、 磁盘存储器、 磁带存储器、 或者能够 用于携带或存储数据的计算机可读的任何其他介质。 0157 以上仅为本发明的优选实施例而已， 并不用于限制本发明， 对于本领域的技术人 员来说， 本发明可以有各种更改和变化。 凡在本发明的精神和原则之内， 所作的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围之内。 说明书 10/10 页 13 CN 110781577 A 13 图1 图2 说明书附图 1/2 页 14 CN 110781577 A 14 图3 图4 说明书附图 2/2 页 15 CN 110781577 A 15 。