三维稳态地温场全信息数值模拟

查看详情 浏览次数:1
单位中南大学有色金属成矿预测与地质环境监测教育部重点实验室;中南大学地球科学与信息物理学院;中国海洋大学海洋地球科学学院;
来源地球物理学报
出版年2025
期号03
摘要
针对常规的数值模拟方法存在计算量大、边界条件近似和难以准确地模拟三维偏微分方程表征的物理信息等问题,本文提出了任意起伏地形下的一种高效、高精度的全信息数值模拟方法.该方法沿水平方向(-∞,+∞)进行二维全信息傅里叶变换,将三维偏微分方程转化为不同波数下的一维常微分方程进行求解,减少了计算量及存储需求,不同波数之间常微分方程相互独立,具有高度并行性.保留垂向为空间域并引入行波分解,消除了上下边界对模拟区域的影响,具有严格的上下边界条件,且每个垂向深度对应的水平网格任意,具有分辨率高,剖分灵活等特点.水平方向采用二维全信息傅里叶变换,获得异常场的正确波数谱信息,再对异常场所有的波数谱信息进行二维全信息逆变换,从而得到空间域的异常场.该方法在水平方向不存在截断效应,相对于标准FFT和Gauss-FFT,没有丢掉任何波数谱信息;每个单元利用一维形函数的高精度解和傅里叶变换的高效性,实现了高效、快速和准确地模拟三维偏微分方程描述的物理信息.本文设计棱柱体,对比数值解与解析解结果,验证了本文算法的正确性,分析异常场的波数谱分布特征,总结波数选取规律,同时对比本文算法与有限单元法、Gauss-FFT的计算精度与效率.设计高、低热导组合模型,进一步讨论边界条件的影响.最后设计任意起伏地形,说明本文算法对任意复杂地形下的适应性,为任意复杂条件大规模、高效、精细化反演成像、定量解释和人机相互解释奠定了基础.

@ 2023 版权所有 中国地质图书馆 (中国地质调查局地学文献中心)

京ICP备 05064591号 京公网安备11010802017129号

建议浏览器: 火狐、谷歌、微软 Edge、不支持 IE