Investigation of MIMO Channel Characteristics in a Two-Section Tunnel at 1 . 4725 GHz

Copyright © 2017 Rongchen Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper presents results from a wide band single-input–single-output (SISO) and 16 × 16 virtual multiple-input–multiple-output (MIMO)measurement campaign at a center frequency of 1.4725GHz in a 100-meter long tunnel laboratorywhich is terminated by a vertical wall with a metallic door.The path loss, root-mean-square delay spread (RMS-DS) characteristics, and power delay profiles (PDPs) are described. In addition, we provide results for the MIMO channel amplitude matrix, which offers a new perspective in understanding MIMO characteristics in tunnel scenarios. Our measurement results are analyzed and compared to ray tracing simulations.The relationships among the angle spread, channelmatrix singular values, andMIMOcapacity at various link distances are illustrated, and these provide insights into MIMO system deployment.


Introduction
To ease the growing population pressure on urban traffic, major cities in China have accelerated subway construction.It is expected that by 2020 more than 45 cities in China will have subways, and all new subways under construction must have a large rush hour capacity.Increasing train running speed and reducing the interval between adjacent trains enable a passenger capacity increase, but at the same time this introduces some security risks such as those due to train malfunction or other disasters (e.g., fires).Thus, these trains require various levels of monitoring to ensure safe and efficient operation.For example, video information should be obtained in real time at various locations.
A wide band and MIMO wireless communication system can satisfy the increasing data transmission requirements for both train control systems and passenger communications.Before any MIMO system is deployed underground, the MIMO channel should be accurately characterized.Much work in the literature has discussed MIMO link performance by measurement and various theoretical methods.
Extensive work has been done on propagation and MIMO channel capacity based on measurements in subway tunnels.Using electromagnetic field theory, for example, [1][2][3][4][5][6], researchers have confirmed that strong wave guiding effects make the tunnel channel characteristics different from those in other scenarios such as indoor scenario.The authors of [7,8] provided some wide band measurement results, including time delay dispersion and Ricean K factors.In [1,2] the authors noted that antenna alignment and the shape of the tunnel cross section have a great influence on the correlations between antennas, the channel matrix singular values, and MIMO capacity.However, to fully explore the advantages of using a large number of antennas in the tunnel, more detailed information regarding medium or large scale MIMO system performance is required.Additionally, dense measurements in the spatial domain also provide information that may be used to improve MIMO capacity when there are limits to antenna alignment precision or antenna selection algorithms [9].
The authors of [10,11] regarded the tunnel as a nonideal waveguide, but this makes it difficult to consider various factors such as the effects of cross section shape, curves, wall materials, and metal reflecting objects in the tunnel.Therefore, some authors employ ray tracing, for example, [12,13].Nonetheless, more work should be done to verify the simulation accuracy through comparison with measurements.This paper focuses on such a comparison of wide band MIMO channel characteristics with ray tracing results.
The organization of the paper is as follows.In Section 2, the measurement campaign is summarized.In Section 3, the measurement results, including path loss, time delay dispersion, and the array amplitude matrix, are analyzed and compared with ray tracing results.The relationships among the angle spread, singular values, and channel capacity are also illustrated.Section 4 presents the conclusion.

Description of the Measurement Campaign
2.1.Measurement Environment.Our measurement campaign was conducted in the Tunnel Laboratory in Zhongtian Technology Group Company, which is located in Nantong City, Jiangsu Province.The laboratory was built to test the performance of the company's leaky cables.The tunnel is made of reinforced concrete and consists of two parts of different cross-sectional shape, rectangular and circular; see Figure 1.The length of each part is 50 m.The interior width and height of the rectangular part are 5 m and 3 m, respectively.The interior radius of the circular part is 2 m and the height is 3 m; the circular section is truncated at its base.The receiver is fixed at the end of the rectangular part and the transmitter was moved to several locations in the tunnel.

Measurement
Equipment.The measurement system is shown in Figure 2. The transmitter (Tx) is a Rohde & Schwarz commercial signal generator SMBV100A.The receiver was made by our research group at Beijing Jiaotong University using several components.This receiver (Rx) equipment can acquire and store the IF signal after downconversion by the RF module.A GPS synchronization signal is transmitted by an optical fiber to both Tx and Rx.The antenna is a biconical antenna.This measurement system is mainly used in virtual MIMO measurements.Table 1 lists the primary measurement system parameters.

