发明名称 双差压节流湿气测量方法
摘要 本发明属于多相流量测量领域,涉及一种利用双节流装置实现的湿气测量方法,首先进行实验室标定,获取双节流装置的简单虚高模型和精确虚高模型,然后再利用两节流差压湿气测量装置进行在线测量,获得两节流装置的差压信号,然后进行多参数迭代计算实现湿气测量。本发明提出的测量方法能够实现在线测量,实施简单、方便,具有较高的测量精度,完全符合工业现场的需要,可广泛应用于石油天然气工业现场。
申请公布号 CN101413817A 申请公布日期 2009.04.22
申请号 CN200810153806.9 申请日期 2008.12.05
申请人 天津大学 发明人 徐英;张强;张涛
分类号 G01F1/88(2006.01)I;G01F7/00(2006.01)I 主分类号 G01F1/88(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 代理人 江镇华
主权项 1.一种湿气测量方法,采用的测量系统包括安装在测量管路里的锥形节流装置和文丘里节流装置,在测量管路里安装有压力变送器,在锥形节流装置的首尾处安装有第一压差变送器,在文丘里节流装置的节流管路内外安装有第二压差变送器,在其测量管路下游内还安装有温度变送器,进行湿气测量的步骤包括:(1)对于锥形节流装置,通过对实验数据进行曲线拟合,求得锥形节流装置的湿气虚高模型式,即其与L-M参数和Froude数之间的关系式:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mo>-</mo><mi>cone</mi></mrow></msub><mo>=</mo><msub><mi>b</mi><mn>1</mn></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mi>X</mi><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mo>&CenterDot;</mo><mi>X</mi><mo>&CenterDot;</mo><msub><mi>F</mi><mi>rg</mi></msub><mo>+</mo><msub><mi>b</mi><mn>4</mn></msub><mo>&CenterDot;</mo><msub><mi>F</mi><mi>rg</mi></msub><mo>+</mo><msub><mi>b</mi><mn>5</mn></msub><mrow><mo>(</mo><mfrac><msub><mi>&rho;</mi><mi>l</mi></msub><msub><mi>&rho;</mi><mi>g</mi></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中的L-M参数与Froude数分别为<maths num="0002"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfrac><msub><mi>W</mi><mi>l</mi></msub><msub><mi>W</mi><mi>g</mi></msub></mfrac><msqrt><mfrac><msub><mi>&rho;</mi><mi>g</mi></msub><msub><mi>&rho;</mi><mi>l</mi></msub></mfrac></msqrt></mrow></math>]]></maths>和<maths num="0003"><![CDATA[<math><mrow><msub><mi>F</mi><mi>rg</mi></msub><mo>=</mo><mfrac><msub><mi>U</mi><mi>sg</mi></msub><mrow><msqrt><mi>gD</mi></msqrt></mrow></mfrac><msqrt><mfrac><msub><mi>&rho;</mi><mi>g</mi></msub><mrow><msub><mi>&rho;</mi><mi>l</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>g</mi></msub></mrow></mfrac></msqrt></mrow></math>]]></maths>,式中,U<sub>sg</sub>为气相表观速度,W<sub>g</sub>为气体质量流量,W<sub>l</sub>为液体质量流量,C为仪表流出系数,ε为气体的可膨胀系数,β为节流比,D为管道直径,ρ<sub>g</sub>为气相流体密度,ρ<sub>l</sub>为液相流体密度,<img file="A200810153806C00024.GIF" wi="67" he="34" />为通过实验测得的锥形节流装置的压差信号;(2)对于文丘里节流装置,通过对实验数据进行曲线拟合,求得文丘里节流装置的湿气虚高模型式,即其与L-M参数和Froude数之间的关系式:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mo>-</mo><mi>venturi</mi></mrow></msub><mo>=</mo><msqrt><mn>1</mn><mo>+</mo><mi>CX</mi><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></msqrt></mrow></math>]]></maths>,式中,<maths num="0005"><![CDATA[<math><mrow><mi>C</mi><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>&rho;</mi><mi>l</mi></msub><msub><mi>&rho;</mi><mi>g</mi></msub></mfrac><mo>)</mo></mrow><mi>n</mi></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>&rho;</mi><mi>g</mi></msub><msub><mi>&rho;</mi><mi>l</mi></msub></mfrac><mo>)</mo></mrow><mi>n</mi></msup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>n</mi><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>F</mi><mi>rg</mi></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><mo>&CenterDot;</mo><msqrt><msub><mi>F</mi><mi>rg</mi></msub></msqrt><mo>+</mo><msub><mi>a</mi><mn>3</mn></msub><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mn>0.5</mn><mo>&le;</mo><msub><mi>Fr</mi><mi>g</mi></msub><mo>&le;</mo><mn>1.5</mn><mo>;</mo></mrow></math>]]></maths>(3)对于两个节流装置,分别通过对实验数据进行曲线拟合,求得锥形和文丘里节流装置的简单虚高模型式分别为:Φ<sub>g-cone</sub>=a+b·X、<maths num="0007"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mo>-</mo><mi>venturi</mi></mrow></msub><mo>=</mo><msqrt><mn>1</mn><mo>+</mo><mi>CX</mi><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></msqrt></mrow></math>]]></maths>(4)由压力变送器和温度变送器获取管道内流体的压力和温度,计算管道内气相流体的密度ρ<sub>g</sub>;(5)利用在锥形节流装置和文丘里节流装置上安装的两个压差变送器获得的压差信号,根据下列公式,分别计算两个节流装置的虚高质量流量W<sub>TP1</sub>和W<sub>TP2</sub>:<maths num="0008"><![CDATA[<math><mrow><msub><mi>W</mi><mi>tpk</mi></msub><mo>=</mo><mfrac><mrow><mi>C</mi><mo>&CenterDot;</mo><mi>&epsiv;</mi></mrow><msqrt><mn>1</mn><mo>-</mo><msup><mi>&beta;</mi><mn>4</mn></msup></msqrt></mfrac><mo>&times;</mo><mfrac><mi>&pi;</mi><mn>4</mn></mfrac><msup><mi>&beta;</mi><mn>2</mn></msup><msup><mi>D</mi><mn>2</mn></msup><mo>&times;</mo><msqrt><mn>2</mn><mi>&Delta;</mi><msub><mi>P</mi><mi>tpk</mi></msub><mo>&CenterDot;</mo><msub><mi>&rho;</mi><mi>g</mi></msub></msqrt><mo>,</mo></mrow></math>]]></maths>式中k=1,2;(6)根据公式Φ<sub>g</sub>=Φ<sub>g</sub>(X),<maths num="0009"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><mfrac><mrow><msub><mi>W</mi><mrow><mi>tp</mi><mn>1</mn></mrow></msub><mo>/</mo><msub><mi>W</mi><mi>g</mi></msub></mrow><mrow><msub><mi>W</mi><mrow><mi>tp</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>W</mi><mi>g</mi></msub></mrow></mfrac><mo>=</mo><mfrac><msub><mi>&Phi;</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><msub><mi>&Phi;</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></math>]]></maths>以及简单虚高模型式Φ<sub>g-cone</sub>=a+b·X、<maths num="0010"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>g</mi><mo>-</mo><mi>venturi</mi></mrow></msub><mo>=</mo><msqrt><mn>1</mn><mo>+</mo><mi>CX</mi><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></msqrt></mrow></math>]]></maths>,获得X参数的初值与气相质量流量的初值X<sub>0</sub>、W<sub>g0</sub>,将气相质量流量的初值W<sub>g0</sub>带入Froude数公式,获得Froude数的迭代初值F<sub>rg0</sub>。;(7)将X参数,Froude数的初值X<sub>0</sub>和F<sub>rg0</sub>带入文丘里节流装置的虚高模型式,实现对X参数的第一次修正,得到X<sub>1</sub>;(8)将修正的X参数即X<sub>1</sub>与Froude数F<sub>rg0</sub>带入到锥形节流装置的虚高模型式中,对气相质量流量进行修正,得到W<sub>g1</sub>,利用W<sub>g1</sub>计算气相表观速度U<sub>sgl</sub>,并代入式Froude数公式,修正F<sub>rg0</sub>,得到F<sub>rgl</sub>;(9)如此重复执行步骤(7)和(8),实现对X参数、Froude数和气相质量流量W<sub>g</sub>的迭代修正,当W<sub>g</sub>的迭代值的相对变化小于设定值,则认为迭代结果收敛,结束迭代。
地址 300072天津市南开区卫津路92号天津大学
您可能感兴趣的专利