Next Article in Journal
Comparison of Damage Indexes for Assessing Seismic Fragility of Bearings in an Offshore Bridge
Previous Article in Journal
Assessment of RANS Turbulence Models in Prediction of the Hydrothermal Plume in the Longqi Hydrothermal Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electromagnetic Field Variation of ELF Near-Region Excited by HED in a Homogeneous Half-Space Model

1
School of Electrical Engineering, Naval University of Engineering, Wuhan 430033, China
2
Institute of Education, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(13), 7499; https://doi.org/10.3390/app13137499
Submission received: 15 May 2023 / Revised: 21 June 2023 / Accepted: 23 June 2023 / Published: 25 June 2023

Abstract

:
Great attention has been paid to the propagation of electromagnetic (EM) waves across the sea surface due to its important applications. Most of the previous research, however, focuses on the half-space model illustrating the deep sea environment. In this paper, EM field distribution in the extremely low frequency (ELF) near-region under horizontal electric dipole (HED) excitation in homogeneous half-space seawater is analyzed based on the general expression of the Sommerfeld integral using the quasistatic approximation method. The focus is on deriving complete and effective solutions in air and seawater regions under the cylindrical coordinates for the EM near-field, which is generated by an HED in a shallow sea. The resulting formulas can be given by a few summands in closed form as the well-known Fourier–Bessel integrals. The analytical approximate expression of ELF Sommerfeld EM field integral excited by the HED in the homogeneous half-space seawater is deduced under the condition that the propagation distance ρ satisfies << 1. To this end, the EM field distribution in the range close to the HED antenna in seawater is simulated, the results have shown that the minimum attenuation value of the vertical electric component Ez is about 15 dB, and that of the radical magnetic components Hφ is about 30 dB, and these values are found to be of greatest potential for the near-field region propagation among the electric and magnetic components. Finally, the correctness of the proposed method is verified by comparison with Pan’s approximation method and Margetis’s exact expression approximation method, which demonstrated the correctness of the proposed method.

1. Introduction

Due to the high salinity and complex temperature and current distribution characteristics of seawater, its conductivity and dielectric constant are higher than those of air. The higher the conductivity of water, the greater the attenuation of EM waves. Therefore, the propagation characteristics of EM waves in seawater are considerably different from those in air. Experiments have shown that the amplitude and phase characteristics of ELF (3–30 Hz) EM waves are stable when they propagate in seawater, with small propagation loss and strong seawater penetration ability. Therefore, ELF EM waves have broad application prospects in geological exploration, submarine exploration, seismic detection, submarine communication, and other forms of geophysical exploration [1,2,3,4,5,6]. Many scholars have conducted in-depth studies on EM fields generated by low-frequency HED located in seawater or soil. In the representative monograph written by A. Banos [7], a complete analytical solution of the EM field excited by a dipole in half-space is given. In the decades that followed, R. W. P. King [8] described in detail the lateral EM waves generated by the excitation of horizontal and vertical electric dipoles (VEDs) in two-layer and three-layer media. It is worth mentioning that the lateral wave term of his solution can be expressed as an approximate expression containing a Fresnel integral. Several years later, Margetis derived the exact solution expression of an EM field excited by electric dipoles at the interface of half-space [9]. However, the transmitting point and receiving point can only be on the ground; off the ground is not accurate. Figure 1 demonstrates the EM wave propagation model in the half-space model in the two-layer medium in King’s monograph, in which L is a lateral wave, L1 is a direct wave, and L2 is a reflected wave.
Generally speaking, the long-term focus of EM research has been on the far-field characteristics of traditional problems such as EM propagation, radiation, and scattering, while the attention and research on the near-field characteristics are still relatively weak [10,11,12]. Sommerfeld was the first to give the approximate integral expression of EM field for electric dipoles in homogeneous half-space (hereinafter referred to as the Sommerfeld integral) [13], but this expression is difficult to solve by simple numerical integration. The key reason is that the Fourier–Bessel expression in the Sommerfeld integral has non-integrable terms at → 0.
To solve this problem, D. Margetis [14] gave the exact near-field expression of the fundamental oscillator in half-space under the condition of >> 1 by assuming that k2/k1 (k2 and k1 represent the wave number in seawater and air, respectively) tends to infinity. However, the half-space model of “air–seawater” cannot satisfy the condition of |k2/k1| << 1. Therefore, the exact formula is not suitable for solving the ELF EM field of the “air–seawater” half-space model. The quasistatic approximate method was first proposed by Wait [15], but it uses Laplace and Poisson equations to calculate the field in the air, so it cannot be applied to the near-field in half-space under the condition of k2 << k1. In recent years, this problem was also revisited by W. Y. Pan and K. Li [16], who addressed the ELF near-field wave propagation for the application of the marine controlled-source electromagnetics method as an extension of the above work. The computational scheme is developed from the approximated formulas in the book [8] by retaining fundamental terms. This motivates us to employ Maclaurin’s expansion near the poles of integrands from the Fourier–Bessel representations. H. L. Xu [17] proposed an approximate quasistatic approximation method for solving the “sea–rock” half-space, but the ELF near-field distribution in the deep-sea half-space was obtained using a simplified Sommerfeld integral in the paper, which was not applicable to the “air–seawater” half-space environment.
In addition, with the rapid development of computer technology and computational EM methods, efficient numerical methods have been developed for solving far-field problems with arbitrary source distribution. However, the calculation of the low-frequency near-field has always been a difficult problem to solve. For example, the method of moment (MOM) [18,19,20] has a singularity in the calculation of the near-field, and it is difficult to deal with it. The finite element method (FEM) [21,22,23] and finite difference time domain (FDTD) [24,25,26] have errors caused by boundary truncation; the former may have pseudo-solutions, and the latter is troubled by the difficulty of accurately fitting the surface boundary. In addition, numerical methods may even have problems such as “low-frequency crash” when the frequency is in the very low ranges. Therefore, in view of the above problems, we still have to use the regression analytic method to calculate the EM field propagation in the near-field region of ELF.
For ELF, the wavelength can be compared with the circumference of the Earth. For example, when the working frequency of ELF is 10 Hz, the wavelength is close to three-quarters of the circumference of the Earth. Thus, when the field and source points are in close proximity, it can be regarded as a near-field propagation problem. The propagation environment can be assumed as an idealized model, i.e., a homogeneous half-space, and the interface as a homogeneous plane. In this study, the EM field distribution in the ELF near-field region excited by an HED in homogeneous half-space seawater is studied based on the Sommerfeld integral. The expression is then simplified using an approximate method under the condition that the propagation distance k0ρ << 1 (k0 is the wave number in air and ρ is the propagation distance) in the near-field region. Then, the analytical approximate expression and spatial distribution of the EM field components in the near-field region are derived via simulation. Finally, we compare the quasistatic approximate calculation results with Pan’s approximation method and Margetis’s exact expression approximation method to verify the effectiveness of the proposed method.

