引用本文
叶凯迪, 秦敏, 方武, 段俊, 唐科, 孟凡昊, 张鹤露, 谢品华, 徐文斌. 深紫外波段苯的差分光学吸收光谱DOAS定量方法研究[J]. 光谱学与光谱分析, 2021,41(10): 3007-3013.
YE Kai-di, QIN Min, FANG Wu, DUAN Jun, TANG Ke, MENG Fan-hao, ZHANG He-lu, XIE Pin-hua, XU Wen-bin. Study on the Measurement of C6H6 in the 195~208 nm Band Based on Differential Absorption Spectroscopy [J]. Spectroscopy and Spectral Analysis, 2021,41(10): 3007-3013.
Doi:10.3964/j.issn.1000-0593(2021)10-3007-07  
Permissions
Copyright©2021, 《光谱学与光谱分析》期刊社
《光谱学与光谱分析》期刊社 所有
深紫外波段苯的差分光学吸收光谱DOAS定量方法研究
叶凯迪1,2, 秦敏1,*, 方武1, 段俊1, 唐科1,2, 孟凡昊1,2, 张鹤露1,2, 谢品华1,2,3, 徐文斌4
1.中国科学院安徽光学精密机械研究所, 中国科学院环境光学与技术重点实验室, 安徽 合肥 230031
2.中国科学技术大学科学岛分院, 安徽 合肥 230026
3.中国科学院区域大气环境研究卓越创新中心, 中国科学院城市环境研究所, 福建 厦门 361021
4.光学辐射重点实验室, 北京 100854
*通讯作者 e-mail: mqin@aiofm.ac.cn

作者简介: 叶凯迪, 1994年生, 中国科学院安徽光学精密机械研究所硕士研究生 e-mail: kdye@aiofm.ac.cn

收稿日期: 2020-03-04 修回日期: 2020-07-21
基金: 国家重点研发计划项目(2017YFC0209900, 2016YFC0201002), 中国科学院重点部署项目(KFZD-SW-320), 安徽省科技重大专项(16030801120)资助
摘要

取代基通过取代苯环上的H原子形成不同苯系物(苯、 甲苯、 二甲苯等), 其共有结构苯环上的不固定π键电子受到激发, 使得苯系物在紫外波段240~280 nm具有明显的特征吸收结构, 鉴于此大气中的苯及相关的苯系物可以通过差分光学吸收光谱(DOAS)方法来进行定量, 但采用该波段测量需要考虑以下问题: 首先是氧气(O2)的吸收干扰问题, 苯(C6H6)在该波段的吸收截面与O2在243~287 nm Herzzberg带相互重叠, 且O2的特征光谱结构随O2的浓度不同而变化, 导致O2的吸收光学密度与O2的浓度不成线性关系。 其次, 苯系物结构上的相似性使其在紫外波段的特征吸收结构差别较小并且相互重叠, 从而对C6H6的拟合产生干扰。 此外, 除了O2和苯系物以外, 还有臭氧(O3)、 二氧化硫(SO2)等干扰。 C6H6在195~208 nm的深紫外波段具有较大的吸收截面(2.417×10-17 cm2·molecule-1), 为240~260 nm处截面大小(2.6×10-18 cm2·molecule-1)的9倍左右, 针对C6H6在深紫外195~208 nm波段的吸收特征, 开展便携式DOAS定量方法研究, 采用该波段进行C6H6的光谱定量分析并应用到实际的外场观测。 通过建立C6H6与干扰气体SO2, 氨(NH3), 二硫化碳(CS2)和一氧化氮(NO)的差分吸收截面的二维相关性矩阵, 获取C6H6光谱定量的最优反演波段。 通过开展实验室条件下C6H6, SO2和NH3不同浓度配比的混气实验对195~208 nm波段反演C6H6的效果进行评估。 实验结果显示, 采用195~208 nm波段进行光谱反演的探测限为17.6 μg·m-3, 光谱反演浓度与理论浓度的相对测量误差均小于5%且RSD(相对标准偏差)小于3%, 同时与240~260 nm波段反演结果进行对比, 相对误差小于5%。 在外场实际情况下, 利用便携式DOAS系统获取190~300 nm的大气测量光谱, 通过DOAS方法解析并结合GPS信息, 获得了某化工园区C6H6的污染浓度分布, 实验结果表明采用195~208 nm深紫外波段同样能适用于对C6H6的光谱定量, 与240~260 nm波段反演结果进行对比, 二者的相关性达到了0.98且相对误差小于10%。

关键词: 差分光学吸收光谱; DOAS; 苯; 深紫外
中图分类号:O433.4 文献标志码:A
Study on the Measurement of C6H6 in the 195~208 nm Band Based on Differential Absorption Spectroscopy
YE Kai-di1,2, QIN Min1,*, FANG Wu1, DUAN Jun1, TANG Ke1,2, MENG Fan-hao1,2, ZHANG He-lu1,2, XIE Pin-hua1,2,3, XU Wen-bin4
1. Key Laboratory of Environment Optics and Technology, Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Hefei 230031, China
2. Science Island Branch of Graduate School, University of Science and Technology of China, Hefei 230026, China
3. CAS Center for Excellence in Regional Atmospheric Environment, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China
4. Science and Technology on Optical Radiation Laboratory, Beijing 100854, China
*Corresponding author
Abstract

Substituents form BTX (benzene, toluene, xylene, etc) by replacing the H atom on the benzene ring. The unfixed π-bond electrons on the benzene ring of the common structure are stimulated, which makes BTX have a distinct characteristic absorption structure in the ultraviolet band of 240~280 nm. BTX in this atmosphere can be quantified by the differential optical absorption spectroscopy (DOAS) method. However, problems need to be taken into consideration when adopting the measurement of this band, such as the non-linear absorption of O2, cross interference between BTX due to the existence of similar absorption structures, and overlapping absorption interference of other gases. Benzene has a large absorption cross-section (2.417×10-17 cm2·molecule-1) in the deep ultraviolet band of 195~208 nm, which is about nine times the cross-section size (2.6×10-18 cm2·molecule-1) at 240~260 nm. Therefore, according to the characteristics of benzene in the deep ultraviolet 195~208 nm band, a portable differential optical absorption spectrum (DOAS) quantitative method is studied. This band is used to carry out the quantitative spectral analysis of benzene for the first time and applied to the actual field observation. The optimal retrieval band of benzene spectral quantification is obtained by establishing the two-dimensional correlation matrix of the differential absorption cross sections of benzene and the interfering gases SO2, NH3, CS2, and NO. Furthermore, the effect of retrieval of benzene in 195~208 nm band was evaluated by carrying out mixed gas experiments with different concentrations of benzene, SO2 and NH3 under laboratory conditions. The experimental results show that the detection limit of the spectral retrieval using the 195~208 nm band reaches 17.6 μg·m-3, the relative measurement error between the measured concentration of the spectral retrieval and the configured concentration is less than 5% and the relative standard deviation (RSD) is less than 3%. Compared with the retrieval results in the 240~260 nm band, the relative error is less than 5%. In the actual field condition, a portable differential optical absorption spectroscopy (DOAS) system was used to obtain the atmospheric measurement spectrum of 190~300 nm, and the DOAS method was applied to analyze and combine with the GPS information to obtain the pollution concentration distribution of benzene in a chemical industrial park. The experimental results show that the deep ultraviolet band of 195~208 nm can also be adapted to the quantitative analysis of benzene. And compared with the retrieval results of the 240~260 nm band, the correlation coefficient of the two retrieval bands has reached 0.98 and the relative error is less than 10%.

Keyword: Differential optical absorption spectroscopy; DOAS; Benzene; Deep ultraviolet
引言

苯(C6H6)是苯系物(BTX)的重要成分之一, 广泛参与大气中的光化学反应。 C6H6在大气中主要通过与OH自由基反应促使光化学烟雾的形成[1]; 同时C6H6也具有一定的毒性, C6H6在1993年即被世界卫生组织( WHO) 确定为致癌物, 中国于2007年也将C6H6列为致癌物质[2]。 因此, 开展C6H6的高灵敏准确定量测量对研究其化学行为过程至关重要。

差分光学吸收光谱(DOAS)技术由Platt等于20世纪70年代末提出[3], 由于具有多组分同步测量、 测量精度高、 可实时在线等优点, 被广泛应用于大气痕量污染物的监测, 大气中苯系物的DOAS测量始于20世纪90年代[4, 5]。 目前DOAS技术主要通过测量C6H6在240~260 nm处的特征吸收来定量C6H6的浓度, 但采用该波段测量需要考虑苯系物之间由于相似的特征吸收结构存在的交叉干扰、 O2的非线性吸收以及其他气体的吸收重叠等问题[6, 7]。 在采用该波段进行苯系物的光谱解析时, 要求光谱仪的最优光学分辨率需满足0.15~0.16 nm[8]。 Dawes等[9]于2017年通过真空紫外光谱技术获得了115~330 nm的C6H6的吸收截面。 相比于DOAS方法通常采用的240~260 nm测量波段, 190~210 nm波段具有更大的吸收截面, 且O2和O3等干扰气体在190~210 nm波段没有窄带吸收特征, 同时对光谱仪的光学分辨率也不如传统测量波段要求严苛。

本文利用便携式DOAS系统开展对C6H6在深紫外波段的定量方法研究, 通过建立二维相关性矩阵选择195~208 nm作为C6H6的反演波段, 通过混气实验评估了195~208 nm波段对C6H6的浓度进行定量测量的可行性, 计算了系统采用这一波段进行光谱反演的检测限。 在此基础上, 利用便携式DOAS系统搭载车载平台, 结合GPS信息, 应用于外场观测, 获取了某化工园区内的C6H6浓度分布。

1 DOAS技术
1.1 DOAS原理

DOAS的基本原理是基于Lambert-Beer定律[3], 光源发出的光经过一定光程后, 会存在许多消光因素的影响, 包括气体分子的瑞利散射[ε R(λ )]、 气溶胶米散射[ε M(λ )]等, 接收端的光强衰减为I(λ )

I(λ)=I0(λ)exp-L∑i=1n(σi(λ)ci)+εR(λ)+εM(λ)(1)

式(1)中: σi(λ )和ci分别代表第i种气体的吸收截面和浓度; L为光程; ε R(λ )和ε M(λ )分别为瑞利散射和米散射造成的光强衰减。 差分吸收光谱技术的思想是将吸收截面分为随波长做快速变化和缓慢变化的两部分, 通过数字滤波将与气体无特征关系的慢变化部分去除得到差分光学密度OD(differential optical density)。

OD=lnI'0(λ)I(λ)=L∑i=1n(σ'i(λ)ci)(2)

最后将OD与各个待测成分的差分吸收截面通过非线性最小二乘拟合得到痕量气体浓度。

1.2 测量系统

用于实验室样气以及外场测量的自研便携式DOAS系统[10]由氘灯光源、 准直镜头、 怀特(White)型多次反射池和光谱仪等单元组成, 图1为测量系统示意图。 系统中的多次反射池基长为39.6 cm, 采用三块严格等焦距的凹球面反射镜共轭放置, 反射镜表面镀紫外增强铝膜, 使其对紫外190~310 nm波段的反射率大于85%, 通过入射光在多反池内部来回反射增加吸收光程。 氘灯光源(185~400 nm)发射的紫外光经准直镜头准直后进入多次反射池, 光束在反射池中经过28次反射后最终在多次反射池的出射口出射, 总光程为11.1 m。 出射光通过汇聚系统耦合至光纤, 经光纤传输至光谱仪(Maya2000 Pro, 分辨率约为0.3 nm), 光谱仪将测量的光谱传输到电脑, 最终通过DOAS算法对光谱数据分析得到气体浓度信息。

图1 测量系统Fig.1 Measurement system

实验室实验中, 在多次反射池的出射光束口和光纤汇聚端的光路上加入一个长度为3 cm的样品池(通光面为JGS1石英材质), 使用多组分动态配气仪(ZTD-003, 中国计量科学研究院)配置不同比例浓度的混合气体(标气浓度C6H6: 2 552 mg· m-3, SO2: 5 497 mg· m-3, NH3: 7 300 mg· m-3, 液化空气气体科技有限公司)通入样品池中, 模拟样气进入多次反射池中的吸收情况, 配气系统如图1中虚线框图所示。 为避免混合样气进入多次反射池干扰测试结果, 实验过程中向多次反射池持续通入高纯氮气; 外场观测中则移除样品池, 通过泵抽取外界大气进入多次反射池内进行测量。

2 C6H6反演波段的选取
2.1 C6H6的吸收截面及干扰气体成分的吸收重叠影响

