楊樹坤， 楊 佩， 李 翔， 郭宏峰

(1. 中海油田服務股份有限公司油田生產(chǎn)事業(yè)部，天津 300459; 2. 西安石油大學電子工程學院，陜西 西安 710065)

0 引 言

隨著社會的不斷發(fā)展，目前石油工業(yè)已經(jīng)成為國民的支柱產(chǎn)業(yè)，并帶動了機械、交通運輸、鋼鐵、化工等行業(yè)的不斷發(fā)展。然而，石油的開采是一道極其復雜的過程，現(xiàn)階段石油的開采主要采用的是高壓注水的方法，基于石油和水之間不同的性質(zhì)，使石油漂浮起來，這種開采方式導致開采出來的原油中含水率較高。因此，在進行石油開采之前應對每個油井的含水率進行檢測，如果產(chǎn)出原油的價值低于開發(fā)成本時，應當停止開發(fā)或者改進開采方案降低成本?？梢?，對油井含水率這一指標進行實時監(jiān)測具有重要的意義。

目前國內(nèi)外含水率的測量方法有很多種，主要分為兩大類[1]：離線測量法和在線測量法。常見的離線測量法[2]有：蒸餾法、離心法、電脫法和卡爾-費休法等；此類方法采用人工取樣方式，測試結(jié)果準確，但缺點是每次采集的樣品隨機性大且不及時，不能有效反應原油含水率的實況。在線測量法能夠?qū)崿F(xiàn)對油井含水率實時而有效的監(jiān)控，常用的在線測量法有：射線法[3-7]、電容法[8-13]、電導法[14]、射頻法[15-17]、紅外光譜法[18]以及微波法[19-26]等。這些方法仍可以被分為兩大類：接觸式測量法和非接觸式測量方法，其中電容法、電導法、射頻法以及微波法中的傳輸法和諧振腔法均屬于接觸式測量法，而射線法、紅外光譜法和微波透射法屬于非接觸式測量法。與接觸式測量法相比，非接觸式測量法不容易受流態(tài)、結(jié)垢、結(jié)蠟的影響，容易實現(xiàn)無損耗、實時、在線快速測量。射線法雖然能夠?qū)崿F(xiàn)較寬范圍內(nèi)的含水率測試，但存在輻射源，使用不安全且價格昂貴；而紅外光譜法屬于新技術(shù)，還處于發(fā)展階段。與之相比，微波透射法技術(shù)成熟，對人體幾乎沒有危害，操作和使用非常安全，因此具有很好的應用前景。

對于微波透射法來說，其微波傳感器被設(shè)置在管道的外壁處，不需要占用管道的中心通道且不與流體接觸，避免了微波傳感器粘油、結(jié)臘等可能造成的測量失效問題。微波透射法的核心是其中的微波傳感器，其性能的優(yōu)劣決定了天線的測試精度和量程。在實際使用中最小的油管的內(nèi)徑僅為40.30 mm，因此微波傳感器的小型化設(shè)計尤為重要。在現(xiàn)有的微波透射法研究中，所應用的微波傳感器主要有喇叭天線[23]、金屬絲螺旋天線[24]、微帶螺旋天線[25]、微帶貼片天線[26]等。如表1所示可以看出，在頻率較低的情況下，微波傳感器輻射口徑已經(jīng)與最小油管的內(nèi)徑基本相等或大于最小油管尺寸，因此上述研究僅能應用于內(nèi)徑較大的油管中；另外在頻率較高情況下，雖然天線的尺寸結(jié)構(gòu)可以減小，但是對微波透射法含水率測量儀的其他模塊提出了較高的要求，而且工作頻率越高器件造價越昂貴，成本較高。

表1 現(xiàn)有研究中微波透射法所使用的微波傳感器

