Ground Resistance Estimation for a DC Converter Station: Theory versus Experiment

: This paper is aimed at proposing a calculation model for the ground resistance of a grounding scheme servicing a high-voltage direct-current converter station. The method is based on the equivalence of current conduction and electric field from the grounding scheme through the surrounding medium. The grounding scheme is composed of three concentric ring electrodes supported by two horizontal conductors and eight vertical rods. The calculated ground resistance is 4.8 Ω against the experimental value of 5 Ω with an error of 4.2%. The calculated ground resistance value agrees reasonably well with that of 4.7 Ω as obtained using CYMGRD software (version 7.0). The calculated surface-potential values over the ground surface agreed reasonably well with those measured experimentally, with an average deviation not exceeding 6.5%. This study is designed to investigate how ground resistance is decreased by the increase in the scheme parameters, including the rods’ diameter and length, as well as the radius of the inner and outer rings. The dependency of the ground resistance on the soil type is also investigated.


Introduction
Due to the intensive use of electrical energy, power transmission over long distances has become necessary in recent years.With the growth in interconnection projects, highvoltage direct current (HVDC) has become suitable as it has several advantages.Among these advantages are greater power transfer per conductor, simpler line construction as each conductor can be operated as an independent circuit, no charging current, no skin effect, lesser corona power loss and radio interference in foul weather conditions, and synchronization is not required [1].HVDC is more cost-effective than high-voltage alternating current for projects with transmission distances exceeding 280 to 300 km [2].
Because of these benefits, several DC interconnection projects have been implemented worldwide.Three different types of HVDC links are used in power transmission, named monopolar, bipolar, and homopolar HVDC transmission systems [1].One of the most significant HVDC projects worldwide is a ±1100 kV transmission line erected in China [3].In addition, several DC interconnection projects have been implemented in the Arabian Golf [4].
Generally speaking, the AC substation grounding scheme is aimed at securing safety for humans [5] and ensuring equipment protection under single-line and double-lineto-ground faults [6].Also, the DC substation grounding scheme is essential during its operation in the monopolar mode, where the current returns through the ground [7].The effectiveness of a grounding scheme is analyzed by assessing the scheme particulars including ground resistance, step voltage, touch voltage, mesh voltage, and surface potential [5] to make sure that these particulars are within the allowable limits.Ground resistance is essential for a grounding scheme in power plants and subconverter stations.A grounding system's effectiveness and safety are evaluated by ground resistance.However, measuring ground resistance is a complex problem because many factors can impact it.It is not only related to the grounding device's size, shape, and ground resistivity [8] but also influenced by the surrounding electromagnetic field, metal substances in the soil, uniformity of earth resistivity, measurement method, electrode arrangement, and seasonal changes [9].
Previously, the grounding scheme design was implemented using empirical methods by assuming uniform soil resistivity.The calculations of ground resistance and voltage gradients are executed using analytical formulas.Previous work used the boundary element method and charge simulation method [10] (point charge) to calculate the ground resistance and surface potentials in different types of soil structures [11,12].By using the uniform soil approximation, where the variation between soil layers' resistivities is enormous, calculation methods may lead to risky or over-designed earthing systems [11,13].These cases may endanger human life by exceeding restrictions on touch and step potentials.Also, it may cause damage to the electrical equipment or unjustified extra costs [14].
Therefore, a well-designed grounding scheme is a primary factor in securing the safe operation of electric systems, which plays a significant role in assuring electromagnetic affinity, human safety, and the devices' reliability.A typical grounding scheme aims to reduce construction costs and time while satisfying the permissible safety parameters.The magnitudes of the safety parameters must be below the maximum permissible values for the grounding system of a substation [9,14], which can demonstrate the application of a multi-in-one station grounding system [15].Modeling of the earth resistance after collecting data for a deep geographic mode was reported [16].For a vertical electrode with a radius of 0.01 m, the influence of the length in the range of 0.6-3 m and the cross-sectional area of the electrode on the ground resistance, current dissipation, and step and touch voltages was investigated [17] Other research works are interested in experimentally studying temperature characteristics of soils with saturated water content (water-saturated soil) and showing its effect on heating the grounding scheme without evaluating other parameters [18].Another experimental work has studied the change in the surface layer and its effect on the step voltage and earth surface potential.No attention was directed to studying the ground resistance and how it will be affected by the surface layer [19], while other research focused only on geographical sites to perform deep surveys and commissioning tests [20].Another work used the IEEE equations with the modified Harris Hawks optimization and a ground grid with unequal conductor spacing to calculate the ground resistance in two-layer soil [21].
Other manuscripts used different numerical methods to model the grounding system.One of them is the boundary element method (BEM) with a uniform soil model.The main results are grid resistance, ground potential rise, mesh voltage, and step potential [11].Another experiment measured the transient impedance using auxiliary circuits displaced vertically in the ground [22].
Research work simulated a grounding system with grid conductors to calculate the ground resistance with the CDEGS software and executed these results experimentally with a ground earth tester [15,23].Another research used the finite element method to calculate ground resistance, step potential, and touch potential by modeling the grounding system with different shapes in multilayer soil [24].Others used the simplest case of the current simulation method (point source) to model the earthing system and calculate the ground resistance, mesh voltage, and step potential in multilayer soil [25].This work was verified with experimental measurements [13,26].Also, experimental work was performed to measure the ground resistance in nonhomogeneous soils using the conventional fall-of-potential method [27,28].Different works have been published using the finite element method via the COMSOL Multiphysics package for modeling an earthing grid, which was provided with or without electrodes buried vertically.These models were simulated in uniform, multilayer, wet, and dry soil [19].In [29], the estimation of ground resistance based on the fall of the potential curve near the current electrode was studied.An introduction of a practical approach for determining the ground resistance of the grounding grid was reported [30] They calculated the ground resistance and earth surface, step, and touch potentials.Part of this work experimentally measured these parameters using the Wenner four-pole [31,32].Another one used the external charges method to model a grounding grid in a multilayer homogenous soil to calculate the earth surface potential [33].
Table 1 reports a review of previous research work on AC system grounding.The grounding of DC power system is reviewed as follows.
An HVDC system can operate in the monopolar mode in cases where the positive and negative poles are locked out.This results in the flow of ground-return current with a subsequent impact on the normal operation of the nearby AC systems.An exact model of the HVDC system, including above-ground and underground parts, was proposed to analyze the operation of the whole system.The equivalent model of the ground current loop during monopolar operation of an HVDC system was proposed for calculating the earth surface potential (ESP) near the ground electrode.The calculated ESP values were compared against those predicted by CDEGS software to verify the proposed model [7].
During monopolar operation of HVDC links or lines, ground return current (GRC) flows into the earth through the grounding neutrals causing a DC offset of the power transformers with adverse effects on the power quality of the AC system.The distribution of GRCs in the earth was simulated based on the GRC distribution model; the penetration depth and the penetration ratio of various earth models were calculated.The penetration ratio expresses the proportion of the current flowing from the earth's surface to a given depth.The simulation results were compared with those obtained using the magnetotelluric method and Wenner method (four-point method) and measured data to investigate the necessity for deep earth resistivity [34].
With the development of the UHV DC transmission lines, the grounding electrode at the converter station has a significant effect on the AC power system.The DC current in the ground results in a DC offset in the power transformers with a subsequent oscillation of the reactive power in the AC power system.With a significant DC offset, permanent damage to the transformer body may occur.The soil resistivity model of the multilayer soil around the grounding pole is simulated by using CDEGS and ANSYS packages for monopolar and bipolar schemes of the HVDC overhead line.The calculated earth potential distribution over a distance of 50 km measured from the converter station was reported.The results obtained by ANSYS are slightly smaller than those predicted by CDEGS, with a difference not exceeding 0.6 V [35].
The ground return current of a long-distance HVDC system may penetrate into AC substations through the neutrals of grounded transformers to bias them.The complex image method was used to calculate the potential distribution in multilayer earth between the ground electrodes of the HVDC system.The impact of deep earth layers on the potential distribution along the earth's surface was investigated.The obtained results showed a negligible impact from the deep earth layers on the potential distribution around the ground electrode.However, these layers have to be considered for the earth's potential rise over many kilometers away from the ground electrode.This sheds some light on the importance of investigating the impact of the ground return current from HVDC ground electrodes on AC systems provided with grounding schemes [36].
The grounding devices in a power system are aimed at ensuring its safe and stable operation.Therefore, accurate estimations of the grounding parameters for these devices are sought.The currently used calculation methods for grounding parameters are based on the complex image method with global sampling, which is not accurate enough for horizontal multilayer soils.A novel segmented sampling numerical calculation method (SSNCM) was proposed for estimating the grounding parameters in these soils.These case studies for typical grounding devices were reported using SSNCM and the CDEGS package.The results showed that the SSNCM has better adaptability for different multilayered soils [37].Changes in soil parameters, due to climatic changes, can affect the distribution of the earth potential around the HVDC ground electrodes.A model was proposed to describe how the surface layer soil moisture is related to surface layer soil resistivity under climatic changes by using image recognition technology.A soil model was also proposed to assess the surface resistivity for soils near the ground electrodes based on the relatively stable resistivity of lower soil layers.A finite element method was applied to calculate the ESP.The image recognition of surface soil resistivity was considered a low-cost, real-time, and online accurate method for assessing ESP and step voltage values around the ground electrode.The higher the resistivity of the soil, the higher the fluctuation in ESP values with subsequent potential safety hazards.To ensure personal safety, the effect of climate changes on soil surface resistivity has to be considered in the design of HVDC ground electrodes [38].
The working group WG B4.61 of CIGRE provided [3] the general guidelines for the design of ground return electrode stations for the HVDC transmission system.The grounding scheme consists of vertical rods distributed around the periphery of two concentric elliptic contours.The rods of each contour are tied together by a cable.The groups are connected together through horizontal conductors as shown in Figure 1.
horizontal multilayer soils.A novel segmented sampling numerical calculation method (SSNCM) was proposed for estimating the grounding parameters in these soils.These case studies for typical grounding devices were reported using SSNCM and the CDEGS package.The results showed that the SSNCM has better adaptability for different multilayered soils [37].
Changes in soil parameters, due to climatic changes, can affect the distribution of the earth potential around the HVDC ground electrodes.A model was proposed to describe how the surface layer soil moisture is related to surface layer soil resistivity under climatic changes by using image recognition technology.A soil model was also proposed to assess the surface resistivity for soils near the ground electrodes based on the relatively stable resistivity of lower soil layers.A finite element method was applied to calculate the ESP.The image recognition of surface soil resistivity was considered a low-cost, real-time, and online accurate method for assessing ESP and step voltage values around the ground electrode.The higher the resistivity of the soil, the higher the fluctuation in ESP values with subsequent potential safety hazards.To ensure personal safety, the effect of climate changes on soil surface resistivity has to be considered in the design of HVDC ground electrodes [38].
The working group WG B4.61 of CIGRE provided [3] the general guidelines for the design of ground return electrode stations for the HVDC transmission system.The grounding scheme consists of vertical rods distributed around the periphery of two concentric elliptic contours.The rods of each contour are tied together by a cable.The groups are connected together through horizontal conductors as shown in Figure 1.This paper is organized into seven sections.Section 2 is devoted to the research gaps in the above literature survey.Section 3 describes the problem formulation.Section 4 presents the method of analysis.Section 5 describes the experimental setup and techniques.Section 6 reports the calculated and measured results and their discussions.Section 7 summarizes the main findings in the present paper.This paper is organized into seven sections.Section 2 is devoted to the research gaps in the above literature survey.Section 3 describes the problem formulation.Section 4 presents the method of analysis.Section 5 describes the experimental setup and techniques.Section 6 reports the calculated and measured results and their discussions.Section 7 summarizes the main findings in the present paper.

Research Gaps
The above survey dictates that research work on grounding in HVDC systems in the literature is very limited when compared with that in the HVAC system.To the authors' knowledge, there is no method supported in the literature to evaluate the ground resistance for a DC converter station based on the well-known charge simulation technique [10].The equivalence concept of current conduction and electric field from the grounding scheme through the surrounding medium is applied for the first time for calculating the ground resistance of the DC converter station.This motivates the authors to apply the charge simulation technique to calculate the ground resistance for a DC converter station.Not only the ground resistance but also the potential distribution over the ground surface is calculated in the present paper for different depth values of the grounding scheme.The proposed grounding scheme for the converter station is aimed at improving the CIGRE scheme [3], where the contours positioning the vertical ground rods are no longer elliptic but circular to ensure a more uniform distribution of the current delivered from the scheme into the surrounding soil.This is in addition to the simplification of the proposed charge simulation technique (CST) due to more symmetry in the proposed scheme.

Problem Formation
This article is aimed at proposing a method for calculating the ground resistance of a grounding scheme serving an HVDC substation.The proposed method is based on the charge simulation technique that utilizes the analogy between the conduction of electric flux lines in dielectric media and current density flow in conducting media.The effects of different parameters on the calculated ground resistance value are investigated.The ground resistance of the proposed grounding scheme and ground surface potentials are measured along a distance of four meters from the scheme's center.

Geometry of the Grounding Scheme
The proposed grounding scheme is an extension to that proposed by CIGRE, as shown in Figure 1.A rodless ring is added to increase the contact surface area of the scheme with the soil.Therefore, the proposed scheme consists of a group of three concentric rings; named rodless, inner, and outer rings, which are located in a horizontal plane.The inner and outer rings are provided by a set of vertical rods uniformly distributed around the periphery of the ring.The inner ring has a radius r i , with four rods extending from it with rod length h i .The outer ring has a radius r o with four rods extending from it with rod length h o .Two perpendicular conductors in the plane of the rings divide them into four quarters, as shown in Figure 2. At each intersection point of the conductors with the inner and outer rings, there is a vertical rod.The radius of the rods and conductors are r d and r c , respectively.
Figure 2 of the proposed grounding scheme is quoted from the grounding of a wind turbine [40] and from the ground return electrode stations for an HVDC transmissions system [3].The vertical rods are tied to the outer ring electrodes.The rodless ring electrode is free from vertical rods because the outer ring electrodes shield the rodless one, so almost no current conduction to ground through the vertical rods occurs when tied to the rodless ring electrode [41].It is the same as in grounding an electric substation using a grid, where vertical rods are tied only to the outer boundary of the grid.This is because vertical rods toward the center of the grid will not carry current to the ground, if any.
distributed around the periphery of the ring.The inner ring has a radius  , with four rods extending from it with rod length ℎ .The outer ring has a radius  with four rods extending from it with rod length ℎ .Two perpendicular conductors in the plane of the rings divide them into four quarters, as shown in Figure 2. At each intersection point of the conductors with the inner and outer rings, there is a vertical rod.The radius of the rods and conductors are  and  , respectively.

Simulation Technique
To evaluate the resistance-to-ground for the proposed grounding scheme, the equivalence of current conduction from the grounding scheme through the surrounding medium and electric flux flow from the scheme in a dielectric medium is utilized.To evaluate the flux flow, the well-known CST is applied to determine the flux density at the surface of the grounding scheme composed of rings, conductors, and rods.With the flow of a fault current in the grounding scheme, a voltage V ap is developed on the scheme.The surface charges on the components of the grounding scheme are simulated by fictitious charges depending on the shape of the scheme components.
On stressing the geometry relative to the ground plane by an applied voltage V ap , the resulting surface charge on the rings is simulated with three-ring charges and on the rods and conductors by finite line charges.For each simulation charge, a boundary point is selected on the electrode surface.A set of boundary points are located on each rod and conductor being modeled as finite-lines charges extending along their axes [10,42].
The equation expressing the potential ∅ i at the ith point due to the simulation charges is expressed as Q j , j = 1, 2, 3, . . . . . ., N C , and is equal to the applied voltage V ap .
where N C is the number of simulation charges and N n is the number of boundary points, which may be equal to or greater than N C .P ij is the potential coefficient calculated at the ith boundary point due to the jth simulation charge.Images of the simulation charges are considered and the potential coefficients for the ring and finite line charges are expressed as follows.
The potential coefficient P at a point (x, y, z) due to a ring charge with r' and z', denoting ring radius and z-coordinate of the ring, is expressed as where ε is the permittivity of the soil and K(k 1 ) and K(k 2 ) are elliptic integrals of the first kind.
The potential coefficient P at a point (x, y, z), due to a finite line simulating vertical rods which extend along the Z-axis from z 1 to z 2 , is expressed as where z 1 and z 2 are the z-coordinates of the start and end of the finite line charge provided that the line charge extends alone the Z-axis.
The potential coefficient P at a point (x, y, z), due to a finite line simulating horizontal conductor which extends along the X-axis from x 1 to x 2 , is expressed as where x 1 and x 2 are the x-coordinates of the start and end of the finite line charge provided that the line charge extends alone the X-axis.
The potential coefficient P at a point (x, y, z), due to a finite line simulating horizontal conductor which extends along the Y-axis from y 1 to y 2 , is expressed as where y 1 and y 2 are the y-coordinates of the start and end of the finite line charge provided that the line charge extends alone the Y-axis.
Equation ( 1) is applied at each boundary point to satisfy a pertinent boundary condition, which is the Dirichlet boundary condition.This formulates a set of equations, whose simultaneous solution determines the unknown simulation charges.
If the number of boundary points N n exceeds the number of simulation charges N C , the method of least squares [43] is applied to formulate the set of describing equations.
To check the accuracy of the proposed charge CST, a set of N n checkpoints is chosen on the electrode surface midway between the boundary points.The potential deviation of the calculated potential from the applied voltage V ap at the checkpoints is a measure of the accuracy of the proposed CST.
The ground resistance (R g ) is obtained from [12] where C is the capacitance of the grounding scheme as expressed by where ρ is the soil resistivity and ε is the soil permittivity.ρ and ε are assumed uniform around the grounding scheme.Figure 3 describes the flow chart of the charge simulation technique for calculating the ground resistance.First, the geometry of the grounding scheme is defined, the simulation charges for the grounding scheme are located, and then the charges are calculated.After this step, the potential error is checked; if this error is acceptable within a predefined tolerance, the electrical capacitance and hence the ground resistance is calculated.If not, the number and location of the fictitious simulation charges are modified, and the calculation procedure is repeated.
where  is the soil resistivity and ε is the soil permittivity. and  are assumed uniform around the grounding scheme.Figure 3 describes the flow chart of the charge simulation technique for calculating the ground resistance.First, the geometry of the grounding scheme is defined, the simulation charges for the grounding scheme are located, and then the charges are calculated.After this step, the potential error is checked; if this error is acceptable within a predefined tolerance, the electrical capacitance and hence the ground resistance is calculated.If not, the number and location of the fictitious simulation charges are modified, and the calculation procedure is repeated.As shown in Figure 2, the total number of rods is eight and each rod is simulated by 18 finite line charges.Each of the two horizontal conductors is simulated by four finite line charges.This is in addition to the three simulation ring charges.This adds up to a total of 155 simulation charges and a total of 155 boundary points.
For a parametric study, the electrical resistivity for soil has a range from 2 to 60 Ω. m.The rods' length has a range from 0.6 m to 2.4 m.The radius of the outer ring has a range from 1.25 m to 3.0 m against 0.2 m to 1.2 m for the inner ring.

Experimental Setup and Technique
In Figure 2, three-ring electrodes are made from bare copper conductors with a diameter of 10 mm.The eight vertical rods are 19 mm in diameter and 1.2 m in length.The depth of the rings of the grounding scheme  is half a meter from ground level.The rings' diameters are one, two, and three meters.Figure 4 shows the schematic diagram for the setup used to test the DC grounding scheme.The full-wave rectifier bridge injects a direct current into the DC grounding scheme for testing purposes.As shown in Figure 2, the total number of rods is eight and each rod is simulated by 18 finite line charges.Each of the two horizontal conductors is simulated by four finite line charges.This is in addition to the three simulation ring charges.This adds up to a total of 155 simulation charges and a total of 155 boundary points.
For a parametric study, the electrical resistivity for soil has a range from 2 to 60 Ω•m.The rods' length has a range from 0.6 m to 2.4 m.The radius of the outer ring has a range from 1.25 m to 3.0 m against 0.2 m to 1.2 m for the inner ring.

Experimental Setup and Technique
In Figure 2, three-ring electrodes are made from bare copper conductors with a diameter of 10 mm.The eight vertical rods are 19 mm in diameter and 1.2 m in length.The depth of the rings of the grounding scheme H r is half a meter from ground level.The rings' diameters are one, two, and three meters.Figure 4 shows the schematic diagram for the setup used to test the DC grounding scheme.The full-wave rectifier bridge injects a direct current into the DC grounding scheme for testing purposes.Figure 5 shows the experimental setup circuit with a Variac transformer 220/(0 − 250) V that controls the AC voltage to the circuit, which is kept at 30 V.An analogue ammeter and voltmeter are used as shown in Figure 5 for the three-electrode method.The earth tester is to measure the surface potential.The three-electrode method is one of the famous methods used to measure ground resistance.This method is applied for measuring ground resistance.One electrode is connected to the DC grounding scheme.The second one (current electrode  ) is connected far away from the grounding scheme at a distance (S1 + S2), and the third one (potential electrode Pe) is connected at a distance S1 equal to 62.5% of the whole distance (S1 + S2), Figure 6.The voltage difference between the first and third rods is measured.
The ground resistance is the potential difference divided by the DC injected current.The values of 5.5 m and 3.6 m are assigned to the distance S1 and S2. Figure 5 shows the experimental setup circuit with a Variac transformer 220/(0 − 250) V that controls the AC voltage to the circuit, which is kept at 30 V.An analogue ammeter and voltmeter are used as shown in Figure 5 for the three-electrode method.The earth tester is to measure the surface potential.Figure 5 shows the experimental setup circuit with a Variac transformer 220/(0 − 250) V that controls the AC voltage to the circuit, which is kept at 30 V.An analogue ammeter and voltmeter are used as shown in Figure 5 for the three-electrode method.The earth tester is to measure the surface potential.The three-electrode method is one of the famous methods used to measure ground resistance.This method is applied for measuring ground resistance.One electrode is connected to the DC grounding scheme.The second one (current electrode  ) is connected far away from the grounding scheme at a distance (S1 + S2), and the third one (potential electrode Pe) is connected at a distance S1 equal to 62.5% of the whole distance (S1 + S2), Figure 6.The voltage difference between the first and third rods is measured.
The ground resistance is the potential difference divided by the DC injected current.The values of 5.5 m and 3.6 m are assigned to the distance S1 and S2.The three-electrode method is one of the famous methods used to measure ground resistance.This method is applied for measuring ground resistance.One electrode is connected to the DC grounding scheme.The second one (current electrode C e ) is connected far away from the grounding scheme at a distance (S1 + S2), and the third one (potential electrode Pe) is connected at a distance S1 equal to 62.5% of the whole distance (S1 + S2), Figure 6.The voltage difference between the first and third rods is measured.The ground resistance is the potential difference divided by the DC injected current.The values of 5.5 m and 3.6 m are assigned to the distance S1 and S2.

Results and Discussion
The satisfaction of the Dirichlet boundary conditions at all boundary points and the solution to the resulting equations determines the unknown simulation charges.The summation of all simulation charges and division of the sum by the applied voltage  determines the capacitance of the ground scheme.

