Study on the Coincident-Loop Transient Electromagnetic Method in Seafloor Exploration—Taking Jiaodong Polymetallic Mine as a Model

Abstract

The transient electromagnetic (TEM) method becomes more urgent than ever for marine exploration due to abundant resource reserves and the increasing undersea engineering construction activities, especially in the offshore exploration of mineral deposits such as Sanshandao gold mine. However, the research and application of TEM method in marine environment are still challenged by many problems. Such contradiction motivates our study on the coincident-loop TEM in seafloor exploration. The TEM response of coincident loops is firstly derived in the integral form, based on the potential functions in Helmholtz equations for a magnetic source locating in the whole-space layered model. The frequency-domain vertical magnetic field is described as the Hankel integral with double first-order Bessel functions of first kind. Secondly, the time-domain induced voltage is obtained by transforming the frequency-domain response through the cosine transform and then taking the derivative of time. To simultaneously solve the Hankel transform and the cosine transform, a novel algorithm is introduced by adapting the fixed-point quadrature and extrapolation via the Shanks transformation. Finally, a typical conductivity model for marine polymetallic deposit is designed to investigate the characteristic of TEM response under various conditions. Numerical results demonstrate that existence of conductive seawater causes the TEM response to increase significantly and decay slower. The air-sea reflected electromagnetic waves lead to a significantly large fake negative response (NR) in shallower seawater with depth less than 300 m. Increase in the height of loops will weaken and delay the anomaly response and shorten the observation time-window. The height of configuration should be no more than 100 m for shallower targets and 50 m for deeper targets, respectively. The observation time-window should cover 10–1 000 ms. Increase in the radius of loops only enhances the TEM response proportionally but hardly improves the relative anomaly. The vertical resolution on the low-resistivity target approximates 20 m for the configuration considered in the study. Decreases in D.C. resistivity and chargeability cause the positive response (PR) to increase significantly and decay more rapidly. Meanwhile, the NR is advanced and enlarged significantly and decays slower compared with the PR. The influence of time constant is not monotony and there exists an optimal value for producing the maximum NR. As the frequency parameter increases, the PR is caused to decay more rapidly without magnitude change and the NR is advanced and decays more rapidly with significant increase in magnitude. The influence of frequency parameter is more pronounced than that of time constant.

Keywords

marine transient electromagnetic method / coincident loops / induced polarization / airwave effects

References

[1] Anderson W L. Computation of Green’s Tensor Integrals for Three-Dimensional Electromagnetic Problems Using Fast Hankel Transforms. Exploration Geophysics, 1984, 15(3): 195-195.

[2] Anderson W L. A Hybrid Fast Hankel Transform Algorithm for Electromagnetic Modeling. Geophysics, 1989, 54(2): 263-266.

[3] Chang Y, Xiao M, Wu Y. One-Dimensional Inversion Strategies for Transient Electromagnetic Data Based on Apparent Vertical Conductance. Oil Geophysical Prospecting, 2010, 45(2): 295-298. (in Chinese with English Abstract)

[4] Chang Y, Zhang G. Comparison among Three Transformation Algorithms of Electromagnetic Field from Frequency Domain to Time Domain. Computing Techniques for Geophysics and Geochemical Exploration, 1995, 3(17): 25-29. (in Chinese with English Abstract)

[5] Chave A D. Numerical Integration of Related Hankel Transforms by Quadrature and Continued Fraction Expansion. Geophysics, 1983, 48(12): 1671-1686.

[6] Chave A D, Cox C S. Controlled Electromagnetic Sources for Measuring Electrical Conductivity beneath the Oceans: 1. Forward Problem and Model Study. Journal of Geophysical Research: Solid Earth, 1982, 87(B7): 5327-5338.

[7] Cheesman S J, Edwards R N, Chave A D. On the Theory of Sea-Floor Conductivity Mapping Using Transient Electromagnetic Systems. Geophysics, 1987, 52(2): 204-217.

[8] Chen S, Li X, Duan W, . Petrological and Geochronological Study of Amphibolite from Jiaobei Terrane. Journal of Earth Science, 2018, 43(3): 716-732. (in Chinese with English Abstract)

[9] Commer M, Newman G A. New Advances in Three-Dimensional Controlled-Source Electromagnetic Inversion. Geophysical Journal International, 2008, 172(2): 513-535.

[10] Constable S C. Marine Electromagnetic Induction Studies. Surveys in Geophysics, 1990, 11(2): 303-327. 3

[11] Constable S C, Orange A S, Hoversten G M, . Marine Magnetotellurics for Petroleum Exploration Part I: A Sea-floor Equipment System. Geophysics, 1998, 63(3): 816-825.

[12] Constable S, Srnka L J, . An Introduction to Marine Controlled-Source Electromagnetic Methods for Hydrocarbon Exploration. Geophysics, 2007, 72(2): WA3-WA12.

