首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
In this paper, the micro-cracks in the brittle rocks are assumed to be penny shaped and evenly distributed; the damage and dilatancy of the brittle rocks is attributed to the growth and expansion of numerous micro-cracks under the local tensile stress. A single crack’s behaviour under the local tensile stress is generalized to all cracks based on the distributed damage mechanics. The relationship between the local tensile stress and the external loading is derived based on the Maxwell model. The damage factor corresponding to the external loading is represented using the p–alpha (pα) model. A dilatancy equation that can build up a link between the external loading and the rock dilatancy is established. A test of dilatancy of a brittle rock under triaxial compression is conducted; the comparison between experimental results and our theoretical results shows good consistency.  相似文献   

2.
Recent estimates of fracture energy G in earthquakes show a power-law dependence with slip u which can be summarized as G u a where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to G f : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G and G f should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s2). The resulting fracture energy G f is similar to the seismological estimates, with G f and G being comparable over most of the slip range. However, G f appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G appears to increase further and surpasses G f at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.  相似文献   

3.
Rock brittleness directly affects reservoir fracturing and its evaluation is essential for establishing fracturing conditions prior to reservoir reforming. Dynamic and static brittleness data were collected from siltstones of the Qingshankou Formation in Songliao Basin. The brittle–plastic transition was investigated based on the stress–strain relation. The results suggest that the brittleness indices calculated by static elastic parameters are negatively correlated with the stress drop coefficient and the brittleness index B2, defined as the average of the normalized Young’s modulus and Poisson’s ratio, is strongly correlated with the stress drop. The brittleness index B2, Young’s modulus, and Poisson’s ratio correlate with the brittle minerals content; that is, quartz, carbonates, and pyrite. We also investigated the correlation between pore fluid and porosity and dynamic brittle characteristic based on index B2. Pore fluid increases the plasticity of rock and reduces brittleness; moreover, with increasing porosity, rock brittleness decreases. The gas-saturated siltstone brittleness index is higher than that in oil- or water-saturated siltstone; the difference in the brittleness indices of oil- and water-saturated siltstone is very small. By comparing the rock mechanics and ultrasonic experiments, we find that the brittleness index obtained from the rock mechanics experiments is smaller than that obtained from the ultrasonic experiments; nevertheless, both decrease with increasing porosity as well as their differences. Ultrasonic waves propagate through the rock specimens without affecting them, whereas rock mechanics experiments are destructive and induce microcracking and porosity increases; consequently, the brittleness of low-porosity rocks is affected by the formation of internal microcrack systems.  相似文献   

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

5.
The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.  相似文献   

6.
We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h Am , which corresponds to the maximum increase, and layer peak height h m F, while h Am is always less than h m F. The difference between h Am and h m F linearly increases with increasing h m F. Whereas the difference is ~38 km for h m F = 280 km, it is ~54 km for h m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.  相似文献   

7.
Based on the known forecast of solar cycle 25 amplitude (Rz max ≈ 50), the first assessments of the shape and amplitude of this cycle in the index of solar activity F10.7 (the magnitude of solar radio flux at the 10.7 cm wavelength) are given. It has been found that (F10.7)max ≈ 115, which means that it is the lowest solar cycle ever encountered in the history of regular ionospheric measurements. For this reason, many ionospheric parameters for cycle 25, including the F2-layer peak height and critical frequency (hmF2 and foF2), will be extremely low. For example, at middle latitudes, typical foF2 values will not exceed 8–10 MHz, which makes ionospheric heating ineffective in the area of upper hybrid resonance at frequencies higher than 10 MHz. The density of the atmosphere will also be extremely low, which significantly extends the lifetime of low-orbit satellites. The probability of F-spread will be increased, especially during night hours.  相似文献   

