合成地震图的广义反射透射系数矩阵方法

简述了在均匀层状介质中，合成地震图的广义反射透射系数矩阵方法，并通过与有人的计算结果对比，检验了经过修改调试后的合成地震图的计算程序。计算与对比表明，经修改调试后的程序对于各种因素，包括震源类型、震源深度、震中、频率和慢度都有较好的普适性，适用于实际的地震波形研究。     

An efficient approach for dynamic impedance of surface footing on layered half-space

A semi-analytical model for the evaluation of dynamic impedance of rigid surface footing bonded to multi-layered subsoil is proposed. The technique is based on the dual vector form of wave motion equation and Green's influence function of subdisk for horizontally layered half-space. The multi-layered half-space is divided into a quite large number of mini-layers and the precise integration method (PIM) is introduced for the numerical implementation. The PIM is highly accurate for solving sets of first-order ordinary differential equations with specified two-end boundary conditions. It can produce numerical results of Green's influence functions up to the precision of computer used. The dual vector form of wave motion equation makes the combination of two adjacent mini-layers/layers very easy. As a result, the computational effort for the evaluation of Green's influence function of the multi-layered half-space is reduced to a great extent. In order to satisfy the mixed boundary condition at the surface, the footing–soil interface is discretized into a number of uniformly spaced subdisk-elements. Comparisons illustrating the efficiency and accuracy of the proposed approach are made with a number of solutions available in the literature.     

Matrix impedance in the problem of the calculation of synthetic seismograms for a layered-homogeneous isotropic elastic medium

A new method is proposed for calculating synthetic seismograms caused by a force in a plane-parallel medium consisting of homogeneous elastic isotropic layers. The matrix impedance, i.e., the matrix function of depth, by which motion vector must be multiplied in order to obtain the stress vector, is introduced for solving a system of ordinary differential equations with respect to the motion-stress vector, which appears during the separation of variables. An independent nonlinear equation is obtained for the impedance. The propagator for the motion vector is constructed with the aid of the impedance. The closed analytical formulas, which do not contain any exponents with positive indices, are obtained both for the impedance and for the motionvector propagator. The algorithm for the calculation of seismograms, free of limitations on the number and thickness of layers, as well as on the frequency range of interest, is constructed on the basis of these formulas. The algorithm is tested with the aid of an analytical solution.     

表面矩形基础阻抗函数的集中参数模型

本文在系统地研究了Ｗｏｌｆ二阶集中参数模型的基础上,完整地给出了表面刚性矩形基础阻抗函数的集中参数模型,同以往的工作相比,本文给出的集中参数模型能在更宽的频段上反映精确解的变化。     

Time domain method for calculating free field motion of a layered half-space subjected to obliquely incident body waves

In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.     

Three-dimensional response of a layered cylindrical valley embedded in a layered half-space

An indirect boundary integral method to obtain the three-dimensional response of an infinitely long, layered, viscoelastic valley of arbitrary cross-section embedded in a layered viscoelastic half-space is presented. The valley is excited by homogeneous plane waves impinging at an oblique angle with respect to the axis of the valley. The method and associated computer programs are tested by comparison with available results in the limiting two-dimensional case of incidence normal to the axis of the valley. Additional comparisons with previous three-dimensional results obtained by a hybrid finite element-boundary integral method for cylindrical valleys subjected to obliquely incident waves show large differences. However, the results obtained here for an infinitely long valley appear to be in some agreement with earlier results for an elongated prolate semi-ellipsoidal valley and with results obtained by a discrete wavenumber boundary element approach. An extensive bibliography on the dynamic response of valleys is also presented.     

A precise integration approach for dynamic impedance of rigid strip footing on arbitrary anisotropic layered half-space

The precise integration method (PIM) is proposed for the dynamic response analysis of rigid strip footing resting on arbitrary anisotropic multi-layered half-space. In the frequency domain, the governing equation of wave motion is converted into dual vector form of first-order ordinary differential equations which is solved by PIM. Each layer is divided into a large number (say, 2^N) of mini-layers of equal thickness, within which characteristic matrices are assumed to vary following the Taylor series expansion to the fourth order. As a result, any desired accuracy of the displacements and stresses can be achieved by PIM. In addition, dual vector form equation makes it quite easily to combine two adjacent mini-layers into a new one. Each pass of combination reduces the total number of mini-layers by a half. The computational effort for the evaluation of the dynamic impedance of rigid strip footing can be reduced to a great extent. Numerical examples are provided to validate the efficiency and accuracy of the proposed approach.     

Stochastic dynamic analysis of a layered half-space

A stochastic thin-layer method is developed for the analysis of wave propagation in a layered half-space. A random field of shear moduli in the layered system is considered in terms of multiple correlated random variables. Expanding the random moduli and uncertain responses by means of Hermite polynomial chaos expansions and applying the Galerkin method in the spatial as well as stochastic domains, stochastic versions of thin-layer methods for a layered half-space in plane strain and antiplane shear are obtained. In order to represent the infinite half-space, continued-fraction absorbing boundary conditions are included in the thin-layer models of the half-space. Using these stochastic methods, dynamic responses of a layered half-space subjected to line loads are examined. Means, coefficients of variance, and probability density functions of the half-space responses with a varying correlation coefficient of the shear moduli are computed and verified by comparison with Monte Carlo simulations. It is demonstrated that accurate probabilistic dynamic analysis is possible using the developed stochastic thin-layer methods for a layered half-space.     

含孔隙介质的分层半空间表面瑞利波的衰减特性

