海洋单阶次表面多次波分离方法研究

海洋地震资料普遍发育强能量表面多次波,传统表面多次波压制技术(SRME)能够预测出所有阶次表面多次波,但是各阶次表面多次波相互混叠.为了能够单独利用不同阶次的表面多次波成像,降低干涉假象对多次波成像的影响,需要将不同阶次的表面多次波分离出来.本文提出一种基于扩展SRME的海洋单阶次表面多次波分离方法.首先,应用SRME...     

一阶多次波聚焦变换成像

将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率.     

表面多次波最小二乘逆时偏移成像

使用相同的炮记录,多次波偏移能提供比反射波偏移更广的地下照明和更多的地下覆盖但是同时产生很多的串声噪声.相比传统逆时偏移,最小二乘逆时偏移反演的反射波成像结果具有更高的分辨率和更均衡的振幅.我们主要利用最小二乘逆时偏移压制多次波偏移产生的串声噪声.多次波最小二乘逆时偏移通常需要一定的迭代次数以较好地消除串声噪声.若提前将一阶多次波从所有阶数的多次波中过滤出来,使用相同的迭代次数,一阶多次波的最小二乘逆时偏移能够得到具有更高信噪比的成像剖面,而且能够提供与多次波最小二乘逆时偏移相似的有效地下结构成像.     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

The influence of errors in the source wavelet on inversion‐based surface‐related multiple attenuation

We investigate the influence of source wavelet errors on inversion‐based, surface‐related multiple attenuation, in order to address how the inverted primary impulse response, estimated primaries, and predicted multiples are affected by the estimated wavelet. In theory, errors in estimated wavelet can lead to errors in the upgoing waves. Because of smoothness and the band‐limitedness characteristics of the estimated wavelet, errors in the upgoing waves are usually not white and random. Theoretical analysis and two synthetic examples demonstrate that (i) when the overall amplitude scalar of the estimated wavelet is underestimated, the inversion of the primary impulse response suffers from instability, which will distort the estimation of primaries, and (ii) when the wavelet is overestimated, the estimated primaries will simply mimic the recorded upgoing waves. Nevertheless, the quality of the estimated primaries in the region above the first‐order, water‐bottom multiples is independent of the estimated wavelet. Synthetic results illustrate that inversion‐based, surface‐related multiple attenuation with a known wavelet is stable, since slight inaccuracy in amplitude spectrum and/or phase spectrum of the given wavelet or the corresponding upgoing waves will not lead to considerable deviation in the waveforms of the inverted results from those of the references. Furthermore, shot‐to‐shot wavelet variations, with maximum amplitude difference of 5% and maximum phase difference of 10°, create just slight artefacts in both the inverted primary impulse response and the estimated primaries. Moreover, the sensitivity test of estimation of primaries by sparse inversion method involving wavelet estimation shows that this method can stably and alternately update the wavelet and the primary impulse response; however, different choices of the initial wavelet can lead to different final inverted results.     

基于稀疏反演的OBS数据多次波压制方法

自由表面多次波压制是海底地震仪（Ocean Bottom Seismometer,OBS）数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计（Estimation of Primaries and Multiples by Sparse Inversion,EPSI）方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME（Surface Related Multiple Elimination）方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

Closed-loop SRME based on 3D L1-norm sparse inversion

In many situations, the quality of seismic imaging is largely determined by a proper multiple attenuation as preprocessing step. Despite the widespread application of surface-related multiple elimination (SRME) and estimation of primaries by sparse inversion (EPSI) for the removal of multiples, there still exist some limitations in the process of prediction and subtraction (SRME) or inversion (EPSI), which make the efficiency of multiple attenuation less satisfactory. To solve these problems, a new fully data-driven method called closed-loop SRME was proposed, which combines the robustness of SRME and the multi-dimensional inversion strategy of EPSI. Due to the selection of inversion approach and constraint, primary estimation by closed-loop SRME may fall into a local optimum during the solving process, which lowers the accuracy of deep information and weakens the continuity of seismic events. To avoid these shortcomings, we first modified the solving method for closed-loop SRME to an L1 norm-based bi-convex optimization method, which stabilizes the solution. Meanwhile, in the L1 norm constraint-based optimization process, the 3D sparsifying transform, being a 2D Curvelet-1D wavelet transform, is brought in as a 3D sparse constraint. In the 3D sparsifying domain, the data become sparser, thus making the result of optimization more accurate, the information of seismic events more continuous and the resolution higher. Examples on both synthetic and field data demonstrate that the method proposed in this paper, compared with the traditional SRME and closed-loop SRME, have an excellent effect on primary estimation and suppress multiples effectively.     

3D targeted multiple attenuation

