Indirect Boundary Element Method applied to fluid–solid interfaces

In this paper scattering of elastic waves in fluid–solid interfaces is investigated. We use the Indirect Boundary Element Method to study this wave propagation phenomenon in 2D models. Three models are analyzed: a first one with an interface between two half-spaces, one fluid on the top part and the other solid in the bottom; a second model including a fluid half-space above a layered solid; and finally, a third model with a fluid layer bounded by two solid half-spaces. The source, represented by Hankel's function of the second kind, is always applied in the fluid. This indirect formulation can give to the analyst a deep physical insight on the generated diffracted waves because it is closer to the physical reality and can be regarded as a realization of Huygens' principle. In any event, mathematically it is fully equivalent to the classical Somigliana's representation theorem. In order to gauge accuracy we test our method by comparing with an analytical solution known as Discrete Wave Number. A near interface pulse generates scattered waves that can be registered by receivers located in the fluid and it is possible to infer wave velocities of solids. Results are presented in both time and frequency domain, where several aspects related to the different wave types that emerge from this kind of problems are pointed out.     

Inversion of Scholte wave dispersion and waveform modeling for shallow structure of the Ninetyeast Ridge

The construction of S-wave velocity models of marine sediments down to hundreds of meters below the seafloor is important in a number of disciplines. One of the most significant trends in marine geophysics is to use interface waves to estimate shallow shear velocities which play an important role in determining the shallow crustal structure. In marine settings, the waves trapped near the fluid–solid interface are called Scholte waves, and this is the subject of the study. In 1998, there were experiments on the Ninetyeast Ridge (Central Indian Ocean) to study the shallow seismic structure at the drilled site. The data were acquired by both ocean bottom seismometer and ocean bottom hydrophone. A new type of seafloor implosion sources has been used in this experiment, which successfully excited fast and high frequency (>500 Hz) body waves and slow, intermediate frequency (<20 Hz) Scholte waves. The fundamental and first higher mode Scholte waves have both been excited by the implosion source. Here, the Scholte waves are investigated with a full waveform modeling and a group velocity inversion approach. Shear wave velocities for the uppermost layers of the region are inferred and results from the different methods are compared. We find that the full waveform modeling is important to understand the intrinsic attenuation of the Scholte waves between 1 and 20 Hz. The modeling shows that the S-wave velocity varies from 195 to 350 m/s in the first 16 m of the uppermost layer. Depths levels of high S-wave impedance contrasts compare well to the layer depth derived from a P-wave analysis as well as from drilling data. As expected, the P- to S-wave velocity ratio is very high in the uppermost 16 m of the seafloor and the Poisson ratio is nearly 0.5. Depth levels of high S-wave impedance contrasts are comparable to the layer depth derived from drilling data.     

软地层横波速度的偶极随钻测井反演

在电缆测井中,可以利用偶极声场中的弯曲波反演软地层的横波速度.然而,在随钻声波测井（LWD）中,钢制钻铤的存在使得井孔结构变得复杂,同时改变了井孔声场,弯曲波也变得难以测定.此外,弯曲波与钻铤波耦合在一起,使得地层横波速度的反演变得困难.本文计算分析了随钻声场的频散曲线和激发曲线,注意到了偶极舒尔特波在较宽频带内速度频散很弱,特别是在本文研究的软地层情况下,偶极舒尔特波速度在3至25 kHz的频率范围内几乎为一个常值,并且该值与地层横波速度存在一一对应关系.舒尔特波速度远小于其他模式波速度,与其他模式波在时域上易分离.相对于其他地层参数,舒尔特波对地层横波速度十分敏感.因此,它可以用来反演软地层的横波速度.

     

Scattering and diffraction of earthquake motions in irregular,elastic layers,II: Rayleigh and body P and SV waves

We present our study of the wave propagation in an irregularly layered, elastic wave-guide excited by incoming Rayleigh surface waves and P and SV body waves. Our aim is to show examples of applying a method that will make it possible to analyze the distribution and amplification of displacements, rotations, curvatures, strains, and stresses on or below the ground surface during passage of strong earthquake ground motion. We employ the weighted-residuals method, which makes it possible to calculate the scattered and diffracted waves, and then we illustrate the amplification of motions in the vicinity of inhomogeneity.     

Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method

In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion.     

Diffraction of seismic waves in an elastic, cracked halfplane using a boundary integral formulation

The presence of subsurface cracks in a halfspace excited by elastic waves may give rise to scattered body and surface waves. For many engineering applications, such as non-destructive testing or oil exploration, the scattered field may yield valuable information to detect cracks and other scatterers. We use the Indirect Boundary Element Method (IBEM) to study the diffraction of P, SV waves with various incidence angles and Rayleigh surface waves. This approximate boundary integral technique is based upon the integral representation for scattered elastic waves using single-layer boundary sources. Our approach is usually called indirect BEM as the sources' strengths should be obtained as an intermediate step. This indirect formulation can give to the analyst a deep physical insight on the generated diffracted waves because it is closer to the physical reality and can be regarded as a realization of Huygens' Principle. In any event, mathematically it is fully equivalent to the classical Somigliana's representation theorem. In order to gauge accuracy we test our method by comparing with previous results in the literature. Various crack configurations, including multiple cracks, are investigated. Results in frequency and time domains are displayed. Under certain conditions the amplitude spectra of those waves clearly show conspicuous resonance peaks.     

3-D wave propagation in fluid-filled irregular boreholes in elastic formations

