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

摘要

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

主权项

1.一种三维双控下的叠前地质统计学反演方法，其特征在于，包括以下步骤：1)目标工区岩石物理分析及井震标定：在测井资料处理及通过岩石物理建模进行横波预测的基础上，对测井进行岩石物理分析；根据测井数据、部分叠加地震数据和地质数据进行井震标定，拾取多井子波，搭建可识别高精度地震地层格架；2)地震速度约束下低频模型的构建：采用地震速度约束下的多井插值法进行低频模型构建，具体过程包括：①首先提取地震速度的平面属性，在井曲线位置处找到井曲线与属性之间的最佳关系；②然后对远离井的井曲线值进行插值/外推，从而创建三维属性数据体；③最后应用简单克里格法消除井点位置的预测误差，得到低频模型；3)基于部分叠加地震数据进行叠前同时反演：对叠加地震数据中的CRP道集进行动校正、去噪、拉平的优化处理，改善道集质量，并进行分角度叠加，得到5个部分叠加地震数据体，并采用Aki‑Richard近似公式进行反演计算，得到弹性参数体，包括纵波阻抗、纵横波速度比和密度；同时，采用频率域合并的方式，将步骤2)得到的低频模型频带范围内的值，合并到步骤3)得到的反演结果中，补充原始地震缺失的低频成分；4)岩性流体概率分析：对测井数据进行交汇分析，建立不同岩性对应步骤3)得到的反演弹性参数体的响应范围及定量的概率密度函数，利用贝叶斯推断将步骤3)得到的反演弹性参数体变为岩性概率体；5)进行三维双控约束下的地质统计学反演，具体过程如下：①进行统计学建模：以精细解释地震层位搭建的精细地层格架、综合地质和测井认识定纵横向变程、以及多变量概率密度函数、岩性概率体作为先验信息进行统计学建模；在这个过程中，利用步骤4)中所得的三维岩性概率体代替常规的零维岩性概率体，约束得到先验概率密度函数；②然后，根据地震数据和测井信息，综合步骤①中的先验概率密度函数得到储层的后验概率分布函数；③基于马尔科夫链‑蒙特卡罗算法对后验概率密度函数进行采样，得到反演弹性参数体；④对步骤3)中的反演结果进行频率合并和降频处理；6)反演结果对比及盲井检验。