耿艷峰 華陳權(quán) 王微微 邢蘭昌

(中國石油大學(xué)(華東)信息與控制工程學(xué)院，山東 青島 266580)

0 引言

凝析天然氣也稱濕氣，是指在一定氣井井口溫度、壓力等條件下，體積含液率小于10%的氣井產(chǎn)出物，其液相成分以C5以下的烷烴類組分為主，是非常寶貴的化工原料?，F(xiàn)有的天然氣計量以計量分離器為主，部分海上氣田應(yīng)用了凝析天然氣流量計，這兩種方式都能達到5% ～10%的生產(chǎn)計量要求，但無法滿足產(chǎn)品交接、氣藏管理等精度要求較高的應(yīng)用需求［1－2］。

隨著海上資源的深入開發(fā)，凝析天然氣的不分離在線計量技術(shù)研究引起了國內(nèi)外研究人員的普遍關(guān)注［3］。采用凝析天然氣流量計代替昂貴、笨重的計量分離器，可簡化海洋油氣生產(chǎn)平臺和陸上油氣輸送工藝設(shè)計，降低生產(chǎn)與建設(shè)成本，從而提高氣田開發(fā)的經(jīng)濟效益。

1 技術(shù)現(xiàn)狀與發(fā)展動態(tài)分析

1.1 濕氣計量技術(shù)現(xiàn)狀

基于節(jié)流元件具有的簡單實用、維護成本低等特點，研究人員采用了各種形式的節(jié)流元件，如槽式孔板、多孔孔板、文丘里管及其改進型、V-Cone等作為流量傳感器，輔以必要的信號處理與流動模型技術(shù)，實現(xiàn)了凝析天然氣的在線不分離計量。

目前，基于節(jié)流元件的凝析天然氣計量技術(shù)可以達到±10%以內(nèi)的氣液相計量精度，基本滿足生產(chǎn)計量的精度要求。但由于節(jié)流元件的輸出受流體溫度、壓力、密度、氣液兩相流流型等參數(shù)的影響，很難進一步提高計量精度。

科氏流量計是目前精度較高的單相質(zhì)量流量計，其計量精度遠高于節(jié)流元件［4-6］。其氣體流量測量精度優(yōu)于0.35%，且基本不受流體物性、溫度、壓力與管道截面流速分布等因素的影響。特別是科氏流量計1000∶1的量程比能夠很好地適應(yīng)氣液兩相流動的瞬態(tài)質(zhì)量流量波動，為氣液兩相流的高精度計量提供了可能性。

1.2 科氏流量計與兩相流

Baker R C和Martin A分別對1994年前后的科氏流量計的研究進展與技術(shù)進步進行了全面綜述［5－6］。圍繞流量計本身的性能，國內(nèi)的太航儀表廠、中國測試技術(shù)研究院、北京航空航天大學(xué)、合肥工業(yè)大學(xué)、中國科技大學(xué)、西北工業(yè)大學(xué)等單位先后開展了大量的研究工作，為我國科氏流量計的研制、生產(chǎn)及應(yīng)用奠定了堅實的基礎(chǔ)［4，7-8］。國外的 Micro Motion、E+H、Krohne、Foxboro等公司則先后實現(xiàn)了科氏流量計的商品化，英國 Oxford、Cranfield，德國 Erlangen-Nuerenburg等大學(xué)為科氏流量計原理作出了很多原創(chuàng)性的貢獻［5，9］。

開展利用科氏流量計測量氣液兩相流研究的國內(nèi)外機構(gòu)有:美國 Micro Motion［9］、Foxboro 公司［10］，英國Oxford［11］、Cranfield 大學(xué)［12］，瑞士 E+H 公司，以及我國的浙江大學(xué)［13］、西安交通大學(xué)等。由于受到二次儀表的驅(qū)動與控制能力、信號檢測靈敏度與信號處理算法等方面的限制，早期認為科氏流量計只能適用于含氣率小于10%的場合，國內(nèi)外絕大多數(shù)研究集中在油水兩相流、低含氣稠油等被測介質(zhì)。本文嘗試將科氏流量計應(yīng)用于含氣率大于90%的凝析天然氣計量技術(shù)研究。

