NIOM Analysis of Earthquake Observation Records in a Rockfill Dam

2022 ◽  
Author(s):  
Hidenori Mogi ◽  
Hideji Kawakami
Keyword(s):  
Propagation Velocity ◽  
Seismic Waves ◽  
Damping Ratio ◽  
Lower Section ◽  
S Wave ◽  
Minimization Method ◽  
Maximum Shear Strain ◽  
P Waves ◽  
Miyagi Prefecture ◽  
Earthquake Observation

ABSTRACT We applied the normalized input–output minimization method (a method developed for the analysis of propagation times in vertical array records) to long-term earthquake observation records from Aratozawa Dam (in Kurihara, Miyagi prefecture, Japan), spanning the period from July 1992 to December 2019 to determine the propagation velocity of seismic waves in the embankment, and investigated changes in soil properties. As a result, we showed that (1) the velocities of S and P waves in the upper section were 449 and 993 m/s, respectively, prior to the strong earthquake motions derived from earthquake records from January 1997 through October 2001, whereas 608 and 1538, respectively, in the lower section, (2) in the Iwate–Miyagi Nairiku earthquake, the S-wave velocity in the upper section decreased to 158 m/s in the principal shock, and (3) in subsequent minor earthquakes the propagation velocity increased more or less in proportion with the logarithm of the number of elapsed days, requiring three years or longer to return to the initial value, (4) although similar changes were observed in the Great East Japan earthquake of 2011, the reduction in propagation velocity that remained after the principal shock was smaller than in the case of the Iwate–Miyagi Nairiku earthquake, and it was judged that there were no large effects on the dam body such as those that occurred in the Iwate–Miyagi Nairiku earthquake, and furthermore (5) in the principal shock of the Iwate–Miyagi Nairiku earthquake, the shear modulus in the upper part of the dam body decreased from 400 to 50 MPa (with a maximum shear strain of 10−3), resulting in more pronounced changes than in the lower section, whereas the damping ratio increased by at least 10% in the lower section during the principal shock of the Iwate–Miyagi Nairiku earthquake, resulting in much greater changes than in the upper section.

scholarly journals Estimation of Geodynamic Properties using Seismic Techniques at a Steel Rolling Factory, Northwestern Gulf of Suez, Egypt

2020 ◽  
Vol 8 (6) ◽  
pp. 1785-1794
Keyword(s):  
Wave Velocity ◽  
Seismic Waves ◽  
Seismic Refraction ◽  
Industrial Area ◽  
P Wave ◽  
Gulf Of Suez ◽  
S Wave ◽  
P Waves ◽  
Steel Rolling ◽  
P Wave Velocity

The objective of the current investigations is to estimate the dynamic geotechnical properties necessary for evaluating the conditions of the subsurface in order to make better decisions for economic and safe designs of the proposed structures at a Steel Rolling Factory, Ataqa Industrial Area, Northwestern Gulf of Suez, Egypt. To achieve this purpose, four seismic refraction profiles were conducted to measure the velocity of primary seismic waves (P-waves) and four profiles were conducted using Multichannel Analysis of Surface Waves (MASW) technique in the same locations of refraction profiles to measure the velocity of shear waves (S-waves). SeisImager/2D Software Package was used in the analysis of the measured data. Data processing and interpretation reflect that the subsurface section in the study area consists of two layers, the first layer is a thin surface layer ranges in thickness from 1 to 4 meters with P-wave velocity ranges from 924 m/s to 1247 m/s and S-wave velocity ranges from 530 m/s to 745 m/s. The second layer has a P-wave velocity ranges from 1277 m/s to 1573 m/s and the S-wave velocity ranges from 684 m/s to 853 m/s. Geotechnical parameters were calculated for both layers. Since elastic moduli such as Poisson’s ratio, shear modulus, Young’s modulus, and bulk’s modulus were calculated. Competence scales such as material index, stress ratio, concentration index, and density gradient were calculated also. In addition, the ultimate and allowable bearing capacities

Corner frequencies of P and S waves and models of earthquake sources

