共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
Ground roll attenuation using the S and x-f-k transforms 总被引：2，自引：0，他引：2
Ground roll, which is characterized by low frequency and high amplitude, is an old seismic data processing problem in land‐based seismic acquisition. Common techniques for ground roll attenuation are frequency filtering, f‐k or velocity filtering and a type of f‐k filtering based on the time‐offset windowed Fourier transform. These techniques assume that the seismic signal is stationary. In this study we utilized the S, x‐f‐k and t‐f‐k transforms as alternative methods to the Fourier transform. The S transform is a type of time‐frequency transform that provides frequency‐dependent resolution while maintaining a direct relationship with the Fourier spectrum. Application of a filter based on the S transform to land seismic shot records attenuates ground roll in a time‐frequency domain. The t‐f‐k and x‐f‐k transforms are approaches to localize the apparent velocity panel of a seismic record in time and offset domains, respectively. These transforms provide a convenient way to define offset or time‐varying reject zones on the separate f‐k panel at different offsets or times. 相似文献
3.
In the present study, a digital waveform dataset of 216 local earthquakes recorded by the Egyptian National Seismic Network (ENSN) was used to estimate the attenuation of seismic wave energy in the greater Cairo region. The quality factor and the frequency dependence for Coda waves and S-waves were estimated and clarified. The Coda waves (Q c) and S-waves (Q d) quality factor were estimated by applying the single scattering model and Coda Normalization method, respectively, to bandpass-filtered seismograms of frequency bands centering at 1.5, 3, 6, 12, 18 and 24?Hz. Lapse time dependence was also studied for the area, with the Coda waves analyzed through four lapse time windows (10, 20, 30 and 40?s). The average quality factor as function of frequency is found to be Q c?=?35?±?9f 0.9±0.02 and Q d?=?10?±?2f 0.9±0.02 for Coda and S-waves, respectively. This behavior is usually correlated with the degree of tectonic complexity and the presence of heterogeneities at several scales. The variation of Q c with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth. In fact, by using the Coda Normalization method we obtained low Q d values as expected for a heterogeneous and active zone. The intrinsic quality factor (Q i ?1 ) was separated from the scattering quality factor (Q s ?1 ) by applying the Multiple Lapse Time Domain Window Analysis (MLTWA) method under the assumption of multiple isotropic scattering with uniform distribution of scatters. The obtained results suggest that the contribution of the intrinsic attenuation (Q i ?1 ) prevails on the scattering attenuation (Q s ?1 ) at frequencies higher than 3?Hz. 相似文献
4.
Giuliana Rossi Davide Gei Gualtiero Böhm Gianni Madrussani José M. Carcione 《Geophysical Prospecting》2007,55(5):655-669
We estimate the quality factor (Q) from seismic reflections by using a tomographic inversion algorithm based on the frequency‐shift method. The algorithm is verified with a synthetic case and is applied to offshore data, acquired at western Svalbard, to detect the presence of bottom‐simulating reflectors (BSR) and gas hydrates. An array of 20 ocean‐bottom seismographs has been used. The combined use of traveltime and attenuation tomography provides a 3D velocity–Q cube, which can be used to map the spatial distribution of the gas‐hydrate concentration and free‐gas saturation. In general, high P‐wave velocity and quality factor indicate the presence of solid hydrates and low P‐wave velocity and quality factor correspond to free‐gas bearing sediments. The Q‐values vary between 200 and 25, with higher values (150–200) above the BSR and lower values below the BSR (25–40). These results seem to confirm that hydrates cement the grains, and attenuation decreases with increasing hydrate concentration. 相似文献
5.
A method for classifying pre‐stack seismic data based on amplitude–frequency attributes and self‐organizing maps 
下载免费PDF全文
下载免费PDF全文 Erik Molino‐Minero‐Re Ernesto Rubio‐Acosta Héctor Benítez‐Pérez Juan Marcos Brandi‐Purata Nora Isabel Pérez‐Quezadas Demetrio Fabián García‐Nocetti 《Geophysical Prospecting》2018,66(4):673-687
6.
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. 相似文献
7.
A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks 总被引：4，自引：0，他引：4
The laboratory ultrasonic pulse‐echo method was used to collect accurate P‐ and S‐wave velocity (±0.3%) and attenuation (±10%) data at differential pressures of 5–50 MPa on water‐saturated core samples of sandstone, limestone and siltstone that were cut parallel and perpendicular to the vertical borehole axis. The results, when expressed in terms of the P‐ and S‐wave velocity and attenuation anisotropy parameters for weakly transversely isotropic media (ɛ, γ, ɛQ, γQ) show complex variations with pressure and lithology. In general, attenuation anisotropy is stronger and more sensitive to pressure changes than velocity anisotropy, regardless of lithology. Anisotropy is greatest (over 20% for velocity, over 70% for attenuation) in rocks with visible clay/organic matter laminations in hand specimens. Pressure sensitivities are attributed to the opening of microcracks with decreasing pressure. Changes in magnitude of velocity and attenuation anisotropy with effective pressure show similar trends, although they can show different signs (positive or negative values of ɛ, ɛQ, γ, γQ). We conclude that attenuation anisotropy in particular could prove useful to seismic monitoring of reservoir pressure changes if frequency‐dependent effects can be quantified and modelled. 相似文献
8.
Attenuation of seismic waves, quantified by the seismic quality factor Q, holds important information for seismic interpretation, due to its sensitivity to rock and fluid properties. A recently published study of Q, based on surface seismic reflection data, used a modified spectral ratio approach (QVO), but both source and receiver responses were treated as isotropic, based on simple raypath arguments. Here, this assumption has been tested by computing apparent attenuation generated by frequency-dependent directivity of typical marine source and receiver arrays and acquisition geometries. Synthetic wavelet spectra were computed for reflected rays, summed over the first Fresnel zone, from the base of a single interval, 50–3000 m thick and velocity 2000 m/s, overlying a 2200 m/s half-space, and for offsets of 71–2071 m. The source and receiver geometry were those of an actual survey. The modelled spectra are clearly affected by directivity, most strongly because of surface ghosts. In general, the strong high-frequency component, produced by the array design, leads to apparently negative attenuation in individual reflection events, though this is dependent on offset and target depth. For shallow targets (less than 400–500 ms two-way traveltime (TWT) depth), apparent Q-values as extreme as ?50 to ?100 were obtained. For deeper target depths, the directivity effect is far smaller. The implications of the model study were tested on real data. QVO was applied to 20 true-spectrum-processed CMPs, in a shallow (405–730 ms TWT) and a deeper (1000–1300 ms TWT) interval, firstly using a measured far-field source signature (effectively isotropic), and secondly using computed directivity effects instead. Mean interval Q?1-values for the deeper interval, 0.029 ± 0.011 and 0.027 ± 0.018 for conventional and directional processing, respectively, suggested no directivity influence on attenuation estimation. For the shallow interval (despite poor spectral signal-to-noise ratios and hence scattered attenuation estimates), directional processing removed directivity-generated irregularities from the spectral ratios, resulting in an improvement from Q?1int = ?0.036 ± 0.130 to a realistic Q?1int = 0.012 ± 0.030: different at 94% confidence level. Equivalent Q-values are: for the deeper interval, 35 and 37 for conventional and directional processing, respectively, and ?28 and 86 for the shallow interval. These results support the conclusions of the model studies, i.e. that source/receiver directivity has a negligible effect except for shallow targets (e.g. TWT depth ≤ 500 ms) imaged with conventional acquisition geometry. In such cases directivity corrections to spectra are strongly recommended. 相似文献
9.
Maksim Bano 《Geophysical Prospecting》2004,52(1):11-26
The attenuation of ground‐penetrating radar (GPR) energy in the subsurface decreases and shifts the amplitude spectrum of the radar pulse to lower frequencies (absorption) with increasing traveltime and causes also a distortion of wavelet phase (dispersion). The attenuation is often expressed by the quality factor Q. For GPR studies, Q can be estimated from the ratio of the real part to the imaginary part of the dielectric permittivity. We consider a complex power function of frequency for the dielectric permittivity, and show that this dielectric response corresponds to a frequency‐independent‐Q or simply a constant‐Q model. The phase velocity (dispersion relationship) and the absorption coefficient of electromagnetic waves also obey a frequency power law. This approach is easy to use in the frequency domain and the wave propagation can be described by two parameters only, for example Q and the phase velocity at an arbitrary reference frequency. This simplicity makes it practical for any inversion technique. Furthermore, by using the Hilbert transform relating the velocity and the absorption coefficient (which obeys a frequency power law), we find the same dispersion relationship for the phase velocity. Both approaches are valid for a constant value of Q over a restricted frequency‐bandwidth, and are applicable in a material that is assumed to have no instantaneous dielectric response. Many GPR profiles acquired in a dry aeolian environment have shown a strong reflectivity inside dunes. Changes in water content are believed to be the origin of this reflectivity. We model the radar reflections from the bottom of a dry aeolian dune using the 1D wavelet modelling method. We discuss the choice of the reference wavelet in this modelling approach. A trial‐and‐error match of modelled and observed data was performed to estimate the optimum set of parameters characterizing the materials composing the site. Additionally, by combining the complex refractive index method (CRIM) and/or Topp equations for the bulk permittivity (dielectric constant) of moist sandy soils with a frequency power law for the dielectric response, we introduce them into the expression for the reflection coefficient. Using this method, we can estimate the water content and explain its effect on the reflection coefficient and on wavelet modelling. 相似文献
10.
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. 相似文献
11.
《Journal of Applied Geophysics》2010,72(4):157-161
An intimate mathematical relation between Hartley and Hilbert transforms is given here in contrast with the well known Fourier and Hilbert transform relations. It is interesting to note that the Fourier–Hilbert and Hartley–Hilbert transforms while possessing the same magnitude differ in phase by 270°. The inverse Hartley–Hilbert transform returns the original function unlike the Fourier–Hilbert transform which results the negative of the original function. Further, it may be realized that the envelope defined here of the analytic signal in both Fourier–Hilbert and Hartley–Hilbert domains numerically remain the same while differing in polarity. The feasibility of Hartley–Hilbert transform for a straight forward interpretation, total magnetic anomaly due to a thin plate from Tejpur, India and self potential data of the Sulleymonkey anomaly in the Ergani Copper district, Turkey are illustrated in contrast with the Fourier–Hilbert transform. This pair of transforms have innumerable geophysical applications. 相似文献
12.
Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real world, seismic events usually have curvature or amplitude variability, which means that their Fourier transforms are no longer strictly linear but rather occupy conic regions of the Fourier domain that are narrow at low frequencies but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised “signal cones”. In this work, we consider a space–time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting, and translation of the kernel for cone‐limited (C‐limited) functions (functions whose Fourier transform lives within a cone) or C‐Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency), which constitutes our dictionary. We find a discrete number of scaling, slanting, and translation parameters from a continuum by optimally matching the data. This is a non‐linear optimisation problem, which we address by a fixed‐point method that utilises a variable projection method with ?1 constraints on the linear parameters and bound constraints on the non‐linear parameters. We observe that slow decay and oscillatory behaviour of the kernel for C‐limited functions constitute bottlenecks for the optimisation problem, which we partially overcome by the C‐Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity‐promoting curvelet decomposition, matching pursuit Fourier interpolation, and sparsity‐promoting plane‐wave decomposition methods. 相似文献
13.
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. 相似文献
14.
Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。 相似文献
15.
We study path effects on prediction equations of pseudo‐velocity response spectra (natural period of 0.1–5.0 s) in northern Japan, where heterogeneous attenuation structure exists. The path effects have been examined by comparing the regression analysis results for two different prediction equations. The first equation consists of a single term of anelastic attenuation conventionally. The second equation consists of two terms of anelastic attenuation in consideration of the heterogeneous attenuation structure. In the second equation, we divide a source‐to‐site distance into two distances at the attenuation boundary beneath the volcanic front. The boundary is considered to separate the relatively high Q fore‐arc side mantle wedge (FAMW) from the low Q back‐arc side mantle wedge (BAMW). Strong motion records (hypocentral distances less than 300 km) from interplate and intraslab events with Mw 5.1–7.3 are used. Regression analysis results show that the standard errors are significantly reduced by the second prediction equation at short periods (0.1–0.5 s), whereas the difference in standard errors from both prediction equations is negligible at intermediate and long periods. The Qs values (quality factor for S‐wave) converted from two anelastic attenuation coefficients for the second prediction equation are remarkably similar to the path‐averaged Qs values for the FAMW and BAMW by other studies using spectral inversion method. From these findings, we conclude that the path effects on the prediction equation of pseudo‐velocity response spectra are satisfactorily accomplished by the second prediction equation. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
16.
宽频带地震观测数据中有效信号和干扰噪声经常发生混频效应,常规的频率域滤波方法很难将二者分离.地震波信号属于时变非平稳信号,时频分析方法能够同时得到地震波信号随着时间和频率变化的振幅和相位特征,S变换是其中较为高效的时频分析工具之一.本文以S变换为例,提出了基于相位叠加的时频域相位滤波方法.与传统叠加方法相比,相位叠加方法对强振幅不敏感,对波形一致性相当敏感,更加利于有效弱信号信息的检测.时频域相位滤波方法滤除与有效信号不相干的背景噪声,保留了相位一致的有效信号成分,显著提高了信噪比.运用理论合成的远震接收函数数据和实际的宽频带地震观测数据检验结果显示该方法较传统的带通滤波方法相比,即使在信噪较低且混频严重条件下,时频域相位滤波方法的滤波效果依然很明显,有助于识别能量较弱的有效信号. 相似文献
17.
Z.A. Der W.D. Rivers T.W. McElfresh A. O&#x;Donnell P.J. Klouda M.E. Marshall 《Physics of the Earth and Planetary Interiors》1982,30(1):12-25
A worldwide study of short-period teleseismic body wave spectra shows that the high frequency falloff rates of spectra are correlated with the tectonic type of the source and receiver regions and with source depth. The data indicate, in a consistent manner, that the main cause for such variations is the lateral variation of Q in the upper mantle as well as change of Q with depth. Using the internal consistency checks provided by redundancies in the data set other effects such as crustal, site dependent distortion of the spectra, source effects and instrument non-linearity can be ruled out as significant factors influencing the estimates obtained. The results indicate high attenuation in the upper mantle under tectonic regions and new oceans. Long-period regional attenuation studies indicate similar variations in mantle Q among the types of regions mentioned but yield significantly lower Q estimates in all areas. The short- and long-period attenuation results can be reconciled only by assuming a frequency dependent Q that increases with frequency along all types of paths, such that the relative differences in Q along various types of paths retain the same sign over the short- and long-period bands. 相似文献
18.
Sediment rating curves are commonly used to estimate the suspended sediment load in rivers and streams under the assumption of a constant relation between discharge (Q) and suspended sediment concentrations (SSC) over time. However, temporal variation in the sediment supply of a watershed results in shifts in this relation by increasing variability and by introducing nonlinearities in the form of hysteresis or a path‐dependent relation. In this study, we used a mixed‐effects linear model to estimate an average SSC–Q relation for different periods of time within the hydrologic cycle while accounting for seasonality and hysteresis. We tested the performance of the mixed‐effects model against the standard rating curve, represented by a generalized least squares regression, by comparing observed and predicted sediment loads for a test case on the Chilliwack River, British Columbia, Canada. In our analyses, the mixed‐effects model reflected more accurate patterns of interpolated SSC from Q data than the rating curve, especially for the low‐flow summer months when the SSC–Q relation is less clear. Akaike information criterion scores were lower for the mixed‐effects model than for the standard model, and the mixed‐effects model explained nearly twice as much variance as the standard model (52% vs 27%). The improved performance was achieved by accounting for variability in the SSC–Q relation within each month and across years for the same month using fixed and random effects, respectively, a characteristic disregarded in the sediment rating curve. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
19.
To better understand (and correct for) the factors affecting the estimation of attenuation (Q), we simulate subsurface wave propagation with the Weyl/Sommerfeld integral. The complete spherical wavefield emanating from a P‐wave point source surrounded by a homogeneous, isotropic and attenuative medium is thus computed. In a resulting synthetic vertical seismic profile, we observe near‐field and far‐field responses and a 90° phase rotation between them. Depth dependence of the magnitude spectra in these two depth regions is distinctly different. The logarithm of the magnitude spectra shows a linear dependence on frequency in the far‐field but not in those depth regions where the near‐field becomes significant. Near‐field effects are one possible explanation for large positive and even negative Q‐factors in the shallow section that may be estimated from real vertical seismic profile data when applying the spectral ratio method. We outline a near‐field compensation technique that can reduce errors in the resultant Q estimates. 相似文献
20.
G.‐Akis Tselentis Paraskevas Paraskevopoulos Nikos Martakis 《Geophysical Prospecting》2010,58(5):845-859
We determine the attenuation structure of a three‐dimensional medium based on first pulse‐width measurements from microearthquake data. Ninety‐five microearthquakes from a seventy stations local network were considered in this study. Measurements of the first half cycle of the wave, the so‐called rise time τ were carried out on high quality velocity seismograms and inverted to estimate the P‐waves intrinsic quality factor Qp. The results of this investigation indicate that first pulse width data from a local microearthquake network permit retrieval with sufficient accuracy of the heterogeneous Qp structure. The inferred attenuation variability corresponds to the known geological formations in the region. 相似文献