[发明专利]多波动态匹配方法及系统

说明书

公开了一种多波动态匹配方法及系统。该方法包括：步骤1：输入初始纵横波速比，获得初始纵横波速比变化率；步骤2：进行时间匹配，获得纵波和匹配后的转换波；步骤3：对纵波和转换波进行反演，根据反演属性值计算基于反演的纵横波速比变化率；步骤4：将基于反演的纵横波速比变化率与初始纵横波速比变化率做差，获得纵横波速比变化率差值；步骤5：根据纵横波速比变化率差值，调整初始纵横波速比，重复步骤1‑4，直到纵横波速比变化率差值最小，将此时的转换波作为最终转换波。本发明通过动态调整纵横波速比，解决了输入波速比值与反演得到的波速比值的不一致问题，实现了纵横波的自动匹配，提高了匹配的可信度。

技术领域

本发明涉及地震勘探领域，更具体地，涉及一种多波动态匹配方法及系统。

背景技术

常规的多波地震勘探是指利用地下的反射纵波和转换横波进行油气藏勘探的一种技术手段。但是由于传播路径以及波场性质的差异，纵波和转换横波在走时、能量和相位都存在差异。因此，在多波勘探过程中，非常关键的一步既是对纵横波进行匹配。常规的匹配技术是估算地下的纵横波速比，再直接将转换波进行压缩，得到匹配结果。这种方法需要对纵横波速比进行估算，实现简单，但是误差比较大；另外一种常见的方法是利用井震标定，分别对纵波和转换波进行层位拾取，之后再根据对应层位进行匹配(付雷，田晓冬，2014.高精度三维多波层位匹配技术在大庆地区的应用.石油地球物理勘探，49(3)，572-577)。该方法相对更加准确，应用最为广泛，但是人工参与度非常高，主观性强，且层位之间无法保证匹配的准确性；更进一步，在解决走时差异之后，有些方法尝试解决纵横波能量和波形差异方面，例如通过统计纵横波的能量比例因子，再进行校正；通过整形滤波的方式，校正波形之间的差异性，以提高后续联合解释的精确度(一种纵横波匹配方法，专利号：CN200910237895)。对于后续的纵横波联合反演而言，这些匹配方法普遍存在的问题是：利用纵横波速比(不管是估算，还是利用人工拾取层位计算)进行匹配之后直接用于反演，而反演本身获得的纵横波速比信息，这两者之间通常存在前后不一致，从而影响了反演结果的可信度。因此，有必要开发一种多波动态匹配方法及系统。

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

发明内容

本发明提出了一种多波动态匹配方法及系统，其能够通过动态调整纵横波速比，解决了输入波速比值与反演得到的波速比值的不一致问题，实现了纵横波的自动匹配，提高了匹配的可信度。

根据本发明的一方面，提出了一种多波动态匹配方法。所述方法可以包括：步骤1：输入初始纵横波速比，进而获得初始纵横波速比变化率；步骤2：通过所述初始纵横波速比进行时间匹配，获得纵波和匹配后的转换波；步骤3：分别对所述纵波和所述转换波进行AVO反演，根据反演属性值计算基于反演的纵横波速比变化率；步骤4：将所述基于反演的纵横波速比变化率与所述初始纵横波速比变化率做差，获得纵横波速比变化率差值：步骤5：根据所述纵横波速比变化率差值，调整所述初始纵横波速比，重复步骤1-4，直到所述纵横波速比变化率差值最小，将此时的转换波作为最终转换波。

优选地，所述反演属性值包括纵波属性参数和转换波属性参数，其中：所述纵波属性参数的值通过对所述纵波进行AVO反演获得；所述转换波属性参数的值通过对所述转换波进行AVO反演获得。

优选地，通过公式(1)对所述纵波进行AVO反演，获得纵波属性参数的值：

R_pp(θ)≈A+Bsin²θ+C(tan²θ-sin²θ) (1)

其中，R_pp(θ)为纵波反射系数，θ为角度，A、B、C为纵波属性参数。

优选地，通过公式(2)对所述纵波进行AVO反演，获得转换波属性参数的值：