2. 吉林大学地球科学学院, 长春 130061;
3. 中兴通讯股份有限公司, 西安 710065
2. College of Earth Sciences, Jilin University, Changchun 130061, China;
3. Zhongxing Telecommunication Equipment Corporation, Xi'an 710065, China
0 前言
2013年12月初, 长征三号乙运载火箭携带中国第一艘月球车从西昌卫星发射中心发射, 实现了中国首次月面软着陆。嫦娥三号(CE3)由着陆器和巡视探测器(“玉兔号”月球车)组成, 其主要任务是进行首次月球软着陆和自动巡视勘察, 获取月球内部的物质成分并进行分析, 将一期工程的“表面探测”引申至内部探测。其中“玉兔号”月球车携带了高能所粒子天体物理中心研制的粒子激发X射线谱仪(APXS), 它是机械臂上唯一的载荷。APXS通过放射源主动激发月岩或月壤中的元素, 使其产生特征X射线, 从而获得月岩和月壤中的化学成分及含量, 其可为月球地质化学过程和形成演化研究提供重要依据[1]。2013年12月25日在机械臂投放过程中, 粒子激发X射线谱仪采用自身的距离感知功能, 成功指导机械臂将探头精确投放到距离月面2~3 cm处, 完成了首次月面元素的原位分析[2]。本文对CE3 “玉兔号”月球车传输回来的APXS数据[3]进行了特征X射线谱数据处理和分析, 以期获得月壤或岩石的化学成分, 为着陆区的地质演化提供技术支持。
1 粒子激发X射线谱线的获得和处理 1.1 特征X射线能谱数据说明CE3玉兔号月球车软着陆后, 共实行了两次APXS测量[1, 3], 第一次于2013年12月25日将数据传回地面应用系统[4-5]。经处理后得到两组数据文件, 分别为“CE3_BMYK_PIXS-E_SCI_N_20131221124501_ 20131223174500_0005_A.2B”和“CE3_BMYK_PIXS-E_SCI_N_20131223174501_ 20131226000000_0006_A.2B”, 均为地面应用系统发布的PDS格式2B级数据。其中, 0005_A测量时间共12 min16 s左右, 0006 _A测量时间共3 h43 min05 s左右。
数据文件每条记录的格式为:
测量时间
工作模式
月球车机械臂的方位角
月球车机械臂的俯视角
月球车机械臂腕的俯视角
温度060
温度001
探测器的俯视位置
探测器的翻滚位置
探测器的偏航位置
特征X射线谱数据第1道
特征X射线谱数据第2道
⋮
特征X射线谱数据第2 048道
仪器状态参数
其中, 谱数据的第1道和最后道数据是无效数据。
1.2 特征X射线谱数据处理X射线能谱仪APXS每8 s记录一次测量结果, 由于时间非常短, 计数率较低, 为了获得较高的分析精度, 采用累加谱记录构成测量点的谱曲线进行X射线能谱分析[3-4]。首先得到CE3 0006测量点的累积X射线能谱曲线(图 1a), 然后归一化为2 h的计数率。从图 1a的谱曲线可以看出, 曲线上毛刺较多, 主要原因是放射性统计涨落, 以及噪声的影响。为了消除这种影响, 采用小波去噪[6]的方法进行了处理, 经处理后得到图 1b。从图 1b中可看出, 去噪后曲线更圆滑, 特征峰更清晰, 这对进一步从谱线特征峰识别元素及分析提供了更高的质量保证。
|
| a.原始测量能谱曲线; b.经小波滤波去噪后的能谱曲线。 图 1 CE3 0006测点累积X射线能谱曲线 Figure 1 Accumulation X-ray energy spectrum of CE3 No.0006 measurement |
|
|
由CE3月球车载的X射线能谱仪APXS所给定的标定能量公式[5]为
(1) 其中:
式中:E为能量; Ch为道号; a1=10.445;b1=7.43×10-3; c1=1.21×10-4; d1=-8.60×10-7; e1=-4.43×10-8; a2=183.36;b2=-0.74;c2=3.27×10-3; d2=-6.12×10-4; T为温度。
式(1)对0006_A文件的累计X射线谱进行了能量标定, 结果见图 2。
|
| ( )内数字为峰编号。 图 2 CE3 0006测量点X射线谱特征元素及能量 Figure 2 Energy and element of X-ray energy spectrum characteristic peak |
|
|
1913年, 莫塞莱(H.G.J.Moseley)发现元素X射线光谱频率的平方根与原子序数Z之间存在着以下简单的线性关系[7]:
(2) 式中:v是X射线的频率; Z为原子序数; a和b均为常数。这就是莫塞莱定律。
由于在轨道上以一定速度运动的电子向心力是电子在核库仑场中所受的静电引力, 而且若系统保持稳定, 则电子运动产生的离心力与电子在核库仑场中所受静电引力必须相等, 与此同时, 还要考虑到在原子中其他核外电子对核电荷有屏蔽作用; 因此第n轨道电子的能态En可表示为
(3) 式中:R为里德伯常数, R=1.097×107 m-1; h为普朗克常数, h=6.625 59×10-34 J·s; C为光速, C=3×108 m/s; n为主量子数; σ为考虑电子对核电荷的屏蔽作用引入屏蔽系数, 其基本与Z无关。
当电子由量子数为n1的能态跃迁到主量子数为n2的能态, 其能量关系为
(4) 式中, ΔE为主量子数n1能态到n2能态的能量差。
式(4)表明, 不同的元素具有不同的能级差, 所以能发射出不同的特征X射线能量, 同一线系的特征X射线能量随原子序数的增大而增加; 所以可以通过不同核素发射出的不同特征X射线能量来识别物质的构成元素。这就是X荧光测量能对物质进行元素成分分析的基本原理。有关不同元素的X射线特征峰能量值见文献[7]附录5。
3.2 月表土壤特征X射线谱的元素识别月球车X射线能谱仪APXS的测量点的坐标为(-19.508 827, 44.119 717), 位于雨海西北扩展开来的一个玄武岩熔岩平原[4], 月海的玄武岩[8-11]主要由辉石、长石、橄榄石和不透明矿物组成, 其中不透明矿物主要是钛铁矿。月球车X射线能谱仪的激发源由4个55Fe(半衰期2.73 a, 活度4×(259.00×107 Bq))和4个109Gd(半衰期1.27 a, 活度4×(9.25×107 Bq))组成, 以Si漂移型半导体为探测器, 2048道幅度分析器, 仪器的测量能量范围为0.39~22.1 keV[3]。
从图 2各峰所对应的能量, 我们可以找到相应被激发元素的特征X射线能量, 从而得到相应能量特征峰可能对应的元素。由于55Fe激发源的能量为5.89 keV和6.49 keV, 而109Gd激发源的能量为22.16 keV和24.94 keV[4], 可激发的元素为Na—Mo, 所以能量高于18 keV的谱峰不进行元素确定。根据元素的特征能量, 本文确定了各峰对应的元素, 将其结果列于表 1中。
| 测点X射线能谱特征 | 相对应元素及理论特征峰能量 | 判别核素 | ||||
| 峰道号 | 峰编号 | E/keV | 元素 | E/keV | ||
| 68 | 1 | 1.251 | Mg | 1.254(Kα) | Mg | |
| 89 | 2 | 1.477 | Al | 1.478(Kα) | Al | |
| 113 | 3 | 1.735 | Si | 1.740(Kα) | Si | |
| 133 | 4 | 1.951 | P | 2.015(Kα) | P | |
| Ca逃逸峰 | ||||||
| 164 | 5 | 2.284 | S | 2.308(Kα) | S | |
| Ca逃逸峰 | ||||||
| 208 | 6 | 2.757 | Ti逃逸峰 | 逃逸峰 | ||
| 259 | 7 | 3.300 | K | 3.313(Kα) | K | |
| 295 | 8 | 3.694 | Ca | 3.689(Kα) | Ca | |
| 324 | 9 | 4.005 | Ca | 4.012(Kβ) | Ca | |
| 370 | 10 | 4.500 | Ti | 4.507(Kα) | Ti | |
| 410 | 11 | 4.930 | Ti | 4.931(Kβ) | Ti | |
| 454 | 12 | 5.400 | Cr | 5.410(Kα) | Cr | |
| 500 | 13 | 5.899 | Mn | 5.893(Kα) | Mn(55Fe激发源能量所致) | |
| 547 | 14 | 6.405 | Fe | 6.401(Kα) | Fe和Mn(55Fe激发源能量所致) | |
| Mn | 6.490(Kβ) | |||||
| 607 | 15 | 7.050 | Fe | 7.057(Kβ) | Fe | |
| 646 | 16 | 7.470 | Ni | 7.469(Kα) | Ni | |
| 698 | 17 | 8.029 | Cu | 8.037(Kα) | Cu(仪器探头校样的外壳材料) | |
| 778 | 18 | 8.890 | Cu | 8.904(Kβ) | Cu(仪器探头校样的外壳材料) | |
| 1261 | 19 | 14.089 | Sr | 14.130(Kα) | Sr | |
| 1334 | 20 | 14.872 | Y | 14.919(Kα) | Y | |
| 1408 | 21 | 15.668 | Zr | 15.774(Kα) | Zr | |
| 1495 | 22 | 16.604 | Nb | 16.567(Kα) | Nb | |
| 1587 | 23 | 17.590 | Mo | 17.425(Kα) | Mo | |
| 注:仪器能量分辨率为135 eV(对于5.9 keV)。Kα为K系α特征X射线, Kβ为K系β特征X射线。 | ||||||
由表 1可知, 测点共有23个特征峰, 各特征峰分别对应相应的元素, 将各元素的符号标于相应的特征峰位置上, 见图 2。其中Cu的特征峰(17峰、18峰)是由于仪器的灵敏探头的校准样的外壳所致[3], 而Mn的特征峰(13峰、14峰)是激发源55Fe的散射影响所致; 所以APXS测量到的元素有Si、Ca、Al、Mg、K、P、S、Fe、Ti、Ni、Cr、Mo、Sr、Y、Zr、Nb, 共16种。这一结果与早期人们的研究是相似的[3-5], 其中Fe、Ti、Ca、Si计数率较高, 而K、Mg、Al计数率较低。这些结果与人们研究的有关月海玄武岩具有高的FeO质量分数(>16%), 低或中等的Mg质量分数, 低Al2O3, 变化范围大的TiO2质量分数以及贫Na和K的规律一致[12]; Ni、Cr、Sr、Y、Zr、Nb、Mo计数率较低, 预示着它们是以微量元素的形式存在。P元素的出现可能反映了克里普岩(KREEP岩)的混杂。这与文献[4]中的认识“嫦娥三号着陆区月壤高Fe和Ti、低Al的特征, 说明下覆的玄武岩是一种新的类型。次要元素K和微量元素Zr、Y和Nb的质量分数, 以及它们之间的比值关系, 表示该玄武岩可能混入了10%~20%的克里普组分”[4]相一致, 而且符合“克里普岩极为富含稀土元素、U、Th、K和P”[9]。
4 结论1) 通过对CE3月球车APXS特征X射线的能谱曲线分析, 获得了月岩或月壤的元素包括Si、Ca、Al、Mg、K、P、S、Fe、Ti、Ni、Cr、Mo、Sr、Y、Zr、Nb。
2) 从CE3月球车X射线能谱仪的特征X射线能谱曲线上还可以看出月壤Fe、Ti、Ca、Si、元素特征X射线计数率较高, 预示着这些元素的含量较高, 这与测点的玄武岩特征是一致的; 但P元素的存在, 反映有克里普岩层的混杂。
| [1] | Chen Dengyi, Wu Jian, Hu Yiming, et al. Design and Experimental Verification of Chang'e-3 Moon-Night Survival Device for APXS[J]. Chinese Astronomy and Astrophysics, 2016, 40: 414-424. DOI:10.1016/j.chinastron.2016.07.001 |
| [2] |
梁晓华, 吴明烨, 王焕玉, 等. 嫦娥三号粒子激发X射线谱仪红外距离感知方[J].
光谱学与光谱分析, 2013, 33(5): 1360-1363.
Liang Xiaohua, Wu Mingye, Wang Huanyu, et al. Near Infrared Distance Method for Chang'e-3 Alpha Particle X-Ray Spectrometer[J]. Spectroscopy and Spectral Analysis, 2013, 33(5): 1360-1363. |
| [3] | Fu Xiaohui, Li Chunlai, Zhang Guangliang, et al. Data Processing for the Active Particle-Induced X-Ray Spectrometer and Initial Scientific Results from Chang'e-3 Mission[J]. Research in Astronomy and Astrophysics, 2014, 14(12): 1595-1606. DOI:10.1088/1674-4527/14/12/008 |
| [4] | Zhang Jinhai, Yang Wei, Hu Sen, et al. Volcanic History of the lmbrium Basin:A Close-up View from the Lunar Rover Yutu[J]. PNAS, 2015, 112(17): 5342-5347. DOI:10.1073/pnas.1503082112 |
| [5] | Ling Zongcheng, Jolliff B L, Wang Alian, et al. Correlated Compositional and Mineralogical Investigatios at the Chang'e-3 Landing Site[J]. Nature Communications, 2015. DOI:10.1038/ncomms9880/www.nature.com/naturecommunications |
| [6] |
刘菁华, 王祝文, 陈圣波, 等. CE-1伽马能谱测量中4种月球岩石的能谱特征分析[J].
地学前缘, 2014, 21(6): 121-128.
Liu Jinghua, Wang Zhuwen, Chen Shengbo, et al. Spectrum Characteristics Analysis of Four Lunar Rocks Based on the CE-1 Gamma-Ray Spectrometry Data[J]. Earth Sicence Frontiers, 2014, 21(6): 121-128. |
| [7] |
刘菁华, 王祝文.
放射性勘探方法[M]. 北京: 地质出版社, 2014.
Liu Jinghua, Wang Zhuwen. Radioactive Exploration Method[M]. Beijing: Geological Publishing House, 2014. |
| [8] | Zhao Jiannan, Huang Jun, Qiao Le, et al. Geologic Characteristics of the Chang'e-3 Exploration Region[J]. Science China Physics, Mechanics & Astronomy March, 2014, 57(3): 569-576. DOI:10.1007/s11433-014-5399-z |
| [9] |
张明皓, 陈超, 兰瑞平, 等. 月球表面多种金属元素的分布特征初探[J].
地学前缘, 2007, 14(5): 277-284.
Zhang Minghao, Chen Chao, Lan Ruiping, et al. APreliminary Analysis:The Distribution of Metallic Elements on the Lunar Surface[J]. Earth Science Frontiers, 2007, 14(5): 277-284. |
| [10] |
李泳泉, 刘建忠, 欧阳自远, 等. 月球表面岩石类型的分布特征:基于Lunar Prospector(LP)伽马射线谱仪探测数据的反演[J].
岩石学报, 2007, 23(5): 1169-1174.
Li Yongquan, Liu Jianzhong, Ouyang Ziyuan, et al. Petrologic Distrbutions on the Moon:Based on the Lunar Prospector(LP) Gamma-Ray Spectrometer Data Transformation[J]. Acta Petrologica Sinica, 2007, 23(5): 1169-1174. |
| [11] |
胥涛, 李春来. 月球表面元素含量的定量分析方法[J].
空间科学学报, 2001, 21(1): 332-340.
Xu Tao, Li Chunlai. Abundances and Distribution of Elements of the Lunar Surface[J]. Chinese Journal of Space Science, 2001, 21(1): 332-340. |
| [12] |
徐义刚. 月海玄武岩与月球演化[J].
地球化学, 2010, 39(1): 50-62.
Xu Yigang. Mare Basalts and Lunar Evolution[J]. Geochimica, 2010, 39(1): 50-62. |