Effect of a large-scale three-dimensional sedimentary basin on Rayleigh wave propagation by using spectral element method combined with frequency-wavenumber method

Abstract

Site effects study has always been a key research topic in earthquake engineering. This study proposes a hybrid method to analyze large-scale three-dimensional sedimentary basin under Rayleigh (R) wave incidence. The proposed hybrid method includes two steps: 1) calculate the free field responses of layered sites subjected to R-wave using the frequency-wavenumber method; 2) Simulate the local site region using spectral element method with the equivalent forces input computed from the free field responses. A comprehensive verification study is conducted demonstrating the accuracy of this method. To investigate the effect of sedimentary basin on R-wave propagation, a parametric study is performed on the medium impedance contrast ratio of sedimentary basins and the incident seismic wave predominant frequency, revealing the scattering patterns of sedimentary basins under R-wave incidence. Finally, a practical case of the Wudu Basin in the Tibetan Plateau region of China is simulated. Results indicate significant amplification of R-wave by sedimentary basin, and the proposed hybrid method could serve as a reliable and efficient approach for large-scale R-wave propagation simulation.

Graphical abstract

Keywords

Rayleigh wave / hybrid method / basin effect / wave propagation / spectral element method

References

[1] Sánchez-Sesma F J, Francisco J, Velázquez S A. On the seismic response of a dipping layer. Wave Motion, 1987, 9(5): 387–391

[2] Aki K. Local site effects on weak and strong ground motion. Tectonophysics, 1993, 218(1–3): 93–111

[3] Kawase H. The cause of the damage belt in Kobe: “The basin-edge effect,” constructive interference of the direct S-wave with the basin-induced diffracted/Rayleigh waves. Seismological Research Letters, 1996, 67(5): 25–34

[4] Kakoty P, Dyaga S M, Molina Hutt C. Impacts of simulated M9 Cascadia subduction zone earthquakes considering amplifications due to the Georgia sedimentary basin on reinforced concrete shear wall buildings. Earthquake Engineering & Structural Dynamics, 2021, 50(1): 237–256

[5] Marafi N A, Eberhard M O, Berman J W, Wirth E A, Frankel A D. Effects of deep basins on structural collapse during large subduction earthquakes. Earthquake Spectra, 2017, 33(3): 963–997

[6] Anderson J G, Bodin P, Brune J N, Prince J, Singh S K, Quaas R, Onate M. Strong ground motion from the Michoacan, Mexico, earthquake. Science, 1986, 233: 1043–1049

[7] Pitarka A, Irikura K, Iwata T, Sekiguchi H. Three-dimensional simulation of the near-fault ground motion for the 1995 Hyogo-ken Nanbu (Kobe), Japan, earthquake. Bulletin of the Seismological Society of America, 1998, 88(2): 428–440

[8] Brissaud Q, Bowden D C, Tsai V C. Extension of the basin Rayleigh-wave amplification theory to include basin-edge effects. Bulletin of the Seismological Society of America, 2020, 110(3): 1305–132

[9] Trifunac M D. Surface motion of a semi-cylindrical alluvial valley for incident plane SH waves. Bulletin of the Seismological Society of America, 1971, 61(6): 1755–1770

[10] Todorovska M I, Lee V W. A note on scattering of Rayleigh waves by shallow circular canyons: Analytical approach. Bulletin of the Indian Society of Earthquake Technology, 1991, 28(2): 1–16

[11] Poursartip B, Fathi A, Tassoulas J L. Large-scale simulation of seismic wave motion: A review. Soil Dynamics and Earthquake Engineering, 2020, 129: 105909

[12] Kawase H, Aki K. A study on the response of a soft basin for incident S, P, and Rayleigh waves with special reference to the long duration observed in Mexico City. Bulletin of the Seismological Society of America, 1989, 79(5): 1361–1382

[13] Semblat J F, Duval A M, Dangla P. Seismic site effects in a deep alluvial basin: Numerical analysis by the boundary element method. Computers and Geotechnics, 2002, 29(7): 573–585

[14] Sánchez-Sesma F J, Ramos-Martínez J, Campillo M. An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves. Earthquake Engineering & Structural Dynamics, 1993, 22(4): 279–295

[15] Sánchez-Sesma F J, Luzon F. Seismic response of three-dimensional alluvial valleys for incident P, S, and Rayleigh waves. Bulletin of the Seismological Society of America, 1995, 85(1): 269–284

[16] Frankel A, Vidale J. A three-dimensional simulation of seismic waves in the Santa Clara Valley, California, from a Loma Prieta aftershock. Bulletin of the Seismological Society of America, 1992, 82(5): 2045–2074

