首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Source parameter scaling relations are examined for microseismic events (–2.4M–0.3) occurring within higyly and moderately stressed and fractured rock masses at Strathcona mine, Sudbury, Canada. Insight into scaling is provided by waveform complexities, calculated rupture velocities, and maximum shear stresses based onin situ and numerical modelling data. The importance of normal stress on the failure process is also considered. Our results show that a strong dependence exists between stress release and seismic moment. An observed positive scaling in excess stress release (/2– a ) is consistent with the concept of overshoot. Rupture velocities ranging from 0.2 to 0.5 and waveform complexities less than 1.5 suggested that overshoot was related to healing behind a slowly advancing rupture front. Scaling in seismic efficiency paralleled that in apparent stress, implying that seismic stress release estimates are quasi-independent of the maximum shear stress. High levels of normal stress further supported the importance of high resisting stress in the observed overshoot behaviour and its role in the failure process.  相似文献   

2.
We estimated the source parameters of 53 local earthquakes (2.0<ML<5.7) of the Friuli-Venezia Giulia (Northeastern Italy) area, recorded by the short-period local seismic network of the Istituto Nazionale di Oceanografia e Geofisica Sperimentale (OGS), in the period 1995-2003. Data were selected on the basis of high quality locations and focal mechanisms. Standard H/V spectral ratios (HVRS) of the three-component stations of the network were performed in order to assess local amplifications, and only stations showing HVRS not exceeding two were considered for the source parameters estimation. Both velocity and acceleration data were used to compute the SH-wave spectra. Observed spectra were corrected for attenuation effects using an independent regional estimate of the quality factor Q and a station dependent estimate of the spectral decay parameter k. Only earthquakes with ML>3.0 recorded with a sampling rate of 125 cps were used to compute k, thus allowing to visualize a linear trend of the high frequency acceleration spectrum up to 40-50 Hz. SH-wave spectra, corrected for attenuation, showed an ω−2 shape allowing a good fit with the Brune model. Seismic moments and Brune radii ranged between 1.5×1012 and 1.1×1017  N m and between 0.1 and 2.7 km respectively. We obtained Mo=1.1×1017  N m for the seismic moment of the Kobarid (SLO) main shock, in good agreement with the Harvard CMT solution (Mo=3.5×1017  N m). Brune stress drops were confined to the range from 0.07 to 5.31 MPa, with an average value of 0.73 MPa and seem to be approximately constant over five orders of magnitude of seismic moment. Radiated seismic energy computed from two nearby stations scales with seismic moment according to , and apparent stress values are between 0.02 and 4.26 MPa. The observed scatter of Brune stress drop data allowed to hypothesize a scaling relation between seismic moment and corner frequency in order to accommodate both Brune stress drop and apparent stress scalings. No systematic differences are evidenced between stress parameters of earthquakes with different focal mechanisms. As a consequence, a relation of the seismic stress release with the strength of rocks can be hypothesized. A high correlation (r>0.9) of Brune stress drop is found with both apparent stress and RMS stress drop, according to and respectively.  相似文献   

3.
TeleseismicP-waves of some large earthquakes that occurred in the eastern Mediterranean region have been analysed by using an iterative maximum entropy technique in order to obtain the independent spectral parameters, the long-period spectral level 0 and the corner frequencyf 0 of the far-field displacement spectra.Based on these parameters, the seismic source parameters seismic momentM 0, source dimensionr, fault lengthl, average displacement u, shear stress drop , radiated energyE s and apparent stressn are calculated for the considered earthquakes by using Brune's and Madariaga's models.The striking feature of the source parameters obtained in this study is the low stress drop value which varies between 5 and 15 bars. If Madariaga's model had been used, higher stress drop values would have been obtained.The low stress drop earthquakes in the eastern Mediterranean region might be interpreted either by the possible presence of low strength material near the source or by the partial stress drop model.  相似文献   

