Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
PADL: A Modeling and Deployment Language for Advanced Analytical Services
Next Article in Special Issue
Design, Implementation, and Empirical Validation of an IoT Smart Irrigation System for Fog Computing Applications Based on LoRa and LoRaWAN Sensor Nodes
Previous Article in Journal
A Numerical Model for Experimental Designs of Vibration-Based Leak Detection and Monitoring of Water Pipes Using Piezoelectric Patches
Previous Article in Special Issue
Fabricating a Portable ECG Device Using AD823X Analog Front-End Microchips and Open-Source Development Validation
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radio Wave Propagation and WSN Deployment in Complex Utility Tunnel Environments †

by
Mikel Celaya-Echarri
1,
Leyre Azpilicueta
1,*,
Peio Lopez-Iturri
2,3,
Imanol Picallo
2,3,
Erik Aguirre
2,
Jose Javier Astrain
4,
Jesús Villadangos
4 and
Francisco Falcone
2,3
1
School of Engineering and Sciences, Tecnologico de Monterrey, 64849 Monterrey, NL, Mexico
2
Electric, Electronic and Communication Engineering Department, Public University of Navarre, 31006 Pamplona, Navarra, Spain
3
Institute of Smart Cities, Public University of Navarre, 31006 Pamplona, Navarra, Spain
4
Mathematical Engineering and Computer Science Department, Public University of Navarre, 31006 Pamplona, Navarra, Spain
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in Celaya-Echarri, M.; Azpilicueta, L.; López-Iturri, P.; Picallo, I.; Aguirre, E.; Astrain, J.J.; Villadangos, J.; Falcone, F. Wireless Channel Characterization and System Analysis of Complex Utility Tunnel Environments. In Proceedings of the 6th International Electronic Conference on Sensors and Applications, 15–30 November 2019; Available online: https://ecsa-6.sciforum.net.
Sensors 2020, 20(23), 6710; https://doi.org/10.3390/s20236710
Submission received: 29 October 2020 / Revised: 17 November 2020 / Accepted: 19 November 2020 / Published: 24 November 2020
Figure 1
<p>Schematic view of the structure of this work. WSN: wireless sensor network and 3D-RL: three-dimensional ray-launching.</p> ">
Figure 2
<p>Principle of operation of the in-house 3D-RL algorithm.</p> ">
Figure 3
<p>Underground utility tunnels under analysis. (<b>a</b>) Tunnel location on the map of the Old Town of Pamplona. (<b>b</b>) Picture of the tunnel where the measurements were taken. (<b>c</b>) Picture of a newer tunnel built in the neighborhood of Lezkairu in Pamplona.</p> ">
Figure 4
<p>Schematic representation of the simulation scenario implemented to analyze the radio propagation losses within the underground utility tunnels.</p> ">
Figure 5
<p>Comparison of different cases of simulation and measurement results for the service tunnel section under analysis.</p> ">
Figure 6
<p>Radiofrequency (RF) power distribution for the XY bidimensional plane at 2-m height for 868 MHz.</p> ">
Figure 7
<p>Power delay profile at a given cuboid located at the center (60 m, 1 m and 2 m) at 868-MHz frequency.</p> ">
Figure 8
<p>Root mean square delay spread (RMS DS) estimation within the utility tunnel for 2-m height for 868-MHz frequency.</p> ">
Figure 9
<p>Coherence bandwidth within the considered utility tunnel for a 2-m height for 868-MHz frequency.</p> ">
Figure 10
<p>Bidimensional plane of the RF power level distribution at a height of 3 m for (<b>a</b>) ZigBee 868 MHz and (<b>b</b>) LoRaWAN.</p> ">
Figure 11
<p>Bidimensional plane of the RF power level distribution at a 1.5-m height for the case without human presence.</p> ">
Figure 12
<p>Bidimensional planes of the RF power level distribution at a 1.5-m height for different human being locations along the tunnel.</p> ">
Figure 13
<p>RF power level distribution differences between the case without a human presence and the cases with human beings at different locations along the tunnel.</p> ">
Figure 14
<p>Inclusion of 3 computational human body models in the considered scenario for the 3D-RL simulation, emulating a maintenance team working in the utility tunnels.</p> ">
Figure 15
<p>RF power level distribution difference between the case without people and the cases with a work team of 3 people at different locations along the tunnel.</p> ">
Figure 16
<p>Bidimensional plane at a 3-m height representation of the RF power level distribution after applying the 20-dB margin due to the human presence for (<b>a</b>) ZigBee 868 MHz and (<b>b</b>) LoRaWAN.</p> ">
Figure 17
<p>Schematic view of the proposed ZigBee 868-MHz WSN deployment.</p> ">
Figure 18
<p>Architecture of the monitoring tool.</p> ">
Figure 19
<p>Time series monitoring by means of Grafana (iFrames).</p> ">
Figure 20
<p>Web-based monitoring tool.</p> ">
Versions Notes

Abstract

:
The significant growth of wireless communications systems in the last years has led to the adoption of a wide range of applications not only for the general public but, also, including utilities and administrative authorities. In this context, the notable expansion of new services for smart cities requires, in some specific cases, the construction of underground tunnels in order to enable the maintenance and operation works of utilities, as well as to reduce the visual impact within the city center. One of the main challenges is that, inherently, underground service tunnels lack coverage from exterior wireless communication systems, which can be potentially dangerous for maintenance personnel working within the tunnels. Accordingly, wireless coverage should be deployed within the underground installation in order to guarantee real-time connectivity for safety maintenance, remote surveillance or monitoring operations. In this work, wireless channel characterization for complex urban tunnel environments was analyzed based on the assessment of LoRaWAN and ZigBee technologies operating at 868 MHz. For that purpose, a real urban utility tunnel was modeled and simulated by means of an in-house three-dimensional ray-launching (3D-RL) code. The utility tunnel scenario is a complex and singular environment in terms of radio wave propagation due to the limited dimensions and metallic elements within it, such as service trays, user pathways or handrails, which were considered in the simulations. The simulated 3D-RL algorithm was calibrated and verified with experimental measurements, after which, the simulation and measurement results showed good agreement. Besides, a complete wireless sensor network (WSN) deployment within the tunnels was presented, providing remote cloud data access applications and services, allowing infrastructure security and safety work conditions. The obtained results provided an adequate radio planning approach for the deployment of wireless systems in complex urban utility scenarios, with optimal coverage and enhanced quality of service.

1. Introduction

