一种注气井井筒压力分布的确定方法与流程

1.本发明涉及天然气注采领域，具体涉及一种注气井井筒压力分布的确定方法。


背景技术：

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


技术实现要素：

5.本发明的目的在于提供一种注气井井筒压力分布的确定方法，该确定方法针对注气井在注入气体、且较高注气量(20-80万方/d)情况下的井筒压力分布的精确计算，是一种计算高压(21-35mpa)注气井井筒压力分布的快速、简便的方法，能够精确、快速计算任意流动状况下气体沿井筒内壁的压力分布。
6.为了实现以上目的，本发明采用以下技术方案：
7.本发明提供一种注气井井筒压力分布的确定方法，该确定方法包括以下步骤：
8.s100、获取注气井井筒的测试参数，
9.s200、取注气井井筒的井口为坐标原点，竖直向下为正，按照预定计算步长将所述注气井井筒划分为多个微元段；
10.s300、根据所述测试参数，以注气井井筒井口处为初始计算点，依次计算各微元段出口处的压力值和温度值；
11.s400、根据所述各微元段出口处的压力值确定所述注气井井筒压力分布。
12.根据本发明的确定方法，优选的，所述测试参数包括油管内径d、井口油压p0、气井深度h、井口温度t0、井底温度te、地温梯度g
t
、注气量q
sc
、天然气与空气的相对密度γ。
13.根据本发明的确定方法，优选的，s300的过程包括：
14.s301、以注气井井筒井口处为初始计算点，任意设定步长dl，l＝0，i＝1，
15.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个微元段的入口温度；
16.s303、根据能量守恒和地层温度变化特性得到该微元段的出口压力p2’
和出口温度t2’
；
17.s304、计算该微元段的平均压力和平均温度：
18.s305、根据该微元段的入口压力和入口温度得到关于压缩因子、温度和压力的简化变量i1；根据该微元段的出口压力和出口温度得到关于压缩因子、温度和压力的简化变量i2；
19.s306、将i1和i2代入井筒压力计算模型式中计算得到dl段出口处的井筒压力p2；
20.s307、计算压力差dp＝p
2-p1；计算天然气在p1和p2压力条件下的真实流速v1和v2，dv＝dv
2-dv1；
21.s308、当|t
2-t2′
|≤0.001，|p
2-p2′
|≤0.001时，结束循环；否则，令t2′
＝t2、p
′2＝p2返回s304再次计算；
22.s309、i++，令p
′1＝p2，t
′1＝t2，l＝l+dl，返回s303继续迭代计算下一微元段的温度、压力；当l≤h时则迭代结束。
23.根据本发明的确定方法，优选的，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]
井筒内传热为稳定传热；
[0042]
井筒内天然气流体符合能量守恒规律；
[0043]
地层温度按线性变化。
[0044]
本发明中的井筒压力计算模型式4)通过以下推导过程得到：
[0045]
取井口为坐标原点，竖直向下为正，在井筒上任意取一段长为dl的微元段，根据能量守恒推导出天然气的稳定流动的运动方程为：
[0046][0047]
式中，dp为dl井段井筒压力变化值，mpa；v为气体真实流速，m/s；ρ为气体的真实密度，kg/m3；g为常数，重力加速度，9.81m/s2；f为天然气的摩擦阻力系数；d为流动油管、套管的内径，m。
[0048]
天然气在一定压力条件下的真实流速：
[0049][0050]
式中，q
sc
为天然气在地面标准状态下的真实流速(井口可以直接读取注气量，根据气体状态方程可以进行换算)，m3/d；z为天然气的压缩因子；p为dl井段筒天然气的压力平均值，mpa；t为dl井段天然气的温度平均值，k。
[0051]
天然气在一定压力条件下的密度：
[0052][0053]
式中，γ为天然气气体与空气的相对密度。
[0054]
将(2)(3)代入(1)得到dl段井筒压力计算模型为：
[0055][0056]
式中，p1为dl段入口处的井筒压力，mpa；p2为dl段出口处的井筒压力，mpa；θ为dl井段与水平方向的夹角。
[0057]
其中：a1至a4为压力模型的回归系数，本发明经过大量现场实际测试资料拟合得到a1至a4的取值，如表1所示：
[0058]
表1提高压力精度的回归系数
[0059]
系数a1a2a3a4系数值2.2950.8699.17
×
10-3
1.853
[0060]
本发明通过合理的假设和必要的简化，提出了一个修正改进的方法，可为储气库井，特别是定向井、高注气量井的工程设计和动态分析提供可靠的理论依据和科学的确定方法。该方法根据气体稳定流动的运动方程，在平均参数法基础上，考虑了动能变化时计算井筒压力分布的数学模型。该方法具有较好的稳定性和收敛性，适合计算机编程。
附图说明
[0061]
图1为实施例中注气井井筒压力分布的确定方法的流程图。
具体实施方式
[0062]
为了更清楚地说明本发明，下面结合优选实施例对本发明做进一步的说明。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本发明的保护范围。
[0063]
如图1所示，本发明实施例提供一种注气井井筒压力分布的确定方法，该确定方法包括以下步骤：
[0064]
s100、获取注气井井筒的测试参数；
[0065]
s200、取注气井井筒的井口为坐标原点，竖直向下为正，按照预定计算步长将所述注气井井筒划分为多个微元段；
[0066]
s300、根据所述测试参数，以注气井井筒井口处为初始计算点，依次计算各微元段出口处的压力值和温度值；
[0067]
s400、根据所述各微元段出口处的压力值确定所述注气井井筒压力分布。
[0068]
具体的，所述测试参数包括油管内径d、井口油压p0、气井深度h、井口温度t0、井底温度te、地温梯度g
t
、注气量q
sc
、天然气与空气的相对密度γ。
[0069]
该确定方法基于以下假设条件：
[0070]
(1)流体流动状态为稳定单向流动，流体为天然气单相流；
[0071]
(2)井筒内传热为稳定传热；
[0072]
(3)井筒内天然气流体符合能量守恒规律；
[0073]
(4)地层温度按线性变化，已知地温梯度；
[0074]
(5)已知井口温度、压力等测试参数。
[0075]
具体的，s300的过程包括：
[0076]
s301、以注气井井筒井口处为初始计算点，任意设定步长dl，l＝0，i＝1。
[0077]
s302、令p
′1＝p0；t
′1＝t0。
[0078]
s303、令
[0079]
s304、计算该微元段的平均压力和平均温度：
[0080]
s305、根据该微元段的入口压力和入口温度得到关于压缩因子、温度和压力的简化变量i1；根据该微元段的出口压力和出口温度得到关于压缩因子、温度和压力的简化变量i2。
[0081]
具体的，通过公式(5)计算得到：
[0082][0083]
其中，p为压力，t为温度，z为压缩因子。
[0084]
s306、将i1和i2代入井筒压力计算模型式中计算得到dl段出口处的井筒压力p2。
[0085]
具体的，所述井筒压力计算模型为式4)：
[0086][0087]
p1为dl段入口处的井筒压力，mpa；p2为dl段出口处的井筒压力，mpa；θ为dl井段与水平方向的夹角；a1至a4为压力模型的回归系数；q
sc
为注气量，m3/d；d为油管内径，m；g为重力加速度，9.81m/s2；γ为天然气与空气的相对密度；f为天然气的摩擦阻力系数。
[0088]
其中：a1至a4为压力模型的回归系数，本发明经过大量现场实际测试资料拟合得到a1至a4的取值，如表1所示：
[0089]
表1提高压力精度的回归系数
[0090]
系数a1a2a3a4系数值2.2950.8699.17
×
10-3
1.853
[0091]
s307、计算压力差dp＝p
2-p1；计算天然气在p1和p2压力条件下的真实流速v1和v2，
dv＝dv
2-dv1。
[0092]
s308、当|t2一t2′
|≤0.001，|p
2-p2′
|≤0.001时，结束循环；否则，令t2′
＝t2、p
′2＝p2返回s304再次计算。
[0093]
s309、i++，令p
′1＝p2，t
′1＝t2，l＝l+dl，返回s303继续迭代计算下一微元段的温度、压力；当l≤h时则迭代结束。
[0094]
本发明实施例通过合理的假设和必要的简化，提出了一个修正改进的算法，可为储气库井，特别是定向井、高注气量井的工程设计和动态分析提供可靠的理论依据和科学的确定方法。该方法根据气体稳定流动的运动方程，在平均参数法基础上，考虑了动能变化时计算井筒压力分布的数学模型。该方法具有较好的稳定性和收敛性，适合计算机编程。
[0095]
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。