首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 421 毫秒
1.
Multifractal of spatial distribution of seismicity in Liaoning area   总被引：1，自引：0，他引：1  
Making use of multifractal theory and corresponding computational method and according to the feature of evolution of spatial distribution with respect to seismicity by earthquake data in Liaoning area, earthquake activity of the area has been studied in detail. The results show that the evolution of increase in seismicity and distributive process in space are a multifractal structure. Whole characteristic of evolution in fractal increasing process of seismicity is described by obvious variation in regard toτ(q)-q curve,f(α) spectrum and other parameters before and after moderate and strong earthquakes.  相似文献   

2.
The geometrical and scaling properties of the energy flux of the turbulent kinetic energy in the solar wind have been studied. Using present experimental technology in solar wind measurements we cannot directly measure the real volumetric dissipation rate, <varepsilon>(t), but are constrained to represent it by its surrogate the energy flux near the dissipation range at the proton gyro scale. There is evidence for the multifractal nature of the so defined dissipation field <varepsilon>(t), a result derived from the scaling exponents of its statistical moments. The generalized dimension D q has been determined and reveals that the dissipation field has a multifractal structure, which is not compatible with a scale-invariant cascade. The related multifractal spectrum f(<alpha>) has been estimated for the first time for MHD turbulence in the solar wind. Its features resemble those obtained for turbulent fluids and other nonlinear multifractal systems. The generalized dimension D q can for turbulence in high-speed streams be fitted well by the functional dependence of the p-model with a comparatively large parameter p 1=0.87, indicating a strongly intermittent multifractal energy cascade. The experimental value for D p/3 used in the scaling exponent s(p) of the velocity structure function gives an exponent that can describe some of the observations. The scaling exponent <mu> of the autocorrelation function of <varepsilon>(t) has also been directly evaluated, being 0.37. Finally, the mean dissipation rate was determined, which could be used in solar wind heating models.  相似文献   

3.
The study of seismic anomalies, related both to the temporal trend of aftershock sequences and to the temporal series of mainshocks, is important for an understanding of the physical processes relating to the existence and the characteristics of seismic precursors. The purpose of this work is to highlight some methodological aspects related to the observation of possible anomalies in the temporal decay of an aftershock sequence. It is realized by means of several parameters. We focused our work on an analysis of the Papua New Guinea seismic sequence that occurred on November 16, 2000. The magnitude of the mainshock is M = 8.2. The observed temporal series of shocks per day can be considered as a sum of a deterministic contribution and a stochastic contribution. If the decay can be modeled as a nonstationary Poisson process where the intensity function is equal to n(t) = K(t + c)−p + K 1, the number of aftershocks in a small time interval Δt is the mean value n(t)Δt, with a standard deviation σ = √n(t)Δt. We observe that there are some variations in seismicity that can be considered as seismic anomalies before the occurrence of a large aftershock. The data, checked according to completeness criteria, come from the website of the USGS NEIC data bank (). The text was submitted by the authors in English.  相似文献   

4.
The purpose of this work is to highlight some methodological aspects related to the observation of possible anomalies in the temporal decay of aftershocks temporal series following a mainshock with magnitude M ≥ 7.0. In this paper we present the results for the Taiwan seismic sequence started on 20 September 1999 (M = 7.7) by tuning some seismic parameters that show considerable variations during the aftershock decay process. In here we also present the results obtained using a fractal approach for the seismic sequence. Earthquakes belong to a class of phenomena known as multifractals. In general it is important to define the fractal dimension D, but sometimes is not useful if we are describing a natural phenomenon; so it is necessary to define D 0 called box-counting dimension and D 2 called correlation dimension, usually D 0 ≥ D 2. In the elaborations of the fractal dimensions, for this sequence, we have obtained values lesser than 1, with a greater tendency of aftershocks to clusterize in time before a large aftershock. This is coherent with the possible existence of seismic anomalies, that could occur before the large aftershock. We also report the results obtained by using the delta/sigma method described firstly in [Caccamo et al., 2005] and later applied to different seismic sequence. The observed temporal series of the aftershocks per day can be considered as a sum of a deterministic and a stochastic contribution. If the decay can be modeled as a non-stationary Poissonian process, the number of aftershocks in a small time interval Δt is the mean value n(t) Δt, with a standard deviation (δ = √n(t)Δt. Investigating both aftershock behavior and a wide spectrum of parameters may find the key to explain better the mechanism of seismicity as a whole.  相似文献   

5.
It is shown that within the framework of the Kolmogorov model the “energy” of the pole E(t) = x 12 + x 22 can be interpreted as a Markovian process. The exact analytical expression has been obtained for the density of the conditional probability of the quantity E(t) and the problem of the first passage time of the process E(t) has been analyzed. It was shown that the available data on the swing of the function E(t) are not at variance with the Kolmogorov model and a short-period drop of the amplitude of the Chandler wobble in the early 20th century fits this model at Q = 50–200 too; values of Q > 350 are less reasonable.  相似文献   

6.
It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).  相似文献   

7.
Abstract

Bayly (1993) introduced and investigated the equation (? t + v·▽-η ▽2)S = RS as a scalar analogue of the magnetic induction equation. Here, S(r, t) is a scalar function and the flow field v(r, t) and “stretching” function R(r, t) are given independently. This equation is much easier to handle than the corresponding vector equation and, although not of much relevance to the (vector) kinematic dynamo problem, it helps to study some features of the fast dynamo problem. In this note the scalar equation is considered for linear flow and a harmonic potential as stretching function. The steady equation separates into one-dimensional equations, which can be completely solved and therefore allow one to monitor the behaviour of the spectrum in the limit of vanishing diffusivity. For more general homogeneous flows a scaling argument is given which ensures fast dynamo action for certain powers of the harmonic potential. Our results stress the singular behaviour of eigenfunctions in the limit of vanishing diffusivity and the importance of stagnation points in the flow for fast dynamo action.  相似文献   

8.
Chronology and geochemistry of the Shangyu gabbro-diorite in western Shandong were studied to understand their petrogenesis and the nature of the Mesozoic lithospheric mantle. The Shangyu intru-sion is mainly composed of a suite of gabbro-diorite. Zircons from the intrusion display eu-hedral-subhedral in shape and have high Th/U ratios (1.23―2.87), implying their magmatic origin. LA-ICP-MS zircon U-Pb dating results for two samples indicate that they were formed in the Early Cre-taceous, yielding weighted mean 206Pb/238U ages of 129±1Ma and 134±2Ma, respectively. Except for early cumulate such as sample QT-19, their SiO2 and MgO contents range from 50.12% to 56.37% and from 3.52% to 6.37%, respectively. Moreover, the gabbro-diorites are characterized by high Mg# (0.54―0.63), enrichment in Na (Na2O/K2O ratios more than 1), Cr (73×10-6―217×10-6) and Ni (34×10-6―241×10-6), and intensive enrichments in light rare earth elements (LREEs) and large ion lithophile elements (LILEs) and depletion in high field strength elements (HFSEs). Their initial 87Sr/86Sr ratios and ε Nd(t) values range from 0.70962 to 0.71081 and from-16.60 to-13.04, respectively. Taken together with the Early Creta-ceous high-Mg diorites and the mantle xenoliths from the Tietonggou and Jinling as well as basalts from the Fangcheng and Feixian, it is suggested that the primary magma for the Shangyu gab-bro-diorites should be derived from the enriched lithospheric mantle intensively modified by conti-nental crust. The Sr-Nd-Pb isotopic compositions for the Early Cretaceous high-Mg diorites in western Shandong display a trend of spatial variations, i.e., initial 87Sr/86Sr, 207Pb/204Pb and 208Pb/204Pb ratios de-creasing and ε Nd(t) values increasing from southeast to northwest in western Shandong, which is con-sistent with the tectonic model that the Yangtze Craton subducted beneath the North China Craton oriented in north-west direction in the Early Mesozoic.  相似文献   

9.

The multifractal fields modeling the spatial distribution of the epicenters and temporal behavior of weak seismicity before the earthquake that occurred on June 17, 2000 in Iceland (M = 6.6) are analyzed. It is shown that the main shock of this event was preceded by the broadening of the f(a) spectra of the proxy fields of both the spatial distribution of the epicenters and the temporal behavior of weak seismicity. The physical interpretation of this effect is presented.

  相似文献   

10.
Seismicity parameters for important urban agglomerations in India   总被引：1，自引：0，他引：1  
India’s urban population has increased in the recent times. An earthquake near an urban agglomeration has the potential to cause severe damage. In this article, seismicity parameters for region surrounding important urban agglomerations in India are estimated. A comprehensive earthquake catalogue for the region (6°E–42°E latitude and 60°N–100°N longitude) including historic and pre-historic events has been compiled from various sources. To estimate the parameters, past earthquake data in a control region of radius 300 km has been assembled to quantify the seismicity around each urban agglomeration. The collected earthquake data is first evaluated for its completeness. From combined (historical and instrumental) data, the seismicity parameters b-value, seismic activity rate, λ and maximum expected magnitude (m max ) have been obtained from the methodology proposed by Kijko and Graham (1998). The obtained activity rates indicate that region surrounding Guwahati urban agglomeration is the most seismically active region followed by Srinagar, Patna, Amritsar and Chandigarh.  相似文献   

11.
ＩｎｔｒｏｄｕｃｔｉｏｎＳｉｎｃｅｔｈｅ１９６０′ｓ,ｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｍｏｄｅｒｎｓｃｉｅｎｃｅａｎｄｔｅｃｈｎｏｌｏｇｙｈａｓｇｒａｄｕａｌｙｍａｄｅｉｔｐｏｓｉｂｌｅｔｏｐｒｅｄｉｃｔｅａｒｔｈｑｕａｋｅｓａｎｄｈａｓｐｒｏｍｏｔ...  相似文献   

12.
水丰水库亚洲公鱼和大银鱼种群生态的比较研究   总被引：4，自引：0，他引：4  
1998年 10月至 2 0 0 0年 11月对水丰水库亚洲公鱼和大银鱼的食性、生长、繁殖、群体结构进行了比较研究 .亚洲公鱼为杂食性鱼类但主要以浮游动物为食 ,繁殖时还摄食同类卵 ,出现频率达 4 2 .9% ,在冬季也摄食大银鱼卵 (2 .8% ) ;体长 (L)、体重 (W)呈异速生长 (W =1.176×10 - 6 L3.5555) ;群体繁殖速度比 1984年以前提高了 4 9.2 % ;群体结构中体长、体重的离散度很小 .大银鱼前期主要以浮游动物为食 ,后期以小鱼小虾为主要食物 ,也捕食公鱼 ,出现频率 6 .9% -9 6 % ;体长、体重呈等速生长 ,W =4 .4 82× 10 - 6 L3.0 0 0 9;群体繁殖速度高出公鱼 1.18倍 ;群体结构中体长、体重的离散度较大 ,表现出了个体间的差异性  相似文献   

13.
The source parameters of the M W = 7.6 Olyutorskii earthquake were estimated using the moments of the slip rate function with degrees 1 and 2. The moments were estimated from broadband P-wave records at 52 stations of the worldwide network. The first step was to find a function S(t) for each station; this function is an apparent source time function, i.e., the P-wave slip as radiated by the source toward a station under consideration. The method of empirical Green’s functions was used to estimate S(t). The next step was to calculate the moments of S(t) of degrees 1 and 2 over time and to set up relevant equations to be solved by least squares for the unknown source moments. The horizontal linear source was used as a nonparametric model for calculating the source moments. Haskell’s parametric model was used for further interpretation of the source moments. The resulting estimates are as follows: the source centroid was 13–25 km southwest of the epicenter, the source was 105–120 km long, the source strike was 222°–228°, the rupture velocity was 2.7–3.0 km/s, and the total radiation duration was 24–27 s. These estimates indicate a bilateral rupture dominated by a southwestward sense of rupture propagation. The source characteristics are consistent with the aftershock area geometry and with the focal mechanism, as well as with surface breakage as observed by geologists in the field.  相似文献   

14.
Ａｌｉｇｎｍｅｎｔｓｉｌｋｗｏｒｍｓａｓｓｅｉｓｍｉｃａｎｉｍａｌａｎｏｍａｌｏｕｓｂｅｈａｖｉｏｒ（ＳＡＡＢ）ａｎｄｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｏｄｅｌｏｆａｆａｕｌｔ：ａｔｈｅｏｒｙａｎｄｌａｂｏｒａｔｏｒｙｅｘｐｅｒｉｍｅｎｔＭＯＴＯ...  相似文献   

15.
In this paper we suggest that conditional estimator/predictor of rockburst probability (and rockburst hazard, P T (t)) can be approximated with the formula P T (t) = P 1(θ 1)…P N (θ N )·P dyn T (t), where P dyn T (t) is a time-dependent probability of rockburst given only the predicted seismic energy parameters, while P i (θ i ) are amplifying coefficients due to local geologic and mining conditions, as defined by the Expert Method of (rockburst) Hazard Evaluation (MRG) known in the Polish mining industry. All the elements of the formula are (approximately) calculable (on-line) and the resulting P T value satisfies inequalities 0 ≤ P T (t) ≤ 1. As a result, the hazard space (0–1) can be always divided into smaller subspaces (e.g., 0–10−5, 10−5–10−4, 10−4–10−3, 10−3–1), possibly named with symbols (e.g., A, B, C, D, …) called “hazard states” — which saves the prediction users from worrying of probabilities. The estimator P T can be interpreted as a formal statement of (reformulated) Comprehensive Method of Rockburst State of Hazard Evaluation, well known in Polish mining industry. The estimator P T is natural, logically consistent and physically interpretable. Due to full formalization, it can be easily generalized, incorporating relevant information from other sources/methods.  相似文献   

