Neural Network-Based Laser Interferometer Compensation for Seismic Signal Detection.
An earthquake is a devastating but common natural disaster around the world. To reduce and predict the damage caused by earthquakes, it is necessary to precisely measure and analyze the seismic signal. Hence, the study of earthquakes has been performed through various seismic analyses and predictions. Methods to analyze earthquakes include spectrum analysis, magnitude measurement, intensity interpretation, and epicenter decision. The short-time average/long-time average (STA/LTA) trigger algorithm is commonly applied in microseismic monitoring systems and has gained significant increased interest. In the STA/LTA algorithm, the most important factor for seismic wave analysis is to identify a seismic event that includes an external disturbance. Widely used seismic measurement devices include seismic velocimeters and accelerometers. These apparatuses detect seismic waves based on force and pressure sensors. They can be classified according to their measurable frequency range as short-period, long-period, and broadband seismographs. These apparatuses have defects such as limited frequency range and low resolution. Therefore, errors exist when calculating the exact P- and S-wave arrival time.
To improve seismic analyses, Colak et al.  proposed a new P- and S-wave arrival time decision method. The discrete wavelet transform was used to analyze the change of frequency and energy in high frequency bands. As a result, P- and S-wave arrival time can be measured using a series of data windows. Hloupis and Vallianatos  suggested an earthquake early warning system, based on a redundant wavelet transform as a processing tool to rapidly estimate earthquake magnitude. The suggested algorithm strengthens the ability to distinguish body waves from noise. Panagiotakis et al.  proposed an automatic identification method of P-phase arrival based on energy and frequency characteristics of the local maxima distribution. Moreover, the method performs a robust detection against the background noise of seismograms. Baziw  proposed a robust seismic event identification algorithm by obtaining the optimal estimate of a possible seismic event based on the Rao-Blackwellized particle filter and hidden Markov model. This enables obtaining the real-time phase estimate and the quantified frequency content of seismic data.
In this paper, we propose an interferometric seismometer that is based on a heterodyne laser interferometer and analyze seismic data in the frequency range to estimate epicenter location. Using the orthogonal characteristic of two beams from a laser head, the seismic signals are expressed as a form of intensity signal, which is represented as a trigonometrical function. Hence, the seismic signal can be detected using the laser interferometer. Each observatory has a different experimental environment since the observatory needs to be located in various regions for the seismic signal measurement. The nonuniformity of measurement environment such as temperature and humidity causes the measurement inaccuracy. To compensate for the environmental error that resides in seismic measurement data using a data fusion method, we measure the displacement with two independent laser heads set separately in opposite directions. Using the neural network, the standard deviation of errors is updated adaptively according to the experiment temperature change. With the measured data sets, the environmental error is reduced by projecting the measured data onto the points on the constraint set based on the ideal relationship between the two interferometers [5, 6]. The compensated data are converted into a signal in the frequency domain to extract the P- and S-wave arrival times, which are used to calculate the STA/LTA ratio. The difference between the two arrival times leads to the determination of the epicenter location.
2. A Heterodyne Laser Interferometric Seismometer
To measure displacement, a heterodyne laser interferometer uses the phase difference caused by the Doppler effect. The heterodyne laser interferometer system consists of a He-Ne laser head, a beam splitter, a polarization beam splitter, photodetectors, a moving mirror, and a fixed mirror. Figure 1 shows the heterodyne laser interferometer seismic system. The electric fields are measured by detector X. Likewise, the electric fields reflected by the fixed mirror and the moving mirror are measured by detector Y. The intensities obtained by detectors X and Y are expressed as follows.
[I.sub.r] [infinity] 1/2 ([A.sup.2] + [B.sup.2]) + AB cos [[DELTA][omega]t + ([[phi].sub.B] - [[phi].sub.A])],
[I.sub.r] [infinity] 1/2 ([A.sup.2] + [B.sup.2]) + AB cos [[DELTA][omega]t + ([[phi].sub.B] - [[phi].sub.A]) + [DELTA][phi]]]. (1)
In (1), A and B are the magnitudes of the electric fields and [[phi].sub.A] [[phi].sub.B] are the initial phases of laser sources. The signal [I.sub.r] measured by detector X is the reference intensity, and the signal [I.sub.m] obtained by detector Y is the measurement intensity. [DELTA][omega] is the difference between the frequencies of the two orthogonal beams [[omega].sub.1] and [[omega].sub.2], and [DELTA][phi] is the phase difference caused by the Doppler effect. We use a high-pass filter to eliminate DC components of the intensity signal [I.sub.m]. The initial phase values [[phi].sub.A] and [[phi].sub.B] are ignored to obtain phase information. We obtain the AC signal of [I.sub.m] as follows.
[I.sub.m,AC] [infinity] AB cos([DELTA][omega]t + [DELTA][phi]). (2)
To obtain the phase value, we use a lock-in amplifier. The phase value is derived by applying the [I.sub.r] and [I.sub.m,AC] signals into a power divider, multiplier, and low pass filter, sequentially. As a result, two orthogonal intensity signals of [I.sub.x] [infinity] [I.sub.m,AC] * [I.sub.r] and [I.sub.y] [infinity] [I.sub.m,AC] * [I.sub.r][e.sup.j[pi]/2] can be represented as a trigonometric function.
[I.sub.x] [infinity] AB/2 cos ([DELTA][phi]),
[I.sub.y] [infinity] AB/2 sin ([DELTA][phi]). (3)
By using an inverse trigonometric function, the phase information [DELTA][phi] related to the displacement is obtained. We calculate the displacement change using [DELTA][phi] = 4[pi][rho][DELTA]L/[lambda]. [rho] is the refractive index, [lambda] is the average wavelength, and [DELTA]L is the difference in optical path length between the fixed mirror and the moving mirror.
3. Data Fusion-Based Error Compensation
The environmental error causes a significant problem of a measurement accuracy in case of the nanoscale instrument . In order to improve measurement precision, we use a data fusion method to reduce the error that exists in the measurement data. Data fusion uses two measurement data sets, which in this study are obtained by two laser systems installed separately in opposite directions. To utilize this measured data as a seismic signal, we compensate the error to improve the interferometer's accuracy .
We use two independent laser heads, which have the same performance characteristics, for the data fusion method. The two separate laser interferometers measure the same moving mirror with a displacement of [[??].sub.1] and [[??].sub.2], respectively. In the ideal case, the displacement [d.sub.1] and [d.sub.2] should satisfy the following constraint, because the two laser interferometers measure the identical moving mirror in the opposite direction. When we measure the displacement, the data includes environmental errors since environmental conditions are not uniform. Therefore, the measurement values ([[??].sub.1], [[??].sub.2]) do not satisfy the ideal constraint of [d.sub.1](k) + [d.sub.2](k) = 0.
We can represent the relationship between [[??].sub.1] and [[??].sub.2] in the Cartesian coordinate, where [[??].sub.1] and [[??].sub.2] are distributed around the constraint of [d.sub.1](k) + [d.sub.2](k) = 0. These patterns show that the measured displacement data contains environmental errors. We project the measured points in the Cartesian coordinate onto the points on the ideal constraint line to compensate for the environmental errors. An objective error function is utilized to minimize the environmental error.
[J.sub.e] = 1/2 [N.summation over (k=1)][parallel][([I.sub.x] (k) - [I*.sub.x] (k)).sup.2] + [([I.sub.y](k) - [I*.sub.y](k)).sup.2][parallel]. (4)
Here, [I.sub.x] and [I.sub.y] are measurement signals that contain environmental error, while [I*.sub.x] and [I*.sub.y] are reference signals without environmental error.
Even though the standard deviation is calculated for each measurement step, a compensation parameter ([zeta]) represents one of the values derived under the fixed environment. To improve measurement accuracy, it is necessary to adaptively modify [zeta] by considering temperature change, which is one of the primary environmental error factors. To estimate the adaptive standard deviations, a multilayer neural network is used. The network is composed of one input layer, two hidden layers, and one output layer. The neurons in the hidden layers and the output layer follow the hyperbolic tangent-sigmoid function [f.sub.1](x) and the log-sigmoid function [f.sub.2](x) as an activation function, respectively. The activation functions enable the neural network to obtain a complex nonlinear model.
[f.sub.1](x) = [e.sup.x] - [e.sup.-x] / [e.sup.x] + [e.sup.-x],
[f.sub.2] (x) = 1/1 + [e.sup.-x]. (5)
The backpropagation algorithm is applied to train and update weights in the network. The steepest gradient descent method is used during the training period. As a target vector, we obtain the error between the reference signals and the measurement signals.
[mathematical expression not reproducible]. (6)
Here, [I.sub.tp] (k) = [tau] is the temperature value during the measurement process. The intensities from (3) are used as the inputs of the neural network during the training period. We define [V.sub.I] = [[[I.sub.x], [I.sub.y], [I.sub.tp]].sup.T] as an input vector and [mathematical expression not reproducible] as an output vector. The objective function uses the errors between the output from the neural network and the target vector, [V.sub.T]:
where m is the amount of sampled data. Weight values are calculated with a steepest gradient descent method and updated until minimizing the objective function.
[w.sup.k+1.sub.j] = [w.sup.k.sub.j] - [eta] [partial derivative]E/[partial derivative]w. (8)
Here, [eta] is a learning rate and -[partial derivative]E/[partial derivative]w is calculated using the chain rule. To guarantee a reliable performance, the error bound of the objective function is set as [10.sup.-9]. The maximum number of the training iterations is limited to 10,000 times for computational efficiency. Finally, we get the neural network modeled in Figure 2.
The estimated error vector [mathematical expression not reproducible] associated with temperature is obtained from the neural network by applying measurement signals ([I.sub.x], [I.sub.y]) and temperature ([I.sub.tp]) as an input vector.
[mathematical expression not reproducible]. (9)
The adaptive standard deviations ([[sigma]*.sub.1], [[sigma]*.sub.2]) are calculated for each measurement step with the neural network output, [mathematical expression not reproducible]. The loci of the ellipsoidal equation with a centroid, ([[??].sub.1] (k), [[??].sub.2] (k)), and the contact points between [x.sub.m2] + [x.sub.m1] = 0 and each ellipse can be represented as follows:
[mathematical expression not reproducible]. (10)
Here, [[xi].sub.k] is the adaptive standard deviation ratio, [[xi].sub.k] = [[sigma]*.sub.2] (k) / [[sigma]*.sub.1](k). Finally, the environmental error is com pensated by using the adaptive data fusion based on an artificial intelligence projection method.
4. Seismic Data Analysis and STA/
To detect the seismic wave, we use the heterodyne laser interferometer system. Using the above equations, we eliminate the environmental error by the data fusion method, enabling more precise measurement of the seismic wave. With the seismic data, we obtain P-S wave arrival time (P-S time), which is an important factor for calculating the epicenter location. The methods that measure P-S time generally analyze seismic data in the time domain. However, if the amplitude of a seismic wave does not change enough during the measurement interval, it is difficult to distinguish between P-wave and S-wave arrival . Hence, we suggest a P-S time measurement method in the frequency domain to more accurately identify the arrival times of P-wave and S-waves.
When we detect the seismic wave with a laser interferometer, the seismic signal can be expressed in terms of the intensity of the laser signal in a trigonometric function. The equation can be represented as follows:
[mathematical expression not reproducible]. (11)
In (11), [[??].sub.x] and [[??].sub.y] are the intensity signals after applying the data fusion method to [I.sub.x] and [I.sub.y] in (3). As the refractive index n and wavelength [lambda] are constants, the change in seismic wave movement, [DELTA]L, is proportional to the phase change, [DELTA][??]. Therefore, the frequency data of the phase signal [DELTA][??] is obtained by the polarity change of the seismic signal, and the [I.sub.x] and [I.sub.y] signals are both suitable for frequency analysis. However, the [I.sub.y] signal is used to analyze the seismic wave as a sinusoidal wave to avoid a latent effect of nonlinearity . The signals are measured in a discrete form. The frequency can be calculated from the signal's polarity change with the following equation.
Sgn[[[??].sub.y] (k)] * sgn [[[??].sub.y] (k +1)] <0. (12)
Here, sgn(*) is the signum function. The final frequency value is represented as [f.sub.c](k) = [N.sub.c](k)/2. [f.sub.c](k) is the frequency of the incoming intensity signal, and [N.sub.c](k) is the number of [[??].sub.y] (k) sign changes in one second. Using this method, we can measure the seismic signal in the frequency domain.
An earthquake's epicenter can be determined using the frequency variation of [[??].sub.y] (k) measured by the heterodyne laser interferometer. We analyze this data with the STA/ LTA ratio to calculate P-S time. With a threshold trigger algorithm, we can enhance the analyzing accuracy of seismic data. The STA/LTA algorithm parameters, such as short-time window size, long-time window size, and threshold value, are determined by the seismic system environment. The STA/LTA algorithm analyzes the data using the difference between STA and LTA characteristics. Finally, we compare the measured STA/LTA ratio with the threshold value. The accuracy of the analysis data depends on the STA and LTA window sizes and threshold value. The STA and LTA are represented as follows [11-13].
[mathematical expression not reproducible]. (13)
Here, [[tau].sub.s] and [[tau].sub.1] are the size of the short- and long-time windows, respectively. We measure the time when the STA/LTA ratio becomes greater than the threshold value. The distance between a seismological observatory and an epicenter is calculated using measured P-S time by the following equation.
D = [V.sub.P] [V.sub.S]/[V.sub.P]-[V.sub.S] [T.sub.ps]. (14)
Here, D is the epicentral distance, [V.sub.P] is the P-wave velocity, [V.sub.S] is the S-wave velocity, and [T.sub.ps] is the P-S time.
5. Performance Evaluation
In this section, we show the performance of the proposed seismic detection system based on a heterodyne laser interferometer. To prove the performance and effectiveness of the interferometric seismometer, we use measured seismic data detected by the laser interferometric seismometer and the accelerometer, respectively. We use the He-Ne laser head with an intermode beat of 632.9912 nm (Wavetronics: WT307B). The laser head emits two-frequency beams of [[lambda].sub.1] = 632.9912576nm for vertical polarization and [[lambda].sub.2] = 632.9912604 nm for horizontal polarization, respectively. The linear stage (PI: M531.5IM) is used as a vibrational source to generate a seismic signal. An accelerometer (Mitutoyo: JEP-8A3) is used as a reference. Figure 3 shows the seismic wave measurement system.
First, we prove the environmental error compensation by data fusion method for the measurement precision enhancement in a laser interferometer. In our simulation, we set the wavelength [lambda] as 632.9912 nm and the refractive index as 1.00000026653516. Figure 4 shows the compensated data by the data fusion algorithm. The thick line and thin line represent the compensated seismic data and the measurement seismic data, respectively. To verify the accuracy of compensated seismic data using the data fusion in Figure 4, we compare the root-mean-square (RMS) of the compensated measurement error with the uncompensated one. The RMS error comparison results appear in Figure 5. The environmental errors that occurred in the seismic signal detection are compensated by data fusion. Figure 6 shows the seismic data measured by the laser interferometric seismometer and the accelerometer, respectively. Figure 7 denotes the difference between the true amplitude and the measurement value at each sampling time. The thick solid line represents the RMS error of the compensated interferometric seismometer amplitude data. The thin dotted line means the RMS error measured from an accelerometer. As shown in this figure, the proposed neural network-based compensated seismic wave is more precise than the measurement of accelerometer.
Figure 8 shows the intensity signal [I.sub.y](k) of the seismic wave and the signal in a time-frequency plane measured by the interferometric seismometer. We obtain the intensity signal ([I.sub.y]) to minimize the nonlinearity error. The displacement change of the seismic wave is measured from the phase difference between [I.sub.r] and [I.sub.m] of the heterodyne interferometer. To analyze the seismic data, we represent the seismic signal in a time-frequency plane and observe sudden frequency changes around 6 and 12 seconds, which are generated by P-wave and S-wave, respectively.
We use the STA/LTA algorithm to measure P-S time for calculating the epicenter location. First, we determine the STA and LTA window size. The LTA window size is set as 14 times the STA. Next, the threshold value used for the P-wave and S-wave arrival time is set as 1.3 times the average of the STA/LTA ratio. The STA/LTA ratio and the threshold value are shown in Figure 9.
From the STA/LTA ratio in Figure 8, we measure P-wave and S-wave arrival time as 5.66 and 12.71 seconds, respectively. With the P-wave velocity and the S-wave velocity as 8 km/s and 4 km/s, respectively, we calculate the epicentral distance as 24.20 km using (14). Figure 10 shows the epicenter location. [n.sub.e] means the maximum absolute value of the environmental noises. We suppose the parameter [n.sub.e] as 0.05 m/[s.sup.2] in Figure 10. The distance from the observatory to the epicenter is measured by the STA/LTA algorithm.
With three distinct sets of observatory data, we can estimate the epicenter location in the overlapping region. Table 1 shows the RMS error comparison of epicenter localization method based on the laser interferometric seismometer and accelerometer. Although the localization accuracy is decreased as the parameter ne increases, we can confirm that the localization result based on the laser interferometric seismometer shows better performance than accelerometer.
In this paper, we proposed a seismic signal detection method using a heterodyne laser interferometer in the frequency domain. We applied the data fusion algorithm to minimize the objective error function. To prove the performance of the interferometric seismometer, we compared seismic data measured by a laser interferometric seismometer with data detected by an accelerometer. Moreover, we measured the seismic signal data in the frequency range with the laser interferometer. Using the STA/LTA algorithm, the P-S time was derived. The epicentral distance was calculated using the measured P-S time. We confirmed the effectiveness of the laser interferometer-based seismometer through simulation and experiment.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper.
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (NRF-2016R1A2B4011369).
 O. Colak, T. Destici, S. Ozen, H. Arman, and O. Cerezci, "Detection of P- and S- wave arrival times using the discrete wavelet transform in real seismograms," The Arabian Journal for Science and Engineering, vol. 34, no. 1A, pp. 79-89, 2009.
 G. Hloupis and F. Vallianatos, "Wavelet-based rapid estimation of earthquake magnitude oriented to early warning," IEEE Geoscience and Remote Sensing Letters, vol. 10, no. 1, pp. 43-47, 2013.
 C. Panagiotakis, E. Kokinou, F. Vallianatos, and P. Automatic, "Automatic P-phase picking based on local-maxima distribution," IEEE Transactions on Geoscience and Remote Sensing, vol. 46, no. 8, pp. 2280-2287, 2008.
 E. Baziw, "Real-time seismic signal enhancement utilizing a hybrid Rao-Blackwellized particle filter and hidden Markov model filter," IEEE Geoscience and Remote Sensing Letters, vol. 2, no. 4, pp. 418-422, 2005.
 A. Effler, R. M. S. Schofield, V. V. Frolov et al., "Environmental influences on the LIGO gravitational wave detectors during the 6th science run," Classical and Quantum Gravity, vol. 32, no. 3, article 035017, 2015.
 L. Yan, B. Chen, E. Zhang, S. Zhang, and Y. Yang, "Precision measurement of refractive index of air based on laser synthetic wavelength interferometry with Edlen equation estimation," Review of Scientific Instruments, vol. 86, no. 8, article 085111, 2015.
 S. G. Gevorgyan, V. S. Gevorgyan, H. G. Shirinyan, G. H. Karapetyan, and A. G. Sarkisyan, "A radically new principle of operation seismic detector of nano-scale vibrations," IEEE Transactions on Applied Superconductivity, vol. 17, no. 2, pp. 629-632, 2007.
 W. Lee, J. Lee, and K. You, "Improved measurement accuracy in heterodyne laser interferometer using WDF," Electronics Letters, vol. 45, no. 21, pp. 1085-1087, 2009.
 S. Gabarda and G. Cristobal, "Detection of events in seismic time series by time-frequency methods," IET Signal Processing, vol. 4, no. 4, pp. 413-420, 2010.
 W. Hou, "Optical parts and the nonlinearity in heterodyne interferometers," Precision Engineering, vol. 30, no. 3, pp. 337-346, 2006.
 B. Sharma, A. Kumar, and V. Murthy, "Evaluation of seismic events detection algorithms," Journal of the Geological Society of India, vol. 75, no. 3, pp. 533-538, 2010.
 A. Velasco, S. Hernandez, T. Parsons, and K. Pankow, "Global ubiquity of dynamic earthquake triggering," Nature Geoscience, vol. 1, no. 6, pp. 375-379, 2008.
 A. Maggi, C. Tape, M. Chen, D. Chao, and J. Tromp, "An automated time-window selection algorithm for seismic tomography," Geophysical Journal International, vol. 178, no. 1, pp. 257-281, 2009.
Kyunghyun Lee, Hyungkwan Kwon, and Kwanho You
Department of Electrical Engineering, Sungkyunkwan University, Suwon, Republic of Korea
Correspondence should be addressed to Kwanho You; email@example.com
Received 19 July 2017; Revised 18 September 2017; Accepted 4 October 2017; Published 22 February 2018
Academic Editor: Jesus Corres
Caption: FIGURE 1: Seismic measurement system based on heterodyne laser interferometer.
Caption: FIGURE 2: Neural network model to compensate for the environmental error.
Caption: FIGURE 3: Laser interferometric seismometer.
Caption: FIGURE 4: Seismic data compensation with data fusion.
Caption: FIGURE 5: RMS of displacement error.
Caption: FIGURE 6: Seismic data detection by interferometric seismometer and accelerometer.
Caption: FIGURE 7: RMS error comparison between a laser interferometer and accelerometer.
Caption: FIGURE 8: Intensity and frequency variation of the seismic signal.
Caption: FIGURE 9: STA/LTA ratio.
Caption: FIGURE 10: Epicenter location.
TABLE 1: RMS error comparison for various environmental errors ([n.sub.e]). RMS error (km) Measurement [n.sub.e] = 0.05 [n.sub.e] = 0.07 Laser interferometric 2.54 2.81 seismometer Accelerometer 4.18 5.22 RMS error (km) Measurement [n.sub.e] = 01 [n.sub.e] = 0.15 Laser interferometric 4.27 5.83 seismometer Accelerometer 7.19 9.45 RMS error (km) Measurement [n.sub.e] = 0.2 Laser interferometric 7.19 seismometer Accelerometer 13.21
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Lee, Kyunghyun; Kwon, Hyungkwan; You, Kwanho|
|Publication:||Journal of Sensors|
|Date:||Jan 1, 2018|
|Previous Article:||Research on Real-Time Supervisory System for Compaction Quality in Face Rockfill Dam Engineering.|
|Next Article:||Secure Data Encryption for Cloud-Based Human Care Services.|