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

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

二维大地电磁正演模拟综述

<正>随着经济的发展,有色金属作为一种重要的物质基础,是关系到一个国家的战略性资源;随着浅层地面的资源日益开采,目前提出了深部找矿的概念,而大地电磁由于其勘探深度大,目前在深部找矿中的应用越来越广泛。本文综述了二维大地电磁正演模拟理论,其方法主要有有限单     

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

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

海洋大地电磁二维正演及结果分析

研究了海洋大地电磁法的二维正演并编制了相应程序 ,程序采用有限单元法求出海底各节点处的电磁场值 ,进一步计算得到视电阻率。为了验证程序的正确性与有效性 ,根据卡尼亚理论推导出计算一维海洋模型的阻抗公式 ,并得到视电阻率的解析解 ,将其与有限元法得出的数值解进行了对比 ,表明了此方法的正确性。通过对正演结果的分析 ,表明海洋大地电磁测深能够成功地反映海底的电性结构     

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

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

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

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

基于ILU-BICGSTAB算法的起伏地形MT二维正演模拟

针对起伏地形对异常体的探测效果的影响,利用有限单元法,基于MTLAB平台编写了大地电磁二维正演程序,为更好地模拟实际地形,采用矩形网格中进一步三角细化的剖分方式,并利用不完全LU分解的预条件BICGSTAB算法提高计算速度与精度,最后通过建立多种地电模型分析讨论了起伏地形对大地电磁场的影响,验证了该方法的实用性与准确性。     

倾斜电各向异性介质大地电磁正演研究

研究地层电各向异性性质能够圈定沉积地层中的油气储层和与地球动力学过程、地震预报有关的深部线性构造。通过得到的N层倾斜电各向异性介质大地电磁响应公式,在不同的电各向异性系数、地层倾角、倾斜层厚度参数下,对电各向异性倾斜地层模型进行了正演计算。结果表明:当电各向异性系数、地层倾角、地层厚度三者中任一变量改变时,视电阻率曲线的形态基本保持不变,但视电阻率曲线振幅值变化较大;视电阻率曲线对地层倾角参数的变化反应最为敏感;电各向异性系数的增加,能突显薄层电阻率异常。这有助于以后的电各向异性介质大地电磁反演和应用研究。     

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

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

大地电磁(MT)自适应有限元各向异性正演

电各向异性在自然界中普遍存在,特别是沉积盆地中的部分岩层经过压实变质作用,表现出很强的各向异性导电性,采用基于各向同性模型对实测资料进行正反演解释必然会造成困难甚至结果的错误,只有基于各向异性理论的正反演解释才更为合理准确。本文通过自适应有限元对大地电磁各向异性进行正演计算,分别模拟水平、垂直和倾斜各向异性介质在不同偏转角和主轴电阻率下的响应结果。结果表明：自适应有限元能够在后验误差的控制下得到合理的网格,使计算结果更加接近解析解;在各向异性介质中,大地电磁TE极化模式的视电阻率和阻抗相位与垂直于层面的电阻率无关;二维电各向异性结构中,大地电磁TM极化模式响应结果总是由主轴上的电阻率在y轴方向上的分量所决定。     

水平层状电各向异性介质大地电磁正演研究

电各向异性现象在地球内部是普遍存在的,通过理论上推导的电各向异性介质的张量电导率的表达式及水平层状电各向异性介质的大地电磁正演公式,对一些地电模型进行数值计算,同时与电各向同性地电模型进行对比分析。结果表明,电各向异性介质的张量电导率的大小不仅与介质的固有电阻率性质有关,还与测量的方向有联系;电各向异性系数和相对测量角的变化都对地层的视电阻率和相位曲线振幅有较大的影响,但其曲线形态基本保持不变;当野外观测方向和电各向异性介质的电性主轴重合时,电各向异性介质模型可以看成两个独立的电各向同性介质模型。这些结果有助于以后的电各向异性介质实验研究和其大地电磁反演研究。     

起伏地形下地面瞬变电磁法三维正演数值模拟研究

