发明名称 两相流分相流速声电双模态测量方法
摘要 本发明提供一种两相流分相流速声电双模态测量方法,所采用的双模态传感器包括一套电学传感器和一套超声多普勒探头,测量方法为:利用电学传感器采集信号并获取水相含率、油相含率和两相流混合密度;确定超声多普勒探头的测量空间,计算测量空间的等效高度;计算超声接收探头获取的超声信号频率与超声发射探头所发射的激励超声信号频率间的频率差,并计算测量空间中的平均流速;建立管道中与距离管壁位置处流体流速的分布关系;建立测量空间中多普勒测量流速与边界层流速间关系;求解计算边界层流速;计算总流速,并分别得到水相流速和油相流速。本发明能够以非侵入方式准确获取两相流相含率。
申请公布号 CN104965104A 申请公布日期 2015.10.07
申请号 CN201510263691.9 申请日期 2015.05.21
申请人 天津大学 发明人 谭超;袁野;董虓霄;董峰
分类号 G01P11/00(2006.01)I;G01P5/00(2006.01)I 主分类号 G01P11/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 一种两相流分相流速声电双模态测量方法,所采用的双模态传感器包括一套电学传感器和一套超声多普勒探头;所述电学传感器用于获取两相流分相含率,超声多普勒探头用于获取两相流的流速;所述电学传感器与超声多普勒探头同时安装于管道中;所述超声多普勒探头包括一个超声发射探头和一个超声接收探头;所述超声发射探头和超声接收探头均与管道间以夹角θ安装,超声发射探头安装于管道顶部,超声接收探头安装于管道底部,并保证超声发射探头、超声接收探头与管道中心处于同一纵向截面内;该测试方法包含如下步骤:1)当两相流的连续相为导电相时,将电学传感器用作电导传感器获取水相含率ε<sub>w</sub>与油相含率ε<sub>o</sub>;当连续相为非导电相时,将电学传感器用作电容传感器获取水相含率ε<sub>w</sub>与油相含率ε<sub>o</sub>2)利用水相含率ε<sub>w</sub>与油相含率ε<sub>o</sub>计算两相流混合密度ρ<sub>m</sub>=ε<sub>w</sub>ρ<sub>w</sub>+ε<sub>o</sub>ρ<sub>o</sub>、粘度μ<sub>m</sub>=ε<sub>w</sub>μ<sub>w</sub>+ε<sub>o</sub>μ<sub>o</sub>与声速c<sub>m</sub>=c<sub>w</sub>ε<sub>w</sub>+c<sub>o</sub>ε<sub>o</sub>,其中ρ<sub>w</sub>与ρ<sub>o</sub>分别为水相密度与油相密度,μ<sub>w</sub>与μ<sub>o</sub>分别为水相粘度与油相粘度,c<sub>w</sub>与c<sub>o</sub>分别为水中声速与油中声速;3)确定超声多普勒探头的测量空间,该测量空间为三维立体结构,其在管道截面方向投影为椭圆形,在管道剖面方向投影为近似菱形的四边形,分散相流经该区域时,其速度信息可以被超声接收探头有效获取,引入水力学直径的概念将其等效为一球形区域,计算测量空间的等效高度<img file="FDA0000721887210000011.GIF" wi="283" he="148" />其中,测量空间理论高度<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>H</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>Len</mi><mi>sin</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>+</mo><mi>&alpha;</mi><mo>)</mo></mrow></mrow><mrow><mi>sin</mi><mrow><mo>(</mo><mi>&pi;</mi><mo>-</mo><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000012.GIF" wi="696" he="158" /></maths>宽度<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>W</mi><mo>=</mo><mi>d</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><mfrac><mrow><mn>4</mn><mi>R</mi><mo>+</mo><mi>d</mi></mrow><mrow><mn>4</mn><mi>sin</mi><mi>&theta;</mi></mrow></mfrac><mo>-</mo><mi>L</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000013.GIF" wi="539" he="146" /></maths>长度<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Len</mi><mo>=</mo><mfrac><mi>d</mi><mrow><mi>sin</mi><mi>&theta;</mi></mrow></mfrac><mo>+</mo><mfrac><mi>d</mi><mrow><mn>2</mn><mi>tan</mi><mi>&theta;</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mo>+</mo><mi>d</mi><mo>-</mo><mn>2</mn><mi>L</mi><mi>sin</mi><mi>&theta;</mi></mrow><mrow><mn>2</mn><mi>tan</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><mrow><mi>R</mi><mo>-</mo><mi>L</mi><mi>sin</mi><mi>&theta;</mi></mrow><mrow><mi>tan</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>+</mo><mi>&alpha;</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000014.GIF" wi="1074" he="146" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mfrac><msup><mi>d</mi><mn>2</mn></msup><mrow><mn>4</mn><mi>&lambda;</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000721887210000015.GIF" wi="167" he="141" /></maths>为超声近场区域长度,d为超声多普勒探头直径,λ为超声波波长,R为管道内半径,α=sin<sup>‑1</sup>(1.22λ/d)为超声发射探头的发射角;4)计算超声接收探头获取的超声信号频率f与超声发射探头所发射的激励超声信号频率f<sub>0</sub>间的频率差f<sub>d</sub>=f‑f<sub>0</sub>,并计算测量空间中的平均流速<img file="FDA0000721887210000016.GIF" wi="384" he="128" />其中<img file="FDA0000721887210000017.GIF" wi="430" he="199" />为平均频率差,S<sub>d</sub>(f<sub>d</sub>)为超声频移的功率谱;5)利用两相流混合密度ρ<sub>m</sub>与混合粘度μ<sub>m</sub>,建立管道中与距离管壁位置z<sub>i</sub>处流体流速u<sub>i</sub>的分布关系<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>i</mi><mo>+</mo></msubsup><mo>=</mo><mn>2.5</mn><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>z</mi><mi>i</mi></msub><msub><mi>&rho;</mi><mi>m</mi></msub><msub><mi>u</mi><mi>&tau;</mi></msub></mrow><msub><mi>&mu;</mi><mi>m</mi></msub></mfrac><mfrac><mrow><mn>1.5</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mi>r</mi><mi>R</mi></mfrac><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mn>2</mn><msup><mrow><mo>(</mo><mfrac><mi>r</mi><mi>R</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><mn>5.5</mn><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000021.GIF" wi="744" he="295" /></maths>其中<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>i</mi><mo>+</mo></msubsup><mo>=</mo><mfrac><mover><msub><mi>u</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><msub><mi>u</mi><mi>&tau;</mi></msub></mfrac></mrow>]]></math><img file="FDA0000721887210000022.GIF" wi="150" he="156" /></maths>为无量纲流速,且<img file="FDA0000721887210000023.GIF" wi="59" he="85" />为位置z<sub>i</sub>处的平均流速,u<sub>τ</sub>为边界流速,r为位置z<sub>i</sub>距离管道中心的距离且r=R‑z<sub>i</sub>;6)根据5)中的流体流速u<sub>i</sub>的分布关系与测量空间等效高度H,建立测量空间中多普勒测量流速与边界层流速u<sub>τ</sub>间关系<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>u</mi><mi>dop</mi></msub><mo>=</mo><msub><mi>u</mi><mi>&tau;</mi></msub><msub><mi>f</mi><mi>H</mi></msub><mo>+</mo><mn>2.5</mn><msub><mi>u</mi><mi>&tau;</mi></msub><mi>ln</mi><mfrac><mrow><mn>1.5</mn><msub><mi>u</mi><mi>&tau;</mi></msub><msub><mi>&rho;</mi><mi>m</mi></msub><mi>R</mi></mrow><msub><mi>&mu;</mi><mi>m</mi></msub></mfrac><mo>+</mo><mn>5.5</mn><msub><mi>u</mi><mi>&tau;</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000024.GIF" wi="856" he="141" /></maths>其中,系数<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>f</mi><mi>H</mi></msub><mo>=</mo><mn>2.5</mn><mo>[</mo><mi>ln</mi><mfrac><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>-</mo><msup><mi>H</mi><mn>2</mn></msup></mrow><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mi>H</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><msup><mi>R</mi><mn>2</mn></msup><msup><mi>H</mi><mn>2</mn></msup></mfrac><mi>ln</mi><mfrac><msup><mi>R</mi><mn>2</mn></msup><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>-</mo><msup><mi>H</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><msup><mi>R</mi><mn>2</mn></msup><mrow><mn>2</mn><msup><mi>H</mi><mn>2</mn></msup></mrow></mfrac><mi>ln</mi><mfrac><msup><mi>R</mi><mn>2</mn></msup><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mi>H</mi><mn>2</mn></msup></mrow></mfrac><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA0000721887210000025.GIF" wi="1263" he="161" /></maths>7)对ln u<sub>τ</sub>进行级数展开<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>ln</mi><msub><mi>u</mi><mi>&tau;</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>u</mi><mi>&tau;</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mo>&infin;</mo></msubsup><mfrac><mrow><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>m</mi></msup><msup><mrow><mo>(</mo><msub><mi>u</mi><mi>&tau;</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>m</mi></msup></mrow><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>&ap;</mo><msub><mi>u</mi><mi>&tau;</mi></msub><mo>-</mo><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000721887210000026.GIF" wi="946" he="142" /></maths>m为级数展开的阶数,参数D为ln u<sub>τ</sub>级数展开的高次项之和,液液两相流下的一般取值范围为2‑5之间,其他应用条件下可通过实验数据以及级数展开方法推算8)将ln u<sub>τ</sub>=u<sub>τ</sub>‑D代入6)关系式中,求解计算边界层流速<img file="FDA0000721887210000027.GIF" wi="454" he="154" />其中系数a=B,b=B+B ln K‑B*D,c=‑u<sub>dop</sub>,A=f<sub>H</sub>+5.5,B=2.5,<img file="FDA0000721887210000028.GIF" wi="292" he="147" />9)通过两相流混合密度ρ<sub>m</sub>、粘度μ<sub>m</sub>以及边界层流速u<sub>τ</sub>,计算总流速<img file="FDA0000721887210000029.GIF" wi="541" he="400" />n为不同流态下的摩擦指数,当水为连续相时n=‑0.25,当油为连续相时n=‑0.2;将水相含率ε<sub>w</sub>与油相含率ε<sub>o</sub>分别代入得到水相流速J<sub>w</sub>=J*ε<sub>w</sub>和油相流速J<sub>o</sub>=J*ε<sub>o</sub>。
地址 300072 天津市南开区卫津路92号