太阳辐射度计算公式

光伏太阳能系统——太阳光源 太阳辐射几何学 根据地球自转和围绕太阳公转的日常以及季节性变化，我们推导了在任意朝向的光伏 阵列上的太阳辐射强度计算方法。为了介绍这个计算方法，需要引入两组描述太阳和地球 相对位置的坐标系。 第一组给出太阳相对于固定在地球上的坐标系的位置，坐标系的一个轴指向天顶，另 一个与之正交的轴指向地平线，即下图中的i, j坐标系，称为地平坐标系。 另一组坐标系也固定在地表的同一个位置，但是一个轴指向极点，即北天极，另一个 与之正交的轴则指向与赤道平行的方向，即下图中的I, J坐标系，称为时角坐标系。 两个坐标系之间有如下关系： JIi )(cos)(sin λλ += JIj )(sin)cos( λλ +−= Kk = 其中的λ是观测点的纬度。 下图给出了正午时分的太阳（在天空的最高点）在地平坐标系中的方位。 其中z称为天顶角而 z a称为太阳方位角。从而，指向太阳的单位矢量可以写成： kazjazizs zz )sin(sin)cos(sin)(cos −+= 第二个坐标系的一个轴始终指向北天极而另一轴位于赤道面上，方向指向南点。因而， 指向正午太阳的单位矢量可以写成： KJIs )sin(cos)cos(cos)(sin ωδωδδ −+= 由于在两组坐标系中的矢量s相同，我们可以得到： ωλδλδ cos))(cos(cos))(sin(sincos +=z δλλωδ sin)(cossin))(cos(cos))(cos(sin −= z az )(sincossin)(sin ωδ= z az 如果一年中的天数确定了，就可以按照库珀公式计算太阳的赤纬角δ：       + ⋅⋅= 365 284 360sin45.23 n δ 如果观测者所在的地理纬度λ也确定了，那么就可以定义太阳的时角ω，它反映了从 日出到日落时分太阳的天顶角在一天中的变化。从而，一天中的时间变化也可以用太阳时 角的变化来表示。在日出或日落时分，太阳的时角正好等于90度或 2 π ，即太阳处在地平 线上。 设cos z = cos 90 = 0，则有: []λδωω tantancos / −== ar sssunsetsunrise 下面再来推导太阳方位角和时角之间的关系。 首先，天顶角的余弦函数cos z反映了太阳的辐射强度投射到地平表面上的光照比。如 果我们把太阳的辐射想象成从太阳流向地球的流体，那么这束流体投射到有某个偏转角的 表面上时，穿过这个表面的流量和偏转角的余弦成正比，如下图所示。 另一种理解，也可以把受到辐射的面积想象成垂直于来流方向，而流体是以某个偏转 角投射到这个表面上的，当然也等于Acos z。 参考前面的内容，cos z的表达式可以写成： ωcoscos 21 CCz += 上式中的余弦函数cos z主要依赖于C 1 的符号，因为C 2 总是正的。 照射在倾斜阵列上的光强 我们总能在所处的地理纬度上找出比倾斜布置好得多的光伏阵列布置方式。为了获得 照射在倾斜阵列上的光强，我们需要在当地坐标系下计算光伏阵列的方位。 在下图所示的坐标系中，垂直于一个倾斜布置的光伏阵列的单位矢量可以写成： kjin c ))(sin(sin))(cos(sin)(cos γβγββ −+= 对于那些没有追踪功能的固定阵列，它们通常是面向正南的。垂直于阵列的单位矢量 可以简化成： jin c )(sin)(cos ββ += 照射到阵列上的光强与指向太阳的单位矢量s和垂直于阵列的单位矢量n c 之间的夹角 的余弦成正比，这个量可以通过求两个矢量的标量积得到。使用上面给出的那些关系式， 这个余弦可以写成： ωβλδβλδ cos)cos(cos)sin(sin),cos( −+−=sn c 这表明，对于指定的赤纬角，地理纬度和光伏阵列的倾斜角是一对余弦变换。 月平均日辐射量和其它参数 假设你处在外太空，并且运行在地球轨道上，而你的太阳光伏阵列是水平的。也就是 说，光伏阵列是平躺在你所作运动的轨道面上。而且在你和太阳之间没有任何东西，哪怕 是一丝以太，来削弱太阳的辐射。此时的太阳辐射强度我们称为大气外界辐射强度，等于： 2 /350.1 mkw ext =φ 在一个昼夜的过程中，投射到水平安置的光伏阵列上的所有太阳辐射等于： extext Hdz ss ≡       ⋅ ∫ ω π φ ω 0 2 24 cos2 使用前面提到的cos z的计算式，可以积分得到： [] ssssssextext H ωωωλδφ π cossincoscos 24 −= 假设你还在地球的大气层以外，同时你的光伏电池阵列相对于水平面倾斜了一个角度 β，但是面向正南，也就是说γ =0。那么，垂直于光伏阵列的单位矢量和指向太阳的 矢量之间的夹角的余弦等于： ωβλδβλδ cos)cos(cos)sin(sin),cos( −+−=sn c 这个余弦当然是太阳时角的函数，或者说是一日当中的时间的函数。 将这个余弦函数对时间积分，从日出积分到日落，可得： [ ] cossin)cos(cos 24 )( ssssssextext H ωωωβλδφ π β −−= 在此，我们假设日落时分的太阳方位角与光伏阵列平行，即假设cos(n c ,s) = cos 90 = 0， 因而有： []λδω tantancos −= ar ss 晴朗指数 晴朗指数 T K是一个经验参数，可以根据在某一固定位置测量得到的日照强度计算得 到。它的定义是地表水平面上的月平均日太阳辐射量与大气层外界水平面上的月平均日太 阳辐射量之比值。 extT HHK /= 下表是一些不同地点的晴朗指数随月份的变化。 倾斜面上的天空散射辐射 有了当地的晴朗指数，以及大气层外界水平面上的辐射强度计算方法，我们可以估算 出地面倾斜布置的光伏阵列上的月平均日辐射总量。这里所用的“水平”晴朗指数一般是 根据当地收集到的太阳辐射数据计算出来的。因而，它包含了天气的反常、太阳的直射辐 射、天空的散射辐射、以及反射的太阳辐射在内的所有辐射在水平面上的总量。如果我们 将光伏电池倾斜一个角度（即前面图示的β角），那么我们必须考虑到天空的散射辐射和 反射的太阳辐射的方式都发生了变化。换言之，我们需要构造一个适用于当地的倾斜表面 系数 t K。 我们在这儿使用的方法是把倾斜布置的光伏阵列上的月平均日辐射总量和倾斜放置的 光伏阵列上的月平均日辐射总量拿一个因子R联系起来。这个R因子的计算方法最早由Liu and Jordan 1 提出，后来经过了Klein 3 的评估，随后又经Klein and Theilacker 4 修改和提炼。 我们在此给出它的一阶近似表达式： HRH ⋅= β 其中       − ⋅+       + ⋅         +⋅               −= 2 cos1 2 cos1 1 β ρ β H H R H H R d b d 这两个表达式中的前一个表示倾斜放置的光伏阵列上的月平均日辐射总量可以用倾 斜布置的光伏阵列上的月平均日辐射总量和一个因子R相乘得到。而这个因子R的计算方法 由第二个式子给出。至此，我们的计算都是稳妥而保守的。 不过，这个R因子的计算还是有趣的。在其计算式中，         H H d 是倾斜布置的光伏阵列上 的月平均日散射辐射量和同一光伏阵列上的月平均日辐射总量的比值。算式中的R b 则是倾 斜布置的光伏阵列上的月平均直射辐射量和水平布置的光伏阵列上的月平均直射辐射量的 比值。第一个比值是根据实验数据构造的，而第二个比值则可以根据我们前面学过的辐射 几何学知识分析得到。参数ρ表示当地的地面反射系数，最好是从当地测量获得。 让我们首先考虑β等于零的情况，即光伏阵列是水平布置的。那么，在我们的计算式 中R＝R b 。这意味着，在计算R的式子的右边第一项，我们“收回了”散射辐射对水平布置 的光伏阵列上的总辐射量的贡献，使第一项只表示直射辐射的贡献。然后，又在第二项中 把散射辐射的贡献加回去。 R b 因而可以用大气层外界的倾斜表面上的日辐射总量（以)(β ext H表示）和水平表面上 的日辐射总量（以 ext H表示）的比值来近似。 [ ] [] ssssss ssssss ext ext b H H R ωωωλδ ωωωβλδβ cossincoscos cossin)cos(cos)( − −− == 上式中的 ss ω是倾斜布置的光伏阵列在日落时分的太阳时角， ss ω则是水平布置的光伏 阵列在日落时分的太阳时角。 我们不得不回过头来认真考虑散射辐射的贡献。如前所述，我们主要依靠实验数据。 目前的文献中至少已经有两条曲线，如下图所示。 参考文献 1. Liu, B.Y.H & Jordan, R.C., “The interrelationship and characteristic distribution of direct, diffuse, and total solar radiation”, Solar Energy 4(3) 1960. 2. Liu, B.Y.H & Jordan, R.C., “Daily insolation on surfaces tilted toward the equator”, Trans. ASHRAE 526-541, 1962. 3. Klein, S.A. “Calculation of Monthly Average Insolation on Tilted Surfaces”, Solar Energy, vol. 19, 325-329, 1977. 4. Klein. S.A. & Theilacker, J.C., “An Algorithm for Calculating Monthly-Average Radiation on Inclined Surfaces”, Jour. of Solar Energy Engineering, vol. 103, 29-33, 1981.