首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 296 毫秒
1.
Conventional f?x empirical mode decomposition (EMD) is an effective random noise attenuation method for use with seismic profiles mainly containing horizontal events. However, when a seismic event is not horizontal, the use of f?x EMD is harmful to most useful signals. Based on the framework of f?x EMD, this study proposes an improved denoising approach that retrieves lost useful signals by detecting effective signal points in a noise section using local similarity and then designing a weighting operator for retrieving signals. Compared with conventional f?x EMD, f?x predictive filtering, and f?x empirical mode decomposition predictive filtering, the new approach can preserve more useful signals and obtain a relatively cleaner denoised image. Synthetic and field data examples are shown as test performances of the proposed approach, thereby verifying the effectiveness of this method.  相似文献   

2.
This study considers the possibility of using the new methods of time-frequency transforms, such as chirplet and warblet transforms, to analyze the digital observational data of geomagnetic pulsations of Pc5 type. For this purpose, necessary algorithms of calculation and appropriate software were developed. The chirplet transform method (CT) is used to analyze signals with a linear frequency modulation. A chirplet variation, the so-called warblet transform, is used to analyze signals with a nonlinear frequency modulation. Since, in studying geomagnetic pulsations, it is difficult to make assumptions on the character of the behavior of the instantaneous frequency of the signal, the special generalized warblet transform (GWT) was used for the analysis. The GWT has a high spatiotemporal resolution and was developed to analyze oscillations both with a periodic and nonperiodic change of the instantaneous frequency. The software developed for GWT calculation was used to study daytime geomagnetic Pc5 pulsations with durations of several hours that were detected via the network of ground-based magnetometers of the Scandinavian IMAGE profile during the magnetic storm of May 29–30, 2003. For the first time, temporal variations of the instantaneous frequency of geomagnetic pulsations are determined and their possible use in studying the fine spatial structure of Pc5 waves is shown.  相似文献   

3.
The effect of baric variations of different origins on characteristics of seismic noise is analyzed in the frequency range 0.03–20 Hz. Long period variations in atmospheric pressure caused by cyclones, whose period T ranges from half a day to a few days, are shown to increase the microseismic background amplitude by two to four times in the frequency range 0.03–1 Hz (the coefficient of linear correlation between time variations in the amplitude and atmospheric pressure is K = 0.65 at a significance level of r = 0.95). Short-period baric variations with T ～ 5–30 min associated with the passage of cold fronts lead to a tenfold increase in the microseismic background amplitude in the frequency range 4–8 Hz (K = 0.67 at r = 0.95). In this case, disturbances of seismic background are recorded for 20–60 min after the passage of an atmospheric front and display an exponential drop in the amplitude. In distinction to cyclones, an atmospheric front increases the number of impulsive microseismic events of the resonance type.  相似文献   

4.
In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q S (f)?=?(147?±?16)f 0.71?±?0.02 and (126?±?12)f 0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).  相似文献   

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

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

7.
  相似文献   

8.
  相似文献   

9.
Attenuation of seismic compression waves leads to the real existence of a fast P1 wave in rocks which are fully saturated with dropping fluid and a slow P2 wave in the rocks containing gas in their pores. This accounts for the seismic blanking zones below the gas horizons for the P1 waves. Oscillations of gaseous inclusions ensure the energy transfer to the dominant frequencies which are different for the cases of passive seismic (few Hz) and active source seismic (10–20 Hz). The intervals of dominant frequencies are determined from the negative attenuation of these low-frequency waves. According to the observations and the suggested equation, random noise amplifies the signal at these frequencies. Thus, the P2 waves at the dominant frequency of the active source seismics are applicable for elaborating on the details of the saturation of the production layer by hydrocarbons. The relation to the AVO method (Amplitude Variation with Offset) and dilatancy effect during the preparation of an earthquake is noted.  相似文献   

10.
Seismic signal windowing is the preliminary step for many analysis procedures in engineering seismology (standard spectral ratio, quality factor, general inversion techniques, etc.). Moreover a noise window is often necessary for the data quality control through the signal-to-noise verification. Selecting the noise window can be challenging when large heterogeneous datasets are considered, especially when they include short pre-event noise signals. This study proposes a fully automatic and configurable (i.e., with default parameters that can also be user-defined) algorithm to windowing the noise and the P, S, coda and full signal once the P-wave (T P ) and S-wave (T S ) first arrivals are known. An application example is given on a KiK-net dataset. A Matlab language implementation of this algorithm is proposed as an online resource.  相似文献   

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

12.
A method for constructing the empirical model of the F2 layer critical frequency (foF2) under magnetically quiet conditions, aimed at analyzing disturbances of any nature, is proposed. This method has been analyzed, and typical features of regular changes in foF2 of the quiet ionosphere and day-to-day foF2 variability are analyzed using the data from Irkutsk and Slough stations as an example. In particular, it has been obtained that this model differs from the international IRI model, and this difference is mainly caused by the fact that the foF2 values in the IRI model do not correspond to quiet conditions. Therefore, this model gives a larger amplitude of the annual and semiannual variations in foF2 than the IRI model. In addition, this model more accurately reproduces the rate of foF2 annual variations at a fixed local time, especially in equinoxes, when foF2 variations can exceed 1 MHz within one month.  相似文献   

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

14.
The existence of strong random noise in surface microseismic data may decrease the utility of these data. Non‐subsampled shearlet transform can effectively suppress noise by properly setting a threshold to the non‐subsampled shearlet transform coefficients. However, when the signal‐to‐noise ratio of data is low, the coefficients related to the noise are very close to the coefficients associated with signals in the non‐subsampled shearlet transform domain that the coefficients related to the noise will be retained and be treated as signals. Therefore, we need to minimise the overlapping coefficients before thresholding. In this paper, a singular value decomposition algorithm is introduced to the non‐subsampled shearlet transform coefficients, and low‐rank approximation reconstructs each non‐subsampled shearlet transform coefficient matrix in the singular value decomposition domain. The non‐subsampled shearlet transform coefficients of signals have bigger singular values than those of the random noise, which implies that the non‐subsampled shearlet transform coefficients can be well estimated by taking only a few largest singular values. Therefore, those properties of singular value decomposition may significantly help minimise overlapping of noise and signals coefficients in the non‐subsampled shearlet transform domain. Finally, the denoised microseismic data are obtained easily by giving a simple threshold to the reconstructed coefficient matrix. The performance of the proposed method is evaluated on both synthetic and field microseismic data. The experimental results illustrate that the proposed method can eliminate random noise and preserve signals of interest more effectively.  相似文献   

15.
The 40-year period of observations of short-term variations (with characteristic times of up to 1–2 days) in the critical frequency of the ionospheric F2 layer (foF2) is analyzed. The continuous (with a step of 1 h) series of fluctuations (F) of the foF2 critical frequency (with eliminated daily variations) has been calculated using the hourly variations in foF2 at Moscow stations. The fractal dimension (FRH) of the fluctuations, characterizing short-term variations in foF2, has been determined and analyzed on a 30-day interval, using the Higuchi method. It has been established that FRH estimates substantially change in time. The 11-year cycle, which is in antiphase with the solar cycle, and the total annual and semiannual variations, similar to the variations observed in the normalized critical frequency of the E region and in the electron density of the D region, are clearly defined in these changes. Thus, the parameters of fast variations in the ionospheric F2 layer are affected by the phase of the 11-year solar cycle and by the position of the Earth in the orbit or seasonal variations in the atmosphere.  相似文献   

16.
Using the foF2 database obtained from satellites and ground-based ionospheric stations, we have constructed a global empirical model of the critical frequency of the ionospheric F2-layer (SDMF2—Satellite and Digisonde Data Model of the F2 layer) for quiet geomagnetic conditions (Kp < 3). The input parameters of this model are the geographical coordinates, UT, day, month, year, and the integral index F10.7 (day, τ = 0.96) of solar activity for a given day. The SDMF2 model was based on the Legendre method for the spatial expansion of foF2 monthly medians to 12 in latitude and 8 in longitude of spherical harmonics. The resulting spatial coefficients have been expanded by the Fourier method in three spherical harmonics with respect to UT. The effect of the saturation of critical frequency of the ionospheric F2-layer at high solar activity was described in the SDMF2 model by foF2 as a logarithmic function of F10.7 (day, τ = 0.96). The difference between the SDMF2 and IRI models is a maximum at low solar activity as well as in the Southern Hemisphere and in the oceans. The testing on the basis of ground-based and satellite data has indicated that the SDMF2 model is more accurate than the IRI model.  相似文献   

17.
The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b m , N, M 1, M 2) that an observed b m estimate, based on a sample of N magnitudes within an [M 1????≤?ΔM/2,?M 2?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.  相似文献   

18.
Floods have changed in a complex manner, triggered by the changing environment (i.e., intensified human activities and global warming). Hence, for better flood control and mitigation in the future, bivariate frequency analysis of flood and extreme precipitation events is of great necessity to be performed within the context of changing environment. Given this, in this paper, the Pettitt test and wavelet coherence transform analysis are used in combination to identify the period with transformed flood-generating mechanism. Subsequently, the primary and secondary return periods of annual maximum flood (AMF) discharge and extreme precipitation (Pr) during the identified period are derived based on the copula. Meanwhile, the conditional probability of occurring different flood discharge magnitudes under various extreme precipitation scenarios are estimated using the joint dependence structure between AMF and Pr. Moreover, Monte Carlo-based algorithm is performed to evaluate the uncertainties of the above copula-based analyses robustly. Two catchments located on the Loess plateau are selected as study regions, which are Weihe River Basin (WRB) and Jinghe River Basin (JRB). Results indicate that: (1) the 1994–2014 and 1981–2014 are identified as periods with transformed flood-generating mechanism in the WRB and JRB, respectively; (2) the primary and secondary return periods for AMF and Pr are examined. Furthermore, chance of occurring different AMF under varying Pr scenarios also be elucidated according to the joint distribution of AMF and Pr. Despite these, one thing to notice is that the associate uncertainties are considerable, thus greatly challenges measures of future flood mitigation. Results of this study offer technical reference for copula-based frequency analysis under changing environment at regional and global scales.  相似文献   

19.
The Aki-Utsu maximum likelihood method is widely used for estimation of the Gutenberg-Richter b-value, but not all authors are conscious of the method’s limitations and implicit requirements. The Aki/Utsu method requires a representative estimate of the population mean magnitude; a requirement seldom satisfied in b-value studies, particularly in those that use data from small geographic and/or time windows, such as b-mapping and b-vs-time studies. Monte Carlo simulation methods are used to determine how large a sample is necessary to achieve representativity, particularly for rounded magnitudes. The size of a representative sample weakly depends on the actual b-value. It is shown that, for commonly used precisions, small samples give meaningless estimations of b. Our results give estimates on the probabilities of getting correct estimates of b for a given desired precision for samples of different sizes. We submit that all published studies reporting b-value estimations should include information about the size of the samples used.  相似文献   

20.
This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f 0(P)/f 0(S). The ratios were computed for each event with corner frequencies f 0 of P- and S-wave, which were obtained from the measured f 0 I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f 0(P)/f 0(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f 0(P)/f 0(S), with the empirical threshold value of the ratio for Israel as 1.48.  相似文献   

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

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