计算碳酸盐岩缝洞型储层充填程度的方法技术

本发明专利技术公开了一种计算碳酸盐岩缝洞型储层充填程度的方法，应用于地震振幅类属性勘探中，包括：对研究区沿地震解释层位提取地震属性平面图；获得碳酸盐岩缝洞型储层内充填物与充填程度信息；选择充填物相似但充填程度不同的样品点并定位其大地坐标；确定有效样品点并计算有效样品点之间的地震属性差异数值；将有效样品点之间的砂质充填物充填量差异数值与地震属性差异数值进行交汇分析得到函数关系式；计算每个有效样品点处的缝洞型储层的充填程度，根据函数关系式及该充填程度，得到地震属性差异数值与砂质充填物充填程度的运算关系式；选取地震剖面上具有串珠状反射的地震反射区，运用运算关系式计算地震反射区缝洞型储层的充填程度。

【技术实现步骤摘要】

本专利技术涉及地球物理
，尤其涉及一种地震振幅类属性勘探中。
技术介绍
地震属性是指那些由叠前或叠后地震数据经过数学变换而导出的有关地震波的几何形态、运动学特征和统计学特征的特殊测量值。根据不同的算法、不同的属性具有不同的地质意义，有些属性值可能擅长于揭示不易探测到的岩性变化，而有些属性则可以直接用于烃类检测。地震属性的研宄和使用始于20世纪70年代，最早主要是以振幅为基础的瞬时属性，用来直接指示油气与地下不同岩层的展布情况。70年代早期，Nigel Anstey发现了含气砂岩波阻抗的异常变化，使用了反射波振幅变化特征“亮点”与“暗点”对含气砂岩储集体进行预测。80年代，地震属性的数量迅速增加，多元属性分析技术首次被使用。80年代晚期，多维属性如倾角和方位角的初步发展，导致了三维连续属性在90年代的广泛应用。宄其原因，它的成功得益于清楚和明确的地质含义。90年代以来，基于精细储层描述以及三维地震数据体解释的需要，地震属性分析技术更是急剧发展。现今，地震属性技术已广泛应用于储层预测、油气藏动态监测、油气藏特征描述等领域，并取得了很好的效果。综合地震属性的数学、物理及地质意义，也可以将地震属性分类两类:第一类是基本属性，包括反射时间、振幅及相位，它们各自具有独立的物理意义和地质意义；第二类是延伸属性，如对多种属性进行处理(叠合、求差、相关、求主成分等)得到的复合属性以及根据单一地震属性再提取出的新的级联属性。其中振幅类属性主要是用于体现地下不同岩性地层的分布特征与展布规律。然而传统的勘探技术不易探测到岩性的变化，且不能直接用于烃类检测，对于寻找剩余油、提高采收率成为油田开发工作的重点及难点，也是稳产的关键。目前国内外已知的碳酸盐岩缝洞型储层中计算孔隙度的方法，总体上还是平均孔隙度法，这种方法没有考虑到不同缝洞型储层中的充填程度的差异性，对实际地下资源量计算造成较大的误差，进而影响到了实际的油田勘探开发工作部署，造成经济损失。
技术实现思路
本专利技术提供了一种，以至少解决目前碳酸盐岩缝洞型储层储量计算工作使用平均孔隙度法，存在较大误差的问题。根据本专利技术的一个方面，提供了一种，应用于地震振幅类属性勘探中，该方法包括:步骤1，利用研宄区的地震数据、测井数据及录井数据对所述研宄区完成地震解释，沿地震解释层位提取地震属性平面图；步骤2，在对所述研宄区进行岩心观察、录井解释报告分析的基础上进行小层对比，在等时格架内结合实际生产数据，获得碳酸盐岩缝洞型储层内部充填物与充填程度的信息；根据对取心井观察与数据统计，选择充填物相似但充填程度不同的样品点；并根据井位坐标、井斜数据与层位深度数据共同定位所述样品点的大地坐标；步骤3，将所述地震属性平面图按数值格式导出，得到所述研宄区的沿地震解释层位的地震属性数值；确定与所述地震属性数值中的大地坐标相同的样品点，作为有效样品点；并计算所述有效样品点之间的地震属性差异数值；步骤4，根据研宄对象精度要求，设定充填量最小计算单元，将所述有效样品点之间的砂质充填物充填量差异数值与所述地震属性差异数值进行交汇分析，得到地震属性差异数值与砂质充填物充填量差异数值的函数关系式；步骤5，根据研宄对象精度要求，设定充填程度最小计算单元，根据所述地震属性数值和地震剖面求取每个有效样品点处的每一缝洞型储层的总体积，计算每一缝洞型储层中砂质充填物充填量与总体积的比值，得到所述缝洞型储层的充填程度；根据所述地震属性差异数值与所述砂质充填物充填量差异数值的函数关系式以及所述缝洞型储层的充填程度，得到所述地震属性差异数值与砂质充填物充填程度的运算关系式；步骤6，选取所述地震剖面上具有串珠状反射的地震反射区，运用所述地震属性差异数值与所述砂质充填物充填程度的运算关系式在选取的所述地震反射区内对沿所述地震解释层位导出的所述地震属性数值进行二次运算，得到所述地震反射区的缝洞型储层的充填程度。在一个实施例中，所述步骤I包括:根据所述测井数据和所述录井数据完成井上层位标定；根据井上层位标定结合地震反射特征完成地震解释；根据所述地震数据进行属性体提取的运算；沿地震解释层位对二次运算后的属性数据体进行切片，得到沿所述地震解释层位的所述地震属性平面图。在一个实施例中，所述地震属性差异数值与砂质充填物充填量差异数值的函数关系式为Y = aXX+b，其中，Y为砂质充填物充填量差异数值，X为地震属性差异数值，a、b为交汇分析所得常量。在一个实施例中，所述地震属性差异数值与砂质充填物充填程度的运算关系式为Z = (cXX+d)/V，其中，Z为砂质充填物充填程度，X为地震属性差异数值，C、d为交汇分析所得运算常量，V为缝洞型储层的体积。在一个实施例中，所述步骤2中，在构造单元中的不同构造部位以及不同构造单元中提取所述样品点。在一个实施例中，所述步骤4和所述步骤5中，设定的所述充填量最小计算单元和所述充填程度最小计算单元均不小于研宄对象要求精度。通过本专利技术的，在地震振幅类属性勘探中，量化碳酸盐岩缝洞型储层内部充填程度在地震振幅类属性数值上的表征，客观体现了储层内部真实孔隙度的数值，避免了使用平均孔隙度法对油气储量估算所导致的误差，从而能明确客观地反应碳酸盐岩缝洞型储层的真实地质储量。【附图说明】此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的限定。在附图中:图1是本专利技术实施例的地震振幅类属性勘探中的流程图。【具体实施方式】下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术的保护范围。本专利技术实施例提供了一种地震振幅类属性勘探中计算碳酸盐岩缝洞型储层充填程度(充填程度也称为充填率)的方法，图1是本专利技术实施例的地震振幅类属性勘探中的流程图。如图1所示，该方法包括:步骤S101，利用研宄区的地震数据、测井数据及录井数据对研宄区完成地震解释，沿地震解释层位提取地震属性平面图。其中，地震数据为segy格式的数据。沿地震解释层位提取地震属性平面图，要求地震层位以及断层解释，然后合成记录制作以及时深转换公式。步骤S102，在对研宄区进行岩心观察、录井解释报告分析的基础上进行小层对比，在等时格架内结合实际生产数据，获得碳酸盐岩缝洞型储层内部充填物与充填程度的信息；根据对取心井观察与数据统计，选择充填物相似但充填程度不同的样品点；并根据井位坐标、井斜数据与层位深度数据共同定位样品点的大地坐标。同时，也可以统计样品点数据，样品点数据包括主要充填物与充填程度。步骤S103，将地震属性平面图按数值格式导出，得到研宄区的沿地震解释层位的地震属性数值；确定与地震属性数值中的大地坐标相同的样品点，作为有效样品点；并计算有效样品点之间的地震属性差异数值。本步骤中，将地震属性数值中的大地坐标与样品点的大地坐标进行对应，要求两个坐标尽可能接近甚至重合，要求平面上距离最短，需考虑样品点坐标中方位角和偏移距能引起误差的各种因素。步骤S104，根据研宄对象精度要求，设定充填量最小计算单元，本文档来自技高网...

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

【技术特征摘要】

【专利技术属性】
技术研发人员：邢翔，张义杰，王召明，段书府，杨海军，贾进华，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11