简介：针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题，本文展开研究．解法方面，传统离散镜像法多采用复数运算、数字滤波等方法计算量大，针对这些问题，提出一种改进的离散镜像方法：基于Gaver-Stehfest概率变换算法将电磁场解式实数化，选用Prony方法对目标核函数进行指数级数逼近，根据离散镜像原理和近似系数闭合式求解瞬变电磁场．通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果，证明该方法有效且具有较好的精度和适用性．继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场，对“地面激发一地层中测量”方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关，感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况，在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。

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

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

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

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

简介：已有的均匀磁化长方体磁场计算方法大都假设观测点位于上半无源空间，对于起伏地形条件而言，这些方法可能存在解析“奇点”。为此，本文基于地磁场基本理论，采用变量替换的积分方法，导出了改进的磁场表达式，详细讨论并有效解决了整个无源空间的所有奇点问题。相比前人的方法，其积分过程更自然、简单，最后的积分结果形式更加统一，并且不需要坐标变换即可求出无源区任意点处的磁场值，从而简化了正演过程。对比模型试验表明，新导出的磁场无解析“奇点”理论表达式是正确的，并能适应地形起伏的情况。

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

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

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

简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：通过最新的全球地磁模型——CALS10K.1b,结合CALS3K.4与IGRF11模型,计算并分析了10000BC~1990AD期间中国大陆及邻近地区非偶极子（ND）磁场Z分量的时空变化。为了深入了解ND场的变化,从场源的角度,对2n（n=2-10）极子ND场及其对应的能量进行了分析。结果显示在研究期间ND场的变化可分为3个阶段。在10000BC-2500BC期间,ND场以正值为主并持续了近7500a,在2500BC-1500AD期间强度转弱为以负值为主并持续至1500AD,自此快速增强为以正值为主。东亚地区ND场异常基本在截断阶数（n）为3时即形成,且该异常区已在大陆地区内形成了封闭的圆形区域,这意味着前3阶的ND场占据了总ND场强度的大部分。ND场在核幔边界（CMB）处衰减较快,在地表处则趋于稳定。

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

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

简介：尾波是持续时间比主要波型长得多的多重散射波，它包含了直达波之外的部分有用信号。在油气田的开采过程中，由储层物性参数的微小变化而导致的储层速度的微小变化对初至的影响很小，无法通过初至变化直接观测，但是因为尾波为多重散射波，故储层速度的微小变化会在尾波的传播过程中被放大，因此利用尾波可以观测到这种变化。本文通过实验与数值模拟研究了尾波干涉方法在储层微小时移差异监测方面的有效性。合成地震数据是基于部分MarmousiII模型采用有限差分波动方程正演方法计算获得，为了模拟时移地震中的储层波速变化，在模型中选取波速发生微小变化的目标区，计算速度扰动前后的合成地震数据，进而观测尾波变化。研究表明利用尾波干涉可检测储层物性的微小变化，为时移地震信号监测提供新的方法和手段，从而提高时移监测的准确度，也为开发生产工作提供可靠的指导。

简介：在CSEM野外工作中,通常把接收机置于离源3～5倍趋肤深度以外的远区。本文从正演和反演两个角度研究了当源和接收机之间存在导电体时,该导电体对远区中目标体响应的影响。2D有限元正演结果表明导电体主要影响中低频的观测结果,导电体的电导率越低,尺寸越大,影响也越大,并且影响程度随着观测频率的降低而逐渐加大。反演结果表明若不考虑源和勘探区间导电体的存在,反演得到的目标体横向位置会向源的方向移动,带来解释上的误差。在CSEM实际工作中,应该是采用多方位的场源作三维采集和三维反演,以减少源和勘探区间低阻体的影响。更多还原

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

简介：阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期（超过100秒）,逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。

简介：提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明：本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。