留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

太阳辐射计的衍射效应修正

刘国栋 方伟 宋宝奇 叶新 王凯

刘国栋, 方伟, 宋宝奇, 叶新, 王凯. 太阳辐射计的衍射效应修正[J]. 中国光学, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
引用本文: 刘国栋, 方伟, 宋宝奇, 叶新, 王凯. 太阳辐射计的衍射效应修正[J]. 中国光学, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
LIU Guo-dong, FANG Wei, SONG Bao-qi, YE Xin, WANG Kai. Diffraction effect correction of solar radiometer[J]. Chinese Optics, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
Citation: LIU Guo-dong, FANG Wei, SONG Bao-qi, YE Xin, WANG Kai. Diffraction effect correction of solar radiometer[J]. Chinese Optics, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851

太阳辐射计的衍射效应修正

doi: 10.3788/CO.20181105.0851
基金项目: 

国家自然科学基金资助项目 41474161

国家高技术研究发展计划(863计划)资助项目 2015AA123703

详细信息
    作者简介:

    刘国栋(1990-), 男, 河南商丘人, 硕士研究生, 主要从事太阳辐射计衍射效应方面的研究。E-mail:Liugd0512@163.com

    方伟(1965—),女,辽宁朝阳人,博士,研究员,博士生导师,主要从事空间辐射计量方面的研究。E-mail:fangw@ciomp.ac.cn

  • 中图分类号: TH753

Diffraction effect correction of solar radiometer

Funds: 

National Natural Science Foundation of China 41474161

National High Technology Research and Development Program of China 2015AA123703