[17] Opršal I, Zahradnik J. Three-dimensional finite difference method and hybrid modeling of earthquake ground motion. Journal of Geophysical Research. Solid Earth, 2002, 107(B8): ESE–2

[18] Wang L, Xu Y, Xia J, Luo Y. Effect of near-surface topography on high-frequency Rayleigh-wave propagation. Journal of Applied Geophysics, 2015, 116: 93–103

[19] Zhang X, Peng X, Li X, Zhou Z, Mebarki A, Dou Z, Nie W. Seismic effects of a small sedimentary basin in the eastern Tibetan plateau based on numerical simulation and ground motion records from aftershocks of the 2008 Mw7.9 Wenchuan, China earthquake. Journal of Asian Earth Sciences, 2020, 192: 104257

[20] Qiang S, Wang H, Wen R, Liu Q, Zhou Y. Investigating the effects of structural parameters on seismic aggravation of two-dimensional sedimentary valleys. Soil Dynamics and Earthquake Engineering, 2023, 171: 107964

[21] Komatitsch D, Vilotte J P. The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures. Bulletin of the Seismological Society of America, 1998, 88(2): 368–392

[22] Komatitsch D, Liu Q, Tromp J, Suss P, Stidham C, Shaw J H. Simulations of ground motion in the Los Angeles basin based upon the spectral-element method. Bulletin of the Seismological Society of America, 2004, 94(1): 187–206

[23] Shen W, Yang D, Xu X, Yang S, Liu S. 3D simulation of ground motion for the 2015 Mw7.8 Gorkha earthquake, Nepal, based on the spectral element method. Natural Hazards, 2022, 112(3): 2853–2871

[24] Touhami S, Gatti F, Lopez-Caballero F, Cottereau R, de Abreu Corrêa L, Aubry L, Clouteau D. SEM3D: A 3D high-fidelity numerical earthquake simulator for broadband (0–10 Hz) seismic response prediction at a regional scale. Geosciences, 2022, 12(3): 112

[25] Ba Z, Zhao J, Zhu Z, Zhou G. 3D physics-based ground motion simulation and topography effects of the 05 September 2022 Mw6.6 Luding earthquake, China. Soil Dynamics and Earthquake Engineering, 2023, 172: 108048

[26] Noureddine M A, de Martin F, Elmeouche R, Ababsa F, Beaufils M, Sammuneh M A. Source to site seismic ground motion prediction on an earth dam. Japanese Geotechnical Society Special Publication, 2024, 10(31): 1188–1193

[27] Castro-Cruz D, Gatti F, Lopez-Caballero F, Hollender F, El Haber E, Causse M. Blind broad-band (0–10 Hz) numerical prediction of the 3-D near field seismic response of an Mw6.0 extended fault scenario: application to the nuclear site of Cadarache. Geophysical Journal International, 2023, 232(1): 581–600

[28] De Martin F, Chaljub E, Thierry P, Sochala P, Dupros F, Maufroy E, Hollender F. Influential parameters on 3D synthetic ground motions in a sedimentary basin derived from global sensitivity analysis. Geophysical Journal International, 2021, 227(3): 1795–1817

[29] Tong P, Chen C W, Komatitsch D, Basini P, Liu Q. High-resolution seismic array imaging based on an SEM-FK hybrid method. Geophysical Journal International, 2014, 197(1): 369–395

[30] Tong P, Komatitsch D, Tseng T L, Hung S H, Chen C W, Basini P, Liu Q A. 3-D spectral-element and frequency-wave number hybrid method for high-resolution seismic array imaging. Geophysical Journal International, 2014, 41(20): 7025–7034

[31] Liang J, Wu M, Ba Z, Liu Y. A hybrid method for modeling broadband seismic wave propagation in 3D localized regions to incident P, SV, and SH waves. International Journal of Applied Mechanics, 2021, 13(10): 2150119

[32] MoncadaV SLopez-caballero F. Coupled approach for the assessment of basin effects on the seismic demand of nonlinear structures. In: Proceedings of 8th International Conference on Earthquake Geotechnical Engineering. Osaka: Japanese Geotechnical Society Special Publication, 2024

[33] SotoVLopez-Caballero F. Integrated assessment of basin effects on seismic damage: A coupled 3D SEM-FEM approach with domain reduction method. 2024 (available at the website of Research Square)

[34] Korres M, Lopez-Caballero F, Alves Fernandes V, Gatti F, Zentner I, Voldoire F, Castro-Cruz D. Enhanced seismic response prediction of critical structures via 3D regional scale physics-based earthquake simulation. Journal of Earthquake Engineering, 2023, 27(3): 546–574

