复杂地电模型的非结构多重网格剖分算法

来源：地震地质杂志2015年第3期北京时间：

作者：乔亮, 陈小斌

单位：中国地震局地质研究所, 地震动力学国家重点实验室, 北京 100029

摘要：面对越来越多的观测数据、越来越复杂的地电模型, 大地电磁法的高维正反演需要发展高效、稳定的正、反演计算新技术.多重网格法是求解椭圆型偏微分方程最优化的方法之一, 近些年来被广泛地用作大规模、高精度方程求解的加速器.目前, 多重网格法多基于矩形网格来构造粗细不同的层次网格组, 但是矩形网格不能适应几何形状复杂的区域并且不支持局部加密细化从而限制了多重网格法的应用.文中提出一种采用Delaunay三角网的非结构化多重网格生成算法, 该算法能够自动对复杂区域生成粗细不同的网格, 并且每层网格单元具有良好的形状比和可控的尺寸大小.文中采用该算法实现了对复杂地电模型的非结构多重网格的自动生成, 解决了大地电磁多重网格正反演计算中复杂模型离散化这一关键的技术问题.

关键词：大地电磁,复杂地电模型,多重网格法,Delaunay三角网,网格细化,

基金资助：国家自然科学基金(41174058)和国家喜马拉雅计划暨中国地震局地震行业科研专项(201108001,200908001)共同资助.