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

Yang Jun,Liu Ying,Wu Xiao-Ping

DOI: https://doi.org/10.6038/cjg20150817

2015-01-01

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.