3D NUMERICAL SIMULATIONS OF EARTHQUAKE 3D NUMERICAL SIMULATIONS OF EARTHQUAKE GROUND MOTION IN SEDIMENTARY BASINS: GROUND MOTION IN SEDIMENTARY BASINS: THE CASES OF GUBBIO AND L’AQUILA, THE CASES OF GUBBIO AND L’AQUILA, CENTRAL ITALY CENTRAL ITALY Roberto Paolucci and Chiara Smerzini Department of Structural Engineering, Politecnico di Milano
29
Embed
3D NUMERICAL SIMULATIONS OF EARTHQUAKE GROUND MOTION IN SEDIMENTARY BASINS: THE CASES OF GUBBIO AND L’AQUILA, CENTRAL ITALY Roberto Paolucci and Chiara.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
3D NUMERICAL SIMULATIONS OF EARTHQUAKE GROUND 3D NUMERICAL SIMULATIONS OF EARTHQUAKE GROUND MOTION IN SEDIMENTARY BASINS: THE CASES OF MOTION IN SEDIMENTARY BASINS: THE CASES OF
GUBBIO AND L’AQUILA, CENTRAL ITALYGUBBIO AND L’AQUILA, CENTRAL ITALY
Roberto Paolucci and Chiara Smerzini
Department of Structural Engineering, Politecnico di Milano
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
2ContentsContents
Motivation for 3D numerical simulations of earthquake ground motion
The spectral element code GeoELSE
Case studies
Seismic response of the Gubbio basin during the 1997 Umbria-Marche earthquake
Modeling of the MW 6.3 2009 L’Aquila earthquake
Conclusions
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
3
To simulate “synthetic earthquakes” as realistic as possible in terms of:
the complexity of the seismic source
the complexity of the geological and morphological environment
PGV maps in the Grenoble Valley due to a Mw6 earthquake along the Belledonne fault. From left to right: neutral, forward, backward directivity conditions with respect to the urban area of Grenoble.
After Stupazzini et al., 2009.
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
5
integration to PSHA, especially for long return periods
- CyberShake (Graves et al., 2010)
- S2 Project DPC-INGV 2007-2009 (Faccioli et al, 2010)
seismic risk assessment of urban areas under scenario
earthquakes
ShakeOut Scenario: Southern California (Tech. report, 2008)
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
7ContentsContents
Motivation for 3D numerical simulations of earthquake ground motion
The spectral element code GeoELSE
Case studies
Seismic response of the Gubbio basin during the 1997 Umbria-Marche earthquake
Modeling of the MW 6.3 2009 L’Aquila earthquake
Conclusions
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
8
Developers
Department of Structural Engineering, Politecnico di Milano E. Faccioli, R. Paolucci, L. Scandella, C. Smerzini, M.Stupazzini, M. Vanini
CRS4 (Center of Advanced Studies, Research and Development in Sardinia) F. Maggio, L. Massidda
Department of Modeling and Scientific Computing (MOX), Politecnico di Milano P. Antonietti, I. Mazzieri, A. Quarteroni, F. Rapetti
Web site: http://geoelse.stru.polimi.it
The Spectral Element code The Spectral Element code GeoELSEGeoELSE
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
9
Main purpose of GeoELSE
Studying 2D/3D linear and non-linear visco-elastic seismic wave propagation in heterogeneous media, including within the same numerical model:- seismic source (extended fault / plane wave with arbitrary incidence angle)- propagation path
- complex geological structures / SSI effects
The Spectral Element code The Spectral Element code GeoELSEGeoELSE
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
10The Spectral Element code The Spectral Element code GeoELSEGeoELSE
Dynamic response of structuresSeismic wave propagation in complex
geological configurations
L’Aquila basin
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
11
Kosloff D, Baysal E. Forward modelling by the Fourier methodGeophysics 1982 47: 1402-1412.
Kosloff D, Kessler D, Filho AQ, Tessmer E, Behle A, Strahilevitz R. Solutions of the equations of dynamics elasticity by a Chebyshev spectral method Geophysics 1990; 55: 748-754.
Faccioli E, Maggio F, Paolucci R, Quarteroni A. 2D and 3D elastic wave propagation by a pseudo-spectral domain decomposition methodJournal of Seismology 1997; 1 237-251.
Komatitsch D, Vilotte J-P. The spectral element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures. Bull. Seism. Soc. Am. 1998; 88: 368-392.
Some “historical” references on spectral approaches for the numerical integration of the wave equation
The Spectral Element code The Spectral Element code GeoELSEGeoELSE
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
slip distribution according to Walters et al. (2009)
Effect of stochastic source parameters
3D numerical simulations of the M3D numerical simulations of the MWW6.3 L6.3 L’’Aquila earthquake Aquila earthquake
AQK
AQV
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
Heterogeneous kinematic parameters, defined by spatially correlated stochastic fields for rise time, rup. velocity and rake angle, with correlation length 4 km
Effect of stochastic source parameters
25
AQK
AQV
3D numerical simulations of the M3D numerical simulations of the MWW6.3 L6.3 L’’Aquila earthquake Aquila earthquake
AQK
AQV
slip rise time rup.vel rake
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
263D numerical simulations of the M3D numerical simulations of the MWW6.3 L6.3 L’’Aquila earthquake Aquila earthquake
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
27
Observed Simulated
Comparison with observed MCS intensity
Model CM1
3D numerical simulations of the M3D numerical simulations of the MWW6.3 L6.3 L’’Aquila earthquake Aquila earthquake
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO
28
3D numerical simulations of earthquake ground motion in near-fault conditions, accounting for complex geological and morphological conditions, may provide realistic seismic scenarios, up to frequencies of 2 – 3 Hz.
The frequency limit is mainly related to insufficient details in the source kinematic models, as well as on the local geology description. A moderate random variability of the kinematic source parameters may significantly improve the high-frequency energy radiation, improving as well the agreement with observed records during L’Aquila earthquake.
The typical features of long period ground motion amplification and propagation of surface waves within sedimentary basins in Central Italy, such as in Gubbio, can be captured well by 3D numerical simulations.
Generation of realistic earthquake ground motion scenarios for future damaging earthquakes within complex tectonic and geological environments is becoming more and more feasible, also for engineering applications.
Conclusions Conclusions
Roberto Paolucci: 3D numerical simulations of earthquake ground motion POLITECNICO DI MILANO