[发明专利]一种深度域地震流体分析方法、装置、系统及存储介质

说明书

本发明公开了一种深度域地震流体分析方法、装置、系统及计算机可读存储介质，该方法包括分别获取在不同入射角下各个采样点在成像过程中的脉冲数据；依据与每个入射角对应的各个脉冲数据，建立在相应的入射角下、与每个所述采样点一一对应的子空间扩散矩阵；依据各个所述子空间扩散矩阵建立深度域的褶积模型，以得到深度域合成地震记录；对所述深度域合成地震记录进行反演分析，得到反演解析结果。本发明基于不同入射角下每个采样点对应的脉冲数据建立不同入射角下每个采样点的子空间扩散矩阵，进而建立深度域的褶积模型，使由该褶积模型得到的深度域合成地震记录能够更好地反应地下介质的实际情况，有利于提高深度域地震流体识别的精确度。

技术领域

本发明实施例涉及油气勘探开发技术领域，特别是涉及一种深度域地震流体分析方法、装置、系统及计算机可读存储介质。

背景技术

近年来，地震反演技术已经广泛应用于石油勘探开发的各个阶段，是储层预测的核心技术，能够将测井纵向上高分辨率和地震横向上密集采样的优势结合起来，估算出地层的岩性、物性特征信息横向上的变化。目前，常规的地震反演是在时间域内进行的，由于测井曲线从深度域转换到时间域损失了宝贵的高频信息，所以时间域地震反演的局限性是对薄层的识别能力较差。随着勘探开发的深入，时间域的地震反演由于自身的局限性经常无法满足储层精细描述的要求，因此深度域地震反演应运而生。

目前，现有的关于深度域地震反演的方法是利用速度替换法完成深度域褶积的反演方法，其主要是利用伪深度域的思想，将深度域资料调整到某一恒定速度的伪深度域来满足线性深不变条件，提取伪深度域子波以合成深度域地震合成记录，然后分析确定深度域反演中的各个参数，得出深度域反演的最终结果。

但是，该方法假设地下介质为水平层状介质，在构造较为复杂、速度变化大的地区，采用现有技术中的方法所提取的深度域子波不能够准确表示地下介质的实际情况，无法满足实际生产的需要，并且在不同域间相互转化会造成累积误差，导致深度域反演的最终结果不准确。

发明内容

本发明实施例的目的是提供一种深度域地震流体分析方法、装置、系统及计算机可读存储介质，在使用过程中能够使反演解析结果更加准确，有利于提高深度域地震流体识别的精确度。

为解决上述技术问题，本发明实施例提供了一种深度域地震流体分析方法，包括：

分别获取在不同入射角下各个采样点在成像过程中的脉冲数据；

依据与每个所述入射角对应的各个脉冲数据，建立在相应的入射角下、与每个所述采样点一一对应的子空间扩散矩阵；

依据各个所述子空间扩散矩阵建立深度域的褶积模型，以得到深度域合成地震记录；

对所述深度域合成地震记录进行反演分析，得到反演解析结果。

可选的，所述依据各个所述子空间扩散矩阵建立深度域的褶积模型的过程为：

依据各个所述子空间扩散矩阵建立深度域的整体空间扩散矩阵；

依据所述整体空间扩散矩阵及固液解耦流体因子反射系数方程建立深度域的褶积模型，其中，所述褶积模型为S(h)＝P(h)A(θ)M(h)，其中：

S(h)为深度域合成地震记录，P(h)为整体空间扩散矩阵，A(θ)为系数矩阵，M(h)为待反演参数反射系数矩阵。

可选的，所述固液解耦流体因子反射系数方程为：

其中：

K_f为固液解耦流体因子，φ为孔隙度，μ为剪切模量，ρ为密度，V_p表示地震波纵波速度，V_s表示地震波横波速度，γ_dry表示干岩纵横波速度比，γ_sat表示湿岩纵横波速度比，θ为入射角，Δ表示梯度；