# An Optimal Electric Dipole Antenna Model and Its Field Propagation

Authors: Yidong Xu; Lili Guo; Wei Xue; Korochentsev Vladimir; Junwei Qi
Journal: International Journal of Antennas and Propagation
 (2016)

---

## 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.


---

## Body

## 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 lengthL, a diameter d, and an embedded depth h, 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
R
=
1
σ
1
2
π
L
ln
⁡
2
L
d
1
+
ln
⁡
(
L
/
2
h
)
ln
⁡
(
2
L
/
d
)
Tagg, Dwight
R
=
1
σ
1
2
π
L
ln
⁡
4
L
d
+
ln
⁡
L
h
-
2
+
2
h
L
-
h
2
L
2
+
h
4
8
L
4
Sunde, Schwarz
R
=
1
σ
1
π
L
ln
⁡
2
L
d
h
-
1
Neutral-point potential method
R
=
1
σ
1
π
L
ln
⁡
2
L
d
h
Average potential method (h
=
0)
R
=
1
σ
1
π
L
ln
⁡
4
L
d
-
1
Neutral-point potential method (h
=
0)
R
=
1
σ
1
π
L
ln
⁡
2
L
dFor 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 impedanceR
d can be expressed as follows:(1)
R
d
=
R
L
n
η
,where R
L is the grounding impedance of a single electrode; n 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 resistanceR
d 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(2)
r
=
L
2
ln
⁡
2
L
/
d
.Forn parallel horizontal grounding electrodes arrayed in-line, the distance between two adjacent electrodes is D. In order to simplify the calculation, a radius ratio α
=
r
/
D is defined. The injection current over each electrode is I
1
,
I
2
,
…
,
I
i
,
…
,
I
n with the the electric potential V. The electric potential of the ith electrode as a reference electrode is given by Φ
i
=
1
/
σ
1
/
2
π
r
I
i, based on the hemispherical grounding equivalent model. The electric potential provided by the i
±
kth electrode at the ith electrode can be expressed as Φ
i
±
k
=
1
/
σ
1
/
2
π
r
α
/
k
I
i
±
k. Then, we obtain the electric potential V
=
∑
i
=
1
n
Φ
i by the superposition principle. By setting each electrode as reference, we have(3)
1
σ
1
2
π
r
1
α
α
2
⋯
α
n
-
1
α
1
α
⋯
α
n
-
2
⋱
α
i
-
1
⋯
1
⋯
α
n
-
i
⋱
α
n
-
1
α
n
-
2
α
n
-
3
⋯
1
I
1
I
2
⋮
I
i
⋮
I
n
=
V
V
⋮
V
⋮
V
R
d
=
V
∑
i
=
1
n
I
i
=
R
L
n
η
α
=
r
D
.Whenn
=
2,3
,
4,5
,
6, the expressions of grounding electrode utilization factor η with argument radius ratio α are obtained and shown in Table 2. The variations of η with α
-
1 are shown in Figure 1.Table 2
Functional relationships betweenη and α.
Numbers of electrodes
Utilization factorη
2
1
1
+
α
3
6
-
7
α
6
+
3
α
-
12
α
2
4
12
-
10
α
12
+
16
α
-
23
α
2
5
360
-
654
α
+
241
α
2
360
+
270
α
-
1495
α
2
+
690
α
3
6
2160
-
2952
α
+
555
α
2
2160
+
3312
α
-
8439
α
2
+
2047
α
3Figure 1
Variation of utilization factorη with α
-
1 under n dipole antennas.From Figure1, it can be seen that the grounding electrode utilization factor η decreases as the parameter α increases for the same number of electrodes n. 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 lengthL and radius b, where b
≪
L, are filled with resistance reducing material with conductivity σ
s. Metal horizontal grounding electrodes with length L and radius a
=
d
/
2 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 b
m. The soil moisture, temperature, and soil gap have a significant influence on the resistance reducing material seepage property, meaning that b
m is variable but always larger than b, as shown in Figure 2.Figure 2
Grounding electrodes surrounded by resistance reducing material with high conductivity.The grounding resistanceR
electrode is equivalent to a series connection circuit of R
c
1, R
e
1, R
c
2, and R
e
2, where R
c
1 is the contact resistance between the resistance reducing material and the electrode, R
e
1 is the resistance reducing material current divergence impedance, R
c
2 is the contact resistance between the resistance reducing material and the earth, and R
e
2 is the earth current divergence impedance:(4)
R
electrode
=
R
c
1
+
R
e
1
+
R
c
2
+
R
e
2
.R
e
1 can be obtained by using(5)
R
e
1
=
1
σ
s
1
π
L
ln
⁡
b
m
a
.It can be simplified by considering the semicylinder for the resistance reducing material and the current divergence impedance.R
e
2 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,(6)
R
e
2
=
1
σ
1
π
L
ln
⁡
L
b
m
.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 impedancesR
c
1 and R
c
2 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 electrodeL is 1.2 m with radius a
=
0.035 m and a grounding hemisphere equivalent radius of r
=
0.17 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 b
=
0.2 m that are filled with resistance reducing material, which gives their equivalent radius b
m
>
b, in this case, we obtain R
e
1
=
0.25 Ω and R
e
2
<
9.5 Ω, meaning that the earth current divergence impedance R
electrode
<
9.75 Ω. 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 isD (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 Table3. According to the utilization factor expression in Table 2, the radius ratio can be written as α
=
2
R
L
/
R
d
-
1 when n
=
2. The radius ratio α
0 is 0.46 with the electrode interval 0.5 m in configuration A when R
L
=
15.30 Ω and R
d
=
21.00 Ω. For the electrode intervals 1m, 1.5 m, and 2 m, the radius ratios are α
0
/
2, α
0
/
3, and α
0
/
4. The radius ratio α
1 is 0.25 in configuration B when R
L
=
13.11 Ω. It is α
1
/
2 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
Intervalm
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±
5
% 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 is40 m but 80 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.
FrequencykHz
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.9As seen in Table4, 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.


## 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 lengthL, a diameter d, and an embedded depth h, 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
R
=
1
σ
1
2
π
L
ln
⁡
2
L
d
1
+
ln
⁡
(
L
/
2
h
)
ln
⁡
(
2
L
/
d
)
Tagg, Dwight
R
=
1
σ
1
2
π
L
ln
⁡
4
L
d
+
ln
⁡
L
h
-
2
+
2
h
L
-
h
2
L
2
+
h
4
8
L
4
Sunde, Schwarz
R
=
1
σ
1
π
L
ln
⁡
2
L
d
h
-
1
Neutral-point potential method
R
=
1
σ
1
π
L
ln
⁡
2
L
d
h
Average potential method (h
=
0)
R
=
1
σ
1
π
L
ln
⁡
4
L
d
-
1
Neutral-point potential method (h
=
0)
R
=
1
σ
1
π
L
ln
⁡
2
L
dFor 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 impedanceR
d can be expressed as follows:(1)
R
d
=
R
L
n
η
,where R
L is the grounding impedance of a single electrode; n 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 resistanceR
d 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(2)
r
=
L
2
ln
⁡
2
L
/
d
.Forn parallel horizontal grounding electrodes arrayed in-line, the distance between two adjacent electrodes is D. In order to simplify the calculation, a radius ratio α
=
r
/
D is defined. The injection current over each electrode is I
1
,
I
2
,
…
,
I
i
,
…
,
I
n with the the electric potential V. The electric potential of the ith electrode as a reference electrode is given by Φ
i
=
1
/
σ
1
/
2
π
r
I
i, based on the hemispherical grounding equivalent model. The electric potential provided by the i
±
kth electrode at the ith electrode can be expressed as Φ
i
±
k
=
1
/
σ
1
/
2
π
r
α
/
k
I
i
±
k. Then, we obtain the electric potential V
=
∑
i
=
1
n
Φ
i by the superposition principle. By setting each electrode as reference, we have(3)
1
σ
1
2
π
r
1
α
α
2
⋯
α
n
-
1
α
1
α
⋯
α
n
-
2
⋱
α
i
-
1
⋯
1
⋯
α
n
-
i
⋱
α
n
-
1
α
n
-
2
α
n
-
3
⋯
1
I
1
I
2
⋮
I
i
⋮
I
n
=
V
V
⋮
V
⋮
V
R
d
=
V
∑
i
=
1
n
I
i
=
R
L
n
η
α
=
r
D
.Whenn
=
2,3
,
4,5
,
6, the expressions of grounding electrode utilization factor η with argument radius ratio α are obtained and shown in Table 2. The variations of η with α
-
1 are shown in Figure 1.Table 2
Functional relationships betweenη and α.
Numbers of electrodes
Utilization factorη
2
1
1
+
α
3
6
-
7
α
6
+
3
α
-
12
α
2
4
12
-
10
α
12
+
16
α
-
23
α
2
5
360
-
654
α
+
241
α
2
360
+
270
α
-
1495
α
2
+
690
α
3
6
2160
-
2952
α
+
555
α
2
2160
+
3312
α
-
8439
α
2
+
2047
α
3Figure 1
Variation of utilization factorη with α
-
1 under n dipole antennas.From Figure1, it can be seen that the grounding electrode utilization factor η decreases as the parameter α increases for the same number of electrodes n. 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 lengthL and radius b, where b
≪
L, are filled with resistance reducing material with conductivity σ
s. Metal horizontal grounding electrodes with length L and radius a
=
d
/
2 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 b
m. The soil moisture, temperature, and soil gap have a significant influence on the resistance reducing material seepage property, meaning that b
m is variable but always larger than b, as shown in Figure 2.Figure 2
Grounding electrodes surrounded by resistance reducing material with high conductivity.The grounding resistanceR
electrode is equivalent to a series connection circuit of R
c
1, R
e
1, R
c
2, and R
e
2, where R
c
1 is the contact resistance between the resistance reducing material and the electrode, R
e
1 is the resistance reducing material current divergence impedance, R
c
2 is the contact resistance between the resistance reducing material and the earth, and R
e
2 is the earth current divergence impedance:(4)
R
electrode
=
R
c
1
+
R
e
1
+
R
c
2
+
R
e
2
.R
e
1 can be obtained by using(5)
R
e
1
=
1
σ
s
1
π
L
ln
⁡
b
m
a
.It can be simplified by considering the semicylinder for the resistance reducing material and the current divergence impedance.R
e
2 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,(6)
R
e
2
=
1
σ
1
π
L
ln
⁡
L
b
m
.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 impedancesR
c
1 and R
c
2 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 electrodeL is 1.2 m with radius a
=
0.035 m and a grounding hemisphere equivalent radius of r
=
0.17 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 b
=
0.2 m that are filled with resistance reducing material, which gives their equivalent radius b
m
>
b, in this case, we obtain R
e
1
=
0.25 Ω and R
e
2
<
9.5 Ω, meaning that the earth current divergence impedance R
electrode
<
9.75 Ω. 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 isD (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 Table3. According to the utilization factor expression in Table 2, the radius ratio can be written as α
=
2
R
L
/
R
d
-
1 when n
=
2. The radius ratio α
0 is 0.46 with the electrode interval 0.5 m in configuration A when R
L
=
15.30 Ω and R
d
=
21.00 Ω. For the electrode intervals 1m, 1.5 m, and 2 m, the radius ratios are α
0
/
2, α
0
/
3, and α
0
/
4. The radius ratio α
1 is 0.25 in configuration B when R
L
=
13.11 Ω. It is α
1
/
2 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
Intervalm
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±
5
% 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 is40 m but 80 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.
FrequencykHz
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.9As seen in Table4, 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 Table5.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
200The full-wave electromagnetic simulation results with the CST STUDIO SUIT electromagnetic tool are presented with a momentp
=
100 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
f
M kHz
1.5
1.0
0.3
−3 dB bandwidth kHz
5.0
3.0
1.0Figure 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 frequencyf
M at each propagation depth. For the field point at position −200 m, f
M is 1.5 kHz with a −3 dB bandwidth 5 kHz. When the field point is located at position −300 m, f
M
=
1.0 kHz and the −3 dB bandwidth decreases to 3 kHz. For the field point at −600 m, f
M is 0.3kHz and the −3 dB bandwidth is 1 kHz. In this instance, f
M 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 f
M is within the -
3 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 Figure4, 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 26
0
°.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 ofp
=
100 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. Figure5 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 Figure5, 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 6
%. 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


## 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 Table5.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
200The full-wave electromagnetic simulation results with the CST STUDIO SUIT electromagnetic tool are presented with a momentp
=
100 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
f
M kHz
1.5
1.0
0.3
−3 dB bandwidth kHz
5.0
3.0
1.0Figure 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 frequencyf
M at each propagation depth. For the field point at position −200 m, f
M is 1.5 kHz with a −3 dB bandwidth 5 kHz. When the field point is located at position −300 m, f
M
=
1.0 kHz and the −3 dB bandwidth decreases to 3 kHz. For the field point at −600 m, f
M is 0.3kHz and the −3 dB bandwidth is 1 kHz. In this instance, f
M 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 f
M is within the -
3 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 Figure4, 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 26
0
°.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 ofp
=
100 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. Figure5 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 Figure5, 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 6
%. 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 f
M, 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. Thef
M in the full-wave simulation model at position −200 m is 1.5 kHz, 1 kHz at −300 m, and 0.3 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 of300 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.



---

*Source: 8601497-2016-05-09.xml*