共查询到20条相似文献，搜索用时 500 毫秒
1.
We propose a method for imaging small‐scale diffraction objects in complex environments in which Kirchhoff‐based approaches may fail. The proposed method is based on a separation between the specular reflection and diffraction components of the total wavefield in the migrated surface angle domain. Reverse‐time migration was utilized to produce the common image gathers. This approach provides stable and robust results in cases of complex velocity models. The separation is based on the fact that, in surface angle common image gathers, reflection events are focused at positions that correspond to the apparent dip angle of the reflectors, whereas diffracted events are distributed over a wide range of angles. The high‐resolution radon‐based procedure is used to efficiently separate the reflection and diffraction wavefields. In this study, we consider poststack diffraction imaging. The advantages of working in the poststack domain are its numerical efficiency and the reduced computational time. The numerical results show that the proposed method is able to image diffraction objects in complex environments. The application of the method to a real seismic dataset illustrates the capability of the approach to extract diffractions. 相似文献
2.
Moshe Reshef 《Geophysical Prospecting》2014,62(6):1453-1467
We present an innovative approach for seismic image enhancement using multi‐parameter angle‐domain characterization of common image gathers. A special subsurface angle‐domain imaging system is used to generate the multi‐parameter common image gathers in a summation‐free image space. The imaged data associated with each common image gathers depth point contain direction‐dependent opening‐angle image contributions from all the available incident and scattered wave‐pairs at this point. Each direction‐dependent opening‐angle data can be differently weighted according to its coherency measure. Once the optimal migration velocity is used, it is assumed that in the actual specular direction, the coherency measure (semblance) along reflection events, from all available opening angles and opening azimuths, is larger than that along non‐specular directions. The computed direction‐dependent semblance attribute is designed to operate as an imaging filter which enhances specular migration contributions and suppresses all others in the final migration image. The ability to analyse the structural properties of the image points by the multi‐parameter common image gather allows us to better handle cases of complicated wave propagation and to improve the image quality at poorly illuminated regions or near complex structures. The proposed method and some of its practical benefits are demonstrated through detailed analysis of synthetic and real data examples. 相似文献
3.
地震绕射波是地下非连续性地质体的地震响应,绕射波成像对地下断层、尖灭和小尺度绕射体的识别具有重要的意义.在倾角域共成像点道集中,反射波同相轴表现为一条下凸曲线,能量主要集中在菲涅耳带内,绕射波能量则比较发散.由于倾角域菲涅耳带随偏移距变化而存在差异,因此本文提出一种在倾角-偏移距域道集中精确估计菲涅耳带的方法,在各偏移距的倾角域共成像点道集中实现菲涅耳带的精确切除,从而压制反射波.在倾角-偏移距域道集中还可以分别实现绕射波增强,绕射波同相轴相位校正,因此能量弱的绕射波可以清晰地成像.在倾角域共成像点道集中,反射波同相轴的最低点对应于菲涅耳带估计所用的倾角,因此本文提出一种在倾角域共成像点道集中直接自动拾取倾角场的方法.理论与实际资料试算验证了本文绕射波成像方法的有效性. 相似文献
4.
Extended common‐image‐point gathers (CIP) constructed by wide‐azimuth TI wave‐equation migration contain all the necessary information for angle decomposition as a function of the reflection and azimuth angles at selected locations in the subsurface. The aperture and azimuth angles are derived from the extended images using analytic relations between the space‐ and time‐lag extensions using information which is already available at the time of migration, i.e. the anisotropic model parameters. CIPs are cheap to compute because they can be distributed in the image at the most relevant positions, as indicated by the geologic structure. If the reflector dip is known at the CIP locations, then the computational cost can be reduced by evaluating only two components of the space‐lag vector. The transformation from extended images to angle gathers is a planar Radon transform which depends on the local medium parameters. This transformation allows us to separate all illumination directions for a given experiment, or between different experiments. We do not need to decompose the reconstructed wavefields or to choose the most energetic directions for decomposition. Applications of the method include illumination studies in complex areas where ray‐based methods fail, and assuming that the subsurface illumination is sufficiently dense, the study of amplitude variation with aperture and azimuth angles. 相似文献
5.
6.
Tariq Alkhalifah 《Geophysical Prospecting》2015,63(2):271-282
Reverse‐time migration can accurately image complex geologic structures in anisotropic media. Extended images at selected locations in the Earth, i.e., at common‐image‐point gathers, carry rich information to characterize the angle‐dependent illumination and to provide measurements for migration velocity analysis. However, characterizing the anisotropy influence on such extended images is a challenge. Extended common‐image‐point gathers are cheap to evaluate since they sample the image at sparse locations indicated by the presence of strong reflectors. Such gathers are also sensitive to velocity error that manifests itself through moveout as a function of space and time lags. Furthermore, inaccurate anisotropy leaves a distinctive signature in common‐image‐point gathers, which can be used to evaluate anisotropy through techniques similar to the ones used in conventional wavefield tomography. It specifically admits a V‐shaped residual moveout with the slope of the “V” flanks depending on the anisotropic parameter η regardless of the complexity of the velocity model. It reflects the fourth‐order nature of the anisotropy influence on moveout as it manifests itself in this distinct signature in extended images after handling the velocity properly in the imaging process. Synthetic and real data observations support this assertion. 相似文献
7.
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. 相似文献
8.
State‐of‐the‐art 3D seismic acquisition geometries have poor sampling along at least one dimension. This results in coherent migration noise that always contaminates pre‐stack migrated data, including high‐fold surveys, if prior‐to‐migration interpolation was not applied. We present a method for effective noise suppression in migrated gathers, competing with data interpolation before pre‐stack migration. The proposed technique is based on a dip decomposition of common‐offset volumes and a semblance‐type measure computation via offset for all constant‐dip gathers. Thus the processing engages six dimensions: offset, inline, crossline, depth, inline dip, and crossline dip. To reduce computational costs, we apply a two‐pass (4D in each pass) noise suppression: inline processing and then crossline processing (or vice versa). Synthetic and real‐data examples verify that the technique preserves signal amplitudes, including amplitude‐versus‐offset dependence, and that faults are not smeared. 相似文献
9.
Interval velocity analysis using post‐stack data has always been a desire, mainly for 3D data sets. In this study we present a method that uses the unique characteristics of migrated diffractions to enable interval velocity analysis from three‐dimensional zero‐offset time data. The idea is to perform a standard three‐dimensional prestack depth migration on stack cubes and generate three‐dimensional common image gathers that show great sensitivity to velocity errors. An efficient ‘top‐down’ scheme for updating the velocity is used to build the model. The effectiveness of the method is related to the incorporation of wave equation based post‐stack datuming in the model building process. The proposed method relies on the ability to identify diffractions along redatumed zero‐offset data and to analyse their flatness in the migrated local angle domain. The method can be considered as an additional tool for a complete, prestack depth migration based interval velocity analysis. 相似文献
10.
Unequal illumination of the subsurface highly impacts the quality of seismic imaging. Different image points receive different folds of reflection‐angle illumination, which can be caused by irregular acquisition or by wave propagation in complex media. Illumination problems can deteriorate amplitudes in migrated images. To address this problem, we present a method of stacking angle‐domain common‐image gathers, in which we use local similarity with soft thresholding to determine the folds of local illumination. Normalization by local similarity regularizes local illumination of reflection angles for each image point of the subsurface model. This approach compensates for irregular illumination by selective stacking in the image space, regardless of the cause of acquisition or propagation irregularities. Additional migration is not required because the methodology is implemented in the reflection angle domain after migration. We use two synthetic examples to demonstrate that our method can normalize migration amplitudes and effectively suppress migration artefacts. 相似文献
11.
C.J. Thomson 《Geophysical Prospecting》2012,60(1):64-92
The reflection operator for a simple flat‐lying interface can be thought of as the set of all its plane‐wave reflection coefficients or as the set of virtual surveys with sources and receivers along the interface. When there is dip, however, it is necessary to include the varying effects of propagation between the virtual‐survey level and the interface. Hence, step one in this paper is to derive the reflection operator for a dipping plane interface as observed at a datum level some distance away. The key assumption is that the aperture at the datum level is sufficient to characterize the reflector properties around a particular point. This translates into an assumption that the dip is moderate, though no explicit small‐angle approximation is required. The second step is to find the apparent reflection operator that would relate data that have been extrapolated from the datum towards and possibly beyond the reflector using an assumed migration velocity. This apparent reflection operator is closely related to extended common‐image gathers. The apparent reflection operator may be analysed asymptotically in terms of rays and other signals, shedding light on the structure of extended image gathers. In keeping with the virtual‐survey idea, the results are considered in a subsurface space‐time or slowness‐time domain at various extrapolation levels around the interface. An important distinction is drawn between using subsurface midpoint‐offset coordinates and the wavefield coordinates of the incident and reflected waves. The latter reveal more clearly the effects of dip, because they lead to a more asymmetric apparent reflection operator. Properties such as an up‐dip shift of a traveltime minimum and its associated curvature theoretically provide information about the reflector location and dip and the migration‐velocity error. The space‐time form of the reflection operator can be highly intricate around the offset‐time origin and it was described for a simple flat interface in a background paper. To avoid a layer of mathematics, the reflection‐operator formulas presented here are in the intermediate space‐frequency domain. They are analysed by considering their stationary‐phase and branch‐point high‐frequency contributions. There is no Born‐like assumption of weak reflector contrast and so wide‐angle, total reflection and head‐wave effects are included. Snell’s law is an explicit part of the theory. It is hoped that the work will therefore be a step towards the goal of unifying amplitude‐versus‐offset, imaging and waveform inversion. 相似文献
12.
Reverse‐time migration has become an industry standard for imaging in complex geological areas. We present an approach for increasing its imaging resolution by employing time‐shift gathers. The method consists of two steps: (i) migrating seismic data with the extended imaging condition to get time‐shift gathers and (ii) accumulating the information from time‐shift gathers after they are transformed to zero‐lag time‐shift by a post‐stack depth migration on a finer grid. The final image is generated on a grid, which is denser than that of the original image, thus improving the resolution of the migrated images. Our method is based on the observation that non‐zero‐lag time‐shift images recorded on the regular computing grid contain the information of zero‐lag time‐shift image on a denser grid, and such information can be continued to zero‐lag time‐shift and refocused at the correct locations on the denser grid. The extra computational cost of the proposed method amounts to the computational cost of zero‐offset migration and is almost negligible compared with the cost of pre‐stack shot‐record reverse‐time migration. Numerical tests on synthetic models demonstrate that the method can effectively improve reverse‐time migration resolution. It can also be regarded as an approach to improve the efficiency of reverse‐time migration by performing wavefield extrapolation on a coarse grid and by generating the final image on the desired fine grid. 相似文献
13.
地震绕射波源于介质非连续性,从地震记录中将绕射波分离出来并进行成像,其结果对研究诸如碳酸盐岩缝洞储层这类复杂非均质储层具有重要意义.对炮集记录进行平面波分解,在地层倾角不大的假设下,反射波和绕射波同相轴在平面波分解剖面上存在较大的倾角差异.基于此,我们提出分步进行绕射波分离的方法:(1)利用局部倾角滤波方法将绕射波的较大倾角信息成分分离出来,此时,余下的部分包含有反射波和残留的低倾角绕射波信息;(2)利用频率-空间域预测反演方法从上述含有反射波和残留的低倾角绕射波信息中分离出残留绕射波成分;(3)将两次分离的绕射波信息相加得到最终的绕射波估计.用该方法能够得到相对完整的绕射波信息,有效地克服了靠单一的倾角差异进行绕射波分离时明显损失低倾角信息,从而影响绕射波成像结果横向分辨率这一问题.理论与实际资料试算验证了该方法的有效性. 相似文献
14.
We present preserved‐amplitude downward continuation migration formulas in the aperture angle domain. Our approach is based on shot‐receiver wavefield continuation. Since source and receiver points are close to the image point, a local homogeneous reference velocity can be approximated after redatuming. We analyse this approach in the framework of linearized inversion of Kirchhoff and Born approximations. From our analysis, preserved‐amplitude Kirchhoff and Born inverse formulas can be derived for the 2D case. They involve slant stacks of filtered subsurface offset domain common image gathers followed by the application of the appropriate weighting factors. For the numerical implementation of these formulas, we develop an algorithm based on the true amplitude version of the one‐way paraxial approximation. Finally, we demonstrate the relevance of our approach with a set of applications on synthetic datasets and compare our results with those obtained on the Marmousi model by multi‐arrival ray‐based preserved‐amplitude migration. While results are similar, we observe that our results are less affected by artefacts. 相似文献
15.
In recent years, wave‐equation imaged data are often presented in common‐image angle‐domain gathers as a decomposition in the scattering angle at the reflector, which provide a natural access to analysing migration velocities and amplitudes. In the case of anisotropic media, the importance of angle gathers is enhanced by the need to properly estimate multiple anisotropic parameters for a proper representation of the medium. We extract angle gathers for each downward‐continuation step from converting offset‐frequency planes into angle‐frequency planes simultaneously with applying the imaging condition in a transversely isotropic with a vertical symmetry axis (VTI) medium. The analytic equations, though cumbersome, are exact within the framework of the acoustic approximation. They are also easily programmable and show that angle gather mapping in the case of anisotropic media differs from its isotropic counterpart, with the difference depending mainly on the strength of anisotropy. Synthetic examples demonstrate the importance of including anisotropy in the angle gather generation as mapping of the energy is negatively altered otherwise. In the case of a titled axis of symmetry (TTI), the same VTI formulation is applicable but requires a rotation of the wavenumbers. 相似文献
16.
Paul Sava 《Geophysical Prospecting》2015,63(5):1086-1096
Waveform inversion is a velocity‐model‐building technique based on full waveforms as the input and seismic wavefields as the information carrier. Conventional waveform inversion is implemented in the data domain. However, similar techniques referred to as image‐domain wavefield tomography can be formulated in the image domain and use a seismic image as the input and seismic wavefields as the information carrier. The objective function for the image‐domain approach is designed to optimize the coherency of reflections in extended common‐image gathers. The function applies a penalty operator to the gathers, thus highlighting image inaccuracies arising from the velocity model error. Minimizing the objective function optimizes the model and improves the image quality. The gradient of the objective function is computed using the adjoint state method in a way similar to that in the analogous data‐domain implementation. We propose an image‐domain velocity‐model building method using extended common‐image‐point space‐ and time‐lag gathers constructed sparsely at reflections in the image. The gathers are effective in reconstructing the velocity model in complex geologic environments and can be used as an economical replacement for conventional common‐image gathers in wave‐equation tomography. A test on the Marmousi model illustrates successful updating of the velocity model using common‐image‐point gathers and resulting improved image quality. 相似文献
17.
Claudio Strobbia Alexander Zarkhidze Roger May Fatma Ibrahim 《Geophysical Prospecting》2014,62(5):1143-1161
Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data. 相似文献
18.
Martin Tygel Bjørn Ursin Einar Iversen Maarten V. de Hoop 《Geophysical Prospecting》2012,60(2):201-216
Starting from a given time‐migrated zero‐offset data volume and time‐migration velocity, recent literature has shown that it is possible to simultaneously trace image rays in depth and reconstruct the depth‐velocity model along them. This, in turn, allows image‐ray migration, namely to map time‐migrated reflections into depth by tracing the image ray until half of the reflection time is consumed. As known since the 1980s, image‐ray migration can be made more complete if, besides reflection time, also estimates of its first and second derivatives with respect to the time‐migration datum coordinates are available. Such information provides, in addition to the location and dip of the reflectors in depth, also an estimation of their curvature. The expressions explicitly relate geological dip and curvature to first and second derivatives of reflection time with respect to time‐migration datum coordinates. Such quantitative relationships can provide useful constraints for improved construction of reflectors at depth in the presence of uncertainty. Furthermore, the results of image‐ray migration can be used to verify and improve time‐migration algorithms and can therefore be considered complementary to those of normal‐ray migration. So far, image‐ray migration algorithms have been restricted to layered models with isotropic smooth velocities within the layers. Using the methodology of surface‐to‐surface paraxial matrices, we obtain a natural extension to smooth or layered anisotropic media. 相似文献
19.
为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果. 相似文献
20.
Wave‐equation migration velocity analysis is a technique designed to extract and update velocity information from migrated images. The velocity model is updated through the process of optimizing the coherence of images migrated with the known background velocity model. The capacity for handling multi‐pathing of the technique makes it appropriate in complex subsurface regions characterized by strong velocity variation. Wave‐equation migration velocity analysis operates by establishing a linear relation between a slowness perturbation and a corresponding image perturbation. The linear relationship and the corresponding linearized operator are derived from conventional extrapolation operators and the linearized operator inherits the main properties of frequency‐domain wavefield extrapolation. A key step in the implementation is to design an appropriate procedure for constructing an image perturbation relative to a reference image that represents the difference between the current image and a true, or more correct image of the subsurface geology. The target of the inversion is to minimize such an image perturbation by optimizing the velocity model. Using time‐shift common‐image gathers, one can characterize the imperfections of migrated images by defining the focusing error as the shift of the focus of reflections along the time‐shift axis. The focusing error is then transformed into an image perturbation by focusing analysis under the linear approximation. As the focusing error is caused by the incorrect velocity model, the resulting image perturbation can be considered as a mapping of the velocity model error in the image space. Such an approach for constructing the image perturbation is computationally efficient and simple to implement. The technique also provides a new alternative for using focusing information in wavefield‐based velocity model building. Synthetic examples demonstrate the successful application of our method to a layered model and a subsalt velocity update problem. 相似文献