基于太阳能光热发电系统效率问题的方案分析——以2023年全国数学建模竞赛A题为例.pdf

Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(4),1778-1799 Published Online April 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.134168 文章引用文章引用:王秉荃,徐洁,邢建民.基于太阳能光热发电系统效率问题的方案分析J.应用数学进展,2024,13(4):1778-1799.DOI:10.12677/aam.2024.134168 基于太阳能光热发电系统效率问题的基于太阳能光热发电系统效率问题的 方案分析方案分析 以以2023年全国数学建模竞赛年全国数学建模竞赛A题为例题为例 王王秉荃秉荃1，徐徐 洁洁1，邢建民邢建民2 1青岛科技大学化学与分子工程学院，山东 青岛 2青岛科技大学数理学院，山东 青岛 收稿日期：2024年3月28日；录用日期：2024年4月23日；发布日期：2024年4月30日 摘摘 要要 本文本文以以2023年高教社杯全国大学生数学建模竞赛年高教社杯全国大学生数学建模竞赛A题为例，分析了题为例，分析了“定日镜场的优化设计定日镜场的优化设计”的可行方案，的可行方案，利用卷积积分和构建模拟镜场等方法，在定日镜工作原理的基础上，建立了优化模型，并通过遗传算法利用卷积积分和构建模拟镜场等方法，在定日镜工作原理的基础上，建立了优化模型，并通过遗传算法和多为粒子群算法对模型进行求解，以解决定日镜场优化设计中的和多为粒子群算法对模型进行求解，以解决定日镜场优化设计中的关键关键问题。本文提供的计算与优化方问题。本文提供的计算与优化方案具有良好的可借鉴性，可为同类问题的分析解决提供参考。由于模型未考虑天气因素、镜面弯曲程度、案具有良好的可借鉴性，可为同类问题的分析解决提供参考。由于模型未考虑天气因素、镜面弯曲程度、灰尘覆盖等复杂因素的影响，对所求结果有一定的影响，有待后续分析探讨。灰尘覆盖等复杂因素的影响，对所求结果有一定的影响，有待后续分析探讨。关键词关键词 模拟镜场，效率云图，卷积积分，模拟镜场，效率云图，卷积积分，HFLCAL模型，多维粒子群算法模型，多维粒子群算法 Efficiency Analysis of the Solar Energy System for Electric Generation via Light-to-Heat Conversion A Case Study of Problem A in the 2023 National Mathematical Modeling Competition Bingquan Wang1,Jie Xu1,Jianmin Xing2 1School of Chemistry&Molecular Engineering,Qingdao University of Science and Technology,Qingdao Shandong 2School of Mathematics&Physics,Qingdao University of Science and Technology,Qingdao Shandong 王秉荃 等 DOI:10.12677/aam.2024.134168 1779 应用数学进展 Received:Mar.28th,2024;accepted:Apr.23rd,2024;published:Apr.30th,2024 Abstract This paper presents a feasible resolution to the optimization design of helioscope field,based on the problem A of 2023 Higher Education Press Cup National Undergraduates Mathematical Mod-eling Competition.Following the working principle of helioscope,an optimization model is estab-lished by means of convolutional integral and construction of simulated optical field,and the model is solved by genetic algorithm and multi-particle swarm optimization algorithm to solve the problems in the optimization design of helioscope field.The calculation and optimization scheme proposed provides a good solution for similar problem.Lacking the effects of some practical pa-rameters,i.e.weather conditions,mirror curvature and covering of dust,which may somehow af-fect the accuracy of the results,further analysis and improvement are still desirable.Keywords Simulated Helioscope Field,Efficiency Cloud Map,Convolution Integral,HFLCAL Model,Particle Swarm Optimization Copyright 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 面对日益紧缺的能源供应与严重的能源问题，寻求环保持久的新能源代替化石燃料越来越受到重视。近年来，新能源技术不断发展，新能源发电系统应用广泛。在各种新兴能源中，太阳能的来源最稳定持续和丰富。塔式太阳能光热发电是一种低碳环保的新型清洁能源技术1，其发电系统的效率优于传统的光热发电效率，对实现我国“碳达峰”和“碳中和”的战略目标具有重要意义1。构建以新能源为主体的新型电力系统，是我国实现“碳达峰”“碳中和”目标的一项重要措施。定日镜是塔式太阳能光热发电站(简称塔式电站)收集太阳能的基本组件，其底座由纵向转轴和水平转轴组成，平面反射镜安装在水平转轴上。纵向转轴的轴线与地面垂直，可以控制反射镜的方位角。水平转轴的轴线与地面平行，可以控制反射镜的俯仰角。两转轴的交点(定日镜中心)离地面的高度称为定日镜的安装高度。塔式电站利用大量的定日镜组成阵列，称为定日镜场。定日镜将太阳光反射汇聚到安装在镜场中吸收塔顶端上的集热器，加热其中的导热介质，并将太阳能以热能形式储存起来，再经过热交换实现由热能向电能的转化。太阳光并非平行光线，而是具有一定锥形角的一束锥形光线，因此太阳入射光线经定日镜任意一点的反射光线也是一束锥形光线2。定日镜在工作时，控制系统根据太阳的位置实时控制定日镜的法向，使得太阳中心点发出的光线经定日镜中心反射后指向集热器中心。集热器中心的离地高度称为吸收塔高度。塔式太阳能光热发电系统的效率与定日镜的太阳光追踪技术和定日镜的位置坐标息息相关，定日镜场的优化设计问题至关重要。2.问题简析问题简析 本文基于定日镜的工作原理，通过卷积积分和构建模拟镜场等方法，建立定日镜的工作优化模型，Open AccessOpen Access王秉荃 等 DOI:10.12677/aam.2024.134168 1780 应用数学进展 同时利用多维粒子群算法和遗传算法对建立模型进行求解。问题问题 1 将吸收塔建于圆形定日镜场中心，定日镜尺寸均为 6 m 6 m，安装高度均为 4 m，且给定所有定日镜中心的位置(以下简称为定日镜位置)，计算该定日镜场的年平均光学效率、年平均输出热功率，以及单位镜面面积年平均输出热功率。解决方案：建立光学效率模型，其模型中包含阴影遮挡效率、余弦效率和集热器截断效率等多个影响因素，利用平面投影法计算阴影遮挡效率，然后通过 MATLAB 编程画出效率云图，再通过云图分析镜场光学效率的分布规律，然后利用模型计算每个定日镜在每一特殊时刻的光学效率，再对所有时间点的数据取平均值，得到年平均光学效率；再建立输出热功率模型。利用定日镜面积、光学效率与 DNI 计算出总输出热功率，对所有时间点的数据进行平均，求得年平均输出热功率；最后利用总输出热功率除以总镜面面积，对所有时间点数据进行平均，得到单位镜面面积年平均输出热功率。问题问题 2 按设计要求，定日镜场的额定年平均输出热功率(简称额定功率)为 60 MW。若所有定日镜尺寸及安装高度相同，设计定日镜场的以下参数：吸收塔的位置坐标、定日镜尺寸、安装高度、定日镜数目、定日镜位置，使得定日镜场在达到额定功率的条件下，单位镜面面积年平均输出热功率尽量大。解决方案：在第一问的基础上，先建立 CAMPO 布置仿真镜场，再用粒子群算法优化模型，目标函数为单位镜面面积年平均输出热功率最大，约束条件为吸收塔周围 100 m 内不安装定日镜，安装高度使镜面绕水平转轴旋转时不触及地面，同时相邻定日镜底座中心间的距离比镜面宽度多 5 m 以上。然后利用问题 1 中建立的模型与优化模型进行优化求解，算得定日镜尺寸为 4.7064 4.5954，定日镜安装高度为 4.2006 m，优化后定日镜数目为 5128 个，优化后定日镜总面积 110907.31 m2。问题问题 3 如果定日镜尺寸可以不同，安装高度也可以不同，额定功率设置同问题 2，重新设计定日镜场的各个参数，使得定日镜场在达到额定功率的条件下，单位镜面面积年平均输出热功率尽量大。解决方案：沿用问题 1 和 2 建立优化模型，将定日镜尺寸与安装高度变为可变参数，在定日镜场达到额定功率 60 MW 条件下，应用遗传算法使其单位镜面面积年平均输出热功率尽量大，然后对模型进行优化求解，算得优化后定日镜总片数 4610 片，单位面积年平均输出功率 0.5318 m2，定日镜总面积 158,620 m2。3.符号说明与模型假设符号说明与模型假设 本文符号及说明如表 1 所示。Table 1.Symbols and illustrations 表表 1.符号及说明 序号 符号 说明 1 s 太阳高度角 2 s 太阳方位角 3 当地纬度 4 太阳时角 5 ST 当地时间 6 太阳赤纬角 7 DNI 法向直接辐射辐照量 8 G0 太阳常数 9 N 定日镜总数 王秉荃 等 DOI:10.12677/aam.2024.134168 1781 应用数学进展 续表 10 sb 阴影遮挡效率 11 at 大气透射率 12 trunc 集热器截断效率 13 H0 吸收塔高度 14 i 入射向量 15 反射向量 为了对模型进行合理的简化，本文做出如下假设：1)假设镜面反射率 ref取为 0.92；2)镜面为理想平面；3)不考虑阴雨天气对太阳辐射的影响；4)假设本题计算中的春分日为 3 月 21 日；5)假设阴影与遮挡只发生在目标定日镜周围的一定范围内。4.问题问题 1 模型的建立与求解模型的建立与求解 4.1.问题分析问题分析 问题 1 要求通过附件所给公式，建立模型计算定日镜安装高度为 4m，定日镜尺寸大小为 6m6m,且已知定日镜中心位置时定日镜场的年平均光学效率、年平均输出热功率及单位镜面面积年平均输出热功率。首先，对太阳高度角和太阳方位角等后续模型所需数据进行计算，其次，建立光学效率模型，分析阴影遮挡效率、余弦效率和集热器截断效率等多个影响因素，然后通过 MATLAB 编程画出效率云图，通过云图分析镜场光学效率的分布规律，利用模型计算每个定日镜在每一特殊时刻的光学效率，对所有时间点的数据取平均值，得到年平均光学效率；再建立输出热功率模型，通过定日镜面积、光学效率与 DNI计算出总输出热功率，对所有时间点的数据进行平均，求得年平均输出热功率；最后利用总输出热功率除以总镜面面积，对所有时间点数据进行平均，得到单位镜面面积年平均输出热功率3。4.2.模型的建立模型的建立 4.2.1.基础数据计算基础数据计算 首先，对后续建立模型经常用到的公式进行处理。太阳角高度 s sincos coscossinsins=+(1)太阳方位角 s sinsinsincoscoscossss=(2)式中，为当地纬度，北纬为正；为太阳时角()1212ST=(3)其中，ST代表当地时间，指太阳赤纬角 王秉荃 等 DOI:10.12677/aam.2024.134168 1782 应用数学进展 22sinsinsin23.45365365D=(4)D 为以春分日作为第 0 天起算的天数。因而可以通过日期和当地时间来推导出太阳时角和太阳赤纬角，再通过推导出的数据利用公式(1)和(2)计算出后续计算所需要的太阳高度角与太阳赤纬角。4.2.2.光学效率模型光学效率模型 定日镜的光学效率 为：cossbattruncref =(5)其中，镜面反射率 ref可取为常数，本题取为 0.92。具体计算如下：1)阴影遮挡效率 sb Figure 1.Illustration of shadow formation by sheltering 图图 1.阴影遮挡示意图 阴影遮挡如图 1 所示。阴影定日镜在入射坐标系中的坐标：()()sincossincossincoscossincossinisossosissossossosissossossosXXYYSYZZYXZ=+=+=+(6)阴影定日镜在入射坐标系中的坐标：cosuiuimuiXXYYZZZ=(7)Zm：阴影定日镜中心在被阴影定日镜面坐标 Z 坐标，由题知 Zm=0。由地面坐标系到镜面坐标系中变换矩阵：sincos00coscoscossinsin0sincossinsincos00001HHHHHHHHHHHHAAEAEAEMEAEAE=(8)王秉荃 等 DOI:10.12677/aam.2024.134168 1783 应用数学进展 由入射坐标系到镜面坐标系中变换矩阵：T432MM M=(9)阴影定日镜的镜面中心沿入射光方向在被阴影镜面上投影坐标：()()()()coscossincossincossincoscoscoscossinsinmuHsmuHHsuHsHsHsuHsHsHsXXAYXEAYEAEZEAE=+(10)根据下式判断是否存在阴影 0mwimlXHZYH且 (11)其中，Hw为镜面宽度，Hl为镜面高度。以镜面中心为坐标系原点，其法向量为 z 轴，镜面两条边建立镜面坐标系，利用矩阵 M，可得镜面追迹点在地面系中坐标。()()()0000sincoscoscossincossin4mgHHHmgHHHmgHmxAxAEyxyAxAEzEy=+=+(12)计算相邻定日镜面之间阴影率，以被遮挡的定日镜中心为原点，建立阴影坐标系，X0指向水平正南，Y0指向正东，Z0与 X0OY0面垂直，阴影定日镜在 X0Y0Z0坐标系中的坐标为(x0，y0，4)，(Xsg，Ysg，Zsg)，(X0g，Y0g，Z0g)表示定日镜及被定日镜在地面坐标系中的坐标。000000ssggssggssggXXXYYYZZZ=(13)以定日镜为中心，入射光线为 Z 轴建立入射坐标系，由地面坐标系到入射坐标系的变换矩阵为 2sincos00sincossinsincos0coscoscossinsin00001ssssssssssssM=(14)()()001tan10038.72138.7sssbsHHH+=集 (15)其中，H0=80 m，为吸收塔的高度，H集=8 m 为集热器高度，s为太阳高度角。通过参考资料及后续模拟仿真，当太阳高度角达到 38.7 度左右时，镜场中的阴影遮挡现象完全消失4。实际计算过程如下：STEP1：将定日镜均匀划分成 m n 网格(图 2)，计算网格上各点映射到定日镜场坐标系中的坐标。STEP2：循环遍历每一块定日镜，对它四周一圈的定日镜进行重合网格点判断并计算重合面积。STEP3：将所有重合的面积和除以总定日镜面积即得阴影遮挡效率。sb=阴影遮挡积镜场总积面面 王秉荃 等 DOI:10.12677/aam.2024.134168 1784 应用数学进展 Figure 2.Illustration of grid-divided helioscopefield 图图 2.定日镜网格表示示意图 2)集热器截断效率 trunc 集热器截断效率表示在太阳能光热系统中，定日镜将太阳光线聚焦到集热器上时光线因为镜子的截断而失去的概率。这个效率取决于镜面的形状和位置，以及集热器的大小和位置。参考文献5采用HFLCAL 模型计算集热器截断效率效率。22221expd d2truncxytottotxyy x+=(16)()22222totsunbqasttrackd=+(17)()222bqs=(18)()220.54ttastHWd+=(19)costlwdHHHwf=(20)cos1slwdWHHwf=(21)上式中，2tot指吸热器表面光斑的总标准差，其受反射距离 d、定日镜焦距 f，太阳形状误差标准差 sun、光束质量误差标准差 bp、像散误差标准差 ast、及跟踪误差标准差 track。s指的是斜率误差标准差；Ht、Ws为吸热器装置上光斑在子午方向与弧矢方向的尺寸大小；Hl、Hw分别为定日镜的高度和宽度；f 是定日镜的轴向焦距；w 为太阳光线入射角。根据搜集数据，sun=2.51 mrad、s=0.94 mrad、track=0.63 mrad。3)余弦效率 cos(太阳入射角余弦值 cos)大部分光学效率损失为镜场的余弦损失，入射光线与镜面反射点处的法线方向存在夹角 (图 3)，其余弦值 cos 即为所求的余弦效率值。利用题中所建立的坐标系，以正东为 x 轴正方向，正北为 y 轴正方向建立镜场坐标系，集热器坐标为(0，0，H0)，镜面 A 坐标为(x，y，h0)。定日镜镜面中心指向集热器的光线的单位法向量n()()0022200,xy HhxyHh=+n (22)入射向量i()coscos,cossin,sinsssss=i (23)反射向量 (),80 xy=(24)王秉荃 等 DOI:10.12677/aam.2024.134168 1785 应用数学进展 Figure 3.Illustration of incident angle 图图 3.入射角示意图 单一定日镜的余弦效率可通过入射光线单位向量与镜面法向量的点积计算：coscos=i n (25)大气透射率 at()820.99321 0.00011761.97 101000atHRHRHRddd=+(26)dHR表示镜面中心到集热器中心的距离，故 22280HRdxy=+(27)带入公式求出大气透射率 at。最终建立模型为()()()002sincoscoscossinsinsinsinsincoscoscos121222sinsinsin23.453653651tan10038.72s.t.138.71exp2sssssbcosattruncrefsssbstruncytotSTDHHH =+=+=集()()()()2220022200822220.993210.00011761.97 10d dcos,coscos,cossin,si0080n10 xtotcossssssatHRHRHRRHdddxyy xxy HhxyHhdxy+=+=+=+i nni (28)王秉荃 等 DOI:10.12677/aam.2024.134168 1786 应用数学进展 4.2.3.平均输出热功率模型平均输出热功率模型 定日镜场的输出热功率 Efield NfieldiiiEDNIA=(29)法向直接辐射辐照量 DNI 为地球中垂直太阳光线的平面单位面积、单位时间内所受到的太阳辐射能量，可用如下公式计算：()()()0222expsin0.42370.00821 60.50550.00595 6.50.2711 0.01858 2.5scDNIGabaHbHcH=+=+=+(30)G0为太阳常数，其值取为 1.366 kW/m2，H 为海拔高度，本题 H=3 km。其中 DNI 为法向直接辐射辐照量，N 为定日镜总数，Ai为第 i 面定日镜采光面积；i为第 i 面镜子的光学效率。最终建立模型为()()()()0222expsin0.42370.00821 60.50550.00595 6.50.27110.01858 2.5sincoscoscossinsinsinsinsincoscoscos1212sinss.t.NfieldiiisssssEDNIAcDNIGabaHbHcHST=+=+=+=+=()()()()002222002220022insin23.453653651tan10038.72138.71expd d2cos,coscos,sbcosattruncrefsssbstruncxytottotcosssDHHHxyy xxy HhxyHh =+=+=+=i nni集()()22282cossin,sin80.99321 0.00011761.97 1001000ssatHRHHsRRHRddddxy=+=+(31)王秉荃 等 DOI:10.12677/aam.2024.134168 1787 应用数学进展 4.2.4.模拟镜场的建立及光学效率分布规律模拟镜场的建立及光学效率分布规律 为了更好表达太阳位置变化对定日镜场中各部分光学效率影响规律，参考文献5，利用 MATLAB 建立针对本题的模拟镜场。其中每一个黑色空心圆点代表一台定日镜，而在圆心处的红色圆点代指吸热塔。镜场半径为 350 m，第一行定日镜距离吸收塔的距离为 100 m。横轴为镜场的东西方向且正半轴为正东，为方便将不同日期和不同时间点的仿真结果做对比，仿真结果的效率云图使用相同的颜色表。其中，红色为高的效率值，蓝色为低的效率值，颜色表的取值范围是 01。通过建立模拟镜场来研究光学效率的分布规律，在特殊日期的时刻对目标镜场进行效率仿真，记录仿真后各个定日镜的效率值，再绘制镜场在各时间点的效率分布云图和日均效率分布云图。分别计算镜场效率的最高值、均值及最差值，并用百分比来表示。通过不同的颜色来区分定日镜的效率值，绘制镜场效率云图的方式展示镜场光学效率的分布规律。4.3.模型的模型的求解求解 4.3.1.效率云图的求解效率云图的求解 从春分、夏至、秋分和冬至的余弦效率分布云图(图 4)可以看出，其在镜场内分布并不均，且镜场中红色区面积和蓝色区面积的占比随着时间而变化。同一天，太阳高度角越大，镜场内红色面积越大，镜场余弦效率越高；总的来看，镜场北半边比南半边的余弦效率高。Figure 4.Cosine efficiency comparisons 图图 4.余弦效率对比图 王秉荃 等 DOI:10.12677/aam.2024.134168 1788 应用数学进展 从图 5 可看出，在圆形镜场中，同一环上的定日镜大气透效率值是相等的。从整体来看，越靠近吸收塔的定日镜的大气透射率越高，而远离的大气透射率越低。边缘定日镜相比中心定日镜的大气衰减效率值会存在 2.2%左右的下降。Figure 5.Shadow-affected efficiency comparisons 图图 5.阴影遮挡效率对比图 王秉荃 等 DOI:10.12677/aam.2024.134168 1789 应用数学进展 Figure 6.Simulation diagrams of atmospheric transmissivity 图图 6.镜场大气透射率仿真图 对比不同四个月的综合效率分布云图(图6)可以发现，镜场效率分布不均其主要是余弦效率的影响，其次为大气透射率。并且南半边镜场的光学效率明显低于北半边镜场。4.3.2.效率云图的求解效率云图的求解 利用建立的模型，通过 MATLAB 进行编程求解，求解出 X 年每月 21 日平均光学效率及输出功率；然后为简化计算，以当地时间每月 21 日 9:00、10:30、12:00、13:30 与 15:00 五个时间点计算出的结果取平均值作为年均指标，所求结果如表 2 与表 3。Table 2.Average optical efficiency and output power of 21st every month in question 1 表表 2.问题 1 每月 21 日平均光学效率及输出功率表 日期 平均光学效率 平均余弦效率 平均阴影遮挡效率 平均截断效率 单位面积镜面 平均输出热功率(kW/m2)1 月 21 日 0.513637083 0.888268603 0.905636502 0.713342735 0.495919819 2 月 21 日 0.517224871 0.888113154 0.911048745 0.715045616 0.505121568 3 月 21 日 0.518505967 0.88806065 0.913369467 0.715795539 0.508618074 王秉荃 等 DOI:10.12677/aam.2024.134168 1790 应用数学进展 续表 4 月 21 日 0.517320867 0.888125011 0.910577331 0.714894578 0.504744835 5 月 21 日 0.513796118 0.888278965 0.905315632 0.713243665 0.495711224 6 月 21 日 0.512020177 0.888373968 0.902702206 0.712452786 0.490949363 7 月 21 日 0.513850017 0.888276358 0.905394701 0.713268016 0.495852741 8 月 21 日 0.517471388 0.888119071 0.91081672 0.7149713 0.505131349 9 月 21 日 0.518446213 0.888060898 0.913362236 0.715793278 0.508552655 10 月 21 日 0.516878708 0.888127066 0.910496011 0.71486863 0.504231793 11 月 21 日 0.513275436 0.88828718 0.905083344 0.713172827 0.494945178 12 月 21 日 0.511726233 0.888373766 0.902705331 0.712453462 0.490671277 Table 3.Average annual optical efficiency and output power in question 1 表表 3.问题 1 年平均光学效率及输出功率表 年平均 光学效率 年平均 余弦效率 年平均 阴影遮挡效率 年平均 截断效率 年平均 输出热功率 单位面积镜面 年平均输出热功率 0.515346 0.88820539 0.90804235 0.71410853 31.4111715 0.50001864 5.问题问题 2 模型的建立与求解模型的建立与求解 5.1.问题分析问题分析 问题二是在定日镜场额定功率为 60 MV、所有定日镜的尺寸大小和安装高度相同的条件下，对定日镜场进行设计优化。因此在第一问的结果基础上建立优化模型。首先用 Campo 布置方法，建立仿真定日镜场6，再用粒子群算法进行优化，先得到定日镜场年工作效率的最大值时，定日镜场里每个镜子的坐标位置，再进行二次优化，将功率限制为 60 MV，得到新功率下的定日镜场镜子数目和每个镜子的坐标位置。5.2.模型的建立模型的建立 5.2.1.Campo 布置方法模型布置方法模型 问题二在问题一模型的基础上，建立了 Campo 布置方法模型：112mdNr=(32)其中，r1为第一行定日镜的半径，dm为定日镜特征圆的直径，N1为第一行定日镜的数目。112arcsin2mdr=(33)其中，1为第一行定日镜之间的方位角。mhddd=+安 (34)其中，dh为定日镜对角线长度，d安为定日镜之间的安全距离。因为每行定日镜的数目是相同的，所以连续行之间的径向间隔也固定不变，且奇数行和偶数行连续交替分布。通过这种最密集的排列方式，得到相邻的行半径的最小变化量rmin。王秉荃 等 DOI:10.12677/aam.2024.134168 1791 应用数学进展 mincos30mrdh=(35)22114mdhrr=(36)因为每行上的定日镜数目是相同的，所以同一行上定日镜之间的间距会随着行数的增加而逐渐变大。规定最后一行上两个相邻的定日镜之间只能放置一个额外的定日镜，然后才能进行下个部分的设置，接下来的这些部分与第一部分的排列原理是相同的。下个部分上的每一行定日镜的数目是上个部分每一行定日镜数目的两倍，所以，在每一部分各行的定日镜数目 Nn和每一部分第一行的行半径1nr能够下面的两个式子计算得出。其中，n 为被设置部分的序号。112nnNN=(37)112nnrr=(38)因为在相同的部分内相邻两行之间的径向间隔固定不变，所以能够通过这种特殊的设置布局计算第n 个部分内能够接受的最大行数 Nmax。11maxmin2nrNr=(39)根据光热发电站的设计需要，先选择合适的第一行定日镜的数目，然后再通过运算预测定日镜场中每个部分的定日镜行数、每一行上的定日镜数目以及每一行上的行半径，得到由 Campo 布置方法模型合成的紧密型定日镜场。再通过软件进行仿真，能够看到紧密型定日镜场具有很低的阴影遮挡效率和很高的其他光学效率。在进行本问的定日镜场的优化布置之前，需先确定定日镜场优化布置的目标函数。国内外学者对定日镜场布置方法，存在较大的差别。但总的来说，塔式光热电站定日镜的优化布置需要遵循一定的设计原则，本问选择以镜场的年均效率最高为设计原则，通过相应的评价准则来判断镜场布置结果的优劣。定日镜场的瞬时光学效率受太阳位置影响，在每个时间点都是不同的。为了综合评价定日镜场的聚光性能，可采用“典型日平均法”或“月平均法”作为定日镜场年均光学效率的评价指标7。“典型日平均法”指分别在典型日(春分日、夏至日、秋分日、冬至日)内的各时间点计算定日镜场的瞬时光学效率再求平均值以近似代表定日镜场的年均光学效率。计算公式如下：()221144,11ddttfielde fhel kkkttttt=(40)k 代表典型日序号，t1与 t2分别代表镜场的开闭场时间。“月平均法”指按一定的方式在典型年的 12 个月中分别选出 1 天作为本月的代表日。在代表日内分别计算各时间点的瞬时光学效率再求平均值以近似代表定日镜场的年均光学效率。计算公式如下：()22111212,11ddttfielde mhel mmmttttt=(41)m 代表月序号，t1与 t2分别代表镜场的开闭场时间。“典型日平均法”和“月平均法”均能代表定日镜场的综合聚光效率，不同之处在于“月平均法”计算的样本天数更大，其计算结果更接近真实情况，因此本文中采用“月平均法”计算定日镜场的年均光学效率。粒子群优化算法(Particle Swarm Optimization,PSO)于1995年由EBERHAR和KENNED提出，粒子群算法要解决的大部分问题就是最优化问题。最优化问题通俗来讲就是指在大量众多甚至无数的可行方案王秉荃 等 DOI:10.12677/aam.2024.134168 1792 应用数学进展 中挑选出一种最优的解决方案，或者说花费最小的代价取得最好的结果8。粒子群优化算法的流程如下：STEP1：根据实际情况确定镜场参数；STEP2：根据步骤一确定的参数，生成初排镜场；STEP3：根据步骤二初排镜场，用粒子群优化算法对镜场中每台定日镜的坐标位置进行优化计算。maxNNfieldiiiiieDNIAA=(42)()()()222828265,lwlwijijlHHhHHxxyyHij ij+为定日镜编号和 (43)其中，Hl为镜面宽度，Hw为镜面高度，h 为安装高度。5.2.2.多维粒子群优化算法多维粒子群优化算法 多维粒子群优化算法流程图如图 7 所示。Figure 7.Flow sheet of multi-dimensional particle swarm optimiza-tion alogrithm 图图 7.多维粒子群优化算法流程图 5.3.模型求解模型求解 在第一问求解的基础上，先用 Campo 布置方法，以简化计算量，建立仿真定日镜场。先得到定日镜场年工作效率的最大值，该值由收敛曲线直观得出。结果如图 8，可见迭代 70 次以上时定日镜场输出热王秉荃 等 DOI:10.12677/aam.2024.134168 1793 应用数学进展 功率趋于稳定，为 81 MW。Figure 8.Convergence curve of particle swarm optimizationalo-grithm 图图 8.粒子群算法收敛曲线 然后使用 MATLAB 做出二次优化前的定日镜单位面积年平均输出功率和定日镜年平均输出功率图，如图 9、图 10。可见南半边输出功率相对北半边低。目前已得到定日镜热输出功率最大时的效率云图，随即将功率限制为 60 MW，得到新的定日镜单位镜面面积年平均输出功率图和定日镜年平均输出功率图，如图 11、图 12。肉眼可见南半边有很多定日镜并不工作。同时可求出定日镜数量，为 8362 片。由于优化后的定日镜总数 8362 片，但实际很多定日镜并不提供功率，因此筛选掉不工作的定日镜数量，实际工作定日镜数量为 5128 片。优化结果见表 4表 6。Figure 9.Annually average output power per unit area of helioscope 图图 9.定日镜单位面积年平均输出功率图 王秉荃 等 DOI:10.12677/aam.2024.134168 1794 应用数学进展 Figure 10.Annually average output power of helioscope 图图 10.定日镜年平均输出功率图 Figure 11.Optimized annually average output power per unit area of helioscope 图图 11.优化后的定日镜年平均输出功率图 Figure 12.Optimized annually average output power of helioscope 图图 12.优化后定日镜单位镜面面积年平均输出功率图 王秉荃 等 DOI:10.12677/aam.2024.134168 1795 应用数学进展 Table 4.Average optical efficiency and output power of 21st every month in question 2 表表 4.问题 2 每月 21 日平均光学效率及输出功率表 日期 平均光学效率 平均余弦效率 平均阴影遮挡效率 平均