基于波射线路径偏移压制多次波

波射线路径压制多次波的反射波成像是在偏移过程去除多次波同时仅对反射波成像.通过在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角计算射线的路径.从炮点入射的射线与从检波点出射的射线的交点形成的走时,若等于观测走时,可以判断此条射线是反射波;反之,若不相等,则是多次波.数值实验表明此方法可以有效地去掉由于多次波能量产生的假成像点和压制多次波,因此界面可以正确归位,同时去掉由于多次波引起的假成像位置.     

波路径偏移压制层间多次波的理论与应用

消除层间多次波是地震勘探资料处理研究领域的难题,尤其对于实际资料的处理,到目前为止还很难找到一种完全有效的方法. 本文给出了仅对一次波成像既波路径偏移方法压制层间多次波方法,在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角,由此计算的角度作为射线追踪的初始角度,计算地震波射线的传播路径. 结合由程函方程计算的走时表,判断偏移范围是反射波还是多次波. 在前期偏移过程压制多次波的理论研究基础上,本文主要研究波路径偏移消除多次波的应用部分. 为了进一步说明效果的有效性,计算了在单炮和共成像点道集压制层间多次波,给出了实际资料的压制多次波的偏移结果.     

Cost-effective 3D one-pass depth migration1

A crucial point in the processing of 3D seismic data is the migration step, both because of its 3D nature and the computational cost involved. The efficiency and accuracy of 3D migration are determined by the wavefield extrapolation technique employed. Wavefield extrapolation based on second-order differential operators of variable-length is very efficient and accurate at the same time. Compared to migration based on the McClellan transform and operator splitting, the use of variable-length second-order differential operators offers significant advantages. The 3D migration operator has an almost perfect circular symmetry. No positioning errors in the 45° azimuth between the in-line and cross-line directions are evident. The method is, in practice, only limited by spatial aliasing and does not require expensive interpolation of data to reduce numerical artifacts. This reduces the computational cost of 3D one-pass depth migration by a large factor.     

Parsimonious 3D post-stack Kirchhoff depth migration

Parsimonious post‐stack migration is extended to three dimensions. By tracing single rays back along each incident wave direction (as determined by a local slant stack at the receivers), the ray tracing can be embedded in the migration. This approach significantly reduces the computer time and disk space needed because it is not necessary to build and save image time maps; 3D migration can be performed on a workstation or personal computer rather than using a supercomputer or cluster. The location of a reflector in the output image is defined by tracing a zero‐offset ray to the one‐way traveltime (the image condition); the orientation of the reflector is defined as a surface perpendicular to the raypath. The migration impulse response operator is confined to the first Fresnel zone around the estimated reflection point, which is much smaller than the large isochronic surface in traditional Kirchhoff depth migration. Additional efficiency is obtained by applying an amplitude threshold to reduce the amount of data to be migrated. Tests on synthetic data show that the proposed implementation of parsimonious 3D post‐stack Kirchhoff depth migration is at least two orders of magnitude faster than traditional Kirchhoff migration, at the expense of slightly degraded migration image coherence. The proposed migration is expected to be a useful complement to conventional time migrations for fast initial imaging of subsurface structures and for real‐time imaging of near‐offset sections during data acquisition for quality control.     

三维偏移距平面波有限差分叠前时间偏移

本文提出了中点-半偏移距域内的三维偏移距平面波(offset plane-wave)方程,并给出了其有限差分解法.偏移距平面波可通过对CMP道集进行平面波分解(倾斜叠加或线性Radon变换)生成,然而这样做会产生严重的噪音干扰.本文提出了局部倾斜叠加方法(local slant-stacking)来消除离散线性Radon变换引入的噪音.针对实际三维数据的不规则性(中点-偏移距域内方位角展布不均匀及偏移距采样不规则),本文还提出了与方位角无关的三维倾斜叠加方法(azimuth-independent 3D slant-stacking),解决了三维平面波分解中存在的问题.使用文中提出的平面波分解方法,可以得到高信噪比的偏移距平面波数据体.同时,三维偏移距平面波偏移可以输出偏移距射线参数域共成像点道集,基于此道集的剩余速度分析方法可以用来更新偏移速度场.偏移距平面波偏移具有很高的计算效率,相较Kirchhoff积分叠前时间偏移有较好的保幅特性,可作为水平地表三维叠前时间偏移的一个很好的解决方案.     

起伏地表采集数据的三维直接叠前时间偏移方法

提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性.     

3D depth migration by rotated McClellan filters1

The application of McClellan transformations considerably reduces the computational cost of 3D wavefield depth extrapolation by explicit convolutional methods. The accuracy of migration methods based on McClellan transformation depends on how well the transformation filter (cos !;κ!;) is approximated; errors in this approximation cause anisotropy in the extrapolation operator and frequency dispersion in the migrated results. The anisotropy can be greatly reduced by rotating the approximate filter by 45° and averaging the rotated filter with the original filter. The application of the rotated filter yields a migration method that correctly images very steep dips, with little or no additional computational cost. McClellan migration with the improved circular response enhances the imaging of synthetic and real data.     

Efficient traveltime compression for 3D prestack Kirchhoff migration

Kirchhoff 3D prestack migration, as part of its execution, usually requires repeated access to a large traveltime table data base. Access to this data base implies either a memory intensive or I/O bounded solution to the storage problem. Proper compression of the traveltime table allows efficient 3D prestack migration without relying on the usually slow access to the computer hard drive. Such compression also allows for faster access to desirable parts of the traveltime table. Compression is applied to the traveltime field for each source location on the surface on a regular grid using 3D Chebyshev polynomial or cosine transforms of the traveltime field represented in the spherical coordinates or the Celerity domain. We obtain practical compression levels up to and exceeding 20 to 1. In fact, because of the smaller size traveltime table, we obtain exceptional traveltime extraction speed during migration that exceeds conventional methods. Additional features of the compression include better interpolation of traveltime tables and more stable estimates of amplitudes from traveltime curvatures. Further compression is achieved using bit encoding, by representing compression parameters values with fewer bits.     

Numerical tests of 3D true-amplitude zero-offset migration1

An amplitude-preserving migration aims at imaging compressional primary (zero-or) non-zero-offset reflections into 3D time or depth-migrated reflections so that the migrated wavefield amplitudes are a measure of angle-dependent reflection coeffcients. The principal objective is the removal of the geometrical-spreading factor of the primary reflections. Various migration/inversion algorithms involving weighted diffraction stacks proposed recently are based on Born or Kirchhoff approximations. Here, a 3D Kirchhoff-type zero-offset migration approach, also known as a diffraction-stack migration, is implemented in the form of a time migration. The primary reflections of the wavefield to be imaged are described a priori by the zero-order ray approximation. The aim of removing the geometrical- spreading loss can, in the zero-offset case, be achieved by not applying weights to the data before stacking them. This case alone has been implemented in this work. Application of the method to 3D synthetic zero-offset data proves that an amplitude-preserving migration can be performed in this way. Various numerical aspects of the true-amplitude zero-offset migration are discussed.     

3D Fourier finite-difference migration by alternating-direction-implicit plus interpolation

Conventional two‐way splitting Fourier finite‐difference migration for 3D complex media yields azimuthal anisotropy where an additional phase correction is needed with much increase of computational cost. We incorporate the alternating‐direction‐implicit plus interpolation scheme into the conventional Fourier finite‐difference method to reduce azimuthal anisotropy. This scheme retains the high‐order remnants ignored by the two‐way splitting in the form of a wavefield interpolation in the wavenumber domain. The wavefield interpolation for each step of downward extrapolation is implemented between the wavefields before and after the conventional Fourier finite‐difference extrapolation. As the Fourier finite‐difference migration is implemented in the space and wavenumber dual space, the Fourier transforms between space and wavenumber domain that were needed for the alternating‐direction‐implicit plus interpolation in frequency domain (FD) migration are saved in Fourier finite‐difference migration. Since the azimuth anisotropy in Fourier finite‐difference is much less than that in FD, the application of the alternating‐direction‐implicit plus interpolation scheme in Fourier finite‐difference migration is superior to that in FD migration in handling complex media with large velocity contrasts and steep dips. Impulse responses show that the presented method reduces the azimuthal anisotropy at almost no extra cost.     

