首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 865 毫秒
1.
Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.  相似文献   

2.
Seismic potential of Southern Italy   总被引：1，自引：2，他引：1  
To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.  相似文献   

3.
V.I. German   《Tectonophysics》2006,424(3-4):167
The paper describes the unified scaling theory for distribution functions of temporal and spatial characteristics in seismology. It is based on the scaling of seismological characteristics calculated for various energy–spatial–temporal intervals. The common mathematical methods for the scaling of distribution functions are developed. The means to test possibility of such scaling are found as well. The relationship between the unified scaling theory and other present scaling approaches is determined. The theory is applied to two characteristics of different seismoactive regions. The first characteristic is the waiting time between earthquakes ΔT, the second one is a new space parameter ΔDmin, which is the minimum distance of a current seismic event to the nearest (in space) neighbor in an energy–spatial–temporal interval. The distribution of the characteristics ΔT and ΔDmin allows estimating the time interval to the next earthquake and the distance of the following earthquake from previous earthquakes. Thus, these characteristics are very important for seismic hazard estimations. Scaling of distributions functions is proven to be successful for ΔDmin in all energy–spatial–temporal intervals and for ΔT with variations of energy/magnitude range. The distribution function of ΔT for various time domains was stable in only 60% of the cases, and near to unstable for spatial variations.  相似文献   

4.
  相似文献   

5.
The EEPAS (“Every Earthquake a Precursor According to Scale”) model is a space–time point-process model based on the precursory scale increase (Ψ) phenomenon and associated predictive scaling relations. It has previously been fitted to the New Zealand earthquake catalogue, and applied successfully in quasi-prospective tests on the CNSS catalogue for California for forecasting earthquakes with magnitudes above 5.75 and on the JMA catalogue of Japan for magnitudes above 6.75. Here we test whether the Ψ scaling relations extend to lower magnitudes, by applying EEPAS to depth-restricted subsets of the NIED catalogue of the Kanto area, central Japan, for magnitudes above 4.75. As in previous studies, the EEPAS model is found to be more informative than a quasi-static baseline model based on proximity to past earthquakes, and much more informative than the stationary uniform Poisson model. The information that it provides is illustrated by maps of the earthquake occurrence rate density, covering magnitudes from 5.0 to 8.0, for the central Japan region as at the beginning of year 2004, using the NIED and JMA catalogues to mid-2003.  相似文献   

6.
L.L. Romashkova   《Tectonophysics》2009,470(3-4):329-344
Following our experience in intermediate-term monitoring of seismicity before large earthquakes worldwide and in some regions we apply a similar approach to the analysis of seismicity on global scale several decades before the December 26, 2004 Sumatra–Andaman mega-earthquake. The Earth lithosphere is being considered as a single whole, representing the ultimate level of the complex Earth's hierarchy. The study brings up the following questions: Are there any anomalies of the global seismicity behaviour observed in advance of the mega-earthquake? Do these anomalies correspond to the global scale seismic patterns similar to those detected on the regional scale before great, major, and strong events? In other words, does the Earth lithosphere considered as a single whole show up an approach of the mega-earthquake in the way typical for events on the smaller levels of hierarchy? The results of our investigation favour the positive answer to the questions. Specifically we found that during the decade before the Sumatra–Andaman mega-earthquake the Earth lithosphere reveals, at least in intermediate-term scale, classical symptoms of instability, which can be depicted by known precursory seismic patterns. These are: (i) transformation of frequency–magnitude distribution, (ii) change in the rate of seismic activity, and (iii) depth redistribution of activity. Moreover changes of dependencies between magnitudes of different types are detected at the same time. The observed global scale patterns of collective behaviour of seismicity may indicate the state of criticality of the Earth lithosphere before the Sumatra–Andaman mega-earthquake.  相似文献   

7.
The Pattern Informatics (PI) technique [Tiampo, K.F., Rundle, J.B., McGinnis, S., Gross, S., Klein, W., 2002. Mean-field threshold systems and phase dynamics: An application to earthquake fault systems, Europhys. Lett., 60, 481–487] is founded on the premise that changes in the seismicity rate are a proxy for changes in the underlying stress. This new approach to the study of seismicity quantifies its local and regional space–time patterns and identifies regions of local quiescence or activation. Here we use a modification of the PI method to quantify localized changes surrounding the epicenters of large earthquakes in California in an attempt to objectively quantify the rupture zones of these upcoming events. We show that this method can be used to forecast the size and magnitude of future earthquakes.  相似文献   

8.
We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆmax, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an MS-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆmax,RP6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.  相似文献   

9.
A probabilistic estimate of seismic hazard can be obtained from the spatial distribution, of earthquake sources, their frequency–magnitude distribution and the rate of attenuation of strong ground motion with distance. We calculate the earthquake perceptibility, i.e. the annual probability that a particular level of ground shaking will be generated by earthquakes of particular magnitude, by weighting frequency–magnitude data with the predicted felt area for a given level of ground shaking at a particular magnitude. This provides an earthquake selection criterion that can be used in the anti-seismic design of non-critical structures. We calculate the perceptibility, at a particular value of isoseismal intensity, peak ground acceleration and velocity, as a function of source magnitude and frequency for the broad Aegean area using local attenuation laws. We use frequency–magnitude distributions that were previously obtained by combining short-term catalogue data with tectonic moment rate data for 14 tectonic zones in Greece with sufficient earthquake data, and where contemporary strain rates are available from satellite data. Many of the zones show a ‘characteristic earthquake’ distribution with the most perceptible earthquake equal to the maximum magnitude earthquake, but a relatively flat perceptibility between magnitudes 6 and 7. The maximum perceptible magnitude is in the fastest-deforming region in the middle of the Aegean sea, and tends to be systematically low on the west in comparison to the east of the Aegean sea. The tectonic data strongly constrain the long-term recurrence rates and lead to low error estimates (±0.2) in the most perceptible magnitudes.  相似文献   

10.
We test the Bowman and King [Bowman, D.D., King, G.C.P., 2001a, Accelerating seismicity and stress accumulation before large earthquakes. Geophys. Res. Lett., 28 (21), 4039–4042, Bowman, D.D., King, G.C.P., 2001b. Stress transfer and seismicity changes before large earthquakes. C. R. Acad. Sci. Paris, 333, 591–599] Stress Accumulation model by examining the evolution of seismicity rates prior to the 1992 Landers, California earthquake. The Stress Accumulation (SA) model was developed to explain observations of accelerating seismicity preceding large earthquakes. The model proposes that accelerating seismicity sequences result from the tectonic loading of large fault structures through aseismic slip in the elasto-plastic lower crust. This loading progressively increases the stress on smaller faults within a critical region around the main structure, thereby causing the observed acceleration of precursory activity. A secondary prediction of the SA model is that the precursory seismicity rates should increase first at the edges of the critical region, with the rates gradually rising over time at closer distances to the main fault. We test this prediction by examining year-long seismicity rates between 1960 and 2004, as a function of distance from the Landers rupture. To quantify the significance of trends in the seismicity rates, we auto-correlate the data, using a range of spatial and temporal lags. We find weak evidence for increased seismicity rates propagating towards the Landers rupture, but cannot conclusively distinguish these results from those obtained for a random earthquake catalog. However, we find a strong indication of periodicity in the rate fluctuations, as well as high correlation between activity 130–170 km from Landers and seismicity rates within 50 km of the Landers rupture temporally offset 1.5–2 years. The implications of this spatio–temporal correlation will be addressed in future studies.  相似文献   

