Comparative in situ study of dynamic load generated by gravel piles measured by a fiber-optic interferometer

Citation: Stolárik, M.; Kepák, S.; Pinka, M.; ˇ Cubík, J.; Nedoma, J. Comparative In Situ Study of Dynamic Load Generated by Gravel Piles Measured by a Fiber-Optic Interferometer. Sensors 2022,22, 5579. https://doi.org/10.3390/s22155579 Academic Editors: Swee Chuan Tjin and Lei Wei Received: 8 July 2022 Accepted: 21 July 2022 Published: 26 July 2022 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). sensors Article Comparative In Situ Study of Dynamic Load Generated by Gravel Piles Measured by a Fiber-Optic Interferometer Martin Stolárik 1, Stanislav Kepák2,* , Miroslav Pinka 1, Jakub ˇ Cubík2and Jan Nedoma 2 1Department of Geotechnics and Underground Engineering, Faculty of Civil Engineering, VSB—Technical University of Ostrava, 708 00 Ostrava, Czech Republic; [email protected] (M.S.); miroslav[email protected] (M.P.) 2 Department of Telecommunications, Faculty of Electrical Engineering and Computer Science, VSB—Technical University of Ostrava, 708 00 Ostrava, Czech Republic; [email protected] (J. ˇ C.); [email protected] (J.N.) *Correspondence: stanislav[email protected] Abstract: Currently, all the technology used for seismic monitoring is based on sensors in the electrical domain. There are, however, other physical principles that may enable and fully replace existing devices in the future. This paper introduces one of these approaches, namely the field of fiber optics, which has great potential to be fully applied in the field of vibration measurement. The proposed solution uses a Michelson fiber-optic interferometer designed without polarization fading and with an operationally passive demodulation technique using three mutually phase-shifted optical outputs. Standard instrumentation commonly used in the field of seismic monitoring in geotechnical engineering was used as a reference. Comparative measurements were carried out during the implementation of gravel piles, which represents a significant source of vibration. For the correlation of the data obtained, the linear dependence previously verified in laboratory measurements was used. The presented results show that the correlation is also highly favorable (correlation coefficient in excess of 0.9) from the values measured in situ, with an average deviation for the oscillation velocity amplitude of the optical sensor not exceeding 0.0052. Keywords: vibration measurement; gravel piles; anthropogenic sources; seismic monitoring; civil engineering; fiber-optic sensor; interferometer 1. Introduction Seismic monitoring is an integral part of the implementation of most geotechnical structures where vibration-generating equipment is used. These vibrations can adversely affect both the surrounding buildings in terms of adverse dynamic effects on the construction of building structures, which is dealt with in technical standards at the national or European level, as well as the comfort of the population, when the generated vibrations have a negative effect on the human body, which is covered by sanitary and government regulations. From the point of view of technical seismicity acting on buildings, we monitor the maximum amplitudes of velocity or acceleration obtained on the basis of in situ seismic measurements in the building and the prevailing frequencies. The magnitude of the amplitude and frequency pattern depend on the source of the vibration, the distance of the object being monitored, and the geological environment through which the vibration propagates through the foundation joint and foundation structure to the upper structure. As the distance from the source of dynamic effects increases, the intensity of these vibrations decreases, so one of the basic tasks of seismic engineering is also to determine the attenuation characteristics of the rock environment. An important source of vibration in geotechnical engineering is mainly equipment designed to modify the properties of the rock environment, such as vibratory packers, vibratory plates, and vibratory rollers, which accelerate the primary consolidation process through vibration, as well as deep foundation technologies in the form of drilled, driven, or Sensors 2022,22, 5579. https://doi.org/10.3390/s22155579 https://www.mdpi.com/journal/sensors Sensors 2022,22, 5579 2 of 20 vibrated piles. In the case of these technologies, we are usually talking about harmonic oscillations with a predominant frequency determined by the manufacturer of the equipment used and with oscillation speed amplitudes in the range of tenths to tens of millimeters per second. For seismic monitoring, i.e., monitoring the effects of technical seismicity, standard seismic instrumentation consisting of a three-component velocity or acceleration sensor operating on the electrical principle and a digital seismic station is used. This measurement system has been used for several decades, but in the developing 21st century, other physical principles are also becoming available that could successfully complement or even replace the established standard with a number of benefits. In the case of successful implementation in seismic engineering, several advantages can be found, be it the design of the sensor itself, resistance to external influences, and even price. The results of previous experimental studies dealing with the use of different physical principles in the monitoring of the effects of engineering seismicity have confirmed the applicability of the alternative devices under development, especially in the frequency domain. The logical follow-up was to focus on the time domain, where a first (laboratory) experiment performed in a homogeneous environment with a calibrated strike yielded interesting and thought-provoking results [ 1 ]. This served as the impetus for a larger-scale follow-up experiment, in which vibrations generated during the implementation of in situ gravel piles were monitored using standard seismic instrumentation, which acted as a reference in the experiment, and by an optical interferometric sensor under development. The article presents the entire experiment carried out on the construction of a road bypass in the Czech Republic, the results of this experiment with emphasis on correlations in the time domain, and the subsequent determination of the attenuation characteristics of the rock environment using data obtained from the alternative monitoring device under development. The fundamental innovation of the presented solution is mainly in applications in the field of anthropogenic vibrations. The research on the following showed basic experimental applications in the field of natural seismicity and vibrations around high-speed railways and tunnels, without comparison to a reference gauge. 2. State of the Art In this part of the article, the problem of seismic monitoring of vibrations from anthropogenic sources using commercially manufactured seismic instrumentation for this purpose is presented, as well as an alternative approach in the form of a physical principle based on optical fibers which can be used for vibration measurement. In conclusion, the actual implementation of the gravel piles whose dynamic effects were monitored is presented. 2.1. Commercially Manufactured Seismic Instrumentation Three-component commercially manufactured seismic stations operating in the electrical domain are commonly used throughout the world in a wide range of applications, be it the monitoring of natural seismicity, induced seismicity, or engineering seismicity due to anthropogenic activity. In the field of technical seismicity, which is dealt with by the authors, there are two major areas of interest. The first is the area of blasting, both in the excavation of underground works, especially shallow below the surface and in the intracity of large agglomerations [ 2 ], and also in the extraction of minerals in surface quarries by blasting, where these quarries are often located close to civil buildings, especially in Central Europe [ 3 ]. In these cases, the results of seismic monitoring are used to optimize the blasting operations: limit and total charge, as well as charge timing, in addition to monitoring the impact of the blasting or mining on structures. The second major area is the monitoring of dynamic expressions of heavy motor vehicles, especially rail transport, which have an undesirable effect on structures [ 4 ], or the monitoring of dynamic phenomena related to technological processes on construction sites, especially in the area of improving foundation conditions by compaction or the Sensors 2022,22, 5579 3 of 20 implementation of deep foundations [ 5 ]. In all these cases, seismic stations are usually installed as solitary temporary meters located in buildings according to applicable national legislation and standards. With the help of profile measurements, when using multiple seismic stations, it is then possible to derive the attenuation of a given environment from the measured values and to predict the range of dynamic loading in the vicinity of the vibration source based on the attenuation curves. Last but not least, seismic measurements are also used to calibrate advanced mathematical models. The common frequency range of seismic stations used is from 2 to 200 Hz, and the nature of the time domain recording depends on the particular source of dynamic loading. The fundamental disadvantage of these devices is the mechanical principle of the measurement itself, which makes the measuring instrumentation highly sensitive to handling. In contrast to newly proposed fiber-optic sensor technology, these commonly used devices are also not immune to electromagnetic influences and cannot withstand extreme climatic conditions in the long term. 2.2. Optical Sensors The authors are building on the results of previous work. Papers [ 6 , 7 ] focus on monitoring vibrations caused by blasting or other sources of dynamic loading. This uses fiber-optic sensors in combination with conventional vibration measurement equipment. As a result, the pilot studies provide a good starting point for the further deployment of these alternative devices. Other achievements include the measurement of vibration during landform compaction with a vibratory roller [ 8 ], the study of harmonic vibration in comparative measurements of alternative equipment and instrumentation for seismic monitoring [ 9 ], and research on vibration associated with works and construction technologies in brownfield areas [ 10 ]. Last but not least, the authors presented a description of new methods for seismic monitoring, namely a laboratory comparison of fiber-optic sensors with pneumatic measurement systems [1]. Vibration measurement with a focus on seismic monitoring (earthquakes, blasting, mine-induced seismicity, construction work, etc.) can be approached in various ways, one of which is fiber optic technology. These are classified according to the application, principle, technology, or design of the sensor itself. For example, a description of these technologies is discussed in [11]. Nowadays, many publications mainly focus on distributed measurements using optical fibers. A specific example is the distributed acoustic sensor (DAS), which evaluates the Rayleigh scattering coming back from a transmitting device. This way, it can turn a single-mode fiber into several thousand acoustic-vibration sensors. This system has been described and used by the authors of papers [ 12 , 13 ]. Distributed measurements using optical fibers are usually implemented on fibers for telecommunication purposes, but their cladding and mechanical robustness can be optimized [ 14 ]. The field of measurements with DAS is diverse: it can be, for example, the study of Earth tide frequencies for the definition of bedrock properties [15] or systems for seismic monitoring [16]. Another fiber-optic technology presented in the field of vibration measurement is fiber Bragg gratings (FBG). This is an inscribed periodic structure that changes the refractive index directly into the optical fiber. This structure can change the reflected wavelength on the mechanical stresses applied to it. This principle can be encountered, for example, in [ 17 ], where a three-dimensional accelerometer consisting of Bragg grating structures embedded in a seven-core optical fiber is presented. The technology of Bragg grates as seismometers, their production, and application are being addressed by various teams of authors, for example, the authors of publications [ 18 , 19 ]. They are analyzing their fiber-optic grating seismometers and improving their accuracy and fabrication methods [ 20 , 21 ]. The FBG measurement system is used in various applications, such as seismic activity monitoring in coal mines [ 22 ], or health structure monitoring [ 23 ]. The most scientifically active areas of FBG sensors are evaluation methods and their encapsulation [24,25]. Sensors 2022,22, 5579 4 of 20 The latest fiber-optic system technology suitable for seismic monitoring is phase sensors. Specifically, these are fiber-optic interferometers that evaluate the phase change in light between the measuring arm and the reference arm. These sensors can come in a variety of configurations and have one common denominator: high sensitivity to vibration. The first reports come from 2014, when the authors [ 26 ] present a fiber-optic Michelson interferometer-based accelerometer. However, the solution only describes the encapsulation of the measuring arm of the interferometer, with the optical fiber forming the outer winding around the transducer and it has only been verified in the laboratory. The sensor design was taken further by the authors of the publication [ 27 ]. The solution describes the housing of the entire sensor head, which now includes both the measuring and reference arms of the Michelson interferometer forming the fiber-optic accelerometer. Numerical simulation of the natural frequency with the result of 258.38 Hz had been carried out using the finite element method (FEM). Sensor performance was evaluated in the laboratory on a shaking table against a reference accelerometer. According to the authors, the sensor is potentially suitable for ocean floor seismic event monitoring; however, no experiments outside the laboratory were presented. For practical use, the encapsulation of the coupling element, whose fibers up to the sensor head are sensitive to vibrations and thus effectively form part of the sensor, will have to be solved. The three-dimensional interferometric system is used by the authors of publications [ 28 , 29 ], for example, to monitor seismic activity on the seabed. This type of sensor can be embedded in multicore fibers for compactness [30,31]. This sensor can be used, for example, as an accelerometer to measure low-frequency vibrations [ 32 , 33 ]. The Michelson fiber-optic interferometer as a seismic sensor can also be found in articles regarding various applications, such as the vibration measurement of oil wells [ 34 ] or general use in seismic applications [ 35 – 38 ]. In [ 35 ], a prototype device for seismological applications was developed. In this case, the light source is fed into the device via optical fiber, but the measurement system is assembled using bulk optics. The measurement of vibration disturbances in a subway tunnel during construction work is the subject of this article [ 37 ]. It is a proof-ofconcept technology without in-depth analysis of measured data and frequency analysis. The authors of the paper [ 38 ] developed a system to measure non-anthropogenic sources of vibration in three directions, who installed their equipment in the seismic observation cave at the Earthquake Administration Bureau. Three-component commercially manufactured seismic stations operating in the electrical domain are commonly used throughout the world in a wide range of applications, be it the monitoring of natural seismicity, induced seismicity, or engineering seismicity due to anthropogenic activity. In the field of technical seismicity, which is dealt with by the authors, there are two major areas of interest. 2.3. Gravel Piles There are two methods of installing gravel piles in the soil. The first is by means of a vibrating or vibroflotation device, and the second method is pre-driving by so-called Franki piles. During installation, vibrations are mainly generated and transmitted through the rock environment [39]. Gravel piles are used to improve the mechanical and physical parameters of cohesive and non-cohesive soils. For example, they can be used to improve the deformation properties of soils in existing road embankments and therefore limit the settlement of any new embankments intended for road widening [ 40 ]. Another application can be to optimize the transfer of stress from the superstructure to the soil. With the use of gravel piles, smaller foundation dimensions can be used [ 41 ] and the construction will have a lower economic and environmental impact. Gravel piles can also be used to modify the strength properties of the soil environment and thus may be used, for example, to stabilize potentially dangerous slopes [42]. The most typical use of gravel piles, however, is under new embankments to limit consolidation while strengthening the subgrade. The consolidation limitation is implemented Sensors 2022,22, 5579 5 of 20 by shortening the drainage path, which is reflected in the design of the pile grid [ 43 ]. For this purpose, gravel piles were also installed at the site where the experimental measurements presented in this article were carried out. 3. Methods This section of the paper focuses on the specific instrumentation used to measure the dynamic response of the subsurface from the effects of hammer-drop impact during the implementation of in situ gravel piles, both commercial seismic equipment and the fiber-optic system under development. 3.1. Compact Seismic Station BRS32 Two BRS32 seismic stations, see Figure 1, were used as the reference instrument at two sites: these had also been used by the authors in previous studies aimed at comparing results obtained from alternative approaches to vibration measurement [ 1 , 6 , 8 ]. The special mobile compact datalogger had originally been developed for the Institute of Rock Structure and Mechanics of the CAS of the Czech Academy of Sciences, v.v.i. and was subsequently put into custom production in the Czech Republic. This station is mainly used in the Visegrad Group countries for all seismic engineering tasks, especially in the fields of monitoring of technical and induced seismicity, and also in geophysics [44–46]. Sensors 2022, 22, x FOR PEER REVIEW 5 of 20 embankments intended for road widening [40]. Another application can be to optimize the transfer of stress from the superstructure to the soil. With the use of gravel piles, smaller foundation dimensions can be used [41] and the construction will have a lower economic and environmental impact. Gravel piles can also be used to modify the strength properties of the soil environment and thus may be used, for example, to stabilize potentially dangerous slopes [42]. The most typical use of gravel piles, however, is under new embankments to limit consolidation while strengthening the subgrade. The consolidation limitation is implemented by shortening the drainage path, which is reflected in the design of the pile grid [43]. For this purpose, gravel piles were also installed at the site where the experimental measurements presented in this article were carried out. 3. Methods This section of the paper focuses on the specific instrumentation used to measure the dynamic response of the subsurface from the effects of hammer-drop impact during the implementation of in situ gravel piles, both commercial seismic equipment and the fiberoptic system under development. 3.1. Compact Seismic Station BRS32 Two BRS32 seismic stations, see Figure 1, were used as the reference instrument at two sites: these had also been used by the authors in previous studies aimed at comparing results obtained from alternative approaches to vibration measurement [1,6,8]. The special mobile compact datalogger had originally been developed for the Institute of Rock Structure and Mechanics of the CAS of the Czech Academy of Sciences, v.v.i. and was subsequently put into custom production in the Czech Republic. This station is mainly used in the Visegrad Group countries for all seismic engineering tasks, especially in the fields of monitoring of technical and induced seismicity, and also in geophysics [44–46]. Figure 1. BRS32 seismic stations used as reference meters. The instrument contains three 31-bit A/D converters that allow the sampling of input signals in the frequency range of 250 Hz to 4 kHz and the storage of the measured data, together with data from a GPS receiver (time and position), which is part of the instrument, on a built-in flash drive with a capacity of up to 64 GB. A trio of SM6-3D seismometers with a frequency range of 4.5 to 100 Hz is housed inside the instrument as standard, Figure 1. BRS32 seismic stations used as reference meters. The instrument contains three 31-bit A/D converters that allow the sampling of input signals in the frequency range of 250 Hz to 4 kHz and the storage of the measured data, together with data from a GPS receiver (time and position), which is part of the instrument, on a built-in flash drive with a capacity of up to 64 GB. A trio of SM6-3D seismometers with a frequency range of 4.5 to 100 Hz is housed inside the instrument as standard, but the sensors can also be connected externally. For data transfer, the device offers a high-speed USB interface and, when connected to a computer, it acts similar to a regular USB drive. The entire seismic station is built into a rugged PeliCASE case that allows the station to be completely buried in the ground. The instrument is powered by an internal battery with a running time of up to 40 h, or from an external source for long-term seismic monitoring. 3.2. Fiber-Optic Measurement System under Development Fiber-optic sensing uses the physical properties of light as it travels along a fiber for a measurement of quantities. One of the physical properties of light is a phase. In an optical Sensors 2022,22, 5579 6 of 20 fiber, changes in external conditions (e.g., vibration) change the optical path length and thus the phase delay of the light passing through the fiber. The phase delay φ can be described by the following term [47], φ=2π λnL =βL, (1) where λ is the wavelength in a vacuum, nis the fiber core refractive index, Lis the fiber length and β is the propagation coefficient of the fiber. The vibrations result in a change in the core ∆ nand cladding refractive indices of the fiber due to photoelastic effect, its length ∆ Ldue to strain effect and, to a lesser extent, the fiber core diameter change ∆ ddue to Poisson effect. The phase change ∆φgenerated in the fiber can be then written as [48]: ∆φ=(∆β)L+β(∆L)=L∂β ∂n∆n+L∂β ∂d∆d+βL∆L L, (2) where dis the core diameter. The elasticity of the fiber can be expressed by the Young’s modulus Yand the Poisson ratio ν representing the transverse expansion (contraction) coefficient. According to the direction of the forces applied to the optical fiber, we distinguish axial strain, lateral radial and lateral unidirectional pressure. When a force is applied to a section of fiber in its axial direction, it will be stretched or compressed and an axial strain occurs. The lateral radial pressure applies the same force to the optical fiber from all directions, so the axial mechanical strains are the same for all directions. The fiber section under radial pressure is much longer than its transverse size, so the axial stress here can be neglected, and the axial strain can be deduced from the transverse stress. Under this condition, stresses and strains only exist in the x–y plane. Under lateral unidirectional pressure, which is the case of the sensor under development, the stress and strain distributions are usually analyzed in a circular disk model, compressed by two equal and opposite pressures applied at the ends of a fiber core diameter. Eventually, (2) can be rewritten as [48]: ∆φ=βLP Y(1−2ν)·1 2n2(p11 +2p12)−1, (3) where p 11 and p 12 are photoelastic coefficients of the optical fiber and the resulting phase shift is proportional to the external pressure Papplied. At a distance, all waves approach plane wave propagation, and the relationship between pressure and particle velocity is directly proportional (this is particularly valid for the sound waves). When transverse waves (typically vibrations) propagate in a medium, pressure remains constant in the direction of wave propagation, but the particles are displaced perpendicular to the direction of the wave propagation, and with that the pressure, which is then recorded by our sensor in the z-axis. Interferometers are then used to measure the phase shift between the light wave traveling through the sensor (or so-called measuring arm) and a reference wave split from the same radiation source. There are several possible configurations of the fiberoptic two-beam interferometer, and Michelson is one of the suitable configurations for vibration measurement. A fiber-optic coupler and a pair of mirrors form a Michelson interferometer. Between the coupler and the mirror, there is a measurement path and a reference path, respectively. For directional coupling and decoupling of light, a circulator or isolator is used. A photodetector is connected to the output, while the input is connected to the radiation source. The proposed system connection diagram of the fiber-optic measurement system is in Figure 2. The radiation source is the distributed feedback laser diode (LD) operating at a wavelength of 1550 nm and with an output power of 3 mW. A three-port circulator is used for directional coupling. The optoelectronic part is further composed of three pieces of photodetectors (PD) with p-i-n diodes made with InGaAs. The ADC used was the NI-9222 module in the cDAQ-9171 chassis made by National Instruments. Digital Signal Processing (DSP) and demodulation were performed on a laptop containing our custom-developed LabVIEW application. Sensors 2022,22, 5579 7 of 20 Sensors 2022, 22, x FOR PEER REVIEW 7 of 20 A fiber-optic coupler and a pair of mirrors form a Michelson interferometer. Between the coupler and the mirror, there is a measurement path and a reference path, respectively. For directional coupling and decoupling of light, a circulator or isolator is used. A photodetector is connected to the output, while the input is connected to the radiation source. The proposed system connection diagram of the fiber-optic measurement system is in Figure 2. The radiation source is the distributed feedback laser diode (LD) operating at a wavelength of 1550 nm and with an output power of 3 mW. A three-port circulator is used for directional coupling. The optoelectronic part is further composed of three pieces of photodetectors (PD) with p-i-n diodes made with InGaAs. The ADC used was the NI9222 module in the cDAQ-9171 chassis made by National Instruments. Digital Signal Processing (DSP) and demodulation were performed on a laptop containing our custom-developed LabVIEW application. Figure 2. Proposed system connection diagram. Interrogation of conventional two-beam interferometers is necessary as the phase shift is wrapped in the following term [46]: 󰇛󰇜=+  cos󰇟∆󰇛󰇜󰇠, (4) where C is the mean value of the optical intensity, A is the amplitude of the variation of the optical intensity and ∆ϕ(t) is the phase shift between the measurement and the reference path. A demodulation technique is essential to unwrap the phase shift and acquire the sensor data. Several different methods can be used for this purpose, one of which is homodyne demodulation [49]. The passive variant works on the principle of symmetry of the 3 × 3 coupler having three mutually phase-shifted outputs [47]: =+  cosΔ󰇛󰇜+2 3+󰇛󰇜, (5) where n = 1,2,3 and δn is 3 × 3 coupler phase asymmetry. Before the demodulation, it is necessary to equalize signal so the C and A values for all channels are Cn = 0 and An = 1. The phase change can be then easily unwrapped using a mathematical combination of two or three signals produced by the 3 × 3 coupler, in this case [50]: Δ󰇛󰇜=arctan󰇩 √ 3󰇛−󰇜 +−2󰇪. (6) The phase shift measurement and dynamic range are given by the photodetector bandwidth and the sampling rate of the ADC used. The sensor consists of a 3 × 3 coupler, and one pair of fiber mirrors with a fiber length of 3 m. The layout of the internal sensor can be seen in Figure 3. The outer case of the sensor consists of a waterproof aluminum box with dimensions 253 × 159 × 72 mm. A Figure 2. Proposed system connection diagram. Interrogation of conventional two-beam interferometers is necessary as the phase shift is wrapped in the following term [46]: I(t) = C+Acos[∆φ(t)], (4) where Cis the mean value of the optical intensity, Ais the amplitude of the variation of the optical intensity and ∆φ (t) is the phase shift between the measurement and the reference path. A demodulation technique is essential to unwrap the phase shift and acquire the sensor data. Several different methods can be used for this purpose, one of which is homodyne demodulation [ 49 ]. The passive variant works on the principle of symmetry of the 3 × 3 coupler having three mutually phase-shifted outputs [47]: un=Cn+Ancos∆φ(t) + 2πn 3+δn(t), (5) where n= 1,2,3 and δn is 3 × 3 coupler phase asymmetry. Before the demodulation, it is necessary to equalize signal so the Cand Avalues for all channels are C n = 0 and A n = 1. The phase change can be then easily unwrapped using a mathematical combination of two or three signals produced by the 3 ×3 coupler, in this case [50]: ∆φ(t) = arctan"√3(u2−u3) u2+u3−2u1#. (6) The phase shift measurement and dynamic range are given by the photodetector bandwidth and the sampling rate of the ADC used. The sensor consists of a 3 × 3 coupler, and one pair of fiber mirrors with a fiber length of 3 m. The layout of the internal sensor can be seen in Figure 3. The outer case of the sensor consists of a waterproof aluminum box with dimensions 253 × 159 × 72 mm. A measuring arm in the form of a fiber mirror was attached to the bottom of the box using epoxy resin. A reference arm of the same optical length (within the coherence length of the LD used) was loosely placed in vibroacoustic insulating foam to minimize the transmission of vibrations to this arm. Optical FC connectors were built into the front panel that represent the input and output ports of the sensor. The sensor was connected to the evaluation site by a 100 m-long multi-fiber-optic cable. The test measurements of a vibration test bench found that in the limited frequency range of 5–50 Hz given by the bench used (Netter NEG 5020), the observed sensor response was almost linear (statistically insignificant difference). The dominant frequency of the hammer impact lies in this frequency range and therefore the following amplitude–frequency analysis was not affected by the frequency dependence of the sensor sensitivity. Sensors 2022,22, 5579 8 of 20 Sensors 2022, 22, x FOR PEER REVIEW 8 of 20 measuring arm in the form of a fiber mirror was attached to the bottom of the box using epoxy resin. A reference arm of the same optical length (within the coherence length of the LD used) was loosely placed in vibroacoustic insulating foam to minimize the transmission of vibrations to this arm. Optical FC connectors were built into the front panel that represent the input and output ports of the sensor. The sensor was connected to the evaluation site by a 100 m-long multi-fiber-optic cable. Figure 3. Sensor internal layout. The test measurements of a vibration test bench found that in the limited frequency range of 5–50 Hz given by the bench used (Netter NEG 5020), the observed sensor response was almost linear (statistically insignificant difference). The dominant frequency of the hammer impact lies in this frequency range and therefore the following amplitude– frequency analysis was not affected by the frequency dependence of the sensor sensitivity. 4. Setup of Experimental Measurement The experimental measurements were carried out during the construction of the structure SO101 as part of the construction of the road 1/11—northern bypass of Opava (49.9562050N, 17.9048656E). The structure will consist of a transition embankment between two bridges that span the Ostrá stream and the local road 0468 towards Poland. During the design of the embankment, unsatisfactory properties of the underlying soil layers were discovered during the exploratory works. One of the biggest problems was the presence of groundwater from a nearby watercourse. As a result, different types of improvement of rock environment properties were proposed. One of the main measures included the construction of a set of gravel piles to strengthen the underlying soil and, at the same time, accelerate the process of consolidation of the underlying layers. Gravel piles are made of crushed or quarried aggregate. They are installed in mixed and fine-grained soils. The installation was carried out with a heavy crawler (Figure 4) and consisted of the drilling of a borehole in the soil, which was then filled with aggregate. The borehole is protected by a steel casing. The outer diameter of the casing is 520 mm. The total length of the borehole was 12 m. The internal void space of the pile is filled with loose aggregate. The aggregate is continuously forced into the soil surroundings, while the casing is simultaneously pulled up. Forcing is performed with a drop hammer. An aggregate of approximately 1.5 m3 is always added, and this is then forced and compacted by a 3.5 t drop hammer falling from a height of 8 m. With each blow, the aggregate is forced into the surrounding area to form a gravel pile with a minimum diameter of 600 mm. The shock generated by the drop hammer was the subject of the experimental measurements. The experimental measurement consisted of measuring the dynamic effects of the impact of the drop hammer on the surface layers of the soil. Measurements were carried Figure 3. Sensor internal layout. 4. Setup of Experimental Measurement The experimental measurements were carried out during the construction of the structure SO101 as part of the construction of the road 1/11—northern bypass of Opava (49.9562050N, 17.9048656E). The structure will consist of a transition embankment between two bridges that span the Ostrástream and the local road 0468 towards Poland. During the design of the embankment, unsatisfactory properties of the underlying soil layers were discovered during the exploratory works. One of the biggest problems was the presence of groundwater from a nearby watercourse. As a result, different types of improvement of rock environment properties were proposed. One of the main measures included the construction of a set of gravel piles to strengthen the underlying soil and, at the same time, accelerate the process of consolidation of the underlying layers. Gravel piles are made of crushed or quarried aggregate. They are installed in mixed and fine-grained soils. The installation was carried out with a heavy crawler (Figure 4) and consisted of the drilling of a borehole in the soil, which was then filled with aggregate. The borehole is protected by a steel casing. The outer diameter of the casing is 520 mm. The total length of the borehole was 12 m. The internal void space of the pile is filled with loose aggregate. The aggregate is continuously forced into the soil surroundings, while the casing is simultaneously pulled up. Forcing is performed with a drop hammer. An aggregate of approximately 1.5 m 3 is always added, and this is then forced and compacted by a 3.5 t drop hammer falling from a height of 8 m. With each blow, the aggregate is forced into the surrounding area to form a gravel pile with a minimum diameter of 600 mm. The shock generated by the drop hammer was the subject of the experimental measurements. The experimental measurement consisted of measuring the dynamic effects of the impact of the drop hammer on the surface layers of the soil. Measurements were carried out at two sites. The entire procedure of the installation of the gravel pile was monitored by continuous measurement. The installation of a pile consisted of the phases of arrival of the machine at a predetermined site, drilling to the required depth and then piling the aggregate and compacting it with a drop hammer. From the evaluation of the entire record, the phase with the largest amplitude was determined to be the hammer impact phase, which was an order of magnitude greater than the other phases described above. As a result, this phase was identified for evaluation and processing. Two machines were installing the gravel piles on the site. One was always performing the installation, while the other was moving to another station. There was also normal construction activity by the construction company at the site, but despite this seismic noise, the dynamic effect of the hammer strike was several times greater and therefore clearly identifiable. Sensors 2022,22, 5579 9 of 20 Sensors 2022, 22, x FOR PEER REVIEW 9 of 20 out at two sites. The entire procedure of the installation of the gravel pile was monitored by continuous measurement. The installation of a pile consisted of the phases of arrival of the machine at a predetermined site, drilling to the required depth and then piling the aggregate and compacting it with a drop hammer. From the evaluation of the entire record, the phase with the largest amplitude was determined to be the hammer impact phase, which was an order of magnitude greater than the other phases described above. As a result, this phase was identified for evaluation and processing. Two machines were installing the gravel piles on the site. One was always performing the installation, while the other was moving to another station. There was also normal construction activity by the construction company at the site, but despite this seismic noise, the dynamic effect of the hammer strike was several times greater and therefore clearly identifiable. Figure 4. Measurement site (red site A)—on the left optical sensor under development, in the middle BRS32. Figure A1a–e shows a schematic of the individual measurements. The measurements of a gravel pile installation were always performed at two sites (site “a”; site “b”). It was measured in five situations. In one situation, two gravel pile installations were always measured. Each installation was performed with a different machine of the same type. For this paper, situations where individual piles (pile 1; pile 2) were installed one after the other and their effects which could be clearly distinguished were always selected. However, there were also cases of simultaneous installation during the measurements, but these were discarded due to ambiguous interpretation. Figure A1 shows the distances (in mm) between the gravel pile installation being monitored and the site for each situation. For the sake of clarity, each distance is described by a name in the Mx.ya/b format. When M indicates a measurement, the first number indicates the situation, the second number indicates the pile’s serial number (1 or 2) and the letter (a or b) indicates if the distance to “a” or “b” is indicated. The example marking M4.1a (whose waveforms are shown in figures below) indicates a measurement at the fourth situation of the first pile at Figure 4. Measurement site (red site A)—on the left optical sensor under development, in the middle BRS32. Figure A1a–e shows a schematic of the individual measurements. The measurements of a gravel pile installation were always performed at two sites (site “a”; site “b”). It was measured in five situations. In one situation, two gravel pile installations were always measured. Each installation was performed with a different machine of the same type. For this paper, situations where individual piles (pile 1; pile 2) were installed one after the other and their effects which could be clearly distinguished were always selected. However, there were also cases of simultaneous installation during the measurements, but these were discarded due to ambiguous interpretation. Figure A1 shows the distances (in mm) between the gravel pile installation being monitored and the site for each situation. For the sake of clarity, each distance is described by a name in the Mx.ya/b format. When M indicates a measurement, the first number indicates the situation, the second number indicates the pile’s serial number (1 or 2) and the letter (a or b) indicates if the distance to “a” or “b” is indicated. The example marking M4.1a (whose waveforms are shown in figures below) indicates a measurement at the fourth situation of the first pile at position “a” with a distance of 9 m (this measurement configuration is marked red, see Figure A1d). 5. Results The following section presents the results of the experimental measurements in both the time and frequency domains from both the commercially produced BRS32 seismic station and the fiber-optic system under development. For further processing, only those data were selected that corresponded to the predriving phase in the Franki pile process, and the highest vibrations are generated in this phase. The set of five consecutive most powerful blows of the hammer was always chosen. Sensors 2022,22, 5579 16 of 20 7. Conclusions The paper presented an interferometric sensor developed for the in situ monitoring of vibrations due to anthropogenic sources and using linear correlation dependence for the basic conversion of measured values to SI units. The fundamental innovation of the presented paper is mainly the application of technical seismicity in the field of geotechnical engineering and comparation with standard instrumentation for seismic monitoring—it was used as a reference instrument. A total of 200 phenomena were evaluated in both the time and frequency domains, and the individual phenomena were similar in nature and character in both the frequency and time records. The correlation coefficient R2was 0.936 and the average deviation from the measured values of the oscillation velocity was 0.0052 mm·s−1, which is very low. The results of this experiment constitute part of the long-term research of the author’s team, the aim of which is a fully functional three-axis interferometric seismic sensor for use in geotechnical engineering. The conclusions from the in situ measurements show that the interferometric sensor (even if only one-axis) has the potential to replace standard seismic instrumentation even in field use, at least in an area where vibrations are generated by anthropogenic sources (harmonic vibrations of technological processes during geotechnical constructions, vibrations generated by traffic especially in cities and blasting in open pits near built-up areas or in excavated shallow tunnels), thus fully exploiting its many benefits. Author Contributions: Conceptualization, S.K. and M.S.; methodology, M.S.; software, S.K.; validation, M.S., J. ˇ C. and S.K.; formal analysis, M.S.; investigation, S.K., J. ˇ C. and M.P.; resources, M.S.; data curation, S.K. and M.P.; writing—original draft preparation, M.S., S.K., M.P. and J. ˇ C.; writing—review and editing, M.S., S.K., J. ˇ C. and J.N.; visualization, M.S., S.K., J. ˇ C. and M.P.; supervision, J.N.; project administration, S.K.; funding acquisition, M.S., S.K., J. ˇ C. and J.N. All authors have read and agreed to the published version of the manuscript. Funding: This article was supported by the Ministry of Education of the Czech Republic (Project No. SP2022/18) and by the project entitled “Monitoring of parameters of rocks disintegration of small and large in scale by fiber-optic sensors”, project number FW03010207 within the Technology Agency of the Czech Republic. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons. Conflicts of Interest: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. Sensors 2022,22, 5579 17 of 20 Appendix A Sensors 2022, 22, x FOR PEER REVIEW 17 of 20 Informed Consent Statement: Not applicable. Data Availability Statement: The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons. Conflicts of Interest: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. Appendix A (a) (b) (c) (d) Figure A1. Cont. Sensors 2022,22, 5579 18 of 20 Sensors 2022, 22, x FOR PEER REVIEW 18 of 20 (e) Figure A1. (a–e) Measuring situations 1–5 with distances between gravel piles and measuring stations marked. References 1. Stolarik, M.; Nedoma, J.; Martinek, R.; Kepak, S.; Hrubesova, E.; Pinka, M.; Kolarik, J. New Methods to Seismic Monitoring: Laboratory Comparative Study of Michelson Fiber-Optic Interferometer and Pneumatic Measurement Systems. Photonics 2021, 8, 147. https://doi.org/10.3390/photonics8050147. 2. Zhou, X.; Zhang, X.; Feng, H.; Zhang, S.; Yang, J.; Mu, J.; Hu, T. Study on Dominant Frequency Attenuation of Blasting Vibration for Ultra-Small-Spacing Tunnel. Appl. Sci. 2022, 12, 1058. https://doi.org/10.3390/app12031058. 3. Feher, J.; Cambal, J.; Pandula, B.; Kondela, J.; Sofranko, M.; Mudarri, T.; Buchla, I. Research of the Technical Seismicity Due to Blasting Works in Quarries and Their Impact on the Environment and Population. Appl. Sci. 2021, 11, 2118. https://doi.org/10.3390/app11052118. 4. Sanayei, M.; Pradeep, M.; Moore, J. Measurement of building foundation and ground-borne vibrations due to surface trains and subways. Eng. Struct. 2013, 53, 102–111. https://doi.org/10.1016/j.engstruct.2013.03.038. 5. Czech, K.; Gosk, W. Analysis of the Vibration Propagation Induced by Pulling out of Sheet Pile Wall in a Close Neighbourhood of Existing Buildings. Procedia Eng. 2013, 143, 1460–1467. https://doi.org/10.1016/j.proeng.2016.06.172. 6. Kepak, S.; Stolarik, M.; Nedoma, J.; Martinek, R.; Kolarik, J.; Pinka, M. Alternative Approaches to Vibration Measurement Due to the Blasting Operation: A Pilot Study. Sensors 2019, 19, 4084. https://doi.org/10.3390/s19194084. 7. Nedoma, J.; Stolarik, M.; Fajkus, M.; Pinka, M.; Hejduk, S. Use of Fiber-Optic Sensors for the Detection of the Rail Vehicles and Monitoring of the Rock Mass Dynamic Response Due to Railway Rolling Stock for the Civil Engineering Needs. Appl. Sci. 2019, 9, 134. https://doi.org/10.3390/app9010134. 8. Nedoma, J.; Stolarik, M.; Kepak, S.; Pinka, M.; Martinek, R.; Frnda, J.; Fridrich, M. Alternative Approaches to Measurement of Ground Vibrations Due to the Vibratory Roller: A Pilot Study. Sensors 2019, 19, 5420. https://doi.org/10.3390/s19245420. 9. Stolarik, M.; Pinka, M.; Zabka, S.; Novak, M. Dynamic effect of harmonic vibrations: Various approaches to monitor. In ElectroOptical and Infrared Systems: Technology and Applications XV; SPIE: Bellingham, WA, USA, 2018; Volume 107950. https://doi.org/1117/12.2325740. 10. Stolarik, M.; Pinka, M.; Nedoma, J. Ground-Borne Vibration Due to Construction Works with Respect to Brownfield Areas. Appl. Sci. 2019, 9, 3766. https://doi.org/10.3390/app9183766. 11. Wentao, Z.; Li, H.; Huang, W.; Li, Z.; Li, L.; Liu, R. Research Progress of Optical Fiber Seismograph. J. Appl. Sci. 2021, 39, 821– 842. https://doi.org/10.3969/j.issn.0255-8297.2021.05.008. 12. Okamoto, T.; Iida, D.; Koshikiya, Y.; Honda, N. Deployment Condition Visualization of Aerial Optical Fiber Cable by Distributed Vibration Sensing Based on Optical Frequency Domain Reflectometry. J. Lightwave Technol. 2021, 39, 6942–6951. https://doi.org/10.1109/JLT.2021.3107855. 13. He, Z.; Liu, Q. Principles and Applications of Optical Fiber Distributed Acoustic Sensors. Laser Optoelectron. Prog. 2021, 58, 1306001. https://doi.org/10.3788/LOP202158.1306001. Figure A1. ( a – e ) Measuring situations 1–5 with distances between gravel piles and measuring stations marked. References 1. Stolarik, M.; Nedoma, J.; Martinek, R.; Kepak, S.; Hrubesova, E.; Pinka, M.; Kolarik, J. New Methods to Seismic Monitoring: Laboratory Comparative Study of Michelson Fiber-Optic Interferometer and Pneumatic Measurement Systems. Photonics 2021 , 8, 147. [CrossRef] 2. Zhou, X.; Zhang, X.; Feng, H.; Zhang, S.; Yang, J.; Mu, J.; Hu, T. Study on Dominant Frequency Attenuation of Blasting Vibration for Ultra-Small-Spacing Tunnel. Appl. Sci. 2022,12, 1058. [CrossRef] 3. Feher, J.; Cambal, J.; Pandula, B.; Kondela, J.; Sofranko, M.; Mudarri, T.; Buchla, I. Research of the Technical Seismicity Due to Blasting Works in Quarries and Their Impact on the Environment and Population. Appl. Sci. 2021,11, 2118. [CrossRef] 4. Sanayei, M.; Pradeep, M.; Moore, J. Measurement of building foundation and ground-borne vibrations due to surface trains and subways. Eng. Struct. 2013,53, 102–111. [CrossRef] 5. Czech, K.; Gosk, W. Analysis of the Vibration Propagation Induced by Pulling out of Sheet Pile Wall in a Close Neighbourhood of Existing Buildings. Procedia Eng. 2013,143, 1460–1467. [CrossRef] 6. Kepak, S.; Stolarik, M.; Nedoma, J.; Martinek, R.; Kolarik, J.; Pinka, M. Alternative Approaches to Vibration Measurement Due to the Blasting Operation: A Pilot Study. Sensors 2019,19, 4084. [CrossRef] 7. Nedoma, J.; Stolarik, M.; Fajkus, M.; Pinka, M.; Hejduk, S. Use of Fiber-Optic Sensors for the Detection of the Rail Vehicles and Monitoring of the Rock Mass Dynamic Response Due to Railway Rolling Stock for the Civil Engineering Needs. Appl. Sci. 2019 , 9, 134. [CrossRef] 8. Nedoma, J.; Stolarik, M.; Kepak, S.; Pinka, M.; Martinek, R.; Frnda, J.; Fridrich, M. Alternative Approaches to Measurement of Ground Vibrations Due to the Vibratory Roller: A Pilot Study. Sensors 2019,19, 5420. [CrossRef] 9. Stolarik, M.; Pinka, M.; Zabka, S.; Novak, M. Dynamic effect of harmonic vibrations: Various approaches to monitor. In Electro-Optical and Infrared Systems: Technology and Applications XV; SPIE: Bellingham, WA, USA, 2018; Volume 107950. [CrossRef] 10. Stolarik, M.; Pinka, M.; Nedoma, J. Ground-Borne Vibration Due to Construction Works with Respect to Brownfield Areas. Appl. Sci. 2019,9, 3766. [CrossRef] 11. Wentao, Z.; Li, H.; Huang, W.; Li, Z.; Li, L.; Liu, R. Research Progress of Optical Fiber Seismograph. J. Appl. Sci. 2021 ,39, 821–842. [CrossRef] 12. Okamoto, T.; Iida, D.; Koshikiya, Y.; Honda, N. Deployment Condition Visualization of Aerial Optical Fiber Cable by Distributed Vibration Sensing Based on Optical Frequency Domain Reflectometry. J. Lightwave Technol. 2021,39, 6942–6951. [CrossRef] 13. He, Z.; Liu, Q. Principles and Applications of Optical Fiber Distributed Acoustic Sensors. Laser Optoelectron. Prog. 2021 ,58, 1306001. [CrossRef] 14. Han, B.; Guan, H.; Yao, J.; Rao, Y.-J.; Ran, Z.; Gong, Y.; Li, Q.; Li, M.; Zhang, R.; An, S.; et al. Distributed Acoustic Sensing with Sensitivity-Enhanced Optical Cable. IEEE Sens. J. 2021,21, 4644–4651. [CrossRef] Sensors 2022,22, 5579 19 of 20 15. Becker, M.W.; Coleman, T.I. Distributed Acoustic Sensing of Strain at Earth Tide Frequencies. Sensors 2019,19, 1975. [CrossRef] 16. Lindsey, N.; Martin, E.R.; Dreger, D.S.; Freifeld, B.; Cole, S.; James, S.R.; Biondi, B.L.; Ajo-Franklin, J.B. Fiber-optic network observations of earthquake wavefields: Fiber-optic earthquake observations. Geophys. Res. Lett. 2017 ,44, 11792–11799. [CrossRef] 17. Zhou, R.; Chen, F.; Li, S.; Wang, R.; Qiao, X. Three-Dimensional Vector Accelerometer Using a Multicore Fiber Inscribed with Three FBGs. J. Lightwave Technol. 2021,39, 3244–3250. [CrossRef] 18. Zhang, W.; Huang, W.; Li, L.; Liu, W.; Li, F. High resolution strain sensor for earthquake precursor observation and earthquake monitoring. In Sixth European Workshop on Optical Fibre Sensors; SPIE: Bellingham, WA, USA, 2016; Volume 9916, pp. 58–61. [CrossRef] 19. Zhang, W.; Huang, W.; Li, L.; Liu, W.; Li, F. High resolution FBG sensor and its applications in Geophysics. In Proceedings of the 2017 16th International Conference on Optical Communications and Networks (ICOCN), Wuzhen, China, 7–10 August 2017; pp. 1–3. [CrossRef] 20. Huang, W.; Zhang, W.; Luo, Y.; Li, L.; Liu, W.; Li, F. Broadband FBG resonator seismometer: Principle, key technique, self-noise, and seismic response analysis. Opt. Express 2018,26, 10705–10715. [CrossRef] 21. Zhang, W.; Huang, W. Applications of Fiber Optics Sensors in Seismology. In Proceedings of the 2018 10th International Conference on Advanced Infocomm Technology, Stockholm, Sweden, 12–15 August 2018; pp. 16–20. [CrossRef] 22. Wang, J.; Hu, B.; Li, W.; Song, G.; Jiang, L.; Liu, T. Design and application of fiber Bragg grating (FBG) geophone for higher sensitivity and wider frequency range. Measurement 2016,79, 228–235. [CrossRef] 23. Saad, S. FBG sensors for seismic control and detection in extradosed bridges. Int. J. Smart Sens. Intell. Syst. 2019 ,14, 1–13. [CrossRef] 24. Fallais, D.J.M.; Henkel, M.; Noppe, N.; Weijtjens, W.; Devriendt, C. Multilevel RTN Removal Tools for Dynamic FBG Strain Measurements Corrupted by Peak-Splitting Artefacts. Sensors 2021,22, 92. [CrossRef] 25. Perez-Alonzo, A.; Sandoval-Romero, G.E. Accelerometer prototype based on enhanced fiber Bragg grating overlapping interrogation method. Optik 2021,242, 167027. [CrossRef] 26. Han, J.; Zhang, W.; Jiang, D.; Wang, Z.; Li, F. Fiber optic 3-component seismometer. Photonic Sens. 2014,4, 102–107. [CrossRef] 27. Chen, J.; Chang, T.; Fu, Q.; Lang, J.; Gao, W.; Wang, Z.; Yu, M.; Zhang, Y.; Cui, H.-L. A Fiber-Optic Interferometric Tri-Component Geophone for Ocean Floor Seismic Monitoring. Sensors 2017,17, 47. [CrossRef] [PubMed] 28. Jena, J.; Wassin, S.; Bezuidenhout, L.; Doucouré, M.; Gibbon, T. Polarization-based optical fiber acoustic sensor for geological applications. J. Opt. Soc. Am. 2020,37, A147–A153. [CrossRef] 29. Chen, J.; Yang, Y.; Gao, W.; Chang, T.; Cui, H. High Sensitivity Optical Fiber Interferometric Accelerometer for Seismic Observation. In Proceedings of the 26th International Conference on Optical Fiber Sensors, Lausanne, Switzerland, 24–28 September 2018; p. ThE103. [CrossRef] 30. Amorebieta, J.; Ortega-Gomez, A.; Durana, G.; Fernández, R.; Antonio-Lopez, E.; Schülzgen, A.; Zubia, J.; Amezcua-Correa, R.; Villatoro, J. Compact omnidirectional multicore fiber-based vector bending sensor. Sci. Rep. 2021,11, 5989. [CrossRef] 31. Azmi, A.; Abdullah, A.; Noor, M.M.; Ibrahim, M.; Ibrahim, R.R.; Tan, T.; Zhang, J. Dynamic bending and rotation sensing based on high coherence interferometry in multicore fiber. Opt. Laser Technol. 2021,135, 106716. [CrossRef] 32. Amorebieta, J.; Ortega-Gomez, A.; Durana, G.; Fernández, R.; Antonio-Lopez, E.; Schülzgen, A.; Zubia, J.; Amezcua-Correa, R.; Villatoro, J. Highly sensitive multicore fiber accelerometer for low frequency vibration sensing. Sci. Rep. 2020 ,10, 16180. [CrossRef] 33. Villatoro, J.; Antonio-Lopez, E.; Schülzgen, A.; Amezcua-Correa, R. Miniature multicore optical fiber vibration sensor. Opt. Lett. 2017,42, 2022–2025. [CrossRef] 34. Yi, D.; Liu, F.; Zhang, M.; He, X.; Zhou, X.; Long, K.; Li, X. Demonstration of Fiber-Optic Seismic Sensor with Improved Dynamic Response in Oilfield Application. IEEE Trans. Instrum. Meas. 2022,71, 1–8. [CrossRef] 35. Ponceau, D.; Millier, P.; Olivier, S. Subnanometric Michelson interferometry for seismological applications. Opt. Sens. 2008 ,7003, 251–260. [CrossRef] 36. Chang, T.; Wang, Z.; Yang, Y.; Zhang, Y.; Zheng, Z.; Cheng, L.; Cui, H.-L. Fiber optic interferometric seismometer with phase feedback control. Opt. Express 2020,28, 6102–6122. [CrossRef] [PubMed] 37. Huang, J.; Zhang, W.; Huang, W.; Huang, W.; Wang, L.; Luo, Y. High-resolution fiber optic seismic sensor array for intrusion detection of subway tunnel. In Proceedings of the Asia Communications and Photonics Conference ACP, Hangzhou, China, 26–29 October 2018. [CrossRef] 38. Yang, Y.; Wang, Z.; Chang, T.; Yu, M.; Chen, J.; Zheng, G.; Cui, H. Seismic observation and analysis based on three-component fiber optic seismometer. IEEE Access 2020,8, 1374–1382. [CrossRef] 39. Jin-Hung, H.; Tung-Yueh, T. Ground Vibration during Gravel Pile Construction. J. Mar. Sci. Technol. 2002,10, 6. [CrossRef] 40. Lv, G.; Cui, W.; Zhang, Q.; Wu, S.; Wang, S. Experimental Study on Embankment Failure Law and Reinforcement Technology for Highway Widening Project over Silt Soils. J. Test. Eval. 2022,50, 1–13. [CrossRef] 41. Li, Q.; Wang, Z.; Yu, Y. Optimization of Area Replacement Ratio of Stone Column Based on Robustness Principle. J. Jilin Univ. Earth Sci. Ed. 2022,52, 171–180. 42. Xiong, S.; Wu, G.; Huang, S.; Ming, W. Method and Parameter Reliability of In-situ Direct Shear Test on Vibroflotation Gravel Piles. IOP Conf. Ser. Earth Environ. Sci. 2021,2021, 861. [CrossRef] Sensors 2022,22, 5579 20 of 20 43. Guo, B.; Yang, Y.; Chen, H.; Du, F. Nonlinear consolidation solution for composite foundation reinforced by stone columns. J. Xi’an Univ. Archit. Technol. 2009,53, 652–658. [CrossRef] 44. Brož, M.; Štrunc, J.; Levý, O.; Haislová, R. Using Multi-Purpose Instruments Brs 32 for Environment Attributes Determination on Seismic Profiles. Explor. Geophys. Remote Sens. Environ. 2018,25, 51–60. [CrossRef] 45. Brož, M.; Štrunc, J. Long-Term On-Line Expedition Monitoring of Technical Seismicity. Explor. Geophys. Remote Sens. Environ. 2019,26, 17–25. [CrossRef] 46. Molinek, O.; Vokurka, M.; Mokrosova, A.; Kotrbanec, J. Utilization of field seismometers for measuring the paraseismic vibrations of mining blasts. IOP Conf. Ser. Earth Environ. Sci. 2021,684, 012003. [CrossRef] 47. Yin, S.; Ruffin, P.B.; Yu, F.T.S. (Eds.) Fiber Optic Sensors, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2008. [CrossRef] 48. Liu, X.; Jin, B.; Bai, Q.; Wang, Y.; Wang, D.; Wang, Y. Distributed Fiber-Optic Sensors for Vibration Detection. Sensors 2016 ,16, 1164. [CrossRef] [PubMed] 49. Zhang, W.; Lu, P.; Qu, Z.; Zhang, J.; Wu, Q.; Liu, D. Passive Homodyne Phase Demodulation Technique Based on LF-TIT-DCM Algorithm for Interferometric Sensors. Sensors 2021,21, 8257. [CrossRef] [PubMed] 50. Todd, M.; Seaver, M.; Bucholtz, F. Improved, operationally-passive interferometric demodulation method using 3 × 3 coupler. Electron. Lett. 2002,38, 784–786. [CrossRef]