摘要
在瞬变电磁测深正演中,时间域电磁场响应通常由频率域响应的逆Laplace变换得到,但常规的逆Laplace变换算法不能保证在大的计算时间范围内满足精度要求,有时甚至不能收敛。为此,对Laplace逆变换的求积公式进行了变换,得到与正弦有关的积分公式,然后对核函数进行离散,利用复数的围道积分方法推导得到快速正弦变换的算法及滤波系数,实现了逆Laplace变换的快速计算。通过具有解析解的算例,总结了多种逆Laplace变换算法的特点及适应性。将上述算法应用于瞬变电磁时间域响应的正演计算,得到了大时间范围的收敛解,解决了常规算法在大时间值时正演计算精度和收敛性差的问题。
Time domain response of transient electromagnetic sounding in forward modeling are often obtained from frequency responses by means of inverse Laplace transform.However,conventional inverse Laplace transform cannot ensure to satisfy the requirement of accuracy in a large range of computation time,sometimes it even cannot converge.Therefore,the inverse Laplace transform was alternated to obtain the integral formula related to sine.Then,the kernel function was discrete,and contour integral method of complex was used to reduce the algorithm of fast Sine transform and filtering coefficient and realize the fast computation of inverse Laplace transform. By the example with analytic solution,the characteristics and adaptability of several inverse Laplace transform methods were concluded. The algorithm was used for forward modeling of the response of transient electromagnetic time domain,and converging of large time range was obtained,which overcame the bad forward computation accuracy and convergence of conventional methods under large time value.
出处
《石油物探》
EI
CSCD
北大核心
2011年第3期213-218,4,共6页
Geophysical Prospecting For Petroleum
基金
国家重点基础研究发展计划(973)项目(2007CB209607)
国家自然科学基金项目(40727001
44704073)
高等学校博士学科点科研专项基金(20070489001)联合资助。
关键词
瞬变电磁响应
正演
逆Laplace变换
正弦变换
数值滤波算法
transient electromagnetic response
forward modeling
inverse Laplace transform
Sine transform
numerical filtering algorithm
