Propagation of magneto-elastic waves in a perfectly conducting plate extending to infinity in vacuum

Summary The effects of a uniform external magnetic field on the propagation of waves in a homogeneous, infinitely conducting flat plate with free boundaries have been studied. It has been found that in general all the three types of waves —P, SV andSH waves—are coupled and the influence may be more pronounced in coupling the symmetric and antisymmetric types of motions in every mode.When the magnetic field is parallel to the plane faces and transverse to the direction of wave propagation, the shear wave polarized parallel to the field is purely elastic whereas the coupledP andS V waves are magnetoelastic and exhibit dispersion strikingly similar to the non-magnetic case, provided the electro-magnetic radiation into the surrounding free space is neglected.The results reported in an earlier communication [1]²) are also confirmed.     

重力波在中层大气温度波导中的传播模式研究

本文给出了重力波在中层大气温度波导中的导制传播模型,并在此模型的基础上详细讨论了重力波部分导制传播下的对称模式与非对称模式,导出了不同模式下相应的特征函数和色散方程,进一步用离散的方法对两类色散方程进行了求解;同时还利用二维全隐欧拉格式（FICE）对重力波在温度波导中的传播进行了模拟,模拟的结果也成功地展现了对称与非对称两种传播模式.研究表明,下边界的扰动能量在向上传播进入波导区域后被俘获,形成导制传播.不同周期的初始扰动,在波导内均会形成对称与非对称形式两种模式的导制传播,由于两者的行进速度不一致,最终会引起两种不同模式的分离.数值模拟中重力波的水平行进速度与线性模型预测值非常接近.波导中不同模式下重力波的水平波长与初始扰动的水平波长非常一致,然而波导区域内重力波的频率与初始扰动的频率无关,频率不同的初始扰动会激发出相同频率的重力波对称与非对称导制传播模式.这表明在确定的温度波导中,水平波数才是决定重力波传播特性的决定因素.进一步的分析显示,初始扰动的水平波数-频率分布越接近完全导制传播的色散关系时,温度波导中更易于生成以该种模式部分导制传播的重力波.     

Propagation of normal modes due to impulsive loading to 3-D medium on a rigid basement

Generally, waves in layered media include several surface wave modes showing dispersive properties. These properties often complicate the understanding of the propagation characteristics of elastic waves. In this paper, not only a method for decomposing transient elastic waves into normal modes is presented but also investigations are applied to decomposed normal modes. This method is developed by means of DWFE and the modal superposition method. In the formulation process, transfer functions describing the relationship between the amplitude of normal modes and that of force density are defined. It was found that the propagation properties of the normal modes were well explained by the dispersion curves and transfer functions. From investigations of the normal mode propagation, properties of each normal mode can be discussed.     

含偏心钻铤充液井孔中正交偶极子声波测井的数值模拟研究

采用三维有限差分方法模拟了正交偶极子声源在含偏心钻铤的充液井孔中激发的声场,研究了钻铤偏心对模式波的种类、激发幅度、以及频散特征的影响.研究结果表明,钻铤偏心导致偶极子声源激发的声场的模式不唯一,观察到了除偶极模式外的单极模式波和四极模式波;钻铤偏心导致偶极模式波出现分裂现象,尤其是快速地层F2模式和慢速地层弯曲波的稍高频率的部分,且快、慢波所对应的两个方位为偏心的方位和与偏心方位垂直的方向;井孔折射横波以及快速地层F2模式的低频部分的慢度基本未受到钻铤偏心的影响,仍然能够正确反映地层的横波慢度及各向异性;对于本文研究的慢速地层井孔模型,当偏心距离l小于等于0.01 m时,弯曲波的慢度和激发幅度受钻铤偏心的影响很小,从快、慢弯曲波中提取的快、慢横波慢度基本能够反映地层的各向异性特征.     

基于黏弹滑移界面理论揭示套管中模式波的传播特征

