共查询到20条相似文献，搜索用时 15 毫秒
1.
岩石破裂滑动过程及其声发射数值模拟 总被引：3，自引：0，他引：3
运用岩石破裂与失稳过程分析RFPA2D系统, 对滑动形成过程进行了数值模拟研究。滑动岩体似为非均匀弹 -脆岩石材料, 模拟结果再现了滑动形成从变形到破坏直至失稳的全过程及其声发射规律。 相似文献
2.
We have used a combination of macroseismic and instrumental techniques to re-evaluate the seismicity of Cyprus and the surrounding region for the period 1890–1900. We identified 56 events which we believe to have a surface-wave magnitude (Ms) of 4.5 or greater. Early events are located by felt effects alone, and macroseismic information is used to improve the location of many later events. We find the broad pattern of seismicity follows that given by standard listings: an area of weak seismicity to the north-west of Cyprus in the Gulf of Antalia includes intermediate-depth events as are found in the Hellenic Arc to the west; a gap of reduced activity separates this activity from that of the main region of Cyprus, where we could find no evidence for sub-crustal foci; there appears to be no present seismological connection along the presumed plate boundary to the north-east to the junction of the East Anatolian fault and the extension of the Dead Sea rift system. We find a much greater level of activity in the first part of this century than is given in the standard listings, such as that of the ISC. Up to 1963 we identified 47 events of Ms > 4.5, compared with only 10 given in the ISC files. Recent improved station coverage locates more events in the region, but these have not been of large magnitude, and the moment release in recent decades has been appreciably less than earlier in the century. For the events in our list the total moment release is 2.7x 1019 Nm with 75% of this originating in the main Cyprus region. We reaffirm the danger of uncritical use of standard earthquake listings for tectonic interpretation and in particular for hazard evaluation in areas such as this, where neither the historical nor the instrumental seismicity record is complete. 相似文献
3.
随着页岩气开采、废水回注和CO2地质封存工程活动的进行,储层应力环境改变诱发地震的问题得到广泛关注,研究储层岩石变形破裂和声发射规律对于理解诱发地震活动具有重要意义。本文选取了页岩、致密砂岩和白云岩开展三轴压缩声发射试验,获取了空间裂缝形态,揭示了不同类型储层岩石的变形破裂过程和声发射特征。结果表明：(1)页岩的扩容应力与峰值应力比最高,其次为白云岩和致密砂岩,表明页岩主要发生脆性破坏。(2)岩石结构对裂缝扩展和强度具有明显影响,页岩层理发育,抗压强度和裂缝形态各向异性显著。当层理角度为0°时,试样发生剪切-拉张复合型破坏。当层理角度为30°和60°时,试样主要发生剪切破坏。当层理角度增大至90°时,试样主要发生拉张破坏。致密砂岩破裂形成剪切主裂缝,白云岩形成两条剪切主裂缝和微裂缝。(3)不同类型储层岩石破裂过程声发射特征差异显著。页岩在扩容应力点附近有少量声发射活动,达到峰值应力时,声发射活动迅速增强。致密砂岩仅在破裂瞬间有少量声发射活动。相比之下,白云岩在裂缝非稳定扩展阶段和峰后阶段声发射活动显著。因此,在工程实践中需要根据微震监测调整施工措施,避免页岩作为储... 相似文献
4.
Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D2 and D15 of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D2 and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture. 相似文献
5.
Cellular automata are simple mathematical idealizations of natural systems and they supply useful models for many investigations in natural science. Examples include sandpile models, forest fire models, and slider block models used in seismology. In the present paper, they have been used for establishing temporal relations between the energy releases of the seismic events that occurred in neighboring parts of the crust. The catalogue is divided into time intervals, and the region is divided into cells which are declared active or inactive by means of a threshold energy release criterion. Thus, a pattern of active and inactive cells which evolves over time is determined. A stochastic cellular automaton is constructed starting with these patterns, in order to simulate their spatio-temporal evolution, by supposing a Moore's neighborhood interaction between the cells. The best model is chosen by maximizing the mutual information between the past and the future states. Finally, a Probabilistic Seismic Hazard Map is given for the different energy releases considered. The method has been applied to the Greece catalogue from 1900 to 1999. The Probabilistic Seismic Hazard Maps for energies corresponding to m = 4 and m = 5 are close to the real seismicity after the data in that area, and they correspond to a background seismicity in the whole area. This background seismicity seems to cover the whole area in periods of around 25–50 years. The optimum cell size is in agreement with other studies; for m > 6 the optimum area increases according to the threshold of clear spatial resolution, and the active cells are not so clustered. The results are coherent with other hazard studies in the zone and with the seismicity recorded after the data set, as well as provide an interaction model which points out the large scale nature of the earthquake occurrence. 相似文献
6.
The shear wave attenuation field in the lithosphere of Eastern Tien Shan has been mapped. The method based on analysis of
the ratio between amplitudes of Sn and Pn waves was used. On aggregate, about 120 seismograms made at the Makanchi station (MKAR), mainly in the period of 2003–2009,
at epicentral distances of about 350–1200 km were analyzed. It was found that shear wave attenuation in the lithosphere of
Eastern Tien Shan is weaker than that in the region of Central Tien Shan. This agrees with the fact that the rate of deformation
of the Earth’s crust in Eastern Tien Shan is lower (based on GPS data), as is the seismicity level, in comparison to Central
Tien Shan. The zones of high attenuation, where strong earthquakes with M > 7.0 have not occurred for the last 200 years, have been identified: first of all, these are the area west of Urumqi and
that of the Lop Nur test site. It is suggested that in the first zone, where an annular seismicity structure has formed over
the last 30 years, a strong earthquake may be being prepared. The second zone is most probably related to the uplift of mantle
fluids resulting from a long-term intensive technogenic effect, analogous to what has occurred in areas of other nuclear test
sites (Nevada and Semipalatinsk). 相似文献
7.
The intra-plate seismicity map for southwest Japan, based on fairly complete historical data for the past four hundred years, reveals an inverse correlation between the seismic activity along the island arc and the slip-rate along the Median Tectonic Line during the Late Quaternary. In the eastern part, the tectonic line is geologically inactive but regional historic seismicity has been high. The intra-plate seismic activity is probably related with the well developed mosaic-like conjugate system of strike-slip faults there. Conversely, the historic seismicity has been low in the western part, especially low in an area along the most geologically active segment of the Median Tectonic Line. Since no creep movement has been found there, energy greater than that of the Mino-Owari earthquake of 1891 (M = 8.0) seems to be stored in this seismicity gap. The difference in seismic released energy between the two regions for the last four hundred years would be balanced by the strain energy accumulated in the seismic gap. The fairly uniform strain release is conformable to the idea, proposed on the basis of the trend of maximum compression axes, that the Philippine Sea plate is dragging southwest Japan southwestward along the Nankai trough. 相似文献
8.
We use a continuum model of reservoir evolution to explore the interaction of coupled thermal, hydraulic and chemical processes that influence the evolution of seismicity within a fractured reservoir from stimulation to production. Events occur from energy release on seeded fractures enabling moment magnitude, frequency and spatial distribution to be determined with time. Event magnitudes vary in the range ?2 to +2 with the largest event size (~2) corresponding to the largest fracture size (~500 m) and a prescribed stress drop of 9 MPa. Modelled b‐values (~0.6–0.7) also correspond to observations (~0.7–0.8) for response in the Cooper Basin (Australia). We track the hydrodynamic and thermal fronts to define causality in the triggering of seismicity. The hydrodynamic front moves twice as fast as the thermal front and envelops the triggered seismicity at early time (days to month) – with higher flow rates correlating with larger magnitude events. For later time (month to years), thermal drawdown and potentially chemical influences principally trigger the seismicity, but result in a reduction in both the number of events and their magnitudes. 相似文献
9.
We have analysed the earthquake sequence occurred at Campi Flegrei during an unrest episode of strong ground uplift and seismicity, occurred in the period 1982–1984. The maximum magnitude of these events was about 4. Both earthquake occurrence and ground deformation have been interpreted in terms of the role played by a ring fault system, inward dipping, related to phenomena of caldera collapse and resurgence. Earthquakes are of mixed strike-slip and normal fault type. They show a dip movement opposite to the static ground deformation. The rising of the internal block with respect to the zone external to the ring fault, as observed by ground deformation, should cause thrust fault slip on the fracture system, whereas a normal fault dip component is observed. The simulation of the stress field generated by overpressure in a magma chamber in presence of lateral discontinuities, as performed by a boundary element method, allows to hypothesise that reverse fault slip on the ring fault is mainly aseismic, and such aseismic movement is able to focus normal fault shear stress along the lateral discontinuities. Aseismic slip on the ring fault in response to static deformation is also supported by the low amount of seismic moment released (M0 ≅ 1015 Nm), about two orders of magnitude lower than expected from the shear slip on the discontinuities needed to accomplish the total static surface deformation (1.8 m). Such results have been compared with observations at Rabaul caldera, during a similar unrest episode. In such area, the seismic moment release is in good agreement with shear slip produced on a system of outward dipping ring faults, and seismicity is much more focused on the fault structures. The comparison between the two areas shed new light about the dynamics of earthquakes in calderas, as due to the role of bordering ring fault systems. 相似文献
10.
An analysis of the vane test using an Arbitrary Lagrangian–Eulerian formulation within a finite element framework is presented. This is suitable for soft clays for which the test is commonly used to measure in situ undrained shear strength. Constitutive laws are expressed in terms of shear stress–shear strain rates, and that permits the study of time effects in a natural manner. An analysis of the shear stress distributions on the failure surface according to the material model is presented. The effect of the constitutive law on the shear band amplitude and on the position of the failure surface is shown. In general, the failure surface is found at 1–1·01 times the vane radius, which is consistent with some experimental results. The problem depends on two dimensionless parameters that represent inertial and viscous forces. For usual vane tests, viscous forces are predominant, and the measured shear strength depends mainly on the angular velocity applied. That can explain some of the comparisons reported when using different vane sizes. Finally, the range of the shear strain rate applied to the soil is shown to be fundamental when comparing experimental results from vane, triaxial and viscosimeter tests. Appart from that, an experimental relation between undrained shear strength and vane angular velocity has been reproduced by this simulation. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
11.
Summary A number of earthquakes have been recorded in strata above underground potash mines in Saskatchewan, Canada. These seismic events are widely understood to be generated in the carbonate Dawson Bay Formation and, possibly, other superincumbent beds, and are believed to be directly related to mining operations. The mechanical response of these higher strata to mining is not observable andin situ instrumentation for measuring post-mining stresses and strains in the carbonate beds is sparse. Numerical models are believed to be the most effective option for investigating the response of these higher strata to potash mining. In this regard a realistic finite element model based on elasto-plastic material behaviour is developed in this paper. Special capabilities added to this general elasto-plastic model permit simulation of the post-failure behaviour of the rock mass. Possible modes of failure in the Dawson Bay are investigated and the brittle failure of intact limestone is affirmed as the possible mode of failure producing major seismicity. Slip along existing discontinuities in overlying strata is also demonstrated to be another possible mechanism, though with lower energy levels. 相似文献
12.
Deformation within the Earth's lithosphere is largely controlled by the rheology of the rock. Fracture and faulting are characterized by elastic rheologies with brittle mechanisms, while folding and flow are characterized by plastic and/or viscous rheologies due to ductile mechanisms. However, it has been recognized that deformation that resembles ductile behavior can be produced within the confines of the brittle lithosphere. Specific examples are folds that form in the shallow crust, steep hinges at subduction zones that are accompanied by seismicity, and large-scale deformation at plate boundaries. In these cases, the brittle lithosphere behaves elastically with fracture and faulting yet produces ductile behavior. In this paper, we attempt to simulate such ductile behavior in elastic materials using continuum damage mechanics. Engineers utilize damage mechanics to model the continuum deformation of brittle materials. We utilize a modified form of damage mechanics that represents a reduction in frictional strength of preexisting fractures and faults. We use this empirical approach to simulate the bending of the lithosphere under the application of a constant moment.We use numerical simulations to obtain elastostatic solutions for plate bending and where the longitudinal stress at a particular node exceeds a yield stress, we apply damage to reduce Young's modulus at the node. Damage is calculated at each time step by a power-law relationship of the ratio of the yield stress to the longitudinal stress and the yield strain to the longitudinal strain. This results in the relaxation of the material due to increasing damage. To test our method, we apply our damage rheology to an infinite plate deforming under a constant bending moment. We simulate a wide range of behaviors from slow relaxation to instantaneous failure, over timescales that span six orders of magnitude. Using this method, stress relaxation produces elastic-perfectly plastic behavior in cases where failure does not occur. For cases of failure, we observe a rapid increase in damage leading to failure, analogous to the acceleration of microcrack formation and acoustic emissions prior to failure. The changes in the rate of damage accumulation in failure cases are similar to the changes in b-values of acoustic emissions observed in triaxial compression tests of fractured rock and b-value changes prior to some large earthquakes. Thus continuum damage mechanics can simulate the phenomenon of ductile behavior due to brittle mechanisms as well as observations of laboratory experiments and seismicity. 相似文献
13.
吉塞尔斯地震活动可能是因蒸汽开发引起的.水在一个承受很高构造剪应力和应变的大范围破裂体中汽化为蒸汽。 汽田地震震源机制解与区域构造应变场几乎一致,并且在该区域范围内汽田地震与别的构造地震不易区分。观测中注意到地震活动与注液历史无关,这表明孔隙水压力增高与注液不可能是吉塞尔斯诱发地震的成因。 相反,所有证据都表明诱发地震与孔隙水压力及温度降低有关。形成机制有两种最大可能:其一,是裂隙排水(汽)导致局部剪应力增加所致,其二,是由稳定滑动转化为不稳定滑动(粘滑)。没有其它记载的诱发地震机制与吉塞尔斯汽田条件相符。 相似文献
14.
Moment tensors and micromechanical models 总被引：4，自引：0，他引：4
A numerical modelling approach that simulates cracking and failure in rock and the associated seismicity is presented and a technique is described for quantifying the seismic source mechanisms of the modelled events. The modelling approach represents rock as an assemblage of circular particles bonded together at points of contact. The connecting bonds can break under applied stress forming cracks and fractures in the modelled rock. If numerical damping is set to reproduce realistic levels of attenuation, then energy is released when the bonds break and seismic source information can be obtained as damage occurs. A technique is described by which moment tensors and moment magnitudes can be calculated for these simulated seismic events. The technique basically involves integrating around the source and summing the components of force change at the surrounding particle contacts to obtain the elements of the moment tensor matrix. The moment magnitude is then calculated from the eigenvalues of the moment tensor. The modelling approach is tested by simulating a well-controlled experiment in which a tunnel is excavated in highly stressed granite while microseismicity is recorded. The seismicity produced by the model is compared to the actual recorded seismicity underground. The model reproduces the spatial and temporal distribution of seismicity observed around the tunnel and also the magnitudes of the events. A direct comparison between the actual and simulated moment tensors is not performed due to the two-dimensional nature of the model, however, qualitative comparisons are presented and it is shown that the model produces intuitively realistic source mechanisms. The ability to obtain seismic source information from the models provides a unique means for model validation through comparison with actual recorded seismicity. Once it is established that the model is performing in a realistic manner, it can then be used to examine the micromechanics of cracking, failure and the associated seismicity and to help resolve the non-uniqueness of the geophysical interpretation. This is demonstrated by examining in detail the mechanics of one of the modelled seismic events by observation of the time dependence of the moment tensor and by direct examination of the particle motions at the simulated source. 相似文献
15.
Short-term earthquake prediction, months in advance, is an elusive goal of earth sciences, of great importance for fundamental science and for disaster preparedness. Here, we describe a methodology for short-term prediction named RTP (Reverse Tracing of Precursors). Using this methodology the San Simeon earthquake in Central California (magnitude 6.5, Dec. 22, 2003) and the Tokachi-Oki earthquake in Northern Japan (magnitude 8.1, Sept. 25, 2003) were predicted 6 and 7 months in advance, respectively. The physical basis of RTP can be summed up as follows: An earthquake is generated by two interacting processes in a fault network: an accumulation of energy that the earthquake will release and a rise of instability triggering this release. Energy is carried by the stress field, instability is carried by the difference between the stress and strength fields. Both processes can be detected and characterized by “premonitory” patterns of seismicity or other relevant fields. Here, we consider an ensemble of premonitory seismicity patterns. RTP methodology is able to reconstruct these patterns by tracing their sequence backwards in time. The principles of RTP are not specific to earthquakes and may be applicable to critical transitions in a wide class of hierarchical non-linear systems. 相似文献
16.
The focal mechanisms for 86 selected earthquakes (3.0 mb 5.5) located in central Alaska have been investigated from P-wave first motions; the data were gathered by local seismic networks. The results show a depth-dependent characteristic to the fault-plane solutions. For earthquakes having focal depths shallower than 60–70 km, the focal mechanisms indicate either strike-slip or normal faults, while for earthquakes with foci at intermediate depths the focal mechanisms correspond to thrust faults. The nature of the seismicity indicates the hinge line of the Pacific lithospheric plate under the study area to be striking N17°E from Cook Inlet towards interior Alaska. The comparison of the focal mechanisms with the seismicity shows that the strike-slip and normal faults are the predominant processes of stress release along the shallow section of the plate. The earthquakes with intermediate foci systematically occur along the inclined section of the plate. If the gently dipping nodal planes for these earthquakes are chosen as the fault planes, the focal mechanisms correspond to underthrust motions at the foci. In these, the slip vectors are oriented either to the west or north with the resultant being in the N30°W direction. The tension axes for the underthrust solutions are also found to be parallel to the local dip of the plate, indicating that the subducted plate in interior Alaska is undergoing gravitational sinking. 相似文献
17.
Consolidation theory and concepts of rock failure can be used to evaluate the probable risk of induced seismicity as a result of filling of reservoirs. This evaluation indicates the safest way to fill a reservoir, and depends only on the geometry of the load, the rate of filling and the geological structures in the area. The stability function is actually a measure of the risk of having failure, with time, for a particular loading history in respect to a plane of weakness.
The stability function is applied to the area of the Itzantun reservoir, which will be in southern Mexico. Drawdowns can increase the risk of triggering earthquakes in this area, which is prone to thrust faulting. It is possible to estimate the stresses after a period during which the water level is maintained and a decrease in stresses with the depth of the observation point.
The estimates of the probable induced seismicity are limited as the residual stress in the area prior to the impounding is unknown. With a measure of the residual tectonic stress it will be possible to determine an optimal filling rate to reduce the probability of induced seismicity. 相似文献
18.
An examination of the Alaskan earthquake catalogs revealed that from 1928 through 1965, twelve earthquakes of magnitudes (M) in the range 5.6 to 7.3 were located in and around Seward Peninsula region of western Alaska. Moreover, a number of earthquakes of M < 5.0 were found to distribute over the same area. Further investigation of the seismicity employing a local seismographic network in the above area showed a higher level of onshore and offshore seismic activity than had been previously recognized. A number of clusters of earthquakes have been identified. Some of them are located in the epicentral areas of past strong earthquakes (M > 5.5) and some are located in areas traversed by mapped faults. The nature of the seismicity as identified with the local network data is primarily crustal over the entire study area. Investigation of focal mechanisms of isolated strong earthquakes or clusters of small earthquakes show normal faulting as the dominant mode of strain energy release in the western part of Alaska. Moreover, in areas lying, approximately, south and north of Kotzebue Sound, the principal components of horizontal stresses tend to orient in the NW-SE and WNW-ESE directions, respectively. 相似文献
19.
V.V. Ruzhich S.G. Psakhie E.N. Chernykh S.A. Bornyakov N.G. Granin 《Russian Geology and Geophysics》2009,50(3):214-221
The mechanics of the ice cover of Lake Baikal has been studied through monitoring of its deformation and seismic effects and full-sizeuniaxial compression and shear tests in 2005–2007. We measured the shear strength of ice specimens and large in situ blocks (σ = 0.2?1.9 MPa) and investigated it as a function of air temperature and ice structure. Deformation was analyzed in terms of various natural controls, such as air temperature and pressure, wind, sub-ice currents, and local earthquakes. Precise strain measurements along ice cracks were used to explore the strain behavior of ice, including the cases of dynamic failure (ice shocks). Measurements by seismic station Baikal-12 were used to monitor diurnal background micro seismicity variations and to record an ice quake with its magnitude (M = 0.3–0.8; E = 104–105 J) comparable to a medium-size rock burst or a small earthquake. Ice quakes were studied in terms of their nucleation, dynamics, and after effects, as well as the strain and seismic responses of the ice, using sub-ice explosions in the latter case. The natural conditions of deformation in the elastoviscoplastic Baikal ice are similar to lithospheric processes and thus can be employed in tectonophysical modeling with scientific and practical implications for hazard mitigation. 相似文献
20.
The Himalayas has experienced varying rates of earthquake occurrence in the past in its seismo-tectonically distinguished segments which may be attributed to different physical processes of accumulation of stress and its release, and due diligence is required for its inclusion for working out the seismic hazard. The present paper intends to revisit the various earthquake occurrence models applied to Himalayas and examines it in the light of recent damaging earthquakes in Himalayan belt. Due to discordant seismicity of Himalayas, three types of regions have been considered to estimate larger return period events. The regions selected are (1) the North-West Himalayan Fold and Thrust Belt which is seismically very active, (2) the Garhwal Himalaya which has never experienced large earthquake although sufficient stress exists and (3) the Nepal region which is very seismically active region due to unlocked rupture and frequently experienced large earthquake events. The seismicity parameters have been revisited using two earthquake recurrence models namely constant seismicity and constant moment release. For constant moment release model, the strain rates have been derived from global strain rate model and are converted into seismic moment of earthquake events considering the geometry of the finite source and the rates being consumed fully by the contemporary seismicity. Probability of earthquake occurrence with time has been estimated for each region using both models and compared assuming Poissonian distribution. The results show that seismicity for North-West region is observed to be relatively less when estimated using constant seismicity model which implies that either the occupied accumulated stress is not being unconfined in the form of earthquakes or the compiled earthquake catalogue is insufficient. Similar trend has been observed for seismic gap area but with lesser difference reported from both methods. However, for the Nepal region, the estimated seismicity by the two methods has been found to be relatively less when estimated using constant moment release model which implies that in the Nepal region, accumulated strain is releasing in the form of large earthquake occurrence event. The partial release in second event of May 2015 of similar size shows that the physical process is trying to release the energy with large earthquake event. If it would have been in other regions like that of seismic gap region, the fault may not have released the energy and may be inviting even bigger event in future. It is, therefore, necessary to look into the seismicity from strain rates also for its due interpretation in terms of predicting the seismic hazard in various segments of Himalayas. 相似文献