退化的Fourier偏移算子及其在复杂断块成像中的应用

波动方程宽角抛物逼近得到的通常是非常系数的单程波传播算子,其系数是速度横向变化的函数,因此需要利用有限差分(FD)进行数值实施. 通过对Lippmann Schwinger单程波动积分方程的退化核逼近,本文研究了一类宽角退化算子的偏移成像. 这种退化偏移算子只用快速Fourier变换进行波场延拓,将常规的Fourier分裂步地震偏移方法(SSF)推广适应强速度横向变化介质和大角度传播波场. 退化的Fourier偏移算子通过在两个分裂步项之间作波数域线性插值来实现波场延拓,每延拓一层需要比常规的SSF地震偏移方法多一次快速Fourier变换(FFT). 通过SEG/EAGE盐丘模型和实际地震资料的应用表明,退化Fourier偏移算子能很好地对盐下的陡倾角断层和实际地震剖面上的复杂小断块和大断裂地质构造成像.     

High-order generalized screen propagator migration based on particle swarm optimization

Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unresolved. To describe the seismic wave propagation in media with lateral velocity variations and to image high-angle structures, we propose the generalized screen propagator based on particle swarm optimization (PSO-GSP), for the precise fitting of the single-square-root operator. We use the 2D SEG/EAGE salt model to test the proposed PSO-GSP migration method to image the faults beneath the salt dome and compare the results to those of the conventional high-order generalized screen propagator (GSP) migration and split-step Fourier (SSF) migration. Moreover, we use 2D marine data from the South China Sea to show that the PSO-GSP migration can better image strong reflectors than conventional imaging methods.     

基于多参量的模拟退火全局优化傅里叶有限差分算子

傅里叶有限差分法(FFD)能够处理复杂地质构造中的波传播问题,但对陡倾角成像仍有明显的误差.优化参数的方法能够在保持计算效率的前提下进一步提高陡倾角的成像精度.本文在有理近似的基础上,将FFD算子展开式中的常系数由两个拓展为四个,然后采用模拟退火算法对这四个参数进行全局优化.本方法除了考虑速度对比度以外,还考虑了频率和延拓步长等参量的影响.理论误差分析和脉冲响应测试均表明该方法能极大地提高FFD算子的精确传播角度.二维SEG/EAGE盐丘模型实验表明本文方法对陡倾角以及盐下构造的成像精度明显高于未优化的FFD法.将本文的方法与交替方向加插值的方法结合应用于三维脉冲响应测试更进一步证实了本文方法的有效性.     

RESIDUAL DEPTH-MIGRATION: ENHANCING DIPS1

Residual migration can correct for migration done with the wrong velocity. The initial migration can be done with a less accurate, fast time-migration algorithm followed by a residual depth-migration operator. This combines the speed and accuracy of different migration schemes. The residual depth-migration equation can be obtained by an appropriate transformation of a typical finite-difference, depth-migration equation. Since the cheaper, initial time-migration has already partially migrated the data, the residual depth-migration has less work to do than if the depth-migration were done in a single step. As a result, the overall process has increased the speed and reduced the dip limitations associated with typical depth-migration operators. Results on a steeply dipping model demonstrate the validity and usefulness of residual depth-migration.     

三维VTI介质中波动方程深度偏移的最优分裂Fourier方法

从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点.     

傅里叶有限差分法三维波动方程正演模拟

傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势,能够处理复杂地质构造中的波传播问题,但在三维情形下,算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项,消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构,明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题,并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比,结果证实了该方法对一次反射波具有较高的模拟精度,在内存需求和计算效率方面则具有更大的优势.     

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.     

NUMERICAL SEISMOGRAM COMPUTATIONS FOR INHOMOGENEOUS MEDIA USING A SHORT,VARIABLE LENGTH CONVOLUTIONAL DIFFERENTIATOR1

A short convolutional differentiator (CD) for computing second spatial derivatives in the acoustic wave equation is presented. This differentiator is obtained by tapering the inverse Fourier transform of the band-limited Fourier spectrum of the second-derivative operator. This new filter has been applied to seismogram computations for inhomogeneous media and results are compared with the conventional high-order finite-difference (FD) and Fourier schemes. The operator can be progressively shortened at the model edges to reduce boundary artefacts. The CD method is superior to the conventional FD operator and comparable with the Fourier method in accuracy but faster to run. A strategy to reduce computation time by 20%, which exploits the localized nature of the operator, is given. The method is illustrated using simple 2D models.     

混合域单程波传播算子及其在偏移成像中的应用

以地震波的单程波传播方程为基础,利用算子近似展开的方法推导出了当前波动方程叠前深度偏移方法研究中广泛使用的裂步Fourier、Fourier有限差分和广义屏传播算子的一般形式及其近似式．讨论了它们间差异、相互关系以及他们的特点,最后给出了基于裂步Fourier、Fourier有限差分和广义屏传播算子的偏移成像方法时Marmousi模型的偏移成像结果,以说明它们间的优劣与计算效率．     

自适应空间分区裂步傅立叶偏移方法

本文提出了一种新的偏移方法——自适应空间分区裂步傅立叶（ASDSSF）偏移方法。该方法将剥层相位移方法的思想推广到裂步傅立叶偏移方法，使之当速度场出现强间断时也能精确而有效地成像。原理上ASDSSF偏移属于多参考慢度（MRS）偏移方法，本文的重点是，在不损失精度的同时选取比同类MRS偏移方法更少的参考慢度。我们根据全局速度函数的变化和误差控制参数来选择参考慢度，同一个参考慢度所对应的速度构成一个分区，每一个空间分区可以由几个空间上不连续的子分区组成，从而有效地减小了参考慢度的个数。每一延拓步的参考慢度以及参考慢度的个数和如何构建分区都是根据速度函数自动生成，因此更为合理。为了消除速度场强间断产生的人为噪音，设计了简单有效的f-k域的光滑滤波。我们对一个生成的二维叠前模型和SEG／EAEG盐丘模型进行了试算。     

