【毕业学位论文】基于ANSYS 的激发极化法中地形影响问题的数值模拟-地球探测与信息技术

分类号 密级 U D C 编号 士学位论文 论文题目 基于 激发极化法中 地形影响问题的数值模拟 学科、专业 地球探测与信息技术 研究生姓名 赵 斌 导师姓名及 专业技术职务 何继善 院士中南大学硕士学位论文 摘 要 要 目前， 电法勘探在我国国民经济和社 会发展各个领域应用非常普遍，遇到的问题也多，特别是地形影响问题。我们知道，虽然地形不平对激电异常的影响比对导电性异常的影响小， 但为了正确推断解释激电异常，仍必须考虑地形影响问题 。目前地形改正的方法很多，本文应用限元进行大量的电阻率法和激电法数值模拟， 然后用比值法进行地形改正。众所周知，有限元法适合于模拟形状不规则的模型，计算精度高，是地球物理数值模拟中较常用的方法。 本文首先回顾了电法勘探地形影响问题的研究现状， 总结了现有激发极化法的数值模拟方法， 简单介绍了在进行激发极化法数值模拟过程中的主要步骤及其注意的关键问题。然后介绍了点源场电位边值、变分问题；通过使用 言编制相应的计算程序，以山脊、山谷地形为例，进行了大量的电阻率法和激发极化法地电模型计算， 并对模拟计算出来的结果加以验证和分析，总结了起伏地形情况下的异常分布规律，得到了一些比较满意的结果。 最后总结全文， 表明了关键字： 激发极化法，有限元，地形影响，数值模拟 中南大学硕士学位论文 is in of We P is on be in to P of A of P by in is a it is of by As we to of It is in In we of in up P. we of P we a of of P by by of as of on of is by we At we be in IP 南大学硕士学位论文 目 录 录 摘 要. . 录. 一章 绪论. 1 11 引言.2 电法勘探地形影响问题国内外研究现状.文的主要工作.2章 激发极化法及其数值模拟.1激发极化法原理 .2 激发极化法数值模拟.21有限差分法.22有限单元法.3 .4 .41生成模型的典型步骤.42 .三章 .1 点源三维电场.11 微分方程 .12 边值问题 .13 变分问题 .2 点源二维电场 .21 微分方程.22 边值问题 .23 变分问题 .3 视电阻率和视极化率的计算 .4 应用.四章 激发极化法有限元模拟与分析.1线电流源的.11线电流源.模型计算.点电流源.模型计算.五章 总结全文.南大学硕士学位论文 目 录 . 56 致 谢. .读硕士学位期间主要研究成果.南大学硕士学位论文 第一章 绪论 1第一章 绪论 11 引言 在地球物理的理论计算中 ,存在着两类基本问题8：正问题和反问题。 给定场源的分布，求解场值的大小，这是正问题，或称正演问题。例如，在磁法勘探中，给定了铁矿体的分布形态、磁化率和磁化强度，计算铁矿体产生的磁异常；在电法勘探中，给定了各种岩石电阻率及岩石的分布形态，计算与各种测量排列装置相应的视电阻率；等等。 长期以来，电法正演问题的计算中用场论中的经典方法解析求解，都只针对于少数简单规则的地电条件。除此以外，对大量的实际地电条件尚未得出解析求解的公式。 解决地球物理正问题的途径7有： （ 1）解析法，即用解析公式表达场值的大小，这是各种地球物理教科书讲的最多的方法。但是，只对少数规则形体，例如球体，板状体和水平层等，才能用解析法导出场值的解析表达式，因而，此方法实用的范围非常有限。 （ 2）模型实验法，这是电法勘探中使用较多的方法。例如，各种水槽、土槽模型实验。但模型的制作比较麻烦，尤其是物性分布复杂的模型很难制作。 （ 3）数值模拟法，即根据地球物理中的偏微分方程和边界条件，用数值方法求场值的近似解。这是一种近似方法，但它实用于复杂物性分布和复杂边界形状的地球物理计算，所以适用范围非常广泛，数值模拟法需要进行大量的计算工作。近三十年来，由于计算机技术的迅速发展，克服了计算工作量浩大的困难，数值模拟法已成为地球物理正演的最主要的方法。 在大多数勘探运用中,由于都在较平地形上进行,故可以不考虑地形影响。 然而在地形起伏较大时,如铁路、公路隧道勘探中,它导致异常曲线畸变严重13。有时淹没地质体的信息,有时 歪曲有用信息,有时呈现假 信息,使人们不能做出正确的推断解释,降低了异常解 释精度,甚至导致错误的解 释。因而,研究地形影响及其改正方法,是提高山区电法勘探地质效果的一个关键。本文通过对有地形影响的异常体模型进行电阻率和极化率的数值模拟， 并对所模拟的结果用比值法进行地形改正，从中找到异常体的空间状态与曲线形态的内在规律性，对野外工作及野外数据的定量解释具有一定的指导性。 12 电法勘探地形影响问题国内外研究现状 地形影响问题是20世纪50年代伴随着直流电法在矿产资源勘查等方面的大规模应用开始受到重视。此后随着大地电磁法(可控源音频大地电磁法中南大学硕士学位论文 第一章 绪论 2(代表的无源和有源电磁法应用的扩展，人们也开始关注电磁法中的地形影响问题。理论与实践表明，地形不但可以引起假异常，而且会掩盖地下由矿体或目标物引起的真异常，如果不能正确的认识和消除，就会导致完全错误的解释结果。 早期，日本、前苏联及我国，乃至西方国家的地球物理学者曾采用导电纸模拟、电阻网络模拟、水槽模拟、土槽模拟等物理模拟方法研究地形影响问题，如按一定比例制作导电纸、电阻网络、水槽、土槽等模型来模拟野外情况，得到地形影响的定性认识，并在野外实际测量中加以识别。有些学者也应用模拟结果，对野外观测数据作些简单的计算和改正，试验将地形影响去掉。这些做法基本属定性解释。其后，也有学者采用一些解析计算方法研究地形影响，如保角变换法和基于保角变换的坐标网法、角域地形叠加法、格林理论、模拟电荷法等。但由于数学公式的复杂性和只能计算简单地形，使应用受到一定限制。 近些年来，随着计算机技术的发展和物探解释方法的进步，国内外学者在水平地形的电法、电磁法资料2究开发出适合于起伏地形的2而使地形问题的研究更加方便可行，也使复杂的2别和校正成为可能。起初是用比值法进行校正，后来发展到带地形进行2得了一定效果，一定程度上提高了电法勘探数据的解释水平。 国内外电阻率、激发极化法、充电法等电法地形影响问题的文献比较多，尤在电阻率2物理模拟方法:日本952,1953)借用物理模拟方法研究电阻率曲线的地形影响。 前苏联959,1960,1965)借助物理模拟研究电阻率测量的20世纪6070年代，我国徐世浙、杨文采、易永森、葛为中、何继善等曾采用导电纸2阻网络2槽2槽3一时期基本是定性解释阶段。 解析计算方法:957)从理论上得到棱形和阶梯形不规则地形的均匀电场的畸变。 980)使用985)和徐世浙等(1988)使用格林理论于电阻率地形模拟，徐世浙等(1988)的研究与985)在结果方程上有所不同，徐世浙等使用立体角方法导出位方程。 刘阳等(1996)将模拟电荷法引人纯地形异常的计算，计算几个典型的2 998)针对埋在均匀起伏地形里具有明显电阻率和极化差异的2 第一章 绪论 3界值问题的一种精确解法，研究起伏地形下深度、隆起地形的宽高、电阻率比等对电场振幅的影响，采用“比值法”进行地形校正。所谓比值法即是用实测视电阻率除以“地形校正因子(地形模型视电阻率/地形模型电阻率)”，从而得到改正后或平地形下的视电阻率。 面积分方程法:1984)对充电法和磁电阻率法(发出面电荷源像的新方法处理空气有限个矩形面元模拟地形表面，面元各自倾向和电荷密度不同。电荷值用迭代解技术得到。面电场根据面电荷分布、电流源和像来计算。磁场由涉及面倾向和电场的面积分计算。展示了充电法和边界元方法11: 徐世浙(1984)兰珍等(1986)介绍点源场电阻率法2计算了山脊、山谷、角域地形上联合剖面法视电阻率异常曲线，展示用比值法对实测曲线进行地改的应用效果。徐世浙和汪晓东(1987)介绍用边界元法计算2世浙等(1989)对2少了计算量，用比值法对对称四极剖面曲线进行了地改。 钱家栋和马钦忠(1992)讨论用边界元法解决层状介质中存在2世浙(1993)介绍边界元法模拟点源2兰珍和田宪谟(1997)介绍边界单元法应用于工程地质勘查电阻率法的地形改正实例48。刘继东(1998)介绍求解异常电位的2世浙等(2002)用2宪漠等(1986)介绍点源电阻率法3并计算了3对相应地形下埋藏着良导电性椭球体的用比值法进行了地形改正。徐世浙等(1988)提出模拟3用格林函数将3用边界元方法解积分方程;给出用比值法电阻率3有限差分法:吴小平等(2001)开展起伏地形条件下有限差分电阻率/激发极化3演技术研究25。 有限元方法20,46:971)针对直立导电极化体，用有限元法计算得到平底、斜边的山谷和山脊21973)用有限元技术模拟比较极化率双偶极和三级装置的地形影响。 977)使用有限元方法模拟电和电磁数据。 980)利用2谷、 山脊)对极化率影响的系统研究。中南大学硕士学位论文 第一章 绪论 4对视电阻率地形影响校正使用有限元法计算校正系数，用比值法进行校正。周熙襄等(1986)应用有限元法计算偶极装置下台阶、山脊、山谷2谷2 999)华和李金铭(1999)以山脊、山谷地形为例，通过点源2用比值法对地形影响进行校正23。001)推出有限元2)，可反演带地形电阻率/ 984)给出3该算法进行多种测试并展示实际模型精确结果，采用比值法校正。 本文主要采用有限元法进行大量的激发极化发地形影响问题的数值模拟。 文的主要工作 本文的工作主要由以下部分组成： 一、 绪论：主要介绍了电法勘探地形影响问题的研究现状。 二、 激发极化法及其数值模拟：首 先简要介绍激发极化法的原理和进行数值模拟的方法。然后对后主要阐述总结了在遵循的一些基本原则。 三、 别简要推导了点源三维电场、点源二维电场问题的基本方程、边值问题及其变分问题，然后分析了如何应用四、 激发极化法有限元模拟与分析 ：对点电流源和线电流源激励的电场进行了大量的有限元带地形的地电模型正演模拟。 其中涉及到了有限元计算中的源、边界条件、计算区域的选择、单元剖分等问题，应用 制 言程序,对激发极化法地形影响问题进行大量的数值模拟并加以分析。 五、 总结全文：总结了应用 现激发极化法地形影响问题数值模拟的一些结论和分析结果，说明 中南大学硕士学位论文 第2章 激发极化法及其数值模拟 5第 2 章 激发极化法及其数值模拟 21 激发极化法原理 激发极化法是50年代和60年代初， 在我国开始实验研究和推广的。 实践证明，它是应用最广的和效果最好的一类电法勘探方法，早期是以直流（时间域）激电法为主，通过长期的应用和研究取得了许多重要的成果。为了提高激电法的抗干扰能力和减轻装备，70年代初又开始研究推广了交流（频率域）激电法，这一时期应用最多的主要是变频法。 中南工业大学提出的双频道激电法是对变频法的发展，为了异常的区分和电磁耦合问题，80年代初又开始对频谱激电法进行研究，取得了不少有价值的成果。 激发极化2(缩写为发生在地质介质中，因外电流激发而引起介质内部出现电荷分离，产生一个附加的“过电位” （一种物理化学现象。这个现象可用图2图 2发极化法现象示意图 -1 P 由供电极观察测量电极果地下介质是均匀非极化的，并且保持供电电流不变，那么在如图2不会随时间变化。如果地下存在可极化的介质，比如像图2供入的电流流过矿体时，矿体将产生电化学反应而被“极化”成原电池。此原电池也将产生电流，俗称“二次电流” 。二次电流在空间也要产生电场，在为二次电位差22随时间变化的，常记为 )(2供电时在 所测到的电位差 )(一次电位差1二次电位差 )(2和，俗称总场电位差，它的涵义是供电电场与极化电场二者在 产生的总电位差。由于 )(2时间而变，是时间的函数，故 )(是时间的函数，上述现象可用图2 中南大学硕士学位论文 第2章 激发极化法及其数值模拟 6图 2时间域观察激发极化法示意图 -2 P in 体（或者别的可极化介质）被极化是一个过程。过程随时间变化的速率先快后慢，最后趋于饱和。饱和值的大小 s 是 缩写） ，以及达到饱和经历的时间般如果在图2会看到，虽然经由 矿体的极化依旧存在，它将继续放电。因切断了激发电流，随着放电的进行，矿体的极化也逐渐减小。相应地，断电后供电引起的电位差1不再存在，但极化引起的电位差 )(2不立刻消失，仍然能观测到，并且先快后慢地逐渐减小，经过数秒甚至到数分钟，才最终归于消失。如图2b。上面观察到的现象有点像电 容器的充放电，故习惯上将供电时 (曲线称为充电曲线，而把断电后测到的二次电位差)(2曲线称为放电曲线。 上述电现象因外电流的激发而引起（外电流接通以前决不会发生） ，外电流切断后却并不立即消失，说明电流的激发导致了介质本身的极化，故称之为激发极化。 22 激发极化法数值模拟 模型和模拟实验（统称物理模拟）方法，虽然比数学解析法更易于解决比较复杂地电条件的问题，但仍有其缺点。比如，对比较复杂的地电条件构建模型也十分复杂，此外，实验结果受电位差观测误差的影响，往往精度不高；最后，物理模拟的效率较低，并且实验结果不便于用于自动化数据处理。现代数值模拟方法可以避免这些缺点。 电场的数值模拟，就其原理来说，是早已解决了的问题；但由于其计算工作量大，所以一直未能付诸实现。近年来，由于电子计算技术的飞速发展，为对电场做快速和精确的数值模拟提供了技术条件，使之得以实现。目前，国内外已经有不少单位和学者，成功地将数值模拟技术应用于研究激电异常。他们研究激电异常的途径也是基于“等效电阻率法”原理1，先模拟计算没有激电效应的一次场，然后将地下各地质体的真电阻率换成相应的等效电阻率 )1/(=，并中南大学硕士学位论文 第2章 激发极化法及其数值模拟 7计算二次场和极化率。所以实际上，这种对激电场的模拟方法还是归结为模拟不包括激电效应场的一次场（即稳定电流场） 。目前，用来对稳定电流场作数值模拟的方法主要有两种，即有限差分法和有限单元法8。 221 有限差分法 有限差分法39是从电场（电位）或电磁场所满足的偏微分方程和边界条件出发，将微分方程转变为差分方程，其研究步骤是：首先将研究区域按一定方式离散化，然后在每个单元内设位、场呈线性变化，电性为均匀的，因而微分方程的微分就可用差分来代替，于是就可以建立一组线性差分方程，最后求解此线性方程组即得相应的位场分布。 通常解决稳定电流场空间分布的方法是根据电磁场理论导出一定的微分方程，并用解析法求解；有限差分法则是将导出的微分方程离散化，得到相应的差分方程以求解。具体来说，在稳定电流场中，电流密度j 、电场强度E、电位 之间有如下关系： = ， =（2 ： = （2电流源以外的点上，根据电荷守恒定律有： 0)( = （2者经过矢量运算，将上式改写为： 0)(222=+ 2式便是在非均匀介质中，稳定电流场电位均匀或分区域均匀介质中，它简化为众所周知的拉普拉斯方程： 02= 。 在用有限差分法求微分方程（2数 值解时，首先要在求解区（或模拟区）建立网格，以便将微分方程离散化。对于模拟区（即网格）边缘上的节点（外边界点） ，其电位值可根据具体的地电条件给出，或者表示为该点坐标的已知函数，或者表示为相邻节点电位的线性函数。这样，网格上的每一个节点都有一个表示该点电位与相邻点电位之间关系得线性方程，总合起来，便得到一个未知数（即各点的电位值）和方程个数等于网格总节点数的高阶线性方程组。现代计算技术有多种方法可以用来解这个方程组。 由此， 便可算出网格中各节点的电位值，即求得模拟区内空间各点的电位分布。 中南大学硕士学位论文 第2章 激发极化法及其数值模拟 8222 有限单元法 有限元法46是从位、场所满足的偏微分方程 出发，根据微分方程解与泛函极小问题的等价性，将微分方程和其边界条件转化为相应泛函的变分问题，其研究的步骤是将研究区域按一定方式离散化，设单元内位、场呈线性变化，电性参数均匀，这时泛函是各节点位、场的二次函数，利用求极小的必要条件，即泛函对各节点位、场的变分为零，二次函数的变分为一次函数，由此得到一个线性方程组，解此线性方程组便可得各节点的位、场值。 有限单元法9是近二十多年来随有限差分法之后 迅速发展起来的又一种有效的数值模拟方法。七十年代初， 后其他一些学者进一步完善了它， 使之成为求解复杂条件下电阻率法和激电法问题的有力工具。 有限单元法在作激电法模拟时，也是分别求解地电断面为真电阻率和有效电阻率分布时的稳定电流场（即一次场和总场） ，然后根据等效电阻率法计算激电异常。这与有限差分法的做法是一样的，所以，下面只讨论求解稳定电流场的电位问题。 有限单元法求解稳定电流场电位 V，是根据变分法原理或最小总电能原理将给定边界条件下求解（ V）微分方程的问题，等价地转换成求 V 的泛函的极值问题；再经过离散化，得到由空间各点未知电位值组成的高阶线性方程组；最后在数字电子计算机上求解该方程组，以确定空间各点的电位。 在包括电流源在内的不均匀三度空间中，稳定电流场的电位 V 满足泊松方程： = （ 2中 为介质电导率；V 为稳定电流场电位；电流源体密度。它们都是空间坐标的函数。在直角坐标系中，公式（2写成： =+)()( （2据变分原理，若在区域 的一部分边界1s 上，其余边界2s 上，自然边界条件” ：2， 其中外法线方向。 则求解微分方程(2等价转换为泛函极值问题，即找出一个函数V，使其泛函中南大学硕士学位论文 第2章 激发极化法及其数值模拟 9的全部边界为以知；极小+=+=212222)()()(21限单元法和有限差分法是根据不同的原理作数值模拟，实践证明，两种方法可以得到彼此相近的模拟结果。 有限差分法的程序编制比较简单， 采用这种方法对平面规则形体 （如板状体、层状体等）进行数值模拟比较方便；但它对内节点、内边界点及两种类型的外边界点要分别用不同的电位计算格式。所以，当不同电导率的介质较多或不均匀体表面积为曲面时，使用这种方法不甚有利。 有限单元法易于模拟具有不规则表面的地电模型；此外，这种方法在整个求解区可采用同一种计算格式， 因而可避免采用多种计算格式时给模拟结果带来的误差。但有限单元法程序编制较复杂，占用计算机的内存量也较大1。 本文所用的地电模型主要是是有地形的模型，故采用有限元法进行实现。 是融结构、热、流体、电磁和声学为一体的大型 用有限元软件，可广泛地应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、电子、土木工程、生物医学等一般工业及科学研究。该软件提供了不断改进的功能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网络划分及利用外，时为用户自己进行二次开发提供了友好的开发环境， 23 件介绍 软件创始人是 士，在 士的领导下，件的研究与开发不断汲取当今计算方法和计算机技术的最新发展，领导着有限元界发展的趋势， 并为全球工业界所广泛接受， 拥有全球最大的用户群。 供给用户的操作方式有图形用户界面 、数化设计编程语言 用户界面设计语言后两者也是 化设计、自适应网格以及二次开发的主要基础。一般而言，： 前处理（ 施加荷载与求解（ 后处理（中南大学硕士学位论文 第2章 激发极化法及其数值模拟 10前处理主要在前处理器（完成，包括有限元模型的创建等。前处理的主要任务是建立结构分析的有限单元模型，它是结构分析的开始，一般步骤为5： （1） 分析准备 （2） 设置单元类型 （3） 设定实常数 （4） 定义材料属性 （5） 创建模型 （6） 划分网格 施加荷载与求解主要在求解器（n，/完成，主要操作包括定义自由度、施加荷载、设定分析类型以及求解等。施加荷载与求解的主要任务是进行结构的计算，一般步骤为： （1） 施加荷载条件 （2） 给定边界条件及求解 后处理在后处理器中完成， 通用后处理器（ 时程后处理器（种。前者主要用于静态分析的结果后处理，而后者主要用于动态分析结果的后处理。后处理的主要任务是进行结构的分析。 应用 行电法有限元正演模拟分析，创建有限元模型是整个正演模拟的基础，模型的建立直接关系到模拟的速度和精度，在建立模型时，我们可以进行合理的简化，尽量使进行的网格划分和加载变得简单。在划分网格时，我们可以在电源点和异常体附近划分较密的网格，随着与电源点的距离的增加，可以把网格划的稀疏一点，这样在不增加总节点数的情况下，可以提高我们所关心区域的网格划分的密度，从而提高精度、加快计算速度。下面着重介绍 模与网格划分。 24 限元建模与网络划分 有限元分析是针对特定的模型而进行的，因此，用户必须建立一个有物理原型的准确的数学模型。一般情况下，有限元模型包括所有的节点、单元、材料属性、实常数、边界条件、荷载条件以及用这些条件来定义的物理场特征。建模与网格划分是有限元分析的基础。 利用以快速、方便地建立实际工程系统的有限元模型。 中南大学硕士学位论文 第2章 激发极化法及其数值模拟 11241 生成模型的典型步骤 通常的建模过程41如下： 设计建模方案。首先确定分析目标，设计模型 采取的基本形式，选择合适的单元类型，并考虑如何建立适当的网格密度； 进入前处理器开始建立模型； 建立工作平面； 利用几何元素和布尔运算操作生成基本的几何形状； 激活适当的坐标系； 生成单元属性表（单元类型、实常数、材料属性和单元坐标系） ； 设置属性指针； 通过对实体模型划分网格来生成节点和单元； 模型数据保存； 退出前处理。 242 限元网格划分 现今所有的有限元分析都是采用实体建模，于在里面填充节点和单元，还可以在几何边界上方便地施加载荷，但是几何实体模型并不参与有限元分析，所有施加在有限元边界上的载荷或约束，最终必须传递到有限元模型上（节点和单元）进行求解。 通常网格划分有三个步骤5： 建立选取单元数据 设定网格建立所需的参数 产生网格与修正网格 第一步是建立单元的数据，这些数据包括单元的种类（，单元的几何常数（R） ，单元的材料性质（，及单元形成时所在的 坐标系，以及不同的单元坐标系统。 第二步是设定网格划分的参数，最主要是定义对象边界单元的大小和数目。网格设定所需的参数，将决定网格的大小和形状，这一步将影响分析结果的正确性和过程的经济性。 完成前两步之后即可进行网格划分，即进行第三步的工作，完成有限元模型的建立。如果不满意网格划分的结果，也可以清楚网格划分，重新定义单元的大小和数目，再进行网格划分，直到得到满意的有限元结果为止。 有限元方法的基本思想是将结构离散化， 用有限个容易分析的单元来表示复中南大学硕士学位论文 第2章 激发极化法及其数值模拟 12杂的对象， 单元之间通过有限个节点相互连接， 然后根据变形协调条件综合求解。由此，有限元网格的划分一方面要考虑对各物体几何形状的准确描述，另一方面也要考虑变形梯度的准确描述。网格划分的密度是个重要的问题，太密则会大大增加计算时间，但计算精度却不会成比例地增加。这样，就存在一个最佳网格密度问题，这个问题往往需要多年工作经验的积累。如果前置处理程序能够自动确定网格密度，对节省机时的意义非常大。对于地球物理模型而言，一般都是一些较复杂的地电模型，因此建立适当的地电模型和划分适当网格应遵循一定原则，建立复杂有限元网格划分的基本原则如下30： 1) 网格的大小 网格大小就是所取泛函积分区域的大小，一般来说网格越大越好。对于微分方程边值问题的求解，只有给出正确的边界条件，才能求出域中比较精确的函数值。例如，对于使用第一类边界条件时，要求给定正确的边界函数值。但是对于不均匀的地电断面，边界函数值和其他边界条件值，特别是地下部分，均无法正确求出。因此往往采用地下边界值为零，或由均匀地电条件给值。这时就要求区域的边界远离不均匀区，即要求网格大，否则就会影响计算的精度，但是，另一方面，如果网格内单元的大小不变的话，那么网格太大，势必要大量的增加节点数，从而需要更多的计算机内存和增加计算工作量，这是因为在 点个数的增加，将使线性方程组的阶数增加，于是刚度系数矩阵的形成，解此方程的时间必然会增加。对于网格内部单元的大小，一般来说单元越小，计算精度越高。我们假定 果单元太大，实际函数便可能不满足这个条件，从而增加计算误差。我们还假定单元内电性是均匀的，即电导率为常数，这也要求单元较小，特别是要拟合复杂的地电断面和地形剖面，更需要划分得细致些，才能满足单元内电性均匀的条件。为了克服网格大小和单元大小选择精度和工作量之间的矛盾，我们在划分网格时采用非均匀的网格，网格的中心部分单元小，节点密，边界单元大，节点稀，由中心到边缘单元逐渐放大， 这样既保证了网格有足够的大小， 又保证地电断面的复杂部位位于网格中心，以满足单元内电性均匀和 2) 网格密度 有限元网格数量的多少将直接影响计算结果的精度和计算规模的大小。 一般来讲，网格数量增加，计算精度会有所提高，但同时计算规模也会增加，所以在确定网格数量时应权衡这两个因素综合考虑。 网格较少时增加网格数量可以使计算精度明显提高，而计算时间不会有大的增加。当网格数量增加到一定程度后，再继续增加网格时精度提高甚微，而计算时间划分的计算结果，如果两次计算结果相差较大，应该继续增加网格，重新计算，直到误差在允许范围之内。在决定中南大学硕士学位论文 第2章 激发极化法及其数值模拟 13网格密度时还应考虑分析类型。对于进行电磁场分析时（对本文而言特指直流电场分析） ，网格类型与源点的位置对计算精度结果是有影响的，在二维计算中，网格形状的不对称，源点位置的不对称，都会引起计算结果的不对称；在地下半空间中，源附近、近地表附近及异常体处的电位的变化比较剧烈，所以该处的网格应较密，其它地方适当可以稀疏一些。 3) 单元阶次 00多种，适合电磁场分析模块的单元类型就有几十种，许多单元都具有线性、二次和三次等形式。其中二次和三次形式的单元称为高阶单元。选用高阶单元可提高计算精度，因为高阶单元的曲线或曲面边界能够更好地逼近结构地曲线和曲面边界， 且高次插值函数可以更高精度地逼近复杂场函数，所以当结构形状不规则、场的分布复杂时可以选用高阶单元。增加网格数量和单元阶次都可以提高计算精度。因此在精度一定的情况下，用高阶单元离散模型时应选择适当的网格数量，太多的网格并不能明显提高计算精度，反而会使计算时间大大增加。为了兼顾计算精度和计算量，同一模型可以采用不同阶次的单元，即精度要求高的重要部位用高阶单元，精度要求低的次要部位用低阶单元。 4) 网格质量 它是指网格几何形状的合理性