CN103500277B - 麦田型塔式太阳能热电系统的镜场优化设计方法 - Google Patents
麦田型塔式太阳能热电系统的镜场优化设计方法 Download PDFInfo
- Publication number
- CN103500277B CN103500277B CN201310456386.2A CN201310456386A CN103500277B CN 103500277 B CN103500277 B CN 103500277B CN 201310456386 A CN201310456386 A CN 201310456386A CN 103500277 B CN103500277 B CN 103500277B
- Authority
- CN
- China
- Prior art keywords
- heliostat
- jing chang
- cornfield
- mirror field
- efficiency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000003287 optical effect Effects 0.000 claims abstract description 31
- 230000005540 biological transmission Effects 0.000 claims description 9
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 238000005286 illumination Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000005619 thermoelectricity Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 241000127225 Enceliopsis nudicaulis Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000004035 construction material Substances 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Landscapes
- Mounting And Adjusting Of Optical Elements (AREA)
- Optical Elements Other Than Lenses (AREA)
Abstract
本发明公开了一种麦田型塔式太阳能热电系统镜场优化设计方法,该方法以镜场的综合光学效率与镜场的土地覆盖率的乘积作为优化目标,采用单纯形算法对镜场中第一排定日镜到吸热器的距离、各定日镜之间的列间距和行间距进行优化,使得定日镜的数量与排布在某个特定的情况下,优化目标取值最大。本发明的镜场优化设计方法解决了定日镜数量与镜场综合光学效率之间的矛盾,且采用单纯形算法具有更好的鲁棒性和稳定性。通过本发明的镜场优化设计方法进行布局的麦田型小镜场的综合光学效率都达到了80%以上,相对于现有的一般镜场提高了10%。
Description
技术领域
本发明涉及太阳能技术领域,尤其涉及麦田型塔式太阳能热电系统的镜场优化设计方法。
背景技术
塔式系统又称集中式系统,主要包括若干定日镜、吸收塔和安装于吸收塔上的吸热器。塔式太阳能热发电技术的工作原理是通过一定数量的定日镜将太阳光汇聚到塔顶的吸热器产生高温,再加热流经吸热器内的介质将光能转化成热能,并产生高温蒸汽,推动汽轮机进行发电,是一种大面积大规模的聚光发电方式。每台定日镜都各自配有跟踪机构准确地将太阳光反射集中到一个高塔顶部的吸热器上。
在整个发电系统中,镜场的排布类型和规模会直接影响镜场的光学效率和发电量,同时也会影响镜场的建造成本,因此如何合理优化镜场布局是镜场设计的关键,而如何快速计算镜场的光学效率是设计优化的前提。镜场年均综合光学效率的指标主要包括余弦效率、阴影遮挡效率、大气透射效率和溢出效率。综合光学效率为余弦效率、阴影遮挡效率、溢出效率、大气透射效率的乘积,即:
ηfield=ηcos×ηS&B×ηint×ηAt.M
其中,ηfield为综合光学效率,ηcos为余弦效率,ηS&B阴影遮挡效率,ηint为溢出效率ηAt.M为大气透射效率。
太阳光照射到定日镜表面时,和定日镜表面产生一定的角度,将入射光线的入射向量和定日镜的法向量的夹角的余弦值定义为余弦效率。如图1所示,太阳2的光线照射到定日镜1表面时,和定日镜1表面产生一定的夹角θ即为入射光线的入射向量和定日镜的法向量的夹角,即余弦角,其中坐标原点为吸收塔在地面的中心,-X表示正东,Y表示正南。余弦效率的计算只需要确定定日镜的坐标、太阳位置和吸热器位置即可求取。
在从定日镜反射至吸热器的过程中,太阳光因在大气中的衰减所导致的能量损失称为大气透射(或大气衰减)损失,而经传播后的反射光线强度与反射初的光线强度的比值即为镜面的大气透射效率。该效率通常与太阳的位置,当地海拔高度以及大气条件(如灰尘、湿气、二氧化碳的含量等)等因素有关。
余弦效率和大气透射效率与定日镜间相对位置无关,容易求取;而阴影遮挡效率和溢出效率需要考虑到镜面与镜面之间、镜面与吸热器之间的位置和角度关系,求取过程较为复杂。
如图2所示,当入射光线照射到定日镜11时,被定日镜12遮挡则造成了阴影M,产生阴影损失,被阴影的光线数量与照到该定日镜的总光线数量之比为阴影效率;同理,当反射光线被定日镜3遮挡时则造成了遮挡N,产生遮挡损失,被遮挡的光线数量与照到该定日镜的总光线数量之比为遮挡损失。一面定日镜可能同时发生被阴影和被遮挡的情况。
当反射光线由于受到吸热器尺寸的限制而没有照进吸热器内,称这些光线为溢出的光线,将溢出的光线与总的光线数之比为溢出效率。
对于镜场优化设计的方案,有最传统的将接受能量和投资成本作为优化目标的,但对于投资成本中包括的占地面积成本、导线成本、接收塔成本等具体参数和建造材料、地域位置等其他相关因素,不易统一估算。但是定日镜数量与镜场综合光学效率存在矛盾关系,若土地面积固定,当定日镜数量较多、排布较密集的时候,会产生较多的阴影遮挡损失,镜场整体的综合光学效率会随之降低;而当定日镜数量较少、排布稀疏的时候,镜场的综合光学效率会相对增加。
麦田型镜场是多面定日镜交叉排布,是目前常用的镜场的排布类型之一。确定了第一排定日镜到吸热器的距离、各定日镜之间的列间距和行间距即可确定麦田型镜场中各定日镜的排布。
发明内容
本发明针对上述现有技术中的不足,提供了一种麦田型塔式太阳能热电系统的镜场优化设计方法。
一种麦田型塔式太阳能热电系统的镜场优化设计方法,包括以下步骤:
(1)确定第一排定日镜到吸热器的距离、各定日镜之间的列间距和行间距为优化变量,并分别确定各优化变量的搜索上下界;
(2)以综合光学效率与土地覆盖率的乘积作为优化目标,根据太阳圆盘模型、镜场参数和各优化变量的搜索上下界,采用单纯形算法进行变量优化,得到各个优化变量的最优值,所述的最优值使优化目标取最大值;
(3)根据步骤(2)得到的各优化变量的最优值,进行麦田型塔式太阳能热电系统的镜场布局。
本发明的麦田型塔式太阳能热电系统的镜场优化设计方法通过将镜场的综合光学效率与镜场的土地覆盖率(即定日镜的总面积与给定土地面积之比)的乘积作为优化目标,采用单纯形算法对镜场中第一排定日镜到吸热器的距离、各定日镜之间的列间距和行间距进行优化,使得定日镜的数量与排布在某个特定的情况下,优化目标取值最大。本发明的镜场优化设计方法将镜场的综合光学效率与镜场的土地覆盖率的乘积作为优化目标,有效解决了定日镜数量与镜场综合光学效率之间的矛盾,从而进行镜场的优化布局。单纯形搜索是不基于梯度搜索的优化算法,能较快收敛到最优设定值而具有较良好的效率和可用性,且单纯形算法又是基于排序的算法,具有更好的鲁棒性和稳定性。
所述步骤(1)中各定日镜之间的列间距和行间距的搜索下界为使定日镜在三维空间内随意转动不发生机械碰撞的安全距离。
所述步骤(1)中各定日镜之间的列间距和行间距的搜索上界保证定日镜之间不发生遮挡。
定日镜在工作过程中需要转动,对自由度要求较高,通过限定优化变量的搜索上下界能够有效地减小变量在优化过程中的搜索范围,提高运算效率。
所述步骤(2)中的太阳圆盘模型采用蒙特卡洛撒点法模拟其能流密度建立得到。
作为优选,所述步骤(2)中的太阳圆盘模型为非平行入射光模型。
太阳是整个热电系统的能量来源,在精确的太阳模型中,光线是发散的,太阳表面各点向任意方向都有光线射出,因此在地球上观察太阳的像为圆形,被称为太阳圆盘(Solar Disk),该像上能量分布并不均匀。若利用蒙特卡洛撒点法来模拟太阳圆盘能量的大小,则越靠近圆心的位置,点数越密集,相应的能量高,靠近圆盘边缘则点数就越稀少,相应的能量低。根据太阳能流密度分布与太阳圆锥角的关系,结合具体时间具体地点的太阳位置,生成非平行入射光的单位向量。太阳光为非平行光,现有技术中建立的太阳模型多为平行入射光模型,本发明通过蒙特卡洛撒点法模拟其能流密度建立的太阳圆盘模型为非平行入射光的模型,更加接近实际的太阳光。根据太阳能流密度分布与太阳圆锥角的关系,结合具体时间具体地点的太阳位置,生成非平行入射光的单位向量,根据入射向量和定日镜的中心坐标可计算出定日镜中心和顶点沿着入射光线在地面上的投影坐标,并可根据入射光向量与每面定日镜的反射向量计算出相应的单位法向量。
所述步骤(2)中采用蒙特卡洛光线追迹法计算镜场的综合光学效率。
蒙特卡洛光线追迹法能有效追踪判断每根光线与每面定日镜的位置关系,同时也能获得吸热器表面的能流密度分布情况,因此该方法能有效计算阴影遮挡效率和溢出效率。
采用蒙特卡洛光线追迹法计算镜场的阴影遮挡效率和溢出效率的主要思想如下:
撒点:利用蒙特卡洛法在确定范围的定日镜场内随机投撒大量的点,由任一点和入射向量组成一根光线。
阴影判断阶段:由随机投撒在地面上的点和定日镜的四个顶点沿着入射光线在地面上的投影坐标依次确定入射光线与每一面定日镜的相交情况,若光点投影不在任何一面定日镜内,则考虑下一个点,该点所对应的光线为无效入射光线;否则判定该根入射光线最终照射到的定日镜,并求取该根入射光线与该定日镜的交点;
遮挡判定阶段:先计算该根入射光线所对应的反射光线,判断反射光线是否被其它定日镜所遮挡,若不被遮挡,则求取反射光线最终与吸热器的交点;
溢出判定阶段:判断该反射光线和吸热器的交点是否在吸热器内,若不在吸热器内则为溢出。
全部光线处理完毕后,即可计算镜场的阴影遮挡效率和溢出效率。
作为优选,所述步骤(2)基于CUDA计算平台进行变量优化。
CUDA计算平台利用GPU的双层并行结构实现其高性能并行运算,让多个线程同时计算和判断,有效地提高了运算速率,为镜场布局的优化设计奠定了基础,且实现简单,成本较低。
本发明的塔式太阳能热电系统的镜场优化设计方法将镜场的综合光学效率与镜场的土地覆盖率(定日镜的总面积与给定土地面积之比)的乘积作为优化目标,使得定日镜的数量与排布在某个特定的情况下,优化目标达到最大值,有效解决了定日镜数量与镜场综合光学效率之间的矛盾,从而完成镜场的优化设计。采用单纯形算法进行优化计算,有效提高了计算效率和可用性,且单纯形算法是基于排序的算法,具有更好的鲁棒性和稳定性。通过本发明的优化设计方法优化得到的小镜场的年均综合光学效率都达到了80%以上,与一般的塔式太阳能热电系统镜场的综合光学效率(70%)相比,提高了10%。
附图说明
图1为塔式太阳能镜场的大地坐标系示意图(-X轴指向正东,Y轴水平指向正南,Z轴指向天顶);
图2为阴影和遮挡示意图;
图3为本发明的麦田型塔式太阳能热电系统的镜场优化设计方法的流程图;
图4为麦田型镜场定日镜之间无遮挡俯视图;
图5为麦田型镜场定日镜之间无遮挡侧视图;
图6为太阳光锥角示意图;
图7为太阳圆盘能流密度示意图;
图8为非平行入射向量组成示意图;
图9为定日镜双轴固定示意图;
图10为定日镜坐标系示意图;
图11为单纯形算法示意图;
图12为单纯形初始化序贯摄动示意图;
图13为麦田型镜场优化示意图。
具体实施方式
下面将结合具体实施例对本发明做进一步说明。
本实施例中麦田型塔式太阳能热电系统位于北半球。本实施例中镜场参数如表1所示:塔高30m,定日镜的尺寸为1m×1m,土地面积为15m×15m,镜场所处纬度为北纬30°。
表1 镜场参数
镜场参数 | 具体取值 |
吸收塔高 | 30m |
定日镜的尺寸() | 1m×1m |
土地面积 | 15m×15m |
镜场所处纬度 | 北纬30° |
本实施例的麦田型塔式太阳能热电系统的镜场优化设计方法,如图3所示,具体包括以下步骤:
(1)确定第一排定日镜到吸热器的距离R0、各定日镜之间的列间距dx和行间距dy为优化变量,并分别确定各优化变量的搜索上下界;
本实施例中各优化变量的搜索上下界如下:
(a)吸热器到第一排定日镜的距离R0的搜索上下界
本实施例中R0的搜索下界和搜索上界分别为0.75Ht和1.5Ht,即搜索上下界为[0.75Ht,1.5Ht],其中Ht为吸热器中心所在的高度。
(b)列间距dx的搜索上下界
定日镜之间的最小列间距(即列间距的搜索下界)应保证定日镜在三维空间内可随意转动而不发生机械碰撞,这里将定日镜之间不发生碰撞的距离定义为安全距离,其表达式为:
其中,lm是定日镜的长,wm是定日镜的宽,dS是安全距离系数,本实施例中dS=0.1m。
为了减少算法的搜索时间,定日镜间的最大行间距和最大列间距(即搜索上界)应保证定日镜之间无遮挡。如图4所示,以定日镜A为例,首先判断可能对其造成遮挡损失的定日镜。在镜场中建立一个矩形K,使定日镜A的中心位于矩形边L1的中点位置,并使另外一边L2与太阳入射光在x-y平面的分量平行。图中的等价圆rA,rB和rC分别表示定日镜A、定日镜B、定日镜C转动时,定日镜的端点形成的轨迹。矩形两组边长的确定方式如下:假设定日镜B的镜面边缘的投影刚好和定日镜A的边缘相切,此时过定日镜B的镜面中心做到L1的垂线,交点到定日镜A的镜面中心点距离为l,当定日镜A、B均和阳光在x-y平面的分量垂直,即阳光分别与定日镜A、B等价圆相切时,l取最大值lm,如果l大于lm,则B任何时刻都不会遮挡镜面A,因此令L1=2lm。
因此,定日镜之间的列间距的搜索上下界是[DM,2lm]。
(c)行距离dy的搜索上下界
定日镜之间的最小行间距(即行间距的搜索下界)为安全距离DM。
接着确定它的搜索上界值:图5为图4的侧视图,假设通过定日镜B的镜面上边缘的太阳入射光刚好和定日镜A的镜面下边缘相切,为求L2的最大值我们假设定日镜A、定日镜B和太阳入射光垂直,即定日镜A的等价圆rA和定日镜B的等价圆rB刚好和太阳入射光在不同侧相切,此时有(0.5wm)/(0.5L2)=sin(hs),即L2=wm/sin(hs),hs是太阳的高度角。
因此,定日镜之间的行间距的搜索上下界是[DM,max(DM,wm/sin(hs))]。
(2)以综合光学效率与土地覆盖率的乘积作为优化目标,根据太阳圆盘模型、镜场参数和各优化变量的搜索上下界,采用单纯形算法进行变量优化,优化目标取最大值时优化变量的取值作为各个优化变量的最优值;
(S1)生成镜场
根据单纯形算法的初始化序贯摄动法,对各个优化变量进行初始化,并根据各个优化变量的初始值和镜场参数(土地面积大小和定日镜的规格)确定定日镜的数量,从而进一步计算每一面定日镜的中心坐标与吸热器的坐标,生成镜场。
本实施例以吸热器所在的地面为中心,从该点向西为X轴正半轴方向,向正南为Y轴正半轴方向,建立右手直角坐标系如图1所示。
根据第一排定日镜与吸热器之间的距离R0、列间距dx和行间距dy的取值,以及给定的土地面积大小和定日镜的尺寸,计算得到镜场中定日镜的数量,镜场中定日镜的列数JMAX和行数IMAX的计算公式如下:
其中,FieldLength为土地的长,FieldWidth为土地的宽,floor是向下取整。
其中,第j行,第i列的定日镜的坐标表达式为:
偶数行(从0行开始计数):
奇数行:
计算得到镜场中所有定日镜的中心坐标O(x,y,z),并生成镜场。
(S2)根据太阳圆盘模型,采用蒙特卡洛光线追迹法计算所生成镜场的综合光学效率
其中,太阳圆盘模型采用蒙特卡洛撒点法模拟其能流密度建立得到,为非平行入射光的太阳圆盘模型。
本实施例中采用如下公式:
建立太阳圆盘模型。其中:
f(α)表示的是太阳圆盘(Solar Disk)上任意一点的能流密度;
λ=0.5138为常量;
S0同f(α)单位一样,为W/m2,其数值取决于观察点到太阳的距离;
αs是太阳光锥角α0的一半,如图6所示,太阳光锥角α0指观察点Q到太阳圆盘的切线P1Q和P2Q的夹角(P1,P2分别表示切点),αs表示观察点Q到太阳圆盘的切线P1Q或P2Q与Solar Disk中心点O1到观察点Q的连线O1Q的夹角。根据几何计算可以得到:αs≈(R1/R2)/2,其中,太阳半径R1=1.39×106km,太阳与地球的距离R2=1.5×108km,代入数值计算可得:αs≈(1.39×106/1.5×108)/2≈4.6mrad≈0.5°。
f(α)表示的是太阳圆盘(Solar Disk)上任意一点P的能流密度,α为连线PQ与连线O1Q的夹角,当α大于αs时,f(α)等于0,即Solar Disk外部是黑体。
根据公式(8),采用轮盘赌原则随机产生点来模拟太阳能流密度,如图7所示,其中αs被均匀等分为100份。则生成的非平行入射光的单位向量如图8所示,表达式为:
其中,向量描述太阳圆盘平面上点的周向位置(0°~360°);向量的含义是当入射光为平行光(即太阳锥角α=0°)时的单位向量,其表达式可由太阳的高度角hs和方位角ps唯一确定:
太阳圆盘模型根据太阳高度角hs和方位角ps(即太阳光线与正南方的夹角,逆时针为负,顺时针为正)唯一确定,整个优化过程中所使用的太阳圆盘模型相同。
(S21)计算镜场的相关参数
(a)计算定日镜的反射向量与平面法向量
定日镜的反射向量是由吸热器中心到定日镜中心的方向向量,在本实施例中统一将定日镜的反射向量与平面法向量化成单位向量,镜面单位法向量可由入射光线和反射光线的向量和表示。则任意一面定日镜的反射向量和法向量分别为:
其中,|R|是定日镜中心到吸热器中心向量的模;|N|为镜面法向量的模;(0,R0,Ht)为吸热器的中心坐标。
(b)计算入射光线的单位反射向量
对任何单位入射向量它所对应的单位反射向量可由该入射向量所在镜面的法向量和表示。其表达式为:
(c)计算定日镜的顶点坐标
定日镜一般为双轴定日镜。如图9所示:一条纵向支撑轴轴A,连接定日镜镜面中心O,将定日镜固定在地面上某一位置,在定日镜背面有一条横向支撑轴轴B,轴B过镜面中心O且与镜面法向量垂直,自旋时,定日镜的镜面绕轴B旋转,且保持镜面中心O相对于地面固定不变。通过由定日镜的双轴构造(轴A和轴B),使得定日镜的与轴B平行的两条边始终与地面保持平行。
在本实施例中,通过坐标旋转的方式来描述顶点坐标。如图10所示,建立定日镜坐标系,以定日镜的中心为原点,以定日镜中平行于东西方向的边为X轴,以定日镜中平行于南北方向的边为Y轴,以定日镜的法线方向为Z轴,X轴正方向为西,Y轴正方向为南,Z轴正方向为垂直镜面向上。则在定日镜坐标系中,根据定日镜的尺寸(即定日镜的长lm和宽wm),确定定日镜的中心坐标O'(x,y,z)'为O'(0,0,0),四个顶点坐标分别为P1'(-lm/2,wm/2,0),P2'(lm/2,wm/2,0),P3'(lm/2,-wm/2,0)和P4'(-lm/2,-wm/2,0)。通过旋转公式将定日镜坐标系下四个顶点的坐标转换到大地坐标系下(如图1所示),则定日镜的第i个顶点在大地坐标系中可表示为:
Pi=(x,y,z)'+Rxyz·Pi',i=1,2,3,4 (11)
其中,Rxyz为旋转矩阵,根据公式:
Rxyz=Rot(-pm,z)·Rot(0,y)·Rot(hm-90,x) (12)
求取,其中hm和pm分别是每面定日镜的高度角和方位角,函数Rot(a,b),表示绕b轴旋转,旋转角度为a,则有:
(d)计算定日镜中心和各顶点沿着入射光线在地面上的投影坐标
在三维空间中,相当于计算直线与平面的交点,则定日镜中心和四个顶点在在地面上的投影坐标分别为(xos,yos)和(xos,yos),其中:
xos=-(a/b)·z+x,yos=-(b/c)·z+y (14)
xpsi=-(a/b)·zpi+xpi,ypsi=-(b/c)·zpi+ypi (15)
其中,i=1,2,3,4分别代表四个顶点序号,xpi,ypi,zpi分别为定日镜的第i个顶点在大地坐标系中的x轴、y轴和z轴的坐标。
(S22)利用蒙特卡洛光线追迹法计算所生成镜场的综合光学效率
将余弦效率、大气透射效率、阴影遮挡效率和溢出效率4者相乘即得到所生成镜场的综合光学效率。
余弦效率和大气透射效率只需要确定镜面中心坐标,太阳位置以及接收塔位置即可求取,具体计算方法可参见公开号为CN102519152A的专利文献。
阴影遮挡效率和溢出效率与定日镜的相对位置有关,利用蒙特卡洛光线追迹法分别计算所生成镜场的阴影遮挡效率和溢出效率,具体计算过程如下:
(a)确定投撒范围,投撒随机光点
动态跟踪定日镜顶点和中心坐标沿着入射光在地面上的投影区域,确定投撒范围,利用蒙特卡洛法随机投撒光点且投撒光点数量应尽量多,多到足可以覆盖镜场的所有镜面在地面上形成的投影区域,在地面上的投影光点的坐标可表示为:
其中,xmin,xmax分别为定日镜顶点坐标沿着入射光在地面上投影的最小和最大X坐标值,ymin,ymax分别为定日镜顶点坐标沿着入射光在地面上投影的最小和最大Y坐标值,N为投撒光点数,rand(N,1)表示随机生成N个(0,1)之间的随机数。
(b)利用光线追迹法,追踪每一根光线的运动轨迹,判断每一根光线和每一面定日镜或吸热器的相交情况,有效计算阴影遮挡效率和溢出效率。
(S3)计算土地覆盖率,并进一步计算得到综合光学效率和土地覆盖率的乘积,即优化目标。
(S4)利用单纯形优化方法重复步骤S1~S3,将得到的优化目标的大小与前一次的进行比较,保留较大值,多次重复,得到最大优化目标,并以优化目标最大时优化变量的取值作为相应优化变量的最优值,即变量优化的结果。
步骤(S4)中多次循环进行步骤S1~S3,在循环过程中每次优化变量的取值不同,每次生成的镜场中定日镜的数量,以及定日镜的列数JMAX和行数IMAX均也在不断优化改变。
采用经过适当修改的单纯形搜索算法进行镜场优化设计,如图11所示,算法的具体思路如下:
(2-1)算法初始化。单纯形计数s=0,单纯形构建算子v=1,初始化单纯形搜索算法参数{α,β,γ,δ},其中α为反射因子,β为收缩因子,γ为膨胀因子,δ为塌缩因子;
(2-2)初始化判定条件:如果当前v≥n+1(n为优化变量数),则转步骤(2-4);否则返回执行步骤(2-3);
(2-3)初始单纯形构建。假设第K个优化变量的搜索上界和搜索下界分别为和(本实施例中优化变量为第一排定日镜到吸热器的距离R0、各定日镜之间的列间距dx和行间距dy为优化变量),则对初始点采用序贯摄动法构建初始单纯形,设τ为摄动因子,若则向右摄动反之,则向左摄动令k=k+1,v=v+1,初始单纯形构建的摄动示意图如图12所示。
(2-4)单纯形顶点排序。将单纯形的顶点所生成的镜场相对应的响应值,即优化目标:镜场的年均光学率和土地覆盖率的乘积,的大小进行排序,V1 s表示最小点,表示距离最大点,表示次大点。令Vs+1=Vs,Fs+1=Fs,s=s+1;
(2-5)反射:根据产生反射工作点其中令k=k+1,如果Yref<F1 s,执行步骤(2-6)执行膨胀操作;如果转步骤(2-7),执行收缩操作;其他情况下,用Vref代替Vn+1,Yref代替返回执行步骤(2-4);
(2-6)膨胀:根据产生膨胀点k=k+1,
若Yexp≤Yref,用代替 返回执行步骤(2-4);
否则,用代替 返回执行步骤(2-4);
(2-7)收缩:根据式产生收缩点其中当否则,令k=k+1,收缩后,比较收缩点和收缩参考点如果用代替 返回执行步骤(2-4);否则,执行步骤(2-8)执行塌缩操作,设v=2;
(2-8)塌缩:执行塌缩操作k=k+1,v=v+1,当v≥n+1,返回执行步骤(2-4);否则继续步骤(2-8)。
通过优化得到各优化变量的最优值,如表2所示。
(3)根据变量优化的结果,进行镜场布局,如图13所示,即令第一排定日镜到吸热器的距离R0=23m、各定日镜之间的列间距dx=1.81m和行间距dy=1.51m。
表2 变量优化结果
优化变量 | 最优值 |
R0 | 23m |
dx | 1.81m |
dy | 1.51m |
结果表明,经过优化后的小镜场的综合光学效率达到了88.8%,与一般高于塔式太阳能热电系统的镜场的综合光学效率(约70%)比较,提高了约18%,且根据本实施例的镜场优化设计方法布局的镜场的土地覆盖率为32%。
Claims (8)
1.一种麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,包括以下步骤:
(1)确定第一排定日镜到吸热器的距离、各定日镜之间的列间距和行间距为优化变量,并分别确定各优化变量的搜索上下界;
(2)以综合光学效率与土地覆盖率的乘积作为优化目标,根据镜场参数、太阳圆盘模型和各优化变量的搜索上下界,采用单纯形算法进行变量优化,得到各个优化变量的最优值,所述的最优值使优化目标取最大值;
所述综合光学效率为余弦效率、阴影遮挡效率、溢出效率和大气透射效率的乘积,所述土地覆盖率为定日镜的总面积与给定土地面积之比;
(3)根据步骤(2)得到的各优化变量的最优值,进行麦田型塔式太阳能热电系统的镜场布局。
2.如权利要求1所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(1)中各定日镜之间的列间距和行间距的搜索下界为使定日镜在三维空间内随意转动不发生机械碰撞的安全距离。
3.如权利要求2所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(1)中各定日镜之间的列间距和行间距的搜索上界保证定日镜之间不发生遮挡。
4.如权利要求3所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(2)中的镜场参数包括定日镜的尺寸和土地面积。
5.如权利要求4所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(2)的太阳圆盘模型采用蒙特卡洛撒点法模拟其能流密度建立得到。
6.如权利要求5所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(2)中的太阳圆盘模型为非平行入射光模型。
7.如权利要求6所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(2)中采用蒙特卡洛光线追迹法计算镜场的综合光学效率。
8.如权利要求7所述的麦田型塔式太阳能热电系统的镜场优化设计方法,其特征在于,所述步骤(2)中基于CUDA计算平台进行变量优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310456386.2A CN103500277B (zh) | 2013-09-29 | 2013-09-29 | 麦田型塔式太阳能热电系统的镜场优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310456386.2A CN103500277B (zh) | 2013-09-29 | 2013-09-29 | 麦田型塔式太阳能热电系统的镜场优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103500277A CN103500277A (zh) | 2014-01-08 |
CN103500277B true CN103500277B (zh) | 2017-01-04 |
Family
ID=49865485
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310456386.2A Active CN103500277B (zh) | 2013-09-29 | 2013-09-29 | 麦田型塔式太阳能热电系统的镜场优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103500277B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105402901A (zh) * | 2015-11-17 | 2016-03-16 | 浙江大学 | 塔式太阳能热电系统镜场在柱状接收器上的成像方法 |
CN106096085B (zh) * | 2016-05-30 | 2019-07-30 | 浙江大学 | 一种塔式太阳能热发电系统的辐射能密度模拟方法 |
CN109102121B (zh) * | 2018-08-15 | 2020-09-11 | 浙江大学 | 塔式太阳能热电站定日镜子镜倾斜角度优化方法 |
US11017561B1 (en) * | 2018-10-09 | 2021-05-25 | Heliogen, Inc. | Heliostat tracking based on circumsolar radiance maps |
CN109446650B (zh) * | 2018-10-29 | 2023-01-31 | 东南大学 | 一种多塔镜场的边界优化方法 |
CN109813754B (zh) * | 2019-02-14 | 2022-06-28 | 浙江可胜技术股份有限公司 | 一种测量与优化吸热器截断效率的系统与方法 |
CN110209207B (zh) * | 2019-05-07 | 2022-05-24 | 中国神华能源股份有限公司 | 确定定日镜的损失面积的方法和装置及机器可读存储介质 |
US11732928B1 (en) | 2019-09-09 | 2023-08-22 | Heliogen Holdings, Inc. | Gas receiver for capturing solar energy |
CN112152559B (zh) * | 2020-09-27 | 2021-12-07 | 浙江中光新能源科技有限公司 | 一种定日镜与光伏结合的镜场系统及其布置方法 |
CN113359868B (zh) * | 2021-06-23 | 2022-11-22 | 中国气象局公共气象服务中心(国家预警信息发布中心) | 塔式发电辐射在定日镜到吸收器间折减计算方法、介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930160A (zh) * | 2012-11-02 | 2013-02-13 | 浙江大学 | 塔式太阳能热电系统镜场光学效率的计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8033110B2 (en) * | 2008-03-16 | 2011-10-11 | Brightsource Industries (Israel) Ltd. | Solar power generation with multiple energy conversion modes |
-
2013
- 2013-09-29 CN CN201310456386.2A patent/CN103500277B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930160A (zh) * | 2012-11-02 | 2013-02-13 | 浙江大学 | 塔式太阳能热电系统镜场光学效率的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103500277A (zh) | 2014-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103500277B (zh) | 麦田型塔式太阳能热电系统的镜场优化设计方法 | |
CN102930160B (zh) | 确定塔式太阳能热电系统镜场光学效率的方法 | |
Li et al. | Optimization of a heliostat field layout using hybrid PSO-GA algorithm | |
CN102519152B (zh) | 基于gpu塔式太阳能热电系统镜场光学效率的计算方法 | |
Collado et al. | A review of optimized design layouts for solar power tower plants with campo code | |
Yu et al. | Analysis and improvement of solar flux distribution inside a cavity receiver based on multi-focal points of heliostat field | |
He et al. | Numerical simulation of solar radiation transmission process for the solar tower power plant: from the heliostat field to the pressurized volumetric receiver | |
Yao et al. | Heliostat field layout methodology in central receiver systems based on efficiency-related distribution | |
CN105469160B (zh) | 塔式太阳能的扇形定日镜场布置方法 | |
CN109102121B (zh) | 塔式太阳能热电站定日镜子镜倾斜角度优化方法 | |
CN104034058B (zh) | 基于gpu的塔式太阳能热电系统镜场的成像方法 | |
CN106369838B (zh) | 一种槽式聚光太阳能集热系统设计方法 | |
CN103530697A (zh) | 辐射型塔式太阳能热电系统的镜场优化设计方法 | |
CN111881576A (zh) | 太阳能塔式光热电站定日镜场优化调度控制方法 | |
CN102877457B (zh) | 太阳能定向线聚光沙漠表面硬化系统 | |
Deng et al. | Rose pattern for heliostat field optimization with a dynamic speciation‐based mutation differential evolution | |
Sánchez-González et al. | FluxSPT: Tool for heliostat field aiming and flux mapping in solar power tower plants | |
CN103530518B (zh) | 太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法 | |
Meng et al. | A vector based freeform approach for reflecting concentrator of solar energy | |
CN202809597U (zh) | 太阳能定向线聚光沙漠表面硬化系统 | |
CN104639025A (zh) | 密集阵列式聚光光伏系统 | |
CN105402901A (zh) | 塔式太阳能热电系统镜场在柱状接收器上的成像方法 | |
CN102877459A (zh) | 偏置式太阳能定向线聚光沙漠表面硬化系统和硬化方法 | |
Bu et al. | Optimization Model of Heliostat Field Layout based on Genetic Algorithm | |
CN110209209A (zh) | 优化太阳能定日镜场的方法和装置及机器可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |