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

doi:10.3969/j.issn.0253-4967.2015.03.014

地震地质 ›› 2015, Vol. 37 ›› Issue (3): 840-850.DOI: 10.3969/j.issn.0253-4967.2015.03.014

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

乔亮, 陈小斌

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

收稿日期:2014-03-26 修回日期:2014-07-27 出版日期:2015-09-20 发布日期:2015-10-20
通讯作者: 陈小斌,研究员,E-mail:cxb@pku.edu.cn
作者简介:乔亮,男,1988年生,2014年于中国地震局地质研究所获固体地球物理学专业硕士学位,主要从事大地电磁测深正反演方法研究与应用等,电话:15652803650,E-mail:qliang@mail3.sysu.edu.cn.
基金资助:
国家自然科学基金(41174058)和国家喜马拉雅计划暨中国地震局地震行业科研专项(201108001,200908001)共同资助.

AN UNSTRUCTURED GEOMETRIC MULTIGRID GENERATION ALGORITHM FOR ARBITRARY GEOELECTRIC MODEL

QIAO Liang, CHEN Xiao-bin

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

摘要/Abstract

摘要：

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

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

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

中图分类号:

P315.72⁺2

乔亮, 陈小斌. 复杂地电模型的非结构多重网格剖分算法[J]. 地震地质, 2015, 37(3): 840-850.

QIAO Liang, CHEN Xiao-bin. AN UNSTRUCTURED GEOMETRIC MULTIGRID GENERATION ALGORITHM FOR ARBITRARY GEOELECTRIC MODEL[J]. SEISMOLOGY AND GEOLOGY, 2015, 37(3): 840-850.

参考文献

