1. School of Geosciences and Infophysics of Central South University, Changsha 410083, China.
2. School of Geophysics and Measurement-control Technology, East China Institute of Technology, Nanchang 330013, China.
Abstract We have developed a hybrid solver that combines the finite-element and integral-equation method for 3D CSEM modeling based on unstructured meshes. To avoid the source singularity, the secondary field is used in the modeling framework. The primary electromagnetic field from an electric dipole source in a layered medium is calculated based on the magnetic vector potential method. The inhomogeneities of the computational region are discretized by a vector-based finite-element mesh with boundaries at finite distance from the inhomogeneities by using the dyadic Green’s function, reducing the truncation boundary effect and the solution region. The electric and magnetic Green’s function is used in data postprocessing to reduce the numerical errors owing to inaccurate gradients because of unstructured meshes; thus, the electromagnetic field is more accurately calculated. Finally, the proposed algorithm is applied to a block and a disc model, and we assess the topography effect on the field components.
This study was supported by the National Nature Science Foundation of China (Nos. 41830107 and 41574120), and Doctoral Student Innovation Program (No. 2016zzts086).
Cite this article:
. A hybrid finite-element and integral-equation method for forward modeling of 3D controlled-source electromagnetic induction[J]. APPLIED GEOPHYSICS, 2018, 15(3-4): 536-544.
[1]
Ansari, S., and Farquharson, C. G., 2014, 3D finite-element forward modeling of electromagnetic data using vector and scalar potentials and unstructured grids: Geophysics, 79(4), E149-E165.
[2]
Best, M. E., Duncan, P., Jacobs, F. J., et al., 1985, Numerical modeling of the electromagnetic response of three-dimensional conductors in a layered earth: Geophysics, 50(4), 665−676.
[3]
Chen, H., Deng, J. Z., Tan, H. D., et al., 2011, Study on divergence correction method in three-dimensional magnetotelluric modeling with staggered-grid finite difference method: Chinese Journal of Geophysics, 54(6), 1649-1659.
[4]
Farquharson, C. G., and Oldenburg D. W., 2002, Chapter 1An integral equation solution to the geophysical electromagnetic forward-modelling problem: Methods in Geochemistry and Geophysics, 35, 3-19.
[5]
Hohmann, G. W., 1975, Three-dimensional induced polarization and electromagnetic modeling: Geophysics, 40(2), 309-324.
[6]
Kalinkin, A., Anders, A., and Anders, R., 2015, Schur complement computations in Intel® Math Kernel Library PARDISO: Applied Mathematics, 6(2), 304-311.
[7]
Lee, K. H., Pridmore, D. F, and Morrison, H. F., 1981, A hybrid three-dimensional electromagnetic modeling scheme: Geophysics, 46(5), 796-805.
[8]
Li, G., Xiao, X., Tang, J. T., et al., 2017, Near-source noise suppression of AMT by compressive sensing and mathematical morphology filtering: Applied Geophysics, 14(4), 581-589.
[9]
Mo, D., Jiang, Q. Y., Li, D. Q., et al., 2017, Controlled-source electromagnetic data processing based on gray system theory and robust estimation: Applied Geophysics, 14(4), 570-580.
[10]
Nabighian, M. N., 1988, Electromagnetic methods in applied geophysics: theory, in Corbett, J., D. Ed., Society of Exploration Geophysicists, USA.
[11]
Ren, Z. Y., Kalscheuer, T., Greenhalgh, S., et al., 2013, A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling: Geophysical Journal International, 194(2), 700-718.
[12]
Scheen, W. L., 1978, ‘EMMMMA’ a computer program for 3D modelling of airborne EM surveys: Proc. of Workshop on Modelling of Electric and EM Methods, Univ. of California, Berkeley, 53-61.
[13]
Si, H., 2015, TetGen, a Delaunay-based quality tetrahedral mesh generator: ACM transactions on Mathematical Software, 41(2), 1-36: https://dl.acm.org/citation.cfm?doid=2629697
[14]
Tan, H. D., Yu, Q. F., Booker, J., et al., 2003, Magnetotelluric three-dimensional modeling using the staggered-grid finite difference method: Chinese Journal Geophysics, 46(5), 705-711.
[15]
Tang, J. T., Zhou, F., Ren, Z. Y., et. al., 2018, Three-dimensional forward modeling of the controlled-source electromagnetic problem based on the integral equation method with an unstructured grid: Chinese Journal of Geophysics, 61(4), 1549-1562.
[16]
Wannamaker, P. E., 1991, Advances in three-dimensional magnetotelluric modeling using integral equations: Geophysics, 56(11), 1716-1728.
[17]
Xie, G. Q., Li, J. H., Majer, E. L., et al., 2000, 3-D electromagnetic modeling and nonlinear inversion: Geophysics, 65(3), 804-822.
[18]
Ye, Y. X., Li, Y. G., Liu, Y., et al., 2016, 3D finite element modeling of marine controlled-source electromagnetic fields using locally refined unstructured meshes: Chinese Journal of Geophysics, 59(12), 4747-4758.
[19]
Yin, C. C., Zhang, B., Liu, Y. H., et al., 2016, A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling: Geophysics, 81(5), E337-E346.
[20]
Zhang, B., Yin, C. C., Liu, Y. H., et al., 2018, 3D forward modeling and response analysis for marine CSEM towed by two ships: Applied Geophysics, 15(1), 11-25.
[21]
Zhdanov, M. S., Lee, S. K., and Yoshioka, K., 2006, Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity: Geophysics, 71(6), G333-G345.
2018, Vol. 15 
