Numerical Analysis of Low-Cost Recognition of Tunnel Cracks with Compressive Sensing along the Railway

: Currently, the use of microseismic detection technology for crack detection and localization in rock masses has great potential in detecting structural damage. As engineering safety has always been a very important issue, this study investigated the problem of multi-crack identiﬁcation in rock masses within the environment of track tunnels using transient waves. A tunnel rock was modeled using MIDAS GTS NX software (2019.v1.2) and a crack transient wave model in the frequency domain was obtained through data analysis and simulation. Then, this was combined with compressive sensing techniques to locate and detect multiple cracks in tunnel rock. The performance of the proposed approach was validated through experimental simulations, which included experiments on differences in the number of cracks, as well as spatial samples. The experimental results indicate that the technique performs well for single-crack localization in tunnel rock mass, where the average localization error is 4 m. Meanwhile, the localization error is larger in multi-crack localization, and the number of spatial sample points set using compressive sensing also has a large impact on the


Introduction
In the field of track tunnel engineering, engineering safety holds exceptional significance.Strengthening the real-time monitoring of tunnel rock mass stress is a common practice in track tunnel construction to prevent accidents, such as tunnel collapse [1].Currently, microseismic detection technology is widely used in tunnel disaster recognition and crack monitoring.When cracks appear in rock, they release strain energy and propagate elastic waves to the surroundings [2][3][4][5].Crack detection based on visual and optical instrument assistance is currently the main means of daily inspection [6,7], but the method is time consuming and labor intensive, and the detection data carry a certain degree of subjectivity and are not sensitive to early cracks.There have been many research works on crack localization and identification; some of them have focused on the study of elastic wave propagation in a single medium [8,9], while others have analyzed the propagation and attenuation characteristics of R-waves in non-uniform media [10,11].Most of them have focused on the study of propagation characteristics of elastic waves and few work on specific crack localization detection techniques.
In recent years, theoretical models for leak detection in pipelines in frequency-domain transient waves (pressure) have been studied [12].Frequency-domain transient waves have great performance in terms of localization detection accuracy.In the research on crack identification, significant progress has been made in various aspects, including (1) analysis using the reverse transient method [13][14][15][16], (2) integration of transient damping techniques [17][18][19], (3) utilization of reflection techniques in the time domain [20][21][22], and (4) identification based on frequency-domain responses [23][24][25][26].In this study, the aforementioned frequency-domain methods of transient waves were applied to crack localization and monitoring in tunnels.By conducting experimental simulations, the effect of localization detection can be simulated, which greatly simplifies the estimation problem of multiple cracks.
Despite the aforementioned work, there are still many unresolved issues in multicrack identification.Solving the identification and localization problem of multiple cracks is challenging; the number of cracks in practical applications was previously unknown.Previous studies have mainly focused on identifying and locating individual cracks to simplify the process.As the number of cracks increases, the computational complexity grows rapidly.As a result, a new low-cost approach was needed to make the estimation process less sensitive to the number of cracks.
With the development of compressive sensing (CS) technology [27], it has been increasingly utilized in various fields.In simple terms, CS leverages the sparsity of signals in a domain to achieve accurate signal reconstruction at low sampling rates, thereby breaking the conventional Nyquist sampling criterion under specific conditions.Currently, CS is used successfully in image processing and communications [28][29][30][31][32][33].In the context of crack detection in tunnel rock mass within the field of railway transportation, there has been limited work.Due to the spatial sparseness of the crack locations in the whole tunnel, the transient wave model for multi-crack identification can be linearly represented in terms of crack size but not in terms of crack location [26].Nonetheless, the actual cracks within the tunnel can be considered sparse; therefore, investigating the problem of multi-crack localization using a sparse structure became a viable approach.
By establishing a transient wave model in the frequency domain for tunnel rock masses, using the geotechnical analysis software MIDAS GTS NX 2019v1.2, a transient wave model suitable for track tunnel rock masses was simulated and obtained the corresponding pressure transfer function through data analysis, which can reflect the propagation relationship.
This work was a CS-based multi-crack identification technique, due to the nonlinearity of the measurement model with respect to crack parameters, directly modeling and estimating the data, which obtained from real-world applications can be highly challenging.To address this issue, spatial sampling techniques were employed to transform the application problem into a defined space, thereby converting the problem of crack localization into a mathematical problem for modeling and solution.
Using MATLAB software (R2023a) with simulation analysis, the facilitation of further evaluation of the performance of this technique can be obtained.For a single-crack case, the average localization error is approximately 4 m.For the localization of two or more cracks, the experimental simulation results show that the average localization error is related to the distance between multiple cracks, and when the distance between multiple crack points is short, larger errors or missed detections will be caused.

Materials and Methods
The pressure-domain transient wave model for tunnel rock masses is elucidated and the presence of multiple cracks within the rock mass can be simulated.Based on this, the mathematical relationship between the tunnel cracks and the transmission of pressuredomain transient waves can be derived.

Transient Wave Model in the Frequency Domain of the Pressure of Tunnel Rock
Considering a tunnel rock mass with dimensions of 100 m in length, 50 m in width, and 50 m in height.Through simulation, modeling of the rock mass and setting its properties was based on real-world mechanical parameters.For illustration purposes, this model simulated a common rock mass with a density of 2730 kg/m 3 , an elastic modulus of 11 GPa, and a Poisson's ratio of 0.21.The simulation was conducted using the Midas GTS NX (2019.v1.2) geotechnical analysis software.A grid system within the tunnel was created, with 50 observation points (c i , where i = 1, 2, . . ., 50) evenly distributed along the Appl.Sci.2023, 13, 13007 3 of 15 100 m length of the rock mass.At one end of the rock mass, dynamic loads in the form of concentrated forces were applied.Table 1 lists some of the simulated data Then, we performed experimental simulations, which consisted of two parts: Simulating the perfect rock mass under dynamic loads of different angular frequencies to obtain the corresponding data Z 0 (c i ; ω); simulating the practical rock mass with cracks under dynamic loads of different angular frequencies to obtain the corresponding data Z(c i ; ω).
During the process of experimental simulation, the model design was based on the geological properties of the tunnel rock mass and the experimental requirements.Figure 1 illustrates the application of dynamic loads at the position and direction indicated by the red arrow in the tunnel, with loads applied at different angular frequencies.Figure 2 demonstrates the introduction of cracks within the rock mass, where a spherical cavity is used to simulate the cracks, representing their location and volume as an indication of their size.concentrated forces were applied.Table 1 lists some of the simulated data Simulating the perfect rock mass under dynamic loads of different angular fr cies to obtain the corresponding data Z c ; ω ; simulating the practical rock ma cracks under dynamic loads of different angular frequencies to obtain the corresp data Z  ;  .
During the process of experimental simulation, the model design was based geological properties of the tunnel rock mass and the experimental requirements.1 illustrates the application of dynamic loads at the position and direction indica the red arrow in the tunnel, with loads applied at different angular frequencies.F demonstrates the introduction of cracks within the rock mass, where a spherical c used to simulate the cracks, representing their location and volume as an indica their size.

Data Analysis and Derivation
Based on the pressure frequency-domain transient wave transmission funct scribed in [26], we adopted this model to establish the transmission law of pressure the rock mass.Reference [12] proposed a pressure function related to position and frequency within the tunnel:

Data Analysis and Derivation
Based on the pressure frequency-domain transient wave transmission function described in [26], we adopted this model to establish the transmission law of pressure within the rock mass.Reference [12] proposed a pressure function related to position and angular frequency within the tunnel: here, q U (ω) represents the head oscillation, Z(ω) is the characteristic impedance, µ(ω) is the propagation function, consisting of two constants with angular frequency as the independent variable, s m is the distance from the loading point to the rock mass, and ω is the angular frequency of the applied dynamic load.This function is position-and angular frequency-dependent, representing the theoretical pressure value in the absence of cracks at the distance s m from the loading point.The specific physical parameters of the formula are provided in Appendix A.
Proposing an expression for the pressure frequency-domain transient wave model of the tunnel rock mass: This expression is derived by simplifying Equation ( 1) and representing the constants as four unknown constants: a, b, c, and d.These four constants are the unknowns, which we aim to determine through experimental data.The observed data sets are used for machine learning training using the least squares method, with 20% of the data set used for testing.By substituting the test set data into Equation (2) for solving, the four unknown constants a, b, c, and d are determined as 1.5724, 3.4860 × 10 −6 , 2.5129, and 9.2144 × 10 −2 , respectively.Testing based on the training set achieved an MSE of 0.00197, indicating good accuracy of the pressure frequency-domain transient wave model for the tunnel rock mass.
Based on the simulated pressure frequency-domain transient wave model of the tunnel rock mass, along with the modeling of the rock mass, frequency-domain transient wave analysis can be performed, providing the transmission relationship of pressure frequency-domain transient waves with respect to position (as given by Equation ( 2)).This provides theoretical support for transforming the pressure difference vector into a vector representation of crack position and size.

Tunnel Rock System Model
A transient wave model was developed for locating and detecting cracks in the tunnel rock in a rail transit environment and a compressive sensing (CS) model for multi-crack detection problems was proposed.
Considering a rock mass (geotechnical properties are set as in Section 2, corresponding to common soft rock tunnels) with N cracks, as shown in Figure 3, where setting the first and last ends to s U and s D , setting x n and v n to be the position and size information of the n-th crack, ∀n = 1 : N, and assuming s A total of M pressure sensors were deployed along the tunnel to collect observational results for crack identification, where s m denotes the location of the mth sensor installation, and the location is a section of the rock body that emits a dynamic load point.

Tunnel Rock System Model
A transient wave model was developed for locating and detecting cracks in the tunnel rock in a rail transit environment and a compressive sensing (CS) model for multi-crack detection problems was proposed.
Considering a rock mass (geotechnical properties are set as in Section 2, corresponding to common soft rock tunnels) with  cracks, as shown in Figure 3, where setting the first and last ends to  and  , setting x and v to be the position and size information of the -th crack, ∀ = 1: , and assuming s < x < x ⋯ < x < s .A total of  pressure sensors were deployed along the tunnel to collect observational results for crack identification, where  denotes the location of the mth sensor installation, and the location is a section of the rock body that emits a dynamic load point.Crack detection employing the frequency-domain transient wave model [12].Based on the description in Section 2, it can be established that the pressure values, at the detection points, vary depending on the position and size of the cracks [26].Thus, the pressure values can serve as the primary observations for determining the location and size of the cracks, which will be utilized in subsequent model development.
The data on pressure is collected by  sensors for subsequent crack detection.  denotes the data collected by sensor  at angular frequency .The set of frequencies used in collecting data at Ω has a total of  angular frequencies Setting the pressure values  (), ∀ ∈ Ω by sensor  , the pressure difference can be expressed as: where  () is the theoretical pressure at  when there are no cracks within the range 0, s , as given by Equation (2).
Expressing d () as a vector  ∈ ℝ : here, setting  and  as the information about the location and size of the cracks in the tunnel, which can be obtained as: where x < ⋯ < x for the position .
From the transient wave model derived in [12], a representation of the pressure difference of the crack related to the crack size and location information is written.
This linearly represents the crack size but not the crack position. represents measurement errors, () ∈ ℝ × is the measurement matrix, and the composition of the measurement matrix is related to the location information of the crack point, which can be expressed as: Crack detection employing the frequency-domain transient wave model [12].Based on the description in Section 2, it can be established that the pressure values, at the detection points, vary depending on the position and size of the cracks [26].Thus, the pressure values can serve as the primary observations for determining the location and size of the cracks, which will be utilized in subsequent model development.
The data on pressure is collected by M sensors for subsequent crack detection.H m (ω) denotes the data collected by sensor s m at angular frequency ω.The set of frequencies used in collecting data at Ω has a total of N Ω angular frequencies Setting the pressure values {H m (ω), ∀ω ∈ Ω} by sensor s m , the pressure difference can be expressed as: where H 0 m (ω) is the theoretical pressure at s m when there are no cracks within the range [0, s m ], as given by Equation (2).
Expressing d m (ω) as a vector d ∈ R N Ω M : here, setting x and v as the information about the location and size of the cracks in the tunnel, which can be obtained as: where From the transient wave model derived in [12], a representation of the pressure difference of the crack related to the crack size and location information is written.
This linearly represents the crack size but not the crack position.represents measurement errors, ρ(x) ∈ R N Ω M×N is the measurement matrix, and the composition of the measurement matrix is related to the location information of the crack point, which can be expressed as: where ϕ(x N ), n = 1 : N, is the measurement function ψ m (x n ; ω) dependent on the crack position x n , given by: and ψ m (x n ; ω) is provided in [26], with the expression: In Equation ( 9), the relevant unknown variables are represented in Equation ( 2).Based on the modeling analysis of the tunnel rock mass in this Section, the transient wave model in the tunnel rock mass allows for compressive sensing (CS) sparse modeling.This transforms the practical application problem into a mathematical problem for solution.

