共查询到20条相似文献，搜索用时 31 毫秒
1.
Seismic body waves in anisotropic media: synthetic seismograms 总被引：5，自引：0，他引：5
Summary. Synthetic seismograms and particle motion diagrams are computed for simple, layered Earth models containing an anisotropic layer. The presence of anisotropy couples the P, SV and SH wave motion so that P waves incident on the anisotropic layer from below produce P, SV and small-amplitude SH waves at the surface both the P velocity and the amplitudes of the converted phases vary with azimuth. Significant SH amplitudes may be generated even when the wavelength of the P wave is much greater than the thickness of the anisotropic layer. Incident SV or SH waves may each generate large amplitudes of both SV and SH motion. This strong coupling is largely independent of the degree of velocity anisotropy of the medium. The arrivals from short-period S waves exhibit S-wave splitting, but arrivals from longer period S waves superpose into a modified waveform. This strong coupling does not allow the arrival of separate phases with pure SV and SH polarization except along directions of symmetry where the motion decouples. 相似文献
2.
R. A. Stephen 《Geophysical Journal International》1977,51(1):169-181
Synthetic seismograms are shown and discussed for the case of the receiver within the medium. Most of the discussion is on the reflectivity method with the receiver within the reflectivity zone, but results using the ray method are shown for comparison. Such synthetic seismograms can be used to interpret data from Oblique Seismic Experiments where shots generated on the surface up to large ranges are recorded in crustal boreholes. 相似文献
3.
Michael S. Thorne Thorne Lay Edward J. Garnero Gunnar Jahnke † Heiner Igel 《Geophysical Journal International》2007,170(2):635-648
We use an axisymmetric, spherical Earth finite difference algorithm to model SH -wave propagation through cross-sections of laterally varying lower mantle models beneath the Cocos Plate derived from recent data analyses. Synthetic seismograms with dominant periods as short as 4 s are computed for several models: (1) a D" reflector 264 km above the core–mantle boundary with laterally varying S -wave velocity increases of 0.9–2.6 per cent, based on localized structures from a 1-D double-array stacking method; (2) an undulating D" reflector with large topography and uniform velocity increase obtained using a 3-D migration method and (3) cross-sections through the 3-D mantle S -wave velocity tomography model TXBW. We apply double-array stacking to assess model predictions of data. Of the models explored, the S -wave tomography model TXBW displays the best overall agreement with data. The undulating reflector produces a double Scd arrival that may be useful in future studies for distinguishing between D" volumetric heterogeneity and D" discontinuity topography. Synthetics for the laterally varying models show waveform variability not observed in 1-D model predictions. It is challenging to predict 3-D structure based on localized 1-D models when lateral structural variations are on the order of a few wavelengths of the energy used, particularly for the grazing geometry of our data. Iterative approaches of computing synthetic seismograms and adjusting model characteristics by considering path integral effects are necessary to accurately model fine-scale D" structure. 相似文献
4.
Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic 总被引：19，自引：0，他引：19
Stuart Crampin 《Geophysical Journal International》1978,53(3):467-496
summary . A new technique is presented for modelling the elastic constants of cracked structures with application to systems with weak concentrations of parallel cracks, and of simple biplanar and triplanar cracks. The velocities and Vp /Vs ratios of these anisotropic structures are used to provide quantitative models for some earthquake precursors. These results indicate the great importance of crack geometry to the behaviour of precursors. The velocities of saturated cracks appear to favour the dilatancy-diffusion model of precursory phenomena. Synthetic seismograms are calculated for propagation through possible dilatancy zones. The seismograms show some characteristic features which may be useful for the investigation of earthquake dilatancy. 相似文献
5.
B. R. Cassell 《Geophysical Journal International》1982,69(2):339-354
Summary An effective algorithm for computing synthetic seismograms in laterally inhomogeneous media has been developed. The method, based on zero-order asymptotic ray theory, is primarily intended for use in refraction and reflection studies and provides an economical means of seismic modelling.
A given smoothed velocity-depth-distance model is divided into small squares with constant seismic parameters and first-order interfaces are represented by an arbitrary number of dipping linear segments. The computation of ray propagation and amplitudes through such a model does not involve complicated analytic expressions and therefore minimizes computer time.
Amplitudes are determined by geometrical spreading of spherical wave-fronts and energy partitioning at interfaces. Synthetic seismograms calculated for laterally homogeneous models are in good agreement with those obtained by the Reflectivity Method. 相似文献
A given smoothed velocity-depth-distance model is divided into small squares with constant seismic parameters and first-order interfaces are represented by an arbitrary number of dipping linear segments. The computation of ray propagation and amplitudes through such a model does not involve complicated analytic expressions and therefore minimizes computer time.
Amplitudes are determined by geometrical spreading of spherical wave-fronts and energy partitioning at interfaces. Synthetic seismograms calculated for laterally homogeneous models are in good agreement with those obtained by the Reflectivity Method. 相似文献
6.
Summary. A lower mantle S-wave triplication detected with short- and long-period WWSSN and CSN recordings indicates a substantial shear velocity discontinuity near 280 km above the core–mantle boundary. The triplication can be observed in rotated SH seismograms from intermediate and deep focus events throughout the distance range from 70° to 95°. Three distinct source region–receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication, with slight systematic shifts in the triplication indicating up to 40 km variations in the depth of the discontinuity. Modelling of the observations with synthetic seismograms produced with the Cagniard de Hoop and reflectivity methods constrains the shear velocity increase to be 235 ± 0.25 per cent, comparable to upper mantle discontinuities. Short-period observations indicate that the velocity increase may be a sharp first-order discontinuity, or may extend over a transition zone no more than 50 km thick. The shear velocity gradient below the discontinuity, within the D" layer, is not well-constrained by the SH data, but slightly positive or near zero velocity gradients are consistent with the long-period amplitude ratios of ScSH/SH . 相似文献
7.
Summary. We report the initial results of our attempts to obtain theoretical seismograms for direct comparison with the experimental time series obtained with the long-period instruments of the WWSSN. The entire theoretical seismogram, including both body waves and surface waves, can be generated for a spherical, anelastic earth by simple inverse Fourier transformation of the sum of the propagating fundamental and higher-mode surface waves. The key to success in reproducing the WWSSN records involves the number of modes, and the minimum period used in these computations; here we use eight modes and a minimum period of 2 s. Efficient computational algorithms make it possible to handle up to 2000 frequency points for each mode; approximately 200 layers are used to model the radial heterogeneity of the earth; attenuation is treated exactly. Examples are given of the SH theoretical seismograms resulting from dislocation sources buried at various depths in the Earth. 相似文献
8.
Summary. Seismic investigations using shear-wave and converted wave techniques show that very often reflected PS - and SS -waves have anomalous polarizations ( accessory components ). This phenomenon cannot be explained in terms of isotropic models with dipping boundaries. Computations of synthetic seismograms of reflected PS - and SS -waves were made for different models of transversely isotropic media with dipping anisotropic symmetry axes not normal to the boundaries. Synthetic seismograms were computed by ray techniques using an optimization algorithm to construct all rays arriving at a given receiver. These computations indicate that accessory components arise when the medium above the boundary is anisotropic, where they are caused by the constructive interference of qSV - and qSH -waves. If a low-velocity layer is present, displacement vectors of both waves have horizontal projections which are approximately orthogonal. The algorithm for wave separation is presented and some results of its use are given. 相似文献
9.
We investigate the importance of wide-band coupling of normal modes due to inner core anisotropy. We compare four different seismic models of inner core anisotropy, which were obtained by others using the splitting of Earth's normal modes. These models have been developed using a self-coupling (SC) approximation, which assumes that coupling between nearby modes through anisotropic inner core structure is negligible. We test the SC approximation by comparing the frequencies and quality factors of 90 inner core sensitive modes, computed for these models using either the SC approximation or full-coupling (FC) among large groups of modes. We find significant shifts in the quality factors and frequencies for some modes. Groups of modes which significantly couple together are constructed for six target modes. These groups are model dependent and in some cases contain large numbers of modes. Synthetic seismograms are calculated to show that the difference between SC and FC is observable on the scale of seismograms and of the same order of magnitude as the difference between synthetic and observed seismograms. Thus, future models of inner core anisotropy should take cross-coupling between large groups of modes into account. 相似文献
10.
Joseph Ha 《Geophysical Journal International》1984,78(1):291-305
Summary. A set of recurrence relations similar to that of Kennett suitable for SH -wave generation in an ( n + l)-layered half-space is presented. The recurrence relations contain no growing terms, thus providing a stable and efficient algorithm for computing complete SH synthetic seismograms. The complete expansion of these recurrence relations gives the explicit form of the transfer function for SH -waves. The transfer function for a point source in layer s of the stratification is a series of 2 n terms in the denominator and a series of 2 n−s +1 terms in the numerator. The result of Wang from ray summation is shown to be a special case of our general result. Numerical comparison of the algorithm of this paper with the propagator matrix method is also made. 相似文献
11.
T. J. Clarke 《Geophysical Journal International》1984,78(3):775-793
Summary. We develop a méthod of reconstructing the elastic paraméters as functions of depth, for a horizontally stratified, isotropic elastic half-space. Unlike previous schemes, which have been able to retrieve the shear wave speed and density from SH seismograms slant stacked at two angles, our méthod makes use of P - SV data at a single stacking paraméter to obtain all three elastic constants. The data required are the elements of the full reflection matrix at the surface, corresponding to measurements of two separate components of the response to two independent sources, one explosive, the other generating shear waves.
In developing this inverse scheme fundamental differences emerge between the acoustic or SH problem, and the coupled P - SV case, the most important being in the nature of the interfacial scattering matrix. We show that it is not possible to make use of the downward reflection data for an interface to determine directly the remaining reflection and transmission coefficients, but that the scattering data may be completed by applying a simple iterative procedure at each interface.
We show the result of applying our inverse scheme to seismograms generated for a six-layered model, including a low-velocity layer. We are able to reconstruct both wave speeds and the density as functions of depth, all quantities being in close agreement with the original model. 相似文献
In developing this inverse scheme fundamental differences emerge between the acoustic or SH problem, and the coupled P - SV case, the most important being in the nature of the interfacial scattering matrix. We show that it is not possible to make use of the downward reflection data for an interface to determine directly the remaining reflection and transmission coefficients, but that the scattering data may be completed by applying a simple iterative procedure at each interface.
We show the result of applying our inverse scheme to seismograms generated for a six-layered model, including a low-velocity layer. We are able to reconstruct both wave speeds and the density as functions of depth, all quantities being in close agreement with the original model. 相似文献
12.
B. G. Mikhailenko 《Geophysical Journal International》1984,79(3):963-986
Summary. An algorithm which is part analytical and part numerical is suggested for the computation of complete synthetic seismograms for complex three-dimensional geological structures with radial symmetry. A partial separation of variables based on the combination of a finite Fourier integral transform with respect to the spatial coordinate z together with the finite difference method is the essence of the algorithm. Upon application of the finite transform the problem reduces to solving a system of equations containing only partial derivatives with respect to one spatial coordinate ( r ) and time. As radial symmetry is assumed, there is no functional dependence on φ in the cylindrical system of coordinates ( r , φ, z ). The coefficients of the transformed equations may contain finite Fourier integrals of the z dependence of the elastic parameters. Several examples of synthetic seismograms computed for both SH - and P – SV -waves propagating in complex subsurface geometries are presented and their interpretation discussed. 相似文献
13.
David W. Eaton 《Geophysical Journal International》2006,166(3):1249-1258
The Kirchhoff (or tangent plane) approximation, derived from the theoretically complete Kirchhoff–Helmholtz integral representation for the seismic wavefield, has been used extensively for the analysis of seismic-wave scattering from irregular interfaces; however, the accuracy of this method for curved interfaces has not been rigorously established. This paper describes an efficient Kirchhoff algorithm to simulate scattered waves from an arbitrarily curved interface in an elastic medium. Synthetic seismograms computed using this algorithm are compared with exact synthetics computed using analytical formulae for scattering of plane P waves by a spherical elastic inclusion. A windowing technique is used to remove strong internal reverberations from the analytical solution. Although the Kirchhoff method tends to underestimate the total scattering intensity, the accuracy of the approximation improves with increasing value of the wavenumber-radius product, kR . The arrival times and pulse shapes of primary reflections from the sphere are well approximated using the Kirchhoff approach regardless of curvature of the scattering surface, but the amplitudes are significantly underestimated for kR ≤ 5. The results of this work provide some new guidelines to assess the accuracy of Kirchhoff-synthetic seismograms for curved interfaces. 相似文献
14.
Joel M. Crichlow 《Geophysical Journal International》1982,70(3):563-575
Summary. Two-dimensional finite element modelling of underground structural anomalies at shallow depths has been done to obtain the response at the ground surface to damped, vertically incident, SH -waves. Power spectral ratios are examined to determine what effects the position, shape, depth and size of the anomaly have on the surface seismograms. Based on the results gathered from a number of models, inferences are made with respect to the inverse problem: given the seismic motion of the ground surface, what can be said about the underground structure? 相似文献
15.
Summary . Synthetic seismograms represented by integrals generally display signals associated with the limits of integration. Sometimes these 'end-point' contributions are spurious (e.g. in the WKBJ seismogram) and sometimes they are the main physical interest (e.g. the Kirchhoff integral for an edge). The end-point contributions may be asymptotically approximated using integration by parts or Laplace's method and it may then be possible to reduce them if desired. We describe examples in the WKBJ seism ogram for reflected or transmitted waves in homogeneous layers and for turning waves. We also study signals due to discontinuities in reflection coefficients, by partitioning the real slowness integral so that the discontinuities lie at end points. Examples are the head wave, which is a physically correct signal, and spurious diffractions caused by using plane-wave coefficients for grazing rays in the WKBJ seismogram. 相似文献
16.
Peter R. Shaw 《Geophysical Journal International》1986,85(2):291-313
Summary . In this paper the accuracy of velocity-depth profiles derived by matching WKBJ seismograms to observations is quantitatively evaluated. Seismograms computed with the WKBJ method are generally quite reliable but possess predictable, systematic inaccuracies in the presence of strong velocity gradients. The effects of these inaccuracies on models derived through WKBJ waveform inversion are studied, using reflectivity seismograms as 'data'. The velocity structure used is an oceanic lithosphere model that contains several transition regions separated by relatively homogeneous layers, producing partially-reflected reverberations in the reflectivity synthetics that are absent from the WKBJ seismograms. The inversion incorporates the 'jumping' strategy to solve for the smoothest models consistent with the data. We find these solutions to be independent of the starting model and to have a stable basic structure that agrees well with the correct model. The differences, everywhere less than a seismic wavelength, depend on the frequency content of the seismograms. Reverberations in the reflectivity seismograms that are well separated from WKBJ arrivals are treated as 'noise' in the inversion. 相似文献
17.
Summary. The Lanczos method of separating exponentials is applied to the Fourier transform of seismograms in order to separate the various modes that contribute to the given portion of the seismograms. Phase velocities and amplitudes are obtained as functions of the frequency. When applying the method to artificial seismograms, which are built as an exact superposition of a number of modes, the separation is very accurate. The method was also applied to the surface-wave portion of numerical seismograms for a vertical point force in a layered medium. The phase velocity and amplitude of the fundamental mode are obtained. These functions were taken as the first guess in the Backus—Gilbert generalized inverse procedure and the process converged very rapidly. When a perturbation of the phases and amplitudes is taken as the first guess the process converges to the true model when enough data are available. 相似文献
18.
Offset-dependent characteristics of seismic scattering are useful for characterizing fractured reservoirs. We use two models that have different background medium properties and different azimuthal AVO responses to study elastic wave propagation and scattering in gas-saturated, heterogeneously fractured reservoirs. Heterogeneous fracture density distributions are built through stochastic modelling. Synthetic seismograms are generated by 3-D finite difference modelling, and waveforms along crack-normal and strike directions are considered in this paper. The multiple signal classification (MUSIC) frequency estimator is used in waveform estimation to provide frequency-domain attributes related to seismic wave scattering by fracture heterogeneity. Our results indicate that the strength of the scattering field is a function of the background medium. The strength also increases with increasing fracture scatterer density and with decreasing correlation length of spatial variations of fracture density. The scattering field is weak at the top of the fractured reservoir. The first-order results are dominated by velocity anisotropy of the mean fracture density field. However, the base of the fractured reservoir corresponds to a strong scattering field on which fracture heterogeneity has a larger effect and is characterized by the loss of coherence. 相似文献
19.
Summary. The study of the asymptotic behaviour of eigenfrequencies of torsional modes of the Earth is of interest in the problem of the existence of discontinuities in the Earth's interior and the determination of their depth and scale, since the solotone effect, which is a persistent oscillatory cohponent in the asymptotic overtone structure, is very sensitive to the presence of discontinuities. The asymptotic behaviour of torsional eigenfrequencies of the second order differential equation for the Earth's free oscillations can be compared with that obtained from eigenfrequencies evaluated from synthetic SH-wave seismograms by Brune's phase correlation method, using various earth models. The solotone effect that appears in the former for discontinuous models can be explained in terms of multiple reflections fram the discontinuities, and can be reconstructed from synthetic SH -wave pulses arising from these reflections. Its properties vary systematically with the depth and the scale of discontinuities and can be superposed for several discontinuities. 相似文献
20.
A fault plane solution using theoretical P seismograms 总被引：1，自引：0，他引：1
We use the method of Hudson and Douglas, Hudson & Blarney to compute seismograms which simulate the codas of 10 short period P -wave seismograms from a shallow earthquake. The polarities and relative amplitudes of P and pP measured from seven of the observed seismograms are used to compute a fault plane solution with confidence limits, assuming that the source radiates as a double couple. This solution is in approximate agreement with that given for the same earthquake by Sykes & Sbar, who used only the onset polarities of short-period P waves. The small difference between the two solutions can be explained by interference between the true first motion of P and microseismic noise at two stations.
The results show that, for some shallow earthquakes, the relative amplitude method has the following advantages over the first motions method. First, a P/pP amplitude ratio (with appropriate confidence limits) can always be measured, even in seismograms which are so noisy that the first motion of P is uncertain. Second, the fault plane solutions obtained from relative amplitudes have known confidence limits. Finally, by using more information from each seismogram, the relative amplitude method requires considerably fewer seismograms than the first motions method. 相似文献
The results show that, for some shallow earthquakes, the relative amplitude method has the following advantages over the first motions method. First, a P/pP amplitude ratio (with appropriate confidence limits) can always be measured, even in seismograms which are so noisy that the first motion of P is uncertain. Second, the fault plane solutions obtained from relative amplitudes have known confidence limits. Finally, by using more information from each seismogram, the relative amplitude method requires considerably fewer seismograms than the first motions method. 相似文献