专利名称：一种高精度地震波层速度反演方法
技术领域：
本发明涉及地震勘探数据处理技术，本发明目的是提供一种利用地震数据的相位谱属性，且具有很强抗噪性的高精度的地震波层速度反演方法。
背景技术：
在地震勘探中，零偏垂直地震剖面(VSP)由于炮点离井口很近，可认为不同深度接收到的下行直达波具有相同的传播路径，能直接利用不同深度地震记录的初至时间反演地震波的层速度，为地面地震处理以及解释提供有力手段。而目前的VSP资料处理中，都是先对地震记录在时间域进行重采样，提高采样率，再进行初至拾取。其具体做法是利用该道地震记录中相邻几个样点的值内插出更多的样点值，比如将I毫秒采样率的信号内插成
O.0001毫秒采样率的信号，然后再进行初至拾取，利用获得的初至，计算相邻道之间的时·差，再用相邻道的深度差除以时间差得到层速度。由于这种提高采样率的方法仅仅只是一种数学算法，其内插出的信号并不是真实信号，不能提高原始数据的分辨率，同时内插结果受内插方法和资料信噪比的影响很大，无法真正起稿初至拾取的精度，导致反演的层速度不够精确，制约了 VSP资料的应用价值。

发明内容
本发明目的是提供一种利用地震数据的相位谱属性，且具有很强抗噪性的高精度地震波层速度反演方法。本发明具体步骤包括I)地面震源激发，井下检波器接收得到VSP数据；2)对地震记录进行傅里叶变换得到振幅谱，将高频部分的振幅谱清零，然后反傅里叶变换，得到只保留低频信号的地震记录，然后拾取每道地震记录的初至；3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱，然后在深度方向再进行一次傅里叶变换得到F-K谱，并将上行波对应的F-K谱清零，然后分别在波数方向和频率方向进行反傅里叶变换得到下行波；4)从每道的下行波中沿初至往后开一时窗，对该时窗内的信号进行傅里叶变换，得到每个频率的相位谱；5)重复步骤4)，得到所有道下行波每个频率的相位谱；6)求取相邻道中相同频率的相位谱之差；7)对步骤6)中的相位谱之差进行最小二乘拟合，得到相邻道初至差；8)相邻道的深度差除以相邻道初至差得到高精度的层速度。步骤6)所述的相邻道相位谱之差计算如下式Si+1(< ) =「Si^t — &ti+x)e-jl7^ftdt =「河、卜dt
J-OOJ-OO则有θ +1(ω) = 2JifAti+1
式中，Si (t_ Λ ti+1)为深度为Iii的检波器得到的下行波为Si⑴，深度为hi+1的检波器得到的下行波；Si+1( )为si+1(t)的傅里叶变换；θ +1(ω)为Si+1(t)的傅里叶变换的相位谱Δ ti+1为第i+Ι道和第i道初至差；j为复数单位。本发明具有很强抗随即干扰能力，对激发子波和底层的吸收衰减情况适应能力强，不仅算法简单还能大大节约工作量，且反演的层速度精度很高。


图I为本发明截取的时窗内信号对应的频谱；图2为相邻道中相同频率的相位谱之差；图3为反演层速度对比。
具体实施例方式以下结合附图详细说明本发明。本发明是一种利用地震数据的相位谱属性，且具有很强抗噪性的高精度层速度反演方法。本发明具体实现步骤如下I)地面震源激发，井下检波器接收得到VSP数据；2)对地震记录进行傅里叶变换得到振幅谱，将高频部分的振幅谱清零，然后反傅里叶变换，得到只保留低频信号的地震记录，然后拾取每道地震记录的初至；3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱，然后在深度方向再进行一次傅里叶变换得到F-K谱，并将上行波对应的F-K谱清零，然后分别在波数方向和频率方向进行反傅里叶变换得到下行波；4)从每道的下行波中沿初至往后开一时窗，对该时窗内的信号进行傅里叶变换，得到每个频率的相位谱；图I截取的该时窗内信号对应的频谱；5)重复步骤4)，得到所有道下行波每个频率的相位谱；6)求取相邻道中相同频率的相位谱之差；图2为相邻道中相同频率的相位谱之差;设深度为Iii的检波器得到的下行波为Si (t)，深度为hi+1的检波器得到的下行波为Si(t_Ati+1)，根据傅里叶变换的性质信号时间域的时移对应其傅里叶变换频率域的相移；有
权利要求
1.一种高精度地震波层速度反演方法，特点是采用以下步骤 1)地面震源激发，井下检波器接收得到VSP数据； 2)对地震记录进行傅里叶变换得到振幅谱，将高频部分的振幅谱清零，然后反傅里叶变换，得到只保留低频信号的地震记录，然后拾取每道地震记录的初至； 3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱，然后在深度方向再进行一次傅里叶变换得到F-K谱，并将上行波对应的F-K谱清零，然后分别在波数方向和频率方向进行反傅里叶变换得到下行波； 4)从每道的下行波中沿初至往后开一时窗，对该时窗内的信号进行傅里叶变换，得到每个频率的相位谱； 5)重复步骤4)，得到所有道下行波每个频率的相位谱； 6)求取相邻道中相同频率的相位谱之差； 7)对步骤6)中的相位谱之差进行最小二乘拟合，得到相邻道初至差； 8)相邻道的深度差除以相邻道初至差得到高精度的层速度。
2.根据权利要求I所述的方法，特点是步骤6)所述的相邻道相位谱之差计算如下式
全文摘要
本发明是勘探数据处理的高精度的地震波层速度反演方法，井下接收数据变换得到振幅谱，将高频振幅谱清零再反傅里叶变换，保留低频信号记录，拾取每道初至，在深度方向进行傅里叶变换得到F-K谱，将上行波对应的F-K谱清零，在波数方向和频率方向反傅里叶变换得到下行波，从每道的下行波中沿初至往后开一时窗，傅里叶变换得到每个频率的相位谱，求取相邻道中相同频率的相位谱之差，进行最小二乘拟合得到相邻道初至差，相邻道的深度差除以初至差得到层速度。本发明具有很强抗随即干扰能力，对激发子波和底层的吸收衰减情况适应能力强，不仅算法简单反演精度高。
文档编号G01V1/28GK102890289SQ201110203589
公开日2013年1月23日 申请日期2011年7月20日 优先权日2011年7月20日
发明者张固澜 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司