一种地外天体表面热惯量反演与验证方法技术

本发明专利技术一种地外天体表面热惯量反演与验证方法，包括：利用目标地物不同深度处的内部测温数据计算目标地物热扩散率，结合目标地物热容的实验测量值，获得目标地物热惯量的计算值；利用热成像仪对目标地物的表面测温数据，根据热惯量理论计算模型，求解目标地物的热惯量；利用由内部测温数据计算得到的热惯量拟合修正热惯量反演模型，提高基于表面测温数据的热惯量测算精度。本发明专利技术旨在解决目前基于热成像仪表面测温数据的热惯量反演方法的测算精度问题，为融合热特性的地外天体表面危险地形识别和可通过性分析提供技术储备。识别和可通过性分析提供技术储备。识别和可通过性分析提供技术储备。

【技术实现步骤摘要】
一种地外天体表面热惯量反演与验证方法


[0001]本专利技术属地外天体表面探测
，尤其涉及一种地外天体表面热惯量反演与验证方法。

技术介绍

[0002]已有地外天体表面探测中，对地形的识别和可通过性分析均基于可见光图像获取的表面几何特征。但单纯的可见光视觉测量易受光照和表面纹理的影响，造成对表面几何特征不明显危险地形(如易打滑/沉陷沙地)可通过性的误判，影响地外星表探测器通行的安全性。
[0003]热惯量是物质热特性的一种综合量度，反映了物质与周围环境能量交换的能力，也代表了表面物质对抗外界温度变化的能力，可以近似认为其代表表面物质的粒度。热惯量较高的物质昼夜温差较小，如石块(热惯量通常大于1200J
·
m
‑2·
K
‑1·
s
‑
1/2
)；热惯量较小的物质昼夜温差较大，如灰尘(热惯量通常小于150J
·
m
‑2·
K
‑1·
s
‑
1/2
)。在特定的表面气压条件下，由颗粒构成的地形的热惯量是颗粒大小、密度以及聚合程度的函数。而火星车行驶过程中，在疏松沙质地形中的可通过性与颗粒大小、密度及聚合程度强相关，表明车轮在行进中的深陷、滑移情况都与地形的热惯量相关。
[0004]利用火星表面热特性反演其内部地质特征，用于火星表面危险地形判断和火星车运动轨迹规划，是近年火星表面探测研究前沿。美国NASA、加州理工等进行了初步研究，但由于热惯量难以直接测量，目前均利用轨道器或车载热成像仪的测温数据计算热惯量，存在的主要问题是缺少对测量和计算效果的评估，热惯量的测算精度不能保证，影响基于热惯量对危险地形及可通过性判断的准确性。

技术实现思路

