砂土液化引起大位移对地下管道影响的非线性分析

地下管线是生命线工程的主要部分,已经成为现代工农业生产和城镇生活的大动脉。已有震害调查表明,饱和砂土液化引起的地基大变形(侧向变形和沉降)是导致强震区生命线工程震害的主要原因。采用三维非线性有限差分分析方法来研究砂土液化引起的大位移对地下管道的破坏特征,分析砂土液化的斜坡变形特征、孔隙水的演化过程。结果表明,砂土液化引起的大位移对地下管道有破坏作用,导致管道变形规律与其斜坡的位移规律相同,地下管线的变形随着振动频率和幅值的增加其非线性增大。     

地震作用下反倾岩质斜坡动力响应规律及频谱特征研究

以四川昌都贡扎滑坡为原型,设计制作了1:1000的振动台试验模型,开展了包含软弱岩性组合和贯通性结构面的反倾岩质斜坡的振动台模型试验.考虑场地工程地质条件和场地特性,将富含多种频率的人工合成基岩地震波作为动力输入条件,探讨了不同概率水准人工合成地震波作用下反倾岩质斜坡的动力响应规律,并与输入正弦波的放大效应规律进行对比...     

Love waves in inhomogeneous and anisotropic earth — I

Summary The frequency equation is derived for the propagation of Love waves in the earth's crust, composed of transversely isotropic layers and overlying anisotropic and inhomogeneous mantle. The exact boundary value problem is solved for a single layer and extended to multilayered media by generalizing theHaskell's technique. In fact the problem of deriving the frequency equation has been reduced to finding out the solution of the equation of motion subject to the appropriate boundary conditions. To illustrate the method, the author has derived frequency equations of Love waves for linear, exponential and generalized power law variation of vertical shear wave velocity with depth in the half space overlain by transversely isotropic inhomogeneous stratum.     

线性波浪加载下海底斜坡失稳机制的数值分析

基于大型有限元软件ABAQUS中的荷载模块,添加一阶波浪力载荷模式,并结合强度折减技术,实现波浪力作用下海底斜坡稳定性与失稳机制的弹塑性有限元数值分析。引入典型算例,利用先前提出的波浪荷载下海底斜坡稳定性的极限分析上限方法开展数值解的对比验证;在此基础上,通过深入地变动参数比较分析,探讨不同波长、波高和水深等波浪参数对计算结果的影响以及波浪力影响下海底斜坡潜在滑动面的变化规律,获得波浪荷载下海底斜坡失稳滑动机制的初步认识。     

Wave-induced steady streaming,mass transport and net sediment transport in rough turbulent ocean bottom boundary layers

The interaction between two important mechanisms which causes streaming has been investigated by numerical simulations of the seabed boundary layer beneath both sinusoidal waves and Stokes second order waves, as well as horizontally uniform bottom boundary layers with asymmetric forcing. These two mechanisms are streaming caused by turbulence asymmetry in successive wave half-cycles (beneath asymmetric forcing), and streaming caused by the presence of a vertical wave velocity within the seabed boundary layer as earlier explained by Longuet-Higgins. The effect of wave asymmetry, wave length to water depth ratio, and bottom roughness have been investigated for realistic physical situations. The streaming induced sediment dynamics near the ocean bottom has been investigated; both the resulting suspended load and bedload are presented. Finally, the mass transport (wave-averaged Lagrangian velocity) has been studied for a range of wave conditions. The streaming velocities beneath sinusoidal waves (Longuet-Higgins streaming) is always in the direction of wave propagation, while the streaming velocities in horizontally uniform boundary layers with asymmetric forcing are always negative. Thus the effect of asymmetry in second order Stokes waves is either to reduce the streaming velocity in the direction of wave propagation, or, for long waves relative to the water depth, to induce a streaming velocity against the direction of wave propagation. It appears that the Longuet-Higgins streaming decreases as the wave length increases for a given water depth, and the effect of wave asymmetry can dominate, leading to a steady streaming against the wave propagation. Furthermore, the asymmetry of second order Stokes waves reduces the mass transport (wave-averaged Lagrangian velocity) as compared with sinusoidal waves. The boundary layer streaming leads to a wave-averaged transport of suspended sediments and bedload in the direction of wave propagation.     

Wave propagation in two layered medium under initial stresses

Summary In this paper, the frequency equation for phase velocity of waves propagated in a laminated medium consisting of two eleastic layers of finite thickness under initial stresses, has been obtained. It has been shown that when wave length becomes very small compared to the thickness of each layer, the wave approaches two Rayleigh waves at the two outer surfaces with the possibility of Stoneley waves at the interface. The propagation ofSH-waves in the composite medium under initial stresses has also been discussed. A particular case has been taken to find the velocity of Love wave in the homogeneous half space under initial compressive stresses.Biot's incremental deformation theory has been used.     

Crustal Structure of the Korean Peninsula Using Surface Wave Dispersion and Numerical Modeling

Surface wave dispersion is studied to obtain the 1-D average velocity structure of the crust in the Korean Peninsula by inverting group- and phase-velocities jointly. Group velocities of short-period Rayleigh and Love waves are obtained from cross-correlations of seismic noise. Multiple-filter analysis is used to extract the group velocities at periods between 0.5 and 20 s. Phase velocities of Rayleigh waves in 10- and 50-s periods are obtained by applying the two-station method to teleseismic data. Dispersion curves of all group and phase velocities are jointly inverted for the 1-D average model of the Korean Peninsula. The resultant model from surface wave analysis can be used as an initial model for numerical modeling of observations of North Korean events for a velocity model appropriated to the Korean Peninsula. The iterative process is focused especially on the surface sedimentary layer in the numerical modeling. The final model, modified by numerical modeling from the initial model, indicates that the crust shear wave velocity increases with depth from 2.16 km/s for a 2-km-thick surface sedimentary layer to 3.79 km/s at a Moho depth of 33 km, and the upper mantle has a velocity of 4.70 km/s.     

基于多波束地形测量的上荆江典型河段沙波形态特征分析

沙波形态影响水流结构、泥沙输移及动床阻力。本研究采用多波束测深系统首次精细测量了上荆江典型河段的床面地形,采用改进后的沙波形态量化算法统计了各类沙波形态参数,分析了不同水流强度下沙波形态的变化特征。计算结果表明:(1)测量河段小型与大型沙波的平均波高分别为0.16～0.81和0.96～2.31 m,波长分别为13～27和16～41m;沙波尺度相较于水深较小,小型与大型沙波的波高分别不超过水深的0.045和0.150倍;(2)沙波背流面坡度基本不超过14°,小于泥沙水下休止角,其与陡度之间的关系可以用线性方程描述;(3)中洪水流量对沙波形态尺度的塑造作用强于枯水流量,且对浅水区大型沙波形态尺度的塑造作用强于深水区。本研究量化了天然河流的沙波形态,较好地反映了沙波形态特征,能为大型冲积河流沙波形态的量化及特征参数的统计分析提供参考。     

来自海底高速层径向波的理论地震图研究

本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.     

A note on specific variability of long surface gravity waves and drag coefficient in coastal upwelling zone

In this paper we solve analytically wave kinematic equations and the wave energy transport equation, for basic long surface gravity wave in the coastal upwelling zone. UsingGent andTaylor's (1978) parameterization of drag coefficient (which includes interaction between long surface waves and the air flow) we find variability of this coefficient due to wave amplification and refraction caused by specific surface water current in the region. The drag coefficient grows towards the shore. The growth is faster for stronger current. When the angle between waves and the current is less than 90° the growth is mainly connected with the waves steepness, but when the angle is larger, it is caused by relative growth of the wave phase velocity.     

