点源场井地电位测量三维有限元模拟.PDF

第 卷 第 期 地 球 物 理 学 报 ， ５３ ３ Ｖｏｌ．５３ Ｎｏ．３ 　 　 　 　 　 年 月 ， ２０１０ ３ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｍａｒ． ２０１０ 　 　 　 李长伟，阮百尧，吕玉增等 点源场井 地电位测量三维有限元模拟 地球物理学报， ， （）： ， ： ／ ．  ． ２０１０５３３ ７２９ ７３６ＤＯＩ１０．３９６９ ． ～ ｊ ｉｓｓｎ．０００１５７３３．２０１０．０３．０２８ ， ， ， ＬｉＣＷ ＲｕａｎＢＹ ＬüＹＺ ｅｔａｌ．ＴｈｒｅｅｄｉｍｅｎｓｉｏｎａｌＦＥＭ ｍｏｄｅｌｉｎ ｏｆ ｏｉｎｔｓｏｕｒｃｅｂｏｒｅｈｏｌｅｒｏｕｎｄｅｌｅｃｔｒｉｃａｌ ｏｔｅｎｔｉａｌ ｇ ｐ ｇ ｐ （ ）， ， （）： ， ： ／ ｍｅａｓｕｒｅｍｅｎｔｓ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２０１０ ５３３ ７２９ ７３６ＤＯＩ１０．３９６９ ．ｉｓｓｎ．０００１５７３３．２０１０．０３．０２８ 狆狔 ～ ｊ 点源场井 地电位测量三维有限元模拟  ， １２ ２ ２ ２ ２ ２ 李长伟 ，阮百尧 ，吕玉增 ，段长生 ，杨庭伟 ，王建历 １中南大学信息物理工程学院，长沙 ４１００８３ 　 ２桂林理工大学地球科学学院，桂林 ５４１００４ 　 摘 要 井地电法是深部地质矿产勘查的重要方法之一，研究高效率、高精度、适用于起伏地形等复杂条件下的正 　 　 演模拟算法具有理论和实际意义 本文研究了点源场井地电位测量的三维有限元正演模拟，采用与井轴一致倾斜 ． 的计算区域，使剖分单元随地形起伏，用预处理共轭梯度法求解线性方程组，利用行压缩存储和改进的行压缩存储 方式来存储刚度矩阵和预处理矩阵，实现了起伏地形和倾斜井情形下的正演计算 此外，本文利用仿射坐标变换技 ． 术，给出了平行五面体单元精确的单元积分公式，这种处理技术同样可推广应用于其他平行多面体单元的单元分 析，相比于采用等参变换和高斯数值积分，能极大地提高计算效率 文中构建了一些模型算例，其计算结果验证了 ． 程序的可靠性和高效性． 关键词 井地电法，有限元，三维数值模拟，起伏地形，倾斜井，仿射变换