复电阻率法二维有限元数值模拟.pdf

第50卷第6期 地球物理学报 Vol50．No．6 OFGEOPHYSICS 2007年11月 CHINESEJOURNAL Nov．．2007 蔡军涛，阮百尧，赵国泽等．复电阻率法二维有限元数值模拟．地球物理学报．2007．S0(6)：1869—1876 B G dTwo-dimensionalof finiteelement CatJ T，RuanY，Zhaoz，n complexresisfvityusing modeling Chinese)．2007．钟(6)：1869—1876 复电阻率法二维有限元数值模拟 蔡军涛1，阮百尧2，赵国泽1，朱光萍2 1中国地震局地质研究所。北京100029 2桂林工学院资源与环境丁程系，佳林541004 捕要伴随着复电阻率法的广泛应用，发展精确和快速的正演和反演算法成为复电阻率法研究的重点．本文采 用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究为了提高计算速度，对无穿远边界进行了 近似处理整个正演计算过程分为两步，首先采用有限单元法计算四个不同频率的视复电阻率数据．然后对前一步 得到的视复电阻率数据采用递推算法计算视Cole．Cole参数采用这种正演算法与一维正演的结果进行丁对比，验 证了本文方洼的正确性设计了两个二维极化模型，散值模拟结果表明视复电阻率和CoI}COh视参数等值线断面 图对于异常目标体都有比较明显的反映． 关键词复电阻率法，有限单元法，递推算法，Cole—Cole参数，数值模拟 文章编号0001—5733C2007)06—1869一嘴中图分类号P631 收稿日期2f07一陀一皿，2007—09—19收修定稿 finite Two·dimensionalof elementmethod modelingcomplexresistivityusing CAI Jun-Ta01，RUANBai—Ya02，ZHAOGuo-Zel，ZHU Guang·Pin92 1 lnsttttae 100029，Ch／na qrGea／ogy．c^im岫mk“m自“∞m，脚w AbstractWithwide of forward andinversion applicationcomplexresistivitymethod，searchingprecise，rapid have a hasbecomethefocusof research．We finiteelement algorithms eomplexresistivity developed algorithm in for the of t