目录总框架光谱反演深度积分模块一月球表面3He丰度 C₀太阳风强度计算月壤成熟度计算TiO2 含量计算模块二月壤剖面3He丰度变化函数 f(z)模块三月壤厚度 d总结总框架光谱反演深度积分核心逻辑利用遥感能够观测的信息逐步逼近无法直接观测的月壤³He储量。对任意一个像元³He 储量计算为M³He 的单位像素储量C₀月壤最表层的³He质量丰度ppb 级f(z)³He丰度随深度的剖面分布函数d月壤厚度A像元面积全球累加得总量C₀ 本身也不是直接测的通常依据太阳风注入条件、月壤成熟度以及钛铁矿TiO₂丰度之间建立的经验关系估算。每个影像因素的发展变化大致如下表格归纳参数/模块早期2007前中期2010s最新2020s–至今太阳风强度 (F FF)宏观纬度函数F ∝ cos ⁡ ϕ F \propto \cos \phiF∝cosϕ结合太阳活动周期的修正模型微观粒子输运模拟 古太阳风记录反演月壤成熟度 (OMAT)简单光谱比值I s / FeO I_s/\text{FeO}Is​/FeOOMAT 算法 神经网络多光谱/热红外联合反演 原位光谱校正钛含量 (TiO 2 \text{TiO}_2TiO2​)Clementine UVVIS 比值法易饱和去成熟度化反演 微波辅助微波穿透探测 嫦娥样品实验室标定剖面丰度 (f ( z ) f(z)f(z))单一指数衰减阿波罗样品拟合多矿物加权物理模型纳米气泡结构约束的分段函数月壤厚度 (d dd)DEM 经验公式精度粗嫦娥一号微波辐射计全球反演嫦娥三/四号 GPR 五/六号 MIMO 雷达原位成像模块一月球表面3He丰度 C₀太阳风强度计算太阳风是³He的来源。由于无法直接测量过去几十亿年的太阳风通量科学家通常通过“代理指标”来推算。早期阶段1990s - 2000s基于纬度与经度的经验模型方法 假设太阳风通量主要受月球纬度影响赤道最强两极较弱并结合经度效应迎风面与背风面。模型 使用简单的球谐函数或三角函数拟合。例如认为通量 F与纬度 ϕ 的余弦成正比 (F∝cosϕ)。局限 忽略了太阳风随时间的长期变化以及局部磁场对太阳风的屏蔽作用精度较粗。中期阶段2010s结合太阳活动周期的修正方法 引入太阳黑子数或太阳风卫星如 ACE, WIND的实测数据建立长期平均模型。进步 区分了“快太阳风”和“慢太阳风”对月球不同区域的轰击差异开始考虑地磁尾对月球的周期性屏蔽效应每月几天。最新阶段2020s - 至今基于古太阳风记录与微观模拟方法 利用嫦娥五号/六号样品中的太阳风注入深度剖面反推古代太阳风通量。同时利用粒子输运模型如 SPIS, Geant4模拟太阳风粒子与月表复杂地形/磁异常的相互作用。进步 从“宏观平均”走向“微观局地”。能够计算出特定撞击坑或磁异常区如 Reiner Gamma实际接收到的有效太阳风剂量而非仅仅是理论值。月壤成熟度计算成熟度反映了月壤暴露在表面的时间长短。暴露越久接收的太阳风越多但同时也伴随着微陨石轰击导致的气体逃逸饱和效应。早期阶段1970s - 2000s基于 Is/FeO的经验参数方法 阿波罗样品分析发现纳米铁含量 (I S ​ I_S​IS​​) 与氧化亚铁 (FeO) 的比值能很好地表征成熟度。遥感应用 早期尝试用 Clementine 紫外/可见光比值来近似但受矿物成分干扰大。关键转折2007 - 2010s光学成熟度参数 (OMAT) 的确立方法 Lucey 等人提出了著名的 OMAT 算法。利用 750nm 反射率和 950nm/750nm 比值在“铁-成熟度”散点图中通过几何距离计算 OMAT。地位 这成为了后续十几年³He估算的标准输入量解决了成分与成熟度的解耦问题。最新阶段2020s多光谱与热惯量联合反演方法 结合 Kaguya、Chandrayaan-1 以及嫦娥系列的高光谱数据引入热红外数据热惯量也能反映表层颗粒大小和孔隙度与成熟度相关。进步 修正了高钛区如静海和低钛区高地使用同一套 OMAT 公式产生的系统误差。特别是针对月背高地建立了新的成熟度标定曲线。TiO2 含量计算钛铁矿是³He的主要“储气罐”。准确获取全球钛含量分布是估算储量的基础。早期阶段1990sClementine UVVIS 比值法方法 利用 415nm/750nm 和 750nm/950nm 的比值关系。局限 在高钛区域 TiO26%容易出现饱和导致高估或低估且容易受到成熟度空间风化的干扰——成熟的土壤会变红干扰钛的光谱特征。改进阶段2000s - 2010sLucey 算法与神经网络方法 Lucey (2000) 改进了算法将成熟度作为变量引入钛的计算中实现了“去成熟度化”的钛含量反演。随后学者们引入人工神经网络ANN利用全波段光谱训练非线性模型。进步 大幅提高了中高钛区域的反演精度生成了更可靠的全球钛图。最新阶段2020s - 至今微波辐射与原位校验方法①微波遥感 利用嫦娥一号/二号的微波辐射计数据。微波穿透深度大不受表面空间风化成熟度影响直接反映月壤介电常数与钛含量强相关。②原位/样品校准 利用玉兔号/玉兔二号巡视器的 Visible-NIR Spectrometer 实测数据以及嫦娥五号/六号返回样品的实验室精确化验数据作为“真值”去校正轨道遥感数据。进步 解决了“表面风化层”与“下层新鲜土壤”钛含量不一致的问题特别是修正了月背南极-艾特肯盆地SPA这种特殊地质单元的钛含量估算偏差。模块二月壤剖面3He丰度变化函数 f(z)f(z) 描述的是³He含量随深度衰减的规律其核心是太阳风粒子注入深度与月壤颗粒捕获能力的耦合关系。早期阶段1990s–2007基于阿波罗样品的指数衰减模型方法 基于阿波罗样品在不同深度的实测数据拟合出经验公式。局限 仅适用于阿波罗着陆区月海未考虑月壤矿物组成、颗粒大小、成熟度对注入深度的影响且未区分正面/背面差异。中期阶段2010s引入物理机制的修正模型方法 结合太阳风粒子输运模拟如 SRIM、Geant4考虑不同矿物如钛铁矿、斜长石对氦离子的阻止本领差异建立分矿物的 f(z) 函数。进步 从“单一指数”走向“多组分加权”能更准确反映月海与高地、成熟与未成熟月壤的差异。最新阶段2020s–至今嫦娥五号/六号样品驱动的“纳米级”剖面重构方法 通过阶段升温提取实验发现³He主要赋存在钛铁矿表层玻璃层中的纳米气泡内其释放温度与深度相关。结合透射电镜TEM观测建立“纳米气泡深度分布模型”。突破 从“宏观平均”走向“微观结构约束” f(z) 不再是一个平滑曲线而是与矿物颗粒表面结构、玻璃化程度强相关的分段函数。模块三月壤厚度 d月壤厚度决定了积分上限直接影响总储量估算。其测量手段从“间接反演”发展到“直接成像”。早期阶段2007前基于DEM的经验公式法方法 利用月球数字高程模型DEM通过撞击坑密度、坡度等地貌参数经验性估算月壤厚度。局限 精度低仅能提供粗略的全球趋势无法反映局部变化。中期阶段2010嫦娥一号微波辐射计物理反演方法 利用多通道微波亮温数据结合三层介质辐射传输模型反演全球月壤厚度。首次给出全球³He总量约 6.6×108kg。进步 首次实现全球覆盖物理模型比经验公式更可靠。同期独立研究2010王振占等团队的微波反演成果 得出月壤厚度多数为 4–6 m5 m 以上区域占全月 43%估算³He约 10.3×108kg。意义 验证了微波反演的可靠性并提供了更精细的区域统计。原位探测阶段2013/2018嫦娥三号/四号测月雷达GPR方法 首次在月面原位获取月壤全厚度剖面分辨率达米级。突破 从“遥感反演”到“实地验证”确认了微波反演结果的准确性并揭示了局部厚度变化。高精度成像阶段2020/2024嫦娥五号/六号MIMO雷达方法 采用多输入多输出MIMO体制获取分辨率优于 3 cm 的月壤精细结构图。意义 能识别月壤内部的层理、石块分布、玻璃质含量等为 f(z) 的局部修正提供直接依据。总结本文系统梳理了利用遥感技术估算月球月壤中³He储量的技术发展脉络。其核心框架可概括为“光谱反演深度积分”即通过遥感数据反演获取月壤表层³He丰度C₀、剖面分布函数f(z)和月壤厚度d三个关键参数再进行积分计算。从早期基于阿波罗样品的简单经验模型到中期引入物理机制和全球遥感数据如嫦娥一号微波数据、Lucey OMAT算法再到最新阶段依托嫦娥五号/六号返回样品、原位雷达GPR/MIMO和微观粒子模拟如Geant4的“纳米级”精细刻画整个技术体系正朝着“微观化”、“原位化”和“多源融合”的方向不断演进。每一次技术进步都使得我们对月球³He这一潜在战略资源的认知和估算精度迈上一个新台阶。未来随着更多月球探测任务的实施和数据的积累这一估算体系将更加精确和可靠为人类未来的月球资源开发利用提供坚实的数据支撑。