Time domain scattering of acoustic plane waves by vertical faults

Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large.     

Scattering of Seismic Waves by Cracks with the Boundary Integral Method

— We develop a new scheme to compute 2-D SH seismograms for media with many flat cracks, based on the boundary integral method. A dry or traction-free boundary condition is applied to crack surfaces although other kinds of cracks such as wet or fluid-saturated cracks can be treated simply by assigning different boundary conditions. While body forces are distributed for cavities or inclusions to express scattered wave, dislocations (or displacement discontinuities between the top and the bottom surfaces of each crack) are used as fictitious sources along crack surfaces. With these dislocations as unknown coefficients, the scattered wave is expressed by the normal derivative of Green's function along the crack surface, which is called “double-layer potentials” in the boundary integral method, while we used “single-layer potentials” for cavities or inclusions. These unknowns are determined so that boundary conditions or crack surfaces are satisfied in the least-squared sense, for example, traction-free for dry cracks. Seismograms with plane-wave incidence are synthesized for homogeneous media with many cracks. First, we check the accuracy of our scheme for a medium with one long crack. All the predicted phases such as reflected wave, diffraction from a crack tip and shadow behind the crack are simulated quite accurately, under the same criterion as in the case for cavities or inclusions. Next, we compute seismograms for 50 randomly distributed cracks and compare them with those for circular cavities. When cracks are randomly oriented, waveforms and the strength of scattering attenuation are similar to the cavity case in a frequency range higher than k d $\simeq$ 2 where the size of scatterers d (i.e., crack length or cavity diameter) is comparable with the wavelength considered (k is the wavenumber). On the other hand, the scattering attenuation for cracks becomes much smaller in a lower frequency range (k d<2) because only the volume but not detail geometry of scatterers becomes important with wavelength much longer than each scatterer. When all the cracks are oriented in a fixed direction, the scattering attenuation depends strongly on the incident angle to the crack surface as frequency increases (k d>2): scattering becomes weak for cracks oriented parallel to the direction of the incident wave, while it gets close to the cavity case for cracks aligned perpendicular to the incident wave.     

Scattering ofP andS waves by a spherically symmetric inclusion

Scattering of an arbitrary elastic wave incident upon a spherically symmetric inclusion is considered and solutions are developed in terms of the spherical vector system of Petrashen, which produces results in terms of displacements rather than displacement potentials and in a form suitable for accurate numerical computations. Analytical expressions for canonical scattering coefficients are obtained for both the cases of incidentP waves and incidentS waves. Calculations of energy flux in the scattered waves lead to elastic optical theorems for bothP andS waves, which relate the scattering cross sections to the amplitude of the scattered fields in the forward direction. The properties of the solutions for a homogeneous elastic sphere, a sphere filled by fluid, and a spherical cavity are illustrated with scattering cross sections that demonstrate important differences between these types of obstacles. A general result is that the frequency dependence of the scattering is defined by the wavelength of the scattered wave rather than the wavelength of the incident wave. This is consistent with the finding that the intensity of thePS scattering is generally much stronger than theSP scattering. When averaged over all scattering angles, the mean intensity of thePS converted waves is2V _p ² /V _s ⁴ times the mean intensity of theSP converted waves, and this ratio is independent of frequency. The exact solutions reduce to simple and easily used expressions in the case of the low frequency (Rayleigh) approximation and the low contrast (Rayleigh-Born) approximation. The case of energy absorbing inclusions can also be obtained by assigning complex values to the elastic parameters, which leads to the result that an increase in attenuation within the inclusion causes an increased scattering cross section with a marked preference for scatteredS waves. The complete generality of the results is demonstrated by showing waves scattered by the earth's core in the time domain, an example of high-frequency scattering that reveals a very complex relationship between geometrical arrivals and diffracted waves.     

Scattering of SH-wave by multiple circular cavities in half space

In this paper, an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space. The scattered wave function used for scattering of SH-waves by multiple circular cavities, which automatically satisfies the stress-free condition at the horizontal surface, is constructed by applying the symmetry of the SH-wave scattering and the method of multi-polar coordinates system. Applying this scattered wave function and method of moving coordinates, the original problem can be transformed to the problem of SH-wave scattering by multiple circular cavities in the full space. Finally, the solution of the problem can be reduced to a series of algebraic equations and solved numerically by truncating the infinite algebraic equations to the finite ones. Numerical examples are provided for case with two cavities to show the effect of wave number, and the distances between the centers of the cavities and from the centers to the ground surface on the dynamic stress concentration around the cavity impacted by incident steady SH-wave.     

