Separation of equation of motion of love waves in curvilinearly anisotropic inhomogeneous medium

Summary The general problem of Love wave propagation, in a medium with cylindrical anisotropy of hexagonal type, is formulated. The method of seperation of variables is applied to examine the possibility of obtaining formal solutions for different types of inhomogeneities present in the medium. It is found that when the elastic parameters (C ₄₄ and ) are functions of bothv and the equation of motion is not separable. The use of the technique is illustrated, by considering radial inhomogeneity in an anisotropic cylindrical crustal layer, for obtaining the characteristic frequency equation of Love waves in such a medium.     

Reflection and transmission of plane waves by a layer of compact inhomogeneities

Reflection and transmission of elastic wave motion by a layer of compact inhomogeneities has been analyzed. For identical inhomogeneities whose geometrical centers are periodically spaced, the problem has been formulated and solved rigorously. The reflected and transmitted longitudinal and transverse wave motions have been expressed as superpositions of wavemodes, where each wavemode has its own cut-off frequency. At its cut-off frequency a mode converts from a standing into a propagating wavemode. The standing wavemodes decay exponentially with distance to the plane of the centers of the inhomogeneities. At small frequencies only the lowest order modes of longitudinal and transverse wave motion are propagating. Reflection and transmission coefficients have been defined in terms of the coefficients of the zeroth-order scattered wavemodes. These coefficients have been computed by a novel application of the Betti-Rayleigh reciprocal theorem. They are expressed as integrals over the surface of a single inhomogeneity, in terms of the displacements and tractions on the surface of the inhomogeneity. The system of singular integral equations for the surface fields has been solved numerically by the boundary integral equation method. Curves show the reflection and transmission coefficients for the reflected and transmitted longitudinal and transverse waves as functions of the frequency. Some results are also presented for planar distributions of cracks whose spacing and size are random variables. Finally, dispersion relations are discussed for solids which are completely filled with periodically spaced inhomogeneities.     

LOVE SEAM-WAVES IN A HORIZONTALLY INHOMOGENEOUS THREE-LAYERED MEDIUM*

Using the WKBJ-method the absorption-dispersion relation and the amplitude functions are derived for Love seam-waves that propagate in a horizontally inhomogeneous three-layered medium. To describe the anelastic friction the constant Q-model is applied. The inhomogeneity that appears in either the elastic moduli or quality factors is assumed to remain weak in the coal as well as in the adjacent layers, which are assumed to have different material properties (asymmetric channel). Using numerical solutions of the dispersion relation, it is shown that the weak horizontal inhomogeneities can be optimally detected using channel-wave constituents of a frequency near to the Airy frequency while inhomogeneities of the adjacent rock can only be detected at frequencies close to, but higher than, the cut-off frequency.     

ELASTIC EXTRAPOLATION OF PRIMARY SEISMIC P- AND S-WAVES1

The elastic Kirchhoff-Helmholtz integral expresses the components of the monochromatic displacement vector at any point A in terms of the displacement field and the stress field at any closed surface surrounding A. By introducing Green's functions for P- and S-waves, the elastic Kirchhoff-Helmholtz integral is modified such that it expresses either the P-wave or the S-wave at A in terms of the elastic wavefield at the closed surface. This modified elastic Kirchhoff-Helmholtz integral is transformed into one-way elastic Rayleigh-type integrals for forward extrapolation of downgoing and upgoing P- and S-waves. We also derive one-way elastic Rayleigh-type integrals for inverse extrapolation of downgoing and upgoing P- and S-waves. The one-way elastic extrapolation operators derived in this paper are the basis for a new prestack migration scheme for elastic data.     

Elastic wave scattering by inhomogeneous and anisotropic bodies in a half space

Scattering of elastic waves by inhomogeneous and anisotropic bodies in a half space is considered. The integral equation method is formulated by using the fundamental solution of a homogeneous isotropic body in elastostatics and regarding the resulting inhomogeneous terms as equivalent body forces. Numerical examples are presented for the wave scattering by inhomogeneous and/or anisotropic alluvial valleys and for the dynamic analysis of an inhomogeneous dam. The effect of inhomogeneities and anisotropy on the dynamic behaviour of alluvial valley and dam is discussed.     

