标量拟解析近似方法模拟立方体的异常电场.pdfVIP

第33卷 第2期 物探化探计算技术 2011年3月 文章编号：1ool一1749(2011)O2__ol47 6 标量拟解析近似方法模拟立方体的异常电场 张金会 ，2 ，孙建国1,2，张亚东4 (1．吉林大学 地球探测科学与技术学院，长春 130026 2．国土资源部 应用地球物理综合解释理论开放实验室，长春 130026； 3．安徽省勘查技术院，蚌埠 233005；4．长庆油田公司勘探部，西安 710021) 摘 要：拟解析近似方法是一种求解积分方程的一种近似方法，它可 以处理强散射或者大扰动 的电磁散射 问题，在计算过程 中避免了传统微分数值方法解决问题时所遇到的大型矩阵或大型 代数方程组的求解。孙建 国[51将其引入直流电场的积分方程 中，并给出了求解异常电场积分方 程的标量拟解析近似公式。在 以前的研究中，已经验证了均匀场中异常球体的拟解析近似解的 精度 ，这里对均匀场 中的立方体异常体进行数值模拟，得到了直流电场 中异常立方体模型的标量 拟解析近似解。由于复杂地 电模型可以用立方体的组合进行模拟，因此对立方体异常电场拟解 析近似解的研究，为三维直流电场 中复杂地 电模型的快速正反演模拟打下了基础。 关键词 ：积分方程 ；拟解析近似 ；直流电场 ；立方体 中图分类号 ：P631．3 文献标识码：A 到下一次迭代中，而积分方程法不存在误差的累积 0 前言 和传递。基于这些优点，前人对积分方程法进行了 很多的研究。TMHabashy等…提出用局部非线 近年来我国资源形势 13益严峻，大规模的经济 性近似方法 ，对 电磁场问题进行研究 ，但这种方法 建设对金属矿的需求越来越大，对于资源的需求应 运算复杂不适于实际应用推广。Zhdanov等 受 该立足于国内，因此要加大 国内金属矿勘探 的力 到TMHabashy的局部非线性近似方法的启发，在 度。在金属矿勘探 中，电法勘探 占有重要地位，研 遇电磁场散射场问题时提出了拟线性近似方法求 究直流电场快速数值模拟，对于推进 电法勘探的发 解散射场的积分方程。这种方法通过假设异常电 展有着重要的意义。三维直流电场数值模拟的方 场 E。(r)与正常电场 (r)之间是 由电反射张量 法很多，大致可以分为微分方法和积分方法。在传 ^ A A(r)联系起来 的拟线性关系，即：E。(r)=A · 统的数值方法中，都会遇到大型矩阵或求解大型线 E (r)。简化 了电磁散射的积分方程，从而可 以提 性代数方程组的问题，这对计算机硬件要求高，计 ^ 高计算速度，电反射张量A(r)由最优化的方法求 算周期长。尤其是在复杂地质条件下，这些问题更 为突出，所以寻求快速精确的数值模拟方法 ，一直 出。Zhdanov 在2000年将拟线性近似方法应用 是地球物理工作者的研究重点。积分方程法在解 于电磁场反演，并用来反演美国新墨西哥州的 CSAMT数据。由于拟线性近似时仍无法避开解线 决三维直流电场问题中有其优点，如积分方程法只 需对异常体进行计算，而不需要对背景区域进行计 性代数方程组，因此Zhdanov经过研究后 ，提