Research and Application of High-Density Resistivity Method in Seawater Intrusion Investigation in Laizhou Bay

The disaster of seawater intrusion seriously affects people's lives and restricts economic development,so the detection and treatment of seawater intrusion is a long-term task.On the basis of field investigation and water quality analysis,according to the change characteristics of apparent resistivity of groundwater after Cl-reaches 250 mg/L,the theoretical basis for the application of high-density resistivity method was determined,and the characteristic values of apparent resistivity for seawater intrusion interfaces in different geological characteristic regions in Laizhou Bay area were determined by typical profile tests.Combined with water quality investigation and other means,profiles for the high-density resistivity method were arranged,and the interfaces between saline and fresh water were accurately divided.

Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data

The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.

Application of 2-D and 3-D Geo-electrical Resistivity Tomography and Geotechnical Soil Evaluation for Engineering Site Investigation:A Case Study of Okerenkoko Primary School,Warri-Southwest,Delta State,Nigeria

In the design of building structures,joint efforts must be decided to resolve the depth of competent layers across the intended site to safeguard the durability of civil engineering structures and to avert the disastrous consequences of structural failure and collapse.In this study,an integrated methodology that employed DC resistivity tomography involving 2-D and 3-D techniques and geotechnical-soil analysis was used to evaluate subsoil conditions for engineering site investigation at Okerenkoko primary school,in the Warri-southwest area of Delta State,to adduce the phenomena responsible for the visible cracks/structural failure observed in the buildings.The results obtained brought to light the geological structure beneath the subsurface,which consists of four geoelectric layers identified as topsoil,dry/lithified upper sandy layer,wet sand(water-saturated)and peat/clay/sandy clayey soil(highly water-saturated).The deeply-seated peat/clay materials(ρ≤20Ωm)were delineated in the study area to the depths of 17.1 m and 19.8 m from 2-D and 3-D tomography respectively.3-D images presented as horizontal depth slices revealed the dominance of very low resistivity materials i.e.peat/clay/sandy clay within the fourth,fifth and sixth layers at depths ranging from 8.68-12.5 m,12.5-16.9 m and 16.9-21.9 m respectively.The dominance of mechanically unstable peat/clay/sandy clay layers beneath the subsurface,which are highly mobile in response to volumetric changes,is responsible for the noticeable cracks/failure detected on structures within the study site.These observations were validated by a geotechnical test of soil samples in the study area.Atterberg’s limits of the samples revealed plasticity indices of zero.Thus,the soil samples within the depth analyzed were representatives of sandy soil that does not possess any plasticity.The methods justifiably provided relevant information on the subsurface geology beneath the study site and should be appropriated as major tools for engineering site assessment/geotechnical projects.

Application of high-density resistivity method to evaluate safety around Minyan Tunnel,Helong City

By determining the distribution and extent of geological structures surrounding the Mingyan Tunnel,Xicheng Town,Helong City,Jilin Province,we can evaluate the stability of the rock mass and assess potential hazards during tunnel construction.We use the high-density resistivity method to analyze the subsurface structure of the study area.Conductive anomalies are likely to represent joint and fissure systems within strongly weathered host rocks,and the bedrock surrounding the tunnel is relatively stable and does not contain well-developed faults.High-density resistivity analysis can provide valuable information in the context of tunnel engineering and safety.

Viability of high-density resistivity method for evaluating mountain slope stability in Erdaojiang District,Tonghua City,China

The authors employ the high-density resistivity method to image the subsurface structure of a mountain in Erdaojiang District,Tonghua City,Jilin Province,China,to evaluate the potential risk of slope failure on surrounding residential areas and infrastructure,and identify a shallow fault that extends across the center of the mountain and is perpendicular to the mountain slope and accurately locate the spatial position and depth of another fault on the southern side of the mountain.The results provide an important basis for evaluating mountain slope stability.This study also demonstrates that the high-density resistivity method is effective for detecting mountain faults.

Application of high-density resistivity method in archaeological investigation of Sumicheng ancient city site

