震源子波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

震源子波matlab,VSP震源子波应用相关推荐

matlab 子波,基于Matlab的Bark子波实现
0引言子波的种类较多,根据是否正交,子波可以分为正交子波.半正交子波.双正交子波和非正交子波.通常,在多分辨率框架下构造子波基的主要优点在于它满足正交性,然而这一点对语音来说并不必要,因为听觉感知本身 ...
震源子波matlab,matlab 我们地震处理相关的一些程序，包括雷克子波，segy的读动校正这些 272万源代码下载- www.pudn.com...
文件名称: matlab下载收藏√ [ 5 4 3 2 1 ] 开发工具: Others 文件大小: 27355 KB 上传时间: 2017-03-14 下载次数: 0 提供者: 王 ...
ricker 子波matlab,毕设求助！！！雷克子波合成地震信号
我在运行下面代码的时候碰到了这个错误 ,请问怎么修改啊未定义与 'double' 类型的输入参数相对应的函数 'wigb'. 出错 monidizhen (line 66) wigb(R) clea ...
matlab求雷克子波相位谱,求雷克子波的振幅谱和相位谱的MATLAB程序，谢谢！！！！！...
%% 自己写了如下程序,望采纳! clc;clear;close all; %Ricker子波 f=40;dt=0.004; nw=6./f/dt; nw=2*floor(nw/2)+1; nc=fl ...
matlab画雷克子波
clear; clcfm = 400000; %主频 dt = 1e-4; %时间域采样间隔 number = 2000; %采样点数 t = linspace(-0.05,0.05,number)' ...
c语言雷克子波函数表达式,雷克子波频率研究
引用本文张莹莹, 王云专, 石颖, 等. 2017. 雷克子波频率研究. 地球物理学进展, 2017, 32(5): 2162-2167, doi: 10.6038/pg20170542 ZHANG ...
用matlab作地震波vsp图,利用VSP提高叠后地面地震资料分辨率
为提高地面地震的分辨率,特别是地面地震中深层的分辨率,人们通常采用反Q滤波.时变谱白化.振幅补偿.DMO(倾角时差).小波高分辨处理等方法[1],而这些方法均依赖于精确的反褶积算子.子波反褶积是地震资 ...
c语言地震子波6,地震子波反褶积程序
<地震子波反褶积程序>由会员分享,可在线阅读,更多相关<地震子波反褶积程序(19页珍藏版)>请在人人文库网上搜索. 1.盲信号实验报告盲解卷积算法姓名: 丁宪成系别: 电 ...
一维地震子波合成记录c语言,地震子波波形显示及一维地震合成记录
<地震子波波形显示及一维地震合成记录>由会员分享,可在线阅读,更多相关<地震子波波形显示及一维地震合成记录(22页珍藏版)>请在人人文库网上搜索. 1.地震勘探原理实验一地震子 ...

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

震源子波matlab,VSP震源子波应用相关推荐

最新文章

热门文章