OBN数据去噪方法及相关装置制造方法及图纸

本发明专利技术公开了一种OBN数据去噪方法及相关装置。所述方法包括：分别对P分量数据和Z分量数据进行双树复小波变换；根据所述P分量数据的振幅谱振幅最大值对所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值；计算得到振幅匹配后的Z分量数据的振幅谱；根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数；根据所述振幅匹配系数，计算得到振幅匹配后的Z分量数据复小波系数；将所述振幅匹配后的Z分量数据复小波系数进行双树复小波逆变换，得到去噪后的Z分量数据。本发明专利技术能够明显压制Z分量数据中的Vz噪声，去噪后Z分量道集数据及叠加剖面的信噪比得到明显提高。噪比得到明显提高。噪比得到明显提高。

【技术实现步骤摘要】
OBN数据去噪方法及相关装置


[0001]本专利技术涉及物探数据处理领域，特别涉及一种OBN数据去噪方法及相关装置。

技术介绍

[0002]海底节点(Ocean bottom nodes，OBN)是一种放置在水底的四分量检波器，震源在海面，各个采集节点之间相互独立。OBN放置在水底，其位置固定，受海况影响相对较小，采集到的地震资料具有全方位、宽频带等特点。OBN节点由四个分量的接收器组成，其中一个分量接收器是水检压力检波器，用来记录标量波场(P分量数据)；另外三个相互垂直分量的接收器是陆检速度检波器，用来记录矢量波场(X、Y、Z分量数据)。由于水、陆双检检波器的工作原理不同，水检检波器对其周围水的压力非常敏感，而陆检检波器对其所在位置波的传播速度非常敏感，因此即使是同一地震波场，二者的响应是不同。P分量数据由于仅记录标量波场，因此其信噪比较高，但由于受到海底地质条件及检波器与海底耦合因素影响，到达检波器的地震波场的传播方向与检波器的接收方向不一致，会导致Z分量数据因接收了泄露的横波能量而混叠入低速、低频的噪声，即Vz噪声，信噪比差。如不对其进行处理，将对后续偏移成像处理产生不利影响。Vz噪声具有形态近似双曲、低频、强振幅且与有效信号的位置及频带皆有混叠的特征。目前去除Vz噪声的方法包括：利用X和Y分量数据对Z分量数据进行自适应相减方法，KL变换去噪方法，分频噪声压制方法。

技术实现思路

