摘(zhai)摘要:四(si)電極外(wai)流式電(dian)磁流量(liang)計 是一(yi)種新型(xíng)的測量(liang)注入剖(pou)面流量(liang)的測井(jing)理想儀(yí)器,廣泛(fan)應用于(yú)油田注(zhù)水井、注(zhù)聚井的(de)流量測(cè)量。目前(qián)四電極(ji)外流式(shì)電磁流(liu)量計的(de)研究主(zhu)要在實(shi)際環境(jìng)中開展(zhan),實驗效(xiào)率低、成(cheng)本高。建(jian)立了四(si)電極外(wài)流式電(diàn)磁流量(liang)計的準(zhǔn)确的有(you)限元模(mo)型,将強(qiang)耦合的(de)方法應(ying)用在電(diàn)磁結構(gòu)流體耦(ou)合.上,并(bìng)在不同(tong)流速下(xia)開展了(le)模型的(de)響應及(ji)誤差分(fèn)析。研究(jiu)表明,該(gāi)有限元(yuan)模型在(zài)一定程(cheng)度上可(kě)用于電(diàn)磁流量(liàng)計的流(liu)場仿真(zhēn)分析。流(liu)速較小(xiǎo)時,有限(xian)元仿真(zhen)結果與(yǔ)實驗誤(wù)差較大(dà);流速較(jiao)大時,流(liu)場趨近(jin)于勻速(su)場,仿真(zhēn)結果與(yǔ)實驗結(jie)果誤差(cha)較小。
　　在(zai)油田三(sān)次采油(yóu)中,注聚(ju)合物驅(qu)油是提(ti)高原油(yóu)采收率(lü)的重要(yao)手段之(zhī)一,它比(bi)水驅效(xiao)果提高(gāo)了20%左右(yòu);現場實(shi)驗表明(míng),過去常(cháng)用的注(zhù)入剖面(miàn)測井儀(yí)器已經(jīng)不适合(he)注聚合(hé)物測井(jing)的剖面(mian)測試的(de)要求。電(dian)磁流量(liàng)計是一(yī)種新型(xing)的測量(liàng)注入剖(pōu)面的儀(yi)器,較好(hǎo)地解決(jue)了聚合(hé)物注入(rù)剖面的(de)測井問(wen)題。
　　四電(dian)極外流(liú)式電磁(ci)流量計(ji)是針對(duì)油田應(yīng)用開發(fa)的一種(zhǒng)特殊電(diàn)磁流量(liang)計,其不(bu)僅具有(yǒu)普通工(gong)業電磁(cí)流量計(ji)無節流(liu)阻流,不(bu)易堵塞(sāi),耐腐蝕(shí)性好,測(cè)量精度(dù)不受被(bei)測介質(zhi)溫度、黏(nián)度、密度(dù)、壓力等(deng)物理參(cān)數的影(yǐng)響且其(qi)示值在(zài)一定的(de)電導率(lǜ)範圍内(nèi)與被标(biāo)定的液(yè)體種類(lei)無關等(deng)特點,還(hai)具有體(ti)積小、耐(nai)高溫高(gāo)壓、流場(chang)不對稱(chēng)對測量(liang)精度影(yǐng)響較小(xiǎo)的優點(dian),可以作(zuò)爲獨立(li)設備進(jin)行井下(xià)測量,也(ye)可以作(zuò)爲複雜(za)智能測(ce)調系統(tǒng)的數據(jù)采集終(zhōng)端。其基(ji)本原理(lǐ)是基于(yu)法拉第(dì)電磁感(gǎn)應定律(lü),即當導(dao)電液體(tǐ)流過磁(ci)場作切(qiē)割磁力(li)線運動(dòng)時,則在(zai)垂直于(yu)流速向(xiang)量和磁(cí)場向量(liàng)的方向(xiang)上會産(chan)生一個(gè)與流量(liang)大小成(chéng)正比的(de)感應電(dian)動勢,其(qi)表達式(shi)爲
 
　　式中(zhōng):Ɛab爲感應(yīng)電動勢(shi);α爲電極(ji)1的位置(zhi)坐标;b爲(wèi)電極2的(de)位置坐(zuò)标;B爲流(liú)體微元(yuán)處的磁(cí)場強度(dù);V爲流體(ti)微元的(de)速度;dl爲(wèi)流體微(wēi)元的長(zhǎng)度。
　　因此(cǐ)可知，通(tong)過測得(de)感應電(dian)動勢的(de)大小，即(jí)可測得(dé)流量大(dà)小。
　　目前(qian),在電磁(ci)流量計(jì)方面的(de)有限元(yuan)建模研(yan)究較少(shǎo)。1996年，MICHALSKI等基(jī)于有限(xian)元建立(li)的不同(tong)形狀和(hé)尺寸的(de)流體管(guan)道數值(zhí)模型對(duì)勵磁線(xiàn)圈的橫(héng)截面形(xíng)狀進行(háng)尋優,以(yǐ)獲得均(jun1)勻的矢(shi)量積3];2002年(nián)，MICHALSKI等用有(you)限元方(fang)法建立(lì)了電磁(ci)流量計(ji)勵磁線(xiàn)圈的3D混(hun)合數學(xué)模型;2009年(nián),金甯德(dé)等用Ansys對(duì)四電極(jí)外流式(shì)電磁流(liú)量計建(jian)立了二(er)維有限(xian)元模型(xing),得出了(le)數值模(mo)拟結果(guo),提出了(le)四電極(ji)外流式(shì)電磁流(liu)量計的(de)理論分(fen)析方法(fǎ)(但這個(gè)模型無(wu)法進行(háng)仿真實(shi)驗);邬惠(huì)峰等建(jiàn)立了普(pu)通工業(ye)内流式(shi)電磁流(liu)量計的(de)二維仿(pang)真模型(xing)°0(内流式(shi)和外流(liu)式因其(qí)應用的(de)場合不(bú)同,整個(gè)流量計(ji)的結構(gou)也不同(tong));2010年，張志(zhì)剛利用(yong)Matlab對四電(dian)極外流(liu)式電磁(ci)流量計(ji)權重函(hán)數分布(bu)情況進(jin)行了理(li)論推導(dao)和仿真(zhēn)計算,爲(wèi)進一步(bù)開展四(sì)電極外(wài)流式電(dian)磁流量(liàng)計的研(yan)究和開(kai)發設計(ji)奠定了(le)理論基(jī)礎”。大量(liang)研究表(biǎo)明，有限(xiàn)元方法(fǎ)是一種(zhong)研究電(diàn)磁流量(liàng)計的有(yǒu)效手段(duàn)。由于四(sì)電極外(wai)流式電(dian)磁流量(liang)計系統(tong)本身受(shòu)結構參(cān)數和電(diàn)氣參數(shù)等衆多(duō)參數的(de)影響,影(ying)響規律(lǜ)複雜,改(gǎi)變某--個(gè)參數就(jiù)需要變(biàn)換硬件(jian),實驗效(xiào)率低而(ér)且成本(ben)高。因此(ci)采用有(you)限元方(fang)法建立(lì)能反映(yìng)其特性(xìng)的多物(wu)理場仿(páng)真模型(xíng),開展電(dian)磁流量(liàng)計勵磁(cí)規律和(hé)三維尺(chi)度下磁(cí)場分布(bu)規律及(jí)影響因(yin)素研究(jiū),可優化(huà)磁場設(shè)計參數(shu),指導傳(chuan)感器的(de)實驗與(yǔ)設計,顯(xiǎn)著降低(dī)成本,提(ti)高開發(fā)準确率(lǜ)及效率(lü)。
1流量計(jì)場路耦(ǒu)合有限(xian)元模型(xíng)的建立(lì)
1.1三維實(shi)體模型(xing)的建立(lì)與簡化(hua)
　　電磁流(liú)量計實(shí)體模型(xing)中不僅(jin)包括線(xiàn)圈、線圈(quan)架、電極(ji)、測量管(guan)、絕緣套(tao)、空氣域(yu)、流場域(yu)等主要(yao)部件,還(hai)包括平(píng)衡柱體(ti)、平衡柱(zhu)套、電纜(lǎn)插頭過(guo)線塞座(zuo)過線塞(sāi)套等輔(fu)助零件(jiàn)。由于輔(fǔ)件對磁(ci)場和電(diàn)極的感(gan)應電動(dòng)勢沒有(you)影響,同(tóng)時各個(gè)主要部(bù)件上都(dōu)加工有(you)裝配特(te)征,且這(zhè)些特征(zheng)都對磁(cí)場和信(xin)号也沒(méi)有影響(xiǎng),因此爲(wei)了提高(gao)計算效(xiao)率，可對(dui)傳感器(qi)模型進(jin)行簡化(hua)。簡化後(hòu)的模型(xing)包括:1)線(xiàn)圈，如圖(tú)1a);2)線圈架(jià)，如圖1b);3)電(diàn)極,如圖(tú)1c);4)空氣域(yù)，如圖1d);5)流(liu)體域,如(ru)圖1e);6)測量(liang)管域,如(ru)圖1f)。
 
