《脑磁感应成像的数学模型及仿真结果》由会员分享，可在线阅读，更多相关《脑磁感应成像的数学模型及仿真结果（6页珍藏版）》请在金锄头文库上搜索。
1、脑磁感应成像的数学模型及仿真结果殷朝庆，董秀珍，刘锐岗，尤富生，付峰，史学涛，王聪，刘程睿【关键词】磁感应；电阻抗；脑磁图描记术；仿真；磁势【Abstrat】AI:Tstudytheinfluenefskullnthebrderinthebrainagnetiindutintgraphy(IT).ETHDS:Theheadassiplifiedandapprxiatedtalayeredandsyetrialediusphere，andtheexitedilandtheindutedilerelatedinthebrderfthesphere.ebtainedthefieldequatinsb
2、ytheethdftieharnieletragnetifieldandthedistributinftheexitedagnetifieldbythefiniteeleentanalysissftare.RESULTS:ebtainedthehangesftheagnetiptentialnthebrderandthedistributinftheexitedilandthetrendftheagnetiptentialnthebrder.Thetrendasthesaeithrithutskullaxistene.NLUTIN:Thehangingtrendftheagnetiptenti
3、alnthebrderuldntbeaffetedbytheskullandthehangeanbedeteted.【Keyrds】agnetiindutin;eletriipedane;agnetenephalgraphy;siulatin;agnetiptential【摘要】目的：讨论脑磁感应断层成像中颅骨存在对边界磁场的影响.方法：将头颅近似为一个分层均匀介质球，鼓励和测量线圈位于球边界均匀放置.采用时谐电磁场分析方法推导出该模型的有限元方程并利用有限元分析软件仿真鼓励磁场分布.结果：获得了该模型的有限元方程、鼓励磁场的分布和边界上磁势的变化.有无颅骨时导体位置的变化均引起边界上磁势发生
4、相应的变化，且两者变化趋势一样.结论：颅骨存在不影响边界上磁势的变化趋势.边界上磁势的变化可以被检测到.【关键词】磁感应；电阻抗；脑磁图描记术；仿真；磁势磁感应电阻抗断层成像是利用电磁感应原理，通过检测导体感应出的涡流的相位进展成像的技术1.该技术通过测量生物组织的感应磁场，根据重构算法利用计算机来表现被测组织电导率分布的一种新兴的成像方法2-3.目前，国内外进展脑磁感应成像仿真研究的有：1997年Krzhenevskii和herepenin使用滤波反投影方法进展了IT(磁感应电阻抗断层成像)的仿真研究；1999年Peytn使用加权反投影进展了IT的仿真研究；1999年Gener和Tek使用N
5、etnRaphsn方法进展了IT的仿真研究；2000年AKrjenesky,Vherepenin和SSapetsky使用滤波反投影方法进展IT的仿真研究；2022年第四军医大学阻抗成像课题组使用反投影方法实现了IT的图像重构，获得了初步结果.2022年HerannSharfetter,Rbertera和KarlPilz等人用GaussNetnnestep方法进展了IT的仿真研究.我们采用仿真软件对脑磁感应阻抗成像的鼓励线圈的鼓励磁场进展了仿真研究，为编程实现IT的正向问题和IT的图像重构提供理论参考.1IT数学模型和变分问题的离散化1.1IT的数学模型研究模型为一个圆面如图1a所示.模型从外向
6、里依次为空气域、颅骨层、脑本质，里面黑色局部表示导体.鼓励线圈位于模型右侧.根据电磁场理论：对于稳态时变场，从axell方程组出发，以矢量位为求解对象，在忽略位移电流的情况下，axell方程可以简化为：=1+2=-J=0(1)相应的场量关系可以表示为：=(2)2=(3)在上面的方程中，1为源电流密度，2为涡流电流密度，为磁导率，为导体的电导率.根据方程=0可以引入矢量位，并引入库仑标准=0.由以上方程组可以得到关于矢量位的微分方程：2=-01+j0边界条件为：=0|l=01.2条件变分问题离散化上述问题的条件变分问题为：F()=12(x)2+(y)2dxdy+j2s2dxdy-l2tdl-s1
7、dxdy=inl=0通过推导，得到有限元方程为：K+T=P式中,是由单元矩阵e,e,e合成的总体矩阵.e,e,e的一般表示是：Kers=Kesr=14e(brbs+rs)Ters=Tesr=j12(1+rs)r,s,l=i,j,Pers=Je13其中，为三角形单元的面积：br,bs,r,s与单元顶点的坐标有关；为角频率；e为单元磁导率；为单元电导率；rs为狄拉克函数，当rs时，rs=1，当r=s时，rs=0；J1为源电流密度.2FELAB仿真及结果分析2.1仿真模型FELAB的电磁场有限元分析是以axell方程组作为根据和出发点，其分析过程为：建立模型；设置模型每局部参数，施加边界条件；将模型
8、划分成有限个单元；求解.利用FELLAB分别建立半径为17.5，15.5，14，3的圆(图1a).每一层分别表示(从外向里)空气域、颅骨层、脑本质，黑色局部表示导体.模型边缘处的矩形为鼓励线圈.以往的研究说明，脑部介质的磁导率和真空磁导率0接近4.因此，把脑模型内每层的相对磁导率设为1，即：r=1.脑模型内每层电导率分别设为：空气域，0；颅骨层，0.0042s/；脑本质，0.75s/5-6导体的电导率设为2s/，鼓励频率为10Hz；电流密度设为200A/2.将模型划分为1079个单元(图1b，图2b)，求解时间为0.42s.图1和图2分别为当导体位于相应位置的磁感应强度的分布情况(采用均匀化)
9、.图3为不存在颅骨的模型、剖分和求解.2.2区域方程及边界条件仿真时，为脑模型内每一层设定相应的参数来确定每个区域的方程.根据前面的推导,空气域方程为：2=01脑模型内每一层和导体的方程为：2=j0为脑模型内每一层导体的电导率.在前面的条件变分问题中，边值问题的第二类边界条件和媒质间的交界面上的条件为自然边界条件，在泛函求极值过程中自动满足.因此，可以不考虑第二类边界条件和媒质的交界面条件.边值问题的第一类边界条件为强加边界条件.这样，第一类边界条件为模型的整个边界.在模型的周边上施加磁场边界条件：=0.2.3结果分析当模型内存在颅骨且导体中心位于坐标原点(0,0)，在y0区域，模型边界上z方
10、向的磁势随x的增大而增大；在y0的区域，模型边界上z方向的磁势随x的增大而减小；并且与y0的区域A：颅骨存在且导体中心位于(0,10)物理模型；B：模型的网格剖分；：模型内鼓励线圈产生的磁场分布.图1颅骨存在且导体中心位于(0,10)的仿真模型A：颅骨存在且导体中心位于(0,0)的物理模型；B：模型的网格剖分；：模型内鼓励线圈产生的磁场分布.图2颅骨存在且导体中心位于(0,0)的仿真模型A：没有颅骨且导体中心位于(0,0)的物理模型；B：模型的网格剖分；：模型内鼓励线圈产生的磁场分布.图3无颅骨存在且导体位于中心(0,0)的仿真模型.完全对称，如图4中间曲线所示.当导体中心位于(0,10)，模
11、型边界上z方向的磁势呈现不对称，偏向y0的区域，如图4上面曲线所示.当导体中心位于(0,-10)，模型边界上z方向的磁势偏向y0的区域，如图4下面曲线所示.图3为没有颅骨时的模型、网格剖分、求解.参数的设置与有颅骨存在设置的不同之处在于将颅骨层的电导率设为脑本质的电导率0.75s/.从图2,图3可以看出：在没有颅骨和存在颅骨的情况下，颅骨附近的磁力线图存在差异.进一步说明：颅骨的存在会影响到模型边界上的磁势.图5给出了没有颅骨时，导体中心分别位于(0,10),(0,0),(0,-10)，边界上z方向磁势的变化.计算说明：导体位于模型中心位置时，边界上磁势的改变最明显为76%图6.3讨论以上根据
12、建立的IT的数学模型，利用时谐电磁场分析方法基于以下假设得到该模型有限元方程：1，媒质为各向同性.2，忽略角频率对场的影响、忽略位移电流.第一个假设确保本构关系=和2=成立；第二个假设确保每层媒质的电导率、磁导率为常数，不受频率的影响；理论说明：位移电流和传导电流都可以产生磁场，但是位移电流相对于传导电流很小，在50Hz时，位移电流相对于传导电流只有110-7.以上的假设使问题简化，但并不失一般性.图6颅骨存在较无颅骨存在时边界磁势变化根据得到的有限元方程，利用有限元软件对脑模型内的鼓励磁场在模型边界上产生的磁势进展了仿真.结果给出了有颅骨存在和没有颅骨情况下，随着导体位置的不同，模型边界上z
13、方向磁势的分布和变化，如图4和图5所示.从边界上磁势的分布曲线可以看出，随着导体位置的不同，有颅骨的头模型和没有颅骨的头模型的边界上的磁势变化趋势一样.可以通过检测边界上的磁势反映出导体的位置.仿真过程中，模型的剖分、每层参数的设置的不同都会给计算结果带来一定的影响.特别是每层的电导率、相对磁导率的设置，这些设置在很大程度上影响到鼓励线圈在模型边界上产生的磁势的分布.计算的结果会影响到整个模型的每一点的磁势，而边界上磁势的变化趋势不会改变.【参考文献】1HerannS,Rbert,KarlP.AnetypefgradieterfrthereEivingiruitfagnetiindutintg
14、raphy(IT)J.Physileas,2022,26(2):S307-S318.2herepeninV,A3Deletrialipedanetgraphy(EIT)systefrbreastanerdetetinJ.Physialeas,2001,22(1):9-18.3KrjenevskyAV.agnetiindutintgraphyneiagingethdinbiediine.Pr，2ndrldngr：IndustrialPressTgraphy.2001,7(10):240-246.4TenfrdTS.InteratinfExtreelyLFrequenyEletriandagnetiFieldsithHuansJ.HealthPhys,1987,53(6):585-605.5冯远明，王明时.脑部磁刺激磁场和感应电场的初步研究J.中国生物医学工程学报，1995,14(1):74-75.6uffinBN.EentriSpheresdelsftheHndJ.IEEETransBiedEng,1991,38(9):871-878.
