首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In this paper, attenuation values are obtained from seismic data using instantaneous-frequency matching and spectral ratios. to obtain differential t * values using instantaneous-frequency matching, a near offset reference pulse is attenuated until the resulting instantaneous frequency matches the observed value at the receiver. Prior to matching, filtering can be applied to each trace in order to reduce the effects of noise on the calculated instantaneous frequencies. In the second method, the spectral ratio between a receiver pulse and a reference pulse is used to obtain differential t * values. to obtain an unbiased estimate, a variable spectral bandwidth is used depending on the noise level of the data. the two methods are tested using synthetic traces and then applied to crustal refraction data from the 1986 PASSCAL Ouachita experiment. Results show that the differential t * values obtained using filtered, instantaneous-frequency matching are consistent with and have less scatter than those obtained from spectral ratios with a variable bandwidth.  相似文献   

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

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

4.
We present velocity constraints for the upper-mantle transition zones beneath Central Siberia based on observations of the 1982 RIFT Deep Seismic Sounding (DSS) profile. The data consist of seismic recordings of a nuclear explosion in north-western Siberia along a 2600 km long seismic profile extending from the Yamal Peninsula to Lake Baikal. We invert seismic data from the mantle transition zones using a non-linear inversion scheme using a genetic algorithm for optimization and the WKBJ method to compute the synthetic seismograms. A statistical error analysis using a graph-binning technique was performed to provide uncertainty values in the velocity models.
Our best model for the upper-mantle velocity discontinuity near 410 km depth has a two-stage velocity-gradient structure, with velocities increasing from 8.70–9.25 km s−1 over a depth range of 400–415 km, a gradient of 0.0433 s−1, and from 9.25–9.60 km s−1 over a depth range of 415–435 km, a gradient of 0.0175 s−1. This derived model is consistent with other seismological observations and mineral-physics models. The model for the velocity discontinuity near 660 km depth is simple, sharp and includes velocities increasing from 10.15 km s−1 at 655 km depth to 10.70 km s−1 at 660 km depth, a gradient of 0.055 s−1.  相似文献   

5.
About 50 000 P and S arrival times and 25 000 values of t * recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P , S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp , Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting.  相似文献   

6.
The frequency-domain version of waveform tomography enables the use of distinct frequency components to adequately reconstruct the subsurface velocity field, and thereby dramatically reduces the input data quantity required for the inversion process. It makes waveform tomography a computationally tractable problem for production uses, but its applicability to real seismic data particularly in the petroleum exploration and development scale needs to be examined. As real data are often band limited with missing low frequencies, a good starting model is necessary for waveform tomography, to fill in the gap of low frequencies before the inversion of available frequencies. In the inversion stage, a group of frequencies should be used simultaneously at each iteration, to suppress the effect of data noise in the frequency domain. Meanwhile, a smoothness constraint on the model must be used in the inversion, to cope the effect of data noise, the effect of non-linearity of the problem, and the effect of strong sensitivities of short wavelength model variations. In this paper we use frequency-domain waveform tomography to provide quantitative velocity images of a crosshole target between boreholes 300 m apart. Due to the complexity of the local geology the velocity variations were extreme (between 3000 and 5500 m s−1), making the inversion problem highly non-linear. Nevertheless, the waveform tomography results correlate well with borehole logs, and provide realistic geological information that can be tracked between the boreholes with confidence.  相似文献   

7.
Effects of fractures on seismic-wave velocity and attenuation   总被引：1，自引：0，他引：1  
The effects of fractures on the seismic velocity and attenuation of a rock are investigated using theoretical results and experimental data. Fractures in a rock mass influence the traveltimes and amplitudes of seismic waves that have propagated through them. The displacement discontinuity model, recently employed in fracture investigations, is modified to describe the effect of fractures on seismic-wave velocity and attenuation. This new model, the modified displacement discontinuity model (MDD), is formulated in a way analogous to transmission-line analysis. The fractures are treated as transmission lines for the passage of seismic waves. The MDD takes into consideration realistic fracture parameters which include the fracture length, the fractional area of a fracture surface in contact, and the nature of the infilling material. A single fracture of varying geometric and material properties is shown to affect dramatically the transmission properties of a propagating waveform, and hence the seismic velocity and attenuation. These effects have been shown to result in a frequency-dependent velocity and attenuation. The sensitivity of the fracture parameters to seismic-wave velocity and attenuation was investigated and interesting results were obtained. Fracture parameters used in designing experimental models consisting of synthetically manufactured cracks were fed into the MDD and a well-known crack model, Hudson's model, for comparison. Velocities as a function of the incident-wave angle were obtained from both numerical models and were compared with the results from the experimental modelling. For P waves, the MDD model results show better agreement with those of the experimental model for all crack densities investigated than those from Hudson's model.  相似文献   

