[发明专利]地球物理勘探密度信息提取的方法和装置

说明书

本发明公开了一种地球物理勘探密度信息提取的方法和装置。其中，该方法包括：获取目标区的目标地震剖面，其中，目标地震剖面用于表征目标区的地层界面密度变化规律；通过反演模型对目标地震剖面进行叠后反演，得到目标区的地层密度剖面，其中，反演模型基于预设的密度测井数据来建立。本发明解决了现有技术中难以获取地层密度信息的技术问题。

技术领域

本发明涉及地球物理勘探数据处理领域，具体而言，涉及一种地球物理勘探密度信息提取。

背景技术

纵波速度(Vp)、横波速度(Vs)、密度(ρ)是表征储层的三个基本地震弹性参数，其中，密度与储层孔隙度，饱和度、泥质含量相关系数最高，能够最直观地反映储层和流体的变化规律。

目前，叠前反演是获取上述参数的最有效手段，其理论基础为Knott-Zoeppritz方程，由于该方程形式复杂，物理意义不明确，学者们进行了大量的近似处理，其中Aki-Richards近似(1980年)运用最为广泛，确切地表达了PP波、PS波反射振幅与弹性参数之间的关系。对于纵波(PP波)而言，纵波AVO特征主要取决于地层Vp、Vs、ρ三者的变化率，当地震偏移距有限(近道、中道)时，密度项对反射系数的贡献权重较小，在实际反演过程中，通常采用速度——密度的拟合关系(如Gardner公式、各工区拟合公式)进行密度项替换或直接省略密度项，将纵波三参数反演降维成两参数反演，这种近似处理明显地改善了中小偏移距道集的AVO反演稳定性，不仅提高了反演效率，而且使反演方程的病态性大大削弱，但直接损失就是无法实现精确的密度参数求解。

实现高精度的密度反演将对叠前地震数据用于勘探阶段储层精细描述、开发阶段油气藏动态监测产生深远的影响，但由于密度反演存在严重的病态性，密度参数一直难以获取。

针对现有技术中难以获取地层密度信息的问题，目前尚未提出有效的解决方案。

发明内容

本发明实施例提供了一种地球物理勘探密度信息提取的方法和装置，以至少解决现有技术中难以获取地层密度信息的技术问题。

根据本发明实施例的一个方面，提供了一种地球物理勘探密度信息提取的方法，包括：获取目标区的目标地震剖面，其中，目标地震剖面用于表征目标区的地层界面密度变化规律；通过反演模型对目标地震剖面进行叠后反演，得到目标区的地层密度剖面，其中，反演模型基于预设的密度测井数据来建立。

进一步地，获取目标区中角道集的最优叠加区间；通过对最优叠加区间的角道集进行叠加，得到目标区的目标地震剖面。

进一步地，获取目标区中角道集的最大叠加区间和最大叠加区间对应的地震信息；根据目标区中角道集的最大叠加区间和最大叠加区间对应的地震信息确定最优叠加区间。

进一步地，获取目标区地震波入射临界角和入射角角度改变量；根据入射临界角和入射角角度改变量确定最大叠加区间。

进一步地，对目标区进行层位标定，得到目标层段；根据目标层段的横波速度、纵波速度以及第一关系，得到地震波的入射临界角，其中，第一关系用于表征在剪切模量项系数为零的情况下，入射角与纵横波速度比的关系。

进一步地，根据目标层段的横波速度、纵波速度以及第二关系，得到入射角角度改变量，其中，入射角角度改变量包括第一入射角角度改变量和第二入射角角度改变量，第二关系为M((θ₀-Δθ₁),γ)+M((θ₀+Δθ₂),γ)＝0，其中，M用于表征剪切模量项系数，θ₁用于表征第一入射角，θ₂用于表征第二入射角，γ用于表征纵波速度与横波速度之比。

进一步地，获取目标区的地震波数据，其中，地震波数据包括：角道集数据、纵波速度、横波速度和测井解释结果；对地震波数据进行层位标定。