三维层状孔隙介质中弹性波的一种积分表达式II：正确性检验和数值模拟实验

基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段.     

地震体波斜入射情形下台阶地形引起的波型转换

本文首次对地震体波斜入射情形下竖直、倾斜台阶地形引起的波型转换进行了数值模拟，结果表明：地震波的斜入射会使台阶上角点引起更强的转换面波；转换Rayleigh面波最大振幅可达弹性半空间表面自由场位移的1.1倍左右。     

各向异性半空间中浅埋孔洞对地表 反平面运动的影响

本文采用复变函数方法研究了稳态水平剪切波(SH波)在各向异性弹性半空间中任意形状孔洞上的散射及其对地面运动的影响.在变换空间中构造出自动满足各向异性半空间水平表面上应力为零的散射波函数,并根据孔洞表面应力为零的边界条件,采用最小二乘法求解散射波函数的系数.用介质的各向异性性质来模拟地质条件,给出了SH波作用下含圆孔、椭圆孔和方孔的各向异性半空间表面位移幅值的数值结果,并分析了介质特性、孔洞的形状、埋深、入射波波数及入射角度等因素对地面运动的影响规律.数值结果表明:介质的各向异性对含有孔洞的半空间表面的地表位移具有显著的影响;沿一定角度的入射波在某一频段内所引起的地表位移幅值比各向同性介质的可能要大,且随着孔洞埋深的增加,地表位移的幅值逐渐减小.      

求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟

间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.     

P-SV波斜入射时成层半空间自由场的时域算法

刘晶波和王艳提出了一种弹性水平成层半空间中平面P-SV波斜入射时平面内自由波场时域计算的一维化有限元方法,该方法采用黏性人工边界条件近似地模拟下部基岩半空间的辐射阻尼,可导致平面P-SV波以大角度入射时自由场计算精度降低.本文提出一种精确模拟基岩半空间辐射阻尼的人工边界条件.由于基岩半空间中外行波是传播方向已知的平面P波和SV波,利用弹性介质的应力-位移本构关系建立了人工边界处应力与速度的阻抗边界条件;采用该人工边界条件替代黏性边界条件,改进了P-SV波斜入射时成层半空间自由场的时域算法.数值试验表明,与采用黏性边界条件的自由场算法相比,改进方法具有更高的计算精度,其计算结果与理论解吻合更好.     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

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

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

地震波传播的哈密顿表述及辛几何算法

地震波传播过程本质上是能量在传播过程中逐步损耗直至殆尽的过程，而在实际应用中，常在无能量损耗假设下，用弹性波动方程或标量波动方程描述它.在哈密顿(Hamilton)体系表述下，地震波传播过程即为一个无限维的哈密顿系统随时间的演化过程.若不计能量损耗，波场演化过程实质上为一个单参数连续的辛变换，因而对应的数值算法应为辛几何算法.本文首先从地震波标量方程出发，给出哈密顿体系下地震波传播的表述，即任意两个时刻的波场是通过辛变换联系起来的.随后，把波场在时间和相空间离散化后，给出了用于波场计算的一些辛格式，如显式辛格式、隐式辛格式和蛙跳辛格式.并进一步讨论了有限差分格式和辛格式的异同.然后，应用显式辛格式和同阶的有限差分方法给出了同一理论速度模型下的波场和Marmousi速度模型下的单炮记录.数值结果表明，辛算法是一类可行的波场模拟的数值算法.在时间步长较小时，有限差分方法是辛算法的一个很好近似.文中的理论和方法，为地震波传播理论及实际应用研究提供了新的途径.     

多层弹性半空间中的地震波（一）

为了了解地震震源和地球介质的性质,很有必要对地震波的辐射、传播和衰减问题作仔细的分析。作为一种近似,可以暂且忽略地球的曲率,把传播地震波的地球介质视为多层半空间。为简便起见,地震波的衰减问题另作考虑。这样,便需要研究多层、均匀、各向同性和完全弹性半空间中地震震源辐射的地震波传播问题。 用哈斯克尔（Haskell）矩阵法解多层介质中弹性波的传播问题是很方便的。如果     

多层弹性半空间中的地震波（二）

一、引言 地震面波的频散性质、地震波辐射的方向性等特性已经广泛地用于地壳和上地幔结构以及震源机制的研究中,并且取得了许多有用的成果.研究地震波如何从震源辐射出来、如何在实际介质中传播和衰减的这一问题,对于利用地震波确定地壳和上地幔结构以及震源的参数,是很有必要的.关于这一问题的研究,已经作过许多工作。已往的工作中,为了分析方便,往往采用简单的地壳-上地幔模型或简单的震源模型,或两者都相     

Scattering and diffraction of plane P-waves in a 2-D elastic half-space Ⅱ:shallow arbitrary shaped canyon

Scattering and Diffraction of elastic in-plane P- and SV- waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P- and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P- and SV- scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.     

A NUMERICAL STUDY OF LAMB'S PROBLEM*

This paper considers propagation of elastodynamic waves in an imperfectly elastic half-space. Two different excitation modes are investigated: a buried source of compressional waves and a vertically directed areal load applied to the surface. Numerical integration of the analytical solution of the wave equation allows study of the vertical and horizontal components of displacement and/or particle velocity anywhere in the half-space. One case of particular interest concerns the examination of particle displacement and velocity at the surface in a circular area above the source. In another application seismograms generated by an explosive buried source are contrasted with seismograms generated by the transient application of a vertically directed load to the free surface. Still another application of considerable practical interest concerns the study of the nongeometrical pS—wave, in particular its characteristics as functions of range and depth. Finally, in the last application the behavior of a rarely observed wave (denoted here by the letter U) is studied in both elastic and visco-elastic half-spaces.     

