复电阻率法三维正反演问题研究

【摘要】：电磁法勘探是一种有效的地球物理手段,随着高分辨率探测技术的发展,对电磁法正反演技术要求也越来越高。复电阻率法可以用多个参数从不同角度评价地质异常体,适合做精细反演解释。目前,复电阻率法正反演方法还不能满足实际工作的要求,三维数值模拟及反演解释仍处于探索阶段。本文研究复电阻率法三维正反演算法,推动复电阻率法正反演进入实用化阶段。 复电阻率法三维正演问题采用交错采样有限差分算法。选取有限长度电偶源,将总场分解成一次场和二次场分别求解。在有限长度电偶源激发下,通过解层状介质的边值问题,采用快速汉克尔变换计算得到一次场的分布。基于Cole-Cole模型建立复电阻率模型,采用交错网格对研究区域剖分,将连续的积分方程离散化,设定边界处二次场值为零,形成线性方程组,解线性方程组可得到二次场。通过与已有算法计算结果对比来检验了三维交错采样有限差分算法的精度。 复电阻率法三维反演的目的是获得Cole-Cole模型参数的分布,本文采用了三维共轭梯度—最小二乘顺序反演方案。为了获得每个网格单元随频率变化的复电阻率值,三维共轭梯度反演应用“拟正演”算法,每进行一次模型参数更新量计算只需要1次正演和2次“拟正演”,避免实际计算雅克比矩阵。计算雅可比矩阵的转置与一个向量的乘积可得到模型的更新方向,雅可比矩阵与另一个向量的乘积可得到模型更新步长。计算中仅需存储上一次反演模型更新方向,有效减小了存储空间的需求。 对三维共轭梯度反演获得的复电阻率模型,利用最小二乘法反演获得每个网格单元Cole-Cole模型参数的分布。使用理论模型合成数据对反演算法的稳定性和有效性进行了检验。