CN112379427B - 面波噪声压制方法及装置 - Google Patents
面波噪声压制方法及装置 Download PDFInfo
- Publication number
- CN112379427B CN112379427B CN202011197554.7A CN202011197554A CN112379427B CN 112379427 B CN112379427 B CN 112379427B CN 202011197554 A CN202011197554 A CN 202011197554A CN 112379427 B CN112379427 B CN 112379427B
- Authority
- CN
- China
- Prior art keywords
- surface wave
- seismic data
- frequency
- prediction result
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 58
- 239000006185 dispersion Substances 0.000 claims abstract description 48
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 230000010363 phase shift Effects 0.000 claims abstract description 23
- 230000001629 suppression Effects 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims description 38
- 238000001914 filtration Methods 0.000 claims description 21
- 238000012163 sequencing technique Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 230000002742 anti-folding effect Effects 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 12
- 238000005070 sampling Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000009466 transformation Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种面波噪声压制方法及装置，其中该方法包括：设定面波时窗，确定待处理数据的频散谱；利用频散谱得到频率速度曲线；对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据。本发明可以解决面波噪声压制中的假频能量问题。
Description
技术领域
本发明涉及地球物理地震勘探技术领域，尤其涉及一种面波噪声压制方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
随着勘探程度的不断深入，复杂地区和复杂油气藏的油气勘探，已经成为地球物理勘探的主要目标。复杂地区和复杂油气藏的油气勘探是一项非常复杂的系统工程，其中提高信噪比是复杂地区地震勘探数据处理重要环节之一，而压制噪声则在地震资料处理中格外重要。在地震资料处理中，面波是一种极为常见的噪音，它表现为低频、低速、强能量的特征，并具有频散性和多模性。为了使反射波和折射波等有效信号更加突显出来，需要有效的压制消除面波噪声。面波去除的好坏直接影响后续的处理，如何有效的去除面波，一直是地震资料处理中的重点问题。
由于面波与有效波在速度、频率、能量以及波长等方面存在明显区别，传统方法主要是对比面波和有效波在速度、频率、能量以及波长其中一个或者几个特征上的区别，来压制面波的干扰，目前比较成熟的做法有：F-K域滤波，区域滤波，区域噪声衰减以及S变换、小波变换等，这些方法在实际应用中都取得了一定的效果，但对于由于空间稀疏采样导致的假频能量问题一直不能得到很好的解决。
发明内容
本发明实施例提供一种面波噪声压制方法，用以解决面波噪声压制中的假频能量问题，该方法包括：
获取采集的地震数据，对地震数据进行预处理；
设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱；
自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线；
对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；
将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；
利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据；
其中，设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱，包括：
设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱；
其中，根据面波的空间分布设定面波时窗，以使利用面波时窗圈定范围中包含面波；
对地震数据进行预处理，包括：
对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序；
将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果后，还包括：
对所选时窗内的每一道都进行如上操作，得到时窗内面波的最佳预测结果。
本发明实施例还提供一种面波噪声压制装置，用以解决面波噪声压制中的假频能量问题，该装置包括：
获取模块，用于获取采集的地震数据，对地震数据进行预处理；
确定模块，用于设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱；
曲线拟合模块，用于自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线；
相移处理模块，用于对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；
时移处理模块，用于将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；
确定模块，还用于利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据；
其中，确定模块具体用于：
设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱；
其中，根据面波的空间分布设定面波时窗，以使利用面波时窗圈定范围中包含面波；
获取模块具体用于：
对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序；
时移处理模块具体用于：
对所选时窗内的每一道都进行如上操作，得到时窗内面波的最佳预测结果。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述面波噪声压制方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述面波噪声压制方法。
本发明实施例中，利用面波的频散特征，通过分析地震数据得到频散谱，自动拾取准确的频率-速度曲线，不用对面波的多模型逐一预测，对于面波的预测和压制更加灵活有效，对假频面波同样具有有效性；在自动拾取频率-速度曲线后，对面波数据进行相移处理及时移处理，修正面波的频散，进而预测面波，根据面波的预测结果来获取压制面波后的地震数据。本发明实施例不需要在数学变换域内对面波数据和有效信号进行分离，从而避免了在空间采样不足时面波信号与有效信号重叠严重、无法分开的问题，所以对空间采样不足导致的假频能量问题同样具有有效性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中一种面波噪声压制方法的流程图；
图2为本发明实施例中一去噪前的原始炮集；
图3为本发明实施例中经过预处理后的原始炮集以及设定的面波时窗范围；
图4为本发明实施例中通过面波时窗选取的待处理数据；
图5为本发明实施例中对图4所示待处理数据经过非线性信号相似性分析得到的频散谱；
图6为本发明实施例中某一道及其相邻道数据经过频散校正后的结果；
图7为本发明实施例中面波时窗范围内的面波预测结果；
图8为本发明实施例中消除面波后的地震数据；
图9为本发明实施例中另一去噪前炮集；
图10为本发明实施例中另一炮集数据的面波预测结果；
图11为本发明实施例中压制面波后的另一炮集的地震数据；
图12为本发明实施例中一种面波噪声压制装置的结构示意图；
图13为本发明实施例中一种计算机设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
本发明实施例提供了一种面波噪声压制方法，使得去噪后能较好的保持原始地震数据的波形特征，且能在很大程度上提高地震数据的信噪比和分辨率，为地震成像、属性提取和油藏开发等提供有利条件。如图1所示，该方法包括步骤101至步骤106：
步骤101、获取采集的地震数据，对地震数据进行预处理。
其中，预处理的过程包括：对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序。
其中，解编、置观测系统、静校正、反褶积和能量补偿处理均是地震数据处理中的常用技术手段，对于其具体实现方法，在此不做赘述。
步骤102、设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱。
具体的，设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱。
其中，根据面波的空间分布设定面波时窗，以使利用面波时窗圈定范围中包含面波。面波时窗和滤波频率范围均由用户设定。
准确获取待处理数据的频散特征，拾取准确的频散谱，对校正面波频散效应、预测面波至关重要。传统方法得到的频散谱在低频部分不如高频部分能量聚焦，在面波存在的低频段难以获取可靠的时移量，从而严重影响面波的预测。本发明实施例使用的非线性信号相似性分析方法得到频散谱在高低频能量聚焦上都很平衡，而且分辨率更高。非线性信号相似性分析方法使用的公式如下所示：
其中，SNLSC(ω,Vph；σ)为规则化后的非线性信号相似度表征；Sπ为背景值；SNL(ω,Vph)为两个接收点处的地震数据对某一频率ω和数据频率在ω的相速度Vph的非线性信号相关；σ为用于调整频散谱分辨率的系数；
其中，
T为信号长度；t为时间；x为信号坐标。
步骤103、自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线。
步骤104、对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果。
本发明实施例中，利用相移校正面波的频散效应。具体的，利用如步骤102中公式计算得到的频散谱SNLSC(ω,Vph；σ)，无需对面波的多模逐一处理，经过自动拾取得到不同频率与其对应的速度值，利用这一频率-速度关系，再加上不同道数据的偏移距信息对相邻道做相移来校正面波的频散效应，这一过程使用的公式为：
其中，Us(ω)为重建后的面波数据；U(ω)为原始数据；Xwin为预测道面波时窗的最大范围；i为虚部。
消除频散特征的相邻道数据中，面波信号展现出真线性的特征，而非原始炮集中的视线性，而有效信号在炮集中的双曲特征变成了随机噪声，通过叠加可以有效消除有效信号的影响，使面波信号得到增强，从而得到预测出来的面波噪声。
在一种实现方式中，在利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理之前，还可以对频率-速度曲线进行平滑处理。
步骤105、将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果。
步骤106、利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据。
本发明对二维或三维的地震数据均可以进行处理，下面将结合具体的附图来介绍上述实现过程。
①、对采集的地震数据进行预处理，如图2所示，将数据按照炮线号、炮点号、检波线号、检波点号的索引排列。
②、设定面波时窗，对炮集上的每条排列的地震数据按照用户设定的滤波频率范围进行滤波处理，得到如图3所示的数据，图中虚线框为根据用户参数选取的面波时窗。
③、对图4中面波时窗内的数据进行非线性信号相似性分析，得到如图5所示的高精度频散谱。参照图5，可知得到的频散谱的分辨率较高。
④、本发明实施例中，不需要人工拾取，也不需要对频散谱中的多阶进行分阶处理，通过自动拾取从该频散谱中得到频率-速度曲线。
⑤、根据该频率-速度曲线，在对时窗内任一道进行预测时，利用这一频散关系加上各道的偏移距信息，对其相邻道进行相移处理。图6所示即为对该数据时窗内第143道进行预测时，对其左右各142道进行了相移处理得到的结果。该处理校正了面波的频散，原始的面波特征标明为进行的线性，通过护理面波的同相轴成为线性，此时有效信号表现为随机噪声。
⑥、对相移后的相邻道进行叠加，有效信号的能量互相抵消，而波能量得到增强，得到用户所选时窗内一道的面波预测结果。对该结果与原始数据进行相位匹配后通过一定的时移就得到该道面波的最佳预测结果。对所选时窗内的所有道都进行如上操作，得到时窗内面波的最佳预测结果，如图7所示。
⑦、从去噪前的原始炮集中消减掉步骤⑥得到的最佳预测结果就可以得到压制面波噪声的地震数据，如图8所示。
图9所示为另一去噪前炮集，该数据道间距为50m，由于空间采样不足，面波在炮集上表现出严重的假频特征。图10所示为通过本发明实施例中方法得到的面波预测结果。图11为消除面波噪声后的结果，可以看出假频面波得到了很好的压制，有效信号得到凸显，信噪比明显得到提升。
本发明实施例中，利用面波的频散特征，通过分析地震数据得到频散谱，自动拾取准确的频率-速度曲线，不用对面波的多模型逐一预测，对于面波的预测和压制更加灵活有效，对假频面波同样具有有效性；在自动拾取频率-速度曲线后，对面波数据进行相移处理及时移处理，修正面波的频散，进而预测面波，根据面波的预测结果来获取压制面波后的地震数据。本发明实施例不需要在数学变换域内对面波数据和有效信号进行分离，从而避免了在空间采样不足时面波信号与有效信号重叠严重、无法分开的问题，所以对空间采样不足导致的假频能量问题同样具有有效性。
本发明实施例中还提供了一种面波噪声压制装置，如下面的实施例所述。由于该装置解决问题的原理与面波噪声压制方法相似，因此该装置的实施可以参见面波噪声压制方法的实施，重复之处不再赘述。
如图12所示，该装置1200包括获取模块1201、确定模块1202、曲线拟合模块1203、相移处理模块1204和时移处理模块1205。
其中，获取模块1201，用于获取采集的地震数据，对地震数据进行预处理；
确定模块1202，用于设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱；
曲线拟合模块1203，用于自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线；
相移处理模块1204，用于对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；
时移处理模块1205，用于将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；
确定模块1202，还用于利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据。
在本发明实施例的一种实现方式中，获取模块1201，用于：
对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序。
在本发明实施例的一种实现方式中，确定模块1202，用于：
设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；
对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱。
在本发明实施例的一种实现方式中，装置1200还包括：
平滑处理模块1206，用于对频率-速度曲线进行平滑处理。
本发明实施例中，利用面波的频散特征，通过分析地震数据得到频散谱，自动拾取准确的频率-速度曲线，不用对面波的多模型逐一预测，对于面波的预测和压制更加灵活有效，对假频面波同样具有有效性；在自动拾取频率-速度曲线后，对面波数据进行相移处理及时移处理，修正面波的频散，进而预测面波，根据面波的预测结果来获取压制面波后的地震数据。本发明实施例不需要在数学变换域内对面波数据和有效信号进行分离，从而避免了在空间采样不足时面波信号与有效信号重叠严重、无法分开的问题，所以对空间采样不足导致的假频能量问题同样具有有效性。
本发明实施例还提供一种计算机设备，图13为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的面波噪声压制方法中全部步骤，该计算机设备具体包括如下内容：
处理器(processor)1301、存储器(memory)1302、通信接口(CommunicationsInterface)1303和通信总线1304；
其中，所述处理器1301、存储器1302、通信接口1303通过所述通信总1304完成相互间的通信；所述通信接口1303用于实现相关设备之间的信息传输；
所述处理器1301用于调用所述存储器1302中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的面波噪声压制方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述面波噪声压制方法的计算机程序。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种面波噪声压制方法，其特征在于，所述方法包括：
获取采集的地震数据，对地震数据进行预处理；
设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱；
自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线；
对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；
将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；
利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据；
其中，设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱，包括：
设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱；
其中，根据面波的空间分布设定面波时窗，以使利用面波时窗圈定范围中包含面波；
对地震数据进行预处理，包括：
对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序；
将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果后，还包括：
对所选时窗内的每一道都进行如上操作，得到时窗内面波的最佳预测结果。
2.根据权利要求1所述的方法，其特征在于，在利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理之前，所述方法还包括：
对频率-速度曲线进行平滑处理。
3.一种面波噪声压制装置，其特征在于，所述装置包括：
获取模块，用于获取采集的地震数据，对地震数据进行预处理；
确定模块，用于设定面波时窗，圈定预处理后的地震数据中的待处理数据，确定待处理数据的频散谱；
曲线拟合模块，用于自动拾取频散谱中各频率值对应的速度最大值，利用各频率值和对应的速度最大值对应的点进行曲线拟合，得到频率速度曲线；
相移处理模块，用于对待处理数据中的每一道的指定道数的相邻道，利用频率速度曲线，以及相邻道各道的偏移距对相邻道各道分别做相移处理，并分别叠加相移处理后的各道得到各道的面波预测结果；
时移处理模块，用于将相邻道每一道的面波预测结果与预处理后的地震数据中相应道做相位匹配，得到匹配最佳的时移量，并利用匹配最佳的时移量对面波预测结果中相应道做时移处理，得到该道面波的最佳预测结果；
确定模块，还用于利用每一道的预处理后的地震数据减去该道面波的最佳预测结果，得到压制面波噪声的地震数据；
其中，确定模块具体用于：
设定面波时窗和滤波频率范围，在一条检波线上的预处理后的地震数据中，利用面波时窗圈定出待处理数据，利用滤波频率范围对待处理数据进行滤波；对滤波后的待处理数据进行非线性信号相似性分析，得到面波的频散谱；
其中，根据面波的空间分布设定面波时窗，以使利用面波时窗圈定范围中包含面波；
获取模块具体用于：
对地震数据进行解编、置观测系统、静校正、反褶积和/或能量补偿处理，并按照炮线号、炮点号、检波线号、检波点号对地震数据进行重排，其中，先按照炮线号的大小顺序对地震数据进行排序；对于相同炮线号的地震数据，按照炮点号的大小顺序对该些地震数据进行排序；对于相同炮点号的地震数据，按照检波线号的大小顺序对该些地震数据进行排序；对于相同检波线号的地震数据，按照检波点号的大小顺序对该些地震数据进行排序；
时移处理模块具体用于：
对所选时窗内的每一道都进行如上操作，得到时窗内面波的最佳预测结果。
4.根据权利要求3所述的装置，其特征在于，所述装置还包括：
平滑处理模块，用于对频率-速度曲线进行平滑处理。
5.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至2任一所述方法。
6.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至2任一所述方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197554.7A CN112379427B (zh) | 2020-10-30 | 2020-10-30 | 面波噪声压制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197554.7A CN112379427B (zh) | 2020-10-30 | 2020-10-30 | 面波噪声压制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112379427A CN112379427A (zh) | 2021-02-19 |
| CN112379427B true CN112379427B (zh) | 2024-03-26 |
Family
ID=74576501
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011197554.7A Active CN112379427B (zh) | 2020-10-30 | 2020-10-30 | 面波噪声压制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379427B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103215A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 三维高密度地震勘探记录叠前面波压制方法 |
| CN103217709A (zh) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | 一种提高地震数据信噪比和分辨率的面波衰减方法 |
| CN104570115A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油天然气集团公司 | 一种面波衰减方法及装置 |
| CN104614769A (zh) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | 一种压制地震面波的聚束滤波方法 |
| CN106707342A (zh) * | 2016-12-12 | 2017-05-24 | 中国石油天然气集团公司 | 共炮点道集多级面波压制方法和装置 |
| WO2020033465A1 (en) * | 2018-08-10 | 2020-02-13 | University Of Houston System | Surface wave estimation and removal from seismic data |
| CN111103621A (zh) * | 2019-12-09 | 2020-05-05 | 东华理工大学 | 一种主动源共成像点叠加多道面波分析方法 |
| CN114428333A (zh) * | 2020-09-24 | 2022-05-03 | 中国石油化工股份有限公司 | 一种基于面波重建的散射面波压制方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7917295B2 (en) * | 2008-04-30 | 2011-03-29 | Westerngeco L.L.C. | Modeling and filtering coherent noise in seismic surveying |
-
2020
- 2020-10-30 CN CN202011197554.7A patent/CN112379427B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103215A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 三维高密度地震勘探记录叠前面波压制方法 |
| CN103217709A (zh) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | 一种提高地震数据信噪比和分辨率的面波衰减方法 |
| CN104570115A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油天然气集团公司 | 一种面波衰减方法及装置 |
| CN104614769A (zh) * | 2015-02-05 | 2015-05-13 | 中铁隧道集团有限公司 | 一种压制地震面波的聚束滤波方法 |
| CN106707342A (zh) * | 2016-12-12 | 2017-05-24 | 中国石油天然气集团公司 | 共炮点道集多级面波压制方法和装置 |
| WO2020033465A1 (en) * | 2018-08-10 | 2020-02-13 | University Of Houston System | Surface wave estimation and removal from seismic data |
| CN111103621A (zh) * | 2019-12-09 | 2020-05-05 | 东华理工大学 | 一种主动源共成像点叠加多道面波分析方法 |
| CN114428333A (zh) * | 2020-09-24 | 2022-05-03 | 中国石油化工股份有限公司 | 一种基于面波重建的散射面波压制方法及系统 |
Non-Patent Citations (4)
| Title |
|---|
| David Li等.Removal of dispersive borehole waves and beam imaging of fractures around the wellbore using sonic logging data .SEG International Exposition and 88th annual Meeting.2018,全文. * |
| Yingcai Zheng and Hao Hu.Nonlinear Signal Comparison and High-Resolution Measurement of Surface-Wave Dispersion.Bulletin of the Seismological Society of America.2017,第1552-1555页. * |
| 伍敦仕 ; 孙成禹 ; 林美言 ; .基于互相关相移的主动源地震面波频散成像方法.地球物理学进展.2017,(第04期),全文. * |
| 董烈乾 ; 张慕刚 ; 张翊孟 ; .预测自适应相减面波压制方法.石油地球物理勘探.2016,(第06期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112379427A (zh) | 2021-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102590862B (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
| CN111190227B (zh) | 基于残差卷积生成对抗模型的低信噪比地震数据消噪方法 | |
| CN109581516B (zh) | 曲波域统计量自适应阈值探地雷达数据去噪方法及系统 | |
| CN102681014A (zh) | 基于多项式拟合的规则线性干扰压制方法 | |
| CN111679328A (zh) | 一种基于变分模态分解的瞬变电磁探测信号的降噪方法 | |
| CN104345341A (zh) | 一种基于区域约束的分频段能量地震面波处理方法 | |
| CN113885077A (zh) | 一种基于深度学习的多震源地震数据分离方法 | |
| CN113962244A (zh) | 瑞雷波地震数据噪声去除方法、存储介质及电子设备 | |
| CN110261910A (zh) | 基于自适应稀疏s变换的地震数据面波去除方法 | |
| CN104133248A (zh) | 一种高保真声波干扰压制方法 | |
| CN111736222A (zh) | 单炮数据信噪比确定方法及装置 | |
| CN113219539B (zh) | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 | |
| Yuan et al. | Application of ICEEMDAN to noise reduction of near-seafloor geomagnetic field survey data | |
| CN112379427B (zh) | 面波噪声压制方法及装置 | |
| Zhang et al. | Segmented noise reduction based on Brillouin-spectrum-partition in Brillouin optical time domain sensors | |
| CN112198547A (zh) | 深层或超深层地震资料处理方法及装置 | |
| CN105093275B (zh) | 基于速度模型的有效信号提取iss波场分离方法 | |
| CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
| CN112213784A (zh) | 复杂地表地震数据一次处理快速静校正方法 | |
| CN114428333A (zh) | 一种基于面波重建的散射面波压制方法及系统 | |
| CN113589384B (zh) | 基于信号随偏移距变化特征的叠前道集保幅去噪方法 | |
| CN113050156B (zh) | 沙漠区地震采集检波点优选方法及装置 | |
| CN114740530B (zh) | 基于双曲时窗约束的中高频拟线性噪声压制方法及装置 | |
| CN113917539B (zh) | 火山岩覆盖区地震数据叠前道集处理方法、系统和装置 | |
| CN113721297B (zh) | 利用速度特征压制沙丘鸣震的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |