透射式油水兩相流超聲檢測(cè)系統(tǒng)設(shè)計(jì)與測(cè)試

張濰，譚超

?

透射式油水兩相流超聲檢測(cè)系統(tǒng)設(shè)計(jì)與測(cè)試

張濰，譚超

(天津大學(xué)電氣與自動(dòng)化工程學(xué)院，天津市過(guò)程檢測(cè)與控制重點(diǎn)實(shí)驗(yàn)室，天津，300072)

為實(shí)現(xiàn)油水兩相流過(guò)程信息的提取和參數(shù)的檢測(cè)，設(shè)計(jì)一種基于超聲透射模式的測(cè)試系統(tǒng)。該系統(tǒng)可實(shí)現(xiàn)超聲信號(hào)發(fā)生、發(fā)射、多通道測(cè)試數(shù)據(jù)采集、測(cè)試數(shù)據(jù)處理、存儲(chǔ)及流動(dòng)信息實(shí)時(shí)顯示等功能。對(duì)不同含油率的油水兩相進(jìn)行流測(cè)試。研究結(jié)果表明，該系統(tǒng)能夠有效獲取油水兩相流相含率的變化信號(hào)。超聲透射測(cè)試信號(hào)衰減率與相含率、混相速度相關(guān)：在不同水流量條件下，超聲透射信號(hào)隨含油率增加呈現(xiàn)先增加后減小的趨勢(shì)；在水包油和油包水過(guò)渡流型中，超聲透射測(cè)試信號(hào)衰減率與混相速度成線性關(guān)系。

油水兩相流；超聲檢測(cè)；透射式；系統(tǒng)設(shè)計(jì)；相含率；混相速度

多相流過(guò)程廣泛存在于能源、動(dòng)力、石油、化工、冶金、醫(yī)藥等工業(yè)過(guò)程中。多相流問(wèn)題的研究對(duì)促進(jìn)相關(guān)工業(yè)過(guò)程工藝、設(shè)備的發(fā)展和創(chuàng)新有重要作用[1]。多相流現(xiàn)象經(jīng)常出現(xiàn)在管道輸送、流化床、電廠以及反應(yīng)堆等生產(chǎn)過(guò)程中，并與安全問(wèn)題、經(jīng)濟(jì)問(wèn)題息息相關(guān)。對(duì)多相流參數(shù)的精確測(cè)量可以更好地指導(dǎo)實(shí)際工業(yè)過(guò)程操作，但多相流參數(shù)測(cè)量的難度給工業(yè)界和科學(xué)研究提出挑戰(zhàn)[2]。油水兩相流為多相流的一種，是石油工業(yè)中的主要流態(tài)，其流動(dòng)狀態(tài)受油水兩相含率與流速影響呈現(xiàn)多種復(fù)雜的流動(dòng)結(jié)構(gòu)。準(zhǔn)確測(cè)量和估計(jì)油水兩相流的流型、相含率、流量等參數(shù)對(duì)于相關(guān)生產(chǎn)工藝設(shè)計(jì)、流動(dòng)狀態(tài)預(yù)測(cè)、生產(chǎn)設(shè)備維護(hù)等有非常重要的價(jià)值。與氣水兩相流相比，油水兩相介質(zhì)的密度相對(duì)接近，流動(dòng)的波動(dòng)程度相對(duì)較小，因而參數(shù)波動(dòng)范圍較小且容易出現(xiàn)連續(xù)相轉(zhuǎn)換現(xiàn)象，對(duì)檢測(cè)系統(tǒng)與方法設(shè)計(jì)提出了更高要求。常用的兩相流參數(shù)檢測(cè)的方法有：介入式檢測(cè)方法，如文丘里、孔板等差壓式測(cè)量方法[3]；非介入方式的檢測(cè)方法，如微波、射線、超聲等[4]。其中，超聲波在不同的介質(zhì)中具有不同的衰減系數(shù)并具有很好的方向性，可用來(lái)獲取兩相流的分相含率[5]。1955年“鳴環(huán)法”超聲波流量計(jì)在美國(guó)研制成功，這也是世界上第1臺(tái)超聲波流量 計(jì)[6]。20世紀(jì)70年代，隨著電子技術(shù)的高速發(fā)展，超聲波檢測(cè)設(shè)備與電子設(shè)備結(jié)合得到進(jìn)一步發(fā)展[7]。近年來(lái)，超聲方法檢測(cè)多相流參數(shù)獲得了越來(lái)越多的重視，檢測(cè)方法主要有反射式、透射式和Doppler 3種。VOLKER等[8]應(yīng)用反射模式檢測(cè)懸濁液并獲得了成功。IHARA等[9]應(yīng)用超聲透射方法檢測(cè)泡狀流紊流條件下速度分布，進(jìn)而解釋和證實(shí)氣泡產(chǎn)生的摩擦阻力減少現(xiàn)象。WILLATZEN等[10]解決了透射聲衰減模式下的聲波傳播非線性問(wèn)題。CAMILLI等[11]利用超聲Doppler方法進(jìn)行了相含率測(cè)量。IHARA等[12]利用Doppler方法測(cè)量頻率變化，對(duì)低流速流體的流速分布做出了描述。隨著數(shù)據(jù)融合技術(shù)的發(fā)展，超聲方法與其他方法進(jìn)行融合，可以得到更好的測(cè)量結(jié)果[13]。近年來(lái)，超聲測(cè)試方法在氣液兩相流動(dòng)過(guò)程的研究中發(fā)展較快，也為油水兩相流的測(cè)量提供了一種解決方 案[14?16]。本文作者利用透射模式超聲衰減效應(yīng)構(gòu)建了一套兩相流相含率超聲檢測(cè)系統(tǒng)，在實(shí)驗(yàn)裝置上進(jìn)行油水兩相流動(dòng)態(tài)實(shí)驗(yàn)對(duì)該系統(tǒng)進(jìn)行了測(cè)試，并利用相含率范圍、兩相流流型以及多重散射理論對(duì)測(cè)試結(jié)果進(jìn)行分析。

1 超聲波多相流檢測(cè)基本原理

在超聲波檢測(cè)技術(shù)中，根據(jù)超聲波在不同媒質(zhì)中的反射、折射、衍射和散射等傳播特性，通過(guò)對(duì)相應(yīng)傳播特性信息的獲取和分析，即可實(shí)現(xiàn)對(duì)被測(cè)物場(chǎng)的特性的檢測(cè)。由于聲波能量與頻率平方成正比，能很好地探測(cè)到不透明物體的內(nèi)部構(gòu)造。

透射模式超聲檢測(cè)方法利用聲波的反射和折射規(guī)律，設(shè)為兩相介質(zhì)分界面處聲阻抗之比：

式中：1和2為相鄰界面的兩相聲阻抗。因此，透射模式的聲壓衰減率T、聲強(qiáng)衰減率T和反射模式的聲壓衰減率R、聲強(qiáng)衰減率R為：

在20℃、垂直入射時(shí)，油水聲壓反射率ow=14.73%，油水聲壓透射率ow=85.27%。在油水兩相流中，利用超聲波在穿透油、水介質(zhì)時(shí)，被吸收能量的不同，通過(guò)測(cè)量聲壓在被測(cè)兩相介質(zhì)場(chǎng)中的衰減，可得到油水介質(zhì)的分布、相含率和流速等信息。

在油水兩相流檢測(cè)過(guò)程中，定義如下：

1) 表觀流速：

式中：o和w為油相和水相的表觀流速；o和w為油相和水相的入口流量；為管道截面積。

2) 相含率：

3) 混相速度：

4) 信號(hào)幅值衰減率：

式中：w為純水時(shí)超聲透射測(cè)試信號(hào)幅值；o/w為油水兩相流時(shí)超聲透射測(cè)試信號(hào)幅值。

2 透射式超聲檢測(cè)系統(tǒng)設(shè)計(jì)

2.1 超聲檢測(cè)系統(tǒng)結(jié)構(gòu)

超聲檢測(cè)系統(tǒng)主要包括5個(gè)主要部分：超聲波傳感器陣列模塊、激勵(lì)信號(hào)發(fā)生模塊、信號(hào)轉(zhuǎn)換及超聲波發(fā)射(或接收)模塊，信號(hào)采集模塊和計(jì)算機(jī)處理系統(tǒng)，其結(jié)構(gòu)如圖1所示。其中，信號(hào)發(fā)生模塊產(chǎn)生一定頻率的激勵(lì)信號(hào)；信號(hào)轉(zhuǎn)換及超聲波波發(fā)射(或接收)模塊完成信號(hào)的濾波、放大、輸出發(fā)射(或接收)等功能；信號(hào)采集和處理模塊，對(duì)接收的信號(hào)進(jìn)行A/D轉(zhuǎn)換和數(shù)據(jù)的預(yù)處理等功能，并將數(shù)據(jù)傳入計(jì)算機(jī)進(jìn)行存儲(chǔ)、計(jì)算和顯示。

圖1 檢測(cè)系統(tǒng)結(jié)構(gòu)圖

2.2 超聲波傳感器陣列

超聲波傳感器選用中心頻率為1 MHz，直徑9 mm的收發(fā)一體超聲波傳感器，分布方式如圖1中超聲傳感器陣列部分所示。超聲波傳感器垂直安裝在管道外壁，相鄰傳感器的夾角相等。每個(gè)傳感器以一定的時(shí)間間隔輪流被觸發(fā)，依次發(fā)射超聲信號(hào)。該信號(hào)經(jīng)過(guò)物場(chǎng)后，由相對(duì)的傳感器探頭進(jìn)行接收。

2.3 激勵(lì)信號(hào)發(fā)生模塊

信號(hào)發(fā)生模塊結(jié)構(gòu)如圖2所示，包括電源轉(zhuǎn)換單元、FPGA(field-programmable gate array)、DDS(direct digital synthesizer)芯片以及濾波放大環(huán)節(jié)。電源轉(zhuǎn)換模塊為VCVS單元，該模塊將輸入電壓轉(zhuǎn)換為15 V電源給FPGA供電和3.3 V電源給DDS芯片供電。FPGA作為一個(gè)控制單元，為DDS芯片提供控制字，調(diào)節(jié)信號(hào)頻率。選用時(shí)鐘頻率為300 MHz的DDS模塊生成1 MHz的正弦信號(hào)。

圖2 超聲波發(fā)生和發(fā)射模塊結(jié)構(gòu)圖

2.4 信號(hào)轉(zhuǎn)換及超聲波收發(fā)模塊

濾波放大環(huán)節(jié)采用四階Butterworth濾波器將DDS生成的正弦信號(hào)進(jìn)行濾波,該環(huán)節(jié)選用Butterworth高通濾波器，截止頻率為800 kHz。放大器將信號(hào)放大到±10 V，激勵(lì)超聲波傳感器。

2.5 信號(hào)采集模塊

超聲信號(hào)經(jīng)過(guò)管道中的流動(dòng)物場(chǎng)，每遇到兩相分界面即產(chǎn)生透射和折射過(guò)程，因此會(huì)產(chǎn)生諸多頻率與原始頻率不相同的噪聲。為減少該噪聲與系統(tǒng)本身噪聲，本模塊中先對(duì)信號(hào)進(jìn)行了濾波處理。信號(hào)采集模塊完成超聲信號(hào)的濾波、數(shù)模轉(zhuǎn)換、數(shù)據(jù)采集等功能。信號(hào)采集模塊結(jié)構(gòu)如圖3所示。

圖3 信號(hào)采集模塊結(jié)構(gòu)圖

2.6 計(jì)算機(jī)處理系統(tǒng)

計(jì)算機(jī)處理系統(tǒng)包括信號(hào)采集模塊的驅(qū)動(dòng)、數(shù)據(jù)采集程序、顯示、存儲(chǔ)等。數(shù)據(jù)采集程序中完成采樣參數(shù)設(shè)置、數(shù)據(jù)存儲(chǔ)管理和系統(tǒng)設(shè)置等功能，流程如圖4所示。數(shù)據(jù)顯示將采集到的多路信號(hào)進(jìn)行實(shí)時(shí)波形顯示。在數(shù)據(jù)存儲(chǔ)管理中，通過(guò)設(shè)置一個(gè)可變的地址，將采集到的數(shù)據(jù)以.txt格式存儲(chǔ)在該地址上。系統(tǒng)設(shè)置完成對(duì)系統(tǒng)的保護(hù)和總體控制，該環(huán)節(jié)創(chuàng)建了參數(shù)可變的采樣數(shù)窗口，當(dāng)實(shí)際采集達(dá)到該參數(shù)時(shí)，程序即跳出循環(huán)停止運(yùn)行。

圖4 數(shù)據(jù)采集程序流程圖

