边界单元法在电法资料解释中的应用

边界单元法应用于地形和地面断面的电位计算，本给出了二维点源场中若干个地形和地电断面模型的三极，联剖ρｓ曲线。应用边界单元法对某金矿０线剖面进行了实际计算和地形改正。     

点源场电阻率法二维地形改正的边界元法

边界元法是一种新的数值计算方法。本文讨论了该法的基本原理。由于边界元法可降低数值模拟的空间维数,因而在地形改正中,边界元法与区域型解法(例如。有限单元法和有限差分法)比较,该法具有基本方程组阶数低、输入数据少、精度高、速度快等优点。本文用解析法验证了角域地形上视电阻率数值解的正确性。实践表明:用边界元法获得的地形改正的视电阻率曲线突出了矿异常;点源场电阻率法2—D地形改正的边界元解法是一种经济有效的方法。     

适合地形起伏的二维有限单元数值模拟技术

在有限单元法矩形～三角形单元剖分的基础上,提出了适合复杂地形条件下的网格剖分方法:平行四边形～三角形剖分方法。计算结果表明,该方法具有快速、高效的特点。只要对原有程序稍加修改,即可实现对复杂地形条件下的复杂地电断面进行自动剖分和网格节点自动编号。因此,对于复杂地形条件下的地电模型,该方法能快速地实现较为精确的数值模拟。     

基于二次场二维起伏地形MT有限元数值模拟

通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。     

电导率分块连续变化的二维MT有限元模拟（Ⅱ）

用有限单元法进行电导率分块连续变化的二维大地电磁场的模拟，本文的第一部分给出有限单元法系数矩阵的计算公式，这部分给出四个模型的计算结果。在有限元中，分别采用双线性播值和双二次插值，计算结果与解析法或有限差分法进行对比，吻合得很好。算例表明，本方法具有计算精度高、网络剖分简单、模拟功能强的优点。     

任意三维地电模型视电阻率异常的边界元法数值解

边界元法(BEM)是一种新的数值方法。由于该法可以降低所研究场问题的维数,因而边界元法较之域型法(FEM和FDM)具有应用简便、数据量少、计算快、精度高等优点。从而,用该法解决了域型法难以实现的三维地电模型视电阻率异常的计算问题。本文论述了用边界无法求解点源场地表水平和起伏下三维地电体位场问题的方法原理和数值处理技术,并给出了若干算例:导电球状矿体上视电阻率数值解与解析解结果;导电球状矿体上不同测线上视电阻率平剖曲线;三维山脊地形下导电椭球状矿体上视电阻率联剖曲线及其地形改正结果等。由本文内视电阻率的边界元法数值解与解析解结果对比的一致性和三维地形校正的显著效果,表明了边界元法是求解任意三维地电模型上位场问题的有效方法。由于用该法实现了对任意三线地质体上空间位场计算,必将推动三维电法勘探工作的开展与深入研究,无疑会对提高电法勘探的地质效果发挥重要作用。     

水文地质计算中的有限单元和边界单元交连耦合法

 本文提出的有限单元和边界单元交连耦合法(CFEBEM)可作为地下水流场和水文地质参数正反演计算的有效工具。该方法改进了Liggett(1983)的代入耦合法,通过内边界节点组和交连方程,实现有限单元和边界单元的多区耦合。CFEBEM方法发挥了有限单元和边界单元各自的特长,具有单元区划分灵活、方程组装简便和计算精度高等特点,可有效地解决结构复杂含水层及大型密集抽水条件下的流场计算和水文地质参数反演。在太原岩溶地下水资源评价和合理开发研究中的应用,证明了这种方法的有效性。     

基于ILU-BICGSTAB算法的起伏地形MT二维正演模拟

针对起伏地形对异常体的探测效果的影响,利用有限单元法,基于MTLAB平台编写了大地电磁二维正演程序,为更好地模拟实际地形,采用矩形网格中进一步三角细化的剖分方式,并利用不完全LU分解的预条件BICGSTAB算法提高计算速度与精度,最后通过建立多种地电模型分析讨论了起伏地形对大地电磁场的影响,验证了该方法的实用性与准确性。     

边界单元法点电源二维地形改正的PC-1500机程序

边界单元法是一种较新的数值模拟方法。目前,在国内外发展较快,在地球物理勘探领域也有日益扩大的势头,理由是明显的。与有限单元法相比,边界单元法具有以下优点:①精度高;②占有计算机内存少,仅需输入二维地形剖分点的坐标就可进行地改;③输入数据准备及操作简单;④可降低数值模拟的空间维数;⑤特別适合于无界域,而很多地球物理正问题都是在无界域中进行。     

二维电法数值模拟中的网格剖分方法

本文提出一种新的有限元及有限差分等方法中的网格剖分及单元电阻率值设置算法。由此算法设计的软件，使用者只需输入地电体的边界顶点坐标及其电阻率值，程序就可自动判断及设置各单元电阻率值，对程序使用者来说极为简单方便。     

