一种转换波角道集抽取方法
【专利摘要】本发明是提高地震处理中计算转换波入射角的速度和精度的转换波角道集抽取方法，利用纵波入射和横波反射转换波数据形成振幅随偏移距变化的正常时差校正后的道集、转换波速度及其各向异性参数，计算纵波速度和纵波各向异性参数、转换波炮检距对应的纵波炮检距和零炮检距的转换波旅行时、纵波旅行时和层间纵波速度，由纵波时差方程和层间纵波速度计算转换波入射角θ，计算每一样点的入射角，将此样点放置在入射角的角度道上，对所有样点重复此过程，得到转换波角道集。本发明计算容易并可提高抽取转换波角道集的速度和精度。
【专利说明】一种转换波角道集抽取方法
【技术领域】
[0001] 本发明涉及石油地球物理勘探领域，具体是一种用于转换波AVA分析、转换波振幅分析和地震属性反演中的转换波角道集抽取方法。
【背景技术】
[0002]在多波多分量地震勘探中，有效利用转换波信息能够获得更加可靠的岩石弹性参数。从转换波AVA分析中能够提取有用的岩性信息，用来预测地层岩性和裂缝等介质信息。在转换波AVA方程中，自变量是转换波入射角。转换波AVA分析中的一个基本步骤就是进行将炮检距域道集转换为角度域道集，即进行转换波入射角的计算。目前，有两种方法可以计算转换波入射角:一是直射线法；二是射线追踪法。
[0003]前者在已知转换波的转换点位置和反射层深度的条件下，用直射线法求取转换波入射角；后者在已知层速度和层厚度的情况下，用射线追踪法求取转换波入射角。
[0004]直射线法的计算速度快，但是其计算的角度总是偏小，误差较大；射线追踪法的计算结果准确，但是其运算量大，计算速度慢。
【发明内容】
[0005]本发明目的在于提供一种提高计算转换波入射角的速度和精度的转换波角道集抽取方法。
[0006]本发明采用以下步骤实现:
[0007]I)在野外利用纵波震源激发地震波并利用三分量检波器记录地震波，按照多波多分量地震资料处理流程对采集的纵波入射和横波反射(PS)转换波数据进行相对振幅保持的高保真处理，形成振幅随偏移距变化(AVO)的正常时差校正(NMO)后的道集；
[0008]2)在横波反射(PS)转换波数据处理过程中通过速度分析得到横波反射(PS)转换波时差校正(NMO)速度Vc2、垂直速度比Y C1、有效速度比Yeff和PS转换波各向异性参数
V.A eff ?
[0009]3)由横波反射(PS)转换波时差校正(NMO)速度Vc2、垂直速度比Y。、有效速度比
Yeff和PS转换波各向异性参数X eff，计算纵波速度Vp2和纵波各向异性参数Il ；
[0010]4)由横波反射(PS)转换波时间tc、横波反射(PS)转换波炮检距X、横波反射(PS)转换波时差校正(NMO)速度Vc2、垂直速度比Y。和有效速度比Yeff，计算与横波反射(PS)转换波炮检距X对应的纵波炮检距X’ ；
[0011]所述的横波反射(PS)转换波炮检距X与纵波炮检距X’具有相同的震源O和相同的转换点水平位置Xe ;
[0012]5)由横波反射(PS)转换波时间t。、横波反射(PS)转换波炮检距X、横波反射(PS)转换波时差校正(NMO)速度Vc2、垂直速度比Ytl、有效速度比Yeff和PS转换波各向异性参数X eff,计算零炮检距的转换波旅行时；
[0013]6)由零炮检距的转换波旅行时tC(l和垂直速度比Y C1计算零炮检距的纵波旅行时lPO ,
[0014]7)由零炮检距的纵波旅行时tP(l、纵波速度Vp2、纵波各向异性参数η和纵波炮检距X’得到纵波时差方程tp ；
[0015]8)由零炮检距的纵波旅行时tP(l和纵波速度Vp2得到层间纵波速度Vint ；
[0016]9)由纵波时差方程tP和层间纵波速度Vint,用射线参数法计算PS转换波入射角Θ ;
[0017]10)计算横波反射(PS)转换波数据中每一样点的入射角，将此样点放置在入射角的角度道上，对所有样点重复此过程，得到横波反射(PS)转换波角道集。
[0018]本发明容易实现，实验结果证明提高了抽取转换波角道集的速度和精度。
【专利附图】
【附图说明】
[0019]图1是转换波入射角计算示意图，图中纵波反射:P波入射，P波反射；PS转换波反射:p波入射，S波反射；
[0020]图2是转换波炮检距道集，图中纵轴为时间，单位秒，横轴为道号，道间距为100米，偏移距为O米；
[0021]图3是不同方法计算的入射角结果对比，根据表1中的转换波速度数据，分别用本发明、射线追踪方法和直射线法计算的转换波炮检距为1000米的不同深度界面的入射角的结果；
[0022]图4是转换波角道集，图中纵轴为时间，单位秒，横轴为角度道号，角度增量5度，第一道为O度。
【具体实施方式】
[0023]本发明是通过转化将转换波炮检距、旅行时和速度映射为相应的纵波炮检距、旅行时和速度，用纵波的射线参数法来求取转换波的入射角。如图1所示，炮检距为X的转换波入射角为Θ等于炮检距为X’的纵波入射角Θ，通过计算得到炮检距为X’的纵波入射角Θ，即得到了炮检距为X的转换波入射角Θ。
[0024]本发明具体实施步骤如下:
[0025]I)在野外利用纵波震源激发地震波并利用三分量检波器记录地震波，按照多波多分量地震资料处理流程对采集的纵波入射和横波反射(PS)转换波数据进行相对振幅保持的高保真处理，形成振幅随偏移距变化(AVO)的正常时差校正(NMO)后的道集(图2)；
[0026]2)在PS转换波数据处理过程中，通过对PS转换波数道集进行速度分析得到PS转换波速度文件(表1)，此文件包括零炮检距的PS转换波时间tra、PS转换波时差校正(NMO)速度Vc2、垂直速度比Ytl、有效速度比Yeff和转换波各向异性参数xeff;
[0027]表1
[0028]
tco(s) Vc2 (m/s)Y οY eff I X eff I
0.6749 2481.9492 1.8414 1.8414 0.0000
【权利要求】
1.一种转换波角道集抽取方法，特点是采用以下步骤实现: 1)在野外利用纵波震源激发地震波并利用三分量检波器记录地震波，按照多波多分量地震资料处理流程对采集的纵波入射和横波反射转换波数据进行相对振幅保持的高保真处理，形成振幅随偏移距变化的正常时差校正后的道集； 2)在横波反射转换波数据处理过程中通过速度分析得到横波反射转换波时差校正速度Vc2、垂直速度比Y C1、有效速度比Yeff和横波反射转换波各向异性参数Xeff; 3)由横波反射转换波时差校正速度Vc2、垂直速度比Y。、有效速度比Yeff和横波反射转换波各向异性参数X rff，计算纵波速度Vp2和纵波各向异性参数η ； 4)由横波反射转换波时间t。、横波反射转换波炮检距X、横波反射转换波时差校正速度Vc2、垂直速度比Ytl和有效速度比Yeff,计算与横波反射转换波炮检距X对应的纵波炮检距X ； 所述的横波反射转换波炮检距X与纵波炮检距X’具有相同的震源O和相同的转换点水平位置Xc; 5)由横波反射转换波时间t。、横波反射转换波炮检距X、横波反射转换波时差校正速度V?、垂直速度比Y O、有效速度比Yeff和横波反射转换波各向异性参数Xeff，计算零炮检距的转换波旅行时te(l; 6)由零炮检距的横波反射转换波旅行时tC(l和垂直速度比Ytl计算零炮检距的纵波旅行时tP(l ； 7)由零炮检距的纵波旅行时tro、纵波速度Vp2、纵波各向异性参数η和纵波炮检距X’得到纵波时差方程tP ； 8) 由零炮检距的纵波旅行时tP(l和纵波速度Vp2得到层间纵波速度Vint； 9)由纵波时差方程tP和层间纵波速度Vint，用射线参数法计算横波反射转换波入射角θ ; 10)计算横波反射转换波数据中每一样点的入射角，将此样点放置在入射角的角度道上，对所有样点重复此过程，得到横波反射转换波角道集。
【文档编号】G01V1/30GK103576197SQ201210272650
【公开日】2014年2月12日 申请日期:2012年8月2日 优先权日:2012年8月2日
【发明者】张铁强, 钱忠平, 孙鹏远, 马光凯 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司