第1章电法勘探数字解释基础08.doc

PAGE  1-PAGE 25 0 绪论 地球物理资料解释的基本问题 两个基本问题： 获得可靠的反映地质情况的地球物理数据及有关物性资料； 对这些地球物理资料进行正确的解释。 定性解释：依据地球物理异常的特征定性判断引起异常的地质原因，大致判断异常源的形状、大小、产状和埋深等。 定量解释：对地球物理异常定量计算异常源的几何参数和物理参数 地球物理解释：依据地球物理异常，应用数学物理方法推断场源体物性参数、埋深、位置、形状、大小、产状等要素。 地质解释：侧重于结合地质和其他资料对地球物理异常及其场源体做出地质解释。 正演问题：已知地质体的形状、埋深及其与围岩的物性参数，求取该种地球物理场的剖面曲线或分布。 反演问题：已知地球物理曲线，求取地质体的形状、埋深、和岩石物性。 电法勘探数据处理的发展历史 70年代，D.P.Ghosh将线性滤波方法引入电测深领域。 国内80年代，电测深一维正反演开发、应用。 三维正演亦始于70年代：J.H.Coggon有限单元法 A.Dty等 积分方程法 边界单元法 有限差分法 常用方法 一维正反演:拟和核函数法;直接拟和视电阻率法;构组Dar-Zarrouk曲线法;其它快速方法。 二维正演：有限单元法、边界单元法、有限差分法、积分方程法。 高密度电阻率法的二维反演。 交流电法的数值计算问题 发展 三维反演 电阻率CT成像 电测深一维数值计算基本理论 §1-1 地下稳定电流场的基本问题 一、稳定电流场基本规律 欧姆定律微分形式: 稳定电流场的连续性: 稳定电流场的势场性: 无旋场： 二、稳定电流场基本方程和边界条件 拉普拉斯方程 边界条件: 第一类边界条件: 第二类边界条件: 第三类边界条件:界面两侧有限: 三、均匀各向同性半空间点源电场 采用各种装置形式的电法勘探，在地下建立起的稳定电流场多为双异极性点源场，因此研究点源场在均匀大地的分布有一定的意义。 1．地表半空间中沟点源电流场 在地面设置一电流强度为I的电流源A(I)，求电阻率为 的均匀各向同性半空间中距点源A为r的M点处电位值U(图1—2)。这种简单地电条件下的电位分布计其问题，可在 给定边界条件下直接求解拉普拉斯方程式，得到均匀半空间条件下电位计算表达式 由此可知，地下任意点处点源场电位U、电流密度j和电场强度E均与电流强度I成正比．而电位U与r成反比，E及j与r的平方成反比。 在地面电法勘探中，最常用的装置是采用双异极性点源场(图1—3M)。即在均匀地下半 空间设置相距2L的异极性电极A(+I)和B(-J)，根据电场叠加原理，可写出地下任意点 见的电位为 （推广到视电阻率公式计算过程） 地下任意一点电流密度的矢量值为 2．偶极电流源 当两个异极性电流源A（+I）和B（-I）的距离与观测点到AB中心的距离相比很小时，可把AB看作为电偶极子，称为偶极源（如下图）。 偶极源置于地面时，地下任意点电位可写成： 式中——偶极距， ——AB间距 空间任一点的电场强度可分解为方向的和垂直方向的 其总场为 四、水平层状大地表面点电场及视电阻率表达式： 水平层状大地电位方程 地电模型：n层参数： ，，………，， h1,h2………hn-1。 各层内电位函数为U1，U2，……… Un 采用柱坐标（电位Z轴对称，U与无关） 拉氏方程： ① 用分离变量法求解： 设 U(r,z)=R(r)Z(z) ② 代入①得： 两边相等，只有都等于一个常数才有可能。 ∴ ③ ④ ③式的解为: ——零阶贝塞尔函数 ④式的解为: , ∴①式的通解为: ⑤ 式中: ,为待定的积分变量的函数 水平层状大地各层电位的表达式 利用边界条件: 时 第一层电位函数: 界面上积累电荷引起的影响电位函数 ∴⑤改为: 根椐地面上任意点 因此 ∴ 韦伯--李普希茨积分 第一层电位为: ⑥ 第二～n-1层: ⑦ 第