不同装置二维电阻率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 电阻率层析成像的有限元算法 　　有限单元法是一种以变分原理和剖分插值为基础的数值计算方法。首先要利用变分原理将给定边值条件下求解电位U的微分方程问题；其次，离散化连续的求解区；第三，在各单元上近似地将变分方程离散化，导出以各节点电位值为变量的高阶线性方程组；最后解此方程组算出各节点的电位值，得到空间场的分布[4-5]。 　　1.3 层析成像的圆滑约束最小二乘法 　　RES2DINV 高密度电阻率二维反演程序是基于圆滑约束最小二乘法的反演计算程序，采用了准牛顿最优化非线性最小二乘法（Loke and Barker 1996），使得大量数据下的计算速度较常规最小二乘法快10倍以上且占用内存较少。圆滑约束最小二乘法基于以下方程： 　　J′J+uF=J′g（4） 　　其中F=fxfxT+fzfzT；fx=水平平滑滤波系数矩阵；fz=垂直平滑滤波系数矩阵；J=偏导数矩阵；J′=J的转置矩阵；u=阻尼系数；d=模型参数修改矢量；g=残差矢量。 　　2 电阻率CT数值模拟 　　本文利用RES2D正反演软件，研究了不同装置类型在纵向水平产状组合模型的异常响应特征，分析了不同装置在两种理论模型上的分辨率。 　　2.1 单一低阻、高阻模型 　　模型目标区长100m，深度为16m，视电阻率的背景值为100Ω?m。低阻体和高阻体的厚度、宽度均相同，为1m、11m。模型参数设置如表1，模型如图1。 　　表1 模型1参数 　　图1 单一低阻、高阻模型 　　在温纳排列、二极排列、偶极排列和三极排列装置中进行正演，得到的视电阻率拟断面图，如图2（按照从上到下顺序依次为：温纳、正向三极、二极、反向三极、偶极）。 　　图2 单一模型正演拟断面图 　　从上述几种装置正演断面图可得出以下结论： 　　（1）各模型均能反应左侧为低阻特征，右侧为高阻特征。 　　（2）二极和偶极装置断面图低阻异常均表现为类抛物线特征，区别在于对高阻体的响应上，二极高阻异常为半圆形，偶极高阻异常仍为类似于低阻的抛物线形。 　　（3）正、反向三极装置对低高阻异常的响应呈现出镜像对称性，形状亦类似于抛物线。 　　（4）温纳排列显现为分别以低、高阻体为中心的发散状。 　　据此，要想准确定位断面中存在的两个高阻体和两个低阻体是十分困难的。所以，通过进行反演模拟看能否达到比较好的效果。 　　通过分析可知：温纳和二极装置只能反映上部异常体，对于下部异常体不能做出准确的响应；偶极和三极装置均能够反映出上、下异常体，但偶极装置的分辨率较三极装置低。 　　2.2 低、高阻混合模型 　　各个参数的设置和单一模型相似，只是改变了下部异常体的类型。通过对比分析得出以下结论： 　　（1）温纳和二极装置的响应同单一模型在形状上没有什么差别，只是下部视电阻率值的差异较单一模型增大了。 　　（2）偶极装置所示断面图能够反映出上、下不同地质异常体的存在，同样呈现抛物线类型。 　　（3）正、反三极装置还是镜像对称，但能够较偶极装置更好的反映出上、下不同地质体的交互。 　　存在的共同问题是均不能准确定位交互异常地质体的位置。同样进行反演模拟，看出，二极装置只能较为准确的反映出上部异常体，下部异常体的响应和实际模型偏离很大，存在严重的错位；温纳装置对交互地质异体有显著地响应，优于二极装置，但对异常体4下部低阻体的反映存在错位现象；偶极装置和正、反