The authors employ the high-density resistivity method during an archaeological investigation of Sumicheng site,an ancient city of the Tang Dynasty,to find evidence of human activities and locate a favorable target for archaeological excavation in the southern part of the outer city.There are two obvious high-resistivity structures,the south wall of the inner city and an ancient building near the south gate along the outer city wall,of which the resistivities are indicative of rammed soil foundations.The south wall of the inner city is continuous but is cut off abruptly to the east,which we suggest it is due to either wall damage or destruction.The resistivity signature of the target area is verified by archaeological excavation,proving the feasibility and effectiveness of implementing the high-density resistivity method for archaeological exploration.

Application of high-density resistivity method for assessing construction safety of Shimodong tunnel in Helong City of Jilin Province

Some unfavorable geological conditions can affect the construction of tunnels.In order to evaluate the damage degree of tunnel construction and determine the surrounding rock grade and stability of the tunnel,the authors used high-density resistivity method to detect the surrounding rocks of Shimodong tunnel in Xicheng Town of Helong City.The underground resistivity structures of the entrance,exit and middle parts of the tunnel are obtained.Through analysis,it is found that there are no bedrock faults near the tunnel,although some joints and fissures are developed in some locations,which are characterized by low-resistivity anomalies.The tunnel structures are stable overall,favorable for safe and efficient construction.The study also proves the good application effect of the high-density resistivity method in tunnel safety detection.

Landslide monitoring in southwestern China via time-lapse electrical resistivity tomography

The dynamic monitoring of landslides in engineering geology has focused on the correlation among landslide stability,rainwater infiltration,and subsurface hydrogeology.However,the understanding of this complicated correlation is still poor and inadequate.Thus,in this study,we investigated a typical landslide in southwestern China via time-lapse electrical resistivity tomography（TLERT） in November 2013 and August 2014.We studied landslide mechanisms based on the spatiotemporal characteristics of surface water infiltration and flow within the landslide body.Combined with borehole data,inverted resistivity models accurately defined the interface between Quaternary sediments and bedrock.Preferential flow pathways attributed to fracture zones and fissures were also delineated.In addition,we found that surface water permeates through these pathways into the slipping mass and drains away as fissure water in the fractured bedrock,probably causing the weakly weathered layer to gradually soften and erode,eventually leading to a landslide.Clearly,TLERT dynamic monitoring can provide precursory information of critical sliding and can be used in landslide stability analysis and prediction.

Investigation of fractured rock aquifer in South China using electrical resistivity tomography and self-potential methods

Assessment of fractured rock aquifers in many parts of the world is complicated given their strong heterogeneity. Delineation of the subsurface geological formation in the weathered terrain is essential for groundwater exploration. To achieve this goal, 2D electrical resistivity tomography(ERT) and self-potential(SP) in combination with joint profile method(JPM) and boreholes have been carried out to delineate the subsurface geological units, detect the fracture/fault zones in hard rock, monitor the groundwater flow, and estimate the groundwater reserves contained within the weathered terrain at a complex heterogeneous site of Huangbu, South Guangdong of China. The integration of resistivity images with the borehole lithology along three profiles delineates three subsurface distinct layers namely topsoil cover, weathered and unweathered layers. The incorporation of ERT and SP with JPM reveal five fractures/faults, i.e., F_1, F_2, F_3, F_4 and F_5. 2D ERT models interpret the less resistive anomalies as the fractures/faults zones, and high resistive anomalies as the fresh bedrock. The inversion program based on the smoothness-constraint is used on the resistivity field data to get more realistic three layered model. SP measurements are obtained along the same electrical profiles which provide the negative anomalies clearly indicating the groundwater preferential flow pathways along the fracture/fault zones. Hydraulic parameters namely hydraulic conductivity and transmissivity were determined to estimate the groundwater resources contained within the fractures/faults. The integrated results suggest that the fractures/faults zones are most appropriate places of drilling for groundwater exploration in the investigated area. Geophysical methods coupled with the upfront borehole data provides better understanding about the conceptual model of the subsurface geological formations. The current investigation demonstrates the importance of the integrated geophysical methods as a complementary approach for groundwater assessment in the hard rock weathered areas.

Nonlinear inversion for electrical resistivity tomography based on chaotic DE-BP algorithm

Nonlinear resistivity inversion requires efficient artificial neural network(ANN)model for better inversion results.An evolutionary BP neural network(BPNN)approach based on differential evolution(DE)algorithm was presented,which was able to improve global search ability for resistivity tomography 2-D nonlinear inversion.In the proposed method,Tent equation was applied to obtain automatic parameter settings in DE and the restricted parameter Fcrit was used to enhance the ability of converging to global optimum.An implementation of proposed DE-BPNN was given,the network had one hidden layer with 52 nodes and it was trained on 36 datasets and tested on another 4 synthetic datasets.Two abnormity models were used to verify the feasibility and effectiveness of the proposed method,the results show that the proposed DE-BP algorithm has better performance than BP,conventional DE-BP and other chaotic DE-BP methods in stability and accuracy,and higher imaging quality than least square inversion.

2. 5-D RESISTIVITY TOMOGRAPHY USINGBOUNDARY INTEGRAL EQUATIONS

DC Resistivity Tomography is a non-linear inversion problem. So far there are mainly two kinds of inversion methods, based on the finite-element method and alpha centers method. In this paper, the disadvantages of these two kinds of methods were analysed,and a new method of forward modeling and inversion (Tomography) based on boundary integral equations was proposed. This scheme successfuly overcomes the difficulties of the two formarly methods. It isn’t necessary to use the linearization approximation and calculate the Jacobi matrix. Numerical modeling results given in this paper showed that the computation speed of our method is fast, and there is not any special requirement for initial model, and satisfying results of tomography can be obtained in the case of great contrast of conductivity. So it has wide applications.

The application of resistivity tomography to hydrogeological exploration

Resistivity tomography, an advanced technique reconstructing resistivity image from a series of potential survey of electrode arrays is quite different from the non linear inversion of wave problems. The main problem in the resistivity reconstruction is how to solve the Poisson′s equation of direct current fields in voluminous media and to complete the inversion iteration efficiently. A mathematical idea of cascade algorithm proposed by Shima (1992), as an example, is introduced briefly in this paper. The emergence of water flood in tunnels in Wennan Coal Mine, Shandong Province, gave us chances to carry out hydrogeological exploration twice using resistivity tomography in 1995～1996. Three profiles with a total length of 5832 m and a maximum depth of 120 m in tomograms are completed. The series of resulting tomograms distinctly reveal the distributions of stratigraphic structure, mined out areas, fracture zones, crack belts and piping water loss zones. These results from tomograms are verified by drill records and then successfully adopted in the followed hydrogeological engineering in the coal mine. Finally, the authors discussed the technical method and existing problems in resistivity tomography.

Contribution of Electrical Resistivity Tomography and Boring Technique in the Realization of Ten (10) Large Boreholes in a Crystalline Basement Rocks in the Centre-West of Benin

In order to ensure access to drinking water for Benin populations by 2021, the Emergency Measure program for the reinforcement of the drinking water supply system of Savalou city was initiated in 2018. This program focuses on densification and extension of hydraulic infrastructures. Therefore, it is prominent to use rigorous approach for implementation and execution of drilling activities. The present work has the advantage of combining the use of electrical resistivity tomography and borehole technique to locate ten high flow drilling in Savalou city. The electrical resistivity tomography (ERT) panels were made based on the dipole-dipole arrays with 48 electrodes with 5 m inter-electrode spacing. The drilling was carried out over ten selected points and in two stages: confirmation test using piezometer and borehole diameter enlargement. Moreover, only piezometers with flow rate greater than 10 m3/h were enlarged. The tomography processing has identified 10 fractured zones that are defined by 250 - 1000 ohm.m resistivity values and a width between 15 - 55 m. The confirmation test carried out over ten piezometers exhibits high flow rates ranging from 9 to 35 m3/h with depths of 30 to 68 m. Nine over the ten boreholes with a flow rate equal or greater than 10 m3/h, have improved their flow rates by 50% to 100% after the boring technique. Thus, the cumulative flow rate has reached 252. 7 m3/h for Savalou city and his surrounding areas.

Geophysical Investigation of a Solid Waste Disposal Site Using Integrated Electrical Resistivity Tomography and Multichannel Analyses of Surface Waves Methods