在理論研究方面，Cranfield的Hemp J提出了科氏流量計應(yīng)用于氣液兩相流測量的第一個理論模型Bubble Model［12］。E+H的Alfred R等進一步將其發(fā)展成為Moving Resonator Model。GiDRA的 Gysling提出了Aeroelastic Model，用于說明由于氣相存在引起的科氏流量計誤差與兩相流體的聲速變化規(guī)律［14］。隨后Hemp J和Goeke J分別從理論分析和試驗研究兩方面對科氏流量計的氣液兩相流測量特性進行了深入討論［15］。上述理論模型較好地預(yù)報了含氣率小于10%的氣液兩相流測量誤差，但由于流動結(jié)構(gòu)、相間作用等方面的顯著差別，科氏流量計仍無法應(yīng)用于含氣率大于90%的凝析天然氣計量。

在實際應(yīng)用方面，Micro Motion與Foxboro公司分別將科氏流量計成功地應(yīng)用于油水兩相流、低含氣稠油的計量;Foxboro公司正在研究將其科氏流量計應(yīng)用于凝析天然氣計量［16］。相比較而言，科氏流量計在油水兩相流、含懸浮物的固液兩相流測量中效果較好。在氣液兩相流測量中，流量計輸出信號波動很大，有時甚至無法正常工作。為此，絕大多數(shù)研究采用了基于數(shù)據(jù)驅(qū)動的神經(jīng)網(wǎng)絡(luò)、數(shù)據(jù)融合等智能算法，以達到較高的測量精度。

現(xiàn)有技術(shù)存在的問題主要有以下幾個方面。

①盡管應(yīng)用了智能算法的網(wǎng)絡(luò)集成、奇異點處理等技術(shù)，仍無法解決其推廣性不足的缺陷。

②由于商品化的科氏流量計適用對象是單相流測量，氣液兩相流引起的內(nèi)部信號波動經(jīng)流量計廠家的單相流數(shù)據(jù)處理算法平滑和濾波后，不可能通過智能算法得到完全補償。

③流量計輸出信號的波動既是由氣液兩相流動的不穩(wěn)定性造成的，也是由單相流信號處理算法處理兩相流測量數(shù)據(jù)造成的，智能算法無法從根本上解決該問題。

綜上所述，只有在充分了解氣液兩相流流經(jīng)科氏流量計的流動規(guī)律以后，才能建立科氏流量計的氣液相流量、激振電流等輸入?yún)?shù)與振動頻率、幅度、阻尼、相位差等輸出參數(shù)之間的理論模型，進而獲得科氏流量計的氣液兩相流測量模型和更為準確的流量計信號處理算法與控制方法，從而利用科氏流量計實現(xiàn)氣液兩相流高精度計量的目的。

2 凝析天然氣的測量模型

為簡化理論分析，同時不失一般性，以單直管科氏流量計為例，根據(jù)參考文獻［17］的研究結(jié)果，本文提出科氏流量計氣液兩相流測量模型［17］，其單相流測量模型如圖1所示。

圖1 科氏流量計的單相流測量模型Fig.1 The model of Coriolis meter for single-phase flow

圖1將科氏流量計等效為三個連接在一起的二階欠阻尼振動系統(tǒng)，三個位置分別對應(yīng)流量計的一個激勵源和兩個檢測點的位置。若僅考慮振動系統(tǒng)的一階模態(tài)，可以由圖1(a)中的cm、km、mm來代替整個科氏流量計的振動狀態(tài)，因為除了振幅不同以外，圖1(a)中三個振動系統(tǒng)的相位、頻率均相同。振動系統(tǒng)的二階模態(tài)可以等效為左右兩個繞中心點的小幅旋轉(zhuǎn)振動，這兩個振動的幅值與頻率相同，但相位相反。這樣，科氏流量計的工作狀態(tài)就可以等效為兩個檢測點處的一階振動模態(tài)與相位相反的二階振動模態(tài)的疊加。

在單相流測量條件下，振動管內(nèi)的流體被認為是不可壓縮的，管內(nèi)被測流體與振動管振動速度相等，故兩個等效的質(zhì)量塊ml=mr，且均正比于振動管與管內(nèi)流體質(zhì)量之和，ki、ci(i=l，r)分別表示振動管扭轉(zhuǎn)彈性系數(shù)與機械阻尼系數(shù)。在正常工作狀態(tài)下，流量計的機械阻尼由電磁激振力完全補償，從而實現(xiàn)理想的等幅振動與高精度流量測量。

