首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (QP and QS) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 MW4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average vP and vS of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high vP/vS ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that QS is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with QS modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.  相似文献   

2.
The single scattering model has been applied for the estimation of codaQ values for local earthquakes that occurred in northern Greece during the period 1983–1989 and recorded by the telemetered network of the Geophysical Laboratory of the University of Thessaloniki. CodaQ estimations were made for four frequency bands centered at 1.5 Hz, 3.0 Hz, 6.0 Hz and 12.0 Hz and for the lapse time windows 10–20 sec, 15–30 sec, 20–45 sec, 30–60 sec and 50–100 sec. The codaQ values obtained show a clear frequency dependence of the formQ c =Q 0 f n , whileQ 0 andn depend on the lapse time window.Q 0 was found equal to 33 andn equal to 1.01 for the time window of 10 to 20 sec, while for the other windowsQ 0 increased from 60 to 129, withn being stable, close to 0.75. This lapse time dependence is interpreted as due to a depth dependent attenuation. The high attenuation and the strong frequency dependence found are characteristic of an area with high seismicity, in agreement with studies in other seismic regions.  相似文献   

3.
We analyzed the local earthquakes waveform recorded on a broadband seismic network in the northwestern Himalayan Region to compute lapse time and frequency dependence of coda Q (Q c). The observed Q c values increase with increasing lapse time at all frequency bands. The increase in Q c values with lapse time is attributed to an increase in Q c with depth. This implies that attenuation decreases with increasing depth. The approximate radius of medium contributing to coda generation varies from 55 to 130 km. By comparing the Q c values with those from other regions of the world, we find that they are similar to those obtained from tectonically active regions. The estimated Q c values show a frequency-dependent relationship, Q c = Q 0 f n , where Q 0 is Q c at 1 Hz and n represents degree of frequency dependence. They represent the level of heterogeneity and tectonic activity in an area. Our results show that northwest Himalayas are highly heterogeneous and tectonically very active. Q 0 increases from 113 ± 7 to 243 ± 10 and n decreases from 1.01 ± 0.05 to 0.85 ± 0.03 when lapse time increases from 30 to 70 s. As larger time window sees the effect of deeper part of the Earth, it is concluded that Q 0 increases and n decreases with increasing depth; i.e., heterogeneity decreases with depth in the study area.  相似文献   

4.
We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q P and Q S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.  相似文献   

5.
The seismic energy attenuation in the frequency range of 1–18 Hz was studied in the two tectonically active zones of Irno Valley (Southern Italy) and Granada Basin (South-East Spain). Data were recorded by short period vertical components seismographs for low-magnitude local earthquakes. The method of coda waves, assuming singleS toS scattering approximation, was used to calculate the quality factorQ from the two data set. Results show a quality factor increasing with frequency, following the empirical lawQ=Q o f n .Q o andn are lower for the Irno Valley than for Granada. This result is interpreted in terms of different scattering environments present in the two investigated areas.  相似文献   

