Uncertainties Analysis of Electrical Resistivity Tomography

Author(s):  
Yonatan Garkebo Doyoro ◽  
Chang Ping-Yu ◽  
Jordi Mahardika Puntu

<p>We examined the uncertainty of the resistivity method in cavity studies using a synthetic cavity model set at six-different depths. Conceptual models were simulated to generate synthetic resistivity data for dipole-dipole, pole-dipole, Wenner-Schlumberger, and pole-pole arrays. The 2D geoelectric models were recovered from the inversion of the synthetically measured resistivity data. The highest anomaly effect (1.46) and variance (24400) in resistivity data were recovered by dipole-dipole array, while the pole-pole array obtained the lowest anomaly effect (0.60) and variance (2401) for the target cavity T<sub>1</sub>. The anomaly effect and variance were linearly associated with the quality of the inverted models. The steeper anomaly gradient of resistivity indicated more distinct cavity boundaries, while the gentler gradient prevents the inference of the cavity boundaries. The recovered model zone above the depth of investigation index of 0.1 has shown relatively higher sensitivity. Modeling for dipole-dipole array provided the highest model resolution and anomaly gradient that shows a relatively distinct geometry of the cavity anomalies. On the contrary, the pole-dipole and Wenner-Schlumberger arrays recovered good model resolutions and moderate anomaly gradient but determining the anomaly geometries is relatively challenging. Whereas, the pole-pole array depicted the lowest model resolution and anomaly gradient with less clear geometry of the cavity anomalies. At deeper depths, the inverted models showed a reduction in model resolutions, overestimation in anomaly sizes, and deviation in anomaly positions, which can create ambiguity in resistivity model interpretations. Despite these uncertainties, our modeling specified that the 2D resistivity imaging is a potential technique to study subsurface cavities.</p>

2021 ◽  
Vol 11 (7) ◽  
pp. 3143
Author(s):  
Yonatan Garkebo Doyoro ◽  
Ping-Yu Chang ◽  
Jordi Mahardika Puntu

We examined the uncertainty of the two-dimensional (2D) resistivity method using conceptual cavity models. The experimental cavity study was conducted to validate numerical model results. Spatial resolution and sensitivity to resistivity perturbations were also assessed using checkerboard tests. Conceptual models were simulated to generate synthetic resistivity data for dipole-dipole (DD), pole-dipole (PD), Wenner–Schlumberger (WS), and pole-pole (PP) arrays. The synthetically measured resistivity data were inverted to obtain the geoelectric models. The highest anomaly effect (1.46) and variance (24,400 Ω·m) in resistivity data were recovered by the DD array, whereas the PP array obtained the lowest anomaly effect (0.60) and variance (2401 Ω·m) for the shallowest target cavity set at 2.2 m depth. The anomaly effect and variance showed direct dependency on the quality of the inverted models. The DD array provided the highest model resolution that shows relatively distinct anomaly geometries. In contrast, the PD and WS arrays recovered good resolutions, but it is challenging to determine the correct anomaly geometries with them. The PP array reproduced the lowest resolution with less precise anomaly geometries. Moreover, all the tested arrays showed high sensitivity to the resistivity contrasts at shallow depth. The DD and WS arrays displayed the higher sensitivity to the resistivity perturbations compared to the PD and PP arrays. The inverted models showed a reduction in sensitivity, model resolution, and accuracy at deeper depths, creating ambiguity in resistivity model interpretations. Despite these uncertainties, our modeling specified that two-dimensional resistivity imaging is a potential technique to study subsurface cavities. We inferred that the DD array is the most appropriate for cavity surveys. The PD and WS arrays are adequate, while the PP array is the least suitable for cavity studies.


Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. E11-E24 ◽  
Author(s):  
Anand Singh ◽  
Shashi Prakash Sharma ◽  
İrfan Akca ◽  
Vikas Chand Baranwal

We evaluate the use of a fuzzy c-means clustering procedure to improve an inverted 2D resistivity model within the iterative error minimization procedure. The algorithm is coded in MATLAB language for the Lp-norm inversion of 2D direct current resistivity data and is referred to as fuzzy constrained inversion (FCI). Two additional input parameters are required to be provided by the interpreter: (1) the number of geologic units in the model (i.e., the number of clusters) and (2) the mean resistivity values of each geologic unit (i.e., cluster center values of the geologic units). The efficacy of our approach is evaluated by tests carried on the synthetic and field electrical resistivity tomography (ERT) data. Inversion results from the FCI algorithm are presented for conventional L1- and L2-norm minimization techniques. FCI indicates improvement over conventional inversion approaches in differentiating the geologic units if a proper number of the geologic units is provided to the algorithm. Inappropriate clustering information will affect the resulting resistivity models, particularly conductive geologic units existing in the model. We also determine that FCI is only effective when the observed ERT data can recognize the particular geologic units.


Author(s):  
Olisah Nzemeka

An investigation using electrical resistivity method was conducted around a solid waste dumpsite at Nsukka in Nsukka L.G.A of Enugu State, Nigeria to investigate the level of groundwater contamination. Electrical Sounding (VES) and 2D resistivity imaging were used with a digital read out resistivity meter (ABEM SAS 1000) to acquire data in the area. A total of eight (8) sounding and six (6) 2D resistivity imagings were carried out in the area. A contaminant leachate plume was delineated in 2D resistivity sections as low resistivity zones while the VES shows the depth of aquifer. In 2D pseudosections where bluish colours with low resistivities (less than 20.80Ωm) with the depth ranging from 1.28m to 17.1m in the Line 1 and 2 are seen as contaminated zones. The rest of the lines are not contaminated because of their high resistivities (greater than 20.80Ωm). The result of the electrical resistivity survey also showed 4 - 5 layers geo-electric sections and an AA and AK type sounding curves. The VES result shows that VES 1A, 1B, 2A and 2B which are carried out on line 1 & 2 of the wenner lines showed signs of contamination with low resistivity values less than 20.80Ωm complementing the wenner results. The contamination has not yet got to where the aquifer is located on the lines. Since the depth to the aquifer ranges from 30.26m to 155.43m while maximum depth of contamination is 17.1m. It is believed that the leachate has not percolated down to the aquiferous zones as such aquifers are presumed to be free.


2014 ◽  
Vol 695 ◽  
pp. 823-827
Author(s):  
Siti Norsalkini Mohd Akip Tan ◽  
Edy Tonnizam Mohamad ◽  
Rosli Saad ◽  
Mustaza M. Nordiana

An assessment of subsurface profile of a study area at Nusajaya. Johor was conducted using 2D resistivity method. The aims of the study are to determine the subsurface profile hence contouring the appearance of outcrop by 2D resistivity imaging. Subsurface structure can be delineated based on the physical and chemical parameter which is apparent resistivity of the medium. A survey line of 80m and minimum 2m electrode spacing was adopted using Pole-dipole array. Data acquisition on the survey line was completed and processed by using Res2Dinv software. The inversion model resistivity shows sandstone contains iron mineral (30-250 ohm-m) and weathered sandstone (500-1000 ohm-m). The lower part of the layer represents sandstone and siltstone exhibit high resistivity values of 1500-5000 ohm-m. The range from 30-1000 ohm-m correspond to the grade IV and V of sandstone in the actual profile. Meanwhile, the high range of the resistivity value matched the material with weathering grade II and III in ground mass. The deterioration of most physical properties was associated with the increase number of weathering grade.


2011 ◽  
Vol 3 (4) ◽  
Author(s):  
Eleni Kokinou ◽  
Apostolos Sarris