數(shù)據(jù)采集系統(tǒng)軟件程序界面由多個(gè)選項(xiàng)組成，包括通道選擇與采樣率控制的設(shè)置窗口；用于觀察觀察波形是否正確的實(shí)時(shí)波形顯示窗口；用于采集前設(shè)定數(shù)據(jù)存儲(chǔ)位置的保存地址設(shè)置窗口；另外還包含系統(tǒng)啟動(dòng)運(yùn)行開關(guān)按鈕，系統(tǒng)狀態(tài)指示信號(hào)顯示等。

3 油水兩相流實(shí)驗(yàn)及結(jié)果分析

3.1 實(shí)驗(yàn)過(guò)程

在水平管路中，針對(duì)不同油水比的兩相流進(jìn)行實(shí)驗(yàn)測(cè)試。實(shí)驗(yàn)管路為水平放置的有機(jī)玻璃透明管道，便于觀察油水兩相流動(dòng)狀態(tài)，管道內(nèi)徑為50 mm。

實(shí)驗(yàn)中改變兩相流流量，調(diào)節(jié)水量這3~8 m3/h，每增加1 m3/h測(cè)1次，共計(jì)6次；調(diào)節(jié)油量為0~14 m3/h，每增加1 m3/h測(cè)1次，共15次；兩兩匹配共90個(gè)測(cè)量值。超聲激勵(lì)信號(hào)為正弦波，幅值10 V、頻率1 MHz。測(cè)試信號(hào)的采樣頻率為10 MHz。實(shí)驗(yàn)測(cè)試時(shí)，每個(gè)測(cè)試點(diǎn)油水配比穩(wěn)定流動(dòng)持續(xù)2 min，每次測(cè)試采樣時(shí)間為20 s。實(shí)驗(yàn)條件對(duì)應(yīng)的流型如圖5所示。

圖5實(shí)驗(yàn)條件對(duì)應(yīng)流型圖

3.2 實(shí)驗(yàn)結(jié)果

圖6所示為水流量為3 m3/h、含油率分別為0和50%時(shí)，超聲傳感器接收的電壓信號(hào)。從圖6可見：當(dāng)油相加入時(shí)，采集所得到的信號(hào)波形、頻率與超聲激勵(lì)信號(hào)一致、沒有發(fā)生變化，但信號(hào)幅值發(fā)生了一定程度的衰減。

含油率/%：1—0；2—50。

對(duì)每個(gè)采樣時(shí)間內(nèi)的超聲透射幅值信號(hào)取正弦波峰值，然后取平均，得到該采樣時(shí)間內(nèi)的平均幅值，不同含油率下的超聲電壓幅值信號(hào)如圖7所示。

水流量/(m3?h?1)：1—3；2—4；3—5；4—6；5—7；6—8。

由式(2)和(3)可得：在超聲波傳播過(guò)程中，多重散射是超聲測(cè)試信號(hào)幅值衰減的主要原因，且衰減量與超聲波在兩相介質(zhì)傳播過(guò)程中遇到兩相分界面數(shù)量近似成正比，也即與流速和含率相關(guān)。對(duì)于不同含率條件的流動(dòng)狀態(tài)，存在2種不同衰減形式：1) 當(dāng)含油率低于30%時(shí)，位于流型圖中油包水流型區(qū)域，此時(shí)水為連續(xù)相、油為離散相，流動(dòng)狀態(tài)為水包油；當(dāng)含油率高于70%時(shí)，位于流型圖中水包油流型區(qū)域，此時(shí)油為連續(xù)相、水為離散相流動(dòng)狀態(tài)為油包水[17]。這2種含率條件下，兩相流的油水界面比較清晰、且分布比較規(guī)律，超聲測(cè)試信號(hào)的衰減主要是透射衰減。超聲測(cè)試信號(hào)遇到分界面的概率隨離散相含率升高而升高，測(cè)試信號(hào)的衰減率是隨著離散相含率增加而增大。2) 當(dāng)含油率在30%~70%時(shí)，位于流型圖中油包水和水包油流型區(qū)域，該區(qū)域?yàn)樗团c油包水流型交替出現(xiàn)的過(guò)渡區(qū)域，測(cè)試信號(hào)的衰減主要與多重散射相關(guān)[18]。隨兩相流混相流速的不同，油水兩相流出現(xiàn)水包油與油包水流型之間的過(guò)渡，此時(shí)油、水作為離散相交替出現(xiàn)，油水界面分布復(fù)雜。當(dāng)流體速度變大時(shí)，油水液滴分散更均勻，超聲測(cè)試信號(hào)在液滴間發(fā)生多重散射，導(dǎo)致信號(hào)的衰減。

