SEISMOLOGY AND GEOLOGY ›› 2015, Vol. 37 ›› Issue (3): 840-850.DOI: 10.3969/j.issn.0253-4967.2015.03.014

Previous Articles     Next Articles

AN UNSTRUCTURED GEOMETRIC MULTIGRID GENERATION ALGORITHM FOR ARBITRARY GEOELECTRIC MODEL

QIAO Liang, CHEN Xiao-bin   

  1. Institute of Geology, China Earthquake Administration, Beijing 100029, China
  • Received:2014-03-26 Revised:2014-07-27 Online:2015-09-20 Published:2015-10-20

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

乔亮, 陈小斌   

  1. 中国地震局地质研究所, 地震动力学国家重点实验室, 北京 100029
  • 通讯作者: 陈小斌,研究员,E-mail:cxb@pku.edu.cn
  • 作者简介:乔亮,男,1988年生,2014年于中国地震局地质研究所获固体地球物理学专业硕士学位,主要从事大地电磁测深正反演方法研究与应用等,电话:15652803650,E-mail:qliang@mail3.sysu.edu.cn.
  • 基金资助:

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

Abstract:

With the growing number of observational data, increasingly complex geoelectric model, and high-dimensional magnetotelluric method, the efficient and stable forward and inversion technologies become more and more required. Multigrid method enables the asymptotically optimal approximate numerical solution to elliptic partial differential equations including the MT 2D forward modeling problem. Currently, the coarsening and refining operations are based on regular grid which is easy to understand and implement. However, when the geometry of a region is complex or the region needs local refinements, regular grid becomes not well applicable. Unstructured triangular grid has better geometrical adaptability than structured grid. The Delaunay triangulation ensures that each triangle is generated to be Delaunay. This property is particularly suitable for numerical interpolation and FEM calculation. Besides, the Delaunay triangulation has a strong mathematical theory foundation which makes it very convenient for local refinement and sparseness. The key to solve the complicated geoelectric model using multigrid method is the automatic generation of unstructured multi-level grids. In this paper, we present an algorithm for generating unstructured grid for Multigrid method based on Delaunay triangulation. The algorithm will automatically generate the coarse and fine grid for complex input region, and all the triangle elements abide by the Delaunay criteria which ensure numerical accuracy and high convergence rate.

Key words: magnetotelluric, arbitrary geoelectric model, multigrid method, Delaunay triangulation, mesh refinement

摘要:

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

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

CLC Number: