两界面之间的多次波的叠后预测方法和装置制造方法及图纸

本发明专利技术提供了一种两界面之间的多次波的叠后预测方法和装置，其中，该方法包括：获取给定的第一时深和第二时深，其中，第二时深大于第一时深，第一时深对应的地下反射界面为第一界面，第二时深对应的地下反射界面为第二界面；以第一界面将叠加剖面分为第一部分数据和第二部分数据；根据第一界面的全程和微屈多次波计算得到第一界面多次波的地震道的能量，根据第二界面的全程和微屈多次波计算得到第二界面多次波的地震道的能量，并计算得到归一化的互相关剖面；根据归一化的相关剖面和第二部分数据，计算得到第一界面和第二界面之间的层间多次波。本发明专利技术上述实施例达到了有效提高预测准确性和预测效率的技术效果。

【技术实现步骤摘要】
两界面之间的多次波的叠后预测方法和装置
本专利技术涉及地球物理勘探
，特别涉及一种两界面之间的多次波的叠后预测方法和装置。
技术介绍
在一些特殊的岩体与围岩之间的界面往往会产生较强的多次波或层间多次波，例如：碳酸盐岩、火成岩和煤层等发育的探区就很容易产生较强的多次波或层间多次波，这会对地震资料的速度分析、成像和反演产生比较大的影响。为了减少这种影响，一般需要对多次波进行压制，目前对于信噪比较高的地震资料，压制多次波一般采用叠前方法，而对于信噪比低的地震资料，压制多次波一般采用叠后方法。进一步的，考虑到海上地震资料信噪比较高，陆上低信噪比地震资料信噪比低，因此压制海上地震资料的多次波一般采用叠前方法，压制陆上地震资料多次波一般采用叠后方法，下面对叠后方法进行具体描述：叠后压制多次波的方法主要采用波场延拓的方法，即给定地下一个或多个产生多次波的界面及对应的叠加速度，然后通过一次下行波波场延拓和一次上行波波场延拓来预测与该界面有关的全程多次波和微屈多次波。对于两个界面之间的层间多次波预测，理论上可以通过六次波场延拓来实现，然而，这种实现方式的缺点是需要两个界面之间的层速度，而由叠加速度转化为层速度的过程会使得层速度的结果存在较大的误差，且，较多次的波场延拓会使得计算效率较低。针对如何高效准确地预测两个界面之间的层间多次波的问题，目前尚未提出有效的解决方案。
技术实现思路
本专利技术实施例提供了一种两界面之间的多次波的叠后预测方法，以解决现有技术中两界面的层间多次波的预测的准确性不高，预测效率低下的技术问题，该方法包括：获取叠加剖面和与所述叠加剖面对应的叠加速度场；获取给定的第一时深和第二时深，其中，第二时深大于第一时深，第一时深对应的地下反射界面为第一界面，第二时深对应的地下反射界面为第二界面；以所述第一界面将所述叠加剖面分为第一部分数据和第二部分数据，其中，所述第一部分数据为所述第一界面以上部分的数据，所述第二部分数据为所述第一界面以下部分的数据；从所述叠加速度场中查找出所述第一时深对应的第一叠加速度和所述第二时深对应的第二叠加速度；根据所述第一时深和所述第一叠加速度计算得到所述第一界面的第一深度，根据所述第二时深和所述第二叠加速度计算得到所述第二界面的第二深度；根据所述第一深度和所述第一叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第一界面，得到第一波场，再采用上行波方程，将所述第一波场从所述第一界面延拓到所述基准面，得到所述第一界面的全程和微屈多次波；根据所述第二深度和所述第二叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第二界面，得到第二波场，再采用上行波方程，将所述第二波场从所述第二界面延拓到所述基准面，得到所述第二界面的全程和微屈多次波；根据所述第一界面的全程和微屈多次波计算得到所述第一界面多次波的地震道的能量，根据所述第二界面的全程和微屈多次波计算得到所述第二界面多次波的地震道的能量；根据所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，计算得到所述第一界面和所述第二界面的多次波之间的互相关剖面；根据所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量和所述互相关剖面，计算得到归一化的互相关剖面；根据所述归一化的相关剖面和所述第二部分数据，计算得到所述第一界面和所述第二界面之间的层间多次波。在一个实施方式中，按照以下公式计算所述第一深度和第二深度：Z1(x,y)＝T1(x,y)V1(x,y)Z2(x,y)＝T2(x,y)V2(x,y)其中，V1(x,y)＝V(x,y,T1(x,y))，V2(x,y)＝V(x,y,T2(x,y))，Z1表示第一深度，Z2表示第二深度，T1表示第一时深，T2表示第二时深，V1表示第一叠加速度，V2表示第二叠加速度，V表示所述叠加剖面的叠加速度场，x表示叠加剖面的纵测线方向坐标，y表示叠加剖面的横测线方向坐标。在一个实施方式中，按照以下公式计算所述第一界面多次波的地震道的能量和所述第二界面多次波的地震道的能量：其中，M1(x,y,t)表示第一界面的全程和微屈多次波，M2(x,y,t)表示第二界面的全程和微屈多次波，E1(x,y)表示第一界面多次波的地震道的能量，E2(x,y)表示第二界面多次波的地震道的能量，x表示叠加剖面的纵测线方向坐标，y表示叠加剖面的横测线方向坐标，t表示时间坐标，tm表示叠加剖面的最大记录时间。在一个实施方式中，按照以下公式计算所述第一界面和所述第二界面的多次波之间的互相关剖面：其中，R(x,y,t)表示所述第一界面和所述第二界面的多次波之间的互相关剖面，M1(x,y,t)表示第一界面的全程和微屈多次波，M2(x,y,t)表示第二界面的全程和微屈多次波，表示互相关运算。在一个实施方式中，按照以下公式计算归一化的互相关剖面：其中，Rn(x,y,t)表示归一化的互相关剖面。在一个实施方式中，按照以下公式计算所述第一界面和所述第二界面之间的层间多次波：M3(x,y,t)＝Rn(x,y,t)*D2(x,y,t)其中，M3(x,y,t)表示所述第一界面和所述第二界面之间的层间多次波，D2(x,y,t)表示所述第二部分数据，*表示褶积运算。本专利技术实施例还提供了一种两界面之间的多次波的叠后预测装置，以解决现有技术中两界面的层间多次波的预测的准确性不高，预测效率低下的技术问题，该装置包括：叠加剖面和叠加速度场获取单元，用于获取叠加剖面和与所述叠加剖面对应的叠加速度场；时深获取单元，用于获取给定的第一时深和第二时深，其中，第二时深大于第一时深，第一时深对应的地下反射界面为第一界面，第二时深对应的地下反射界面为第二界面；叠加剖面划分单元，用于以所述第一界面将所述叠加剖面分为第一部分数据和第二部分数据，其中，所述第一部分数据为所述第一界面以上部分的数据，所述第二部分数据为所述第一界面以下部分的数据；查找单元，用于从所述叠加速度场中查找出所述第一时深对应的第一叠加速度和所述第二时深对应的第二叠加速度；深度计算单元，用于根据所述第一时深和所述第一叠加速度计算得到所述第一界面的第一深度，根据所述第二时深和所述第二叠加速度计算得到所述第二界面的第二深度；第一延拓单元，用于根据所述第一深度和所述第一叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第一界面，得到第一波场，再采用上行波方程，将所述第一波场从所述第一界面延拓到所述基准面，得到所述第一界面的全程和微屈多次波；第二延拓单元，用于根据所述第二深度和所述第二叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第二界面，得到第二波场，再采用上行波方程，将所述第二波场从所述第二界面延拓到所述基准面，得到所述第二界面的全程和微屈多次波；地震道的能量计算单元，用于根据所述第一界面的全程和微屈多次波计算得到所述第一界面多次波的地震道的能量，根据所述第二界面的全程和微屈多次波计算得到所述第二界面多次波的地震道的能量；互相关剖面计算单元，用于根据所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，计算得到所述第一界面和所述第二界面的多次波之间的互相关剖面；归一化单元，用于根据所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量和所述互相关剖面本文档来自技高网...