The evolution of wireless communication has led to their adoption in a wide range of applications, including utilities and administrative authorities. In addition, the growth in population, as well as in the services offered in cities, has led, in some cases, to the construction of underground tunnels in order to reduce the visual impact within the city center, as well as facilitating maintenance and operation works of utilities. In recent years, the underground sector has seen the use of wireless networking for certain projects, i.e., tunneling construction [1], but they do not make full use of the benefits offered by wireless communications. The utilization of wireless networking in the underground environment offers the potential for major benefits in productivity and safety. Currently, most of the communications infrastructure in operational tunnels such as transport tunnels are provided by cables. However, if it is an infrequently accessed tunnel, a wired communication infrastructure will not generally be available; therefore, the maintenance, upgrade and inspection of those tunnels will rely on wireless communications [2]. Nevertheless, one of the issues that local authorities have considered is the fact that, inherently, underground service tunnels lack coverage from exterior wireless systems, such as mobile networks or municipal wireless local area networks (WLAN), which can lead to potentially dangerous situations for maintenance personnel working within the tunnel. Therefore, wireless coverage should be deployed within the underground installation in order to guarantee real-time connectivity.
In emergency situations, such as fires, explosions or earthquakes, the normal wired communication infrastructure of an underground environment could be inoperable. Due to that, sometimes, it is necessary for the emergency services to have a wireless communication infrastructure following a serious incident. Rescue services normally are equipped with very high-frequency (VHF) or ultra-high-frequency (UHF) radio equipment, which is effective aboveground, but such radios have a very poor range in tunnels. The use of leaky feeder cables [3] is usually used to extend the coverage of VHF and UHF systems. However, the use of microwave frequencies provides a much greater range, without expensive leaky feeder systems, in such environments, and, also, they offer greater bandwidths. In addition, leaky feeders are also susceptible to damage during a fall or explosion. In [4], a system including both wired and wireless communications was described that enabled voice communication, as well as the exchange of photographs, video streams or other data, inside the tunnel. In [5], a wireless-based personnel tracking system to track people in an underground environment was presented, showing that the system can reduce the time taken for rescuers to reach any casualties.
The characterization of wireless channel behavior within complex heterogenous underground environments, such as utility tunnel systems, enables to perform link-balance estimations, which lead to coverage/capacity estimations. In this sense, an adequate estimation of losses is relevant in order to characterize the complex operation conditions, such as those given by high-obstacle density, with a high percentage of non-line-of-sight conditions and multipath propagation buildup owing to scatterers distributions within the scenarios under test or the effect of diffuse scattering, among others. This information is relevant in order to perform radio-planning tasks, which are required in order to design the network topology, node configuration and device-related parameters, such as antenna location/type. A radio-planning analysis was studied for different types of wireless networks (e.g., broadcasting systems, trunk radios, public land mobile systems, microwave transport links, etc.), considering aspects such as multi-variable optimization (energy consumption, node density and transmission rate as a function of the coverage area under analysis) and the impact of elements such as user density variation, estimation of cell edge conditions or interference background levels [6,7,8,9,10,11].
With the progressive adoption of wireless sensor networks, radio-planning tasks have been adapted to optimize network configuration [12,13,14,15,16,17,18]. In this case, the specific considerations given by the wireless sensor network (WSN) operation have to be taken into account, such as the random node location, the existence of predominant non-line-of-sight operation conditions or the inherent energy consumption and processing constraints, among others. In this way, a wireless channel analysis can be effectively employed in order to assess elements such as cooperative MAC schemes, optimized routing protocols or specific data aggregator policies in order to optimize the overall system operation. One of the main challenges in this case is to consider accurate radio channel conditions, given the large variability of node location and configuration within the wireless network layout, usually leading to the use of empirical/analytical-based channel estimation models. Radio propagation has been studied for different wireless systems within the interiors of galleries and tunnels [19,20,21,22], as well as for different tunnel geometries [23]. Other approaches have developed diffraction-based models, such as the alternate direction implicit parabolic equation method [24], boundary element method [25], numerical calculation of a vector parabolic equation [26] or application of the perturbation theory, among others. Another approach is based on the use of deterministic methods, such as ray launching (RL) or ray tracing (RT). One of the main issues of these approaches is the inherent divergence within the RL calculation as the distance increases between the transmit and receive antenna. In this case, solutions in scaled frequency models have been given by reconstructing the initial wave front of the source, achieved by ray density normalization [27]. One of the approaches adopted to perform wireless characterization within tunnels has been the adoption of dielectric waveguide models, in which the tunnel is emulated by a waveguide with dielectric walls [28]. These models, however, do not take into account the additional complexities, which appear in underground service tunnels, such as ducts, walkways or distribution trays, among others. Complete modeling can be achieved by employing full-wave electromagnetic simulation techniques, such as the finite difference time domain (FDTD) or frequency domain finite element techniques [29,30]. However, the computational cost of these techniques can be very high, due to the large size of the tunnel structures as compared to the wavelengths considered in conventional wireless communication systems. Moreover, the tunnels analyzed are usually devoted to road traffic or trains, with inherently larger cross-sections and lower structural complexity as compared with utility tunnels, with dense obstacle distributions.
In this paper, an in-house developed three-dimensional ray-launching (3D-RL) code is employed to analyze the wireless propagation in a real urban utility tunnel scenario at 868-MHz operation frequency. The aim of the article is a thorough analysis of the real results of measurements in the tunnels in order to obtain the optimal parameters for the simulation’s 3D-RL algorithm, so that, after their application, the simulated and experimental results show good agreement. By adequately modeling the tunnel gallery structure, a precise characterization for the complete extension of the tunnel system can be obtained, which can aid in the deployment of wireless communication and monitoring systems. Furthermore, a complete radio wave propagation analysis is performed, and a WSN deployment is proposed for the considered real underground utility tunnel system gallery. Finally, an infrastructure security and a work safety monitoring application prototype is developed in order to warrant work safety legislation.
The paper is organized as follows: the 3D-RL deterministic employed technique, as well as the considered real and simulated underground scenario are described in Section 2. Section 3 presents the radio channel analysis and the required convergence analysis due to the particular characteristics of the scenario. This analysis was performed compared with experimental measurements in the real utility tunnel in order to obtain the optimal simulation parameters. Once the optimal parameters were obtained, large- and small-scale simulation results within the tunnel were presented in order to characterize the radio wave propagation within this particular environment. In Section 4, the WSN deployment issues are discussed, and the infrastructure security and work safety monitoring application prototype are presented. Finally, the conclusions are given in Section 5. A schematic view of the structure of this paper is presented in Figure 1.

2. Materials and Methods

In this section, the in-house deterministic technique based on the 3D-RL approach is described, as well as the considered real and simulated underground utility tunnel scenario used for this analysis.

2.1. 3D Ray-Launching Technique

The in-house implemented RL algorithm, based on MATLAB programming environment, was used to analyze radio wave propagation in the considered underground tunnel scenario. The main difference between RL and classical RT methods is that, in the RL technique, the principle is that rays are launched in a solid angle and the rays arriving at the receiver are considered the true path. However, in classical RT techniques, the image of the transmitter or of the receiver is calculated, to be used to find the reflected and refracted paths.
The RL algorithm is based on geometrical optics (GO) and the geometrical theory of diffraction (GTD). To avoid the nonfield regions predicted by GO in the shadow areas, the diffracted rays are considered with the GTD and its uniform extension, called the uniform GTD (UTD). These diffracted rays are introduced to avoid the abrupt field changes and consider the proper field corrections. Electromagnetic phenomena such as reflection, refraction and diffraction are considered. The principle of the algorithm is that hundreds or thousands of rays are launched for a fixed position of the transmitter, which is an input parameter of the algorithm. The radiation pattern of the transmitter and receiver antennas can be introduced, with the aid of the elevation angle θ and the azimuth angle φ, considering the spherical coordinate system. Each ray has a different path in the considered scenario, and all the parameters related with its path are stored in a matrix. All the rays that achieve the receiver are considered in the calculation of the electromagnetic field at that point. The whole scenario is divided into a three-dimensional array of cuboids, with a different spatial resolution in the x-axis, y-axis and z-axis, depending on the specific characteristics and dimensions of the considered scenario. The basis is that, when a ray impacts an obstacle or a wall in the scenario under analysis, reflected and refracted rays are created, following Snell’s law. When a ray hits an edge, a new family of diffracted rays are created, according to the UTD. The methodology is shown in Figure 2.
As input parameters of the algorithm, the frequency of operation, number of multipath reflections and spatial and angular resolution are considered. In the simulations, the material properties of all the obstacles within the scenario are also considered, taking into account the conductivity and relative permittivity at the frequency range of operation of the system under analysis. The in-house implemented 3D-RL algorithm has already been validated in the literature for different complex scenarios for frequencies up to 6 GHz. Some scenarios that have been considered have been typical office scenarios [31,32], sports environments [33], vehicular communications [34] or WSNs in smart cities [35]. The detailed description of the algorithm can be found in [36].

2.2. Scenario Description