复杂地形条件下激发极化有限单元法三维数值模拟

野外激电勘查过程中,据激电异常进行钻探验证时,常发生验证目标矿体擦边而过现象.究其原因是在实际勘查中,地形的起伏变化改变了地下电场的分布规律,导致激电异常发生位移及幅值改变.为了改善、克服这种现象,必须从理论上进行研究.运用四面体单元三维有限单元数值正演计算方法对此进行了理论研究,选用典型山谷地形水平层状极化体模型进行数值模拟计算.研究结果表明:①起伏地形不产生假的激电异常,仅对激电异常的形态和幅值产生影响;②采用四面体剖分单元有限单元法研究复杂地形条件对激电异常的影响是行之有效的.     

电法测井的三维有限元模拟

研究了不同测量方式(井-地,地-井,井-井)下点源场井中电法的三维有限元数值模拟。考虑到深度方向上大范围的网格剖分和井眼的影响及井-井测量等因素,采用放射状三棱柱单元的网格剖分方式,以提高网格质量,减少剖分单元数;给出了三棱柱单元的坐标变换公式,进行精确的单元积分,减少了单元积分时间;结合非结构化网格技术,实现了复杂模型的模拟;开发出相应的程序实现了复杂条件下(如考虑井眼影响、井井测量、倾斜井情形、地形起伏等)电法测井的三维有限元模拟,数值算例验证了方法的可靠性及计算效率,并对不同情形下的异常响应进行了分析,为进一步的反演工作奠定了基础。     

电阻率层析成像及演中若干问题的讨论

以双异性点电流源电流场二维地电断面为例，着重讨论了电阻率展析成像反演中有限单元法数值模拟边界条件的选择、网格剖分以及模拟精度等问题，同时也讨论了成像反演中最优化及演方法及拟合参数的评价等。文中提出的新的网格剖分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。     

任意四边形剖分二维MT有限元正演模拟

为了模拟实际的起伏地形,在前人研究的基础上,单元网格设计为任意四边形网格,单元内的场值双线性变化,采用了高斯数值积分法计算单元系数矩阵,并给出二维大地电场的辅助场表达式。通过模型算例表明,模拟结果与解析解的均方误差在1%以内,地形模型与前人的模拟结果相吻合。分析对比了两个不同坡度模型对视电阻率和相位的影响。采用任意四边形网格剖分法,降低了编程难度,可以方便地适应野外地形的起伏变化。     

起伏地形条件下二维大地电磁各向异性正演

在大地电磁资料处理和解释中,大地电磁的各向异性正演一直是国内、外研究的前沿课题。这里首先从麦克斯韦方程组出发,得出各向异性介质二维大地电磁场的边值问题以及等价的变分问题,进而对计算区域进行完全非结构化三角形剖分,在单元内采用线性插值,将变分方程转化成线性方程组,求解出有限单元法数值解,获得各向异性介质下的电磁场值。讨论了电性主轴与笛卡尔坐标轴之间的夹角以及地形的起伏变化对视电阻率和相位值的影响,结果表明:起伏地形条件下,TE和TM两种极化模式的视电阻率和相位值都会受到影响,且TM模式受地形影响较大;由于介质的各向异性,视电阻率和相位断面图出现明显不同于各向同性的形态,导致横向分界面模糊。     

起伏地形二维位场上延与换算的边界单元法

本文用边界单元法进行起伏地形二维位场向上延拓与换算,介绍了方法的原理和具体的计算公式。以地形上的磁异常Z为例,换算出地形上的H,延拓出地形上部的Z及导数δz/δx、δx/δy。计算结果与精确解符合得很好。在有效异常范围内,误差基本小于1％。     

三维地形大地电磁场的边界元模拟方法

提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果.      

二维电法数值模拟中的网格部分方法

本提出一种新的有限元及有限差分等方法中的网格剖分及单元电阻率设置算法，由此算法设计的软件，使用只需输入地电体的边界顶点坐标及其电阻率值，程序就可自动判断及设置各单元电阻率值，对程序使用来说极为简单方便。     

水下直流电阻率法数值模拟

利用有限单元法对水下电阻率测深进行数值模拟。与地面测量不同的是，水下电极所处的深度、水体本身和水底地形都会给电阻率数据造成较大影响，给数据的解释带来困难。若先通过计算，分析不同电阻率装置在水下测量时的视电阻率与水深的关系，再用有限单元法对水下的三维地电断面进行数值模拟，讨论水底地形对视电阻率的影响．即可用比较法消除地形影响。     

三维地面断面电阻率测深有限元数值模拟

用有限单元法进行了电导率分块均匀的三维点源电场电阻率测探的数值模拟。首先给出了三维构造中点源电场的边值问题、变分问题；然后，用有限单元法求解变分问题，将区域剖分成六面体单元，在单元中进行三线性函数插值，将变分方程化为线性代数方程组；最后解方程组，得各节点的电位值，进而计算出地表的视电阻率。对几例较典型的地电模型进行试算，结果表明本方法是行之有效的。