Parallelized 3D CSEM modeling using edge-based finite element with total field formulation and unstructured mesh

Cai, Hongzhu; Hu, Xiangyun; Li, Jianhui; Endo, Mitsuru; Xiong, Bin

doi:10.1016/j.cageo.2016.11.009

Cited by 78 publications

(27 citation statements)

References 55 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…; Cai et al . ; Dunham, Ansari and Farquharson ; Zeng et al . ), as well as full transmitting‐current waveforms (Sun et al .…”

Section: Introductionmentioning

confidence: 99%

A finite‐element time‐domain forward‐modelling algorithm for transient electromagnetics excited by grounded‐wire sources

Cai

et al. 2020

Geophysical Prospecting

Self Cite

View full text Add to dashboard Cite

A B S T R A C TA three-dimensional finite-element time-domain forward-modelling algorithm is developed to simulate transient electromagnetics excited by grounded-wire sources. The main advantage of this finite-element time-domain algorithm is that full transmitting-current waveforms and complex-shaped sources resulting from topography can be directly dealt with in this algorithm. The models used to test this algorithm include a homogeneous half-space model, a stratified-medium model, the model of a complex conductor at a vertical contact and the Ovoid Zone massive sulfide deposit at Voisey's Bay, Canada. The homogeneous half-space model is used to determine the truncation boundary for a computational domain, and to compare with the electromagnetic responses excited by step-off, step-on and direct current waveforms. For the stratified-medium model, results demonstrate that full transmitting waveforms have strong effects on the observed electromagnetic responses. The model of a complex conductor at a vertical contact is designed for the grounded electrical source airborne transient electromagnetic method and is also used to examine the effectiveness of the broadside and inline configurations for such a vertical, thin plate embedded in the subsurface. The area of the Ovoid Zone massive sulfide deposit possesses non-negligible topography, the effects of which on the shapes of the grounded-wire sources must be taken into account when implementing the finite-element time-domain solver. The results show that both the broadside and inline electromagnetic responses are strongly affected by the massive conductive ore body.

show abstract

“…; Cai et al . ; Dunham, Ansari and Farquharson ; Zeng et al . ), as well as full transmitting‐current waveforms (Sun et al .…”

Section: Introductionmentioning

confidence: 99%

A finite‐element time‐domain forward‐modelling algorithm for transient electromagnetics excited by grounded‐wire sources

Cai

et al. 2020

Geophysical Prospecting

Self Cite

View full text Add to dashboard Cite

show abstract

“…Thus, both preconditioned iterative solvers and sparse direct solvers are commonly considered. Although the use of unstructured tetrahedral grids [17] or hexahedral grids with hanging nodes [18] reduces the size of the discrete problem, an efficient solver still will be required.…”

Section: Introductionmentioning

confidence: 99%

Parallel Simulation of Audio- and Radio-Magnetotelluric Data

Yavich

Malovichko

Shlykov³

2019

Minerals

View full text Add to dashboard Cite

This paper presents a novel numerical method for simulation controlled-source audio-magnetotellurics (CSAMT) and radio-magnetotellurics (CSRMT) data. These methods are widely used in mineral exploration. Interpretation of the CSAMT and CSRMT data collected over an area with the complex geology requires application of effective methods of numerical modeling capable to represent the geoelectrical model of a deposit well. In this paper, we considered an approach to 3D electromagnetic (EM) modeling based on new types of preconditioned iterative solvers for finite-difference (FD) EM simulation. The first preconditioner used fast direct inversion of the layered Earth FD matrix (Green’s function preconditioner). The other combined the first with a contraction operator transformation. To illustrate the effectiveness of the developed numerical modeling methods, a 3D resistivity model of Aleksandrovka study area in Kaluga Region, Russia, was prepared based on drilling data, AMT, and a detailed CSRMT survey. We conducted parallel EM simulation of the full CSRMT survey. Our results indicated that the developed methods can be effectively used for modeling EM responses over a realistic complex geoelectrical model for a controlled source EM survey with hundreds of receiver stations. The contraction-operator preconditioner outperformed the Green’s function preconditioner by factor of 7–10, both with respect to run-time and iteration count, and even more at higher frequencies.

show abstract

“…They have been frequently applied for the modeling of 3D structures from EM data (e.g. Ahagon and Kameari 2017;Ahagon et al 2018;Cai et al 2017;Camargos et al 2017;Liu et al 2010Liu et al , 2018Puzyrev et al 2013;Ren et al 2013;Schwarzbach et al 2011;Weiss 2017). On the other hand, hexahedron mesh geometries are also commonly applied (Cai et al 2014;Farquharson and Miensopust 2011;Mitsuhata and Uchida 2004;Mogi 1996;Nam et al 2007;Shi et al 2004;Xiao et al 2018;Xiong et al 2018;Zyserman and Santos 2000).…”

Section: Introductionmentioning

confidence: 99%

“…it is advantageous to extend the computational capabilities of modeling algorithms toward the handling of full 3D anisotropic structures. Recent implementations of anisotropy in VFEM have been applied with either linear-order tetrahedral elements (Cai et al 2017;Liu et al 2018) or using linear-order hexahedral elements (Xiao et al 2018).…”

Section: Introductionmentioning

confidence: 99%

Multi-order vector finite element modeling of 3D magnetotelluric data including complex geometry and anisotropy

Rivera-Rios

Zhou

Heinson

et al. 2019

Earth Planets Space

View full text Add to dashboard Cite

We introduce MoVFEM, a computational algorithm for the modeling of three-dimensional magnetotelluric (MT) data using a vector finite element method of specific order from multiple elements' orders. Our algorithm allows complex geometries, topography, and anisotropic resistivity structures. The software calculates secondary electric and magnetic fields for a plane-wave primary magnetic field. Accurate calculation of fields in the boundary regions of the computational domain are ensured by the implementation of the Generalized Perfect Matched Layers method. We validate the MoVFEM algorithm by applications to various scenarios, which allow a comparison with analytical or accepted numerical solutions where available. The respective results of our algorithm are in good agreement with existing solutions. which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

show abstract

Parallelized 3D CSEM modeling using edge-based finite element with total field formulation and unstructured mesh

Cited by 78 publications

References 55 publications

A finite‐element time‐domain forward‐modelling algorithm for transient electromagnetics excited by grounded‐wire sources

A finite‐element time‐domain forward‐modelling algorithm for transient electromagnetics excited by grounded‐wire sources

Parallel Simulation of Audio- and Radio-Magnetotelluric Data

Multi-order vector finite element modeling of 3D magnetotelluric data including complex geometry and anisotropy

Contact Info

Product

Resources

About