CN107832524B - 一种确定储层改造体积的方法、装置及系统 - Google Patents
一种确定储层改造体积的方法、装置及系统 Download PDFInfo
- Publication number
- CN107832524B CN107832524B CN201711083854.0A CN201711083854A CN107832524B CN 107832524 B CN107832524 B CN 107832524B CN 201711083854 A CN201711083854 A CN 201711083854A CN 107832524 B CN107832524 B CN 107832524B
- Authority
- CN
- China
- Prior art keywords
- crack
- elliptical surface
- parameters
- axis
- coordinates
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000012545 processing Methods 0.000 claims abstract description 30
- 230000015654 memory Effects 0.000 claims description 14
- 230000004048 modification Effects 0.000 claims description 13
- 238000012986 modification Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 6
- 230000009466 transformation Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 14
- 238000012544 monitoring process Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 239000003245 coal Substances 0.000 description 5
- 230000004075 alteration Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Remote Sensing (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本申请实施例公开了一种确定储层改造体积的方法、装置及系统，所述方法包括获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；将裂缝等效为椭圆面，根据所述裂缝参数获取裂缝的轮廓坐标；根据所述裂缝的轮廓坐标构建凸包；获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。利用本申请各个实施例，可以提高储层改造体积计算的精确度，有利于更好的认识和优化压裂设计。
Description
技术领域
本发明涉及油气开发技术领域，特别地，涉及一种确定储层改造体积的方法、装置及系统。
背景技术
储层改造体积(Stimulated Reservoir Volume,简称SRV)用以描述非常规储层中水力压裂形成的裂缝网络的形状和尺寸。目前，非常规储层中油气井的产能主要依赖于SRV及其毗邻区域的渗透率，因此，SRV是非常规储层中分析油气井产能的一个重要参数，SRV可用于评价储层改造效果以及预测产能。
现有技术中，应用最为广泛的SRV计算方法是基于微地震裂缝监测，通常采用包含微地震事件点的长方体叠加计算SRV，但是微地震监测不能有效区分裂缝扩展与应力扰动诱发的微地震事件，这样会导致得到的储层改造体积不准确。进一步的，公开号为CN105447238A的中国发明专利公开了一种页岩气水平井压裂改造体积的估算方法，该方法通过试井双对数拟合解释裂缝半长，从而计算SRV；但其在计算SRV时对于整个水平井段采用相同的缝长和缝高，而实际的多级压裂水平井在各压裂段的裂缝方位及尺寸都是不同的，从而影响了储层改造体积估计的精确度。因此，业内需要一种更加精确的估计非常规储层的储层改造体积的方法。
发明内容
本申请实施例的目的在于提供一种确定储层改造体积的方法、装置及系统，可以更加精确直观的确定储层改造体积的尺寸和形状，从而提高储层改造效果评价和产能预测的精确度。
本申请提供的一种储层含气识别的方法、装置及系统是通过包括以下方式实现的：
一种确定储层改造体积的方法，所述方法包括：
获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
根据所述裂缝参数获取裂缝的轮廓坐标；
整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
本申请实施例的确定储层改造体积的方法，所述对所述地面倾角数据进行处理获得裂缝参数，包括：
对所述地面倾角进行处理，获取裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标。
本申请实施例的确定储层改造体积的方法，所述根据所述裂缝参数获取裂缝的轮廓坐标，包括：
将裂缝等效为椭圆面，以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零；
将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；
将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；
将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；
获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；
其中，所述旋转以逆时针旋转为正。
本申请实施例的确定储层改造体积的方法，所述对所述地面倾角数据进行处理获得裂缝参数，包括：
对所述地面倾角数据进行反演获得裂缝参数。
本申请实施例的确定储层改造体积的方法，所述对所述地面倾角数据进行反演获得裂缝参数，包括：
根据预先建立的裂缝模型模拟获得理论地面倾角数据，对理论地面倾角数据与实际采集的地面倾角数据进行拟合，根据拟合结果确定裂缝参数。
另一方面，本申请实施例还提供一种确定储层改造体积的装置，包括：
裂缝参数获取模块，用于获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
裂缝坐标获取模块，用于根据所述裂缝参数获取裂缝的轮廓坐标；
凸包构建模块，用于根据所述裂缝参数获取裂缝的轮廓坐标，整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
储层改造体积确定模块，用于获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
本申请实施例的确定储层改造体积的装置，所述裂缝参数获取模块，包括：
裂缝参数获取单元，用于对所述地面倾角进行处理，获取裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标。
本申请实施例的确定储层改造体积的装置，所述裂缝坐标获取模块，包括：
将裂缝等效为椭圆面，以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零；
将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；
将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；
将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；
获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；
其中，所述旋转以逆时针旋转为正。
本申请实施例的确定储层改造体积的装置，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
根据所述裂缝参数获取裂缝的轮廓坐标；
整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
另一方面，本申请实施例还提供一种确定储层改造体积的系统，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意实施例所述方法的步骤。
本说明书一个或多个实施例提供的一种确定储层改造体积的方法、装置及系统，可以通过对地面倾角数据进行处理获得裂缝参数数据，根据裂缝参数构建凸包，从而利用凸包表征储层改造体积。利用本申请实施例提供的方案，定量的确定储层改造体积的形状及尺寸，提高储层改造体积估计的精确度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本说明书提供的一种确定储层改造体积的方法实施例的流程示意图；
图2为本说明书提供的一个具体实例中H2页岩气井组水平井分段压裂所形成三维裂缝网络示意图；
图3为本说明书提供的一个具体实例中H2-1水平井压裂储层改造体积示意图；
图4为本说明书提供的一个具体实例中H2-2水平井压裂储层改造体积示意图；
图5为本说明书提供的一个具体实例中H2-3水平井压裂储层改造体积示意图；
图6为本说明书提供的一个具体实例中H2-4水平井压裂储层改造体积示意图；
图7为本说明书提供的另一个具体实例中H3页岩气井组水平井分段压裂所形成三维裂缝网络示意图；
图8为本说明书提供的另一个具体实例中H3-1水平井压裂储层改造体积示意图；
图9为本说明书提供的另一个具体实例中H3-2水平井压裂储层改造体积示意图；
图10为本说明书提供的另一个具体实例中H3-3水平井压裂储层改造体积示意图；
图11为本说明书提供的另一个具体实例中Q28煤层气井组多井压裂储层改造体积示意图；
图12为本说明书提供的一种确定储层改造体积的装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。
储层改造体积可用于描述非常规储层中水力压裂形成的裂缝网络的形状和尺寸。在压裂过程中，压裂裂缝可以引起岩石变形，岩石变形形成的变形场向各个方向辐射，可以引起地面的倾角变化，形成地面变形场。压裂形成的变形场是可测且唯一的，并且与储层内裂缝特征相关，通过对地面变形场进行分析，可以较为准确的确定裂缝的几何参数、方位等裂缝参数。本说明书中提供的实施方案，可以通过获取压裂裂缝形成的地面倾角数据，对地面倾角数据进行处理获得精确的压裂裂缝几何参数、方位等裂缝参数信息；然后，可以根据裂缝参数信息构建凸包，利用凸包来定量表征储层改造体积。通过本说明书中提供的方案，可以提高储层改造体积确定的精确度；并且可以通过凸包来定量表征储层改造体积，更加直观地确定出压裂储层改造体积的三维形状和尺寸，从而有利于储层改造效果的精确评价，以及更好的认识和优化压裂设计。
图1是本说明书提供的所述一种确定储层改造体积的方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
具体的一个实施例如图1所示，本说明书提供的一种确定储层改造体积的方法的一个实施例中，所述方法可以包括：
S2、获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
S4、根据所述裂缝参数获取裂缝的轮廓坐标；
S6、整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
S8、根据所述凸包的几何参数确定储层改造体积。
本实施例中，可以通过测斜仪微形变技术获取压裂裂缝形成的地面倾角数据。本说明书的一个实施例中，可以在压裂井周围布置精度极高的一组测斜仪来测量压裂形成的地面倾角数据。其中，测斜仪裂缝监测的具体方案可以根据压裂井的轨迹参数、压裂设计参数来确定。然后根据所述测斜仪裂缝监测方案布置测斜仪监测孔，在所述测斜仪监测孔中下入测斜仪，以采集压裂中裂缝引起的地面倾角数据。
本说明书的一个实施例中，可以对所述地面倾角数据进行处理获得裂缝参数。本说明书的一个或者多个实施例中，所述裂缝参数可以包括裂缝起裂点的三维空间坐标、裂缝的方位角、倾角、缝长、缝高。
本实施例中，可以对采集到的地面倾角数据进行处理获得裂缝参数，例如，可以采用数据拟合的方法获得裂缝参数。本说明书一个实施例中，可以对采集到的地面倾角数据进行反演，获得各压裂层/段裂缝的方位和尺寸。本说明书一个或者多个实施例中，可以预先建立裂缝模型，以模拟不同的裂缝理论上引起的地面变形场；然后将理论上的地面变形场与实际采集到的地面变形场进行比较，通过寻找最佳拟合结果来确定各压裂层/段裂缝参数。
本说明书的一个实施例中，可以建立三维坐标系，根据裂缝参数获得裂缝轮廓坐标。本说明书的一个或者多个实施例中，可以以裂缝的缝长为长轴，以裂缝的缝高为短轴，将裂缝等效为椭圆面，然后根据裂缝参数获得裂缝轮廓坐标。当然，在本说明书的一些实施例中，也可以将裂缝等效成圆形或者矩形等其他形状，这里不做限定。因裂缝的自然形状类似于椭圆形，因此，本说明书实施例中优选的采用将裂缝等效为椭圆形。本说明书的一个或者多个实施例中，可以设定水平井与垂直井的交接点为原点，以原点为起点向正东方向为X轴正向，向正北方向为Y轴正向，垂直井朝向地层表面方向为Z轴正向，建立三维坐标系。以原点为椭圆的中心，裂缝的缝长沿Y轴，缝高沿X轴，Z轴坐标为零，将裂缝等效为椭圆面。
本说明书的一个实施例中，可以根据裂缝参数将裂缝对应的椭圆面进行旋转或者平移，然后，获取旋转或者平移后的椭圆面对应的轮廓坐标，将椭圆面对应的轮廓坐标确定该裂缝的轮廓坐标。本说明书的一个或者多个实施例中，可以将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面，将所述第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面，将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面。然后，获取第三椭圆面的轮廓线的三维坐标，将所述第三椭圆面的轮廓线的三维坐标确定为所述裂缝的轮廓坐标。其中，所述旋转以逆时针旋转为正。
利用上述实施例提供的方法，获得压裂井段中每个裂缝的轮廓坐标。整合该压裂井段中每个裂缝的轮廓坐标，形成对应于该压裂井段的三维压裂缝网。将压裂井段的三维压裂缝网中裂缝轮廓坐标集标记为三维空间点阵集合P，计算得到三维点阵集合P所构成的凸包。其中，所述凸包可以包括在一个实数向量空间V中，对于给定集合P，所有包含P的凸集的交集S被称为P的凸包。本说明书的一个实施例中，对于上述三维空间点阵集合P，可以采用德洛涅三角剖分函数(Delaunay Triangulation)构造三维点阵所建立的德洛涅(Delaunay)三角形面的集合DT；然后，采用凸包函数(convex Hull)计算Delaunay三角形面的集合DT所形成的凸包K，所述凸包K即为三维空间点阵集合P的凸包。
本说明书的一个实施例中，可以获取凸包的几何参数信息，将所述凸包的几何参数信息确定为所述储层改造体积的几何参数信息，从而通过凸包来定量表征储层改造体积。本说明书提供的一个或者多个实施例中，还可以利用MATLAB函数将所述凸包进行可视化输出，从而进行直观的观测储层改造体积的尺寸及形状。
本说明书上述实施例根据地面倾角数据获得裂缝参数，然后，根据裂缝参数构建凸包，从而利用凸包定量的表征储层改造体积，提高储层改造体积估计的精确度。进一步的，还可以可视化凸包，直观地观测压裂储层改造体积的三维形状及尺寸，从而可以更加精确直观的分析储层改造的效果，有利于更好的认识和优化压裂设计。
为了使得本说明书提供的实施例中的方案更加清楚，本说明书还提供了应用上述方案的实际待测区域的具体实例，如表1以及图2至图10所示。表1表示利用本说明书实施例提供的方法确定的页岩储层和煤岩储层压裂改造所形成的储层改造体积数据。
图2表示H2页岩气井组水平井分段压裂经测斜仪监测反演所形成的三维裂缝网络；图3至图6分别表示H2页岩气井组中各水平井储层改造体积，其中，图3表示H2-1水平井储层改造体积，图4表示H2-2水平井储层改造体积，图5表示H2-3水平井储层改造体积，图6表示H2-4水平井储层改造体积。图7表示H3页岩气井组水平井分段压裂经测斜仪监测反演所形成的三维裂缝网络；图8至图10分别表示H3页岩气井组中各水平井储层改造体积，其中，图8表示H3-1水平井储层改造体积，图9表示H3-2水平井储层改造体积，图10表示H3-3水平井储层改造体积。H2井组和H3井组在空间距离上较近，施工规模相当，但由于水平井组钻井方位与水平主应力夹角不同，H3井组钻井方向更接近于最小水平主应力方向，易于形成横切缝并造成缝网。分析表1中的储层改造体积大小以及图3至图6、图8至图9中的储层改造体积形状，可以精确且直观的得出，H3井组的储层改造体积远大于H2井组，因此H3井组的改造效果优于H2井组。
表1页岩和煤岩储层压裂改造所形成的储层改造体积数据表
图11表示利用本说明书上述实施例提供的方法获得的Q28煤层气井组压裂储层改造体积。煤层割理比较发育，水力压裂容易形成复杂裂缝网络，各压裂层经测斜仪微形变监测反演得到水平、垂直相交裂缝网络，其改造体积也因施工规模和地质条件的不同而各异。并且，通过图11中的凸包可以直观的分析各压裂层的储层改造体积，有利于更好的认识和优化对煤层压裂设计。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。
本说明书一个或多个实施例提供的一种确定储层改造体积的方法，可以通过对地面倾角数据进行处理获得裂缝参数数据，根据裂缝参数构建凸包，从而利用凸包表征储层改造体积。利用本申请实施例提供的方案，定量的确定储层改造体积的形状及尺寸，提高储层改造体积估计的精确度。
基于上述所述的确定储层改造体积的方法，本说明书一个或多个实施例还提供一种确定储层改造体积的装置。所述的装置可以包括使用了本说明书实施例所述方法的系统、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图12是本说明书提供的一种确定储层改造体积的装置实施例的模块结构示意图，如图12所示，所述装置可以包括：
裂缝参数获取模块102，可以用于获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
裂缝坐标获取模块104，可以根据所述裂缝参数获取裂缝的轮廓坐标；
凸包构建模块106，可以用于根据所述裂缝参数获取裂缝的轮廓坐标，整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
储层改造体积确定模块108，可以用于获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
当然，参照前述方法实施例描述，所述装置的其他实施例中，所述裂缝参数获取模块102可以包括：
裂缝参数获取单元，用于对所述地面倾角进行处理，获取裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标。
本说明书提供的所述装置的另一个实施例中，所述裂缝坐标获取模块104可以包括：
以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零，将裂缝等效为椭圆面；
将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；
将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；
将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；
获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；
其中，所述旋转以逆时针旋转为正。
本说明书提供的所述装置的另一个实施例中，所述裂缝参数获取模块102可以包括：
对所述地面倾角数据进行反演获得裂缝参数。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书一个或多个实施例提供的一种确定储层改造体积的装置，可以通过对地面倾角数据进行处理获得裂缝参数数据，根据裂缝参数构建凸包，从而利用凸包表征储层改造体积。利用本申请实施例提供的方案，定量的确定储层改造体积的形状及尺寸，提高储层改造体积估计的精确度。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种确定储层改造体积的装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数；
根据所述裂缝参数获取裂缝的轮廓坐标；
整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种确定储层改造体积的装置，可以通过对地面倾角数据进行处理获得裂缝参数数据，根据裂缝参数构建凸包，从而利用凸包表征储层改造体积。利用本申请实施例提供的方案，定量的确定储层改造体积的形状及尺寸，提高储层改造体积估计的精确度。
本说明书还提供一种确定储层改造体积的系统，所述系统可以为单独的确定储层改造体积的系统，也可以应用在多种类型的压裂设计或者评价系统中。所述的系统可以为单独的服务器，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的服务器集群、系统(包括分布式系统)、软件(应用)、实际操作装置、逻辑门电路装置、量子计算机等并结合必要的实施硬件的终端装置。所述确定储层改造体积的系统可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意实施例所述方法的步骤。
需要说明的，上述所述的系统根据方法或者装置实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种确定储层改造体积的系统，可以通过对地面倾角数据进行处理获得裂缝参数数据，根据裂缝参数构建凸包，从而利用凸包表征储层改造体积。利用本申请实施例提供的方案，定量的确定储层改造体积的形状及尺寸，提高储层改造体积估计的精确度。
需要说明的是，本说明书上述所述的装置或者系统根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类、存储介质+程序实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (6)
1.一种确定储层改造体积的方法，其特征在于，包括：
获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数，其中，所述裂缝参数包括裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标；
根据所述裂缝参数获取裂缝的轮廓坐标，包括：将裂缝等效为椭圆面，以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零；将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；其中，所述旋转以逆时针旋转为正；
整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
2.根据权利要求1所述的确定储层改造体积的方法，其特征在于，所述对所述地面倾角数据进行处理获得裂缝参数，包括：
对所述地面倾角数据进行反演获得裂缝参数。
3.根据权利要求2所述的确定储层改造体积的方法，其特征在于，所述对所述地面倾角数据进行反演获得裂缝参数，包括：
根据预先建立的裂缝模型模拟获得理论地面倾角数据，对理论地面倾角数据与实际采集的地面倾角数据进行拟合，根据拟合结果确定裂缝参数。
4.一种确定储层改造体积的装置，其特征在于，包括：
裂缝参数获取模块，用于获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数，其中，所述裂缝参数包括裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标；
裂缝坐标获取模块，用于根据所述裂缝参数获取裂缝的轮廓坐标，包括：将裂缝等效为椭圆面，以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零；将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；其中，所述旋转以逆时针旋转为正；
凸包构建模块，用于根据所述裂缝参数获取裂缝的轮廓坐标，整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
储层改造体积确定模块，用于获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
5.一种确定储层改造体积的装置，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取压裂裂缝形成的地面倾角数据，对所述地面倾角数据进行处理获得裂缝参数，其中，所述裂缝参数包括裂缝的方位角、倾角、缝长、缝高、以及裂缝起裂点的三维空间坐标；
根据所述裂缝参数获取裂缝的轮廓坐标，包括：将裂缝等效为椭圆面，以原点为椭圆的中心，以缝长为长轴，沿Y轴，以缝高为短轴，沿X轴，Z轴坐标为零；将所述椭圆面沿Y轴旋转对应的倾角获得第一椭圆面；将第一椭圆面沿Z轴旋转对应的方位角获得第二椭圆面；将所述第二椭圆面的中心点从原点平移至裂缝的起点，获得第三椭圆面；获取第三椭圆面轮廓线的三维坐标，将所述第三椭圆面轮廓线的三维坐标确定为裂缝的轮廓坐标；其中，所述旋转以逆时针旋转为正；
整合压裂井段裂缝的轮廓坐标获得三维空间点阵集合，计算所述三维空间点阵集合所构成的凸包；
获取所述凸包的几何参数，根据所述凸包的几何参数确定储层改造体积。
6.一种确定储层改造体积的系统，其特征在于，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1-3中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711083854.0A CN107832524B (zh) | 2017-11-07 | 2017-11-07 | 一种确定储层改造体积的方法、装置及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711083854.0A CN107832524B (zh) | 2017-11-07 | 2017-11-07 | 一种确定储层改造体积的方法、装置及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107832524A CN107832524A (zh) | 2018-03-23 |
| CN107832524B true CN107832524B (zh) | 2021-01-01 |
Family
ID=61653780
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711083854.0A Active CN107832524B (zh) | 2017-11-07 | 2017-11-07 | 一种确定储层改造体积的方法、装置及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107832524B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111553504B (zh) * | 2019-01-23 | 2022-06-03 | 中国石油化工股份有限公司 | 一种计算压裂改造油气藏体积的方法及系统 |
| CN113803043B (zh) * | 2020-06-17 | 2023-07-07 | 中国石油化工股份有限公司 | 用于页岩气储层的压裂改造体积预测方法及装置 |
| CN111859260B (zh) * | 2020-07-03 | 2024-07-16 | 中国科学技术大学 | 一种利用压裂施工数据获得srv区域的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
| CN106226813A (zh) * | 2016-09-08 | 2016-12-14 | 南京特雷西能源科技有限公司 | 基于微地震的压裂缝网重构方法及装置 |
| CN106408208A (zh) * | 2016-10-14 | 2017-02-15 | 东北石油大学 | 体积压裂改造效果评价方法 |
| CN106869896A (zh) * | 2017-04-26 | 2017-06-20 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140058713A1 (en) * | 2012-08-21 | 2014-02-27 | Schlumberger Technology Corporation | Seismic modeling system and method |
-
2017
- 2017-11-07 CN CN201711083854.0A patent/CN107832524B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
| CN106226813A (zh) * | 2016-09-08 | 2016-12-14 | 南京特雷西能源科技有限公司 | 基于微地震的压裂缝网重构方法及装置 |
| CN106408208A (zh) * | 2016-10-14 | 2017-02-15 | 东北石油大学 | 体积压裂改造效果评价方法 |
| CN106869896A (zh) * | 2017-04-26 | 2017-06-20 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| "储层改造体积预测模型的研究与应用";翁定为等;《石油钻探技术》;20160131;第44卷(第1期);第95-100页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107832524A (zh) | 2018-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shakiba et al. | Using embedded discrete fracture model (EDFM) and microseismic monitoring data to characterize the complex hydraulic fracture networks | |
| CN105550448B (zh) | 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置 | |
| CN107832524B (zh) | 一种确定储层改造体积的方法、装置及系统 | |
| US9792726B2 (en) | Determining and displaying reserve estimates for a reservoir | |
| US10853533B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
| EP2924471A1 (en) | Systems and methods for modeling fracture networks in reservoir volumes from microseismic events | |
| CN108254784B (zh) | 一种基于二维地震数据的断层识别方法、装置及系统 | |
| CN108828669B (zh) | 一种二维相交测线静校正处理方法、装置及系统 | |
| US20170323477A1 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
| Yangsheng et al. | Three-dimensional fractal distribution of the number of rock-mass fracture surfaces and its simulation technology | |
| CN106597542B (zh) | 一种储层特征参数的预测方法及装置 | |
| Huan et al. | A new statistical parameter for determining joint roughness coefficient (JRC) considering the shear direction and contribution of different protrusions | |
| Díaz-Fernández et al. | Computation of influence functions for automatic mining subsidence prediction | |
| Chang et al. | Strength anisotropy of jointed rock slope considering mining damage: a numerical analysis method | |
| CN109358364B (zh) | 一种建立地下暗河储集体地质模型的方法、装置及系统 | |
| CN109001804B (zh) | 一种基于三维地震数据确定有效应力的方法、装置及系统 | |
| He et al. | A procedure to detect, construct, quantify, numerically simulate and validate fracture networks in coal blocks | |
| Zheng et al. | Stability analysis method for rock blocks formed by curved fractures | |
| CN104121864A (zh) | 一种岩石错动裂缝面的间隙评价方法 | |
| Scott et al. | Using field data to constrain a numerical kinematic model for ridge-transform deformation in the Troodos ophiolite, Cyprus | |
| CN111815769B (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
| Ma et al. | A rational and realistic rock mass modelling strategy for the stability analysis of blocky rock mass | |
| CN108874735B (zh) | 一种确定沉积盆地古流体压力的方法及装置 | |
| CN106291673B (zh) | 基于横波双折射的裂缝属性因子提取方法和装置 | |
| CN110907987A (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 |