CS Sparse Model
In this Section, the process of transforming rock transient wave modeling to coefficient modeling is described, transforming the application problem into the problem of finding a sparse recovery signal.The corresponding column vectors of the basis matrix of the spatial sample points α i are: where ψ m (α i ; ω) depends on α i and is given by Equation ( 8) with α i substituted.Let G(α) ∈ R N Ω M×I be the base matrix, defined as: where N Ω represents the number of frequencies utilized.Consequently, the measured value d: The error s in this is due to the fact that the spatial sample points set during the spatial sample modeling cannot be aligned with the sensor sampling points.The vector t represents the sparse crack size signal, where its non-zero elements correspond to the actual crack sizes, and its support corresponds to the spatial samples closest to the actual crack positions.When setting up the spatial sample points, a larger number of spatial sample points are generally set up to obtain better results, so a small number of cracks in the tunnel rock mass can be regarded as sparse for the spatial sample point.As shown in Figure 4, this illustrates the graphical representation of sparse modeling for crack detection in tunnel rock mass, where d = G(α)t + + s .From the figure, it can be observed that the measurement value d can be expressed as the product of G(α) and the estimated value t, where each column of G(α) corresponds to φ(α i ).The spatial sample points corresponding to the non-zero values in t represent possible crack locations in the rock mass, and the magnitude of the non-zero values represents the crack sizes.

Estimating Crack Positions 𝒙 and Sizes 𝒗
In the expression of Equation ( 12), by incorporating the known transient wave model, it is possible to represent the tunnel cracks by the recovered vector, which simultaneously provides information about the crack locations and sizes.Within the recovered vector, all non-zero entries correspond to crack points, with the position information reflected by the index of the entry, and the value at that position indicating the size of the respective crack; therefore, the focus is to estimate  from , translated into questions 0.Problem 0 is explained by the following lemma, where  ( , ) =  ∈ ℂ : ∥  ∥ ≤ 2 represents the set of sparse signals:

Lemma 2. [For the error analysis of a single crack] There exists an optimal solution to the problem for a single crack based on a specific basis matrix 𝑮(𝜶) under the transient wave model in the frequency domain, which has the following constraints:
The proof of Lemma 2 is provided in Appendix C. The optimal solution to the problem 0 is the one that is closest to the real signal with enough spatial sample points.By modeling the crack detection problem in tunnel rock mass using CS, the problem of crack detection is transformed into seeking the optimal solution of 0.

Crack Detection Method Based on CS
In this Section, the logic of the sparse recovery algorithm and analysis of the algorithm are described.Then, based on the estimated value , the crack position  and crack size  are determined.

