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

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

A numerical calculation method of wellbore temperature field in the process of acid thickening

【技术实现步骤摘要】
一种稠化酸酸化过程中井筒温度场数值计算方法
本专利技术涉及石油工程领域，尤其是深层碳酸盐岩储层稠化酸酸化过程中考虑酸液流变性影响的井筒温度场的数值计算方法。采用数值计算方法计算酸化过程中井筒流动酸液温度变化，提高深层碳酸盐岩储层酸化设计的准确性，实现深层碳酸盐岩储层高效酸化改造。
技术介绍
深层海相碳酸盐岩层系是实现四川盆地天然气效益勘探开发的主力层系，储层具有埋藏深(>5000m)、温度高(>140℃)等特征(马新华,杨雨,文龙,罗冰.四川盆地海相碳酸盐岩大中型气田分布规律及勘探方向[J].石油勘探与开发,2019,46(1):1-13)。该类储层开发难以自然获得工业油气流，酸化技术是区块储层主要的增产改造手段之一。酸化是一种通过向地层注入一定类型、浓度的酸液和添加剂组成的配方酸液，溶蚀地层岩石部分矿物或孔隙、裂缝内的堵塞物，提高地层或裂缝渗透性，改善渗流条件，达到恢复或提高油气井产能目的的油气田增产技术(李颖川.采油工程[M].石油工业出版社,2009)。稠化酸是目前广泛应用于海相碳酸盐岩储层酸化改造的主要酸液体系之一(韩慧芬.四川盆地上震旦统灯四段气藏提高单井产量的技术措施[J].天然气工业,2017,37(08):40-47)。酸化过程中，酸液在井筒中流动会受到高温地层加热而持续升温，酸液温度直接影响酸液性能以及酸与岩石的反应速率，进而影响储层酸化效果。因此，准确预测井筒温度场有利于优化储层酸化工艺与设计。目前常用的井筒温度计算方法主要关注井筒与地层的传热过程分析(ZhangZ,XiongYM,GaoY,etal.Wellboretemperaturedistributionduringcirculationstagewhenwell-kickoccursinacontinuousformationfromthebottom-hole.Energy,164(2018):964-977,2018)，却忽略了酸化过程中酸液流变性对于井筒-地层传热效果的影响。因此，针对四川盆地高温深井的特点，研究酸化过程井筒温度变化规律，考虑酸液流变性对于井筒内酸液温度分布的影响，对于四川海相碳酸盐岩储层酸化改造设计至关重要。
技术实现思路
本专利技术的目的在于提供一种稠化酸酸化过程中井筒温度场数值计算方法，该方法重点研究酸化注液过程中酸液流变性对于井筒-地层传热效率的影响，采用数值方法计算考虑上述影响的井筒内酸液温度分布。该方法原理可靠，计算高效，有利于计算高温深井井筒内酸液温度，对于四川海相碳酸盐岩储层酸化改造设计具有指导意义。为达到以上技术目的，本专利技术采用以下技术方案。首先，利用流变仪测定不同温度下现场常用稠化酸液样品的表观粘度；其次，基于非牛顿流动流态判别方法判定井筒内稠化酸流动流态，并建立不同流态下的井筒内稠化酸流动剪切速率模型；再者，结合稠化酸流变模型以及流动剪切速率模型求解不同流态下的稠化酸对流换热系数；然后，结合对流换热系数，建立考虑酸液流变性影响的井筒温度场数值计算模型；最后，基于井筒温度场模型计算井筒内酸液温度分布分布。一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型；(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，基于解析公式以及ANSYS软件，计算不同流态下的稠化酸井筒流动剪切速率；(3)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数；(4)建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；(5)基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，计算不同注入条件下的井筒内酸液温度分布。在本专利技术中，所述步骤(1)中，利用流变仪测定稠化酸样品不同温度不同剪切速率下的表观粘度，建立稠化酸流变模型，过程如下：选取现场常用的聚丙烯酰胺稠化酸体系，配制等量稠化酸作为流变实验样品，利用流变仪开展不同实验温度下液体变剪切表观粘度测试，获取酸液流变曲线，采用最小二乘法拟合酸液流变曲线，建立考虑温度和流动剪切速率共同影响的酸液流变模型：μa＝keaTγn-1(1)式中：μa——酸液表观粘度，Pa·s；k——稠度系数，Pa·sn；a——常数，无因次；T——酸液温度，℃；γ——剪切速率，s-1；n——流态指数，无因次。在本专利技术中，所述步骤(2)中，基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，基于解析公式以及ANSYS软件，计算不同流态下的稠化酸井筒流动剪切速率，过程如下：1)基于刘崇建等推导的非牛顿流体雷诺数计算公式(刘崇建,刘孝良,柳世杰.非牛顿流体流态判别方法的研究[J].天然气工业,2001(04):49-52)，计算不同注入条件下井筒内流体流动雷诺数：式中：Re——流动雷诺数，无因次；D——井筒直径，m；v——酸液注入速度，m/s；ρ——酸液密度，kg/m3。2)基于刘乃震等推到的紊流流动临界雷诺数公式(刘乃震,王廷瑞,刘孝良,等.非牛顿流体的稳定性及其流态判别[J].天然气工业,2003(01):53-57)，计算不同注入条件下稠化酸井筒流动临界雷诺数：式中：Rec——临界雷诺数，无因次。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；⑥在“CellZoneConditions”模块中设置模型边界条件：“inlet”设置为“velocity-inlet”并赋予入本文档来自技高网...

【技术保护点】
1.一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：/n(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型；/n(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，计算不同流态下的稠化酸井筒流动剪切速率；/n(3)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数；/n(4)建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；/n(5)基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，计算不同注入条件下的井筒内酸液温度分布。/n

【技术特征摘要】
1.一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：
(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型；
(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，计算不同流态下的稠化酸井筒流动剪切速率；
(3)将步骤(2)计算的稠化酸井筒流动剪切速率带入步骤(1)建立的流变模型，得到不同注入条件下的酸液表观粘度，结合不同流态下的对流换热系数公式，计算不同注入条件下的稠化酸对流换热系数；
(4)建立考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型；
(5)基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，计算不同注入条件下的井筒内酸液温度分布。


2.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(1)过程如下：
利用流变仪开展不同实验温度下液体变剪切表观粘度测试，获取酸液流变曲线，采用最小二乘法拟合酸液流变曲线，建立考虑温度和流动剪切速率共同影响的酸液流变模型：
μa＝keaTγn-1
式中：μa——酸液表观粘度，Pa·s；
k——稠度系数，Pa·sn；
a——常数，无因次；
T——酸液温度，℃；
γ——剪切速率，s-1；
n——流态指数，无因次。


3.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(2)过程如下：
计算不同注入条件下井筒内流体流动雷诺数Re：



式中：D——井筒直径，m；
v——酸液注入速度，m/s；
ρ——酸液密度，kg/m3；
计算不同注入条件下稠化酸井筒流动临界雷诺数Rec：



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



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



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


4.如权利要求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·℃)。


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


...

【专利技术属性】
技术研发人员：郭建春，刘壮，路千里，李鸣，任冀川，陶亮，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川;51