CN113743037A - 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 - Google Patents
一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 Download PDFInfo
- Publication number
- CN113743037A CN113743037A CN202111081319.8A CN202111081319A CN113743037A CN 113743037 A CN113743037 A CN 113743037A CN 202111081319 A CN202111081319 A CN 202111081319A CN 113743037 A CN113743037 A CN 113743037A
- Authority
- CN
- China
- Prior art keywords
- fracture
- variable
- crack
- pressure
- water injection
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 91
- 238000002347 injection Methods 0.000 title claims abstract description 70
- 239000007924 injection Substances 0.000 title claims abstract description 70
- 238000004364 calculation method Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000003860 storage Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 18
- 230000000694 effects Effects 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000007704 transition Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 9
- 230000035699 permeability Effects 0.000 claims description 22
- 239000000243 solution Substances 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 210000003491 skin Anatomy 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 210000002615 epidermis Anatomy 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 4
- 206010017076 Fracture Diseases 0.000 description 124
- 208000010392 Bone Fractures Diseases 0.000 description 108
- 238000011161 development Methods 0.000 description 11
- 239000003129 oil well Substances 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000033228 biological regulation Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 208000013201 Stress fracture Diseases 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000004200 deflagration Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- UUQHAAWMCLITRZ-KEOOTSPTSA-N 2-[(3S,6S,12S,20R,23S)-20-carbamoyl-12-[3-(diaminomethylideneamino)propyl]-3-(1H-indol-3-ylmethyl)-2,5,8,11,14,22-hexaoxo-17,18-dithia-1,4,7,10,13,21-hexazabicyclo[21.3.0]hexacosan-6-yl]acetic acid Chemical compound N1C(=O)[C@H](CC(O)=O)NC(=O)CNC(=O)[C@H](CCCN=C(N)N)NC(=O)CCSSC[C@@H](C(N)=O)NC(=O)[C@@H]2CCCN2C(=O)[C@@H]1CC1=CNC2=CC=CC=C12 UUQHAAWMCLITRZ-KEOOTSPTSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000009958 sewing Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010200 validation analysis 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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- 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)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明具体涉及一种针对低渗透油藏的注水诱导动态裂缝变导流能力计算方法。一种低渗透油藏注水诱导动态裂缝变导流能力计算方法,包括步骤1:建立物理模型;步骤2:将整个***分为井筒、裂缝和基质3部分;裂缝‑变缝长&变导流能力模型典型曲线包括5个流动阶段:井筒储集段;过渡段;变裂缝储集段;第二过渡段;受裂缝属性参数变化影响的线性流段。本发明针对低渗透油藏注水诱导动态裂缝问题,揭示了低于破裂压力条件下注水诱导裂缝的机理,建立了考虑关井测压过程中裂缝属性参数变化的动态诱导裂缝压力分析模型,揭示了井筒存储和由于裂缝闭合导致的“裂缝存储”效应,可快速准确解释注水诱发裂缝长度及导流能力等重要参数。
Description
技术领域
本发明属于石油开采技术领域,具体涉及一种针对低渗透油藏的注水诱导动态裂缝变导流能力计算方法。
背景技术
裂缝普遍存在于低渗透砂岩中,国内大量低渗透砂岩油藏的共同特征是,油层致密,水动力联系差,油井无自然产能,需经压裂改造后才具有产能,油藏压力低,注水保持压力成为主要开发方式。但在注水开发过程中,低渗透砂岩油藏面临的突出问题是,局部裂缝发育区注水时,注入水沿裂缝突进、引起水窜和水淹,而裂缝不发育区储层物性差,导致注入压力高、水井欠注,水驱波及体积小,储量动用程度低,裂缝对油田注水开发的影响越来越大,控制水驱开发效果及剩余油分布规律的关键因素,通过对井网与注水井动态裂缝进行匹配优化,可以最大限度地发挥油藏潜力。因此正确认识动态裂缝特征成为这类油田科学开发的重要依据。
低渗透油藏在长期注水开发过程中,由于注水井近井地带憋压,当井底压力超过岩层破裂、延伸压力,导致岩层破裂,或原始状态下闭合、充填的天然裂缝被激动、复活,产生了新生、有效裂缝通道。这些裂缝受现今地应力场控制,随着注水量的增长和井底压力的升高,不断向油井方向延展,直至与油井压裂缝连通,这些新生、有效裂缝称之为动态裂缝,动态裂缝按成因分成3种类型:天然潜在闭合型动态裂缝、超破裂压力型动态裂缝以及人工措施诱导型动态裂缝。动态裂缝的开启是受现今水平地应力和天然裂缝双重因素控制的,注水开发中裂缝开启顺序取决于天然裂缝走向与现今最大水平主应力方向的夹角,夹角越小,天然裂缝开启压力就越小,注入水将首先沿现今最大水平主应力方向的天然裂缝推进,随着注水压力的增大,其它方向的天然裂缝将依次开启,导致油藏单方向或多方向水淹,油井产量递减幅度大,油田开发效果变差。
低渗透油藏注水开发过程中,注水井和采油井的流量变化差异大,注水井虽未压裂或只进行了小型压裂,但注水量大;采油井即使进行了大规模压裂,但产液量仍然很低,同时油井见水具有单方向性,这些现象与注水诱导微细裂缝形成水淹通道密切相关。国内目前对压力的解释模型,主要使用复合模型、均质模型来解释,不能用于注水诱导裂缝情况。
发明内容
本发明旨在针对上述问题,提出一种针对低渗透油藏的注水诱导动态裂缝变导流能力计算方法。
本发明的技术方案在于:
(一)注水诱导动态裂缝机理
(1)天然裂缝被激活
随着注水的进行,地层压力逐渐升高,当达到天然微裂缝的开启压力时,将诱导原先呈闭合状态的天然裂缝,不同走向的裂缝开启压力与裂缝产状与最大主应力方向上的夹角有关,若天然裂缝为垂直缝,缝面在现今应力条件下所受到的应力为:
式中:σx和σy为平行、垂直于裂缝面所受到的应力,MPa;σmax和σmin为最大、最小水平主应力,MPa;β为裂缝走向与最大主应力方向的夹角,θ;τxy为裂缝尖端的切应力,MPa;天然裂缝的开启主要是由垂直于裂缝面上的张应力引起的,当缝内压力超过最大张应力时,裂缝将会开启,注水诱导动态裂缝沿着最大主应力方向延伸。
(2)天然裂缝被激活
低渗油藏由于储层渗透能力差,吸水能力差,为增强注入能力,注水井采用***压裂、爆燃及复合射孔投注方式,在近井地带形成多条径向、短小的裂缝。裂缝在注水过程中相互沟通和延伸,最终形成动态裂缝渗流通道。裂缝的延伸压力可由下式计算:
式中:ptip为裂缝延伸压力,MPa;U为缝面能,J;ν为泊松比,无因次;E为弹性模量,MPa;xf为裂缝半长,cm。
(3)地层破裂形成新缝由于注水压力超过地层破裂压力而形成的注水诱导裂缝,主要适用于天然裂缝不发育的储层,注水诱导动态裂缝的条件为:
ptip=3σmin-σmax+σt-pp (5)
式中:σt为岩石抗张强度,MPa;pp为地层孔隙压力;MPa;
若储层中天然微裂缝不发育,未经人工爆燃、***压裂,注入水仅能依靠基质储层吸收,由于低渗透储层吸水能力弱,持续的注水使得井底附近压力不断升高,当储层压力达到岩石破裂压力时,将会破裂地层产生新缝。说明低渗油藏诱导裂缝的主要机理不是地层破裂形成新缝,而是激活了原先呈闭合状态的微细缝,然后不断延伸形成裂缝条带,因此将地层破裂压力作为注水压力的上限不足以防止诱导裂缝的产生。
(二)注水诱导动态裂缝矿场特征
注水诱导裂缝是一个动态的形成过程,它是由于长期注水过程中,由于注水井近井地带憋压,诱导了原先呈闭合状态的微裂缝、或破裂地层产生了新缝。注水诱导裂缝的形成有以下几个矿场特征;
(1)油井含水快速上升、对应的水井注入压力对注水量不敏感;
(2)注水指示曲线存在明显拐点;
(3)注水井吸水剖面呈指状/尖峰状吸水特征;
(4)注水井试井曲线表现出裂缝特征、解释的有效渗透率远大于基质渗透率;
(5)油井见效及水淹呈明显的方向性特征;
(6)示踪剂、或水驱前缘等测试结果显示方向性连通。
(三)裂缝变缝长/变导流能力模型
步骤1:建立物理模型
(1)储层均质、等厚、无限大,上下边界封闭,忽略温度的影响;
(2)仅考虑单相水的渗流,假设其满足达西规律;
(3)仅考虑一条注水诱导缝的影响,假设其贯穿整个油层,并沿井筒对称;
(4)考虑关井测压阶段由于裂缝闭合引起的裂缝导流系数及裂缝半长的变化;
(5)考虑裂缝为有限导流,由于低渗储层渗透性差,假设水体在储层中渗流满足双线性流态;忽略重力的作用;
将整个***分为井筒、裂缝和基质3部分;
裂缝-变缝长&变导流能力模型典型曲线包括5个流动阶段:井筒储集段;过渡段;变裂缝储集段;第二过渡段;受裂缝属性参数变化影响的线性流段;
井底压力解pFD的计算过程为:
式中:pFD为裂缝处压力,无因次;
γD为渗透率变化系数,无因次;
ξFD为考虑裂缝导流系数变化的实空间解,无因次;(一)其中,γD的计算过程为:
式中:qw为注入量,cm3/d;
μ为流体粘度,cP;B为体积系数,cm3/cm3;
km为基质渗透率,无因次;
h为油藏厚度,cm;
γ为渗透率变化系数,MPa-1;
考虑关井时注水诱导裂缝渗透率的变化满足以下规律:
式中:kF某一压力下的渗透率,mD;
pF为裂缝处压力,MPa;
Ki为裂缝初始渗透率,mD;
Pi为原始地层压力,MPa;
PF为某时刻地层压力,MPa;
(二)其中,ξFD的计算过程为:
s为表皮系数,无因次;
CFD0为初始裂缝导流系数,无因次;
ε为中间变量,无因次;
其中:
式中:ηfD为裂缝导压系数,无因次;
φ为孔隙度,无因次;
ct为综合压缩系数,atm-1;
k为渗透率,mD;
kf为裂缝渗透率,mD;
w为裂缝宽度,cm;
km为基质渗透率,mD;
Lf0裂缝与井筒距离,cm;
考虑裂缝壁面的压力降,有:
ξwD2(tD)=ξFD'(tD)+Sfs (11)
式中:ξwD2为考虑裂缝壁面压力降的拉式空间解,无因次;
Sfs为裂缝面表皮系数,无因次;
关井后注水诱导动态裂缝闭合时,会将原先存储于裂缝中的水挤入地层,形成井筒续流现象,即为“裂缝储集”,则有
式中:ξwD3为考虑因裂缝闭合引起变流量作用的拉式空间解,无因次;
CfD为裂缝储集系数,无因次;
τD为无因次积分时间变量,无因次;
由于裂缝的储集系数不是常数,而是压力的函数,考虑裂缝-井筒间的额外压力损失,可得:
ξwD4(tD)=ξwD3(tD)+Schoke (13)
式中:ξwD4为考虑裂缝-井筒间额外压力损失的拉式空间解,无因次;
Schoke为井筒表皮系数,无因次;
根据杜哈美原理,考虑井筒储集的影响,则有:
式中:ξwD5为井筒储集压力,MPa;
CwD为井筒储集系数,无因次;
求得的ξwD5即为ξFD;
将γD及ξFD带入公式(10)计算得到井底压力解pFD,得到裂缝-变缝长&变导流能力曲线,并从该裂缝-变缝长&变导流能力曲线解释参数。
本发明的技术效果在于:
本发明针对低渗透油藏注水诱导动态裂缝问题,揭示了低于破裂压力条件下注水诱导裂缝的机理,建立了考虑关井测压过程中裂缝属性参数变化的动态诱导裂缝压力分析模型,揭示了井筒存储和由于裂缝闭合导致的“裂缝存储”效应,可快速准确解释注水诱发裂缝长度及导流能力等重要参数,为精细水驱实时调控提供科学依据。
附图说明
图1为裂缝-变缝长&变导流能力模型典型曲线。
图2为本模型与Spivey&Lee模型计算结果对比图。
图3为本模型与Saphir(Ecrin 2011)数值模型计算结果对比图。
图4为X1井双对数拟合对比图。
图5为X2井Saphir与新模型拟合对比图。
具体实施方式
实施例1-模型验证Spivey&Lee
输入参数如下表1所示
表1解析/数值验证的输入参数
(1)解析法验证:本专利提出的“裂缝变缝长/变导流能力模型”与Spivey&Lee模型作对比,忽略了由于裂缝长度变化引起的“变裂缝储集”效应的影响,并假设裂缝导流系数为定值。利用表1中的基础数据,通过本发明提供的“裂缝变缝长/变导流能力模型”计算得到井底压力解,形成裂缝-变缝长&变导流能力曲线;与Spivey&Lee模型进行结果对比,如图2所示;从图上可知,两种方法的结果表现出较好的一致性,验证了“裂缝变缝长/变导流能力模型”的准确性。
(2)数值法验证:Saphir(Ecrin 2011)
利用Saphir软件作试井设计,采用非结构化网格,模拟压裂井的井底压力响应特征,裂缝半长取300m,为了满足无限大油藏的特征,油藏大小取3000×2000m,采用局部网格加密,裂缝两端为长方形网格来模拟地层线性流。由于Saphir无法模拟动态裂缝的影响,因此将新模型中的裂缝设置为静态的,忽略裂缝闭合的影响,即不考虑裂缝长度和导流系数的变化,利用表1中的基础数据,通过本发明提供的“裂缝变缝长/变导流能力模型”计算得到井底压力解,形成裂缝-变缝长&变导流能力曲线;与Saphir模型进行结果对比,如图3所示;从图上可知,两种方法的结果表现出较好的一致性,验证了“裂缝变缝长/变导流能力模型”的准确性。
本专利建立了“裂缝-变缝长&变导流能力”数学模型,通过对模型数值求解得到井底压力解,绘制出不同时刻的压力随时间的变化关系曲线,得到了裂缝-变缝长&变导流能力模型典型曲线,该典型曲线包括5个流动阶段:井筒储集段、过渡段、变裂缝储集段、第二过渡段、受裂缝属性参数变化影响的线性流段。
由图3中可以看出,典型曲线具有两个储集阶段:一是真实的井筒储集,二是由于裂缝的闭合导致的裂缝储集效应,两者可用井筒-裂缝间表皮Schoke表征。当Schoke=0时,前三段(井筒储集、过渡段、裂缝储集段)将合并为一段,由于裂缝半长的逐渐减小,使得裂缝的储集系数减小,因此第三段会出现变裂缝储集的典型特征,表现为压力导数斜率大于1,同时动态裂缝的闭合会导致导流系数的减小,第五线性段则表现为压力导数斜率大于1/2的特征。该典型曲线创新的解释了注水动态裂缝的变化引起的裂缝储集效应,有效地获取动态裂缝参数,为精细注水调控提供依据。
具体应用例1
X1井是延长西部油田长6层位的一口注水井,日注水量30m3/d,采取井口关井压力降落测试,测得的井底压力解的裂缝-变缝长&变导流能力曲线如图4所示;与图1裂缝-变缝长&变导流能力模型典型曲线相比,X1井只测到了前四个流动阶段:井筒储集阶段、过渡流阶段、裂缝储集阶段和第二过渡流阶段。且与压力历史曲线拟合程度较好;可以明显地看到压力导数有两个斜率为1的线,这是注水井典型的“双储集”特征。即:井筒储集和裂缝储集;
通过拟合出的裂缝储集系数,得到动态裂缝的裂缝半长,为后续精细注水调控提供重要的依据。本应用例中,解释出裂缝储集系数为8.94m3/MPa,动态裂缝的裂缝半长为49.8m,解释结果与该井的生产动态比较符合,认为拟合分析结果是可靠的,正确的。
具体应用例2
X2井该井关井前的日注水量为25m3/d,同样采取井口关井压力降落的测试方式。其压力特征表现为,初期压力和压力导数相互重合呈斜率为1的线,随后压力导数上翘,超过了对应的压力曲线,然后压力导数下降,紧接着又发生明显的上翘。
首先利用Saphir软件,选择变井储、一条封闭边界模型进行解释,拟合结果如下图5(a)所示,可以明显看出,理论曲线与实测压力数据的拟合度较低,变储集现象和后期的压力导数上翘并没拟合上;同时解释出的储集系数大,为9.3m3/MPa,表皮系数为-9.74,表明注水井附近可能诱导了动态裂缝。
同时利用“裂缝-变缝长&变导流能力”模型解释该井,如下图5(b)所示,与Saphir软件结果对比可以看出,新模型能较好拟合变井储效应和后期受裂缝闭合影响的线性流,解释出裂缝储集系数29.17m3/MPa,得到裂缝半长约为100m,试井拟合分析法得到的裂缝半长约为110m,两者结果基本一致,验证了新模型结果的可靠性、正确性。说明了注水井由于诱导裂缝闭合的“变井储效应”和后期的压力导数上翘,这个变储集效应是由于关井测压过程中裂缝半长的“缩短”而导致裂缝储集系数减小,而后期压力导数上翘是因为裂缝闭合引起其导流系数减小。
Claims (1)
1.一种低渗透油藏注水诱导动态裂缝变导流能力计算方法,其特征在于:
步骤1:建立物理模型
(1)储层均质、等厚、无限大,上下边界封闭,忽略温度的影响;
(2)仅考虑单相水的渗流,假设其满足达西规律;
(3)仅考虑一条注水诱导缝的影响,假设其贯穿整个油层,并沿井筒对称;
(4)考虑关井测压阶段由于裂缝闭合引起的裂缝导流系数及裂缝半长的变化;
(5)考虑裂缝为有限导流,由于低渗储层渗透性差,假设水体在储层中渗流满足双线性流态;忽略重力的作用;
步骤2:
将整个***分为井筒、裂缝和基质3部分;
裂缝-变缝长&变导流能力模型典型曲线包括5个流动阶段:井筒储集段;过渡段;变裂缝储集段;第二过渡段;受裂缝属性参数变化影响的线性流段;
井底压力解pFD的计算过程为:
式中:pFD为裂缝处压力,无因次;
γD为渗透率变化系数,无因次;
ξFD为考虑裂缝导流系数变化的实空间解,无因次;
(一)其中,γD的计算过程为:
式中:qw为注入量,cm3/d;
μ为流体粘度,cP;B为体积系数,cm3/cm3;
km为基质渗透率,无因次;
h为油藏厚度,cm;
γ为渗透率变化系数,MPa-1;
考虑关井时注水诱导裂缝渗透率的变化满足以下规律:
式中:kF为某一压力下的渗透率,mD;
pF为裂缝处压力,MPa;
Ki为裂缝初始渗透率,mD;
Pi为原始地层压力,MPa;
(二)其中,ξFD的计算过程为:
s为表皮系数,无因次;
CFD0为初始裂缝导流系数,无因次;
ε为中间变量,无因次;
其中:
式中:ηfD为裂缝导压系数,无因次;
φ为孔隙度,无因次;
ct为综合压缩系数,atm-1;
k为渗透率,mD;
kf为裂缝渗透率,mD;
w为裂缝宽度,cm;
km为基质渗透率,mD;
Lf0裂缝与井筒距离,cm;
ξwD2(tD)=ξFD'(tD)+Sfs (11)
式中:ξwD2为考虑裂缝壁面压力降的拉式空间解,无因次;
Sfs为裂缝面表皮系数,无因次;
关井后注水诱导动态裂缝闭合时,会将原先存储于裂缝中的水挤入地层,形成井筒续流现象,即为“裂缝储集”,则有
式中:ξwD3为考虑因裂缝闭合引起变流量作用的拉式空间解,无因次;
CfD为裂缝储集系数,无因次;
τD为无因次积分时间变量,无因次;
由于裂缝的储集系数不是常数,而是压力的函数,考虑裂缝-井筒间的额外压力损失,可得:
ξwD4(tD)=ξwD3(tD)+Schoke (13)
式中:ξwD4为考虑裂缝-井筒间额外压力损失的拉式空间解,无因次;
Schoke为井筒表皮系数,无因次;
根据杜哈美原理,考虑井筒储集的影响,则有:
式中:ξwD5为井筒储集压力,MPa;
CwD为井筒储集系数,无因次;
求得的ξwD5即为ξFD;
将γD及ξFD带入公式(10)计算得到井底压力解pFD,得到裂缝-变缝长&变导流能力曲线,并从该裂缝-变缝长&变导流能力曲线解释参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111081319.8A CN113743037B (zh) | 2021-09-15 | 2021-09-15 | 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111081319.8A CN113743037B (zh) | 2021-09-15 | 2021-09-15 | 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113743037A true CN113743037A (zh) | 2021-12-03 |
CN113743037B CN113743037B (zh) | 2023-07-25 |
Family
ID=78739084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111081319.8A Active CN113743037B (zh) | 2021-09-15 | 2021-09-15 | 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113743037B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114834A (zh) * | 2022-07-29 | 2022-09-27 | 西南石油大学 | 一种复杂情况下的压裂井试井模拟方法 |
CN117307152A (zh) * | 2023-11-28 | 2023-12-29 | 四川省能源地质调查研究所 | 页岩气井压裂裂缝导流能力动态评估方法及其装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150211339A1 (en) * | 2014-01-24 | 2015-07-30 | King Fahd University Of Petroleum And Minerals | Quantification of skin in hydraulic fracturing of low and tight reservoirs |
CN107480383A (zh) * | 2017-08-21 | 2017-12-15 | 中国石油大学(北京) | 一种通过测压数据监测注水动态裂缝的方法 |
CA3020545A1 (en) * | 2017-10-13 | 2019-04-13 | Uti Limited Partnership | Completions for inducing fracture network complexity |
CN110159260A (zh) * | 2019-04-02 | 2019-08-23 | 中国石油大学(北京) | 用于裂缝部分闭合压裂直井主要来水方向判别方法及装置 |
CN111236908A (zh) * | 2020-01-09 | 2020-06-05 | 西南石油大学 | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 |
-
2021
- 2021-09-15 CN CN202111081319.8A patent/CN113743037B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150211339A1 (en) * | 2014-01-24 | 2015-07-30 | King Fahd University Of Petroleum And Minerals | Quantification of skin in hydraulic fracturing of low and tight reservoirs |
CN107480383A (zh) * | 2017-08-21 | 2017-12-15 | 中国石油大学(北京) | 一种通过测压数据监测注水动态裂缝的方法 |
CA3020545A1 (en) * | 2017-10-13 | 2019-04-13 | Uti Limited Partnership | Completions for inducing fracture network complexity |
CN110159260A (zh) * | 2019-04-02 | 2019-08-23 | 中国石油大学(北京) | 用于裂缝部分闭合压裂直井主要来水方向判别方法及装置 |
CN111236908A (zh) * | 2020-01-09 | 2020-06-05 | 西南石油大学 | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 |
Non-Patent Citations (1)
Title |
---|
赵向原;曾联波;靳宝光;王策;李恕军;: "低渗透油藏注水诱导裂缝特征及形成机理――以鄂尔多斯盆地安塞油田长6油藏为例", 石油与天然气地质, no. 04 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114834A (zh) * | 2022-07-29 | 2022-09-27 | 西南石油大学 | 一种复杂情况下的压裂井试井模拟方法 |
CN115114834B (zh) * | 2022-07-29 | 2024-02-23 | 西南石油大学 | 一种复杂情况下的压裂井试井模拟方法 |
CN117307152A (zh) * | 2023-11-28 | 2023-12-29 | 四川省能源地质调查研究所 | 页岩气井压裂裂缝导流能力动态评估方法及其装置 |
CN117307152B (zh) * | 2023-11-28 | 2024-02-09 | 四川省能源地质调查研究所 | 页岩气井压裂裂缝导流能力动态评估方法及其装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113743037B (zh) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | An experimental study on the mechanism of degradable fiber-assisted diverting fracturing and its influencing factors | |
Zeng et al. | Breakthrough in staged fracturing technology for deep shale gas reservoirs in SE Sichuan Basin and its implications | |
Zou et al. | 3-D numerical simulation of hydraulic fracturing in a CBM reservoir | |
US11408264B2 (en) | Volumetric fracturing method of temporarily plugging and diverting through functional slick water with oil displacement agent injected simultaneously | |
CN110344799B (zh) | 一种提高裂缝复杂性的临界砂堵压裂方法 | |
Shiozawa et al. | EGS designs with horizontal wells, multiple stages, and proppant | |
Zhu et al. | Coupled flow-stress-damage simulation of deviated-wellbore fracturing in hard-rock | |
CN107965305A (zh) | 一种分层重复压裂方法 | |
CN103089224A (zh) | 一种综合控制裂缝高度的压裂方法 | |
CN113743037A (zh) | 一种低渗透油藏注水诱导动态裂缝变导流能力计算方法 | |
WO2019014090A2 (en) | METHODS AND SYSTEMS FOR HYDRAULIC BALLOON FRACTURES AND COMPLEX END FLOODING | |
CN110259421A (zh) | 一种裂缝性的致密油藏注水补充能量方法 | |
CN112324412A (zh) | 一种体积压裂形成复杂缝网的方法 | |
Nguyen et al. | Fracture height containment by creating an artificial barrier with a new additive | |
Lu et al. | Numerical simulation on the initiation and propagation of synchronous perforating fractures in horizontal well clusters | |
US20150369032A1 (en) | Method for selecting the location of a simulating geothermal well | |
CN113743036B (zh) | 一种低渗-致密油藏注水诱导动态裂缝条带模型计算方法 | |
CN108121844A (zh) | 水力波及半径的获得方法 | |
CN115081352B (zh) | 深水高温高压疏松砂岩端部脱砂工艺设计优化方法及装置 | |
Huang et al. | Study on the coalbed methane development under high In Situ stress, large buried depth, and low permeability reservoir in the Libi Block, Qinshui Basin, China | |
Smith et al. | Enhanced 2d proppant-transport simulation: The key to understanding proppant flowback and post-frac productivity | |
Chen et al. | Chemical Agents Diversion with Microseismic Monitoring–New Prospects of Refracturing for Open-Hole Horizontal Well in Tight Oil Reservoir of Junggar Basin, China | |
CN114198077A (zh) | 一种水平井分段压裂提高远井裂缝复杂性的方法及其应用 | |
Di et al. | Study on the production mode of horizontal well from water huff and puff to segmented injection and production in tight reservoir | |
Bao et al. | An evaluation of enhanced oil recovery strategies for extra heavy oil reservoir after cold production without sand in Orinoco, Venezuela |
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 |