目前固井质量评价主要是利用套管中模式波的幅度或衰减变化判断套管与水泥间的胶结状况.在超声测量频段（70~500 kHz）,现有声波测井仪器工作时在套管中激发的模式波类似薄板中的对称兰姆波（拉伸波）或反对称兰姆波（弯曲型兰姆波）.本文将黏弹滑移界面理论应用到套管和水泥的耦合界面,探讨了套管中的准兰姆波（包括拉伸波和弯曲型兰姆波）与套管后物质的耦合方式,拉伸波与套后物质主要以剪切刚度耦合,即便是0.01 mm的微环也会使得拉伸波的衰减明显降低,接近套后是流体时的响应特征;而弯曲型兰姆波与套管后物质主要以法向刚度耦合,这使得在套后胶结轻质水泥时弯曲型兰姆波对微环不敏感,在套后耦合常规水泥时微环的存在又会使得弯曲型兰姆波的衰减明显增强,衰减的增加说明其与套后物质的声耦合更紧密,也表明微环与常规水泥的等效声阻抗适当降低,更有利于弯曲型兰姆波传播时与套后介质实现法向刚度的匹配.拉伸波和弯曲型兰姆波对同一胶结状况的不同响应特征也为固井质量综合评价时区分不同的水泥胶结状况提供了可能.     

Quadratic normal moveouts of symmetric reflections in elastic media: A quick tutorial

The design of reflection traveltime approximations for optimal stacking and inversion has always been a subject of much interest in seismic processing. A most prominent role is played by quadratic normal moveouts, namely reflection traveltimes around zero-offset computed as second-order Taylor expansions in midpoint and offset coordinates. Quadratic normal moveouts are best employed to model symmetric reflections, for which the ray code in the downgoing direction coincides with the ray code in the upgoing direction in reverse order. Besides pure (non-converted) primaries, many multiply reflected and converted waves give rise to symmetric reflections. We show that the quadratic normal moveout of a symmetric reflection admits a natural decomposition into a midpoint term and an offset term. These, in turn, can be be formulated as the traveltimes of the one-way normal (N) and normal-incidence-point (NIP) waves, respectively. With the help of this decomposition, which is valid for propagation in isotropic and anisotropic elastic media, we are able to derive, in a simple and didactic way, a unified expression for the quadratic normal moveout of a symmetric reflection in its most general form in 3D. The obtained expression allows for a direct interpretation of its various terms and fully encompasses the effects of velocity gradients and Earth surface topography.     

Surface waves in isotropic,laterally inhomogeneous media

Summary Two different viewpoints of the phase velocities of the elastic surface waves in isotropic, laterally inhomogeneous media have led to inconsistent results. Arguments in terms of surface wave modes give the conclusion that the phase velocity is independent of the propagation direction, while the outcome of calculations based on a constructive interference of body waves in a surface layer is that the phase velocity is dependent on the propagation direction. Both arguments are summarized and an error in the calculations giving dependence is pointed out. The calculations and observations of surface wave amplitude changes in laterally inhomogeneous media are also summarized.     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

REFLECTION OF ELASTIC WAVES FROM PERIODICALLY STRATIFIED MEDIA WITH INTERFACIAL SLIP*

A periodically stratified elastic medium can be replaced by an equivalent homogeneous transverse isotropic medium in the long wavelength limit. The case of a homogeneous medium with equally spaced parallel interfaces along which there is imperfect bonding is a special instance of such a medium. Slowness surfaces are derived for all plane wave modes through the equivalent medium and reflection coefficients for a half-space of such a medium are found. The slowness surface for the SH mode is an ellipsoid. The exact solution for the reflection of SH-waves from a half-space with parallel slip interfaces is found following the matrix method of K. Gilbert applied to elastic waves. Explicit results are derived and in the long wavelength limit, shown to approach the results for waves in the equivalent homogeneous medium. Under certain conditions, a half-space of a medium with parallel slip interfaces has a reflection coefficient independent of the angle of incidence and thus acts like an acoustic reducing mirror. The method for the reflection of P- and SV-waves is fully outlined, and reflection coefficients are shown for a particular example. The solution requires finding the eigenvalues of a 4 × 4 transfer matrix, each eigenvalue being associated with a particular wave. At higher frequencies, unexpected eigenvalues are found corresponding to refracted waves for which shear and compressional parameters are completely coupled. The two eigenvalues corresponding to the transmitted wavefield give amplitude decay perpendicular to the stratification along with up- and downgoing phase propagation in some other direction. Much of this work was performed while the author was at the Department of Geophysics and Planetary Sciences, Tel-Aviv University, Ramat-Aviv, Israel. The author is grateful for illuminating discussions with K. Helbig and K. Gilbert.     

Analysis of effects of active sources on observed phase velocity based on the thin layer method

