[发明专利]一种地震波精确反射系数解析方法及装置

说明书

本发明公开了一种地震波精确反射系数解析方法及装置，其中该方法包括：获取纵波入射角θ_p1、纵波透射角θ_p2、横波反射角θ_s1、横波透射角θ_s2、纵波速度v_p、横波速度v_s、纵波反射系数rp、横波反射系数rs及密度反射系数rd；按照下列公式对地震波反射系数进行解析：按照本申请公开的技术方案实现地震波反射系数的解析具有精度较高的特点，进一步，利用起来建立AVO正演模拟的过程能够避免对于常规AVO三参数反演的应用范围的限制，从而解决了现有技术中存在的问题。

技术领域

本发明涉及叠后反演技术领域，更具体地说，涉及一种地震波精确反射系数解析方法及装置。

背景技术

叠后反演利用叠后地震资料，主要反演地层纵波信息。叠前地震资料包含了许多叠后资料所不具备的信息，除反演纵波信息外，还可估计地层横波、密度、岩石模量、流体敏感参数等信息，所以叠前地震反演被广泛应于储层预测和流体识别。

常规AVO三参数反演通常利用Zoeppritz方程的近似式(Aki-Richards近似和Fatti近似等)来建立AVO正演模拟的过程，然而在P波入射角过临界角和弹性参数变化剧烈的情况下，近似公式精度有限，限制了常规AVO三参数反演的应用范围。

综上所述，如何提供一种精度较高的技术方案来建立AVO正演模拟的过程，是目前本领域技术人员亟待解决的问题。

发明内容

本发明的目的是提供一种地震波精确反射系数解析方法及装置，以提高地震波反射系数解析的精度。

为了实现上述目的，本发明提供如下技术方案：

一种地震波精确反射系数解析方法，包括：

获取纵波入射角θ_p1、纵波透射角θ_p2、横波反射角θ_s1、横波透射角θ_s2、纵波速度v_p、横波速度v_s、纵波反射系数rp、横波反射系数rs及密度反射系数rd；

按照下列公式对地震波反射系数进行解析：

其中，

优选的，获取纵波速度v_p及横波速度v_s，包括：

分别获取第一介质的纵波速度v_p1、横波速度v_s1，第二介质的纵波速度v_p2、横波速度v_s2，及第一介质和第二介质之间的纵波速度差Δv_p、横波速度差Δv_s；

按照下列公式分别计算纵波速度v_p及横波速度v_s：

优选的，获取纵波反射系数rp、横波反射系数rs及密度反射系数rd，包括：

分别获取所述第一介质的密度ρ₁、第二介质的密度ρ₁和第一介质和第二介质之间的密度差Δρ；

按照下列公式分别计算纵波反射系数rp、横波反射系数rs及密度反射系数rd：