Estimating Crack Positions x and Sizes v
In the expression of Equation ( 12), by incorporating the known transient wave model, it is possible to represent the tunnel cracks by the recovered vector, which simultaneously provides information about the crack locations and sizes.Within the recovered vector, all non-zero entries correspond to crack points, with the position information reflected by the index of the entry, and the value at that position indicating the size of the respective crack; therefore, the focus is to estimate t from d, translated into questions P0.

(P0)
Here, ||t|| 0 is the 0 -norm, is the number of non-zero terms in the recovery vector, d − G(α)t 2  2 is the 2 -norm, and d th is The set error threshold, d th = (1 + ς z )σ 2 , where ς z ∈ [0.1, 0.3] and σ 2 is the variance of the measurement error .
Problem P0 is explained by the following lemma, where D (I,2) = µ ∈ C I : µ 0 ≤ 2 represents the set of sparse signals: the value of δ (2) is defined as (λ max − λ min )/(λ max + λ min ), and C (2) is calculated as (λ max + λ min ) /2, where λ max and λ min are obtained from Appendix B. The proof of Lemma 1 can be found in Appendix B.

Lemma 2. [For the error analysis of a single crack]
There exists an optimal solution to the problem for a single crack based on a specific basis matrix G(α) under the transient wave model in the frequency domain, which has the following constraints: The proof of Lemma 2 is provided in Appendix C. The optimal solution to the problem P0 is the one that is closest to the real signal with enough spatial sample points.By modeling the crack detection problem in tunnel rock mass using CS, the problem of crack detection is transformed into seeking the optimal solution of P0.

Crack Detection Method Based on CS
In this Section, the logic of the sparse recovery algorithm and analysis of the algorithm are described.Then, based on the estimated value d, the crack position x and crack size v are determined.

Reconstruction of Sparse Signals
Problem P0 is a standard sparse recovery problem.For the selection of sparse recovery algorithms, the Orthogonal Matching Pursuit (OMP) algorithm is employed.According to Equation (12), given G(α) and d, the goal is to recover t from these values.It is evident that we need to fully utilize the sparsity of u.Based on the knowledge of linear algebra, d can be represented as a linear combination of the column vectors of matrix G(α), where t acts as a weight for the column vectors of G(α).As t is sparse, only a few column vectors in G(α) significantly contribute to d.The objective is to identify the column vectors that have a significant impact on d.Simultaneously, based on the positions of the column vectors in G(α), the positions of the non-zero elements in t can be determined.This is the fundamental application of the OMP algorithm in the crack estimation.Summarizing the pseudocode for the OMP algorithm to recover t, the specific algorithm steps are as Algorithm 1 below:

1.
Set algorithm parameters, such as I and the setting of d th 2.
Initialize the index set S = ∅ and the set of spatial samples A = {1, . . . ,I} 3.
Set the measured value d = d 4.
Loop through the following steps until d 2 ≤ d th

Determination of Crack Position and Size
Given t, the crack position x and size v are determined.The computational steps used are to first set up an index set of spatial samples, expressed using the following Equation ( 14), which needs to be taken into account when setting up the index set, with respect to the previous values of the error thresholds: R{•} is used to represent the part of the real part of the complex number; subsequent problems on crack estimation are analyzed based on the part of the real part, and the estimate obtained by the specific algorithm can be expressed as follows: ti represents the ith estimated value of t.Due to the error in the measurement estimation, there will be a projection of the error in the OMP basis matrix.So when setting the threshold t th , the value should be set to be larger than this projection of the error, but should not be too large in order to produce false detections and missed detections in the test.If t th is set too large, it may lead to missed detections where small cracks may go undetected.On the other hand, if t th is too small, it may cause false detections.
The recovery vector calculated using the algorithm is able to reflect the real crack situation to a certain extent; specifically, the estimated vector can be correlated with the real crack sparse signal throughout the transformation.For the algorithm, the signal with good recovery effect, it corresponds to the non-zero spatial sample point, which is actually very close to the real crack.This is also related to the setting of the algorithm parameter, as well as the setting of the spatial sample point; its ith element t th will be equal to v n , otherwise it will be zero for ∀i ∈ S. The OMP algorithm essentially comes up with a search for column vectors on G(α S ) that contribute more to the measurements, while at the same time, based on the position of the column vectors in G(α S ), the position of the non-zero elements in t can be determined.

