[发明专利]地震数据处理方法和装置

说明书

技术领域

本发明涉及地震勘探资料处理领域，更具体地，涉及一种地震数据处理方法和一种地震数据处理装置。

背景技术

在中国南方和西部的勘探区块存在地表地形地貌、近地表岩性条件和复杂地下地质构造同时并存的工区，地形高差变化剧烈、低降速层厚度变化大、高陡岩层出露地表、地下逆掩推覆构造多，这些复杂的地震地质条件对现有的地震采集和处理技术提出了更大的要求。如此复杂的地质情况，采集到的地震信号信噪比非常低，常规的去噪方法和技术并不完全适用，具有很大的局限性和不适应性。

τ-p变换滤波技术是一种常用的地震数据处理手段，但是，在复杂山前带探区，有效信号往往被淹没在噪声中，看不到有效同相轴，或仅能看到极少数同相轴，信噪比极低(一般信噪比低于0.5)，利用现有τ-p变换滤波技术处理后得到的数据仍存在假同相轴及混叠现象，无法满足地震勘探的需要。

发明内容

本发明提出了一种加权τ-p变换滤波技术，其能够显著提高资料的信噪比。

根据本发明的一方面，提出了一种地震数据处理方法，该方法包括：得到t-x域的地震道集d(t,x)，其中，t表示双程旅行时，x表示偏移距；基于下式得到经加权τ-p变换后的τ-p域数据m(τ,p)：

其中，N表示所述地震道集中的道数，τ表示截距时间，p表示射线参数，x_n表示所述地震道集中第n道的偏移距，Δx表示相邻道的偏移距的差。

根据本发明的另一方面，提出了一种地震数据处理装置，该装置包括：地震道集获取单元，用于得到t-x域的地震道集d(t,x)，其中，t表示双程旅行时，x表示偏移距；加权τ-p变换单元，用于基于下式得到经加权τ-p变换后的τ-p域数据m(τ,p)：

其中，N表示所述地震道集中的道数，τ表示截距时间，p表示射线参数，x_n表示所述地震道集中第n道的偏移距，Δx表示相邻道的偏移距的差。

本发明中所采用的加权后的τ-p变换滤波技术能够大大降低低信噪比资料中存在的随机噪声，显著提高资料的信噪比，消除现有τ-p变换处理中存在的假同相轴及混叠现象，能够明显改善后续的偏移成像结果。本发明尤其适用于信噪比非常低的地震资料，一般指信噪比低于0.5，现有商业处理模块对此几乎无能为力，而应用本发明，能从这类极低信噪比资料中搜索到有效地震信号。

附图说明

通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。

图1示出了根据本发明的一个实施例的地震数据处理方法的流程图。

图2(a)为信噪比为0.6时应用现有τ-p变换得到的滤波结果的速度剖面和时间剖面；图2(b)为对相同资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图3(a)为信噪比为0.2时应用现有τ-p变换得到的滤波结果的速度剖面和时间剖面；图3(b)为对相同资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图4(a)为信噪比为1的原始地震资料的速度剖面和时间剖面；图4(b)为对图4(a)所示的资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图5(a)为信噪比为0.6的原始地震资料的速度剖面和时间剖面；图5(b)为对图5(a)所示的资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图6(a)为信噪比为0.4的原始地震资料的速度剖面和时间剖面；图6(b)为对图6(a)所示的资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图7(a)为信噪比为0.2的原始地震资料的速度剖面和时间剖面；图7(b)为对图7(a)所示的资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图8(a)为信噪比为0.1的原始地震资料的速度剖面和时间剖面；图8(b)为对图8(a)所示的资料应用本发明得到的滤波结果的速度剖面和时间剖面。

图9(a)、(b)、(c)、(d)依次为原始CMP(共中心点)地震道集、应用本发明得到的滤波后的数据、被滤除的噪声的速度剖面、被滤除的噪声的时间剖面。

图10(a)、(b)、(c)为图9(a)、(b)、(d)对应的频谱。

图11(a)为去噪前库东资料平原区CMP道集；图11(b)为去噪后库东资料平原区CMP道集。