Abstract:Nodal finite element (FE) has become a popular numerical method in computational geo-electromagnetics since 1970 s. However, the zero-divergence property of the electromagnetic field and the discontinuity of the normal component of the electrical field harass the elegant implementation of nodal FE. The former brings in the spurious solution in the simulation while the latter makes it impossible for regular nodal FE to solve electromagnetic field directly. This paper focuses on 3D forward modeling of marine controlled-source electromagnetics (CSEM) with an electrical dipole source, using vector finite element method on unstructured tetrahedral grids which can simulate complicated geometry and submarine topography without spurious solution. The curl-curl electrical field full-wave equation, which is derived from Maxwell's equations, is used as the governing equation in our marine CSEM simulation. Somerfeld radiation operator is chosen as the boundary condition on the truncated boundary. By applying the generalized variational principle to the boundary value problem of the electrical field, the equivalent variational problem is obtained. Unstructured tetrahedral grids are employed for spatial discretization in the 3D calculation domain using a non-commercial mesh generator Gmsh. Nedelec type vector basis functions are adopted for the interpolation of the electrical field within each tetrahedral element. These basis functions are divergence-free intrinsically and they guarantee the continuity of the tangential field component when crossing the interface between two different media. Unlike the nodal FE, the unknowns are assigned to the edges of the elements rather than the nodes. Substituting the electrical field by the interpolation, a system of linear equations is constructed. We solve the system of linear equations with a direct solver. A 1D stratified air-sea-sediment model is designed for the validation of our 3D vector FE algorithm. We output the numerical solution of the electrical field of 0.01 Hz, 0. 1 Hz and 1. 0 Hz with three measurement configurations: inline, broadside and a 45-degree angle between the survey line and the dipole source direction. The comparison between the analytical and numerical solution shows that the relative error are less than 1% for most cases when the observation points are a few hundred meters away from the source. The CSEM response of a 2D air-sea-sediment-basement model with a horst type interface between the sediment and basement layers is then compared with the numerical solution of the well-known 2D FE algorithm MARE2DCSEM and it shows high consistency between the 2D and 3D numerical solutions. The influence of submarine topography to CSEM response is also studied. The simulation indicates that submarine topography has significant impact on the electrical field in marine CSEM. The existence of submarine topography may cover the anomalous signal induced by oil reservoir. Finally, the CSEM response to a 3D reservoir model covered by complicated submarine topography and fluctuant stratum interface is simulated and the results show that the complicated topography makes unsmooth changes to the MVO (Magnitude Versus Offset) curve. We develop a 3D marine CSEM simulation algorithm using vector FE on unstructured meshes. The main advantages of vector FE we use in marine CSEM over regular nodal FE are: (1) no divergence correction is needed because of the vector basis functions are divergence-free; (2) the electrical field is solved directly without worrying about the discontinuity of the normal component of the electrical field. On the other hand, unstructured tetrahedral grids can easily handle complicated submarine topography and refined reservoir structures. The comparisons with 1D analytical solution and 2D numerical solution show high accuracy of our algorithm. The simulations also show that submarine topography may be an important contributor to the uncertainties of CSEM MVO anomalies.

Efficient Simulation of Marine Controlled Source Electromagnetic Responses for Axisymmetric Scatter by Using Numerical Mode Matching Approach

An Efficient Marching-on-in-Degree Solution of Transient Multiscale EM Scattering Problems

Hybrid Time-Domain Adaptive Integration Method for Simulating Transient Electromagnetic Responses of an Arbitrary Structure with Large Scale

3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method

3-D Marine Controlled-Source Electromagnetic Modeling in Electrically Anisotropic Formations Using Scattered Scalar–Vector Potentials

3-D Marine CSEM Forward Modeling in Anisotropic Media Using a Space-Wavenumber Domain Integral Equation Method

Special Section — Marine Controlled-Source Electromagnetic Methods 2D marine controlled-source electromagnetic modeling: Part 1—A nadaptive finite-element algorithm

Finite Element Simulation of 3-D Marine Controlled Source Electromagnetic Fields in Anisotropic Media with Unstructured Tetrahedral Grids

Study on 1d Forward Modeling of Marine Csem Survey

3D simulation of marine CSEM using vector finite element method on unstructured grids

2D Marine Controlled-Source Electromagnetic Modeling: Part 1 — an Adaptive Finite-Element Algorithm

Sensitivity Study of Three-Dimensional Marine Controlled-Source Electromagnetic Method

MCSEM 3D modeling for arbitrarily anisotropic media

3D Finite Element Modeling of Marine Controlled-Source Electromagnetic Fields Using Locally Refined Unstructured Meshes

Simulation and analysis on the prospecting capability of marine controlled-source electromagnetic methods to hydrocarbon reservoirs

Controlled-source Electromagnetic Modeling Using a Versatile Secondary Electric-Field Formulation and Efficient Multigrid-Based Preconditioner

Goal-Oriented 3-D Time-Domain Marine CSEM Modeling With Anisotropy and Topography

3-D Marine CSEM Modeling in General Anisotropic Media by Using an Adaptive Finite Element Approach Based on the Vector-Scalar Potential

Three-dimensional Modeling of Marine Controlled Source Electromagnetic Field Using a Space-Wavenumber Domain Method

3-D finite-volume forward modeling of wide-field EM using scattered potentials

Analyzing the Effect of Dipping Anisotropy on 3D Marine Controlled-Source Electromagnetic Response Based on a Goal-Oriented Adaptive Finite Element Method