共查询到20条相似文献，搜索用时 31 毫秒
1.
不连续地质体（如断层）的自动检测一直以来都是叠后地震数据解释中的关键问题之一,尤其在三维情况中尤为重要.然而,大多数边缘检测和相干算法都对随机噪声很敏感,随机噪声衰减是叠后地震数据解释的另一个主要问题.针对构造保护去噪和断层检测问题,本文基于非平稳相似性系数完善一种构造导向滤波方法并且提出一种自动断层检测方法,形成了一套匹配的处理技术.该构造导向滤波既能够有效地衰减随机噪声又可以很好地保护地震资料中的断层等信息不被破坏,增强地震剖面中弯曲、倾斜同相轴的连续性.根据地震数据局部倾角走向,利用相邻道构建当前地震道的预测,通过预测道的叠加得到参考道,计算预测道与参考道之间的非平稳相似性系数可以设计出数据驱动的加权中值滤波.另一方面,预测道与原始道之间的非平稳相似性系数能够用于带有断层指示性的相干分析.这两种方法都基于构造预测和非平稳相似性系数,但是使用不同的调节参数和处理方案.理论模型和实际数据的处理结果证明了本文提出构造导向滤波和断层检测方法的有效性. 相似文献
2.
随着地震勘探精度的不断提高,常规速度分析方法已不能满足精细构造分析、储层描述及各向异性研究等对速度精度的要求.研究应用高密度叠前速度求取技术,自动拾取地震记录规则波的走时、振幅、波场相干性及炮点射线参数、检波点射线参数等信息,求取不同偏移距道上波峰时间处的反射波速度,并针对地质条件,采用空变时变速度滤波、频率滤波、倾角参数滤波等手段提高速度精度.将该技术应用于塔河油田西南部TFT地区1200km2叠前三维地震资料的精细速度分析,落实了主要目的层低幅度构造,确定了构造形态和要素;通过差异速度分析、AVO反射能量包络差及速度变化规律研究,对泥盆系东河砂岩进行了储层异常区带划分. 相似文献
3.
Traditional coherence algorithms are often based on the assumption that seismic traces are stationary and Gaussian. However, seismic traces are actually non-stationary and non-Gaussian. A constant time window and the canonical correlation analysis in traditional coherence algorithms are not optimal for non-stationary seismic traces and cannot describe the similarity between adjacent seismic traces in detail. To overcome this problem, a new coherence algorithm using the high-resolution time–time transform and the feature matrix is designed. The high-resolution time–time transform used to replace the constant time window can produce a frequency-dependent time local series to analyse non-stationary seismic traces. The feature matrix, constructed by the frequency-dependent time local series and the related local gradients, defines a new correlation metric that enhances more details of the geological discontinuities in seismic images than does the canonical correlation analysis. Additionally, the Riemannian metric is introduced for related calculations because the feature matrices are not defined in a Euclidean space but rather in a manifold space. Application to field data illustrates that the proposed method reveals more details of structural and stratigraphic features. 相似文献
4.
Several papers presented at the last SEG Convention in Houston by Schneider, Backus et al have shown how important and fruitful it was to obtain a continuous knowledge of the velocity functions and they have solved their problem by a Dynamic Correlation Analysis. Our purpose is to introduce here a method based on the best summation of a set of traces instead of the best correlation. Practically, this approach has several advantages: 1) Two traces only can be correlated at each step whereas the summation can bear on any number of them; 2) Optimizing the summation is actually what we are looking for since, at the long end, the success of the improvement is evaluated from the compositing of several traces either weighted or not. On the other hand, an advantage of correlation is the possibility of adding correlations obtained at several places in a same neighbourhood in order to improve the results. With the summation method this is feasible only when dips are inexistent: we shall see that the difficulty due to the dip effect can be turned around. The basic principle of the method can be summed up as follows: traces relating to a same reflection point are considered; several composites are made, each after applying different move out corrections ranging widely around an estimated adequate velocity function. At each time coordinate, the best adapted velocity function, i.e. the one that yields the best phase relation between reflected events, corresponds to the composite trace the average amplitude of which is the largest. This way, the velocity function corresponding to primary reflections as well as those corresponding to multiple reflections can be established accurately. Some examples are shown. 相似文献
5.
In this work, we propose the application of the wavelet transform analysis in well-logs (radioactivity, resistivity and sonic) to identify facies. The wavelet transform is applied to a set of well-log data for identifying correlations between wavelet coefficients and lithofacies sequences. Our results indicate that the scales, in a multiscale analysis, are related to the rock thickness and depending on the scale used it is possible to identify other particular or general sequences. The results obtained are compared and corroborated by standard geological procedures for lithological characterization, indicating that the wavelet analysis provides qualitative guides for the identification of lithological properties in wells. All our analyses are based on a siliciclastic oil field that belongs to Chicontepec Formation of the Tampico–Misantla basin in Mexico. 相似文献
6.
Implication of seismic noise for determining the structure of the upper Earth Rock Mass 总被引：1，自引：0，他引：1
T. Yu. Koroleva T. B. Yanovskaya S. S. Patrusheva 《Izvestiya Physics of the Solid Earth》2009,45(5):369-380
An assertion that the cross-correlation function of seismic noise, considered as a result of the superposition of the surface waves, excited by the sources, randomly distributed over the Earth’s surface, determines the Green function of the surface wave is verified by numerical modeling. The maximum wave periods, for which this assertion is correct, are estimated and the errors in determination of the phase and group velocity of the surface waves are evaluated. The procedure for the determination of the correlation function and estimation from it of the group velocity are tested thoroughly based on the example of the pair of the BJT and TLY stations in Asia. This procedure is used for obtaining the group and phase velocities of the Rayleigh waves on the traces between the OBN-ARU and PUL-ARU stations. The velocity sections of the transverse waves are built based on the dispersion curves of the phase and group velocities of averages along these traces. The region of the lowered velocity in the upper mantle at depths of 150–300 km is revealed on both traces. From the analysis of correlation functions, which are subjected to narrow-band filtering, it is shown that the frequency composition of noise varies from the East and from the West from the profiles between the stations: in the East (Siberia) the noise has an appreciably lower-frequency than in the West (Western Europe). 相似文献
7.
Spectral decomposition is a widely used technique in analysis and interpretation of seismic data. According to the uncertainty principle, there exists a lower bound for the joint time–frequency resolution of seismic signals. The highest temporal resolution is achieved by a matching pursuit approach which uses waveforms from a dictionary of functions (atoms). This method, in its pure mathematical form can result in atoms whose shape and phase have no relation to the seismic trace. The high‐definition frequency decomposition algorithm presented in this paper interleaves iterations of atom matching and optimization. It divides the seismic trace into independent sections delineated by envelope troughs, and simultaneously matches atoms to all peaks. Co‐optimization of overlapping atoms ensures that the effects of interference between them are minimized. Finally, a second atom matching and optimization phase is performed in order to minimize the difference between the original and the reconstructed trace. The fully reconstructed traces can be used as inputs for a frequency‐based reconstruction and red–green–blue colour blending. Comparison with the results of the original matching pursuit frequency decomposition illustrates that high‐definition frequency decomposition based colour blends provide a very high temporal resolution, even in the low‐energy parts of the seismic data, enabling a precise analysis of geometrical variations of geological features. 相似文献
8.
提出一种自适应协方差的时频域极化滤波方法。该方法在广义S变换时频方法的基础上,构造时频域自适应协方差矩阵,通过特征分析计算时频域瞬时极化参数,设计极化滤波器,实现多分量地震极化分析和滤波。其优势在于协方差矩阵的分析时窗的长度由多分量地震数据的瞬时频率确定,可以自适应于有效信号的周期,在每个时频点计算极化参数不需要进行插值处理;结合时间频率信息,解决在时间域或频率域波形或频率重叠的信号具有明显的直观性。模型数据及实际三分量台站地震数据处理结果表明,该极化滤波方法在台站地震资料分析和处理方面具有很好的直观性和较高的分辨率。 相似文献
9.
The widespread use of common depth point techniques has emphasized the need for accurate static corrections. Manual interpretation methods can give excellent results, but a computer technique is desirable because of the great volumn of data recorded in common depth point shooting. The redundancy inherent in common depth point data may be used to compute a statistical estimate of the static corrections. The corrections are assumed to be time-invarient, surface-consistent, and independent of frequency. Surface consistency implies that all traces from a particular shot will receive the same shot static correction and all traces from a particular receiver position will receive the same receiver correction. Time shifts are computed for all input traces using crosscorrelation functions between common depth point traces. The time shift for each trace is composed of a shot static, a receiver static, residual normal moveout if present, and noise. Estimates of the shot and receiver static corrections are obtained by averaging different sets of the measured time shifts. Time shifts which are greatly in error are detected and removed from the computations. The method is useful for data which has a moderate to good signal to noise ratio. Residual normal moveout should be corrected before estimating the statics. The program estimates the statics for correctly stacking common depth point traces but it is not sensitive to constant or very slowly changing static errors. 相似文献
10.
Proper stacking of three-dimensional seismic CDP-data generally requires the knowledge of normal moveout velocities in all source-receiver directions contributing to a CDP-gather. The azimuthal variation of the stacking velocities mainly depends on the dip of the seismic interfaces. For a single dipping plane a simple relation exists between the dip and the azimuthal variation of NMO-velocity. Varying strike and dip of subsequent reflectors, however, result in a complex dependency of the seismic parameters. Reliable information on the spatial distribution of the normal moveout (NMO)-velocity can be derived from a wavefront curvature estimation using a 3-D ray-tracing technique. These procedures require additional information, e.g. reflection time gradients or depth maps to show interval velocities between leading interfaces. Moreover, their application to an extended 3-D data volume is restricted by high costs. The need for a routine 3-D procedure resulted in a special data selection to create pseudo 2-D profiles and to apply existing velocity estimation routines to these profiles. At least three estimates in different directions are necessary to derive the full azimuthal velocity variation, characterized by the large and the small main axis and the orientation of the velocity ellipse. Errors are estimated by means of computer models. Stacking velocities obtained by mathematical routines (least-squares fit) and by seismic standard routines (NMO-correction and correlation) are compared. Finally, a general 3-D velocity procedure using cross-correlation of preliminarily NMO-corrected traces is proposed. 相似文献
11.
Harmonic or subharmonic noise is often present in vibroseis data as reverberation‐like, laterally coherent bands occurring parallel to and before or after, the main events. Such periodic noise is typically generated during the standard correlation process when the actual source signal travelling through the subsurface is, for whatever reason, different from the desired source signal, i.e., the pilot‐sweep controlling the baseplate and used for correlation. A typical cause can be that harmonic or subharmonic frequency partials are generated in addition to the vibroseis sweep's desired fundamental frequencies. These harmonics produce strong ‘ghost events’ during correlation of the geophone trace with the pilot‐sweep, originating from additional correlations between the fundamental and harmonic frequencies. Especially subharmonic ‘ghosts’ will overlap with ‘good’ fundamental signals, since for typically used up‐sweeps they are folded to later traveltimes, where the signal/noise‐ratio is already lower, thus aggravating or preventing a reliable interpretation of possible later reflections. Here, a method is introduced to remove these unwanted noise trains (with only negligible impact on the fundamental signal) by transforming the seismogram traces into a so‐called ‘(sub)harmonic domain’. In this domain, the respective harmonic noise portions are focused and separated from the fundamental signals, enabling easier detection and appropriate suppression. After back‐transformation to the x‐T domain, the records are free from the corresponding harmonic contamination and can then be processed as usual. The method operates in a data‐driven fashion, i.e., the traces are not uniformly processed but are processed depending upon their actual (sub)harmonic content. The decontamination procedure can be applied universally, i.e., to uncorrelated/correlated and/or vertically unstacked/stacked data either in a manual, semiautomated or fully automated manner. The method works perfectly for synthetic vibroseis traces with or without harmonic/subharmonic portions. The application to real, crustal‐scale vibroseis records that were acquired in 2006 in the Dead Sea region, Israel and that were severely contaminated by subharmonic ground‐roll ghosts covering reflectivity from the basement to the Moho, shows the robustness and success of the presented method. 相似文献
12.
对稀疏/非规则采样或者低信噪比数据,射线束提取困难并伴随有假频产生,对叠加剖面和道集造成严重干扰.为了提升射线束偏移在稀疏和低信噪比地震数据采集中的成像效果,本文提出基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法.射线束偏移首先将地震数据划分为超道集,经过部分NMO后转化为以射线束中心定义的共偏移距数据,倾斜叠加和反假频操作均在局部共中心点坐标上实现.时间域倾斜叠加是对地震数据的时移累加操作,三角低通滤波同样可以在时间域完成,在对地震数据进行因果和反因果积分后,亦为地震数据的时移累加.因此,三角低通滤波与倾斜叠加可在时间域结合同时完成,避免了频域滤波的正反傅里叶变换.本文在反假频公式中加入权重系数,用以对反假频的程度进行控制,达到分辨率和噪声压制的最佳折衷.以某海上三维实际数据为例,文中展示了反假频射线束形成对偏移叠加剖面和共成像点偏移距道集中的噪声进行了有效压制. 相似文献
13.
首先讨论倾角测井资料自适应窗长人机交互式对比方法．依据地质背景资料对倾角测井电导率曲线进行粗略分层，然后电脑进行精细对比，最后人的知识和经验对对比结果进行确认与解释，从而保证每个资料点对比结果的可靠与客观性．同时讨论了智能型知识库支持下的倾角测井资料交互式构造解释方法，采用人机交互技术识别地层倾角矢量模式组合。然后完成交互式地质解释过程．实践表明，上述方法在复杂地下地质条件和井眼条件下是一种有效而快速的方法，同时可以克服倾角测井资料地质响应的多解性． 相似文献
14.
V. BARDAN 《Geophysical Prospecting》1987,35(4):343-358
Trace interpolation of spatially aliased seismic data is used to improve the quality of multi-trace processing, especially that of migration. We consider non-rectangular sampling lattices. After stating the real sampling requirements for (two-dimensional or 3-D) digitized signals, the paper explains geometrically the effect of migration and of the 2-D dip filtering upon the aliased events, as well as our ability to interpolate traces from undersampled seismic data. 相似文献
15.
Nonlinear sweeps have often successfully been employed in the 1960s. However, this area of sweep technology has been neglected since the introduction of digital recording techniques in the Vibroseis system. Now the advent of computerized recording instruments yields a new economical possibility of forming approximately nonlinear sweeps by combining several linear sweeps with or without time gaps to a “Combisweep”. The total duration of a Combisweep may be as long as the maximum available recording time, for example 32 s. Beside the attenuation of correlation noise, the new method has further merits, such as the weighting of predetermined frequency ranges, in order to effect a certain kind of optimum filtering on the emitter side, or in order to compensate to some degree for frequency dependent absorption. In all these applications the Combisweep is considered as one signal in the correlation process. But by correlating with the individual sweeps or a partial combination of them and by applying automatic switching at predetermined times within the gaps between the individual sweeps additional possibilities arise, such as obtaining in one run with a twenty-four channel recording unit twenty-four traces with small distances between vibrators and geophones for shallow reflections and another twenty-four traces with larger distances for deeper reflections. Various Combisweeps and their applications are presented. 相似文献
16.
D. WALTHAM 《Geophysical Prospecting》1992,40(3):343-358
A method is presented for developing and/or evaluating 2D filters applied to seismic data. The approach used is to express linear 2D filtering operations in the space-frequency (x–ω) domain. Correction filters are then determined using plane-wave constraints. For example, requiring a vertically propagating plane wave to be unaffected by migration necessitates application of a half-derivative correction in Kirchhoff migration. The same approach allows determination of the region of time-offset space where half-derivative corrections are correct in x–t domain dip moveout. Finally, an x–ω domain dip filter is derived using the constraint that a plane wave be attenuated as its dip increases. This filter has the advantage that it is significantly faster than f–k domain dip filtering and can be used on irregularly spaced data. This latter property also allows the filter to be used for interpolation of irregular data onto a regular grid. 相似文献
17.
Phase identification and attribute analysis of broadband seismograms at far-regional distances 总被引：1，自引：0，他引：1
A trained analyst can frequently provide a rapid assessment of a seismic record and provide identification for many seismic phases. For digital data a challenge is to find methods (or combinations of methods) which can provide equivalent levels of phase identification and attribute analysis. Until now, there have been few discussions on phase attribute analysis for broadband records, even though the character of the major phases has been recognised several decades ago. We introduce a combination of four simple methods into the analysis of broadband seismograms so as to provide a means of improving phase recognition and the full use of broadband information for far-regional distances where the seismograms are particularly complex (because of the influence of the upper mantle discontinuities). For arrival detection we can employ the energy ratios of the short term behaviour to the long-term trend, using the vertical component and horizontal components of unrotated seismic records. We also use auto-regressive analysis to endeavour to separate broadband records into three parts: the seismic signal, microseismic noise and white noise. The higher order auto- and cross-correlation coefficient (representing the similarity of waveform) can be used to identify the presence of seismic phases, by avoiding the influence of the relatively low order correlation of microseismic noise. For each broadband 3-component record a set of complex traces are constructed and then a variety of definitions of instantaneous phase and frequency can be exploited to separate the behaviour of signal and noise. The complex traces can also be used for polarisation analysis. The changes in the character of the eigenvectors are particularly helpful for recognising the phases of broadband records in the far-regional range. The individual methods are quite powerful but when used in combination can provide a very effective means of phase characterisation. 相似文献
18.
Extended reflectivity method for modelling the propagation of diffusive–viscous wave in dip‐layered media 
下载免费PDF全文
下载免费PDF全文 The reflectivity method plays an important role in seismic modelling. It has been used to model different types of waves propagating in elastic and anelastic media. The diffusive–viscous wave equation was proposed to investigate the relationship between frequency dependence of reflections and fluid saturation. It is also used to describe the attenuation property of seismic wave in a fluid‐saturated medium. The attenuation of diffusive–viscous wave is mainly characterised by the effective attenuation parameters in the equation. Thus, it is essential to obtain those parameters and further characterise the features of the diffusive–viscous wave. In this work, we use inversion method to obtain the effective attenuation parameters through quality factor to investigate the characteristics of diffusive–viscous wave by comparing with those of the viscoacoustic wave. Then, the reflection/transmission coefficients in a dip plane‐layered medium are studied through coordinate transform and plane‐wave theory. Consequently, the reflectivity method is extended to compute seismograms of diffusive–viscous wave in a dip plane multi‐layered medium. Finally, we present two models to simulate the propagation of diffusive–viscous wave in a dip plane multi‐layered medium by comparing the results with those in a viscoacoustic medium. The numerical results demonstrate the validity of our extension of reflectivity method to the diffusive–viscous medium. The numerical examples in both time domain and time–frequency domain show that the reflections from a dip plane interface have significant phase shift and amplitude change compared with the results of horizontal plane interface due to the differences in reflection/transmission coefficients. Moreover, the modelling results show strong attenuation and phase shift in the diffusive–viscous wave compared to those of the viscoacoustic wave. 相似文献
19.
T. B. Yanovskaya V. M. Kozhevnikov O. A. Solovei K. R. Akchurin 《Izvestiya Physics of the Solid Earth》2008,44(8):622-630
Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region. 相似文献
20.
M.T. Carrozzo R. de Franco L. De Luca D. Luzio R. Primiceri T. Quarta M. Vitale 《Geophysical Prospecting》2002,50(6):547-564
A new filtering technique for single‐fold wide‐angle reflection/refraction seismic data is presented. The technique is based on the wavelet decomposition of a set of adjacent traces followed by coherence analysis. The filtering procedure consists of three steps. In the first, a wavelet decomposition of traces into different detail levels is performed. In the second, the coherence attributes for each level are evaluated by calculating cross‐correlation functions of detail portions contained in a space–time moving window. Finally, the filtered traces are obtained as a weighted reconstruction of the trace details. Each weight is obtained from the coherence‐attributes distribution estimated in a proper interval. A sequence of tests is then conducted in order to select possible optimum or unsuitable wavelet bases. The efficiency of the filter proposed was assessed by calculating some properly designed parameters in order to compare it with other standard de‐noising techniques. The proposed method produced a clear signal enhancement in high‐density wide‐angle seismic data, thus proving that it is a useful processing tool for a reliable correlation of seismic phases. 相似文献