在凝析天然氣測量條件下，水平安裝的振動管內(nèi)氣液兩相流的流型如圖2所示，分別是分層光滑流、分層波浪流、段塞流與環(huán)狀流四種流型。

圖2 凝析天然氣的常見流型Fig.2 Flow regimes of wet gas

由于氣相的存在，科氏流量計的測量模型不再等效為圖1(b)的簡單形式。首先可能出現(xiàn)ml≠mr的現(xiàn)象;其次是等效的ci、ki不再是常數(shù)，其數(shù)值大小受氣液兩相流相關(guān)參數(shù)的影響;第三，由于氣體的可壓縮性，氣液相的振動速度與振動管的振動速度之間出現(xiàn)了差異，也會影響等效科氏力與阻尼系數(shù)的大小。

參考流體力學(xué)中 Sloshing Dynamics的內(nèi)容［18］，本文提出了科氏流量計測量凝析天然氣模型，其示意圖如圖3所示。

圖3 凝析天然氣測量模型Fig.3 The model of wet gas flow measurement

首先，將振動管的振動與被測流體的振動分開考慮，兩者之間通過一個二階振動系統(tǒng)進行耦合;其次將不同的凝析天然氣流型等效為不同的振動參數(shù)。為簡化起見，圖3(a)只畫出了圖1(b)左邊振動系統(tǒng)的等效圖，右邊類同。圖3(a)中，c1、k1、m1分別表示振動管本身的阻尼系統(tǒng)、彈性系數(shù)和質(zhì)量;c2、k2、m2分別表示由管內(nèi)兩相流體等效的阻尼、彈性系數(shù)和兩相質(zhì)量之和;f表示流量計的激振電磁力;x1、x2分別表示振動管與管內(nèi)流體的位移。流量計的振動方程如下:

圖3(a)是一個二自由度振動系統(tǒng)，振動參數(shù)除了與振動管本身的結(jié)構(gòu)、材料有關(guān)以外，還與凝析天然氣流型、氣液相截面含率、管道壓力、流體黏度等流動參數(shù)有關(guān)。其中，由兩相流體與振動管的相對運動引起的阻尼系數(shù)c2和由氣相的可壓縮性引起的彈性系數(shù)k2無法忽略。

在正常工作條件下，流量計兩個振動系統(tǒng)的阻尼系數(shù)由激振電磁力完全補償，振動管工作在等幅振動狀態(tài)，流量計振動模型轉(zhuǎn)變?yōu)槿鐖D3(b)所示的模型。此時的振動方程式由(1)簡化為:

式(2)仍然是一個二自由度振動系統(tǒng)，對式(2)進行求解，即可得到振動系統(tǒng)的兩個共振頻率

設(shè)兩個共振頻率ω1＜ω2。當流量計正常工作時，ω1小于相同壓力的單相氣體共振頻率，但大于單相液體共振頻率，其表示振動管與兩相流體作同方向運動，ω1是流量計的工作頻率;ω2大于相同壓力的單相氣體的共振頻率，表示振動管與兩相流體之間作相對運動。

當c2→0、k2→∞時，即忽略流體與振動管管壁的相對運動、假設(shè)管內(nèi)流體是不可壓縮流體，上述模型就是文獻［17］給出的科氏流量計單相流測量模型。

3 試驗研究

3.1 試驗裝置簡介

試驗研究工作是在英國Cranfield大學(xué)的多相流試驗裝置上進行的，試驗裝置的組成和基本功能介紹參見文獻［19］。

試驗測試數(shù)據(jù)分布如圖4所示。

圖4 試驗數(shù)據(jù)分布Fig.4 The experimental data set

試驗選用了兩臺E+H公司的 Promass 83F微彎科氏流量計(研究工作獲得了E+H公司的技術(shù)與資金支持)，采集并記錄了流量計內(nèi)部的振動阻尼、激振電流等信號。試驗管路直徑為50.8 mm，試驗的氣相流量為200 ～500 m3·h－1;液相流量為 0.11 ～1.5 kg·s－1;壓力為100～300 kPa;折算的氣相體積含率為88%～99%，屬于凝析天然氣的參數(shù)范圍。兩臺質(zhì)量流量計順序串聯(lián)，分別為垂直安裝和水平安裝。

3.2 試驗結(jié)果分析

為便于敘述，以下稱垂直安裝的科氏流量計為CMF01，稱水平安裝的科氏流量計為CMF02。除專門對比兩種安裝方式的測量結(jié)果以外，主要介紹CMF02的凝析天然氣測量特性。

①流量計參數(shù)的波動特性

氣液兩相流測量條件下，科氏流量計參數(shù)波動很大。流量計輸出頻率(fCMF01、fCMF02)、激振電流(ICMF01、ICMF02)和阻尼系數(shù)(CCMF01、CCMF02)三個參數(shù)在三種不同測量條件下的實測值如表1所示。測試過程中，單相氣體流量為150～700 m3·h－1;壓力為300 kPa、溫度為室溫;單相液體流量為 0.2 ～7.1kg·s－1;氣液兩相流由350 m3·h－1的空氣流量和 0.45 kg·s－1的水流量組成，氣相含率為97.7%。

表1 測量方式Tab.1 Data comparision of flow meter

表1中的數(shù)據(jù)說明:科氏流量計可以在凝析天然氣測量條件下正常工作，但流量計輸出參數(shù)與內(nèi)部參數(shù)波動較大。

相關(guān)數(shù)據(jù)表明，兩相流測量時振動頻率的波動約大于氣相測量的四倍，約為液相測量的100倍，這就必然要求對流量計的采集與控制電路、控制算法與數(shù)據(jù)處理算法作出相應(yīng)的改進。

②水平安裝與垂直安裝的測量結(jié)果比較

在單相流測量條件下，流量計水平安裝和垂直安裝區(qū)別不大。對于氣液兩相流計量，以上兩種安裝方式有很大差別。

由于垂直管道中氣液兩相的截面含率分布均勻，流型種類少，使得密度測量值相對穩(wěn)定可靠，因此，利用伽馬射線、電容/電阻等傳感器進行含率測量的多相流量計更傾向于垂直安裝;利用節(jié)流元件的多相流測量系統(tǒng)傾向于水平安裝，回避垂直安裝時重力壓降與節(jié)流壓降的難以區(qū)分的問題。

兩臺流量計不同安裝形式的測量結(jié)果如圖5所示。其中:圖5(a)是科氏流量計用于單相液體的測量結(jié)果，測量過程中液體流量從 25 kg·s－1變化到45 kg·s－1;圖5(b)是凝析天然氣的測量結(jié)果，氣相流量為 200 m3·h－1，液相流量為 0.11kg·s－1。

由圖5可以看出，在單相流測量中，兩種安裝方式的測量結(jié)果基本相同。由于CMF01安裝在CMF02前，圖5(a)中CMF01的流量變化超前于 CMF02，流量計輸出信號波動是由于單相流管線中混入了少量的氣體;圖5(b)的兩相流測量中，水平安裝的流量計CMF02輸出數(shù)據(jù)大于垂直安裝的 CMF01，且其波動較大。

圖5 兩種安裝形式的測量結(jié)果對比Fig.5 Comparison results of two installation methods

③兩相混合物密度測量結(jié)果

流量計輸出的氣液兩相流體密度測量值與科氏流量計的振動頻率有關(guān)。

在不同氣相流量下，CMF02的密度輸出結(jié)果與由單相氣液流量換算得到的兩相真實密度之間的關(guān)系如圖6所示。

由圖6可知，在氣相流量一定的情況下，真實密度與流量計讀數(shù)之間近似呈線性關(guān)系。

圖6 流量計輸出密度與真實密度Fig.6 Density output and true density

④兩相混合物質(zhì)量流量的測量結(jié)果

在整個測量參數(shù)范圍內(nèi)，流量計輸出的混合物質(zhì)量流量與由氣液單相計量結(jié)果換算所得的兩相真實質(zhì)量流量之間的關(guān)系如圖7所示。

圖7 輸出質(zhì)量流量與真實質(zhì)量流量Fig.7 Mass flowrate output and true mass flowrate

由圖7可知，流量計輸出值總體上低于兩相真實質(zhì)量流量的理論計算值，但兩者之間近似呈線性關(guān)系。這種簡單的線性關(guān)系說明，對流量計輸出結(jié)果進行校正是可行的，并有望達到較高的校正精度，對此需要進一步研究影響線性關(guān)系斜率的主要因素。

4 結(jié)束語

