基于格林函数理论的波场预测和鬼波压制方法

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

Evolution of deghosting process for single‐sensor streamer data from 2D to 3D

In marine acquisition, reflections of sound energy from the water–air interface result in ghosts in the seismic data, both in the source side and the receiver side. Ghosts limit the bandwidth of the useful signal and blur the final image. The process to separate the ghost and primary signals, called the deghosting process, can fill the ghost notch, broaden the frequency band, and help achieve high‐resolution images. Low‐signal‐to‐noise ratio near the notch frequencies and 3D effects are two challenges that the deghosting process has to face. In this paper, starting from an introduction to the deghosting process, we present and compare two strategies to solve the latter. The first is an adaptive mechanism that adjusts the deghosting operator to compensate for 3D effects or errors in source/receiver depth measurement. This method does not include explicitly the crossline slowness component and is not affected by the sparse sampling in the same direction. The second method is an inversion‐type approach that does include the crossline slowness component in the algorithm and handles the 3D effects explicitly. Both synthetic and field data examples in wide azimuth acquisition settings are shown to compare the two strategies. Both methods provide satisfactory results.     

Joint up/down decomposition and reconstruction using three‐component streamers with or without ghost model: the sampling theory

This paper addresses two artefacts inherent to marine towed‐streamer surveys: 1) ghost reflections and 2) too sparse a sampling in the crossline direction. A ghost reflection is generated when an upcoming reflection bounces off the sea surface back into the sensors and can, in principle, be removed by decomposing the measured wavefield into its up‐ and downgoing constituents. This process requires a dense sampling of the wavefield in both directions along and perpendicular to the streamers. A dense sampling in the latter direction is, however, often impossible due to economical and operational constraints. Recent multi‐component streamers have been designed to record the spatial gradients on top of the pressure, which not only benefits the wavefield decomposition but also facilitates a lower‐than‐Nyquist sampling rate of the pressure. In this paper, wavefield reconstruction and deghosting are posed as a joint inverse problem. We present two approaches to establish a system matrix that embeds both a deghosting and an interpolation operator. The first approach is derived with a ghost model, whereas the second approach is derived without a ghost model. The embodiment of a ghost model leads to an even lower sampling rate but relies on a more restrictive assumption on the sea surface.     

Three‐dimensional receiver deghosting of seismic streamer data using L1 inversion and redundant extended radon dictionary

In this paper, we propose a novel three‐dimensional receiver deghosting algorithm that is capable of deghosting both horizontal and slanted streamer data in a theoretically consistent manner. Our algorithm honours wave propagation phenomena in a true three‐dimensional sense and frames the three‐dimensional receiver deghosting problem as a Lasso problem. The ultimate goal is to minimise the mismatch between the actual measurements and the simulated wavefield with an L1 constraint applied in the extended Radon space to handle the underdetermined nature of this problem. We successfully demonstrate our algorithm on a modified three‐dimensional EAGE/SEG Overthrust model and a Red Sea marine dataset.     

基于波动方程的上下缆地震数据鬼波压制方法研究

本文发展基于波动方程的上下缆鬼波压制方法,推导了上下缆地震波场频率波数域波动方程延拓合并公式.基于Fourier变换的波场解析延拓确保上下缆资料振幅相位的一致性,消除了长拖缆远偏移距信号的计算误差,同时具有较高的计算效率;上下缆地震波场的波动方程法合并有效解偶鬼波干涉,实现综合利用上下缆地震数据压制鬼波.理论模型数据和实际采集地震数据的测试表明了方法的有效性.     

频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

海上倾斜缆采集技术具有多样的陷波特征,通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题,本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法,该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明,频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波,能达到拓宽地震记录频带目的.     

Noise transfer in variable‐depth streamer deghosting

