[发明专利]高角度断裂约束的叠前裂缝预测方法、装置及存储介质

说明书

本发明涉及一种高角度断裂约束的叠前裂缝预测方法，包括如下步骤：利用高分辨率Radon变换技术，从叠后地震数据中提取地震剖面所反映的高角度断裂地震反射信息；采用三维高斯平滑与归一化处理将上述高分辨率三维Radon变换提取的高角度断裂估计转换为三维地震数据体中各采样点参与叠前裂缝预测的先验权重；建立方位傅里叶反射系数公式的矩阵‑向量化方程；将三维地震数据中各采样点参与叠前裂缝预测的先验权重改写为矩阵‑向量系统，进一步加入上述方位傅里叶反射系数公式的矩阵‑向量化方程中。

技术领域

本发明涉及油气田开发技术领域，具体是关于一种高角度断裂约束的叠前裂缝预测方法、装置及存储介质。

背景技术

裂缝发育程度是决定潜山类油气储层品质最核心的因素之一，现有裂缝预测方法可以分为依靠叠前地震资料的叠前裂缝预测与依靠叠后地震资料的叠后裂缝预测。叠前裂缝预测方法具有预测精度更高，同时能够预测裂缝走向的优势，而叠后裂缝预测方法则具有计算效率更快，同时预测结果与地下宏观断裂发育规律更为吻合的优势。上述叠前、叠后裂缝预测方法从不同尺度上表征着地下裂缝的发育规律与特征，单独依靠任何一种方法的裂缝预测结果均存在不确定性。

然而，现阶段缺乏能够有效结合叠前、叠后裂缝预测两者优势的裂缝预测方法和手段，这制约了裂缝预测结果的可靠性，进一步导致裂缝预测成果在实际生产项目中应用受限，而这种局限是由叠前预测方法或叠后预测方法本身的理论局限所导致的，并不能通过改进叠前或叠后裂缝预测方法来解决。

发明内容

针对上述问题，本发明的目的是提供一种高角度断裂约束的叠前裂缝预测方法、装置及存储介质，以解决上述技术问题。

为实现上述目的，本发明采取以下技术方案：

本发明所述的高角度断裂约束的叠前裂缝预测方法，包括如下步骤：

利用高分辨率Radon变换技术，从叠后地震数据中提取地震剖面所反映的高角度断裂地震反射信息；

采用三维高斯平滑与归一化处理将上述高分辨率三维Radon变换提取的高角度断裂估计转换为三维地震数据体中各采样点参与叠前裂缝预测的先验权重；

建立方位傅里叶反射系数公式的矩阵-向量化方程；

将三维地震数据中各采样点参与叠前裂缝预测的先验权重改写为矩阵-向量系统，进一步加入上述方位傅里叶反射系数公式的矩阵-向量化方程中。

所述的叠前裂缝预测方法，优选地，还包括如下步骤：

将带有先验权重约束的方位傅里叶反射系数公式的矩阵-向量化方程改写成最小二乘反演目标函数形式。

所述的叠前裂缝预测方法，优选地，还包括如下步骤：

将叠前宽方位地震道集数据与对应的叠后三维地震数据中各采样点参与叠前裂缝预测的先验权重，输入到上述带有先验权重约束的方位傅里叶反射系数裂缝预测目标函数中，基于最优化算法求解得到裂缝密度与裂缝走向估计。

所述的叠前裂缝预测方法，优选地，所述从叠后地震数据中提取地震剖面所反映的高角度断裂地震反射信息的具体方法为：

通过三维高分辨率Radon正变换沿着特定的空间路径将叠后地震数据变换到τ-p_x-p_y域中，再通过在τ-p_x-p_y域中压制非高角度反射信号，然后将只保留高角度反射信号的τ-p_x-p_y域阈值结果再通过反Radon变换至t-x-y域，从而得到地震剖面中包含的高角度断裂地震反射信息；

具体地，高分辨率三维Radon变换技术实质是对输入数据沿特定空间路径进行积分，从而将具有规律排列的信号分解为Radon域内的稀疏散点，以此实现信号识别与分离，三维线性Radon正变换可表达为：