多次波压制技术与应用(1-2)

② 速度谱上识别：多次波的速度一般都小于相同位置处的一次 反射的速度；
③ 速度扫描识别：用一系列速度对叠加剖面进行常速扫描也是 识别多次波的一种有效手段。如果存在多次波，那么多次波 在剖面上相同位置，不同速度下的成像和形态将会有较大的 差别；
④ 频谱识别方法：多次波，特别是鸣震，某些单频能量很强；
2019/12/30
② 如果多次波在时差上满足周期性，理论上来说，用预 测反褶积是可以精确压制多次波能量的。叠前道集， 多次波时差随偏移距增大变小，步长选取以远道的自 相关次极值为好。
③ 预测反褶积是子波截尾反褶积，当预测步长小于子波 长度时，无论是一次波还是多次波都会子波都会被切 除。从这个角度来看，预测反褶积去多次波，对一次 波不是无损伤的，更谈不上是保幅。
NMO之前的一次波 与多次波几何形态
NMO之前的一次波 双曲时差曲线
NMO之后的一次波 与多次波几何形态
NMO之后的一次波时差曲 线，被拉平为一直线
沿抛物线计算一次 波与多次波模型
仅在这一抛物线范 围内计算一次波模 型
2019/12/30
NMO之前的多次波 双曲时差曲线
NMO之后的多次波时差曲 线，为抛物线形态

慢度: u 1/ v z
2019/12/30
下行波旅行时间
t(z) zq z cos / v z 1 p2v2 / v
Tau-p域P道集周期性推导
一次反射波
t p (z) 2zq 2z cos / v 2z 1 p2v2 / v t0 1 p2v2
一阶多次波 tm(1) (z) 4zq 2t0 1 p2v2
在该抛物线 范围内计算 多次波模型
Radon法压制多次波原理
NMO校正
2019/12/30
v
sin2
sin2
t0

sin m sin
t0
当地层倾角较小时: t0 mt0 周期性、可预测性
水平地层、近炮检距或叠加剖面
2019/12/30
预测反褶积压制多次波的适用条件
自相关函数
反褶积算子长度 n 预测步长 m
白噪系数
xx (0)


