Search Results
Comparison of novel semi-airborne electromagnetic data with multi-scale geophysical, petrophysical and geological data from Schleiz, Germany
In the framework of the Deep Electromagnetic Sounding for Mineral EXploration (DESMEX) project, we carried out multiple geophysical surveys from regional to local scales in a former mining area in the state of Thuringia, Germany. We prove the applicability of newly developed semi-airborne electromagnetic (EM) systems for mineral exploration by cross-validating inversion results with those of established airborne and ground-based investigation techniques. In addition, supporting petrophysical and geological information to our geophysical measurements allowed the synthesis of all datasets over multiple scales. An initial regional-scale reconnaissance survey was performed with BGR's standard helicopter-borne geophysical system deployed with frequency-domain electromagnetic (HEM), magnetic and radiometric sensors. In addition to geological considerations, the HEM results served as base-line information for the selection of an optimal location for the intermediate-scale semi-airborne EM experiments. The semi-airborne surveys utilized long grounded transmitters and two independent airborne receiver instruments: induction coil magnetometers and SQUID sensors. Due to the limited investigation depth of the HEM method, local-scale electrical resistivity tomography (ERT) and long-offset transient electromagnetic (LOTEM) measurements were carried out on a reference profile, enabling the validation of inversion results at greater depths. The comparison of all inversion results provided a consistent overall resistivity distribution. It further confirmed that both semi-airborne receiver instruments achieve the bandwidth and sensitivity required for the investigation of the resistivity structure down to 1 km depth and therewith the detection of deeply seated earth resources. A 3D geological model, lithological and geophysical borehole logs as well as petrophysical investigations were integrated to interpret of the geophysical results. Distinct highly-conductive anomalies with resistivities of less than 10 Om were identified as alum shales over all scales. Apart from that, the petrophysical investigations exhibited that correlating geophysical and geological information using only one single parameter, such as the electrical resistivity, is hardly possible. Therefore, we developed a first approach based on clustering methods and self-organizing maps (SOMs) that allowed us to assign geological units at the surface to a given combination of geophysical and petrophysical parameters, obtained on different scales. © 2020 The Authors
Towards an open-source landscape for 3-D CSEM modelling
Large-scale modelling of 3-D controlled-source electromagnetic (CSEM) surveys used to be feasible only for large companies and research consortia. This has changed over the last few years, and today there exists a selection of different open-source codes available to everyone. Using four different codes in the Python ecosystem, we perform simulations for increasingly complex models in a shallow marine setting. We first verify the computed fields with semi-analytical solutions for a simple layered model. Then we validate the responses of a more complex block model by comparing results obtained from each code. Finally, we compare the responses of a real-world model with results from the industry. On the one hand, these validations show that the open-source codes are able to compute comparable CSEM responses for challenging, large-scale models. On the other hand, they show many general and method-dependent problems that need to be faced for obtaining accurate results. Our comparison includes finite-element and finite-volume codes using structured rectilinear and octree meshes as well as unstructured tetrahedral meshes. Accurate responses can be obtained independently of the chosen method and the chosen mesh type. The runtime and memory requirements vary greatly based on the choice of iterative or direct solvers. However, we have found that much more time was spent on designing the mesh and setting up the simulations than running the actual computation. The challenging task is, irrespective of the chosen code, to appropriately discretize the model. We provide three models, each with their corresponding discretization and responses of four codes, which can be used for validation of new and existing codes. The collaboration of four code maintainers trying to achieve the same task brought in the end all four codes a significant step further. This includes improved meshing and interpolation capabilities, resulting in shorter runtimes for the same accuracy. We hope that these results may be useful for the CSEM community at large and that we can build over time a suite of benchmarks that will help to increase the confidence in existing and new 3-D CSEM codes.
custEM: Customizable finite-element simulation of complex controlled-source electromagnetic data
We have developed the open-source toolbox custEM (customizable electromagnetic modeling) for the simulation of complex 3D controlled-source electromagnetic (CSEM) problems. It is based on the open-source finite-element library FEniCS, which supports tetrahedral meshes, multiprocessing, higher order polynomials, and anisotropy. We use multiple finite-element approaches to solve the time-harmonic Maxwell equations, which are based on total or secondary electric field and gauged potential formulations. In addition, we develop a secondary magnetic field formulation, showing superior performance if only magnetic fields are required. Using Nédélec basis functions, we robustly incorporate the current density on the edges of the mesh for the total field formulations. The latter enable modeling of CSEM problems taking topography into account. We evaluate semianalytical 1D layered-earth solutions with the pyhed library, supporting arbitrary configurations of dipole or loop sources for secondary field calculations. All system matrices have been modified to be symmetric and solved in parallel with the direct solver MUMPS. Aside from the finite-element kernel, mesh generation, interpolation, and visualization modules have been implemented to simplify and automate the modeling workflow. We prove the capability of custEM, including validation against analytic-solutions, crossvalidation of all implemented approaches, and results for a model with 3D topography with four examples. The object-oriented implementation allows for customizable modifications and additions or to use only submodules designed for special tasks, such as mesh generation or matrix assembly. Therefore, the toolbox is suitable for crossvalidation with other codes and as the basis for developing 3D inversion routines.