1. 引言
我国正努力实现“双碳”目标,其中大力推进具有环保和低碳排放特点的新型清洁能源技术一一塔式太阳能光热发电是一项核心举措,而定日镜作为塔式电站光热发电系统中的重要构件,其聚焦零落的太阳光线并将其传导到吸收塔中的功能,大大提高了太阳能光热系统的光学功率和能源密度。因此,研究定日镜的性质和配置排列方式有助于实现能源高效转换与充分利用,并为能源的可持续发展做出贡献。
本文将基于所获取的数据,综合运用反射定律、坐标转化、能流密度分布等数学物理知识,建立空间模型以求解定日镜场的光学效率和输出热功率,并从定日镜的排布方式、吸收塔的阴影域、光学效率的呈现分布三个角度考虑,利用遗传算法优化定日镜场的排列设计,使得达到额定功率的情况下,单位镜面面积年平均输出热功率最大。
其中,阴影挡光损失的计算通常涉及到复杂的几何关系,例如太阳位置、定日镜的几何形状等,本文通过坐标转换,将坐标系转换为统一的标准,对不同定日镜上的光线遮挡情况进行有效的刻画和评估。此外,本文使用的遗传算法在优化模型中具有较高的适配性和解决能力,能够较好地处理高维参数且出现多种局部最优解的搜索空间,且在搜索过程中同时考虑多个解,增加算法对噪声和随机扰动的鲁棒性。
2. 定日镜场的概念模型建立
定日镜是塔式太阳能光热发电站(以下简称塔式电站)收集太阳能的基本组件,其底座由纵向转轴和水平转轴组成,平面反射镜安装在水平转轴上。纵向转轴的轴线与地面垂直,可以控制反射镜的方位角。水平转轴的轴线与地面平行,可以控制反射镜的俯仰角,两转轴的交点(也是定日镜中心)离地面的高度称为定日镜的安装高度。定日镜将太阳光反射汇聚到安装在镜场中吸收塔顶端上的集热器,加热其中的导热介质,并将太阳能以热能形式储存起来,再经过热交换实现由热能向电能的转化。太阳光并非平行光线,而是具有一定锥形角的一束锥形光线,因此太阳入射光线经定日镜任意一点的反射光线也是一束锥形光线 [1] 。定日镜及底座示意图如图1:
Figure 1. Schematic diagram of the heliostat field
图1. 定日镜场示意图
本文假设在中心位于东经98.5˚,北纬39.4˚,海拔3000 m,半径350 m的圆形区域内建设一个圆形定日镜场。以圆形区域中心为原点,正东方向为𝑥轴正向,正北方向为𝑦轴正向,垂直于地面向上方向为z轴正向建立坐标系,称为镜场坐标系。规划的吸收塔高度为80 m,集热器采用高8 m、直径7 m的圆柱形外表受光式集热器。吸收塔周围100 m范围内不安装定日镜,留出空地建造厂房,用于安装发电、储能、控制等设备。此外,镜面边长在2 m 至8 m之间,安装高度在2 m至6 m之间,安装高度保证镜面在绕水平转轴旋转时不会触及地面。为简化计算,本文中所有“年均”指标的计算时点均为当地时间每月21日9:00、10:30、12:00、13:30、15:00。
3. 基于光学性能理论的镜场输出热功率测算
3.1. 坐标系的建立与变换
为了研究太阳入射光线在每个平面定日镜上的反射过程,本文需要建立定日镜和采光点的局部坐标系,以便后续通过向量的运算确定阴影的遮挡损失情况等。
由线性代数知识,在给定的线性空间中,任意两个笛卡尔坐标系之间只需通过平移和旋转即可完成转换。以地面坐标系作为参考坐标系
,其原点位于东经98.5˚,北纬39.4˚,海拔3000 m处。镜场中任一定日镜的局部坐标系
的3个坐标轴在
中的方向余弦分别为
、
和
。因此,假设镜面坐标系中某离散点的坐标为
,则它在地面坐标系
中的坐标为
,满足:
(1)
其中,
是该镜面坐标系
的原点在地面坐标系
中的坐标。上式也可表示成:
(2)
此处,
是两坐标系之间的旋转矩阵。由旋动理论的旋转变换关系Rodrigues等式 [2] ,存在一个旋转向量
,由其构成的关系矩阵
为反向斜对称矩阵,其结构为:
(3)
且旋转矩阵
与
之间满足下述关系式:
(4)
考虑到
是奇异矩阵,本文可以任取空间中的两个向量,已知它们在地面坐标系
和镜面坐标系
下的坐标,通过构造增广矩阵方程即可求得两坐标系之间的旋转向量
。此外,还需将旋转矩阵
中的3个列向量归一化,使之成为正交矩阵 [3] ,具体示意如图2:
Figure 2. Schematic diagram of coordinate system transformation
图2. 坐标系转化示意图
3.2. 镜场年平均光学效率的推算
不妨假设,该吸收塔位于圆形定日镜场的中心,定日镜尺寸均为6 m × 6 m,安装高度均为4 m,且给定所有定日镜中心的位置,则该定日镜场的年平均光学效率
满足下式:
(5)
其中,对任一特定时刻(每月21日9:00、10:30、12:00、13:30、15:00),每个定日镜的光学效率表示为:
(6)
它包括了阴影挡光效率
、余弦效率
、大气透射率
、集热器的截断效率
和镜面反射率
。下文将逐一对其进行分析。
其中,定日镜的光学效率示意图如图3:
Figure 3. Schematic diagram of the optical efficiency of a heliostat
图3. 定日镜的光学效率示意图
3.2.1. 阴影挡光效率
通过分析可知其基本表达式为:
(7)
阴影遮挡损失是指:在太阳能镜场设置中,若一个定日镜的阴影投射到另一定日镜的反射面,或者其反射的光束照射到另一定日镜的背部,就会导致整个镜场接收到的太阳光线减少 [4] 。而一般定义下,阴影遮挡损失应考虑吸收塔对镜场造成的阴影损失、前排定日镜对后排定日镜的入射太阳光线造成的阴影损失、后排定日镜的出射光线被前排定日镜遮挡造成的挡光损失共三部分。
由于吸收塔周围100 m范围内不安装定日镜,因此需要由计算吸收塔在地面的投影面积。已知太阳高度角
的计算公式 [5] :
(8)
其中,
为当地纬度(北纬为正),
为太阳时角,满足:
(9)
其中,
为当地时间,
为太阳赤纬角,满足:
(10)
其中,
为以春分作为第0天起算的天数。
将太阳高度角
代入下式:
(11)
其中,
为吸收塔高度,
为其影子长度。
计算可得,恰好覆盖到定日镜安置范围100 m处的临界值
约为0.6247。即当太阳高度角小于
时,部分定日镜会受到吸收塔的阴影遮挡,通过等比例放缩可求得吸收塔对镜场造成的阴影损失面积满足:
(12)
其中,
为阴影损失面积,
为吸收塔的直径。
在此基础上,还需计算定日镜重合分布造成的相邻镜面间的遮挡损失,如图4所示:
Figure 4. Schematic diagram of the shadow blocking of the front and rear heliostats
图4. 前后定日镜的阴影挡光示意图
探究阴影挡光损失的过程,实则是计算通过定日镜A上任一点的入射光线或反射光线,沿着相反的方向是否会落入其相邻定日镜B所在的区域,并利用坐标变换求出该光线在B镜面坐标系下的坐标 [1] 。
假设某光线(入射或反射)在地面坐标系
下的坐标为
,该光线在镜面局部坐标系下的坐标为
,则利用两坐标系间的旋转矩阵
有:
(13)
若该光线在定日镜A坐标系上的采光点坐标为
,设其在地面坐标系
下的向量坐标为
,转换到定日镜B坐标系上的坐标为
,则根据两点一线原理可求得该光线与定日镜B之间的交点。具体步骤如下:
(14)
其中,
为镜A的坐标原点在
下的向量坐标,同理有:
(15)
其中,
为镜B的坐标原点在
下的向量坐标。
将
坐标系下的光线向量
转换到镜面坐标系B,得到坐标
表示为:
(16)
联立直线方程即可得到落入定日镜B中的向量坐标
。
解:
(17)
得:
(18)
由此,即可判断
是否落入镜面B的镜面区域。
3.2.2. 余弦效率
由于镜面法线与太阳入射光线有一定角度
,所以会造成余弦损失,定义余弦效率
,其物理意义为表示定日镜面实际接受光照的面积和总面积的比值,对于单个定日镜而言,余弦效率随角度
的减小而增大 [6] 。其具体计算方法如下:
对于定日镜面法线与太阳入射光线夹角
的求解,根据几何光学的知识,反射角与入射角是相等的,其值均为
,因此可以将问题转化为研究入射光线和反射光线。
首先,太阳入射光线的方向向量
可以表示为
(19)
其中
为太阳高度角,
为太阳方位角,且已知太阳方位角余弦值的计算公式为
(20)
其中
为太阳赤纬角,
为当地纬度。
对于反射光线,其方向由定日镜面中心指向集热器中心,而定日镜面中心坐标为
,集热器中心坐标为
,因此反射光线的方向向量
可以表示为
(21)
其中
为吸收塔高度,
为定日镜安装高度。于是根据
(22)
即得余弦效率的值。通过计算得到年平均余弦效率为0.7730。
因为太阳相对于各面定日镜的位置在实时发生变化,且定日镜面的朝向各异,因此余弦效率的值也在实时进行改变,而余弦效率在较大程度上反映了定日镜场的采光情况,所以对平均光学效率有着较大影响。此外,还可以根据入射光线和反射光线的方向向量求出定日镜面的法线
,利用反射定律 [7] ,有:
(23)
镜面法线的求解对于研究定日镜的朝向以及建立镜面坐标系有着直接联系。
3.2.3. 大气透射率
大气透射率是指大气对可见光的透过程度,用来刻画太阳辐射经过大气层时的折损程度。太阳能传输损失是由大气对太阳辐射的吸收和散射导致的,透射率的范围通常在0到1之间,越接近1表示大气透明度越高,光线传播损失越小。其基本表达式 [8] 为:
(24)
(25)
其中,
为定日镜的安装高度,
为定日镜在镜场的分布半径,
为镜面中心到集热器中心的直线距离。由此可知,当外界环境条件一定时,太阳辐射传输折损与定日镜到集热器中心的直线距离有关,定日镜距集热器越远,太阳辐射传输折损越大。因此,镜场的规模越大,大气透射率越低 [9] 。
3.2.4. 集热器的截断效率
通过分析可得其基本表达式为:
(26)
分析可知,它表示的是集热器截获的能量占镜场汇聚的能量的百分比,而镜场汇聚的总能量应考虑的是去除阴影挡光后能反射出的太阳光线。
由HFLCAL的假设,利用中心极限定理可证明,定日镜在集热器开口平面上的能流密度分布接近圆形高斯分布,且该定理指出两个或多个分布函数(日形误差和镜面误差等)的卷积向高斯分布收敛 [10] 。由此,截断效率
可表示为:
(27)
其中,
为反射光线的整体误差。由Rabl (1985),通量分布的离散度可以计算为:
(28)
其中,
为太阳形状的标准误差,
为定日镜的面形误差,
为定日镜的跟踪误差。而根据Neumann等人(2002)收集的DLR太阳形状的平均剖面和测量结果,可以独立于DNI使用
为2.09 mrad。
此外,本文采用定日镜的平面面形误差 [11]
为1.801 mrad,而根据Sener定日镜的测试结果 [12] ,其跟踪误差
取值为0.63 mrad。
因此,将上述参数代入公式,可得整体误差
为3.033,平均截断效率
为0.9255。
3.2.5. 平均光学效率
基于上文,可推得各月特定时刻的阴影遮挡效率、余弦效率,各定日镜的大气透射率以及集热器的截断效率,不妨取镜面反射率为常数
,可以算出各月特定时刻的光学效率,由此,月平均光学效率和年平均光学效率均可求出,其中各项的月平均值如图5所示:
Figure 5. Line chart for each indicator over month
图5. 各指标随月份变化折线图
从图5可以看出,在各项效率中,月平均光学效率、平均余弦效率、平均阴影遮挡效率、平均截断效率的波动较为稳定,而单位面积镜面平均输出热功率呈现出先增大后减少的趋势。
3.3. 定日镜场的输出热功率
法向直接辐射辐照度DNI (单位:kW/m2)是指地球上垂直于太阳光线的平面单位面积上、单位时间内接收到的太阳辐射能量,可按以下公式近似计算 [7] :
(29)
其中
为圆形定日镜场的海拔高度,计算得
,
,
,由:
(30)
其中
,为太阳常数,
为太阳高度角,根据得到的每一特定时刻的太阳高度角正弦值,可以计算出对应的每一特定时刻的DNI,再根据:
(31)
其中
为第i个定日镜面积,
为其对应的光学效率,可以求出各时刻的输出热功率。
最后,得到年平均输出热功率为38.16 MW,单位面积定日镜面年平均输出热功率为0.6075KW。
3.4. 结果分析
现将年平均光学效率、输出功率等数据结果汇总如表1:
Table 1. Average optical efficiency and output power meter over the year
表1. 年平均的光学效率及输出功率表
从另一角度进行分析,对于每一面定日镜,同样可以求出一年内的平均光学效率,为进一步考察平均值的区域分布特征,本文将这一结果做了可视化处理,具体效果如图6:
Figure 6. Heat map of the annual average optical efficiency of a heliostat under coordinate distribution
图6. 坐标分布下定日镜年平均光学效率热力图
基于以上结果图,可以作出以下分析:该圆形定日镜场在一年中的平均光学效率大致呈现出东北高西南低的分布,这与该地处于北半球相关,太阳光线来自于偏南方位置,因此从镜面朝向来看,应是位于北面的定日镜采光效果更好。上述结果有利于进一步分析决策过程中吸收塔位置坐标的确定,以及定日镜的数量与位置排布等优化配置。
4. 基于遗传算法的定日镜场优化设计
4.1. 模型的构建
针对定日镜场的优化配置问题,本文首先考虑定日镜的排布方式,需要尽量避免定日镜的光线遮挡以获取效率最大化,按照个体效率最大的原则将每个定日镜放置在其效率最优的位置,但这种基于贪婪选择的方式往往计算成本巨大且整体效率较低。因此,基于整体效益最大化原则和光照均匀性理论,本文设立某种具有特定几何特征如圆形径向交错的排列规则,而每两排之间的径向间距相等,这种规则采用较少参数来表征整个镜场的布置,更加切合实际计算要求且效率较高。
本文设计的阴影分布及径向交错排列示意图如图7:
Figure 7. Shadow distribution and radial staggered arrangement
图7. 阴影分布与径向交错排列图
其次,考虑吸收塔的阴影域对定日镜是否有影响,经过计算,吸收塔(塔高80 m)形成的阴影长度超出100 m (即留出空地建造厂房不安装定日镜的范围),因而会对定日镜造成阴影覆盖,本文计算出1~3月、9~12月的阴影角度以及长度进行计算,利用几何知识计算出相应扇形点位并进行去除,对剩余坐标进行遍历计算。
最后,结合上文中考虑的光学效率问题,鉴于平均光学效率大致呈现出东北高西南低的分布,本文将吸收塔的位置约束于西南部分光学效率较低的位置,可以最大化定日镜的光学效率以提高热功率。
基于此,本文建立下列定日镜场的目标优化模型:
决策变量:设立吸收塔位置坐标为
,定日镜的尺寸(其中高为
,宽为
),定日镜的数目为
,定日镜位置
。
约束条件:
a) 定日镜场总的年平均输出热功率达到额定值60 MW
(32)
b) 定日镜安装尺寸限制,如下:
① 镜面边长范围
。
② 安装高度范围
。
③ 保证镜面在绕水平转轴旋转时不会触及地面即
。
④ 相邻间隔
。
a) 定日镜安装位置限制,如下:
考虑到吸热器100 m内不安装定日镜以及镜场的最大径向为350 m,其圆形区域边界限制为:
(33)
b) 吸收塔位置约束域,如下:
目标函数:单位镜面面积平均输出热功率最大化,即
,其中
(34)
综上,目标优化模型为:
(35)
4.2. 模型求解与分析
(1) 算法构建——遗传算法(GA)
针对定日镜场目标优化问题,本文使用遗传算法进行求解。在参数设定时,定日镜的位置、尺寸、数目等决策变量构成了一个高维且存在多个局部最优解的搜索空间,同时受到各种条件约束的限制。而遗传算法通过选择、交叉和变异操作,在整个搜索空间中通过不断的迭代和演化,逐步提升种群中的个体质量,遍历问题中存在的多个局部最优解,最终找到最佳的参数组合。同时,遗传算法具有较高的并行处理能力,能同时对多个个体进行评估,以加快算法的收敛速度和搜索效率。
(2) 定日镜场目标优化模型的求解与分析
对相关指标进行月度分析,发现平均光学效率、平均余弦效率、平均阴影遮挡效率等波动仍较为稳定,而单位面积镜面平均输出热功率仍呈现出先增大后减少的趋势,这是由于每个定日镜尺寸高度等参数保持一致,所以分布规律相似。其中,年平均光学效率、输出功率等数据结果如表2:
Table 2. Optimized annual average optical efficiency and output power meters
表2. 优化后年平均的光学效率及输出功率表
基于遗传算法,利用Python对整个镜场坐标进行遍历搜索,计算出吸收塔的大致位置信息与定日镜的交错散落分布,其中阴影区域是不考虑的定日镜位置范围。每个定日镜的尺寸、安装高度是相同的,因此,定日镜的分布较为均匀、错落有致,且避开了太阳辐射吸收较差的影带,空间分布示意如图8,吸收塔位置坐标与定日镜各个参数见表3。
通过深入分析,可以发现各项指标都在优化布局后基本都得到了提升,此外,通过计算,本文得到了年平均光学效率的优化率为16.65%,年平均输出热功率的优化率为58.20%,单位镜面面积的年平均输出热功率的优化率为12.18%,有力地说明了本文对定日镜场的优化设计是有效且值得应用推广的。
Table 3. Absorption tower location coordinates and heliostat parameters
表3. 吸收塔位置坐标与定日镜各个参数表
Figure 8. Schematic diagram of the location distribution of the absorption tower and heliostat
图8. 吸收塔与定日镜的位置分布示意图
由于本文所研究的定日镜场地处北半球,东北方位的定日镜朝向偏南,能够更好地有效采光,尽可能地减少了采光时造成的损失,所以该区域的排布更集中,且定日镜的安装高度为2.75 m,这可能和阴影遮挡效率、大气透射率有关,因为设置较高高度可能会造成较大的阴影域,遮挡部分的定日镜,而设置过低的高度与吸收塔的距离较远,太阳辐射传输折损较多,降低大气透射率影响太阳辐射量,从这个角度来看,各个数据合理且符合预期效果。
5. 结论
本文应用了数学物理相关理论,建立了定日镜场的光热转化模型,并通过遗传算法对定日镜场进行优化配置设计,使得达到额定功率60 MW的情况下,单位镜面面积年平均输出热功率尽可能大,得到了年平均光学效率为0.6929,较优化前的0.5940提升了16.65%,单位镜面输出热功率为0.6815 kW/m2,较优化前的0.6075 kW/m2提升了12.18%。同时,本文所研究的模型可以推广到太阳能、建筑设计、温室农业等领域,对定日镜进行调节提高能源利用效率和可再生能源的使用比例,还可以与其他能源系统相结合,如光伏发电、风能等,形成综合的能源利用系统。
NOTES
*通讯作者。