xx
(1)
xx (2)
Tau-p 反褶积前叠加剖面
2019/12/30
Tau-p 反褶积后叠加剖面
① 利用多次波可预测性压制多次波的方法对短周期的全程多 次波压制效果较好，但对长周期的全程多次波和层间多次 波压制却难以奏效。有可能在记录时间长度内没有足够的 多次波满足周期性的要求；另一方面，长周期的多次波需 要长的算子，而在一个长的时窗范围内，甚至一次波也有 可能呈周期性，这样的算子就潜伏着去多次的同时也除去 了一次波的风险。
f（t）
y(t)
y(t) f (t)* x(t)
预测反褶积
f（t）
x(t m) f (t)* x(t)
预测反褶积与数据匹配
x(t m)
确定性反褶积
①水平常速介质；②法向入射；③稳态子波
x(t) w(t) r(t) n(t)
④没有噪声
x(t) w(t) r(t) ( f ) W ( f )R( f )
a)
P
Twoway time (t)
M1
Offset (x)
t-x域非周期性
b)
Horizontal slowness (p)
Tau-p域周期性
P
Inter cept time (t)
M1
t( p)
M2
M2
2019/12/30
Tau-p域预测反褶积的核心思路
垂直慢度:
q cos / v
水平慢度: p sin / v
2019/12/30
基于波动方程的预测减去方法
算法
基本物理单元
附加信息
波场外推法
SRME 反散射级数法
水层 海底
自由界面 层界面
自由界面 点散射
水层深度
数据规则化 数据规则化 背景速度场
2019/12/30
P1
S
P2
B
R
P1 P2 M1-P1 M2-P1 M3-P1
M4-P1
P1 P2 M1-P1
⑤ 分偏移距叠加识别：多次波在不同偏移距上具有不同的动态 时差，因此不同偏移距的叠加剖面形态差别很大，一般在小 偏移距剖面上表现出多次波的形态，而在大偏移距剖面上为 一次反射的形态，据此可以识别多次波；
⑥ 自相关识别多次波：根据多次波的周期性，即多次波在自相 关上表现为子波的续至相位，可以通过自相关来识别和选取 压制多次波的方法和参数；
⑤已知子波
r(t) w*(t) x(t) E( f ) W 1( f )( f )
⑥最小相位子波 子波的逆稳定和最小相位
预测反褶积与数据匹配
2019/12/30
预测反褶积 ①水平常速介质；②法向入射；③稳态子波
x(t) w(t) r(t) n(t)
④没有噪声 x(t) w(t) r(t) ⑤随机反射系数（子波与数据自相关相等）
多次波压制技术与应用
目录
1. 多次波及压制技术简介 2. 基于信号分析的多次波压制技术 3. 基于波动理论的多次波压制技术
2019/12/30
直达波
2019/12/30
直达虚反射
虚反射
一次波
多次波
2019/12/30
按照反射层的位置，多次波可以分为：
一次波
全程多次波
2019/12/30
全程 多次波
2019/12/30
④ 预测算子长度至少要大于多次波的一个周期。算子长 度取小了，会带进假的能量。因此算子应该适当取长 一点。
⑤ 级联预测反褶积效果有待进一步探讨。一般而言，先 做较大步长多次波去噪，再做较小步长子波压缩提高 分辨率是比较妥当的；先小步长、再做大步长，意义 不大。但预测反褶积步长本身不能一分为二，严格来 说，要想通过反褶积串联来实现多步长多次波的去噪 是不可能的。
M2-P1 M1-P2
M3-P1
M2-P2 M4-P1
P1
P2 M1-P1
PL1M1 M2-P1 M1-P2 PL1M2 M3-P1
PL1M3 M2-P2 M4-P1
P1
P2 M1-P1
PL1M1 M2-P1 M1-P2 PL2M2 PL1M2 M3-P1 PL2M2 PL1M3 M2-P2 M4-P1
按照反射层的位置，多次波可以分为：
层间多次波
层间多次波
2019/12/30
按照反射层的位置，多次波可以分为：
微屈多次波
2019/12/30
海底一次波 海底多次波 微屈多次波 微屈多次波
按照反射层的位置，多次波可以分为：
2019/12/30
鸣震
按照延迟时间的长短，多次波可以分为：
一次波
长程多次波
2019/12/30
2.1 预测反褶积压制多次波技术 2.2 Tau-p域反褶积压制多次波技术 2.3 F-K滤波压制多次波方法 2.4 Radon变换压制多次波技术 2.5 聚束滤波压制多次波技术 2.6 内切与加权叠加
2019/12/30
x(t) x(t)
2019/12/30
数据滤波和数据匹配
2019/12/30
2.1 预测反褶积压制多次波技术 2.2 Tau-p域反褶积压制多次波技术 2.3 F-K滤波压制多次波方法 2.4 Radon变换压制多次波技术 2.5 聚束滤波压制多次波技术 2.6 内切与加权叠加
2019/12/30
排列内Βιβλιοθήκη 地层水平ox
x
p
2019/12/30
t
Tau-p域预测反褶积压制多次波的适用条件
J ( (t) w(t)*w(t)1)2 dt min
2019/12/30
预测反褶积与数据匹配
输入
期望输出
求解
数据匹配
数据x
数据y
匹配算子
预测反褶积 压制虚反射 压制虚反射
输入数据或子 脉冲（反射系 反子波算子
波
数）
子波+虚反射 脉冲（反射系 反子波和虚反
数）
射算子
子波
虚反射
预测算子
⑦ Radon域：一次波和多次波的聚焦点不同； ⑧ VSP记录：一次波和多次波的运动轨迹有差别。
2019/12/30
2019/12/30
动校正道集
速度谱
速度扫描V=2400米/秒
2019/12/30
速度扫描V=4000米/秒
0-4500米炮检距叠加
2019/12/30
1500米-4500米炮检距叠加
2019/12/30
反射波与多次波的NMO特征
2019/12/30
2019/12/30
F-K滤波压制多次波原理
反射波
多次波
2019/12/30
NMO校正后的反射波与多次波F-K谱
2019/12/30
叠加速度对F-K滤波压制多次波的影响
2019/12/30
NMO校正后
NMO校正+FK滤波后
2019/12/30
② Tau-p域的预测反褶积效果优于 t-x域。 ③ Tau-p域的预测反褶积的关键是要有足够的转换精度。
2019/12/30
2.1 预测反褶积压制多次波技术 2.2 Tau-p域反褶积压制多次波技术 2.3 F-K滤波压制多次波方法 2.4 Radon变换压制多次波技术 2.5 聚束滤波压制多次波技术 2.6 内切与加权叠加
2019/12/30
多次波
按照延迟时间的长短，多次波可以分为：
短程多次波