In the wave field induced by active sources, the observed phase velocity of surface waves is influenced by both mode incompatibility (i.e. non-planar spread of surface waves is idealized as plane waves) and body waves. Effects of sources are usually investigated based on numerical simulations and physical models. Several methods have been proposed to mitigate the effects. In application, however, these methods may also have difficulties since the energy of the body waves depends on soil stratification and parameters. There are multiple modes of surface waves in layered media, among which the higher modes dominate the wave field for soils with the irregular shear velocity profiles. Considering the mode incompatibility and the higher modes, we derive analytical expressions for the effective phase velocity of the surface waves based on the thin layer stiffness method, and investigate the effects of the body waves on the observed phase velocity through the phase analysis of the vibrations of both the surface waves and the body waves. The results indicate that the effective phase velocity of the surface waves in layered media varies with the frequency and the spread distance, and is underestimated compared to that of the plane surface waves in the spread range less than about one wavelength. The oscillations that appeared in the observed phase velocity are due to the involvement of the body waves. The mode incompatibility can be ignored in the range beyond one wavelength, while the influence range of the body waves is far beyond one wavelength. The body waves have a significant influence on the observed phase velocity of the surface waves in soils with a soft layer trapped between the first and the second layers because of strong reflections.     

Wave dispersion studies in granular media by analytical and analytical–numerical methods

The phenomenon of wave dispersion in dry sand is studied both by purely analytical studies and by analytical–numerical experiments on the basis of gradient elastic and viscoelastic material models. These material models are employed in order to simulate the microstructural characteristics of dry sand. The analytical studies treat the material body as a one-dimensional (for the viscoelastic case) and three-dimensional (for the gradient elastic case) and for both material models provide explicit expressions for the velocity of propagation of harmonic compressional (P) and shear (S) waves. These velocities are found to be functions of frequency, i.e., dispersive. The analytical–numerical studies treat the material as a one-dimensional one and try to simulate P and S wave propagation along the axial direction of cylindrical dry sand specimens. Thus, a sinusoidal pulse with a specific frequency is applied at one end of the specimen and the response is determined at some other point by solving a transient dynamic boundary value problem with the aid of a numerical Laplace transform. This analytical–numerical experiment is repeated for various frequencies. Thus, one determines the velocities of P and S waves as functions of frequency, thereby proving again that wave propagation in dry sand is dispersive.     

套后介质声学性质对套管中三种模式波影响的对比分析

由于套管的厚度小于目前声波测井仪器工作频率下的波长,声波仪器工作时可在套管中激发Lamb波或SH模式波.目前评价水泥胶结质量的声波测井仪器主要是利用套管中的对称Lamb波(拉伸波,S0模式)或反对称Lamb波(泄漏弯曲波,A0模式)的相对幅度或衰减评价套管与水泥之间的胶结程度.Lamb波可看作是纵波和横波在套管边界的耦合形成的,在套管中还存在另一种非耦合的SH横波,偏振方向与套管轴向平行沿着周向传播.本文基于定向声源辐射技术,实现了对称Lamb波、反对称Lamb波的分波模拟技术,数值研究了套管后耦合不同声阻抗介质时,对称Lamb波、反对称Lamb波以及SH横波的测井响应,通过频散衰减曲线以及模式波的波形响应幅度的直观显示表明,由于套管中的SH波仅向固体水泥中辐射剪切波,这一传播特征与套管中的Lamb波相比可以极大提高套管-水泥界面胶结状况的测量灵敏度.Lamb波的固有缺陷是套管与两侧介质的较大声阻抗反差使得在井中测量套管后面水泥界面的灵敏度受到相当大的限制.     

Alfven波在低纬电离层中的传播研究

Alfven波在低纬地区电离层的传播有其特殊性,一方面,低纬地区同样存在Alfven速度梯度的巨大变化,导致电离层Alfven谐振器(Ionospheric Alfven resonator, IAR)的形成;另一方面,由于在低纬地区磁倾角很小,所以剪切Alfven波在传播的过程中纬度方向跨度很大,不同纬度电离层参数将共同对其产生影响;并且,由于电离层水平分层,故磁力线与电离层不正交.本文选取双流体力学模型,在忽略场向电场的条件下,利用非正交坐标系,结合IRI07模型与MSISE00模型模拟低纬地区Alfven波的传播,得到其反射及耦合特性.结果表明,低纬地区同样存在电离层Alfven谐振现象,由耦合产生的压缩模有向磁赤道方向传播的趋势,夜间电离层状态相对于白天更适合IAR的形成,谐振频率沿磁力线L值增大单调递增.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