4.
Intrinsic and scatteringS-wave quality factors (Q ) were estimated using the Multiple Lapse Time Window Analysis (MLTWA) for microseismic events (M<–1) with source-sensor distances of 45 to 120 m, associated with an excavation at 630 m depth in Strathcona Mine, Sudbury, Canada. Additional information on the rock mass was provided by underground structural mapping data. IntrinsicQ values, at 800 Hz, were on the order of 140, similar to quality factor values obtained in previous studies using Spectral Decay and Coda-Q methods (120 to 170). The scattering quality factor at this frequency was about 520. An observed frequency dependence of the scattering attenuation suggested that a decrease in the density of scatterers, with scale lengths on the order of 2 m, exists at the site. Characteristic fracture scale lengths were considered to range from 4 to 6 m as identified in the mapping data. These observations were supported by the increase in scattering found for seismic waves with frequencies less than 1000 Hz. By assuming that the identified scatters are characteristic faults, these scatterers can then be considered to increase nonsimilar behavior in source scaling. Overall, our results suggest that MLTWA provides a practical method for remotely characterizing the quality of a rock mass when visual observations are not attainable.  相似文献   

5.
The physical meaning of the characteristic displacement that has been observed in velocity-stepping friction experiments was investigated based on the micromechanics of asperity contact. It has been empirically found for bare rock surfaces that the magnitude of the characteristic displacement is dependent only on surface roughness and insensitive to both slip velocity and normal stress. Thus the characteristic displacement has been interpreted as the displacement required to change the population of contact points completely. Here arises a question about the physical mechanism by which the contact population changes. Because individual asperity contacts form, grow and are eliminated with displacement, there are at least two possible interpretations for the characteristic displacement: (1) it is the distance over which the contacts existing at the moment of the velocity change all fade away, being replaced by new asperity contacts, or (2) it is the distance required for a complete replacement in the real contact area that existed at the moment of the velocity change. In order to test these possibilities, theoretical models were developed based on the statistics of distributed asperity summits. A computer simulation was also performed to check the validity of the theoretical models using three-dimensional surface topography data with various surface roughnesses. The deformation was assumed to be elastic at each asperity contact. The results of both the simulation and the theoretical models show that the characteristic displacement in (1) is about three times longer than that in (2). Comparison of the results with the experimental observations obtained by others indicates that the possibility (2) is the correct interpretation. This means that the state in the rate and state variable friction law is memorized in a very confined area of real contact. Further, our results explain why the characteristic displacement is insensitive to normal stress: this comes from the fact that the microscopic properties such as the mean contact diameter are insensitive to normal stress. The approach based on the micromechanics of asperity contact is useful to investigate the underlying mechanism of various phenomena in rock friction.  相似文献   

6.
A new model for Plinian eruption columns is derived from first principles and investigated numerically. The dynamics particular to the momentum-driven basal gas-thrust region and the upper buoyancy-driven convective region are treated separately. The thermal interactions in the column are modelled by the steady-flow-energy equation. The main results of the present paper are that: (1) the basal gas-thrust region model predicts a very rapid initial expansion of the material on leaving the vent; (2) the gas-thrust region height decreases with initial temperature, inital gas content of the erupted material and initial velocity, but increases with vent radius; (3) the total column height increases with initial temperature, initial velocity and vent radius, but decreases with initial gas content; (4) column collapse occurs for initial velocities of the order of 100 m/s; the precise value increases as the initial gas content in the erupted material decreases; (5) for large vent radii or low initial gas content of the erupted material, the velocity in the column can increase with height once in the buoyancy-driven region instead of decaying to zero monotonically; (6) the interaction of the potential energy with the enthalpy is found to be the dominant thermal interaction in the upper part of the column. Previous models of eruption columns involve inconsistencies and simplifications; these are shown to lead to significant differences in the results in comparison to the present model.  相似文献   

