脑磁感应成像数学模型及仿真结果.docVIP

脑磁感应成像数学模型及仿真结果 　 作者： 殷朝庆，董秀珍，刘锐岗，尤富生，付峰，史学涛，王聪，刘程睿 【关键词】 磁感应；电阻抗；脑磁图描记术；仿真；磁势 　　 【Abstract】 AIM: To study the influence of skull on the border in the brain magnetic induction tomography(MIT). METHODS: The head was simplified and approximated to a layered and symmetrical medium sphere，and the excited coil and the inducted coil were located in the border of the sphere. We obtained the field equations by the method of time harmonic electromagnetic field and the distribution of the excited magnetic field by the finite element analysis software. RESULTS: We obtained the changes of the magnetic potential on the border and the distribution of the excited coil and the trend of the magnetic potential on the border. The trend was the same with or without skull axistence. CONCLUTION: The changing trend of the magnetic potential on the border wouldnt be affected by the skull and the change can be detected. 【Keywords】 magnetic induction; electric impedance; magnetoencephalography; simulation; magnetic potential 【摘要】目的： 探讨脑磁感应断层成像中颅骨存在对边界磁场的影响. 方法： 将头颅近似为一个分层均匀介质球，激励和测量线圈位于球边界均匀放置. 采用时谐电磁场分析方法推导出该模型的有限元方程并利用有限元分析软件仿真激励磁场分布. 结果： 获得了该模型的有限元方程、激励磁场的分布和边界上磁势的变化. 有无颅骨时导体位置的变化均引起边界上磁势发生相应的变化，且两者变化趋势相同. 结论： 颅骨存在不影响边界上磁势的变化趋势. 边界上磁势的变化可以被检测到. 【关键词】 磁感应；电阻抗；脑磁图描记术；仿真；磁势 磁感应电阻抗断层成像是利用电磁感应原理，通过检测导体感应出的涡流的相位进行成像的技术［1］ . 该技术通过测量生物组织的感应磁场，根据重构算法利用计算机来表现被测组织电导率分布的一种新兴的成像方法［2-3］. 目前，国内外进行脑磁感应成像仿真研究的有：1997年Korzhenevskii和Cherepenin使用滤波反投影方法进行了MIT(磁感应电阻抗断层成像)的仿真研究；1999年Peyton使用加权反投影进行了MIT的仿真研究；1999年Gencer 和Tek使用NewtonRaphson方法进行了MIT的仿真研究；2000年A Korjenesky, V Cherepenin和S Sapetsky使用滤波反投影方法进行MIT的仿真研究；2003年第四军医大学阻抗成像课题组使用反投影方法实现了MIT的图像重构，获得了初步结果. 2004年Hermann Scharfetter, Robert Merwa 和Karl Pilz等人用GaussNewtononestep方法进行了MIT的仿真研究. 我们采用仿真软件对脑磁感应阻抗成像的激励线圈的激励磁场进行了仿真研究，为编程实现MIT的正向问题和MIT的图像重构提供理论参考. 1MIT数学模型和变分问题的离散化 1.1MIT的数学模型研究模型为一个圆面如图1a所示. 模型从外向里依次为空气域、颅骨层、脑实质，里面黑色部分表示导体. 激励线圈位于模型右侧. 根据电磁场理论：对于稳态时变场，从Maxwell方程组出发，以矢量位为求解对象，在忽略位移电流的情况下，Maxwell方程可以简化为： ×=1+2 ×=-Jω #12539;=0(1) 相应的场量关系可以表示为：=