Short-period multiple reflections pose a particular problem in the North Sea where predictive deconvolution is often only partially successful. The targeted multiple attenuation (TMA) algorithm comprises computation of the covariance matrix of preflattened prestack or post-stack seismic data, the determination of the dominating eigenvectors of the covariance matrix, and subtraction of the related eigenimages followed by reverse flattening. The main assumption made is that the flattened multiple reflections may be represented by the first eigenimage(s) which implies that the spatial amplitude variations of primaries and associated multiples are similar. This assumption usually limits the method to short-period multiple reflections. TMA is applicable post-stack or prestack to common-offset gathers. It is computationally fast, robust towards random noise, irregular geometry and spatial aliasing, and it preserves the amplitudes of primaries provided they are not parallel to the targeted multiples. Application of TMA to 3D wavefields is preferable because this allows a better discrimination between primaries and multiples. Real data examples show that the danger of partially removing primary energy can be reduced by improving the raw multiple model that is based on eigenimages, for example by prediction filtering.     

Estimation and application of near‐surface full waveform redatuming operators

Land seismic data quality can be severely affected by near‐surface anomalies. The imprint of a complex near‐surface can be removed by redatuming the data to a level below the surface, from where the subsurface structures are assumed to be relatively smooth. However, to derive a velocity‐depth model that explains the propagation effects of the near‐surface is a non‐trivial task. Therefore, an alternative approach has been proposed, where the redatuming operators are obtained in a data‐driven manner from the reflection event related to the datum. In the current implementation, the estimation of these redatuming operators is done in terms of traveltimes only, based on a high‐frequency approximation. The accompanying amplitudes are usually derived from a local homogeneous medium, which is obviously a simplification of reality. Such parametrization has produced encouraging results in the past but cannot completely remove the near‐surface complexities, leaving artefacts in the redatumed results. In this paper we propose a method that estimates the redatuming operators directly from the data, i.e., without using a velocity model, in a full waveform manner, such that detailed amplitude and phase variations are included. The method directly outputs the inverse propagation operators that are needed for true‐amplitude redatuming. Based on 2D synthetic data it is demonstrated that the resulting redatuming quality is improved and artefacts are reduced.     

Repeatability enhancement in deep-water permanent seismic installations: a dynamic correction for seawater velocity variations

With the increasing use of permanently installed seismic installations, many of the issues in time‐lapse seismic caused by the lack of repeatability can be reduced. However, a number of parameters still influence the degree of reliability of 4D seismic data. In this paper, the specific impact of seawater velocity variations on time‐lapse repeatability is investigated in a synthetic study. A zero‐lag time‐lapse seabed experiment with no change in the subsurface but with velocity changes in the water column is simulated. The velocity model in the water column is constant for the baseline survey while the model for the repeat survey is heterogeneous, designed from sea salinity and temperature measurements in the West of Shetlands. The difference section shows up to 80% of residual amplitude, which highlights the poor repeatability. A new dynamic correction which removes the effect of seawater velocity variations specifically for permanent installations is developed. When applied to the synthetic data, it reduces the difference residual amplitude to about 3%. This technique shows substantial improvement in repeatability beyond conventional time‐lapse cross‐equalization.     

Interferometric interpolation of sparse marine data

We present the theory and numerical results for interferometrically interpolating 2D and 3D marine surface seismic profiles data. For the interpolation of seismic data we use the combination of a recorded Green's function and a model‐based Green's function for a water‐layer model. Synthetic (2D and 3D) and field (2D) results show that the seismic data with sparse receiver intervals can be accurately interpolated to smaller intervals using multiples in the data. An up‐ and downgoing separation of both recorded and model‐based Green's functions can help in minimizing artefacts in a virtual shot gather. If the up‐ and downgoing separation is not possible, noticeable artefacts will be generated in the virtual shot gather. As a partial remedy we iteratively use a non‐stationary 1D multi‐channel matching filter with the interpolated data. Results suggest that a sparse marine seismic survey can yield more information about reflectors if traces are interpolated by interferometry. Comparing our results to those of f‐k interpolation shows that the synthetic example gives comparable results while the field example shows better interpolation quality for the interferometric method.     

基于相关迭代的非因果匹配滤波器多次波压制方法

针对常规多次波压制(SRME)方法存在的缺陷进行了改进.首先在多次波模型道预测阶段,采用数据相关和迭代更新的策略预测多次波模型道,该方法降低了SRME方法采用数据空间褶积对输入数据要求严格的限制,提高了对近偏移距缺失和空间假频数据的适应性.在匹配相减阶段,本文设计了一种非因果非平稳的匹配滤波器,该滤波器可以对整道多次波模型进行处理,而且当预测多次波模型道滞后于实际多次波时,也能够对多次波模型道进行很好的匹配.模型数据和实际数据试算证明该方法在多次波模型道预测和匹配相减阶段的优越性.     

基于Marchenko理论压制自由表面多次波方法研究

常规虚源点Marchenko自聚焦多次波预测方法只适用于预测不含自由表面的多次波模型,局限于压制层间多次波,该方法在构建上下行格林函数场前,必须从反射响应中去除所有与表面相关的多次波.本文对构建上下行Marchenko格林函数方程进行改进,得到了包含一次波、层间多次波和自由表面多次波的格林函数,利用改进的Marchenko自聚焦预测方法预测自由表面多次波.本文利用水平层状模型数据及SMARRT模型数据证明,改进后的Marchenko法预测海底相关的自由表面多次波效果较为理想,该方法避免了常规SRME自由表面多次波预测方法需要近道重构的缺陷,能够有效提高地震资料的信噪比和分辨率.     

Reverse time migration using water-bottom-related multiples

Reverse time migration of multiples can be used to construct subsurface structures where primaries cannot illuminate well. However, the images generated using multiples suffer from severe artefacts due to the cross-talks created by interference among unrelated multiples. We developed a migration approach using water-bottom-related multiples to reduce these cross-talk artefacts. This approach first isolates primaries from the original data and predicts water-column primaries. The nth-order water-column multiples can be obtained by auto-convolving the water-column primaries n times, followed by convolving the nth-order water-column multiples with the primaries to extract the (n+1)th-order water-bottom-related multiples. The approach takes the nth-order water-column multiples as the secondary source and regards the (n+1)th-order water-bottom-related multiples as the receiver wavefield, followed by a cross-correlation imaging condition. Numerical examples from synthetic and field data sets reveal that our approach can provide images with substantially fewer cross-talk artefacts than conventional reverse time migration using multiples, as well as greatly improving shallow imaging compared with reverse time migration of primaries.     

3‐D shallow‐water seismic survey evaluation and design using the focal‐beam method: a case study offshore Abu Dhabi

For 3‐D shallow‐water seismic surveys offshore Abu Dhabi, imaging the target reflectors requires high resolution. Characterization and monitoring of hydrocarbon reservoirs by seismic amplitude‐versus‐offset techniques demands high pre‐stack amplitude fidelity. In this region, however, it still was not clear how the survey parameters should be chosen to satisfy the required data quality. To answer this question, we applied the focal‐beam method to survey evaluation and design. This subsurface‐ and target‐oriented approach enables quantitative analysis of attributes such as the best achievable resolution and pre‐stack amplitude fidelity at a fixed grid point in the subsurface for a given acquisition geometry at the surface. This method offers an efficient way to optimize the acquisition geometry for maximum resolution and minimum amplitude‐versus‐offset imprint. We applied it to several acquisition geometries in order to understand the effects of survey parameters such as the four spatial sampling intervals and apertures of the template geometry. The results led to a good understanding of the relationship between the survey parameters and the resulting data quality and identification of the survey parameters for reflection imaging and amplitude‐versus‐offset applications.     

The impact of surface‐wave attenuation on 3‐D seismic survey design

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.     

在成像空间中衰减多次波方法研究

在偏移后成像空间中的共成像点道集中可以对多次波进行衰减,对于给定的偏移速度模型,一次波与多次波在叠前偏移后的共成像点道集中具有不同的动校时差,这样我们就可以使用类似于偏移前衰减多次波的方法将一次波和多次波进行分离.本文在成像空间中应用抛物Radon变换分离多次波和有效波,由于每个共成像点道集都包含了复杂三维波场传播效应,所以本文方法具有处理三维数据和复杂地下构造的能力.相比于SRME以及传统Radon变换衰减多次波方法,本文方法能够在保持较小的计算量的同时,保证了衰减多次波的准确性.通过对模型数据试算和对实际数据的处理验证了本文方法在叠前时间偏移后衰减多次波的能力,并取得了很好的成像效果.     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

地表数据驱动的与层相关的层间多次波消除方法及应用

基于反馈模型的地表相关多次波消除方法SRME(Surface-Related Multiple Elimination)近年来已得到了广泛的应用.利用共聚焦点CFP(Common Focus Point)道集代替炮集,可以将该方法扩展至层间多次波的消除.地表数据驱动的层间多次波消除方法直接利用地表观测数据进行层间多次波预测,避免了构建CFP道集所需的聚焦运算,特别是与层相关的层间多次波消除方法有效提高了多次波预测的计算效率.但地表数据驱动的与层相关的层间多次波消除方法并没有从理论上被严格地推导证明,其与CFP方法之间的关系亦未被讨论.本文在CFP方法的理论基础上推导了地表数据驱动的与层相关的层间多次波消除方法,阐明了CFP方法与地表数据驱动方法之间的内在联系.并将该方法应用于模型数据和野外实际数据,应用实例表明了所提方法的有效性.