[发明专利]地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法

说明书

技术领域

本发明涉及地球物理技术领域，尤其涉及一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法。

背景技术

地震反演技术是充分利用钻井、测井、地质资料所提供的有关构造、岩性、层位等方面的丰富信息，由常规地震剖面来反演推导出地下地层的密度、波阻抗、速度、渗透率、孔隙度、压力、砂泥岩百分比等信息。波阻抗反演是指利用地震资料反演推导出地层波阻抗(或速度)的地震特殊处理解释技术。与地震模式识别预测油气、神经网络预测地层参数、振幅拟合预测储集层厚度等统计性方法相比，波阻抗反演处理解释技术具有十分明确的物理意义，是储集层岩性预测、油藏特征描述的确定性方法，在实际应用中取得了显著的地质效果，因此地震反演通常特指波阻抗反演。

目前碳酸盐岩缝洞型储层储量计算工作主要使用平均孔隙度这一参数，平均孔隙度法在计算缝洞型油气藏的过程中存在较大的误差，对实际的勘探开发生产工作造成了一定的不利影响。针对上述问题，目前尚未提出有效的解决方案。

发明内容

本发明提供了一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法，以至少解决目前碳酸盐岩缝洞型储层储量计算工作使用平均孔隙度法，存在较大误差的问题。

根据本发明的一个方面，提供了一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法，包括：步骤1，利用研究区的地震数据、测井数据及录井数据对所述研究区完成地震解释与地震反演，沿地震解释层位在反演结果上提取地震反演波阻抗平面图；步骤2，在对所述研究区进行岩心观察、录井解释报告分析的基础上进行小层对比，在等时格架内结合实际生产数据，获得碳酸盐岩缝洞型储层内部充填物与充填程度的信息；根据对取心井观察与数据统计，选择充填物相似但充填程度不同的样品点；并根据井位坐标、井斜数据与层位深度数据共同定位所述样品点的大地坐标；步骤3，将所述地震反演波阻抗平面图按数值格式导出，得到所述研究区的沿地震解释层位的地震波阻抗数值；确定与所述地震波阻抗数值中的大地坐标相同的样品点，作为有效样品点；并计算所述有效样品点之间的地震波阻抗差异数值；步骤4，根据研究对象精度要求，设定充填量最小计算单元，将所述有效样品点之间的砂质充填物充填量差异数值与所述地震波阻抗差异数值进行交汇分析，得到地震波阻抗差异数值与砂质充填物充填量差异数值的函数关系式；步骤5，根据研究对象精度要求，设定充填程度最小计算单元，根据所述地震波阻抗数值和地震剖面求取每个有效样品点处的每一缝洞型储层的总体积，计算每一缝洞型储层中砂质充填物充填量与总体积的比值，得到所述缝洞型储层的充填程度；根据所述地震波阻抗差异数值与所述砂质充填物充填量差异数值的函数关系式以及所述缝洞型储层的充填程度，得到所述地震波阻抗差异数值与砂质充填物充填程度的运算关系式；步骤6，选取所述地震剖面上具有串珠状反射的地震反射区，运用所述地震波阻抗差异数值与所述砂质充填物充填程度的运算关系式在选取的所述地震反射区内对沿所述地震解释层位导出的所述地震波阻抗数值进行运算，得到所述地震反射区的缝洞型储层的充填程度。

在一个实施例中，所述步骤1包括：根据所述测井数据和所述录井数据完成井上层位标定；根据井上层位标定结合地震反射特征完成地震解释；根据所述地震数据进行地震反演的二次运算；沿地震解释层位对二次运算后的反演波阻抗数据体进行切片，得到沿所述地震解释层位的所述地震反演波阻抗平面图。

在一个实施例中，所述地震波阻抗差异数值与砂质充填物充填量差异数值的函数关系式为Y＝a×X+b，其中，Y为砂质充填物充填量差异数值，X为地震波阻抗差异数值，a、b为交汇分析所得常量。

在一个实施例中，所述地震波阻抗差异数值与砂质充填物充填程度的运算关系式为Z＝(c×X+d)/V，其中，Z为砂质充填物充填程度，X为地震波阻抗差异数值，c、d为交汇分析所得运算常量，V为缝洞型储层的体积。

在一个实施例中，所述步骤2中，在构造单元中的不同构造部位以及不同构造单元中提取所述样品点。

在一个实施例中，所述步骤4和所述步骤5中，设定的所述充填量最小计算单元和所述充填程度最小计算单元均不小于研究对象要求精度。

通过本发明的计算碳酸盐岩缝洞型储层充填程度的方法，在地震反演过程中，量化碳酸盐岩缝洞型储层内部充填程度在地震反演波阻抗数值上的表征，客观体现了储层内部真实孔隙度的数值，避免了使用平均孔隙度法对油气储量估算所导致的误差，从而能明确客观地反应碳酸盐岩缝洞型储层的真实地质储量。

附图说明