弹性层状半空间中凸起地形对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决凸起地形对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中凸起地形对入射平面SH波的放大作用问题。本文模型的显著特点之一是考虑了层状半空间的动力特性以及层状半空间和凸起地形的阻尼;特点之二是计算精度高。文中以基岩上单一土层中半圆凸起地形对入射平面SH波的放大作用为例进行了数值计算分析。研究表明,基岩上单一土层中凸起地形对入射平面SH波放大作用和均匀半空间中凸起地形有着本质的差别;土层动力特性不仅影响凸起地形地表位移的幅值,还会影响地表位移的频谱;阻尼会显著降低凸起地形对高频波的放大作用。     

弹性层状半空间中沉积谷地对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决沉积谷地对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中沉积谷地对入射平面SH波的放大作用问题。通过自由场反应分析,求得假想边界上各点的位移和各单元的应力反应。施加虚拟均布荷载,求得位移和应力的格林函数。根据应力和位移连续边界条件确定虚拟分布荷载,将自由场反应和虚拟分布荷载产生的反应叠加起来,即得到问题的解答。最后,以基岩上单一土层中沉积谷地对入射平面SH波的放大作用为例进行了数值计算分析,结果表明,层状半空间中沉积谷地和均匀半空间中沉积谷地对波的放大作用存在显著差别。文中分析了造成差别的原因,并讨论了覆盖层厚度和刚度对放大作用的影响。     

瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明: 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响.      

无破碎带断层场地对SH波的散射

在层状半空间精确动力刚度矩阵和斜线荷载动力格林函数的基础上建立间接边界元方法,在频域内求解无破碎带断层场地对入射平面SH波的散射。为方便求解,将总波场分解为自由波场和散射波场,自由波场由直接刚度法求得,断层两侧的散射波场通过在断层面上分别对两侧施加均布斜线荷载产生的动力响应来模拟,虚拟荷载的密度可通过引入断层表面的边界条件确定,最后叠加自由波场和散射波场求得总波场。以有落差断层和无落差断层模型为例进行数值计算,分析断层落差、断层倾角以及断层两侧介质的刚度比对散射效应的影响。研究表明,断层落差与波长相当时,断层对SH波的放大作用最大;地表位移幅值随着断层倾角的增大逐渐增大;若断层无落差且其两侧刚度不同时,一般刚度较小一侧地表位移幅值较大且振荡更为剧烈,波从刚度较小一侧入射时位移幅值放大尤为显著。     

Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

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.     

3D scattering of obliquely incident plane SV waves by an alluvial valley embedded in a fluid-saturated,poroelastic layered half-space

The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obliquely incident SV waves. A wave-number transform is first applied along the valley’s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is confirmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are discussed in detail resulting in some conclusions.     

平面SV波在层状半空间中沉积谷地周围的散射

采用间接边界元方法求解了入射平面SV波在层状半空间中沉积谷地周围的散射问题.问题的解答包含自由场和散射场两部分.自由场可由直接刚度法求得,散射场由层状半空间中斜线荷载动力格林函数来模拟.文中以入射平面SV波在基岩上单一土层中沉积谷地周围散射为例研究了土层和沉积谷地周围的影响.结果表明,由于考虑了土层的动力特性,平面SV...     

平面SV波在饱和半空间中沉积谷地周围的散射

采用一种特殊的间接边界积分方程法,求解了平面SV波在饱和半空间中任意形状沉积谷地周围的二维散射问题。结合饱和半空间中膨胀波源和剪切波源格林函数,由分布在沉积和半空间交界面附近两虚拟波源面上的波源分别构造沉积内外的散射波场,由交界面连续条件建立方程并求解确定虚拟波源密度,总波场反应即可由自由波场和散射波场叠加而得。然后通过边界条件验算、退化解答与现有结果的比较以及稳定性检验,验证了方法的计算精度。通过一组典型算例,研究了平面SV波在饱和半空间中沉积谷地周围散射的基本规律,详细给出了不同参数情况沉积谷地附近地表位移幅值和孔隙水压,着重分析了入射SV波频率和角度、边界渗透条件、沉积孔隙率等因素对场地反应的影响,得出了一些有益的结论。     

2.5D scattering of incident plane SV waves by a canyon in layered half-space

This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green’s functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fic...     

Scattering of plane sh waves by a semi-cylindrical hill

A new analytical method is presented to study the problem of scattering of plane SH waves by a semi-cylindrical hill in an otherwise homogeneous, isotropic and elastic two-dimensional half-space, using the series of wave functions and a new expansion technique. The results obtained show that (1) the hill has quite considerable effects on ground motions for both the points on the hill and its environs, (2) these effects depend mainly on the frequency, the angle of wave incidence and the ratio of radius of hill to one-half the wave length of incident waves, and (3) prominent and depression topographies having equal form and radius for the same incident waves would cause quite different mechanisms of wave propagation in both response performance and magnitude. The results, doubtless, would be useful for further deepening knowledge of the effects of prominent topography on seismic ground motion and for testing the accuracy of various approximate methods so far available.     

