CN112329187B - 断层封闭性分析方法、装置及系统 - Google Patents
断层封闭性分析方法、装置及系统 Download PDFInfo
- Publication number
- CN112329187B CN112329187B CN201910649401.2A CN201910649401A CN112329187B CN 112329187 B CN112329187 B CN 112329187B CN 201910649401 A CN201910649401 A CN 201910649401A CN 112329187 B CN112329187 B CN 112329187B
- Authority
- CN
- China
- Prior art keywords
- core
- fracture zone
- pores
- fault
- seepage
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种断层封闭性分析方法、装置及系统,该方法包括:获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。本发明可以定量分析断层封闭性,准确度高。
Description
技术领域
本申请涉及油气勘探开发技术领域,特别涉及一种断层封闭性分析方法、装置及系统。
背景技术
断层封闭性是形成油气藏及控制油气藏规模的重要因素,它的研究对于了解断层在油气成藏和石油勘探开发中具有重要意义。因此,对断层封闭性的研究显得很重要。同时,影响断层封闭性的因素有很多,对这些因素的分析也是十分必要的。断层在油气运聚过程中究竟起何种作用,关键取决于其封闭性。所谓断层封闭性是指断层面或断裂带对地层流体封堵并阻止流体渗流的能力。断层封闭性的研究,包括断层封闭性的定性和定量(或半定量)研究,一直以来都是油气勘探领域研究的热点问题,目前已成为涉及到与断层有关的油气成藏条件研究和减少勘探风险的一个重要指标。尤其是近几年,海域勘探的高风险性促进了断层封闭性定量评价的步伐。断层封闭性的研究自20世纪50~60年代起步以来,已有半个多世纪之久,取得了一系列进展和可喜的成果。国内外对断层封闭性的研究,主要侧重于断层的几何学、形态学、断层面的物质涂抹及动力学、运动学等方面,偏向于宏观定性研究较多,准确度不高。
发明内容
本发明实施例提出一种断层封闭性分析方法,用以定量分析断层封闭性,准确度高,该方法包括:
获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
根据孔隙数据体,获得岩心断裂带渗流场;
计算岩心断裂带渗流场在断裂带处的强度;
根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
本发明实施例提出一种断层封闭性分析装置,用以定量分析断层封闭性,准确度高,该装置包括:
第一数据获得模块,用于获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
第二数据获得模块,用于对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
第三数据获得模块,用于对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
岩心断裂带渗流场获得模块,用于根据孔隙数据体,获得岩心断裂带渗流场;
计算模块,用于计算岩心断裂带渗流场在断裂带处的强度;
分析模块,用于根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
本发明实施例提出一种断层封闭性分析系统,用以定量分析断层封闭性,准确度高,该系统包括:上述断层封闭性分析装置,扫描单元,其中,
所述扫描单元,用于:
对含有断层的目标层段岩心模型进行扫描,获得含有孔隙的第一岩心数据体;
将含有孔隙的第一岩心数据体发送至第一数据获得模块。
本发明实施例还提出了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述断层封闭性分析方法。
本发明实施例还提出了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述断层封闭性分析方法的计算机程序。
在本发明实施例中,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。在上述过程中,获得了孔隙数据体,构造了岩心断裂带渗流场,并计算了岩心断裂带渗流场在断裂带处的强度,然后根据岩心断裂带渗流场在断裂带处的强度的具体数据,分析断层封闭性的过程即为定量分析断层封闭性的过程,相对于定性分析,准确度高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中断层封闭性分析方法的流程图;
图2为本发明实施例提出的断层封闭性分析方法的详细流程图;
图3为本发明实施例中含有孔隙的第一岩心数据体的示意图;
图4为本发明实施例中实心的第二岩心数据体的示意图;
图5为本发明实施例中孔隙数据体的示意图;
图6为本发明实施例中渗流联通通道的示意图;
图7为本发明实施例中渗流线的示意图;
图8为本发明实施例中岩心断裂带渗流场在断裂带处的强度的示意图;
图9为本发明实施例中断层封闭性分析装置的示意图;
图10为本发明实施例中断层封闭性分析系统的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
在现有技术中,国内外对断层封闭性的研究,主要侧重于断层的几何学、形态学、断层面的物质涂抹及动力学、运动学等方面,偏向于宏观定性研究较多,准确度不高。针对上述问题,发明人考虑将数字岩心技术应用到断层封闭性分析中,数字岩心技术是近年兴起的岩心分析的有效方法,在常规砂岩和碳酸盐岩等岩心分析领域应用广泛,获得了极大的成功,基本原理是基于二维扫描电镜图像或三维CT扫描图像,运用计算机图像处理技术,通过一定的算法完成数字岩心重构。数字岩心建模均借助高倍光学显微镜、扫描电镜或CT成像仪等高精度仪器获取岩心的平面图像,之后对平面图像进行三维重建即可得到数字岩心;数值重建方法则借助岩心平面图像等少量资料,通过图像分析提取建模信息,之后采用某种数学方法建立数字岩心。早期的数字岩心是基于二维扫描电镜图片,通过数值算法实现三维重构。数值算法包括模拟退火算法、过程法、马尔科夫-蒙特卡洛等,模拟退火算法首先随机产生孔隙度为φ的随机多孔介质,通过不断调整孔隙和骨架的位置,产生符合条件的多孔介质。过程法模拟了真实岩心形成过程,包括沉积、压实和成岩作用,重构的数字岩心有较好的连通性。马尔科夫-蒙特卡洛方法首先通过邻域模板对原始图片进行遍历,获取条件概率函数,然后用蒙特卡洛算法确定重构图像的每一点的状态。技术上述数字岩心技术,发明提出了一种断层封闭性分析方法。
图1为本发明实施例中断层封闭性分析方法的流程图,如图1所示,该方法包括:
步骤101,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
步骤102,对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
步骤103,对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
步骤104,根据孔隙数据体,获得岩心断裂带渗流场;
步骤105,计算岩心断裂带渗流场在断裂带处的强度;
步骤106,根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
在本发明实施例中,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。在上述过程中,获得了孔隙数据体,构造了岩心断裂带渗流场,并计算了岩心断裂带渗流场在断裂带处的强度,然后根据岩心断裂带渗流场在断裂带处的强度的具体数据,分析断层封闭性的过程即为定量分析断层封闭性的过程,相对于定性分析,准确度高。
具体实施时,首先获得含有断层的目标层段岩心模型,含有断层的目标层段岩心模型包括泥岩涂抹类型断层,砂泥岩对接型断层等多种类型。含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的,扫描方式可以有多种,例如,可以采用高精度工业CT扫描,根据具体需求设置合理扫描步长和范围,扫描范围需包含断裂带结构,保证断裂带附近清晰扫描,当然,还可以采用其他三维扫描方式进行扫描。实心的第二岩心数据体不含空隙,因此,对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得的即为岩心数据体中的孔隙,称之为空隙数据体。根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度,该强度为一个具体的定量的数值,根据该强度,即可定量分析断层的封闭性。
在一实施例中,含有孔隙的第一岩心数据体包括岩心骨架体积数据和/或孔隙体积数据。
具体实施时,对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体的方法有多种,下面给出其中一个实施例。
在一实施例中,对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体,包括:
对含有孔隙的第一岩心数据体中的孔隙进行重建填充,获得实心的第二岩心数据体。
在上述实施例中,含有孔隙的第一岩心数据体为一个三维结构,三维重建即对孔隙进行重建填充,获得实心的第二岩心数据体。
具体实施时,根据孔隙数据体,获得岩心断裂带渗流场的方法有多种,下面给出其中一个实施例。
在一实施例中,根据孔隙数据体,获得岩心断裂带渗流场,包括:
获得孔隙数据体的孔隙三维成像数据;
根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙;
确定渗流联通通道的起点和终点;
根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道;
根据所述多个渗流联通通道,构造岩心断裂带渗流场。
在上述实施例中,孔隙三维成像数据可以通过将孔隙数据体加载进专业的成像软件获得,例如,VG或者Simpleware软件,渗流联通通道的起点一般为断层一侧,渗流联通通道的终点为断层的另一侧,然后连通这些联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道,将形成的垂直通过断层方向的渗流联通通道按照顺序进行标记,通道N1、N2、N3……,便于后续步骤分析。
具体实施时,根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙的方法有多种,下面给出其中一个实施例。
在一实施例中,根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙,包括:
根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙,包括:
确定孔隙尺寸和孔隙间的距离;
根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙。
在上述实施例中,可根据空隙联通原理,设定空隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙。
具体实施时,计算岩心断裂带渗流场在断裂带处的强度的方法有多种,下面给出其中一个实施例。
在一实施例中,计算岩心断裂带渗流场在断裂带处的强度,包括:
获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线;
根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度。
在上述实施例中,对每一渗流联通通道,从该渗流联通通道的起点到渗流联通通道的终点,获得该渗流联通通道的渗流线L1、L2、L3……,其中,渗流线的疏密表示渗流场的强度,渗流线从起点穿过断层到终点的方向为标定的渗流场的方向。
在一实施例中,采用如下公式,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度:
W=L/S
W为岩心断裂带渗流场在断裂带处的强度;
L为渗流线的数量;
S为断层的断面面积。
在上述实施例中,岩心断裂带渗流场在断裂带附近强度系数W越强的地方断层封闭性越弱,岩心断裂带渗流场在断裂带强度系数W越弱的地方或者无渗流线的地方断层封闭性越强。
根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性的方法有多种,下面给出其中一个实施例。
在一实施例中,根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性,包括:
根据钻井数据,确定岩心断裂带渗流场的强度阈值;
将岩心断裂带渗流场在断裂带处的强度与所述强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭。
在上述实施例中,将岩心断裂带渗流场在断裂带处的强度W与强度阈值W1相比较,W大于W1处表示渗流能力强,断层对流体阻碍能力不强,流体易通过,确定岩心断裂带渗流场在断裂带处不封闭,W小于W1处表示渗流能力弱,断层对流体阻碍能力强,流体不易通过,确定岩心断裂带渗流场在断裂带处封闭。
在一实施例中,可根据目标断层的空间断面中各处的封闭性,确定是否对所述目标断层控制下的断块圈闭进行钻井。
基于上述实施例,本发明提出如下一个实施例来说明断层封闭性分析方法的详细流程,图2为本发明实施例提出的断层封闭性分析方法的详细流程图,如图2所示,在一实施例中,断层封闭性分析方法的详细流程包括:
步骤201,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
步骤202,对含有孔隙的第一岩心数据体中的孔隙进行重建填充,获得实心的第二岩心数据体;
步骤203,对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
步骤204,获得孔隙数据体的孔隙三维成像数据;
步骤205,确定孔隙尺寸和孔隙间的距离;
步骤206,根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙;
步骤207,确定渗流联通通道的起点和终点;
步骤208,根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道;
步骤209,根据所述多个渗流联通通道,构造岩心断裂带渗流场;
步骤210,获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线;
步骤211,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度;
步骤212,根据钻井数据,确定岩心断裂带渗流场的强度阈值;
步骤213,将岩心断裂带渗流场在断裂带处的强度与所述强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭。
当然,可以理解的是,上述断层封闭性分析方法的详细流程还可以有其他变化例,相关变化例均应落入本发明的保护范围。
下面给出一具体实施例,说明本发明提出的断层封闭性分析方法的具体应用。
获得含有目标断层A、目标断层B的目标层段岩心,长度为80mm,岩心变径为79mm,为了保证扫描效果以及满足扫描平台尺寸要求,对目标层段岩心做必要的处理,处理为长80mm,宽高分别为40mm的柱状体,形成目标层段岩心模型,其中,切去外部多余砂岩,可使扫描射线更容易穿透外表到达内部断层部位。初步观察识别出目标层段岩心中的断层为逆断层,岩性对接方式为砂泥互接。
将准备好的含有断层A和断层B的目标层段岩心模型放置于扫描平台上,采用高精度工业CT扫描仪进行数据采集,设置扫描方式为三维体积扫描,扫描长度为80mm,扫描宽度为50mm,为了保证扫描效果,设置合理扫毛电压和电流分别为420kv和1.2A,通过扫描获得的含有孔隙的第一岩心数据体中,包含了岩心骨架体积和孔隙体积,断裂带附近数据特征为泥岩致密层。图3为本发明实施例中含有孔隙的第一岩心数据体的示意图,对含有孔隙的第一岩心数据体利用专业三维重建VG软件进行三维重建,获得实心的第二岩心数据体,图4为本发明实施例中实心的第二岩心数据体的示意图,对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体,图5为本发明实施例中孔隙数据体的示意图。
将所述孔隙数据体加载进VG等专业成像软件,得到孔隙三维成像数据,确定孔隙尺寸和孔隙间的距离;根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙;确定渗流联通通道的起点和终点;根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道,图6为本发明实施例中渗流联通通道的示意图,根据所述多个渗流联通通道,构造岩心断裂带渗流场;获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线,图7为本发明实施例中渗流线的示意图,渗流线的疏密表示渗流场的强度,渗流线从起点穿过断层到终点的方向为标定的渗流场的方向,渗流线的疏密表示渗流场的强度,渗流线从起点穿过断层到终点的方向为标定的渗流场的方向,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度,图8为本发明实施例中岩心断裂带渗流场在断裂带处的强度的示意图。
根据钻井数据,确定岩心断裂带渗流场的强度阈值为2.8,将岩心断裂带渗流场在断裂带处的强度与所述强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭,表1为本发明实施例中岩心断裂带渗流场在断裂带处封闭的分析结果,从表1中可以看出,岩心断裂带渗流场在31个断裂带处封闭的分析结果。
表1岩心断裂带渗流场在断裂带处封闭的分析结果
在本发明实施例提出的方法中,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。在上述过程中,获得了孔隙数据体,构造了岩心断裂带渗流场,并计算了岩心断裂带渗流场在断裂带处的强度,然后根据岩心断裂带渗流场在断裂带处的强度的具体数据,分析断层封闭性的过程即为定量分析断层封闭性的过程,相对于定性分析,准确度高。
基于同样的发明构思,本发明实施例还提供了一种断层封闭性分析装置,如下面的实施例所述。由于这些解决问题的原理与断层封闭性分析方法相似,因此装置的实施可以参见方法的实施,重复之处不在赘述。
图9为本发明实施例中断层封闭性分析装置的示意图,如图9所示,该装置包括:
第一数据获得模块901,用于获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
第二数据获得模块902,用于对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
第三数据获得模块903,用于对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
岩心断裂带渗流场获得模块904,用于根据孔隙数据体,获得岩心断裂带渗流场;
计算模块905,用于计算岩心断裂带渗流场在断裂带处的强度;
分析模块906,用于根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
在一实施例中,第二数据获得模块902具体用于:
对含有孔隙的第一岩心数据体中的孔隙进行重建填充,获得实心的第二岩心数据体。
在一实施例中,岩心断裂带渗流场获得模块904具体用于:
获得孔隙数据体的孔隙三维成像数据;
根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙;
确定渗流联通通道的起点和终点;
根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道;
根据所述多个渗流联通通道,构造岩心断裂带渗流场。
在一实施例中,岩心断裂带渗流场获得模块904具体用于:
确定孔隙尺寸和孔隙间的距离;
根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙。
在一实施例中,计算模块905具体用于:
获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线;
根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度。
在一实施例中,分析模块906具体用于:
根据钻井数据,确定岩心断裂带渗流场的强度阈值;
将岩心断裂带渗流场在断裂带处的强度与所述强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭。
在一实施例中,计算模块905具体用于:
采用如下公式,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度:
W=L/S
W为岩心断裂带渗流场在断裂带处的强度;
L为渗流线的数量;
S为断层的断面面积。
在一实施例中,含有孔隙的第一岩心数据体包括岩心骨架体积数据和/或孔隙体积数据。
在本发明实施例提出的装置中,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。在上述过程中,获得了孔隙数据体,构造了岩心断裂带渗流场,并计算了岩心断裂带渗流场在断裂带处的强度,然后根据岩心断裂带渗流场在断裂带处的强度的具体数据,分析断层封闭性的过程即为定量分析断层封闭性的过程,相对于定性分析,准确度高。
本发明实施例提出一种断层封闭性分析系统,图10为本发明实施例中断层封闭性分析系统的示意图,该系统包括:
上述断层封闭性分析装置1001,扫描单元1002,其中,
所述扫描单元1002,用于:
对含有断层的目标层段岩心模型进行扫描,获得含有孔隙的第一岩心数据体;
将含有孔隙的第一岩心数据体发送至第一数据获得模块901。
在本发明实施例提出的系统中,获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;根据孔隙数据体,获得岩心断裂带渗流场;计算岩心断裂带渗流场在断裂带处的强度;根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。在上述过程中,获得了孔隙数据体,构造了岩心断裂带渗流场,并计算了岩心断裂带渗流场在断裂带处的强度,然后根据岩心断裂带渗流场在断裂带处的强度的具体数据,分析断层封闭性的过程即为定量分析断层封闭性的过程,相对于定性分析,准确度高。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种断层封闭性分析方法,其特征在于,包括:
获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
根据孔隙数据体,获得岩心断裂带渗流场;
计算岩心断裂带渗流场在断裂带处的强度;
根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
2.如权利要求1所述的断层封闭性分析方法,其特征在于,对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体,包括:
对含有孔隙的第一岩心数据体中的孔隙进行重建填充,获得实心的第二岩心数据体。
3.如权利要求1所述的断层封闭性分析方法,其特征在于,根据孔隙数据体,获得岩心断裂带渗流场,包括:
获得孔隙数据体的孔隙三维成像数据;
根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙;
确定渗流联通通道的起点和终点;
根据渗流联通通道的起点和终点,联通孔隙和能够在设定压力下形成联通的孔隙,形成多个垂直通过断层方向的渗流联通通道;
根据所述多个渗流联通通道,构造岩心断裂带渗流场。
4.如权利要求3所述的断层封闭性分析方法,其特征在于,根据孔隙三维成像数据,获得联通孔隙和能够在设定压力下形成联通的孔隙,包括:
确定孔隙尺寸和孔隙间的距离;
根据孔隙尺寸和孔隙间的距离,从孔隙三维成像数据中,查找联通孔隙和能够在设定压力下形成联通的孔隙。
5.如权利要求1所述的断层封闭性分析方法,其特征在于,计算岩心断裂带渗流场在断裂带处的强度,包括:
获得岩心断裂带渗流场中多个渗流联通通道对应的渗流线;
根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度。
6.如权利要求1所述的断层封闭性分析方法,其特征在于,根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性,包括:
根据钻井数据,确定岩心断裂带渗流场的强度阈值;
将岩心断裂带渗流场在断裂带处的强度与所述强度阈值进行比较,若岩心断裂带渗流场在断裂带处的强度大于所述强度阈值,则岩心断裂带渗流场在断裂带处不封闭,否则,岩心断裂带渗流场在断裂带处封闭。
7.如权利要求5所述的断层封闭性分析方法,其特征在于,采用如下公式,根据渗流线的数量和断层的断面面积,获得岩心断裂带渗流场在断裂带处的强度:
W=L/S
W为岩心断裂带渗流场在断裂带处的强度;
L为渗流线的数量;
S为断层的断面面积。
8.如权利要求1所述的断层封闭性分析方法,其特征在于,含有孔隙的第一岩心数据体包括岩心骨架体积数据和/或孔隙体积数据。
9.一种断层封闭性分析装置,其特征在于,包括:
第一数据获得模块,用于获得含有孔隙的第一岩心数据体,所述含有孔隙的第一岩心数据体是对含有断层的目标层段岩心模型进行扫描获得的;
第二数据获得模块,用于对含有孔隙的第一岩心数据体进行三维重建,获得实心的第二岩心数据体;
第三数据获得模块,用于对实心的第二岩心数据体和含有孔隙的第一岩心数据体作差,获得孔隙数据体;
岩心断裂带渗流场获得模块,用于根据孔隙数据体,获得岩心断裂带渗流场;
计算模块,用于计算岩心断裂带渗流场在断裂带处的强度;
分析模块,用于根据岩心断裂带渗流场在断裂带处的强度,分析断层封闭性。
10.一种断层封闭性分析系统,其特征在于,包括:权利要求9所述的断层封闭性分析装置,扫描单元,其中,
所述扫描单元,用于:
对含有断层的目标层段岩心模型进行扫描,获得含有孔隙的第一岩心数据体;
将含有孔隙的第一岩心数据体发送至第一数据获得模块。
11.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8任一项所述方法。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至8任一项所述方法的计算机程序。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910649401.2A CN112329187B (zh) | 2019-07-18 | 2019-07-18 | 断层封闭性分析方法、装置及系统 |
GB2112593.5A GB2595809B (en) | 2019-07-18 | 2020-08-26 | Method, apparatus, and system for quantitative analysis of seal of trap faults |
PCT/CN2020/111223 WO2021008630A1 (zh) | 2019-07-18 | 2020-08-26 | 圈闭断层的封闭性定量分析方法、装置及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910649401.2A CN112329187B (zh) | 2019-07-18 | 2019-07-18 | 断层封闭性分析方法、装置及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112329187A CN112329187A (zh) | 2021-02-05 |
CN112329187B true CN112329187B (zh) | 2022-10-04 |
Family
ID=74319457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910649401.2A Active CN112329187B (zh) | 2019-07-18 | 2019-07-18 | 断层封闭性分析方法、装置及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112329187B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3023954A1 (zh) * | 2014-07-15 | 2016-01-22 | Schlumberger Services Petrol | |
CN105626053A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 通过取心井钻取断裂带开展断裂带研究的方法 |
CN108682020A (zh) * | 2018-04-28 | 2018-10-19 | 中国石油大学(华东) | 岩心微米ct孔隙结构重构方法 |
CN109164026A (zh) * | 2018-07-25 | 2019-01-08 | 中国石油天然气股份有限公司 | 岩石渗流能力评价方法及装置 |
-
2019
- 2019-07-18 CN CN201910649401.2A patent/CN112329187B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3023954A1 (zh) * | 2014-07-15 | 2016-01-22 | Schlumberger Services Petrol | |
CN105626053A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 通过取心井钻取断裂带开展断裂带研究的方法 |
CN108682020A (zh) * | 2018-04-28 | 2018-10-19 | 中国石油大学(华东) | 岩心微米ct孔隙结构重构方法 |
CN109164026A (zh) * | 2018-07-25 | 2019-01-08 | 中国石油天然气股份有限公司 | 岩石渗流能力评价方法及装置 |
Non-Patent Citations (2)
Title |
---|
利用断储排替压力差定量评价断层封闭性――以海拉尔盆地贝尔凹陷贝301区块为例;黄劲松等;《科学技术与工程》;20110508(第13期);全文 * |
断层封闭性研究现状;赵密福;《新疆石油地质》;20040601(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112329187A (zh) | 2021-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3144468B1 (fr) | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter | |
EP2560030B1 (fr) | Méthode pour construire un modèle géologique comportant une mise en position de depot d'unites stratigraphiques | |
EP2103964B1 (fr) | Procédé d'estimation de la densité de fractures dans un milieu rocheux | |
CA2713603C (fr) | Procede pour exploiter un reservoir petrolier a partir d'un calage d'historique optimise | |
US20110310101A1 (en) | Pillar grid conversion | |
US20140114632A1 (en) | Method for modeling a reservoir using 3d multiple-point simulations with 2d training images | |
Ishutov et al. | Resurrection of a reservoir sandstone from tomographic data using three-dimensional printing | |
EA032063B1 (ru) | Системы и способы усовершенствования прямого численного моделирования свойств материала из образцов горной породы и выявления неопределенности в свойствах материала | |
Abadpour et al. | 4D seismic history matching with ensemble Kalman filter-assimilation on Hausdorff distance to saturation front | |
CN111739149B (zh) | 一种岩石ct扫描图像的油水分布连续性修复方法 | |
Dubetz et al. | Characterization of rock pore-size distribution and its effects on wormhole propagation | |
Theologou et al. | An MICP-based pore typing workflow–core scale to log scale | |
CN112329187B (zh) | 断层封闭性分析方法、装置及系统 | |
CN110967738A (zh) | 基于深度学习的地震数据插值处理方法及系统 | |
CN112634428A (zh) | 一种基于双向循环生成网络的多孔介质三维图像重建方法 | |
CN111077053A (zh) | 一种裂缝型多孔介质渗流分析可视化方法 | |
Kamal et al. | Numerical well testing: a method to use transient testing results in reservoir simulation | |
CN112818511B (zh) | 一种基于分形蒙特卡罗算法的多孔材料重构方法 | |
WO2009090522A1 (fr) | Procede, programme et systeme informatique de construction d'un modele geologique 3d | |
CN113640879A (zh) | 基于双网络的储层时移参数预测方法和系统 | |
Franc et al. | Image-based effective medium approximation for fast permeability evaluation of porous media core samples | |
Branets et al. | Capturing geologic complexity in a simulation grid | |
CN107507276A (zh) | 用于任意方向的岩心切片显示的三维数字岩心存储方法 | |
Hu et al. | Correlating recovery efficiency to pore throat characteristics using digital rock analysis | |
Fadakar Alghalandis | Stochastic modelling of fractures in rock masses. |
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 |