諧振腔微擾法測(cè)量原油含水率仿真研究

摘 要： 與傳統(tǒng)在線測(cè)量原油含水率的方法相比，高頻微波諧振腔微擾法具有受外界環(huán)境影響小，不受礦化度影響，測(cè)量精度高的特點(diǎn)。由于原油和水的介電常數(shù)相差較大，微小的含水率變化會(huì)引起介電常數(shù)較大變化，使得不同含水率原油對(duì)微波場(chǎng)的損耗不同。通過(guò)軟件仿真的方法可以測(cè)量出含微擾源（樣品原油）的微波諧振腔的諧振頻率和品質(zhì)因數(shù)變化情況，根據(jù)矩形諧振腔微擾公式可得出樣品原油的介電常數(shù)，然后將所得結(jié)果代入原油含水率和介電常數(shù)關(guān)系公式求得原油含水率。并且將鋁制邊界條件和銅制邊界條件下所測(cè)得的原油介電常數(shù)數(shù)據(jù)與理論值進(jìn)行分析，求出其相對(duì)百分比誤差，得出一種較為簡(jiǎn)單、實(shí)用、快速準(zhǔn)確的測(cè)量中低含水率原油的方法。

關(guān)鍵詞： 諧振腔； 諧振頻率； 品質(zhì)因數(shù)； 含水率； 微擾法

中圖分類號(hào)： TN98?34 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2016）20?0004?04

Abstract： In comparison with the traditional on?line detecting method of moisture content of crude oil， the high?frequency microwave resonant cavity perturbation method is seldom influenced by the external environment and never influenced by the mineralization， and has high detecting accuracy. Since the permittivity of water and crude oil is significantly different， and the small change of water content can change the permittivity greatly， the crude oil with different moisture contents may result in the different losses of microwave field. The change of resonant frequency and quality factor of the microwave resonant cavity containing perturbation source （sample crude oil） are detected by means of software simulation. The permittivity of sample crude oil can be obtained according the perturbation formula of rectangle resonant cavity. And then the obtained result is substituted into the relation formula reflecting crude oil moisture content and permittivity to get the permittivity of the sample moisture content of crude oil. The crude oil permittivity data detected in aluminum and copper boundary conditions is analyzed with the theoretical value， and the relative percentage error is found out. A simple， practical， fast and accurate method to detect the crude oil with medium and low moisture content is obtained.

Keywords： resonant cavity； resonant frequency； quality factor； water content； perturbation method

0 引 言