Single‐component towed‐streamer marine data acquisition records the pressure variations of the upgoing compressional waves followed by the polarity‐reversed pressure variations of downgoing waves, creating sea‐surface ghost events in the data. The sea‐surface ghost for constant‐depth towed‐streamer marine data acquisition is usually characterised by a ghost operator acting on the upgoing waves, which can be formulated as a filtering process in the frequency–wavenumber domain. The deghosting operation, usually via the application of the inverse Wiener filter related to the ghost operator, acts on the signal as well as the noise. The noise power transfer into the deghosted data is proportional to the power spectrum of the inverse Wiener filter and is amplifying the noise strongly at the notch wavenumbers and frequencies of the ghost operator. For variable‐depth streamer acquisition, the sea‐surface ghost cannot be described any longer as a wavenumber–frequency operator but as a linear relationship between the wavenumber–frequency representation of the upgoing waves at the sea surface and the data in the space–frequency domain. In this article, we investigate how the application of the inverse process acts on noise. It turns out that the noise magnification is less severe with variable‐depth streamer data, as opposed to constant depth, and is inversely proportional to the local slant of the streamer. We support this statement via application of the deghosting process to real and numerical random noise. We also propose a more general concept of a wavenumber–frequency ghost power transfer function, applicable for variable‐depth streamer acquisition, and demonstrate that the inverse of the proposed variable‐depth ghost power transfer function can be used to approximately quantify the action of the variable‐depth streamer deghosting process on noise.     

A probabilistic model for ghost delay time estimation based on recording geometry

In order to deconvolve the ghost response from marine seismic data, an estimate of the ghost operator is required. Typically, this estimate is made using a model of in‐plane propagation, i.e., the ray path at the receiver falls in the vertical plane defined by the source and receiver locations. Unfortunately, this model breaks down when the source is in a crossline position relative to the receiver spread. In this situation, in‐plane signals can only exist in a small region of the signal cone. In this paper, we use Bayes' theory to model the posterior probability distribution functions for the vertical component of the ray vector given the known source–receiver azimuth and the measured inline component of the ray vector. This provides a model for the ghost delay time based on the acquisition geometry and the dip of the wave in the plane of the streamer. The model is fairly robust with regard to the prior assumptions and controlled by a single parameter that is related to the likelihood of in‐plane propagation. The expected values of the resulting distributions are consistent with the deterministic in‐plane model when in‐plane likelihood is high but valid everywhere in the signal cone. Relaxing the in‐plane likelihood to a reasonable degree radically simplifies the shape of the expected‐value surface, lending itself for use in deghosting algorithms. The model can also be extended to other plane‐wave processing problems such as interpolation.     

Marine,seabed, and land seismic equipment for broadband acquisition: a review

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.     

Robust dual‐sensor data sea‐surface ghost attenuation for quality control purposes

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.     

最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用

海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的.     

Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.     

基于编码解码思想的鬼波预测与压制方法

高分辨率的宽带地震勘探技术是最近几年海上油气地震勘探的热点问题,鬼波压制是其中的核心议题.石油工业界,除了提出各种压制鬼波的采集方式外,资料处理过程中压制鬼波的方法也在不断地推陈出新.源、检鬼波的存在使得反射子波的有效频带变窄,成像分辨率降低,且干扰后续的自由表面多次波压制、FWI等.不同于常规的鬼波压制模型,本文基于编码与解码理论框架,用编码建立起鬼波预测模型,在Bayes反演框架下建立起解码估计一次波、从而压制鬼波的方法.基于此鬼波预测模型及相应反演理论的鬼波压制方法对鬼波的压制更为彻底.理论模型数据与实际资料测试结果验证了本文提出的理论框架和方法技术的有效性与优越性.     

Deghosting towed streamer data in τ/p domain based on rough sea surface reflectivity

Currently, the deghosting of towed streamer seismic data assumes a flat sea level and a sea-surface reflection coefficient of ?1; this decreases the precision of deghosting. A new method that considers the rough sea surface is proposed to suppress ghost reflections. The proposed deghosting method obtains the rough sea surface reflection coefficient using Gaussian statistics, and calculates the optimized deghosting operator in the τ/p domain. The proposed method is closer to the actual sea conditions, offers an improved deghosting operator, removes the ghost reflections from marine towed seismic data, widens the bandwidth and restores the low-frequency information, and finally improves the signal-tonoise ratio and resolution of the seismic data.     

Theoretical study on multi‐level source design

This paper studies aspects that influence the de‐ghosting performance of marine multi‐level sources based on a modified Johnson model. The normalized squared error between actual signature and its corresponding ghost‐free signature is introduced to evaluate the multi‐level source design. The results show that optimum depth combinations and volume combinations exist in the design. However, there is also some flexibility in the volume combination which makes it possible to optimize with respect to residual bubble oscillation. By considering both operational aspects and performance, we propose that three or four levels in a multi‐level source are reasonable. Compared to a horizontal source, a multi‐level source can be designed to reduce the notch effect, strengthen the down‐going energy and improve the energy transmission directivity. Studies of the influence of depth and firing time deviations indicate that a multi‐level source is more stable than a normal horizontal source in an operational environment.     

Source location using time‐reverse imaging

We present the chain of time‐reverse modeling, image space wavefield decomposition and several imaging conditions as a migration‐like algorithm called time‐reverse imaging. The algorithm locates subsurface sources in passive seismic data and diffractors in active data. We use elastic propagators to capitalize on the full waveforms available in multicomponent data, although an acoustic example is presented as well. For the elastic case, we perform wavefield decomposition in the image domain with spatial derivatives to calculate P and S potentials. To locate sources, the time axis is collapsed by extracting the zero‐lag of auto and cross‐correlations to return images in physical space. The impulse response of the algorithm is very dependent on acquisition geometry and needs to be evaluated with point sources before processing field data. Band‐limited data processed with these techniques image the radiation pattern of the source rather than just the location. We present several imaging conditions but we imagine others could be designed to investigate specific hypotheses concerning the nature of the source mechanism. We illustrate the flexible technique with synthetic 2D passive data examples and surface acquisition geometry specifically designed to investigate tremor type signals that are not easily identified or interpreted in the time domain.     

Extending illumination using all multiples: application to 3D acquisition geometry analysis

Recent advances in survey design have led to conventional common‐midpoint‐based analysis being replaced by subsurface‐based seismic acquisition analysis, with emphasis on advanced techniques of illumination analysis. Among them is the so‐called focal beam method, which is a wave‐equation‐based seismic illumination analysis method. The objective of the focal beam method is to provide a quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on the resolution and angle‐dependent amplitude fidelity of the image. The method distinguishes between illumination and sensing capability of a particular acquisition geometry by computing the focal source beam and the focal detector beam, respectively. Sensing is related to the detection properties of a detector configuration, whereas illumination is related to the emission properties of a source configuration. The focal source beam analyses the incident wavefield at a specific subsurface grid point from all available sources, whereas the focal detector beam analyses the sensing wavefield reaching at the detector locations from the same subsurface grid point. In the past, this method could only address illumination by primary reflections. In this paper, we will extend the concept of the focal beam method to incorporate the illumination due to the surface and internal multiples. This in fact complies with the trend of including multiples in the imaging process. Multiple reflections can illuminate a target location from other angles compared with primary reflections, resulting in a higher resolution and an improved illumination. We demonstrate how an acquisition‐related footprint can be corrected using both the surface and the internal multiples.     

How rough sea affects marine seismic data and deghosting procedures