Simulation Results
Simulating the leakage detection method based on compressive sensing (CS) to validate the experimental performance of the proposed scheme.

Simulation Setup
Considering the overall properties of the tunnel rock mass, a rectangular model with a length of 100 m and both width and height of 50 m was constructed.A scenario that involves a maximum of two cracks was built, with crack sizes ranging from 5 × 10 −4 m 2 to 10 × 10 −4 m 2 .For this purpose, we adopted that I = 1000 uniformly distributed spatial samples, and equally spaced M = 50 sensors along the 100 m length of the tunnel rock mass for data acquisition.The frequency set Ω was defined as The transmission relationship of frequency-domain transient waves in tunnel rock mass is explained in Section 2 of this study, based on the same material properties and relevant parameters, with the transfer function given by Equation (2).

Analysis of the Crack Detection Results of the OMP Algorithm
Analyzing the results of the experimental simulation in terms of a single crack as ll as two cracks.

Analysis of the Results for a Single Crack
The localization of individual cracks was simulated and analyzed.The simulated estimation of the leak location is shown in Figure 5, where the number of spatial sample points determines the localization accuracy.The average error in crack location analysis for a single crack was 4.3547 m.The method was not applied in tunnel crack detection, and a comparison of the average errors in [12] shows that the simulation error is decreased, and the cumulative distribution functions (CDFs) of the estimated crack position errors are shown in Figure 6.
As shown in Figure 5, the simulation results obtained using the CS-based Orthogonal Matching Pursuit (OMP) algorithm are presented.The widely spaced points in the figure represent the measurement nodes, which are the locations where the sensor detection data are obtained.The closely spaced nodes represent the spatial sample points, which correspond to the potential crack locations.
As shown in Figure 6, the CDF of the error obtained from experimental simulation for a single crack indicates the distribution of the error.The error, denoted as |x − x|, represents the distance between the estimated position and the actual location of the crack, reflecting the accuracy of the estimation.timation of the leak location is shown in Figure 5, where the number of spatial sample points determines the localization accuracy.The average error in crack location analysis for a single crack was 4.3547 m.The method was not applied in tunnel crack detection, and a comparison of the average errors in [12] shows that the simulation error is decreased, and the cumulative distribution functions (CDFs) of the estimated crack position errors are shown in Figure 6.As shown in Figure 6, the CDF of the error obtained from experimental simulation for a single crack indicates the distribution of the error.The error, denoted as |x − x|, represents the distance between the estimated position and the actual location of the crack, reflecting the accuracy of the estimation.

Analysis of Results for Two Cracks
The simulated estimation of crack positions for two cracks is shown in Figure 7.It can be observed that when the distance between the two cracks is close, the localization error increases.This indicates that the resolution for the performance of crack identification is degraded, which may be related to the wavelength of the transient wave, i.e., it is difficult to identify two cracks if the distance between them is smaller than the wavelength of the transient wave.The cumulative distribution functions (CDFs) of the estimated crack position errors are shown in Figure 8.

Analysis of Results for Two Cracks
The simulated estimation of crack positions for two cracks is shown in Figure 7.It can be observed that when the distance between the two cracks is close, the localization error increases.This indicates that the resolution for the performance of crack identification is degraded, which may be related to the wavelength of the transient wave, i.e., it is difficult to identify two cracks if the distance between them is smaller than the wavelength of the transient wave.The cumulative distribution functions (CDFs) of the estimated crack position errors are shown in Figure 8.As shown in Figure 6, the CDF of the error obtained from experimental simulation for a single crack indicates the distribution of the error.The error, denoted as |x − x|, represents the distance between the estimated position and the actual location of the crack, reflecting the accuracy of the estimation.

Analysis of Results for Two Cracks
The simulated estimation of crack positions for two cracks is shown in Figure 7.It can be observed that when the distance between the two cracks is close, the localization error increases.This indicates that the resolution for the performance of crack identification is degraded, which may be related to the wavelength of the transient wave, i.e., it is difficult to identify two cracks if the distance between them is smaller than the wavelength of the transient wave.The cumulative distribution functions (CDFs) of the estimated crack position errors are shown in Figure 8.As shown in Figure 8, the CDF analysis of the error using the CS-based OMP algorithm for two cracks reveals that the detection error for two cracks is obviously greater than that for a single crack.During the simulation process, the distance between the two crack points has a significant impact on the error.The closer the crack points are, the larger the error will be, and it may even fail to detect the two cracks.The detection error for multiple cracks is greater than that for a single crack.The chosen spatial sample size I also has a substantial impact on the simulation results, as increasing the number of spatial samples can improve the accuracy of localization detection simulation to some extent.

Discussion and Conclusions
A compression-aware tunnel crack localization technique using the OMP recovery algorithm in the field of railway transportation is proposed.This CS-based crack detection method relies on the pressure difference at a point in the rock mass, considering both the As shown in Figure 8, the CDF analysis of the error using the CS-based OMP algorithm for two cracks reveals that the detection error for two cracks is obviously greater than that for a single crack.During the simulation process, the distance between the two crack points has a significant impact on the error.The closer the crack points are, the larger the error will be, and it may even fail to detect the two cracks.The detection error for multiple cracks is greater than that for a single crack.The chosen spatial sample size I also has a substantial impact on the simulation results, as increasing the number of spatial samples can improve the accuracy of localization detection simulation to some extent.

Discussion and Conclusions
A compression-aware tunnel crack localization technique using the OMP recovery algorithm in the field of railway transportation is proposed.This CS-based crack detection method relies on the pressure difference at a point in the rock mass, considering both the The detection error for multiple cracks is greater than that for a single crack.The chosen spatial sample size I also has a substantial impact on the simulation results, as increasing the number of spatial samples can improve the accuracy of localization detection simulation to some extent.

Discussion and Conclusions
A compression-aware tunnel crack localization technique using the OMP recovery algorithm in the field of railway transportation is proposed.This CS-based crack detection method relies on the pressure difference at a point in the rock mass, considering both the presence and absence of cracks.As the linear representation of crack positions is not feasible, CS for modeling analysis, utilizing the OMP algorithm for signal recovery, is employed, followed by experimental simulations to analyze the results.The specific conclusions are as follows: 1.
The existing frequency-domain transient wave analysis method is applicable to tunnel rock mass environments.Through data analysis and simulation, the existing frequency-domain transient wave model is applied to the rail transportation field to solve the crack detection technology of tunnel rock, which has scalability.

2.
According to the proposed crack detection method, a more efficient scheme can be proposed for sensor deployment based on the results obtained from the simulations.
Sensor placement in practical engineering can be guided by the expected error range, providing practical guidance for engineering applications.

3.
The accuracy of CS-based crack detection primarily depends on the number of spatial samples.These simulation results demonstrate that increasing the number of spatial sample points can reduce the error.4.
The resolution of the proposed crack identification technique is related to the wavelength of the transient wave.If the distance between the two cracks is smaller than the wavelength of the transient wave, it becomes challenging to distinguish between the two crack points.5.
In the experimental simulations for localization detection, the employed physical formulas and corresponding relationships are currently based on theoretical derivations.
To improve the accuracy of localization, further research should focus on enhancing the content related to these aspects.
The CS-based framework utilizing the transient wave model shows potential for crack detection in tunnel rock masses.The findings of this study provide insights into the application of CS-based methods; different models should be used to study specific rock properties and future works will analyze the theoretical performance of CS-based tunnel crack detection.

