不同装置二维电阻率CT数值模拟研究.docVIP

不同装置二维电阻率CT数值模拟研究 　　【摘 要】本文通过对比分析相对应的视电阻率断面图，总结出了三极装置的数值模拟结果与所建地质模型最符合，分辨率和准确度较其它装置高，从而对以后探测不同地质异常体应采用何种装置提供理论依据。 　　【关键词】电阻率CT；RES2DMOD正演；RES2DINV反演；分辨率 　　0 引言 　　本研究利用Loke的二维正反演软件包，对温纳（Wenner）、二极（pole-pole）、偶极（dipole-dipole）、三极（pole-dipole）装置在纵向水平产状模型对比分析，总结出了不同装置电极排列的优缺点和适用性。正演模拟采用的数值方法为二维有限元法，正演模拟的响应视电阻率用做反演的输入数据，建立的正演模型均具有电阻率突变界面的性质，采用的反演方法为圆滑约束最小二乘法。 　　1 二维电阻率层析成像的正反演计算 　　1.1 二维电阻率法的理论基础 　　电阻率层析成像的理论基础是欧姆定律和电流连续性方程。欧姆定律： 　　J=■E=-■■U（1） 　　电流连续性方程：divJ=0（2） 　　求解位场分布不同边界条件下的拉普拉斯方程： 　　div gradU=■2U=0（3） 　　式中：J为电流密度（A/m2）；ρ为电阻率（Ω?m2）电场强度（v/m）；U为电位（V）。 　　1.2 电阻率层析成像的有限元算法 　　有限单