Most seismic processing algorithms generally consider the sea surface as a flat reflector. However, acquisition of marine seismic data often takes place in weather conditions where this approximation is inaccurate. The distortion in the seismic wavelet introduced by the rough sea may influence (for example) deghosting results, as deghosting operators are typically recursive and sensitive to the changes in the seismic signal. In this paper, we study the effect of sea surface roughness on conventional (5–160 Hz) and ultra‐high‐resolution (200–3500 Hz) single‐component towed‐streamer data. To this end, we numerically simulate reflections from a rough sea surface using the Kirchhoff approximation. Our modelling demonstrates that for conventional seismic frequency band sea roughness can distort results of standard one‐dimensional and two‐dimensional deterministic deghosting. To mitigate this effect, we introduce regularisation and optimisation based on the minimum‐energy criterion and show that this improves the processing output significantly. Analysis of ultra‐high‐resolution field data in conjunction with modelling shows that even relatively calm sea state (i.e., 15 cm wave height) introduces significant changes in the seismic signal for ultra‐high‐frequency band. These changes in amplitude and arrival time may degrade the results of deghosting. Using the field dataset, we show how the minimum‐energy optimisation of deghosting parameters improves the processing result.     

Local outlier factor as part of a workflow for detecting and attenuating blending noise in simultaneously acquired data

A number of deblending methods and workflows have been reported in the past decades to eliminate the source interference noise recorded during a simultaneous shooting acquisition. It is common that denoising algorithms focusing on optimizing coherency and weighting down/ignoring outliers can be considered as deblending tools. Such algorithms are not only enforcing coherency but also handling outliers either explicitly or implicitly. In this paper, we present a novel approach based on detecting amplitude outliers and its application on deblending based on a local outlier factor that assigns an outlier-ness (i.e. a degree of being an outlier) to each sample of the data. A local outlier factor algorithm quantifies outlier-ness for an object in a data set based on the degree of isolation compared with its locally neighbouring objects. Assuming that the seismic pre-stack data acquired by simultaneous shooting are composed of a set of non-outliers and outliers, the local outlier factor algorithm evaluates the outlier-ness of each object. Therefore, we can separate the data set into blending noise (i.e. outlier) and signal (i.e. non-outlier) components. By applying a proper threshold, objects having high local outlier factors are labelled as outlier/blending noise, and the corresponding data sample could be replaced by zero or a statistically adequate value. Beginning with an explanation of parameter definitions and properties of local outlier factor, we investigate the feasibility of a local outlier factor application on seismic deblending by analysing the parameters of local outlier factor and suggesting specific deblending strategies. Field data examples recorded during simultaneous shooting acquisition show that the local outlier factor algorithm combined with a thresholding can detect and attenuate blending noise. Although the local outlier factor application on deblending shows a few shortcomings, it is consequently noted that the local outlier factor application in this paper obviously achieves benefits in terms of detecting and attenuating blending noise and paves the way for further geophysical applications.     

Acquisition design for direct reflectivity and velocity estimation from blended and irregularly sampled data

Blended acquisition along with efficient spatial sampling is capable of providing high-quality seismic data in a cost-effective and productive manner. While deblending and data reconstruction conventionally accompany this way of data acquisition, the recorded data can be processed directly to estimate subsurface properties. We establish a workflow to design survey parameters that account for the source blending as well as the spatial sampling of sources and detectors. The proposed method involves an iterative scheme to derive the survey design leading to optimum reflectivity and velocity estimation via joint migration inversion. In the workflow, we extend the standard implementation of joint migration inversion to cope with the data acquired in a blended fashion along with irregular detector and source geometries. This makes a direct estimation of reflectivity and velocity models feasible without the need of deblending or data reconstruction. During the iterations, the errors in reflectivity and velocity estimates are used to update the survey parameters by integrating a genetic algorithm and a convolutional neural network. Bio-inspired operators enable the simultaneous update of the blending and sampling operators. To relate the choice of survey parameters to the performance of joint migration inversion, we utilize a convolutional neural network. The applied network architecture discards suboptimal solutions among newly generated ones. Conversely, it carries optimal ones to the subsequent step, which improves the efficiency of the proposed approach. The resultant acquisition scenario yields a notable enhancement in both reflectivity and velocity estimation attributable to the choice of survey parameters.