日散射辐射模型经验系数确定方法

曹其梦,杨 柳,于 瑛

(1.西部绿色建筑国家重点实验室(西安建筑科技大学)，西安 710055； 2.西安建筑科技大学 建筑学院,西安 710055； 3.西安建筑科技大学 机电工程学院,西安 710055)

到达地面的太阳辐射分为直射辐射和散射辐射，散射辐射主要受大气成分、云量的影响；准确可靠的散射辐射数据是建筑能耗模拟、负荷计算、遮阳设计、太阳能利用的基础，计算斜面的太阳辐射数据也需要水平面的直射、散射辐射数据做支撑. 使用高精度的测量设备实时采集是获得准确、足量数据的最有效的方法，但限于资金和维护等问题，中国辐射观测台站较少. 1993年，辐射观测站调整为98个，其中仅有17个台站记录了水平面散射辐射数据，并且呈现了东多稀少、平原多山地少的分布特征[1]；在已有的观测台站中，连续记录超过10年的数据大多集中在20世纪80～90年代初，近20年的直、散辐射数据严重匮乏. 因此，亟需选择适当的模型计算中国不同地区散射辐射数据.

散射辐射模型的研究起源于20世纪60年代，发展至今，出现了上百篇相关计算方法的研究，总结起来，主要为建立散射比(水平面日散射辐射与水平面日总辐射之比)或散射率(定义为水平面散射辐射与水平面天文辐射之比)与晴空指数、日照时数比、云量等气象因子之间的相关关系，包括单个因子的线性、多项式、指数、对数、幂函数形式以及考虑多因子综合作用的多项式形式. 根据计算尺度的不同，可分为月均值与日值，其中计算月均值的模型较多，日值模型包含在月均值模型中，详见表1.

表1 日散射辐射模型类别及表达式

注：H为总辐射,MJ·m-2;H0为天文辐射,MJ·m-2;Hd为散射辐射,MJ·m-2;S为日照时数，h；S0为可照时数，h；Tavg为平均温度，℃；RH为相对湿度，%；δ为太阳赤纬角，°；φ为纬度，°；晴空指数(kt)为总辐射与天文辐射的比值；日照时数比(S/S0)为日照时数与可照时数的比值；*为用来计算日值.

这些研究多以实测数据进行模型的建立或对模型进行修正，采用误差分析方法评价模型的优劣，进而得到适合这一地区的最优模型，包括模型形式及经验系数. 由于地区气候差异性，没有得到统一的模型；例如，Boukelia等[17]基于阿尔及利亚地区，发现以晴空指数为输入参数的二次或三次多项式形式最准确；Jamil 等[18]认为以晴空指数、日照时数比为输入参数的一次函数适合印度地区；Sabzpooshani 等[19]发现基于日照时数比的三次多项式最适合伊朗伊斯法罕地区月均散射辐射的计算. 李峥嵘等[20]研究得到了适合上海地区的日散射辐射模型；马年骏等[21]研究表明，考虑晴空指数、日照百分率的模型估算精度优于单一晴空指数的分段函数模型；冯巍等[22]以北京地区数据为基础，验证了Erbs、张晴原等10种模型的计算精度，得出晴空指数是影响散射辐射估算精度的最主要因素，考虑晴空指数为输入参数的分段函数最准确. 相关研究[17, 23]证实，这类模型的经验系数与地理位置有很大的关系，具有显著的地区差异性，对于缺少观测台站的西部及偏远的农村地区，在应用此类模型计算日散射辐射时，其经验系数取多少合适，尚有待研究. 文献[23-25]在计算总辐射时，在温度、日照时数、相对湿度等气象参数的基础上，加入了纬度、海拔高度等地理参数的修正，结果表明，修正后的模型准确度有所降低，但可以满足计算精度的要求，实现了总辐射数据的“从无到有”；文献[17]根据阿尔及利亚6个城镇的散射辐射、气象数据，拟合了散射比与日照时数比、晴空指数的关系，但这是一种平均化状态，弱化了纬度、海拔高度对太阳辐射的影响. 因此，本文根据观测台站的辐射数据建立模型经验系数与地理参数关系，以计算未观测站点的散射辐射数据. 本文以实测数据进行不同模型在不同地区的适应性分析，以准确度相对较高、形式简单的模型为基础，建立考虑纬度、海拔高度等地理参数修正的模型，对比修正前后模型的准确度，验证此类方法的可行性.

