多源条件下直流电阻率法有限元三维数值模拟中一种-地球物理学报.PDF

第 卷 第 期 地 球 物 理 学 报 ， ５９ ９ Ｖｏｌ．５９ Ｎｏ．９ 　 　 　 　 　 年 月 ， ２０１６ ９ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｓｅ．２０１６ 　 　 　 ｐ 张钱江，戴世坤，陈龙伟等 多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件 地球物理学报， ．２０１６． ． ５９ （）： ， ： ／ ９ ３４４８３４５８ｄｏｉ１０．６０３８ ｊｇ ， ， ， Ｚｈａｎ ＱＪ ＤａｉＳＫ ＣｈｅｎＬＷ ｅｔａｌ．２０１６．Ａｎａ ｒｏｘｉｍａｔｅｂｏｕｎｄａｒ ｃｏｎｄｉｔｉｏｎｆｏｒＦＥＭｂａｓｅｄ３Ｄｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ ｇ ｐｐ ｙ （ ）， （）： ， ： ／ ｗｉｔｈｍｕｌｔｉｓｏｕｒｃｅｄｉｒｅｃｔｃｕｒｒｅｎｔｒｅｓｉｓｔｉｖｉｔ ｍｅｔｈｏｄ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ５９ ９ ３４４８３４５８ｄｏｉ１０．６０３８ ｙ 狆狔 ｊｇ 多源条件下直流电阻率法有限元三维数值 模拟中一种近似边界条件 ，， ， ， ， ， １２３ １２ ３ １２ １２ １２ 张钱江 ，戴世坤 ，陈龙伟 ，强建科 ，李昆 ，赵东东 １中南大学有色金属成矿预测与地质环境监测教育部重点实验室，长沙　４１００８３ ２中南大学地球科学与信息物理学院，长沙 ４１００８３ 　 ３桂林理工大学地球科学学院，桂林 ５４１００４ 　 摘要 在多源直流电阻率法有限元三维数值模拟中，传统混合边界条件由于刚度矩阵与源点位置相关，无法实现 　 线性方程组的快速回代求解，目前常用齐次边界条件或无限元边界进行替代，虽然实现了快速回代求解，但同时也 降低了数值模拟的精度 为了实现快速回代求解，并确保数值模拟的计算精度，本文提出了一种近似边界条件方 ． 法 其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来，使得分离后的刚度矩阵与场源位置 ． 无关 并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中，当场源位置发生改变时，只需 ． 要重新计算右端源项就可以实现快速回代求解 理论模型数值计算表明，在水平地形条件下，本文边界条件数值精 ． 度优于混合边界条件；在起伏地