《如何计算太阳的方位角》由会员分享，可在线阅读，更多相关《如何计算太阳的方位角（4页珍藏版）》请在金锄头文库上搜索。
1、太阳辐射怎么计算在太阳能利用工作中，太阳辐射计算十分重要。为了帮助读者掌握太阳辐射计算方法，我们请长期从事太阳辐射研究工作的中国气象科学研究院王炳忠研究员编写了太阳辐射计算讲座 ，供大家学习、参考。1日地距离地球绕太阳公转的轨道是椭圆形的，太阳位于椭圆两焦点中的一个。发自太阳到达地球表面的辐射能量与日地间距离的平方成反比，因此，一个准确的日地距离值 R 就变得十分重要了。日地平均距离 R0，又称天文单位，1 天文单位=1.496108km 或者，更准确地讲等于 149597890500km。日地距离的最小值（或称近日点）为0.983 天文单位，其日期大约在 1 月 3 日；而其最大值（或称远日
2、点）为 1.017 天文单位，日期大约在 7 月 4 日。地球处于日地平均距离的日期为 4 月 4 日和 10 月 5 日。 由于日地距离对于任何一年的任何一天都是精确已知的，所以这个距离可用一个数学表达式表述。为了避免日地距离用具体长度计量单位表示过于冗长，一般均以其与日地平均距离比值的平方表示，即 ER=（rr0） 2，也有的表达式用的是其倒数，即 r0r，这并无实质区别，只是在使用时，需要注意，不可混淆。我们得到的数学表达式为 ER=1.0004230.032359sin0.000086sin2 0.008349cos0.000115cos2 （1）式中 称日角，即=2t365.2422
3、 （2） 这里 t 又由两部分组成，即 t=NN0 （3） 式中 N 为积日，所谓积日，就是日期在年内的顺序号，例如，1 月 1 日其积日为 1，平年 12 月 31 日的积日为 365，闰年则为 366，等等。N0=79.67640.2422 （年份1985）INT（年份1985 ）4 （4） 2太阳赤纬角 地球绕太阳公转的轨道平面称黄道面，而地球的自转轴称极轴。极轴与黄道面不是垂直相交，而是呈 66.5角，并且这个角度在公转中始终维持不变。正是由于这一原因形成了每日中午时刻太阳高度的不同，以及随之而来的四季的变迁。太阳高度的变化可以从图1 中形象地看到。图中日地中心的连线与赤道面间的夹角每
4、天（实际上是每一瞬间）均处在变化之中，这个角度称为太阳赤纬角。它在春分和秋分时刻等于零，而在夏至和冬至时刻有极值，分别为正负 23.442。 图 1 地球绕太阳运行轨迹 由于太阳赤纬角在周年运动中任何时刻的具体值都是严格已知的，所以它（ED）也可以用与式（1）相类似的表达式表述，即：ED=0.372323.2567sin0.1149sin20.1712sin30.758cos0.3656cos20.0201cos3 （5） 式中 的含义与式（1 ）中的相同。3时差 真正的太阳在黄道上的运动不是匀速的，而是时快时慢，因此，真太阳日的长短也就各不相同。但人们的实际生活需要一种均匀不变的时间单位，这
5、就需要寻找一个假想的太阳，它以均匀的速度在运行。这个假想的太阳就称为平太阳，其周日的持续时间称平太阳日，由此而来的小时称为平太阳时。 平太阳时 S 是基本均匀的时间计量系统，与人们的生活息息相关。由于平太阳是假想的，因而无法实际观测它，但它可以间接地从真太阳时 S求得，反之，也可以由平太阳时来求真太阳时。为此，需要一个差值来表达二者的关系，这个差值就是时差，以 Et 表示，即 S=S Et （6）由于真太阳的周年视运动是不均匀的，因此，时差也随时都在变化着，但与地点无关，一年当中有 4 次为零，并有 4 次达到极大。时差也可以以式（1）相似的表达式表示： Et=0.00281.9857sin
6、9.9059sin27.0924cos0.6882cos2（7） 上面，我们给出了 3 个计算式，从形式上讲，它们与一般书籍中给出的并无不同。我们之所以又重新研究它，是因为以往的公式存在以下的通病：对平年和闰年不加区分，一方面，这对闰年就不好处理，另一方面，闰年的影响有累计效应，会逐步增长；即使是从当年天文年历查到的数值，也是格林尼治经度处 0 点时刻的数值，而我们所需要的数值，会因所在地点的地理经度以及具体时刻与表值有异而不同。具体地讲，一般要进行如下 3 项订正：2 如何计算太阳的方位角?（1）年度订正：除非我们只用当年的天文年历值，此外均需使用此项订正，引入此项订正的原因就是一回归年的实
7、际长度不是 365 日，而是 365.2422 日，但日历上只有整日，不可能有小数日。假定我们选用的是 1981 年的表值，1982 年再用时，就要加上0.2 （0.2422 ）日的订正了。这个订正到了 1983 年为0.51（0.4844）日，1984年为0.7（0.7266）日，但此年为闰年，多了 1 日，实际订正应为0.7 1=0.3（0.2734 ）日， 1985 年为 0.0（0.0312）日，等等，余类推。（2）经度订正：即使我们查阅的是当年的天文年历，也需此项订正。在我国的地理经度范围内，各地的订正值是90E 0.2 日 90E2thenC=31.890G=INT（30.6*YC
8、0.5）R100 L=（JDJF60）15110 H=S8 F 60 120 N=G（HL ）24130 T=（N N0）K140ER=1.0004230.032359sinT0.000086sin2T0.008349cosT0.000115cos2T150ED=0.372323.2567sinT 0.1149sin2T0.1712sin3 T0.758cos T0.3656cos2T0.0201cos3T160Et=0.00281.9857sinT9.9059sin2T7.0924cosT0.6882cos2T170print“Er=”；Er；“Ed=”；Ed， “Et=”；Et 180in
9、put“ 是否仍要计算 yn？”，W0 190ifW=“Y”orW=“y”then10else200200end 程序中 5090 各句的目的在于计算当天的积日，100 句是经度订正，110 句是时刻订正，130 句包含了年度订正的内容。 在太阳能利用中，最常见的是要求计算太阳高度和太阳方位。太阳高度（h）的计算公式为sinh=sinsin coscoscos （8）式中， 就是太阳赤纬角，即式（5）中的 Ed， 为当地的地理纬度， 为当时的太阳时角。 值不难获得，且一旦确定，不会改变。 值的计算可以从前述程序中得到。唯一需要说明的是太阳时角 的计算。其计算式为 （9）这里时 S 和分 F 的
10、符号均加上了下标，表示是真太阳时，为了从北京时求出真太阳时，需要两个步骤：首先，将北京时换成地方时 Sd： SdS+（地理经度-120）*4/60 （10）式中，120 是北京时的标准经度，乘 4 是将角度转化成时间，即每度相当于 4 分钟，除 60 是将分钟化成小时。 其次，进行时差订正，即S=SdEt60 （11）这里应该指出的是，时角是以太阳正午时刻为 0 点的，顺时针方向（下午）为正，反之为负。太阳方位角的计算式为cosA=（sinh sinsin ）coshcos （12）由此可求出二个 A 值，第一个 A 值是午后的太阳方位， 当 cosA0 时 90A180当 cosA0 时 0
11、A90第 2 个 A 值为午前的太阳方位，取 360A 。 实例：计算东经 110北回归线上 1999 年 6 月 23 日北京时 1242 的太阳高度角及当日的日落时的方位角。 计算：将 JD=110，JF=0，NF=1999 ，Y=6，R=23，S=12，F=42 ，各参数输入运行中的程序；屏幕上立即显示：Er=1.0330，Ed=23.438 ，Et= 1.84 将北京时 1242 换算成东经 110的地方时，利用式（10） ，可得 Sd=1202加当日时差 Et2，得此时当地的 S=1200，将其代入式（ 9）得 =0，北回归线处 =23.442 最后根据式（8）求得 h=89.966读者可能产生疑问，为何在北回归线上，夏至日的中午时刻的太阳高度不等于 90，大家不妨变换 NF 的输入值，看一看结果不仅都不等于 90，且各年之间还略有差异。之所以会如此，是因为夏至不仅有日期，还有时刻，很难遇到夏至时刻在正午是 12 时的。在计算日落时的方位角时，由于此时 h=0，所以式（12）的形式有所变化：cosA=sin cos （13）将已知参数代入，得 cosA=0.3977 依照判据 90A180，故 A=113.44