[0005]本专利技术的技术解决问题：克服现有技术的不足，提出了一种地外天体表面热惯量反演与验证方法，通过构建实验获取内部温度测量数据，并计算得到目标地物的热惯量，与利用热惯量理论计算模型和热成像仪表面测温数据的热惯量反演值进行对比，评估基于热成像仪测量数据的热惯量反演效果，并对反演模型进行校正，保证热惯量的计算精度，为地外天体表面危险地形识别及可通过性判断提供技术储备。
[0006]本专利技术的技术解决方案是：
[0007]一种地外天体表面热惯量反演与验证方法，包括：
[0008]选取N种不同材质地物作为目标，利用目标地物测点不同深度处的内部测温数据计算目标地物热扩散率D
h
，结合目标地物热容的实验测量值C
h
，获得每种目标地物热惯量的计算值P
T
；
[0009]利用目标地物的表面测温数据，获得每种目标地物的热惯量的估测值P
C
；
[0010]利用每种目标地物的热惯量计算值P
T
和每种目标地物的热惯量的估测值P
C
，拟合修正热惯量反演模型，获得修正后的热惯量反演模型。
[0011]优选地，所述获得目标地物热惯量的计算值P
T
的方法，具体为：
[0012]选取地表下方两处不同位置的深度(z1、z2)，获得两处位置在一个光周期t
p
内的i个温度值；
[0013]利用两处位置在一个光周期t
p
内的i个温度值，拟合获得目标地物热扩散率D
h
；
[0014]通过体积热容的测量实验，获得目标地物的热容C
h
；
[0015]利用获得的目标地物热扩散率D
h
和获得的目标地物热容C
h
，得到目标地物热惯量的计算值P
T
。
[0016]优选地，利用T型热电偶作为温度测量装置，获得两处位置在一个光周期t
p
内的i个温度值；
[0017]两处位置之间的深度差值不小于0.1米，且两处位置中深度最浅的位置深度不大于0.1米；或者，一个光周期t
p
内，两处位置的温差最大值不小于10℃。
[0018]优选地，拟合获得目标地物热扩散率D
h
的方法，具体为：
[0019]获得两处位置在一个光周期t
p
内的24个温度值，分别记温度测量数据为和
[0020]将温度测量数据和分为6组，每组对应z1和z2位置的4个时刻的温度测量数据；
[0021]则目标地物热扩散率D
h
具体为：
[0022][0023][0024]式中，ω＝2π/t
p
，k＝0,1,2,3,4,5。
[0025]优选地，获得目标地物的热容C
h
的方法，具体为：
[0026]选取体积为V的密闭绝热容器，该容器内装满目标地物，并且缠绕加热丝，通过功率为P
C
的电加热使容器内的目标地物升温，测量加热时间Δt内目标地物的升温值ΔT
C
，并计算目标地物的热容C
h
如下：
[0027][0028]优选地，所述目标地物热惯量的计算值P
T
，具体为：
[0029][0030]优选地，获得目标地物的热惯量的估测值P
C
的方法，具体为：
[0031]21)计算热惯量理论计算模型中傅立叶级数的系数(A1，A2)；
[0032]22)根据表面温度测量数据，计算热惯量迭代求解的初值B0；
[0033]23)求解目标地物的热惯量的估测值P
C
。
[0034]优选地，热惯量理论计算模型中傅立叶级数的系数(A1，A2)，具体为：
[0035][0036][0037]ψ＝arccos(tanδ
·
tanα)
[0038]其中，α为当地纬度，δ为太阳赤纬。
[0039]优选地，获取目标地物测点表面位置在一个光周期t
p
内的表面温度测量数据T
i
,i＝0,1,2,23，筛选该光周期内的最高温度值得到最高温时刻t
max
；
[0040]从一个光周期t
p
内的表面温度测量数据T
i
,i＝0,1,2,
…
23中，任意选取两个不相同的温度测量值，获得两个不相同的温度测量值对应的时刻为t1,t2，计算得到两个不相同的温度测量值对应温度差的绝对值
[0041]计算热惯量迭代求解的初值B0，具体为：
[0042][0043]ω＝2π/t
p
[0044][0045][0046]式中，A为目标测点处的地表反照率，C
t
为大气透过率，S0为太阳常数；
[0047]根据如下以P
C
,B为未知变量的非线性方程，利用热惯量迭代求解的初值B0，采用牛顿迭代法求解得到目标地物热惯本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种地外天体表面热惯量反演与验证方法，其特征在于，包括：选取N种不同材质地物作为目标，利用目标地物测点不同深度处的内部测温数据计算目标地物热扩散率D
h
，结合目标地物热容的实验测量值C
h
，获得每种目标地物热惯量的计算值P
T
；利用目标地物的表面测温数据，获得每种目标地物的热惯量的估测值P
C
；利用每种目标地物的热惯量计算值P
T
和每种目标地物的热惯量的估测值P
C
，拟合修正热惯量反演模型，获得修正后的热惯量反演模型。2.根据权利要求1所述的一种地外天体表面热惯量反演与验证方法，其特征在于，所述获得目标地物热惯量的计算值P
T
的方法，具体为：选取地表下方两处不同位置的深度(z1、z2)，获得两处位置在一个光周期t
p
内的i个温度值；利用两处位置在一个光周期t
p
内的i个温度值，拟合获得目标地物热扩散率D
h
；通过体积热容的测量实验，获得目标地物的热容C
h
；利用获得的目标地物热扩散率D
h
和获得的目标地物热容C
h
，得到目标地物热惯量的计算值P
T
。3.根据权利要求2所述的一种地外天体表面热惯量反演与验证方法，其特征在于，利用T型热电偶作为温度测量装置，获得两处位置在一个光周期t
p
内的i个温度值；两处位置之间的深度差值不小于0.1米，且两处位置中深度最浅的位置深度不大于0.1米；或者，一个光周期t
p
内，两处位置的温差最大值不小于10℃。4.根据权利要求2所述的一种地外天体表面热惯量反演与验证方法，其特征在于，拟合获得目标地物热扩散率D
h
的方法，具体为：获得两处位置在一个光周期t
p
内的24个温度值，分别记温度测量数据为和将温度测量数据和分为6组，每组对应z1和z2位置的4个时刻的温度测量数据；则目标地物热扩散率D
h
具体为：具体为：式中，ω＝2π/t
p
，k＝0,1,2,3,4,5。5.根据权利要求2所述的一种地外天体表面热惯量反演与验证方法，其特征在于，获得目标地物的热容C
h
的方法，具体为：选取体积为V的密闭绝热容器，该容器内装满目标地物，并且缠绕加热丝，通过功率为P
C
的电加热使容器内的目标地物升温，测量加热时间Δt内目标地物的升温值ΔT
C
，...

【专利技术属性】
技术研发人员：邢琰，姜甜甜，贾永，魏春岭，胡海霞，胡海涛，
申请(专利权)人：北京控制工程研究所，
类型：发明
国别省市：