文章快速检索  
  高级检索
月球风化层钻取采样过程密实度分类研究
郑燕红1, 邓湘金1, 庞勇2, 金晟毅1, 姚猛1, 赵志晖1     
1. 北京空间飞行器总体设计部, 北京 100094;
2. 北京卫星制造厂, 北京 100094
摘要: 钻取采样是月球风化层土壤样品获取的重要方式,密实度是重要的风化层月壤原位特性,对钻进过程中的策略制定有重要影响。本文结合钻取采样过程特点,提出了通过采样机构的力、速度、电流、温度等传感器获取的瞬时信息感知月壤密实度的方法,利用深度学习方法构建一类适应于可变长度序列数据的门控型循环神经网络,实现钻进过程月壤密实度在线分类。研究表明,该分类方法在风化层钻进过程中月壤密实度感知滞后时间约为33 s,对未知序列数据识别正确率大于89.36%,具有较高的分类精度和泛化能力。
关键词: 月球    钻取采样    月壤密实度    循环神经网络    在线分类    
Research on classification of relative density in lunar regolith drilling
ZHENG Yanhong1, DENG Xiangjin1, PANG Yong2, JIN Shengyi1, YAO Meng1, ZHAO Zhihui1     
1. Beijing Institute of Space System and Engineering, Beijing 100094, China;
2. Beijing Spacecrafts, Beijing 100094, China
Abstract: Drilling is an important approach to acquire samples from lunar regolith. The relative density is an in situ characteristic of lunar regolith and will influence the strategy of drilling process. According to the working principle of drilling, the perception scheme of lunar soil relative density is proposed with the instantaneous information derived from sensors of force, velocity, current, and temperature. Based on the deep learning method, a class of gated recurrent unit neural network is constructed for variable length data sequence. The network can achieve the online classification of relative density of lunar soil. The delay time for relative density perception is about 33 s and the correct ratio is higher than 89.36% to unknown data sequence, indicating that the method has high classification precision and generality.
Keywords: moon    drilling    lunar soil relative density    recurrent neural network    online classification    

中国月球探测按照“绕、落、回”三步走战略规划,正在开展第三步“回”的研制工作,将发射嫦娥五号探测器实施无人采样返回任务[1-2],通过钻取采样获取月球风化层月壤样品是任务的重要环节之一。早在20世纪70年代,钻取采样在Apollo任务、Luna任务中就作为一种重要的采样方式获取了大量月球样品[3-5]。随着月球南极成为月球探测新的热点,美国、俄罗斯新发布的深空探测器计划[6-7],均将月球采样作为重返月球的重要标志之一,钻取采样必将发挥重要作用。

钻取采样的对象是月球风化层的月壤,不同于地球土壤,其不受空气、流动水、微生物等因素影响,大多数位置的岩土特性处于较为稳定状态,然而由于小天体撞击等影响,月壤密实度在不同位置、不同深度则呈现较大差异。月壤密实度尤其是钻进过程中的在线密实状态,是钻取采样策略制定的重要依据。目前月壤密实度信息仍主要以离线方式获得,测定途径主要利用月壤样品最大最小密度、宇航员脚印、月球车轨迹、月面漂砾轨迹、月壤耐穿透性等。这类测量结果可为钻取采样产品研制提供边界,但对月球地质构造的不可预知性、不同采样位置的随机性反映不足,Luna 16、Apollo 15任务实施过程中,均出现过密实月壤卡滞钻杆的情况[8]。因此更快捷地获取钻取采样过程中的月壤密实度状态对钻取工作的顺利进行具有重要意义。文献[9]综合考虑月壤、岩石情况,分析了钻进扭矩、钻压力的时域与频域特征,利用支持向量机(Support Vector Mechine, SVM)方法实现了高准确率的可钻性等级辨识,但可钻性对差异较小的风化层月壤刻画略有不足。文献[10]通过加速度计或振动仪测量钻取采样机构在多个方向的频率响应,通过神经网络自主分析钻取采样机构健康状态,在地面试验中取得了较好的效果,但其关注点主要集中在机构上,在对象状态分析方面偏少。

本文针对月球风化层无人自动钻取采样过程,通过合理配置轻小型月壤状态感知传感器,获取月壤状态瞬时信息,构造一类适用于月壤密实度分类的门控型循环神经网络(GRU),对长度可变的信息序列进行自主分类识别,通过地面试验数据建立分类模型,并验证其泛化能力,为月壤原位特性分析与在轨钻进策略制定提供技术支持。

1 钻取采样过程

