共炮点道集多级面波压制方法和装置制造方法及图纸

本发明专利技术提供了一种共炮点道集多级面波压制方法和装置，其中，该方法包括：在共炮点道集上确定出面波区域和不含面波的反射波区域；确定面波区域中面波的起始频率、截至频率和最低视速度；在面波的起始频率至截至频率范围内对面波区域的数据，以多道统计值为期望振幅谱，进行零相位滤波处理，得到初步压制面波的结果；利用面波的截至频率、最低视速度、空间采样间隔和炮检距计算线性相移算子；根据线性相移算子对初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果。本发明专利技术解决了面波压制方法中所存在的处理精度低、计算效率不高以及面波假频难以消除的技术问题，达到了方便有效地实现共炮点道集多级面波压制的目的。

Method and device for suppressing common shot point set multi-stage surface wave

The invention provides a common shot gathers pressing method and apparatus, multistage surface wave, the method includes determining a surface wave reflection region and free surface wave wave region in the shot gathers; determine the surface wave in the region as the starting frequency, surface wave frequency and the lowest apparent velocity; in the initial frequency of the surface wave as to the frequency range of surface wave data of the regional statistics is expected to pass, amplitude spectrum, zero phase filtering, get surface wave suppression preliminary results; using surface wave as the lowest frequency, apparent velocity, spatial sampling interval and offset calculation of linear phase shift operator; according to the linear phase shift operator on the initial surface wave suppression results of linear phase shift, phase filtering and anti linear phase shift processing, get the results of surface wave suppression. The invention solves the technical problems existing in the processing precision of surface wave suppression in the low calculation efficiency and high frequency surface wave is difficult to eliminate the false, to realize common shot gathers to multistage surface wave suppression.

【技术实现步骤摘要】
共炮点道集多级面波压制方法和装置
本专利技术涉及地质勘探
，特别涉及一种共炮点道集多级面波压制方法和装置。
技术介绍
面波是陆上勘探中常见的干扰波，主要分布在近炮检距，呈“扫帚状”，具有低频、低速、强振幅和频散等特点。地震记录中面波与反射波的同相轴交织在一起，严重影响着地震记录的信噪比，给叠前振幅信息的有效利用带来极大的干扰，从而影响了地震资料处理的质量。合理地压制面波，有利于改善地震资料的品质，提高地震数据的信噪比，对于后续的叠前属性分析、波形反演以及改善深部地质目标的成像质量都是至关重要的。目前，主要有以下几种面波压制方法：1)基于传统滤波的面波压制方法在数据采集阶段可以通过检波器的组合来压制面波，但这种方法使得反射波的接收频带变窄，而且不能对已经接收到的数据进行处理，缺乏灵活性。随着地震数据处理技术的发展，出现了多种滤波方法来衰减面波干扰，主要包括基于频域滤波类的高通滤波、带通滤波以及基于视速度特征差异的FK滤波等。在频域中反射波与面波的频谱不是完全分开的，部分区域存在重叠。因而这些方法都存在缺陷，例如，高通滤波在去除低频面波时也压制了反射波的低频成分，而且截频滤波会引起波形的失真，造成视觉上的假象。为了克服这些缺点，许多滤波方法利用反射波与面波在几何特性上的差异进行面波压制。例如：Randon变换将时空域的数据变换到τ-p域中，以直线形式分布的面波在τ-p域中是以一系列能量点的形式存在，切除这些能量点，反变换到时空域就能够达到压制面波的目的。但由于“端点效应”的影响，使得面波在τ-p域中不能完全集中，从而导致面波压制不彻底。2)基于多尺度分析的面波压制方法地震信号是非平稳信号，各种多尺度分析工具(例如：S变换，Wavelet、Ridgelet、Curvelet变换等)也广泛应用于地震数据的面波衰减中。这类方法基本上都能取得一定的处理效果，但又存在各自的不足。例如：在曲波域中可以综合利用面波与反射波在空间、尺度、方向上的差异来压制面波，但随着尺度的降低，波原子的局域性和方向性都降低，对于面波中低频低波数成分，曲波低尺度波原子难以提供信噪分离所需的精度，给曲波系数的处理带来困难；Askari等基于S变换对FK滤波法进行扩展，通过设计随时间或炮检距变化的视速度滤波器来压制面波，对于面波受到地表非均质体散射或者频散严重的地震记录，能够取得有效的滤波效果，但该方法运算量较大，非常耗时。3)基于稀疏表征策略的面波压制方法随着稀疏表征技术的发展出现了多种地震信号建模的策略，Karsli等利用线性或非线性的扫频信号对面波进行参数化建模，通过反褶积的思路来预测并消除面波干扰。Trad等将单炮记录描述为具有线性同相轴的面波和具有双曲同向轴的反射波的叠加，利用线性Radon变换和双曲Radon变换，结合稀疏约束可以实现面波与反射波波场的分离。此外，利用地震记录中面波与体波信号波形结构的差异性，利用平稳小波变换和局部离散余弦变换作稀疏表示字典，将面波与体波分离问题转化为稀疏优化问题来求解。由于近地表结构变化多样，该类方法的处理精度依赖于所建模型特征是否与实际情况相匹配。由上述描述可知，现有的面波压制方法主要存在以下问题：1)由于近地表特征复杂，依据视速度、频率、线性特征等难以完全地区分面波和反射波；2)多尺度几何分析工具具有局域性表征能力，但随着频率的降低，局域性越来越差，难以精细地刻画低频面波与反射波的特征，给变换系数的处理带来困难，导致处理精度不高，并且计算效率低。3)在现有三维地震勘探观测系统采样间隔条件下，当面波视速度较低时，往往存在假频，现有方法不能有效的压制面波假频能量。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本专利技术实施例提供了一种共炮点道集多级面波压制方法，以解决现有的面波压制方法中所存在的处理精度低、计算效率不高以及面波假频难以消除的技术问题，该方法包括：在共炮点道集上确定出面波区域和不含面波的反射波区域；确定所述面波区域中面波的起始频率、截至频率和最低视速度；在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值；在面波的起始频率至截至频率范围内对所述面波区域的数据，以所述多道统计值为期望振幅谱，进行零相位滤波处理，得到初步压制面波的结果；利用面波的截至频率、最低视速度、空间采样间隔和炮检距计算线性相移算子；利用所述线性相移算子对所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果。在一个实施方式中，在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值，包括：将所述反射波区域的数据变换至频率域，得到频率空间域的数据体；在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值和标准差；将所述几何平均值和所述标准差相加，得到面波的起始频率至截至频率范围内的每个频率上振幅的多道统计值。在一个实施方式中，在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值，包括：对所述频率空间域数据体取自然对数；在面波的起始频率至截至频率范围内的每个频率上，对取自然对数后的结果计算几何平均值；对所述几何平均值取做e指数处理，得到在面波的起始频率至截至频率范围内的所述频率空间域数据体的每个频率上振幅的多道统计几何平均值。在一个实施方式中，利用面波的截至频率、最小视速度、空间采样间隔和炮检距计算线性相移算子，包括：利用面波的截至频率和空间采样间隔计算满足无面波假频的最小视速度；确定所述满足无面波假频的最小视速度与面波的最小视速度之间的大小关系；如果所述满足无面波假频的最小视速度小于面波的最小视速度，则确定线性相移系数为零；如果所述满足无面波假频的最小视速度大于面波的最小视速度，则将所述满足无面波假频的最小视速度的倒数与面波最小视速度的倒数之间的差值作为线性相移系数；将所述线性相移系数、频率和炮检距三者相乘得到线性相移算子。在一个实施方式中，利用所述线性相移算子将所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果，包括：计算所述初步压制面波结果的相位谱；对所述初步压制面波结果的相位谱加上所述线性相移算子得到线性相移后的相位谱；对所述线性相移后的相位谱，取e的复指数运算得到包含相位信息的信号；对所述包含相位信息的信号进行FKK域的滤波处理，消除相位信息中的面波信息，并做去e的复指数处理，得到相位滤波后的相位谱；对所述相位滤波后的相位谱减去线性相移因子，得到反线性相移后的相位谱；对所述反线性相移后的相位谱与所述初步压制面波后的振幅谱结合反变换到时间域得到压制面波的结果。在一个实施方式中，将所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果之后，所述方法还包括：确定得到的所述压制面波后的结果中面波是否被完全压制；如果没有，则重复执行以下操作，直至面波被完全压制：将所述压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到再次面波压制的结果。本专利技术实施例提供了一种共炮点道集多级面波压制装置，以解决现有的面波压制方法中所存在的处理精度低、计算效率不高以及面波假频难以消除的技术问题，该装置包括本文档来自技高网...

