1887
Volume 62, Issue 6
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

A fully three‐dimensional finite‐element algorithm has been developed for simulating controlled‐source electromagnetic surveys. To exploit the advantages of geometric flexibility, frequency‐domain Maxwell's equations of the secondary electric field were discretised using edge‐based finite elements while the primary field was calculated analytically for a horizontally layered‐earth model. The resulting system of equations for the secondary field was solved using a parallel version of direct solvers. The accuracy of the algorithm was successfully verified by comparisons with integral‐equations and iterative solutions, and the applicability to models containing large conductivity contrasts was verified against published data. The advantages of geometry‐conforming meshes have been demonstrated by comparing different mesh systems to simulate an inclined sheet model. A comparison of the performance between direct and iterative solvers demonstrated the superior efficiency of direct solvers, particularly for multisource problems.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12132
2014-06-05
2020-04-06
Loading full text...

Full text loading...

References

  1. AmestoyP.R., GuermoucheA., L'ExcellentJ.‐Y. and PraletS.2006. Hybrid scheduling for the parallel solution of linear systems. Parallel Computing32, 136–156.
    [Google Scholar]
  2. AvdeevD.B., KuvshinovA.V., PankratovO.V. and NewmanG.A.2002. Three‐dimensional induction logging problems, Part I: An integral equation solution and model comparisons. Geophysics67, 413–426.
    [Google Scholar]
  3. BadeaE.A., EverettM.E., NewmanG.A. and BiroO.2001. Finite‐element analysis of controlled‐source electromagnetic induction using Coulomb‐gauged potentials. Geophysics66, 786–799.
    [Google Scholar]
  4. da SilvaN.V., MorganJ.V., MacGregorL. and WarnerM.2012. A finite element multifrontal method for 3D CSEM modeling in the frequency domain. Geophysics77, E101–E115.
    [Google Scholar]
  5. ErnstO.G. and GanderM.J.2012. Why it is Difficult to Solve Helmholtz Problems with Classical Iterative Methods. Numerical Analysis of Multiscale Problems83, 325–363.
    [Google Scholar]
  6. FarquharsonC.G., DuckworthK. and OldenburgD.W.2006. Comparison of integral equation and physical scale modeling of the electromagnetic responses of models with large conductivity contrasts. Geophysics71, G169–G177.
    [Google Scholar]
  7. FarquharsonC.G. and MiensopustM.P.2011. Three‐dimensional finite‐element modelling of magnetotelluric data with a divergence correction. Journal of Applied Geophysics75, 699–710.
    [Google Scholar]
  8. GrayverA. and StreichR.2012. Comparison of iterative and direct solvers for 3D CSEM modeling. 82nd SEG Annual Meeting, Expanded Abstracts.
  9. GribenkoA. and ZhdanovM.2007. Rigorous 3D inversion of marine CSEM data based on the integral equation method. Geophysics72, WA73–WA84.
    [Google Scholar]
  10. GuptaA., 2000, WSMP: Watson Sparse Matrix Package, Part I – direct solution of symmetric sparse systems.
  11. HaberE., AscherU.M., AruliahD.A. and OldenburgD.W.2000. Fast Simulation of 3D Electromagnetic Problems Using Potentials. Journal of Computational Physics163, 150–171.
    [Google Scholar]
  12. Intel
    Intel , 2009, Intel Math Kernel Library Reference Manual: Intel Corporation.
  13. JinJ.2002. The finite element method in electromagnetics, 2nd ed. John Wiley & Sons.
    [Google Scholar]
  14. KangS., SeolS.J., ChungY. and KwonH.‐S.2013. Pitfalls of 1D inversion of small‐loop electromagnetic data for detecting man‐made objects. Journal of Applied Geophysics90, 96–109.
    [Google Scholar]
  15. KeyK. and OvallJ.2011. A parallel goal‐oriented adaptive finite element method for 2.5‐D electromagnetic modelling. Geophysical Journal International186, 137–154.
    [Google Scholar]
  16. KimH.J., SongY. and LeeK.H.1997. High‐frequency electromagnetic inversion for a dispersive layered earth. Journal of geomagnetism and geoelectricity49, 1439–1450.
    [Google Scholar]
  17. MitsuhataY. and UchidaT.2004. 3D magnetotelluric modeling using the T‐Ω finite‐element method. Geophysics69, 108–119.
    [Google Scholar]
  18. MiyataK.2006. Magnetic field analysis by the edge element FEM. Electrical Engineering in Japan154, 59–65.
    [Google Scholar]
  19. MukherjeeS. and EverettM.2011. 3D controlled‐source electromagnetic edge‐based finite element modeling of conductive and permeable heterogeneities. Geophysics76, F215–F226.
    [Google Scholar]
  20. NamM.J., KimH.J., SongY., LeeT.J., SonJ.‐S. and SuhJ.H.2007. 3D magnetotelluric modelling including surface topography. Geophysical Prospecting55, 277–287.
    [Google Scholar]
  21. NewmanG.A. and AlumbaughD.L.2002. Three‐dimensional induction logging problems, Part 2: A finite‐difference solution. Geophysics67, 484–491.
    [Google Scholar]
  22. PuzyrevV., KoldanJ., de la PuenteJ., HouzeauxG., VázquezM. and CelaJ.M.2013. A parallel finite‐element method for three‐dimensional controlled‐source electromagnetic forward modelling. Geophysical Journal International193, 678–693.
    [Google Scholar]
  23. RenZ., KalscheuerT., GreenhalghS. and MaurerH.2013. A goal‐oriented adaptive finite‐element approach for plane wave 3‐D electromagnetic modelling. Geophysical Journal International194, 700–718.
    [Google Scholar]
  24. SaadY.2003. Iterative methods for sparse linear systems. Society for Industrial Mathematics.
  25. SasakiY.2001. Full 3‐D inversion of electromagnetic data on PC. Journal of Applied Geophysics46, 45–54.
    [Google Scholar]
  26. SasakiY.2004. Three‐dimensional inversion of static‐shifted magnetotelluric data. Earth, Planets and Space56, 239–248.
    [Google Scholar]
  27. SasakiY. and MejuM.A.2006. A multidimensional horizontal‐loop controlled‐source electromagnetic inversion method and its use to characterize heterogeneity in aquiferous fractured crystalline rocks. Geophysical Journal International166, 59–66.
    [Google Scholar]
  28. SchenkO. and GärtnerK.2004. Solving unsymmetric sparse systems of linear equations with PARDISO. Future Generation Computer Systems20, 475–487.
    [Google Scholar]
  29. SchenkO., GärtnerK., FichtnerW. and StrickerA.2001. PARDISO: a high‐performance serial and parallel sparse linear solver in semiconductor device simulation. Future Generation Computer Systems18, 69–78.
    [Google Scholar]
  30. SchwarzbachC., BörnerR.‐U. and SpitzerK.2011. Three‐dimensional adaptive higher order finite element simulation for geo‐electromagnetics—a marine CSEM example. Geophysical Journal International187, 63–74.
    [Google Scholar]
  31. SmithJ.T.1996. Conservative modeling of 3‐D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator. Geophysics61, 1319–1324.
    [Google Scholar]
  32. SongY., KimH.J. and LeeK.H.2002. An integral equation representation of wide‐band electromagnetic scattering by thin sheets. Geophysics67, 746–754.
    [Google Scholar]
  33. StreichR.2009. 3D finite‐difference frequency‐domain modeling of controlled‐source electromagnetic data: Direct solution and optimization for high accuracy. Geophysics74, F95–F105.
    [Google Scholar]
  34. StreichR., SchwarzbachC., BeckenM. and SpitzerK.2010. Controlled‐source electromagnetic modeling studies: utility of auxiliary potentials for low frequency stabilization. 72nd EAGE Conference & Exhibition, Expanded Abstracts.
  35. TongX.‐Z., LiuJ.‐X., XieW., XuL.‐H., GuoR.‐W. and ChengY.‐T.2009. Three‐dimensional forward modeling for magnetotelluric sounding by finite element method. Journal of Central South University of Technology16, 136–142.
    [Google Scholar]
  36. UmE.S., CommerM. and NewmanG.A.2013. Efficient pre‐conditioned iterative solution strategies for the electromagnetic diffusion in the Earth: finite‐element frequency‐domain approach. Geophysical Journal International193, 1460–1473.
    [Google Scholar]
  37. UmE.S., HarrisJ.M. and AlumbaughD.L.2010. 3D time‐domain simulation of electromagnetic diffusion phenomena: A finite‐element electric‐field approach. Geophysics75, F115–F126.
    [Google Scholar]
  38. WeissC.J. and ConstableS.2006. Mapping thin resistors and hydrocarbons with marine EM methods, Part II — Modeling and analysis in 3D. Geophysics71, G321–G332.
    [Google Scholar]
  39. WeissC.J. and NewmanG.A.2002. Electromagnetic induction in a fully 3‐D anisotropic earth. Geophysics67, 1104–1114.
    [Google Scholar]
  40. WonI.J.2003. Small frequency‐domain electromagnetic induction sensors. The Leading Edge22, 320–322.
    [Google Scholar]
  41. XiongZ.1992. Electromagnetic modeling of 3‐D structures by the method of system iteration using integral equations. Geophysics57, 1556–1561.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12132
Loading
/content/journals/10.1111/1365-2478.12132
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): 3D , CSEM , Direct solver , FEM and Modelling
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error