[13] Descloitres M, Guérin R, Albouy Y, . Improvement in TDEM Sounding Interpretation in Presence of Induced Polarization. a Case Study in Resistive Rocks of the Fogo Volcano, Cape Verde Islands. Journal of Applied Geophysics, 2000, 45(1): 1-18.

[14] El-Kaliouby H M, Hussain S A, Bayoumi A E R, . Effect of Clayey Media Parameters on the Negative Response of a Coincident Loop1. Geophysical Prospecting, 1995, 43(5): 595-603.

[15] Ellis M, Evans R L, Hutchinson D, . Electromagnetic Surveying of Seafloor Mounds in the Northern Gulf of Mexico. Marine and Petroleum Geology, 2008, 25(9): 960-968.

[16] Evans R L, Everett M E. Discrimination of Hydrothermal Mound Structures Using Transient Electromagnetic Methods. Geophysical Research Letters, 1994, 21(6): 501-504.

[17] Everett M E. Transient Electromagnetic Response of a Loop Source over a Rough Geological Medium. Geophysical Journal International, 2009, 177(2): 421-429.

[18] Fiandaca G, Ramm J, Binley A, . Resolving Spectral Information from Time Domain Induced Polarization Data through 2-D Inversion. Geophysical Journal International, 2012, 192(2): 631-646.

[19] Flis M F, Newman G A, Hohmann G W. Induced-Polarization Effects in Time-Domain Electromagnetic Measurements. Geophysics, 1989, 54(4): 514-523.

[20] Flosadóttir Á H H, Constable S. Marine Controlled-Source Electromagnetic Sounding: 1. Modeling and Experimental Design. Journal of Geophysical Research: Solid Earth, 1996, 101(B3): 5507-5517.

[21] Ghosh D P. The Application of Linear Filter Theory to the Direct Interpretation of Geoelectrical Resistivity Sounding Measurements. Geophysical Prospecting, 1971, 19(2): 192-217.

[22] Goldfarb R J, Santosh M. The Dilemma of the Jiaodong Gold Deposits: Are They Unique?. Geoscience Frontiers, 2014, 5(2): 139-153.

[23] Guo T, Lü G X. Measurement of Metallogenic Depth and Deep Second Enrichment Belt Forecasting in Jiaodong Gold Deposit. Earth Science—Journal of China University of Geosciences, 2018, 6(2): 102-107. (in Chinese with English Abstract)

[24] Guptasarma D. Computation of the Time-Domain Response of a Polarizable Ground. Geophysics, 1982, 47(11): 1574-1576.

[25] Guptasarma D, Singh B. New Digital Linear Filters for Hankel J₀ and J₁ Transforms. Geophysical Prospecting, 1997, 45(5): 745-762.

[26] Hohmann G W, Kintzinger P R, Van Voorhis G D, . Evaluation of the Measurement of Induced Electrical Polarization with an Inductive System. Geophysics, 1970, 35(5): 901-915.

[27] Hohmann G W, Newman G A. Transient Electromagnetic Responses of Surficial, Polarizable Patches. Geophysics, 1990, 55(8): 1098-1100.

[28] Hu H L, Liu S L, Fan H R, . Structural Networks Constraints on Alteration and Mineralization Processes in the Jiaojia Gold Deposit, Jiaodong Peninsula, China. Journal of Earth Science, 2020, 31(3): 500-513.

[29] Hu J, Chang Y, Lei S, . 1D forward of TEM of Central Loop Configuration on Seafloor and Calculation of All-Time Apparent Resistivity. Geophysical & Geochemical Exploration, 2013, 37(6): 1137-1140. (in Chinese with English Abstract)

[30] Hu J, Chang Y, Zuo Q, . The application of Hankel Transform Based on QWE in Electromagnetic Modeling. Oil Geophysical Prospecting, 2015, 50(2): 370-375. (in Chinese)

[31] Jang H, Kim H J. Mapping Deep-Sea Hydrothermal Deposits with an In-Loop Transient Electromagnetic Method: Insights from 1D Forward and Inverse Modeling. Journal of Applied Geophysics, 2015, 123: 170-176.

[32] Johansen H, Sørensen K. Fast Hankel Transforms. Geophysical Prospecting, 1979, 27(4): 876-901.

[33] Keller G V, Kaufman A A. Frequency and Transient Sounding, 1987, Beijing: Geology Press, 134-167.

[34] Key K. 1D Inversion of Multicomponent, Multifrequency Marine CSEM Data: Methodology and Synthetic Studies for Resolving Thin Resistive Layers. Geophysics, 2009, 74(2): F9-F20.

[35] Key K. Is the Fast Hankel Transform Faster than Quadrature?. Geophysics, 2012, 77(3): F21-F30.

