地震数据处理方法和装置制造方法及图纸

本申请实施例提供了一种地震数据处理方法和装置，其中，该方法包括：获取目标区域的地震数据；从地震数据中提取第一地震子波；在时间域中对第一地震子波进行压缩处理，以获取第二地震子波；根据地震数据、第一地震子波、第二地震子波建立反演目标函数；根据反演目标函数，确定处理后的地震数据，由于该方案是基于反演的原理通过建立并根据反演目标函数来获取具有较高分辨率的处理后的地震数据，从而避免了由于使用反子波算子造成的数据处理不稳定，解决了现有方法中存在的地震数据处理的准确度差、可靠性低的技术问题，达到了稳定获取保幅、高分辨率的地震数据的技术效果。

【技术实现步骤摘要】
地震数据处理方法和装置
本申请涉及地球物理勘探
，特别涉及一种地震数据处理方法和装置。
技术介绍
在地球物理勘探领域，常常会采集目标区域的地震数据，以便后续可以利用上述地震数据对目标区域进行储层预测等具体的地震勘探。但是，由于在采集地震数据的过程中，地震波会受到大地滤波的作用，例如大地或其他传播介质会对地震波产生一种低通滤波效应，使得所采集得到的地震数据往往是低分辨率的波形数据。此外，又由于地下还存在各类的噪声干扰，会降低地震波的反射波的信噪比，进一步对所采集的地震数据的质量造成影响。目前，为了提高所采集的地震数据的分辨率，大多采用反褶积方法进行地震数据的处理，通常会先确定反子波算子，再利用反子波算子和地震数据做褶积处理。但是，对于上述方法，数据处理的效果在很大程度上依赖于所确定的反子波算子。具体的，例如反子波算子的微小扰动容易导致反褶积结果发生大幅度的抖动，从而影响基于反褶积方法的获取的处理后的地震数据的准确度。即，现有方法具体实施往往存在地震数据处理的准确度差、可靠性低的技术问题。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本申请实施例提供了一种地震数据处理方法和装置，以解决现有方法中存在的地震数据处理的准确度差、可靠性低的技术问题，达到了稳定获取保幅、高分辨率的地震数据的技术效果。本申请实施例提供了一种地震数据处理方法，包括：获取目标区域的地震数据；从所述地震数据中提取第一地震子波；在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波；根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，其中，所述反演目标函数为基于反演机理建立的用于进行高分辨率地震数据处理的函数；根据所述反演目标函数，确定处理后的地震数据。在一个实施方式中，根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，包括：按照以下公式，建立所述反演目标函数：其中，minimize表示求取最小值，表示求取L2范数，y表示反演目标函数的解析解，W为基于第一地震子波的第一地震子波矩阵，Wh为基于第二地震子波的第二地震子波矩阵，x为地震数据的列向量形式，μ为正则化参数。在一个实施方式中，根据所述第二地震子波，按照以下公式确定所述基于第二地震子波的第二地震子波矩阵：其中，Wh为基于第二地震子波的第二地震子波矩阵，wh(i)为第二地震子波中编号为i的时间点的振幅值，n为第二地震子波时间点的总数。在一个实施方式中，根据所述反演目标函数，确定处理后的地震数据，包括：根据所述反演目标函数，求取所述反演目标函数的导数为0时所对应的解析解，并将所述解析解作为所述处理后的地震数据。在一个实施方式中，按照以下公式确定所述反演目标函数的导数为0时所对应的解析解：y′＝(WTW+μI)-1WTd其中，(·)-1表示矩阵求逆，y′为反演目标函数的导数为0时所对应的解析解，W为基于第一地震子波的第一地震子波矩阵，WT为基于第一地震子波的第一地震子波矩阵转置，I为单位矩阵，d为模型向量，μ为正则化参数。在一个实施方式中，按照以下方式确定所述模型向量：d＝Wh·x其中，d为模型向量，Wh为基于第二地震子波的第二地震子波矩阵，x为地震数据的列向量。本申请实施例还提供了一种地震数据处理装置，包括：获取模块，用于获取目标区域的地震数据；提取模块，用于从所述地震数据中提取第一地震子波；压缩处理模块，用于在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波；建立模块，用于根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，其中，所述反演目标函数为基于反演机理建立的用于进行高分辨率地震数据处理的函数；求解模块，用于根据所述反演目标函数，确定处理后的地震数据。在本申请实施例中，由于是基于反演的原理通过建立并根据反演目标函数来获取具有较高分辨率的处理后的地震数据，从而解决了现有方法中存在的地震数据处理的准确度差、可靠性低的技术问题，达到能够稳定地获取保幅、高分辨率的地震数据的技术效果。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1是根据本申请实施方式提供的地震数据处理方法的处理流程图；图2是根据本申请实施方式提供的地震数据处理装置的组成结构图；图3是基于本申请实施例提供的地震数据处理方法的电子设备组成结构示意图；图4是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的某油田处理前的地震记录的示意图；图5是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的某油田处理后的地震记录的示意图；图6是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的油田A处理前的叠加地震记录的示意图；图7是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的油田A处理后的叠加地震记录的示意图；图8是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的油田C处理前的地震数据的示意图；图9是在一个场景示例中应用本申请实施例提供的地震数据处理方法和装置获得的油田C处理后的地震数据的示意图。具体实施方式为了使本
的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。考虑到现有的地震数据处理方法大多是采用反褶积方法进行具体的地震数据处理，即：先确定反子波算子，再利用反子波算子和地震数据做褶积处理。但是，上述方法的处理效果过于依赖于所确定的反子波算子，且处理效果对于反子波算子非常敏感，例如反子波算子出现微小的扰动就很容易导致反褶积结果发生大较幅度的抖动，进而影响基于反褶积方法所获取的处理后的地震数据的准确度，导致现有方法具体实施时往往存在地震数据处理的准确度差、可靠性低的技术问题。针对产生上述技术问题的根本原因，本申请考虑可以不采用需要依赖反子波算子的反褶积方法，而是从地震数据的反演机理出发，设计一套基于反演的数据处理方法，即：可以基于反演的原理通过建立并根据反演目标函数来获取具有较高分辨率的处理后的地震数据，从而避免了由于使用反子波算子造成的数据处理不稳定，解决了现有方法中存在的地震数据处理的准确度差、可靠性低的技术问题，达到了稳定获取保幅、高分辨率的地震数据的技术效果。基于上述思考思路，本申请实施例提供了一种地震数据处理方法。具体请参阅图1所示的根据本申请实施方式提供的地震数据处理方法的处理流程图。本申请实施例提供的地震数据处理方法，具体实施时，可以包括以下步骤。S11：获取目标区域的地震数据。在本实施方式中，通过深入分析现有的地震数据处理方法所采用的反褶积方法的实现机理，即：先确定反子波算子作为一种反子波滤波器，再利用这反子波滤波器和地震数据做褶积处理，以消除地震子波对地震数本文档来自技高网...

