首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 546 毫秒
1.
This paper provides a review of array-based imaging techniques that use converted and scattered teleseismic waves. It addresses various aspects of the imaging process, from the preprocessing of the data to the application of the imaging algorithms. The reviewed techniques form a continuum with respect to the level of complexity adopted in the treatment of the scattering problem. On one end of the spectrum, images may be produced by simple stacking of normalized P-to-S conversion records (i.e., receiver functions), which are binned according to station or common conversion points (CCP) and mapped to depth. Finer resolution can be achieved through the stacking of singly scattered wavefields along diffraction hyperbolae to recover relative scattering intensity/potential at individual points through a 2-D or 3-D model space. Moving to higher levels of complexity, we find methods that involve inversion/backprojection of scattered teleseismic wavefield to recover estimates of localized material property perturbations with respect to an a priori background model.  相似文献   

2.
Over the past decade, the H-κ stacking technique of Zhu and Kanamori (J Geophys Res 105:2969–2980, 2000) has become a standard tool to determine the crustal thickness H and the bulk crustal vP/vS ratio κ from teleseismic receiver functions. It is obvious that unfavorable noise conditions as well as a complex 3D velocity structure can severely hamper the interpretation of receiver-function data. However, we observe that ambiguities can even arise from a simple 1D layered velocity structure which raises a high potential for misinterpretations. To analyze the feasibility and basic limitations of the H-κ stacking method, we conduct a series of tests based on synthetic data. The impact of different given elementary parameters, related either to the velocity structure or to the data processing, is evaluated in a series of eight individual tests. We deliberately exclude complications such as 3D structural variations and/or noise to show that even a simple 1D velocity structure, involving, e.g., an additional inter-crustal discontinuity, can have significant consequences for the interpretation of the results. However, our modeling suggests that more complex crustal structures may lead to even less reliable results. Additionally, our tests illustrate that time shifts of the maxima in the H-κ domain due to the superposition and merging of individual phases can lead to significantly overestimated vP/vS ratios. In general, the depth to the Moho (or other discontinuities of interest) is less significantly affected. Our tests indicate the necessity to accurately check delay times derived from the maxima of the H-κ stacks against corresponding phases in the receiver functions. Repeating the stacking with varied weighting factors and filter ranges can help to reduce the ambiguities and to avoid possible misinterpretation.  相似文献   

3.
Common-reflection-surface (CRS) stack for common offset   总被引：8，自引：0，他引：8  
We provide a data-driven macro-model-independent stacking technique that migrates 2D prestack multicoverage data into a common-offset (CO) section. We call this new process the CO common-reflection-surface (CRS) stack. It can be viewed as the generalization of the zero-offset (ZO) CRS stack, by which 2D multicoverage data are stacked into a well-simulated ZO section. The CO CRS stack formula can be tailored to stack P-P, S-S reflections as well as P-S or S-P converted reflections. We point out some potential applications of the five kinematic data-derived attributes obtained by the CO CRS stack for each stack value. These include (i) the determination of the geometrical spreading factor for reflections, which plays an important role in the construction of the true-amplitude CO section, and (ii) the separation of the diffractions from reflection events. As a by-product of formulating the CO CRS stack formula, we have also derived a formula to perform a data-driven prestack time migration.  相似文献   

4.
The P receiver function includes P-to-SV converted phases and multiple reverberations of the discontinuities in the crust and mantle. The time of these phases is related to the crustal thickness and vP/vS ratio, and the amplitude of these phases is mainly controlled by the velocity and density contrast of interfaces. By using H-κ stacking method, this work estimated the crustal thickness and vP/vS ratio beneath the stations in the Guangdong province of South China. The velocity and density contrast (δβ-δρ) scanning stacking algorithm of the receiver function is applied to constrain the velocity and density contrast of the Moho in Guangdong province. This work analyzed the results of the crustal thickness, vP/vS ratio, and the velocity and density contrasts of Moho. The results indicate that the velocity contrast is higher beneath Yangjiang area in western Guangdong province and Nanao area in eastern Guangdong, which has a strong correlation with the distribution of geothermal springs in local areas and the characteristics of high heat flow. The velocity contrast of Moho has also a good correlation with the vP/vS ratio and the crustal thickness, which indicates that there is a strong material composition contrasts of the Moho in the study area. Velocity and density contrasts of Moho in some local area (such as western Guangdong) are somewhat consistent with the seismic activities.  相似文献   

