地震第3章_反褶积（一） .pptVIP

Chapter3 反褶积; 反褶积处理是褶积处理的反过程，因而称为反褶积。在前面第一章第一节中，我们曾讲过，一个滤波器的滤波过程在时间域的输出是输入信号与滤波器滤波因子的褶积。因此，时间域的褶积处理就相当于一个滤波过程。而反褶积则相当于时间域的一个反滤波过程。地震记录可以看作是地震子波与地层脉冲响应的褶积，即;(3一1)式可视为一个滤波过程，如图3-1所示。; 如果我们设计一个滤波器，其滤波因子 具有与滤波器 约恰好相反的性质，即当输入为地震道记录 约时，其输出为 地层脉冲响应 。我们称这个反过程为反滤波或反褶积，如 图3-2所示。;二、褶积模型; (3-1)式所表示的地震道记录 是不包含随机环境噪声的。合成地震记录，则需要在(3-l)式中加入适量的随机环境噪声 。 为了得到更真实的人工，这时，人工合成地震记录褶积公式为;式中 分别为 和 的频谱。;式中。 —震源脉冲值，为一常数; —反射界面的反射系数。 但是，由于地层介质具有滤波作用，这种大地的滤波作用相当于一个滤波器。因此，由震源发出的尖脉冲经过大地滤波器的滤波作用后，变成一个具有一定时间延续的波形 ，通常叫作地震子波(图3一6)。这时，地震记录是许多反射波叠加的结果，即地震记录 是地震子波 与反射系数 的褶积;实际的地震记录城 除了(3一4)式所表示的一系列反射波 而外，还存在着干扰波 ，因此，地震记录双 的一般模型可以写为; 在普通的地震记录上，一个界面的反射波一般是一个延续时间为几十毫秒的波形。由于地下反射界面一般是相距为几米至几十米的密集层，他们的到达时间差仅为几毫秒到几十毫秒，因此，在反射地震记录上它们彼此干涉，难于区分开来。为了提高反射地震记录的分辨能力，我们希望在所得到的地震记录上，每个界面的反射波表现为一个窄脉冲，每个脉冲的强弱与界面的反射系数的大小成正比，而脉冲的极性反映界面反射系数的符号。那么，怎样把延续几十毫秒的地震子波 ，压缩成为一个反映反射系数 的窄脉冲呢?这就是反褶积所要解决的问题。 如果，地震记录是(3-4)式所表示的地震子波 与反射系数 的褶积，即地震记录中只有反射波 ，而没有干扰波 。这时反褶积问题很简单。 根据(3一4)式，在频率域相应有;式中， 分别是地震记录 、地震子波 和反射系数 的频谱。 显然;在时间域，得到; 因而，所谓的反褶积或反滤波实际就是一个滤波过程，只不过是这种滤波过程其作用恰好与某个滤波过程的作用相反。; 1.直接观测法 这种方法是用专门布置在震源附近的检波器直接记录地震子波 ，此方法只适用于海上地震勘探。 在某些地区的海上地震勘探中，在地震记录上海底反射波到达之前曾记录到一个地震波。经过分析知道这是由于海水含盐量有分层性所形成的。由于海水的含盐量有分层性使海水明显地分成上下两层。下层的含盐量较上层含盐量高，形成了一个较为清楚的界面。由震源出发的地震波到达这个界面引起反射返回到海面下的检波器，被记录下来。由于这个波没有与其他波干涉，所以可以作为地震子波 。使用这样求取的地震子波进行反褶积，得到了良好的效果。;假设反射系数 为白噪声，其z变换为 ，则 自相关 的Z变换;将(3-13)式代入上式，得到;因而，由(3一1)式得到;地震子波， 还满足稳定性条件，即;因而得到;其中，M要求取的足够大。即可求出地震子波 。 如果地震子波 不是最小相位的，而是零相位的。假设反射系数 为白噪声，则其振幅谱;再计算其频谱;上式也可由(3一1）式中，令 ，直接得出。 将所得到的地震子波的频谱 进行傅里叶反变换，得到;假设反射系数 白噪声，即;即地震记录 的自相关函数与地震子波 的自相关函数相同;将上式两边同乘以 ，得到;则(3-35)式可以写成;由此得出最小相位的地震子波;其中 ，得到随反射时间变化的声阻抗曲线 。然后利用反射系数公式;最后，对 进行傅里叶反变换，就得到地震子波; 为了将地震子波 和反射系数 从地震记录 中分离开来，对上式等式两端取对数得到; 上式表明，地震记录