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

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

横观各向同性层状半空间中凹陷地形对平面SH波的散射

沉积风化作用下土体表现出明显的成层特性和横观各向同性(TI)特性。将更为符合实际的TI介质模型引入到局部地形对地震波的散射问题研究有着重要的理论意义和工程实用价值。本文在TI土层平面外动力刚度矩阵的基础上,推导了斜线均布荷载动力格林函数,进而以该格林函数为基本解建立间接边界元方法研究了TI层状半空间中凹陷地形对平面SH波的散射问题。文中验证了方法的正确性,并分别在频域和时域内进行了数值计算分析。研究表明,TI介质中凹陷地形对SH波的散射与各向同性介质情况存在较大差异,凹陷附近地表位移幅值依赖于SH波的入射角度、入射频率和竖向与水平剪切模量比值;另外土体竖向与水平剪切模量比值对基岩上单一土层场地动力特性有着显著影响,进而使得TI层状场地凹陷附近位移幅值与各向同性情况差异变得更为复杂。     

斜入射平面SH波在层状半空间中沉积谷地周围的三维散射

采用直接刚度法计算自由场动力响应,以层状半空间中移动斜线均布荷载动力格林函数模拟散射波场,采用间接边界元方法求解了层状半空间中沉积谷地对斜入射平面SH波的三维散射问题.由于文中采用的层状场地三维动力刚度矩阵是精确的,且用于模拟散射波场的均布移动斜线荷载可以直接施加在沉积交界面处而不存在奇异性,所以本文方法具有很高精度.文中以均匀半空间和基岩上单一土层中沉积谷地对入射SH波的散射为例进行了数值计算,研究表明,沉积谷地对地震波的三维散射与二维散射之间存在本质差别;层状半空间中沉积谷地与均匀半空间中沉积谷地附近地表位移存在显著差异.     

