首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 250 毫秒
1.
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 2 /V s 4 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.  相似文献   

2.
We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ –1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ –1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ –1 is proportional tok –1 independently of model parameters and the wave modes. If the seismological observation thatQ –1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.  相似文献   

3.
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 2. 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.  相似文献   

4.
Numerical modelling techniques are now becoming common for understanding the complicated nature of seismic wave propagation in fractured rock. Here the Indirect Boundary Element Method (IBEM) is applied to study scattering of elastic waves by cracks. The problem addressed in this paper is the diffraction of P and S waves by open 3-D cracks of arbitrary shape embedded in a homogeneous isotropic medium. The IBEM yields the value of the jump of displacements between opposite surfaces of the crack, often called Crack Opening Displacement (COD). This is used to evaluate the solution away from the crack. We use a multi-regional approach which consists of splitting a surface S into two identical surfaces S+ and S− chosen such that the crack lies at the interface. The resulting integral equations are not hyper-singular and wave propagation within media that contain open cracks can be rigorously solved. In order to validate the method, we compare results of displacements of a penny-shaped crack for a vertical incident P-wave with the classic results by Mal (1970) obtaining excellent agreement. This comparison gives us confidence to study cases where no analytic solutions exist. Some examples of incidence of P or S waves upon cracks with various shapes are depicted and the salient aspects of the method are also discussed. Both frequency and time-domain results are included.  相似文献   

5.
Guided seismic waves are generated at lateral discontinuities in coal seams as an integral part of the mechanics of scattering downgoing body waves from a surface energy source.Numerical and laboratory model studies undertaken for both thin and thick coal seam waveguides show that small faults generate recognizable diffraction patterns by convertingP andS waves into channel wave arrivals at in-seam geophones on a walkawayVSP profile.These results suggest that in-seam geophones can be used during subsequent surface reflection surveys so as to improve mine profitability by optimising longwall layout in faulted ground.  相似文献   

6.
Parameters of seismic waves from clusters of local weak earthquakes that occurred at the upper boundary of the seismofocal zone in depth intervals of 40–60 and 70–90 km along eastern Hokkaido are investigated for the period 1998–2003, including the strong (M = 8.0) Tokachi-Oki earthquake of September 26, 2003. Analysis of data indicates that the distribution of anisotropic properties along Hokkaido is inhomogeneous and parameters of split waves (the azimuth of the fast S wave and the time delay between split S waves) are sensitive to variations in the stress-strain state of the medium. Unstable behavior of split wave parameters and increased values of the ratio V P /V S for clusters of events in the areas of the Hidaka Mountains and Nemuro Peninsula imply that the medium is in a mechanically weakened state (in the regime of intense dilatant deformation). On the contrary, the regions beneath the Tokachi and Kushiro plains are more rigid and are characterized by lower V P /V S values and comparatively stable behavior of wave parameters. Anomalous parameters of split waves from events of clusters in areas of different stations correlate with each other and are related to occurrence times of large earthquakes around Hokkaido, which may point to a redistribution of stresses and strains and fluid migration in the subduction zone.  相似文献   

7.
A single scattering model was used to analyse the temporary changes in the mean density of scattered waves in a discrete random medium. The model of the mean energy density, originally proposed bySato (1977) for spherical radiation and isotropic scattering, has been modified and applied to a medium in which the scatterers are confined to a specified volume. The time variation of the early part of the mean energy density function for the different source durations was investigated. The dominant effect on the theoretical mean energy density is caused by the specified volume containing scatterers. The duration of the source pulse influences the early part of the coda fort/t 0<1.2, wheret is the lapse time measured from the source origin time, andt 0is arrival time of the body wave.The analysis of the coda signal of micro-events occurring immediately in front of the face enables us to estimate the size of the fracture zone induced by the stope. The model of the mean energy density of coda for a medium containing scatterers close to the seismic source was used to analyse a large number of events recorded close to an advancing mine face in a deep level gold mine in South Africa. The coda decay rate has two trends: the first, with a steep decay of coda, is produced by a larger deviation of rock parameters and/or larger size of the scatterers; the second trend, which decays more slowly, has the corresponding mean-free path ranging from 20 m to 200 m. The analysis indicates that the rock mass about 15–20 m from the stope contains a large proportion of fractured and blocked rock, which is the source of scattering. The scattering of theS-wave was much stronger and more stable, with the mean-free path varying from 11 m to 45 m. This is due to the shorter wavelength of theS wave in comparison with theP wave. The quality factor for theP coda wave varies from 30 to 100 in the fracture zone of stope and outside this zone it has a value of 300. The quality factor of theS wave varies from 20 to 78 in the equivalent volume. For rock surrounding the stope the ratioQ sp –1 /Q ss –1 varied from 0.31 to 0.69. This suggests that the radii of scatterers are smaller than 3.5 m.  相似文献   