2. Propagation Model of ELF Near-Field Region Excited by HED

The EM field model in the near-field region of ELF excited by an HED in an unbounded, homogeneous half-space medium is shown in Figure 2. The homogeneous half-space medium model can be an air–seawater or air–lithosphere model. This paper focuses on the air–seawater medium. The sea surface was taken as the z = 0 plane; the sea region and air region were marked as regions 1 and 2, respectively.
Assuming that the HED antenna is located in the positive direction of the z-axis, parallel to the x-axis direction, the distance from the xoy plane, i.e., the sea level, is set as d meters; the electric moment is Idl; and the distance from the observation point P to the sea level is z meters. In a cylindrical coordinate system, the coordinates of the observation point P can be expressed as P(ρ, φ, z). The components of each electromagnetic field at point P are shown in Figure 2. Region 1 (z ≥ 0) was filled with seawater. µ0, σ1, and εr1 are the permeability, conductivity, and relative dielectric constant, respectively, and the wave number can be expressed as k 1 = ε r 1 ω / c , where ω is the angular frequency and c is the speed of light. Region 2 (z ≤ 0) is an air medium with a permeability, conductivity, relative permittivity, and wave number of μ0, σ2, εr2 = 1, and k2 = k0 = ω/c, respectively. ρ is the propagation distance between the projection point of the electric dipole source and the observation point P, and r1 is the geometric distance between the electric dipole source point and the observation point P in seawater. For low-frequency electric dipole excitation, which varies with the time harmonic factor e−jωt, the integral expression of the EM field generated in region 1 (z ≥ 0) is given in A. N. Sommerfeld’s work [13] and will not be repeated here. The approximate analytical method is used to simplify the Sommerfeld integral of the ELF near-field region EM fields in homogeneous half-space seawater and derive its approximate analytical expression for the EM fields.

3. Analytical Approximate Expression of the Field Component

3.1. Quasistatic Approximate Propagation Conditions

Under the condition that the propagation distance satisfies k2ρ << 1, the Fourier–Bessel function integral in the expression of the ELF electromagnetic field at or near the boundary of the two layers of media has discrete terms. The distribution of divergent poles is shown in Figure 3.
In order to evaluate the half-space Sommerfeld integral using the complex function method, the integrand must be divided into single-valued branches; that is, several slits must be made in the complex plane of λ, and the multivalued functions must be specified clearly and reasonably. For this reason, we take the radical function γ1 = (k12λ2)1/2 as an example and first analyze its amplitude change on the complex plane; if k1 = β1 + 1, λ = ξ + , then
γ 1 = ( β 1 2 α 1 2 ) ξ 2 η 2 + 2 j ( α 1 β 1 ξ η )
On the isoaxial hyperbola ξη = α1β1, the amplitude angle of γ1 in Equation (1) is divided into three parts, as shown in Figure 3c: the amplitude angle of γ1 near the virtual axis is 0, the amplitude angle of γ1 on the dotted line is π/2, and the amplitude angle of γ1 away from the virtual axis is π. If two symmetric slits are made over the entire λ complex plane, then γ1 in the λ complex plane is singly analytic.
Under the same condition, these integrals diverge in the traditional sense. For example, the multivalued function D(λ) = k12γ2 + k22γ1 included in the Sommerfeld integral has four zeros on the Riemann surface. γ1 is the propagation parameter in the seawater medium of region 1, and γ2 is the propagation parameter in the air medium of region 2. Therefore, two sets of symmetrical slits are created on the λ complex plane, as shown in Figure 3a. Here, slits parallel to the imaginary axis were used to simplify the calculation. According to the loop integral theory of the complex variable function [27], the integral path along the real axis of the Sommerfeld integral can be converted into an integral along the slits C1 and C2.

3.1.1. Approximation of the Propagation Parameter γ2

