P. Rulff
info
Please Note
<p>This page displays the records of the person named above and is not linked to a unique person identifier. This record may need to be merged to a profile.</p>
14 records found
1
Computational electromagnetic geophysics for groundwater system studies
A review of established practices and recent advances
Journal article
(2026)
-
Paula Rulff, Octavio Castillo-Reyes, Wouter Deleersnyder, Maria Carrizo Mascarell, Burke J. Minsley, Jude King
Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies.As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.
...
Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies.As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.
Journal article
(2026)
-
David Bruhn, Hemmo A. Abels, Patrick Fulton, Virginie Harcouët-Menou Harcouët-Menou, Ernst Huenges, Stefan Jansen, Alexis Koulidis, Susanne Laumann, Haiyan Lei, Joseph Moore, Paula Rulff, Thorben Schöfisch, Auke Barnhoorn, Evert Slob, Philip J. Vardon, Liliana Vargas Meleza, Denis Voskov, Claire Bossennec, Aoife K. Braiden, Maren Brehme, Romain Chassagne, Alexandros Daniilidis, Mathieu Darnet, Guy Drijkoningen
Low-enthalpy geothermal heat production is becoming increasingly common, which leads to the potentially competitive use of the available subsurface space, especially in densely populated urban areas. A specific challenge presented by the high density of different geothermal systems is understanding the details of convective and conductive heat flow processes and detailed monitoring of properties and processes in the subsurface.
On the TU Delft campus, we aim to drill a borehole of around 4.5 km depth to be used for the exploration, observation, and monitoring of subsurface processes that will be part of a larger research infrastructure under development. This so-called urban energy laboratory includes – in addition to the deep multi-use borehole – a well-instrumented geothermal doublet drilled in 2023, reaching to a depth of 2.2 km; a local seismic monitoring system (installed in 2022); an ultra-sensitive portable seismic monitoring array; and a high-temperature aquifer heat storage system (HT-ATES), for which a pilot well was drilled in 2024. With this urban energy laboratory, we want to tackle problems and better understand processes related to multiple and/or competing subsurface uses in urban environments. The deep exploration and monitoring borehole is designed specifically to monitor fluid and/or flux movement in 3D with unprecedented precision, aiming to understand the propagation of the geothermal cold front and reservoir pressures.
During the 3 d International Continental Scientific Drilling Program (ICDP)-sponsored UrbEnLab workshop, 75 scientists from 17 countries met in Delft, the Netherlands, in June 2024 to prioritize the scientific ambitions of the deep exploration and monitoring borehole and to discuss potential techniques that could be applied to tackle them. Assessing the life cycle of a geothermal system situated in a complex heterogeneous sedimentary system was defined as the broad aim, with revealing the detailed flow field established being a key priority. ...
On the TU Delft campus, we aim to drill a borehole of around 4.5 km depth to be used for the exploration, observation, and monitoring of subsurface processes that will be part of a larger research infrastructure under development. This so-called urban energy laboratory includes – in addition to the deep multi-use borehole – a well-instrumented geothermal doublet drilled in 2023, reaching to a depth of 2.2 km; a local seismic monitoring system (installed in 2022); an ultra-sensitive portable seismic monitoring array; and a high-temperature aquifer heat storage system (HT-ATES), for which a pilot well was drilled in 2024. With this urban energy laboratory, we want to tackle problems and better understand processes related to multiple and/or competing subsurface uses in urban environments. The deep exploration and monitoring borehole is designed specifically to monitor fluid and/or flux movement in 3D with unprecedented precision, aiming to understand the propagation of the geothermal cold front and reservoir pressures.
During the 3 d International Continental Scientific Drilling Program (ICDP)-sponsored UrbEnLab workshop, 75 scientists from 17 countries met in Delft, the Netherlands, in June 2024 to prioritize the scientific ambitions of the deep exploration and monitoring borehole and to discuss potential techniques that could be applied to tackle them. Assessing the life cycle of a geothermal system situated in a complex heterogeneous sedimentary system was defined as the broad aim, with revealing the detailed flow field established being a key priority. ...
Low-enthalpy geothermal heat production is becoming increasingly common, which leads to the potentially competitive use of the available subsurface space, especially in densely populated urban areas. A specific challenge presented by the high density of different geothermal systems is understanding the details of convective and conductive heat flow processes and detailed monitoring of properties and processes in the subsurface.
On the TU Delft campus, we aim to drill a borehole of around 4.5 km depth to be used for the exploration, observation, and monitoring of subsurface processes that will be part of a larger research infrastructure under development. This so-called urban energy laboratory includes – in addition to the deep multi-use borehole – a well-instrumented geothermal doublet drilled in 2023, reaching to a depth of 2.2 km; a local seismic monitoring system (installed in 2022); an ultra-sensitive portable seismic monitoring array; and a high-temperature aquifer heat storage system (HT-ATES), for which a pilot well was drilled in 2024. With this urban energy laboratory, we want to tackle problems and better understand processes related to multiple and/or competing subsurface uses in urban environments. The deep exploration and monitoring borehole is designed specifically to monitor fluid and/or flux movement in 3D with unprecedented precision, aiming to understand the propagation of the geothermal cold front and reservoir pressures.
During the 3 d International Continental Scientific Drilling Program (ICDP)-sponsored UrbEnLab workshop, 75 scientists from 17 countries met in Delft, the Netherlands, in June 2024 to prioritize the scientific ambitions of the deep exploration and monitoring borehole and to discuss potential techniques that could be applied to tackle them. Assessing the life cycle of a geothermal system situated in a complex heterogeneous sedimentary system was defined as the broad aim, with revealing the detailed flow field established being a key priority.
On the TU Delft campus, we aim to drill a borehole of around 4.5 km depth to be used for the exploration, observation, and monitoring of subsurface processes that will be part of a larger research infrastructure under development. This so-called urban energy laboratory includes – in addition to the deep multi-use borehole – a well-instrumented geothermal doublet drilled in 2023, reaching to a depth of 2.2 km; a local seismic monitoring system (installed in 2022); an ultra-sensitive portable seismic monitoring array; and a high-temperature aquifer heat storage system (HT-ATES), for which a pilot well was drilled in 2024. With this urban energy laboratory, we want to tackle problems and better understand processes related to multiple and/or competing subsurface uses in urban environments. The deep exploration and monitoring borehole is designed specifically to monitor fluid and/or flux movement in 3D with unprecedented precision, aiming to understand the propagation of the geothermal cold front and reservoir pressures.
During the 3 d International Continental Scientific Drilling Program (ICDP)-sponsored UrbEnLab workshop, 75 scientists from 17 countries met in Delft, the Netherlands, in June 2024 to prioritize the scientific ambitions of the deep exploration and monitoring borehole and to discuss potential techniques that could be applied to tackle them. Assessing the life cycle of a geothermal system situated in a complex heterogeneous sedimentary system was defined as the broad aim, with revealing the detailed flow field established being a key priority.
Journal article
(2025)
-
Paula Rulff, Ute Weckmann, Thomas Kalscheuer, Gregor Willkommen, Laura Maria Buntin, Anna Platz
The pathways of fluids and mantle-originated carbon dioxide in the seismically active Ohře (Eger) Rift system appearing as mofettes at the surface are currently subject to investigation, especially by the International Continental Scientific Drilling Program “Drilling the Eger Rift”. If the aquifers show significant contrast in electrical resistivity to the host rocks, they can be investigated with geo-electromagnetic methods. However, imaging complex fluid and CO
2 pathways in detail in near-surface structures is challenging, because, in contrast to the background stratigraphy, they are often oriented in near-vertical directions. Therefore, we aim to investigate how the shallow aquifer structures can be examined best with an inductive electromagnetic method. For this purpose, we collected radio-magnetotelluric data in the Hartoušov mofette field and evaluated them by two- and three-dimensional inversions. Data from a nearby magnetotelluric station, drill hole data, gas flux measurements and electrical resistivity tomography models were used to assess the reliability and robustness of our inversion results. We concluded that the near-surface fluid reservoirs are adequately depictable, while the migration paths of gaseous CO
2 cannot be traced properly due to a lack of resistivity contrast. Our model analyses suggest that imaging the given geological setting with fluids and gases ascending in anastomosing pathways benefits from a fine-scale three-dimensional inversion approach because the fluids mostly appear as local conductive reservoir-like anomalies, which can be falsely projected onto the profiles during inversion in two dimensions. The resistivity models contribute with detailed images of the near-surface aquifers to the geodynamic model of the Ohře Rift.
...
The pathways of fluids and mantle-originated carbon dioxide in the seismically active Ohře (Eger) Rift system appearing as mofettes at the surface are currently subject to investigation, especially by the International Continental Scientific Drilling Program “Drilling the Eger Rift”. If the aquifers show significant contrast in electrical resistivity to the host rocks, they can be investigated with geo-electromagnetic methods. However, imaging complex fluid and CO 2 pathways in detail in near-surface structures is challenging, because, in contrast to the background stratigraphy, they are often oriented in near-vertical directions. Therefore, we aim to investigate how the shallow aquifer structures can be examined best with an inductive electromagnetic method. For this purpose, we collected radio-magnetotelluric data in the Hartoušov mofette field and evaluated them by two- and three-dimensional inversions. Data from a nearby magnetotelluric station, drill hole data, gas flux measurements and electrical resistivity tomography models were used to assess the reliability and robustness of our inversion results. We concluded that the near-surface fluid reservoirs are adequately depictable, while the migration paths of gaseous CO 2 cannot be traced properly due to a lack of resistivity contrast. Our model analyses suggest that imaging the given geological setting with fluids and gases ascending in anastomosing pathways benefits from a fine-scale three-dimensional inversion approach because the fluids mostly appear as local conductive reservoir-like anomalies, which can be falsely projected onto the profiles during inversion in two dimensions. The resistivity models contribute with detailed images of the near-surface aquifers to the geodynamic model of the Ohře Rift.
elfe3D v1.0.1
Modelling with the total electric field approach using finite elements in 3D
The controlled-source electromagnetic method is a geophysical technique that detects variations in electric and magnetic material properties in the subsurface. The method is valuable for both environmental assessments and resource exploration. To evaluate electromagnetic data obtained over complex subsurface structures, precise three-dimensional numerical modeling software is required. This is particularly important for tasks related to survey design and for inversion routines that estimate subsurface models based on field data inputs. The forward modelling program elfe3D simulates synthetic data of frequency-domain controlled-source electromagnetic configurations for three-dimensional subsurface models. To optimise the balance between the size of the problems and the accuracy of the solutions, adaptive mesh refinement strategies are employed. The elfe3D program allows the user to define an arbitrary survey setup: The sensor locations, where the electric and magnetic field responses are calculated, can be placed at or below the subsurface, or in the air. This flexibility facilitates realistic modelling of surface-based, borehole and airborne controlled-source electromagnetic surveys.
...
The controlled-source electromagnetic method is a geophysical technique that detects variations in electric and magnetic material properties in the subsurface. The method is valuable for both environmental assessments and resource exploration. To evaluate electromagnetic data obtained over complex subsurface structures, precise three-dimensional numerical modeling software is required. This is particularly important for tasks related to survey design and for inversion routines that estimate subsurface models based on field data inputs. The forward modelling program elfe3D simulates synthetic data of frequency-domain controlled-source electromagnetic configurations for three-dimensional subsurface models. To optimise the balance between the size of the problems and the accuracy of the solutions, adaptive mesh refinement strategies are employed. The elfe3D program allows the user to define an arbitrary survey setup: The sensor locations, where the electric and magnetic field responses are calculated, can be placed at or below the subsurface, or in the air. This flexibility facilitates realistic modelling of surface-based, borehole and airborne controlled-source electromagnetic surveys.
Geoelectrical and electromagnetic imaging methods applied to groundwater systems
Recent advances and future potentials
Abstract
(2024)
-
Paula Rulff, Octavio Castillo-Reyes, Philipp Koyan, Tina Martin, Wouter Deleersnyder, Maria Carrizo Mascarell
The impacts of climate change, combined with population growth, necessitate practical and effective solutions for locating groundwater resources and ensuring drinking water quality. Our contribution explores recent advances in geoelectrical and electromagnetic imaging methods applied to investigate groundwater systems. Geoelectrical and electromagnetic imaging techniques are popular methods for characterising subsurface properties, such as electrical resistivity or dielectric permittivity. These electrical properties are strongly related to the hydrogeological characteristics of the subsurface. Therefore, geoelectrical and electromagnetic investigations can provide valuable insights into finding groundwater resources, assessing the water quality in terms of contaminations and conducting effective groundwater management.
Our study examines state-of-the-art approaches in modelling and instrumentation of induced polarisation and electrical resistivity tomography, as well as time- and frequency-domain electromagnetics and ground-penetrating radar methods. We review recent impactful and innovative groundwater case studies where the above-mentioned methods were applied and further developed. Emphasising the combination of geoelectrical and electromagnetic methods, the studies provide insights into the variation of electrical subsurface properties at different scales, contributing to an improved understanding of the hydrological dynamics in the studied areas. Furthermore, we provide an outlook on the potential for applying geoelectrical and electromagnetic imaging techniques for large-scale groundwater investigations in the exascale computing area. ...
Our study examines state-of-the-art approaches in modelling and instrumentation of induced polarisation and electrical resistivity tomography, as well as time- and frequency-domain electromagnetics and ground-penetrating radar methods. We review recent impactful and innovative groundwater case studies where the above-mentioned methods were applied and further developed. Emphasising the combination of geoelectrical and electromagnetic methods, the studies provide insights into the variation of electrical subsurface properties at different scales, contributing to an improved understanding of the hydrological dynamics in the studied areas. Furthermore, we provide an outlook on the potential for applying geoelectrical and electromagnetic imaging techniques for large-scale groundwater investigations in the exascale computing area. ...
The impacts of climate change, combined with population growth, necessitate practical and effective solutions for locating groundwater resources and ensuring drinking water quality. Our contribution explores recent advances in geoelectrical and electromagnetic imaging methods applied to investigate groundwater systems. Geoelectrical and electromagnetic imaging techniques are popular methods for characterising subsurface properties, such as electrical resistivity or dielectric permittivity. These electrical properties are strongly related to the hydrogeological characteristics of the subsurface. Therefore, geoelectrical and electromagnetic investigations can provide valuable insights into finding groundwater resources, assessing the water quality in terms of contaminations and conducting effective groundwater management.
Our study examines state-of-the-art approaches in modelling and instrumentation of induced polarisation and electrical resistivity tomography, as well as time- and frequency-domain electromagnetics and ground-penetrating radar methods. We review recent impactful and innovative groundwater case studies where the above-mentioned methods were applied and further developed. Emphasising the combination of geoelectrical and electromagnetic methods, the studies provide insights into the variation of electrical subsurface properties at different scales, contributing to an improved understanding of the hydrological dynamics in the studied areas. Furthermore, we provide an outlook on the potential for applying geoelectrical and electromagnetic imaging techniques for large-scale groundwater investigations in the exascale computing area.
Our study examines state-of-the-art approaches in modelling and instrumentation of induced polarisation and electrical resistivity tomography, as well as time- and frequency-domain electromagnetics and ground-penetrating radar methods. We review recent impactful and innovative groundwater case studies where the above-mentioned methods were applied and further developed. Emphasising the combination of geoelectrical and electromagnetic methods, the studies provide insights into the variation of electrical subsurface properties at different scales, contributing to an improved understanding of the hydrological dynamics in the studied areas. Furthermore, we provide an outlook on the potential for applying geoelectrical and electromagnetic imaging techniques for large-scale groundwater investigations in the exascale computing area.
We develop a three-dimensional inversion code to image the resistivity distribution of the subsurface from frequency-domain controlled-source electromagnetic data. Controlled-source electromagnetic investigations play an important role in many different geophysical prospecting applications. To evaluate controlled-source electromagnetic data collected with complex measurement setups, advanced three-dimensional modelling and inversion tools are required.
We adopt a preconditioned non-linear conjugate gradient algorithm to enable three-dimensional inversion of impedance tensor and vertical magnetic transfer function data produced by multiple sets of two independent active sources. Forward simulations are performed with a finite-element solver. Increased sensitivities at source locations can optionally be counteracted with a weighting function in the regularization term to reduce source-related anomalies in the resistivity model. We investigate the capabilities of the inversion code using one synthetic and one field example. The results demonstrate that we can produce reliable subsurface models, although data sets from single pairs of independent sources remain challenging. ...
We adopt a preconditioned non-linear conjugate gradient algorithm to enable three-dimensional inversion of impedance tensor and vertical magnetic transfer function data produced by multiple sets of two independent active sources. Forward simulations are performed with a finite-element solver. Increased sensitivities at source locations can optionally be counteracted with a weighting function in the regularization term to reduce source-related anomalies in the resistivity model. We investigate the capabilities of the inversion code using one synthetic and one field example. The results demonstrate that we can produce reliable subsurface models, although data sets from single pairs of independent sources remain challenging. ...
We develop a three-dimensional inversion code to image the resistivity distribution of the subsurface from frequency-domain controlled-source electromagnetic data. Controlled-source electromagnetic investigations play an important role in many different geophysical prospecting applications. To evaluate controlled-source electromagnetic data collected with complex measurement setups, advanced three-dimensional modelling and inversion tools are required.
We adopt a preconditioned non-linear conjugate gradient algorithm to enable three-dimensional inversion of impedance tensor and vertical magnetic transfer function data produced by multiple sets of two independent active sources. Forward simulations are performed with a finite-element solver. Increased sensitivities at source locations can optionally be counteracted with a weighting function in the regularization term to reduce source-related anomalies in the resistivity model. We investigate the capabilities of the inversion code using one synthetic and one field example. The results demonstrate that we can produce reliable subsurface models, although data sets from single pairs of independent sources remain challenging.
We adopt a preconditioned non-linear conjugate gradient algorithm to enable three-dimensional inversion of impedance tensor and vertical magnetic transfer function data produced by multiple sets of two independent active sources. Forward simulations are performed with a finite-element solver. Increased sensitivities at source locations can optionally be counteracted with a weighting function in the regularization term to reduce source-related anomalies in the resistivity model. We investigate the capabilities of the inversion code using one synthetic and one field example. The results demonstrate that we can produce reliable subsurface models, although data sets from single pairs of independent sources remain challenging.
Monitoring temperature changes in geothermal applications is crucial to ensure sustainable heat production and storage operations. This work focuses on a geothermal project situated on the campus of Delft University of Technology in the Netherlands. In connection to deep low-enthalpy geothermal reservoir exploration, an aquifer thermal energy storage system for the purpose of seasonal shallow heat storage is planned. To enable monitoring changes in the electrical resistivity distribution due to heat injection and extraction operations in the shallow subsurface, a new 480 m deep borehole will be equipped with an electrode setup. Measuring the vertical component of the electric field in the borehole using a frequency-domain surface-to-borehole controlled-source electromagnetic setup is particular effective for monitoring reservoir changes. Due to corrosion effects, conventional electrodes have rather limited lifespans, which may not be sufficient for the multi-decade operational plan for this geothermal application.
We use a new approach integrating capacitive electrodes in composite borehole casings. Tests in shallow boreholes have shown comparable results to standard electrodes. Integrating the capacitive sensors in the composite borehole casing is rather time- and cost-intense requiring pre-drilling installation and specially designed electronics. Therefore, we want to optimise the electrode placement along the borehole trajectory. We simulate vertical electric fields at closely spaced receivers along the borehole trajectory for different subsurface scenarios originating from resistivity changes introduced by injecting and extracting hot fluids. Applying the Ramer–Douglas–Peucker algorithm, we determine which electrode locations and combinations are optimal to record resulting variations in the vertical electric field component. The proposed methodology for optimal electrode placement is promising to improve monitoring efficiency in geothermal applications, ensuring sustainable and effective operation over extended periods. ...
We use a new approach integrating capacitive electrodes in composite borehole casings. Tests in shallow boreholes have shown comparable results to standard electrodes. Integrating the capacitive sensors in the composite borehole casing is rather time- and cost-intense requiring pre-drilling installation and specially designed electronics. Therefore, we want to optimise the electrode placement along the borehole trajectory. We simulate vertical electric fields at closely spaced receivers along the borehole trajectory for different subsurface scenarios originating from resistivity changes introduced by injecting and extracting hot fluids. Applying the Ramer–Douglas–Peucker algorithm, we determine which electrode locations and combinations are optimal to record resulting variations in the vertical electric field component. The proposed methodology for optimal electrode placement is promising to improve monitoring efficiency in geothermal applications, ensuring sustainable and effective operation over extended periods. ...
Monitoring temperature changes in geothermal applications is crucial to ensure sustainable heat production and storage operations. This work focuses on a geothermal project situated on the campus of Delft University of Technology in the Netherlands. In connection to deep low-enthalpy geothermal reservoir exploration, an aquifer thermal energy storage system for the purpose of seasonal shallow heat storage is planned. To enable monitoring changes in the electrical resistivity distribution due to heat injection and extraction operations in the shallow subsurface, a new 480 m deep borehole will be equipped with an electrode setup. Measuring the vertical component of the electric field in the borehole using a frequency-domain surface-to-borehole controlled-source electromagnetic setup is particular effective for monitoring reservoir changes. Due to corrosion effects, conventional electrodes have rather limited lifespans, which may not be sufficient for the multi-decade operational plan for this geothermal application.
We use a new approach integrating capacitive electrodes in composite borehole casings. Tests in shallow boreholes have shown comparable results to standard electrodes. Integrating the capacitive sensors in the composite borehole casing is rather time- and cost-intense requiring pre-drilling installation and specially designed electronics. Therefore, we want to optimise the electrode placement along the borehole trajectory. We simulate vertical electric fields at closely spaced receivers along the borehole trajectory for different subsurface scenarios originating from resistivity changes introduced by injecting and extracting hot fluids. Applying the Ramer–Douglas–Peucker algorithm, we determine which electrode locations and combinations are optimal to record resulting variations in the vertical electric field component. The proposed methodology for optimal electrode placement is promising to improve monitoring efficiency in geothermal applications, ensuring sustainable and effective operation over extended periods.
We use a new approach integrating capacitive electrodes in composite borehole casings. Tests in shallow boreholes have shown comparable results to standard electrodes. Integrating the capacitive sensors in the composite borehole casing is rather time- and cost-intense requiring pre-drilling installation and specially designed electronics. Therefore, we want to optimise the electrode placement along the borehole trajectory. We simulate vertical electric fields at closely spaced receivers along the borehole trajectory for different subsurface scenarios originating from resistivity changes introduced by injecting and extracting hot fluids. Applying the Ramer–Douglas–Peucker algorithm, we determine which electrode locations and combinations are optimal to record resulting variations in the vertical electric field component. The proposed methodology for optimal electrode placement is promising to improve monitoring efficiency in geothermal applications, ensuring sustainable and effective operation over extended periods.
Journal article
(2024)
-
Paula Rulff, Thomas Kalscheuer
Controlled-source electromagnetic methods are applied to survey the electrical resistivity distribution of the subsurface. This work compares normalized electromagnetic field components and transfer functions such as impedance tensors and vertical magnetic transfer functions generated by two independent source polarizations as input data for three-dimensional inversion. As most other available inversion codes allow for inverting only one of the mentioned input data types, it is unclear which data type is preferable for controlled-source electromagnetic inversion. Our three-dimensional non-linear conjugate gradient inversion code can handle both input data types, facilitating a comparison of normalized electromagnetic field components and transfer functions inversion. Examining inversion results for a three-dimensional synthetic model with two anomalies, we infer that the transfer functions inversion is favourable for recovering the overall resistivity distribution below the receiver sites in fewer iterations. The inversion of normalized electromagnetic field components produces a sharper image of the anomalies and may be capable of detecting the resistivity distribution below the extended sources, which comes at the price of introducing a more heterogeneous background resistivity in the model.
...
Controlled-source electromagnetic methods are applied to survey the electrical resistivity distribution of the subsurface. This work compares normalized electromagnetic field components and transfer functions such as impedance tensors and vertical magnetic transfer functions generated by two independent source polarizations as input data for three-dimensional inversion. As most other available inversion codes allow for inverting only one of the mentioned input data types, it is unclear which data type is preferable for controlled-source electromagnetic inversion. Our three-dimensional non-linear conjugate gradient inversion code can handle both input data types, facilitating a comparison of normalized electromagnetic field components and transfer functions inversion. Examining inversion results for a three-dimensional synthetic model with two anomalies, we infer that the transfer functions inversion is favourable for recovering the overall resistivity distribution below the receiver sites in fewer iterations. The inversion of normalized electromagnetic field components produces a sharper image of the anomalies and may be capable of detecting the resistivity distribution below the extended sources, which comes at the price of introducing a more heterogeneous background resistivity in the model.
Journal article
(2023)
-
Octavio Castillo-Reyes, Paula Rulff, Evan Schankee Um, Adrian Amor-Martin
In 3D geo-electromagnetic modeling, an adequate discretisation of the modeling domain is crucial to obtain accurate forward responses and reliable inversion results while reducing the computational cost. This paper investigates the mesh design for subsurface models, including steel-cased wells, which is relevant for many exploration settings but still remains a numerically challenging task. Applying a goal-oriented mesh refinement technique and subsequent calculations with the high-order edge finite element method, simulations of 3D controlled-source electromagnetic models in the presence of metallic infrastructure are performed. Two test models are considered, each needing a distinct version of approximation methods to incorporate the conductive steel casings of the included wells. The influence of mesh quality, goal-oriented meshing, and high-order approximations on problem sizes, computational cost, and accuracy of electromagnetic responses is investigated. The main insights of our work are: (a) the applied numerical schemes can mitigate the computational burden of geo-electromagnetic modeling in the presence of steel artifacts; (b) investigating the processes driving the meshing of models with embedded metallic infrastructures can lead to adequate strategies to deal with the inversion of such electromagnetic data sets. Based on the modeling results and analyses conducted, general recommendations for modeling strategies are proposed when performing simulations for challenging steel infrastructure scenarios.
...
In 3D geo-electromagnetic modeling, an adequate discretisation of the modeling domain is crucial to obtain accurate forward responses and reliable inversion results while reducing the computational cost. This paper investigates the mesh design for subsurface models, including steel-cased wells, which is relevant for many exploration settings but still remains a numerically challenging task. Applying a goal-oriented mesh refinement technique and subsequent calculations with the high-order edge finite element method, simulations of 3D controlled-source electromagnetic models in the presence of metallic infrastructure are performed. Two test models are considered, each needing a distinct version of approximation methods to incorporate the conductive steel casings of the included wells. The influence of mesh quality, goal-oriented meshing, and high-order approximations on problem sizes, computational cost, and accuracy of electromagnetic responses is investigated. The main insights of our work are: (a) the applied numerical schemes can mitigate the computational burden of geo-electromagnetic modeling in the presence of steel artifacts; (b) investigating the processes driving the meshing of models with embedded metallic infrastructures can lead to adequate strategies to deal with the inversion of such electromagnetic data sets. Based on the modeling results and analyses conducted, general recommendations for modeling strategies are proposed when performing simulations for challenging steel infrastructure scenarios.
Abstract
(2023)
-
Paula Rulff, Thomas Kalscheuer, Dominik Zbinden
We present a synthetic inversion study illustrating two approaches which help to deal with heterogeneous sensitivities in 3D frequency-domain controlled-source electromagnetic inverse problems.
Using edge-based vector finite-element approximations on tetrahedral meshes and a preconditioned non-linear conjugate gradient algorithm, we invert for impedance tensor elements generated by a set of two coincident perpendicularly oriented horizontal electric or horizontal magnetic dipole sources. Depending on the number and locations of sources and the choice of impedance tensor components used for inversion, the sensitivity patterns can differ significantly. Measurement setups with a small number of sources, but many receiver stations at the surface covering near-field, transition zone and far-field, are often deployed for land-based controlled-source electromagnetic measurements. Such a setup can result in accumulated sensitivities close to the sources and receivers, which implies strongest model updates in these regions and can mislead the inverse algorithm to a search direction, where no physically meaningful model can be produced nor the data are fitted.
In order to mitigate the influence of strong sensitivities near sources and receivers on the inversion process, we apply an efficient preconditioner and customisable weights in the model regularisation matrix. The preconditioner is updated with the Broyden-Fletcher-Goldfarb-Shanno algorithm using the diagonal of the approximate Hessian matrix as start preconditioner. The latter is computationally expensive to obtain, but aims at finding a favourable search direction for the inverse algorithm already in early iterations and distributing the model update more evenly in the domain. To account for the sensitivity loss with depth, we implemented a depth weighting functional in the model regularisation term. The approach is based on counteracting the exponential and geometrical decay of the electromagnetic fields with depth and distance from the sources. In practical, we increase the smoothing in the shallow part of the model close to the source locations, where no structure is expected. We present synthetic examples indicating that this approach is an efficient way of helping the inversion to converge, obtaining a reliable model and resolving structure at depth. ...
Using edge-based vector finite-element approximations on tetrahedral meshes and a preconditioned non-linear conjugate gradient algorithm, we invert for impedance tensor elements generated by a set of two coincident perpendicularly oriented horizontal electric or horizontal magnetic dipole sources. Depending on the number and locations of sources and the choice of impedance tensor components used for inversion, the sensitivity patterns can differ significantly. Measurement setups with a small number of sources, but many receiver stations at the surface covering near-field, transition zone and far-field, are often deployed for land-based controlled-source electromagnetic measurements. Such a setup can result in accumulated sensitivities close to the sources and receivers, which implies strongest model updates in these regions and can mislead the inverse algorithm to a search direction, where no physically meaningful model can be produced nor the data are fitted.
In order to mitigate the influence of strong sensitivities near sources and receivers on the inversion process, we apply an efficient preconditioner and customisable weights in the model regularisation matrix. The preconditioner is updated with the Broyden-Fletcher-Goldfarb-Shanno algorithm using the diagonal of the approximate Hessian matrix as start preconditioner. The latter is computationally expensive to obtain, but aims at finding a favourable search direction for the inverse algorithm already in early iterations and distributing the model update more evenly in the domain. To account for the sensitivity loss with depth, we implemented a depth weighting functional in the model regularisation term. The approach is based on counteracting the exponential and geometrical decay of the electromagnetic fields with depth and distance from the sources. In practical, we increase the smoothing in the shallow part of the model close to the source locations, where no structure is expected. We present synthetic examples indicating that this approach is an efficient way of helping the inversion to converge, obtaining a reliable model and resolving structure at depth. ...
We present a synthetic inversion study illustrating two approaches which help to deal with heterogeneous sensitivities in 3D frequency-domain controlled-source electromagnetic inverse problems.
Using edge-based vector finite-element approximations on tetrahedral meshes and a preconditioned non-linear conjugate gradient algorithm, we invert for impedance tensor elements generated by a set of two coincident perpendicularly oriented horizontal electric or horizontal magnetic dipole sources. Depending on the number and locations of sources and the choice of impedance tensor components used for inversion, the sensitivity patterns can differ significantly. Measurement setups with a small number of sources, but many receiver stations at the surface covering near-field, transition zone and far-field, are often deployed for land-based controlled-source electromagnetic measurements. Such a setup can result in accumulated sensitivities close to the sources and receivers, which implies strongest model updates in these regions and can mislead the inverse algorithm to a search direction, where no physically meaningful model can be produced nor the data are fitted.
In order to mitigate the influence of strong sensitivities near sources and receivers on the inversion process, we apply an efficient preconditioner and customisable weights in the model regularisation matrix. The preconditioner is updated with the Broyden-Fletcher-Goldfarb-Shanno algorithm using the diagonal of the approximate Hessian matrix as start preconditioner. The latter is computationally expensive to obtain, but aims at finding a favourable search direction for the inverse algorithm already in early iterations and distributing the model update more evenly in the domain. To account for the sensitivity loss with depth, we implemented a depth weighting functional in the model regularisation term. The approach is based on counteracting the exponential and geometrical decay of the electromagnetic fields with depth and distance from the sources. In practical, we increase the smoothing in the shallow part of the model close to the source locations, where no structure is expected. We present synthetic examples indicating that this approach is an efficient way of helping the inversion to converge, obtaining a reliable model and resolving structure at depth.
Using edge-based vector finite-element approximations on tetrahedral meshes and a preconditioned non-linear conjugate gradient algorithm, we invert for impedance tensor elements generated by a set of two coincident perpendicularly oriented horizontal electric or horizontal magnetic dipole sources. Depending on the number and locations of sources and the choice of impedance tensor components used for inversion, the sensitivity patterns can differ significantly. Measurement setups with a small number of sources, but many receiver stations at the surface covering near-field, transition zone and far-field, are often deployed for land-based controlled-source electromagnetic measurements. Such a setup can result in accumulated sensitivities close to the sources and receivers, which implies strongest model updates in these regions and can mislead the inverse algorithm to a search direction, where no physically meaningful model can be produced nor the data are fitted.
In order to mitigate the influence of strong sensitivities near sources and receivers on the inversion process, we apply an efficient preconditioner and customisable weights in the model regularisation matrix. The preconditioner is updated with the Broyden-Fletcher-Goldfarb-Shanno algorithm using the diagonal of the approximate Hessian matrix as start preconditioner. The latter is computationally expensive to obtain, but aims at finding a favourable search direction for the inverse algorithm already in early iterations and distributing the model update more evenly in the domain. To account for the sensitivity loss with depth, we implemented a depth weighting functional in the model regularisation term. The approach is based on counteracting the exponential and geometrical decay of the electromagnetic fields with depth and distance from the sources. In practical, we increase the smoothing in the shallow part of the model close to the source locations, where no structure is expected. We present synthetic examples indicating that this approach is an efficient way of helping the inversion to converge, obtaining a reliable model and resolving structure at depth.
Abstract
(2021)
-
Michael Weiss, Paula Rulff, Thomas Kalscheuer
We developed two forward modelling approaches to simulate 3-dimensional land-based controlled-source electromagnetic (CSEM) problems in frequency domain with hexahedral spectral-element meshes and tetrahedral finite-element meshes. In recent years, the geo-electromagnetic community made a lot of progress in modelling and inversion of EM data in three dimensions using a variety of approaches. The available software is used to verify the accuracy of newly developed codes, which apply e.g. different element shapes or interpolation schemes. However, a direct comparison in terms of advantages and disadvantages of different modelling strategies, especially discretisation methods in 3D, is often not focused on in publications.
Having two modelling codes and their developers available at the same place, gives us the unique opportunity to compare the approaches in a very detailed way. Our spectral-element as well as our finite-element solution is based on Galerkin’s weighted residual method and we solve the electromagnetic diffusion equations for the total electric field on the element edges.
The main differences between both codes are the choice and order of the interpolation functions and the discretisation of the modelling domain employing hexahedral and tetrahedral elements. While the tetrahedral meshes used in our finite-element approach are known for being able to properly resolve complex structures in the subsurface, this issue is addressed in the spectral-element method by utilising curvilinear instead of orthogonal hexahedral elements.
In this contribution, we focus on the comparison of both approaches for a simple 1D model and a complex 3D model in terms of accuracy, effort in mesh generation and computational resources such as simulation time and memory requirement. Moreover, we contrast the influence of mesh discretisation on the solution for the two methods as well as the order of approximation. A preliminary test simulation of a model consisting of a conductive body buried within a resistive background covered by a thin conductive layer yielded comparable results in terms of accuracy. It also revealed significant differences concerning the mesh discretisation meaning the solution's dependency on the meshing of the model domain. ...
Having two modelling codes and their developers available at the same place, gives us the unique opportunity to compare the approaches in a very detailed way. Our spectral-element as well as our finite-element solution is based on Galerkin’s weighted residual method and we solve the electromagnetic diffusion equations for the total electric field on the element edges.
The main differences between both codes are the choice and order of the interpolation functions and the discretisation of the modelling domain employing hexahedral and tetrahedral elements. While the tetrahedral meshes used in our finite-element approach are known for being able to properly resolve complex structures in the subsurface, this issue is addressed in the spectral-element method by utilising curvilinear instead of orthogonal hexahedral elements.
In this contribution, we focus on the comparison of both approaches for a simple 1D model and a complex 3D model in terms of accuracy, effort in mesh generation and computational resources such as simulation time and memory requirement. Moreover, we contrast the influence of mesh discretisation on the solution for the two methods as well as the order of approximation. A preliminary test simulation of a model consisting of a conductive body buried within a resistive background covered by a thin conductive layer yielded comparable results in terms of accuracy. It also revealed significant differences concerning the mesh discretisation meaning the solution's dependency on the meshing of the model domain. ...
We developed two forward modelling approaches to simulate 3-dimensional land-based controlled-source electromagnetic (CSEM) problems in frequency domain with hexahedral spectral-element meshes and tetrahedral finite-element meshes. In recent years, the geo-electromagnetic community made a lot of progress in modelling and inversion of EM data in three dimensions using a variety of approaches. The available software is used to verify the accuracy of newly developed codes, which apply e.g. different element shapes or interpolation schemes. However, a direct comparison in terms of advantages and disadvantages of different modelling strategies, especially discretisation methods in 3D, is often not focused on in publications.
Having two modelling codes and their developers available at the same place, gives us the unique opportunity to compare the approaches in a very detailed way. Our spectral-element as well as our finite-element solution is based on Galerkin’s weighted residual method and we solve the electromagnetic diffusion equations for the total electric field on the element edges.
The main differences between both codes are the choice and order of the interpolation functions and the discretisation of the modelling domain employing hexahedral and tetrahedral elements. While the tetrahedral meshes used in our finite-element approach are known for being able to properly resolve complex structures in the subsurface, this issue is addressed in the spectral-element method by utilising curvilinear instead of orthogonal hexahedral elements.
In this contribution, we focus on the comparison of both approaches for a simple 1D model and a complex 3D model in terms of accuracy, effort in mesh generation and computational resources such as simulation time and memory requirement. Moreover, we contrast the influence of mesh discretisation on the solution for the two methods as well as the order of approximation. A preliminary test simulation of a model consisting of a conductive body buried within a resistive background covered by a thin conductive layer yielded comparable results in terms of accuracy. It also revealed significant differences concerning the mesh discretisation meaning the solution's dependency on the meshing of the model domain.
Having two modelling codes and their developers available at the same place, gives us the unique opportunity to compare the approaches in a very detailed way. Our spectral-element as well as our finite-element solution is based on Galerkin’s weighted residual method and we solve the electromagnetic diffusion equations for the total electric field on the element edges.
The main differences between both codes are the choice and order of the interpolation functions and the discretisation of the modelling domain employing hexahedral and tetrahedral elements. While the tetrahedral meshes used in our finite-element approach are known for being able to properly resolve complex structures in the subsurface, this issue is addressed in the spectral-element method by utilising curvilinear instead of orthogonal hexahedral elements.
In this contribution, we focus on the comparison of both approaches for a simple 1D model and a complex 3D model in terms of accuracy, effort in mesh generation and computational resources such as simulation time and memory requirement. Moreover, we contrast the influence of mesh discretisation on the solution for the two methods as well as the order of approximation. A preliminary test simulation of a model consisting of a conductive body buried within a resistive background covered by a thin conductive layer yielded comparable results in terms of accuracy. It also revealed significant differences concerning the mesh discretisation meaning the solution's dependency on the meshing of the model domain.
Magnetotellurics in the Eger Rift
An overview of subsurface imaging of different tectonic features
Abstract
(2021)
-
Ute Weckmann, Anna Platz, Basel Aleid, Gregor Willkommen, Johannes Mair, Radek Klanica, Svetlana Kovacikova, Paula Rulff
The Bohemian Massif represents the easternmost part of the geodynamically active European Cenozoic Rift System. This region hosts the contact between three tectonic units of the Variscan Belt, the NE-SW trending Eger Rift and the NNW-SSE striking Marianské Lázne fault. It is characterised by ongoing magmatic processes in the intra-continental lithospheric mantle, repeated earthquake swarms, extensive CO2 degassing in mineral springs and mofettes and the presence of Quaternary volcanoes. While the ICDP drilling programme utilizes information gathered within shallow boreholes in the region, we applied the Magnetotelluric (MT) method to obtain site characterizations in the vicinity of the proposed drill sites. The electrical conductivity has proven to be an important parameter to image the above-mentioned tectonic from the lower crust to the shallow subsurface as well as on a regional and a local scale. Here, we present 2D and 3D inversion models of different MT and Radio-MT (RMT) experiments to study e.g. the Hartouŝov mofette fields, the Quaternary scoria cones, the regional faults and their interplay. Thereby the experiments were designed that we can use lower frequency data from MT to support shallow 3D inversions of e.g. the scoria cones in the regions. The most prominent large-scale conductivity features map channels from the lower crust to the surface possibly forming pathways for fluids into the region of earthquake swarms, mofette fields and know spas. However, the locations of the scoria cones seem to be bound to regional fault zones.
...
The Bohemian Massif represents the easternmost part of the geodynamically active European Cenozoic Rift System. This region hosts the contact between three tectonic units of the Variscan Belt, the NE-SW trending Eger Rift and the NNW-SSE striking Marianské Lázne fault. It is characterised by ongoing magmatic processes in the intra-continental lithospheric mantle, repeated earthquake swarms, extensive CO2 degassing in mineral springs and mofettes and the presence of Quaternary volcanoes. While the ICDP drilling programme utilizes information gathered within shallow boreholes in the region, we applied the Magnetotelluric (MT) method to obtain site characterizations in the vicinity of the proposed drill sites. The electrical conductivity has proven to be an important parameter to image the above-mentioned tectonic from the lower crust to the shallow subsurface as well as on a regional and a local scale. Here, we present 2D and 3D inversion models of different MT and Radio-MT (RMT) experiments to study e.g. the Hartouŝov mofette fields, the Quaternary scoria cones, the regional faults and their interplay. Thereby the experiments were designed that we can use lower frequency data from MT to support shallow 3D inversions of e.g. the scoria cones in the regions. The most prominent large-scale conductivity features map channels from the lower crust to the surface possibly forming pathways for fluids into the region of earthquake swarms, mofette fields and know spas. However, the locations of the scoria cones seem to be bound to regional fault zones.
Journal article
(2021)
-
Paula Rulff, Laura M. Buntin, Thomas Kalscheuer
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.
...
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.
3D CSEM Forward Modelling
Testing Adaptive Mesh Refinement Approaches on an Ore Body Model
Abstract
(2020)
-
Paula Rulff, Thomas Kalscheuer, Laura Maria Schmidt, Mehrdad Bastani
We develop a forward modelling code to simulate 3D land-based controlled-source electromagnetic (CSEM) problems in frequency domain with unstructured tetrahedral meshes. The algorithm accounts for isotropic electrical resistivity and magnetic permeability variations in the subsurface. The latest addition to the software is a goal-oriented adaptive mesh refinement strategy driven by error estimators based on “face-jumps” of current density and magnetic flux density. In this study, we demonstrate that the goal-oriented adaptive refinement approaches are suitable to design a problem-specific mesh, which helps to solve 3D CSEM forward problems efficiently and accurately.
In mineral exploration, ore bodies often exhibit a strong resistivity contrast and sometimes a non-negligible contrast in magnetic permeability to their host rock. Accurate 3D modelling of electromagnetic measurement setups is therefore needed for feasibility studies and incorporation of the forward modelling in inversion approaches. To obtain sufficiently accurate solutions in time- and memory efficient computations, one option is to employ guided mesh refinement strategies.
The so called goal-oriented adaptive mesh refinement method aims at designing a mesh, which is fine where necessary and coarse where discretisation errors do not influence the accuracy of the solution at the points of interest, typically the receiver sites. We apply the total electric field approach and first order Nédélec basis functions as interpolation functions defined on the edges of the finite elements to solve the electromagnetic diffusion equations. Thus, we achieve continuity of the electric and magnetic fields inside the elements and tangential to the edges and faces. However, the continuity of the normal components of current density and magnetic flux density across element interfaces cannot be ensured, resulting in small errors in the solution. We calculate these so called “face-jumps” and use them in combination with the elemental residuals and the dual solution of the problem to obtain error estimators that guide our adaptive refinement approach. The dual problem simulates influence sources at the receiver sites to weight the elemental error estimators with their influence to the solution accuracy at the receivers.
We utilise a model of an iron ore body in central Sweden with a known magnetic permeability contrast and unknown electrical resistivity to study the behaviour of our implemented adaptive mesh refinement approaches. This is combined with a feasibility study to investigate the detectability of the ore body with CSEM.
From literature examples on magnetotelluric forward modelling we know, that the error estimator based on the continuity of the normal current density shows robust performance, when modelling for electrical resistivity. We observe the same behaviour after adapting it to the controlled-source problem. The error estimator using the continuity of the magnetic flux density seems mathematically most promising to improve the mesh, when variations in magnetic permeability are significant. Numerical experiments with the ore body model indicate, that best results can be achieved, when mesh refinement guided by both error estimators is applied. ...
In mineral exploration, ore bodies often exhibit a strong resistivity contrast and sometimes a non-negligible contrast in magnetic permeability to their host rock. Accurate 3D modelling of electromagnetic measurement setups is therefore needed for feasibility studies and incorporation of the forward modelling in inversion approaches. To obtain sufficiently accurate solutions in time- and memory efficient computations, one option is to employ guided mesh refinement strategies.
The so called goal-oriented adaptive mesh refinement method aims at designing a mesh, which is fine where necessary and coarse where discretisation errors do not influence the accuracy of the solution at the points of interest, typically the receiver sites. We apply the total electric field approach and first order Nédélec basis functions as interpolation functions defined on the edges of the finite elements to solve the electromagnetic diffusion equations. Thus, we achieve continuity of the electric and magnetic fields inside the elements and tangential to the edges and faces. However, the continuity of the normal components of current density and magnetic flux density across element interfaces cannot be ensured, resulting in small errors in the solution. We calculate these so called “face-jumps” and use them in combination with the elemental residuals and the dual solution of the problem to obtain error estimators that guide our adaptive refinement approach. The dual problem simulates influence sources at the receiver sites to weight the elemental error estimators with their influence to the solution accuracy at the receivers.
We utilise a model of an iron ore body in central Sweden with a known magnetic permeability contrast and unknown electrical resistivity to study the behaviour of our implemented adaptive mesh refinement approaches. This is combined with a feasibility study to investigate the detectability of the ore body with CSEM.
From literature examples on magnetotelluric forward modelling we know, that the error estimator based on the continuity of the normal current density shows robust performance, when modelling for electrical resistivity. We observe the same behaviour after adapting it to the controlled-source problem. The error estimator using the continuity of the magnetic flux density seems mathematically most promising to improve the mesh, when variations in magnetic permeability are significant. Numerical experiments with the ore body model indicate, that best results can be achieved, when mesh refinement guided by both error estimators is applied. ...
We develop a forward modelling code to simulate 3D land-based controlled-source electromagnetic (CSEM) problems in frequency domain with unstructured tetrahedral meshes. The algorithm accounts for isotropic electrical resistivity and magnetic permeability variations in the subsurface. The latest addition to the software is a goal-oriented adaptive mesh refinement strategy driven by error estimators based on “face-jumps” of current density and magnetic flux density. In this study, we demonstrate that the goal-oriented adaptive refinement approaches are suitable to design a problem-specific mesh, which helps to solve 3D CSEM forward problems efficiently and accurately.
In mineral exploration, ore bodies often exhibit a strong resistivity contrast and sometimes a non-negligible contrast in magnetic permeability to their host rock. Accurate 3D modelling of electromagnetic measurement setups is therefore needed for feasibility studies and incorporation of the forward modelling in inversion approaches. To obtain sufficiently accurate solutions in time- and memory efficient computations, one option is to employ guided mesh refinement strategies.
The so called goal-oriented adaptive mesh refinement method aims at designing a mesh, which is fine where necessary and coarse where discretisation errors do not influence the accuracy of the solution at the points of interest, typically the receiver sites. We apply the total electric field approach and first order Nédélec basis functions as interpolation functions defined on the edges of the finite elements to solve the electromagnetic diffusion equations. Thus, we achieve continuity of the electric and magnetic fields inside the elements and tangential to the edges and faces. However, the continuity of the normal components of current density and magnetic flux density across element interfaces cannot be ensured, resulting in small errors in the solution. We calculate these so called “face-jumps” and use them in combination with the elemental residuals and the dual solution of the problem to obtain error estimators that guide our adaptive refinement approach. The dual problem simulates influence sources at the receiver sites to weight the elemental error estimators with their influence to the solution accuracy at the receivers.
We utilise a model of an iron ore body in central Sweden with a known magnetic permeability contrast and unknown electrical resistivity to study the behaviour of our implemented adaptive mesh refinement approaches. This is combined with a feasibility study to investigate the detectability of the ore body with CSEM.
From literature examples on magnetotelluric forward modelling we know, that the error estimator based on the continuity of the normal current density shows robust performance, when modelling for electrical resistivity. We observe the same behaviour after adapting it to the controlled-source problem. The error estimator using the continuity of the magnetic flux density seems mathematically most promising to improve the mesh, when variations in magnetic permeability are significant. Numerical experiments with the ore body model indicate, that best results can be achieved, when mesh refinement guided by both error estimators is applied.
In mineral exploration, ore bodies often exhibit a strong resistivity contrast and sometimes a non-negligible contrast in magnetic permeability to their host rock. Accurate 3D modelling of electromagnetic measurement setups is therefore needed for feasibility studies and incorporation of the forward modelling in inversion approaches. To obtain sufficiently accurate solutions in time- and memory efficient computations, one option is to employ guided mesh refinement strategies.
The so called goal-oriented adaptive mesh refinement method aims at designing a mesh, which is fine where necessary and coarse where discretisation errors do not influence the accuracy of the solution at the points of interest, typically the receiver sites. We apply the total electric field approach and first order Nédélec basis functions as interpolation functions defined on the edges of the finite elements to solve the electromagnetic diffusion equations. Thus, we achieve continuity of the electric and magnetic fields inside the elements and tangential to the edges and faces. However, the continuity of the normal components of current density and magnetic flux density across element interfaces cannot be ensured, resulting in small errors in the solution. We calculate these so called “face-jumps” and use them in combination with the elemental residuals and the dual solution of the problem to obtain error estimators that guide our adaptive refinement approach. The dual problem simulates influence sources at the receiver sites to weight the elemental error estimators with their influence to the solution accuracy at the receivers.
We utilise a model of an iron ore body in central Sweden with a known magnetic permeability contrast and unknown electrical resistivity to study the behaviour of our implemented adaptive mesh refinement approaches. This is combined with a feasibility study to investigate the detectability of the ore body with CSEM.
From literature examples on magnetotelluric forward modelling we know, that the error estimator based on the continuity of the normal current density shows robust performance, when modelling for electrical resistivity. We observe the same behaviour after adapting it to the controlled-source problem. The error estimator using the continuity of the magnetic flux density seems mathematically most promising to improve the mesh, when variations in magnetic permeability are significant. Numerical experiments with the ore body model indicate, that best results can be achieved, when mesh refinement guided by both error estimators is applied.