摘摘要:四電極外流式電磁流量計是一種新型的測量注入剖面流量的測井理想儀器,廣泛應(yīng)用于油田注水井、注聚井的流量測量。目前四電極外流式電磁流量計的研究主要在實(shí)際環(huán)境中開展,實(shí)驗(yàn)效率低、成本高。建立了四電極外流式電磁流量計的準(zhǔn)確的有限元模型,將強(qiáng)耦合的方法應(yīng)用在電磁結(jié)構(gòu)流體耦合.上,并在不同流速下開展了模型的響應(yīng)及誤差分析。研究表明,該有限元模型在一定程度上可用于電磁流量計的流場仿真分析。流速較小時,有限元仿真結(jié)果與實(shí)驗(yàn)誤差較大;流速較大時,流場趨近于勻速場,仿真結(jié)果與實(shí)驗(yàn)結(jié)果誤差較小。
　　在油田三次采油中,注聚合物驅(qū)油是提高原油采收率的重要手段之一,它比水驅(qū)效果提高了20%左右;現(xiàn)場實(shí)驗(yàn)表明,過去常用的注入剖面測井儀器已經(jīng)不適合注聚合物測井的剖面測試的要求。電磁流量計是一種新型的測量注入剖面的儀器,較好地解決了聚合物注入剖面的測井問題。
　　四電極外流式電磁流量計是針對油田應(yīng)用開發(fā)的一種特殊電磁流量計,其不僅具有普通工業(yè)電磁流量計無節(jié)流阻流,不易堵塞,耐腐蝕性好,測量精度不受被測介質(zhì)溫度、黏度、密度、壓力等物理參數(shù)的影響且其示值在一定的電導(dǎo)率范圍內(nèi)與被標(biāo)定的液體種類無關(guān)等特點(diǎn),還具有體積小、耐高溫高壓、流場不對稱對測量精度影響較小的優(yōu)點(diǎn),可以作為獨(dú)立設(shè)備進(jìn)行井下測量,也可以作為復(fù)雜智能測調(diào)系統(tǒng)的數(shù)據(jù)采集終端。其基本原理是基于法拉第電磁感應(yīng)定律,即當(dāng)導(dǎo)電液體流過磁場作切割磁力線運(yùn)動時,則在垂直于流速向量和磁場向量的方向上會產(chǎn)生一個與流量大小成正比的感應(yīng)電動勢,其表達(dá)式為
 
　　式中:Ɛab為感應(yīng)電動勢;α為電極1的位置坐標(biāo);b為電極2的位置坐標(biāo);B為流體微元處的磁場強(qiáng)度;V為流體微元的速度;dl為流體微元的長度。
　　因此可知，通過測得感應(yīng)電動勢的大小，即可測得流量大小。
　　目前,在電磁流量計方面的有限元建模研究較少。1996年，MICHALSKI等基于有限元建立的不同形狀和尺寸的流體管道數(shù)值模型對勵磁線圈的橫截面形狀進(jìn)行尋優(yōu),以獲得均勻的矢量積3];2002年，MICHALSKI等用有限元方法建立了電磁流量計勵磁線圈的3D混合數(shù)學(xué)模型;2009年,金寧德等用Ansys對四電極外流式電磁流量計建立了二維有限元模型,得出了數(shù)值模擬結(jié)果,提出了四電極外流式電磁流量計的理論分析方法(但這個模型無法進(jìn)行仿真實(shí)驗(yàn));鄔惠峰等建立了普通工業(yè)內(nèi)流式電磁流量計的二維仿真模型°0(內(nèi)流式和外流式因其應(yīng)用的場合不同,整個流量計的結(jié)構(gòu)也不同);2010年，張志剛利用Matlab對四電極外流式電磁流量計權(quán)重函數(shù)分布情況進(jìn)行了理論推導(dǎo)和仿真計算,為進(jìn)一步開展四電極外流式電磁流量計的研究和開發(fā)設(shè)計奠定了理論基礎(chǔ)”。大量研究表明，有限元方法是一種研究電磁流量計的有效手段。由于四電極外流式電磁流量計系統(tǒng)本身受結(jié)構(gòu)參數(shù)和電氣參數(shù)等眾多參數(shù)的影響,影響規(guī)律復(fù)雜,改變某--個參數(shù)就需要變換硬件,實(shí)驗(yàn)效率低而且成本高。因此采用有限元方法建立能反映其特性的多物理場仿真模型,開展電磁流量計勵磁規(guī)律和三維尺度下磁場分布規(guī)律及影響因素研究,可優(yōu)化磁場設(shè)計參數(shù),指導(dǎo)傳感器的實(shí)驗(yàn)與設(shè)計,顯著降低成本,提高開發(fā)準(zhǔn)確率及效率。
1流量計場路耦合有限元模型的建立
1.1三維實(shí)體模型的建立與簡化
　　電磁流量計實(shí)體模型中不僅包括線圈、線圈架、電極、測量管、絕緣套、空氣域、流場域等主要部件,還包括平衡柱體、平衡柱套、電纜插頭過線塞座過線塞套等輔助零件。由于輔件對磁場和電極的感應(yīng)電動勢沒有影響,同時各個主要部件上都加工有裝配特征,且這些特征都對磁場和信號也沒有影響,因此為了提高計算效率，可對傳感器模型進(jìn)行簡化。簡化后的模型包括:1)線圈，如圖1a);2)線圈架，如圖1b);3)電極,如圖1c);4)空氣域，如圖1d);5)流體域,如圖1e);6)測量管域,如圖1f)。
 
　　在Solidworks中建立了簡化的傳感器實(shí)體模型,然后將其導(dǎo)入強(qiáng)大的網(wǎng)格劃分軟件HyperMesh中進(jìn)行布爾運(yùn)算和網(wǎng)格劃分，由于實(shí)體模型導(dǎo)入后會丟失體信息，因此模型導(dǎo)入后要重新利用各個實(shí)體的面重新生成體。
1.2有限元模型的前處理及設(shè)置
　　Ansys在工程領(lǐng)域強(qiáng)大的求解能力眾所周知(8],故采用Ansys軟件作為電磁場求解軟件。由于整個有限元模型中的各個部件都是三維實(shí)體,模型尺寸頗大,在進(jìn)行網(wǎng)格劃分時會有大量網(wǎng)格產(chǎn)生,增大計算量,而該模型中除了流體域和電極是計算域外,其他部分都不需要參與計算,因此將線圈、線圈架、電極的網(wǎng)格大小設(shè)置為2mm,空氣域的網(wǎng)格大小設(shè)置為3mm,流體域的網(wǎng)格大小設(shè)置為1mm。
　　有限元網(wǎng)格的質(zhì)量直接影響計算精度,采用自動網(wǎng)格劃分,單元形狀為四面體,粗網(wǎng)格和細(xì)網(wǎng)格之間過渡并不光滑,因此將流體域和空氣域之間的測量管域的網(wǎng)格單元大小設(shè)置為2mm。網(wǎng)格劃分后,導(dǎo)入Anrsys中進(jìn)行單元類型、材料、實(shí)常數(shù)、載荷、邊界條件和場路耦合單元設(shè)置。線圈用銅線實(shí)現(xiàn),匝數(shù)共6500匝,其截面積為2.72×10-4:mm²,體積為1.49×10-5mm3;線圈坐標(biāo)系單獨(dú)定義為局部柱坐標(biāo)系，軸向?yàn)檎齓方向,其余部件的坐標(biāo)系使用全局笛卡爾坐標(biāo)系，軸向?yàn)檎齓方向。各個部件的材料參數(shù)設(shè)置見表1。
 
　　為了實(shí)現(xiàn)勵磁方式的可編程，需要把線圈單元耦合到電路,因此建立2個Circu124分別實(shí)現(xiàn)獨(dú)立電壓源單元和耦合單元,V;節(jié)點(diǎn)的電位定義為0,然后將線圈單元的任意-一個節(jié)點(diǎn)定義為耦合單元的K節(jié)點(diǎn)以實(shí)現(xiàn)耦合,具體如圖2所示。.
 
2模型校驗(yàn)
　　為了保證模型的正確率,對建立的四電極外流式電磁流量傳感器的三維有限元模型,從2個方面進(jìn)行了校驗(yàn):首先,給有限元模型施加恒值電流激勵,選用靜態(tài)求解類型,將模型最外圈節(jié)點(diǎn)的Ax,Ay,Ax自由度均設(shè)為0,選擇所有單元后進(jìn)行求解，然后在后處理器中讀入結(jié)果，畫出電極附近的磁場;分布,如圖3所示,磁場分布符合金寧德等數(shù)值分析的結(jié)果'5],如圖4所示;其次,在現(xiàn)有模型基礎(chǔ)上加密網(wǎng)格單元，所得感應(yīng)電動勢大小前后誤差小于5%,從而保證有限元計算結(jié)果不受網(wǎng)格質(zhì)量變化的影響。綜上所述,該有限元模型是準(zhǔn)確的,可用來進(jìn)行仿真研究。
 
3不同流速下模型的響應(yīng)及誤差分析
　　耦合分析分2種方法:強(qiáng)耦合(或稱緊耦合)和弱耦合(或稱松耦合)。強(qiáng)耦合通過單元矩陣或荷載向量把耦合作用構(gòu)造到控制方程中,然后對控制方程直接求解,其缺點(diǎn)是在構(gòu)造控制方程過程中常常不得不對問題進(jìn)行某些簡化,有時候計算準(zhǔn)確程度較難保證。弱耦合是在每一步內(nèi)分別對每一種場方程進(jìn)行一次求解,通過把第1個物理場的結(jié)果作為外荷載加于第2個物理場來實(shí)現(xiàn)2個場的耦合。其優(yōu)點(diǎn)是可以利用現(xiàn)有的通用流場和電磁場軟件,并且可以分別對每--個軟件單獨(dú)地制定合適的求解方法;缺點(diǎn)是計算過程比較復(fù)雜。強(qiáng)耦合通常適合于對耦合場的理論分析,弱耦合適用于對耦合場的數(shù)值計算。
　　仿真對象的外徑尺寸是38mm,其工作的管道內(nèi)徑為46mm,根據(jù)截面積相等的原則,其等效管徑為26mm。當(dāng)雷諾數(shù)Re<2000時,管道內(nèi)流動狀態(tài)為層流;當(dāng)4000>Re>2000時,管道內(nèi)流動狀態(tài)不確定;當(dāng)Re>4000時,管道內(nèi)流動狀態(tài)為湍流。當(dāng)流動狀態(tài)為湍流時,由可計算出紊流流動對應(yīng)的最小平均流速V=0.092m/s.
　　式中:V為平均流速;D為圓管直徑,取26mm;ʋ為運(yùn)動黏度,取0.6×10-6m2/s。
　　因此,當(dāng)管道內(nèi)平均流速V>0.092m/s時,管道內(nèi)的流動狀態(tài)為紊流;事實(shí)上，四電極外流式電磁流量計在工作的時候，管道內(nèi)大多數(shù)的流動速度都大于這個值。當(dāng)管道內(nèi)的流動狀態(tài)為紊流時,用CFD軟件進(jìn)行流場分析、計算,通過CFD模擬，可以分析并且顯示流體流動過程中發(fā)生的現(xiàn)象,及時預(yù)測流體在模擬區(qū)域的流動性能[10],用有限元軟件Ansys中的FLOTRANCFD模塊對其流場進(jìn)行仿真分析,計算結(jié)果如圖5所示。
 
　　在近壁0.2mm處速度較小，其余位置都接近平均速度。基于此,可以將流體等效為一個勻速導(dǎo)體,用強(qiáng)耦合的方法進(jìn)行電磁流場耦合的分析。
　　在紊流場共選定了10個不同的流量值,獨(dú)立電壓源編程為兩值矩形波,幅值為15V,頻率為1Hz,對流量數(shù)據(jù)進(jìn)行了仿真計算，并在實(shí)驗(yàn)臺上得出了實(shí)驗(yàn)數(shù)據(jù),實(shí)驗(yàn)臺采用精度為0.5%的電磁流量計讀取流量值,用信號處理電路采集四電極外流式電磁流量計的感應(yīng)電動勢信號,通過串口輸入到計算機(jī)顯示,實(shí)驗(yàn)臺原理圖如圖6所示，最后對這2種數(shù)據(jù)進(jìn)行了誤差分析,結(jié)果見表2。
 
　　從實(shí)驗(yàn)結(jié)果和仿真結(jié)果的誤差來看,流速較小的時候誤差非常大,隨著流速的加快，誤差逐漸減小。這是因?yàn)榱魉僭酱?流場就越趨近于勻速場,仿真計算的方法越接近真實(shí)情況。誤差一方面是由仿真模型的簡化引起的,另一方面是由信號處理電路引起的,仿真模型反映了實(shí)際的情況,可以用于勵磁技術(shù)實(shí)驗(yàn)等的理論分析。
4結(jié)語
　　通過有限元方法建立了四電極電磁流量計的仿真模型,從2個方面對模型進(jìn)行了校驗(yàn),驗(yàn)證了模型的正確率。在不同平均流速下,用強(qiáng)耦合的方法仿真計算了模型的響應(yīng),并計算了誤差。研究表明,該有限元模型在一定程度上可用于電磁結(jié)構(gòu)流場的仿真分析,流速較小時,有限元仿真結(jié)果與實(shí)驗(yàn)誤差較大;流速較大時,流場趨近于勻速場,仿真結(jié)果與實(shí)驗(yàn)結(jié)果誤差較小。

本文來源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！