CN110378004B - 一种用于微地震解释压裂裂缝参数结果的校正方法 - Google Patents

一种用于微地震解释压裂裂缝参数结果的校正方法 Download PDF

Info

Publication number
CN110378004B
CN110378004B CN201910626869.XA CN201910626869A CN110378004B CN 110378004 B CN110378004 B CN 110378004B CN 201910626869 A CN201910626869 A CN 201910626869A CN 110378004 B CN110378004 B CN 110378004B
Authority
CN
China
Prior art keywords
fracturing
fracture
well
adding amount
microseism
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
Application number
CN201910626869.XA
Other languages
English (en)
Other versions
CN110378004A (zh
Inventor
侯冰
肖月
金衍
曹博文
陈勉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201910626869.XA priority Critical patent/CN110378004B/zh
Publication of CN110378004A publication Critical patent/CN110378004A/zh
Priority to US17/593,927 priority patent/US20220206174A1/en
Priority to PCT/CN2020/100548 priority patent/WO2021004442A1/zh
Application granted granted Critical
Publication of CN110378004B publication Critical patent/CN110378004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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
    • 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 DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP 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 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
    • 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. analysis, for interpretation, for correction
    • 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 DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP 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为：
Figure DDA0002126341410000011
其中，[L5/4H]W为利用微地震监测解释结果的计算值；[L5/4H]T为利用现场施工参数得到的理论计算值；C、根据计算得到的各所述压裂井的单井体积校正系数，计算各类压裂井的类体积校正系数。本发明结合现场的施工参数能够对微地震解释结果进行校正，从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态，进而可以辅助微地震位置精确反演，有效指导水力压裂施工作业。

Description

一种用于微地震解释压裂裂缝参数结果的校正方法
技术领域
本发明涉及微地震监测解释技术领域，具体涉及一种用于微地震解释压裂裂缝参数结果的校正方法。
背景技术
随着矿藏开采技术的不断开发，水力裂缝遇到天然裂缝的行为受到越来越多人的关注。水力压裂是公认的在致密裂隙储层中提高采收率的主要增产技术。水力压裂产生的水力裂缝可激活天然裂缝，形成复杂的弥散性裂缝缝网，起到增渗、增产的效果，因此对水力裂缝与天然裂缝相交模式及作用效果的研究具有极高的应用价值。
微地震监测是解释压裂时人工裂缝(即水力裂缝)缝长、缝高及方位的有效方法之一。水力压裂微地震裂缝监测技术可以用于指示裂缝位置、分析裂缝发育情况，并辅助微地震位置精确反演，指导水力压裂施工作业。水力压裂微地震裂缝监测技术是通过监测水力压裂过程中产生的微地震事件，来解释压后人工裂缝形态的参数。
目前研究重点大多是在改进微地震重构理论算法上，微地震监测过程中由于噪点影响导致微地震解释裂缝参数失真，微地震解释裂缝形态与现场压裂过程中的加砂量和液量时常不具备对应关系，导致微地震反演出的裂缝参数不具备参考价值，严重影响了对压后裂缝形态以及压裂效果的评价。
因此有必要发明一种可以用于微地震解释压裂裂缝参数结果的校正方法，以能够更准确真实地推导出裂缝参数，指示裂缝位置、分析裂缝发育情况，辅助微地震位置精确反演，指导水力压裂施工作业。
发明内容
高排量压裂施工时会产生大量微地震信号事件点，利用微地震监测数据难以真实解释出压裂后的有效裂缝形态，因此有必要结合现场的施工参数对对微地震解释结果进行校正。本发明旨在提供一种用于微地震解释压裂裂缝参数结果的校正方法。
本发明提供一种用于微地震解释压裂裂缝参数结果的校正方法，所述方法包括：
A、根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井(部分取样，作为学习样本)的沟通模式进行分类；
B、对分类后的各所述压裂井单井进行体积系数校正，各所述压裂井的体积校正系数BV为：
Figure BDA0002126341390000021
其中，[L5/4H]W为利用微地震监测解释结果的计算值；[L5/4H]T为利用现场施工参数得到的理论计算值，
Figure BDA0002126341390000022
而Vi＝Vfp+Vfs，
Figure BDA0002126341390000023
上式中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
其中，
Figure BDA0002126341390000031
Figure BDA0002126341390000032
为所述压裂井的类体积校正系数。
上述的用于微地震解释压裂裂缝参数结果的校正方法，根据所述体积校正系数得到缝高校正系数BH：
BH＝a
其中，
Figure BDA0002126341390000033
Figure BDA0002126341390000034
为所述压裂井的类体积校正系数。
上述的用于微地震解释压裂裂缝参数结果的校正方法，所述压裂井的裂缝为对称的双翼垂直裂缝，裂缝井的水平截面形状为椭圆形。
本发明结合现场的施工参数能够对微地震解释结果进行校正，从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态，进而可以辅助微地震位置精确反演，有效指导水力压裂施工作业。
附图说明
图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。
根据公式，计算得到第一主成分贡献率，看是否可以满足要求。
如果可以满足要求，即可得到第一主成分对应的特征向量为
Figure BDA0002126341390000051
Figure BDA0002126341390000052
故第一主成分可以表示为
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的概念，其表示前置液与携砂液的总的造缝效率。
Figure BDA0002126341390000071
Ef根据上式中第二个等式以及现场压裂施工参数可以直接计算得到。
式中：
Vfrac-单侧裂缝体积，m3；
根据上述PKN模型，最大裂缝宽度Wmax可由下式计算得到。
Figure BDA0002126341390000072
式中：
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)联立，得到
Figure BDA0002126341390000073
联立公式(6)和公式(4)得到
Figure BDA0002126341390000081
将L5/4H与微地震监测结果做比值，即可得到体积校正系数Bv，即：
Figure BDA0002126341390000082
式中：
[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。
上述推导得到的校正系数为单井的校正系数，而为了实际描述各类压裂井的裂缝形态，需要得到各类的校正系数，得到各类的校正系数，才能够对裂缝形态进行校正。
计算各类压裂井的类体积校正系数，各类的类体积校正系数取该类压裂井中单井体积校正系数计算结果的平均值。则计算各类的缝宽与缝高校正系数时，
Figure BDA0002126341390000083
Figure BDA0002126341390000084
为所述压裂井的类体积校正系数。
本发明结合现场的施工参数能够对微地震解释结果进行校正，从而能够利用微地震监测数据较真实地解释出压裂后的有效裂缝形态，进而可以辅助微地震位置精确反演，有效指导水力压裂施工作业。
实施例
对切六区块36口压裂井的加砂量以及微地震解释裂缝长度进行统计，得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵，随后对观测矩阵进行计算得到均值向量为：
Figure BDA0002126341390000091
根据公式计算得到切六区块二维随机变量样本的A矩阵为：
Figure BDA0002126341390000092
根据数理统计相关知识，用样本的无偏估计来代替总体协方差。
Figure BDA0002126341390000093
由样本的协方差矩阵可以发现，其数值相差较大，为防止大数吃小数的现象，在求解切六区块的主成分的时候可以基于相关系数矩阵求解，根据上述的样本协方差矩阵，求解得到切六区块分类过程中二维随机变量的相关系数矩阵。
Figure BDA0002126341390000094
根据相关系数矩阵，计算得到了切六区块分类相关系数矩阵的特征向量。
λ1＝1.58
λ2＝0.42
根据公式，计算得到第一主成分贡献率为0.79，可以满足要求。
随后得到第一主成分对应的特征向量为
Figure BDA0002126341390000095
故第一主成分可以表示为
F1＝0.5x1+0.87x2
x1为加砂量，x2为微地震解释裂缝长度(双翼缝长之和)。
由于切六区块压裂井较多，可按照水力裂缝沟通天然裂缝的复杂程度将二类及三类分别再细分成两个亚类，随后确定各类的临界值。
(1)基于主成分分析的分类标准
显然，主成分得分越高表示施工加砂量以及裂缝长度越大，则形成天然裂缝水力裂缝沟通模式也更为复杂，结合现场分析需要给出基于主成分得分对沟通模式进行分类的分类标准，分类标准如表1.1所示。
表1.1基于主成分得分的分类标准表
Figure BDA0002126341390000101
(2)基于加砂量的分类标准
根据现场施工需要，认为当裂缝长度与加砂量比值小于10时，认为现场实际的施工加砂量与微地震解释缝长不能很好对应，微地震解释结果不具有参考价值，微地震解释裂缝长度不能作为分类标准，而以加砂量作为沟通模式分类的标准，其标准如表1.2所示。
表1.2切六区块基于加砂量的分类标准表
Figure BDA0002126341390000102
将待判的36口井参数逐一输入所编制的软件，可得到切六区块待判的36口井的归类结果，如表1.3所示。
表1.3切六区块压裂井分类结果表
Figure BDA0002126341390000111
根据切六区块压裂井的分类结果，分别根据本发明的推导思路得到一类井、二类井和三类井及其亚类的体积校正数据BV，类的体积校正系数取该类井中单井体积校正系数计算结果的平均值。在这里，取切六区块岩石的弹性模量为336GPa，泊松比为0.2，裂缝半缝长取裂缝微地震监测结果两个半翼缝长的平均值。切六区块各类体积校正系数计算结果如表2.1所示。
表2.1切六区块体积校正系数计算结果表
Figure BDA0002126341390000112
根据切六区块的体积校正系数计算结果，依据公式(9)、(10)、(11)得到了各类井及其亚类的裂缝半缝长校正系数BL和缝高校正系数BH，其统计结果见表2.2。
表2.2切六区块压裂井各校正系数统计表
Figure BDA0002126341390000121
分别以切六区块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 (8)

1.一种用于微地震解释压裂裂缝参数结果的校正方法，所述方法包括：
A、根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类；根据微地震监测结果与加砂量的吻合程度对待校正区块取样压裂井的沟通模式进行分类的标准为：当微地震解释的水力裂缝总的长度与加砂量的比值<10时，将微地震解释的水力裂缝与对应的加砂量作为分类的依据；当微地震解释的水力裂缝长度与加砂量的比值>＝10时，将加砂量作为分类的依据：
当将微地震解释的水力裂缝与对应的加砂量作为分类的依据时，利用数理分析计算微地震解释的水力裂缝长度与加砂量的主成分，将该主成分值作为分类的依据，微地震解释裂缝长度与加砂量的主成分得分越高表示施工加砂量及所述裂缝长度越大，反之亦然；
B、对分类后的各所述压裂井单井进行体积系数校正，各所述压裂井的体积校正系数BV为：
Figure FDA0002826313100000011
其中，[L5/4H]W为利用微地震监测解释结果的计算值；[L5/4H]T为利用现场施工参数得到的理论计算值，
Figure FDA0002826313100000012
而Vi＝Vfp+Vfs，
Figure FDA0002826313100000013
上式中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所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，当所述主成分值<130时，将所述压裂井分类为第一类；当所述主成分值为130-160时，将所述压裂井分类为第二类第I型；当所述主成分值为160-200，将所述压裂井分类为第二类第II型；当所述主成分值为200-250时，将所述压裂井分类为第三类第I型；当所述主成分值>250时，将所述压裂井分类为第三类第II型。
3.根据权利要求2所述的用于微地震解释压裂裂缝参数结果的校正方法，进行主成分计算时，首先对带校正区块部分样本压裂井的加砂量以及微地震解释裂缝长度进行统计，得到由加砂量和裂缝长度组成的二维随机变量的观测矩阵，随后对观测矩阵进行计算得到均值向量；之后计算得到该区块二维随机变量样本的A矩阵；再用样本的无偏估计来代替总体协方差S*；再计算得到区块分类相关系数矩阵的特征向量；如果随后计算得到第一主成分贡献率不能满足要求，则添加第二主成分；如果随后计算得到第一主成分贡献率满足要求，则以该第一主成分值作为依据进行分类。
4.根据权利要求3所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，如果所述协方差的矩阵数值相差较大，在求解所述区块的取样压裂井的主成分的时候基于相关系数矩阵求解，根据样本协方差矩阵，求解得到所述区块的所述二维随机变量的相关系数矩阵；再根据所述相关系数矩阵，计算得到所述区块的相关系数矩阵的所述特征向量。
5.根据权利要求1所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，当所述加砂量<10m3时，将所述压裂井分类为第一类；当所述加砂量为10-13m3时，将所述压裂井分类为第二类第I型；当所述加砂量为13-16m3时，将所述压裂井分类为第二类第II型；当所述加砂量为16-20m3时，将所述压裂井分类为第三类第I型；当所述加砂量>20m3时，将所述压裂井分类为第三类第II型。
6.根据权利要求1至5任一项所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，根据各类所述压裂井的类体积校正系数得到各类所述压裂井的类裂缝半缝长校正系数BL，BL＝a5/4
其中，
Figure FDA0002826313100000021
Figure FDA0002826313100000023
为所述压裂井的类体积校正系数。
7.根据权利要求1至5任一项所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，根据所述体积校正系数得到缝高校正系数BH，
BH＝a，其中，
Figure FDA0002826313100000022
Figure FDA0002826313100000024
为所述压裂井的类体积校正系数。
8.根据权利要求7所述的用于微地震解释压裂裂缝参数结果的校正方法，其特征在于，所述压裂井的裂缝为对称的双翼垂直裂缝，裂缝井的水平截面形状为椭圆形。
CN201910626869.XA 2019-07-11 2019-07-11 一种用于微地震解释压裂裂缝参数结果的校正方法 Active CN110378004B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201910626869.XA CN110378004B (zh) 2019-07-11 2019-07-11 一种用于微地震解释压裂裂缝参数结果的校正方法
US17/593,927 US20220206174A1 (en) 2019-07-11 2020-07-07 Correction method for micro seismic interpretation fracturing fracture parameter result
PCT/CN2020/100548 WO2021004442A1 (zh) 2019-07-11 2020-07-07 一种用于微地震解释压裂裂缝参数结果的校正方法

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 CN110378004A (zh) 2019-10-25
CN110378004B true 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)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110378004B (zh) * 2019-07-11 2021-05-07 中国石油大学(北京) 一种用于微地震解释压裂裂缝参数结果的校正方法
CN110967748A (zh) * 2019-12-02 2020-04-07 中国矿业大学(北京) 微地震成像裂缝提取方法、装置、设备和存储介质
CN114718536B (zh) * 2021-01-05 2024-03-26 中国石油天然气股份有限公司 压裂工艺调整方法
CN113283197B (zh) * 2021-06-10 2022-04-05 西南石油大学 一种基于复杂裂缝缝宽分布的加砂参数设计方法
CN114542040B (zh) * 2022-02-25 2024-04-26 山西蓝焰煤层气集团有限责任公司 一种连续脉冲水力压裂系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297783A (zh) * 2014-06-30 2015-01-21 中国石油天然气集团公司 一种水力压裂微地震事件的解释方法以及系统
CN105672973A (zh) * 2016-01-26 2016-06-15 西南石油大学 一种煤层气藏整体体积压裂优化设计方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100256964A1 (en) * 2009-04-07 2010-10-07 Schlumberger Technology Corporation System and technique to quantify a fracture system
US9618652B2 (en) * 2011-11-04 2017-04-11 Schlumberger Technology Corporation Method of calibrating fracture geometry to microseismic events
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
CN104459775B (zh) * 2014-11-28 2017-02-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于微地震监测数据的页岩气藏裂缝建模方法
CN104730596B (zh) * 2015-01-25 2017-11-07 中国石油大学(华东) 一种基于多尺度因素约束的离散裂缝建模方法
CN105629295B (zh) * 2015-12-29 2019-11-29 四川圣诺油气工程技术服务有限公司 一种页岩气体积压裂微地震监测方法
CN106408208B (zh) * 2016-10-14 2019-05-28 东北石油大学 体积压裂改造效果评价方法
CN108984839B (zh) * 2018-06-15 2022-02-18 中国石油化工股份有限公司 一种基于微地震的页岩改造体积快速判断方法
CN110378004B (zh) * 2019-07-11 2021-05-07 中国石油大学(北京) 一种用于微地震解释压裂裂缝参数结果的校正方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297783A (zh) * 2014-06-30 2015-01-21 中国石油天然气集团公司 一种水力压裂微地震事件的解释方法以及系统
CN105672973A (zh) * 2016-01-26 2016-06-15 西南石油大学 一种煤层气藏整体体积压裂优化设计方法