【技术保护点】
一种共炮点道集多级面波压制方法，其特征在于，包括：在共炮点道集上确定出面波区域和不含面波的反射波区域；确定所述面波区域中面波的起始频率、截至频率和最低视速度；在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值；在面波的起始频率至截至频率范围内对所述面波区域的数据，以所述多道统计值为期望振幅谱，进行零相位滤波处理，得到初步压制面波的结果；利用面波的截至频率、最低视速度、空间采样间隔和炮检距计算线性相移算子；利用所述线性相移算子对所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果。

【技术特征摘要】
1.一种共炮点道集多级面波压制方法，其特征在于，包括：在共炮点道集上确定出面波区域和不含面波的反射波区域；确定所述面波区域中面波的起始频率、截至频率和最低视速度；在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值；在面波的起始频率至截至频率范围内对所述面波区域的数据，以所述多道统计值为期望振幅谱，进行零相位滤波处理，得到初步压制面波的结果；利用面波的截至频率、最低视速度、空间采样间隔和炮检距计算线性相移算子；利用所述线性相移算子对所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果。2.如权利要求1所述的方法，其特征在于，在面波的起始频率至截至频率范围内计算所述反射波区域中的每个频率上振幅的多道统计值，包括：将所述反射波区域的数据变换至频率域，得到频率空间域的数据体；在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值和标准差；将所述几何平均值和所述标准差相加，得到面波的起始频率至截至频率范围内的每个频率上振幅的多道统计值。3.如权利要求2所述的方法，其特征在于，在面波的起始频率至截至频率范围内计算所述频率空间域数据体的每个频率上振幅的多道统计几何平均值，包括：对所述频率空间域数据体取自然对数；在面波的起始频率至截至频率范围内的每个频率上，对取自然对数后的结果计算几何平均值；对所述几何平均值取做e指数处理，得到在面波的起始频率至截至频率范围内的所述频率空间域数据体的每个频率上振幅的多道统计几何平均值。4.如权利要求1所述的方法，其特征在于，利用面波的截至频率、最小视速度、空间采样间隔和炮检距计算线性相移算子，包括：利用面波的截至频率和空间采样间隔计算满足无面波假频的最小视速度；确定所述满足无面波假频的最小视速度与面波的最小视速度之间的大小关系；如果所述满足无面波假频的最小视速度小于面波的最小视速度，则确定线性相移系数为零；如果所述满足无面波假频的最小视速度大于面波的最小视速度，则将所述满足无面波假频的最小视速度的倒数与面波最小视速度的倒数之间的差值作为线性相移系数；将所述线性相移系数、频率和炮检距三者相乘得到线性相移算子。5.如权利要求1所述的方法，其特征在于，利用所述线性相移算子将所述初步压制面波的结果进行线性相移、相位滤波和反线性相移处理，得到压制面波的结果，包括：计算所述初步压制面波结果的相位谱；对所述初步压制面波结果的相位谱加上所述线性相移算子得到线性相移后的相位谱；对所述线性相移后的相位谱，取e的复指数运算得到包含相位信息的信号；对所述包含相位信息的信号进行FKK域的滤波处理，消除相位信息中的面波信息，并做去e的复指数处理，得到相位滤波后的相位谱；对所述相位滤波后的相位谱减去线性相移因子，得到反线性...

【专利技术属性】
技术研发人员：王德营，凌云，王永军，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11