3D MCSEM modeling based on domain-decomposition finite-element method

摘要

Domain decomposition approach is effective in converting large modeling problems into many small ones, so that the parallel computation can be easily facilitated and the memory consumption can be largely reduced. In this paper, we developed a domain-decomposition method for 3D time-domain marine controlled-source (MCSEM) modeling based on Dual-Primal Finite-Element Tearing and Interconnecting technique (FETI-DP). We use tetrahedral grids to do the spatial discretization first and then partition the mesh into subdomains. In each subdomain, the finite-element equations are independently constructed using the first-order vector shape functions and boundary conditions. To build the interconnections among these subdomains, we establish the interface equation for the Lagrange multipliers that are used as the boundary conditions for the subdomains. After solving the equations, the electric field for all subdomains can be obtained from the Lagrange multipliers. For the time discretization, we use unconditionally stable backward Euler scheme. We validate our algorithm by comparing with a semi-analytical solutions for a uniform half-space model under the ocean. Numerical experiments show that compared to the conventional finite-element method, our domain-decomposed method can save large amount of memory, render it possible to run 3D MCSEM modeling in small computational equipment.

关键词