首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
The popularly used viscoelastic models have some shortcomings in describing relationship between quality factor (Q) and frequency, which is not consistent with the observation data. Based on the theory of viscoelasticity, a new approach to construct constant-Q viscoelastic model in given frequency band with three parameters is developed. The designed model describes the frequency-independence feature of quality factor very well, and the effect of viscoelasticity on seismic wave field can be studied relatively accurate in theory with this model. Furthermore, the number of required parameters in this model has been reduced fewer than that of other constant-Q models, this can simplify the solution of the viscoelastic problems to some extent. At last, the accuracy and application range have been analyzed through numerical tests. The effect of viscoelasticity on wave propagation has been briefly illustrated through the change of frequency spectra and waveform in several different viscoelastic models.  相似文献   

2.
The boundary value problems most frequently encountered in studies of elastic wave propagation in stratified media can be formulated in terms of a finite number of linear, first order and ordinary differential equations with variable coefficients. Volterra (1887) has shown that solutions to such a system of equations are conveniently represented by the product integral, or propagator, of the matrix of coefficients. In this paper we summarize some of the better known properties of propagators plus numerica methods for their computation. When the dispersion relation is somem th order minor of the integral matrix it is possible to deal withm th minor propagators so that the dispersion relation is a single element of them th minor integral matrix. In this way one of the major sources of loss of numerical accuracy in computing the dispersion relation is avoided. Propagator equations forSH and forP-SV waves are given for both isotropic and transversely isotropic media. In addition, the second minor propagator equations forP-SV waves are given. Matrix polynomial approximations to the propagators, obtained from the method of mean coefficients by the Cayley-Hamilton theorem and the Lagrange-Sylvester, interpolation formula, are derived.  相似文献   

3.
Accurate Q parameter is hard to be obtained, but there is great difference between Q measurements from different measurement methods in seismic physical modelling. The influence factors, stability and accuracy of different methods are analyzed through standard sample experiment and the seismic physical modelling. Based on this, we proposed an improved method for improving accuracy of pulse transmission method, in which the samples with similar acoustic properties to the test sample are selected as the reference samples. We assess the stability and accuracy of the pulse transmission, pulse transmission insertion, and reflection wave methods for obtaining the quality factor Q using standard and reference samples and seismic physical modeling. The results suggest that the Q-values obtained by the pulse transmission method are strongly affected by diffraction and the error is 50% or greater, whereas the relative error of the improved pulse transmission method is about 10%. By using a theoretical diffraction correction method and the improved measurement method, the differences among the Q-measuring methods can be limited to within 10%.  相似文献   

4.
A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q m started to appear from Nov. 9, 1999. The events giving the deduced high Q m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q m versus distance observation obviously deviates from the normal trend of Q m linearly increasing with distance. (2) The average Q m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q m after the M=5.1 event is 204.  相似文献   

5.
Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M L ?≤?4.4. Up to five seismograph stations were utilised to determine Q c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q 0 and α were determined for the attenuation relation Q c(f)?=?Q 0 f α . The result was Q 0?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t s???t 0) (time from origin time t 0 to the start of coda analysis window is 1.9 times the S-travel time, t s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q c was not calculated. The derived Q c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q S was Q 0?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q c result.  相似文献   

6.
Shear and compressional wave velocities, coupled with other petrophysical data, are very important for hydrocarbon reservoir characterization. In situ shear wave velocity (Vs) is measured by some sonic logging tools. Shear velocity coupled with compressional velocity is vitally important in determining geomechanical parameters, identifying the lithology, mud weight design, hydraulic fracturing, geophysical studies such as VSP, etc. In this paper, a correlation between compressional and shear wave velocity is obtained for Gachsaran formation in Maroon oil field. Real data were used to examine the accuracy of the prediction equation. Moreover, the genetic algorithm was used to obtain the optimal value for constants of the suggested equation. Furthermore, artificial neural network was used to inspect the reliability of this method. These investigations verify the notion that the suggested equation could be considered as an efficient, fast, and cost-effective method for predicting Vs from Vp.  相似文献   

7.
The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer. We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation. In this study, we drilled number of receiver holes around the source hole, each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys. We also propose a novel tomographic inversion of the Q factor without the effect of the source signature, and examine its stability and reliability using synthetic data. We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield. The results show that seismic absorption in the near-surface layer is much greater than that in the subsurface strata. Thus, it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption. In addition, we derive different Q factors from two frequency bands, which can be treated, to some extent, as evidence of a frequency-dependent Q.  相似文献   

8.
Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R 1.0 as the geometrical spreading function in the NMSZ, the Q P and Q S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q P?=?(115.80?±?1.36) f (0.495?±?0.129) and Q S?=?(161.34?±?1.73) f (0.613?±?0.067), respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q S/Q P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.  相似文献   

9.
We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ c and P d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.  相似文献   

10.
The Q-factor estimates of the Earth’s crust and upper mantle as the functions of frequency (Q(f)) are obtained for the seismic S-waves at frequencies up to ~35 Hz. The estimates are based on the data for ~40 earthquakes recorded by the Kislovodsk seismic station since 2000. The magnitudes of these events are MW > 3.8, the sources are located in the depth interval from 1 to 165 km, and the epicentral distances range from ~100 to 300 km. The Q-factor estimates are obtained by the methods developed by Aki and Rautian et al., which employ the suppression of the effects of the source radiation spectrum and local site responses in the S-wave spectra by the coda waves measured at a fixed lapse time (time from the first arrival). The radiation pattern effects are cancelled by averaging over many events whose sources are distributed in a wide azimuthal sector centered at the receiving site. The geometrical spreading was specified in the form of a piecewise-continuous function of distance which behaves as 1/R at the distances from 1 to 50 km from the source, has a plateau at 1/50 in the interval from 50–70 km to 130–150 km, and decays as \({\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 {\sqrt R }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\sqrt R }$}}\) beyond 130–150 km. For this geometrical spreading model and some of its modifications, the following Q-factor estimates are obtained: Q(f) ~ 85f0.9 at the frequencies ranging from ~1 to 20 Hz and Q(f) ~ 75f1.0 at the frequencies ranging from ~1 to 35 Hz.  相似文献   

11.
The relationship between the characteristics of seismic waves in the Western Caucasus and the geological-tectonic structure of the region is studied for identifying the specificity of seismic propagation in the mountainous regions with a complicated geological structure and forecasting the characteristics of the propagation from the geological and tectonic data. The interpretation is presented for the estimates of the Q-factor of the medium (Q(f) ~ 55f0.9 in the region of Sochi and Q(f) ~ 90f0.7 in the region of Anapa), seismic wave enhancement in the upper crustal layers (A(f) ~ 1), and peak ground acceleration residuals, which were previously determined from the records of the local earthquakes and show the distributions of local variations in the parameters of seismic wave radiation and propagation. The obtained characteristics are interpreted in the context of the up-to-date information about the tectonic, geological, and deep structure of the epicentral zones in the Western Caucasus and neighboring territory of the Black Sea. The discrepancies revealed in the low-frequency behavior of the Q-factor in the vicinities of Sochi and Anapa is accounted for by the spatial scale and character of tectonic dislocations of the rocks in these regions. The local variations in the parameters of seismic radiation and propagation are probably related to the geological features of the region such as the fault structures, including the thrusts, shatter zones, oblique seismic boundaries, variations in the thickness and consolidation of the sedimentary cover, as well as the peculiarities in the structure and material composition of the basement.  相似文献   

12.
In log time-frequency spectra, the nonstationary convolution model is a linear equation and thus we improved the Gabor deconvolution by employing a log hyperbolic smoothing scheme which can be implemented as an iteration process. Numerical tests and practical applications demonstrate that improved Gabor deconvolution can further broaden frequency bandwidth with less computational expenses than the ordinary method. Moreover, we attempt to enlarge this method’s application value by addressing nonstationary and evaluating Q values. In fact, energy relationship of each hyperbolic bin (i.e., attenuation curve) can be taken as a quantitative indicator in balancing nonstationarity and conditioning seismic traces to the assumption of unchanging wavelet, which resultantly reveals more useful information for constrained reflectivity inversion. Meanwhile, a statistical method on Q-value estimation is also proposed by utilizing this linear model’s gradient. In practice, not only estimations well agree with geologic settings, but also applications on Q-compensation migration are favorable in characterizing deep geologic structures, such as the pinch-out boundary and water channel.  相似文献   

13.
Coalbed methane can be detected employing the amplitudevariation-with-offset technique. However, there are two issues in applying this technique to a coalbed: strong azimuthal anisotropy resulting from high-density fractures, and the seismic response being composed of many or several individual reflections within the coalbed. To overcome these difficulties, we present an exact solution for reflections in extensive dilatancy anisotropy media. First, we build a three-layer model and simulate the wave propagation in this model. Then we derive an exact P- and converted S-wave reflection coefficient equation based on boundary conditions. Finally, substituting given model parameters into the exact equation, we obtain the variation in the reflection coefficient with incidence angle. The results show that the fracture factors, wavelet frequency and thickness of the coalbed have different effects on the reflection coefficient. Furthermore, we create a synthetic seismogram by forward calculation, and the result fits well with results of the exact equation.  相似文献   

14.
A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used to record vibration graphs. An expression of higher degree of strain (deformation force) is introduced. It is testified theoretically that the displacement, velocity and acceleration cannot be obtained by simple integral acceleration and differential velocity when long displacement and high strain exist, namely when the pile phase generates a whole slip relative to the soil body. That is to say that there are non-linear relations between them. It is educed accordingly that the force P and displacement S are calculated from the amplitude of wave train and (dynamic) P-S curve is drew so as to determine the yield points. Further, a method of determining the vertical bearing capacity for single-pile is discussed. A static load test is utilized to check the result of dynamic test and determine the correlative constants of dynamic-static P(Q)-S curve.  相似文献   

15.
Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F 2 region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F 2 region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.  相似文献   

16.
Based on the Anapa (ANN) seismic station records of ~40 earthquakes (MW > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.  相似文献   

17.
Theory of wave boundary layers (WBLs) developed by Reznik (J Mar Res 71: 253–288, 2013, J Fluid Mech 747: 605–634, 2014, J Fluid Mech 833: 512–537, 2017) is extended to a rotating stratified fluid. In this case, the WBLs arise in the field of near-inertial oscillations (NIOs) driven by a tangential wind stress of finite duration. Near-surface Ekman layer is specified in the most general form; tangential stresses are zero at the lower boundary of Ekman layer and viscosity is neglected below the boundary. After the wind ceases, the Ekman pumping at the boundary becomes a linear superposition of inertial oscillations with coefficients dependent on the horizontal coordinates. The solution under the Ekman layer is obtained in the form of expansions in the vertical wave modes. We separate from the solution a part representing NIO and demonstrate development of a WBL near the Ekman layer boundary. With increasing time t, the WBL width decays inversely proportional to \( \sqrt{t} \) and gradients of fields in the WBL grow proportionally to \( \sqrt{t} \); the most part of NIO is concentrated in the WBL. Structure of the WBL depends strongly on its horizontal scale L determined by scale of the wind stress. The shorter the NIO is, the thinner and sharper the WBL is; the short-wave NIO with L smaller than the baroclinic Rossby scale LR does not penetrate deep into the ocean. On the contrary, for L?≥?LR, the WBL has a smoother vertical structure; a significant long-wave NIO signal is able to reach the oceanic bottom. An asymptotic theory of the WBL in rotating stratified fluid is suggested.  相似文献   

18.
In this paper, the relationships of the plunges and azimuths of T and P axes versus the strikes, dips, and rakes of two seismic nodal planes were derived to provide reference for earthquake researchers. The independence of the plunges and azimuths of T, B, and P axes in focal mechanism solution was discussed, and it was concluded that three parameters, i.e., the azimuths of T, B and P axes, are completely independent. The focal mechanism solution representation based on Euler rotation was introduced, using three Euler angles in place of the plunges and azimuths of T, B, and P axes, and three focal mechanism solution representations were briefly compared and analyzed in respect of accuracy on the basis of the assumption of rounding; it was concluded that the Euler angle representation has better accuracy, compared with the azimuth representation and the traditional representation with T, B, and P axes.  相似文献   

19.
Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s???1 (for Rayleigh waves) and 8 ×10???3 s???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.  相似文献   

20.
This paper presents a stochastic model for multicomponent competitive monovalent cation exchange in hierarchical porous media. Reactive transport in porous media is highly sensitive to heterogeneities in physical and chemical properties, such as hydraulic conductivity (K), and cation exchange capacity (CEC). We use a conceptual model for multimodal reactive mineral facies and develop a Eulerian-based stochastic theory to analyze the transport of multiple cations in heterogeneous media with a hierarchical organization of reactive minerals. Numerical examples investigate the retardation factors and dispersivities in a chemical system made of three monovalent cations (Na+, K+, and Cs+). The results demonstrate how heterogeneity influences the transport of competitive monovalent cations, and highlight the importance of correlations between K and CEC. Further sensitivity analyses are presented investigating how the dispersion and retardation of each cation are affected by the means, variances, and integral scales of K and CEC. The volume fraction of organic matter is shown to be another important parameter. The Eulerian stochastic framework presented in this work clarifies the importance of each system parameters on the migration of cation plumes in formations with hierarchical organization of facies types. Our stochastic approach could be used as an alternative to numerical simulations for 3D reactive transport in hierarchical porous media, which become prohibitively expensive for the multicomponent applications considered in this work.  相似文献   

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

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