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

权利要求书

1.一种深度域地震子波的确定方法，其特征在于，包括：

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

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

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

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

根据所述深度域反射系数序列，确定深度域反射系数序列的自相关参数；

根据所述深度域反射系数序列和所述井位置的深度域成像数据，确定深度域反射系数序列和深度域成像数据的互相关参数；

根据所述自相关参数和所述互相关参数，确定匹配算子，并将所述匹配算子作为所述深度域地震子波，其中，深度域成像数据为地震成像数据。

2.根据权利要求1所述的方法，其特征在于，根据声波测井数据，确定深度域反射系数序列，包括：

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

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

3.根据权利要求2所述的方法，其特征在于，根据所述速度场数据和所述密度场数据，确定所述深度域反射系数序列，包括：

按照以下公式确定所述深度域反射系数序列：

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

4.根据权利要求1所述的方法，其特征在于，根据所述深度域反射系数序列，确定深度域反射系数序列的自相关参数，包括：按照以下公式确定所述自相关参数：

其中，为第i个采样点的自相关参数，R^d(n)为深度为n时的深度域反射系数，R^d(n+z_i)为深度为第i个采样点的深度和n之和时的深度域反射系数，n为深度变量，取值为1至l范围内的自然数，z_i为测井数据中的第i个采样点的深度值，l为算子的长度。

5.根据权利要求1所述的方法，其特征在于，根据所述深度域反射系数序列和所述井位置的深度域成像数据，确定深度域反射系数序列和深度域成像数据的互相关参数，包括：按照以下公式确定所述互相关参数：

其中，μ_Rx(z_i)为第i个采样点的互相关参数，x^d(n)为深度为n时的成像数据，R^d(n+z)为深度为第i个采样点的深度和n之和时的深度域反射系数，n为深度变量，取值为1至l范围内的自然数，z_i为测井数据中的第i个采样点的深度值，l为算子的长度。

6.根据权利要求1所述的方法，其特征在于，根据所述自相关参数和所述互相关参数，确定匹配算子，包括：

根据所述自相关参数和所述互相关参数，构建求解矩阵方程；

通过求解所述求解矩阵方程，确定所述匹配算子。

7.根据权利要求6所述的方法，其特征在于，根据所述自相关参数和所述互相关参数，构建求解矩阵方程，包括：

根据所述自相关参数和所述互相关参数，构建以下方程作为所述求解矩阵方程：

其中，为深度为z时的自相关参数，μ_Rx(z)为深度为z时的互相关参数，w^d(z)为深度为z时的匹配算子，l为匹配算子的长度，z为深度，z的取值为0至l范围内的自然数。