3D angle-domain common-image gathers for migration velocity analysis

Angle‐domain common‐image gathers (ADCIGs) are an essential tool for migration velocity analysis (MVA). We present a method for computing ADCIGs in 3D from the results of wavefield‐continuation migration. The proposed methodology can be applied before or after the imaging step in a migration procedure. When computed before imaging, 3D ADCIGs are functions of the offset ray parameters (p, p) ; we derive the geometric relationship that links the offset ray parameters to the aperture angle γ and the reflection azimuth φ. When computed after imaging, 3D ADCIGs are directly produced as functions of γ and φ. The mapping of the offset ray parameters (p, p) into the angles (γ, φ) depends on both the local dips and the local interval velocity; therefore, the transformation of ADCIGs computed before imaging into ADCIGs that are functions of the actual angles is difficult in complex structure. By contrast, the computation of ADCIGs after imaging is efficient and accurate even in the presence of complex structure and a heterogeneous velocity function. On the other hand, the estimation of the offset ray parameters (p, p) is less sensitive to velocity errors than the estimation of the angles (γ, φ). When ADCIGs that are functions of the offset ray parameters (p, p) are adequate for the application of interest (e.g. ray‐based tomography), the computation of ADCIGs before imaging might be preferable. Errors in the migration velocity cause the image point in the angle domain to shift along the normal to the apparent geological dip. By assuming stationary rays (i.e. small velocity errors), we derive a quantitative relationship between this normal shift and the traveltime perturbation caused by velocity errors. This relationship can be directly used in an MVA procedure to invert depth errors measured from ADCIGs into migration velocity updates. In this paper, we use it to derive an approximate 3D residual moveout (RMO) function for measuring inconsistencies between the migrated images at different γ and φ. We tested the accuracy of our kinematic analysis on a 3D synthetic data set with steeply dipping reflectors and a vertically varying propagation velocity. The tests confirm the accuracy of our analysis and illustrate the limitations of the straight‐ray approximation underlying our derivation of the 3D RMO function.     

On the aperture effect in 3D Kirchhoff-type migration

It is well known that the migrated image given by a Kirchhoff-type (diffraction-stack) migration with limited aperture is always accompanied by some events which depend on the migration aperture. Although these events may severely affect the quality of migration, they have been studied only in 2D cases. Here, the events due to the migration aperture in 3D situations are investigated using a new method of analysing the reconstructed wavefield. It is found that a finite migration aperture results in a reconstructed wavefield with two components. One comes from the tangent points and curves between the traveltime surfaces of reflected and point-diffracted rays and is independent of the migration aperture, and the other is from the boundary of the migration aperture and depends strongly on the location and size as well as on the shape of the migration aperture. It is this last component that describes the aperture effect in migration. If the migration aperture is not sufficiently large, and if the input for migration is not zero on the boundary of the migration aperture, the boundary component may partially or totally cancel the migration signal. Furthermore, for synthetic data, the aperture effect cannot be eliminated by enlarging the migration aperture because, except for the common-shotpoint data, the aperture effect always exists however large the migration aperture becomes. This leads to the conclusion that the published Kirchhoff-type operators are not the exact inverse operators of the Fresnel–Kirchhoff integral if the input data are synthetic.     

三维最小二乘弹性高斯束叠前深度偏移

与声波高斯束成像相比,弹性波高斯束偏移更适用于复杂油气藏多波多分量地震勘探.但是由于观测系统的局限性和深部构造的复杂性,该方法同样存在成像分辨率低、照明不均衡等问题.本文结合最小二乘偏移和弹性高斯束偏移的优势,提出了一种通过弹性高斯束叠加构建Born正演（反偏移）算子和偏移算子的三维最小二乘叠前深度偏移方法.依据最小二乘反演理论,建立基于反偏移数据与实际观测数据残差的目标函数,采用共轭梯度算法迭代更新来建立地下真实的反射率.与传统弹性高斯束偏移方法相比,该方法不仅提高了成像分辨率,而且使复杂构造特别是陡倾角地层的成像照明也得到了补偿.理论模型测试结果证明了本文方法的可行性和有效性.

     

Eliminating 3D pre‐stack migration artefacts by 6D filtering

State‐of‐the‐art 3D seismic acquisition geometries have poor sampling along at least one dimension. This results in coherent migration noise that always contaminates pre‐stack migrated data, including high‐fold surveys, if prior‐to‐migration interpolation was not applied. We present a method for effective noise suppression in migrated gathers, competing with data interpolation before pre‐stack migration. The proposed technique is based on a dip decomposition of common‐offset volumes and a semblance‐type measure computation via offset for all constant‐dip gathers. Thus the processing engages six dimensions: offset, inline, crossline, depth, inline dip, and crossline dip. To reduce computational costs, we apply a two‐pass (4D in each pass) noise suppression: inline processing and then crossline processing (or vice versa). Synthetic and real‐data examples verify that the technique preserves signal amplitudes, including amplitude‐versus‐offset dependence, and that faults are not smeared.     

基于Born波路径的高斯束初至波波形反演

为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低.     

基于纵横波解耦的三维弹性波逆时偏移

纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性.

     

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀疏采样情况下可实现陡倾角构造正确成像的反假频波场延拓算法.SEG二维Hess VTI各向异性介质理论模型数据及野外数据集观测系统下的三维脉冲响应计算表明,本文提出了一种具有工业应用潜力的三维各向异性波动方程叠前深度方法.     

三维VSP数据高效偏移成像的超道集方法

当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当.     

三维电磁偏移数值滤波器实现及参数分析

基于满足扩散方程的电磁偏移场解析式,推导了能够应用于实际数据处理的电磁偏移数值滤波器.为了合理选取滤波器参数,本文基于偏移微分方程,构建了满足偏移微分方程的解析测试场,通过测试场与偏移数值滤波器处理结果比较,重点研究了滤波器窗口大小、采样距离、偏移深度等滤波器参数对数值偏移滤波效果的影响,确定了数值滤波器参数选择原则.考虑到实际应用,在计算偏移电导率时,本文利用背景场定义了反射函数.建立了三维典型地电模型做数值计算,利用电磁偏移数值滤波器处理模拟数据,得到了合理的目标体位置和形态.测试场比较和数值试验证明,电磁偏移成像是一种可靠稳定的电磁资料解释技术.     

三维共偏移距叠前地震数据反假频射线束形成偏移成像

对稀疏/非规则采样或者低信噪比数据,射线束提取困难并伴随有假频产生,对叠加剖面和道集造成严重干扰.为了提升射线束偏移在稀疏和低信噪比地震数据采集中的成像效果,本文提出基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法.射线束偏移首先将地震数据划分为超道集,经过部分NMO后转化为以射线束中心定义的共偏移距数据,倾斜叠加和反假频操作均在局部共中心点坐标上实现.时间域倾斜叠加是对地震数据的时移累加操作,三角低通滤波同样可以在时间域完成,在对地震数据进行因果和反因果积分后,亦为地震数据的时移累加.因此,三角低通滤波与倾斜叠加可在时间域结合同时完成,避免了频域滤波的正反傅里叶变换.本文在反假频公式中加入权重系数,用以对反假频的程度进行控制,达到分辨率和噪声压制的最佳折衷.以某海上三维实际数据为例,文中展示了反假频射线束形成对偏移叠加剖面和共成像点偏移距道集中的噪声进行了有效压制.