In order to perform the radio wave propagation analysis within the tunnel, two scenarios were considered for the 3D-RL simulations: a real underground utility tunnel scenario located in the Old Town of Pamplona (Pamplona, Navarre, Spain) and a tunnel recreation of one of its gallery sections. The distribution of the considered experimental underground utility tunnels gallery is presented in Figure 3a with green lines on the map. As it can be seen from the figure, the tunnels have different sections, where the most remarkable characteristics are the different straight section lengths connected to one section with some intersections and slight curves. The dashed yellow line in the figure represents the straight tunnel section of 80-m-long where the experimental measurements were performed in order to calibrate the 3D-RL algorithm, as it is explained in Section 3. These tunnels were built in the early 2000s, under some specific streets, with the aim of easing maintenance operations. The tunnels consist of an outer layer of connected cement rings, with the inner part allocating several metallic trays for the corresponding utilities. The longest straight tunnel section is 327 m. A metallic hallway was also installed in order to provide access to the maintenance personnel, as can be seen in Figure 3b. Currently, they provide a route for the electrical connections, telephone connections, waste disposal and water distribution systems of the neighborhood. It is worth noting that slightly different utility tunnels have been built in newer neighborhoods of Pamplona, as can be seen in Figure 3c. Specifically, this picture corresponds to a recently repaired tunnel located in the suburb of Lezkairu, which began being built in 2008, and it has a total length of tunnels of 5600 m, plus 2300 m of sub-tunnels.
Regarding security and safety, the utility tunnels include common intrusion detection and antifire systems, with their corresponding local alarms. However, one of the main problems of these types of tunnels is the lack of wireless coverage within the tunnel, which initially was planned to be provided by external infrastructures, such as municipal Wi-Fi or mobile networks (3G/4G). However, in reality, these technologies are unable to provide coverage to all tunnel sections, as signal attenuation can be very high in some sections due to the different tunnel depths. In this sense, one of the main concerns of municipal authorities has been this lack of wireless coverage within the tunnels, as coverage is needed in order to provide security and safety to the maintenance staff. Thus, a specific indoor coverage system is required, and therefore, a complete radio planning analysis before the deployment of the wireless system is needed in order to secure coverage and provide good performance at an affordable cost.
As stated before, in order to characterize the variations of received power within the proposed utility tunnel, simulations by means of a 3D-RL code were performed. Material parameters were considered; basically, cement walls, aluminum trays and PVC tubes of different dimensions. A section of the simulation scenario is depicted in Figure 4, which corresponds with the same section where the experimental measurements were performed (see yellow dashed line in Figure 3a). The dimensions of the tunnel section are the following: width 2 m, height 3 m and length 80 m. Simulation results can be extended in length; however, due to measurement limitations, they will be considered up to the previously stated distance. The complete topological and morphological characteristics of this scenario were taken into account, prior visual inspection and measurement of dimensions within the tunnels. A schematic view of the employed simulation scenario is shown in Figure 4.
One of the particular aspects of this scenario is the fact that the tunnel is not a conventional tunnel (i.e., free hallway with walls that can be smooth or can present a certain degree of surface roughness and, hence, diffuse scattering), due to the fact that there is a set of service trays placed inside the tunnel. This leads to new interactions with the electromagnetic waves launched by the source, with strong multipath component generation given the large ratio between the length of the tunnel in comparison with the width of the tunnels, tray dimension and inter-tray distance. This condition implies the need for including more precise topological information of the scenario in order to consider the interaction, specifically for launched rays within a RL algorithm, such as the one employed and described in next subsection.

3. Radio Channel Analysis

3.1. 3D-RL Calibration and Validation

The underground tunnel scenario has distinctive characteristics that convert it in a particular scenario for radio wave propagation purposes. These peculiar characteristics are that it is an area that is considerably long (80-m-long, in this specific case) and extremely narrow (2-m width and 3-m height). It can be found in the literature [36,37], the convergence analysis of the in-house RL algorithm to obtain the optimal input parameters for the best performance of the algorithm in typical indoor complex scenarios. However, for this specific environment, those parameters are not the optimal, because the RL algorithm diverges due to the extremely narrow width of the scenario. The divergence of the algorithm means that there are a lot of cuboids of the scenario that do not have any ray if the same simulation parameters are used for this specific scenario. Due to that, the 3D-RL algorithm needs to be calibrated according to the experimental measurements to validate the results accuracy with the optimal input parameters. Therefore, a convergence analysis of the algorithm was done for this underground tunnel environment in order to characterize the 3D-RL convergence as a function of multiple parameters, such as cuboid dimensions, horizontal angular resolution (Δφ), vertical angular resolution (Δθ) or the maximum number of reflections. Due to the specific characteristics of the utility tunnel environment (i.e., strong reflections due to the partial conductor material presence and large obstacle density), the maximum number of reflections and the angular resolution of the launching rays was analyzed. For that purpose, the different number of multipath reflections and different number of launching rays were considered. Measurement results were compared with the different simulation results to find which simulation results achieved the best trade-off between precision and simulation computational cost, thus obtaining a calibrated 3D-RL technique that could be used for a radio-planning analysis in this type of complex underground environment. The considered frequency of operation was 868 MHz, emulating the operation of a typical LoRaWAN or ZigBee wireless communication system. The different cases for the convergence analysis, which was considered for the simulations of the underground tunnel scenario, are shown in Table 1.
Once the simulation results for the different cases were obtained, the validation and the corresponding calibration was obtained by performing radio channel measurements inside the tunnel section. The scenario in which the measurements were performed is a section of the utility underground tunnels in the city of Pamplona, which is depicted in Figure 3b. Measurements were carried out with the aid of an Agilent CSA N1996A to generate the radiofrequency (RF) signal and an Agilent FieldFox N9912A portable spectrum analyzer, both from Keysight Technologies. The transmitter was connected to a 5-dBi monopole antenna, whereas the receiver was connected to small vertical 0-dBi monopole antenna. The considered frequency of the operation was 868 MHz, with a transmitted power of −10 dBm. The transmitter antenna was placed at the beginning of the tunnel, located at the center with a height of 2 m, whereas the receiver points were placed at different distances in front of the transmitter along the tunnel. Measurements were performed with no users crossing the tunnels or other potential time-varying scatterers within the scenario. Results, compared with the different cases performed (shown in Table 1) in the 3D-RL simulations from the implemented scenario, are depicted in Figure 5. It can be seen from the picture that Cases 1, 3, 4 and 5 did not have good agreement with the measurements after 50-m distance approximately. This is because Cases 3, 4 and 5 considered larger angular resolutions than the other cases (one degree and 0.5 degrees of angular resolution) in comparison with 0.2 degrees in the other cases. In Case 1, the angular resolution was high, but the number of reflections was the lowest. It was shown that Case 1 did not have good agreement with the measurements. It was observed that Case 2 and Case 6 (both with angular resolutions of 0.2 degrees but different numbers of reflections considered) both had a good match with real measurements.
Table 2 shows the mean error and standard deviation of the comparison of the different simulation cases with the real measurements. It is also shown the simulation computational time of each one of the considered cases. It can be seen that Case 2 and Case 6 have the lowest mean error, as stated before. In Cases 1, 3, 4 and 5, the mean error increases considerably. With regards to the simulation computational time, it is observed that Case 6 has an increase of 186,861 s more than Case 2, with a slight decrease in the mean error. So, it can be concluded that, for the underground tunnel scenario, the optimal simulation parameters are Case 2, achieving a good trade-off between accuracy results and simulation computational burden.

3.2. Simulation Results

