The seismic safety assessment of structures and components of nuclear power plants is based on a well-established methodology consisting in the (i) establishment of the site-specific seismic hazard, the (ii) computation of the dynamic response of structures interacting with the surrounding soil to compute seismic demand, the iii) evaluation of the seismic fragility of structures and components and iv) the convolution of hazard and fragility for safety assessment. The development of best-estimate approaches together with comprehensive consideration of uncertainties as well as the development of integrated approaches from hazard to structural response allow for more accurate risk estimates (e.g. SINAPS@1 and METIS2 research projects). Recently, an integrated methodology was proposed by Korres et al. [1] accounting for the seismic wave propagation, from the source, through the complex geological medium in the regional and local scale of the site under consideration, and up to the structure’s foundation and its transmission to the sensitive equipment of interest (e.g., electro-mechanical devices, piping systems, etc.). However, the link between physics-based simulation and regulatory approaches where seismic load is defined by response spectra still needs to be fully established. On the other hand, source-to-structure 3D simulation remain computationally expensive and requires the development of surrogate models in particular in the light of uncertainty propagation and sensitivity analysis. In this context, the main objective of this work is to develop a global methodology for seismic assessment of structures integrating: (i) a source-to-structure wave propagation modeling conditioned on the site-specific seismic hazard defined on the outcropping rock, and (ii) uncertainties at the different scales, while propagating them through the different steps of the procedure.
The proposed strategy for this PhD work is based on ideas synthesized in the following:
Source-to-structure 3D physics-based simulations (PBS):
In order to appropriately model the seismic wave propagation at the different scales of analysis (regional scale ~km, site/structure scale ~m), a spectral element (SEM) - finite element method (FEM) coupling based on the domain reduction method [2, 3] already implemented in [1] is used as a basis for the numerical simulations. In this context, important components of the source-to- structure wave propagation [4, 5, 6] can be explicitly taken into account. More precisely, such numerical approach requires : i) the definition of a realistic 3D geology (regional/local scale) and the mechanical properties of the structure [7, 8], ii) the characterization and modeling of the excitation source on the regional scale [9, 10], iii) the definition of the ground motion on the scale of the site (hazard analysis - [11, 12], iv) the constitutive model for the description of soil’s or structure’s behavior (structural analysis - [4]), and v) the way that the dynamic excitation is imposed as an input to the soil-structure interaction model [5, 13];
Uncertainty quantification and propagation:
Uncertainties in the source-to-structure wave propagation modeling need also to be considered [14]. In this context, several studies have been proposed in the literature in order to account for these uncertainties : i) random fields for material properties [15, 16], ii) source modelling variability [17, 18], iii) sensitivity analysis for constitutive law parameters [19], iv) influence of the input motion excitation [3, 1]. Since hundreds to thousands of simulations are necessary for uncertainty propagation analyses in Monte-Carlo or more advances sampling strategies, surrogate models (Gaussian processes, Kriging, Polynomial chaos, etc.) can be constructed to obtain robust estimates of intensity measures (IM) [15] based on a reduced set of numerically costly high-fidelity simulations.
Site-specific hazard definition on the outcropping rock:
Finally, to account for the site-specific hazard definition, ground motion generated from the 3D PBS simulation need to comply with the seismic hazard defined at the outcropping rock. Hazard is generally determined using a probabilistic seismic hazard analysis (PSHA), where the seismic ground motion is defined by seismic response spectra The latter are an input for our analysis. An ultimate objective of this work is to create a link between the ground motion defined at the outcropping rock, and the response spectrum at the structural level, used for the verification of structural and non-structural components.