[36] Knight J H, Raiche A P. Transient Electromagnetic Calculations Using the Gaver-Stehfest Inverse Laplace Transform Method. Geophysics, 1982, 47(1): 47-50.

[37] Kong F N. Hankel Transform Filters for Dipole Antenna Radiation in a Conductive Medium. Geophysical Prospecting, 2007, 55(1): 83-89.

[38] Kozhevnikov N O, Antonov E Y. Fast-Decaying IP in Frozen Unconsolidated Rocks and Potentialities for Its Use in Permafrost-Related TEM Studies. Geophysical Prospecting, 2006, 54(4): 383-397.

[39] Lee T. Transient Electromagnetic Response of a Polarizable Ground. Geophysics, 1981, 46(7): 1037-1041.

[40] Li R X, Wang H, Xi Z Z, . The 3d Transient Electromagnetic forward Modeling of Volcanogenic Massive Sulfide Ore Deposits. Chinese Journal of Geophysics, 2016, 59(6): 725-733.

[41] Li Y G, Key K. 2D Marine Controlled-Source Electromagnetic Modeling: Part 1—An Adaptive Finite-Element Algorithm. Geophysics, 2007, 72(2): WA51-WA62.

[42] Li Y, Constable S. Transient Electromagnetic in Shallow Water: Insights from 1D Modeling. Chinese Journal of Geophysics, 2010, 3(53): 737-742. (in Chinese with English Abstract)

[43] Liu C S, Lin J. Transient Electromagnetic Response Modeling of Magnetic Source on Seafloor and the Analysis of Seawater Effect. Chinese Journal of Geophysics, 2006, 49(6): 1726-1736.

[44] Liu C S, Lin J, Zhou F D. Transient Response Characteristics of Central Loop Configuration on Seafloor. Dianbo Kexue Xuebao/Chinese Journal of Radio Science, 2010, 25: 195-200. (in Chinese with English Abstract)

[45] Løseth L O, Pedersen H M, Ursin B, . Low-Frequency Electromagnetic Fields in Applied Geophysics: Waves or Diffusion?. Geophysics, 2006, 71(4): W29-W40.

[46] Lucas S K, Stone H A. Evaluating Infinite Integrals Involving Bessel Functions of Arbitrary Order. Journal of Computational and Applied Mathematics, 1995, 64(3): 217-231.

[47] Luo Y Z, Chang Y J. A Rapid Algorithm for G-S Transform. Chinese Journal of Geophysics, 2000, 43(5): 724-730. in Chinese with English Abstract)

[48] Michalski K A. Extrapolation Methods for Sommerfeld Integral Tails. IEEE Transactions on Antennas and Propagation, 1998, 46(10): 1405-1418.

[49] Newman G A, Alumbaugh D L. Frequency-Domain Modelling of Airborne Electromagnetic Responses Using Staggered Finite Differences. Geophysical Prospecting, 1995, 43(8): 1021-1042.

[50] Newman G A, Hohmann G W, Anderson W L. Transient Electromagnetic Response of a Three-Dimensional Body in a Layered Earth. Geophysics, 1986, 51(8): 1608-1627.

[51] Nie J Z, Liu Y C, Yang Y. Phase Equilibria Modeling and P-T Evolution of the Mafic Lower-Crustal Xenoliths from the Southeastern Margin of the North China Craton. Journal of Earth Science, 2018, 29(5): 1236-1253.

[52] Orange A, Key K, Constable S. The Feasibility of Reservoir Monitoring Using Time-Lapse Marine CSEM. Geophysics, 2009, 74(2): F21-F29.

[53] Pelton W H, Ward S H, Hallof P G, . Mineral Discrimination and Removal of Inductive Coupling with Multifrequency IP. Geophysics, 1978, 43(3): 588-609.

[54] Qi Y Z, El-Kaliouby H, Revil A, . Three-Dimensional Modeling of Frequency- and Time-Domain Electromagnetic Methods with Induced Polarization Effects. Computers & Geosciences, 2019, 124: 85-92.

[55] Qi Y Z, Huang L, Wu X, . Effect of Loop Geometry on TEM Response over Layered Earth. Pure and Applied Geophysics, 2014, 171(9): 2407-2415.

[56] Raiche A P. Transient Electromagnetic Field Computations for Polygonal Loops on Layered Earths. Geophysics, 1987, 52(6): 785-793.

[57] Raiche A P, Bennett L A, Clark P J, . The Use of Cole-Cole Impedances to Interpret the TEM Response of Layered Earths. Exploration Geophysics, 1985, 16(2): 271-273. 3

[58] Raiche A. A Flow-through Hankel Transform Technique for Rapid, Accurate Green’s Function Computation. Radio Science, 1999, 34(2): 549-555.

