一种叠后地震数据处理方法、装置、电子设备和介质制造方法及图纸

本发明专利技术实施例提供了一种叠后地震数据处理方法、装置、电子设备和介质，所述方法包括：获取待处理数据；所述待处理数据包括与至少一个地层对应的叠后地震数据；确定目标地层；沿所述目标地层进行层拉平处理，得到拉平体数据；将所述拉平体数据进行连续小波变换拓频，得到拓频数据；对所述拓频数据进行反层拉平处理以及奇异值分解去噪处理，得到去噪数据；针对所述去噪数据进行各项异性拉普拉斯滤波处理，得到优化数据。本发明专利技术实施例能够对叠后地震数据进行优化处理，以提高数据的信噪比和分辨率。辨率。辨率。

【技术实现步骤摘要】
一种叠后地震数据处理方法、装置、电子设备和介质


[0001]本专利技术涉及石油天然气勘探开发
，特别是涉及一种叠后地震数据处理方法、装置、电子设备和介质。

技术介绍

[0002]对地区进行地震数据分析，有利于推进石油等资源勘探。但是部分工区地表类型复杂，影响了地震采集资料的整体信噪比和分辨率，另外，在斜坡区，地层产状较陡，较陡的产状会影响资料成像效果。
[0003]现有技术中，针对工区特殊的地表环境以及复杂的地质构造特征，从以下两个方面来提高地震资料的分辨率和信噪比：
①
在采集上，采用井炮、气枪、震源联合高密度采集的方式；
②
在处理上，通过子波一致性处理消除子波在频率、波形、相位的差异，改善子波的一致性，来保证地震资料的分辨率、信噪比。按照上述方法虽然对采集处理联合攻关后的地震资料品质有所改善，但由于原始地震资料受复杂采集地表激发接收条件以及地层构造的影响，整体信噪比、分辨率较低，不能满足开展精细构造解释及储层预测工作。
[0004]因此，在一定程度上使上述方法仍具备一定局限性，不适于普遍使用。

技术实现思路

[0005]鉴于上述问题，提出了本专利技术实施例以便提供一种克服上述问题或者至少部分地解决上述问题的一种叠后地震数据处理方法、装置、电子设备和介质。
[0006]为了解决上述问题，本专利技术实施例公开了一种叠后地震数据处理方法，包括：
[0007]获取待处理数据；所述待处理数据包括与至少一个地层对应的叠后地震数据；
>[0008]确定目标地层；
[0009]沿所述目标地层进行层拉平处理，得到拉平体数据；
[0010]基于连续小波变换将所述拉平体数据进行拓频，得到拓频数据；
[0011]对所述拓频数据进行反拉平处理以及奇异值分解去噪处理，得到去噪数据；
[0012]针对所述去噪数据进行各项异性拉普拉斯滤波处理，得到优化数据。
[0013]可选地，还包括：
[0014]基于所述优化数据进行反演，得到反演数据；
[0015]对比所述反演数据和预设的参照数据，确定所述优化数据的符合率。
[0016]可选地，所述将所述拉平体数据进行拓频，得到拓频数据包括：
[0017]对所述拉平体数据进行连续小波变换分解，得到分解数据；所述分解数据包括与多个频率对应的地震信号；
[0018]确定拓频参数；所述拓频参数包括拓频范围以及基准频率；
[0019]基于拓频参数和所述地震信号进行重新构建，得到拓频数据。
[0020]可选地，所述对所述拓频数据进行反拉平处理以及去噪处理，得到去噪数据包括：
[0021]对拓频数据进行反拉平处理，得到反拉平数据；
[0022]对所述反拉平数据进行奇异值分解去噪处理，得到去噪数据。
[0023]可选地，所述对所述反拉平数据进行奇异值分解去噪处理，得到去噪数据包括：
[0024]生成与所述反拉平数据匹配的Hankel矩阵；
[0025]对所述Hankel矩阵进行奇异值分解，得到一个或多个相互对应的奇异值和奇异值矢量；
[0026]按照预设规则确定目标奇异值；
[0027]采用所述目标奇异值以及所述目标奇异值对应的奇异值矢量，生成去噪数据。
[0028]可选地，所述针对所述去噪数据进行滤波处理，得到优化数据包括：
[0029]基于所述去噪数据，计算梯度结构张量；
[0030]依据所述梯度结构张量，确定地震剖面置信度信息；
[0031]依据所述梯度结构张量，地层倾角；
[0032]依据所述置信度信息和所述地层倾角，构建各向异性拉普拉斯滤波器；
[0033]采用所述各向异性拉普拉斯滤波器对所述去噪数据进行处理，得到优化数据；
[0034]其中，所述各向异性拉普拉斯滤波器中的滤波窗口的尺寸和形状由所述置信度信息确定，所述各向异性拉普拉斯滤波器中的窗户方向由所述地层倾角确定；所述置信度信息包括线型结构特征置信度和横向不连续性结构特征置信度。
[0035]本专利技术实施例还公开了一种叠后地震数据处理装置，包括：
[0036]待处理数据获取模块，用于获取待处理数据；所述待处理数据包括与至少一个地层对应的叠后地震数据；
[0037]目标地层确定模块，用于确定目标地层；
[0038]拉平体数据生成模块，用于沿所述目标地层进行层拉平处理，得到拉平体数据；
[0039]拓频数据生成模块，用于将所述拉平体数据进行连续小波变换拓频，得到拓频数据；
[0040]去噪数据生成模块，用于对所述拓频数据进行反拉平处理以及奇异值分解去噪处理，得到去噪数据；
[0041]滤波处理模块，用于针对所述去噪数据进行各项异性拉普拉斯滤波处理，得到优化数据。
[0042]可选地，还包括：
[0043]反演模块，用于基于所述优化数据进行反演，得到反演数据；
[0044]对比模块，用于对比所述反演数据和预设的参照数据，确定所述优化数据的符合率。
[0045]本专利技术实施例还公开了一种电子设备，包括：处理器、存储器及存储在所述存储器上并能够在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如上所述的叠后地震数据处理方法的步骤。
[0046]本专利技术实施例还公开了一种计算机可读存储介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被处理器执行时实现如上所述的叠后地震数据处理方法的步骤。
[0047]本专利技术实施例包括以下优点：
[0048]通过对叠后地震数据沿目标地层拉平后进行连续小波变换拓频处理得到拓频数
据，以拓宽地震数据频带；对拓频数据进行反拉平和奇异值分解去噪处理，以对拓频数据有效去除随机噪声得到去噪数据；针对去噪数数据进行各项异性拉普拉斯滤波处理，有效衰减随机噪声、增强了反射同相轴的横向一致性的同时，还能够保持剖面内的地层边缘和细微结构，通过上述处理后得到的优化数据，相对于初始的叠后地震数据，具备更好的信噪比和分辨率。
[0049]进一步的，通过依据优化数据进行反演处理，能够得出更准确的解释性数据，提高对探测区域地层预测准确性。
附图说明
[0050]图1是本专利技术的一种叠后地震数据处理方法实施例的步骤流程图；
[0051]图2是本专利技术实施例提供的一种降噪过程示意图；
[0052]图3a是本专利技术实施例提供的部分频段拓频前的剖面示意图；
[0053]图3b是本专利技术实施例提供的部分频段拓频后的剖面示意图；
[0054]图4a是本专利技术实施例提供的去噪处理和滤波处理前的剖面示意图；
[0055]图4b是本专利技术实施例提供的去噪处理和滤波处理后的剖面示意图；
[0056]图5a是本专利技术实施例提供的待处理数据的反演剖面示意图；
[0057]图5b是本专利技术实施例提供的优化数据的反演剖面示意图；...

