瞬变电磁正演之带激励源三维曲线交错网格时域有限差分方法.pdfVIP

·240· 中国地球物理2011 7 电磁方法研究与应用 瞬变电磁正演之带激励源的三维曲线 交错网格时域有限差分方法 A3Dcurvilinear finitedifferencetime—domainmethod staggered-grid v能ithall fortransient method source electromagnetic impulse 李展辉+黄清华 Li Zhanhui QinghuaHuang 北京大学地球与空间科学学院地球物理学系北京100871 瞬变电磁法 TEM 是一种在勘探与工程地球物理等领域得到广泛应用的人工源电磁探测方法。20 际TEM测量常在山区进行，地表是起伏不平的，源也不是偶极子，大多是矩形框。如要使正演更加真实， 地形因素以及源的形状也需要考虑。因此，Wang和Hohmann所阐述的方法会碰到这样两个问题： 1 向上延拓的方法不再适用，空气层需要考虑，时间步长会因此大大缩短，计算时间难以承受。 2 有地形情况下，不管是偶极子源还是矩形框发射源，初始值很难计算。 针对以上两个问题，尤其是第一个问题，一些学者进行了尝试【2J，但还是不能很好的解决。我们针对以 上两个问题，在使用曲线交错网格的情况下进行了探索研究。 1．引入介电常数可变化的空气层 上空气层并不需要进行低频近似。因此在对地下介质使用Dufort格式的同时，可以让空气中的介电常数随 地下介质介电常数共同变化。根据菲涅尔反射公式，在地下电磁场满足低频近似的情况下，空气中的介电 常数同地下介电常数一起变化对空气一地面边界的反射率和折射率的改变也是非常微小的。因此，这种变化 对地下电磁场的影响也是非常微小的。在此基础上，空气层不会再对Dufort格式的时间步长造成大的影响， 可以置于计算空间内直接计算。 2．引入脉冲激励源替代初始值场 初始值场在带地形的情况下很难给出，需要其他形式的源。加入空气层之后，便可以直接使用脉冲激励 源。TEM正演理想源是6函数，但是在有限差分中是不可能实现的。理论分析和正演中，所考虑的只是低频 近似的情况，那么6函数中不满足低频近似的高频部分可以不用考虑。这表明可以使用其他脉冲函数代替6 函数，高频部分可能不完整但是低频部分和6函数几乎一致，比如高斯函数。TEM野外测量普遍采用矩形 线框，即使在平坦地表，这种源的初始值计算也是比较耗时的。而使用有限差分对TEM进行正演很容易构 造矩形线框激励源，不再需要初始值的计算。因此我们引入矩形线框产生高斯脉冲激励源替代初始值场。这 种源不像初始值场一样严重受地形的限制，对地形有比较强的适应性，更加适合在有起伏地表的情况下使用。 3．方法验证 方法验证是数值模拟中很重要的一方面。我们采用两种验证方式： 1 半空间情况下数值结果和理论值比较。 2 起伏地表情况下，因为没有理论解，我们首先令空气中和地下介质中的介电常数不随时间变化， 时间步长也不变化，然后以此结果作为参考值和我们的方法进行对比。 两种验证结果都呈现了很好的一致性。 参考文献 HohmarmG彤A solutionfor3-dimensional [11 Wangr finite-difference，time·domain electromagnetic 1993，58：797-809． Endo K．Three—dimensional the forthecaseofthetime—domain [21 M，Noguchi modelingconsideringtopography electromagnetic Three-dimensional ofthesecondinternational in method,in Geochemistry eleetromagnetics：proceedingssymposium[J]．Methods and Science，2002，35：85～107． Geophysics，Elsevier