陈小斌. 1999. 有限元直接迭代算法［J］. 物探化探计算技术, 2: 70—76.
CHEN Xiao-bin. 1999. Direct iterative algorithm of the finite element［J］. Computing Techniques for Geophysical and Geochemical Exploration, 2: 70—76(in Chinese).
陈小斌, 胡文宝. 2002. 有限元直接迭代算法及其在线源频率域电磁响应计算中的应用［J］. 地球物理学报, 1: 119—130.
CHEN Xiao-bin, HU Wen-bao. 2002. Direct iterative finite element(DIFE)algorithm and its application to electromagnetic response modeling of line current source［J］. Chinese Journal of Geophysics, 1: 119—130(in Chinese).
陈小斌, 赵国泽. 2009. 自动构建大地电磁二维反演的测点中心网格［J］. 地球物理学报, 6: 1564—1572.
CHEN Xiao-bin, ZHAO Guo-ze. 2009. Automatic construction of a site-centered grid(SCG)for 2D magnetotelluric inversion［J］. Chinese Journal of Geophysics, 6: 1564—1572(in Chinese).
陈小宏, 牟永光. 1995. 地震波动方程反演的多重网格方法［A］. 见: 1995年中国地球物理学会第十一届学术年会论文集.
CHEN Xiao-hong, MOU Yong-guang. 1995. Multigrid method inversion of seismic wave equation［A］. In: Proceedings of Eleventh Annual Conference of the Chinese Geophysical Society, 1995(in Chinese).
许和勇, 叶正寅. 2011. 三维非结构自适应多重网格技术［J］. 空气动力学学报, 29(3): 365—369.
XU He-yong, YE Zheng-yin. 2011. A technique of three dimensional unstructured adaptive multigrid［J］. Acta Aerodynamica Sinica, 29(3): 365—369(in Chinese).
杨钦. 2005. 限定Delaunay三角网格剖分技术［M］. 北京: 电子工业出版社.
YANG Qin. 2005. Constrained Delaunay Triangular Mesh Generation Techniques［M］. Electronic Industry Press, Beijing(in Chinese).
叶涛, 陈小斌, 严良俊. 2013. 大地电磁资料精细处理和二维反演解释技术研究(三): 构建二维反演初始模型的印模法［J］. 地球物理学报, 10: 3596—3606.
YE Tao, CHEN Xiao-bin, YAN Liang-jun. 2013. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅲ): Using the impressing method to construct starting model of 2D magnetotelluric inversion［J］. Chinese Journal of Geophysics, 10: 3596—3606(in Chinese).
赵国泽, 陈小斌, 汤吉. 2007. 中国地球电磁法新进展和发展趋势［J］. 地球物理学进展, 4: 1171—1180.
ZHAO Guo-ze, CHEN Xiao-bin, TANG Ji. 2007. Advanced geo-electromagnetic methods in China［J］. Progress in Geophysics, 4: 1171—1180(in Chinese).
赵国泽, 詹艳, 王立凤, 等. 2010. 鄂尔多斯断块地壳电性结构［J］. 地震地质, 32(3): 345—359.
ZHAO Guo-ze, ZHAN Yan, WANG Li-feng, et al. 2010. Electric structure of the crust beneath the Ordos fault block［J］. Seismology and Geology, 32(3): 345—359(in Chinese).
赵延文, 聂在平, 卢涛, 等. 2002. 基于变形玻恩迭代法的多重网格反演方法［J］. 电子科技大学学报, 31(4): 340—344.
ZHAO Yan-wen, NIE Zai-ping, LU Tao, et al. 2002. Multi-grid inversion method based on distorted Born iterative method［J］. Journal of UEST of China, 31(4): 340—344(in Chinese).
Achi Brandt. 1997. Multi-level adaptive solutions to boundary-value problems［J］. Mathematics of Computation, 31: 333—90.
Bowyer A. 1981. Computing dirichlet tessellations［J］. Comput J, 24(2): 162—166. doi: 10.1093/comjnl/24.2.162.
Cheng S W, Dey T K, Shewchuk J. 2012. Delaunay Mesh Generation［M］. Taylor & Francis.
Chew P L. 1989. Constrained Delaunay triangulations［J］. Algorithmica, 4(1-4): 97—108.
Delaunay B. 1934. Sur la sphère vide. A la mémoire de Georges Voronoï, Bulletin de l'Académie des Sciences de l'URSS［J］. Classe des sciences mathématiques et naturelles, 6: 793—800.
George P L, Borouchaki H. 1998. Delaunay Triangulation and Meshing: Application to Finite Elements［M］. Butterworth-Heinemann.
Mark de Berg, Marc van Kreveld, et al. 2008. Computational Geometry: Algorithms and Applications［M］. Springer-Verlag, ISBN 978-3-540-77973—5.
Miller G L, Talmor D, Teng S H. 1999. Optimal coarsening of unstructured meshes［J］. Journal of Algorithms, 31(1): 29—65.
Rand A. 2009. Delaunay Refinement Algorithms for Numerical Methods［M］. Carnegie Mellon University.
Rebay S. 1993. Efficient unstructured mesh generation by means of Delaunay triangulation and Bowyer-Watson algorithm［J］. Journal of Computational Physics, 106(1): 125—138.
Rineau L, Yvinec M. 2007. A generic software design for Delaunay refinement meshing［J］. Computational Geometry, 38(1-2): 100—110.
Ruppert J. 1995. A Delaunay refinement algorithm for quality 2-dimensional mesh generation［J］. Journal of Algorithms, 18(3): 548—585.
Shewchuk J R. 2002. Delaunay refinement algorithms for triangular mesh generation［J］. Computational Geometry, 22(1-3): 21—74.
Trottenberg U, Oosterlee C W, Schüller A. 2001. Multigrid［M］. Academic Press.
Watson, David F. 1981. Computing the n-dimensional Delaunay tessellation with application to Voronoi polytopes［J］. Comput J, 24(2): 167—172. doi: 10.1093/comjnl/24.2.162.
Wesseling P. 2004. An Introduction to Multigrid Methods［M］. R T Edwards, Inc.
Zhu Y. 2006. Multigrid Finite Element Methods for Electromagnetic Field Modeling［M］. Cangellaris, A.C.