瞬变电磁法在地形条件复杂的地区应用时,勘探效果往往不理想,究其原因,主要是忽略地形因素影响致使解释结果偏差较大.笔者采用三维时域有限差分法模拟起伏地形下瞬变电磁响应,通过麦克斯韦方程直接导出瞬变电磁扩散方程,引入虚拟位移电流项构建时域有限差分方程.首先,针对均匀半空间模型的计算结果与解析解进行对比,检验本文算法精度;然后,模拟分析了山峰、沟谷模型的地形响应特征,得出早期信号表现为两种特征,沟谷地形,早期“前段”较正常地形感应电动势值减弱,导致“假高阻异常”,早期“后段”较正常地形感应电动势值值增强,导致“假低阻异常”,山峰地形特征恰好相反,至晚期地形影响呈逐渐减弱的趋势;最后,通过实测数据证实正演模拟结果的正确性.研究结果对地形效应识别及校正工作具有较好的参考意义.     

考虑激电效应的二维大地电磁测深正演

通常情况下地下岩矿石在通电的情况下,不仅表现出电磁效应,还表现出激电效应,在激电效应强烈的地区两者难以分离。笔者通过对均匀半空间中两种极化体的正演模拟,分析Cole--Cole模型参数对极化体的影响规律。模型正演结果表明:激电效应的存在会使视电阻率值增大,不同激电参数对视电阻率产生的影响不同,其中极化率对视电阻率的影响最大;极化率越大,视电阻率越大。     

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

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

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

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

基于三次插值的大地电磁自适应有限元二维正演模拟

大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。     

大地电磁三维矢量有限元正演

采用广义变分原理,基于矢量基函数详细推导了大地电磁三维矢量有限元方程。为了提高计算精度和效率,应用直接法强加边界条件改善总体系数矩阵的条件数,同时使用SSOR(symmetric successive over relaxation)预处理的双共轭稳定梯度法求解复对称大型稀疏线性方程组。并利用国际标准模型与相关参考文献的结果进行了对比,验证了算法的准确性。对一个典型的三维低阻体模型进行正演,得到了不同测线的视电阻率和相位断面图,并与二维正演结果进行对比分析。结果表明:在x方向测线上,ρ_(yx)变化幅度较ρ_(xy)小,中心测线上的ρ_(yx)和ρ_(xy)响应均与二维TM模式条件下的响应特征相似。     

起伏地形下可控源音频大地电磁三维数值模拟

起伏地形对可控源音频大地电磁(CSAMT)响应具有强烈的影响,因此在CSAMT数据处理解释时需要考虑地形。同时,实际的地下地质情况和地表的地形情况通常比较复杂,地质结构和地形大部分情况下都是三维的。在水平地表三维有限差分CSAMT数值模拟算法的基础上,推导了利用地下交错网格采样点处的总磁场计算起伏地形下空气-地下介质分界面处的总电场和总磁场的表达式,从而实现了起伏地形下三维CSAMT数值模拟算法。在算法实现过程中,采用伪δ函数代替麦克斯韦方程中的场源项和直接计算总场的策略,避免了原有的将总场分离成背景场和二次场的策略在复杂地质条件下难以选择合适背景电阻率的问题。为了直接模拟总场,起伏地形下三维CSAMT数值模拟算法给出了新的三维正演方程的边界条件。将模拟水平地表三维异常体和三维山峰地形两个理论模型得到的响应结果与前人算法的计算结果进行对比,验证了所实现算法的正确性和有效性。     

起伏地形重磁三维快速正演计算

正演是反演技术的基础,正演速度和求解反演问题的系数矩阵存放一直是起伏地形下重、磁三维反演的关键技术问题。这里提出了一种起伏地形下重磁快速正演计算方法,其计算原理是根据反演在垂向的剖分层数,利用水平地形正演计算形成二个不同大小刻度标尺矩阵,然后在模型空间,使用分段线性插值的方式,直接计算出起伏地形观测点的正演值。该方法的主要特点是在保持很高的计算精度下计算速度可提高二倍,且节省计算内存,适合起伏地形下重磁三维反演技术研究。     

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

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