简介：传统上,时间域航空电磁数据通过拟合迭代反演计算得到大地模型,然而,由于航空电磁数据道间的较强相关性,导致病态反演,并引起超定问题;同时电磁数据的相关性使其与模型参数的映射关系复杂,增加了反演的复杂度。采用主成分分析法将航空电磁数据变换为正交的较少数量的主成分,不仅降低了数据道间的相关性,减小了数据量,同时压制了数据的不相关噪声。本文利用人工神经网络(ANN)逼近主成分与大地模型参数间的映射关系,避免了传统反演算法中雅克比矩阵的复杂计算。层状模型的主成分神经网络与数据神经网络的反演结果对比显示,主成分神经网络反演方法网络结构简单,训练步数少,反演结果好,特别是对于含噪数据。准二维模型的主成分ANN、数据ANN以及Zhody方法的反演结果显示了主成分神经网络具有更接近真实模型的反演效果,进一步证明了主成分神经网络反演方法适合海量航空电磁探测数据反演。

简介：本文在Oristaglio等(1984)和Adhidjaja等(1985)工作基础上,给出线源二维时间域瞬变电磁二次场的DuFort-Frankel有限差分数值解,有效避免了在总场求解法中场源附近的奇异问题,并对地-空边界电导率的处理、归一化感应电动势偏导数的计算、推进时间步的确定,提出了改进方法;吸取前人成就中二次场地-空边界向上延拓和零值边界处理技术,从而简化了计算方法;通过对均匀大地、水平层状大地模型的计算,二次场求解法与解析法的最大相对误差小于0.01%,计算速度比总场求解法提高了约3倍;模拟计算不同时刻瞬变电磁场在地下的分布形态,描绘出感应涡流向下向外的传播特征,以及与地下异常体相互作用的物理过程。

简介：波场重构反演是一种改进的全波形反演理论。该反演方法通过将波动方程引入目标函数中拓宽了解的寻找空间，通过重构真实波场来计算模型梯度，大大提高了计算效率的同时还减弱了局部极小值的影响。但目前该理论基本在频率域进行，而频率域反演对计算内存的需求太高，并且很难应用到实际生产中。因此，本文将波场重构反演拓展到时间域，推导了时间域波场重构的增广方程，结合模型试算结果对波场重构的模型梯度进行了修改。数值实验表明，时间域波场重构反演准确性较高并且对低频信息具有良好的重建能力。

简介：本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.

简介：无网格法形函数构造不依赖预定义的单元，具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法（EFGM）、点插值法（PIM）与径向基点插值法（RPIM）三种全域弱式无网格法的近似原理及特点；以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演，给出了RPIM形状参数的推荐值；分析了均匀介质模型大地电磁（MT）二维正演无网格法边界条件直接加载与罚函数法加载的精度差异，结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势，提出了EFG—PIM及EFG．RPIM耦合算法，数值计算结果验证了耦合算法的有效性。研究发现：无网格法及其耦合方法适用于电磁法数值模拟；支持域无量纲尺寸取1．0时无网格法精度与效率高，场节点与背景网格重合时计算效果佳；泊松方程求解PIM及RPIM精度较EFGM低，计算均匀介质MT响应精度较EFGM高；RPIM改善了PIM计算涉及的奇异性问题，对应支持域无量纲尺寸选择空间大。

简介：海洋天然气水合物是一种潜在的巨大能源,而地球物理探测技术是勘探天然气水合物资源的重要方法。本文介绍了利用可控源电磁法探测海洋天然气水合物的可行性,在参考了大洋钻探计划(ODP)164航次对天然气水合物的采样数据的基础上建立了几种不同水合物含量的海底天然气水合物一维地电模型,并且利用建立的模型讨论了天然气水合物的频率域电磁响应。文章就天然气水合物的电场振幅值和相位随着收发距和频率等参数的变化关系等做了相应研究,同时就水合物厚度不同时其相应的电磁响应特征做了初步探讨,为海洋天然气水合物的勘探和资源评价提供了参考。

简介：频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有�

简介：本文是中国陆域航磁地质构造研究最新进展的综述性文章，对基于航磁数据的中国陆域地质构造编图最新成果和取得的新认识进行了论述。该项研究利用我国50多年的海量实测航空磁测数据，采用最新的数字化编图技术及先进实用的位场数据处理方法编制了1：5000000万和1：1000000万中国陆域航磁系列图件，以板块构造和大陆动力学理论为指导，以岩石物性研究与航磁异常特征分析为基础，定性、定量与综合解释相结合，编制了中国陆域磁性基底深度图及主要地质构造单元图件，并依据航磁数据对新划分的地质构造进行了论述。

简介：本文介绍了一种改进的角度域VSP逆时偏移方法。对VSP逆时偏移中的逆推公式进行了改进,为方便数值计算出相对保幅的角度域共成像点道集（ADCIGs）。此外VSP记录到的波场信息丰富,包括上行波场、下行波场和直达波场等,本文分析了这些波场的响应特征,发现直达波和下行波在角度域共成像点道集（ADCIGs）上都产生了成像噪音,直达波产生的噪音尤为严重。把该方法用于我国西部地区实际观测的VSP资料,不仅获得相对保幅角度域共成像点道集（ADCIGs）,而且压制了成像噪音。通过数值模型试算,实际资料的应用验证了该方法的实用性与有效性,从而为VSP偏移速度分析、VSPAVA/AVO分析和反演等提供可靠的基础资料。

