光阵列面倾角倾斜面辐射量的计算.doc

光整列面倾角及辐射量的计算直接辅射 不同倾斜面上的直射辐照度可利用下式求出：S（，）= Smcos 式中是太阳光线对倾斜面的入射角，可由下式得出： cos=cosSinh+Sincoshcos(-) 式中是倾斜面与水平面间的夹角，h是太阳高度角，是太阳的方位角，是倾斜面的方位角，方位角从正南算起，向西为正，向东为负。对于水平面来说，由于=0，所以cos=Sinh，因此：S（0，0）= SmSinh 设KS=S（，）/S（0，0），将前面的公式代入，则有：KS=cos/Sinh=cos+Sincos(-) /tanhKS称为换算系数。有了KS值，根据水平面上的辐射值很容易求出倾斜面的辐射值。对于不同时段的曝辐射量，也是如此。只时求算KS时，、h等值要代入相应时段的平均值。当计算较长时段内的曝辐射量时，如日总量，使用换算系数也很方便，只是这时的KS值应从实测值中得出，而不能用上述几何关系计算出来。对于实用来说，用月平均日总量的KS值最方便，它比个别日子的KS值对云量和透明状况的依赖性更少。其他影响KS的因子是地点的纬度、倾斜面的朝向和月份等。表13给出了不同纬度三种倾斜角度月平均日总量的KS值。散射辐射 计算不同朝向倾斜面上的散射辐照度，困难要大得多。通常的解决办法是假定辐射是各向同性的，即呈均匀分布。这样，散射辐照度Ed和反射辐照度Er可按下列公式计算。 Ed（，）= Ed(1+ Cos)/2 Er（，）= Er(1- Cos)/2式中Ed和Er是水面上的散射和反射辐照度。不过，用下式根据水平面上的散射辐照度计算倾斜面上的散射辐照度，要比利用各向同性的假设更准确此。 Ed（，）+ Er（，）=K(Ed+ Er)Ed换算系数K（Ed+Er）是在各种太阳高度角和方位角下，用总辐射表对各种倾斜表面上的散射辐照度和反射辐照度进行实测的结果确定的。表14给出了不同混浊情况下不同的K（Ed+Er）值。总辅射 在各向同性的前提下，倾斜面上的总辐射可用下式算出：Eg（，）=KsSm+ Ed(1+ Cos)/2+ Er(1- Cos)/2不过，对于大多数用户来说，通过换算系数Kg直接从水平面的总辐射求出Eg(，)更方便，即Eg（，）=KgEg表15 是国外发表的在一些情况下总辐射月平均日总量的Kg值。 表14 晴天落日垂直面上的散射辐照度相对值(下垫面草地,反射系数约20%)(-)度h(度)2510305070B=0.0505.83.62.62.22.01.8403.52.21.71.81.91.8901.71.10.91.31.61.71301.00.70.71.11.41.61800.70.40.50.91.21.4B=0.1004.43.12.62.21.71.5402.72.22.01.91.51.4901.61.41.21.21.21.31301.11.00.80.91.01.21800.70.70.60.70.91.1B=0.1503.93.32.61.91.51.2402.62.42.01.61.31.1901.41.41.21.11.01.01300.90.80.70.80.90.91800.60.60.50.60.70.8B=0.2004.53.32.41.61.31.0403.12.52.01.41.11.0901.61.41.10.90.90.81300.60.60.50.70.70.71800.30.40.40.50.60.7计算实例：例1已知Eg=845W/M2，Ed=845W/M2，反射系数=0.20，计算北纬45某地6月16日真太阳时11时朝南30倾斜面上的直射、散射、反射和总辐射照度。水平面上的直射辐照度Sm=Eg- Ed=705W/ M2太阳高度角h和方位角可由下式得出Sinh=Sin Sin+CosCosCos =0.7070.396+0.7070.9180.966 =0.907 h=65.1Sin=CosSin/Cosh=-0.238/0.421=-0.566因为是午前11时，太阳偏东=-34.5 -=-34.5-0=-34.5 Ks= cos+Sincos(-) /tanh =0.866+0.5000.824/2.154 =1.06倾斜面上的直射辐照度 S（，）=Ks Sm =1.06705W/M2 =747W/M2 倾斜面上的散射辐照度 Ed（，）= Ed(1+ Cos)/2=140W/M20.933 =131W/ M2 倾斜面上的反射辐照度 Er（，）= Er(1- Cos)/2 Er= Eg=0.2845=169 W/ M2 ErT（，）=1690.067=11 W/ M2倾斜面上的总辐射辐照度 Eg（，）=S（，）+ Ed（，）+ Er（，）=747+131+11=889W/M2表15 晴天不同纬度各种倾斜面上总辐射日总量的相对值Kg(月平均)倾斜角度纬度北半球,朝南的4567830100.950.951.02501.011.081.23601.061.201.4560300.700.730.80400.750.780.92500.800.901.10600.901.081.259000.200.190.190.210.21100.250.200.200.220.25200.300.220.220.240.30300.400.250.250.280.40400.600.350.350.360.50500.950.500.500.530.72600.620.620.680.96南半球,朝北的10111212例2.条件仍如例1，求向东倾斜30平面上的总辐射辐照度。因为倾斜平面朝东，=-90-=-34.5-（-90）=55.5Cos(-)=0.566这时，KS=0.866+0.50.566/2.145=0.9981倾斜面上的直射辐照度Sm(30,-90)=1705=705 W/M2由于各向同性的假定，所以散射和反射辐照度没有什么变化，于是倾斜面上的总辐射辐照度Eg（30，-90）=705+131+11=847 W/M2例3已知多年实测的平均日总量为Hg=31MJ/ M2dHd=5.3MJ/ M2d地表反射系数=0.23，求北纬41某地6月某一晴天朝南倾斜30平面上的总辐射和直接辐射的多年平均日总量。具体的计算步骤如下：(1)水平面上的直射辐照度Hs= Hg- Hd=31-5.3=25.7 MJ/ M2d(2)从表13中内插出北纬41朝南倾斜30平面上的换算系数Ks=0.93(3)倾斜面上的直射曝辐射量Hs（30，0）=0.9325.7=23.9 MJ/ M2d(4)从表15中内插求出北纬41朝南倾斜30平面上的总辐射换算系数 Kg=0.956(5)倾斜面上的总辐射曝辐射量Hg（30，0）=0.95631=29.6 MJ/ M2d当然，也可按各向同性的假设去求算，即（1）确定散射和反射的换算系数散射（1+Cos30）/2=1.866/2=0.933反射（1-Cos30）/2=0.134/2=0.067（2）根据反射系数计算反射辐射量Hr=Hg=0.2331=7.1 MJ/ M2d（2）计算倾斜面上的总辐射曝辐射量Hg（30，0）=23.9+0.9335.3+0.0677.1=29.3MJ/ M2d用两种方法所行结果是相当接近的。例4已知南纬2231某地日射站2月份中午的多年月平均总辐射辐照度为1117W/M2，散射辐照度为70W/M2，=0.22，求算当地向南、向北倾斜60平面上的中午月平均总辐射辐照度。具体步骤如下：（1）确定SmSm= Eg- Ed=1117-70=1047W/M2(2)确定h，即2月15日中午的hSinh=Sin Sin+CosCosCos =0.985h=80(3)确定方位南半球2月份太阳位于太空的北半部=180向南倾斜面的方位=0向北倾斜面的方位=180向南倾斜面方位与太阳方位夹角-=180向北倾斜面方位与太阳方位夹角-=0(4)确定直射换算系数 Ks=Cos+ SinCos（-）/tanh= Cos60+ Sin60Cos（-）/tan80向南的Ks=0.5+0.866(-1)/5.6713=0.345向北的Ks=0.5+0.8661/5.6713=0.65(5)确定倾斜面上的直射辐照度向南的S(60,0)=0.3451047=361W/M2向北的S(60,180)=0.6551047=685W/M2(6)确定倾斜面上的散射辐照度,向南与向北的无差异,均为 Ed（60，0或180）= Ed(1+ Cos)/2=70W/M20.75 =52W/ M2 (7)确定倾斜面上的反射辐照度Er（60，0或180）=Er(1- Cos)/2=0.2211170.25=61 W/ M2(8)确定倾斜面上的总辐射辐照度向南的Eg（60，0）=361+52+61=474 W/M2向北的Eg（60，180）=685+52+61=798 W/M2在太阳能利用工作中，太阳辐射计算十分重要。为了帮助读者掌握太阳辐射计算方法，我们请长期从事太阳辐射研究工作的中国气象科学研究院王炳忠研究员编写了太阳辐射计算讲座，供大家学习、参考。 编者 在太阳能利用中，为了多获得能量，一般将太阳能装置朝南倾斜放置。若想了解倾斜状态下，地表的太阳能量，最直接、也是最准确的方法是利用仪器进行直接测量。在无法进行直接测量的情况下，如果当地的气象站有日射观测项目，则可利用其测得的水平面上的数据进行换算。(注：气象站测量的都是水平面上的数据) 具体的换算方法如下(文中符号的物理意义与前几讲一样)： 参照图1并结合第四讲中介绍的内容可知： En=Eocosz(1) En,=Eo,coso(2) 因此， rb=En,/En=coso/cosz(3) 同理，对于曝辐量来说，可将上一讲中的计算Ho和H的公式代入，经过简化可得出： (4) 以上各式都是针对地外，即没有大气情况下得到的，当我们感兴趣的是实际地表(即存在大气)时的情况，又当如何呢?我们说，情况是相同的，因为rb和Rb是个比值，影响它们的只是入射角度的不同，至于经过大气衰减，还是未经衰减，对于比值来说，可以相抵消，而实际上没有什么影响。 以上的讨论还只限于直射部分。对于斜面上的辐射来说，除了直射外，还有天空中的散射以及地面反射部分。为了简化问题，我们只能近似地假定散射和反射特性都是完美的，即各向同性。实际上，问题当然不会这么简单，比如说天空中影响散射最明显的云量，其分布完全是随机的，当然不会是均匀的，但不做如此假设，无法处理这么复杂的问题。 在各向同性的假定下，对于反射辐射，则有：对于散射辐射，则有 (7) (8) 式(5)和(6)中的为地表的反射比。 如果在当地的气象站有直射和散射的数据，就可直接代入上式进行计算，但是可提供直射的气象站较少。如果有总辐射和散射数据，则可以EEd或HHd作为水平面上的直射量。这样一来，斜面上的总辐射便可以写成 E=(EEd)rbErEd(9) H=(HHd)RbHrHd(10) 计算实例：设位于北纬4530的某地，6月4日全天的总辐射日曝辐量为29.71MJ/m2，散射曝辐量为4.73MJ/m2，地表反射比=0.20，求朝南倾斜50表面上的日曝辐量? 解：根据式(4)，计算得出 Rb=0.77 直射部分： Hb=(HHd).Rb=(29.714.73)0.77=19.23MJ/m2 反射部分： 29.710.2(1cos50)=1.06MJ/m2 散射部分： (1cos50)=3.89MJ/m2 图1 斜面上的总辐射：H=HbHrHd=24.18MJ/m2 本讲座至此共进行了5讲，其实太阳辐射需要计算的项目还有很多，如光谱量、大气浑浊度、太阳能资源等等。考虑到这些内容涉及专业知识较多，一方面枯躁乏味，另一方面偏离科普宗旨，同时应用面也较窄，所以本讲座对这些内容不再介绍。如读者对某项内容感兴趣，可通过本刊编辑部与作者联系，作者当尽力协助解决。最后还想提醒几点在计算中应注意的事项： (1)由于日地距离订正系数存在着r20/r2与r2/r20两种形式，同时还有究竟是由日地平均距离处向具体日期的订正，还是相反，由具体某一日期向日地平均距离处订正的问题。虽然订正不外乎乘、除两种方式，但面临具体问题时，却容易搞混，从而导致错误。为了避免该情况的发生，请大家切不可死记某一公式，而应牢记如下方法：从日地平均距离处向某一具体日期订正，并规定此日期为1月1日。如果乘(或除)以该系数后，其值大于原值，则表明做对了，因为1月1日地球处于近日点附近，其值理应大于日地平均距离处的，以此可以类推其它日期的做法，结果亦应如此；反这，如果数值小于原值，则应改乘(除)为除(乘)。 (2)地方时、世界时、北京时在运用中也是容易混淆的。例如在计算中使用的太阳赤纬都是以天文年历为准的，而天文年历所给出的参数都是世界时0时的值，但时角又是以地方时为