反褶积2016课件.pptVIP

Seismic Processing 地震数据处理 任课教师：李宇 职 称：讲师 所属院系：地测学院地球物理系 电 话邮 箱：liyupa@126.com = * f A 1 f A = f A 反褶积 第四章 反褶积 = * f A = f A f A 1 第四章 反褶积 f A f 反褶积 第四章 反褶积 = * f A = f A f A 1 求反射系数，即滤波算子对地震记录的滤波，相当于地震子波对反射系数滤波的反过程（反滤波）。 滤波算子是地震子波的逆。 滤波算子对地震记录进行滤波，得到反射系数。 地震子波对反射系数进行滤波，得到地震记录。 第四章 反褶积 反滤波是一种反褶积方法； 反滤波是一种确定性反褶积方法。 反滤波已地震子波已知为条件； 反滤波 地震子波对反射系数滤波的反过程(反滤波)。 = * 第四章 反褶积 地震子波的z变换是一个z的多项式； 变量z的幂次等于子波时间序列单位时间的延迟量； z的系数是子波时间序列的样点值。 0延迟 1延迟 1 -1/2 幂次 系数 1 -1/2 反滤波/Z变换 第四章 反褶积 反滤波/Z变换 时间序列f(1,-1,1/2)与时间序列g(1,1/2,-1/2)褶积 1 -1 1/2 延迟 输出 -1/2 1/2 1 -2 1 -1/2 1/2 1 -1 -1/2 -1/2 1/2 1 0 -1/2 -1/2 1/2 1 1 3/4 -1/2 1/2 1 2 -1/4 第四章 反褶积 反滤波/Z变换 0延迟 1延迟 2延迟 1 -1 1/2 幂次 系数 1 -1 1/2 0延迟 1延迟 2延迟 1 1/2 -1/2 幂次 系数 1 1/2 -1/2 第四章 反褶积 反滤波/Z变换 第四章 反褶积 滤波算子 无限长，截短滤波算子。 反滤波/Z变换 第四章 反褶积 反滤波/Z变换 地震记录x(1,-1/2)与滤波算子f(1,1/2)褶积 1 -1/2 延迟 输出 期望输出 1/2 1 -2 1 1 1/2 1 -1 0 0 1/2 1 -0 -1/4 0 地震记录x(1,-1/2)与滤波算子f(1,1/2,1/4)褶积 1 -1/2 延迟 输出 期望输出 1/4 1/2 1 -3 1 1 1/4 1/2 1 -2 0 0 1/4 1/2 1 -1 0 0 1/4 1/2 1 0 -1/8 0 第四章 反褶积 反滤波/Z变换 子波 (1, -1/2) 滤波输出 期望输出 误差 子波的逆 (1, 1/2, 1/4, 1/8······) 滤波算子 (1, 1/2) (1, 0, -1/4) (1, 0, 0) 0.063 滤波算子 (1, 1/2,1/4) (1, 0, 0, -1/8) (1, 0, 0, 0) 0.016 滤波算子系数越多，输出越接近反射系数 1。 误差L定义为滤波与期望输出的系数之差的平方和： 第四章 反褶积 Z 页岩 2.4 2.3 含气砂岩 2.2 1.8 含水砂岩 2.5 2.3 页岩 2.5 2.4 反射系数R与透射系数T 第四章 反褶积 Z 页岩 2.4 2.3 5.52 含气砂岩 2.2 1.8 3.96 含水砂岩 2.5 2.3 5.75 页岩 2.5 2.4 6.00 反射系数R与透射系数T 第四章 反褶积 1. ，产生反射。即反射波条件：波阻抗差异 2. ，差值越大，反射波振幅越大 3. ，反射系数为正，反射波与入射波相位相同 ，反射波与入射波相位相反 4. 恒成立，透射波与反射波相位相同 入射波:Incident wave. 反射波:Reflection wave. 透射波:Transmission wave. 反射系数:Reflection coefficient. 透射系数:Transmission coefficient. 反射系数R与透射系数T 第四章 反褶积 反射系数序列 反射系数序列：反射系数有序的排列。 有序排列：时间域按照时间排序，深度域按照深度排序。 第四章 反褶积 Ricker wavelet /雷克子波 地震子波(Seismic Wavelet) 40ms 第四章 反褶积 40ms 地震子波(Seismic Wavelet) 第四章 反褶积 40ms 64ms 地震子波(Seismic Wavelet) 第四章 反褶积 40ms 20ms 地震子波(Seismic Wavelet) 第四章 反褶积 反褶积是通过压缩地震记录中的地震子波