首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data.  相似文献   

2.
波场延拓深度滤波方法   总被引:5,自引:1,他引:4       下载免费PDF全文
消除面波是地震数据处理中的一个重要内容.本文提出了基于15°波动方程的深度滤波方法.由于面波与有效反射波具有不同的传播深度,可利用波场向下延拓方法将二者进行波场分离.把向下延拓后的波场中集中在地表附近的面波能量切除后,再将波场重新延拓回原始的观测面,达到去除干扰的目的.实际资料处理显示:方法计算稳定,消除面波能力强,能更好地保持波场的有效成份和幅值,符合波场的实际传播状态,表明该方法正确可行.  相似文献   

3.
The broadband capabilities of marine, seabed, and land seismic equipment are reviewed with respect to both the source and the receiver sides. In marine acquisition, the main issue at both ends of the spectrum relates to ghosts occurring at the sea surface. Broadband deghosting requires towing at variable depth to introduce notch diversity or using new equipment like multi‐component and/or low‐noise streamers. As a result, a doubling of the bandwidth from about three to six octaves (2.5–200 Hz) has been achieved. Such improvement is not yet observed for seabed surveys in spite of deghosting being a standard process on the receiver side. One issue may be related to the coupling of the particle motion sensor, particularly at high frequencies. For land acquisition, progress came from the vibrators. New shakers and control electronics using broadband sweeps made it possible to add two more octaves to the low‐frequency signal (from 8 to 2 Hz). Whereas conventional 10 Hz geophones are still able to record such low frequencies, 5 Hz high gain geophones or digital accelerometers enhance them to keep the signal above the noise floor. On the high end of the bandwidth, progress is not limited by equipment specifications. Here, the issue is related to a low signal‐to‐noise ratio due to the strong absorption that occurs during signal propagation. To succeed in enlarging the bandwidth, these improved equipment and sweeps must be complemented by a denser spatial sampling of the wavefield by point–source and point–receiver acquisition.  相似文献   

4.
P-alert台网实时数据对地震预警及烈度速报和工程地震研究都是重要的补充,处理分析这些数据对客观衡量P-alert台网数据质量和数据用途有重要意义。对2016年2月6日台湾美浓ML6.4地震P-alert台网获取的记录进行了处理和初步分析,统计分析显示100gal以上加速度记录有112条,200gal以上加速度记录有33条,400gal以上加速度记录有7条,东西向最大峰值加速度为466.4gal,南北向最大峰值加速度为498.4gal,竖直向最大峰值加速度为258.6gal,最大仪器地震烈度为9.5度。竖直向峰值加速度和峰值速度比水平向峰值加速度和峰值速度衰减快。峰值加速度比峰值速度衰减快,观测峰值加速度和峰值速度与台湾西南地区峰值加速度和峰值速度衰减公式比较一致。计算得到了近场台站的永久位移,显示P-alert台网绝大多数永久位移在1cm到5cm之间,最大永久位移达8cm。  相似文献   

5.
基于连续小波变换的自适应面波压制方法   总被引:5,自引:3,他引:2       下载免费PDF全文
面波干扰的抑制是陆上地震资料处理的主要问题之一.本文根据炮集记录中面波与反射波主要能量在小波域分布区域的不同,以及面波干扰的影响随炮检距变化等特点,提出了一种具有时变、空变特性的自适应面波衰减方法.文中将该方法用于模型及实际炮集资料的处理,并与常用的高通滤波方法进行对比,结果表明,该方法在衰减面波干扰的同时,能更好地保持反射波的振幅及相位信息.  相似文献   

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

7.
In the field of seismic exploration, ground roll seriously affects the deep effective reflections from subsurface deep structures. Traditional curvelet transform cannot provide an adaptive basis function to achieve a suboptimal denoised result. In this paper, we propose a method based on empirical curvelet transform (ECT) for ground roll attenuation. Unlike the traditional curvelet transform, this method not only decomposes seismic data into multiscale and multi-directional components, but also provides an adaptive filter bank according to frequency content of seismic data itself. So, ground roll can be separated by using this method. However, as the frequency of reflection and ground roll components are close, we apply singular value decomposition (SVD) in the curvelet domain to differentiate the ground roll and reflection better. Examples of synthetic and field seismic data reveal that the proposed method based ECT performs better than the traditional curvelet method in terms of the suppression of ground roll.  相似文献   

8.
Faithful recording of the elastic wavefield at the sea‐bed is required for quantitative applications of 4C seismic. The accuracy of the recorded vectorial wavefield depends on factors that vary from deployment to deployment. This paper focuses on one such factor: the interaction of the acquisition system with the sea‐bed, which is referred to here as coupling. We show, using multi‐azimuth data recorded with a cable‐based sea‐bed acquisition system, whose sensor housing is cylindrically shaped and with the in‐line geophone fixed to the cable, that coupling depends on the propagation direction and wave type (P‐ or S‐waves) of the incident wavefield. We show that coupling is more critical for S‐waves than for P‐waves. Detection of inconsistent coupling using both P‐ and S‐waves is therefore mandatory. A data‐driven processing method to compensate for the frequency‐dependent coupling response of the cross‐line geophone is derived. Its application to field data verifies the effectiveness of the method.  相似文献   

9.
Ground roll is an interference wave that severely degrades the signal-to-noise ratio of seismic data and affects its subsequent processing and interpretation. In this study, according to differences in morphological characteristics between ground roll and reflected waves, we use morphological component analysis based on two-dimensional dictionaries to separate ground roll and reflected waves. Because ground roll is characterized by low-frequency, low-velocity, and dispersion, we select two-dimensional undecimated discrete wavelet transform as a sparse representation dictionary of ground roll. Because of a strong local correlation of the reflected wave, we select two-dimensional local discrete cosine transform as the sparse representation dictionary of reflected waves. A sparse representation model of seismic data is constructed based on a two-dimensional joint dictionary then a block coordinate relaxation algorithm is used to solve the model and decompose seismic record into reflected wave part and ground roll part.The good effects for the synthetic seismic data and application of real seismic data indicate that when using the model, strong-energy ground roll is considerably suppressed and the waveform of the reflected wave is effectively protected.  相似文献   

10.
The attenuation of seismic waves propagating in reservoirs can be obtained accurately from the data analysis of vertical seismic profile in terms of the quality-factor Q. The common methods usually use the downgoing wavefields in vertical seismic profile data. However, the downgoing wavefields consist of more than 90% energy of the spectrum of the vertical seismic profile data, making it difficult to estimate the viscoacoustic parameters accurately. Thus, a joint viscoacoustic waveform inversion of velocity and quality-factor is proposed based on the multi-objective functions and analysis of the difference between the results inverted from the separated upgoing and downgoing wavefields. A simple separating step is accomplished by the reflectivity method to obtain the individual wavefields in vertical seismic profile data, and then a joint inversion is carried out to make full use of the information of the individual wavefields and improve the convergence of viscoacoustic full-waveform inversion. The sensitivity analysis of the different wavefields to the velocity and quality-factor shows that the upgoing and downgoing wavefields contribute differently to the viscoacoustic parameters. A numerical example validates our method can improve the accuracy of viscoacoustic parameters compared with the direct inversion using full wavefield and the separate inversion using upgoing or downgoing wavefield. The application on real field data indicates our method can recover a reliable viscoacoustic model, which helps reservoir appraisal.  相似文献   

11.
In the past few years, distributed acoustic sensing has gained great interest in geophysics. This acquisition technology offers immense improvement in terms of efficiency when compared with current geophysical acquisition methods. However, the fundamentals of the measurement are still not fully understood because direct comparisons of fibre data with conventional geophysical sensors are difficult during field tests. We present downscaled laboratory experiments that enabled us to characterise the relationship between the signals recorded by conventional seismic point receivers and by distributed fibre optic sensors. Interrogation of the distributed optical fibre sensor was performed with a Michelson interferometer because this system is suited to compact test configurations, and it requires only a very simple data processing workflow for extracting the signal outputs. We show acoustic data that were recorded simultaneously by both the fibre optical interferometer and conventional three‐component accelerometers, thus enabling the comparison of sensor performance. We present results focused on the directionality of fibre measurements, on the amplitude variation with angle of incidence, and on the transfer function that allows accelerometer signals to be transformed into optical fibre signals. We conclude that the optical fibre response matches with the array of the displacement differences of the inline accelerometers deployed along the fibre length. Moreover, we also analysed the influence of various types of coupling and fibre cable coating on the signal responses, emphasising the importance of these parameters for field seismic acquisitions when using the distributed fibre optic technology.  相似文献   

