考虑关断时间的回线源激发TEM三维时域有限差分正演

被引:120
作者
孙怀凤
李貅
李术才
戚志鹏
王祎鹏
苏茂鑫
薛翊国
刘斌
机构
[1] 山东大学岩土与结构工程研究中心
[2] 长安大学地质工程与测绘学院
基金
国家自然科学基金重点项目;
关键词
瞬变电磁; 三维正演; FDTD; 关断时间;
D O I
暂无
中图分类号
P631.325 [];
学科分类号
摘要
从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.
引用
收藏
页码:1049 / 1064
页数:16
相关论文
共 23 条
  • [11] 瞬变电磁法三维正演研究[D]. 陈丹丹.中国地质大学(北京). 2008
  • [12] 带地形的瞬变电磁法三维数值模拟[D]. 肖怀宇.中国地质大学(北京). 2006
  • [13] 并行时域有限差分[M]. 中国传媒大学出版社 , 余文华[等著], 2005
  • [14] 时间域电磁法原理[M]. 中南工业大学出版社 , 牛之琏编著, 1992
  • [15] 电磁测深法原理[M]. 地质出版社 , 朴化荣编著, 1990
  • [16] Inversion of time domain three-dimensional electromagnetic data
    Haber, Eldad
    Oldenburg, Douglas W.
    Shekhtman, R.
    [J]. GEOPHYSICAL JOURNAL INTERNATIONAL, 2007, 171 (02) : 550 - 564
  • [17] Three-dimensional inversion of airborne time-domain electromagnetic data with applications to a porphyry deposit .2 Yang D K,Oldenburg D W. Geophysics . 2012
  • [18] 3Dinversion of airborne electromagnetic data using a moving footprint .2 Cox L H,Wilson G A,Zhdanov M S. Exploration Geophysics . 2010
  • [19] 瞬变电磁法中两种关断电流对响应函数的影响及其应对策略
    白登海
    Maxwell Meju
    [J]. 地震地质, 2001, (02) : 245 - 251
  • [20] 3D瞬变电磁场的有限差分正演计算
    宋维琪
    仝兆歧
    [J]. 石油地球物理勘探, 2000, (06) : 751 - 756+820