16.
Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains εxx (r, t), εyy (r, t) and εzz (r, t) and the bulk strain θ (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (α, β, γ) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.  相似文献   

17.
Tonalites from the island arc rock assemblage in the Zêtang segment of the Yarlung Zangbo suture zone were analyzed for major, trace elements (including REE) and Sr-Nd isotope. The experimental datademonstrate that the tonalites have the adakite-like characteristics, including high SiO2 (58%-63%),Al2O3 (18.4%-22.4%), Sr (810×10-6-940×10-6), Sr/Y (77-106), low HREE (Y=9×10-6-11×10-6, Yb=1×10-6-1.3×10-6), with LREE enrichment and faint Eu positive anomaly. Isr (0.70421-0.70487) is relatively low whereas 143Nd/144Nd (0.512896-0.512929) and εNd(t) values ( 6.7- 7.3) are high. These feainvolvement of a small amount of oceanic sediments. The identification of Z(e)tang adakites, derived from slab melting, presents new evidence for the intra-Tethyan subduction and the previous suggestion about the existence of intra-oceanic island arc within Tethys.  相似文献   

18.
Two empirical magnitude scaling relationships, predominant period (tpmax{\tau _p^{\max}}) and peak ground displacement (Pd) magnitudes, were investigated for the first 3 s after P-wave arrivals using 1,412 vertical waveforms recorded by the Korea National Seismic Network (KNSN) between 2001 and 2007. To evaluate the accuracy of the derived magnitude relationships, we simulated off-line ElarmS tests using 65 events occurring inside the KNSN. While the average magnitude error was ∼0.70 magnitude units when using only the closest station to the epicentre, the error dropped to ∼0.62 and ∼0.42 magnitude units when using the closest two and closest four stations, respectively. For events M L ≥ 3.0, the average magnitude error was ∼0.33 and showed stable values when the closest four stations were available. Our magnitude scaling relationships may be useful for initial work in developing an earthquake early warning system in South Korea.  相似文献   

19.
Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S i i (t)>, and the plume centroid covariance, R i i (t) (i=1,2), were simulated for the variance of Y=log T, σ Y 2=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ Y 2=0.1, all simulated <S i i (t)>−S i i (0) and R i i (t) agree well with the first-order theoretical values, where S i i (0) are the initial values of S i i (t). For σ Y 2=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S 11(t)>−S 11(0) and R 11(t), agree well with the first-order theoretical results but the simulated transverse moments <S 22(t)>−S 22(0) and R 22(t) are significantly larger than the first-order values. For the same two larger values of σ Y 2 but the line sources parallel to μ, the simulated <S 11(t)>−S 11(0) are larger than but the simulated R 11 are smaller than the first-order values, and both simulated <S 22(t)>−S 22(0) and R 22(t) stay larger than the first-order values. For a fixed value of σ Y 2, the summations of <S i i (t)>−S i i (0) and R i i , i.e., X i i (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X 11 are in good agreement with the first-order theory but the simulated X 22 are significantly larger than the first-order values. The simulated X 22, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S 22(t)>−S 22(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ Y 2 is small.  相似文献   

20.
 Simulated gas-driven eruptions using CO2–water-polymer systems are reported. Eruptions are initiated by rapidly decompressing CO2–saturated water containing up to 1.0 wt.% CO2. Both cylindrical test cells and a flask test cell were used to examine the effect of magma chamber/conduit geometry on eruption dynamics. Bubble-growth kinetics are examined quantitatively in experiments using cylindrical test cells. Uninhibited bubble growth can be roughly expressed as dr/dt≈λD(β-1)/(γt 1/3) for a CO2–water-polymer system at 0–22  °C and with viscosities up to 5 Pa·s, where r is the radius of bubbles, λ and D are the Ostwald solubility coefficient and diffusivity of the gas in the liquid, β is the degree of saturation (decompression ratio), and γ characterizes how the boundary layer thickness increases with time and is roughly 1.0×10–5 m/s1/3 in this system. Unlike the radius of cylindrical test cells, which does not affect the eruption threshold and dynamics, the shape of the test cells (flask vs cylindrical) affects the dynamics but not the threshold of eruptions. For cylindrical test cells, the front motion is characterized by constant acceleration with both Δh (the height increase) and ΔV (the volume increase) being proportional to t 2; for the flask test cell, however, neither Δh nor ΔV is proportional to t 2 as the conduit radius varies. Test-cell geometry also affects foam stability. In the flask test cell, as it moves from the wider base chamber into the narrower conduit, the bubbly flow becomes fragmented, affecting the eruption dynamics. The fragmentation may be caused by a sudden increase in acceleration induced by conduit-shape change, or by the presence of obstacles to the bubbly flow. This result may help explain the range in vesicularities of pumice and reticulite. Received: 16 May 1997 / Accepted: 11 October 1997  相似文献   

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

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