地震波信号的处理方法、装置、设备及存储介质制造方法及图纸

本发明专利技术提供一种地震波信号的处理方法、装置、设备及存储介质，属于地震勘探技术领域。方法包括：获取多次波图像；对动校正前的CMP道集和多次波图像进行二值化处理；对二值化图像进行连通域分析；对动校正前的CMP道集以及动校正后的CMP道集进行多次波衰减。本发明专利技术可以实现对动校正前以及动校正后的CMP道集的多次波进行衰减处理，继而可以降低多次波对地震勘探的影响，有效地提高地震勘探的准确性。

Processing method, device, equipment and storage medium of seismic wave signal

【技术实现步骤摘要】
地震波信号的处理方法、装置、设备及存储介质
本专利技术涉及地震勘探
，具体而言，涉及一种地震波信号的处理方法、装置、设备及存储介质。
技术介绍
在地震勘探中，多次地震波是一种能量很强的干扰信号，会影响地震成像，使反射波形态发生畸变，从而降低地震资料的成像效果，甚至导致对地震资料的解释和地震介质构造的错误认识。为了提高地震勘探的准确性，降低多次波对地震勘探的影响，通常会对多次地震波进行衰减。目前一般通过拉东变换(RadonTransformation)，对道集的多次波信号进行处理，以实现对多次地震波信号的衰减处理。然而，多次波通常在不同炮检距上具有不同时差，炮检距越小，多次波和一次波的时差越小，拉东变换衰减多次波的能力越弱，因此近炮检距地震道中的多次波同相轴难以消除。同时，远炮检距地震道往往存在道集未拉平和拉伸效应，也影响了多次波衰减效果，使得地震波信号成像效果较差。
技术实现思路
本专利技术的目的在于提供一种地震波信号的处理方法、装置、设备及存储介质，可以提高对多次波的衰减效果，提高基于衰减后的地震波的成像效果。本专利技术的实施例是这样实现的：本专利技术实施例的一方面，提供一种地震波信号的处理方法，包括：对动校正后的共中心点CMP道集进行处理，获取多次波图像；分别对动校正前的CMP道集和多次波图像进行二值化处理，得到动校正前的CMP道集的第一二值化图像，和多次波图像的第二二值化图像；对第一二值化图像、第二二值化图像进行连通域分析，确定第二二值化图像中多次波连通域位置；根据第二二值化图像中多次波连通域位置，确定第一二值化图像中多次波连通域位置；根据第一二值化图像中多次波连通域位置，对动校正前、后的CMP道集进行多次波衰减。可选地，对动校正后的共中心点CMP道集进行处理，获取多次波图像，包括：采用高分辨率抛物拉东变换对动校正后的CMP道集进行处理，得到多次波图像。可选地，分别对动校正前的CMP道集和多次波图像进行二值化处理，得到动校正前的CMP道集的第一二值化图像，和多次波图像的第二二值化图像，包括：将动校正前的CMP道集中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将动校正前的CMP道集中幅度值小于预设阈值的像素点的灰度值修改为最小灰度值，得到第一二值化图像；将多次波图像中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将多次波图像中幅度值小于预设阈值的像素点的灰度值修改为最小灰度值，得到第二二值化图像。可选地，对第一二值化图像、第二二值化图像进行连通域分析，确定第二二值化图像中多次波连通域位置，包括：分别对第一二值化图像和第二二值化图像进行连通域分析，确定第一二值化图像的多个连通域和第二二值化图像的多个连通域；确定第一二值化图像中的多个连通域中，第一个一次波的采样时间；根据采样时间，确定第二二值化图像中的多次波连通域。可选地，分别对第一二值化图像和第二二值化图像进行连通域分析，确定第一二值化图像的多个连通域和第二二值化图像的多个连通域包括：对第一二值化图像进行连通域分析，确定第一二值化图像的所有连通域以及每个连通域的面积，并将第一二值化图像的所有连通域中面积小于或等于预设面积的连通域进行过滤，得到第一二值化图像的多个连通域；对第二二值化图像进行连通域分析，确定第二二值化图像的所有连通域以及每个连通域的面积，并将第二二值化图像的所有连通域中面积小于或等于预设面积的连通域进行过滤，得到第二二值化图像的多个连通域。可选地，对第一二值化图像进行连通域分析，确定第一二值化图像的所有连通域以及每个连通域的面积，包括：对第一二值化图像进行连通域分析，确定第一二值化图像的所有连通域；对第一二值化图像的每个连通域中的像素点进行标记，并根据第一二值化图像的每个连通域中标记的像素点的总面积，确定第一二值化图像的每个连通域的面积。相应的，对第二二值化图像进行连通域分析，确定第二二值化图像的所有连通域以及每个连通域的面积，包括：对第二二值化图像进行连通域分析，确定第二二值化图像的所有连通域；对第二二值化图像的每个连通域中的像素点进行标记，并根据第二二值化图像的每个连通域中标记的像素点的总面积，确定第二二值化图像的每个连通域的面积。可选地，该方法还包括：根据第一二值化图像的每个连通域中像素点的标记个数，确定第一二值化图像的每个连通域的像素个数，并根据第一二值化图像的每个连通域中标记的像素点的位置，确定第一二值化图像的每个连通域的位置；根据第二二值化图像的每个连通域中像素点的标记个数，确定第二二值化图像的每个连通域的像素个数，并根据第二二值化图像的每个连通域中标记的像素点的位置，确定第二二值化图像的每个连通域的位置。本专利技术实施例的另一方面，提供一种地震波信号的处理装置，包括：第一处理模块，用于对动校正后的共中心点CMP道集进行处理，获取多次波图像。第二处理模块，用于分别对动校正前的CMP道集和多次波图像进行二值化处理，得到动校正前的CMP道集的第一二值化图像，和多次波图像的第二二值化图像。连通域分析模块，用于对第一二值化图像、第二二值化图像进行连通域分析，确定第二二值化图像中多次波连通域位置，还用于根据第二二值化图像中多次波连通域位置，确定第一二值化图像中多次波连通域位置。衰减模块，用于根据第一二值化图像中多次波连通域位置，对动校正前、后的CMP道集进行多次波衰减。可选地，第一处理模块，具体用于采用高分辨率抛物拉东变换对动校正后的CMP道集进行处理，得到多次波图像。可选地，第二处理模块，具体用于将动校正前的CMP道集中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将动校正前的CMP道集中幅度值小于预设阈值的像素点的灰度值修改为最小灰度值，得到第一二值化图像。可选地，第二处理模块，具体用于将多次波图像中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将多次波图像中幅度值小于预设阈值的像素点的灰度值修改为最小灰度值，得到第二二值化图像。可选地，连通域分析模块，具体用于分别对第一二值化图像和第二二值化图像进行连通域分析，确定第一二值化图像的多个连通域和第二二值化图像的多个连通域。可选地，连通域分析模块，具体用于确定第一二值化图像中的多个连通域中，第一个一次波的采样时间。可选地，连通域分析模块，具体用于根据采样时间，确定第二二值化图像中的多次波连通域位置。可选地，连通域分析模块，具体用于对第一二值化图像进行连通域分析，确定第一二值化图像的所有连通域以及每个连通域的面积，并将第一二值化图像的所有连通域中面积小于或等于预设面积的连通域进行过滤，得到第一二值化图像的多个连通域。可选地，连通域分析模块，具体用于对第二二值化图像进行连通域分析，确定第二二值化图像的所有连通域以及每个连通域的面积，本文档来自技高网...

