1. Introduction
The track quality not only directly affects the safe operation of the train [
1] but also its interaction performance with infrastructures [
2]. Irregularities of the track lead to an increased force input for both the vehicle and the track itself, are in this way self-reinforcing and therefore increase without appropriate measures. On the infrastructure side, the interactions are determined by the amplitude and the wavelength of the track discontinuity [
3]. Effects with short wavelengths can cause rail fractures and are therefore relevant to safety [
4]; at the same time, the elimination of the effects is cost-relevant, as Muhamedsalih demonstrates on squats [
5]. Noise problems caused by short-wave effects, as discussed by Li [
6], can reduce the attractiveness of the railway system for passengers and residents. In addition, according to Chiengson, the increased force effect due to the high-frequency excitation of vehicles leads to increased deterioration processes of the entire track [
7].
To make short-wave effects measurable, measurement systems with high-frequency sampling rates have to be used. Historically, measurement systems used by most infrastructure operators evaluate the infrastructure in longer wavelength ranges. The European standard EN 13848-1 covers three wavelength ranges: D1 (3–25 m), D2 (25–70 m) and D3 (70–150 m) [
8]. D3 is only relevant for speeds above 250 km/h, while D1 and D2 must be used obligatorily for asset management. The annex of the standard also refers to possible measurement methods for short-wave effects, primarily by means of axle box accelerations. However, this part of the standard is declared to be purely informative and is not specified as a minimum requirement.
A requirement that comes with the handling of short-wave and thus punctual faults is exact positioning. Since the timely development of a quality parameter must be considered for a meaningful condition assessment, various measurement runs must be positioned to each other in sufficient quality. Historically, location deviations of measurement runs in the meter scale were justifiable, as characteristic values were calculated over a length of 100 m and more [
9]. In addition, today the condition of the track geometry of open track is successfully assessed in this way. Even if the longitudinal level itself is used as a quality index [
10], positioning accuracies of 25 cm are sufficient, since this parameter is usually not stored in a higher resolution anyway. To evaluate short-wave effects, we need relative positioning accuracies in the centimeter to millimeter scale. The sector’s experience in describing short-wave effects with measuring data is relatively low due to these issues of missing precise position. Nevertheless, also today these effects are assessed by different methods.
The most common way of evaluating short-wave effects (the term “short-wave effect” is defined differently in various works. For this paper, the term is used for all effects of the railway track whose characteristics have a wavelength of less than 1 m. This includes any kind of corrugations, squats, welding joints, insulated joints, switch components, rail breaks, skid marks, indentations, break-outs at the running edge and machining errors) is still visual in-track inspection. Due to the high competence of the responsible track staff, visual inspection still delivers the most reliable results. Nevertheless, significant disadvantages of visual inspection occur. Regular route inspections are carried out at walking pace and are therefore very slow, resulting in high costs on the one hand and low frequencies on the other one. As the trend moves towards less available staff, track condition assessment can no longer be carried by track inspections only in the future. In addition, as the number of trains increases, entering the danger zone must be reduced to a minimum.
One very broad investigated method for the evaluation of short-wave effects are axle box accelerations. Li gives an example for the detection of squats [
11], and Molodova describes the possibility for evaluating the health condition of insulated rail joints [
12]. Axle box accelerators can be applied at very high speeds and still provide high sampling rates. Additional track unavailability can be avoided by mounting the devices on axles of in-service trains or regular measuring cars. Axle box accelerations are a direct consequence of short-wave effects in the wheel–rail contact point. On the other hand, accelerations are not only triggered by rail condition but also by the condition of the wheels. While a wheel that is flat produces periodic peaks that can be handled relatively well, the influence of wheel out-of-roundness and indentations is rather unpredictable. In addition, the speed of the vehicle has an influence on the amplitude of the measurements making it difficult to compare the signal characteristics [
13].
High-resolution images of rails are another method for rail defect detection with great potential, but has not yet been successfully implemented beyond trial operations. This procedure needs intelligent image processing algorithms. De Ruvo and Hovad give an overview on the principles of this method [
14,
15]. Despite a great deal of research, no system has been fully proven in operation yet.
It is also possible to evaluate short-wave effects using lasers. There are different possibilities of laser arrangement and subsequent data processing. The measurement system of the rail surface described in detail in the method section corresponds to one of these designs.
Vidovic could show that it is also possible to use Fibre Optic Sensing for the detection of short-wave defects [
16]. In his thesis, Vidovic demonstrated that the system provides information on the condition of insulated rail joints and turnout crossings. Research on distributed acoustic sensing is ongoing.
There are several other concepts conceivable. Papaelias gives an overview on measuring principles like ultrasonic transducers, magnetic induction, eddy current technologies, radiography, long-range ultrasonics and laser ultrasonics [
17].
Although all the described methods are fundamentally suitable for the assessment of short-wave effects, they have not been able to establish themselves in practice for a variety of reasons. This is partly due to the fact that the achievable measurement speeds are clearly too low, which makes a network-wide implementation unfeasible. This applies, for example, to many of the methods described by Papaelias. For various technical reasons, fiber optic sensing and image recognition tools are not yet developed far enough to make reliable statements. It is also not expected that all relevant effects can be evaluated with these systems. The same applies to axle box accelerations. These have already been intensively researched and provide reliable results for certain scenarios. However, they do not only react to the effects of the rail, but also to vehicle properties and stiffness changes of the track. A reliable classification of effects only with axle box accelerations is not considered possible, at least not at present. In addition, many of the methods raise the question of the relative and absolute positioning of the data.
Exact positioning of the measurement data is required for all measurement principles. Wang distinguishes absolute positioning errors and relative positioning errors [
18]. Since absolute positioning errors describe the deviation between the positioning and the absolute position of a point on the track, they are particularly important for planning maintenance activities. Relative errors, on the other hand, describe the deviation of the positioning between different measurement runs and are relevant for the investigation of time series. For short-wave effects, relative errors are primarily relevant, but the absolute error can also be very important. Additionally, Wang describes a way to significantly decrease the relative positioning errors using big-data fusion and incremental learning [
18].
Khosravi also deals with the positioning of track geometry measurements [
19], using different alignment methods including cross-correlation function, recursive alignment by fast Fourier transformation, dynamic time warping, correlation optimized warping and a combined method. The combined method has proven to be the most reliable, but also the most complex.
The data sources of the underlying research of this paper are measurement systems of the EM 250, the main track measurement car of OeBB-Infrastruktur AG (Austrian Federal Railways). Fellinger has developed an algorithm allowing for a price post-positioning of the majority of the measured area with the help of signal characteristics [
20]. The rail surface signal dealt with in this paper was not included in the method due to the significantly higher sampling rate and the requirement of an even more precise positioning. Nevertheless, Fellinger’s method is adapted to the data format of the measuring car and therefore forms the basis for the method described later on.
2. Positioning Method
For the detection and evaluation of short-wave effects, we use the original output of an available measurement system. This data is obtained from the rail surface measurement system, a chord-based system in which three optical distance sensors are mounted in a row.
Figure 1 illustrates the basic setup of the laser sensors.
The first and the last laser build a virtual chord. The actual measurement value describes the distance between the rail surface of the middle laser and the chord surface at the same position. The measurement system can be mounted onto vehicles with speeds up to 250 km/h and still returns a data point every 5 mm. The recorded wavelength range is between 0.01 and 1 m. The combination of a sampling rate of 5 mm and the high frequency spectrum of the system makes us confident for possible use cases, as other measurement systems (in Austria) cannot provide this.
The asymmetrical chord principle with three measuring points leads to very stable measuring results. Whereas a direct measuring method would lead to a lack of accuracy if vehicle movements occur, the chord principle can handle these very well. Linear correlated shifts influence all three outputs in the same way and are thus automatically compensated. Small pitch motions lead to linear correlated but opposite shifts of the first and last laser and therefore do not influence the output of the chord value. The chord measurement system’s output-signal reflects the corrugation of the railhead with distorted amplitudes. The amount of distortion depends on the wavelengths of the output signal characteristics and can be extracted from the transfer function of the measuring system. If the actual space curve of the rail corrugation is needed, relatively complicated algorithms have to be used for the recalculation of the space curve. These algorithms can only deliver an approximated space curve and thus a lower data quality. As for many evaluations, it is sufficient to work with chord-values and signal characteristics of chord values; pros and cons of a recalculation of the space curve must be weighed against each other.
Although, the measurement system has been mounted onto the EM250 since 2005 and some modified parameters have been partly in use since then, the original output signal was not investigated up to now. The main reason might be the already mentioned precondition of an exact positioning process. Without re-positioning, the potential of the small sampling rate cannot be exploited.
For re-positioning we use an algorithm shifting signals of different measurement runs synchronously to each other. For this purpose, signal characteristics are used and a synchronicity criterion is defined. We chose the sum of the Euclidean distance as it delivers the most reliable results. To determine the best possible relative position of two measurement runs, the signal of one measurement run is shifted sequentially by a certain distance and the similarity measure to the second run is calculated anew each time. In contrast to Fellinger, the signal is shifted by only 25 mm per iteration. The relative displacement with the smallest summed Euclidean distance is therefore carried out for all measurement runs with synchronous measurements as result. To demonstrate synchronicity,
Figure 2 shows input and output of the computing process.
The different colors of the signals in
Figure 2 reflect different runs of the measuring car. In this case five measurement runs are processed corresponding to about 1.5 years. The lines in the upper half of the figure represent the filtered signal without exact positioning. While the signal characteristics appear to match each other, the different runs are displaced to each other. The lower half of the figure shows the same signals after precise positioning. Here, all signal characteristics match each other.
Although the described method works in principle, it is computationally expensive and time-consuming due to the large amount of data (one data point every 5 mm) and the high number of iterations (to catch outliers, all potential displacements ±25 m must be covered leading to 2000 iterations per measurement run). For this reason, a process step is carried out upstream guaranteeing a positioning accuracy of ±1 m and reducing the computing time by a multiple.