频散介质中地质雷达波传播的数值模拟

地质雷达所探测的地球介质常常具有频散性．为了研究地质雷达在频散介质中的探测能力,提出了频散介质中时间域有限差分法计算麦克斯韦方程的方法,给出了满足Debye关系的频散介质中的电位移和磁场的迭代算法,以及由电位移计算电场的算法．只有在电场计算时才用到介质的物性参数．提出一种新的吸收边界条件的算法,通过增加假想的介电常数和磁导率,实现了吸收层中波的无反射衰减,克服了以往Berenger完全匹配层计算时对场进行分裂带来的麻烦,从而提高了计算效率．计算实例表明,频散介质中电磁波的衰减更快,测量信号变得很弱．     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

地中海直布罗陀海峡附近内孤立波的地球物理特征

本文主要利用地震海洋学方法研究地中海直布罗陀海峡附近内孤立波的结构特征,此处内孤立波为第一模态下沉型,为中幅度和大幅度内孤立波,垂向振幅最大可达74.5 m,振幅随深度增加呈增大趋势,传播速度随振幅增大而增大,可以确定"真"最大振幅位置位于密跃层附近.由于类多普勒效应和孤立波与测量船之间存在夹角的原因,从地震剖面上得到的为视半高宽参数,需要进行校正后才能得到比较真实的半高宽参数,校正后半高宽最高可达到1721.8 m,但是校正后的半高宽与理论结果有些差距,这可能与内孤立波的发育稳定程度有关.随着内孤立波包不断向东运动,整体波宽变大,垂向速度变小.本文将地震海洋学方法拓展应用于地中海区域内孤立波分析,进一步证明了利用地震海洋学方法研究海水运动的可行性.     

切变基本纬向流中非线性赤道Rossby长波

为了解决观测和理论研究中的一些问题以及更好地了解热带大气动力学 ,有必要进一步研究基本气流的变化对大气中赤道Rossby波动的影响 .本文研究分析基本气流对赤道Rossby长波的影响 ,利用一个简单赤道 β平面浅水模式和摄动法 ,研究纬向基本气流切变中非线性赤道Rossby波 ,推导出在切变基本纬向流中赤道Rossby长波振幅演变所满足的非线性KdV方程并得到其孤立波解 .分析表明 ,孤立波存在的必要条件是基本气流有切变 ,而且基流切变不能太强 ,否则将产生正压不稳定 .     

固体介质中弹性波传播机制的细胞自动机有限深势阱模型

分析了采用细胞自动机研究波动问题的建模方法，针对一维、均匀、各向同性固体介质中弹性纵波的微观机制，借用一经典弹簧振子模型、细胞自动机格子气模型，以及量子力学中的无限深势阱模型，建立了一个细胞自动机有限深势阱模型，从量子力学角度出发，基于介观物理和纳米概念，以微观精子的德布罗意假设为基础，利用薛定谔方程，讨论了该模型中粒子（分子组）的振动速度与粒子物质波波速之间的联系，给出了模型中的波动方程，得出ζ＝Vp（ζ为粒子振动速度，Vp为物质波纵波波速）。同时还讨论了模型中粒子的大小和能量传递问题，引入引力场，得出了能量及引力势的量子化条件，另外，对声波速度、格子气粒子振动速度和本文模型中分子组振动速度进行了比较；还对本文模型中的粒子能量分布作了分析。     

二维横各向同性弹性随机介质中的波场特征

本文通过交错网格有限差分正演．模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录．为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值)．这样,对应每一个横各向同性弹性随机介质模型．均可计算得到15个不同的波场特征量．我们通过在二维横各向同性弹性随机介质中的正演模拟．研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点．证实了在随机介质模型中．各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论．     

安达曼海南部内孤立波的生成与传播的数值模拟

大量的内孤立波在安达曼南部被观测到.文章利用二维的MITgcm模式来研究该海域内孤立波的动力机制并探索底地形及潮流驱动力对内孤立波生成及传播的影响.结果 表明,大振幅下凹型内孤立波主要由潮流越过格雷特海峡附近的海脊(标记为海脊A)生成,生成机制为背风波机制.无畏舰岸坡(Dreadnought Bank)本身不能激发生成...     

双相各向异性介质中弹性波传播特征研究

随着地震工程和能源地震勘探的深入发展,人们所遇到的地下介质愈来愈复杂．常规的各向异性介质理论或双相各向同性介质理论难以精确描述含流体的各向异性介质,如裂缝性气藏、含水页岩等．本文以Biot双相各向异性介质理论为基础,利用弹性平面波方程,推导出了任意双相各向异性介质中弹性波的Christoffel方程．根据Christoffel方程,计算并分析了频率对双相横向各向同性介质中弹性波的相速度、衰减、双相振幅比和偏振特征的影响．结果表明,在4类波(快纵波、慢纵波、快横波和慢横波)中,频率对慢纵波影响最大;当耗散很大时,快纵波、快横波和慢横波的流固相振幅比值近似为1．对偏振特征分析的结果表明,在双相各向异性介质中,弹性波的固相位移偏振方向与流相位移偏振方向将不再保持同向或反向,而是呈不同大小的夹角．     

Modelling Internal Solitary Waves in the Coastal Ocean

In the coastal oceans, the interaction of currents (such as the barotropic tide) with topography can generate large-amplitude, horizontally propagating internal solitary waves. These waves often occur in regions where the waveguide properties vary in the direction of propagation. We consider the modelling of these waves by nonlinear evolution equations of the Korteweg–de Vries type with variable coefficients, and we describe how these models are used to describe the shoaling of internal solitary waves over the continental shelf and slope. The theories are compared with various numerical simulations.     

Solitary Rossby Waves in Zonally Varying Jet Flows

The dynamics of solitary Rossby waves (SRWs) embedded in a meridionally sheared, zonally varying background flow are examined using a non-divergent barotropic model centered on a midlatitude g -plane. The zonally varying background flow, which is produced by an external potential vorticity (PV) forcing, yields a modified Korteweg-de Vries (K-dV) equation that governs the spatial-temporal evolution of a disturbance field that contains both Rossby wave packets and SRWs. The modified K-dV equation differs from the classical equation in that the zonally varying background flow, which varies on the same scale as the disturbance field, directly affects the disturbance linear translation speed and linear growth characteristics. In the limit of a locally parallel background flow, equations governing the amplitude and propagation characteristics of SRWs are derived analytically. These equations show, for example, that a sufficiently large (small) translation speed and/or a sufficiently weak (strong) background zonal shear favor transmission (reflection) of the SRW through (from) the jet. Conservation equations are derived showing that time changes in the domain averaged amplitude ("mass") or squared amplitude ("momentum") are due to zonal variation in both the linear, long-wave phase speed and linear growth; dispersion and nonlinearity do not affect the "mass" or "momentum". Provided (1) the background PV forcing is sufficiently small, or (2) the background PV forcing is meridionally symmetric and the disturbance is a SRW, the dynamics of the disturbance field is Hamiltonian and mass and energy are thus conserved. Numerical solutions of the K-dV equation show that the zonally varying background flow yields three general classes of behavior: reflection, transmission, or trapping. Within each class there exists SRWs and Rossby wave packets. SRWs that become trapped within the zonally localized jet region may exhibit the following behaviors: (1) an oscillatory decay to a steady state at the jet center, (2) the creation of additional SRWs within the jet region, or (3) a steady-state wherein the solution has a smoothed step-like structure located downstream along the jet axis.     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

Resonant forcing of coastally trapped waves in a continuously stratified ocean

In a previous paper (Grimshaw, 1987) the resonant forcing of coastally trapped waves was discussed in the barotropic case. In order to extend that theory to more realistic situations, we have considered the analogous theory whereby a longshore current interacts with a longshore topographic feature, or the forcing is due to longshore wind stress, for the case of the continuously stratified ocean. As in the previous theory, near resonance, when a long-wave phase speed is close to zero (in the reference frame of the forcing), the wave motion is governed by a forced evolution equation of the KdV-type. The behaviour of the wave field is characterized by three parameters representing the bandwidth for resonance, the forcing amplitude and the dissipation. We have evaluated these parameters in various practical cases, and found that the bandwidths, which scale with ^1/2 when the forcing has dimensionless amplitude , can often be quite broad. Typically the second, third, or higher, modes may be resonant. Concurrently, the dissipation is also usually significant, leading to a steady state balance between the forcing, dissipation and nonlinear terms.     

