基于二次场算法的大地电磁二维有限单元法正演

推导了同时考虑电阻率与磁导率变化的大地电磁二维方程,并应用有限单元法进行数值模拟.为了提高计算精度与效率、简化计算节点生成,采用格林定理处理二次场方程源项,并设计实现了一种基于二叉树结构的收缩网格剖分算法;采用基于最少填入元思想的稀疏矩阵符号分析方法,实现了稀疏线性方程组的LDLT求解.利用二次场算法进行模型试算,结果表明所采用的新计算方法大大减少了计算单元数量,提高了计算精度与效率.     

基于二次场的二维大地电磁有限元法数值模拟

引入基于计算二次场响应的方法,将其用于二维大地电磁数值模拟.导出了二维大地电磁二次场的表达式,在使用有限元法解偏微分方程时,利用二次场偏微分方程与总场偏微分方程的相似性,直接采用总场法中的系数矩阵作为二次场的合成矩阵.将根据此原理所编制的二次场有限元法模拟软件,对几个典型模型进行了试算,并与传统总场模拟法作了比较分析.结果表明,直接计算二次场的结果精度更高,更接近真实解.     

不完全LU分解预处理的BICGSTAB算法在大地电磁二维正演模拟中的应用

基于双二次插值的有限单元法求解大地电磁二维正演问题,以不均匀网格剖分为基础,推导出大地电磁响应的计算公式。针对有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称正定复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难等问题,采用不完全LU分解（即上三角与下三角分解）处理的稳定双共轭梯度算法（BICGSTAB算法）求解该线性方程组,通过对层状介质和二维模型电磁响应进行计算,获得二维大地电磁的视电阻率曲线和阻抗相位曲线。研究结果表明,BICGSTAB算法具有速度快、精度高和稳定性好等优点。     

基于可控源音频大地电磁二维反演的MPI并行算法研究

可控源音频大地电磁(CSAMT)是以岩石的电性差异为基础和前提的勘探方法。该方法采用了人工场源,弥补了天然场源信号微弱、不易观测等缺点,已被广泛应用到矿产普查、油气勘探、水文环境等方面,并发挥了巨大的作用。基于CSAMT根据各频率依次独立处理数据的特点,采用了MPI并行运算方法,用多个进程同时来计算各频率数据,最后再将数据进行收集,得到最后的计算结果。通过对正演和反演的计算结果图的比较,验证了程序的正确性。对并行算法的效率进行了统计,进程数为2~5时,加速比能达到1.9~4.5,极大的提高了计算效率。     

基于二次场算法的时间域激电2.5D有限元正演

二次场算法是提高时域激电正演计算精度的有效方法之一,能很好处理场源奇异问题。研究工作首先通过引入Cole-Cole模型,推导了在此模型下复电阻率二次场满足的边值问题;并利用有限单元法进行数值计算。然后通过Guptasarma滤波算法,实现了频率域激电到时域激电的转换,改变了传统求解时域激电的方法。为了提高计算效率,首先对边界采取近似处理。经处理的刚度矩阵只与供电频率及波数有关,而与供电点位置无关。其次设计了基于图论理论的矩阵重排与填入元分析算法,实现了矩阵的LDLT分解,从而实现了基于二次场算法的时间域2.5维激电正演的快速算法。利用快速算法,计算并分析了视电阻率、视极化率异常特征;其次通过计算频散率异常特征,验证了极化率与频散率对异常体响应具有相似性。     

大地电磁二维正演响应特征试算

采用基于有限差分的大地电磁二维正演程序,分别对单个低阻异常体、高低阻复合异常体和断层两侧高阻异常体模型进行了模拟计算,对比TE、TM视电阻率和相位特征以及TP振幅的响应规律进行模拟归纳并且总结出大地电磁二维正演的响应特征,为大地电磁方法的设计和数据解释提供参考与借鉴.通过计算结果表明:在二维大地电磁测深正演响应特征中,TE、TM视电阻率和相位以及TP幅值展现出不同的形态,其中TE模式可以更好地分辨纵向地质体异常,TM模式可以更好地分辨横向地质体异常.     

正六边形二维场域电位的有限差分算法

以正六边形二维场域为例,建立了正三角形网格划分的有限差分法的计算模型,并通过MATLAB编程得到场域内电位分布的数值解和相应图形,对有限差分法的计算处理进行了讨论和分析.     

基于互相关算法的CSAMT接收信号处理研究

针对可控源音频大地电磁法（CSAMT: Controlled Source Audio-frequency Magnetotelluric）实际勘探中, 因深部目标体响应信号微弱, 测量环境噪声严重等特殊情况而导致接收信号的清晰度大大降低的问题, 提出一种基于互相关算法的CSAMT接收数据去噪处理的新方法。该方法利用发射信号与接收信号之间相关性强, 与随机噪声相关性弱的特点, 设计发射信号电流波形记录器, 将接收信号和发射信号作互相关处理, 滤掉噪声等干扰信号。实验结果表明, 该方法可压制噪声, 可提高CSAMT测量的准确性。     

基于非重叠区域分解算法的大地电磁法二维正演模拟

