[发明专利]盲反演方法及装置

说明书

本发明公开了一种盲反演方法及装置，该方法包括：利用稀疏脉冲反演确定地震记录的反射系数矩阵，将反射系数矩阵转换为反射系数向量，将构建的时间维度、线维度及道维度上的一阶差分矩阵和二阶差分矩阵分别转换为总变分正则化矩阵和吉洪诺夫正则化矩阵，根据反射系数向量和总变分正则化矩阵，及吉洪诺夫正则化矩阵构建波阻抗的目标函数，波阻抗包括横向突变波阻抗和横向渐变波阻抗；根据目标函数取得最小值时的横向突变波阻抗和横向渐变波阻抗确定地震记录的波阻抗。本发明利用总变分正则化和吉洪诺夫正则化分别对横向突变波阻抗和横向渐变波阻抗进行正则化约束，利用构建的目标函数求解波阻抗，可以提高波阻抗横向能量的连续性。

技术领域

本发明涉及油气田勘探技术领域，尤其涉及盲反演方法及装置。

背景技术

本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

波阻抗反演实际上是从地震剖面上消除子波的影响，留下反射系数，再由反射系数计算出能够反映地层物性变化的物理参数波阻抗。波阻抗反演在地震勘探领域非常重要，广泛用于地震解释及储层预测，同时在岩石物理方面也得到了很好的应用。

波阻抗反演问题通常会涉及到反卷积的问题，通过从地震记录中提取地震子波求解出反射系数，然后根据递推原理得到波阻抗。但是，由于子波的不确定性，导致反射系数在横向上连续性较差，同时递推的过程对于噪声干扰往往非常敏感，使得反演后的波阻抗在横向上存在能量变化大、连续性差的问题，在地震波阻抗剖面上就会出现“挂面条”现象，不利于储层预测分析。

因此，现有的波阻抗反演存在波阻抗横向能量连续性差的问题。

发明内容

本发明实施例提供一种盲反演方法，用以提高波阻抗横向能量的连续性，该方法包括：

利用稀疏脉冲反演确定地震记录的反射系数矩阵，将反射系数矩阵按照时间维度、线维度及道维度的顺序转换为反射系数向量；

将构建的时间维度、线维度及道维度上的一阶差分矩阵分别转换为时间维度、线维度及道维度上的总变分正则化矩阵；

将构建的时间维度、线维度及道维度上的二阶差分矩阵分别转换为时间维度、线维度及道维度上的吉洪诺夫正则化矩阵；

根据反射系数向量和总变分正则化矩阵，及吉洪诺夫正则化矩阵构建波阻抗的目标函数；所述波阻抗包括横向突变波阻抗和横向渐变波阻抗；

根据波阻抗的目标函数取得最小值时的横向突变波阻抗和横向渐变波阻抗确定所述地震记录的波阻抗；

根据时间维度、线维度及道维度上的一阶差分矩阵和时间维度、线维度及道维度上的单位矩阵，将时间维度、线维度及道维度上的一阶差分矩阵分别转换为时间维度、线维度及道维度上的总变分正则化矩阵，包括：

根据道维度上的单位矩阵和线维度上的单位矩阵，及时间维度上的一阶差分矩阵，确定时间维度上的总变分正则化矩阵；

根据道维度上的单位矩阵和线维度上的一阶差分矩阵，及时间维度上的单位矩阵，确定线维度上的总变分正则化矩阵；

根据道维度上的一阶差分矩阵和线维度上的单位矩阵，及时间维度上的单位矩阵，确定道维度上的总变分正则化矩阵。

本发明实施例还提供一种盲反演装置，用以提高波阻抗横向能量的连续性，该装置包括：

反射系数确定模块，用于利用稀疏脉冲反演确定地震记录的反射系数矩阵，将反射系数矩阵按照时间维度、线维度及道维度的顺序转换为反射系数向量；

总变分正则化转换模块，用于将构建的时间维度、线维度及道维度上的一阶差分矩阵分别转换为时间维度、线维度及道维度上的总变分正则化矩阵；