CN115903026A - 复合砂体构型解析方法、设备及介质 - Google Patents

复合砂体构型解析方法、设备及介质 Download PDF

Info

Publication number
CN115903026A
CN115903026A CN202310024380.1A CN202310024380A CN115903026A CN 115903026 A CN115903026 A CN 115903026A CN 202310024380 A CN202310024380 A CN 202310024380A CN 115903026 A CN115903026 A CN 115903026A
Authority
CN
China
Prior art keywords
curve
configuration
composite sand
sand body
mode
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.)
Granted
Application number
CN202310024380.1A
Other languages
English (en)
Other versions
CN115903026B (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.)
Sanya Offshore Oil And Gas Research Institute Of Northeast Petroleum University
Original Assignee
Sanya Offshore Oil And Gas Research Institute Of Northeast Petroleum University
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 Sanya Offshore Oil And Gas Research Institute Of Northeast Petroleum University filed Critical Sanya Offshore Oil And Gas Research Institute Of Northeast Petroleum University
Priority to CN202310024380.1A priority Critical patent/CN115903026B/zh
Publication of CN115903026A publication Critical patent/CN115903026A/zh
Application granted granted Critical
Publication of CN115903026B publication Critical patent/CN115903026B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请公开了一种复合砂体构型解析方法、设备及介质,涉及石油开发技术领域,该方法包括:利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线;利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体;基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。能够获取更加具体的复合砂体构型。

Description

复合砂体构型解析方法、设备及介质
技术领域
本发明涉及石油开发技术领域,特别涉及复合砂体构型解析方法、设备及介质。
背景技术
以潮汐为主要营力的潮控河口湾沉积体系是最主要的油气储集层类型之一,其主要发育潮汐河道、潮汐沙坝、沙坪、泥坪等沉积微相。受潮汐和河流共同作用,该类沉积单元内砂体演化复杂,垂向及侧向切叠严重。因此单纯依靠井点信息或地震资料难以对潮控河口湾复合砂体进行准确解析。在油田开发的中后期,石油产量减少,大量的石油以剩余油的形式存在于储层之中,而储层砂体的叠置关系是制约剩余油分布的重要因素。复合砂体构型是指空间上多个单一成因砂体(构型单元)的叠置关系,包括平面的及垂向上的复合。不同的叠置关系控制着不同构型单元间的流通性,如两个构型单元构型特征为空间上相对独立,则每个构型单元内的流体也相对独立,在油田的注采开发过程中应将两个构型单元采用不同的注采方案,如果两个构型单元构型特征为空间上叠合在一起,则需要考虑其叠合界面对流体的阻隔能力。如过界面阻隔能力强,将流体完全阻隔,则参照独立构型单元开发,如果界面不能完全阻隔,两个构型单元间流体具有流通性,则将两个构型单元视为同一储集单元开采。因此如何将岩心资料、稀井网资料、三维地震资料有机结合进行潮控河口湾复合砂体解析是油田地质研究的重点研究内容,对于油田持续稳产,提高采收率具有重要意义。
目前油田对于储层构型研究程度非常高,领先于国际开发技术,老油田由于多年开采,在开发区的井控程度高,井网密度大,井距可达到200m,小于复合砂体规模,应用钻井、测井资料,通过多井组合,井网控制就可以在空间上解析地下储集层的复合砂体叠置关系。但对于井网密度小、井距大的研究区,密井网条件下复合砂体构型解析的方法无法使用。潮控河口湾类型的沉积环境中,复合砂体的构型研究成熟程度相较河流相复合砂体构型研究还处在最初级的阶段,构型单元类型不明确。
综上可见,如何获取更加具体的复合砂体构型是本领域有待解决的问题。
发明内容
有鉴于此,本发明的目的在于提供一种复合砂体构型解析方法、设备及介质,能够获取更加具体的复合砂体构型。其具体方案如下:
第一方面,本申请公开了一种复合砂体构型解析方法,包括:
利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线;
利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体;
基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。
可选的,所述利用测井曲线建立构型要素类别曲线,包括:
利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线;其中,所述构型元素类型包括潮汐河道、潮汐沙坝、潮汐泥坪、潮汐沙坪中任意一种或几种。
可选的,所述利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线,包括:
利用岩心标定测井曲线,并进行井震标定,以便利用预设地震数据和预设测井数据确定目标复合砂体的顶底界面;
利用所述岩心确定所述目标复合砂体的构型元素类型,并基于所述目标复合砂体的顶底界面逐井识别所述目标复合砂体的内部的隔层顶底;
基于所述隔层顶底,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线。
可选的,所述利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线,包括:
利用傅里叶变换方法对所述测井曲线进行傅里叶变换,得到所述测井曲线的第一频谱,并在频域下利用所述第一频谱融合所述构型要素类别曲线和所述测井曲线,以得到第一融合后曲线,利用所述傅里叶变换方法对所述第一融合后曲线进行傅里叶反变换,以得到岩性模式特征曲线;
利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,得到所述孔渗曲线的第二频谱,并在频域下利用所述第二频谱融合所述构型要素类别曲线和所述孔渗曲线,以得到第二融合后曲线,利用所述傅里叶变换方法对所述第二融合后曲线进行傅里叶反变换,以得到孔渗模式特征曲线。
可选的,所述利用傅里叶变换方法对所述测井曲线进行傅里叶变换,包括:
对所述测井曲线进行归一化处理,以得到归一化后的测井曲线,并利用傅里叶变换方法对所述归一化后的测井曲线进行傅里叶变换;
相应的,所述利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,包括:
对孔渗曲线进行归一化处理,以得到归一化后的孔渗曲线,并利用傅里叶变换方法对所述归一化后的孔渗曲线进行傅里叶变换。
可选的,所述利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,包括:
利用预设地震数据和所述岩性模式特征曲线确定第一波形指示的低频特征和高频特征,以便利用所述第一波形指示进行模拟计算,得到目标复合砂体的岩性模式反演体;
相应的,所述基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体,包括:
基于所述岩性模式反演体确定第二波形指示的低频特征,并利用所述孔渗模式特征曲线确定所述第二波形指示的高频特征,以便利用所述第二波形指示进行模拟计算,得到所述目标复合砂体的孔渗模式反演体。
可选的,所述基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征,包括:
确定过井剖面,调整与所述岩性模式反演体对应的色标显示范围,以便确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征;
迁移所述构型界面信息至所述孔渗模式反演体,以便获取所述目标复合砂体的物性特征。
第二方面,本申请公开了一种电子设备,包括:
存储器,用于保存计算机程序;
处理器,用于执行所述计算机程序,以实现前述公开的复合砂体构型解析方法的步骤。
第三方面,本申请公开了一种计算机可读存储介质,用于存储计算机程序;其中,所述计算机程序被处理器执行时实现前述公开的复合砂体构型解析方法的步骤。
本申请有益效果为:利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线;利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体;基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。由此可见,本申请利用述构型要素类别曲线和测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线,并基于岩性模式特征曲线和孔渗模式特征曲线得到岩性模式反演体和孔渗模式反演体,进而可以获取在潮控河口湾类型的沉积环境中,更加具体的目标复合砂体构型信息,即目标复合砂体的构型界面信息、构型特征以及物性特征,以便后续基于目标复合砂体的构型特征以及目标复合砂体的物性特征,可以更加容易的进行石油开发相关研究。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本申请公开的一种复合砂体构型解析方法流程图;
图2为本申请公开的一种具体的复合砂体构型解析方法流程图;
图3为本申请公开的另一种具体的复合砂体构型解析方法流程图;
图4为本申请公开的一种具体的构型界面解释结果示意图;
图5为本申请公开的一种具体的砂体垂向叠置类型模式关系示意图;
图6为本申请公开的一种电子设备结构图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前油田对于储层构型研究程度非常高,领先于国际开发技术,老油田由于多年开采,在开发区的井控程度高,井网密度大,井距可达到200m,小于复合砂体规模,应用钻井、测井资料,通过多井组合,井网控制就可以在空间上解析地下储集层的复合砂体叠置关系。但对于井网密度小、井距大的研究区,密井网条件下复合砂体构型解析的方法无法使用。潮控河口湾类型的沉积环境中,复合砂体的构型研究成熟程度相较河流相复合砂体构型研究还处在最初级的阶段,构型单元类型不明确。
为此本申请相应的提供了一种复合砂体构型解析方案,能够获取更加具体的复合砂体构型。
参见图1所示,本申请实施例公开了一种复合砂体构型解析方法,包括:
步骤S11:利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线。
在井点处识别目标复合砂体顶底,并在复合砂体内部识别隔层顶底,建立单井复合砂体基础框架,并创建构型要素类别曲线。可以理解的是,输入测井曲线,则输出目标复合砂体的顶部与目标复合砂体的底部之间的深度值、目标复合砂体内部隔层深度值、构型要素类别曲线。
在单井复合砂体基础框架控制下,使用傅里叶变换方法融合构型要素类别曲线和测井曲线,构建复合砂体特征曲线,即分别为岩性模式特征曲线和孔渗模式特征曲线。可以理解的是,输入测井曲线和构型要素曲线,则输出岩性模式特征曲线,孔渗模式特征曲线。
步骤S12:利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体。
本实施例中,所述利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,具体包括:利用预设地震数据和所述岩性模式特征曲线确定第一波形指示的低频特征和高频特征,以便利用所述第一波形指示进行模拟计算,得到目标复合砂体的岩性模式反演体。使用复合砂体顶底建立构造模型,预设地震数据作为第一波形指示的低频特征,岩性模式特征曲线作为第一波形指示的高频特征,使用第一波形指示进行模拟计算,以得到复合砂体的岩性模式反演体。岩性模式反演体是岩性体,体现空间上的岩性变化。
步骤S13:基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。
本实施例中,所述基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体,包括:基于所述岩性模式反演体确定第二波形指示的低频特征,并利用所述孔渗模式特征曲线确定所述第二波形指示的高频特征,以便利用所述第二波形指示进行模拟计算,得到所述目标复合砂体的孔渗模式反演体。使用复合砂体顶底建立构造模型,岩性模式反演体作为第二波形指示的低频特征,孔渗特征曲线作为第二波形指示的,使用第二波形指示进行模拟计算,以得到复合砂体的孔渗模式反演体。孔渗模式反演体是体现孔隙度和渗透率三维空间变化的反演体。
本实施例中,所述基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征,包括:确定过井剖面,调整与所述岩性模式反演体对应的色标显示范围,以便确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征;迁移所述构型界面信息至所述孔渗模式反演体,以便获取所述目标复合砂体的物性特征。首先在岩性模式反演体上解释构型界面,具体方法是选择过井剖面,调整反演体色标显示范围,凸显构型界面的特征,逐个过井剖面识别并追踪砂体的顶底界面,反演体中解释构型界面要确保每个解释剖面都有井或者过井解释结果作为参考。解释完成后,迁移岩性模式反演体上解释的构型界面到孔渗模式特征反演体中;然后,获取垂直流向的解释剖面,分析各构型单元间的垂向叠置关系,得到研究目标的砂体垂向叠置类型模式,制作延构型界面提取的岩性模式反演平面属性和孔渗模式反演平面属性,用以反应单砂体间的渗流能力。
可见,本申请利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线;利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体;基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。由此可见,本申请利用述构型要素类别曲线和测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线,并基于岩性模式特征曲线和孔渗模式特征曲线得到岩性模式反演体和孔渗模式反演体,进而可以获取在潮控河口湾类型的沉积环境中,更加具体的目标复合砂体构型信息,即目标复合砂体的构型界面信息、构型特征以及物性特征,以便后续基于目标复合砂体的构型特征以及目标复合砂体的物性特征,可以更加容易的进行石油开发相关研究。
参见图2所示,本申请实施例公开了一种具体的复合砂体构型解析方法,包括:
步骤S21:利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线;其中,所述构型元素类型包括潮汐河道、潮汐沙坝、潮汐泥坪、潮汐沙坪中任意一种或几种。
本实施例中,所述利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线,包括:利用岩心标定测井曲线,并进行井震标定,以便利用预设地震数据和预设测井数据确定目标复合砂体的顶底界面;利用所述岩心确定所述目标复合砂体的构型元素类型,并基于所述目标复合砂体的顶底界面逐井识别所述目标复合砂体的内部的隔层顶底;基于所述隔层顶底,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线。
先使用有岩心的井,利用岩心标定测井曲线,并完成井震标定,确保复合砂体顶底地震层位和测井分层严格对应,完成研究区内高精度层序地层格架建立,在地震数据和测井数据上确定目标复合砂体的顶底界面。再根据岩心确定目标段复合砂体的构型元素类型。潮控河口湾沉积体系中复合砂体构型要素包括:潮汐河道、潮汐沙坝、潮坪,其中潮坪分为潮汐泥坪和潮汐沙坪。为有效反演各构型要素之间的空间接触关系,定义隔夹层要素,其岩性为泥岩。最后在复合砂体顶底界面的控制下,逐井识别复合砂体内部的隔层顶底,并制作表达复合砂体各个构型要素的构型要素类别曲线,例如构型要素类别标签分别是:潮汐河道类别定义为5,潮汐沙坝类别定义为4,潮汐沙坪类别定义为3,潮汐泥坪和隔层类别定义为2。
步骤S22:利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线。
步骤S23:利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体。
步骤S24:基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。
由此可见,本申请通过综合利用测井数据、地震数据、反演数据,获取在潮控河口湾类型的沉积环境中,更加具体的目标复合砂体构型信息,即获取目标复合砂体的构型界面信息、构型特征以及物性特征,实现潮控河口湾复合砂体构型深度解析,为油田油气勘探开发、剩余油挖潜提供有力的技术支撑。
参见图3所示,本申请实施例公开了另一种具体的复合砂体构型解析方法,包括:
步骤S31:利用测井曲线建立构型要素类别曲线。
步骤S32:利用傅里叶变换方法对所述测井曲线进行傅里叶变换,得到所述测井曲线的第一频谱,并在频域下利用所述第一频谱融合所述构型要素类别曲线和所述测井曲线,以得到第一融合后曲线,利用所述傅里叶变换方法对所述第一融合后曲线进行傅里叶反变换,以得到岩性模式特征曲线。
本实施例中,所述利用傅里叶变换方法对所述测井曲线进行傅里叶变换,包括:对所述测井曲线进行归一化处理,以得到归一化后的测井曲线,并利用傅里叶变换方法对所述归一化后的测井曲线进行傅里叶变换。
使用测井曲线交会图法优选砂泥可分曲线,确定使用GR在测井曲线上区分砂泥岩;使用傅里叶变换获取构型要素类别曲线的第一频谱,如下所示:
式中,表示测井曲线、GR曲线、孔渗曲线,表示测井曲线GR曲线和孔渗曲线的均值,表示测井曲线GR曲线和孔渗曲线的方差。
岩性模式特征定义为复合砂体垂向上砂泥相对组分变化曲线,通过在频域上融合GR曲线和构型要素类别曲线构建岩性模式特征曲线,具体步骤为:首先对测井曲线进行归一化处理,以得到归一化后的测井曲线,并利用傅里叶变换方法对归一化后的测井曲线进行傅里叶变换,如下所示:
式中,表示深度域的测井曲线信号,如GR曲线、孔渗曲线,表示是频域的测井曲线的频谱,表示复数单位,表示测井曲线信号的基频。
在频域融合构型要素曲线和GR曲线,如下所示:
式中,表示构型要素类别曲线和测井曲线匹配频域融合后的频谱,表示GR、孔渗曲线在不同频率上的振幅,表示构型要素类别曲线在不同频率上的振幅,表示构型要素类别曲线的频谱,表示GR、孔渗曲线的频谱,表示频谱主频带掩码,表示最大频率取值。
最后进行反傅里叶变换得到岩性模式特征曲线,如下所示:
式中,表示深度域的测井曲线信号,如GR曲线、孔渗曲线,表示频域的测井曲线的频谱,表示复数单位,表示测井曲线信号的基频,表示最大频率取值。
步骤S33:利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,得到所述孔渗曲线的第二频谱,并在频域下利用所述第二频谱融合所述构型要素类别曲线和所述孔渗曲线,以得到第二融合后曲线,利用所述傅里叶变换方法对所述第二融合后曲线进行傅里叶反变换,以得到孔渗模式特征曲线。
本实施例中,所述利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,包括:对孔渗曲线进行归一化处理,以得到归一化后的孔渗曲线,并利用傅里叶变换方法对所述归一化后的孔渗曲线进行傅里叶变换。
孔渗模式特征定义为复合砂体垂向上孔渗相对变化曲线,通过在频域上融合孔渗曲线和构型要素类别曲线构建孔渗模式特征曲线,具体步骤为:首先对孔渗曲线进行归一化后傅里叶变换,得到孔渗曲线的第二频谱,然后在频域融合构型要素曲线和孔渗曲线,以得到第二融合后曲线,最后利用傅里叶变换方法对第二融合后曲线进行傅里叶反变换,得到孔渗模式特征曲线。
步骤S34:利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体。
步骤S35:基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。
由此可见,本申请利用述构型要素类别曲线和测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线,并基于岩性模式特征曲线和孔渗模式特征曲线得到岩性模式反演体和孔渗模式反演体,进而可以获取在潮控河口湾类型的沉积环境中,更加具体的目标复合砂体构型信息,即目标复合砂体的构型界面信息、构型特征以及物性特征,实现潮控河口湾复合砂体构型深度解析。
下面对本申请中的技术方案进行说明。使用沉积背景确定为潮控河口湾的Mariann油田的LU层段为研究目标,实施过程中使用的数据包括地震数据、测井数据。在井点处识别Mariann油田LU层段复合砂体顶底,并在复合砂体内部识别隔层顶底,建立单井复合砂体基础框架,并创建构型要素类别曲线。建立Mariann油田的LU层段高精度层序地层格架,在地震数据和测井数据上确定LU层段复合砂体顶底界面LU顶和LU底。确定Mariann油田LU层段确定复合砂体构型要素包括:潮汐河道、潮汐沙坝、潮汐泥坪、潮汐沙坪。逐井识别复合砂体内部的隔层顶底,并制作表达复合砂体各个构型要素的构型要素类别曲线。构型要素类别标签分别是:潮汐河道类别定义为5,潮汐沙坝类别定义为4,潮汐沙坪类别定义为3,潮汐泥坪和隔层类别定义为2。
在单井复合砂体基础框架控制下,使用频域融合构型要素类别曲线和GR曲线,构建岩性模式特征曲线,使用频域融合构型要素类别曲线和孔渗曲线,构建孔渗模式特征曲线。使用测井曲线交会图法优选LU层段复合砂体砂泥可分曲线是GR,使用傅里叶变换获取LU层段内每口井的构型要素类别曲线的频谱。计算岩性模式特征曲线,首先对LU层段内GR曲线进行归一化并进行傅里叶变换,得到GR曲线的频谱,然后在频域融合构型要素曲线和GR曲线,最后反傅里叶变换得到岩性模式特征曲线。计算孔渗模式特征曲线,首先对LU层段内孔渗曲线进行归一化并进行傅里叶变换,得到孔渗曲线的频谱,然后在频域融合构型要素曲线和孔渗曲线,最后反傅里叶变换,得到孔渗模式特征曲线。
迭代反演获取LU层段复合砂体内部单砂体和隔层展布特征。先使用岩性模式特征曲线制作复合砂体岩性模式反演体,在复合砂体岩性模式反演体的基础上,使用孔渗模式特征曲线制作复合砂体孔渗模式反演体。使用LU顶和LU低建立构造模型,地震数据作为第一波形指示模拟低频特征,岩性模式曲线作为高频特征,确定高通频率为300,高截频率为450,使用第一波形指示模拟计算复合砂体岩性模式反演体。再LU顶和LU低建立构造模型,岩性模式反演体作为第二波形指示模拟低频特征,孔渗特征作为高频特征,确定高通频率为250,高截频率为300,使用第二波形指示模拟计算复合砂体孔渗模式反演体。
在LU层段岩性模式反演体上解释构型界面,迁移岩性模式反演体上解释的构型界面到LU层段孔渗模式特征反演体中,完成复合砂体平面和垂向上物性特征解析。首先在LU层段岩性模式反演体上解释构型界面,例如图4所示的一种具体的构型界面解释结果示意图。选择过井剖面,调整反演体色标显示范围为深红、红、黄、蓝的值分别对应5、4、3、2。逐个过井剖面识别并追踪砂体的顶底界面。解释完成后,迁移岩性模式反演体上解释的构型界面到LU层段孔渗模式特征反演体中。然后,制作东-西向LU层段岩性模式反演体解释剖面,分析各构型单元间的垂向叠置关系,得到研究目标的砂体垂向叠置类型模式,例如图5所示的一种具体的砂体垂向叠置类型模式关系示意图。
延构型界面提取孔渗模式反演平面属性和岩性模式反演平面属性,综合两种属性制作LU层段内隔层分布图。其中,I型隔层对应砂体为垂向上砂体分离,表现为复合砂体构型被厚层隔层分开。II型隔层对应砂体为垂向上砂体叠置,表现为复合砂体构型单元间隔层厚度很小。III型隔层对应砂体为垂向上砂体切叠。为复合砂体构型单元垂向严重的叠置。
进一步的,本申请实施例还提供了一种电子设备。图6是根据一示例性实施例示出的电子设备20结构图,图中的内容不能认为是对本申请的使用范围的任何限制。
图6为本申请实施例提供的一种电子设备的结构示意图。具体可以包括:至少一个处理器21、至少一个存储器22、电源23、通信接口24、输入输出接口25和通信总线26。其中,所述存储器22用于存储计算机程序,所述计算机程序由所述处理器21加载并执行,以实现前述任一实施例公开的由电子设备执行的复合砂体构型解析方法中的相关步骤。
本实施例中,电源23用于为电子设备上的各硬件设备提供工作电压;通信接口24能够为电子设备创建与外界设备之间的数据传输通道,其所遵循的通信协议是能够适用于本申请技术方案的任意通信协议,在此不对其进行具体限定;输入输出接口25,用于获取外界输入数据或向外界输出数据,其具体的接口类型可以根据具体应用需要进行选取,在此不进行具体限定。
其中,处理器21可以包括一个或多个处理核心,比如4核心处理器、8核心处理器等。处理器21可以采用DSP(Digital Signal Processing,数字信号处理)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)、PLA(Programmable Logic Array,可编程逻辑阵列)中的至少一种硬件形式来实现。处理器21也可以包括主处理器和协处理器,主处理器是用于对在唤醒状态下的数据进行处理的处理器,也称CPU(Central ProcessingUnit,中央处理器);协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中,处理器21可以在集成有GPU(Graphics Processing Unit,图像处理器),GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中,处理器21还可以包括AI(Artificial Intelligence,人工智能)处理器,该AI处理器用于处理有关机器学习的计算操作。
另外,存储器22作为资源存储的载体,可以是只读存储器、随机存储器、磁盘或者光盘等,其上所存储的资源包括操作系统221、计算机程序222及数据223等,存储方式可以是短暂存储或者永久存储。
其中,操作系统221用于管理与控制电子设备上的各硬件设备以及计算机程序222,以实现处理器21对存储器22中海量数据223的运算与处理,其可以是Windows、Unix、Linux等。计算机程序222除了包括能够用于完成前述任一实施例公开的由电子设备执行的复合砂体构型解析方法的计算机程序之外,还可以进一步包括能够用于完成其他特定工作的计算机程序。数据223除了可以包括电子设备接收到的由外部设备传输进来的数据,也可以包括由自身输入输出接口25采集到的数据等。
进一步的,本申请实施例还公开了一种计算机可读存储介质,所述存储介质中存储有计算机程序,所述计算机程序被处理器加载并执行时,实现前述任一实施例公开的由复合砂体构型解析过程中执行的方法步骤。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明所提供的一种复合砂体构型解析方法、设备及介质进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种复合砂体构型解析方法,其特征在于,包括:
利用测井曲线建立构型要素类别曲线,并利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线;
利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体;
基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征。
2.根据权利要求1所述的复合砂体构型解析方法,其特征在于,所述利用测井曲线建立构型要素类别曲线,包括:
利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线;其中,所述构型元素类型包括潮汐河道、潮汐沙坝、潮汐泥坪、潮汐沙坪中任意一种或几种。
3.根据权利要求2所述的复合砂体构型解析方法,其特征在于,所述利用岩心标定测井曲线以及确定目标复合砂体的构型元素类型,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线,包括:
利用岩心标定测井曲线,并进行井震标定,以便利用预设地震数据和预设测井数据确定目标复合砂体的顶底界面;
利用所述岩心确定所述目标复合砂体的构型元素类型,并基于所述目标复合砂体的顶底界面逐井识别所述目标复合砂体的内部的隔层顶底;
基于所述隔层顶底,并利用所述测井曲线建立与所述构型元素类型分别对应的构型要素类别曲线。
4.根据权利要求1所述的复合砂体构型解析方法,其特征在于,所述利用傅里叶变换方法、所述构型要素类别曲线和所述测井曲线,以得到岩性模式特征曲线和孔渗模式特征曲线,包括:
利用傅里叶变换方法对所述测井曲线进行傅里叶变换,得到所述测井曲线的第一频谱,并在频域下利用所述第一频谱融合所述构型要素类别曲线和所述测井曲线,以得到第一融合后曲线,利用所述傅里叶变换方法对所述第一融合后曲线进行傅里叶反变换,以得到岩性模式特征曲线;
利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,得到所述孔渗曲线的第二频谱,并在频域下利用所述第二频谱融合所述构型要素类别曲线和所述孔渗曲线,以得到第二融合后曲线,利用所述傅里叶变换方法对所述第二融合后曲线进行傅里叶反变换,以得到孔渗模式特征曲线。
5.根据权利要求4所述的复合砂体构型解析方法,其特征在于,所述利用傅里叶变换方法对所述测井曲线进行傅里叶变换,包括:
对所述测井曲线进行归一化处理,以得到归一化后的测井曲线,并利用傅里叶变换方法对所述归一化后的测井曲线进行傅里叶变换;
相应的,所述利用所述傅里叶变换方法对孔渗曲线进行傅里叶变换,包括:
对孔渗曲线进行归一化处理,以得到归一化后的孔渗曲线,并利用傅里叶变换方法对所述归一化后的孔渗曲线进行傅里叶变换。
6.根据权利要求1至5任一项所述的复合砂体构型解析方法,其特征在于,所述利用所述岩性模式特征曲线获取目标复合砂体的岩性模式反演体,包括:
利用预设地震数据和所述岩性模式特征曲线确定第一波形指示的低频特征和高频特征,以便利用所述第一波形指示进行模拟计算,得到目标复合砂体的岩性模式反演体;
相应的,所述基于所述岩性模式反演体,并利用所述孔渗模式特征曲线制作所述目标复合砂体的孔渗模式反演体,包括:
基于所述岩性模式反演体确定第二波形指示的低频特征,并利用所述孔渗模式特征曲线确定所述第二波形指示的高频特征,以便利用所述第二波形指示进行模拟计算,得到所述目标复合砂体的孔渗模式反演体。
7.根据权利要求1所述的复合砂体构型解析方法,其特征在于,所述基于所述岩性模式反演体确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征,基于所述构型界面信息和所述孔渗模式反演体获取所述目标复合砂体的物性特征,包括:
确定过井剖面,调整与所述岩性模式反演体对应的色标显示范围,以便确定所述目标复合砂体的构型界面信息,并获取所述目标复合砂体的构型特征;
迁移所述构型界面信息至所述孔渗模式反演体,以便获取所述目标复合砂体的物性特征。
8.一种电子设备,其特征在于,包括:
存储器,用于保存计算机程序;
处理器,用于执行所述计算机程序,以实现如权利要求1至7任一项所述的复合砂体构型解析方法的步骤。
9.一种计算机可读存储介质,其特征在于,用于存储计算机程序;其中,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的复合砂体构型解析方法的步骤。
CN202310024380.1A 2023-01-09 2023-01-09 复合砂体构型解析方法、设备及介质 Active CN115903026B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310024380.1A CN115903026B (zh) 2023-01-09 2023-01-09 复合砂体构型解析方法、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310024380.1A CN115903026B (zh) 2023-01-09 2023-01-09 复合砂体构型解析方法、设备及介质

