摘要 在对浙江省某铜矿进行瞬变电磁勘查时，由于感应激发极化效应的影响，使得瞬变电磁晚期测道的数据发生倒转。采用常规的瞬变电磁数据处理方法难以处理，造成晚期测道数据不可用，影响TEM的探测深度及精度。针对这个情况，采用Cole-Cole模型分析了均匀半空间模型中直流电阻率、充电率、时间常数以及频率相关系数对磁性源瞬变电磁响应的影响规律。利用奇异值分解法(The Singular Value Decomposition,简称SVD)对实测TEM数据进行分析及反演，并且从瞬变电磁响应数据中分离出勘探区瞬变电磁测量数据Cole-Cole模型的各参数，将分离出来的参数用于探测结果的辅助解释。经过分析得出充电率和频率相关系数对瞬变电磁响应的影响较大,而直流电阻率和时间常数对其影响较小的结论。通过常规方法和奇异值分解法对实测数据分别进行处理，发现后者视电阻率断面图的异常更加突出，从而较准确的圈定了铜矿体的范围，此结论得到了钻孔的验证，与实际地质情况吻合较好。
Abstract：
During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the induced polarization (IP) effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition (SVD) is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.
YU Chuan-Tao,LIU Hong-Fu,ZHANG Xin-Jun et al. The analysis on IP signals in TEM response based on SVD[J]. APPLIED GEOPHYSICS, 2013, 10(1): 79-87.
[1]
El-Kaliouby, H. M., El-Diwany, E. A., Hussain, S. A., et al., 1997, Optimum negative response of a coincident-loop electromagnetic system above a polarizable half-space: Geophysics, 62, 75 - 79.
[2]
Fitterman, D. V., and Anderson, W. L., 1987, Effect of transmitter turn-off time on transient soundings: Geoexploration, 24, 131 - 146.
[3]
Flis, M. F., Newman, G. A., and Hohmann, G. W., 1989, Induced-polarization effects in time-domain electromagnetic measurements,Geophysics, 54(4), 514 - 523.
[4]
Flores, C., and Peralta-Ortega, S. A., 2009, Induced polarization with in-loop transient electromagnetic soundings: A case study of mineral discrimination at El Arco porphyry copper, Mexico: Journal of Applied Geophysics, 68, 423 - 436.
[5]
Ghorbani, A., Camerlynck, C., Florsch, N., Cosenza, P., and Revil, A., 2007, Bayesian inference of the Cole-Cole parameters from time- and frequency-domain induced polarization: Geophysical Prospecting, 55, 589 - 605.
[6]
Guptasarma, D., 1982, Computation of time-domain response of a polarizable ground: Geophysics, 47(11), 1574 - 1576.
[7]
Han, Y. L., 2006, Application of transient electromagnetic induced IP: Chinese Journal of Engineering Geophysics, 3(5), 366 - 369.
[8]
Jiang, R. Z., Yao Y. P., et al., 2005, Improved algorithm for singular value decomposition inversion of T2 spectrum in nuclear magnetic resonance: Acta Petrolei Sinica, 11(6), 27 - 59.
[9]
Jupp, D. L. B., and Vozof, K. F., 1975, Stable iterative methods for the inversion of geophysical data: Geophysical Journal of the Royal Astronomical Society, 42(3), 957 - 976.
[10]
Kozhevnikov, N. O., and Antonov, E. Y., 2009, Joint inversion of IP-affected TEM data: Russian Geology and Geophysics, 50, 136 - 142.
[11]
Kozhevnikov, N. O., and Antonov, E. Y., 2008, Inversion of TEM data affected by fast-decaying induced polarization: Numerical simulation experiment with homogeneous half-space: Journal of Applied Geophysics, 66, 31 - 43.
[12]
Morrison, H. F., 1969, Quantitave interpretation of transient electromagnetic fields over a layered half space: Geophysical Prospecting, 17, 82 - 101.
[13]
Niu, Z. L., 2007, Theory of time 2 domain electromagnetic methods: Changsha Central South University of Technology Press (in Chinese).
[14]
Pelton, W. H., Ward, S. H., Hallof, P. G., Sill, W. R., and Nelson, P. H., 1978, Mineral discrimination and removal of inductive coupling with multi-frequency IP: Geophysics, 43, 588 - 609.
[15]
Reninge, P. A., Martelet, G., Deparis, J., Perrin ,J., Chen, Y.,2011, Singular value decomposition as a denoising tool for airborne time domain electromagnetic data: Journal of Applied Geophysics, 75, 264 - 276.
[16]
Sidorov, V. A., and Yakhin, A. M., 1979, Induced polarization of rocks in TEM soundings: Izv. An SSSR. Ser. Fizika Zemli, 11, 46 - 52.
[17]
Spies, B. R., 1980, A field occurrence of sign reversals with the transient electromagnetic method: Geophysical Prospecting, 28, 620 - 632.
[18]
Stoyer, C. H., 1990, Efficient computation of transient sounding curves for wire segments of finite length using an equivalent dipole approximation: Geophysical Prospecting, 38, 87 - 100.
[19]
Wang, L. P., and Wen, P. L., 1998, IP effects in TEM response: Journal of Central South University Technology, 29(3), 208 - 210.
Yin, C. C., and Liu, B., 1994, Research on the 3D TDEM modeling and IP effect: Chinese J. Geophys, (in Chinese), 37(Suppl. 2), 486 - 4927.
[22]
Yue, J. H., Yang, H. Y., 2008, Research on response of transient electromagnetic field in underground mine with boundary condition of laneway: Journal of China University of Mining & Technology, 37(2), 152 - 156.
2013, Vol. 10 