　　在Solidworks中(zhōng)建立了(le)簡化的(de)傳感器(qì)實體模(mo)型,然後(hou)将其導(dǎo)入強大(da)的網格(gé)劃分軟(ruan)件HyperMesh中進(jin)行布爾(ěr)運算和(he)網格劃(huà)分，由于(yu)實體模(mo)型導入(ru)後會丢(diū)失體信(xin)息，因此(cǐ)模型導(dao)入後要(yào)重新利(lì)用各個(ge)實體的(de)面重新(xin)生成體(tǐ)。
1.2有限元(yuan)模型的(de)前處理(lǐ)及設置(zhi)
　　Ansys在工程(chéng)領域強(qiang)大的求(qiú)解能力(lì)衆所周(zhou)知(8],故采(cai)用Ansys軟件(jiàn)作爲電(diàn)磁場求(qiú)解軟件(jian)。由于整(zheng)個有限(xiàn)元模型(xíng)中的各(gè)個部件(jiàn)都是三(san)維實體(ti),模型尺(chǐ)寸頗大(dà),在進行(hang)網格劃(hua)分時會(hui)有大量(liàng)網格産(chan)生,增大(dà)計算量(liàng),而該模(mó)型中除(chú)了流體(tǐ)域和電(diàn)極是計(jì)算域外(wài),其他部(bù)分都不(bu)需要參(cān)與計算(suàn),因此将(jiang)線圈、線(xiàn)圈架、電(diàn)極的網(wǎng)格大小(xiǎo)設置爲(wèi)2mm,空氣域(yu)的網格(gé)大小設(shè)置爲3mm,流(liu)體域的(de)網格大(da)小設置(zhi)爲1mm。
　　有限(xiàn)元網格(gé)的質量(liang)直接影(ying)響計算(suan)精度,采(cǎi)用自動(dòng)網格劃(huà)分,單元(yuan)形狀爲(wèi)四面體(ti),粗網格(ge)和細網(wǎng)格之間(jian)過渡并(bìng)不光滑(huá),因此将(jiāng)流體域(yu)和空氣(qi)域之間(jiān)的測量(liàng)管域的(de)網格單(dān)元大小(xiao)設置爲(wei)2mm。網格劃(huà)分後,導(dǎo)入Anrsys中進(jin)行單元(yuán)類型、材(cái)料、實常(cháng)數、載荷(he)、邊界條(tiáo)件和場(chǎng)路耦合(he)單元設(shè)置。線圈(quan)用銅線(xian)實現,匝(zā)數共6500匝(za),其截面(miàn)積爲2.72×10-4:mm²,體(ti)積爲1.49×10-5mm3;線(xian)圈坐标(biāo)系單獨(du)定義爲(wei)局部柱(zhù)坐标系(xi)，軸向爲(wèi)正Y方向(xiang),其餘部(bù)件的坐(zuo)标系使(shǐ)用全局(ju)笛卡爾(er)坐标系(xi)，軸向爲(wèi)正Y方向(xiang)。各個部(bu)件的材(cai)料參數(shù)設置見(jiàn)表1。
 
　　爲了(le)實現勵(lì)磁方式(shì)的可編(bian)程，需要(yao)把線圈(quan)單元耦(ǒu)合到電(dian)路,因此(ci)建立2個(ge)Circu124分别實(shi)現獨立(lì)電壓源(yuan)單元和(he)耦合單(dān)元,V;節點(diǎn)的電位(wei)定義爲(wei)0,然後将(jiang)線圈單(dān)元的任(rèn)意-一個(ge)節點定(ding)義爲耦(ou)合單元(yuan)的K節點(diǎn)以實現(xian)耦合,具(jù)體如圖(tu)2所示。.
 
