1. Introduction
The background of this research project is that after disasters such as underground space collapse, conventionally installed underground communication systems are damaged, making it difficult for ground search-and-rescue personnel to contact trapped personnel, resulting in low search and rescue efficiency. Highly reliable communication technology equipment is crucial in the emergency rescue process.
Ground-penetrating communication technology is a wireless communication technology that uses the Earth’s strata as a communication channel. Therefore, this technology is not affected by disasters such as underground engineering collapses and can achieve two-way communication, representing an important communication method for emergency rescue in the case of underground engineering disasters such as those in tunnels and mines. However, in such situations, the achievement of bidirectional signal transmission is a key bottleneck problem that this technology faces [
1,
2].
At present, there are three main types of ground-penetrating communication technologies studied domestically and internationally: magnetic induction near-field electromagnetic wave ground-penetrating communication, ground electrode current field ground-penetrating communication, and mechanical vibration wave ground-penetrating communication. Near-field electromagnetic wave and elastic wave transmission through the ground based on magnetic induction antennae have been widely discussed by researchers at home and abroad; however, owing to the need for larger circular antennae and greater power, the transmission attenuation is considerable, making the technology susceptible to the influence of geological media, resulting in poor applicability. Ground electrode current field communication has the advantages of small antenna size, a simple structure, and flexibility, enabling the achievement of stable communication and long-distance transmission [
3,
4].
In 2012, Spanish scholar V. Bataller proposed a current field through ground communication based on the injection of ground electrode current and the detection of potential differences between receiving electrodes, and proposed a circuit-equivalent model for electrode impedance. In 2015, L. Van and C. Sunderman comprehensively considered the noise level and signal attenuation characteristics in soil, and derived an analytical solution for the field distribution of an electrode-based Through the Earth(TTE)communication system in a uniform half space. In 2016, researchers from Beijing Jiaotong University studied the propagation and attenuation characteristics of electromagnetic waves in layered media channels and established a layered media model and calculation model. In 2017, a research team from Beijing University of Information Technology designed an anti-multipath interference information system based on a direct-sequence spread spectrum. By using highly random logistic spread spectrum codes, the anti-multipath interference performance of the tunnel engineering geological exploration information system was effectively improved. In 2017, researchers from China University of Mining and Technology designed a wireless underground magnetic induction information transmission system to achieve short-distance wireless data transmission in underground environments. In 2017, Ralchenko et al. from Carlton University in Canada discovered the presence of low-frequency signals on two slender conductors near a transmitter, which were found to considerably increase the signal transmission distance [
5,
6,
7,
8,
9,
10,
11].
In 2021, Dongkuo University in Seoul, South Korea developed a practical multi hop magnetic induction wireless underground sensor network deployment model and evaluated the transmission and reception performance of the network. These studies provide us with important references and references, and have important research significance for promoting the development and application of ground penetrating communication technology. In 2021, Liu Baoheng and others from Naval Engineering University constructed a three-dimensional model of an underground electric field, providing complete explicit solutions for three electric field components. Based on the model and expressions, the effects of the offset, relative orientation, length, and operating frequency of the transmitting and receiving antennas on the electric field distribution were simulated and analyzed. In 2023, Zhao Jun and others from Shanxi Mechanical and Electrical Vocational and Technical College used the finite element method to solve the normal potential and abnormal potential of the space field. Based on Comsol software, a homogeneous three-dimensional geoelectric detection model was constructed, and tetrahedral mesh adaptive algorithm was used for dissection. The numerical solution accuracy was compared and analyzed. Study the distribution and variation of the focusing current field, and determine the range of the focusing effect current ratio coefficient [
12,
13,
14].
The principle of ground penetrating communication based on the current field of grounding electrodes is to apply low-frequency electrical signals to two electrodes hit on a rock or soil layer, forming a current field in the rock or soil layer. The signal voltage is detected between the two receiving electrodes at a distance, thereby achieving long-distance transmission of signals. This article studies the principle of this technology and the signal design, signal detection and recognition of ground electrode arrays, providing theoretical guidance and experimental basis for further development in this field.
2. Technical Principles
Earth electrode current injection refers to the injection of a certain frequency of current into the ground through one or more earth electrodes below the surface, resulting in a charged volume and the formation of a current field within that volume. The current field can cause the movement of underground charges, resulting in the propagation of electromagnetic waves underground.
2.1. Mechanism of Earth Electrode Current Injection
(1) The Basic Steps of the Mechanism of Earth Electrode Current Injection
The process of injecting a ground electrode current can be divided into the following basic steps:
- a)
Inject current into the ground through a flexibly grounded electrode;
- b)
Current moves in the direction of the current through underground charges;
- c)
Electric current forms a charged volume underground;
- d)
A current field is formed inside the volume, causing the movement of underground charges;
- e)
The movement of charges forms electromagnetic waves, thereby achieving ground-penetrating communication.
(2) Factors Influencing the Mechanism of Earth Electrode Current Injection
The implementation of the earth electrode current injection mechanism is influenced by multiple factors, among which the most important factors include the following:
- a)
Electrode shape and material: During the process of current injection into the ground electrode, the shape of the electrode and the quality of the material determine the effect of current injection into the ground and the intensity of the generated current field.
- b)
Frequency and intensity of current: The frequency and intensity of the current have significant impacts on the movement of charges injected into the ground. The frequency and intensity of the current should not be too high, as excessive current can cause serious electromagnetic interference.
- c)
Geological conditions and terrain: Geological conditions and terrain have a significant impact on the propagation rate and amplitude of electromagnetic waves. Some special underground rock formations or terrain conditions can produce phenomena such as reflection and refraction, which affect the propagation rate and amplitude of electromagnetic waves.
(3) Application of the Earth Electrode Current Injection Mechanism
Ground electrode current injection is widely used in the fields of earth-penetrating communication, geophysics, underground mineral exploration, etc. Ground electrode current injection is one of the key technical links to achieve ground-penetrating communication, playing an important role in the propagation and reception of underground electromagnetic waves.
2.2. The Principle of Ground Electrode Current Field through Ground Communication Technology
The ground electrode current field penetrate ground communication technology is a communication method that utilizes the ground current field to penetrate the earth medium to transmit information. The principle of this method is that the transmitting device injects a high-frequency current into a pair of electrodes (i.e., ground electrodes) and forms a ground current field signal through the underground ground layer medium. The receiving device detects, identifies, processes, and demodulates the signal by receiving the ground current field signal, thereby achieving information transmission.
The transmitting device injects a high-frequency current into the ground through a ground electrode, and the current field strength decreases with depth. At the same time, reflection, refraction, and other phenomena occur at different depths, resulting in different intensities and phases of the injected current. These currents propagate underground, forming a current field.
During the transmission process, the conductivity of underground rock layers can have an impact on the transmission effect, and the shape, size, depth, and contact with the ground of the ground electrode can also lead to changes in the current field.
The frequency used to inject the ground electrode current also affects the transmission effect. Generally speaking, the delay of high-frequency signals underground is shorter than that of low-frequency signals. Therefore, selecting an appropriate frequency can improve the transmission efficiency of the ground current field through ground communication technology.
3. Signal Design
Due to the special nature of the geological environment, the terrestrial channel is very complex, and the attenuation of high-frequency electromagnetic waves is very severe. In through-the-ground communication, only lower-frequency waves can be selected for transmission, thereby limiting the bandwidth of the channel. Therefore, the optimal working mode of the excitation source is studied to adapt the signal transmission mode, working frequency, and the condition of the collapsed medium. Through waveform design, energy is concentrated in the main conduction direction, which is conducive to signal filtering at the receiving end and achieving high-reliability data transmission [
15]. Signal design includes constant-envelope continuous phase synchronization signal design and OFDM transmission signal design. The details are as follows:
3.1. Design of a Constant-Envelope Continuous Phase Synchronous Signal
Constant-envelope continuous signals have good correlation performance, which is conducive to achieving signal synchronization detection. A Minimum-Shift Keying (
MSK) signal is used as the synchronization head; the expression of this signal is shown in Equation (1):
where
is the carrier angular frequency;
is the input code unit;
is the code width; and
is the initial phase of the m-th code unit, which is invariant across a code unit width. According to the continuity of the phase, that is, the total phase of the previous code unit, the terminal is equal to the initial total phase of the next code unit (
), which must satisfy Equation (2):
Therefore, Equation (3) is obtained.
MSK is a special type of binary-frequency shift keying with the following characteristics [
16]:
- a)
A constant-envelope that can be amplified using nonlinear amplitude saturation devices;
- b)
At the time of symbol conversion, the phase of the signal is continuous, which means that the waveform of the signal does not undergo sudden changes;
- c)
The instantaneous frequency is always one of two values, as shown in Equation (4)
The modulation index is .
The synchronous head design adopts an orthogonal modulation method. First, the data to be transmitted are differentially encoded, then serially converted into I-channel and Q-channel signals. Then, each channel is multiplied by the weighting functions
and
and carriers
and
. Finally, the two channels are summed to obtain a complete MSK signal. The signal modulation box is shown in
Figure 1 [
17].
The waveform of two signals with different frequencies in the MSK signal, that is, when
is +1, the frequency is
; when
is −1, the frequency is
. The final designed waveform is shown in
Figure 2, with
at 10 Hz,
at 5 Hz, the carrier at 7.5 Hz, the bit cycle at 0.1 s, and the sampling rate at 500 Hz, generating a constant-envelope signal curve with 1600 points [
18].
The constant-envelope waveform generated by the MSK modulation method has a constant amplitude and a continuous phase, and the maximum and average values of the waveform envelope are equal everywhere. Therefore, the peak-to-average ratio of this waveform is 0, which has strong stability and anti-interference ability [
19].
3.2. OFDM Transmission Signal Design
Orthogonal frequency division multiplexing (OFDM) technology is applied to ground-penetrating communication. Essentially, it converts signals from high-speed serial signals to low-speed parallel signals through serial-to-parallel conversion, with each parallel signal corresponding to a subcarrier with a different frequency. However, the attenuation amplitude of signals at different frequencies varies in terrestrial channels, so OFDM not only exhibits multicarrier but also multichannel performance in ground-penetrating communication. After the serial parallel conversion, the symbol period in each subchannel is increased relative to the original symbol, which reduces the impact of time dispersion caused by the multipath delay expansion of the terrestrial wireless channel on the entire ground-penetrating system. During current field OFDM modulation, appropriate protection intervals are inserted between OFDM symbols, and the protection interval is greater than the maximum delay extension of the terrestrial wireless channel, thereby minimizing intersymbol interference caused by multipath interference in the terrestrial channel. In waveform design, the system uses a cyclic prefix as the protection interval for the ground-penetrating communication system [
20,
21].
Figure 3 shows the OFDM signal workflow.
The principle of OFDM is to convert signal symbols from high-speed serial signals to low-speed parallel signals through serial–parallel conversion. Each parallel signal corresponds to subcarriers with different frequencies. The earth channel has different attenuation amplitudes for signals of different frequencies, so OFDM is both multicarrier and multi-channel in TTE communication.
In signal design, cyclic prefix (CP) is used as the protection interval for the TTE system. In addition, OFDM embeds the transmission information into the carrier coefficients. Its carrier has orthogonality, and the spectrum before the carrier can overlap, which improves spectrum utilization.
The first step in generating OFDM signals is to perform serial-to-parallel conversion. Before OFDM modulation, the ground-penetrating communication information is a serial data stream, and the source information is continuously transmitted as in ordinary communication systems. The spectrum of each source data symbol occupies the entire channel bandwidth. After converting the serial data stream into a parallel signal, multiple data symbols can be transmitted simultaneously, and different modulation methods can be used on each subcarrier. In order to improve anti-interference performance, the system uniformly selects QPSK as the modulation method for the subcarriers.
Figure 4 shows a constellation diagram of the binary bitstream at the transmitter after QPSK modulation [
22].
Assuming that the data symbol of each subcarrier is
,
is the number of subcarriers, and
i is the subcarrier sequence number, the OFDM symbol can be expressed as Equation (5):
Equation (5) and the Inverse Discrete Fourier Transform (IDFT) (Equation (6)) are consistent.
In the OFDM expression, can be regarded as the frequency-domain information () in the IDFT formula. Therefore, the generation of signals transmitted through the ground in OFDM is equivalent to performing IDFT operations on , and Inverse Fast Fourier Transform(IFFT) is a fast algorithm for IDFT. Therefore, IFFT is used to generate OFDM signals.
At present, the sampling rate (
) of the current field transmitter in the system is 1000 Hz, while the number of IFFT points is generally 2
N, and
is the number of subcarriers. According to the relationship between the minimum subcarrier interval and the sampling rate, as well as
:
To generate low-frequency ground-penetrating OFDM signals, it is necessary to reduce the subcarrier frequency interval and select 1024 points as the IFFT calculation length: .
Placing QPSK data on subcarriers 1–10 in the frequency domain ensures that the maximum frequency of the generated OFDM symbol is 9.766 Hz, not exceeding 10 Hz, achieving low-frequency OFDM signal design. The conjugate data of modulation information (
) are placed in a symmetrical position in the frequency domain to ensure that the generated signal after IFFT is a real-number signal.
Figure 5 shows the frequency-domain information of low-frequency ground-penetrating OFDM.
After performing the IFFT operation on frequency domain data, low-frequency OFDM symbols are generated.
Figure 6 shows a time-domain waveform carrying a 20 bit ground-penetrating low-frequency OFDM symbol.
Due to the multipath effect in the terrestrial channel, signals generate Inter Symbol Interference(ISI) and Inter Channel Interference (ICI) during transmission in the terrestrial channel. Therefore, the period of OFDM symbols is extended to achieve the goal of eliminating ISI and ICI.
Therefore, the tail of the quarter length of the low-frequency OFDM symbol through the ground is selected as the cyclic prefix, as shown in
Figure 7, which is the OFDM time-domain symbol with the added cyclic prefix.
4. Signal Detection and Recognition
Due to the special nature of information transmission in closed collapsed-earth electrodes, the received voltage in the information transmission system based on the current field of the earth electrode rapidly decays with increased penetration distance. When the penetration distance exceeds 100 m, the voltage value detected by the receiving end is extremely low. The purpose of weak signal detection at μV or even nV levels is to extract signals submerged in noise.
In order to verify the correctness of the signal design proposed in this article, we conducted on-site testing on a tunnel under construction in Shunyi District, Beijing, using the unfinished part in the middle. The tunnel is 8 m wide and 6 m high, and the unconnected part is 63 m long. The tunnel is 20 m under the ground, and the ground electrodes are 3 m apart. The ground electrode is 1.5 m deep into the ground.
A schematic diagram of the testing structure is shown in
Figure 8, and the testing site is shown in
Figure 9.
We verified the transmission feasibility and reliability of the signal design in actual tunnel application sites through testing experiments in tunnels under construction.
In the information transmission system of the earth electrode current field, due to the complexity of the electrical conductor of the collapsed medium, power-frequency interference and the earth current field form noise during signal reception, interfering with signal detection. Therefore, external electromagnetic environment testing was conducted to analyze interference characteristics and verify anti-interference technology.
Figure 10 shows the time-domain waveform of background noise when no signal is sent.
By performing an FFT, the frequency domain information of the noise signal can be obtained [
23], as shown in
Figure 11. It can be seen that the interference noise mainly comes from 50 Hz power-frequency interference and signal components near 0 Hz.
The designed constant-envelope synchronization curve is used as the transmitter signal input, and nine synchronization signals are sent at fixed intervals. The maximum voltage of the transmitter is 128 V, and the sampling rate of the receiver is 500 Hz. The standard output signal and the signal received by the receiver are shown in
Figure 12.
A frequency-domain FFT is performed on the received signal, as shown in
Figure 13. It can be seen that the frequency band of the transmitted signal is around 5–10 Hz, while the interference signals comprise DC interference and 50 Hz power-frequency interference, which are consistent with the background noise test.
By combining time-domain and frequency-domain information analysis, it can be concluded that the presence of DC interference and power-frequency interference leads to the submergence of useful signals.
Therefore, for the received signal, filtering processing is required, using a high-pass filter to filter out the DC component [
24] and a low-pass filter to filter out the power-frequency interference. The filtering effect is shown in
Figure 14. The above figure of
Figure 14 shows the time-domain waveform after passing through a high-pass filter, and it can be seen that the signal mean approaches 0. The below figure of
Figure 14 shows the time-domain waveform after passing through a low-pass filter, indicating that the high-frequency part of the signal has been filtered out.
Correlation calculations are performed between the synchronization curve and the intercepted receiving curve of the same length, with point-by-point sliding, and the correlation value of the point with the highest correlation value is recorded, obtaining the relevant results presented in
Figure 15.
Nine related peaks are identified, consistent with the number of times the sender sends a signal.
The waveform of the signal is shown in
Figure 16, with a constant-envelope signal as the synchronization head, taking into account the sudden change effect caused by the capacitance effect during transmission and adding a 200 point presignal. Finally, the 1024-point OFDM symbol, combined with a quarter-length cyclic prefix and the constant-envelope synchronization head, forms a 3080-point signal.
In the field experiment conducted on the simulated test site, the transmitter sent the signal at a sampling rate of 1000 Hz above the OFDM waveform five times, and the receiver received the signal at a sampling rate of 500 Hz. The received time-domain waveform is shown in
Figure 17.
After transforming the received signal into the frequency domain, as shown in
Figure 18, it can be seen that the interference is also DC interference, with 50 Hz power-frequency interference.
A high-pass filter is used to filter out the DC component, and a low-pass filter is applied to filter out power-frequency interference. The filtered signal is shown in
Figure 19, clearly identifying the five sent curves.
Similarly, using the correlation-matching algorithm, the synchronization head is slid point-by-point, and five obvious correlation peaks can be identified in
Figure 20.
After filtering and constant-envelope correlation detection matching at the receiving end, all five waveforms can be found, consistent with the number of transmitted synchronization signals. This indicates the synchronization effect of the MSK constant-envelope synchronization signal in the ground electrode current field, as shown in
Figure 21.
For the extracted OFDM signal, the FFT operation, which is completely opposite to the IFFT modulation method, is used for demodulation.
The demodulated frequency domain information (
) is shown in Equation (8), as shown in
Figure 22, indicating that there is no significant deviation in the constellation diagram.
QPSK demodulation is performed on
, and the complex information is mapped to binary bit stream information. Compared with the sent information, this information is completely correct, as shown in
Figure 23.
Due to the known data of the transmitted signal, the OFDM system can be represented as Equation (9):
where
is the channel response,
is a known transmission signal,
is the received pilot signal, and
is the noise vector. Therefore, least-squares (LS) channel estimation can be used to estimate parameter
and minimize the function.
The channel estimation value of the LS algorithm is expressed as Equation (10):
An IFFT operation is performed on the channel estimation value to obtain the estimated channel impulse response [
25], as shown in
Figure 24.
In
Figure 24, the presence of multiple pulses indicates the presence of delay extension in the channel. The root-mean-squared delay spread (
) can be expressed as Equation (11):
where
is 0.063 s. It is generally believed that the multipath fading of the channel follows a negative exponential distribution. Maximum delay extension (
) is four to six times that of
; therefore
is 0.252 s. When designing OFDM signals, the length of the cyclic prefix should be greater than this delay.
From the above tests, it can be seen that the OFDM signal modulated by 10 Hz QPSK achieves satisfactory penetration performance in outdoor ground-penetrating environments. The influence of ground electric field noise and power-frequency interference noise can be filtered out by FIR filters. The ground channel has an impact on the amplitude and phase of the OFDM subcarriers, which can cause misjudgment of the QPSK demodulation. Using known symbols for channel estimation and equalization can compensate for errors to a certain extent; however, residual bias will remain, and after compensating for residual phase tracking, the rotation effect of the constellation can be restored, enabling the received signal to undergo normal QPSK demodulation.
5. Conclusions
This article mainly focuses on the principle of ground electrode current field communication, signal design, signal detection, and recognition in the ground electrode current field penetration communication technology. Through reviewing relevant literature and analyzing experimental data, we have reached the following conclusions.
- (1)
The ground electrode current injection technology has good transmission performance through the ground and can achieve reliable underground communication. Its advantage is that it can still achieve good communication performance in harsh underground road conditions.
- (2)
Due to the particularity of the geological environment, the earth’s geological channels are very complex, and the attenuation of high-frequency electromagnetic waves is very severe. In system design, only lower frequency (10-20Hz) waves can be selected for transmission.
- (3)
Using the minimum frequency shift keying MSK signal as the synchronization head is conducive to achieving signal synchronization detection. Using OFDM signals as transmission signals is beneficial for reducing the impact of time dispersion caused by multipath delay expansion in terrestrial wireless channels on the entire ground penetrating system and eliminating inter symbol interference caused by multipath in terrestrial channels.
The underground engineering underground emergency rescue communication technology is of great significance. Firstly, this technology is widely used and can be applied in scenarios such as tunnels, coal mines, non coal mines, and large and medium-sized underground engineering. Secondly, the application of this technology can shorten rescue time, achieve scientific and precise rescue, and reduce rescue risks. Once again, this technology can also be applied in fields such as settlement monitoring, geological exploration, and groundwater exploration in underground engineering construction in the future. After years of continuous research, the project team has applied relevant technologies to develop and produce engineering prototypes. The communication frequency used is 10Hz, and the maximum voltage injected into the grounding electrode is 250V. The engineering prototype has achieved a maximum transmission distance of 1000 meters and a maximum communication capability of 24 characters of self generated information length. Based on existing research results and future application needs, we have reason to believe that the grounding electrode current field through ground communication technology will be widely applied and promoted in future research and practice, and contribute to human development.
Figure 1.
MSK signal modulation process.
Figure 1.
MSK signal modulation process.
Figure 2.
Constant-envelope synchronous head signal.
Figure 2.
Constant-envelope synchronous head signal.
Figure 3.
OFDM signal workflow.
Figure 3.
OFDM signal workflow.
Figure 4.
QPSK-modulated constellation diagram.
Figure 4.
QPSK-modulated constellation diagram.
Figure 5.
Frequency domain information of low-frequency ground-penetrating OFDM signals.
Figure 5.
Frequency domain information of low-frequency ground-penetrating OFDM signals.
Figure 6.
Time-domain waveform of a ground-penetrating low-frequency OFDM symbol.
Figure 6.
Time-domain waveform of a ground-penetrating low-frequency OFDM symbol.
Figure 7.
A low-frequency ground-penetrating OFDM symbol with a cyclic prefix.
Figure 7.
A low-frequency ground-penetrating OFDM symbol with a cyclic prefix.
Figure 8.
Schematic diagram of the testing structure.
Figure 8.
Schematic diagram of the testing structure.
Figure 10.
Background noise of simulated testing ground.
Figure 10.
Background noise of simulated testing ground.
Figure 11.
Background noise frequency-domain information.
Figure 11.
Background noise frequency-domain information.
Figure 12.
The standard signal sent by transmitter and the signal received by the receiver.
Figure 12.
The standard signal sent by transmitter and the signal received by the receiver.
Figure 13.
Frequency domain of the received signal.
Figure 13.
Frequency domain of the received signal.
Figure 14.
Filter results.
Figure 14.
Filter results.
Figure 15.
The result is consistent with the number of transmitter transmissions.
Figure 15.
The result is consistent with the number of transmitter transmissions.
Figure 16.
OFDM signal with constant-envelope synchronization head.
Figure 16.
OFDM signal with constant-envelope synchronization head.
Figure 17.
Time-domain waveform of OFDM signal reception.
Figure 17.
Time-domain waveform of OFDM signal reception.
Figure 18.
Frequency domain of the received OFDM signal.
Figure 18.
Frequency domain of the received OFDM signal.
Figure 19.
Time-domain waveform of the OFDM signal after filtering.
Figure 19.
Time-domain waveform of the OFDM signal after filtering.
Figure 20.
The result is consistent with the number of transmitter transmissions.
Figure 20.
The result is consistent with the number of transmitter transmissions.
Figure 21.
Extracting OFDM signals.
Figure 21.
Extracting OFDM signals.
Figure 22.
Receiver OFDM signal constellation diagram.
Figure 22.
Receiver OFDM signal constellation diagram.
Figure 23.
OFDM signal transmission bit and reception bit.
Figure 23.
OFDM signal transmission bit and reception bit.
Figure 24.
Channel impact response.
Figure 24.
Channel impact response.