首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
Incorrect results including false low frequency components may appear during the recovery of ground motion in the frequency domain. In order to check the quality of restoration in the frequency domain of CDSN broad-band seismograms, a digital model of the seismograph systems, based on modern control engineering, is presented and improvement of the method of restoration in the frequency domain is described. Examples shown in this paper demonstrate the effectiveness of the proposed approach in improving the restoration of true ground motion. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 90–99, 1992. The work is supported by the Chinese Joint Seismological Science Foundation.  相似文献   

2.
A method is presented to derive approximate versions of the wave equation which allow finite-difference migration for very steep dips (> 50°). It is shown that for conventional finite-difference schemes, in addition to the dip limitation, the maximum acceptable frequency should be specified. A finite-difference migration technique is proposed in the frequency domain. It is derived that finite-difference wave field extrapolation in the frequency domain consists of a space-variant convolution procedure for each frequency component, the space-variance being defined by the lateral variation in the velocity. Finally it is shown that with finite-difference migration, particle velocity data can be easily obtained from pressure data.  相似文献   

3.
Summary The reasoning that joule heating in multiple-stroke lightning flashes associated with an active tornado constitutes the major energy source for itsdrive is questioned. An electrical discharge is proposed as a mechanism for the formation of the initial vortexsink, a high specific-energy source being supplied by Joule heating in a leader-stroke coronal discharge of several hundred amperes and milliseconds duration. Magnetic pinch effects are invoked as a means of preventing large lateral heat losses from the ionised column and a resulting degradation of the energy source. Applying the equation for magneto-hydrostatic balance, it is shown that magnetic pinch exerted on the column by the current-induced magnetic field alone is insufficient to achieve balance and thus prevent energy losses. Qualitative arguments are used to demonstrate that the pinch effect can be enhanced through the interaction with the external, geomagnetic field. Observational evidence is presented to support the hypothesis of geomagnetic control for tornado initiation on a global scale. It is seen in large differences of tornado frequency between regions on the surface of the Earth featured by comparable pre-requisite thunderstorm frequency but significant differences in geomagnetic intensity which qualitatively agree with the predicted effect.  相似文献   

4.
Attenuation in seismic wave propagation is a common cause for poor illumination of subsurface structures. Attempts to compensate for amplitude loss in seismic images by amplifying the wavefield may boost high‐frequency components, such as noise, and create undesirable imaging artefacts. In this paper, rather than amplifying the wavefield directly, we develop a stable compensation operator using stable division. The operator relies on a constant‐Q wave equation with decoupled fractional Laplacians and compensates for the full attenuation phenomena by performing wave extrapolation twice. This leads to two new imaging conditions to compensate for attenuation in reverse‐time migration. A time‐dependent imaging condition is derived by applying Q‐compensation in the frequency domain, whereas a time‐independent imaging condition is formed in the image space by calculating image normalisation weights. We demonstrate the feasibility and robustness of the proposed methods using three synthetic examples. We found that the proposed methods are capable of properly compensating for attenuation without amplifying high‐frequency noise in the data.  相似文献   

5.
We have developed a novel method for missing seismic data interpolation using f‐x‐domain regularised nonstationary autoregression. f‐x regularised nonstationary autoregression interpolation can deal with the events that have space‐varying dips. We assume that the coefficients of f‐x regularised nonstationary autoregression are smoothly varying along the space axis. This method includes two steps: the estimation of the coefficients and the interpolation of missing traces using estimated coefficients. We estimate the f‐x regularised nonstationary autoregression coefficients for the completed data using weighted nonstationary autoregression equations with smoothing constraints. For regularly missing data, similar to Spitz f‐x interpolation, we use autoregression coefficients estimated from low‐frequency components without aliasing to obtain autoregression coefficients of high‐frequency components with aliasing. For irregularly missing or gapped data, we use known traces to establish nonstationary autoregression equations with regularisation to estimate the f‐x autoregression coefficients of the complete data. We implement the algorithm by iterated scheme using a frequency‐domain conjugate gradient method with shaping regularisation. The proposed method improves the calculation efficiency by applying shaping regularisation and implementation in the frequency domain. The applicability and effectiveness of the proposed method are examined by synthetic and field data examples.  相似文献   