【技术保护点】
1.一种地震数据处理方法，其特征在于，包括：获取目标区域的地震数据；从所述地震数据中提取第一地震子波；在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波；根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，其中，所述反演目标函数为基于反演机理建立的用于进行高分辨率地震数据处理的函数；根据所述反演目标函数，确定处理后的地震数据。

【技术特征摘要】
1.一种地震数据处理方法，其特征在于，包括：获取目标区域的地震数据；从所述地震数据中提取第一地震子波；在时间域中对所述第一地震子波进行压缩处理，以获取第二地震子波；根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，其中，所述反演目标函数为基于反演机理建立的用于进行高分辨率地震数据处理的函数；根据所述反演目标函数，确定处理后的地震数据。2.根据权利要求1所述的方法，其特征在于，根据所述地震数据、所述第一地震子波、所述第二地震子波建立反演目标函数，包括：按照以下公式，建立所述反演目标函数：其中，minimize表示求取最小值，表示求取L2范数，y表示反演目标函数的解析解，W为基于第一地震子波的第一地震子波矩阵，Wh为基于第二地震子波的第二地震子波矩阵，x为地震数据的列向量形式，μ为正则化参数。3.根据权利要求2所述的方法，其特征在于，根据所述第二地震子波，按照以下公式确定所述基于第二地震子波的第二地震子波矩阵：其中，Wh为基于第二地震子波的第二地震子波矩阵，wh(i)为第二地震子波中编号为i的时间点的振幅值，n为第二地震子波时间点的总数。4.根据权利要求2所述的方法，其特征在于，根据所述反演目标函数，确定处理后的地震数据，包括：根据所述反演目标函数，求取所述反演目标函数的导数为0时所对应的...

【专利技术属性】
技术研发人员：杜昕，李皓，李国发，周振晓，张文波，唐博文，王伟，
申请(专利权)人：中国石油天然气集团有限公司，中国石油集团东方地球物理勘探有限责任公司，中国石油大学北京，
类型：发明
国别省市：北京,11