Analysis of an Underground Horizontal Electrically Small Antenna

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 fieldmagnitude, 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.


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][5][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 quasistatic 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.

Mathematical Procedure
Assume that the earth is homogeneous medium; the magnetic vector potential A satisfies (1) obtained from the Maxwell's equations: where J is the current density of an impressed source,  is the wave number given by  2 = −( + ), , , 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   is assumed and suppressed in this equation and the following analysis.As shown in Figure 1, 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  direction and located at the coordinates (0, 0, −ℎ), while the coordinates of its image are (0, 0, ℎ), where ℎ > 0.
Considering that the dipole is represented by the current density J = (r − r  )e  , the magnetic vector potentials in two media satisfy the following two equations, respectively: From the boundary conditions, there are two components of the magnetic vector potential oriented in  and  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  component of the vector potentials in two media satisfies (4) according to the image theory [16]: where  and  are undetermined coefficients determined by the boundary conditions and  0 is the zero-order Bessel function of the first kind.Here some variables are given by and , where  0 and  1 are, respectively, the wave numbers of the air and the ground.
The boundary conditions are given by ( 5), from which the relationship equations between  and  can be derived: Further, the coefficients  and  are obtained as the following expressions: International Journal of Antennas and Propagation 3 So far, the  component of the vector potentials in two media can be expressed as the following expressions: where  = √ 2 + ( + ℎ) 2 and The following work is to solve the  component of the vector potentials.Since it is independent of the source, the wave equations ( 11) can be got: According to the knowledge of partial differential equations, the general solutions of (11) in a cylindrical coordinate system (, , ) are as follows: Because of lim  → ∞   0  = 0 and lim  → −∞  − 1  = 0, in order to ensure that  0 and  1 are finite, the coefficients  0 and  1 must be 0.Meanwhile, the vector potential satisfies  0 (, , ) =  0 (, −, ) and  1 (, , ) =  1 (, −, ), so  0 and  1 are expressed by The boundary conditions related to  0 and  1 are as follows: lim By using the boundary conditions, the  component of the vector potentials in two media is got, respectively: The magnetic vector potential underground is obtained so far, and it is In the following analysis, what we are concerned about is the vector potential underground presented as A 1 .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 To solve the electric field E, ∇ ⋅ A should be solved first.Since   = cos   and   = − sin   , the expression of ∇ ⋅ A can be derived in cylindrical coordinate system (, , ): By using the recurrence formula (/) 0 (  ) = −   1 (  ), the following result is acquired: Then International Journal of Antennas and Propagation Therefore, the three components of the electric field can be obtained in our model: 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
Utilizing the Sommerfeld identity [21] Then 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, (, , ) and (  ,   ,   ).The horizontal antenna is oriented along -axis, and the coordinates of its center point are (0, 0,   ).Considering that the current is distributed on the antenna surface only, the  component of the electric field generated by current is E  while the  component of the electric field generated by other sources is E   .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: Since the antenna radius  and the wavelength  meet the condition  ≪ , 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 -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: where 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   : where   is the voltage received on the ground and   is the distance of the th point.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. Figure 2 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.

Numerical Results and Discussion
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.Figure 3 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.
The curves in Figure 3(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 Figure 5.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 6 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.
The curve of the voltage received on the ground with the radial distance is shown in Figure 7.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.

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 Figure 8, 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  = 50 m,  = 150 m, and  = 250 m 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.

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 Section 3.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.

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 International Journal of Antennas and Propagation electric field.In Figure 12, 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 This shows that the radial distance of the mutation point increases with the increasing frequency.

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 points  = 50 m,  = 150 m, and  = 250 m 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     The effect of the relative permittivity on current and electric field is studied later.Curve in Figure 16 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  = 50 m 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.
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 Figure 18 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.

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.

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  = 150 m 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  = 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

Figure 1 :
Figure 1: Schematic of the horizontal electric dipole.

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  = 6 m while the radius is  = 0.05 m.

Figure 2 :
Figure 2: Current distribution on the wire antenna surface.

Figure 3 :
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.

Figure 4 :
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.

Figure 5 :Figure 6 :
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.

Figure 7 :
Figure 7: Curve of the voltage received on the ground with the radial distance.

Figure 8 :Figure 9 :
Figure 8: Curve of the maximum current value with the buried depth.

Figure 10 :
Figure 10: Curve of the maximum current value with the antenna length.

Figure 11 :
Figure 11: Curves of the radial component of the electric field magnitude with the antenna length.

Figure 12 :
Figure 12: Curve of the maximum current value with the frequency.

Figure 13 :Figure 14 :Figure 15 :
Figure 13: Curves of the radial component of the electric field magnitude with the frequency.

Figure 16 :
Figure 16: Curve of the maximum current value with the relative permittivity.

Figure 17 :Figure 18 :
Figure 17: Curve of the radial component of the electric field magnitude with the relative permittivity.

Figure 19 :Figure 20 :
Figure 19: Current distribution on the antenna surface by HFSS.