Abstract

Radio wave propagation in confined spaces is consequent upon the reflections of boundaries; thus, the radiation characteristics of the antenna have significant influence on the field coverage in the confined space. This paper investigates the effects of antenna parameters on field coverage characteristics in a tunnel environment. A modified modal method is proposed to analyse the wave propagation properties along the tunnel. The relationships between the amplitudes of modes and the antenna parameters, including the beam width, beam direction, and antenna location, are analysed. The results indicate that by properly selecting the antenna parameters, optimum field coverage in tunnel environments can be realized.

1. Introduction

Effective and reliable electromagnetic wave coverage in confined environments is important to the quality of wireless communications. In general, the discrete antenna system and the leaky coaxial cable (LCX) system are two types of infrastructure widely used to generate wave coverage in confined spaces [1, 2]. The discrete antenna system is more flexible and efficient and is less expensive compared with the LCX system. However, most available discrete antennas suffer from unsatisfying performance, which may result in the interruption of communications in long and narrow confined areas, such as tunnels [3]. To overcome this shortcoming, the characteristics of the wave coverage generated by discrete antennas in tunnels should be investigated.

Radio propagation in tunnels has been studied for decades. Different types of methods have been employed to predict the propagation characteristics. Due to the bounding effect of the tunnel walls, the propagation characteristics in tunnels are different from those in open space. In [4, 5], the theoretical modal analysis is applied to build the propagation model, which has high accuracy when the receiver is far away from the transmitter. Nevertheless, the method has difficulty in taking the antenna parameters into account. The ray tracing method [6, 7] is another powerful way to address radio propagation in complex environments. The antenna parameters can be considered in this method, but it is difficult to explain the physical mechanism of the propagation procedure directly. Full wave calculation methods [8, 9], including finite difference time domain (FDTD), finite element method (FEM), and integral equations, can contain the effect of different antenna parameters, and precise results can be obtained. However, these methods are time and memory consuming for large-scale problems. Recently, the parabolic equation method [10, 11] has become popular in solving the radio propagation in tunnels as it can determine the path loss and fading characteristics in long tunnel environments in an acceptable amount of time. However, in applying this method, it is difficult to accurately model high order modes, especially in the vicinity of the transmitting antenna.

By using these reported methods, channel parameters, which are essential in wireless communication applications, including path loss, time delay, and break point, can be studied. In [4, 5], the field composition and propagation constants of the modes in a hollow rectangular dielectric waveguide were investigated. The propagation characteristics in the tunnel and the relationship between the beak point and the operating frequency, antenna gain, and dimensions of the tunnel were analysed in [12]. However, most available studies generally consider the antenna and the wave propagation independently. According to the results in [1315], the wave propagation is influenced by the properties of the antennas used in a confined space [1315].

In this paper, the field coverage instead of wave propagation is used to describe the wave characteristics in tunnels where the antenna is located. As described in [16], the ray tracing method is mathematically equivalent to the modal method in the simulation of radio propagation in straight rectangular tunnels. However, the modal method in [16] is derived based on an isotropic source, which means that it is not suitable for the case of directional antennas. Hence, a modified modal method is proposed in this paper. By applying this method, the effects of antenna parameters (beam width, beam direction, and location of antennas) on the field coverage in tunnels can be successfully studied. The results achieved in this paper are valuable for the design of antennas to be used in confined spaces. Additionally, the amplitudes of the modes can be obtained, which is helpful in analysing the mechanism of field coverage.

The rest of this paper is organized as follows: in Section 2, the novel modal method is depicted in detail; in Section 3, the effect of antenna parameters on field coverage is presented; and, in Section 4, the conclusions are summarized.

2. Methodology