弹性半空间位错内源的数值解

本文首先将一般剪切位错点源表示成四种基本点位错的组合,断层面倾角和滑动系数仅与组合系数有关.利用将三维动力学问题化为两维的处理方法,将每一基本点源问题化成求解位移函数的轴对称问题,最后将人工透射边界的处理方法同有限差分结合起来,提出一种简便、节省的计算理论地震图的方法.利用这种方法获得了四种基本点位错源的典型近场理论地震图,并初步讨论了近场地面运动的若干特征,表明浅源地震近场地面位移中震相不能分离,且点位错源在近场产生大小相当的水平与垂直位移.本文的方法可用于研究复杂的近场地震波传播问题.      

圆形隧洞对地下地震动的影响:SH波入射

地下结构对地震波的散射改变了场地的动力特性,无论是地上建筑还是地下结构的安全评价和抗震设计中,目前均没有很好的考虑由于地下结构的存在对原场地地震动尤其是地下地震动的影响。基于弹性波动理论,运用波函数展开法和镜像原理,分析了弹性半空间中圆形隧洞对柱面SH波的散射问题,得到了含圆形隧洞的弹性半空间位移解析解。通过数值算例分析了圆形隧洞对原场地地震动的影响,重点考察了隧洞埋深、隧洞半径和围岩衬砌模量比等参数的影响规律。结果表明,地下结构对沿线场地的动力特性有着显著的影响,对其自身以及沿线工程结构的抗震设计提供一定参考价值。     

MODAL ANALYSIS OF SOFT-SOIL SITES INCLUDING RADIATION DAMPING

For the one-dimensional analysis of soft-soil layers on an elastic half-space, a general form of analytical solution is developed for converting radiation damping due to energy leaking back to the half-space into equivalent modal damping, allowing the modal analysis technique to be extended to a site where radiation damping has to be accounted for. Closed-form solutions for equivalent modal damping ratios and effective modal participation factors are developed for a single layer with a shear wave velocity distribution varying from constant to linearly increasing with depth. Compact and recursive forms of solutions for equivalent modal damping ratios are developed for a system with an arbitrary number of homogeneous layers on an elastic half-space. Comparisons with numerical solutions show that the modal solutions are accurate. The nominal frequency of a site, i.e. the inverse of four times the total shear wave travel time through the layers, is an important parameter for estimating the high mode frequencies. A parameter study shows that for the same impedance ratio of the bottom layer to the elastic half-space, a system of soil layers with an increasing soil rigidity with depth has, in general, larger peak modal amplifications at the ground surface than does a single homogeneous layer on an elastic half-space, while a system with a decreasing soil rigidity with depth has smaller modal peak amplifications. © 1997 by John Wiley & Sons, Ltd.     

Two-dimensional scattering and diffraction of P- and SV-waves around a semi-circular canyon in an elastic half-space: An analytic solution via a stress-free wave function

A well-defined boundary-valued problem of wave scattering and diffraction in elastic half-space should have closed-form analytic solutions. This two-dimensional (2-D) scattering around a semi-circular canyon in elastic half-space subjected to seismic plane and cylindrical waves has long been a challenging boundary-value problem. In all cases, the diffracted waves will consist of both longitudinal (P-) and shear (S-) rotational waves. Together at the half-space surface, these in-plane longitudinal P- and shear SV-waves are not orthogonal over the infinite half-space flat-plane boundary. Thus, to simultaneously satisfy both the zero normal and shear stresses at the flat-plane boundary, some approximation of the geometry and/or wave functions often has to be made, or in some cases, relaxed (disregarded). This paper re-examines this two-dimensional (2-D) boundary-value problem from an applied mathematics points of view and redefines the proper form of the orthogonal cylindrical-wave functions for both the longitudinal P- and shear SV-waves so that they can together simultaneously satisfy the zero-stress boundary conditions at the half-space surface. With the zero-stress boundary conditions satisfied at the half-space surface, the most difficult part of the problem will be solved, and the remaining boundary conditions at the finite-canyon surface are then comparatively less complicated to solve. This is now a closed-form analytic solution of the 2-D boundary-valued problem satisfying the half-space zero-stress boundary conditions exactly.     

Weighted residual method for diffraction of plane P-waves in a 2-D elastic half-space Ⅲ: on an almost circular arbitrary-shaped alluvial valley

Scattering and diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong motion seismologists for over forty years. The case of out-ofplane SH-waves on the same elastic canyon that is semicircular in shape on the half-space surface is the first such problem that was solved by analytic closed-form solutions over forty years ago by Trifunac. The corresponding case of in-plane P-and SV-waves on the same circular canyon is a much more complicated problem because the in-plane P-and SV-scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by Lee and Liu. This paper uses their technique of defining these stress-free scattered waves, which Brandow and Lee previously used to solve the problem of the scattering and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape, to the study of the scattering and diffraction of these in-plane waves on an almost circular arbitrary-shaped alluvial valley.     

Weighted residual method for diffraction of plane P-waves in a 2-D elastic half-space III: on an almost circular arbitrary-shaped alluvial valley

Scattering and diffraction of elastic in-plane P- and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong motion seismologists for over forty years. The case of out-of-plane SH-waves on the same elastic canyon that is semicircular in shape on the half-space surface is the first such problem that was solved by analytic closed-form solutions over forty years ago by Trifunac. The corresponding case of in-plane P- and SV-waves on the same circular canyon is a much more complicated problem because the in-plane P- and SV- scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by Lee and Liu. This paper uses their technique of defining these stress-free scattered waves, which Brandow and Lee previously used to solve the problem of the scattering and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape, to the study of the scattering and diffraction of these in-plane waves on an almost circular arbitrary-shaped alluvial valley.