探月雷达方法技术研究及嫦娥三号雷达数据解释

探月雷达方法技术研究及嫦娥三号雷达数据解释我国“嫦娥三号”(Chang’e 3,CE-3)和“嫦娥四号”(Chang’e 4,CE-4)搭载的探月雷达(Lunar Penetrating Radar,LPR)是一种装载在月球巡视器上对月表进行实地登陆观测的探测方式。其理论基础是利用天线向探测目标发射高频电磁波(1MHz~10GHz)进行探测,可以高分辨率地探测月球次表层地质结构与月壤厚度。根据探月雷达测量的数据,可以实现以下科学和工程研究目标:获得沿巡视器路径的月壤的厚度分布,并可以通过反演描述月壤和月球浅层岩石的电磁参数,为月壤层矿产资源的估算提供依据;更好地了解月球地形特征与地质构造的关系,为月球地形和地质研究提供科学数据;为详细研究月球的形成和演化提供科学依据;为未来选取月球基地的位置提供必要的科学数据。

本文针对月壤层的模型建立,电磁响应模拟,探月雷达数据预处理,月壤参数估计,近月表地质解释等方面进行研究。开展了月球浅表层介质的建模及其电磁响应计算与分析,为探月雷达探测机制和数据解释提供方法技术支撑;研发了嫦娥三号探月雷达数据的解算、校正和初步处理的完整流程,对雷达数据中的反射信号进行分析;设计实现了月壤、月岩物理参数的高精度估计、反演和分析,量化潜在可开发的月球资源,帮助分析资源分布规律;解释了近月表地质界面,有利于确定界面的地质意义和介质组成,为加深对浅月表层的形成演化认识提供基础,也为进一步的月球以及深空探测提供支持。(1)建立了月壤综合模型。

通过将分形技术应用于地形建模,准确地描述了小区域的月球地形,引入磨粒技术来模拟不规则的异常体,应用随机等效介质描述随机介质。同时也考虑了介电常数随深度的变化。月壤综合模型的建立加深了对月壤形成和演化的认识。

根据时间域有限差分法对建立的月壤综合模型进行电磁响应模拟。与简单模

型的仿真结果相比,月壤综合模型的仿真结果包含更多细节,更加真实。(2)处理了嫦娥三号探月雷达数据。

介绍嫦娥三号着落点的区域信息,研究我国探月工程所采集的探月雷达数据体,结合在月球表面采集数据的复杂环境以及雷达仪器本身的影响,分别针对第一通道(CH-1,60MHz)和第二通道(CH-2,500MHz)进行处理流程设计,研究包括带通滤波以及增益等常规方法在内的处理方法对探月雷达数据处理的效果,最后得到了清晰的探月雷达图像并分析了数据存在的问题。(3)研究了嫦娥三号着陆区月壤层的结构和演变过程。根据正演模拟的结果,探讨了月壤与基岩接触面的形态。

随后,利用F-K(频率-波数)滤波突出了月壤和基岩的接触面。根据探月雷达第二通道数据能量的分布对着陆区月壤层进行分层。通过综合考虑月球的演变历史,区域地质情况,特别是探月雷达数据结果,推断了嫦娥三号着陆区月壤层的演变过程。

(4)估计了月壤的电性参数和铁钛含量。从理论上证明了跟据两组收发距不同的雷达数据,利用两者对同一个异常体反射时间的不同,从而推断介质介电常数的方法具有可行性。针对一个简单模型和一个月壤复杂模型,利用该参数估计方法,反推出各异常体的位置和模型的介质参数,得到结果与实际模型一致。

利用探月雷达第二通道的两套收发距不同的数据,得到了嫦娥三号着陆点月壤的介电常数。Apollo计划的成果给出,介电常数与月壤的密度和铁钛含量存在内在联系,结合估计得到的月壤介电常数,可以估算得到月壤的铁钛含量。(5)定量分析了月壤中孤石的位置和数量。

提出了一种基于局部相似性约束的自适应石块定位方法,该方法利用到了探

月雷达第二通道两套数据之间的相似性。该方法的实现首先采用基于f-x EMD

的倾角滤波器来提取两个数据中含有信号的波数分量。然后,计算滤波后的两套数据之间的局部相似谱。

通过利用软阈值函数,认为相似谱中的局部最大值可作为每个孤石的位置。最后,根据提取结果,一方面获得了月壤的深度,另一方面也揭示了月壤中孤石随路径和深度变化的分布信息。进一步为嫦娥五号月壤采集提供了先验方法。

(6)分析了嫦娥三号着陆点浅月表地质结构。嫦娥三号第一通道数据探测到了虹湾地区的浅表层地质结构。由于数据信噪比较小,为了解决噪声的限制,本文提出了一种基于shearlet变换的去噪策略。

在该策略中,shearlet变换不仅是一种数据分析的工具,还可以提高CH-1数据的信噪比并突出同相轴。处理后质量提高的CH-1数据有助于地质分层和解释。在分别获得浅部和深部的处理数据之后,跟踪处理结果中的同相轴,得12个可能的分界面。

在评估每个同相轴后,根据可靠性对每个同相轴进行判断和分类。利用探月雷达数据分析浅月表地质结构提供了有价值的地质信息,并为未来的月球探测任务提供了参考。