CN109611067A - 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 - Google Patents
深层石灰岩储层酸压酸液有效作用距离的数值计算方法 Download PDFInfo
- Publication number
- CN109611067A CN109611067A CN201910009184.0A CN201910009184A CN109611067A CN 109611067 A CN109611067 A CN 109611067A CN 201910009184 A CN201910009184 A CN 201910009184A CN 109611067 A CN109611067 A CN 109611067A
- Authority
- CN
- China
- Prior art keywords
- acid
- unit
- temperature
- rock
- reaction heat
- 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 178
- 238000000034 method Methods 0.000 title claims abstract description 38
- 235000019738 Limestone Nutrition 0.000 title claims abstract description 19
- 239000006028 limestone Substances 0.000 title claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 64
- 239000011435 rock Substances 0.000 claims abstract description 63
- 238000010276 construction Methods 0.000 claims abstract description 26
- 238000002513 implantation Methods 0.000 claims abstract description 6
- 238000012546 transfer Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 239000007788 liquid Substances 0.000 claims description 5
- 239000007789 gas Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- GPRLSGONYQIRFK-UHFFFAOYSA-N hydron Chemical compound [H+] GPRLSGONYQIRFK-UHFFFAOYSA-N 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 9
- 239000000446 fuel Substances 0.000 abstract description 4
- 239000003208 petroleum Substances 0.000 description 9
- 238000009795 derivation Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005530 etching Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000007787 solid Substances 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
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000008733 Citrus aurantifolia Nutrition 0.000 description 1
- 235000011941 Tilia x europaea Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000004571 lime Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction 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/16—Enhanced recovery methods for obtaining hydrocarbons
-
- 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/06—Measuring temperature or pressure
-
- 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/06—Measuring temperature or pressure
- E21B47/07—Temperature
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (AREA)
- Geophysics (AREA)
- Primary Health Care (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明公开了深层石灰岩储层酸压酸液有效作用距离的数值计算方法，包括：（1）计算施工结束时的井底温度；（2）计算施工结束时的裂缝缝长，在缝长方向划分网格单元并计算各单元的平均缝宽，得到各单元酸液流速；（3）以井底温度为初始条件，计算各单元不考虑酸岩反应热的酸液温度以及酸液浓度；（4）计算各单元的酸岩摩尔反应热；（5）计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，建立酸液浓度分布曲线，定义酸液初始注入浓度的10%作为残酸极限浓度，以残酸极限浓度对应的缝长为酸液有效作用距离。本发明考虑了酸压过程中酸液受到的热效应，提高了酸液有效作用距离计算的准确性，对深层石灰岩储层酸压设计具有指导意义。
Description
技术领域
本发明涉及石油工程领域，尤其是深层石灰岩储层酸压过程中酸液有效作用距离的数值计算方法。采用数值计算方法计算高温条件下的酸液有效作用距离，提高深层石灰岩储层酸压设计的准确性，实现深层石灰岩储层高效酸压改造。
背景技术
川西地区二叠统深层碳酸盐岩储层展现出良好的勘探潜力，已成为四川盆地现阶段最现实的接替层系之一。储层以石灰岩、白云岩为主，具有埋藏深(>5000m)、温度高(>145℃)等特征。该地区开发难以自然获得工业油气流，酸压裂技术是区块储层主要的增产改造手段之一(齐宝权,赵佐安,贺洪举等.应用测井储层品质综合评价指数来判别储层的有效性——以四川盆地西部地区二叠系栖霞组为例[J].天然气工业,2018,38(02):25-32)。
酸化压裂(简称酸压)是指在高于储层破裂压力或天然裂缝的闭合压力下，将酸液挤入储层，在储层中形成裂缝，同时酸液与裂缝壁面岩石发生化学反应，非均匀刻蚀裂缝壁岩石，形成沟槽状或凹凸不平的刻蚀裂缝，施工结束后裂缝不完全闭合，最终形成具有一定几何尺寸和导流能力的酸蚀裂缝，实现油气井增产。酸液有效作用距离是评价酸压效果的重要指标之一，它主要受温度、酸液浓度等因素影响。其中，温度是影响石灰岩储层酸压效果的重要因素(李颖川.采油工程[M].石油工业出版社,2009)。
目前常用的酸液有效作用距离计算方法主要关注残酸极限浓度的计算(李建辉,慕立俊,赵振峰,等.一种基于残酸极限评价酸液有效作用距离的方法,CN105044288A[P].2015)，很少综合考虑压裂过程中酸液的温度变化对于酸压改造效果的影响。因此，针对川西地区储层温度高的特点，研究酸压过程酸液温度变化规律，明确温度对于酸压效果(酸液有效作用距离)的影响，对于川西地区储层酸压改造设计至关重要。
发明内容
本发明的目的在于提供深层石灰岩储层酸压酸液有效作用距离的数值计算方法，该方法重点研究酸压过程中酸液受到的热效应(井筒传热、酸液流动换热、酸岩反应放热)，采用数值方法计算考虑以上热效应的酸液温度与浓度，最终基于酸液浓度确定酸液有效作用距离。该方法原理可靠，有利于计算高温地层的酸液温度及浓度，提高酸液有效作用距离计算的准确性，对深层石灰岩储层酸压设计具有指导意义。
为达到以上技术目的，本发明提供以下技术方案。
首先，建立井筒温度场模型模拟井筒传热并计算井底温度；其次，采用PKN公式计算酸压裂缝的缝长和缝宽，根据裂缝缝宽分布求解酸液流速；再者，以井底温度为初始条件，根据流速求解酸液温度和浓度；然后，通过酸液温度和浓度计算酸岩反应热；最后，计算考虑酸岩反应热后的酸液温度和浓度，建立酸液浓度与裂缝长度的关系曲线，确定酸液有效作用距离。
深层石灰岩储层酸压酸液有效作用距离的数值计算方法，依次包括以下步骤：
(1)采用Eickmeier井筒温度场模型模拟井筒温度场并计算施工结束时的井底温度(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998)；
(2)采用PKN公式计算施工结束时的裂缝缝长(李颖川.采油工程[M].北京:石油工业出版社,2001)，在缝长方向划分网格单元并计算各单元的平均缝宽，再计算各单元酸液流速；
(3)以步骤(1)计算的井底温度为初始条件，通过步骤(2)计算的酸液流速计算各单元不考虑酸岩反应热的酸液温度以及酸液浓度；
(4)通过步骤(3)中计算的酸液温度和浓度计算各单元的酸岩摩尔反应热；
(5)将步骤(4)计算的酸岩摩尔反应热代入步骤(3)中，计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，建立酸液浓度分布曲线。定义酸液初始注入浓度的10％作为残酸极限浓度，根据分布曲线以残酸极限浓度对应的缝长为酸液有效作用距离。
在本发明中，所述步骤(1)中，采用Eickmeier井筒温度场模型模拟井筒温度场并计算施工结束时的井底温度，过程如下：
采用Eickmeier井筒温度场模型建模方法建立目标井井筒温度场模型并划分网格单元。以油管中心为轴：径向上，定义ri为网格单元半径，划分N个网格单元(i＝0,1,2,…,N)。其中模型边界点rN点的选择应满足该点处温度始终等于地层温度；纵向上，设目的层深度为H，单元深度为ΔH，则纵向上可划分为M个单元(j＝0,1,2…,M)，边界点M为(H/ΔH)取整。定义施工总时间为tT，时间步长为Δt，时间步为n(n＝0,1,2…,tT)。结合目标井地质工程参数，利用下式计算n时刻单元(i,j)的温度最终得到施工结束时的井底温度
式中：Q——注酸排量，m3/min；
ρi——网格单元i的密度，kg/m3；
Ci——网格单元i的比热容，J/(kg·℃)；
ΔHj——网格单元j的高度，m；
ri——模型径向网格半径，m；
λi——网格单元i的热传导系数，W/(m·℃)；
Δt——时间步长，min；
tT——施工总时间，min；
n——时间步，min；
——n时刻网格单元(i,j)的温度(为施工结束时的井底温度)，℃。
在本发明中，所述步骤(2)中，采用PKN公式计算施工结束时的裂缝缝长，在缝长方向划分网格单元并计算各单元的平均缝宽，再计算各单元酸液流速，过程如下：
1)通过计算施工结束即n＝tT时的裂缝缝长，定义x为缝长方向并划分网格单元，计算各单元对应的平均缝宽w：
式中：L——施工结束时的缝长，m；
c——酸液滤失系数，m/min0.5；
μf——酸液粘度，Pa·min；
h——裂缝高度，m；
E——储层岩石杨氏模量，Pa
ε——储层岩石泊松比，无因次；
x——缝长方向上单元坐标，m；
w——施工结束时各单元的平均缝宽，m。
2)通过下式得到各单元酸液流速vx：
式中：vl——各单元滤失速度，m/s；
vx——各单元x方向流速，m/s。
公式(3)的推导过程如下：
酸液连续性方程为：
假设裂缝内酸液流速考虑为平均流速，缝内压力均匀分布，在裂缝轴线附近(即y＝0处)的压力梯度为0，所以有vy|y＝0＝0。另外，裂缝壁面上酸液流速等于滤失速率，即因此将式(4)由裂缝轴线(y＝0)至裂缝壁面(y＝w/2)积分，即可推导公式(3)。
在本发明中，所述步骤(3)中，以步骤(1)计算的井底温度为初始条件，通过步骤(2)计算的酸液流速计算各单元不考虑酸岩反应热的酸液温度以及酸液浓度，过程如下：
1)假设裂缝入口温度等于步骤(1)计算的井底温度，缝长方向各单元温度考虑为平均温度，将步骤(2)求解的酸液流速代入下式求解各单元酸液温度与裂缝壁面温度：
其中，对流换热系数可由努谢尔特准则求解：
hT＝NNuKhf/w (7)
式中：T——各单元酸液平均温度，℃；
Tw——各单元裂缝壁面温度，℃；
hT——对流换热系数，J/(m2·min·℃)；
NNu——努谢尔特数，无因次；
Khf——酸液的导热系数，J/(m·min·℃)；
ρf——储层岩石密度，kg/m3；
Cf——储层岩石比热容，J/(kg·℃)。
2)假设酸液对流换热热量hT(Tw-T)主要由地层传热热量以及酸岩反应热热量提供：
hT(Tw-T)＝q(tT)-kRCmΔrQm(Tw,p) (8)
式中：kR——酸岩反应速率常数，(kg/m3)-m·mol/(m2·s)；
C——各单元酸液浓度，kg/m3；
ΔrQm(Tw,p)——酸岩摩尔反应热，kJ/mol；
m——反应级数，无因次；
q(tT)——施工结束地层传热热量，kJ。
其中，热流函数公式为(Whitsitt,Dysart.The effect of temperature onstimulation design[J].SPE 2497,1970)：
式中：Mma——地层岩石体积热容，J/m3·℃；
Khr——地层岩石导热系数，W/(m·℃)；
——孔隙度，％。
由于酸岩反应热是关于酸液温度和浓度的函数，因此在步骤(3)中先不考虑酸岩反应热(设其值为0)，将式(9)、式(8)代入式(6)即可求出各单元不考虑酸岩反应热的酸液温度。
式(6)的推导过程如下：
假设酸液温度场为稳态。各单元可建立能量守恒方程：
同公式(5)的处理方法，对公式(10)由裂缝轴线(y＝0)至裂缝壁面(y＝w/2)积分：
由于裂缝轴线(y＝0)速度梯度为0，因此温度梯度也为0。各单元温度考虑为平均温度，将式(11)处理为常微分方程：
对流换热系数为表征流体与固体表面之间换热能力的参数，定义为单位面积交换热量与流体固体温度差的比值(何潮洪，冯霄.化工原理[M].北京：科学出版社，2001:190)。因此，定义hT为对流换热系数：
需要说明的是，式(7)是基于努谢尔数定义推导的求解公式，是hT另一种表现形式。将式(3)、式(13)代入式(12)，即可推导式(6)。
3)假设各单元酸液浓度为平均浓度，将步骤(2)求解的酸液流速代入下式得到各单元酸液浓度与裂缝壁面浓度：
式中：kg——对流传质系数，m/min；
Cw——各单元裂缝壁面酸液浓度，kg/m3；
假设裂缝壁面的酸液与裂缝壁面岩石完全反应，因此有：
kg(C-Cw)＝kRCw m (15)
式中，由舍伍德数定义，对流传质系数可以表示为：
kg＝DeNSh/w (16)
酸岩反应速率常数KR由阿伦尼乌斯公式求解(李颖川.采油工程[M].北京:石油工业出版社,2001)：
式中：De——酸液氢离子有效传质系数，m2/s；
Nsh——舍伍德数，无因次；
k0——频率因子，(kg/m3)-m·mol/(m2·s)；
Ea——反应活化能，J/mol；
R——气体常数，8.314J/(mol·℃)。
因此，将步骤2)求解的各单元酸液温度代入式(17)求解得到KR，而后联立求解式(14)、式(15)即可求得各单元酸液浓度。
式(14)的推导过程如下：
各单元建立酸液浓度物质平衡方程：
同公式(5)的处理方法，对公式(18)由裂缝轴线(y＝0)至裂缝壁面(y＝w/2)积分：
由于裂缝轴线(y＝0)速度梯度为0，因此浓度梯度也为0。各单元浓度考虑为平均浓度，将式(19)处理为常微分方程：
由对流传质系数定义(刘瑾,王豪,王治红.传热与传质[M].西南石油大学,2007)：
需要说明的是，式(16)是基于舍伍德数定义推导的求解公式，是kg另一种表现形式。将式(3)、式(21)代入式(20)，即可推导式(14)。综上，由步骤(3)可求得各单元不考虑酸岩反应热的酸液温度和浓度。
在本发明中，所述步骤(4)中，结合步骤(3)中计算的酸液温度和浓度计算各单元的酸岩摩尔反应热，过程如下：
1)考虑目标工区储层岩性为石灰岩，通过查阅物理化学手册(孙艳辉,何广平,马国正等.物理化学实用手册[M].北京:化学工业出版社.2016)获得反应体系中各组分的标准摩尔反应热，并计算反应体系的摩尔反应热
式中：υi——反应体系中i组分的化学计量系数，无因次；
——i组分的标准摩尔反应热，kJ/mol；
——反应体系标准摩尔反应热，kJ/mol。
2)结合步骤1)的计算结果，计算高温高压条件酸岩摩尔反应热ΔrHm(Tw,p)(GuoJ,Liu H,Zhu Y,et al.Effects of acid–rock reaction heat on fluid temperatureprofile in fracture during acid fracturing in carbonate reservoirs[J].Journalof petroleum science&engineering,2014,122:31-37)：
其中，fg为游离态CO2摩尔分数：
式中：ΔrHm(Tw)——温度Tw下的摩尔反应热，kJ/mol；
p——地层压力，MPa；
Cp,m(i,T)——组分i在温度T下的常压摩尔热容，J/(℃·mol)；
fg——CO2摩尔分数，无因次；
Vco2——CO2摩尔体积，m3/mol；
Vacid——残酸体积，m3；
Sco2——残酸中CO2溶解度，m3/m3。
ΔrHm(Tw,p)——温度Tw压力p下的摩尔反应热，kJ/mol。
3)酸岩反应过程中，反应体系对CO2做体积功。因此，综合考虑温度、压力和CO2的体积功的影响，得到地层条件下酸岩摩尔反应热ΔrQm(Tw,p)：
在本发明中，所述步骤(5)中，将步骤(4)计算的酸岩摩尔反应热代入步骤(3)中，计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，建立酸液浓度分布曲线。定义酸液初始注入浓度的10％作为残酸极限浓度，根据分布曲线以残酸极限浓度对应的缝长为酸液有效作用距离，过程如下：
1)将步骤(4)计算的酸岩摩尔反应热代入公式(8)中计算各单元酸岩反应热，而后由步骤(3)计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，并作缝长方向各单元酸液浓度分布曲线；
2)定义酸液初始注入浓度的10％作为残酸极限浓度，参考酸液浓度分布曲线，取残酸极限浓度对应的裂缝缝长为酸液有效作用距离。
与现有技术相比，本发明提供了一种针对深层石灰岩储层的酸液有效作用距离的数值计算方法。本方法在计算过程中耦合了井筒传热、酸液对流换热、酸岩反应热三种热效应：①首先通过井筒温度场模型计算施工结束井底温度；②利用PKN公式计算的缝宽求解缝长方向各单元酸液流速，以井底温度为初始条件，通过酸液流速求解各单元酸液温度以及浓度，并实现井筒传热与对流换热的耦合；③通过各单元酸液温度、浓度求解酸岩反应热，耦合酸岩反应热后计算各单元酸液浓度得到酸液浓度分布曲线，确定酸液有效作用距离。本方法克服了现有方法未能综合考虑压裂过程中多种热效应对于酸压改造效果影响的局限。
附图说明
图1是本发明中井筒温度场模型网格划分示意图。
图2是本发明中井底温度变化模拟结果图。
图3是本发明中酸岩反应热计算结果。
图4是本发明中酸液有效作用距离计算结果图。
具体实施方式
下面结合附图和应用实例进一步说明本发明，旨在对本发明做示意性说明和解释，并不限定本发明的范围。具体如下：
1.基于目标工区某探井X井酸压设计基础参数(表1)，建立X井井筒温度场物理模型并划分网格(见图1)，采用步骤(1)建立的井筒温度场模型计算该井施工过程中动态井底温度，得到施工结束时的井底温度。对比实测井底温度曲线，可以证明该模型可靠(见图2)。
表1 X井酸压设计部分基础参数
2.以施工结束时的井底温度为初始条件，采用数值方法计算步骤(2)和步骤(3)中的裂缝酸液温度和浓度方程，得到施工结束时各单元酸液温度和浓度。
3.基于具体实施方式中的步骤2计算得到的各单元的酸液温度和浓度计算各单元的酸岩反应热，如图3所示。
4.将具体实施方式中的步骤3计算的各单元酸岩反应热代入步骤2中，计算考虑酸岩反应热的酸液温度以及酸液浓度，得到缝长方向各单元酸液浓度分布曲线，取残酸极限浓度对应的裂缝缝长为酸液有效作用距离，如图4所示。
Claims (6)
1.深层石灰岩储层酸压酸液有效作用距离的数值计算方法，依次包括以下步骤：
(1)采用Eickmeier井筒温度场模型模拟井筒温度场并计算施工结束时的井底温度；
(2)采用PKN公式计算施工结束时的裂缝缝长，在缝长方向划分网格单元并计算各单元的平均缝宽，再计算各单元酸液流速；
(3)以步骤(1)计算的井底温度为初始条件，通过步骤(2)计算的酸液流速计算各单元不考虑酸岩反应热的酸液温度以及酸液浓度；
(4)通过步骤(3)中计算的酸液温度和浓度计算各单元的酸岩摩尔反应热；
(5)将步骤(4)计算的酸岩摩尔反应热代入步骤(3)中，计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，建立酸液浓度分布曲线，定义酸液初始注入浓度的10％作为残酸极限浓度，根据分布曲线以残酸极限浓度对应的缝长为酸液有效作用距离。
2.如权利要求1所述的深层石灰岩储层酸压酸液有效作用距离的数值计算方法，其特征在于，所述步骤(1)过程如下：
建立目标井井筒温度场模型并划分网格单元，以油管中心为轴：径向上，定义ri为网格单元半径，划分N个网格单元(i＝0,1,2,…,N)，其中模型边界点rN点的选择应满足该点处温度始终等于地层温度；纵向上，设目的层深度为H，单元深度为ΔH，则纵向上可划分为M个单元(j＝0,1,2…,M)，边界点M为(H/ΔH)取整；定义施工总时间为tT，时间步长为Δt，时间步为n(n＝0,1,2…,tT)，利用下式计算n时刻单元(i,j)的温度最终得到施工结束时的井底温度
式中：Q——注酸排量，m3/min；
ρi——网格单元i的密度，kg/m3；
Ci——网格单元i的比热容，J/(kg·℃)；
ΔHj——网格单元j的高度，m；
ri——模型径向网格半径，m；
λi——网格单元i的热传导系数，W/(m·℃)；
Δt——时间步长，min；
tT——施工总时间，min；
n——时间步，min；
——n时刻网格单元(i,j)的温度，℃。
3.如权利要求1所述的深层石灰岩储层酸压酸液有效作用距离的数值计算方法，其特征在于，所述步骤(2)过程如下：
计算施工结束即n＝tT时的裂缝缝长，定义x为缝长方向并划分网格单元，计算各单元对应的平均缝宽w：
式中：L——施工结束时的缝长，m；
c——酸液滤失系数，m/min0.5；
μf——酸液粘度，Pa·min；
h——裂缝高度，m；
E——储层岩石杨氏模量，Pa
ε——储层岩石泊松比，无因次；
x——缝长方向上单元坐标，m；
w——施工结束时各单元的平均缝宽，m；
从而得到各单元酸液流速vx：
式中：vl——各单元滤失速度，m/s；
vx——各单元x方向流速，m/s。
4.如权利要求1所述的深层石灰岩储层酸压酸液有效作用距离的数值计算方法，其特征在于，所述步骤(3)过程如下：
假设裂缝入口温度等于井底温度，设ΔrQm(Tw,p)为0，求出各单元不考虑酸岩反应热的酸液温度与裂缝壁面温度：
hT＝NNuKhf/w
hT(Tw-T)＝q(tT)-kRCmΔrQm(Tw,p)
式中：T——各单元酸液平均温度，℃；
Tw——各单元裂缝壁面温度，℃；
hT——对流换热系数，J/(m2·min·℃)；
NNu——努谢尔特数，无因次；
Khf——酸液的导热系数，J/(m·min·℃)；
ρf——储层岩石密度，kg/m3；
Cf——储层岩石比热容，J/(kg·℃)；
kR——酸岩反应速率常数，(kg/m3)-m·mol/(m2·s)；
C——各单元酸液浓度，kg/m3；
ΔrQm(Tw,p)——酸岩摩尔反应热，kJ/mol；
m——反应级数，无因次；
q(tT)——施工结束地层传热热量，kJ；
Mma——地层岩石体积热容，J/m3·℃；
Khr——地层岩石导热系数，W/(m·℃)；
——孔隙度，％；
通过下式得到各单元酸液浓度与裂缝壁面浓度：
kg(C-Cw)＝kRCw m
kg＝DeNSh/w
式中：kg——对流传质系数，m/min；
Cw——各单元裂缝壁面酸液浓度，kg/m3；
De——酸液氢离子有效传质系数，m2/s；
Nsh——舍伍德数，无因次；
k0——频率因子，(kg/m3)-m·mol/(m2·s)；
Ea——反应活化能，J/mol；
R——气体常数，8.314J/(mol·℃)。
5.如权利要求1所述的深层石灰岩储层酸压酸液有效作用距离的数值计算方法，其特征在于，所述步骤(4)过程如下：
首先计算高温高压条件酸岩摩尔反应热ΔrHm(Tw,p)：
式中：ΔrHm(Tw)——温度Tw下的摩尔反应热，kJ/mol；
p——地层压力，MPa；
Cp,m(i,T)——组分i在温度T下的常压摩尔热容，J/(℃·mol)；
fg——CO2摩尔分数，无因次；
Vco2——CO2摩尔体积，m3/mol；
Vacid——残酸体积，m3；
Sco2——残酸中CO2溶解度，m3/m3；
ΔrHm(Tw,p)——温度Tw压力p下的摩尔反应热，kJ/mol；
从而得到地层条件下酸岩摩尔反应热ΔrQm(Tw,p)：
6.如权利要求1所述的深层石灰岩储层酸压酸液有效作用距离的数值计算方法，其特征在于，所述步骤(5)过程如下：
1)通过步骤(4)计算的酸岩摩尔反应热，由步骤(3)计算考虑酸岩反应热的各单元酸液温度以及酸液浓度，并作缝长方向各单元酸液浓度分布曲线；
2)定义酸液初始注入浓度的10％作为残酸极限浓度，取残酸极限浓度对应的裂缝缝长为酸液有效作用距离。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910009184.0A CN109611067B (zh) | 2019-01-04 | 2019-01-04 | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910009184.0A CN109611067B (zh) | 2019-01-04 | 2019-01-04 | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109611067A true CN109611067A (zh) | 2019-04-12 |
| CN109611067B CN109611067B (zh) | 2020-10-09 |
Family
ID=66015624
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910009184.0A Active CN109611067B (zh) | 2019-01-04 | 2019-01-04 | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109611067B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
| CN110929447A (zh) * | 2019-12-20 | 2020-03-27 | 西南石油大学 | 一种稠化酸酸化过程中井筒温度场数值计算方法 |
| CN111197476A (zh) * | 2020-01-08 | 2020-05-26 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
| CN112345436A (zh) * | 2020-09-22 | 2021-02-09 | 中国石油天然气股份有限公司 | 酸蚀裂缝导流能力分布测试方法 |
| CN113671148A (zh) * | 2020-05-14 | 2021-11-19 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN116838309A (zh) * | 2023-06-15 | 2023-10-03 | 长江大学 | 一种碳酸盐岩储层酸压裂缝有效长度测定方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102022105A (zh) * | 2010-10-21 | 2011-04-20 | 中国石油化工股份有限公司 | 缝洞型碳酸盐岩储层大型复合酸压方法 |
| CN104179488A (zh) * | 2013-05-27 | 2014-12-03 | 中国石油化工股份有限公司 | 一种提高开发低渗透碳酸盐岩稠油油藏效果的方法 |
| EP1904601B1 (en) * | 2005-06-13 | 2015-05-06 | Sun Drilling Products Corp. | Thermoset particles with enhanced crosslinking, processing for their production, and their use in oil and natural gas drilling applications |
| CN105156081A (zh) * | 2014-06-10 | 2015-12-16 | 中国石油化工股份有限公司 | 一种碳酸盐岩稠油油藏酸化模拟评价方法 |
| CN107066769A (zh) * | 2017-06-09 | 2017-08-18 | 西南石油大学 | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 |
| US20180320054A1 (en) * | 2014-04-17 | 2018-11-08 | Saudi Arabian Oil Company | Compositions and Methods For Enhanced Fracture Cleanup Using Redox Treatment |
-
2019
- 2019-01-04 CN CN201910009184.0A patent/CN109611067B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1904601B1 (en) * | 2005-06-13 | 2015-05-06 | Sun Drilling Products Corp. | Thermoset particles with enhanced crosslinking, processing for their production, and their use in oil and natural gas drilling applications |
| CN102022105A (zh) * | 2010-10-21 | 2011-04-20 | 中国石油化工股份有限公司 | 缝洞型碳酸盐岩储层大型复合酸压方法 |
| CN104179488A (zh) * | 2013-05-27 | 2014-12-03 | 中国石油化工股份有限公司 | 一种提高开发低渗透碳酸盐岩稠油油藏效果的方法 |
| US20180320054A1 (en) * | 2014-04-17 | 2018-11-08 | Saudi Arabian Oil Company | Compositions and Methods For Enhanced Fracture Cleanup Using Redox Treatment |
| CN105156081A (zh) * | 2014-06-10 | 2015-12-16 | 中国石油化工股份有限公司 | 一种碳酸盐岩稠油油藏酸化模拟评价方法 |
| CN107066769A (zh) * | 2017-06-09 | 2017-08-18 | 西南石油大学 | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 |
Non-Patent Citations (1)
| Title |
|---|
| 刘壮: "提升高温气井酸压有效缝长方法_以川西地区栖霞组为例", 《天然气地球科学》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
| CN110609974B (zh) * | 2019-09-09 | 2022-09-16 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
| CN110929447A (zh) * | 2019-12-20 | 2020-03-27 | 西南石油大学 | 一种稠化酸酸化过程中井筒温度场数值计算方法 |
| CN111197476A (zh) * | 2020-01-08 | 2020-05-26 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
| CN111197476B (zh) * | 2020-01-08 | 2020-10-09 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
| CN113671148A (zh) * | 2020-05-14 | 2021-11-19 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN113671148B (zh) * | 2020-05-14 | 2024-05-28 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN112345436A (zh) * | 2020-09-22 | 2021-02-09 | 中国石油天然气股份有限公司 | 酸蚀裂缝导流能力分布测试方法 |
| CN112345436B (zh) * | 2020-09-22 | 2023-07-25 | 中国石油天然气股份有限公司 | 酸蚀裂缝导流能力分布测试方法 |
| CN116838309A (zh) * | 2023-06-15 | 2023-10-03 | 长江大学 | 一种碳酸盐岩储层酸压裂缝有效长度测定方法 |
| CN116838309B (zh) * | 2023-06-15 | 2024-03-19 | 长江大学 | 一种碳酸盐岩储层酸压裂缝有效长度测定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109611067B (zh) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109611067A (zh) | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 | |
| CN110334431B (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
| CN104989341B (zh) | 一种确定低渗透油藏有效驱替注采井距的方法 | |
| CN105354639B (zh) | 致密油多重介质耦合渗流的全周期产能预测方法及其装置 | |
| CN108868748A (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
| CN110188503B (zh) | 一种致密油藏产能评价方法 | |
| CN105631754B (zh) | 一种确定海外油田的产量剖面数据的方法和装置 | |
| CN107829718A (zh) | 基于均衡水驱理念的油藏井网及注采方案优化设计方法 | |
| CN109577929B (zh) | 一种超低渗致密油藏水平井建立有效驱替的定量评价方法 | |
| CN103670369A (zh) | 一种注采井间连通状况的判别方法及装置 | |
| CN106547973A (zh) | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 | |
| CN104975827B (zh) | 预测二氧化碳驱油藏指标的物质平衡方法 | |
| CN106156439A (zh) | 一种具有泡沫油现象的稠油油藏溶解气驱数值模拟方法 | |
| CN103161436A (zh) | 一种稠油热采水平井试井解释方法 | |
| CN102536223B (zh) | 煤层裂隙孔隙度快速计算方法 | |
| CN111950112B (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
| CN110173252A (zh) | 一种水平井开发气藏井网的加密方法及其评价方法 | |
| CN107288595A (zh) | 一种注水利用率的评价方法 | |
| CN106203699A (zh) | 一种粗糙裂缝初始导流能力的预测方法 | |
| Ma et al. | Production performance optimization of water-enhanced geothermal system based on effective time-averaged volume | |
| Meng et al. | Experimental research on three dimensional flow characteristics of multiple horizontal fractures utilizing an innovative experiment apparatus | |
| CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
| Rodriguez | Inferences of two dynamic processes on recovery factor and well spacing for a shale oil reservoir | |
| CN104747154A (zh) | 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法 | |
| CN109653717A (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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200611 Address after: No. 8 Road, Xindu Xindu District of Chengdu city of Sichuan Province in 610000 Applicant after: SOUTHWEST PETROLEUM University Applicant after: CHINA PETROLEUM AND NATURAL GAS Co.,Ltd. SOUTHWEST OIL AND GAS FIELD BRANCH CHUANZHONG OIL AND GAS MINE Address before: No. 8 Road, Xindu Xindu District of Chengdu city of Sichuan Province in 610500 Applicant before: SOUTHWEST PETROLEUM University |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |