One-dimensional Laterally Constrained Joint Anisotropic Inversion of CSRMT and ERT Data

2021 ◽  
Vol 26 (1) ◽  
pp. 35-48
Author(s):  
Arseny Shlykov ◽  
Alexander Saraev ◽  
Sudha Agrahari ◽  
Bülent Tezkan ◽  
Akarsh Singh
Keyword(s):  
High Frequency ◽  
Electrical Resistivity Tomography ◽  
Joint Inversion ◽  
Inversion Algorithm ◽  
Borehole Data ◽  
Resistivity Tomography ◽  
One Dimensional ◽  
Controlled Source ◽  
High Anisotropy ◽  
And Inversion

In this paper, we discuss several approaches for a joint inversion of controlled source radiomagnetotelluric (CSRMT) and electrical resistivity tomography (ERT) data observed over anisotropic media. We compare results of 2D isotropic joint inversion with results of a newly developed joint 1D anisotropic inversion algorithm. The developed algorithm involves the full controlled source high frequency forward and inversion formulations without the plane wave assumption. We demonstrate that for measurements on an anisotropic subsurface the isotropic joint inversion cannot fit both datasets properly due to a high anisotropy of shallow horizons of quaternary sands and loams. The joint anisotropic inversion helps to solve this problem and highlights the advantages of a joint inversion of CSRMT and ERT data. We also demonstrate application of the laterally constrain algorithm for the anisotropic inversion. Results of the joint 1D anisotropic inversion of CSRMT and ERT data were successfully compared with existing borehole data.

scholarly journals Identifying Changes in Sediment Texture along an Ephemeral Gravel-Bed Stream Using Electrical Resistivity Tomography 2D and 3D

2021 ◽  
Vol 11 (7) ◽  
pp. 3030
Author(s):  
Marcos A. Martínez-Segura ◽  
Carmelo Conesa-García ◽  
Pedro Pérez-Cutillas ◽  
Pedro Martínez-Pagán ◽  
Marco D. Vásconez-Maza
Keyword(s):  
Electrical Resistivity ◽  
Cross Sections ◽  
Electrical Resistivity Tomography ◽  
Environmental Changes ◽  
Subsurface Layer ◽  
Sediment Texture ◽  
Borehole Data ◽  
Resistivity Tomography ◽  
Gravel Bed ◽  
Basin Scale

Differences in deposit geometry and texture with depth along ephemeral gravel-bed streams strongly reflect fluctuations in bedload which are due to environmental changes at the basin scale and to morphological channel adjustments. This study combines electrical resistivity tomography (ERT) with datasets from borehole logs to analyse the internal geometry of channel cross-sections in a gravel-bed ephemeral stream (southeast Spain). The survey was performed through longitudinal and transverse profiles in the upper channel stretch, of 14 to 30 m in length and 3 to 6 m in depth, approximately. ERT values were correlated with data on sediment texture as grain size distribution, effective grain sizes, sorting, and particle shape (Zingg’s classification). The alluvial channel-fills showed the superposition of four layers with uneven thickness and arrangement: (1) the softer rocky substrate (<1000 Ω.m); (2) a thicker intermediate layer (1000 to 2000 Ω.m); and (3) an upper set composed of coarse gravel and supported matrix, ranging above 2000 Ω.m, and a narrow subsurface layer, which is the most resistive (>5000 Ω.m), corresponding to the most recent armoured deposits (gravel and pebbles). The ERT results coupled with borehole data allowed for determining the horizontal and vertical behaviour of the materials in a 3D model, facilitating the layer identification.

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 ◽  
...  
Keyword(s):  
Electrical Resistivity ◽  
Electrical Resistivity Tomography ◽  
Geophysical Research ◽  
Resistivity Tomography ◽  
Subsurface Structures ◽  
Data Inversion ◽  
Scripting Language ◽  
And Inversion ◽  
Geoelectric Data ◽  
Insight Into

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.

scholarly journals Characterization of aquifer system using electrical resistivity tomography (ERT) and induced polarisation (IP) techniques

2021 ◽  
Vol 880 (1) ◽  
pp. 012025
Author(s):  
N Akhtar ◽  
M S Mislan ◽  
M I Syakir ◽  
M T Anees ◽  
M S M Yusuff
Keyword(s):  
Electrical Resistivity ◽  
Electrical Resistivity Tomography ◽  
Imaging System ◽  
Borehole Data ◽  
Resistivity Tomography ◽  
Aquifer System ◽  
Saturation Zone ◽  
Average Percentage ◽  
Aquifer Systems ◽  
Freshwater Resource

Abstract Groundwater plays a major role as an alternative freshwater resource for irrigation and industrial purposes. This study aimed to characterize the subsurface of aquifer systems in TelukIntan district, Perak, Malaysia using Electrical Resistivity Tomography (ERT) and Induced Polarization (IP) methods. The horizontal profiling (TL1 and TL2) was conducted at length of 400 m. The estimated depth is 150 m below ground level (b.g.l.). An ABEM SAS 4000 Terrameter and ABEM LUND ES464 Imaging System were applied to create a resistivity pseudo-section using polar-dipole configuration. The collected geo-electrical data was interpreted using RESIST software with partial curve matching and computer iteration. ERT and IP survey profile results were validated with in-situ borehole data from borehole 2 (B2). Eleven samples of soil profile were collected at depth from 5.6 m to 61.2 m, with average percentage of sand, silt and clay are 93.77 %, 5.78 % and 0.02 %, respectively. The geology of subsurface settings is the key factor in determining the aquifer system characterized by interlayer sand-silt sequence indicating the saturation zone of aquifer underlain by shale at the bottom. Further study on hydraulics perspective is important to understand the overall capacity of the aquifer.

Measurement and inversion schemes for single borehole-to-surface electrical resistivity tomography surveys

2011 ◽  
Vol 8 (4) ◽  
pp. 487-497 ◽  
Author(s):  
Panagiotis Tsourlos ◽  
Richard Ogilvy ◽  
Constantinos Papazachos ◽  
Philip Meldrum
Keyword(s):  
Electrical Resistivity ◽  
Electrical Resistivity Tomography ◽  
Resistivity Tomography ◽  
And Inversion

scholarly journals Three-Dimensional Induced Polarization Parallel Inversion Using Nonlinear Conjugate Gradients Method

2015 ◽  
Vol 2015 ◽  
pp. 1-12 ◽  
Author(s):  
Huan Ma ◽  
Handong Tan ◽  
Yue Guo
Keyword(s):  
Message Passing ◽  
Message Passing Interface ◽  
Joint Inversion ◽  
Synthetic Data ◽  
Induced Polarization ◽  
Conjugate Gradients ◽  
Processing Unit ◽  
Inversion Algorithm ◽  
Borehole Data ◽  
Surface Borehole

Four kinds of array of induced polarization (IP) methods (surface, borehole-surface, surface-borehole, and borehole-borehole) are widely used in resource exploration. However, due to the presence of large amounts of the sources, it will take much time to complete the inversion. In the paper, a new parallel algorithm is described which uses message passing interface (MPI) and graphics processing unit (GPU) to accelerate 3D inversion of these four methods. The forward finite differential equation is solved by ILU0 preconditioner and the conjugate gradient (CG) solver. The inverse problem is solved by nonlinear conjugate gradients (NLCG) iteration which is used to calculate one forward and two “pseudo-forward” modelings and update the direction, space, and model in turn. Because each source is independent in forward and “pseudo-forward” modelings, multiprocess modes are opened by calling MPI library. The iterative matrix solver within CULA is called in each process. Some tables and synthetic data examples illustrate that this parallel inversion algorithm is effective. Furthermore, we demonstrate that the joint inversion of surface and borehole data produces resistivity and chargeability results are superior to those obtained from inversions of individual surface data.

Joint inversion of surface and three‐component borehole magnetic data

Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 540-552 ◽  
Author(s):  
Yaoguo Li ◽  
Douglas W. Oldenburg
Keyword(s):  
Field Data ◽  
Joint Inversion ◽  
Magnetic Data ◽  
Data Sets ◽  
Source Depth ◽  
Inversion Algorithm ◽  
Borehole Data ◽  
Weighting Functions ◽  
Surface Data ◽  
Geometric Decay

The inversion of magnetic data is inherently nonunique with respect to the distance between the source and observation locations. This manifests itself as an ambiguity in the source depth when surface data are inverted and as an ambiguity in the distance between the source and boreholes if borehole data are inverted. Joint inversion of surface and borehole data can help to reduce this nonuniqueness. To achieve this, we develop an algorithm for inverting data sets that have arbitrary observation locations in boreholes and above the surface. The algorithm depends upon weighting functions that counteract the geometric decay of magnetic kernels with distance from the observer. We apply these weighting functions to the inversion of three‐component magnetic data collected in boreholes and then to the joint inversion of surface and borehole data. Both synthetic and field data sets are used to illustrate the new inversion algorithm. When borehole data are inverted directly, three‐component data are far more useful in constructing good susceptibility models than are single‐component data. However, either can be used effectively in a joint inversion with surface data to produce models that are superior to those obtained by inversion of surface data alone.

scholarly journals MONITORING RECYCLED WATER INJECTION INTO A CONFINED AQUIFER IN SINDOS (THESSALONIKI) USING ELECTRICAL RESISTIVITY TOMOGRAPHY (ERT): INSTALLATION AND PRELIMINARY RESULTS

2007 ◽  
Vol 40 (2) ◽  
pp. 580 ◽  
Author(s):  
P. Tsourlos ◽  
G. Vargemezis ◽  
C. Voudouris ◽  
T. Spachos ◽  
A. Stampolidis
Keyword(s):  
Electrical Resistivity ◽  
Electrical Resistivity Tomography ◽  
Low Cost ◽  
Water Injection ◽  
General Setting ◽  
Confined Aquifer ◽  
Inversion Algorithm ◽  
Recycled Water ◽  
Resistivity Tomography ◽  
Spatial Coverage

This work describes the installation and preliminary measurements of an electrical resistivity tomography (ERT) system to monitor the injection of recycled water into a confined aquifer in the area of Sindos. The aim is to provide, through time-lapse electrical resistivity tomography (ERT) measurements and processing, geoelectrical images of rather increased volumetric sampling around and between the holes and to obtain improved understanding of the flow and transport of the injected water. The details about the general setting and the design of the injection utility are presented and the construction and installation of the ERT cables into the boreholes are explained in full. Preliminary measurements involving single-hole ERT measurements were obtained and processed with a 2D inversion algorithm to produce images of the subsurface. Results depict a very good correlation between ERT images and the lithology logs an indication of the reliability of the approach. This images offer increased resolution and spatial coverage compared to traditional approaches. The entire ERT permanent installation is of low cost, easy to implement and can be used to understand and evaluate the effects of the water injection process.

Sign in / Sign up

Export Citation Format

Share Document