Air and seawater have different electric parameters in the “air–seawater” homogeneous half-space model; therefore, the ELF near-field EM fields of the HED near the interface of the two media conforms to the quasistatic assumption: ω → 0, k2ρ << 1, k2 << k1. Based on this assumption, the simplified root function can be approximated as follows:
γ 2 = k 2 2 λ 2 = lim ω 0 ( ε r 2 1 / 2 c ω ) 2 λ 2 j λ
where ω is the angular frequency, λ is the wavelength, and j is an imaginary unit.

3.1.2. Approximation When |k2| << |k1|

Considering that the wave number k2 of the air medium is considerably lower than that of the seawater medium k1, i.e., |k2| << |k1|, the corresponding |k22γ1|/|k12γ2| approaches zero and the higher-order term can be ignored:
k 2 2 γ 1 k 1 2 γ 2 k 2 2 k 1 2 k 1 2 λ 2 j λ k 2 k 1 k 2 λ 2 1 < < 1
According to the aforementioned assumptions, the integration path in Figure 3b can replace the original integration path in Figure 3a; the wave number k2 is very close to the origin of the coordinate system. Figure 3c shows the geometric diagram of γ1.

3.1.3. Approximation d << ρ and z << ρ

Referring to (1), under the quasistatic assumption, the slot position of the Sommerfeld integral fulcrum k2 is close to the zero point of the imaginary axis. Because the ELF near-field conditions in seawater meet the relationship equations d << ρ and z << ρ, the multivalue function of the Sommerfeld integral in seawater can be approximated as follows:
k 1 2 γ 2 k 2 2 γ 1 k 1 2 γ 2 + k 2 2 γ 1 = 1 k 2 2 γ 1 / k 1 1 γ 2 1 + k 2 2 γ 1 / k 1 1 γ 2 1
Substituting the Equations (2)~(4) into the Sommerfeld integral expressions, the Sommerfeld integral expression of each electromagnetic field component at any observation point P in the homogeneous seawater medium can be simplified to
E 1 ρ = ω μ 0 I d l 4 π k 1 2 cos ϕ k 1 2 F 1 ( k 1 , d 0 ) 1 2 F 2 ( k 1 , d 0 ) + 1 2 G 4 ( k 1 , d 1 ) k 1 2 2 H 1 ( d 1 )
E 1 ϕ = ω μ 0 I d l 4 π k 1 2 sin ϕ k 1 2 F 1 ( k 1 , d 0 ) 1 2 F 2 ( k 1 , d 0 ) + 1 2 G 5 ( k 1 , d 1 ) k 1 2 2 H 2 ( d 1 )
E 1 z = j ω μ 0 I d l 4 π k 1 2 cos ϕ ± F 3 ( k 1 , d 0 ) + G 1 ( k 1 , d 1 )
B 1 ρ = μ 0 I d l 4 π sin ϕ ± F 4 ( k 1 , d 0 ) + 1 2 G 2 ( k 1 , d 1 ) 1 2 H 3 ( d 1 )
B 1 ϕ = μ 0 I d l 4 π cos ϕ ± F 4 ( k 1 , d 0 ) + 1 2 G 3 ( k 1 , d 1 ) 1 2 H 4 ( d 1 )
B 1 z = j μ 0 I d l 4 π sin ϕ F 5 ( k 1 , d 0 ) H 5 ( d 1 )
where F1 (k1, d0)~F5 (k1, d0) represents the direct wave from the field source point (0, 0, d) to any near-field observation point P(ρ, φ, z) in the lower half-space, G1 (k1, d1)~G5 (k1, d1) represents the reflected wave from the field source image point (0, 0, −d) to any near-field observation point (ρ, φ, z) in the lower half-space, and H1 (d1)~H5 (d1) represents the side wave affected by the interface of the two media. φ is azimuth coordinates. r0 is the geometric distance from the source point to the receiving point of the seawater medium, r1 is the geometric distance from the image point of the field source to the receiving point of the air medium, d0 is the vertical height between the source and receiving points of the electric dipole in the seawater medium, and d1 is the vertical height between the image point of the field source and receiving point of the air medium. r0, r1 and d0, d1 are defined as follows:
d 0 = z d ,   d 1 = z + d
r 0 = d 0 2 + ρ 2 ,   r 1 = d 1 2 + ρ 2

3.2. Bessel Integral Solution