以分层半空间内部含有一层孔隙介质为物理模型进行数值计算,研究半空间表面瑞利波的传播和衰减特性.为更加接近实际,结合瑞利波的激发特性,确定了瑞利波的主衰减曲线,并主要以此进行规律分析.针对速度递增和含低速层这两种典型的地质模型,讨论了瑞利波的传播衰减在不同地质模型下的特性,并分析了各自的规律.结果表明,在这两种模型下瑞利波的主衰减曲线都受孔隙介质所处空间位置影响产生比较明显的变化,但衰减系数极大值对应的波长与模型的表层厚度存在较明显的线性对应关系,利用这一关系,可以在实际勘探中快速得到表层介质厚度.另外,通过对比分析还可以看到,瑞利波主衰减曲线随孔隙介质的孔隙度和渗透率的变化都强于主频散曲线的变化,表明衰减曲线对孔隙度和渗透率的变化更加敏感,理论上更加适合进行介质参数反演工作.综合对比结果,我们认为瑞利波主衰减曲线中包含了更丰富的介质参数信息,如果能够有效利用,将可以提高瑞利波勘探的准确性和应用范围.

     

Group theoretic reduction of the electromagnetic impedance matrix for large-contrast symmetric prisms in a layered earth

Calculation using integral equations of the electromagnetic response of a geologic body which is much more conductive than the surrounding media requires the use of both current pulse and current tube basis functions. The impedance matrices for such cases can be large and expensive to form, factor, and solve. However, if the surrounding media is layered and the scatterer is symmetric under symmetry operations which preserve the depth of transformed points, then we can apply group representation theory to drastically reduce storage and computation requirements. I discuss this application of group representation theory in detail, using the symmetry groupC ₂ for purposes of illustration. In a sample calculation for a body which is invariant under the symmetry operations of the groupC ₂ , storage was reduced by a factor of 4, matrix formation time was reduced by a factor of 2, and the sum of matrix factorization and solution times was reduced by a factor of 10.     

Seismic response of foundations embedded in a layered half-space

An integral equation technique to determine the response of foundations embedded in a layered viscoelastic half-space when subjected to various types of seismic waves is presented. The technique is validated by comparison with previous results for rigid hemispherical and cylindrical foundations embedded in a uniform half-space. Illustrative results for rigid cylindrical foundations embedded in layered media are also presented.     

Three-dimensional response of a cylindrical canyon in a layered half-space

An indirect boundary-method formulation to obtain the three-dimensional response of an infinitely long canyon of uniform but arbitrary cross-section cut in a layered viscoelastic half-space is presented. Seismic excitation in the form of plane elastic waves acting at an arbitrary angle with respect to the axis of the canyon is considered. Numerical results for SH-, SV- and P-wave excitation of a circular canyon and of a canyon with a topography similar to that in the vicinity of Pacoima Dam are discussed in some detail.     

Dynamic interaction between rigid foundations in a layered half-space

A computationally efficient boundary integral equation technique to calculate the dynamic response of a group of rigid surface foundations bonded to a layered viscoelastic half-space and subjected to external forces and seismic waves is presented. The technique relies on an iterative scheme which minimizes in-core memory requirements and takes advantage of any geometrical symmetry of the foundations. Extensive results for the case of two rigid square foundations placed at different separations and bonded to a viscoelastic half-space are presented. It was found that the choice of discretization of the foundations has a marked effect on the calculated impedance functions for extremely small separations. Illustrative results for a case of several closely-spaced foundations bonded to a layered half-space are also presented.     

Disturbance due to a point source in a layered half-space

Summary In this paper the displacement components due to a compressional point source in a layered half-space consisting of a liquid layer of finite thickness overlying a semi-infinite solid homogeneous medium have been deduced. Two different cases have been considered. In the first case the source is taken in the liquid layer and in the second case the source is taken in the solid layer. The displacement components in cylindrical coordinates have been obtained in closed form in detail at the interface for all values of distance from the source for the second case.     

Propagation of anSH-torque pulse in a three layered solid half-space

Summary The potential functions caused due to anSH-torque pulse from a point source in any of the layers of a three-layered solid half-space have been evaluated in terms of definite integrals, usingCagniard's method. The solution is given as a double infinite series, of which only a finite number of terms contribute to the disturbance at a definite time, and the number of such terms increases with time.     

Dynamic soil-tunnel interaction in layered half-space for incident P-and SV-waves

The dynamic soil-tunnel interaction is studied by the model of a rigid tunnel embedded in layered half-space,which is simplified as a single soil layer on elastic bedrock to the excitation of P- and SV-waves.The indirect boundary element method is used,combined with the Green's function of distributed loads acting on inclined lines.It is shown that the dynamic characteristics of soil-tunnel interaction in layered half-space are different much from that in homogeneous half-space,and that the mechanism of soil-tunnel interaction is also different much from that of soil-foundation-superstructure interaction.For oblique incidence,the tunnel response for in-plane incident SV-waves is completely different from that for incident SH-waves,while the tunnel response for vertically incident SV-wave is very similar to that of vertically incident SH-wave.     

层状半空间中洞室对入射平面P波的放大作用

文中利用间接边界元法,在频域内求解了层状弹性半空间中洞室对入射平面P波的放大作用问题。通过自由场反应分析,求得假想洞室边界上各单元的应力响应。在洞室边界各个单元上施加虚拟分布荷载,求得位移和应力的格林函数。根据应力边界条件确定虚拟分布荷载,将自由场位移响应和虚拟分布荷载产生的位移响应叠加起来,即得到问题的解答。比较了层状半空间和均匀半空间中洞室对入射平面P波的放大作用,结果表明,层状半空间情况有可能导致较大的地表位移幅值,尤其是对于较低频率入射波。     

层状半空间中洞室对平面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.