CN113803043B - 用于页岩气储层的压裂改造体积预测方法及装置 - Google Patents
用于页岩气储层的压裂改造体积预测方法及装置 Download PDFInfo
- Publication number
- CN113803043B CN113803043B CN202010551931.6A CN202010551931A CN113803043B CN 113803043 B CN113803043 B CN 113803043B CN 202010551931 A CN202010551931 A CN 202010551931A CN 113803043 B CN113803043 B CN 113803043B
- Authority
- CN
- China
- Prior art keywords
- well section
- fracturing
- similarity
- parameters
- sample well
- 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
- 230000009466 transformation Effects 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000010276 construction Methods 0.000 claims abstract description 56
- 238000012216 screening Methods 0.000 claims abstract description 31
- 238000012986 modification Methods 0.000 claims abstract description 15
- 230000004048 modification Effects 0.000 claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 230000035945 sensitivity Effects 0.000 claims abstract description 11
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000011002 quantification Methods 0.000 claims description 8
- 238000000556 factor analysis Methods 0.000 claims description 7
- 238000000691 measurement method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000004927 clay Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 24
- 238000011161 development Methods 0.000 abstract description 12
- 238000005516 engineering process Methods 0.000 abstract description 12
- 238000005457 optimization Methods 0.000 abstract description 5
- 238000012544 monitoring process Methods 0.000 description 15
- 238000012360 testing method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 239000000700 radioactive tracer Substances 0.000 description 5
- 238000005056 compaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000007906 compression Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 235000015110 jellies Nutrition 0.000 description 1
- 239000008274 jelly Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种用于页岩气储层的压裂改造体积预测方法,其包含:将区块内具有微地震检测数据的井段作为样本井段,收集样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入基础数据库;以基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从工程地质参数以及压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数;计算目标井段与基础数据库中的样本井段之间高敏感地质参数的相似度以及高敏感压裂施工参数的相似度,将筛选出与目标井段相似度最高的样本井段的压裂改造体积作为所述目标井段的压裂改造体积。本发明低成本、客观、快速评价压后压裂改造体积评价技术,为压后评价和开发方案优化提供指导。
Description
技术领域
本发明涉及油气田开发技术领域,具体地说,涉及一种用于页岩气储层的压裂改造体积预测方法及装置。
背景技术
在常规油气藏和致密砂岩气藏中,常用单层裂缝的半长和导流能力作为增产压裂的主要描述参数,但是在页岩气藏中这些参数并不足以较好的表述增产特征。页岩气井实施压裂改造措施后,需用有效的方法确定压裂作业效果,量化表征改造体积,于是工程技术人员提出改造体积这一参数,即产生的复杂裂缝体积的大小可以近似用一个三维的体积(改造体积,SRV)表征。改造体积是增产压裂设计中的一个重要参数,能够对分析验证压裂设计参数、反馈施工质量、反演地层物性参数和岩石力学参数、研究正常的裂缝延伸规律和评价压裂效果、改善增产作业效果以及气井产能等提供必要的技术手段。改造体积可以通过微地震图测得,其与泵入流体体积和井的生产动态有关。
目前对于改造体积的评价方法主要依靠微地震监测、压裂过程示踪剂监测、水平井压后试井评价技术、水平井压后生产动态分析技术。
微地震压裂监测技术就是通过观测、分析由压裂过程中岩石破裂或错断所产生的微小地震事件来监测地下状态的地球物理技术,该方法是目前公认的最科学有效的评价改造体积的方法,缺点在于成本高,无法在每一个口井上都开展监测。
压后返排示踪剂监测技术是将用于监测多级压裂返排液的示踪剂分层段加入到实施多级压裂措施的油井中,在压裂液返排过程中,定期对压裂返排液进行取样,对所取样品中的多种示踪元素络合物的含量同时进行检测,根据检测到的各层段的示踪元素络合物的含量,实现压裂改造体积的评价,该方法解释结果量化程度低、对改造体积的判定不确定程度大。
试井是油气藏动态描述及动态监测的重要手段之一,在基于页岩气的渗流模型的理论计算,通过理论技术和实际测试的井底压力拟合的方法确定改造体积。水平井压后生产动态分析技术和试井技术的理论基础一致,在单井或井组理论产量计算和实际生产动态的拟合过程中得到页岩气储层的改造体积。
试井和产能动态分析的方法由于目前在理论模型方面尚不成熟,因此,所给出的分析结果与实际地下的情况误差较大。
虽然针对页岩气储层改造体积的评价技术越来越受到重视,但目前仍然缺乏一种科学、有效、低成本的改造体积的评价方法。针对页岩气的开发和储层改造压后评价的技术需求,需要一种可快速评价页岩气储层水平井压裂后改造体积的评价方法,为页岩气压后评价和开发方案优化提供指导。
因此,本发明提供了一种用于页岩气储层的压裂改造体积预测方法及装置。
发明内容
为解决上述问题,本发明提供了一种用于页岩气储层的压裂改造体积预测方法,所述方法包含以下步骤:
步骤一:建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集所述样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入所述基础数据库;
步骤二:利用多因素分析方法,以所述基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从所述工程地质参数以及所述压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数;
步骤三:基于相似度计量方法,计算目标井段与所述基础数据库中的样本井段之间所述高敏感地质参数的相似度以及所述高敏感压裂施工参数的相似度,将筛选出与所述目标井段相似度最高的样本井段的压裂改造体积作为所述目标井段的压裂改造体积。
根据本发明的一个实施例,所述步骤一中具体包含以下步骤:
收集所述样本井段的微地震事件三维坐标数据,利用体积包络法对所述样本井段的压裂改造体积进行计算,得到所述样本井段的压裂改造体积。
根据本发明的一个实施例,通过三维凸包增量法或二维凸包算法计算得到所述样本井段的压裂改造体积。
根据本发明的一个实施例,所述步骤二中具体包含以下步骤:
选取所述样本井段的压裂改造体积作为参考序列,所述工程地质参数以及所述压裂施工参数分别作为第一比较序列组以及第二比较序列组;
对所述参考序列、所述第一比较序列组以及所述第二比较序列组进行无量纲化处理,得到标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组。
根据本发明的一个实施例,所述步骤二中具体包含以下步骤:
基于所述标准化参考序列、所述标准化第一比较序列组以及所述标准化第二比较序列组进行关联系数计算,得到关联系数组;
通过关联度计算公式结合所述关联系数组计算关联度,筛选得到所述高敏感工程地质参数以及所述高敏感压裂施工参数。
根据本发明的一个实施例,所述步骤三中具体包含以下步骤:
基于相似度定量方法依次类比所述目标井段与所述样本井段的高敏感工程地质参数的相似度,筛选得到与所述目标井段相似度大于预设值的相似样本井段。
根据本发明的一个实施例,所述步骤三中具体包含以下步骤:
对于所述相似样本井段,基于所述相似度定量方法依次类比所述目标井段与所述相似样本井段的高敏感压裂施工参数的相似度,筛选得到与所述目标井段相似度最高的的样本井段。
根据本发明的一个实施例,所述工程地质参数包含:目标层深度、层位、储层厚度、地层倾角、杨氏模量、泊松比、孔隙度、测井渗透率、垂向应力、最大水平主应力、最小水平主应力、水平主应力差异、水平主应力差异系数、长英质含量、碳酸盐岩含量、黏土含量、裂缝指数。
根据本发明的一个实施例,所述压裂施工参数包含:段间距、簇间距、射孔密度、射孔孔径、射孔簇数、净压力指数、破裂压力、延伸压力、停泵压力、滑溜水用量、线性胶用量、冻胶用量、排量、100目支撑剂用量、40/70目支撑剂用量、30/50目支撑剂用量、20/40目支撑剂用量。
根据本发明的另一个方面,还提供了一种用于页岩气储层的压裂改造体积预测装置,所述装置包含:
基础数据库模块,其用于建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集所述样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入所述基础数据库;
参数筛选模块,其用于利用多因素分析方法,以所述基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从所述工程地质参数以及所述压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数;
相似度比较模块,其用于基于相似度计量方法,计算目标井段与所述基础数据库中的样本井段之间所述高敏感地质参数的相似度以及所述高敏感压裂施工参数的相似度,将筛选出与所述目标井段相似度最高的样本井段的压裂改造体积作为所述目标井段的压裂改造体积。
本发明提供的用于页岩气储层的压裂改造体积预测方法及装置首次运用类比评价的方法预测压裂改造体积,针对具有微地震等监测结果的已压裂页岩气井建立地质工程参数、压裂施工参数与微地震监测数据等样本数据集,通过类比样本井段与目标井段之间的相似度来快速获取目标井段压裂后的改造体积;本发明为油气田开发现场提供一种低成本、客观、快速评价压后压裂改造体积评价技术,为压后评价和开发方案优化提供指导。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1显示了根据本发明的一个实施例的用于页岩气储层的压裂改造体积预测方法流程图;
图2显示了根据本发明的一个实施例的凸包搜索示意图;
图3显示了根据本发明的一个实施例的基于凸包搜索的压裂改造体积搜索效果示意图;
图4显示了根据本发明的一个实施例的高敏感工程地质参数以及高敏感压裂施工参数筛选流程图;
图5显示了根据本发明的一个实施例的相似度最高的样本井段筛选流程图;以及
图6显示了根据本发明的一个实施例的用于页岩气储层的压裂改造体积预测装置结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下结合附图对本发明实施例作进一步地详细说明。
图1显示了根据本发明的一个实施例的用于页岩气储层的压裂改造体积预测方法流程图。
本发明涉及到非常规油气储层直井、斜井或水平井改造体积的压裂改造体积评价方法。
如图1,在步骤S101中,建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入基础数据库。
具体来说,工程地质参数包括:目标层深度(m)、层位、储层厚度(m)、地层倾角(°)、杨氏模量(GPa)、泊松比(无量纲)、孔隙度(%)、测井渗透率(mD)、垂向应力(MPa)、最大水平主应力(MPa)、最小水平主应力(MPa)、水平主应力差异(MPa)、水平主应力差异系数(无量纲)、长英质含量(%)、碳酸盐岩含量(%)、黏土含量(%)、裂缝指数(无量纲)等。
压裂施工参数包括:段间距(m)、簇间距(m)、射孔密度(孔/m)、射孔孔径(mm)、射孔簇数(簇)、净压力指数(无量纲)、破裂压力(MPa)、延伸压力(MPa)、停泵压力(MPa)、滑溜水用量(m3)、线性胶用量(m3)、冻胶用量(m3)、排量(m3/min)、100目支撑剂用量、40/70目支撑剂用量、30/50目支撑剂用量、20/40目支撑剂用量等。
具体来说,收集样本井段的微地震事件三维坐标数据,利用体积包络法对样本井段的压裂改造体积进行计算,得到样本井段的压裂改造体积。微地震事件三维坐标数据即采用微地震压裂监测技术得到的三维坐标数据。
进一步地,通过三维凸包增量法或二维凸包算法计算得到样本井段的压裂改造体积。在进行样本井段的压裂改造体积计算时,首先采用三维凸包增量法在三维坐标数据中构建一个四面体,如果无法构建四面体,则转而采用二维凸包算法计算样本井段的压裂改造体积。
图2显示了根据本发明的一个实施例的凸包搜索示意图。具体来说,利用三维凸包增量算法求出三维凸包的每个面,再用点到面的距离,暴力找出最小的距离。
采用三维凸包增量法计算样本井段的压裂改造体积时:
首先,在三维坐标数据中寻找四个点,构建一个四面体,可以先找两个不同点P1,P2,寻找和它们不共线的第三个点P3,再找不共面的第四个点P4。如果找不到即无法构建一个四面体,则调用二维凸包算法计算样本井段的压裂改造体积。
接下来计算除P1、P2、P3、P4外剩下点的随机排列。每次加一个点,有两种情况:
情况1:新点在当前凸包内部,只需简单地忽略该点,如图2中a所示。
情况2:新点在当前凸包外部,需要计算并的凸包,在这种情况下,首先需要计算原凸包相对于Pr的水平面,即Pr可以看到的封闭区域,如图2中b所示。
将P点能看到的所有平面删去,同时求出相对于Pr的水平面,将水平面上的边与Pr相连组成面,如图2中c所示。
之后,判断一个点P在凸包内还是凸包外,可利用有向体积的方法。在存储面时,保证面的法线方向朝向凸包外部,若存在某一平面和点P所组成的四面体的有向体积为正,则P点在凸包外部,并且此面可被P点看见。此时,只要将此面删去并用同样的方法判断与它相邻的其他平面是否可被P点看见,可以使用深度优先搜索(如图3所示)。
如图1,在步骤S102中,利用多因素分析方法,以基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从工程地质参数以及压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数。
具体来说,采用灰色关联法等多因素分析方法,分析确定影响SRV的工程地质参数(脆性、应力差异系数、层理发育等),压裂施工参数(排量、液量、施工压力等),用于后期的类比分析。
图4显示了根据本发明的一个实施例的高敏感工程地质参数以及高敏感压裂施工参数筛选流程图。
如图4,在步骤S401中,选取样本井段的压裂改造体积作为参考序列,工程地质参数以及压裂施工参数分别作为第一比较序列组以及第二比较序列组。
具体来说,选取压裂改造体积作为参考序列,该序列是反映系统行为特征的数据序列,亦即因变量;选取所有工程地质参数和所有压裂工程参数作为比较序列(包含第一比较序列组以及第二比较序列组),该序列由影响系统行为的因素的数据组成,即自变量。设参考序列为{Y(i)|i=1,2,3,…,n},比较序列为{Xj(i)|j=1,2,3,…,m,i=1,2,3,…,n}。
如图4,在步骤S402中,对参考序列、第一比较序列组以及第二比较序列组进行无量纲化处理,得到标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组。
具体来说,由于各因素序列的数据量纲可能不同,因此彼此之间就不便于比较或者比较之后难以得出正确的结论,因而必须进行数据的无量纲化处理。首先根据式(1)求出参考序列和比较序列的平均值和标准差。
式(1)中:μ表示所求序列的平均值,σ表示所求序列的标准差。然后利用式(2)即原始数据减去相应的平均值再除以标准差,由此得到的新数据序列便是标准化序列,其量纲为1,均值为0,方差为1。
如图4,在步骤S403中,基于标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组进行关联系数计算,得到关联系数组。
具体来说,对数据进行无量纲化处理之后,就可以进行关联系数的计算,y(i)与xj(i)的关联系数为下式(3):
式中,εj(i)为比较序列与参考序列的关联系数,分别为j个比较序列在各期的绝对差值中的最小者和最大者。式中ρ称为分辨系数。ρ越小,分辨力越大,一般ρ的取值区间为(0,1),具体取值可视情况而定。当ρ≤0.5463时,分辨力最好,通常取ρ=0.5。
如图4,在步骤S404中,通过关联度计算公式结合关联系数组计算关联度,筛选得到高敏感工程地质参数以及高敏感压裂施工参数。
具体来说,由于比较序列与参考序列的关系反映在n个关联系数上,关联信息比较分散,不利于整体比较,因此采用各比较序列的关联系数的平均值作为比较序列和参考序列间关联程度的数值表示,关联度ri公式如下式(4):
如图1,在步骤S103中,基于相似度计量方法,计算目标井段与基础数据库中的样本井段之间高敏感地质参数的相似度以及高敏感压裂施工参数的相似度,将筛选出与目标井段相似度最高的样本井段的压裂改造体积作为目标井段的压裂改造体积。
具体来说,计算目标井段与样本数据集中的样本井段之间的高敏感工程地质参数的相似度,筛选出高敏感工程地质参数相似度大于预设值(80%)的样本。在筛选基础上,再进行高敏感压裂施工参数的相似度类比,选择相似度最高的样本井段压裂改造体积近似为目标井段的压裂改造体积。
图5显示了根据本发明的一个实施例的相似度最高的样本井段筛选流程图。
如图5所示,在步骤S501中,基于相似度定量方法依次类比目标井段与样本井段的高敏感工程地质参数的相似度,筛选得到与目标井段相似度大于预设值的相似样本井段。
如图5所示,在步骤S502中,对于相似样本井段,基于相似度定量方法依次类比目标井段与相似样本井段的高敏感压裂施工参数的相似度,筛选得到与目标井段相似度最高的的样本井段。
具体做法:根据相似性科学理论,对归一化处理后的参数采用相似度定量方法类比井段的参数评价值。数据之间的相似性有两种表现形式,一种是以数据之间的差异度来衡量;另一种则是以其相似度来衡量。相似性度量函数sim(X,Y),其表达式为:
其中,X=(x1,…,xd)和Y=(y1,…,yd)是d维空间中的两个向量,表示第i维的区间长度,即数据在第i维上的最大距离,目的是使相似性度量不仅依赖于X和Y,而且也依赖于全体数据;是一个很小的常数,可以取值为0.001,以保证算法不被0除。根据该函数的特点,只要X和Y在某些维上的数值比较接近,它们就会表现出一定的相似性,而且随着数值接近的维数越多,它们之间的相似性程度也就越高。
一般相似性度量函数不满足三角不等式,所以不能将原来基于距离的聚类算法照搬到基于相似度的聚类算法中。算法参考凝聚层次聚类算法的基本过程,采用自底向上方法进行聚类。为了比较簇与对象之间的相似度,给出相似度定义为:
定义1对于一个对象集D={Xi|1=<i<=N}及其对象两两之间的相似度s,簇C和对象y的相似度定义为s(C,y)=min{s(x,y)|x∈C}。
取最小值目的是当一个对象想归并到一个簇里面时就必须与这个簇里面的所有对象做一一比较,看看相似度是否满足大于或等于阈值,满足条件才能并归进去。
然后,假设有N个d维数据Xi(1=<i<=n),算法第一步运用sim函数计算高维对象两两之间的相似度s(Xi,Xj),形成一个N×N的相似度矩阵,显然该矩阵为对称矩阵。
得出相似度矩阵后,算法将所有对象排成一列,然后把第一个对象归入到一个新簇中,再从头到尾扫描其它对象,检查未聚类对象和新簇对象的相似度是否大于或等于阈值,如果是就把它归并进去,这样对未聚类对象扫描一次以后就会得到一个新簇,同时保证簇里面对象两两之间的相似度至少达到指定的阈值。接下来算法再把剩下的未聚类对象排成一列,按照上述方法进行聚类,直到所有对象都聚类为止。
在一个实施例中,以南方某页岩区块为选定区块进行实验:
(1)建立南方某页岩气区块基础数据库,将具有微地震监测数据的56个井段作为样本井段,收集该批井段的工程地质参数和压裂施工参数基础数据,总计可参与类比的参数为67项;
(2)利用体积包络法对微地震监测数据进行解释,获取压裂后各井段的压裂改造体积,输入基础数据库;
(3)以压裂改造体积为目标,开展样本数据敏感性参数分析,确定主要影响该页岩气区块压后改造体积的高敏感工程地质参数及高敏感压裂施工参数各10项;
(4)以新钻水平井的某一段为目标井段,基于相似度计量方法计算目标井段与样本数据集中的样本井段之间的高敏感工程地质参数的相似度,筛选出高敏感工程地质参数相似度大于80%的样本,结果为33个样本类似;
(5)在筛选基础上,再进行高敏感压裂施工参数的相似度类比,类比得到与其相似度大于80%的段有3段,其中相似度最高为86%。相似度最高的样本井段改造体积近似为目标井段改造体积。
总结来说,本发明属于油气田勘探开发中的水力压裂技术领域,具体而言,涉及一种基于类比评价的压裂改造体积预测方法,包括以下几个步骤:
A、收集类比基础参数,包括工程地质参数与压裂井段的压裂施工参数;
B、收集微地震井解释井段事件三维坐标数据,利用包络法对样本井段SRV体积进行估算;
C、建立微地震井解释井段关键参数的样本数据集;
D、对目标井段和样本井段的参数进行标准化处理;
E、通过体积敏感性参数分析,确定影响某开发区块压后改造体积的高敏感工程地质及高敏感压裂施工参数各10项;
F、基于相似度计量方法,计算目标井段与样本数据集中的样本井段之间的相似度;
G、样本井段与目标井段的定量类比评价值计算;
H、通过参数集类比计算,筛选高相似度的样本井段的压裂改造体积作为目标井段的压裂改造体积,从而获取目标井段的改造体积。
本发明可以为现场提供一种低成本、客观、快速评价压后改造体积评价技术,为压后评价和开发方案优化提供指导。
图6显示了根据本发明的一个实施例的用于页岩气储层的压裂改造体积预测装置结构框图。
如图6,压裂改造体积预测装置600包含基础数据库模块601、参数筛选模块602以及相似度比较模块603。其中,基础数据库模块601包含三维凸包单元6011以及二维凸包单元6012;参数筛选模块602包含序列生成单元6021、标准化单元6022、关联系数单元6023以及关联度单元6024;相似度比较模块603包含第一筛选单元6031以及第二筛选单元6032。
基础数据库模块601用于建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入基础数据库。
具体来说,三维凸包单元6011以及二维凸包单元6012分别通过三维凸包增量法以及二维凸包算法计算得到样本井段的压裂改造体积。
参数筛选模块602用于利用多因素分析方法,以基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从工程地质参数以及压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数。
具体来说,序列生成单元6021用于选取样本井段的压裂改造体积作为参考序列,工程地质参数以及压裂施工参数分别作为第一比较序列组以及第二比较序列组。
具体来说,标准化单元6022用于对参考序列、第一比较序列组以及第二比较序列组进行无量纲化处理,得到标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组。
具体来说,关联系数单元6023用于基于标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组进行关联系数计算,得到关联系数组。
具体来说,关联度单元6024用于通过关联度计算公式结合关联系数组计算关联度,筛选得到高敏感工程地质参数以及高敏感压裂施工参数。
相似度比较模块603用于基于相似度计量方法,计算目标井段与基础数据库中的样本井段之间高敏感地质参数的相似度以及高敏感压裂施工参数的相似度,将筛选出与目标井段相似度最高的样本井段的压裂改造体积作为目标井段的压裂改造体积。
具体来说,第一筛选单元6031用于基于相似度定量方法依次类比目标井段与样本井段的高敏感工程地质参数的相似度,筛选得到与目标井段相似度大于预设值的相似样本井段。
具体来说,第二筛选单元6032对于相似样本井段,基于相似度定量方法依次类比目标井段与相似样本井段的高敏感压裂施工参数的相似度,筛选得到与目标井段相似度最高的的样本井段。
综上,本发明提供的用于页岩气储层的压裂改造体积预测方法及装置首次运用类比评价的方法预测压裂改造体积,针对具有微地震等监测结果的已压裂页岩气井建立地质工程参数、压裂施工参数与微地震监测数据等样本数据集,通过类比样本井段与目标井段之间的相似度来快速获取目标井段压裂后的改造体积;本发明为油气田开发现场提供一种低成本、客观、快速评价压后压裂改造体积评价技术,为压后评价和开发方案优化提供指导。
应该理解的是,本发明所公开的实施例不限于这里所公开的特定结构、处理步骤或材料,而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是,在此使用的术语仅用于描述特定实施例的目的,而并不意味着限制。
说明书中提到的“一个实施例”或“实施例”意指结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中。因此,说明书通篇各个地方出现的短语“一个实施例”或“实施例”并不一定均指同一个实施例。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (6)
1.一种用于页岩气储层的压裂改造体积预测方法,其特征在于,所述方法包含以下步骤:
步骤一:建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集所述样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入所述基础数据库;
步骤二:利用多因素分析方法,以所述基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从所述工程地质参数以及所述压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数;
步骤三:基于相似度计量方法,计算目标井段与所述基础数据库中的样本井段之间所述高敏感工程地质参数的相似度以及所述高敏感压裂施工参数的相似度,将筛选出与所述目标井段相似度最高的样本井段的压裂改造体积作为所述目标井段的压裂改造体积;
所述步骤二中具体包含以下步骤:选取所述样本井段的压裂改造体积作为参考序列,所述工程地质参数以及所述压裂施工参数分别作为第一比较序列组以及第二比较序列组;对所述参考序列、所述第一比较序列组以及所述第二比较序列组进行无量纲化处理,得到标准化参考序列、标准化第一比较序列组以及标准化第二比较序列组;
所述步骤二中具体包含以下步骤:基于所述标准化参考序列、所述标准化第一比较序列组以及所述标准化第二比较序列组进行关联系数计算,得到关联系数组;通过关联度计算公式结合所述关联系数组计算关联度,筛选得到所述高敏感工程地质参数以及所述高敏感压裂施工参数;
所述步骤三中具体包含以下步骤:基于相似度定量方法依次类比所述目标井段与所述样本井段的高敏感工程地质参数的相似度,筛选得到与所述目标井段相似度大于预设值的相似样本井段;
所述步骤三中具体包含以下步骤:对于所述相似样本井段,基于所述相似度定量方法依次类比所述目标井段与所述相似样本井段的高敏感压裂施工参数的相似度,筛选得到与所述目标井段相似度最高的样本井段。
2.如权利要求1所述的方法,其特征在于,所述步骤一中具体包含以下步骤:
收集所述样本井段的微地震事件三维坐标数据,利用体积包络法对所述样本井段的压裂改造体积进行计算,得到所述样本井段的压裂改造体积。
3.如权利要求2所述的方法,其特征在于,通过三维凸包增量法或二维凸包算法计算得到所述样本井段的压裂改造体积。
4.如权利要求1所述的方法,其特征在于,所述工程地质参数包含:目标层深度、层位、储层厚度、地层倾角、杨氏模量、泊松比、孔隙度、测井渗透率、垂向应力、最大水平主应力、最小水平主应力、水平主应力差异、水平主应力差异系数、长英质含量、碳酸盐岩含量、黏土含量、裂缝指数。
5.如权利要求1所述的方法,其特征在于,所述压裂施工参数包含:段间距、簇间距、射孔密度、射孔孔径、射孔簇数、净压力指数、破裂压力、延伸压力、停泵压力、滑溜水用量、线性胶用量、冻胶用量、排量、100目支撑剂用量、40/70目支撑剂用量、30/50目支撑剂用量、20/40目支撑剂用量。
6.一种用于页岩气储层的压裂改造体积预测装置,其特征在于,执行如权利要求1-5中任一项所述的方法,所述装置包含:
基础数据库模块,其用于建立目标井段所在区块的基础数据库,将区块内具有微地震检测数据的井段作为样本井段,收集所述样本井段的工程地质参数、压裂施工参数以及压裂改造体积录入所述基础数据库;
参数筛选模块,其用于利用多因素分析方法,以所述基础数据库中的压裂改造体积为目标,开展敏感性参数分析,从所述工程地质参数以及所述压裂施工参数中筛选出影响压裂改造体积的高敏感工程地质参数以及高敏感压裂施工参数;
相似度比较模块,其用于基于相似度计量方法,计算目标井段与所述基础数据库中的样本井段之间所述高敏感工程地质参数的相似度以及所述高敏感压裂施工参数的相似度,将筛选出与所述目标井段相似度最高的样本井段的压裂改造体积作为所述目标井段的压裂改造体积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010551931.6A CN113803043B (zh) | 2020-06-17 | 2020-06-17 | 用于页岩气储层的压裂改造体积预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010551931.6A CN113803043B (zh) | 2020-06-17 | 2020-06-17 | 用于页岩气储层的压裂改造体积预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113803043A CN113803043A (zh) | 2021-12-17 |
CN113803043B true CN113803043B (zh) | 2023-07-07 |
Family
ID=78943225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010551931.6A Active CN113803043B (zh) | 2020-06-17 | 2020-06-17 | 用于页岩气储层的压裂改造体积预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113803043B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104018831A (zh) * | 2014-06-24 | 2014-09-03 | 西南石油大学 | 一种压裂井储层评价方法 |
CN107025384A (zh) * | 2015-10-15 | 2017-08-08 | 赵乐平 | 一种复杂数据预测模型的构建方法 |
CN107742020A (zh) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层压裂改造体积预测方法 |
CN107832524A (zh) * | 2017-11-07 | 2018-03-23 | 中国石油天然气股份有限公司 | 一种确定储层改造体积的方法、装置及系统 |
CN108009716A (zh) * | 2017-11-28 | 2018-05-08 | 西南石油大学 | 一种水平井体积压裂效果影响因素多层次评价方法 |
CN108446797A (zh) * | 2018-03-06 | 2018-08-24 | 西南石油大学 | 一种致密油藏水平井体积压裂初期产能预测方法 |
CN108960651A (zh) * | 2018-07-11 | 2018-12-07 | 西南石油大学 | 一种致密油气藏多级压裂水平井完井效率的综合评价方法 |
CN109859245A (zh) * | 2019-01-22 | 2019-06-07 | 深圳大学 | 视频目标的多目标跟踪方法、装置及存储介质 |
CN110322363A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 页岩气储层改造体积计算方法及系统 |
-
2020
- 2020-06-17 CN CN202010551931.6A patent/CN113803043B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104018831A (zh) * | 2014-06-24 | 2014-09-03 | 西南石油大学 | 一种压裂井储层评价方法 |
CN107025384A (zh) * | 2015-10-15 | 2017-08-08 | 赵乐平 | 一种复杂数据预测模型的构建方法 |
CN107742020A (zh) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层压裂改造体积预测方法 |
CN107832524A (zh) * | 2017-11-07 | 2018-03-23 | 中国石油天然气股份有限公司 | 一种确定储层改造体积的方法、装置及系统 |
CN108009716A (zh) * | 2017-11-28 | 2018-05-08 | 西南石油大学 | 一种水平井体积压裂效果影响因素多层次评价方法 |
CN108446797A (zh) * | 2018-03-06 | 2018-08-24 | 西南石油大学 | 一种致密油藏水平井体积压裂初期产能预测方法 |
CN110322363A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 页岩气储层改造体积计算方法及系统 |
CN108960651A (zh) * | 2018-07-11 | 2018-12-07 | 西南石油大学 | 一种致密油气藏多级压裂水平井完井效率的综合评价方法 |
CN109859245A (zh) * | 2019-01-22 | 2019-06-07 | 深圳大学 | 视频目标的多目标跟踪方法、装置及存储介质 |
Non-Patent Citations (4)
Title |
---|
国外页岩气资源储量评价方法分析;张金华等;《中外能源》(第09期);全文 * |
基于GA-BP神经网络的长宁地区页岩气水平井产能预测技术;陈娟等;《科学技术与工程》(第05期);全文 * |
基于微地震事件点的SRV拟合方法比较研究;邵媛媛等;《西南石油大学学报(自然科学版)》(第04期);全文 * |
微地震监测技术在川东南页岩气压裂中的应用;王瑜;《中国石油和化工标准与质量》(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113803043A (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10607043B2 (en) | Subsurface reservoir model with 3D natural fractures prediction | |
CA2725923C (en) | Heterogeneous earth models for a reservoir field | |
US20130297272A1 (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
US20130046524A1 (en) | Method for modeling a reservoir basin | |
CN103026202A (zh) | 获取多孔介质的相容和综合物理性质的方法 | |
US9176255B2 (en) | Permeability prediction systems and methods using quadratic discriminant analysis | |
EP2279437A1 (en) | Continuous measurement of heterogeneity of geomaterials | |
US20140297186A1 (en) | Rock Classification Based on Texture and Composition | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
CN113803043B (zh) | 用于页岩气储层的压裂改造体积预测方法及装置 | |
CN110322363A (zh) | 页岩气储层改造体积计算方法及系统 | |
CN115905917A (zh) | 综合静动态数据构建海域低渗气层分类曲线的方法 | |
Sun et al. | Evaluation of Hydraulic Conductivity Estimates from Various Approaches with Groundwater Flow Models | |
CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
CN112983394A (zh) | 基于测井资料的曲线构建方法、装置和存储介质 | |
CN113075748A (zh) | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 | |
Tavoosi Iraj et al. | Integrated Petrophysical and heterogeneity assessment of the karstified fahliyan formation in the abadan plain, Iran | |
CN114510808A (zh) | 一种循缝找洞酸压数值模拟缝洞属性参数精细表征方法 | |
US11953647B2 (en) | System and method for radioactivity prediction to evaluate formation productivity | |
CN114624784B (zh) | 灰岩储层的状态确定方法、装置及可读存储介质 | |
Deng et al. | An Automated Data-Driven Workflow for Identifying Fractured Horizontal Well Sweet Spots in Shale Reservoirs | |
US20240036022A1 (en) | Workflow to predict source rock richness and net thickness using integrated inorganic, pyrolysis, and wireline data | |
CN116050656A (zh) | 基于主成分分析与神经网络的页岩脆性指数预测方法 | |
Ma et al. | A Novel Semi-Supervised Clustering Method for Carbonate Pore Size Distribution in Middle East Using Earth Mover Distance | |
CN116906033A (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 |