简介：传统叠前深度偏移只能够提供地下的构造信息，但工业界在需要构造信息的同时还要与地下界面反射系数成比例的振幅信息。最近几年，基于单程波方程的保幅叠前深度偏移算法有了一定的发展，但是，基于炮域、单程波的保幅型叠前深度偏移必须应用反褶积型的成像条件，这种成像条件在构造复杂、速度变化剧烈的地区会出现不稳定现象。基于角度域的保幅深度偏移克服了这一不稳定性缺点的同时，还域的保幅深度偏移，模型和实际资料的试算分析验证该思路方法的正确性和有效性。

简介：F-X域经验模态分解去噪方法在处理非稳态地震数据时存在两个局限,一是单纯剔除第一个固有模态分量将导致有效信号缺失及去噪能力偏弱问题,二是分解复信号时对实部和虚部分别分解存在分解数目不一致的风险。本文对上述两个方面进行了改进,提出了一种新的F-X域投影法复数经验模态分解预测滤波方法,首先采用基于空间投影的复数经验模态分解将F-X域地震数据直接分解为不同的复固有模态分量,然后再对这些分量分别进行F-X域预测滤波。合成记录及实际资料测试表明,本文的新方法能更好地衰减随机噪声,更有效地保持地震信号。

简介：高频电磁测深法是一种利用天然高频电磁场作为场源的电磁勘探方法,具有高分辨率和相对直流电法较深的勘探深度的特点,特别适合中浅层地热和地下水勘探。本文介绍了一个应用高频电磁测深评价地下热水资源的一个成功实例。野外数据采集采用了高频电磁测深系统MT-U5A,采用远参考技术,以保证数据采集质量。高频电磁测深资料二维反演成像处理结果较好地描述了地下热水的空间分布范围,为查明该区的地下热水的来源提供了较有利的证据。

简介：在对浙江省某铜矿进行瞬变电磁勘查时，由于感应激发极化效应的影响，造成瞬变电磁晚期测道的数据发生倒转。采用常规的瞬变电磁数据处理方法难以处理，造成晚期测道数据不可用，影响TEM的探测深度及精度。针对这个情况，采用Cole—Cole模型分析了均匀半空间模型中直流电阻率、充电率、时间常数以及频率相关系数对磁性源瞬变电磁响应的影响规律。利用奇异值分解法（TheSingularValueDecomposition，简称SVD）对实测TEM数据进行分析及反演，并且从瞬变电磁响应数据中分离出勘探区瞬变电磁测量数据Cole-Cole模型的各参数，将分离出来的参数用于探测结果的辅助解释。经过分析得出充电率和频率相关系数对瞬变电磁响应的影响较大，而直流电阻率和时间常数对其影响较小的结论。通过常规方法和奇异值分解法对实测数据分别进行处理，发现后者视电阻率断面图的异常更加突出，从而较准确的圈定了铜矿体地范围，此结论得到了钻孔的验证，与实际地质情况吻合较好。

简介：用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深（LOTEM）的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。

简介：电磁全息测量数据包含两种模态，因此重建流动图像也需采用“双模”融合的敏感场。首先，通过电磁全息探测物理场分析，结合层析成像数学理论Radon反变换，从定解问题推导了全息测量敏感场函数表达式。其次，将有限元仿真计算得到的全息测量敏感场应用于模拟流动试验和全息成像。结果表明，基于复电位φ关于极径r的偏导数的全息测量敏感场契合了Radon反变换的数学表达，且充分体现了幅度、相位测量敏感性；反演所得流动图像更加吻合实际流型。全息测量敏感场构建有效克服了传统单模敏感场在计算精度和计算效率等方面的局限性。

简介：频率域全波形反演充分利用全波场的振幅、相位以及频率信息，采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础，对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中，针对截断边界问题，并考虑多频率联合反演中计算区域采用同一套剖分网格的需求，提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子，通过一定方式调节衰减因子使得声波在衰减层中充分衰减，达到压制截断边界影响的目的。根据指数衰减规律，文中推导出了一种新的衰减因子计算公式，并给出了不同频率条件下衰减层厚度计算公式；在反演算法中，采用共轭梯度法求解高斯牛顿反演迭代方程组，避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量，并采用相同的反演模型，对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明，本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。

简介：应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。

简介：为了进一步提高大地电磁三维快速松弛反演的计算效率，在深入分析大地电磁三维快速松弛反演算法的基础上，结合MPI自身的优越性，确定了并行计算的思路，实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算，分析对比了在多种情况下程序的执行效率。测试结果表明，所实现的三维快速松弛反演并行程序运行结果正确，效率提高明显。此成果在普通微机上实现，推动了大地电磁三维反演技术的实用化，可为其它地球物理三维正反演研究所借鉴。

简介：研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解，由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件，有效地消除了低波数干扰，同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析，在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像，利用校正因子对各角度的像进行校正，进而实现对成像结果的校正，从而实现了一种高效的倾角域幅度校正方法。最后通过SEG／EAGE模型进行数值计算验证了文中所述的计算方法。

简介：全波形反演是一种高精度的地震成像方法，可以对地下介质物性参数模型进行准确的重构。然而在实际应用中，尤其是在三维复杂介质反演中，计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本，提高反演效率。但是使用震源编码技术也带来了两个问题：一方面，参与编码的各个震源之间会产生“串扰噪声”，导致反演结果中出现假象；另一方面，基于震源编码的频率域全波形反演方法周围噪声较为敏感，使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”，并基于震源编码技术提出一种频率域自适应全波形反演方法，通过一个与频率相关的自适应选择机制，将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来，在保证反演质量的同时也最大程度地提高了反演效率。