[发明专利]一种地震数据的噪声去除方法及装置

说明书

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

技术领域

本发明涉及地震勘探技术领域，具体而言，涉及一种地震数据的噪声去除方法及装置。

背景技术

随着技术形势的发展，勘探目标逐渐向深层和地表复杂地区转移，给地震资料处理工作提出了更高的要求，并且地震数据处理一直以高信噪比、高分辨率、高保真为目标，因此更为精确的计算地震数据中的噪声成为地震数据处理的关键。

现有技术中，传统的地震数据噪声去除方法计算复杂程度较高，对地震数据中噪声的计算准确度不高且计算速度较慢。

发明内容

有鉴于此，本发明的目的在于提供一种地震数据的噪声去除方法和装置，以提高计算的准确度以及运算处理速度。

第一方面，本发明实施例提供了一种地震数据的噪声去除方法，包括：

根据获取的待处理地震数据，通过C-C算法计算得到相空间重构参数：嵌入维度和延迟时间；

根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构；

计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值；

基于所述邻域点的均值进行局部投影，得到所述待处理地震数据的噪声分量。

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所述方法，还包括：

基于所述噪声分量和所述待处理地震数据计算得到修正值，对所述待处理地震数据进行处理，以去除噪声。

结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述根据所述嵌入维数和所述延迟时间对所述待处理地震数据进行相空间重构，具体包括：

利用所述嵌入维度m和所述延迟时间τ，将所述待处理地震数据一维的时间序列重构到多维的相空间中，得到一个由N个相点组成的一个m维相空间。

结合第一方面，本发明实施例提供了第一方面的第三种可能的实施方式，其中，所述计算相空间每个相点的邻域以及每个邻域内所有邻域点的均值，具体包括：

计算所述待处理地震数据的方差，得到初始邻域半径值r；

根据所述的嵌入维度m进行迭代，判断是否满足迭代终止条件n≥m，其中n表示邻域内的相点数，m表示邻域内相点最小值；

如果满足条件，则终止迭代计算，得到该次迭代对应的邻域半径值r；

根据得到的所述邻域半径值，计算得到每个相点的邻域U_i＝{X_j；||X_j-X_i||<r}，X_j表示第j个相点，X_i表示第i个相点；

根据所述每个相点的邻域计算得到所有邻域点的均值

结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，所述基于所述邻域点的均值进行局部投影，计算得到所述待处理地震数据的噪声分量，具体包括：

根据预先设定的矩阵和由该矩阵确定的对角权重矩阵R确定协方差矩阵C；