Estimation of Soil Electrical Properties in a Multilayer Earth Model with Boundary Element Formulation

1 Department of Civil and Structural Engineering, Faculty of Engineering & Built Environment, Universiti Kebangsaan Malaysia, Bangi 43600, Malaysia 2 Department of Electrical and Electronics Engineering, Faculty of Computer Science and Engineering, PSTU, Patuakhali 8602, Bangladesh 3 Department of Electrical and Electronics & System Engineering, Faculty of Engineering & Built Environment, Universiti Kebangsaan Malaysia, Bangi 43600, Malaysia


Introduction
Soil monitoring using electrical resistivity has been widely applied to many geotechnical and other engineering problems to investigate near surface soil profile 1 .In particular, direct current DC resistivity monitoring has been actively used in geotechnical investigations, since the resistivity of subsurface material is easily affected by conductive or resistive fluid injection 2-4 .

Mathematical Problems in Engineering
Near surface soil characterizations and soil strength determinations are prerequisite in highway and road engineering including construction of highway embankments, earth dams, geotechnical engineering, and other divisions of civil engineering.The soil electrical resistivity and electric potential differences can be used for determining the specific soil characteristics 5-7 corresponding to the depth of soil.The model of soil electric properties with multilayer earth structure will be able to perform specified soil characterizations through laver-by-layer earth analysis in geotechnical investigations.This study will encourage innovations in multilayer soil resistivity profile to scrutinize arrangement of specific characteristics of near surface soil such that the objectives of the geotechnical investigation can be realized.
A two-layer soil model is generally used for nonhomogeneous soil characterizations in geoelectric engineering.The soil resistivity measurements are commonly done by four-probe Wenner method.Current is passed through two points current sources, I, at the surface of the earth.Measurements of voltage, V , between two points at the surface of the earth are taken between these points as shown in Figure 1.The electrical resistance, R, is obtained by dividing V by I according to the Ohm's law.The basic principle of the soil resistivity, ρ measurement system is that when a constant voltage is applied to one of the two probes placed in the soil the current that flows between the probes is inversely proportional to the resistance of the soil 8 .
Though the apparent resistivity estimation using the conventional Wenner fourelectrode measurements is relatively simple and its usefulness has always been recognized, it is not a method explicitly directed at multilayer soil structures.The resulting data still shows lack of accuracy for two-layer soil model.The limitations are found especially in terms of the interrelationships between soil apparent electrical resistivity ρ and several soil physical properties.In response to this situation, we have calculated the true resistivity in a horizontally multilayer earth model using the four-probes method.
This study includes resistivity curve features for a variety of earth parameters and has made possible the investigations of techniques for estimating these earth parameters.
Soil electric potentiality and thickness are also considered to obtain efficient outcomes with multilayer analysis.The ability to perform multilayer soil resistivity profile is useful to get bulk density of compacted soil, to estimate the type of soil and rocks, to mark granular soil and groundwater surface, to detect anomalous materials, and to estimate the depths of bedrock surfaces in geotechnical engineering.
In this work, boundary element method BEM is implemented on soil electric potential for electric field analyses which is the physical quantity of interest to analysts 9, 10 .Numerical analysis of multilayer model also includes the boundary conditions to show multilayer model.Differentiating of the voltage also demonstrates how to obtain the electric field intensity for multilayer soil profile.This paper first derives the theoretical equations for calculating earth resistivity in order to prepare the resistivity curves for layer-by-layer analysis.It also describes the features of resistivity curves in a multilayer structure and some numerical analysis for estimating earth parameters.Finally, the model for estimating earth parameters is presented through potential differences and resistivity ratio in a multilayer earth structure.
The model based on soil resistivity ratio and voltage difference of multilayered earth structure is able to define soil characteristics with more specifications in geotechnical profile.Including potential differences and thickness in the model shows the reliable and efficient outcomes in near surface soil investigations.The aim of the research is demonstration of the model for obtaining more specific soil characterizations through multilayer soil resistivity profile.

Wenner Apparent Resistivity
The basic equation for the electric potentiality and soil resistivity 11 is also included in the study.The equation of the electric filed is given corresponding to the gradient of scalar potential as, E −grad ψ.

2.1
Another basic equation of current density, J is that div J 0.

