Investigation of Soil Ionization Propagation in Two- Layer Soil Samples A. Elzowawi Cardiff University, UK

Similar documents
Heat Transfer in Evacuated Tubular Solar Collectors

Effect of Foundation Embedment of PWR Structures in SSI Analysis considering Spatial Incoherence of Ground Motions

RESPONSE OF ANCHOR IN TWO-PHASE MATERIAL UNDER UPLIFT

TRANSMISSIVITY BEHAVIOR OF SHREDDED SCRAP TIRE DRAINAGE LAYER IN LANDFILL COVER SYSTEM *

Modeling of Ceiling Fan Based on Velocity Measurement for CFD Simulation of Airflow in Large Room

Analysis of Pullout Resistance of Soil-Nailing in Lateritic Soil

STACK EFFECT IN LIGHT WELL OF HIGH RISE APARTMENT BUILDING

FLIGHT UNLOADING IN ROTARY SUGAR DRYERS. P.F. BRITTON, P.A. SCHNEIDER and M.E. SHEEHAN. James Cook University

Heat Transfer Enhancement using Herringbone wavy & Smooth Wavy fin Heat Exchanger for Hydraulic Oil Cooling

Experimental Study of Direct Contact Condensation of Steam on Water Droplets

INFLUENCE OF SOLAR RADIATION AND VENTILATION CONDITIONS ON HEAT BALANCE AND THERMAL COMFORT CONDITIONS IN LIVING-ROOMS

Numerical study of heat pipe application in heat recovery systems

EFFECT OF BOLT CONNECTION OF SQUARE-SHAPED GEOCELL MODEL ON PULLOUT TEST RESULTS

Further Development of an Ion current Combustion Control System (ICCS)

Open and Closed Door Moisture Transport and Corresponding Energy Consumption in Household Refrigerator

Numerical Analysis of the Bearing Capacity of Strip Footing Adjacent to Slope

Simulation Of Pneumatic Drying: Influence Of Particle Diameter And Solid Loading Ratio

Numerical Simulation of Fluid Flow and Heat Transfer in a Water Heater

Stability analysis of slopes with surcharge by LEM and FEM

SIMULATION ANALYSIS ON THE FRESH AIR HANDLING UNIT WITH LIQUID DESICCANT TOTAL HEAT RECOVERY

Modeling and Simulation of Axial Fan Using CFD Hemant Kumawat

Consolidation Stress Effect On Strength Of Lime Stabilized Soil

A Method of Heat Exchange Structure Optimization of the Cricoid Plastic Parts

EFFECT OF CENTRAL PILE IN INCREASING THE BEARING CAPACITY OF BORED PILE GROUPS

Effect of Placement of Footing on Stability of Slope

Piles subject to excavation-induced soil movement in clay

OPTIMIZATION OF VENTILATION MODE OF SMOKE CONTROL SYSTEM IN HIGH-RISE BUILDING FIRE

The Research of Performance Comparison of Displacement and Mixing Ventilation System in Catering Kitchen *

Smoldering Propagation Characteristics of Flexible Polyurethane Foam under Different Air Flow Rates

ASSESSMENT OF FIRE BEHAVIOUR OF TIMBER PARTITION MATERIALS WITH A ROOM CALORIMETER

Investigation of moisture content fluctuation in mixedflow

Numerical Stability Analysis of a Natural Circulation Steam Generator with a Non-uniform Heating Profile over the tube length

NUMERICAL SIMULATION OF THE NEW SOUTH WALES FIRE BRIGADE COMPARTMENT FIRE BEHAVIOUR TRAINING TEST CELL

Study and Design Considerations of HRSG Evaporators in Fast Start Combined Cycle Plants. Govind Rengarajan, P.E.CEM, Dan Taylor

A Numerical study of the Fire-extinguishing Performance of Water Mist in an Opening Machinery Space

The Effect of Quantity of Salt on the Drying Characteristics of Fresh Noodles

American International Journal of Research in Science, Technology, Engineering & Mathematics

Dynamic Simulation of Double Pipe Heat Exchanger using MATLAB simulink

International Forum on Energy, Environment Science and Materials (IFEESM 2015)

APPLICATION OF CFD SIMULATIONS FOR LOCATING OUTDOOR UNITS ON HIGH RISE TOWER Basant Kumar Gupta 1, Anand Patil 2

Experimental Study on Match for Indoor and Outdoor Heat Exchanger of Residential Airconditioner

Experimental Research On Gas Injection High Temperature Heat Pump With An Economizer

NUMERICAL STUDIES ON BARE CABIN FIRES WITH OPERATION OF SMOKE EXTRACTION SYSTEM

FIRE DYNAMICS IN FAÇADE FIRE TESTS: Measurement, modeling and repeatability

Maximum Obtainable Efficiency For Engines And Refrigerators Based On The Stirling Cycle

Dynamic Simulation of Liquid Chillers

Behaviour of Black Cotton Soil Reinforced with Sisal Fibre

Analysis of freeze protection methods for recuperators used in energy recovery from exhaust air

EXPERIMENTAL PERFORMANCE INVESTIGATION OF SWIRLING FLOW ENHANCEMENT ON FLUIDIZED BED DRYER

CFD Analysis of Fire Characteristics on Subway Junction Station

Purdue e-pubs. Purdue University

Equipment Required. WaveRunner Zi series Oscilloscope 10:1 High impedance passive probe. 3. Turn off channel 2.

Characterization of a Heat Sink with Embedded Heat Pipe with Variable Heat Dissipating Source Placement for Power Electronics Applications

Impact of Air Distribution on Heat Transfer during Night-Time Ventilation

EFFECTS OF VARIABLES ON NATURAL CONVECTIVE HEAT TRANSFER THROUGH V-CORRUGATED VERTICAL PLATES

Simulation study of evacuation in high-rise buildings

The Low-Temperature Radiant Floor Heating System Design and Experimental Study on its Thermal Comfort

A Study on Cooling Rate with Blade and Sound Fire Extinguisher

NUMERICAL ANALYSIS OF VERTICAL UPLIFT RESISTANCE OF HORIZONTAL STRIP ANCHOR EMBEDDED IN COHESIVE FRICTIONAL WEIGHTLESS SOIL

Performance Enhancement of Inclined Bubbling Fluidized Bed Paddy Dryer by Design Modification

A Study on Soil Stabilization using Cement and Coir Fibres

FROSTING BEHAVIOR OF FIN-TUBE HEAT EXCHANGER ACCORDING TO REFRIGERANT FLOW TYPE

PULLOUT CAPACITY OF HORIZONTAL AND INCLINED PLATE ANCHORS IN CLAYEY SOILS

MODELLING AND SIMULATION OF A ROOM WITH A RADIANT COOLING CEILING. Technicka 4, Prague 6, Czech Republic

Experimental Research on Monomial Cooling Measure of Greenhouse in Summer

Fire Dynamics Simulation and Evacuation for a Large Shopping Center (Mall), Part II, Evacuation Scenarios

MINIATURE HEAT TRANSPORT SYSTEMS WITH LOOP HEAT PIPES

EXPERIMENTAL STUDY ON PULL-OUT CAPACITY OF HELICAL PILE IN CLAYEY SOIL

CHOOSING A FIRE VENTILATION STRATEGY FOR AN UNDERGROUND METRO STATION

Thermal comfort investigation on a naturally ventilated two- storey residential house in Malaysia

Long-range Ultrasonic Non-destructive Testing of Fuel Tanks

A TRNSYS Simulation Case Study on Utilization of Heat Pump For both Heating and Cooling

Passive cooling of telecommunication outdoor cabinets for mobile base station

Development Of 2-Cylinder Rotary Compressor Series For Light Commercial Use With R410A

Evaluation on Combustion Characteristics of Finishing Materials for Exterior Walls

Performance Characteristics of Air-Conditioner Under Tropical Ambient Condition

Study of Soil Cement with Admixture Stabilization for Road Sub-Grade

CFD Model of a Specific Fire Scenario

Performance Comparison of Ejector Expansion Refrigeration Cycle with Throttled Expansion Cycle Using R-170 as Refrigerant

Effects of Non-Uniform Refrigerant and Air Flow Distributions on Finned- Tube Evaporator Performance

Developing a Fire Test Strategy for Storage Protection Under Sloped Ceilings

APPLICATION OF MINI HEAT PIPES FOR THERMAL MANAGEMENT OF OPTO-ELECTRONIC INSTRUMENTS

EFFECT OF COMPACTION ON THE UNSATURATED SHEAR STRENGTH OF A COMPACTED TILL

Load-Carrying Capacity of Stone Column Encased with Geotextile. Anil Kumar Sahu 1 and Ishan Shankar 2

Exercise 2-4. Heat Exchangers (Optional Exercise) EXERCISE OBJECTIVE DISCUSSION OUTLINE. Description of a brazed plate heat exchanger DISCUSSION

COOLING ENERGY SAVING SYSTEM USING A WATER MIST DEVICE ON WINDOWS

ISO AIR FILTERS FOR GENERAL VENTILATION: DETERMINING FRACTIONAL EFFICIENCY

A Method for Predicting the Matric Suction of Unsaturated Soils with Soil Color Recognition

IEC Type Test Report Report No. EU1120-HR Type PDV 100 Polymer Arrester 10,000 A Line Discharge Class 1

Performance Analysis of Electronic Expansion Valve in 1 TR Window Air Conditioner using Various Refrigerants

The Effects of Woody Vegetation on Levees

PERFORMANCE OF GEOSYNTHETICS IN THE FILTRATION OF HIGH WATER CONTENT WASTE MATERIAL

for family composition, time allocation of people, ownership and specification of home appliances. This paper first presents the simulation procedure

Effect of Magnetic Field on Convective Drying Of Ceramic Clay

Performance of Geosynthetics in the Filtration of High Water Content Waste Material

EFFECT OF RELICT JOINTS IN RAIN INDUCED SLOPE FAILURES IN RESIDUAL SOIL

Indoor Climate Control Effect of AAC Panel Heat Capacity Experimental rooms and simulations with three structural materials

Jülich Solar Power Tower System Behavior During Downtime

Improvement of Granular Subgrade Soil by Using Geotextile and Jute Fiber

Effect of characteristics of unsaturated soils on the stability of slopes subject to rainfall

Transcription:

nvestigation of Soil onization Propagation in Two- Layer Soil Samples A. Elzowawi ElzowawiA@cf.ac.uk A. Haddad Haddad@cf.ac.uk H. Griffiths GriffithsH@cf.ac.uk D. Clark ClarkD@cf.ac.uk Abstract- High current lightning strikes into earthing systems can result in ionization in the soil surrounding the earthing electrode. Most of the published studies investigating this phenomenon have assumed uniform one-layer soil, but soil ionization propagation in a multilayered soil sample has not been extensively addressed. Practical soils may consist of several layers with different water contents, and hence soil resistivity will vary continuously with depth. This investigation considers several sand samples, consisting of two layers with different water contents subjected to standard lightning impulse voltages. A rod-plane electrode configuration was constructed inside a cylindrical plastic test rig, in order to house both wet and dry soil test samples. n order to quantify the propagation of ionization inside the test sample, voltage probes were installed along the tube at specific positions. Localized changes in the ionization zone potential could, therefore, be monitored in real time. ndex Terms Lightning impulse, soil breakdown, soil ionisation, soil ionisation propagation, two-layer soil. NTRODUCTON Multilayer soils with variable resistivity values are a common feature of real earthing systems, due to variations in soil composition with depth, localized inhomogeneity and numerous hydrological and geological processes. The resulting distribution of soil resistivity can thus be highly non-uniform. Therefore, in the literature, several research investigations have been reported in which two layer soil scenarios were studied for steady state performance and resistivity measurements [1-5]. However, the high impulse current performance and the soil ionization phenomenon have not been extensively investigated with two-layer soils. While soil ionization in earthing systems has been intensively studied [-8], the majority of these investigations consider a uniform soil in the vicinity of the earthing electrode, which may not always be fully representative of practical soils. Variations in weather conditions lead to changes in water content within the soil strata. This means that the upper and lower layers may have either high or low resistivity, thus affecting soil ionization initiation and propagation in these layers. This, in turn, affects the localized soil resistivity and the earth potential rise (EPR) due to current injection at the earth electrode []. The simplification proposed in [1], considering a bulk resistivity for all layers rather than each layer individually, could be acceptable where high water contents (more than 1%) in all soil layers keep the resistivity low. n the case of poor soil conditions (high resistivity soil), this simplification may not be applicable. The resistivity and thickness of the layers are significant factors that influence the impulse behavior of earthing systems in two-layer soils. Taking advantage of the current s tendency to flow in the lower resistivity soil layer will help understand the behaviour of the current dissipation in each layer. For better electrical safety, it is preferable to have the lower resistivity layer just below the upper surface layer, with an electrode long enough to reach this layer to allow the majority of the current to be dissipated in this layer rather than in the upper surface layer. Soil ionization is thus encouraged to initiate and propagate to a greater depth, reducing the EPR at the ground surface.[] n this paper, a variety of sand samples representing common soil configurations were considered to investigate the initiation and propagation of soil ionization in two-layer soils with different water contents. Two voltage probes installed inside the test sand sample at selected positions were used to measure dynamic changes in the ionization zone potentials. These measurements allow determination of the position along the sample column to which the soil ionization has propagated. From this measurement, an estimate of the velocity of the ionization propagation between the two probes may be determined.. TEST ARRANGEMENT A. Test Setup: A test rig with consisting of a rod-plane electrode configuration was used in these series of tests connected with the test circuit shown in Fig. 1. The two-layer sample was placed in a vertical plastic tube between the two electrodes. A four-stage 4 kv Haefely impulse voltage generator was used to generate the lightning impulse 1.2/5 wave shape. A capacitive voltage divider with ratio 2731:1 was used to measure the applied voltage, and two other dividers with ratios 2:1 and 1:1 were used to measure the voltages inside the sample during the discharge (see tem 7 in Fig. 1 depicting a voltage probe installed in the tube and connected to the voltage divider. A current transformer with a sensitivity.1 V/A was used to measure the current flowing through the sample. The voltages and current signals were captured and recorded on a LeCroy digital oscilloscope. B. Sample Preparation Medium grain size (.25-.) mm sand was used in all the samples. Tap water was also used to make the wetted sand layer, where the water content (wc) was calculated as a percentage of the mass of the dry sand as stated in [11]. The sand and water were thoroughly mixed so that a good moisture distribution among the sand grains is achieved. 78-1-473-82-/15/$31. 215 EEE

Fig. 1. Test circuit This wet sand mixture is placed in the tube and then moderate pressure was applied to the sample to ensure adequate compaction and uniformity of distribution within the tube; this is thought to prevent formation of gaps inside the sample. A similar procedure was followed when setting each layer. at the instant when the ionisation reaches the dry-wet boundary. This indicates a sharp drop in resistivity caused by the soil ionisation in the dry sand layer. Also visible from Fig. 3 theree are two different delay times (t 1, t 2 ), which represent the propagation times for the ionisation expansion in each probe interval. The delay time (t 1 ) represents not only the propagation, but also the initiation time of the ionization at the electrode surface. The delay time for current initiation (t 1 +t 2 ) was found to depend mainly on the applied voltage and the thickness of the dry layer. The above delay times indicate that the propagation of the ionization region in the dry layer has various successive stages. The first stage is the initiation of the ionisation phenomenon around the electrode; the second stage being the initiation of a streamer discharge from the electrode, and the third stage is the propagation of the discharge through the sand with what appears to be to a variable velocity which may be attributable to the driving local electric field. When the dry layer breaks down, nearly the full applied voltage will appear across the wet layer. This explains why the applied voltage did not show any sign of breakdown. After about 3.88µs from the initiation of the current flow, a full breakdown of the sample occurred, causing the collapse of the three voltages due to the breakdown of the wet layer as well.. PROPAGATON OF ONZATON ZONE N DRY AND WET LAYER SAMPLESS Soil ionization may have different initiation and propagation processes in dry and wet conditions; hence, it has different effects on the soil in both cases. n this test, the electrode e is embedded in the upper layer only. Therefore, by changing the positions of the voltage probes in the upper dry and lower wet layers, the propagation of the ionization zone can be tracked in both layers. A. Propagation in Dry Soil n this experiment, the sample consists of a 1 cm column of dry sand above 15 cm of wet sand with 5 % water content (wc). The first voltage divider is connected to the active electrode and the second is connected to the probe installed in the middle of the dry column. The last divider is connected to the probe installed at the interface between the dry and wet layers, as depicted in Fig. 2. The voltage probes are so located to measure the development of the ionisation propagation in the upper dry layer. An ncreasing lightning impulse voltage was applied in steps to the point where a current was detected by the CT. From the waveforms in Fig. 3, it can be seen that the current did not flow until the ionisation has crossed the whole dry layer up to the second voltage probe between the two layers. Hence, the ionisation discharge current flowed from the highest field region at the HV electrode, crossing the dry layer to reach the wet layer. The voltage (V 1 ) can be seen to rise as the ionisation extends to the probe in the middle of the dry layer. Then, the current and the voltage (V 2 ) start to increase at the same instant when the ionisation discharge arrives at the second probe, as seen in Fig. 3. The total voltage developed across the dry layer (V d = V t V 2 ) collapses Fig. 2. Sample and voltage probe arrangement for -A Fig. 3. Voltages and current traces in dry sand (1cm) above wet sand (15cm, 5% wc)

B. Propagation in Wet Soil Fig. 4 shows the configuration of the sample and the voltage probes, where the upper dry layer has a height of 1 cm as in the previous test, but the lower 2cm wetted layer contains higher wc (1%). A voltage probe is installed at the interface between the dry and wet layers, with a second installed at the middle of the wetted layer. Thus, these two voltage probes will measure the potentials in the expected ionisation zone in the wetted layer. When the applied voltage is sufficient to break down the upper dry layer, soil ionisation is initiated in the lower wetted layer and the current exhibited the second peak. The voltage (V 1 ) can be seen to rise very quickly at the instant when the streamer reaches the dry-wet boundary. This voltage then decreases, indicating the initiation of ionisation at the top of the wet layer. The voltage (V 2 ) increased at the same time as (V 1 ) but it continued to increase as the current increased (indicating a linear resistive behaviour) up to the instant when the ionisation propagation reached the second probe, then V 2 started to decrease, as can be seen in Fig. 5. The propagation time (t p ) between the two peaks of (V 1 and V 2 ) represents the propagation velocity of the ionisation zone boundary between the two probes in the lower layer. This test shows that, despite the distance between the active electrode and the wet layer (1 cm dry sand layer), soil ionization is initiated in the lower layer, which may be driven by the local electric field magnitude provided by the discharge to further ionise the soil, allowing the ionisation to continue its propagation from the upper dry layer to the wetted lower layer. n a similar test, wheree 1% wc was used instead of 1% wc, soil ionisation did not initiate in the bottom layer. Therefore, it can be inferred that soil ionisation could be initiated in areas well below the electrode given the right conditions of field and water content. trace as can be seen in Fig. 7. t is obvious from the voltage trace (V 1 ) that the soil ionization only extended from probe (V 1 ) to the lower layer but it did not propagat te to the second probe, where V 1 started with similar behaviour to the current until the ionization reached probe (V 1 ), then started to decline. As can be seen on the figure, (V 2 ) exhibits the same shape as the current trace, indicating a linear behaviour (no ionization). At kv, the ionization zone has propagated further down to the second probe (V 2 ) as shown in Fig. 8. n a similar test, the height of the upper layer was increased to 15cm and the lower layer height was reduced to 15cm, as shown in Fig b. n this case, the threshold voltage to initiate the ionization increased to 2kV instead of 54kV. Thus, the thickness of the upper layer has a major effect on the initiation of the soil ionization phenomenon in two layer soil. Fig. 4. Sample and voltage probes arrangements for -B V. PROPAGATON OF ONZATON ZONE N TWO WETTED SAND LAYERS WTH DFFERENT WATER CONTENTS This investigation examine the scenario of two wetted layers with various water content in each layer, the water contents considered in this test are 1% and 1%. These percentages have been chosen to obtain a clear difference in the performance of the two layers under the applied voltage. n separate one layer soil sample tests, it was found that 1% wc was not enough to help initi iate soil ionization in this test configuration, whilst 1% was very appropriate for the ionization. Therefore, these series of tests will investigate the initiation and propagation of soil ionization in two layer samples, where one layer contains (1% wc) and the other (1% wc). The thickness of these two layers will also be changed as in the following sections. A. Lower Layer of Higher Water Content The upper layer contains 1% wc with 1 cm height, and the bottom layer contains 1% wc with 2 cm height. The first voltage probe (V 1 ) was installed at the interface between the two layers, and the second voltage probe (V 2 ) was installed in the middle of the lower layer as in Fig.. When the applied voltage was around 42 kv, a linear behaviour was dominating the current conduction. However at 54kV, a second peak started to appear in the current Fig. 5. Voltages and current traces in dry sand (1cm) over wet sand (2cm, 1% wc) Fig.. Sample and voltage probe arrangements for V-A

18 15 3 5-3 - - V V2 2 3 4 Fig. 7. Voltages and current traces in 1% wet sand (11cm) above 1% wet sand (2cm), propagation only to probe 7 5 4 3 2 1-1 -1-2 -3 V22 28 24 2 1 8 4-4 -8-5 Fig. shows the currents waveforms obtaiined under different voltage levels when both layers had a 15cm thiickness. The increase of the second peak is thought to represent the growth of the mple as the applied ionization zone further down in the sand sam voltage increased. t is also noticeable thatt the soil ionization initiation time decreases with voltage increase. 7-4kV -52kV -2kV -4kV -8kV 5 4 3 2 1 2 3 4 5 7 8 Fig.. Current traces for ionization propagation in 1% wet sand (15cm) above 1% wet sand (15cm) B. Upper Layer of Higher Water Content The configuration of the sample and the probbes are shown in Fig. 1. The upper layer has a 1cm column with 1% wc, and the lower layer has 1% wc with 2 cm height. The probee () is placed at the Fig. 1. Sample and voltage prrobes arrangements for V-B 8 7 5 4 3 2 1-1 -1 V V2 24 21 18 15 3 7-3 3 Fig. 8. Voltages and current traces in 1% wet sand (11cm) above 1% wet sand (2cm), propagation up to V2 probe p -1-1 interface between the two layers, and (V2) is at the middle of the bottom layer. Fig. 11 shows an exxample of this discharge, where a voltage impulse of 7kV peak waas applied, but no second current peak observed. However, the longeer rise time of the current than the applied voltage may indicate that there was a weak ionization has only initiated around the electrode,, but it did not propagate far from the electrode and did not cross to thhe bottom layer. n a similar test, where the upperr and lower layers have the same height (15 cm) as shown in Fig. 1 b. A small second current peak started to emerge at kv, as can be b seen in Fig.. However, soil ionization has initiated and propagaated only in the upper layer, and it did not cross to the bottom layer. The T second current peaks obtained in this test were lower than those obbtained in the previous tests, even at voltage near the breakdown levell. (µs) 3 Time 5 Fig. 11. Voltages and current traces in 1% wet sand (1cm) above 1% wet sand (22cm) 1 14 1 8 4 2-2 -1-4 -4kV -52kV -kv -4kV -72kV Currents (A) 5 4 3 2 1-1 -1-2 -3 2 3 4 5 7 Tim me (µs) p in 1% wet sand (15cm) Fig.. Current traces for ionization propagation above 1% wet sand (15cm)

This test may indicate that the higher resistivity layer being at the bottom of a two-layered soil is not desirable for the earthing systems under transient conditions. This scenario gives less resistivity reduction, limited ionisation propagation and lower dissipating currents. V. CONCLUSON nvestigation of soil ionization initiation and propagation in several two layer sand samples with various moisture contents was conducted. Soil ionization is thought to cause the breakdown of the dry sand after few propagation stages. The propagated ionization discharge in dry soil could initiate the soil ionization in the wet sand placed underneath the dry layer if the wet layer has sufficient moisture content. This means that the electric field at the tip of the discharge after the breakdown of the dry layer can be high enough to initiate the ionization in the wet layer. The presence of a high resistivity layer on top of a lower resistivity layer offers better resistivity reduction, deeper ionization propagation and higher dissipated current values than the other scenario with the high resistivity layer is at the bottom of the sample. Furthermore, soil ionization does not tend to propagate from the low resistivity to the high resistivity layers, which may be due to the tendency of the current to flow in the less resistivity soil. Given the findings in this investigation, soil ionization should be considered when designing earthing systems in two-layer soil, where the resistivity and the thickness of each layer should be taken in account. REFERENCES [1] C. J. Blattner, "Analysis of Soil Resistivity Test Methods in Two- Layer Earth," EEE Transactions on Power Apparatus and Systems, vol. PAS-14, pp. 33-38, 85. [2] J. L. del Alamo, "A comparison among eight different techniques to achieve an optimum estimation of electrical grounding parameters in two-layered earth," EEE Transactions on Power Delivery, vol. 8, pp. 18-18, 3. [3] J. Ma, F. P. Dawalibi, and R. D. Southey, "On the equivalence of uniform and two-layer soils to multilayer soils in the analysis of grounding systems," EE Proceedings-Generation, Transmission and Distribution, vol. 143, pp. 4-55,. [4] M. Nayel, L. Boyang, T. Yu, and Z. Yingzhen, "Study of Soil Resistivity Measurements in Vertical Two-Layer Soil Model," in Power and Energy Engineering Conference (APPEEC), 2 Asia-Pacific, Shanghai, 2, pp. 1-5. [5] "EEE Guide for Safety in AC Substation Grounding," ANS/EEE Std 8-8, 8. [] A. C. Liew and M. Darveniza, "Dynamic model of impulse characteristics of concentrated earths," Proceedings of the nstitution of Electrical Engineers, vol. 1, pp. 3-135, 74. [7] A. M. Mousa, "The soil ionization gradient associated with discharge of high currents into concentrated electrodes," EEE Transactions on Power Delivery, vol., pp. 1-177, 4. [8] A. Habjanic and M. Trlep, "The simulation of the soil ionization phenomenon around the grounding system by the finite element method," EEE Transactions on Magnetics,, vol. 42, pp. 87-87, 2. [] Y. Liu, N. Theethayi, R. Thottappillil, R. M. Gonzalez, and M. Zitnik, "An improved model for soil ionization around grounding system and its application to stratified soil," Journal of Electrostatics, vol., pp. 23-2, 24. [1] K. J. Nixon,. R. Jandrell, and A. J. Phillips, "A simplified model of the lightning performance of a driven rod earth electrode in multi-layer soil that includes the effect of soil ionisation," in ndustry Applications Conference, 2. 41st AS Annual Meeting. Conference Record of the 2 EEE, Tampa, FL, 2, pp. 1821-1825. [11] "Methods of test for soils for civil engineering purposes Part 1: General requirements and sample preparation," British standard BS 1377-1,.