The model and parameters of the straight hollow rectangular tunnel considered in this paper are shown in Figure 1. The tunnel is directed along the -axis, and the origin of the coordinates is located in the centre of the tunnel cross section. Based on the ray tracing theory, the electric field at the point (, , ) in the tunnel can be expressed by the summation of the rays radiated by all the images of the source [17]: In the above equation, is the magnitude of the transmitted electric field, is the wave number, is the distance between the receiver and the image , and is expressed bywhere and . and are the coordinates of the centre of the transmitting antenna in the plane. It should be noted that only the case with vertical polarization is considered in this paper. in (1) represents the reflection coefficients on the perpendicular and parallel walls, respectively, and can be calculated approximately by [17]where and are the reflection angles on the perpendicular and parallel walls, respectively; is the quantity related to the effective surface impedance defined by , ; and are the complex relative permittivities of the vertical and horizontal walls, respectively.

in (1) represents the antenna radiation pattern. To simulate the radiation properties of the antenna, the Gaussian pattern is used as the excitation in this paper, which is given by where and represent the beam angles in the azimuth and elevation directions, respectively, and represent the beam widths in the azimuth and elevation directions, respectively, and is a constant equal to 0.7213. Because the tunnel is set along the direction, the range of and in the forward direction is (, ), and when the tunnel is sufficiently long, they can be approximated asThen, (1) can be rewritten asSimilar to that in [16], (6) can be written aswhere , , , and . and for even numbers of and , respectively, and and for odd numbers of and , respectively. For the case in this paper, is expressed asThe definitions of the variables in (8) are the same as those defined in [16].

The Poisson Summation Formula (PSF) and saddle-point method can be utilized to convert (7) into the mode formula [16, 18], and then (7) can be expressed as whereand is the wave number of free space. and are mode numbers corresponding to the horizontal and vertical directions, respectively. As derived in the Appendix, (9) can be simplified toNote that (11) is exactly the weighted summation of the eigenfunctions of the propagation modes, whereare the classic modal attenuation and propagation constants [4, 5]. The weight coefficients in (11) can be expressed asFrom (14), it is observed that the weight of each eigenfunction is exactly the mode amplitude in the excitation plane, and it is related to the beam width of the antenna, azimuth and elevation angles of the beam, and the location of the antenna. Hence, the effect of the antenna parameters on the radio propagation and field coverage characteristics can be analysed by modes. In actuality, all the antenna radiation patterns can be taken into account in rectangular tunnels by utilizing this method, indicating that the radiation pattern can be analytically expressed. The excitation mode amplitude can be expressed aswhere represents and represents .

In this paper, the Gaussian radiation pattern is used as a special case to demonstrate the efficiency of the method. Figures 2(a) and 2(b) give the comparison between the received powers obtained with different beam tilt angles and frequencies using the modified modal method proposed in this paper and the ray tracing method, respectively. In the calculation, the transmitting antenna is located in the centre of the cross section of the tunnel with dimensions of 6 m × 6 m × 1000 m and an operating frequency of 900 MHz. The received power is calculated by [19] where is the received field by a vertical polarized half wave dipole, is the wave impedance of free space, and is the wavelength. To ensure the accuracy of the approximation, the maximum orders of and in the modal method are set to 40. As shown in Figure 2, the result obtained from the new modal method agrees well with that of the ray tracing method. To further clarify the accuracy of the proposed modal method, the measured result in [20] is compared with that of this method. As shown in Figure 2(c), the result of the new modal method agrees well with the result measured in the Massif Central tunnel. The dimensions of the equivalent rectangle, which are suitable for simulating this tunnel, are listed in Table 1.

3. Results and Analysis

In this section, the effect of antenna parameters on the received power along a rectangular tunnel is analysed by utilizing the above modified modal method. The rectangular tunnels in this paper are set along the direction, and the parameters of the tunnels are listed in Table 1. The relative permittivity and conductivity of the tunnel wall are shown in Table 1, which imitate the concrete. It is assumed that the antenna efficiency is 100% and the operating frequency is 900 MHz. The transmitting and receiving antennas are all vertically polarized.

3.1. Effect of Tunnel Dimensions

The curves of received power in tunnels A, B, and C are shown in Figure 3. To obtain the general rule of the effect of the tunnel dimension parameters on the field coverage characteristics, the transmitter and receiver are located in the centre of the cross section of the tunnels. From Figure 3, it is observed that the variations of the received power in the tunnels with the same width are similar. This similarity is because in the vertical polarization case, the dominant propagation modes are modes, and the phase constants of those modes are related to the tunnel width. Thus, if the tunnel width is fixed, the superposed manner of the modes is the same. If the tunnel widths are different, the positions of peaks and nulls on the variation curves of the received power are changed, and the pseudo-periodicities become larger.

