一种随机噪声衰减的方法及装置制造方法及图纸

本申请实施例提供了一种随机噪声衰减的方法及装置，其中方法包括：获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录；根据所述频域内的单炮记录，构建汉克尔矩阵集；对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解；根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵的前K个部分和，所述预设数值K不小于所述单炮记录中的倾角个数；将每一个所述前K个部分和矩阵沿其反对角线平均，得到平均后的矩阵；用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值；将替换后的频域内的单炮记录进行傅里叶反变换。本申请实施例可以实现高保真度的压制随机噪声。

【技术实现步骤摘要】

本申请涉及石油勘探地震资料处理
，尤其是涉及一种随机噪声衰减的方法及装置。
技术介绍
在地震传播过程中，不仅激发了发射波，也激发了面波等干扰波，这就导致地震检波器既检测到了反射信号，又检测到了噪声。地震勘探中的噪声主要有两种：随机噪声和规则噪声，所谓随机噪声就是没有固定频率，在地震整个记录中随机出现，频带很宽，视速度不确定，无一定传播方向的一种噪声，因而很难利用随机噪声和有效波之间频率上的差异或传播方向上的差异对其进行压制。在强噪声下，即使较强的反射信号也可能不明显，甚至淹没在噪声中，地震勘探难度增大。因此，压制随机噪声，提高地震资料信噪比是当前勘探的研究热点之一。维纳滤波是地震勘探随机噪声压制最早也是最经典的方法，但是它要求已知信号或噪声的相关函数或功率谱，实际上，这很难满足，因此维纳滤波的效果不理想。另一方面在精确地震勘探中，要求地震资料处理满足高保真度、高信噪比和高分辨率的原则，现有技术在压制随机噪声过程中，在提高信噪比的同时或多或少都压制了部分有效信号，这就导致无法实现高保真度的要求。
技术实现思路
本申请实施例的目的在于提供一种随机噪声衰减的方法及装置，可以实现高保真度的压制随机噪声。为达到上述目的，本申请实施例提供了一种随机噪声衰减的方法，该方法包括：获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录；根据所述频域内的单炮记录，构建汉克尔矩阵集，所述汉克尔矩阵集中的每一个汉克尔矩阵的元素为时间频率变换值，所述时间频率变换值为所述单炮记录在预设频率范围内每一个频率切片所对应的时间频率变换值；对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解；根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵的前K个部分和，获取所述每一个汉克尔矩阵的前K个部分和矩阵，所述预设数值K不小于所述单炮记录中的倾角个数；将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵；将所述平均后的矩阵进行傅里叶反变化，获取傅里叶反变换后的矩阵；用所述傅里叶反变换后的矩阵的值代替所述单炮记录中的对应数值。将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵；用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值；本申请实施例还提供了一种随机噪声衰减的装置，用以实现高保真度的压制随机噪声，该装置包括：获取模块，用于获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录；构建模块，用于根据频域内的单炮记录，构建汉克尔矩阵集，所述汉克尔矩阵集中的每一个汉克尔矩阵的元素为时间频率变换值，所述时间频率变换值为所述单炮记录在预设频率范围内每一个频率切片所对应的时间频率变换值；分解模块，用于对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解；计算模块，用于根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵前K个部分和，获取每一个汉克尔矩阵的前K个部分和矩阵，所述预设数值K不小于所述单炮记录中的倾角个数；平均模块，用于将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵；替换模块，用于用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值；反变换模块，用于将替换后的频域内的单炮记录进行傅里叶反变换。由以上本实施例提供的技术方案可见，本实施例在频域范围内进行了压制随机噪声处理，首先通过将频域内的单炮记录中各个检波器在每一频率切片对应的值构建成汉克尔矩阵集，然后对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解，取前K个部分和代替原来的汉克尔矩阵，然后将汉克尔矩阵还原到时间域，所述K值等于所述单炮记录中的倾角个数。由于地震数据的有效秩数不超过数据的倾角个数，取前K个部分和矩阵代替原有的矩阵，可以保留下所有的地震信号，且前K个部分和都是奇异值较大的部分，奇异值越大，对地震信号贡献越大，因此取前K个部分和可以实现高保真的压制随机噪声。附图说明此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的一部分，并不构成对本申请实施例的限定。在附图中：图1为本申请实施例的一种随机噪声衰减的方法流程示意图；图2为本申请实施例的无噪声的剖面图；图3为本申请实施例的K＝1时得到的还原剖面；图4为本申请实施例的K＝2时得到的还原剖面图5为本申请实施例的K＝3时得到的还原剖面图6为本申请实施例的K＝1时得到的还原剖面与无噪声的剖面图想减结果；图7为本申请实施例的K＝2时得到的还原剖面与无噪声的剖面图想减结果；图8为本申请实施例的K＝3时得到的还原剖面与无噪声的剖面图想减结果；图9为本申请实施例的含随机噪声的剖面图；图10为本申请实施例的含随机噪声的剖面处理后的效果示意图；图11为本申请实施例的一种随机噪声衰减的模块示意图。具体实施方式为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于解释本申请实施例，但并不作为对本申请实施例的限定。下面结合附图，对本申请实施例的具体实施方式作进一步的详细说明。参考图1所示，本申请实施例一种随机噪声衰减的方法可以包括以下步骤：S101，获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录。S102，根据所述频域内的单炮记录，构建汉克尔矩阵集。所述汉克尔矩阵集中的每一个汉克尔矩阵的元素为所述单炮记录在预设频率范围内每一个频率切片的时间频率变换值，即预设频率范围内每一个频率切片在各个地震道的时间频率变换值。所述各个地震道与各个检波器相对应。具体的，各个地震道的数据对应各个检波器的数据。每一个频率切片的频率数值根据采样定理，结合采样时间间隔决定。汉克尔矩阵为每一条副对角线上元素都相等的方阵。这里构建汉克尔矩阵的目的就是为了估算地震数据的有效秩数。S103，对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解。S104，根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵的前K个部分和，获取每一个汉克尔矩阵的前K个部分和矩阵，又称为降秩近似和，所述预设数值K不小于所述单炮记录中的倾角个数。对于时间—频域的地震数据，国外学者发现数据的有效秩数不超过数据的倾角个数，因此本申请的实施例首先通过汉克尔矩阵矩阵估算地震数据的有效秩数，再计算汉克尔矩阵的前K个部分和矩阵。按照有效秩数不超过数据的倾角个数的原则，理论上，当K等于数据倾角个数时，前K个部分和矩阵中就会包含全部的地震数据。S105，将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵。由于噪声的影响，使得前K个部分和矩阵可能不再满足汉克尔矩阵的形式要求，为了恢复其汉克尔矩阵的形式，因此这里对前K个部分和矩阵进行反对角线平均处理。S106，用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值。平均后的矩阵与汉克尔矩阵集中的汉克尔矩阵是一一对应的，且平均后的矩阵也具备汉克尔矩阵的形式要求，即每一条副对角线上元素都相等。用平均后的矩阵的元素对应替换单炮记录在预设频率范围内的每一频率切片中各个检波器的频率值，得到替换后的频域内的单炮记录。S107，将替换后的频域内的单炮记录进行傅里叶反变换。由图1所示的流程图可知，本实本文档来自技高网...

【技术保护点】
一种随机噪声衰减的方法，其特征在于，所述包括以下步骤：获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录；根据所述频域内的单炮记录，构建汉克尔矩阵集，所述汉克尔矩阵集中的每一个汉克尔矩阵的元素为时间频率变换值，所述时间频率变换值为所述单炮记录在预设频率范围内每一个频率切片所对应的时间频率变换值；对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解；根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵的前K个部分和，获取所述每一个汉克尔矩阵的前K个部分和矩阵，所述预设数值K不小于所述单炮记录中的倾角个数；将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵；用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值；将替换后的频域内的单炮记录进行傅里叶反变换。

【技术特征摘要】
1.一种随机噪声衰减的方法，其特征在于，所述包括以下步骤：获取单炮记录，并对其进行傅里叶变换，得到频域内的单炮记录；根据所述频域内的单炮记录，构建汉克尔矩阵集，所述汉克尔矩阵集中的每一个汉克尔矩阵的元素为时间频率变换值，所述时间频率变换值为所述单炮记录在预设频率范围内每一个频率切片所对应的时间频率变换值；对所述汉克尔矩阵集中的每一个汉克尔矩阵进行奇异值分解；根据所述奇异值分解的结果，并结合预设数值K，计算所述每一个汉克尔矩阵的前K个部分和，获取所述每一个汉克尔矩阵的前K个部分和矩阵，所述预设数值K不小于所述单炮记录中的倾角个数；将每一个所述前K个部分和矩阵沿其反对角线平均，得到对应的平均后的矩阵；用每一个所述平均后的矩阵的元素值替换频域内的单炮记录的对应数值；将替换后的频域内的单炮记录进行傅里叶反变换。2.如权利要求1所述的方法，其特征在于，所述将所述每一个汉克尔矩阵的前K个部分和矩阵沿反对角线平均包括：用所述前K个部分和矩阵中的每一条反对角线元素的算术平均值代替相应的反对角线上的所有元素。3.如权利要求1所述的方法，其特征在于，所述计算所述每一个汉克尔矩阵的前K个部分和为：FK＝I1+I2+…IK式中，FK表示所述前K个部分和，Ii表示第i个加权特征图像。4.如权利要求1所述的方法，其特征在于，所述预设频率范围为所述单炮记录中的地震数据频率范围。5.如权利要求1所述的方法，其特征在于，所述预设数值K大于所述单炮记录中的倾角个数。6.一种随机噪声衰减的装置，其特征在于，所...

【专利技术属性】
技术研发人员：王万里，杨午阳，魏新建，禄娟，何欣，陈德武，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11