月球钻取采样可分为硬管取芯与软管取芯2种方式,软管取芯方式可同步完成月壤钻进与样品收纳,并可保持风化层月壤的层序信息,相对于硬管取芯其具有质量轻、易回收等优点,有利于无人自动月球采样返回任务的实施[11]。钻取采样获取样品是通过机构与月壤的相互作用,使局部月壤与相邻区域月壤发生分离的过程。采样机构对原位月壤施加纵向下压力、横向切削力[12],对可能存在的硬物质提供冲击破碎力,对软管取芯还提供一定的软管保持力,通常包含加载机构、钻进机构、提芯机构,如图 1所示。加载机构利用电机卷绕绳系带动钻进机构以一定的速度向下运动提供加载力;钻进机构通过电机驱动钻杆以一定转速或冲击频率对风化层月壤进行回转或冲击破碎,提供切削力和冲击力;提芯机构通过拉绳与内嵌软管连接,提供软管保持力,确保软管可用容积不断增加,降低月壤进入阻力。

图 1 钻取采样机构示意图 Fig. 1 Sketch of drilling mechanism

月球风化层钻取采样对象是具有一定深度剖面的月壤,由不同的层构成,各层月壤的粒径分布、密实度等存在随机性,层与层沿纵向变化的不连续性,采样作用对原位月壤引入的扰动性,导致螺旋钻杆在各层的负载变得不确定[13-14],对钻进机构、加载机构工作策略适应性带来了挑战。对采样机构而言,不同月壤状态将导致机构作用力、力矩、速度、温度发生变化,因此可设置相应的传感器对这类状态信号进行获取,从而感知风化层月壤状态信息。从图 1可见,钻取采样纵向运动的钻进速度v(可通过位移s差分)、加载力ft可通过加载机构进行控制,径向运动的回转转速nr、回转扭矩(与电流i成正比)可通过钻进机构进行控制,同时提芯拉绳的提拉力fk可对芯管收纳月壤的阻力进行感知,而加载机构、钻进机构的电机随着工作时间、负载的变化,其工作温度tttr也随之变化。因此对加载机构可配置位移、力、温度传感器分别对钻进速度、加载力、加载电机温度变化情况进行测量;对钻进机构配置速度、电流、温度传感器对回转转速、回转扭矩、回转电机温度变化情况进行测量;对提芯机构配置力传感器对提供的提芯力进行测量。

图 2所示,在钻取采样过程中,通过上述传感器获取的反映风化层月壤状态的瞬时信息x,可表示为

图 2 钻取采样月壤状态感知示意图 Fig. 2 Sketch of state perception in drilling lunar soil
$ \boldsymbol{x}=\left[\begin{array}{lllllll} v & f_{\mathrm{t}} & n_{\mathrm{r}} & i & f_{\mathrm{k}} & t_{\mathrm{t}} & t_{\mathrm{r}} \end{array}\right]^{\mathrm{T}} $ (1)
 

以一定时间段内的瞬时信息构成序列输入,通过有效的状态感知方法,生成表征局部月壤特性的输出信息c,为钻取采样策略生成与控制提供依据信息,形成器上自主或天地协同回路。

2 月壤密实度

根据第1节的描述,钻取采样过程是在加载机构、钻进机构、提芯机构的作用下,通过剪切破碎使局部区域月壤发生原位分离的过程,局部区域月壤的抗剪性是影响钻具与月壤相互作用的重要因素,而风化层月壤的抗剪性与月壤密实度正相关,因此钻取采样过程中风化层月壤密实度可作为状态感知的特征输出。密实度是反映颗粒土壤聚集情况的重要参数,对于颗粒型土壤,密实度dr的表达式为

$ d_{\mathrm{r}}=\frac{\rho_{\max }}{\rho} \times \frac{\rho-\rho_{\min }}{\rho_{\max }-\rho_{\min }} \times 100 \% $ (2)
 

式中:ρmaxρmin分别为最大、最小密度;ρ为当前状态密度。根据密实度定义,风化层月壤可根据密实程度分为5类[15],如表 1所示。

表 1 月壤密实度分类 Table 1 Relative density types of lunar soil
月壤密实度 状态分类
0~15% 极疏
15%~35% 疏松
35%~65% 中密
65%~85% 密实
85%~100% 极密

Apollo任务证实,月壤与地球土壤的密实度存在较大差异,对浅层月壤密实度最佳估计表明,月表15 cm以内的月壤密实度约65%,当深度超过30 cm时密实度将急剧上升至92%左右[8]。因此,月球风化层钻取采样过程中的月壤密实度主要集中在中密、密实与极密范围内。

月球探测器采用的钻取采样机构往往具有较为严格的重量约束,仅能通过轻小型传感器获取力、位移、转速、电流、温度等状态量,难以配置复杂的传感器直接获取表征月壤状态的瞬时密实度。另一方面,钻取采样控制系统为确保系统的鲁棒性,通常具有较好的扰动抑制能力[16],过于精细的月壤密实度信息将增大控制系统设计难度。因此,采用既具有一定区分度,又具有一定容忍度的月壤密实度刻画可与后端策略生成与控制形成良好的匹配。

中国探月工程实施过程中,对模拟月壤进行了广泛的研究[17-18],采样模拟月壤与月面软着陆、月面巡视探测研究的模拟月壤有所不同,其对月壤密实度的关注度更高。其从物质成份、颗粒形态、颗粒密度、水份含量等方面进行基础原料控制,按照一定的粒径级配进行混杂,通过控制如图 3所示振动台的振动量级、振动时间获取所需的模拟月壤当前状态密度ρ,从而根据式(2)制备不同密实程度的模拟月壤分层,然后将各个分层进行拼接组成具有一定深度的钻取采样模拟月壤。利用该方法,在地面试验过程中可制备中密、密实、极密3类不同密实度状态的模拟月壤,开展大量的钻取采样试验验证,试验装置如图 4所示,获取各传感器的特征数据,建立传感器状态量与模拟月壤密实度的感知模型。

图 3 模拟月壤制备示意图 Fig. 3 Scheme of lunar soil simulant fabricate
图 4 钻取采样试验装置 Fig. 4 Test facility for ground drilling

钻取采样在轨实施过程中,可利用获取的感知模型,完成风化层月壤密实度的分类预测,为在轨钻进策略制定提供支持。

3 感知分类

由于采集到的各传感器数据为基于时间的序列信息,门控型循环神经网络对该类数据处理具有广泛的应用,并取得了良好的效果[19-22]。针对风化层月壤钻取采样过程,设计如图 5所示密实度分类网络。

图 5 月壤密实度分类网络 Fig. 5 Classification network for relative density of lunar soil

月球风化层钻取采样过程中,探测器以一定的时间间隔下传传感器数据,每帧下传数据可作为一个时间步数据x(k)Rp×1(k=1, 2,…, n),其中p为传感器特征数据维数。设用于密实度分类的数据序列长度为n,可围绕标称值上下波动,每个时间步单元Tk具有相同的门控结构,重置门、更新门输出取值范围为0~1,具有门控特性,重置门输出r(k)决定前一时间步信息量h(k-1)的遗忘因子,更新门输出z(k)决定传递至后一时间步的信息量h(k)的组合因子。门控单元Tk可描述为

$ \left\{\begin{aligned} \boldsymbol{r}^{(k)} &=\sigma\left(\boldsymbol{u}_{r}^{(k)}\right)=\sigma\left(\boldsymbol{W}_{rh} \boldsymbol{h}^{(k-1)}+\boldsymbol{W}_{r x} \boldsymbol{x}^{(k)}+\boldsymbol{b}_{r}\right) \\ \boldsymbol{z}^{(k)} &=\sigma\left(\boldsymbol{u}_{z}^{(k)}\right)=\sigma\left(\boldsymbol{W}_{z h} \boldsymbol{h}^{(k-1)}+\boldsymbol{W}_{z x} \boldsymbol{x}^{(k)}+\boldsymbol{b}_{z}\right) \\ \boldsymbol{g}^{(k)} &=\tanh \left(\boldsymbol{u}_{g}^{(k)}\right)=\tanh \left(\boldsymbol{W}_{g r h} \boldsymbol{r}^{(k)} \cdot \boldsymbol{h}^{(k-1)}+\right.\\ &\left.\boldsymbol{W}_{g x} \boldsymbol{x}^{(k)}+\boldsymbol{b}_{g}\right) \\ \boldsymbol{h}^{(k)} &=\left(1-\boldsymbol{z}^{(k)}\right) \cdot \boldsymbol{h}^{(k-1)}+\boldsymbol{z}^{(k)} \cdot \boldsymbol{g}^{(k)} \end{aligned}\right. $ (3)
 

式中:g(k)为组合信息;ur(k)uz(k)ug(k)为净输入。设门控单元内部维数为q,则g(k)r(k)z(k)h(k)h(k-1)Rq×1,权值矩阵WrhWzhWgrhRq×qWrxWzxWgxRq×p,偏置矩阵brbzbgRq×1。门控单元采用了sigmoid和双曲正切激活函数:

$ \left\{\begin{array}{l} \sigma(x)=\frac{1}{1+\exp (-x)} \\ \tanh (x)=2 \sigma(2 x)-1 \end{array}\right. $ (4)
 

风化层月壤钻取采样过程中的传感器数据序列经门控单元循环处理,不断提取数据特征,为区分月壤的密实度情况,以最后一个门控单元输出h(n)为输入,通过全连接层与Softmax层对当前钻进月壤密实度进行分类。设Softmax层输出为c类,全连接层神经元数量为m,则其数学描述为

$ \left\{\begin{array}{l} \boldsymbol{Z}_{\mathrm{Fc}}=\boldsymbol{W}_{\mathrm{Fe}} \boldsymbol{h}^{(n)}+\boldsymbol{B}_{\mathrm{Fc}} \\ \boldsymbol{Y}_{\mathrm{Fc}}=\operatorname{relu}\left(\boldsymbol{Z}_{\mathrm{Fc}}\right) \\ \boldsymbol{Y}_{\mathrm{s}}=\operatorname{softmax}\left(\boldsymbol{W}_{\mathrm{s}} \boldsymbol{Y}_{\mathrm{Fc}}+\boldsymbol{B}_{\mathrm{s}}\right) \end{array}\right. $ (5)
 

式中:WFcRm×qWsRc×m为权值矩阵,BFcRm×1BsRc×1为偏置矩阵,ZFcRm×1为全连接层加权输入,YFcRm×1YsRc×1分别为全连接层、Softmax层输出。全连接层、Softmax层激活函数描述为

$ \left\{\begin{array}{ll} \operatorname{relu}(x)=\left\{\begin{array}{ll} x & x>0 \\ 0 & x \leqslant 0 \end{array}\right. \\ \operatorname{softmax}(x)= \frac{\exp (x)}{\sum\limits_{i=1}^{c} \exp \left(x_{i}\right)} \end{array}\right. $ (6)
 

为建立传感器特征数据与月壤密实度的映射关系,选取如式(7)所示交叉熵函数为目标函数L,进行监督学习。

$ L=-\sum\limits_{o=1}^{c} \boldsymbol{d}(o) \ln \boldsymbol{Y}_{\mathrm{s}}(o)+\frac{1}{2} \lambda \sum\limits_{i, j} W_{*}(i, j)^{2} $ (7)
 

式中:d为训练样本目标输出;o为输出类别序号;λ为L2正则系数;W*表示各层权值矩阵;ij表示其元素序号。根据反向传播原理[23-24],目标函数对Softmax层、全连接层权值偏导数为

$ \left\{\begin{array}{l} \frac{\partial L}{\partial \boldsymbol{W}_{\mathrm{s}}}=e \cdot \boldsymbol{Y}_{\mathrm{Fc}}+\lambda \boldsymbol{W}_{\mathrm{s}} \\ \frac{\partial L}{\partial \boldsymbol{W}_{\mathrm{Fc}}}=\boldsymbol{W}_{\mathrm{s}}^{\mathrm{T}} \boldsymbol{e} \cdot \operatorname{relu}^{\prime}\left(\boldsymbol{Z}_{\mathrm{Fc}}\right)\left(\boldsymbol{h}^{(n)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{\mathrm{Fc}} \end{array}\right. $ (8)
 

式中:输出偏差e=Ys-d;relu′(·)为relu函数导函数。同理可得目标函数对偏置矩阵的偏导数为

$ \left\{\begin{array}{l} \frac{\partial L}{\partial \boldsymbol{B}_{\mathrm{s}}}=e \\ \frac{\partial L}{\partial \boldsymbol{B}_{\mathrm{Fc}}}=\boldsymbol{W}_{\mathrm{s}}^{\mathrm{T}} \boldsymbol{e} \cdot \operatorname{relu}^{\prime}\left(\boldsymbol{Z}_{\mathrm{Fc}}\right) \end{array}\right. $ (9)
 

随着误差信号反向传播至各时间步单元,可得其权值梯度和偏置梯度分别为

$ \left\{\begin{array}{l} \frac{\partial L}{\partial \boldsymbol{W}_{g x}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{g}^{(k)}\left(\boldsymbol{x}^{(k)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{g x} \\ \frac{\partial L}{\partial \boldsymbol{W}_{g r h}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{g}^{(k)}\left(\boldsymbol{h}^{(k-1)} \cdot \boldsymbol{r}^{(k)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{g r h} \\ \frac{\partial L}{\partial \boldsymbol{W}_{z x}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{z}^{(k)}\left(\boldsymbol{x}^{(k)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{z x} \\ \frac{\partial L}{\partial \boldsymbol{W}_{z h}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{z}^{(k)}\left(\boldsymbol{h}^{(k-1)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{z h} \\ \frac{\partial L}{\partial \boldsymbol{W}_{r x}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{r}^{(k)}\left(\boldsymbol{x}^{(k)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{r x} \\ \frac{\partial L}{\partial \boldsymbol{W}_{rh}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{r}^{(k)}\left(\boldsymbol{h}^{(k-1)}\right)^{\mathrm{T}}+\lambda \boldsymbol{W}_{rh} \end{array}\right. $ (10)
 
$ \left\{\begin{array}{l} \frac{\partial L}{\partial \boldsymbol{b}_{r}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{r}^{(k)} \\ \frac{\partial L}{\partial \boldsymbol{b}_{z}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{z}^{(k)} \\ \frac{\partial L}{\partial \boldsymbol{b}_{g}}=\sum\limits_{k=1}^{n} \boldsymbol{\delta}_{g}^{(k)} \end{array}\right. $ (11)
 

式中:组合信息灵敏度δg(k)=∂L/∂ug(k),更新门灵敏度δz(k)=∂L/∂uz(k),重置门灵敏度δr(k)=∂L/ur(k),其递推关系为

$ \left\{\begin{array}{ll} \boldsymbol{\delta}_{h}^{(k)}=\left\{\begin{array}{ll} \boldsymbol{W}_{\mathrm{Fc}}^{\mathrm{T}} \boldsymbol{W}_{\mathrm{s}}^{\mathrm{T}} \boldsymbol{e} \cdot \operatorname{relu}^{\prime}\left(\boldsymbol{Z}_{\mathrm{Fc}}\right) & k=n \\ \boldsymbol{\delta}_{h}^{(k+1)}\left(1-\boldsymbol{z}^{(k+1)}\right)+\boldsymbol{W}_{zh}^{\mathrm{T}} \boldsymbol{\delta}_{z}^{(k+1)}+ & \\ \boldsymbol{W}_{{grh }}^{(k+1)} \boldsymbol{\delta}_{g}^{(k+1)} \boldsymbol{r}^{(k+1)}+\boldsymbol{W}_{rh}^{(k+1)} \boldsymbol{\delta}_{r}^{(k+1)} & \\ & k < n \end{array}\right.\\ \boldsymbol{\delta}_{g}^{(k)}=\boldsymbol{\delta}_{h}^{(k)} \cdot \boldsymbol{z}^{(k)} \cdot\left(\bf{1}-\boldsymbol{g}^{(k)} \cdot \boldsymbol{g}^{(k)}\right)\\ \boldsymbol{\delta}_{r}^{(k)}=\boldsymbol{W}_{g r h}^{\mathrm{T}} \boldsymbol{\delta}_{g}^{(k)} \cdot \boldsymbol{h}^{(k-1)} \cdot \boldsymbol{r}^{(k)} \cdot\left(\bf{1}-\boldsymbol{r}^{(k)}\right)\\ \boldsymbol{\delta}_{z}^{(k)}=\boldsymbol{\delta}_{h}^{(k)} \cdot\left(\boldsymbol{g}^{(k)}-\boldsymbol{h}^{(k-1)}\right) \cdot \boldsymbol{z}^{(k)} \cdot\left(\bf{1}-\boldsymbol{z}^{(k)}\right) \end{array}\right. $ (12)
 

式中: δh(k)为目标函数对各隐层输出h(k)的偏导。

为实现各层、各门控单元的权值、偏置迭代更新,对获取的传感器数据分小批量训练,每一批进行一次权值及偏置更新,即

$ \boldsymbol{\theta}_{*}^{b+1}=\boldsymbol{\theta}_{*}^{b}-\eta \frac{1}{N} \sum\limits_{l=1}^{N}\left(\frac{\partial L}{\partial \boldsymbol{\theta}_{*}}\right)^{(l)} $ (13)
 

式中:θ*表示各权值或偏置;b为样本批号;N为每批样本量;l为样本序号;η为学习率。

训练过程中,对目标函数及样本分类正确率进行交叉验证,建立月球风化层密实度分类模型;预测过程中,利用获得的模型权值及偏置进行正向计算,实现月球风化层钻取采样过程密实度快速感知分类。

4 验证与分析

为验证密实度感知分类算法的正确性,以不同密实状态模拟月壤钻取采样试验数据基础,将获取的传感器数据分为2组(对应的密实度状态已知),第1组用于网络的训练与验证,获取分类算法映射参数(权值及偏置),第2组数据用于分类预测,确认密实度分类算法的泛化能力。

考虑到月壤粒径级配存在一定波动,钻取采样试验中采用的颗粒直径在0~41 mm范围内,形状为棱角状或次棱角状,按照如图 6所示的4种粒径级配方案进行制备,其中小于1 mm颗粒占主导地位,最小占比不少于60%,大于10 mm颗粒最大占比不超过8.7%,大于20 mm颗粒最大占比不超过6.3%,制备过程中各种粒径颗粒摆放位置未进行特殊设计,采用了随机排布方式。结合中密、密实及极密状态相对密实度要求,利用图 3所示振动台获得了如表 2所示的7种模拟月壤,通过如图 4所示钻取采样过程试验,获取分类算法验证所需的2组数据。

图 6 模拟月壤粒径分布 Fig. 6 Particle distribution of simulant lunar soil
表 2 钻取采样过程使用的模拟月壤状态 Table 2 Simulant lunar soil state in ground drilling process
模拟月壤标号 密实度 状态分类 测试编号
Ⅰ-B 60%±5% 中密 1
Ⅰ-C 60%±5% 中密 1
Ⅰ-D 60%±5% 中密 1
Ⅱ-A 87%±2% 密实 1,2
Ⅱ-B 87%±2% 密实 1,2
Ⅲ-B 97%±2% 极密 1,2
Ⅲ-C 97%±2% 极密 1

图 7所示选取模拟月壤Ⅰ-C-1、Ⅰ-D-1、Ⅲ-C-1、Ⅲ-B-1、Ⅱ-A-1、Ⅱ-B-1 (图中从上至下按此顺序)试验中获取的传感器数据为第1组。图中Ⅱ-B-1、Ⅲ-B-1转速、温度、电流存在短暂下落是由于临时停机所致。从图 7(c)图 7(d)可见,钻杆未与模拟月壤接触时,回转电机停止、加载电机运动,此时为无月壤状态,当钻杆与模拟月壤接触后,两电机均运动,此时为月壤负载状态,因此将网络分类c设置为无月壤、中密、密实、极密4类。训练前,将第1组数据分别按序列长度n∈[25, 40]、n∈[22, 44]进行随机分割,结合已知密实度状态,分别生成训练数据与验证数据。从图 7可见,钻取采样过程中根据钻进状态对回转速度、钻进速度进行了适应性调整(不作调整可能发生卡钻或样品量下降),加载力、提芯力等随之变化,同时不同粒径颗粒在空间分布的随机性也对部分传感器数据带来了影响,因此,依靠单个传感器数据、多传感器单帧数据都难以对密实度状态进行正确地识别。

图 7 第1组传感器数据 Fig. 7 Data of the first group sensors

根据前述网络结构,设置门控单元内神经元数量为100,全连接层神经元为40,学习率η=0.01,正则系数λ=0.000 5。根据式(3)、式(5)、式(8)、式(10)、式(11)进行正反向传播,并按式(13)进行网络更新,每5次训练迭代状态如图 8所示。

图 8 训练过程变化曲线 Fig. 8 Curves of training process

图 8(b)整体来看,目标函数值逐步下降,经285次迭代后趋于稳态,对验证数据的分类正确率为91.27%。整个网络36 604个训练参数数量较多,利用图 9进行图示,各层权值、偏置范围见图中标注。

图 9 网络参数训练结果 Fig. 9 Training results of network parameters

为确认密实度分类算法的泛化能力,选取未参与训练的模拟月壤Ⅰ-B-1、Ⅱ-A-2、Ⅱ-B-2、Ⅲ-B-2(图 10中中从上至下按此顺序)试验中的传感器数据为第2组,其输入序列长度先后按照n∈[22, 27]、[28, 33]、[34, 38]、[39, 44]进行随机分割,结合已知密实度状态,生成预测数据。

图 10 第2组传感器数据曲线 Fig. 10 Data of the second group sensors

对各预测数据序列进行归一化处理后,利用如图 9所示映射参数,进行密实度状态分类,通过与真实模拟月壤状态进行比对,其预测正确率不低于89.36%,详细结果如表 3所示。

表 3 第2组数据预测结果 Table 3 Predicted results of the second group data
月壤类型 输入序列数 正确预测数 正确率/%
无月壤 170 152 89.41
中密 142 129 90.85
密实 235 210 89.36
极密 134 122 91.04

表 3可见,密实度分类算法对不同状态下的密实度状态,具有较好的泛化能力。结合图 6表 2及钻进初始的空行程运动可知,训练过程中验证数据包含了4种粒径级配,预测过程中测试数据包含了2种粒径级配,网络对训练数据、预测数据的总密实度分类正确率优于90%,验证数据、测试数据对应的7种模拟月壤,通过对典型粒径、密实度交叉组合,有效降低了粒径级配对网络提取密实度状态特征的影响。随着钻取采样过程中传感器特征数据随遥测不断下传,试验过程中以一定时间间隔(代替序列长度)组织各时间步输入,利用密实度分类算法实现小滞后密实度状态感知,最小时间间隔约为22 s(试验中遥测下传时间间隔为1 s),平均时间间隔约为33 s,可为工作策略的维持或调整提供依据。

5 结论

1) 提出了钻取采样机构传感器特征数据感知月球风化层月壤状态方法,实现钻取采样数据信息对就位月壤密实度的感知,并支持钻取采样过程自主策略调整。

2) 构建了一类适用于月球风化层月壤密实度感知分类的深度神经网络,推导了其正、反向传播关系,可适应于变长度的传感器特征数据序列。

3) 利用不同状态的模拟月壤开展钻取采样试验,获取了传感器特征数据与模拟月壤密实度间的感知分类模型,对钻取采样过程中未知传感器数据序列,密实度感知分类正确率大于89.36%,平均感知滞后时间约为33 s(遥测下传时间间隔为1 s),具有较好的泛化能力和原位特性在线识别能力。

4) 月球风化层密实度是重要的原位特性,该方法对密实度所处范围具有较好的区分度,后续可结合钻取采样机构工作特性、模拟月壤颗粒分布随机性,进行更为精细的密实度刻画与反演研究。

参考文献
[1] 彭兢, 黄昊, 向开恒, 等. 月球无人采样返回任务概念设想[J]. 航天器工程, 2010, 19(5): 99-104.
PENG J, HUANG H, XIANG K H, et al. Conception design of a lunar robotic sampling and return mission[J]. Spacecraft Engineering, 2010, 19(5): 99-104. (in Chinese)
Cited By in Cnki (22) | Click to display the text
[2] LI F, MAO Y, YAN J G, et al. A simulation of the four-way lunar lander-orbiter tracking mode for the Chang'E-5 mission[J]. Advances in Space Research, 2016, 2016(57): 2376-2384.
Click to display the text
[3] KELLY T J. A review of the Apollo lunar module program and its lessons for future space missions[C]//AIAA Space Programs and Technologies Conference.Reston, VA: AIAA, 1990: 1-7.
[4] ROBINSON M S, PLESCIA J B, JOLLIFF B L, et al. Soviet lunar sample return missions:Landing site identification and geologic context[J]. Planetary and Space Science, 2012, 69(2012): 76-88.
Click to display the text
[5] DJACHKOVA M V, LITVAK M L, MITROFANOV I G, et al. Selection of Luna-25 landing sites in the south polar region of the moon[J]. Solar System Research, 2017, 51(3): 185-195.
Click to display the text
[6] 郑燕红, 邓湘金, 赵志晖, 等. 地外天体采样任务特点及关键技术发展建议[J]. 探矿工程(岩土钻掘工程), 2014, 41(9): 71-74.
ZHENG Y H, DENG X J, ZHAO Z H, et al. The character and key technology suggestion of extraterrestrial sampling mission[J]. Exploration Engineering (Rock&Soil Drilling and Tunneling), 2014, 41(9): 71-74. (in Chinese)
Cited By in Cnki (12) | Click to display the text
[7] NASA. National space exploration campaign report.[R].Washington, D.C.: NASA, 2018.
[8] HEIKEN G H, VANIMAN D T, FRENCH B M. Lunar sourcebook a user's guide to the moon[M]. New York: Cambrige University Press, 1991.
[9] 唐钧跃.基于可钻性在线辨识的月面钻取采样钻进规程研究[D].哈尔滨: 哈尔滨工业大学, 2014.
TANG J Y. Research on drilling strategy of lunar drilling and coring based on real-time recognition of drillability[D].Harbin: Harbin Institute of Technology, 2014(in Chinese).
Cited By in Cnki | Click to display the text
[10] STATHAM S, HANAGUD S, RUZZENE M, et al. Design and validation of an LDV-Based structural health monitoring in drilling automation for mars exploration[C]//48th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Confe4rence. Reston, VA: AIAA, 2007: 1-13.
[11] 王丽丽, 刘志全, 吴伟仁, 等. 月球钻取采样机构的钻杆结构与运动参数分析[J]. 航空学报, 2016, 37(2): 738-748.
WANG L L, LIU Z Q, WU W R, et al. Analysis of drill stem structural and kinematic parameters of lunar drilling sampling mechanism[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(2): 738-748. (in Chinese)
Cited By in Cnki (5) | Click to display the text
[12] 邓宗全, 田野, 唐德威, 等. 用于地外星体探测的一种新结构取心钻头研究[J]. 机械工程学报, 2013, 49(19): 104-110.
DENG Z Q, TIAN Y, TANG D W, et al. Research on new structure coring bit for extraterrestrial bodies exploration[J]. Journal of Mechanical Engineering, 2013, 49(19): 104-110. (in Chinese)
Cited By in Cnki (21) | Click to display the text
[13] 赖小明.月壤剖面模拟及其钻进负载特性研究[D].长沙: 国防科技大学, 2016: 68-69.
LAI X M. The cross-section simulant of lunar regolith and research on drilling load feature[D].Changsha: National University of Defense Technology, 2016: 68-69(in Chinese).
[14] 庞勇, 冯亚杰, 孙启臣, 等. 月壤大颗粒对钻进力载影响的仿真及实验研究[J]. 北京大学学报(自然科学版), 2019, 55(3): 397-403.
PANG Y, FENG Y J, SUN Q C, et al. Simulation and experimental study on the effect of large granular rocks in lunar soil on drilling load[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2019, 55(3): 397-403. (in Chinese)
Cited By in Cnki | Click to display the text
[15] LAMBE T W, WHITMAN R V. Soil mechanics[M]. New York: John Wiley&Sons, 1969.
[16] 郑燕红, 邓湘金, 姚猛, 等. 一类月面钻进采样机构的鲁棒控制[J]. 空间控制技术与应用, 2014, 40(5): 36-41.
ZHENG Y H, DENG X J, YAO M, et al. Robust control for a class of lunar drill sampling mechanism[J]. Aero-space Control and Application, 2014, 40(5): 36-41. (in Chinese)
Cited By in Cnki | Click to display the text
[17] 贾阳, 申振荣, 党兆龙, 等. 模拟月壤研究及其在月球探测工程中的应用[J]. 航天器环境工程, 2014, 31(3): 241-247.
JIA Y, SHENG Z R, DANG Z L, et al. Lunar soil simulant and its engineering application in lunar exploration program[J]. Spacecraft Environment Engineering, 2014, 31(3): 241-247. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[18] 邹猛, 李建桥, 何玲, 等. 不同粒径分布模拟月壤承压特性试验研究[J]. 航空学报, 2012, 33(12): 2338-2345.
ZOU M, LI J Q, HE L, et al. Experiment study on the pressure-sinkage characteristic of the simulant lunar regolith with different particle size distribution[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(12): 2338-2345. (in Chinese)
Cited By in Cnki (16) | Click to display the text
[19] GRAVES A. Supervised sequence labeling with recurrent neural networks[M]. Berlin Heidelberg: Springer-Verlag, 2012.
[20] 王鑫, 吴际, 刘超, 等. 基于LSTM循环神经网络的故障时间序列预测[J]. 北京航空航天大学学报, 2018, 44(4): 772-784.
WANG X, WU J, LIU C, et al. Exploring LSTM based recurrent neural network for failure time series prediction[J]. Journal of Beijing University of Aeronautics and Astronautics, 2018, 44(4): 772-784. (in Chinese)
Cited By in Cnki (8) | Click to display the text
[21] LIPTON Z C, BERKOWITZ J. A critical review of recurrent neural networks for sequence learning[J]. Computer Science, 2015, 23-30.
Click to display the text
[22] DEY R, SALEM F M. Gate-variants of gated recurrent unit (GRU) neural networks[C]//Circuits and Systems, 2017 IEEE 60th International Midwest Symposium. Piscataway, NJ: IEEE Press, 2017.
[23] GOODFELLOW I, BENGIO Y, COURVILLE A. Deep learning[M]. Massachusetts: MIT Press, 2016: 203-218.
[24] 石叶楠, 郑国磊. 三种用于加工特征识别的神经网络方法综述[J]. 航空学报, 2019, 40(9): 022840.
SHI Y N, ZHENG G L. A review of three neural network methods for manufacturing feature recognition[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(9): 022840. (in Chinese)
Cited By in Cnki | Click to display the text
http://dx.doi.org/10.7527/S1000-6893.2019.23391
中国航空学会和北京航空航天大学主办。
0

文章信息

郑燕红, 邓湘金, 庞勇, 金晟毅, 姚猛, 赵志晖
ZHENG Yanhong, DENG Xiangjin, PANG Yong, JIN Shengyi, YAO Meng, ZHAO Zhihui
月球风化层钻取采样过程密实度分类研究
Research on classification of relative density in lunar regolith drilling
航空学报, 2020, 41(4): 223391.
Acta Aeronautica et Astronautica Sinica, 2020, 41(4): 223391.
http://dx.doi.org/10.7527/S1000-6893.2019.23391

文章历史

收稿日期: 2019-08-19
退修日期: 2019-09-09
录用日期: 2019-09-26
网络出版时间: 2019-10-22 09:58

相关文章

工作空间