首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 984 毫秒
1.
Finite-frequency sensitivity kernels for head waves   总被引：2，自引：0，他引：2  
Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.  相似文献   

2.
We investigate the sensitivity of finite-frequency body-wave observables to mantle anisotropy based upon kernels calculated by combining adjoint methods and spectral-element modelling of seismic wave propagation. Anisotropy is described by 21 density-normalized elastic parameters naturally involved in asymptotic wave propagation in weakly anisotropic media. In a 1-D reference model, body-wave sensitivity to anisotropy is characterized by 'banana–doughnut' kernels which exhibit large, path-dependent variations and even sign changes. P -wave traveltimes appear much more sensitive to certain azimuthally anisotropic parameters than to the usual isotropic parameters, suggesting that isotropic P -wave tomography could be significantly biased by coherent anisotropic structures, such as slabs. Because of shear-wave splitting, the common cross-correlation traveltime anomaly is not an appropriate observable for S waves propagating in anisotropic media. We propose two new observables for shear waves. The first observable is a generalized cross-correlation traveltime anomaly, and the second a generalized 'splitting intensity'. Like P waves, S waves analysed based upon these observables are generally sensitive to a large number of the 21 anisotropic parameters and show significant path-dependent variations. The specific path-geometry of SKS waves results in favourable properties for imaging based upon the splitting intensity, because it is sensitive to a smaller number of anisotropic parameters, and the region which is sampled is mainly limited to the upper mantle beneath the receiver.  相似文献   

3.
In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.  相似文献   

4.
We investigate the reconstruction of a conductive target using crosswell time-domain electromagnetic tomography in the diffusive limit. The work is a natural extension of our ongoing research in the modification of time-domain methods for the rugged marine mid-ocean-ridge environment, an environment characterized by extreme topography and pronounced variations in crustal conductivity on all scales. We have proved both in theory and in practice that 'traveltime', the time taken for an electromagnetic signal to be identified at a receiver following a change of current in the transmitter, is an excellent, robust estimator of average conductivity on a path between transmitter and receiver. A simple estimate of the traveltime for a parallel electric dipole-dipole system is the time at which the derivative of the electric field with respect to logarithmic time at the receiver reaches its maximum. We have derived the fundamental relationship between the traveltime and the conductivity of the medium for a uniform whole-space. We have applied the concept of the traveltime inversion to the related crosswell problem and demonstrated reconstructions of finite targets based on tomographic analyses. Results show that the crosswell time-domain electromagnetic tomography can supply useful information, such as the location and shape of a conductive target.  相似文献   

5.
A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area.  相似文献   

6.
Velocity estimation remains one of the main problems when imaging the subsurface with seismic reflection data. Traveltime inversion enables us to obtain large-scale structures of the velocity field and the position of seismic reflectors. However, as the media currently under study are becoming more and more complex, we need to know the finer-scale structures. The problem is that below a certain range of velocity heterogeneities, deterministic methods become difficult to use, so we turn to a probabilistic approach. With this in view, we characterize the velocity field as a random field defined by its first and second statistical moments. Usually, a seismic random medium is defined as a homogeneous velocity background perturbed by a small random field that is assumed to be stationary. Thus, we make a link between such a random velocity medium (together with a simple reflector) and seismic reflection traveltimes. Assuming that the traveltimes are ergodic, we use 2-D seismic reflection geometry to study the decrease in the statistical traveltime fluctuations as a function of the offset (the source–receiver distance). Our formulae are based on the Rytov approximation and the parabolic approximation for acoustic waves. The validity and the limits are established for both of these approximations in statistically anisotropic random media. Finally, theoretical inversion procedures are developed for the horizontal correlation structure of the velocity heterogeneities for the simplest case of a horizontal reflector. Synthetic seismograms are then computed (on particular realizations of random media) by simulating scalar wave propagation via finite difference algorithms. There is good agreement between the theoretical and experimental results.  相似文献   

