scholarly journals A review of open software resources in python for electrical resistivity modelling

2022 ◽  
Vol 9 (1) ◽  
Author(s):  
Yonatan Garkebo Doyoro ◽  
Ping-Yu Chang ◽  
Jordi Mahardika Puntu ◽  
Ding-Jiun Lin ◽  
Tran Van Huu ◽  
...  

AbstractGeophysical modelling performs to obtain subsurface structures in agreement with measured data. Freeware algorithms for geoelectrical data inversion have not been widely used in geophysical communities; however, different open-source modelling/inversion algorithms were developed in recent years. In this study, we review the structures and applications of openly Python-based inversion packages, such as pyGIMLi (Python Library for Inversion and Modelling in Geophysics), BERT (Boundless Electrical Resistivity Tomography), ResIPy (Resistivity and Induced Polarization with Python), pyres (Python wrapper for electrical resistivity modelling), and SimPEG (Simulation and Parameter Estimation in Geophysics). In addition, we examine the recovering ability of pyGIMLi, BERT, ResIPy, and SimPEG freeware through inversion of the same synthetic model forward responses. A versatile pyGIMLi freeware is highly suitable for various geophysical data inversion. The SimPEG framework is developed to allow the user to explore, experiment with, and iterate over multiple approaches to the inverse problem. In contrast, BERT, pyres, and ResIPy are exclusively designed for geoelectric data inversion. BERT and pyGIMLi codes can be easily modified for the intended applications. Both pyres and ResIPy use the same mesh designs and inversion algorithms, but pyres uses scripting language, while ResIPy uses a graphical user interface (GUI) that removes the need for text inputs. Our numerical modelling shows that all the tested inversion freeware could be effective for relatively larger targets. pyGIMLi and BERT could also obtain reasonable model resolutions and anomaly accuracies for small-sized subsurface structures. Based on the heterogeneous layered model and experimental target scenario results, the geoelectrical data inversion could be more effective in pyGIMLi, BERT, and SimPEG freeware packages. Moreover, this study can provide insight into implementing suitable inversion freeware for reproducible geophysical research, mainly for geoelectrical modelling.

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. EN85-EN96 ◽  
Author(s):  
Timothy C. Johnson ◽  
Roelof J. Versteeg ◽  
Mark Rockhold ◽  
Lee D. Slater ◽  
Dimitrios Ntarlagiannis ◽  
...  

Continuing advancements in subsurface electrical resistivity tomography (ERT) are increasing its capabilities for understanding shallow subsurface properties and processes. The inability of ERT imaging data to resolve unique subsurface structures and the corresponding need to include constraining information remains one of the greatest limitations, yet provides one of the greatest opportunities for further advancing the utility of the method. We propose a new method of incorporating constraining information into an ERT imaging algorithm in the form of discontinuous boundaries, known values, and spatial covariance information. We demonstrated the approach by imaging a uranium-contaminated wellfield at the Hanford Site in southeastern Washington State, USA. We incorporate into the algorithm known boundary information and spatial covariance structures derived from the highly resolved near-borehole regions of a regularized ERT inversion. The resulting inversion provides a solution which fits the ERT data (given the estimated noise level), honors the spatial covariance structure throughout the model, and is consistent with known bulk-conductivity discontinuities. The results are validated with core-scale measurements, indicating a significant improvement in accuracy over the standard regularized inversion and revealing important subsurface structure known to influence flow and transport at the site.


2011 ◽  
Vol 8 (4) ◽  
pp. 487-497 ◽  
Author(s):  
Panagiotis Tsourlos ◽  
Richard Ogilvy ◽  
Constantinos Papazachos ◽  
Philip Meldrum

2014 ◽  
Vol 661 ◽  
pp. 73-80 ◽  
Author(s):  
Andy Anderson Bery ◽  
Rosli Saad

This paper presents 4-D high resolution electrical resistivity tomography study for slope monitoring using two optimized (modified) arrays of the Wenner-Schlumberger and Pole-Dipole. These optimized resistivity arrays give a total of 2052 number of datum points from each data acquisition set. These two optimized were used because they have more datum points compared to their original arrays. Perhaps these optimized arrays were able to resolve the subsurface structures from surface areas. Inversion results from the computer (mathematical models) suggested that these optimized arrays were able to give in imaging the Earth’s subsurface structures and characterization at different period of time. Furthermore, in this paper, we are presented merging data levels in order to give high resolution in electrical resistivity tomography. Even though the time taken for data acquisition using this technique is twice, the outcome is compromised and reliable which helpful in the Earth’s subsurface interpretation.


Author(s):  
Peter Milo ◽  
◽  
Tomáš Tencer ◽  
Michal Vágner ◽  
Michaela Prišťáková ◽  
...  

The aim of this paper is to present the results of geophysical surveys at the Staré zámky site near Brno-Líšeň, which were carried out in 2019. Electrical resistivity tomography, georadar survey and large-scale magnetic prospection were all carried out there. The primary task of the first two methods was to investigate the inner structure of the still-existing ramparts and to identify their individual structural elements. The magnetic survey focused on the identification of areas where potential archaeological features can be found – together with the identification of previously-unknown fortifications. The surveys were successful: we have found numerous settlement features from different phases of settlement, an early medieval cemetery and fortifications of various types, sizes and state of preservation. The results of previous archaeological research of the site played an important role for interpretation of the geophysical data. Together these results provide important insights into the study of the complex fortified settlement of Staré zámky. Despite our results, some questions which cannot be answered by geophysical research alone remain unanswered.


2019 ◽  
Vol 65 (1) ◽  
pp. 92-104 ◽  
Author(s):  
V. V. Olenchenko ◽  
L. V. Tsibizov ◽  
A. A. Kartoziya ◽  
E. I. Esin

