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

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

有限单元数值模拟在张家界风景区地质调查中的应用

对张家界风景区进行1/25000遥感地质调查,解译出1031条线性构造,经野外验证可靠性达85％,新确认该区存在北北东向褶皱,提出了该区经历了两次以上的褶皱构造运动的新认识。认为北北东和北东东向构造是该区的基本构造骨架,设计了区域性印支期和张家界风景区印支期、燕山早期、燕山晚期四种应力场模型,采用有限单元法进行构造应力场反演计算,所获结果与野外调查结果吻合。认为砂岩峰林地貌是泥盆纪石英砂岩处于宽缓向斜扬起端及后期背斜叠加使地层趋于水平的部位才利于砂岩峰柱的稳定形成,多次构造应力作用产生的高角度裂隙有利于峰林地貌的造型,挽近时期地壳抬升是峰林地貌形成的动力。     

迁安盆地地下水资源有限元数值模拟与评价

依据迁安盆地的水文地质模型,建立数学模型。通过有限单元法对迁安盆地内水源地的地下水资源评价,并对其开采量和开采后的水位进行预测,同时与水均衡法求解的结果进行比较,以论证该地区地下水资源开采的保证程度,并对中长期取水要求提出合理的方案。     

双频激电法有限元数值模拟研究

利用有限单元法和Cole-Cole模型对双频激电法进行数值模拟研究。首先简要介绍了双频激电法的基本原理,在忽略电磁效应的情况下,给出了电场的边值问题和变分方程,将模型中大地和极化异常体的电阻率用Cole-Cole模型的参数替换,通过解方程,最后得到模型节点上的复电位值。通过改变地电模型的参数,得到了在不同参数下视幅频率和绝对相位的响应曲线。模拟结果与实际情况符合,表明该方法是正确和适用的。     

淄博市王旺庄水源地地下水开采预测的数值模拟

王旺庄水源地是淄博市晴纶工程的后备水源地,设计供水能力20 000 m3/d.利用Ritz有限元数值分折方法,结合单纯性线性规划,在完成了对山东省淄博市临淄区王旺庄-朱台地段的淄博市腈纶工程水源地的水文地质模型建造的基础上,成功地对水源地进行地下水开采模拟,同时对该水源地的未来10年开采作了两种大气降水条件下的地下水动态预测.为腈纶工程水源地地下水开采提供了充分的设计依据.     

频率域激发极化法有限元数值模拟

在频率较低和忽略电磁效应的情况下,利用有限单元方法和Cole-Cole模型对频率域激发极化法进行数值模拟.首先在三维地电条件下,给出电场的边值问题和变分问题,将Cole-Cole模型的频率响应引入到地电模型中,运用有限单元法对模型进行单元剖分、插值、积分和总体合成等,通过解方程最后得到表征频率域极化强度的参数幅频率.通过改变地电模型的参数,得到了不同的幅频率响应曲线.模拟结果与实际情况符合,表明该方法是正确和适用的.     

探地雷达有限元数值模拟及衰减特性探讨

从Maxwell方程组出发,推导了带衰减项与不带衰减项的探地雷达(GPR)有限元波动方程。利用matlab平台编译了二维GPR有限元正演模拟程序。通过对典型模型例子的模拟分析,给出了GPR电磁波衰减比的概念及其计算公式,加深理解了GPR电磁波在不同介质中的传播特性。在GPR正演模拟中,必须考虑衰减项的吸收作用,得到的模拟结果更符合GPR电磁波在地下介质中的传播规律,更有助于理解和掌握地下介质中探地雷达波的传播规律,对探地雷达资料的地质解释具有更好地指导意义。     

软土地下工程抗震数值模拟的若干关键问题

针对我国地下工程建设项目抗震安全性能设计和评估的需求,并基于目前缺少完善的地下结构抗震分析方法和专门的地下结构抗震设计规范的现状,重点阐述了采用有限元方法进行地下工程数值模拟和分析时需要特别关注的三方面问题：计算范围、土一结构动力相互作用和边界条件。着重分析了有限元网格的划分、土一结非线性接触问题和不同边界条件的适用性,以及这些建模条件对计算结果可靠性的影响。对基于数值手段开展地下工程抗震安全性能的计算、分析和评价,制定地下结构抗震设计规范或规程,可提供有益参考。     

适合地形起伏的二维有限单元数值模拟技术

在有限单元法矩形～三角形单元剖分的基础上,提出了适合复杂地形条件下的网格剖分方法:平行四边形～三角形剖分方法。计算结果表明,该方法具有快速、高效的特点。只要对原有程序稍加修改,即可实现对复杂地形条件下的复杂地电断面进行自动剖分和网格节点自动编号。因此,对于复杂地形条件下的地电模型,该方法能快速地实现较为精确的数值模拟。     

数值模拟软件FLAC及其在地学应用简介

