International Journal of Antennas and Propagation
Volume 2016 (2016), Article ID 8601497, 8 pages
http://dx.doi.org/10.1155/2016/8601497

Research Article

An Optimal Electric Dipole Antenna Model and Its Field Propagation

Yidong Xu,¹ Lili Guo,¹ Wei Xue,¹ Korochentsev Vladimir,² and Junwei Qi¹

¹College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
²School of Engineering, Far Eastern Federal University, Vladivostok 690950, Russia

Received 4 February 2016; Revised 5 March 2016; Accepted 17 April 2016

Academic Editor: Christoph F. Mecklenbräuker

Copyright © 2016 Yidong Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An optimal electric dipole antennas model is presented and analyzed, based on the hemispherical grounding equivalent model and the superposition principle. The paper also presents a full-wave electromagnetic simulation for the electromagnetic field propagation in layered conducting medium, which is excited by the horizontal electric dipole antennas. Optimum frequency for field transmission in different depth is carried out and verified by the experimental results in comparison with previously reported simulation over a digital wireless Through-The-Earth communication system. The experimental results demonstrate that the dipole antenna grounding impedance and the output power can be efficiently reduced by using the optimal electric dipole antenna model and operating at the optimum frequency in a vertical transmission depth up to 300 m beneath the surface of the earth.

1. Introduction

Electric dipole antennas have been widely used in Through-The-Earth (TTE) communication via a current injected into the earth at the power stage [1–3]. The Coal Mine Security Standard has many restrictions on the permissible transmission power of TTE communication systems. There are many methods which have been proposed and investigated to reduce the output power by changing the shape and material, increasing the interelectrode distance [4, 5], or using only one electrode with different kinds of material, such as parabolt, copper braid extended aluminum foil, and steel rod [3]. However, these methods have not discussed the grounding impedance model of several electrodes in-array and still cannot reduce the electrodes grounding impedance to a level less than 300 Ω. Many reports take the contact impedance between the electrode and the earth into account [2] and suggest lowering it with resistance reducing material but have not given the quantitative relationship between them. Some models have been proposed and investigated the dipole antenna in two- or three-layered conducting medium [6–8]. However, most of them discuss the propagation from surface to surface or for far field transmission and have not given the optimal operating frequency to reduce the TTE communication power in a stratified dissipative medium vertically. In [9], the author gives the electromagnetic wave radiation propagate in four-layered medium model in the VHF/UHF bands which is a great limitation for TTE communication. References [7, 10] have proposed the models of extremely low frequency electromagnetic field generated by submerged horizontal electric dipole (HED) in three-layer medium, where the field point and the source point are in the same layer. These models cannot be applied in practical TTE communication because the source point and the field point are not in the same layer.

To reduce the output power, a model for the parallel grounding impedance of several HED antennas arrayed in-line is presented and the quantitative relationship between the resistance reducing material and the impedance is proposed in this paper. It is found that the grounding impedance model of several electrodes in-array has good match with the experimental results, and the optimal electric dipole antenna can reduce the grounding impedance to less than 5 Ω. In order to further lower the output power and to improve communication quality, a full-wave electromagnetic simulation model for the electromagnetic field propagation in layered conducting medium is proposed, which is excited by the HED antennas. An optimum frequency for field transmission at different propagation depth is carried out by this model. The receiver can get a maximum field when the TTE communication transmitter works at this optimum frequency, which is numerically and experimentally verified over a digital wireless system.

2. Modeling and Optimization of an Electrical Dipole Antenna

2.1. HED Antenna Grounding Model

Electric dipole antennas are used to inject electrical power into the earth by establishing galvanic coupling between the earth and the power stage of the transmitter. The frequency-dependent grounding impedance includes the wire impedance, the contact impedance between the wire and the electrode, the impedance of the electrode, the earth current divergence impedance, and the contact impedance between the earth and the electrode [2, 11, 12]. However, when the grounding impedance is in the frequency range of lower than 50 kHz, the grounding electrode exhibits the grounding impedance close to the grounding resistance [13]. In this paper, we investigate the earth current divergence impedance.

For a horizontal grounding electrode, with a length , a diameter , and an embedded depth , six possible earth current divergence impedance models are presented in Table 1 [12].

Table 1: Horizontal grounding electrode grounding resistance models.

Model
Formula

Rudenberg, Zingraff

Tagg, Dwight

Sunde, Schwarz

Neutral-point potential method

Average potential method ()

Neutral-point potential method ()

For the sake of simplicity, the resistance of a horizontal grounding electrode with an embedded depth of zero can be expressed by the neutral-point potential method. This states that increasing the earth conductivity and the size of electrode can reduce its current divergence resistance. By placing several electrodes in parallel as one new electrode, the grounding impedance can be further reduced. The new grounding impedance can be expressed as follows:

where is the grounding impedance of a single electrode; is the number of electrodes; is the grounding electrode utilization factor, which is determined by the shape and number of electrodes and their relative positions.

The parallel resistance can be similarly calculated by introducing a hemispherical grounding equivalent model [12]. For a horizontal grounding electrode buried at the subsurface level, its equivalent radius expression can be obtained based on image theory and the neutral-point potential method as

For parallel horizontal grounding electrodes arrayed in-line, the distance between two adjacent electrodes is . In order to simplify the calculation, a radius ratio is defined. The injection current over each electrode is with the the electric potential . The electric potential of the th electrode as a reference electrode is given by , based on the hemispherical grounding equivalent model. The electric potential provided by the th electrode at the th electrode can be expressed as . Then, we obtain the electric potential by the superposition principle. By setting each electrode as reference, we have

When , the expressions of grounding electrode utilization factor with argument radius ratio are obtained and shown in Table 2. The variations of with are shown in Figure 1.

Table 2: Functional relationships between and .

Numbers of electrodes
Utilization factor

2

3

4

5

6

Figure 1: Variation of utilization factor with under dipole antennas.

From Figure 1, it can be seen that the grounding electrode utilization factor decreases as the parameter increases for the same number of electrodes . The main reason is that the shielding effect increases as the contiguous electrodes distance decreases, restricting the current injection into the earth. When there is no shielding effect between electrodes, becomes infinitely small and tends towards 1.

The grounding impedance can be reduced when the grounding electrodes are surrounded by high conductivity resistance reducing material [14].

Semicylinder pits with length and radius , where , are filled with resistance reducing material with conductivity . Metal horizontal grounding electrodes with length and radius are buried at the pit center. The transitional area is the area between the resistance reducing material and the surrounding soil which depends on the seepage property of the material, and it expands the pits to an equivalent radius . The soil moisture, temperature, and soil gap have a significant influence on the resistance reducing material seepage property, meaning that is variable but always larger than , as shown in Figure 2.

Figure 2: Grounding electrodes surrounded by resistance reducing material with high conductivity.

The grounding resistance is equivalent to a series connection circuit of , , , and , where is the contact resistance between the resistance reducing material and the electrode, is the resistance reducing material current divergence impedance, is the contact resistance between the resistance reducing material and the earth, and is the earth current divergence impedance:

can be obtained by using

It can be simplified by considering the semicylinder for the resistance reducing material and the current divergence impedance. of a horizontal semicylindrical electrode can be derived from a model of a vertical electrode grounding impedance and the principle of mirror image; that is,

It can be seen that the resistance of the horizontal semicylindrical grounding electrode buried under the subsurface of the earth is equal to that derived by the midpoint potential method. The resistance reducing material has a gel-like property that keeps it in a certain shape and makes good electrical contact with the medium surface. Therefore, the impedances and can be negligible.

For demonstration purposes, it can be assumed that the conductivity of the moist soil is 0.05 S/m, and the conductivity of the resistance reducing material is 2 S/m. The length of each grounding electrode is 1.2 m with radius  m and a grounding hemisphere equivalent radius of  m. The earth current divergence impedance is 18.6 Ω when the horizontal grounding electrodes are buried at the subsurface of the earth directly without the resistance reducing material. If the horizontal grounding electrodes are buried in semicylindrical pits of radius  m that are filled with resistance reducing material, which gives their equivalent radius , in this case, we obtain  Ω and  Ω, meaning that the earth current divergence impedance  Ω. It can be seen that the grounding impedance would be effectively reduced by filling the pit with the resistance reducing material. In order to reduce the grounding impedance to 9.75 Ω by using the grounding method without the resistance reducing material, there should be at least two electrodes placed parallelly with 3.4 m electrode interval, which would be an inefficient use of space and electrode material.

2.2. Grounding Impedance Physical Experiment Simulation and Verification

One set of experiments has been carried out on a sandy beach in Harbin. A stainless steel tube has been used as the test electrode, with length 1.2 m and diameter 0.14 m. An Agilent U1733C LCR meter has been used to measure the impedance.

There are two configurations for placing the antennas: several electrodes are placed in parallel and the distance between any two contiguous electrodes is (configuration A), and the electrodes are placed in parallel in a head to tail configuration (configuration B).

The measurement results of the electrodes placed in configuration A and configuration B are shown in Table 3. According to the utilization factor expression in Table 2, the radius ratio can be written as when . The radius ratio is 0.46 with the electrode interval 0.5 m in configuration A when  Ω and  Ω. For the electrode intervals m, 1.5 m, and 2 m, the radius ratios are , , and . The radius ratio is 0.25 in configuration B when  Ω. It is when the interval is 1.7 m. Substituting radius ratios into utilization factor models in Table 2, the theoretical utilization factor can be derived. We can obtain the practical utilization factor based on (1).

Table 3: Grounding impedance of electrodes and comparison between practical and theoretical utilization factor.

Electrodes configuration
Interval
Number of electrodes

Measured impedance

Utilization factor deviation

A
0.5
1
ł
21.00
ł
ł
ł

0.5
2
ł
15.30
ł
ł
ł

0.5
3
0.46
11.82
0.59
0.58
2.8%

0.5
4
0.46
10.65
0.49
0.51
−3.7%

0.5
5
0.46
9.36
0.45
0.47
−4.6%

1.0
2
0.23
13.14
0.80
0.81
−1.2%

1.0
3
0.23
10.20
0.69
0.72
−4.2%

1.5
2
0.15
12.83
0.82
0.86
−4.7%

2.0
2
0.12
12.39
0.85
0.89
−4.5%

B
1.7
1
ł
21.0
ł
ł
ł

1.7
2
ł
13.11
ł
ł
ł

1.7
3
0.25
9.57
0.73
0.71
2.8%

3.4
2
0.13
11.94
0.88
0.89
−1.1%

The theoretical utilization factor matches the practical utilization factor well within a deviation range. Compared with the electrodes placed in parallel, the grounding impedance of the electrodes arranged in parallel in a head to tail configuration is much smaller when they have the same gap. The arrangement of electrodes in parallel in a head to tail configuration increases the distance between each electrode because of its length, causing different grounding impedance. For instance, the interval in configuration A is 0.5 m with the gap 0.5 m but 1.7 m in configuration B, accounting 1.2 m in length of the electrode itself, and the grounding impedance is 11.82 Ω in configuration A and 9.57 Ω in configuration B. It has been proven that an impedance model of the horizontal grounding electrodes buried at the subsurface of the earth could be similarly calculated using the hemisphere grounding model. The electrodes arranged in parallel in a head to tail configuration are more suitable for applications such as in mines where there is a width restriction.

The experiment to optimize the antenna has been investigated at Xinqiang Mine in Qitaihe, Heilongjiang. The electrodes are buried in semicylindrical pits filled with resistance reducing material under the subsurface of the earth and in the mine and are parallelly arranged in a head to tail configuration. Each side of the antenna is composed of three electrodes, which are 1.3 m in length, 0.07 m in diameter, and at an interval of 1.5 m. The antennas are buried in pits that are 1.3 m in length and 0.3 m in diameter. The distance between two sides of the antenna on the ground is  m but  m within the mine. Moist clay is the predominant soil on the ground, which has a conductivity within the range of 0.01~0.05 S/m. The conductivity in mine is in the range of 0.005~0.02 S/m, due to the fact that the predominant soil is rock. The impedance of each antenna is measured by an Agilent U1733C LCR meter, and the results are shown in Table 4.

Table 4: Grounding impedance of electrodes on the ground and in the mine.

Frequency
kHz
Grounding impedance of electrodes on the ground
Grounding impedance of electrodes in the mine

0.1
5.4
12.7

1
4.5
14.0

10
4.4
13.9

As seen in Table 4, the impedance of the dipole antenna on the ground is about 5 Ω, indicating that the grounding impedance of each side of the antenna is less than 2.5 Ω. The conductivity in the mine tunnel is lower than that on the ground as the geological material within the mine is mainly rock. The impedance of the dipole antenna in the mine is 13 Ω and 6.5 Ω at each side of the antenna.

The output current ability of the TTE communication system and its range can be improved by placing the antennas in parallel in a head to tail configuration and filling the pits with resistance reducing material. This grounding model meets the requirement to reduce the underground communication systems’ power.

3. Propagation of Dipole Antennas in a Semi-Infinite Layered Conducting Medium

3.1. Simulation Model of Dipole Antenna in Layered Conducting Medium

A numerical calculation of the electric field has been employed successfully in many applications, but it cannot describe the propagation successfully for complicated structures, such as layered conducting mediums [15]. To investigate the propagation of dipole antenna in layered conducting medium, the full-wave electromagnetic simulation model with the CST STUDIO SUIT electromagnetic tool is set up to accurately investigate the propagation.

According to a mine regional stratigraphic structural table, the Quaternary, Cretaceous, and Jurassic compositions are the main structural components of mines. The surface is a Quaternary system composed of alluvial sand. The middle layer is a Cretaceous system dominated by purple sandstone and grey-greenish mudstone on the upper layer and purplish-red conglomerate with mudstone in the lower layer. The bottom layer is a Jurassic system, with grey-greenish sandstone, conglomerate, mudstone, and marl on the upper level, volcanic rock, rhyolite, and mudstone in the middle level, and sandstone and conglomerate in the lower level. In order to reduce the influence of boundary conditions, the thickness of the air on the ground is set at 50 m, and the simulation area is 2000 m in length and width. The electrical parameters of these strata are shown in Table 5.

Table 5: The electrical parameters of strata.

Strata
Attribute
Conductivity
Thickness

S/m
m

Quaternary
Alluvial sand, clay
0.01
20

Cretaceous
Sandstone
0.006
180

Coal
0.01
2

Conglomerate
0.01
200

Coal
0.01
3

Conglomerate mixed with mudstone
0.008
160

Jurassic
Marl
0.003
260

Volcanic rock
0.001
240

Conglomerate
0.01
200

The full-wave electromagnetic simulation results with the CST STUDIO SUIT electromagnetic tool are presented with a moment  A·m on the dipole antenna. To understand the full-wave propagation characteristics of the dipole in the layered conducting earth, the dipole antenna working at operating frequency 1.5 kHz is simulated and shown in Figure 3. It is found that there is a huge change at the interface between every two mediums, making the propagation complicated. As the points investigated are far enough from the boundary, the boundary condition has little influence on the simulation results. The full-wave electromagnetic simulation results of the field point at positions of −200 m, −300 m, and −600 m are shown in Table 6.

Table 6: Electric field of the field point.

Frequency
Field strength at
Field strength at
Field strength at

kHz
−200 m μV/m
−300 m μV/m
−600 m μV/m

0.1
174.0
56.0
8.8

0.3
187.5
65.2
12.2

1
204.3
66.6
9.0

1.5
204.5
61.9
6.3

2
199.7
56.2
4.5

3
183.7
45.3
2.4

5
147.4
29.0
0.84

10
81.8
10.4
0.096

20
28.2
1.9
0.003

30
11.2
0.5
0.0002

kHz
1.5
1.0
0.3

−3 dB bandwidth kHz
5.0
3.0
1.0

Figure 3: Propagation of the dipole in layered conducting earth at operating frequency 1.5 kHz.

As shown in the simulation results, there is an optimum frequency at each propagation depth. For the field point at position −200 m, is 1.5 kHz with a −3 dB bandwidth 5 kHz. When the field point is located at position −300 m,  kHz and the −3 dB bandwidth decreases to 3 kHz. For the field point at −600 m, is kHz and the −3 dB bandwidth is 1 kHz. In this instance, and the −3 dB bandwidth decrease as the distance between the source point and the field point increases. For extremely low frequency, the electric field at the field point shows characteristic of a static field, and the bandwidth lower than is within the  dB bandwidth. A further simulation result has concluded that when the communication distance is within a range of 200~600 m vertically and the conductivity of the earth is within 0.001~0.05 S/m, the most favorable operating frequencies are below 10 kHz.

3.2. Field Experiment for the Propagation of a Dipole Antenna

As shown in Figure 4, a set of experiments has been carried out at Xinqiang Mine in Heilongjiang Province. The test point at the mine is 300 m from the test point on the ground vertically and 300 m in straight-line distance from the wellhead horizontally. The test point is not connected with the lift shaft directly so that the electromagnetic energy cannot propagate along that. The interelectrode distance for the antenna in the mine is 80 m and for the antenna on the ground is 40 m. The dipole antenna on the ground is parallel to that in the mine, with their azimuth at .

Figure 4: TTE communication experiment in the mine.

For the propagation measurement, the signal source, power amplifier, low noise linear amplitude to amplify the voltage on the dipole antenna, and scope to store the wave data are used. The power signal output by the power amplifier is loaded onto the dipole antenna, with an operating frequency in the range of 1~15 kHz. The power stage and the dipole antenna are connected by a multifiber copper wire with a length in the range of 20~40 m. There is an additional inductive-resistive component on the total impedance. Therefore, the power amplifier output should be regulated so that the moment on the antenna is a fixed value of  A·m.

The full-wave electromagnetic simulation model with the CST STUDIO SUIT electromagnetic tool is set, where the thickness of the air on the ground is 50 m, and the simulation area is 2000 m in length and width. The upper layer of the earth is alluvial sand and sandstone with the thickness of about 180 m and conductivity 0.007 S/m. The next layer is conglomerate and coal with the thickness of about 200 m and conductivity 0.005 S/m. The bottom layer is mudstone with the thickness of about 160 m and conductivity 0.008 S/m. Figure 5 shows a comparison of the relative magnitude of the full-wave electromagnetic simulation results and field experimental results.

Figure 5: The full-wave electromagnetic simulation results and the field experimental results curve.

From Figure 5, the magnitude decreases as the frequency increases, with a slope of about −1.2 dB/kHz for the frequency higher than 3 kHz, but for lower frequency the slope tends to be stable. It can be seen from Table 7 that most of the linear errors between the simulation results and the experimental results are less than . Four of the linear errors deviate from normal value because of the tone interference at certain frequency. But the experimental results are close to the simulation results in Figure 5, indicating that the electromagnetic tool can represent the propagation of a dipole antenna in a layered conducting medium accurately. For the mine studied, the optimal frequency at a depth of −300 m is 1 kHz.

Table 7: The full region electromagnetic simulation results and the field experimental results.

Frequency
Simulation results
Experimental results
Linear error

kHz
dB μV/m
dB μV/m
%

0.3
39.72
ł
ł

1
40.01
40.01
0.00

2
39.65
39.72
0.81

3
38.67
39.02
4.11

4
37.48
37.86
4.47

5
36.21
36.12
−1.03

6
34.87
34.95
0.93

7
33.58
32.14
−15.28

8
32.29
32.18
−1.26

9
31.04
30.54
−5.59

10
29.80
29.85
0.58

11
28.56
28.80
2.80

12
27.37
26.87
−5.59

13
26.19
27.48
16.01

14
25.09
25.80
8.52

15
23.94
22.66
−13.70

4. TTE Communication for Mine

A TTE communication system for mines has been designed, which operates at frequencies of 2 kHz or 8 kHz, with a message data rate of 25 bps and 50 bps for digital voice information. It works in the half duplex, using 2DPSK modulation and demodulation technology. A high-performance speech compression MDV chip is used, with an ultralow voice coding rate of 600 bps. Digital voice data are error correction coded, stored in a buffer, and then transmitted with a low data rate.

The TTE communication experiment has been performed in Xinqiang Coal Mine, Heilongjiang Province. A dipole antenna is placed in the pits with radius 0.2 m filled with a resistance reducing material. Each side of the antenna includes 3 electrodes that are 1.2 m in length and 0.7 m in diameter, parallelly arranged in a head to tail configuration, with a gap of 1 m. The antenna impedance on the ground is reduced to 8 Ω, with a 40 m electrode interval. The impedance of the antenna in the mine is reduced to 13 Ω, with an 80 m electrode interval. Thus, the grounding impedance is effectively reduced compared with that in [3]. For an efficient transmitter, the transformer turn ratios 1.4 : 1 and 1.8 : 1 are fixed in order to match the ideal power stage load value 4 Ω [3]. The vertical distance from the ground to the tunnel of the mine is 300 m. The results show that the short message and digital voice can be effectively received at 8 kHz, when the output power is greater than 180 W. However, the system output power is only 30 W at the operating frequency 2 kHz with the same communication performance. Thus, the dipole antenna can be optimized using the optimal electric dipole antenna model and the communication output power can be reduced by transporting information at the carrier frequency , which further proves the results in this paper.

5. Conclusions

In this paper, a grounding impedance model is presented for several horizontal electrodes placed in parallel. The analysis and contrast experiments have shown that the impedance of the dipole antenna on the ground is reduced to 5 Ω and 13 Ω in the mine. Compared with the grounding impedance 300 Ω in [3], the optimal grounding impedance model can effectively reduce the grounding impedance.

