CN110929447B - 一种稠化酸酸化过程中井筒温度场数值计算方法 - Google Patents

Abstract

本发明涉及一种稠化酸酸化过程中井筒温度场数值计算方法，包括：（1）利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型；（2）判定不同注入条件下的稠化酸流动流态，计算不同流态下的稠化酸井筒流动剪切速率；（3）将稠化酸井筒流动剪切速率带入流变模型，得到不同注入条件下的酸液表观粘度，计算不同注入条件下的稠化酸对流换热系数；（4）建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；（5）计算不同注入条件下的井筒内酸液温度分布。本发明考虑酸化过程中酸液流变性对于井筒‑地层传热效率的影响，采用数值方法计算井筒内酸液温度分布，对于海相碳酸盐岩储层酸化改造设计具有指导意义。

Description

一种稠化酸酸化过程中井筒温度场数值计算方法

技术领域

本发明涉及石油工程领域，尤其是深层碳酸盐岩储层稠化酸酸化过程中考虑酸液流变性影响的井筒温度场的数值计算方法。采用数值计算方法计算酸化过程中井筒流动酸液温度变化，提高深层碳酸盐岩储层酸化设计的准确性，实现深层碳酸盐岩储层高效酸化改造。

背景技术

深层海相碳酸盐岩层系是实现四川盆地天然气效益勘探开发的主力层系，储层具有埋藏深(>5000m)、温度高(>140℃)等特征(马新华,杨雨,文龙,罗冰.四川盆地海相碳酸盐岩大中型气田分布规律及勘探方向[J].石油勘探与开发,2019,46(1):1-13)。该类储层开发难以自然获得工业油气流，酸化技术是区块储层主要的增产改造手段之一。

酸化是一种通过向地层注入一定类型、浓度的酸液和添加剂组成的配方酸液，溶蚀地层岩石部分矿物或孔隙、裂缝内的堵塞物，提高地层或裂缝渗透性，改善渗流条件，达到恢复或提高油气井产能目的的油气田增产技术(李颖川.采油工程[M].石油工业出版社,2009)。稠化酸是目前广泛应用于海相碳酸盐岩储层酸化改造的主要酸液体系之一(韩慧芬.四川盆地上震旦统灯四段气藏提高单井产量的技术措施[J].天然气工业,2017,37(08):40-47)。酸化过程中，酸液在井筒中流动会受到高温地层加热而持续升温，酸液温度直接影响酸液性能以及酸与岩石的反应速率，进而影响储层酸化效果。因此，准确预测井筒温度场有利于优化储层酸化工艺与设计。

目前常用的井筒温度计算方法主要关注井筒与地层的传热过程分析(Zhang Z,Xiong Y M,Gao Y,et al.Wellbore temperature distribution during circulationstage when well-kick occurs in a continuous formation from the bottom-hole.Energy,164(2018):964-977,2018)，却忽略了酸化过程中酸液流变性对于井筒-地层传热效果的影响。因此，针对四川盆地高温深井的特点，研究酸化过程井筒温度变化规律，考虑酸液流变性对于井筒内酸液温度分布的影响，对于四川海相碳酸盐岩储层酸化改造设计至关重要。

发明内容

本发明的目的在于提供一种稠化酸酸化过程中井筒温度场数值计算方法，该方法重点研究酸化注液过程中酸液流变性对于井筒-地层传热效率的影响，采用数值方法计算考虑上述影响的井筒内酸液温度分布。该方法原理可靠，计算高效，有利于计算高温深井井筒内酸液温度，对于四川海相碳酸盐岩储层酸化改造设计具有指导意义。

为达到以上技术目的，本发明采用以下技术方案。

首先，利用流变仪测定不同温度下现场常用稠化酸液样品的表观粘度；其次，基于非牛顿流动流态判别方法判定井筒内稠化酸流动流态，并建立不同流态下的井筒内稠化酸流动剪切速率模型；再者，结合稠化酸流变模型以及流动剪切速率模型求解不同流态下的稠化酸对流换热系数；然后，结合对流换热系数，建立考虑酸液流变性影响的井筒温度场数值计算模型；最后，基于井筒温度场模型计算井筒内酸液温度分布分布。

一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：

(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型；

(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，基于解析公式以及ANSYS软件，计算不同流态下的稠化酸井筒流动剪切速率；

(3)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数；

(4)建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；

(5)基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，计算不同注入条件下的井筒内酸液温度分布。

在本发明中，所述步骤(1)中，利用流变仪测定稠化酸样品不同温度不同剪切速率下的表观粘度，建立稠化酸流变模型，过程如下：

选取现场常用的聚丙烯酰胺稠化酸体系，配制等量稠化酸作为流变实验样品，利用流变仪开展不同实验温度下液体变剪切表观粘度测试，获取酸液流变曲线，采用最小二乘法拟合酸液流变曲线，建立考虑温度和流动剪切速率共同影响的酸液流变模型：

μ_a＝ke^aTγ^n-1 (1)

式中：μ_a——酸液表观粘度，Pa·s；

k——稠度系数，Pa·sⁿ；

a——常数，无因次；

T——酸液温度，℃；

γ——剪切速率，s^-1；

n——流态指数，无因次。

在本发明中，所述步骤(2)中，基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，基于解析公式以及ANSYS软件，计算不同流态下的稠化酸井筒流动剪切速率，过程如下：

1)基于刘崇建等推导的非牛顿流体雷诺数计算公式(刘崇建,刘孝良,柳世杰.非牛顿流体流态判别方法的研究[J].天然气工业,2001(04):49-52)，计算不同注入条件下井筒内流体流动雷诺数：

式中：Re——流动雷诺数，无因次；

D——井筒直径，m；

v——酸液注入速度，m/s；

ρ——酸液密度，kg/m³。

2)基于刘乃震等推到的紊流流动临界雷诺数公式(刘乃震,王廷瑞,刘孝良,等.非牛顿流体的稳定性及其流态判别[J].天然气工业,2003(01):53-57)，计算不同注入条件下稠化酸井筒流动临界雷诺数：

式中：Re_c——临界雷诺数，无因次。

3)对比步骤1)和步骤2)计算的雷诺数与临界雷诺数，判断井筒内酸液流动流态：

4)基于步骤3)判定的结果，若井筒内酸液流态为层流，采用刘崇建等推导的层流剪切速率解析公式计算稠化酸井筒流动剪切速率γ：

5)基于步骤3)判定的结果，若井筒内酸液流态为紊流，利用ANSYS软件建立井筒流体流动模型并计算稠化酸井筒流动剪切速率γ。其具体操作步骤为：①利用ANSYS-ICEM软件建立长为L，直径为D的井筒流动物理模型(其中L和D参考实际井筒尺寸，要求设置L>50D以保证流体流动充分发展)；②对物理模型自动进行网格划分，井筒长度方向采用矩形网格，井筒径向方向采用“Y”型网格，并将划分好网格的物理模型保存为“.mesh”文件；③将“.mesh”文件导入ANSYS-Fluent软件，并在“General”模块对导入的网格进行网格质量检测，若网格质量低于0.3，须返回步骤①重新划分网格；若高于0.3，则继续下一步骤；④在“Models”模块选择“K-epsilon(2eqn)”紊流计算模型；⑤在“Materials”模块选中“non-newtonian-power-law”类型流体，并基于步骤(1)计算的流变参数赋予该流体稠度系数k以及流态指数n；⑥在“Cell Zone Conditions”模块中设置模型边界条件：“inlet”设置为“velocity-inlet”并赋予入口注入速度v，“outlet”设置为“pressure-outlet”并赋予出口压力(默认为大气压)；⑦在“Reference Value”模块中将步骤(1)中的稠化酸实验参数赋为流体参数；⑧在“Run Calculation”中设置迭代次数(默认为500次)，模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ。

所述模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ，过程如下：

基于模拟的稠化酸流动剪切速率曲线，拟合得到不同注入条件下的紊流剪切速率计算公式，通过该公式计算紊流状态下的酸液剪切速率：

式中：c₁，c₂——拟合经验参数，无因次。

在本发明中，所述步骤(3)中，将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数，过程如下：

1)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，可以计算得到不同注入条件下(温度，剪切速率)的酸液表观粘度；

2)若井筒内酸液流态为层流，将步骤1)计算的表观粘度代入Rohsenow等给出的层流状态下对流换热系数计算公式(Rohsenow W,Harnett J,Cho Y.Handbook of heattransfer(3rd ed).New York:McGraw Hill book company,1998)，计算层流状态下稠化酸井筒流动对流换热系数：

h＝1.86λ(RePrD/L)^0.33(μ_at/μ_aw)^0.14/D (7)

式中：h——对流换热系数，W/(m²·℃)；

λ——酸液导热系数，W/(m·℃)；

L——井筒长度，m；

Pr——Prandtl数，无因次；

μ_at——油管中(井筒内)流动酸液的表观粘度，Pa·s；

μ_aw——油管壁面温度下的酸液表观粘度，Pa·s；

3)若井筒内酸液流态为紊流，将步骤1)计算的表观粘度代入Gniellinski等建立的紊流状态下对流换热系数计算公式(Gnielinski V.New equations for heat and masstransfer in turbulent pipe and channel flows.International ChemicalEngineering,1976.16(2):359-368)，计算得到紊流状态下稠化酸井筒流动对流换热系数：

公式(7)、公式(8)中Pr数的定义为：

Pr＝μ_atC/λ (9)

C——酸液比热容，J/(kg·℃)。

在本发明中，所述步骤(4)中，建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型，过程如下：

1)建立酸化井筒温度场模型(廖毅.非均质碳酸盐岩转向酸化模型研究[D].西南石油大学,2017)：

式中：Q——注酸排量，m³/s；

T₀——油管内酸液温度，℃；

z——井筒垂向深度，m；

T₁——油管温度，℃；

T₂——环空液温度，℃；

T₃——套管温度，℃；

T_i——第i层管柱温度，℃；

r₀——油管内壁半径，m；

r₁——油管半径，m；

r₂——套管内壁半径，m；

r₃——套管半径，m；

r_i——第i层管柱半径，m；

h_a,h_w——酸液，环空液对流换热系数，W/(m²·℃)；

ρ_t——油管密度，kg/m³；

C_t——油管比热容，J/(kg·℃)；

ρ_w——环空液密度，kg/m³；

C_w——环空液比热容，J/(kg·℃)；

λ_i——第i层管柱导热系数，W/(m·℃)；

ρ_i——第i层管柱密度，kg/m³；

C_i——第i层管柱比热容，J/(kg·℃)。

需要说明的是，上述模型中酸液对流换热系数由恒定酸液粘度计算所得，未考虑酸液流变性影响，计算过程中酸液对流换热系数为常值。

2)将步骤(3)得到的稠化酸井筒流动对流换热系数代入步骤1)建立的酸化井筒温度场模型中，得到考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型：

若井筒酸液流态为层流(Re≤Re_c)：

若井筒酸液流态为紊流(Re＞Re_c)：

在本发明中，所述步骤(5)中，基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场模型，计算不同注入条件下的井筒内酸液温度分布，内容如下：

基于步骤(4)建立的稠化酸酸化井筒温度场模型，带入不同的注入参数即可求解不同注入条件下的井筒垂直深度z对应的油管内酸液温度T₀。

与现有技术相比，本方法具有以下有益效果：

①本发明基于实验方法建立了一套稠化酸流变模型，模型简单，理论可靠，能够计算指定注入条件下的酸液表观粘度；

②本发明建立了考虑稠化酸流变性影响的酸化井筒温度场模型，克服了采用经验常数判断流态以及采用实验测得的稳态酸液粘度计算井筒传热的误差；

③本发明考虑了稠化酸流变性对于稠化酸对流传热效率的影响，参数来源真实，预测结果更加符合现场实际情况。

附图说明

图1是本发明中不同温度下稠化酸流变曲线。

图2是本发明中井筒流体流动物理模型。

图3是本发明中稠化酸紊流流动(酸液温度50℃)不同注入速率条件下的剪切速率。

图4是本发明中稠化酸酸化井筒垂向深度—酸液温度分布图。

具体实施方式

下面结合附图和应用实例进一步说明本发明，旨在对本发明做示意性说明和解释，并不限定本发明的范围。

实施例1

一种稠化酸酸化过程中井筒温度场数值计算方法，具体如下：

1.基于川西海相碳酸盐岩气藏某井现场常用的一套稠化酸体系配制实验样品，现场酸化注酸速率为2m³/min，地温梯度为0.22℃/m，井筒直径为90mm，井深为4800m。利用流变仪测试不同温度下的稠化酸流变曲线(见图1)，采用步骤(1)所述方法，计算式(1)中各参数为：k＝3.1，a＝－0.003，n＝0.29。建立不同温度下的稠化酸流变模型。

2.将步骤(1)获得的相关参数带入步骤(2)中，可以计算当前注入条件下的稠化酸流动雷诺数为7571.2，计算流动临界雷诺数为3692.5。

3.由步骤2可知当前注入条件下稠化酸井筒流动为紊流流态，基于步骤(3)所述方法以及现场参数，利用ANSYS软件建立井筒流体流动物理模型(见图2)，模拟不同注入条件下稠化酸井筒紊流流动的剪切速率(见图3)。代入步骤2所得参数，计算当前注入条件下流动剪切速率为320s^-1，式(6)各参数为：c₁＝52.3，c₂＝0.96。

4.将步骤3建立的剪切速率模型引入步骤1建立的稠化酸流变模型可以计算特定注入条件下酸液的表观粘度，带入式(7)可以计算紊流状态下的稠化酸对流换热系数。

5.将步骤4计算的对流换热系数带入步骤(4)所述模型中替换原有对流换热系数，建立稠化酸酸化井筒温度场模型。

6.基于步骤5建立的稠化酸酸化井筒温度场模型，模拟当前注入条件下井筒内酸液温度分布(见图4)。

Claims (6)

1.一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：

