一种多尺度非均匀月壤层内离散随机介质建模方法
【技术领域】
[0001] 本发明涉及月球行星科学探测领域，尤其涉及一种多尺度非均匀月壤层内离散随 机介质建模方法。
【背景技术】
[0002] 20世纪60年代人类第一次登上月球之后，人类对月球的内部构造探索就从未停 止。根据天然和人工月震所提供的资料表明，月球跟地球一样，也可以分成月壳、月幔和月 核等层次。最外层的月壳厚约60~65km，它最上部的1~2km主要是月壤和岩石碎块。通 过对Apol 1〇样品的分析表明，月壤主要由矿物和岩石碎肩、角砾碎肩、各种玻璃物质、陨石 碎片等组成，月壤密度与深度成双曲线关系，介电常数为2. 3~3. 5,损耗角正切为0. 005~ 0.009;月球高地岩石主要由斜长岩构成，而月海则主要由玄武岩组成，月海玄武岩介电常 数为6. 6~8. 6,损耗角正切为0. 009~0. 016。
[0003]1972年美国Apollo 17搭载的ALSE雷达探测仪，用于探测月球次表层地质结构。 2007年日本SELENE搭载的LRS，其原理与ALSE相同。LRS探测了大部分月海地区的次表面 分层结构。2013年12月2日，我国成功发射了嫦娥三号，其由着陆器与巡视器组成，在巡视 器上搭载了一部测月雷达（LPR，Lunar Penetrating Radar)，其是一种高分辨率月球表面 穿透成像雷达。LPR由两个通道组成：一通道中心频率为60MHz，厚度分辨率为米级，用于探 测巡视路上次表层结构；二通道中心频率为500MHz，厚度分辨率小于30cm，用于探测巡视 路上月壤结构以及厚度。测月雷达回波信号图反映了月壤内部介质的阻抗差异，是介电常 数、电导率和磁导率差异的综合体现，另外还由于介质的非均匀性、多次反射、环境反射、衍 射等的影响，不可避免的存在多解性与复杂性。在地球物理方法中，通常的研宄思路是建立 正演模型，进行数值模拟计算，用于解释测月雷达回波信号的地质结构信息。
[0004] 然而，在月球探测过程中，反演月壤厚度时，建立的月壤模型都是均匀且理想的结 构模型，因此，现有技术的月壤模型难以解释测月雷达回波信号图所反映的地质结构信息。
【发明内容】
[0005] 本发明目的在于提出一种多尺度非均匀月壤层内离散随机介质建模方法，该方法 针对测月雷达回波信号图所反映的月壤内部介质的阻抗差异以及介质的非均匀性、多次反 射、环境反射、衍射等的影响，建立能够解释测月雷达回波信号图所反映的地质结构信息的 月壤模型。
[0006] 根据本发明的一方面，一种多尺度非均匀月壤层内离散随机介质建模方法包括以 下步骤：步骤1，确定月壤层内随机介质的各参数值，并构建自相关函数；步骤2,通过对所 构建的自相关函数进行傅里叶变换，求出空间随机扰动函数的功率谱密度函数；步骤3,针 对所求出的功率谱密度函数加入随机相位函数，计算随机功率谱函数；步骤4,通过对计算 出的随机功率谱函数做逆傅里叶变换，得到随机介质的空间随机扰动函数；步骤5,规范化 空间随机扰动函数，得到离散随机介质模型；以及步骤6,根据月岩采样样品数据以及月壤 层内部月岩颗粒分布特征，在离散随机介质模型内添加月岩颗粒，从而建立多尺度非均匀 月壤层内离散随机介质模型。
[0007] 优选地，在假设月壤层内的随机介质模型为平稳随机过程的前提下，多尺度非均 匀月壤层内随机介质模型的介电常数由式（*)表述：
[0008]
(*)
[0009] 其中，F = A +A，F是矢径，x，z为空间点坐标值，ex，％分别是沿x， z坐标轴 正向的单位矢量；％为背景大尺度扰动，以均值表示；为模型标准差；为小尺度 随机扰动，其空间分布特征服从自相关函数。
[0010] 优选地，在步骤1中，构建自相关函数式（1)，
[0011]
(1)
[0012] 其中，a, b分别为随机介质水平自相关长度与垂直自相关长度；0为自相关角度； r为粗糙度因子，0 < r < 1。
[0013] 优选地，当r = 0时式（1)即为高斯自相关函数，当r= 1时式（1)即为指数型自 相关函数，当〇〈r〈l时式（1)即为混合型自相关函数。
[0014] 优选地，自相关函数为椭圆自相关函数。
[0015] 优选地，在步骤2中，求出空间随机扰动函数的功率谱密度函数R(kx，k z)，其 等于自相关函数的傅里叶变换：
[0016]
(2)
[0017] 其中，左+ ，是是矢社，kx，k z为空间点坐标值，e x，ez是沿k x，匕坐标轴正 向的单位矢量。
[0018] 优选地，在步骤3中，针对功率谱密度函数:加入随机相位函数舛F)，g卩，在值 域[0,2 31)区间内独立、均匀分布的二维随机序列，并由式（3)计算随机功率谱函数
[0019]
(3)
[0020] 其中，^二f是矢径，kx，kz为空间点坐标值，e x，ez分别是沿k x，kz坐标 轴正向的单位矢量。
[0021] 优选地，在步骤4中，对随机功率谱函数r(F).进行逆傅里叶变换，即可得到随机介 质的空间随机扰动函数心:
[0022]
(4)〇
[0023]优选地，在步骤5中，规范化空间随机扰动函数，并将式⑷带入式（*)，便得 到了离散随机介质模型。
[0024]优选地，在步骤6中，月岩采样样品数据以及月壤层内部月岩颗粒分布特征采用 Apol 1〇登月探测获取的月岩采样样品数据以及月壤层内部月岩颗粒分布特征。
[0025] 本发明基于统计学概念，运用随机等效介质模型理论，建立了多尺度非均匀月壤 层内离散随机介质建模方法，这更加符合月壤内部地质结构的模型。该方法针对测月雷达 回波信号图所反映的月壤内部介质的阻抗差异以及介质的非均匀性、多次反射、环境反射、 衍射等的影响，建立能够解释测月雷达回波信号图所反映的地质结构信息的月壤模型。
【附图说明】
[0026] 图1是根据本发明的一种多尺度非均匀月壤层内离散随机介质建模方法的流程 图。
[0027] 图2是根据本发明的多尺度非均匀月壤层内离散随机介质建模方法的优选实施 例中的自相关函数的示例性示意图。
[0028] 图3是根据本发明的多尺度非均匀月壤层内离散随机介质建模方法的优选实施 例中的均匀随机数的示例性示意图。
[0029] 图4是根据本发明的多尺度非均匀月壤层内离散随机介质建模方法的优选实施 例中的空间随机扰动函数的示例性示意图。
[0030] 图5是根据本发明的多尺度非均匀月壤层内离散随机介质建模方法的优选实施 例中的随机介质模型的示例性示意图。
[0031] 图6是根据本发明的多尺度非均匀月壤层内离散随机介质建模方法建立的多尺 度非均匀月壤层内离散随机介质模型的一个优选实施例。
[0032] 图7是依据本发明的多尺度非均匀月壤层内离散随机介质建模方法的多尺度非 均匀月壤层内离散随机介质建模流程的一个优选实施例。
【具体实施方式】
[0033] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照 附图，对本发明进一步详细说明。
[0034] 图1是根据本发明的一种多尺度非均匀月壤层内离散随机介质建模方法的流程 图。如图1所示，本发明的多尺度非均匀月壤层内离散随机介质建模方法包括以下步骤：
[0035] 步骤1，确定月壤层内随机介质的各参数值，并构建自相关函数；
[0036] 步骤2,通过对所构建的自相关函数进行傅里叶变换，求出空间随机扰动函数的功 率谱密度函数；
[0037] 步骤3,针对所求出的功率谱密度函数加入随机相位函数，计算随机功率谱函数；
[0038] 步骤4,通过对计算出的随机功率谱函数做逆傅里叶变换，得到随机介质的空间随 机扰动函数；
[0039] 步骤5,规范化空间随机扰动函数得到离散随机介质模型；以及
[0040] 步骤6，根据月岩采样样品数据以及月壤层内部月岩颗粒分布特征，在离散随机介 质模型内添加月岩颗粒，从而建立多尺度非均匀月壤层内离散随机介质模型。
[0041] 在本发明的多尺度非均匀月壤层内离散随机介质建模方法中，在假设月壤层内的 随机介质模型为平稳随机过程的前提下，多尺度非均匀月壤层内随机介质模型的介电常数 由式（*)表述：
[0042]
(*)
[0043] 其中，r =X4 + K，7是矢径，x，z为空间点坐标值，ex，别是沿x，z坐标轴 正向的单位矢量；％为背景大尺度扰动，以均值表示；A3:为模型标准差为小尺度 随机扰动（空间随机扰动函数），其空间分布特征服从自相关函数。
[0044] 在实际建模过程中，计算机无法模拟无限大的区域，需要在有限的范围内将以上 连续随机介质模型进行离散化，并由离散间隔确定模型的精细程度。因此，可以通过根据本 发明的多尺度非均匀月壤层内离散随机介质建模方法来进行建模。
[0045] 在步骤1中，确定月壤层内随机介质的各参数值，并构建自相关函数。例如，将 随机介质模型区域设为2mX2m，均匀网格划分间隔为A = 0.01m，即划分点数为MXN = 200X200;r = 0.5，a = b = 0.01，0 =5°。构建自相关函数式（1)，如图2和图7所示；
[0046]
(1)
[0047] 其中，a，b分别为随机介质水平自相关长度与垂直自相关长度；0为自相关角度； r为粗糙度因