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

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

第二类非线性Fredholm型积分方程数值解

配置法研究了地球物理中常见的第二类非线性 Fredholm 型积分方程的数值解法,将第二类非线性 Fredholm 型积分方程转化为非线性代数方程组进行求解,采用高斯数值积分公式,给出了数值计算的具体实例.利用Matlab软件的符号运算功能编程计算,克服了非线性方程难于变成求解的困难,数值例子表明该方法编程简便有效.对非线性积分方程和非线性代数方程组的求解都有重要价值.     

点源三维地形视电阻率的计算

利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。     

地下多个异常体电磁散射的数值模型

本文针对地下两孔电磁法的实际,建立多异常体(包括非均匀异常体)散射的数值模型。给出一种使计算机程序具有一定通用性的处理分块序号与坐标关系的方法。用矩量法解体积分方程,并采用合适的近似公式,计算出不同工作频率下各种异常体散射的场分布。由计算结果可得出一些对解释实测数据有参考价值的结论。     

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

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

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

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

积分方程法求解三度极化体的激发极化效应

最近十多年来,在电法勘探的正演理论中引入了一种新的求解场的边值问题的方法——积分方程法。在求解积分方程时采用了将积分方程化成N个未知量的线性代数方程组的矩量法,本文较详细介绍这种求解方法的基本原理     

计算任意形体磁异常的边界元法

本文提出一种新的数值方法,计算任意形状磁性体的磁异常。首先将磁位的边值问题转变为磁性体表面的积分方程。然后用边界单元法解积分方程,求得磁性体表面的磁位及其法向导数。根据这些数值,即可计算空间任一点的磁场强度。与有限单元法相比,本方法所需的计算机内存少。原始数据输入较简单,算例表明,计算精度高。本方法可在微型或小型计算机上实行三维运算。     

圆柱体点源场有限单元模拟

我们对有限长度园柱形电性不均匀体的点源场的有限单元法数值模拟进行了研究和讨论。利用“镜像法”原理将三维数值模拟问题转化成二维问题,从理论上严格地推导出了求“旁测剖面”电位的公式。依据此算法编制出的程序对各种模型进行了计算,结果令人满意。     

层状地基中单桩负摩擦问题积分方程解法

利用Biot固结理论和积分方程方法研究了表面有堆载的层状地基中单桩负摩擦问题。根据层状饱和土的圆形载荷基本解得出了单桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化。再结合传递和刚度矩阵传递到各个层中去,对变换域内的积分方程采用Schapery 逆变换方法得到时域内单桩的近似积分方程。求解积分方程组并进行相应的数值逆变换,就可得出层状地基中的单桩在表面圆形均布载荷作用下的位移、轴力、孔压和桩侧摩阻力随时间的变化情况。计算结果表明,桩侧剪力和孔压分层明显。     

点源三维层状介质边界元解

本文以二层介质中存在球体的模型为例，讨论了用边界元方法解决层状介质中存在三维不均匀体情况下点源场视电阻率的计算问题。文中通过对基本解的选取，避免了对水平界面的积分计算，并且在积分方程中消去了位的法向导数，使得积分方程的形式更为简单，也大大地减少了所需的边界剖分节点数和未知量的个数，节省了计算机内存，提高了计算精度。计算的结果表明，边界元法对本问题的研究是成功的。     

标量拟解析近似方法模拟立方体的异常电场

拟解析近似方法是一种求解积分方程的一种近似方法,它可以处理强散射或者大扰动的电磁散射问题,在计算过程中避免了传统微分数值方法解决问题时所遇到的大型矩阵或大型代数方程组的求解。孙建国[5]将其引入直流电场的积分方程中,并给出了求解异常电场积分方程的标量拟解析近似公式。在以前的研究中,已经验证了均匀场中异常球体的拟解析近似解的精度,这里对均匀场中的立方体异常体进行数值模拟,得到了直流电场中异常立方体模型的标量拟解析近似解。由于复杂地电模型可以用立方体的组合进行模拟,因此对立方体异常电场拟解析近似解的研究,为三维直流电场中复杂地电模型的快速正反演模拟打下了基础。     

越流含水层井群系统地下水流的边界元数值模拟

本文提出一种求解越流含水层井群系统的边界元方法。该法把抽水井井壁作为内边界,通过积分处理,不需要对井壁内边界进行剖分,通过推导求得水头H及其法向导数的边界积分方程,离散化后得到求解该积分方程的一组线性代数方程组、该法用于实例计算表明,计算精度较高。     

混合法计算三维非均匀介质中的电磁场

介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。     

有耗媒质中线天线激励下三维异常体的散射

本文研究均匀有耗媒质中,在半波长偶极天线的场激励下,多种形状和参数的异常体的散射特性。用矩量法分别求解良导电球体散射的磁场积分方程和有限导电的长方形异常体散射的电场积分方程。讨论了与数值法有关的近似、误差和检验等问题。通过数值计算得出媒质电参数改变、散射体大小、形状以及空间位置改变时散射场变化的规律。水中金属球和介质块散射的测量结果与数值计算结果吻合较好,验证了数值模型的正确性。     

基于三角元法计算三度体球冠模型的重力异常

提出利用高斯公式将体积分变成包围该积分域的全表面积,应用一系列三角形拟合计算表面,然后通过格林公式把对每个三角形的面积分变成对三角形每边的线积分之和的三角元法来计算三度体球冠模型重力异常,同时给出了合理的计算方案。计算了均匀密度球冠模型和非均匀密度球冠模型在水平地表面的重力异常,得到了2种模型的平面异常图。最后通过选择均匀密度球体作为误差对比的对象进行误差分析,验证了该方法的可行性以及计算结果的可靠性。     

工程中渗流问题的边界元算法与分析

针对河堤、水坝、高速公路档水墙等存在多孔单元介质中的定常渗流问题,将渗流Laplace方程,通过数学变换转化为边界积分方程,进而离散为线性代数方程组。采用边界元编写FORTRAN程序。计算结果表明:该法与有限元法相比,占用CPU时间短,不需要网格化分,使繁杂问题简单化,同时很好地切近实际。      

一种新的二维电阻率模拟方法

目前，对地下介质作电阻率成像需要将研究区域网格化。本文以格林函数为数学工具导出了网格化地电模型二维问题的边界积分方程，该方程适合于电阻率成像。用边界单元法求解该方程可以得出空间电位分布，从而进行视电阻率模拟。计算与实验结果的对比表明，这种方法是正确、可行的。     

复合单元法在求解有自由面渗流问题中的应用

论文指出了在有限元分析过程中地下水自由面上、下的区别,提出可采用复合单元法处理常应变交截单元,并以四面体单元为例,推导出了自由面与单元交截的截面方程,详细说明了复合单元法的原理、系数矩阵的求解方法、推导出体积分和曲面积分的求解公式。     

层状介质中含水岩溶的数值模拟

针对用均匀高阻介质中的低阻体来模拟含水岩溶构造的缺陷,用积分方程法对两层介质中的岩溶构造进行了数值模拟,并分析了计算曲线的特征,以及表层、岩溶与围岩的电阻率差异对探测结果的影响。计算结果对野外施工和资料解释具有指导意义。