CN109670248B - 塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 - Google Patents
塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 Download PDFInfo
- Publication number
- CN109670248B CN109670248B CN201811578452.2A CN201811578452A CN109670248B CN 109670248 B CN109670248 B CN 109670248B CN 201811578452 A CN201811578452 A CN 201811578452A CN 109670248 B CN109670248 B CN 109670248B
- Authority
- CN
- China
- Prior art keywords
- heliostat
- energy density
- density distribution
- radiant energy
- power generation
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000005855 radiation Effects 0.000 title claims abstract description 37
- 238000004088 simulation Methods 0.000 title claims abstract description 35
- 238000010248 power generation Methods 0.000 title claims abstract description 20
- 230000014509 gene expression Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 6
- 238000002310 reflectometry Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims 2
- 241000764238 Isis Species 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 9
- 230000000052 comparative effect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 230000004907 flux Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 241000321453 Paranthias colonus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，涉及塔式太阳能热发电系统模拟技术领域，包括以下步骤：(1)在考虑阴影和遮挡的影响的情况下，提取定日镜表面的有效反射区域；(2)推导和建立定日镜局部坐标系下定义的虚拟辐射能密度分布标量场函数的解析模型；(3)将虚拟的所述虚拟辐射能密度分布标量场函数沿着定日镜反射的方向斜平行投影到接收平面上，完成接收器上接收到的辐射能密度分布建模；本发明考虑了太阳光方向、太阳表面的辐射能分布、定日镜的位置、大小、旋转、定日镜微表面以及阴影和遮挡等因素，实验和对比结果表明本发明提出的解析模型比已有的解析方法精度都要高并且可以快速的计算。
Description
技术领域
本发明涉及塔式太阳能热发电系统模拟技术领域，特别涉及塔式太阳能热发电中光斑的辐射能密度分布的模拟方法。
背景技术
塔式太阳能热发电是一种既环保又可持续的应对当代世界能源危机的新兴科技(Romero M，Steinfeld A.Concentrating solar thermal power and thermochemicalfuels[J].Energy&Environmental Science，2012， 5(11)：9234-9245.)。在该类型的系统中，操控策略，尤其是当系统运行不稳定的时候，需要大量的及时辐射能仿真模拟。除了发电效率的考虑，接收器表面的辐射能密度分布关系到系统的安全运营(Roldán M I，Monterreal R.Heat flux and temperature prediction on a volumetric receiverinstalled in a solar furnace[J].Applied Energy，2014， 120(3)：65-74.)，因此，对接受器表面接收到的来自于全定日镜场的辐射能分布进行精确快速的仿真模拟就变得非常重要和意义重大，但前提是要对一面定日镜反射的辐射能密度进行建模。
但是，对定日镜反射的辐射能光斑进行精确的模拟并不是一件容易的事情，因为这个光斑受很多因素的影响，比如定日镜位置，大小，朝向以及定日镜微表面，太阳方向以及太阳表面的辐射能密度分布。另外，场地上定日镜之间的阴影和遮挡损失也需要重点考虑。现有的仿真方法可以分为两类，基于光线跟踪的方法和解析的方法(Garcia P，Ferriere A， Bezian J J.Codes for solar flux calculation dedicated to centralreceiver system applications：A comparative review[J].Solar Energy， 2008，82(3)：189-197.)。光线跟踪方法通过离散地采样并跟踪大量的光线来模拟光能在太阳、定日镜和接收器之间的传输。这种方法可以得到相对精确的辐射能密度分布的预测结果，伴随的代价是计算复杂度比较高。另一方面，解析方法通过对辐射能光斑进行整体的刻画和近似，从而极大地降低计算时间。一般而言，解析方法刻画分布模型通过卷积的方式，关联因素包括太阳、定日镜、接收器。但是，这种卷积的结果到如今依旧没有一个解析表达式。卷积模型通常通过数值方法或者对模型简化进行求解，前一种方式比较耗时，后一种方式仿真精度低。
模拟定日镜反射到接收器上的辐射能密度分布，解析方法是一种有效途径，涵盖多种影响因素，比如太阳参数，定日镜属性和大气衰减等。Biggs 和Vittitoe(1976)提出通过卷积的方式在统计学意义上综合多个非确定因素的影响(太阳表面形状和定日镜镜片误差)，这个卷积利用傅里叶变换进行数值求解(Vittitoe C N，Biggs F.The HELIOS modelfor the optical behavior of reflecting solar concentrators[J].Sandia NationalLaboratories Report No.SAND76-0347，1976.)。Walzel等人在1977年利用二维的诶米特多项式来预测一个平面定日镜反射到接收器上的辐射能密度分布(Walzel M D，Lipps FW，Vant-Hull L L.A solar flux density calculation for a solar towerconcentrator using a two-dimensional hermite function expansion[J].SolarEnergy，1977， 19(3)：239-253.)。他们的结果跟之前的解析模型结果很接近但是只需要1/10的计算时间。Lipps和Walzel在1978年提出了一个针对圆形的定日镜的解析模型并且也是通过数值的方法刻画由任意形状的定日镜反射的辐射能密度分布，同时考虑了阴影和遮挡的影响。作者指出他们的方法比诶米特函数的方法精度更高但计算速度更慢。Hennet和Abatut在1984 年给出了一个复杂的辐射能密度分布解析模型，计算太阳圆盘和定日镜在接收面上的投影轮廓的卷积(Hennet J C，Abatut J L.An analytical method forreflected flux density calculations[J].Solar energy， 1984，32(3)：357-363.)。阴影的处理方式很直接，根据被遮挡的面积占比统一均匀降低反射的辐射能。显然，这种做法对于辐射能密度分布的预测不是很准确。Collado等人在1986年以积分的形式给出了UNIZAR模型来刻画聚焦型定日镜在接收器上形成的辐射能密度分布，本质上就是几个高斯分布的卷积(Collado F J，Gomez A，Turégano J A.An analytic function for the fluxdensity due to sunlight reflected from a heliostat[J].Solar Energy，1986，37(3)：215-234.)。他们通过比较模拟的辐射能密度分布等高线、总能量与真实测量的结果来验证提出的解析模型。总之，一个普适的具有闭合形式的基于卷积积分方式的辐射能密度分布模拟方案还至今未知。
Elsayed等人在1995年设计了一个仪器测量了一个平面定日镜反射在平面接收面板上的辐射能密度分布(Elsayed M M，Fathalah K A， A1-Rabghi O M.Measurements ofsolar flux density distribution on a plane receiver due to a flat heliostat[J].Solar energy，1995， 54(6)：403-411.)。根据测量的数据，他们识别出了两个图像，也就是可见光斑区域和定日镜反射的辐射能密度分布区，结果分布图跟本文提出的解析模型结果高度一致。太阳辐射能沿着定日镜在接收面上的投影轮廓的主轴有规律的分布。他们还给出了一个一维的解析模型来逼近测量的结果 (Elsayed M M，Fathalah K A.Solarflux density distribution using a separation of variables/superpositiontechnique[J].Renewable energy，1994，4(1)：77-87.)。
作为另一种方式，基于卷积的解析模型可以简化为一些具有闭合形式的简单形式表达式，代价是牺牲一点精度。Schwarzbozl等人在2009年利用一个各向同性的高斯分布模型来逼近接收器表面的辐射能密度分布，也就是HFLCAL模型(
P，Pitz-PaalR，Schmitz M.Visual HFLCAL-A Software Tool for Layout and Optimisation ofHeliostat Fields[C]//Proceedings of 15th International SolarPACES Symposium，Berlin，September.2009：15-18.)。太阳模型、定日镜微表面以及跟踪误差等因素综合为高斯模型的参数来刻画，这个模型之后跟测量的数据进行了验证(Collado F J.One-pointfitting of the flux density produced by a heliostat[J].Solar Energy，2010，84(4)：673-684.)。之后， HFLCAL模型被应用在了定日镜聚焦策略优化等相关工作(SaloméA，Chhel F，Flamant G，et al.Control of the flux distribution on a solar towerreceiver using an optimized aiming point strategy： Application to THEMIS solar tower[J].Solar Energy，2013，94： 352-366.)。Garcia等人在2015年对原始的HFLCAL模型进行了改进，与之前的不同，他们提出首先在定日镜平面上建立一个各向同性的高斯分布辐射能密度标量场，然后将这个标量场通过单应性变换(也就是斜平行投影)投影到平面接收面上(García L，Burisch M，Sanchez M.Spillage Estimation in aHeliostats Field for Solar Field Optimization[J]. Energy Procedia，2015，69(12)：1269-1276.)。这种投影变换方式更好地逼近了真实场景中接收面上由于斜投影变形了的辐射能光斑因此刻画得更精确。Huang和Sun在2016年提出用一个椭圆高斯模型来建模由聚焦行定日镜形成的光斑的辐射能密度分布，但是没有考虑阴影和遮挡的影响 (HuangW，Sun L.Solar flux density calculation for a heliostat with an ellipt icalGaussian di stribution source[J].Applied Energy，2016， 182：434-441.)。基于改进的HFLCAL模型，He等人在2017年提出了一个快速辐射能密度分布的计算框架，通过挖掘现代GPU绘制流水线的强大并行计算能力，并且考虑了阴影和遮挡(He C，Feng J，ZhaoY.Fast flux density distribution simulation of central receiver system on GPU[J].Solar Energy，2017，144：424-435.)。局限于采用的辐射能密度分布模型以及处理阴影和遮挡采用的z-buffering算法，光斑上被遮挡区域的轮廓边缘的辐射能密度分布呈现出尖锐的过渡，显然是跟真实的情况不符。
发明内容
本发明提供了塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，提高模拟精度和效率。
塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，包括以下步骤：
(1)在考虑阴影和遮挡的影响的情况下，提取定日镜表面的有效反射区域；
(2)推导和建立定日镜局部坐标系下定义的虚拟辐射能密度分布标量场函数的解析模型；
(3)将虚拟的所述虚拟辐射能密度分布标量场函数沿着定日镜反射的方向斜平行投影到接收平面上，完成接收器上接收到的辐射能密度分布建模。
为了进一步提高模拟精度和效率，优选的，步骤(1)中，在考虑阴影和遮挡的影响的情况下，提取定日镜表面的有效反射区域的具体过程如下：
1-1、确定定日镜场上所有定日镜之间的阴影遮挡关系；
1-2、会产生阴影和遮挡的定日镜分别沿着太阳能光线方向和反射方向的逆方向投影到当前考察的定日镜平面上，得到阴影和遮挡镜面的投影轮廓；
1-3、将步骤1-2得到的投影轮廓裁剪掉，得到了定日镜表面上除去被遮挡的区域，定义为定日镜有效反射表面。
为了进一步提高模拟精度和效率，优选的，步骤1-1，采用在GPU上实现的光柱遍历均匀网格算法(He C，Feng J，Zhao Y.Fast flux density distribution simulation ofcentral receiver system on GPU[J].Solar Energy，2017，144：424-435.)快速并行确定定日镜场上所有定日镜之间的阴影遮挡关系。
为了进一步提高模拟精度和效率，优选的，步骤1-3中，采用 Weiler-Atherton裁剪算法将步骤1-2得到的投影轮廓裁剪掉。
为了进一步提高模拟精度和效率，优选的，步骤(2)中，推导和建立定日镜局部坐标系下定义的虚拟辐射能密度分布标量场函数的解析模型的具体过程如下：
2-1、在定日镜局部坐标系o-xy下，定日镜有效反射表面S1是个平面多边形，顶点为
对应的边为
每一条边从vi(xi，yi)到vi+1(xi+1，yi+1)，可以表示成参数形式：
其中：
ai＝xi+1-xi，bi＝yi+1-yi，υn(xn，yn)＝υ0(x0，y0)；
2-2、对于定日镜二维平面上的任意一点p(x，y)，p(x，y)对应的虚拟辐射能密度通过均匀的IFDSF函数H(q)和归一化的准Cauchy核函数C(r)在定日镜平面上作卷积得到：
其中：
Ph表示定日镜反射的总辐射能，等于太阳参数ID、定日镜反射表面面积SH、定日镜余弦效应cosω和定日镜反射率的乘积；
ηaa是大气衰减因子；
SH′是定日镜反射表面沿着反射方向斜平行投影到接收面上的区域的面积；
s是准Cauchy核函数的参数；控制准Cauchy核的半径的大小，从而控制最终光斑的扩散范围。
因此：
结合格林公式，令：
则：
结合以上，得到：
C是有效反射表面S1的边界，也就是线段集合
令Ii和fi(t)分别表示：
2-3、求解fi(t)，得到：
其中，
得到：
2-4、注意到，归一化的准柯西核函数
是定义在整个定日镜平面上的，而接收面区域是有界的。考虑到能量守恒，对步骤3-3得到的辐射能密度分布表达式进行归一化操作，除以核函数在接收面投影到镜面上的区域的积分，记为Cg，得到点p(x，y)上对应的虚拟辐射能密度值表示为：
其中，
2-5、根据步骤2-4的解析表达式计算得到的虚拟辐射能密度分布标量场。
其中，s＝3～3.5。
为了进一步提高模拟精度和效率，优选的，步骤(3)中，将虚拟的所述虚拟辐射能密度分布标量场函数沿着定日镜反射的方向斜平行投影到接收平面上的具体过程如下：
在全局坐标系O-XYZ中，令R＝(Rx，Ry，Rz)为接收面上的一点， H＝(Hx，Hy，Hz)为在平行投影操作下在定日镜上对应R的一点，那么根据线面求交方程推导，有：
其中：
r＝(rx，ry，rz)是考察定日镜的反射光线方向；
T表示向量的转置操作；
N＝(Nx，Ny，Nz)是定日镜的法向；
o＝(ox，oy，oz)是定日镜的中心；
E是单位矩阵。
本发明的有益效果：
本发明的模拟方法给出了具有闭合表达式的解析模型来刻画平面定日镜反射到平面接收面上的辐射能密度分布，考虑了太阳光方向、太阳表面的辐射能分布、定日镜的位置、大小、旋转、定日镜微表面以及阴影和遮挡等因素，实验和对比结果表明本发明提出的解析模型比已有的解析方法精度都要高并且可以快速的计算。
附图说明
图1为本发明的模拟方法中建立的坐标系。
图2为定日镜表面虚拟的辐射能密度分布标量场的建模原理示意图，其中，图2(a)为均匀分布的密度函数示意图，图2(b)为标准准柯西核函数做卷积示意图，图2(c)为得到的光滑分布的函数示意图。
图3为定日镜场上阴影和遮挡现象的可视化示意图。
图4为在定日镜局部坐标系下剔除阴影和遮挡区域的定日镜有效反射镜面S1的示意图。
图5为定日镜局部坐标系下的虚拟辐射能密度标量场示意图。
图6为本发明为验证提出的模型的准确性设计的辐射型交错分布的大型定日镜场的示意图。
图7为中午时刻例子中光线跟踪方法(第一列)和本实施例的方法与改进的HFLCAL模型的对比等高线图(第二列)对实验定日镜反射的光斑进行模拟的结果对比示意图。
图8为傍晚时刻例子中本实施例的方法(第一列)和本实施例的方法与改进的HFLCAL模型的对比等高线图(第二列)对实验定日镜反射的光斑进行模拟的结果对比示意图。
具体实施方式
如图1～6所示，本实施例的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法分为三步，首先，考虑阴影和遮挡的影响，提取定日镜表面的的有效反射区域；接着，推导和建立定日镜局部坐标系下定义的虚拟辐射能密度分布标量场函数的解析模型；最后，这个虚拟的标量场函数沿着定日镜反射的方向斜平行投影到接收平面上，从而完成接收器上接收到的辐射能密度分布建模，其中图1为本发明的模拟方法中建立的坐标系，分别是定日镜场全局坐标系O-XYZ、定日镜局部坐标系o-xy、接收面坐标系c-uv，以及在平行投影下他们之间的关系。图2为定日镜表面虚拟的辐射能密度分布标量场的建模原理示意图，通过一个如图2(a)的均匀分布的密度函数与如图2(b)的标准准柯西核函数做卷积，从而得到一个光滑分布的函数如图2(c)。图3为定日镜场上阴影和遮挡现象的可视化示意图，红色标示的入射光柱和黄色表示的反射光柱分别被不同的相邻定日镜遮挡了一部分，从而在考察定日镜上产生了阴影遮挡现象。具体过程如下：
(1)定日镜有效反射区域的提取：
首先，定日镜场上所有定日镜之间的阴影遮挡关系通过在GPU上实现的光柱遍历均匀网格算法(He C，Feng J，Zhao Y.Fast flux density distribution simulation ofcentral receiver system on GPU[J].Solar Energy，2017，144：424-435.)快速并行确定。然后，那些会产生阴影和遮挡的定日镜分别沿着太阳能光线方向和反射方向的逆方向投影到当前考察的定日镜平面上，从而得到阴影和遮挡镜面的投影轮廓，应用 Weiler-Atherton裁剪算法就可以将这些遮挡区域裁剪掉。这样就得到了定日镜表面上除去被遮挡的区域，也称为定日镜有效反射表面。
(2)定日镜上虚拟辐射能密度标量场的解析建模：
在这一小节里，推导一个刻画虚拟辐射能密度分布的解析函数F(x，y)，它定义在定日镜局部坐标系上，它沿着定日镜反射方向平行投影到接收器表面得到的就是最终要求的辐射能密度分布结果。
如图3所示，在定日镜局部坐标系o-xy下，定日镜有效反射表面S1是个平面多边形，它的顶点为
对应的边为
每一条边从 vi(xi，yi)到vi+1(xi+1，yi+1)，可以表示成参数形式：
其中：
ai＝xi+1-xi，bi＝yi+1-yi，υn(xn，yn)＝υ0(x0，y0)
基于分析，如图3所示，对于定日镜二维平面上的任意一点p(x，y)，它对应的虚拟辐射能密度可以通过均匀的IFDSF函数H(q)和一个归一化的准Cauchy核函数C(r)在定日镜平面上作卷积得到：
其中：
这里Ph表示定日镜反射的总辐射能，它等于太阳参数DNI(ID)，定日镜反射表面面积(SH)，定日镜余弦效应(cosω)和定日镜反射率的乘积。ηaa是大气衰减因子，SH′是定日镜反射表面沿着反射方向斜平行投影到接收面上的区域的面积。
因此，
结合格林公式，令：
则：
结合以上，有：
C是有效反射表面S1的边界，也就是线段集合
令Ii和fi(t)分别表示：
求解fi(t)，得到：
其中：
所以得到：
归一化的准柯西核函数C(r)是定义在整个定日镜平面上的，而接收面区域是有界的。考虑到能量守恒，上述的辐射能密度分布表达式应该进行归一化操作，除以核函数在接收面投影到镜面上的区域的积分，记为Cg。最终，点p(x，y)上对应的虚拟辐射能密度值可以表示为：
其中：
根据上式解析表达式计算得到的虚拟辐射能密度分布标量场样例参见图5，考虑了阴影和遮挡的影响，定日镜的尺寸是2m*2m，A＝600W/m2以及s＝3.2。
(3)斜平行投影到接收面：
接收器的接收平面上的最终辐射能密度分布是将定日镜上的虚拟辐射能标量场沿着定日镜反射方向平行投射到接收器上得到的，这个操作通过斜平行投影实现。
在全局坐标系O-XYZ中，令R＝(Rx，Ry，Rz)为接收面上的一点， H＝(Hx，Hy，Hz)为在平行投影操作下在定日镜上对应R的一点，那么根据线面求交方程推导，有：
其中r＝(rx，ry，rz)是考察定日镜的反射光线方向，N＝(Nx，Ny，Nz)是定日镜的法向，o＝(ox，oy，oz)是定日镜的中心。E是单位矩阵。
本实施例方法和现有技术的对比结果如下：
如图7和8所示，通过大量的实验和跟以往相关工作的对比，证明本实施例提出的解析模型的精确性。
根据Lipps和Vant-Hull在1978年提出的算法，本实施例设计了一个辐射型交错分布的定日镜场地用于实验，如图6所示，实验定日镜对象用他们的编号标示，并且选择了两个具有代表性的太阳方位角和高度角组合(表格1)。一个对应正午时刻，另一个是下午太阳快下山高度角很小的时候，这时定日镜场地上的阴影和遮挡现象很常见。定日镜的尺寸是3.2m*2.2m，反射率是0.88。定日镜的微表面倾斜误差的标准差设为 2.0mrad。随机选取了50面定日镜作为实验，他们到接收器的距离从100m 到900m。这些实验定日镜的参数如表格2。所有定日镜都对准接收面 (12m*12m)的中心(0.0，137.0，1.0)，接收面的法向是(0.0，0.0，1.0)。在光线跟踪模拟方法和基于本文提出的解析模型方法模拟的实验中，接收器的细分粒度都一样，都是5*5cm2。
表格1实验设计的两个时刻太阳属性
| 太阳高度角(度) | 太阳方位角(度) | DNI(W/m<sup>2</sup>) | |
| 正午时刻 | 83.15 | 88.0 | 1000.0 |
| 傍晚 | 10.0 | 162.0 | 600.0 |
表格2用于对比实验的定日镜参数
三种模拟方法对于选择的实验定日镜在接收器上形成的辐射能密度分布的模拟结果和比较，分别是光线跟踪方法、本实施例方法和改进的 HFLCAL模型。
结论：本实施例提出的模型考虑了定日镜的位置、几何尺寸、微表面、旋转角、以及阴影和遮挡，因为本模型的本质是对提出阴影和遮挡区域的有效反射面做卷积得到的标量场。另一方面，HFLCAL模型用一个简单的高斯分布对辐射能光斑进行刻画，不考虑定日镜的大小以及阴影遮挡，因此，基于本文模型预测出来的辐射能光斑跟真实值更接近，无论从光斑形状还是密度分布的角度。另外，在本模型中，沿着定日镜有效反射轮廓的辐射能密度分布呈现出了从亮到暗的光滑过度的现象，这跟光线跟踪的结果以及(Elsayed M M，Fatha]ah K A，Al-Rabghi O M.Measurements of solar flux density distribution on aplane receiver due to a flat heliostat[J].Solar energy，1995，54(6)：403-411.)的报告一致。相比较之下，HFLCAL模型的结果的辐射能密度分布偏差较大，尤其是在被阴影和被遮挡区域，而且正因为没有考虑阴影和遮挡，它预测出来的接收面上接收到的总能量也偏大。
综上所述，本实施例的模拟方法给出了具有闭合表达式的解析模型来刻画平面定日镜反射到平面接收面上的辐射能密度分布，考虑了太阳光方向、太阳表面的辐射能分布、定日镜的位置、大小、旋转、定日镜微表面以及阴影和遮挡等因素，实验和对比结果表明本发明提出的解析模型比已有的解析方法精度都要高并且可以快速的计算。
Claims (9)
1.塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，包括以下步骤：
(1)在考虑阴影和遮挡的影响的情况下，提取定日镜表面的有效反射区域；
(2)推导和建立定日镜局部坐标系下定义的虚拟辐射能密度分布标量场函数的解析模型，具体过程如下：
2-1、在定日镜局部坐标系o-xy下，定日镜有效反射表面S1是个平面多边形，顶点为
对应的边为
2-2、对于定日镜二维平面上的任意一点p(x，y)，p(x，y)对应的虚拟辐射能密度通过均匀的IFDSF函数H(q)和归一化的准Cauchy核函数C(r)在定日镜平面上作卷积，获得虚拟辐射能密度的卷积积分；
2-3、对步骤2-2的卷积积分进行求解，获得虚拟辐射能密度分布；
2-4、对步骤2-3得到的虚拟辐射能密度分布进行归一化操作，获得p(x，y)对应的虚拟辐射能密度表达式；
2-5、根据步骤2-4得到的虚拟辐射能密度表达式，计算得到虚拟辐射能密度分布标量场；
(3)将虚拟的所述虚拟辐射能密度分布标量场函数沿着定日镜反射的方向斜平行投影到接收平面上，完成接收器上接收到的辐射能密度分布建模。
2.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，步骤(1)中，在考虑阴影和遮挡的影响的情况下，提取定日镜表面的有效反射区域的具体过程如下：
1-1、确定定日镜场上所有定日镜之间的阴影遮挡关系；
1-2、会产生阴影和遮挡的定日镜分别沿着太阳能光线方向和反射方向的逆方向投影到当前考察的定日镜平面上，得到阴影和遮挡镜面的投影轮廓；
1-3、将步骤1-2得到的投影轮廓裁剪掉，得到了定日镜表面上除去被遮挡的区域，定义为定日镜有效反射表面。
3.如权利要求2所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，步骤1-1，采用在GPU上实现的光柱遍历均匀网格算法快速并行确定定日镜场上所有定日镜之间的阴影遮挡关系。
4.如权利要求2所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，步骤1-3中，采用Weiler-Atherton裁剪算法将步骤1-2得到的投影轮廓裁剪掉。
5.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，所述步骤2-1中的边
是从vi(xi，yi)到vi+1(xi+1，yi+1)构成，具体参数形式如下：
其中：
ai＝xi+1-xi，bi＝yi+1-yi，vn(xn，yn)＝v0(x0，y0)。
6.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，所述步骤2-2的具体过程如下：
其中：
Ph表示定日镜反射的总辐射能，等于太阳参数ID、定日镜反射表面面积SH、定日镜余弦效应cosω和定日镜反射率ρ的乘积；
ηaa是大气衰减因子；
SH′是定日镜反射表面沿着反射方向斜平行投影到接收面上的区域的面积；
s是准Cauchy核函数的参数；
因此：
结合格林公式，令：
则：
结合以上，得到：
C是有效反射表面S1的边界，也就是线段集合
令Ii和fi(t)分别表示：
7.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，所述步骤2-3的具体过程如下：
其中，
得到辐射能密度分布的表达式：
s是准Cauchy核函数的参数。
8.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，所述步骤2-4的归一化操作是将获得的辐射能密度分布除以核函数在接收面投影到镜面上的区域的积分，具体表达式如下：
Cg是核函数在接收面投影到镜面上的区域的积分。
9.如权利要求1所述的塔式太阳能热发电中光斑的辐射能密度分布的模拟方法，其特征在于，步骤(3)中，将虚拟的所述虚拟辐射能密度分布标量场函数沿着定日镜反射的方向斜平行投影到接收平面上的具体过程如下：
在全局坐标系O-XYZ中，令R＝(Rx，Ry，Rz)为接收面上的一点，H＝(Hx，Hy，Hz)为在平行投影操作下在定日镜上对应R的一点，那么根据线面求交方程推导，有：
其中：
r＝(rx，ry，rz)是考察定日镜的反射光线方向；
T表示向量的转置操作；
N＝(Nx，Ny，Nz)是定日镜的法向；
o＝(ox，oy，oz)是定日镜的中心；
E是单位矩阵。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811578452.2A CN109670248B (zh) | 2018-12-21 | 2018-12-21 | 塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811578452.2A CN109670248B (zh) | 2018-12-21 | 2018-12-21 | 塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109670248A CN109670248A (zh) | 2019-04-23 |
| CN109670248B true CN109670248B (zh) | 2022-09-30 |
Family
ID=66145959
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811578452.2A Active CN109670248B (zh) | 2018-12-21 | 2018-12-21 | 塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109670248B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110209207B (zh) * | 2019-05-07 | 2022-05-24 | 中国神华能源股份有限公司 | 确定定日镜的损失面积的方法和装置及机器可读存储介质 |
| CN110414058B (zh) * | 2019-06-28 | 2020-12-01 | 浙江大学 | 一种塔式太阳能热电站中卷积核函数的拟合方法 |
| CN110414059B (zh) * | 2019-06-28 | 2020-12-01 | 浙江大学 | 塔式太阳能热电站中平面定日镜的辐射能密度模拟方法 |
| CN110414057B (zh) * | 2019-06-28 | 2020-12-11 | 浙江大学 | 塔式太阳能热电站中聚焦型定日镜的辐射能密度模拟方法 |
| CN110705077B (zh) * | 2019-09-25 | 2021-06-08 | 河海大学 | 一种塔式太阳能吸热器聚焦光斑能流密度分布的计算方法 |
| CN111625957B (zh) * | 2020-05-25 | 2022-06-10 | 浙江大学 | 塔式太阳能镜场接收器平面的辐射能密度分布模拟方法 |
| CN115082504B (zh) * | 2022-08-23 | 2022-11-11 | 菏泽学院 | 一种太阳能光伏板光斑识别方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530518A (zh) * | 2013-10-14 | 2014-01-22 | 中国科学院电工研究所 | 太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法 |
| CN106056670A (zh) * | 2016-05-30 | 2016-10-26 | 浙江大学 | 塔式太阳能热发电系统中剔除遮挡的辐射能密度模拟方法 |
| CN106096085A (zh) * | 2016-05-30 | 2016-11-09 | 浙江大学 | 一种塔式太阳能热发电系统的辐射能密度模拟方法 |
| WO2017131544A1 (en) * | 2016-01-25 | 2017-08-03 | Universidade De Évora | Compact linear fresnel reflective solar concentrator designed for direct molten salt operation as heat transfer fluid in evacuated tubes |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120304981A1 (en) * | 2011-06-01 | 2012-12-06 | Rolf Miles Olsen | Dynamic distributed tower receiver system for collecting, aiming and receiving solar radiation |
-
2018
- 2018-12-21 CN CN201811578452.2A patent/CN109670248B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530518A (zh) * | 2013-10-14 | 2014-01-22 | 中国科学院电工研究所 | 太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法 |
| WO2017131544A1 (en) * | 2016-01-25 | 2017-08-03 | Universidade De Évora | Compact linear fresnel reflective solar concentrator designed for direct molten salt operation as heat transfer fluid in evacuated tubes |
| CN106056670A (zh) * | 2016-05-30 | 2016-10-26 | 浙江大学 | 塔式太阳能热发电系统中剔除遮挡的辐射能密度模拟方法 |
| CN106096085A (zh) * | 2016-05-30 | 2016-11-09 | 浙江大学 | 一种塔式太阳能热发电系统的辐射能密度模拟方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于整数规划的塔式太阳能镜场能量调度方法;刘紫军等;《信息与控制》;20161015(第05期);第615-620 * |
| 塔式太阳能电站定日镜场的建模与仿真;余强等;《中国电机工程学报》;20120815(第23期);第90-98 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109670248A (zh) | 2019-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109670248B (zh) | 塔式太阳能热发电中光斑的辐射能密度分布的模拟方法 | |
| Sánchez-González et al. | Solar flux distribution on central receivers: A projection method from analytic function | |
| Jafrancesco et al. | Optical simulation of a central receiver system: Comparison of different software tools | |
| Wagner et al. | SolarPILOT: A power tower solar field layout and characterization tool | |
| CN109697315B (zh) | 辐射能光斑解析模型参数的优化方法 | |
| Yao et al. | Heliostat field layout methodology in central receiver systems based on efficiency-related distribution | |
| Oh et al. | A new algorithm using a pyramid dataset for calculating shadowing in solar potential mapping | |
| He et al. | An analytical flux density distribution model with a closed-form expression for a flat heliostat | |
| Gebreiter et al. | sbpRAY–A fast and versatile tool for the simulation of large scale CSP plants | |
| Lee et al. | Numerical study on optical performances of the first central-receiver solar thermal power plant in Korea | |
| Huang et al. | Development of an analytical method and its quick algorithm to calculate the solar energy collected by a heliostat field in a year | |
| de Almeida Rocha et al. | A pixel counting technique for sun patch assessment within building enclosures | |
| Huang et al. | Development of a new flux density function for a focusing heliostat | |
| He et al. | An improved flux density distribution model for a flat heliostat (iHFLCAL) compared with HFLCAL | |
| Lin et al. | GPU-based Monte Carlo ray tracing simulation considering refraction for central receiver system | |
| Wang et al. | A novel and efficient method for calculating beam shadows on exterior surfaces of buildings in dense urban contexts | |
| JP5137039B2 (ja) | 熱放射エネルギーのシミュレーション装置、及び、方法 | |
| Huang et al. | Gauss–Legendre integration of an analytical function to calculate the optical efficiency of a heliostat | |
| CN110414057B (zh) | 塔式太阳能热电站中聚焦型定日镜的辐射能密度模拟方法 | |
| CN111625957B (zh) | 塔式太阳能镜场接收器平面的辐射能密度分布模拟方法 | |
| Hu et al. | A comprehensive mathematical approach and optimization principle for solar flux distribution and optical efficiency in a solar tower | |
| Arrif et al. | Shadowing and blocking factors in heliostats: Comparison between parallel and oblique projections | |
| Richter et al. | Computationally fast analytical ray-tracer for central receiver systems | |
| He et al. | Analytical radiative flux model via convolution integral and image plane mapping | |
| Martínez-Hernández et al. | Advanced surface reconstruction method for solar reflective concentrators by flux mapping |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |