# Analysis of an Underground Horizontal Electrically Small Antenna

Authors: Shuwei Dong; Aiguo Yao; Fanhe Meng
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

This work is about a horizontal electrically small wire antenna located underground, which transmits electromagnetic signals to the ground. To solve this problem, the expressions of the magnetic vector potential and the electric field are derived. Further, a quasi-static situation is considered in the condition of extremely low frequency (ELF) or super low frequency (SLF) to make an approximation on Sommerfeld integral for easy calculation. The method of moments (MOM) is used to solve the current distribution along the antenna surface at different frequencies, which lays a good foundation for obtaining the electric field of the antenna. Then the three components of the electric field along the radial distance at different polar angles on the ground are investigated, as well as the voltage received on the ground. Furthermore, some influence factors of the antenna are analyzed in order to know the operating characteristics of the antenna better. The results indicate that the antenna length and relative permittivity have a positive correlation with the radial component of the electric field magnitude, while the buried depth, frequency, and conductivity have a negative correlation with that. The influence of these factors on electric field is obvious except the relative permittivity.


---

## Body

## 1. Introduction

Electromagnetic underground communication plays an important role in the mine production. Through the stratum, the communication system can propagate signals from underground to ground, the kernel technology of which is the use of underground antenna. Unlike conventional communication, the underground communication faces a series of difficulties because the transmission medium is the lossy formation, which will cause attenuation of electromagnetic waves. The issue of electromagnetic propagation through different media has been studied first by Sommerfeld [1] in 1909. Then Sommerfeld [2] established a vertical dipole model placed at the surface of ground and used the evaluation of Fourier-Bessel integrals to solve it. Moore and Blair [3] gave the field expressions of a dipole near a conducting half space in 1961. Bannister [4–6] did a series of studies about quasi-static field of antennas at or above ground using complex image theory and made a conclusion of his achievements in 1979. Chang and Wait [7] did the study of a horizontal wire antenna which is located above or buried in the earth with extremely low frequency and derived explicit expressions for the propagation constant in 1974. King [8] did the study about a dipole over an imperfectly conducting half space as early as 1990. In 1997, King et al. [9] studied a very low frequency antenna in the sea and gave the results of electromagnetic fields. Tai and Collin [10] researched the radiation of a dipole immersed in a dissipative medium in 2000. In 2004, Thottappillil et al. [11] gave the formulas of a semi-infinite antenna above a conducting plane. Vong et al. [12] modeled an electromagnetic telemetry system and gave a portion of results in 2005. Xu and Huang [13] used an efficient computation to solve the problem of an electric dipole above a lossy half space in 2007. Khalatpour et al. [14] used numerical methods to study a wire antenna above a lossy half space in 2010. Parise [15] did a research on quasi-static theory and obtained the harmonic electromagnetic field components excited by a vertical electric dipole lying on the surface of a flat and homogeneous lossy half space in 2013. In addition, there are many scholars who have done researches on related problems. However, the detailed study of a horizontal electrically small antenna in a lossy medium has been not common in previous researches. A horizontal electrically small antenna with ELF or SLF is the best choice for underground electromagnetic communication generator, which can be used effectively in mine communication. This kind of antenna can play a good role in other areas such as energy development, trenchless technology, measurement while drilling, and lateral wells butt.In this paper, an accurate model about a horizontal electrically small antenna located underground with ELF or SLF has been proposed. Analytical solutions of the field excited by underground antenna have been obtained to solve this problem. Sommerfeld integral in the expressions can cause some troubles on calculation, so quasi-static approximation has been used to overcome the difficulties in this low frequency problem. The main work is to get the induction field of the electrically small antenna; for this purpose, numerical methods including the MOM have been used. Numerical results such as the current distribution and the electric field have been obtained in order to describe the antenna character. As a further study, the effect of inside and outside parameters on induction field has been considered and the variation curves of the current and the electric field with the influence factors have been got for better understanding of the antenna.


## 2. Mathematical Procedure

Assume that the earth is homogeneous medium; the magnetic vector potentialA satisfies (1) obtained from the Maxwell’s equations:(1)∇2A+k2A=-μJ,where J is the current density of an impressed source, k is the wave number given by k2=-iωμσ+iωε, ε, μ, and σ are, respectively, the permittivity, permeability, and conductivity of the medium, and ω is the angular frequency. What should be noted is that the time dependence eiωt is assumed and suppressed in this equation and the following analysis.As shown in Figure1, the air occupies the upper half space referred to as medium 0 whereas the ground occupies the lower one referred to as medium 1, and the coordinate origin is at the interface. The horizontal electric dipole is oriented in the x direction and located at the coordinates 0,0,-h, while the coordinates of its image are 0,0,h, where h>0. Considering that the dipole is represented by the current density J=Ilδr-r′ex, the magnetic vector potentials in two media satisfy the following two equations, respectively:(2)∇2A0+k02A0=0,(3)∇2A1+k12A1=-μIlδr-r′ex.Figure 1
Schematic of the horizontal electric dipole.From the boundary conditions, there are two components of the magnetic vector potential oriented inx and z direction, respectively. The field generated by the horizontal dipole is the superposition of the primary field and the reflected field of the source, which can be derived by the image theory. According to what we have known of the magnetic vector potentials in free-space, the x component of the vector potentials in two media satisfies (4) according to the image theory [16]:(4)A0x=μIl4π∫0∞Pkρξ0e-ξ0z+hJ0kρρdkρ,A1x=μIl4π∫0∞kρξ1e-ξ1z+h+Qkρξ1e-ξ1z-hJ0kρρdkρ,where P and Q are undetermined coefficients determined by the boundary conditions and J0 is the zero-order Bessel function of the first kind. Here some variables are given by ρ2=x2+y2, kρ2=kx2+ky2, ξ02=kρ2-k02, and ξ12=kρ2-k12, where k0 and k1 are, respectively, the wave numbers of the air and the ground.The boundary conditions are given by (5), from which the relationship equations between P and Q can be derived:(5)limz→0+⁡A0x=limz→0-⁡A1x,limz→0+⁡∂A0x∂z=limz→0-⁡∂A1x∂z.Further, the coefficientsP and Q are obtained as the following expressions:(6)P=2ξ0ξ1+ξ0eξ0-ξ1h,Q=2ξ1ξ1+ξ0-1.So far, thex component of the vector potentials in two media can be expressed as the following expressions:(7)A0x=μIl4πI0x,(8)A1x=μIl4πe-ik1RR-e-ik1RiRi+I1x,(9)I0x=∫0∞2kρξ1+ξ0e-ξ0z-ξ1hJ0kρρdkρ,(10)I1x=∫0∞2kρξ1+ξ0eξ1z-hJ0kρρdkρ,where R=ρ2+z+h2 and Ri=ρ2+z-h2.The following work is to solve thez component of the vector potentials. Since it is independent of the source, the wave equations (11) can be got:(11)∇2A0z+k02A0z=0,∇2A1z+k12A1z=0.According to the knowledge of partial differential equations, the general solutions of (11) in a cylindrical coordinate system ρ,θ,z are as follows:(12)A0z=μIl4π∑n=-∞∞einθ∫0∞a0ne-ξ0z+b0neξ0zJnkρρdkρ,A1z=μIl4π∑n=-∞∞einθ∫0∞a1ne-ξ1z+b1neξ1zJnkρρdkρ.Because oflimz→∞eξ0z=0 and limz→-∞e-ξ1z=0, in order to ensure that A0z and A1z are finite, the coefficients b0n and a1n must be 0. Meanwhile, the vector potential satisfies A0zρ,θ,z=A0zρ,-θ,z and A1zρ,θ,z=A1zρ,-θ,z, so A0z and A1z are expressed by (13)A0z=μIl4π∑n=-∞∞cos⁡⁡θ∫0∞ane-ξ0zJnkρρdkρ,A1z=μIl4π∑n=-∞∞cos⁡⁡θ∫0∞bneξ1zJnkρρdkρ.The boundary conditions related toA0z and A1z are as follows:(14)limz→0+⁡A0z=limz→0-⁡A1z,(15)limz→0+⁡1k02∂A0x∂x+∂A0z∂z=limz→0-⁡1k12∂A1x∂x+∂A1z∂z.By using the boundary conditions, thez component of the vector potentials in two media is got, respectively:(16)A1z=μIlcos⁡⁡θ4πI1z,A0z=μIlcos⁡⁡θ4πI0z,I0z=∫0∞2k02-k12kρ2ξ1+ξ0k02ξ1+k12ξ0e-ξ1h-ξ0zJ1kρρdkρ,I1z=∫0∞2k02-k12kρ2ξ1+ξ0k02ξ1+k12ξ0eξ1z-hJ1kρρdkρ.The magnetic vector potential underground is obtained so far, and it is(17)A1=μIl4πe-ik1RR-e-ik1RiRi+I1xex+cos⁡⁡θI1zez.In the following analysis, what we are concerned about is the vector potential underground presented asA1. For simplicity, the subscript 1 will be omitted. From Maxwell’s equations, the relationship expression between the electric field E and the vector potential A is shown as(18)E=-iωA+1μσ+iωμε∇∇·A.To solve the electric fieldE, ∇·A should be solved first. Since Aρ=cos⁡θAx and Aθ=-sinθAx, the expression of ∇·A can be derived in cylindrical coordinate system ρ,θ,z:(19)∇·A=1ρ∂ρAρ∂ρ+1ρ∂Aθ∂θ+∂Az∂z=∂Aρ∂ρ+∂Az∂z.By using the recurrence formula(∂/∂ρ)J0kρρ=-kρJ1kρρ, the following result is acquired:(20)∂Az∂z=-μIlcos⁡⁡θ4π∂∂ρI1z′,I1z′=∫0∞2k02-k12kρξ1ξ1+ξ0k02ξ1+k12ξ0eξ1z-hJ0kρρdkρ.Then(21)∇·A=μIl4πcos⁡⁡θ∂∂ρe-ik1RR-e-ik1RiRi+IA,IA=∫0∞2k12kρk02ξ1+k12ξ0eξ1z-hJ0kρρdkρ.Therefore, the three components of the electric field can be obtained in our model:(22)Eρ=cos⁡⁡θIl4πσ1+iωε1k12e-ik1RR-e-ik1RiRi+I1x+∂2∂ρ2e-ik1RR-e-ik1RiRi+IA,Eθ=-sin⁡θIl4πσ1+iωε1k12e-ik1RR-e-ik1RiRi+I1x+1ρ∂∂ρe-ik1RR-e-ik1RiRi+IA,Ez=cos⁡⁡θIl4πσ1+iωε1∂2∂z∂ρe-ik1RR+e-ik1RiRi-k02k12IA.An important problem is how to solve the Sommerfeld integrals in the field expressions. Under the condition of ELF or SLF, the field is in quasi-static situation. In order to facilitate the calculation, the quasi-static approximations can be used. In this case, Wait [17] has shown that k02ξ1+k12ξ0≈k12ξ0, ξ0=kρ2-k02≈kρ, and these approximations have been verified available by many scholars, such as Cooray [18]. In addition, a Wait and Spies [19] image theory approximation is utilized here:(23)ξ1-kρξ1+kρ≈e-kρd,where d=δ1-i and δ is the skin depth given by δ=1/πfμσ1.Thus, by using the quasi-static approximations, it is possible to obtain the approximate expressions of Sommerfeld integralsI1x and IA as follows:(24)I1x≈∫0∞2kρξ1+kρeξ1z-hJ0kρρdkρ≈∫0∞1-e-kρdeξ1z-hJ0kρρdkρ=A-B,A=∫0∞eξ1z-hJ0kρρdkρ,B=∫0∞e-kρdeξ1z-hJ0kρρdkρ,IA≈2∫0∞eξ1z-hJ0kρρdkρ.Use the approximation that Bannister and Dube [5] have verified as eξ1z-h≈eik1mz-hekρnz-h, where the constants m and n get different values in different cases. When ρ2+z-h2≪δ (δ is the skin depth), m=0 and n=1. When ρ2+z-h2 and δ are of the same order of magnitude, m=0.4 and n=0.96 for ρ2+z-h2/δ<1 and m=0.96 and n=0.4 for 1<ρ2+z-h2/δ<10. This approximation has a good match presented in the paper written by Gavriloska et al. [20].Utilizing the Sommerfeld identity [21](25)∫0∞e-ξz-hξJ0kρρkρdkρ=e-ikRR.Then(26)A≈eik1mz-h∫0∞ekρnz-hJ0kρρdkρ=e-ik1m-z+hρ2+n-z+h2,B≈∫0∞e-kρdeik1mz-hekρnz-hJ0kρρdkρ=e-ik1m-z+hρ2+n-z+h+d2,IA≈2eik1mz-h∫0∞ekρnz-hJ0kρρdkρ=2e-ik1m-z+hρ2+n-z+h2.Thus, the three components of the electric field can be defined by more easily calculated expressions with our efforts, which eliminate the difficulties in computing the Sommerfeld integrals. To obtain the electric field, which should be done first, is to calculate the current distribution along the antenna. Take the field point coordinates and the source point coordinates which are, respectively,x,y,z and x′,y′,z′. The horizontal antenna is oriented along x-axis, and the coordinates of its center point are 0,0,z′. Considering that the current is distributed on the antenna surface only, the x component of the electric field generated by current is Ex while the x component of the electric field generated by other sources is Exi. Because the antenna surface is considered to be an ideal conductor, according to the boundary conditions of the ideal conductor, the following equation should be satisfied:(27)Exs+Exis=0.Since the antenna radiusa and the wavelength λ meet the condition a≪λ, the assumption that current is uniformly distributed along the surface is defined. Then the effect of the current along the antenna surface is equivalent to a line current along x-axis.The MOM [22] is used to solve the current distribution, and in this case the MOM formula of the horizontal antenna is as shown in (28). In order to get the numerical result, pulse base and point matching method are chosen in this section:(28)Exi=-Il4πσ1+iωε1∫lIx′·k12e-ik1RR-e-ik1RiRi+I1x+∂2∂x2e-ik1RR-e-ik1RiRi+IAdx′,where R=a2+x-x′2 and Ri=a2+x+x′2.Already knowing the current distribution, the three components of the electric field are easily computed by utilizing formulas (22). The voltage is a parameter that can be directly measured, so we should obtain the voltage received on the ground. Formula (29) shows the expression of the voltage received on the ground by making use of Eρ:(29)Vr=∑iNEρρi+Eρρi+12ρi+1-ρi,where Vr is the voltage received on the ground and ρi is the distance of the ith point.


## 3. Numerical Results and Discussion

### 3.1. Field of the Antenna

Consider an underground horizontal wire antenna which is located 500 m down from the ground surface as the signal generator. The ground is seen as homogenous medium, and the characters areσ=0.075 S/m, ε=20ε0=1.7708×10-10 F/m, and μ=4π×10-7 H/m. The length of the antenna is l=6m while the radius is a=0.05 m. As the feed terminal, there is a 1 V δ voltage source excitation added on the center, where the δ gap is 0.1 m. In the following work, we take the modulus of the complex numbers as the values of the current.The current distribution of the horizontal antenna is calculated by MOM. Figure2 shows the current distribution on the wire antenna surface when the frequency is 10 Hz. The current distribution is approximate triangular distribution, which well matches the known current distribution of electrically small antenna.Figure 2
Current distribution on the wire antenna surface.The electric field on the ground plane is obtained afterwards. The magnitude and the phase of the field are both analyzed to know the field more clearly. Figure3 shows the curves of the three components of the electric field magnitude along the radial distance while Figure 4 shows the curves of the magnitude at different radial distances with the polar angle in polar coordinates.Figure 3
Curves of (a) the radial component of the electric field magnitude, (b) the tangential component of the electric field magnitude, and (c) the axial component of the electric field magnitude with the radial distance.
(a)
(b)
(c)Figure 4
Curves of (a) the radial component of the electric field magnitude, (b) the tangential component of the electric field magnitude, and (c) the axial component of the electric field magnitude with the polar angle.
(a)
(b)
(c)The curves in Figure3(a) represent the variation of the radial component of the electric field along the radial distance in different polar angles. The magnitude of the radial component of the electric field goes through a process that gets smaller to a minimum value from the maximum value first and then gets bigger and afterwards gets smaller gently again along the radial distance in all polar angles. When the polar angle θ=0°, the magnitude is larger than the corresponding ones in the same radial distance when the polar angle takes other values. Curves on Figure 3(b) show the variation of the tangential component of the electric field, from which we can see the magnitude gets smaller with the increasing radial distance, and the magnitude reaches maximum when the polar angel θ=90°. The axial component of the electric field magnitude becomes larger with the increasing distance first and then becomes smaller after a maximum value, which is shown in Figure 3(c). The trends of the axial component of the electric field with the radial distance at different polar angles are the same, but the value at θ=0° is larger than the others. Compared with the other two components, the axial component of the electric field is several orders of magnitude smaller. The curves of the electric field magnitude with the polar angle are shown in Figure 4; they are symmetrical in the polar coordinates and shaped almost like circles in corresponding half areas.The phase of the electric field is also studied as shown in Figure5. The phase of the radial component of the electric field undergoes a mutation at the radial distance between 150 and 200 m, and the other two components get to gradually increase with the increasing radial distance. The phase mutation point is exactly the minimum point of the magnitude, which illustrates the direction of the radial component of the electric field changes in this point. In addition, the trends of the phase variation with the radial distance at different polar angles are similar.Figure 5
Curves of (a) the radial component of the electric field phase, (b) the tangential component of the electric field phase, and (c) the axial component of the electric field phase with the radial distance.
(a)
(b)
(c)Figure6 shows the curves about the variation of the total electric field with the radial distance on the ground. The curve generally presents a tendency which the field magnitude decreases with the increasing distance.Figure 6
Curves of the total electric field magnitude with the radial distance.The curve of the voltage received on the ground with the radial distance is shown in Figure7. Here the original point on the ground is taken as the relative voltage zero point, and the voltage value is gained at the polar angle θ=0°. The received voltage value tends to increase first and gets to decrease after a maximum value. The maximum value appears at the point that the direction of the radial component of the electric field begins to change.Figure 7
Curve of the voltage received on the ground with the radial distance.


