CN103439694A - 一种验证探地雷达在月球次表层探测深度和分辨率的方法 - Google Patents
一种验证探地雷达在月球次表层探测深度和分辨率的方法 Download PDFInfo
- Publication number
- CN103439694A CN103439694A CN2013103847689A CN201310384768A CN103439694A CN 103439694 A CN103439694 A CN 103439694A CN 2013103847689 A CN2013103847689 A CN 2013103847689A CN 201310384768 A CN201310384768 A CN 201310384768A CN 103439694 A CN103439694 A CN 103439694A
- Authority
- CN
- China
- Prior art keywords
- ice
- rock
- eta
- moon
- regolith
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种验证探地雷达在月球次表层探测深度和分辨率的方法,该方法包括:选取地球表面接近月球表面电磁特性的位置,获取测月雷达在该位置的回波数据,得到测月雷达在该位置的冰层以及冰石混合层的最大穿透深度Dice_1+Dice_rock+Dice_2,以及测月雷达在冰层中的最小可分辨层厚Tice,并据此获得测月雷达在月球次表层的探测深度Drock和测月雷达在月壳浅层岩石中的厚度分辨率Trock。
Description
技术领域
本发明涉及探地雷达技术领域,特别是涉及一种验证探地雷达在月球次表层探测深度和分辨率的方法。
背景技术
测月雷达是人类首次采用表面探地雷达的方式进行月球次表层结构探测,其探测目的是获取巡视路线上的月壤厚度和月壳浅层结构。测月雷达采用超宽带无载频脉冲探地雷达技术,对巡视路线上月壳浅层岩石分层结构进行探测,要求探测深度大于100m,厚度分辨率米级。
在巡视器行进过程中,测月雷达经过发射天线向月面辐射超宽带电磁脉冲信号。信号在月壤和月壳岩石的传播过程中,如果遇到不同阻抗的介质交界面,将产生反射和散射。测月雷达接收天线接收到该反射和散射信号后,经过接收机放大、采样后获得探测数据。
测月雷达的最终性能指标将体现为探测深度和厚度分辨率两方面。为了验证测月雷达在月球次表层探测中能否满足指标要求,需要通过地面试验,评估测月雷达在月球次表层的探测深度和厚度分辨率。
由于测月雷达是我国首台月球探地雷达,也是世界首台搭载在月球表面巡视器上的探地雷达,没有成熟的地面验证方法可供采用。
发明内容
本发明的目的在于,针对月球表面探地雷达地面验证工作,提出一种验证其在月球次表层探测深度和分辨率的方法,并成功应用在了测月雷达的地面验证试验中,能够为今后的月球表面探地雷达的地面验证提供依据。
本发明提出的验证探地雷达在月球次表层探测深度和分辨率的方法包括步骤:选取地球表面接近月球表面电磁特性的位置,获取测月雷达在该位置的回波数据,得到测月雷达在该位置的冰层以及冰石混合层的最大穿透深度Dice_1+Dice_rock+Dice_2,以及测月雷达在冰层中的最小可分辨层厚Tice,测月雷达在月球次表层的探测深度Drock与在地球获得的冰层厚度Dice_1、Dice_2以及冰石混合层的厚度Dice_rock之间关系为:
其中ηair、ηice、ηice_rock、ηbase_rock、ηvacuum、ηregolith、ηrock和ηnext_layer分别表示空气、冰、冰石混合层、冰川基岩、真空、月壤、月海玄武岩和月岩中漂石或分层的波阻抗;αice、αice_rock、αregolith和αrock分别表示冰、冰石混合层、月壤和月海玄武岩的衰减系数,Dice_1表示冰川表面到冰石混合层上界面的厚度,Dice_rock表示冰石混合层厚度,Dice_2表示冰石混合层下界面到基岩的厚度,dregolith表示月壤厚度,Drock表示雷达在月球次表层的探测深度,其中Drock为待求量,Dice_1、Dice_rock和Dice_2是测量值,其他参数的为介质特性,是已知量。
优选地,
umedium,εmedium,tanδmedium分别为介质的磁导率、介电常数和损耗角正切,ω为雷达波角频率,然后代入测得的Dice_1、Dice_2和Dice_rock,以及各种介质的电磁参数,可获得测月雷达在月球次表层的探测深度Drock。
优选地,根据所述最小可分辨层厚Tice,确定测月雷达在月壳浅层岩石中的厚度分辨率Trock,二者的关系为:
其中,Trock为测月雷达在月壳浅层岩石中的厚度分辨率,Tice为测月雷达在上述地球位置的冰川中的厚度分辨率,εice为上述地球位置的旱地冰川探测区域介电常数,εrock为月海玄武岩的介电常数。
本发明提出一种验证探地雷达在月球次表层探测深度和分辨率的方法,并成功应用在了测月雷达的地面验证试验中。能够为今后的月球表面探地雷达的地面验证提供依据。
附图说明
图1为测月雷达在老虎沟12号冰川的回波图像;
图2为测月雷达在老虎沟12号冰川的回波图像,放大区域114m~156m;
图3为月球次表层分层结构模型;
图4为冰川内部分层结构模型。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明提出了一种验证探地雷达在月球次表层探测深度和分辨率的方法,并成功应用在了测月雷达的地面验证试验中。能够为今后的月球表面探地雷达的地面验证提供依据。
由探地雷达原理可知,其探测深度和厚度分辨率除了与测月雷达系统自身的技术参数相关外,还取决于月壤和月壳浅层岩石的电磁特性。因此本发明在地球上选取与虹湾地区月球次表层电磁特性接近的地方来开展验证试验,有利于探测深度和厚度分辨率的评估。
根据已有的探测结果,月球跟地球一样,也可以分成月壳、月幔和月核等层次。最外层的月壳厚约60-65km,它最上部的1-2km主要是月壤和岩石碎块。月壤主要由矿物和岩石碎屑、角砾碎屑、各种玻璃物质、陨石碎片等组成,虹湾地区月壤平均厚度约为4-5m,介电常数2.3-3.5,损耗角正切0.005-0.009;月岩介电常数6.6-8.6,损耗角正切0.009-0.016。
从上述电磁特性可以得出,虹湾地区处于干燥的环境,电磁波衰减非常小,因此本发明选取地球上旱地冰川区域来进行地面验证试验,并建立相应的分析评价模型,来评估测月雷达在月球次表层的探测深度和厚度分辨率。
本发明提出的验证探地雷达在月球次表层探测深度和分辨率的方法,包括以下步骤:首先,选取地球上接近月球特性的位置(例如甘肃省酒泉市肃北蒙古族自治县老虎沟12号冰川)进行地面验证试验,获取测月雷达在冰川中的回波数据,可以得到测月雷达在冰层以及冰石混合层的最大穿透深度Dice_1+Dice_rock+Dice_2,其中Dice_1表示冰川表面到冰石混合层上界面的厚度,Dice_rock表示冰石混合层厚度,Dice_2表示冰石混合层下界面到基岩的厚度,如图1所示,以及获取测月雷达在冰层中的最小可分辨层厚Tice,如图2所示。
其中探测深度的评估通过以下方式来实现:根据已有探测结果,建立月球次表层分层结构模型如图3所示,分为真空-月壤-月海玄武岩三层,雷达信号到达最大可探测深度处的漂石或分层时产生回波,被接收天线捕获。月壤厚度dregolith,雷达在月球次表层的探测深度Drock。根据钻取冰芯分析结果,建立冰川内部分层结构模型如图4所示,分为空气-冰-冰石混合层-冰-基岩五层,Dice_1表示冰川表面到冰石混合层上界面的厚度,Dice_rock表示冰石混合层厚度,Dice_2表示冰石混合层下界面到基岩的厚度。分析中假设各层为均匀介质,雷达波为平面波,且各层的入射角为0。
在上述模型中,通过计算雷达波传播过程中的衰减,就可以得到测月雷达在月球次表层的深度探测能力Drock与在地球旱地冰川试验中获得的冰层厚度Dice_1、Dice_2以及冰石混合层的厚度Dice_rock之间关系为:
其中,ηair、ηice、ηice_rock、ηbase_rock、ηvacuum、ηregolith、ηrock和ηnext_layer分别表示空气、冰、冰石混合层、冰川基岩、真空、月壤、月海玄武岩和月岩中漂石或分层的波阻抗;αice、αice_rock、αregolith和αrock分别表示冰、冰石混合层、月壤和月海玄武岩的衰减系数。Dice_1表示冰川表面到冰石混合层上界面的厚度,Dice_rock表示冰石混合层厚度,Dice_2表示冰石混合层下界面到基岩的厚度。dregolith表示月壤厚度,Drock表示雷达在月球次表层的探测深度。其中Drock为待求量,Dice_1、Dice_rock和Dice_2是在冰川试验中的测量值。其他的为介质特性,已知量。
其中αmedium为雷达波在该种介质中的衰减系数(例如αregolith表示在月壤中的衰减系数),在地球和月球中都是采用下面公式计算,如αregolith就是在月壤中的衰减系数,αice就是在冰中的衰减系数。
同样的ηmedium是雷达波在该种介质中的波阻抗
umedium,εmedium,tanδmedium分别为该介质的磁导率、介电常数和损耗角正切,ω为雷达波角频率。
代入测得的测月雷达在冰层以及冰石混合层的最大穿透深度Dice_1、Dice_2和Dice_rock,以及各种介质的电磁参数(参见表1),可以评估测月雷达在月球次表层的探测深度Drock。
表1:介质电磁参数
其中厚度分辨率的评估通过以下方式来实现:根据测月雷达在冰层中的最小可分辨层厚Tice,可以评估其在月壳浅层岩石中的厚度分辨率,两者的关系为:
其中,Trock为测月雷达在月壳浅层岩石中的厚度分辨率,Tice为测月雷达在冰川中的厚度分辨率,εice为冰川介电常数,εrock为月海玄武岩的介电常数。各参数的取值参见表1。如图1所示,可以得到测月雷达在冰层以及冰石混合层的最大穿透深度Dice_1=91m、Dice_2=37m、Dice_rock=24m,则将表1的参数代入式(1),可评估测月雷达在月壳浅层岩石的深度探测能力为:27-117m。
如图2所示,测月雷达在冰层中的最小可分辨层厚Tice=1.5m,根据表1的评价参数,考虑最恶劣的情况,取εrock=6.6,代入式(2)可得,测月雷达在月壳浅层岩石中的厚度分辨率不大于1m。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种验证探地雷达在月球次表层探测深度和分辨率的方法,该方法包括:选取地球表面接近月球表面电磁特性的位置,获取测月雷达在该位置的回波数据,得到测月雷达在该位置的冰层以及冰石混合层的最大穿透深度Dice_1+Dice_rock+Dice_2,以及测月雷达在冰层中的最小可分辨层厚Tice,测月雷达在月球次表层的探测深度Drock与在地球获得的冰层厚度Dice_1、Dice_2以及冰石混合层的厚度Dice_rock之间关系为:
其中ηair、ηice、ηice_rock、ηbase_rock、ηvacuum、ηregolith、ηrock和ηnext_layer分别表示空气、冰、冰石混合层、冰川基岩、真空、月壤、月海玄武岩和月岩中漂石或分层的波阻抗;αice、αice_rock、αregolith和αrock分别表示冰、冰石混合层、月壤和月海玄武岩的衰减系数,Dice_1表示冰川表面到冰石混合层上界面的厚度,Dice_rock表示冰石混合层厚度,Dice_2表示冰石混合层下界面到基岩的厚度,dregolith表示月壤厚度,Drock表示雷达在月球次表层的探测深度,其中Drock为待求量,Dice_1、Dice_rock和Dice_2是测量值,其他参数的为介质特性,是已知量,εice和εice_rock分别表示冰和冰石混合层的介电常数。
2.根据权利要求1所述的方法,其特征在于,
umedium,εmedium,tanδmedium分别为介质的磁导率、介电常数和损耗角正切,ω为雷达波角频率,然后代入测得的Dice_1、Dice_2和Dice_rock,以及各种介质的电磁参数,可获得测月雷达在月球次表层的探测深度Drock。
3.根据权利要求1所述的方法,其特征在于,根据所述最小可分辨层厚Tice,确定测月雷达在月壳浅层岩石中的厚度分辨率Trock,二者的关系为:
其中,Trock为测月雷达在月壳浅层岩石中的厚度分辨率,Tice为测月雷达在上述地球位置的冰川中的厚度分辨率,εice为上述地球位置的旱地冰川探测区域介电常数,εrock为月海玄武岩的介电常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384768.9A CN103439694B (zh) | 2013-08-29 | 2013-08-29 | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384768.9A CN103439694B (zh) | 2013-08-29 | 2013-08-29 | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103439694A true CN103439694A (zh) | 2013-12-11 |
CN103439694B CN103439694B (zh) | 2015-06-10 |
Family
ID=49693399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310384768.9A Active CN103439694B (zh) | 2013-08-29 | 2013-08-29 | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103439694B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044696A (zh) * | 2015-07-30 | 2015-11-11 | 中国科学院国家天文台 | 一种基于相关系数分析法计算测月雷达探测深度的方法 |
CN105158744A (zh) * | 2015-07-30 | 2015-12-16 | 中国科学院国家天文台 | 一种测月雷达探测深度计算方法 |
CN106019272A (zh) * | 2016-06-13 | 2016-10-12 | 中国科学院电子学研究所 | 一种采用超宽带雷达技术探测月壤分层结构的系统及方法 |
CN107153190A (zh) * | 2017-05-12 | 2017-09-12 | 中国科学院电子学研究所 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
CN113624144A (zh) * | 2021-08-06 | 2021-11-09 | 内蒙古工业大学 | 一种基于多源遥感数据估算月海玄武岩厚度的方法及系统 |
CN113884351A (zh) * | 2021-09-03 | 2022-01-04 | 中国矿业大学 | 一种基于生物固化技术的磁敏性模拟月壤的制备方法 |
CN116381821A (zh) * | 2023-06-05 | 2023-07-04 | 中国科学院地质与地球物理研究所 | 室内验证探地雷达对复杂地层探测分辨率的装置及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08105957A (ja) * | 1994-10-03 | 1996-04-23 | Natl Space Dev Agency Japan<Nasda> | レーダ装置の観測方法 |
JPH08105953A (ja) * | 1994-10-03 | 1996-04-23 | Natl Space Dev Agency Japan<Nasda> | レーダ装置 |
CN101083020A (zh) * | 2007-07-05 | 2007-12-05 | 中国航天科技集团公司第五研究院第五一○研究所 | 星球着陆探测器地面模拟试验场的建立方法 |
CN101968542A (zh) * | 2010-09-29 | 2011-02-09 | 中国科学院国家天文台 | 一种利用地球站对月球探测器进行跟踪的方法 |
CN102830392A (zh) * | 2012-08-29 | 2012-12-19 | 南京大学 | 一种月球光学遥感数据相对辐射与光度校正法 |
-
2013
- 2013-08-29 CN CN201310384768.9A patent/CN103439694B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08105957A (ja) * | 1994-10-03 | 1996-04-23 | Natl Space Dev Agency Japan<Nasda> | レーダ装置の観測方法 |
JPH08105953A (ja) * | 1994-10-03 | 1996-04-23 | Natl Space Dev Agency Japan<Nasda> | レーダ装置 |
CN101083020A (zh) * | 2007-07-05 | 2007-12-05 | 中国航天科技集团公司第五研究院第五一○研究所 | 星球着陆探测器地面模拟试验场的建立方法 |
CN101968542A (zh) * | 2010-09-29 | 2011-02-09 | 中国科学院国家天文台 | 一种利用地球站对月球探测器进行跟踪的方法 |
CN102830392A (zh) * | 2012-08-29 | 2012-12-19 | 南京大学 | 一种月球光学遥感数据相对辐射与光度校正法 |
Non-Patent Citations (1)
Title |
---|
邹立星等: "距离-多普勒法对月面目标定位精度的研究", 《测绘科学》, vol. 31, no. 2, 31 March 2006 (2006-03-31), pages 19 - 21 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044696A (zh) * | 2015-07-30 | 2015-11-11 | 中国科学院国家天文台 | 一种基于相关系数分析法计算测月雷达探测深度的方法 |
CN105158744A (zh) * | 2015-07-30 | 2015-12-16 | 中国科学院国家天文台 | 一种测月雷达探测深度计算方法 |
CN105158744B (zh) * | 2015-07-30 | 2017-08-25 | 中国科学院国家天文台 | 一种测月雷达探测深度计算方法 |
CN106019272A (zh) * | 2016-06-13 | 2016-10-12 | 中国科学院电子学研究所 | 一种采用超宽带雷达技术探测月壤分层结构的系统及方法 |
CN106019272B (zh) * | 2016-06-13 | 2019-10-22 | 中国科学院电子学研究所 | 一种采用超宽带雷达技术探测月壤分层结构的系统及方法 |
CN107153190A (zh) * | 2017-05-12 | 2017-09-12 | 中国科学院电子学研究所 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
CN107153190B (zh) * | 2017-05-12 | 2020-06-26 | 中国科学院电子学研究所 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
CN113624144A (zh) * | 2021-08-06 | 2021-11-09 | 内蒙古工业大学 | 一种基于多源遥感数据估算月海玄武岩厚度的方法及系统 |
CN113884351A (zh) * | 2021-09-03 | 2022-01-04 | 中国矿业大学 | 一种基于生物固化技术的磁敏性模拟月壤的制备方法 |
CN113884351B (zh) * | 2021-09-03 | 2023-10-27 | 中国矿业大学 | 一种基于生物固化技术的磁敏性模拟月壤的制备方法 |
CN116381821A (zh) * | 2023-06-05 | 2023-07-04 | 中国科学院地质与地球物理研究所 | 室内验证探地雷达对复杂地层探测分辨率的装置及方法 |
CN116381821B (zh) * | 2023-06-05 | 2023-08-08 | 中国科学院地质与地球物理研究所 | 室内验证探地雷达对复杂地层探测分辨率的装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103439694B (zh) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103439694B (zh) | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 | |
Hong et al. | Analyses of GPR signals for characterization of ground conditions in urban areas | |
Olsson et al. | Borehole radar applied to the characterization of hydraulically conductive fracture zones in crystalline rock 1 | |
Jeannin et al. | Multiconfiguration GPR measurements for geometric fracture characterization in limestone cliffs (Alps) | |
US8680865B2 (en) | Single well reservoir imaging apparatus and methods | |
DanieIs | Fundamentals of ground penetrating radar | |
RU2461026C1 (ru) | Способ определения геометрических характеристик трещины гидроразрыва | |
CN109521479A (zh) | 一种地下连续墙渗漏水的非开挖检测方法 | |
CN104656157B (zh) | 一种识别页岩气甜点区的方法及装置 | |
CN102759491B (zh) | 堆石体碾压密度测定方法及装置 | |
CN102012463A (zh) | 一种地层相对介电常数现场测定法 | |
CN112857698B (zh) | 一种基于声表面波检测墙体渗漏的方法 | |
Liu et al. | Asphalt pavement characterization by GPR using an air-coupled antenna array | |
Zou et al. | Study of wavelet entropy for airport pavement inspection using a multistatic ground-penetrating radar system | |
Müller et al. | Quantitative moisture measurement of road pavements using 3D noise-modulated GPR | |
Wu et al. | Mapping vertical profile of discontinuous permafrost with ground penetrating radar at Nalaikh depression, Mongolia | |
Liu | Fracture characterization using borehole radar: Numerical modeling | |
CN106198741B (zh) | 隧洞管片回填灌浆质量无损检测的工程b超方法 | |
CN113050085A (zh) | 一种超前地质预报的方法 | |
Lv et al. | Dielectric characteristics of unsaturated loess and the safety detection of the road subgrade based on GPR | |
Zhang et al. | Application of cross-borehole integrated geophysical methods for the detailed investigation of karst in urban metro construction | |
Łój et al. | Geophysical surveys and modelling for recognizing of gypsum karst | |
Aktürk et al. | Integration of electrical resistivity imaging (ERI) and ground-penetrating radar (GPR) methods to identify soil profile around Necatibey Subway Station, Ankara, Turkey | |
Liu et al. | A system for inspecting karst voids during construction of cast-in-place pile foundations | |
Silvast et al. | Onkalo edz-measurements using ground penetrating radar (gpr) method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |