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

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

一种含约束条件的二维自适应三角剖分新算法

为了对二维复杂地电断面进行三角网格剖分,在已有算法的基础上,提出了一种含约束条件的二维自适应三角剖分新算法,并编制了该算法的通用程序,可对任意含约束条件的二维区域进行自动剖分.通过引入节点剖分尺度函数,算法允许在指定位置设置网格剖分尺寸,自动生成渐变网格.实际算例表明,形态因子大于0.9的网格占95%以上,能很好地对带地形、含断层或地质体等二维复杂地电断面进行有限元三角网格剖分.     

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

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

频率域航空电磁法人机交互式二维正演研究

研究了频率域航空电磁法的电磁场理论、边界条件、有限元单元法求解等内容,重点介绍了网格剖分原理,实现了人机交互网格剖分软件以及二维有限单元数值模拟算法。通过二维数值模拟计算,了解了圆柱体模型的响应曲线特征,提高了频率域航空电磁数据的解释水平。     

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

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

无单元Galerkin法大地电磁三维正演模拟

无单元Galerkin法(EFGM)作为一种相对成熟的无网格方法,避免了网格剖分,其精度高,适用于复杂电导率分布和复杂边界形状的计算。本文将EFGM用于大地电磁三维正演,详述了三维EFGM形函数的构造过程,从大地电磁三维变分问题出发,利用Galerkin法结合高斯积分公式推导了相应的系统矩阵离散表达式,简述了边界条件的加载技术,研究了支持域尺寸对EFGM三维正演计算精度的影响,最后通过数值计算验证了EFGM三维算法的正确性。     

基于非结构化网格的瞬变电磁2．5维有限元正演模拟

利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。     

二度体的重力张量有限元正演模拟

介绍了重力梯度张量,并将有限单元法应用于二维重力梯度张量的正演计算。为了验证有限元正演方法的精度,对截面为矩形的两个二度体组合模型进行有限元正演模拟,结果表明正演曲线与理论曲线形态一致,拟合情况好。通过对截面形状不规则、密度分块均匀的二度体进行正演模拟,说明有限元法可通过网格剖分来逼近不规则目标体的边界,并对剖分单元赋予不同的密度值来实现对复杂二度体的重力张量的正演模拟。     

基于异常复电位2.5维CR有限元数值模拟

论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交流双点电源复电阻法的变分问题;然后利用有限单元法求解变分问题,在矩形网格中增加2条对角线形成三角网格对研究区域进行剖分,在单元中对复电位和岩、矿石复电导率进行线性插值,将变分方程化为线性代数方程组;最后,求解线性方程组,得到各节点的异常复电位值,并根据选定的测量装置计算出地表的视复电阻率,进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法的有效性后构建模型算例,给出了其不同测量装置下的计算结果。     

电阻率/激电测深二维人机交互正演模拟

在有限元点源二维电阻率正演基础上,针对传统对称四极装置其电极距跨度大且不等间隔分布等特点,提出将整个正演计算区域划分为目标区、观测区和扩展区三个区域,在观测区域先采用等间隔剖分,然后在区域中插入实际电极作为正演计算节点的剖分方式,同时将无穷远边界条件简化技术引入到正演计算中,并结合人机交互与可视化技术实现了大极距电阻率/激电测深的二维正演模拟。最后通过对相关模型的试算结果表明,作者在本文中提出的针对传统对称四极装置电阻率/激电测深的二维正演模拟技术是有效可行的。     

基于非结构化网格的大规模磁测数据正则化反演研究

磁测数据解析正演的内存需求与模型规模、观测点数成正比,以其为依托开展大规模磁测数据反演时,计算成本较大。文章利用非结构化网格有限单元法进行磁测正演,分析了地表观测网格、边界节点数和局部网格密度等因素对正演计算精度的影响;在此基础上,基于灵敏度矩阵隐式存储的高斯-牛顿法求解正则化反演目标函数,实现了大规模磁测数据三维反演。理论模型试算和铀矿区域实测数据案例表明,文章所述方法可保证正演计算的精确度,提高了计算效率、减少了反演的内存消耗,可用于复杂条件下大规模磁测数据的三维反演。     

网格尺寸及边界对大地电磁有限元正演精度的影响

大地电磁正演的精度对后续反演意义重大。决定大地电磁有限元正演精度的因素,主要有插值函数、方程组的求解误差,以及网格的剖分程度。而网格剖分的合适与否,为有限元的求解提供了先决条件。因此,从边界条件和横向网格以及纵向网格三方面,讨论了网格剖分对大地电磁正演精度的影响。经研究表明,对于一维和左右对称的二维介质,左右和下边界都无需放置在无穷远处,并且横向网格对精度影响不大。而由于纵向网格与近地表的最小趋肤深度密切相关,所以地表网格不得大于1/4个趋肤深度。     

有限单元法求解水中二维体问题正演研究

潜艇在现代海军中占有举足轻重的地位,所以反潜技术也越来越重要,相对于传统的声波、磁法等反潜技术,电磁法探测是近几年才开始起步的。这里研究了电磁法探测水下二维体的问题,基于天然电磁场的有限单元正演,同时考虑电阻率与磁导率,采用三角剖分,通过模型试算验证了天然电磁场的有限元算法在水下二维体探测中,具有较高的精度和适用性,模拟结果 TM模式优于TE模式,可以为电磁法探测水下潜艇打下坚实基础。     

复电阻率法二维正演并行计算研究

复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。     

无网格局部Petrov-Galerkin法大地电磁二维正演

将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。      

一种适用于电磁法数值模拟的二维建模与网格剖分方法

为满足电磁法二维数值模拟解释对交互建模与可视化网格剖分的需求,利用计算机图形学、人机交互、拓扑关系学等技术,设计了基于测线、测点等信息的二维矢量建模、四边形和三角形网格剖分方法与流程,并开发形成了软件模块。该方法与软件可提供起伏地形下大地电磁、可控源音频大地电磁法二维数值模拟的三角形、四边形网格(下边界随地形变化或平地形)剖分实用化工具。将该方法应用到多种模型试验中,取得了良好的建模与网格剖分效果。     

大地电磁二维正演中的无网格局部径向基点插值法

无单元Galerkin法作为较成熟的一种无网格方法,已成功应用于有限元法触及的领域,还解决了如大变形、裂纹扩展及高速冲击等网格方法较难处理的问题,但其最大的缺陷在于系统方程的离散需借助背景网格,因此该方法并非真正意义上的无网格方法。无网格局部径向基点插值法采用子域法构造系统方程,加权残量只要求在局部积分域消除,大大降低了对背景网格的依赖,向真正的无网格迈进了一大步.这里将此方法用于大地电磁二维正演,介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法推导了与之对应的无网格局部弱式系统方程,并用高斯积分将其离散化;论述了局部径向基点插值法较无单元Galerkin法及有限元法的优缺点;最后通过二维模型的计算验证了算法的有效性。     

基于泊松方程的空间波数混合域二度体磁异常数值模拟

正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。     

基于非结构网格的带地形MT二维Occam反演及应用

本文对基于自适应非结构三角形网格的带地形MT数据二维Occam反演进行应用研究.自适应非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,正演网格由有限元解的后验误差估计指导自动细化,保证了模型响应的准确性,反演网格在反演目标区域采用精细网格剖分,在模型边界区域采用粗网格剖分,在满足反演精度的前提下减少了不必要的反演...     

起伏地形下偶极-偶极激电测深二维反演软件开发及应用

在有限元点源二维正演模拟中,采用第二类齐次边界条件结合相应的网格剖分技术,并将二维正演模拟程序引入到偶极-偶极激电测深二维反演中,编制相应反演程序并进行集成,形成集数据输入、反演计算、数据成图和结果输出为一体的二维反演软件。通过对起伏地形下的两个理论模型合成数据及一例实测数据的反演,反演结果表明反演软件有效。