共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
Paul Sava 《Geophysical Prospecting》2015,63(5):1050-1069
Imaging the change in physical parameters in the subsurface requires an estimate of the long wavelength components of the same parameters in order to reconstruct the kinematics of the waves propagating in the subsurface. One can reconstruct the model by matching the recorded data with modeled waveforms extrapolated in a trial model of the medium. Alternatively, assuming a trial model, one can obtain a set of images of the reflectors from a number of seismic experiments and match the locations of the imaged interfaces. Apparent displacements between migrated images contain information about the velocity model and can be used for velocity analysis. A number of methods are available to characterize the displacement between images; in this paper, we compare shot‐domain differential semblance (image difference), penalized local correlations, and image‐warping. We show that the image‐warping vector field is a more reliable tool for estimating displacements between migrated images and leads to a more robust velocity analysis procedure. By using image‐warping, we can redefine the differential semblance optimization problem with an objective function that is more robust against cycle‐skipping than the direct image difference. We propose an approach that has straightforward implementation and reduced computational cost compared with the conventional adjoint‐state method calculations. We also discuss the weakness of migration velocity analysis in the migrated‐shot domain in the case of highly refractive media, when the Born modelling operator is far from being unitary and thus its adjoint (migration) operator poorly approximates the inverse. 相似文献
3.
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. 相似文献
4.
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. 相似文献
5.
The estimation of a velocity model from seismic data is a crucial step for obtaining a high‐quality image of the subsurface. Velocity estimation is usually formulated as an optimization problem where an objective function measures the mismatch between synthetic and recorded wavefields and its gradient is used to update the model. The objective function can be defined in the data‐space (as in full‐waveform inversion) or in the image space (as in migration velocity analysis). In general, the latter leads to smooth objective functions, which are monomodal in a wider basin about the global minimum compared to the objective functions defined in the data‐space. Nonetheless, migration velocity analysis requires construction of common‐image gathers at fixed spatial locations and subsampling of the image in order to assess the consistency between the trial velocity model and the observed data. We present an objective function that extracts the velocity error information directly in the image domain without analysing the information in common‐image gathers. In order to include the full complexity of the wavefield in the velocity estimation algorithm, we consider a two‐way (as opposed to one‐way) wave operator, we do not linearize the imaging operator with respect to the model parameters (as in linearized wave‐equation migration velocity analysis) and compute the gradient of the objective function using the adjoint‐state method. We illustrate our methodology with a few synthetic examples and test it on a real 2D marine streamer data set. 相似文献
6.
Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix–vector formulation from which different image gathers can be formed and with which amplitude‐versus‐angle and wave‐equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two‐way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave‐equation‐based migration‐velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips. 相似文献
7.
Migration velocity analysis and waveform inversion 总被引：3，自引：0，他引：3
William W. Symes 《Geophysical Prospecting》2008,56(6):765-790
Least‐squares inversion of seismic reflection waveform data can reconstruct remarkably detailed models of subsurface structure and take into account essentially any physics of seismic wave propagation that can be modelled. However, the waveform inversion objective has many spurious local minima, hence convergence of descent methods (mandatory because of problem size) to useful Earth models requires accurate initial estimates of long‐scale velocity structure. Migration velocity analysis, on the other hand, is capable of correcting substantially erroneous initial estimates of velocity at long scales. Migration velocity analysis is based on prestack depth migration, which is in turn based on linearized acoustic modelling (Born or single‐scattering approximation). Two major variants of prestack depth migration, using binning of surface data and Claerbout's survey‐sinking concept respectively, are in widespread use. Each type of prestack migration produces an image volume depending on redundant parameters and supplies a condition on the image volume, which expresses consistency between data and velocity model and is hence a basis for velocity analysis. The survey‐sinking (depth‐oriented) approach to prestack migration is less subject to kinematic artefacts than is the binning‐based (surface‐oriented) approach. Because kinematic artefacts strongly violate the consistency or semblance conditions, this observation suggests that velocity analysis based on depth‐oriented prestack migration may be more appropriate in kinematically complex areas. Appropriate choice of objective (differential semblance) turns either form of migration velocity analysis into an optimization problem, for which Newton‐like methods exhibit little tendency to stagnate at nonglobal minima. The extended modelling concept links migration velocity analysis to the apparently unrelated waveform inversion approach to estimation of Earth structure: from this point of view, migration velocity analysis is a solution method for the linearized waveform inversion problem. Extended modelling also provides a basis for a nonlinear generalization of migration velocity analysis. Preliminary numerical evidence suggests a new approach to nonlinear waveform inversion, which may combine the global convergence of velocity analysis with the physical fidelity of model‐based data fitting. 相似文献
8.
We develop a new time‐domain reverse‐time migration method called double plane‐wave reverse‐time migration that uses plane‐wave transformed gathers. Original shot gathers with appropriate data acquisition geometry are double slant stacked into the double plane‐wave domain with minimal slant stacking artefacts. The range of plane‐wave components needed for migration can be determined by estimating the maximum time dips present in shot gathers. This reduces the total number of input traces for migration and increases migration efficiency. Unlike the pre‐stack shot‐profile reverse‐time migration where the number of forward propagations is proportional to the number of shots, the number of forward propagations needed for the proposed method remains constant and is relatively small even for large seismic datasets. Therefore, the proposed method can improve the efficiency of the migration and be suitable for migrating large datasets. Double plane‐wave reverse‐time migration can be performed for selected plane‐wave components to obtain subsurface interfaces with different dips, which makes the migration method target oriented. This feature also makes the method a useful tool for migration velocity analysis. For example, we are able to promptly obtain trial images with nearly horizontal interfaces and adjust velocity models according to common image gathers. Seismic signal coming from steeply dipping interfaces can be included into the migration to build images with more detailed structures and higher spatial resolution as better velocity models become available. Illumination compensation imaging conditions for the proposed method are also introduced to obtain images with balanced amplitudes. 相似文献
9.
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. 相似文献
10.
Using both image and data domains to perform velocity inversion can help us resolve the long and short wavelength components of the velocity model, usually in that order. This translates to integrating migration velocity analysis into full waveform inversion. The migration velocity analysis part of the inversion often requires computing extended images, which is expensive when using conventional methods. As a result, we use pre‐stack wavefield (the double‐square‐root formulation) extrapolation, which includes the extended information (subsurface offsets) naturally, to make the process far more efficient and stable. The combination of the forward and adjoint pre‐stack wavefields provides us with update options that can be easily conditioned to improve convergence. We specifically use a modified differential semblance operator to split the extended image into a residual part for classic differential semblance operator updates and the image (Born) modelling part, which provides reflections for higher resolution information. In our implementation, we invert for the velocity and the image simultaneously through a dual objective function. Applications to synthetic examples demonstrate the features of the approach. 相似文献
11.
方位角度域共成像点道集能够客观反映地下介质的速度、各向异性参数异常以及振幅随角度变化（AVA）和裂缝信息。传统Kirchhoff PSTM通常输出偏移距域共成像点道集,对于速度分析、各向异性分析、AVA分析、裂缝识别等均存在诸多不便。本文提出了基于走时梯度的Kirchhoff叠前时间偏移全方位角度集输出方法并提出工业上切实可行的实现方案。通过走时场梯度计算波场传播方向矢量,形成能够反映观测系统参数和波场传播情况的全方位角度域共成像点道集。为了在大规模地震数据Kirchhoff积分叠前时间偏移中输出全方位角度道集,本文给出基于输入道方式的偏移实现方法,采用逐条inline线进行线偏移成像,从而大大降低了全方位角度道集输出对计算机内存的压力,显著提高了Kirchhoff积分时间偏移输出全方位角度道集的可行性。三维盐丘模型测试和海上某区块三维实际资料试验证明了本文方法的正确性。 相似文献
12.
Pre-stack Kirchhoff depth migration local space-shift imaging condition: synthetic and data examples
Extracting accurate common image angle gathers from pre-stack depth migrations is important in the generation of any incremental uplift to the amplitude versus angle attributes and seismic inversions that can lead to significant impacts in exploration and development success. The commonly used Kirchhoff migration outputs surface common offset image gathers that require a transformation to angle gathers for amplitude versus angle analysis. The accuracy of this transformation is one of the factors that determine the robustness of the amplitude versus angle measurements. Here, we investigate the possibility of implementing an extended imaging condition, focusing on the space-lag condition, for generating subsurface reflection angle gathers within a Kirchhoff migration. The objective is to determine if exploiting the spatial local shift imaging condition can provide any increase in angle gather fidelity relative to the common offset image gathers. The same restrictions with a ray-based approach will apply using the extended imaging condition as both the offset and extended imaging condition method use travel times derived from solutions to an Eikonal equation. The aims are to offer an alternative ray-based method to generate subsurface angle gathers and to understand the impact on the amplitude versus angle response. To this end, the implementation of the space-shift imaging condition is discussed and results of three different data sets are presented. A layered three-dimensional model and a complex two-dimensional model are used to assess the space shift image gathers output from such a migration scheme and to evaluate the seismic attributes relative to the traditional surface offset common image gathers. The synthetic results show that the extended imaging condition clearly provides an uplift in the measured amplitude versus angle over the surface offset migration. The noise profile post-migration is also improved for the space-lag migration due to the double summation inside the migration. Finally, we show an example of a space-lag gather from deep marine data and compare the resultant angle gathers with those generated from an offset migration and a time-shift imaging condition Kirchhoff migration. The comparison of the real data with a well log shows that the space-lag result is a better match to the well compared to the time-lag extended imaging condition and the common offset Kirchhoff migration. Overall, 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 with an incremental improvement in amplitude versus angle fidelity and lower noise but comes at a higher computational cost. 相似文献
13.
Wavefield‐based migration velocity analysis using the semblance principle requires computation of images in an extended space in which we can evaluate the imaging consistency as a function of overlapping experiments. Usual industry practice is to assemble those seismic images in common‐image gathers that represent reflectivity as a function of depth and extensions, e.g., reflection angles. We introduce extended common‐image point (CIP) gathers constructed only as a function of the space‐ and time‐lag extensions at sparse and irregularly distributed points in the image. Semblance analysis using CIP's constructed by this procedure is advantageous because we do not need to compute gathers at regular surface locations and we do not need to compute extensions at all depth levels. The CIP's also give us the flexibility to distribute them in the image at irregular locations aligned with the geologic structure. Furthermore, the CIP's remove the depth bias of common‐image gathers constructed as a function of the depth axis. An interpretation of the CIP's using the scattering theory shows that they are scattered wavefields associated with sources and receivers inside the subsurface. Thus, when the surface wavefields are correctly reconstructed, the extended CIP's are characterized by focused energy at the origin of the space‐ and time‐lag axes. Otherwise, the energy defocuses from the origin of the lag axes proportionally with the cumulative velocity error in the overburden. This information can be used for wavefield‐based tomographic updates of the velocity model, and if the velocity used for imaging is correct, the coordinate‐independent CIP's can be a decomposed as a function of the angles of incidence. 相似文献
14.
Seismic velocity analysis in the scattering-angle/azimuth domain 总被引：2，自引：0，他引：2
Migration velocity analysis is carried out by analysing the residual moveout and amplitude variations in common image point gathers (CIGs) parametrized by scattering angle and azimuth. The misfit criterion in the analysis is of the differential-semblance type. By using angles to parametrize the imaging we are able to handle and exploit data with multiple arrivals, although artefacts may occur in the CIGs and need to be suppressed. The CIGs are generated by angle migration, an approach based on the generalized Radon transform (GRT) inversion, and they provide multiple images of reflectors in the subsurface for a range of scattering angles and azimuths. Within the differential semblance applied to these CIGs, we compensate for amplitude versus angle (AVA) effects. Thus, using a correct background velocity model, the CIGs should have no residual moveout nor amplitude variation with angles, and the differential semblance should vanish. If the velocity model is incorrect, however, the events in the CIGs will appear at different depths for different angles and the amplitude along the events will be non-uniform. A standard, gradient-based optimization scheme is employed to develop a velocity updating procedure. The model update is formed by backprojecting the differential semblance misfits through ray perturbation kernels, within a GRT inverse. The GRT inverse acts on the data, subject to a shift in accordance with ray perturbation theory. The performance of our algorithm is demonstrated with two synthetic data examples using isotropic elastic models. The first one allows velocity variation with depth only. In the second one, we reconstruct a low-velocity lens in the model that gives rise to multipathing. The velocity model parametrization is based upon the eigentensor decomposition of the stiffness tensor and makes use of B-splines. 相似文献
15.
16.
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. 相似文献
17.
18.
提出三维Kirchhoff叠前时间偏移角度域共像点道集的改进算法,克服传统角度求取算法局限,可计算相对倾斜地层法线入射角;与Kirchhoff直射线叠前时间偏移求角度算法相比,本文方法考虑射线弯曲效应,包含层速度,角度范围加大,更接近真实入射角;计算走时采取弯曲射线或者适应线性横向变速介质的非对称走时等算法,角度道集在大角度处得到拉平;采用相对保幅的权因子以及覆盖次数校正技术,有利于叠前AVA反演.模型测试结果表明：叠前时间偏移角度道集,相对CMP、CRP所转化角度道集,更准确反应AVA效应;实际三维数据测试表明本文方法可以提供品质优良的角度道集,适用于AVA分析、反演,提高叠前反演分辨率. 相似文献
19.
逆时偏移(RTM)是复杂介质条件下地震成像的重要手段.因受观测系统限制、上覆地层影响以及波场带宽有限等因素的影响,现行的常规RTM所采用的互相关成像条件通常对地下构造进行模糊成像.最小二乘逆时偏移(LSRTM)通过最小化线性Born近似正演数据和采集数据之间的波形差异,采用梯度类反演算法优化反射系数模型,获得的成像结果具有更高的分辨率和更可靠的振幅保真度.然而,基于波形拟合的LSRTM对背景速度模型的依赖性很强.误差太大的速度模型容易产生周波跳跃现象,导致LSRTM难以获得全局最优解.为了克服这一问题,本文基于扩展模型的思想,在线性Born近似下,推导得到RTM扩展成像条件.并基于最小二乘反演理论,提出扩展成像条件下的LSRTM方法.理论模型试算表明,本文方法不仅可以提供分辨率更高、振幅属性更为可靠的成像结果,而且能够在一定程度上消除速度误差对反演成像的影响. 相似文献
20.
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. 相似文献