复杂电性结构大地电磁二维响应特征

为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。     

可控源音频大地电磁法方位非各向同性层状介质Occam反演

目前,可控源音频大地电磁法测深资料反演主要基于各向同性介质,在复杂的地电体中,由于岩石裂隙或矿物形状等原因使得在各个方向的电阻率不同。本文主要探讨更为符合实际地质体情况的方位非各向同性层状介质在正交水平电偶源激发条件下的可控源音频大地电磁数据的反演。运用Occam反演方法,实现了可控源音频大地电磁法方位非各向同性层状介质反演。构建理论地电模型,通过理论模型合成数据的反演算例检验了所实现的可控源音频大地电磁法方位非各向同性层状介质反演算法的有效性和稳定性。同时,不同数据的反演结果表明：方位非各向同性介质的反演应同时使用视电阻率ρsxy、ρsyx和相位?xy、?yx数据。     

大地电磁二维反演方法对比研究

由反演方程出发, 从构造目标函数和方程求解方面, 比较了国际上通用的OCCAM (奥可姆法)、RRI(快速松弛法)和REBOCC(简化基奥可姆法), 并由两个合成模型分析比较了OCCAM、RRI和REBOCC的反演结果。试验表明, REBOCC是继RRI之后又一种较快的反演方法, 可为电法勘探资料处理提供新的途径。而在RRI和REBOCC反演时, 网格的划分对结果很有影响, 这是值得注意的地方。     

起伏地形大地电磁二维反演

为了适应地形起伏的实际地质情况,开展带地形的最小二乘二维反演研究。鉴于大地电磁（ MT）反演的不适定问题,引入Tikhonov的正则化方法,从而获得关于总目标函数的方程,利用光滑约束最小二乘法求解总目标函数方程。由于正则化因子值与反演精度以及稳定性相关,采用主动约束平衡方法获取最优化的正则化因子,以确保反演精度和稳定性都达到最佳。与此同时,利用电磁场互易定理以节省反演迭代过程求解雅可比矩阵的计算时间。构建了若干地质构造模型进行试算,分别讨论TE、TM模式以及二者联合模式的反演结果,并与前人研究工作对比以说明本文方法的反演效果。     

用阻抗相位识别二维音频大地电磁垂向异常

从原理上分析了音频大地电磁阻抗相位在静位移和垂向电性异常体影响下的变化特性,并采用二维有限元方法进行了正演模拟。考虑到阻抗相位在静位移和真实垂向异常中的响应差异,认为利用阻抗相位开展综合解释对于识别二维垂向电性异常具有十分积极的意义。理论模型计算和工程实例表明,基于阻抗相位的综合解释方法对于识别二维垂向电性异常具有较好的应用效果。     

二维大地电磁主轴各向异性理论模型反演

系统阐述了广泛认可的二维有限差分正演算法理论,再考虑一种特殊情况,将任意各向异性情况的电导率简化为只有对角线的主轴电导率,得到电场E x与磁场H x解耦的方程组,并利用TE模式只与σxx有关、TM模式只与σyy和σzz有关的特点,使用最小二乘法将TM模式单独反演计算.在y轴方向上,反演结果较为准确地还原了各向异性异常体...     

二维大地电磁正则化共轭梯度法反演算法

针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。     

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

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

二维最平模型的大地电磁快速反演

在已有原理的基础上，本文给出了一种新的大地电磁二维快速反演方法，该方法在合理的模型网格下，求模型的纵向和横向一阶偏导数平方最小，反演得到最平滑模型，有效地消除了二维近似反演中易出现的多构造和数据过拟合现象。同时，引入总体控制函数以及用牛顿单变量寻根法求拉格朗日因子，使反演快速、稳定地收敛。理论模型试算表明，该方法比近来发展的同类反演方法效果更好，另外，野外资料的反演也取得了比较满意的结果。     

电性可控源音频大地电磁法二维反演方法技术研究

在多年来电性可控源音频大地电磁法方法实际应用的基础上,对电性可控源音频大地电磁法二维反演方法技术进行研究。用美国ZONGE公司的SCS2D软件对一个应用实例资料的处理,给出不同的计算结果以及一般的反演计算工作流程。     