Seismic energy partitioning and scattering in laterally heterogeneous ocean crust

We present finite difference forward models of elastic wave propagation through laterally heterogeneous upper oceanic crust. The finite difference formulation is a 2-D solution to the elastic wave equation for heterogeneous media and implicitly calculatesP andSV propagation, compressional to shear conversion, interference effects and interface phenomena. Random velocity perturbations with Gaussian and self-similar autocorrelation functions and different correlation lengths (a) are presented which show different characteristics of secondary scattering. Heterogeneities scatter primary energy into secondary body waves and secondary Stoneley waves along the water-solid interface. The presence of a water-solid interface in the model allows for the existence of secondary Stoneley waves which account for much of the seafloor noise seen in the synthetic seismograms for the laterally heterogeneous models.Random incoherent secondary scattering generally increases aska (wavenumber,k, and correlation length,a) approaches one. Deterministic secondary scattering from larger heterogeneities is the dominant effect in the models aska increases above one. Secondary scattering also shows up as incoherence in the primary traces of the seisograms when compared to the laterally homogeneous case. Cross-correlation analysis of the initialP-diving wave arrival shows that, in general, the correlation between traces decreases aska approaches one. Also, because many different wave types exist for these marine models, the correlation between traces is range dependent, even for the laterally homogeneous case.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Harmonic wave diffraction by two circular cavities in a poroelastic formation

This study investigates the dynamic interaction of time harmonic plane waves with a pair of parallel circular cylindrical cavities of infinite length buried in a boundless porous elastic fluid-saturated medium. The novel features of Biot dynamic theory of poroelasticity along with the appropriate wave field expansions, the pertinent boundary conditions, and the translational addition theorems for cylindrical wave functions are employed to develop a closed-form solution in the form of infinite series. The analytical results are illustrated with numerical examples in which two empty cavities are insonified by a fast compressional or a shear wave at end-on incidence. The basic dynamic field quantities such as the hoop stress amplitude and the radial displacement of the elastic frame are evaluated and discussed for representative values of the parameters characterizing the system. The effects of the proximity of the two cavities, the incident wave frequency and type are examined. Particular attention has been focused on multiple scattering interactions in addition to the slow wave coupling effects which is known to be the primary distinction of the scattering phenomenon in poroelasticity from the classical elastic case. Limiting case involving two empty cylindrical cavities in an elastic solid is considered and excellent agreement with a well-known solution is established.     

Model studies on electro-telluric interpretation problems of some geological inhomogeneities