(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型，过程如下：

利用流变仪开展不同实验温度下液体变剪切表观粘度测试，获取酸液流变曲线，采用最小二乘法拟合酸液流变曲线，建立考虑温度和流动剪切速率共同影响的酸液流变模型：

μ_a＝ke^aTγ^n-1

式中：μ_a——酸液表观粘度，Pa·s；

k——稠度系数，Pa·sⁿ；

a——常数，无因次；

T——酸液温度，℃；

γ——剪切速率，s^-1；

n——流态指数，无因次；

(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，计算不同流态下的稠化酸井筒流动剪切速率，过程如下：

计算不同注入条件下井筒内流体流动雷诺数Re：

式中：D——井筒直径，m；

v——酸液注入速度，m/s；

ρ——酸液密度，kg/m³；

计算不同注入条件下稠化酸井筒流动临界雷诺数Re_c：

判断井筒内酸液流动流态：

若井筒内酸液流态为层流，利用下式计算稠化酸井筒流动剪切速率γ：

若井筒内酸液流态为紊流，利用ANSYS软件建立井筒流体流动模型并计算稠化酸井筒流动剪切速率γ；

(3)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数；

(4)建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；

(5)基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，计算不同注入条件下的井筒内酸液温度分布。

2.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(3)过程如下：

1)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，计算得到不同注入条件下的酸液表观粘度；

2)若井筒内酸液流态为层流，利用下式计算层流状态下稠化酸井筒流动对流换热系数：

h＝1.86λ(RePrD/L)^0.33(μ_at/μ_aw)^0.14/D

若井筒内酸液流态为紊流，利用下式计算紊流状态下稠化酸井筒流动对流换热系数：

Pr＝μ_atC/λ

式中：λ——酸液导热系数，W/(m·℃)；

L——井筒长度，m；

Pr——Prandtl数，无因次；

μ_at——油管中(井筒内)流动酸液的表观粘度，Pa·s；

μ_aw——油管壁面温度下的酸液表观粘度，Pa·s；

C——酸液比热容，J/(kg·℃)。

3.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(4)过程如下：

1)建立酸化井筒温度场模型：

式中：Q——注酸排量，m³/s；

T₀——油管内酸液温度，℃；

z——井筒垂向深度，m；

T₁——油管温度，℃；

T₂——环空液温度，℃；

T₃——套管温度，℃；

T_i——第i层管柱温度，℃；

r₀——油管内壁半径，m；

r₁——油管半径，m；

r₂——套管内壁半径，m；

r₃——套管半径，m；

r_i——第i层管柱半径，m；

h_a,h_w——酸液，环空液对流换热系数，W/(m²·℃)；

ρ_t——油管密度，kg/m³；

C_t——油管比热容，J/(kg·℃)；

ρ_w——环空液密度，kg/m³；

C_w——环空液比热容，J/(kg·℃)；

λ_i——第i层管柱导热系数，W/(m·℃)；

ρ_i——第i层管柱密度，kg/m³；

C_i——第i层管柱比热容，J/(kg·℃)；

2)将步骤(3)得到的稠化酸井筒流动对流换热系数代入酸化井筒温度场模型中，得到考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型：

若井筒酸液流态为层流：

若井筒酸液流态为紊流：

式中：h——对流换热系数，W/(m²·℃)；

λ——酸液导热系数，W/(m·℃)；

L——井筒长度，m；

C——酸液比热容，J/(kg·℃)。

4.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(5)内容如下：

基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，带入不同的注入参数即可求解不同注入条件下的井筒垂直深度z对应的油管内酸液温度T₀。

5.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，若井筒内酸液流态为紊流，利用ANSYS软件建立井筒流体流动模型并计算稠化酸井筒流动剪切速率γ，具体步骤为：

①利用ANSYS-ICEM软件建立长为L，直径为D的井筒流动物理模型；②对物理模型自动进行网格划分，井筒长度方向采用矩形网格，井筒径向方向采用“Y”型网格，并将划分好网格的物理模型保存为“.mesh”文件；③将“.mesh”文件导入ANSYS-Fluent软件，并在“General”模块对导入的网格进行网格质量检测，若网格质量低于0.3，须返回步骤①重新划分网格；④在“Models”模块选择“K-epsilon(2eqn)”紊流计算模型；⑤在“Materials”模块选中“non-newtonian-power-law”类型流体；⑥在“Cell Zone Conditions”模块中设置模型边界条件：“inlet”设置为“velocity-inlet”并赋予入口注入速度v，“outlet”设置为“pressure-outlet”并赋予出口压力；⑦在“Reference Value”模块中将稠化酸实验参数赋为流体参数；⑧在“Run Calculation”中设置迭代次数，模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ。

6.如权利要求5所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ，过程如下：

基于模拟的稠化酸流动剪切速率曲线，拟合得到不同注入条件下的紊流剪切速率计算公式，计算紊流状态下的酸液剪切速率：

式中：c₁，c₂——拟合经验参数，无因次。

Images