8.
Investigation of teleseismic P -wave recordings at a temporary network in the Eastern Carpathians, equipped with predominantly short-period sensors, is compared with synthetic modelling of anelastic attenuation of teleseismic waves in the upper mantle. Using the t * approach, we examine variations of amplitude decrease over frequency for teleseismic recordings in the frequency band 0.5–1.5 Hz. The results reveal a consistent pattern of increased t * values in the centre of the network, in the Vrancea region at the bend of the Carpathian Arc, although the magnitude of the observed variation in t * is much higher than expected. Synthetic t * parameter computations for the same event-receiver configurations reproduce the observed pattern in terms of relative variations. However, the amplitude of the synthetic t * values explains only 10–20 per cent of the observed variation in t *. t * is not a direct measure for anelastic attenuation but rather for a combination of anelastic and other attenuating effects such as scattering and amplitude fluctuation related to velocity inhomogeneities. If regional amplitude variations are solely attributed to anelastic attenuation, all other effects are mapped into Q . We discuss the role of anelastic attenuation and other effects in the case of the Eastern Carpathians and conclude that t * is an unsuitable parameter to characterize anelastic attenuation in the Eastern Carpathians.  相似文献   

9.
Summary. In this paper we discuss some aspects of estimating t * from short-period body waves and present some limits on t* (f) models for the central and south-western United States (CUS and SWUS). We find that for short-period data, with frequencies above 1 or 2 Hz, while the average spectral shape is stable, the smaller details of the spectra are not; thus, only an average t *, and not a frequency-dependent t *, can be derived from such information. Also, amplitudes are extremely variable for short-period data, and thus a great deal of data from many stations and azimuths must be used when amplitudes are included in attenuation studies.
The predictions of three pairs of models for t* (f) in the central and south-western United States are compared with time domain observations of amplitudes and waveforms and frequency domain observations of spectral slopes to put bounds on the attenuation under the different parts of the country. A model with the t * values of the CUS and SWUS converging at low frequencies and differing slightly at high frequencies matches the spectral domain characteristics, but not the time domain amplitudes and waveforms of short-period body waves. A model with t * curves converging at low frequencies, but diverging strongly at high frequencies, matches the time domain observations, but not the spectral shapes. A model with nearly-parallel t* (f) curves for the central and south-western United States satisfies both the time and frequency domain observations.
We conclude that use of both time and frequency domain information is essential in determining t* (f) models. For the central and south-western United States, a model with nearly-parallel t* (f) curves, where Δ t *∼ 0.2 s, satisfies both kinds of data in the 0.3–2 Hz frequency range.  相似文献   

10.
The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic.  相似文献   

11.
We report results from the Seismic Wide-Angle and Broadband Survey carried out over the Mid North Sea High. This paper focuses on integrating the information from a conventional deep multichannel reflection profile and a coincident wide-angle profile obtained by recording the same shots on a set of ocean bottom hydrophones (OBH). To achieve this integration, a new traveltime inversion scheme was developed (reported elsewhere) that was used to invert traveltime information from both the wide-angle OBH records and the reflection profile simultaneously. Results from the inversion were evaluated by producing synthetic seismograms from the final inversion model and comparing them with the observed wide-angle data, and an excellent match was obtained. It was possible to fine-tune velocities in less well-resolved parts of the model by considering the critical distance for the Moho reflection. The seismic velocity model was checked for compatibility with the gravity field, and used to migrate and depth-convert the reflection profile. The unreflective upper crust is characterized by a high velocity gradient, whilst the highly reflective lower crust is associated with a low velocity gradient. At the base of the crust there are several subhorizontal reflectors, a few kilometres apart in depth, and correlatable laterally for several tens of kilometres. These reflectors are interpreted as representing a strike section through northward-dipping reflectors at the base of the crust, identified on orthogonal profiles by Freeman et al. (1988) as being slivers of subducted and imbricated oceanic crust, relics of the mid-Palaeozoic Iapetus Ocean.  相似文献   

12.
A large data set of amplitude measurements of minor and major arc Rayleigh waves in the period range 73–171 s is collected. By comparing these amplitudes with the amplitudes of synthetic waveforms calculated by mode summation, maps of lateral variations in the apparent attenuation structure of the Earth are constructed. An existing formalism for predicting the effects of focusing is employed to calculate amplitude perturbations for the same data set. These perturbations are used to construct 'pseudo‐attenuation' maps and these results are compared with the apparent attenuation maps calculated from the data. It is shown that variations in Rayleigh wave amplitude perturbations in the Earth are dominated by attenuation at long wavelengths (below about degree 8) and by elastic structure at shorter wavelengths. It is also shown that the linear approximation for focusing is successful at predicting Rayleigh wave amplitudes using existing phase velocity maps. These results indicate that future attempts to model the velocity structure of the Earth would be assisted by incorporating amplitude data and by jointly inverting for Q structure.  相似文献   