Inhomogeneous Harmonic Plane Waves in Viscoelastic Anisotropic Media

In viscoelastic media, the slowness vector p of plane waves is complex-valued, p = P + iA. The real-valued vectors P and A are usually called the propagation and the attenuation vector, repectively. For P and A nonparallel, the plane wave is called inhomogeneousThree basic approaches to the determination of the slowness vector of an inhomogeneous plane wave propagating in a homogeneous viscoelastic anisotropic medium are discussed. They differ in the specification of the mathematical form of the slowness vector p. We speak of directional specification, componental specification and mixed specification of the slowness vector. Individual specifications lead to the eigenvalue problems for 3 × 3 or 6 × 6 complex-valued matrices.In the directional specification of the slowness vector, the real-valued unit vectors N and M in the direction of P and A are assumed to be known. This has been the most common specification of the slowness vector used in the seismological literature. In the componental specification, the real-valued unit vectors N and M are not known in advance. Instead, the complex-valued vactorial component p of slowness vector p into an arbitrary plane with unit normal n is assumed to be known. Finally, the mixed specification is a special case of the componental specification with p purely imaginary. In the mixed specification, plane represents the plane of constant phase, so that N = ±n. Consequently, unit vector N is known, similarly as in the directional specification. Instead of unit vector M, however, the vectorial component d of the attenuation vector in the plane of constant phase is known.The simplest, most straightforward and transparent algorithms to determine the phase velocities and slowness vectors of inhomogeneous plane waves propagating in viscoelastic anisotropic media are obtained, if the mixed specification of the slowness vector is used. These algorithms are based on the solution of a conventional eigenvalue problem for 6 × 6 complex-valued matrices. The derived equations are quite general and universal. They can be used both for homogeneous and inhomogeneous plane waves, propagating in elastic or viscoelastic, isotropic or anisotropic media. Contrary to the mixed specififcation, the directional specification can hardly be used to determine the slowness vector of inhomogeneous plane waves propagating in viscoelastic anisotropic media. Although the procedure is based on 3 × 3 complex-valued matrices, it yields a cumbersome system of two coupled equations.     

SH波在表面多层介质中传播的精确模拟

针对地震横波在地表低速层内的振幅放大效应问题，提出了一种模拟SH波在地表层状介质中传播的递推算法，并用它模拟了新西兰Alfredton盆地A10场址的SH波地震动响应特性。这个方法适用于具线性吸收性质的粘弹性介质。由于方法不受介质层厚薄制约，层厚可以无限薄化，实践上可以用许多薄层逼近的办法来模拟纵向上任意变化的连续介质。通过求取不同频率不同波数平面简谐波解并按实际问题的加权迭加可求解具特定波形和传播方向组合的任意SH波场。此方法在计算上具有解析解特有的精确性，稳定性和方便性，特别适用于模拟薄层介层，次波长现象及需要进行大量而又精确模拟计算的情形。     

Asymptotic Ray Theory for Seismic Surface Waves in Laterally Inhomogeneous Media

A recurrence procedure is outlined for constructing asymptotic series for surface wave field in a half-space with weak lateral heterogeneity. Both horizontal variations of the elastic parameters and of the wave field are assumed small on the distances comparable with the wavelength. This is equivalent to the condition that the frequency is large. The Surface Wave Asymptotic Ray Theory (SWART) is an analog of the asymptotic ray theory (ART) for body waves. However the case of surface waves presents additional difficulty: the rate of amplitude variation is different in vertical and horizontal directions. In vertical direction it is proportional to the large parameter . To overcome this difficulty the transformation equalizing vertical and horizontal coordinated is suggested, Z = z. In the coordinates x,y,Z the wave field is represented as an asymptotic series in inverse powers of . The amplitudes of successive terms of the series are determined from a recurrent system of equations. Attention is paid to similarity and difference of the procedures for constructing the ray series in SWART and ART. Applications of SWART to interpretation of seismological observations are discussed.