基于异常复电位2.5维CR有限元数值模拟.pdf

基于异常复电位2.5维CR有限元数值模拟 李勇1林品荣1李桐林2王有学3尚晓1 1.中国地质科学院地球物理地球化学勘查研究所，河北廊坊065000 2.吉林大学地球探测科学与技术学院，长春 130026 3.桂林理工大学资源与环境工程系，广西桂林541004 摘要：论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交 流双点电源复电阻法的变分问题；然后利用有限单元法求解变分问题，在矩形网格中增加2条对角线形成 三角网格对研究区域进行剖分，在单元中对复电位和岩、矿石复电导率进行线性插值，将变分方程化为线 性代数方程组；最后，求解线性方程组，得到各节点的异常复电位值，并根据选定的测量装置计算出地表的 视复电阻率，进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法 的有效性后构建模型算例，给出了其不同测量装置下的计算结果。 有限单元法；复电阻率法；Cole-Cole模型；异常复电位 P631.3 A 1671-5888(2011)05-1596-09 FiniteElementMethodforSolvingAnomalous ComplexPotentialof2.5-DComplexResistivity LIYongLINPin-rongLITong-linWANG You-xueSHANG Xiao 2010-11-22 国家“863”计划重点项目(2007AA060601)；中国地质调查局项目(1212011120193)；中国地质科学院地球物 理地球化学勘查研究所基本科研业务费专项资金(AS2010J01，AS2007J14) 李勇(1980-)，男，云南宾川人，硕士，主要从事地球物理数值模拟与反演成像研究，E-mail:liyong@igge. cn。 幅频谱(a) Fig.3P (12) (24) mericalsolution delshownonFig.5 @@[1]罗延钟，熊宗厚，崔先文.同时存在激电和电磁效应的 频谱激电异常性质[J].物探与化探，1987，11(2):106 -114 LUOYan-zhong,XIONGZong-hou,CUIXianwen. PropertiesofthespectralIPanomalyincoexistenceof IPandelecteomagneticeffects[J]. Geophysical GeochemicalExploration,1987,11(2):106-114. @@[2]刘崧谱激电法[M].武汉：中国地质大学出版社， 1998. LIUSong.Spectralinducedpolarization[M].Wuhan: ChinaUniversityofGeosciencesPress,1998. GeologyExploration,2006,34(1):70-73. @@[3]席振铢，张友山，张宪润.运用对称四极测深研究层状 @@[7]张辉，李桐林，董瑞霞.体积分方程法模拟电偶源三维 极化介质频谱特征[J].中南工业大学学报：自然科学 电磁响应[J].地球物理学进展，2006，21(2):386- 版，2003，34(1):8-10. 390. XIZhen-zhu,ZHANGYou-shan,ZHANGXian-run. ZHANGHui,LITong-lin,DONGRui-xia.Modeling SpectralcharacteristicsofmultifrequencyIPsounding 3-Delectromagneticresponsesoftheelectricdipole with schlumberger array over a multi-layered usingvolumeintegraleq