首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio (i.e., ellipticity), and receiver function data to better resolve 1D crustal shear wave velocity (v S) structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensitivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast (e.g., due to the existence of crustal interfaces) and v P/v S ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the v S model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute v S model and then incorporate receiver function data in the joint inversion to obtain a finer v S model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal v S structures and with little initial model dependency.  相似文献   

2.
Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C v value and to explore its physical meaning experimentally. The result shows that the change of C v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C v value increases firstly and then decreases before the instability, and the instability occurs when C v value decreases to the level before increasing. This indicates that C v value may be a useful parameter for earthquake prediction.  相似文献   

3.
The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K sat and K sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V P/V S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K sat–μ, EEI, V P/V S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.  相似文献   

4.
Based on a comparison of the cases of a decrease in the ratio of A n/A μ (where A n and A μ are the amplitudes of the diurnal variations of the neutron and hard cosmic ray components) to the instants of the Earth crossing the neutral IMF, it has been indicated that the process of such crossing is most effective for stimulating large destructive earthquakes with a magnitude of M ≥ 6. The 11-year period in the cyclicity of the occurrence probability of the above earthquakes has been revealed.  相似文献   

5.
The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.  相似文献   

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

7.
The study of rock samples from the Upper Permian Khei-Yaga River section revealed an r-n-r-n-r magnetic polarity succession based on the prefolding characteristic component of natural remanent magnetization. With account for stratigraphic and previous magnetostratigraphic data on Lower Triassic rocks from the Khei-Yaga River section, the examined strata of the Pechora Group (Silova Formation) may be compared with magnetic zones R1P2u and N1P3t in the magnetostratigraphic scale of European Russia. The gap in the paleomagnetic record, which corresponds in the examined section to the interval of the mid-Severodvinian Stage to basal Induan Stage, Zone N1T included, is estimated to be 10 Ma long. It is assumed that this hiatus represents one of the local signs of the global Permian-Triassic crisis.  相似文献   

8.
Recent estimates of fracture energy G ′ in earthquakes show a power-law dependence with slip u which can be summarized as G ′ ∝ u a where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to G f : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G ′ and G f should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s2). The resulting fracture energy G f is similar to the seismological estimates, with G f and G ′ being comparable over most of the slip range. However, G f appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G ′ appears to increase further and surpasses G f at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.  相似文献   

9.
Results of investigation of the lithosphere in the Kamchatka seismic focal zone from dynamic characteristics of earthquake records obtained at regional stations are presented. It is assumed that the specificity of the source zone can be estimated by the relation Cr = K P ? bK S ? c characterizing relative energies (energy classes, according to [Fedotov, 1972]) of short period transverse and longitudinal waves in the source. Azimuthal, spatial, and temporal variations in Cr and their relation to focal mechanisms are examined. Spatiotemporal variations in this parameter are shown to be caused by the influence of variations in the conditions in the source zone (its substance or process) on the radiation of P and S waves.  相似文献   

10.
Lower Cretaceous C-isotope records show intermittent negative/positive spikes, and consistent patterns of coeval chemostratigraphic curves thus document shifts that signal simultaneous responses of temporal changes in the global carbon reservoir. The standard pattern registered by the δ 13Corg and δ 13Ccarb in Lower Aptian sediments includes distinct isotope segments C1 to C8 (Menegatti et al., 1998). In the El Pui section, Organyà Basin, Spain, C-isotope segment C2 is the longest interval preceding segments C3–C6 associated with oceanic anoxic event 1a (OAE 1a), and reveals a distinct negative shift of ~1.8‰ to ~2.23‰ defining the C-isotope pattern within that interval. Total inorganic carbon (TIC), total organic carbon (TOC), δ 13Corg, microfacies, n-alkanes show no difference before, during, or after the negative inflection. The biomarkers indicate that organic matter (OM) mainly originates from algal/microbial sources because short-chain length homologues (≤nC19) dominate. nC20 through nC25 indicate some contribution from aquatic vegetation, but little from higher plants (>nC25), as also suggested by the terrestrial/aquatic ratio of n-alkanes or (TAR) = [(nC27+nC29+nC31)/(nC15+nC17+nC19)] (averages 0.085). We suggest that conjoint pulses of contemporaneous LIPs (Ontong Java) and massive explosive volcanism in northeast Asia, the Songliao Basin (SB-V), best conform to plausible causes of the negative intra-C2 carbon isotopic excursion (CIE) at that time. Because of its apparent common occurrence the intra-C2 inflection could be a useful marker harbinger to the more pronounced CIE C3, the hallmark of OAE1a.  相似文献   

11.
We made an attempt to assess the shear wave velocity values V S and, to a lesser extent, the V P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.  相似文献   

12.
An alternative model for the nonlinear interaction term Snl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of Snl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of Snl, and to compare it to the existing models for Snl in WAVEWATCH; the DIA and WRT. Some differences between the different models for Snl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of Snl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.  相似文献   

13.
We describe an algorithm for rapidly computing the surface displacements induced by a general polygonal load on a layered, isotropic, elastic half-space. The arbitrary surface pressure field is discretized using a large number, n, of equally-sized circular loading elements. The problem is to compute the displacement at a large number, m, of points (or stations) distributed over the surface. The essence of our technique is to reorganize all but a computationally insignificant part of this calculation into an equivalent problem: compute the displacements due to a single circular loading element at a total of m n stations (where m n is the product m × n). We solve this “parallel” problem at high computational speed by utilizing the sparse evaluation and massive interpolation (SEMI) method. Because the product m n that arises in our parallel problem is normally very large, we take maximum possible advantage of the acceleration achieved by the SEMI algorithm.  相似文献   

14.
Aftershock hazard maps contain the essential information for search and rescue process, and re-occupation after a main-shock. Accordingly, the main purposes of this article are to study the aftershock decay parameters and to estimate the expected high-frequency ground motions (i.e., Peak Ground Acceleration (PGA)) for recent large earthquakes in the Iranian plateau. For this aim, the Ahar-Varzaghan doublet earthquake (August 11, 2012; M N =6.5, M N =6.3), and the Ilam (Murmuri) earthquake (August 18, 2014 ; M N =6.2) have been selected. The earthquake catalogue has been collected based on the Gardner and Knopoff (Bull Seismol Soc Am 64(5), 1363-1367, 1974) temporal and spatial windowing technique. The magnitude of completeness and the seismicity parameters (a,??b) and the modified Omori law parameters (P,??K,??C) have been determined for these two earthquakes in the 14, 30, and 60 days after the mainshocks. Also, the temporal changes of parameters (a,??b,??P,??K,??C) have been studied. The aftershock hazard maps for the probability of exceedance (33%) have been computed in the time periods of 14, 30, and 60 days after the Ahar-Varzaghan and Ilam (Murmuri) earthquakes. For calculating the expected PGA of aftershocks, the regional and global ground motion prediction equations have been utilized. Amplification factor based on the site classes has also been implied in the calculation of PGA. These aftershock hazard maps show an agreement between the PGAs of large aftershocks and the forecasted PGAs. Also, the significant role of b parameter in the Ilam (Murmuri) probabilistic aftershock hazard maps has been investigated.  相似文献   

15.
We study the frictional and viscous effects on earthquake nucleation, especially for the nucleation phase, based on a one-degree-of-freedom spring-slider model with friction and viscosity. The frictional and viscous effects are specified by the characteristic displacement, U c, and viscosity coefficient, η, respectively. Simulation results show that friction and viscosity can both lengthen the natural period of the system and viscosity increases the duration time of motion of the slider. Higher viscosity causes a smaller amplitude of lower velocity motion than lower viscosity. A change of either U c (under large η) or η (under large U c) from a large value (U ch for U c and η h for η) to a small one (U cl for U c and η l for η) in two stages during sliding can result in a clear nucleation phase prior to the P-wave. The differences δU c = U ch ? U cl and δη = η h ? η l are two important factors in producing a nucleation phase. The difference between the nucleation phase and the P-wave increases with either δU c or δη. Like seismic observations, the peak amplitude of P-wave, which is associated with the earthquake magnitude, is independent upon the duration time of nucleation phase. A mechanism specified with a change of either η or U c from a larger value to a smaller one due to temporal variations in pore fluid pressure and temperature in the fault zone based on radiation efficiency is proposed to explain the simulation results and observations.  相似文献   

16.
The composition characteristics of light hydrocarbons from crude oil, chloroform bitumen A, saturated hydrocarbon fraction, aromatic hydrocarbon fraction, and asphaltene fraction during cracking have been studied systematically. The results revealed that the content of n-alkanes, branched alkanes and cycloalkanes in light hydrocarbons from the samples gradually decreased as the simulation temperature increased, and finally almost depleted completely, while the abundance of methane, benzene and its homologues increased obviously and became the main products. The ratios of benzene/ n-hexane and toluene/n-heptane can be used as measures for oil cracking levels. Variation characteristics of maturity parameters of light hydrocarbons, for example, iC4/nC4, iC5/nC5, isoheptane value, 2,2-DMC4/nC6, and 2-MC6+3-MC6/nC7 for different samples with increasing pyrolysis temperature, are consistent with those in petroleum reservoirs, indicating that these parameters may be efficient maturity index.  相似文献   

17.
This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f 0(P)/f 0(S). The ratios were computed for each event with corner frequencies f 0 of P- and S-wave, which were obtained from the measured f 0 I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f 0(P)/f 0(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f 0(P)/f 0(S), with the empirical threshold value of the ratio for Israel as 1.48.  相似文献   

18.
Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν j is computed at each layer. For a given phase velocity (c), ν j becomes complex depending on the transversely isotropic parameters. When ν j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α H , β V , ξ, ? and η and find their limits beyond which ν j becomes complex. It is seen that ν j depends on ? and η, but does not depend on ξ. The complex ν j occurs when ? is small and η is large. For a given c/β V , the region of complex ν j in a ? -η plane increases with the increase of α H /β V . Further, for a given α H /β V , the complex region of ν j increases significantly with the decrease of c/β V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.  相似文献   

19.
Statistical tests have been used to adjust the Zemmouri seismic data using a distribution function. The Pareto law has been used and the probabilities of various expected earthquakes were computed. A mathematical expression giving the quantiles was established. The extreme values limiting law confirmed the accuracy of the adjustment method. Using the moment magnitude scale, a probabilistic model was made to predict the occurrences of strong earthquakes. The seismic structure has been characterized by the slope of the recurrence plot γ, fractal dimension D, concentration parameter Ksr, Hurst exponents Hr and Ht. The values of D, γ, Ksr, Hr, and Ht diminished many months before the principal seismic shock (M = 6.9) of the studied seismoactive zone has occurred. Three stages of the deformation of the geophysical medium are manifested in the variation of the coefficient G% of the clustering of minor seismic events.  相似文献   

20.
We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V P values measured for the soils. The downtown part of Baku city is characterized by low VS30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.  相似文献   

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

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