井地电阻率法及双频激电三维数值模拟与反演研究

【摘要】： 本文主要研究点状电流源(单点或多点)及任意线电流源供电时,在井中供电—地面接收方式下的三维电阻率正、反演问题。通过引入Cole-Cole复电阻率模型,解决了相应的频率域井地双频激电的正、反演问题。 在三维电阻率正演方面,通过对地下点电流源的线积分,计算出任意线电源供电时,在均匀介质情况下的空间任意点的理论电位。三维电阻率的有限差分法正演计算使用ICCG算法对正演方程组求解,可以大大加快正演计算的速度。任意线电源供电时,三维离散化网格在边界节点采用异常电位为零的强制边界条件,在不均匀网格剖分(边界网格长度远大于中间网格长度)的情况下,正演计算仍能取得非常高的精度,简化了正演计算。视电阻率幅值响应与电阻率异常体的空间关系的数值模拟试验,为野外观测系统的设计提供了非常有意义的指导。 在三维电阻率反演方面,针对不同的实际地电模型的反演需求,采用了多种不同的计算Jacobi矩阵元素的方法。从Jacobi矩阵元素的物理意义上认识到了在垂直方向上,较深部电阻率异常体被反演计算过程重构成浅部的异常块体的原因,即反演的异常体重心上飘现象。研究证明:三维电阻率反演的水平分辨率较高,垂向分辨率在很大程度上取决于最佳Lagrange乘子λ。 以Cole-Cole模型来模拟地下电性单元的复电阻率,双频激电的数值模拟结果与物理模拟结果非常相似,证明了Cole-Cole模型的正确性,也证明了0.nHz～nHz频率的双频电流对极化体反映的灵敏性。适用于井地观测方式的二极法观测系统的数值模拟试验显示了井地双频激电对电性异常体、极化异常体的反映十分灵敏,其视谱剖面的梯度带能非常好的勾绘出异常体的水平投影。不同观测方式的零频电阻率的反演结果说明了对于双频激电,用低频视振幅(或视实分量)剖面曲线或图像来反演零频电阻率的可行性;极化率的反演结果说明用高频视相位(或视幅频率、视虚分量)剖面曲线或图像来反演极化率的合理性;在固定反演区地下网格单元的时间常数和频率相关系数的情况下,反演出的极化率并不代表异常体真实的极化率,应该是极化体激发极化效应的综合反映,其空间参数的意义大于其幅值的意义。 最后,通过数值模拟,对井地观测方式的伪随机多频波供电勘探方法进行了可行性分析。