### 3.2. Influence of the Buried Depth on Field

The buried depth of the antenna can change the electric field distribution on the ground. In this part, the maximum value of the current on the antenna surface is chosen to study the influence of the buried depth on current distribution. In Figure8, it is shown that the current value tends to increase with the deepening depth, but the influence is not obvious numerically. The radial component of the electric field is of our most concern, which directly affects the voltage received, so the radial component of the electric field can be a representative to study the influence on electric field. Three points corresponding to the radial distance R=50m, R=150m, and R=250m at the polar angle θ=0° are taken as the study objects. As shown in Figure 9, the radial component of the electric field magnitude becomes smaller when the buried depth becomes deeper at different radial distances, which indicates that the deeper the buried depth, the smaller the electric field magnitude in the underground communication system.Figure 8
Curve of the maximum current value with the buried depth.Figure 9
Curves of the radial component of the electric field magnitude with the buried depth.


### 3.3. Influence of the Antenna Length on Field

The antenna size is a factor that can influence the electric field. Since the model is a line antenna model, the radius is much smaller than the length, so the study aims at the influence of the antenna length on the electric field. The same study objects as the study on Section3.2 are used in this case. Figure 10 shows the curve of the maximum current value with the antenna length from 5 m to 7 m, which demonstrates that the current value has an obvious increase when the antenna length gets longer. Figure 11 represents the curves of the radial component of the electric field magnitude with the antenna length at different radial distance. The data indicates that the field magnitude gets larger at any radial distance with the increasing antenna length. So we can obtain the conclusion that longer antenna will produce better results.Figure 10
Curve of the maximum current value with the antenna length.Figure 11
Curves of the radial component of the electric field magnitude with the antenna length.


### 3.4. Influence of the Frequency on Field

In lossy medium, the attenuation of the electric field increases as the frequency gets higher, so the frequency is an important parameter that can influence the electric field on the ground. Similarly, the maximum current value on the antenna surface and the radial component of the electric field magnitude can be the representative to study the influence of the frequency on electric field. In Figure12, it is shown that the current value tends to decrease with the increasing frequency, but the influence is too tiny to distinguish the size of the value. As shown in Figure 13, the field value at the points whose radial distances R=50m and R=250m decreases when the frequency increases, but the field value at the point R=150m gets a slight increase when the frequency increases. The point R=150m is just in the vicinity of the point where the radial component of the electric field magnitude takes the minimum value and the field phase undergoes a mutation. This shows that the radial distance of the mutation point increases with the increasing frequency.Figure 12
Curve of the maximum current value with the frequency.Figure 13
Curves of the radial component of the electric field magnitude with the frequency.


### 3.5. Influence of the Stratum Parameters on Field

The characteristic of the formation affects the field as well, two of which are the stratum conductivity and the relative permittivity. Similar with the previous study, we choose the maximum current value on antenna surface as the object to study the influence on current and use the radial component of the electric field magnitude corresponding to the pointsR=50m, R=150m, and R=250m at the polar angle θ=0° to study the influence on electric field. Figure 14 shows the curve of current value varying with the conductivity, which indicates that the current value has a positive correlation with the conductivity. Curves in Figure 15 demonstrate the variation of the radial component of the electric field magnitude with the conductivity from 0.05 to 0.1 S/m at the three points selected. With the increasing conductivity, the field magnitudes at the three points all get smaller. There are obvious variations at the points R=50m and R=250m and unconspicuous variation at the point R=150m.Figure 14
Curve of the maximum current value with the conductivity.Figure 15
Curves of the radial component of the electric field magnitude with the conductivity.The effect of the relative permittivity on current and electric field is studied later. Curve in Figure16 represents the variation of the current value on antenna surface with the relative permittivity. The current value increases when the relative permittivity gets larger, but the influence is so tiny that the change is not obvious. Figure 17 represents the dependence of the radial component of the electric field with the relative permittivity corresponding to the point R=50m at the angle θ=0°. The value of the radial component of the electric field is positively correlated with the relative permittivity, and the variation is too tiny to recognize the difference of the value as well.Figure 16
Curve of the maximum current value with the relative permittivity.Figure 17
Curve of the radial component of the electric field magnitude with the relative permittivity.


