Finite element modelling of marine controlled-source electromagnetic responses in two-dimensional dipping anisotropic conductivity structures

SUMMARY We present an adaptive finite element (FE) solution for the marine controlled-source electromagnetic (CSEM) forward problem in 2-D dipping anisotropic conductivity structures. Our code is implemented on an unstructured triangular mesh, which allows for arbitrary model geometries including bathymetry and dipping layers. We have verified the FE code using a layered 1-D model with anisotropy. For this model, the FE algorithm provides accurate results. The relative amplitude error between the analytical solution and numerical results for all electromagnetic components except Hz for large offsets is less than 1 per cent and the error in phase is less than 0.5°. We simulate the marine CSEM responses of a 2-D anisotropic model, and the numerical results show that the CSEM fields are affected considerably by anisotropy in both the reservoir target and the surrounding sediment, but to different extent.

[1]  Chuan Yi Tang,et al.  A 2.|E|-Bit Distributed Algorithm for the Directed Euler Trail Problem , 1993, Inf. Process. Lett..

[2]  Roland W. Freund,et al.  An Implementation of the QMR Method Based on Coupled Two-Term Recurrences , 1994, SIAM J. Sci. Comput..

[3]  James Demmel,et al.  A Supernodal Approach to Sparse Partial Pivoting , 1999, SIAM J. Matrix Anal. Appl..

[4]  Chester J. Weiss,et al.  A matrix‐free approach to solving the fully 3D electromagnetic induction problem , 2001 .

[5]  Lucy MacGregor,et al.  Sea Bed Logging (SBL), a new method for remote and direct identification of hydrocarbon filled layers in deepwater areas , 2002 .

[6]  Yuguo Li.  A finite-element algorithm for electromagnetic induction in two-dimensional anisotropic conductivity structures , 2002 .

[7]  Jonathan Richard Shewchuk,et al.  Delaunay refinement algorithms for triangular mesh generation , 2002, Comput. Geom..

[8]  L. Lines,et al.  Influence of seismic anisotropy on prestack depth migration , 2004 .

[9]  Chester J. Weiss,et al.  Adaptive finite-element modeling using unstructured grids: The 2D magnetotelluric example , 2005 .

[10]  Michael J. Tompkins,et al.  The Role of Vertical Anisotropy In Interpreting Marine Controlled-Source Electromagnetic Data , 2005 .

[11]  Jeffrey S. Ovall.  Asymptotically exact functional error estimators based on superconvergent gradient recovery , 2006, Numerische Mathematik.

[12]  Kerry Key,et al.  First results from a marine controlled‐source electromagnetic survey to detect gas hydrates offshore Oregon , 2006 .

[13]  B. Ursin,et al.  Electromagnetic fields in planarly layered anisotropic media , 2007 .

[14]  Kerry Key,et al.  2D marine controlled-source electromagnetic modeling: Part 1 — An adaptive finite-element algorithm , 2007 .

[15]  Tage Røsten,et al.  A 2.5D finite-element-modeling difference method for marine CSEM modeling in stratified anisotropic media , 2008 .

[16]  Josef Pek,et al.  Adaptive finite element modelling of two-dimensional magnetotelluric fields in general anisotropic media , 2008 .