长偏移距瞬变电磁测深甚晚期响应及视电阻率数值计算.pdfVIP

第 25 卷 　第 4 期 地　震　地　质 Vol. 25 , No. 4 2003 年 12 月 SEISMOLO GY AND GEOLO GY Dec. , 2003 长偏移距瞬变电磁测深甚晚期响应 及视电阻率的数值计算 翁爱华　王雪秋 (吉林大学地球探测科学与技术学院 ,长春 　130026) 摘 　要 　　文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计 算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域 甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬 变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的 ( ) 方法可有效地提高长偏移距瞬变电磁测深 甚 晚期数据的模拟与解释精度。 关键词 　　长偏移距瞬变电磁测深 　甚晚期 　数值积分 　全区视电阻率 　连分式 + ( ) 中图分类号: P631. 3 25 　　　文献标识码 : A 　　　　文章编号 : 0253 - 4967 2003 04 - 0664 - 07 0 　引言 ( ) ( 长偏移距瞬变电磁测深 LO TEM 是研究深部电性结构的一个有效方法 Strack , 1992 ;严 良俊等,2001) 。目前 ,一般观测小于 10s 的瞬变响应。但当存在较厚的导电性较好的地层时 , 为了观测到较深处地层的响应 ,需要观测更晚的瞬变电磁信号。然而 ,数值结果发现 ,在更晚的 ( ) 时间里 , 目前理论计算的LO TEM 响应具有振荡特征 Gubatenko et al . , 1999 ,原因是频率域 响应在低频时的计算精度不够。这是由于频率域电磁响应是以 Hankel 积分形式给出的 ,而目 ( ) 前的计算方法主要是数字滤波方法 Ghosh ,1971 ; Anderson , 1979 。但数字滤波本身计算精度 ( 不够高 ,且滤波系数计算精度取决于具体的计算方法 ,只有在合适的条件下才有足够的精度 朴 化荣 ,1990) 。低精度的频率响应会导致变换得到的瞬变响应尤其是晚期响应出现振荡 ,相应的 ( ) 全区视电阻率曲线也出现振荡与畸变 翁爱华 ,1996 。显然 ,沿用先计算频率响应再经变换获 得时间域响应的计算思路 ,为提高层状模型 LO TEM 晚期响应的精度 ,将晚期计算时间加以延 长 ,需要提高频率域响应的计算精度。 基于高精度的 LO TEM 电磁场分量可得到定性反映地电模型的全区视电阻率。全区视电 ( ) 阻率最成功的计算方法是多项式逼近算法 长谷川健 ,1985 ,但多项式逼近需要采用合适的多 (