三维french模型地震正演问题的全波有限差分方法.docxVIP

三维french模型地震正演问题的全波有限差分方法

1全波正演方法

波动方程的正转模拟在反演、解释和观测系统设计中起着重要作用。采用最常用的波动方程的方法是基于有限差、有限差和伪谱法等全波方程，从而模拟各种波动现象，并可应用于复杂介质模型。然而，计算和存储成本相对较高。在3d问题上，当3d波方程模拟时，当前波场的波场从整个三维模型的第一个（或更多）分钟的波场值中得到。相比之下，波动方程将3d模型划分为一系列二维平面成像单元。当每个薄膜层的波场计算时，只知道前层的薄面，而不是整个模型的波场值。因此，波动波模拟的效率高，对存储的要求也显著下降。它比全波方程的模拟方法更适合于当前的三维地震勘探模型计算。此外，在当前的地震勘探中，大多数地震勘探是基于反射波的有效信息，而直接波和多个波被认为是噪声信息。全波方程不仅可以保证反射波的模拟精度，还可以控制直接波和多个波的生成，这有利于提高正演数据集的信噪比。

到目前为止,已经发展了许多单程波正演方法,如隐式有限差分法、Fourier法、傅里叶有限差分(FFD)法和短算子方法等.隐式有限差分法在地震偏移领域已经得到了长足的发展和广泛的应用.它的优点是差分格式稳定性强,能适用于速度任意变化的介质.但是,由于数值频散和传播角度的制约,该方法只适用于传播角度和模型网格都较小的情况.Fourier方法几乎无空间数值频散,计算网格可以很大,具有很高的执行效率,在均匀介质中其传播角度可以精确到90°.但是,Fourier法只能在弱横向非均匀介质中才能得到精确的大角度传播解.作为Fourier法和隐式有限差分的有机结合,FFD法在相移延拓和时移校正的基础上增加了一项具有四阶空间导数精度的有限差分校正项,既保持了隐式有限差分法适应速度横向变化的能力,又继承了相位屏法稳定高效的优势,可以使用相对较大的网格,数值频散也比较弱,在诸如盐丘模型、高陡倾角等复杂构造中取得了很好的应用效果.

将FFD法拓展到三维勘探领域的主要障碍就是三维隐式有限差分所带来的庞大计算量.在实际处理中,通常采用双向分裂方法(Two-waysplitting,也称为交替方向法)将三维隐式有限差分沿着x和y方向进行分裂得到两个串联的二维FFD算子.尽管双向分裂法具有很高的执行效率,它同时也会导致十分明显的分裂误差(也称为方位各向异性,azimuthalanisotropy),严重地制约了FFD法在三维介质中的应用.

为了消除双向分裂引起的误差,人们提出了很多应对措施,但是,这些方法都要以成倍增加的计算量为代价.受到交替方向加插值(ADIPI)方法的启发,我们保留了常规双向分裂中忽略的交叉项,并通过波数域的波场插值消除了常规双向分裂引起的分裂误差.有限差分是在频率空间域实现的,而ADIPI格式中的插值部分却需要在波数域实现,因而利用FFD法在双域实现的结构特点可以进一步改进相应的算法.

基于ADIPI格式的FFD法既消除了分裂误差,又不影响x方向和y方向的精度,而且保持了双向分裂FFD法的高效性,因此它可以胜任高陡倾角、强横向变化介质的三维正演模拟,成为精度高、效率高、适应性强的三维单程波正演方法.

2基于api格式的ffd方法

2.1基于相移算子的相位校正算法

横向非均匀介质中单程波方程的三维FFD算子为

kz≈k0z+ωΔs+b(?2?x2+?2?y2)1+a(?2?x2+?2?y2)?(1)

其中ω是角频率,k0z=√ω2/v20+?2/?x2+?2/?y2是背景速度中的垂直波数,Δs=1/v-1/v0是慢度扰动,v为波场传播速度,v0为均匀背景介质的参考速度,a=0.25(v2+vv0+v20)/ω2,b=0.5(v-v0)/ω.(1)式右端的第一项是相移算子,在波数域完成参考速度的相移;第二项是时移算子,用于校正主传播方向上由于引入参考速度引起的慢度扰动误差,在空间域进行;第三项是有限差分校正算子,用于校正垂直于主传播方向上的相位误差,也是在空间域进行.前两项构成了著名的相位屏算子,可以完成大部分的相位校正,在小倾角和弱横向变化介质中得到了广泛应用.

在某一延拓步长Δz范围内,先将波场P(x,y,z;ω)通过快速傅里叶变换(FFT)变换到波数域,记为～Ρ(kx?ky?z;ω),然后利用(2)式在波数域完成相移:

～Ρ′(kx?ky?z+Δz；ω)=e±ik0zΔz～Ρ(kx?ky?z；ω)?(2)

式中符号±分别对应反向和正向传播.将延拓后的瞬时波场～Ρ′(kx?ky?z+Δz;ω)反变换回空间域,记为P′(x,y,z+Δz;ω),并将其作为时移校正项的输入波场,可以得到

Ρ(x?y?z+Δz；ω)=e±iωΔsΔzΡ′(x?y?z+Δz；ω).(3)

最后,将瞬时波场P(x,y,z+Δz;ω)代入(4)式进行有限差分校正,