1973 ◽  
Vol 63 (6-1) ◽  
pp. 2091-2104 ◽  
Author(s):  
Peter Molnar ◽  
Brian E. Tucker ◽  
James N. Brune
Keyword(s):  
Seismic Waves ◽  
Wave Spectra ◽  
S Wave ◽  
P Waves ◽  
S Waves ◽  
Earthquake Sources ◽  
Dislocation Models ◽  
San Fernando

Abstract P- and S-wave spectra of 144 aftershocks 12≦M≦412 of the February 9, 1971 San Fernando earthquake corroborate previous work showing that the corner frequencies for P waves in general are greater than those for S waves. This observation is consistent not only with models that treat earthquakes as volume sources, but also with physically reasonable dislocation models for which (1) the source is approximately equidimensional, (2) both the duration of slip at each point on the fault and the time for the ruptured area to develop are not long compared with the time for seismic waves to cross the ruptured area, and (3) much of the source radiates essentially simultaneously. There may be other physically reasonable dislocation models compatible with the observations. Savage's calculations indicate that models that involve propagating dislocations on long thin faults are not adequate for describing most moderate and small earthquakes studied.

Geometrical spreading in a layered transversely isotropic medium with vertical symmetry axis

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad
Keyword(s):  
Ray Tracing ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Analytical Approximation ◽  
Geometrical Spreading ◽  
S Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
Amplitude Versus

Compensation for geometrical spreading is important in prestack Kirchhoff migration and in amplitude versus offset/amplitude versus angle (AVO/AVA) analysis of seismic data. We present equations for the relative geometrical spreading of reflected and transmitted P‐ and S‐wave in horizontally layered transversely isotropic media with vertical symmetry axis (VTI). We show that relatively simple expressions are obtained when the geometrical spreading is expressed in terms of group velocities. In weakly anisotropic media, we obtain simple expressions also in terms of phase velocities. Also, we derive analytical equations for geometrical spreading based on the nonhyperbolic traveltime formula of Tsvankin and Thomsen, such that the geometrical spreading can be expressed in terms of the parameters used in time processing of seismic data. Comparison with numerical ray tracing demonstrates that the weak anisotropy approximation to geometrical spreading is accurate for P‐waves. It is less accurate for SV‐waves, but has qualitatively the correct form. For P waves, the nonhyperbolic equation for geometrical spreading compares favorably with ray‐tracing results for offset‐depth ratios less than five. For SV‐waves, the analytical approximation is accurate only at small offsets, and breaks down at offset‐depth ratios less than unity. The numerical results are in agreement with the range of validity for the nonhyperbolic traveltime equations.

A periodic seismic isolation foundation with an extremely broad low-frequency attenuation zone: Theoretical analysis and experimental verification

2021 ◽  
pp. 136943322110646
Author(s):  
Peng Zhou ◽  
Shui Wan ◽  
Xiao Wang ◽  
Yingbo Zhu ◽  
Muyun Huang
Keyword(s):  
Seismic Isolation ◽  
Seismic Waves ◽  
Periodic Structures ◽  
Finite Size ◽  
Low Frequency ◽  
Bridge Structure ◽  
Engineering Structures ◽  
P Waves ◽  
New Type ◽  
Dominant Frequencies

The attenuation zones (AZs) of periodic structures can be used for seismic isolation design. To cover the dominant frequencies of more seismic waves, this paper proposes a new type of periodic isolation foundation (PIF) with an extremely wide low-frequency AZ of 3.31 Hz–17.01 Hz composed of optimized unit A with a wide AZ and optimized unit B with a low-frequency AZ. The two kinds of optimized units are obtained by topology optimization on the smallest periodic unit with the coupled finite element-genetic algorithm (GA) methodology. The transmission spectra of shear waves and P-waves through the proposed PIF of finite size are calculated, and the results show that the AZ of the PIF is approximately the superposition of the AZs of the two kinds of optimized units. Additionally, shake tests on a scale PIF specimen are performed to verify the attenuation performance for elastic waves within the designed AZs. Furthermore, numerical simulations show that the acceleration responses of the bridge structure with the proposed PIF are attenuated significantly compared to those with a concrete foundation under the action of different seismic waves. Therefore, the newly proposed PIF is a promising option for the reduction of seismic effects in engineering structures.

Diffraction of elastic waves by a fluid-filled crack in a fluid-saturated poroelastic half-space

2021 ◽  
Author(s):  
Zhongxian Liu ◽  
Jiaqiao Liu ◽  
Sibo Meng ◽  
Xiaojian Sun
Keyword(s):  
Half Space ◽  
Seismic Waves ◽  
Incident Angle ◽  
Excitation Frequency ◽  
Vertical Displacement ◽  
Crack Model ◽  
P Waves ◽  
Amplification Effect ◽  
Resonance Frequencies ◽  
Angle Crack

Summary An indirect boundary element method (IBEM) is developed to model the two-dimensional (2D) diffraction of seismic waves by a fluid-filled crack in a fluid-saturated poroelastic half-space, using Green's functions computed considering the distributed loads, flow, and fluid characteristics. The influence of the fluid-filled crack on the diffraction characteristics is investigated by analyzing key parameters, such as the excitation frequency, incident angle, crack width and depth, and medium porosity. The results for the fluid-filled crack model are compared to those for the fluid-free crack model under the same conditions. The numerical results demonstrate that the fluid-filled crack has a significant amplification effect on the surface displacements, and that the effect of the depth of the fluid-filled crack is more complex compared to the influence of other parameters. The resonance diffraction generates an amplification effect in the case of normally incident P waves. Furthermore, the horizontal and vertical displacement amplitudes reach 4.2 and 14.1, respectively. In the corresponding case of the fluid-free crack, the vertical displacement amplitude is only equal to 4.1, indicating the amplification effect of the fluid in the crack. Conversely, for normally incident SV waves at certain resonance frequencies, the displacement amplitudes above a fluid-filled crack may be lower than the displacement amplitudes observed in the corresponding case of a fluid-free crack.

scholarly journals Simulation of seismic waves at the earth's crust (brittle–ductile transition) based on the Burgers model

Solid Earth ◽  
2014 ◽  
Vol 5 (2) ◽  
pp. 1001-1010 ◽  
Author(s):  
J. M. Carcione ◽  
F. Poletto ◽  
B. Farina ◽  
A. Craglietto
Keyword(s):  
Seismic Waves ◽  
In Situ Stress ◽  
Seismic Attenuation ◽  
Earth’S Crust ◽  
S Wave ◽  
Fourier Pseudospectral Method ◽  
Stress Criterion ◽  
Brittle Ductile Transition ◽  
Earth's Crust

Abstract. The earth's crust presents two dissimilar rheological behaviors depending on the in situ stress-temperature conditions. The upper, cooler part is brittle, while deeper zones are ductile. Seismic waves may reveal the presence of the transition but a proper characterization is required. We first obtain a stress–strain relation, including the effects of shear seismic attenuation and ductility due to shear deformations and plastic flow. The anelastic behavior is based on the Burgers mechanical model to describe the effects of seismic attenuation and steady-state creep flow. The shear Lamé constant of the brittle and ductile media depends on the in situ stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P and S wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the P−S and SH equations of motion recast in the velocity-stress formulation, including memory variables to avoid the computation of time convolutions. The equations correspond to isotropic anelastic and inhomogeneous media and are solved by a direct grid method based on the Runge–Kutta time stepping technique and the Fourier pseudospectral method. The algorithm is tested with success against known analytical solutions for different shear viscosities. A realistic example illustrates the computation of surface and reverse-VSP synthetic seismograms in the presence of an abrupt brittle–ductile transition.

Numerical simulation of azimuthal acoustic logging in a borehole penetrating a rock formation boundary

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  
Keyword(s):  
Domain Model ◽  
P Wave ◽  
S Wave ◽  
Angle Variation ◽  
P Waves ◽  
Acoustic Logging ◽  
S Waves ◽  
Rock Formation ◽  
Logging Tool ◽  
Difference Time

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.

Elastic Parameters from Compressional and Shear Wave Tomographic Survey: A Case Study from Kuala Lumpur, Malaysia