Electrical resistivity tomography survey was deployed at a solid waste landfill in southwest Missouri USA with the intent to map variations in moisture content through the solid waste and underlying subsurface, and to map the top of bedrock. Multichannel analyses of surface waves survey was also deployed to map variations in engineering properties of the solid waste and underlying subsurface, and to constrain the interpretations of top of bedrock. The 2-D resistivity images through the waste suggest rainwater seeps through the cap cover system of the solid waste landfill, and moisture content within the solid waste increases with solid waste burial depth. The resistivity anomalies displayed by the soil and bedrock directly underneath the solid waste suggests a lateral component to moisture infiltrating at the toe of the landfill, which is flowing inward to the base of solid waste structural low. The 1-D shear wave velocity profiles obtained from the multichannel analyses of surface waves survey helped interpret the top of bedrock underneath the solid waste, where top of bedrock is difficult to map using electrical resistivity tomography, as shallow fractured bedrock is moist and displays comparable resistivity values to that of overlying soil. Not surprisingly, the top of bedrock is readily identified on the electrical resistivity tomography profiles in places where subsurface is relatively dry. The deployment of the combined non- invasive, cost and time effective geophysical surveys, along with engineering judgement on available site history data, has reasonably identified potential landfill seepage pathways. The methodology presented could be used in similar site investigation settings.

Comparison of smoothness-constrained and geostatistically based cross-borehole electrical resistivity tomography for characterization of solute tracer plumes

Experiments using electrical resistivity tomography(ERT) have shown promising results in reducing the uncertainty of solute plume characteristics related to estimates based on the analysis of local point measurements only.To explore the similarities and differences between two cross-borehole ERT inversion approaches for characterizing salt tracer plumes,namely the classical smoothness-constrained inversion and a geostatistically based approach,we performed two-dimensional synthetic experiments.Simplifying assumptions about the solute transport model and the electrical forward and inverse model allowed us to study the sensitivity of the ERT inversion approaches towards a variety of basic conditions,including the number of boreholes,measurement schemes,contrast between the plume and background electrical conductivity,use of a priori knowledge,and point conditioning.The results show that geostatistically based and smoothness-constrained inversions of electrical resistance data yield plume characteristics of similar quality,which can be further improved when point measurements are incorporated and advantageous measurement schemes are chosen.As expected,an increased number of boreholes included in the ERT measurement layout can highly improve the quality of inferred plume characteristics,while in this case the benefits of point conditioning and advantageous measurement schemes diminish.Both ERT inversion approaches are similarly sensitive to the noise level of the data and the contrast between the solute plume and background electrical conductivity,and robust with regard to biased input parameters,such as mean concentration,variance,and correlation length of the plume.Although sophisticated inversion schemes have recently become available,in which flow and transport as well as electrical forward models are coupled,these schemes effectively rely on a relatively simple geometrical parameterization of the hydrogeological model.Therefore,we believe that standard uncoupled ERT inverse approaches,like the ones discussed and assessed in this paper,will continue to be important to the imaging and characterization of solute plumes in many real-world applications.

Resistivity tomography study on samples with water-bearing structure

The apparent resistivity of the samples with water-bearing configuration was measured by an electrode-array and 2-D resistivity images of these samples were reconstructed then. The obtained series of tomograms reveal the dis tribution and its variation of true resistivity within the samples caused by the changes of crack and liquid distribu tion. Applying this method to the simulation experiment on the electrical properties of rocks, the fracturing and water filling process, which produces the electrical changes, can be brought to light clearly.

Hydrate formation and distribution within unconsolidated sediment:Insights from laboratory electrical resistivity tomography

Laboratory visual detection on the hydrate accumulation process provides an effective and low-cost method to uncover hydrate accumulation mechanisms in nature.However,the spatial hydrate distribution and its dynamic evolutionary behaviors are still not fully understood due to the lack of methods and experimental systems.Toward this goal,we built a two-dimensional electrical resistivity tomography(ERT)apparatus capable of measuring spatial and temporal characteristics of hydrate-bearing porous media.Beach sand(0.05–0.85 mm)was used to form artificial methane hydrate-bearing sediment.The experiments were conducted at 1°C under excess water conditions and the ERT data were acquired and analyzed.This study demonstrates the utility of the ERT method for hydrate mapping in laboratory-scale.The results indicate that the average electrical conductivity decreases nonlinearly with the formation of the hydrate.At some special time-intervals,the average conductivity fluctuates within a certain scope.The plane conductivity fields evolve heterogeneously and the local preferential hydrate-forming positions alternate throughout the experimental duration.We speculate that the combination of hydrate formation itself and salt-removal effect plays a dominant role in the spatial and temporal hydrate distribution,as well as geophysical parameters changing behaviors during hydrate accumulation.