Piezomagnetic fields produced by dislocation sources

Tectonomagnetic modeling based on the linear piezomagnetic effect is reviewed with special attention to dislocation models. Stacey's scheme was the prototype for such modeling, as proposed in his first seismomagnetic calculations in 1964. The linear piezomagnetic law is presented, in which the stress-induced magnetization is expressed as a linear combination of stress components. The Gauss law for magnetic field and the Cauchy-Navier equation for static elastic equilibrium are combined through linear piezomagnetism and the Hooke law to yield the basic equation for piezomagnetic potential. A representation theorem for its solution is given by surface integrals of the displacement and its normal derivative over the strained body.A Green's function method is developed to compute the piezomagnetic field produced by a dislocation surface in an elastic half-space. Volterra's formula for piezomagnetic potential is derived by modifying Stacey's scheme for tectonomagnetic modeling. The Green's functions for the problem are called elementary piezomagnetic potentials, which are defined as potentials produced by elementary dislocations. Special consideration is required to construct the elementary piezomagnetic potentials, because the stress field around a point dislocation has a singularity of orderr ^–3. The integral representing elementary piezomagnetic potentials is not uniformly convergent. Owing to inappropriate convergency, the Green's functions obtained in an earlier study led to a puzzling outcome. Revised Green's functions give consistent results with those obtained so far by numerical integrations. Generally the piezomagnetic field produced by dislocation sources is weak in the case of a homogeneous earth model. Two enhancement effects for piezomagnetic signals are suggested: one due to inhomogeneous magnetization and the other via bore-hole observations.     

Dynamics of rigid strip foundations embedded in orthotropic elastic soils

This study is concerned with the dynamic response of an arbitrary shaped rigid strip foundation embedded in an orthotropic elastic soil. The foundation is subjected to time-harmonic vertical, horizontal and moment loadings. The boundary-value problem related to an embedded foundation is analysed by using the indirect boundary integral equation method. The kernel functions of the integral equations are displacement and traction Green's functions of an anisotropic elastic half plane. Exact analytical solutions are used for the Green's functions. The boundary integral equation is solved by using numerical techniques. Selected numerical results are presented for the impedances of rectangular and semi-circular rigid strip foundations embedded in four types of anisotropic soils. A discussion on the influence of soil anisotropy and frequency of excitation on the impedances is presented. The versatility of the analysis is demonstrated by considering the through soil interaction between two semi-circular strip foundations.     

两层大地中三维体的激发极化与电阻率响应的积分方程模拟

用散射、叠加方法推导出两层大地的并矢格林函数。使用这些函数,含三维异常体的二层大地的边值问题转变为积分方程。使用矩量法,可解此积分方程。 使用偶极-偶极装置计算了激发极化和电阻率响应的几个数值结果,并对这一算法作了几方面检验。     

Torsional vibration of an orthotropic Cylindrical Shell

Summary This paper is concerned with the determination of torsional vibration of an orthotropic Cylindrical Shell assuming the elastic constants to be proportional to thenth power of the distance from the axis of the shell. The solution is obtained in terms of Bessel functions.c/o Dr. A. K. Dutta, P. 7. A H. B. Town, Block A, P. O. Sodepur, Dt. 24-Paraganas, West Bengal, India.     

Propagation of Rayleigh waves in an axially symmetric heterogeneous layer lying between two homogeneous halfspaces

Summary The propagation of Rayleigh type waves in an axially symmetric inhomogeneous layer lying between two halfspaces is studied. The halfspaces are supposed to be identical in their elastic properties. The variation of the parameters in the layer is assumed to be of the form where is a constant andz is the distance measured from one interface into the layer. With this assumption, the vector wave equation for the layer is separable. The solution is obtained in terms ofWhittaker's functions and the frequency equation of Rayleigh type waves is derived.     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.     

Nature of the scattered seismic response from zones of random clusters of cavities and fractures in a massive rock