11.
Both the genesis and rates of activity of shallow intraplate seismic activity in central Chile are poorly understood, mainly because of the lack of association of seismicity with recognizable fault features at the surface and a poor record of seismic activity. The goal of this work is to detail the characteristics of seismicity that takes place in the western flank of the Andes in central Chile. This region, located less than 100 km from Santiago, has been the site of earthquakes with magnitudes up to 6.9, including several 5+ magnitude shocks in recent years. Because most of the events lie outside the Central Chile Seismic Network, at distances up to 60 km to the east, it is essential to have adequate knowledge of the velocity structure in the Andean region to produce the highest possible quality of epicentral locations. For this, a N–S refraction line, using mining blasts of the Disputada de Las Condes open pit mine, has been acquired. These blasts were detected and recorded as far as 180 km south of the mine. Interpretation of the travel times indicates an upper crustal model consisting of three layers: 2.2-, 6.7-, and 6.1-km thick, overlying a half space; their associated P wave velocities are 4.75–5.0 (gradient), 5.8–6.0 (gradient), 6.2, and 6.6 km/s, respectively.Hypocentral relocation of earthquakes in 1986–2001, using the newly developed velocity model, reveals several regions of concentrated seismicity. One clearly delineates the fault zone and extensions of the strike-slip earthquake that took place in September 1987 at the source of the Cachapoal River. Other regions of activity are near the San José volcano, the source of the Maipo River, and two previously recognized lineaments that correspond to the southern extension of the Pocuro fault and Olivares River. A temporary array of seismographs, installed in the high Maipo River (1996) and San José volcano (1997) regions, established the hypocentral location of events with errors of less than 1 km. These events are clustered along no particular lineament approximately 25 km away from the San José and Maipo volcanoes. Recurrence intervals, based on a frequency magnitude relationship for lower magnitude events, indicate that earthquakes with magnitudes of 4.7 and 7 have a repeat time of 1 and 1200 years, respectively. Focal mechanisms of the two largest events indicate horizontal maximum and minimum compressive stresses with σ1 varying from a NW–SE orientation in the north to E–W at the southern extreme.  相似文献   

12.
In this study, we analyze the recent (1990–1997) seismicity that affected the northern sector (Sannio–Benevento area) of the Southern Apennines chain. We applied the Best Estimate Method (BEM), which collapses hypocentral clouds, to the events of low energy (Md max=4.1) seismic sequences in order to constrain the location and geometry of the seismogenetic structures. The results indicate that earthquakes aligned along three main structures: two sub-parallel structures striking NW–SE (1990–1992, Benevento sequence) and one structure striking NE–SW (1997, Sannio sequence). The southernmost NW–SE structure, which dips towards NE, overlies the fault that is likely to be responsible for a larger historical earthquake (Io max=XI MCS, 1688 earthquake). The northernmost NW–SE striking structure dips towards SW. The NE–SW striking structure is sub-vertical and it is located at the northern tip of the fault segment supposed to be responsible for the 1688 earthquake. The spatio-temporal evolution of the 1990–1997 seismicity indicates a progressive migration from SE (Benevento) to NW (Sannio) associated to a deepening of hypocenters (i.e., from about 5 to 12 km). Hypocenters cluster at the interface between the major structural discontinuities (e.g., pre-existing thrust surfaces) or within higher rigidity layers (e.g., the Apulia carbonates). Available focal mechanisms from earthquakes occurred on the recognized NW–SE and NE–SW faults are consistent with dip-slip normal solutions. This evidences the occurrence of coexisting NW–SE and NE–SW extensions in Southern Apennines.  相似文献   

13.
Observations indicate that earthquake faults occur in topologically complex, multi-scale networks driven by plate tectonic forces. We present realistic numerical simulations, involving data-mining, pattern recognition, theoretical analyses and ensemble forecasting techniques, to understand how the observable space–time earthquake patterns are related to the fundamentally inaccessible and unobservable dynamics. Numerical simulations can also help us to understand how the different scales involved in earthquake physics interact and influence the resulting dynamics. Our simulations indicate that elastic interactions (stress transfer) combined with the nonlinearity in the frictional failure threshold law lead to the self-organization of the statistical dynamics, producing 1) statistical distributions for magnitudes and frequencies of earthquakes that have characteristics similar to those possessed by the Gutenberg–Richter magnitude–frequency distributions observed in nature; and 2) clear examples of stress transfer among fault activity described by stress shadows, in which an earthquake on one group of faults reduces the Coulomb failure stress on other faults, thereby delaying activity on those faults. In this paper, we describe the current state of modeling and simulation efforts for Virtual California, a model for all the major active strike slip faults in California. Noting that the Working Group on California Earthquake Probabilities (WGCEP) uses statistical distributions to produce earthquake forecast probabilities, we demonstrate that Virtual California provides a powerful tool for testing the applicability and reliability of the WGCEP statistical methods. Furthermore, we show how the simulations can be used to develop statistical earthquake forecasting techniques that are complementary to the methods used by the WGCEP, but improve upon those methods in a number of important ways. In doing so, we distinguish between the “official” forecasts of the WGCEP, and the “research-quality” forecasts that we discuss here. Finally, we provide a brief discussion of future problems and issues related to the development of ensemble earthquake hazard estimation and forecasting techniques.  相似文献   

14.
The aim of the present work is to compile and update a catalogue of the instrumentally recorded earthquakes in Egypt, with uniform and homogeneous source parameters as required for the analysis of seismicity and seismic hazard assessment. This in turn requires a detailed analysis and comparison of the properties of different available sources, including the distribution of events with time, the magnitude completeness, and the scaling relations between different kinds of magnitude reported by different agencies. The observational data cover the time interval 1900–2004 and an area between 22°–33.5° N and 25°–36° E. The linear regressions between various magnitude types have been evaluated for different magnitude ranges. Using the best linear relationship determined for each available pair of magnitudes, as well as those identified between the magnitudes and the seismic moment, we convert the different magnitude types into moment magnitudes M W, through a multi-step conversion process. Analysis of the catalogue completeness, based on the M W thus estimated, allows us to identify two different time intervals with homogeneous properties. The first one (1900–1984) appears to be complete for M W ≥ 4.5, while the second one (1985–2004) can be considered complete for magnitudes M W ≥ 3.  相似文献   

15.
Southern Italy is dominated by extensional tectonics that in the Calabrian arc and Eastern Sicily produced the development of the Siculo–Calabrian Rift Zone (SCRZ). This zone is represented by a ≈ 370 km-long fault belt consisting of 10 to 50 km long distinct fault segments which extend both offshore and on land being also responsible of the crustal seismicity of this region. The geological and morphological observations indicate that the active normal faults of the SCRZ are characterized by throw-rates ranging from 0.7 to 3.1 mm/a. They accommodate an almost uniform horizontal extension-rate of about 3.0 mm/a along a WNW–ESE regional extension direction. Based on our field observations and following empirical relationships between magnitude and surface rupture length connections between large crustal earthquakes and distinct fault segments of the SCRZ have been also tentatively tested. Our data indicate moreover that the magnitudes (M) of the historical and instrumental earthquakes are consistent with the estimated values and that the geometry and kinematics of the fault segments and the related different crustal features of the SCRZ control the different seismic behaviours of adjacent portions of the active rift zone.  相似文献   

16.
Distributions of time between consecutive earthquakes verify an approximately universal scaling law for stationary seismicity. The shape of these distributions is shown to arise as a mixture of one distribution for short‐distance events and an exponential distribution for far‐off events, the distinction from short and long distances being relative to the size of the region studied. The distributions of consecutive distances show a double power law decay and verify an approximate scaling law which guarantees the simultaneous fulfillment of the scaling laws for time. The interplay between space and time can be seen as well by looking at the distribution of distances for a fixed time separation. These results suggest that seismicity can be understood as a series of intertwined independent continuous‐time random walks, with power law‐distributed waiting times and Lévy‐flight jumps. However, a simple model based on these ideas does not capture the invariance of seismicity under renormalization.  相似文献   

17.
Hot dry rock is becoming an important clean energy source. Enhanced geothermal systems (EGS) hold great promise for the potential to make a contribution to the energy inventory. However, one controversial issue associated with EGS is the impact of induced seismicity. In August 2019, a hydraulic stimulation experiment took place at the hot dry rock site of the Gonghe Basin in Qinghai, China. Earthquakes of different magnitudes of 2 or less occurred during the hydraulic stimulation. Correlations between hydraulic stimulation and seismic risk are still under discussion. Here, we analyze the hydraulic stimulation test and microseismic activity. We quantify the evolution of several parameters to explore the correlations between hydraulic stimulation and induced seismicity, including hydraulic parameters, microseismic events, b-value and statistical forecasting of event magnitudes. The results show that large-magnitude microseismic events have an upward trend with an increase of the total fluid volume. The variation of the b-value with time indicates that the stimulation experiment induces small amounts of seismicity. Forecasted magnitudes of events can guide operational decisions with respect to induced seismicity during hydraulic fracturing operations, thus providing the basis for risk assessment of hot dry rock exploitation.  相似文献   

18.
We present a revision and a seismotectonic interpretation of deep crust strike–slip earthquake sequences that occurred in 1990–1991 in the Southern Apennines (Potenza area). The revision is motivated by: i) the striking similarity to a seismic sequence that occurred in 2002  140 km NNW, in an analogous tectonic context (Molise area), suggesting a common seismotectonic environment of regional importance; ii) the close proximity of such deep strike–slip seismicity with shallow extensional seismicity (Apennine area); and iii) the lack of knowledge about the mechanical properties of the crust that might justify the observed crustal seismicity. A comparison between the revised 1990–1991 earthquakes and the 2002 earthquakes, as well as the integration of seismological data with a rheological analysis offer new constraints on the regional seismotectonic context of crustal seismicity in the Southern Apennines. The seismological revision consists of a relocation of the aftershock sequences based on newly constrained velocity models. New focal mechanisms of the aftershocks are computed and the active state of stress is constrained via the use of a stress inversion technique. The relationships among the observed seismicity, the crustal structure of the Southern Apennines, and the rheological layering are analysed along a crustal section crossing southern Italy, by computing geotherms and two-mechanism (brittle frictional vs. ductile plastic strength) rheological profiles. The 1990–1991 seismicity is concentrated in a well-defined depth range (mostly between 15 and 23 km depths). This depth range corresponds to the upper pat of the middle crust underlying the Apulian sedimentary cover, in the footwall of the easternmost Apennine thrust system. The 3D distribution of the aftershocks, the fault kinematics, and the stress inversion indicate the activation of a right-lateral strike–slip fault striking N100°E under a stress field characterized by a sub-horizontal N142°-trending σ1 and a sub-horizontal N232°-trending σ3, very similar to the known stress field of the Gargano seismic zone in the Apulian foreland. The apparent anomalous depths of the earthquakes (> 15 km) and the confinement within a relatively narrow depth range are explained by the crustal rheology, which consists of a strong brittle layer at mid crustal depths sandwiched between two plastic horizons. This articulated rheological stratification is typical of the central part of the Southern Apennine crust, where the Apulian crust is overthrusted by Apennine units. Both the Potenza 1990–1991 and the Molise 2002 seismic sequences can be interpreted to be due to crustal E–W fault zones within the Apulian crust inherited from previous tectonic phases and overthrusted by Apennine units during the Late Pliocene–Middle Pleistocene. The present strike–slip tectonic regime reactivated these fault zones and caused them to move with an uneven mechanical behaviour; brittle seismogenic faulting is confined to the strong brittle part of the middle crust. This strong brittle layer might also act as a stress guide able to laterally transmit the deviatoric stresses responsible for the strike–slip regime in the Apulian crust and may explain the close proximity (nearly overlapping) of the strike–slip and normal faulting regimes in the Southern Apennines. From a methodological point of view, it seems that rather simple two-mechanism rheological profiles, though affected by uncertainties, are still a useful tool for estimating the rheological properties and likely seismogenic behaviour of the crust.  相似文献   

19.
The study computes time-dependant earthquake probabilities on the basis of seismicity data mainly deriving from historic records. It provides a methodological approach useful for those countries where the scarcity of instrumental data and/or paleoseismological evidences requires that historical information shall be stressed. Thus, the conditional probability that damaging earthquakes (M ≥ 6) may occur in Italy in the next 30 years is shown, and the potential for the main worldwide known Italian cities with a cultural heritage is outlined. Earthquake probabilities are computed referring to the application of renewal processes, where the periodicity is analytically modelled by means of the Brownian Passage Time function; an estimate of the dispersion (i.e., uncertainty) introduced on probabilities is provided making use of Monte Carlo simulations. The computed probabilities refer to seismic source zones deriving from the spatial clustering of the historically documented seismicity. The computation of probabilities based on the interaction of earthquakes occurring in nearby zones, has been also attempted for a test area to explore the influence exerted by the stress transfer effect. The main findings of this study are that (1) seismic source zones in Southern Italy are the most prone to experience damaging earthquakes in the next 30-years, with conditional probabilities a large as 10%; and (2) the influence exerted by the earthquake interaction in increasing such probabilities, doesn’t seem to be relevant, because the mean recurrence times of large earthquakes (above the threshold magnitude of six chosen in this study) are in general much longer than the time shortening produced by the stress transfer.  相似文献   

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

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

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