一种叠前深度偏移方法及装置制造方法及图纸

本申请实施例公开了一种叠前深度偏移方法及装置。所述方法包括：获取目的工区的时间域的叠前地震数据，基于叠前地震数据，分别确定目的工区中各个成像点位置的旅行时信息和补偿旅行时信息；对叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据；分别确定第一地震道数据对应的炮点位置到目的工区中第一成像点位置的总旅行时和总补偿旅行时，以及分别确定第一地震道数据对应的检波点位置到目的工区中第一成像点位置的总旅行时和总补偿旅行时；基于总旅行时和总补偿旅行时，对初始补偿地震道数据进行插值处理，得到第一地震道数据对应的目标补偿地震道数据。可以提高叠前深度偏移的运算效率。

【技术实现步骤摘要】
一种叠前深度偏移方法及装置
本申请涉及石油地球物理勘探地震数据处理
，特别涉及一种叠前深度偏移方法及装置。
技术介绍
地下地层介质可能存在非完全弹性介质，其对地震波具有吸收衰减作用，因此，在叠前深度偏移过程中考虑这种介质的吸收衰减作用，可以提高叠前深度偏移的成像结果的准确性、精度或分辨率。在叠前深度偏移过程中考虑非完全弹性介质对地震波的吸收衰减作用的方法常称为粘声介质叠前深度偏移方法。现有技术中常用的粘声介质叠前深度偏移方法的主要过程是：获取目的工区的叠前地震数据；根据叠前地震数据确定目的工区的速度场数据和品质因子场数据；根据品质因子场数据，计算目的工区中各个成像点处射线路径补偿旅行时；获取叠前地震数据对应的旅行时表和补偿旅行时表，以及叠前地震数据中的第一地震道数据；根据旅行时表和补偿旅行时表，对第一地震道数据进行补偿处理；根据补偿处理后的第一地震道数据，确定目的工区中第一成像点的第一成像结果；对于第一成像点对应的其他地震道数据重复上述步骤，以确定目的工区中该成像点对应的由不同地震道数据对应的成像结果生成的目标成像结果。对于一个成像点，其可以对应多个炮点或检波点，那么不同炮点或检波点对应的旅行时不相同，同时，不同炮点或检波点对应的补偿旅行时间也不同。对于叠前地震数据中的某一地震道数据，该地震道数据中不同采样时间对应的相位和振幅补偿因子也不同，因此，现有技术中常用的粘声介质叠前深度偏移方法需要先将第一地震道数据从时间域变换到频率域，然后计算第一地震道数据对应的所有频点的补偿结果，再将这些补偿结果进行叠加处理，得到第一成像点的第一成像结果，通过对叠前地震数据中所有地震道数据循环进行上述补偿过程，得到第一成像点的目标成像结果。专利技术人发现现有技术中至少存在如下问题：由于需要计算叠前地震数据中每一个地震道数据对应的各个频点的补偿结果，以及将这些补偿结果进行叠加处理，导致现有技术中常用的粘声介质叠前深度偏移方法的计算量较大，从而导致采用该方法进行叠前深度偏移的运算效率较低。
技术实现思路
本申请实施例的目的是提供一种叠前深度偏移方法及装置，以提高叠前深度偏移的运算效率。为解决上述技术问题，本申请实施例提供一种叠前深度偏移方法及装置是这样实现的：一种叠前深度偏移方法，包括：获取目的工区的时间域的叠前地震数据，基于所述叠前地震数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息；根据所述补偿旅行时信息和预设间隔时间，对所述叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据；根据所述旅行时信息和补偿旅行时信息，分别确定第一地震道数据对应的炮点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时，以及分别确定第一地震道数据对应的检波点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时；基于所述总旅行时和所述总补偿旅行时，对所述初始补偿地震道数据进行插值处理，得到第一地震道数据对应的目标补偿地震道数据。优选方案中，所述基于叠前地震数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息，包括：基于所述叠前地震数据，确定所述目的工区的速度场数据和品质因子场数据；根据所述速度场数据和所述品质因子场数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息。优选方案中，所述根据速度场数据和品质因子场数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息，包括：根据所述速度场数据，分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的旅行时；根据所述品质因子场数据，分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的补偿旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的补偿旅行时；基于所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的旅行时和补偿旅行时，以及所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的旅行时和补偿旅行时，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息。优选方案中，采用下述公式分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的补偿旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的补偿旅行时：公式中，T*(x)表示所述目的工区中一个炮点位置到所述目的工区中一个成像点位置x的补偿旅行时，或所述目的工区中一个检波点位置到所述目的工区中一个成像点位置x的补偿旅行时，v0表示地震波在所述目的工区中的传播速度，Q表示所述品质因子场中一个成像点位置处的品质因子，x表示所述目的工区中一个成像点位置，s表示所述目的工区中一个炮点位置分别到所述目的工区中一个成像点位置的射线路径，或所述目的工区中一个检波点位置分别到所述目的工区中一个成像点位置的射线路径。优选方案中，所述根据补偿旅行时信息和预设间隔时间，对所述叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据，包括：通过傅立叶变换将第一地震道数据从时间域转换到频率域，得到频率域的第一地震道数据；获取所述补偿旅行时信息中的最大补偿旅行时；利用所述预设间隔时间将所述最大补偿旅行时划分为等间隔的补偿旅行时；基于所述频率域的第一地震道数据和所述等间隔的补偿旅行时，确定第一地震道数据对应的初始补偿地震道数据。优选方案中，采用下述公式确定第一地震道数据对应的初始补偿地震道数据：公式中，表示第一地震道数据对应的初始补偿地震道数据，表示预设补偿滤波函数，i表示虚数单位，ω表示第一地震道数据的角频率，u(ω)表示所述频率域的第一地震道数据，为所述等间隔的补偿旅行时，t表示所述时间域中第一地震道数据的采样时间。优选方案中，所述根据旅行时信息和补偿旅行时信息，分别确定第一地震道数据对应的炮点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时，以及分别确定第一地震道数据对应的检波点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时，包括：根据所述旅行时信息和补偿旅行时信息，分别获取第一地震道数据对应的炮点位置和检波点位置到所述目的工区中第一成像点位置的旅行时，以及分别获取第一地震道数据对应的炮点位置和检波点位置到所述目的工区中第一成像点位置的补偿旅行时；将所述第一地震道数据对应的炮点位置到第一成像点位置的旅行时与第一地震道数据对应的检波点位置到第一成像点位置的旅行时相加，得到第一地震道数据对应的炮点位置和检波点位置到所述目的工区中第一成像点位置的总旅行时；将所述第一地震道数据对应的炮点位置到第一成像点位置的补偿旅行时与第一地震道数据对应的检波点位置到第一成像点位置的补偿旅行时相加，得到第一地震道数据对应的炮点位置和检波点位置到所述目的工区中第一成像点位置的总补偿旅行时。优选方案中，采用下述公式进行所述插值处理，得到第一地震道数据对应的目标补偿地震道数据：公式中，T和T*分别表示所述总旅行时和所述总补偿旅行时，和分别表示与所述总补偿旅行时T*紧邻的前后两个等间隔的补偿旅行时，表示第一地震道数据对应的目标补偿地震道数据，本文档来自技高网...