The full-wave electromagnetic simulation model with the CST STUDIO SUIT electromagnetic tool can describe the propagation property of dipole in layered conducting medium and verified by field experiment. There is a single operating frequency for each field point at each depth of the layered earth that can achieve maximum field strength. The in the full-wave simulation model at position −200 m is 1.5 kHz, 1 kHz at −300 m, and  kHz at −600 m. The operating frequency decreases as the communication distance or the earth conductivity increases. Thus, TTE communication is suitable for short range communication in the vertical direction.

It is shown that the TTE communication method could be successfully used for communication across a distance of  m vertically from the ground into the mine, with an output power reduced from 180 W at 8 kHz to 30 W at 2 kHz and an information data rate of 50 bps. The communication system can be placed in a fixed position and can be used for transmitting low data rate information, such as monitoring information or a distress signal in emergency situation.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This research has been supported by International Science & Technology Cooperation Program of China (2014 DFR10240), Harbin Science and Technology Research Projects (2013AE1BE003), Science Foundation of Heilongjiang Province QC2015075, and the Fundamental Research Funds for the Central Universities (China) [HEUCF1508].

References

H. Wang, K. Zheng, K. Yang, and Y. Ma, “Electromagnetic field in air produced by a horizontal magnetic dipole immersed in sea: theoretical analysis and experimental results,” IEEE Transactions on Antennas& Propagation, vol. 62, no. 9, pp. 4647–4655, 2014.

V. Bataller, A. Muñoz, P. M. Gaudó, A. Mediano, J. A. Cuchí, and J. L. Villarroel, “Earth impedance model for through-the-earth communication applications with electrodes,” Radio Science, vol. 45, no. 6, Article ID RS6015, pp. 2970–2980, 2010.

V. Bataller, A. Muñoz, P. M. Gaudó, A. Mediano, J. A. Cuchí, and J. L. Villarroel, “Improving medium access in through-the-earth VLF-LF communications,” Journal of Communications, vol. 4, no. 4, pp. 284–294, 2009.

A. A. Abdou, A. Shaw, A. Mason, A. Al-Shamma'a, J. Cullen, and S. Wylie, “Electromagnetic (EM) wave propagation for the development of an underwater Wireless Sensor Network (WSN),” in Proceedings of the IEEE Sensors, pp. 1571–1574, IEEE, Limerick, Ireland, October 2011.

G. Adw, Channel Characterisation and System Design for Sub Surface Communications, University of Leeds, Leeds, UK, 2003.

M. B. Kraichman, Handbook of Electromagnetic Propagation in Conducting Media, University of Michigan Library, 1970.

U. M. Cella, R. Johnstone, and N. Shuley, “Electromagnetic wave wireless communication in shallow water coastal environment: theoretical analysis and experimental results,” in Proceedings of the 4th ACM International Workshop on UnderWater Networks (WUWNet '09), pp. 1–8, Berkeley, Calif, USA, November 2009.

P. R. Bannister, “Quasi-static fields of dipole antennas at the earth’s surface,” Radio Science, vol. 1, no. 11, pp. 1321–1332, 1966.

Y. H. Xu, K. Li, and L. Liu, “Electromagnetic field of a horizontal electric dipole in the presence of a four-layered region,” Progress in Electromagnetics Research, vol. 81, pp. 371–391, 2008.

C. Liu, L.-G. Zheng, and Y.-P. Li, “Study of ELF electromagnetic fields from a submerged horizontal electric dipole positioned in a sea of finite depth,” in Proceedings of the 3rd IEEE International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications (MAPE '09), pp. 152–157, Beijing, China, October 2009.

S. Bourg, B. Sacepe, and T. Debu, “Deep earth electrodes in highly resistive ground: frequency behaviour,” in Proceedings of the IEEE International Symposium on Electromagnetic Compatibility, pp. 584–589, Atlanta, Ga, USA, August 1995.

J. He and R. Zeng, Power System Grounding, Science Press, 2007.

K. S. Lee and Y. K. Chung, “Apparatus and method for calculating length of carbon grounding electrode module based on two-layered distributed constant circuit,” Patent WO 2012050338 A3, 2013

Z. Feng, L. Lu, and J. Feng, “Research on reducing grounding resistance of transmission line tower grounding grid,” in Proceedings of the 2nd Annual Conference on Electrical and Control Engineering (ICECE '11), pp. 1216–1219, Yichang, China, September 2011.

K. Li, Fields and Waves in Layered Media, Zhejiang University Press, Zhejiang, China, 2010.