[0003]专利技术人发现，现有技术手段难以达到对Z分量数据有效去噪，利用X和Y分量数据对Z分量数据进行自适应相减方法的去噪效果不明显，KL变换去噪方法仅能去除线性干扰噪声，分频噪声压制方法为单道处理手段，对信号和噪声的分解性较差，去噪效果有限。为了至少部分地解决现有技术存在的技术问题，专利技术人做出本专利技术，通过具体实施方式，提供一种OBN数据去噪方法及相关装置。
[0004]第一方面，本专利技术实施例提供一种OBN数据去噪方法，包括：
[0005]分别对P分量数据和Z分量数据进行双树复小波变换，分别得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱；
[0006]确定所述P分量数据的振幅谱振幅最大值及最大值在双树复小波域中对应的位置；
[0007]将所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，与所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值；
[0008]根据所述P分量数据的振幅谱振幅最大值、与所述P分量匹配的Z分量数据的振幅谱振幅值和所述P分量数据的振幅谱，计算得到振幅匹配后的Z分量数据的振幅谱；
[0009]根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数；
[0010]根据所述振幅匹配系数，计算得到振幅匹配后的Z分量数据复小波系数；
[0011]将所述振幅匹配后的Z分量数据复小波系数进行双树复小波逆变换，得到去噪后
的Z分量数据。
[0012]在一些可选的实施例中，所述分别对P分量数据和Z分量数据进行双树复小波变换之前，包括：
[0013]获取OBN的水检压力检波器在激发震源后收集并记录的P分量数据；以及OBN的陆检速度检波器在激发震源后收集并记录的Z分量数据。
[0014]在一些可选的实施例中，所述分别对P分量数据和Z分量数据进行双树复小波变换，分别得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱，包括：
[0015]按照预设的分级级数，分别对所述P分量数据和Z分量数据进行对应分级级数的双树复小波变换，得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱。
[0016]在一些可选的实施例中，所述将所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，与所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值，包括：
[0017]根据所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，确定相同位置的Z分量数据的振幅谱振幅值，作为与所述P分量匹配的Z分量数据的振幅谱振幅值。
[0018]在一些可选的实施例中，根据所述P分量数据的振幅谱振幅最大值、与所述P分量匹配的Z分量数据的振幅谱振幅值和所述P分量数据的振幅谱，计算得到振幅匹配后的Z分量数据的振幅谱，包括：
[0019]将与所述P分量匹配的Z分量数据的振幅谱振幅值，除以所述P分量数据的振幅谱振幅最大值，乘以所述P分量数据的振幅谱，得到振幅匹配后的Z分量数据的振幅谱。
[0020]在一些可选的实施例中，所述根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数，包括：
[0021]将所述振幅匹配后的Z分量数据的振幅谱，除以相同位置的所述Z分量数据的振幅谱，得到所述振幅匹配后的Z分量数据的振幅谱中各位置的振幅匹配系数。
[0022]在一些可选的实施例中，所述根据所述振幅匹配系数，计算得到振幅匹配后的Z分量数据复小波系数，包括：
[0023]将所述振幅匹配后的Z分量数据的振幅谱每一位置的振幅匹配系数，分别乘以对应位置的所述Z分量数据复小波系数的实部和虚部，分别得到振幅匹配后的Z分量数据复小波系数的实部和虚部。
[0024]在一些可选的实施例中，所述将所述振幅匹配后的Z分量数据复小波系数进行双树复小波逆变换，得到去噪后的Z分量数据,包括：
[0025]按照所述Z分量数据双树复小波变换的分级级数，将所述振幅匹配后的Z分量数据复小波系数逐级进行双树复小波逆变换，得到去噪后的Z分量数据。
[0026]第二方面，本专利技术实施例提供一种OBN数据去噪装置，包括：
[0027]正变换模块，用于分别对P分量数据和Z分量数据进行双树复小波变换，分别得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱；
[0028]振幅匹配模块，用于确定所述P分量数据的振幅谱振幅最大值及最大值在双树复小波域中对应的位置；将所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，与所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值；根据所述P分量数据的振幅谱振幅最大值、与所述P分量匹配的Z分量数据的振幅谱振幅
值和所述P分量数据的振幅谱，计算得到振幅匹配后的Z分量数据的振幅谱；根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数；根据所述振幅匹配系数，计算得到振幅匹配后的Z分量数据复小波系数；
[0029]逆变换模块，用于将所述振幅匹配后的Z分量数据复小波系数进行双树复小波逆变换，得到去噪后的Z分量数据。
[0030]在一些可选的实施例中，还包括：
[0031]P和Z分量数据收集模块，用于获取OBN的水检压力检波器在激发震源后收集并记录的P分量数据；以及OBN的陆检速度检波器在激发震源后收集并记录的Z分量数据。
[0032]基于同一专利技术构思，本专利技术实施例还提供一种计算机存储介质，所述计算机存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器执行时实本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种OBN数据去噪方法，其特征在于，包括：分别对P分量数据和Z分量数据进行双树复小波变换，分别得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱；确定所述P分量数据的振幅谱振幅最大值及最大值在双树复小波域中对应的位置；将所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，与所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值；根据所述P分量数据的振幅谱振幅最大值、与所述P分量匹配的Z分量数据的振幅谱振幅值和所述P分量数据的振幅谱，计算得到振幅匹配后的Z分量数据的振幅谱；根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数；根据所述振幅匹配系数，计算得到振幅匹配后的Z分量数据复小波系数；将所述振幅匹配后的Z分量数据复小波系数进行双树复小波逆变换，得到去噪后的Z分量数据。2.如权利要求1所述的方法，其特征在于，所述分别对P分量数据和Z分量数据进行双树复小波变换之前，包括：获取OBN的水检压力检波器在激发震源后收集并记录的P分量数据；以及OBN的陆检速度检波器在激发震源后收集并记录的Z分量数据。3.如权利要求1所述的方法，其特征在于，所述分别对P分量数据和Z分量数据进行双树复小波变换，分别得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱，包括：按照预设的分级级数，分别对所述P分量数据和Z分量数据进行对应分级级数的双树复小波变换，得到所述P分量数据和Z分量数据在双树复小波域中的振幅谱。4.如权利要求1所述的方法，其特征在于，所述将所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，与所述Z分量数据的振幅谱进行匹配，得到与所述P分量匹配的Z分量数据的振幅谱振幅值，包括：根据所述P分量数据的振幅谱振幅最大值在双树复小波域中对应的位置，确定相同位置的Z分量数据的振幅谱振幅值，作为与所述P分量匹配的Z分量数据的振幅谱振幅值。5.如权利要求1所述的方法，其特征在于，根据所述P分量数据的振幅谱振幅最大值、与所述P分量匹配的Z分量数据的振幅谱振幅值和所述P分量数据的振幅谱，计算得到振幅匹配后的Z分量数据的振幅谱，包括：将与所述P分量匹配的Z分量数据的振幅谱振幅值，除以所述P分量数据的振幅谱振幅最大值，乘以所述P分量数据的振幅谱，得到振幅匹配后的Z分量数据的振幅谱。6.如权利要求1所述的方法，其特征在于，所述根据所述振幅匹配后的Z分量数据的振幅谱，计算得到振幅匹配系数，包括：将所述振幅匹配后...

【专利技术属性】
技术研发人员：张立彬，王伟，胡鑫，王克斌，王永明，张胜，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：