[发明专利]叠前各向异性特征参数反演方法及计算机可读存储介质

说明书

公开了一种叠前各向异性特征参数反演方法及计算机可读存储介质。该方法可以包括：基于叠前地震道集，获得反射系数；针对叠前地震数据进行叠前同时反演，分别获得上层介质与下层介质的纵波速度、横波速度与密度；基于反射系数、上层介质与下层介质的纵波速度、横波速度与密度，根据各向异性参数反演公式进行反演，获得上层介质与下层介质的各向异性参数。在Ruger公式的基础上做了变换，以反映地层各向异性参数与裂缝发育之间的潜在联系。

技术领域

本发明涉及地震资料解释领域，更具体地，涉及一种叠前各向异性特征参数反演方法及计算机可读存储介质。

背景技术

基于纵波数据的方位各向异性AVO特征来预测裂缝是目前较为经济的、可行的方法之一。Zoeppritz方程作为AVO技术的理论基础，给出了水平层状介质中平面波传播的反射和折射规律，但由于其形式复杂，难以给出确定的物理意义，因此，难以直观表述振幅系数与介质参数之间的关系。因此，许多专家学者对Zoeppritz方程进行近似简化，得到了Aki、Shuey、Hilterman等多种近似公式，但这些都是以各向同性介质为前提。随着地层各向异性问题越来越受到人们的重视，讨论各向异性参数对反射系数的影响成为了必然的趋势。

1977年，Dely和Hron讨论了双层横向各向同性介质间界面处反射系数与透射系数的精确表达式。1987年，Wright分析了横向各向同性介质中各向异性对分界面处反射和折射所产生的影响。同年，Banik讨论了横向各向同性介质中各向异性参数的物理意义以及各向异性参数对分界面处反射和透射系数的影响。1993年，Kim建立了一个上层为各向异性页岩、下层为各向同性含气砂岩的模型，分析了分界面处各向异性对P波反射系数的影响。Ruger于1996年推导出了弱各向异性的横向各向同性介质中纵波反射系数的近似表达式。1998年，Vavrycuk V和Psencik I给出了任意各向异性介质中纵波反射系数随入射角、方位角变化的近似表达式。以往的研究都表明，介质的各向异性会严重影响各向同性介质模型假设下的AVO分析。

根据方位各向异性AVO特性有助于预测裂缝走向。1998年，Ruger给出了基于Ruger公式反演拟各向异性参数和裂缝方位角的推导，用拟各向异性参数表征裂缝发育密度，从而实现裂缝预测(Rüger A.Variation of P-wave reflectivity with offset andazimuth in anisotropic media[J].Geophysics,1998,63(3):935-947)。2004年，Gray给出了根据叠前方位角道集数据进行裂缝预测的应用实例，提供了利用地震方位各向异性AVO特征预测裂缝的思路和技术流程。2005年，朱兆林等人研究了裂缝介质的纵波方位各向异性AVO的反演方法，推导了基于Ruger公式的反演公式并用模型进行验证。2009年，肖鹏飞等人探讨了TTI介质中裂缝密度公式，进一步分析了裂缝倾角对密度反演的影响。2010年，Downton分析了当时方位AVO技术的局限性，并提出了一种新的方位各向异性弹性反演方法来克服这些局限性。在1998年Ruger给出的反演方法中，其反演结果为拟各向异性参数和方位角，而不是反演得到Thomsen各向异性参数。因此，有必要开发一种叠前地震各向异性特征参数反演方法，以反映地层各向异性参数与裂缝发育之间的潜在联系。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种叠前各向异性特征参数反演方法及计算机可读存储介质，其能够在Ruger公式的基础上做了变换，以反映地层各向异性参数与裂缝发育之间的潜在联系。