备战2024全国大学生数学建模竞赛A题第一问思路-定日镜场的优化设计模型

摘 要

定日镜场能将太阳光反射汇聚到安装在镜场中的吸收塔集热器上，实现能量转换不同的定日镜场参数会影响定日镜场的效率。本文通过光线反射、入射等模型，计算给定参数下的定日镜场效率，并建立定日镜场的效率优化模型，研究定日镜场参数变化对单位镜面面积平均输出热功率的影响。

针对问题一，本文建立了定日镜场的效率计算模型。首先，建立镜场大地坐标系镜面坐标系和光锥束坐标系，描述定日镜场的参数。太阳光束的入射光线和反射光线考虑为光锥束，计算太阳光的入射和反射，考虑阴影遮挡损失、余弦损失、集热器截断效率等因素。通过坐标变换和方程联立，判断某入射光线及对应的反射光线是否造成阴影遮挡损失所涉及的塔身阴影损失、入射镜面遮挡损失、反射镜面遮挡损失，再判断成功反射出的光线是否被集热器接收。最后得到光学效率和输出热功率的表达式进一步计算年平均光学效率、年平均热功率、单位镜面面积年平均输出热功率等。在某个定日镜上取点步长t为1m，在光锥束上取6=0.002rad，6=0,m/2,m,3x/2，进行网格化构建入射光线束和反射光线束，计算得到表1和表2。该定日镜场的年平均输出光学效率0.536230167、年平均输出热功率32.76117051MW、单位镜面面积年平均输出热功率 0.521508604kW/m?，最后，进行敏感性分析，分析步长，对结果的影响。

针对问题二，本文提出了单位镜面面积年平均输出热功率的优化模型。问题二是在定日镜场的额定年平均输出热功率为60MW的条件下，设计定日镜场的参数，使得单位镜面面积年平均输出热功率尽量大。决策变量包括吸收塔的位置坐标、定日镜的尺寸(相同)、安装高度、数量和位置。目标函数是单位镜面面积年平均输出热功率的最大值，约束条件包括镜面边长在2m至8m之间、安装高度在2m至 6m 之间、相邻定日镜底座中心距离比镜面宽度多5m等。先取w,v相等，用蜂窝排列法确定圆环内能排列的定日镜数目和位置。然后将蜂窝进行绕原点旋转，用蒙特卡洛模拟法对部分定日镜进行随机抽样，得到粗决策变量，从而进一步计算光学效率和单位面积输出热功率等。最后遍历“w=v,h,X,”寻找决策变量的最优解。得到该定日镜场的年平均输出光学效率 0.591643667、年平均输出热功率68.24427914MW、单位镜面面积年平均输出热功率0.572538333kW/m2。

针对问题三，本文提出了单位镜面面积年平均输出热功率的优化模型。问题三是在定日镜尺寸和安装高度可以不同的情况下，设计定日镜场的参数，使得单位镜面面积年平均输出热功率尽量大。以单位镜面面积年平均输出热功率为目标函数，以吸收塔位置坐标、定日镜尺寸、安装高度、定日镜数目、定日镜位置为决策变量，以额定功率、相邻定日镜距离、圆形区域半径等为约束条件，建立优化模型。考虑通过第二问求得结果进一步优化，将最外围的定日镜安装高度升高为6m，其他不变，得到该定日镜场的年平均输出光学效率0.496428083、年平均输出热功率60.336111MW、单位镜面面积年平均输出热功率0.506192417kWW/m2。

关键词:坐标旋转 光学效率 输出热功率 锥形光束 蜂窝排列法

一、问题分析

1.1 问题一的分析

问题一是在定日镜场参数确定时，计算定日镜场的效率。同时，题中简化计算“年均”指标只需计算每月21日5个时间点等。为了得到最大的效率，我们认为经过定日镜中心反射后的反射光线通过集热器中心，那么其他入射光线也能大部分被集热器接受。因此，太阳位置确定后，经过定日镜中心的入射光线和反射光线确定，可根据向量法则计算定日镜的法向量，从而确定每个定日镜场的俯仰角和方位角。阴影遮挡损伤考虑塔身的阴影损伤、后排镜面接受的太阳光线被前镜阻挡的阴影损伤和后排镜面反射的太阳光线被前镜阻挡的遮挡损伤。余弦损伤考虑为入射光线和镜面法线不平行的能量损伤。最后，建立模型表示集热器的接收能量，计算定日镜场的相应量。

1.2 问题二的分析

优化模型下，吸收塔的位置坐标、定日镜尺寸、安装高度、定日镜数量、定日镜位置为决策变量，单位面积年平均输出热功率尽量大为目标函数，每个定日镜尺寸和安装高度相同，在满足约束条件的情况下，使得目标的数尽量大。进行逐步寻优，先确定定日镜尺寸和安装高度，之后逐步寻找其他决策变量的优解，使得目标函数达到局部最优解即可。

1.3 问题三的分析

每个定日镜尺寸和安装高度不同，确定定日镜场的其他参数，在满足约束条件的情况下，使得单位面积年平均输出热功率尽量大。根据问题二的寻优过程，适当调整尺寸和高度。

二、模型假设

为了对模型进行合理简化，我们建立了以下的模型假设:
1、假设太阳光直接照射到集热器上，不计入我们模型的计算。
2、假设相对于3000m的海拔，定日镜作为垂直于太阳光线的平面，计算DNI时的海拔可以忽略定日镜的高度。
3、假设吸收塔的直径和集热器的直径相同，为7m。
4、假设安装定日镜场时不考虑人工成本、土地利用率等。

 注:其他符号将在下文给出具体说明

四、问题一的模型建立与求解

4.1 问题一的模型准备

4.1.1 镜场大地坐标系的建立

在计算定日镜场的效率时，由于一年中太阳是运动的，因此需要确定太阳的位置才能计算定日镜场的效率。根据题目描述，为了简化计算，“平均”指标只需在每月21日的9:00、10:30、12:00、13:30、15:00计算，因此当太阳位置位于每月21日的5个时间点时，太阳位置固定，可计算镜场效率。首先，以形区域中心为原点，正东方向为x轴方向，正北方向为y轴方向，垂直地面方向为z轴方向建立定日镜场的大地坐标系。

 在该坐标系中，第i个定日镜中心位置坐标d(x,”,z)，集热器中心坐标J(0,0,H)

4.1.2 太阳主光线入射、反射以及定日镜法向

太阳光并非平行光线，而是具有一定锥形角的一束锥形光线。因此，先对太阳的一条主光线进行描述。进行在计算定日镜场效率时，需要首先需要对定日镜的入射主光线、反射主光线、法向、俯仰角和方位角进行建模。首先，根据太阳高度角和太阳方位角的定义和公式，结合下图，则可得到太阳光的入射主光线方向向量如下:

 其中，根据太阳方位角和太阳高度角的公式，ST为当地时间，8为太阳赤纬角，D为以春分作为第0天起算的天数，@为太阳时角，p为当地纬度:

 而正常来说，对于中国区域，早上太阳光从东边射来，中午太阳光从南边射来，傍晚太阳光从西边射来。早上的太阳方位角在90”左右(但一年当中，有一定的角度范围变化)，正中午的太阳方位角在180°(正南方)，傍晚的太阳方位角在270°左右(但一年当中，有一定的角度范围变化)。因此，对于太阳方位角的正弦值需要进行分类讨论，如下:

我们认为，控制器控制定日镜的法向(即控制定日镜的俯仰角和方位角)使得经过定日镜中心的反射主光线刚好经过集热器中心，即为定日镜法向对不同太阳位置的最佳调节。因此，经过定日镜中心的反射主光线单位向量为:

 根据平行四边形定理，从而得到该时刻太阳位置下，第i个定日镜的法向量为(以下简称定日镜法向量):

 通过镜场大地坐标系z轴的方向向量Z=(0,0,1)，计算定日镜的俯仰角B，即定日镜的法向量与镜场大地坐标系投影轴的夹角，如下图:

图3:定日镜的方位角(a)和仰角()计算

 将定日镜的法向量投影于大地坐标系的xoy平面，得到定日镜法向量的投影向量元，计算元与大地坐标系x轴的方向向量S的夹角为方位角 afa，如下:

其中s=(1,0,0),n=(n(1),n(2),0)。光学效率包含了阴影遮挡效率、余弦效率、大气透射率、集热器截断效率和镜面反射率。因此确定太阳的反射和入射及定日镜的法向调节后，需表示太阳光线入射反射造成的阴影遮挡效率、集热器接收情况和余弦效率

4.1.3 定日镜面坐标系的建立
为了方便表示和计算定日镜面上的坐标，我们需要将定日镜面的坐标系转换到镜场大地坐标系或者将镜场大地坐标系转换到定日镜面坐标系。因此，以定日镜中心为原点，以经过原点d并垂直镜面的法向量为z轴。以经过原点平行镜面的宽为x轴，以经过原点平行镜面的长为y轴，建立镜面坐标系(见下图):

 4.1.4光锥坐标系的建立
太阳光并非平行光线，而是具有一定锥形角的一束锥形光线，因此太阳入射光线经定日镜任意一点的反射光线也是一束锥形光线。在考虑阴影遮挡和集热器吸收能量时，需要考虑锥形光束的吸收和遮挡。因此，我们需要建立起以主入射光线和主反射光线为锥体中心线的光锥坐标系。大地坐标系下的反射光锥束的中心主光线L.=(m,",p)，以主光线为z轴，则z轴的方向向量为了=(0,0,1)，则锥体光线束的俯仰角B为:

 

 

 

 

 

 

 

 

通过距离d判断是否相邻，考虑题中要求相邻定日镜底座中心之间的距离比镜面宽度邻接矩阵如下;多5m。先计算距离，得到d'矩阵，观察距离大小，给定d=20， 

 其中，I(i,j)=1表示第i块定日镜和第j块定日镜相邻，否则不相邻。根据此减少a定日镜的阴影遮挡损失的定日面遍历数，更加精确地定位问题定日镜。

4.3.2 网格化镜面和取点求解

在某个定日镜上取点，步长1.5m，即镜上取25个点，在光锥束上取6=0.002rad，0,=0,π/2,π,3x/2，即光锥束取4个点。首先，计算以3月21日为起止的天数D:

 根据不同月份、时间内，计算所对应太阳位置。而后需要计算定日镜的仰角和方位角。而后，需要遍历每一个入射光线圆锥的光线，根据入射主光线，计算主光线锥体系的旋转矩阵，再通过旋转矩阵将不同坐标系的坐标进行转换。
接着判断入射光线(锥束)是否被塔身遮挡，若被塔身遮挡则造成阴影遮挡损失然后，判断入射光线(锥束)是否被其他光镜遮挡。根据几何关系，计算所对应的反射光线(锥束)，判断是否被其他镜遮挡。若未造成遮挡，则计算集热器是否接受。之后计算定日镜场的计算。

4.3.3 结果计算
根据定日镜场效率计算的公式，求解5个时刻的平均值，可得到表1:表1:问题1的每月21日平均光学效率及输出功率

 

 

 

参考资料：

2023全国大学生数学建模竞赛A题思路-定日镜场的优化设计模型 | 2024国赛数学建模_CSDN

完整版本可以通过下方名片加群获取：