Influence of Sand Production on Critical Liquid Carrying Flow in Gas Wells

Guijiu Zhou, Xiuwu Wang*, Yu Han, Jie Liu, Ruiquan Liao

1Faculty of Chemical Engineering, Kunming University of Science and Technology, Kunming Yunnan

2The Fourth Exploit Factory, Huabei Oilfield Company, Langfang Hebei

3Department of Petroleum Engineering, Yangtze University, Wuhan Hubei

4Multiphase Flow Laboratory of Gas Lift Test Base of China National Petroleum Corporation, Wuhan Hubei

Abstract

Keywords

1. 引言

对于生产中的气井，实时有效地判断气井积液情况，是保证气井能够稳定生产的前提。通常，气井积液判断是通过临界携液流量模型进行。不同地区气井生产情况不同，气井临界携液流量存在一定差异。

较早发展的模型是垂直气井方面的临界携液理论，这其中比较有代表性的是不同学者提出了一系列临界携液模型，如Duggan 模型[1]、Turner 模型[2]、Meshack 模型[3]、Coleman 模型[4]、Nosseir 模型[5]、杨川东模型[6]、李闽模型[7] [8]、吴志均模型[9]、Guo 模型[10]、王毅忠模型[11]以及魏纳模型[12]。针对不同实际情况，近几年的临界携液理论研究很多都是基于这三种模型的实验或者现场修正。由于液膜模型中的计算参数难以确定，现场很难推广应用，因此选用应用较为广泛的液滴模型[2]或BF 模型[13]，但在气流携液中完全没有考虑到气井出砂的影响，所以需对其进行修正。

针对实际产水气井积液判断模型的确定，一般需要结合室内气水流动模拟实验，进行实验理论研究得出气井气流临界携液的修正模型，再通过现场实际情况进行验证后用于气井积液判断。对于气井存在出砂的情况，本文通过室内气流携液携砂的模拟实验，建立了出砂情况下气井临界携液流量新模型，并对模型准确性进行了检验。

2. 实验研究

由于现场气井的生产情况无法进行可视化，而气井临界携液方面的理论分析研究仍然以理论与实验相结合的方法[13]。目前，这方面的实验研究中国内外主要参考Awolusi 和魏纳[14] [15]最为常见，针对本研究，自主设计了一套完整的气流携液携砂实验装置，该装置类似于国内实用新型专利[16]中的实验装置，但可以进行不同出砂速度下的气流携液模拟实验。

为了模拟该气田生产现场的实际情况，研究过程中对该气田地层水水质报告进行了计算统计，明确了实际地层水的水型和相关物性特征。该气田地层水的矿化度平均值为120,581 mg/L，结合室内配伍性原则[17]，配制与此矿化度值相接近的水作为实验液体介质。在标况下经测定配制后的模拟地层水，其主要物性参数如下：水型——氯化钙型CaCl2；密度——1.08 mg/L；黏度——17.4 mPa∙s；表面张力——53.37 mN。

2.1. 常规临界携液模型误差分析

首先结合常用的气井积液判断方法[18]，综合分析实际气井生产数据，统计积液的气井。以某气田实际No.1 气井为例分析气井井筒积液和确定临界积液点。

Figure 1. Gas well production dynamic curve 图1. 气井生产动态曲线

根据实际生产数据对No.1 井进行生产动态分析，由图1 中可以看出在2010 年4 月之后出现明显的“剪刀差”现象，初步判定该井可能在2010 年4 月积液。然后将计算井底流压和实测井底流压进行对比分析，结果由表1 可以看出，在2013 年12 月之前，实测值与计算值很接近，2013 年12 月之后，测压值都大于计算井底流压，认为积液时间在2012 年7 月至2013 年12 月之间，所以，该气井井筒开始产生积液的时间应该为2010 年4 月。

由于液膜模型中的计算参数难以确定，现场很难推广应用，因此选用应用较为广泛的液滴模型[2]或BF 模型[13]进行积液判断并优选。根据90˚倾角条件下的实验数据进行两个模型的计算值与实测值的误差分析，结果见图2 和图3。

Table 1. Comparative analysis of measured bottomhole pressure and calculated bottomhole pressure in Well No.1 表1. No.1 井实测井底流压与计算井底流压的对比分析

Figure 2. Error analysis of calculated and measured values of the model 图2. 模型计算值与实测值的误差分析

Figure 3. Correction of laboratory experiment data of Befoird model 图3. Befoird 模型室内实验数据修正

由图2 误差分析表明，采用Befoird 模型[19]计算的误差率较小，因此优选Befoird 模型[19]作为该气井积液判断数学模型。

2.2. 临界携液模型实验修正

采用Befoird 模型[13]对实验数据进行计算，将计算结果与实际数据进行对比，修正系数情况如图3 所示。将修正系数进行平均，对Befoird 模型[19]进行修正后的计算公式如下：

