一种基于模型统计拾取的深水OBC资料多次波压制方法技术

本发明专利技术涉及一种基于模型统计拾取的深水OBC资料多次波压制方法，包括以下步骤：1)原始地震资料预处理；2)预处理后的地震资料道集分选；3)确定需要压制的多次波同相轴；4)读取一个道集，在该道集上拾取出需要压制的多次波同相轴；5)平滑拾取的多次波同相轴；6)以平滑后的多次波同相轴为中线，对共偏移距道集上的各道做静态时移，将平滑后的多次波同相轴拉平，再沿多次波同相轴分时窗统计时窗内各道多次波的振幅特征，建立多次波模型；7)从拉平的地震资料中减去多次波模型；8)将减去多次波模型后得到的一次波数据沿步骤6)中的多次波同相轴中线反向静态时移；9)以步骤5)得到的平滑后的多次波同相轴为参考，在下一个道集上重复步骤4)～9)。

【技术实现步骤摘要】
一种基于模型统计拾取的深水OBC资料多次波压制方法
本专利技术涉及石油天然气地震勘探领域，特别是关于一种基于模型统计拾取的深水OBC资料多次波压制方法。
技术介绍
多次波是海上地震勘探中最常见的一种规则干扰波，它源于震源激发的地震波在几个强反射界面之间的多次反射。多次波容易导致在接收的地震剖面上，来自地下深部目的层的有效反射信号被多次波强能量所掩盖，进而使叠加速度场求取不准确，在偏移剖面上形成很多假象。为此，行业内先后发展了预测反褶积、二维滤波、SRME(Surface-RelatedMultipleElimination,表面相关多次波去除)等一系列压制多次波的技术方法。这些方法都有各自的基本假设和适用条件，在一定范围内取得了良好的应用效果。然而，在对深水OBC(OceanBottomCable，海底电缆)采集的多波多分量地震资料压制多次波时，上述这些方法却不能得到满意的多次波压制效果。例如预测反褶积压制鸣震，通过设计一个反滤波算子压缩地震子波，对多次波的周期性要求很高，且只能适用于浅水近道的情况。二维滤波根据动校正后的剩余时差分离一次波和多次波，先要拾取均方根速度，再到变换域中滤除多次波，可是OBC资料的放炮点和接收点不在同一个基准面上，当海水较深时，静校正误差很大，如图1所示，从(d)中可以看出OBC资料与拖缆采集资料的时差，既随偏移距变化，又随反射层深度变化，埋深越浅，变化越大，最大时差为海水层的t0/2(t0是海水层的自激自收时间)。所以，通过简单的高程静校正，将OBC资料的接收点校正到海面，会产生很大的误差，且海水越深，误差越大。这会导致均方根速度拾取不准确，从而影响对动校剩余时差的区分，而且无论F-K变换还是Radon变换都会引入假频，降低资料的信噪比。且由于近道的一次波与多次波剩余时差很小，所以该方法对近道效果不好，容易伤害到有效波。此外，对OBC横波分量而言，纵波激发横波接收，其均方根速度拾取本身也是十分困难的。SRME是目前海上拖缆采集资料去多次波效果最好，应用最广的一项技术，它最大的特点是无需任何地下信息，仅利用地震数据自身进行时空褶积来预测多次波。实际应用中，要注意数据规则化和提供近道数据。可是，将SRME用于深水OBC资料同样得不到满意的多次波压制效果。这是因为，SRME将多次波看作是由两段在“地表反射点”相连的一次反射路径组成的(“地表反射点”位置射线传播遵循Snell定律)。而OBC数据是在海面放炮，海底接收，缺少了从接收点到炮点的射线路径。虽然可以通过静校正将海底电缆校正到炮点所在的基准面上，但当海水较深时，静校正误差很大，且不满足地震波传播规律，这使得SRME对OBC资料的中远道多次波压制效果很差，如图2所示，随偏移距增加，预测的多次波误差迅速增大。由此可见，目前生产中还缺乏有效的针对OBC资料的多次波压制技术。将常规多次波压制技术应用到OBC资料的主要困难是：①炮点和接收点不在同一个基准面上；②静校正量时变空变；③射线路径不对称；④转换波速度拾取困难；⑤多次波能量越强压制越困难。现今主流的多次波压制技术包都含多次波预测和多次波相减两个步骤。而其中的多次波预测的准确性常常是方法成功的关键。
技术实现思路
针对上述问题，本专利技术的目的是提供一种能够有效压制深水OBC地震资料中多次波干扰的基于模型统计拾取的深水OBC资料多次波压制方法。为实现上述目的，本专利技术采取以下技术方案：一种基于模型统计拾取的深水OBC资料多次波压制方法，包括以下步骤：1)对原始地震资料进行预处理；2)对预处理后的地震资料进行道集分选，使得多次波在分选后的地震道集上呈现周期性特征或在速度谱上呈现低速特征，其中，对于2-D地震资料，按带符号的共偏移距剖面进行道集分选；对于3-D地震资料，按CDP进行道集分选；3)通过2-D共偏移距道集剖面上反射同相轴的周期性特征，或通过3-D共CDP道集剖面结合速度谱能量团信息识别出低速的多次波强能量，确定地震资料中所有需要压制的多次波同相轴；4)读取其中一个道集，在该道集上拾取出需要压制的多次波同相轴；5)平滑拾取的多次波同相轴，滤除多次波同相轴中偏离的局部极值；6)以平滑后的多次波同相轴为中线，对共偏移距道集上的各道做静态时移，将平滑后的多次波同相轴拉平，再沿多次波同相轴分时窗统计时窗内各道多次波的振幅特征，建立多次波模型；7)从拉平的地震资料中减去步骤6)所建立的多次波模型，得到压制多次波后的一次波数据；8)将步骤7)得到的一次波数据沿步骤6)中的多次波同相轴中线反向静态时移，校正回原来的位置，得到压制多次波后的地震道集；9)以步骤5)得到的平滑后的多次波同相轴为参考，在下一个道集上重复步骤4)～步骤8)，直至完成地震资料全部道集的多次波压制处理，得到压制多次波后的地震资料。所述步骤1)中的预处理包括去除坏炮、坏道，地震道极性反转校正，去除随机噪音，球面扩散补偿，炮集间的能量均衡。所述步骤4)在该道集上拾取出需要压制的多次波同相轴，根据多次波同相轴特点分为以下两种情况：①针对海底一次反射波能量强、同相轴连续的海底多次波，首先，使用计算机自动拾取海底一次反射波；然后，根据海底一次反射波的到达时间T0和当前地震道集的偏移距ofs计算出海底一阶、二阶多次波的达到时间T1和T2，其中，v代表声波在海水中的传播速度；最后，分别沿T1和T2在上下各1/4子波波长范围内根据多次波同相轴的子波极性和振幅相位特征，自动拾取海底一阶、二阶多次波同相轴；②针对局部强能量多次波，首先，在多次波同相轴上按相等的水平间隔人工拾取若干参考点，以及人工拾取多次波同相轴的左右端点，人工拾取参考点时，在其它强能量反射轴斜交所拾取的多次波同相轴的交点附近适当添加参考点；然后，采用样条插值计算参考基线，并根据多次波同相轴的子波极性、振幅相位特征或时窗内相邻道的相关分析信息，沿着参考基线逐点追踪出要拾取的多次波同相轴。所述步骤5)中平滑拾取的多次波同相轴时，采用中值滤波。所述步骤6)中统计时窗内各道多次波的振幅特征采用平均值或中值叠加方法。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术包括多次波模型建立和相减，即首先根据道集上多次波的走时特征识别并解释多次波，再沿同相轴统计多次波振幅特征，采用中值叠加或统计方法建立多次波模型，然后从原始地震资料中减去多次波，从而达到多次波压制的目的，可以处理放炮点位置与接收点位置不在同一个水平面上的情况，避开了二维滤波和SRME等常规方法处理深水OBC资料时需要面对的静校正和速度拾取两个难题，不仅适用于深水OBC地震资料，对常规海上拖缆采集地震资料也同样有效。如图3所示，对Marmousi-II模型OBC纵波资料压制海底多次波处理，该模型海水水深450m，海底水平，成像剖面中的海底多次波得到了有效压制，此时，诸如预测反褶积、二维滤波和SRME等常用方法都不满足假设条件，不能取得很好的多次波压制效果。2、本专利技术采用不同的方式对地震资料进行分选，对于2-D地震资料，采用Soffset(带符号的共偏移距剖面)分选地震数据，对于3-D资料地震资料，采用CDP(共深度点道集)进行道集分选，有助于多次波的识别和追踪。3、本专利技术根据多次波同相轴特点，采用不同方法拾取需要压制本文档来自技高网...

【技术保护点】
一种基于模型统计拾取的深水OBC资料多次波压制方法，包括以下步骤：1）对原始地震资料进行预处理；2）对预处理后的地震资料进行道集分选，使得多次波在分选后的地震道集上呈现周期性特征或在速度谱上呈现低速特征，其中，对于2‑D地震资料，按带符号的共偏移距剖面进行道集分选；对于3‑D地震资料，按CDP进行道集分选；3）通过2‑D共偏移距道集剖面上反射同相轴的周期性特征，或通过3‑D共CDP道集剖面结合速度谱能量团信息识别出低速的多次波强能量，确定地震资料中所有需要压制的多次波同相轴；4）读取其中一个道集，在该道集上拾取出需要压制的多次波同相轴；5）平滑拾取的多次波同相轴，滤除多次波同相轴中偏离的局部极值；6）以平滑后的多次波同相轴为中线，对共偏移距道集上的各道做静态时移，将平滑后的多次波同相轴拉平，再沿多次波同相轴分时窗统计时窗内各道多次波的振幅特征，建立多次波模型；7）从拉平的地震资料中减去步骤6）所建立的多次波模型，得到压制多次波后的一次波数据；8）将步骤7）得到的一次波数据沿步骤6）中的多次波同相轴中线反向静态时移，校正回原来的位置，得到压制多次波后的地震道集；9）以步骤5）得到的平滑后的多次波同相轴为参考，在下一个道集上重复步骤4）～步骤9），直至完成地震资料全部道集的多次波压制处理，得到压制多次波后的地震资料。...

