共查询到20条相似文献，搜索用时 718 毫秒
1.
地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声. 相似文献
2.
《中国科学:地球科学(英文版)》2015,(3)
The South China Sea where water depth is up to 5000 m is the most promising oil and gas exploration area in China in the future.The seismic data acquired in the South China Sea contain various types of multiples that need to be removed before imaging can be developed.However,compared with the conventional reflection migration,multiples carry more information of the underground structure that helps provide better subsurface imaging.This paper presents a method to modify the conventional reverse time migration so that multiple reflections can migrate to their correct locations in the subsurface.This approach replaces the numerical impulsive source with the recorded data including primaries and multiples on the surface,and replaces the recorded primary reflection data with multiples.In the reverse time migration process,multiples recorded on the surface are extrapolated backward in time to each depth level,while primaries and multiples recorded on the surface are extrapolated forward in time to the same depth levels.By matching the difference between the primary and multiple images using an objective function,this algorithm improves the primary resultant image.Synthetic tests on Sigsbee2 B show that the proposed method can obtain a greater range and better underground illumination.Images of deep water in the South China Sea are obtained using multiples and their matching with primaries.They demonstrate that multiples can make up for the reflection illumination and the migration of multiples is an important research direction in the future. 相似文献
3.
VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性. 相似文献
4.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
5.
Extending illumination using all multiples: application to 3D acquisition geometry analysis 
下载免费PDF全文
下载免费PDF全文 Recent advances in survey design have led to conventional common‐midpoint‐based analysis being replaced by subsurface‐based seismic acquisition analysis, with emphasis on advanced techniques of illumination analysis. Among them is the so‐called focal beam method, which is a wave‐equation‐based seismic illumination analysis method. The objective of the focal beam method is to provide a quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on the resolution and angle‐dependent amplitude fidelity of the image. The method distinguishes between illumination and sensing capability of a particular acquisition geometry by computing the focal source beam and the focal detector beam, respectively. Sensing is related to the detection properties of a detector configuration, whereas illumination is related to the emission properties of a source configuration. The focal source beam analyses the incident wavefield at a specific subsurface grid point from all available sources, whereas the focal detector beam analyses the sensing wavefield reaching at the detector locations from the same subsurface grid point. In the past, this method could only address illumination by primary reflections. In this paper, we will extend the concept of the focal beam method to incorporate the illumination due to the surface and internal multiples. This in fact complies with the trend of including multiples in the imaging process. Multiple reflections can illuminate a target location from other angles compared with primary reflections, resulting in a higher resolution and an improved illumination. We demonstrate how an acquisition‐related footprint can be corrected using both the surface and the internal multiples. 相似文献
6.
A three-dimensional (3-D) kinematic migration algorithm for media in which migration velocity varies linearly with depth is developed, implemented and tested. The algorithm is based on the concept that a single reflection or diffraction in a (zero- or finite-offset) trace may have originated at any point on a constant traveltime surface within the Earth defined by the observed two-way traveltime. The envelope of all such constant time surfaces, for all observed reflections and diffractions produced by one reflector, is the desired migrated 3-D image. The optimal envelope position in depth is determined, beneath each point on a regular grid, by a statistical imaging condition; an incremental function of depth containing the number of constant time surfaces passing through that depth increment is cross-correlated with a Gaussian function whose width is chosen to correspond to the vertical scale of the features of interest. The numerical procedures are based on the observation that, in a medium in which velocity varies linearly with depth, ray segments are circular so traveltimes can be computed analytically. Also, traveltimes are independent of azimuth so the 3-D problem can be collapsed into an equivalent 2-D problem. The algorithm is illustrated and tested by application to synthetic data and to scale-model data from the Seismic Acoustics Laboratory at the University of Houston. 相似文献
7.
使用相同的炮记录,多次波偏移能提供比反射波偏移更广的地下照明和更多的地下覆盖但是同时产生很多的串声噪声.相比传统逆时偏移,最小二乘逆时偏移反演的反射波成像结果具有更高的分辨率和更均衡的振幅.我们主要利用最小二乘逆时偏移压制多次波偏移产生的串声噪声.多次波最小二乘逆时偏移通常需要一定的迭代次数以较好地消除串声噪声.若提前将一阶多次波从所有阶数的多次波中过滤出来,使用相同的迭代次数,一阶多次波的最小二乘逆时偏移能够得到具有更高信噪比的成像剖面,而且能够提供与多次波最小二乘逆时偏移相似的有效地下结构成像. 相似文献
8.
在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义. 相似文献
9.
研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献
10.
在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善. 相似文献
11.
局域化相空间中的VSP偏移成像方法(英文) 总被引：1，自引：1，他引：0
对VSP资料进行偏移成像可提高井附近地下结构的成像分辨率。本文给出了一种基于局域化相空间波场分解的VSP偏移成像方法。此方法采用了基于Gabor-Daubechies紧标架的延拓算子(G-D延拓算子)及其高频渐近形式对相空间波场进行延拓;基于局部平面假设,提出了一种局部角度域相关成像条件。合成和实际VSP资料的偏移成像结果表明,在满足渐近展开的条件下利用G-D延拓算子的高频近似式能够有效的减少计算时间;采用局部角度域相关成像条件能够在不增加计算量的同时,有效减弱VSP成像剖面上的偏移假象。 相似文献
12.
The imaging of offset VSP data in local phase space can improve the image of the subsurface structure near the well. In this
paper, we present a migration scheme for imaging VSP data in a local phase space, which uses the Gabor-Daubechies tight framebased
extrapolator (G-D extrapolator) and its high-frequency asymptotic expansion to extrapolate wavefields and also delineates
an improved correlation imaging condition in the local angle domain. The results for migrating synthetic and real VSP data
demonstrate that the application of the high-frequency G-D extrapolator asymptotic expansion can effectively decrease computational
complexity. The local angle domain correlation imaging condition can be used to weaken migration artifacts without increasing
computation. 相似文献
13.
14.
Clement Fleury 《Geophysical Prospecting》2013,61(5):891-906
Reverse‐time migration is a two‐way time‐domain finite‐frequency technique that accurately handles the propagation of complex scattered waves and produces a band‐limited representation of the subsurface structure that is conventionally assumed to be linear in the contrasts in model parameters. Because of this underlying linear single‐scattering assumption, most implementations of this method do not satisfy the energy conservation principle and do not optimally use illumination and model sensitivity of multiply scattered waves. Migrating multiply scattered waves requires preserving the non‐linear relation between the image and perturbation of model parameters. I modify the extrapolation of source and receiver wavefields to more accurately handle multiply scattered waves. I extend the concept of the imaging condition in order to map into the subsurface structurally coherent seismic events that correspond to the interaction of both singly and multiply scattered waves. This results in an imaging process referred to here as non‐linear reverse‐time migration. It includes a strategy that analyses separated contributions of singly and multiply scattered waves to a final non‐linear image. The goal is to provide a tool suitable for seismic interpretation and potentially migration velocity analysis that benefits from increased illumination and sensitivity from multiply scattered seismic waves. It is noteworthy that this method can migrate internal multiples, a clear advantage for imaging challenging complex subsurface features, e.g., in salt and basalt environments. The results of synthetic seismic imaging experiments, including a subsalt imaging example, illustrate the technique. 相似文献
15.
对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善. 相似文献
16.
Two-dimensional VSP surveys are often conducted to provide structural illumination of the subsurface away from the borehole. The illumination is achieved through offsetting the source with respect to the downhole geophone. This inevitably gives rise to mode-conversions in both downgoing and upgoing wavefields. Migration of mixed-mode wavefields is complex because the velocity profile used for wavefield extrapolation is valid only for a particular propagation mode; the other mode always propagates at a different velocity. It is therefore advisable to separate the wave-types (P-wave and SV-wave) prior to migration. This may be achieved through wavemode filtering, a multichannel process which exploits the relation between propagation velocity, slowness of events at the recording array and particle motion. The necessary information about particle motion is available only if the VSP data are acquired with a three-component downhole geophone assembly. The wavemode filter partitions wave-types at the recording array; it provides no information about the various changes of propagation mode experienced by the energy as it travels from source to geophone. For the purpose of migration, the intermediate modes of propagation must be deduced. Much of the energy arriving at the receivers is P-wave which has followed the P-wave velocity profile from the source. It can therefore be imaged by conventional (Kirchhoff) migration. As an example of SV-wave imaging, a common mode-code is P-wave from source to reflector and SV-wave from reflector to geophone. Migration of such data calls for back-propagation of the geophone array wavefield, at SV-wave velocity, to the point in the subsurface where it is time-coincident with the forward propagated downwave, at P-wave velocity. 相似文献
17.
COMPUTATION OF ZERO-OFFSET VERTICAL SEISMIC PROFILES INCLUDING GEOMETRICAL SPREADING AND ABSORPTION*
Synthetic vertical seismic profiles (VSP) provide a useful tool in the interpretation of VSP data, allowing the interpreter to analyze the propagation of seismic waves in the different layers. A zero-offset VSP modeling program can also be used as part of an inversion program for estimating the parameters in a layered model of the subsurface. Proposed methods for computing synthetic VSP are mostly based on plane waves in a horizontally layered elastic or anelastic medium. In order to compare these synthetic VSP with real data a common method is to scale the data with the spherical spreading factor of the primary reflections. This will in most cases lead to artificial enhancement of multiple reflections. We apply the ray series method to the equations of motion for a linear viscoelastic medium after having done a Fourier transformation with respect to the time variable. This results in a complex eikonal equation which, in general, appears to be difficult to solve. For vertically traveling waves in a horizontally layered viscoelastic medium the solution is easily found to be the integral along the ray of the inverse of the complex propagation velocity. The spherical spreading due to a point source is also complex, and it is equal to the integral along the ray of the complex propagation velocity. Synthetic data examples illustrate the differences between spherical, cylindrical, and plane waves in elastic and viscoelastic layered media. 相似文献
18.
Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated. 相似文献
19.
In this paper, we image the subsurface reflectors by interferometric imaging using primary and downgoing first-order free-surface
related multiple reflections in walkaway VSP data. By analyzing the stack fold distribution, we find that primary-direct interferometric
imaging has a smaller image range, but its stack fold is higher near the well while ghost-direct interferometric imaging is
the opposite. We try to solve this problem by the joint interferometric imaging of walkaway VSP data, combining primary-direct
interferometric imaging with ghost-direct interferometric imaging. In this way, we can effectively widen the imaging range,
simultaneously increase the fold (especially near the well), suppress spurious interference, and improve the image SNR, so
that we can get a more credible image. Test results on synthetic walkaway VSP data and field data show that joint interferometric
imaging is very effective. 相似文献
20.
Walkaway VSP数据的联合相干成像(英文) 总被引：1，自引：1，他引：0
本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。 相似文献