5.
—Shear (S) waves differ from compressional (P) waves because of their lower propagation velocities, their lower frequencies and due to the different character of their particle motion. The move-out of travel-time branches of S-wave reflections is different from P waves owing to the difference in the propagation velocities. To distinguish between P and S waves requires broadband-frequency acquisition, long receiver arrays and three-component recording. S-wave generation at the source and P-to-S-wave conversion at crustal interfaces can be very efficient, implying that there is a real danger of misinterpreting signals if only vertical components are used. On the other hand, integrated P- and S-wave studies promise to provide very efficient lithological discriminators in the crystalline crust, in particular concerning the quartz content, and indicators for rock anisotropy, which can be interpreted for the existence of fine layering, the direction of the recent stress regime (alignments of micro-fractures) or for the direction of palaeo-stress (alignments of minerals).  相似文献   

6.
The North China Craton (NCC) is one of the oldest cratons on earth. Several important tectonic transformations of Mesozoic-Cenozoic tectonic regime led to the destruction of the North China craton. The knowledge of crustal structure can provide important constraints for the formation and evolution of cratons. New maps of sediment thickness, crustal thickness (H) and vP/vS (κ) in the central and western NCC were obtained using sequential H-κ stacking. P-wave receiver functions are calculated using teleseismic waveform data recorded by 405 stations from ChinArray project. Benefiting from the densely distribution of temporary seismic stations, our results reveal details of the crustal structure in the study area. The thickness of sedimentary layer in North China ranges from 0–6.4 km, and the thickest sedimentary layer is in Ordos block and its surroundings (about 2.8–6 km); The thickness of sedimentary layer in the Mongolia fold belt and Yinshan orogenic belt is relatively thin (less than 1 km). The crustal thickness of the study area varies between 27–48 km, of which the crust of the North China Plain is about 30–33 km, the central NCC is about 33–40 km, and the Ordos block is 40–48 km thick. The average vP/vS ratios in the study area is mostly between 1.66 and 1.90, and that in the Yanshan-Taihang mountain fold belt is between 1.70 and 1.85, and that in the Ordos block is between 1.65 and 1.90, with an average value of 1.77, indicating the absence of a thick basaltic lower crust. The obvious negative correlation between crustal thickness and average vP/vS ratio within Ordos and Central Asia orogenic belt may be related to magmatic underplating during the crustal formation. There is no significant correlation between the crustal thickness and the vP/vS ratio in the Lüliang-Taihang mountain fold belt, which may be related to the multiple geological processes such as underplating and crustal extension and thinning in this area. The lack of correlation between crust thickness and topography in the central orogenic belt and the North China Basin indicates the topography of these areas are controlled not only by crustal isostatic adjustment but also by the lithospheric mantle processes.  相似文献   

7.
—Reflectivity synthetic seismograms demonstrate that the type, layering and orientation of 1-D anisotropy influences strongly the coda of teleseismic P waves at periods T > 1 sec, particularly P-SH converted waves. We assume the simplest form of anisotropy described by an elastic tensor with a symmetry axis ? of arbitrary orientation. The resulting phase velocities vary as cos 2ξ with respect to that axis. Using three families of simple crustal models, we compare the effects of an anisotropic surface layer with reverberations caused by both "thick" and "thin" layers of anisotropy at depth. If anisotropy in the surface layer is significant, the polarization of direct P can be distorted to generate a transverse component, followed by Ps and a prominent shear reverberation converted from direct P at the free surface. If the anisotropic layer is buried, the first, and often the most prominent, arrival on the transverse component is the P-to-SH conversion at its upper surface. If the anisotropic layer is sufficiently thin, P-to-SH conversions from its boundaries interfere to form a derivative pulse shape on the transverse component, which could be mistaken as the signature of shear-wave splitting. If ? is horizontal, compressional (P) and shear (S) anisotropy both produce similar waveform perturbations with four-lobed azimuthal patterns, suggesting that a weighted stack of P coda from different back-azimuths would improve signal-to-noise. For ? tilted between the horizontal and vertical, however, the effects of P- and S-anisotropy differ greatly. The influence of P-anisotropy on P-to-S conversion is greatest for a symmetry axis tilted at 45° to the vertical, where its azimuthal pattern has two lobes, rather than four. Combinations of P- and S-anisotropy typically lead to a composite azimuthal dependence in the P-coda reverberations.  相似文献   

