太阳辐射度计算公式

光伏太阳能系统——太阳光源 太阳辐射几何学 根 据 地 球 自 转 和 围 绕 太 阳 公 转 的 日 常 以 及 季 节 性 变 化 ， 我 们 推 导 了 在 任 意 朝 向 的 光 伏 阵 列 上 的 太 阳 辐 射 强 度 计 算 方 法 。 为 了 介 绍 这 个 计 算 方 法 ， 需 要 引 入 两 组 描 述 太 阳 和 地 球 相 对 位 置 的 坐 标 系 。 第一组给出太阳相对于固定在地球上的坐标系的位置，坐标系的一个轴指向天顶，另 一个与之正交的轴指向地平线，即下图中的i, j坐标系，称为地平坐标系。 另一组坐标系也固定在地表的同一个位置，但是一个轴指向极点，即北天极，另一个 与之正交的轴则指向与赤道平行的方向，即下图中的I, J坐标系，称为时角坐标系。 两个坐标系之间有如下关系： i  (sin )I  (cos ) J j  (cos )I  (sin )J k  K 其中的  是观测点的纬度。 下图给出了正午时分的太阳（在天空的最高点）在地平坐标系中的方位。   其中 z 称为天顶角而a z 称为太阳方位角。从而，指向太阳的单位矢量可以写成： s  (cos z)i  (sin z cos az ) j  (sin z sin az )k 第二个坐标系的一个轴始终指向北天极而另一轴位于赤道面上，方向指向南点。因而， 指向正午太阳的单位矢量可以写成： s  (sin )I  (cos cos)J  (cos sin)K 由于在两组坐标系中的矢量s相同，我们可以得到： cos z  (sin )(sin  )  (cos )(cos  ) cos (sin z)(cos az )  (cos )(cos) sin   (cos  ) sin (sin z) sin az  cos (sin) 如果一年中的天数确定了，就可以按照库珀公式计算太阳的赤纬角  ：   23.45 sin 360  284  n   365  如果观测者所在的地理纬度  也确定了，那么就可以定义太阳的时角 ，它反映了从日 出到日落时分太阳的天顶角在一天中的变化。从而，一天中的时间变化也可以用太阳时角 的变化来表示。在日出或日落时分，太阳的时角正好等于 90 度或  2 ，即太阳处在地平线 上。 设cos z = cos 90 = 0，则有: sunrise / sunset  ss  ar cos tan tan   下面再来推导太阳方位角和时角之间的关系。 首先，天顶角的余弦函数cos z反映了太阳的辐射强度投射到地平表面上的光照比。如 果我们把太阳的辐射想象成从太阳流向地球的流体，那么这束流体投射到有某个偏转角的 表面上时，穿过这个表面的流量和偏转角的余弦成正比，如下图所示。 另一种理解，也可以把受到辐射的面积想象成垂直于来流方向，而流体是以某个偏转 角投射到这个表面上的，当然也等于Acos z。 参考前面的内容，cos z的表达式可以写成： cos z  C1  C2 cos 上式中的余弦函数cos z主要依赖于C 1的符号，因为 C2总是正的。 照射在倾斜阵列上的光强 我们总能在所处的地理纬度上找出比倾斜布置好得多的光伏阵列布置方式。为了获得 照射在倾斜阵列上的光强，我们需要在当地坐标系下计算光伏阵列的方位。 在下图所示的坐标系中，垂直于一个倾斜布置的光伏阵列的单位矢量可以写成： nc  (cos  )i  (sin  )(cos ) j  (sin  )(sin )k ext 对于那些没有追踪功能的固定阵列，它们通常是面向正南的。垂直于阵列的单位矢量 可以简化成： nc  (cos  )i  (sin  ) j 照射到阵列上的光强与指向太阳的单位矢量s和垂直于阵列的单位矢量n c之间的夹角 的 余 弦 成 正 比 ， 这 个 量 可 以 通 过 求 两 个 矢 量 的 标 量 积 得 到 。 使 用 上 面 给 出 的 那 些 关 系 式 ， 这 个 余 弦 可 以 写 成 ： cos(nc , s)  sin sin(   )  cos cos(   ) cos 这表明，对于指定的赤纬角，地理纬度和光伏阵列的倾斜角是一对余弦变换。 月平均日辐射量和其它参数 假设你处在外太空，并且运行在地球轨道上，而你的太阳光伏阵列是水平的。也就是 说，光伏阵列是平躺在你所作运动的轨道面上。而且在你和太阳之间没有任何东西，哪怕 是 一 丝 以 太 ， 来 削 弱 太 阳 的 辐 射 。 此 时 的 太 阳 辐 射 强 度 我 们 称 为 大 气 外 界 辐 射 强 度 ， 等 于 ：   1.350kw / m 2 在一个昼夜的过程中，投射到水平安置的光伏阵列上的所有太阳辐射等于： ss 24  2ext cos z  2 d  H ext 0   使用前面提到的cos z的计算式，可以积分得到： Hext  24  ext cos cossinss ss cosss  ssss ss 假设你还在地球的大气层以外，同时你的光伏电池阵列相对于水平面倾斜了一个角度  ，但是面向正南，也就是说 =0。那么，垂直于光伏阵列的单位矢量和指向太阳的单位 矢量之间的夹角的余弦等于： cos(nc , s)  sin sin(   )  cos cos(   ) cos 这个余弦当然是太阳时角的函数，或者说是一日当中的时间的函数。 将这个余弦函数对时间积分，从日出积分到日落，可得： Hext ( )  24   ext cos cos(   )sin cos  在 此 ， 我 们 假 设 日 落 时 分 的 太 阳 方 位 角 与 光 伏 阵 列 平 行 ， 即 假 设 cos(nc,s) = cos 90 = 0， 因而有：  ar cos tan  tan   晴朗指数 晴朗指数 K T 是一个经验参数，可以根据在某一固定位置测量得到的日照强度计算得 到。它的定义是地表水平面上的月平均日太阳辐射量与大气层外界水平面上的月平均日太 阳辐射量之比值。 KT  H / H ext 下表是一些不同地点的晴朗指数随月份的变化。  ss  倾斜面上的天空散射辐射 有 了 当 地 的 晴 朗 指 数 ， 以 及 大 气 层 外 界 水 平 面 上 的 辐 射 强 度 计 算 方 法 ， 我 们 可 以 估 算 出 地 面 倾 斜 布 置 的 光 伏 阵 列 上 的 月 平 均 日 辐 射 总 量 。 这 里 所 用 的 “水 平 ”晴 朗 指 数 一 般 是 根 据 当 地 收 集 到 的 太 阳 辐 射 数 据 计 算 出 来 的 。 因 而 ， 它 包 含 了 天 气 的 反 常 、 太 阳 的 直 射 辐 射 、 天 空 的 散 射 辐 射 、 以 及 反 射 的 太 阳 辐 射 在 内 的 所 有 辐 射 在 水 平 面 上 的 总 量 。 如 果 我 们 将 光 伏 电 池 倾 斜 一 个 角 度 （ 即 前 面 图 示 的  角），那么我们必须考虑到天空的散射辐射和 反 射 的 太 阳 辐 射 的 方 式 都 发 生 了 变 化 。 换 言 之 ， 我 们 需 要 构 造 一 个 适 用 于 当 地 的 倾 斜 表 面 系 数 Kt 。 我们在这儿使用的方法是把倾斜布置的光伏阵列上的月平均日辐射总量和倾斜放置的 光伏阵列上的月平均日辐射总量拿一个因子R 联系起来。这个R因子的计算方法最早由Liu and Jordan1提出，后来经过了Klein 3的评估，随后又经 Klein and Theilacker4修改和提炼。 我们在此给出它的一阶近似表达式： H   R  H 其中 R    Hd   Hd   1  cos    1  cos   1   H   Rb   H      2 2         这两个表达式中的前一个表示倾斜放置的光伏阵列上的月平均日辐射总量可以用倾 斜 布 置 的 光 伏 阵 列 上 的 月 平 均 日 辐 射 总 量 和 一 个 因 子 R相 乘 得 到 。 而 这 个 因 子 R的计算方法  ss 由第二个式子给出。至此，我们的计算都是稳妥而保守的。 不过，这个R因子的计算还是有趣的。在其计算式中，  Hd  是倾斜布置的光伏阵列上   H  的月平均日散射辐射量和同一光伏阵列上的月平均日辐射总量的比值。算式中的R b则是倾 斜布置的光伏阵列上的月平均直射辐射量和水平布置的光伏阵列上的月平均直射辐射量的 比值。第一个比值是根据实验数据构造的，而第二个比值则可以根据我们前面学过的辐射 几何学知识分析得到。参数 表示当地的地面反射系数，最好是从当地测量获得。 让 我 们 首 先 考 虑  等于零的情况，即光伏阵列是水平布置的。那么，在我们的计算式 中 R＝ Rb。 这 意 味 着 ， 在 计 算 R的 式 子 的 右 边 第 一 项 ， 我 们 “收 回 了 ”散 射 辐 射 对 水 平 布 置 的 光 伏 阵 列 上 的 总 辐 射 量 的 贡 献 ， 使 第 一 项 只 表 示 直 射 辐 射 的 贡 献 。 然 后 ， 又 在 第 二 项 中 把 散 射 辐 射 的 贡 献 加 回 去 。 Rb因 而 可 以 用 大 气 层 外 界 的 倾 斜 表 面 上 的 日 辐 射 总 量 （以 H ext ( ) 表示） 和水平表面上 的日辐射总量（以 H ext 表示）的比值来近似。 H ( ) cos cos(   )sin   cos  R  ext b ss   ss ssH ext cos cos sin ss  ss cos ss 上式中的  是倾斜布置的光伏阵列在日落时分的太阳时角，  则是水平布置的光伏 ss 阵列在日落时分的太阳时角。 我们不得不回过头来认真考虑散射辐射的贡献。如前所述，我们主要依靠实验数据。 目前的文献中至少已经有两条曲线，如下图所示。