Publications (2)

Publication Number Publication Date
CN115903026A true CN115903026A (zh) 2023-04-04
CN115903026B CN115903026B (zh) 2023-05-09

Family

ID=85746690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310024380.1A Active CN115903026B (zh) 2023-01-09 2023-01-09 复合砂体构型解析方法、设备及介质

Country Status (1)

Country Link
CN (1) CN115903026B (zh)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103454685A (zh) * 2013-08-09 2013-12-18 中国石油天然气股份有限公司 利用测井约束波阻抗反演预测砂体厚度的方法和装置
CN103792573A (zh) * 2012-10-26 2014-05-14 中国石油化工股份有限公司 一种基于频谱融合的地震波阻抗反演方法
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
CN105445800A (zh) * 2015-11-13 2016-03-30 中国石油化工股份有限公司 一种厚层砂体顶部分异岩性油藏的识别方法
CN106094019A (zh) * 2016-04-12 2016-11-09 中国石油化工股份有限公司 基于地质信息映射的深度域地层结构反演方法
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
US20180246254A1 (en) * 2015-11-02 2018-08-30 Petrochina Company Limited Method of acquiring rock component content of stratum
CN108680955A (zh) * 2017-12-12 2018-10-19 中国地质大学(武汉) 一种岩性识别方法及识别系统
CN108802812A (zh) * 2017-04-28 2018-11-13 中国石油天然气股份有限公司 一种井震融合的地层岩性反演方法
CN110082835A (zh) * 2019-03-25 2019-08-02 中国石油化工股份有限公司 曲流河点坝疏松砂岩储层物性建模新方法
CN111045114A (zh) * 2019-12-31 2020-04-21 核工业北京地质研究院 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法
CN111856566A (zh) * 2019-04-28 2020-10-30 中国石油天然气股份有限公司 湖相滩坝砂体中薄储层预测方法及装置
CN112034527A (zh) * 2020-09-01 2020-12-04 中国石油大学(华东) 基于岩相组合模型及地震频谱特征的砂泥岩薄互层沉积相分析方法
US20210041596A1 (en) * 2019-08-06 2021-02-11 Exxonmobil Upstream Research Company Petrophysical Inversion With Machine Learning-Based Geologic Priors
CN114152985A (zh) * 2021-12-14 2022-03-08 中国地质大学(北京) 一种确定地下古河道边界及其内部薄砂体厚度的方法
US20220221614A1 (en) * 2021-01-11 2022-07-14 Shandong University Of Science And Technology Analysis method, system and storage media of lithological and oil and gas containing properties of reservoirs

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103792573A (zh) * 2012-10-26 2014-05-14 中国石油化工股份有限公司 一种基于频谱融合的地震波阻抗反演方法
CN103454685A (zh) * 2013-08-09 2013-12-18 中国石油天然气股份有限公司 利用测井约束波阻抗反演预测砂体厚度的方法和装置
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
US20180246254A1 (en) * 2015-11-02 2018-08-30 Petrochina Company Limited Method of acquiring rock component content of stratum
CN105445800A (zh) * 2015-11-13 2016-03-30 中国石油化工股份有限公司 一种厚层砂体顶部分异岩性油藏的识别方法
CN106094019A (zh) * 2016-04-12 2016-11-09 中国石油化工股份有限公司 基于地质信息映射的深度域地层结构反演方法
CN108802812A (zh) * 2017-04-28 2018-11-13 中国石油天然气股份有限公司 一种井震融合的地层岩性反演方法
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN108680955A (zh) * 2017-12-12 2018-10-19 中国地质大学(武汉) 一种岩性识别方法及识别系统
CN110082835A (zh) * 2019-03-25 2019-08-02 中国石油化工股份有限公司 曲流河点坝疏松砂岩储层物性建模新方法
CN111856566A (zh) * 2019-04-28 2020-10-30 中国石油天然气股份有限公司 湖相滩坝砂体中薄储层预测方法及装置
US20210041596A1 (en) * 2019-08-06 2021-02-11 Exxonmobil Upstream Research Company Petrophysical Inversion With Machine Learning-Based Geologic Priors
CN111045114A (zh) * 2019-12-31 2020-04-21 核工业北京地质研究院 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法
CN112034527A (zh) * 2020-09-01 2020-12-04 中国石油大学(华东) 基于岩相组合模型及地震频谱特征的砂泥岩薄互层沉积相分析方法
US20220221614A1 (en) * 2021-01-11 2022-07-14 Shandong University Of Science And Technology Analysis method, system and storage media of lithological and oil and gas containing properties of reservoirs
CN114152985A (zh) * 2021-12-14 2022-03-08 中国地质大学(北京) 一种确定地下古河道边界及其内部薄砂体厚度的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵继龙: "塔北西部玉东地区古近系薄层砂岩储层地球物理特征与预测", 《中国博士学位论文全文数据库 基础科学辑》 *
魏立花等: "测井约束岩性反演关键技术分析", 《天然气地球科学》 *

Also Published As

Publication number Publication date
CN115903026B (zh) 2023-05-09

Similar Documents

Publication Publication Date Title
Chen et al. History matching of the Norne full-field model with an iterative ensemble smoother
CN102681013B (zh) 碳酸盐岩储层空间模型的建立方法及装置
CN105093313B (zh) 一种岩溶型油藏单井油气产能预测方法及装置
CN109388817A (zh) 一种储层裂缝三维建模方法
Dubois et al. Multiscale geologic and petrophysical modeling of the giant Hugoton gas field (Permian), Kansas and Oklahoma, USA
CN110309597A (zh) 基于构型界面的阻流带模型确定方法、装置及存储介质
Singleton Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas
CN115903026B (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
Lyu et al. Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China
Aghchelou et al. Lithofacies estimation by multi-resolution graph-based clustering of petrophysical well logs: Case study of south pars gas field of iran
Sukmono Application of multi-attribute analysis in mapping lithology and porosity in the Pematang-Sihapas groups of Central Sumatra Basin, Indonesia
Bigoni et al. A multi-scale path for the characterization of heterogeneous karst carbonates: How log-to-seismic machine learning can optimize hydrocarbon production
Abdolahi et al. Improving the 3D facies model with the seismic-derived log volumes: a case study from the Asmari Formation in the Hendijan Field, southwest Iran
Limin et al. Case study of fine geological modeling for Mishrif carbonate reservoir of H oilfield in Iraq
Fluckiger* et al. Predicting Reservoir Heterogeneity in the Upper Cretaceous Frontier Formation in the Western Powder River Basin-An Integrated Stratigraphic, Sedimentologic, Petrophysical and Geophysical Study
Abdy et al. An integrated and effective method for well placement using a numerical simulation model: a North Kuwait field case study
Carrasquilla et al. Using facies, data mining and artificial intelligence concepts in the evaluation of a carbonate reservoir in Campos basin, Southeastern Brazil
McGeer et al. Dynamically Conditioned Modeling to Address Development Challenges in a Highly Complex Fractured Basement Reservoir, Yemen
Freeman et al. A New Facies Classification Scheme Using Gamma Ray and Bulk Density Logs, With Multiple Practical Applications in North Kuwait Heavy Oil Fields
Stephens et al. Detailed Sector Model to Reconcile Surveillance Results and High Permeability, Geologic, Features; Example from Kashagan Field, Kazakhstan
CN117345208B (zh) 一种压裂优势区定量表征方法、装置、电子设备及介质
US20230127237A1 (en) Integrating wells into adaptive multi-scale geological modeling
CN111396026B (zh) 基于水平井测井数据的碳酸盐岩高渗层识别方法及装置
Wang et al. Super-large scale horizontal well position optimization method and infilling practice in high water-cut large complex fluvial reservoir based on multidisciplinary innovative techniques

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Sun Yu

Inventor after: Yao Xiutian

Inventor after: Yan Baiquan

Inventor after: Wang Suibao

Inventor after: Han Liyang

Inventor after: Zhang Kexin

Inventor after: Zhang Chaoqian

Inventor before: Sun Yu

Inventor before: Yan Baiquan

Inventor before: Wang Suibao

Inventor before: Han Liyang

Inventor before: Zhang Kexin

Inventor before: Zhang Chaoqian