The scattering of elastic energy by random clusters of fractures and/or cavities in a massive rock is studied. The interpretation of the scattered seismic response reveals crucial information about the clusters of inhomogeneities (fractures/cavities), which may correspond to reservoirs. The study is based on a new two‐dimensional numerical‐modelling method that relaxes the constraints on the location and orientation of the inhomogeneities, accounts for inhomogeneities that have almost no volume but a finite surface area (fractures) and improves the accuracy of the calculation when the size of the inhomogeneities is comparable to the mesh size. It is shown that the nature of the seismic response of zones of diffuse fracturing and/or cavities is associated with the non‐uniformity of micro‐inhomogeneities in such zones; accumulations of these micro‐inhomogeneities are known as clusters. The relationship between the non‐uniformity of micro‐inhomogeneities and the strength of the seismic response has been established and measured. Considerable differences in the structure of the seismic response of zones of diffuse fracturing and diffuse cavities have been identified. Converted PS‐waves dominate in the scattered wavefield associated with fractures. This is explained, as the modelling results show, by a greater transparency of fluid‐filled fractures, which reduces the reflected energy of compressional waves. The wavefield associated with cavities is characterized by the predominance (in terms of strength) of compressional PP‐waves. The strength of converted PS‐waves in the scattered wavefields for both media is approximately the same. On the whole, according to the results of the modelling, the energy of the scattered response of fractured reservoirs is considerably less (about two times) than that of cavernous reservoirs.     

Rayleigh wave dispersion equation for a layered spherical earth with exponential function solutions in each shell

We consider the second-order differential equations ofP-SV motion in an isotropic elastic medium with spherical coordinates. We assume that in the medium Lamé's parameters , r ^p and compressional and shear-wave velocities , r, wherer is radial distance. With this regular heterogeneity both the radial functions appearing in displacement components satisfy a fourth-order differential equation which provides solutions in terms of exponential functions. We then consider a layered spherical earth in which each layer has heterogeneity as specified above. The dispersion equation of the Rayleigh wave is obtained using the Thomson-Haskel method. Due to exponential function solutions in each layer, the dispersion equation has similar simplicity, as in a flat-layered earth. The dispersion equation is further simplified, whenp=–2. We obtain numerical results which agree with results obtained by other methods.     

Scattering of elastic waves by a general anisotropic basin. Part 2: a 3D model

Scattering of plane harmonic waves by a three‐dimensional basin of arbitrary shape embedded within elastic half‐space is investigated by using an indirect boundary integral equation approach. The materials of the basin and the half‐space are assumed to be the most general anisotropic, homogeneous, linearly elastic solids without any material symmetry (i.e. triclinic). The unknown scattered waves are expressed in terms of three‐dimensional triclinic time harmonic full‐space Green's functions. The results have been tested by comparing the surface response of semi spherical isotropic and transversely isotropic basins for which the numerical solutions are available. Surface displacements are presented for a semicircular basin subjected to a vertical incident plane harmonic pseudo‐P‐, SV‐, or SH‐wave. These results are compared with the motion obtained for the corresponding equivalent isotropic models. The results show that presence of the basin may cause significant amplification of ground motion when compared to the free‐field displacements. The peak amplitude of the predominant component of surface motion is smaller for the anisotropic basin than for the corresponding isotropic one. Anisotropic response may be asymmetric even for symmetric geometry and incidence. Anisotropic surface displacement generally includes all three components of motion which may not be the case for the isotropic results. Furthermore, anisotropic response strongly depends upon the nature of the incident wave, degree of material anisotropy and the azimuthal orientation of the observation station. These results clearly demonstrate the importance of anisotropy in amplification of surface ground motion. Copyright © 2003 John Wiley & Sons, Ltd.     

Thermo-elastic interactions in an infinite elastic solid due to a concentrated transient heat source

Summary The present paper is concerned with the determination of thermo-elastic stress and temperature distribution in an infinite elastic solid when it is subjected to a concentrated transient heat source taking into account the effect of coupling. The fundamental partial differential equation for the thermo-elastic potential is solved by means of the operational method. The solutions are obtained for small values of time in terms of known functions such as complementary error functions and associated complementary error functions and have the form of power series with respect to the parameter of coupling . Since this parameter is generally very small for most of the metals, the expressions for stress and temperature distribution have been obtained by considering terms upto linear .     

A Second-Order Approach for P-Wave Elastic Impedance Technology

