天空分布模型

本文讨论了确定天穹顶上日光和太阳辐射的发光分布所涉及的过程，以便准确解释建筑工地内部和周围的阴影和反射的影响。它着眼于所使用的数学模型及其实现，以及与在较长时间内聚合多个天空相关的一些问题。

介绍

从地面观察时，天空穹顶上的光线分布很少均匀。来自太阳的太阳辐射穿过地球大气层，可以被气体和尘埃颗粒以及水蒸气云吸收、反射和散射到各个方向。浑浊和浑浊也是高度动态的，当受到高层大气中的风的驱动时，可以在几分钟内改变整个天空的特征。图 1 显示了一些常见的不同类型的真实天空条件示例——在这种情况下是阴天、部分多云和晴朗的天空。

阴天（左）、部分多云（中）和晴（右）（
来源：http://www.philohome.com/skycollec/skycollec.htm）

图1 - 使用180°半球形“鱼眼镜头”记录的不同天空的一些真实示例。

CIE S 011/E：2003 和 ISO 15469：2004（E） 中描述的 CIE 标准通用天空提供了一个数学模型，该模型使用五 （5） 个系数和太阳在天空中的位置的组合来描述整个天空穹顶上日光的空间分布。它包括 15 种不同天空类型的一组系数，这些系数大致代表了具有不同云量和浑浊度的常见天空条件。但是，也可以直接从测量的瞬时天空照度或辐照度中得出系数值，或者独立调整它们以匹配任何天空条件。图 2 显示了如何使用特定系数值对较早使用的真实天空示例的辐射分布进行非常紧密的建模。

（来源： http://andrewmarsh.com/）

图2 - 与图 1 中的相同示例天空，使用 CIE 标准通用天空系数建模。

CIE 标准包括 15 种一般天空类型的系数设置，这些类型代表了一系列不同的天空条件 - 从阴天到万里无云，以及不同程度的阳光直射。还给出了额外的第 16 种类型，它基于 CIE 在当前标准之前定义和使用的更简单的阴天数学模型。图 3 显示了这 16 种天空类型中每一种的瞬时分布示例。

（来源： http://andrewmarsh.com/）

图3 - 16 种 CIE 标准通用天空类型的示例。

使用 CIE 标准 General Sky

CIE Standard General Sky 模型的核心是基于三个关键特征的变化;天空晴朗度、天空亮度和太阳高度。清晰度和亮度可以通过比较直接法线和漫反射水平辐照度和/或照度得出。太阳高度可以针对特定日期和时间进行计算，也可以作为水平面上方的测量角度给出。然后可以直接从这三个值计算五个天空系数。

使用计算出的天空系数，可以相对于天顶照度计算出天穹上任何点的亮度。鉴于 EnergyPlus 天气文件格式 （EPW） 包括天顶照度的每小时值，人们可能会认为可以直接且相对简单地使用此数据。然而，在许多EPW文件中，天顶照度的值是不可靠的，在太阳高度相对较高的中午被记录为零。图 4 显示了纽约市、新加坡和澳大利亚珀斯的 EPW 数据中此问题的示例。图表中心的蓝色/紫色区域显示与同一数据集中的全局水平照度值相比，天顶照度应真正达到峰值的零值区域。

图片 4 - 从 EnergyPlus 天气数据文件中获取的年度每小时天顶照度示例图，在中午显示不正确的零值。

值得庆幸的是，EPW 文件确实包含更可靠的年度每小时直接法线和漫射水平照度和辐照度数据，可以使用这些数据来代替。有趣的是，根据Perez et.al 的说法，天空晴朗度和亮度的计算实际上取决于漫反射水平辐照度而不是照度。1990 （参考文献： 1），从照度转换为辐照度有时可能涉及少量迭代以收敛到正确的值。因此，使用来自天气数据的辐照度数据实际上使亮度计算速度略快。

此外，需要注意的是，CIE Standard General Sky 仅模拟日光的空间分布，而不是阳光。日光是指从整个天穹反射和折射的全向漫射光，而太阳光仅指直接从太阳直接穿过大气层的高定向太阳辐射光线。

获取绝对值

CIE Standard General Sky 模型计算的值始终相对于天空天顶的照度给出。这意味着分布在天空中的值的尺度和范围是非常任意的，并且会随着不同的天空类型和太阳高度而发生显着变化。要将这些相对值转换为照度或辐照度的绝对值，有必要将整个天幕的总贡献归一化，使其等于 1，然后再乘以漫反射水平辐照度或照度的参考值。

这通常意味着将天穹细分为一系列离散的天空斑块，然后将它们的所有值进行积分以获得总贡献。由于参考值是漫反射水平辐照度或照度，因此在确定每个斑块的贡献时，积分需要考虑朗伯余弦定律和水平面上每个天空斑块的光入射角。此外，许多天空细分技术不会产生大小完全相同的天空斑块。因此，每个贴片的贡献也必须精确地进行区域加权，如下式所示。

哪里：

• F规范= 应用于所有修补程序的归一化因子。
• n = 天穹被划分为的天空斑块的总数。
• Lp= 与天顶相比，每个天空斑块的相对贡献。
• 一个p= 每个天空斑块的相对面积占天空穹顶总面积的一小部分。
• δ = 天空斑块中心与天空天顶之间的角度。
• L禅= 天空天顶的相对贡献。

一旦每个天空斑块的值被归一化和面积加权，就可以通过简单地乘以相应的瞬时漫反射水平辐照度或照度值来确定它们的绝对辐射度或亮度值。这通常作为最后的单独步骤完成，以允许缓存和重用计算出的天空补丁数据。

需要注意的是，当太阳在日出和日落时或接近天空时非常低时，CIE 标准通用天空模型的准确性会降低。当太阳在地平线上方的高度小于2度时，大部分亮度发生在接近水平面的掠入射处。这可能会导致非常小的归一化因子来测试 Javascript 引擎的浮点精度。然而，一旦理解了这一点，就可以相对容易地通过在较低的太阳高度进行值缩放来适应，并且任何残余效应都变得相对微不足道，因为漫射水平辐照度和照度在这些时间左右都收敛到零。

聚集和积累天空

为了计算长时间的值，可以生成许多连续的瞬时天空分布，对每个分布运行计算，然后对结果求和。根据计算的性质，有时生成单个天空分布（即许多连续瞬时天空分布的累积总和）会更有效，然后使用该累积天空运行一次所需的太阳计算。

生成和聚合多个天空分布是一项相对简单的任务，尤其是当它们都基于相同的天空斑块时。然而，使用每小时的天气数据这样做可能会导致重大的潜在误差，因为太阳在天空中的实际路径是平滑和连续的，而每小时采样会导致太阳位置之间的巨大差距，如图5所示。

（来源： http://andrewmarsh.com/）

图片 5 - 两个示例天空分布的动画，相隔一小时。

亚小时采样是解决方案，但这需要谨慎地对亚小时天气数据值进行插值。

在天空条件之间插值

亚小时采样的明显方法是简单地计算两个小时的天空条件，然后在每个天空段中的两个值之间线性插值。然而，如图 6 所示，这种方法只会导致太阳从前一小时的位置淡出，并在下一小时的新位置淡入。在两个小时之间插值的天空分布并不能准确表示当时实际发生的真实天空条件。

（来源： http://andrewmarsh.com/）

图6 - 插值各个天空段值只会使太阳从前一小时淡出，并在下一小时淡入。

另一种方法是在每个小时边界处对每小时直射太阳辐照度值和漫射太阳辐照度值进行线性插值，并得出每个过渡时间步的天空系数。然而，从直接和漫射太阳辐照度得出的天空系数是基于表格数据范围的。这意味着，在一个小时内从一组漫射/直接值过渡到另一组漫射/直接值，有时会在天空分布的特征中引入明显的步进或跳跃，从而导致包含可能在记录数据的一小时内实际没有发生的中间天空条件。

另一种方法

图 7 展示了一个示例，清楚地显示了在表格数据范围内越过边界阈值时发生的突然阶跃变化。在所示示例中，漫反射水平照度在 6000 和 20000 lux 之间插值，并显示在 10000 和 15000 lux 处发生的大步进。在这些阈值下，您可以清楚地看到一些天空系数滑块突然跳到明显不同的值。

（来源： http://andrewmarsh.com/）

图片 7 - 在不同的直接/漫射值之间进行插值时可能导致的天空条件突然发生阶跃变化的示例。

正确的方法

为了避免上述两个问题，作者找到的最佳解决方案是计算每个小时天空条件的结果天空系数，然后在两组计算的系数之间线性插值 - 在每个过渡时间步生成亚小时太阳位置和结果天空条件。这种方法可以在任何两个已知的天空条件之间实现最平滑和最真实的过渡。图 8 显示了使用此方法在两个小时天空条件之间产生的插值示例。

（来源： http://andrewmarsh.com/）

图8 - 通过插值两组每小时的天空系数并生成临时天空，在两个已知的每小时天空条件之间产生亚小时过渡。

分配阳光直射

为了正确模拟真实的天空，模型中还必须包含直接太阳分量（通常称为太阳光）。起初，将阳光直射分配给天空分布似乎很简单，但是一旦将天空划分为离散的斑块，就会变得复杂。

从文献中，有两种常用的方法：

1. 太阳圆盘法：
   仅将阳光分布到直径为0.533°的太阳圆盘实际重叠的天空斑块上，按受影响的斑块数量及其相对重叠面积进行分配。

2. 最近邻法：
   将太阳光分布到四 （4） 个最近的相邻天空斑块，根据每个斑块中心到太阳中心的相对距离进行分配。

虽然这两种方法的绝对精度主要取决于天空细分的分辨率，但第一种方法中太阳盘的小尺寸有时会导致值在斑块之间转换时出现明显的跳跃，尤其是在使用少量大型天空斑块时。这种情况会导致一段时间的静态天空状况，因为当太阳缓慢穿过它时，所有的阳光都被分配到一个斑块上，当它移动到下一个斑块时，会穿插突然的变化。图 9 显示了两种方法在同一时期和同一天的可视化比较示例。

方法一：太阳盘。

方法2：最近邻。

图片 9 - GIF动画，显示每种方法之间差异的示例。这里使用了 5° 等面积的天空细分来说明这种效果。

最近邻方法可能看起来不太准确，但实际上，当太阳在天空中移动时，它提供了更平滑和更连续的值过渡，并更好地模拟了实际的太阳位置。当对窗口和光圈使用预先计算的阴影遮罩时，这一点很重要，在这些遮光罩中，一个天空斑块可能可见，但另一个紧邻的天空斑块可能不可见。使用太阳盘方法将导致在单个时间步长的空间内在无阳光直射和完全阳光直射之间突然跳跃。使用最近邻方法将导致内部光照条件发生更现实的渐进变化。

日光分布方法的选择也会影响累积天空计算中的聚合值。图 10 显示了使用两种阳光直射分布方法对美国纽约的年累积天空的比较。每种方法的整体模式非常相似，但是最近邻方法的拖尾效果是显而易见的。两种方法的总体天空总辐射度相同，但比例略有不同，因为任何特定天空段的最大值都会受到所选方法的影响。最好使用哪种方法是一个公开的争论。

方法一：太阳盘。

方法2：最近邻。

图片 10 - 比较两种方法的纽约年度累计总天空量。

低于小时时间步长

下一个显而易见的问题是，应该使用什么亚小时时间增量来生成累积天空分布。如图 11 所示，答案绝对不是 60 分钟，甚至不是 30 分钟。为了获得平滑准确的累积天空分布，需要 10-12 分钟或更短的亚小时时间增量。

（来源： http://andrewmarsh.com/）

图片 11 - 不同次小时采样时间增量对所得累积天空分布的影响。

在大约 13 分钟时，由于太阳位置的跳跃而导致的变化开始在视觉上变得明显，随着时间增量的增加，服务器性也会增加。此外，1 分钟和 5 分钟的增量天空分布之间几乎没有差异，这表明 5 到 10 分钟之间的某个地方是计算速度和精度之间的最佳权衡。

---

引用

1. Perez R、Ineichen P、Seals R、Michalsky J 和 Stewart R. （1990），从直接辐照度和全球辐照度模拟日光可用性和辐照度分量，太阳能 44（5）：271-289。