Rayleigh wave and detection of low-velocity layers in a stratified half-space

The excitation and propagation of the guided waves in a stratified half-space and a Rayleigh wave exploration method in shallow engineering seismic exploration are studied in this paper. All the modes of the guided waves are calculated by the bisection method in the case where the low velocity layers are contained in a stratified half-space. Cases when the formation shear wave velocity gradually decreases from the top to the bottom layers are also studied. The dispersion curves obtained in actual Rayleigh wave exploration are usually noncontinual zigzag curves, but the dispersion curves given by the elastic theory for given modes of the guided waves are smooth and continual curves. In this paper, the mechanism of zigzag dispersion curves in Rayleigh wave exploration is investigated and analyzed thoroughly. The zigzag dispersion curves can give not only the possible positions of the low-velocity layers but also the other information on the formation structure (fractures, oil, gas, etc.). It is found that the zigzag dispersion curves of the Rayleigh wave are the result of the leap of the modes and the existence of low velocity layers in a stratified half-space. The effects of the compressional wave velocity, shear wave velocity, and density of each layer on zigzag dispersion curves and the relationship of the low velocity layers to zigzag dispersion curves are also investigated in detail. Finally, the exploration depth of the Rayleigh wave is discussed. The exploration depth of the Rayleigh wave is equal to the wavelength multiplied by a coefficient that is variable and usually given by the work experience and the formation properties of the local work area.     

Perturbation Analysis on Guided Waves in a Fluid-Filled Borehole Surrounded by a Cubic Crystal Anisotropic Medium

The perturbation method is employed to analyse the guided waves in a borehole surrounded by a cubic crystal medium for the first time. The cubic crystal medium is regarded as a reference unperturbed isotropic state added to the perturbation. The dispersion characteristics of Stoneley wave, pseudo-Rayleigh wave, flexural wave, and screw wave are investigated in detail. It is found that dispersion of the guided waves excited by monopole and dipole sources does not depend on the azimuth of the source, whereas the dispersion of screw wave excited by quadrupole source is significantly related to the azimuth of the source. Screw waves propagated along different azimuth in the borehole can be split. This is different from screw waves in transversely isotropic media （hexagonal crystal）, which have been widely studied.     

分层固体板中导波的激发与频散特性

针对无限大弹性分层固体板,研究了结构中导波的频散和激发特性。首先使用传递矩阵法推导分层板模型中导波的频散方程,然后用二分法求取导波各模式的频散曲线,进而分析结构中导波的频散特性。结果表明:在速度递增或递减的分层板中,基阶模式和高阶模式的高频极限分别等于低速层的瑞利波速和横波波速。对于含低速夹层的分层板,所有模式的高频极限都等于低速层的横波速度。在导波激发特性方面,研究了在具有一定宽度的法向力源作用下的分层板中导波各模式在结构中的法向位移谱。发现在速度递增的分层板结构中基阶模式是主导模式,而对于速度递减和含低速夹层模型,主导模式在不同的频段范围内对应不同的导波模式。      

Guided waves in a stratified half-space

The dispersion and excitation mechanisms and the energy distribution of guided waves in a stratified half-space are studied. All possible guided waves excited by a symmetric point source in two or three-layer medium models and their relation to the medium parameters are analyzed in detail. The excitation and propagation characteristics, as well as the energy distribution along the depth direction, of all modes of the surface waves and trapped waves are numerically investigated and analyzed thoroughly not only in the case when the shear wave velocity increases from up to down layers but also when a low-velocity layer is contained in halfspace, especially when the shear wave velocity decreases from up to down layers. It is found that there exist many guided wave modes in the case where the shear wave velocity of each layer increases from up to down layers. However, there is less than one guided wave mode in the case where the shear wave velocity of each layer decreases from up to down layers. The trapped waves exist and propagate along the low-velocity structure in the stratified half-space. It is also found that the characteristic of a mode is related to the source frequency. It is possible that a surface wave at one value of frequency is like a trapped wave at another value of frequency. Finally, the relation of the characteristics of all guided waves (surface waves and trapped waves) to the parameters of media is studied.     

无限流体中孔隙介质圆柱周向导波的传播特性

为研究无限大流体约束的孔隙圆柱中周向导波的传播规律,分析孔隙参数对导波传播特性的影响,建立了无限流体中孔隙介质圆柱的理论模型,利用孔隙介质弹性波动理论,建立了周向导波频散方程,通过数值模拟计算得到无限流体中孔隙介质圆柱的频散曲线,探讨了圆柱半径和孔隙参数对导波传播特性的影响,并对导波的衰减特性进行了分析;通过数值计算,得到了周向导波的时域波形,讨论了孔隙参数对波形的影响.结果表明,孔隙介质圆柱半径的改变影响圆柱尺度,孔隙度的改变影响孔隙介质中体声波的波速,都对周向导波频散曲线产生一定的影响,所得到的频散曲线特征及衰减曲线与时域波形吻合.研究结果对开展无限流体中孔隙介质圆柱的超声无损评价提供了一定的理论参考.     

北极冰-水耦合系统中声波传播特性

利用固体和流体介质中波传播理论,导出了冰-水两层复合结构中导波频散方程.进一步,利用二分法对频散方程进行了数值求解,得到了 w-k频散曲线(w与k分别为圆频率和波数),以及相速度和群速度频散曲线.结果表明:冰-水两层复合结构中导波由具有相同厚度水层和冰层中导波耦合而成,但与水层和冰层中导波频散曲线相比,复合结构中导波频...     

Guided wave dispersion curves for a bar with an arbitrary cross-section,a rod and rail example

Theoretical and experimental issues of acquiring dispersion curves for bars of arbitrary cross-section are discussed. Since a guided wave can propagate over long distances in a structure, guided waves have great potential for being applied to the rapid non-destructive evaluation of large structures such as rails in the railroad industry. Such fundamental data as phase velocity, group velocity, and wave structure for each guided wave mode is presented for structures with complicated cross-sectional geometries as rail. Phase velocity and group velocity dispersion curves are obtained for bars with an arbitrary cross-section using a semi-analytical finite element method. Since a large number of propagating modes with close phase velocities exist, dispersion curves consisting of only dominant modes are obtained by calculating the displacement at a received point for each mode. These theoretical dispersion curves agree in characteristic parts with the experimental dispersion curves obtained by a two-dimensional Fourier transform technique.     

Research on the propagation of guided waves in multi-layered plates

The dispersion characteristics and excitation mechanisms of the guided waves in multilayered plates are studied in this paper. Firstly, the dispersion equation is obtained by the propagator matrix method. Then, the bisection technique is employed to find all the roots of the dispersion equation. The dispersion characteristics of the guided waves are investigated and analyzed. For the multilayered plates in which the S-wave velocity increases or decreases from top to bottom, it is found that the velocity limits in high frequency of the first and high modes are equal to the Rayleigh wave and S-wave velocities of the low-velocity layer, respectively. It is also found that the velocity limits in the high frequency of all modes are equal to the S-wave velocity of the low-velocity layer for the plate with a low-velocity middle layer. The normal displacement spectrum of all the modes excited by the normal force source with a definite width on the surface of the plate is also investigated. It is proved that the dominant mode is the first mode when the S-wave velocity increases from top to bottom layer and the dominant mode is different in different frequency range for the plate with a low-velocity middle layer.     

Ultrasonic guided wave focusing by a generalized phased array

Ultrasonic guided wave focusing by a generalized phased array is studied based on dispersion curves in a multi-layered medium. The different phase of the guided waves with different frequency is added on the excitation signal on each element of the transducer array for focusing. This can be realized by designing a proper excitation pulse based on the dispersion curves of the guided waves for each of the transducer array elements. The numerical simulation results show that the guided waves with different modes, different frequency components, and from different elements of the transducer array can all be focused at the target and focusing is achieved.     

Propagation characteristics of guided waves in a rod surrounded by an infinite solid medium

The dispersion and excitation characteristics of the guided waves in a rod surrounded by an infinite solid medium (cladding) are investigated. First, the bisection technique is employed to find all the roots of the dispersion function on the basis of theoretical analysis and to obtain the complex phase and group velocity dispersion curves of the guided modes. Second, according to their different dispersion characteristics, the guided modes are divided into two categories: normal modes and Stoneley modes. And it is concluded that the normal modes merely exist in the “hard cladding” model in which the cladding’s shear velocity is larger than the rod’s; while the Stoneley modes in cylindrical interface are highly dispersive and merely exist in the model whose acoustical parameters satisfied the existence condition of the Stoneley waves. Third, the seldom discussed issue, the excitation mechanisms of the guided waves, excited by three source models: symmetric point source, axial and radial force sources, are simulated respectively. Attention is paid on the dominant mode which has better excitation sensitivity and the suitable excitation frequency range. Moreover, the propagation characteristics of the Stoneley modes, ignored in previous references, are analyzed and compared with those of the normal modes.     