1 数据与方法

1.1 数据来源

文中所用到的数据均来自中国国家气象信息中心日值气象、辐射数据集，建模所涉及到的数据包括水平面日总辐射、水平面日散射辐射及日照时数等气象要素. 选择北京等17个台站2000—2017年观测数据，其中，每个台站前10年数据用来建模，后4年数据用来验证，各个台站经纬度等基本信息及建模、验证数据时间段见表2.

表2 17个台站的基本信息

1.2 数据质量控制

文中以北京(54511)为例，按照上述方法分别计算得到上、下包络线，将包络线与散点图叠加，如图1所示，上、下包络线之外的数据则为异常数据. 将17个台站依次完成质量控制，将异常数据删去，其余通过质量控制的数据将用于后续建模与模型验证.

图1 散射辐射数据质量控制图

1.3 模型选择

文献[27-28]证实散射率模型的准确度低于散射比模型，加入日平均温度和相对湿度的多因子的散射比模型并不能显著提高模型准确度. 因此本文不再考虑散射率模型，而是选择散射比模型，模型的输入参数为单一的晴空指数、日照时数比以及晴空指数、日照时数比共同作用的模型，模型形式为多项式函数、对数形式、指数形式，共12个模型，具体见表3.

表3 日散射辐射推测模型汇总表

1.4 模型评价方法

为了比较模型的优劣，引入了3个评价指标，分别为平均绝对误差百分比(MAE)、均方根误差百分比(RMSE)和决定系数(R2)，计算如下：

MAE、RMSE越小越好.R2数值范围介于-∞～1之间，R2越接近于1越好；R2接近于0表示计算值接近实测平均值，总体结果可信，但估算误差较大；R2<0，表示计算值小于实测值，结果不可信.

2 结果与讨论

2.1 单站点模型

以上述17个台站质量控制后的数据进行模型适应性研究，选择相应的数据建模，拟合17个模型的经验系数，如图2所示；依据相应时间段的数据进行模型验证，依次计算模型的R2、MAE和RMSE值，按晴空指数由低到高排序，绘制不同模型的MAE、RMSE及R2，结果如图3～5所示.

以系数a、b、c为例分析模型经验系数的差异，如图2所示，模型的经验系数具有台站差异性，尤其以包含晴空指数为变量模型M1、M2、M3、M4、M12最为明显，这与文献[23]的研究结果是一致的.

图2 M1～M12的经验系数

分析图3(a)～图3(c)可知：3类模型的R2，以晴空指数、日照时数比为输入参数的模型最大、波动小，其次是以晴空指数为输入参数的模型，最后是以日照时数比为输入参数的模型，MAE、RMSE的数值与R2有相反的变化趋势，表明以晴空指数、日照时数比为输入参数的模型准确度最高，其次是以晴空指数为输入参数的模型，以日照时数比为输入参数的模型准确度最低.文献[18]以印度地区2013—2016年间的数据分析以晴空指数、日照时数比为输入参数的模型的准确度，得到了相同的结果.

图3 12个模型的误差评价指标