AbstractThe present survey aimed to image the subsurface structure, including karstic voids, and to evaluate the extent of the heterogeneities that can result in potentially dangerous collapse of road segments overlying these features. A multidisciplinary geophysical approach (seismic refraction, frequency domain electromagnetic and ground penetrating radar) in combination with a detailed geological survey indicated the presence of tectonic faults as well as velocity and conductivity anomalies along an old road within the area of Akrotiri at Chania (Crete). Due to the presence of subsurface fuel pipes, perpendicular to the direction of the road, 2D resistivity imaging was excluded from the applied geophysical methods.Interpretation of the geophysical data revealed that the section of the road investigated overlies prominent voids attributed mostly to karst features. The conductivity and velocity anomalies are interpreted to indicate an area where the host limestone rock has been downthrown by faulting and associated karstification. The continuation of this fault zone was observed in the slope of the road during later excavations. Interpretation, using geographic information systems (GIS) to integrate data, allowed these controls and relationships to be understood and monitored. The above methodology was proved successful for areas where the application of resistivity method is not possible.


Author(s):  
Muhammad Sehah ◽  
Sukmaji Anom Raharjo ◽  
Fajar Destiani

Interpretation of 2D-subsurface rock resistivity data has been carried out in the iron ore prospect area of Eastern Binangun Coastal in ​​Cilacap Regency, Central Java. The background of this research is the potential for abundant iron sand in this area that prospects to be exploited. The research was conducted using a magnetic method in 2017 to map the distribution patterns of the local magnetic anomalies that were interpreted to originate from the distribution of iron ore in the subsurface. In 2018, the research continued using the 2D-resistivity method to find out the lithology section in the subsurface of research area. 2D-resistivity data acquisition is carried out on four tracks consisting of Bng-01 to Bng-04. The resistivity data modeling have produced the true resistivity value for each track in the form of the subsurface resistvity section, which including the Bng-01 track is 2.27 – 44.1 Ωm; the Bng-02 track is 4.5 – 58.6 Ωm; the Bng-03 track is 6.37 – 63.4 Ωm; and the Bng-04 track of 4.98 – 83.3 Ωm. After interpretation process, some models of subsurface rocks lithology section is obtained under the four trajectories. The rocks resulted from interpretation process consists of sand which inserted with gravel (> 58.6 Ωm); sand containing iron ore grains (28.2 – 83.3 Ωm), clayey sand (11.1 – 32.9 Ωm), sandy clay (4.98 – 13.5 Ωm), and fine sand which intruded by saltwater (<6.49 Ωm). Sand containing iron ore grains is main part of the coastal aquifer so that exploitation of iron sand has the potential to reduce aquifer function in storing and flowing of groundwater and causing of abrasion in the Eastern Binangun Coastal area.


2020 ◽  
pp. 1345-1352
Author(s):  
Amina M. Salman ◽  
Jassim M. Thabit ◽  
Ali M. Abed

The 2D resistivity imaging technique was applied in an engineering study for the investigation of subsurface weakness zones within University of Anbar, western Iraq. The survey was carried out using Dipole-dipole array with an n-factor of 6 and a-spacing values of 2 m and 5 m. The inverse models of the 2D electrical imaging clearly show the resistivity contrast between the anomalous parts of the weakness zones and the background resistivity distribution. The thickness and shape of the subsurface weakness zones were well defined from the 2D imaging using Dipole-dipole array of 2 m a-spacing. The thickness of the weakness zone ranges between 9.5 m to 11.5 m. Whereas the Dipole-dipole array with a-spacing of 5 m and n-factor of 6 allocated the geoelectrical stratigraphic layers sequence in low-accuracy of weakness zones, but deeper than the inverse model of 2 m a-spacing. This survey was made to explain the correlation between the weakness zone and the deeper layers in the study area. It points out that the deeper layers were not affected in the weakness zones. The inverse model was produced using the Standard Least-Squares Inversion Method and the Robust Inversion Model Constraints Method. The first method had a gradational boundary of the weakness zones and the second had sharper and straighter boundaries of fractures and voids within the weakness zones.


Sign in / Sign up

Export Citation Format

Share Document