2模(mo)型校驗(yàn)
　　爲了保(bǎo)證模型(xíng)的正确(què)率,對建(jian)立的四(sì)電極外(wai)流式電(diàn)磁流量(liang)傳感器(qi)的三維(wei)有限元(yuán)模型,從(cong)2個方面(miàn)進行了(le)校驗:首(shǒu)先,給有(you)限元模(mó)型施加(jiā)恒值電(diàn)流激勵(li),選用靜(jìng)态求解(jiě)類型,将(jiang)模型最(zui)外圈節(jiē)點的Ax,Ay,Ax自(zì)由度均(jun1)設爲0,選(xuan)擇所有(yǒu)單元後(hou)進行求(qiu)解，然後(hòu)在後處(chu)理器中(zhong)讀入結(jie)果，畫出(chu)電極附(fù)近的磁(cí)場;分布(bù),如圖3所(suǒ)示,磁場(chǎng)分布符(fú)合金甯(ning)德等數(shù)值分析(xi)的結果(guo)'5],如圖4所(suo)示;其次(ci),在現有(yǒu)模型基(jī)礎上加(jiā)密網格(ge)單元，所(suo)得感應(ying)電動勢(shì)大小前(qián)後誤差(cha)小于5%,從(cong)而保證(zhèng)有限元(yuan)計算結(jié)果不受(shòu)網格質(zhì)量變化(hua)的影響(xiǎng)。綜上所(suo)述,該有(you)限元模(mo)型是準(zhǔn)确的,可(kě)用來進(jin)行仿真(zhen)研究。
 
3不(bú)同流速(su)下模型(xíng)的響應(yīng)及誤差(chà)分析
　　耦(ǒu)合分析(xī)分2種方(fang)法:強耦(ǒu)合(或稱(cheng)緊耦合(he))和弱耦(ou)合(或稱(cheng)松耦合(hé))。強耦合(hé)通過單(dān)元矩陣(zhèn)或荷載(zǎi)向量把(ba)耦合作(zuò)用構造(zào)到控制(zhi)方程中(zhong),然後對(duì)控制方(fang)程直接(jiē)求解,其(qí)缺點是(shì)在構造(zao)控制方(fang)程過程(chéng)中常常(chang)不得不(bú)對問題(ti)進行某(mǒu)些簡化(hua),有時候(hòu)計算準(zhun)确程度(du)較難保(bao)證。弱耦(ǒu)合是在(zai)每一步(bu)内分别(bié)對每一(yi)種場方(fang)程進行(hang)一次求(qiú)解,通過(guò)把第1個(ge)物理場(chang)的結果(guǒ)作爲外(wài)荷載加(jiā)于第2個(ge)物理場(chang)來實現(xian)2個場的(de)耦合。其(qí)優點是(shi)可以利(li)用現有(yǒu)的通用(yòng)流場和(he)電磁場(chǎng)軟件,并(bing)且可以(yi)分别對(duì)每--個軟(ruǎn)件單獨(dú)地制定(ding)合适的(de)求解方(fāng)法;缺點(diǎn)是計算(suàn)過程比(bi)較複雜(za)。強耦合(he)通常适(shi)合于對(duì)耦合場(chǎng)的理論(lùn)分析,弱(ruo)耦合适(shi)用于對(duì)耦合場(chang)的數值(zhí)計算。
　　仿(pang)真對象(xiàng)的外徑(jìng)尺寸是(shi)38mm,其工作(zuo)的管道(dào)内徑爲(wei)46mm,根據截(jié)面積相(xiang)等的原(yuán)則,其等(děng)效管徑(jing)爲26mm。當雷(lei)諾數Re<2000時(shi),管道内(nei)流動狀(zhuàng)态爲層(ceng)流;當4000>Re>2000時(shi),管道内(nèi)流動狀(zhuàng)态不确(que)定;當Re>4000時(shí),管道内(nei)流動狀(zhuang)态爲湍(tuān)流。當流(liú)動狀态(tài)爲湍流(liu)時,由可(ke)計算出(chu)紊流流(liú)動對應(ying)的最小(xiao)平均流(liu)速V=0.092m/s.
　　式中(zhōng):V爲平均(jun)流速;D爲(wei)圓管直(zhí)徑,取26mm;ʋ爲(wei)運動黏(nian)度,取0.6×10-6m2/s。
　　因(yīn)此,當管(guan)道内平(ping)均流速(su)V>0.092m/s時,管道(dao)内的流(liu)動狀态(tài)爲紊流(liu);事實上(shang)，四電極(ji)外流式(shì)電磁流(liú)量計在(zai)工作的(de)時候，管(guan)道内大(da)多數的(de)流動速(sù)度都大(dà)于這個(ge)值。當管(guǎn)道内的(de)流動狀(zhuang)态爲紊(wen)流時,用(yòng)CFD軟件進(jin)行流場(chang)分析、計(jì)算,通過(guò)CFD模拟，可(ke)以分析(xi)并且顯(xian)示流體(tǐ)流動過(guò)程中發(fā)生的現(xian)象,及時(shi)預測流(liú)體在模(mo)拟區域(yù)的流動(dòng)性能[10],用(yong)有限元(yuán)軟件Ansys中(zhōng)的FLOTRANCFD模塊(kuai)對其流(liu)場進行(háng)仿真分(fen)析,計算(suàn)結果如(ru)圖5所示(shi)。
 