【技术保护点】
1.一种地震波信号的处理方法，其特征在于，包括：/n对动校正后的共中心点CMP道集进行处理，获取多次波图像；/n分别对动校正前的CMP道集和所述多次波图像进行二值化处理，得到所述动校正前的CMP道集的第一二值化图像，和所述多次波图像的第二二值化图像；/n对所述第一二值化图像、所述第二二值化图像进行连通域分析，确定所述第二二值化图像中多次波连通域位置；/n根据所述第二二值化图像中多次波连通域位置，确定所述第一二值化图像中多次波连通域位置；/n根据所述第一二值化图像中多次波连通域位置，对所述动校正前、后的CMP道集进行多次波衰减。/n

【技术特征摘要】
1.一种地震波信号的处理方法，其特征在于，包括：
对动校正后的共中心点CMP道集进行处理，获取多次波图像；
分别对动校正前的CMP道集和所述多次波图像进行二值化处理，得到所述动校正前的CMP道集的第一二值化图像，和所述多次波图像的第二二值化图像；
对所述第一二值化图像、所述第二二值化图像进行连通域分析，确定所述第二二值化图像中多次波连通域位置；
根据所述第二二值化图像中多次波连通域位置，确定所述第一二值化图像中多次波连通域位置；
根据所述第一二值化图像中多次波连通域位置，对所述动校正前、后的CMP道集进行多次波衰减。


2.如权利要求1所述的方法，其特征在于，所述对动校正后的共中心点CMP道集进行处理，获取多次波图像，包括：
采用高分辨率抛物拉东变换对所述动校正后的CMP道集进行处理，得到所述多次波图像。


3.如权利要求1所述的方法，其特征在于，所述分别对动校正前的CMP道集和所述多次波图像进行二值化处理，得到所述动校正前的CMP道集的第一二值化图像，和所述多次波图像的第二二值化图像，包括：
将所述动校正前的CMP道集中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将所述动校正前的CMP道集中幅度值小于所述预设阈值的像素点的灰度值修改为最小灰度值，得到所述第一二值化图像；
将所述多次波图像中幅度值大于或等于预设阈值的像素点的灰度值修改为最大灰度值，并将所述多次波图像中幅度值小于所述预设阈值的像素点的灰度值修改为最小灰度值，得到所述第二二值化图像。


4.如权利要求1所述的方法，其特征在于，所述对所述第一二值化图像、所述第二二值化图像进行连通域分析，确定所述第二二值化图像中多次波连通域位置，包括：
分别对所述第一二值化图像和所述第二二值化图像进行连通域分析，确定所述第一二值化图像的多个连通域和所述第二二值化图像的多个连通域；
确定所述第一二值化图像中的多个连通域中，第一个一次波的采样时间；
根据所述采样时间，确定所述第二二值化图像中的多次波连通域。


5.如权利要求4所述的方法，其特征在于，所述分别对所述第一二值化图像和所述第二二值化图像进行连通域分析，确定所述第一二值化图像的多个连通域和所述第二二值化图像的多个连通域包括：
对所述第一二值化图像进行连通域分析，确定所述第一二值化图像的所有连通域以及每个连通域的面积，并将所述第一二值化图像的所有连通域中面积小于或等于预设面积的连通域进行过滤，得到所述第一二值化图像的多个连通域；
对所述第二二值化图像进行连通域分析，确定所述第二二值化图像的所有连通域以及每个连通域的面积，并将所述第二二值化图像的所有连通域中面...

【专利技术属性】
技术研发人员：毛伟建，宋欢，唐欢欢，
申请(专利权)人：中国科学院测量与地球物理研究所，
类型：发明
国别省市：湖北;42