【摘 要】
：
海洋可控源电磁法数值模拟策略分为背景/异常场和总场算法.其中,总场算法需要通过靠近源的网格细化方法减少源的奇异性的影响,会显著增加电磁问题未知数的数量和解决时间.因此,如何通过增加适当的网格来适应高度复杂的地下电性结构是提高计算效率的关键.本文实现了一种三维复杂结构的面向目标自适应海洋可控源电磁有限元方法。通过自适应网格细化来提高有限元解决方案的准确性，对于大规模三维MCSEM数据，使用了METI
【机 构】
：
河南理工大学物理与电子信息学院,焦作、河南、中国;中石油东方地球物理公司综合物化探处,涿州、河北、中国;成都理工大学“地球勘探与信息技术”教育部重点实验室,成都、四川、中国
【出 处】
：
第十三届中国国际地球电磁学术讨论会
论文部分内容阅读
海洋可控源电磁法数值模拟策略分为背景/异常场和总场算法.其中,总场算法需要通过靠近源的网格细化方法减少源的奇异性的影响,会显著增加电磁问题未知数的数量和解决时间.因此,如何通过增加适当的网格来适应高度复杂的地下电性结构是提高计算效率的关键.本文实现了一种三维复杂结构的面向目标自适应海洋可控源电磁有限元方法。通过自适应网格细化来提高有限元解决方案的准确性，对于大规模三维MCSEM数据，使用了METIS来进行网格计算任务量的划分，形成的大型复数线性方程组分解为等价的实数形式，并利用Auxiliary Space Maxwell(AMS)预条件的共轭梯度法求解，使并行计算达到均衡负载的效果，通过对比一维解析解与三维自适应MCSEM计算结果，验证了程序的正确性。三维复杂模型在不同迭代条件下的网格划分和误差估计子的计算，证明了程序的鲁棒性和可扩展性。
其他文献
大地电磁测深是研究深部电性结构各向异性的主要方法,完善的各向异性正反演算法及参数灵敏度计算方法对于提高其解释可靠性是不可或缺的.模型参数灵敏度的计算分析对评价地球物理观测的地质构造分辨能力和实现地球物理线性化反演算法都极为重要,因而很有必要研究实用、高效而又准确的灵敏度计算方法.尽管已有相关研究报道,如Pek等(2003)展示了对离散方程求导,采用直接法或互易原理计算二维模型参数灵敏度的过程;利用
目前勘探的很多区域,地震资料对中浅层有很高的分辨率,但在复杂区对深层结构揭示不清,而电磁资料可以为地震提供重要的补充.同时,利用已知的地震、地质、电测井等资料进行电法约束反演是提高电磁勘探分辨率的重要途径.SEG年会也为此开辟了多种数据约束与联合反演的Workshop,并且在其他很多专题报告中,约束联合反演已成为一项提高分辨率的重要技术.本文在Hu et a1．(2015)工作的基础上，进一步发展
经过几十年的发展,海洋可控源电磁法(MCSEM)逐渐被纳入了海洋油气勘探的基本流程,并能极大的提高海洋勘探作业的油气储层探测能力,降低勘探成本.海水能够屏蔽空间电磁噪声,有利于应用MCSEM探测海底高阻异常体,但是MCSEM信号依然会受到随机噪声、空气波及海水扰动噪声的影响,导致信号失真,直接影响电磁数据的定性显示分析及定量反演解释效果.因此对MCSEM数据进行去噪方法研究,是后期数据处理和资料解
时域地空电磁信号易受到天电噪声、飞行载体发动机电磁噪声的影响,导致实测数据的信噪比降低,成像结果出现假异常,干扰对大地电性结构的准确识别.因此,研究时域地空电磁噪声抑制方法,提高电磁数据成像质量是必要的.本文基于平稳小波变换对时域地空电磁数据进行分解，并采用衰减斜率法识别细节分量中的电磁噪声，对相应采样点上的细节分量处理后，最终通过平稳小波逆变换实现电磁噪声抑制的目的。
在大地电磁观测中,由于天然电磁源的复杂性,数据误差之间存在着频率相关性和空间相关性.研究表明,在固定维数的大地电磁反演中,忽略误差的频率和空间相关性会导致模型不确定度的低估,并且会使参数估计产生偏差.但是,在跨维大地电磁反演中,数据误差相关性会对模型参数化的选择和模型参数值及其不确定估计产生什么影响尚未可知,本文采用高效的跨维贝叶斯反演方法，研究了数据误差的相关性对地电构造反演影响。结果表明忽略相
大地电磁阻抗张量畸变校正是大地电磁资料处理的重要内容,目前阻抗张量分解的方法主要有GB分解(Groom&Bailey,1989)及其衍生算法(多测点多频点GB分解(Joneset al,2001);分解参数在有噪声存在的情况下不稳定，尤其是电性主轴方向，而电性主轴方向是反映区域构造信息的重要参数，提出了一种基于相位张量分解估算稳定电性主轴方向的方法，并通过二次方程对视电阻率曲线畸变进行校正。本文对
逆时偏移一直以来是探地雷达(Ground Penetrating Radar,GPR)理论和应用研究的重点,通过逆时偏移成像可获得地下复杂结构精细分布.与基于射线理论的波动方程积分偏移和波动理论的微分波动方程单程波偏移相比,基于全波方程的逆时偏移具有精度高、相位准确和不受横向变速和高陡倾角的影响,能够精确的处理沿任意方向传播的波场的优点,可更精准地将复杂电磁波场归位到其真实的地下空间位置.GPR逆
空中电场电磁法(Air electric field electromagnetic method,AEMT),使地表的交变电场在空气中观测成为可能.本文根据这种方法将EH-4的一对接地电极换成电容进行测量,其它装置保持不变,并与传统接地方法在野外进行实验比对.用WINGLINK对其所测数据进行二维反演,其结果有效反映了地下的地质结构,验证了这种方法的可行性.
大地电磁(MT)反演方法最主要的问题是解的非唯一性,如何解决不适定问题是一项巨大的挑战.传统反演方法中,解决此问题的绝大多数方法都是基于Tikhonov正则化理论,最光滑约束是最常用的一种方法，基于模糊聚类算法，应用先验信息的约束反演方法能够对得到较好的地电模型，然而，此方法要依赖于先验信息。如果先验信息未知或者不准确，得到的反演结果并不理想。本文所提出的基于模糊类隶属度控制的MT反演方法，能够得
地震的发生受构造运动控制,构造运动会改变局部岩石圈的物理特性,导致岩石圈磁场发生变化.川滇菱形地块是南北地震带南段的重要地质构造单元,汇集了丰富的断裂构造分布.在该区域开展地磁场时空变化观测,研究地下磁性结构时空分布,可以为研究地震孕育过程中地下介质变化过程提供磁性结构分析依据,为川滇菱形地块北端的地壳物质动力学过程提供磁学印证.