CN108868731A - 一种裂缝性储层酸压动态综合滤失系数的计算方法 - Google Patents
一种裂缝性储层酸压动态综合滤失系数的计算方法 Download PDFInfo
- Publication number
- CN108868731A CN108868731A CN201810704079.4A CN201810704079A CN108868731A CN 108868731 A CN108868731 A CN 108868731A CN 201810704079 A CN201810704079 A CN 201810704079A CN 108868731 A CN108868731 A CN 108868731A
- Authority
- CN
- China
- Prior art keywords
- time step
- intrinsic fracture
- fracture
- acid
- leakoff coefficent
- 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
- 239000002253 acid Substances 0.000 title claims abstract description 136
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 14
- 238000005260 corrosion Methods 0.000 claims abstract description 14
- 230000007797 corrosion Effects 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 11
- 238000013461 design Methods 0.000 claims abstract description 5
- 239000011435 rock Substances 0.000 claims description 22
- 239000012530 fluid Substances 0.000 claims description 18
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 12
- 241000361919 Metaphire sieboldi Species 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 4
- 238000007789 sealing Methods 0.000 claims description 4
- QRDZSRWEULKVNW-UHFFFAOYSA-N 6-hydroxy-2-oxo-1h-quinoline-4-carboxylic acid Chemical compound C1=C(O)C=C2C(C(=O)O)=CC(=O)NC2=C1 QRDZSRWEULKVNW-UHFFFAOYSA-N 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 3
- 208000005156 Dehydration Diseases 0.000 description 17
- 238000009792 diffusion process Methods 0.000 description 7
- 238000005530 etching Methods 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 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 description 1
- 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
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006557 surface reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Examining Or Testing Airtightness (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种裂缝性储层酸压动态综合滤失系数的计算方法,依次包括以下步骤:(1)计算考虑酸蚀蚓孔的基质滤失系数C mt;(2)计算初始状态第i条天然裂缝滤失系数(C ft 0) i 及初始状态酸压综合滤失系数C t 0;(3)计算Tn时间步第i条天然裂缝内酸液有效作用距离(L ef n) i ;(4)计算Tn时间步第i条天然裂缝宽度(w n) i ;(5)计算Tn时间步第i条天然裂缝滤失系数(C ft n) i ;(6)计算Tn时间步酸压综合滤失系数C t n。本发明考虑了酸液溶蚀天然裂缝壁面导致天然裂缝宽度变化,原理可靠,计算简便,可以处理随机复杂分布的天然裂缝,能够准确快捷地计算出酸压过程中酸液的滤失速度,进而为酸压方案设计提供有效指导。
Description
技术领域
本发明涉及石油天然气勘探开发领域,尤其是酸化压裂领域中一种裂缝性储层酸压动态综合滤失系数的计算方法。
背景技术
裂缝性油气藏储量丰富,是油气勘探开发领域的重点。酸压是此类储层改造的重要手段之一。酸压是利用地面高压泵组,以高于储层吸收能力的排量将酸液泵入地层以产生水力裂缝。同时通过酸液与裂缝壁面岩石发生反应,非均匀刻蚀缝壁岩石,形成沟槽状或凹凸不平的刻蚀裂缝。施工结束后裂缝不完全闭合,最终在地层中形成具有一定几何尺寸和流动能力的人工裂缝,改善油气井的流动状况,从而使油气井获得增产。
酸液的滤失是酸压过程中必然存在的现象,具体是指酸液在流动过程会沿裂缝壁面及天然裂缝流入地层。滤失会导致酸液在水力裂缝中的质量损失,进而影响酸压裂缝的几何尺寸及导流能力,影响改造效果。因此酸压设计过程中需要准确预测酸液的滤失特征。
酸压过程中的滤失特征与水力压裂存在较大的差别。酸液作为反应性流体会溶蚀地层岩石,增加储层的渗透率,从而增加酸液的滤失量。当水力裂缝在延伸过程中与天然裂缝相交时,大量的酸液会通过天然裂缝滤失到地层,同时溶蚀天然裂缝壁面岩石,增加天然裂缝宽度。天然裂缝宽度的增加会显著提高天然裂缝的流动能力,进而增加酸液的滤失量。可以看出,酸压过程中滤失系数是动态变化的,与压裂液滤失过程中的静态滤失系数有较大差异。因此,准确计算酸压过程中的动态滤失系数,尤其是考虑天然裂缝宽度变化的滤失系数非常重要。
关于酸压过程中的动态滤失系数计算方法,国内外学者已经做了许多相关的研究。但大多数都针对于酸液在基质溶蚀形成蚓孔开展,少有考虑天然裂缝在酸压过程中的动态变化。因此这些方法都不能很好地满足裂缝性储层酸压过程中动态滤失系数计算的需求。
发明内容
本发明目的在于提供一种裂缝性储层酸压动态综合滤失系数的计算方法,该方法考虑了酸液溶蚀天然裂缝壁面导致天然裂缝宽度变化,原理可靠,计算简便,可以处理随机复杂分布的天然裂缝,能够准确快捷地计算出酸压过程中酸液的滤失速度,进而为酸压方案设计提供有效指导。
为达到以上技术目的,本发明提供以下技术方案。
一种裂缝性储层酸压动态综合滤失系数的计算方法,依次包括以下步骤:
(1)计算考虑酸蚀蚓孔的基质滤失系数Cmt;
(2)计算初始状态第i条天然裂缝滤失系数(Cft 0)i及初始状态酸压综合滤失系数Ct 0;
(3)计算Tn时间步第i条天然裂缝内酸液有效作用距离(Lef n)i;
(4)计算Tn时间步第i条天然裂缝宽度(wn)i;
(5)计算Tn时间步第i条天然裂缝滤失系数(Cft n)i;
(6)计算Tn时间步酸压综合滤失系数Ct n。
本发明中,为了计算考虑酸液溶蚀天然裂缝壁面导致天然裂缝宽度变化的酸压综合滤失系数,将总注酸时间T化为间隔相同的l个时间节点。相邻时间节点的间隔为:Δt =T / l,时间步定义为Tn,n=0,1,2,3……l,其中T0代表初始状态,即初始时间步。当时间步长较短时,可忽略时间步内缝宽的微小变化对酸液流动的影响,故假设时间步内天然裂缝宽度不随酸液溶蚀而改变。则当天然裂缝表面随机分布有m条天然裂缝时,令单条天然裂缝的编号为i,i=1,2,3……m。
本发明中,所述步骤(1)计算考虑酸蚀蚓孔的基质滤失系数Cmt,具体过程如下:
为了计算酸压过程中于裂缝壁面岩石及形成的酸蚀蚓孔造成的滤失,采用公式(1)计算考虑酸蚀蚓孔的基质滤失系数Cmt(Hill,A.D.,Zhu,D.,and Wang,Y.M.,1995.TheEffect of Wormholing on the Fluid-Loss Coefficient in Acid Fracturing.SPEProduction&Facilities,November,257-263)。此滤失系数能够表征酸压过程中考虑酸蚀蚓孔发展的滤失量变化,故只需要在步骤(1)中计算一次:
式中:Cmt—考虑酸蚀蚓孔的基质滤失系数,
Cc—非蚓孔控制的滤失系数,
Cwh—蚓孔控制的滤失系数,
φ—岩石孔隙度,无因次;
ct—地层流体压缩系数,Pa-1;
k—岩石渗透率,m2;
Pf—酸压裂缝内流体压力,Pa;
Pr—地层流体压力,一般取原始地层压力,Pa;
μ—地层流体粘度,Pa·s;
Qwh—酸液突破体积,一般通过酸液驱替实验获得,无因次。
本发明中,所述步骤(2)计算初始状态第i条天然裂缝滤失系数(Cft 0)i及初始状态酸压综合滤失系数Ct 0,具体过程如下:
当n=0,即为初始状态时,裂缝内酸岩反应尚未发生,则可直接采用式(4)、式(5)计算获得第i条天然裂缝初始的天然裂缝滤失系数(Cft 0)i及初始入口面积(Af 0)后,带入式(6)计算获得初始状态酸压综合滤失系数Ct 0(Liu Y,Guo J,Chen Z.Leakoffcharacteristics and an equivalent leakoff coefficient in fractured tight gasreservoirs[J].Journal of Natural Gas Science&Engineering,2016,31:603-611):
式中:(w0)i—第i条天然裂缝初始宽度,由储层地质资料获得,m;
—第i条天然裂缝初始天然裂缝滤失系数,
(H)i—第i条天然裂缝高度,由储层地质资料获得,m;
—第i条天然裂缝初始入口面积,m2;
Am—单翼单面裂缝面积,由酸压设计资料中获得,m2。
本发明中,所述步骤(3)计算Tn时间步第i条天然裂缝内酸液有效作用距离(Lef n)i,具体过程如下:
本发明中,为了求取随时间步数n变化的酸压综合滤失系数,需在初始状态参数的基础上,按时间步依次进行循环计算。例如针对Tn时间步的计算,所需要输入的参数则均为Tn-1时间步计算所得的结果或能根据Tn-1时间步结果计算。同时本发明中计算过程应从T1时间步开始,基于初始状态参数开展循环计算。
Tn时间步第i条天然裂缝内酸液有效作用距离(Lef n)i可由式(7)计算。此处(Lef n)i实际为Tn-1到Tn时间步所经历的时间中酸液在天然裂缝内的酸液有效作用距离,根据本发明中忽略时间步内缝宽微小变化对酸液流动影响的假设,(Lef n)i的值由Tn-1时间步获得的已知参数计算:
式中:—第i条天然裂缝Tn时间步酸液有效作用距离,m;
(wn-1)i—第i条天然裂缝Tn-1时间步宽度,m;
—第i条天然裂缝Tn-1时间步内平均裂缝滤失速度,m/s;
—第i条天然裂缝Tn-1时间步酸液传质系数,m/s;
—第i条天然裂缝Tn-1时间步天然裂缝滤失系数,
kc—反应速率常数,通过酸岩反应动力学实验获得,m/s;
Δt—间时步步长,s。
由于天然裂缝内流速及缝宽均较小,故流态一般为层流。式(7)中Tn-1时间步第i条天然裂缝内酸液传质系数(kg n-1)i可采用式(9)计算获得(C.Y.Chang,J.A.Guin,andL.D.Roberts.Surface Reaction With Combined Forced and Free Convection[J].AIChEJ,1976,22,No.2:252-259):
式(7)推导过程如下:
假设Tn-1到Tn时间步内天然裂缝内酸液流动反应为稳态二维均匀流动扩散反应,根据平衡原理可写出二维对流扩散方程:
式中:—酸浓度在缝长方向的梯度,mol/(m3·m);
—酸浓度在缝宽方向的梯度,mol/(m3·m);
—H+扩散系数,m2/s。
根据假设,对流项(式10左边)通量面积为天然裂缝入口面积,扩散项(式10右边)通量面积为酸液有效作用距离乘以天然裂缝高度。故对式(10)在反应体积内进行积分,并通过高斯定理展开有:
式中:Ce—反应控制体积出口边界酸浓度,mol/m3;
Cw—反应控制体积入口边界酸浓度,mol/m3;
—第i条天然裂缝Tn-1时间步入口面积,m2。
式(11)中为扩散边界项,根据酸压裂缝传质-反应耦合边界(式12,式13)对该项进行处理(Roberts,L.D.,&GUIN,J.A.(1975,August 1).A New Method forPredicting Acid Penetration Distance.Society of Petroleum Engineers.doi:10.2118/5155-PA):
式中:CP—反应控制体积中心酸浓度,mol/m3;
Cwall—天然裂缝壁面酸浓度,mol/m3。
将式(13)带入式(12),消去边界上未知酸液浓度Cwall,则可得式(11)中扩散项边界条件可表示为:
另假设进入天然裂缝的酸液在反应控制体积内完全反应,则反应出口边界e处酸液浓度Ce=0。入口处H+通量等于单位时间内进入天然裂缝的H+量,则式(11)左边第二项可表示为:
式中:(Bn-1)i—Tn-1时间步进入第i条天然裂缝的酸液流量分量,无因次;
—Tn-1时间步酸液平均综合滤失速度,m/s;
Cacid—天然裂缝缝口H+浓度,mol/m3。
则根据线性插值原理,式(11)中反应控制体积中心酸浓度CP可以根据线性插值理论根据式(16)进行计算:
将式(14),式(15)、式(16)带入式(11)得式(17):
将式(17)方程两边进行消参,则可Tn时间步第i条天然裂缝内酸液有效作用距离计算公式:
本发明中,所述步骤(4)计算Tn时间步第i条天然裂缝宽度(wn)i,具体过程如下:
Tn时间步第i条天然裂缝宽度(wn)i根据式(18)计算
式中:M—H+摩尔质量,kg/mol;
β—酸液对岩石的溶解能力,kg/kg;
ρ—岩石密度,kg/m3;
n—时间步数,无因次;
—Tn-1时间步酸液平均综合滤失速度,m/s;
(wn)i—第i条天然裂缝Tn时间步宽度,m;
(wn-1)i—第i条天然裂缝Tn-1时间步宽度,m;
(Δwn-1)i—Tn-1时间步第i条天然裂缝酸液溶蚀产生的宽度变化,m;
Cacid—天然裂缝缝口H+浓度,mol/m3;
(Bn-1)i—Tn-1时间步进入第i条天然裂缝的酸液流量分量,无因次;
—Tn-1时间步酸压综合滤失系数,由上个时间步的计算结果获取,
Tn-1时间步进入第i条天然裂缝的酸液流量分量(Bn-1)i,可表达为进入第i条天然裂缝的酸液流量与滤失入基质的酸液流量以及进入所有天然裂缝的酸液流量加和之比。
式中:—Tn-1时间步平均基质滤失速度,m/s;
—第i条天然裂缝Tn-1时间步入口面积,m2;
—第i条天然裂缝Tn-1时间步内平均裂缝滤失速度,m/s。
式(18)中Tn-1时间步第i条天然裂缝酸液溶蚀产生的宽度变化推导过程如下:
步骤(3)中假设进入裂缝的酸液在酸蚀有效作用距离之内完全反应。故Tn-1时间步第i条天然裂缝内参加反应的H+摩尔量为则反应掉的岩石体积可写为同时,参加反应的天然裂缝壁面面积可写为故Tn-1时间步第i条天然裂缝酸液溶蚀产生的宽度变化可写为:
本发明中,所述步骤(5)计算Tn时间步第i条天然裂缝滤失系数(Cft n)i,具体过程如下:
Tn时间步第i条天然裂缝的表观滤失系数(Cf n)i可根据式(24)进行计算:
式中:—Tn时间步第i条天然裂缝的表观滤失系数,
天然裂缝在酸压过程中会因为酸液在缝口处的大量溶蚀,而变成“喇叭状”。因此,入口处的流动能力会远远大于天然裂缝内部未被大量溶蚀的部分。而由于流动的连续性,天然裂缝的整体流动能力应由天然裂缝的平均宽度决定,直接采用缝口宽度计算会导致较大误差。因此对长度为Ln的天然裂缝根据式(25)在天然裂缝长度方向进行Tn-1时间步的天然裂缝滤失系数与Tn时间步的表观滤失系数的加权平均以获得Tn时间步第i条天然裂缝滤失系数(Cft n)i:
式中:(Ln)i—第i条天然裂缝长度,由油田地质资料获得,m;
当大于(Ln)i时,则取
本发明中,所述步骤(6)计算Tn时间步酸压综合滤失系数Ct n,具体过程如下:
根据式(26)可计算Tn时间步酸压综合滤失系数:
附图说明
图1是酸压综合动态滤失系数Ct n计算结果展示图。
具体实施方式
下面结合附图和实施例进一步说明本发明。需要指出的是,该实施方式旨在对本发明做示意性说明和解释,并不限定本发明的范围。
实施例
假设裂缝壁面存在两条天然裂缝,输入参数如表1所示。
表1本井的储层物性资料、流体参数资料、井身结构资料
参数 | 数值 | 参数 | 数值 |
酸压裂缝内流体压力,106Pa | 70 | 第一条天然裂缝初始宽度,m | 0.00012 |
地层流体压力,106Pa | 60 | 第一条天然裂缝高度,m | 13 |
岩石孔隙度,% | 8 | 第二条天然裂缝初始宽度,m | 0.00018 |
岩石渗透率,10-15m2 | 1 | 第二条天然裂缝高度,m | 19 |
地层流体压缩系数,106Pa-1 | 5×10-5 | 酸岩反应常数,m/s | 1.5×10-5 |
地层流体粘度,Pa·s | 0.001 | 酸液对岩石的溶解能力,kg/kg | 1.37 |
单翼单面裂缝面积,m2 | 3000 | H+扩散系数,m2/s | 5×10-5 |
H+摩尔质量,kg/mol; | 36.5 | 酸压裂缝内酸液浓度,mol/m3 | 6020 |
岩石密度,kg/m3; | 2600 | 酸压时间,s | 3600 |
时间步数,无因次 | 60 |
A基于表1中的数据,根据步骤(1)中式(1)~(3)计算,本实例中考虑酸蚀蚓孔的基质滤失系数
B基于表1中的数据,根据步骤(2)中式(4)~(5)计算出两条天然裂缝T0时间步的天然裂缝滤失系数。计算结果为
C基于表1中的数据及B中的计算结果,根据步骤(2)中式(6)计算T0时间步酸压综合滤失系数
D基于表1中的数据及B、C中的计算结果,根据步骤(3)中式(7)~(9)计算两条天然裂缝内T1时间步酸液有效作用距离。计算结果为:
E基于表1中的数据及D中的计算结果,根据步骤(4)中式(18)~(22)计算得两条天然裂缝T1时间步的裂缝宽度。计算结果为:(w1)1=8.1×10-6m以及(w1)2=11.8×10-6m;
F基于表1中的数据及E中的计算结果,根据步骤(5)中式(24)~(25)计算两条天然裂缝的天然裂缝滤失系数。此时,由于及均大于Ln,则式(25)计算式取计算得T1时间步两条天然裂缝的天然裂缝滤失系数分别为:
G基于表1中的数据及E、F中的计算结果,根据步骤(6)中式(26)计算T1时间步酸压综合滤失系数Ct 1,计算结果为
H将E、F、G中计算所得结果带入步骤(3),重复步骤(3)~(6),即可获得随时间变化的酸压综合滤失系数Ct n,计算结果及与不考虑酸液溶蚀天然裂缝的综合滤失系数的对比展示于图1。
Claims (9)
1.一种裂缝性储层酸压动态综合滤失系数的计算方法,依次包括以下步骤:
(1)计算考虑酸蚀蚓孔的基质滤失系数Cmt;
(2)计算初始状态第i条天然裂缝滤失系数(Cft 0)i及初始状态酸压综合滤失系数Ct 0;
(3)计算Tn时间步第i条天然裂缝内酸液有效作用距离(Lef n)i;
(4)计算Tn时间步第i条天然裂缝宽度(wn)i;
(5)计算Tn时间步第i条天然裂缝滤失系数(Cft n)i;
(6)计算Tn时间步酸压综合滤失系数Ct n。
2.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(1)根据下式计算考虑酸蚀蚓孔的基质滤失系数Cmt:
式中:Cc—非蚓孔控制的滤失系数,
Cwh—蚓孔控制的滤失系数,
φ—岩石孔隙度,无因次;
ct—地层流体压缩系数,Pa-1;
k—岩石渗透率,m2;
Pf—酸压裂缝内流体压力,Pa;
Pr—地层流体压力,一般取原始地层压力,Pa;
μ—地层流体粘度,Pa·s;
Qwh—酸液突破体积,一般通过酸液驱替实验获得,无因次。
3.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(2)计算初始状态第i条天然裂缝滤失系数(Cft 0)i及初始状态酸压综合滤失系数Ct 0,过程如下:
式中:(w0)i—第i条天然裂缝初始宽度,由储层地质资料获得,m;
—第i条天然裂缝初始天然裂缝滤失系数,
(H)i—第i条天然裂缝高度,由储层地质资料获得,m;
—第i条天然裂缝初始入口面积,m2;
Am—单翼单面裂缝面积,由酸压设计资料中获得,m2。
4.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(3)通过下式计算Tn时间步第i条天然裂缝内酸液有效作用距离(Lef n)i:
式中:(wn-1)i—第i条天然裂缝Tn-1时间步宽度,m;
—第i条天然裂缝Tn-1时间步内平均裂缝滤失速度,m/s;
—第i条天然裂缝Tn-1时间步酸液传质系数,m/s;
—第i条天然裂缝Tn-1时间步天然裂缝滤失系数,
kc—反应速率常数,通过酸岩反应动力学实验获得,m/s;
Δt—间时步步长,s。
5.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(4)根据下式计算Tn时间步第i条天然裂缝宽度(wn)i:
式中:M—H+摩尔质量,kg/mol;
β—酸液对岩石的溶解能力,kg/kg;
ρ—岩石密度,kg/m3;
n—时间步数,无因次;
—Tn-1时间步酸液平均综合滤失速度,m/s;
(wn)i—第i条天然裂缝Tn时间步宽度,m;
(wn-1)i—第i条天然裂缝Tn-1时间步宽度,m;
(Δwn-1)i—Tn-1时间步第i条天然裂缝酸液溶蚀产生的宽度变化,m;
Cacid—天然裂缝缝口H+浓度,mol/m3;
(Bn-1)i—Tn-1时间步进入第i条天然裂缝的酸液流量分量,无因次;
—Tn-1时间步酸压综合滤失系数,由上个时间步的计算结果获取,
6.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(5)计算Tn时间步第i条天然裂缝滤失系数(Cft n)i,过程如下:
式中:(Ln)i—第i条天然裂缝长度,由油田地质资料获得,m;
当大于(Ln)i时,则取
—Tn时间步第i条天然裂缝的表观滤失系数,
7.如权利要求1所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述步骤(6)通过下式计算Tn时间步酸压综合滤失系数Ct n:
8.如权利要求4所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述Tn-1时间步第i条天然裂缝内酸液传质系数(kg n-1)i通过下式计算获得:
9.如权利要求5所述的一种裂缝性储层酸压动态综合滤失系数的计算方法,其特征在于,所述Tn-1时间步进入第i条天然裂缝的酸液流量分量(Bn-1)i通过下式计算获得:
式中:—Tn-1时间步平均基质滤失速度,m/s;
—第i条天然裂缝Tn-1时间步入口面积,m2;
—第i条天然裂缝Tn-1时间步内平均裂缝滤失速度,m/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810704079.4A CN108868731B (zh) | 2018-06-30 | 2018-06-30 | 一种裂缝性储层酸压动态综合滤失系数的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810704079.4A CN108868731B (zh) | 2018-06-30 | 2018-06-30 | 一种裂缝性储层酸压动态综合滤失系数的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108868731A true CN108868731A (zh) | 2018-11-23 |
CN108868731B CN108868731B (zh) | 2020-05-01 |
Family
ID=64297799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810704079.4A Active CN108868731B (zh) | 2018-06-30 | 2018-06-30 | 一种裂缝性储层酸压动态综合滤失系数的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108868731B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
CN111197476A (zh) * | 2020-01-08 | 2020-05-26 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
CN111219175A (zh) * | 2020-01-08 | 2020-06-02 | 西南石油大学 | 考虑应力敏感的裂缝性碳酸盐岩酸压裂缝匹配性优化方法 |
CN112989716A (zh) * | 2021-02-05 | 2021-06-18 | 成都北方石油勘探开发技术有限公司 | 一种岩板人工裂缝酸刻蚀测试参数计算方法 |
CN113969773A (zh) * | 2021-10-21 | 2022-01-25 | 中国石油化工股份有限公司 | 一种粗糙天然裂缝压裂液滤失测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5005643A (en) * | 1990-05-11 | 1991-04-09 | Halliburton Company | Method of determining fracture parameters for heterogenous formations |
CN106777663A (zh) * | 2016-12-12 | 2017-05-31 | 西南石油大学 | 一种考虑天然裂缝的压裂液滤失速度计算方法 |
CN107992690A (zh) * | 2017-12-06 | 2018-05-04 | 西南石油大学 | 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法 |
-
2018
- 2018-06-30 CN CN201810704079.4A patent/CN108868731B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5005643A (en) * | 1990-05-11 | 1991-04-09 | Halliburton Company | Method of determining fracture parameters for heterogenous formations |
CN106777663A (zh) * | 2016-12-12 | 2017-05-31 | 西南石油大学 | 一种考虑天然裂缝的压裂液滤失速度计算方法 |
CN107992690A (zh) * | 2017-12-06 | 2018-05-04 | 西南石油大学 | 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法 |
Non-Patent Citations (3)
Title |
---|
A.D.HILL 等: "The Effect of Wormholing on the Fluid-Loss Coefficient in Acid Fracturing", 《SPE PRODUCTION & FACILITIES》 * |
B.BAZIN 等: "Acid Filtration in Dynamic Conditions to Mimic Fluid Loss in Acid Fracturing", 《SPE》 * |
YE ZHANG 等: "Wormhole Propagation Behavior and Its Effect on Acid Leakoff under In Situ Conditions in Acid Fracturing", 《TRANS POROUS MED》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
CN110609974B (zh) * | 2019-09-09 | 2022-09-16 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
CN111197476A (zh) * | 2020-01-08 | 2020-05-26 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
CN111219175A (zh) * | 2020-01-08 | 2020-06-02 | 西南石油大学 | 考虑应力敏感的裂缝性碳酸盐岩酸压裂缝匹配性优化方法 |
CN111197476B (zh) * | 2020-01-08 | 2020-10-09 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
CN112989716A (zh) * | 2021-02-05 | 2021-06-18 | 成都北方石油勘探开发技术有限公司 | 一种岩板人工裂缝酸刻蚀测试参数计算方法 |
CN112989716B (zh) * | 2021-02-05 | 2023-03-24 | 成都北方石油勘探开发技术有限公司 | 一种岩板人工裂缝酸刻蚀测试参数计算方法 |
CN113969773A (zh) * | 2021-10-21 | 2022-01-25 | 中国石油化工股份有限公司 | 一种粗糙天然裂缝压裂液滤失测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108868731B (zh) | 2020-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113826099B (zh) | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 | |
CN108868731A (zh) | 一种裂缝性储层酸压动态综合滤失系数的计算方法 | |
Hu | An improved analytical model for vertical borehole ground heat exchanger with multiple-layer substrates and groundwater flow | |
Li et al. | Experimental and numerical studies on heat transfer characteristics of vertical deep-buried U-bend pipe to supply heat in buildings with geothermal energy | |
Liu et al. | Modeling and simulation of wormhole formation during acidization of fractured carbonate rocks | |
CN105840187B (zh) | 致密性油藏水平井分段压裂产能计算方法 | |
Xu et al. | Numerical modeling of gas extraction from coal seam combined with a dual-porosity model: Finite difference solution and multi-factor analysis | |
Perera et al. | Numerical simulation of gas flow through porous sandstone and its experimental validation | |
CN106295095B (zh) | 基于常规测井资料预测低渗透砂岩储层产能的方法 | |
Liu et al. | Fractal analysis of directional permeability of gas shale fracture networks: A numerical study | |
Xu et al. | Study on relative permeability characteristics affected by displacement pressure gradient: Experimental study and numerical simulation | |
Yu et al. | The anisotropic seepage analysis of water-sealed underground oil storage caverns | |
CN109611067B (zh) | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 | |
CN106951649B (zh) | 一种测定水平井sagd蒸汽腔扩展速度的方法 | |
CN103590824A (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
CN104297126A (zh) | 低渗透储层气体渗流启动压力梯度测量装置及测量方法 | |
Zhang et al. | Efficient flow rate profiling for multiphase flow in horizontal wells using downhole temperature measurement | |
CN103726815B (zh) | 一种co2驱采出井井筒流态确定及参数优化方法 | |
CN112304842B (zh) | 一种页岩油co2/n2交替驱替注入量模拟分析方法 | |
Guo et al. | Analysis of 2D flow and heat transfer modeling in fracture of porous media | |
CN111734394B (zh) | 一种确定致密油藏压裂井不定常流井底压力的方法 | |
CN114201932A (zh) | 一种复杂情况下的致密油藏压裂井试井模拟方法 | |
Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
Xu et al. | A fractal physics-based data-driven model for water-flooding reservoir (FlowNet-fractal) | |
Wang et al. | Comparative study on effects of macroscopic and microscopic fracture structures on the performance of enhanced geothermal systems |
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 |