2.3. 临界携液模型现场校正

对该气田的气井进行积液统计，将积液气井的实际气体流量与上述修正模型计算结果进行对比，结果如表2。

Table 2. Model correction of critical fluid accumulation point of gas well 表2. 气井临界积液点对模型的校正

由于Befoird 模型[13]建模过程中忽略一些影响因素，这些影响因素很大程度上受气体流速的影响。根据实测/模型的平均值，对上述修正模型进行进一步校正，得到的临界携液流量的计算公式为：

2.4. 现场验证

统计2011 年~2017 年间该气田23 口井的液面测试数据对校正的模型进行验证，结果如表3 所示。虽然通过液面测试数据不能获得气井临界携液流量的大小，但可以准确判断气井是否积液[20]，进而可以准确地对修正模型进行的验证。

根据测试液面时对应的油压、井口温度、井底温度等参数可以应用上述校正模型及常用几种模型计算该状态下的临界携液流量，与该井实际产气量进行对比分析，以此判断该井的积液状况。

Table 3. Comparison and analysis of diagnosis results of wellbore effusion with different calculation formulas 表3. 不同计算公式井筒积液诊断结果对比分析表

由表3 可知，Turner 模型[2]、Coleman 模型[4]、LiMin 模型[7] [8]现场判断井筒积液的准确率分别为74%、74%、48%，本课题校正的临界携液模型判断准确率达到83%。说明校正后的模型能够准确地判断该气田气井的积液情况。

3. 出砂对临界携液模型的影响

3.1. 出砂时临界携液模型的建立

由于该气田出砂量较小，砂砾主要影响流体密度及摩阻，所以将含有砂砾的液体看作悬浮液体，根据固相颗粒存在的液体进行密度的修正。液体密度为考虑固体颗粒存在时的修正密度，建立密度修正计算公式如下：

将密度和粘度的修正结果纳入到上述校正的临界携液模型中，有：

3.2. 出砂时临界携液实验及数据的分析

由于该气田砂粒的粒径在600 目~1000 目之间，根据粒径中值及筛网目数与粒径之间的对应关系，实验过程选取600 目、800 目、1000 目及1200 目四种砂粒进行模拟实验。

进行了不同直径颗粒、不同出砂速度条件下的气流临界携液实验，共完成实验数据240 组，实验数据计算结果绝对误差如图4 所示。

Figure 4. Comparison of absolute error with and without sand model 图4. 考虑与不考虑砂模型绝对误差对比

由图4 可以看出，砂粒对气井临界携液流量存在一定影响，进行统计分析发现，用建立的出砂时临界携液模型计算结果的误差仅为4.5%，用上述校正的临界携液模型计算结果的误差为6.2%。

虽然室内采用四种不同直径的砂粒进行了实验，但出砂时临界携液模型的计算结果与不同粒径的室内实验结果十分接近，说明在砂粒流速一定的条件下砂粒粒径对临界携液流量影响不大，这是因为该气藏储层颗粒直径较小，对液体的密度和粘度影响较小，从而对校正的临界携液模型影响不大。

Figure 5. Impact analysis of sand flow and liquid flow 图5. 砂粒流量、液体流量的影响分析

出砂速度对临界携液气量的影响从图5 可以看出，在井底压力、温度一定的条件下，液体流量一定时临界携液流量随着出砂速度的增大而增大，但是增大幅度随着液体流量的增大而减小；出砂一定时，临界携液流量随着液体流量的减小而减小，但减小幅度降低。从总体上看，出砂对气井临界携液有一定影响。

3.3. 出砂时临界携液流量模型现场验证

将出砂时临界携液流量模型计算公式(4)与上述校正的临界携液模型计算公式(2)对气井积液的诊断结果进行对比分析，统计结果见表4。可以看出，与校正的临界携液模型分析结果相比，出砂时临界携液流量模型的精度提高了6.15%。

Table 4. Statistics of diagnostic results with and without considering sand particles affecting wellbore effusion 表4. 考虑和不考虑砂粒影响井筒积液的诊断结果统计

Continued

Continued

4. 结论

通过理论与室内实验研究发现：

1) 现场常用的临界携液流量模型并不一定适合所有气田及气井，针对不同区块气井需要结合实际情况进行模型优选及修正，以满足现场精度的要求。

2) 将含有较小粒径砂粒的液体看作悬浮液，基于这一思想建立了考虑砂粒影响的临界携液流量计算模型，出砂时临界携液模型的计算精度提高了2 个百分点，积液井的诊断准确率提高了6 个百分点。

3) 砂粒对气流临界携液流量存在一定的影响，但出砂速度较小的小颗粒对临界携液流量影响不大，大颗粒粒径的需要进一步进行理论与实验分析。