在太阳能利用工作中

在太阳能利用工作中,阳太射计算辐太阳辐射计算十分重要。为讲了帮助读者掌握太阳辐射计算方法,我们请长期从事太阳辐射研究工作的中国气象科学研究院王炳忠研究员编写了《太阳辐射计算讲座》供大家学习、。,参考座—编者—第一讲太阳能中天文参数的计算王炳忠地球绕太阳公转的轨道是椭圆形的,太阳位于椭圆两焦点中的一个。发自太阳到达地球表面的辐射能量与日地间距离的平方成反比,因此,一个准确的日地距离值R就变得十分重要了。日地平均距离R0,又称天文单位,81天文单位=11496×km10或者,更准确地讲等于149597890±500km。日地距离的最小值(或称近日点)为01983天文单位,其日期大约在1月3日;而其最大值(或称远日点)为11017天文单位,日期大约在7月4日。地球处于日地平均距离的日期为4月4日和10月5日。由于日地距离对于任何一年的任何一天都是精确已知的,所以这个距离可用一个数学表达式表述。为了避免日地距离用具体长度计量单位表示过于冗长,一般均以其与日地平均距离比值的平方表示,即ER=(r/r0)2,也有的表达式用的是其倒数,即r0/r,这并无实质区别,只是在使用时,需要注意,不可混淆。我们得到的数学表达式为θθER=11000423+01032359sin+01000086sin2θθ-01008349cos+01000115cos2(1)式中θ称日角,即θ=2πt/36512422(2)这里t又由两部分组成,即t=N-N0式中N为积日,所谓积日,就是日期在年内的顺序号,例如,1月1日其积日为1,平年12月31日的积日为365,闰年则为366,等等。811日地距离(3)地球绕太阳公转的轨道平面称黄道面,而地球的自转轴称极轴。极轴与黄道面不是垂直相交,而是呈6615°,并且这个角度在公转中始终维持不变。正是角由于这一原因形成了每日中午时刻太阳高度的不同,以及随之而来的四季的变迁。太阳高度的变化可以从图1中形象地看到。图中日地中心的连线与赤道面间的夹角每天(实际上是每一瞬间)均处在变化之中,这个角度称为太阳赤纬角。它在春分和秋分时刻等于零,而在夏至和冬至时刻有极值,分别为正负231442°。由于太阳赤纬角在周年运动中任何时刻的具体值都是严格已知的,所以它(ED)也可以用与式(1)相类似的表达式表述,即:θθED=013723+2312567sin+011149sin2θθθ-011712sin3-01758cos+013656cos2θ(5)+010201cos3式中θ的含义与式(1)中的相同。真正的太阳在黄道上的运动不是匀速的,而是时快时慢,因此,真太阳日的长短也就各不相同。但人们的实际生活需要一种均匀不变的时间单位,这就需要寻找一个假想的太阳,它以均匀的速度在运行。这个假想的太阳就称为平太阳,其周日的持续时间称平太阳日,由此而来的小时称为平太阳时。平太阳时S是基本均匀的时间计量系统,与人们的生活息息相关。由于平太阳是假想的,因而无31时差21太阳赤纬角N0=7916764+012422×(年份-1985)-INT〔(年份-1985)/4〕(4)图1地球绕太阳运行轨迹实际观测它,但它可以间接地从真太阳时S⊙求得,反之,也可以由平太阳时来求真太阳时。为此,需要一个差值来表达二者的关系,这个差值就是时差,以Et表示,即S⊙=S+Et(6)由于真太阳的周年视运动是不均匀的,因此,时差也随时都在变化着,但与地点无关,一年当中有4次为零,并有4次达到极大。时差也可以以式(1)相似的表达式表示:θθEt=010028-119857sin+919059sin2θθ(7)-710924cos-016882cos2上面,我们给出了3个计算式,从形式上讲,它们与一般书籍中给出的并无不同。我们之所以又重新研究它,是因为以往的公式存在以下的通病:①对平年和闰年不加区分,一方面,这对闰年就不好处理,另一方面,闰年的影响有累计效应,会逐步增长;②即使是从当年天文年历查到的数值,也是格林尼治经度处0点时刻的数值,而我们所需要的数值,会因所在地点的地理经度以及具体时刻与表值有异而不同。具体地讲,一般要进行如下3项订正:(1)年度订正:除非我们只用当年的天文年历值,此外均需使用此项订正,引入此项订正的原因就是一回归年的实际长度不是365日,而是36512422日,但日历上只有整日,不可能有小数日。假定我们选用的是1981年的表值,1982年再用时,就要加上012(-012422)日的订正了。这个订正到了1983年为-0151(-014844)日,1984年为-017(017266)日,但此年为闰年,多了1日,实际订正应为-017+1=013(012734)日,1985年为010(010312)日,等等,余类推。(2)经度订正:即使我们查阅的是当年的天文年的订正值是:002424484812时段336-6006-88-1010-13日+012+013+014+015时段1312-15361536-18001800-2024日+016+017+018由于我国普遍采用的是北京时,它与格林尼治的地方时相差8小时,故具体到我国情况:时段(北京时)200-424424-648648-912912-1136订正值(日)-012-0110011时段1136-14001400-16241624-18481848-2112订正值012013014015前面3个计算式,项数多计算麻烦,后面多项订正,更显繁琐。为了方便实际应用,特编制如下仅含20句的BASIC语言程序,供使用:10input经度,经分和年份”D,JF,NF“,J20A=NF/4:K=23311415926#/3651242230N0=7916764+0124223(NF-1985)-INT((NF-1985)/4)40input月,日,时,分(按北京时)”Y,R,S,F“,50B=A-INT(A)6070if80if90C=3218100L=(JD+JF/60)/15H=S-8+F/60=(N-N0)/KN=G+(H-L)/24110120130140式(1)150式(5)160式(7)G=INT(30163Y-C+015)+RB=0andY2thenC=3118Y≤thenC=30162历,也需此项订正。在我国的地理经度范围内,各地的订正值是≤°90E-012日90°～128°EE-013日≥°128E-014日(3)时刻订正:要求同前一项。即使在格林尼治当地,不同时刻也需加以订正。各时段9170printEr=”Er“Ed=”Ed,“Et=”Et“;;;;180input是否仍要计算y/n?”“,W190ifW200end=Y”rW“o=y”then10else200“顺时针方向(下午)为正,反之为负。太阳方位角的计算式为δcosA=(sinh⊙sin-sin)/cosh⊙cosφφ(12)程序中50-90各句的目的在于计算当天的积日,100句是经度订正,110句是时刻订正,130句包含3年度订正的内容。在太阳能利用中,最常见的是要求计算太阳高度和太阳方位。太阳高度(h⊙)的计算公式为δφδφτ(8)sinh⊙=sinsin+coscoscosδ就是太阳赤纬角,即式(5)中的Ed,为当地的式中,φ地理纬度,τ为当时的太阳时角。值不难获得,且一φδ旦确定,不会改变。值的计算可以从前述程序中得到。唯一需要说明的是太阳时角的计算。其计算式为τ=(S⊙+F⊙60这里时S和分F的符号均加上了⊙下标,表示是真太阳时,为了从北京时求出真太阳时,需要两个步骤:首先,将北京时换成地方时Sd:JF)]×60Sd=S+{F-[120°(JD+4}/60式中,120°是北京时的标准经度,乘4是将角度转化成时间,即每度相当于4分钟,除60是将分钟化成小时。其次,进行时差订正,即S⊙=Sd+Et/60这里应该指出的是,时角是以太阳正午时刻为0点的,15.-12)×(9)(10)(11)10由此可求出二个A值,第一个A值是午后的太阳方位,当cosA≤时090°A≤°≤180当cosA≥时00≤A≤°90第2个A值为午前的太阳方位,取360°A。实例:计算东经110°北回归线上1999年6月23日北京时12∶42的太阳高度角及当日的日落时的方位角。计算:将JD=110,JF=0,NF=1999,Y=6,R=23,S=12,F=42,各参数输入运行中的程序;屏幕上立即显示:Er=110330,Ed=231438,Et=-1184将北京时12∶换算成东经110°42的地方时,利用(10),可得Sd=12∶式02加当日时差Et≈-2,得此时当地的S⊙=12∶,将00其代入式(9)得τ=0°北回归线处φ=231442°,最后根据式(8)求得h⊙=891966°读者可能产生疑问,为何在北回归线上,夏至日的中午时刻的太阳高度不等于90°大家不妨变换NF的,输入值,看一看结果不仅都不等于90°且各年之间还,略有差异。之所以会如此,是因为夏至不仅有日期,还有时刻,很难遇到夏至时刻在正午是12时的。在计算日落时的方位角时,由于此时h⊙=0,所以式(12)的形式有所变化:δφ(13)cosA=-sin/cos将已知参数代入,得cosA=-013977依照判据90°A≤°故A=113.44°≤180,