CN112001134A - 一种非常规气藏流场结构识别方法及装置 - Google Patents

一种非常规气藏流场结构识别方法及装置 Download PDF

Info

Publication number
CN112001134A
CN112001134A CN202010920207.6A CN202010920207A CN112001134A CN 112001134 A CN112001134 A CN 112001134A CN 202010920207 A CN202010920207 A CN 202010920207A CN 112001134 A CN112001134 A CN 112001134A
Authority
CN
China
Prior art keywords
flow field
field structure
stage
reservoir
production rate
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
Application number
CN202010920207.6A
Other languages
English (en)
Other versions
CN112001134B (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202010920207.6A priority Critical patent/CN112001134B/zh
Priority to US17/027,951 priority patent/US20210003004A1/en
Priority to AU2020102430A priority patent/AU2020102430A4/en
Publication of CN112001134A publication Critical patent/CN112001134A/zh
Application granted granted Critical
Publication of CN112001134B publication Critical patent/CN112001134B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • E21B49/008—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 by injection test; by analysing pressure variations in an injection or production test, e.g. for estimating the skin factor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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/30—Specific pattern of wells, e.g. optimising the spacing of wells
    • E21B43/305—Specific pattern of wells, e.g. optimising the spacing of wells comprising at least one inclined or horizontal well
    • 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
    • E21B47/00—Survey of boreholes or wells
    • E21B47/10—Locating fluid leaks, intrusions or movements
    • 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
    • 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/2605—Methods for stimulating production by forming crevices or fractures using gas or liquefied gas
    • 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
    • E21B47/00—Survey of boreholes or wells
    • E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • E21B47/138—Devices entrained in the flow of well-bore fluid for transmitting data, control or actuation signals
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本公开提供一种非常规气藏流场结构识别方法和装置。非常规气藏包括至少一口水平井，在忽略井间干扰的情况下，对于所述至少一口水平井中的任一口水平井，非常规气藏流场结构识别方法包括：在压裂液返排阶段，获取第一生产数据，并根据第一生产数据绘制压裂液返排阶段的实际生产曲线图；获取第二生产数据和不同缝网形态对应的储层物性参数，并对压裂液返排阶段的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制压裂液返排阶段不同缝网形态对应的流场结构版图；将实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构。本公开能够实现对压裂后的储层流场结构的识别。

Description

一种非常规气藏流场结构识别方法及装置
技术领域
本公开涉及非常规天然气开发技术领域，尤其涉及一种非常规气藏流场结构识别方法及装置。
背景技术
非常规油气藏是指油气藏特征、成藏机理及开采技术有别于常规油气藏的石油天然气矿藏。其中，非常规天然气主要包括致密气、页岩气、煤层气等。非常规天然气开发过程中，通常采用对储层进行大规模体积压裂的方式形成裂缝网络，以获得更好的油气流动通道，从而得到更大的油气产量。压裂施工后，每口水平井周围储层所形成的裂缝网络的几何形态及对应的流场结构关系到后期生产能力的评价以及开发方案的设计与调整，因此有必要针对生产井对应的储层进行清晰的流场结构形态描述。
发明内容
本公开的一些实施例提供一种非常规气藏流场结构识别方法及装置，以实现对压裂后的储层流场结构的识别。
一方面，提供一种非常规气藏流场结构识别方法，非常规气藏包括至少一口水平井，在忽略井间干扰的情况下，对于所述至少一口水平井中的任一口水平井，非常规气藏流场结构识别方法包括：在压裂液返排阶段，获取第一生产数据，并根据第一生产数据绘制压裂液返排阶段的实际生产曲线图；获取第二生产数据和不同缝网形态对应的储层物性参数，并对压裂液返排阶段的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制压裂液返排阶段不同缝网形态对应的流场结构版图；将实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构。
在本公开的至少一个实施例中，获取第二生产数据和不同缝网形态对应的储层物性参数，并对水平井的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制不同缝网形态对应的压裂液返排阶段的流场结构版图，包括：预设压裂后的储层缝网形态为羽状缝、网状缝、簇状缝或树状缝，获取第二生产数据，并分别获取羽状缝、网状缝、簇状缝、树状缝四种缝网形态对应的储层物性参数；根据第二生产数据和羽状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的羽状缝流场结构版图；根据第二生产数据和网状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的网状缝流场结构版图；根据第二生产数据和簇状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的簇状缝流场结构版图；根据第二生产数据和树状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的树状缝流场结构版图。
在本公开的至少一个实施例中，网状缝对应的物性参数包括网状缝对应的等效渗透率，其中，网状缝对应的等效渗透率公式为：
Figure BDA0002666470090000021
式中，Ke1为网状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数。
在本公开的至少一个实施例中，簇状缝对应的物性参数包括簇状缝对应的等效渗透率，其中，簇状缝对应的等效渗透率公式为：
Figure BDA0002666470090000022
式中，Ke2为簇状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角。
在本公开的至少一个实施例中，羽状缝对应的物性参数包括羽状缝对应的等效渗透率，其中，羽状缝对应的等效渗透率公式为：
Figure BDA0002666470090000031
式中，Ke3为羽状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角。
在本公开的至少一个实施例中，树状缝对应的物性参数包括树状缝对应的等效渗透率，其中，树状缝对应的等效渗透率公式为：
Figure BDA0002666470090000032
式中，Ke4为树状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角；
l0为0级分叉长度；
D为裂缝分形维数；
Dτ为迂曲度维数；
dmax为初级裂缝最大开度；
m为裂缝最大分叉级数；
rc为裂缝最大延伸长度；
rw为井筒半径。
在本公开的至少一个实施例中，将实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构，包括：将实际生产曲线图分别与羽状缝流场结构版图、网状缝流场结构版图、簇状缝流场结构版图以及树状缝流场结构版图进行对比，根据四种流场结构版图呈现的生产曲线的变化趋势，识别压裂后的储层流场结构。
在本公开的至少一个实施例中，识别压裂后的储层流场结构的分类依据包括：当实际生产曲线图不包括纯产水阶段，无产水速率上升阶段，且套管压力一直下降，则压裂后的储层流场结构为网状缝所形成的流场结构；当实际生产曲线图包括纯产水阶段，且具有产水速率上升阶段与套管压力上升阶段，则压裂后的储层流场结构为簇状缝所形成流场结构；当实际生产曲线图包括纯产水阶段，无产水速率上升阶段，且产气速率维持至少240小时不下降，则压裂后的储层流场结构为树状缝所形成的流场结构；当实际生产曲线图包括两个生产阶段，具有纯产水阶段，无产水速率上升阶段，且产气速率上升后快速下降，则压裂后的储层流场结构为羽状缝所形成的流场结构。
在本公开的至少一个实施例中，第一生产数据包括产水速率、产水量、产气速率、产气量以及套管压力。
另一方面，还提供一种非常规气藏流场结构识别装置，装置包括处理器和存储器，存储器中存储有适于处理器执行的计算机程序指令，计算机程序指令被处理器运行时执行如上述任一实施例所述的非常规气藏流场结构识别方法中的一个或多个步骤。
附图说明
附图示出了本公开的示例性实施方式，并与其说明一起用于解释本公开的原理，其中包括了这些附图以提供对本公开的进一步理解，并且附图包括在本说明书中并构成本说明书的一部分。
图1为根据一些实施例的一种网状缝缝网形态示意图；
图2为根据一些实施例的一种簇状缝缝网形态示意图；
图3为根据一些实施例的一种羽状缝缝网形态示意图；
图4为根据一些实施例的一种树状缝缝网形态示意图；
图5为根据一些实施例的一种非常规气藏流场结构识别方法的实际生产曲线图；
图6为根据一些实施例的一种非常规气藏流场结构识别方法的网状缝流场结构版图；
图7为根据一些实施例的一种非常规气藏流场结构识别方法的簇状缝流场结构版图；
图8为根据一些实施例的一种非常规气藏流场结构识别方法的树状缝流场结构版图；
图9为根据一些实施例的一种非常规气藏流场结构识别方法的羽状缝流场结构版图；
图10为根据一些实施例的一种簇状缝流场结构模拟产气量曲线与实际产气量曲线对比图。
具体实施方式
下面结合附图和实施方式对本公开作进一步的详细说明。可以理解的是，此处所描述的具体实施方式仅用于解释相关内容，而非对本公开的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本公开相关的部分。
需要说明的是，在不冲突的情况下，本公开中的实施方式及实施方式中的特征可以相互组合。下面将参考附图并结合实施方式来详细说明本公开。
需要说明的是，文中的步骤编号，仅为了方便具体实施例的解释，不作为限定步骤执行先后顺序的作用。
本公开一些实施例提供的方法可以由相关的处理器执行，且下文均以处理器作为执行主体为例进行说明。其中，执行主体可以根据具体案例进行调整，如服务器、电子设备、计算机等。
在相关技术中，大多采用微地震监测技术对压裂裂缝的网络进行描述。微地震监测技术通常利用压裂过程中地层压力升高而诱发的微地震，采用地面监测***收集地震波数据，并进行数据的信号识别与处理，从而记录、定位各个微地震源，各个微地震源的分布可以反映地层中的裂缝轮廓。
通过微地震监测技术得到的裂缝形态数据存在一些不足之处。首先，由于目前监测与数据解释技术的发展水平的限制，微地震监测技术的精度有限，往往难以反映地层中真实的裂缝条件。其次，施行微地震等工程需要较高的经济成本。最后，通过微地震监测技术得到的裂缝数据难以应用到油气藏的数值模拟计算过程中，从而无法利用微地震监测技术得到的裂缝数据进行清晰地流场结构识别。
针对目前利用微地震监测技术描述裂缝网络的方式存在的不足，本公开一些实施例提供一种非常规气藏流场结构识别方法，以期更加精确、更加合理的进行非常规气藏流场结构识别。非常规气藏包括至少一口水平井，在忽略井间干扰的情况下，对于所述至少一口水平井中的任一口水平井，非常规气藏流场结构识别方法包括S1～S3。
S1，在压裂液返排阶段，获取第一生产数据，并根据第一生产数据绘制压裂液返排阶段的实际生产曲线图。
在非常规气藏的开采过程中，对储层进行压裂后，即可对储层中的压裂液进行返排，以防止注入的压裂液堵塞储层，并为油气提供流出通道。压裂液返排阶段持续的时间较短，通过对该阶段生产参数的分析可以对储层的压裂效果进行早期判断。压裂液返排后，即可按需进行后期生产。
在一些实施例中，第一生产数据包括产水速率、产水量、产气速率、产气量以及套管压力。
可选的，所述水平井的井口设置有与处理器连接的至少一个压力计，通过上述至少一个压力计可以测量套管内的流体压力(即套管压力)。处理器还与多相分离装置连接。流体流出井口后，通过管道进入多相分离装置。通过对时间与流体体积的计量，可以得到不同流体相(例如气体和液体)的累计产量和速率。处理器采集来自压力计和多相分离装置的数据信号后，对包括产水速率、产水量、产气速率、产气量以及套管压力的数据信号进行处理，从而形成压裂液返排阶段的实际生产曲线图。处理器还与显示器连接，显示器被配置为显示处理器产生的实际生产曲线图。
S2，获取第二生产数据和不同缝网形态对应的储层物性参数，并对压裂液返排阶段的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制压裂液返排阶段不同缝网形态对应的流场结构版图。
可选的，第二生产数据包括焖井时间、喷嘴直径、返排速度、返排时间等参数。储层物性参数包括渗透率、孔隙度、储层厚度、储层杨氏模量、泊松比、储层密度等储层参数，以及裂缝开度、压裂段数、裂缝间距等压裂相关参数，其中，压裂相关参数可以根据实际压裂施工的参数获得。
可选的，采用的数值模拟软件例如为Eclipse、CMG等油气藏数值模拟软件。
可以理解的是，进行数值模拟时选择的生产模式与压裂液返排阶段的实际生产模式相对应。例如，实际压裂液返排阶段没有进行焖井而直接进行返排，则数值模拟过程中即对应没有进行焖井而直接进行返排的情况。又例如，实际压裂液返排阶段对返排速度进行了控制，则数值模拟过程中即对应控制返排速度的情况。是否进行焖井、焖井时间以及控制不同的返排速度等不同的生产模式，均会对形成的流场结构版图中的生产曲线的变化趋势产生影响。
S3，将实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构。
本公开一些实施例提供的非常规气藏流场结构识别方法，通过绘制压裂液返排阶段不同缝网形态对应的流场结构版图，并将压裂液返排阶段的实际生产曲线图与不同缝网形态对应的流场结构版图进行综合对比，能够更加准确地识别压裂后的储层流场结构，从而更加准确地对后期产量进行预测和计算，同时为后期的生产方案的设计与调整提供依据，进而达到提高采收率的目的。
在一些实施例中，S2包括S21～S25。
S21，预设压裂后的储层缝网形态为羽状缝、网状缝、簇状缝或树状缝，获取第二生产数据，并分别获取羽状缝、网状缝、簇状缝、树状缝四种缝网形态对应的储层物性参数。
发明人在研究中发现，非常规气藏压裂后，地层中的缝网形态主要包括网状缝(如图1所示)、簇状缝(如图2所示)、羽状缝(如图3所示)以及树状缝(如图4所示)四种形态。通过将压裂后的储层缝网形态预设为以上四种，能够在保证数值模拟的精确度的同时，简化模拟过程、提高模拟效率。
S22，根据第二生产数据和网状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的网状缝流场结构版图。
在一些实施例中，网状缝对应的物性参数包括网状缝对应的等效渗透率，其中，网状缝对应的等效渗透率公式为：
Figure BDA0002666470090000081
式中，Ke1为网状缝裂缝-基质***整体渗透率，m2；
Km为基质***渗透率，m2；
W为裂缝开度，m；
X为裂缝间距，m；
n为裂缝段数。
S23，根据第二生产数据和簇状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的簇状缝流场结构版图。
在一些实施例中，簇状缝对应的物性参数包括簇状缝对应的等效渗透率，其中，簇状缝对应的等效渗透率公式为：
Figure BDA0002666470090000082
式中，Ke2为簇状缝裂缝-基质***整体渗透率，m2；
Km为基质***渗透率，m2；
W为裂缝开度，m；
X为裂缝间距，m；
n为裂缝段数；
γ为裂缝之间的夹角，rad。
S24，根据第二生产数据和树状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的树状缝流场结构版图。
在一些实施例中，树状缝对应的物性参数包括树状缝对应的等效渗透率，其中，树状缝对应的等效渗透率公式为：
Figure BDA0002666470090000091
式中，Ke4为树状缝裂缝-基质***整体渗透率，m2；
Km为基质***渗透率，m2；
W为裂缝开度，m；
X为裂缝间距，m；
n为裂缝段数；
γ为裂缝之间的夹角，rad；
l0为0级分叉长度，m；
D为裂缝分形维数；
Dτ为迂曲度维数；
dmax为初级裂缝最大开度，m；
m为裂缝最大分叉级数；
rc为裂缝最大延伸长度，m；
rw为井筒半径，m。
S25，根据第二生产数据和羽状缝对应的储层物性参数，对水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的羽状缝流场结构版图。
在一些实施例中，羽状缝对应的物性参数包括羽状缝对应的等效渗透率，其中，羽状缝对应的等效渗透率公式为：
Figure BDA0002666470090000101
式中，Ke3为羽状缝裂缝-基质***整体渗透率，m2；
Km为基质***渗透率，m2；
W为裂缝开度，m；
X为裂缝间距，m；
n为裂缝段数；
γ为裂缝之间的夹角，rad。
通过公式(1)～(4)计算得出的四种缝网形态对应的等效渗透率能够更加准确地反映四种不同缝网形态对应的渗透特性。将计算得出的四种缝网形态对应的等效渗透率分别应用于数值模拟过程中，模拟出的生产曲线的变化趋势更加贴合实际生产规律，从而使获得的压裂液返排阶段不同缝网形态对应的流场结构版图更加准确。
在一些实施例中，S3包括：将实际生产曲线图分别与羽状缝流场结构版图、网状缝流场结构版图、簇状缝流场结构版图以及树状缝流场结构版图进行对比，根据四种流场结构版图呈现的生产曲线的变化趋势，识别压裂后的储层流场结构是羽状缝所形成的流场结构、网状缝所形成的流场结构、簇状缝所形成的流场结构，还是树状缝所形成的流场结构。
下面以某气田中的以后压裂后的水平井A为例，详细介绍一下本公开一些实施例中的非常规气藏流场结构识别方法。
对水平井A所在储层进行压裂施工之后，立即对储层中的压裂液进行返排，即直接进入压裂液返排阶段。在压裂液返排阶段，获取水平井A产水速率、产水量、产气速率、产气量以及套管压力数据，并根据这些数据绘制压裂液返排阶段水平井A的实际生产曲线图，如图5所示。
获取第二生产数据和储层物性参数。焖井时间为0(即不进行焖井直接返排)，返排速度由喷嘴直径控制，取喷嘴直径为10mm，返排时间为200h；储层孔隙度为0.2，储层厚度为5m，储层岩石杨氏模量为4530MPa，泊松比为0.23，储层密度为2350kg/m3；压裂范围为100m，压裂段数为17，平均压裂段间距76.47m，裂缝平均开度为5.42mm。网状缝对应的等效渗透率通过公式(1)计算得出，簇状缝对应的等效渗透率通过公式(2)计算得出，树状缝对应的等效渗透率通过公式(3)计算得出，羽状缝对应的等效渗透率通过公式(4)计算得出。
根据第二生产数据和网状缝对应的储层物性参数，对水平井A的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的网状缝流场结构版图，如图6所示。
图6中，根据网状缝流场结构版图呈现的生产曲线的变化趋势可以看出，生产曲线(或返排曲线)包括一个生产阶段，没有纯产水阶段，无产水速率上升阶段，且套管压力一直下降。因此，当实际生产曲线图不包括纯产水阶段，无产水速率上升阶段，且套管压力一直下降，则压裂后的实际储层流场结构与网状缝所形成的流场结构相符。
根据第二生产数据和簇状缝对应的储层物性参数，对水平井A的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的簇状缝流场结构版图，如图7所示。
图7中，根据簇状缝流场结构版图呈现的生产曲线的变化趋势可以看出，生产曲线可以包括三个生产阶段，即纯产水阶段(Ⅰ1阶段)、产水速率上升与套管压力上升阶段(II1阶段)、两相排出阶段(Ⅲ1阶段)。因此，当实际生产曲线图包括纯产水阶段，且具有产水速率上升与套管压力上升阶段，则压裂后的实际储层流场结构与簇状缝所形成的流场结构相符。
根据第二生产数据和树状缝对应的储层物性参数，对水平井A的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的树状缝流场结构版图，如图8所示。
图8中，根据树状缝流场结构版图呈现的生产曲线的变化趋势可以看出，生产曲线可以包括两个生产阶段，即纯产水阶段(Ⅰ2阶段)和产气速率维持阶段(II2阶段)，并且其不包括产水速率上升阶段。因此，当实际生产曲线图包括纯产水阶段，无产水速率上升阶段，且产气速率维持至少240小时不下降，则压裂后的实际储层流场结构与树状缝所形成的流场结构相符。
根据第二生产数据和羽状缝对应的储层物性参数，对水平井A的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的羽状缝流场结构版图，如图9所示。
图9中，根据羽状缝流场结构版图呈现的生产曲线的变化趋势可以看出，生产曲线可以被划分为两个生产阶段，即纯产水阶段(Ⅰ3阶段)和产气速率上升后快速下降的阶段(II3阶段)，无产水速率上升阶段。因此，当实际生产曲线图包括纯产水阶段，无产水速率上升阶段，且产气速率上升后快速下降，则压裂后的实际储层流场结构与羽状缝所形成的流场结构相符。
基于此，由图5中所示的实际生产曲线图可以看出，实际生产曲线图明显具有纯产水阶段(Ⅰ1阶段)，有产水速率上升与套管压力上升阶段(II1阶段)，以及两相排出阶段(Ⅲ1阶段)，符合簇状缝流场结构版图中生产曲线的变化趋势，故水平井A对应的压裂后储层流场结构被识别为簇状缝流场结构。
在确定水平井A对应的压裂后储层流场结构为簇状缝流场结构的基础上，对水平井A进行后期生产过程模拟。选取岩石基质为黑色页岩，选取缝网形态为簇状缝，压裂范围为100m，等效渗透率为0.135μm2，压裂段数为17，平均压裂段间距76.47m，裂缝平均开度5.42mm，模拟生产时间为850天。产气量数值模拟结果曲线与实际产气量曲线对比如图10所示。从图10中可以看出，数值模拟得到的产气量曲线与实际产气量曲线基本保持一致，表明水平井A对应的压裂后储层流场结构为簇状缝流场结构的结果判断准确。在此基础上对后期的生产方案做出的设计与调整能够符合实际开发需求，达到提高采收率的目的。
本公开的一些实施例还提供一种非常规气藏流场结构识别装置，装置包括处理器和存储器。
处理器用于支持非常规气藏流场结构识别装置执行上述任一实施例所述的非常规气藏流场结构识别方法中的一个或多个步骤。处理器可以是中央处理单元(CentralProcessing Unit，简称CPU)，还可以是其他通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。其中，通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器中存储有适于所述处理器执行的计算机程序指令，所述计算机程序指令被所述处理器运行时执行上述任一实施例所述的非常规气藏流场结构识别方法中的一个或多个步骤。
存储器可以是只读存储器(Read-Only Memory，ROM)或可存储静态信息和指令的其他类型的静态存储设备，随机存取存储器(Random Access Memory，RAM)或者可存储信息和指令的其他类型的动态存储设备，也可以是电可擦可编程只读存储器(ElectricallyErasable Programmable Read-Only Memory，EEPROM)、只读光盘(Compact Disc Read-Only Memory，CD-ROM)或其他光盘存储、光碟存储(包括压缩光碟、激光碟、光碟、数字通用光碟、蓝光光碟等)、磁盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质，但不限于此。存储器可以是独立存在，通过通信总线与处理器相连接。存储器也可以和处理器集成在一起。
在本说明书的描述中，参考术语“一个实施例/方式”、“一些实施例/方式”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例/方式或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例/方式或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例/方式或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例/方式或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例/方式或示例以及不同实施例/方式或示例的特征进行结合和组合。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本公开的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。同时，在本公开的描述中，除非另有明确的规定和限定，术语“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电性连接；可以是直接相连，也可以通过中间媒介间接相连。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本公开中的具体含义。
本领域的技术人员应当理解，上述实施方式仅仅是为了清楚地说明本公开，而并非是对本公开的范围进行限定。对于所属领域的技术人员而言，在上述公开的基础上还可以做出其它变化或变型，并且这些变化或变型仍处于本公开的范围内。

Claims (10)

1.一种非常规气藏流场结构识别方法，其特征在于，所述非常规气藏包括至少一口水平井，在忽略井间干扰的情况下，对于所述至少一口水平井中的任一口水平井，所述非常规气藏流场结构识别方法包括：
在压裂液返排阶段，获取第一生产数据，并根据所述第一生产数据绘制压裂液返排阶段的实际生产曲线图；
获取第二生产数据和不同缝网形态对应的储层物性参数，并对压裂液返排阶段的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制压裂液返排阶段不同缝网形态对应的流场结构版图；
将所述实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构。
2.根据权利要求1所述的非常规气藏流场结构识别方法，其特征在于，所述获取第二生产数据和不同缝网形态对应的储层物性参数，并对所述水平井的产水速率、产气速率以及套管压力进行数值模拟，根据数值模拟的结果，绘制压裂液返排阶段不同缝网形态对应的流场结构版图，包括：
预设压裂后的储层缝网形态为羽状缝、网状缝、簇状缝或树状缝，获取第二生产数据，并分别获取羽状缝、网状缝、簇状缝、树状缝四种缝网形态对应的储层物性参数；
根据所述第二生产数据和网状缝对应的储层物性参数，对所述水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的网状缝流场结构版图；
根据所述第二生产数据和簇状缝对应的储层物性参数，对所述水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的簇状缝流场结构版图；
根据所述第二生产数据和树状缝对应的储层物性参数，对所述水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的树状缝流场结构版图；
根据所述第二生产数据和羽状缝对应的储层物性参数，对所述水平井的产水速率、产气速率以及套管压力进行数值模拟，并根据数值模拟结果，绘制压裂液返排阶段的羽状缝流场结构版图。
3.根据权利要求2所述的非常规气藏流场结构识别方法，其特征在于，网状缝对应的物性参数包括网状缝对应的等效渗透率，其中，网状缝对应的等效渗透率公式为：
Figure FDA0002666470080000021
式中，Ke1为网状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数。
4.根据权利要求2所述的非常规气藏流场结构识别方法，其特征在于，簇状缝对应的物性参数包括簇状缝对应的等效渗透率，其中，簇状缝对应的等效渗透率公式为：
Figure FDA0002666470080000022
式中，Ke2为簇状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角。
5.根据权利要求2所述的非常规气藏流场结构识别方法，其特征在于，羽状缝对应的物性参数包括羽状缝对应的等效渗透率，其中，羽状缝对应的等效渗透率公式为：
Figure FDA0002666470080000031
式中，Ke3为羽状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角。
6.根据权利要求2所述的非常规气藏流场结构识别方法，其特征在于，树状缝对应的物性参数包括树状缝对应的等效渗透率，其中，树状缝对应的等效渗透率公式为：
Figure FDA0002666470080000032
式中，Ke4为树状缝裂缝-基质***整体渗透率；
Km为基质***渗透率；
W为裂缝开度；
X为裂缝间距；
n为裂缝段数；
γ为裂缝之间的夹角；
l0为0级分叉长度；
D为裂缝分形维数；
Dτ为迂曲度维数；
dmax为初级裂缝最大开度；
m为裂缝最大分叉级数；
rc为裂缝最大延伸长度；
rw为井筒半径。
7.根据权利要求2所述的非常规气藏流场结构识别方法，其特征在于，所述将所述实际生产曲线图与不同缝网形态对应的流场结构版图进行对比，识别压裂后的储层流场结构，包括：
将所述实际生产曲线图分别与羽状缝流场结构版图、网状缝流场结构版图、簇状缝流场结构版图以及树状缝流场结构版图进行对比，根据四种流场结构版图呈现的生产曲线的变化趋势，识别压裂后的储层流场结构。
8.根据权利要求7所述的非常规气藏流场结构识别方法，其特征在于，识别压裂后的储层流场结构的分类依据包括：
当所述实际生产曲线图不包括纯产水阶段，无产水速率上升阶段，且套管压力一直下降，则压裂后的储层流场结构为网状缝所形成的流场结构；
当所述实际生产曲线图包括纯产水阶段，且具有产水速率上升阶段与套管压力上升阶段，则压裂后的储层流场结构为簇状缝所形成流场结构；
当所述实际生产曲线图包括纯产水阶段，无产水速率上升阶段，且产气速率维持至少240小时不下降，则压裂后的储层流场结构为树状缝所形成的流场结构；
当所述实际生产曲线图包括纯产水阶段，无产水速率上升阶段，且产气速率上升后快速下降，则压裂后的储层流场结构为羽状缝所形成的流场结构。
9.根据权利要求1～8任一项所述的非常规气藏流场结构识别方法，其特征在于，所述第一生产数据包括产水速率、产水量、产气速率、产气量以及套管压力。
10.一种非常规气藏流场结构识别装置，其特征在于，所述装置包括处理器和存储器，所述存储器中存储有适于所述处理器执行的计算机程序指令，所述计算机程序指令被所述处理器运行时执行如权利要求1～9中任一项所述的非常规气藏流场结构识别方法中的一个或多个步骤。
CN202010920207.6A 2020-09-04 2020-09-04 一种非常规气藏流场结构识别方法及装置 Active CN112001134B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010920207.6A CN112001134B (zh) 2020-09-04 2020-09-04 一种非常规气藏流场结构识别方法及装置
US17/027,951 US20210003004A1 (en) 2020-09-04 2020-09-22 Method and device for identifying structure of flow field in unconventional gas reservoir
AU2020102430A AU2020102430A4 (en) 2020-09-04 2020-09-25 Method and device for identifying structure of flow field in unconventional gas reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010920207.6A CN112001134B (zh) 2020-09-04 2020-09-04 一种非常规气藏流场结构识别方法及装置

Publications (2)

Publication Number Publication Date
CN112001134A true CN112001134A (zh) 2020-11-27
CN112001134B CN112001134B (zh) 2021-09-07

Family

ID=73016676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010920207.6A Active CN112001134B (zh) 2020-09-04 2020-09-04 一种非常规气藏流场结构识别方法及装置

Country Status (3)

Country Link
US (1) US20210003004A1 (zh)
CN (1) CN112001134B (zh)
AU (1) AU2020102430A4 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113779764A (zh) * 2021-08-12 2021-12-10 中国石油大学(北京) 页岩气藏压裂裂缝网络参数的定量评价方法和装置

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113111607B (zh) * 2021-04-15 2022-04-15 西南石油大学 一种油藏流动全耦合的压采一体化数值模拟方法
CN113586969B (zh) * 2021-07-22 2022-11-25 杭州电子科技大学 一种基于准瞬态压力信号的爆管侦测方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105484717A (zh) * 2015-11-24 2016-04-13 西南石油大学 一种提高富有机质页岩基块渗透率的方法
US20180010429A1 (en) * 2015-01-23 2018-01-11 Schlumberger Technology Corporation Control system and method of flowback operations for shale reservoirs
CN108518212A (zh) * 2018-04-09 2018-09-11 西南石油大学 一种计算页岩气藏复杂裂缝网络非稳态产量的方法
CN109101773A (zh) * 2018-09-27 2018-12-28 北京科技大学 致密油水平井立体缝网簇网压裂优化方法
CN109268004A (zh) * 2018-09-21 2019-01-25 北京科技大学 一种页岩气藏介质结构耦合及缝网形态判识方法
CN110210144A (zh) * 2019-06-05 2019-09-06 西南石油大学 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法
CN110219630A (zh) * 2019-06-04 2019-09-10 西南石油大学 一种裂缝性砂岩气藏压裂水平井的压裂液返排量计算方法
US20200032623A1 (en) * 2018-03-21 2020-01-30 ResFrac Corporation Systems and methods for hydraulic fracture treatment and earth engineering for production
CN110924933A (zh) * 2019-11-18 2020-03-27 中国石油集团川庆钻探工程有限公司 一种动态模拟页岩压裂缝网可视化实验方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180010429A1 (en) * 2015-01-23 2018-01-11 Schlumberger Technology Corporation Control system and method of flowback operations for shale reservoirs
CN105484717A (zh) * 2015-11-24 2016-04-13 西南石油大学 一种提高富有机质页岩基块渗透率的方法
US20200032623A1 (en) * 2018-03-21 2020-01-30 ResFrac Corporation Systems and methods for hydraulic fracture treatment and earth engineering for production
CN108518212A (zh) * 2018-04-09 2018-09-11 西南石油大学 一种计算页岩气藏复杂裂缝网络非稳态产量的方法
CN109268004A (zh) * 2018-09-21 2019-01-25 北京科技大学 一种页岩气藏介质结构耦合及缝网形态判识方法
CN109101773A (zh) * 2018-09-27 2018-12-28 北京科技大学 致密油水平井立体缝网簇网压裂优化方法
CN110219630A (zh) * 2019-06-04 2019-09-10 西南石油大学 一种裂缝性砂岩气藏压裂水平井的压裂液返排量计算方法
CN110210144A (zh) * 2019-06-05 2019-09-06 西南石油大学 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法
CN110924933A (zh) * 2019-11-18 2020-03-27 中国石油集团川庆钻探工程有限公司 一种动态模拟页岩压裂缝网可视化实验方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHU WEIYAO，QI QIAN，MA QIAN，DENG JIA，YUE MING， LIU YUZHANG: "Unstable seepage modeling and pressure propagation of shale gas reservoirs", 《PETROLEUM EXPLORATION AND DEVELOPMENT》 *
朱维耀，马东旭，亓倩，宋智勇，岳明: "复杂缝网页岩压裂水平井多区耦合产能分析", 《天然气工业》 *
杨柳，石富坤，赵逸清: "复杂结构井在页岩气开发中的应用进展", 《科学技术与工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113779764A (zh) * 2021-08-12 2021-12-10 中国石油大学(北京) 页岩气藏压裂裂缝网络参数的定量评价方法和装置
CN113779764B (zh) * 2021-08-12 2022-07-29 中国石油大学(北京) 页岩气藏压裂裂缝网络参数的定量评价方法和装置

Also Published As

Publication number Publication date
AU2020102430A4 (en) 2020-11-05
US20210003004A1 (en) 2021-01-07
CN112001134B (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
CN112001134B (zh) 一种非常规气藏流场结构识别方法及装置
EP1774136B1 (en) Well characterisation method
US11236596B2 (en) Real-time diversion control for stimulation treatments using fiber optics with fully-coupled diversion models
CN103413030B (zh) 一种缝洞型碳酸盐岩气藏动态分析方法及***
EP3152391B1 (en) Method of estimating well productivity along a section of a wellbore
US20020177955A1 (en) Completions architecture
US11499392B2 (en) Designing a wellbore cement sheath in compacting or subsiding formations
AU2020481721A1 (en) Real-time parameter adjustment in wellbore drilling operations
CN110439519A (zh) 一种基于极限限流设计的压裂方法及***
US11867034B2 (en) Systems and methods for automated gas lift monitoring
Rodriguez Inferences of two dynamic processes on recovery factor and well spacing for a shale oil reservoir
WO2021126963A1 (en) Method of predicting and preventing an event of fracture hit
US20190186237A1 (en) Method and control system for allocating steam to multiple wells in steam assisted gravity drainage (sagd) resource production
CN109726450B (zh) 确定页岩气藏水平井重复压裂时机的方法及设备
CN112345429B (zh) 非常规油气藏非均匀地质参数的获取方法及装置
KR101920689B1 (ko) 수평정의 수직 위치 결정방법
CN113107474A (zh) 利用短期关井数据预测低渗气藏的储层参数的方法及装置
Celier et al. Zuidwal: a gas field development with horizontal wells
Yeleussinov et al. Evaluation of vertical fracture propagation using combined spectral noise logging and high precision temperature logging data in a low permeability sandstone reservoir, Western Kaz
Kremenetsky et al. New possibilities of well testing and production logging in horizontal wells with non-uniform inflow profile
Wilson Wireless Downhole Gauges Help Maximize Value of Appraisal Test in Abandoned Well
CN116556934A (zh) 一种含水致密气藏优质储量区域的开发方法、装置和设备
CN114526042A (zh) 一种针对长井段裸眼井的分段设计方法及***
WO2021230959A1 (en) Methods and systems for selecting inflow control device design simulations based on case selection factor determinations
CN115704307A (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