8.
Summary Analytic expressions for the far-field displacement both in the time and the frequency domains are obtained for the kinematic earthquake model of a rupture which nucleates at a point, grows in the form of a rectangle and stops at a later time.The spectra ofS andP waves are analyzed with the object of determining fault dimensions from the corner frequencies. The approach by which the solution of the present problem is obtained can be extended to get closed form solutions to the problem of rupture front expanding in the form of a polygon and terminating over a fault area which has the form of another polygon.  相似文献   

9.
Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ –1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ –1 and intrinsicQ –1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves.  相似文献   

10.
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]2) are also confirmed.  相似文献   

11.
Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic  相似文献   

12.
We demonstrate how multiples, generated at the interfaces of plane parallel beds, modify the propagation characteristics of an originally coherent seismic wave. For waves propagating at an angle to the bedding plane we find that theSV andP-waves couple so that neither is a pure mode. TheSH-wave, while modified in its propagation characteristics by multiples, remains a pure mode. The coupling ofSV-multiples into the quasi-P-mode appears weaker than the coupling ofP-wave multiples into the quasi-SV mode; at least this is so for the two simple cases of (a) density fluctuations only and (b) correlatedV p andV s fluctuations which conserve Poisson's ratio.We also find that the coupling is sensitive to both the angle of propagation and frequency. In addition there is a cut-off angle forP-wave multiples influencing the quasi-SV mode. Propagation angles larger than the cut-off permit theP-multiples to modify the phase of the quasi-SV mode, but not its effective attenuation. No such cut-off effect is found for SV-multiples influencing the quasi-P mode, whose angle-dependent and frequency-dependent phase distortion and effective attenuation are influenced both byP-wave multiples andSV-multiples.In view of the mathematical complexity of the expressions describing the phase, and effective attenuation of modes when allowance is made forP-andS-wave multiples, we strongly advocate numerical coding of the major mathematical formulae. By so doing a systematic study can be undertaken of the frequency and offset dependence of seismic waves as a function of seismic source input and power spectral behavior of the fluctuations in density and elastic constants of beds. It is our opinion that the full mathematical expressions are too involved to permit an analytic, systematic investigation to be given of the phase and attenuation of seismic waves with any degree of sophistication or generality.  相似文献   

13.
  相似文献   

14.
Interpretation ofP/S-wave successions is used in induced or passive microseismicity. It makes the location of microseismic events possible when the triangulation technique cannot be used. To improve the reliability of the method, we propose a technique that identifies theP/S-wave successions among recorded wave successions. A polarization software is used to verify the orthogonality between theP andS polarization axes. The polarization parameters are computed all along the 3-component acoustic signal. Then the algorithm detects time windows within which the signal polarization axis is perpendicular to the polarization axis of the wave in the reference time window (representative of theP wave). The technique is demonstrated for a synthetic event, and three application cases are presented. The first one corresponds to a calibration shot within which the arrivals of perpendicularly polarized waves are correctly detected in spite of their moderate amplitude. The second example presents a microseismic event recorded during gas withdrawal from an underground gas storage reservoir. The last example is chosen as a counter-example, concerning a microseismic event recorded during a hydraulic fracturing job. The detection algorithm reveals that, in this case, the wave succession does not correspond to aP/S one. This implies that such an event must not be located by the method based on the interpretation of aP/S-wave succession as no such a succession is confirmed.  相似文献   

15.
There remains much uncertainty on the absolute elastic wave energy released by fault rupture. Few direct estimates of the partition of seismic wave energy in ground shaking have been made. In this work, ground particle velocities from integrated accelerograms are used to compute the kinetic energy crossing unit area per unit time. Simplified theory for the near-field strong-motion case would appear to give a valid lower energy bound; the wave attenuation does not present a major problem. The partition of energy in predominantly P, S, and surface wave portions, for given frequency windows, is tabulated using strong-motion accelerograms recorded at different azimuths to the fault-sources of six California earthquakes (5.5<M L<7.2). Checks against earlier magnitudeM L and momentM 0 correlations indicate significantly higher overall wave energy outputs than expected, but further calibration is needed.The study demonstrates that stable estimates of frequency-dependent seismic wave energies in the nearfield can be obtained from strong-motion records. Hence, energy flux may have a wider application to risk mapping than previously thought. In particular, a shift from peak acceleration scaling to (kinetic) energy inputs for engineering design appears to involve only routine processing.  相似文献   