　　在近壁(bì)0.2mm處速度(du)較小，其(qi)餘位置(zhi)都接近(jìn)平均速(sù)度。基于(yú)此,可以(yi)将流體(tǐ)等效爲(wei)一個勻(yún)速導體(ti),用強耦(ǒu)合的方(fang)法進行(hang)電磁流(liu)場耦合(he)的分析(xī)。
　　在紊流(liu)場共選(xuan)定了10個(ge)不同的(de)流量值(zhí),獨立電(diàn)壓源編(biān)程爲兩(liǎng)值矩形(xíng)波,幅值(zhí)爲15V,頻率(lǜ)爲1Hz,對流(liú)量數據(ju)進行了(le)仿真計(ji)算，并在(zài)實驗台(tai)上得出(chū)了實驗(yàn)數據,實(shí)驗台采(cǎi)用精度(du)爲0.5%的電(diàn)磁流量(liàng)計讀取(qǔ)流量值(zhi),用信号(hao)處理電(diàn)路采集(jí)四電極(jí)外流式(shì)電磁流(liú)量計的(de)感應電(dian)動勢信(xìn)号,通過(guo)串口輸(shū)入到計(ji)算機顯(xian)示,實驗(yàn)台原理(li)圖如圖(tú)6所示，最(zuì)後對這(zhe)2種數據(ju)進行了(le)誤差分(fen)析,結果(guo)見表2。
 
　　從(cong)實驗結(jié)果和仿(pang)真結果(guǒ)的誤差(cha)來看,流(liú)速較小(xiao)的時候(hou)誤差非(fēi)常大,随(suí)着流速(sù)的加快(kuai)，誤差逐(zhu)漸減小(xiao)。這是因(yin)爲流速(sù)越大,流(liu)場就越(yuè)趨近于(yú)勻速場(chang),仿真計(jì)算的方(fang)法越接(jiē)近真實(shí)情況。誤(wù)差一方(fang)面是由(yóu)仿真模(mo)型的簡(jian)化引起(qi)的,另一(yi)方面是(shi)由信号(hào)處理電(diàn)路引起(qǐ)的,仿真(zhēn)模型反(fǎn)映了實(shí)際的情(qíng)況,可以(yǐ)用于勵(li)磁技術(shu)實驗等(deng)的理論(lun)分析。
4結(jie)語
　　通過(guò)有限元(yuán)方法建(jiàn)立了四(si)電極電(dian)磁流量(liang)計的仿(pang)真模型(xing),從2個方(fang)面對模(mo)型進行(háng)了校驗(yan),驗證了(le)模型的(de)正确率(lǜ)。在不同(tóng)平均流(liú)速下,用(yong)強耦合(he)的方法(fa)仿真計(jì)算了模(mó)型的響(xiǎng)應,并計(jì)算了誤(wu)差。研究(jiu)表明,該(gāi)有限元(yuán)模型在(zài)一定程(chéng)度上可(ke)用于電(dian)磁結構(gòu)流場的(de)仿真分(fèn)析,流速(sù)較小時(shí),有限元(yuán)仿真結(jié)果與實(shi)驗誤差(chà)較大;流(liú)速較大(dà)時,流場(chǎng)趨近于(yu)勻速場(chang),仿真結(jie)果與實(shí)驗結果(guo)誤差較(jiao)小。

本文(wen)來源于(yu)網絡,如(ru)有侵權(quán)聯系即(ji)删除！