2017 ◽  
Vol 22 (4) ◽  
pp. 427-434
Author(s):  
Julius K. von Ketelhodt ◽  
Thomas Fechner ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim
Keyword(s):  
Wave Train ◽  
Kuala Lumpur ◽  
S Wave ◽  
Elastic Parameters ◽  
High Quality ◽  
P Waves ◽  
Receiver Array ◽  
Wave Tomography ◽  
Karstic Limestone

An integrated P- and S-wave cross-borehole tomographic survey was performed in the city center of Kuala Lumpur, Malaysia, with the aim of exploring a karstic limestone area near an area that previously encountered cavities. Horizontally polarized shear waves were generated with two opposing, perpendicular strike directions and recorded with a multi-level, three-component receiver array. This allowed a high quality picking of the traveltimes, whereby the wave train reverses at the time of the S-wave arrival. In addition, high quality sparker generated P-waves were recorded. The P- and S-wave traveltimes were used to invert for two co-located tomograms. These tomograms enabled a better interpretation capability than a P- or S-wave tomogram on its own. The tomograms enabled the calculation of the elastic parameters, i.e., P- to S-wave velocity (Vp/Vs) ratio, Poisson's ratio, bulk modulus, Young's modulus and the shear modulus, on a 2D surface between the boreholes. This further aided the interpretation, as areas with limited traveltime accuracy and thus, an increase in tomographic error, could be easily identified, and the extent of a large cavity could be estimated. The interpretation of the tomograms was constrained by two additional boreholes, which provided more confidence on the delineation and location of cavities at depths. The survey shows the benefit of co-locating P- and S-wave tomography surveys.

On the Green function of the Lord–Shulman thermoelasticity equations

2019 ◽  
Vol 220 (1) ◽  
pp. 393-403 ◽  
Author(s):  
Zhi-Wei Wang ◽  
Li-Yun Fu ◽  
Jia Wei ◽  
Wanting Hou ◽  
Jing Ba ◽  
...  
Keyword(s):  
Green’S Function ◽  
Green's Function ◽  
Parabolic Type ◽  
Second Order ◽  
P Wave ◽  
Thermal Mode ◽  
S Wave ◽  
Order Tensor ◽  
P Waves ◽  
Thermal Conductivities

SUMMARY Thermoelasticity extends the classical elastic theory by coupling the fields of particle displacement and temperature. The classical theory of thermoelasticity, based on a parabolic-type heat-conduction equation, is characteristic of an unphysical behaviour of thermoelastic waves with discontinuities and infinite velocities as a function of frequency. A better physical system of equations incorporates a relaxation term into the heat equation; the equations predict three propagation modes, namely, a fast P wave (E wave), a slow thermal P wave (T wave), and a shear wave (S wave). We formulate a second-order tensor Green's function based on the Fourier transform of the thermodynamic equations. It is the displacement–temperature solution to a point (elastic or heat) source. The snapshots, obtained with the derived second-order tensor Green's function, show that the elastic and thermal P modes are dispersive and lossy, which is confirmed by a plane-wave analysis. These modes have similar characteristics of the fast and slow P waves of poroelasticity. Particularly, the thermal mode is diffusive at low thermal conductivities and becomes wave-like for high thermal conductivities.

Vertical Vibration Control of Series Isolation System under Near-Fault Earthquake

2014 ◽  
Vol 638-640 ◽  
pp. 1952-1955
Author(s):  
Q. Rong ◽  
Yan Sheng ◽  
Shi Xin Liu
Keyword(s):  
Seismic Waves ◽  
Damping Ratio ◽  
Vertical Vibration ◽  
Isolation Effect ◽  
Earthquake Intensity ◽  
Isolation System ◽  
Vertical Stiffness ◽  
Isolation Layer ◽  
Near Fault ◽  
Composite Disk

Series isolation system consists of rubber isolation bearings and composite disk springs, determination method of vertical stiffness and vertical damping of isolation layer is given. Entering the near-fault vertical seismic waves, the affect of isolation layer parameters and earthquake intensity on the isolation effect is studied. Studies have shown that the vertical isolation effect increases with the increase of vertical damping ratio. When the damping ratio reaches a certain value, the isolated effect leveles off. When calculating model is adopted as the hierarchical model, vertical isolation effect has nothing to do with the increases of earthquake intensity.

Sign in / Sign up

Export Citation Format

Share Document