一种基于统计的逐时太阳辐射数据计算方法2.docVIP

一种基于统计的逐时太阳辐射数据计算方法2行政论文范文大全 一种基于统计的逐时太阳辐射数据计算方法 尽管在各个云量水平下，总辐射/天文辐射和直接辐射/总辐射的比例系数分布都满足正态分布，但是由于数据分布范围较广，其中还存在一些极端值，因此笔者认为，如果直接采用它们的算术平均值或中位值可能会带来较大的误差。因此，笔者采用了hampels redescending m-estimator最大似然均值估计算法，这种算法在数据分布符合正态分布规律，但是分布范围较广，以及有少数极端值存在的情况下，可以得到比算术平均更加接近真值的平均值。表1中列出了根据典型气象年逐时数据，在各个云量水平下，采用hampels redescending m-estimator最大似然均值估计算法得到的总辐射/天文辐射比例系数和直接辐射/总辐射比例系数。 表1 不同云量下总辐射/天文辐射及直接辐射/总辐射比例系数的最大似然均值 云量 总辐射/天文辐射 直接辐射/总辐射 0 .201996 .629582 1 .225952 .707603 2 .279112 .794809 3 .286788 .821102 4 .273717 .821468 5 .276085 .796764 6 .238903 .723629 7 .199643 .612589 8 .146616 .399326 9 .089103 .134332 10 .028377 .006261 从表中数据可以看出，(1) 最多只有不到30%的太阳天文辐射能量到达地面；(2) 当云量为2、3和4时，到达地面的太阳辐射能量最多；(3) 当云量为3和4时，在到达地面的太阳辐射能量中，直接辐射所占的比例最大，超过了80%。 根据表1中的数据，可以得到两个拟合公式： (1) (2) 式中：k1 --- 总辐射/天文辐射 k2 --- 直接辐射/总辐射 cc --- 云量 公式(1)的r2值为0.9899，公式(2)的r2值为0.9873，拟合曲线见图5。 图5 总辐射/天文辐射和直接辐射/总辐射拟合曲线 在得到上述拟合公式以后，就可以计算实际气象年的逐时太阳辐射数据了。具体步骤如下： 1) 利用天文学公式计算出当地全年逐时太阳天文辐射值，并剔除那些辐射值为0（太阳位于地平线以下）的时段。 2) 根据实际气象年的逐时天气状况，对云量值进行量化，并对变化剧烈的时段进行人工平滑处理，得到逐时云量值。 3) 由1) 和2) 的结果，应用公式(1) 或直接采用表1中的数据，计算得到当地实际气象年逐时水平面太阳总辐射量。 4) 由3) 的结果，应用公式(2) 或直接采用表1中的数据，计算得到当地实际气象年逐时法线方向太阳直射辐射量。 由于计算全年数据时的数据量很大，因此在计算2004年上海逐时太阳辐射数据时，所有的计算均由一笔者自编的计算机程序完成。由此得到的太阳辐射数据已经作为2004年上海逐时气象参数的一部分，应用于上海某高层建筑计算机能耗模拟与分析项目。 3 结论 1. 根据影响因子分析的结果，在忽略大气污染影响的前提下，水平面太阳辐射总量与大气层外水平面太阳辐射总量之间的关系只与云量相关。同样，法线方向太阳直射辐射量与水平面太阳辐射总量的关系也只与云量相关。 2. 尽管在各个不同的云量水平下，水平面太阳辐射总量与大气层外水平面太阳辐射总量的比例系数，以及法线方向太阳直射辐射量与水平面太阳辐射总量的比例系数，都分布在一个相当宽的范围内，但是这些数据都很好地服从正态分布。 3. 上述两个与不同云量水平相对应的比例系数，可以采用hampels redescending m-estimator最大似然均值估计算法计算得到，这种算法能得到比算术平均更加接近真值的平均值。这两个比例系数的最大似然均值与云量的对应关系，可以很好地分别用一个二次三项式来拟合。 4. 结合采用天文学公式计算得到的当地大气层外水平面逐时太阳辐射总量数据，以及从当地实际气象年的逐时天气状况中得到的云量数据，应用上述两个比例系数，就可以得到实际气象年的逐时太阳辐射数据。 5. 由于本计算方法的基础是当地的典型气象年逐时太阳辐射数据和当地实际气象年的逐时云量数据，因此本方法可以适用于任何地点和任何年份的逐时太阳辐射数据计算。 参考文献 1. 郎四维. 建筑能耗分析逐时气象资料的开发研究. 暖通空调. 2002年第4期. 2. 张素宁 田胜元. 太阳辐射逐时模型的建立. 太阳能学报. 1997年第3期. 3. 王炳忠. 太阳辐射计算讲座 第一讲太阳能中天文参数的计算. 太阳能. 1999年第2期. 4. 北京市气象局气候资料室编著. 北京城市气候. 气象