根据Lambert-Beer定律, 要获得待测气体的浓度信息, 需要已知该气体分子的吸收截面。 190~300 nm处苯系物以及部分干扰气体的吸收截面[4]如图2a—j所示, 此部分讨论的气体成分是在该波段存在特征窄带吸收的成分, 对于O2, 尽管在195~241 nm波段也存在吸收会造成短波紫外的光强衰减[11], 但由于没有明显的窄带吸收特征, DOAS方法可通过高通滤波去除此类气体的影响。 在190~210 nm处除了C6H6的吸收外, 还存在SO2和NH3等气体成分的吸收。 在240~260 nm处由于苯系物共有结构苯环上的不固定π 键电子受到激发使得苯系物具有相似吸收结构, 除了C6H6, 甲苯(Toluene)、 二甲苯(M-xylene, P-xylene, O-xylene)等苯系物的吸收也集中在这一波段; 此外O2和O3等成分由于在一些特定环境大气中的含量较高, 通常成为这一波段主要的吸收结构而干扰包括C6H6在内的苯系物的准确定量, 因此在光谱反演中需要特别注意。

图2 190~300 nm波段原始吸收截面(a—j)以及C6H6的差分吸收截面(k)Fig.2 Original absorption cross section at 190~300 nm (a—j) and differential absorption cross section (k) of C6H6

C6H6在298 K下从190~300 nm波段的吸收截面[9]如图2中a所示。 根据DOAS算法将C6H6的原始吸收截面与仪器函数卷积获得与实验测量系统相同分辨率下的吸收截面, 通过高通数字滤波, 得到C6H6的差分吸收截面, 如图2中k所示, 其主要的吸收波段集中在190~210 nm波段和240~260 nm波段。 通过计算得到, 190~210 nm处的差分截面大小为8.628× 10-8(峰峰值), 约为240~260 nm处(240~260 nm处的截面大小为2.181× 10-8)的4倍, 采用190~210 nm光谱波段对苯进行定量解析时需要考虑扣除同一波段内其他痕量气体的吸收。

2.2 反演波段的选择

190~210 nm波段反演C6H6的干扰因素主要来自于拟合波段内其他气体的吸收, 如图2中b—d所示, SO2和NH3等气体的吸收与C6H6的吸收在该波段有部分重叠。 为了优化苯的反演波段, 避免190~210 nm波段内其他气体(SO2, NH3等)对C6H6光谱准确定量的干扰, 对该波段内C6H6与其他干扰气体的吸收截面的相关性进行研究。 通过建立C6H6与其他气体的差分吸收截面的二维相关性矩阵, 评估在同一波段不同气体吸收截面相互间的相关性, 如两者之间相关性较弱, 可认为对C6H6反演的干扰较小。 如图3所示, 横坐标为起始波段, 纵坐标为截止波段, 图中每一个点都对应一个反演波段。

图3 C6H6与其他干扰气体吸收截面的相关性矩阵Fig.3 Correlation matrices between C6H6 and other interference gas absorption cross sections

C6H6和SO2, CS2的相关性均接近0, 但是C6H6和NH3, NO在部分区域有一定相关性。 图3(a)中C6H6与NO在起始波长196~199 nm, 截止波长206~210 nm波段(区域1), 图3(c)中C6H6与NH3在起始波长196~199 nm, 截止波长206~210 nm波段(区域2)两个区域具有一定相关性, 存在一定的交叉干扰。 因此在选择C6H6的反演波段时需要避开这两个区域。

综合分析, 考虑各种干扰气体的影响, 在截面相关系数小的波段范围内, 选择相对较宽的波段作为反演波段, 最终采用195~208 nm作为C6H6的反演波段。

3 苯的实验室样气测量

为了评估在多种干扰气体环境下195~208 nm波段反演C6H6的效果, 使用配气系统(ZTD-003)在实验室进行混气实验。 根据实际应用环境不同[10, 12], 模拟城区交通环境以及化工园区环境。

由于短波紫外光源的稳定性和探测器的响应、 镜面反射率等因素的影响, 为获取实验的最佳测试参数, 实验首先将高纯度氮气经过PFA管导入样品池和多次反射池中进行吹扫, 采集无被测气体吸收的背景光谱。 为保证195~208和240~260 nm能够同时获得足够的光强信号, 选择80 ms的光谱积分时间(在240~260 nm不饱和, 195~208 nm波段获得相对较大的光强情况下的最大积分时间)。 为提高光谱信噪比, 测量光谱采用500次叠加后平均作为单条光谱进行后续光谱分析。 将参考光谱扣除背景后, 如图4(a)所示, 195~208 nm处的平均光强为11 400 counts, 而240~260 nm处的光强达到了50 500 counts, 两个波段的平均光强相差约4倍, 相对而言, 后者的信噪比要高于前者。 采用相同的积分时间依次测量两条灯谱, 分别扣除背景后将两条谱相除, 作高通, 处理结果如图4(b)。 根据下式确定两个波段所对应的最低检测浓度cmin

