Research Article  Open Access
Erlei Yao, Yu Miao, Guobo Wang, "Synthesis of Spatially Correlated Earthquake Ground Motions Based on Hilbert Transform", Modelling and Simulation in Engineering, vol. 2017, Article ID 2614769, 9 pages, 2017. https://doi.org/10.1155/2017/2614769
Synthesis of Spatially Correlated Earthquake Ground Motions Based on Hilbert Transform
Abstract
A simplified method for synthesizing spatially correlated earthquake ground motions is developed based on Hilbert transform and a reference earthquake record. In this method, one reference earthquake record is treated as the original ground motion, based on a series of generated ground motions. This procedure uses the instantaneous amplitude and the instantaneous phase of the record obtained using Hilbert transform to achieve the nonstationarity of ground motion. To establish the coherency between generated ground motions, an incoherence model is employed to describe the relation between the instantaneous phase at the present station and the instantaneous phases at previous stations. This type of phase is defined as the instantaneous coherence phase. In addition, time lag is included in the instantaneous coherence phase to prescribe the wave passage effect. The proposed Hilberttransformbased method is efficient and avoids cumbersome parameter estimations as well as other drawbacks involved in some traditional synthesizing methods. Applications of this method demonstrate that the generated ground motions are statistically analogous with the reference record.
1. Introduction
In the past decades, the methods of generating spatially variable ground motions were studied thoroughly. A number of studies demonstrate that the effect of spatially variable ground motions on the responses of long structures is not negligible [1–7]. Nowadays, the synthesizing methods [8–11] based on stochastic process theory are popular and are implemented as a tool in most of the finiteelement packages [12]. Stochastic procedures usually include power spectral matrix/incoherence matrix decomposition [13, 14] or spectral factorization [15], which involves massive calculation, thus decreasing the synthesizing efficiency. In addition, the nonstationarity of ground motion is a significant factor in the process of generating ground motions, which includes intensity nonstationarity and frequency nonstationarity [16]. Intensity nonstationarity is usually achieved by multiplying the generated stationary process using an envelope function. However, this method may not be applicable for simulations based on a recorded accelerogram, because the applied envelop function may not completely present the intensity nonstationarity of the original ground motion. Moreover, studies by Ohsaki [17] showed that the seismic waveform was governed by the distribution of the phase difference spectrum; thus, phase difference spectrum could be used to achieve intensity nonstationarity. Zhu and Feng [18, 19] studied the distribution characteristic of the phase difference spectrum and proposed a method of generating a random phase. This method can achieve full nonstationarity; however, the effect of the randomness of the generated phase on the ground motions is not confirmed. Several methods have been proposed for frequency nonstationarity [20–26]. These methods involve complex processes and massive calculation. Conditional simulation is an alternative and several contributions have been made to the study of ground motion simulation [27, 28].
In the present paper, a simplified conditional simulation method of synthesizing spatially correlated earthquake ground motions is proposed based on Hilbert transform and earthquake record. In this method, one reference earthquake record is treated as the original ground motion, based on which a series of ground motions are generated. By performing Hilbert transform on the known earthquake record, the instantaneous amplitude and the instantaneous phase of the earthquake record can be obtained. The instantaneous amplitude is utilized as an envelope function to achieve the intensity nonstationarity of each simulated ground motion. To establish the coherency between generated ground motions, a coherence model is employed to describe the relation between the instantaneous phase at the present station and the instantaneous phases at previous stations. The instantaneous coherence phases at different stations are all statistically analogous with that of the known record. Thus every generated ground motion can present similar frequency nonstationarity. This type of phase is defined as an instantaneous coherence phase. In addition, time lag is included in the instantaneous coherence phase to prescribe the wave passage effect. The proposed Hilberttransformbased method can efficiently achieve intensity nonstationarity and frequency nonstationarity, and this approach avoids cumbersome parameter estimations, as well as other drawbacks involved in some traditional synthesizing methods. The application of this method demonstrates its validity and practical value.
2. Hilbert Transform
For a realvalued function , its Hilbert transform is defined aswhere denotes the Cauchy principal value of the integral. Thus, by using the Hilbert transform, an analytic signal which is a complexvalued function can be obtained as follows [29]:which can be further expressed as where .
Thus, the original function can be expressed as where the independent variable denotes time and time functions and are called the instantaneous amplitude and the instantaneous phase function or, by Bendat and Piersol [30], the envelop signal and the instantaneous phase signal of , respectively. In addition, represents the amplitude modulation, and represents the frequency modulation mechanisms contained in the original signal and is between –π and π. In other words, governs the intensity nonstationarity (or the temporal variation of amplitude), whereas dominates the frequency nonstationarity (or the temporal variation of frequency content) of the signal .
3. Proposed Method
In this study, the spatial variation of ground motion is prescribed in terms of wave scattering as well as wave passage effects. The soil condition and geology of the field of interest were assumed to be uniform. The wellknown northsouth component of the natural ground motion recorded at the El Centro station during the 1940 earthquake in Imperial Valley, California, which had a magnitude of 6.95 Mw, is chosen as the reference accelerogram, as shown in Figure 1(a). The recorded ground motion exhibits a sampling rate of 50 Hz and the anterior 40.94 s time history of the record is extracted to be studied with a total of 2048 values. First, the instantaneous amplitude and the instantaneous phase functions can be obtained by performing the Hilbert transform on the known earthquake record, as indicated in Figures 1(b) and 1(c), respectively. The instantaneous amplitude is treated by the method as an envelop function, with a value that will be preserved into the synthetic samples.
(a)
(b)
(c)
According to (5), the original earthquake ground motion can be expressed, and the acceleration process at station can be given aswhere is the time lag between the two stations given by [13]where denotes the separation distance between the two stations projected parallel to the dominant wave propagation direction and is assumed to be 200 m uniformly; represents the apparent seismic wave velocity in the medium, which is 500 m/s; and indicates the coherence function between two instantaneous phases. The instantaneous phase at the present station is assumed to be only affected by these phases at previous stations. Thus, the effect of phases at previous stations on the present phase can be expressed by . The relevant verification is shown in appendix.
Many coherency models were proposed, such as the models proposed by Loh and Lin [31], Feng and Hu [32], Loh and Yeh [33], Hao et al. [34], and Harichandran and Vanmarcke [35]. In this study, the Sobczyk model [36] is selected to describe the coherency loss between the ground motions at points and as follows:which reflects the level of coherency loss; is used in the present paper, which indicates highly correlated motions [37]; is the incident angle of the incoming wave to the site, and it is assumed to be 0°.
4. Verification of the Proposed Method
To verify the proposed method, a set of comparisons are conducted, including comparison of peak ground accelerations (PGA), peak ground velocity (PGV), and peak ground displacement (PGD); Fourier amplitude and response spectra amplitude; and coherence function of simulated motions to those of the known record.
4.1. Comparison of PGA, PGV and PGD
A series of accelerograms at three stations arranged at 200 m, 400 m, and 600 m from the original record are simulated using the proposed methodology. The acceleration, velocity, and displacement time histories at each point are presented in Figures 2–4, respectively. The waveform of time histories of acceleration, velocity, and displacement at each point are all consistent with the original record. The PGAs, PGVs, and PGDs of the simulated motions are compared with those of the original ground motion, which are presented in Table 1. The PGA, PGV, and PGD at every station are close to those of the known record. Besides, from those time histories figures, it can be seen that the generated ground motions have the identical intensity nonstationarity and frequency nonstationarity with the original record. These comparison results verify that the simulated ground motions are rational.

(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
4.2. Comparison of Fourier Amplitude and Response Spectra
Figure 5 provides the results of Fourier amplitude of simulated motions compared with those of the original record. Figure 6 presents the response spectra amplitude of simulated motions, as well as that of the original record. The spectral values are close to that of the known record. The comparison of the timefrequency spectra of the generated ground motions to that of the known record is shown in Figure 7. The timefrequency spectra were obtained by using the short time Fourier transform. The figure shows that both spectral amplitude and tendency are similar to that of the known record.
(a)
(b)
(c)
(d)
4.3. Comparison of Coherence
Figure 8 indicates the coherency values between the generated time histories and the known record as well as the corresponding Sobczyk model. A good match can also be observed except for and in the high frequency range. However, this result is expected because the cross correlation between the generated motions or their coherency values decrease rapidly with frequency as the separation increases [37]. A previous study [38] indicated that the coherency value of approximately 0.3–0.4 is the threshold of the cross correlation between two signals because numerical calculations of the coherency function between any two white noise series result in a value of roughly 0.3–0.4. Thus the calculated coherency function between two simulated time histories is smaller than the target value.
(a)
(b)
(c)
5. Conclusion
In this study, an approach was developed to simulate spatially correlated earthquake ground motions based on Hilbert transform and earthquake records where instantaneous coherence phase is defined. In this method, one recorded accelerogram is treated as the original ground motion, based on which a series of ground motions are generated. By performing Hilbert transform on the known earthquake record, the instantaneous amplitude and the instantaneous phase of the earthquake record can be obtained. The instantaneous amplitude is utilized as an envelop function to achieve the intensity nonstationarity of each simulated ground motion. To establish the coherency between generated ground motions, the Sobczyk model is employed in describing the relation between the instantaneous phase at the present station and the instantaneous phases at previous stations. The instantaneous coherence phases at different stations are all statistically analogous with that of the known record. Thus every generated ground motion can present similar frequency nonstationarity. In addition, time lag is included in the instantaneous coherence phase to prescribe the wave passage effect. The proposed Hilberttransformbased method is concise and avoids many drawbacks involved in some traditional synthesizing methods. The power spectra density of every simulation point is considered to be identical in the proposed method, so it is applicable to the relatively small field only. Because for largescale field, the power spectra density of each simulation point might be different. Thus, it has been considered that the local site effect should be incorporated in the synthesis method in the future.
Appendix
Verification for the Instantaneous Coherency Phase
As mentioned above, it is assumed that the soil condition and geology of the field of interest were uniform, so the power spectra density at each point is uniform and the coherency function between th and th point can be given as follows:where denotes a time lag. By using Hilbert transform, and can be expressed as
Substituting (A.3) into (A.2), can be expressed asand then
When , , (A.6) is valid. Herein, is set to be 0 and then (A.8) can be expressed as
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
Financial support for the project from the National Key Research and Development Program of China (no. 2016YFC0800206), the National Natural Science Foundation of China (no. 51378234 and no. 51678465), and the Fundamental Research Funds for the Central Universities (no. 2015MS060) is acknowledged.
References
 A. Hindy and M. Novak, “Pipeline response to random ground motion,” Journal of Engineering Mechanics ASCE, vol. 106, no. 2, pp. 339–360, 1980. View at: Google Scholar
 W. L. Liu, X. G. Wu, L. M. Zhang, J. Zheng, and J. Teng, “Global sensitivity analysis of tunnelinduced building movements by a precise metamodel,” Journal of Computing in Civil Engineering, vol. 31, no. 5, p. 04017037, 2017. View at: Publisher Site  Google Scholar
 A. S. Nazmy and A. M. Abdel‐Ghaffar, “Effects of ground motion spatial variability on the response of cable‐stayed bridges,” Earthquake Engineering & Structural Dynamics, vol. 21, no. 1, pp. 1–20, 1992. View at: Publisher Site  Google Scholar
 R. S. Harichandran, A. Hawwari, and B. N. Sweidan, “Response of longspan bridges to spatially varying ground motion a,” Journal of Structural Engineering, vol. 122, no. 5, pp. 476–484, 1996. View at: Publisher Site  Google Scholar
 M.T. Chen and R. S. Harichandran, “Response of an earth dam to spatially varying earthquake ground motion,” Journal of Engineering Mechanics, vol. 127, no. 9, pp. 932–939, 2001. View at: Publisher Site  Google Scholar
 A. A. Dumanogluid and K. Soyluk, “A stochastic analysis of long span structures subjected to spatially varying ground motions including the siteresponse effect,” Engineering Structures, vol. 25, no. 10, pp. 1301–1310, 2003. View at: Publisher Site  Google Scholar
 K. Soyluk, “Comparison of random vibration methods for multisupport seismic excitation analysis of longspan bridges,” Engineering Structures, vol. 26, no. 11, pp. 1573–1583, 2004. View at: Publisher Site  Google Scholar
 M. Shinozuka and G. Deodatis, “Stochastic process models for earthquake ground motion,” Probabilistic Engineering Mechanics, vol. 3, no. 3, pp. 114–123, 1988. View at: Publisher Site  Google Scholar
 A. Zerva, “Seismic ground motion simulations from a class of spatial variability models,” Earthquake Engineering & Structural Dynamics, vol. 21, no. 4, pp. 351–361, 1992. View at: Publisher Site  Google Scholar
 A. Zerva and M. Shinozuka, “Stochastic differential ground motion,” Structural Safety, vol. 10, no. 1–3, pp. 129–143, 1991. View at: Publisher Site  Google Scholar
 G. Deodatis, “Nonstationary stochastic vector processes: seismic ground motion applications,” Probabilistic Engineering Mechanics, vol. 11, no. 3, pp. 149–168, 1996. View at: Publisher Site  Google Scholar
 A. A. Shama, “Simplified procedure for simulating spatially correlated earthquake ground motions,” Engineering Structures, vol. 29, no. 2, pp. 248–258, 2007. View at: Publisher Site  Google Scholar
 Y.F. Gao, Y.X. Wu, and B. Li, “Simplified method for simulation of ergodic spatially correlated seismic ground motion,” Applied Mathematics and Mechanics. English Edition, vol. 32, no. 10, pp. 1297–1314, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Wu, Y. Gao, D. Li, C. Xu, and A. H. Mahfouz, “Study of the approximate approaches to the POD based spectral representation method,” Science China Technological Sciences, vol. 56, no. 4, pp. 970–979, 2013. View at: Publisher Site  Google Scholar
 M. Shrikhande and V. K. Gupta, “Synthesizing ensembles of spatially correlated accelerograms,” Journal of Engineering Mechanics, vol. 124, no. 11, pp. 1185–1192, 1998. View at: Publisher Site  Google Scholar
 S. Rezaeian and A. Der Kiureghian, “A stochastic ground motion model with separable temporal and spectral nonstationarities,” Earthquake Engineering and Structural Dynamics, vol. 37, no. 13, pp. 1565–1584, 2008. View at: Publisher Site  Google Scholar
 Y. Ohsaki, “On the significance of phase content in earthquake ground motions,” Earthquake Engineering & Structural Dynamics, vol. 7, no. 5, pp. 427–439, 1979. View at: Publisher Site  Google Scholar
 Zhu Y. and Feng Q. M., “Digital characteristics of phase difference spectrum distribution of earthquake accelerogram,” Earthquake Engineering & Engineering Vibration, vol. 12, no. 1, pp. 37–44, 1992. View at: Google Scholar
 Zhu Y. and Feng Q. M., “Distribution characteristic of phase difference spectrum and artificial accelerogram,” Earthquake Engineering & Engineering Vibration, vol. 13, no. 2, pp. 30–37, 1993. View at: Google Scholar
 M. B. Priestley, SpEctral Analysis and Time Series, Academic Press, 1981. View at: MathSciNet
 K. Sarkar and V. K. Gupta, “Waveletbased simulation of spatially correlated and spectrumcompatible accelerograms,” in Proceedings of the national symposium on structural dynamics, random vibrations earthquake engineering (NSSD '05), pp. 69–78, Department of Civil Engineering, Indian Institute of Science, Bangalore, 2005. View at: Google Scholar
 P. Cacciola and G. Deodatis, “A method for generating fully nonstationary and spectrumcompatible ground motion vector processes,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 3, pp. 351–360, 2011. View at: Publisher Site  Google Scholar
 T. Kubo and J. Penzien, “Simulation of three‐dimensional strong ground motions along principal axes, San Fernando earthquake,” Earthquake Engineering & Structural Dynamics, vol. 7, no. 3, pp. 279–294, 1979. View at: Publisher Site  Google Scholar
 Y. K. Lin and Y. Yong, “Evolutionary KanaiTajimi earthquake models,” Journal of Engineering Mechanics, vol. 113, no. 8, pp. 1119–1137, 1987. View at: Publisher Site  Google Scholar
 A. Der Kiureghian and J. Crempien, “An evolutionary model for earthquake ground motion,” Structural Safety, vol. 6, no. 2, pp. 235–246, 1989. View at: Publisher Site  Google Scholar
 J. P. Conte, K. S. Pister, and S. A. Mahin, “Nonstationary ARMA modeling of seismic motions,” Soil Dynamics and Earthquake Engineering, vol. 11, no. 7, pp. 411–426, 1992. View at: Publisher Site  Google Scholar
 L. Hu, Y. L. Xu, and Y. Zheng, “Conditional simulation of spatially variable seismic ground motions based on evolutionary spectra,” Earthquake Engineering and Structural Dynamics, vol. 41, no. 15, pp. 2125–2139, 2012. View at: Publisher Site  Google Scholar
 Y. Wu, Y. Gao, and D. Li, “Error assessment of multivariate random processes simulated by a conditionalsimulation method,” Journal of Engineering Mechanics, vol. 141, no. 5, Article ID 04014155, 2015. View at: Publisher Site  Google Scholar
 L. Cohen, TimeFrequency Analysis: Theory and Applications, PrenticeHall, Inc, Englewood Cliffs, New Jersey, 1995.
 J. S. Bendat and A. G. Piersol, Random Data: Analysis and Measurement Procedures, John & Wileys, New York, NY, USA, 1986. View at: MathSciNet
 C.H. Loh and S.G. Lin, “Directionality and simulation in spatial variation of seismic waves,” Engineering Structures, vol. 12, no. 2, pp. 134–143, 1990. View at: Publisher Site  Google Scholar
 Q. M. Feng and Y. X. Hu, “Spatial Correlation of Earthquake Motion and Its Effect on Structural Response, Bilateral Workshop on Earthquake Engineering, Beijing,” 1982. View at: Google Scholar
 C. Loh and Y. Yeh, “Spatial variation and stochastic modelling of seismic differential ground movement,” Earthquake Engineering & Structural Dynamics, vol. 16, no. 4, pp. 583–596, 1988. View at: Publisher Site  Google Scholar
 H. Hao, C. S. Oliveira, and J. Penzien, “Multiplestation ground motion processing and simulation based on smart1 array data,” Nuclear Engineering and Design, vol. 111, no. 3, pp. 293–310, 1989. View at: Publisher Site  Google Scholar
 R. S. Harichandran and E. H. Vanmarcke, “Stochastic variation of earthquake ground motion in space and time,” Journal of Engineering Mechanics, vol. 112, no. 2, pp. 154–174, 1986. View at: Publisher Site  Google Scholar
 K. Sobczyk, StochasTic Wave Propagation, Kluwer Academic Publishers, Netherlands, 1991.
 K. Bi and H. Hao, “Modelling and simulation of spatially varying earthquake ground motions at sites with varying conditions,” Probabilistic Engineering Mechanics, vol. 29, pp. 92–104, 2012. View at: Publisher Site  Google Scholar
 H. Hao, EffEcts of Spatial Variation of Ground Motions on Large Multiple Supported Structures UCB/EERC8906, University of California, Berkeley, CA, USA, 1989. View at: MathSciNet
Copyright
Copyright © 2017 Erlei Yao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.