【技术特征摘要】
1.一种基于模型统计拾取的深水OBC资料多次波压制方法，包括以下步骤：1)对原始地震资料进行预处理；2)对预处理后的地震资料进行道集分选，使得多次波在分选后的地震道集上呈现周期性特征或在速度谱上呈现低速特征，其中，对于2-D地震资料，按带符号的共偏移距剖面进行道集分选；对于3-D地震资料，按CDP进行道集分选；3)通过2-D共偏移距道集剖面上反射同相轴的周期性特征，或通过3-D共CDP道集剖面结合速度谱能量团信息识别出低速的多次波强能量，确定地震资料中所有需要压制的多次波同相轴；4)读取其中一个道集，在该道集上拾取出需要压制的多次波同相轴；5)平滑拾取的多次波同相轴，滤除多次波同相轴中偏离的局部极值；6)以平滑后的多次波同相轴为中线，对共偏移距道集上的各道做静态时移，将平滑后的多次波同相轴拉平，再沿多次波同相轴分时窗统计时窗内各道多次波的振幅特征，建立多次波模型；7)从拉平的地震资料中减去步骤6)所建立的多次波模型，得到压制多次波后的一次波数据；8)将步骤7)得到的一次波数据沿步骤6)中的多次波同相轴中线反向静态时移，校正回原来的位置，得到压制多次波后的地震道集；9)以步骤5)得到的平滑后的多次波同相轴为参考，在下一个道集上重复步骤4)～步骤8)，直至完成地震资料全部道集的多次波压制处理，得到压制多次波后的地震资料。2.如权利要求1所述的一种基于模型统计拾取的深水OBC资料多次波压制方法，其特征在于：所述步骤1)中的预处理包括去除坏炮、坏道，地震道极性反转校正，去除随机噪音，球面扩散补偿，炮集间的能量均衡。3.如权利要求1所述的一种基于模型统计拾取的深水OBC资料多次波压制方法，其特征在于：所述步骤4)在该道集上拾取出需要压制的多次波同相轴，根据多次波同相轴特点分为以下两种情况：①针对海底一次反射波能量强、同相轴连续的海底多次波，首先，使用计算机自动拾取海底一次反射波；然后，根据海底一次反射波的到达时间T0和当前地震道集的偏移距ofs计算出海底一阶、二阶多次波的达到时间T1和T2，其中，v代表...

【专利技术属性】
技术研发人员：李维新，薛东川，张云鹏，江南森，黄小刚，杨俊，桑淑云，王艳冬，王小六，王清振，
申请(专利权)人：中国海洋石油总公司，中海油研究总院，
类型：发明
国别省市：北京;11