8.
It is noticed that few geophysical studies have been carried out to decipher the crustal structure of southwestern part of the Northeast India comprising of Tripura fold belt and Bengal basin as compared to the Shillong plateau and the Brahmaputra basin. This region has a long history of seismicity that is still continuing. We have determined first-order crustal features in terms of Moho depths (H) and average VP/VS ratios (κ) using H-κ stacking technique. The inversion of receiver functions data yields near surface thick sedimentary layer in the Bengal basin, which is nearly absent in the Shillong plateau and Tripura fold belt. Our result suggests that the crust is thicker (38–45 km) in the Tripura fold belt region with higher shear-wave velocity in the lower crust than the Shillong plateau. The distribution of VP/VS ratio indicates heterogeneity throughout the whole region. While low to medium value of Poisson’s ratio (1.69–1.75) indicates the presence of felsic crust in the Shillong plateau of the extended Indian Archean crust. The medium to high values of VP/VS ratio (> 1.780) in the Bengal basin and the Tripura fold belt region represent mafic crust during the formation of the Bengal delta and the Tripura fold belt creation in the Precambrian to the Permian age. The depth of the sediments in the Bengal basin is up to 8 km on its eastern margin, which get shallower toward its northeastern and southeastern margins.  相似文献   

9.
The crustal structure in Myanmar can provide valuable information for the eastern margin of the ongoing Indo-Eurasian collision system. We successively performed H–k stacking of the receiver function and joint inversion of the receiver function and surface wave dispersion to invert the crustal thickness (H), shear wave velocity (VS), and the VP/VS ratio (k) beneath nine permanent seismic stations in Myanmar. H was found to increase from 26 ?km in the south and east of the study area to 51 ?km in the north and west, and the VP/VS ratio was complex and high. Striking differences in the crust were observed for different tectonic areas. In the Indo-Burma Range, the thick crust (H ?～ ?51 ?km) and lower velocities may be related to the accretionary wedge from the Indian Plate. In the Central Myanmar Basin, the thin crust (H ?= ?26.9–35.5 ?km) and complex VP/VS ratio and VS suggest extensional tectonics. In the Eastern Shan Plateau, the relatively thick crust and normal VP/VS ratio are consistent with its location along the western edge of the rigid Sunda Block.  相似文献   

10.
A Seismic Model of Casing Failure in Oil Fields   总被引：1，自引：0，他引：1  
—We develop a seismic model that characterises the sudden tensional failure of oil-well casings. The energy released by the rupture of a well casing is transformed into heat and seismic energy. The upper bound of the seismic efficiency of this process is estimated at about 3%. The static situation at the completion of a casing failure episode is modelled by calculating the static displacement field generated by two opposing forces separated by an arm. The azimuthal patterns of these displacements and the change in the strain and stress fields caused by the force couple are described. The dynamics of the failure episode are modelled as a dipole with a seismic moment equivalent to the product of the average drop in shear stress, the failure surface, and an arm. The radiated P and S waves have mean-square radiation pattern coefficients of 1/5 for P waves and 2/15 for S waves. The displacement field as a function of time during rupture and the spectral properties in the far field are derived. The most promising seismic parameters that can be used for distinguishing between casing failure events and other possible events are polarisation properties of S waves and S/P amplitude ratios. S-wave polarisation distinguishes between shear events and casing failure events. S/P amplitude ratios distin guish between tensile events and casing failure events.  相似文献   

11.
Several parameters are needed to describe the converted-wave （C-wave） moveout in processing multi-component seismic data, because of asymmetric raypaths and anisotropy. As the number of parameters increases, the converted wave data processing and analysis becomes more complex. This paper develops a new moveout equation with two parameters for C-waves in vertical transverse isotropy （VTI） media. The two parameters are the C-wave stacking velocity （Vc2） and the squared velocity ratio （7v,i） between the horizontal P-wave velocity and C-wave stacking velocity. The new equation has fewer parameters, but retains the same applicability as previous ones. The applicability of the new equation and the accuracy of the parameter estimation are checked using model and real data. The form of the new equation is the same as that for layered isotropic media. The new equation can simplify the procedure for C-wave processing and parameter estimation in VTI media, and can be applied to real C-wave processing and interpretation. Accurate Vc2 and Yvti can be deduced from C-wave data alone using the double-scanning method, and the velocity ratio model is suitable for event matching between P- and C-wave data.  相似文献   

12.
Based on the experimental and calculated data, the model of depth variations in density, as well as the velocities of longitudinal (P) and transverse (S) waves, within the upper, middle, and lower crusts of the Kola-Norwegian block down to a depth of about 40 km is suggested. The variations in density and the velocities of the P- and S-waves are primarily caused by the changes in the mineral composition of the rocks. The relative reduction in the velocities of the P- and S-waves under the action of the increasing pressure and temperature in the depth interval from 5 to 37 km is estimated at ～2%.  相似文献   

13.
Molecular dynamics (MD) modeling of the 10-Å phase, Mg3Si4O10(OH)2·xH2O, with x=2/3, 1.0 and 2.0 shows complex structural changes with pressure, temperature and water content and provides new insight into the structures and stabilization of these phases under subduction zone conditions. The structure(s) of this phase and its role as a reservoir of water in the mantle have been controversial, and these calculations provide specific predictions that can be tested by in situ diffraction studies. At ambient conditions, the computed structures of talc (x=0) and the 10-Å phases with x=2/3 and 1.0 are stable over the 350-ps period of the MD simulations. Under these conditions, the 10-Å phases show phlogopite-like layer stacking in good agreement with previously published structures based on powder X-ray diffraction data for samples quenched from high-pressure and high-temperature experiments. The calculations show that the 10-Å phase with x=2.0 is unstable at ambient conditions. The computed structures at P=5.5 GPa and T=750 K, well within the known stability field of the 10-Å phase, change significantly with water content, reflecting changing H-bonding configurations. For x=2/3, the layer stacking is talc-like, and for x=1.0, it is phlogopite-like. The calculations show that transformation between these two stackings occurs readily, and that the talc-like stacking for the x=2/3 composition is unlikely to be quenchable to ambient conditions. For x=2.0, the layer stacking at P=5.5 GPa and T=750 K is different than any previously proposed structure for a 10-Å phase. In this structure, the neighboring basal oxygens of adjacent magnesium silicate layers are displaced by b/3 (about 3 Å) resulting in the Si atoms of one siloxane sheet being located above the center of the six-member ring across the interlayer. The water molecules are located 1.2 Å above the center of all six-member rings and accept H-bonds from the OH groups located below the rings. The b/3-displaced structure does not readily transform to either the talc-like or phlogopite-like structure, because neither of these stackings can accommodate two water molecules per formula unit. There is likely to be a compositional discontinuity and phase transition between the b/3-displaced phase and the phase with phlogopite-like stacking. The simulations reported here are the first to use the recently developed CLAYFF force field to calculate mineral structures at elevated pressures and temperatures.  相似文献   

14.
Coalbed methane can be detected employing the amplitudevariation-with-offset technique. However, there are two issues in applying this technique to a coalbed: strong azimuthal anisotropy resulting from high-density fractures, and the seismic response being composed of many or several individual reflections within the coalbed. To overcome these difficulties, we present an exact solution for reflections in extensive dilatancy anisotropy media. First, we build a three-layer model and simulate the wave propagation in this model. Then we derive an exact P- and converted S-wave reflection coefficient equation based on boundary conditions. Finally, substituting given model parameters into the exact equation, we obtain the variation in the reflection coefficient with incidence angle. The results show that the fracture factors, wavelet frequency and thickness of the coalbed have different effects on the reflection coefficient. Furthermore, we create a synthetic seismogram by forward calculation, and the result fits well with results of the exact equation.  相似文献   

15.
The primary objective of this study is for presenting some simple-to-use expressions relating the shear and dilatational wave velocities (VS and VP) to some physical and constitutive parameters of unsaturated soils. To this purpose, a simplified formulation is developed using the theory of linear poroelasticity in conjunction with some constitutive parameters widely used in geotechnical engineering. The derived expressions are of practical interest in view of the fact that they could be employed for evaluating the involved soil parameters from VS and VP measurements by in-situ or laboratory geophysical tests.  相似文献   

16.
—Anisotropy in the subcontinental lithosphere becomes increasingly important, because it is observed in many seismic studies especially for P n -waves. Typical rocks of the uppermost mantle are peridotites, which predominantly exhibit a pronounced elastic anisotropy. This anisotropy is mainly caused by the anisotropic elastic properties and the lattice preferred orientation (here referred to as texture) of olivine. To evaluate the elastic anisotropy of peridotites from the subcontinental lithosphere, specimens of the Northern Hessian Depression (Germany) and the Balmuccia Ultramafic Massif (Northern Italy) have been used. They comprise four olivine texture types, which are characteristic for olivine textures observed worldwide. The bulk rock elastic properties have been calculated using olivine and orthopyroxene textures, their single-crystal elastic constants at ambient pressure/temperature conditions and their volume fraction. Clinopyroxene and spinel are assumed to be randomly distributed. The effect of four different orientations of the foliation within the uppermost mantle has been evaluated, since this orientation is usually unknown.¶Two of the olivine textures have a pronounced azimuthal dependence of compressional waves when a horizontal foliation within the uppermost mantle is presumed. These variations cause significant azimuthal variations of the P-wave reflections coefficients at the Moho. Primarily, we predict a significant azimuthal dependence of the critical points where the reflected amplitude increases from approximately 15% to 95%. Possibly, these azimuthal variations can be detected by seismic reflection measurements carried out at earth surface.¶The remaining two texture types only manifest a small directional dependence. When anisotropy of compressional waves is observed in seismic studies, these latter types can only be of subordinate importance. However, all of the peridotites investigated are able to explain the seismically observed azimuthal variations of compressional waves when a vertical foliation is proposed. This ambiguity can be substantially reduced when shear waves (S-waves) are considered. The directional distribution of S-wave velocities and of the S-wave splitting exhibits characteristic patterns for the different olivine texture types. This could be used to discriminate between different texture types and orientations of the foliation within the uppermost mantle. A fundamental requirement for a more comprehensive interpretation is the availability of detailed S-wave observations. The maximum S-wave splitting in the peridotites investigated coincides with the maximum of the faster (leading) S-wave. This may be of importance to detect S-wave splitting in future seismic studies.  相似文献   

17.
Ray equation migration of wide-angle reflections in 2-D medium is one kind of Kirchhoff prestack depth migration method. Based on ray theory, this method can be used extensively in 2-D inhomogeneous medium, and shows its advantage in wide-angle reflection study. After calculating ray fields, we can get the wave fields of sources and receivers by interpolation, and the intensity (or amplitude) A p corresponding to the depth point P at travel-time t P =t SP +t PR can be found in seismograms. Ray equation migration is completed by transforming A P according to the specific relation, migrating it to the depth point, then calculating the ray fields of lots of sources and receivers in the same way and finally stacking the outcomes. Numerical calculation has yielded satisfactory results. Data processing of the Zhuangmu-Yuexi-Huangmei-Zhanggondu profile passing through the Dabie orogenic zone provides a structural form of M-discontinuity along the profile. The result shows that the high pressure metamorphic rock zone of South Dabie rock mass is related to the uplifting of M-discontinuity here. "Mountain root" exists under North Dabie rock mass, which conforms with gravity isostasy theory. The Xiaotian-Mozitan fault is a suture zone as a result of the collision of North China plate with Yangtze plate, and extends through M-discontinuity. The abyssal fault near Liu’an is the extended eastern section of the Luonan-Minggang fault, which is also confirmed here. Contribution No. RCEG 96012, Research Center of Exploration Geophysics, SSB, China. This project was supported by the National Natural Science Foundation of China and the Chinese Joint Seismological Science Foundation.  相似文献   

18.
Microearthquake spectra from the Shillong region are analyzed to observe the effect of attenuation and site on these spectra. The spectral ratio method is utilized to estimate the Q values for both P- and S-waves in the subsurface layer, wherein the ratio of spectral amplitudes at lower and higher frequencies are taken into consideration for three stations at varying epicentral distances. Average estimates of Q P and Q S are 178 and 195. The ratio of Q S to Q P is estimated to be greater than 1 in major parts of the Shillong area, which can be related to the dry crust prevailing in the Shillong region. Typically, the variation in corner frequencies for these spectra is inferred to be characteristic of the site. Simultaneously, observations from spectral content of local earthquakes recorded at two different stations with respect to the reference site yield greater amplification of incoming seismic signals in the frequency range of 2–5 Hz, which is found to be well supported by the existing local lithology pertinent to that region.  相似文献   

19.
Records of deep-focus Hindu Kush earthquakes in the depth ranges 70–110 and 190–230 km made by 45 digital and analogue seismic stations were analyzed to study the attenuation field of short period seismic waves in the lithosphere of central Tien Shan. The dynamic characteristics studied include the ratio of peak amplitudes in S and P waves (S/P) and the ratio of the S-wave maximum to the coda level in the range t = 400 ± 5 s, where t is the lapse time (S/c400) for 1.25 Hz. Comparatively high values of S/P are shown to prevail in most of the area, corresponding to lower S-wave attenuation. Upon this background is a band of high and intermediate attenuation in the west of the area extending along the Talas-Fergana fault in the south and afterwards turning north-northeast. The rupture areas of the two largest (M ≥ 7.0) earthquakes which have occurred in Tien Shan during the last 25 years are confined to this band. Abnormally high values of S/c400 were obtained for stations situated in the rupture zone of the August 19, 1992, magnitude 7.3 Suusamyr earthquake and around it. For two of the stations we found considerable time variations in the coda envelope before the earthquake. The effective Q was derived from compressional and shear wave data for the entire area, as well as for the band of high attenuation. Comparison with previous data shows that the attenuation field in the area has changed appreciably during 20–25 years, which can only be due to a rearrangement of the fluid field in the crust and uppermost mantle. It is hypothesized that a large earthquake is very likely to occur in the northern part of the attenuating band.  相似文献   

20.
— A P-wave tomographic method for 3-D complex media (3-D distribution of elastic parameters and curved interfaces) with orthorhombic symmetry is presented in this paper. The technique uses an iterative linear approach to the nonlinear travel-time inversion problem. The hypothesis of orthorhombic anisotropy and 3-D inhomogeneity increases the set of parameters describing the model dramatically compared to the isotropic case. Assuming a Factorized Anisotropic Inhomogeneous (FAI) medium and weak anisotropy, we solve the forward problem by a perturbation approach. We use a finite element approach in which the FAI medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, analytical expressions for rays and travel times, valid to first-order, are given for P waves in orthorhombic inhomogeneous media. More complex media can be modeled by introducing interfaces separating FAI media with different elastic properties. Simple formulae are given for the Fréchet derivatives of the travel time with respect to the elastic parameters and the interface parameters. In the weak anisotropy hypothesis the P-wave travel times are sensitive only to a subset of the orthorhombic parameters: the six P-wave elastic parameters and the three Euler angles defining the orientation of the mirror planes of symmetry. The P-wave travel times are inverted by minimizing in terms of least-squares the misfit between the observed and calculated travel times. The solution is approached using a Singular Value Decomposition (SVD). The stability of the inversion is ensured by making use of suitable a priori information and/or by applying regularization. The technique is applied to two synthetic data sets, simulating simple Vertical Seismic Profile (VSP) experiments. The examples demonstrate the necessity of good 3-D ray coverage when considering complex anisotropic symmetry.  相似文献   

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

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