区域分解算法采用分而治之的思想,将大规模问题转化为若干个小问题进行求解,已成为大规模数值计算领域的常用算法之一。将非重叠区域分解算法引入到大地电磁法二维正演模拟中。首先将整个求解区域分解为多个互不重叠的子域,子域之间共享边界元素;然后对每个子域采用有限差分进行离散,采用Schur补偿算法解耦得到共享边界节点上的未知数,并作为子域问题的边界条件得到关于子域内部节点的线性方程组;最后,利用直接求解算法对上述方程组进行求解,实现了大地电磁法二维正演。该算法的准确性和可行性通过多个地电模型的对比试算得到了验证。此外,还统计分析了采用不同子域分区方式和分解个数时的计算耗时,结果表明子域分区的方式对计算效率影响不大,但子域分解个数的影响则较大,进行区域分解时需要选择合适的子域个数。     

用有限单元法模拟各向异性介质中二维电性异常体的大地电磁响应

为了研究电导率各向异性介质中二维大电磁场的变化规律,首先从其满足的变分问题出发,利用有限单元法计算出大地电磁场,借助辅助场获得视电阻率,其次讨论TEz模式下各向异性介质和各向同性介质中视电阻率响应曲线的区别,接着探讨不同层面倾角、不同频率下各向异性异常体引起大地电磁场的变化,另外就高、低阻各向异性介质对大地电磁视电阻率影响的差异进行讨论,最后讨论平行层面和垂直层面的电阻率对各向异性介质平均电阻率的影响,从而获得各向异性大地电磁场响应特征。     

有限元法在表面波等离子体中的数值模拟应用

工业用平板型表面波等离子体( SWP)源具有大面积、高密度、高活性等优良性能,其建造之前必须进行结构优化,而对表面电磁波的数值模拟,可以辅助其优化设计.基于有限元法( FEM),采用Ansoft HFSS模拟软件,数值仿真了大面积矩形SWP源的电磁场空间分布.研究了开槽天线激发表面波场的特性,得出了表面波场的空间分布,...     

打结钢丝绳有限元数值模拟方法研究

采用线性强化模型得到钢丝绳等效的弹性模量和屈服应力,将其应用在打结钢丝绳的有限元分析中。结果表明,该方法应用于打结钢丝绳的数值模拟切实可行。     

基于有限元法的高压挤压砾石充填数值模拟

从疏松砂岩的岩石特性以及施工工艺角度分析了影响高压挤压砾石充填效果的各种可能因素。在此基础上,采用有限元分析软件针对射孔井进行了数值模拟,模拟结果显示射孔沿着最大水平主应力方向,有助于裂缝的产生。在数值模拟的基础上分析了高压挤压过程中疏松砂岩的启裂机理,认为压密砂层的破裂压力即为高压挤压砾石充填的启裂压力,其裂缝形态,直接影响着近井地带的裂缝形态。     

用有限元法解热水驱油数值模拟问题

在热水驱油中，地下油藏的温度场，油、水的压力场以及它们的饱和度场均随时 间变化，并相互影响。为有效的对场变量进行数值计算。具体地建立了热水驱油的数 学模型并进行了数学处理，提出了适于迭代计算的泛函，并作了有限元离散和构造了 算法，编制了计算程序。通过一定的数值计算，说明有限元法能很好的解决油藏数值 模拟问题。     

基于动力有限元的路基动力响应数值模拟分析

为了研究路基动力响应特性,本文基于动力有限元原理,采用考虑Rayleigh阻尼的空间三维非均质线弹性体系作为结构模型,通过Newmark法求解在FWD脉冲荷载作用下结构的动力响应,并用Fortran90语言编制了相应的正演计算程序,并与实测数据及大型有限元程序计算结果进行了对比与分析,结果表明:该程序计算结果与实际结果...     

波浪绕射问题的无限相似单元数值模拟

采用无限相似单元方法与有限元法（Garlerkin）相结合对波浪绕射进行了数值模拟，求得任意截面柱状物体上的波浪荷载．模型中，对内外域交界面的选择并无限制，极大减少了内域单元的个数；外域组合刚度矩阵的阶数仅与内外域交界面上的结点个数有关，较大程度上减少了内外域拼接后总刚度矩阵的阶数．数值模拟结果和解析解吻合较好，且该方法简单快捷、收敛性好．     

利用有限元法对板坯电磁搅拌中电磁场的数值模拟

基于Maxwell方程,建立了板坯内部电磁场的数学模型.为了满足计算机求解的需要,对数学模型进行了离散化和简化处理.使用有限元法对简化后的数学模型进行了仿真,得到了电磁场等势线图和磁力分布图.通过对磁场分布进行定性分析,得到了在其他条件不变的情况下,电磁力和电流强度、电磁力和电流频率的变化关系.     

Burgers方程基于混合有限元法的差分格式及其数值模拟

给出了Burgers方程的一种基于混合有限元的最低阶的差分格式，并给出了数值解的例子，与以往的处理Burgers方程的有限差分法不同之处是该方法能同时求出速度和流通量的近似解，而且得到的数值解具有很好的稳定性。     

注塑制品翘曲变形有限元数值模拟的误差分析

从数值分析理论、力学理论和注塑成型CAE模拟的积累误差三个方面 ,系统地分析注塑制品翘曲变形有限元数值模拟中产生误差的主要原因 ,以求找到有效的改进方法和更深入研究的可行性方向 ,从而能提高模拟结果的可靠性 .