[35] Miller G F, Pursey H, Bullard E C. On the partition of energy between elastic waves in a semi-infinite solid. Mathematical and Physical Sciences, 1955, 233(1192): 55–69

[36] Wong H L. Effect of surface topography on the diffraction of P, SV, and Rayleigh waves. Bulletin of the Seismological Society of America, 1982, 72(4): 1167–1183

[37] Yang Y, Yu H, Yuan Y, Lu D, Huang Q. A numerical framework for underground structures in layered ground under inclined P-SV waves using stiffness matrix and domain reduction methods. Frontiers of Structural and Civil Engineering, 2023, 17(1): 10–24

[38] YangYZhou ZWangXZhangXWangZ. 2.5-dimension soil seismic response to oblique incident waves based on exact free-field solution. Frontiers of Structural and Civil Engineering, 2024, 18(2): 216-235

[39] Zhang J, Lan J, Li H, Ba Z, Zhuang H, Xu Z. Seismic wave-field modeling of 3D irregularities in layered half-space under P-, S-, and Rayleigh waves with arbitrary incident azimuth. Soil Dynamics and Earthquake Engineering, 2024, 182: 108669

[40] Nguyen K T, Kusanovic D S, Asimaki D. Three-dimensional nonlinear soil–structure interaction for Rayleigh wave incidence in layered soils. Earthquake Engineering & Structural Dynamics, 2022, 51(11): 2752–2770

[41] Bouchon M, Aki K. Discrete wave-number representation of seismic-source wave fields. Bulletin of the Seismological Society of America, 1977, 67(2): 259–277

[42] HaskellN A. Vincit Veritas: A Portrait of the Life and Work of Norman Abraham Haskell. Washington, D.C.: American Geophysical Union, 1990

[43] Zhao X, Zhu J B, Zhao J, Cai J. Study of wave attenuation across parallel fractures using propagator matrix method. International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(10): 1264–1279

[44] Singla V K, Gupta V K. Surface rotations due to kinematic shear dislocation point source in a multilayered elastic medium. Bulletin of the Seismological Society of America, 2019, 109(1): 433–447

[45] Kausel E. Generalized stiffness matrix method for layered soils. Soil Dynamics and Earthquake Engineering, 2018, 115: 663–672

[46] Kausel E, Roësset J M. Stiffness matrices for layered soils. Bulletin of the Seismological Society of America, 1981, 71(6): 1743–1761

[47] Cao Z, Tao X, Tao Z, Tang A. Kinematic source modeling for the synthesis of broadband ground motion using the F-K approach. Bulletin of the Seismological Society of America, 2019, 109(5): 1738–1757

[48] WolfJ P. Dynamic Soil–Structure-Interaction. Englewood Cliffs, NJ: Prentice-Hall, Inc., 1985

[49] Ba Z, Liang J, Lee V W. Wave propagation of buried spherical SH-, P1-, P2-and SV-waves in a layered poroelastic half-space. Soil Dynamics and Earthquake Engineering, 2016, 88: 237–255

[50] Ba Z, Niu J, Liu Y, Liang J. Dynamic response of a multi-scale layered saturated porous half-space due to seismic dislocation source by using a revised dynamic stiffness matrix method. Applied Mathematical Modelling, 2023, 120: 217–245

[51] Ba Z, Liu Y, Liang J, Sang Q, Wu M, Zhang Y. The dynamic stiffness matrix method for seismograms synthesis for layered transversely isotropic half-space. Applied Mathematical Modelling, 2022, 104: 205–227

[52] Ba Z, Sang Q, Wu M, Liang J. The revised direct stiffness matrix method for seismogram synthesis due to dislocations: from crustal to geotechnical scale. Geophysical Journal International, 2021, 227(1): 717–734

[53] Komatitsch D, Tromp J. Spectral-element simulations of global seismic wave propagation—I. Validation. Geophysical Journal International, 2002, 149(2): 390–412

[54] Komatitsch D, Tromp J. Spectral-element simulations of global seismic wave propagation—II. Three-dimensional models, oceans, rotation and self-gravitation. Geophysical Journal International, 2002, 150(1): 303–318

[55] Clayton R, Engquist B. Absorbing boundary conditions for acoustic and elastic wave equations. Bulletin of the Seismological Society of America, 1977, 67(6): 1529–1540

[56] Mossessian T K, Dravinski M. Amplification of elastic waves by a three dimensional valley. Part 1: Steady state response. Earthquake Engineering & Structural Dynamics, 1990, 19(5): 667–680

RIGHTS & PERMISSIONS

Higher Education Press