基于以上基礎(chǔ)，考慮到微波傳感器小型化這一技術(shù)難題，一款Vivaldi天線被設(shè)計作為微波傳感器。Vivaldi天線的輻射口徑面遠遠小于油管的內(nèi)徑，同時具有工作頻帶寬、輻射增益穩(wěn)定等優(yōu)點，非常符合微波透射法含水率測量系統(tǒng)的設(shè)計要求。因此，在油管實際尺寸的約束條件下，由于Vivaldi天線的輻射口徑寬度遠小于油管內(nèi)徑，因此不會發(fā)生天線體積過大而超出油管內(nèi)徑限制的問題。最后通過CST仿真和矢量網(wǎng)絡(luò)分析儀測試，研究如何使用Vivaldi天線作為微波傳感器實現(xiàn)原油含水率的測試。

1 天線設(shè)計

1.1 Vivaldi天線的理論分析與設(shè)計

Vivaldi天線的結(jié)構(gòu)如圖1所示，該天線在介質(zhì)基板兩側(cè)的輻射貼片均采用相同的結(jié)構(gòu)，包括一條微帶線和一個漸變槽線結(jié)構(gòu)，并在末端進行阻抗匹配后通過50 Ω的同軸接口進行激勵。該天線被印制在相對介電常數(shù)為εr= 4.40、損耗角正切為tanδ =0.02且厚度為Hs= 2 mm 的FR-4介質(zhì)基板上。如圖1（a）所示，天線內(nèi)指數(shù)漸變槽線和外指數(shù)漸變槽線可以通過式（1）和（2）表示為

圖1 Vivaldi天線的詳細結(jié)構(gòu)圖

式中：(xi, yi)——內(nèi)指數(shù)漸變槽線的坐標；

C1i，C2i，F(xiàn)i，Si——內(nèi)指數(shù)漸變槽線表達式（1）的系數(shù)；

(xo, yo)——外指數(shù)漸變槽線的坐標；

C1o，C2o，F(xiàn)o，So——內(nèi)指數(shù)漸變槽線表達式（2）的系數(shù)。

通過采用CST仿真軟件對天線的指數(shù)槽線和微帶線的相應參數(shù)進行優(yōu)化，可以得出整個天線的各參數(shù)尺寸如表2所示。

表2 Vivaldi天線的結(jié)構(gòu)尺寸

1.2 Vivaldi天線的仿真結(jié)果

對于一個天線來說，其最重要的性能指標為反射系數(shù)S11、輻射方向圖、增益以及工作帶寬，其中反射系數(shù)S11是衡量天線阻抗匹配和工作狀態(tài)優(yōu)劣的重要指標，通常認為如果反射系數(shù)S11低于-10 dB則天線處于完美的工作狀態(tài)。圖2和圖3給出了該Vivaldi天線的輻射性能和帶寬性能。圖2給出了該天線在中心頻率（5 GHz）處的3D輻射方向圖，可以看出該天線實現(xiàn)了完美的前向輻射，其輻射增益為6.70 dB。與此同時由圖3（a）中可以看出，該天線在2～8 GHz的頻帶內(nèi)滿足反射系數(shù)S11低于-10 dB，說明天線在這個頻段范圍內(nèi)阻抗匹配和工作狀態(tài)良好，且其工作帶寬為120%，實現(xiàn)了超寬帶工作。另外還可以看出，隨著頻率的增加天線的輻射增益也隨之增加，如圖3（b）所示，當頻率從2 GHz變化至8 GHz時，增益逐漸增大并伴隨著一些波動。通過以上結(jié)果可以看出，該天線在2～8 GHz具有良好的輻射性能，工作頻帶寬、輻射增益穩(wěn)定。

圖2 Vivaldi天線的輻射方向圖

圖3 Vivaldi天線的帶寬性能

1.3 Vivaldi天線的測試結(jié)果

為實現(xiàn)微波透射法原油含水率測量儀系統(tǒng)，加工了兩個相同的Vivaldi天線作為發(fā)射天線和接收天線，如圖4所示?？紤]到反射系數(shù)S11是衡量天線阻抗匹配和工作狀態(tài)優(yōu)劣的重要指標，先通過矢量網(wǎng)絡(luò)分析儀對兩個天線的反射系數(shù)S11進行測試，測試結(jié)果如圖5所示?？梢钥闯觯@兩個天線在2.5～8 GHz的頻帶內(nèi)基本滿足反射系數(shù)S11低于-10 dB，在5.10～5.40 GHz反射系數(shù)S11出現(xiàn)了微小的抬升，且在5.32 GHz處反射系數(shù)S11達到最大值-7.38 dB。當反射系數(shù)為S11=-7.38 dB時天線輻射的能量能夠被計算為，可見當反射系數(shù)S11等于-7.38 dB，超過82%的能量被輻射出去，這說明反射系數(shù)的微小變化并不影響天線的整體性能，因此所設(shè)計的該Vivaldi天線可用于原油含水率測量儀中作為微波傳感器。

圖4 Vivaldi天線的加工圖

圖5 Vivaldi天線的反射系數(shù)實測

1.4 Vivaldi天線的安裝方式

在微波傳感器的設(shè)計中，這個Vivaldi天線的輻射口徑為70 mm×2 mm，其在油管上的安裝方式如圖6所示?？梢钥闯?，油管在x方向上的尺寸為油管內(nèi)徑的尺寸，在y方向上的尺寸為任意尺寸。由于2 mm的尺寸遠小于油管內(nèi)徑，因此該Vivaldi天線能夠適用于任意尺寸的油管。

圖6 Vivaldi天線的安裝方式

在現(xiàn)有的研究，所設(shè)計的微波傳感器基本上為正方形結(jié)構(gòu)，容易受到油管尺寸的限制，無法自由應用于內(nèi)徑較小的油管上，因此需要在油管之間增加一個內(nèi)徑較大的管子作為測試區(qū)域，這種突變將會導致油管內(nèi)的壓差較大。而Vivaldi天線將輻射口徑中一個維度的尺寸大大縮減，實現(xiàn)了微波傳感器的小型化設(shè)計。

2 微波透射法原油含水率測量儀的設(shè)計

2.1 微波透射法原油含水率測量儀系統(tǒng)

圖7給出了微波透射法原油含水率測量儀的系統(tǒng)框圖。微波透射法測試是根據(jù)不同含水率的原油對微波傳輸能量有不同的衰減和反射，使得接收到的功率發(fā)生變化，再經(jīng)過一系列的信號處理就能得到被測物值的含水率情況。由圖7可以看出，該測量系統(tǒng)的核心是微波傳感器（收、發(fā)天線），天線的性能好壞直接影響著測量的精度和量程?？紤]到含水率的不同將會導致相對介電常數(shù)和損耗角正切的不同，而介電常數(shù)的變化也會導致反射系數(shù)S11的偏移。因此，盡管測量頻率不變，但含水率的變化也會導致反射系數(shù)S11的變化，由于天線的正常工作要求S11低于-10 dB，因此一個寬帶化的天線被需要，保證在含水率變化時天線的正常工作，Vivaldi天線作為一個超寬帶天線剛好符合該系統(tǒng)的設(shè)計需求；與此同時，考慮到油管內(nèi)徑尺寸的限制，天線的小型化設(shè)計尤為重要，由于Vivaldi天線輻射口徑面的寬度僅為介質(zhì)基板的厚度，因此不會發(fā)生天線體積過大而超出油管內(nèi)徑限制的技術(shù)問題。

圖7 微波透射法原油含水率測量儀的系統(tǒng)框圖

2.2 微波透射法原油含水率測量儀的仿真模型

圖8給出了微波透射法原油含水率測量儀在CST中的數(shù)值仿真模型。在本次仿真中，油水混合液處于一個高為150 mm，內(nèi)徑為50.30 mm，管壁厚度為10 mm的油管中。由于油管是金屬材質(zhì)，為了使天線發(fā)射的電磁波穿透油管，在油管兩邊分別開有30 mm × 90 mm的長方形開孔，并加入耐高溫高壓、且透波性好的材料進行填充，避免油水混合液的泄露，然后在兩邊搭載前面設(shè)計的Vivaldi天線作為發(fā)射天線和接收天線。至此，完成了微波透射法原油含水率測量儀的數(shù)值仿真模型建立。

