共查询到20条相似文献，搜索用时 15 毫秒
1.
苏北大陆科学钻探靶区深反射地震的叠前深度偏移 总被引：4，自引：2，他引：2
由于深反射地震数据具有信噪比低和记录长度长等特点,叠前深度偏移方法的应用有许多困难．为此,我们研究了一种适合于深反射地震的叠前深度偏移方法;包括：用逆风有限差分方法计算程函方程;在常规速度扫描的基础上,用协方差控制提高速度分析精度;用联合反演算法计算层速度,再插值后得到初始速度模型;用Ｋｉｒｃｈｈｏｆｆ法作为偏移速度分析工具,求得最终的速度模型;最终的速度模型作为有限差分深度偏移的输入,求得最终的偏移结果．用该方法对“中国大陆科学深钻工程”东海二维深反射地震数据ＤＨ－４线进行了叠前深度偏移,取得了良好的效果。 相似文献
2.
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. 相似文献
3.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
4.
基于炮集数据的常规叠前深度偏移波场外推的计算量非常大,控制照明叠前深度偏移是基于快速的组合炮技术,通过在目的层上控制震源波场,高效实现对地质目标的高质量成像. 组合炮所需的合成算子,由在目的层定义的震源波场反向外推到地表建立. 本文提出一种更快速的地表旋转控制照明方法,直接在地表对震源波场进行面向目标的照明控制,实现目标区域高质量成像. 该方法只需对某一核合成算子进行简单的几何旋转建立合成算子函数,避免了复杂的波场外推. 这种新的控制照明方法与快速的一般组合炮偏移效率相同,但对复杂目标区域的成像能力显著增强. 应用本文方法对Marmousi模型进行计算,得到深层复杂构造非常理想的成像结果. 相似文献
5.
基于单程波深度延拓算法,得到震源波场在成像点的入射角度,结合对地层倾角的估计,获得入射地震波与界面法线的夹角.通过运用"保幅"的反褶积成像条件和考虑累加的炮数,解决了炮点覆盖不均匀导致的成像幅值误差问题,进而建议了炮域波动方程叠前深度偏移直接产生角道集的方法和流程.与基于空间移动或时移成像条件的波动方程叠前深度偏移提取角道集的方法相比,本文建议的方法只需少量额外的存储空间,又可补偿观测系统非均匀覆盖对成像幅值的影响;其增加的计算量与炮域偏移算法相比几乎可以忽略.文中算例表明,本文方法提取的角道集可为叠前反演提供较精确的AVO振幅特性.此外,就改善地震成像效果本身而言,提取角道集使得可在波动方程叠前深度偏移中应用剩余动校和拉伸切除技术,从而可更好地保持高频成分并提高成像的信噪比. 相似文献
6.
Diffraction imaging in depth 总被引：3，自引：0，他引：3
High resolution imaging is of great value to an interpreter, for instance to enable identification of small scale faults, and to locate formation pinch-out positions. Standard approaches to obtain high-resolution information, such as coherency analysis and structure-oriented filters, derive attributes from stacked, migrated images. Since they are image-driven, these techniques are sensitive to artifacts due to an inadequate migration velocity; in fact the attribute derivation is not based on the physics of wave propagation. Diffracted waves on the other hand have been recognized as physically reliable carriers of high- or even super-resolution structural information. However, high-resolution information, encoded in diffractions, is generally lost during the conventional processing sequence, indeed migration kernels in current migration algorithms are biased against diffractions. We propose here methods for a diffraction-based, data-oriented approach to image resolution. We also demonstrate the different behaviour of diffractions compared to specular reflections and how this can be leveraged to assess characteristics of subsurface features. In this way a rough surface such as a fault plane or unconformity may be distinguishable on a diffraction image and not on a traditional reflection image.
We outline some characteristic properties of diffractions and diffraction imaging, and present two novel approaches to diffraction imaging in the depth domain. The first technique is based on reflection focusing in the depth domain and subsequent filtering of reflections from prestack data. The second technique modifies the migration kernel and consists of a reverse application of stationary-phase migration to suppress contributions from specular reflections to the diffraction image. Both techniques are proposed as a complement to conventional full-wave pre-stack depth migration, and both assume the existence of an accurate migration velocity. 相似文献
We outline some characteristic properties of diffractions and diffraction imaging, and present two novel approaches to diffraction imaging in the depth domain. The first technique is based on reflection focusing in the depth domain and subsequent filtering of reflections from prestack data. The second technique modifies the migration kernel and consists of a reverse application of stationary-phase migration to suppress contributions from specular reflections to the diffraction image. Both techniques are proposed as a complement to conventional full-wave pre-stack depth migration, and both assume the existence of an accurate migration velocity. 相似文献
7.
8.
9.
三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节．由于其技术实现不仅涉及波场延拓理论的创新,而且需要大规模计算,因而研究难度较大．本文以实验效果的取得为目的,完整地实现了SEG／EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算．文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力,具体包括：1)对于盐丘模型,文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较;2)对推覆体模型,文中进行了脉冲响应测试;3)由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力。 相似文献
10.
Tariq Alkhalifah 《Geophysical Prospecting》2005,53(5):643-653
Imaging pre‐salt reflections for data acquired from the coastal region of the Red Sea is a task that requires prestack migration velocity analysis. Conventional post‐stack time processing lacks the lateral inhomogeneity capability, necessary for such a problem. Prestack migration velocity analysis in the vertical time domain reduces the velocity–depth ambiguity that usually hampers the performance of prestack depth‐migration velocity analysis. In prestack τ‐migration velocity analysis, the interval velocity model and the output images are defined in τ (i.e. vertical time). As a result, we avoid placing reflectors at erroneous depths during the velocity analysis process and thus avoid inaccurately altering the shape of the velocity model, which in turn speeds up the convergence to the true model. Using a 1D velocity update scheme, the prestack τ‐migration velocity analysis produces good images of data from the Midyan region of the Red Sea. For the first seismic line from this region, only three prestack τ‐migration velocity analysis iterations were required to focus pre‐salt reflections in τ. However, the second line, which crosses the first line, is slightly more complicated and thus required five iterations to reach the final, reasonably focused, τ‐image. After mapping the images for the two crossing lines to depth, using the final velocity models, the placements of reflectors in the two 2D lines were consistent at their crossing point. Some errors occurred due to the influence of out‐of‐plane reflections on 2D imaging. However, such errors are identifiable and are generally small. 相似文献
11.
盐丘模型 推覆体模型 波动方程 叠前深度偏移 总被引：4，自引：0，他引：4
三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节.由于其技术实现不仅涉及波场延拓理论的创新，而且需要大规模计算，因而研究难度较大. 本文以实验效果的取得为目的，完整地实现了SEG/EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算. 文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力. 具体包括：1）对于盐丘模型，文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较；2）对推覆体模型，文中进行了脉冲响应测试；3）由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力. 相似文献
12.
Kai Yang Bao‐shu Chen Xiao‐jiang Wang Xiao‐chun Yang Jun‐rong Liu 《Geophysical Prospecting》2012,60(2):255-269
In the application of a conventional common‐reflection‐surface (CRS) stack, it is well‐known that only one optimum stacking operator is determined for each zero‐offset sample to be simulated. As a result, the conflicting dip situations are not taken into account and only the most prominent event contributes to any a particular stack sample. In this paper, we name this phenomenon caused by conflicting dip problems as ‘dip discrimination phenomenon’. This phenomenon is not welcome because it not only leads to the loss of weak reflections and tips of diffractions in the final zero‐offset‐CRS stacked section but also to a deteriorated quality in subsequent migration. The common‐reflection‐surface stack with the output imaging scheme (CRS‐OIS) is a novel technique to implement a CRS stack based on a unified Kirchhoff imaging approach. As far as dealing with conflicting dip problems is concerned, the CRS‐OIS is a better option than a conventional CRS stack. However, we think the CRS‐OIS can do more in this aspect. In this paper, we propose a workflow to handle the dip discrimination phenomenon based on a cascaded implementation of prestack time migration, CRS‐OIS and prestack time demigration. Firstly, a common offset prestack time migration is implemented. Then, a CRS‐OIS is applied to the time‐migrated common offset gather. Afterwards, a prestack time demigration is performed to reconstruct each unmigrated common offset gather with its reflections being greatly enhanced and diffractions being well preserved. Compared with existing techniques dealing with conflicting dip problems, the technique presented in this paper preserves most of the diffractions and accounts for reflections from all possible dips properly. More importantly, both the post‐stacked data set and prestacked data set can be of much better quality after the implementation of the presented scheme. It serves as a promising alternative to other techniques except that it cannot provide the typical CRS wavefield attributes. The numerical tests on a synthetic Marmousi data set and a real 2D marine data set demonstrated its effectiveness and robustness. 相似文献
13.
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. 相似文献
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.
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. 相似文献
16.
Seismic imaging is an important step for imaging the subsurface structures of the Earth. One of the attractive domains for seismic imaging is explicit frequency–space (f – x) prestack depth migration. So far, this domain focused on migrating seismic data in acoustic media, but very little work assumed visco‐acoustic media. In reality, seismic exploration data amplitudes suffer from attenuation. To tackle the problem of attenuation, new operators are required, which compensates for it. We propose the weighted L 1 ‐error minimisation technique to design visco‐acoustic f – x wavefield extrapolators. The L 1 ‐error wavenumber responses provide superior extrapolator designs as compared with the previously designed equiripple L 4 ‐norm and L∞‐norm extrapolation wavenumber responses. To verify the new compensating designs, prestack depth migration is performed on the challenging Marmousi model dataset. A reference migrated section is obtained using non‐compensating f – x extrapolators on an acoustic dataset. Then, both compensating and non‐compensating extrapolators are applied to a visco‐acoustic dataset, and both migrated sections are then compared. The final images show that the proposed weighted L 1 ‐error method enhances the resolution and results in practically stable images. 相似文献
17.
不同于常规广义屏传播算子的推导中使用散射理论,本文利用单平方根算子的渐近展开,推导出了单程波方程广义屏传播算子的高阶表达式.高阶广义屏传播算子不仅可提高常规广义屏传播算子的计算精度,而且还能改善广义屏传播算子对速度强横向变化介质的适应性.把高阶广义屏传播算子应用于波动方程叠前深度偏移,可得到比常规广义屏传播算子更好的效果.高阶广义屏传播算子的阶数越高,计算精度越高,但计算量也越多.以SEG EAGE二维盐丘模型数据的波动方程叠前深度偏移为例,二阶广义屏传播算子相对于常规（一阶）广义屏传播算子增加了30%的计算量.高阶广义屏传播算子是常规广义屏传播算子理论的发展和完善. 相似文献
18.
Jiangtao Hu Huazhong Wang Zhongyu Fang Tiancai Li Jiannan Zhang 《Geophysical Prospecting》2016,64(6):1483-1497
Least‐squares reverse time migration provides better imaging result than conventional reverse time migration by reducing the migration artefacts, improving the resolution of the image and balancing the amplitudes of the reflectors. However, it is computationally intensive. To reduce its computational cost, we propose an efficient amplitude encoding least‐squares reverse time migration scheme in the time domain. Although the encoding scheme is effective in increasing the computational efficiency, it also introduces the well‐known crosstalk noise in the gradient that degrades the quality of the imaging result. We analyse the cause of the crosstalk noise using an encoding correlation matrix and then develop two numerical schemes to suppress the crosstalk noise during the inversion process. We test the proposed method with synthetic and field data. Numerical examples show that the proposed scheme can provide better imaging result than reverse time migration, and it also generates images comparable with those from common shot least‐squares reverse time migration but with less computational cost. 相似文献
19.
各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献
20.
发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据（Marmousi模型）和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法. 相似文献