本实用新型属于天然气地下储气库井口注采双向计量技术领域领域，更具体地，涉及一种利用节流阀节流原理进行计量的装置，主要用于天然气地下储气库在注气和采气过程中，实现单井的实时、在线、连续、气液分相的双向计量功能，属于一种内部计量器具。
背景技术：
天然气地下储气库(简称“储气库”)是确保向天然气用户平稳供气的最佳选择，具有储存量大、安全性高、使用年限长、经济易行等优点，是目前较为普遍采用的一种天然气存储设施。储气库地面井组中各单井的注气计量数据和采气计量数据是评价该区域地层天然气储存能力和采出能力的重要依据。储气库注采单井要求双向计量、承载高低压交替变化，而且在注采气成分、压力、流量、相态、流态变化较大条件下仍然能保证较高的计量精度及可靠度等特点，对计量仪表提出了较为苛刻的要求，因此选用适合储气库运行工况特点的计量装置至关重要。
目前储气库井口注采双向计量多采用靶式流量计和超声流量计。靶式流量计根据伯努利方程，当介质在测量管中流动时，因其自身的动能与靶板产生压差，而产生对靶板的作用力，使靶板产生微量的位移，其作用力的大小与介质流速的平方成正比。这样，传感器获得的信号可换算为流体的流量。靶式流量计结构简单，没有滑动、转动可动部件，主要用于测量高粘度、低雷诺数的流体。当靶式流量计用于储气库井口双向计量时，由于靶板安装于管道中，需要与气流接触，在注采气过程中，固体颗粒、油、水等固液态杂质对靶板形成交变载荷，容易受气流流速、气流冲蚀、杂质干扰、采出液波动等复杂工况的影响，造成计量误差偏大。
超声流量计主要通过信号传输时差校正原理来测量介质流量，其换能器在逆流和顺流方向同时发射和接收信号，当介质流动时，信号在介质中传播会跟随流动产生位移，其传播的时间在顺流方向小于逆流方向，而此时间差与平均介质流动的速度成正比例关系，通过测量此时间差，即可计算出介质流量。超声换能器与气体直接接触时，称为插入式超声流量计；超声换能器不与气体直接接触时，称为外夹式超声流量计。超声波流量计具有测量精度高、可靠性强、无扰动部件、压力损失小等特点，尤其适用于大管径气体流量的高精度计量(管径可达1.6m、精度0.5％)，但在储气库的实际应用过程中，天然气中的杂质及液滴、换能器定位误差、声学噪声干扰等因素会对测量结果造成误差。根据《用气体超声流量计测量天然气流量(GB/T 18604—2014)》国家标准，在正常输气工作条件下，在流量计表体内的附着物(如凝析液或带有加工杂质的油品残留物、灰和砂等) 会减少流量计的流通面积而影响计量准确度，同时附着物还会阻碍或衰减超声换能器发射和接收超声信号，或者影响超声信号在流量计表体内壁的反射，因此对流量计应定期检查清洗。
以上靶式流量计和超声流量计(含插入式超声流量计和外夹式超声流量计)的性能对比见表1。从现场应用情况来看，注气阶段的计量难度不大，出现的问题较少；在采气初期，天然气气体中携带油分或者凝析地层水及岩屑杂质等，影响流量计的测量精度，容易造成气量计量误差，当流量计承载静压差较大时，极易导致流量计阻流件卡阻、损坏或脱落。此外，靶式流量计和超声流量计均不能对注采井的凝析液或地层水进行计量。
表1储气库双向流量计对比表
技术实现要素：
针对上述问题，本实用新型的目的提供一种具有数据自动采集功能的计量装置，为建设数字化储气库创造条件。实现在复杂气流环境中单井注采同管的双向、实时、在线、连续、气液分相计量，既可满足储气库关于内部计量的精度及可靠性要求，又能大幅度降低计量成本。
为实现上述目的，本实用新型提供一种天然气地下储气库单井注采同管双向计量装置，所述装置包括：
依次连接的第一温压测取直管、节流阀和第二温压测取直管，所述第一温压测取直管的一端连通于所述储气库，另一端连通于所述节流阀，所述第二温压测取直管的一端伸出所述储气库，另一端连通于所述节流阀；
压力传感器，所述压力传感器分别设置在所述第一温压测取直管和第二温压测取直管上；
温度传感器，所述温度传感器分别设置在所述第一温压测取直管和第二温压测取直管上；
计量处理装置，所述计量处理装置通信连接于所述节流阀、所述压力传感器及所述温度传感器，用于计算气体流量。
进一步地，所述第一温压测取直管和第二温压测取直管垂直设置，所述节流阀设置在拐角处。
进一步地，还包括标定接头直管，所述标定接头直管连接于所述第二温压测取直管的一端。
进一步地，标定接头直管与所述第二温压测取直管为可拆卸的安装。
进一步地，所述第一温压测取直管、所述节流阀、所述第二温压测取直管、所述标定接头直管依次通过法兰严密连接而成。
进一步地，所述计量处理装置为远程计算机工作站或单片机。
进一步地，所述第一温压测取直管的一端连通于储气库作为采气入口和注气出口，所述第二温压测取直管的一端伸出储气库作为注气入口和采气出口。
进一步地，所述节流阀设置有开度传感器。
本实用新型的装置结构简单，工作范围宽广且性能稳定，不受流态、相态、凝析液、固体杂质等因素的干扰。设置了标定接头，可方便地采用移动计量装置在不停产的条件下对本装置进行在线调试和定期标定，也可以采用下游生产分离器在启停井期间对本装置进行校准。
而且，本实用新型的装置可同时为储气库(群)的多口井分别工作，进一步降低了计量成本，提高了运维管理效率。
本实用新型的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本实用新型示例性实施方式进行更详细的描述，本实用新型的上述以及其它目的、特征和优势将变得更加明显，其中，在本实用新型示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本实用新型的计量装置结构示意图。
图2示出了根据本实用新型的一个实施例的计量方法流程图。
图3示出了根据本实用新型的一个实施例的节流过程示意图。
附图标记说明：
1、第一温压测取直管；2、节流阀；3、第二温压测取直管；4、标定接头直管；5、计量处理装置；D、管道内直径；P1、P2、测取点压力；T1、 T2、测取点温度；CD、节流阀开度；Qg、Ql、气相、液相的流量。
具体实施方式
下面将更详细地描述本实用新型的优选实施方式。虽然以下描述了本实用新型的优选实施方式，然而应该理解，可以以各种形式实现本实用新型而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本实用新型更加透彻和完整，并且能够将本实用新型的范围完整地传达给本领域的技术人员。
下面结合附图和实施例对本实用新型进行详细的描述。
如图1所示，本实用新型提供的一种天然气地下储气库单井注采同管双向计量装置，包括第一温压测取直管1、节流阀2(例如，角式针阀、笼套阀等)、第二温压测取直管3、标定接头直管4、计量处理装置5(例如，远程计算机工作站或安装在设备上的单片机等)，及计量处理装置5安装有计量应用程序。
第一温压测取直管1、节流阀2、第二温压测取直管3、标定接头直管4依次通过法兰严密连接而成。标定接头直管4用连接移动计量装置(如车载旋风分离计量撬)的软管接头。
压力传感器分别设置在第一温压测取直管1和第二温压测取直管3上， P1和P2分别为第一温压测取直管1和第二温压测取直管3上压力传感器的压力值。温度传感器分别设置在第一温压测取直管1和第二温压测取直管3 上，T1和T2分别为第一温压测取直管1和第二温压测取直管3上温度传感器的温度值。第一温压测取直管1的一端连通于储气库作为采气入口(注气出口)，另一端连通于节流阀2，第二温压测取直管3的一端伸出储气库作为注气入口(采气出口)，另一端连通于节流阀2。节流阀2设置有开度传感器。
计量处理装置5通信连接于温压传感器、节流阀开度传感器并保持通讯状态，通过计量处理装置5的计量应用程序实时演算完成计量任务。本实用新型的装置可自动识别并完成储气库在注气、采气阶段，天然气不同流向、流量、流态、相态条件下的单井计量功能。
通过计量处理装置5实时在线采集第一温压测取直管1、第二温压测取直管3及节流阀2的开度传感器的数据，导入计量应用程序在线演算完成。
根据本实用新型的另一实施方式，提供一种天然气地下储气库单井注采同管双向计量方法，所述方法包括：
采集压力、温度及节流阀开度数据；
根据天然气组分以及采集的温度、压力数据的平均值，计算当前条件下天然气混合物的气相质量分数、液相质量分数、含水率及物性参数；
判别当前的作业流程；
如果当前为采气作业流程，则计算单井的气相流量和液相流量；
如果当前为注气作业流程，则计算单井的气体流量。
接下来，参照图2描述本实用新型的具体实施例。参照图2，计量方法的实施流程如下：
(1)输入天然气的组分(CH4、C2H6、C3H8、C4H10、H2O等)及其摩尔分数，并进行定期监测。本步骤并不是每次计量都需要输入的，如果组分发生较大变化，及时更新即可。
(2)采集井口节流阀前后的压力、温度及节流阀开度等实时数据。例如，可以基于OPC通讯协议，借助SCADA/DCS控制系统采集实时数据。
(3)根据天然气组分以及实时采集的温度压力数据的平均值，利用相平衡热力学计算模型计算当前条件下，天然气混合物的气、液相质量分数、含水率及物性参数等。物性参数包括密度、摩尔质量、比体积、比焓、比熵、比热、比内能、压缩因子、JT系数等参数。相平衡热力学计算模型选用本领域公认成熟的SRK模型、PR模型等，在此无需赘述。例如，可以通过计量应用程序中的“相平衡热力学计算模块”进行上述的计算过程。
(4)由节流阀前后压力P1、P2的相对大小，判别当前的作业流程。若 P1>P2,则为采气作业；反之则为注气作业。
(5)如果当前为采气作业状态，天然气中夹带有一定量的凝析液或地层水，计算单井的气相、液相流量，其理论雏形来源于Perkins模型，考虑了气液相之间的滑脱影响。例如，可以通过计量应用程序中的“多相流动力学计算模块”进行上述的计算过程。图3示出了节流过程示意图，参照图3，计算过程如下：
引入气液相之间的滑脱比R：
式中，ug、uL分别为气、液相的当地流速。
对亚临界流，滑脱比可用Grolmes&Leung模型描述：
式中，αi(i＝0,1,2,3)为拟合常数；xg为气相的质量分数；ρg、ρL分别为气相、液相的密度；μg、μL分别为气相、液相的动力粘度。
考虑气液相之间的滑脱后，其混合密度ρm为：
假设气体在节流过程中作等熵膨胀运动，则对气液相混合物定义等熵膨胀指数：
式中，k为气体绝热指数；Cvg、CL分别为气、液相的定容比热容。
根据热力学等熵运动的状态方程，则有：
npdv＝-vdp (5)
式中，p为压力；v为比体积。
积分上式，可得如下平衡关系：
根据Perkins模型，忽略气体内能影响，节流过程的气液多相流动的能量平衡方程可表述为：
式中，为气液相的混合质量流量。
将式(3)、式(6)代入式(7)，可得：
式中，vg、vL分别为气、液相的比体积。
定义压力比r为：
将混合密度表达式(3)代入式(8)得：
同时，定义α：
根据式(6)可得：
将式(11)、式(12)代入式(10)，有：
进一步整理式(13)后，得：
当时，式(14)可进一步化简为：
当流态为临界流时，应满足以下微分方程：
积分上式，化简得：
当时，式(17)可进一步简化为：
对式(17)或式(18)，采用迭代法即可解出rc。
当r≥rc时，流态为亚临界流，利用式(17)或式(18)计算质量流量当r＜rc时，流态为临界流，其质量流量采用以下方法计算：
如果可用下式进行简化计算：
计算出质量流量后，气相、液相的体积流量由以下方程解得：
式中，Qgsc为气体在标况下的体积流量，ρgsc为气体在标况下的密度，C 为分配系数，为质量流量，xg为气相的质量分数，QL为液相的体积流量，ρL为液相的密度。
(6)如果当前为注气作业状态，节流阀一般处于全开状态，即CD＝1.0，注入的天然气为纯度较高的干气，节流压降较小，计算单井的气体流量(例如，可以通过采用计量应用程序中的“单相流动力学计算模块”进行计算)，其计算过程如下：
根据热力学原理，临界压力比为：
当r≤rc时，流态为临界流；否则为亚临界流。
根据气体节流的等熵原理，对于亚临界流状态，流量与压力比的关系可表示为：
对于临界流动，最大气流量与p2无关，按下式进行计算：
式中，Qgsc为气体在标况下的体积流量，d为节流阀孔眼直径，Z为气体偏差系数，γg为天然气的相对密度，P1、P2为节流阀前后的压力，r为压力比，k为气体绝热指数，T1为节流阀前的测取点温度。
根据上面的计算结果，显示和存储流量计量结果。
此外，为确保计量始终处于较高精度水平，需定期对装置进行标定。根据计量性能的实际表现，标定周期一般为1～3年。两种可选的计量标定数据来源如下：
1、拆下标定接头4，连接移动计量装置(如车载旋风分离计量撬)的软管接头，该计量装置和本实用新型的装置同时计量该井的注、采气量。
2、关井，测量稳定后下游生产分离器(多井共用)的流量；然后启井，测量稳定后下游生产分离器(多井共用)的流量。启井、停井期间下游分离器的流量差值即为该井的实测产量，启井阶段本实用新型装置同时计量该井采气量。
上述实施例仅用于说明本实用新型，其中通过不同节流方式或节流级数实现储气库单井天然气注采计量的各部件的结构、连接方式和制作工艺等，以及计量软件所采用的计算模型都是可以有所变化的，凡是在本实用新型技术方案的基础上进行的等同变换和改进，均不应排除在本实用新型的保护范围之外。
以上已经描述了本实用新型的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。