一种地震数据的噪声去除方法及装置制造方法及图纸

本发明专利技术提供了一种地震数据的噪声去除方法及装置，其中，方法包括：根据获取的待处理地震数据，通过C‑C算法计算得到相空间重构参数：嵌入维度和延迟时间；根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构；计算相空间每个相点的邻域以及每个领域内所有邻域点的均值；基于所述领域点的均值进行局部投影，得到所述待处理地震数据的噪声分量。本发明专利技术通过在时间域内对地震信号进行处理，无需考虑噪声类型，具有运算方便和处理速度较快的有益效果。

【技术实现步骤摘要】
一种地震数据的噪声去除方法及装置
本专利技术涉及地震勘探
，具体而言，涉及一种地震数据的噪声去除方法及装置。
技术介绍
随着技术形势的发展，勘探目标逐渐向深层和地表复杂地区转移，给地震资料处理工作提出了更高的要求，并且地震数据处理一直以高信噪比、高分辨率、高保真为目标，因此更为精确的计算地震数据中的噪声成为地震数据处理的关键。现有技术中，传统的地震数据噪声去除方法计算复杂程度较高，对地震数据中噪声的计算准确度不高且计算速度较慢。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种地震数据的噪声去除方法和装置，以提高计算的准确度以及运算处理速度。第一方面，本专利技术实施例提供了一种地震数据的噪声去除方法，包括：根据获取的待处理地震数据，通过C-C算法计算得到相空间重构参数：嵌入维度和延迟时间；根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构；计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值；基于所述邻域点的均值进行局部投影，得到所述待处理地震数据的噪声分量。结合第一方面，本专利技术实施例提供了第一方面的第一种可能的实施方式，其中，所述方法，还包括：基于所述噪声分量和所述待处理地震数据计算得到修正值，对所述待处理地震数据进行处理，以去除噪声。结合第一方面，本专利技术实施例提供了第一方面的第二种可能的实施方式，其中，所述根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构，具体包括：利用所述嵌入维度m和所述延迟时间τ，将所述待处理地震数据一维的时间序列重构到多维的相空间中，得到一个由N个相点组成的一个m维相空间。结合第一方面，本专利技术实施例提供了第一方面的第三种可能的实施方式，其中，所述计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值，具体包括：计算所述待处理地震数据的方差，得到初始邻域半径值r；根据所述的嵌入维度m进行迭代，判断是否满足迭代终止条件n≥m，其中n表示邻域内的相点数，m表示邻域内相点最小值；如果满足条件，则终止迭代计算，得到该次迭代对应的邻域半径值r；根据得到的所述邻域半径值，计算得到每个相点的邻域Ui＝{Xj；||Xj-Xi||<r}，Xj表示第j个相点，Xi表示第i个相点；根据所述每个相点的邻域计算得到所有邻域点的均值结合第一方面，本专利技术实施例提供了第一方面的第四种可能的实施方式，其中，所述基于所述邻域点的均值进行局部投影，计算得到所述待处理地震数据的噪声分量，具体包括：根据预先设定的矩阵和由该矩阵确定的对角权重矩阵R确定协方差矩阵C；计算所述协方差矩阵C的特征向量aq以及特征值λq；将所述相空间重构的各个相点投影到由所述协方差矩阵的特征向量组成的相空间中，得到噪声分量θ。第二方面，本专利技术实施例还提供了一种地震数据的噪声去除装置，包括：第一计算模块，用于根据获取的待处理地震数据通过C-C算法计算得到相空间重构参数：嵌入维数和延迟时间；相空间重构模块，用于根据所述嵌入维度和所述延迟时间对所述待处理地震数据进行相空间重构；第二计算模块，用于计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值；投影模块，用于基于所述邻域点的均值进行局部投影，计算得到所述待处理地震数据的噪声分量。结合第二方面，本专利技术实施例提供了第二方面的第一种可能的实施方式，其中，所述装置，还包括：第三计算模块，用于基于所述噪声分量和所述待处理地震数据计算得到修正值，对所述待处理地震数据进行处理，以去除噪声。结合第二方面，本专利技术实施例提供了第二方面的第二种可能的实施方式，其中，所述相空间重构模块，包括：利用所述嵌入维度m和所述延迟时间τ，将所述待处理地震数据一维的时间序列重构到多维的相空间中，得到一个由N个相点组成的一个m维相空间。结合第二方面，本专利技术实施例提供了第二方面的第三种可能的实施方式，其中，所述第二计算模块，包括：第一计算单元，用于计算所述待处理地震数据的方差，得到初始邻域半径值r；判断单元，根据所述的嵌入维度m进行迭代，判断是否满足迭代终止条件n≥m，其中n表示邻域内的相点数，m表示邻域内相点最小值，如果满足条件，则终止迭代计算，得到该次迭代对应的邻域半径值r；第二计算单元，根据所述邻域半径值，计算得到每个相点的邻域Ui＝{Xj；||Xj-Xi||<r}，Xj表示第j个相点，Xi表示第i个相点；根据所述每个相点的邻域计算得到所有邻域点的均值。结合第二方面，本专利技术实施例提供了第二方面的第四种可能的实施方式，其中，所述投影模块，包括：设定单元，用于根据预先设定的矩阵和由该矩阵确定的对角权重矩阵R确定协方差矩阵C；计算单元，用于计算所述协方差矩阵的特征向量aq以及特征值λq；投影单元，用于将所述相空间重构的各个相点投影到由所述协方差矩阵的特征向量组成的相空间中，得到噪声分量θ。本专利技术实施例提供的一种地震数据的噪声去除方法及装置，采用相空间重构地震数据，对噪声分量进行计算，与现有技术中相比，其只需在时间域内对地震信号进行处理，运算方便且处理速度较快。为使本专利技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本专利技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。图1示出了本专利技术一个实施例所提供的一种地震数据的噪声去除方法的流程示意图；图2示出了本专利技术的一个实施例所提供的一种相空间重构参数计算流程示意图；图3示出了本专利技术一个实施例所提供的一种计算每个相点邻域以及邻域点均值的计算流程示意图；图4示出了本专利技术一个实施例所提供的一种对邻域点的均值进行局部投影得到噪声分量的计算流程示意图；图5a示出了本专利技术一个实施例所提供的合成地震数据示意图；图5b示出了本专利技术一个实施例所提供的经过地震数据的噪声去除方法去除噪声后的地震数据示意图；图6a示出了本专利技术一个实施例所提供的合成地震数据的相空间重构结果示意图；图6b示出了本专利技术一个实施例所提供的经过地震数据的噪声去除方法去除噪声后的地震数据相空间重构结果示意图；图7示出了本专利技术一个实施例所提供的一种地震数据的噪声去除装置结构示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本专利技术实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本专利技术的实施例的详细描述并非旨在限制要求保护的本专利技术的范围，而是仅仅表示本专利技术的选定实施例。基于本专利技术的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。考虑到随着技术水平的提高，勘探目标逐渐向深层和地表复杂地区转移，给地震资料处理工作提出了更高的要求，传统地震数据中噪声去除方法复杂程度高，处理速度慢。基于此，本专利技术实施例提供了一种地震数据的噪声去除方法及装置，下面通过实施例进行描述。参照图1所示的实施本文档来自技高网...

【技术保护点】
1.一种地震数据的噪声去除方法，其特征在于，所述方法包括：根据获取的待处理地震数据，通过C‑C算法计算得到相空间重构参数：嵌入维度和延迟时间；根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构；计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值；基于所述邻域点的均值进行局部投影，得到所述待处理地震数据的噪声分量。

【技术特征摘要】
1.一种地震数据的噪声去除方法，其特征在于，所述方法包括：根据获取的待处理地震数据，通过C-C算法计算得到相空间重构参数：嵌入维度和延迟时间；根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构；计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值；基于所述邻域点的均值进行局部投影，得到所述待处理地震数据的噪声分量。2.根据权利要求1所述的方法，其特征在于，还包括：基于所述噪声分量和所述待处理地震数据计算得到修正值，对所述待处理地震数据进行处理，以去除噪声。3.根据权利要求1所述的方法，其特征在于，所述根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构，具体包括：利用所述嵌入维度m和所述延迟时间τ，将所述待处理地震数据一维的时间序列重构到多维的相空间中，得到一个由N个相点组成的一个m维相空间。4.根据权利要求1所述的方法，其特征在于，所述计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值，具体包括：计算所述待处理地震数据的方差，得到初始邻域半径值r；根据所述的嵌入维度m进行迭代，判断是否满足迭代终止条件n≥m，其中n表示邻域内的相点数，m表示邻域内相点最小值；如果满足条件，则终止迭代计算，得到该次迭代对应的邻域半径值r；根据得到的所述邻域半径值，计算得到每个相点的邻域Ui＝{Xj；||Xj-Xi||<r}，Xj表示第j个相点，Xi表示第i个相点；根据所述每个相点的邻域计算得到所有邻域点的均值5.根据权利要求1所述的方法，其特征在于，所述基于所述邻域点的均值进行局部投影，计算得到所述待处理地震数据的噪声分量，具体包括：根据预先设定的矩阵和由该矩阵确定的对角权重矩阵R确定协方差矩阵C；计算所述协方差矩阵C的特征向量aq以及特征值λq；将所述相空间重构的各个相点投影到由所述协方差矩阵的特征向量组成的相空间中，得到噪声分量θ。6...

【专利技术属性】
技术研发人员：钟庆良，陈双全，莫莉，张颖燕，曾晓平，孟会杰，
申请(专利权)人：中国石油化工股份有限公司江汉油田分公司物探研究院，
类型：发明
国别省市：湖北,42