共查询到20条相似文献，搜索用时 184 毫秒
1.
非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理. 相似文献
2.
苏北大陆科学钻探靶区深反射地震的叠前深度偏移 总被引：4，自引：2，他引：2
由于深反射地震数据具有信噪比低和记录长度长等特点,叠前深度偏移方法的应用有许多困难．为此,我们研究了一种适合于深反射地震的叠前深度偏移方法;包括：用逆风有限差分方法计算程函方程;在常规速度扫描的基础上,用协方差控制提高速度分析精度;用联合反演算法计算层速度,再插值后得到初始速度模型;用Ｋｉｒｃｈｈｏｆｆ法作为偏移速度分析工具,求得最终的速度模型;最终的速度模型作为有限差分深度偏移的输入,求得最终的偏移结果．用该方法对“中国大陆科学深钻工程”东海二维深反射地震数据ＤＨ－４线进行了叠前深度偏移,取得了良好的效果。 相似文献
3.
Seismic tomography is a well‐established approach to invert smooth macro‐velocity models from kinematic parameters, such as traveltimes and their derivatives, which can be directly estimated from data. Tomographic methods differ more with respect to data domains than in the specifications of inverse‐problem solving schemes. Typical examples are stereotomography, which is applied to prestack data and Normal‐Incidence‐Point‐wave tomography, which is applied to common midpoint stacked data. One of the main challenges within the tomographic approach is the reliable estimation of the kinematic attributes from the data that are used in the inversion process. Estimations in the prestack domain (weak and noisy signals), as well as in the post‐stack domain (occurrence of triplications and diffractions leading to numerous conflicting dip situations) may lead to parameter inaccuracies that will adversely impact the resulting velocity models. To overcome the above limitations, a new tomographic procedure applied in the time‐migrated domain is proposed. We call this method Image‐Incident‐Point‐wave tomography. The new scheme can be seen as an alternative to Normal‐Incidence‐Point‐wave tomography. The latter method is based on traveltime attributes associated with normal rays, whereas the Image‐Incidence‐Point‐wave technique is based on the corresponding quantities for the image rays. Compared to Normal‐Incidence‐Point‐wave tomography the proposed method eases the selection of the tomography attributes, which is shown by synthetic and field data examples. Moreover, the method provides a direct way to convert time‐migration velocities into depth‐migration velocities without the need of any Dix‐style inversion. 相似文献
4.
Soazig Le Bégat Hervé Chauris † Vincent Devaux Sylvain Nguyen Mark Noble 《Geophysical Prospecting》2004,52(5):427-438
The estimation of velocity models is still crucial in seismic reflection imaging as it controls the quality of the depth‐migrated image, which is the basis of geological interpretation. Among the numerous existing methods for velocity determination, tomographic methods are very attractive for their efficiency and ability to retrieve heterogeneities of the medium. We present three tomographic methods in order to estimate heterogeneous velocity models from 2D prestack PP reflection data: a traveltime tomography in the time‐migrated domain, a traveltime and slope tomography in the non‐migrated time domain, and a slope tomography in the depth‐migrated domain. The first method (traveltime tomography in the time domain) is based on continuous picked events, whereas the two slope tomographic methods, one in the time domain and the other in the depth domain, are based on locally coherent events, with no assumptions about reflector geometry or the unknown velocity field. The purpose of this paper is not to describe in detail the theoretical basis and implementation of the methods, but to apply and compare their output using the same marine real data set. Based on the estimated velocity models, the migrated images and the common‐image gathers from the three processing routes, the relative strengths and weaknesses of the methods are discussed. Finally, similarities are indicated and potential alternative approaches are proposed. 相似文献
5.
Sylvain Nguyen † Reda Baina Mathias Alerini ‡ Gilles Lambaré Vincent Devaux Mark Noble 《Geophysical Prospecting》2008,56(5):613-625
Depth velocity model building remains a difficult step within the seismic depth imaging sequence. Stereotomography provides an efficient solution to this problem but was limited until now to a picking of seismic data in the prestack time un-migrated domain. We propose here a method for stereotomographic data picking in the depth migrated domain. Picking in the depth migrated domain exhibits the advantage of a better signal-to-noise ratio and of a more regular distribution of picked events in the model, leading to a better constrained tomographic inverse problem. Moreover, any improvement on the velocity model will improve the migrated results, again leading to improved picking. Our strategy for obtaining a stereotomographic dataset from a prestack depth migration is based on migration of attributes (and not on a kinematic demigration approach!). For any locally coherent event in the migrated image, migration of attributes allows one to compute ray parameter attributes corresponding to the specular reflection angle and dip. For application to stereotomography, the necessary attributes are the source/receiver locations, the traveltime and the data slopes. For the data slope, when the migration velocity model is erroneous, some additional corrections have to be applied to the result of migration of the attributes. Applying these corrections, our picking method is theoretically valid whatever the quality of the migration velocity model. We first present the theoretical aspects of the method and then validate it on 2D synthetic and real seismic reflection data sets. 相似文献
6.
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. 相似文献
7.
8.
Full‐waveform inversion is re‐emerging as a powerful data‐fitting procedure for quantitative seismic imaging of the subsurface from wide‐azimuth seismic data. This method is suitable to build high‐resolution velocity models provided that the targeted area is sampled by both diving waves and reflected waves. However, the conventional formulation of full‐waveform inversion prevents the reconstruction of the small wavenumber components of the velocity model when the subsurface is sampled by reflected waves only. This typically occurs as the depth becomes significant with respect to the length of the receiver array. This study first aims to highlight the limits of the conventional form of full‐waveform inversion when applied to seismic reflection data, through a simple canonical example of seismic imaging and to propose a new inversion workflow that overcomes these limitations. The governing idea is to decompose the subsurface model as a background part, which we seek to update and a singular part that corresponds to some prior knowledge of the reflectivity. Forcing this scale uncoupling in the full‐waveform inversion formalism brings out the transmitted wavepaths that connect the sources and receivers to the reflectors in the sensitivity kernel of the full‐waveform inversion, which is otherwise dominated by the migration impulse responses formed by the correlation of the downgoing direct wavefields coming from the shot and receiver positions. This transmission regime makes full‐waveform inversion amenable to the update of the long‐to‐intermediate wavelengths of the background model from the wide scattering‐angle information. However, we show that this prior knowledge of the reflectivity does not prevent the use of a suitable misfit measurement based on cross‐correlation, to avoid cycle‐skipping issues as well as a suitable inversion domain as the pseudo‐depth domain that allows us to preserve the invariant property of the zero‐offset time. This latter feature is useful to avoid updating the reflectivity information at each non‐linear iteration of the full‐waveform inversion, hence considerably reducing the computational cost of the entire workflow. Prior information of the reflectivity in the full‐waveform inversion formalism, a robust misfit function that prevents cycle‐skipping issues and a suitable inversion domain that preserves the seismic invariant are the three key ingredients that should ensure well‐posedness and computational efficiency of full‐waveform inversion algorithms for seismic reflection data. 相似文献
9.
本文通过选取合适的叠前时间偏移软件，对两块三维地震资料进行偏移成像试验，验证叠前时间偏移中影响偏移成像效果的几个主要因素.该软件偏移算法的核心技术是弯曲射线偏移处理，这不同于工业界常用的直射线假设.偏移速度是偏移成像好坏的主要因素，通过迭代进行偏移、速度分析，使共成像点道集拉平，从而实现构造的准确成像；偏移孔径也是影响偏移成像的一个关键参数，其选取与成像目标层的倾斜角、深度、速度等有关；反假频参数对偏移成像效果有一定影响，是偏移中需要考虑的因素之一. 相似文献
10.
描述地震波衰减特征的品质因子Q对地震数据处理和油藏描述非常重要,在地震勘探领域,Q值一般通过垂直地震剖面(VSP)数据或地面地震数据得到.由于叠前地面地震数据具有复杂的射线路径且存在噪声、调谐干涉效应等影响,从叠前地震数据中准确估计Q值相对困难.本文以地震波射线传播为基础,根据同相轴局部斜率和射线参数的映射关系,将多射线波形频谱同时带入谱比法联合反演估计Q值,提出了基于多射线联合反演的速度无关叠前Q值估计方法.该方法通过局部斜率属性避开了速度对Q值估计的影响,局部斜率携带地震波传播的速度信息,具有相同局部斜率的地震反射波具有相同的传播射线参数.同相轴局部斜率是地震数据域的属性,而速度是模型域的参数,在估计Q值中采用数据域的属性参数可以直接应用于数据的联合反演,而不需要通过速度对其做进一步的转化,从而提高了Q值估计的精度.同时,本方法采用预测映射(predictive mapping)技术将非零炮检距反射信息映射到零炮检距处,从而获得零偏移距走时对应的Q值.模拟和实际算例验证了本文方法的有效性. 相似文献
11.
速度是地震偏移成像准确与否的关键所在.全波形反演综合利用地震波场运动学和动力学信息,能够得到相比传统速度建模方法更高频的成分.全波形反演的理论比较成熟,但实际应用成功的例子相对较少,特别是对于陆上地震资料.塔里木盆地地震地质条件复杂,为了实现缝洞型储层的准确成像,本文开展了针对目标靶区的全波形反演精细速度建场研究.采用一种时间域分层多尺度全波形反演流程:首先通过层析成像建立初始速度模型;其次利用折射波反演浅层速度模型;最后利用反射波反演中深层速度模型.偏移成像结果表明基于全波形反演的速度建模技术能有效改善火成岩下伏构造的成像精度,显示了全波形反演在常规陆上采集资料的应用潜力. 相似文献
12.
与其他偏移方法相比,逆时偏移基于精确的波动方程而不是对其近似,用时间外推来代替深度外推.因此,它具有良好的精度,不受地下构造倾角和介质横向速度变化的限制.激发时间成像条件的求取是叠前逆时偏移的难点之一,本文采用求解程函方程的方法得到地下各点的初至波走时,以此作为叠前逆时偏移的成像条件.基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法精度较高并适用于强纵横向变速的复杂介质.试算结果表明,在复杂介质模型中利用叠前逆时深度偏移收到了很好的成像效果. 相似文献
13.
In this study we present the workflow and results of 2D frequency domain waveform tomography applied to the global‐offset seismic data acquired in central Poland along a 50‐km long profile during the GRUNDY 2003 experiment. The waveform tomography method allows full exploitation of the wide‐aperture content of these data and produces in a semi‐automatic way both the detailed P‐wave velocity model and the structural image (i.e., perturbations in respect to the starting model). Thirteen frequencies ranging from 4 to 16 Hz were inverted sequentially, gradually introducing higher wavenumbers and more details into the velocity models. Although the data were characterised by relatively large shot spacings (1.5 km), we obtained clear images both of the Mesozoic and Permian sedimentary cover. Velocity patterns indicated facies changes within the Jurassic and Zechstein strata. A high velocity layer (ca. 5500 m/s) was found near the base of Triassic (Scythian), which made the imaging of a deeper layer difficult. Nevertheless, we were able to delineate the base of the Permian (i.e., base of the Rotliegend), which was not possible to derive from conventional common‐depth‐point processing, as well as some deeper events, attributed to the Carboniferous. The sub‐Permian events formed a syn‐form which favoured our previous interpretation of a depression filled with Upper Carboniferous molasse. The validity of the waveform tomography‐derived model was confirmed by well‐log data. Forward ray‐tracing modelling and synthetic seismograms calculations provided another justification for the key structures present in the waveform tomography model. 相似文献
14.
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. 相似文献
15.
偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势. 相似文献
16.
Wave-equation migration velocity analysis. I. Theory 总被引：2，自引：0，他引：2
We present a migration velocity analysis (MVA) method based on wavefield extrapolation. Similarly to conventional MVA, our method aims at iteratively improving the quality of the migrated image, as measured by the flatness of angle‐domain common‐image gathers (ADCIGs) over the aperture‐angle axis. However, instead of inverting the depth errors measured in ADCIGs using ray‐based tomography, we invert ‘image perturbations’ using a linearized wave‐equation operator. This operator relates perturbations of the migrated image to perturbations of the migration velocity. We use prestack Stolt residual migration to define the image perturbations that maximize the focusing and flatness of ADCIGs. Our linearized operator relates slowness perturbations to image perturbations, based on a truncation of the Born scattering series to the first‐order term. To avoid divergence of the inversion procedure when the velocity perturbations are too large for Born linearization of the wave equation, we do not invert directly the image perturbations obtained by residual migration, but a linearized version of the image perturbations. The linearized image perturbations are computed by a linearized prestack residual migration operator applied to the background image. We use numerical examples to illustrate how the backprojection of the linearized image perturbations, i.e. the gradient of our objective function, is well behaved, even in cases when backprojection of the original image perturbations would mislead the inversion and take it in the wrong direction. We demonstrate with simple synthetic examples that our method converges even when the initial velocity model is far from correct. In a companion paper, we illustrate the full potential of our method for estimating velocity anomalies under complex salt bodies. 相似文献
17.
Double-square-root one-way wave equation prestack tau migration in heterogeneous media 总被引：1，自引：0，他引：1
In this paper, source‐receiver migration based on the double‐square‐root one‐way wave equation is modified to operate in the two‐way vertical traveltime (τ) domain. This tau migration method includes reasonable treatment for media with lateral inhomogeneity. It is implemented by recursive wavefield extrapolation with a frequency‐wavenumber domain phase shift in a constant background medium, followed by a phase correction in the frequency‐space domain, which accommodates moderate lateral velocity variations. More advanced τ‐domain double‐square‐root wave propagators have been conceptually discussed in this paper for migration in media with stronger lateral velocity variations. To address the problems that the full 3D double‐square‐root equation prestack tau migration could meet in practical applications, we present a method for downward continuing common‐azimuth data, which is based on a stationary‐phase approximation of the full 3D migration operator in the theoretical frame of prestack tau migration of cross‐line constant offset data. Migrations of synthetic data sets show that our tau migration approach has good performance in strong contrast media. The real data example demonstrates that common‐azimuth prestack tau migration has improved the delineation of the geological structures and stratigraphic configurations in a complex fault area. Prestack tau migration has some inherent robust characteristics usually associated with prestack time migration. It follows a velocity‐independent anti‐aliasing criterion that generally leads to reduction of the computation cost for typical vertical velocity variations. Moreover, this τ‐domain source‐receiver migration method has features that could be of help to speed up the convergence of the velocity estimation. 相似文献
18.
The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. VP, VS and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data. 相似文献
19.
地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性. 相似文献
20.
立体层析成像是一种新的地震反射波层析成像方法,能为叠前深度偏移提供较为精确的宏观速度模型。本文研究了立体层析成像的实现方法,包括斜率与走时数据的拾取、离散速度模型构建和初始化、射线参数的确定、斜率和走时及射线计算以及反演问题解法等,建立立体层析成像的算法流程。并通过对Marmousi模型试验,对立体层析成像运行所需的主要参数,如初始速度模型、拾取数据量、离散网格尺寸、速度平滑权重等进行测试和分析,总结这些不同参数对立体层析反演结果的影响规律,用以指导生产实践。 相似文献