[发明专利]基于叠前反演的“甜点”概率直接估算方法及系统

权利要求书

1.基于叠前反演的“甜点”概率直接估算方法，其特征是，包括：

利用测井数据得到反演参数之间的联合概率分布，即纵波模量与密度的联合概率分布；

在反演参数初始模型的基础上，利用概率估计方法得到初始“甜点”相概率；

在初始“甜点”相概率的约束下，计算“甜点”敏感参数纵波模量及剪切模量的条件概率分布，并从条件概率分布中采样得到新的纵波模量与剪切模量模型；

然后根据新的纵波模量模型，结合纵波模量与密度的联合概率分布，计算密度的条件概率分布，并从条件概率分布中采样得到新的密度模型；

最后，计算新的“甜点”相概率，并利用接受概率公式来判断是否更新反演参数及“甜点”相概率模型，通过不断的迭代，即可得到有关反演参数及“甜点”相概率后验概率分布的一系列样本，最终实现“甜点”相概率的直接估算；

初始“甜点”相概率为：

式中，a，b为权重系数，且a+b＝1；

对于纵波模量与剪切模量，假设其在“甜点”相为k的条件概率p(Ι＝k|M)，p(Ι＝k|u)分布均服从高斯分布，即：

式中，是纵波模量与剪切模量在“甜点”相为k时的均值序列，与为其对应的方差序列。

2.如权利要求1所述的基于叠前反演的“甜点”概率直接估算方法，其特征是，纵波模量与密度的联合概率分布即是考虑两个参数之间相关性的特征,在已知一个参数的情况下，得到另一个参数的值域范围及概率分布特征。

3.如权利要求1所述的基于叠前反演的“甜点”概率直接估算方法，其特征是，纵波模量及剪切模量的条件概率分布p(M|I′)，p(μ|I′)是通过井数据统计得到，就是在甜点状态下，反演参数的先验分布特征，并将其定义为下个状态的参数采样空间。

4.如权利要求1所述的基于叠前反演的“甜点”概率直接估算方法，其特征是，密度的条件概率分布是拟合在给定纵波模量值时，所对应得密度取值得到。

5.如权利要求1所述的基于叠前反演的“甜点”概率直接估算方法，其特征是，在密度的条件概率分布中采样即可得到新的密度模型。

6.如权利要求1所述的基于叠前反演的“甜点”概率直接估算方法，其特征是，反演参数之间的联合概率分布的获得为：

包含纵波模量、剪切模量及密度的Zeoppritz近似公式表示为PP波反射系数；

将PP波反射系数与子波的褶积表示为地震响应，从而得到纵横波模量及密度与地震记录之间的正演关系，从而构建反演所需要的似然函数表达式。

7.基于叠前反演的“甜点”概率直接估算系统，其特征是，包括：

纵波模量与密度的联合概率分布计算模块：利用测井数据得到反演参数之间的联合概率分布，即纵波模量与密度的联合概率分布；

初始“甜点”相概率计算模块：在反演参数初始模型的基础上，利用概率估计方法得到初始“甜点”相概率；

新的纵波模量与剪切模量模型获得模块：在初始“甜点”相概率的约束下，计算“甜点”敏感参数纵波模量及剪切模量的条件概率分布，并从条件概率分布中采样得到新的纵波模量与剪切模量模型；

新的密度模型获得模块：然后根据新的纵波模量模型，结合纵波模量与密度的联合概率分布，计算密度的条件概率分布，并从条件概率分布中采样得到新的密度模型；

新的“甜点”相概率计算模块：计算新的“甜点”相概率，并利用接受概率公式来判断是否更新反演参数及“甜点”相概率模型，通过不断的迭代，即可得到有关反演参数及“甜点”相概率后验概率分布的一系列样本，最终实现“甜点”相概率的直接估算；

初始“甜点”相概率为：

式中，a，b为权重系数，且a+b＝1；

对于纵波模量与剪切模量，假设其在“甜点”相为k的条件概率p(Ι＝k|M)，p(Ι＝k|u)分布均服从高斯分布，即：

式中，是纵波模量与剪切模量在“甜点”相为k时的均值序列，与为其对应的方差序列。