12.
在合肥地震台和安徽省地震局仪器检测室均架设MI 3000和GL-P2C地震烈度仪,分别进行同台址观测,对记录的背景噪声进行RMS值和噪声功率谱计算,对比分析2种烈度仪性能。研究结果显示:①自噪声水平:2种烈度仪水平向自噪声水平均约为垂直向的50%,且MI 3000地震烈度仪约为GL-P2C地震烈度仪的3倍;②频带范围:GL-P2C地震烈度仪比MI 3000地震烈度仪频带范围更宽,且在0.1—40 Hz范围内,比MI 3000地震烈度仪记录信号的噪声功率谱降低约10 dB。分析认为,相对于MI 3000地震烈度仪,GL-P2C地震烈度仪性能更佳。  相似文献   

13.
14.
Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface.  相似文献   

15.
作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.  相似文献   

16.
利用伪谱和有限差分混合方法对兰州盆地进行二维强地面运动模拟.假设3个不同的震源深度5、10、20 km的同一地震下,以一个5层的二维剖面作为模型,以此来研究地震波传播过程并分析复杂的二维非均质路径效应,以及局部地下速度构造对强地面运动振幅加乘的影响,从而了解路径效应对强地面运动最大峰值位移的影响程度.研究结果显示:伪谱和有限差分混合方法模拟地震波场结合了有限差分法和伪谱法的优点,弥补了二者的不足,能较好地处理介质不连续面的计算,同时保证了和伪谱法相当的计算精度.通过地震波场模拟可知沉积盆地的强地面运动较基岩相比具有放大效应,震源深度的不同对兰州盆地垂直分量地面最大峰值位移影响较大,在震源深度为5 km时最小,为0.06 cm;震源深度10 km时在水平剖面30~55 km范围内峰值位移最大,达到0.14 cm;但对水平分量的峰值位移影响较小.  相似文献   

17.
Three‐dimensional seismic survey design should provide an acquisition geometry that enables imaging and amplitude‐versus‐offset applications of target reflectors with sufficient data quality under given economical and operational constraints. However, in land or shallow‐water environments, surface waves are often dominant in the seismic data. The effectiveness of surface‐wave separation or attenuation significantly affects the quality of the final result. Therefore, the need for surface‐wave attenuation imposes additional constraints on the acquisition geometry. Recently, we have proposed a method for surface‐wave attenuation that can better deal with aliased seismic data than classic methods such as slowness/velocity‐based filtering. Here, we investigate how surface‐wave attenuation affects the selection of survey parameters and the resulting data quality. To quantify the latter, we introduce a measure that represents the estimated signal‐to‐noise ratio between the desired subsurface signal and the surface waves that are deemed to be noise. In a case study, we applied surface‐wave attenuation and signal‐to‐noise ratio estimation to several data sets with different survey parameters. The spatial sampling intervals of the basic subset are the survey parameters that affect the performance of surface‐wave attenuation methods the most. Finer spatial sampling will reduce aliasing and make surface‐wave attenuation easier, resulting in better data quality until no further improvement is obtained. We observed this behaviour as a main trend that levels off at increasingly denser sampling. With our method, this trend curve lies at a considerably higher signal‐to‐noise ratio than with a classic filtering method. This means that we can obtain a much better data quality for given survey effort or the same data quality as with a conventional method at a lower cost.  相似文献   

18.
Wave field reconstruction – the estimation of a three‐dimensional (3D) wave field representing upgoing, downgoing or the combined total pressure at an arbitrary point within a marine streamer array – is enabled by simultaneous measurements of the crossline and vertical components of particle acceleration in addition to pressure in a multicomponent marine streamer. We examine a repeated sail line of North Sea data acquired by a prototype multicomponent towed‐streamer array for both wave field reconstruction fidelity (or accuracy) and reconstruction repeatability. Data from six cables, finely sampled in‐line but spaced at 75 m crossline, are reconstructed and placed on a rectangular data grid uniformly spaced at 6.25 m in‐line and crossline. Benchmarks are generated using recorded pressure data and compared with wave fields reconstructed from pressure alone, and from combinations of pressure, crossline acceleration and vertical acceleration. We find that reconstruction using pressure and both crossline and vertical acceleration has excellent fidelity, recapturing highly aliased diffractions that are lost by interpolation of pressure‐only data. We model wave field reconstruction error as a linear function of distance from the nearest physical sensor and find, for this data set with some mismatched shot positions, that the reconstructed wave field error sensitivity to sensor mispositioning is one‐third that of the recorded wave field sensitivity. Multicomponent reconstruction is also more repeatable, outperforming single‐component reconstruction in which wave field mismatch correlates with geometry mismatch. We find that adequate repeatability may mask poor reconstruction fidelity and that aliased reconstructions will repeat if the survey geometry repeats. Although the multicomponent 3D data have only 500 m in‐line aperture, limiting the attenuation of non‐repeating multiples, the level of repeatability achieved is extremely encouraging compared to full‐aperture, pressure‐only, time‐lapse data sets at an equivalent stage of processing.  相似文献   

19.
We modified the common-offset–common-reflection-surface (COCRS) method to attenuate ground roll, the coherent noise typically generated by a low-velocity, low-frequency, and high-amplitude Rayleigh wave. The COCRS operator is based on hyperbolas, thus it fits events with hyperbolic traveltimes such as reflection events in prestack data. Conversely, ground roll is linear in the common-midpoint (CMP) and common-shot gathers and can be distinguished and attenuated by the COCRS operator. Thus, we search for the dip and curvature of the reflections in the common-shot gathers prior to the common-offset section. Because it is desirable to minimize the damage to the reflection amplitudes, we only stack the multicoverage data in the ground-roll areas. Searching the CS gathers before the CO section is another modification of the conventional COCRS stacking. We tested the proposed method using synthetic and real data sets from western Iran. The results of the ground-roll attenuation with the proposed method were compared with results of the f–k filtering and conventional COCRS stacking after f–k filtering. The results show that the proposed method attenuates the aliased and nonaliased ground roll better than the f–k filtering and conventional CRS stacking. However, the computation time was higher than other common methods such as f–k filtering.  相似文献   

20.
In previous publications, we presented a waveform-inversion algorithm for attenuation analysis in heterogeneous anisotropic media. However, waveform inversion requires an accurate estimate of the source wavelet, which is often difficult to obtain from field data. To address this problem, here we adopt a source-independent waveform-inversion algorithm that obviates the need for joint estimation of the source signal and attenuation coefficients. The key operations in that algorithm are the convolutions (1) of the observed wavefield with a reference trace from the modelled data and (2) of the modelled wavefield with a reference trace from the observed data. The influence of the source signature on attenuation estimation is mitigated by defining the objective function as the ℓ2-norm of the difference between the two convolved data sets. The inversion gradients for the medium parameters are similar to those for conventional waveform-inversion techniques, with the exception of the adjoint sources computed by convolution and cross-correlation operations. To make the source-independent inversion methodology more stable in the presence of velocity errors, we combine it with the local-similarity technique. The proposed algorithm is validated using transmission tests for a homogeneous transversely isotropic model with a vertical symmetry axis that contains a Gaussian anomaly in the shear-wave vertical attenuation coefficient. Then the method is applied to the inversion of reflection data for a modified transversely isotropic model from Hess. It should be noted that due to the increased nonlinearity of the inverse problem, the source-independent algorithm requires a more accurate initial model to obtain inversion results comparable to those produced by conventional waveform inversion with the actual wavelet.  相似文献   

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

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