Figure 1 .
Figure 1.Dynamic load applied to a perfect rock mass in normal conditions.

Figure 1 .
Figure 1.Dynamic load applied to a perfect rock mass in normal conditions.

Figure 2 .
Figure 2. Dynamic load applied to a perfect rock mass with cracks.

Figure 2 .
Figure 2. Dynamic load applied to a perfect rock mass with cracks.

Figure 3 .
Figure 3.A tunnel rock mass with cracks and sensors.

Figure 3 .
Figure 3.A tunnel rock mass with cracks and sensors.

4. 1 .
Sparse Modeling of Pressure Differences The purpose is to estimate the location and size of the crack based on the collected data d.Directly estimating x and v is challenging due to the nonlinear relationship between d and x and v. Proposing to use spatial sampling to transform a portion of the linear model into a sparse representation simplifies the associated signal processing.Consider I spatial samples of crack positions α = vec[α i | ∀i = 1 : I], where α i ∈ [s U , s D ].These sample positions α represent potential crack positions.

Algorithm 1 :
The OMP-Based Sparse Signal Recovery Input: Measurement values of pressure difference d.

6 .||ν i || 2 7 .
Find i such that i = arg max i∈A S Gain the index set S as S = i ∪ S 8. Update the estimated value t as t = G(α)) † d , where G(α)) † is the pseudo-inverse of G(α) Calculation result of d = d − G(α S ) t Output: Estimated value t

Figure 5 .
Figure 5. Monitoring results for a single crack.

Figure 5 .Figure 6 .
Figure 5. Monitoring results for a single crack.As shown in Figure 5, the simulation results obtained using the CS-based Orthogonal Matching Pursuit (OMP) algorithm are presented.The widely spaced points in the figure represent the measurement nodes, which are the locations where the sensor detection data are obtained.The closely spaced nodes represent the spatial sample points, which correspond to the potential crack locations.

Figure 7 .
Figure 7. Monitoring results for two cracks.

Figure 6 .
Figure 6.CDF of estimated error in single-crack position.

15 Figure 5 .Figure 6 .
Figure 5. Monitoring results for a single crack.As shown in Figure 5, the simulation results obtained using the CS-based Orthogonal Matching Pursuit (OMP) algorithm are presented.The widely spaced points in the figure represent the measurement nodes, which are the locations where the sensor detection data are obtained.The closely spaced nodes represent the spatial sample points, which correspond to the potential crack locations.

Figure 7 .
Figure 7. Monitoring results for two cracks.

Figure 7 .
Figure 7. Monitoring results for two cracks.

Figure 8 .
Figure 8. Cumulative distribution function of estimated error in two crack positions.6.2.3.Analysis of Results for Different Spatial Sample Sizes (I) Then, the performance of CS-based crack detection with different spatial sample sizes was analyzed.The range of values for the spatial sample number I is from 100 to 3000.The resulting crack position and size estimation errors relative to I are shown in Figure 9.The simulation results demonstrate that the crack detection performance improves as the number of spatial samples increases.

Figure 9 .
Figure 9. Relationship between crack detection error and spatial sample size.

Figure 8 .
Figure 8. Cumulative distribution function of estimated error in two crack positions.

6. 2 . 3 .Figure 8 .
Figure 8. Cumulative distribution function of estimated error in two crack positions.6.2.3.Analysis of Results for Different Spatial Sample Sizes (I)Then, the performance of CS-based crack detection with different spatial sample sizes was analyzed.The range of values for the spatial sample number I is from 100 to 3000.The resulting crack position and size estimation errors relative to I are shown in Figure9.The simulation results demonstrate that the crack detection performance improves as the number of spatial samples increases.

Figure 9 .
Figure 9. Relationship between crack detection error and spatial sample size.

Figure 9 .
Figure 9. Relationship between crack detection error and spatial sample size.

Table 1 .
Position-dependent pressure data at specified angular frequency (200π as an example).

Table 1 .
Position-dependent pressure data at specified angular frequency (200π as an exam