首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 641 毫秒
1.
  相似文献   

2.
The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log10E2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log10E2/3, and log10N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log10E2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log10Ni and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log10Ni time series in the different magnitude ranges (M ≥ Mci) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔMi) in the range of low magnitudes are known. It is concluded that using the parameter log10N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.  相似文献   

3.
Aftershock sequences of some strong earthquakes of Kamchatka, the Kurile Islands, and Japan are examined. Such source parameters as the length L, along-dip width W, motion on fault D, and stress drop Δσ are determined from the aftershock sequences considered. The values of these parameters were obtained by the formal estimation of linear source parameters (lower bound estimates) and visually (upper bound estimates). The correlation dependences of the obtained parameters on the surface wave (M S ) and seismic moment (M W ) magnitudes are calculated.  相似文献   

4.
During the ruptures of an earthquake,the strain energy.△E,.will be transferred into,at least,three parts,i.e..the seismic radiation energy(E_s),fracture energy(E_g),and frictional energy(E_f),that is,△E = E_s + E_g + E_f.Friction,which is represented by a velocity- and state-dependent friction law by some researchers,controls the three parts.One of the main parameters of the law is the characteristic slip displacement.D_c.It is significant and necessary to evaluate the reliable value of D_c from observed and inverted seismic data.Since D_c controls the radiation efficiency.η_R = E_s/(E_s+ E_g),the value of η_r is a good constraint of estimating D_c.Integrating observed data and inverted results of source parameters from recorded seismograms.the values of E_s and E_g of an earthquake can be measured,thus leading to the value of η_R.The constraint used to estimate the reliable value of D_c will be described in this work.An example of estimates of D_c.based on the observed and inverted values of source parameters of the September 20,1999 M_S 7.6 Chi-Chi(Ji-Ji).Taiwan region,earthquake will be presented.  相似文献   

5.
In this study, we aim to improve the scaling between the moment magnitude (M W), local magnitude (M L), and the duration magnitude (M D) for 162 earthquakes in Shillong-Mikir plateau and its adjoining region of northeast India by extending the M W estimates to lower magnitude earthquakes using spectral analysis of P-waves from vertical component seismograms. The M W-M L and M W-M D relationships are determined by linear regression analysis. It is found that, M W values can be considered consistent with M L and M D, within 0.1 and 0.2 magnitude units respectively, in 90 % of the cases. The scaling relationships investigated comply well with similar relationships in other regions in the world and in other seismogenic areas in the northeast India region.  相似文献   

6.
Based on the example of the Vrancea zone of concentrated seismicity, it is shown how the stress-strain state of the medium responds to a disturbance of the geomagnetic field. Geomagnetic conditions are examined in relation to earthquakes in the Vrancea zone in the period 1988–1996. It is established that the seismic energy release in the Vrancea zone is associated with differences (“gradients”) in the H component of the geomagnetic field. Such a gradient preceding earthquakes is shown to be the midnight polar substorm and the degree of its mid-latitude effect. The time interval from the maximum of the substorm development to a shock (τ, h) is directly related to the focal depth. The seismic characteristics K en and h (km) are demonstrated to be related to morphological features of the substorm development, namely, its duration T (min), intensity, and background. Differences in the duration of polar substorms before crustal (shallow) and deep earthquakes are revealed. Morphological features of the spectrum of geomagnetic variations preceding the seismic energy release are established.  相似文献   

7.
We continue applying the general concept of seismic risk analysis in a number of seismic regions worldwide by constructing regional seismic hazard maps based on morphostructural analysis, pattern recognition, and the Unified Scaling Law for Earthquakes (USLE), which generalizes the Gutenberg-Richter relationship making use of naturally fractal distribution of earthquake sources of different size in a seismic region. The USLE stands for an empirical relationship log10N(M, L)?=?A?+?B·(5 – M)?+?C·log10L, where N(M, L) is the expected annual number of earthquakes of a certain magnitude M within a seismically prone area of linear dimension L. We use parameters A, B, and C of USLE to estimate, first, the expected maximum magnitude in a time interval at seismically prone nodes of the morphostructural scheme of the region under study, then map the corresponding expected ground shaking parameters (e.g., peak ground acceleration, PGA, or macro-seismic intensity). After a rigorous verification against the available seismic evidences in the past (usually, the observed instrumental PGA or the historically reported macro-seismic intensity), such a seismic hazard map is used to generate maps of specific earthquake risks for population, cities, and infrastructures (e.g., those based on census of population, buildings inventory). The methodology of seismic hazard and risk assessment is illustrated by application to the territory of Greater Caucasus and Crimea.  相似文献   

8.
Based on the Anapa (ANN) seismic station records of ~40 earthquakes (MW > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.  相似文献   

9.
Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M L ?≤?4.4. Up to five seismograph stations were utilised to determine Q c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q 0 and α were determined for the attenuation relation Q c(f)?=?Q 0 f α . The result was Q 0?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t s???t 0) (time from origin time t 0 to the start of coda analysis window is 1.9 times the S-travel time, t s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q c was not calculated. The derived Q c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q S was Q 0?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q c result.  相似文献   

10.
We analyzed the most relevant seismic sequences that occurred from 1977 to 2007 in the Friuli-Venezia Giulia region (northeastern Italy) and western Slovenia. The eight aftershock sequences were triggered by low- to moderate-magnitude earthquakes with mainshock duration magnitude ranging from 3.7 to 5.6. The b-value of the Gutenberg–Richter law varies from 0.8 to 1.1. The modified Omori’s modeling of the sequences evidences values of the p exponent ranging from 0.8 to 1.0. Using the Reasenberg and Jones (Science 243:1173–1176, 1989; Science 265:1251–1252, 1994) approach, we computed the probabilistic estimate of the aftershock rates and the largest aftershock in given time intervals. The difference in magnitude between the mainshock and the largest aftershock is calculated according to the modified Båth law and using an approach that considers the partitioning of the radiated seismic energy between mainshock and aftershocks. The partitioning of the radiated seismic energy appears to play a significant role in the evolution of the sequences. We define the parameter R ES as the ratio between the radiated seismic energy of the mainshock and the summation of the seismic energy radiated by the aftershocks. The difference in magnitude between the mainshock and the largest aftershock, calculated with the parameter R ES, agrees well with the observed difference. In most sequences, the parameter R ES decreases very quickly until the occurrence of the largest aftershock and then becomes constant. By analyzing the values of R ES during the early hours following the mainshock, we found that the R ES values after 24 h are well related to the final ones, calculated on the whole sequence, and to the differences in magnitude between the mainshock and the largest aftershock.  相似文献   

11.
In the present study, the level of the largest earthquake hazard is assessed in 28 seismic zones of the NW Himalaya and its vicinity, which is a highly seismically active region of the world. Gumbel’s third asymptotic distribution (hereafter as GIII) is adopted for the evaluation of the largest earthquake magnitudes in these seismic zones. Instead of taking in account any type of Mmax, in the present study we consider the ω value which is the largest earthquake magnitude that a region can experience according to the GIII statistics. A function of the form Θ(ω, RP6.0) is providing in this way a relatively largest earthquake hazard scale defined by the letter K (K index). The return periods for the ω values (earthquake magnitudes) 6 or larger (RP6.0) are also calculated. According to this index, the investigated seismic zones are classified into five groups and it is shown that seismic zones 3 (Quetta of Pakistan), 11 (Hindukush), 15 (northern Pamirs), and 23 (Kangra, Himachal Pradesh of India) correspond to a “very high” K index which is 6.  相似文献   

12.
The purpose of this work is to define a seismic regionalization of Mexico for seismic hazard and risk analyses. This seismic regionalization is based on seismic, geologic, and tectonic characteristics. To this end, a seismic catalog was compiled using the more reliable sources available. The catalog was made homogeneous in magnitude in order to avoid the differences in the way this parameter is reported by various agencies. Instead of using a linear regression to converts from m b and M d to M s or M w , using only events for which estimates of both magnitudes are available (i.e., paired data), we used the frequency-magnitude relations relying on the a and b values of the Gutenberg-Richter relation. The seismic regions are divided into three main categories: seismicity associated with the subduction process along the Pacific coast of Mexico, in-slab events within the down-going COC and RIV plates, and crustal seismicity associated to various geologic and tectonic regions. In total, 18 seismic regions were identified and delimited. For each, the a and b values of the Gutenberg-Richter relation were determined using a maximum likelihood estimation. The a and b parameters were repeatedly estimated as a function of time for each region, in order to confirm their reliability and stability. The recurrence times predicted by the resulting Gutenberg-Richter relations obtained are compared with the observed recurrence times of the larger events in each region of both historical and instrumental earthquakes.  相似文献   

13.
Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.  相似文献   

14.
The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M0?=?1.12×1017 Nm, M w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω?2 source model) by attenuation parameters Q(f)?=?500f0.9, κ?=?0.04 s, and fmax?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.  相似文献   

15.
Based on the analysis of the world’s earthquakes with magnitudes M ≥ 6.5 for 1960–2013, it is shown that they cause global-scale coherent seismic oscillations which most distinctly manifest themselves in the period interval of 4–6 min during 1–3 days after the event. After these earthquakes, a repeated shock has an increased probability to occur in different seismically active regions located as far away as a few thousand km from the previous event, i.e., a remote interaction of seismic events takes place. The number of the repeated shocks N(t) decreases with time, which characterizes the memory of the lithosphere about the impact that has occurred. The time decay N(t) can be approximated by the linear, exponential, and powerlaw dependences. No distinct correlation between the spatial locations of the initial and repeated earthquakes is revealed. The probable triggering mechanisms of the remote interaction between the earthquakes are discussed. Surface seismic waves traveling several times around the Earth’s, coherent oscillations, and global source are the most preferable candidates. This may lead to the accumulation and coalescence of ruptures in the highly stressed or weakened domains of a seismically active region, which increases the probability of a repeated earthquake.  相似文献   

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

17.
Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T p , respectively, versus M are inferred from compiled data. The difference between the values of T and T p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.  相似文献   

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

19.
The effect of baric variations of different origins on characteristics of seismic noise is analyzed in the frequency range 0.03–20 Hz. Long period variations in atmospheric pressure caused by cyclones, whose period T ranges from half a day to a few days, are shown to increase the microseismic background amplitude by two to four times in the frequency range 0.03–1 Hz (the coefficient of linear correlation between time variations in the amplitude and atmospheric pressure is K = 0.65 at a significance level of r = 0.95). Short-period baric variations with T ～ 5–30 min associated with the passage of cold fronts lead to a tenfold increase in the microseismic background amplitude in the frequency range 4–8 Hz (K = 0.67 at r = 0.95). In this case, disturbances of seismic background are recorded for 20–60 min after the passage of an atmospheric front and display an exponential drop in the amplitude. In distinction to cyclones, an atmospheric front increases the number of impulsive microseismic events of the resonance type.  相似文献   

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

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

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