| 1. |
- Buntin, Laura M.
(författare)
-
Efficient Electromagnetic Induction Modelling : Adaptive mesh optimisation, advanced boundary methods and iterative solution techniques
- 2023
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Forward modelling of electromagnetic induction data simulates the electric and magnetic fields within a computational domain for a given distribution of electromagnetic material properties and a given source of the electromagnetic field. The quantities of interest are the fields at receiver locations at the Earth's surface. Reliable results require high accuracy solutions at the receivers. First and foremost, numerical computations need to be accurate, but ideally they are also resource efficient, i.e., as fast and cheap as possible. Run time and memory demand mainly depend on the size of the numerical problem to be solved. This thesis addresses specific steps within the forward modelling procedure of electromagnetic induction data in order to improve the solution accuracy of forward modelling as well as to reduce computational resources. The solution accuracy is strongly influenced by the spatial discretisation of the computational domain, which directly correlates with the numerical problem size. To optimise the solution accuracy while keeping the numerical problem size as small as possible, a goal-oriented adaptive mesh refinement scheme for three-dimensional controlled-source electromagnetic models is developed. In addition, this thesis investigates the influence of different types of boundary methods on the solution accuracy at the receivers. To replace inhomogeneous boundary conditions in magnetotelluric total-field modelling by perfectly-matched layers (PML), a domain decomposition approach (the total and scattered field decomposition) is adapted for Earth models. By reducing boundary effects, the approach yields superior solution accuracy for specific types of models. The fastest and most memory-efficient way to solve large numerical problems are iterative solution methods. Iterative solvers, however, work poorly for numerical systems arising from domains bounded by PML. In this thesis a preconditioned iterative solution framework that efficiently solves PML-bounded magnetotelluric models is proposed.
|
|
| 2. |
- Buntin, Laura M., et al.
(författare)
-
Improved accuracy of plane-wave electromagnetic modelling by application of the total and scattered field decomposition and perfectly matched layers
- 2023
-
Ingår i: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 235:2, s. 1201-1217
-
Tidskriftsartikel (refereegranskat)abstract
- In 2-D magnetotelluric modelling, the standard application of Dirichlet boundary conditions (BC) may severely diminish the solution accuracy, because the unknown scattered part of the electromagnetic field is erroneously reflected at the domain boundary. Therefore, we adapt the total and scattered field decomposition (TSFD) to geophysical modelling, enabling the application of fully absorbing boundary methods, here perfectly matched layers (PML), to the scattered field. Our novel TSFD divides the modelling domain into two regions. In the total-field region containing the area of interest, the solution is computed for the total field. In the scattered-field region containing the boundaries, the solution is obtained for the scattered field, which is fully absorbed by PML at the boundaries. The plane-wave source is excited at the TSFD interface between both regions. Thus, boundary reflections are eradicated leading to superior solution accuracy, and boundaries can be placed closer to the receivers, shrinking the computational problem. Especially for challenging models with strong lateral changes, the solution accuracy of the TSFD is superior to that of the standard Dirichlet approach. Owing to the linearity of Maxwell's equations, the inaccuracy introduced to the electric and magnetic fields by using Dirichlet BC can be expected to partly cancel out in the magnetotelluric transfer functions, for example the impedance tensor. In this work, we quantify this cancellation effect. The inaccuracy is less than typical measurement errors in the vast majority of apparent resistivity and phase data, even, when the primary fields are strongly inaccurate.
|
|
| 3. |
|
|
| 4. |
- Rulff, Paula, et al.
(författare)
-
Efficient goal-oriented mesh refinement in 3-D finite-element modelling adapted for controlled source electromagnetic surveys
- 2021
-
Ingår i: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 227:3, s. 1624-1645
-
Tidskriftsartikel (refereegranskat)abstract
- We developed a 3-D forward modelling code, which simulates controlled source electromagnetic problems in frequency domain using edge-based finite elements and a total electric field approach. To evaluate electromagnetic data acquired across complex subsurface structures, software performing accurate 3-D modelling is required, especially for incorporation in inversion approaches. Our modelling code aims at finding a good compromise between the necessary solution accuracy at the points of interest and the general problem size by using a goal-oriented mesh refinement strategy designed for models of variable electric conductivity and magnetic permeability. To formulate an improved error estimator suitable for controlled source electromagnetic problems, we developed literature approaches of mesh refinement further targeting three aspects. First, to generate a roughly homogeneously fine mesh discretization around all receiver sites, our new error estimator weights the adjoint source term by the approximate decay of the electric field with increasing distance from the primal source using the expression for a homogeneous half-space. This causes almost no additional computational cost. Second, the error estimator employed in the refinement approach can be optimized for models with pronounced conductivity and magnetic permeability contrasts as often encountered in, for example, mineral prospecting scenarios by optionally including terms that measure the continuity of the normal component of current flow and the tangential component of the magnetic field across interfaces of abutting elements. Third, to avoid amplitude-dependent over-refining of the mesh, we formulate our element-wise error estimators relative to the local amplitude of the electromagnetic field. In this work, we evaluate the implemented adaptive mesh refinement approach and its solution accuracy comparing our solutions for simple 1-D models and a model with 3-D anomalies to semi-analytic 1-D solutions and a second-order finite-element code, respectively. Furthermore, a feasibility study for controlled-source electromagnetic measurements across ferrous mineral deposits is conducted. The numerical experiments demonstrate that our new refinement procedure generates problem-specific finite-element meshes and yields accurate solutions for both simple synthetic models and realistic survey scenarios. Especially for the latter, characteristics of our code, such as the possibility of modelling extended sources as well as including arbitrary receiver distributions and detailed subsurface anomalies, are beneficial.
|
|
