CN113111607A - 一种油藏流动全耦合的压采一体化数值模拟方法 - Google Patents
一种油藏流动全耦合的压采一体化数值模拟方法 Download PDFInfo
- Publication number
- CN113111607A CN113111607A CN202110422759.9A CN202110422759A CN113111607A CN 113111607 A CN113111607 A CN 113111607A CN 202110422759 A CN202110422759 A CN 202110422759A CN 113111607 A CN113111607 A CN 113111607A
- Authority
- CN
- China
- Prior art keywords
- fracture
- artificial
- matrix
- grid
- reservoir
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 77
- 238000004088 simulation Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000010168 coupling process Methods 0.000 title claims abstract description 9
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims abstract description 102
- 239000007791 liquid phase Substances 0.000 claims abstract description 64
- 239000012071 phase Substances 0.000 claims abstract description 58
- 239000012530 fluid Substances 0.000 claims abstract description 38
- 239000003129 oil well Substances 0.000 claims abstract description 34
- 238000010276 construction Methods 0.000 claims abstract description 33
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000008878 coupling Effects 0.000 claims abstract description 6
- 230000000977 initiatory effect Effects 0.000 claims abstract description 4
- 206010017076 Fracture Diseases 0.000 claims description 159
- 208000010392 Bone Fractures Diseases 0.000 claims description 131
- 208000013201 Stress fracture Diseases 0.000 claims description 22
- 230000035699 permeability Effects 0.000 claims description 13
- 239000007788 liquid Substances 0.000 claims description 12
- 239000011435 rock Substances 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000005465 channeling Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000005086 pumping Methods 0.000 claims description 2
- 239000000758 substrate Substances 0.000 claims description 2
- 230000008595 infiltration Effects 0.000 claims 1
- 238000001764 infiltration Methods 0.000 claims 1
- 238000007670 refining Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract 1
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/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
-
- 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)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种油藏流动全耦合的压采一体化数值模拟方法,依次包括以下步骤:(1)建立结构化储层网格和裂缝初始单元;(2)计算压裂过程中的人工裂缝半长、人工裂缝宽度分布、人工裂缝内流体压力分布、基质油相压力分布和基质液相饱和度分布;(3)计算油井返排‑生产过程中储层压力分布和储层含水饱和度分布;(4)计算油井的产能。本发明采用结构化网格和嵌入式离散裂缝模型,使用结构化网格能够有利于油藏边界问题的处理,嵌入式离散裂缝模型具有网格划分简单、计算快速、易于处理复杂裂缝等优点。本发明用于模拟油藏开发的压裂、返排和生产过程,通过一体化结合得到更贴近实际情况的油井生产动态,从而指导压裂施工和生产制度设计。
Description
技术领域
本发明属于油气田开发领域,具体涉及一种基于嵌入式离散裂缝模型的油藏流动全耦合的压采一体化数值模拟方法。
背景技术
随着油气资源的不断开发,数值模拟技术在油气田开发中得到了极大的应用。在水力压裂、压后返排、油井生产等领域,通过数值模拟方法能够对施工效果和油井产能进行预测分析,从而对实际的施工、生产提供可靠的指导依据。
目前油气藏开发过程的数值模拟中,压裂施工的数值模拟与油井返排生产过程的数值模拟相对独立,很少有油藏压裂-返排-生产相互结合的一体化数值模拟方法。在水力压裂施工过程中,由于压裂液的滤失和裂缝扩展,油藏的压力场、含水饱和度场均会发生变化,裂缝周围一定范围内的储层流体压力和含水饱和度会得到提升。因此压裂施工结束时与原始初始状态下的储层压力场、饱和度场有着较大的差异,若想要获得与实际生产情况贴切的油井返排-生产模拟结果,就应当运用压裂模拟结束时的储层状态参数进行返排-生产模拟。由此可知,油藏压采一体化数值模拟方法的提出有着重要的意义,油藏流动全耦合的压采一体化模拟可以获得更符合油气井实际生产情况的模拟结果,为油藏开发提供可靠依据。
目前,油藏压采一体化数值模拟方面的研究还比较少。现有的研究主要基于三角形非结构化网格提出了一种耦合页岩气藏水力压裂返排生产过程的一体化模拟方法(唐慧莹,梁海鹏,张东旭,张烈辉.一种耦合页岩气藏水力压裂返排生产过程一体化模拟方法[P].中国专利:CN110863810A,2020.03.06)。但该方法运用非结构化网格进行模拟降低了网格利用率,对运算设备也有着较高的要求,同时无法实现复杂裂缝的模拟。该方法虽然给出了压裂-生产过程模拟的思路,但模型中未考虑储层流动的全耦合,没有记录水力裂缝扩展过程中的压力变化历史,未能将压裂过程与返排生产过程较好结合,无法模拟超低含水饱和度油藏中的压裂液水相圈闭损害。
因此,本发明采用了结构化网格和嵌入式离散裂缝模型进行模拟,建立一种油藏流动全耦合的压采一体化数值模拟方法。使用结构化网格能够有利于油藏边界问题的处理,提升计算效率;同时嵌入式离散裂缝模型还具有网格划分简单、计算快速、易于处理复杂裂缝等优点。
发明内容
本发明的目的在于提供一种基于嵌入式离散裂缝模型的油藏流动全耦合压采一体化数值模拟方法,用于模拟油藏开发中的压裂、返排和生产过程,通过一体化结合能够模拟得到更贴近实际情况的油井生产动态,从而指导压裂施工和生产制度设计,弥补了现有技术的不足。
一种油藏流动全耦合的压采一体化数值模拟方法,依次包括以下步骤:
(1)建立结构化储层网格和裂缝初始单元。
首先,在x-y直角坐标系下建立储层的结构化网格,将储层长度Lx划分成ni段,储层宽度Ly划分为nj段,整个储层便被划分成了ni×nj的结构化网格。xi,j和yi,j分别代表着每个网格的长度和宽度,下标i和j代表着每一个网格在储层中的位置。
随后,在储层网格基础上,建立人工裂缝的初始单元。对压裂延伸位置处一定范围内的储层网格进行细化,并在细化后的网格上添加一段有3个网格长度的较小人工裂缝,裂缝延伸方向为y方向,将其作为人工裂缝的初始单元。在压裂阶段裂缝会沿着裂缝初始单元的方向进行延伸。由于本发明选用嵌入式离散裂缝模型来建立人工裂缝,储层网格中人工裂缝会被基质网格切分为具有一定长度的线段。每一个线段代表一个人工裂缝单元,因此人工裂缝的初始单元数为nL0,人工裂缝单元总数表示为nL,任意裂缝单元编号为L,每个人工裂缝单元长度为ξL,随着裂缝的延伸扩展,nL会逐渐增大,直至压裂结束。
(2)计算压裂过程中的人工裂缝半长、人工裂缝宽度分布、人工裂缝内流体压力分布、基质油相压力分布和基质液相饱和度分布。
在已建立的油藏网格的基础上,运用水力压裂模型进行数值计算,获得压裂施工过程中的人工裂缝单元总数、裂缝单元的宽度、裂缝单元内的流体压力、人工裂缝半长、基质网格的油相压力以及基质网格的液相饱和度。
水力压裂模型中人工裂缝内为单相流动。水力压裂模型如下:
a)考虑流体滤失的裂缝宽度模型(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998):
式中,W(y)——人工裂缝任意位置处宽度,m;
H——人工裂缝高度,m;
v——储层岩样的泊松比,无因次;
E——储层岩样的杨氏模量,MPa;
σh——水平最小主应力,MPa;
Pm(y)——人工裂缝任意位置处的基质流体压力,MPa。
μ——压裂液的黏度,mPa·s;
vl——压裂液滤失速度,m/s;
t——时间,s;
kmf——人工裂缝与基质的平均渗透率,mD;
b)压裂过程中的基质渗流模型
Pmc=Pmo-Pml (6)
式中,Amf——裂缝与基质的接触面积,m;
KF——人工裂缝的渗透率,D;
VF——人工裂缝单元的体积,m3;
KFrl——人工裂缝液相相对渗透率,无量纲;
PF——人工裂缝流体压力,MPa;
φF——人工裂缝的孔隙度,无因次;
φm——储层基质孔隙度,无因次;
Km——基质渗透率,mD;
Kmrl——基质中液相的相对渗透率,无因次;
Kmro——基质中油相的相对渗透率,无因次;
Sml——基质中的液相饱和度,无因次;
Vb——基质单元体积,m3;
μl——基质中的液相黏度,mPa·s;
μo——基质中的油相黏度,mPa·s;
Bl——基质中液相的体积系数,无因次;
Bo——基质中油相的体积系数,无因次;
Pml、Pmo——基质中液相、油相的压力,MPa;
Pmc——基质中的毛管压力,MPa;
δm——基质网格是否含有人工裂缝的判断参数,由于有裂缝穿过的储层基质网格会发生基质-裂缝的流体交换,因此当储层基质网格有裂缝穿过时δm=1;当储层基质网格无裂缝穿过时δm=0。
c)油藏渗流初始条件为:
Pmo(i,j,t)|t=0=Pe (7)
式中,Pe——油藏原始地层压力,MPa;
i,j——网格的位置坐标。
d)裂缝延伸边界条件为:
PFtL=1,t=PZ (9)
式中,G——储层岩样的体积模量,MPa;
PFL=1,t——压裂施工中的井底流体压力,即第1段人工裂缝单元内的流体压力,MPa;
PZ——压裂施工中的井底泵注压力,MPa;
nL——压裂施工时间t下的人工裂缝单元的总数;
ξL——人工裂缝单元的长度,m。
e)油藏基质渗流边界条件为:
式中,Lx,Ly——分别表示储层长度和储层宽度,m;
对上述压裂模型进行数值求解可得到任意时间下的人工裂缝任意位置处的宽度W(y)、人工裂缝任意位置处的流体压力PF、基质油相压力Pmo、基质液相饱和度Sml。由此可以得到,压裂施工时间t下的每个裂缝单元的宽度WL,t、每个裂缝单元的流体压力PFL,t、每个基质网格的油相压力Pmoi,j,t和每个基质网格的液相饱和度Smli,j,t,以及压裂施工时间t下的尖端裂缝单元内的流体压力
将计算得到的代入下式,则可以得到压裂施工时间t下的裂缝尖端应力强度因子KIi,t。(Abbas S,Gordeliy E,Peirce A,et al.Limited Height Growth and ReducedOpening of Hydraulic Fractures due to Fracture Offsets:An XFEM Application[C].SPE Hydraulic Fracturing Technology Conference.Society of PetroleumEngineers,2014)
式中,KIi,t——压裂施工时间t下的裂缝尖端应力强度因子,MPa·m1/2。
将计算得到的裂缝尖端应力强度因子KIi,t与储层岩石的断裂韧性KIC(取KIC=2MPa·m1/2)进行比较。当KIi,t≤KIC时,压裂施工时间t下的人工裂缝不发生裂缝扩展,裂缝长度不变,重复步骤(2)进行时间t+1下的计算,时间t+1下的人工裂缝单元总数仍为nL;当KIi,t>KIC时,人工裂缝沿y方向发生扩展,裂缝长度增加,重复步骤(2)进行时间t+1下的计算,时间t+1下的人工裂缝单元总数为nL=nL+1。
计算直到压裂施工完成时t=tend结束,此时可获得时间t=tend时的人工裂缝单元总数nL、每个裂缝单元的宽度每个裂缝单元内的流体压力人工裂缝半长nL×ξL、每个基质网格的油相压力和每个基质网格的液相饱和度
(3)计算油井返排-生产过程中储层压力分布和储层含水饱和度分布。
将已计算获得的人工裂缝半长、基质网格的油相压力分布、基质网格的液相饱和度分布等参数当作油井返排-生产过程中的初始参数来进行模拟。
本发明中的油井的返排-生产模型运用双重介质模型和嵌入式离散裂缝模型进行建立,考虑了人工裂缝、微裂缝、基质、人工裂缝-微裂缝之间和基质-微裂缝之间的流体渗流。模型中储层网格包含了基质和微裂缝,同时人工裂缝中忽略毛管力作用,液相压力与油相压力相同,且均等于人工裂缝单元内的流体压力PF。油井返排-生产模型如下(Lee S H,Lough M.F,Jensen C.L.Hierarchical modeling of flow in naturally fracturedformations with multiple length scales[J].Water Resources Research,2001,37(3):443-455):
a)油藏中油-水两相渗流微分方程:
Pfc=Pfo-Pfl (20)
式中:KFro——人工裂缝液相和油相的相对渗透率,无量纲;
Kfrl、Kfro——微裂缝液相和油相的相对渗透率,无量纲;
Kf——微裂缝的渗透率,D;
Vf——微裂缝网格的体积,m3;
qFl、qFo——人工裂缝中液相和油相的源汇项,m3/s;
SFl、Sfl——人工裂缝和微裂缝中液相饱和度,无量纲;
φF、φf、φm——人工裂缝、微裂缝和基质的孔隙度,无量纲;
Pfl、Pfo——微裂缝的液相压力和油相压力,MPa;
QfFl、QfFo——主裂缝的液相、油相窜流量,m3/s;
δf——微裂缝网格是否含有人工裂缝的判断参数,由于有人工裂缝穿过的微裂缝网格会发生微裂缝-人工裂缝的流体交换,因此当微裂缝网格有裂缝穿过时δ=1;当微裂缝网格无裂缝穿过时δ=0。
t2——油井返排-生产的时间,s;
β——单位转换系数,取β=0.001;
b)初始条件包括初始压力及初始饱和度的分布,即:
c)内边界条件为:
PF(xw,yw,t2)=Pwf(t2) (23)
式中,xw、yw——油井所在网格单元的横、纵坐标值,m;
Pwf——井底流压,MPa。
d)外边界条件为:
将步骤(2)中计算得到的参数当作上述油井返排-生产模型的初始条件,可以数值求解得到不同时间下的人工裂缝流体压力PF、微裂缝油相压力Pfo和基质油相压力Pmo;以及不同时间下的人工裂缝液相饱和度SFl、微裂缝液相饱和度Sfl和基质液相饱和度Sml。由此可以获得时间t2下基质网格的液相饱和度Sml(i,j,t2)、基质网格的油相压力Pmo(i,j,t2)、微裂缝网格的液相饱和度Sfl(i,j,t2)、微裂缝网格的油相压力Pfo(i,j,t2)、人工裂缝单元的液相饱和度SFl(L,t2)、人工裂缝单元的油相压力PF(L,t2)。
(4)计算油井的产能。
将时间t2=0和时间t2下基质网格液相饱和度、微裂缝网格液相饱和度、人工裂缝单元液相饱和度代入下式,可计算得到开井生产至时间t2时的油井累积产量Q。
式中,Q——开井生产至时间t2时的油井累积产量,m3;
ni,nj——基质和微裂缝网格中x方向上和y方向上的网格总数;
nL——人工裂缝单元的总数;
ξL——每个人工裂缝单元的长度,m;
xi,j、yi,j——i,j位置处基质和微裂缝网格的长度和宽度,m;
Sml(i,j,0)——油井返排-生产开始时i,j位置处基质网格的初始液相饱和度;
Sfl(i,j,0)——油井返排-生产开始时i,j位置处微裂缝网格的初始液相饱和度;
SFl(L,0)——油井返排-生产开始时第L段人工裂缝单元的初始液相饱和度;
Sml(i,j,t2)——开井生产至时间t2时i,j位置处基质网格的液相饱和度;
Sfl(i,j,t2)——开井生产至时间t2时i,j位置处微裂缝网格的液相饱和度;
SFl(L,t2)——开井生产至时间t2时第L段人工裂缝单元的液相饱和度。
附图说明
图1为G12井压裂施工模拟结束时的裂缝周围油相压力分布图。
图2为G12井生产60天后的储层油相压力分布图。
图3为G12井生产80天后的储层油相压力分布图。
图4为G12井累积产量随时间变化的曲线图。
具体实施方式
下面根据附图和实例进一步说明本发明,以便于本技术领域的技术人员理解本发明。但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,均在保护之列。
实施例
以华北地区某区块的油井(G12井)为例运用本发明的压采一体化数值模拟方法进行压裂、返排、生产过程的模拟,并绘制模拟过程中的裂缝扩展压力云图、返排-生产过程中的储层压力云图用以直观展示模拟效果。对G12井的计算步骤如下。
步骤1,建立结构化储层网格和裂缝初始单元。首先,在x-y直角坐标系下建立储层的结构化网格,模拟的储层长度Lx=400m,且划分成ni=200段;储层宽度Ly=400m,划分为nj=200段,整个储层便被划分成了100×100的结构化网格。随后,在储层网格基础上,建立人工裂缝的初始单元。以裂缝延伸位置为中心,取50m×100m的矩形范围内的储层网格进行细化,细化后的网格长度和宽度均为0.05m。之后在裂缝延伸位置添加一段延伸y方向有3个网格长度的线段作为人工裂缝的初始单元。此时,人工裂缝单元的总数nL=3,每个人工裂缝单元长度ξL=0.05m。
步骤2,基于已建立的网格,运用压裂模型公式(1)~(10)数值计算得到第1时间步下的人工裂缝单元的宽度WL,1、裂缝单元内的流体压力PFL,1、基质网格的油相压力Pmoi,j,1、基质网格的液相饱和度Smli,j,1,以及第1时间步下的尖端裂缝单元内的流体压力之后通过公式(11)可计算获得第1时间步下的裂缝尖端应力强度因子KIi,1,将其与储层岩石的断裂韧性KIC进行比较,若KIi,1≤KIC,则在第1时间步下,人工裂缝不发生裂缝扩展,裂缝长度不变,重复步骤2进行第2时间步的计算,第2时间步的人工裂缝单元总数仍为nL;若KIi,1>KIC,则在第1时间步下人工裂缝会沿y方向发生扩展,增加1个人工裂缝单元数,重复步骤2进行第2时间步的计算,此时人工裂缝单元总数为nL=nL+1。
计算直到压裂施工完成时t=tend结束,此时可获得时间t=tend时的人工裂缝单元总数nL、每个裂缝单元的宽度每个裂缝单元内的流体压力人工裂缝半长nL×ξL、每个基质网格的油相压力和每个基质网格的液相饱和度
图1为裂缝扩展模拟结束时的裂缝周围压力分布图。
步骤3,将步骤(2)中计算得到压裂施工完成时的参数当作油井返排-生产模型的初始条件。运用公式(12)~(24)可计算得到,不同时间步下的基质网格的液相饱和度Sml(i,j,t2)、基质网格的油相压力Pmo(i,j,t2)、微裂缝网格的液相饱和度Sfl(i,j,t2)、微裂缝网格的油相压力Pfo(i,j,t2)、人工裂缝单元的液相饱和度SFl(L,t2)、人工裂缝单元的流体压力PF(L,t2)。
图2、图3分别是油井返排-生产模型中不同时间下的压力分布图。
步骤4,将步骤(3)中得到的不同时间下的基质、微裂缝、人工裂缝的液相饱和度代入公式(25),便可计算得到不同时间下的油井累积产量Q。由此可以绘制累积产量随时间变化的曲线图(见图4)。
Claims (7)
1.一种油藏流动全耦合的压采一体化数值模拟方法,依次包括以下步骤:
(1)建立结构化储层网格和裂缝初始单元;
(2)计算压裂过程中的人工裂缝半长、人工裂缝宽度分布、人工裂缝内流体压力分布、基质油相压力分布和基质液相饱和度分布;
(3)计算油井返排-生产过程中储层压力分布和储层含水饱和度分布;
(4)计算油井的产能。
2.如权利要求1所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述步骤(1)过程如下:在x-y直角坐标系下建立储层的结构化网格,将储层长度Lx划分成ni段,储层宽度Ly划分为nj段,整个储层便被划分成了ni×nj的结构化网格,xi,j和yi,j分别代表着每个网格的长度和宽度,下标i和j代表着每一个网格在储层中的位置;在储层网格基础上,建立人工裂缝的初始单元,对压裂延伸位置处一定范围内的储层网格进行细化,并在细化后的网格上添加一段有3个网格长度的人工裂缝,裂缝延伸方向为y方向,将其作为人工裂缝初始单元,储层网格中人工裂缝被基质网格切分为具有一定长度的线段,每一个线段代表一个人工裂缝单元,人工裂缝的初始单元数为nL0,人工裂缝单元总数表示为nL,任意裂缝单元编号为L,每个人工裂缝单元长度为ξL,随着裂缝的延伸扩展,nL会逐渐增大,直至压裂结束。
3.如权利要求1所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述步骤(2)过程如下:
水力压裂模型中人工裂缝内为单相流动,水力压裂模型如下:
a)考虑流体滤失的裂缝宽度模型:
式中,W(y)——人工裂缝任意位置处宽度,m;
H——人工裂缝高度,m;
v——储层岩样的泊松比,无因次;
E——储层岩样的杨氏模量,MPa;
σh——水平最小主应力,MPa;
Pm(y)——人工裂缝任意位置处的基质流体压力,MPa。
μ——压裂液的黏度,mPa·s;
vl——压裂液滤失速度,m/s;
t——时间,s;
kmf——人工裂缝与基质的平均渗透率,mD;
b)压裂过程中的基质渗流模型
Pmc=Pmo-Pml
式中,Amf——裂缝与基质的接触面积,m;
KF——人工裂缝的渗透率,D;
VF——人工裂缝单元的体积,m3;
KFrl——人工裂缝液相相对渗透率,无量纲;
PF——人工裂缝流体压力,MPa;
φF——人工裂缝的孔隙度,无因次;
φm——储层基质孔隙度,无因次;
Km——基质渗透率,mD;
Kmrl——基质中液相的相对渗透率,无因次;
Kmro——基质中油相的相对渗透率,无因次;
Sml——基质中的液相饱和度,无因次;
Vb——基质单元体积,m3;
μl——基质中的液相黏度,mPa·s;
μo——基质中的油相黏度,mPa·s;
Bl——基质中液相的体积系数,无因次;
Bo——基质中油相的体积系数,无因次;
Pml、Pmo——基质中液相、油相的压力,MPa;
Pmc——基质中的毛管压力,MPa;
δm——基质网格是否含有人工裂缝的判断参数,当储层基质网格有裂缝穿过时δm=1;当储层基质网格无裂缝穿过时δm=0;
c)油藏渗流初始条件为:
Pmo(i,j,t)|t=0=Pe
式中,Pe——油藏原始地层压力,MPa;
i,j——网格的位置坐标;
d)裂缝延伸边界条件为:
PFtL=1,t=PZ
式中,G——储层岩样的体积模量,MPa;
PFL=1,t——压裂施工中的井底流体压力,即第1段人工裂缝单元内的流体压力,MPa;
PZ——压裂施工中的井底泵注压力,MPa;
nL——压裂施工时间t下的人工裂缝单元的总数;
ξL——人工裂缝单元的长度,m;
e)油藏基质渗流边界条件为:
式中,Lx,Ly——分别表示储层长度和储层宽度,m;
对水力压裂模型进行数值求解,得到任意时间下的人工裂缝任意位置处的宽度W(y)、人工裂缝任意位置处的流体压力PF、基质油相压力Pmo、基质液相饱和度Sml,由此得到压裂施工时间t下的每个裂缝单元的宽度WL,t、每个裂缝单元的流体压力PFL,t、每个基质网格的油相压力Pmoi,j,t和每个基质网格的液相饱和度Smli,j,t,以及压裂施工时间t下的尖端裂缝单元内的流体压力将代入下式,得到压裂施工时间t下的裂缝尖端应力强度因子KIi,t:
将裂缝尖端应力强度因子KIi,t与储层岩石的断裂韧性KIC进行比较,当KIi,t≤KIC时,压裂施工时间t下的人工裂缝不发生裂缝扩展,裂缝长度不变,重复步骤(2)进行时间t+1下的计算,时间t+1下的人工裂缝单元总数仍为nL;当KIi,t>KIC时,人工裂缝沿y方向发生扩展,裂缝长度增加,重复步骤(2)进行时间t+1下的计算,时间t+1下的人工裂缝单元总数为nL=nL+1;
4.如权利要求1所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述步骤(3)过程如下:
油井返排-生产模型如下:
a)油藏中油-水两相渗流微分方程:
Pfc=Pfo-Pfl
式中:KFro——人工裂缝液相和油相的相对渗透率,无量纲;
Kfrl、Kfro——微裂缝液相和油相的相对渗透率,无量纲;
Kf——微裂缝的渗透率,D;
Vf——微裂缝网格的体积,m3;
qFl、qFo——人工裂缝中液相和油相的源汇项,m3/s;
SFl、Sfl——人工裂缝和微裂缝中液相饱和度,无量纲;
φF、φf、φm——人工裂缝、微裂缝和基质的孔隙度,无量纲;
Pfl、Pfo——微裂缝的液相压力和油相压力,MPa;
QfFl、QfFo——主裂缝的液相、油相窜流量,m3/s;
δf——微裂缝网格是否含有人工裂缝的判断参数,当微裂缝网格有裂缝穿过时δ=1;当微裂缝网格无裂缝穿过时δ=0;
t2——油井返排-生产的时间,s;
β——单位转换系数,取β=0.001;
b)初始条件包括初始压力及初始饱和度的分布:
c)内边界条件为:
PF(xw,yw,t2)=Pwf(t2)
式中,xw、yw——油井所在网格单元的横、纵坐标值,m;
Pwf——井底流压,MPa;
d)外边界条件为:
将步骤(2)中计算得到的参数当作油井返排-生产模型的初始条件,求解得到不同时间下的人工裂缝流体压力PF、微裂缝油相压力Pfo和基质油相压力Pmo;以及不同时间下的人工裂缝液相饱和度SFl、微裂缝液相饱和度Sfl和基质液相饱和度Sml,由此获得时间t2下基质网格的液相饱和度Sml(i,j,t2)、基质网格的油相压力Pmo(i,j,t2)、微裂缝网格的液相饱和度Sfl(i,j,t2)、微裂缝网格的油相压力Pfo(i,j,t2)、人工裂缝单元的液相饱和度SFl(L,t2)、人工裂缝单元的油相压力PF(L,t2)。
5.如权利要求1所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述步骤(4)过程如下:
将时间t2=0和时间t2下基质网格液相饱和度、微裂缝网格液相饱和度、人工裂缝单元液相饱和度代入下式,计算得到开井生产至时间t2时的油井累积产量Q:
式中,Q——开井生产至时间t2时的油井累积产量,m3;
ni,nj——基质和微裂缝网格中x方向上和y方向上的网格总数;
nL——人工裂缝单元的总数;
ξL——每个人工裂缝单元的长度,m;
xi,j、yi,j——i,j位置处基质和微裂缝网格的长度和宽度,m;
Sml(i,j,0)——油井返排-生产开始时i,j位置处基质网格的初始液相饱和度;
Sfl(i,j,0)——油井返排-生产开始时i,j位置处微裂缝网格的初始液相饱和度;
SFl(L,0)——油井返排-生产开始时第L段人工裂缝单元的初始液相饱和度;
Sml(i,j,t2)——开井生产至时间t2时i,j位置处基质网格的液相饱和度;
Sfl(i,j,t2)——开井生产至时间t2时i,j位置处微裂缝网格的液相饱和度;
SFl(L,t2)——开井生产至时间t2时第L段人工裂缝单元的液相饱和度。
6.如权利要求3所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述储层岩石的断裂韧性KIC=2MPa·m1/2。
7.如权利要求4所述的一种油藏流动全耦合的压采一体化数值模拟方法,其特征在于,所述油井返排-生产模型运用双重介质模型和嵌入式离散裂缝模型进行建立,考虑人工裂缝、微裂缝、基质、人工裂缝-微裂缝之间和基质-微裂缝之间的流体渗流,人工裂缝中忽略毛管力作用,液相压力与油相压力相同,且均等于人工裂缝单元内的流体压力PF。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110422759.9A CN113111607B (zh) | 2021-04-15 | 2021-04-15 | 一种油藏流动全耦合的压采一体化数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110422759.9A CN113111607B (zh) | 2021-04-15 | 2021-04-15 | 一种油藏流动全耦合的压采一体化数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113111607A true CN113111607A (zh) | 2021-07-13 |
CN113111607B CN113111607B (zh) | 2022-04-15 |
Family
ID=76718641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110422759.9A Active CN113111607B (zh) | 2021-04-15 | 2021-04-15 | 一种油藏流动全耦合的压采一体化数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113111607B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114547998A (zh) * | 2022-02-28 | 2022-05-27 | 西南石油大学 | 一种通过耦合储层流动确定水平井压裂改造体积的方法 |
CN116084906A (zh) * | 2023-03-21 | 2023-05-09 | 西南石油大学 | 一种小规模砂体油藏增能提采工艺 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103953323A (zh) * | 2014-05-08 | 2014-07-30 | 西南石油大学 | 一种水平井产生缝网的水力压裂工艺 |
CN105243182A (zh) * | 2015-09-08 | 2016-01-13 | 中国石油天然气股份有限公司 | 致密油压裂水平井的动态储量计算方法 |
CN110863810A (zh) * | 2019-11-21 | 2020-03-06 | 西南石油大学 | 一种耦合页岩气藏水力压裂返排生产过程一体化模拟方法 |
CN111651912A (zh) * | 2020-05-09 | 2020-09-11 | 大连理工大学 | 一种基于断裂力学准则的水力裂缝扩展路径数值计算方法 |
AU2020102430A4 (en) * | 2020-09-04 | 2020-11-05 | University Of Science And Technology Beijing | Method and device for identifying structure of flow field in unconventional gas reservoir |
US20200387650A1 (en) * | 2019-06-10 | 2020-12-10 | Southwest Petroleum University | Fracturing fluid flow-back simulation method for fractured horizontal well in shale gas reservoir |
-
2021
- 2021-04-15 CN CN202110422759.9A patent/CN113111607B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103953323A (zh) * | 2014-05-08 | 2014-07-30 | 西南石油大学 | 一种水平井产生缝网的水力压裂工艺 |
CN105243182A (zh) * | 2015-09-08 | 2016-01-13 | 中国石油天然气股份有限公司 | 致密油压裂水平井的动态储量计算方法 |
US20200387650A1 (en) * | 2019-06-10 | 2020-12-10 | Southwest Petroleum University | Fracturing fluid flow-back simulation method for fractured horizontal well in shale gas reservoir |
CN110863810A (zh) * | 2019-11-21 | 2020-03-06 | 西南石油大学 | 一种耦合页岩气藏水力压裂返排生产过程一体化模拟方法 |
CN111651912A (zh) * | 2020-05-09 | 2020-09-11 | 大连理工大学 | 一种基于断裂力学准则的水力裂缝扩展路径数值计算方法 |
AU2020102430A4 (en) * | 2020-09-04 | 2020-11-05 | University Of Science And Technology Beijing | Method and device for identifying structure of flow field in unconventional gas reservoir |
Non-Patent Citations (5)
Title |
---|
DAPENG GAO 等: "Coupled numerical simulation of multi-layer reservoir developed by lean-stratified water injection", 《JOURNAL OF PETROLEUM EXPLORATION AND PRODUCTION TECHNOLOGY》 * |
任冬梅: "考虑井眼周围应力变化的油藏—井筒耦合的一体化模拟研究", 《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅰ辑》 * |
吕坐彬 等: "基于非结构化网格的离散裂缝建模技术研究", 《重庆科技学院学报(自然科学版)》 * |
张建光 等: "页岩气藏水力压裂中应力-流压耦合效应及人工裂缝扩展规律", 《中国石油大学学报(自然科学版)》 * |
郭小哲 等: "页岩气藏压裂水平井物质平衡模型的建立", 《西南石油大学学报(自然科学版)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114547998A (zh) * | 2022-02-28 | 2022-05-27 | 西南石油大学 | 一种通过耦合储层流动确定水平井压裂改造体积的方法 |
CN114547998B (zh) * | 2022-02-28 | 2024-03-22 | 西南石油大学 | 一种通过耦合储层流动确定水平井压裂改造体积的方法 |
CN116084906A (zh) * | 2023-03-21 | 2023-05-09 | 西南石油大学 | 一种小规模砂体油藏增能提采工艺 |
CN116084906B (zh) * | 2023-03-21 | 2023-10-31 | 西南石油大学 | 一种小规模砂体油藏增能提采工艺 |
Also Published As
Publication number | Publication date |
---|---|
CN113111607B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111322050B (zh) | 一种页岩水平井段内密切割暂堵压裂施工优化方法 | |
RU2723769C1 (ru) | Способ расчета объема обратного потока жидкости для гидроразрыва пласта при гидроразрыве в горизонтальных скважинах в залежах газа трещиноватых песчаников | |
CN107044277B (zh) | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 | |
CN112561144B (zh) | 一种致密油压裂水平井产能主控因素评判与产能预测方法 | |
CN113111607B (zh) | 一种油藏流动全耦合的压采一体化数值模拟方法 | |
CN105952432B (zh) | 一种超低渗致密油藏体积压裂水平井准自然能量开发交错布井方法 | |
CN107506948B (zh) | 一种基于动态泄流体积的页岩油气综合产量分析方法 | |
CN105354639B (zh) | 致密油多重介质耦合渗流的全周期产能预测方法及其装置 | |
CN111980654B (zh) | 一种非均匀页岩油藏分段压裂水平井产能计算方法 | |
CN110685657B (zh) | 一种转向压裂用暂堵颗粒用量计算方法 | |
CN107066769B (zh) | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 | |
CN104594872A (zh) | 一种优化致密气藏压裂水平井裂缝导流能力的方法 | |
CN107630686B (zh) | 水平井段间交错驱替和渗吸驱油的致密油能量补充方法 | |
CN109408859A (zh) | 页岩气藏压裂水平井二维三重介质数值模型建立方法 | |
CN108240214A (zh) | 一种页岩气藏水平井分段压裂压后产能预测方法 | |
CN110259421B (zh) | 一种裂缝性的致密油藏注水补充能量方法 | |
CN112541287A (zh) | 疏松砂岩压裂充填防砂增产调剖一体化设计方法 | |
CN116894572B (zh) | 一种超深井考虑岩崩后出砂的合理配产方法 | |
CN108316915B (zh) | 一种确定油气井致密储层中纤维暂堵转向液最优用量的方法 | |
CN116108572A (zh) | 页岩凝析气井体积压裂外区产能贡献分析方法 | |
CN107247816B (zh) | 一种稠油油藏开展冷采化学降粘最远施工半径的判定方法 | |
CN114547998B (zh) | 一种通过耦合储层流动确定水平井压裂改造体积的方法 | |
CN110671088A (zh) | 一种考虑固相封堵主裂缝的新裂缝起裂压力预测方法 | |
CN113803042B (zh) | 一种单段单簇密集压裂方法及*** | |
CN113417616B (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 |