在水包油與油包水流型過(guò)渡區(qū)間內(nèi)，相同相含率條件下，超聲波測(cè)試信號(hào)衰減率與混相速度的關(guān)系，如圖8和圖9所示。由圖8可知：超聲波測(cè)試信號(hào)衰減率與含油率變化趨勢(shì)成反比，同一含油率的超聲衰減隨水流量的增加而減小。由圖9可知：超聲波測(cè)試信號(hào)衰減率隨著油水兩相流混相速度增加而減小，呈近似線性關(guān)系。

上述實(shí)驗(yàn)結(jié)果說(shuō)明，測(cè)試信號(hào)的衰減與油水兩相流含率和速度的變化同時(shí)相關(guān)。當(dāng)含油率低于30%或高于70%時(shí)，在管道內(nèi)形成水包油和油包水流型，測(cè)試信號(hào)的衰減以透射衰減為主，衰減率隨著離散相的含率增加而增大；當(dāng)含油率為30%~70%時(shí)，管道內(nèi)形成水包油與油包水流型的過(guò)渡時(shí)，測(cè)試信號(hào)的衰減以多重散射衰減為主，需通過(guò)結(jié)合流動(dòng)狀態(tài)和流速的變化，進(jìn)一步建立超聲測(cè)試的模型，實(shí)現(xiàn)油水兩相流相含率等參數(shù)的估計(jì)和表征。

4 結(jié)論

1) 基于透射模式的超聲波測(cè)試信號(hào)衰減率與油水兩相流的相含率、混相速度相關(guān)；在不同水流量條件下，測(cè)試信號(hào)衰減隨油相含率增加呈現(xiàn)先增加后減小的趨勢(shì)。

2) 當(dāng)含油率低于30%(水包油流型)或高于70%(油包水流型)時(shí)，超聲波測(cè)試信號(hào)衰減率與相含率和混相速度呈線性關(guān)系；當(dāng)含油率為30%~70%時(shí)，超聲波衰減以多重散射效應(yīng)為主，其衰減率同時(shí)受相含率與流速影響。

3) 該系統(tǒng)可廣泛應(yīng)用于油水兩相流工業(yè)測(cè)量中。由于在特定含率條件下，超聲波測(cè)試信號(hào)衰減與油含率呈現(xiàn)線性趨勢(shì)，可以實(shí)時(shí)估測(cè)管道中油的含率，亦可為進(jìn)一步實(shí)現(xiàn)油水兩相流參數(shù)的準(zhǔn)確獲取，提供測(cè)試工具和方法基礎(chǔ)。

[1] 林宗虎, 王棟, 王樹眾, 等. 多相流的近期工程應(yīng)用趨向[J]. 西安交通大學(xué)學(xué)報(bào), 2001, 35(9): 886?890. LIN Zonghu, WANG Dong, WANG Shuzhong, et al. Recent trend towards engineering application of multi-phase flow[J]. Journal of Xi’an Jiaotong University, 2001, 35(9): 886?890.

[2] 譚超, 董峰. 多相流過(guò)程參數(shù)檢測(cè)技術(shù)綜述[J]. 自動(dòng)化學(xué)報(bào), 2013, 39(11):1923?1932.TAN Chao, DONG Feng. Parameters measurement for multiphase flow process[J]. Acta Automatica Sinica, 2013, 39(11):1923?1932.

[3] THORN R, JOHANSEN G A, HJERTAKER B T. Three-phase flow measurement in the petroleum industry[J]. Measurement Science & Technology, 2012, 24(1):12003?12019.

[4] THORN R, JOHANSEN, G A, HAMMER E A. Recent developments in three-phase measurement[J]. Measurement Science and Technology, 1997, 8(7): 691?701.

[5] ZHAO An, HAN Yunfeng, REN Yingyu, et al. Ultrasonic method for measuring water holdup of low velocity and high-water-cut oil-water two-phase flow[J]. Applied Geophysics, 2016, 13(1):179?193.