2.2
The conductivity, σ z , varies according to the depth, z.Hence, the partial differential equation of electrical potentiality is expressed as, According to cylindrical coordinates r, z , we obtain 1 r ∂ψ ∂r where r x 2 y 2 .Here, ∂ψ/∂r ∂σ/∂r in 2.4 vanishes as consideration of conductivity, σ z , varies according to depth z.
Due to the azimuthal symmetry, we can separate 2.4 using separation of variables as ψ r, z R r Z z , which is a product of a function of r and a function of z.Thus, where R is the vector of cylindrical coordinates.And where λ is the separation constant.Thus, a general solution on electric potential can be written as 2.7 In addition, the electric potential and electric intensity being varied according to depth is considered as a function of z.Now, the potentiality equation can be derived from 2.7 with Bessel function theory for layer 1 and 2 as, where, ρ is the soil resistivity.Thus, f 1 λ and g 1 λ , f 2 λ and g 2 λ can be solved with following boundary conditions: Now, f 1 λ and g 1 λ are obtained according to the boundary condition: Mathematical Problems in Engineering where reflection coefficient, K, is defined as In the Wenner array, four electrodes are arranged in a straight line with equal spacing of a 12 .For the two-layer model, the ratio of apparent resistivity, ρ a , and upper layer resistivity are estimated according to the relative spacing a/h including various values of the factor, K, as shown in Figure 2. Curves of previously used two-layer model include upper layer resistivity, ρ 1 , and apparent resistivity, ρ a , of having electrical properties in near surface soil investigations.

Methodology
The In the study, the parameters of soil resistivity, potential differences, and thickness of each layer are considered for the implementation of multilayer earth model as Figure 3. Four-probe resistivity method is used to estimate soil electric properties in geotechnical investigations 13 .When current is injected in the soil, electric field is generated in the near surface soil.The potentiality and electric field also vary according to changing of depth in soil 14 .
The n potential functions can be written as n layers model for normal surface in multilayer earth structure.The derivation of potential function of the uppermost layer can be revealed with Bessel function: B 0 m J 0 mr e mz dm .

3.1
The potentiality of ith layer model is shown with 3.2 : Moreover, the potentiality of nth layer is written as A n m J 0 mr e −mz dm .

3.3
A set of n equations is obtained when that potential is continuous across each boundary.The potential equation for multilayer soil profile also shows the derivation of resistivity with Bessel function as, A n 1 m J 0 mr e −mz n dm .

3.4
In addition, a second set of n equations is formed by equating the components of current density normal to the boundary for each boundary: mA n 1 m J 0 mr e −mz n dm.

3.5
All the integrations are carried out through the same limits.Thus, 3.4 can be written as A n m e −mz n B n m e mz n ρ n − A n 1 m e −mz n ρ n 1 J 0 mr dm 0.

3.6
And, we get from 3.5 as A n 1 m e −mz n J 0 mr mdm 0.

3.7
These equations hold for all points along various planes separating the layers

Mathematical Problems in Engineering
These sets of equations are solved directly for all the A and B parameters by working from the above equations: 3.9 The kernel function of the integration, A 1 m B 1 m for the multilayer soil structure, is obtained which is related to the resistivity depth function ρ z .The multilayer reflection coefficient is shown as With subtraction and summation of the resistivity, the multilayer reflection coefficient can be represented as Here, multilayer reflection coefficient, K 1 , is derived to get the function of resistivity ration for the adjacent two layers in multilayer resistivity model.The algorithm involving iterative search for multilayer layer earth structure is given below 1 Estimate soil resistance, R, resistivity, ρ, voltage, v, and thickness, t with four-probe soil resistivity measurements technique.
2 Obtain n layers according to resistivity and potential data in soil electric field.
3 Apply boundary element formulations with k nodes for each layer and iteration are implemented using array technique.
where ψ i is the electric potential in volts at coordinates x, y , ρ is resistivity and Q is the current source.BEM is also used in the proposed multilayer model to derive electric potential in terms of potential functions ψ i based on cylindrical coordinate system.

