深水陡坡带绕射多次波压制方法研究.docxVIP

深水陡坡带绕射多次波压制方法研究 1 层间绕射多次波的压制 自由地表多个波的抑制已经成为海洋地震资料处理的传统处理过程。 绕射多次波通常是经海底一次绕射以后导致波的传播路径发生畸变, 畸变后的情况可以分为三种:第一种, 畸变后波的传播方向仍然在原先设定的电缆长度的记录范围之内, 经自由表面反射产生的多次波;第二种, 畸变后波的传播方向超出了原先设定的电缆长度的记录范围, 但是仍然由自由表面反射产生的多次波;第三种, 无论是否超出了预先设定的电缆长度的记录范围, 畸变后波的传播方向发生变化, 在海底地层发生多次反射形成层间绕射多次波;第一种情况下的绕射多次波可以归纳到传统的SRME预测和消除多次波方法当中予以压制;第三种情况因为属于层间绕射多次波, 不在本文方法的讨论范围之列;但是对于第二种情况, 常规SRME方法由于电缆长度的限制, 在深水陡坡带和海底构造异常复杂的情况下局限性尤为明显, 绕射多次波的压制也愈加困难.本文借助炮点、检波点互换原理, 将邻炮的炮记录依照互换原理引入到反馈迭代循环当中, 补充绕射多次波的近偏移距信息不足和负偏移距信息缺失等问题, 修正深水陡坡带因为电缆长度限制而造成的信息缺失, 最大限度地准确预测深水陡坡带的与自由表面有关的绕射多次波, 也可以给类似多次波偏移成像方法 2 反馈迭代模型 在2D情况下, 从波动方程出发可以推导出Kirchhoff积分, 由Kirchhoff积分可以推导出Rayleigh积分, 这里由第二型Rayleigh积分可以知道地震波的传播过程对应空间褶积模型的具体表达式如下 式中, p (x, z 由 (2) 式可知, 地震波的传播过程可以看作是一个大地滤波过程, 表示的是地震信号在地层中的传播和反射效应.由反馈迭代模型 上式的物理解释可以用图1说明. 在z 联合 (2) 式和 (3) 式有: 与自由表面相关的绕射多次波遵循多次波传播的一般规律, 由于单边海上拖缆数据采集的拖缆长度限制, 在复杂海底地质构造情况下, 炮记录上记录到的绕射多次波常常超出x 针对绕射多次波的产生机理, 在现有单边海上数据采集基础上, 通过将邻炮的炮记录依照炮、检互换原理引入到反馈迭代循环当中 (如图3中虚线部分地震记录信息) , 补充绕射多次波的近偏移距信息不足和负偏移距缺失, 则对公式 (6) 、 (7) 做如下改进; 其中, S 3 应用分析的例子 3.1 正演模拟记录 图5a为SMARRT发布的Sigsbee2b模型, Sigsbee2b是墨西哥湾深水模型, 由于该模型构造复杂、绕射点多, 故而绕射多次波异常发育.该模型划分为3200×1200个网格, 水平方向和垂直吹响网格间距均为7.6m.该模型的观测系统设计有500炮的正演模拟, 最小检波距为0, 最大炮检距约为3900m, 炮距为检波距2倍, 前348炮为每炮348道, 往后每炮依次减少2道; 图5b为该模型的第251炮正演记录, 由该正演模型记录可以发现时距曲线为标准双曲类型的自由表面多次波相当发育, 但由于盐丘和绕射点的原因, 时距曲线为非双曲类型的多次波也相当发育, 而且能量相当强;采用同样的均衡多道二范数匹配相减 3.2 深水斜坡带绕射多次波压制 该数据为南海深水盆地某二维海洋地震叠加剖面, 其海水深度变化大, 构造复杂且陡坡带绕射多次波较为明显, 常规的自由表面多次波消除方法不能较好的解决深水陡坡带的绕射多次波压制;采用改进的反馈迭代方法以后, 可以更为准确地预测和匹配压制深水陡坡带绕射多次波, 通过均衡多道匹配滤波以后可以得到如下所示结果, 图6a为自由表面绕射多次波压制前的叠加剖面, 图6b为自由表面绕射多次波压制后的叠加剖面, 从两个椭圆标注出的范围不难发现由于陡坡带和崎岖海底产生的绕射多次波得到了较好的压制, 被多次波掩盖的一次波能量得到显现. 4 相邻炮记录的对比分析 深水陡坡带和复杂海底的绕射多次波受到采集拖缆的长度限制, 由陡坡带和复杂海底产生的绕射多次波压制方法是海上资料处理的一个难题, 本文在反馈迭代模型的基础上, 依据炮点检波点互换原理, 引入相邻炮记录的相关信息, 补充了由于绕射而导致炮记录的近偏移距信息不足和负偏移距信息缺失问题, 通过Sigsbee2b模型和实际资料的压制结果对比分析, 较好地证明了该方法的有效性.另外, 海上拖缆采集的方向也和自由表面绕射多次波的预测也有较大影响, 采用不同方向采集的拖缆数据综合处理应该会效果更好. 在频率域将 (4) 式进一步简化有: 其中, 式中, x 该式的物理意义可以理解为叠前炮记录上的自由表面相关多次波可以由炮点检波点之间的自由表面多次波和当前炮的最小最大偏移距之外的绕射多次波之和, λ为叠加系数因子, λ的选择需要依据实际资料而定.如图4所示, P (