圖8 微波透射法原油含水率測量儀在CST中的仿真模型

對油水混合液而言，含水率不同，其介電常數(shù)也不同。其中，油水混合液的相對介電常數(shù)以及損耗角正切均與含水率有關(guān)，具體函數(shù)關(guān)系如下：

式中： εmix、 εwater和 εoil——油水混合液、純水和純油的相對介電常數(shù)；

(tanδ)mix、 ( tanδ)water和——油水混合液、純水和純油的損耗角正切；

φ——油的含水率，單位%。

根據(jù)上式，可以計算得出油水混合液的相對介電常數(shù)和損耗角正切。在常溫工作狀態(tài)下，純油的相對介電常數(shù)為2.33，損耗角正切為0.000 1，而純水的相對介電常數(shù)為78.00，損耗角正切為0.120 0，因此油水混合液的相對介電常數(shù)和損耗角正切則可以被確定，如表3所示。

表3 油水混合液介質(zhì)參數(shù)與含水率的關(guān)系

3 仿真和實測分析

3.1 仿真分析

在CST求解過程中，設(shè)置求解器為時域求解器，邊界條件為開放邊界，能量收斂為-40 dB、工作頻帶設(shè)置為2～8 GHz，左邊的天線設(shè)置為發(fā)射天線并激勵，右邊的天線設(shè)置為接收天線，研究不同含水率下傳輸系數(shù)S21的變化，其中傳輸系數(shù)S21指的是發(fā)射天線發(fā)出的能量經(jīng)過油水混合液衰減后被接收天線接收到的能量，根據(jù)傳輸系數(shù)S21經(jīng)過不同油水混合液后的衰減值來描述含水率的變化。

圖9給出了傳輸系數(shù)S21隨含水率的變化曲線，可以看出，當天線工作頻率在4.39 GHz、含水率在0～100%的變化范圍內(nèi)時，天線的傳輸系數(shù)S21隨著含水率的增加逐漸減小。當油水混合液的含水率為0%，即為純油時，電磁波在純油中的傳輸系數(shù)S21為-10.58 dB；當原油含水率為100%，即為純水時，電磁波在純水中的傳輸系數(shù)S21為-26.77 dB。根據(jù)仿真結(jié)果，可以發(fā)現(xiàn)曲線基本上呈現(xiàn)線性的變化趨勢。為了更具體的描述0～100%范圍內(nèi)傳輸系數(shù)S21和含水率的關(guān)系，對該變化曲線進行了擬合，擬合函數(shù)為

圖9 S21隨含水率變化的仿真和擬合結(jié)果

式中：x（%）——含水率；

y（dB）——傳輸系數(shù)S21。

從圖9（a）中可以發(fā)現(xiàn)擬合函數(shù)得到的擬合曲線和仿真曲線在高含水率時具有較好的一致性，在低含水率時一致性稍差。將油水混合液的仿真結(jié)果帶入擬合關(guān)系式中，得出相應的含水率，其結(jié)果如表4所示。通過數(shù)據(jù)處理可以看出，當含水率在0～100%范圍內(nèi)變化時，仿真得到的絕對誤差1小于8%。

表4 含水率的仿真結(jié)果及絕對誤差

另外，考慮到6%和8%兩個狀態(tài)下傳輸系數(shù)產(chǎn)生了較為明顯的突變，這是造成絕對誤差增大的主要原因。因此如果不考慮0～6%范圍的數(shù)據(jù)，僅對8%～100%范圍內(nèi)傳輸系數(shù)S21和含水率之間的關(guān)系進行擬合，則可以得到擬合關(guān)系式為

基于上述關(guān)系式，通過數(shù)據(jù)處理可以得出絕對誤差2小于4%。因此可以看出，微波透射法仍存在一定局限性，在低含水率時存在較大的誤差。這是由于此時油水混合液的介電常數(shù)較小，導致該油水混合液的相對波長較大，將會導致微波傳感器處于感應近場區(qū)，而非遠場輻射區(qū)。

