# CSEM forward modelling & Edge elements formulation¶

The last decade has been a period of rapid growth for electromagnetic methods (EM) in geophysics, mostly because of their industrial adoption. In particular, the 3D marine controlled-source electromagnetic (3D CSEM) method has become an important technique for reducing ambiguities in data interpretation in hydrocarbon exploration. In order to be able to predict the EM signature of a given geological structure, modelling tools provide us with synthetic results which we can then compare to measured data. In particular, if the geology is structurally complex, one might need to use methods able to cope with such complexity in a natural way by means of, e.g., an unstructured mesh representing its geometry. Among the modelling methods for EM based upon 3D unstructured meshes, the Nédélec Edge Finite Element Method (EFEM) offers a good trade-off between accuracy and number of degrees of freedom, e.g. size of the problem. Furthermore, its divergence-free basis is very well suited for solving Maxwell’s equation. On top of that, we choose to support tetrahedral meshes, as these are the easiest to use for very large domains or complex geometries. We present the numerical formulation and results of 3D CSEM forward modelling (FM) using tetrahedral EFEM on unstructured meshes.

## CSEM problem¶

3D CSEM FM is typically solved in frequency domain, which involves the numerical solution of Maxwell’s equations in stationary regimes for heterogeneous anisotropic electrically conductive domains. Furthermore, CSEM surveys generally work with low frequency electromagnetic fields ( 1 Hz to 3 Hz) because the electric conductivity of the geological structures is much larger than their dielectric permittivity. As a consequence, in an unbound domain , the electric field can be obtained by solving Maxwell’s equations in their diffusive form:

(1)¶

where the harmonic time dependence is omitted, with is the angular frequency, the free space magnetic permeability, the distribution of source current, the induced current in the conductive Earth and the electrical conductivity which is assumed isotropic for simplicity.

In numerical approximations of EM fields there are two main drawbacks. The first one is the inevitable spatial singularity at the source. The second is the grid refinement requirements in order to capture the rapid change of the primary field [1]. In order to mitigate these issues, PETGEM used a secondary field approach where the total electric field is obtained as:

(2)¶

where subscripts and represent a primary field and secondary field respectively. For a general layered Earth model, can be computed semi-analytically by using Hankel transform filters. Based on this decomposition and following the work by [3] the equation system to solve is:

(3)¶

where the electrical conductivity is a function of position that is allowed to vary in 3D, whereas the vacuum permeability is set to the free space value . Homogeneous Dirichlet boundary conditions, on , are applied. The range of applicability of this conditions can be determined based on the skin depth of the electric field [4].

## Edge finite element formulation¶

For the computation of , PETGEM implemented the Nédélec EFEM which uses vector basis functions defined on the edges of the corresponding elements. Its basis functions are divergence-free but not curl-free [2]. Thus, EFEM naturally ensures tangential continuity and allows normal discontinuity of at material interfaces. PETGEM used unstructured tetrahedral meshes because of their ability to represent complex geological structures such as bathymetry or reservoirs as well as the local refinement capability in order to improve the solution accuracy. Figure 4.1 shows the tetrahedral Nédélec elements (lowest order) together with their node and edge indexing.

In PETGEM workflow, the tangential component of the secondary electric field is assigned to the edges in the mesh. Therefore, all components of the electric field at a point located inside a tetrahedral element can be obtained as follows:

(4)¶

where are the vector basis functions associated to each edge and their degrees of freedom. Considering the node and edge indexing in Figure 4.1, the vector basis functions can be expressed as follows:

(5)¶

where subscripts and are the first and second nodes linked to the -th edge, are the linear nodal basis functions, and is the length of the -th edge of the element .

By substituting equation (4) into (3), and using Galerkin’s approach, the weak form of the original differential equation becomes:

(6)¶

The compact discretized form of (6) is obtained after applying the Green’s theorem:

(7)¶

where and are the elemental stiffness and mass matrices which can be calculated analytically or numerically [2], and is the right hand side which requires numerical integration.

[1] | Cai, H., Xiong, B., Han, M. and Zhdanov, M. (2014). 3D controlled-source electromagnetic modeling in anisotropic medium using edge-based finite element method. Computers & Geosciences, 73, 164–176. |

[2] | (1, 2) Jin, J. (2002). The Finite Element Method in Electromagnetics. Wiley, New York, second edn. |

[3] | Newman, G.A. and Alumbaugh, D.L. (2002). Three-dimensional induction logging problems, Part 2: A finite difference solution. Geophysics, 67(2), 484–491. |

[4] | Puzyrev, V., Koldan, J., de la Puente, J., Houzeaux, G., Vázquez, M. and Cela, J.M. (2013). A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling. Geophysical Journal International, ggt027. |