16.
Results from plasma wave experiments in spacecrafts give support to nonlinear interactions involving Langmuir, electromagnetic, and ion-acoustic waves in association with type III solar radio bursts. Starting from a general form of Zakharov equation (Zakharov, V.E., 1985. Collapse and self-focusing of Langmuir waves. Hand-book of Plasma Physics Cap.2, 81–121) the equations for electric fields and density fluctuations (density gratings) induced by a pair of counterpropagating Langmuir waves are obtained. We consider the coupling of four triplets. Each two triplets have in common the Langmuir pump wave (forward or backward wave) and a pair of independent density gratings. We numerically solve the dispersion relation for the system, extending the work of (Alves, M.V., Chian, A.C.L., Moraes, M.A.E., Abalde, J.R., Rizzato, F.B., 2002. A theory of the fundamental plasma emission of type- III solar radio bursts. Astronomy and Astrophysics 390, 351–357). The ratio of anti-Stokes (AS) (ω0+ω) to Stokes (S) (ω0-ω*) electromagnetic mode amplitudes is obtained as a function of the pump wave frequency, wave number, and energy. We notice that the simultaneous excitation of AS and S distinguishable modes, i.e., with Re{ω}=ωr≠0, only occurs when the ratio between the pump wave amplitudes, r is ≠1 and the pump wave vector k0 is , W0 being the forward pump wave energy. We also observe that the S mode always receives more energy.  相似文献   

17.
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.  相似文献   

18.
The Alaska Tsunami Warning Center has the responsibility of providing timely tsunami warning services for Alaska and the west coasts of Canada and the United States. Recently, the ATWC implemented a new microcomputer system which is used for both automatic and interactive earthquake processing, and for disseminating critical information to the Tsunami Warning System recipients.Real-time seismic wave form data from 23 short-period and 9 long-period sites in Alaska, the lower 48 States, and Hawaii, are continually computer-monitored for the occurrence of an earthquake. Once detected from the short-period wave form data, pre- and post-earthquake data are displayed on a graphics terminal along with an indicator to identify the time of the onset of theP waves (P-picks). TheP-picks can easily be changed during or after data collection via a mouse. Magnitudes (M b ,M l ,M B ,M S ) are automatically computed from appropriate short- and long-period wave form data concurrently with the above processing. A second graphics terminal displays cycle-by-cycle long-period wave form data that was used to compute an earthquake'sM B andM S magnitudes.An earthquake's parametric data and other information are available and printed within tens of seconds after theP wave arrivals are recorded at the first 5 sites, then 7 sites, 9 sites, and a final parametric computation using all collected data. Three video display monitors are used for displaying the parameters, procedural aids, and a map showing the epicenter. Additionally, selected event parameters are immediately transmitted by VHF radio to alphanumeric beepers which are carried by standby duty personnel during those times that the Center is not manned.Using a dedicated video display terminal and printer, the interactive system can use data and parameters resulting from the automatic processes for concurrent parameter recomputations; perform additional computations; disseminate critical information; and generate procedural aids for duty geophysicists to facilitate an earthquake/tsunami investigation.  相似文献   

19.
In the present article, the dependencies of the acoustic signal total energy and the energies of the wave packets of different types of the waves on the elastic parameters and permeability of rocks have been studied. We have considered traditional logging tools containing acoustical monopole source. Calculations were performed in a frequency range of dozens of kilohertz, typical for acoustic well logging. It was shown that in a typical high-velocity formation (vs > vf, where vs and vf are the velocities of the shear wave in the rock and of the compressional wave in the borehole fluid, respectively), the pseudo-Rayleigh waves, whose elastic properties depend mainly on the shear modulus of the rock, contributed significant energy to the total signal energy in the borehole. The energies of different wave packets depend on the permeability in different ways. The greatest sensitivity to permeability change has been shown by the acoustic signal total energy and the energy of the low-velocity part of the pseudo-Rayleigh wave packet. The theoretical analysis was illustrated by real sonic log data.  相似文献   

20.
Acoustic full waveforms recorded in wells are the simplest way to get the velocity of P, S, and Stoneley waves in situ. Processing and interpretation of acoustic full waveforms in hard formations does not generate problems with identification packets of waves and calculation of their slowness and arrivals, and determination of the elastic parameter of rocks. But in shallow intervals of wells, in soft formations, some difficulties arise with proper evaluation of the S-wave velocity due to the lack of refracted S wave in case when its velocity is lower than the velocity of mud. Dynamic approach to selection of a proper value of semblance to determine the proper slowness and arrival is presented. Correlation between the results obtained from the proposed approach and the theoretical modeling is a measure of the correctness of the method.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号