1. Introduction
In recent years, particular attention has been paid to the local seismic amplification phenomena, intensifying the geological studies of the territories and the soils’ seismic characteristics to understand the causes that can trigger seismic waves’ surface focalizations [
1].
The numerical estimates of the surface seismic amplification are influenced by various factors, including the selected analysis method, the dynamic properties of the soil, the input motion, and the shear wave velocity profiles
. In this paper, in order to introduce a stochastic approach, the
parameter was selected as a random parameter [
2]. The
can be evaluated with various geophysical methods (active and/or passive) or through empirical correlations [
3]. Currently, engineering design codes have few guidelines regarding the inclusion of
numerical uncertainty values for performing site response analyses. The uncertainties in the definition of the shear wave velocity can generate a numerical seismic response that is not always representative of the real surface amplifications. A common practice is the single deterministic analysis, using an average
value for each, eventually multiple, layers. However, this kind of approach does not commonly guarantee a correct and, as far as possible, complete estimate of the implications due to the variability of the selected geotechnical parameters. Accordingly, in this paper, the probabilistic method, for which each physical parameter can be assumed as a function of a pseudo-random variable
was assumed. All compilers (Fortran, C++, Python, etc.) use default algorithms which, based on an internal clock and on an extremely long series of numbers, even with 10
36 elements, are able to select numerical values that are not completely random, and for this reason, they are called “pseudo-random” numbers.
Currently, a large international bibliography describes in detail the characteristics of the stochastic method applied to seismic amplification. An attempt to couple Probabilistic Seismic Risk Analysis and seismic micro-zonation, two of the most important components of seismic risk mitigation strategies, has been discussed in [
4]. In the article [
2], for the stochastic study of seismic risk, the variability of the
was considered, and in [
5], the effects resulting from the choice of three different ways of considering the
were discussed. In [
6], the Monte Carlo method was applied. The parameters that substantially influence the surface seismic amplification depend on the deformation in a non-linear way [
7] (among many others). Accordingly, a technique to avoid the complexity of full non-linear numerical approaches has been proposed in [
8] called “The Equivalent-Linear Method”. By this approach, a pseudo iterative linear analysis is performed, with some initial values assumed for damping ratio and shear modulus. The maximum cyclic shear strain is recorded for each element and used to determine new values for damping and modulus by reference to laboratory-derived curves that relate damping ratio and secant modulus to the amplitude of cycling shear strain.
Some empirical scaling factor is usually used when relating laboratory strains to model strains. The new values of damping ratio and shear modulus are then used in a new updated numerical run. The whole process is repeated several times until there are no further changes in properties. At this point, it is said that “strain-compatible” values of damping and modulus have been found, and the simulation using these values is representative of the response of the real site.
It is important to note that, to introduce stochastic variability in numerical modeling, there are two main commonly used techniques, the Spectral Approach [
9] and the Monte Carlo method, which is more time-consuming, but more realistic.
In the practice of engineering modelling, but also in research activities, a commonly internationally used software is the Quad4M [
10], based on the 2D finite element method. The non-linearity of the soil is considered by “The Equivalent-Linear Method”. However, the parameters are entered in a deterministic way. To overcome this aspect, VisualQ4M, a new software (
www.visualq4m.com, “URL (accessed on 16 January 2023)” [
11], was developed, including a graphical user interface (GUI) for the pre- and post-processor with functionality and visual command able to generate, with excellent accuracy, geometric models suitable for the Quad4M program. The new software, selected to perform simulations discussed in this paper, manages all the seismic analyses in a single integrated interface to perform a complete 2D Quad4M simulation, producing various color-map results and spectral responses, as well. In addition, the selected software allows for the automatic performance of the Monte Carlo approach, even with a distribution of random values defined by the user, for example, variables with depth.
Accordingly, the Monte Carlo approach was applied to analyze different possible realizations of the system under consideration, resulting from different, but possible, statistical distributions of the parameter, selected as a random variable.
The spatial domain was discretized by geometric elements of varying size depending on location, but whose minimum dimension was no less than one meter. Therefore, in order to simplify the approach, we assumed that the autocorrelation and cross-correlation scales were, for the system under consideration, below 1.0 m, thus avoiding the elaboration of a spatially random distribution of geotechnical parameters using approaches such as those developed by Metropolis–Hastings [
12] and others [
13,
14]. This option will be implemented in future activities.
For the test case, a site located in the city of L’Aquila (Italy) was selected [
15]. The actual topographic shape is the result of the anthropic excavation realized at the beginning of the 2000s. The natural geometry was modified to realize a private construction in the valley area. About 8000 m
3 of soil was removed for the new project, and a prefabricated retaining wall with a height of about 7.50 m was realized. After the main earthquake of April 2009, the upstream structure has been completely damaged. In addition to the classic X-shaped pattern cracks (typical of the seismic action), other types of fissures also appeared, resulting from a sliding movement and minimal vertical deformation of the building. The reduced restriction downstream of the structure generated a relaxation of the horizontal stress component and an increase in the horizontal surface accelerations by topographical amplification. In this case, the simplified approach (1D model) [
16] was not correct to simulate multiple reflections of seismic waves on the surface. Thus, 2D simulations were correctly selected.
Moreover, only the fundamental frequencies of the ground have been considered and not those of the structures. Moreover, it should be noted that the Monte Carlo approach, based on the generation of random numbers for the determination of shear wave velocity, introduces frequencies that are not included in the natural frequencies of the selected site [
17]. Accordingly, not all seismic frequencies coming from the Fourier analyses of the seismic spectrum, selected as featuring the large scale territory within which is located the particular site under study, could impact the small scale of the site under consideration.
Accordingly, many (pseudo) random, parametric numerical analyses were performed to simulate the seismic wave propagations. In addition, an approach of this type was also useful to understand the impact of a modification of the topography on local seismic analysis after an excavation, including uncertainties considerations as well. Another useful result, based on the statistical standard deviation variability, was the experimental numerical test of how many Monte Carlo simulations were necessary to be performed.
4. Discussion and Conclusions
In this paper, the dynamic characterization of the seismic response analysis was based on the shear wave profile. The uncertainties in the definition of the shear wave velocity can generate seismic responses with results that may not be representative of the surface amplifications. Thus, in order to estimate a reasonable seismic response, despite the uncertainties, the Monte Carlo approach was selected. The outcomes of this kind of method were pseudo-random distributions of numerical values according to variability intervals. Sequences of pseudo-random distributions with two geometric configurations were calculated. For each realization of the system under consideration, seven natural waveforms were considered, and for each of them, 25 distinct simulations were performed. The spectrum response was calculated as the average of all 175 resulting spectra. The analysis was carried out according to the equivalent linear model. The calculations were performed with two distinct approaches: deterministic, for which a unique initial
wave velocity was considered for each layer, and probabilistic, for which the
wave velocity was calculated as a pseudo-random variable for each element. In the last case, the pseudo-random distribution model was re-calculated at each simulation. Subsequently, only for the probabilistic approach with the actual geometrical profile, the results were compared with the seismic site surveys to select only the realizations compatible with the real conditions. This kind of selection was performed through a band-pass filter (BPF) focused on the fundamental frequency of the site. Comparing the results of the two geometries with the deterministic approach, it is evident that the impact of the topographic modification was important. In fact, from
Table 6, it can be observed that the inclusion of the excavation in the numerical simulations caused an increase in the maximum spectral acceleration, raising it from 0.978 g to 1.051, 7% higher, as also displayed in
Figure 14, while the period, 0.36 s, was not affected. This aspect is more marked in the higher frequencies.
The topographic impact is more evident in the models analyzed with a probabilistic approach. The outcomes of the model, resulting from the inclusion of pseudo-random distribution parameters (shear wave) and the selection of the seven reference accelerograms), show, as can be observed from
Table 9 and also from
Figure 15, an increase in the maximum spectral acceleration from 0.974 g before excavation to 1.096 after excavation, or 12.5% higher. Unlike the deterministic case, with the probabilistic approach, the period related to the max value of the acceleration spectrum also changes, decreasing from 0.28 s to 0.24s, 14.4% lower (
Table 9). Accordingly, the numerical simulations carried out with the probabilistic approach show a period decrease compared to the values obtained with the deterministic approach of 22.2% and 33.3% pre-excavation and post-excavation, respectively. In order to highlight further important differences emerging from the choice of the methodological approach,
Figure 16 displays the comparison between the numerical outcomes resulting from the simulations based on probabilistic and deterministic approaches, both related to the geometry after excavation and the application of the band-pass filter (BPF) centered on the fundamental frequency of the site.
The models based on pseudo-random distributions generated responses with greater values of the amplification factors.
Of course, these differences are also a function of the amount of available data, the laboratory tests, and mainly the quality of the site investigations. Furthermore, in this work, only the (and therefore only the variability of the shear modulus) was considered as a reference parameter, but the surface seismic response could also be a function of the variability influence of other reference parameters or a combination of the same. Using numerical analyses with a pseudo-random approach certainly allows the reduction in the impact due to uncertainties on the determination of seismic parameters and, accordingly, on the estimation of as realistic as possible surface seismic amplification. Moreover, it is believed that not all pseudo-random combinations are representative of the site (especially if a variability on several parameters is used). The comparison with the field survey data (direct and/or indirect tests) is, however, fundamental for the validation of the numerical results.