Different seismic testing techniques rely on the propagation of acoustic waves in fluid-filled boreholes from sources placed within the borehole and in the solid media. The interpretation of the signals recorded relies on understanding how waves propagate in the borehole and its immediate vicinity. It is known that very complex wave patterns can arise, depending on the distance between the source and the receiver, and their placement and orientation relative to the axis of a circular borehole. The problem becomes more complex if the cross-section is not circular, conditions for which analytical solutions are not known. In this work, the Boundary Element Method (BEM) is used to evaluate the three-dimensional wave field elicited by monopole sources in the vicinity of a fluid-filled borehole. This model is used to assess the effects of the receiver position on the propagation of both axisymmetric and non-axisymmetric wave modes when different borehole cross-sections are used. Both frequency vs. axial-wave number responses and time-domain responses are calculated.     

Oblique wave diffraction by a flexible floating structure in the presence of a submerged flexible structure

Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering.     

Surface waves in a micropolar elastic solid containing voids

The present paper investigates the effect of voids on the propagation of surface waves in a homogeneous micropolar elastic solid medium which contains a distribution of vacuous pores (voids). The general theory for surface wave propagation in micropolar elastic media containing voids has been presented. Particular cases of surface waves (Rayleigh’s, Love’s and Stoneley’s) in micropolar media which contain vacuous pores have been deduced from the above general theory. Discussions have been made in each case to highlight the effect of voids and micropolar character of the material medium separately. Their joint effect has also been studied in details. Modulation of Rayleigh wave velocity has been studied numerically. It is observed that Love waves are not affected by the presence of voids.     

Propagation of guided waves in a fluid layer bounded by two viscoelastic transversely isotropic solids

The analysis of Stoneley wave propagation in a fracture is essential for the identification and evaluation of fracture parameters from the borehole Stoneley wave. Also, it is important for many geophysics considerations, e.g. for tremor and long-period events observed in volcanoes and geothermal areas. In this paper, we investigate the guided waves propagation in a fluid layer lying between two viscoelastic vertically transversely isotropic media. The viscoelastic mechanism models the attenuation due to the presence of fluid saturation in the rock. A model based on the superposition of three inhomogeneous partial plane waves: one in the fluid and two heterogeneous waves in the solid is developed. The dispersion equation is obtained for this case. A numerical solution is carried out to obtain the guided wave velocity and attenuation coefficient. The results of this investigation show that there is a strong correlation between the velocity dispersion and attenuation of Stoneley wave and the anisotropic parameters of the medium especially in a sandstone (fast) medium.     

瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明： 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响.     

交错网格高阶差分法三维弹性波数值模拟

本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。     

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.     

三维层状孔隙介质中弹性波的一种积分表达式I：理论

孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.     

饱和多孔介质近场波动分析的一种黏弹性人工边界

采用有限模型数值方法求解能量开放系统中的波动问题时,虚拟人工边界的处理方式对计算结果的准确性和精度具有重要的影响.本文针对无限域饱和多孔介质中波传播问题的人工边界处理方式进行了研究,提出了饱和多孔介质近场波动分析的一种黏弹性人工边界处理方法.在考虑多孔介质中固相和液相的相互作用的情况下,通过在人工边界处分别施加反映固相和液相介质波传播效应的弹簧及阻尼来模拟饱和多孔介质中波的能量辐射效应影响.算例表明,本文建议的黏弹性人工边界具有较好的模拟效果.     

Rectangular footing on soil with depth-degrading stiffness: Vertical and rocking impedances under conditional existence of surface waves

The response of rectangular rigid footings resting on an elastic soil of shear modulus decreasing monotonically with depth is studied. Such profiles are typically encountered after ground improvement. The propagation characteristics of SV/P surface waves are investigated, showing the appearance of cut-off frequencies above which surface waves do not exist. The semi-analytical method of the subdivision of the footing/soil contact area is then used for solving the boundary value problem, whereby the influence functions for the sub-regions are determined by integration of the corresponding surface-to-surface Green׳s functions. Impedance functions are presented over a wide range of frequencies for typical values of the non-homogeneity parameters, the Poisson׳s ratio and the foundation geometry. The salient features that are associated with the non-homogeneity and the appearance of cut-off frequencies are elucidated.     

On the modified Taylor constraint

A dynamo driven by motions unaffected by viscous forces is termed magnetostrophic. Although such a model might describe magnetic field generation in Earth’s core well, a magnetostrophic dynamo has not yet been found even though Taylor [Proc. R. Soc. Lond. A 1963, 274, 274–283] devised an apparently viable method of finding one. His method for determining the fluid velocity from the magnetic field and the energy source involved only the evaluation of integrals along lines parallel to the Earth’s axis of rotation and the solution of a second-order ordinary differential equation. It is demonstrated below that an approximate solution of this equation for a broad family of magnetic fields is immediate. Furthermore inertia, which was neglected in Taylor’s theory, is restored here, so that the modified theory includes torsional waves, whose existence in the Earth’s core has been inferred from observations of the length of day. Their theory is reconsidered.     

On the completeness of inertial wave modes in rotating annular channels

An important unanswered mathematical question in the theory of rotating fluids has been the completeness of the inviscid eigenfunctions which are usually referred to as inertial waves or inertial modes. We provide for the first time a mathematical proof for the completeness of the inertial modes in a rotating annular channel by establishing the completeness relation, or Parseval’s equality, for any piecewise continuous, differentiable velocity of an incompressible fluid.     

各向异性双相介质弹性波场褶积算法数值模拟

从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到.     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函...