Once the calibration and verification of the 3D-RL algorithm for this singular scenario was achieved, the input optimal parameters, in terms of angular resolution and number of reflections, were applied in the simulations of the underground utility scenario. Large and small-scale radio wave propagation analyses were performed with the aid of simulations in order to provide a general signal propagation understanding within this complex and singular environment. The tunnel cross-sections have an important impact on the signal attenuation in this type of environment, as it is stated in [38]. In this sense, typical tunnel cross-sections can be classified in rectangular or circular shapes, where the tunnel width and height need to be considered when calculating the radio signal attenuation. Nevertheless, it has been shown in the literature [39,40] that cross-section dimensions within tunnels (i.e., different heights and widths of the tunnels) have a higher impact on the signal attenuation than distinct shapes, such as rectangular or circular. Therefore, the methodology, signal propagation analysis and WSN deployment presented in this work could potentially also be applied to cylinder tunnels, as long as tunnel cross-sections are of similar dimensions than the ones presented in this work. In fact, the cylindrical shape of the tunnel will generate slightly lower multipath propagation within the environment, as the number of reflections and diffractions within the surface of the tunnel will be lower.
Large scale propagation analysis was performed by means of the RF power distribution for the whole volume of the scenario. For that purpose, the transmitter antenna was placed at the beginning of the tunnel model, at a height of 2 m (represented by a red triangle in Figure 6). The same simulation parameters as in the experimental measurements were considered for the operational frequency and the transmitted power: 868 MHz and −10 dBm, respectively, emulating a ZigBee or LoRaWAN system at this specific frequency. Figure 6 shows the obtained RF power distribution estimations for the XY bidimensional plane at a height of 2 m (the same height as the transmitter antenna). It can be seen, the variability of received power along the distance, showing some spatial points with lower received power in the final part of the tunnel. The service trays on the sides of the tunnel generated a nonhomogeneous variation of the signal distribution levels along the tunnel distance.
To gain insight into the multipath behavior of the considered scenario, due to the rough surfaces of the walls and the service trays located in the tunnel, the power delay profile (PDP) is shown in Figure 7 for a specific point of the scenario, showing the high multipath trajectories within a time span of 200 and 2500 ns. The high multipath propagation within the environment can lead to signal fading in some spatial locations, which can lead to areas with a lack of coverage within the tunnel if the network design is not adequately performed. In addition, another parameter that can grossly quantify the multipath channel is the root mean square delay spread (RMS DS), referred also as σ τ , which shows the effects of dispersion and is determined from the PDP. The RMS DS for the considered utility tunnel environment is depicted in Figure 8. It is the square root of the second central moment of the PDP, defined as [41]:
σ τ = τ 2 ¯ ( τ ¯ ) 2
where
τ 2 ¯ = k a k   2 τ k 2 k a k 2 = k P ( τ k ) τ k 2 k P ( τ k )
and
τ ¯ = k a k 2   τ k k a k 2 =   k P (   τ k ) τ k k P ( τ k )
It must be pointed out that values for τ ¯ , τ 2 ¯ and σ τ depend on the choice of the noise threshold to process the received power at each delay P ( τ ) . The noise threshold used to depict the RMS DS shown in Figure 8 was −120 dBm, chosen to emulate the typical sensitivity of a ZigBee or LoRaWAN WSN that can be deployed in the tunnel environment.
It is observed that the RMS DS is higher in some areas closest to the center and the final part of the tunnel. This is because of the distribution of the service trays and tubes in the utility tunnel. Each spatial sample of the RMS DS is associated with a PDP. In this specific case, a multipath propagation contribution is higher in the central and final part of the tunnel, but it is necessary to consider all the spatial samples of the environment due to the influence of the topology and morphology of the specific scenario under consideration. Therefore, all the obstacles within the tunnel and the metallic structure of the service trays play an important role in the multipath propagation.
In order to have knowledge of the frequency domain channel description, the coherence bandwidth (CB) was obtained from the RMS DS for the considered scenario. The CB can be defined as the bandwidth over which the frequency correlation function is above 0.9; then, the CB is approximately B c 1 / 50 σ τ ; see [41].
Figure 9 shows the CB for the considered utility tunnel showing values that vary from 0.01 to 0.6 MHz. It can be seen that the lowest values of RMS DS cause higher values of CB. Therefore, a receiver in those areas can encounter higher multipath interference.

4. WSN Deployment

In this section, a generic WSN deployment case for complex heterogeneous utility tunnels is presented considering environmental and propagation conditions by means of a calibrated 3D-RL simulation tool. A complete WSN deployment feasibility study is introduced focusing on the reliability of wireless communications technologies (Zigbee and LoRaWAN) within the considered scenario and, particularly, in terms of security management, assessing the human presence propagation impact over the system. Finally, an optimized WSN deployment is presented considering the morphology and topology of the real case study underground tunnel system environment and an infrastructure and work safety monitoring application prototype.

4.1. Safety Management and Occupational Risk Control

One of the main problems of working in tunnels and service galleries, which is aggravated by their length and narrowness, is the safety management and the occupational risks control. The health of the working environment of the maintenance staff must always be guaranteed. This includes monitoring temperature and humidity levels but, also, monitoring CO and CO2 levels; the presence of flammable or toxic gases and compounds and the presence of animals (rats, snakes, etc.), among others. In addition, mobility within the tunnel, the absence of obstacles or elements that hinder movement must be guaranteed. It should be borne in mind that employees, on many occasions, must move tools to carry out their maintenance, repair and installation tasks.
Current regulations related to labor risk prevention state that work must be done at least in pairs and desirable in work teams to minimize risks and prevent accidents. In addition, workers must notify the coordination center of the task to be performed; the starting time; the estimated duration; the working conditions (lighting, power supply cut, etc.); the possible risks and any other information that may help to ensure safety at work during the job performance. Hence, personnel monitoring implies that the staff will carry their own sensors (located in the inhalation zone and at head height) and must be able to interconnect to the monitoring network. Thus, the use of a wireless network is particularly advantageous in the latter case, as it offers a greater ease of work (hands-free) and greater mobility for workers. However, the presence of human beings within the tunnels can affect the radio propagation, even more considering their narrowness. Therefore, an adequate RF characterization and a complete risk control assessment must be performed in order to warrant compliance with occupational safety management within the utility tunnels.

4.2. Technology Analysis: ZigBee 868 MHz vs. LoRaWAN Range within the Tunnels

3D-RL simulations were shown to be accurate for an optimized radio-planning analysis in this complex and particular underground environment. Based on these preliminary results, the coverage for both Zigbee and LoRaWAN wireless devices can be estimated (i.e., the number of wireless nodes required for the whole tunnel) considering their corresponding device sensitivities, which are around −100 dBm and −148 dBm, respectively. In order to provide clear insight, range estimations within the tunnel for both wireless technologies were obtained. For that purpose, a transmitter was placed on the ceiling at the beginning of the tunnel with a length equal to the longest tunnel section of the scenario under analysis (327 m, as it is presented in Figure 3a). The transmitted powers were set to 13 dBm and 14 dBm for ZigBee and LoRaWAN, respectively. As can be concluded from the results depicted in Figure 10, both wireless technologies covered the entire length of the tunnel, although LoRaWAN presented better results due to its better sensitivity (−148 dBm) compared to ZigBee (−100 dBm). Thus, in terms of wireless link distance within a straight tunnel of these characteristics, both wireless technologies could be employed.
Nevertheless, in the real environment, different phenomena can lead to channel impairments, such as the presence of people within the tunnels or worsening the radio propagation conditions. Therefore, a further study is required in order to provide safe, robust, reliable and efficient WSN deployment.

4.3. Scatterers Analysis: Human Presence within the Tunnel