隨著社會(huì)經(jīng)濟(jì)對(duì)石油能源的需求量逐年提升，工業(yè)設(shè)備和石化采集、運(yùn)輸系統(tǒng)對(duì)原油含水率信息的要求越來(lái)越精細(xì)，使得原油含水率的測(cè)量技術(shù)直接影響著石油工業(yè)的發(fā)展，主要表現(xiàn)在原油的開(kāi)采、脫水、集輸計(jì)量、儲(chǔ)運(yùn)銷售和石化煉制等重要環(huán)節(jié)，而精確及時(shí)的測(cè)量含水率十分困難，所以尋找一種及時(shí)準(zhǔn)確的含水率測(cè)量方法顯得尤為重要。及時(shí)準(zhǔn)確地獲得含水率信息不僅能夠?yàn)榉€(wěn)定地控制石油生產(chǎn)的過(guò)程提供理論依據(jù)，而且能夠保證成品原油的生產(chǎn)，進(jìn)而降低生產(chǎn)成本，減少能源消耗。傳統(tǒng)的測(cè)量方法（如電容法、電導(dǎo)法等）只能滿足較低準(zhǔn)確度，即使能夠達(dá)到較高準(zhǔn)確度也需要耗費(fèi)大量的時(shí)間和人力（如電化學(xué)方法中的卡爾?費(fèi)休法、蒸餾法和電脫法等）。在低含水（<30%）時(shí)，水的礦化度對(duì)于水的介電常數(shù)影響較小，當(dāng)對(duì)高含水（>60%）時(shí)，水的礦化度會(huì)較大地影響電容式測(cè)量、電導(dǎo)式等在線測(cè)量方法的精度。據(jù)此提出了一種測(cè)量精確度高、測(cè)量時(shí)間短、并且不受原油礦化度影響的方法，即微波諧振腔微擾法測(cè)量原油含水率[1?2]。微波具有穿透能力強(qiáng)，不僅能夠測(cè)量樣品原油表層的水分，還可以檢測(cè)樣品原油內(nèi)部（油包水）的水分，并且具有測(cè)量時(shí)間短、效率高、不受礦化度影響等顯著優(yōu)點(diǎn)[3?4]。諧振腔微擾法根據(jù)引入微擾源中的油、水介電常數(shù)相差較大，即水分子電場(chǎng)極化強(qiáng)度遠(yuǎn)遠(yuǎn)高于油。原油樣品中含水量的多少對(duì)于電磁場(chǎng)的幅值損耗和相位延遲產(chǎn)生影響，使得高頻電磁波在諧振腔內(nèi)產(chǎn)生根據(jù)含水量變化的高頻振蕩諧振，諧振腔中的諧振頻率和品質(zhì)因數(shù)會(huì)隨著發(fā)生變化。根據(jù)諧振腔內(nèi)電磁諧振頻率和品質(zhì)因數(shù)可以得出樣品原油的介電常數(shù)的變化，然后由原油含水率和介電常數(shù)的關(guān)系公式可求得樣品原油的含水率。

1 原油含水率與介電常數(shù)關(guān)系

介電常數(shù)又稱為電容率，表示介質(zhì)在外電場(chǎng)作用下極化程度的物理量，反映了介質(zhì)存儲(chǔ)電荷的能力，是相對(duì)介電常數(shù)和真空絕對(duì)介電常數(shù)的乘積，常溫常壓下，一般水的相對(duì)介電常數(shù)為80左右，而純油（如原油、柴油和煤油等）相對(duì)介電常數(shù)均為2.3左右。油水混合物的相對(duì)介電常數(shù)與其含水率密切相關(guān)，水分子作為極性分子，而油分子作為非極性分子，由于極性分子與非極性分子在電極化場(chǎng)的作用下極化機(jī)理不同，因此，油水混合物的相對(duì)介電常數(shù)值不等于它們的算術(shù)平均值。

在實(shí)際生產(chǎn)中原油和水的混合存在三種取向極化狀態(tài)[5]：

本文為了符合實(shí)際情況，仿真設(shè)計(jì)中采取第三種取向極化狀態(tài)的表達(dá)式，根據(jù)極化式（3）可求得不同含水率原油的混合介電常數(shù)的變化情況。

2 微擾法原理

微波諧振腔主要分為圓柱形諧振腔和矩形諧振腔。由于矩形諧振腔有著廣泛的理論基礎(chǔ)和測(cè)試研究，有較為準(zhǔn)確的測(cè)試表達(dá)公式，并且具有品質(zhì)因數(shù)高和傳輸損耗小的特點(diǎn)，所以在此只討論矩形諧振腔的微擾法[6?8]，其結(jié)構(gòu)如圖1所示。 由于水的介電常數(shù)遠(yuǎn)遠(yuǎn)大于原油的介電常數(shù)，所含水分對(duì)微波場(chǎng)的損耗比純凈原油大的多，將含水原油樣品放到諧振腔的適當(dāng)位置，微波通過(guò)諧振腔中的油體樣品時(shí)，其諧振頻率和品質(zhì)因數(shù)會(huì)發(fā)生相應(yīng)變化，對(duì)于原油含水率測(cè)量來(lái)說(shuō)，主要關(guān)心的是由于含水量不同所引起的介電常數(shù)變化，所以將樣品放在諧振腔微波電場(chǎng)最大，而微波磁場(chǎng)為零處，當(dāng)樣品[dh?110] 很小時(shí)，可以將[f-f0] 看成微擾[9]。