Previous formulation for P-wave elastic impedance (EI) technology considers only first-order effects in isotropic reflectivity. In this paper, Wang's pseudo-quadratic approximation for PP-wave reflection (R_PP) coefficients is used, in order to incorporate nonlinear effects into EI equation. In comparison with coefficients computed with the conventional linear approximation, Wang's pseudo-quadratic formula shows higher accuracy at far incidence. A further nonlinear component in the intermediate region of incidence is responsible for the high accuracy achieved with the pseudo-quadratic Rpp coefficient formula. By applying the same procedures of previous linear formulation to Wang's pseudo-quadratic Rpp coefficient, a second-order approach for EI equation is obtained. This novel approach is formed by multiplication of two terms. The first term represents the previous linear approach for EI equation. As for the second term, it is interpreted as the correction of first-order EI formula to second-order effects. As expected, specialization of the second-order EI equation to normal incidence results in the well-known acoustic impedance (AI). Assumption of invariability in fundamental elastic properties leads to simplification of mathematical procedures. However, high contrasts possibly found within the log window under investigation may corrupt the computation of EI logs by introducing numerical errors. Although two procedures are proposed to cope with numerical errors, modeling shows that the second-order approach for EI is robust enough to handle high contrasts in elastic parameters. Actual well logs are used to verify performance of the novel EI equation in reproducing the amplitude-versus-offset (AVO) response of a mature, oil-bearing sandstone resevoir. As a result, influence of nonlinear effects, which are incorporated into EI equation, is observed on amplitudes and on frequency bandwidth of synthetic seismograms generated at a high angle of incidence. Further experiments with actual well data focus on crossplotting EI logs against fundamental elastic parameters. In terms of accuracy, the outcomes reveal that lithofacies classification can benefit from using the elaboration of EI technology derived in this work.     

Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements

The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500 km², as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.We derive a set of local Green's functions defined globally on a 1° × 1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0 mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.     

Analytical representation of Green’s functions of the Biot equations for elastic waves in a homogeneous two-phase medium

The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride’s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green’s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green’s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green’s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type.

     

An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P,S and Rayleigh waves

A boundary integral formulation is presented and applied to model the ground motion on alluvial valleys under incident P, S and Rayleigh waves. It is based on integral representations for the diffracted and the refracted elastic waves using single-layer boundary sources. This approach is called indirect BEM in the literature as the sources' strengths should be obtained as an intermediate step. Boundary conditions lead to a system of integral equations for boundary sources. A discretization scheme based on the numerical and analytical integration of exact Green's functions for displacements and tractions is used. Various examples are given for two-dimensional problems of diffraction of elastic waves by soft elastic inclusion models of alluvial deposits in an elastic half-space. Results are displayed in both frequency and time domains. These results show the significant influence of locally generated surface waves in seismic response and suggest approximations of practical interest. For shallow alluvial valleys the response and its resonant frequencies are controlled by a coupling mechanism that involves both the simple one-dimensional shear beam model and the propagation of surface waves.     

Using the pseudospectral method on curved grids for 2D elastic forward modelling1

When applying the conventional Fourier pseudospectral method (FSM) on a Cartesian grid that has a sufficient size to propagate a pulse, spurious diffractions from the staircase representation of the curved interfaces appear in the wavefield. It is demonstrated that these non-physical diffractions can be eliminated by using curved grids that conform to all the interfaces of the subsurface. Methods for solving the 2D acoustic wave equation using such curved grids have been published previously by the authors. Here the extensions to the full 2D elastic wave equations are presented. The curved grids are generated by using the so-called multiblock strategy which is a well-known concept in computational fluid dynamics. In principle the sub-surface is divided into a number of contiguous subdomains. A separate grid is generated for each subdomain patching the grid lines across domain boundaries to obtain a globally continuous grid. Using this approach, even configurations with pinch outs can be handled. The curved grid is taken to constitute a generalized curvilinear coordinate system. Thus, the elastic equations have to be written in a curvilinear frame before applying the numerical scheme. The method implies that twice the number of spatial derivatives have to be evaluated compared to the conventional FSM on a Cartesian grid. However, it is demonstrated that the extra terms are more than compensated for by the fewer grid points needed in the curved approach.