为表征冻土地层天然存在的三维空间变异性并阐明冻结帷幕温度特征值演化过程与统计规律，将冻土地层热学参数建模为三维空间随机场，提出用于不规则三维空间随机场离散的四面体单元局部平均法；引进体积坐标变换和高斯数值积分呈现了不规则三维四面体随机场单元协方差矩阵的解析计算方法和数值计算方法；构建考虑冻土地层三维空间变异性的冻结帷幕水热耦合随机分析模型，获得温度特征值动态演化过程与统计规律，讨论表征三维空间变异性的随机场相关结构模式对冻结帷幕温度特征值的影响。结果表明：提出的三维随机场四面体单元离散法能与有限元四面体单元离散法完美结合，随机场网格与有限元网格的对应关系清晰，易于随机有限元程序的编制；冻土地层热学参数的空间变异性对冻结帷幕温度特征值有显著影响，其中导热系数空间变异性影响最为强烈；不同三维随机场相关结构模式对冻结帷幕温度特征值影响程度存在差异，其中三角型模式影响最大，高斯型为中间值，指数型影响最小。

来源平台：岩石力学与工程学报