共查询到20条相似文献，搜索用时 31 毫秒
1.
当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当. 相似文献
2.
Xiaoxue Jiang Fan Zheng Haiqing Jia Jun Lin Hongyuan Yang 《Studia Geophysica et Geodaetica》2016,60(1):91-111
A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Radon panel due to their differences in slowness, and an inverse transform of only part of the data leads to separated wavefields. A problem of hyperbolic Radon transform is that it works in the time domain entailing the inversion of large operators which is prohibitively time-consuming. By applying the conjugate gradient algorithm during the inversion of hyperbolic Radon transform, the computational cost can be kept reasonably low for practical application. Synthetic data examples prove that P-P and P-SV wavefield separation by hyperbolic Radon transform produces more accurate separated wavefields compared with separation by high-resolution parabolic Radon transform, and the feasibility of the proposed separation scheme is also verified by a real field data example. 相似文献
3.
Paul Sava 《Geophysical Prospecting》2012,60(1):29-48
Seismic waves propagate through the earth as a superposition of different wave modes. Seismic imaging in areas characterized by complex geology requires techniques based on accurate reconstruction of the seismic wavefields. A crucial component of the methods in this category, collectively known as wave‐equation migration, is the imaging condition that extracts information about the discontinuities of physical properties from the reconstructed wavefields at every location in space. Conventional acoustic migration techniques image a scalar wavefield representing the P‐wave mode, in contrast to elastic migration techniques, which image a vector wavefield representing both the P‐ and S‐waves. For elastic imaging, it is desirable that the reconstructed vector fields are decomposed into pure wave modes, such that the imaging condition produces interpretable images, characterizing, for example, PP or PS reflectivity. In anisotropic media, wave mode separation can be achieved by projection of the reconstructed vector fields on the polarization vectors characterizing various wave modes. For heterogeneous media, because polarization directions change with position, wave mode separation needs to be implemented using space‐domain filters. For transversely isotropic media with a tilted symmetry axis, the polarization vectors depend on the elastic material parameters, including the tilt angles. Using these parameters, we separate the wave modes by constructing nine filters corresponding to the nine Cartesian components of the three polarization directions at every grid point. Since the S polarization vectors in transverse isotropic media are not defined in the singular directions, e.g., along the symmetry axes, we construct these vectors by exploiting the orthogonality between the SV and SH polarization vectors, as well as their orthogonality with the P polarization vector. This procedure allows one to separate all three modes, with better preserved P‐wave amplitudes than S‐wave amplitudes. Realistic synthetic examples show that this wave mode separation is effective for both 2D and 3D models with strong heterogeneity and anisotropy. 相似文献
4.
C. P. A. WAPENAAR P. HERRMANN D. J. VERSCHUUR A. J. BERKHOUT 《Geophysical Prospecting》1990,38(6):633-661
Inversion of multicomponent seismic data can be subdivided in three main processes: (1) Surface-related preprocessing (decomposition of the multicomponent data into ‘primary’ P-and S-wave responses). (2) Prestack migration of the primary P- and S-wave responses, yielding the (angle-dependent) P-P, P-S, S-P and S-S reflectivity of the subsurface. (3) Target-related post-processing (transformation of the reflectivity into the rock and pore parameters in the target). This paper deals with the theoretical aspects of surface-related preprocessing. In a multicomponent seismic data set the P- and S-wave responses of the subsurface are distorted by two main causes: (1) The seismic vibrators always radiate a mixture of P- and S-waves into the subsurface. Similarly, the geophones always measure a mixture of P- and S-waves. (2) The free surface reflects any upgoing wave fully back into the subsurface. This gives rise to strong multiple reflections, including conversions. Therefore, surface-related preprocessing consists of two steps: (1)Decomposition of the multicomponent data (pseudo P- and S-wave responses) into true P- and S-wave responses. In practice this procedure involves (a) decomposition per common shot record of the particle velocity vector into scalar upgoing P- and S-waves, followed by (b) decomposition per common receiver record of the traction vector into scalar downgoing P- and S-waves. (2) Elimination of the surface-related multiple reflections and conversions. In this procedure the free surface is replaced by a reflection-free surface. The effect is that we obtain ‘primary’ P-and S-wave responses, that contain internal multiples only. An interesting aspect of the procedure is that no knowledge of the subsurface is required. In fact, the subsurface may have any degree of complexity. Both the decomposition step and the multiple elimination step are fully determined by the medium parameters at the free surface only. After surface-related preprocessing, the scalar P- and S-wave responses can be further processed independently by existing scalar algorithms. 相似文献
5.
将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率. 相似文献
6.
7.
Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration 
下载免费PDF全文
下载免费PDF全文 The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections. 相似文献
8.
We introduce a new method for prestack depth migration of seismic common-shot gathers. The computational procedure follows
standard steps of the reverse-time migration, i.e., downward continuation of the source and the receiver wavefields, followed
by application of an imaging condition (e.g. zero-lag cross-correlation of these fields). In our method we first find a sparse
data representation with a small number of Gaussian wave packets. We then approximate the downward wavefield propagation (for
the source and the receiver fields) by a rigid flow of these wave packets along seismic rays. In this case, the wave packets
are simply translated and rotated according to the ray geometry. One advantage of using Gaussian wave packets is that analytic
formulas can be used for translation, rotation, and the application of the cross-correlation imaging condition. Moreover,
they allow more sparse representations than competing methods. Finally we formulate a computationally and memory efficient
migration procedure, as only few rays have to be traced, and since it is cheap to compute the cross-correlation for the intersecting
rays. 相似文献
9.
The conventional reverse time migration of ground-penetrating radar data is implemented with the two-way wave equation. The cross-correlation result contains low-frequency noise and false images caused by improper wave paths. To eliminate low-frequency noise and improve the quality of the migration image, we propose to separate the left-up-going, left-down-going, right-up-going and right-down-going wavefield components in the forward- and backward-propagated wavefields based on the Hilbert transform. By applying the reverse time migration of ground-penetrating radar data with full wavefield decomposition based on the Hilbert transform, we obtain the reverse time migration images of different wavefield components and combine correct imaging conditions to generate complete migration images. The proposed method is tested on the synthetic ground-penetrating radar data of a tilt-interface model and a complex model. The migration results show that the imaging condition of different wavefield components can highlight the desired structures. We further discuss the reasons for incomplete images by reverse time migration with partial wavefields. Compared with the conventional reverse time migration methods for ground-penetrating radar data, low-frequency noise can be eliminated in images generated by the reverse time migration method with full wavefield decomposition based on the Hilbert transform. 相似文献
10.
Bergler Steffen Duveneck Eric Höcht German Zhang Yonghai Hubral Peter 《Studia Geophysica et Geodaetica》2002,46(2):165-175
The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections. 相似文献
11.
The polarization direction or 'sign’ of reflected converted P–S waves depends upon the angle of incidence of the incident P-wave. Sign reversal due to reversal of the angle of incidence is often encountered and is an impediment to P–S wave processing and imaging, because when P–S events or P-S migrated images with mixed signs are stacked, destructive interference occurs. We have created and demonstrated a means of correcting for this reversal. To do this, a P-wave angle of incidence is calculated for every point in the image space. This is done by calculating a P–S reflected waveform for every point, by extrapolating the reflected S-wavefield backwards from the receiver line, and then cross-correlating this waveform with the S-wave reflections observed at the receiver line. A multiplier, (sgn α) is assigned to each point in the image space, where α is the angle of incidence of the P-wave. The multiplier was applied to a set of prestack reverse time migration images derived from a cross-borehole physical elastic model data set. The improvement in the stacked image when the sign correction is applied is spectacular. The P-S image quality is comparable to, or better than, stacked migrated P-P images. The method appears to be applicable to all reflection modes and to all recording geometries. 相似文献
12.
Reverse-time migration has attracted more and more attention owing to the advantages of high imaging accuracy, no dip restriction, and adaptation to complex velocity models. Cross-correlation imaging method is typically used in conventional reverse-time migration that produces images with strong low-frequency noise. Wavefield decomposition imaging can suppress such noise; however, some residual noise persists in the imaging results. We propose a 2D multidirectional wavefield decomposition method based on the traditional wavefield decomposition method. First, source wavefields and receiver wavefields are separated into eight subwavefields, respectively. Second, cross-correlation imaging is applied to selected subwavefields to produce subimages. Finally, the subimages are stacked to generate the final image. Numerical examples suggest that the proposed method can eliminate the low-frequency noise effectively and produce high-quality imaging profiles. 相似文献
13.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
14.
M. NOVOTNÝ 《Geophysical Prospecting》1990,38(8):833-851
The slant-stack migration formula based on the Radon transform is studied with respect to the depth step Δz of wavefield extrapolation. It can be viewed as a generalized trace-interpolation procedure including wave extrapolation with an arbitrary step Δz. For Δz= 0 the formula yields the familiar plane-wave decomposition, while for Δz > 0 it provides a robust tool for migration transformation of spatially undersampled wavefields. Using the stationary phase method, it is shown that the slant-stack migration formula degenerates into the Rayleigh-Sommerfeld integral in the far-field approximation. Consequently, even a narrow slant-stack gather applied before the diffraction stack can significantly improve the representation of noisy data in the wavefield extrapolation process. The theory is applied to synthetic and field data to perform trace interpolation and dip reject filtration. The data examples presented prove that the Radon interpolator works well in the dip range, including waves with mutual stepouts smaller than half the dominant period. 相似文献
15.
在各向异性地震波场中,qP波与qS波常常是耦合在一起的.多分量地震数据处理中一个关键环节就是波型分离(即模式解耦),以纵波成分为主的常规单分量地震数据的成像则需要合理描述标量qP波的传播算子.本文作者曾构建了在运动学上同弹性波动方程等价,动力学上突出标量qP波的伪纯模式波动方程.为了彻底消除qS波残余,本文根据波矢量与qP波偏振矢量之间的偏差,提出从伪纯模式波场提取纯模式标量qP波的方法.数值分析展示了投影偏差算子在波数域和空间域的特征.基于不同复杂程度理论模型的试验结果表明,联合"伪纯模式传播算子"与"投影偏差校正"可为各向异性介质分离模式波场传播过程提供一种简便的描述工具. 相似文献
16.
VSP 转换波跟VSP 纵波或常规地面转换波相比,具有较高的分辨率和信噪比,但传统的VSP成像方法只利用了反射P波信息,而把转换波（反射S波、透射S波）以及透射P波当作影响成像质量的噪音.本文给出了一种VSP共炮点道集多分量地震资料克希霍夫法偏移成像的方法.本方法充分利用了多波（反射P波、反射S波、透射P波、透射S波）信息,根据转换点处四种波同时起跳,能量叠加最大的原理,从接收点分别用向绕射点延拓它们的能量,并将其叠加起来,求得的和最大的一点即反射点.通过模型试算和实际资料处理表明,此法成像精度高,信噪比高,且有利于改善剖面的频率特性. 相似文献
17.
Extrapolating wavefields and imaging at each depth during three‐dimensional recursive wave‐equation migration is a time‐consuming endeavor. For efficiency, most commercial techniques extrapolate wavefields through thick slabs followed by wavefield interpolation within each thick slab. In this article, we develop this strategy by associating more efficient interpolators with a Fourier‐transform‐related wavefield extrapolation method. First, we formulate a three‐dimensional first‐order separation‐of‐variables screen propagator for large‐step wavefield extrapolation, which allows for wide‐angle propagations in highly contrasting media. This propagator significantly improves the performance of the split‐step Fourier method in dealing with significant lateral heterogeneities at the cost of only one more fast Fourier transform in each thick slab. We then extend the two‐dimensional Kirchhoff and Born–Kirchhoff local wavefield interpolators to three‐dimensional cases for each slab. The three‐dimensional Kirchhoff interpolator is based on the traditional Kirchhoff formula and applies to moderate lateral velocity variations, whereas the three‐dimensional Born–Kirchhoff interpolator is derived from the Lippmann–Schwinger integral equation under the Born approximation and is adapted to highly laterally varying media. Numerical examples on the three‐dimensional salt model of the Society of Exploration Geophysicists/European Association of Geoscientists demonstrate that three‐dimensional first‐order separation‐of‐variables screen propagator Born–Kirchhoff depth migration using thick‐slab wavefield extrapolation plus thin‐slab interpolation tolerates a considerable depth‐step size of up to 72 ms, eventually resulting in an efficiency improvement of nearly 80% without obvious loss of imaging accuracy. Although the proposed three‐dimensional interpolators are presented with one‐way Fourier extrapolation methods, they can be extended for applications to general migration methods. 相似文献
18.
双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象. 相似文献
19.
Xiang-Yang Li 《Geophysical Prospecting》1997,45(1):39-64
The characteristic seismic response to an aligned-fracture system is shear-wave splitting, where the polarizations, time-delays and amplitudes of the split shear waves are related to the orientation and intensity of the fracture system. This offers the possibility of delineating fractured reservoirs and optimizing the development of the reservoirs using shear-wave data. However, such applications require carefully controlled amplitude processing to recover properly and preserve the reflections from the target zone. Here, an approach to this problem is suggested and is illustrated with field data. The proposed amplitude processing sequence contains a combination of conventional and specific shear-wave processing procedures. Assuming a four-component recording (two orthogonal horizontal sources recorded by two orthogonal horizontal receivers), the split shear waves can be simulated by an effective eigensystem, and a linear-transform technique (LTT) can be used to separate the recorded vector wavefield into two principal scalar wavefields representing the fast and slow split shear waves. Conventional scalar processing methods, designed for processing P-waves, including noise reduction and stacking procedures may be adapted to process the separated scalar wavefields. An overburden operator is then derived from and applied to the post-stacked scalar wavefields. A four-component seismic survey with three horizontal wells drilled nearby was selected to illustrate the processing sequence. The field data show that vector wavefield decomposition and overburden correction are essential for recovering the reflection amplitude information in the target zone. The variations in oil production in the three horizontal wells can be correlated with the variations in shear-wave time-delays and amplitudes, and with the variations in the azimuth angle between the horizontal well and the shear-wave polarization. Dim spots in amplitude variations can be correlated with local fracture swarms encountered by the horizontal wells. This reveals the potential of shear waves for fractured reservoir delineation. 相似文献
20.
The state-of-the-art joint migration inversion faces the so-called amplitude-versus-offset challenge, due to adopting over-simplified one-way propagation, reflection and transmission operators to avoid over-parameterization in the inversion process. To overcome this challenge, we apply joint migration inversion to horizontally layered media (or 1.5-dimensional media) and parameterize the solution space via density and velocity models. In this scenario, one-way propagation, reflection and transmission operators required by the joint migration inversion process can be analytically and correctly derived from the subsurface models, so the amplitude-versus-offset challenge is successfully overcome. We introduce a new concept, which is named ‘inverse propagation’, into our 1.5-dimensional amplitude-versus-offset joint migration inversion. It can correctly reconstruct subsurface wavefields by using a surface-recorded receiver wavefield with all the influence of transmission, reflection and multiples accounted for. A synthetic example is used to demonstrate the correctness of the inverse propagation. This work is the foundation to further develop the 1.5-dimensional amplitude-versus-offset joint migration inversion technology. 相似文献