电法资料处理与解释电法资料处理与解释 视极化率曲线正演程序设计视极化率曲线正演程序设计姓姓 名：名： 丁钰峰丁钰峰 学学 号：号： 2012100389920121003899 专专 业：业： 地球物理学地球物理学 班班 学学 号：号： 061121061121 院院（（系系）） ：： 地空学院地空学院 指导教师：指导教师： 师学明师学明 2015 年 6 月目录目录一、基本原理.2二、程序设计.3三、程序代码.41、前期准备 .42、输入数据 .5（1）公式参数 .5（2）层参数设定 .53、求等效电阻率 .64、计算视极化率，调用两个子程序 .65、视极化率曲线成图 .6附：子程序 7（1）wjhxysdzl（无极化效应视电阻率） 7（2）jhxysdzl(极化效应视电阻率) .8四、相关成果.9五、心得体会.9一、基本原理一、基本原理激电二次场的形成和衰减相对于电磁感应现象来说是一个比较缓慢的过程 在时间域，充放电过程大约为数秒到数百秒；在频率域，激电效应主要发生在 0.0110Hz 的低频段。
对于这样缓慢变化的激电场，可以近似按照稳定电流𝐻𝑧~场方法进行计算，其电位分布仍满足拉普拉斯方程：（3-10）∇2𝑈 = 0因此，计算激发极化场的基本方法就是解拉普拉斯方程A．等效电阻率 对于电阻率为 的均匀介质，当不存在激电效应时，在地面上国内采用任意𝜌装置进行观测，并按下式计算电阻率：（3-11）𝜌 = 𝐾∆𝑈1𝐼若介质存在激电效应，此时按上式计算的电阻率为：（3-12）𝜌*= 𝐾∆𝑈𝐼式中，为总场电位差∆𝑈由于，故可见介质的激发极化效应等效于介质电阻率∆𝑈 ∆𝑈1𝜌* 𝜌的增大，故为等效电阻率，在长时间供电的情况下，极化二次场达到饱和，𝜌*此时η =∆𝑈2∆𝑈=∆𝑈 - ∆𝑈1∆𝑈故（3-13）η =𝜌*- 𝜌𝜌*B．体极化的边界条件 在体极化条件下，极化单元分布在极化体内故在极化体与围岩的界面上 没有点位跃变，即极化总场点位连续3-14）𝑈(1)= 𝑈(2)由于长时间充电达到饱和时的总场可视为稳定电流场，故在界面上总场电流密 度的法向分量连续，即（3-15）𝐽(1)𝑛= 𝐽(2)𝑛或（3-16）１𝜌∗ １∂𝑈(1)∂ｎ＝１𝜌∗ ２∂𝑈(2)∂ｎ式（3-16）采用了与总场电流密度相对应的等效电阻率及。
𝜌∗ １𝜌∗ ２C．体极化场的计算将式（3-15） 、式（3-16）与无激电效应的一次场边界条件对比，可以发现 它们在形式上完全相同而且一次场电位与总场电位都同样满足拉普拉斯方程， 故它们的解在形式上也完全相同由此得出结论：将一次电位表达式中各介质的电阻率换成相应的等效电阻率，便可得到极化总场的电位表达式从一𝜌𝑖𝜌* 𝑖次场和总场的电位表达式可以很容易得到激电二次场以及视极化率的表达式二、程序设计二、程序设计根据实验原理的内容可以制作出如上图所示的流程图，运用两次循环计算出相关参数即可实现正演目的，这也是本次编程的理论公式的图像化成果具体化到程序，相关的流程及参数设置如下图所示：三、程序代码三、程序代码1、前期准备2、输入数据（1）公式参数（2）层参数设定3、求等效电阻率4、计算视极化率，调用两个子程序5、视极化率曲线成图附：子程序（1）wjhxysdzl（无极化效应视电阻率）（2）jhxysdzl(极化效应视电阻率)四、相关成果四、相关成果根据反演出图可以看出，从第一层的数据和最后一层的数据可以看出来， 成图的结果和原先的设定是吻合的，说明本次程序编辑是成功的五、心得体会五、心得体会本次程序设计原理与一维直流电测深几乎相同，所以在原理推导及程序流 程设计上面不需要花费太多的精力，而在设计程序上面因为有了前期的设计程 序经验以及相关的知识积累，所以显得不是特别地困难。
本次实验利用相关的视电阻率推导结论，用无极化效应视电阻率数据直接 替代极化效应视电阻率以得到视极化率，便可以方便快捷地进行视电阻率的计 算我们在接下来的学习中，也可以主动去探索相关的简便算法以及数学结论， 这样将会对我们的学习大有益处。