An important factor that has to be assessed in terms of radio propagation in small cross-section utility tunnels is the presence of human beings. It is a well-known fact that human bodies create a shadowing effect [42], and in such a narrow environment like the underground tunnels presented in this work, this effect was expected to be especially significant. In order to have insight into this phenomenon, different simulations were performed using the calibrated 3D-RL simulation tool. Human body models are included in the 3D-RL software, emulating occupational personnel at the location of interest. This human body model was implemented for the simulation, considering the dimensions of a standard person and all the dispersive material properties of its different body parts (i.e., skin, bones, muscles, blood, different organs, etc.). A further description of the considered human body model, as well as a signal attenuation analysis for different parts of the human body, can be found in reference [43]. By including computational human body models within the considered scenario, the results of the human presence effect were obtained, allowing analysis and comparison. In the context of the utility tunnel, teams of up to three workers are generally sent into the tunnels for performing different security and maintenance tasks, in order to minimize the risks and prevent accidents due to municipal authorities’ requirements regarding occupational risk prevention legislation. Therefore, a complete analysis of the human presence within the tunnel was carried out considering different densities of human body models in different locations.
First of all, the received power results of a specific scenario with no humans is presented in Figure 11, considering a transmitter located on the tunnel ceiling at x-axis = 20 m for reference and comparison purposes.
Then, a single human body model was included in the scenario at different x-axis locations in order to analyze the influence of the surrounding environment and the transmitter distance along the tunnel. The considered locations were the following: 20 m, 27 m, 34 m, 41 m and 48 m. In Figure 12, the corresponding RF power levels distribution results at 1.5-m height are presented. As expected, the major RF distribution changes are mainly in the area behind the location of the human body, due to the shielding effect.
In order to gain more insights into these results, Figure 13 shows the differences between the case without human presence (i.e., Figure 11) and the different case studies considering a single human body model at distinct locations (i.e., Figure 12). As can be clearly observed from Figure 13, the presence of a single worker within the tunnel in different locations leads to RF power level differences up to 18 dB. Although this difference is lower the higher the distance from the transmitter is, it has to be taken into account for an accurate range estimation. Thus, a margin of about 20 extra dBs should be added to the wireless communication range estimated in the previous subsection both for ZigBee and LoRaWAN technologies.
Finally, extra simulations were carried out in order to assess if the obtained 20 dBs are enough when work teams of three people are within the utility tunnel. In Figure 14, a rendered view of the inclusion of the three human body models within the created tunnel scenario for the 3D-RL simulation is presented.
Figure 15 presents the obtained estimation results in terms of the RF power level distribution differences between the case without a human presence and the case study with a three-worker team at different locations along the tunnel. From the results, a mild growth variation can be observed when comparing the three-worker team case study with the single-worker case study, and therefore, it can be stated that the previous estimated 20-dB margin is enough for a proper RF planification when considering teams of up to three workers within the tunnel.
In general, the obtained results regarding the study of the human presence within the tunnel show that a margin of about approximately 20 extra dBs is required in order to preserve the received RF power level above the sensitivity level, thus ensuring the wireless communication links between nodes, even when up to a three-worker team is considered within the tunnel. This extra margin corresponds to no need of extra range reduction estimated in the previous subsection analysis, as can be seen in Figure 16. However, since the sensitivity value of ZigBee is −100 dBm, in this case, a reduction of the distance between nodes is desirable in order to assure reliable wireless communication at every circumstance. This restriction was applied for the design of the WSN presented in the following subsection.

4.4. Proposed WSN Deployment

In this section, an optimized WSN deployment is proposed considering the morphology and topology of the real underground utility tunnel system. In general, the network topology gains a great importance for an adequate deployment and infrastructure design in this type of complex underground environments, since their particular morphology (for instance, tunnels in the presented case) could limit the use of the star topology (the typical one for LoRaWAN). From the previous results, it can be concluded that LoRaWAN technology presents a much higher range in terms of distance and better robustness with human presence due to its better sensitivity. Besides, it provides a high tolerance for interference, which is highly relevant in a confined scenario with the presence of tubes, grids and metal trays. However, considering the wireless system technologies’ capabilities and constraints, and the real utility tunnel system conditions in terms of connections, intersections and infrastructure location possibilities, there are several key factors that lead to the election of ZigBee technology to be deployed within the considered scenario. The first LoRaWAN limitation is the reduced transfer capacity (up to 242 bytes), but the most critical one is the network topology and its deployment in underground conditions. While the mesh topology provided by ZigBee gives the option of deploying one single gateway, the star topology of LoRaWAN leads to the deployment of several gateways due to the morphology of the tunnels. Furthermore, taking into account that the considered tunnels lack of communication systems infrastructure and coverage of external systems (such as the mobile network), the deployment of gateways presents serious difficulties, since they need external communication links in order to send the data of the underground tunnel network to the cloud. Therefore, the possibility of deploying a unique gateway is very advantageous.
For the proposed WSN deployment, a Digi International’s XBee Zigbee Gateway was selected, allowing GPRS communications. The gateway was designed to be installed in the center of the tunnel system for network performance optimization and located under a sewer, with the aim of facilitating the wireless communication with the mobile network. An efficient node of deployment was implemented with specific nodes located in key corners and intersections, considering the particular morphology of the tunnels and the maximum distance between the nodes obtained in the previous subsection analysis. Digi’s XBee SX 868 RF Modules were employed for the sensor nodes. A schematic view of the WSN deployment is depicted in Figure 17 representing the location of the ZigBee nodes and the gateway over the real utility tunnel system map.

4.5. Infrastructure and Work Safety Monitoring Application

In order to guarantee the security and protection of the infrastructure against unauthorized access, vandalism or sabotage, as well as to ensure occupational safety at work, a surveillance and monitoring system prototype was developed to allow infrastructure and work safety remote monitoring. In general, these types of systems usually include WSN, closed-circuit television (CCTV) and an antifire system, each of them with its own independent communication network. Although they could share the same communication network, for security and fault tolerance reasons, normally, three different communication networks are deployed, both wired and wireless. Adequate precautions, such as the use of shielded cables or fiber optics, should be incorporated when installing the complete wiring system to ensure the immunity to the effects of electromagnetic interference. However, sometimes, due to cost or deployment reasons, a wired network is not available, requiring a wireless communication network. In the analyzed real case, a fiber-optic network was available for the city’s security cameras, although no cameras were included inside the underground utility tunnel systems. If they were included, they could be connected to the fiber-optic network without further complications. Nevertheless, either special infrared cameras or visible spectrum cameras, in conjunction with presence sensors that activate the lighting system, should be included as well. This second option is more efficient in terms of infrastructure installation costs due to the difference in price of the cameras, but it has a higher probability of failure because it involves more elements (sensors and actuators). While, currently, no CCTV cameras are available in the galleries of the underground utility tunnel system, a wireless infrastructure and work safety monitoring application is proposed.
Regarding infrastructure security and work safety monitoring, WSN are devoted to monitoring possible leaks, zones or areas and personnel. When monitoring leaks (mainly water and gases), sensors should be installed as close as possible of the possible leakage in order to facilitate that the detection of the leak, if any, be as quick as possible. When monitoring zones/areas, detectors are regularly distributed throughout the area according to the potential risk (and its assessment) of leakage. On the one hand, some combustible gases, such as hydrogen or methane, are lighter than air, so the detector should be placed above the possible leak origin point or on the ceiling, where there is a high probability of accumulations of these gases. On the other hand, other gases such as propane or butane tend to move through the lower parts and can cover large areas. Therefore, in order to detect this type of gas in the fastest possible way, it is necessary to place the detectors in the lower parts of the installation, avoiding the contact with water. Finally, personnel monitoring is pivotal in order to provide work safety conditions, especially in underground environments. By means of a wireless network, surveillance and maintenance teamwork can be controlled and monitored along the galleries using wearables or worn body sensors, allowing free mobility and work service.
The development of an infrastructure-monitoring tool implies an adequate analysis and design of the whole system but, also, an adequate choice of the types of sensors, their locations, their configurations and their communication frequencies. In addition, the tool must allow the definition of alarm levels and the processes and procedures to be activated in case of alarm confirmation. Infrastructure monitoring will focus on both temperature and humidity control to prevent mold and bacteria buildup and to control potential sources of contamination. To do this, the humidity must be kept below 50% (desired values 30–50%). The same applies to dangerously high levels of carbon monoxide or dioxide in the air or nitrogen compounds (monoxide, dioxide and other nitrogen oxides) to prevent poisoning of the maintenance staff.
Considering the real case utility tunnel system and the proposed WSN deployment presented in the previous subsection, an infrastructure and work safety monitoring application prototype was developed in order to collect the data obtained by the sensors located in the gallery in real time and store them in a database allowing cloud computing. A schematic view of the monitoring tool architecture is depicted in Figure 18. The prototype was developed using open-source tools such as the Bootstrap framework [44] for the development of the web application Grafana [45] for the analytics and interactive visualization and the D3.js [46] JavaScript library for the production of dynamic and interactive data visualizations in web browsers. The system also includes a web service in charge of data collection, developed in Java, and a relational database (MySQL) [47].
The design premises are security and use simplicity. The monitoring tool is responsive, which allows its use on mobile devices (smartphones and tablets), workstations or video walls. It embeds Grafana’s panels using iframes. On the one hand, there is a Grafana dashboard for the direct monitoring of the gallery’s health parameters, which allows the visualization of the data provided by the tunnel’s sensor network, which is depicted in Figure 19. On the other hand, the monitoring tool (see Figure 20), in addition to integrating the information provided by Grafana, includes the monitoring of the activities performed inside the gallery. These activities include dead man functionalities, monitoring of the task progress, access and exit points of the gallery, evacuation protocol and instructions in case of emergency, emergency acoustic and light signals activation, etc. The tool includes tree-based and sunburst visualization to describe the workflows.
The data collected by the WSN is stored, processed and analyzed in a private cloud. For that purpose, an Apache CloudStack was implemented over hypervisors VMware ESXI managed by a VMware vSphere. The resulting infrastructure allows the agile and easy deployment and management of the services. Besides, it allows computing orchestration, user and account management and provides data storage and data analytics. The proposed computing cloud is easily scalable and allows hot migrations according to the availability of the servers.

