CN115831243B - 一种结垢临界温度压力预测方法 - Google Patents
一种结垢临界温度压力预测方法 Download PDFInfo
- Publication number
- CN115831243B CN115831243B CN202310140970.0A CN202310140970A CN115831243B CN 115831243 B CN115831243 B CN 115831243B CN 202310140970 A CN202310140970 A CN 202310140970A CN 115831243 B CN115831243 B CN 115831243B
- Authority
- CN
- China
- Prior art keywords
- dimensionless
- pressure
- scaling
- model
- temperature
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 150000002500 ions Chemical class 0.000 claims abstract description 100
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 60
- 230000000694 effects Effects 0.000 claims abstract description 43
- 239000000126 substance Substances 0.000 claims abstract description 41
- 239000008151 electrolyte solution Substances 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000013178 mathematical model Methods 0.000 claims abstract description 17
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims abstract description 16
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000000903 blocking effect Effects 0.000 claims abstract description 5
- 239000000243 solution Substances 0.000 claims description 60
- 239000012071 phase Substances 0.000 claims description 47
- 239000007789 gas Substances 0.000 claims description 42
- 239000007788 liquid Substances 0.000 claims description 31
- 239000003792 electrolyte Substances 0.000 claims description 27
- 230000015572 biosynthetic process Effects 0.000 claims description 24
- 239000012530 fluid Substances 0.000 claims description 23
- 230000014509 gene expression Effects 0.000 claims description 23
- 239000007791 liquid phase Substances 0.000 claims description 21
- 239000000203 mixture Substances 0.000 claims description 19
- 230000009471 action Effects 0.000 claims description 17
- 239000002245 particle Substances 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 14
- 238000004090 dissolution Methods 0.000 claims description 12
- 230000008020 evaporation Effects 0.000 claims description 11
- 238000001704 evaporation Methods 0.000 claims description 11
- 150000001875 compounds Chemical class 0.000 claims description 10
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 239000011575 calcium Substances 0.000 claims description 5
- 229910017053 inorganic salt Inorganic materials 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 239000007790 solid phase Substances 0.000 claims description 5
- 150000001450 anions Chemical class 0.000 claims description 4
- 229910000019 calcium carbonate Inorganic materials 0.000 claims description 4
- 239000008398 formation water Substances 0.000 claims description 3
- 239000003345 natural gas Substances 0.000 claims description 3
- 235000013619 trace mineral Nutrition 0.000 claims description 2
- 239000011573 trace mineral Substances 0.000 claims description 2
- 230000009467 reduction Effects 0.000 abstract description 6
- 229940021013 electrolyte solution Drugs 0.000 description 22
- 239000000047 product Substances 0.000 description 20
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000005684 electric field Effects 0.000 description 4
- 150000001768 cations Chemical class 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000003595 mist Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000004836 empirical method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000002332 oil field water Substances 0.000 description 2
- 239000002244 precipitate Substances 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- RECVMTHOQWMYFX-UHFFFAOYSA-N oxygen(1+) dihydride Chemical compound [OH2+] RECVMTHOQWMYFX-UHFFFAOYSA-N 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000009491 slugging Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Landscapes
- Physical Or Chemical Processes And Apparatus (AREA)
Abstract
本发明公开了一种结垢临界温度压力预测方法,包括以下步骤:S1、构建井筒结垢堵塞物理模型和结垢热力学平衡数学模型;S2、构建水蒸气含量相态计算模型同时构建混合电解质溶液中离子活度系数模型;S3、构建硫酸盐与碳酸盐结垢预测数学模型;S4、构建井筒压力温度梯度模型,通过迭代循环求解模型参数,得到井筒温度压力关系曲线;S5、求解结垢预测模型,获得结垢生成线;S6、根据现场资料进行实例分析。本发明考虑高温下压力降低使得水蒸发与溶解在水中的CO2逸出的因素,基于相平衡、化学平衡、活度系数及溶度积规则,建立了结垢预测模型,获得了结垢临界温度压力判别图版。
Description
技术领域
本发明涉及油气开发领域,具体涉及一种结垢临界温度压力预测方法。
背景技术
随着油田开发逐渐趋于成熟,结垢问题已经成为很多油田在生产运行中面临的主要问题。地层结垢则会引起储层渗透率降低,井筒或管线结垢会导致管柱有效内径减小,流体输送效率降低,开发成本增加,严重影响油田的正常生产。
目前,国内外已有众多学者对结垢预测模型进行研究,Langelier(1936)提出的饱和指数法未考虑压力、含盐量及其他碱度等因素的影响,只能表示某种水中的产生碳酸钙垢的趋势,不能计算结垢量;Ryzner(1944)提出的稳定指数法能够较为准确的判断矿化度较高和pH值较高地方的碳酸钙结垢趋势,但只考虑了CaCO3这一种结垢物质的溶解平衡;Stiff-Davis(1952)基于反应过程中的热力学平衡,提出了半经验性的饱和指数法,对部分特定油藏预测结果与实际偏差较大;Skillman(1969)提出了热力学溶解度法,基于热力学理论以及溶解平衡原理建立预测模型,可用于粗略计算油田水CaSO4结垢趋势,未考虑化学平衡;Oddo JE(1982)考虑了不同的温度压力条件下的结垢Is指数,未考虑化学动力学因素;Oddo-Tomosn(1994)考虑了热力学平衡原理、溶解平衡原理二氧化碳逸度和弱酸等因素,属于半经验性预测模型;罗明良(2003)基于经典溶液理论、离子缔合理论、饱和指数以及对实验数据进行多元非线性回归,建立了基于硫酸盐沉淀的结垢趋势预测模型;李丽(2012)基于物料守衡、相态变化、电解质溶液中的化学平衡建立了适用于实际气藏流体中气-液(液烃)-液(电解质溶液)-固(无机盐垢)多相平衡的无机盐结垢。
但是,目前的结垢趋势预测模型还存在以下问题:(1)未考虑高温下压力降低使得水蒸发加剧,引起成垢离子含量增加,同时未考虑溶解在水中的CO2逸出,其分压降低,容易形成CaCO3垢;(2)目前的结垢趋势预测模型都是属于半经验性的油田水饱和指数方程,未综合考虑相平衡、化学平衡、离子缔合理论、活度系数及溶度积规则;(3)目前的结垢预测模型都不能得到结垢生成线,不能获取结垢临界温度压力判别图版。
发明内容
本发明的目的在于提供一种结垢临界温度压力预测方法,考虑高温下压力降低使得水蒸发与溶解在水中的CO2逸出的因素,基于相平衡、化学平衡、活度系数及溶度积规则,建立了结垢预测模型,获得了结垢临界温度压力判别图版。
为实现上述目的,本发明提供了以下技术方案:
本发明提供的一种结垢临界温度压力预测方法,包括以下步骤:
S1、构建井筒结垢堵塞物理模型和结垢热力学平衡数学模型;
S2、构建水蒸气含量相态计算模型,考虑水蒸发对离子浓度的影响,同时构建混合电解质溶液中离子活度系数模型;
S3、构建硫酸盐与碳酸盐结垢预测数学模型,考虑高温下节流降压导致CO2逸出对碳酸盐结垢的影响;
S4、构建井筒压力温度梯度模型,通过迭代循环求解模型参数,得到井筒温度压力关系曲线;
S5、满足相平衡、化学平衡以及溶解平衡时输出饱和指数,求解结垢预测模型,获得结垢生成线;
S6、根据现场资料进行实例分析,结合结垢生成线与井筒温度压力曲线关系得到结垢临界温度压力判别图版。
进一步的,步骤S1中,所述井筒结垢堵塞物理模型的假设条件为:
(1)井筒分为长度为dL的等距微元体,每个微元体中垢物类型主要为碳酸盐与硫酸盐;
(2)不考虑动力学因素,高压气井井筒中流体的无机盐结垢是综合考虑各组分间相平衡、化学平衡以及溶解平衡的复杂过程;
(3)平衡条件下,天然气、地层水、电解质溶液应满足物料守恒与热力学平衡。
进一步的,步骤S1中,所述结垢热力学平衡数学模型满足液-固相平衡理论、原子个数守恒定律和电荷守恒定律,其中,
液-固相平衡理论的表达式为
原子个数守恒,实质上就是物质平衡,其表达式为
进一步的,步骤S2中,所述水蒸气含量相态计算模型的表达式如下:
基于相平衡原理,不同压力温度下的水蒸气含量可通过闪蒸计算确定,根据稳定状态吉布斯自由能最小的原则,气液闪蒸计算的相态判定准则为:
当体系处于气液相平衡状态时,由物料守恒可得:
式中:xi为液相中i组分的摩尔分数,无量纲;yi为气相中i组分的摩尔分数,无量纲;由上式可知:
进一步的,步骤S2中,所述混合电解质溶液中离子活度系数模型的表达式如下:
式中:为长程作用项;与为二粒子作用项;和为三粒子作用项;是正离子的电荷数,无量纲;ma为负离子浓度,mol/kg;mi为离子浓度,mol/kg;zi为电荷数,无量纲;mc为正离子浓度,mol/kg;、分别为正离子与负离子的电荷数,无量纲;是电解质电离出的正离子的物质的量,mol;是电解质电离出的负离子的物质的量,mol;是电解质电离出总的离子物质的量,mol;c表示正离子个数,无量纲;a表示负离子个数,无量纲。
进一步的,步骤S3中,所述硫酸盐结垢预测数学模型的表达式如下:
进一步的,步骤S3中,所述碳酸盐结垢预测数学模型的表达式如下:
溶液中CaCO3垢的饱和指数可表示为:
式中:为碳酸钙的饱和指数,无量纲;为Ca2+的活度,mol/kg;为的活度,mol/kg;K2为的溶度积常数,无量纲;为CO2分压,bar;为CO2的逸度系数,无量纲;Ksp为CaCO3的溶度积常数,无量纲;K1为H2CO3的溶度积常数,无量纲;为CO2的亨利常数,无量纲;
溶液中结垢趋势判别方法:当饱和指数等于0时,溶液饱和无结垢;当饱和指数大于1时,溶液过饱和有结垢;当饱和指数小于1时,溶液未饱和无结垢。
进一步的,步骤S4中,所述井筒压力温度梯度模型包括井筒压力梯度模型和井筒温度梯度模型;
井筒压力梯度模型的表达式为:
式中:P为压力,Pa;z为井筒单位微元体长度,m;ρm为气液混合物密度,kg/m3;g为重力加速度,m/s2;θ为倾斜角,度;fm为混合物摩阻系数,无量纲;vm为气液混合物流速,m/s;D为油管内径,m;vsg为气体表观流速,m/s。
进一步的,所述井筒温度梯度模型的表达式为:
式中:为井筒温度,K;为气体的比热,J/(kg·K);Z为偏差因子,无量纲;T为体系温度,K;P为压力,Pa;g为重力加速度,m/s2;v为流体流速,m/s;D为油管内径,m;为气体密度,kg/m3;为流体质量流量,kg/s;q为消耗热量,J/(m·s);z为井筒单位微元体长度,m;为流体的比热,J/(kg·K)。
基于上述技术方案,本发明实施例至少可以产生如下技术效果:
本发明提供的结垢临界温度压力预测方法,在结垢预测模型中,温度和压力是影响结垢趋势的两个关键因素。油田现场可通过测试得到井口压力与地层温度,利用井筒压力温度梯度模型计算得到全井筒的温度与压力分布曲线。此外,结垢预测模型中还考虑了不同温度压力下水的相态变化与CO2在水中的溶解与逸出,使得模型更符合油田现场的结垢环境,同时提出了结垢临界温度压力生成线。为方便油田现场使用结垢预测模型,更加直观的判断地层或井筒的结垢情况,将地层温度压力、井筒温度压力关系曲线与结垢生成线组合得到结垢判别图版,使得现场人员可直观分析地层或井筒的结垢情况。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1是本发明井筒球座附近结垢堵塞机理示意图;
图2是本发明水蒸气对结垢的影响示意图;
图3是本发明结垢临界温度压力生成线求解流程图;
图4是本发明A井CaSO4饱和指数与压力的关系曲线;
图5是本发明A井CaSO4饱和指数与温度的关系曲线;
图6是本发明A井CaSO4垢生成线与井筒温度压力线对比;
图7是本发明 A井CaCO3饱和指数与压力的关系曲线;
图8是本发明A井CaCO3饱和指数与温度的关系曲线;
图9是本发明A井CaCO3垢生成线与井筒温度压力线对比;
图10是本发明A井结垢判别图版对应温度压力。
具体实施方式
本发明提供一种结垢临界温度压力预测方法,考虑高温下压力降低使得水蒸发与溶解在水中的CO2逸出的因素,基于相平衡、化学平衡、活度系数及溶度积规则,建立了结垢预测模型,获得了结垢临界温度压力判别图版。本发明采用下述的技术方案:
S1、构建井筒结垢堵塞物理模型和结垢热力学平衡数学模型
1.井筒球座附近结垢堵塞机理示意图如图1所示,结垢类型为碳酸盐与硫酸盐结垢。
井筒中结垢堵塞物理模型的假设条件为:
(1)井筒分为长度为dL的等距微元体,每个微元体中垢物类型主要为碳酸盐与硫酸盐;
(2)不考虑动力学因素,高压气井井筒中流体的无机盐结垢是综合考虑各组分间相平衡、化学平衡以及溶解平衡的复杂过程;
(3)平衡条件下,天然气、地层水、电解质溶液应满足物料守恒与热力学平衡。
2.构建结垢热力学平衡数学模型
随着油田现场生产不断进行,流体在井筒内不同位置处的温度压力不同,因此化合物的沉淀溶解过程是一个动态平衡,满足液-固相平衡理论,其表达式为
对于电解质溶液而言,溶液中存在着不同离子的电离平衡,电离平衡过程满足原子个数守恒,实质上就是物质平衡,其表达式为
在溶液的化学平衡过程中,电荷守恒也必须满足,在式(2)中表现为。同时溶液中气液相平衡也必须满足,在井筒中主要体现为不同温度压力下水蒸气与液态水的相互转化以及CO2在溶液中的溶解与逸出,水蒸气与CO2在溶液中的相平衡将在下文中详细阐述。
S2、构建水蒸气含量相态计算模型,考虑水蒸发对离子浓度的影响,同时构建混合电解质溶液中离子活度系数模型
1.水蒸气含量相态计算模型
井筒内水的存在状态为水蒸气和液态水,不同井深处的水蒸气含量不一样,不同井深处的液态水量也不一样,液态水量将影响离子浓度、离子强度以及离子活度,进一步影响结垢趋势与结垢量。因此,水蒸气的含量计算将影响结垢量的多少,水蒸气对结垢的影响如图2所示。
井底附近的高温使水蒸发较快,气组分中的水蒸气含量较高,水中离子浓度增加会加剧垢的形成;井口附近的温度较低,水蒸发较慢,气组分中的水蒸气含量较低,且井底的高温水蒸气在井口附近的低温条件下冷凝为液态水。基于相平衡原理,不同压力温度下的水蒸气含量可通过闪蒸计算确定。
根据稳定状态吉布斯自由能最小的原则,气液闪蒸计算的相态判定准则为:
采用威尔逊(Wilson)公式来估算平衡常数的初值,其表达式为
当体系处于气液相平衡状态时,由物料守恒可得通过式(6)得到xi和yi的表达式为
式中:xi为液相中i组分的摩尔分数,无量纲;yi为气相中i组分的摩尔分数,无量纲。
通过式(6)得到xi和yi的表达式为
由上式可知
2.混合电解质溶液活度系数模型
对于实际井筒中的流体可以看作一种结垢反应溶液体系,结垢反应的关键是计算溶液中的各种垢的溶解度,溶解度的大小则与化合物及离子的活度系数直接相关。活度系数的计算有经验方法与热力学方法,其中经验法由于误差太大且没有严密的理论基础,因此目前普遍采用热力学方法。
Debye-Hückel方程是最早提出的热力学活度系数模型,但其适用离子强度较低的电解质溶液,对于高离子强度的电解质溶液计算误差较大。Pitzer则在Debye-Hückel方程基础上,考虑了长程作用相、水的介电常数、二粒子作用项及三粒子作用项,其中长程作用相为
式中:NA为阿伏伽德罗常数,mol-1;为水的密度,kg/m3;e是元电荷,C;D为水的介电常数,无量纲;k是Boltzmann常数,1.3806×10-23J/K;T是温度,K;是离子最近距离,m;I是离子强度,mol/kg;mi为离子浓度,mol/kg;zi为电荷数,无量纲,不同温度下的值如下表1所示。
水在外加电场时会产生感应电荷而削弱电场,外加电场与介质中电场的比值即为介电常数。水的介电常数计算相关参数见表2所示,其表达式为
式中:D为水的介电常数,无量纲;U1~U9是计算介电常数参数,无量纲;P是压力,Pa;T是温度,K。
表2 水的介电常数计算相关参数
对于电解质溶液中的单个离子,Pitzer活度系数模型中的二粒子作用项与三粒子作用项分别为
式中:、、为二粒子作用项相关参数,kg·mol-1;为1.4;I是离子强度,mol/kg;e是元电荷,C;是参数,与电解质在水中的缔合程度有关,无量纲;为12;是正离子的电荷数,无量纲;是负离子的电荷数,无量纲;是电解质电离出的正离子的物质的量,mol;是电解质电离出的负离子的物质的量,mol;、是三粒子作用系数,无量纲。
式中:、是二粒子作用项相关参数,kg·mol-1;是正离子的电荷数,无量纲;是负离子的电荷数,无量纲;是电解质电离出的正离子的物质的量,mol;是电解质电离出的负离子的物质的量,mol;、是三粒子作用系数,无量纲。
式中:为长程作用项;与为二粒子作用项;和为三粒子作用项;是正离子的电荷数,无量纲;m a为负离子浓度,mol/kg;m i为离子浓度,mol/kg;z i为电荷数,无量纲;m c为负离子浓度,mol/kg;、分别为正离子与负离子的电荷数,无量纲;是电解质电离出的正离子的物质的量,mol;是电解质电离出的负离子的物质的量,mol ;是电解质电离出总的离子物质的量,mol。
S3、构建硫酸盐与碳酸盐结垢预测数学模型,考虑高温下节流降压导致CO2逸出对碳酸盐结垢的影响
1.硫酸盐结垢预测数学模型
在一定的温度压力下,难溶电解质在溶液中有如下化学平衡
溶度积规则实际上是化学平衡移动规则在多相离子平衡中的应用,K sp即为成垢沉淀物质A m B n的热力学溶度积,其表达式为
溶度积常数与温度压力相关,与温度相关的表达式为
表3 溶度积常数相关参数值
溶度积常数与压力相关的表达式为
式中:Ksp为与压力相关的热力学溶度积,无量纲;R i 表示、、,无量纲;是参考压力下,平衡时成垢阴、阳离子的浓度乘积,无量纲;为化合物的平均活度系数,无量纲;P是压力,Pa;是参考压力,取1bar;I是离子强度,mol/kg;T是温度,K;A、B、C、D为常数,无量纲。
溶液中的离子能否形成沉淀与饱和指数直接相关,可用来判断结垢趋势,其表达式为
2.碳酸盐结垢预测数学模型
在油田现场,地层或井筒中常出现碳酸盐垢,经过油田现场垢物成分检验为CaCO3垢。在地层或井筒中的高温条件下,压力降低使得水蒸发与CO2逸出,会加剧碳酸盐垢的形成。
电解质溶液中CaCO3的化学分解反应为
则CaCO3溶度积可表示为
式中:为Ca2+的活度,mol/kg;为的活度,mol/kg;为CaCO3的活度,mol/kg;为溶液中的标准化学势,kJ/mol;为溶液中Ca2+的标准化学势,kJ/mol;为CaCO3的标准化学势,kJ/mol。
电解质溶液中H2CO3的化学分解反应为
根据实验数据回归得到关联了温度与压力的溶度积常数,其表达式为
式中:K为溶度积常数,无量纲;P为体系压力,bar;T为温度,K;Ps表示不同条件下的压力,T<373.15K时,P s=1bar,T>373.15K时,P s取水的饱和蒸汽压,bar;a 1~a 11为计算溶度积常数所需参数,其值见下表4所示;K1为H2CO3的溶度积常数,无量纲;K2为的溶度积常数,无量纲;
表4 电离平衡常数的相关参数
油气井在生产过程中,井筒中的气体组分含有CO2,随着温度压力的不断变化,CO2气体会在水中不断的溶解与逸出,当井流物处于井筒下部球座位置时,球座缩径会产生节流降压效应,溶解在水中的CO2逸出,其分压降低,易形成CaCO3垢。电解质溶液中CO2的溶解与逸出为
式中:CO2(g)为气相中的CO2,无量纲;CO2(aq)为溶液中的CO2,无量纲。
CO2在溶液体系中的浓度为
式中:为CO2在溶液体系中的质量摩尔浓度,mol/kg;为CO2的摩尔分数,无量纲;为CO2的逸度系数,无量纲;为CO2在溶液中标准化学势,无量纲;R为气体常数,取8.314J/mol/K;T为温度,K;、是分别为CO2与溶液中正、负离子的交互作用系数,无量纲;是为溶液中CO2-正离子-负离子的交互作用系数,无量纲。
电解质溶液中的CO2活度系数为
电解质溶液中的CO2活度为
电解质溶液达到平衡时,CO2的亨利常数可表示为
CaCO3垢的结垢机理方程为
溶液中CaCO3垢的饱和指数可表示为
式中:为饱和指数,无量纲;为Ca2+的活度,mol/kg;为的活度,mol/kg;K2为的溶度积常数,无量纲;为CO2分压,bar;为CO2的逸度系数,无量纲;Ksp为CaCO3的溶度积常数,无量纲;K1为H2CO3的溶度积常数,无量纲;为CO2的亨利常数,无量纲。
溶液中结垢趋势判别方法:当饱和指数等于0时,溶液饱和无结垢;当饱和指数大于1时,溶液过饱和有结垢;当饱和指数小于1时,溶液未饱和无结垢。
S4、构建井筒压力温度梯度模型,通过迭代循环求解模型参数,得到井筒温度压力关系曲线
Mukherjee-Brill在Beggs-Brill基础上提出了M-B模型,其压力梯度方程为
式中:P为压力,Pa;z为井筒单位微元体长度,m;ρ m为气液混合物密度,kg/m3;g为重力加速度,m/s2;θ为倾斜角,度;f m为混合物摩阻系数,无量纲;v m为气液混合物流速,m/s;D为油管内径,m;v sg为气体表观流速,m/s。
M-B模型持液率为
式中:H l 为持液率,无量纲;c 1、c 2、c 3、c 4、c 5、c 6为持液率公式回归参数,无量纲;θ为倾斜角,度;N l为液相粘度数,无量纲;N gv为气相速度数,无量纲;N lv为液相速度数,无量纲。
M-B模型区分泡状流、段塞流与环雾流,其判别式为:
如果N gv≥N gvsm,则为环雾流,否则为泡状流或者段塞流。
泡状流或段塞流:两相摩阻系数为无滑脱摩阻系数,由Jain公式得出。
无滑脱雷诺数为
环雾流:摩阻系数考虑相对持液率和无滑脱摩阻系数的函数,相对持液率为
混合物摩阻系数为
H R与f R的关系见下表5。
表5 HR与fR的关系
对于气体,对外做功=0,则总的压降梯度为
重力压降为
式中:dL为单位井筒段长度,m。
摩擦压降为
其中
气体节流压降(动能压降)dP u为
对于液体的动能变化为
2.井筒温度梯度模型
根据能量守恒原理,井筒中流体能量平衡方程为
式中:为流体质量流量,kg/s;为流体混合物的比焓,J/kg;v为流体流速,m/s; g为重力加速度,m/s2;为生成焓,J/mol;M h为平均摩尔质量,kg/mol;q为消耗热量,J/(m·s)。
比焓梯度用温度和压力梯度表示
井筒内温度梯度表达式为
井筒流体到地层中传热方程为
式中:T ei为地层温度,K;A′为松弛距离。
松弛距离可表示为
式中:U to为地层综合传热系数。
U to可表示为
式中:为油管外半径,m;为油管内半径,m;为井眼半径,m;为套管外半径,m;为套管内半径,m;为油管导热系数,W/(m·K);为套管导热系数,W/(m·K);为水泥环导热系数,W/(m·K);h c为环空流体热对流传热系数,W/(m·K);h r为环空流体辐射系数,W/(m·K)。
S5、满足相平衡、化学平衡以及溶解平衡时输出饱和指数,求解结垢预测模型,获得结垢生成线
由结垢预测模型可知,饱和指数是温度压力的函数,即Is=f(T,P)。当饱和指数Is=0时,对应的温度压力为结垢临界温度压力生成线,即结垢生成线,其计算框如图3所示。
S6、根据现场资料进行实例分析,结合结垢生成线与井筒温度压力关系曲线得到结垢临界温度压力判别图版
根据A井现场资料可得到模型预测的基础数据为:日产气量65×104m3/d,日产水量7.5t/d,井口压力50MPa,井底温度为135℃。A井水离子分析结果见下表所示。
表6 A井取样水离子分析结果
由结垢趋势预测模型可知,饱和指数是井筒温度压力的函数,即I s=f(T,P)。以A井为例,CaSO4与CaCO3的饱和指数与井筒压力温度的关系如图4、图5、图7、图8所示。当饱和指数F s=0时,对应的压力温度为生成垢的临界压力温度,即可得到不同压力温度下结垢判别图版如图6、图9所示。由此可知,在结垢生成线右下区域的温度压力满足I s>0,为垢生成区域;在结垢生成线左上区域的温度压力满足I s<0,不会生成垢。
由此可知,本模型预测在108.85℃以上结垢,108.85℃以下不结垢,与现场实际井筒结垢情况吻合。通过本模型获得A井的结垢判别图版如图10所示,其对应温度压力见表7所示。将A井的结垢判别图版分为三个区:地层高压不结垢区、井筒结垢区和完全不结垢区。在地层温度条件下,地层压力大于73.8MPa,地层不结垢,小于73.8MPa,地层结垢;在目前生产条件下,井筒温度大于119.8℃,井深大于4230m,井筒结垢;针对产出水成分,温度<108.85℃,压力<49.67MPa,产出水完全不结垢。
表7 A井结垢判别图版对应温度压力
上面对本发明的较佳实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (9)
1.一种结垢临界温度压力预测方法,其特征在于,包括以下步骤:
S1、构建井筒结垢堵塞物理模型和结垢热力学平衡数学模型;
S2、构建水蒸气含量相态计算模型,考虑水蒸发对离子浓度的影响,同时构建混合电解质溶液中离子活度系数模型;
S3、构建硫酸盐与碳酸盐结垢预测数学模型,考虑高温下节流降压导致CO2逸出对碳酸盐结垢的影响;
S4、构建井筒压力温度梯度模型,通过迭代循环求解模型参数,得到井筒温度压力关系曲线;
S5、满足相平衡、化学平衡以及溶解平衡时输出饱和指数,求解结垢预测模型,获得结垢生成线;
S6、根据现场资料进行实例分析,结合结垢生成线与井筒温度压力关系曲线得到结垢临界温度压力判别图版。
2.根据权利要求1所述的结垢临界温度压力预测方法,其特征在于,步骤S1中,所述井筒结垢堵塞物理模型的假设条件为:
(1)井筒分为长度为dL的等距微元体,每个微元体中垢物类型主要为碳酸盐与硫酸盐;
(2)不考虑动力学因素,高压气井井筒中流体的无机盐结垢是综合考虑各组分间相平衡、化学平衡以及溶解平衡的复杂过程;
(3)平衡条件下,天然气、地层水、电解质溶液应满足物料守恒与热力学平衡。
3.根据权利要求1所述的结垢临界温度压力预测方法,其特征在于,步骤S1中,所述结垢热力学平衡数学模型满足液-固相平衡理论、原子个数守恒定律和电荷守恒定律,其中,
液-固相平衡理论的表达式为
原子个数守恒,实质上就是物质平衡,其表达式为
4.根据权利要求1所述的结垢临界温度压力预测方法,其特征在于,步骤S2中,所述水蒸气含量相态计算模型的表达式如下:
基于相平衡原理,不同压力温度下的水蒸气含量可通过闪蒸计算确定,根据稳定状态吉布斯自由能最小的原则,气液闪蒸计算的相态判定准则为:
当体系处于气液相平衡状态时,由物料守恒可得:
式中:xi为液相中i组分的摩尔分数,无量纲;yi为气相中i组分的摩尔分数,无量纲;由上式可知:
5.根据权利要求1所述的结垢临界温度压力预测方法,其特征在于,步骤S2中,所述混合电解质溶液中离子活度系数模型的表达式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310140970.0A CN115831243B (zh) | 2023-02-21 | 2023-02-21 | 一种结垢临界温度压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310140970.0A CN115831243B (zh) | 2023-02-21 | 2023-02-21 | 一种结垢临界温度压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115831243A CN115831243A (zh) | 2023-03-21 |
CN115831243B true CN115831243B (zh) | 2023-05-02 |
Family
ID=85522002
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310140970.0A Active CN115831243B (zh) | 2023-02-21 | 2023-02-21 | 一种结垢临界温度压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115831243B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116205162B (zh) * | 2023-04-24 | 2023-07-25 | 西南石油大学 | 一种高温高压气藏结垢动态沉积堵塞预测模型构建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104343427A (zh) * | 2013-07-24 | 2015-02-11 | 中国石油化工股份有限公司 | 一种预测co2驱油藏无机结垢趋势的方法 |
CN105388083A (zh) * | 2015-10-28 | 2016-03-09 | 陕西省石油化工研究设计院 | 一种油田结垢趋势预测方法 |
CN111382486A (zh) * | 2018-12-27 | 2020-07-07 | 中国石油天然气股份有限公司 | 一种基于稳定指数预测co2气驱井筒结垢趋势的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL1016306C2 (nl) * | 2000-10-02 | 2002-08-06 | Kiwa Nv | Werkwijze voor het vroegtijdig signaleren van het optreden van scaling bij de zuivering van water. |
CN108804862A (zh) * | 2017-05-02 | 2018-11-13 | 中国石油化工股份有限公司 | 一种用于碳酸钙结垢趋势的预测方法 |
KR102479346B1 (ko) * | 2018-11-12 | 2022-12-19 | 주식회사 엘지화학 | 폴리올레핀 제조 공정에서의 파울링 예측 방법 |
CN111101925A (zh) * | 2019-11-26 | 2020-05-05 | 西南石油大学 | 一种评价注水井结垢趋势的方法 |
CN112100850B (zh) * | 2020-09-16 | 2021-04-13 | 西南交通大学 | 基于系统动力学的高浓度有机废水集输管路结垢预测方法 |
-
2023
- 2023-02-21 CN CN202310140970.0A patent/CN115831243B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104343427A (zh) * | 2013-07-24 | 2015-02-11 | 中国石油化工股份有限公司 | 一种预测co2驱油藏无机结垢趋势的方法 |
CN105388083A (zh) * | 2015-10-28 | 2016-03-09 | 陕西省石油化工研究设计院 | 一种油田结垢趋势预测方法 |
CN111382486A (zh) * | 2018-12-27 | 2020-07-07 | 中国石油天然气股份有限公司 | 一种基于稳定指数预测co2气驱井筒结垢趋势的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115831243A (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cui et al. | Geothermal energy exploitation from depleted high-temperature gas reservoirs by recycling CO2: The superiority and existing problems | |
Jia et al. | Modeling and analysis of carbonate matrix acidizing using a new two-scale continuum model | |
Shen et al. | Fouling of enhanced tubes for condensers used in cooling tower systems: A literature review | |
CN115831243B (zh) | 一种结垢临界温度压力预测方法 | |
CN106761733A (zh) | 一种稠油油藏水平井蒸汽吞吐初期产能预测方法 | |
Navarrete et al. | Laboratory and theoretical studies for acid fracture stimulation optimization | |
CA2656593C (en) | A method of detecting dryness of wet steam of a once-through steam injection boiler and a detecting apparatus therefor | |
CN104895560A (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
Liu et al. | Maximum gas production rate for salt cavern gas storages | |
Ogden et al. | Desalination of hypersaline brines via Joule-heating: Experimental investigations and comparison of results to existing models | |
CN111101925A (zh) | 一种评价注水井结垢趋势的方法 | |
Yang et al. | Prediction of CaCO3 scaling in water injection wellbore | |
Luo et al. | Experimental study of gas-liquid two-phase flow for high velocity in inclined medium size tube and verification of pressure calculation methods | |
Addison et al. | Brine silica management at mighty river power, New Zealand | |
Von Hirtz | Silica scale control in geothermal plants—Historical perspective and current technology | |
Wei et al. | Experimental research on brine crystallization mechanism in solution mining for salt cavern energy storage | |
CN115936258B (zh) | 一种井筒结垢动态沉积堵塞预测模型构建方法 | |
Lu et al. | Experimental investigation and numerical modelling of transient two-phase flow in a geysering geothermal well | |
CN116738113A (zh) | 预测地热井筒中结垢位置的数学模型的构建方法 | |
Lai et al. | Calcium carbonate scaling kinetics in oilfield gathering pipelines by using a 1D axial dispersion model | |
Cao et al. | Modeling of carbon dioxide dissolution in an injection well for geologic sequestration in aquifers | |
Francke | Thermo-hydraulic model of the two-phase flow in the brine circuit of a geothermal power plant | |
Tjelta et al. | Corrosion, scaling and material selection in deep geothermal wells-Application to IDDP-2 | |
Champel | Discrepancies in brine density databases at geothermal conditions | |
Garber et al. | Internal corrosion rate prediction in pipelines and flowlines using a computer model |
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 |