The integral expressions of each field component from Equation (5)~(10) can be reduced to infinite integrals containing Bessel functions. The expressions of direct and reflected waves can be obtained by strictly solving the corresponding integral. According to the general solution of the integral formula given by King, the analytical expressions of the integrals F1 (k1, d0)~F5 (k1, d0), G1 (k1, d1)~G5 (k1, d1) can be obtained [8], and we have given the detailed solution in the Appendix A. Therefore, we only need to find the solution of the integral formulas H1 (d1)~H5 (d1) to obtain the final expression of the electromagnetic field.
As the integral formulas H1 (d1)~H5 (d1) of the side wave component contain a slow decay oscillation function, which cannot be directly obtained analytically, we use the approximate method for asymptotic estimation. H1 (d1)~H5 (d1) are expressed as follows [13]:
H 1 ( d 1 ) = 0 γ 2 γ 1 γ 2 + γ 1 [ J 0 ( λ ρ ) + J 2 ( λ ρ ) ] γ 1 1 e j γ 1 d 1 λ d λ
H 2 ( d 1 ) = 0 γ 2 γ 1 γ 2 + γ 1 [ J 0 ( λ ρ ) J 2 ( λ ρ ) ] γ 1 1 e j γ 1 d 1 λ d λ
H 3 ( d 1 ) = 0 γ 2 γ 1 γ 2 + γ 1 [ J 0 ( λ ρ ) J 2 ( λ ρ ) ] e j γ 1 z + d 1 λ d λ
H 4 ( d 1 ) = 0 γ 2 γ 1 γ 2 + γ 1 [ J 0 ( λ ρ ) + J 2 ( λ ρ ) ] e j γ 1 z + d 1 λ d λ
H 5 ( d 1 ) = 0 γ 2 γ 1 γ 2 + γ 1 J 1 ( λ ρ ) e j γ 1 d 1 γ 1 1 λ 2 d λ
where Jn (λρ), n = 0, 1, 2, is the expression of the n-order Bessel function.
In Equations (13)~(17), the integrand function contains the multivalued function term A(λ) = (γ2 − γ1)/(γ2 + γ1). As it is known from the form text that the intermediate variable γ1 can be expressed by γ1jλ, the absolute values for γ1, γ2, and k1 are approximately defined by a triangular relationship, as shown in Figure 3d. According to the quasistatic approximation (2),
0 < γ 2 γ 1 j λ k 1 2 λ 2 < λ 2 + k 1 2 λ 2 = k 1
Based on Equation (18), the absolute value of vector γ2 − γ1 is proportional to the absolute value of vector k1, which can be deduced as follows:
γ 2 γ 1 k 1 d 1 j λ k 1 2 λ 2 k 1 d 1 ξ k 1 d 1 < < 1
where |ξ|≤1. Therefore, for infinitely small imaginary numbers [17],
j γ 2 γ 1 k 1 d 1 e j ( γ 2 γ 1 k 1 ) d 1 1
In addition, the McLaughlin expansion of the function e j ( γ 2 γ 1 k 1 ) d m is reserved for the first two orders:
e j ( γ 2 γ 1 k 1 ) d 1 1 + j ( γ 2 γ 1 k 1 ) d 1 ( ( γ 2 γ 1 k 1 ) d 1 ) 2 2 ! , < 0 <
Therefore,
( γ 2 γ 1 ) 2 = 2 d 1 2 2 e j ( γ 2 γ 1 k 1 ) d 1 d 1 2 + 2 j ( γ 2 γ 1 k 1 ) d 1 d 1 2 k 1 2 + 2 k 1 ( γ 2 γ 1 )
Substituting Equations (20) into (22), and simplifying Equation (22) further, we obtain the following:
( γ 2 γ 1 ) 2 2 j k 1 e j ( γ 2 γ 1 k 1 ) d 1 1 d 1 + k 1 2
Combining the multivalued function A(λ) = (γ2 − γ1)/(γ2 + γ1) and Equation (23), we obtain the following:
A ( λ ) = γ 2 γ 1 γ 2 + γ 1 = ( γ 2 γ 1 ) 2 k 2 2 k 1 2 2 j k 1 e j ( γ 2 γ 1 k 1 ) d 1 1 ( k 2 2 k 1 2 ) d 1 + k 1 2 k 2 2 k 1 2
According to the property of the Bessel function, H1 (d1)~H5 (d1) along the real axis can be converted into the integral of slits C1 and C2 in Figure 3b, and the approximate expressions of H1 (d1)~H5 (d1) can be obtained as follows:
H 1 ( d 1 ) = 0 γ 2 γ 1 2 k 2 2 k 1 2 J 0 λ ρ + J 2 λ ρ γ 1 1 e j γ 1 d 1 λ d λ 2 k 2 2 k 1 2 d 1 ρ 2 2 j + k 1 d 1 e j k 1 d 1 e j k 1 r 1 + 2 j e j k 2 d 1 e j k 1 r 1
H 2 ( d 1 ) = 0 γ 2 γ 1 2 k 2 2 k 1 2 J 0 λ ρ J 2 λ ρ γ 1 1 e j γ 1 d 1 λ d λ 2 k 2 2 k 1 2 d 1 2 j k 1 k 2 e j k 2 d 1 ρ 2 + j k 2 r 1 1 ρ 2 e j k 2 r 1 2 j k 1 d 1 e j k 2 d 1 ρ 2 + j k 2 r 1 1 ρ 2 e j k 1 r 1
H 3 ( d 1 ) = 0 γ 2 γ 1 2 k 2 2 k 1 2 J 0 λ ρ J 2 λ ρ e j γ 1 z + d 1 λ d λ 2 k 1 k 2 2 k 1 2 d 1 2 j k 1 d 1 e j k 1 d 1 ρ 2 + d 1 e j k 1 r 1 r 1 j k 1 r 1 1 r 1 2 1 ρ 2 2 j e j k 2 d 1 ρ 2 + d 1 e j k 2 r 1 r 1 j k 2 r 1 1 r 1 2 1 ρ 2
H 4 ( d 1 ) = 0 γ 2 γ 1 2 k 2 2 k 1 2 J 0 λ ρ + J 2 λ ρ e j γ 1 z + d 1 λ d λ 2 k 1 k 2 2 k 1 2 d 1 2 j + k 1 d 1 e j k 1 d 1 ρ 2 d 1 e j k 1 r 1 r 1 ρ 2 + 2 j e j k 2 d 1 ρ 2 d 1 e j k 2 r 1 r 1 ρ 2
H 5 ( d 1 ) = 0 γ 2 γ 1 2 k 2 2 k 1 2 J 1 λ ρ γ 1 1 e j γ 1 d 1 λ 2 d λ k 1 ρ k 2 2 k 1 2 d 1 2 k 2 e j k 2 r 1 j r 1 2 1 k 2 r 1 3 e j k 1 r 1 j k 1 2 d 1 + 2 k 1 j r 1 2 1 k 1 r 1 3

