一种基于喷嘴的低产气井湿气两相流量计量方法与流程

1.本发明涉及天然气计量技术领域，特别涉及一种基于喷嘴的低产气井湿气两相流量计量方法。


背景技术：

2.湿天然气井的气液两相流动十分复杂，两相流形成机理和运动规律尚未完全被认识，两相流参数的测量难度较大，需要建立合适的数学和物理模型，选择特征参数，研发特种仪器以及采用先进的数据处理方法。
3.现有气液两相流的参数测量多采用分离法，但是该方法所用分离设备体积庞大，占地面积大且费用较高。节流式差压流量计的结构简单、容易替换且可靠性较高，因此基于节流式差压流量计的两相流测量模型的研究一直受到大量研究者的重视。murdock、bizon和collins、v.c.ting等，对多相流相关模型做了大量的实验研究，提出很多修正公式。贺登辉、白博峰等人分析了两相流的质量流量系数k的影响因素，并提出了一种新的两相流测量模型。lockhart和martinelli根据两相流在管内的摩阻压降规律定义了与两相流的液相含率相关的参数，洛玛参数(xlm)；克列姆莱夫斯基用不同节流比孔板流量计在假设两相流为均相流的基础上，利用大干度的水和水蒸气混合的流体进行了大量实验，通过实验得到两相流的压降关系式，并得到按绝热膨胀计算流量时的修正因子k。r.v.smith和j.t.leang认为气液两相的压降高于单相气的压降是由于液相对气相的阻塞，因此，提出用阻塞系数来修正两相流的计量流量。chisholm利用孔板流量计进行水和水蒸气的两相流实验，推导得到滑动比与密度比的关系，证明了两相流的压降大小与含气率和工作压力有关；林宗虎基于分相流动假设提出文丘里管的气液两相流流量计算相关式，并通过改变修正系数的方法将适用于孔板流量计的压降模型应用到文丘里流量计中。岳伟挺、张宏建等人认为两相流的压差波动与两相流的空隙度有关，并进行相关实验。研究结果表明，空隙度从0到1的过程中，差压波动先升高后降低，波动大小主要取决于空隙率，同时也受流量、压力的影响。
4.近些年国内外能源和仪表类企业也开始投入研究，希望能够打造计量精准的多相流产品，有了很多的研究成果，但多数只能在一定工况条件下计量准确，且各方法均存在一定的短板。比如放射法计量相对精准，但造价高、结构复杂、存在安全顾虑；文丘里节流法受硬件结构限制，节流内径不能太小，因此对于低产气井无法测量，无法形成有效差压进行计量。


技术实现要素：

5.针对现有技术中的问题，本发明的目的在于提供一种基于喷嘴的低产气井湿气两相流量计量方法，实现低产气井中气液两相流量的准确计量。
6.为了实现上述目的，本发明采用以下技术方案予以实现：
7.一种基于喷嘴的低产气井湿气两相流量计量方法，包括以下步骤：
8.(1)采集管道内待计量湿气的温度t和压力p，利用经验公式计算待计量湿气中气
相密度ρg和粘度μg；
9.ρg＝f(t，p)；
10.μg＝f(t，p)；
11.(2)以一定时间为一个区间，在该一定时间内采集1个待计量湿气的温度值t，30个压力值p[30]，3000个差压值
△
p[3000]；
[0012]
对以上数据进行滤波处理，剔除掉过高过低的异常值，然后进行均值计算得到压力p
tp
，差压
△
p
tp
；
[0013]
p
tp
＝average(p[30])
[0014]
△
p
tp
＝average(
△
p[3000])
[0015]
(3)使用差压
△
p
tp
计算产气量qg：
[0016][0017]
其中，qg是产气量；c是喷嘴流出系数；β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；ρg是待计量湿气中气相密度；
[0018]
(4)计算
△
p[3000]的标准差
△
p
st
：
[0019]
△
p
st
＝stdev(
△
p[3000])
[0020]
(5)设气相单独流过节流件的差压为
△
pg；建立
△
pg与
△
p
tp
、
△
p
st
的关系式：
[0021][0022][0023]
δpg＝δp
tp-θδp
st
[0024]
δpg＝δp
tp-(az2 bz c)δp
st
[0025]
其中，a，b，c均为实验拟合系数；
[0026]
ρg是待计量湿气中气相密度，ρ
l
是待计量湿气中液相密度；z是设定的中间参数，无量纲变量；
[0027]
(6)修正后的气相流量为：
[0028][0029]
其中，q
gx
是修正后的气相流量；c是喷嘴流出系数；β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；
[0030]
(7)气相弗劳德数表示重力对流动影响的程度，与表观速度、密度相关，即受流量与压力的影响，气相弗劳德函数定义为气相表观惯性力与液相重力之比的均方根；
[0031][0032]
[0033]
其中，v
sg
为气相表观速度；g为重力加速度；d为管道内径；s为设备管道横截面积。
[0034]qgx
是修正后的气相流量；frg是气相弗劳德数；
[0035]
(8)设液相单独流过节流件的差压为
△
p
l
；建立
△
p
l
与
△
p
tp
的模型：
[0036]
δp
l
＝(d*frg2 f*frg g)δp
tp
[0037][0038]
其中，d，f，g为实验拟合系数；
[0039]ql
是液相流量，c是喷嘴流出系数，β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；
[0040]
(9)设质量含气率为：
[0041][0042]
其中，x
hg
是质量含气率，q
gx
是修正后的气相流量，q
l
是液相流量；
[0043]
(10)湿气两相体积流量分别为：
[0044][0045][0046]
其中，q
气
为标况气相体积流量，q
液
为液相体积流量，q
gx
是修正后的气相流量，q
l
是液相流量，ρ
l
待计量湿气中液相密度，ρ
gb
为标况下天然气密度。
[0047]
与现有技术相比，本发明具有以下技术效果：
[0048]
1、本发明提供的基于喷嘴的湿气两相流量计量方法，通过利用高频差压进行实时数据采集，使得数据更加反映正式流动状态，计算更加可靠；利用差压均值和标准差，间接计算单相气和单相液流经节流件的差压，进而计算气液两相的流量，计算方法新颖独特，同时确保了对气液两相流量计算的准确性。
[0049]
2、本发明提供基于喷嘴的湿气两相流量计量方法特别适用于低产气井的湿气两相流量的计量。
附图说明
[0050]
图1示出为根据本发明具体实施方式提供的一种基于喷嘴的低产气井湿气两相流量计量方法流程图；
[0051]
图2示出为根据本发明具体实施方式提供的一种基于喷嘴的低产气井湿气两相流量计量装置的结构示意图；
[0052]
图3示出为图2中计量装置的侧视图；
[0053]
图4示出为图2中k1-k1截面的剖视图；
[0054]
图5示出为图2中计量装置的俯视图；
[0055]
图6示出为图2中k2-k2截面的剖视图；
[0056]
图中标号说明：10、设备管道；11、连接法兰；12、排污口；13、排污阀；14、护罩外壳；20、压力传感器；30、温度传感器；40、差压传感器组；50、主电路组件；60、电源电路板；70、液
晶显示器。
具体实施方式
[0057]
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体附图和实施例，进一步阐明本发明。
[0058]
现有技术中，油气开发现场仍以分离式设备进行湿气两相流量的计量，但分离式设备体积大、成本高，更多的适用于场站计量，而对于单气井计量，现主要有文丘里节流法和射线加节流法两种；其中，文丘里节流尺寸不能过小，对日产量在1.5万以下的气井几乎无法计量，且文丘里节流设备长度较长，现场要求条件高；射线法存在放射性，不符合现场安全规范，且射线法成本较高，对于低产量气井成本效益低。
[0059]
基于此，结合图1所示，本技术的发明人提出了一种基于喷嘴的低产气井湿气两相流量计量方法，包括以下步骤：
[0060]
s1：采集管道内待计量湿气的温度t和压力p，利用经验公式计算待计量湿气中气相密度ρg和粘度μg；
[0061]
ρg＝f(t，p)
[0062]
μg＝f(t，p)
[0063]
s2：以一定时间为一个区间，在该一定时间内采集1个待计量湿气的温度值t，30个压力值p[30]，3000个差压值
△
p[3000]；对以上数据进行滤波处理，剔除掉过高过低的异常值，然后进行均值计算得到压力p
tp
，差压
△
p
tp
；
[0064]
p
tp
＝average(p[30])
[0065]
△
p
tp
＝average(
△
p[3000])；
[0066]
根据本发明提供的计量方法，本发明中，所述的一定时间可以在一定范围内进行选择，作为优选的，所述一定时间为10～30s，更优选为15s。
[0067]
s3：使用差压
△
p
tp
计算产气量qg：
[0068][0069]
其中，qg是产气量；c是喷嘴流出系数；β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；ρg是待计量湿气中气相密度；
[0070]
s4：计算
△
p[3000]的标准差
△
p
st
：
[0071]
△
p
st
＝stdev(
△
p[3000])
[0072]
s5：设气相单独流过节流件的差压为
△
pg；建立
△
pg与
△
p
tp
、
△
p
st
的关系式：
[0073][0074][0075]
δpg＝δp
tp-θδp
st
[0076]
δpg＝δp
tp-(az2 bz c)δp
st
[0077]
其中，a，b，c均为实验拟合系数；
[0078]
ρg是待计量湿气中气相密度，ρ
l
待计量湿气中液相密度；z是设定的中间参数，无量纲变量；
[0079]
s6：修正后的气相流量为：
[0080][0081]
其中，q
gx
是修正后的气相流量；c是喷嘴流出系数；β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；
[0082]
s7：气相弗劳德数表示重力对流动影响的程度，与表观速度、密度相关，即受流量与压力的影响，气相弗劳德函数定义为气相表观惯性力与液相重力之比的均方根；
[0083][0084][0085]
其中，v
sg
为气相表观速度；g为重力加速度；d为管道内径；s为设备管道横截面积。
[0086]qgx
是修正后的气相流量；frg是气相弗劳德数；
[0087]
s8：设液相单独流过节流件的差压为
△
p
l
；建立
△
p
l
与
△
p
tp
的模型：
[0088]
δp
l
＝(d*frg2 f*frg g)δp
tp
[0089][0090]
其中，d，f，g为实验拟合系数；
[0091]ql
是液相流量，c是喷嘴流出系数，β＝d/d，其中d为喷嘴内径，d为管道内径；ε是可膨胀系数；
[0092]
s9：计算质量含气率为：
[0093][0094]
其中，x
hg
是质量含气率，q
gx
是修正后的气相流量，q
l
是液相流量；
[0095]
s10：计算湿气两相体积流量分别为：
[0096][0097][0098]
其中，q
气
为标况气相体积流量，q
液
为液相体积流量，q
gx
是修正后的气相流量，q
l
是液相流量，ρ
l
待计量湿气中液相密度，ρ
gb
为标况下天然气密度。
[0099]
以实际应用举例，仪表管道内径为dn50，喷嘴节流内径为dn12，采集某气井的实际数据，计算分析如下：
[0100]
(1)以15秒为周期采集管道内待计量湿气的温度t，30个压力值p[30]，3000个差压值
△
p[3000]，整理数据得到：
[0101]
温度：t＝12.85℃
[0102]
平均压力：p
tp
＝average(p[30])＝1.33mpa
[0103]
平均差压：
△
p
tp
＝average(
△
p[3000])＝19.81kpa
[0104]
差压标准差：
△
p
st
＝stdev(
△
p[3000])＝3.46
[0105]
(2)利用经验公式计算待计量湿气中气相密度ρg和粘度μg；
[0106]
ρg＝f(t，p)＝11.06kg/m3[0107]
μg＝f(t，p)＝0.22mpa.s
[0108]
(3)使用差压
△
p
tp
计算产气量qg：
[0109][0110]
其中，qg是气相流量；c是喷嘴流出系数＝0.98；β＝d/d＝0.24，其中d为喷嘴内径＝12mm，d为管道内径＝50mm；ε是可膨胀系数＝0.99；ρg是待计量湿气中气相密度＝11.06kg/m3；
[0111]
(4)求气相单独流过节流件的差压为
△
pg；
[0112][0113]
θ＝f(z)＝az2 bz c＝0.86
[0114]
δpg＝δp
tp-θδp
st
＝19.81-0.86*3.46＝16.83kpa
[0115]
(5)求修正后的气相流量为：
[0116][0117]
其中，q
gx
是修正后的气相流量；c是喷嘴流出系数＝0.98；β＝d/d＝0.24，其中d为喷嘴内径＝12mm，d为管道内径＝50mm；ε是可膨胀系数＝0.99；
[0118]
(6)气相弗劳德数表示重力对流动影响的程度，与表观速度、密度相关，即受流量与压力的影响，气相弗劳德函数定义为气相表观惯性力与液相重力之比的均方根；
[0119][0120][0121]
其中，v
sg
为气相表观速度；g为重力加速度；d为管道内径；s为设备管道横截面积。q
gx
是修正后的气相流量；frg是气相弗劳德数；
[0122]
(7)求液相单独流过节流件的差压为
△
p
l
；
[0123]
δp
l
＝(d*frg2 f*frg g)δp
tp
＝0.0056*19.81＝0.11kpa
[0124][0125]ql
是液相流量，c是喷嘴流出系数＝0.98；β＝d/d＝0.24，其中d为喷嘴内径＝12mm，d为管道内径＝50mm；ε是可膨胀系数＝0.99；
[0126]
(8)求质量含气率为：
[0127][0128]
(9)求湿气两相体积流量分别为：
[0129][0130][0131]
其中，ρ
gb
为标况下天然气密度＝0.75kg/m3。
[0132]
(10)油田现场分离设备计量该气井的产气量为8012.7nm3/d，产液量为4.28m3/d。与计量设备对比分析：
[0133]
气相误差＝(7689.6-8012.7)/8012.7＝-4.03％
[0134]
液相误差＝(4.5-4.28)/4.28＝5.14％
[0135]
对于气井井口计量，气相和液相计量偏差控制在10％以内，可接受。
[0136]
结合图2-6所示，本发明还提供了一种基于喷嘴的湿气两相在线计量装置，所述计量装置包括设备管道10、压力传感器20、温度传感器30、差压传感器组40和主电路组件50。
[0137]
其中，所述设备管道10的两端通过连接法兰11串联至待计量湿气两相的管道中，所述设备管道10中安装有喷嘴节流件；所述压力传感器20设置在所述设备管道10中，用于测量进入喷嘴节流件前的待计量湿气两相的压力；所述温度传感器30设置在所述设备管道10中，用于测量进入喷嘴节流件前的待计量湿气两相的温度；所述差压传感器组40设置在所述设备管道10中，用于测量待计量湿气进入喷嘴节流件前后的差压；所述主电路组件50设置在所述设备管道10外部，用于数据采集、参数解析、液晶显示和数据传输。
[0138]
进一步地，根据本发明提供的计量装置，所述压力传感器20、温度传感器30、喷嘴节流件和差压传感器组40沿湿气两相的流动方向依次布置。也即，压力传感器20设置在最上游，温度传感器30设置在压力传感器20的下游，喷嘴节流件的上游；差压传感器组40设置在下游，如此，避免了温度传感器30对压力传感器20的测量产生影响。
[0139]
作为优选地，所述压力传感器20在设备管道10上水平安装，所述温度传感器30在设备管道10上垂直安装，所述差压传感器组40在设备管道10上呈45
°
安装，如此，确保了整个计量装置结构紧凑，满足多数现场井口的安装要求。
[0140]
进一步的，本发明提供的计量装置还包括有电源电路板60用于供电。
[0141]
所述主电路组件50包括主电路板和计算电路板，其中，所述主电路板与所述压力传感器20、温度传感器30和差压传感器组40电性连接，用于数据采集、参数解析、液晶显示和数据传输；
[0142]
所述计算电路板与所述主电路板电性连接，用于接收采集到的温度数据和压力数据，经计算得到湿气两相流量后反馈至主电路板，进行存储、显示和传输。
[0143]
进一步地，根据本发明提供的计量装置，所述设备管道10上开设有排污口12，所述排污口12上设有排污阀13用于清理设备管道10内杂物。作为优选地，所述排污口12和排污阀13设置在设备管道10的垂直下方，便于后续维修维护。
[0144]
本发明中，所述计量装置还包括护罩外壳14，所述护罩外壳14罩射在所述设备管道10外部，用于防护主电路组件50。
[0145]
本发明中，所述计量装置还包括液晶显示器70，所述液晶显示器70用于实时显示待计量湿气的采集参数和计量结果。作为优选地，该液晶显示器70嵌入都护罩外壳14中，省去了仪表表头，美观且保护了显示模块。
[0146]
本发明中，所述温度传感器30可选用本领域人员所常用的，具体如铂电阻温度传感器。
[0147]
本发明提供的计量装置，在具体使用时，将气井井口的工作管道切割，焊接配对法兰，通过配套螺栓即可实现连接，简单方便；并且计量装置整体结构紧凑，高度集成，在本发明的一个具体实施例中，设备全长仅为360mm，进而对井口空间的要求较低，可满足多数现场井口的安装要求。
[0148]
从井场电源引线24v连接设备，所述线路按防爆要求封装，在设备通电后，配置参数，检查设备，开井运行，温度、压力、差压采集实时数据，计算流量，液晶显示；并且，本发明提供的计量装置还可依据现场要求，可应用有线或无线技术将数据远传。
[0149]
本发明中，通过采用喷嘴作为节流件，在工艺上可以将喷嘴尺寸制造的更小，即节流件与管道内径的比值范围更宽，可以实现对日产量低于1.5万方的气井进行计量，也即，本发明特别适用于低产气井的湿气两相流量的计量。此外，喷嘴节流件拆卸灵活，前后直管段更短，使得计量设备整体结构紧凑，对现场井口安装空间要求低。
[0150]
本发明提供的计量系统属于节流法，相较于射线法设备的成本可控，更符合低气井的现场成本效益要求。本发明提供的方法中通过采用高频差压采集方式，快速采集差压数据，更好捕捉管道内湿气的流动和波动情况，使得数据更能体现真实流体变化，流量计算更加精准。以采集的差压数据和计算得到的差压标准差为基础，建立气相差压和液相差压修正模型，进而进行气液两相流量计算，方法新颖独特，实践证明具有较高的稳定性、可靠性和精准度。
[0151]
以上显示和描述了本发明的基本原理、主要特征和本发明的特点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明的范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。