共查询到20条相似文献，搜索用时 31 毫秒
1.
In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity. 相似文献
2.
In many land seismic situations, the complex seismic wave propagation effects in the near‐surface area, due to its unconsolidated character, deteriorate the image quality. Although several methods have been proposed to address this problem, the negative impact of 3D complex near‐surface structures is still unsolved to a large extent. This paper presents a complete 3D data‐driven solution for the near‐surface problem based on 3D one‐way traveltime operators, which extends our previous attempts that were limited to a 2D situation. Our solution is composed of four steps: 1) seismic wave propagation from the surface to a suitable datum reflector is described by parametrized one‐way propagation operators, with all the parameters estimated by a new genetic algorithm, the self‐adjustable input genetic algorithm, in an automatic and purely data‐driven way; 2) surface‐consistent residual static corrections are estimated to accommodate the fast variations in the near‐surface area; 3) a replacement velocity model based on the traveltime operators in the good data area (without the near‐surface problem) is estimated; 4) data interpolation and surface layer replacement based on the estimated traveltime operators and the replacement velocity model are carried out in an interweaved manner in order to both remove the near‐surface imprints in the original data and keep the valuable geological information above the datum. Our method is demonstrated on a subset of a 3D field data set from the Middle East yielding encouraging results. 相似文献
3.
4.
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. 相似文献
5.
Land seismic data quality can be severely affected by near‐surface anomalies. The imprint of a complex near‐surface can be removed by redatuming the data to a level below the surface, from where the subsurface structures are assumed to be relatively smooth. However, to derive a velocity‐depth model that explains the propagation effects of the near‐surface is a non‐trivial task. Therefore, an alternative approach has been proposed, where the redatuming operators are obtained in a data‐driven manner from the reflection event related to the datum. In the current implementation, the estimation of these redatuming operators is done in terms of traveltimes only, based on a high‐frequency approximation. The accompanying amplitudes are usually derived from a local homogeneous medium, which is obviously a simplification of reality. Such parametrization has produced encouraging results in the past but cannot completely remove the near‐surface complexities, leaving artefacts in the redatumed results. In this paper we propose a method that estimates the redatuming operators directly from the data, i.e., without using a velocity model, in a full waveform manner, such that detailed amplitude and phase variations are included. The method directly outputs the inverse propagation operators that are needed for true‐amplitude redatuming. Based on 2D synthetic data it is demonstrated that the resulting redatuming quality is improved and artefacts are reduced. 相似文献
6.
Quantitative imaging of complex structures from dense wide-aperture seismic data by multiscale traveltime and waveform inversions: a case study 总被引：3，自引：0，他引：3
S. Operto C. Ravaut L. Improta J. Virieux A. Herrero P. Dell'Aversana 《Geophysical Prospecting》2004,52(6):625-651
An integrated multiscale seismic imaging flow is applied to dense onshore wide‐aperture seismic data recorded in a complex geological setting (thrust belt). An initial P‐wave velocity macromodel is first developed by first‐arrival traveltime tomography. This model is used as an initial guess for subsequent full‐waveform tomography, which leads to greatly improved spatial resolution of the P‐wave velocity model. However, the application of full‐waveform tomography to the high‐frequency part of the source bandwidth is difficult, due to the non‐linearity of this kind of method. Moreover, it is computationally expensive at high frequencies since a finite‐difference method is used to model the wave propagation. Hence, full‐waveform tomography was complemented by asymptotic prestack depth migration to process the full‐source bandwidth and develop a sharp image of the short wavelengths. The final traveltime tomography model and two smoothed versions of the final full‐waveform tomography model were used as a macromodel for the prestack depth migration. In this study, wide‐aperture multifold seismic data are used. After specific preprocessing of the data, 16 frequency components ranging from 5.4 Hz to 20 Hz were inverted in cascade by the full‐waveform tomography algorithm. The full‐waveform tomography successfully imaged SW‐dipping structures previously identified as high‐resistivity bodies. The relevance of the full‐waveform tomography models is demonstrated locally by comparison with a coincident vertical seismic profiling (VSP) log available on the profile. The prestack depth‐migrated images, inferred from the traveltime, and the smoothed full‐waveform tomography macromodels are shown to be, on the whole, consistent with the final full‐waveform tomography model. A more detailed analysis, based on common‐image gather computations, and local comparison with the VSP log revealed that the most accurate migrated sections are those obtained from the full‐waveform tomography macromodels. A resolution analysis suggests that the asymptotic prestack depth migration successfully migrated the wide‐aperture components of the data, allowing medium wavelengths in addition to the short wavelengths of the structure to be imaged. The processing flow that we applied to dense wide‐aperture seismic data is shown to provide a promising approach, complementary to more classical seismic reflection data processing, to quantitative imaging of complex geological structures. 相似文献
7.
地下复杂介质地震处理中的CFP技术 总被引：1，自引：3，他引：1
要简要介绍CFP(Common Focus Point)方法技术的基本原理和主要应用。CFP是复杂介质地震处理中的一项新技术，它把叠前偏移分成两个独立的步骤：首先对检波点(炮点)进行聚焦处理，产生共聚焦点道集(CFP道集)，然后再对炮点(检波点)进行聚焦，产生叠前偏移的输出。两个步骤中间的CFP道集则可以进行其它处理。如果是为了寻找构造信息，那么可以应用共焦点CFP偏移；如果是为了寻找岩石、孔隙或流体的信息，则要应用双焦点CFP偏移。目前，该技术主要应用于：(1)CFP两步聚焦法偏移；(2)叠前深度偏移速度模型的建立；(3)试图通过算于而不是速度来解决复杂地表的静校正；(4)消除全程或层间多次波；(5)CFP方法基准面的延拓和盐下成像；(6)多分量地震资料的偏移成像。 相似文献
8.
While seismic imaging for crustal and mantle structures has traditionally relied on surface wave and refraction data, the use of reflection data for crustal-scale targets has been largely limited to the common midpoint (CMP) stack techniques. The rapid increase in the number of seismograph array deployments in recent years in crustal and mantle seismology has reached a level such that a re-examination of the imaging techniques is becoming necessary. In this paper we show the advantage of prestack depth imaging for crustal reflection studies, based on data from two reflection surveys of the Los Angeles Regional Seismic Experiment (LARSE) to map faults and crustal-scale structures. Our analysis indicates that the quality of the previous images of these surveys is limited by the CMP stack technique. For comparison, we present here depth images of the same LARSE data using wave equation prestack depth imaging and a tomographic velocity model based on first arrivals of the LARSE surveys and local earthquakes. Our new images are considerably improved over previous images in terms of resolution and reflector continuity. The new images show reflectors throughout the crust and suggest truncations in the Moho associated with the San Andreas Fault. A series of bright reflector segments, which are associated with the San Gabriel and San Andreas faults have been identified and might represent reflections from the fault zones. Our results suggest that the presence of high noise level, strong lateral velocity heterogeneity and wide angle geometry argue for, rather than against, the use of prestack depth imaging over the simple CMP stack techniques. As demonstrated in this study, it is now viable to conduct prestack depth imaging of crustal reflection data using a velocity model based on earthquake first arrivals thanks to the dense acquisition deployment. 相似文献
9.
A single‐layer model of the near surface throughout the Kingdom of Saudi Arabia is available. While this simple model suffices for most areas and large subsurface structures, it fails in situations where the surface topography is complex, the base of weathering is below the datum, or where the time structural closure is less than the uncertainty in the static correction. In such cases, multiple‐layered models that incorporate velocities derived from analysis of first arrivals picked from seismic shot records have proved to be successful in defining the lateral heterogeneity of the near surface. The additional velocity information obtained from this first‐arrival analysis (direct as well as refracted arrivals) vastly improves the velocity–depth model of the near surface, regardless of the topography. Static corrections computed from these detailed near‐surface velocity models have significantly enhanced subsurface image focusing, thereby reducing the uncertainty in the closure of target structures. Other non‐seismic methods have been used either to confirm qualitatively or to enhance the layer models previously mentioned. Gravity data may be particularly useful in sandy areas to confirm general structure, while geostatistical modelling of vibrator base‐plate attributes has yielded information that enhances the velocity field. In the global context, exploration targets of the oil and gas industry are seeking smaller and lower relief‐time structures. Thus, near‐surface models will need to enhance and integrate these methods, particularly in areas where the assumption of flat‐lying near‐surface layers cannot be met. 相似文献
10.
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. 相似文献
11.
Wavefield extrapolation and prestack depth migration in anelastic inhomogeneous media 总被引：6，自引：0，他引：6
Wavefield depth extrapolation and prestack depth migration in complex anelastic media are studied. Kjartansson's frequency‐independent Q law is used to describe the absorption of seismic energy. The macromodel used is analogous to the macromodel used for current migration schemes except that an additional frequency‐independent Q macromodel needs to be provided. Absorption in the forward one‐way propagator is introduced by assuming a complex phase velocity, and the inverse one‐way propagator is obtained using the reciprocity theorem for one‐way wavefields in dissipative media. The stability of the inverse propagator is achieved by limiting the angle of propagation of wavefields. A table‐driven explicit operator scheme for imaging complex 2D anelastic media is presented. High‐accuracy, short convolution operators are designed by the weighted least‐squares method, and two kinds of imaging conditions are proposed. Numerical examples of depth extrapolation in laterally varying media, the migration of a spatial impulse with dispersion as well as shot record depth migration demonstrate the potential of the proposed explicit forward operator, the explicit inverse operator and the prestack depth migration scheme, respectively. 相似文献
12.
复杂构造地震成像主要取决于叠前地震数据品质、偏移速度可靠性和偏移算子成像精度. 库车坳陷异常复杂的近地表条件导致极低信噪比的地震采集数据. 该区逆冲推覆高陡构造刺穿盐体大面积分布, 盐层厚度变化大、顶底面形态复杂, 盐下断裂带破碎、小断块发育, 形成异常复杂的地震成像问题. 本文重点研究三个关键环节:(1)精细的叠前地震预处理研究: 根据该区地震地质复杂性和地震资料特征, 采用一些新的方法技术和技术组合从振幅与时移的大、中、小尺度变化三个层次来解决资料信噪比问题, 重建深部反射信号; (2)三级偏移速度分析研究:利用库车坳陷盐刺穿逆冲推覆构造建模理论及变速成图配套技术解决叠前时间偏移速度场时深转换问题,利用井约束低频速度地震迭代反演技术解决连井层速度场与偏移速度场的融合问题,实现从DMO速度分析、叠前时间偏移速度分析到叠前深度偏移速度分析的有机衔接,建立拓扑结构相对保持的叠前深度偏移速度模型;(3)基于退化Fourier偏移算子的半解析波动方程叠前时间和深度偏移研究, 极大地改善了地震偏移过程中高波数波的成像问题. 通过对库车坳陷大北、博孜、却勒、西秋4和西秋10等复杂高陡构造的叠前时间和深度偏移地震成像处理,取得了较好的应用效果. 相似文献
13.
Shear‐wave statics in marine seismic exploration data are routinely too large to be estimated using conventional techniques. Near‐surface unconsolidated sediments are often characterized by low values of Vs and steep velocity gradients. Minor variations in sediment properties at these depths correspond to variations in the shear‐wave velocity and will produce significant static shifts. It is suggested that a significant proportion of the shear‐wave statics solution can be estimated by performing a separate high‐resolution survey to target near‐surface unconsolidated sediments. Love‐wave, shear‐wave refraction and geotechnical measurements were individually used to form high‐resolution near‐surface shear‐wave velocity models to estimate the shear‐wave statics for a designated survey line. Comparisons with predicted statics revealed that shear‐wave statics could not be estimated using a velocity model predicted by substituting geotechnical measurements into empirical relationships. Empirical relationships represent a vast simplification of the factors that control Vs and are therefore not sufficiently sensitive to estimate shear‐wave statics. Refraction measurements are potentially sensitive to short‐wavelength variations in sediment properties when combined with accurate navigational data. Statics estimated from Love‐wave data are less sensitive, and sometimes smoothed in appearance, since interpreted velocity values represent an average both laterally and vertically over the receiver array and the frequency–depth sensitivity range, respectively. For the survey site, statics estimated from near‐surface irregularities using shear‐wave refraction measurements represent almost half the total statics solution. More often, this proportion will be greater when bedrock relief is less. 相似文献
14.
针对鄂尔多斯盆地西缘黄土塬区复杂地表和复杂地下构造导致难以准确成像问题,采用浅层潜水波层析反演（DWT）速度建模技术,同时辅以中深层反射波层析成像技术,形成一套实用叠前深度偏移速度建模方法。首先生成一个基于钻井和解释信息的起始近地表速度,其次利用潜水波层析反演建立初始近地表模型,将其与常规处理获得的中深层速度模型进行匹配拼接,建立起初始的起伏地表全速度模型,然后在此基础上利用基于反射波的网格层析进行中深层速度建模,经过多轮次迭代,最终获得可靠的高精度速度模型。鄂尔多斯盆地西缘MJT工区地震资料的成像处理验证了这一套速度建模技术的有效性,地下构造成像更合理也更精确。 相似文献
15.
Ralph Bridle Nikolai Barsoukov Mohammad Al-Homaili Robert Ley Ameerah Al-Mustafa 《Geophysical Prospecting》2006,54(6):667-680
We present a seismic Test Line, provided by Saudi Aramco for various research teams, to highlight a few major challenges in land data processing due to near‐surface anomalies. We discuss state‐of‐the‐art methods used to compensate for shallow distortions, including single‐layer, multilayer, plus/minus, refraction and tomostatics methods. They are a starting point for the new technologies presented in other papers, all dealing with the same challenging data described here. The difficulties on the Test Line are mostly due to the assumption of vertical raypaths, inherent in classical applications of near‐surface correction statics. Even the most detailed velocity/depth model presents difficulties, due to the compleX‐raypath. There is a need for methods which are based on more complex models andtheories. 相似文献
16.
Yonghai Zhang 《Geophysical Prospecting》2010,58(6):915-924
Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface. 相似文献
17.
复杂地表探区,尤其是盆山过渡区的油气勘探是我国也是世界上油气勘探的重点区域,但是此类区域油气地震勘探中满足精确地震成像的速度建模一直是个没有很好解决的问题.本文提出了一种综合性的数据域初至波走时与成像域反射波走时联合层析复杂地表浅中深层速度建模方法,并针对联合层析速度反演解的非唯一性问题,深入地分析了层析反演中正则化的本质意义,指出了建立构造特征正则化方法的具体技术路线,提出了联合层析的实现流程及策略.理论和实际数据试验表明,本文提出的数据域初至走时与成像域反射走时联合层析浅中深层速度建模技术避免了常规建模方法中浅层速度模型与中深层速度模型的融合问题,较好地解决了传统成像域反射层析对近地表模型的不可控更新问题,整体提升了深度域浅中深层速度模型的建模精度,进而提高了复杂地表、复杂构造区的地震成像质量. 相似文献
18.
19.
Survey sinking migration downward continues the entire surface observed multi‐shot data to the subsurface step by step recursively. Reflected energy from reflectors at current depth appear at zero time and zero offset in the extrapolated wavefield. The data (seismic records) of t > 0 at this depth are equivalent to the data acquired by a survey system deployed at this depth. This is the reason to name the process ‘survey sinking’. The records of negative time need not to be further propagated since they carry no information to image structures beneath the new survey system. In this paper, we combine survey sinking with dreamlet migration. The dreamlet migration method decomposes the seismic wavefield and one‐way wave propagator by complete time‐space localized bases. The localization on time gives flexibility on time‐varying operations during depth extrapolation. In dreamlet survey sinking migration, it only keeps the data for imaging the structures beneath the sunk survey system and gets rid of the data already used to image structures above it. The deeper the depth is, the shorter is the valid time records of the remaining data and less computation is needed for one depth step continuation. For data decomposition, in addition to time axis, dreamlet survey sinking also decomposes the data for source and receiver gathers, which is a fully localized decomposition of prestack seismic data. A three‐scatter model is first used to demonstrate the computational feature and principle of this method. Tests on the two‐dimensional SEG/EAGE salt model show that with reduced data sets the proposed method can still obtain good imaging quality on complex geology structures and a strong velocity contrast environment. 相似文献
20.
研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献