8.
The rupture processes of any heterogeneous material constitute a complex physical problem. Earthquake aftershocks show temporal and spatial behaviors which are consequence of the heterogeneous stress distribution and multiple rupturing following the main shock. This process is difficult to model deterministically due to the number of parameters and physical conditions, which are largely unknown. In order to shed light on the minimum requirements for the generation of aftershock clusters, in this study, we perform a simulation of the main features of such a complex process by means of a fiber bundle (FB) type model. The FB model has been widely used to analyze the fracture process in heterogeneous materials. It is a simple but powerful tool that allows modeling the main characteristics of a medium such as the brittle shallow crust of the earth. In this work, we incorporate spatial properties, such as the Coulomb stress change pattern, which help simulate observed characteristics of aftershock sequences. In particular, we introduce a parameter (P) that controls the probability of spatial distribution of initial loads. Also, we use a “conservation” parameter (π), which accounts for the load dissipation of the system, and demonstrate its influence on the simulated spatio-temporal patterns. Based on numerical results, we find that P has to be in the range 0.06 < P < 0.30, whilst π needs to be limited by a very narrow range (0.60 < π < 0.66) in order to reproduce aftershocks pattern characteristics which resemble those of observed sequences. This means that the system requires a small difference in the spatial distribution of initial stress, and a very particular fraction of load transfer in order to generate realistic aftershocks.  相似文献   

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

10.
The consideration of the relation between the daytime and nighttime values of the critical frequency F2, foF2 of the ionospheric F2 layer, started in the previous publication of the authors, is continued. The main regularities in variations in the correlation coefficient R(foF2) characterizing this relation are confirmed using larger statistical material (more ionospheric stations and longer observational series). Long-term trends in the R(foF2) value are found: at all stations the negative value of R(foF2) increases with time after 1980.  相似文献   

11.
High-pressure (HP) and ultrahigh-pressure (UHP) eclogites exposed in collisional orogens are widely regarded to record the history of crustal rocks that were subducted to mantle depths and exhumed back to the surface. Insight into subduction and exhumation processes plays an important role in understanding the nature and evolution of subduction zones, geodynamics and plate tectonics. In contrast to continental-type HP to UHP metamorphic rocks that are dominated by felsic lithology, oceanic-type HP to UHP metamorphic rocks are dominated by mafic eclogites and thus have greater density, and their exhumation needs to overcome large barriers and may involve complicated tectonic processes. The exhumation of HP to UHP rocks is mainly influenced by the internal buoyancy, however, the external tectonic forces (such as channel flow) also act as effective exhumation drivers; in addition, effects of tectonic settings (such as slab rollback and breakoff) should take into account. The HP-UHP metamorphic terrane in Southwestern Tianshan, which mainly comprises of metasediments with interlayered metamafic lenses and blocks, represents a typical accretionary mélange associated with deep subduction of oceanic crust. However, the exhumation mechanisms of these once deeply buried HP-UHP rocks are still under discussion. Based on the field occurrences, petrographic features, peak metamorphic P-T conditions and ages of the eclogites/blueschists and their metasedimentary country rocks, a “sediment-type subduction channel” model is advocated in this study to appraise/decipher the formation and evolution of Southwestern Tianshan HP-UHP metamorphic belt. Poly-cyclic metamorphic eclogites record the multistage burial-exhumation cycling manipulated by convective flow in a channel-like interface between the plates, giving robust evidence for the presence of a subduction channel. In addition, this study summarizes some remaining geotectonic problems and research perspectives concerning the Southwestern Tianshan HP-UHP metamorphic belt.  相似文献   

12.
This paper investigated the effects of basin geometry and material property on the response of 2D trapezoidal sediment-filled basin to incident plane SH waves. Ten basin configurations with different geometries were developed, and then their seismic responses to both Ricker wavelets and seismic records were simulated by using an explicit finite difference scheme. The definition of deep/shallow basin, the precondition for the observation of prominent surface waves and the influential area of edge effects of the shallow basin were discussed quantitatively in this study. The followings were concluded: in the common velocity contrast range (v s1/v s2 < 10), the fundamental frequency a basin with W/H > 3.0 can be estimated approximately by 1D theory. The complexity of peak ground acceleration distribution pattern, the width of the most affected section as well as the amplitude of ground motion in the Edge Region increase with incident frequency. Prominent surface waves can only be observed when the incident wavelength is shorter than the critical wavelength λ c . The interaction between incident wave and basin dynamic property plays a dominant role on the peak ground acceleration amplitude while the interaction between incident wave and geometry plays a more significant role on the peak ground acceleration distribution. For very shallow basin, different areas along the basin width are affected to different extents. Only a limited area close to the basin edge is influenced significantly. It is more feasible to propose spectral aggravation factor for different surface zones respectively than a uniform constant as a tool to calibrate the 1D-based design spectrum so as to take the basin effects into account.  相似文献   

