[发明专利]一种含气饱和度反演方法

摘要

本发明涉及地质勘探技术领域，具体的讲是一种含气饱和度反演方法。本发明方法通过加权统计的概念以及解析化的目标反演函数表达式，利用测井资料统计分析与统计岩石物理建模相结合的思路，推导了较传统概率统计类方法效率更高、不受弹性参数精度及含气敏感性差异影响的目标反演函数及高效求解策略，方法具有较高的实用性，能有效地运用于三维实际资料，识别出高含气饱和度气藏的位置。

主权项

1.一种含气饱和度反演方法，其特征在于包括，根据含气饱和度参数R_i先验分布特征，得到含气饱和度先验概率密度函数P(Ri)=Σn=1NrαnN(Ri;μRin,ΣRin);]]>上式中，N(·)表示高斯分布概率密度函数，N_r,α_n,和分别为与饱和度参数R_i有关的混合高斯分布部件数、权值、均值与方差；设各弹性参数变量m^k间相互条件独立，则似然函数P(m_j|R_i)为：P(mj1,mj2,...mjNm|Ri)=Πk=1NmP(mjk|Ri);]]>根据训练样本统计出所述似然函数P(m^k|R_i)的分布信息，得到所述似然函数P(m^k|R_i)的解析表达式：P(mk|Ri)=Σπ=1NcαπN(Ri;μmk|Riπ,Σmk|Riπ),]]>上式中，N(·)表示高斯分布概率密度函数，Nc,απ,和分别为与m^k、R_i有关的混合高斯分布部件数、权值、均值与方差；根据所述含气饱和度参数先验分布特征和所述似然函数P(m^k|R_i)的解析表达式以及引入权重系数W_k取对数运算得到反演目标函数的解析化表达式：R~=lnΣn=1NrαnN(Ri;μRin,ΣRin)+Σk=1NmWklnΣπ=1NcαπN(Ri;μmk|Riπ,Σmk|Ri1π);]]>上式中，ln表示取自然对数；输入弹性参数数据，取所述反演目标函数最大后验概率密度解作为反演结果R~final=argMaxRi∈R{lnΣn=1NrαnN(Ri;μRin,ΣRin)+Σk=1NmWklnΣπ=1NcαπN(Ri;μmkRiπ,ΣmkRiπ)},]]>上式中，表示返回极大值对应的R_i。