In the tunnels, the number of modes, which significantly affects the received power, is an important parameter for interpreting the field coverage characteristics. The number of “active” modes, (10%), is defined as the number of modes with power equal to at least 10% of the power in the most powerful mode [21]. The variations of versus distance when transmitter is in the centre of the different tunnels are shown in Figure 4. The number of active modes helps for interpreting the fast variations of the received power around its mean value. In a small tunnel, the number of active modes is much smaller than in a large tunnel (as shown in Figure 4), so less modes can be reserved in the propagation process and arrive at the far region. Therefore, due to the superposition of different modes, the received power in the small tunnel suffers from less fluctuation as shown in Figure 3.

3.2. Effect of the Beam Width on Field Coverage

Considering that a rectangular tunnel with a height of 6 mm and width of 8 mm is common in practice, tunnel B is taken as an example to analyse the effect of antenna parameters on the field coverage characteristics. To evaluate the effect of the antenna beam width on the received power, the antenna is located in the centre of the cross section of the tunnel and points to the axial direction of the tunnel. The beam width of the transmitting antenna is set to 15°, 30°, 45°, and 60°, respectively. Figure 5 shows the received powers along the tunnel axis for different beam widths. As shown in Figure 5, the curves of the received powers generally have similar shapes, but, in the near region, the received power fluctuates significantly, leading to large differences in terms of the shapes of the curves, whereas, in the far region, the fluctuations become smaller. It is also shown that a narrower beam width of the transmitting antenna leads to smaller fluctuations in the near region and higher received power in the far region. In addition, the received power levels for different beam widths are similar to each other in the near region but significantly different in the far region. This phenomenon can be explained with the help of the proposed modal method. Figure 6 shows the variation of the number of active modes versus distance to the transmitter, which is excited by different antenna beam widths. It is shown that, for the narrower beam, the number of active modes is less, and the power is mainly transferred into the lower order modes. For the wider beam, the number of active modes is much higher than that of the narrower beam and more power is transferred into the higher order modes. From (12), it is obvious that the higher order modes have larger attenuation constants, so they decrease rapidly in the near region. Alternatively, for lower order modes, the attenuation constants are small, so they decline slowly. As a result, the received power for the wider beam presents rapid decay and deep fluctuation in the near region and becomes stable in the far region. This phenomenon can also be interpreted from Figure 6 in which the curve of the antenna with a beam width of 15° declines slowly, whereas the curve of the antenna with a beam width of 60° decreases rapidly. Thus, different mode compositions cause different received power distributions, which cannot be explained simply by the antenna gain. It is also noticed that increasing the number of active modes can improve the capacity of the diversity antenna systems, but it also reduces the received power level. Therefore, there is a trade-off between the number of active modes and the received power level.

In formula (14), there are two parameters, and , related to the transmitting antennas. Because it is obvious that the antenna beam pointing to the tunnel axis exhibits the highest level, installing an antenna with a slightly tilted angle is a common method. Thus, and are usually fixed to 0°.

3.3. Effect of Antenna Locations on Field Coverage

Many studies have investigated the effect of antenna position on the path loss characteristics from the angles of measurement or simulation [22, 23]. However, few studies consider this issue theoretically. Figure 7 shows the received power along the tunnel axis for different transverse locations of transmitting antenna in the tunnel cross section. It can be observed from Figure 7 that when the antenna moves vertically from the centre to the ceiling, the received power level reduces markedly, but the variation characteristic along the tunnel axis changes slightly. When the antenna moves horizontally from the centre to the side wall, the variation characteristic of the received power along the tunnel axis changes significantly, and the positions of peaks and nulls markedly shift. When the transmitter is off centred, a large number of high order modes are excited and carry a large amount of energy. is not the mode with highest energy compared with these modes, different from the case of excitation by a centred transmitter. The high order modes rapidly attenuate in the process of propagation in the tunnel and take away a lot of energy, which is the reason why the received power level reduces markedly when the transmitter is off centre as shown in Figure 7.

