CN110929447B - 一种稠化酸酸化过程中井筒温度场数值计算方法 - Google Patents
一种稠化酸酸化过程中井筒温度场数值计算方法 Download PDFInfo
- Publication number
- CN110929447B CN110929447B CN201911326643.4A CN201911326643A CN110929447B CN 110929447 B CN110929447 B CN 110929447B CN 201911326643 A CN201911326643 A CN 201911326643A CN 110929447 B CN110929447 B CN 110929447B
- Authority
- CN
- China
- Prior art keywords
- acid
- wellbore
- flow
- thickening
- under different
- 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.)
- Active
Links
- 239000002253 acid Substances 0.000 title claims abstract description 183
- 230000008719 thickening Effects 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 title claims abstract description 36
- 230000008569 process Effects 0.000 title claims abstract description 28
- 230000020477 pH reduction Effects 0.000 title claims abstract description 23
- 238000002347 injection Methods 0.000 claims abstract description 48
- 239000007924 injection Substances 0.000 claims abstract description 48
- 239000007788 liquid Substances 0.000 claims abstract description 34
- 238000012546 transfer Methods 0.000 claims abstract description 33
- 239000000243 solution Substances 0.000 claims abstract description 22
- 238000000518 rheometry Methods 0.000 claims abstract description 17
- 239000012530 fluid Substances 0.000 claims description 36
- 238000012850 discrimination method Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 230000008676 import Effects 0.000 claims 1
- 230000009916 joint effect Effects 0.000 claims 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 abstract description 11
- 238000013461 design Methods 0.000 abstract description 5
- 239000011435 rock Substances 0.000 abstract description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 7
- 239000003921 oil Substances 0.000 description 7
- 239000007789 gas Substances 0.000 description 6
- 230000000638 stimulation Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 229920002401 polyacrylamide Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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＝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；⑥在“Cell Zone Conditions”模块中设置模型边界条件：“inlet”设置为“velocity-inlet”并赋予入口注入速度v，“outlet”设置为“pressure-outlet”并赋予出口压力(默认为大气压)；⑦在“Reference Value”模块中将步骤(1)中的稠化酸实验参数赋为流体参数；⑧在“Run Calculation”中设置迭代次数(默认为500次)，模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ。
所述模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ，过程如下：
基于模拟的稠化酸流动剪切速率曲线，拟合得到不同注入条件下的紊流剪切速率计算公式，通过该公式计算紊流状态下的酸液剪切速率：
式中：c1，c2——拟合经验参数，无因次。
在本发明中，所述步骤(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/(m2·℃)；
λ——酸液导热系数，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——注酸排量，m3/s；
T0——油管内酸液温度，℃；
z——井筒垂向深度，m；
T1——油管温度，℃；
T2——环空液温度，℃；
T3——套管温度，℃；
Ti——第i层管柱温度，℃；
r0——油管内壁半径，m；
r1——油管半径，m；
r2——套管内壁半径，m；
r3——套管半径，m；
ri——第i层管柱半径，m；
ha,hw——酸液，环空液对流换热系数，W/(m2·℃)；
ρt——油管密度，kg/m3；
Ct——油管比热容，J/(kg·℃)；
ρw——环空液密度，kg/m3；
Cw——环空液比热容，J/(kg·℃)；
λi——第i层管柱导热系数，W/(m·℃)；
ρi——第i层管柱密度，kg/m3；
Ci——第i层管柱比热容，J/(kg·℃)。
需要说明的是，上述模型中酸液对流换热系数由恒定酸液粘度计算所得，未考虑酸液流变性影响，计算过程中酸液对流换热系数为常值。
2)将步骤(3)得到的稠化酸井筒流动对流换热系数代入步骤1)建立的酸化井筒温度场模型中，得到考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型：
若井筒酸液流态为层流(Re≤Rec)：
若井筒酸液流态为紊流(Re＞Rec)：
在本发明中，所述步骤(5)中，基于步骤(3)计算得到的对流换热系数，基于步骤(4)建立的稠化酸酸化井筒温度场模型，计算不同注入条件下的井筒内酸液温度分布，内容如下：
基于步骤(4)建立的稠化酸酸化井筒温度场模型，带入不同的注入参数即可求解不同注入条件下的井筒垂直深度z对应的油管内酸液温度T0。
与现有技术相比，本方法具有以下有益效果：
①本发明基于实验方法建立了一套稠化酸流变模型，模型简单，理论可靠，能够计算指定注入条件下的酸液表观粘度；
②本发明建立了考虑稠化酸流变性影响的酸化井筒温度场模型，克服了采用经验常数判断流态以及采用实验测得的稳态酸液粘度计算井筒传热的误差；
③本发明考虑了稠化酸流变性对于稠化酸对流传热效率的影响，参数来源真实，预测结果更加符合现场实际情况。
附图说明
图1是本发明中不同温度下稠化酸流变曲线。
图2是本发明中井筒流体流动物理模型。
图3是本发明中稠化酸紊流流动(酸液温度50℃)不同注入速率条件下的剪切速率。
图4是本发明中稠化酸酸化井筒垂向深度—酸液温度分布图。
具体实施方式
下面结合附图和应用实例进一步说明本发明，旨在对本发明做示意性说明和解释，并不限定本发明的范围。
实施例1
一种稠化酸酸化过程中井筒温度场数值计算方法，具体如下：
1.基于川西海相碳酸盐岩气藏某井现场常用的一套稠化酸体系配制实验样品，现场酸化注酸速率为2m3/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)各参数为：c1＝52.3，c2＝0.96。
4.将步骤3建立的剪切速率模型引入步骤1建立的稠化酸流变模型可以计算特定注入条件下酸液的表观粘度，带入式(7)可以计算紊流状态下的稠化酸对流换热系数。
5.将步骤4计算的对流换热系数带入步骤(4)所述模型中替换原有对流换热系数，建立稠化酸酸化井筒温度场模型。
6.基于步骤5建立的稠化酸酸化井筒温度场模型，模拟当前注入条件下井筒内酸液温度分布(见图4)。
Claims (6)
1.一种稠化酸酸化过程中井筒温度场数值计算方法，依次包括以下步骤：
(1)利用流变仪测定稠化酸样品不同温度、不同剪切速率下的表观粘度，建立稠化酸流变模型，过程如下：
利用流变仪开展不同实验温度下液体变剪切表观粘度测试，获取酸液流变曲线，采用最小二乘法拟合酸液流变曲线，建立考虑温度和流动剪切速率共同影响的酸液流变模型：
μa＝keaTγn-1
式中：μa——酸液表观粘度，Pa·s；
k——稠度系数，Pa·sn；
a——常数，无因次；
T——酸液温度，℃；
γ——剪切速率，s-1；
n——流态指数，无因次；
(2)基于非牛顿流体流动流态判别方法，判定不同注入条件下的稠化酸流动流态，计算不同流态下的稠化酸井筒流动剪切速率，过程如下：
计算不同注入条件下井筒内流体流动雷诺数Re：
式中：D——井筒直径，m；
v——酸液注入速度，m/s；
ρ——酸液密度，kg/m3；
计算不同注入条件下稠化酸井筒流动临界雷诺数Rec：
判断井筒内酸液流动流态：
若井筒内酸液流态为层流，利用下式计算稠化酸井筒流动剪切速率γ：
若井筒内酸液流态为紊流，利用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——注酸排量，m3/s；
T0——油管内酸液温度，℃；
z——井筒垂向深度，m；
T1——油管温度，℃；
T2——环空液温度，℃；
T3——套管温度，℃；
Ti——第i层管柱温度，℃；
r0——油管内壁半径，m；
r1——油管半径，m；
r2——套管内壁半径，m；
r3——套管半径，m；
ri——第i层管柱半径，m；
ha,hw——酸液，环空液对流换热系数，W/(m2·℃)；
ρt——油管密度，kg/m3；
Ct——油管比热容，J/(kg·℃)；
ρw——环空液密度，kg/m3；
Cw——环空液比热容，J/(kg·℃)；
λi——第i层管柱导热系数，W/(m·℃)；
ρi——第i层管柱密度，kg/m3；
Ci——第i层管柱比热容，J/(kg·℃)；
2)将步骤(3)得到的稠化酸井筒流动对流换热系数代入酸化井筒温度场模型中，得到考虑酸液流变性影响的稠化酸酸化井筒温度场数值计算模型：
若井筒酸液流态为层流：
若井筒酸液流态为紊流：
式中：h——对流换热系数，W/(m2·℃)；
λ——酸液导热系数，W/(m·℃)；
L——井筒长度，m；
C——酸液比热容，J/(kg·℃)。
4.如权利要求1所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述步骤(5)内容如下：
基于步骤(4)建立的稠化酸酸化井筒温度场数值计算模型，带入不同的注入参数即可求解不同注入条件下的井筒垂直深度z对应的油管内酸液温度T0。
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所述的一种稠化酸酸化过程中井筒温度场数值计算方法，其特征在于，所述模拟得到不同注入条件下的稠化酸流动剪切速率曲线，以此计算稠化酸井筒流动剪切速率γ，过程如下：
基于模拟的稠化酸流动剪切速率曲线，拟合得到不同注入条件下的紊流剪切速率计算公式，计算紊流状态下的酸液剪切速率：
式中：c1，c2——拟合经验参数，无因次。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911326643.4A CN110929447B (zh) | 2019-12-20 | 2019-12-20 | 一种稠化酸酸化过程中井筒温度场数值计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911326643.4A CN110929447B (zh) | 2019-12-20 | 2019-12-20 | 一种稠化酸酸化过程中井筒温度场数值计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110929447A CN110929447A (zh) | 2020-03-27 |
| CN110929447B true CN110929447B (zh) | 2020-11-03 |
Family
ID=69864508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911326643.4A Active CN110929447B (zh) | 2019-12-20 | 2019-12-20 | 一种稠化酸酸化过程中井筒温度场数值计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110929447B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112198092A (zh) * | 2020-08-24 | 2021-01-08 | 浙江南都电源动力股份有限公司 | 测试电池浆料稳定性的方法 |
| CN112502685B (zh) * | 2020-12-03 | 2022-03-11 | 西南石油大学 | 一种考虑热效应的碳酸盐岩储层交替酸压级数优选方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105154038A (zh) * | 2009-02-16 | 2015-12-16 | 切弗朗菲利浦化学公司 | 包括阳离子聚合物的井筒应用液及其使用方法 |
| CN104806230B (zh) * | 2015-02-16 | 2018-01-05 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
| CN108316251A (zh) * | 2018-02-06 | 2018-07-24 | 浙江大学 | 一种采用超高韧性水泥基复合材料加固堆石坝混凝土面板的方法 |
| CN110159239A (zh) * | 2019-05-06 | 2019-08-23 | 中国石油天然气股份有限公司 | 一种直井大规模水力压裂油套同注压裂方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA3069503A1 (en) * | 2017-07-10 | 2019-01-17 | Exxonmobil Upstream Research Company | Methods for deep reservoir stimulation using acid-forming fluids |
| US10647903B2 (en) * | 2017-08-15 | 2020-05-12 | Saudi Arabian Oil Company | Oil-based drilling fluid compositions which include layered double hydroxides as rheology modifiers and amino amides as emulsifiers |
| CN109057762B (zh) * | 2018-07-23 | 2019-08-23 | 中国石油大学(北京) | 一种碳酸盐岩油气藏的酸化方法 |
| CN109611067B (zh) * | 2019-01-04 | 2020-10-09 | 西南石油大学 | 深层石灰岩储层酸压酸液有效作用距离的数值计算方法 |
-
2019
- 2019-12-20 CN CN201911326643.4A patent/CN110929447B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105154038A (zh) * | 2009-02-16 | 2015-12-16 | 切弗朗菲利浦化学公司 | 包括阳离子聚合物的井筒应用液及其使用方法 |
| CN104806230B (zh) * | 2015-02-16 | 2018-01-05 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
| CN108316251A (zh) * | 2018-02-06 | 2018-07-24 | 浙江大学 | 一种采用超高韧性水泥基复合材料加固堆石坝混凝土面板的方法 |
| CN110159239A (zh) * | 2019-05-06 | 2019-08-23 | 中国石油天然气股份有限公司 | 一种直井大规模水力压裂油套同注压裂方法 |
Non-Patent Citations (4)
| Title |
|---|
| "An optimal design of network-fracture acidification for ultra-deep gas wells in the Lower Permian strata of the western Sichuan Basin";Guo Jianchun et al.;《Natural Gas Industry B》;20171102;第4卷(第6期);415-422 * |
| "Wellbore temperature distribution during circulation stage when well-kick occurs in a continuous formation from the bottom-hole";Zheng Zhang et al.;《Energy》;20180910;第164卷;964-977 * |
| "复杂碳酸盐岩储层水平井高效酸化数值模拟研究";薛衡;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20180515(第05期);B019-7 * |
| "非均质碳酸盐岩转向酸化模型研究";廖毅;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20171115(第11期);B019-130 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110929447A (zh) | 2020-03-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111400893B (zh) | 一种套管式地埋管换热器流体温度场分析方法 | |
| CN101560880B (zh) | 超临界井筒多相流动实验装置 | |
| CN107045671A (zh) | 产水气井积液风险预测方法 | |
| Zhang et al. | Analysis of wellbore temperature distribution and influencing factors during drilling horizontal wells | |
| CN104832131B (zh) | 一种基于机采和集油系统整体的集油参数优选方法 | |
| CN108827830B (zh) | 一种高温高压钻井液流动性能测试装置及方法 | |
| CN104807850B (zh) | 一种测量油气井井筒流体、油井管热力学参数的实验装置及方法 | |
| CN110929447B (zh) | 一种稠化酸酸化过程中井筒温度场数值计算方法 | |
| CN106777709A (zh) | 基于地面实测电参数的潜油电泵井动液面计算的方法 | |
| CN109948182A (zh) | 一种针对中深层地热井井间距的计算方法 | |
| CN104389589B (zh) | 基于空心杆确定井筒温度场分布的方法以及系统 | |
| Lin et al. | Numerical model for geothermal energy utilization from double pipe heat exchanger in abandoned oil wells | |
| CN103439356A (zh) | 用于超临界二氧化碳定压比热测量的实验装置与方法 | |
| CN111679061A (zh) | 冻土地层融化沉降试验模拟装置及方法 | |
| CN105910951A (zh) | 油气井井筒流体中气体溶解度的测量装置及测量方法 | |
| Guo et al. | Study of wellbore heat transfer considering fluid rheological effects in deep well acidizing | |
| CN105443052A (zh) | 一种油井洗井用空心抽油杆下入深度的确定方法 | |
| CN205157527U (zh) | 钻井液性能与油气含量相关性实验装置 | |
| CN204536240U (zh) | 一种测量油气井井筒流体、油井管热力学参数的实验装置 | |
| CN207660565U (zh) | 一种定向井钻井液携岩效率评价实验装置 | |
| Nie | Circulating preheating model of full-length horizontal wellbore in heavy oil reservoirs with multiple thermal fluid injection | |
| CN116629057A (zh) | 海上油气田高含co2原油井筒流动温度分布分析方法 | |
| CN104019919A (zh) | 一种驱替实验过程中测量岩心内部温度的装置及方法 | |
| CN212255154U (zh) | 一种低温地热套管的检测装置 | |
| CN204832945U (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |