CN110378004A - 一种用于微地震解释压裂裂缝参数结果的校正方法 - Google Patents
一种用于微地震解释压裂裂缝参数结果的校正方法 Download PDFInfo
- Publication number
- CN110378004A CN110378004A CN201910626869.XA CN201910626869A CN110378004A CN 110378004 A CN110378004 A CN 110378004A CN 201910626869 A CN201910626869 A CN 201910626869A CN 110378004 A CN110378004 A CN 110378004A
- Authority
- CN
- China
- Prior art keywords
- microseism
- fracture
- well
- feeding amount
- fractured
- 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
Links
- 238000012937 correction Methods 0.000 claims abstract description 82
- 239000004576 sand Substances 0.000 claims abstract description 68
- 238000012544 monitoring process Methods 0.000 claims abstract description 27
- 238000010276 construction Methods 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 102000003712 Complement factor B Human genes 0.000 claims abstract 2
- 108090000056 Complement factor B Proteins 0.000 claims abstract 2
- 239000011159 matrix material Substances 0.000 claims description 37
- 239000012530 fluid Substances 0.000 claims description 13
- 239000007788 liquid Substances 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 208000010392 Bone Fractures Diseases 0.000 description 112
- 206010017076 Fracture Diseases 0.000 description 112
- 238000010586 diagram Methods 0.000 description 17
- 238000000034 method Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000513 principal component analysis Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2431—Multiple classes
-
- 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 OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Data Mining & Analysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Emergency Management (AREA)
- Business, Economics & Management (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种用于微地震解释压裂裂缝参数结果的校正方法,所述方法包括:A、根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类;B、对分类后的各所述压裂井单井进行体积系数校正,各所述压裂井的体积校正系数BV为:其中,[L5/4H]W为利用微地震监测解释结果的计算值;[L5/4H]T为利用现场施工参数得到的理论计算值;C、根据计算得到的各所述压裂井的单井体积校正系数,计算各类压裂井的类体积校正系数。本发明结合现场的施工参数能够对微地震解释结果进行校正,从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态,进而可以辅助微地震位置精确反演,有效指导水力压裂施工作业。
Description
技术领域
本发明涉及微地震监测解释技术领域,具体涉及一种用于微地震解释压裂裂缝参数结果的校正方法。
背景技术
随着矿藏开采技术的不断开发,水力裂缝遇到天然裂缝的行为受到越来越多人的关注。水力压裂是公认的在致密裂隙储层中提高采收率的主要增产技术。水力压裂产生的水力裂缝可激活天然裂缝,形成复杂的弥散性裂缝缝网,起到增渗、增产的效果,因此对水力裂缝与天然裂缝相交模式及作用效果的研究具有极高的应用价值。
微地震监测是解释压裂时人工裂缝(即水力裂缝)缝长、缝高及方位的有效方法之一。水力压裂微地震裂缝监测技术可以用于指示裂缝位置、分析裂缝发育情况,并辅助微地震位置精确反演,指导水力压裂施工作业。水力压裂微地震裂缝监测技术是通过监测水力压裂过程中产生的微地震事件,来解释压后人工裂缝形态的参数。
目前研究重点大多是在改进微地震重构理论算法上,微地震监测过程中由于噪点影响导致微地震解释裂缝参数失真,微地震解释裂缝形态与现场压裂过程中的加砂量和液量时常不具备对应关系,导致微地震反演出的裂缝参数不具备参考价值,严重影响了对压后裂缝形态以及压裂效果的评价。
因此有必要发明一种可以用于微地震解释压裂裂缝参数结果的校正方法,以能够更准确真实地推导出裂缝参数,指示裂缝位置、分析裂缝发育情况,辅助微地震位置精确反演,指导水力压裂施工作业。
发明内容
高排量压裂施工时会产生大量微地震信号事件点,利用微地震监测数据难以真实解释出压裂后的有效裂缝形态,因此有必要结合现场的施工参数对对微地震解释结果进行校正。本发明旨在提供一种用于微地震解释压裂裂缝参数结果的校正方法。
本发明提供一种用于微地震解释压裂裂缝参数结果的校正方法,所述方法包括:
A、根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井(部分取样,作为学习样本)的沟通模式进行分类;
B、对分类后的各所述压裂井单井进行体积系数校正,各所述压裂井的体积校正系数BV为:
其中,[L5/4H]W为利用微地震监测解释结果的计算值;[L5/4H]T为利用现场施工参数得到的理论计算值,而Vi=Vfp+Vfs,
上式中L为裂缝半长,单位为m;H为裂缝高度,单位为m;Ef为总造缝效率;Vi为有效造缝液量,单位为m3;Vfp为前置液量,单位为m3;Vfs为携砂液量,单位为m3;μ为压裂液粘度,单位为Pa.s;Q为排量,单位为m3/min;υ为岩石的泊松比;E为弹性模量,单位是GPa;当Q取地面总排量时,α=1.26;当Q取地面总排量之半时,α=1.5;
C、根据计算得到的各所述压裂井的单井体积校正系数,计算各类压裂井的类体积校正系数,各类的类体积校正系数取该类压裂井中单井体积校正系数计算结果的平均值。
上述的用于微地震解释压裂裂缝参数结果的校正方法,根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类的标准为:当微地震解释的水力裂缝长度(双翼缝长之和)与加砂量的比值<10时,将微地震解释的水力裂缝与对应的加砂量作为分类的依据;当微地震解释的水力裂缝长度与加砂量的比值>=10时,将加砂量作为分类的依据。
上述的用于微地震解释压裂裂缝参数结果的校正方法,当将微地震解释的水力裂缝与对应的加砂量作为分类的依据时,利用数理分析计算微地震解释的水力裂缝长度与加砂量的主成分,将该主成分值作为分类的依据,微地震解释裂缝长度与加砂量的主成分得分越高表示施工加砂量及所述裂缝长度越大,反之亦然。
上述的用于微地震解释压裂裂缝参数结果的校正方法,当所述主成分值<130时,将所述压裂井分类为第一类;当所述主成分值为130-160时,将所述压裂井分类为第二类第I型;当所述主成分值为160-200,将所述压裂井分类为第二类第II型;当所述主成分值为200-250时,将所述压裂井分类为第三类第I型;当所述主成分值>250时,将所述压裂井分类为第三类第II型。
上述的用于微地震解释压裂裂缝参数结果的校正方法,进行主成分计算时,首先对待校正区块部分样本压裂井的加砂量以及微地震解释裂缝长度进行统计,得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵,随后对观测矩阵进行计算得到均值向量;之后计算得到该区块二维随机变量样本的A矩阵;再用样本的无偏估计来代替总体协方差S*;再计算得到区块分类相关系数矩阵的特征向量;如果随后计算得到的第一主成分贡献率不能满足要求,则添加第二主成分;如果随后计算得到的第一主成分贡献率满足要求,则以该第一主成分值作为依据进行分类。
上述的用于微地震解释压裂裂缝参数结果的校正方法,如果所述协方差的矩阵数值相差较大,在求解所述区块的取样压裂井的主成分的时候应基于相关系数矩阵求解,根据上述的样本协方差矩阵,求解得到所述区块的所述二维随机变量的相关系数矩阵;再根据所述相关系数矩阵,计算得到所述区块的相关系数矩阵的所述特征向量。
上述的用于微地震解释压裂裂缝参数结果的校正方法,当所述加砂量<10m3时,将所述压裂井分类为第一类;当所述加砂量为10-13m3时,将所述压裂井分类为第二类第I型;当所述加砂量为13-16m3时,将所述压裂井分类为第二类第II型;当所述加砂量为16-20m3时,将所述压裂井分类为第三类第I型;当所述加砂量>20m3时,将所述压裂井分类为第三类第II型。
上述的用于微地震解释压裂裂缝参数结果的校正方法,根据各类所述压裂井的类体积校正系数得到各类所述压裂井的类裂缝半缝长校正系数BL,
BL=a5/4
其中, 为所述压裂井的类体积校正系数。
上述的用于微地震解释压裂裂缝参数结果的校正方法,根据所述体积校正系数得到缝高校正系数BH:
BH=a
其中, 为所述压裂井的类体积校正系数。
上述的用于微地震解释压裂裂缝参数结果的校正方法,所述压裂井的裂缝为对称的双翼垂直裂缝,裂缝井的水平截面形状为椭圆形。
本发明结合现场的施工参数能够对微地震解释结果进行校正,从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态,进而可以辅助微地震位置精确反演,有效指导水力压裂施工作业。
附图说明
图1a为本发明第一类压裂井的示意图;
图1b为本发明第二类压裂井的示意图;
图1c为本发明第三类压裂井的示意图;
图2为本发明压裂井分类流程示意图;
图3为本发明一实施例压裂井分类输出结果示意图;
图4为本发明裂缝PKN模型示意图;
图5a为本发明一实施例10-12井微地震裂缝回归结果示意图;
图5b为本发明一实施例10-12井微地震裂缝校正结果示意图;
图6a为本发明一实施例5-17井微地震裂缝回归结果示意图;
图6b为本发明一实施例5-17井微地震裂缝校正结果示意图;
图7a为本发明一实施例4-13井微地震裂缝回归结果示意图;
图7b为本发明一实施例4-13井微地震裂缝校正结果示意图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明,以便更清楚理解本发明的目的、特点和优点。应理解的是,附图所示的实施例并不是对本发明范围的限制,而只是为了说明本发明技术方案的实质精神。
高排量压裂施工时会产生大量微地震信号事件点,利用微地震监测数据难以真实解释出压裂后的有效裂缝形态,因此有必要结合现场的施工参数对对微地震解释结果进行校正。
本发明对裂缝的沟通模式进行了划分,本实施方式以青海油田36口压裂井为背景,对水力裂缝与天然裂缝的沟通模式进行定义并给出以裂缝沟通模式为依据对压裂井的分类方法;基于微地震监测结果以及现场施工参数(加砂量)对三类裂缝沟通模式的压裂井微地震解释结果的校正系数(体积校正系数、缝长校正系数和缝高校正系数)进行了理论推导和实际计算,为现场实际施工提出依据。
储层天然裂缝发育是出现复杂压后裂缝形态的客观条件,当进行压裂施工时,人工裂缝沟通天然裂缝会形成复杂裂缝系统,反之,会形成单一的两翼主裂缝,在进行天然裂缝与水力裂缝组合分类之前,需要对研究区块的压裂井段的人工裂缝发育情况进行表征,以评价研究区块是否具有能形成复杂缝网的可能性。因此需要分类表征压裂井段的人工裂缝与天然裂缝的沟通后的形态,为微地震监测有效信号识别以及后续校正系数推导提供依据。
根据现场施工情况,将压后人工裂缝与水力裂缝沟通模式分为三类,如图1a、图1b、图1c所示。图1a为第一类:天然裂缝不发育,水力裂缝主缝为单一横切缝没有沟通天然裂缝;图1b为第二类:天然裂缝较为发育,水力裂缝形态较为单一并且沟通部分天然裂缝;图1c为第三类:天然裂缝发育,水力裂缝为多组不同产状的裂缝,沟通天然裂缝,形成复杂组合形式。其中,图中的各参数代表的意义如下:
σv-垂向正应力;
σH-垂向正应力;
σh-垂向正应力;
HF-水力裂缝(HF1表示第一条水力裂缝,依次类推);
NF-天然裂缝(NF1表示第一条水力裂缝,依次类推)。
在进行水力裂缝沟通天然裂缝模式进行分类分析时,用到的现场定量数据主要是微地震监测裂缝长度解释结果以及现场的加砂量。处理过程可分为两类:
(1)当微地震监测结果与加砂量出现矛盾,即微地震解释长水力裂缝但现场加砂量少,此时可认为微地震解释结果可行度较低,主要以实际施工加砂量作为分类依据;
(2)当微地震监测结果与加砂量吻合程度较好,即微地震解释的长水力裂缝对应的加砂量大,反之亦然。此时由于两者相关性较强,可考虑利用数理分析的主成分分析理论将二者定义成一个新的主成分以减少多元变量维数,简化分类流程。
具体分类流程如图2所示,首先对某一区块部分样本压裂井的加砂量以及微地震解释裂缝长度进行统计,得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵,随后对观测矩阵进行计算得到均值向量为;之后根据公式计算得到该区块二维随机变量样本的A矩阵。
根据数理统计相关知识,用样本的无偏估计来代替总体协方差S*。
如果由样本的协方差矩阵发现,其数值相差较大,为防止大数吃小数的现象,在求解该区块样本的主成分的时候可以基于相关系数矩阵求解,根据上述的样本协方差矩阵,求解得到该区块样本分类过程中二维随机变量的相关系数矩阵ρ。
根据相关系数矩阵,计算得到该区块分类相关系数矩阵的特征向量λ1、λ2。
根据公式,计算得到第一主成分贡献率,看是否可以满足要求。
如果可以满足要求,即可得到第一主成分对应的特征向量为
故第一主成分可以表示为
F1=c1x1+c2x2
其中x1为加砂量,x2为微地震解释裂缝长度(双翼缝长之和)。
而如果计算得到第一主成分贡献率不能满足要求,则需要计算第二主成分,使得第一第二主成分的累计贡献率达到工程需求。
以上计算都基于现有数理统计方法,此处不再详细给出计算公式。
基于主成分分析的分类方法,本实施例以visual basic6.0为平台开发了归类的软件,该软件支持用户自定义的分类标准,输出结果主要包括样本压裂井归属的类比,分类基于的参数等,其运行界面如图3所示。
油田现场水力压裂前期设计常采用直观快捷的PKN模型,该模型假设裂缝空间分布为一均匀的椭球体,因而本发明的基本假设如下:
(1)支撑剂颗粒在高粘压裂液中呈悬浮状态,携砂液在裂缝中近似以柱塞状向前推进;
(2)砂比的大小对携砂液的滤失特性没有影响;
(3)裂缝为对称的双翼垂直裂缝,裂缝水平截面形状为椭圆面。
如图4所示,为了简化计算,我们将计算的裂缝模型考虑为PKN模型。
可按照水力裂缝沟通天然裂缝的复杂程度进行分类,确定各类的临界值。
根据现场施工需要,按照裂缝长度与加砂量比值落入的范围,区分现场实际的施工加砂量与微地震解释缝长是否能很好对应。具体为:
(1)基于主成分分析的分类标准
根据现场施工需要,认为当裂缝长度与加砂量比值大于等于一数值时比如10时,认为现场实际的施工加砂量与微地震解释缝能较好对应,微地震解释裂缝长度与加砂量共同作为分类标准。显然,微地震解释裂缝长度与加砂量的主成分得分越高表示施工加砂量以及裂缝长度越大,则形成天然裂缝水力裂缝沟通模式也更为复杂,结合现场分析需要给出基于主成分得分对沟通模式进行分类的分类标准。
(2)基于加砂量的分类标准
根据现场施工需要,认为当裂缝长度与加砂量比值小于一数值时比如10时,认为现场实际的施工加砂量与微地震解释缝长不能很好对应,微地震解释结果不具有参考价值,微地震解释裂缝长度不能作为分类标准,而以加砂量作为沟通模式分类的标准。
根据上述分类标准得到某区块待判样本井的归类结果。
再根据区块压裂井的分类结果,得到各类井的体积校正数据BV,类的体积校正系数取该类井中单井体积校正系数计算结果的平均值。再根据区块的体积校正系数计算结果,依据公式得到各类井的裂缝半缝长校正系数BL和缝高校正系数BH。裂缝半缝长又称裂缝半径。煤储层压裂强化后,裂缝从井筒沿径向向煤储层延伸的距离,一般指水平裂缝的缝长。裂缝半长是裂缝尺寸要素之一。
本发明各参数的校正推导如下。
设Vi为有效造缝液量,其与前置液量Vfp和携砂液Vfs有如下关系
Vi=Vfp+Vfs (1)
前置液量Vfp和携砂液Vfs根据现场压裂施工可以直接测得。
在这里引用诺尔特(K G.N ol t e)提出的泵注砂比连续变化公式中总造缝效率Ef的概念,其表示前置液与携砂液的总的造缝效率。
Ef根据上式中第二个等式以及现场压裂施工参数可以直接计算得到。
式中:
Vfrac-单侧裂缝体积,m3;
根据上述PKN模型,最大裂缝宽度Wmax可由下式计算得到。
式中:
Wmax-牛顿液体在层流条件下裂缝的最大缝宽,m;
μ-压裂液粘度,Pa.s;
Q-排量,m3/min;
L-裂缝半长,m
υ-岩石泊松比
E-弹性模量,GPa
当Q取地面总排量的时候,α=1.26;当Q取地面总排量之半的时候,α=1.5;在这里由于是取地面总排量,故α=1.26。Q的取值可根据个人计算需要而定。
推导校正系数的时候,可首先根据现场的携砂液量和前置液量,由公式(1)计算得到有效造缝液量Vi,随后根据公式(2)计算得到总造缝效率Ef,最后计算得到裂缝体积如公式(4)所示:
Vfrac=EfVi/2 (4)
另一方面
Vfrac=πWmaxLH (5)
式中:
Wmax-牛顿液体在层流条件下裂缝的最大缝宽,m;
L-裂缝半长,m;
H-裂缝高度,m;
上式中L、H为理论值。
将公式(5)与公式(3)联立,得到
联立公式(6)和公式(4)得到
将L5/4H与微地震监测结果做比值,即可得到体积校正系数Bv,即:
式中:
[L5/4H]W-利用微地震监测解释结果的计算值(L和H利用微地震解释结果得到);
[L5/4H]T-利用现场施工参数得到的理论计算值(依据7式的计算公式,根据现场实际施工参数得到计算值)。
将所得到的体积校正系数BV反演到微地震解释结果以得到缝长和缝高校正系数,推导过程如下
BV[L5/4H]W=[L5/4H]T (9)
另
(La)5/4Ha=[L5/4H]T (10)
a5/4=BL
a=BH
式中:
BL-裂缝半缝长校正系数;
BH-缝高校正系数。
则
a9/4=BV (11)
由公式(11)可以反解出中间变量a,校正系数从而求解出裂缝半缝长校正系数BL和缝高校正系数BH。
上述推导得到的校正系数为单井的校正系数,而为了实际描述各类压裂井的裂缝形态,需要得到各类的校正系数,得到各类的校正系数,才能够对裂缝形态进行校正。
计算各类压裂井的类体积校正系数,各类的类体积校正系数取该类压裂井中单井体积校正系数计算结果的平均值。则计算各类的缝宽与缝高校正系数时, 为所述压裂井的类体积校正系数。
本发明结合现场的施工参数能够对微地震解释结果进行校正,从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态,进而可以辅助微地震位置精确反演,有效指导水力压裂施工作业。
实施例
对切六区块36口压裂井的加砂量以及微地震解释裂缝长度进行统计,得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵,随后对观测矩阵进行计算得到均值向量为:
根据公式计算得到切六区块二维随机变量样本的A矩阵为:
根据数理统计相关知识,用样本的无偏估计来代替总体协方差。
由样本的协方差矩阵可以发现,其数值相差较大,为防止大数吃小数的现象,在求解切六区块的主成分的时候可以基于相关系数矩阵求解,根据上述的样本协方差矩阵,求解得到切六区块分类过程中二维随机变量的相关系数矩阵。
根据相关系数矩阵,计算得到了切六区块分类相关系数矩阵的特征向量。
λ1=1.58
λ2=0.42
根据公式,计算得到第一主成分贡献率为0.79,可以满足要求。
随后得到第一主成分对应的特征向量为
故第一主成分可以表示为
F1=0.5x1+0.87x2
x1为加砂量,x2为微地震解释裂缝长度(双翼缝长之和)。
由于切六区块压裂井较多,可按照水力裂缝沟通天然裂缝的复杂程度将二类及三类分别再细分成两个亚类,随后确定各类的临界值。
(1)基于主成分分析的分类标准
显然,主成分得分越高表示施工加砂量以及裂缝长度越大,则形成天然裂缝水力裂缝沟通模式也更为复杂,结合现场分析需要给出基于主成分得分对沟通模式进行分类的分类标准,分类标准如表1.1所示。
表1.1基于主成分得分的分类标准表
(2)基于加砂量的分类标准
根据现场施工需要,认为当裂缝长度与加砂量比值小于10时,认为现场实际的施工加砂量与微地震解释缝长不能很好对应,微地震解释结果不具有参考价值,微地震解释裂缝长度不能作为分类标准,而以加砂量作为沟通模式分类的标准,其标准如表1.2所示。
表1.2切六区块基于加砂量的分类标准表
将待判的36口井参数逐一输入所编制的软件,可得到切六区块待判的36口井的归类结果,如表1.3所示。
表1.3切六区块压裂井分类结果表
根据切六区块压裂井的分类结果,分别根据本发明的推导思路得到一类井、二类井和三类井及其亚类的体积校正数据BV,类的体积校正系数取该类井中单井体积校正系数计算结果的平均值。在这里,取切六区块岩石的弹性模量为336GPa,泊松比为0.2,裂缝半缝长取裂缝微地震监测结果两个半翼缝长的平均值。切六区块各类体积校正系数计算结果如表2.1所示。
表2.1切六区块体积校正系数计算结果表
根据切六区块的体积校正系数计算结果,依据公式(9)、(10)、(11)得到了各类井及其亚类的裂缝半缝长校正系数BL和缝高校正系数BH,其统计结果见表2.2。
表2.2切六区块压裂井各校正系数统计表
分别以切六区块10-12井,5-17井,4-13井作为三类典型井,利用微地震原始数据,基于八叉树理论算法对人工裂缝进行重构,并考虑校正系数对人工裂缝重构结果进行校正,三类井回归与校正结果分别如图5a的10-12井微地震裂缝回归结果示意图、5b的10-12井微地震裂缝校正结果示意图、图6a的5-17井微地震裂缝回归结果示意图、图6b的5-17井微地震裂缝校正结果示意图、图7a的4-13井微地震裂缝回归结果示意图、图7b的4-13井微地震裂缝校正结果示意图所示。
由以上对比分析看出,根据本发明的校正方法重构的裂缝与最初校正前的裂缝存在差异,不同类别的井裂缝校正之后的结果与最初的结果差异不同,一类井校正裂缝形态参数小于校正之前的裂缝参数,二类校正前后基本一致,三类井校正裂缝形态参数大于校正之前的裂缝参数。
上述各实施例仅用于说明本发明,其中实施例的各零部件、装置、方法步骤都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (10)
1.一种用于微地震解释压裂裂缝参数结果的校正方法,所述方法包括:
A、根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类;
B、对分类后的各所述压裂井单井进行体积系数校正,各所述压裂井的体积校正系数BV为:
其中,[L5/4H]W为利用微地震监测解释结果的计算值;[L5/4H]T为利用现场施工参数得到的理论计算值,而Vi=Vfp+Vfs,
上式中L为裂缝半长,单位为m;H为裂缝高度,单位为m;Ef为总造缝效率;Vi为有效造缝液量,单位为m3;Vfp为前置液量,单位为m3;Vfs为携砂液量,单位为m3;μ为压裂液粘度,单位为Pa.s;Q为排量,单位为m3/min;υ为岩石泊松比;E为弹性模量,单位为GPa;当Q取地面总排量时,α=1.26;当Q取地面总排量之半时,α=1.5;
C、根据计算得到的各所述压裂井的单井体积校正系数,计算各类压裂井的类体积校正系数,各类的类体积校正系数取该类压裂井中单井体积校正系数计算结果的平均值。
2.根据权利要求1所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类的标准为:当微地震解释的水力裂缝总的长度与加砂量的比值<10时,将微地震解释的水力裂缝与对应的加砂量作为分类的依据;当微地震解释的水力裂缝长度与加砂量的比值>=10时,将加砂量作为分类的依据。
3.根据权利要求2所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,当将微地震解释的水力裂缝与对应的加砂量作为分类的依据时,利用数理分析计算微地震解释的水力裂缝长度与加砂量的主成分,将该主成分值作为分类的依据,微地震解释裂缝长度与加砂量的主成分得分越高表示施工加砂量及所述裂缝长度越大,反之亦然。
4.根据权利要求3所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,当所述主成分值<130时,将所述压裂井分类为第一类;当所述主成分值为130-160时,将所述压裂井分类为第二类第I型;当所述主成分值为160-200,将所述压裂井分类为第二类第II型;当所述主成分值为200-250时,将所述压裂井分类为第三类第I型;当所述主成分值>250时,将所述压裂井分类为第三类第II型。
5.根据权利要求4所述的用于微地震解释压裂裂缝参数结果的校正方法,进行主成分计算时,首先对带校正区块部分样本压裂井的加砂量以及微地震解释裂缝长度进行统计,得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵,随后对观测矩阵进行计算得到均值向量;之后计算得到该区块二维随机变量样本的A矩阵;再用样本的无偏估计来代替总体协方差S*;再计算得到区块分类相关系数矩阵的特征向量;如果随后计算得到第一主成分贡献率不能满足要求,则添加第二主成分;如果随后计算得到第一主成分贡献率满足要求,则以该第一主成分值作为依据进行分类。
6.根据权利要求5所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,如果所述协方差的矩阵数值相差较大,在求解所述区块的取样压裂井的主成分的时候基于相关系数矩阵求解,根据上述的样本协方差矩阵,求解得到所述区块的所述二维随机变量的相关系数矩阵;再根据所述相关系数矩阵,计算得到所述区块的相关系数矩阵的所述特征向量。
7.根据权利要求2所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,当所述加砂量<10m3时,将所述压裂井分类为第一类;当所述加砂量为10-13m3时,将所述压裂井分类为第二类第I型;当所述加砂量为13-16m3时,将所述压裂井分类为第二类第II型;当所述加砂量为16-20m3时,将所述压裂井分类为第三类第I型;当所述加砂量>20m3时,将所述压裂井分类为第三类第II型。
8.根据权利要求1至5任一项所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,根据各类所述压裂井的类体积校正系数得到各类所述压裂井的类裂缝半缝长校正系数BL,BL=a5/4
其中, 为所述压裂井的类体积校正系数。
9.根据权利要求1至5任一项所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,根据所述体积校正系数得到缝高校正系数BH,
BH=a,其中, 为所述压裂井的类体积校正系数。
10.根据权利要求9所述的用于微地震解释压裂裂缝参数结果的校正方法,其特征在于,所述压裂井的裂缝为对称的双翼垂直裂缝,裂缝井的水平截面形状为椭圆形。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910626869.XA CN110378004B (zh) | 2019-07-11 | 2019-07-11 | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
PCT/CN2020/100548 WO2021004442A1 (zh) | 2019-07-11 | 2020-07-07 | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
US17/593,927 US20220206174A1 (en) | 2019-07-11 | 2020-07-07 | Correction method for micro seismic interpretation fracturing fracture parameter result |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910626869.XA CN110378004B (zh) | 2019-07-11 | 2019-07-11 | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110378004A true CN110378004A (zh) | 2019-10-25 |
CN110378004B CN110378004B (zh) | 2021-05-07 |
Family
ID=68252773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910626869.XA Active CN110378004B (zh) | 2019-07-11 | 2019-07-11 | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20220206174A1 (zh) |
CN (1) | CN110378004B (zh) |
WO (1) | WO2021004442A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110967748A (zh) * | 2019-12-02 | 2020-04-07 | 中国矿业大学(北京) | 微地震成像裂缝提取方法、装置、设备和存储介质 |
WO2021004442A1 (zh) * | 2019-07-11 | 2021-01-14 | 中国石油大学(北京) | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
CN113283197A (zh) * | 2021-06-10 | 2021-08-20 | 西南石油大学 | 一种基于复杂裂缝缝宽分布的加砂参数设计方法 |
CN114718536A (zh) * | 2021-01-05 | 2022-07-08 | 中国石油天然气股份有限公司 | 压裂工艺调整方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114542040B (zh) * | 2022-02-25 | 2024-04-26 | 山西蓝焰煤层气集团有限责任公司 | 一种连续脉冲水力压裂系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010117584A1 (en) * | 2009-04-07 | 2010-10-14 | Schlumberger Canada Limited | System and technique to quantify a fracture system |
WO2014105659A1 (en) * | 2012-12-27 | 2014-07-03 | Schlumberger Canada Limited | Method of calibrating fracture geometry to microseismic events |
CN104297783A (zh) * | 2014-06-30 | 2015-01-21 | 中国石油天然气集团公司 | 一种水力压裂微地震事件的解释方法以及系统 |
CN104459775A (zh) * | 2014-11-28 | 2015-03-25 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于微地震监测数据的页岩气藏裂缝建模方法 |
CN105629295A (zh) * | 2015-12-29 | 2016-06-01 | 四川圣诺油气工程技术服务有限公司 | 一种页岩气体积压裂微地震监测方法 |
CN105672973A (zh) * | 2016-01-26 | 2016-06-15 | 西南石油大学 | 一种煤层气藏整体体积压裂优化设计方法 |
CN106408208A (zh) * | 2016-10-14 | 2017-02-15 | 东北石油大学 | 体积压裂改造效果评价方法 |
CN108984839A (zh) * | 2018-06-15 | 2018-12-11 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种基于微地震的页岩改造体积快速判断方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2506039A3 (en) * | 2011-03-28 | 2013-08-14 | Conocophillips Company | Methods for Seismic Fracture Parameter Estimation and Gas Filled Fracture Identification From Vertical Well Log Data |
CN104730596B (zh) * | 2015-01-25 | 2017-11-07 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
CN110378004B (zh) * | 2019-07-11 | 2021-05-07 | 中国石油大学(北京) | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
-
2019
- 2019-07-11 CN CN201910626869.XA patent/CN110378004B/zh active Active
-
2020
- 2020-07-07 US US17/593,927 patent/US20220206174A1/en active Pending
- 2020-07-07 WO PCT/CN2020/100548 patent/WO2021004442A1/zh active Application Filing
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010117584A1 (en) * | 2009-04-07 | 2010-10-14 | Schlumberger Canada Limited | System and technique to quantify a fracture system |
WO2014105659A1 (en) * | 2012-12-27 | 2014-07-03 | Schlumberger Canada Limited | Method of calibrating fracture geometry to microseismic events |
CN104297783A (zh) * | 2014-06-30 | 2015-01-21 | 中国石油天然气集团公司 | 一种水力压裂微地震事件的解释方法以及系统 |
CN104459775A (zh) * | 2014-11-28 | 2015-03-25 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于微地震监测数据的页岩气藏裂缝建模方法 |
CN105629295A (zh) * | 2015-12-29 | 2016-06-01 | 四川圣诺油气工程技术服务有限公司 | 一种页岩气体积压裂微地震监测方法 |
CN105672973A (zh) * | 2016-01-26 | 2016-06-15 | 西南石油大学 | 一种煤层气藏整体体积压裂优化设计方法 |
CN106408208A (zh) * | 2016-10-14 | 2017-02-15 | 东北石油大学 | 体积压裂改造效果评价方法 |
CN108984839A (zh) * | 2018-06-15 | 2018-12-11 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种基于微地震的页岩改造体积快速判断方法 |
Non-Patent Citations (3)
Title |
---|
BINGHOU 等: "Parameter simulation and optimization in channel fracturing", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 * |
侯冰 等: "裂缝性页岩储层多级水力裂缝扩展规律研究", 《岩土工程学报》 * |
刘星 等: "利用微地震事件重构三维缝网", 《石油地球物理勘探》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021004442A1 (zh) * | 2019-07-11 | 2021-01-14 | 中国石油大学(北京) | 一种用于微地震解释压裂裂缝参数结果的校正方法 |
CN110967748A (zh) * | 2019-12-02 | 2020-04-07 | 中国矿业大学(北京) | 微地震成像裂缝提取方法、装置、设备和存储介质 |
CN114718536A (zh) * | 2021-01-05 | 2022-07-08 | 中国石油天然气股份有限公司 | 压裂工艺调整方法 |
CN114718536B (zh) * | 2021-01-05 | 2024-03-26 | 中国石油天然气股份有限公司 | 压裂工艺调整方法 |
CN113283197A (zh) * | 2021-06-10 | 2021-08-20 | 西南石油大学 | 一种基于复杂裂缝缝宽分布的加砂参数设计方法 |
CN113283197B (zh) * | 2021-06-10 | 2022-04-05 | 西南石油大学 | 一种基于复杂裂缝缝宽分布的加砂参数设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110378004B (zh) | 2021-05-07 |
WO2021004442A1 (zh) | 2021-01-14 |
US20220206174A1 (en) | 2022-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110378004A (zh) | 一种用于微地震解释压裂裂缝参数结果的校正方法 | |
Clarkson et al. | Incorporating geomechanical and dynamic hydraulic-fracture-property changes into rate-transient analysis: example from the haynesville shale | |
DE60131181T2 (de) | Untersuchung von mehrschichtigen lagerstätten | |
Khanna et al. | Stimulation of the natural fracture system by graded proppant injection | |
El-Rabba et al. | New perforation pressure-loss correlations for limited-entry fracturing treatments | |
Anderson | Evolution of lumpy glacial landscapes | |
CA3063831A1 (en) | Integrating contextual information into workflow for wellbore operations | |
CN107939365B (zh) | 前置液增能压裂方案生成方法及装置 | |
CN106795748A (zh) | 基于与机械比能的变量相对应的钻探数据来确定完井设计的一个或多个参数 | |
Sato et al. | Post-frac analyses indicating multiple fractures created in a volcanic formation | |
CN110532707A (zh) | 焖井时间确定方法及设备 | |
CN110347931A (zh) | 文章新章节的检测方法及装置 | |
Weiss | Stone tool analysis and context of a new late Middle Paleolithic site in western central Europe–Pouch-Terrassenpfeiler, Ldkr. Anhalt-Bitterfeld, Germany: Eine neue spätmittelpaläolithische Fundstelle im westlichen Mitteleuropa–Pouch-Terrassenpfeiler, Ldkr. Anhalt-Bitterfeld, Deutschland. Steinartefaktanalyse und mitteldeutscher Kontext | |
Konstantinou et al. | Fluid injection through artificially reconstituted bio-cemented sands | |
Kori et al. | Optimum strategies of groundwater pumping regime under scavenger tubewells in lower Indus Basin, Sindh, Pakistan | |
WO2019112469A1 (en) | Method for reservoir stimulation analysis and design based on lagrangian approach | |
Mogensen et al. | Case study-Evaluation of horizontal well multi-stage fracturing in the Viking oil formation | |
Rahim et al. | The effects of mechanical properties and selection of completion interval upon the created and propped fracture dimensions in layered reservoirs | |
Todd et al. | Late Archaic bison hunters in Northern Colorado: 1997-1999 excavations at the Kaplan-Hoover bison bonebed (5LR3953) | |
Mohammed et al. | A production-based model for a fractured well in unconventional reservoirs | |
Hou et al. | Numerical investigation and optimization of multiple fractures in tight gas reservoirs | |
Romberg et al. | Predicting Trouble Stages with Geomechanical Measurements and Machine Learning: A Case Study on Southern Midland Basin Horizontal Completions | |
CN115126459B (zh) | 水力裂缝高度的处理方法和装置 | |
Smallwood | Context and spatial organization of the Clovis assemblage from the Topper site, South Carolina | |
Abrahamson | Is VS30 an effective parameter for site characterization? |
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 |