基于八叉树网格的直流电阻率法有限元数值模拟  

Direct Current Resistivity Finite Element Simulation Based on Octree Mesh

作　　者：邹桂红[1] 梁华庆[1] 耿敏[1] 

机构地区：[1]中国石油大学(北京)地球物理与信息工程学院,北京102249

出　　处：《科学技术与工程》2013年第34期10116-10121,10128,共7页Science Technology and Engineering

摘　　要：在VC环境下,使用C语言,实现了基于八叉树网格的三维直流电阻率法有限元数值模拟程序。三个典型模型的算例用来证明八叉树网格方案适合于复杂地电模型的数值模拟。数值模拟了球体模型。随着网格对球面边界的不断逼近,数值模拟的精度也逐步提高。计算了覆盖层下倾斜板状高阻体模型,在偶极-偶极视电阻率断面上,能清晰地看见高阻体模型的形态、埋深。仿真了联合剖面法测量三维地形下埋藏低阻椭球体的情况。由于地形的影响,在视电阻率曲线上不能发现低阻异常的存在,利用比较法进行地形改正后,低阻异常得到凸显。Geoelectric numerical modeling program based on octree mesh is developed using C language under VC development environment.Three experiments are used to prove the scheme based on octree mesh benefits complex geoelectric model simulations.The first experiment about the sphere model shows that less geometric error,higher simulation accuracy.The second experiment is about the calculation of the model with inclined high-resistivity plate-like body buried,and the contour and the burial-depth of the high-resistivity plate-like body are show clearly on dipole-dipole apparent resistivity pseudosection.The third experiment is to simulate the measurement of a lowresistivity ore-body buried under three-dimensional terrain using composite profiling method.Because of topography influence,low resistivity anomaly can not be found from apparent resistivity curves.After eliminating the effect of the terrain by comparative method,low resistivity anomaly is highlighted.

关 键 词：八叉树网格 有限元法 直流电阻率法 

分 类 号：P631.322[天文地球—地质矿产勘探] O245[天文地球—地质学]