7.
We present a technique based on the single-scattering approximation that relates time-lapse localized changes in the propagation velocity to changes in the traveltime of singly scattered waves. We describe wave propagation in a random medium with homogeneous statistical properties as a single-scattering process where the fluctuations of the velocity with respect to the background velocity are assumed to be weak. This corresponds to one of two end-member regimes of wave propagation in a random medium, the first being single scattering, and the second multiple scattering. We present a formulation that relates the change in the traveltime of the scattered waves to a localized change in the propagation velocity by means of the Born approximation for the scattered wavefield. We validate the methodology with synthetic seismograms calculated with finite differences for 2-D acoustic waves. Potential applications of this technique include non-destructive evaluation of heterogeneous materials and time-lapse monitoring of heterogeneous reservoirs.  相似文献   

8.
Amplitude measurements of the transverse component of SKS waves, the so-called splitting intensity, can be used to formulate a non-linear inverse problem to image the 3-D variations of upper mantle anisotropy. Assuming transverse isotropy (or hexagonal symmetry), one can parametrize anisotropy by two anisotropic parameters and two angles describing the orientation of the symmetry axis. These can also be written as two collinear pseudo-vectors. The tomographic process consists of retrieving the spatial distribution of these pseudo-vectors, and thus resembles surface wave vectorial tomography. Spatial resolution results from the sensitivity of low-frequency SKS waves to seismic anisotropy off the ray path. The expressions for the 3-D sensitivity kernels for splitting intensity are derived, including the near-field contributions, and validated by comparison with a full wave equation solution based upon the finite element method. These sensitivity kernels are valid for any orientation of the symmetry axis, and thus generalize previous results that were only valid for a horizontal symmetry axis. It is shown that both lateral and vertical subwavelength variations of anisotropy can be retrieved with a dense array of broad-band stations, even in the case of vertically propagating SKS waves.  相似文献   

9.
Volcano seismology often deals with rather shallow seismic sources and seismic stations deployed in their near field. The complex stratigraphy on volcanoes and near‐field source effects have a strong impact on the seismic wavefield, complicating the interpretation techniques that are usually employed in earthquake seismology. In addition, as most volcanoes have a pronounced topography, the interference of the seismic wavefield with the stress‐free surface results in severe waveform perturbations that affect seismic interpretation methods. In this study we deal predominantly with the surface effects, but take into account the impact of a typical volcano stratigraphy as well as near‐field source effects. We derive a correction term for plane seismic waves and a plane‐free surface such that for smooth topographies the effect of the free surface can be totally removed. Seismo‐volcanic sources radiate energy in a broad frequency range with a correspondingly wide range of different Fresnel zones. A 2‐D boundary element method is employed to study how the size of the Fresnel zone is dependent on source depth, dominant wavelength and topography in order to estimate the limits of the plane wave approximation. This approximation remains valid if the dominant wavelength does not exceed twice the source depth. Further aspects of this study concern particle motion analysis to locate point sources and the influence of the stratigraphy on particle motions. Furthermore, the deployment strategy of seismic instruments on volcanoes, as well as the direct interpretation of the broad‐band waveforms in terms of pressure fluctuations in the volcanic plumbing system, are discussed.  相似文献   

10.
Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 'summary' traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km2 area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.  相似文献   

11.
A simple modification of the waveform inversion formula, based on the normal mode perturbation theory, is shown to lead to a formula for traveltime anomalies. The kernel which is derived can be used for traveltime inversion with automatic inclusion of finite frequency effects. Inversion for Earth structure with such kernels will lead to better resolution estimates than ray-theoretical traveltime inversion. Examples of kernels for transverse component seismograms are shown for direct S waves, ScS , Love waves and diffracted S waves. A measure of finite frequency effects is also proposed by comparing our formula with the one from ray theory. A quantity which should be 1 in the case of ray theory is computed for the finite frequency kernels and is shown to have deviations up to about 30 per cent from 1. Therefore, the use of ray theory for long-period body waves applies incorrect weight along a ray path and may introduce a small bias to an earth model.  相似文献   