Seismic anisotropy beneath Tibet: evidence for eastward extrusion of the Tibetan lithosphere?

Strong seismic anisotropy beneath Tibet has recently been reported from the study of SKS shear wave splitting. The fast split waves are generally polarized in an easterly direction, close to the present day direction of motion of the Tibetan crust relative to stable Eurasia, as deduced from Holocene slip rates on the major active faults in and around Tibet. This correlation may be taken to suggest that the whole Tibetan lithosphere is being extruded in front of indenting India and that the anisotropic layer is the deforming asthenosphere, that accommodates the motion of the Tibetan lithosphere relative to the fixed mantle at depth. Uncertainties about this motion are at present too large to bring unambiguous support to that view. Assuming that this view is correct however, a simple forward model is used to compute theoretical delay times as a function of the thickness of the anisotropic layer. The observed delay times would require a 50–100 km thick anisotropic layer beneath south-central Tibet and an over 200 km thick layer beneath north-central Tibet, where particularly hot asthenosphere has been inferred. This study suggests that the asthenospheric anisotropy due to present absolute block motion might be dominant under actively deforming continents.     

Simplification of strong ground motion considering inelastic responses of structures

Simplification of strong ground motions to 1 cycle sine waves was investigated from the elastic and inelastic earthquake response analyses and response analyses under sine wave input using single‐degree‐of‐freedom systems. Strong ground motions could be simplified to 1 cycle sine waves if large plastic deformations, with ductility factor more than 2, were assumed. This is because the approximate maximum responses from input sine waves are determined by the initial response cycle, due to period elongation and plastic energy dissipation of the systems. A sine wave whose acceleration amplitude is the peak ground acceleration (PGA) and whose period is that of an equivalent 1 cycle sine wave is proposed. The period of an equivalent sine wave is easily obtained from the elastic response acceleration spectrum of a seismic record. This means that the inelastic responses are approximately determined by the PGA and an equivalent 1 cycle sine wave period. Therefore, an equivalent 1 cycle sine wave period provides a single index to express the frequency characteristics of a strong ground motion. Copyright © 2000 John Wiley & Sons, Ltd.     

Double boundary layers in standing surface waves

Summary The double boundary-layer theory ofStuart (1963, 1966) andRiley (1965, 1967) is employed to investigate the mass transport velocity due to two-dimensional standing waves in a homogeneous fluid. The double boundary-layer structure occurs only at the bottom, and the existence of jet-like motions is predicted, due to the nature of the mean flows in the bottom boundary layers. A vector equation governing mean motions in the secondary layer is derived for threedimensional wave motions.     

多层介质中利用sPn与Pn到时差确定震源深度的研究

为利用sPn与Pn波的走时差测定震源深度,进一步提高地震震源深度的测定精度,推导多层介质下sPn与Pn波的走时差与震源深度的关系,发现走时差与震中距无关,只与震源深度及区域地壳速度模型有关。震源在同一层中,走时差曲线的斜率不变,而当震源位于不同层中时,sPn-Pn走时差曲线的斜率不同,并呈分段直线的走时差曲线形态。地壳速度结构纵向越不均匀,多层和单层介质下利用sPn-Pn走时差计算的震源深度误差就越大,走时曲线的各分段直线斜率相差越大;探讨地壳中sPn与Pn波传播路径相同但波速不同的单层地壳速度模型,发现单层介质下波速越大,测定的震源深度越大;对于同一地区相同的地壳分层结构,测得的震源深度随着泊松比的增大而减小。基于前人给出的中国五个典型块体多层平均地壳模型,给出sPn-Pn走时差与震源深度计算公式速查表。     

Some hydrodynamic problems for a nepheloid zone

