物理03反褶积课件.ppt

混和相位子波 反射系数：受子波相位和滤波门的影响，反射系数反极性并有延迟 例1(续) 只有一个非零尖脉冲时，V＝1，达到最大值； 具有两个相等的非零尖脉冲时，V＝0.50； 具有3个相等的非零尖脉冲时，V＝0.33； 具有10个相等的非零尖脉冲，V＝0.10。 因此，记录道中相等的非零尖脉冲的数量越少，V的数值就越大。当记录道段中含有多个不相等的非零尖脉冲时，这些尖脉冲的振幅相差越大，V的数值也就越大。因而，使V值达到最大，有着简化记录道段外形的作用。V值的大小不受尖脉冲的极性及其所在位置的影响 反褶积－预测反滤波 四、利用预测反滤波消除海上鸣震 1．海上鸣震的数学模型 所谓海上鸣震，就是由海底面下面的反射界面产生的反射波返回到海表面后，在海表面和海底面之间产生的多次反射波。这种由一次反射波s(t)(可看作新波源)再在海表面和海底面之间多次(往返)反射的多次波就叫做海上一级鸣震(reverberation)，又叫做混响 反褶积－预测反滤波 设：海水深度为H，波速度为v0海底面的反射系数为R，海表面的反射系数为-1。假设各次多次反射波的波形与一次反射波s(t)的波形相同它们之间只有时间延迟。这样就可以得到含有一次波和各次多次波的地震记录。 是地震波在海水中垂直往返的旅行时间。 FFT 反褶积－预测反滤波 海上鸣震相当于滤波作用 滤波器的振幅频率特性曲线 反褶积－预测反滤波 反褶积－预测反滤波 FFT 反褶积－预测反滤波 2．利用反滤波消除鸣震 但是海水底面的反射系数R往往是不知道的，实际应用应用受限 反褶积－预测反滤波 3．用预测反滤波消除海上鸣震 是地震一次反射波 地震反射波s(t)可以看做是地震子波b(t)与反射系数 的褶积 IFFT 反褶积－预测反滤波 取预测步长 ，有： 在 ，有： 反褶积－预测反滤波 预测值 预测误差 用预测反滤波确实可以消除海上鸣震等多次波干扰。 一次波 四、参数选择 预测反褶积程序中的几个主要参数选择如下: 预测步长 ；是预测反褶积的关键参数，预测反褶积的效果与预测步长 的选择有很大的关系。 为了选择最佳预测步长，可以将叠加前地震记录作为输入，进行预测反褶积试验，在其他参数（例如预测因子长度和预白化量）不变的情况下，改变预测步长 ，从不同预测步长的预测反褶积输出记录直接确定最佳预测步长 。 选择最佳预测步长 ，也可以用人工合成的地震记录作为作为输入，进行预测反褶积试验，同样，在其他参数不变的情况下，改变预测步长，从预测反褶积输出结果确定最佳预测步长 。 如图10所示，预测因子长度 m 不变为128ms，预白化量e为0%，采样间隔为2ms，取预测步长 分别为2ms，8ms，20ms，32ms，44ms，和88ms。图中（a）为脉冲相应，（b）为合成地震记录，地震子波为已知最小相位子波。从图中可以看出预测步长 = 2ms时，预测反褶积效果最好。 图10 用合成地震记录试验预测步长 （2）预测因子长度m；对预测反褶积的结果也有明显的影响。总的来说，预测因子不能太小。长度太小预测反褶积波形尾部出现明显波动，长度增大，这种尾部波动减弱，反褶积效果较好。达到一定长度时，反褶积效果将趋于稳定。预测因子长度也不能太大，长度太大不仅增加计算量，而且会增加相邻反射及噪声的影响。 图11所示为用合成地震记录确定预测因子长度的试验，预测步长 不变为2ms，预白化量 为0%，取预测因子长度分别为20ms，44ms，94ms，128ms，192ms，292ms，图中（a）为脉冲响应，（b）为合成地震记录，地震子波为已知最小相位子波。从图中可以看出。当预测因子长度太小（例如20ms和44ms）时，预测反褶积波形尾部波动明显，当预测因子长度达到 =128ms时，预测反褶积效果已经较好。 （3）预白化量 ：在进行预测反褶积计算预测因子时，要根据（87）式求解地震记录相关矩阵。为了求解矩阵方程的稳定性，需要根据随机噪声干扰水平确定随机噪声自相关值e或者它与地震子波自相关的百分比。 图12所示为合成地震记录确定预白化量得试验，预测步长a为2ms，预测因子长度m为128ms，预白化量分别为0,0.1%，1%，5%，10%和20%，图中（a）为脉冲响应，（b）为合成地震记录，地震子波为最新相位子波。从图中可以看出，当预白化量 =0.1%时，预测反褶积的输出结果基本不变，而当预白化量很大（ 5%），预测反褶积波形尾部波动明显增大，预测反褶积效果明显