Probe Signal and Data
Processing.Zadoff-Chu (ZC) sequences are often used as synchronization signals due to their good orthogonality and low peak-to-average power ratio properties.The complex value of each sequence is given by [14] where [  ] is received signal at time   ,  is the multipath delay,  and  are time and delay indices, and (  ) is the locally generated ZC sequence.
As noted, the length of the ZC sequence used here is 2047.Hence the maximum delay that can be measured is (1/91 MHz) × 2047 = 2.25 × 10 −5 s.This value is much larger than the delay where any echoes can incur within the tunnel.So the receiver noise threshold can be estimated by the last (long-delay) part of the PDP where no echoes can exist.Since the noise threshold varies with Tx-Rx distance, a fixed threshold −40 dB relative to the power delay profile (PDP) peak

SISO Measurement
3.1.1.Path Loss.To measure the path loss (PL) in the tunnel, the Tx was moved from one end to the other in 1 m intervals.
The PL is given by where  Tx is transmit power and  Rx is received power, which is given by [15,16] where  is the number of valid multipath components (MPCs).
The path loss results obtained from the SISO measurement and those from the free space model and ray tracing simulations, all versus linear distance, are shown in Figure 4.The ray tracing tool we used here is Wireless InSite.X3D model was chosen in the simulation and diffractions are included.The accuracy in the ray tracing simulation depends on the reconstruction of the measurement environment and parameter settings established prior to simulation.Thus the ray tracing result is only approximate due to imperfect representation of the tunnel structure and electrical properties.The relative permittivity and conductivity values are 5 and 0.001 S/m, respectively.The number of reflections accounted for in the ray tracer is 8, and no more than 300 rays are employed in the simulation.
The path loss can be divided into two sections: before and after a distance of approximately 45 m.For the short range region, the measured path loss is up to 70 dB.In this short range section, the measured path loss is similar to, but slightly less than, the path loss of free space.After reaching 70 dB at 45 m, the measured path loss fluctuates around a value of approximately 60 dB.The drop at 45 m hence appears to come from the reduction of the cross section at the middle of the tunnel (transition from rectangular to circular), where the waveguide effect evidently increases; this waveguide effect is also expected to be stronger at longer link distances.The steady trend after 45 m is attributable to low order reflected multipath components (MPCs) adding to the LOS component; these MPCs are strongly affected by waveguide effect.The ray tracing results exhibit the same basic trend as the measurements but show more fluctuations in both sections.We point out that these path loss results did not average out small scale fading effects; hence the results are not used for path loss modeling but rather to observe overall trends in this two-section tunnel.

