一种基于ansys的瞬变电磁正演研究.docxVIP

一种基于ansys的瞬变电磁正演研究 摘要：瞬变电磁正演是研究瞬变电磁响应规律的有效途径，瞬变电磁的高精度、快速三维正演技术，是当前瞬变电磁正演研究的热点。ANSYS软件中的Multiphysics模块也适用于瞬变电磁法的正演模拟问题。本文通过一些模型算例介绍了ANSYS软件在瞬变电磁场三维正演中的应用方法，包括建模、加载、求解设置、后处理的一般过程。数值计算结果表明，利用ANSYS软件进行瞬变电磁正演模拟具有很高的计算精度。合理使用ANSYS软件进行瞬变电磁响应模拟可以提高正演模拟计算的效率，为分析、理解瞬变电磁响应规律提供有力的支撑。 0 引言瞬变电磁法是应用最广泛的电磁探测方法之一，近年来在金属矿勘查、煤田地质勘探、隧道超前预报等领域瞬变电磁的高精度、快速三维正演技术已成为当前瞬变电磁正演研究的热点，国内外很多学者都针对这一问题展开了研究。目前，瞬变电磁的三维正演技术主要可分为频时转换法和直接时域法两大类。频时转换法首先计算三维频率域电磁响应，然后采用余弦变换、G-S变换等频时技术将频率域响应转换到时间域，进而获得瞬变电磁响应。时域法一般通过时间离散来直接求取时域响应，主要包括时域有限差分法、有限体积法、有限单元法等。时间域有限差分法基于交错网格和显式时间差分格式随着数值计算技术的发展，瞬变电磁三维正演的速度和精度均有了很大程度的提高，但仍然缺乏模型建立和剖分简便、交互便利、便于推广使用的成熟瞬变电磁三维正演软件。在目前已有的商业软件中，ANSYS是比较流行的大型有限元分析软件之一1 ANSYS电磁仿真理论基础电磁感应现象可由麦克斯韦方程组进行描述。麦克斯韦方程组实际上由4个电磁学定律组成，即安培环路定律、法拉第电磁感应定律、高斯电通定律和高斯磁通定律，其微分形式为： 式中： 使得耦合方程(1)转化为关于电场和磁场的两个偏微分方程 式中：实际上，求解偏微分方程(3)时，必须要先给定初始条件和边界条件，才能形成一个完整的初边问题。一般地，利用ANSYS进行瞬变电磁问题模拟时，初始条件往往要通过一个静态分析获得；边界条件一般使用Dirichlet边界条件(即直接给定磁矢势和电势在边界上的值)或Neumann边界条件(即通量平行或垂直条件)。2 ANSYS电磁仿真基本流程ANSYS作为一种有限元仿真工具，求解电磁问题的一般思路就是依照前节给出的理论基础，将所处理区域划分为有限个单元，然后将磁矢势和电势所满足的变分问题极小化以获得节点上的磁矢势和电势值，继而导出计算区域的其他待求量。其基本流程为：1) 定义物理环境，主要包括单元类型选取、材料定义、坐标系选用等；2) 进行几何建模，并分配几何体材料特性、单元类型等属性，对求解区域按照一定规则进行剖分；3) 根据问题的特征，施加边界条件和载荷；4) 求解和后处理，利用ANSYS求解器求解磁矢势和电势，然后利用其后处理器获得与问题相关的其他量。3 基于ANSYS的瞬变电磁三维正演算例为分析和验证利用ANSYS软件进行瞬变电磁三维正演的正确性和精度，分别设计了均匀半空间模型、层状地层模型和单个三维异常体模型进行了仿真模拟。 均匀半空间模型3.1.1 模型描述如图1所示，仿真模型为回线源激发的均匀半空间，回线半径为100 m, 均匀半空间电阻率值为100 Ω·m。剖分区域分别由线圈区、空气区和地层区三部分组成， 其外边界为一个半径为10 000 m的球面。在3.1.2 几何模型几何模型按照物理模型及各分析对象的特性建立，使用直角坐标系，其3.1.3 有限元模型1)单元类型。ANSYS软件为三维电磁场分析提供了solid97、solid236、solid237等单元类型，其中solid236、solid237为矢量单元类型，基于矢量单元的有限元分析将由ANSYS软件自动采用矢量有限元法进行分析。考虑到基于标量有限单元法求解三维电磁场模拟问题时的“伪解”现象，本次分析使用矢量单元solid236,按照单元选项不同分为2类(表1)。2)材料属性。本次分析地层、线圈、空气3种介质。对应这三种介质，选取3种材料属性，如表2所示。3)网格控制。网格剖分的质量直接影响着计算结果的精度。Solid236单元为一阶单元，其四面体结构属于常应变单元，六面体结构属于梯度单元，当剖分粒度相同时，六面体网格相比于四面体网格具有更高的计算精度4)加载与求解。在线圈中加载20 A/m5) 后处理及分析ANSYS直接求解的是