震源子波matlab,VSP震源子波应用

VSP资料[1,2]采集时,每激发一炮只能测量一个深度段,并在距炮点很近的地方埋置检波器记录震源子波。目前震源子波信息主要用途是利用子波初至起跳时间来监控激发井深变化及子波一致性评价[3,4],并没有充分利用震源子波的信息。目前的VSP资料处理中,当炮与炮间存在激发能量差异时,没有合理的方法用于激发能量消除。提高分辨率处理的反Q滤波[5]所需的Q值难以求准,大大影响实际应用效果;反褶积[69]过程中,由于受地层的吸收造成子波时变,不满足反褶积的假设前提,处理效果也很难尽如人意。本文利用震源子波消除记录中的激发能量差异;同时按照波动方程所描述的地震波传播规律,结合广义S变换[1014],利用震源子波作为参考来消除地层吸收衰减的影响,使所采集资料获得更高的分辨率。1利用震源子波消除激发能量差异零偏VSP资料采集中,震源子波监控点一般位于震源510m处。因此,可忽略地层吸收衰减差异,认为监控子波能量即为激发能量。对监控子波能量进行归一化,即可消除VSP记录中炮与炮间的激发能量差异。1.1方法原理设井下检波器深度为第i米时记录的地震记录为hi(t),监控的震源子波为wi(t)。则激发能量可表示为Ei=Lt=0wi2(t)(1)式中,L的选择视资料的品质而定,一般选择1000个样点即可。对监控子波wi(t)的能量进行归一化,得wi(t)=wi(t)/Ei(2)消除激发能量差异后的地震记录可表示为hi(t)=hi(t)/Ei(3)1.2实际资料应用图1为监控的震源子波,图2为原始零偏VSP记录。图1和图2中,在相同的深度,能量具有明显的一致性。求取每一个深度处震源子波能量(图3),结合图2可以看出:图2中的红色方框对应的区域具有68西南石油大学学报(自然科学版)2011年图1震源子波Fig.2TheSourcewavelet图2原始记录Fig.2Rawdata明显的激发能量差异,最大能量差异可达34倍。对监控子波进行能量归一化(图4),并按照式(3)对零偏记录进行激发能量处理后,即可消除激发造成的记录中的能量差异(图5)。对比图2和图5,除红色区域内地震记录振幅有明显变化外,其他区域没有发生明显变化;较图2来讲,图5的资料品质得到明显改善,消除了由于激发能量差异带来的解释陷阱,有利于后期的处理和解释。图3各道监控子波能量Fig.3Theenergyofwaveletineachtrace图4能量归一化后的监控子波Fig.4Theenergyofwaveletafternormalization图5消除激发能量差异后的零偏VSP记录Fig.5VSPdataafterenergynormalization2利用震源子波消除地层吸收衰减的影响将监控子波作为参考信号,利用广义S变换将地震记录的分辨率提高到子波的分辨率。2.1吸收衰减补偿原理基于GST变换的吸收衰减补偿[1518]基本思想是:如果没有地层吸收,深、浅层反射波具有相同频第6期张固澜,等:VSP震源子波应用69宽的振幅谱,相位谱仅相差一个线性相位;如果把地震记录分成不同的频率,所对应时间的能量分布关系具有相似性。但是,由于地层的吸收作用,造成各频率能量对时间的分布不同,给地层乘以时变因子使它们相同,就起到了对地层吸收的补偿作用。假定散射、干涉、多次反射造成的振幅衰减都已得到了合理的补偿,只需考虑地层吸收、频散以及球面扩散。设地层反射系数为rk(k=1,2,,K),由于子波时变,每个反射系数会受不同的Q滤波,对应时间Tk,经Q滤波后的振幅响应为X(Tk,)=rkp(0,)exp(Tk2QeqTk)exp(j