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.

2D marine controlled-source electromagnetic modeling: Part 2 — The effect of bathymetry

Special Section — Marine Control-Source Electromagnetic Methods 2D Marine Controlled-Source Electromagnetic Modeling: Part 2 — the Effect of Bathymetry

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

Bathymetric analysis and corrections for 3-D marine controlled-source electromagnetic field

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

Finite Element Modelling of Marine Controlled-Source Electromagnetic Responses in Two-Dimensional Dipping Anisotropic Conductivity Structures

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

Finite-Element Scheme With Total-Field for 3-D Marine Controlled-Source Electromagnetic Modeling With Bathymetry and Multisources

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

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

3-D Adaptive Finite-Element Modeling of Marine Controlled-Source Electromagnetics with Seafloor Topography Based on Secondary Potentials

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

2.5D Marine CSEM Modeling in the Frequency-Domain Based on an Improved Interpolation Scheme at Receiver Positions

Effects of the electric anisotropy on marine controlled-source electromagnetic responses

2.5D finite element modeling of marine controlled-source electromagnetic fields

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

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

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

Impacts of Complex Geological Structures on Marine Controlled-Source Electromagnetic Responses