【技术保护点】

【技术特征摘要】
1.一种叠后地震数据处理方法，其特征在于，包括：获取待处理数据；所述待处理数据包括与至少一个地层对应的叠后地震数据；确定目标地层；沿所述目标地层进行层拉平处理，得到拉平体数据；将所述拉平体数据进行连续小波变换拓频，得到拓频数据；对所述拓频数据进行反拉平处理以及奇异值分解去噪处理，得到去噪数据；针对所述去噪数据进行各项异性拉普拉斯滤波处理，得到优化数据。2.根据权利要求1所述的方法，其特征在于，还包括：基于所述优化数据进行反演，得到反演数据；对比所述反演数据和预设的参照数据，确定所述优化数据的符合率。3.根据权利要求1或2所述的方法，其特征在于，所述将所述拉平体数据进行连续小波变换拓频，得到拓频数据包括：对所述拉平体数据进行连续小波变换分解，得到分解数据；所述分解数据包括与多个频率对应的地震信号；确定拓频参数；所述拓频参数包括拓频范围以及基准频率；基于拓频参数和所述地震信号进行重新构建，得到拓频数据。4.根据权利要求1或2所述的方法，其特征在于，所述对所述拓频数据进行反拉平处理以及奇异值分解去噪处理，得到去噪数据包括：对拓频数据进行反拉平处理，得到反拉平数据；对所述反拉平数据进行奇异值分解去噪处理，得到去噪数据。5.根据权利要求4所述的方法，其特征在于，所述对所述反拉平数据进行奇异值分解去噪处理，得到去噪数据包括：生成与所述反拉平数据匹配的汉克尔Hankel矩阵；对所述Hankel矩阵进行奇异值分解，得到一个或多个相互对应的奇异值和奇异值矢量；按照预设规则确定目标奇异值；采用所述目标奇异值以及所述目标奇异值对应的奇异值矢量，生成去噪数据。6.根据权利要求1或2所述的方法，其特征在于，所述针对所述去噪数据进行各项异性拉普拉斯滤波处理，得到优化数据包括：基于所述去噪数据，计算梯度结构张量；依据所述梯度...

【专利技术属性】
技术研发人员：张珊，于宝利，邓勇，钟厚财，朱哲，刘鑫铭，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：