3.3. Simplification of Field Component Expression

Let the electric moment be Idl = 1 A‧m; substituting Equations (25) to (29) into (5) to (10), respectively, we obtain the approximate expression of each EM field component. Taking the EM field components Ez as an example, this paper gives the process of simplifying the near-field approximate expression of the vertical magnetic field component. For k2ρ << 1, the expressions of Ez can be given as follows:
E 1 z = j ω μ 0 ρ 4 π k 1 cos ϕ d 0 r 0 k 1 r 0 2 + 3 j r 0 3 3 k 1 r 0 4 e j k 1 r 0 d 1 r 1 k 1 r 1 2 + 3 j r 1 3 3 k 1 r 1 4 e j k 1 r 1
In the equation number, “+” is taken when zd, and “−” is taken when 0 ≤ zd. According to the quasistatic assumption, the distance relationship between the field source point and observation point is ρ >> d, ρ >> z, and r0r1ρ. Therefore, the approximate expression of the vertical EM field components Ez satisfying the quasistatic assumption can be obtained as follows:
E 1 z = j ω μ 0 ρ 4 π k 1 cos ϕ k 1 ρ 2 + 3 j ρ 3 3 k 1 ρ 4 e j k 1 ρ d 0 d 1
The second approximation of e j k 1 ρ was expanded using McLaughlin’s progressive formula. In the process of simplification, the attenuation rate of EM waves in air was considerably less than that in seawater; thus, it meets the requirements of Im(k2) << Im(k1). The quasistatic approximate expression can be further simplified as follows:
E 1 z = j ω μ 0 4 π k 1 cos ϕ k 1 2 ρ 2 + j k 1 2 2 ρ d 0 + d 1
In the same way, the near-field approximate expressions of the EM field components in the remaining “air–seawater” half-space can be obtained.

4. Spatial Distribution of EM Field Intensity in the Near-Field Region

4.1. Spatial Distribution of EM Field Intensity

To further study the ELF near-field distribution of the HED antenna in seawater, we selected the model shown in Figure 2 and used the algorithm to evaluate the distribution of the EM field intensity of the observation plane on the square grid 100 m × 100 m at the plane of z = 1 m.
Assuming that the distance from the horizontal antenna to the sea surface is d = 10 m, the electric moment is Idl = 1A‧m, φ = π/4, the working frequency is 10 Hz, the propagation distance is ρ = 5000 m, and the height from the observation point to the sea surface is z = 1 m; μ0 and ε0 are the permeability and dielectric constant in a vacuum, respectively. The conductivity and relative permittivity of seawater are σ1 = 4 S/m and ε1 = 80, respectively. The EM field components changed with the height of the observation point. In the model shown in Figure 2, the coordinate origin is taken as the center, and a square area with a length and width of 100 m is obtained at the plane of z = 1 m. The amplitude distributions of each field component were simulated using Delphi software and are shown in Figure 4.
In the model shown in Figure 2, an HED antenna with a working frequency of 10 Hz is placed underwater at 1 m away from the sea surface. A voltage source is fed at the center point of the antenna. It is considered that the HED antenna is excited by an equal amplitude and phase source in order to analyze the law of field distribution in the underwater HED antenna. The radius of the thin wire is 2 cm. The wavelength of an EM wave at 10 Hz in seawater is approximately equal to 0.5 m for seawater conductivity equal to 4 S/m, and it is about 3 × 107 m in air.
The results for electric fields are shown in Figure 4a–c, which include four color intensity plots for the three electric field components and the magnitude of the electric field. Eρ exhibits maximum field intensity along the x-axis, and minimum field intensity along the y-axis, while Ez clearly shows that the shape and polarized direction for the HED antenna are placed along the y-axis; Eφ demonstrates that there is a strong electric field intensity in the y-axis direction and the fastest attenuation in the x-axis direction. In addition, the vertical electric component Ez is much stronger than Eρ and Eφ for near-field region distance propagation in any direction, and the minimum attenuation is 15 dB. At the same time, we should also focus on the electric field strength in the x-axis direction of Eρ and the y-axis direction of Eφ. Similar plots for magnetic fields are shown in Figure 4d–f. Figure 4d,e show the field distribution in air; the field distributions of Hρ and Hφ at 10 Hz are drawn on the z = 1 m plane above the air–sea interface. It can be seen that Hρ and Hφ have an obvious attenuation mutation point at 40 m, which may be caused by changes in the interface between air and seawater. Meanwhile, we can see that Hφ receives the strongest magnetic field strength in the x-axis direction, the decay rate is the slowest, and the lowest attenuation value is 30 dB, so we should pay attention to the magnetic field strength in this direction. Figure 4f shows that the maximum field intensity of Hz is along the x-axis, and the minimum field intensity is along the y-axis, while the radical magnetic component Hz is a little stronger than Hρ and Hφ for near-field region distance in any direction, and the attenuation speed is relatively stable, and the minimum attenuation rate is maintained at 40 dB, which is suitable for near-field region distance propagation.