理論分析與試驗研究均表明，利用科氏流量計實現(xiàn)凝析天然氣的高精度計量既是可能的，也是可行的。本文提出的科氏流量計二自由度振動模型綜合考慮了氣液兩相流流型和科氏流量計的工作原理［20－23］，為流量計工作特性與控制方法研究提供了依據(jù);氣液兩相的質(zhì)量流量測量值與真實值之間存在較好的線性關(guān)系，氣液兩相的密度測量值能夠線性地跟蹤真實值的變化，可簡化測量結(jié)果的修正算法，奠定了流量計實際應(yīng)用基礎(chǔ)。

［1］江懷友，趙文智，張東曉.世界天然氣資源及勘探現(xiàn)狀研究［J］.天然氣工業(yè)，2008，28(7):21 －28.

［2］Bieker P，Slupphaug O.Real-time production optimization of oil and gas production systems:a technology survey［J］.SPE Production ＆Operations，2007，22(4):382 －391.

［3］Geng Yanfeng，Zheng Jinwu，Shi Tianming.Study on the metering characteristics of slotted orifice for wet gas flow［J］.FlowMeasurement and Instrumentation，2006，17(2):123 －128.

［4］樊尚春，張軒.科氏流量計的非線性影響機理分析與抑制［J］.科學(xué)技術(shù)與工程，2005，5(9):569 －571.

［5］Martin A，Wolfgang D.Coriolis mass flow meters:overview of current state of the art and latest research［J］.Flow Measurement and Instrumentation，2006，17(6):317 －323.

［6］Weinstein A，Kassoy R.Experimental study of oscillatory motion of particles and bubbles with applications to Coriolis flow meters［J］.Physics of Fluids，2008，20(10):306 －315.

［7］倪偉，徐科軍.科里奧利流量計數(shù)字信號處理方法研究［D］.合肥:合肥工業(yè)大學(xué)，2004.

［8］仝猛，陳明.科氏流量計理論與應(yīng)用研究［D］.西安:西北工業(yè)大學(xué)，2003.

［9］Baker R.Coriolis flow meters:industrial practice and published information［J］.Flow Measurement and Instrumentation，1994，5(4):229－246.

［10］Seeger M.Coriolis flow measurement in two phase flow［J］.Computing and Control Engineering，2005，16(3):10 －16.

［11］Henry M，Tombs M.Two-phase flow metering of heavy oil using a Coriolis meter:a case study［J］.Flow Measurement and Instrumentation，2006，17(6):399－413.

［12］Hemp J，Kutin J.Theory of errors in Coriolis flow meter reading due to compressibility of the fluid being metered［J］.Flow Measurement and Instrumentation，2006，17(6):359 －369.

［13］馬龍博，張宏建，周洪亮，等.基于Coriolis流量計和SVM的油水兩相流質(zhì)量流量測量的研究［J］.高校化學(xué)工程學(xué)報，2007，21(2):200－205.

［14］Gysling L.An aeroelastic model of Coriolis mass and density meters operating on aerated mixtures［J］.Flow Measurementand Instrumentation，2007，18(2):69 －77.

［15］Goeke J，Stock M.New parameters for analysis of distortions of a Coriolis mass flow meter［J］.Technisches Messen，2007，74(11):577 －588.

［16］Lansangan R，Skinner J.Coriolis mass flow metering for wet gas［J］.Measurement and Control，2008，41(7):213 －216.

［17］Storm R.Coriolis mass flowmeter and method for operating a Coriolis mass flow meters［P］.US 7136777，2006.

［18］Takashi I，Shin M.Auto parametric resonances in a structure/fluid interaction system carrying a cylindrical liquid tank［J］.Journal of Sound and Vibration，2005(285):517 －546.

［19］Cao Y，Hoi Y.The instruction of multiphase flow loop in Cranfield University［D］.Shrivenham:Cranfield University，2004.

［20］Chisholm D.Two-phase Flow Through Sharp-Edged Orifices［J］.Journal of Mechanical Engineering Science，1977，19(3):128 －130.

［21］Lin Z H.Two-phase Flow Measurement with Sharp-edged Orifices［J］.International Journal of Multiphase Flow，1982，8(6):683 －693.

［22］Huang Zhiyao，Xie Dailiang，Zhang Hongjian，et al.Gas-oil twophase flow measurement using an ecectrical capacitance tomography system and a Venturi meter［J］.Flow Measurement and Instrument，2005，16(2):177 －182.

［23］林宗虎.氣液固多相流測量［M］.北京:中國計量出版社，1988:82－86.