[发明专利]深度域地震子波的确定方法和装置

说明书

本申请实施方式提供了一种深度域地震子波的确定方法和装置，其中，该方法包括：获取声波测井数据，根据声波测井数据，分别确定深度域反射系数序列和井位置的深度域成像数据；根据所述深度域反射系数序列和所述井位置的深度域成像数据，确定深度域地震子波。由于该方案通过根据深度域反射系数序列和深度域成像数据，构建矩阵方程，求取以深度域成像数据为期望输出，基于深度域反射系数序列的最佳匹配算子作为深度域地震子波。因而解决了现有的深度域地震子波的确定方法存在的无法直接求取深度域地震子波、确定过程复杂、所确定的深度域地震子波准确度差的技术问题，达到了准确、简便地确定深度域地震子波的技术效果。

技术领域

本申请涉及油气勘探技术领域，特别涉及一种深度域地震子波的确定方法和装置。

背景技术

在油气勘探中，地震子波，尤其是深度域地震子波是一种重要的地震数据。通常，地震子波指的是一段具有确定的起始时间，能量有限，且具有一定延续长度的信号。它是地震记录中的基本单元。一般认为，地震震源激发时所产生的地震波仅是一个延续时间极短的尖脉冲，随着尖脉冲在粘弹性介质中传播，尖脉冲的高频成分会很快衰减，波形会随之增长，进而形成了地震子波，然后以地震子波的形式在地下传播。由于地震子波的上述特点，在具体实施时，经常需要使用深度域地震子波。例如，在正演问题的分析中，需要通过波动方程或褶积模型结合深度域地震子波产生模拟数据，进而可以得到合成记录，以合成记录作为指导，进行具体的油气勘探或者储层预测。因此，如何准确地确定深度域的地震子波一直是人们关注的一个问题。

现有的深度域地震子波的确定方法通常不能直接地提取得到深度域地震子波，常常是先获取时间域地震子波，再根据褶积定理将时间域地震子波转化为深度域地震子波。但是，具体实施时，由于深度域地震子波与反射系数的褶积并不满足褶积定理，即不能够满足子波随深度递变的特性。因此，现有的深度域地震子波的确定方法具体实施时，往往存在无法直接确定深度域地震子波、所确定的深度域地震子波准确度差、确定过程复杂的技术问题。

针对上述问题，目前尚未提出有效的解决方案。

发明内容

本申请实施方式提供了一种深度域地震子波的确定方法和装置，以解决现有的深度域地震子波确定方法中存在的无法直接确定深度域地震子波、所确定的深度域地震子波准确度差、确定过程复杂的技术问题。

本申请实施方式提供了一种深度域地震子波的确定方法，包括：

获取待测区域的声波测井数据；

根据所述声波测井数据，确定深度域反射系数序列和井位置的深度域成像数据；

根据所述深度域反射系数序列和所述井位置的深度域成像数据，确定深度域地震子波。

在一个实施方式中，根据声波测井数据，确定深度域反射系数序列，包括：

根据所述声波测井数据，获取速度场数据和密度场数据；

根据所述速度场数据和所述密度场数据，确定所述深度域反射系数序列。

在一个实施方式中，根据所述速度场数据和所述密度场数据，确定所述深度域反射系数序列，包括：按照以下公式确定所述深度域反射系数序列：

其中，R^d(z_i)为深度域反射系数序列中的第i个采样点的深度域反射系数，υ(z_i)为速度场数据中的第i个采样点的速度场数据，ρ(z_i)为密度场数据中的第i个采样点的密度场数据，z_i为测井数据中的第i个采样点的深度值，i为测井数据中深度方向上的采样点的编号。

在一个实施方式中，根据所述深度域反射系数序列和所述井位置的深度域成像数据，确定深度域地震子波，包括：