Accuracy of Simulation Technique
Figure 7 shows the percent potential error on the rods, which reaches a maximum value of 7.8 × 10 , with checkpoints on the rods equal to 144 (8 rods;

Results and Discussion
The satisfaction of the Dirichlet boundary conditions at all boundary points and the solution to the resulting equations determines the unknown simulation charges.The summation of all simulation charges and division of the sum by the applied voltage V ap determines the capacitance of the ground scheme.
Table 2 lists the basic parameters used for the experimental setup including the values of soil permittivity and resistivity.These values represent approximately the mean values of those used for the parametric study.

Accuracy of Simulation Technique
Figure 7 shows the percent potential error on the rods, which reaches a maximum value of 7.8 × 10 −4 , with checkpoints on the rods equal to 144 (8 rods; 18 checkpoints per rod) as recorded on the horizontal axis of the figure and distributed over the eight rods in the sequence 1-18, 19-38, 37-    Figure 9 shows the percentage potential error on the rings, which reaches a maximum value of 1.1 × 10 with checkpoints numbered 153 on the intersection point of the rodless ring with the x-axis, 154 on the inner ring but shifted 120° from the x-axis, and 155 on the outer ring but shifted 240° from the x-axis.It is worth mentioning that the percent error in Figure 9 over the ring electrodes is calculated at specific checkpoints, where the error records a maximum value.Figure 8 shows the percent potential error on the conductors, which reaches a maximum value of 10 −3 with eight checkpoints on the rods numbered 145-152 with four points numbered 145-148 at a radial distance equal to r i /2 and four points with numbers 149-152 at a radial distance equal to (r i + r o )/2. Figure 8 shows the percent potential error on the conductors, which reaches a maximum value of 10 with eight checkpoints on the rods numbered 145-152 with four points numbered 145-148 at a radial distance equal to  /2 and four points with numbers 149 − 152 at a radial distance equal to ( +  )/2. Figure 9 shows the percentage potential error on the rings, which reaches a maximum value of 1.1 × 10 with checkpoints numbered 153 on the intersection point of the rodless ring with the x-axis, 154 on the inner ring but shifted 120° from the x-axis, and 155 on the outer ring but shifted 240° from the x-axis.It is worth mentioning that the percent error in Figure 9 over the ring electrodes is calculated at specific checkpoints, where the error records a maximum value.Figure 9 shows the percentage potential error on the rings, which reaches a maximum value of 1.1 × 10 −3 with checkpoints numbered 153 on the intersection point of the rodless ring with the x-axis, 154 on the inner ring but shifted 120 • from the x-axis, and 155 on the outer ring but shifted 240 • from the x-axis.It is worth mentioning that the percent error in Figure 9 over the ring electrodes is calculated at specific checkpoints, where the error records a maximum value.

Ground Resistance as Influenced by Grounding-Scheme Parameters
Figure 10 shows the effects of changing the rod diameter on the calculated parameter ground resistance values.The relationship is inversely proportional between the rod diameter and the ground resistance, the larger the rod diameter, the more the contact area of the rods with the soil.This is reflected in the decreasing  with increasing rod diameter, as indicated in Figure 10.By applying Equation ( 6), the ground resistance  is calculated to be equal to 4.8 Ω, which is very close to the measured value of 5 Ω with an error value of 4.2%.By using the CYMGRD, the ground resistance of the proposed grounding scheme records 4.7 Ω.This indicates the accuracy of the proposed simulation of the grounding scheme and enhances the confidence in the proposed method of analysis.

Ground Resistance as Influenced by Grounding-Scheme Parameters
Figure 10 shows the effects of changing the rod diameter on the calculated parameter ground resistance values.The relationship is inversely proportional between the rod diameter and the ground resistance, the larger the rod diameter, the more the contact area of the rods with the soil.This is reflected in the decreasing R g with increasing rod diameter, as indicated in Figure 10.

Ground Resistance as Influenced by Grounding-Scheme Parameters
Figure 10 shows the effects of changing the rod diameter on the calculated parameter ground resistance values.The relationship is inversely proportional between the rod diameter and the ground resistance, the larger the rod diameter, the more the contact area of the rods with the soil.This is reflected in the decreasing  with increasing rod diameter, as indicated in Figure 10.By applying Equation ( 6), the ground resistance  is calculated to be equal to 4.8 Ω, which is very close to the measured value of 5 Ω with an error value of 4.2%.By using the CYMGRD, the ground resistance of the proposed grounding scheme records 4.7 Ω.This indicates the accuracy of the proposed simulation of the grounding scheme and enhances the confidence in the proposed method of analysis.By applying Equation ( 6), the ground resistance R g is calculated to be equal to 4.8 Ω, which is very close to the measured value of 5 Ω with an error value of 4.2%.By using the CYMGRD, the ground resistance of the proposed grounding scheme records 4.7 Ω.This indicates the accuracy of the proposed simulation of the grounding scheme and enhances the confidence in the proposed method of analysis.
The measurements of the ground resistance do not fit well with the numerical calculations, which assume uniform soil resistivity.These calculations predict values that vary depending on the used formula or the adopted technique for evaluating the ground resistance [44].As an example, the ground resistance of a single rod (one of those used in the proposed grounding scheme of Figure 2) of 1.2 m in length which is flush with the surface of the soil with 33 Ω• m resistivity has values of (i) 24.96 Ω and 23.57Ω with a percent decrease of 5.58% on increasing the rod diameter from 16 mm to 22 mm by using the Rudenberg formula [45], (ii) 23.62 Ω and 22.227 Ω with a percent decrease of 5.9% by using the Dwight-Sunde formula [46], (iii) 21.95 Ω and 20.576 Ω with a percent decrease of 6.25% by using Liew-Darveniza [47], and (iv) 25.2 Ω and 23.45 Ω with a percent decrease of 6.95% by using CYMGRD software.
Figure 11 presents the influence of varying the rod length on the ground resistance with h i = h o .The relationship is inversely proportional between the rod length and the ground resistance.The same as in Figure 10, increasing the rod length increases the contact area of the grounding scheme with the surrounding soil with a subsequent decrease in R g with an increase in rod length h i or h o .
Energies 2024, 17, x FOR PEER REVIEW 15 of 22 The measurements of the ground resistance do not fit well with the numerical calculations, which assume uniform soil resistivity.These calculations predict values that vary depending on the used formula or the adopted technique for evaluating the ground resistance [44].As an example, the ground resistance of a single rod (one of those used in the proposed grounding scheme of Figure 2) of 1.2 m in length which is flush with the surface of the soil with 33 Ω. m resistivity has values of (i) 24.96 Ω and 23.57Ω with a percent decrease of 5.58% on increasing the rod diameter from 16 mm to 22 mm by using the Rudenberg formula [45], (ii) 23.62 Ω and 22.227 Ω with a percent decrease of 5.9% by using the Dwight-Sunde formula [46], (iii) 21.95 Ω and 20.576 Ω with a percent decrease of 6.25% by using Liew-Darveniza [47], and (iv) 25.2 Ω and 23.45 Ω with a percent decrease of 6.95% by using CYMGRD software.
Figure 11 presents the influence of varying the rod length on the ground resistance with ℎ = ℎ .The relationship is inversely proportional between the rod length and the ground resistance.The same as in Figure 10, increasing the rod length increases the contact area of the grounding scheme with the surrounding soil with a subsequent decrease in  with an increase in rod length ℎ or ℎ .Figure 12 shows how the ground resistance decreases with an increase in the inner radius r i in the range from 0.2 m to 1.2 m.Of course, the increase in r i of the inner ring results in increasing the contact area of the grounding scheme with the soil with a subsequent decrease in R g , as shown in Figure 11.The minimum resistance value corresponds to a radius equal to 1.2 m.
Figure 13 shows the effect of increasing the radius of the outer ring on the ground resistance.It is clear that the lowest R g value equals 3.9 Ω at a radius equal to three meters.Of course, an increase in r o of the outer ring results in increasing the contact area of the grounding scheme with the soil with a subsequent decrease in R g , as shown in Figure 13.
Figure 14 presents the calculated ground resistance values for different types of soil with resistivity values given in Table 3.Sea water is the best-conducting medium to achieve a minimum ground resistance value because sea water has the lowest resistivity value.Figure 13 shows the effect of increasing the radius of the outer ring on the ground resistance.It is clear that the lowest  value equals 3.9 Ω at a radius equal to three meters.Of course, an increase in  of the outer ring results in increasing the contact area of the grounding scheme with the soil with a subsequent decrease in  , as shown in Figure 13.   3. Sea water is the best-conducting medium to achieve a minimum ground resistance value because sea water has the lowest resistivity value.

Wet clay 8-70 33
Ground well and spring water 10-150 50  Figure 13 shows the effect of increasing the radius of the outer ring on the ground resistance.It is clear that the lowest  value equals 3.9 Ω at a radius equal to three meters.Of course, an increase in  of the outer ring results in increasing the contact area of the grounding scheme with the soil with a subsequent decrease in  , as shown in Figure 13.   3. Sea water is the best-conducting medium to achieve a minimum ground resistance value because sea water has the lowest resistivity value.Clay and sand mix.

4-300 100
Shale, slates and sandstone 10 -1000 120 Figure 15 shows the decrease in the ground resistance  with the increase in the length of the outer rod.This is the same as discussed for the effect of the increase in the rod length in Figure 11, and the effect of the increase in the inner and outer rings in Figures 12 and 13 on decreasing the ground  ; it also decreases with the increase in the length of the outer rods.Figure 15 shows the decrease in the ground resistance R g with the increase in the length of the outer rod.This is the same as discussed for the effect of the increase in the rod length in Figure 11, and the effect of the increase in the inner and outer rings in Figures 12 and 13 on decreasing the ground R g ; it also decreases with the increase in the length of the outer rods.Figure 15 shows the decrease in the ground resistance  with the increase in the length of the outer rod.This is the same as discussed for the effect of the increase in the rod length in Figure 11, and the effect of the increase in the inner and outer rings in Figures 12 and 13 on decreasing the ground  ; it also decreases with the increase in the length of the outer rods.

Surface Potential
Figure 16 shows the measured and calculated potential values at the ground surface starting from the center of the grounding scheme, provided that the DC voltage is applied at the center of the scheme.The maximum value at the center of the electrode is 12 V and decreases in the direction away from the grounding scheme and reaches 1.2 V at a distance

Surface Potential
Figure 16 shows the measured and calculated potential values at the ground surface starting from the center of the grounding scheme, provided that the DC voltage is applied at the center of the scheme.The maximum value at the center of the electrode is 12 V and decreases in the direction away from the grounding scheme and reaches 1.2 V at a distance of 4 m.The deviation of the calculated surface-potential value from those measured experimentally has an average value of 6.5%.
Figure 17 shows the calculated surface-potential distribution over the ground plane at different depths H r values of the grounding scheme into the soil.At a depth equal to 0.4 m, the maximum surface-potential value is 11.7 V, and decreases to 9.65 V on increasing the depth to 0.5 m against 8.47 V at depth of 0.6 m.This is because the increase in depth results in increasing the spacing between the simulation charges of the grounding scheme and the ground surface.The larger the spacing, the smaller the potential values at the ground surface, as shown in Figure 17.
of 4 m .The deviation of the calculated surface-potential value from those measured experimentally has an average value of 6.5%. Figure 17 shows the calculated surface-potential distribution over the ground plane at different depths  values of the grounding scheme into the soil.At a depth equal to 0.4 m , the maximum surface-potential value is 11.7 V , and decreases to 9.65 V on increasing the depth to 0.5 m against 8.47 V at depth of 0.6 m.This is because the increase in depth results in increasing the spacing between the simulation charges of the grounding scheme and the ground surface.The larger the spacing, the smaller the potential values at the ground surface, as shown in Figure 17.   Figure 17 shows the calculated surface-potential distribution over the ground plane at different depths  values of the grounding scheme into the soil.At a depth equal to 0.4 m , the maximum surface-potential value is 11.7 V , and decreases to 9.65 V on increasing the depth to 0.5 m against 8.47 V at depth of 0.6 m.This is because the increase in depth results in increasing the spacing between the simulation charges of the grounding scheme and the ground surface.The larger the spacing, the smaller the potential values at the ground surface, as shown in Figure 17.

Conclusions
This paper presents a method for calculating the ground resistance of a grounding scheme servicing a DC converter station.The method is based on the charge simulation technique and is considered a valuable approach for engineers in charge of designing and analyzing DC grounding schemes.Lower ground resistance values can be achieved by increasing the number of rods, conductors, and rings.Therefore, the present paper improves the grounding scheme proposed by CIGRE [3], seeking a decrease in the grounding resistance of the converter station in DC systems.The improvement is achieved by adding the rodless ring electrode to increase the contact surface area of the grounding scheme with the surrounding soil.
This paper studied the effects of different parameters affecting the ground resistance, such as the rod diameter, radius of the ring, etc., which express the geometry of the grounding station.The obtained results recommend a suitable depth and diameter of the grounding rods to achieve a design value of the ground resistance.Generally speaking, the larger the contact area of the grounding scheme with the soil, the lower the ground resistance of the scheme.
The calculated ground resistance agreed reasonably well with that obtained experimentally with an error value of 4.2%.The calculated surface-potential values agreed with those measured experimentally with an average percent error not exceeding 6.5%.
The maximum percent potential error on the rods, conductors, and ring electrodes recorded values of 7.8 × 10 −4 , 10 −3 , and 1.1 × 10 −3 , respectively.This is reflected in the accuracy of the adopted charge simulation technique.
An increase in rod diameter from 16 mm to 22 mm resulted in a decrease in the ground resistance by about 13%, and an increase in rod length from 0.6 m to 2.4 m resulted in a decrease in the ground resistance by about 22%.
An increase in the radius of the inner ring from 0.2 m to 1.2 m resulted in a decrease in the ground resistance by 9%, and an increase in the radius of the outer ring from 1.25 m to 3 m resulted in a decrease in the ground resistance by 25%.

Figure 3 .
Figure 3. Flow chart for ground resistance calculation using the CST.

Figure 3 .
Figure 3. Flow chart for ground resistance calculation using the CST.

Figure 6 .
Figure 6.Schematic diagram for the three-electrode method.
Figure 7 shows the percent potential error on the rods, which reaches a maximum value of 7.8 × 10 , with checkpoints on the rods equal to 144 (8 rods; 18 checkpoints per rod) as recorded on the horizontal axis of the figure and distributed over the eight rods in the sequence 1 − 18, 19 − 38, 37 − 54, 55 − 72, 73 − 90, 91 − 108, 109 − 126 and 127 − 144.The checkpoints on each rod start from top to bottom.

Figure 6 .
Figure 6.Schematic diagram for the three-electrode method.

Figure 8
Figure8shows the percent potential error on the conductors, which reaches a maximum value of 10 with eight checkpoints on the rods numbered 145-152 with four points numbered 145-148 at a radial distance equal to  /2 and four points with numbers 149 − 152 at a radial distance equal to ( +  )/2.

Figure 7 .
Figure 7. Potential error at rods at the different checkpoints.(h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 8 .
Figure 8. Potential error at the conductors at the different checkpoints.(h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 9 .
Figure 9. Potential error at the rings at the different checkpoints.(h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 10 .
Figure 10.Effect of rods' diameter on the ground resistance.(h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r c = 0.005 m, H r = 0.5 m).

Figure 12
Figure12shows how the ground resistance decreases with an increase in the inner radius  in the range from 0.2 m to 1.2 m.Of course, the increase in  of the inner ring results in increasing the contact area of the grounding scheme with the soil with a subsequent decrease in  , as shown in Figure11.The minimum resistance value corresponds to a radius equal to 1.2 m.

Figure 11 .
Figure 11.Effect of the rods' length on the ground resistance.(h i = h o , r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 14
Figure 14 presents the calculated ground resistance values for different types of soil with resistivity values given in Table3.Sea water is the best-conducting medium to achieve a minimum ground resistance value because sea water has the lowest resistivity value.

Figure 12 .
Figure 12.Effect of varying the inner ring radius on the ground resistance.(h i = h o = 1.2 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 14
Figure 14 presents the calculated ground resistance values for different types of soil with resistivity values given in Table3.Sea water is the best-conducting medium to achieve a minimum ground resistance value because sea water has the lowest resistivity value.

Figure 13 .
Figure 13.Effect of varying the outer ring radius on the ground resistance.(h i = h o = 1.2 m, r i = 1 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 14 .
Figure 14.Ground resistance for different soil types (h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 15 .
Figure 15.Ground resistance value versus the length of the rods for outer ring (h i = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 16 .
Figure 16.Measured and calculated surface potential versus the distance from the center of the grounding scheme (V ap = 30 V, h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m, H r = 0.5 m).

Figure 17 .
Figure 17.Surface-potential values over the ground plane at different depth  values of the scheme rings.( = 30 V, ℎ = ℎ = 1.2 m,  = 1 m,  = 1.5 m,  = 0.0095 m,  = 0.005 m).Figure 17.Surface-potential values over the ground plane at different depth H r values of the scheme rings.V ap = 30 V, h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m .

Figure 17 .
Figure 17.Surface-potential values over the ground plane at different depth  values of the scheme rings.( = 30 V, ℎ = ℎ = 1.2 m,  = 1 m,  = 1.5 m,  = 0.0095 m,  = 0.005 m).Figure 17.Surface-potential values over the ground plane at different depth H r values of the scheme rings.V ap = 30 V, h i = h o = 1.2 m, r i = 1 m, r o = 1.5 m, r d = 0.0095 m, r c = 0.005 m .

Table 1 .
A brief literature review of grounding schemes in AC power systems.

Table 2
lists the basic parameters used for the experimental setup including the values of soil permittivity and resistivity.These values represent approximately the mean values of those used for the parametric study.

Table 2 .
Basic design parameters.

Table 2 .
Basic design parameters.