Summary The rotating nature of telluric field is simulated in an electrolytic model tank. Variations of this rotating field due to the presence of a few two and three dimensional geological models, simulated in the tank, are studied. The telluric ellipses, recorded on a C.R.O. Screen, are found to be useful for qualitative and quantitative interpretation of the non-linearly polarised telluric field data. Relative usefulnesses of the various parameters, e.g., (i)M=ratio of the semimajor and semiminor axis, (ii)K = square root of normalized surface areas, (iii) = linear eccentricity ((a ²/b ²-1) and (iv)e = conventional eccentricity ((1-b ²/a ²)), of a telluric ellipse are tested. Linear eccentricity appears to be the most sensitive parameter. Methods for determining depth and horizontal extent of a structure from a telluric map are suggested. Near surface inhomogeneities, as revealed from the experiment are unlikely to affect the telluric map due to a basement structure.     

Diffraction of SV waves by underground, circular, cylindrical cavities

The scattering and diffraction of plane SV waves underground, circular, cylindrical cavities at various depths in an elastic half space is studied in this paper. The cavities, studied here, are at depths of two to five cavity radii, measured from the surface to the center of the cavity. Fourier-Bessel series are used to satisfy the wave equation and the boundary conditions. When the angle of incidence of the plane SV wave exceeds the critical angle, surface waves are generated, which are expanded in terms of Fourier series, which also involve Bessel functions. The surface displacement amplitudes and phases that are presented show that the results depend on the following parameters: (1) The angle of incidence, θ_β; (2) the ratio cavity depth to the cavity radius, h/a; (3) the dimensionless frequency of the incident SV wave, η; and (4) Poisson's ratio, v. The presence of the cavity in the half space results in significant deviation of both the displacement amplitudes and phases on the nearby half space surface from that of a uniform half space.     

The secular variation of the geomagnetic field in Egypt in the last 5000 years

The palaeo-intensities (F _a) of the geomagnetic field in Egypt at some ages are determined by archaeomagnetic measurements and found to be:F _a=36.2 T at 3100 B.C., F_a=46.8 T at 3000 B.C.,F _a=36.5 T at 2780 B.C., 49.0 T at 2500 B.C., 36.4 T at 2200 B.C., 57.5 T at 1990 B.C., 62.1 T atca 1400 B.C., 61.5 T at 1400 B.C., 69.9 T at 600 B.C., 59.3 T at 550 B.C., 79.9 T at 460 B.C., 73.7 T at 450 B.C., 69.7 T at 320 B.C., 56.2 T at A.D. 50, 64.9 T, at A.D. 400, 54.4 T at A.D. 300, 57.5 T at A.D. 700 and 43.0 T at A.D. 1975.The palaeo-inclinations (I _a) at some ages are found to be:I _a=24.2° at 420 B.C., 44° at A.D. 50, 60.7° at A.D. 703 and 42° at A.D. 1795.The measured values ofF _a are affected by the anisotropy of magnetic susceptibility of the samples by 13% to 20% of the expected correct value. The suitable correction of this effect is by multiplyingF by 1/((1+0.2(/90)) andF by 1/((1–0.13 (/90)), whereF andF are the resultant values ofF _a if the laboratory field is perpendicular or parallel to the wall of the sample during the Thelliers' experiments, respectively, and is the angle between the direction of natural remnant magnetization of the sample and the direction of the laboratory field.The results of this paper, together with the previous results for Egypt and the neighbourhoods, lead to the production of the secular variation curve of the geomagnetic field in Egypt for the last 5000 years. The intensity of the field shows a periodicity of about 400 years with multiples.     

Laboratory measurements of spatial fluctuation and attenuation of elastic waves by scattering due to random heterogeneities

To study the effects of strong scattering on elastic waves, spatial fluctuation and scattering attenuation ofP waves were examined by laboratory experiments for 2-D models of random media approximately characterized by a triangular correlation function in the range of 2<ka<33, wherek is the wave number anda is the correlation distance of the heterogeneities, i.e., the heterogeneity size. The results obtained are as follows: (1) Forka>10, both the intensity and the correlation distance of the amplitude fluctuation are approximate for any phase of theP-wave train. The correlation distance nearly agrees with the heterogeneity size. These fluctuation properties are quite consistent with the theoretical prediction by the forward-scattering approximation. (2) For 3<ka<6, the fluctuation intensity becomes stronger in later phases of theP-wave train. This shows that scattering is approximately isotropic, and therefore, the scattered energy increases with time within theP-wave train. The correlation distance of the amplitude fluctuation disagrees with the heterogeneity size, and it shows a frequency-dependent property decreasing from 7a to 4a with the increase ofka from 3 to 6. These properties for 3<ka<6 have not yet been predicted theoretically. (3) Forka<3, though the fluctuation is considerably smaller compared with that ofka>10 and 3<ka<6, the fluctuation property is considered similar to that of 3<ka<6. (4) The observed scattering attenuation,Q ^–1, increases withka forka<3, has a peak aroundka=35, and then decreases withka. (5) When _min = 15° and = 0.075, the theoreticalQ ^–1 curve, predicted by the approximate theory of Wu, roughly matches the observedQ ^–1 values, where _min is the minimum scattering angle measured from the propagation direction of theP waves and is the rms of fractional velocity fluctuation. This suggests that the energy scattered in the range of >15° is lost from theP waves, while the energy scattered in the range of <15° is retained; and that the approximate theory overestimates by about three times the value of the model media used owing to the neglect of multiple scattering. (6) When the size of velocity heterogeneities responsible for forward scattering at 3<ka<6 is estimated from the _min value of 15° on the basis of Wu's theory, it nearly agrees with the correlation distance for the initial phase of theP-wave train.     

An approximate method for soil-structure interaction for sh-waves—the born approximation

An approximate method is proposed for the scattering of SH-waves by foundations of irregular shape and the resulting soil-structure interaction problems. The scattering of elastic waves by the rigid foundation embedded in half-space is solved approximately by using integral representation of the wave equation. The procedure is the Born approximation which has been widely used in quantum mechanics for collision and scattering theory though not well-known in elastodynamics. This paper extends the previous work of the authors on the scattering of waves to account for soil-structure interaction. The motion of the foundation is evaluated by the balance of momentum under stresses due to the incident waves as well as the waves generated by its own motion and the forces coming from the superstructure. The model investigated consists of an infinitely long elastic shear wall of height H and thickness h erected on a rigid infinitely long foundation. Results are presented for the cases with circular, elliptical and rectangular foundations. For a circular foundation, excellent agreement is found with the exact solutions for the foundation displacement and the relative displacement between the top and bottom of the structure for the entire range of wave numbers. For an elliptical foundation, accuracy decreases with increasing wave numbers. Foundation displacements are compared for foundation shapes that are shallow elliptical, deep elliptical, rectangular and circular. It is observed that foundation displacements are dependent on the angle of incidence except for a semi-circle. The results on the details of the scattered field are, however, not as accurate.     

Lateral variations of coda Q and attenuation of seismic waves in Southwest Anatolia

The seismic quality factor (Q _c) and the attenuation coefficient (δ) in the earth’s crust in southwest (SW) Anatolia are estimated by using the coda wave method based on the decrease of coda wave amplitude by time on the seismogram. The quality factor Q _o, the value of Q _c at 1 Hz, and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves. δ is determined from the observations of amplitude variations of seismic waves. In applying the coda wave method, firstly, a type curve representing the average pattern of the individual coda decay curves for 0.75, 1.5, 3.0, 6.0, 12.0, and 24.0 Hz values was estimated. Secondly, lateral variation of coda Q and the attenuation coefficients for three main tectonic patterns are estimated. The shape of the type curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The Q _o and η values vary from 30 to 180 and from 0.55 to 1.25, respectively for SW Anatolia. In SW Anatolia, coda Q–f relation is described by and δ = 0.008 km⁻¹. These results are expected to help in understanding the degree of tectonic complexity of the crust in SW Anatolia.     

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.     

各向异性半空间中浅埋孔洞对地表 反平面运动的影响

本文采用复变函数方法研究了稳态水平剪切波(SH波)在各向异性弹性半空间中任意形状孔洞上的散射及其对地面运动的影响.在变换空间中构造出自动满足各向异性半空间水平表面上应力为零的散射波函数,并根据孔洞表面应力为零的边界条件,采用最小二乘法求解散射波函数的系数.用介质的各向异性性质来模拟地质条件,给出了SH波作用下含圆孔、椭圆孔和方孔的各向异性半空间表面位移幅值的数值结果,并分析了介质特性、孔洞的形状、埋深、入射波波数及入射角度等因素对地面运动的影响规律.数值结果表明:介质的各向异性对含有孔洞的半空间表面的地表位移具有显著的影响;沿一定角度的入射波在某一频段内所引起的地表位移幅值比各向同性介质的可能要大,且随着孔洞埋深的增加,地表位移的幅值逐渐减小.      

Determination of source duration and moment of earthquakes in the Indian Peninsula using love waveforms

Love waves recorded by a long-period seismograph at New Delhi (NDI) from seven earthquakes of magnitude 4.3 to 5.2 in Koyna and Bhatsa on the western coast and one earthquake in Ongole on the eastern coast of the Indian Peninsula have been used to determine the seismic moment for each of the earthquakes by waveform modeling. Transverse component of the synthetic seismogram shows that the maximum amplitude of waveform decreases with an increase of source duration. Thus for an evaluation of the seismic moment by equating the amplitude level of the observed and synthetic waveforms, we must know the source duration. The synthetic seismogram also indicates that a short source duration gives rise to a small but sharp pulse and this pulse is interpreted as anLg wave. Comparison of the observed and synthetic waveforms has been used for a simultaneous evaluation of the source duration and seismic moment. The source durations are found to vary between 2.2 and 4.4 s; for earthquakes with a magnitude range between 4.3 and 5.2 these durations are slightly higher than normal. We obtain moment (M ₀) of Ongole earthquake (M _L =5.1)as 1.7×10²⁴ dyne-cm; moments of Koyna and Bhatsa earthquakes (4.3M _L 5.2) on the western coast lie between 0.7×10²³ and 3.6×10²³ dyne-cm. Moment (M ₀)-magnitude (M _L ) relation logM ₀=1.5M _L +16.0 for the western United States region agrees as well, in general, with the results for the earthquakes in the Indian Peninsula.     

A review of numerical experiments on seismic wave scattering

This paper reviews applications of the finite-difference and finite-element methods to the study of seismic wave scattering in both simple and complex velocity models. These numerical simulations have improved our understanding of seismic scattering in portions of the earth where there is significant lateral heterogeneity, such as the crust. The methods propagate complete seismic wavefields through highly complex media and include multiply scattered waves and converted phases (e.g.,P toSV, SV toP, body wave to surface wave). The numerical methods have been especially useful in cases of moderate and strong scattering in complex media where multiple scattering becomes important. Progress has been made with numerical methods in understanding how near-surface, low-velocity basin structures scatter surface waves and vertically-incident body waves. The numerical methods have proven useful in evaluating scattering of surface waves and body waves from topography of both the free surface and interfaces buried at depth. Numerical studies have demonstrated the importance of conversions from body waves to surface waves (andvice versa) when lateral heterogeneities and topographic relief are present in the uppermost crust. Recently, several investigations have applied numerical methods to study seismic wave propagation in velocity models which vary randomly in space. This stochastic approach seeks to understand the effects of small-scale complexity in the earth which cannot be resolved deterministically. These experiments have quantified the relationships between the statistical properties of the random heterogeneity and the measurable properties of high-frequency (1 Hz) seismograms. These simulations have been applied to the study of many features observed in actual high-frequency seismic waves, including: the amplitude and time decay of seismic coda, the apparent attenuation from scattering, the dispersion of waveforms, and the travel time and waveform variations across arrays of receivers.     

空间尘埃等离子体电磁散射特性研究

基于尘埃等离子体Mie-Debye散射模型,研究平衡态下尘埃等离子体中带电尘埃的电磁散射特性.根据输运理论计算高层大气环境中尘埃等离子体层的电磁散射场,给出电磁波垂直入射时尘埃等离子体层的反射函数,同时分析不同浓度、不同粒径分布状态下尘埃等离子体层对电磁波传播的影响.结果表明:尘埃粒子尺度对其电磁散射特性影响非常大;在对数正态分布情况下,当入射波长远大于Debye半径时,尘埃等离子体的散射主要表现为Debye散射,而入射波长远小于Debye半径时,Mie散射占主要部分;并且尘埃等离子体层对电磁波的衰减随尘埃粒子浓度以及尘埃粒子半径的增大而明显增强.     

Electromagnetic scattering by a perfectly conducting half-plane in a conductive half-space

FollowingDmitriev (1960) a rigorous theoretical solution for the problem of scattering by a perfectly conducting inclined half-plane buried in a uniform conductive half-space has been obtained for plane wave excitation. The resultant integral equation for the Laplace transform of scattering current in the half-plane is solved numerically by the method of successive approximation. The scattered fields at the surface of the half-space are found by integrating the half-space Green's function over the transform of the scattering current.The effects of depth of burial and inclination, of the half-plane on the scattered fields are studied in detail. An increase in the depth of burial leads to attenuation of the fields. Inclination introduces asymmetry in the field profiles beside affecting its magnitude. Depth of exploration is greater for quadrature component. An interpretation scheme based on a phasor diagram is presented for the VLF-EM method of exploration for rich vein deposits in a conductive terrain.List of symbols x, y, z Space co-ordinates - Half-space conductivity - ₀ Free-space permeability - Excitation frequency (angular) - T Time - h Depth of the half-plane - a Inclination of the half-plane - E _x x-Directed total electric field - E _x ^p x-Directed primary electric field - E _xo ^p x-Directed primary electric field atz=0 directly over the half-plane - H _y y-Component of total magnetic field - H _y ^p y-Component of primary magnetic field - H _y0 ^p y-Component of primary magnetic field atz=0 directly over the half-plane - H _z z-Component of total magnetic field - H _z ^p z-Component of primary magnetic field - J _x Surface density ofx-directed scattering current - G Green's function - k ₀,K Wave numbers - u,u ₀,u ₁,u ₂ Functions - Space co-ordinate - s Variable in transform domain - Variable of integration - Normalized scattering current - Laplace transform of - N Normalized - , ₀, ₁, ₂ Functions - t Variable of integration - Skin depth - H Total magnetic field - H ^p Primary magnetic field - H ₀ ^p Primary magnetic field atz=0 directly over the half-plane - M,Q,R,S,U,V Functions - N ₁,N ₂ Functions