低信噪比转换波地震资料静校正(英文)

Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points （CCP）. Some statics methods, which work well for P-wave data, can＇t be effectively used for solving convertedwave statics problems. This has become the main obstacle to breakthroughs in convertedwave data processing. To improve converted-wave static corrections, first, a statics method based on the common-receiver-point （CRP） stack is used for the initial receiver static correction to enhance the coherency of the CRP stack. Second, a stack-power-maximization static correction which improves the continuity of the CCP stack is used for detailed receiver statics. Finally, a non-surface-consistent residual moveout correction of the CCP gathers is used to enhance the stack power of reflection signals from different depths. Converted-wave statics are solved by the joint use of the three correction methods.     

基于偏移成像道集的剩余静校正方法

针对陆上地震资料处理的静校正问题,提出了一种基于偏移成像道集的剩余静校正方法.与传统的由动校正后的CMP道集中拾取剩余时差不同,本文基于偏移成像道集求取剩余时差,避免了复杂情况下同相轴归位不准确导致的剩余时差拾取误差.通过生成随炮点和检波点位置变化的偏移道集,实现了由偏移道集中直接拾取炮、检点的地表一致性剩余时差;该炮、检点偏移道集只在指定的局部时窗生成,并不增加大的计算量.二维和三维实际数据测试表明了该方法的有效性和实用性.     

三种算法联合迭代反演求取最佳剩余静校正量

估算剩余静校正量可看作是以激发点及接收点剩余静校正量为模型参数的、以叠加能量为目标函数的全局最优值的搜索问题．当地震记录的剩余静校正量大而信噪比较低时，应用相关统计法不能有效地拾取静校正量．本文分别采用了ＣＭＰ道集内道间互相关法、图像边缘提取的Ｌａｐｌａｃｉａｎ算法及ＤＦＰ算法，三种不同算法联合迭代反演可求取最佳的模型参数，使目标函数值达到或接近于全局最优．理论试算及实际资料处理效果表明，这种方法能够解决低信噪比及大静校正量造成的周期跳跃等静校正问题．     

起伏地表煤田地震资料静校正

由于地表起伏和近地表结构变化产生的静校正问题严重影响了煤田地震资料的成像质量.为此,首先利用低速带分片拟合的广义线性反演技术进行折射波静校正,解决长波长静校正问题和部分短波长静校正问题,然后,利用叠加能量最大静校正技术进一步解决剩余静校正问题,最后,利用非地表一致性剩余时差校正技术,解决速度和射线等误差引起的非地表一致性剩余时差问题.实验结果表明,在以串连的方式应用了三种校正方法之后,在共炮点道集上,折射渡同相轴的线性形态得到了恢复;在动校正后的共中心点道集上,煤层反射的双曲线同相轴被拉平;在叠加剖面上,煤层反射的信噪比得到了改善.     

Determination of a shallow velocity–depth model from seismic refraction data by coherence inversion1

Seismic refractions have different applications in seismic prospecting. The travel- times of refracted waves can be observed as first breaks on shot records and used for field static calculation. A new method for constructing a near-surface model from refraction events is described. It does not require event picking on prestack records and is not based on any approximation of arrival times. It consists of the maximization of the semblance coherence measure computed using shot gathers in a time window along refraction traveltimes. Time curves are generated by ray tracing through the model. The initial model for the inversion was constructed by the intercept-time method. Apparent velocities and intercept times were taken from a refraction stacked section. Such a section can be obtained by appling linea moveout corrections to common-shot records. The technique is tested successfully on synthetic and real data. An important application of the proposed method for solving the statics problem is demonstrated.     

3D seismic residual statics solutions derived from refraction interferometry

We apply interferometric theory to solve a three‐dimensional seismic residual statics problem to improve reflection imaging. The approach calculates the static solutions without picking the first arrivals from the shot or receiver gathers. The static correction accuracy can be significantly improved by utilising stacked virtual refraction gathers in the calculations. Shots and receivers may be placed at any position in a three‐dimensional seismic land survey. Therefore, it is difficult to determine stationary shots and receivers to form the virtual refraction traces that have identical arrival times, as in a two‐dimensional scenario. To overcome this problem, we use a three‐dimensional super‐virtual interferometry method for residual static calculations. The virtual refraction for a stationary shot/receiver pair is obtained via an integral along the receiver/shot lines, which does not require knowledge of the stationary locations. We pick the maximum energy times on the interferometric stacks and solve a set of linear equations to derive reliable residual static solutions. We further apply the approach to both synthetic and real data.     

基于共姿态道集的静校正方法

