[论文解读] Continuity, the Bloch-Torrey equation, and Diffusion MRI

Matt G. Hall|arXiv (Cornell University)|Aug 9, 2016

Advanced Neuroimaging Techniques and Applications参考文献 20被引用 4

一句话总结

本文提出了一种具有灵活通量项的广义Bloch-Torrey方程，通过统一多种扩散磁共振成像（dMRI）模型于单一连续介质框架之下。通过将不同的dMRI技术解释为通量的不同选择（例如，Fickian扩散、相干流动），该工作揭示了其背后的数学等价性，并为自旋动力学建模开辟了新途径，包括相位敏感的信号分析以及动量或能量等非传统守恒量的建模。

ABSTRACT

The Bloch equation describes the evolution of classical particles tagged with a magnetisation vector in a strong magnetic field and is fundamental to many NMR and MRI contrast methods. The equation can be generalised to include the effects of spin motion by including a spin flux, which typically contains a Fickian diffusive term and/or a coherent velocity term. This form is known as the Bloch-Torrey equation, and is fundamental to MR modalities which are sensitive to spin dynamics such as diffusion MRI. Such modalities have received a great deal of interest in the research literature over the last few years, resulting in a huge range of models and methods. In this work we make make use of a more general Bloch-Torrey equation with a generalised flux term. We show that many commonly employed approaches in Diffusion MRI may be viewed as different choices for the flux terms in this equation. This viewpoint, although obvious theoretically, is not usually emphasised in the diffusion MR literature and points to interesting new directions for methods development. In the interests of completeness, we describe all the necessary mathematical and physical framework -- including the method used to solve the various Bloch-Torrey equations. We then review the relevant portions of the diffusion MRI literature and show how different flux choices lead to different models. All the models considered here are continuum-based but we also discuss their relationship to so-called microstructure methods and other more advanced ideas such as time-dependent diffusivities.

研究动机与目标

• 基于连续性方程与广义通量，将广泛多样的扩散磁共振成像模型统一于单一理论框架之下。
• 通过强调其源自输运过程的推导，阐明Bloch-Torrey方程的数学与物理基础。
• 证明不同的扩散磁共振成像技术（例如，DTI、DKI、HARDI）对应于广义方程中特定的通量选择。
• 探索新的建模方向，包括相位敏感的信号分析、时变扩散系数，以及动量或能量的守恒。
• 挑战将实部与虚部信号组合为仅幅度数据的常规做法，建议保留相位与频率信息可提升模型精度。

提出的方法

• 从带有矢量通量项的连续性方程推导广义Bloch-Torrey方程，允许同时包含扩散与相干运动的贡献。
• 引入通量的一般形式：∇·(D∇M) + ∇·(uM)，其中D为扩散系数，u为速度，从而实现对多种自旋动力学的建模。
• 使用标准偏微分方程（PDE）技术求解广义方程，对磁化强度矢量M的空间与时间导数进行显式处理。
• 将既有的扩散磁共振成像模型（如DTI、DKI、HCP）映射到特定的通量假设上，证明其在广义框架下的等价性。
• 提出将框架扩展至包含动量守恒（类似Navier-Stokes方程）与能量守恒，以支持新型成像模态。
• 主张分别对磁共振信号的实部与虚部进行建模，以保留高斯噪声特性，避免Rician失真。

实验结果

研究问题

• RQ1 如何将广泛多样的扩散磁共振成像模型统一于单一数学框架之下？
• RQ2 广义Bloch-Torrey方程中通量项的物理与数学意义是什么？
• RQ3 不同的通量假设（如Fickian扩散、相干流动）如何对应于DTI或DKI等现有扩散磁共振成像技术？
• RQ4 对完整复信号（实部与虚部）进行建模而非仅幅度数据，其影响是什么？
• RQ5 广义框架是否能够支持基于动量或能量守恒的新型模型，用于磁共振信号形成？

主要发现

• 具有灵活通量项的广义Bloch-Torrey方程为扩散磁共振成像模型提供了统一框架，表明不同技术源于通量项的不同假设。
• 常见的扩散磁共振成像方法，如扩散张量成像（DTI）与扩散峰度成像（DKI），被证明对应于广义方程中通量项的特定选择。
• 该框架揭示了复磁共振信号中的相位与频率信息目前被严重低估，因为标准的仅幅度处理会引入Rician噪声并扭曲统计特性。
• 分别对实部与虚部信号进行建模可保留高斯噪声特性，可能提升拟合精度，但代价是信噪比降低。
• 该方法开启了新的建模方向，例如用于动量守恒系统的类似Navier-Stokes方程，以及将热力学与弛豫、扩散关联的能量守恒模型。
• 该框架支持对复杂微结构（包括具有非平凡拓扑结构的系统，如二元分支网络）的新型建模，可通过拉普拉斯算子的特征值分析实现。