模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre(GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.     

饱和层状场地中任意凹陷地形对入射SV波的散射

利用饱和土层的精确动力刚度矩阵和动力格林(Green)函数,采用间接边界元法,在频域内求解了层状饱和场地中任意凹陷地形对入射SV波的散射问题.通过自由场反应分析,求得凹陷地形表面各点的位移和各单元的应力响应;同样计算了虚拟分布荷载的格林影响函数,求得相应的位移和应力响应;根据边界条件确定虚拟分布荷载,将自由场位移响应和...     

流体饱和半空间二维地形三分量弹性波散射间接边界元模拟

无限长局部地形地震波斜入射响应问题称为二维三分量问题,在计算量远小于三维的情况下,一定程度上反映了近地表场地的三维动力响应特征.基于天然土体的成层性及固液两相耦合特性,以层状多孔介质内部移动线荷载(孔隙水压)动力格林函数作为基本解,开展流体饱和半空间二维地形三分量弹性波散射的2.5维间接边界元模拟研究.总场响应由自由波场和散射波场叠加构成,前者可由直接刚度法求得,后者则通过施加移动虚拟均布荷载和移动虚拟孔隙水压所产生的动力响应来模拟.该方法优势在于离散仅限于地形底边界(无须离散自由地表),格林函数计算不存在奇异性(荷载可直接加在边界上),容易控制计算精度,对复杂边界条件具有很强的适应性.在退化验证和精度比较的基础上,以梯形凹陷和半椭圆沉积地形为例,模拟了时域和频域的流体饱和半空间三维弹性波散射响应.研究表明:局部地形的地震动响应依赖于入射频率、入射角度、边界透水条件、土层刚度和土层厚度等,入射波、反射波和散射波相互干涉,极大延长了位移的振动持续时间.     

复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值...     

时间域起伏自由地表正演模拟综述

起伏地表对地震波有着非常大的影响,在起伏地表条件下,地震波的解析解不存在的,所以只能用数值方法来模拟地震波在地下的传播形式,起伏地形使得各种特殊的波现象出现,散射波、波的干扰、地滚波、以及体波和面波之间的转换,都给地震勘探带来了巨大的困难,因此,国内外学者对起伏自由地表的正演模拟做了大量的研究.本文从起伏地表的数值模拟方法、网格离散形式和自由边界实现方法三个方面系统的概括了时间域起伏自由地表正演模拟的研究进展.     

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

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

近表面周期裂纹阵对SH波散射问题的边界元解

针对半空间中一近表面周期裂纹对SH波的散射问题,本文采用边界元法进行研究。利用周期性的条件,将问题归结为一个典型条带中含单个裂纹的问题,并推导出了以裂纹张开位移为未知函数的边界积分方程。为了解决高阶奇异的边界积分方程的数值计算问题,本文采用解析的方法进行降阶。计算了半空间表面的位移,与其他方法求解结果的比较证明了数值计算的可靠性。     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography.The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

Evaluation of seismic ground motion induced by topographic irregularity

Results of an extensive numerical study on the 2D scattering of seismic waves by local topography are presented. The investigation has been conducted using the direct boundary element method. Several types of topography (slopes, canyons and ridges) are considered. The influences of some key parameters, such as exciting frequency and geometry of the irregular feature, on surface ground motion are studied in detail. It is found that local topographic conditions play an important role in the modification of seismic ground motion at the irregular feature itself and its neighbourhood. The present results can be considered to be useful from the viewpoint of earthquake engineering and seismology.     

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

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

求解二维弹性半空间散射问题的波函数组合法

提出了一种新的波函数组合法,将体波与面波特征函数组合求解不规则河谷地形在不同波型以不同角度入射的非一致位移输入问题.新方法严格满足不规则河谷散射表面自由条件,克服了目前波函数函数展开法不能完备表达二维不规则地形散射问题的缺点.数值算例表明波函数组合方法具有精确、稳定、计算量小的特点,可以应用于结构-地基动力相互作用的非一致输入问题.     

Review on localized boundary integral equation:Discrete wavenumber method for 2D irregular layers

The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Larner.Since that many numerical methods on the subject have been developed,such as pure numerical techniques,ray method and boundary method.The boundary method based on boundary integral equation is a semi-analytical method which is suitable to modeling wave field induced by irregular border.According to the property of the applied Green's function the boundary methods can be sorted into space domain boundary method and wavenumber domain boundary method.For both of them it is necessary to solve a large equation,which means much computation is needed.Thus,it is difficult for the boundary methods to be applied in simulating wave field with high frequency or in large range.To develop a new method with less computation is meaningful.For this purpose,localized boundary integral equation,i.e.,discrete wavenumber method is proposed.It is rooted in the Bouchon-Campillo method,an important wavenumber domain boundary method.Firstly the force on interface is separated into two parts:one is on flat part and the other on irregular part of the interface.Then Fourier transform is applied to identify their relation,the unknown distributes only on irregular part.Consequently computation efficiency is dramatically improved.Importantly its accuracy is the same as that of Bouchon-Campillo.     

二维频率测深边界单元法正演计算

在二维层状介质中含有横向非均匀体的地电断面情况下，对电偶源的偶极长度作了延长，并选取了适当的基本解及无穷远边界条件，因而将相应的含源谐变电磁场的边值问题转化为较简单的边界积分方程．用边界单元法求出了相应的电磁场分量，进而计算阻抗视电阻率，并对这一算法作了相应的检验．     

Reverse‐time migration from rugged topography using irregular,unstructured mesh

We developed a reverse‐time migration scheme that can image regions with rugged topography without requiring any approximations by adopting an irregular, unstructured‐grid modelling scheme. This grid, which can accurately describe surface topography and interfaces between high‐velocity‐contrast regions, is generated by Delaunay triangulation combined with the centroidal Voronoi tessellation method. The grid sizes vary according to the migration velocities, resulting in significant reduction of the number of discretized nodes compared with the number of nodes in the conventional regular‐grid scheme, particularly in the case wherein high near‐surface velocities exist. Moreover, the time sampling rate can be reduced substantially. The grid method, together with the irregular perfectly matched layer absorbing boundary condition, enables the proposed scheme to image regions of interest using curved artificial boundaries with fewer discretized nodes. We tested the proposed scheme using the 2D SEG Foothill synthetic dataset.     

Uniform half-plane elastodynamic problems by an approximate boundary element method

In this paper the applicability of an approximate Boundary Element Method to uniform half-plane elastodynamic problems is investigated. This method employs the concept of images to construct approximate fundamental solutions for the half-plane and does not require any half-plane surface discretization. The method is formulated in the frequency domain for the case of harmonic disturbances or the Laplace transform domain for the case of transient disturbances. In the latter case a numerical inversion of the transformed solution is necessary to obtain the time domain response. The proposed method can be used as an alternative to boundary element methods that either utilize the infinite plane fundamental solution and thus require a half-plane surface discretization, or employ the exact half-plane fundamental solution, which even though leads to no surface discretization, is of a very lengthy and complicated form. Two characteristics numerical examples are used to illustrate the proposed method and study its advantages and disadvantages.