静校正问题是地震勘探的关键问题,直接影响地震勘探精度和准确性.实际地震采集过程中,当在相同接收点位置上不同时间内插拔布设了不同的检波器时,对于目前基于地表一致性理论假设的基准面静校正和剩余静校正,以及非地表一致性剩余静校正都不具备适用条件.为解决这一问题,本文提出了基于共姿态道集的静校正方法,将相同接收点位置上不同时间布设的检波点所接收的地震数据抽成不同的共姿态道集,在共姿态道集内实施地表一致性静校正;当某接收点位置上具有若干个共姿态道集时,该接收点位置上可能会存在多个检波点静校量;炮点静校问题仍然采用地表一致性静校正方法解决.该方法解决了同一接收点位置上不同共姿态道集之间的非地表一致性静校正问题,同时也解决了全区的检波点和炮点的地表一致性静校正问题,在实际数据应用效果明显.     

EFFICIENT MULTICHANNEL FILTERING OF SEISMIC DATA1

Multichannel filters are used to eliminate coherent noise from surface seismic data, for wavefield separation from VSP stacks, and for signal enhancement. Their success generally depends on the choice of the filter parameters and the domain of application. Multichannel filters can be applied to shots (monitors), common-receiver traces, CDP traces and stacked sections. Cascaded applications in these domains are currently performed in the seismic industry for better noise suppression and for signal enhancement. One-step shot-domain filtering is adequate for some applications. However, in practice, cascaded applications in shot-and common-receiver domains usually give better results when the S/N ratio is low. Multichannel filtering after stacking (especially after repeated applications in shot and/or receiver domains) may create undesirable results such as artificial continuations, or smearing and smoothing of small features such as small throw faults and fine stratigraphic details. Consequently, multichannel filtering after stacking must be undertaken with the utmost care and occasionally only as a last resort. Multichannel filters with fan-shaped responses (linear moveout filters) should be applied after NMO correction. These are the filters commonly used in the seismic industry where they have such names as velocity filters, moveout filters, f-k filters and coherency filters. Filtering before NMO correction may result in break-up and flattening especially of those shallow reflection events with relatively higher curvatures and diffractions. NMO correction is needed prior to wavefield separation from VSP stacks for the same practical reasons outlined above whenever source-receiver offsets are involved. Creation of artificial lineup and smearing at the outputs of multichannel filters is presently the common practical concern. Optimum multichannel filters with well-defined pass, reject and transition bands overcome the latter problems when applied before stacking and after NMO correction. The trace dimension of these filters must be kept small to avoid such lineups and the smoothing of small structures. Good results can be obtained with only five traces, but seven traces seems to be a better compromise both in surface and well seismic applications. The so-called f-k filtering and τ-p domain filtering are no exceptions to the above practical considerations. Residual static computations after multichannel filtering also need special consideration. Since multichannel filtering improves spatial continuity, residual static algorithms using local correlation, i.e. nonsurface-consistent algorithms, may be impractical especially after multichannel filtering.     

CRS-stack-based residual static correction

In the case of onshore data sets, the acquired reflection events can be strongly impaired due to rough top‐surface topography and inhomogeneities in the uppermost low‐velocity layer, the so‐called weathering layer. Without accounting for these influences, the poor data quality will make data processing very difficult. Usually, the correction for the top‐surface topography is not perfect. The residuals from this correction and the influence of the weathering layers lead to small distortions along the reflection events. We integrated a residual static correction method into our data‐driven common‐reflection‐surface‐stack‐based imaging workflow to further eliminate such distortions. The moveout‐corrected traces and the stacked pilot trace are cross‐correlated to determine a final estimate of the surface‐consistent residual statics in an iterative manner. As the handling of top‐surface topography within the common‐reflection‐surface stack is discussed in a separate paper in this special issue, the corresponding residual static correction will be explained in more detail. For this purpose, the results obtained with a data set from the Arabian Peninsula will be presented.     

HORIZONTAL STACKING AND MULTICHANNEL FILTERING APPLIED TO COMMON DEPTH POINT SEISMIC DATA*

The common depth point method of shooting in oil exploration provides a series of seismic traces which yield information about the substrata layers at one location. After normal moveout and static corrections have been applied, the traces are combined by horizontal stacking, or linear multichannel filtering, into a single record in which the primary reflections have been enhanced relative to the multiple reflections and random noise. The criterion used in optimum horizontal stacking is to maximize the signal to noise power ratio, where signal refers to the primary reflection sequence and noise includes the multiple reflections. It is shown when this criterion is equivalent to minimizing the mean square difference between the desired signal (primary reflection sequence) and the weighted horizontally stacked traces. If the seismic traces are combined by multichannel linear filtering, the primary reflection sequence will have undergone some phase and frequency distortion on the resulting record. The signal to noise power ratio then becomes less meaningful a criterion for designing the optimum linear multichannel filter, and the mean square criterion is adopted. In general, however, since more a priori information about the seismic traces is required to design the optimum linear multichannel filter than required for the optimum set of weights of the horizontal stacking process, the former will be an improvement over the latter. It becomes evident that optimum horizontal stacking is a restricted form of linear multichannel filtering.     

HORIZONTAL STACKING AND MULTICHANNEL FILTERING APPLIED TO COMMON DEPTH POINT SEISMIC DATA*

The common depth point method of shooting in oil exploration provides a series of seismic traces which yield information about the substrata layers at one location. After normal moveout and static corrections have been applied, the traces are combined by horizontal stacking, or linear multichannel filtering, into a single record in which the primary reflections have been enhanced relative to the multiple reflections and random noise. The criterion used in optimum horizontal stacking is to maximize the signal to noise power ratio, where signal refers to the primary reflection sequence and noise includes the multiple reflections. It is shown when this criterion is equivalent to minimizing the mean square difference between the desired signal (primary reflection sequence) and the weighted horizontally stacked traces. If the seismic traces are combined by multichannel linear filtering, the primary reflection sequence will have undergone some phase and frequency distortion on the resulting record. The signal to noise power ratio then becomes less meaningful a criterion for designing the optimum linear multichannel filter, and the mean square criterion is adopted. In general, however, since more a priori information about the seismic traces is required to design the optimum linear multichannel filter than required for the optimum set of weights of the horizontal stacking process, the former will be an improvement over the latter. It becomes evident that optimum horizontal stacking is a restricted form of linear multichannel filtering.     

CRS strategies for solving severe static and imaging issues in seismic data from Saudi Arabia

Static shifts from near‐surface inhomogeneities very often represent the key problem in the processing of seismic data from arid regions. In this case study, the deep bottom fill of a wadi strongly degrades the image quality of a 2D seismic data set. The resulting static and dynamic problems are solved by both conventional and common‐reflection‐surface (CRS) processing. A straightforward approach derives conventional refraction statics from picked first breaks and then goes through several iterations of manual velocity picking and residual statics calculation. The surface‐induced static and dynamic inhomogeneities, however, are not completely solved by these conventional methods. In CRS processing, the local adaptation of the CRS stacking parameters results in very detailed dynamic corrections. They resolve the local inhomogeneities that were not detected by manual picking of stacking velocities and largely compensate for the surface‐induced deterioration in the stack. The subsequent CRS residual statics calculations benefit greatly from the large CRS stacking fold which increases the numbers of estimates for single static shifts. This improves the surface‐consistent averaging of static shifts and the convergence of the static solution which removes the remaining static shifts in the 2D seismic data. The large CRS stacking fold also increases the signal‐to‐noise ratio in the final CRS stack.     

厚风化层覆盖区转换波静校正方法

P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.     

Dip‐adaptive singular‐value decomposition filtering for seismic reflection enhancement

We present a singular value decomposition (SVD) filtering method for the enhancement of coherent reflections and for attenuation of noise. The method is applied in two steps. First normal move‐out (NMO) correction is applied to shot or CMP records, with the purpose of flattening the reflections. We use a spatial SVD filter with a short sliding window to enhance coherent horizontal events. Then the data are sorted in common‐offset panels and the local dip is estimated for each panel. The next SVD filtering is performed on a small number of traces and a small number of time samples centred around the output sample position. Data in a local window are corrected for linear moveout corresponding to the dips before SVD. At the central time sample position, we sum over the dominant eigenimages of a few traces, corresponding to SVD dip filtering. We illustrate the method using land seismic data from the Tacutu basin, located in the north‐east of Brazil. The results show that the proposed method is effective and is able to reveal reflections masked by ground‐roll and other types of noise.     

三维地震勘探数据的剩余静校正方法研究与软件实现

剩余静校正的目的是消除地震道中由近地表因素造成的静态剩余时差的影响,通常以时差的四因子分解模型为基础建立剩余静校正方程组,并采用迭代求解方法获得剩余静校正量。但实际上,由于地层的倾斜,时差中还应包含横向和纵向倾角项,因此本文首先对剩余校正的理论模型进行了改进,提出了以扩展面元为计算单元的时差分解六因子模型。又由于地层的走向具有局部线性性,因此剩余动校正项,横向及纵向倾角项系数仅在一些给定的面元上计     

CRS-stack-based seismic imaging considering top-surface topography

In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity.     

提高浅层地震反射静校正精度的一种方法与数值模拟

浅层地震反射法是一种常用的勘探方法.在浅层地震资料处理中,静校正的精度直接影响速度反演的结果和叠加剖面的质量,在地形平缓时,固定基准面静校正可以满足勘探精度的要求,但在复杂地形条件下,其存在较大误差,即使采用浮动基准面,仍会由于地表一致性假设而残余静校正量,不能消除地形起伏引起的影响,为了提高浅层地震反射静校正的精度必须在常规静校正后进行一次剩余静校正,本文给出起伏地形条件下,滑动基准面(过共中心点的水平面即为该共中心点的滑动基准面)的剩余静校正量,该校正量与炮检距、反射层埋深、地层波速以及炮点和接收点高程有关,适用于单一介质和层状介质情况,本文通过对典型地形起伏的3个水平均匀层状介质理论地质模型的速度谱计算和分析,阐明在复杂地形条件下,应用本文提出的剩余静校正方法可以消除地形起伏的影响,提高静校正精度,在此基础上做动校正可以得到高质量的水平叠加剖面.     

Reverse time migration of CMP-gathers an effective tool for the determination of interval velocities

One of the most important steps in the conventional processing of reflection seismic data is common midpoint (CMP) stacking. However, this step has considerable deficiencies. For instance the reflection or diffraction time curves used for normal moveout corrections must be hyperbolae. Furthermore, undesirable frequency changes by stretching are produced on account of the dependence of the normal moveout corrections on reflection times. Still other drawbacks of conventional CMP stacking could be listed.One possibility to avoid these disadvantages is to replace conventional CMP stacking by a process of migration to be discussed in this paper. For this purpose the Sherwood-Loewenthal model of the exploding reflector has to be extended to an exploding point model with symmetry to the lineP _EX M whereP _EX is the exploding point, alias common reflection point, andM the common midpoint of receiver and source pairs.Kirchhoff summation is that kind of migration which is practically identical with conventional CMP stacking with the exception that Kirchhoff summation provides more than one resulting trace.In this paper reverse time migration (RTM) was adopted as a tool to replace conventional CMP stacking. This method has the merit that it uses the full wave equation and that a direct depth migration is obtained, the velocityv can be any function of the local coordinatesx, y, z. Since the quality of the reverse time migration is highly dependent on the correct choice of interval velocities such interval velocities can be determined stepwise from layer to layer, and there is no need to compute interval velocities from normal moveout velocities by sophisticated mathematics or time consuming modelling. It will be shown that curve velocity interfaces do not impair the correct determination of interval velocities and that more precise velocity values are obtained by avoiding or restricting muting due to non-hyperbolic normal moveout curves.Finally it is discussed how in the case of complicated structures the reverse time migration of CMP gathers can be modified in such a manner that the combination of all reverse time migrated CMP gathers yields a correct depth migrated section. This presupposes, however, a preliminary data processing and interpretation.     

LONG PERIOD STATICS DETERMINATION BY INVERSE FILTERING*

The relative time shift of CDP reflections, caused by long period statics, is treated as the result of a filtering process. The shape of the filter is that of a comb having the magnitude and spacing of its teeth determined by the CDP spread geometry. Convolution of a near surface anomaly with a comb filter representing the CDP spread geometry produces the same time fluctuations at CDP reflections (CDP statics) as the stacked traces across the anomaly. Conversely, the near surface time anomaly is given by the convolution of the CDP statics with the inverse of the comb filter. Provided that CDP statics can be separated from noise, dip, and residual normal moveout, it is pcissible to determine long period statics with a relative wavelength of between half a spreadlength and five spreadlengths. The object of this paper is to present the theory and practical applications of the method. Several examples based on synthetic and real data will also be discussed.     

SUPERSTACK—AN ITERATIVE STACKING ALGORITHM*

An algorithm for separation of signals according to their coherency is proposed. The algorithm, called Superstack, is used on common depth point data after normal moveout corrections have been applied. The algorithm can be regarded as an iterative stacking procedure. After each stack, input values are changed depending on the consistency of the output of the previous stack. The Superstack algorithm is able to provide better separation of signals showing a different degree of horizontal consistency (coherency) than the normal horizontal stack. All noise signals in a CDP-gather may be assumed to show less horizontal consistency than the primary reflections. On this basis the algorithm will give an improved primary/multiple and an improved primary/noise ratio. The Superstack is very general, should not be especially expensive in terms of computer time, and can easily be implemented in routine processing.