5. Conclusions

In this work, the radio wave propagation assessment and wireless system deployment for real urban utility tunnels was described. The impact of large obstacle density, given by the inclusion of elements such as service trays or handrails, as well as the presence of occupational staff within the underground utility tunnels, was analyzed by means of an in-house 3D-RL code. The calibration and validation of the 3D-RL algorithm for this singular scenario were presented by comparison with experimental measurements in the real scenario, showing good agreement between them. Once the optimal simulation parameters were obtained from the experimental measurements, the large and small-scale simulation results were presented to characterize the radio wave propagation within the tunnel. A complete volumetric analysis was performed in order to consider the location of wireless transceivers and enable them to be embedded within the elements of the tunnels. The system-level analysis considering the conventional transceiver parameters (ZigBee and LoRaWAN) was performed for an 868-MHz operational frequency, obtaining the initial deployment planning results. The use of higher frequencies could lead to a higher signal attenuation, and then, the coverage could be affected within the underground tunnel. Additionally, an infrastructure security and work safety monitoring application were developed in which the information was stored in a private Cloud. The proposed methodology enables to consider specific propagation losses in complex urban utility tunnels in order to deploy wireless communication systems with optimal coverage, minimum interference and an enhanced quality of service.

Author Contributions

Conceptualization, L.A. and P.L.-I.; methodology, L.A.; software, M.C.-E. and I.P.; validation, M.C.-E.; formal analysis, M.C.-E., L.A. and P.L.-I.; investigation, M.C.-E. and L.A.; data curation, M.C.-E. and I.P.; writing—original draft preparation, L.A. and M.C.-E.; writing—review and editing, P.L.-I., J.J.A. and F.F.; visualization, E.A. and M.C.-E.; supervision, F.F.; wireless network planning, P.L.-I. and I.P.; cloud design and development, J.J.A., monitoring system design and development J.J.A. and J.V. and funding acquisition, F.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the School of Engineering and Sciences at Tecnologico de Monterrey, in part by Project RTI2018-095499-B-C31, funded by Ministerio de Ciencia, Innovación y Universidades, Gobierno de España (MCIU/AEI/FEDER, UE) and, in part, by the European Union’s Horizon 2020 research and Innovation program under grant agreement N°774094 (Stardust-Holistic and Integrated Urban Model for Smart Cities).

Acknowledgments

The authors would like to acknowledge the support and collaboration of the School of Engineering and Sciences and the Focus Research Group of Telecommunications and Networks at Tecnologico de Monterrey.

Conflicts of Interest

The authors declare no conflict of interest. The statements made herein are solely the responsibility of the authors.

References

  1. Mahmoud, S.F. Wireless Transmission in Tunnels with Non-Circular Cross Section. IEEE Trans. Antennas Propag. 2010, 58, 613–616. [Google Scholar] [CrossRef]
  2. Hrovat, A.; Kandus, G.; Javornik, T. A Survey of Radio Propagation Modeling for Tunnels. IEEE Commun. Surv. Tutor. 2014, 16, 658–669. [Google Scholar] [CrossRef]
  3. Wu, Y.; Zheng, G.; Saleem, A.; Zhang, Y.P. An Experimental Study of MIMO Performance Using Leaky Coaxial Cables in a Tunnel. IEEE Antennas Wirel. Propag. Lett. 2017, 16, 1663–1666. [Google Scholar] [CrossRef]
  4. Forooshani, A.E.; Bashir, S.; Michelson, D.G.; Noghanian, S. A Survey of Wireless Communications and Propagation Modeling in Underground Mines. IEEE Commun. Surv. Tutor. 2013, 15, 1524–1545. [Google Scholar] [CrossRef]
  5. Moridi, M.A.; Sharifzadeh, M.; Kawamura, Y.; Jang, H.D. Development of wireless sensor networks for underground communication and monitoring systems (the cases of underground mine environments). Tunn. Undergr. Space Technol. 2018, 73, 127–138. [Google Scholar] [CrossRef] [Green Version]
  6. Athanasiadou, G.E.; Tsoulos, G.V.; Zarbouti, D.A.; Valavanis, I.K. Optimizing Radio Network Planning Evolution towards Microcellular Systems. Wirel. Pers. Commun. 2019, 106, 521–534. [Google Scholar] [CrossRef]
  7. Xu, B.; Li, J.; Yang, Y.; Postolache, O.; Wu, H. Robust modeling and planning of radio-frequency identification network in logistics under uncertainties. Int. J. Distrib. Sens. Netw. 2018, 14, 1550147718769781. [Google Scholar] [CrossRef]
  8. Athanasiadou, G.E.; Fytampanis, P.; Zarbouti, D.A.; Tsoulos, G.V.; Gkonis, P.K.; Kaklamani, D.I. Radio Network Planning towards 5G mmWave Standalone Small-Cell Architectures. Electronics 2020, 9, 339. [Google Scholar] [CrossRef] [Green Version]
  9. González González, D.; Mutafungwa, E.; Haile, B.; Hämäläinen, I.; Poveda, H. A Planning and Optimization Framework for Ultra Dense Cellular Deployments. Mobile Inf. Syst. 2017, 2017, 9242058. [Google Scholar] [CrossRef]
  10. Hoppe, R.; Wölfle, G.; Jakobus, U. Wave Propagation and Radio Network Planning Software WinProp added to the Electromagnetic Solver Package FEKO. In Proceedings of the 2017 International Applied Computational Electromagnetics Society Symposium—Italy (ACES), Firenze, Italy, 26–30 March 2017. [Google Scholar]
  11. Araujo, W.; Fogarolli, R.; Seruffo, M.; Cardoso, D. Deployment of small cells and a transport infrastructure concurrently for next-generation mobile access networks. PLoS ONE 2018, 13, e0207330. [Google Scholar] [CrossRef]
  12. Iqbal, M.; Naeem, M.; Anpalagan, A.; Ahmed, A.; Azam, M. Wireless Sensor Network Optimization: Multi-Objective Paradigm. Sensors 2015, 15, 17572–17620. [Google Scholar] [CrossRef] [PubMed]
  13. Sanchez-Iborra, R.; Sanchez-Gomez, J.; Ballesta-Viñas, J.; Cano, M.-D.; Skarmeta, A.F. Performance Evaluation of LoRa Considering Scenario Conditions. Sensors 2018, 18, 772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sulaiman BenSaleh, M.; Saida, R.; Hadj Kacem, Y.; Abid, M. Wireless Sensor Network Design Methodologies: A Survey. J. Sens. 2020. [Google Scholar] [CrossRef]
  15. Priyadarshi, R.; Gupta, B.; Anurag, A. Deployment techniques in wireless sensor networks: A survey, classification, challenges, and future research. J. Supercomput. 2020, 76, 7333–7373. [Google Scholar] [CrossRef]
  16. Gao, Y.; Wang, J.; Wu, W.; Sangaiah, A.K.; Lim, S.-J. Travel Route Planning with Optimal Coverage in Difficult Wireless Sensor Network Environment. Sensors 2019, 19, 1838. [Google Scholar] [CrossRef] [Green Version]
  17. Iliodromitis, A.; Pantazis, G.; Vescoukis, V. 2D wireless sensor network deployment based on Centroidal Voronoi Tessellation. AIP Conf. Proc. 2017, 1836, 020010. [Google Scholar]
  18. Kanaris, L.; Sergiou, C.; Kokkinis, A.; Pafitis, A.; Antoniou, N.; Stavrou, S. On the Realistic Radio and Network Planning of IoT Sensor Networks. Sensors 2019, 19, 3264. [Google Scholar] [CrossRef] [Green Version]
  19. Delogne, P. EM Propagation in Tunnels. IEEE Trans. Antennas Propag. 1991, 39, 401–406. [Google Scholar] [CrossRef]
  20. Mariage, P.; Lienard, M.; Degauque, P. Theoretical and Experimental Approach of the Propagation of High Frequency Waves in Road Tunnels. IEEE Trans. Antennas Propag. 1994, 42, 75–81. [Google Scholar] [CrossRef]
  21. Chen, S.; Jeng, S. SBR Image Approach for Radiowave Propagation in Tunnels with and without Traffic. IEEE Trans. Antennas Propag. 1996, 45, 570–578. [Google Scholar]
  22. Celaya-Echarri, M.; Azpilicueta, L.; López-Iturri, P.; Picallo, I.; Aguirre, E.; Astrain, J.J.; Villadangos, J.; Falcone, F. Wireless Channel Characterization and System Analysis of Complex Utility Tunnel Environments. In Proceedings of the 6th International Electronic Conference on Sensors and Applications, 15–30 November 2019; Available online: https://www.mdpi.com/2504-3900/42/1/53 (accessed on 23 November 2020).
  23. Sun, Z.; Akyildiz, I. Channel modeling and analysis for wireless networks in underground mines and road tunnels. IEEE Trans. Commun. 2010, 58, 1758–1768. [Google Scholar] [CrossRef]
  24. Martelly, R.; Janaswamy, R. An ADI-PE Approach for Modeling Radio Transmission Loss in Tunnels. IEEE Trans. Antennas Propag. 2009, 57, 1759–1770. [Google Scholar] [CrossRef]
  25. Sakai, K.; Koshiba, M. Analysis of electromagnetic field distribution in tunnels by the boundary-element method. IEE Proc. 1990, 137, 202–208. [Google Scholar] [CrossRef]
  26. Bernardi, P.; Caratelli, D.; Cicchetti, R.; Schena, V.; Testa, O. A Numerical Scheme for the Solution of the Vector Parabolic Equation Governing the Radio Wave Propagation in Straight and Curved Rectangular Tunnels. IEEE Trans. Antennas Propag. 2009, 57, 3249–3257. [Google Scholar] [CrossRef]
  27. Didascalou, D.; Schäfer, T.M.; Weinmann, F.; Wiesbeck, W. Ray-Density Normalization for Ray-Optical Wave Propagation Modeling in Arbitrarily Shaped Tunnels. IEEE Trans. Antennas Propag. 2000, 48, 1316–1325. [Google Scholar] [CrossRef]
  28. Zhang, Y.P.; Hwang, Y. Enhancement of rectangular tunnel waveguide model. In Proceedings of the Asia Pacific Microwave Conference, Hong Kong, China, 2–5 December 1997. [Google Scholar]
  29. Mohammed, F.H.; Samir, F.M. A High-Order Compact-FDTD Algorithm for Electrically Large Waveguide Analysis. IEEE Trans. Antennas Propag. 2008, 56, 2589–2598. [Google Scholar]
  30. Rana, M.; Sanagavarapu, A. Segmented-Locally-One-Dimensional-FDTD Method for EM Propagation Inside Large Complex Tunnel Environments. IEEE Trans. Magn. 2012, 48, 223–226. [Google Scholar] [CrossRef]
  31. Lopez-Iturri, P.; Aguirre, E.; Azpilicueta, L.; Astrain, J.J.; Villadangos, J.; Falcone, F. Challenges in Wireless System Integration as Enablers for Indoor Context Aware Environments. Sensors 2017, 17, 1616. [Google Scholar] [CrossRef] [Green Version]
  32. Azpilicueta, L.; Astrain, J.J.; Lopez-Iturri, P.; Granda, F.; Vargas-Rosales, C.; Villadangos, J.; Perallos, A.; Bahillo, A.; Falcone, F. Optimization and Design of Wireless Systems for the Implementation of Context Aware Scenarios in Railway Passenger Vehicles. IEEE Trans. Intell. Transp. Syst. 2017, 18, 2838–2850. [Google Scholar] [CrossRef]
  33. López-Iturri, P.; Aguirre, E.; Azpilicueta, L.; Astrain, J.J.; Villadangos, J.; Falcone, F. Implementation and Analysis of ISM 2.4GHz Wireless Sensor Network System in Judo Training Venues. Sensors 2016, 16, 1247. [Google Scholar] [CrossRef] [Green Version]
  34. Córdoba, A.; Astrain, J.J.; Villadangos, J.; Azpilicueta, L.; López-Iturri, P.; Aguirre, E.; Falcone, F. SesToCross: Semantic Expert System to Manage Single-Lane Road Crossing. IEEE Trans. Intell. Transp. Syst. 2017, 18, 1221–1233. [Google Scholar] [CrossRef]
  35. Azpilicueta, L.; López-Iturri, P.; Aguirre, E.; Martínez, C.; Astráin, J.J.; Villadangos, J.; Falcone, F. Evaluation of Deployment Challenges of Wireless Sensor Networks at Signalized Intersections. Sensors 2016, 16, 1140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Azpilicueta, L.; Rawat, M.; Rawat, K.; Ghannouchi, F.; Falcone, F. Convergence Analysis in Deterministic 3D Ray Launching Radio Channel Estimation in Complex Environments. ACES J. 2014, 29, 256–271. [Google Scholar]
  37. Azpilicueta, L.; López-Iturri, P.; Aguirre, E.; Falcone, F. An accurate UTD Extension to a Ray Launching Algorithm for the Analysis of Complex Indoor Radio Environments. J. Electromagn. Waves Appl. 2016, 30, 43–60. [Google Scholar] [CrossRef]
  38. Changsen, Z.; Yan, M. Effects of cross section of mine tunnel on the propagation characteristics of UHF radio wave. In Proceedings of the 7th International Symposium on Antennas, Propagation EM Theory ISAPE ’06, Guilin, China, 26–29 October 2006; pp. 1–5. [Google Scholar]
  39. Wang, S. Radio wave attenuation character in the confined environments of rectangular mine tunnel. Mod. Appl. Sci. 2010, 4, 65–70. [Google Scholar] [CrossRef] [Green Version]
  40. Cheng, L.; Zhang, P. Influence of dimension change on radio wave propagation in rectangular tunnels. In Proceedings of the 5th International Conference on Wireless Communications, Networking and Mobile Computing, Beijing, China, 24–26 September 2009; pp. 1–3. [Google Scholar]
  41. Rappaport, T. Wireless Communications: Principles and Practice, 2nd ed.; Prentice Hall PTR: Upper Saddle River, NJ, USA, 2002. [Google Scholar]
  42. Obayashi, S.; Zander, J. A body-shadowing model for indoor radio communication environments. IEEE Trans. Antennas Propag. 1998, 46, 920–927. [Google Scholar] [CrossRef]
  43. Aguirre, E.; Arpon, J.; Azpilicueta, L.; de Miguel Bilbao, S.; Ramos, V.; Falcone, F.J. Evaluation of electromagnetic dosimetry of wireless systems in complex indoor scenarios with human body interaction. Prog. Electromagn. Res. B 2012, 43, 189–209. [Google Scholar] [CrossRef] [Green Version]
  44. Bootstrap. Available online: https://getbootstrap.com/ (accessed on 30 July 2020).
  45. Grafana. Available online: https://grafana.com/ (accessed on 30 July 2020).
  46. D3 Data-Driven Documents. Available online: https://d3js.org/ (accessed on 30 July 2020).
  47. MySQL. Available online: https://www.mysql.com/ (accessed on 30 July 2020).