4.2. Comparative Validation Analysis

In order to verify the effectiveness of the method presented in this paper, Pan’s approximation [27] and Margetis’s exact expression [9] are respectively obtained by MacClaurin’s progressive simplification, and the obtained results are compared with the method presented in this paper. If both the field and source points are in seawater, the seawater and air conductivity are σ1 = 4 S/m and σ2 = 1.8 × 10−4 S/m, respectively, and the relative dielectric constants are εr1 = 80 and εr2 = 1, respectively. When calculated at 3 Hz and 30 Hz, the electric field component in the ρ direction changes with the distance, as shown in Figure 5.
As can be seen from Figure 5, the electric field components in the ρ direction obtained by different methods decrease with the increase in radial distance, and the variation trend is basically consistent. When the frequency is the same, the amplitude of E1ρ obtained by Pan approximation decreases the most with distance, while the amplitude of E1ρ obtained by Margetis’s exact expression is in the middle of the two approximation methods, and the amplitude obtained by this method decreases at the slowest rate, which also verifies that the method proposed in this paper has a good accuracy in calculating the electric field component in the near-field region. It can be seen from Figure 5 that at the same site, the E1ρ amplitude decreases with the increase in operating frequency, and the amplitude decreases greatly in the ELF band high-frequency band. When the horizontal distance is less than 5000 m, the E1ρ amplitude decreases significantly with the increase in the horizontal distance, and when the horizontal distance is greater than 5000 m, the decline trend tends to be gentle. As shown in Figure 5a, when the electric field component E1ρ is operating at 3 Hz, the field strength decreases with the increase in the horizontal distance in the direction of ρ. If the horizontal distance is greater than 10 km, the seawater has little influence on the E1ρ amplitude of the field point when it is higher than 500 m.

5. Conclusions

In this research, we employed the “air–seawater” homogeneous half-space model to obtain EM field integral expression in ELF near-field region under HED antenna excitation based on the general expression of the Sommerfeld integral using the quasistatic approximation method. The EM field distribution of an underwater dipole was simulated to calculate the field distribution in shallow seawater. After the calculation, the EM field distribution law in the vicinity of the HED antenna in seawater was analyzed. The simulated results show different attenuations with distance for each component. The vertical electric component Ez and radical magnetic components Hφ are found to be of the greatest potential for long-distance propagation among the electric and magnetic components, respectively. Finally, we verified the correctness of the proposed method by comparing it with Pan’s approximation method and Margetis’s exact expression approximation method, exhibiting the great potential of the proposed method for various kinds of applications in shallow sea environments. The proposed method could be also applied to the excitation of vertical and horizontal magnetic dipoles (VMDs/HMDs) and vertical electric dipoles (VEDs).

Author Contributions

Conceptualization, T.D.; methodology, H.X.; software, S.H.; validation H.X. and S.H.; formal analysis, H.X. and S.H.; investigation, H.X., S.H., and T.D.; resources, S.H. and T.D.; data curation, T.D.; writing—original draft preparation, S.H.; writing—review and editing, H.X. and S.H.; visualization, S.H.; supervision, H.X.; project administration, H.X.; funding acquisition, H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author on reasonable request.

Acknowledgments

We would like to thank Jianchen Feng for his support of the article and for providing some help during the preparation of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The strict solution of F1 (k1, d0)~F5 (k1, d0) and G1 (k1, d1)~G5 (k1, d1) in the Sommerfeld integral is as follows:
F 1 ( k 1 , d 0 ) = 0 e j γ 1 d 0 γ 1 J 0 ( λ ρ ) λ d λ = j e j k 1 r 0 r 0
F 2 ( k 1 , d 0 ) = 0 e j γ 1 d 0 γ 1 J 0 ( λ ρ ) J 2 ( λ ρ ) λ 3 d λ = 2 k 1 j k 1 r 0 2 r 0 2 2 j k 1 r 0 3 d 0 2 r 0 2 ( j k 1 r 0 3 r 0 2 3 j k 1 r 0 3 ) e j k 1 r 0
F 3 ( k 1 , d 0 ) = 0 e j γ 1 d 0 J 1 ( λ ρ ) λ 2 d λ = ρ k 1 d 0 r 0 k 1 r 0 2 + 3 i r 0 3 3 k 1 r 0 4 e j k 1 r 0
F 4 ( k 1 , d 0 ) = 0 e j γ 1 d 0 J 0 ( λ ρ ) λ d λ = d 0 r 0 ( j k 1 r 0 2 1 r 0 2 ) e j k 1 r 0
F 5 ( k 1 , d 0 ) = 0 e j γ 1 d 0 γ 1 J 1 ( λ ρ ) λ 2 d λ = ρ k 1 ( 1 r 0 2 + 1 k r 0 3 ) e j k 1 r 0
G 1 ( k 1 , d 1 ) = 0 e j γ 1 d 1 J 1 ( λ ρ ) λ 2 d λ = ρ k 1 d 1 r 1 k 1 r 1 2 + 3 j r 1 3 3 k 1 r 1 4 e j k 1 r 1
G 2 ( k 1 , d 1 ) = 0 e j γ 1 d 1 J 0 ( λ ρ ) + J 2 ( λ ρ ) λ d λ = 2 e j k 1 d 1 ρ 2 d 1 r 1 e j k 1 r 1 ρ 2
G 3 ( k 1 , d 1 ) = 0 e j γ 1 d 1 J 0 ( λ ρ ) J 2 ( λ ρ ) λ d λ = 2 ( 1 ρ 2 + d 1 r 1 j k 1 r 1 1 r 1 2 1 ρ 2 ) e j k 1 r 1
G 4 ( k 1 , d 1 ) = 0 γ 1 e j γ 1 d 1 [ J 0 ( λ ρ ) J 2 ( λ ρ ) ] λ d λ = 2 k 1 e j k 1 d 1 ρ 2 + 1 r 1 2 + j k 1 r 1 3 + d 1 2 r 1 2 j k 1 r 1 3 r 1 2 3 j k 1 r 1 3 1 ρ 2 e j k 1 r 1
G 5 ( k 1 , d 1 ) = 0 γ 1 e j γ 1 d 1 [ J 0 ( λ ρ ) + J 2 ( λ ρ ) ] λ d λ = 2 k 1 e j k 1 d 1 ρ 2 + j e j k 1 r 1 k 1 r 1 3 1 + j k 1 r 1 d 1 2 ρ 2
where Jn (λρ), n = 0, 1, 2, is the expression of n-order Bessel function, λ is the wavelength, and ρ is propagation distance. γ1 is the propagation parameter in the seawater medium of region 1, and γ2 is the propagation parameter in the air medium of region 2. r0 is the geometric distance from the source point to the receiving point of the seawater medium, r1 is the geometric distance from the image point of the field source to the receiving point of the air medium, d0 is the vertical height between the source and receiving points of the electric dipole in the seawater medium, and d1 is the vertical height between the image point of the field source and receiving point of the air medium. k1 is the wave number in the seawater medium of region 1. j is an imaginary unit.