[1]	赵凌强, 胡亚轩, 王庆良, 祝意青, 操聪, 李仲巍, 綦伟, 汶宇龙. 吉林龙岗火山区深部电性结构特征分析[J]. 地震地质, 2022, 44(4): 845-858.
[2]	董泽义, 汤吉, 赵国泽, 陈小斌, 崔腾发, 韩冰, 姜峰, 王立凤. 首都圈极低频电磁台网区地下电性结构探测[J]. 地震地质, 2022, 44(3): 649-668.
[3]	赵凌强, 詹艳, 王庆良, 孙翔宇, 韩静, 操聪, 张松, 蔡妍. 大地电磁数据揭示的1303年洪洞8级地震区精细结构[J]. 地震地质, 2022, 44(3): 686-700.
[4]	韩静, 詹艳, 孙翔宇, 赵国泽, 刘雪华, 包雨鑫, 孙建宝, 彭远黔. 强电磁干扰环境下的大地电磁数据特征及处理[J]. 地震地质, 2022, 44(3): 736-752.
[5]	张赟昀, 王培杰, 陈小斌, 詹艳, 韩冰, 王立凤, 赵国泽. 强干扰环境下的大地电磁时间序列处理过程[J]. 地震地质, 2022, 44(3): 786-801.
[6]	刘钟尹, 陈小斌, 蔡军涛, 崔腾发, 赵国泽, 汤吉, 欧阳飚. 大地电磁三维反演云计算系统toPeak的设计与实现[J]. 地震地质, 2022, 44(3): 802-820.
[7]	阮帅, 汤吉, 董泽义, 王立凤, 邓琰, 韩冰. 基于三维大地电磁AR-QN反演的长白山天池火山区电性结构[J]. 地震地质, 2020, 42(6): 1282-1300.
[8]	孙翔宇, 詹艳, 赵国泽, 赵凌强, 邓琰, 胡亚轩, 胡久常, 向小娟. 琼东北马鞍岭-雷虎岭火山区深部岩浆系统大地电磁三维探测[J]. 地震地质, 2020, 42(3): 640-653.
[9]	孙翔宇, 詹艳, 赵凌强, 陈小斌, 李陈侠, 孙建宝, 韩静, 崔腾发. 东昆仑断裂带东端和2017年九寨沟7.0级地震区深部电性结构探测[J]. 地震地质, 2020, 42(1): 182-197.
[10]	姜峰, 陈小斌, 董泽义, 崔腾发, 刘钟尹, 王培杰. 利用单剖面大地电磁三维反演识别沙德和玉农希断裂[J]. 地震地质, 2019, 41(6): 1444-1463.
[11]	张继红, 赵国泽, 董泽义, 王立凤, 韩冰, 王庆林, 唐廷梅, 王梅. 郯庐断裂带安丘、莒县电磁台地壳电性结构研究[J]. 地震地质, 2019, 41(5): 1239-1253.
[12]	赵凌强, 詹艳, 王庆良, 孙翔宇, 杨皓, 陈小斌. 1954年甘肃民勤7级地震区深部电性结构特征及地震构造环境研究[J]. 地震地质, 2018, 40(3): 552-565.
[13]	韩江涛, 王天琪, 刘文玉, 刘国兴, 韩松, 刘立家. 阿尔山火山群深部“拱桥式”岩浆系统及其稳定性分析[J]. 地震地质, 2018, 40(3): 590-610.
[14]	周耀明, 朱文斌, 陈正乐, 朱炳玉, 薛峰. 准噶尔盆地克-百断裂带火山岩分布特征的重磁资料解释[J]. 地震地质, 2018, 40(3): 641-655.
[15]	翁爱华, 李建平, 范小平, 李斯睿, 韩江涛, 李大俊, 李亚彬, 赵祥阳, 唐裕. 大地电磁测深揭示的1668年郯城8.5级地震震中地壳精细结构[J]. 地震地质, 2018, 40(2): 396-409.

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

AN UNSTRUCTURED GEOMETRIC MULTIGRID GENERATION ALGORITHM FOR ARBITRARY GEOELECTRIC MODEL

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价