共查询到20条相似文献，搜索用时 561 毫秒
1.
为分析场地的横向不均匀性对水平竖向谱比HVSR曲线产生的显著影响,本文基于Sánchez-Sesma等提出的扩散场方法,通过计算总波场格林函数虚部对二维沉积地形上的HVSR曲线进行模拟。格林函数虚部则通过刚度矩阵和平面内斜线格林函数采用间接边界元方法进行求解。对二维沉积地形和相应的一维层状半空间的HVSR曲线进行了参数分析,着重讨论了沉积地形的形状、相对计算点位置等因素对HVSR曲线的影响规律。结果表明:沉积地形内外材料阻抗比对HVSR曲线的影响最为显著;随着沉积地形内外材料阻抗差异和沉积侧界面坡度的增大,沉积地形上HVSR曲线的第一峰值点的频率显著增大至相应层状半空间结果的3.3倍,同时HVSR曲线的形态呈现出平台现象;随着计算点到沉积边界距离的减小,HVSR曲线高频段幅值相对较大。根据本文得到的局部地形对HVSR曲线的影响规律,在进行场地勘探时可采用HVSR方法初步确定局部地形的分布位置以降低勘探成本。 相似文献
2.
The simulation of the seismic response of heterogeneous sedimentary basins under incident plane waves is computed using the Indirect Boundary Element Method (IBEM). To deal with these kinds of basins we used approximate analytical expressions for the two-dimensional Greens functions of a medium with constant-gradient wave propagation velocity. On the other hand, for the homogeneous half space underlying the sedimentary basin, the full space Greens functions were used. The response of semi-circular heterogeneous basins under incident SH waves is explored by means of the displacements in the frequency-space diagrams and synthetic seismograms. Moreover, we compared these results with those obtained for other homogeneous semi-circular models. The principal differences among them are pointed out. This simulation provided interesting results that displayed a complex amplification pattern in a rich spectrum of frequencies and locations. The maximum amplitudes levels were found around the edges of the heterogeneous sedimentary basins. In time domain some features characterize the seismic response of the basin which include enhanced generation and trapping of surface waves inside the sediments, and the reduced emission of seismic energy to the hard rock. In the heterogeneous models the lateral reflections of surface waves greatly influence the total displacements at the free surface in comparison with the homogeneous models where the displacements have a shared influence among both vertical and lateral reflections. 相似文献
3.
The dynamic response of infinitely long lined tunnels with a uniform cross-section buried into an elastic or viscoelastic half-space to body and surface harmonic seismic waves is numerically determined by a special direct boundary element method in the frequency domain. The waves have an arbitrary direction of propagation with respect to the axis of the tunnel and this renders the problem three-dimensional. However, this problem is effectively reduced to a two-dimensional one by a coordinate transformation and appropriate integration of the full-space dynamic fundamental solution along the direction of the tunnel axis. Quadratic isoparametric boundary line elements and advanced numerical integration techniques for the treatment of singular integrals produce results of high accuracy. Numerical results are presented for the case of an infinitely long lined tunnel of circular cross-section and compared against those of a full three-dimensional boundary element analysis, as well as those of other methods. Thus the proposed method is illustrated and its merits demonstrated. 相似文献
4.
《地震科学（英文版）》2019,(Z1)
A two and a half dimensional(2.5 D) multidomain indirect boundary element method(IBEM) is developed to study the wave scattering of obliquely incident P-, SVand SH-waves by a hill-valley staggered topography in a multi-layered half-space. The IBEM algorithm includes using2.5 D full-space and half-space Green's functions to construct scattered fields in decomposed closed and opened half-space regions, respectively, and using the dynamic stiffness method to solve the free fields. All regions are finally integrated by introducing the compatibility conditions to obtain the total wave fields. The proposed 2.5 D IBEM has the flexibility in dealing with complex boundaries by directly applying the fictitious loads on the regions' boundaries, with a less storage requirement compared to the full 3 D models. Besides, by combining the specific advantages of the two kinds of Green's functions, the method is well suitable for handling coupled topographies with high accuracy. The method is validated by comparison with published results for a single valley as well as a single hill topography. The effects of height-to-width ratio of hill and layering on dynamic responses are further parametrically investigated by numerical implementations in frequency domain. Results show that the interaction between valley and hills can lead to a more significant amplification within the valley region, and dynamic responses are deeply influenced by the height-towidth of hill, simultaneously depending on incident angle and frequency. Besides, the site effects become more complex when the stratification feature is taken into account. 相似文献
5.
A. Rodríguez-Castellanos F.J. Sánchez-Sesma C. Ortiz-Alemán M. Orozco-del-Castillo 《Soil Dynamics and Earthquake Engineering》2011
In earthquake engineering and seismology it is of interest to know the surface motion at a given site due to the incoming and scattered seismic waves by surface geology. This can be formulated in terms of diffraction of elastic waves and then the indirect boundary element method (IBEM) for dynamic elasticity is used. It is based on the explicit construction of diffracted waves at the boundaries from which they radiate. This provides the analyst with insight on the physics of diffraction. The IBEM has been applied to study the amplification of elastic waves in irregular soil profiles. From the strong or weak satisfaction of boundary conditions and a simple analytical discretization scheme a linear system of equations for the boundary sources is obtained. Here, we explore the use of a weak discretization strategy with more collocation points than force densities. The least squares enforcement of boundary conditions leads to a system with reduced number of unknowns. This approach naturally allows one to use both coarser and finer boundary discretizations for smooth and rapidly varying profiles, respectively. A well studied semicircular canyon under incident P or SV in-plane waves is used to calibrate this method. Several benefits are obtained using mixed meshing that leads to the least squares condensation of the IBEM. 相似文献
6.
《地震科学（英文版）》2019,(Z1)
The indirect boundary element method(IBEM) is applied to investigate the scattering of elastic waves around a 3-D sedimentary basin filled with fluidsaturated poroelastic medium. Based on this method, the free field and scattered field can be solved according to the boundary conditions. And the numerical accuracy has been verified. The effects of parameters on elastic wave scattering are studied, such as boundary condition, incident frequency,incident angle and porosity of medium. Numerical results illustrate that the amplification effect of surface displacement near poroelastic sedimentary basin is notable. In addition, for the case of large porosity the drainage condition has a significant impact on the response amplitude. Due to the fluid exchange at the interface under the drained condition, the displacement amplitude can be much larger than that under the undrained condition in present study. The study can provide a theoretical basis for the anti-seismic design of engineering structures located in sedimentary basin. 相似文献
7.
Elastic wave propagation in an irregularly layered medium 总被引：1，自引：0，他引：1
Rossana Vai Jos Manuel Castillo-Covarrubias Francisco J. Snchez-Sesma Dimitri Komatitsch Jean-Pierre Vilotte 《Soil Dynamics and Earthquake Engineering》1999,18(1):385
The indirect boundary element method (IBEM) is used to simulate wave propagation in two-dimensional irregularly layered elastic media for internal line sources. The method is based on the integral representation for scattered elastic waves using single layer boundary sources. Fulfillment of the boundary conditions leads to a system of integral equations. Results are obtained in the frequency domain and seismograins are computed through Fourier synthesis. In order to test and validate the method we present various comparisons between our results and the time series obtained analytically for a buried line source in a half-space and by using the recently developed spectral element method (SEM). 相似文献
8.
Michel Bouchon 《Pure and Applied Geophysics》1996,148(1-2):3-20
We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering. 相似文献
9.
Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ) 
下载免费PDF全文
下载免费PDF全文 Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and El Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered. 相似文献
10.
基于一种高精度间接边界元法(IBEM), 实现了沉积盆地三维地震响应的频域、 时域精细求解, 并以半空间中椭球形沉积盆地对平面P波和SV波的散射为例, 着重探讨了入射角度、 入射波型、 入射频率、 盆地长宽比和深宽比对沉积盆地地震动放大效应的影响规律． 结果表明: 盆地形状对地震波的放大效应和空间分布状态具有显著影响, 且具体规律受控于入射波频段． ① 随着盆地深度增大, 盆地边缘面波发育更为充分, 在较宽频段内均会出现显著的地震动放大效应, 且深盆地的放大区域集中于盆地中部． ② 圆形盆地对地震波的汇聚效应最为显著, 而狭长盆地对地震波的汇聚作用相对较弱, 高频情况下可在盆地内部形成多个聚焦区域． ③ 不同波型入射下, 盆地对地震动放大效应的机制有所差异: P波入射下, 竖向位移放大主要是由于盆地边缘面波由四周向中部汇聚所致; SV波入射下, 边缘面波汇聚效应相对较弱, 而当盆地较深时, 底部透射体波和边缘面波易形成同相干涉从而显著放大地震动． 按盆地内外介质波速比为1/2, P波和SV波垂直入射下频域最大放大倍数分别为25和15, 时域放大倍数约为4.0和3.7(雷克子波)． ④ 低频波入射下, 位移从盆地中部向边缘逐渐减小, 且浅层沉积盆地对地表位移幅值的放大作用不明显． ⑤ P波和SV波的入射角度对盆地地震动放大幅值及空间分布特征也具有显著影响． 相似文献
11.
针对三维沉积盆地对球面波的散射问题,发展一种快速宽频间接边界元方法(IBEM)。利用ANSYS建立求解模型,基于Intel-Fortran编译器编译相应的计算程序,对基岩半空间三维半椭球形盆地对球面波的散射进行了数值分析,着重探讨入射波频率、波源埋深、波源与不规则地形(沉积盆地)距离等参数对地震动特性的影响规律,计算方法适用于求解任意复杂形状的局部场地,实现对沉积盆地三维地震响应的频域精细求解。研究结果表明:随着入射波频率增大,沉积盆地相对半空间的放大效应越显著,且干涉效应越强烈;高频波入射时,由于盆地底部透射的体波和盆地内部由体波转换为的面波以及从盆地边缘处透射的体波相互叠加,次方向位移出现一定的边缘效应。位移频谱分析表明,沉积盆地在较低频率时,地表不同点位的主方向位移幅值基本相同;在较高频率域内,地表不同点位的位移频谱特性差别很大;球面波入射与平面波入射相比之下,位移幅值均有所降低。整体来看,随着膨胀波源与盆地水平距离增大,地表主方向最大位移幅值降低,地表次方向位移聚焦区域增加且分布更分散。 相似文献
12.
针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明: 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响. 相似文献
13.
The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories. 相似文献
14.
Francisco Luzón Francisco J. Sánchez-Sesma J. Alfonso Pérez-Ruiz Leonardo Ramírez-Guzmán Andrés Pech 《Soil Dynamics and Earthquake Engineering》2009
The propagation of seismic P and SV waves within inhomogeneous alluvial valleys has been investigated using the indirect boundary element method (IBEM). An improvement on the formulae published early, for the 2D Green's functions in an inhomogeneous medium is presented in this work. A modification has been done over these functions in view of its connection to the ray theory. An accuracy analysis validates these modified Green's functions computing the relative error in frequency domain for the SH case, and with a quantitative analysis obtaining envelope and phase misfits of the solution in time domain, for the P–SV case. 相似文献
15.
Seismic response of multi-layered basins with velocity gradients upon incidence of plane shear waves
A boundary integral scheme based on boundary-integral discrete-wave-number approach has been developed to compute the seismic response of two-dimensional irregular-shaped basins with horizontal soil layers. Each layer exhibits a linear gradient of shear wave velocity with depth. The approach combines the boundary-integral representation of the seismic wave field outside the basin with plane wave representation of the seismic wave field inside the basin. The propagation throughout the layers is performed by matrix propagators in which the effect of the vertical variation of the velocity is incorporated by using confluent hyper-geometric functions of the Whittaker type. Our method is tested against otherwell-accepted solutions for the case of a circular basin with excellent agreement. Test of the ground response for a semi-circular basin with radius a shows that stable solutions are obtained if the chosen model parameters satisfy following conditions: (1) the distance from the sources to the interface is greater than 0·1a; (2) the distance between the sources is smaller than a quarter of the incident wavelength; and (3) the discrete wave-number step is smaller than 2π/4a. The computation of ground response of basins with a sharp interface and several horizontal deposits leads to the following main results: (1) the amplification of a basin with velocity gradients is larger than that of a basin with homogeneous layers; (2) the frequencies of the second- and third-order harmonics for a basin with velocity gradients are lower than those of a basin with homogeneous layers; and (3) the response amplitude of the basin with velocity gradients attenuates more slowly in time domain than when layers are homogeneous. Since these results have been obtained for realistic values of basin geometrical and mechanical consideration, they should find some interest in earthquake engineering or seismic microzonation studies. © 1998 John Wiley & Sons, Ltd. 相似文献
16.
2008年11月10日在青海柴达木盆地北缘发生了大柴旦M_W6.3地震,为了研究该地震的区域地震波传播与地面运动特征,本文利用地质资料和地壳速度结构研究成果,构建了柴达木盆地及周边区域三维传播介质模型,采用有限差分方法模拟了大柴旦地震波场传播过程以及地面运动分布特征.结果表明,柴达木盆地对波场传播有明显影响,表现为地震波传入盆地后在边界产生次生面波,盆地沉积物对地震波具有围陷作用,地震地面运动在盆地内振幅增大、持时延长.模拟结果给出的地震地面运动峰值速度分布以及理论地震图均和观测结果符合较好,反映数值模拟较好地给出了观测地面运动的主要特征以及传播介质模型的合理性. 相似文献
17.
实际工程场地中的断层通常具有三维尺寸,应进行三维分析。本文考虑均匀弹性半空间内存在不同长度、深度的三维断层,由下向上垂直入射一个脉冲波,主要采用三维时域显式有限元结合透射人工边界的方法及相应编写的三维显示有限元波动程序DSI3,计算断层附近地震地面运动的时程反应,研究含有三维断层场地地面运动的一些特点,并与自由场反应对比,以说明三维竖向断层对地震波传播的特点和对断层附近地面运动的影响。 相似文献
18.
在层状半空间精确动力刚度矩阵和斜线荷载动力格林函数的基础上建立间接边界元方法,在频域内求解无破碎带断层场地对入射平面SH波的散射。为方便求解,将总波场分解为自由波场和散射波场,自由波场由直接刚度法求得,断层两侧的散射波场通过在断层面上分别对两侧施加均布斜线荷载产生的动力响应来模拟,虚拟荷载的密度可通过引入断层表面的边界条件确定,最后叠加自由波场和散射波场求得总波场。以有落差断层和无落差断层模型为例进行数值计算,分析断层落差、断层倾角以及断层两侧介质的刚度比对散射效应的影响。研究表明,断层落差与波长相当时,断层对SH波的放大作用最大;地表位移幅值随着断层倾角的增大逐渐增大;若断层无落差且其两侧刚度不同时,一般刚度较小一侧地表位移幅值较大且振荡更为剧烈,波从刚度较小一侧入射时位移幅值放大尤为显著。 相似文献
19.
Scattering of incident plane harmonic pseudo P‐, SH‐, and SV‐waves by a two‐dimensional basin of arbitrary shape is investigated by using an indirect boundary integral equation approach. The basin and surrounding half‐space are assumed to be generally anisotropic, homogeneous, linearly elastic solids. No material symmetries are assumed. The unknown scattered waves are expressed as linear combinations of full‐space time‐harmonic two‐dimensional Green functions. Using the Radon transform, the Green functions are obtained in the form of finite integrals over a unit circle. An algorithm for the accurate and efficient numerical evaluation of the Green functions is discussed. A detailed convergence and parametric analysis of the problem is presented. Excellent agreement is obtained with isotropic results available in the literature. Steady‐state surface ground motion is presented for semi‐circular basins with generally anisotropic material properties. The results show that surface motion strongly depends upon the material properties of the basin as well as the angle of incidence and frequency of the incident wave. Significant mode conversion can be observed for general triclinic materials which are not present in isotropic models. Comparison with an isotropic basin response demonstrates that anisotropy is very important for assessing the nature of surface motion atop basins. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
20.
结合单相介质动力格林函数和流体域格林函数,将间接边界元方法拓展到含流体层河谷对地震波散射的求解,并结合具体算例进行大量参数分析。研究结果表明,含流体层河谷地形对平面P波、SV波入射时的地震响应受控于入射波频率、入射波角度及流体深度等多种因素。总体来看:① 在低频域内,含流体河谷底部及附近地表的频谱特性与不含流体的河谷反应基本一致;② P波入射时在水层体系共振频率处,河谷底部位移缩小效应显著,而此频率处流体表面位移达到最大;③ 流体层具有吸收地震波能量的作用,流体深度越大,河谷表面及附近地表的地震动位移越小。研究成果可在一定程度上为河谷地形附近地震动效应的评估及防震减灾工作提供理论依据。 相似文献