逆Laplace变换新算法及其在时间域电磁响应计算中的应用

时间域电磁响应的正演计算多是由频率域响应经逆Laplace变换而得到.逆Laplace变换的计算精度和效率是时间域电磁响应计算中方法选择的重要指标.论文分析了几种逆Laplace变换的算法机制,并优选出Talbot算法计算了水平电偶源层状模型的时间域电磁响应.逆Laplace变换常用的算法有折线法、数字滤波算法和Gaver-Stehfest算法(简称G-S算法).折线法需要精细地确定分割步长以提高精度,数字滤波算法系数很多,适应频率范围受计算问题所限,而G-S算法受计算机字长和问题对象的影响大.本文在64位计算平台中计算比较了G-S算法、Euler算法和Talbot算法的节点数对于精度的影响,发现Talbot算法受节点数影响小,计算精度高,适应频率范围宽.最后利用21点Talbot算法计算了水平电偶源轴向偶极装置均匀大地模型径向电场的阶跃响应和冲激响应,计算精度及响应时间范围均优于G-S算法.计算了水平电偶源赤道偶极装置均匀大地模型垂直磁场的阶跃响应和冲激响应,冲激响应峰值时刻对于电阻率的变化响应灵敏,与轴向偶极径向电场响应能力相当,但垂直磁场随收发距增大,衰减较快.根据层状模型阶跃响应晚期渐近值计算的视电阻率,水平电偶源轴向偶极径向电场有能力发现大埋深高阻或低阻薄层,收发距应大于中间目标层埋深的5~6倍方可完整探测,类似的,采用水平电偶源赤道偶极装置测量垂直磁场也能达到与之相当的探测能力.计算结果证实了21点Talbot算法适应不同地电模型、不同观测方式的时间域电磁响应计算.