排序方式: 共有31条查询结果,搜索用时 796 毫秒
21.
采用规则六面体单元和并矢Green函数奇异积分等效积分技术,已有的大地电磁积分正演方法具有不能有效模拟地下复杂地质体和计算精度偏低的缺点.本文提出了一种新的三维大地电磁积分方程正演技术,即采用四面体单元、解析的并矢Green函数奇异积分表达式,达到既能模拟地下复杂异常体,又能有效提高已有积分方程法计算精度的目的.首先,采用四面体网格技术离散地下复杂异常体,获得四面体单元上的大地电磁积分方程.然后,利用针对四面体单元开发的新的奇异值积分的解析表达式,准确计算线性方程中的并矢Green函数的奇异积分,从而获得精确的线性方程.借助于PARDISO高性能并行直接求解器,实现了三维大地电磁问题的高精度求解.最后,基于国际标准3D-1模型和六棱柱模型,通过与其他方法结果的对比分析,验证了本文方法的正确性、处理高电导率对比度的能力(1000:1)和处理复杂模型的能力. 相似文献
22.
音频大地电磁(Audio Magnetotelluric,AMT)信号常常受到持续性人文噪声影响,这类噪声使用远参考法和Robust阻抗估计等常规方法往往效果不佳.为此,本文从噪声的规律与特征出发,提出一种新的AMT数据处理方法.首先通过字典学习方法从观测数据中自主学习到人文噪声的特征结构,构建冗余字典,然后利用学习到的冗余字典,分离出AMT数据中的人文噪声.为验证方法的有效性,首先进行了合成数据的仿真试验,然后在四川凉山进行了针对性的野外试验研究,最后将本文方法应用于庐枞矿集区实测数据的处理.结果表明,本文方法能够快速、准确地分离出AMT信号中的人文干扰,保留有用信号,显著改善AMT数据质量. 相似文献
23.
提出了一种基于非结构化四面体以及带地形模型的自适应多层快速多极大规模磁法快速正演算法.该算法弥补了传统积分方法采用FFT加速计算时不能采用非结构化网格的缺陷;同时采用自适应快速多极算法突破积分求和方法求解大规模磁法问题耗时长的突出问题.首先,采用非结构化的四面体网格剖分技术能够更好的模拟复杂模型以及带地形模型,实现磁法模型的高精度模拟;其次,采用一种自适应多层快速多极(AMFM)算法实现大规模磁法正演求解.通过将计算区域划分为近区和远区,对近区采用解析计算高精度求解,对远区采用自适应多层快速多极算法进行加速计算,假设有M个观测点,N个四面体源单元,可将计算复杂度由传统积分求和法的O(MN)减少到O(Mlog N).本文设计了组合体模型以及安徽怀宁地区的实际地形模型,模型计算结果体现了采用该方法进行大规模复杂模型三维磁法正演模拟的高效性和准确性. 相似文献
24.
重力梯度张量曲率目前广泛用于重力数据的处理和解释中.为了拓宽重力梯度张量曲率的应用,本文回顾了重力梯度张量曲率的定义,从等位面的曲率定义出发,讨论了正确计算曲率的测量参考系及局部旋转的相关理论,并以单个球体和棱柱体为例来说明曲率的正确计算方式.然后,在正确计算重力梯度张量曲率的基础上,将重力梯度张量曲率应用到重力数据的边界识别中,通过理论模型和实际数据详细分析和比较了各种曲率在重力边界识别中的应用效果.结果表明:基于等位面的局部旋转坐标系是各种曲率正确计算的先决条件,纠正了曲率计算的误区;在边界识别中,局部坐标系下所计算的高斯曲率进行边界识别能够较好的圈定地下地质体的边界. 相似文献
26.
传统的大地电磁法及音频大地电磁法都假设平面电磁波垂直入射地球表面,且忽略位移电流的影响.实际中,电磁波并不总是垂直入射地球表面,在某些情况下,位移电流也不可忽略.本文以新的方式推导了全电流条件下平面电磁波倾斜入射层状介质表面时的波阻抗公式,分析了位移电流、入射角、磁导率以及介电常数对TE和TM模式视电阻率之间相对差别的影响.结果表明:电磁波垂直入射地球表面时,TE和TM模式视电阻率相等,该特性与位移电流大小、磁导率和介电常数变化无关.当位移电流可忽略时,TE和TM模式视电阻率相等,该特性与入射角、磁导率和介电常数变化无关.当电磁波倾斜入射且位移电流不可忽略时,TE和TM模式视电阻率不相等,两者的差别随位移电流和入射角的增大而增大且与磁导率和介电常数有关.一般地,当位移电流超过传导电流的5%时,其影响就不可忽略.因此,实际数据处理及反演解释中,对特殊地质条件或频率很高时,应考虑位移电流、入射角、磁导率以及介电常数对TE和TM两种视电阻率模式之间差异的影响,否则可能得不到正确的结论. 相似文献
27.
射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性. 相似文献
28.
大地电磁信号具有非线性、非平稳、非最小相位特征,不符合以Fourier变换为基础的传统功率谱估计的基本要求. Hilbert-Huang变换是近年发展起来的处理非线性、非平稳信号的完全局部时频分析方法. 本文在简要介绍Hilbert-Huang变换基本原理与算法基础上,以实际数据分析为例,探讨了它在大地电磁信号处理及噪声压制中的应用. 提出利用Hilbert时-频能量谱对大地电磁信号进行时段筛选,以提高信号品质,增强数据处理的质量和资料的可解释性. 利用经验模态分解方法及其多尺度滤波特征,可以有效地分析MT信号中的噪声分布特征,并进行干扰压制. 相似文献
29.
边界识别是位场数据处理解释中的重要环节,传统边界识别方法通常不能均衡深、浅部地质体边界.基于尺度空间技术和归一化的Facet模型检测算子,本文开发了一种带通空间滤波和边缘检测相结合的边界识别方法,有效地提高位场数据边界识别的精度和可靠性.为了验证本文算法的有效性和稳定性,分析了不同尺度空间函数和检测算子对算法的影响,并且对比了传统边界识别方法的效果.理论模拟和实际数据分析表明,利用位场垂向二阶导数进行的基于尺度空间技术的归一化Facet模型边界识别方法不仅算法的稳定性强,而且可以避免高阶导数对噪声干扰放大作用,同时均衡深部和浅部地质体边界,从而可以更精确地识别地质体的形态. 相似文献
30.
为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法. 相似文献