6.
The use of relaxation mechanisms has recently made it possible to simulate viscoelastic (Q) effects accurately in time-domain numerical computations of seismic responses. As a result, seismograms may now be synthesized for models with arbitrary spatial variations in compressional- and shear-wave quality factors (Q9, and Qs, as well as in density (ρ) and compressional- and shear-wave velocities (Vp, and Vs). Reflections produced by Q contrasts alone may have amplitudes as large as those produced by velocity contrasts. Q effects, including their interaction with Vp, Vs and p, contribute significantly to the seismic response of reservoirs. For band-limited data at typical seismic frequencies, the effects of Q on reflectivity and attenuation are more visible than those on dispersion. Synthetic examples include practical applications to reservoir exploration, evaluation and monitoring. Q effects are clearly visible in both surface and offset vertical seismic profile data. Thus, AVO analyses that neglect Q may produce erroneous conclusions.  相似文献   

7.
An approach to extraction of viscoelastic parameters from seismic data is implemented and succesfully tested. Viscoelastic inversion is performed using adaptive damping factors to control the sensitivity of the viscoelastic parameters in relation to the τ-p seismic data. A priori information is incorporated through the damping factors as standard deviations of the data and of the viscoelastic model parameters. The stability of the inversion process is controlled by the variation of the damping factors as a function of the residual errors and parameter updates at each iteration. Tests on synthetic and real data show that P- and S-wave quality factors, Qp and Qs, in addition to P- and S-wave velocities and density Cp, Cs and p, can be extracted from τ-p seismic information. Singular value decomposition analysis demonstrates that estimated Qp and Qs values are more affected by the presence of data inaccuracies and noise than are those of Cp and p. Cs and Qs are not uniquely recovered due to the limited contribution of P-S converted waves. Knowledge of the viscoelastic parameters is of particular importance in accurately describing petrophysical properties of rocks and pore fluids existing in the subsurface; this is demonstrated with real data from the Gulf of Mexico.  相似文献   

8.
We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective QS/QP ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.  相似文献   

9.
Body-wave Attenuation in the Region of Garda, Italy   总被引：1，自引：0，他引：1  
We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (ML 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q P = 65 f 0.9 and Q S = 160 f 0.6 , respectively. At 2 Hz the Q S /Q P ratio reaches the highest value of 2.8. At higher frequencies Q S /Q P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q S /Q P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.  相似文献   

10.
The time‐invariant gain‐limit‐constrained inverse Q‐filter can control the numerical instability of the inverse Q‐filter, but it often suppresses the high frequencies at later times and reduces the seismic resolution. To improve the seismic resolution and obtain high‐quality seismic data, we propose a self‐adaptive approach to optimize the Q value for the inverse Q‐filter amplitude compensation. The optimized Q value is self‐adaptive to the cutoff frequency of the effective frequency band for the seismic data, the gain limit of the inverse Q‐filter amplitude compensation, the inverse Q‐filter amplitude compensation function, and the medium quality factor. In the processing of the inverse Q‐filter amplitude compensation, the optimized Q value, corresponding gain limit, and amplitude compensation function are used simultaneously; then, the energy in the effective frequency band for the seismic data can be recovered, and the seismic resolution can be enhanced at all times. Furthermore, the small gain limit or time‐variant bandpass filter after the inverse Q‐filter amplitude compensation is considered to control the signal‐to‐noise ratio, and the time‐variant bandpass filter is based on the cutoff frequency of the effective frequency band for the seismic data. Synthetic and real data examples demonstrate that the self‐adaptive approach for Q value optimization is efficient, and the inverse Q‐filter amplitude compensation with the optimized Q value produces high‐resolution and low‐noise seismic data.  相似文献   

11.
Summary P n velocities determined from seismic refraction measurements, show significant differences between Southern Finland (7.96 km/sec) and Southeastern Norway (8.20 km/sec). TheP n/Sn velocity ratios (k) were determined from earthquake and explosion data, and the observed variation ofk indicates lateral variations in theP n and/orS n velocities in Fennoscandia.  相似文献   

12.
On 26 January, 2001 (03:46:55,UT) a devastating intraplate earthquake of Mw 7.7 occurred in a region about 5 km NW of Bhachau, Gujarat (23.42°N, 70.23°E). The epicentral distribution of aftershocks defines a marked concentration along an E-W trending and southerly dipping (45°) zone covering an area of (60 × 40) km2. The presence of high seismicity including two earthquakes of magnitudes exceeding 7.7 in the 200 years is presumed to have caused a higher level of shallow crustal heterogeneity in the Kutch area; a site lying in the seismic zone V (zone of the highest seismicity for potentially M8 earthquakes) on the seismic zoning map of India. Attenuation property of the medium around the epicentral area of the Bhuj earthquake covering a circular area of 61,500 km2 with a radius of 140 km is studied by estimating the coda-Qc from 200 local earthquakes of magnitudes varying from 3.0–4.6. The estimated Q0 values at locations in the aftershock zone (high seismicity) are found to be low in comparison to areas at a distance from it. This can be attributed to the fact that seismic waves are highly scattered for paths through the seismically active and fractured zone but they are well behaved outside the aftershock zone. Distribution of Q0 values suggests that the local variation in Q0 values is probably controlled by local geology. The estimated Q0 values at different stations suggest a low value of Q=(102 ± 0.80)*f(0.98 ± 0.02) indicating an attenuative crust beneath the entire region. The frequency-dependent relation indicates a relatively low Qc at lower frequencies (1–3 Hz) that can be attributed to the loss of energy due to scattering attenuation associated with heterogeneities and/or intrinsic attenuation due to fluid movement in the fault zone and fluid-filled cracks. The large Qc at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneity is expected. Based on the attenuation curve estimated for Q0=102, the ground acceleration at 240 km distance is 13% of 1 g i.e., 0.13 g agreeing well with the ground acceleration recorded by an accelerograph at Ahmedabad (0.11 g). Hence, it is inferred that the Q0 value obtained from this study seems to be apt for prediction of ground motion for the region.  相似文献   

13.
The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The amplitude spectral ratio method is used, assuming no frequency dependence of the Q P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor is observed. Lateral variations of the attenuation are obtained, with a very low Q P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different azimuth angles in relation to the epicentral area, Q P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located immediately beneath the crust, but limited in depth to at least 100 km.  相似文献   

14.
Abstract

The normal mode instability of harmonic waves in an ideal incompressible fluid on a rotating sphere is analytically studied. By the harmonic wave is meant a Legendrepolynomial flow αPn(μ) (n ≥ 1) and steady Rossby-Haurwitz wave of set F 1 ⊕ Hn where Hn is the subspace of homogeneous spherical polynomials of the degree n(n ≥ 2), and F 1 is the one-dimensional subspace generated by the Legendre-polynomial P1(μ). A necessary condition for the normal mode instability of the harmonic wave is obtained. By this condition, Fjörtoft's (1953) average spectral number of the amplitude of each unstable mode must be equal to . It is noted that flow αPn (μ) is Liapunov (and hence, exponentially and algebraically) stable to all the disturbances whose zonal wavenumber m satisfies condition |m| ≥ n. The bounds of the growth rate of unstable normal modes are estimated as well. It is also shown that the amplitude of each unstable, decaying or non-stationary mode is orthogonal to the harmonic wave.

The new instability condition can be useful in the search of unstable perturbations to a harmonic wave and on trials of numerical stability study algorithms. For a Legendre-polynomial flow, it complements Kuo's (1949) condition in the sense that while the latter is related to the basic flow structure; the former characterizes the structure of a growing perturbation.  相似文献   

15.
Viscoelastic modelling reveals that the interaction of compressional-wave velocity Cp, compressional-wave quality factor Qp, shear-wave velocity Cs, shear-wave quality factor Qs and Poisson's ratio as a function of time intercept τ and ray parameter p, is complicated; however, distinct, potentially diagnostic behaviours are seen for different combinations of viscoelastic parameters. Synthetic seismograms for three viscoelastic reservoir models show that variations in the Poisson's ratio produce visible differences when compared to the corresponding elastic synthetic seismograms; these differences are attributable to interaction of the elastic parameters with Qp and Qs. When the P-wave acoustic impedance contrast is small, viscoelastic effects become more apparent and more useful for interpretation purposes. The corresponding amplitude and net phase spectra reveal significant differences between the elastic and the viscoelastic responses. When P-wave reflectivities are large, they tend to dominate the total response and to mask the Q reflectivity effects. The attenuation effects are manifested as an amplitude decay that increases with both time and ray parameter. The sensitivity of the computed seismic responses for various combinations of viscoelastic parameters suggests the opportunity for diagnostic interpretation of τ-p seismic data. The interpretation of the viscoelastic parameters can permit a better understanding of the rock types and pore fluid distribution existing in the subsurface.  相似文献   

16.
Attenuation of seismic waves is very essential for the study of earthquake source parameters and also for ground-motion simulations, and this is important for the seismic hazard estimation of a region. The digital data acquired by 16 short-period seismic stations of the Delhi Telemetric Network for 55 earthquakes of magnitude 1.5 to 4.2, which occurred within an epicentral distance of 100 km in an area around Delhi, have been used to estimate the coda attenuation Qc. Using the Single Backscattering Model, the seismograms have been analyzed at 10 central frequencies. The frequency dependence average attenuation relationship Qc = 142f 1.04 has been attained. Four Lapse-Time windows from 20 to 50 seconds duration with a difference of 10 seconds have been analyzed to study the lapse time dependence of Qc. The Qc values show that frequency dependence (exponent n) remains similar at all the lapse time window lengths. While the change in Q0 values is significant, change in Q0 with larger lapsetime reflects the rate of homogeneity at the depth. The variation of Qc indicates a definitive trend from west to east in accordance with the geology of the region.  相似文献   

17.
Summary Earthquake parameters for the forty aftershocks of the main Koyna earthquake of 10 December, 1967, have been determined. Depths of the foci of the earthquakes have been found to vary between 2 to 17 km. The velocities for the phasesP g ,P *,P n have been observed to be 5.78±0.00, 6.58±0.04, 8.19±0.02 km/sec, and forS g ,S *,S n to be 3.42±0.00, 3.92±0.01 and 4.62±0.01 km/sec respectively. A two-layered crustal model has been interpreted for the Peninsular shield with the average thickness of the granitic layer as 20 km and that for the basaltic layer as 18.7 km. A plot of the epicenters suggests a NNE to SSW orientation of the fault.  相似文献   

18.
The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q P) and S (Q S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q P and Q S are: Q P = (105 ± 2) f 0.82 ± 0.01, Q S = (74 ± 2) f 1.06 ± 0.01 for Kachchh region; Q P = (148 ± 2) f 0.92 ± 0.01, Q S = (149 ± 14) f 1.43 ± 0.05 for Saurashtra region and Q P = (163 ± 7) f 0.77 ± 0.03, Q S = (118 ± 34) f 0.65 ± 0.14 for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q C > Q S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q S/Q P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.  相似文献   

19.
Wide‐angle seismic data acquired by use of air‐guns and ocean bottom seismometers (OBS) contain strong direct water arrivals and multiples, generally considered as noise and thus not included in the modelling. However, a recent study showed that standard ray‐tracing modelling of the water multiples recorded off the Bear Island, North Atlantic, provided a reliable estimate of the velocity distribution in the water layer. Here, we demonstrate that including the amplitudes in the modelling provide valuable information about the VP contrast at the seafloor, as well as the VP/VS ratio and attenuation (QP) of the uppermost sediments. The VP contrast at the seafloor is estimated at about 250 m/s, within a precision of approximately ±30 m/s. The VP/VS ratio in the uppermost sedimentary layer is modelled in the range 2.25–2.50 and the QP factor is estimated at 1000 for the water, 30–50 for the uppermost layer and 40–50 for the second sedimentary layer. The values obtained for the sediments suggest a lithology dominated by silty clays, with porosity below average.  相似文献   

20.
A synthetic seismogram that closely resembles a seismic trace recorded at a well may not be at all reliable for, say, stratigraphic interpretation around the well. The most accurate synthetic seismogram is, in general, not the one that displays the smallest errors of fit to the trace but the one that best estimates the noise on the trace. If the match is confined to a short interval of interest or if the seismic reflection wavelet is allowed to be unduly long, there is considerable danger of forcing a spurious fit that treats the noise on the trace as part of the seismic reflection signal instead of making a genuine match with the signal itself. This paper outlines tests that allow an objective and quantitative evaluation of the accuracy of any match and illustrates their application with practical examples. The accuracy of estimation is summarized by the normalized mean square error (NMSE) in the estimated reflection signal, which is shown to be (/n)(PN/PS) where PS/PN is the signal-to-noise power ratio and n is the spectral smoothing factor. That is, the accuracy varies directly with the ratio of the power in the signal (taken to be the synthetic) to that in the noise on the seismic trace, and the smoothing acts to improve the accuracy of the predicted signal. The construction of confidence intervals for the NMSE is discussed. Guidelines for the choice of the spectral smoothing factor n are given. The variation of wavelet shape due to different realizations of the noise component is illustrated, and the use of confidence intervals on wavelet phase is recommended. Tests are described for examining the normality and stationarity of the errors of fit and their independence of the estimated reflection signal.  相似文献   

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

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