References

  1. Ozaki, M.; Yagitani, S.; Nagano, I.; Miyamura, K. Ionospheric penetration characteristics of ELF waves radiated from a current source in the lithosphere related to seismic activity. Radio Sci. 2009, 44, 1–12. [Google Scholar] [CrossRef]
  2. Zhao, G.Z.; Wang, L.F.; Ji, T.; Chen, X.B.; Zhan, Y.; Xiao, Q.B.; Wang, J.J.; Cai, J.T.; Xu, G.J.; Wan, Z.S.; et al. New experiments of CSELF electromagnetic method for earthquake monitoring. Chin. J. Geophys. 2010, 53, 479–486. [Google Scholar]
  3. Liu, Y.; Kong, L.; Ding, K.; Zhuo, X.J. Wireless electromagnetic method and its application in discovery of marine rource. Shipbuiding China 2019, 60, 21–27. [Google Scholar]
  4. Zhang, P. Application of seismic ELF electromagnetic observation in Yunnan, China. Recent Dev. World Seismol. 2019, 8, 158–159. [Google Scholar]
  5. Sainson, S. Electromagnetic Seabed Logging; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  6. Yusof, A.B.M.; Kabir, S. An overview of sonar and electromagnetic waves for underwater communication. IETE Tech. Rev. 2012, 29, 307–317. [Google Scholar] [CrossRef]
  7. Banos, A. Dipole Radiation in the Presence of a Conducting Half-Space; Pergamon: Oxford, UK, 1966. [Google Scholar]
  8. King, R.W.P.; Owens, M.; Wu, T.T. Lateral Electromagnetic Waves; Springer: New York, NY, USA, 1992; pp. 151–170. [Google Scholar]
  9. Margetis, D.; Wu, T.T. Exactly calculable field components of electric dipoles in planar boundary. J. Math. Phys. 2001, 42, 713–745. [Google Scholar] [CrossRef]
  10. He, T.; Zhang, X.W.; Pan, W.Y.; Li, K. Near-field of a VLF electric dipole in an anisotropic plasma. IEEE Trans. Antennas Propag. 2019, 67, 4040–4048. [Google Scholar] [CrossRef]
  11. Zhima, Z.; Huang, J.; Shen, X.; Xia, Z.; Chen, L.; Piersanti, M.; Yang, Y.; Wang, Q.; Zeng, L.; Lei, J.; et al. Simultaneous observations of ELF/VLF Rising-Tone quasiperiodic waves and energetic electron precipitations in the High-Latitude upper ionosphere. JGR Space Phys. 2020, 125, 1–21. [Google Scholar] [CrossRef]
  12. Ghosal, S.; De, A.; Chakrabarty, A.; Shubair, R.M. Vertical dipole above the lossy dielectric half-space—A Characteristic Mode Analysis. IEEE Trans. Electromagn. Compat. 2020, 62, 2832–2841. [Google Scholar] [CrossRef]
  13. Sommerfeld, A.N. Propagation of waves in wireless telegraphy. Ann. Phys. 1909, 26, 665–736. [Google Scholar] [CrossRef] [Green Version]
  14. Margetis, D. Radiation of horizontal electric dipole on large dielectric sphere. J. Math. Phys. 2002, 43, 3162–3201. [Google Scholar] [CrossRef] [Green Version]
  15. Wait, J.R. Electromagnetic Waves in Stratified Media, 2nd ed.; Pergamon: New York, NY, USA, 1970. [Google Scholar]
  16. Pan, W.Y.; Li, K. Propagation of SLF/ELF Electromagnetic Waves, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 161–219. [Google Scholar]
  17. Xu, H.L.; Gu, T.T.; Li, K. Approximated solutions for ELF Near-Field propagation due to a horizontal electric dipole excitation near the Sea-rock boundary. IEEE Trans. Antennas Propag. 2018, 66, 2471–2480. [Google Scholar] [CrossRef]
  18. Tasic, M.S.; Kolundzija, B.M. Method of Moment weighted domain decomposition method for scattering from large platforms. IEEE Trans. Antennas Propag. 2018, 66, 3577–3589. [Google Scholar] [CrossRef]
  19. Xu, L.Y.; Yang, H.; Shen, R.J. Scattering analysis of multiobject electromagnetic systems using stepwise Method of Moment. IEEE Trans. Antennas Propag. 2018, 11, 1–8. [Google Scholar] [CrossRef]
  20. Xu, X.; Xie, H.; Luo, Z.; Wan, L. Analysis of the electrical performance of multituned VLF thirteen-tower umbrella antenna. Int. J. Antennas Propag. 2022, 2022, 1–10. [Google Scholar] [CrossRef]
  21. Cai, H.Z.; Xiong, B.; Zhdanov, M. Three-dimensional marine controlled-source electromagnetic modelling in anisotropic medium using finite element method. Chin. J. Geophys. 2015, 58, 557–563. [Google Scholar]
  22. Cao, X.Y.; Yin, C.C.; Zhang, B. A goal-oriented adaptive finite-element method for 3D MT anisotropic modeling with topography. Chin. J. Geophys. 2018, 61, 2618–2628. [Google Scholar]
  23. Gao, J.; Wang, J.; Yi, Z. Simulation of detecting small underwater insulative target through electromagnetic field based on finite element method. Spaceraft Environ. Eng. 2022, 39, 293–297. [Google Scholar]
  24. Marshall, R.A.; Wallace, T.; Turbe, M. Finite-Difference Modeling of Very-Low-Frequency Propagation in the Earth-Ionosphere Waveguide. IEEE Trans. Antennas Propag. 2017, 65, 7185–7197. [Google Scholar] [CrossRef]
  25. Berenger, J.-P. FDTD propagation of VLF-LF waves in the presence of ions in the earth-ionosphere waveduide. Ann. Telecommun 2020, 75, 437–446. [Google Scholar] [CrossRef]
  26. Shahmohamadi, B.; Sarraf Shirazi, R.; Moradi, G. Analysis of dipole plasma antenna using kinetic method and FDTD numerical approach. J. Electron. Commun. 2022, 145, 154006. [Google Scholar] [CrossRef]
  27. Pan, W.Y. LF VLF ELF Wave Propagation; The University of Electronic Science and Technology Press: Chengdu, China, 2004; pp. 46–48. [Google Scholar]