式中：[f0]為諧振腔諧振頻率；[f]為放入樣品后諧振頻率；[Δ1QL]為放入樣品后諧振腔的有載品質(zhì)因數(shù)[QL] 的倒數(shù)變化；[μ]為樣品的張量磁導(dǎo)率；[I]為單位張量；[εr]為樣品的相對(duì)介電常數(shù)；[E0]和[H0]為空腔中的電磁場(chǎng)；[E]和[H]為放入樣品后腔中的電磁場(chǎng)；[VS]為樣品的體積；[V0]為諧振腔體積（[VS?V0]）。

為了分析簡(jiǎn)便，將樣品放入諧振腔電場(chǎng)最大而磁場(chǎng)為零處時(shí)，仿真結(jié)果和理論相一致，即矩形諧振腔在本征模式條件下存在兩處電場(chǎng)最強(qiáng)而磁場(chǎng)強(qiáng)度為零的點(diǎn)，[x=a2，y=l4]，空腔仿真圖形，如圖2所示。

由式（7）可知，通過(guò)測(cè)量諧振腔諧振頻率和品質(zhì)因數(shù)的變化，可以將復(fù)介電常數(shù)中的[ε′]和[ε″]求出[8]。其中[ε″]是復(fù)介電常數(shù)中表示電磁波在諧振腔中介質(zhì)的電極化損耗，而[ε′]是介電常數(shù)部分[9]。

3 仿真及結(jié)果分析

本文為了方便計(jì)算和仿真，采用國(guó)標(biāo)BJ?9型矩形諧振腔，其起始頻率為0.76 GHz，截止頻率為1.15 GHz，內(nèi)截面尺寸[247.6 mm×123.8 mm，]外截面尺寸[253.6 mm×129.8 mm，]內(nèi)部圓角[r=1.2 mm，]諧振腔長(zhǎng)度選擇為[400 mm]，厚度為[3 mm]。由于本文只對(duì)腔體中所存在的電磁場(chǎng)結(jié)構(gòu)和對(duì)應(yīng)的頻率值進(jìn)行數(shù)值上的求解，所以采用諧振腔本征模的數(shù)值求解方法。其設(shè)置最小求解頻率為0.8 GHz，最大迭代次數(shù)3次，收斂誤差10%，其求解模式數(shù)有5個(gè)，其中模式1為最低次模。仿真計(jì)算所得到的頻率具有虛部和實(shí)部，其中實(shí)部是諧振頻率，而諧振頻率的虛部和各種損耗有關(guān)，在此不考慮損耗，所以忽略虛部的影響。

將空腔的面電場(chǎng)圖和放入純油的面電場(chǎng)圖進(jìn)行比較發(fā)現(xiàn)，矩形諧振腔中電場(chǎng)強(qiáng)度最大的兩個(gè)點(diǎn)的值減少了約0.027 V/m，滿足微擾的條件要求，如圖3所示。

原油樣品含水率取12組油水兩相不同配比的樣品，純油、1%含水率的油、10%含水率的油、20%含水率的油、30%含水率的油、40%含水率的油、50%含水率的油、50%含水率的油、60%含水率的油、70%含水率的油、80%含水率的油、90%含水率的油和純水。矩形諧振腔內(nèi)電場(chǎng)強(qiáng)度隨著微擾樣品中原油含水率的增加而減少，如圖4所示，其諧振腔的諧振頻率也依次降低，如圖5所示，能夠較好地反應(yīng)出原油含水率的變化情況。

相對(duì)介電常數(shù)為第三種既有串聯(lián)又有并聯(lián)的形式，經(jīng)計(jì)算所得理論介電常數(shù)和仿真計(jì)算所得介電常數(shù)對(duì)比如表1所示。

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