### 3.6. Study on the Asymmetric Horizontal Antenna

In the practical application, such as trenchless technology, the antenna used in the underground communication system may be asymmetric. So the asymmetric horizontal antenna is an important subject needed to study. Take one antenna arm length as 3 m, 4.16 m, 5.32 m, and 6.48 m, respectively, and the other antenna arm length as fixed 3 m, and the current distribution on the antenna surface is calculated. The curves in Figure18 describe the variation of the current value on the antenna surface with different antennas. The current distribution is nearly triangular distribution no matter whether the antenna is symmetrical, and the distribution is asymmetric on asymmetric antenna. The current value gets bigger when the antenna length gets longer, and the change is relatively obvious.Figure 18
Current distribution on the asymmetric antenna surface.


### 3.7. Verification

In order to verify the numerical results, a simulation work has been done by using the commercial software HFSS. The model simulated is the same as the previous study with 10 Hz frequency. The ground characters areσ=0.075 S/m, ε=20ε0=1.7708×10-10 F/m, and μ=4π×10-7 H/m. The current distribution on antenna surface is presented in Figure 19, which matches well with the previous result in Figure 2. Curve on Figure 20 shows the variation of the total field along the radial distance. The electric field value in this case is a little smaller than the result in Figure 6. This consequence may be due to the fact that more factors are concerned in HFSS and the model in this work is too large to calculate. Overall, this HFSS case is largely consistent with the previous results as verification.Figure 19
Current distribution on the antenna surface by HFSS.Figure 20
Curve of the total electric field with the radial distance by HFSS.


## 3.1. Field of the Antenna

Consider an underground horizontal wire antenna which is located 500 m down from the ground surface as the signal generator. The ground is seen as homogenous medium, and the characters areσ=0.075 S/m, ε=20ε0=1.7708×10-10 F/m, and μ=4π×10-7 H/m. The length of the antenna is l=6m while the radius is a=0.05 m. As the feed terminal, there is a 1 V δ voltage source excitation added on the center, where the δ gap is 0.1 m. In the following work, we take the modulus of the complex numbers as the values of the current.The current distribution of the horizontal antenna is calculated by MOM. Figure2 shows the current distribution on the wire antenna surface when the frequency is 10 Hz. The current distribution is approximate triangular distribution, which well matches the known current distribution of electrically small antenna.Figure 2
Current distribution on the wire antenna surface.The electric field on the ground plane is obtained afterwards. The magnitude and the phase of the field are both analyzed to know the field more clearly. Figure3 shows the curves of the three components of the electric field magnitude along the radial distance while Figure 4 shows the curves of the magnitude at different radial distances with the polar angle in polar coordinates.Figure 3
Curves of (a) the radial component of the electric field magnitude, (b) the tangential component of the electric field magnitude, and (c) the axial component of the electric field magnitude with the radial distance.
(a)
(b)
(c)Figure 4
Curves of (a) the radial component of the electric field magnitude, (b) the tangential component of the electric field magnitude, and (c) the axial component of the electric field magnitude with the polar angle.
(a)
(b)
(c)The curves in Figure3(a) represent the variation of the radial component of the electric field along the radial distance in different polar angles. The magnitude of the radial component of the electric field goes through a process that gets smaller to a minimum value from the maximum value first and then gets bigger and afterwards gets smaller gently again along the radial distance in all polar angles. When the polar angle θ=0°, the magnitude is larger than the corresponding ones in the same radial distance when the polar angle takes other values. Curves on Figure 3(b) show the variation of the tangential component of the electric field, from which we can see the magnitude gets smaller with the increasing radial distance, and the magnitude reaches maximum when the polar angel θ=90°. The axial component of the electric field magnitude becomes larger with the increasing distance first and then becomes smaller after a maximum value, which is shown in Figure 3(c). The trends of the axial component of the electric field with the radial distance at different polar angles are the same, but the value at θ=0° is larger than the others. Compared with the other two components, the axial component of the electric field is several orders of magnitude smaller. The curves of the electric field magnitude with the polar angle are shown in Figure 4; they are symmetrical in the polar coordinates and shaped almost like circles in corresponding half areas.The phase of the electric field is also studied as shown in Figure5. The phase of the radial component of the electric field undergoes a mutation at the radial distance between 150 and 200 m, and the other two components get to gradually increase with the increasing radial distance. The phase mutation point is exactly the minimum point of the magnitude, which illustrates the direction of the radial component of the electric field changes in this point. In addition, the trends of the phase variation with the radial distance at different polar angles are similar.Figure 5
Curves of (a) the radial component of the electric field phase, (b) the tangential component of the electric field phase, and (c) the axial component of the electric field phase with the radial distance.
(a)
(b)
(c)Figure6 shows the curves about the variation of the total electric field with the radial distance on the ground. The curve generally presents a tendency which the field magnitude decreases with the increasing distance.Figure 6
Curves of the total electric field magnitude with the radial distance.The curve of the voltage received on the ground with the radial distance is shown in Figure7. Here the original point on the ground is taken as the relative voltage zero point, and the voltage value is gained at the polar angle θ=0°. The received voltage value tends to increase first and gets to decrease after a maximum value. The maximum value appears at the point that the direction of the radial component of the electric field begins to change.Figure 7
Curve of the voltage received on the ground with the radial distance.


## 3.2. Influence of the Buried Depth on Field

The buried depth of the antenna can change the electric field distribution on the ground. In this part, the maximum value of the current on the antenna surface is chosen to study the influence of the buried depth on current distribution. In Figure8, it is shown that the current value tends to increase with the deepening depth, but the influence is not obvious numerically. The radial component of the electric field is of our most concern, which directly affects the voltage received, so the radial component of the electric field can be a representative to study the influence on electric field. Three points corresponding to the radial distance R=50m, R=150m, and R=250m at the polar angle θ=0° are taken as the study objects. As shown in Figure 9, the radial component of the electric field magnitude becomes smaller when the buried depth becomes deeper at different radial distances, which indicates that the deeper the buried depth, the smaller the electric field magnitude in the underground communication system.Figure 8
Curve of the maximum current value with the buried depth.Figure 9
Curves of the radial component of the electric field magnitude with the buried depth.


## 3.3. Influence of the Antenna Length on Field

The antenna size is a factor that can influence the electric field. Since the model is a line antenna model, the radius is much smaller than the length, so the study aims at the influence of the antenna length on the electric field. The same study objects as the study on Section3.2 are used in this case. Figure 10 shows the curve of the maximum current value with the antenna length from 5 m to 7 m, which demonstrates that the current value has an obvious increase when the antenna length gets longer. Figure 11 represents the curves of the radial component of the electric field magnitude with the antenna length at different radial distance. The data indicates that the field magnitude gets larger at any radial distance with the increasing antenna length. So we can obtain the conclusion that longer antenna will produce better results.Figure 10
Curve of the maximum current value with the antenna length.Figure 11
Curves of the radial component of the electric field magnitude with the antenna length.


## 3.4. Influence of the Frequency on Field

In lossy medium, the attenuation of the electric field increases as the frequency gets higher, so the frequency is an important parameter that can influence the electric field on the ground. Similarly, the maximum current value on the antenna surface and the radial component of the electric field magnitude can be the representative to study the influence of the frequency on electric field. In Figure12, it is shown that the current value tends to decrease with the increasing frequency, but the influence is too tiny to distinguish the size of the value. As shown in Figure 13, the field value at the points whose radial distances R=50m and R=250m decreases when the frequency increases, but the field value at the point R=150m gets a slight increase when the frequency increases. The point R=150m is just in the vicinity of the point where the radial component of the electric field magnitude takes the minimum value and the field phase undergoes a mutation. This shows that the radial distance of the mutation point increases with the increasing frequency.Figure 12
Curve of the maximum current value with the frequency.Figure 13
Curves of the radial component of the electric field magnitude with the frequency.


## 3.5. Influence of the Stratum Parameters on Field

The characteristic of the formation affects the field as well, two of which are the stratum conductivity and the relative permittivity. Similar with the previous study, we choose the maximum current value on antenna surface as the object to study the influence on current and use the radial component of the electric field magnitude corresponding to the pointsR=50m, R=150m, and R=250m at the polar angle θ=0° to study the influence on electric field. Figure 14 shows the curve of current value varying with the conductivity, which indicates that the current value has a positive correlation with the conductivity. Curves in Figure 15 demonstrate the variation of the radial component of the electric field magnitude with the conductivity from 0.05 to 0.1 S/m at the three points selected. With the increasing conductivity, the field magnitudes at the three points all get smaller. There are obvious variations at the points R=50m and R=250m and unconspicuous variation at the point R=150m.Figure 14
Curve of the maximum current value with the conductivity.Figure 15
Curves of the radial component of the electric field magnitude with the conductivity.The effect of the relative permittivity on current and electric field is studied later. Curve in Figure16 represents the variation of the current value on antenna surface with the relative permittivity. The current value increases when the relative permittivity gets larger, but the influence is so tiny that the change is not obvious. Figure 17 represents the dependence of the radial component of the electric field with the relative permittivity corresponding to the point R=50m at the angle θ=0°. The value of the radial component of the electric field is positively correlated with the relative permittivity, and the variation is too tiny to recognize the difference of the value as well.Figure 16
Curve of the maximum current value with the relative permittivity.Figure 17
Curve of the radial component of the electric field magnitude with the relative permittivity.


## 3.6. Study on the Asymmetric Horizontal Antenna

In the practical application, such as trenchless technology, the antenna used in the underground communication system may be asymmetric. So the asymmetric horizontal antenna is an important subject needed to study. Take one antenna arm length as 3 m, 4.16 m, 5.32 m, and 6.48 m, respectively, and the other antenna arm length as fixed 3 m, and the current distribution on the antenna surface is calculated. The curves in Figure18 describe the variation of the current value on the antenna surface with different antennas. The current distribution is nearly triangular distribution no matter whether the antenna is symmetrical, and the distribution is asymmetric on asymmetric antenna. The current value gets bigger when the antenna length gets longer, and the change is relatively obvious.Figure 18
Current distribution on the asymmetric antenna surface.


## 3.7. Verification

In order to verify the numerical results, a simulation work has been done by using the commercial software HFSS. The model simulated is the same as the previous study with 10 Hz frequency. The ground characters areσ=0.075 S/m, ε=20ε0=1.7708×10-10 F/m, and μ=4π×10-7 H/m. The current distribution on antenna surface is presented in Figure 19, which matches well with the previous result in Figure 2. Curve on Figure 20 shows the variation of the total field along the radial distance. The electric field value in this case is a little smaller than the result in Figure 6. This consequence may be due to the fact that more factors are concerned in HFSS and the model in this work is too large to calculate. Overall, this HFSS case is largely consistent with the previous results as verification.Figure 19
Current distribution on the antenna surface by HFSS.Figure 20
Curve of the total electric field with the radial distance by HFSS.


## 4. Conclusion

In this work, a more accurate model of horizontal electrically small antenna is established to solve the problem of through the earth communication. The expressions of the electric field generated by an underground horizontal antenna are derived in forms which can be easily solved in this mathematical model. Based on this work, numerical method is used to obtain the current distribution and the electric field on the ground in the condition of 10 Hz. In these numerical results, the variations of three components of the electric field with the radial distance at different polar angles are investigated, respectively. The radial and axial components of the electric field get the maximum value at the polar angleθ=0°, and the tangential component of the electric field gets the maximum value at the polar angle θ=90°. With the increasing radial distance, the radial component of the electric field at the polar angle θ=0° decreases first and achieves the minimum value at the phase mutation point R=150m and then makes a change of getting larger and afterwards smaller again. The tangential component of the electric field at the polar angle θ=90° diminishes with the radial distance, and the axial component of the electric field at the polar angle θ=0° appears as a trend that increases first and decreases afterwards with the increasing radial distance. In addition, the voltage received on the ground is obtained. At the polar angle θ=0°, the maximum value of the voltage is approximately at the radial distance R=150 m.In order to have a relatively overall understanding of the antenna, the influence of the buried depth, antenna length, frequency, and stratum parameters on field is studied. The consequence shows that the buried depth, antenna length, conductivity, and relative permittivity are positively related to the current value while the frequency is negatively related to the current. The influence of the antenna length and conductivity on the current value is obvious. The study of these influence factors on the electric field shows that the antenna length and relative permittivity have a positive correlation with the radial component of the electric field magnitude, and the buried depth, frequency, and conductivity have a negative correlation with that. The influence of the relative permittivity is very slight. Additionally, the current distribution on the asymmetric horizontal antenna surface is studied.



---

*Source: 531729-2015-05-03.xml*