三维时间域航空电磁有理Krylov正演研究

作者： 邱长凯   殷长春   刘云鹤   张博   任秀艳   齐彦福   蔡晶   来源： 地球物理学报 年份： 2020 文献类型 ： 期刊 关键词： 时间域   有理Krylov方法   航空电磁   有限元法   三维正演  

描述： 常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.

三维时间域航空电磁有理Krylov正演研究

作者： 邱长凯   殷长春   刘云鹤   张博   任秀艳   齐彦福   蔡晶   来源： 地球物理学报 年份： 2020 文献类型 ： 期刊 关键词： 时间域   有理Krylov方法   航空电磁   有限元法   三维正演  

描述： 常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.

航空瞬变电磁数据背景场去除研究

作者： 黄威   贲放   李军峰   殷长春   胥值礼   刘俊杰   来源： 物探与化探 年份： 2020 文献类型 ： 期刊 关键词： 数据预处理   线性插值法   航空瞬变电磁法   背景场去除  

描述： 与地面瞬变电磁法相比,航空瞬变电磁法受不同的噪声类型影响较大,如天电噪声、运动噪声等,并且背景场对数据的影响也是不可忽略的。由于背景场中含有一次场信息,使得背景场幅值远大于含有地质信息的二次场信号。为较清晰地获得地下电性结构,需找到可以高效去除背景场的方法。本文提出利用线性插值方法来对背景场问题进行研究。首先,获得前、后高空飞行的测线数据,经过叠加、平均计算得到叠加周期内相同采样时间道的响应值,然后,利用插值原理针对每个采样时间道建立一个插值函数,最后,根据飞行时间来计算工区每条测线上采集数据的背景场值,从而达到去除背景场的目的。实测数据背景场去除的算例表明,该方法可较好地去除背景场对二次场异常的掩盖,去除后可较清晰地判断异常出现的位置。时间域航空电磁数据预处理包含较多处理模块,该方法的提出为高效的数据处理奠定了基础。

基于谱元法三维航空电磁电各向异性模拟及识别研究（英文）

作者： 黄鑫   殷长春   曹晓月   刘云鹤   张博   蔡晶   来源： Applied Geophysics 年份： 2018 文献类型 ： 期刊 关键词： 正演模拟   频率域航空电磁   各向异性   谱元法   GLL插值基函数  

描述： 航空电磁法具有高采样率探测特征,传统数值模拟方法为保证航电响应计算精度,特别是针对复杂地质体(如电各向异性),需要对模拟域进行精细剖分,导致计算工作量庞大。因此本文采用谱元法对航空电磁各向异性响应进行三维正演模拟,该方法结合谱方法和有限元法双重优势,它与谱方法类似采样高阶基函数代替有限元中的线性插值,同时具有有限元拟合边界的灵活性。因此该方法能改善离散网格内部的数值模拟精度,减弱数值模拟结果对网格的依赖性,实现航空电磁各向异性响应的高精度计算。首先,本文将旋转张量的各向异性电导率引入Maxwell方程,并给出GLL谱元基函数给出电场表示形式,采用伽辽金加权余量法形成基于谱元法大型线性方程组获得电磁响应,以VMD发射源为例进行航空电磁响应模拟。其次,本文采用粗物理网格四阶谱元法对任意各向异性半空间的航空电磁响应进行计算,通过与一维半解析解对比验证了该算法的高精度性。再次,本文采用不同物理网格和谱元阶数讨论不同各向异性条件下三维异常体的航空电磁响应,得出了不同各向异性异常体的电磁响应的收敛条件。最后,本文首次基于谱元算法,采用Hr/Hz方式给出了三维各向异性异常体,各向异性围岩以及围岩和异常体同时存在各向异性时航空电磁各向异性的识别方式,这航空电磁各向异性的三维反演和各向异性地质区域测量具有指导意义。

时间域航空电磁数据变维数贝叶斯反演

作者： 高宗慧   殷长春   齐彦福   张博   任秀艳   卢永超   来源： Applied Geophysics 年份： 2018 文献类型 ： 期刊 关键词： 贝叶斯反演   时间域航空电磁   加权系数   反褶积   变维数  

描述： 为改善反演效果,获得全局最小解,减小反演结果对初始模型的依赖程度,本文将变维数贝叶斯反演应用于时间域航空电磁数据反演。变维数贝叶斯反演方法在贝叶斯方法基础上利用可逆跳跃马尔科夫链蒙特卡洛方法(RJMCMC)实现反演模型层数的变化。这种方法根据建议分布并利用蒙特卡洛方法充分搜索模型空间进行随机采样。只统计同时满足数据拟合要求和接受概率的候选模型,受初始模型影响小,收敛稳定,反演结果可靠,最终可获得反演模型的概率分布和不确定度信息。由于实际飞行中发射源高度很难精确测量,因此本文在反演过程中将发射源高度分为不变和发射源高度变化两种情况。同时本文在电阻率先验概率密度函数中引入加权系数以调整对反演模型的约束强度,可有效地解决电阻率断面中间层反演效果不理想的问题。本文通过反演中心回线装置的H型和分离装置K型、HK型断面添加高斯噪声后的仿真数据以及实测数据,验证了变维数贝叶斯方法反演时间域航空电磁数据的有效性。

时间域航空电磁各向异性大地三维自适应有限元正演研究

作者： 齐彦福   李貅   殷长春   戚志鹏   周建美   孙乃泉   刘云鹤   张博   来源： 地球物理学报 年份： 2021 文献类型 ： 期刊 关键词： 时间域航空电磁   异常识别   各向异性   自适应有限元   三维正演  

描述： 各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.

时间域航空电磁各向异性大地三维自适应有限元正演研究

作者： 齐彦福   李貅   殷长春   戚志鹏   周建美   孙乃泉   刘云鹤   张博   来源： 地球物理学报 年份： 2020 文献类型 ： 期刊 关键词： 时间域航空电磁   异常识别   各向异性   自适应有限元   三维正演  

描述： 各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.

三维时间域航空电磁任意各向异性正演模拟（英文）

作者： 黄威   贲放   殷长春   孟庆敏   李文杰   廖桂香   吴珊   西永在   来源： Applied Geophysics 年份： 2018 文献类型 ： 期刊 关键词： 任意各向异性   矢量有限元法   三维时间域航空电磁法  

描述： 在自然界中地层会呈现出各向异性特征,各向异性介质的研究有助于资料解释和认识地球动力学的过程,但当前数据处理中多数使用各向同性理论解释各向异性构造,因而对电磁资料解释造成了一定的偏差。当前,时间域航空电磁法关于各向异性的解释仍局限于一维情况,且相应的三维数值模拟也正处于研究中,因此本文针对任意各向异性介质中三维时间域航空电磁法进行正演模拟。将各向同性介质的电导率经过坐标旋转后,得到任意各向异性介质的电导率张量,将其引入到麦克斯韦方程组中,采用矢量有限元法的规则网格对求解区域进行剖分,利用直接求解器对总体合成的线性方程组进行求解。通过与解析解对比,验证了算法的精确性和可行性;同时与各向同性情况进行对比,通过响应的分布特征及幅值变化分析任意各向异性对三维时间域航空电磁响应的影响,得出在不同情况下各向异性对时间域航空电磁响应的影响情况存在较大差异。