拉格朗日元法是一种常见的数值计算方法,以有限差分程序FLAC软件为例,简单论述其原理,分析了它的优缺点并叙述其工作流程,表明其非常适合解决非线性、大变形问题.同时通过介绍了在国内外地学方面的一些应用实例,对数值模拟技术在国内地学方面的应用和研究提出一些建议.     

三维地电断面激发极化法有限元数值模拟

用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.      

伪随机激电法的有限单元法数值模拟研究

利用有限单元法和Cole-Cole模型对伪随机激电法进行数值模拟研究。首先简要介绍了伪随机激电法的基本原理,然后在忽略电磁效应的情况下,给出电场的边值问题和变分方程,将模型中大地和极化异常体的电阻率用Cole-Cole模型的参数替换,通过解方程最后得到模型节点上的复电位值。通过改变地电模型的参数,得到了不同参数下幅频率的模型相应曲线。模拟结果与实际情况符合,表明运用有限单元法模拟伪随机激电法中的参数是正确和适用的。     

基于区域分解法的地下水有限元与边界元耦合模型—淄博市王旺庄水源地地下水数值模拟

区域分解法（DDM）是20世纪90年代兴起的一种求解偏微分方程的新方法，方法本身独到的耦合思想和高效的并行计算机理，对于求解复杂的、大型的地下水问题具有相当的优势和广阔的应用前景。本文以淄博市王旺庄水源地地下水流模型为例，应用重叠型区域分解法（DDM）构造了边界单元法（BEM）与有限单元法（FEM）耦合模型，在两种数值方法各自优点的基础上，更形象地再现了实际水文地质原型，有效地消除了人为边界造成的流场失真。     

有限单元法与有限差分法在MT一维正演模拟中的对比分析

首先从电磁场所满足的麦克斯韦方程组出发,介绍了大地电磁测深正演的基本理论,并针对一维大地电磁模型加以讨论.运用有限单元法及有限差分法分别推导了大地电磁测深一维正演算法,并运用Matlab 7.0软件编写了相应的程序.为了检验这两种一维正演算法的准确性,设计了均匀半空间模型和层状介质模型,并给出了由本文程序、解析解得到的相应结果和图件,从而对正演结果进行对比分析.结果表明:两种方法的正演结果均真实地反映了模型的地电参数.     

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

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

物探方法在区域地质调查中的应用

物探方法是利用地下岩土体的物理性质差异,使用仪器测试地下介质的物性数据,利用计算技术,反演推断出地下不同岩土体分布的一种无损快速的探测方法。不同的岩土体具有不同的密度、导电性、磁性等物理性质,根据不同的物理性质、采用不同仪器设备和方法进行测试,可以解决不同的地质问题。本文论述了高密度电法、重力、磁法在区域地质调查中解决基岩埋深、地下水情况、隐伏断裂位置及特征、矿产分布等地质问题的情况,选择适当的物探方法可以取得良好的地质效果。     

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

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

起伏地表复电阻率法有限元数值模拟及异常特征

为了解起伏地表对复电阻率法数据的影响,发展了起伏地表条件下复电阻率法2.5维有限元数值模拟方法。经数值模拟结果与解析解对比表明,该模拟方法在计算精度上得到了较满意的结果;几种常见纯地形异常的数值结果表明:起伏地形对复电阻率法数据的畸变特征与地形坡度大小以及地形分布范围有很强的相关性,坡度越大,地形影响越为剧烈;坡度角为45°的山脊地形、山谷地形和斜坡地形所引起的复电阻率法实分量异常分别达到背景值的21.2倍、3.9倍和3.4倍,虚分量异常分别达到背景值的2.3倍、2.2倍和2.1倍,振幅异常分别达到背景值的21.2倍、3.9倍和3.4倍;而相位受地形影响非常微弱。这些工作为复电阻率法的地形改正和反演奠定了基础。     

时间域激发极化法非结构化三角网格有限元正演模拟

这里分析了传统网格剖分方式的不足,提出采用非结构化渐变三角网格对求解区域进行离散,采用固定带宽LDLT分解法求解线性方程组。设计了非结构化渐变三角网格自动剖分程序,并针对固定点源测深装置的激发极化2.5维有限元正演模拟程序。通过理论模型试算,验证了程序的正确性与有效性,并对正演模拟的速度和精度进行了初步评价。     

广域电磁法三维有限单元法模拟研究

本文以Maxwell方程组为理论基础,从矢量位和标量位方程出发,采用有限单元法对广域电磁法进行了三维数值模拟研究.在模型算例中,通过具有解析解的均匀半空间模型对本文算法的正确性进行了验证,在此基础上,设计了两种不同深度的低阻异常体模型,对广域电磁法的探测深度进行了模拟研究,数值计算结果表明:广域电磁法不受场源中近区的影响,能够在全区范围内进行观测;此外,还具有探测深度大的特点,为深地资源探测提供了理论与方法支撑.