Figure 1. Half-space model of EM wave propagation in two-layer media.
Figure 1. Half-space model of EM wave propagation in two-layer media.
Applsci 13 07499 g001
Figure 2. Model diagram of the ELF near-field region excited by HED.
Figure 2. Model diagram of the ELF near-field region excited by HED.
Applsci 13 07499 g002
Figure 3. The geometric shape of γ1 (related to the Sommerfeld integral) due to the excitation of the horizontal electric dipole: (a) corresponding pole and integral path of the Sommerfeld integral; (b) integral path when k2 << k1; (c) geometric diagram of γ1; (d) triangular quantitative relationship.
Figure 3. The geometric shape of γ1 (related to the Sommerfeld integral) due to the excitation of the horizontal electric dipole: (a) corresponding pole and integral path of the Sommerfeld integral; (b) integral path when k2 << k1; (c) geometric diagram of γ1; (d) triangular quantitative relationship.
Applsci 13 07499 g003
Figure 4. Distribution diagram of each field component of the electric and magnetic field intensity at the plane of z = 1 m. (a) Eρ (dB); (b) Eφ (dB); (c) Ez (dB); (d) Hρ (dB); (e) Hφ (dB); (f) Hz (dB).
Figure 4. Distribution diagram of each field component of the electric and magnetic field intensity at the plane of z = 1 m. (a) Eρ (dB); (b) Eφ (dB); (c) Ez (dB); (d) Hρ (dB); (e) Hφ (dB); (f) Hz (dB).
Applsci 13 07499 g004
Figure 5. Magnitudes of E1z versus the directional distance in near field (φ = π/2). (a) Magnitudes of E1z for 3 Hz; (b) magnitudes of E1z for 30 Hz.
Figure 5. Magnitudes of E1z versus the directional distance in near field (φ = π/2). (a) Magnitudes of E1z for 3 Hz; (b) magnitudes of E1z for 30 Hz.
Applsci 13 07499 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, S.; Xie, H.; Ding, T. Electromagnetic Field Variation of ELF Near-Region Excited by HED in a Homogeneous Half-Space Model. Appl. Sci. 2023, 13, 7499. https://doi.org/10.3390/app13137499

AMA Style

Hu S, Xie H, Ding T. Electromagnetic Field Variation of ELF Near-Region Excited by HED in a Homogeneous Half-Space Model. Applied Sciences. 2023; 13(13):7499. https://doi.org/10.3390/app13137499

Chicago/Turabian Style

Hu, Sumou, Hui Xie, and Tao Ding. 2023. "Electromagnetic Field Variation of ELF Near-Region Excited by HED in a Homogeneous Half-Space Model" Applied Sciences 13, no. 13: 7499. https://doi.org/10.3390/app13137499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop