电导率分块线性变化的二维高频大地电磁法有限元正演模拟

为了适应电性参数在水平方向和垂直方向是连续变化的实际需要,这里采用矩形剖分,线性插值的有限元法研究了电导率分块线性变化的高频大地电磁高精度,快速正演模拟。首先给出了二维高频大地电磁测深的变分问题以及电导率分块线性变化时的有限元数值解法,编制了相应的有限元程序。利用该程序对层状模型进行了计算,并与解析解的结果做了对比分析,证明了程序的正确性和有效性,然后对典型电导率随深度线性变化的模型和典型地堑模型进行了正演模拟,结果表明能够有效地解决电导率在水平和垂直方向分块线性变化的高频大地电磁正演问题。     

MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用

在矩形有限元和三角形有限元方法大地电磁测深（ＭＴ）二维正演计算的基础上，设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序，使得ＭＴ二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深（ＣＳＡＭＴ）ＴＥ极化方式的二维正演计算。文章给出了在不同模型上计算的ＭＴ和线源ＣＳＡＭＴ的响应，分析了计算ＣＳＡＭＴ响应时，计算点离源的距离对视电阻率和相位值计算的影响，比较了ＭＴ和ＣＳＡＭＴ网格剖分的特点，提出了计算ＣＳＡＭＴ响应的改进的边界条件设计     

基于二次插值的线源可控源有限元数值模拟

在准静态近似条件下,采用矩形网格单元和双二次函数插值就频率域二维线源边值问题进行了有限元数值模拟。在二维地电条件下,给出了边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成,最后通过求解复系数方程组得到了地表视电阻率响应。引入伪delta函数模拟线源,消弱了源带来的奇异性。通过与均匀大地以及层状介质模型的解析解对比,平均相对误差分别为0.71%和1.12%。建立了两个异常体模型,数值模拟表明异常响应比较明显,为进一步实现三维可控源电磁法有限元数值模拟提供了基础。     

频率域航空电磁三维矢量有限元正演模拟

目前有限元技术的开发及在电磁勘查技术中的应用已日趋成熟.然而,有限元正演模拟主要集中于地面和海洋电磁,航空电磁三维有限元模拟尚没有受到足够重视.以前人工作为基础,利用结构化网格实现了航空电磁系统的三维矢量有限元正演模拟.从二次场双旋度矢量非齐次亥姆霍兹方程出发,应用广义变分原理推导出变分方程,并采用六面体单元剖分,将场置于单元棱边上,对每个单元应用线性插值,最后合成含有稀疏矩阵的线性方程组.针对航空电磁多源性问题,利用MUMPS(multifrontal massively parallel sparse direct solver)直接求解器进行求解,在保证精度的前提下大幅度提高计算速度.利用单个异常体模型检验算法的精度和软件的稳定性,进而通过典型地电模型的模拟验证算法的有效性.对不同地下电性结构正演模拟结果进行对比分析,进一步研究了覆盖层和垂直接触带等典型构造对航空电磁响应的影响特征.      

电导率连续变化二维海洋大地电磁场自适应有限元正演研究

岩矿石在形成过程中受湿度、含水率、温度、压力等多种外界因素的影响,这些因素会导致岩矿石的电导率连续变化,于是电导率连续变化地电模型可能更符合实际地质情况。对电导率连续变化二维介质大地电磁场分布规律,进行了有限元正演数值模拟研究。基于自适应有限元模拟技术和非结构网格剖分技术,使得新的算法可以真实地模拟电导率连续变化的复杂二维地电模型,并能够提供高精度的数值解。算例表明,海底近地表垂向电导率发生变化能够对高频大地电磁场响应产生较大影响。     

基于非结构化网格的大地电磁2D自适应有限元正演模拟

有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。     

起伏地形频率域可控源电磁二维快速正反演

起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。     

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础      

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深 CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

基于cole-cole复电阻率模型的线源可控源有限元数值模拟

