共查询到20条相似文献，搜索用时 31 毫秒
1.
用Q值刻画的地震衰减在地震信号处理和解释中具有很广泛的应用。利用反射地震资料进行Q值估计需要解决地震子波和反射系数序列耦合的问题。从反射地震资料中去除反射系数序列的影响,这个过程称为频谱校正。本文提出了一种基于子波估计的求取Q值的方法,进而设计了一个反Q滤波器。该方法利用反射地震资料的高阶统计量进行子波估计,并利用所估计子波实现频谱校正。我们利用合成数据实验给出了质心频移法与频谱比法这两种常用的Q值估计方法在不同参数设置下的性能。人工合成数据和实际数据处理表明,利用本文提出的方法进行频谱校正后,可以得到可靠的Q值估计。经过反Q滤波,地震数据的高频部分得到了有效地恢复。 相似文献
2.
A. Pirogova R. Pevzner B. Gurevich S. Glubokovskikh K. Tertyshnikov 《Geophysical Prospecting》2019,67(7):1778-1797
At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO2-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO2 plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Qscat ≈ 50–65. 相似文献
3.
In this paper, the "spectral amplitude ratio method" (SAR), "energy method" (EN) and "coda wave method" (CW) are used to calculate
theQ value variations of gneiss in the preparing rupture process. The obtained results show that the variation state ofQ values by SAR features the shape of relative stability—gradual increment to the maximum—then decrement and final rupture.
The variation state ofQ values by EN is just contrary to that by SAR, i. e. with the shape of stability—decrement—increment—and final rupture. The
varation state ofQ values by CW is similar to that by EN, its main frequency features the shape of relatively high value—decrement to the minimum—increment—and
final rupture. But to the high frequency (higher than the main frequency), the variation state ofQ values features the shape of the stable value-increment to the maximum-decrement-and final rupture. At the same time, the
results by coda wave amplitude spectrum show that, when stress reaches 70% of rupture stress, the high frequency component
of S wave rapidly reduces (Q
c
increasing); at the time of impending the main rupture, the main frequency component reduces with a large scale (Q
c
increasing again), this may be the reason which causes the different variation states of two codaQ values. The result of amplitude spectra of P, S (initial wave) waves also show that with the appearance of microcracks the
frequency band of S wave turn to be narrow, the high frequency component is reduced quickly, i. e. the S wave spectra have
different variation states with different frequency components. That is why theQ
s
obtained by different methods have different variation characteristics. 相似文献
4.
本文利用云南地区永胜台观测的地磁南北和东西分量开展舒曼谐振的背景变化特征分析.通过提取舒曼谐振各阶频率和功率谱密度的小时均值,分析了前三阶舒曼谐振频率和功率谱密度在分点和至点前后的周日变化特征.舒曼谐振功率谱密度的周日变化与亚洲、非洲和美洲三大闪电活动中心的活跃时段以及观测站相对于三大闪电活动中心的方位密切相关.舒曼谐振频率的周日变化特征更复杂.各阶功率谱密度和频率在夏至和秋分前后的变化幅度比春分和冬至前后大.从2011年舒曼谐振频率和功率谱密度日中值的年变化图中发现,谐振频率随季节变化的特征不明显,而功率谱密度的年度变化曲线呈半周期正弦波形态,以7月份为轴对称分布.功率谱密度的季节变化特征与闪电活动的季节变化特征相一致.南北和东西分量得到的前三阶谐振频率,第一阶约稳定在7.5 Hz.而随着阶数增加,南北分量得到的谐振频率比东西分量约大0.5 Hz.高阶谐振频率发生偏移的原因目前还不清楚. 相似文献
5.
We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q
p
), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q
p
with an interface at 660 km depth and take Q
p
constant with period. Our data set is too small to solve for a possible frequency dependence of Q
p
. The upper mantle Q
1 is 476 [299–1176] and the lower mantle Q
2 is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q
p
–1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135. 相似文献
6.
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. 相似文献
7.
RAINER TONN 《Geophysical Prospecting》1991,39(1):1-27
Ten methods for the computation of attenuation have been investigated, namely: amplitude decay, analytical signal, wavelet modelling, phase modelling, frequency modelling, rise-time, pulse amplitude, matching technique, spectral modelling and spectral ratio. In particular, we have studied the reliability of each of these methods in estimating correct values of Q using three synthetic VSP seismograms for plane P-waves with different noise contents. The investigations proved that no single method is generally superior. Rather, some methods are more suitable than others in specific situations depending on recording, noise or geology. The analytical signal method has been demonstrated to be superior if true amplitude recordings are available. Otherwise spectral modelling or, in the ‘ noise-free’ case the spectral ratio method, is optimal. Finally, two field VSPs in sediments are investigated. Only in the case of the highest quality VSP can significant information be deduced from the computed attenuation. 相似文献
8.
V. V. Pchelkin M. I. Beloglazov A. N. Vasiliev A. I. Voronin 《Geomagnetism and Aeronomy》2010,50(5):623-631
The observation results of Q-type bursts in the measurements of the horizontal component of the noise magnetic field in the range of the first Schumann
resonance in polar regions (Lovozero high-latitude observatory) are presented. Automatic selection of Q-type bursts from the experimental data series is implemented on the basis of a waveform recognition algorithm. The resonant
nature of Q-bursts is shown. The possibility of selecting such events in magnetic excitation conditions is highlighted. The global resonator
quality upon decreasing the selected waveform amplitudes is estimated. The data obtained by this method are compared with
estimates on the basis of Fourier analysis and values known from the world literature. The possible reasons for disagreements
of the estimates are analyzed, including the problem of selecting the function approximating the spectrum, the problem of
accounting for the background, and the possible irregularity of the spectrum. It is shown that Q-type bursts, besides the quality, allow estimating the resonant frequency of the first oscillation mode; however, the accuracy
of such estimation is lower as compared to the results of Fourier analysis methods. 相似文献
9.
10.
The quality factor Q is a vital parameter for quantitatively describing the attenuation information of underground reservoirs, which is of great significance for hydrocarbon detection and reservoir characterization. A frequency-weighted-exponential (FWE) method utilizing the symmetry index and the characteristic frequency can obtain this parameter. Unfortunately, the constant symmetry index assumption of it reduces the accuracy of Q values under the non-standard FWE shape. Selecting an optimal symmetry index is also a problem for this method. Hence, the basic idea of a novel Q estimation method is to substitute the symmetry index with the standard deviation of the source and attenuated wavelet spectrums. The added standard deviation varies with the degree of attenuation under different wavelet shapes, which can reduce the effects of the constant assumption. Meanwhile, it is simple to calculate this parameter directly from the spectrums. In this way, the proposed method has wider applicability to various wavelets. The synthetic records show the better performance of the novel method in improving the accuracy of Q values and the resistance to random noise than the FWE method. Furthermore, the results of matching wavelets exhibit the applicability of the proposed method in real data, and the quadratic spectrum simulation method improves the stability of the spectrums. Finally, real data experiments indicate the effectiveness of the proposed method after the above two processing means. 相似文献
11.
Coda site amplification factors are used to eliminate the site effect from records of three circum-PacificT phases recorded by the Hawaiian Volcano Observatory (HVO) network on the island of Hawaii. ObservedT-phase amplitudes across the island generally decrease with increasing distance from the conversion point where acoustic waves
in the SOFAR channel become seismic waves propagating through the crust. However, the decay of the observedT-phase signal across the island is not regular in regions of dense station coverage, in particular, the Kilauea caldera region.
We divide the maximum observedT-phase amplitudes at a given station by the coda site amplification factor obtained for the same station and frequency band
(3.0Hz); the distribution of these amplitudes reveals a smooth pattern over the entire island. The distance over which the
site effect-correctedT-phase amplitude decreases by one-half, combined with the apparent velocity of propagation ofT phases across the island, allows for an approximate determination of near-surfaceQ over much of the island of Hawaii. We found a region of lowQ in the Kilauea summit area (Q≈30) and east rift zone (Q≈60) with considerably higherQ in the Kaoiki and northern portions of the island (Q≈150 to 200). The lowQ values obtained in the Kilauea summit region and east rift zone are significantly lower than estimates of codaQ in the same region, suggestingT phases may be sampling the earth's near-surface properties. 相似文献
12.
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. 相似文献
13.
3-D S-waveQ structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in
this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal
faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen
earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude
spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time
overQ, provided that earthquake sources conform to Brune’s disk dislocation model. Finally, by 3-D ray tracing method, theoretical
travel time overQ is made to fit observed travel time overQ with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel
time tomography, remains unchanged, while onlyQ model is modified. When fitting came to the given accuracy, the ultimateQ model is obtained. The result shows that an NE-trending lowQ zone exists at the depths of 10–18 km, and an NW-trending lowQ zone exists at the depths of 12–18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones
revealed by other scientists. The difference is that the lowQ zones have a wider range than the low velocity zones.
Foundation item: Joint Seismological Science Foundation of China (957-07-414) and State Key Basic Research Development and Programming Project
(95-13-02-02).
Contribution No. RCEG200105, Research Center of Exploration Geophysics, China Seismological Bureau. 相似文献
14.
The attenuation of P- and S-waves in Southeastern Sicily was estimated by applying two different methods in time and frequency domains. We analyzed waveforms from about 290 local events (0.6≤ML≤4.6) recorded at a three-component digital network.By applying the pulse broadening method to the first P-wave pulse, we found an average Qp value of ca. 140. The application of the frequency decay method provided a Qp value of ca. 120, in the low-frequency band (3-9 Hz). Conversely, in the high frequency range (16-27 Hz) the average Qp is significantly larger (ca. 640). The frequency decay method was also applied to S-waves spectra. In the low frequency range (2-5 Hz) the estimated average Qs is ca. 190. As for Qp, also Qs, in the high frequency range (16-27 Hz), is larger (ca. 700). These results evidenced a frequency dependence of both the quality factors Qp and Qs, as commonly observed in tectonically active zones characterized by high degree of heterogeneity. 相似文献
15.
Compressional-wave Q estimation from full-waveform sonic data 总被引：1，自引：0，他引：1
There is significant evidence that the anelastic loss of seismic energy is linked to petrophysical properties such as porosity, permeability and clay content. Thus, reliable estimation of anelastic attenuation from seismic data can lead to improved methods for the prediction of petrophysical properties. This paper is concerned with methods for the estimation of attenuation at sonic frequencies (5–30 KHz) from in situ data. Two independent methods have been developed and tested for estimating compressional‐wave attenuation from full‐waveform sonic data. A well‐established technique, the logarithm spectral ratio (LSR) method, is compared with a new technique, the instantaneous frequency (IF) method. The LSR method uses the whole spectrum of the seismic pulse whilst the IF method uses a carefully estimated value of instantaneous frequency which is representative of the centre frequency of the pulse. In the former case, attenuation estimation is based on the relative variation of amplitudes at different frequencies, whilst in the latter case it is based on the shift of the centre frequency of the pulse to lower values during anelastic wave propagation. The IF method does not assume frequency independence of Q which is a necessary assumption for the LSR method, and it provides a stable frequency log, the peak instantaneous frequency (PIF) log, which may be used as an indicator for attenuation under certain limitations. The development and implementation of the two methods is aimed at minimizing the effect of secondary arrivals, such as leaky modes, and involved a series of parameter tests. Testing of the two methods using full‐waveform sonic data of variable quality, obtained from a gas‐bearing sandstone reservoir, showed that the IF method is in general more stable and suitable for full‐waveform sonic data compared with the LSR method. This was evident especially in data sets with high background noise levels and wave‐interference effects. For good quality data, the two methods gave results that showed good agreement, whilst comparison with other log types further increased confidence in the results obtained. A significant decrease (approximately 5 KHz) in the PIF values was observed in the transition from an evaporite/shale sequence to the gas‐bearing sandstone. Average Q values of 54 and 51 were obtained using good quality data from a test region within the gas‐saturated sandstone reservoir, using the LSR and IF methods, respectively. 相似文献
16.
Xintao Chai Ronghua Peng Genyang Tang Wei Chen Jingnan Li 《Geophysical Prospecting》2019,67(8):2003-2021
The subsurface media are not perfectly elastic, thus anelastic absorption, attenuation and dispersion (aka Q filtering) effects occur during wave propagation, diminishing seismic resolution. Compensating for anelastic effects is imperative for resolution enhancement. Q values are required for most of conventional Q-compensation methods, and the source wavelet is additionally required for some of them. Based on the previous work of non-stationary sparse reflectivity inversion, we evaluate a series of methods for Q-compensation with/without knowing Q and with/without knowing wavelet. We demonstrate that if Q-compensation takes the wavelet into account, it generates better results for the severely attenuated components, benefiting from the sparsity promotion. We then evaluate a two-phase Q-compensation method in the frequency domain to eliminate Q requirement. In phase 1, the observed seismogram is disintegrated into the least number of Q-filtered wavelets chosen from a dictionary by optimizing a basis pursuit denoising problem, where the dictionary is composed of the known wavelet with different propagation times, each filtered with a range of possible values. The elements of the dictionary are weighted by the infinity norm of the corresponding column and further preconditioned to provide wavelets of different values and different propagation times equal probability to entry into the solution space. In phase 2, we derive analytic solutions for estimates of reflectivity and Q and solve an over-determined equation to obtain the final reflectivity series and Q values, where both the amplitude and phase information are utilized to estimate the Q values. The evaluated inversion-based Q estimation method handles the wave-interference effects better than conventional spectral-ratio-based methods. For Q-compensation, we investigate why sparsity promoting does matter. Numerical and field data experiments indicate the feasibility of the evaluated method of Q-compensation without knowing Q but with wavelet given. 相似文献
17.
Nicolas Florsch Hilaire Legros Jacques Hinderer 《Physics of the Earth and Planetary Interiors》1995,90(3-4)
When considering the search for discovery or amplitude estimation of a spectral line with a probabilistic approach, great attention must be paid to the meaning of each step. We give the probability law for the amplitude of a spectral peak in the presence of random noise appearing in a periodogram and discuss the effective probability of the existence of the corresponding wave. We find that the estimated amplitude of a spectral peak is biased and should be corrected when the signal-to-noise ratio is small. As a first application to gravity data, it results in a re-estimation of the gravimetric amplitude factors (delta factors) provided by least-squares tidal analysis. We also estimate the probability of observing a spectral line above a given level in the spectrum of a purely random noise. This allows us to compute for given spectrum the number of peaks expected to overcross the classical levels used in statistical analysis (like nσ, where σ is the standard deviation of the temporal noise distribution and n is an integer with typical values equal to 2 or 3). A specific application to real data is investigating the gravity spectrum derived from a 5 year record of the French superconducting gravimeter and we show that the predicted statistics are indeed in agreement with the observations. We also show the statistical consequence of using longer observing periods to obtain the spectral estimations. The problem of detecting translational motion of the Earth's solid inner core (Slichter modes) in a gravity spectrum is analyzed and the probabilities of having a triplet of random peaks thresholding specific levels in a given frequency window are computed. We show that, in the case of a typical gravity spectrum (1 year of hourly data and a frequency window of 0.03 cycle h−1), the probability of having a random set of three peaks exceeding a level of 3 σ, is very high. This emphasizes the need for a very careful analysis of spectral lines before inferring the existence of a true physical signal. 相似文献
18.
We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective QS/QP ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small. 相似文献
19.
Marian Ivan 《Journal of Seismology》2007,11(1):73-85
The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during
the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations
are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The
amplitude spectral ratio method is used, assuming no frequency dependence of the Q
P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor
is observed. Lateral variations of the attenuation are obtained, with a very low Q
P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different
azimuth angles in relation to the epicentral area, Q
P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located
immediately beneath the crust, but limited in depth to at least 100 km. 相似文献
20.
Using simulated coda waves, the resolution of the single-scattering model to extract codaQ (Q
c
) and its power law frequency dependence was tested. The back-scattering model ofAki andChouet (1975) and the single isotropic-scattering model ofSato (1977) were examined. The results indicate that: (1) The inputQ
c
models are reasonably well approximated by the two methods; (2) almost equalQ
c
values are recovered when the techniques sample the same coda windows; (3) lowQ
c
models are well estimated in the frequency domain from the early and late part of the coda; and (4) models with highQ
c
values are more accurately extracted from late code measurements. 相似文献