Three-dimensional finite difference time domain simulation for down-hole transient electromagnetic method
XU Xin-Gang1, XU Zheng-Yu2
1. School of Resources and Geosciences, China University of Mining & Technology,Xuzhou 221116,China;; 2. School of Nuclear Engineering and Geophysics, East China Institute of Technology, Nanchang 330013,China
Abstract:As one of the effective methods used in geological prospecting, the down-hole transient electromagnetic method is a geophysical survey method in search for mineral resources, which measures the transient electromagnetic field in a borehole by placing the receiving coils in the borehole. In this paper, based on time domain finite difference algorithm, the authors used a loop source and non-uniform gird, developed a 3D simulation program, and simulated the response of the uniform medium half space, inclined conductor, and inclined conductor with overburden. According to the comparative analysis, the maximum value of response decreases with delay time in uniform half space; in the early period of transient electromagnetic survey, it mainly reflects the response of the surrounding rock; in the late period, abnormal response of the conductor can be obviously displayed. When the source is in the central position, the abnormal response conductor can be reflected, and when the source is in the right position, the need to extend the observation time can be observed in response to abnormal conductor. The impact on overburden conductor is reflected in the entire observation time, making the abnormal response of the conductor ambiguously displayed. The authors studied the work and provide the reference for the down-hole TEM abnormal response characteristics and data interpretation.
[1] 蒋邦远.实用近区磁源瞬变电磁法勘探[M]. 北京: 地质出版社,1998. [2] 牛之琏.时间域电磁法勘探[M]. 长沙:中南大学出版社,2007. [3] Dyck A V, West G F. The role of simple computer models in interpretations of wide-band, drill-hole electromagnetic surveys in mineral exploration[J]. Geophysics, 1984, 49(7):957-980. [4] Macnae J, Staltari G. Classification of sign changes in Borehole TEM decays[J]. Exploration Geophysics, 1987, 18(3): 331-339. [5] Eaton P A, Hohmann G W. The influence of a conductive host on two-dimensional borehole transient electromagnetic responses [J]. Geophysics, 1984, 49(7):861-869. [6] West R C, Ward S H. The borehole transient electromagnetic response of a three-dimensional fracture zone in a conductive half-space[J]. Geophysics, 1988, 53(11):1469-1478. [7] 张杰,邓晓红,郭鑫.等. 地—井TEM在危机矿山深部找矿中的应用实例[J]. 物探与化探,2013,37(1):30-34. [8] 张杰. 地—井瞬变电磁异常特征分析及矢量交会解释方法[D]. 北京:中国地质大学(北京),2009. [9] 戴雪平.地—井瞬变电磁法三维响应特征研究[D]. 北京:中国地质大学(北京),2013. [10] 孟庆鑫, 潘和平. 地—井瞬变电磁响应特征数值模拟分析[J]. 地球物理学报, 2012, 55(3):1046-1053. [11] 宋维琪. 3D瞬变电磁场有限差分正演计算[J].石油地球物理勘探, 2000,35(6):751-756. [12] 闫述,陈明生,傅君眉.瞬变电磁场的直接时域数值分析[J].地球物理学报, 2002,45(2): 275-284. [13] 杨海燕,岳建华. 磁偶源2.5维瞬变电磁场全空间FDTD数值模拟[J]. 物探与化探, 2008, 32(3):326-330. [14] 岳建华,杨海燕,胡搏.矿井瞬变电磁法三维时域有限差分数值模拟[J].地球物理学进展,2007, 22(6):1904-1909. [15] 董浩,魏文博,叶高峰,等.基于有限差分正演的带地形三维大地电磁反演方法[J].地球物理学报,2014,57(3):939-952. [16] 孙怀风,李貅,李术才,等.考虑关断时间的回线源激发TEM三维时域有限差分正演[J].地球物理学报,2013,56(3):1049-1064. [17] 谭捍东,余钦范,John Booker,等. 大地电磁法三维交错采样有限差分数值模拟[J].地球物理学报,2003,46(5):705-711. [18] 杨海燕,邓居智,张华,等. 矿井瞬变电磁法全空间视电阻率解释方法研究[J]. 地球物理学报,2010,53(3):651-656. [19] 米萨克 N.纳比吉安.勘查地球物理电磁法:第1卷[M].赵经祥,王艳君,译.北京:地质出版社,1992. [20] 杨海燕.矿用多匝小回线源瞬变电磁场数值模拟与分布规律研究[D].徐州:中国矿业大学,2009. [21] 廖振鹏,周正华,张艳红.波动数值模拟中透射边界的稳定实现[J].地球物理学报,2002,45(4):533-545.