根據(jù)所求得的諧振頻率代入式（7）中，可以算出所測(cè)油水混合物的介電常數(shù)，對(duì)于本文仿真分析來(lái)說(shuō)，取原油樣品充分混合，對(duì)于樣品即油水混合物相對(duì)介電常數(shù)為第三種既有串聯(lián)又有并聯(lián)的形式取原油樣品充分混合，對(duì)于樣品即油水混合物將兩組數(shù)據(jù)進(jìn)行數(shù)據(jù)處理并畫(huà)出圖形進(jìn)行對(duì)比，如圖6所示，當(dāng)含水率低于70%時(shí)，仿真測(cè)量相對(duì)介電常數(shù)和理論值非常接近，說(shuō)明仿真的準(zhǔn)確性；當(dāng)含水率高于70%時(shí)，測(cè)量值略低于理論計(jì)算值，仿真出現(xiàn)較大誤差。

如圖7所示，當(dāng)含水率小于70%的時(shí)候，樣品原油介電常數(shù)測(cè)量值和理論真實(shí)值的相對(duì)百分比誤差小于5%，和真實(shí)值較為接近，可以滿足實(shí)際測(cè)量的需要（銅制邊界條件下的測(cè)量結(jié)果和鋁制邊界結(jié)果基本相同，故只列出鋁制邊界條件下的仿真結(jié)果）。

5 結(jié) 語(yǔ)

本文提出一種測(cè)量精確度高、測(cè)量時(shí)間短、并且不受原油礦化度影響的方法，即微波諧振腔微擾法測(cè)量原油含水率，并進(jìn)行仿真，結(jié)果表明該方法可以滿足實(shí)際測(cè)量的需要。

參考文獻(xiàn)

[1] 張亞萍，閆向宏.水礦化度及有機(jī)物含量對(duì)原油含水率的影響[J].油氣儲(chǔ)運(yùn)，2002（2）：32?34.

[2] 黃正華，陳世廉.水的礦化度對(duì)測(cè)量原油含水率的影響[J].中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1993（2）：116?120.

[3] 李志茂，黃志堯，王保良，等.微波透射法測(cè)量含水率的研究[J].傳感技術(shù)學(xué)報(bào)，2006，19（6）：2544?2546.

[4] 李志茂.基于微波透射法測(cè)量油水兩相流分相含率的實(shí)驗(yàn)研究[D].杭州：浙江大學(xué)，2006.

[5] 黃正華.油水混合介質(zhì)相對(duì)介電常數(shù)的研究[J].油氣田地面工程，2000，19（2）：8?9.

[6] 張夢(mèng)君.小型高Q值介質(zhì)填充諧振腔微擾法測(cè)量微波復(fù)介電常數(shù)研究[D].南京：南京郵電大學(xué)，2013.

[7] 劉澄宇.微波諧振腔微擾法測(cè)量植物葉片介電常數(shù)[J].長(zhǎng)春師范學(xué)院學(xué)報(bào)，2007，26（6）：44?46.

[8] 唐宗熙，張其劭.改進(jìn)矩形腔微擾法測(cè)試微波電介質(zhì)復(fù)介電常數(shù)[J].通信學(xué)報(bào)，1996（5）：56?61.

[9] 王建華，丁言鎂，左春英.微波諧振腔微擾法測(cè)量油體含水量的研究[J].沈陽(yáng)師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，24（4）：433?435.

[10] 陳波，楊德強(qiáng)，潘錦.同軸諧振腔及微擾法測(cè)量介質(zhì)介電常數(shù)[J].實(shí)驗(yàn)科學(xué)與技術(shù)，2009（4）：4?6.

[11] 張緒峰，劉宇，張婷.微波諧振腔諧振頻率實(shí)時(shí)測(cè)量系統(tǒng)的實(shí)現(xiàn)[J].傳感技術(shù)學(xué)報(bào)，2009，22（10）：1451?1454.