[6] 蘇茜, 董峰. 油水兩相流超聲波衰減測(cè)試方法[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 47(2):647?653.SU Qian, DONG Feng. Ultrasound attenuation measurement in oil-water two-phase flow[J]. Journal of Central South University (Science and Technology), 2016, 47(2):647?653.

[7] DONG Xiaoxiao, TAN Chao, YE Yuan, et al. Measuring oil–water two-phase flow velocity with continuous-wave ultrasound doppler sensor and drift-flux model[J]. IEEE Transactions on Instrumentation & Measurement, 2016, 65(5):1098?1107.

[8] VOLKER A, de KROON M. Characterization of suspensions by ultrasonic reflection measurements[J]. Ultrasonics, 1998, 36(1): 283?289(7).

[9] MURAI Y, FUJI H, TASAKA Y. Turbulent bubbly channel flow investigated by ultrasonic velocity profiler[J]. Journal of Fluid Science and Technology, 2006, 1(1): 12?23.

[10] WILLATZEN M, KAMATH H. Nonlinearities in ultrasonic flow measurement[J]. Flow Measurement and Instrumentation, 2007, 19(2): 79?84.

[11] CAMILLI R, DI I D, BOWEN A, et al. Acoustic measurement of the deepwater horizon macondo well flow rate[J]. Proceedings of the National Academy of the Sciences of the United States of America(PNAS), 2012, 109(50): 20235?20239.

[12] IHARA T, KIKURA H, TAKEDA Y. Ultrasonic velocity profiler for very low velocity field[J]. Flow Measurement and Instrumentation, 2013, 34(6): 127?133.

[13] ZHAI Lusheng, JIN Ningde, GAO Zhongke. The ultrasonic measurement of high water volume fraction in dispersed oil-in-water flows[J]. Chemical Engineering Science, 2013, 94(5): 271?283.

[14] CHAUDHURI A, SINHA D N, ZALTE A. Mass fraction measurements in controlled oil-water flows using noninvasive ultrasonic sensors[J]. Journal of Fluids Engineering, 2014, 136(3): 1?8.

[15] COUTINHO F R, OFUCHI C Y, ARRUDA. L V R. A new method for ultrasound detection of interfacial position in gas-liquid two-phase flow[J]. Sensors, 2014, 14(5): 9093?9116.

[16] ABBASZADEH J, RAHIM H B A, RAHIM R B A. Frequency analysis of ultrasonic wave propagation on metal pipe in ultrasonic tomography system[J]. Sensor Review, 2014, 34(1): 13?23.

[17] TRALLERO J L, SARICA C, BRILL J P. A study of oil/water flow patterns in horizontal pipes[J]. SPE Production and Facilities, 1997, 12(3): 165?172.

[18] MCCLEMENTS D J, POVEY M J W. Scattering of ultrasound by emulsions[J]. Journal of Physics D: Applied Physics, 1989, 22(1): 38?47.

Design and testing of measurement system for oil-water two-phase flow based on ultrasonic transmission mode

ZHANG Wei, TAN Chao

(Tianjin Key Laboratory of Process Measurement and Control, School of Electrical Engineering and Automation, Tianjin University, Tianjin 300072, China)

A transmission mode acoustic system was designed for oil-water two phase flow measurement. The system can be used for signal generation, signal emission, multi-channel data acquisition, data procession, data saving and data visualization and etc. The experiment of oil-water two-phase flow with different oil contents was carried out. The results show that the system can successfully acquire changing signal of oil-water two-phase flow. By analyzing the experiment data, the attenuation of the amplitude of transmission mode ultrasound test voltage is found to have relationship with phase volume fraction and phase velocity. For different quantities of water, ultrasound test voltages increase first and then decrease with the increase of oil contents. In water in oil and oil in water flow pattern transition, the attenuation has a linear relationship with mixed phase velocity.

oil-water two-phase flow; ultrasonic testing; transmission mode; system design; phase volume fraction； mixed phase velocity

10.11817/j.issn.1672-7207.2016.09.046

TP216，TH814+.92

A

1672?7207(2016)09?3252?06

2015?05?15；

2015?09?25

國(guó)家自然科學(xué)基金資助項(xiàng)目(61227006，61473206)；天津科技創(chuàng)新體系及平臺(tái)建設(shè)計(jì)劃資助項(xiàng)目(13TXSYJC40200) (Projects(61227006, 61473206) supported by the National Natural Science Foundation of China? Project(13TXSYJC40200) supported by the Science and Technology Innovation Plan of Tianjin)

譚超，博士，副教授，從事多相流測(cè)試技術(shù)與系統(tǒng)研究；E-mail: tanchao@tju.edu.cn

(編輯 趙俊)