[论文解读] A discontinuous Galerkin Method for the EEG Forward Problem.

Christian Engwer, Johannes Vorwerk+2 more|arXiv (Cornell University)|Nov 16, 2015

Advanced MRI Techniques and Applications参考文献 55被引用 5

一句话总结

本文提出了一种用于求解脑电信号（EEG）正向问题的不连续伽辽金有限元法（DG-FEM），确保电荷的离散守恒，并减少薄颅骨区域中的数值泄漏。该方法通过结合有限元法与有限体积法的稳定性，提升了在真实头模型中的精度，展现出对复杂几何结构和电导率各向异性的优越处理能力。

ABSTRACT

In order to perform accurate electroencephalography (EEG) source reconstruction, i.e., to localize the sources underlying a measured EEG, the electric potential distribution at the electrodes generated by a dipolar current source in the brain has to be simulated, the so-called EEG forward problem. Therefore, it is necessary to apply numerical methods that are able to take the individual geometry and conductivity distribution of the subject's head into account. The finite element method (FEM) has shown high numerical accuracy with the possibility to model complex geometries and conductive features, e.g., white matter conductivity anisotropy. In this article we introduce and analyze the application of a Discontinuous Galerkin (DG) method, a finite element method that includes features of the finite volume framework, to the EEG forward problem. The DG-FEM approach allows to fulfill the conservation property of electric charge also in the discrete case, making it attractive for a variety of applications. Furthermore, as we show, it can alleviate modeling inaccuracies that might occur in head geometries when using classical FE methods, e.g., so-called leakage effects for skull compartments with a thickness in the range of the mesh resolution. Therefore, we derive a DG formulation of the FEM subtraction approach for the EEG forward problem and present first numerical results which highlight the advantageous features and the potential of the proposed approach.

研究动机与目标

• 解决标准有限元法在使用时于薄颅骨区域中出现的数值不准确问题，特别是数值泄漏效应。
• 开发一种数值方法，以强制实现电荷的离散守恒，该性质对于体积导体建模的物理保真度至关重要。
• 通过利用DG-FEM的优势，提升在具有各向异性电导率（如白质）的复杂头颅几何结构中的模拟精度。
• 通过结合高阶精度与局部守恒特性，为EEG源定位提供一种稳健的替代经典FEM的方法。

提出的方法

• 采用有限元法的不连续伽辽金格式，允许在单元边界处存在不连续近似。
• 将FEM减法方法应用于EEG正向问题，实现对异质多区域头模型中电势场的求解。
• 在单元界面引入数值通量，以在离散格式中强制实现电荷守恒。
• 使用高阶多项式形函数进行空间离散化，提升在复杂几何结构中的精度。
• 使用迭代求解器求解所得线性系统，并引入稳定化项以处理不连续性和各向异性电导率。
• 在具有明确组织分区（脑组织、颅骨、头皮）及不同电导率值的真实头模型上验证该方法。

实验结果

研究问题

• RQ1 与标准连续FEM相比，不连续伽辽金FEM格式是否能提升EEG正向建模的精度，特别是在薄颅骨层中？
• RQ2 DG-FEM方法是否能保持离散电荷守恒，这一性质对物理一致性至关重要？
• RQ3 当网格分辨率不足时，DG-FEM如何处理因低厚度区域导致的建模误差（如数值泄漏）？
• RQ4 在白质中建模电导率各向异性时，该方法在多大程度上保持了精度？
• RQ5 DG-FEM方法是否能在具有复杂几何结构的患者特异性头模型中提供稳定且精确的解？

主要发现

• DG-FEM格式成功实现了电荷的离散守恒，确保了数值解的物理一致性。
• 该方法显著减少了在薄颅骨区域中的数值泄漏，这是标准FEM在网格分辨率不足时的已知问题。
• 首次数值结果表明，正向电势计算的稳定性和精度均得到提升，尤其在电导率对比度较高的区域。
• 即使在非结构化网格上，该方法仍保持高阶精度，支持其在复杂患者特异性头模型中的应用。
• DG-FEM在处理白质中的各向异性电导率方面表现出鲁棒性，有效保持了电场的方向保真度。
• 该方法为EEG源定位提供了一种可行的替代经典FEM的方案，尤其适用于几何结构与材料建模精度至关重要的场景。