A Simplified Radiative Transfer Calculation Scheme for Transmittance in Satellite Thermal Infrared Remote Sensing
-
摘要:
热红外大气辐射传输会受到各种因素的影响，大气透过率是其中一个重要参量。目前，对于热红外谱段（8~14 μm）上大气透过率的研究仍然面临着背景参数输入过多且不易获取的问题。基于此，提出了一种应用于卫星热红外遥感的大气透过率简化辐射传输计算方案。利用MODTRAN模型对大气透过率的影响参数进行整体定量模拟分析，从而得到影响大气透过率的关键参数，由此形成一套基于MODTRAN模型的参数简化方案。在MERSI-Ⅱ两个热红外通道的中心波长处，对MODTRAN模型和参数简化方案计算出的大气透过率做对比验证，R2均大于0.99，RMSE最大为0.007458，具有较高精度。参数简化方案仅依赖于观测天顶角和水汽柱浓度，不需要CO2、O3和气溶胶光学厚度。相比于MODTRAN模型，参数简化方案输入参数减少了，计算效率也随之提高。本研究为热红外数据大气校正提供了理论支撑。
Abstract:Atmospheric radiative transfer in the thermal infrared spectrum is influenced by various factors, and atmospheric transmittance is a critical parameter. Currently, researching atmospheric transmittance within 8-14 μm thermal infrared spectrum remains challenging because of numerous elusive input background parameters. Consequently, a simplified parameterization scheme for thermal infrared radiative transfer in satellite remote sensing is proposed. Using the moderate resolution atmospheric transmission (MODTRAN) model, the impact parameters of atmospheric transmittance were quantitatively simulated and analyzed, to identify key parameters, leading to the development of a MODTRAN-model-based simplified scheme. Atmospheric transmittance values calculated by the MODTRAN model and simplified scheme were compared and validated at the central wavelengths of the two thermal infrared channels of a medium resolution spectral imager-LL (MERSI-LL). The R2 value exceeded 0.99, and the root mean-square error (RMSE) was below 0.007458, indicating high accuracy. The simplified scheme relies solely on the view zenith angle and water vapor column, eliminating the need for CO2, O3, and aerosol optical depth. Compared with the MODTRAN model, the simplified scheme reduces the number of input parameters and enhances computational efficiency. This study offers theoretical support for atmospheric correction of thermal infrared data.
-
0. 引言
有机电致发光器件（Organic Light Emitting Device，OLED）具有发光亮度高、响应时间短、可视范围大和可柔性化等优点，被称为“梦幻般的显示器”，被视为液晶显示后的下一代主流显示器，并初步应用于装饰和室内照明[1-6]。近年来，高性能顶发射器件逐渐成为研究热点，诸多科研工作者投身于实现高性能器件的研究中，目前主要从两个方面入手：一是新材料的研发，如新型有机发光分子材料[7]；二是新结构的开发，如超薄结构[8]、量子阱结构[9]和和微腔结构[10]等。在微腔结构方面，主要是通过理论计算改变有机结构层厚度，进而调节器件的微腔长度，获得不同模数的微腔，使器件处于不同微腔加强区，从而提升器件性能。
光学微腔是一种光学微型谐振腔，尺寸在光波长量级。有机微腔电致发光器件最早是日本九州大学在1993年完成的[11]。当前关于有机微腔发光的大部分研究致力于提升器件效率[12-14]，而对具有微腔效应顶发射器件的色纯度及稳定性的研究存在不足。因此，本文在现有器件研究的基础上，通过引入二阶微腔结构[15-16]，制备了一系列顶发射微型器件，验证二阶微腔长度范围内器件的光电性能，最终获得优化后的稳定绿光顶发射器件，实现标准绿光显示。
1. 实验
本文所制备的顶发射器件，微腔结构为简单的FP（Fabry-Perot）微腔结构[17-19]，底部全反射电极采用Ag，顶部光出射端采用半透明的金属阴极Mg/Ag作为半反射镜。器件各膜层通过蒸镀设备依次完成，主要膜层及所用材料见表 1，其中阳极为ITO，空穴注入层（Hole Injection Layer, HIL）为有机材料F16CuPc和NPB，F16CuPc为掺杂料；空穴传输层（Hole Transport Layer, HTL）为有机材料NPB；电子阻挡层（Electron Blocking Layer, EBL）为有机材料TCTA；有机发光层（Emitting Layer, EML）为有机材料mCP和Ir(ppy)3，mCP为绿色发光基质，Ir(ppy)3掺杂料；电子传输层（Electron Transport Layer, ETL）为有机材料Bphen和Liq，Liq为掺杂料；光输出耦合层（Capping Layer, CPL）为有机材料Alq3。器件中涉及的有机材料分子结构如图 1所示。
表 1 器件主要膜层及所用材料Table 1. Layers and materials of deviceLayer Material anode ITO HIL Copper(II)1, 2, 3, 4, 8, 9, 10, 11, 15, 16, 17, 18, 22, 23, 24, 25-hexadecafluoro-29H, 31H-phthalocyanine（F16CuPc）
N, N'-Di-[(1-naphthyl)-N, N'-diphenyl]-1, 1'-biphenyl)-4, 4'-diamine (NPB)HTL N, N'-Di-[(1-naphthyl)-N, N'-diphenyl]-1, 1'-biphenyl)-4, 4'-diamine (NPB) EBL 4, 4', 4''-tris(carbazol-9-yl)-triphenylamine (TCTA) EML 1, 3-bis(9-carbazolyl)benzene（mCP）
Iridium, tris[2-(2-pyridinyl-kN)phenyl-kC]（Ir(ppy)3）ETL 4, 7-Diphenyl-1, 10-phenanthroline（Bphen）
8-hydroxyquinoline lithium（Liq）cathode Mg/Ag CPL 8-Hydroxyquinoline aluminum salt（Alq3） ![]()
图 1 器件中涉及的有机材料分子结构Figure 1. Molecular structures of the materials in the OLED devices该器件采用云南北方奥雷德光电股份有限公司开发的硅基CMOS基板作为器件衬底，依次蒸镀各层有机材料，蒸发速率保持在0.1 nm/s，真空度保持在2×10-4 Pa。器件的亮度及光谱通过PR-655测量，电流和电压采用搭载Keithley 2400测试仪的测试系统进行测量。
2. 结果讨论
2.1 微腔长度对色纯度影响
一般来说，顶发射器件都存在微腔效应，器件发出的光谱强度I(λ)如式(1)[20]：
$$ I\left( \lambda \right) = \frac{{\left( {1 + {R_{\text{h}}}} \right)\left[ {1 + {R_{\text{f}}} + 2\sqrt {{R_{\text{f}}}} \cos \left( {\frac{{4{\rm{ \mathsf{ π} }}Z}}{\lambda }} \right)} \right]}}{{1 + {R_{\text{f}}}{R_{\text{h}}} - 2\sqrt {{R_{\text{f}}}{R_{\text{h}}}} \cos \left( {\frac{{4{\rm{ \mathsf{ π} }}L}}{\lambda }} \right)}}{I_0}\left( \lambda \right) $$ (1) 式中：Rf为全反射镜的反射率；Rh为半透明反射镜的反射率；I0(λ)为自由空间的光谱强度；L为器件微腔光学长度；Z为全反射镜与有机发光层之间的距离。其中，微腔的光学长度L计算式为：
$$ L = \sum {{n_{\text{m}}}{d_{\text{m}}}} + {n_{{\text{ITO}}}}{d_{{\text{ITO}}}} + \left| {\frac{{{\lambda _q}}}{{4{\rm{ \mathsf{ π} }}}}\sum\limits_i {{\phi _i}\left( \lambda \right)} } \right| = q\frac{{{\lambda _q}}}{2} $$ (2) 式中：nm、dm分别为有机材料的折射率和厚度；nITO、dITO分别为ITO的折射率和厚度；q（1, 2, 3, 4, …）是发射模的模（阶）数；λq是模（阶）数为q的共振发射波长；ϕt(λ)为光在有机界面/金属镜面之间的相移，i为阳极/有机界面或阴极/有机界面。由式(1)、(2)可知，通过调节有机材料膜层厚度，可以改变器件微腔长度，使腔模q的位置产生移动，从而改变微腔器件的出射光波长。为了使器件微腔的谐振波长与发光层电致发光谱的峰值波长相匹配以实现增益，利用公式(2)计算得到一阶腔长对应的有机层总厚度约为100 nm，二阶腔长对应的有机层总厚度约为250 nm。
通过调整空穴传输层和电子阻挡层厚度，实验中制作了5种不同微腔长度的器件A~E，如图 2所示。其结构为：Si Substrate/Ag/ITO/ NPB: F16CuPc(10 nm, 3%)/NPB(x nm)/TCTA(y nm)/ mCP: Ir(ppy)3(40 nm, 6%)/ Bphen: Liq(30 nm, 40%)/ Mg/Ag(12 nm)/Alq3(35 nm)，x表示空穴传输层（NPB）的膜层厚度，y表示电子阻挡层（TCTA）的膜层厚度。其中x分别为30、30、60、20、120，y分别为20、15、20、15、40，器件有机层厚度依次为130 nm、125 nm、160 nm、115 nm、240 nm。
![]()
图 2 5种不同微腔长度器件结构图Figure 2. Schematics of device structure with five microcavity lengths图 3为不同腔长器件EL光谱。器件A、B、C、D在524 nm处有一强峰，556 nm、552 nm、560 nm、560 nm处出现一弱峰，器件E为520 nm处唯一单峰。从图中可以看出，器件C→A→B→D→E长波一侧出现明显的窄化趋势，向短波一侧移动，出现蓝移，560 nm处的肩峰逐渐减弱至消失。这一现象是器件微腔效应导致的，根据腔量子电动力学效应，腔内光场的模式密度受到调制，在谐振波长处得到增强，而在其他波长处的受到抑制，光谱得到窄化[21]。微腔效应的强弱常通过半高宽（FWHM, full width at half maximum）来衡量，计算得到器件C→A→B→D→E半高宽从84 nm减小到33 nm，微腔效应逐渐增强。
不同腔长器件的发光性能如表 2所示。在A~E中，D在亮度、电流效率与外量子效率等方面表现较佳，B次之，C表现最差，而E色坐标偏移最小。这主要是因为，D位于一阶加强区，E位于二阶加强区，C远离加强区。可以看出，当器件腔长位于一阶加强区时，器件的光电效率会得到加强；当位于二阶加强区时，器件效率会低于一阶加强区[22-23]，但器件色纯度明显高于一阶加强区，说明处于二阶加强区对器件的色纯度有显著的提升作用。
表 2 不同腔长器件的光电特性Table 2. Optoectronic performance of device with different cavity lengthsDevice Luminance/(cd/m2) Current efficiency/(cd/A) Peak wavelength/nm FWHM/nm External quantum efficiency/% CIEx, y Color shift[CIE 1931] A 6330 33.80 524 73 9.19% (0.3713, 0.6019) (0.1613, 0.1081) B 7439 39.73 524 70 10.59% (0.3601, 0.6110) (0.1501, 0.0990) C 2198 11.74 524 84 3.39% (0.3959, 0.5821) (0.1859, 0.1279) D 9123 48.72 524 66 12.75% (0.3436, 0.6243) (0.1336, 0.0857) E 5477 29.25 520 33 7.67% (0.2092, 0.7167) (0.0008, 0.0067) 通过进一步的测试发现，制作得到的器件色坐标都具有很好的稳定性，如图 4所示。A~E色坐标CIEx，CIEy在低电压阶段经过短暂上升，电压达到2.8 V后，色坐标保持平稳。从整个变化情况来看，器件E色坐标出现了明显的突变，CIEx骤降到0.2左右，CIEy骤升到0.71左右，出现该现象的原因是器件A~D分别在556 nm、552 nm、560 nm、560 nm处存在一弱峰，导致色坐标产生偏离，发光时表现出黄绿光，而器件E为唯一单峰，在器件正常启亮后就表现出近乎接近标准绿光(0.21, 0.71)显示，如图 4(c)所示。这一结果也再次表明微腔长度处于二阶加强区，对器件发光色纯度有明显的提升作用。
![]()
图 4 不同腔长器件色坐标变化Figure 4. Color coordinate variation of device with different cavity lengths2.2 空穴传输层和电子阻挡层厚度对腔长影响
前述结果表明，当器件微腔长度位于二阶加强区时，器件的色纯度会得到明显提升。为了验证器件处于二阶加强区时，空穴传输层和电子阻挡层厚度是否对微腔长度改变起同等作用，制作了器件E1。在其他条件保持不变的情况下，空穴传输层厚度为40 nm，电子阻挡层厚度为120 nm。从表 3可以看出，E、E1在亮度、电流效率、外量子效率等性能方面表现相当，差异很小。通过光谱图（图 5）和色坐标（图 6）也可以看出，两者EL光谱基本重合，且CIEx、CIEy未发生较大改变。这一结果表明，空穴传输层与电子传输层厚度在微腔长度改变中作用相同，均能有效调节色纯度。
表 3 不同HTL & EBL厚度器件的光电特性Table 3. Optoectronic performance of device with different HTL & EBL thicknessDevice Luminance/(cd/m2) Current efficiency/(cd/A) Peak wavelength/nm FWHM/nm External quantum efficiency/% CIEx, y Color shift[CIE 1931] E 5477 29.25 520 33 7.67 (0.2092, 0.7167) (0.0008, 0.0067) E1 5261 28.09 520 32 7.58 (0.2079, 0.7173) (0.0021, 0.0073) ![]()
图 6 HTL&EBL厚度对色坐标影响Figure 6. Color coordinate variation of device with different HTL&EBL thickness3. 结论
研究发现器件结构为Si Substrate/Ag/ITO/ NPB: F16CuPc(10 nm, 3%)/NPB(x nm)/TCTA(y nm)/ mCP: Ir(ppy)3(40 nm, 6%)/Bphen: Liq(30 nm, 40%)/ Mg/ Ag(12 nm)/Alq3(35 nm)的顶发射绿光器件，通过调节器件空穴传输层和电子阻挡层的厚度使器件处于第二阶微腔加强区，可以使光谱明显窄化，器件色纯度得到极大提升，进一步研究发现，空穴传输层与电子阻挡层在微腔长度改变中作用相同，均能有效调节色纯度。器件在腔长为240 nm时，能实现稳定的高色纯度绿光显示，正向出射绿光的色坐标达到了(0.2092，0.7167)，接近标准绿光(0.21, 0.71)，该结果对二阶腔长绿光器件的应用有较好的参考意义。
-
![]()
![]()
图 4 CO2混合比影响下的平均AT以及误差带
Figure 4. Average AT and error bands at different CO2 mixture ratio
![]()
图 6 O3柱浓度影响下的平均AT以及误差带
Figure 6. Average AT and error bands under different O3 column concentration
![]()
图 8 水汽柱浓度影响下的平均AT以及误差带
Figure 8. Average AT and error bands under different water vapor column
![]()
图 12 不同大气模式下的平均RPC和平均AT以及误差带
Figure 12. Average RPC and average AT and error bands under different atmospheric modes
![]()
图 14 不同气溶胶模式下的平均RPC和平均AT以及误差带
Figure 14. Average RPC and average AT and error bands under different aerosol modes
![]()
图 15 MODTRAN模型与参数简化方案在MERSI-Ⅱ 2个TIR通道的AT比较
Figure 15. Comparison of AT between MODTRAN model and simplified scheme in two MERSI-Ⅱ TIR channels
表 1 辐射传输计算的参数设置
Table 1 Parameter settings for radiative transfer calculation
Parameters Settings Atmospheric profile Mid-latitude summer Aerosol extinction URBAN extinction CO2 mixture ratio 420 ppmv Water vapor column 3.0 g/cm2 O3 column concentrations 7.0 g/m2 AOD 0.4 VZA 0°~60°，interval 5° Sensor altitude 700 km DEM 2 m Wavenumber range 714~1250 cm－1 
下载: 导出CSV
-
[1] 赵帅阳, 景欣, 张飞舟, 等. 中波红外大气透过率特性分析[J]. 激光与红外, 2017, 47(5): 568-574. DOI: 10.3969/j.issn.1001-5078.2017.05.009 ZHAO Shuaiyang, JING Xin, ZHANG Feizhou, et al. Analysis on atmospheric transmittance characteristics in middle infrared (3~5 μm) spectrum[J]. Laser & Infrared, 2017, 47(5): 568-574. DOI: 10.3969/j.issn.1001-5078.2017.05.009
[2] 张安琪, 邰会强, 高佳星, 等. 基于大气参数修正的时变大气透过率测试方法研究[J]. 燃气涡轮试验与研究, 2023, 36(2): 39-44. DOI: 10.3969/j.issn.1672-2620.2023.02.007 ZHANG Anqi, TAI Huiqiang, GAO Jiaxing, et al. Time-varying atmospheric transmittance measurement method based on atmospheric parameter correction[J]. Gas Turbine Experiment and Research, 2023, 36(2): 39-44. DOI: 10.3969/j.issn.1672-2620.2023.02.007
[3] Kim Jongchan, Kim Sungho. Temperature estimation adaptive to variables over distance using infrared–LiDAR[J]. Applied Sciences, 2021, 11(9): 4063. DOI: 10.3390/app11094063
[4] 马刚, 郭杨, 苏婧, 等. 卫星探测通道晴空大气吸收透过率快速计算研究进展[J]. 地球物理学报, 2023, 66(6): 2275-2291. MA Gang, GUO Yang, SU Jing, et al. An overview of fast calculation to clear atmospheric absorption transmittance of satellite channel[J]. Chinese Journal of Geophysics, 2023, 66(6): 2275-2291.
[5] 何光聪, 胡月明, 肖志峰. 基于TM8卫星热红外数据地表温度反演及其应用[J]. 广东农业科学, 2015, 42(6): 146-153, 193. HE Guangcong, HU Yueming, XIAO Zhifeng. Land surface temperature retrieval based on TM8 thermal infrared data and its application[J]. Guangdong Agricultural Sciences, 2015, 42(6): 146-153, 193.
[6] 陈秀红, 魏合理, 吕炜煜, 等. CART软件计算的红外大气透过率和实测值比较[J]. 激光与红外, 2009, 39(4): 403-406. DOI: 10.3969/j.issn.1001-5078.2009.04.015 CHEN Xiuhong, WEI Heli, LYU Weiyu, et al. Comparison of infrared atmospheric transmittance calculated by CART software with measured values[J]. Laser & Infrared, 2009, 39(4): 403-406. DOI: 10.3969/j.issn.1001-5078.2009.04.015
[7] 李晋华, 王志斌, 陈媛媛, 等. 近红外氧气A带大气透过率的计算[J]. 激光与红外, 2013, 43(10): 1142-1145. DOI: 10.3969/j.issn.1001-5078.2013.10.13 LI Jinhua, WANG Zhibin, CHEN Yuanyuan, et al. Calculation of near infrared oxygen A-band atmospheric transmittance[J]. Laser & Infrared, 2013, 43(10): 1142-1145. DOI: 10.3969/j.issn.1001-5078.2013.10.13
[8] 王充, 汪卫华. 红外辐射大气透过率研究综述[J]. 装备环境工程, 2011, 8(4): 73-76. WANG Chong, WANG Weihua. Summarization on atmospheric transmittance of infrared radiation[J]. Equipment Environmental Engineering, 2011, 8(4): 73-76.
[9] 赵军, 张建, 杜翠兰. 红外辐射大气透过率修正函数[J]. 激光与红外, 2006, 36(9): 866-867. DOI: 10.3969/j.issn.1001-5078.2006.09.012 ZHAO Jun, ZHANG Jian, DU Cuilan. Revisable function of IR-radiation atmospheric transmissivity[J]. Laser & Infrared, 2006, 36(9): 866-867. DOI: 10.3969/j.issn.1001-5078.2006.09.012
[10] 雷忠腾. 基于MODIS热红外分裂窗波段反演晴空大气水汽的研究[D]. 青岛: 山东科技大学, 2020. LEI Zhongteng. Atmosphere water vapor retrieval over clear sky on land based on modis thermal infrared splitting-window band[D]. Qingdao: Shandong University of Science and Technology, 2020.
[11] 李云红, 张龙, 王延年. 红外热像仪外场测温的大气透过率二次标定[J]. 光学精密工程, 2010, 18(10): 2143-2149. LI Yunhong, ZHANG Long, WANG Yannian. Second calibration of at mospheric trans mission coefficients on temperature measurement of infrared thermal imager in fields[J]. Optics and Precision Engineering, 2010, 18(10): 2143-2149.
[12] 王文慧, 王高, 张猛. 红外辐射大气透过率影响因素的研究[J]. 激光与红外, 2017, 47(3): 308-312. WANG Wenhui, WANG Gao, ZHANG Meng. Study on influence factors of atmospheric transmittance of infrared radiation[J]. Laser & Infrared, 2017, 47(3): 308-312.
[13] ZENG Hui, REN Huazhong, NIE Jing, et al. Land surface temperature and emissivity retrieval from nighttime middle and thermal infrared images of chinese fengyun-3D MERSI-Ⅱ[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 7724-7733. DOI: 10.1109/JSTARS.2021.3098579
[14] 吴骅, 李秀娟, 李召良, 等. 高光谱热红外遥感: 现状与展望[J]. 遥感学报, 2021, 25(8): 1567-1590. WU Hua, LI Xiujuan, LI Zhaoliang, et al. Hyperspectral thermal infrared remote sensing: current status and perspectives[J]. National Remote Sensing Bulletin, 2021, 25(8): 1567-1590.
[15] 孟鹏, 胡勇, 巩彩兰, 等. 热红外遥感地表温度反演研究现状与发展趋势[J]. 遥感信息, 2012, 27(6): 118-123, 132. MENG Peng, HU Yong, GONG Cailan, et al. Thermal infrared remote sensing of surface temperature inversion: current status and future prospects[J]. Remote Sensing Information, 2012, 27(6): 118-123, 132.
[16] 李召良, 段四波, 唐伯惠, 等. 热红外地表温度遥感反演方法研究进展[J]. 遥感学报, 2016, 20(5): 899-920. LI Zhaoliang, DUAN Sibo, TANG Bohui, et al. Review of methods for land surface temperature derived from thermal infrared remotely sensed data[J]. National Remote Sensing Bulletin, 2016, 20(5): 899-920.
[17] 姚前, 许华, 樊程, 等. 一种应用于卫星中红外遥感的辐射传输简化计算方案[J]. 红外与毫米波学报, 2024, 43(3): 375-390. YAO Qian, XU Hua, FAN Cheng, et al. A simplified parameterization scheme of mid-infrared radiative transfer for satellite remote sensing[J]. Journal of Infrared And Millimeter Waves, 2024, 43(3): 375-390.
-
期刊类型引用(4)
1. 吕伽奇，丁帅，庞静珠，许小进. 基于改进LeNet-5网络的堆芯燃料组件编码识别. 东华大学学报(自然科学版). 2024(02): 121-128 . 
百度学术
2. 毛羽，郑怀华，李隆，张傲. 基于热红外图像的光伏板热斑检测方法研究. 自动化仪表. 2024(05): 25-29+34 . 百度学术
3. 王晓君，孙梓林，王雁. 基于AMP架构的青霉素结晶与发酵检测系统设计. 仪表技术与传感器. 2024(05): 66-73 . 百度学术
4. 赵兴文. 机器学习在信用贷款评分中的应用. 福建电脑. 2023(02): 31-34 . 百度学术
其他类型引用(15)
计量
- 文章访问数: 21
- HTML全文浏览量: 4
- PDF下载量: 6
- 被引次数: 19
