[发明专利]一种基于泊松比分解的K值鲁棒YPD叠前同时反演方法

权利要求书

1.一种基于泊松比分解的K值鲁棒YPD叠前同时反演方法，其特征在于，包括：

一、测井曲线地震尺度匹配处理的步骤；

二、实际叠前道集时深标定的步骤；

三、实际叠前道集保AVO处理的步骤：

四、K值扰动AVA鲁棒性分析的步骤，具体包括：

首先，选取井上目的层段的弹性参数构建双层介质模型，模拟对比Aki-Richards近似方程、Zong-YPD近似方程和基于泊松比分解的YPD近似方程的AVA曲线；

然后，扰动K值模拟AVA曲线的变化，以此分析各个近似方程对K值的确定度和鲁棒性，具体包括以下步骤：

①与纵、横波速度和地层密度变化相关的Aki-Richard反射系数近似方程为：

其中，θ为纵波入射角；纵波速度反射系数项为横波速度反射系数项为密度反射系数项为为相邻地层的纵波速度的变化量；△V_S为相邻地层的横波速度的变化量；△ρ为相邻地层的密度的变化量；为相邻地层的纵波速度的平均值；为相邻地层的横波速度的平均值；为相邻地层的地层密度的平均值；K为地层横纵波速度比的平方；

②与井上杨氏模量、泊松比和地层密度有关的Zong-YPD反射系数近似方程为：

其中，井上杨氏模量反射系数项为泊松比反射系数项为地层密度反射系数项为△E为相邻地层的杨氏模量的变化量；△υ为相邻地层的泊松比的变化量；为相邻地层的杨氏模量的平均值；为相邻地层的泊松比的平均值；

③将泊松比υ分解为两个新的泊松比变量：υ₁＝1-υ，

④纵波速度V_P由井上杨氏模量E和泊松比υ表示为代入分解后的υ₁和υ₂则表示为：

其中，ρ为地层密度；

⑤同样地，横波速度V_S由井上杨氏模量E和泊松比υ表示为代入分解后的泊松比υ₁和υ₂则表示为：

⑥将公式(6)和(7)代入公式(4)，得到基于泊松比分解的YPD反射系数近似方程为：

其中，杨氏模量反射系数项为泊松比υ₁的反射系数项为泊松比υ₂的反射系数项为地层密度反射系数项为为相邻地层的分解后的泊松比υ₁的平均值；相邻地层的分解后的泊松比υ₂的平均值；

⑦根据公式(4)、(5)和(8)，在前三步骤的数据基础上，即可进行AVO模拟，给定不同K值进行AVA扰动分析，作为后续YPD参数反演稳定性质控；

五、基于泊松比分解的K值鲁棒YPD叠前同时反演的步骤。

2.根据权利要求1所述的基于泊松比分解的K值鲁棒YPD叠前同时反演方法，其特征在于，所述测井曲线地震尺度匹配处理的步骤具体包括：

①针对井上目的层段的测井数据进行预处理，具体为对测井曲线进行环境校正，剔除泥岩或含煤地层井眼垮塌引起的测量异常值，其中测井曲线包括地层密度曲线、纵波时差曲线和横波时差曲线；

②当缺少横波测井资料时，使用环境校正后的测井曲线进行预测横波速度；

③使用Backus平均及滑动窗函数法地层块化平均，进行测井曲线尺度向地震尺度的匹配处理，确保Backus处理后的地层密度、声波资料能够有效反映地震尺度响应，且消除K值中的高频异常对后续计算分析的影响；

④在尺度匹配处理后测井曲线的数据基础上，通过公式(1)-(3)分别求取出井上杨氏模量、泊松比和K值曲线，并生成相应的测井曲线，用于后续分析和反演：

式中，E为井上杨氏模量；υ为泊松比；ρ为地层密度；V_S为横波速度；V_P为纵波速度。

3.根据权利要求1所述的基于泊松比分解的K值鲁棒YPD叠前同时反演方法，其特征在于，所述实际叠前道集时深标定的步骤为：

结合校验炮或零偏VSP测井的时深关系，进行叠后井震标定，以此时深关系对实际叠前道集进行井震标定，以避免二类AVO引起叠后剖面中“暗点”的漏标，并确定是否需要进行剩余时差校正：当叠前道集存在显著未拉平或抖动时，则需要进行剩余时差校正。