依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法技术

依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法是一种油气地震勘探数据处理与解释技术，它实现了利用交会分析来直接指导依赖频率的AVO频散参数反演中敏感的参考频率的选择方法。首先需要在地震有效频带范围内，以全部频率对叠前角道集进行纵横波频散属性参数反演，充分利用已知的测井解释数据进行约束，将含流体储层段的全部频度的纵横波频散属性投到频率—频散属性交会平面图上，最后通过交会分析，确定含流体储层的纵横波频散属性参数极大值对应的频率以及不同流体储层的纵横波频散属性差异最大时对应的频率或频段，此即为敏感的特征频率，以此频率作为敏感的最优参考频率对叠前角道集进行依赖频率的AVO纵横波频散属性参数反演，则可最大限度突出含流体储层的地震频散异常，且使不同流体储层的频散差异达到最大，从而实现不同流体的区分和识别。

【技术实现步骤摘要】
依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法
本专利技术涉及油气地震勘探数据处理与解释领域，是一种在依赖频率的AVO频散参数反演中，通过交会分析方法来确定反演的参考频率，从而以此频率作为敏感的参考频率，对地震叠前角道集进行频散参数反演，使储层在含不同流体时的地震纵横波频散差异最大，最大限度地区分不同流体储层，突出油气储层的地震频散异常的技术。
技术介绍
AVO(AmplitudeVersusOffset)是利用Zoeppritz方程描述地震波的反射振幅随入射角(或炮检距)变化的规律，用于提取储层岩性、物理参数及孔隙流体的信息。传统的AVO分析基于Zoeppritz方程，讨论反射系数与界面上下层纵横波速度及密度之间的关系，通过岩层弹性参数变化对地震振幅的影响，实现有利的油气储层预测。然而，在不同角度的地震道集中，含流体层的地震反射不仅在振幅上存在差异(即AVO现象)，而且，地震反射中不同频率分量振幅的变化特征也存在差异，即对于同一类型的AVO，当流体特征(如饱和度等参数)不同时，不同频率的地震反射振幅的大小及变化趋势随入射角的变化规律是不同的，而这种差别与储层中流体引起的频散和衰减的共同作用关系密切(Chapman&Odebeatu，2005；2006)。然而，常规的AVO技术没有考虑地震反射中不同频率分量的特征参数随入射角变化的差异，忽略了频率的因素。事实上，地震波的反射系数与频率之间的关系非常密切，地震波经过油气储层后表现出来的强吸收衰减、速度频散等异常现象，反映为反射系数的频率依赖特性。针对上述问题，将常规AVO分析技术进行拓展，发展了新的依赖频率的AVO分析技术，基于地震反射系数会随入射角度(或炮检距)和频率共同变化的特性，利用依赖频率的AVO纵横波频散参数反演，可以检测与油气饱和度有关的地震频散与衰减异常，进一步提高油气检测的精度。由于地震反射系数存在频率依赖特性，当储层含不同流体时，地震反射信号的不同频率分量对于不同流体响应的频散特性存在差异，因此，在依赖频率的AVO纵横波频散参数反演中，如何确定参考频率十分重要，它关系到利用纵横波频散异常识别和区分不同流体、实现油气检测的敏感程度和可靠性。然而，目前在利用地震叠前角道集进行依赖频率的AVO频散反演过程中，选择参考频率的方法通常为“人工观察法”：首先对工区的叠后地震资料进瞬时频谱分解，得到若干共频率剖面，人工观察共频率剖面中已知油气储层的瞬时谱能量随不同频率变化的特征，一般选择接近地震资料主频、油气储层瞬时谱能量相对较强时对应的频率作为参考频率，这种确定方法容易存在较强的人为主观性；另外，叠前角道集中反映的不同流体的频散异常信息在叠后地震资料中易畸变或被削弱，故利用叠后资料也不利于准确确定叠前角道集地震数据中对油气更敏感的参考频率，限制了依赖频率的AVO频散参数反演的精度和可靠性。
技术实现思路
本专利技术是要提供一种依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，它可以直接指导依赖频率的AVO频散参数反演中参考频率的选择，尤其是需要利用频散属性区分储层的流体性质时，本专利技术的方法可以搜索并确定不同流体储层的频散属性差异最大时所对应的频率或频段，这样，利用该频率或频段作为敏感的参考频率进行反演，则可尽量突出不同流体的频散差异，直接指导储层流体类型、饱和度等信息的识别。本专利技术的依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，在确定敏感的参考频率时，充分利用了地震有效频带范围内的全部频率，以及测井解释数据进行约束，对叠前角道集进行纵横波频散参数反演，获得不同流体储层的纵横波频散属性异常最显著、差异最大时对应的频率或频段。本专利技术的依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，是在频率—频散属性平面进行交会分析，能够快速可靠地拾取对应纵横波频散属性异常最为敏感的参考频率。本专利技术的依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，具有如下优越性：(1)在地震有效频带范围内利用地震信号的全部频率成分进行最优参考频率的搜索，不同于人工观察法，有利于精细地反映不同流体储层纵横波频散异常最显著的频率或频带。(2)可充分利用已知的测井解释数据进行约束，在频率—频散属性平面直接对含流体储层顶底地震反射的纵横波频散属性进行交会分析，快速准确地拾取最敏感的参考频率进行依赖频率的AVO频散参数反演，可有效提高依赖频率的AVO频散参数反演的可靠性以及油气检测的精度，计算效率高，算法简洁。本专利技术的具体实现原理如下：在Smith和Gidlow(1987)的AVO二项近似表达式的基础上，由于地震反射存在依赖频率的速度频散，因此引入频率维，在入射角度—频率域建立依赖频率的AVO反射系数分布R(θ,f)，使AVO反射系数分布构成了入射角度和频率的函数，形式如下：令其中，θ是入射角度，f为频率，vp和vs分别为反射界面上下层的平均纵横波速度，△vp和△vs分别为反射界面上下层的纵波速度差和横波速度差。对上述公式中的和关于频率fi进行泰勒展开，并舍去高阶项，可以得到下式：其中分别为纵横波频散属性参数，它们可通过基于上式的最小二乘反演，即对叠前地震角道集进行依赖频率的AVO纵横波频散属性参数反演获得。在频率—频散属性平面，对含流体储层顶底地震反射的纵横波频散属性进行交会分析，快速准确地拾取最敏感的参考频率的过程为：(1)在资料的有效频带范围内，以不同的频率作为参考频率，利用工区所有已知井的井旁叠前角道集，以及各井的测井解释已知的含流体类型，分别计算目的储层顶底(基于井上层位标定确定的目的层段及其顶底位置)的纵横波频散属性，并分别投到频率—频散属性平面上；(2)通过上述步骤，所有已知井目的储层段的纵、横波频散属性均投到了频率—纵波频散属性平面、频率—横波频散属性平面，形成了交会图，据此，即可直接拾取目的储层段纵、横波频散属性趋于最大值时对应的频率或频段，对于工区已知井的目的层段存在不同流体类型或饱和度的情况，则在交会图中分别用不同的颜色或符号表示，这样，可以拾取不同流体或饱和度的储层的纵横波频散属性差异最大时对应的频率或频段；(3)以上述步骤确定的敏感频率或频带作为敏感的最优参考频率，对全区的三维叠前地震角道集执行依赖频率的AVO纵横波频散属性参数反演，从而可以得到最能突出流体异常、区分不同流体类型或饱和度的纵横波频散属性数据体，指导储层流体识别。附图说明图1是设计的5层地质模型(模型中第三层为含两种不同流体的砂岩，)，第一层为泥岩，厚度为300米，纵波速度为1810米/秒，横波速度为1038米/秒，密度为2.16克/立方厘米；第二层为泥岩，厚度为254米，纵波速度为2037米/秒，横波速度为1237米/秒，密度为2.139克/立方厘米；第三层为含流体砂岩，其中两侧为含水砂岩，中间为含气砂岩，厚度为226米，纵波速度为2771米/秒，横波速度为1499米/秒，密度为2.08克/立方厘米；第四层为泥岩，厚度为220米，纵波速度为2037米/秒，横波速度为1237米/秒，密度为2.139克/立方厘米；第五层为泥岩，厚度为200米，纵波速度为2310米/秒，横波速度为1438米/秒，密度为2.20克/立方厘米。图2是与图1的地质模型对应本文档来自技高网...

【技术保护点】
一种依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，其特征在于采用以下具体步骤：⑴导入工区所有已知井的井旁叠前角道集，以及各井测井解释的目的层段及其位置，目的层段的流体类型数据；⑵在地震有效频带范围内，利用井旁叠前角道集的全部频率成分在目的储层段执行依赖频率的AVO纵横波频散参数反演；⑶将目的层段的每个频率成分对应的纵横波频散属性参数，根据测井解释的流体类型用不同的符号或颜色，投到频率—频散属性交会平面图上；⑷对频率—频散属性交会平面图进行交会分析，搜索并拾取纵横波频散属性极大值对应的频率以及不同流体储层的纵横波频散属性差异最大时对应的频率或频段，此即为最敏感的特征频率；⑸以步骤⑷的最敏感的特征频率作为最优参考频率，对工区的全部三维叠前角道集进行依赖频率的AVO频散参数反演，获取最能突出储层流体异常，进而区分不同流体类型的频散属性参数，指导储层油气识别。

【技术特征摘要】
1.一种依赖频率的AVO频散参数反演中敏感参考频率的交会确定方法，其特征在于采用以下具体步骤：导入工区所有已知井的井旁叠前角道集，以及各井测井解释的目的层段及其位置，目的层段的流体类型数据；在地震有效频带范围内，利用井旁叠前角道集的全部频率成分在目的储层段执行依赖频率的AVO纵横波频散参数反演；将目的层段的每个频率成分对应的纵横波频散属性参数，根据测井解释的流体类型用不同的符号或颜色，投到频率—频散属性交会平面图上；对频率—频散属性交会平面图进行交会分析，搜索并拾取纵横波频散属性极大值对应的频率以及不同流体储层的纵横波频散属性差异最大时对应的频率或频段，此即为最敏感的特征频率；以步骤的最敏感的特征频率作为最优参考频率，对工区的全部三维叠前角道集进行依赖频率...

【专利技术属性】
技术研发人员：陈学华，罗鑫，许迪，陈辉，齐迎凯，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川;51