7.
— We have used numerical simulations with the boundary integral equation method to investigate a mechanism to excite super-shear rupture velocities in a homogeneous stress field including an asperity of increased initial stress. When the rupture, with the slip-weakening distance selected to generate sub-Rayleigh speed, encounters the asperity it either accelerates to super-shear velocities or maintains the sub-Rayleigh speed, dependent on the size and amplitude of the asperity. Three classes of rupture propagation are identified: the velocity (a) for the most narrow asperities increases slowly towards the Rayleigh wave speed, (b) for intermediate width of the asperities jumps to super-shear values for a short distance but then decreases to sub-Rayleigh wave speeds, and (c) for the widest asperities jumps to super-shear values and pertains to values between the S- and P-wave velocities. The transitions between the three classes of rupture propagation are characterized by very narrow (critical) ranges of rupture resistance. If the size of the initial asperity is smaller than critical, it becomes difficult for rupture to propagate with super-shear velocities even if the initial stress level is high. Our results suggest that stress variation along the rupture path helps homogenize the rupture velocity and propagate with sub-Rayleigh wave speeds.  相似文献   

8.
The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.  相似文献   

9.
We characterize the heterogeneous source slip model of intraslab earthquakes to compare source scaling properties with those of inland crustal and subduction-zone plate-boundary earthquakes. We extracted rupture area (S), total area of asperity (S a), average slip (D) and average slip on asperity (D a) of eleven intraslab earthquakes following the procedure proposed by Somerville et al. (Seism Res Lett 70:59?C80, 1999) and proposed the empirical scaling relationship formula of S, S a, and D for intraslab earthquakes. Under the same seismic moment, an intraslab earthquake has a smaller rupture area and total area of asperity, and smaller average slip than an inland crustal earthquake. The area ratio of asperity area and total rupture area of intraslab earthquakes are similar to those of inland crustal earthquakes. The strong motion generation area (SMGA) scaling of intraslab earthquakes appears self-similar, and those results support the idea the characterized source model of intraslab earthquakes can be modeled in a manner similar to that of inland crustal earthquakes.  相似文献   

10.
Strong motion (SM) data of six Mexican subduction zone earthquakes (6.4M S8.1) recorded near the epicentral zone are analyzed to estimate their far-field source acceleration spectra at higher frequencies (f0.3 Hz). Apart from the usual corrections such as geometrical spreading (1/R), average radiation pattern (0.6), free surface amplification (a factor of 2), and equal partitioning of the energy into two orthogonal horizontal components (a factor of 1/ ), the observed spectra are corrected for a frequency dependentQ(Q=100f), a site dependent filter (e –kf ), and amplification ofS waves near the surface (a factor of about 2 atf2Hz). We takeR as the average distance from the rupture area to the site. If we model the high frequency plateau (f1 Hz) of the source spectra, by a point source –2-model, and interpret them in terms of Brune's model we obtain between 50 and 100 bars for all earthquakes. The low-frequency broadband teleseismicP wave spectra, corrected witht *=1.0 s, agrees within a factor of two with SM source spectra near 1 Hz. The –2-model is inadequate to explain the observed source spectra in a broad frequency range; these resemble spectra given byGusev (1983) with some differences.SM source acceleration spectra require significant corrections to explain observed spectra and RMS acceleration (arms) (a) at farther coastal sites for extended sources due to directivity effect and (b) at inland sites (100R200 km) because of unaccounted path and site amplification and/or invalidity of body-wave approximation. The observed spectra and arms at these sites are significantly greater than the predicted values from the estimated source spectra.  相似文献   

11.
Excavation induced seismic events with moment magnitudesM<0 are examined in an attempt to determine the role geology, excavation geometry, and stress have on scaling relations. Correlations are established based on accurate measurements of excavation geometry and methodology, stress regime, rock mass structure, local tectonics, and seismic locations. Scaling relations incorporated seismic moments and source radii obtained by spectral analysis, accounting for source, propagation, and site effects, and using Madariaga's dynamic circular fault model. Observations suggest that the interaction of stresses with pre-existing fractures, fracture complexity and depth of events are the main factors influencing source characteristics and scaling behaviour. Self-similar relationships were found for events at similar depths or for weakly structured rock masses with reduced clamping stresses, whereas a non-similar behaviour was found for events with increasing depth or for heavily fractured zones under stress confinement. Additionally, the scaling behaviour for combined data sets tended to mask the non-similar trends. Overall, depth and fracture complexity, initially thought as second order effects, appear to significantly influence source characteristics of seismic events withM<0 and consequently favour a non-similar earthquake generation process.  相似文献   

12.
The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.  相似文献   

13.
A model has been developed to simulate the statistical and mechanical nature of rupture on a heterogeneous strike-slip fault. The model is based on the progressive failure of circular asperities of varying sizes and strengths along a fault plane subjected to a constant far-field shear displacement rate. The basis of the model is a deformation and stress intensity factory solution for a single circular asperity under a unidirectional shear stress. The individual asperities are unified through the fault stiffness and the far-field stress and displacement. During fault deformation asperities can fail and reheal, resulting in changes in the local stresses in the asperities, stress drops, and changes in the stiffness of the fault. Depending on how the stress is redistributed following asperity failure and on the strenghts of the neighboring asperities an earthquake event can be the failure of one or more asperities. Following an earthquake event seismic source parameters such as the stress drop, energy change, and moment magnitude are calculated. Results from the model show a very realistic pattern of earthquake rupture, with reasonable source parameters, the proper magnitude-frequency behavior, and the development of characteristic earthquakes. Also the progression ofb-values in the model gives some insight into the phenomenon of self-organized criticality.  相似文献   

14.
A two-stream approximation is applied to the equation of infrared (IR) radiative transfer in order to compute the detailed structure of fluxes and cooling/heating rates in inhomogeneous, cloudy atmospheres. The spectrum between 4 and 400 m is split into 50 bands, and absorption of water vapour, uniformly mixed gases, ozone, water vapour polymers and water droplets are taken into account. Allowance is made for scattering by water droplets in the atmospheric window region (8–12.8 m) where the delta function approximation is used. The backward scattering coefficient is calculated from the asymmetry parameter by means of a new simple formula. The model is compared to other radiative transfer schemes and some applications are presented.  相似文献   

15.
Seismological models for mining-induced seismic events   总被引：1，自引：0，他引：1  
  相似文献   

16.
The paper is devoted to estimating the dynamical parameters of 14 earthquakes with intermediate magnitudes (energy class 11 to 14), which occurred in the Northern Tien Shan. For obtaining the estimates of these parameters, including the stress drop, which could be then applied in crustal stress reconstruction by the technique suggested by Yu.L. Rebetsky (Schmidt Institute of Physics of the Earth, Russian Academy of Sciences), we have improved the algorithms and programs for calculating the spectra of the seismograms. The updated products allow for the site responses and spectral transformations during the propagation of seismic waves through the medium (the effect of finite Q-factor). By applying the new approach to the analysis of seismograms recorded by the seismic KNET network, we calculated the radii of the sources (Brune radius), scalar seismic moment, and stress drop (release) for the studied 14 earthquakes. The analysis revealed a scatter in the source radii and stress drop even among the earthquakes that have almost identical energy classes. The stress drop by different earthquakes ranges from one to 75 bar. We have also determined the focal mechanisms and stress regime of the Earth’s crust. It is worth noting that during the considered period, strong seismic events with energy class above 14 were absent within the segment covered by the KNET stations.  相似文献   

17.
—The values of uniformly estimated apparent stress σa and Brune’s stress drop, taken as a measure of static stress drop, from 850 seismic events, with moment magnitude ranging from ?3.6 to 3.6, induced at the Underground Research Laboratory (URL) in Canada, Western Deep Levels (WDL) gold mine in South Africa, and two coal and two copper mines in Poland, were collected to study the stress release mode in various mining environments. For this, the quantity epsilon, ?=Δσ/(σa + Δσ/2), where Δσ is the static stress drop, proposed by Zúñiga (1993) as an indicator of stress-drop mechanism was used. The events induced at the URL are characterized by low values of epsilon corresponding to a partial stress drop mechanism, whereas all the events at WDL display a frictional overshoot mechanism in which final stress reaches a lower value than that of frictional stress. The events at Polish coal and copper mines, on the other hand, are in good agreement with the well-known Orowan’s condition such that the final stress is equal to the dynamic frictional stress.¶The Brune stress drop, however, is heavily model dependent through the source radius-corner frequency relation. The Orowan’s condition for the events from the URL would be met if a constant in the source radius-corner frequency relation is equal to 1.82 ± 0.12, and for the events from WDL if it is equal to 3.92 ± 0.40, in contrast to Brune’s constant of 2.34.¶The smoothed values of epsilon displayed as a function of time, represented by the consecutive event numbers, for selected sets of events imply that the largest seismic events in a given set occur when the epsilon is low and a partial stress drop mechanism is dominant. The large events are then followed by high epsilon values when a frictional overshoot mechanism begins to dominate.  相似文献   

18.
Forty-six mining-induced seismic events with moment magnitude between ?1.2 and 2.1 that possibly caused damage were studied. The events occurred between 2008 and 2013 at mining level 850–1350 m in the Kiirunavaara Mine (Sweden). Hypocenter locations were refined using from 6 to 130 sensors at distances of up to 1400 m. The source parameters of the events were re-estimated using spectral analysis with a standard Brune model (slope ?2). The radiated energy for the studied events varied from 4.7 × 10?1 to 3.8 × 107 J, the source radii from 4 to 110 m, the apparent stress from 6.2 × 102 to 1.1 × 106 Pa, energy ratio (E s/E p) from 1.2 to 126, and apparent volume from 1.8 × 103 to 1.1 × 107 m3. 90% of the events were located in the footwall, close to the ore contact. The events were classified as shear/fault slip (FS) or non-shear (NS) based on the E s/E p ratio (>10 or <10). Out of 46 events 15 events were classified as NS located almost in the whole range between 840 and 1360 m, including many events below the production. The rest 31 FS events were concentrated mostly around the production levels and slightly below them. The relationships between some source parameters and seismic moment/moment magnitude showed dependence on the type of the source mechanism. The energy and the apparent stress were found to be three times larger for FS events than for NS events.  相似文献   

19.
On 29–30 September 1996, an earthquake sequence occurred in the Füzesgyarmat region in eastern Hungary. The main shock had a magnitude of ML = 3.2 and was felt with a maximum intensity 4 MSK. It was preceded by a foreshock with a magnitude of ML = 2.8 and was followed, within six hours, by five aftershocks with magnitudes 2.1ML 3.1. The dynamic source parameters of the Füzesgyarmat earthquake sequence have been derived from P-wave spectra of the Hungarian seismograph stations. The average of the obtained values at different stations shows that the main shock occurred on a fault length of 610 m, with relative displacement of 1.13 cm, stress drop of 7 bar and seismic moment of 3.96*1021 dyne.cm. The main shock was small to yield data for a full mechanism solution and no reliable single fault plane solution could be obtained due to the low signal to noise ratio at the recording stations. The parameters of the foreshock are fault length of 560m, seismic moment of 2.09*1021 dyne.cm, stress drop of 5.53 bar and relative displacement of 0.73 cm. The five aftershocks show source parameters similar to the foreshock stress drops (5.26 5.76 bar), fault lengths (415 L 600 m), seismic moments (8.36*1020 Mo 2.31*1021 dyne.cm) and relative displacements (0.52 0.91 cm).  相似文献   

20.
Among the second-order effects on friction the most important are those of variable normal stress and of slip velocity. Velocity weakening, which is usually considered the source of the stick-slip instability in rock friction, has been observed in velocity stepping experiments with Westerly granite. The friction change, , was –0.01 to –0.008 for a tenfold velocity increase. Using normal closure measurements, we observed dilation upon each increase in sliding rate. We also observed, for the first time, time-dependent closure between surfaces during static loading. The dilation that occurred during the velocity stepping experiment was found to be that expected from the static time-dependent closure phenomenon. This change in closure was used to predict friction change with an elastic contact model. The calculated friction change which results from a change in contact area and asperity interlocking, is in good agreement with the observed velocity dependence of steady-state friction. Variable normal stress during sliding has two effects, first in creating new partial slip contacts and locking some existing fully sliding contacts and second in increasing interlocking, for instance when normal load is suddenly increased. As a result, a transient change in friction occurs upon a sudden change in normal load.  相似文献   

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

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