Three-dimensional tomography using high-power induced polarization with the similar central gradient array
Wang Jun-Lu1, Lin Pin-Rong1, Wang Meng2, Li Dang1, and Li Jian-Hua1
1. Institute of Geophysical and Geochemical Exploration, CAGS, Langfang 065000, China.
2. China Aero Geophysics Survey & Remote Sensing Center foe Land and Resources, Beijing 100083, China.
Abstract Induced polarization (IP) 3D tomography with the similar central gradient array combines IP sounding and IP profiling to retrieve 3D resistivity and polarization data rapidly. The method is characterized by high spatial resolution and large probing depth. We discuss data acquisition and 3D IP imaging procedures using the central gradient array with variable electrode distances. A 3D geoelectric model was constructed and then numerically modeled. The data modeling results suggest that this method can capture the features of real geoelectric models. The method was applied to a polymetallic mine in Gansu Province. The results suggest that IP 3D tomography captures the distribution of resistivity and polarization of subsurface media, delineating the extension of abrupt interfaces, and identifies mineralization.
The project was funded jointly by the National High Technology Research and Development Program (863 Program: No. 2014AA06A610), special funds for basic scientific research business expenses of the Chinese Academy of Geological Sciences (No. YYWF201632), and the National Major Scientific Instruments and Equipment Development Projects (No. 2011YQ050060).
Cite this article:
. Three-dimensional tomography using high-power induced polarization with the similar central gradient array[J]. APPLIED GEOPHYSICS, 2017, 14(2): 291-300.
[1]
Cardarelli, E., Cercato, M., and Di Filippo, G., 2007, Assessing foundation stability and soil structure interaction through integrated geophysical techniques: a case history in Rome (Italy): Near Surface Geophysics, 5, 141-147.
[2]
Chen, D. P., Dai, Q. W., Liu, H. F., and Feng, D. S., 2014, 3-D forward modeling and analysis of borehole-surface potential with gradient source: Journal of Central South University (Science and Technology) (inChinese), 45(1), 150-156.
[3]
Dahlin, T., and Bernstone, C., 1997, A roll-along technique for 3-D resistivity data acquisition with multi-electrode arrays, Prccs. SAGEEP’97 (Symposium on the Application of Geophysics to Engineering and Environmental Problems), Reno, Nevada, 23-26 March 1997, 2, 927-935.
[4]
Eaton, P., Anderson, B., Queen, S., et al., 2010, NEWDAS-the Newmont distributed IP data acquisition system: 80th SEG Annual Meeting, Expanded Abstracts, 1768-1772.
[5]
Johnson, T. C., Versteeg, R. J., Ward, A., Day-Lewis, F. D., and Revil, A., 2010, Improved hydrogeophysical characterization and monitoring through high performance electrical geophysical modelling and inversion: Geophysics, 75(4), WA27-WA41.
[6]
Karaoulis, M., Revil, A., Tsourlos, P., Werkema, D. D., and Minsley, B.J., 2013, IP4DI: A software for time-lapse 2D/3-D DC-resistivity and induced polarization tomography: Computers & Geosciences, 54, 164-170.
[7]
Kim, J. H., Yi, M. J., Park, S. G., and Kim, J. G., 2009, 4-D inversion of DC resistivity monitoring data acquired over a dynamically changing earth model: Journal of Applied Geophysics, 68(4), 522-532.
[8]
LaBrecque, D. J., and Yang, X., 2001, Difference inversion of ERT data: a fast inversion method for 3-D in situ monitoring: Journal of Environmental and Engineering Geophysics, 5, 83-90.
[9]
Li, Y., and Oldenburg, D., 2000, 3-D inversion of induced polarization data: Geophysics, 65, 1931-1945.
[10]
Lin, P. R., Guo, P., Shi, F. S., Zheng, C. J., Li, Y., Li, J. H., and Xu, B. L., 2010, A study of the techniques for large depth and multi-function electromagnetic survey: Acta Geoscientica Sinica (inChinese), 31(2),149-154.
[11]
Lin, P. R., Zheng, C. J., Shi, F. S., Guo, P., Xu, B. L., and Zhao, Z. Y., 2006, The research of integrated electromagnetic method system: Acta Geologica Sinica (in Chinese), 80(10),1539-1548.
[12]
Loke, M. H., and Barker, R. D., 1996a, Practical techniques for 3-D resistivity surveys and data inversion: Geophysical Prospecting, 44, 499-523.
[13]
Loke, M. H., and Barker, R. D., 1996b, Rapid least-squares inversion of apparent resistivity pseudosections using a quasi-Newton method: Geophysical Prospecting, 44, 131-152.
[14]
Loke, M. H., 2002, RES3-DINV - Rapid 3-D Resistivity and IP Inversion using the least squares method: Geoelectrical Imaging 2D & 3-D: Computer software manual, Geotomo Software, Malaysia.
[15]
Loke, M. H., and Dahlin, T., 2002, A comparison of the Gauss-Newton and quasi-Newton methods in resistivity imaging inversion: Journal of Applied Geophysics, 49, 149-162.
[16]
Ma, W., and Lin, X., 2015, The application of 3D time domain induced polarization to Kalagailei Copper-Gold mine prospecting in Xinjiang: Chinese Journal of Engineering Geophysics (in Chinese), 12(2),171-175.
[17]
Miller, C. R., Routh, P. S., Brosten, T. R., and McNamara, J. P., 2008, Application of time lapse ERT imaging to watershed characterization: Geophysics, 73, G7-G17.
[18]
Oldenburg, D., and Li, Y., 1994, Inversion of induced polarization data: Geophysics 59, 1327-1341.
[19]
Pidlisecky, A., Haber, E., and Knight, R., 2007, RESINVM3-D: a 3-D resistivity inversion package: Geophysics 72(2), H1-H10.
[20]
Pollock, D., and Cirpka, O. A., 2012, Fully coupled hydrogeophysical inversion of a laboratory salt tracer experiment monitored by electrical resistivity tomography: Water Resources Research, 48, W01505.
[21]
Ruan, B. Y., and Xiong, B., 2002, A finite element modeling of three-dimensional resistivity sounding with continuous conductivity: Chinese Journal of Geophysics (inChinese), 45(1), 131-138.
[22]
Santarato, G., Ranieri, G., Occhi, M., Morelli, G., Fischanger, F. and Gualerzi, D., 2001, Three-dimensional Electrical Resistivity Tomography to control the injection of expanding resins for the treatment and stabilization of foundation soils: Engineering Geology, 119, 18-30.
[23]
Sun, J. J., Li, Y. G. and Nabighian, M., 2012, Lithology differentiation based on inversion of full waveform induced polarization data from Newmont Distributed IP Data Acquisition System (NEWDAS): 82nd SEG Annual Meeting, Expanded Abstracts, 1-5.
[24]
Webb, D., Rowston, P., and McNeill, G., 2003, A Comparison of 2D and 3-D IP from Copper Hill NSW: ASEG 16th Geophysical Conference and Exhibition, Adelaide, February.
[25]
White, M., and Gordon, N., 2003, Deep imaging: New technology lowers cost of discovery: Canadian Mining Journal, 124(3), 27-28.
[26]
White, R. M. S., Collins, S., Denne, R., Hee, R., and Brown, P., 2001, A new survey design for 3-D IP inversion modelling at Copper Hill: Exploration Geophysic, 32, 152-155.
[27]
Xu, S. Z., 1994, The Finite Element Method in Geophysics (in Chinese): Science Press, Beijing, 178-188.
[28]
Xu, S. Z., Liu, B., and Ruan, B. Y., 1994, The finite element method for solving anomalous potential for resistivity surveys: Chinese J. Geophys. (in Chinese), 37(S2), 511-515.
[29]
Yang, J. G., Zhai, J. Y., Yang, H. W., Wang, X. H., Xie, C. L., Wang, X. A., and Ren, B. C., 2010a, Metallotectonics and Prospection of the Huaniushan Exhalogene Gold Silver Lead Zinc Deposit in Beishan, Gansu Province: Geotectonica et Metallogenia (in Chinese), 34(2), 246-254.
[30]
Yang, J. G., Zhai, J. Y., Yang, H. W., Wang, C. F., Xie, C. L., Wang, X. H., and Lei, Y. X., 2010b, LA-ICP-MS zircon U-Pb dating of basalt and its geological significance in Huaniushan Pb-Zn deposit, Beishan area, Gansu, China: Geological Bulletin of China (in Chinese), 29(7), 1017-1023.
[31]
Zhou, B., and Greenhalgh, S. A., 2001, Finite element three-dimensional direct current resistivity modeling: accuracy and efficiency considerations: Geophysical Journal International, 145, 679-688.
2017, Vol. 14 
