共查询到20条相似文献，搜索用时 31 毫秒
1.
Emmanuel Spadavecchia Vincenzo Lipari Nicola Bienati Giuseppe Drufuca 《Geophysical Prospecting》2013,61(4):725-734
Despite being less general than 3D surface‐related multiple elimination (3D‐SRME), multiple prediction based on wavefield extrapolation can still be of interest, because it is less CPU and I/O demanding than 3D‐SRME and moreover it does not require any prior data regularization. Here we propose a fast implementation of water‐bottom multiple prediction that uses the Kirchhoff formulation of wavefield extrapolation. With wavefield extrapolation multiple prediction is usually obtained through the cascade of two extrapolation steps. Actually by applying the Fermat’s principle (i.e., minimum reflection traveltime) we show that the cascade of two operators can be replaced by a single approximated extrapolation step. The approximation holds as long as the water bottom is not too complex. Indeed the proposed approach has proved to work well on synthetic and field data when the water bottom is such that wavefront triplications are negligible, as happens in many practical situations. 相似文献
2.
Extrapolating wavefields and imaging at each depth during three‐dimensional recursive wave‐equation migration is a time‐consuming endeavor. For efficiency, most commercial techniques extrapolate wavefields through thick slabs followed by wavefield interpolation within each thick slab. In this article, we develop this strategy by associating more efficient interpolators with a Fourier‐transform‐related wavefield extrapolation method. First, we formulate a three‐dimensional first‐order separation‐of‐variables screen propagator for large‐step wavefield extrapolation, which allows for wide‐angle propagations in highly contrasting media. This propagator significantly improves the performance of the split‐step Fourier method in dealing with significant lateral heterogeneities at the cost of only one more fast Fourier transform in each thick slab. We then extend the two‐dimensional Kirchhoff and Born–Kirchhoff local wavefield interpolators to three‐dimensional cases for each slab. The three‐dimensional Kirchhoff interpolator is based on the traditional Kirchhoff formula and applies to moderate lateral velocity variations, whereas the three‐dimensional Born–Kirchhoff interpolator is derived from the Lippmann–Schwinger integral equation under the Born approximation and is adapted to highly laterally varying media. Numerical examples on the three‐dimensional salt model of the Society of Exploration Geophysicists/European Association of Geoscientists demonstrate that three‐dimensional first‐order separation‐of‐variables screen propagator Born–Kirchhoff depth migration using thick‐slab wavefield extrapolation plus thin‐slab interpolation tolerates a considerable depth‐step size of up to 72 ms, eventually resulting in an efficiency improvement of nearly 80% without obvious loss of imaging accuracy. Although the proposed three‐dimensional interpolators are presented with one‐way Fourier extrapolation methods, they can be extended for applications to general migration methods. 相似文献
3.
A new method for predicting different kinds of multiples and peg-leg reflections in unstacked seismic data is discussed. The basis for this method is the fact that kinematic properties of multiples can be represented as a combination of kinematic properties of primary reflections. The prediction is made using a two-step process. In the first step, the values for the angle of emergence and radius of curvature of the wavefront for primary reflections from ‘multiple-generating’ interfaces are obtained. These parameters are estimated directly from unstacked data for every source point using the homeomorphic-imaging technique. The second step consists of prediction of multiples from primary reflections that satisfy a so-called ‘multiple condition’. This condition is the equality of the absolute values of the angles of emergence calculated from the first step. This method is effective even in complex media and information on the subsurface geology is not required. The parameters are estimated directly from the unstacked data and do not require any computational efforts such as in wavefield extrapolation of data. 相似文献
4.
Antoine Guitton 《Geophysical Prospecting》2006,54(2):135-152
Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries. 相似文献
5.
Flavio Poletto Cinzia Bellezza Kees Wapenaar Aronne Craglietto 《Geophysical Prospecting》2012,60(5):838-854
We formulate the Kirchhoff‐Helmholtz representation theory for the combination of seismic interferometry signals synthesized by cross‐correlation and by cross‐convolution in acoustic media. The approach estimates the phase of the virtual reflections from the boundary encompassing a volume of interest and subtracts these virtual reflections from the total seismic‐interferometry wavefield. The reliability of the combination result, relevant for seismic exploration, depends on the stationary‐phase and local completeness in partial coverage regions. The analysis shows the differences in the phase of the corresponding seismic interferometry (by cross‐correlation) and virtual reflector (by cross‐convolution) signals obtained by 2D and 3D formulations, with synthetic examples performed to remove water layer multiples in ocean bottom seismic (OBS) acoustic data. 相似文献
6.
Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm 
下载免费PDF全文
下载免费PDF全文 Gareth Shane O'Brien Sean Joseph Delaney Michael Igoe John Doherty Andrew Colhoun 《Geophysical Prospecting》2018,66(4):688-706
Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study. 相似文献
7.
地下低速夹层的存在导致地震数据中包含较强能量的层间多次波,有效识别和预测深部储层上覆地层产生的层间多次波是提高深部储层解释精度的重要环节,而准确模拟层间多次波是辅助识别地震数据中层间多次波的一种非常有效的方法.本文提出了一种基于自适应变步长波场延拓的可控地层分阶层间多次波模拟方法,该方法基于自适应变步长波场延拓,以递归循环的方式实现分阶层间多次波的模拟.通过对模型添加双重层位约束,可以模拟指定地层产生的各阶层间多次波.利用二维反周期延拓方法压制波场延拓的边界反射优于传统方法,例如吸收边界法.提出自适应变步长波场延拓技术,大大提升了波场模拟的效率.理论和数值例子表明,本文方法模拟的一次波和各阶层间多次波与常用的有限差分方法模拟结果具有很好的一致性,且克服了有限差分方法无法分阶模拟波场的不足,显著提升了层间多次波识别的效率. 相似文献
8.
ARNIM B. HAASE 《Geophysical Prospecting》1992,40(4):403-427
Methods for predicting and attenuating water-bottom multiples by wavefield extrapolation have been discussed by several investigators. Because these prediction methods operate on shot records, boundary conditions must be specified for every shot record. The approach presented operates in the common-offset plane; a model of expected water-bottom multiples is generated from the observed surface wavefield using a finite-difference wave-equation migration algorithm with an offset term. An accurate water-depth profile is required, but there is no restriction on the shape of the water bottom other than a dip limit of approximately 18–20°. In generating a multiple model, the water-bottom primary and each water-bottom multiple reflection of the observed surface wavefield are extrapolated to a higher order. Thus, the extrapolated water-bottom primary of the model is lined up with a water-bottom multiple in the data and each multiple in the model is lined up with a higher-order (or later) multiple in the data. Prestack multiple attenuation is achieved, for one offset at a time, by first adapting the model of expected multiples to the observed data and then subtracting the predicted multiple energy. An error-constrained adaptation algorithm is proposed in order to control instabilities. No assumptions are made about primary reflections and no subwater-bottom velocities are required. Computational efficiency of modelling and adaptation can be improved by applying this method only to near and intermediate offsets as the stacking process usually provides sufficient multiple attenuation at far offsets. A field data example demonstrates the potential of the proposed method for improving the primary-to-multiple ratio in prestack and post-stack data. 相似文献
9.
Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated. 相似文献
10.
Two-dimensional VSP surveys are often conducted to provide structural illumination of the subsurface away from the borehole. The illumination is achieved through offsetting the source with respect to the downhole geophone. This inevitably gives rise to mode-conversions in both downgoing and upgoing wavefields. Migration of mixed-mode wavefields is complex because the velocity profile used for wavefield extrapolation is valid only for a particular propagation mode; the other mode always propagates at a different velocity. It is therefore advisable to separate the wave-types (P-wave and SV-wave) prior to migration. This may be achieved through wavemode filtering, a multichannel process which exploits the relation between propagation velocity, slowness of events at the recording array and particle motion. The necessary information about particle motion is available only if the VSP data are acquired with a three-component downhole geophone assembly. The wavemode filter partitions wave-types at the recording array; it provides no information about the various changes of propagation mode experienced by the energy as it travels from source to geophone. For the purpose of migration, the intermediate modes of propagation must be deduced. Much of the energy arriving at the receivers is P-wave which has followed the P-wave velocity profile from the source. It can therefore be imaged by conventional (Kirchhoff) migration. As an example of SV-wave imaging, a common mode-code is P-wave from source to reflector and SV-wave from reflector to geophone. Migration of such data calls for back-propagation of the geophone array wavefield, at SV-wave velocity, to the point in the subsurface where it is time-coincident with the forward propagated downwave, at P-wave velocity. 相似文献
11.
12.
A.J. Berkhout 《Geophysical Prospecting》2017,65(1):106-145
Surface removal and internal multiple removal are explained by recursively separating the primary and multiple responses at each depth level with the aid of wavefield prediction error filtering. This causal removal process is referred to as “data linearization.” The linearized output (primaries only) is suitable for linear migration algorithms. Next, a summary is given on the migration of full wavefields (primaries + multiples) by using the concept of secondary sources in each subsurface gridpoint. These secondary sources are two‐way and contain the gridpoint reflection and the gridpoint transmission properties. In full wavefield migration, a local inversion process replaces the traditional linear imaging conditions. Finally, Marchenko redatuming is explained by iteratively separating the full wavefield response from above a new datum and the full wavefield response from below a new datum. The redatuming output is available for linear migration (Marchenko imaging) or, even better, for full wavefield migration. Linear migration, full wavefield migration, and Marchenko imaging are compared with each other. The principal conclusion of this essay is that multiples should not be removed, but they should be utilized, yielding two major advantages: (i) illumination is enhanced, particularly in the situation of low signal‐to‐noise primaries; and (ii) both the upper side and the lower side of reflectors are imaged. It is also concluded that multiple scattering algorithms are more transparent if they are formulated in a recursive depth manner. In addition to transparency, a recursive depth algorithm has the flexibility to enrich the imaging process by inserting prior geological knowledge or by removing numerical artefacts at each depth level. Finally, it is concluded that nonlinear migration algorithms must have a closed‐loop architecture to allow successful imaging of incomplete seismic data volumes (reality of field data). 相似文献
13.
For increasing time values, isochrons can be regarded as expanding wavefronts and their perpendicular lines as the associated orthogonal isochron rays. The speed of the isochron movement depends on the medium velocity and the source-receiver position. We introduce the term equivalent-velocity to refer to the speed of isochron movement. In the particular case of zero-offset data, the equivalent velocity is half of the medium velocity. We use the concepts of orthogonal isochron-rays and equivalent velocity to extend the application of the exploding reflector model to non-zero offset imaging problems. In particular, we employ these concepts to extend the use of zero-offset wave-equation algorithms for modelling and imaging common-offset sections. In our imaging approach, the common-offset migration is implemented as a trace-by-trace algorithm in three steps: equivalent velocity computation, data conditioning for zero-offset migration and zero-offset wave-equation migration. We apply this methodology for modelling and imaging synthetic common-offset sections using two kinds of algorithms: finite-difference and split-step wavefield extrapolation. We also illustrate the isochron-ray imaging methodology with a field-data example and compare the results with conventional common-offset Kirchhoff migration. This methodology is attractive because it permits depth migration of common-offset sections or just pieces of that by using wave-equation algorithms, it extends the use of robust zero-offset algorithms, it presents favourable features for parallel processing, it permits the creation of hybrid migration algorithms and it is appropriate for migration velocity analysis. 相似文献
14.
综述了当前多次波问题的研究进展,指出了对多次波问题,有两种主要的处理思路：(1)因为后续偏移成像的需要,把多次波看成相干噪声,根据多次波特性设计衰减方法将其消除掉,将其称之为多次波衰减方法.多次波衰减方法大体上可以分为两类：基于信号分析的滤波方法和基于波动方程的预测相减法.预测相减法不需要或较少需要有关地下的先验或后验信息,能适应复杂地下构造.本文对预测相减法中的波场延拓法、反馈迭代法、逆散射级数法和恒定内插法等三种主要方法,进行了较详细的对比分析和介绍.(2)将多次波看成有效信号,对其成像,将其称之为多次波成像方法.这种方法提取了多次波中包含的地下界面信息,可改善地震资料的成像质量. 相似文献
15.
Tariq Alkhalifah 《Geophysical Prospecting》2015,63(5):1126-1141
Wavefield extrapolation operators for elliptically anisotropic media offer significant cost reduction compared with that for the transversely isotropic case, particularly when the axis of symmetry exhibits tilt (from the vertical). However, elliptical anisotropy does not provide accurate wavefield representation or imaging for transversely isotropic media. Therefore, we propose effective elliptically anisotropic models that correctly capture the kinematic behaviour of wavefields for transversely isotropic media. Specifically, we compute source‐dependent effective velocities for the elliptic medium using kinematic high‐frequency representation of the transversely isotropic wavefield. The effective model allows us to use cheaper elliptic wave extrapolation operators. Despite the fact that the effective models are obtained by matching kinematics using high‐frequency asymptotic, the resulting wavefield contains most of the critical wavefield components, including frequency dependency and caustics, if present, with reasonable accuracy. The methodology developed here offers a much better cost versus accuracy trade‐off for wavefield computations in transversely isotropic media, particularly for media of low to moderate complexity. In addition, the wavefield solution is free from shear‐wave artefacts as opposed to the conventional finite‐difference‐based transversely isotropic wave extrapolation scheme. We demonstrate these assertions through numerical tests on synthetic tilted transversely isotropic models. 相似文献
16.
通常工业界实现逆时偏移算法时采用有限差分数值方法模拟地震波场,波场模拟常常受稳定性条件限制,且易产生数值频散,成像精度降低.本文引入了一步法波场延拓方法,首先构建声波传播算子,借助Chebyshev多项式和Jacobi-Anger展开式近似传播算子中的e指数项,进而实现波场递推,该方法时间步长的选取不受稳定性条件限制而且不存在空间频散现象.本文将一步法波场延拓方法用于逆时偏移成像的波场模拟,并提出双缓冲区存储策略,在不增加计算量的前提下,大幅降低了逆时偏移方法的波场存储量.波场模拟和逆时偏移成像测试表明,本文提出的一步法波场延拓方法模拟地震波场精度高,消除了频散影响,可在较大时间步长的情况下实现高精度波场模拟;提出的基于一步法波场延拓的逆时偏移方法成像质量好;基于双缓冲区存储策略的逆时偏移成像方法存储成本低. 相似文献
17.
Riemannian wavefield extrapolation is a technique for one‐way extrapolation of acoustic waves. Riemannian wavefield extrapolation generalizes wavefield extrapolation by downward continuation by considering coordinate systems different from conventional Cartesian ones. Coordinate systems can conform with the extrapolated wavefield, with the velocity model or with the acquisition geometry. When coordinate systems conform with the propagated wavefield, extrapolation can be done accurately using low‐order kernels. However, in complex media or in cases where the coordinate systems do not conform with the propagating wavefields, low order kernels are not accurate enough and need to be replaced by more accurate, higher‐order kernels. Since Riemannian wavefield extrapolation is based on factorization of an acoustic wave‐equation, higher‐order kernels can be constructed using methods analogous to the one employed for factorization of the acoustic wave‐equation in Cartesian coordinates. Thus, we can construct space‐domain finite‐differences as well as mixed‐domain techniques for extrapolation. High‐order Riemannian wavefield extrapolation kernels improve the accuracy of extrapolation, particularly when the Riemannian coordinate systems does not closely match the general direction of wave propagation. 相似文献
18.
VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性. 相似文献
19.
Source–receiver two‐way wave extrapolation for prestack exploding‐reflector modelling and migration 下载免费PDF全文
下载免费PDF全文 Most modern seismic imaging methods separate input data into parts (shot gathers). We develop a formulation that is able to incorporate all available data at once while numerically propagating the recorded multidimensional wavefield forward or backward in time. This approach has the potential for generating accurate images free of artiefacts associated with conventional approaches. We derive novel high‐order partial differential equations in the source–receiver time domain. The fourth‐order nature of the extrapolation in time leads to four solutions, two of which correspond to the incoming and outgoing P‐waves and reduce to the zero‐offset exploding‐reflector solutions when the source coincides with the receiver. A challenge for implementing two‐way time extrapolation is an essential singularity for horizontally travelling waves. This singularity can be avoided by limiting the range of wavenumbers treated in a spectral‐based extrapolation. Using spectral methods based on the low‐rank approximation of the propagation symbol, we extrapolate only the desired solutions in an accurate and efficient manner with reduced dispersion artiefacts. Applications to synthetic data demonstrate the accuracy of the new prestack modelling and migration approach. 相似文献
20.
J.F.B. Bolte 《Geophysical Prospecting》2004,52(6):523-534
The common focal point (CFP) method and the common reflection surface (CRS) stack method are compared. The CRS method is a fast, highly automated procedure that provides high S/N ratio simulation of zero‐offset (ZO) images by combining, per image point, the reflection energy of an arc segment that is tangential to the reflector. It uses smooth parametrized two‐way stacking operators, based on a data‐driven triplet of attributes in 2D (eight parameters in 3D). As a spin‐off, the attributes can be used for several applications, such as the determination of the geometrical spreading factor, multiple prediction, and tomographic inversion into a smooth background velocity model. The CFP method aims at decomposing two‐way seismic reflection data into two full‐aperture one‐way propagation operators. By applying an iterative updating procedure in a half‐migrated domain, it provides non‐smooth focusing operators for prestack imaging using only the energy from one focal point at the reflector. The data‐driven operators inhibit all propagation effects of the overburden. The CFP method provides several spin‐offs, amongst which is the CFP matrix related to one focal point, which displays the reflection amplitudes as measured at the surface for each source–receiver pair. The CFP matrix can be used to determine the specular reflection source–receiver pairs and the Fresnel zone at the surface for reflection in one single focal point. Other spin‐offs are the prediction of internal multiples, the determination of reflectivity effects, velocity‐independent redatuming and tomographic inversion to obtain a velocity–depth model. The CFP method is less fast and less automated than the CRS method. From a pointwise comparison of features it is concluded that one method is not a subset of the other, but that both methods can be regarded as being to some extent complementary. 相似文献