Calculation of leaky Lamb waves with a semi-analytical finite element method

A semi-analytical finite element method (SAFE) has been widely used for calculating dispersion curves and mode shapes of guided waves as well as transient waves in a bar like structures. Although guided wave inspection is often conducted for water-loaded plates and pipes, most of the SAFE techniques have not been extended to a plate with leaky media. This study describes leaky Lamb wave calculation with the SAFE. We formulated a new solution using a feature that a single Lamb wave mode generates a harmonic plane wave in leaky media. Dispersion curves obtained with the SAFE agreed well with the previous theoretical studies, which represents that the SAFE calculation was conducted with sufficient accuracy. Moreover, we discussed dispersion curves, attenuation curves, and displacement distributions for total transmission modes and leaky plate modes in a single side and both two side water-loaded plate.     

转移矩阵法在负折射率介质材料平板波导中的应用研究

利用严格电磁理论,推导出了适用于负折射率介质材料光波导的转移矩阵,分析讨论了转移矩阵的性质和应用.利用转移矩阵方法,推导出导波层为负折射率介质材料、覆盖层和衬底为右手材料的三层对称介质光波导的本征色散方程.用图解法研究了负折射率介质波导中TE波的异常色散特性.在负折射材料介质波导中没有零阶模,最低阶为1阶模,并且有截止频率,只有波导参量满足一定条件的时候才会存在,导模的横向波数可以为实数和纯虚数,而正折射率介质波导导模的横向波数只能为实数.     

Propagation of guided elastic waves in 2D phononic crystals

The phononic band structure of two-dimensional phononic guides is numerically studied. A plane wave expansion method is used to calculate the dispersion relations of guided elastic waves in these periodic media, including 2D phononic plates and thin layered periodic arrangements. We show that, for any guided elastic wave, Lamb or generalised Lamb modes, stop bands appear in the dispersion curves, displaying a phononic band structure in both cases.     

沿均匀无限介质中固体杆传播的导波特性研究

针对均匀无限固体介质包围的无限长固体杆,研究了沿固体杆传播的导波的频散和激发特性。在理论分析的基础上,用二分法快速求取了导波所有模式的频散曲线,在4种介质参数下,对斯东利模式和简正模式进行了深入的分析和总结;研究了固体杆中导波的激发强度特性,分别在对称点源、轴向力源和径向力源激励下,研究了导波诸模式的激发机理及频谱特性,分析了不同频段导波的主导模式,研究了主导模式的频响特性及位移径向分布关系。      

三层介质超声谐振模式随材料和界面粘接性能变化的演变规律

在用超声波谐振对粘接材料的粘接强度进行无损评估时, 不同模式对粘接强度的敏感程度受到众多因素和参数的影响, 对检测结果的可靠性至关重要. 基于多层介质中声传播和界面弱粘接边界条件的理论模型, 将一个上下非对称的金属-粘接剂-金属三层结构的平面波反射系数函数中的谐振模式看作是上下铝金属层各自的Lamb波频散模式通过夹心粘接剂层相互耦合后叠加组成. 改变影响结构粘接强度的因素, 即粘接剂的性能参数(声阻抗、密度、厚度)和界面切向劲度系数k_t来分析三层结构谐振模式耦合方式的变化,得出结论: 粘接结构粘接性能的变化基本上不改变与被粘铝层相关的固有部分的Lamb波模式, 而它们的耦合模式则在谐振频率上产生平移并会与固有模式进行交换和替代; 不同参数的变化引起的模式演变有各自的规律, 大多可彼此区分.     

Extraction of dispersion curves for waves propagating in free complex waveguides by standard finite element codes

This paper presents a fast and reliable method, for obtaining all the range of dispersion curves for wave propagation usually used in practice, by numerical simulation only, via common commercial finite element codes. Essentially, the method is based on a simple and robust approach, consisting in a few series of modal analyses for a representative part of the inspected structure. In this way, for different wave lengths, one can find the mode shapes and corresponding natural frequencies by solving some real, symmetric and well numerically conditioned eigenvalue problems. The method allows the extraction of propagating modes only and, in spite of not producing continuous dispersion curves, it is not susceptible to aliasing effects, as some similar methods are. Additionally, complete graphical representations of guided waves are possible with some minor calculus effort.     