RMS Delay Spread. The RMS-DS is given by [15]
where  2 and  can be written as where () is the power measured at delay  (the PDP) and the mean excess delay  is the first moment of the PDP.
Figure 5 shows the RMS-DS results versus link distance from measurements and ray tracing with and without the reflection from the rear door (wall at 100 m in Figure 3); each point represents a RMS-DS result obtained from an individual PDP.In Figure 5(a), both curves show the same general trend of gradual rising and then falling back gradually; this agrees with intuition in such a "closed" structure.The maximum RMS-DS is approximately 150 to 160 ns at the middle of the tunnel and after a link distance of 60 to 70 m, the RMS-DS rapidly declines to about 20 ns. Figure 5(b) illustrates that after removing the rear door reflection, the RMS-DS stays at a low value of approximately 20 ns.Clearly any such long reflections as those from the rear door will have a strong influence on RMS-DS.As we can see, in both Figures 5(a we found that 8 reflections are appropriate for reaching a reasonable match between measurement and simulation.Also it is worth pointing out that the tunnel walls are not perfectly smooth, and since we did not have detailed information on the tunnel construction, we could not model the detailed material structure of the walls (e.g., the possible presence of reinforcing metal rods, so-called "rebar").As seen in comparing Figure 5(a) with Figure 5(b), removing the back wall reflections does improve the agreement between measurements and simulations.The normalized RMS error εRMSE between measurement and ray tracing results is given by where σrmsM and σrmsR are RMS-DS values that are normalized by the mean RMS-DS in each case.The normalized RMS errors in case of Figures 5(a) and 5(b) are 0.70 and 0.42, respectively.

Virtual MIMO Measurement.
As Figure 6 shows, the channel was measured at six positions along lines orthogonal to the tunnel longitudinal axis, at several link distances.At each position, 16 by 16 MIMO are used and at link distances of 25 m and 75 m, three Tx antenna heights (2.3 m, 2.4 m, and 2.5 m) were used.

PDP Results
. Figure 7 is a measurement and ray tracing PDP example when the Tx was placed at 25 m.Because the tunnel is closed, the echoes from the walls of the end of the tunnel can be found in the measured CIRs.As Figure 7 shows, part A is the useful signal, the echoes of part B are from the junction of the two tunnels, and part C is the reflection from the rear-wall and metal door at the end of tunnel II.It is obvious that in the ray tracing PDP, echoes in parts B and C can be distinguished from those in part A easily, according to the time delay.However, as link distance increases, the MPCs of parts B and C come gradually closer to those of part A, which makes it difficult to identify the valid MPCs of part A from the PDP.Thus, in cases like the PDP at 25 m in Figure 7, MPCs of parts B and C can be removed completely, but in cases like the PDP at 85 m (shown in Figure 13), we can only remove the signal under the noise floor, in which case some reflections from the rear door cannot be removed completely, and this will have some influence on the subsequent parameter analysis.This discussion also pertains to the RMS-DS results in Figure 5.
The measured PDPs for link distances from 15 m to 85 m appear in Figures 8-13.All the PDPs in these figures pertain to the case when both the Tx and Rx antennas are located near the tunnel center.Here number 8 transmit antenna to number 8 receiving antenna PDPs are shown.
In Figures 8-13, it can be observed that most of the useful signal has delay less than 200 ns, and as the distance increases, the path length difference (relative delay) from the first to last strong multipath component becomes smaller, as expected from the geometry.In each figure,  rms is the RMS-DS of the whole PDP, and  rmsA is the RMS-DS of part A in each PDP.
The PDPs show that the reflections from the rear door are strong at each position and can be separated from the LOS component.It also can be seen that in tunnel I the duration of reflected multipath from the rear door, as part C in Figure 7 shows, is less than approximately 100 ns, which corresponds to a 30 m path length difference.Since the tunnel height and width are 3 m and 5 m, this means that there can be 6 to 10 reflections when the signal propagates from transmitter to receiver from the rear-wall reflections.Thus we note again that 8 reflections were chosen as a reference when we set the ray tracing simulation parameters.

Amplitude Matrix.
In general, in small scale MIMO systems, antenna correlation and the MIMO capacity are of primary concern.But in medium or large scale MIMO systems, the received signal strength is also of interest, because when the number of antennas increases, distinct fading across the array may appear in some scenarios.In addition, in confined spaces like indoors and in tunnels, the capacity of any MIMO system is always influenced by specific characteristics of the antenna deployment.Thus novel metrics may be of use to help us understand how the signal propagates; ideally these metrics would be straightforward to compute and intuitive.An amplitude matrix (AM) is introduced here to show the effects of varying the antenna position.
Table 2 lists the number of multipath components using our fixed threshold and the resolved number of multipath   components extracted by SAGE.The comparison indicates that SAGE algorithm can slightly improve the resolution tunnel scenario.For simplicity multipath components were obtained by fixed threshold method.After removing the noise components, the multipath complex components above the noise floor are combined in the delay domain as where   is the amplitude of the channel between the  th Tx antenna and  th Rx antenna,  = 1, . . ., ,  = 1, . . ., .So the MIMO amplitude matrix Ã is given by

International Journal of Antennas and Propagation
For purposes of comparison, we use the normalized MIMO amplitude matrix, which is written as The normalized MIMO amplitude matrices are displayed from 15 m to 85 m in Figure 14.These show the normalized amplitude matrices of the virtual MIMO array at the various link distances.An interesting phenomenon can be observed:  most of the amplitude matrix plots indicate that the signal is transmitted mainly in a diagonal way, except for the case in Figure 14(d).The corresponding simulation results in Figures 15(a)-15(f) show the same phenomenon even more clearly than in Figures 14(a)-14(f).Figures 16(a the antenna position has a great impact on the received signal strength in these tunnel scenarios.To utilize these propagation characteristics, if the MIMO system is deployed in the tunnel, a suitable antenna selection algorithm may help the communication system reach its maximum capacity.
Ray tracing results for the amplitude matrices are shown in Figure 15.These figures reveal that ray tracing results can only roughly agree with the measurement results.Results in Figure 15 are, by visual inspection, the closest to measurements we have been able to obtain by carefully varying the representation of the tunnel structure and the parameters of the tunnel walls, as previously noted in discussing the PDPs and delay spread.In spite of that, the principal diagonal transmission character at each position is even more strongly present in the simulations.Furthermore, as the link distance increases, the number of diagonal "bands" decreases.We hypothesize that this is because the angular spread decreases with distance.
To help us judge the simulation accuracy, we have computed two metrics to compare the agreement between amplitude matrices of simulations and measurements: the root-mean-square error (RMSE) (11) and vector correlation (12).In (11) and ( 12)  and  are the vectorized versions of the measured amplitude matrix Α and simulated amplitude matrix Α, where  AM is total number of elements in the amplitude matrix.
As an absolute error, it is difficult to judge the simulation accuracy solely from the RMSE values, so the correlation is used here as our primary metric.For these metrics, good agreement between simulation and measurements means that  RMSE is small, and  corr is large (near unity).Table 2 shows the RMSE and correlation results, and the latter are always larger than 0.8, indicating at least reasonable simulation accuracy.
The results listed in Table 3 show the relative agreement between the measurement and simulation at fixed points in space.Since small scale effects may increase the actual error in such a confined space, an average of the amplitude  17(a)-17(f) show the measurement and simulated amplitude average results over 4 antennas (1.5 wavelengths) when averaging is done at the Tx and Rx and both at Tx and Rx antennas at 25 m.From these figures and the compilation of the estimation error metrics in Table 4 we observe that the local area amplitude matrices match even better than those for the individual array elements, especially after averaging at both Tx and Rx.

AOA and AOD.
The azimuth angle characteristics of the channel can be extracted from our linear array measurements [17].The angle of arrival (AOA) and angle of departure (AOD) are estimated using the space-alternating generalized expectation-maximization (SAGE) algorithm [18], and the results are shown in Table 5.As can be seen, when the Tx and Rx are both in the rectangular tunnel section (at 15 m, 25 m, and 35 m), the range of the AOA and AOD is about 60 degrees.At the longer link distances (65 m, 75 m, and 85 m), the range of the AOD decreases and is less than that of the AOA.Interestingly, the range of AOA and AOD both decreases with increasing link distance, as expected, but the AOA is larger than the AOD (except for the shortest link distance of 15 m).The RMS angle spread (RMS-AS) is used to quantify the change of the AOA and AOD.It is given by [19]  where  2 and  can be written as where   and (  ) are the angle and power of the  th multipath component.Figure 18 shows the AS versus link distance and summarizes the results of Table 5.We point out that the received signal is somewhat sensitive to the activities of the testers (i.e., the presence and motion of the humans in the environment): this caused a sudden change of the RMS-AS at 75 m, which also can be observed in the prior RMS-DS results and the subsequently provided singular value results.When the Tx and Rx are both located in the rectangular tunnel, the RMS-AS of AOA and AOD is approximately 13 to 16 degrees; when the Tx and Rx are located in differently shaped sections of the tunnel (and link distance is larger), the RMS-AS of AOA is ∼8 to 10 degrees and the RMS-AS of AOA and AOD is ∼4 to 7 degrees.

MIMO Capacity.
The capacity of a wireless communication link depends on the properties of the complex channel matrix , which is an  ×  dimensional matrix: where h is the complex channel coefficient that can be obtained by In general, a MIMO system's capacity is inversely proportional to the channel correlation between antennas.A low correlation between channels can be achieved by rich  scattering from a large number of reflectors; increasing physical separation between antennas also generally decreases correlation.In tunnel scenarios though, due to the waveguide effect, attenuation of high order modes will reduce the number of the propagating modes, which can be interpreted as being somewhat similar to the key hole effect, and this also reduces the rank of the channel matrix [1] .The singular  values of  are often used to quantify the key hole effect, and these will also have direct impact on the MIMO capacity.
According to [20], an  ×  matrix can be put into the form where  × and  × are two distinct unitary matrices and Σ × is a positive semidefinite diagonal matrix.The diagonal elements of Σ are the singular values of , which are represented by   , where  = 1 ⋅ ⋅ ⋅ , and  is equal to the minimum of  and .
Curves in Figure 19 show the normalized singular values, which are obtained by   / max(  ).This figure shows that distance has no obvious or monotonic influence on the normalized singular value, but rather it appears that the crosssectional area has the primary influence: the normalized singular value in the circular tunnel (longer distances) is less than that in the rectangular tunnel (shorter distances).
The maximum capacity of a memoryless  ×  MIMO channel with additive white Gaussian noise can be written as [21] where  is the average signal-to-noise ratio (SNR),   is the identity matrix of size ,  is the channel matrix, the superscript () Hermitian conjugate, and det(⋅) denotes the determinant of a matrix.To illustrate behavior, the ergodic capacity for a SNR of 10 dB is listed in Table 6,  6 and 7 we can observe that the rank of the channel matrix decreased more in circular tunnel (4 to 7) than in rectangular tunnel (2 to 3).Furthermore, it can be hypothesized that increasing the cross-sectional area appears to improve the MIMO link performance.

Conclusion
In this article, results of a wide band SISO and virtual MIMO measurement in a tunnel laboratory at a frequency of 1.4725 GHz have been presented.In the SISO case, fairly good agreement was found between the measurement and ray tracing results for path loss and RMS-DS characteristics.
The path loss presents strong waveguide effects for distances larger than approximately 45 m.The RMS-DS is near 20 ns without the rear door reflection and such remote reflectors can obviously increase the RMS-DS in tunnel scenarios.
In the MIMO case, the measurement was conducted at six link distances.The simulated PDPs show fairly good agreement with the measurements; a slight vertical shift of one antenna does not affect the PDP significantly.Amplitude matrices were introduced, and these show strong regularity characteristics at all of the six positions; specifically a strong diagonal transmission characteristic was found.This observation held more strongly when local averaging was applied to the amplitude matrices, and this may be of use in increasing the MIMO system performance.The simulation accuracy of the amplitude matrices based on ray tracing was analyzed as well, and the result also showed acceptable agreement between measurements and simulations.Angular characteristics and channel matrix singular values were also illustrated, and these were in accordance.Finally, although there are rich reflections in the tunnel, the MIMO capacity shows some resemblance to key hole effects according to the narrow angular characteristics and channel matrix singular values.

Figure 4 :
Figure 4: Path loss versus distance for measurement, ray tracing, and free space results.

Figure 5 : 15 m
Figure 5: RMS-DS versus distance for measurement and ray tracing results.

Figure 16 :
Figure 16: Illustrations of the two propagation modes at 25 m and 35 m based upon the amplitude matrices.

Figure 17 :
Figure 17: Measurement (a)-(c) and simulated (d)-(f) amplitude average results over 4 antennas on Tx and Rx and both on Tx and Rx antennas at 25 m.

Figure 18 :
Figure 18: RMS-AS of AOA and AOD versus link distance.

Figure 19 :
Figure 19: Normalized singular values at each value of link distance.

Table 2 :
Number of multipath components using fixed threshold and SAGE.

Table 3 :
Estimation error metric values versus link distance (m).

Table 4 :
Estimation error metrics after averaging.

Table 5 :
AOA and AOD at each position.

Table 6 :
Ergodic capacity versus link distance.

Table 7 :
Rayleigh channel from 8 × 8 to 16 × 16 MIMO capacity. is the capacity for an outage probability 0.5.As a comparison, the Rayleigh channel capacities for MIMO array sizes from 8 × 8 to 16 × 16 are given in Table7as well.As with the singular value results, the capacities in the circular tunnel section (longer distances 65-85 m) are smaller than those in the rectangular tunnel (shorter distances), with the exception of the 85 m value.As expected, the capacity results have strong positive correlation with the singular values.And also, by comparing Tables which