Figure 1. Schematic view of the structure of this work. WSN: wireless sensor network and 3D-RL: three-dimensional ray-launching.
Figure 1. Schematic view of the structure of this work. WSN: wireless sensor network and 3D-RL: three-dimensional ray-launching.
Sensors 20 06710 g001
Figure 2. Principle of operation of the in-house 3D-RL algorithm.
Figure 2. Principle of operation of the in-house 3D-RL algorithm.
Sensors 20 06710 g002
Figure 3. Underground utility tunnels under analysis. (a) Tunnel location on the map of the Old Town of Pamplona. (b) Picture of the tunnel where the measurements were taken. (c) Picture of a newer tunnel built in the neighborhood of Lezkairu in Pamplona.
Figure 3. Underground utility tunnels under analysis. (a) Tunnel location on the map of the Old Town of Pamplona. (b) Picture of the tunnel where the measurements were taken. (c) Picture of a newer tunnel built in the neighborhood of Lezkairu in Pamplona.
Sensors 20 06710 g003
Figure 4. Schematic representation of the simulation scenario implemented to analyze the radio propagation losses within the underground utility tunnels.
Figure 4. Schematic representation of the simulation scenario implemented to analyze the radio propagation losses within the underground utility tunnels.
Sensors 20 06710 g004
Figure 5. Comparison of different cases of simulation and measurement results for the service tunnel section under analysis.
Figure 5. Comparison of different cases of simulation and measurement results for the service tunnel section under analysis.
Sensors 20 06710 g005
Figure 6. Radiofrequency (RF) power distribution for the XY bidimensional plane at 2-m height for 868 MHz.
Figure 6. Radiofrequency (RF) power distribution for the XY bidimensional plane at 2-m height for 868 MHz.
Sensors 20 06710 g006
Figure 7. Power delay profile at a given cuboid located at the center (60 m, 1 m and 2 m) at 868-MHz frequency.
Figure 7. Power delay profile at a given cuboid located at the center (60 m, 1 m and 2 m) at 868-MHz frequency.
Sensors 20 06710 g007
Figure 8. Root mean square delay spread (RMS DS) estimation within the utility tunnel for 2-m height for 868-MHz frequency.
Figure 8. Root mean square delay spread (RMS DS) estimation within the utility tunnel for 2-m height for 868-MHz frequency.
Sensors 20 06710 g008
Figure 9. Coherence bandwidth within the considered utility tunnel for a 2-m height for 868-MHz frequency.
Figure 9. Coherence bandwidth within the considered utility tunnel for a 2-m height for 868-MHz frequency.
Sensors 20 06710 g009
Figure 10. Bidimensional plane of the RF power level distribution at a height of 3 m for (a) ZigBee 868 MHz and (b) LoRaWAN.
Figure 10. Bidimensional plane of the RF power level distribution at a height of 3 m for (a) ZigBee 868 MHz and (b) LoRaWAN.
Sensors 20 06710 g010
Figure 11. Bidimensional plane of the RF power level distribution at a 1.5-m height for the case without human presence.
Figure 11. Bidimensional plane of the RF power level distribution at a 1.5-m height for the case without human presence.
Sensors 20 06710 g011
Figure 12. Bidimensional planes of the RF power level distribution at a 1.5-m height for different human being locations along the tunnel.
Figure 12. Bidimensional planes of the RF power level distribution at a 1.5-m height for different human being locations along the tunnel.
Sensors 20 06710 g012
Figure 13. RF power level distribution differences between the case without a human presence and the cases with human beings at different locations along the tunnel.
Figure 13. RF power level distribution differences between the case without a human presence and the cases with human beings at different locations along the tunnel.
Sensors 20 06710 g013
Figure 14. Inclusion of 3 computational human body models in the considered scenario for the 3D-RL simulation, emulating a maintenance team working in the utility tunnels.
Figure 14. Inclusion of 3 computational human body models in the considered scenario for the 3D-RL simulation, emulating a maintenance team working in the utility tunnels.
Sensors 20 06710 g014
Figure 15. RF power level distribution difference between the case without people and the cases with a work team of 3 people at different locations along the tunnel.
Figure 15. RF power level distribution difference between the case without people and the cases with a work team of 3 people at different locations along the tunnel.
Sensors 20 06710 g015
Figure 16. Bidimensional plane at a 3-m height representation of the RF power level distribution after applying the 20-dB margin due to the human presence for (a) ZigBee 868 MHz and (b) LoRaWAN.
Figure 16. Bidimensional plane at a 3-m height representation of the RF power level distribution after applying the 20-dB margin due to the human presence for (a) ZigBee 868 MHz and (b) LoRaWAN.
Sensors 20 06710 g016
Figure 17. Schematic view of the proposed ZigBee 868-MHz WSN deployment.
Figure 17. Schematic view of the proposed ZigBee 868-MHz WSN deployment.
Sensors 20 06710 g017
Figure 18. Architecture of the monitoring tool.
Figure 18. Architecture of the monitoring tool.
Sensors 20 06710 g018
Figure 19. Time series monitoring by means of Grafana (iFrames).
Figure 19. Time series monitoring by means of Grafana (iFrames).
Sensors 20 06710 g019
Figure 20. Web-based monitoring tool.
Figure 20. Web-based monitoring tool.
Sensors 20 06710 g020
Table 1. Different cases for the convergence analysis of the algorithm.
Table 1. Different cases for the convergence analysis of the algorithm.
Δφ (rad)Δθ (rad)Cuboids Size (m)Number of ReflectionsNumber of Rays
Case 1 π / 900 π / 900 0.15810,000
Case 2 π / 900 π / 900 0.17810,000
Case 3 π / 360 π / 360 0.17129,600
Case 4 π / 180 π / 180 0.1732,400
Case 5 π / 180 π / 180 0.11332,400
Case 6 π / 900 π / 900 0.113810,000
Table 2. Comparison between different simulation cases.
Table 2. Comparison between different simulation cases.
Case 1Case 2Case 3Case 4Case 5Case 6
Mean Error (dB)2.951.342.745.245.241.06
Standard Deviation (dB)16.145.9016.5614.2914.296.14
Computational Time (s)52,66478,45948528831656265,320
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Celaya-Echarri, M.; Azpilicueta, L.; Lopez-Iturri, P.; Picallo, I.; Aguirre, E.; Astrain, J.J.; Villadangos, J.; Falcone, F. Radio Wave Propagation and WSN Deployment in Complex Utility Tunnel Environments. Sensors 2020, 20, 6710. https://doi.org/10.3390/s20236710

AMA Style

Celaya-Echarri M, Azpilicueta L, Lopez-Iturri P, Picallo I, Aguirre E, Astrain JJ, Villadangos J, Falcone F. Radio Wave Propagation and WSN Deployment in Complex Utility Tunnel Environments. Sensors. 2020; 20(23):6710. https://doi.org/10.3390/s20236710

Chicago/Turabian Style

Celaya-Echarri, Mikel, Leyre Azpilicueta, Peio Lopez-Iturri, Imanol Picallo, Erik Aguirre, Jose Javier Astrain, Jesús Villadangos, and Francisco Falcone. 2020. "Radio Wave Propagation and WSN Deployment in Complex Utility Tunnel Environments" Sensors 20, no. 23: 6710. https://doi.org/10.3390/s20236710

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