6.
It is shown how the empirical equations for scaling the Fourier amplitude spectra in the frequency band from ～0.1 to 25 Hz can be extended to describe the strong motion amplitudes in a much broader frequency range. At long periods, the proposed equations are in excellent agreement with (1) the seismological and field estimates of permanent ground displacement (near field) and (2) the independent estimates of seismic moment (far field). At high frequencies, f ≥ 25 Hz, the spectral amplitudes can be described by exp (? πkf), where k ranges from 0·02 (near source) to about 0·06 at an epicentral distance of about 200 km. It is also shown how amplification by local soil and geological site conditions can be defined to apply in the same broad frequency range.  相似文献   

7.
In this paper two causal models that approximate the nearly frequency‐independent cyclic behaviour of soils are analysed in detail. The study was motivated by the need to conduct time‐domain viscoelastic analysis on soil structures without adopting the ad hoc assumption of Rayleigh damping. First, the causal hysteretic model is introduced in which its imaginary part is frequency independent the same way that is the imaginary part of the popular non‐causal constant hysteretic model. The adoption of an imaginary part that is frequency independent even at the zero‐frequency limit, in conjunction with the condition that the proposed model should be causal, yields a real part that is frequency dependent and singular at zero frequency. The paper shows that the causal hysteretic model, although pathological at the static limit, is the mathematical connection between the non‐causal constant hysteretic model and the physically realizable Biot model. The mathematical structure of the two causal models is examined and it is shown that the causal hysteretic model is precisely the high‐frequency limit of the Biot model. Although both models have a closed‐form time‐domain representation, only the Biot model is suitable for a time‐domain viscoelastic analysis with commercially available computer software. The paper demonstrates that the simplest, causal and physically realizable linear hysteretic model that can approximate the cyclic behaviour of soil is the Biot model. The proposed study elucidates how the dynamic analysis of soil structures can be conducted rigorously in terms of the viscoelastic properties of the soil material and not with the ad hoc Rayleigh damping approach which occasionally has been criticized that tends to overdamp the higher vibration modes. The study concludes that under pulse‐type motions the Rayleigh damping approximation tends to overestimate displacements because of the inappropriate viscous type of dissipation that is imposed. Under longer motions that induce several cycles, the concept of equivalent viscous damping is more appropriate and the Rayleigh damping approximation results to a response that is comparable to the response computed with a rigorous time‐domain viscoelastic finite element analysis. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

8.
The time‐invariant gain‐limit‐constrained inverse Q‐filter can control the numerical instability of the inverse Q‐filter, but it often suppresses the high frequencies at later times and reduces the seismic resolution. To improve the seismic resolution and obtain high‐quality seismic data, we propose a self‐adaptive approach to optimize the Q value for the inverse Q‐filter amplitude compensation. The optimized Q value is self‐adaptive to the cutoff frequency of the effective frequency band for the seismic data, the gain limit of the inverse Q‐filter amplitude compensation, the inverse Q‐filter amplitude compensation function, and the medium quality factor. In the processing of the inverse Q‐filter amplitude compensation, the optimized Q value, corresponding gain limit, and amplitude compensation function are used simultaneously; then, the energy in the effective frequency band for the seismic data can be recovered, and the seismic resolution can be enhanced at all times. Furthermore, the small gain limit or time‐variant bandpass filter after the inverse Q‐filter amplitude compensation is considered to control the signal‐to‐noise ratio, and the time‐variant bandpass filter is based on the cutoff frequency of the effective frequency band for the seismic data. Synthetic and real data examples demonstrate that the self‐adaptive approach for Q value optimization is efficient, and the inverse Q‐filter amplitude compensation with the optimized Q value produces high‐resolution and low‐noise seismic data.  相似文献   

9.
Three‐dimensional receiver ghost attenuation (deghosting) of dual‐sensor towed‐streamer data is straightforward, in principle. In its simplest form, it requires applying a three‐dimensional frequency–wavenumber filter to the vertical component of the particle motion data to correct for the amplitude reduction on the vertical component of non‐normal incidence plane waves before combining with the pressure data. More elaborate techniques use three‐dimensional filters to both components before summation, for example, for ghost wavelet dephasing and mitigation of noise of different strengths on the individual components in optimum deghosting. The problem with all these techniques is, of course, that it is usually impossible to transform the data into the crossline wavenumber domain because of aliasing. Hence, usually, a two‐dimensional version of deghosting is applied to the data in the frequency–inline wavenumber domain. We investigate going down the “dimensionality ladder” one more step to a one‐dimensional weighted summation of the records of the collocated sensors to create an approximate deghosting procedure. We specifically consider amplitude‐balancing weights computed via a standard automatic gain control before summation, reminiscent of a diversity stack of the dual‐sensor recordings. This technique is independent of the actual streamer depth and insensitive to variations in the sea‐surface reflection coefficient. The automatic gain control weights serve two purposes: (i) to approximately correct for the geometric amplitude loss of the Z data and (ii) to mitigate noise strength variations on the two components. Here, Z denotes the vertical component of the velocity of particle motion scaled by the seismic impedance of the near‐sensor water volume. The weights are time‐varying and can also be made frequency‐band dependent, adapting better to frequency variations of the noise. The investigated process is a very robust, almost fully hands‐off, approximate three‐dimensional deghosting step for dual‐sensor data, requiring no spatial filtering and no explicit estimates of noise power. We argue that this technique performs well in terms of ghost attenuation (albeit, not exact ghost removal) and balancing the signal‐to‐noise ratio in the output data. For instances where full three‐dimensional receiver deghosting is the final product, the proposed technique is appropriate for efficient quality control of the data acquired and in aiding the parameterisation of the subsequent deghosting processing.  相似文献   

10.
Two types of implicit time-stepping algorithms have been proposed recently for pseudodynamic tests. The first type consists of an algorithm which relies on Newton iterations to satisfy the equations of motion. The second type consists of an algorithm which is based on the Operator-Splitting technique and does not require any numerical iteration. While one or the other has been preferred by some researchers, these time-stepping algorithms have not been analysed and compared under a uniform setting. In this paper, a concise summary of these schemes is presented, and they are evaluated in a consistent manner in terms of numerical dissipation, frequency distortion and experimental errors. The analytical results are validated by numerical simulations as well as experimental results. It is shown that the algorithm based on Newton iterations can control experimental error effects effectively by means of an error-correction procedure. The algorithm based on the Operator-Splitting technique demonstrates similar performance provided the I-Modification is adopted.  相似文献   

11.
A set of monoial scaling relations to parametrize several measures of strong motion (peak velocity, peak starting-phase acceleration, peak stopping acceleration) is proposed. Dynamic solutions are obtained for a 2-D (P-SV) stress-drop model of faulting, and ground motion from these calculations is used to calibrate the scaling relations. Geometrical spreading, radiation patterns, low frequency near-field radiation, and free surface response are analysed, and introduced as corrections. The calculational finite difference method is sound within its frequency range of validity, which is found to be about 0 to 5–6 Hz for the chosen grid steps. A strong difference is obtained between theP andSV motion scalings, mainly with the source rupture velocity. Also noted are significant differences between the starting and stopping accelerations due to different frequency content and the influence of low frequency near-field radiation. To test the estimated scaling relationships, some synthetic predictions of the kinematic parameters are made, with quite good agreement when compared with dynamic computations (errors within 30%). The results emphasize some features of the near-fieldP-SV radiation and allow a kinematic prediction for a simple and smooth source model, but show the limited reliability of such predictions, arising from the problem's complexity.  相似文献   

12.
Summary Variation ofb with depth has been studied using most likelihood method of estimatingb proposed byUtsu and it has been found thatb remains practicall constant upto a depth of 400 km, at which there is a sharp increase in the value ofb which decreases with depth. The distribution of Maximum magnitude with depth shows that the maximum magnitude decreases upto 300 km and below this level there is an increase in the maximum magnitude which again decreases with depth. The frequency of distribution ofb values is found to cluster around 0.8.  相似文献   

13.
The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (Mw) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (Mw) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.  相似文献   

14.
Wavefield depth extrapolation and prestack depth migration in complex anelastic media are studied. Kjartansson's frequency‐independent Q law is used to describe the absorption of seismic energy. The macromodel used is analogous to the macromodel used for current migration schemes except that an additional frequency‐independent Q macromodel needs to be provided. Absorption in the forward one‐way propagator is introduced by assuming a complex phase velocity, and the inverse one‐way propagator is obtained using the reciprocity theorem for one‐way wavefields in dissipative media. The stability of the inverse propagator is achieved by limiting the angle of propagation of wavefields. A table‐driven explicit operator scheme for imaging complex 2D anelastic media is presented. High‐accuracy, short convolution operators are designed by the weighted least‐squares method, and two kinds of imaging conditions are proposed. Numerical examples of depth extrapolation in laterally varying media, the migration of a spatial impulse with dispersion as well as shot record depth migration demonstrate the potential of the proposed explicit forward operator, the explicit inverse operator and the prestack depth migration scheme, respectively.  相似文献   

15.
In this paper, a new index is proposed for the selection of the best regional frequency analysis method. First, based on the theory of reliability, the new selective index is developed. The variances of three regional T‐year event estimators are then derived. The proposed methodology is applied to an actual watershed. For each regional method, the reliability of various T‐year regional estimates is computed. Finally, the reliability‐based selective index graph is constructed from which the best regional method can be determined. In addition, the selection result is compared with that based on the traditional index, root mean square error. The proposed new index is recommended as an alternative to the existing indices such as root mean square error, because the influence of uncertainty and the accuracy of estimates are considered. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

16.
Empirical frequency distributions of multiplicative cascade weights, or breakdown coefficients, at small timescales are analyzed for 5-min precipitation time series from four gauges in Germany. It is shown that histograms of the weights, W, are strongly deformed by the recording precision of rainfall amounts. A randomization procedure is proposed to statistically remove the artifacts due to precision errors in the original series. Evolution of the probability distributions of W from beta-like for large timescales to combined beta-normal distribution with a pronounced peak at W ≈ 0.5 for small timescales is observed. A new 3N-B distribution built from 3 separate normal, N, distributions and one beta, B, distribution is proposed for reproduction of the empirical histograms of W at small timescales. Parameters of the 3N-B distributions are fitted for all gauges and analyzed timescales. Microcanonical cascades models with a generator based on 3N-B distributions are developed and their performance at disaggregating precipitation at 1280-min intervals down to 5-min intervals is evaluated.  相似文献   

17.
A simple theory is presented to derive explicit recursive migration schemes in three dimensions which are suitable for situations with lateral velocity variations. The proposed schemes can be realized by one-dimensional convolutions along the x-axis and y-axis. An investigation is made on the errors of each scheme. For each frequency component the amplitude error and the phase error are computed as a function of the immergence angle (dip). Based on the error analysis a critical discussion is given on the stability of explicit schemes.  相似文献   

18.
A new method is proposed for the quantitative interpretation of SP field data produced by a polarized ore body and simulated by an inclined sheet. The theoretical concept is based on the study of the amplitude Spectrum. It is shown that the SP amplitude Spectrum is not continuous at zero frequency; this leads to the dip angle determination of the inclined sheet. It is also shown that the SP amplitude spectrum is practically nullified at a characteristic amortization frequency that depends on the depth of the polarized body. The maximum amplitude Spectrum value of the SP gradient is used to estimate the depth to the bottom of the polarized body. Thus, the geometrical parameters h and H, the depths to the top and bottom, respectively, as well as the dip angle of the inclined sheet, can be satisfactorily determined. Some problems may arise in the determination of these parameters, affecting their accuracy, whenever unwanted frequency noise is present.  相似文献   

19.
利用山东台网记录的长岛震群2017年2月14日—9月1日期间的波形与震相资料研究长岛地区非弹性衰减系数,得到该地区介质平均Q值与频率f的关系式为Q(f)=363.9f1.374 1。采用Moya等[1]提出的利用遗传算法联合反演得到长岛周边台站的场地响应,根据Brune模型震源参数计算公式求解长岛震群序列地震震源参数。结果显示,各个震源参数之间均存在一定的相关关系,地震矩随ML震级的增大而增加,地震矩与破裂半径R之间存在半对数关系,拐角频率fc随地震矩的增大而减少;长岛地震序列的应力降数值普遍偏小,最大不超过0.9MPa,这意味着长岛震源区整体构造应力较低,也可能指示长岛震群为低摩擦应力的断层作用;震源参数随时间的变化方面,整体而言,长岛震群地震应力降变化起伏很大,在M4.1地震发生前,拐角频率与应力降均发生快速下降后随即翻转上升的现象,证明在M4.1地震发生前震源区整体应力的挤压逐渐增强。  相似文献   

20.
Assessing the probability of collapse is a computationally demanding component of performance‐based earthquake engineering. This paper examines various aspects involved in the computation of the mean annual frequency of collapse (λc) and proposes an efficient method for estimating the sidesway collapse risk of structures in seismic regions. By deaggregating the mean annual frequency of collapse, it is shown that the mean annual frequency of collapse is typically dominated by earthquake ground motion intensities corresponding to the lower half of the collapse fragility curve. Uncertainty in the collapse fragility curve and mean annual frequency of collapse as a function of the number of ground motions used in calculations is also quantified, and it is shown that using a small number of ground motions can lead to unreliable estimates of a structure's collapse risk. The proposed method is shown to significantly reduce the computational effort and uncertainty in the estimate. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

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

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