13.
Summary. The general problem of inverting Love-wave dispersion and amplitude data to obtain a velocity and Qs structure is considered. A formulation is used which incorporates attenuation into the Haskell-Thompson matrix method in an exact manner and thus retains the inherent non-linearity in the anelasticity. The resulting exact inversion kernels allow simultaneous inversion for velocity and intrinsic attenuation parameters. The method is applied to synthetic data which allows a comparison to be made with inexact kernels. The results indicate that the use of inexact kernels may introduce spurious oscillations into the Qs structure and that a simultaneous inversion can be more stable than inverting for velocity alone.  相似文献   

14.
We present new methods for the interpretation of 3-D seismic wide-angle reflection and refraction data with application to data acquired during the experiments CELEBRATION, 2000 and ALP 2002 in the area of the Eastern Alps and their transition to the surrounding tectonic provinces (Bohemian Massif, Carpathians, Pannonian domain, Dinarides). Data was acquired on a net of arbitrarily oriented seismic lines by simultaneous recording on all lines of seismic waves from the shots, which allows 2-D and 3-D interpretations. Much (80%) of the data set consists of crossline traces. Low signal to noise (S/N) ratio in the area of the young orogens decreases the quality of travel time picks. In these seismically heterogeneous areas it is difficult to assign clearly defined arrivals to the seismic phases, in particular on crossline record sections.
In order to enhance the S/N ratio, signal detection and stacking techniques have been applied to enhance the Pg -, Pn - and PmP phases. Further, inversion methods have been developed for the interpretation of WAR/R-data, based on automated 1-D inversion ( Pg ) and the application of the delay time concept ( Pn ). The results include a 3-D velocity model of the crust based on Pg waves, time and depth maps of the Moho and a Pn -velocity map. The models based on stacked data are robust and provide a larger coverage, than models based on travel time picks from single-fold (unstacked) traces, but have relatively low resolution, especially near the surface. They were used as the basis for constructing models with improved resolution by the inversion of picks from single-fold data. The results correlate well with geological structures and show new prominent features in the Eastern Alps area and their surrounds. The velocity distribution in the crust has strong lateral variations and the Moho in the investigation area appears to be fragmented into three parts.  相似文献   

15.
  相似文献   

16.
Summary. The structure of the upper lithosphere beneath southern Germany, northern Switzerland and west-central Utah (U.S.A.) has been investigated in detail by various geophysical methods. A synoptic interpretation of travel time and amplitude data obtained in seismic refraction and wide-angle reflection surveys, combined with near-normal incidence reflection observations, now permits the elucidation of the fine structure in a more quantitative and unified manner. With this scheme it is possible to unambiguously identify low-velocity zones and to deduce velocity gradients if reliable amplitude information is included in the inversion process.  相似文献   

17.
A crustal seismic velocity model for the UK, Ireland and surrounding seas   总被引：1，自引：0，他引：1  
A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.  相似文献   

18.
Summary. Studies of teleseismic P -and S -wave amplitudes and spectra in the 0.5–4 Hz band show large variations in the attenuative properties of the upper mantle under the United States. The data indicate that attenuation is greatest under the south-western United States including, but not confined to, the Basin and Range province. The lowest attenuation prevails under the north central shield regions. The north-eastern part of the country, consisting of New England and possibly including a larger area along the eastern seaboard, is characterized by moderate attenuation in the mantle.
The level of the high-frequency energy in short-period seismic waves and the differences between Q values derived from short-and long-period data indicate that Q is frequency dependent. The form of frequency dependence of t * compatible with the data in the 0.5–4 Hz range does not allow a rapid decrease of t * with increasing frequency. Rather it supports a gradual decrease covering the broader 0.1–4 Hz range. The curves of t * versus frequency, for shield-to-shield and mixed shield-to-western United States type paths are parallel with an average difference of 0.2 s in t * in the short-period band, but may diverge towards the long-period band. For both curves t *p is below 1 s. For shield-to-shield paths t *p must be below 0.5 s at 1 Hz.  相似文献   

19.
In an accompanying paper, we used waveform tomography to obtain a velocity model between two boreholes from a real crosshole seismic experiment. As for all inversions of geophysical data, it is important to make an assessment of the final model, to determine which parts of the model are well-resolved and can confidently be used for geological interpretation. In this paper we use checkerboard tests to provide a quantitative estimate of the performance of the inversion and the reliability of the final velocity model. We use the output from the checkerboard tests to determine resolvability across the velocity model. Such tests can act as good guides for designing appropriate inversion strategies. Here we discovered that, by including both reference-model and smoothing constraints in initial inversions, and then relaxing the smoothing constraint for later inversions, an optimum velocity image was obtained. Additionally, we noticed that the performance of the inversion was dependent on a relationship between velocity perturbation and checkerboard grid-size: larger velocity perturbations were better-resolved when the grid-size was also increased. Our results suggest that model assessment is an essential step prior to interpreting features in waveform tomographic images.  相似文献   

20.
Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= tu/Qm where t u is the travel time and Qm the quality factor in the absorption band, and τ−1 m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ m. For very small τ m , this model is shown to be very similar to previously proposed attenuation models. As τ m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.  相似文献   

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

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