CN104834007B - 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 - Google Patents
地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 Download PDFInfo
- Publication number
- CN104834007B CN104834007B CN201510219252.8A CN201510219252A CN104834007B CN 104834007 B CN104834007 B CN 104834007B CN 201510219252 A CN201510219252 A CN 201510219252A CN 104834007 B CN104834007 B CN 104834007B
- Authority
- CN
- China
- Prior art keywords
- seismic
- disparity
- filling operation
- charges
- chiltern
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法,该方法包括:对研究区沿地震解释层位提取地震反演波阻抗平面图;获得碳酸盐岩缝洞型储层内充填物与充填程度信息;选择充填物相似但充填程度不同的样品点并定位其大地坐标;确定有效样品点并计算有效样品点之间的地震波阻抗差异数值;将有效样品点之间的砂质充填物充填量差异数值与地震波阻抗差异数值进行交汇分析得到函数关系式;计算每个有效样品点处的缝洞型储层的充填程度,根据函数关系式及该充填程度,得到地震波阻抗差异数值与砂质充填物充填程度的运算关系式;选取地震剖面上具有串珠状反射的地震反射区,运用运算关系式计算地震反射区缝洞型储层的充填程度。
Description
技术领域
本发明涉及地球物理技术领域,尤其涉及一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法。
背景技术
地震反演技术是充分利用钻井、测井、地质资料所提供的有关构造、岩性、层位等方面的丰富信息,由常规地震剖面来反演推导出地下地层的密度、波阻抗、速度、渗透率、孔隙度、压力、砂泥岩百分比等信息。波阻抗反演是指利用地震资料反演推导出地层波阻抗(或速度)的地震特殊处理解释技术。与地震模式识别预测油气、神经网络预测地层参数、振幅拟合预测储集层厚度等统计性方法相比,波阻抗反演处理解释技术具有十分明确的物理意义,是储集层岩性预测、油藏特征描述的确定性方法,在实际应用中取得了显著的地质效果,因此地震反演通常特指波阻抗反演。
目前碳酸盐岩缝洞型储层储量计算工作主要使用平均孔隙度这一参数,平均孔隙度法在计算缝洞型油气藏的过程中存在较大的误差,对实际的勘探开发生产工作造成了一定的不利影响。针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明提供了一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法,以至少解决目前碳酸盐岩缝洞型储层储量计算工作使用平均孔隙度法,存在较大误差的问题。
根据本发明的一个方面,提供了一种地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法,包括:步骤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所得的运算关系式对地震波阻抗数值进行再运算,将地震波阻抗数值转换为充填物数值。
通过上述方法,量化碳酸盐岩缝洞型储层内部充填程度在地震反演波阻抗数值上的表征,客观体现了储层内部真实孔隙度的数值,避免了使用平均孔隙度法对油气储量估算所导致的误差,从而能明确客观地反应碳酸盐岩缝洞型储层的真实地质储量。不同的主要充填物以及充填率的地震反演波阻抗数据是有差异的,因此,本发明分析不同充填程度的缝洞的地球物理响应并利用上述差异在平面上展开推出无井区缝洞型储层的充填程度。
在一个实施例中,步骤S101可以包括:根据测井数据和录井数据完成井上层位标定;根据井上层位标定结合地震反射特征完成地震解释;根据地震数据进行地震反演的二次运算;沿地震解释层位对二次运算后的反演波阻抗数据体进行切片,得到沿地震解释层位的地震反演波阻抗平面图。
在一个实施例中,地震波阻抗差异数值与砂质充填物充填量差异数值的函数关系式为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中,设定的所述充填量最小计算单元和所述充填程度最小计算单元均不小于研究对象要求精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219252.8A CN104834007B (zh) | 2015-05-04 | 2015-05-04 | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219252.8A CN104834007B (zh) | 2015-05-04 | 2015-05-04 | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104834007A CN104834007A (zh) | 2015-08-12 |
CN104834007B true CN104834007B (zh) | 2017-09-26 |
Family
ID=53812003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510219252.8A Active CN104834007B (zh) | 2015-05-04 | 2015-05-04 | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104834007B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572746A (zh) * | 2015-12-11 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种确定钻井位置的方法及装置 |
CN109143397B (zh) * | 2017-06-28 | 2020-05-19 | 中国石油化工股份有限公司 | 碳酸盐岩储层缝洞充填物识别方法及*** |
CN110244363B (zh) * | 2018-03-09 | 2020-10-23 | 中国石油化工股份有限公司 | 一种用于预测缝洞型储集层资源量的方法 |
CN109441541B (zh) * | 2018-11-06 | 2020-01-03 | 中国矿业大学 | 一种煤矿采空区充填体承载压缩率监测***及其监测方法 |
CN111596351B (zh) * | 2020-04-28 | 2023-04-25 | 中国石油天然气股份有限公司 | 碳酸盐岩输导体系定量评价方法、***、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
CN103064118A (zh) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | 一种利用声波测井定量洞穴充填程度的方法 |
CN103116192A (zh) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏储集体建模方法 |
CN103529475A (zh) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
-
2015
- 2015-05-04 CN CN201510219252.8A patent/CN104834007B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
CN103116192A (zh) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏储集体建模方法 |
CN103064118A (zh) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | 一种利用声波测井定量洞穴充填程度的方法 |
CN103529475A (zh) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
Non-Patent Citations (2)
Title |
---|
地震资料解释与盆地模拟在油气成藏研究中的联合运用;苏劲 等;《石油地球物理勘探》;20100228;第45卷(第1期);第72-78页 * |
碳酸盐岩洞穴充填物及其充填程度的测井判别方法;赵军 等;《中国岩溶》;20130630;第32卷(第2期);第225-230,237页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104834007A (zh) | 2015-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2895549C (en) | Fracturing and reactivated fracture volumes | |
US9910938B2 (en) | Shale gas production forecasting | |
CN104834007B (zh) | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 | |
EP3298437B1 (en) | Inversion for tectonic stress | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN104992468A (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
CN105653815A (zh) | 一种基于岩石物理模型理论的定量解释油藏流体分布方法 | |
CN109388817A (zh) | 一种储层裂缝三维建模方法 | |
CN104834008B (zh) | 计算碳酸盐岩缝洞型储层充填程度的方法 | |
CN105986815A (zh) | 一种用于识别页岩地层地质甜点的方法 | |
CN105388525B (zh) | 储层预测方法及装置 | |
CN106875471A (zh) | 煤系含或隔水层三维可视化建模方法 | |
US11249208B2 (en) | Geologic structural model generation | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN104914470B (zh) | 碳酸盐岩缝洞储层储量校正方法 | |
CN105929452B (zh) | 基于地震数据预测地下裂缝空间展布的方法及装置 | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN104063626A (zh) | 一种它源圈闭充满度的预测方法 | |
CN106249300B (zh) | 基于井震结合确定源岩toc含量和非均质性方法和装置 | |
CN104153768A (zh) | 一种评价花岗岩储层储集性能的方法 | |
CN103790579B (zh) | 随钻地质导向中确定钻头与地层界面距离的方法及装置 | |
Qian et al. | Prediction and modeling of petrophysical parameters of deep-buried, low permeability glutenite reservoirs in Yubei area, Turpan-Hami Basin, China | |
CN108665545A (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 |