基于收缩二叉树结构网格剖分的大地电磁二维有限单元法正演研究.pdf

年 月 第 卷 第 期 ２０１３ ６ ４８ ３ 　 　 　 ·非地震 · 文章编号： （ ） １０００７２１０２０１３０３０４８２０６ 基于收缩二叉树结构网格剖分的大地电磁 二维有限单元法正演研究 张志勇① 刘庆成② 　 （ 中南大学地球科学与信息物理学院，湖南长沙 ４１００８３； 东华理工大学核资源与环境教育部重点实验室，江西南昌 ３３００１３） ① ② 张志勇，刘庆成 基于收缩二叉树结构网格剖分的大地电磁二维有限单元法正演研究 石油地球物理勘探， ． ． ， （）： ２０１３４８３ ４８２ ４８７ ～ 摘要 在大地电磁的正演计算中，通常采用不断增加垂向步长以减少垂向剖分节点数，以提高计算效率。本文 　 针对规则网格剖分无法实现随深度增加横向变步长剖分的实际情况，提出了一种基于二叉树结构的收缩网格 剖分算法。该算法利用二叉树结构在实现四边形剖分的基础上，通过四边形三角形化实现三角剖分；利用对边 的访问实现快速的双二次插值；利用带二次插值节点的三角形再剖分实现网格细化，同时引入缓冲层解决二叉 树收缩过快问题。试算结果表明，采用这种收缩结构的二维大地电磁正演精度与规则剖分下的计算精度相当， 同时大大减少了节点数，提高了计算效率。 关键词 二叉树 有限单元法 大地电磁 数值模拟 　 　 　 　 中图分类号： 文献标识码： Ｐ６３１ Ａ 　　 ［］ ３ 等 采用了对四边形单元进行二次剖分为三角形的 ［］ ［］ ４ ５ １ 引言 方案；为提高计算精度，史明娟等 和徐世浙 采 　 ［］ ６ 用了二次插值四边形网格；刘云等 利用等参单元 ［］ 自Ｃｏ ｏｎ１ 将有限单元法引入地球物理邻域 变换，采用四边形网格实现了带地形的正演问题； ｇｇ ［］