一种高精度瞬变电磁响应正演的数值滤波算法

摘要

在瞬变电磁测深正演中,时间域电磁场响应通常由频率域响应的逆Laplace变换得到,但常规的逆Laplace 变换算法不能保证在大的计算时间范围内满足精度要求,有时甚至不能收敛.为此,对Laplace逆变换的求积公式进行了变换,得到与正弦有关的积分公式,然后对核函数进行离散,利用复数的围道积分方法推导得到快速正弦变换的算法及滤波系数,实现了逆Laplace变换的快速计算.通过具有解析解的算例,总结了多种逆Laplace 变换算法的特点及适应性.将上述算法应用于瞬变电磁时间域响应的正演计算,得到了大时间范围的收敛解,解决了常规算法在大时间值时正演计算精度和收敛性差的问题.