电导率连续变化2.5维直流电阻率法有限元数值模拟  

2.5-D numerical modeling on DC resistivity by FEM with continuous variation of conductivity

作　　者：刘云[1] 宋滔[1] 王赟[1] 

机构地区：[1]中国科学院地球化学研究所矿床地球化学国家重点实验室,贵阳550002

出　　处：《地球物理学进展》2014年第3期1194-1200,共7页Progress in Geophysics

基　　金：国家科技973项目"华南大面积低温成矿作用"中"大型-超大型低温矿床成矿规律与找矿预测"课题(2014CB440905);矿床地球化学国家重点实验室"十二五"项目群(SKLODG-ZY125-01)联合资助

摘　　要：在阮百尧、徐世浙(1998,2001)工作基础上,给出四边形内三角网格单元剖分,电导率连续变化2.5维直流电阻率数值模拟方法.在二维起伏地形情况下,有限单元设计为四边形内剖分四个三角网格单元,单元内的电位值和电导率设计为线性变化,将研究区域远边界的第三类边界条件简化为第二类齐次边界条件,重新给出一组适用性较好的Fourier反变换系数,并采用分离式Cholesky分解法快速求解线性方程组.通过对3个地电模型的计算与对比分析,其结果与解析法的均方根误差小于0.2%,计算速度比较常规方法提高了5倍,能有效模拟起伏地形和复杂形态的地电体模型.Based on the study of Ruan and Xu（1998,2001）,the triangle division method in the quadrangle and the thinking of continuous variation of conductivity within each element are presented for numerical modeling on 2.5dimensional（2.5-D）direct-current（DC）resistivity.Under the condition of twodimensional topography,the finite element is designed to the method of four triangles division in a quadrangle,then potential and conductivity are also designed to bilinear variation within each triangle element.In the paper,the third class boundary condition is simplified for the second class of homogeneous boundary condition at far boundaries of the study area,a set of suitability better inverse Fourier transform coefficients is given,and the separate Cholesky decomposition method is used for quickly solving linear equations.By calculating and comparative analysis on three geoelectric models lastly,the result of our method show a high accuracy（the mean square error is less than 0.2%）,the calculation speed of the numerical method can be increased at least 5times than general methods,and the method can model arbitrarily complicated terrain and geoelectric bodies preferably.

关 键 词：电导率连续变化 2 5维 直流电阻率法 有限元 数值模拟 

分 类 号：P631[天文地球—地质矿产勘探]