12.
Anisotropy in multi-offset deep-crustal seismic experiments   总被引：1，自引：0，他引：1  
Modelling of deep-seismic wide-angle data commonly assumes that the Earth is heterogeneous and isotropic. It is important to know the magnitudes of errors that may be introduced by isotropic-based wide-angle models when the Earth is anisotropic. It is equally important to find ways of detecting anisotropy and determining its properties.
  This paper explores the errors introduced by interpreting anisotropic seismic data with isotropic models. Errors in P -wave reflector depths are dependent on the magnitude of the velocity anisotropy and the direction of the fast axis. The interpreted, isotropic, model velocity function is found to correspond closely to the horizontal velocity of the anisotropic medium. An additional observed parameter is the time mismatch , which we define to be the difference between the vertical two-way traveltime to a reflector and the time-converted wide-angle position of the reflector. The magnitude of the time mismatch is typically <1.0  s (when the whole crust is anisotropic) and is found to be closely related to the magnitude and sign of the anisotropic anellipticity. The relationships are extendible to more complicated models, including those with vertical velocity gradients, crustal zonation, and lower symmetry orders.
  A time mismatch may be symptomatic of the presence of anisotropy. We illustrate the observation of a time mismatch for a real multi-offset seismic data set collected north of Scotland and discuss the implications for crustal anisotropy in that region.  相似文献   

13.
  相似文献   

14.
Summary. Factors influencing the seismic response characteristics of marine profiling systems are reviewed. The single frequency case is used to illustrate the influence of different frequencies on the response, as well as the towing depths of the source and receiver, and the geometry of a linear receiving array. The more realistic case of band-limited source waveforms is considered, using frequency spectra calculated from theoretically derived airgun signals. The results show that the number and shape of sidelobes of the profiling system response, as well as the filtering characteristics for reflections arising from reflectors in the vertical plane perpendicular to the axis of the receiver array are determined by the depths of the source and receiver and the relative amplitudes of the frequencies in the source waveform. These factors, along with the configuration of the hydrophone elements in the receiver array, determine the frequency and amplitude attenuation of reflections in the vertical plane containing the receiver array.
The filtering characteristics of the system both in and out of the vertical plane containing the receiver array are discussed, with implications for discriminating between off-axis and in-plane reflections. A plan view of the response of the system is constructed in the time domain for various profiling configurations and sources of different frequency content at a given time. This example shows how useful the resulting pictures are for optimizing acquisition parameters in profiling experiments.  相似文献   

15.
We present the results of Rayleigh wave and Love wave phase velocity tomography in the western United States using ambient seismic noise observed at over 250 broad-band stations from the EarthScope/USArray Transportable Array and regional networks. All available three-component time-series for the 12-month span between 2005 November 1 and 2006 October 31 have been cross-correlated to yield estimated empirical Rayleigh and Love wave Green's functions. The Love wave signals were observed with higher average signal-to-noise ratio (SNR) than Rayleigh wave signals and hence cannot be fully explained by the scattering of Rayleigh waves. Phase velocity dispersion curves for both Rayleigh and Love waves between 5 and 40 speriod were measured for each interstation path by applying frequency–time analysis. The average uncertainty and systematic bias of the measurements are estimated using a method based on analysing thousands of nearly linearly aligned station-triplets. We find that empirical Green's functions can be estimated accurately from the negative time derivative of the symmetric component ambient noise cross-correlation without explicit knowledge of the source distribution. The average traveltime uncertainty is less than 1 s at periods shorter than 24 s. We present Rayleigh and Love wave phase speed maps at periods of 8, 12, 16,and 20 s. The maps show clear correlations with major geological structures and qualitative agreement with previous results based on Rayleigh wave group speeds.  相似文献   

16.
The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.  相似文献   