Evolution of permafrost under thermokarst lakes is an actual question in the light of such problems of cryolythic zone research as greenhouse gas emission, permafrost degradation and cryovolcanism. Recently drained thermokarst lake provide an opportunity to study under-lake permafrost state with ground geophysical methods. This lake located on Kurungnakh island (composed of Yedoma ice complex deposits) in the Lena delta was studied with electrical resistivity tomography. Local low-resistivity anomaly in the central part of the lake was found during previous geophysical research. Main goal of this work is detection of residual thermal effect from a frozen under-lake talik in an electric field. Satellite images of different years show that the drainage has taken place about 30 years ago. The area of the lake was covered in 2016 by high-resolution aerial imagery (5 cm/pixel), then digital elevation model was built using photogrammetry. These data were used for geomorphological description of the alas (depression in permafrost after lake drainage). The alas depth reaches 8 m, its size is about 300 × 500 m. It was formed probably on the last stage of Holocene thermokarst activity and it is relatively shallow in comparison to other typical alases on the island. A number of baijarachs (as a result of polygonal ice wedges thawing) were observed on the alas bottom. A line of the steepest slope marks a coast line of the lake, which allows to estimate a volume of water, which was contained in it earlier. Electrical resistivity tomography was implemented on 8 parallel profiles of 235 m. Measurements were conducted with dipole-dipole array. Basing on 3-dimensional inversion results a 3-dimensional resistivity model of under-lake deposits up to 40 m deep was made. Relatively low resistivity area (16–25 kOhm·m in comparison to 50–100 kOhm·m) was registered at the depth of 15–35 m. It is probably linked to a temperature anomaly (–3...–5 °С in comparison with –8.5 °С average value of the region). Therefore the under lake talik was fully frozen but the rest of temperature anomaly is still observable. 3-dimensional finite-element modeling of talik propagation (500 years) and refreezing (30 years) was done in axisymmetric setting taking phase transition into account. Temperature anomaly up to 0 °С in its center at the depth of 35 m was obtained as a result of the modeling. It qualitatively confirms the interpretation of electrical resistivity tomography data on the residual temperature anomaly below the basin of the drained lake.


2018 ◽  
Vol 10 (4) ◽  
pp. 19
Author(s):  
Peter K. Nyabeze ◽  
Oswald Gwavava ◽  
Matome Sekiba

Electrical resistivity tomography (ERT) was used to delineate subsurface structures at Siloam hot spring in the Soutpansberg Basin of South Africa. The spring water was reported to be hottest in South Africa. Groundwater with temperature ranging between 63.3 °C and 67.5°C was observed flowing from the spring. Relatively wet areas were noted to the south and north of the spring. Electrical resistivity arrays with Dipole-dipole and Wenner-Schlumberger configurations were used. Current injection and potential dipole separations of 5 m and 10 m were applied across the spring. Inversion results for arrays with the shorter dipole separation of 5 m had better resolution. Resistivity sections were generated from inversion of ERT data. The resistivity values for the low and high resistive zones were ρa < 10 and ρa > 70 Ωm respectively. Sub-horizontal resistivity zones were obtained from the inversion of Wenner-Schlumberger array survey data. High resistivity zones associated with dykes or sills were delineated on Dipole-dipole derived depth sections. Near vertical structures associated with water bearing faults or fractures were delineated at the Siloam hot spring.


2021 ◽  
Vol 33 (1) ◽  
pp. 73-84
Author(s):  
Olenchenko V. ◽  
◽  
Shnaider S. ◽  
Osipova P. ◽  
Kalganov A. ◽  
...  

Abstract: Ust’-Biyke-I may become one of the most important archaeological complexes in the northern Altai region. Previous research established that the site was first occupied in the developed Mesolithic. Traces of subsequent occupations up to the Bronze Age were later discovered, in addition to burial structures. A loss of archaeological deposits at the mouth of the Biyke River occurred due to natural and anthropogenic activity. In order to identify the most promising areas for future excavations, geophysical surveys of the site were conducted. An area of 31.5 x 24 m was examined with electrical resistivity tomography (ERT) and ground-penetrating radar (GPR). Results revealed an object with a diameter of about 4 m that was not visible on the ground surface. In its center at a depth of 1.5 m, a local anomaly was observed. Another isometric anomaly was observed in the southwestern part at a depth of 0.75 m. Other anomalies were also noted in the surveyed area, some of which are visible in the modern terrain as small depressions. Future excavations will clarify the geophysical survey. Keywords: Altai, archaeological site, geophysical research, electrical resistivity tomography, groundpenetrating radar, settlement Acknowledgements: Geophysical study at the site was conducted with the financial support of the Russian Foundation for Basic Research (Project № 17–29–04122), archaeological work was carried out at the expense of the project NGS-67706R-20.


Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. A25-A28 ◽  
Author(s):  
Kamini Singha ◽  
Stephen Moysey

A number of issues impact electrical resistivity tomography (ERT) inversions: how ERT measurements sample the subsurface, the nature of subsurface heterogeneity, the geometry selected for data collection, the choice of data-misfit criteria, and regularization of the inverse problem. Lab-scale rock-physics models, typically used to estimate solute concentration from ERT, do not accommodate or account for these issues and therefore produce inaccurate geophysical estimates of solute concentrations. In contrast, the influence of measurement sensitivity and inversion artifacts can be captured by pixel-based rock-physics relationships, determined using numerical analogs that recreate the field-scale ERT experiment based on flow and transport modeling and a priori data. In the 2D synthetic example presented here, where ERT is used to monitor the transport of a saline tracer through the subsurface, improved estimates of concentration are obtained when field-scale rock-physics relationships based on numerical analogs are used.


Sign in / Sign up

Export Citation Format

Share Document