地球磁层顶的漂移动力学Alfvn波不稳定性及其反常输运效应

本文用回旋动力学研究了磁层顶等离子体低频漂移动力学不稳定性.在β≥1附近发现两支不稳定的漂移动力学Alfvèn模（DKA）.它们可在▽T（∥-▽n）≠0时激发,速度剪切提供主要的自由能源.两支DKA均具有非零的（δE_y,δE_∥）和（δB_x,δB_∥）.在湍动的非线性饱和状态下,δB_x的起伏可导致很强的反常输运,当|δB_x|=1nT时,D_⊥可达到10⁹m²/s的数量级.因此,DKA可能在太阳风-磁层耦合过程中起重要作用.     

Spatial dispersion of surface waves at the free surface of a general anisotropic elastic medium

A new technique relates the wave velocity of the surface waves in anisotropic elastic medium to its elastic constants. Anisotropic propagation of surface waves is studied in a half-space occupied by a general anisotropic elastic solid. The phase velocity expressions of quasi-waves, in three-dimensional space, are used to derive the secular equation of surface waves. The complex secular equation is resolved, analytically, into real and imaginary parts and is then solved, numerically, for phase velocity along a given phase direction on the surface. The complete procedure is thus analogous to the one used for conventional Rayleigh waves in isotropic medium. A non-linear equation relates the ray direction of the surface waves to its phase direction on the (plane) surface of the medium. The analytical differentiation of secular equation yields the directional derivative of phase velocity. This derivative is used to calculate the wave velocity of surface waves. Spatial variations of phase velocity, wave velocity and ray direction over the free plane surface are plotted for the numerical models of crustal rocks with orthorhombic, monoclinic and triclinic anisotropies.     

横向各向同性地层中随钻声波测井模式波分析

针对横向各向同性地层随钻声波测井模型,通过模式分析的方法,考察了快速地层和慢速地层井孔内随钻单极子、偶极子和四极子声源激发的斯通利波、弯曲波和螺旋波的相速度频散和激发强度特征,计算了这些模式波对于地层弹性常数的灵敏度,并与电缆测井中的情况进行了比较.结果表明:随钻斯通利波在低频时对地层弹性常数中c₆₆的灵敏度较电缆测井中有了很大提高,可用于反演地层水平向横波速度;随钻偶极子最低阶弯曲波在低频时不能用于直接获取地层横波信息,但在慢速地层中频率较高(例如6 kHz)时却可以间接得到地层垂直向横波速度;随钻四极子螺旋波的特征与电缆测井中的类似,可用于获取地层垂直向横波速度.     

Propagation of plane waves in an unbounded elastic medium with Cauchy's initial stress

Summary Cauchy theory of initial stress has been applied to investigate the influence of uniform initial stress on the propagation of plane waves in an unbounded elastic medium. It is found that if the initial principal stresses are unequal, the velocity depends on the direction in which the wave propagates, whereas when the initial principal stresses are equal, the velocity is independent of the direction of wave propagation.     

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.     

In-plane free-field response of actual sites

To study the characteristic features of the in-plane free-field response, two actual sites of nuclear power plants, a soft and a rock site, are analysed, by varying the location of the control point and the nature of the wave pattern. Harmonic and transient seismic excitations are examined. The conclusions reached in the vast parametric study of Reference 6 are confirmed for the actual sites. These apply to the range of possible apparent velocities and the associated motions, the spatial variations with depth and in the horizontal direction. If only one component of the control motion, e.g. the horizontal, is matched, then it can be associated either with a body wave or with a surface wave. In the latter case, a specific mode is used up to the frequency at which the next higher mode starts, since for a given frequency, the higher modes attenuate less. The other component of the motion follows. If both components are prescribed, the motion can be interpreted as arising from a combination of a P- and an SV-wave (with a common apparent velocity). Surface waves alone cannot be used to match both components. A body wave has to be included, at least up to the frequency at which the second mode starts. For the soft site, the surface waves decay significantly, especially in the range of higher frequencies, where the apparent velocity is considerably smaller than the shear-wave velocity of the rock. For this site, it seems sufficient to examine only (extremely shallow) body waves. For a rock site, however, Rayleigh-waves exist which attenuate little, leading to smaller apparent velocities than extremely shallow body waves.