17.
Large scale seismic anisotropy in the Earth's mantle is likely dynamically supported by the mantle's deformation; therefore, tomographic imaging of 3-D anisotropic mantle seismic velocity structure is an important tool to understand the dynamics of the mantle. While many previous studies have focused on special cases of symmetry of the elastic properties, it would be desirable for evaluation of dynamic models to allow more general axis orientation. In this study, we derive 3-D finite-frequency surface wave sensitivity kernels based on the Born approximation using a general expression for a hexagonal medium with an arbitrarily oriented symmetry axis. This results in kernels for two isotropic elastic coefficients, three coefficients that define the strength of anisotropy, and two angles that define the symmetry axis. The particular parametrization is chosen to allow for a physically meaningful method for reducing the number of parameters considered in an inversion, while allowing for straightforward integration with existing approaches for modelling body wave splitting intensity measurements. Example kernels calculated with this method reveal physical interpretations of how surface waveforms are affected by 3-D velocity perturbations, while also demonstrating the non-linearity of the problem as a function of symmetry axis orientation. The expressions are numerically validated using the spectral element method. While challenges remain in determining the best inversion scheme to appropriately handle the non-linearity, the approach derived here has great promise in allowing large scale models with resolution of both the strength and orientation of anisotropy.  相似文献   

18.
Multicomponent near-surface correction for land VSP data   总被引：1，自引：0，他引：1  
Multicomponent seismic data collected using directional sources are degraded by the wave excitation process due to inaccurate control of the ground motion. unequal activation strengths or ground couplings between differently oriented sources, and misalignment of the pad. These acquisition uncertainties are exacerbated by the complicated near-surface scattering present in most seismic areas. Neither group of effects should be neglected in multicomponent analyses that make use of relative wavefield attributes derived from compressional and shear waves. These effects prevent analysis of the direct and reflected waves using procedures based on standard scalar techniques or a prima facia interpretation of the vector wavefield properties, even for the seemingly straightforward case of a near-offset vertical seismic profile (VSP). Near-surface correction, using a simple matrix operator designed from the shallowest recordings, alleviates many of these interpretational difficulties in near-offset VSP data. Results from application of this technique to direct waves from a nine-component VSP shot at the Conoco test-site facility, Oklahoma, are encouraging. The technique corrects for unexpected compressional-wave energy from shear-wave vibrators and collapses near-surface multiples, thus facilitating further processing for the upgoing wavefield. The method provides a simple and effective processing step for routine application to near-offset VSP analyses.  相似文献   

19.
We develop a systematic approach to the phase identification of late-arriving groups in 2-D seismic data. Waveforms in the same traveltime branch are grouped, and synthetic traveltimes for all phases are calculated using an initial approximation to the 2-D structure. For each group, we identify the two synthetic phases providing the smallest RMS residuals. If their ratio is less than some predetermined threshold, then the group's phase is ambiguous and both assignments must be tested by traveltime inversion. If there are n unidentified groups, we construct 2 n phase tables and perform a traveltime inversion on every plausible phase assignment. The phase table that provides the highest value of the posterior probability density is taken as correct, and a 2-D velocity model is constructed from the data. This approach is shown to be effective and efficient on both simulated and real data. In addition, the residuals associated with late-arriving groups provide a means of identifying deficiencies in the initial model.  相似文献   

20.
We invert differential SdS-SS traveltime residuals measured from stacked waveforms and finite-frequency sensitivity kernels for topography on the 410- and 660-km discontinuities. This approach yields higher resolution images of transition zone thickness than previous stacking methods, which simply average/smooth over topographic features. Apparent structure measured using simple stacking is highly dependent upon the bin size of each stack. By inverting for discontinuity topography with a variety of bin sizes, we can more accurately calculate the true structure. The inverted transition zone model is similar to simple stack models with an average thickness of 242 km, but the lateral variations in thickness are larger in amplitude and smaller in scale. Fast seismic velocities in 3-D mantle models such as SB4L18 correlate with areas of thicker transition zone. The elongated curvilinear regions of thickened transition zone that occur near subduction zones are narrow and high amplitude, which suggests relatively little lateral spreading and warming of subducted lithosphere within the transition zone. The anomalously thin transition zone regions are laterally narrow, and not broadly continuous. If these variations in transition zone thickness are interpreted as thermal in nature, then this model suggests significant temperature variations on small lateral scales.  相似文献   

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

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