比较同类模型之间不同形式对准确度的影响. 对比以晴空指数为输入参数的模型(M1，M2，M3，M4，M5)可知：各个台站的R2由小到大的顺序依次是对数模型、指数模型、多项式模型，MAE、RMSE变化趋势相反；在多项式形式中，二次、三次的MAE、RMSE低于一次函数形式；二次、三次的MAE、RMSE、R2基本重合. 对比以日照时数比为输入参数以及以晴空指数、日照时数比为输入参数的双变量模型发现，多项式模型MAE、RMSE要低于其他形式，而在多项式模型中，二次、三次之间MAE、RMSE的差异较小.上述分析表明：多项式形式要优于对数、指数等形式；在多项式模型中，增加自变量的次数可以提高模型的准确度，但二次到三次不明显.

2.2 修正模型

由上述分析可知，模型的准确度随地区而变化，此外，模型的经验系数具有地区差异性，对于无辐射观测台站地区，其经验系数多少合适，是未知的，能否寻找经验系数与地理参数的关系，以解决无观测台站的辐射数据. 相关研究证实散射辐射的变化与纬度、经度、海拔高度密切相关，挑选单站点模型中准确度相对较高、形式相对简单的模型，即M1、M6、M10，利用多元回归分析建立模型经验系数与纬度(φ)、经度(β)、海拔高度(H)关系式，本文称为修正模型，以M13～M15表示，模型经验系数见表4.为便于和单站点模型对比，计算修正后的模型与单站点模型(即M13与M1、M14与M6、M15与M10)误差指标的差值以及在17个台站的平均值，结果如图4、5所示,以17个台站后4年的数据进行验证，计算3个模型的误差指标，如图6所示.

M13、M14、M15之间对比可知，除57494、59287、58362、57083台站外，M15的MAE、RMSE明显高于M13、M14，R2低于M13、M14，如图4所示；结合误差指标平均值进行分析(图6)，M15的MAE、RMSE最大，R2最小；结果表明，基于晴空指数、日照时数比模型的准确度低于基于晴空指数、日照时数比的单一输入参数模型.

计算修正前后模型之间误差指标(即M13与M1、M14与M6、M15与M10)的差值，并做17个台站的统计图，如图5所示，对于3类模型，大多数台站的单站点模型的MAE、RMSE要低于修正模型；基于晴空指数模型(M1、M13)的MAE、RMSE的差异最小，分别分布在-0.12%～0.08%、-0.1%～0.1%之间；基于晴空指数、日照时数比双变量模型(M10、M15)的MAE、RMSE的差异最大，分别分布在-10%～40%、-10%～50%之间；将空间扩展模型与单站点模型在17个台站的误差指标值进行平均，如图6所示，发现基于晴空指数修正模型(M13)的MAE、RMSE均低于单站点模型，R2高于单站点模型(M1)；而另两类修正模型的MABE、RMSE均高于单站点模型，分布高出0.48%、7.98%、4.80%、9.37%；上述对比说明，基于晴空指数修正模型的准确度稍高于单站点模型，其他修正模型准确度低于单站点模型.

表4 修正模型经验系数

图4 修正模型误差指标

图5 修正模型与单站点模型误差指标差值的统计图

Fig.5 Statistical chart of the difference of error index between modified and individual model

图6 修正模型与单站点模型误差指标的平均值

3 结 论

1)以中国17个辐射观测台站2000年—2017年的实测数据为基础，建立12个单站点模型，采用误差分析的方法对比不同模型的准确度；分析表明：多项式模型的准确度高于对数、指数形式；一次函数的准确度低于二次、三次函数形式，二次与三次函数无明显差异；以晴空指数、日照时数比为输入参数的模型最准确，基于晴空指数的模型次之，基于日照时数比的模型最差.

2)以单站点模型中表现较好形式简单的模型(M1、M6、M10)为基础，建立3个包含地理信息(纬度、经度、海拔高度)的修正模型，数据分析表明，基于晴空指数的模型最准确，基于日照时数比模型次之，基于晴空指数、日照时数比模型最差.

3)建议用本文提出的方法确定散射辐射模型的经验系数，采用基于晴空指数的一次函数模型计算散射辐射.