3.2 實測分析

在上一節(jié)主要從仿真角度分析了電磁波在不同含水率下的傳輸特性，可以看出在4.39 GHz工作頻率下，實現(xiàn)了0～100%范圍內(nèi)的含水率與傳輸系數(shù)S21之間的線性變化，同時在0～6%的低含水率范圍內(nèi)存在較大的絕對誤差。下面主要通過實驗測量電磁波在不同含水率下的傳輸特性，驗證Vivaldi天線實現(xiàn)含水率測量的可行性。實驗測試系統(tǒng)如圖10所示，被測介質(zhì)被放在直徑為50.30 mm的油管內(nèi)，左邊的天線為發(fā)射天線，右邊的天線為接收天線，通過矢量網(wǎng)絡(luò)分析儀實時分析兩個天線之間傳輸系數(shù)S21的變化，從而獲得含水率與傳輸系數(shù)之間的關(guān)系。

圖10 實驗測試系統(tǒng)

圖11給出了傳輸系數(shù)S21隨含水率的變化曲線，可以看出，當天線工作頻率在4.39 GHz、含水率在0～100%的變化范圍內(nèi)時，天線的傳輸系數(shù)S21隨著含水率的增加逐漸減小，但是在0～8%的低含水率產(chǎn)生了非單調(diào)性變化。下面對該變化曲線進行了擬合，擬合函數(shù)為

圖11 S21隨含水率變化的測試和擬合結(jié)果

將油水混合液的測試結(jié)果帶入擬合關(guān)系式中，得出相應的含水率，其結(jié)果如表5所示。通過數(shù)據(jù)處理可以看出，當含水率在0～100%范圍內(nèi)變化時，實測得到的絕對誤差3小于7%。

表5 含水率的測試結(jié)果及絕對誤差

另外，如果不考慮0～8%范圍的非單調(diào)性變化的數(shù)據(jù)，僅對10%～100%范圍內(nèi)傳輸系數(shù)S21和含水率之間的關(guān)系進行擬合，則可以得到擬合關(guān)系式為

基于上述關(guān)系式，通過數(shù)據(jù)處理可以得出絕對誤差4小于5.20%，且在10%～85%范圍內(nèi)絕對誤差小于5%。因此可以看出，微波透射法高含水率時具有較大誤差，這是由于兩種混合液的體積相差較大，不易實現(xiàn)均勻混合，因此導致測試誤差增大。而對于低含水率，一個原因是混合液體積相差較大，另一個原因則是天線可能工作在感應近場區(qū)。

4 結(jié)束語

在本文中，提出了一種基于Vivaldi天線的微波透射法原油含水率測量儀，通過CST全波仿真分析和實驗測試，得到如下結(jié)論：

1）設(shè)計了一種Vivaldi天線，該天線的輻射口徑為70 mm×2 mm，其2 mm的厚度遠小于油管內(nèi)徑，可應用于任意尺寸的油管作為微波傳感器，同時該天線具有較寬的工作帶寬，非常符合微波透射法含水率的測試需求；

2）將Vivaldi天線作為微波傳感器，在CST全波仿真軟件中進行仿真，可以看出所提出的Vivaldi天線能夠?qū)崿F(xiàn)0～100%范圍內(nèi)的含水率與傳輸系數(shù)S21的線性變化。在0～100%范圍內(nèi)含水率的真實值與擬合值之間的絕對誤差小于8%；而在8%～100%范圍內(nèi)含水率的真實值與擬合值之間的絕對誤差小于4%；

3）搭建微波透射法原油含水率測量系統(tǒng)，通過矢量網(wǎng)絡(luò)分析儀對原油含水率與傳輸系數(shù)S21之間的關(guān)系進行測試，測試結(jié)果表明：Vivaldi天線能夠?qū)崿F(xiàn)0～100%范圍內(nèi)的含水率測試，絕對誤差小于7%；而在10%～85%范圍內(nèi)含水率的真實值與擬合值之間的絕對誤差小于5%。