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

权利要求书

1.一种深度域地震流体分析方法，其特征在于，包括：

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

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

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

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

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

依据各个所述子空间扩散矩阵建立褶积模型中的整体空间扩散矩阵；

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

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

2.根据权利要求1所述的深度域地震流体分析方法，其特征在于，所述固液解耦流体因子反射系数方程为：

其中：

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

所述S(h)＝P(h)A(θ)M(h)为：

其中：

S(θ_i)＝[S(θ_i)₁ S(θ_i)₂ … S(θ_i)_N]^T，m_i＝[m_i(h₁) m_i(h₂) … m_i(h_N)]^T，

其中：

N为采样点个数，θ_i为第i个入射角，S(θ_i)_n为入射角为θ_i时第n个采样点的合成地震记录；m_i(h_n)为第n个采样点的第i个待反演参数反射系数，p(θ_i)_n为入射角为θ_i时第n个采样点的子空间扩散矩阵，a(θ_i)_n为入射角为θ_i时第n个采样点的系数，i＝1,2,3,4，n＝1,2,3,...,N。

3.根据权利要求1所述的深度域地震流体分析方法，其特征在于，所述对所述深度域合成地震记录进行反演分析，得到反演解析结果的过程为：

建立与所述深度域合成地震记录对应的目标函数；

对所述目标函数进行参数反演求解，得到与待反演参数对应的解析结果。

4.根据权利要求3所述的深度域地震流体分析方法，其特征在于，所述建立与所述深度域合成地震记录对应的目标函数的过程为：

采用贝叶斯反演方法建立与对所述深度域合成地震记录对应的目标函数。

5.根据权利要求4所述的深度域地震流体分析方法，其特征在于，所述对所述目标函数进行参数反演求解，得到与待反演参数对应的解析结果的过程为：

采用反复加权最小二乘法对所述目标函数进行参数反演求解，得到与待反演参数对应的解析结构。