Broadband constant-coefficient propagators

The phase error between the real phase shift and the Gazdag background phase shift, due to lateral velocity variations about a reference velocity, can be decomposed into axial and paraxial phase errors. The axial phase error depends only on velocity perturbations and hence can be completely removed by the split‐step Fourier method. The paraxial phase error is a cross function of velocity perturbations and propagation angles. The cross function can be approximated with various differential operators by allowing the coefficients to vary with velocity perturbations and propagation angles. These variable‐coefficient operators require finite‐difference numerical implementation. Broadband constant‐coefficient operators may provide an efficient alternative that approximates the cross function within the split‐step framework and allows implementation using Fourier transforms alone. The resulting migration accuracy depends on the localization of the constant‐coefficient operators. A simple broadband constant‐coefficient operator has been designed and is tested with the SEG/EAEG salt model. Compared with the split‐step Fourier method that applies to either weak‐contrast media or at small propagation angles, this operator improves wavefield extrapolation for large to strong lateral heterogeneities, except within the weak‐contrast region. Incorporating the split‐step Fourier operator into a hybrid implementation can eliminate the poor performance of the broadband constant‐coefficient operator in the weak‐contrast region. This study may indicate a direction of improving the split‐step Fourier method, with little loss of efficiency, while allowing it to remain faster than more precise methods such as the Fourier finite‐difference method.     

串联F-K加相移插值法偏移

本文介绍的偏移成像方法为递归F-K法、相移插值法及串联偏移技术的有效结合,简称CFPI法.该方法中,波场外推采用相移插值法中递归向下延拓方式,而每一延拓层间的成像则用文中给出的常速F-K加剩余速度F-K偏移,通过与横向速度变化有关的插值来实现,该过程用FFT完成.该方法可采用较大的延拓步长速归向下外推与成像.理论分析及实际计算结果比较表明,该算法效率远高于相移插值法偏移.     

A self-adaptive algorithm for choosing reference velocities in the presence of lateral velocity variations

Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective.     

自适应横向变速参考速度选择方法

基于扰动理论的混合域波动方程波场外推算子,具有一定介质横向速度变化适应能力,是反射地震学中研究较多的成像方法。此类波场外推算子沿深度层进行波场外推,都存在参考速度选择问题。单参考速度波场外推算子,适应地下介质弱横向变速,而多参考速度波场外推算子可以提高横向变速的适应能力和成像精度,但要以大量计算为代价。本文提出的自适应多参考速度选择策略,根据外推层地质构造的复杂度和给定的速度门槛值自动选择参考速度个数,利用该策略构造混合域SSF、FFD、WXFD和GSP等多参考速度波场外推成像算法。盐丘模型理论数据测试结果表明,自适应多参考速度波场外推算法具有强横向变速适应能力和较高成像精度。     

单向传播子的局部平面波分解以及在定向照明分析和目标特征成像中的应用

本文基于波场的局部平面波分解,得到空间和方向都具有局部性的局部平面波小波束,并由此对Fourier有限差分传播算子进行局部平面波分解.这种用小波束分解的单向传播算子,可被用来对采集系统进行定向照明分析和采集倾角响应估算,以分析采集系统布局和上覆结构对面向目标的叠前偏移成像质量的影响.进而,利用局部平面波成像矩阵及其成像相册,根据目标结构的主倾角方向,抽取并叠加部分局部平面波像,可以进行目标结构定向化特征成像.最后,对具有速度强横向变化上覆结构的陡倾角目标断层的复杂二维SEG-EAGE盐丘模型数据,进行采集系统目标照明分析和目标结构定向化特征成像.     

FINITE DIFFERENCE AND WAVE NUMBER MIGRATION *

Finite difference migration has been developed and popularized by J. F. Claerbout of Stanford University and is now widely used in seismic processing. For most sections finite difference migration gives results comparable to those obtained by conventional Kirchhoff migration and, where events are not dipping too much, a cleaner appearance is often apparent. However, there are two practical limitations to the method, and these occur in regions of very steep dip and where there is a large variation of the velocity in the lateral direction. It is possible to develop successively more accurate equations to deal with the steep dip problem, but above third order these schemes become prohibitively expensive to implement. The finite difference method itself introduces errors and so imposes further limitations on the angle of dip. For the effective treatment of steeply dipping beds there appears to be no method available in the time domain which does not suffer from dispersion inaccuracies. However, by developing wavenumber migration, an exact one-way wave equation can be used, and this eliminates any error except that caused by finite sampling. The other difficulty with wave migration is the correct migration in regions with lateral velocity variation. A number of approaches are possible of which three are discussed here. The first uses an exact theory, the second is based on the deviation from a depth stratified model, and the third uses a transformation to a depth co-ordinate system. All methods are discussed with their advantages and limitations. Finally, some examples are shown of wave migration applied to synthetic and real data.     

基于全波形反演的探地雷达数据逆时偏移成像

逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果.     

Wavefield interpolation in 3D large‐step Fourier wavefield extrapolation

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.     

基于波场垂向和水平方向外推的高陡构造偏移

常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势.     

数据自相关多次波偏移成像

在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义.