CN104834008B - 计算碳酸盐岩缝洞型储层充填程度的方法 - Google Patents

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

Info

Publication number
CN104834008B
CN104834008B CN201510220776.9A CN201510220776A CN104834008B CN 104834008 B CN104834008 B CN 104834008B CN 201510220776 A CN201510220776 A CN 201510220776A CN 104834008 B CN104834008 B CN 104834008B
Authority
CN
China
Prior art keywords
seismic
charges
reservoir
disparity
filling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510220776.9A
Other languages
English (en)
Other versions
CN104834008A (zh
Inventor
邢翔
张义杰
王召明
段书府
杨海军
贾进华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510220776.9A priority Critical patent/CN104834008B/zh
Publication of CN104834008A publication Critical patent/CN104834008A/zh
Application granted granted Critical
Publication of CN104834008B publication Critical patent/CN104834008B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

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

Description

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

Claims (6)

1.一种计算碳酸盐岩缝洞型储层充填程度的方法,应用于地震振幅类属性勘探中,其特征在于,包括:
步骤1,利用研究区的地震数据、测井数据及录井数据对所述研究区完成地震解释,沿地震解释层位提取地震属性平面图;
步骤2,在对所述研究区进行岩心观察、录井解释报告分析的基础上进行小层对比,在等时格架内结合实际生产数据,获得碳酸盐岩缝洞型储层内部充填物与充填程度的信息;根据对取心井观察与数据统计,选择充填物相似但充填程度不同的样品点;并根据井位坐标、井斜数据与层位深度数据共同定位所述样品点的大地坐标;
步骤3,将所述地震属性平面图按数值格式导出,得到所述研究区的沿地震解释层位的地震属性数值;确定与所述地震属性数值中的大地坐标相同的样品点,作为有效样品点;并计算所述有效样品点之间的地震属性差异数值;
步骤4,根据研究对象精度要求,设定充填量最小计算单元,将所述有效样品点之间的砂质充填物充填量差异数值与所述地震属性差异数值进行交汇分析,得到地震属性差异数值与砂质充填物充填量差异数值的函数关系式;
步骤5,根据研究对象精度要求,设定充填程度最小计算单元,根据所述地震属性数值和地震剖面求取每个有效样品点处的每一缝洞型储层的总体积,计算每一缝洞型储层中砂质充填物充填量与总体积的比值,得到所述缝洞型储层的充填程度;根据所述地震属性差异数值与所述砂质充填物充填量差异数值的函数关系式以及所述缝洞型储层的充填程度,得到所述地震属性差异数值与砂质充填物充填程度的运算关系式;
步骤6,选取所述地震剖面上具有串珠状反射的地震反射区,运用所述地震属性差异数值与所述砂质充填物充填程度的运算关系式在选取的所述地震反射区内对沿所述地震解释层位导出的所述地震属性数值进行运算,得到所述地震反射区的缝洞型储层的充填程度。
2.根据权利要求1所述的方法,其特征在于,所述步骤1包括:
根据所述测井数据和所述录井数据完成井上层位标定;
根据井上层位标定结合地震反射特征完成地震解释;
根据所述地震数据进行属性体提取的二次运算;
沿地震解释层位对二次运算后的属性体进行切片,得到沿所述地震解释层位的所述地震属性平面图。
3.根据权利要求1所述的方法,其特征在于,所述地震属性差异数值与砂质充填物充填量差异数值的函数关系式为Y=a×X+b,其中,Y为砂质充填物充填量差异数值,X为地震属性差异数值,a、b为交汇分析所得常量。
4.根据权利要求1所述的方法,其特征在于,所述地震属性差异数值与砂质充填物充填程度的运算关系式为Z=(c×X+d)/V,其中,Z为砂质充填物充填程度,X为地震属性差异数值,c、d为交汇分析所得运算常量,V为缝洞型储层的体积。
5.根据权利要求1所述的方法,其特征在于,所述步骤2中,在构造单元中的不同构造部位以及不同构造单元中提取所述样品点。
6.根据权利要求1所述的方法,其特征在于,所述步骤4和所述步骤5中,设定的所述充填量最小计算单元和所述充填程度最小计算单元均不小于研究对象要求精度。
CN201510220776.9A 2015-05-04 2015-05-04 计算碳酸盐岩缝洞型储层充填程度的方法 Active CN104834008B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510220776.9A CN104834008B (zh) 2015-05-04 2015-05-04 计算碳酸盐岩缝洞型储层充填程度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510220776.9A CN104834008B (zh) 2015-05-04 2015-05-04 计算碳酸盐岩缝洞型储层充填程度的方法

Publications (2)

Publication Number Publication Date
CN104834008A CN104834008A (zh) 2015-08-12
CN104834008B true CN104834008B (zh) 2017-10-17

Family

ID=53812004

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510220776.9A Active CN104834008B (zh) 2015-05-04 2015-05-04 计算碳酸盐岩缝洞型储层充填程度的方法

Country Status (1)

Country Link
CN (1) CN104834008B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107558999B (zh) * 2017-07-03 2021-02-23 恒泰艾普集团股份有限公司 一种孤立缝洞体剩余油定量预测方法
CN109387870B (zh) * 2017-08-14 2020-06-23 中国石油化工股份有限公司 碳酸盐岩表层缝洞储集体刻画方法
CN108020863A (zh) * 2017-11-28 2018-05-11 河海大学 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法
CN109061763B (zh) * 2018-09-04 2020-08-28 中国地质大学(北京) 碳酸盐岩断溶体油藏洞穴测井综合评价方法
CN110941027B (zh) * 2018-09-21 2022-03-22 潜能恒信能源技术股份有限公司 一种计算碳酸盐岩溶蚀孔洞型地热能储量的方法和系统
CN110941028B (zh) * 2018-09-21 2022-03-25 潜能恒信能源技术股份有限公司 一种对碳酸盐岩溶蚀孔洞型地热能储层位置进行定位的方法和系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064118A (zh) * 2013-01-08 2013-04-24 西南石油大学 一种利用声波测井定量洞穴充填程度的方法
CN103487830B (zh) * 2013-09-17 2016-07-06 长江大学 一种地震振幅类属性勘探中补偿差异压实的方法

Also Published As

Publication number Publication date
CN104834008A (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
CN104834008B (zh) 计算碳酸盐岩缝洞型储层充填程度的方法
CN105510993B (zh) 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法
EP2846175B1 (en) Seismic survey analysis
EP3298437B1 (en) Inversion for tectonic stress
CN106951660A (zh) 一种海相碎屑岩水平井储层测井解释方法及装置
US10359529B2 (en) Singularity spectrum analysis of microseismic data
CN109388817A (zh) 一种储层裂缝三维建模方法
CN106780730A (zh) 三维地质模型的构建方法和装置
CN105986815A (zh) 一种用于识别页岩地层地质甜点的方法
CN106875471A (zh) 煤系含或隔水层三维可视化建模方法
CN104636980B (zh) 针对河道砂油藏类型油气汇集条件的地球物理表征方法
CN104155701B (zh) 一种利用叠前地震资料及井信息的多尺度裂缝预测方法
CN105676294B (zh) 一种碳酸盐岩缝洞型储集体定量描述方法
CN104678434A (zh) 一种预测储层裂缝发育参数的方法
CN104834007B (zh) 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法
CN103487830B (zh) 一种地震振幅类属性勘探中补偿差异压实的方法
CN103529473B (zh) 一种地震反演过程中差异压实的补偿方法
CN103278852B (zh) 利用地震数据体波形结构特征模型预测油气的方法
CN102221707B (zh) 含有逆掩断层的地质构造的地震属性提取方法
CN104820237B (zh) 计算地层的方差体的方法
CN112528106A (zh) 火山岩岩性识别方法
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
Wilson et al. Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site
Lyu et al. Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China
CN103217718B (zh) 一种补充井下缺失地层的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant