散射问题中柱面坐标波函数的全域变换公式

为满足偏心源和多体散射理论研究发展的需求,阐明了柱面坐标下波函数变换的基本类型,指出了现有变换公式的局限性,提出了全域变换公式,给出了应用实例并讨论了应用前景。研究表明,不同柱面坐标下波函数变换可以分为内域和外域两种基本类型,现有Graf加法公式仅适用于内域变换问题。提出的新公式突破了现有公式的限制,可完成弹性波散射问题中不同圆柱坐标系下波函数全空间域上的变换。应用新公式,构造了圆形空穴的外线源散射和内线源散射的解答,通过实例与大圆弧近似方法进行了对比,表明了新公式的优势和有效性。波函数变换基本类型的划分、全域变换概念的引入和公式构建的原理具有普适性,可以推广到其它坐标系下波函数的坐标变换公式的建立中。     

圆弧状多层沉积谷地在平面P波入射下稳态响应的解析解

规则的沉积谷地对地震波散射问题解析解的研究是地球物理学、地震学与地震工程学等学科领域重要的研究课题之一.Trifunac最先利用波函数的Fourier-Bessel级数展开在极坐标系下给出了半圆形沉积谷地对平面SH波散射问题的精确解析解.其后,Wong和Trifunac利用波函数的Mathieu级数展开在椭圆坐标系下给出了半椭圆形沉积谷地对平面SH波散射问题的解析解.     

任意圆弧形凸起地形中隧洞对入射平面SH波的影响

借助辅助函数的思想， 采用波函数展开法给出了圆弧形凸起地形中隧洞对入射平面SH波影响问题的一个解析解．数值结果表明，凸起地形中隧洞的存在以及隧洞大小，对凸起地形表面运动，以及凸起地形中的隧洞的动应力集中均具有显著影响．      

层状各向异性地层中含环状天线槽随钻方位电磁波测井几何因子算法

本文基于感应测井中的几何因子与一阶Born近似理论,通过柱坐标系下耦合势Helmholtz方程三维有限体积法研究建立了各向异性地层中含环状天线槽的随钻方位电磁波测井几何因子.首先,通过引入环状电流源与径向磁偶极子源电磁场Green函数,并结合Green第二积分公式,推导出随钻方位电磁波测井仪器中环形电场和横向磁场微小变化与各向异性地层中水平和垂直电导率相对摄动量之间的关系,得到轴向和横向感应电动势微小变化量计算方法与柱坐标系中轴向和横向分量空间灵敏度函数(三维几何因子)表达式.在此基础上,进一步给出径向和纵向微分几何因子计算公式.然后,应用耦合势三维有限体积法确定发射线圈产生的感应电磁场以及环状电流源与径向磁偶极子源电磁场Green函数的数值解.最后,给出不同频率、不同倾角、不同各向异性系数等多种情况下该仪器的空间灵敏度函数与微分几何因子的数值结果,用于分析考查仪器的响应特征与空间探测能力.     

分层坐标变换法起伏自由地表弹性波叠前逆时偏移

传统有限差分方法在处理起伏地表时存在一些困难,而坐标变换法可将起伏地表映射为水平地表以克服此缺点.但同时,地下构造被变换得更加复杂,导致了波传播和成像的不准确.本文提出了一种分层的坐标变换方法,并应用到了弹性波逆时偏移中,此方法既可以克服起伏地表的影响,又可以不破坏地下构造.波场正向延拓、逆时延拓和分离是在辅助坐标系下完成的,而成像是在笛卡尔坐标系下完成的.通过对简单起伏模型和中原起伏模型的试算证明了本文提出方法的准确性.同时,对两种极端起伏地层高程不准确的情况进行测试可以看出:分层坐标变换逆时偏移方法的成像效果远好于传统坐标变换方法.     

圆弧状凹陷地形对平面SV波散射问题的宽频带解析解

利用波函数的Fourier-Bessel级数展开法,推导了具有不同深宽比的圆弧状凹陷地形对入射平面SV波二维散射问题的解析解。区别于现有其他解析解,利用柱函数的渐近性质,使得散射波的待定系数得以直接确定,避免了线性方程组的求解以及相应的高频波入射下的数值计算问题,从而拓展了解析解适用的频带范围。通过与已有解析解的比较论证了该解析解的正确性,并在一个较宽的频带范围内研究了圆弧状凹陷地形对入射平面SV波的散射效应。     