受激发极化效应的影响,大地中岩(矿)石电阻率是一个与频率、电化学性质等因素有关的复电阻率。接收到的电磁场信号中,包含了电磁感应效应以及激发极化效应双重因素的影响。将cole-cole复电阻率模型引入二维线源可控源正演理论,加载伪delta函数,实现了含激发极化效应的线源可控源有限元数值模拟。计算了不同极化参数的二维地质体线源可控源的影响特征。结果表明:线源可控源响应对IP效应反映明显,视电阻率出现增高的特征;低阻高极化相对于高阻高极化模型受IP效应影响更大,从低阻中提取IP效应更加乐观。该研究为实现从复杂地质条件下CSAMT信号提取IP效应奠定了基础。     

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

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

电导率分块线性变化二维地电断面电阻率测深有限元数值模拟

在徐世浙等（１９９５）工作的基础上,给出了二维地电断面电阻率测深的变分问题以及电导率分块线性变化时的有限元方法数值解,并计算了电导率随深度变化时的地电断面的测深曲线．结果与解析解相比,精度很高．     

五极纵轴激电测深三维有限元正演模拟

本文研究了五极纵轴激电测深三维有限元正演模拟方法。首先,从三个点电流源总电位的边值问题出发,导出了异常电位的边值问题,证明了与异常电位的边值问题对应的变分问题。然后,基于三维连续电性介质模型,推导了有限元法求解变分问题的计算过程,编制了五极纵轴激电测深的三维模拟程序。最后,通过模拟水平层状介质模型的五极纵轴激电测深曲线,发现最大相对误差小于0.25%,耗费时间为18 s;从精度和效率的角度验证了算法和程序是正确的,为后续模拟复杂地电模型五极纵轴激电测深曲线、分析曲线异常特征规律以及开展反演工作奠定了基础。     

二维输运方程高精度数值模拟

采用剖开算子法,把二维输运问题剖分为两个子初值问题(对流分步、扩散分步)。在任意三角形网格中,分别对不同性质的算子采用各自适合的算法,即采用特征线法求解对流分步,采用半隐式有限元法求解扩散分步。重点探讨了对流插值问题,给出了一种完全对称三次插值模式,有效地减少了数值阻尼。为了克服高阶插值数值震荡问题,计算中保证了函数及其一阶偏导数连续。算例表明,数值方法模拟结果与精确解吻合较好。该算法在求解输运方程(包括纯对流输运方程)时,既能有效减少数值阻尼,也能保证计算中不出现数值震荡。     

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

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

基于MeshPy的3DGIS与三维有限元数值计算无缝耦合方法

3DGIS与有限元数值计算耦合方法及其系统开发,是目前GIS与岩土工程交叉学科的重要研究内容。通过对MeshPy、GRASS GIS、vtkPython以及Scipy等多个开源库的探索与编程实践,以Python为“黏合剂”,构建了GIS与三维有限元数值计算无缝耦合的平台框架。在满足几何形状、Delaunay特性、角度、面积以及体积等约束条件下,通过MeshPy包含的Triangle与TetGen库获得能用于有限元数值计算的三角形与四面体网格,并阐述了限定Delaunay三角化（CDT）节点插入算法与四面体网格质量控制方法。以排土场为例,针对CDT网格建立的地表模型,结合vtkPython封装的模块与算法,提出了根据节点分布位置综合采用线性与三角形重心坐标的插值方法,有效消除了其他插值手段导致的部分地表高程畸变的现象。借助Scipy提供的矩阵运算模块,开发了邓肯-张E-? 模型,给出了采用中点增量法进行非线性求解的算法,并通过对室内大型三轴剪切试验的数值模拟验证了程序的正确性。最后,对准东露天煤矿北部排土场堆排过程开展三维有限元数值计算,结合GRASS GIS对排弃物料高度与厚度空间分布规律的分析结果,对排土场边坡位移变化规律进行初步研究,结果表明,排土场垂直位移与排弃料高度及厚度相关,而最大水平位移则出现在各台阶坡体边缘位置,与实际情况一致。     

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

用有限元法进行电导率分块连续变化的二维大地电磁的模拟。本文第一部分给了有限元法系数矩阵的计算公式，可用于编制有限元法的计算程度。模型上的计算结果和计算精度的分析交城本文的第二部分介绍。     

基于边界元法的边坡矢量和稳定分析

矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。