地震数字处理-1de wolf、薄板和屏近似.pptVIP

*利用DeWolf近似和相屏传播算子计算反向散射波基于DeWolf近似和相屏单程传播，给出了一种在地面反射地震观测中计算一次反射波或反向散射波的合成地震记录的方法。具体实现时有两种不同的算法：一种是直接利用MFSB近似，其中计算反向散射波的步长等于网格间隔，但计算正向传播的步长使用远远大于网格间隔的屏间隔；另一种是对正向和反向散射波的屏近似，这样两种步长都能采用屏间隔，从而大大提高了计算效率。在复杂不均匀介质（尤其是3-D介质）中，快速模拟方法和算法是发展复杂构造成像和反演方法、地震方法解释和用于复杂构造的关键。有限差分和有限元法从原理上能模拟在任意不均匀介质中波的传播，但计算工作量大、效率低。在大型3-D弹性波问题的情况下更是不可想象。下面我们给出了一种基于多次正向散射-单次反向散射（Multiple-forescatteringsingle-backscatteringMFSB）近似，即DeWolf近似，在地面反射地震观测中计算反向散射场的新方法。为快速实现该方法，我们推出了双域算法。当不均匀体尺度大于地震主波长时，这一理论可进一步通过薄板近似和屏散射近似得到合理近似和简化。屏近似能减少计算工作量、提高计算效率。数值模型试算验证了所提方法的有效性。一．多次正向散射-单次反向散射（MFSB）近似---DeWolf近似当介质体内的不连续性不是很强或不均匀体的参数扰动不太大时，在不均匀体和谐振散射之间的混响通常可以忽略不计。但正向散射的累加效应一般不能忽略掉。实际上，对于大型不均匀体介质或长传播距离，多次散射对正演模拟和反问题来讲是很重要的。在这种情况下，Born近似不满足，但可应用DeWolf近似。1．Lipmann-Schwinger方程从标量波动方程出发（1）定义（2）得到（3）由地震叠加原理，令（4）那么，基于Green’s函数，有（5）其中，是背景介质中的Green’s函数。上式就是Lipmann-Schwinger方程。2．散射序列的重新归一化和DeWolf近似Lipmann-Schwinger方程有一多次散射Born序列的标准解。广泛应用的Born近似是该序列的第一项。Born近似仅当非匀质体较弱和传播距离较短时适用。在对多次散射序列重新归一化后，DeWolf推出了MFSB近似：（6）其中，和分别是重新归一化的多次正向散射场和Green’s函数。利用相屏算子可计算和。值得注意的是：DeWolf近似用重新归一化的多次正向散射近似代替总的实际波场和自由空间的Green’s函数，优于由Wu和Huang给出的近似，它是仅近似总的实际波场，而自由空间的Green’s函数保持不变。MFSB近似适用的条件是反向散射场远远小于正向散射场。3．利用相屏传播算子实现MFSB合成记录利用（6）式可计算反向散射场（7）其中，是在深度z处接收器的水平位置。方程（7）利用相屏传播算子可由数值法实现。为提高反向散射场的计算效率，在一薄板内可使用局部Born近似。这意味着在薄板内正向散射场可保持不扰动，可由常速介质中的Green’s函数来替代。对方程（7）相对做Fourier变换，得到