分离式双圆形隧道衬砌对平面SH波的散射

采用波函数展开法,将入射、散射和折射SH波的位移势函数展开成Fourier-Bessel函数的级数形式,根据Graf加法定理,得到了任一局部圆柱坐标系下的波场势函数的表达式.根据隧道衬砌与周围岩体都在界面处应力和位移连续及衬砌内边界完全自由的边界条件,得到了分离式双隧道衬砌对入射平面,SH波的多重散射问题的理论解.通过数值计算分析得出:入射频率、衬砌间距、衬砌与围岩的模量比及衬砌的厚度等是影响衬砌内侧的无量纲位移和动应力集中因子分布的重要因素,为双衬砌或多衬砌的地震评价提供了理论依据.     

三维地质体磁场梯度计算理论与方法

系统地推导建立了均匀磁化三维多面体磁场梯度理论公式及计算方法,解决了任意三度体磁场梯度正演问题,得出了在坐标系旋转情况下,磁场梯度的转换与磁场的转换具有本质的区别,并且在二维坐标系下,具有统一的任意阶梯度转换公式,在三维坐标系下,没有统一表达式的重要结论.模型计算表明.所建立的磁场梯度计算理论及方法是正确的.     

利用伪谱法模拟横向各向同性介质中的波

介质的弹性常数为三维四阶张量的分量,共有８１个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为２１个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有５个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的（即对称轴可以不平等于铅直方向）,在此情况下,需要进行坐标变换,如果已知介质在某一坐标系（其坐标轴平行或垂直于介质的对称轴）中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。     

横向各向同性地层斜井中正交偶极子激发声场的数值模拟

本文采用三维应力-速度有限差分(SV-FD)方法,数值模拟了横向各向同性(TI)地层对称主轴与井轴斜交情况下正交偶极子声源激发的井孔声场.主要解决了与倾斜角有关的三维空间弹性模量矩阵的推导,柱坐标系下应力-速度有限差分方程组的建立,井轴上场点奇异性与内边界处理等几个关键问题,提高了计算精度.在横向各向同性地层对称主轴与井轴平行的情况下,与实轴积分法所得结果进行了对比,验证了本文方法的正确性.计算了不同倾角情况下xx和yy两分量的弯曲波,并用频域加权相似法提取了弯曲波频散曲线,结果显示了横向各向同性介质中不同方向偏振的弯曲波传播的分裂现象,其频散曲线在低频段分裂,随频率增大而逐渐重合.弯曲波低频截止频率处速度与理论公式得到的横波速度基本符合.     

Scattering of SH waves induced by a symmetrical V-shaped canyon: a unified analytical solution

This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.     

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.     

球坐标系下谱元法三维地震波场模拟

谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系下求解三维地震波方程的谱元法.另一方面,计算Fréchet敏感核是进行全波形反演的关键,本文借助伴随原理,推导了全波走时层析成像三维Fréchet敏感核表达式.为了验证球坐标系下谱元法的精度,我们将数值模拟结果与normal mode方法得到的解析解在1-D PREM模型下进行了对比.同时,我们将此方法应用到华北克拉通区域,以期获得地球内部结构精确成像.基于3-D全球径向各向异性地幔模型S362ANI和3-D地壳模型Crust1.0,我们建立了华北克拉通初始3-D背景模型,并将数值模拟结果与实际观测台站记录波形资料进行对比分析,利用互相关方法提取走时残差,最后给出了Fréchet敏感核在3-D空间中的分布,这些工作为下一步开展球坐标系下三维大尺度全波形反演奠定了基础.     

二维地表相邻多个半圆弧沟谷对SH波的散射

本文由解析法研究了二维地表相信多个圆弧沟谷对稳态ＳＨ波的散射问题。文中把相邻多个圆弧沟谷对波的多重散射表达为各局部坐标的Ｆｏｕｒｉｅｒ－Ｂｅｓｓｅｌ级数之和，再利用Ｇｒａｆ加法公式在某一局部坐标中表达为双重级数形式，将问题归结为对一组线性方程的求解。最后，作为例子，给出了相信两圆弧沟谷对波散射的数值结果，讨论了沟谷对波的屏蔽、隔振作用及沟谷之间的相互作用。     