【技术保护点】
一种两界面之间的多次波的叠后预测方法，其特征在于，包括：获取叠加剖面和与所述叠加剖面对应的叠加速度场；获取给定的第一时深和第二时深，其中，第二时深大于第一时深，第一时深对应的地下反射界面为第一界面，第二时深对应的地下反射界面为第二界面；以所述第一界面将所述叠加剖面分为第一部分数据和第二部分数据，其中，所述第一部分数据为所述第一界面以上部分的数据，所述第二部分数据为所述第一界面以下部分的数据；从所述叠加速度场中查找出所述第一时深对应的第一叠加速度和所述第二时深对应的第二叠加速度；根据所述第一时深和所述第一叠加速度计算得到所述第一界面的第一深度，根据所述第二时深和所述第二叠加速度计算得到所述第二界面的第二深度；根据所述第一深度和所述第一叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第一界面，得到第一波场，再采用上行波方程，将所述第一波场从所述第一界面延拓到所述基准面，得到所述第一界面的全程和微屈多次波；根据所述第二深度和所述第二叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第二界面，得到第二波场，再采用上行波方程，将所述第二波场从所述第二界面延拓到所述基准面，得到所述第二界面的全程和微屈多次波；根据所述第一界面的全程和微屈多次波计算得到所述第一界面多次波的地震道的能量，根据所述第二界面的全程和微屈多次波计算得到所述第二界面多次波的地震道的能量；根据所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，计算得到所述第一界面和所述第二界面的多次波之间的互相关剖面；根据所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量和所述互相关剖面，计算得到归一化的互相关剖面；根据所述归一化的相关剖面和所述第二部分数据，计算得到所述第一界面和所述第二界面之间的层间多次波。...

【技术特征摘要】
1.一种两界面之间的多次波的叠后预测方法，其特征在于，包括：获取叠加剖面和与所述叠加剖面对应的叠加速度场；获取给定的第一时深和第二时深，其中，第二时深大于第一时深，第一时深对应的地下反射界面为第一界面，第二时深对应的地下反射界面为第二界面；以所述第一界面将所述叠加剖面分为第一部分数据和第二部分数据，其中，所述第一部分数据为所述第一界面以上部分的数据，所述第二部分数据为所述第一界面以下部分的数据；从所述叠加速度场中查找出所述第一时深对应的第一叠加速度和所述第二时深对应的第二叠加速度；根据所述第一时深和所述第一叠加速度计算得到所述第一界面的第一深度，根据所述第二时深和所述第二叠加速度计算得到所述第二界面的第二深度；根据所述第一深度和所述第一叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第一界面，得到第一波场，再采用上行波方程，将所述第一波场从所述第一界面延拓到所述基准面，得到所述第一界面的全程和微屈多次波；根据所述第二深度和所述第二叠加速度，采用下行波方程将所述叠加剖面从基准面延拓到所述第二界面，得到第二波场，再采用上行波方程，将所述第二波场从所述第二界面延拓到所述基准面，得到所述第二界面的全程和微屈多次波；根据所述第一界面的全程和微屈多次波计算得到所述第一界面多次波的地震道的能量，根据所述第二界面的全程和微屈多次波计算得到所述第二界面多次波的地震道的能量；根据所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，计算得到所述第一界面和所述第二界面的多次波之间的互相关剖面；根据所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量和所述互相关剖面，计算得到归一化的互相关剖面；根据所述归一化的相关剖面和所述第二部分数据，计算得到所述第一界面和所述第二界面之间的层间多次波。2.如权利要求1所述的方法，其特征在于，按照以下公式计算所述第一深度和第二深度：Z1(x,y)＝T1(x,y)V1(x,y)Z2(x,y)＝T2(x,y)V2(x,y)其中，V1(x,y)＝V(x,y,T1(x,y))，V2(x,y)＝V(x,y,T2(x,y))，Z1表示第一深度，Z2表示第二深度，T1表示第一时深，T2表示第二时深，V1表示第一叠加速度，V2表示第二叠加速度，V表示所述叠加剖面的叠加速度场，x表示叠加剖面的纵测线方向坐标，y表示叠加剖面的横测线方向坐标。3.如权利要求1或2所述的方法，其特征在于，按照以下公式计算所述第一界面多次波的地震道的能量和所述第二界面多次波的地震道的能量：1其中，M1(x,y,t)表示第一界面的全程和微屈多次波，M2(x,y,t)表示第二界面的全程和微屈多次波，E1(x,y)表示第一界面多次波的地震道的能量，E2(x,y)表示第二界面多次波的地震道的能量，x表示叠加剖面的纵测线方向坐标，y表示叠加剖面的横测线方向坐标，t表示时间坐标，tm表示叠加剖面的最大记录时间。4.如权利要求3所述的方法，其特征在于，按照以下公式计算所述第一界面和所述第二界面的多次波之间的互相关剖面：其中，R(x,y,t)表示所述第一界面和所述第二界面的多次波之间的互相关剖面，M1(x,y,t)表示第一界面的全程和微屈多次波，M2(x,y,t)表示第二界面的全程和微屈多次波，表示互相关运算。5.如权利要求4所述的方法，其特征在于，按照以下公式计算归一化的互相关剖面：其中，Rn(x,y,t)表示归一化的互相关剖面。6.如权利要求5所述的方法，其特征在于，按照以下公式计算所述第一界面和所述第二界面之间的层间多次波：M3(x,y,t)＝Rn(x,y,t)*D2(x,y,t)其中，M3(x,y,t)表示所述第一...

【专利技术属性】
技术研发人员：柯本喜，吴艳辉，李鹏，方云峰，杨继云，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，中油油气勘探软件国家工程研究中心有限公司，
类型：发明
国别省市：河北;13