应用有限元法计算磁各向异性体的磁异常

用常规的虚磁荷法计算磁各向异性体的磁异常是十分困难的,甚至是不可能的。王书惠最近提出用有限元法计算磁各向异性体的磁异常的方法。有限元法不仅考虑到磁各向异性效应,而且还考虑到退磁效应和剩余磁性对感应磁性的影响,因而有限元法给出更准确的结果。本文简化了求解磁各向异性体的磁位的数学提法及相应的变分理论,并给出有限元法的具体计算公式。此外,我们在单元中采用二次函数插值,有限元法的运算结果可直接给出各     

边界单元法原理及其在地球物理勘探中的应用简介

边界单元法比有限单元法有重大的优点:1)方程组的规模小得多,2)向计算机输入数据少得多;3)适用于无界区域。这些优点使边界单元法更便于解决二维和三维地球物理正演问题。本文简单介绍边界单元法的原理及其在地球物理勘探中的应用。     

水平地形三维电场的边界单元解法

用格林公式将水平地形三维电场的边值问题转变为不均匀体表面的边界积分方程,然后用三角单元对不均匀体的表面进行剖分,并用高斯求积公式将积分方程转变为线性代数方程组。解方程组,得不均匀体表面的电位及其法向导数。根据电位和法向导数,可用数值方法计算地面的电位和电场强度。     

磁场垂直分量Z与水平分量H互相换算的有限单元法

由地形线上的磁场垂直分量Z计算水平分量H,对山区磁测资料的解释是很有意义的.鉴于目前有限单元法在各工程领域得到广泛应用,成为一种高效、通用的计算方法,本文提出二维磁场Z、H分量互相换算的有限单元法. 基本原理首先阐明由Z求H的基本原理. 在无源空间Ω,磁场水平分量H满足拉普拉斯方程??~2H=0,如果已知H的第二类边界条件??=f,其中n为区域Ω的边界??的外法向,f为已知函数,则求解Ω和??上的H,就等于解下列边值     

二维各向异性地电断面大地电磁场的有限元法解法

假定各向异性主轴之一与层面垂直,另一主轴与构造走向平行,本文给出这种条件下,二维各向异性地电断面的大地电磁场的变分方程。用有限单元法解变分方程。对区域进行三角单元剖分,在单元内进行二次函数插值,将变分方程转变为线性代数方程组。解方程组得各节点的场值,并可计算场的导数。由这些值可计算地表的视电阻率。 算例表明,对于H_x型波(x平行构造走向),各向异性体与各向同性体的大地电磁场有很大区别。      

有限单元法计算良导嵌入体二维模型的大地电磁测深曲线

以具有良导嵌入体的二维模型为例,研究用有限单元法求取二维介质大地电磁测深理论曲线问题,讨论了计算中的区域选取,边界条件,单元部分等问题。给出了几个算例,说明此方法可用于二维模型的理论研究和反演解释中。     

地形对大地电磁勘探的影响

本文讨论了地形对大地电磁场的影响。用有限元方法计算了几种地形的大地电磁影响曲线。计算表明,地形对E_X型波视电阻率的影响,主要出现在地形起伏幅度与电磁波在地下介质中的波长可以比较时。H_X型波的地形影响比E_X波严重,随着周期的增长,H_X型波的地形影响变得稳定。地形凸起处的影响比下凹处严重。