3-D shell analysis of cylindrical underground structures under seismic shear (S) wave action

The 3-D shell theory is employed in order to provide a new perspective to earthquake-induced strains in long cylindrical underground structures, when soil-structure interaction can be ignored. In this way, it is possible to derive analytical expressions for the distribution along the cross-section of axial, hoop and shear strains and also proceed to their consistent superposition in order to obtain the corresponding principal and von Mises strains. The resulting analytical solutions are verified against the results of 3-D dynamic FEM analyses. Seismic design strains are consequently established after optimization of the analytical solutions against the random angles which define the direction of wave propagation relative to the longitudinal structure axis, the direction of particle motion and the location on the structure cross-section. The basic approach is demonstrated herein for harmonic shear (S) waves with plane front, propagating in a homogeneous half-space or in a two layer profile, where soft soil overlays the bedrock.     

Broadband Analytical Solution to the Scattering of Plane SV Waves by Circular Cylindrical Canyons

Based on Fourier-Bessel series expansion of wave functions,an analytical solution to 2-D scattering ofincident plane SV waves by circular cylindrical canyons with variable depthto-width ratios is deduced in this paper. Unlike other analytical solutions,this paper uses the asymptotic behavior of the cylindrical function to directly define the undetermined coefficients of scattered waves,thus,avoiding solving linear equation systems and the related numerical computation problems under high-frequency incident waves,thereby broadening the applicable frequency range of analytical solutions. Through comparison with existing analytical solutions,the correctness of this solution is demonstrated. Finally, the incident plane SV wave scattering effect under circular cylindrical canyons in wider frequency bands is explored.     

Application of a second-order paraxial boundary condition to problems of dynamics of circular foundations on a porous layered half-space

A half-space finite element and a consistent transmitting boundary in a cylindrical coordinate system are developed for analysis of rigid circular (or cylindrical) foundations in a water-saturated porous layered half-space. By means of second-order paraxial approximations of the exact dynamic stiffness for a half-space in plane-strain and antiplane-shear conditions, the corresponding approximation for general three-dimensional wave motion in a Cartesian coordinate system is obtained and transformed in terms of cylindrical coordinates. Using the paraxial approximations, the half-space finite element and consistent transmitting boundary are formulated in a cylindrical coordinate system. The development is verified by comparison of dynamic compliances of rigid circular foundations with available published results. Examination of the advantage of the paraxial condition vis-á-vis the fixed condition shows that the former achieves substantial gain in computational effort. The developed half-space finite element and transmitting boundary can be employed for accurate and effective analysis of foundation dynamics and soil–structure interaction in a porous layered half-space.     

地震波场的复射线展开

利用复射线理论将惠更斯原理推广到复空间,发展出一种简便的波场计算方法.本文以均匀柱面波在平界面上的反射场为例,用复射线展开法计算了时谐地震波及瞬变地震波在部分反射、全反射和临界角过渡区的场分布,并与平面波谱积分结果进行了比较.结果表明,用复射线展开法可以自动平滑过渡区反射场的奇异性,因而该方法可直接处理焦点、焦散线和全反射过渡区等复杂的波场分析与综合问题,并可简化数值计算过程.     

Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm

We present a parallel hybrid algorithm based on pseudospectral method (PSM) and finite difference method (FDM) for two-dimensional (2-D) global SH-wavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corresponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.     

Vibration of surface foundations of arbitrary shapes

Presented is a systematic procedure for generating impedance (or compliance) matrices for foundations with arbitrary shapes, resting on an elastic half-space medium. A technique to decompose prescribed harmonic tractions on the half-space medium is employed to solve analytically the differential wave equations in cylindrical coordinates. However, the interaction stresses due to the vibration of a foundation with arbitrary shape are described in rectangular coordinates, and assumed to be piecewise constant in the region of the arbitrary shape. A coordinate transformation matrix is introduced for the piecewise constant tractions in order to use the solution of the differential wave equations in cylindrical coordinates. Finite element modelling is assumed in rectangular coordinates for the foundation itself. The impedance matrix is then obtained for the finite element model, using a variational principle and the reciprocal theorem. A simple example of a rigid square plate resting on a half-space medium and subjected to vertical excitation is used to demonstrate the efficiency and effectiveness of the procedure. Some numerical aspects are investigated and some possible extensions of the procedure are also discussed.