13.
Month-to-month changes in the statistical characteristics of the ionospheric E layer peak electron density NmE at medium and low geomagnetic latitudes under daytime geomagnetically quiet conditions are investigated. Critical frequencies of the ionospheric E layer measured by the middle latitude ionosonde Boulder and low latitude ionosondes Huancayo and Jicamarca at low solar activity from 1957 to 2015 have been used in the conducted statistical analysis. The mathematical expectation of NmE, standard deviation of NmE from the expectation of NmE, and NmE variation coefficient have been calculated for each month of the year. The months of the formation of extrema of these statistical parameters of NmE were found.  相似文献   

14.
An important task in seismic hazard assessment is estimation of the intensity and frequency of extremely strong earthquake effects, in particular, peak ground velocities (PGV). Earlier, a method was proposed to evaluate PGV values based on the magnitude of displacements of rock blocks (Rodkin et al., 2012). In this study, this method is used to analyze field data on the source zones of the August 19, 1992, MS = 7.3 Susamyr earthquake and the January 3, 1911, Mw = 7.9 Kemin earthquake, and estimate maximum ground shaking at the upper construction site of the Upper Naryn series of hydropower plants, Kyrgyz Republic. It is shown that the resulting estimates are consistent with data obtained through other techniques. Therefore, the new approach can be recommended to estimate earthquake effects.  相似文献   

15.
The earthquake inter-event time distribution is studied, using catalogs for different recent aftershock sequences. For aftershock sequences following the Modified Omori’s Formula (MOF) it seems clear that the inter-event distribution is a power law. The parameters of this law are defined and they prove to be higher than the calculated value (2–1/p). Based on the analysis of the catalogs, it is determined that the probability densities of the inter-event time distribution collapse into a single master curve when the data is rescaled with instantaneous intensity, R(t; M th ), defined by MOF. The curve is approximated by a gamma distribution. The collapse of the data provides a clear view of aftershock-occurrence self-similarity.  相似文献   

16.
A search for trends k(foE) in the critical frequency of the ionospheric E layer at Juliusruh and Slough stations is performed by the method often used by the authors to analyze trends in the F2-layer parameters. It is found that k(foE) could differ in both magnitude and even sign within different time intervals. However, the k(foE) trends have been stably negative over the last two decades for both stations and all months of the year. The k(foE) values averaged over a year are ?0.012 and ?0.005 MHz per year for Juliusruh and Slough stations, respectively. The method used in the recent paper by La?tovi?ka et al. (2016) to determine foE trends is analyzed, and it is shown that the difference in linear approximation of the dependence of the observed foE values on F10.7 within different time intervals could be interpreted not as the presence of a different foE dependence on the F10.7 index within these intervals but as the presence within them of foE trends that change the slope of the linear approximation.  相似文献   

17.
We performed an analysis of mean daily variations in the ΔEz atmospheric electric field at the Hornsund (located near the polar cap boundary) and Kakioka (located at near-equatorial latitudes) observatories under magnetically quiet and weakly disturbed conditions. At both observatories, the mean daily variations in ΔEz were found to be mainly controlled by the location of the observation point with respect to the focuses of the convective vortices of the DP 0 system. The substorm evolution in the nightside of the magnetosphere (a sharp burst in the AE index) was shown to lead to negative variations in ΔEz in the dayside sector at polar latitudes (the Hornsund observatory) and positive deviations in ΔEz at premidnight time at equatorial latitudes (the Kakioka observatory). It is concluded that variations in ΔEz at the Kakioka observatory are largely controlled by the equatorial electrojet, which is maximal during day-time hours, and at the Hornsund observatory these variations are controlled by the auroral electrojet, which is maximal at night-time and early morning hours of local time.  相似文献   

