CN107844637B - 一种油藏数值确定方法及装置 - Google Patents
一种油藏数值确定方法及装置 Download PDFInfo
- Publication number
- CN107844637B CN107844637B CN201710968034.3A CN201710968034A CN107844637B CN 107844637 B CN107844637 B CN 107844637B CN 201710968034 A CN201710968034 A CN 201710968034A CN 107844637 B CN107844637 B CN 107844637B
- Authority
- CN
- China
- Prior art keywords
- grid
- fracture
- matrix
- representing
- grids
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
本申请实施例公开了一种油藏数值确定方法及装置。所述方法包括:基于基质网格与裂缝网格之间的几何连接系数,确定目的油藏中的第一窜流量;其中,第一窜流量表示基质网格与裂缝网格之间的窜流量;基于裂缝网格的裂缝开度和裂缝间距,确定目的油藏中的第二窜流量;其中,第二窜流量表示裂缝网格集合中两个裂缝网格之间的窜流量;分别建立基质网格集合对应的基质流体渗流模型和裂缝网格集合对应的裂缝流体渗流模型;基于基质流体渗流模型和裂缝流体渗流模型,以及第一窜流量和第二窜流量,确定基质网格内的压力和流量,以及裂缝网格内的压力和流量。本申请实施例提供的技术方案,可以提高油藏数值的计算效率。
Description
技术领域
本申请涉及油气田开发技术领域,特别涉及一种油藏数值确定方法及装置。
背景技术
与常规油藏相比,多尺度的裂缝性油藏具有较强的非均质性和储集空间多样化的特征,这些特征造成了油藏复杂的渗流机理和相对较低的采收率,为了高效开发此类油藏,需要采用数值模拟手段对其注水开发过程进行精确模拟,设计合理的开发方案,达到提高采收率的目的。
目前,在多尺度的裂缝性油藏数值模拟中,通常采用常规连续性介质模型表征此类油藏,例如,双重孔隙介质模型、离散裂缝模型等。双重孔隙介质模型将裂缝性油藏的基质系统与裂缝系统看作两个相互交错的连续性多孔介质系统,两个系统间通过串流量系数进行联系,该模型适用于裂缝发育程度较高,连通性较好的网状裂缝油藏,对于裂缝发育程度不高或存在主导流体流动的大尺度裂缝时,数值模拟结果与实际相差较大。离散裂缝模型主要是基于匹配型网格,即将裂缝作为内边界,并以此作为约束面来进行网格剖分。由于裂缝几何形态的复杂性且随机分布的特征,需要采用非结构化网格剖分方法,然而,其剖分过程比较复杂和繁琐,尤其当裂缝间距离或夹角较小时,往往会因网格划分质量较差而导致油藏数值的计算比较困难,可能很难建立真实条件下的离散裂缝模型。
发明内容
本申请实施例的目的是提供一种油藏数值确定方法及装置,以在保证所确定的油藏数值的准确度的前提下,提高油藏数值的计算效率。
为解决上述技术问题,本申请实施例提供一种油藏数值确定方法及装置是这样实现的:
一种油藏数值确定方法,提供有目的油藏中的基质网格集合和裂缝网格集合;其中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格;所述方法包括:
基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;
基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;
分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;
基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。
优选方案中,采用下述公式确定所述目的油藏中的第一窜流量:
其中,ψmf表示所述第一窜流量,CI表示所述基质网格与所述裂缝网格之间的几何连接系数,Km表示所述基质网格的渗透率,μ表示所述目的油藏中流体的粘度,pf表示所述裂缝网格内的压力,pm表示所述基质网格内的压力;Am-f表示所述裂缝网格与所述基质网格之间的接触面的面积,dm-f表示所述基质网格到所述裂缝网格的距离。
优选方案中,采用下述公式确定所述目的油藏中的第二窜流量:
其中,qff表示所述第二窜流量,Tff表示所述裂缝网格集合内两个裂缝网格之间的传导率,和分别表示所述裂缝网格集合内第i个裂缝网格内的压力和第j个裂缝网格内的压力;和分别表示所述第i个裂缝网格和所述第j个裂缝网格的裂缝开度,和分别表示所述第i个裂缝网格和所述第j个裂缝网格的渗透率,μ表示所述目的油藏中流体的粘度,和分别表示所述第i个裂缝网格和所述第j个裂缝网格的裂缝间距。
优选方案中,采用下述公式建立所述基质网格集合对应的基质流体渗流模型:
其中,Km表示所述基质网格的渗透率,μ表示所述目的油藏中流体的粘度,pm表示所述基质网格内的压力,▽表示哈密顿算子,ψmf表示所述第一窜流量,qm表示所述基质网格内的流量。
优选方案中,采用下述公式建立所述裂缝网格集合对应的裂缝流体渗流模型:
其中,Kf表示所述裂缝网格的渗透率,ψmf表示所述第一窜流量,qff表示所述第二窜流量,μ表示所述目的油藏中流体的粘度,pf表示所述裂缝网格内的压力,▽表示哈密顿算子,qm表示所述裂缝网格内的流量。
优选方案中,所述基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量,包括:
基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,建立所述基质网格集合对应的离散基质渗流通用模型和所述裂缝网格集合对应的离散裂缝通用渗流模型;
基于预设边界条件、所述离散基质渗流通用模型和所述离散裂缝通用渗流模型,确定所述基质网格内的压力和所述裂缝网格内的压力;
基于所述基质网格内的压力和所述裂缝网格内的压力,对所述基质网格和所述裂缝网格进行线性耦合处理,得到所述基质网格内的流量和所述裂缝网格内的流量。
优选方案中,采用下述公式建立所述基质网格集合对应的离散基质渗流通用模型:
其中,和分别表示所述基质网格集合中指定地层深度位置处的二维网格平面中,沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格,以及沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格之间的渗透率;其中,所述二维网格平面中x方向与y方向相互垂直;表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格的渗透率;和分别表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格内的压力,以及所述沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格内的压力;μ表示所述目的油藏中流体的粘度,ΔAx和ΔAy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的交界面的面积,Δx和Δy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的间距,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格相连通,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的流量,CIi',j',k表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格之间的几何连接系数,表示所述第k个裂缝网格内的压力。
优选方案中,采用下述公式建立所述裂缝网格集合对应的离散裂缝通用渗流模型:
其中,和分别表示所述裂缝网格集合中第k个裂缝网格与第k-1个裂缝网格、第k+1个裂缝网格之间的渗透率;和分别表示所述第k-1个裂缝网格和所述第k+1个裂缝网格内的压力;表示所述第k个裂缝网格内的流量,表示所述第k个裂缝网格的裂缝开度,和分别表示第k个裂缝网格与第k-1个裂缝网格之间的间距,以及第k个裂缝网格与第k+1个裂缝网格之间的间距。
优选方案中,
当所述基质网格的数量小于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
其中,表示所述基质网格集合中第l个基质网格与相邻的基质网格之间的传导率,l表示所述基质网格集合中基质网格的数量,表示所述裂缝网格集合中第n个裂缝网格与相连通的基质网格之间的传导率,表示所述裂缝网格集合中第n个裂缝网格与相邻的裂缝网格之间的传导率,n表示所述裂缝网格集合中裂缝网格的数量;和分别表示所述第l个基质网格内的压力和流量;和分别表示所述第n个裂缝网格内的压力和流量;
当所述基质网格的数量等于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
当所述基质网格的数量大于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
其中,表示所述基质网格集合中第n个基质网格与相邻的基质网格之间的传导率,n表示所述基质网格集合中基质网格的数量,表示所述裂缝网格集合中第l个裂缝网格与相连通的基质网格之间的传导率,表示所述裂缝网格集合中第l个裂缝网格与相邻的裂缝网格之间的传导率,l表示所述裂缝网格集合中裂缝网格的数量;和分别表示所述第n个基质网格内的压力和流量;和分别表示所述第l个裂缝网格内的压力和流量。
一种油藏数值确定装置,所述装置提供目的油藏中的基质网格集合和裂缝网格集合;其中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格;所述装置包括:第一窜流量确定模块、第二窜流量确定模块、流体渗流模型建立模块和网格数值确定模块;其中,
所述第一窜流量确定模块,用于基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;
所述第二窜流量确定模块,用于基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;
所述流体渗流模型建立模块,用于分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;
所述网格数值确定模块,用于基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。
本申请实施例提供了一种油藏数值确定方法及装置,可以基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;可以基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;可以分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;可以基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。如此,可以在保证所确定的油藏数值的准确度的前提下,提高油藏数值的计算效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种油藏数值确定方法实施例的流程图;
图2是本申请实施例中分别采用现有技术方法和本申请方法确定的油藏数值的对比结果示意图;
图3是本申请一种油藏数值确定装置实施例的组成结构图。
具体实施方式
本申请实施例提供一种油藏数值确定方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
本申请实施例提供了一种油藏数值确定方法。所述油藏数值确定方法提供有目的油藏中的基质网格集合和裂缝网格集合。
在本实施方式中,所述目的油藏可以是指油藏数值尚未确定的油藏。所述目的油藏可以是多尺度裂缝性油藏。
在本实施方式中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格。其中,所述基质网格和所述裂缝网格之间相互连通。所述基质网格和裂缝网格可以是三维网格,也可以是二维网格。
在本实施方式中,可以对所述目的油藏的基质和裂缝进行结构化网格划分处理,得到所述基质网格集合和裂缝网格集合。具体地,首先,可以采用正交网格划分的方法对所述目的油藏的基质进行结构化网格划分处理,得到所述基质网格集合。然后,可以根据所述目的油藏的裂缝与所述基质网格的交点,以及所述目的油藏中两个裂缝之间的交点,将所述裂缝离散划分为所述裂缝网格集合。
图1是本申请一种油藏数值确定方法实施例的流程图。如图1所示,所述油藏数值确定方法,包括以下步骤。
步骤S101:基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量。
在本实施方式中,基于所述基质网格与所述裂缝网格之间的几何连接系数,可以采用下述公式确定所述目的油藏中的第一窜流量:
其中,ψmf表示所述第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;CI表示所述基质网格与所述裂缝网格之间的几何连接系数,Km表示所述基质网格的渗透率,μ表示所述目的油藏中流体的粘度,pf表示所述裂缝网格内的压力,pm表示所述基质网格内的压力;Am-f表示所述裂缝网格与所述基质网格之间的接触面的面积,dm-f表示所述基质网格到所述裂缝网格的距离。
步骤S102:基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量。
在本实施方式中,基于所述裂缝网格的裂缝开度和裂缝间距,可以采用下述公式确定所述目的油藏中的第二窜流量:
其中,qff表示所述第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;Tff表示所述裂缝网格集合内两个裂缝网格之间的传导率,和分别表示所述裂缝网格集合内第i个裂缝网格内的压力和第j个裂缝网格内的压力;和分别表示所述第i个裂缝网格和所述第j个裂缝网格的裂缝开度,和分别表示所述第i个裂缝网格和所述第j个裂缝网格的渗透率,μ表示所述目的油藏中流体的粘度,和分别表示所述第i个裂缝网格和所述第j个裂缝网格的裂缝间距。
步骤S103:分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型。
在本实施方式中,可以采用下述公式建立所述基质网格集合对应的基质流体渗流模型:
其中,Km表示所述基质网格的渗透率,μ表示所述目的油藏中流体的粘度,pm表示所述基质网格内的压力,▽表示哈密顿算子,ψmf表示所述第一窜流量,qm表示所述基质网格内的流量。
在本实施方式中,可以采用下述公式建立所述裂缝网格集合对应的裂缝流体渗流模型:
其中,Kf表示所述裂缝网格的渗透率,ψmf表示所述第一窜流量,qff表示所述第二窜流量,μ表示所述目的油藏中流体的粘度,pf表示所述裂缝网格内的压力,▽表示哈密顿算子,qm表示所述裂缝网格内的流量。
步骤S104:基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。
在本实施方式中,基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量,具体可以包括以下步骤:
(1)基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,建立所述基质网格集合对应的离散基质渗流通用模型和所述裂缝网格集合对应的离散裂缝通用渗流模型;
(2)基于预设边界条件、所述离散基质渗流通用模型和所述离散裂缝通用渗流模型,确定所述基质网格内的压力和所述裂缝网格内的压力;其中,所述预设边界条件具体可以包括:所述目的油藏中的流体的预设流动特征、所述目的油藏中沿x方向两侧的预设边界压力以及所述基质网格和所述裂缝网格内的预设初始压力。所述预设流动特征可以是所述目的油藏中的流体沿x方向流动,且沿y方向不流动。所述基质网格和所述裂缝网格内预设初始压力可以相同。
(3)基于所述基质网格内的压力和所述裂缝网格内的压力,对所述基质网格和所述裂缝网格进行线性耦合处理,得到所述基质网格内的流量和所述裂缝网格内的流量。
在本实施方式中,可以采用下述公式建立所述基质网格集合对应的离散基质渗流通用模型:
其中,和分别表示所述基质网格集合中指定地层深度位置处的二维网格平面中,沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格,以及沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格之间的渗透率;其中,所述二维网格平面中x方向与y方向相互垂直;表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格的渗透率;和分别表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格内的压力,以及所述沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格内的压力;μ表示所述目的油藏中流体的粘度,ΔAx和ΔAy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的交界面的面积,Δx和Δy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的间距,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格相连通,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的流量,CIi',j',k表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格之间的几何连接系数,表示所述第k个裂缝网格内的压力。
在本实施方式中,可以采用下述公式建立所述裂缝网格集合对应的离散裂缝通用渗流模型:
其中,和分别表示所述裂缝网格集合中第k个裂缝网格与第k-1个裂缝网格、第k+1个裂缝网格之间的渗透率;和分别表示所述第k-1个裂缝网格和所述第k+1个裂缝网格内的压力;表示所述第k个裂缝网格内的流量,表示所述第k个裂缝网格的裂缝开度,和分别表示第k个裂缝网格与第k-1个裂缝网格之间的间距,以及第k个裂缝网格与第k+1个裂缝网格之间的间距。
在本实施方式中,当所述基质网格的数量小于所述裂缝网格的数量时,可以采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
其中,表示所述基质网格集合中第l个基质网格与相邻的基质网格之间的传导率,l表示所述基质网格集合中基质网格的数量,表示所述裂缝网格集合中第n个裂缝网格与相连通的基质网格之间的传导率,表示所述裂缝网格集合中第n个裂缝网格与相邻的裂缝网格之间的传导率,n表示所述裂缝网格集合中裂缝网格的数量;和分别表示所述第l个基质网格内的压力和流量;和分别表示所述第n个裂缝网格内的压力和流量。或者,
当所述基质网格的数量等于所述裂缝网格的数量时,可以采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
当所述基质网格的数量大于所述裂缝网格的数量时,可以采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
其中,表示所述基质网格集合中第n个基质网格与相邻的基质网格之间的传导率,n表示所述基质网格集合中基质网格的数量,表示所述裂缝网格集合中第l个裂缝网格与相连通的基质网格之间的传导率,表示所述裂缝网格集合中第l个裂缝网格与相邻的裂缝网格之间的传导率,l表示所述裂缝网格集合中裂缝网格的数量;和分别表示所述第n个基质网格内的压力和流量;和分别表示所述第l个裂缝网格内的压力和流量。
例如,图2是本申请实施例中分别采用现有技术方法和本申请方法确定的油藏数值的对比结果示意图。其中,现有技术的方法为离散裂缝模型方法。图2中的(a)为分别采用现有技术方法和本申请方法确定的目的油藏中的基质网格内的压力对比结果示意图。图2中的(b)为分别采用现有技术方法和本申请方法确定的目的油藏中的裂缝网格内的压力对比结果示意图。图2中的(a)和(b)中的横坐标和纵坐标分别表示所述目的油藏沿x方向的位置和该位置处的压力,压力的单位为兆帕。如图2所示,采用本申请的方法所确定的油藏数值结果与采用离散裂缝模型的方法所确定的油藏数值结果比较一致,表明采用本申请的方法可以保证所确定的油藏数值结果的准确度。不仅如此,由于本申请的方法采用结构化网格划分方法,例如,正交网格划分方法,对所述目的油藏中的基质和裂缝进行网格划分,所划分的基质网格和裂缝网格数量较少,从而可以提高油藏数值计算的效率。
所述油藏数值确定方法实施例,可以基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;可以基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;可以分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;可以基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。如此,可以在保证所确定的油藏数值的准确度的前提下,提高油藏数值的计算效率。
图3是本申请一种油藏数值确定装置实施例的组成结构图。所述油藏数值确定装置提供目的油藏中的基质网格集合和裂缝网格集合;其中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格。如图3所示,所述油藏数值确定装置可以包括:第一窜流量确定模块100、第二窜流量确定模块200、流体渗流模型建立模块300和网格数值确定模块400。
所述第一窜流量确定模块100,可以用于基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量。
所述第二窜流量确定模块200,可以用于基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量。
所述流体渗流模型建立模块300,可以用于分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型。
所述网格数值确定模块400,可以用于基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量。
所述油藏数值确定装置实施例与所述油藏数值确定方法实施例相对应,可以实现油藏数值确定方法实施例的技术方案,并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的装置、模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (8)
1.一种油藏数值确定方法,其特征在于,提供有目的油藏中的基质网格集合和裂缝网格集合;其中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格;所述方法包括:
基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;
基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;
分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;
基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量;其中,包括:基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,建立所述基质网格集合对应的离散基质渗流通用模型和所述裂缝网格集合对应的离散裂缝通用渗流模型;基于预设边界条件、所述离散基质渗流通用模型和所述离散裂缝通用渗流模型,确定所述基质网格内的压力和所述裂缝网格内的压力;基于所述基质网格内的压力和所述裂缝网格内的压力,对所述基质网格和所述裂缝网格进行线性耦合处理,得到所述基质网格内的流量和所述裂缝网格内的流量;
其中,采用下述公式建立所述基质网格集合对应的离散基质渗流通用模型:
其中,和分别表示所述基质网格集合中指定地层深度位置处的二维网格平面中,沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格,以及沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格之间的渗透率;其中,所述二维网格平面中x方向与y方向相互垂直;表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格的渗透率;和分别表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格内的压力,以及所述沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格内的压力;μ表示所述目的油藏中流体的粘度,ΔAx和ΔAy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的交界面的面积,Δx和Δy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的间距,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格相连通,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的流量,CIi',j',k表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格之间的几何连接系数,表示所述第k个裂缝网格内的压力。
7.根据权利要求1所述的一种油藏数值确定方法,其特征在于,
当所述基质网格的数量小于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
其中,表示所述基质网格集合中第l个基质网格与相邻的基质网格之间的传导率,l表示所述基质网格集合中基质网格的数量,表示所述裂缝网格集合中第n个裂缝网格与相连通的基质网格之间的传导率,表示所述裂缝网格集合中第n个裂缝网格与相邻的裂缝网格之间的传导率,n表示所述裂缝网格集合中裂缝网格的数量;和分别表示所述第l个基质网格内的压力和流量;和分别表示所述第n个裂缝网格内的压力和流量;
当所述基质网格的数量等于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
当所述基质网格的数量大于所述裂缝网格的数量时,采用下述公式对所述基质网格和所述裂缝网格进行线性耦合处理:
8.一种油藏数值确定装置,其特征在于,所述装置提供目的油藏中的基质网格集合和裂缝网格集合;其中,所述基质网格集合中包括多个基质网格,所述裂缝网格集合中包括多个裂缝网格;所述装置包括:第一窜流量确定模块、第二窜流量确定模块、流体渗流模型建立模块和网格数值确定模块;其中,
所述第一窜流量确定模块,用于基于所述基质网格与所述裂缝网格之间的几何连接系数,确定所述目的油藏中的第一窜流量;其中,所述第一窜流量表示所述基质网格与所述裂缝网格之间的窜流量;
所述第二窜流量确定模块,用于基于所述裂缝网格的裂缝开度和裂缝间距,确定所述目的油藏中的第二窜流量;其中,所述第二窜流量表示所述裂缝网格集合中两个裂缝网格之间的窜流量;
所述流体渗流模型建立模块,用于分别建立所述基质网格集合对应的基质流体渗流模型和所述裂缝网格集合对应的裂缝流体渗流模型;
所述网格数值确定模块,用于基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,确定所述基质网格内的压力和流量,以及所述裂缝网格内的压力和流量;其中,包括:基于所述基质流体渗流模型和所述裂缝流体渗流模型,以及所述第一窜流量和所述第二窜流量,建立所述基质网格集合对应的离散基质渗流通用模型和所述裂缝网格集合对应的离散裂缝通用渗流模型;基于预设边界条件、所述离散基质渗流通用模型和所述离散裂缝通用渗流模型,确定所述基质网格内的压力和所述裂缝网格内的压力;基于所述基质网格内的压力和所述裂缝网格内的压力,对所述基质网格和所述裂缝网格进行线性耦合处理,得到所述基质网格内的流量和所述裂缝网格内的流量;
其中,采用下述公式建立所述基质网格集合对应的离散基质渗流通用模型:
其中,和分别表示所述基质网格集合中指定地层深度位置处的二维网格平面中,沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格、沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格,以及沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格之间的渗透率;其中,所述二维网格平面中x方向与y方向相互垂直;表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格的渗透率;和分别表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'-1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'+1且沿y方向的网格号为j'的基质网格内的压力、所述沿x方向的网格号为i'且沿y方向的网格号为j'-1的基质网格内的压力,以及所述沿x方向的网格号为i'且沿y方向的网格号为j'+1的基质网格内的压力;μ表示所述目的油藏中流体的粘度,ΔAx和ΔAy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的交界面的面积,Δx和Δy分别表示所述基质网格集合中相邻两个基质网格沿x方向和沿y方向的间距,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格相连通,表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格内的流量,CIi',j',k表示所述沿x方向的网格号为i'且沿y方向的网格号为j'的基质网格与第k个裂缝网格之间的几何连接系数,表示所述第k个裂缝网格内的压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710968034.3A CN107844637B (zh) | 2017-10-18 | 2017-10-18 | 一种油藏数值确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710968034.3A CN107844637B (zh) | 2017-10-18 | 2017-10-18 | 一种油藏数值确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107844637A CN107844637A (zh) | 2018-03-27 |
CN107844637B true CN107844637B (zh) | 2021-03-30 |
Family
ID=61662485
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710968034.3A Active CN107844637B (zh) | 2017-10-18 | 2017-10-18 | 一种油藏数值确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107844637B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111079335B (zh) * | 2019-12-20 | 2020-10-30 | 西南石油大学 | 一种嵌入式离散裂缝模型的数学推导和数值计算方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9068448B2 (en) * | 2008-12-03 | 2015-06-30 | Chevron U.S.A. Inc. | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
CN102339326B (zh) * | 2010-07-16 | 2014-01-15 | 中国石油化工股份有限公司 | 一种分析模拟缝洞型油藏流体流动的方法 |
CN103399970B (zh) * | 2013-04-16 | 2016-02-03 | 西南石油大学 | 用离散裂缝线处理进行数模测算油藏流动状况的方法 |
US9416630B2 (en) * | 2013-06-18 | 2016-08-16 | Baker Hughes Incorporated | Grid calculation for fluid flow modeling |
CN104392109B (zh) * | 2014-11-04 | 2017-10-17 | 中国石油天然气股份有限公司 | 一种油藏的动态离散裂缝模拟方法和装置 |
CN105653746B (zh) * | 2014-11-10 | 2019-01-18 | 中国石油化工股份有限公司 | 基于嵌入式离散裂缝模型的压裂井建模及模拟方法 |
CN105089595B (zh) * | 2015-05-27 | 2019-04-12 | 中国石油天然气股份有限公司 | 水平压裂裂缝导流作用下的油藏数值模拟方法及装置 |
CN106326517B (zh) * | 2015-07-03 | 2019-07-30 | 中国石油化工股份有限公司 | 层状裂缝-基质杂交网格建模方法和装置 |
CN107060746B (zh) * | 2017-04-27 | 2019-03-26 | 中国石油大学(华东) | 一种复杂裂缝性油藏流动模拟的方法 |
-
2017
- 2017-10-18 CN CN201710968034.3A patent/CN107844637B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107844637A (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108106964B (zh) | 一种确定页岩含气量的方法及装置 | |
CN108019205B (zh) | 一种确定储层可压裂性指数的方法及装置 | |
US9361414B2 (en) | Hybrid local nonmatching method for multiphase flow simulations in heterogeneous fractured media | |
EP2791841A1 (en) | Relaxed constraint delaunay method for discretizing fractured media | |
CN107657075B (zh) | 模拟地下水介质交界面处达西速度的区域分解有限元法 | |
Sun et al. | Anisotropic and dynamic mesh adaptation for discontinuous Galerkin methods applied to reactive transport | |
CN108227032B (zh) | 一种地应力的确定方法及装置 | |
CN106845786B (zh) | 注入井和生产井间储层连通性评价方法和装置 | |
CN107844637B (zh) | 一种油藏数值确定方法及装置 | |
CN111025409B (zh) | 一种水淹层评价方法、装置及存储介质 | |
Zhang et al. | An efficient mimetic finite difference method for multiphase flow in fractured reservoirs | |
Majumder et al. | Vectorized simulation of groundwater flow and contaminant transport using analytic element method and random walk particle tracking | |
Li et al. | Tracer flowback based fracture network characterization in hydraulic fracturing | |
CN109612900B (zh) | 一种储层岩石渗透率预测方法、装置及存储介质 | |
CN108595782B (zh) | 一种离散裂缝中基质与裂缝间的传质计算方法 | |
Croucher | Recent developments in the PyTOUGH scripting library for TOUGH2 simulations | |
CN108505991B (zh) | 一种确定岩心中油的采出程度的方法及装置 | |
CN108267561B (zh) | 室内恒速实验注入速度的确定方法及装置 | |
Gao et al. | Hybrid C-and L-Moment–based Hermite transformation models for non-Gaussian processes | |
Rasoulzadeh et al. | Hydrodynamic driven dissolution in porous media with embedded cavities | |
CN102254093A (zh) | 基于泰森多边形的连通域统计相关算法 | |
CN107563617B (zh) | 一种确定油气运移路径上的油气量的方法及装置 | |
Broadbent et al. | Mathematical Proceedings of the Cambridge Philosophical Society | |
CN103279354A (zh) | 设计规则检查文件的自动化产生方法及产生器 | |
CN107605472B (zh) | 一种确定油藏采收率的方法及装置 |
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 |