Also Published As

Publication number Publication date
CN110378004A (zh) 2019-10-25
US20220206174A1 (en) 2022-06-30
WO2021004442A1 (zh) 2021-01-14

Similar Documents

Publication Publication Date Title
CN110378004B (zh) 一种用于微地震解释压裂裂缝参数结果的校正方法
US10242130B2 (en) System, method and computer program product for wellbore event modeling using rimlier data
CN104126052A (zh) 用于确定裂隙间距的方法和使用所述方法的井压裂
CN112502701B (zh) 一种低渗透储层综合地质-工程的分类评价方法
CN105525909A (zh) 分析油藏储层非均质性质的方法
CN103590824A (zh) 经过多段压裂改造后的致密气藏水平井的产能计算方法
CN109063403B (zh) 一种滑溜水压裂优化设计方法
CN109236273A (zh) 油田开发生产动态数据处理方法
CN106651158B (zh) 一种超低渗致密油藏水平井注水开发见效程度定量评价方法
Bin et al. Intelligent identification and real-time warning method of diverse complex events in horizontal well fracturing
CN106777663B (zh) 一种考虑天然裂缝的压裂液滤失速度计算方法
CN112699554B (zh) 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法
Qu et al. Evaluation of hydraulic fracturing of horizontal wells in tight reservoirs based on the deep neural network with physical constraints
CN111553067B (zh) 页岩油藏压裂水平井不停井试井测试和解释方法、装置
CN109426689B (zh) 水平井压裂裂缝的评价方法及系统
Sesetty et al. Modeling dense-arrays of hydraulic fracture clusters: Fracture complexity, net pressure and model calibration
Jiang et al. Multiple fracturing parameters optimization for horizontal gas well using a novel hybrid method
CN109339760B (zh) 一种水平井一段多簇压裂裂缝条数诊断方法
CN116362121A (zh) 水平井压裂的裂缝参数确定方法及装置
CN111311119A (zh) 一种海相页岩气储层缝网改造能力的地质评价方法
CN114169261B (zh) 一种基于g函数曲线分析的多裂缝参数反演方法和装置
CN107832900A (zh) 砂砾岩油藏注水效果评价方法
CN110067555A (zh) 碳酸盐岩油井的最小动态储量的确定方法和装置
US20240076973A1 (en) System and methods for determining the effect of fracture interference on shale well performance
CN114065529B (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