点源二维势场问题的边界元法中点源处理方法

用傅氏变换将点源二维地电断面的三维边值问题,变换成二维问题.提出了波数选取的方法与确定傅氏反变换加权系数的解析方法.研究了提高计算精度和计算速度的方法技术.方法简单,精度较高.     

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

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

点源二维地电模型反演方法的研究

为了提高电法勘探数据解释的能力,寻找较好的二维反演方法,本文在研究改进的阻尼最小二乘法在二维反演中应用的同时,成功地引进和发展了Box-Kanemasu反演法。在微机上实现了电法勘探数据的二维反演解释。理论和实际地电模型的二维反演算例表明,该方法的理论是正确的,方法是切实可行的。     

用边界元法进行二维重力异常反演

重力异常正演计算的边界元方法是利用格林公式将通常的沿整个密度体的体积分转变为沿密度体边界的积分,这样使问题的维数降低了一维,从而加快了计算速度、提高了计算精度。本文利用边界积分公式解反问题,这样减少了求解参数,提高了稳定性。理论模型计算表明,该法适用于形状复杂均匀密度体的反演,特别是它能用于地下起伏密度界面的确定。     

点源二维电阻率测深边界单元法

本文重点阐述了点源二维电阻率测深边界单元方法原理、方法技术。针对电阻率测深装置特点。提出的区间截断剖分法,有利于节省计算机内存,提高计算速度,降低成本。文中用边界元法计算了若干地电形体上对称四极测深和等比装置的ρs测深曲线,获得了一些有益的结果。     

水平地形二维地电断面点源电场的正演

本文给出水平地形下二维地电断面的电阻率及面激发极化异常的数值计算方法。首先用镜象法将半空间问题变为全空间问题,然后用付氏变换将点源二维地电断面的三维电场边值问题变换成二维边值问题。再用格林公式将二维边值问题转变成积分方程。用边界单元法解积分方程,并由付氏反变换,获得三维空间中的电位。用本方法计算了几种模型的电阻率异常及激发极化异常,并与模型实验对比,效果是良好的。与起伏地形二维地电断面的边界单元解法比较,本文所需的计算机内存减少很多,计算时间节省几倍。     

点源二维起伏界面二层介质边界元解

本文以二层介质下界面上隆模型为例,讨论了用边界元法在解决层状介质中存在二维异常体情况下点源场视电阻率计算问题。文中对基本解的选取避免了对水平层状界面的积分计算,从而大大减小了所需的边界剖分点和未知量的数目,大量地节省了计算机内存,改善了计算精度,使其从解决二层介质中横向不均匀性问题推广到解决三层以上介质横向不均匀性的问题。结果表明,边界元法对本问题的研究是成功的。     

用边界元法计算矿体充电电场

理论计算研究表明,对于复杂条件下任意形状的充电矿体的充电场的研究,边界单元法是一种经济有效的数值计算方法。本文从位函数满足的点电源位于地下矿体上所建立的稳定电流场的微分方程出发,论述边界单元法数值计算矿体充电电场的方法原理和技术,并给出部分算例。     

利用有限元及边界元法计算复杂磁性体磁场

近年来,利用有限元及边界元法对具有复杂边界条件的电磁场问题的计算有所发展,可使我们沿此途径对复杂磁性体的磁场特征的认识进一步深化,国内学者在这方面曾做过工作。在前人基础上,本文作了系统的推导和计算。下面先介绍有限元法的应用,再以此为基础,介绍结合边界元法求解磁性体磁场的计算。     

点源二维电阻率法有限差分法正演计算

在寻找金属和非金属矿产,解决水文,工程地质问题和地热资源的勘探中,电阻率法是一种重要的物探方法,但由于该方法计算复杂,正演、反演问题往往停留在定性阶段,从而大大限制了该方法的应用效果和解决问题的能力,我们有必要逐步实现电阻率法资料解释的定量化。如果我们已知一个地区的地电资料,对给定的装置和供电电流值计算电异常值或视电阻率值,这个问题即是电法勘探的正演问题,它是电阻率法定量解释的基础。由于地电条件的复     

不同地电条件下电阻率法的边界元法数值模拟

不同地电条件下电阻率法的边界元法数值模拟田宪谟，黄兰珍（成都理工学院）边界元法用于不同地电条件下电阻率法的正演模拟具有许多优点：该法可以降低位场边值问题的维数，因而其与区域型解法相比较，具有离散性方程组的阶数低、原始数据量少、准备工作简便、工作效率高...     

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

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

边界元法计算曲化平精度的改进

解算“曲化平”问题,以二维情形为例Δμ(x,y)=0 μ|Ω=f(x,y)|Ω这里Ω为一段起伏不平的地形线,r=((x~2)+(y~2))~(1/2)。用边界元法求解参见。本文提出如下的改进意见: 一、边界元计算“曲化平”此类问题多出现在边界测点上精度不高,而远离边界测点精度较高。这是因为地表面异常值的测点边界往往还不是真正的边界,在测点边界以外,较小的异常值还会在一个相当的范围内存在,从而影响着边界附近测点的异常值。把它们考虑进去,会改善边界附近的地改值以及异常值较小的测点上的地改值。     

用边界元法求积分方程数值

本文给出了用边界元法求积分方程的算法。通过算例验算表明，这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。