谱有限元数值分析各向异性复合板中导波色散特性

提出谱有限元方法研究层状各向异性复合板中导波的色散特性和波结构。基于三维弹性动力学方程,用有限元方法离散波导截面,波传播方向的位移用简谐波表示,得到了导波色散的特征方程。分析了单层和双层复合板中导波沿不同方向传播的色散特性和波结构,讨论了双层复合板中层厚比对相速度的影响。数值研究结果表明:导波的对称模态沿纤维方向传播时在较宽的频率范围内保持弱色散状态。双层复合板中导波基本模态的相速度在低频时受层厚比的影响较明显,随着频率的增加趋向于相速度较低的材料。数值模拟结果为导波用于复合材料定量无损检测和性能评价提供理论依据。      

Wave structure analysis of guided waves in a bar with an arbitrary cross-section

Both dispersion curves and wave structures, which are displacement distributions on a bar cross-section, are essential for guided wave NDEs. Theoretical dispersion curves and wave structures for a bar with an arbitrary cross-section are derived in this paper using a special modeling technique called the semi-analytical finite element method (SAFEM). The guidelines for guided wave NDEs of bar-like structures are also shown based on wave structure and modal analysis. First, the relationship between the dispersion curves and their corresponding wave structures were obtained for a square rod. Modes with longitudinal vibration have higher group velocities and torsional modes have constant phase and group velocities. Next, the relationship between the dispersion curves and wave structures for a rail are detailed. The rail is used to represent a bar with a complex cross-section. Similar to the square rod results, the rail's longitudinal modes have higher group velocities. However, the rail contains modes with local vibration. Finally, single mode detection and excitation techniques are introduced. A single mode can be obtained by detecting and exciting with a weighted function that corresponds to a specific mode's wave structure.     

Comparison between the dispersion curves calculated in complex frequency and the minima of the reflection coefficients for an embedded layer

Analytical solutions of Lamb functions for symmetric and antisymmetric elastodynamic modes propagating within a solid layer embedded in an infinite medium are presented. Alternative theoretical analyses of such modes are performed, first in terms of the usual approach of harmonic heterogeneous plane waves (real frequency and complex slowness) and then in terms of transient homogeneous plane waves (complex frequency and real slowness). An example structure of a 0.1-mm-thick "alpha case" (an oxygen-rich phase of titanium that is relatively stiff) plate embedded in titanium is used for the study. A large difference between the usual dispersion curves calculated in real frequency and complex slowness and those calculated in complex frequency and real slowness is shown. Thus the choice between a spatial and a temporal parameter to describe the imaginary part of the guided waves is shown to be significant. The minima and the zeros of the longitudinal and shear plane-wave reflection coefficients are calculated and are compared with the dispersion curves. It is found that they do not match with the dispersion curves for complex slowness, but they do agree quite well with the dispersion curves for complex frequency. This implies that the complex frequency approach is better suited for the comparison of the modal properties with near-field reflection measurements.     

Piezoacoustic wave spectra using improved surface impedance matrix: application to high impedance-contrast layered plates

Starting from the general modal solutions for a homogeneous layer of arbitrary material and crystalline symmetry, a matrix formalism is developed to establish the semianalytical expressions of the surface impedance matrices (SIM) for a single piezoelectric layer. By applying the electrical boundary conditions, the layer impedance matrix is reduced to a unified elastic form whether the material is piezoelectric or not. The characteristic equation for the dispersion curves is derived in both forms of a three-dimensional acoustic SIM and of an electrical scalar function. The same approach is extended to multilayered structures such as a piezoelectric layer sandwiched in between two metallic electrodes, a Bragg coupler, and a semi-infinite substrate as well. The effectiveness of the approach is numerically demonstrated by its ability to determine the full spectra of guided modes, even at extremely high frequencies, in layered plates comprising up to four layers and three materials. Negative slope in f-k curve for some modes, asymptotic behavior at short wavelength regime, as well as wave confinement phenomena made evident by the numerical results are analyzed and interpreted in terms of the surface acoustic waves and of the interfacial waves in connection with the bulk waves in massive materials.