基于单元细分H-自适应有限元全张量重力梯度正演.pdfVIP

第 卷 第 期 地 球 物 理 学 进 展 ， ２５ ３ Ｖｏｌ．２５ Ｎｏ．３ 　 　 　 　 　 　 年 月（页码： ） ２０１０ ６ １０１５ １０２３ ～ ＰＲＯＧＲＥＳＳ ＩＮ ＧＥＯＰＨＹＳＩＣＳ Ｊｕｎｅ２０１０ 　 　 曹书锦，朱自强，鲁光银，等 基于单元细分 自适应有限元全张量重力梯度正演 地球物理学进展， ， （）： ， ． Ｈ ． ２０１０２５３ １０１５ １０２３ ～ ： ／ ＤＯＩ１０．３９６９ ．ｉｓｓｎ．１００４２９０３．２０１０．０３．０４１． ｊ ， ， ， ＣａｏＳＪ ＺｈｕＺＱ ＬｕＧＹ 犲狋犪犾．Ｆｏｒｗａｒｄｍｏｄｅｌｌｉｎ ｏｆｆｕｌｌｒａｖｉｔ ｒａｄｉｅｎｔｔｅｎｓｏｒｓｂａｓｅｄＨＡｄａｔｉｖｅｍｅｓｈｒｅｆｉｎｅｍｅｎｔ． ｇ ｇ ｙｇ ｐ （ ）， ， （）： ， ： ／ 犘狉狅狉犲狊狊犻狀犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２０１０ ２５３ １０１５ １０２３ＤＯＩ１０．３９６９ ．ｉｓｓｎ．１００４２９０３．２０１０．０３．０４１． 犵 狆狔 ～ ｊ 基于单元细分 自适应有限元全张量重力梯度正演 犎 曹书锦， 朱自强， 鲁光银， 曾思红， 郭文斌 　 　 　 　 （中南大学信息物理工程学院，长沙 ４１００８３） 摘 要 在重力梯度各阶张量正演中，难于求取复杂地质体的解析式，转而利用已知解析解的简单规则几何形体近 　 　 似剖分地质体，其会引起较大的几何拟合误差 本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元 ． 分析给出了满足地下复杂地质体正演的型函数，引入体积分的误差指示器用以估计全张量重力梯度误差，提出了适 用于全张量重力梯度正演的自适应计算的策略及迭代算法，数值算例证明了本文方法的正确性和有效性． 关键词 全张量重力梯度正演，非结构化网格，自适应有限元，复杂形体 　 ： ／ 中图分类号 ， 文献标识码 ＤＯＩ１０．３９６９ ．ｉｓｓｎ．１００４２９０３．２０１０．０３．０４１ Ｐ６３１Ｐ３１２ Ａ ｊ 　　　　 　 　　　　　 　 　　　　 　 犉狅狉狑犪狉犱犿狅犱犲犾犾犻狀 狅犳犳狌犾犾狉犪狏犻狋 狉犪犱犻犲狀狋狋犲狀狊狅狉狊犫犪狊犲犱 犵 犵 狔犵 犎犃犱犪狋犻狏犲犿犲狊犺狉犲犳犻狀犲犿犲狀狋 狆 ，