[发明专利]一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法

权利要求书

1.一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法，其特征在于，包括如下步骤：

步骤100、基于纵波阻抗曲线和测井铀曲线计算拟阻抗曲线，计算公式为：

norm(PIp)＝a[norm(Ip)]+b[norm(-U)]

式中：norm(x)代表施加在x上的标准化函数；a和b为权系数，且a+b＝1；Ip和-U分别代表纵波阻抗曲线和测井铀曲线；PIp代表拟阻抗曲线，计算拟阻抗曲线的具体步骤为：

步骤101、将纵波阻抗曲线和测井铀曲线分别归一化至[0，1]；

步骤102、根据norm(PIp)＝a[norm(Ip)]+b[norm(-U)]测试确定权系数a、b，并获得归一化的拟阻抗曲线，所述权系数的确定原则为：

建立的拟阻抗曲线能够鉴别页岩品质；

在地震频带范围内，拟阻抗曲线和纵波阻抗曲线形成的地震响应相似；

步骤103、将归一化的拟阻抗曲线反归化到纵波阻抗曲线的范围内再获得拟阻抗曲线；

步骤200、对地震资料进行保幅处理，抽取高品质共反射点道集，并采用线性叠前弹性反演去除振幅随偏移距变化而变化因素的干扰获得纯纵波数据；

步骤300、结合拟阻抗曲线和纯纵波数据进行叠前反演获得拟阻抗数据体，并构建密度约束条件，基于去除了AVO效应的纯纵波数据，依照纵波阻抗和构建的拟阻抗曲线分别进行拟阻抗反演，获得稳定的三维空间拟阻抗数据，为叠前参数反演提供三维基础约束数据，并基于三维基础约束数据建立拟阻抗约束条件；

步骤400、基于提取的高品质共反射点道集和约束条件进行叠前参数的反演。

2.根据权利要求1所述的一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法，其特征在于，在步骤300中还包括建立纵波速度和入射角之间的关系式并进行线性化，并将各个曲线对应的参数作为加权项，在线性化的基础上确定各加权项的权系数，并基于各加权项和权系数获得弹性参数。

3.根据权利要求1所述的一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法，其特征在于，纯纵波的计算方法为：

叠前纵波反射系数公式为：

设定界面上下介质的纵横波速度和密度分别为α₁、α₂、β₁、β₂、ρ₁、ρ₂，θ₁，θ₂分别是反射角和透射角，则Δα＝α₂-α₁，Δβ＝β₂-β₁，Δρ＝ρ₂-ρ₁，α＝(α₁+α₂)/2，β＝(β₁+β₂)/2，θ＝(θ₁+θ₂)/2，ρ＝(ρ₁+ρ₂)/2，其中，R_p，R_s，R_d分别为纵波、横波反射系数和密度梯度；

设定

计算得到：

4.根据权利要求3所述的一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法，其特征在于，对地震资料进行保幅处理和抽取高品质共反射点道集后，设定叠前的道集为S(θ)，将S(θ)直接视为R_pp(θ)，此时已知泥石关系且道集个数大于三的情况下即可求得准确的Sp，其中，Sp为真正的纵波剖面。

5.根据权利要求4所述的一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法，其特征在于，对求取的Sp结合获得的弹性参数去除不同入射角度子波的影响。