大地电磁测深二维地形影响及其校正方法研究

采用有限元法实现ＭＴ二维模型的正演模拟,将地形条件下的二维地电模型的ＭＴ响应分解为地形引起的畸变场和与测点高程位置有关的稳定场两部分进行分析。用有限元法模拟二维地形条件下地电结构的ＭＴ 响应和纯地形的ＭＴ响应,采用比值法消除地形引起的干扰场,然后,将消除地形影响后的视电阻率值通过延拓方法获水平地形条件下二维地电结构的视电阻率分布,完善大地电磁测深视电阻率曲线的地形校正。     

二维大地电磁正演中的无网格算法

无网格法是近几年来发展的一种新的基于变分原理的数值计算方法,由于在计算形函数中不需要划分网格.在力学、电磁学等领域得到了广泛的研究.基于无网格法在大地电磁勘探正演中的应用进行了研究,首先对无网格法的基本原理进行了阐述,并利用广义变分原理推导出了相应的离散方程,编制了相应的程序,最后通过两个理论模型的计算结果检验了算法的正确性.     

大地电磁资料的分形特征

分形理论是研究一类无规则、非线性但具有自相似性的无序系统。这个无序系统的维数一般是分数。分数维的大小可以反映系统的复杂性，揭示系统的本质规律。本文利用功率谱分析法来揭示大地电磁资料的分形特征，并简单分析分数维在大地电磁中的地质意义。发现当地下电性分布均匀时，电、磁场分量的自功率谱都服从对数坐标轴上的线性分布；而遇到电性不均匀体时，功率谱的线性发生畸变。     

二维大地电磁资料的人机联作解释

本文介绍了在PC微机上实现的大地电磁资料的二维人机联作解释软件(简称MTI2D)。该软件由建模、二维正演和绘图三部分组成。建模、绘图部分采用TURBO BASIC编程,二维正演部分采用有限元解法,用FORTRAN77写成。三者相互独立,通过SHELL功能将三者连在主莱单下。此外,还说明了在建模过程中涉及的几个问题。 经理论验证和实测资料解释表明,该软件使用方便,操作简单,运行稳定可靠,具有良好的人机对话功能。     

大地电磁测控深二维地形影响及其校正方法研究

采用有限元法实现ＭＴ二维模型的正演模拟,将地形条件下的二维地电模型的ＭＴ响应分解为地形引起的畸变场和与测点高程位置有关的稳定场两部进行分析。用有限元法模拟二维地形条件下地电结构的ＭＴ响应和地形的ＭＴ响应,采用幽会法消除地形引起的干扰场,然后,将消除地形影响后的视电阻率值通过延拓方法获水平地形条件下二维地电结构的视电阻率分布,完善在地电磁测深视电阻率的地形校正。     

带地形的可控源音频大地电磁法二维正演

利用二次场算法研究了可控源音频大地电磁法二维正演问题。采用有限单元法进行正演模拟,将矩阵压缩存储和共轭梯度解方程方法应用到正演算法中,加快了正演算法的速度,并且将地形因素考虑到正演算法中。通过不同的模型验算,检验了算法的精度。     

基于球体层状介质模型的大地电磁正演

目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。     

基于有限单元法的二维／三维大地电磁正演模拟策略

对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。     

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

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

大地电磁三维模型二维反演计算数值模拟分析

设计了多个1D/3D、2D/3D模型进行数值模拟分析,从模型三维响应和二维响应的差异大小和反演结果对比两个方面,分析了利用二维反演解释方法来逼近解释三维地质模型的适用性问题.分析结果表明:在异常体沿走向长度逐渐缩短的过程中,模型的三维特性逐渐增强,采用TM极化模式的二维反演结果基本保持不变,TE极化模式的二维反演结果逐渐产生畸变,表明TM极化模式对模型二维性的要求远远低于TE极化模式,TE/TM联合模式二维反演效果通常介于前两者之间,有时优于TM极化模式的二维反演结果.综合得出:一般条件下采用TM极化模式进行二维反演通常可以较好地重建原始模型,TE/TM联合二维反演结果在与TM极化模式反演结果整体电性特征相似的情况下,其对模型电阻率值的反演优于TM极化模式.