SRME原理.docVIP

4 自由界面多次波预测与压制(SRME原理) 利用f-x域波场外推理论，研究在不考虑自由界面反射的情况下，有效波波场的传播及反射规律；然后，在考虑自由界面反射的情况下，详尽推导了自由界面多次波产生的数学模型，给出了任意阶自由界面多次波波场模拟的级数形式和迭代形式，推导作为自由界面多次波产生的逆过程的自由界面多次波消除的数学模型，深刻剖析自由界面多次波的产生和消除的物理机制。 4.1 有效波波场（无自由表面多次波）的模拟 根据第三部分的理论推导，地震波在介质中的传播和反射过程可以用f-x域波场延拓来表示[82-84]，其几何表示如图4-1所示。 图4-1 地震波的传播和反射的物理模型 用系统理论可以表示为图4-2。它把地震波场分解为五个过程:发射，下行传播，反射，上行传播，接收。 图4-2 地面波场的正演模拟（无自由界面的反射） 自由表面 根据Berkhout的研究，对于地下的某一个反射界面m，其对应的有效波波场（f-x域）为[108] （4-1） 如果忽略自由反射界面，则。 无自由表面接收到的总的有效波波场[111]为 代入式（4-1），得 （4-2） 进而 （4-3） （4-4） 式中，算子为自由界面无反射的地下介质的脉冲响应； 算子为总的上行波场； 算子为总的下行波场； 算子为在假设自由界面不反射情况下，在自由界面处的上行波场； 算子为在自由界面处，对于下行入射波场的反射算子； 算子为在自由界面处的下行震源波场，包含激发方式和震源的虚反射效应； 算子为上、下行传播矩阵，每一列对应于地表上一个点的传播算子，整个矩阵描述了波场从到的传播。对于算子的每一元素可以表示为 式中，限带微分算子，为从炮点到反射点的旅行时，为从反射点到接收点的旅行时，球面扩散因子； 算子为界面处的下行入射波场的反射算子。 无自由表面多次波的时间域表达式为[112] （4-5） 式中为kronecker内积。 4.2 自由界面多次波正演模拟 图4-3 有效波场和多次波场的正演模拟，包含自由界面的反射 当不考虑自由界面的反射时，在界面处的上行波场为 （4-6） 当考虑自由界面的反射时，上行波场遇自由界面后发生反射，变为下行波场。因此总的下行波场不仅包含下行的震源波场，而且还应包含被自由界面反射回地下的波场:，如图4-3所示。为自由界面对于上行波场的反射算子，为自由界面处的上行波场。因此自由界面处的总的上行波场为: (4-7) 在实际地震资料中，还应该考虑检波器的特性(例如检波器组合模式、虚反射特性等)，将检波器属性定义为算子，则在自由界面处波场为: (4-8) 将代入式（4-8），得 (4-9) 可以将自由表面接收到的波场分解为有效波场和（自由界面）多次波波场，即 (4-10) (4-11) 由式（4-7）-（4-11），得 (4-12a) (4-12b) Berkhout称为自由界面算子。 由代数运算可得 （4-13） 用Neumann级数展开得[112] ＝ (4-14) 由通项可知： 当n＝0时，，即为1阶自由界面多次波； 当n＝n时，，即为n阶自由界面多次波。 自由界面多次波波场模拟的方法可以分为以下两种。 截断离散级数方法 (4-15) 级数的次数就是多次波的阶数。首先模拟出有效波场，然后叠加上多次波场，就得到总的波场。具体展开可得 当n＝1时， 当n＝2时，+ 当n＝3时，++ …… 当n＝n时，++ 在实际的多次波模拟中，选取有限项。 迭代方法 由，其迭代表达式为 （4-16） 其中，n为迭代次数，是第n+1次迭代得结果，是有效波波场。 从本质上讲，迭代算法和级数算法是一致的，只不过两者具体的实现形式不一样；迭代形式的表达更加简洁，物理含义更加明确(为自由界面多次波的预测算子，为自由界面算子)，在计算机上更易实现，实现起来也能节约计算机资源。 4.3 自由界面多次波的压制 （1）多次波迭代压制理论基础[106] 自由界面多次波的消除过程是自由界面多次波产生过程的反过程。 (4-17a) (4-17b) 经过矩阵变换： = （4-18） 由于原始数据中