Figure 8 shows three cases of the variation of the active mode number when the transmitting antenna is located in the centre, close to the side wall and to the ceiling, respectively. It can be found that when the transmitter is off centre many modes can be reserved in the tunnel compared with the case that transmitter is in the centre, which will result in the larger fluctuation of the received power.

To further understand the excited modes when the transmitter is located in different positions, Table 2 shows the order of the mode with the maximum amplitude in the planes of and 1000. As Table 2 shows, when the antenna moves vertically to the ceiling, modes related to index become significant. Because the high order modes with large have larger attenuation constants, they rapidly decay along the tunnel, and finally only the fundamental modes remain in the far region. In contrast, when the antenna moves horizontally to the side wall, modes related to index become significant, which have smaller attenuation constants, so many can be reserved in the propagation process and arrive at the far region. The normalized amplitude distributions of the active modes at the plane of  m tested with different antenna locations are shown in Figure 9. When the transmitter moves horizontally, , , and are excited and can propagate for long distance, which cannot be excited when the transmitter move vertically. Therefore, due to the superposition of different modes, the received power for the case of the horizontally moving antenna suffers from larger fluctuation and deeper nulls, as shown in Figure 7. In addition, because the phase variations of modes for antennas moving horizontally are substantially different from those of antennas moving vertically, the locations of peaks and nulls change significantly.

From the above analysis, the effect of the antenna parameters (beam width and location) on the received power characteristics can be clearly observed. However, it is worth mentioning that the effects of these parameters are coupled together, and the effects should be considered simultaneously. According to the results, for the vertical polarization case, increasing the amplitudes of fundamental and lower order modes related to index can extend the communication region in the tunnel. Because of the condition of practical application, the transmitting antenna has to be installed on the side wall. Thus, for the vertical polarization case, the narrow beam width antenna with no tilt or slight tilt is preferable.

From the above analysis, it can be found that by only changing the location of the transmitter it is hard to realize the optimum field coverage in the whole tunnel because deep fading along the tunnel still exists. Multitransmitters could realize the optimum field coverage in a line along the tunnel, but not in the whole tunnel. Unlike multitransmitter systems, the spatial diversity system with multireceivers is an effective and flexible way to achieve uniform field coverage because it could utilize the irrelevance of the different positions of the tunnel to obtain an average level. Improving the diversity of antenna systems in tunnels, which has been a hot research topic in recent years [24, 25], is an efficient way to determine the optimal locations of transmitters and receivers to minimize the bit error rate and achieve higher capacity.

4. Conclusions

This paper investigates the field coverage characteristics in different tunnel environments for different antenna parameters, including beam width and antenna location. The analysis is based on a newly proposed, modified modal method, which is mathematically equivalent to the ray tracing method, for the simulation of the wave propagation in a tunnel environment. The relationship between the mode amplitudes and the antenna parameters is obtained and analysed. The results indicate that the antenna parameters have a significant effect on field coverage. The investigation in this paper offers insight into understanding and interpreting field coverage characteristics in tunnels.

Appendix

The last summation in (9) can be expanded asfrom which (12) and (13) can be obtained. The remaining summation formula in (A.1) can be further simplified by the Euler formula as follows:where the second term on the right side can be expressed bySubstituting (A.3) into (9) and considering the summation of indexes and asthen the complex exponential functions corresponding to summation of index in (A.4) can be converted to sinusoidal functions based on Euler’s formula asSimilar to (A.5), the summation corresponding to index can be expressed asBecause of the symmetry of the images, and can be exchanged. By substituting (A.4), (A.5), and (A.6) into (11), then formula (14) can be obtained.

When the antenna pattern is unknown, (A.1) can been expressed asSimilar to (A.5) and (A.6), the final summation for the unknown radiation pattern case can be given asBecause of the symmetry of the images, and can be exchanged, and then (15) can be obtained.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by NSFC Projects under Grant no. 61331002 and in part by the National Key Basic Research Project under Grant no. 2013CB328903.