TEM中心回线法瞬变磁场求取和全区视电阻率计算

1、 第 6期 许建荣等 : TEM 中心回线法瞬变磁场求取和全区视电阻率计算 全区视电阻率。不同关断时间的中晚期道的视电阻 率几乎相等 , 差别在早期道 , 对应三个关断时间的第 一道视 电 阻 率 分 别 为 87 . 5 80. 0 m, 83 . 1 m, m, 关断时间越长, 第一道视电阻率与模型 = S! 2 SB 13( toff + 2 t I! 1 /2 当 t/ 0 时, B / I !/ 2 ( ( /S , / 0 , 所以远区 成份多, 全区视电阻率会降低。 第一层电阻率相差越大。其原因是早期道的首支视 电阻率因关断延时不同程度地含有第二层的低阻信 息 , 使早期道的视电

2、阻率曲线变平滑。算例说明虽 然全区视电阻率考虑了关断时间, 但计算的早期道 视电阻率仍然受到关断时间的影响。 4 结语 1 文章依据 d B / dt 的衰减规律给出了由实测 的 dB / dt计算瞬变磁场 B 的方法, 通过在半空间条 件下的试算, 完成满足精度要求。但对于实测资料, 因晚期数据受干扰很大, 必须要根据变化趋势认真 编辑后才能计算。积分 常数 C 的 误差主要影响晚 期道视电阻率的计算精度。 2 计算 全区视电 阻率对 TEM 资料 的定 性分 析, 半定量解释及一维反演等很有意义 , 文章所介绍 的全区视电阻率计算方法, 尽量用解析关系计算反 演过程中的积分和求导 , 这样

3、不仅提高了计算精度, 也提高了计算速度 , 对于含有概率积分的函数, 采用 建立数据库的方法 , 使 Z 和逻辑记录顺序建立起关 系, 可直接读出 Z 所在的区间, 进行内差, 这也提高 计算速度的关键, 所以本章的方法应该是比较实用 图 8 不同关断时间全区视电阻率曲线 ( 发射框 : 200m ( 200m 和易实现可操作的。 3 由 B 计算的全区视电阻率 曲线和由 d B / dt 计算的全区视电阻率 ! 曲线都能反映地下电性变化 特征 , 但有区别的。相同衰减延时时, 反映的深度 曲线上 大于 ! 的, 且因 B 对电阻率是个单调函数, 不会出 现 Overshoo t.和 U nd

4、ershoot. 现象, 但因 比较平滑 , 可能会模糊一些细节。由算例表明, 在早 期道 , 受关断时间的影响大于 ! ,而 ! 受 发射框大 小的影响大于 。 参考文献 1 王庆 乙 . 近 区 瞬变 场的 频率 效应 与 大深 度 探测 能 力的 研究 J . 地质与勘探 , 2007 , 43 ( 1 : 79 - 83. 2 罗国平 , 魏凤英 . 瞬变电磁法磁场及其视电阻 率研究 J . 中国 煤田地质 , 2000 , 12( 2 : 60- 62 . 3 白登海 , M axw ellM eju. 瞬 变电磁 法中两 种关断 电流对 响应函 数的影响及其 对应 策 略 J .

5、地 震地 质 , 2001 , 23 ( 2 : 245 251. 4 苏朱刘 , 张交东 , 胡文 宝 . 瞬变电 磁法中 磁场的 特性及 测量方 法 J . 工程地球物理学报 , 2004, 1 ( 1 : 70 - 73. 5 蒋邦远主编 . 实用近区磁源瞬变电 磁法勘探 M . 北京 : 地质 出版社 , 1998. 6 R aiche A P. Th e u se and m isuse of apparen t res is tiv ity in electro m agnet ic m ethod s J . G eophys ics , 1984 , 15 : 37 - 41

6、. 7 陈晓东 , 赵毅 , 王赤 军 , 等 . 高温 超导 磁强计 的研 制及 在 TE M 上的应用研究 J . 地质与勘探 , 2002, 38( 增刊 : 18 - 22 . 8 牛之琏 . 时间域电磁法原理 M . 长沙 : 中南 工业大学出版社 , 图 9 不同发射框大小全区视电阻率曲 线 ( 关断时间 : 0. 3m s 3 不同发射框大小的全区视电阻率比较 对同一五层 模型在关断时间 为 0 . 3m s 的条件 下 , 发射框大小分别为 200 m ( 200m 和 400 m ( 400m 时 , 计算的全区视电阻率见图 9 , 差别也主要是在早 期道上 , 第一 道 视

7、电 阻 率 分别 为 76 . 08 m和 71. 78 m。发射框越大 , 实测数据中远区场成份 越多, 在早期道, 视电阻率 和 B 有如下关系: 67 地质与勘探 1992 . 9 10 K au f m an A A. H ar m on ic and transient fields on the su rface of a t wo layer m ed ium J. G eophysics , 1979, 44 : 1208 - 1217 . F itter m an D V and A nderson W L. E ffect of tran s m it ter tu rn

8、off ti m e on tran sient sound ing J . G eoexp loration , 1987 , 24: 131 11 146. 2008 年 R a iche A P. C omp arison of apparent res ist ivitys functions for trans ien t electrom agneticm ethods J . G eophysics , 1983 , 48 : 787 - 789. CALCULATION OF TRANSIENT M AGNETIC FIELD AND ALL TIM E APPARENT RE

9、SISTIV ITY BASED ON CENTRAL TEM LOOPS M ETHOD XU Jian- rong1, 2 , L I A i- yong1, 2, YANG Sheng3 ( 1. School of Info - phy sics E ngineer ing, Cen tral Sou th U niver sity, Chang sha 2. E astern Ch ina Bureau of N onf errous G eolog ical Exp loration, N anj ing 3. Ch ina N onferrous M etals R esourc

10、e G eolog ical Survey, B eij ing 410083; 210007; 100012 Abstract : A ccorded on the attenuation reg i m e of cen tral loop un it TEM field , th e exp licit calcu lation m ethod of trans ien t m agnet ic f ield B , from dB / d t m easurem ent w as g iven. The m ethod , w h ich inverses the all ti m e

11、 apparen t resist ivity from B iteration , under the s lan ting w ave tran s m iss ion curren t w as deduced. The in tegral and derivat ion resolu tion expression w as d edu ced, in order to prom ote th e calcu lat ion precis ion and speed , as w ell as the calcu lation of p robab il ity in tegral function u sing stan dard databasem eth od . Then through the sam ple calcu