本申请提供的一种时深转换曲线估计方法、装置、存储介质和程序产品，通过获取初始时深转换曲线、井旁道地震信号和深度域测井波阻抗，并根据初始时深转换曲线和深度域测井波阻抗，得到初始时间域测井反射系数，然后根据初始时间域测井反射系数、井旁道地震信号和预设的陷波点算子，得到初始时间域测井反射系数对应的第一陷波点图谱和井旁道地震信号第二陷波点图谱。然后对第一陷波点图谱和第二陷波点图谱进行时间方向上的匹配，得到陷波点匹配结果，根据该陷波点匹配结果，得到时深转换曲线，不需要估计地震子波，有效提升时深转换曲线估计的准确性。计的准确性。计的准确性。
【技术实现步骤摘要】
时深转换曲线估计方法、装置、存储介质和程序产品
[0001]本申请涉及地震资料处理
，尤其涉及一种时深转换曲线估计方法、装置、存储介质和程序产品。
技术介绍
[0002]地震资料处理是利用数字计算机对野外地震勘探所获得的原始资料进行加工和改造，可以为地质解释提供高质量且直观可靠的地质信息。由于地震资料横向覆盖范围广，纵向分辨率低且有效频带有限，而测井资料纵向分辨率高，频带宽，但空间上为单点数据，所以将二者结合对高精度的地震资料处理具有重要意义。时深转换是地震资料与测井资料匹配的一种方法，即将深度域的测井资料转换到时间域，然后再进行地震资料与测井资料匹配。
[0003]现有技术中，通常使用井旁道地震信号对时深转换曲线进行估计，先利用测井反射系数和井旁道地震信号估计地震子波，得到对应的人工合成地震信号，然后对井旁道和人工合成的地震信号进行匹配，最后利用匹配结果估计出时深转换曲线。
[0004]但是，地震子波和时深转换曲线的估计是相互影响的，所以时深转换曲线的估计结果一般存在较大误差。
技术实现思路
[0005]本申请提供一种时深转换曲线估计方法、装置、存储介质和程序产品，用以解决现有技术中时深转换曲线的估计结果不准确的问题。
[0006]第一方面，本申请提供一种时深转换曲线估计方法，包括：
[0007]获取初始时深转换曲线、井旁道地震信号和深度域测井波阻抗；
[0008]根据所述初始时深转换曲线和所述深度域测井波阻抗，得到初始时间域测井反射系数；
[0009]根据所述初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到所述初始时间域测井反射系数对应的第一陷波点图谱和所述井旁道地震信号对应的第二陷波点图谱；
[0010]将所述第一陷波点图谱和所述第二陷波点图谱进行时间方向的匹配，得到陷波点匹配结果；
[0011]根据所述陷波点匹配结果，得到时深转换曲线。
[0012]可选的，所述根据所述初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到所述初始时间域测井反射系数对应的第一陷波点图谱和所述井旁道地震信号对应的第二陷波点图谱，包括：
[0013]分别对所述初始时间域测井反射系数和所述井旁道地震信号进行短时傅里叶变换，得到所述初始时间域测井反射系数对应的第一时频振幅谱和所述井旁道地震信号对应的第二时频振幅谱；
[0014]利用所述井旁道地震信号对应的时频振幅谱估计有效频带，获取所述有效频带内的第一时频振幅谱和第二时频振幅谱；
[0015]将所述第一时频振幅谱中的每个元素代入所述预设的陷波点算子的定义式中，得到所述第一陷波点图谱；
[0016]将所述第二时频振幅谱中的每个元素代入所述预设的陷波点算子的定义式中，得到所述第二陷波点图谱。
[0017]可选的，所述预设的陷波点算子的定义式如下：
[0018][0019]其中，所述陷波点算子的窗口大小为(2w1+1)(2w2+1)，a
i，j
为所述第一时频振幅谱或者所述第二时频振幅谱中的元素，w1和w2为表征窗口大小的参数，表示时频振幅谱中的元素a
i，j
映射后得到的值。
[0020]可选的，所述将所述第一陷波点图谱和所述第二陷波点图谱进行时间方向的匹配，得到陷波点匹配结果，包括：
[0021]根据如下公式计算成本矩阵：
[0022][0023]其中，C为大小为n
×
n的成本矩阵，c
i，j
为所述成本矩阵中的元素，其中，w为控制窗口大小的参数，是所述第一陷波点图谱中的多个时刻对应的列向量组成的二维图像块，为所述第二陷波点图谱中的多个时刻对应的列向量组成的二维图像块，表示两个陷波点图谱中对应的二维图像块的结构相似性；
[0024]根据所述成本矩阵，计算得到距离矩阵，根据所述距离矩阵得到最短路径；
[0025]根据所述最短路径，得到所述陷波点匹配结果，所述陷波点匹配结果为新的陷波点图谱。
[0026]可选的，所述根据所述陷波点匹配结果，得到时深转换曲线，包括：
[0027]将所述初始时深转换曲线中的时间序列更新为所述新的陷波点图谱中的新时间序列，根据所述新时间序列和所述新时间序列对应的深度，得到所述时深转换曲线。
[0028]可选的，所述根据所述初始时深转换曲线和所述深度域测井波阻抗，得到初始时间域测井反射系数，包括：
[0029]根据所述初始时深转换曲线和所述深度域测井波阻抗，得到时间域测井波阻抗；
[0030]根据所述时间域测井波阻抗，得到所述初始时间域测井反射系数。
[0031]第二方面，本申请提供一种时深转换曲线估计装置，包括：
[0032]获取模块，用于获取初始时深转换曲线、井旁道地震信号和深度域测井波阻抗；
[0033]第一处理模块，用于根据所述初始时深转换曲线和所述深度域测井波阻抗，得到初始时间域测井反射系数；
[0034]第二处理模块，用于根据所述初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到所述初始时间域测井反射系数对应的第一陷波点图谱和所述井旁道地震信号对应的第二陷波点图谱；
[0035]匹配模块，用于将所述第一陷波点图谱和所述第二陷波点图谱进行时间方向的匹配，得到陷波点匹配结果；
[0036]第三处理模块，用于根据所述陷波点匹配结果，得到时深转换曲线。
[0037]第三方面，本申请提供一种电子设备，包括：处理器，以及与所述处理器通信连接的存储器；
[0038]所述存储器存储计算机执行指令；
[0039]所述处理器执行所述存储器存储的计算机执行指令，以实现第一方面所述的方法。
[0040]第四方面，本申请提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现第一方面所述的一种时深转换曲线估计方法。
[0041]第五方面，本申请提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时第一方面所述的方法。
[0042]本申请提供的一种时深转换曲线估计方法、装置、存储介质和程序产品，通过获取初始时深转换曲线、井旁道地震信号和深度域测井波阻抗，并根据初始时深转换曲线和深度域测井波阻抗，得到初始时间域测井反射系数，然后根据初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到初始时间域测井反射系数对应的第一陷波点图谱和井旁道地震信号第二陷波点图谱。然后对第一陷波点图谱和第二陷波点图谱进行时间方向上的匹配，得到陷波点匹配结果，根据该陷波点匹配结果，得到时深转换曲线，不需要估计地震子波，有效提升时深转换曲线估计的准确性。
附图说明
[0043]此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。
[0044]图1为本申请实本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种时深转换曲线估计方法，其特征在于，包括：获取初始时深转换曲线、井旁道地震信号和深度域测井波阻抗；根据所述初始时深转换曲线和所述深度域测井波阻抗，得到初始时间域测井反射系数；根据所述初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到所述初始时间域测井反射系数对应的第一陷波点图谱和所述井旁道地震信号对应的第二陷波点图谱；将所述第一陷波点图谱和所述第二陷波点图谱进行时间方向的匹配，得到陷波点匹配结果；根据所述陷波点匹配结果，得到时深转换曲线。2.根据权利要求1所述的方法，其特征在于，所述根据所述初始时间域测井反射系数、所述井旁道地震信号和预设的陷波点算子，得到所述初始时间域测井反射系数对应的第一陷波点图谱和所述井旁道地震信号对应的第二陷波点图谱，包括：分别对所述初始时间域测井反射系数和所述井旁道地震信号进行短时傅里叶变换，得到所述初始时间域测井反射系数对应的第一时频振幅谱和所述井旁道地震信号对应的第二时频振幅谱；利用所述井旁道地震信号对应的时频振幅谱估计有效频带，获取所述有效频带内的第一时频振幅谱和第二时频振幅谱；将所述第一时频振幅谱中的每个元素代入所述预设的陷波点算子的定义式中，得到所述第一陷波点图谱；将所述第二时频振幅谱中的每个元素代入所述预设的陷波点算子的定义式中，得到所述第二陷波点图谱。3.根据权利要求2所述的方法，其特征在于，所述预设的陷波点算子的定义式如下：其中，所述陷波点算子的窗口大小为(2w1+1)(2w2+1)，a
i，j
为所述第一时频振幅谱或者所述第二时频振幅谱中的元素，w1和w2为表征窗口大小的参数，表示时频振幅谱中的元素a
i，j
映射后得到的值。4.根据权利要求3所述的方法，其特征在于，所述将所述第一陷波点图谱和所述第二陷波点图谱进行时间方向的匹配，得到陷波点匹配结果，包括：根据如下公式计算成本矩阵：其中，C为大小为n
×
n的成本矩阵，c
i，j
为所述成本矩阵中的元素，其中，w为控制窗口大小的参数，是所述第一陷波点图谱中的多个时刻对应的列向量组成的二维图像块，为所述第
【专利技术属性】
技术研发人员：陆文凯，王佳林，
申请(专利权)人：清华大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。