Detecting the Resistivity Distribution of Carbon Fiber Reinforced Concrete by Electrical Resistance Tomography Method

According to the principle of electrical resistance tomography （ ERT）, the resistivity distribution of the carbon fiber reinforced concrete （CFRC） in the sensing field can be measured by injecting exciting current and measuring the voltage on the sensor electrode arrays installed on the surface of the object. Therefore, measurement of the resistivity distribution of CFRC is divided into first measuring the boundary conditions and then inversely computing the resistivity distribution. To reach this goal, an ERT system was constructed, which is composed of a sensor array unit, a data acquisition unit and an image reconstruction unit. Simulations of static ERT was performed on set-ups with many objects spread in a homogeneous background, and a simulation of dynamic ERT was also done on a rectangular board, the resistivity of which was changed within a small domain of it. Then, the resistivity distribution of a CFRC sample with a circlar hole as the target was detected by the ERT system. Simulation and experimental results show that the reconstructed ERT image reflects the resistivity distribution or the resistivity change of CFRC structure well. Especially, a small change in resistivity can be identified from the reconstructed images in the simulation of dynamic ERT images.

Characterisation of Fractures and Fracture Zones in a Carbonate Aquifer Using Electrical Resistivity Tomography and Pricking Probe Methodes

Position, width and fragmentation level of fracture zones and position, sig-nificance and characteristic distance of fractures were aimed to determine in a carbonate aquifer. These are fundamental parameters, e.g. in hydrogeological modelling of aquifers, due to their role in subsurface water movements. The description of small scale fracture systems is however a challenging task. In the test area (Kádárta, Bakony Mts, Hungary), two methods proved to be applicable to get reasonable information about the fractures: Electrical Resistivity Tomography (ERT) and Pricking-Probe (PriP). PriP is a simple mechanical tool which has been successfully applied in archaeological investigations. ERT results demonstrated its applicability in this small scale fracture study. PriP proved to be a good verification tool both for fracture zone mapping and detecting fractures, but in certain areas, it produced different results than the ERT. The applicability of this method has therefore to be tested yet, although its problems most probably origin from human activity which reorganises the near-surface debris distribution. In the test site, both methods displayed fracture zones including a very characteristic one and a number of individual fractures and determined their characteristic distance and significance. Both methods prove to be able to produce hydrogeologically important parameters even individually, but their simultaneous application is recommended to decrease the possible discrepancies.

Electrical Resistivity Tomography and TDEM Applied to Hydrogeological Study in TaubatéBasin, Brazil

This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m of maximum sediments thickness. The study area is in Taubaté city, São Paulo State, Brazil, where the Taubaté aquifer is an important water source. The Taubaté Group is the main sedimentary package of the basin;it is formed mainly by shales that form aquicludes, and thin layers of sandstones that form the aquifer. There are 40 groundwater exploration wells in Taubaté city that provide important information. The study purpose is to characterize the geoelectrical stratigraphy of the subsurface to locate the contact between the Quaternary and Tertiary sediments and to identify the Taubaté aquifer. The ERT is used for shallow investigations (tens of meters) and the TDEM can reach a great investigation depth (hundreds of meters). Therefore, these geophysical methods are complementary. The ERT data were acquired with the pole-dipole array with 20 m of electrodes spacing and 400 m length, and the TDEM data with the central-loop array with a 200 × 200 m transmitter loop. The results permit to define the contact between the Quaternary and Tertiary sediments around 15 m depth, the Pindamonhangaba Formation between 15 m and 30 m depth and the Taubate Group between 30 m and 300 m depth. The TDEM method defined the Taubaté Group as a single geoelectric layer because the shale and the sandstone layers are all very conductive. The basement is formed by gneiss, which is a very resistive rock. The TDEM method is not able to identify a high conductor/resistor contrast. Overall, the results are consistent with the known geology and the wells information.