On rigid foundations subjected to seismic waves

In the analysis of structural foundations for seismic loads, it is customary to distinguish two types of soil-structure interaction effect: kinematic interaction (or wave passage), and inertial interaction. The former refers to the phenomenon of wave scattering, which occurs because the foundation is much stiffer than the surrounding soil and cannot accommodate to its distortions. Inertial interaction, on the other hand, is caused by feedback of kinetic energy of the structure into the soil. This paper is concerned only with the first phenomenon. The rigorous analysis of rigid, embedded foundations subjected to seismic disturbances requires, in general, substantial computational effort. Indeed, a typical analysis would normally require models with finite elements and/or boundary elements. Although such methods may be used to find an accurate solution to the problem of kinematic interaction, their use is not always warranted, given the many uncertainties involved and the multitude of assumptions that must be considered. Hence, approximate solutions are attractive for this problem. One such approximate method is the remarkably simple algorithm proposed by Iguchi.³ This paper presents first an appraisal of this method by way of a comparison with accurate numerical solutions for cylindrical foundations; next the algorithm is applied to rectangular (prismatic) foundations. It is found that Iguchi's method gives results that are adequate for engineering purposes, even if not entirely accurate.     

Sensitivity of time-lapse seismic to reservoir stress path

The change in reservoir pore pressure due to the production of hydrocarbons leads to anisotropic changes in the stress field acting on the reservoir. Reservoir stress path is defined as the ratio of the change in effective horizontal stress to the change in effective vertical stress from the initial reservoir conditions, and strongly influences the depletion‐induced compaction behaviour of the reservoir. Seismic velocities in sandstones vary with stress due to the presence of stress‐sensitive regions within the rock, such as grain boundaries, microcracks, fractures, etc. Since the response of any microcracks and grain boundaries to a change in stress depends on their orientation relative to the principal stress axes, elastic‐wave velocities are sensitive to reservoir stress path. The vertical P‐ and S‐wave velocities, the small‐offset P‐ and SV‐wave normal‐moveout (NMO) velocities, and the P‐wave amplitude‐versus‐offset (AVO) are sensitive to different combinations of vertical and horizontal stress. The relationships between these quantities and the change in stress can be calibrated using a repeat seismic, sonic log, checkshot or vertical seismic profile (VSP) at the location of a well at which the change in reservoir pressure has been measured. Alternatively, the variation of velocity with azimuth and distance from the borehole, obtained by dipole radial profiling, can be used. Having calibrated these relationships, the theory allows the reservoir stress path to be monitored using time‐lapse seismic by combining changes in the vertical P‐wave impedance, changes in the P‐wave NMO and AVO behaviour, and changes in the S‐wave impedance.     

Sensitivity of seismic hazard evaluations to uncertainties determined from seismic source characterization

The sensitivity and overall uncertainty in peak ground acceleration (PGA)estimates have been calculated for the city of Tabriz, northwestern Iran byusing a specific randomized blocks design. Eight seismic hazard models andparameters with randomly selected uncertainties at two levels have beenconsidered and then a linear model between predicted PGA at a givenprobability level and the uncertainties has been performed. The inputmodels and parameters are those related to the attenuation, magnituderupture-length and recurrence relationships with their uncertainties.Application of this procedure to the studied area indicates that effects ofthe simultaneous variation of all eight input models and parameters on thesensitivity of the seismic hazard can be investigated with a decreasingnumber of computations for all possible combinations at a fixed annualprobability. The results show that the choice of a mathematical model ofthe source mechanism, attenuation relationships and the definition ofseismic parameters are most critical in estimating the sensitivity of seismichazard evaluation, in particular at low levels of probability of exceedance.The overall uncertainty in the expected PGA for an annual probability of0.0021 (10% exceedence in 50 yr) is expressed by a coefficient ofvariation (CV) of about 34% at 68% confidence level for a distance ofabout 5km from the field of the major faults. The CV will decrease withincreasing site-source distance and remains constant, CV = 15%, fordistances larger than 15 km. Finally, treating alternative models on theoverall uncertainty are investigated by additional outliers in input decision.     

Studies on seismic waves

The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous me-dia, anisotropie and porous media, surface wave and seismic wave inversion, and seismic wave study in prospect-ing and logging problems. Important projects in the current studies on seismic wave is suggested as the develop-ment of high efficient numerical methods, and applying them to the studies of excitation and propagation of seis-mic waves in complex media and strong ground motion, which will form a foundation for refined earthquake haz-ard analysis and prediction.     

Studies on seismic waves

The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous media, anisotropic and porous media, surface wave and seismic wave inversion, and seismic wave study in prospecting and logging problems. Important projects in the current studies on seismic wave is suggested as the development of high efficient numerical methods, and applying them to the studies of excitation and propagation of seismic waves in complex media and strong ground motion, which will form a foundation for refined earthquake hazard analysis and prediction.     

Amplitudes of explosion-generated seismic waves

Summary Records of underwater explosions in Iceland in 1959 and 1960, mainly performed for crustal studies, have been investigated with regard to amplitudes. The amplitudes ofP2 in the 1959 explosions and of the first arrivingP waves in the 1960 explosions were found to be proportional to the first power of the charge weight, whereas amplitudes ofP1 andS1 in the 1959 explosions were proportional to the 3/4-power of the charge up to 200 kg, but increased much slower for larger charges. The influence of water depth of shot point on the charge weights, required to obtain a certain amplitude, has been determined and it was found that the logarithm of the charge weight has a linear relation to water depth down to about 8 meters. The amplitudes of first arrivingP waves decrease as the inverse 2.2-power of the distance up to about 30 km. For greater distances an exponential decrease of the form (const./) exp (–) is valid. ForP2 waves with a frequency of 10 cps we found =0.027±0.003 km^–1 for a profile across central Iceland and =0.009±0.004 km^–1 for a profile in the western part of Iceland.

---

Zusammenfassung In den Jahren 1959 und 1960 wurden auf Island Unterwasser-Explosionen ausgelöst, deren Hauptziel Untersuchung der Erdkruste waren. Die Aufzeichnungen der Explosionen werden hier in Hinsicht auf die Amplituden der Wellen untersucht. DieP2-Amplituden in den Explosionen von 1959 und die Amplituden der zuerst ankommendenP-Wellen in den Explosionen von 1960 waren proportional der ersten Potenz des Ladungsgewichtes, während diePl-undSl-Amplituden in den Explosionen von 1959 proportional der 3/4-Potenz des Ladungsgewichtes bis 200 kg waren. Für höhere Ladungsgewichte wuchsen die Amplituden langsamer an. Es wurde der Einfluss der Tiefe des Schusspunktes unterhalb des Wasserspiegels auf das für eine gewisse Amplitude benötigte Ladungsgewicht bestimmt. Es ergab sich dabei, dass der Logarithmus des Ladungsgewichtes in einer linearen Beziehung zur Tiefe, bis etwa 8 m, steht. Die Amplituden der zuerst ankommendenP-Wellen sind umgekehrt proportional der 2.2-Potenz der Entfernung, bis etwa 30 km. Für grössere Entfernungen ist die Abnahme der Amplituden exponential und kann durch (Konst./) exp (–) dargestellt werden. FürP2 mit einer Frequenz von 10 cps haben wir =0.027±0.003 km^–1 für ein Profil in Zentral-Island und =0.009 ± 0.004 km^–1 für ein Profil in West-Island gefunden.

     

Sensitivity of seismic structural response to interpretation of soils data

This paper investigates the sensitivity of the predicted seismic response of buildings in a PWR nuclear power station to the potential changes in the techniques and methods of interpreting soil data that have occurred over the last decade. The investigation is based on the soil-structure interaction (SSI) response of a typical PWR reactor building on a soft site during a seismic event. The current techniques and methods of interpretation of soils data tend to lead to a stiffer site with lower soil material damping than the earlier techniques. This leads to an increase in the SSI natural frequencies of typical buildings and an increase in its seismic response. This increase in the seismic response could put into question any seismic design based on seismic loads derived using the previously accepted generic soil data. The paper concludes with a recommendation for further consideration of the proposed departure from the previously accepted soil data.     

层状围岩隧道交叉结构地震方向敏感性分析

依托某层状围岩深埋铁路隧道工程,以主线单线隧道与横通道正交结构为研究对象,通过ABAQUS建立有限元分析模型.选取12条地震波数据,模拟基岩SV波地震作用,调整地震动峰值加速度和地震动入射方向,对研究结构随地震动输入方向的敏感性进行分析.研究分析得到:二次衬砌在地震入射角度为30°～45°时对地震动方向变化敏感性较低,...     

Torsional response of structures to obliquely incident seismic sh waves

A study is made of the torsional response of an elastic structure placed on a rigid circular foundation supported on an elastic half-space and subjected to the action of obliquely incident plane SH waves. The problem is solved by considering first the steady-state response of a massless rigid foundation excited externally by a harmonic torque and through the soil by an obliquely incident plane SH wave. In a second stage the coupling between the structure and the soil is considered to obtain the torsional response at the base and top of the superstructure. The results obtained indicate a range of conditions under which the torsional effects will be most pronounced.     

论可能存在的地幔柱地震导波

假想的地幔热地幔柱预期地震波速低,因而可支持类似于断层带地震导波、光纤波和海洋声音固定和测距通道声波的弹性导波传播。地幔柱导波在波导中对几何复杂性并不敏感,其频散会使其在地震记录中与众不同并提供可以补充地震层析成像波导结构的有关信息。检测这类波将成为地幔柱存在新一类强有力的证据。嵌于无限介质中的圆柱形通道支持两类轴对称弹性波振型,即扭转振型和纵向—径向振型。扭转振型有与柱体表面横切的直线质点运动。纵向—径向振型在含有圆柱体轴的平面上有椭圆形质点运动,在轴附近有反向运动。椭圆形质点运动的方向随距离与轴反向:基本振型为一倍,第一谐波为两倍,依此类推。每种振型只有在超过其截止频率时才存在,相速度和群速度在无限介质中等于剪切波速。在高频率时两种速度在通道中接近剪切波速。所有振型在它们的群速度中均有最小值,在地震记录上产生艾里相位。对于剪切波速有百分之几的反差尽管对于地球的热地幔柱而言是现实的,但最大的信号是反向频散的,并有约0.1～1Hz和15～30s持续时间的优势频率。可观测到的地幔柱波至少有两个可能的来源:(1)深地幔中地幔柱与高振幅核震相焦散线的相交;(2)向下传播的导波在核幔边界ScS状的反射。最近广泛部署的宽带地震仪使搜索这些波成为可能。     

Dynamic response of rectangular foundations to obliquely incident seismic waves

A study is made of the harmonic response of a rigid massless rectangular foundation bonded to an elastic half-space and subjected to the action of both external forces and obliquely incident plane seismic waves. The associated mixed boundary value problem is discretized and solved numerically. The results obtained indicate that the angle of incidence of the seismic wave has a marked effect on the nature and magnitude of the foundation response.     

On surface seismic waves formation

Summary The interpretation of surface seismic waves records is rather complicated as they include a superposition of oscillations of the fundamental mode and higher modes. Besides recorded oscillations depend on spectral characteristics of motions in earthquakes sources. The consideration of these problems is based on results of surface waves two-dimensional modelling [1]³. Some physical ideas about their formation deals with the change of the nature of the oscillation propagating with dispersion. This report represents a condensate of several independent works. , , . , . , . () . .Scientific communication presented to the IASPEI Assembly, Madrid, 1969.     

An attenuating structure beneath the Aso Caldera determined from the propagation of seismic waves

The attenuation of amplitude is seen in seismic waves which pass through the central region of the Aso caldera, in Kyushu, Japan. It is also recognized from spectral analysis of seismic waves that the higher frequencies of the P-wave are reduced in the waves which pass through the central region of the caldera. It is shown that the relative attenuation increases remarkably for the frequency range of 5 to 10 Hz. The specific attenuation factor Q of the P-wave train is about 100. From the surface projection of the ray paths with low Q values through the Aso caldera to each station, the attenuating region is located beneath the center of the caldera, extending to the north of the central cones. In conjunction with the low Q value of the P-wave and the decreases of S-wave amplitudes, the relative P-wave residual times have comparatively large values for seismic waves passing through the central region beneath the caldera. In order to attempt to provide additional information on the depth configuration of the attenuating material, the ray paths of P-wave's first arrivals are located in three-dimensional space. It indicates that the low-velocity material is located beneath the center of the caldera at depths of about 6 to 9 km. However, lowvelocity anomalies above the depth of 6 km and below the depth of 15 km were not able to be detected, because most of the available seismic ray paths had crossed the caldera at depths of about 6 to 15 km. Furthermore, the relative residual times have numerous errors resulting from incorrect hypocenter locations, origin times, inhomogeneities in the structure and uncertainty of the velocity structure. At shallow depths in the Aso caldera, refraction or reflection studies are required for an accurate estimate of the structure and more detailed properties of the attenuating material.     

Continental dispersion of seismic surface waves

Summary Continental dispersion curves for Love and Rayleigh waves have been determined by means of the Uppsala seismograph records for three earthquakes near the south coast of Turkey in April, 1957. To the authors' knowledge this is the first time dispersion curves have been determined for a relatively homogenous path across the whole of continental Europe and covering a good range of periods, namely from 3 to 86 sec for Love waves and from 5 to 61 sec for Rayleigh waves. The dispersion curves obtained are relevant to determinations of the average structure along the path traversed in particular, and to determinations of crustal structure of the European continent in general.

---

Zusammenfassung Kontinentale Dispersionskurven für Love-und Rayleigh-Wellen wurden vermittels Registrierungen an der Station Uppsala von drei Erdbeben in der Nähe der Südküste der Türkei im April 1957 bestimmt. Soweit es den Autoren bekannt ist, wurden damit erstmalig Dispersionskurven für einen relativ homogenen Weg über das gesamte kontinentale Europa mit einem beträchtlichen Periodenbereich, der 3 bis 86 Sek für Love-Wellen und 5 bis 61 Sek für Rayleigh-Wellen umfasst, ausgearbeitet. Die erhaltenen Dispersionskurven haben ihre Bedeutung in der Festlegung der mittleren Krustenstruktur entlang dem durchlaufenen Wellenweg und überhaupt in der Festlegung der Krustenstruktur im Bereich des europäischen Kontinents.

     

Soil amplification of plane seismic waves

Amplitudes of surface particle velocities are calculated when time-harmonic seismic waves of uniform amplitude are incident upon an arbitrary stratified elastic soil layer from the underlying bedrock. Whereas previous workers have mainly treated normally incident S waves, we allow the waves to be of SV, P, or SH types and to have arbitrary angles of incidence. Following standard practice the problem is set up as a matrix differential system, but in such a way that incident SV and P waves may be treated together (the system for SH decoupling). Though complicated, the 4 × 4 SVP system has considerable structure which is elucidated in Appendices 1 and 2. These results, though not altogether new, are of independent interest, and are gathered together in concise form for reference. The theory for low- and for high-frequency approximations is given. The main results of the work are illustrated by two numerical examples: Model 1 where the soil layer is homogeneous; and Model 2 where the soil layer has a linear velocity profile.     

Sensitivity of reflection seismic data to oil-column height in high-porosity sandstones

Tuning is the effect of interference between the reflections from the top and bottom of a thin layer on the amplitude of the composite reflection. For a homogeneous sandstone reservoir containing an oil column overlying brine, interference between the reflection from the top reservoir and the oil/water contact is a function of the height of the oil column. If the properties of the sandstone do not vary across the oil/water contact, the SS, PS and SP reflection coefficients from the oil/water contact are small in comparison to the PP reflection coefficient. This allows analytic expressions for the effective PP and PS reflection coefficients from the reservoir to be derived that include all P‐wave multiples within the oil column. For a given source/receiver offset, the component of the wavevector inside the oil column normal to the interface is larger for the PPPP reflection than for the PPPS reflection, due to the asymmetry in the raypath for the PPPS reflection. The PPPS reflection is therefore useful for determining oil‐column heights larger than that discriminated by the PPPP reflection, especially when used at wider offsets. A convenient classification of the AVO response of hydrocarbon‐bearing sandstone reservoirs overlain by shale is the scheme of Rutherford and Williams. Class 1 sands have higher acoustic impedance for normal incidence than the overlying shale, Class 2 sands have nearly the same acoustic impedance as the shale and Class 3 sands have lower acoustic impedance. Synthetic shot gathers calculated for these three classes as a function of oil‐column height show that a combination of the PPPP and the PPPS amplitudes can be plotted as a tuning trajectory, which can be used to determine the oil‐column height. This method is most sensitive for reservoirs that belong to AVO classes 1 and 2, and therefore may be useful in AVO analysis of Class 1 and 2 reservoirs where the traditional AVO indicators (developed for Class 3 reservoirs) do not work very well. These results demonstrate the usefulness of shear waves recorded in the marine environment at wide offsets.     

地震波作用下的埋地燃气管道动力响应研究

埋地管道是现代城市的重要生命线工程,因其埋于地下,一旦遭遇地震,管道就有可能发生破坏,导致燃气泄露,可能引发火灾、爆炸等次生灾害,对人们的生命和财产造成严重威胁,因此,如何提高埋地燃气管道的抗震性能,是当前需要迫切解决的问题.认为在管道轴向上和横向上均与土体之间存在相对位移,结合拟静力法进行了地震波作用下埋地管道的动力...     

Simplified longitudinal seismic response of tunnels linings subjected to surface waves

This paper presents an engineering approach for analysing the longitudinal behaviour of tunnels subjected to earthquakes. The tunnel is modelled as a Timoshenko beam connected to the far soil by means of continuous elastic support (Winkler model). Seismic free-field inputs, such as those caused by surface waves travelling parallel to the tunnel axis, were imposed at the base of the springs of the Winkler model, generating bending moments and shear forces on the cross-sections of the tunnel. Closed-form expressions of the tunnel displacements, shear forces, and bending moments were determined at any tunnel section in terms of the seismic excitation, tunnel geometry and material properties, and subgrade reaction modulus of the soil. A dimensional analysis was carried out to ascertain directly the maximum tunnel displacement, bending moment, and shear force.     

Spectra of seismic waves of Matsushiro micro-earthquakes

Summary The spectra of Matsushiro micro-earthquakes were studied using a spectrum analyser and discussed in relation to their magnitude and the epicentral distance. It was found that the spectra of seismic waves by Matsushiro micro-earthquakes did not depend upon the epicentral distance since the sonagrams showed a peak around 5 cps at Matsushiro as well as Suzaka. The spectrum also failed to show any distinct peak for micro-earthquakes of magnitude lower than 1.0 due to the decrease of signal to noise ratio.