Results and Discussions
Recent advancements in electronics and electrical engineering have improved the ability to obtain near surface soil characteristics with electrical properties in geotechnical engineering.Four-probe Wenner method is commonly used with two-layer earth model for soil characterizations through electric field generation in soil.Multilayer earth structure reveals the profile of more specific soil properties with different soil electric parameters and has made possible the investigation of techniques for estimating these parameters.Soil resistivity and soil potentiality are obtained through the creation of soil electric field in near surface soil.The study of the soil electric properties for different layer of soil is needed to obtain reliable outcomes in soil characterization.
Figure 4 shows the simulation results for multilayer soil electric field in near surface soil profile.In the present work, we focus our attention on multilayer soil electric model of the electric potentiality and soil resistivity in grounding system.The numerical approach based on the BEM for layered soil model is presented here.Potential measurements on the near surface soil between the metallic probes are taken when the current is injected into the soil.
In the numerical analysis, the soil is modeled in multilayer according to the soil characterization method.The bottom layer is considered as a half-space domain and bounded with the implementation of boundary condition.The parameters of thickness of soil layers are considered for multilayer model including the interface of soil layer and probes.
Soil resistivity is obtained with the consideration of voltage differences between two probes and injected current using four-probe electric method.The potentiality and electric field strength also decrease according to the increasing of the depth of soil profile shown in 2.8 .The soil profile is modeled in multilayer for numerical analysis according to the collected soil resistivity in soil electric field investigations.in methodology of the study as the resistivity ratio of two adjacent layers in multilayer soil electric model.The resistivity ratio and potential differences for the adjacent two layers are considered in the new proposed model of multilayer soil structure.The variation of reflection coefficient with depth is shown in Figure 6.
In addition, the resistivity ratio is considered for the implementation of proposed multilayer soil electric model.The results of numerical analysis including resistivity ratio are revealed through Figure 7 for multilayer earth structure on soil characterizations.In the half space of the multilayer soil electric model, the thickness is considered as infinity for near surface soil profile.The resistivity ratio of the half space is not considerable to get a simple calculation in multilayer soil investigations.As an example, the resistivity ratio, ρ 2 /ρ 1 is demonstrated as 0.73, where resistivity, ρ 1 considered as 0.38 Kilo Ohm-m and ρ 2 is as 0.28 Kilo Ohm-m for easy analyzing in the multi-layer earth structure.The resistivity of third layerρ 3 is taken as 0.22 Kilo Ohm-m that manifests the resistivity ratio, ρ 3 /ρ 2 as 0.78.
However, there are the effects of soil type, water contents, bulk density of soil on resistivity profile of soil.Consideration of only resistivity profile for soil characterizations can vary for the presence of water contents.Thus, the resistivity ratio for multilayer is considered in proposed model to obtain reliable and specific near surface soil characteristics with electrical properties.
Seedher and Arora 16 shows two layer model for measurements of near surface soil resistivity.This two-layer model using four-probe Wenner method also considers the apparent resistivity which can show unreliable results of near surface soil profile.Moreover, there are no specific characterizations for different layers of near surface soil with two-layer soil resistivity model.
There is the implementation of BEM by De Lacerda et al. 17 for the analysis of cathodic protection systems of buried slender structures.In that paper, the soil resistivity properties are measured along the depth where two-layer apparent resistivity model has been used.In addition, Binley and Kemna 18 shows the implementation of DC resistivity for soil characterizations.They demonstrated the surface imaging with soil apparent resistivity using Wenner method.Though apparent resistivity can be used to obtain deeper soil profile, there are limitations in obtaining accurate results for empirical relationship between probe distances and depth of corresponding soil profile.In addition, considering soil apparent resistivity can affect the results of near surface soil investigations.
Therefore, there is concentration on multilayer true resistivity model to obtain accurate results in soil investigations.The theoretical derivations and numerical analysis of multilayer true resistivity model are demonstrated for the application of near surface soil characterizations.Moreover, as the experience of using this model, it is easy to set up in field as there is no need for excavation of hole in surface soil in geotechnical investigations which is tedious and time consuming.

Conclusion
Multilayer soil electric model is demonstrated for soil electrical characterizations with robust and more reliable performance in near surface soil profile.The theoretical equations with numerical analysis on soil electric field are derived for multilayer soil investigations.The feasibility of implementation of resistivity ratio and potential difference for robust multilayer soil profile is shown with numerical analysis of electrical properties.The boundary element formulations are considered in the numerical analysis on layered soil models.The electrical properties with soil multilayer profile are crucial to obtain specific and reliable outcomes for getting depth corresponding soil characteristics in geoelectric engineering.This multilayer soil electric model would be the improvement in obtaining bulk density of compacted soil and detecting hazardous organic waste in geotechnical engineering.

Figure 1 :
Figure 1: Soil resistivity measurements in near surface profile.

1 Figure 2 :
Figure 2: Resistivity ratio for two-layer model using Wenner method.
mz 1 B 0 m e mz 1 A 1 m e −mz 1 − B 1 m e mz 1 0, ρ n A n e −mz n ρ n B n m e mz n − ρ n 1 A n 1 m e −mz n 0, −A n e −mz n B n m e mz n A n 1 m e −mz n 0.

Figure 5
shows the variation of resistivity with depth in a multilayer soil profile.The theoretical derivations for multilayer soil resistivity profile are shown based on the wave propagation solution of electrical engineering.Reflection coefficient K 1 is derived Mathematical Problems in Engineering f a c e , Y ( m ) S u r fa c e , X ( m )

Figure 4 :
Figure 4: 3D profile of soil electric potential in near surface soil.

Figure 5 :
Figure 5: Multiple-layer considerations according to collected soil resistivity data.

Figure 7 :
Figure 7: Depth corresponding resistivity ratio in soil investigations.
research study on the soil characterizations with multilayer electric model is conducted by geoelectric research group at Geotechnical Laboratory, Faculty of Engineering and Built Environment, University Kebangsaan Malaysia UKM .The analysis for multilayer soil characterizations also includes BEM numerical analysis.Simulations for multilayer layer earth structures are conducted with partial differential equation PDE tools in Matlab 2009 including boundary element formulation.Simulations with PDE tools are done for the conditions of conductive media DC for electric field generations in soil.The data collection for obtaining multilayer electrical properties with four-probe Wenner method is carried out.
BEM is applied to obtain multilayer model with boundary conditions in soil investigations.BEM criteria 15, 16 is applied to the electric potential of the soil as the solution of Poisson equation is shown in 3.14