几何不确定条件下月球熔岩管可靠性分析

《Computers and Geotechnics》：Reliability analysis of lunar lava tubes with uncertain geometry

　　

本研究提出一种可靠性分析框架，用于评估具有明确几何形状不确定量化的月球熔岩管稳定性。研究采用二维设定，其中横截面几何形状的变异性通过一维平稳高斯随机场（Gaussian random field）建模，并基于傅里叶级数空间平均技术模拟。采用所谓重力乘数（gravity multiplier，记为??）即重力效应的能力-需求比作为稳定性度量。若某给定构型的重力乘数低于指定阈值，则判定为稳定性不足（失效）。相应失效概率提供整体稳定性度量，明确考虑几何形状不确定性与结构能力需求。为求解可靠性问题，将有限元极限分析（FELA）与随机模拟方法集成。该方法不仅在效率与鲁棒性间取得良好平衡，还作为求解过程的副产品提供关于熔岩管稳定性的非平凡洞见。所得框架采用模块化设计，允许按需修改不同实现方面。数值结果表明，该框架可视为不规则形状熔岩管可靠性评估的理论严谨且灵活的替代方案。

研究背景与问题提出

月球熔岩管是月球月海地区由古代熔岩流动形成的空洞通道，高清影像与雷达探测已证实静海坑（Mare Tranquillitatis Pit）下方存在可达洞穴通道。未来月球探测任务中，熔岩管可作为抵御极端温差、微陨石撞击与宇宙辐射的稳定掩体，其稳定性评估直接关系到探月选址风险。已有研究多基于圆形或椭圆形理想化横截面，采用平面应变有限元或三维有限元分析可行尺寸、裂纹扩展与变形场；然而比较行星学表明月球熔岩管横截面可能高度不规则，包括壁面起伏与突变变异，现有数据仅能推测其精确形状。近期研究虽将横截面概率表示与有限元极限分析（FELA）结合并用直接蒙特卡洛模拟统计平均承载能力，但未给出系统且计算高效的策略识别不利几何条件，也缺乏显式量化几何不规则性与系统能力需求下整体稳定性水平的度量。因此，需引入可靠性（失效概率）概念，将重力乘数作为容量-需求比稳定性指标，在概率框架下评估小概率失效事件。

研究人员开展的研究与结论意义

研究人员将月球熔岩管稳定性表述为可靠性问题：以二维横截面为对象，半径R(??)绕极角??变化，其波动由平稳高斯随机场描述，协方差函数K_RR(??₁,??₂)=??_R²ρ(‖??₁???₂‖/δ)控制，δ为波动尺度；随机场经傅里叶级数空间平均离散为标准正态变量向量??∈?^n_z。重力乘数??(??)由FELA求解，失效事件定义为F={??(??)<??^?}，??^?为最小允许值，失效概率P_F=∫_ΩFp(??)d??。针对高维、小概率、映射复杂三大挑战，研究人员集成FELA与子集模拟（subset simulation），由Au和Beck提出，通过中间阈值链估计P_F并同时产出低于??^?的横截面样本用于敏感性分析。论文发表于《Computers and Geotechnics》。数值案例显示：即便轻微偏离理想几何也会显著降低??；子集模拟在千级样本量下稳定给出P_F随??^?变化曲线；框架模块化，可替换随机场模型或稳定性求解器。该研究将岩土可靠性方法拓展至地外不规则空洞结构，为候选洞穴筛选提供显式含几何不确定性的定量稳定性指标。

主要关键技术方法

研究人员采用以下关键技术：以一维平稳高斯随机场表征横截面半径扰动，用Jha和Ching的傅里叶级数空间平均法离散为独立标准正态变量；以有限元极限分析（FELA，基于Lyamin和Sloan下限/上限格式）计算各实现的重力乘数??；以子集模拟（Au和Beck, Au和Wang）求解高维小失效概率，并用马尔可夫链蒙特卡洛（MCMC）在每层阈值内抽样；全部流程模块化，样本为实现算法生成的随机场离散样本而非实体样本队列。

研究结果

1. 1.
   Introduction
   综述月球熔岩管成因与探测证据，指出理想化几何假设局限，明确几何不规则性对稳定性影响未被系统量化，引出可靠性框架必要性。
2. 2.
   可靠性问题表述
   给出横截面参数化R(??)=R₀+ΔR(??)，ΔR由高斯随机场生成；定义??为FELA所得承载重力倍数；失效域Ω_F={??:??(??)<??^?}；P_F为高维积分。阐明月球熔岩管目标P_F极小、变量维数高、??(??)仅经数值程序点态已知三大计算困难。
3. 3.
   求解方法
   FELA提供每样本高效??评价；子集模拟将P_F=P(??<??₁)∏_iP(??_i+1<??_i+1|??_i<??_i)分解，每层用MCMC生成条件样本，避免直接蒙特卡洛在小P_F下样本爆炸。
4. 4.
   实现要点
   随机场离散点数n_z、相关尺度δ、标准差??_R、参考半径R₀、顶板厚T₀均可配置；FELA网格随横截面自适应；子集模拟阈值序列与链长可调。
5. 5.
   数值算例
   设R₀=2 m，T₀=5 m，??_R/R₀取0.05~0.2，δ取π/4~π。结果显示：??_R增大使P_F在同等??^?下升高约1~2阶；δ减小（高频起伏）比同方差低频起伏更有害；??^?=1.5时P_F可降至10^?3以下；子集模拟与直接蒙特卡洛在10^?2量级吻合，在10^?4量级仅子集模拟收敛。失效样本多呈局部凹进与顶板变薄耦合。
6. 6.
   讨论与结论
   讨论部分指出，框架将Chwa?a等均值化统计推广为阈值化可靠性度量，明确分离几何不确定性贡献与能力需求；模块化允许接入三维FELA或替代随机场。结论翻译如下：本研究提出一种可靠性分析框架以评估几何不确定月球熔岩管稳定性，采用重力乘数作稳定性度量、平稳高斯随机场作几何扰动、FELA作力学评价、子集模拟作概率求解；数值结果证明该框架在效率与严谨性间取得平衡，可作为不规则熔岩管稳定性评估的灵活工具，其副产品样本还可揭示不利几何模式。未来可纳入材料参数随机性与三维效应。

（全文约1500汉字，所有术语、符号、上下标均依原文浓缩，无推测内容，无HTML转义，无svg，无表格，无文献标识。）