18.
ADCP and temperature chain measurements have been used to estimate the rate of energy input by wind stress to the water surface in the south basin of Windermere. The energy input from the atmosphere was found to increase markedly as the lake stratified in spring. The efficiency of energy transfer (Eff), defined as the ratio of the rate of working in near-surface waters (RW) to that above the lake surface (P 10), increased from ~0.0013 in vertically homogenous conditions to ~0.0064 in the first 40 days of the stratified regime. A maximum value of Eff~0.01 was observed when, with increasing stratification, the first mode internal seiche period decreased to match the diurnal wind period of 24 h. The increase in energy input, following the onset of stratification was reflected in enhancement of the mean depth-varying kinetic energy without a corresponding increase in wind forcing. Parallel estimates of energy dissipation in the bottom boundary layer, based on determination of the structure function show that it accounts for ~15% of RW in stratified conditions. The evolution of stratification in the lake conforms to a heating stirring model which indicates that mixing accounts for ~21% of RW. Taken together, these estimates of key energetic parameters point the way to the development of full energy budgets for lakes and shallow seas.  相似文献   

19.
An experiment on evapotranspiration from citrus trees under irrigation with saline water was carried out for 4 months. Two lysimeters planted with a citrus tree in the green house were used. One lysimeter was irrigated with saline water (NaCl and CaCl2 of 2000 mg/L equivalence,EC = 3.8 dS/m, SAR = 5.9) and the other was irrigated with freshwater using drip irrigation. The applied irrigation water was 1.2 times that of the evapotranspiration on the previous day. Evapotranspiration was calculated as the change in lysimeter weight recorded every 30 minutes. The lysimeters were filled with soil with 95.8% sand. The results of the experiment were as follows. (i) The evapotranspiration from citrus tree was reduced after irrigation with saline water. The evapotranspiration returns to normal after leaching. However it takes months to exhaust the salt from the tree. (ii) To estimate the impact of irrigation with saline water on the evapotranspiration from citrus trees, the reduction coefficient due to salt stress (Ks) was used in this experiment. Evapotranspiration under irrigation with saline water (ET s ) can be calculated from evapotranspiration under irrigation with freshwater (ET) by the equationET s =K s × ET. Ks can be expressed as a function ofEC sw . (iii) The critical soil-water electrical conductivity (EC sw ) is 9.5 dS/m, beyond which adverse effects on evapotranspiration begin to appear. IfEC sw can be controlled at below 9.5 dS/m, saline water can be safely used for irrigation.  相似文献   

20.
We study the frictional and viscous effects on earthquake nucleation, especially for the nucleation phase, based on a one-degree-of-freedom spring-slider model with friction and viscosity. The frictional and viscous effects are specified by the characteristic displacement, U c, and viscosity coefficient, η, respectively. Simulation results show that friction and viscosity can both lengthen the natural period of the system and viscosity increases the duration time of motion of the slider. Higher viscosity causes a smaller amplitude of lower velocity motion than lower viscosity. A change of either U c (under large η) or η (under large U c) from a large value (U ch for U c and η h for η) to a small one (U cl for U c and η l for η) in two stages during sliding can result in a clear nucleation phase prior to the P-wave. The differences δU c = U ch ? U cl and δη = η h ? η l are two important factors in producing a nucleation phase. The difference between the nucleation phase and the P-wave increases with either δU c or δη. Like seismic observations, the peak amplitude of P-wave, which is associated with the earthquake magnitude, is independent upon the duration time of nucleation phase. A mechanism specified with a change of either η or U c from a larger value to a smaller one due to temporal variations in pore fluid pressure and temperature in the fault zone based on radiation efficiency is proposed to explain the simulation results and observations.  相似文献   

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

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