A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling

Puzyrev, Vladimir; Koldan, Jelena; Puente, Josep de la; Houzeaux, Guillaume; Vázquez, Mariano; María, José

doi:10.1093/gji/ggt027

Cited by 136 publications

(65 citation statements)

References 60 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…For example, the performance of several tested approximate inverse preconditioners with static and dynamic sparsity pattern strategies for the systems resulting from the FD discretization was unsatisfactory in some cases. However, for the formulation in potentials, the truncated approximate inverse (TAI) preconditioner from Puzyrev et al (2013) and its variants based on topological and geometric approaches for an a priori chosen sparsity pattern performed reliably. Hence, we use this preconditioner with the block methods in the FE code.…”

Section: Preconditioningmentioning

confidence: 99%

A review of block Krylov subspace methods for multisource electromagnetic modelling

Puzyrev¹,

María

2015

Geophys. J. Int.

Self Cite

View full text Add to dashboard Cite

SUMMARYPractical applications of controlled-source electromagnetic modeling require solutions for multiple sources at several frequencies, thus leading to a dramatic increase of the computational cost. In this paper we present an approach using block Krylov subspace solvers that are iterative methods especially designed for problems with multiple right-hand-sides. Their main advantage is the shared subspace for approximate solutions, hence, these methods are expected to converge in less iterations than the corresponding standard solver applied to each linear system. Block solvers also share the same preconditioner, which is constructed only once. Simultaneously computed block operations have better utilization of cache due to the less frequent access to the system matrix. In this paper we implement two different block solvers for sparse matrices resulting from the finitedifference and the finite-element discretizations, discuss the computational cost of the algorithms and study their dependence on the number of right-hand sides given at once. The effectiveness of the proposed methods is demonstrated on two electromagnetic survey scenarios, including a large marine model. As the results of the simulations show, when a powerful preconditioning is employed, block methods are faster than standard iterative techniques in terms of both iterations and time.

show abstract

Section: Preconditioningmentioning

confidence: 99%

A review of block Krylov subspace methods for multisource electromagnetic modelling

Puzyrev¹,

María

2015

Geophys. J. Int.

Self Cite

View full text Add to dashboard Cite

show abstract

“…Simulations for CSEM synthetic models are performed using the finite-element code by Puzyrev et al (2013). The use of completely unstructured tetrahedral meshes allows us to simulate complicated 3D reservoirs more accurately and with modest computational cost.…”

Section: Design Of the Modelmentioning

confidence: 99%

CSEM Monitoring of a CO2 Reservoir Imaged by MT

Puzyrev¹,

Queralt²,

Ledo³

et al. 2014

Proceedings

View full text Add to dashboard Cite

SUMMARYDuring the last years, different studies based on numerical simulations have shown the potential of CSEM for CO2 monitoring. In this study, we carried out a set of simulations for CSEM monitoring of CO2 realistic deep saline reservoir in a 3D dome anticline structure. We perform numerical simulations in different scenarios (emission frequencies and surface-to-borehole configurations) in order to investigate the effect on the resolution when the simulations are done considering a baseline geoelectric model resulting from the inversion of MT synthetic data.

show abstract

“…Direct matrix solvers for large geophysical 3D forward problems [33,28,30,26,24,14,19,8,35,2,15] can be computationally expensive, both in terms of memory and CPU time. A way to overcome this limitation is by reducing the dimensionality of the problem, namely from 3D modeling to a 2.5D [23], 2D [4] and/or 1D [20] approximations that are only suitable for particular geometries.…”

Section: Introductionmentioning

confidence: 99%

A multi-domain decomposition-based Fourier finite element method for the simulation of 3D marine CSEM measurements

Bakr

Pardo

2017

Comput Geosci

View full text Add to dashboard Cite

We introduce a multi-domain decomposition Fourier finite element (MDDFFE) method for the simulation of three-dimensional (3D) marine controlled source electromagnetic (CSEM) measurements. The method combines a 2D finite element (FE) method in two spatial dimensions with a hybrid discretization based on a Fourier FE method along the third dimension. The method employs a secondary field formulation rather than the full field formulation. By using the secondary field formulation, we avoid to numerically model the source singularities and reduce the effect of the air layer. We apply the MDDFFE method to several synthetic marine CSEM examples exhibiting bathymetry and/or multiple 3D subdomains. Numerical results show that the use of the MDDFFE method reduces the problem size by as much as 87% in terms of the number of unknowns, without any sacrifice in accuracy.

show abstract

A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling

Cited by 136 publications

References 60 publications

A review of block Krylov subspace methods for multisource electromagnetic modelling

A review of block Krylov subspace methods for multisource electromagnetic modelling

CSEM Monitoring of a CO2 Reservoir Imaged by MT

A multi-domain decomposition-based Fourier finite element method for the simulation of 3D marine CSEM measurements

Contact Info

Product

Resources

About