Sensitivity Analysis for Finite‐Difference Seismic Basin Modeling: A Case Study for Kinburn Basin, Ottawa, Canada

2019 ◽  
Vol 109 (6) ◽  
pp. 2305-2324 ◽  
Author(s):  
Amin Esmaeilzadeh ◽  
Dariush Motazedian

Abstract We used a finite‐difference modeling method to investigate the sensitivity of the ground‐motion simulation results to the main input parameters, including the source model, regional path properties, and local site conditions. We used a spectral frequency range of 0.1–1 Hz for the Kinburn bedrock topographic basin, Ottawa, Canada, for the Ladysmith earthquake (Mw 4.7). Some findings are known facts; however, the unique geophysical conditions in the Ottawa area, such as the high contrast between the shear‐wave velocities of the bedrock and the shear‐wave velocity of the soil, were the reason for a comprehensive sensitivity analysis. Using a Gaussian source function with a short half‐duration increased the peak ground velocities (PGVs) and the amplitude of the velocity Fourier spectrum. Relaxation times and relaxation coefficients for the viscoelastic simulation significantly increased the amplitude of later arrivals at the soil site, which, consequently, led to an increase in PGV, the amplitude of the pseudospectral acceleration (PSA) ratio, and the velocity Fourier spectrum for a small earthquake. Employing a small soil Q model damped a significant amount of energy of the waves in the basin; thus, PGV, the PSA of soil to rock ratios, and the velocity Fourier spectrum were dependent on the soil Q model. Also, using a high‐velocity contrast between soil and rock increased PGVs and the amplitude of the PSA of the soil to rock ratios, whereas the frequency content of the waves shifted toward lower frequencies. Using a finite‐fault source model for a scenario large earthquake (Mw 7) significantly reduced the PGV values relative to a point‐source model. Using nonlinear‐viscoelastic simulation for a large earthquake (Mw 7) reduced the amplitude of the later arrivals and the amplitude of the PSA of the soil to rock ratios, and shifted the frequency content of waves toward lower frequency.

1992 ◽  
Vol 82 (6) ◽  
pp. 2308-2327
Author(s):  
Stephen H. Hartzell

Abstract Aftershocks of the 1989 Loma Prieta, California, earthquake are used to estimate site response along the San Francisco Peninsula. A total of 215 shear-wave records from 24 sources and 21 sites are used in a linear inversion for source and site response spectra. The methodology makes no assumptions about the shape of the source spectrum. However, to obtain a stable, unique inverse a Q model and geometrical spreading factor are assumed, as well as a constraint on site response that sets the site response averaged over two specific stations to 1.0. Site responses calculated by this formulation of the problem are compared with other studies in the same region that use different methodologies and / or data. The shear-wave site responses compare favorably with estimates based on an ω2-constrained source model. Comparison with coda amplification factors is not as close, but still favorable considering that the coda values were determined for nearby locations with similar geology, and not the same sites. The degree of agreement between the three methods is encouraging considering the very different assumptions and data used.


1963 ◽  
Vol 53 (5) ◽  
pp. 965-978 ◽  
Author(s):  
David E. Willis

Abstract A comparison of the seismic waves generated by a nuclear explosion and an earthquake is discussed. The epicenter of the earthquake was located within the Nevada Test Site. Both events were recorded at the same station with the same type of equipment. The earthquake waves contained slightly lower frequency than the waves generated by the nuclear shot. The early P phases of the shot had larger amplitudes while the phases after Pg for the earthquake were larger. Seismic waves from collapses were generally found to be composed of lower frequencies than the waves from the original shot. Aftershocks of the Hebgen Lake earthquake were found to generate seismic waves whose frequency content was related to the magnitude of the aftershock. Spectral differences in quarry shot recordings that correlate with source duration times are also discussed.


Author(s):  
Reza Kamgar ◽  
Noorollah Majidi ◽  
Ali Heidari

The nonlinear dynamic analysis provides a more accurate simulation of the structural behavior against earthquakes. On the other hand, this analysis method is time-consuming since the time-step integration schemes are used to calculate the responses of the structure. Wavelet transform is also considered as one of the strong computing tools in studying the properties of the waves. The continuous wavelet transform is a time-frequency study and examines the frequency content of the waves while, the discrete wavelet transform is used to reduce sampling data and also to eliminate the noise of the waves. In this paper, the discrete and continuous wavelet transforms are used to reduce the wave sampling and therefore to reduce the required time for analysis. In this regard, eight near- and far- field earthquakes are studied. The frequency content of the earthquake is investigated by the Fourier spectrum and the continuous wavelet transform. The results show that the first five frequencies for the main earthquakes are similar to those values of earthquakes obtained by wavelet transform. Besides, it is shown that using wavelet transform for the main and decomposed earthquakes indicates that the duration of strong ground motion and the time of dominant frequency occur approximately in the same domain. Finally, it is concluded that the required calculation time reduces to about 80 % with an error less than 6 % when the main earthquake is decomposed by wavelet transform and the approximation waves are used in the nonlinear dynamic analysis.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. R199-R206 ◽  
Author(s):  
Wansoo Ha ◽  
Changsoo Shin

The lack of the low-frequency information in field data prohibits the time- or frequency-domain waveform inversions from recovering large-scale background velocity models. On the other hand, Laplace-domain waveform inversion is less sensitive to the lack of the low frequencies than conventional inversions. In theory, frequency filtering of the seismic signal in the time domain is equivalent to a constant multiplication of the wavefield in the Laplace domain. Because the constant can be retrieved using the source estimation process, the frequency content of the seismic data does not affect the gradient direction of the Laplace-domain waveform inversion. We obtained inversion results of the frequency-filtered field data acquired in the Gulf of Mexico and two synthetic data sets obtained using a first-derivative Gaussian source wavelet and a single-frequency causal sine function. They demonstrated that Laplace-domain inversion yielded consistent results regardless of the frequency content within the seismic data.


Author(s):  
Eben C. Cobb ◽  
Tsu-Chien Cheu ◽  
Jay Hoffman

This paper presents a design methodology to determine the optimal circumferential placement of cylindrical probes upstream of a turbine stage for reduced excitation forces. The potential flow forcing function generated by the probes is characterized by means of a Fourier analysis. A finite difference formulation is used to evaluate the sensitivity of the forcing function to the probe positions. An optimization scheme, based on the linear programming method, uses the sensitivity analysis results to reposition the probes such that the Fourier amplitudes of critical excitation orders are reduced. The results for an example design situation are presented.


Sign in / Sign up

Export Citation Format

Share Document