【技术保护点】
一种叠前深度偏移方法，其特征在于，包括：获取目的工区的时间域的叠前地震数据，基于所述叠前地震数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息；根据所述补偿旅行时信息和预设间隔时间，对所述叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据；根据所述旅行时信息和补偿旅行时信息，分别确定第一地震道数据对应的炮点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时，以及分别确定第一地震道数据对应的检波点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时；基于所述总旅行时和所述总补偿旅行时，对所述初始补偿地震道数据进行插值处理，得到第一地震道数据对应的目标补偿地震道数据。

【技术特征摘要】
1.一种叠前深度偏移方法，其特征在于，包括：获取目的工区的时间域的叠前地震数据，基于所述叠前地震数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息；根据所述补偿旅行时信息和预设间隔时间，对所述叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据；根据所述旅行时信息和补偿旅行时信息，分别确定第一地震道数据对应的炮点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时，以及分别确定第一地震道数据对应的检波点位置到所述目的工区中第一成像点位置的总旅行时和总补偿旅行时；基于所述总旅行时和所述总补偿旅行时，对所述初始补偿地震道数据进行插值处理，得到第一地震道数据对应的目标补偿地震道数据。2.根据权利要求1所述的一种叠前深度偏移方法，其特征在于，所述基于叠前地震数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息，包括：基于所述叠前地震数据，确定所述目的工区的速度场数据和品质因子场数据；根据所述速度场数据和所述品质因子场数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息。3.根据权利要求2所述的一种叠前深度偏移方法，其特征在于，所述根据速度场数据和品质因子场数据，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息，包括：根据所述速度场数据，分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的旅行时；根据所述品质因子场数据，分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的补偿旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的补偿旅行时；基于所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的旅行时和补偿旅行时，以及所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的旅行时和补偿旅行时，分别确定所述目的工区中各个成像点位置的旅行时信息和补偿旅行时信息。4.根据权利要求3所述的一种叠前深度偏移方法，其特征在于，采用下述公式分别计算所述目的工区中各个炮点位置到所述目的工区中各个成像点位置的补偿旅行时，以及分别计算所述目的工区中各个检波点位置到所述目的工区中各个成像点位置的补偿旅行时：公式中，T*(x)表示所述目的工区中一个炮点位置到所述目的工区中一个成像点位置x的补偿旅行时，或所述目的工区中一个检波点位置到所述目的工区中一个成像点位置x的补偿旅行时，v0表示地震波在所述目的工区中的传播速度，Q表示所述品质因子场中一个成像点位置处的品质因子，x表示所述目的工区中一个成像点位置，s表示所述目的工区中一个炮点位置分别到所述目的工区中一个成像点位置的射线路径，或所述目的工区中一个检波点位置分别到所述目的工区中一个成像点位置的射线路径。5.根据权利要求1所述的一种叠前深度偏移方法，其特征在于，所述根据补偿旅行时信息和预设间隔时间，对所述叠前地震数据中第一地震道数据进行补偿预处理，得到第一地震道数据对应的初始补偿地震道数据，包括：通过傅立叶变换将第一地震道数据从时间域转换到频率域，得到频率域的第一地震道数据；获取所述补偿旅行时信息中的最大补偿旅行时；利用所述预设间隔时间将所述最大补偿旅行时划分为等间隔的补偿旅行时；基于所述频率域的第一地震道数据和所述等间隔的补偿旅行时，确定第一地震道数据对应的初始补偿地震道数据。6.根据权利要求5所述的一种叠前深度偏移方法，其特征在于，采用下述公式确定第一地震道数据对应的初始补偿地震道数据：公式中，表示第一地震道数据对应的初始补偿地震道数据，表示预设补偿滤波函数，i表示虚数单位，ω表示第一地震道数据的角频率，u(ω)表示所述频率域的第一地震道数据，为所述等间隔的补偿旅行时，t表示所述时间...

【专利技术属性】
技术研发人员：张巍毅，张建磊，王成祥，罗国安，王硕，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11