More Information
  • 摘要: 太阳辐射测量是研究太阳活动与地球气候演变的重要方式之一,对人类社会的可持续发展具有重要意义。衍射效应作为测量过程中系统误差的主要来源之一,有必要进行精确的修正,从而提高测量数据的精度。首先,对衍射效应理论进行研究,从Kirchhoff衍射理论出发,在高斯光学近似下,逐步确定点与点,点与面,面与面之间的能量传输关系,推导出了衍射效应的一般公式;接着,根据衍射效应的渐近性质,得到了一种简化的计算方法;然后,用简化的方法计算太阳辐照绝对辐射计(SIAR)的衍射效应以及衍射修正因子,最后,根据衍射修正结果,计算相对于世界辐射基准(WRR)的定标系数。结果显示:SIAR的衍射效应以及衍射修正因子分别约为1.002 742和0.997 265。经过衍射修正后,SIAR对WRR的定标系数更接近于1,表明衍射修正降低了系统误差,提高了辐射测量的准确度。
  • 图  1  点光源衍射讨论中的参量

    Figure  1.  Parameters in the discussion of point source diffraction

    图  2  SAD结构

    Figure  2.  Structure of SAD

    图  3  SIAR剖面示意图

    Figure  3.  Profile of SIAR

    图  4  光阑的几何分布

    Figure  4.  Positions of the apertures

    表  1  几种太阳辐射计的衍射效应

    Table  1.   Diffraction effects in some radiometers

    Instrument Ra/mm dd/mm Rd/mm Correction
    PMO6V 4.25 95.4 2.5 1.001 280
    DIARAD 6.52 144 4.001 5 1.000 833
    ERBE 12.09 100.8 4.039 1.000 209
    ACRIM
    Baf1 6.654 8 150.469 6 3.987 8 1.000 828
    Baf2 6.311 9 76.352 4 3.987 8 1.000 466
    Total 1.001 295
    TIM 3.989 4 101.6 7.62 0.999 582
    下载: 导出CSV

    表  2  采用简化方法的衍射效应值

    Table  2.   Diffraction effects by using simplified method

    太阳辐射计 衍射效应
    F Cdiff=F-1 δ
    RMO6V 1.001 229 0.001 299 +0.000 010
    DIARAD 1.000 807 0.000 807 -0.000 026
    ERBE 1.000 204 0.000 204 -0.000 005
    ACRIM Baf1 1.000 793 0.000 793 -0.000 035
    Baf2 1.000 450 0.000 450 -0.000 016
    Baf3 1.001 242 0.00 1242 -0.000 053
    TIM 0.999 577 -0.000 423 -0.000 005
    下载: 导出CSV

    表  3  SIAR的SAD参数

    Table  3.   SAD parameters of the SIAR

    光阑 光阑半径R/mm 光阑与探测器距离dd/mm 探测器(主光阑)半径rd/mm
    Ap1 5.75 100 4.00
    Ap2 5.40 80 4.00
    Ap3 5.05 60 4.00
    Ap4 4.70 40 4.00
    Ap5 4.35 20 4.00
    下载: 导出CSV

    表  4  SIAR的衍射效应

    Table  4.   Diffraction effects in SIAR

    光阑 衍射效应F ε=F-1
    Ap1 1.000 785 0.000 785
    Ap2 1.000 692 0.000 692
    Ap3 1.000 577 0.000 577
    Ap4 1.000 435 0.000 435
    Ap5 1.000 252 0.000 252
    Total 1.002 742 0.002 742
    下载: 导出CSV

    表  5  衍射效应未修正的WRR Factor

    Table  5.   Diffraction effects uncorrected WRR Factor

    SIAR型号 定标系数WRR Factor |WRR Faotor-1|
    SIAR-2C 0.998 949 0.001 051
    SIAR-4A 0.998 445 0.001 555
    SIAR-4B 0.996 734 0.003 266
    SIAR-4D 0.995 917 0.004 083
    下载: 导出CSV

    表  6  衍射修正过的WRR Factor

    Table  6.   Diffraction effects corrected WRR Factor

    SIAR型号 定标系数WRR Factor |WRR Factor-1|
    SIAR-2C 1.001 688 0.001 688
    SIAR-4A 1.001 182 0.001 128
    SIAR-4B 0.999 467 0.000 533
    SIAR-4D 0.998 648 0.001 352
    下载: 导出CSV
  • [1] 王玉鹏, 方伟, 弓成虎, 等.双锥腔互补偿型绝对辐射计[J].光学精密工程, 2007, 15(11):1662-1667. doi:  10.3321/j.issn:1004-924x.2007.11.005

    WANG Y P, FANG W, GONG CH H, et al.. Dual cavity inter-compensating absolute radiometer[J]. Opt. Precision Eng., 2007, 15(11):1662-1667.(in Chinese) doi:  10.3321/j.issn:1004-924x.2007.11.005
    [2] 唐潇, 方伟, 骆杨, 等.低温绝对辐射计ARCPR热沉的热特性分析[J].光学学报, 2016, 36(10):1012004. http://cdmd.cnki.com.cn/Article/CDMD-80139-1017064021.htm

    TANG X, FANG W, LUO Y, et al.. Thermal analysis of an absolute cryogenic radiometer heat sink[J]. Acta Optica Sinica, 2016, 36(10):1012004.(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-80139-1017064021.htm
    [3] 李同保.光辐射测量中的衍射修正问题[J].国外计量, 1979, 5:42-46. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000003320094

    LI T B. Diffraction Correction in Optical Radiation Measurement[J]. Foreign Measurement, 1979, 5:42-46.(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000003320094
    [4] BORN M, WOLF E. Principles of Optics[M]. Oxford, UK:Pergamon Press Ltd., 1965.
    [5] WOLF E. Light distribution near focus in an error-free diffraction image[J]. Proc. Roy. Soc., 1951, A204:533-548. http://www.jstor.org/stable/98739
    [6] SHIRLEY E L. Revised formulas for diffraction effects with point and extended sources[J]. Appl. Opt., 1998, 37:6581-6590. doi:  10.1364/AO.37.006581
    [7] EDWARDS P, MCCALL M. Diffraction loss in radiometry[J]. Appl. Opt., 2003, 42:5024-5032. doi:  10.1364/AO.42.005024
    [8] SHIRLEY E L. Diffraction corrections in radiometry:Spectral and total power, and asymptotic properties[J]. J. Opt. Soc. Am., 2004, A21:1895-1906.
    [9] FOCKE J. Total illumination in an aberration-free diffraction image[J]. Optica Acta, 1956, 3:161-163. doi:  10.1080/713823678
    [10] OLVER F W J. Asymptotics and Special Functions[M].Wellesley, USA:A.K.Peters, Ltd., 1997.
    [11] 方伟, 禹秉熙, 王玉鹏, 等.太阳辐照绝对辐射计及其在航天器上的太阳辐照度测量[J].中国光学, 2009, 2(1):23-28. doi:  10.3969/j.issn.2095-1531.2009.01.004

    FANG W, YU B X, WANG Y P, et al.. Solar irradiance absolute radiometers and solar irradiance measurement on spacecraft[J]. Chinese Journal of Optics, 2009, 2(1):23-28.(in Chinese) doi:  10.3969/j.issn.2095-1531.2009.01.004
    [12] BUTLER J J, JOHNSON B C, RICE J P, et al.. Sources of differences in on-orbital total solar irradiance measurements and description of a proposed laboratory intercomparison[J]. J. Res. Natl. Inst. Stand. Technol., 2008, 113:187-203. doi:  10.6028/jres
    [13] 王玉鹏, 方伟, 禹秉熙.宽视场绝对辐射计的杂散光特性[J].光学精密工程, 2008, 16(3):432-436. http://d.old.wanfangdata.com.cn/Periodical/gxjmgc200803011

    WANG Y P, FANG W, YU B X. Stray light characteristics of WFOV absolute radiometer[J]. Opt. Precision Eng., 2008, 16(3):432-436.(in Chinese) http://d.old.wanfangdata.com.cn/Periodical/gxjmgc200803011
    [14] KOPP G, LEAN J L. A new, lower value of total solar irradiance:Evidence and climate significance[J]. Geophysical Research Letters, 2011, 38:L01706. doi:  10.1029-2010GL045777/
    [15] MARKUS S. Office of Scientific & Technical Information Technical Reports[R]. Advances in Solar Radiometry, 2014.
    [16] FINSTERLE W. International pyrheliometer comparison[R]. IPC-XⅡ, 2015.
  • [1] 吕强, 王玮, 刘兆武, 宋莹, 姜珊, 刘林, 巴音贺希格, 李文昊.  五维自由度衍射光栅精密测量系统 . 中国光学, 2020, 13(1): 189-202. doi: 10.3788/CO.20201301.0189
    [2] 赵维, 刘华, 陆子凤, 卢振武, 王新.  大口径衍射望远系统初始结构研究 . 中国光学, 2019, 12(6): 1395-1402. doi: 10.3788/CO.20191206.1395
    [3] 胡琪, 王喆, 刘洪顺, 薛智文, 邓家春.  基于单次傅里叶变换的分段衍射算法 . 中国光学, 2018, 11(4): 568-575. doi: 10.3788/CO.20181104.0568
    [4] 马国鹭, 刘丽贤, 杨贵洋, 赵斌.  参数可控长距无衍射光束的生成方法研究 . 中国光学, 2018, 11(1): 100-107. doi: 10.3788/CO.20181101.0100
    [5] 吕强, 李文昊, 巴音贺希格, 柏杨, 刘兆武, 王玮.  基于衍射光栅的干涉式精密位移测量系统 . 中国光学, 2017, 10(1): 39-50. doi: 10.3788/CO.20171001.0039
    [6] 李瑶, 杨甬英, 王晨, 陈元恺, 陈晓钰.  点衍射干涉检测技术 . 中国光学, 2017, 10(4): 391-414. doi: 10.3788/CO.20171004.0391
    [7] 代晓珂, 金春水, 于杰.  点衍射干涉仪波面参考源误差及公差分析 . 中国光学, 2014, 7(5): 855-862. doi: 10.3788/CO.20140705.0855
    [8] 杨虹, 黄远辉, 龚昌妹, 吴腾飞, 邵晓鹏.  散射介质超衍射极限技术研究进展 . 中国光学, 2014, 7(1): 1-25. doi: 10.3788/CO.20140701.001
    [9] 张梅, 耿征.  应用折-衍射元件校正人眼色差 . 中国光学, 2012, 5(5): 525-530. doi: 10.3788/CO.20120505.0525
    [10] 王红睿, 王玉鹏.  用于测量太阳总辐射的空间辐射计 . 中国光学, 2012, 5(6): 555-565. doi: 10.3788/CO.20120506.0555
    [11] 李淼, 高慧斌.  标校经纬仪航向标定修正方法 . 中国光学, 2011, 4(6): 583-587.
    [12] 王红睿, 方伟.  自动寻日的太阳辐照绝对辐射计 . 中国光学, 2011, 4(3): 252-258.
    [13] 张海涛.  基于光学设计软件的相移点衍射干涉仪建模 . 中国光学, 2010, 3(6): 616-622.
    [14] 曲艺.  真空紫外波段光栅二级衍射效率测量装置 . 中国光学, 2010, 3(6): 566-571.
    [15] 喻波.  Mo/Si多层膜小角X射线衍射结构表征 . 中国光学, 2010, 3(6): 623-629.
    [16] 曹艳波, 艾华.  亚波长闪耀光栅矢量衍射效率计算 . 中国光学, 2010, 3(6): 679-683.
    [17] 吕银环, 张涛, 童广辉, 崔维鑫.  基于衍射光学元件的红外弱目标探测系统设计 . 中国光学, 2009, 2(6): 543-549.
    [18] 黄启禄, 吴逢铁, 卢文和.  被动方式产生近似无衍射贝塞尔光技术 . 中国光学, 2009, 2(4): 289-295.
    [19] 方伟, 禹秉熙, 王玉鹏, 弓成虎, 杨东军, 叶新, .  太阳辐照绝对辐射计及其在航天器上的太阳辐照度测量 . 中国光学, 2009, 2(1): 23-28.
    [20] 徐福全, 金陆, 李文昊, 裴舒, 巴音贺希格, 齐向东.  曝光系统离焦对平面全息光栅衍射波前的影响 . 中国光学, 2008, 1(1): 57-61.
  • 加载中
图(4) / 表ll (6)
计量
  • 文章访问数:  478
  • HTML全文浏览量:  119
  • PDF下载量:  147
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-18
  • 修回日期:  2018-02-21
  • 刊出日期:  2018-10-01

太阳辐射计的衍射效应修正

doi: 10.3788/CO.20181105.0851
    基金项目:

    国家自然科学基金资助项目 41474161

    国家高技术研究发展计划(863计划)资助项目 2015AA123703

    作者简介:

    刘国栋(1990-), 男, 河南商丘人, 硕士研究生, 主要从事太阳辐射计衍射效应方面的研究。E-mail:Liugd0512@163.com

    方伟(1965—),女,辽宁朝阳人,博士,研究员,博士生导师,主要从事空间辐射计量方面的研究。E-mail:fangw@ciomp.ac.cn

  • 中图分类号: TH753

摘要: 太阳辐射测量是研究太阳活动与地球气候演变的重要方式之一,对人类社会的可持续发展具有重要意义。衍射效应作为测量过程中系统误差的主要来源之一,有必要进行精确的修正,从而提高测量数据的精度。首先,对衍射效应理论进行研究,从Kirchhoff衍射理论出发,在高斯光学近似下,逐步确定点与点,点与面,面与面之间的能量传输关系,推导出了衍射效应的一般公式;接着,根据衍射效应的渐近性质,得到了一种简化的计算方法;然后,用简化的方法计算太阳辐照绝对辐射计(SIAR)的衍射效应以及衍射修正因子,最后,根据衍射修正结果,计算相对于世界辐射基准(WRR)的定标系数。结果显示:SIAR的衍射效应以及衍射修正因子分别约为1.002 742和0.997 265。经过衍射修正后,SIAR对WRR的定标系数更接近于1,表明衍射修正降低了系统误差,提高了辐射测量的准确度。

English Abstract

刘国栋, 方伟, 宋宝奇, 叶新, 王凯. 太阳辐射计的衍射效应修正[J]. 中国光学, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
引用本文: 刘国栋, 方伟, 宋宝奇, 叶新, 王凯. 太阳辐射计的衍射效应修正[J]. 中国光学, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
LIU Guo-dong, FANG Wei, SONG Bao-qi, YE Xin, WANG Kai. Diffraction effect correction of solar radiometer[J]. Chinese Optics, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
Citation: LIU Guo-dong, FANG Wei, SONG Bao-qi, YE Xin, WANG Kai. Diffraction effect correction of solar radiometer[J]. Chinese Optics, 2018, 11(5): 851-859. doi: 10.3788/CO.20181105.0851
    • 太阳辐射是地球最主要的能量来源,除了核能,火山爆发等,地球上大部分的能量都直接或间接地来自太阳辐射,太阳辐射对地球的生态和环境具有深刻的影响[1]。为了自身的长久可持续发展,人类很早就开始了对太阳辐射的监测与研究。太阳辐射的测量主要是指太阳总辐照度(Total Solar Irradiance, TSI)测量和太阳光谱辐照度(Solar Spectral Irradiance, SSI)测量。这些太阳辐射的监测数据,可以应用到很多科学领域,如气候模型的推演、太阳物理研究、地球辐射收支平衡的研究等等,随着研究的深入,这些领域对太阳辐射数据的测量精度要求越来越高。

      太阳绝对辐射计是一种对太阳总辐照度进行监测的仪器[2],其测量精度直接影响太阳辐射数据的可信度。衍射效应是太阳绝对辐射计系统误差的主要来源之一,要获得精确的测量数据,必须对其进行修正。但是在辐射计量领域,衍射效应的修正一直是一个难以解决的问题。一方面是因为衍射效应广泛存在于光学系统中,几乎所有对光起限制作用的器件,都会引起衍射效应;另一方面,衍射效应难以用实验测量,即使能测量,实验中的不确定度也超过能容忍的范围。因此,对衍射效应的研究主要侧重于数学计算上。

      辐射计量中衍射效应数学计算在国外有较多的研究。1962年,C.L.Sanders和O.C.Jones在用铂点黑体实现光单位时,用夫琅和费衍射理论证明了接收像面中心照度比通过几何光学计算结果低大约0.25%,证明了衍射效应对精密辐射测量的显著性影响。W.R.Blevin在实测Stefan-Boltzman常数的工作中发现这项误差甚至可达到0.5%。此后,衍射效应成为黑体辐射测量中不可忽视的因素[3]。近几十年来,随着人类活动走向太空,发射了各种太阳辐射监测仪器,由于对其精度的要求逐渐提高,衍射效应的数学计算,逐渐开展起来。目前在国外,衍射效应的修正已经开始应用于辐射计的测量结果,获得较高的辐射测量精度,然而国内却尚未开展该方面的工作,严重制约了我国在辐射计量领域的发展。本文根据对衍射效应问题的研究,提出一种实用化的计算方法,精度高,计算简便,实现了对太阳辐照绝对辐射计(Solar Irradiance Absolute Radiometer, SIAR)的衍射效应修正。

    • 为了便于论述,本文做如下定义:衍射效应是指由于光的衍射,探测器接收的辐射通量与几何光学期望值的比,是一个与光的衍射有关的“物理量”,用F表示。

      根据Kirchhoff衍射理论和傍轴光学近似,点光源通过光阑衍射后的电场强度为:

      (1)

      其中,C=(d-1+d-1)/2是与结构有关的系数,相关参数如图 1所示,它表示了x′与x电场强度的相互关系。考虑到能量密度W=C|U(x, t)|2,并做一定的坐标变换,可以得到:

      (2)

      图  1  点光源衍射讨论中的参量

      Figure 1.  Parameters in the discussion of point source diffraction

      其中,R表示光阑的半径,r表示x′到光轴的距离,ρ=r″/Ru=。2qCR2v=|2qCRr|。

      式(2)表示了x′与x之间能量密度的关系。而在几何光学中,x′与x之间能量密度的关系如下:

      (3)

      式中,Θ(x)是Heaviside阶跃函数。

      因此,|S(r′)/SG(0)|2就表示点x′处的衍射效应,其表达式由Lommel函数[4]给出:

      (4)
      (5)

      式中,Un(u, v)和Vn(u, v)可以用n阶贝塞尔函数来表示

      (6)
      (7)

      式(4)和式(5)是等价的,当r′>R时,式(4)更简便,而当r′>R时,式(5)更简便。

      实际的辐射测量,是对整个探测平面的能量密度进行积分,为此,引入一个新的量L(u, v0),表示整个探测平面的衍射效应[5],它相当于|S(r′)/SG(0)|2在整个探测平面的积分,类似于|S(r′)/SG(0)|2的推导,Wolf给出了L(u, v0)的表达式:

      v0>|u|时,

      (8)

      v0<|u|时,

      (9)

      式中,Q2s(v0)和Yn(u, v0)表示如下:

      设探测器平面上所有点的集合为{x′},探测器边缘上的点满足v=v0,探测器内部的点满足vv0。显然L(u, v)表示以探测平面中心为原点,以v/(2qCR)为半径的圆形区域的衍射效应。

      实物光源常常可以看作无数个相互独立的点光源的集合,在式(8)和式(9)的基础上,E.L.Shirley推导出了适用于扩展光源的公式[6]:

      (10)

      式中,D=4π3R4rs2rd2/[ds2dd2(λvm)2],Φλ(λ)表示探测器接收的总能量,Lλ(λ)表示光源的亮度。

      若仅仅在几何光学范畴考虑,则可以很容易求得探测器接收的总能量ΦG(λ),其表达式跟具体的光学结构有关。因此,总的衍射效应就可以表示为:

      (11)

      现在通过图 2对式(10)的参数进行说明,可以看出,该光构中包含了扩展型朗伯光源(Source),光阑(Aperture)和探测器(Detector)三个要素,Rrsrddsdd表示3个要素构成的几何参数。事实上,对该结构的研究是辐射计量中研究最多的问题之一,称为SAD问题。

      图  2  SAD结构

      Figure 2.  Structure of SAD

      在SAD问题中常常引入如下的参数:u= vm=max(vs, vd), ,使用这5个参数,可以将SAD的衍射问题分成3种,其衍射效应值分别用F1F2F3表示。下式给出了F1F2F3所对应的条件[7]

      (12)

      对于F1情况,光阑会限制光的几何传输路径,而在F2中则不会。F3情况会导致探测器边缘位于光束与阴影区域的交界附近,衍射效应分析复杂,是结构设计中需要避免的情况。式(8)和式(9)是L(u, v)分别在F1F2情况下的表达式,只不过取v=vm(1+σx)。在几何光学范畴内考虑,ΦG(λ)在F1F2情况下可以表示为

      (13)

      综合式(10)、式(11)和式(13)可得总的衍射效应为:

      (14)

      该式表示了SAD结构对波长为λ的光的衍射效应,然而太阳辐射计中所测光源是太阳,是一个宽光谱的普朗克光源。实际操作中,可以将太阳辐射的等效波长λeff≈902.792 042 6 nm应用于式(14)[8]

    • 由于L(u, v)包含复杂的n阶贝塞尔函数,对式(14)直接积分是困难的。实际操作中,常常根据贝塞尔函数的渐近性质对L(u, v)进行简化。

      v>|u|时,Focke在前人的工作基础上给出了L(u, v)的渐近形式[9]

      (15)

      v<|u|时,则L(u, v)的表达式依赖于Q2s(v)以及Yn(u, v)的渐近形式,具体如下[10]

      (16)
      (17)
      (18)

      其中σk是由简单函数Wk(x)决定的变量

      (19)
      (20)

      通过上述渐近处理,可以将L(u, v)表示为简单函数,则式(14)就转化为简单函数的线性积分,大大简化了计算量。

      为了验证上述简化方法的准确性,本文参考国外一些已经进行过衍射修正的辐射计结构,将使用本文方法计算的结果与其报告的结果进行对比。表 1是国外几个辐射计的SAD结构参数以及衍射效应的计算结果[11]。计算过程是对公式(14)的数值计算,该过程采用精确的L(u, v)(含有n阶贝塞尔函数),计算结果较为精确。可以看出太阳总辐照度测量中衍射效应在10-3~10-4之间,已经成为精确测量中不可忽视的重要因素。

      表 1  几种太阳辐射计的衍射效应

      Table 1.  Diffraction effects in some radiometers

      Instrument Ra/mm dd/mm Rd/mm Correction
      PMO6V 4.25 95.4 2.5 1.001 280
      DIARAD 6.52 144 4.001 5 1.000 833
      ERBE 12.09 100.8 4.039 1.000 209
      ACRIM
      Baf1 6.654 8 150.469 6 3.987 8 1.000 828
      Baf2 6.311 9 76.352 4 3.987 8 1.000 466
      Total 1.001 295
      TIM 3.989 4 101.6 7.62 0.999 582

      表 2是采用本文简化方法计算得到的衍射效应,也就是考虑L(u, v)的渐近性质,进行较为简便的线性积分。可以看出,计算结果与表 1的偏差在10-5~10-6之间,远小于衍射效应值,从而证明本文的方法在简化运算的同时兼顾了计算结果的准确性。

      表 2  采用简化方法的衍射效应值

      Table 2.  Diffraction effects by using simplified method

      太阳辐射计 衍射效应
      F Cdiff=F-1 δ
      RMO6V 1.001 229 0.001 299 +0.000 010
      DIARAD 1.000 807 0.000 807 -0.000 026
      ERBE 1.000 204 0.000 204 -0.000 005
      ACRIM Baf1 1.000 793 0.000 793 -0.000 035
      Baf2 1.000 450 0.000 450 -0.000 016
      Baf3 1.001 242 0.00 1242 -0.000 053
      TIM 0.999 577 -0.000 423 -0.000 005
    • SIAR是一种利用电替代原理对太阳辐射进行测量的仪器,其构造如图 3所示,主要结构包括锥腔,光阑结构,热沉以及相关的电路模块[12]

      图  3  SIAR剖面示意图

      Figure 3.  Profile of SIAR

      SIAR的光阑有6个:主光阑,视场光阑,4个消杂光光阑。视场光阑在前,主光阑在后并靠近锥腔,半径为4 mm。整体结构视场半角为1°,相邻光阑距离为20 mm光阑分布如图 4所示。

      图  4  光阑的几何分布

      Figure 4.  Positions of the apertures

      应用上述SAD结构,对SIAR的衍射效应进行分析。首先由于主光阑靠近锥腔,认为通过主光阑的光都能到达锥腔,则主光阑开口相当于SAD的探测器,且rd=4 mm;太阳作为SAD结构中的扩展型朗伯光源,rs≈6.75×1011 mm, 距离光阑距离ds≈1.5×1014 mm, 由于太阳到地球表面的视场半角约为32′, 位置在前的光阑不会遮挡后面光阑的几何光学路线[13],因此视场光阑和4个消杂光光阑可以作为SAD结构中的非限制光阑,引起F2类型的衍射效应。SIAR的衍射效应问题可以分解为5个SAD问题。视场光阑与随后的4个消杂光光阑用Apx(x=1, 2, 3, 4, 5)表示,具体参数如表 3所示。

      表 3  SIAR的SAD参数

      Table 3.  SAD parameters of the SIAR

      光阑 光阑半径R/mm 光阑与探测器距离dd/mm 探测器(主光阑)半径rd/mm
      Ap1 5.75 100 4.00
      Ap2 5.40 80 4.00
      Ap3 5.05 60 4.00
      Ap4 4.70 40 4.00
      Ap5 4.35 20 4.00

      根据L(u, v)在v<|u|时的渐近性质,利用式(14)进行数值积分得到每个光阑引起的衍射效应以及总的衍射效应,如表 4所示。

      表 4  SIAR的衍射效应

      Table 4.  Diffraction effects in SIAR

      光阑 衍射效应F ε=F-1
      Ap1 1.000 785 0.000 785
      Ap2 1.000 692 0.000 692
      Ap3 1.000 577 0.000 577
      Ap4 1.000 435 0.000 435
      Ap5 1.000 252 0.000 252
      Total 1.002 742 0.002 742

      表 4可以看出,SIAR总的衍射效应达到0.274%,而根据90年代给出的太阳总辐照度为(1 365.4±1.3) W/m2,不确定度占0.095%,由此可见,由衍射效应引起的误差是相当显著的[14],是精确测量中不可忽视的重要因素。另外由于前一级光阑不会遮挡后一级光阑的几何传输路径,衍射效应会随消杂光光阑个数的增加而逐级累加。

      辐射计的衍射修正因子Cdiff定义为[15]

      (21)

      则SIAR的衍射修正因子Cdiff=≈0.997 265,辐射计测量结果乘以Cdiff就得到经过衍射修正后的太阳辐射数据。

      太阳辐射计所测辐射数据,需要定标到国际辐射基准(WRR)。具体做法为:在相同时间,相同实验环境下,SIAR与WRR同时对太阳辐射进行监测,得到的结果分别为TSI|SIAR, TSI|WRR,则SIAR的定标系数为

      (21)

      WRR Factor越接近1,则表明辐射计的测量准确度越高。表 5是未经过衍射修正的SIAR的定标系数[16]

      表 5  衍射效应未修正的WRR Factor

      Table 5.  Diffraction effects uncorrected WRR Factor

      SIAR型号 定标系数WRR Factor |WRR Faotor-1|
      SIAR-2C 0.998 949 0.001 051
      SIAR-4A 0.998 445 0.001 555
      SIAR-4B 0.996 734 0.003 266
      SIAR-4D 0.995 917 0.004 083

      对SIAR进行衍射效应修正,则定标系数为

      (22)

      表 6是经过衍射修正后SIAR的定标系数,可以看出经过衍射修正后,SIAR-4B和SIAR-4D的定标系数都明显地趋近于1。而对于SIAR-2C和SIAR-4A,衍射效应造成的系统误差,使得定标系数从小于1变成大于1,之所以没能更靠近1,是由其他的系统误差导致的。因此,衍射修正降低了辐射计的一部分系统误差,提高了SIAR测量数据的准确度。

      表 6  衍射修正过的WRR Factor

      Table 6.  Diffraction effects corrected WRR Factor

      SIAR型号 定标系数WRR Factor |WRR Factor-1|
      SIAR-2C 1.001 688 0.001 688
      SIAR-4A 1.001 182 0.001 128
      SIAR-4B 0.999 467 0.000 533
      SIAR-4D 0.998 648 0.001 352
    • 本文从Kirchhoff衍射理论出发,根据高斯光学近似,逐步推导出了衍射效应下点与点、点与面、面与面的能量传输关系,最终给出了衍射效应计算的一般公式。为了便于计算,构建了用于衍射分析的SAD模型,重点分析了F1F2两种衍射类型,最后根据渐近性质,提出了简便的计算方法。根据SIAR的光阑结构特点,将其衍射效应分解成5个SAD问题,使用简化的方法分别计算衍射效应,最后按照累加原则,得到总的衍射效应。由于衍射效应与辐射计的衍射修正因子互为倒数关系,从而给出了衍射修正因子,这是国内辐射计首次实现对衍射效应的修正。根据衍射修正结果,计算4个型号的SIAR相对WRR的定标系数,结果表明四个SIAR中衍射效应引起的系统误差被修正,提高了测量结果的准确度。

参考文献 (16)

目录

    /

    返回文章
    返回