太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法
【专利摘要】一种太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法，对于单侧布置定日镜场对应的吸热器，将吸热器的采光面做矩形网格划分，形成矩形平面网格或圆柱面网格，同时确定采光面的“田”字形网格边界；采用逆向光线追迹法，计算定日镜场在吸热器采光面“田”字形网格边界上的聚光能流密度值；把吸热器采光面网格展开成X-Y平面网格，由“田”字形网格边界上的聚光能流密度的计算值，插值重建吸热器采光面所有网格点上的能流密度值。对于周向布置定日镜场，可对每个圆柱面采光面子区域，先计算其“田”字形网格边界上的能流密度值，再通过插值重建出对应采光面子区域上的聚光能流密度分布。
【专利说明】太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法
【技术领域】
[0001]本发明属于太阳能热利用领域，涉及一种计算太阳能塔式电站吸热器采光面上聚光能流密度分布的方法。
【背景技术】
[0002]太阳能塔式热发电站中，吸热器一般放置在中央吸热塔的顶端，太阳塔周围的地面上是众多的定日镜，形成聚光场。这些定日镜均能绕双轴旋转，自动跟踪太阳，把入射太阳光反射向吸热器的米光面。定日镜聚光场相对于吸热塔，一般有北向单侧布置和周向布置两种模式。对于北向单侧布置的定日镜场，塔顶吸热器的一侧是采光面，向着聚光场，且采光面多为平面或圆柱面(圆柱面的一部分，即槽形圆柱面)。对于周向布置的定日镜场，塔顶吸热器的采光面多为圆柱面(完整的圆柱面，即筒形圆柱面)。
[0003]吸热器采光面是太阳能塔式热发电站聚光能量流的一个天然界面，既是聚光场能量流的出口，也是吸热器能量流的入口，因而计算和分析吸热器采光面上的聚光能流密度分布对评价聚光场的光学性能以及分析吸热器的吸热效率都很重要。
[0004]计算的吸热器采光面上的能流密度分布，一般以能流密度值数表的离散形式表示，即以某种网格间距对采光面作矩形网格划分，并计算出采光面上所有网格点的聚光能流密度值。平面采光面对应平面矩形网格划分，圆柱面采光面对应柱面矩形网格划分。
[0005]计算定日镜场在吸热器采光面上的聚光能流密度分布主要有两种方式，一种是解析计算，另一种是采用光线追迹法。不管采用哪种计算方式，都是先计算聚光场中单台定曰镜在吸热器采光面上的能流密度分布，然后把所有定日镜的聚光能流密度分布叠加起来，便得到整个聚光场在吸热器采光面上的聚光能流密度分布。解析的方式是在一些简化假设条件下计算聚光能流密度分布，精度不高，适用范围有限，计算的灵活性差。随着计算机性能的快速提高，以光线追迹法计算聚光能流密度分布越来越受青睐。光线追迹法的好处是计算灵活、方便，几乎适用于所有的计算场合，且计算精度也较高，但计算量较大。光线追迹法又分为蒙特卡洛光线追迹法、基于等效太阳光锥的正向光线追迹法和基于等效太阳光锥的逆向光线追迹法。
[0006]太阳光不是严格的平行光。太阳轮廓相对镜面上一点形成半张角是16'的圆椎面，从而形成太阳光锥。太阳光锥的空间能量分布，对应于太阳圆盘的亮度分布。表示太阳亮度分布的模型比较多，有圆形均匀分布，圆形高斯分布，Kuiper亮度分布模型，LBL实验数据模型，Buie亮度分布模型，等等。由于镜面光学误差的存在，入射的太阳光锥，经镜面反射后形成扩散的太阳光锥。为了在光学计算过程中不用考虑镜面误差的光学扩散，可以把扩散了的太阳光锥作为等效太阳光锥，即等效太阳光锥的空间能量分布是入射太阳光锥的能量分布与镜面误差分布的复合。基于等效太阳光锥的概念，可以用一束光线(由多条光线组成)代替等效太阳光锥，以正向(从光源到采光面)或逆向(从采光面到光源)光线追迹的方式计算聚光能流密度分布。[0007]无论正向还是逆向光线追迹，都需要把反射镜面均匀地划分为网格，把镜面网格单元的中点作为该镜面单元的代表点，把镜面网格单元看作小平面，镜面单元的法向即为镜面网格单元中点处的法向。定日镜在吸热器采光面上的聚光能流密度等于定日镜所有镜面单元聚光能流密度的累加。
[0008]总之，要计算吸热器采光面上的能流密度分布，需要计算采光面上所有网格点处的聚光能流密度值，需要的计算量很大。而且商业化太阳能塔式热发电站总的反射镜面积在几万平米以上，即便是计算采光面上一点的能流密度值，计算量也已不小。因此，需要发展一种快速的吸热器采光面能流密度分布计算方法。
【发明内容】
[0009]本发明的目的是克服上述现有聚光能流密度分布计算方法的不足，针对太阳能塔式热发电站吸热器采光面上的聚光能流密度分布，提供一种简单、快速的计算方法。
[0010]本发明针对太阳能塔式热发电站的定点聚光的特点，即对于单侧布置的聚光场，大量的定日镜同时向吸热器采光面中心点聚光，在吸热器采光面上的聚光能流分布是中央强边缘弱，呈现接近圆形的椭圆形高斯分布特征。本发明利用二元高斯分布密度函数的性质，只需通过基于等效光锥的逆向光线追迹法，计算吸热器采光面“田”字形网格边界上的聚光能流密度值，再插值重建整个吸热器采光面上的聚光能流密度分布。采光面“田”字形网格边界，是指由米光面网格点阵列的最上、中间、最下3行网格点和最左、中间、最右3列网格点组成，是采光面网格点阵列中的一部分典型网格点。
[0011]记二元高斯分布密度函数为f(x，y)，则f(x，y)在矩形区域上的插值重建的数学原理如下:
[0012]设矩形区域田字形网格边界的最左、中间、最右3条纵向网格线的χ值分别是X0> XmidI和xm，最上、中间、最下3条横向网格线的I值分别是yN、ymid2和10，f (x, y)限制在田字形网格上的6条函数边界分别是f (χ。, y)，yQ≤y≤yN，f (xmid1.y).y0≤y≤yN，f (χΜ, y)，y0 ≤ y ≤ yN，f (χ, yN)，χ0 ≤ χ ≤ χΜ，f (χ, ymid2)，χ0 ≤ χ ≤ χΜ 和 f (χ, y0)，χ0 ≤ χ ≤ χΜ。令S (χ, y) =In (f (X，y)), S (X，y)限制在田字形网格上的6条函数边界分别是S (x0, y)，y0 ≤ y ≤ yN，S (xmidl, y)，y0 ≤ y ≤ yN，S (xM, y)，y0 ≤ y ≤ yN，S (x, yN)，χ0 ≤ χ ≤ χΜ,S (χ, ymid2)，χ。≤ χ ≤ χΜ 和 S (x, yQ)，xQ ≤ χ ≤ χΜ,则有:
【权利要求】
1.一种太阳能塔式发电吸热器采光面聚光能流密度分布的计算方法，其特征在于，所述的计算方法为:对于与单侧布置定日镜场对应的吸热器，将吸热器的采光面以一定网格间距做矩形网格划分，形成矩形平面网格或槽形圆柱面网格，同时确定采光面的“田”字形网格边界；所述的“田”字形网格边界由采光面网格点阵列中的最上、中间、最下3行网格点和最左、中间、最右3列网格点组成；再采用逆向光线追迹法，只计算定日镜场在吸热器采光面“田”字形网格边界上的聚光能流密度值；然后把吸热器采光面网格展开成X-Y平面网格，由计算得到的“田”字形网格边界上的聚光能流密度值，插值计算吸热器采光面上其余网格点上的能流密度值；对于计算周向布置定日镜场在吸热器采光面上的聚光能流密度分布，分解为计算多个单侧布置定日镜场在各自的采光面区域上聚光能流密度分布，即对每个圆柱面采光面子区域，先计算其“田”字形网格边界上的能流密度值，再插值重建出对应采光面子区域上的聚光能流密度分布。
2.根据权利要求1所述的太阳能塔式热发电吸热器采光面聚光能流密度分布的计算方法，其特征在于，插值重建所述的吸热器采光面上聚光能流密度分布的步骤是:首先根据吸热器采光面“田”字形网格边界上的聚光能流密度计算值，通过计算式(1)-(5)插值得到吸热器采光面区域内各网格点上的聚光能流密度值；然后根据显示或分析的需要进一步插值细化吸热器采光面上能流密度分布:
【文档编号】G06F19/00GK103530518SQ201310479508
【公开日】2014年1月22日 申请日期:2013年10月14日 优先权日:2013年10月14日
【发明者】郭明焕, 王志峰 申请人:中国科学院电工研究所