感應電流eit中激勵線圈的設計

1、<p>　　感應電流EIT中激勵線圈的設計</p><p>　　作者：向海燕　董秀珍　秦明新　尤富生　史學濤　付峰　劉銳崗　馬建英 </p><p>　　【關鍵詞】 電阻抗 </p><p>　　關鍵詞: 感應電流;電阻抗;體層攝影術;激勵線圈 </p><p>　　摘 要:目的 研究在感應電流電阻抗斷層成像（induced

2、 current electrical impedance tomography，ICEIT）中，線圈的形狀、規(guī)格和位置對整體性能的影響. 方法 研究了用于I-CEIT系統(tǒng)的兩種線圈.以敏感矩陣的條件數為原則，對圓形線圈的配置進行了研究.計算了另一種有軸線圈內感應電流分布及邊界電壓，探討了限制分辨能力的因素.結果 找到了選用圓形線圈時對線圈大小和位置的要求，證實了有軸線圈能增大目標中心區(qū)域的電流. 結論 用恰當的配置，使用圓形線圈

3、可使敏感矩陣的條件數較小，而用有軸線圈可以提高中心區(qū)域的敏感度和徑向分辨能力. </p><p>　　Keywords:induced current;electrical impedance;tomogra-phy;exciting coil </p><p>　　Abstract:AIM To research the effect of the shape，size and conf

4、iguration of the exciting coils on the perform ance of the system in induced current electrical impedance tomography. </p><p>　　METHODS Two kinds of exciting coils were studied.Circu-lar coil’s configuration

5、 was studied on the basis of the condi-tion number of the sensitive matrix.For the axis coil，the currents in the object and the potential of the boundary were calculated.Furthermore，the factors confining the distin-guish

6、ability were discussed.RESULTS The appropriate size and configuration of cirlular coils was found and axis coil was able to enhance the currents in the center of the object.CON┐CLUSION Circular coil of</p><p&g

7、t;<b>　　0 引言 </b></p><p>　　感應電流電阻抗斷層成像（induced current elec-trical impedance tomography，ICEIT）是一種新的EIT技術，它在被測目標的外圍放置若干個激勵線圈，對其施加時變電流，在空間產生交變磁場，從而在被測目標內激勵出感應電流.測量目標表面相鄰電極的電壓差，并用此數據進行目標區(qū)域電導率圖像重建［1-5

8、］ .與傳統(tǒng)的注入式EIT（applied current elec-trical impedance tomography，ACEIT）相比，ICEIT具有以下優(yōu)勢 </p><p>　?。?，3，5］ :①成像目標內的電流不受電極處的電流密度的限制，因而有可能使用更大的電流密度以提高信噪比;②由于周圍電極僅測量輸出電壓，不用于電流驅動，所以可以優(yōu)化電極設計;③通過改變線圈的形狀和位置，使空間磁場發(fā)生改變，從而

9、改變目標內的電流分布，提取某一部分的細節(jié);④當成像目標外有屏蔽層時，選擇適當頻率的驅動電流，就可使屏蔽層對感應電流密度分布影響不大，從而可能得到比較理想的成像結果.因此，ICEIT具有一系列注入式EIT所不具有的優(yōu)點，具有良好的應用前景.激勵磁場分布由線圈配置決定，它對電導率變化的敏感度和圖像重建中的病態(tài)問題都有重要的影響［2，3，6］ .本文中，我們首先討論了圓形線圈的大小和位置對敏感矩陣條件數的影響，接著計算了使用有軸線圈驅動時電導

10、率擾動區(qū)域為同心圓時的邊界電壓和電流分布，并對線圈的分辨能力進行了探討. </p><p>　　1 兩種可用于ICEIT的線圈 </p><p>　　1.1 圓形線圈 圓形線圈是ICEIT中最早也是最常用的線圈，如Fig1.3個半徑相等的線圈等角度放置，線圈的中心位于和目標同中心的半徑為RS 的圓周上，線圈與線圈之間的夾角為120°.目標表面聯(lián)接E個電極，對每一種激勵磁場，有（

11、E-1）個獨立測量數據.若線圈數為P，那么就有P種磁場分布，所以總的獨立測量數就為（E-1）［3］ .圓形線圈的大小和位置對敏感矩陣的條件數有很大的影響［6］ .使用3個半徑均為0.36cm（R </p><p>　　C =0.36cm）的線圈，對半徑為0.12cm的目標進行仿真實驗.所用剖分模型為541節(jié)點，1016單元.測量電極數為16，則獨立測量數為45.當RS 取不同的值時，計算敏感矩陣，并且對敏感矩陣進

12、行奇異值分解（SVD）.當RS =0時，因為3個線圈重合，所以只得到15個非零奇異值.當RS 逐漸增加到0.04，0.08，0.12，0.16和0.20cm時，它們對應的條件數迅速減少到4483，1845，899，438和252. </p><p>　　當線圈和目標的最小距離（穿過線圈圓心和目標圓心的直線被截在線圈和目標外緣之間的部分）固定時，考察線圈半徑變化時條件數的變化.固定線圈和目標間的最小距離為0.01c

13、m.當RC 分別取0.15，0.18，0.26，0.32和0.36cm時，條件數分別為308，218，192，187和185.因為大線圈的遠端對目標內的感應電流的影響越來越弱，所以當RC 增大到一定程度（此處為0.36cm），使用更大的線圈并不能明顯改善奇異值變化曲線.于是，可以得到以下結論，如果使用等角度放置的圓形線圈，應該使線圈半徑盡量大，并且使目標和線圈間的距離盡可能小，這樣才能使敏感矩陣的條件數盡可能小. </p>

14、<p>　　1.2 有軸的圓形線圈［6］ 用圓形線圈進行激勵時，邊緣處的電流密度較大，而中央區(qū)域相對較弱，如Fig3.為了增強中心區(qū)域電流，可以采用如Fig2所示的線圈，我們稱為有軸線圈.在線圈軸附近，磁場的方向相反，所以感應電流的方向也相反.中心的磁場變化率加大，相應就增強了中央區(qū)域的感應電流.而且，在線圈軸附近，感應電流的徑向分量得到了明顯增強.通過旋轉線圈，或者在目標周圍不同角度放置相類似的線圈，就可得到目標內不同

15、的感應電流分布，從而從目標表面采集到不同的數據用于成像. </p><p>　　2 對有軸線圈的進一步分析 </p><p>　　2.1 計算邊界電壓 同心圓電導率擾動問題的解在 兩個假設前提下可由解析方法得到.這兩個假設，一是假設總磁場主要由線圈電流產生，感應電流對總磁場的貢獻可以忽略;第二，在EIT的驅動頻率范圍內，位移電流就可以忽略［2，3，7，8］ </p>&l

16、t;p>　　在擾動區(qū)域為同心圓時，可以用傅立葉展開法（series expansion method，SEM）來求目標內的標量電位V.標量電位V在兩個區(qū)域可分別表示為:V=V1 0r R1 和V=V2 R1 r R2 兩個區(qū)域的解可各自獨立表述為［2，7］ :V1 （rθ）=Σ∞m=1 rm ［em cos（mθ）+fm sin（mθ）］V2 （r，θ）=Σ∞m=1 （rm am +r-m bm ）cos（mθ）+（rm cm

17、 +r-m dm ）sin（mθ）目標邊界條件和介質分界面處的條件為:V2 r n （R2 ）（σ2 V </p><p>　　2 r-V1 r） r=R 1 -（σ2 -σ 1 ）ωAn （R1 ）V1 （R1 ）=V1 （R1 ）系數am ，b m ，cm ，dm ，em ，f </p><p>　　m 與R1 ，R2 ，R3 ，σ1 ，σ2 ，ω和I都有關，可通過計算得到［2，7］

18、 .當m取到5時，也就是對傅立葉序列的前五項求和，就能夠達到優(yōu)于99%的精度［7］ . </p><p>　　2.2 計算目標內的電流分布 因為E=- V-JωA，所以計算出 V和A，就得到了電場的分布.而V=［Vx Vy ］，其中:Vx =V x=V r&#12539;r x+Vθ&#12539;θx Vy =V y=V r&#12539;r y+Vθ&#12539;θy矢量磁

19、位A=μ0I4π ∫θ2θ1 →dl R，它的兩個分量可按下式計算［2］ :Ax =μ 0 4πΣN n=1 -sin（n2πN）R R3 2πN A y =μ0 4πΣN n=1 cos（n2πN）R R3 2πN最后由J =σE，得到電流密度的分布. </p><p>　　Fig3為偏心放置的圓形線圈在均勻目標內的感應電流分布，Fig4為有軸線圈在均勻目標內的感應電流分布.由Fig3可見，目標中心的電

20、流較弱，而邊緣處的電流較大，且電流沿圓周切向方向的分量較大.從Fig4可以看到在線圈軸附近，電流密度較大，中心區(qū)域的電流得到了增強，且電流的徑向分量得到增強. </p><p>　　圖1 － 圖4 略 </p><p>　　2.3 有軸線圈的分辨能力 設背景電導率σ2 =S&#12539;m-1 ，目標半徑R2 =1m.當σ1 =σ2 時，也就是均勻背景時，目標外圍

21、R2 =1m處的電位記為V1 .當擾動區(qū)域電導率σ1 =10σ2 ，半徑為R1 時，測得R2 處電位為V </p><p>　　2 .Isaacson定義分辨能力為V1 和V2 的差Ud 的范數‖Ud ‖［9］ .當電壓差Ud 大于測量精度時，不均勻區(qū)域就可以分辨出來.反之，信號將被噪聲淹沒，這就決定了在一定的測量精度條件下能檢測出的最小目標的大小. </p><p>　　圖5 － 圖8

22、 略 </p><p>　　Fig5顯示的是擾動半徑R1 從0.05m到1m之間變化時Ud 的取值.Fig6顯示的是擾動區(qū)域電導率不變而擾動區(qū)域半徑變化時‖Ud ‖的變化.從圖中可以看出，感應電流EIT具有注入EIT所不具有的一個特性.當電導率擾動區(qū)域半徑增大時，注入EIT能夠容易的檢測出來.但是，在感應電流EIT中，最大分辨能力并不出現在R1 =1m處，而是出現在R1 =0.6m和R1 =0.7m之間，當R1

23、=1m時減小到零. </p><p>　　2.4 線圈半徑對邊界電壓和電場的影響 設R1 =0.6m，R2 =1m，σ1 =10S&#12539;m-1 ，σ2 =1S&#12539;m-1 .當R3 取1.1，1.2，1.3，1.4，1.5，1.6，1.7，1.8，1.9，2.0，2.2，2.4，2.6，2.8，3.0，3.5和4.0m時，求‖Ud ‖和 Ud / V . </p>

24、;<p>　　Fig7表明當線圈半徑增大時，‖Ud ‖的值逐漸變小.Fig8顯示的是在不同的線圈半徑下 Ud / V 的變化.當半徑在1.1～2.0之間變化時， Ud / V 集中在1%～4%之間變化. </p><p><b>　　3 討論 </b></p><p>　　我們探討了在感應電流EIT中兩種不同的線圈配置.使用圓形線圈時，3個等半徑的圓

25、形線圈等角度放置.為了得到較小的條件數，線圈要盡量靠近成像目標，線圈半徑要足夠大.但是，當半徑超過一定值時（視目標大小而定），條件性能的改變就不明顯了.采用有軸線圈可以給中央區(qū)域引入較大的電流，從而提高測量對中心區(qū)域擾動的敏感性.在線圈軸附近的大的徑向電流，也可以提高徑向分辨率. </p><p>　　更進一步的，我們可以采用多個弧形線圈來驅動，每個線圈內的電流可以各不相同.這樣，就可以像注入式EIT一樣來優(yōu)化驅

26、動模式.未來的工作應致力于研究最優(yōu)的線圈數目和最優(yōu)的驅動電流模式. </p><p><b>　　參考文獻: </b></p><p>　　［1］Purvis WR，Tozer RC，Anderson DK，Freeston IL，Induced current impedance imaging ［J］.IEE Proc A Sci Meas Technol，1

27、993;140（2）:135-141. </p><p>　?。?］Gencer NG，Kuzuoglu M，Ider YZ.Electrical impedance tomog-raphy using induced currents ［J］.IEEE Trans Med Imaging，1994;13（2）:338-350. </p><p>　?。?］Gencer NG，Ider Y

28、Z，Willamson SJ.Electrical impedance to-mography:induced-current imaging achieved with a multiple coil system ［J］.IEEE Trans Biomed Eng，1996;43（2）:139-149. </p><p>　　［4］Freeston IL，Tozer RC.Impedance imaging

29、using induced cur-rents ［J］.Physiol Meas，1995;16（Suppl3A）:257-266. </p><p>　?。?］Wenxin R，Robert G，Andy A.Experimental evaluation of two iterative reconstruction methods for induced current electrical impedan

30、ce tomography ［J］.IEEE Trans Med Imaging，1996;15（2）:180-187. </p><p>　　［6］Gencer NG，Ider YZ.A comparative study of several exciting magnetic fields for induced current EIT ［J］.Physiol Meas，1995;15（suppl2A）:5

31、1-57. </p><p>　?。?］Eyuboglu BM，Koksal A，Demirbilek M.Distinguishability anl-ysis of an induced EIT system using discrete coils ［J］.Phys Med Biol，2000;45（7）:1997-2009. </p><p>　?。?］Qin MX，Dong XZ

32、，Han T，Liu RG，Tang MX，Fu F，You FS，Shi XT.Simulating study of forward problem of electrical impedance tomography in brain ［J］.Di-si Junyi Daxue Xuebao（J Fourth Mil Med University），2001;22（5）:468-470. </p><p>