[发明专利]一种三维双控下的叠前地质统计学反演方法

说明书

本发明涉及一种三维双控下的叠前地质统计学反演方法，包括以下步骤：1)目标工区岩石物理分析及井震标定，2)地震速度约束下低频模型的构建，3)基于部分叠加地震数据进行叠前同时反演，4)岩性流体概率分析；5)进行三维双控约束下的地质统计学反演。本发明在求取地质统计学反演先验信息岩性概率这个参数上，首先进行确定性反演，并基于贝叶斯判别原理，将反演得到的弹性数据转化为三维岩性概率体，并将其作为先验信息取代传统的零维岩性概率约束地质统计学反演，并对最终反演结果进行频率合并和降频处理，从而在提高地质统计学反演纵向分辨率的同时，依然能保持原始地震结构特征，降低反演不确定性，使其依然具有较好的预测能力。

技术领域

本发明涉及一种油气储层预测过程中的地震资料反演方法，特别是关于一种适用于纵波阻抗叠置、非均质性强储层的叠前地质统计学反演技术。

背景技术

相对于叠后反演只能得到纵波阻抗信息，难以适用于纵波阻抗叠置储层预测，叠前反演技术基于地震反射波振幅与不用入射角反射系数有关的理论，利用多个不同角度部分叠加地震数据反演各种弹性参数，如纵横波阻抗、密度、泊松比等，进而能够预测储层岩性、物性及流体。目前，叠前地质统计学反演技术已发展成为叠前高分辨率储层预测的有效方法，该技术可利用已知的地震信息、测井信息、地质信息，根据后验概率公式，求解空间中储层的分布特征。其综合测井垂向和地震横向分辨率高的优势，能进一步提高储层反演的空间分辨率，可有效识别薄储层，且可获得多个概率实现，提高储层预测的精度。

虽然地质统计学反演能得到高于地震分辨率的反演结果，但其反演过程中的纵横向变程、概率密度函数等关键参数，往往依靠经验，采用参数测试的方式选取，参数精度低、缺乏地质含义，反演结果不确定较大，预测性较低，难以精确预测横向非均质性强的储层展布，因此需要探索能进一步提升地质统计学反演预测性的方法。

发明内容

针对上述问题，本发明的目的是提供一种三维双控下的叠前地质统计学反演方法，在提高地质统计学反演纵向分辨率的同时，依然使其保持原始地震结构特征，具有较高预测性。

为实现上述目的，本发明采取以下技术方案，一种三维双控下的叠前地质统计学反演方法，其特征在于，包括以下步骤：

1)目标工区岩石物理分析及井震标定：

在测井资料处理及通过岩石物理建模进行横波预测的基础上，对测井进行岩石物理分析；

根据测井数据、部分叠加地震数据和地质数据进行井震标定，拾取多井子波，搭建可识别高精度地震地层格架；

2)地震速度约束下低频模型的构建：

采用地震速度约束下的多井插值法进行低频模型构建，具体过程包括：

①首先提取地震速度的平面属性，在井曲线位置处找到井曲线与属性之间的最佳关系；

②然后对远离井的井曲线值进行插值/外推，从而创建三维属性数据体；

③最后应用简单克里格法消除井点位置的预测误差，得到低频模型；

3)基于部分叠加地震数据进行叠前同时反演：

对叠加地震数据中的CRP道集进行动校正、去噪、拉平的优化处理，改善道集质量，并进行分角度叠加，得到5个部分叠加地震数据体，并采用Aki-Richard近似公式进行反演计算，得到弹性参数体，包括纵波阻抗、纵横波速度比和密度；

同时，采用频率域合并的方式，将步骤2)得到的低频模型频带范围内的值，合并到步骤3)得到的反演结果中，补充原始地震缺失的低频成分；

4)岩性流体概率分析：

对测井数据进行交汇分析，建立不同岩性对应步骤3)得到的反演弹性参数体的响应范围及定量的概率密度函数，利用贝叶斯推断将步骤3)得到的反演弹性参数体变为岩性概率体；