不规则场地上基础的动阻抗及其在入射波作用下的输入运动

本文根据边界元方法建立了位不规则场上刚体的动阻抗和在入射平面波作用下的有效输入运动的分析模型，分析模型考虑了不规则场地和基础对入射波的散射作用以及土与基础的相互作用，通过验证确认了本方法的正确性，文中计算了凹陷，高地和盆地三种不规则场地土不同条件基础的动阻和有效输入的运动，并与半空间地基上相应基础的情况作了对比，计算表明，当基础尺寸与不规则场地范围可比时有必要用本文模型分析不规则场地的影响和土一结     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

二维平面内考虑横向不均匀体散射作用的水平竖向谱比模拟与应用

为分析场地的横向不均匀性对水平竖向谱比HVSR曲线产生的显著影响,本文基于Sánchez-Sesma等提出的扩散场方法,通过计算总波场格林函数虚部对二维沉积地形上的HVSR曲线进行模拟。格林函数虚部则通过刚度矩阵和平面内斜线格林函数采用间接边界元方法进行求解。对二维沉积地形和相应的一维层状半空间的HVSR曲线进行了参数分析,着重讨论了沉积地形的形状、相对计算点位置等因素对HVSR曲线的影响规律。结果表明:沉积地形内外材料阻抗比对HVSR曲线的影响最为显著;随着沉积地形内外材料阻抗差异和沉积侧界面坡度的增大,沉积地形上HVSR曲线的第一峰值点的频率显著增大至相应层状半空间结果的3.3倍,同时HVSR曲线的形态呈现出平台现象;随着计算点到沉积边界距离的减小,HVSR曲线高频段幅值相对较大。根据本文得到的局部地形对HVSR曲线的影响规律,在进行场地勘探时可采用HVSR方法初步确定局部地形的分布位置以降低勘探成本。      

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

The dynamic soil-tunnel interaction is studied by indirect boundary element method (IBEM), using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

3D magnetotelluric modelling including surface topography

An edge finite‐element method has been applied to compute magnetotelluric (MT) responses to three‐dimensional (3D) earth topography. The finite‐element algorithm uses a single edge shape function at each edge of hexahedral elements, guaranteeing the continuity of the tangential electric field while conserving the continuity of magnetic flux at boundaries. We solve the resulting system of equations using the biconjugate gradient method with a Jacobian preconditioner. The solution gives electric fields parallel to the slope of a surface relief that is often encountered in MT surveys. The algorithm is successfully verified by comparison with other numerical solutions for a 3D‐2 model for comparison of modelling methods for EM induction and a ridge model. We use a 3D trapezoidal‐hill model to investigate 3D topographic effects, which are caused mainly by galvanic effects, not only in the Zxy mode but also in the Zyx mode. If a 3D topography were approximated by a two‐dimensional topography therefore errors occurring in the transverse electric mode would be more serious than those in the transverse magnetic mode.     

Surface motion of a semi-cylindrical hill for incident plane SV waves： Analytical solution

Introduction Scattering of incident waves by local topographies is one of the most attractive topics in the field of engineering seismology. They can be resolved by either a numerical method or an analyti-cal method. The numerical methods include finite difference method, finite element method, boundary element method, etc. The analytical method is the wave function expansion method. Al-though numerical methods can be used for arbitrary-shaped topographies, the analytical solutions are still …     

局部地形下入射波散射效应对场地地震响应的影响

局部地形下入射波散射效应对场地地震响应的影响是目前众多学者所关注的问题.其影响主要由入射波频率、入射角度、地形几何形状、介质性质等几个因素所制约,所反应出来的散射效应主要体现在地表位移变化上.本文归纳了凹陷地形、凸起地形、沉积谷地以及复合地形这四类局部地形入射波散射效应对场地地震响应的影响,对目前研究成果进行了评述,并针对数学技术、工程应用、模型建立等7个方面指出了现今存在的问题和今后发展的方向.     

Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method

In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion.     

Propagation and amplification of seismic waves in 2-D valleys excited by obliquely incident P- and SV-waves

The seismic motion in sediment-filled valleys due to incident SH-waves has been studied exhaustively. However, the response of such geologic structures to incident SV- and P-waves has not been studied as thoroughly. The response of a 2-D model of the valley of Caracas, Venezuela—a NS cross-section through the Palos Grandes district—to incident plane SV- and P-waves is investigated using the discrete wave number boundary element method. It is observed that the differences in the predictions of the 1-D and 2-D models are more pronounced for SV-waves than for SH-waves, especially when SV-waves are incident at (or near) the critical angle i_c. The valley responds very strongly to the horizontally propagating P-wave (SP-wave) which is induced when SV-waves, incident at the critical angle, interact with the free surface of the half-space. However, the SP-wave, being a wave diffracted at a boundary, is likely to be sensitive to impedance contrasts, to the presence of other interfaces in the medium, and to the topography surrounding the valley. These aspects of the problem need further investigation.