一种注气井井筒压力分布的确定方法技术

本发明专利技术公开了一种注气井井筒压力分布的确定方法。该确定方法包括以下步骤：获取注气井井筒的测试参数；取注气井井筒的井口为坐标原点，竖直向下为正，按照预定计算步长将所述注气井井筒划分为多个微元段；根据所述测试参数，以注气井井筒井口处为初始计算点，依次计算各微元段出口处的压力值和温度值；根据所述各微元段出口处的压力值确定所述注气井井筒压力分布。该方法根据气体稳定流动的运动方程，在平均参数法基础上，考虑了动能变化时计算井筒压力分布的数学模型。该方法具有较好的稳定性和收敛性，适合计算机编程。适合计算机编程。适合计算机编程。

【技术实现步骤摘要】
一种注气井井筒压力分布的确定方法


[0001]本专利技术涉及天然气注采领域，具体涉及一种注气井井筒压力分布的确定方法。

技术介绍

[0002]注气井井筒压力分布是储气库建井、油管设计、井口选择的重要参数，是注气井生产完井管柱动态分析核心参数，因此研究计算注气井井筒内的压力分布十分必要。
[0003]计算气井单相流体井筒压力的常规算法均采用能量方程，考虑气体稳态流动时的传热稳定，在能量守恒的前提下导出气体流动方程的积分表达式，目前常用的是采用Moody、AGA与Cullender和Smith方法等迭代法求解。
[0004]以上许多研究者在井筒压力计算的各个方面已经做出了大量的研究，并且给出了计算气井井筒压力分布的方程式。但是在计算过程中，需要多次迭代，计算结果针对较高注气量的储气库注气井来说误差仍然是无法接受的；因此本专利技术结合大量现场实际测试资料，采用统计试验中的因子设计分析方法，回归了一种可以匹配注气量在80
×
104m3/d以内的注气井的稳定的井筒内压力数据的代数公式，该公式可准确预测注气井井筒内流动压力。

技术实现思路

[0005]本专利技术的目的在于提供一种注气井井筒压力分布的确定方法，该确定方法针对注气井在注入气体、且较高注气量(20-80万方/d)情况下的井筒压力分布的精确计算，是一种计算高压(21-35MPa)注气井井筒压力分布的快速、简便的方法，能够精确、快速计算任意流动状况下气体沿井筒内壁的压力分布。
[0006]为了实现以上目的，本专利技术采用以下技术方案：
[0007]本专利技术提供一种注气井井筒压力分布的确定方法，该确定方法包括以下步骤：
[0008]S100、获取注气井井筒的测试参数，
[0009]S200、取注气井井筒的井口为坐标原点，竖直向下为正，按照预定计算步长将所述注气井井筒划分为多个微元段；
[0010]S300、根据所述测试参数，以注气井井筒井口处为初始计算点，依次计算各微元段出口处的压力值和温度值；
[0011]S400、根据所述各微元段出口处的压力值确定所述注气井井筒压力分布。
[0012]根据本专利技术的确定方法，优选的，所述测试参数包括油管内径D、井口油压P0、气井深度H、井口温度T0、井底温度Te、地温梯度g
T
、注气量Q
sc
、天然气与空气的相对密度γ。
[0013]根据本专利技术的确定方法，优选的，S300的过程包括：
[0014]S301、以注气井井筒井口处为初始计算点，任意设定步长dl，l＝0，i＝1，
[0015]S302、令P
′1＝P0；T
′1＝T0；其中P0为井口油压，T0为井口温度，P
i
’
为第i个微元段的入口压力，P
i+1
’
为第i个微元段的出口压力，即第i+1个微元段的入口压力，T
i
’
为第i个微元段的入口温度，T
i+1
’
为第i个微元段的出口温度，即第i+1个微元段的入口温度；
[0016]S303、根据能量守恒和地层温度变化特性得到该微元段的出口压力P2’
和出口温度T2’
；
[0017]S304、计算该微元段的平均压力和平均温度：
[0018]S305、根据该微元段的入口压力和入口温度得到关于压缩因子、温度和压力的简化变量I1；根据该微元段的出口压力和出口温度得到关于压缩因子、温度和压力的简化变量I2；
[0019]S306、将I1和I2代入井筒压力计算模型式中计算得到dl段出口处的井筒压力P2；
[0020]S307、计算压力差dP＝P
2-P1；计算天然气在P1和P2压力条件下的真实流速v1和v2，dv＝dv
2-dv1；
[0021]S308、当|T
2-T2′
|≤0.001，|P
2-P2′
|≤0.001时，结束循环；否则，令T2′
＝T2、P
′2＝P2返回S304再次计算；
[0022]S309、i++，令P
′1＝P2，T
′1＝T2，l＝l+dl，返回S303继续迭代计算下一微元段的温度、压力；当l≤H时则迭代结束。
[0023]根据本专利技术的确定方法，优选的，S303中，所述该微元段的出口压力和出口温度为：
[0024][0025]其中P0为井口油压，T0为井口温度，P
i
’
为第i个微元段的入口压力，P
i+1
’
为第i个微元段的出口压力，即第i+1个微元段的入口压力，T
i
’
为第i个微元段的入口温度，T
i+1
’
为第i个微元段的出口温度，即第i+1个微元段的入口温度；g
T
为地温梯度。
[0026]根据本专利技术的确定方法，优选的，S305中所述关于压缩因子、温度和压力的简化变量I1和I2分别通过公式(5)计算得到：
[0027][0028]其中，P为压力，T为温度，Z为压缩因子。
[0029]根据本专利技术的确定方法，优选的，S306中所述井筒压力计算模型为式4)：
[0030][0031]P1为dl段入口处的井筒压力，MPa；P2为dl段出口处的井筒压力，MPa；θ为dl井段与水平方向的夹角；A1至A4为压力模型的回归系数；Q
sc
为注气量，m3/d；D为油管内径，m；g为重
力加速度，9.81m/s2；γ为天然气与空气的相对密度；f为天然气的摩擦阻力系数。
[0032]本专利技术推导出储气库高压注气井在单相流条件下，在井口流压和温度已知的条件下，沿井深分布的压力计算方程，通过该方程可快速计算注天然气井的井筒压力分布。而且本专利技术所涉及的确定方法具有良好的计算精度，计算速度较快。
[0033]根据本专利技术的确定方法，优选的，提高压力精度的回归系数A1至A4分别取值为：2.295、0.869、9.17
×
10-3
和1.853。
[0034]根据本专利技术的确定方法，优选的，S307中所述计算天然气在P1和P2压力条件下的真实流速v1和v2，通过下式(2)计算得到：
[0035]天然气在一定压力条件下的真实流速：
[0036][0037]式中，Q
sc
为天然气在地面标准状态下的真实流速(即为注气量，在井口可以直接读取)，m3/d；Z为天然气的压缩因子；P为所对应的一定压力，MPa；T为温度，K。
[0038]根据本专利技术的确定方法，优选的，该确定方法适用的注气量为20-80万方/d。
[0039]根据本专利技术的确定方法，优选的，该确定方法包括以下假设条件：
[0040]流体流动状态为稳定单向流动，流体为天然气单相流；
[0041]井筒内传热为稳定本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种注气井井筒压力分布的确定方法，其特征在于，该确定方法包括以下步骤：S100、获取注气井井筒的测试参数；S200、取注气井井筒的井口为坐标原点，竖直向下为正，按照预定计算步长将所述注气井井筒划分为多个微元段；S300、根据所述测试参数，以注气井井筒井口处为初始计算点，依次计算各微元段出口处的压力值和温度值；S400、根据所述各微元段出口处的压力值确定所述注气井井筒压力分布。2.根据权利要求1所述注气井井筒压力分布的确定方法，其特征在于，所述测试参数包括油管内径、井口油压、气井深度、井口温度、井底温度、地温梯度、注气量、天然气与空气的相对密度。3.根据权利要求2所述注气井井筒压力分布的确定方法，其特征在于，S300的过程包括：S301、以注气井井筒井口处为初始计算点，任意设定步长dl，l＝0，i＝1；S302、令P1′
＝P0；T1′
＝T0；其中P0为井口油压，T0为井口温度，P
i
’
为第i个微元段的入口压力，P
i+1
’
为第i个微元段的出口压力，即第i+1个微元段的入口压力，T
i
’
为第i个微元段的入口温度，T
i+1
’
为第i个微元段的出口温度，即第i+1个微元段的入口温度；S303、根据能量守恒和地层温度变化特性得到该微元段的出口压力P2’
和出口温度T2’
；S304、计算该微元段的平均压力和平均温度：S305、根据该微元段的入口压力和入口温度得到关于压缩因子、温度和压力的简化变量I1；根据该微元段的出口压力和出口温度得到关于压缩因子、温度和压力的简化变量I2；S306、将I1和I2代入井筒压力计算模型式中计算得到dl段出口处的井筒压力P2；S307、计算压力差dP＝P
2-P1；计算天然气在P1和P2压力条件下的真实流速v1和v2，dv＝dv
2-dv1；S308、当|T
2-T2′
|≤0.001，|P
2-P2′
|≤0.001时，结束循环；否则，令T2′
＝T2、P2′
＝P2返回S304再次计算；S309、i++，令P1′
＝P2，T1′
＝T2，l＝l+dl，返回S303继续迭代计算下一微元段的温度...

【专利技术属性】
技术研发人员：贺梦琦，王玲，陈小凯，秦辉，梁策，张绍辉，孙振宇，宋阳，刘鹍澎，贾纯真，栾睿智，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：