cmin=τResσ(λ)LN(3)

式(3)中, σ(λ )为1 ppb(10-9)痕量气体的差分吸收截面, L为光程(11.1 m), τ Res为剩余结构光学密度, 它体现了干扰因素对系统的影响, 由图4(b)得到195~208 nm波段的τ Res为1.06‰ , 240~260 nm波段的τ Res为0.51‰ 。 N为实际光谱拟合范围内苯吸收峰的个数。 通过计算得到195~208 nm波段C6H6的最低检测限约为17.6 μg· m-3, 240~260 nm的C6H6的最低检测限48.8 μg· m-3。 通过对比195~208 nm波段与240~260 nm波段的计算结果, 195~208 nm波段对C6H6具有更低的检测下限。

图4 参考光谱(a)以及剩余噪声(b)Fig.4 Reference spectrum (a) and residual noise (b)

由于系统测量时的光程为11.1 m, 为样品池长度的370倍, 因此, 样气理论浓度为系统配置浓度的1/370。 分别配置4:1:1, 3:1:1, 2:1:1三个不同浓度比例的C6H6, SO2和NH3的混合气体(C6H6和NH3的系统配置浓度为37 mg· m-3, 对应的理论浓度为100 μg· m-3, SO2的系统配置浓度分别为148, 111和74 mg· m-3, 对应的理论浓度分别为400, 300和200 μg· m-3, NH3由于气路管线以及样品池腔壁的吸附损耗导致反演浓度比实际配气浓度低)。 图5展示了混气实验期间(SO2:C6H6:NH3=2:1:1)测量谱经过处理后得到的195~210 nm波段差分吸收光谱示例。 从图5中可以看出, 经过高通滤波处理后, 能够较好的去除光谱中的慢变化部分, 最终通过最小二乘拟合得到C6H6, SO2和NH3三种气体的浓度。

图5 混气实验测量的差分吸收光谱示例
配气浓度比SO2:C6H6:NH3=2:1:1Fig.5 Example of differential absorption spectrum measured in mixed gas experiment
Ratio of gas concentration SO2:C6H6:NH3=2:1:1

在通入气体浓度稳定后连续测量10分钟, 实验结果如图6所示, SO2的平均浓度分别为399.8, 306.4和205.6 μg· m-3, NH3由于吸附损耗测量得到的实际平均浓度为69 μg· m-3。

图6 混合气体的反演结果Fig.6 Retrieval results of mixed gas

对三种配气比例下的C6H6的理论浓度与反演浓度的平均值进行对比, 计算其相对标准偏差(RSD)以及相对误差。 由表1可知, 在不同比例浓度的SO2和NH3的影响下, C6H6的反演浓度基本保持稳定, RSD小于3%, 且与理论浓度的相对误差均在5%以内。 同时使用C6H6在240~260 nm波段的截面对测量光谱进行反演, 反演结果与前者进行对比分析, 两波段反演误差均在5%以内, 说明两个波段的反演结果具有较好的一致性。

表1 C6H6的测量结果 Table 1 Measurement result of C6H6
4 外场实验及实验结果

选择天津市某化工园区, 将便携DOAS系统放置于车内开展走航观测实验。 实验过程中, 在仪器的进气处使用过滤器进行过滤, 以防止多次反射池被污染。 同时, 为减少车辆排放对实验的影响, 将采样管口放置于车顶上方约10 cm处的位置。 实验期间系统实际采样抽气率约为20 L· min-1; 车辆行驶速度约为20 km· h-1; 实验过程中由车载系统配备的GPS接收机记录经纬度信息。

图7展示了195~208 nm以及240~260 nm波段反演C6H6的光谱示例(2019年1月14日外场实验期间测量得到的大气谱), 考虑园区的生产工艺以及典型污染物排放特征, 在光谱反演中考虑SO2和NH3等作为干扰气体成分。图7中的(a)和(f)是195~208 nm以及240~260 nm波段对应的测量光谱。 (b), (c), (d)和(g), (h)是测量光谱经过高通滤波后的差分吸收截面(灰线)和拟合得到的差分吸收结构(黑线)的叠加, 拟合的波段范围分别是195~208 nm以及240~260 nm, 从图中可以看出SO2, NH3和C6H6的拟合结果较好。 (e)和(i)是测量光谱扣除了气体吸收结构后的拟合残差(黑线), 195~208 nm波段拟合残差的峰峰值为3× 10-3(240~260 nm波段拟合残差的峰峰值为1× 10-3), 主要来源为未知成分的吸收和系统噪声。 由于实际大气测量时, 短波紫外波段易受到瑞利散射和O2的吸收造成短波紫外的衰减[11], 且光谱的信噪比会受到仪器本身和测量环境的影响, 此时根据式(3)计算出两个波段的检测限分别为52.7和96.6 μg· m-3, 其稍差于实验室测试条件下的检出限。

图7 195~208 nm波段以及240~260 nm波段C6H6的拟合示例
(2019年1月14日, 13:28测量的光谱)Fig.7 Example for C6H6 fitting at 195~208 and 240~260 nm
(Spectrum measured on January 14, 2019, 13:28)

图8为天津市某化工园区内走航路径上C6H6的浓度分布情况(环境中同时存在SO2和NH3等气体)以及195~208与240~260 nm波段反演结果的相关性分析, 红色圆圈为化工园区所在位置, 红色箭头代表了观测期间的主导风向。 C6H6在整个实验期间的平均浓度为65 μg· m-3, 在园区东南侧出现高浓度的C6H6, 最高达到了392 μg· m-3, 在主导风向的传输下, 污染气团有可能对下游区域的空气质量造成影响。

图8 某化工园区C6H6的浓度分布Fig.8 Concentration distribution of C6H6in a chemical park

5 结论

利用便携式DOAS系统开展了深紫外波段苯的定量方法研究, 通过建立吸收截面间的二维相关性矩阵, 选择交叉干扰最小的波段, 确定了195~208 nm作为苯的反演波段。 结合车载GPS数据, 实现了某化工园区苯浓度的快速获取。 本文的研究结果表明, 在保证光谱信噪比的基础上, 可通过对195~208 nm波段的光谱解析实现对苯的精确测量, 与240~260 nm波段相比具有更低的检测限, 为苯的光谱定量提供一种新的解决思路。

参考文献
[1] Stockwell W R, Middleton P, Chang J S. Journal of Geophysical Research, 1990, 95(D10): 16343. [本文引用:1]
[2] PEI Bing-an(裴炳安). Petroleum Refinery Engineering(炼油技术与工程), 2012, 42(11): 62. [本文引用:1]
[3] Platt U, Perner D, Pat H. Journal of Geophysical Research, 1979, 84(10): 6329. [本文引用:2]
[4] Lofgren L. International Journal of Environmental Analytical Chemistry, 1992, 47(1): 69. [本文引用:2]
[5] Schaefer K, Emeis S M, Stockhause M. Environmental Monitoring and Remediation Technologies, 1999, 3534: 212. [本文引用:1]
[6] Volkamer R, Etzkorn T, Geyer A. Atmospheric Environment, 1998, 32(21): 3731. [本文引用:1]
[7] Trost B, Stutz J, Platt U. Atmospheric Environment, 1997, 31(23): 3999. [本文引用:1]
[8] Peng Fumin, Xie Pinhua, Zhang Yinghua, et al. Journal of Environmental Sciences, 2008, 20(5): 632. [本文引用:1]
[9] Dawes A, Pascual N, Hoffmann S V, et al. Physical Chemistry Chemical Physics, 2017, 19(40): 27544. [本文引用:2]
[10] SHEN Lan-lan, QIN Min, SUN Wei, et al(沈兰兰, 秦敏, 孙伟, 等). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2016, 36(6): 1936. [本文引用:2]
[11] Cheung S C, Yoshino K, Parkinson W H. Planetary and Space Science, 1987, 35: 397. [本文引用:2]
[12] TAO Shuang-cheng, DENG Shun-xi, HAO Yan-zhao, et al(陶双成, 邓顺熙, 郝艳召, 等). China Environmental Science(中国环境科学), 2019, 39(2): 542. [本文引用:1]