Summary A layer of a few hundred meters thickness with suspended matter (a nepheloid zone) was discovered byEwing andThorndike [4]³⁾ near the bottom on the continental slope of the North Atlantic. A downward pressure gradient is produced in this layer due to increment of water density with suspensoid. When only the Coriolis force balances with this pressure gradient, a bottom nepheloid current flows southwestward parallel to the depth contours with a velocity of about 10 (cm/sec) for a slope of one degree. The pressure gradient for fluid with locally variable density above a sloping bottom is treated and an extra term due to density gradient along the slope is derived. The vertical profiles of the nepheloid current with an effect on the vertical eddy viscosity are computed. Two kinds of vertical distributions of eddy viscosity are determined from the observed nepheloid distributions and used in the calculations: constant but different values at two layers and those increasing with height. The effect of the change of density along the bottom is treated by introducing dimensionless variables. Rossby number of the nepheloid current becomes about 10^–2 indicating inertia terms to be negligible. Rossby number of turbidity currents ranges from 2 (in a decaying area) to 5 (developing area), suggesting that inertia terms are more important than Coriolis terms. The trajectories of turbidity currents are computed from motion of a mass of mud under the Coriolis force and friction, and the results are applied to those inferred byHand andEmery [6] in the San Diego Through off California.LGO Contribution Number 925.     

Effect of point source and heterogeneity on the propagation of SH-Waves in a viscoelastic layer over a viscoelastic half space

The present paper is concerned with the propagation of shear waves in a homogeneous viscoelastic isotropic layer lying over a semi-infinite heterogeneous viscoelastic isotropic half-space due to point source. The inhomogeneity parameters associated to rigidity, internal friction and density are assumed to be functions of depth. The dispersion equation of shear waves has been obtained using Green’s function technique. The dimensionless angular frequency has been plotted against dimensionless wave number for different values of inhomogeneity parameters. The effects of inhomogeneity have been shown in the dispersion curves. graphical user interface (GUI) software in MATLAB has been developed to show the effect of various inhomogeneity parameters on angular frequency. The topic can be of interest for geophysical applications in propagation of shear waves on the Earth’s crust.     

Propagation of shear waves in viscoelastic medium at irregular boundaries

The aim of the paper is to study the shear wave propagation in a viscoelastic layer over a semi-infinite viscoelastic half space due to irregularity in the viscoelastic layer. It is of great interest to study the propaga-tion of shear waves in the assumed medium having a non planar boundary due to its similarity to most of the real situations. The perturbation method is applied to find the displacement field. The effect of complex wave number on dissipation factor is analysed. Finally, as an application, the result obtained has been used to get the reflected field in viscoelastic layer when the shear wave is incident on an irregular boundary in the shape of parabolic irregularity as well as triangular notch. It is observed that the amplitude of this reflected wave decreases with increasing length of the notch, and increases with increasing depth of the irregularity.     

天然地震PS转换波动力学特征的初步研究

本文根据合成理论地震图的广义射线理论,并应用自行编写的FORTRAN计算程序,研究了单一界面、薄层和迭层界面上形成的PS转换波的动力学特征,讨论了沉积盖层及上部地壳介质吸收对转换波的影响。     

Source parameters for the January 1975 Brawley — Imperial valley earthquake swarm

The source parameters, moment, stress drop and source dimension are estimated for 61 events from the January 1975 Brawley earthquake swarm. Earthquakes studied range in local magnitude from 1.0 to 4.7. Stress drops range from 1 to 636 bars and increase with source depth. It is estimated that the sedimentary structure of the Imperial Valley amplifies shear waves by a factor of 2 to 3 in addition to the free surface amplification of 2. Estimates of moment from 10 sec surface waves are 4 to 6 times larger than the moment estimated from the relatively flat part of the local body wave spectrum at 1 sec. This may be due to after-slip on the fault, a long thin fault, or partial stress drop. It is shown that the experimentally determined ratio of stress drop to apparent stress should be approximately 4.0 when spectrum integration is used to obtainS-wave energy and theP-wave energy is 1/3 theS-wave energy.