[59] Seidel M, Tezkan B. 1D Cole-Cole Inversion of TEM Transients Influenced by Induced Polarization. Journal of Applied Geophysics, 2017, 138: 220-232.

[60] Shanks D. Non-Linear Transformations of Divergent and Slowly Convergent Sequences. Journal of Mathematics and Physics, 1955, 34(1): 1-42. 2/3/4

[61] Singh S, Singh R. On the Use of Shank’s Transform to Accelerate the Summation of Slowly Converging Series. IEEE Transactions on Microwave Theory and Techniques, 1991, 39(3): 608-610.

[62] Smith R S, Klein J. A Special Circumstance of Airborne Induced-Polarization Measurements. Geophysics, 1996, 61(1): 66-73.

[63] Smith R S, West G F. An Explanation of Abnormal TEM Responses: Coincident-Loop Negatives, and the Loop Effect. Exploration Geophysics, 1988, 19(3): 435-446.

[64] Song M C, Li S Z, Santosh M, . Types, Characteristics and Metallogenesis of Gold Deposits in the Jiaodong Peninsula, Eastern North China Craton. Ore Geology Reviews, 2015, 65: 612-625.

[65] Song Z J, Liu H M, Meng F X, . Zircon U-Pb Ages and Hf Isotopes of Neoproterozoic Meta-Igneous Rocks in the Liansandao Area, Northern Sulu Orogen, Eastern China, and the Tectonic Implications. Journal of Earth Science, 2019, 30(6): 1230-1242.

[66] Spies B R. A Field Occurrence of Sign Reversals with the Transient Electromagnetic Method. Geophysical Prospecting, 1980, 28(4): 620-632.

[67] Stoyer C. Electromagnetic Fields of Dipoles in Stratified Media. IEEE Transactions on Antennas and Propagation, 1977, 25(4): 547-552.

[68] Swidinsky A, Hölz S, Jegen M. On Mapping Seafloor Mineral Deposits with Central Loop Transient Electromagnetics. Geophysics, 2012, 77(3): E171-E184.

[69] Ward S H, Hohmann G W. Nabighian M N. Electromagnetic Theory for Geophysical Applications. Electromagnetic Methods in Applied Geophysics, 1988, Huston: Society of Exploration Geophysicist, 131-311.

[70] Weidelt P. Response Characteristics of Coincident Loop Transient Electromagnetic Systems. Geophysics, 1982, 47(9): 1325-1330.

[71] Weiss C J. The Fallacy of the “Shallow-Water Problem” in Marine CSEM Exploration. Geophysics, 2007, 72(6): A93-A97.

[72] Weitemeyer K A, Constable S C, Key K W, . First Results from a Marine Controlled-Source Electromagnetic Survey to Detect Gas Hydrates Offshore Oregon. Geophysical Research Letters, 2006, 33 3 L03304

[73] Weniger E J. Nonlinear Sequence Transformations for the Acceleration of Convergence and the Summation of Divergent Series. Computer Physics Reports, 1989, 10(5): 189-371. 6

[74] Xiong Z H. Electromagnetic Fields of Electric Dipoles Embedded in a Stratified Anisotropic Earth. Geophysics, 1989, 54(12): 1643-1646.

[75] Yang Q Y, Santosh M. Early Cretaceous Magma Flare-up and Its Implications on Gold Mineralization in the Jiaodong Peninsula, China. Ore Geology Reviews, 2015, 65: 626-642.

[76] Zach J J, Brauti K. Methane Hydrates in Controlled-Source Electromagnetic Surveys—Analysis of a Recent Data Example. Geophysical Prospecting, 2009, 57(4): 601-614.

[77] Zeng Q D, Liu T B, Shen Y C. The Tanlu Fault Zone and Gold Ore Metallogenesis in Eastern China. International Geology Review, 2001, 43(2): 176-190.

[78] Zhang Y, Ling W, Zhang J, . Zircon U-Pb Geochronology of the Mesozoic Volcanic Rocks from Qingshan Group in the Eastern Shandong. Earth Science, 2019, 44(1): 344-354. (in Chinese with English Abstract)

[79] Zhao Y W, Zhu Z Q, Lu G Y, . The Optimal Digital Filters of Sine and Cosine Transforms for Geophysical Transient Electromagnetic Method. Journal of Applied Geophysics, 2018, 150: 267-277.

[80] Zhou S, Song G, Huang L, . Marine Towed Transient Electromagnetic System Working in Deep-Sea below 6 000 Meters and Its Application. Chinese Journal of Geophysics, 2017, 60(11): 4294-4301.

[81] Zhou S, Xi Z Z, Song G, . Responses of the Towed Transient Electromagnetic Sounding on Deep Seafloor. Journal of Central South University (Science and Technology), 2012, 43(2): 605-610.