Chapter 2: Fractional Models in Biology and Medicine

1. 作者

Kevin Burrage（通信作者），澳大利亚 Queensland University of Technology (QUT) 数学科学学院教授，Email: kevin.burrage@qut.edu.au
Pamela M. Burrage，QUT 数学科学学院，Email: pamela.burrage@qut.edu.au
Alfonso Bueno-Orovio，英国牛津大学（University of Oxford）计算机科学系，Email: alfonso.bueno@cs.ox.ac.uk

三人在分数阶数值方法、计算心脏电生理学、随机建模与化学动力学领域有长期合作。Bueno-Orovio 在分数阶扩散方程求解（特别是用于分数阶 Laplacian 的矩阵转移技术）方面有奠基性贡献。

2. 内容概述

本章是全书唯一聚焦生物医学应用的章节，系统讨论分数阶 PDE 在 4 个应用中的建模与计算：(1) 分数阶反应扩散（空间/时间分数阶）：涵盖 Allen-Cahn、Fisher-KPP、FitzHugh-Nagumo、Gray-Scott 等模式形成模型；(2) 分数阶化学动力学：在异质（拥挤）环境中的随机模拟——给出 ASSA 算法、分数阶化学主方程（FCME）的 Mittag-Leffler 解及 Cauchy 轮廓法；(3) 心脏电生理异质组织建模：通过分数阶 monodomain 方程捕捉多尺度微结构；(4) 分数阶心脏模型的 DTI 校准：基于修正的分数阶 Bloch-Torrey 方程，引入"stretched Mittag-Leffler 模型"从扩散 MRI 数据反演分数阶指数 $\alpha, \beta$。

3. 核心方程与概念

3.1 分数阶反应扩散方程

一般形式（时间-空间分数阶，单变量）： $$\frac{\partial^\beta u}{\partial t^\beta} = -D(-\Delta)^{\alpha/2}u + f(t,u)$$ 其中时间分数阶用 Caputo 导数、空间分数阶用 Riesz 形式的 $(-\Delta)^{\alpha/2}$。当 $D_1=D_2=D$、$\beta=1$、$\alpha=2$ 时退化为标准反应扩散。

特殊应用： - 珊瑚生长模型（Somathilake & Burrage 2018, [X1]）：$f = 1-u-\theta_1 u^2 v$，$g = -\theta_2 v + \theta_1 uv^2$，描述营养物与碳酸钙之间的生物量转化。 - Allen-Cahn（相场、合金动力学）：$f = u-u^3$。分数阶扩散下界面增厚且不稳定界面持续时间更长 [X2]。 - 分数阶 Fisher-KPP：Engler [X3] 证明波前传播速度无界，被 Burrage et al. [X2] 数值验证。 - FitzHugh-Nagumo（神经与心脏电传导）：$D_2 = 0$，$f = u(1-u)(1-a)-v$，$g = \epsilon(bu-\gamma v-\delta)$。随 $\alpha$ 减小，传播速度与波长都下降 [X2]。 - 分数阶 Gray-Scott（形态发生）：随 $\alpha$ 减小，模式形成从曲率驱动演化为核化（边界处小结构）。

Turing 不稳定性：通过线性化得到关于 $k^2$ 的二次色散方程，系数为 $\alpha, D_1, D_2$ 的非线性函数。分数阶设定下 Turing 条件需重新分析。

边界条件难题：分数阶 Laplacian 是非局部算子，定义在无穷域上。处理有限域时需谨慎。Ilić 矩阵转移技术（无证明）和 Cusimano 反射矩阵（已证明，1D）：在 $[0,L]$ 上，$N+1$ 阶反射矩阵的本征值为 $$\frac{2}{h}\left(\frac{\sin\frac{(j-1)\pi}{2N}}{\cos\frac{\pi}{2\alpha}}\right)\cos\left(\frac{2(j-1)+\alpha(N+1-j)\pi}{2N}\right)$$ 当 $h\to 0$ 时收敛到连续算子 $(j-1)\pi/L$ 的 $\alpha$ 次幂 [X4]。

时间分数阶的记忆问题：Lubich 高阶近似需要追溯到初始条件（O($h^p$) 收敛）。可用时间窗或分级网格缓解。

记忆归属争议：Caputo 导数作用于 LHS（Lindenberg & Yuste 视角，反应项也"记旧"）vs 仅作用于扩散项（Sung et al. 视角）——取决于具体应用场景。

3.2 分数阶化学动力学

均相环境（良好搅拌）：标准 SSA（随机模拟算法）→ CLE（化学 Langevin 方程）→ ODE 三个层次。

异质（拥挤）环境：平均平方位移 $\sim t^\alpha$，$\alpha \in (0,1]$。细胞质膜上 40% 障碍密度对应 $\alpha \approx 0.69$ [X5]。

分数阶化学主方程（FCME）： $$_0D_t^\alpha p(t) = Ap(t), \quad p(0) = p_0$$ 解为 Mittag-Leffler 矩阵函数： $$p(t) = E_\alpha(At^\alpha)p_0$$ 其中 $E_\alpha(z) = \sum_{j=0}^\infty \frac{z^j}{\Gamma(\alpha j+1)}$。

计算难题：直接级数展开在 $\alpha$ 远离 1 时收敛慢。Weideman-Trefethen 技巧 [X6]：用 Bromwich 围道积分和抛物线/双曲线轮廓，避开 $A$ 本征值对应的支点奇点。梯形规则配合 $e^{zt}$ 指数衰减因子实现快速收敛。

ASSA 算法（MacNamara, Henry, McLean 2017 [X7]）：把 SSA 的指数等待时间替换为重尾的 Mittag-Leffler 等待时间 $$\tau = -\frac{1}{\lambda}\left(\frac{\sin(\pi\alpha)}{\tan(\pi\alpha(1-u_1))} - \cos(\pi\alpha)\ln(u_2)\right)$$ 当 $\alpha=1$ 时退化为标准指数等待时间。

特殊分数阶 $\alpha = (1/2)^{k+1}$ 的优雅表示（Orsingher-Beghin [X8]）：时间分数阶扩散方程的解是迭代布朗运动（iterated Brownian motion）的密度 $$Z(t) = B_1(|B_2(\cdots|B_{k+1}(t)|\cdots)|), \quad B_j \sim N(0,2)$$ 对 $\alpha = 1/4$，$Z(t) = B_1(|B_2(t)|)$。对一般 $\alpha$ 也有"时间变换"解（虽然不那么简洁）。

单分子反应的精确 CME 解（Jahnke-Huisinga [X9]）：Birth-Death 过程（eq. 19）的解是 Binomial 与 Poisson 的卷积 (eq. 20) $$p(t,x) = \sum_{k=0}^{\min(0,X_0)} \binom{X_0}{k} q(t)^k(1-q(t))^{X_0-k} \frac{\lambda(t)^{X-k}e^{-\lambda(t)}}{(X-k)!}$$ 其中 $q(t) = e^{-k_2 t}$，$\lambda(t) = (k_1/k_2)(1-e^{-k_2 t})$。该结构是分数阶情形下用迭代布朗路径算法生成样本的基础。

3.3 心脏组织电传导的分数阶建模

标准 monodomain 方程（电缆方程同质化）： $$\nabla\cdot(\sigma\nabla V_m) = \chi(C_m\partial_t V_m + I_{ion} - I_{stim})$$ 耦合细胞内 ODE（Beeler-Reuter 1977 [X10] / Tomek 等人的人心室 AP 模型 [X11]）。

标准 bidomain 方程（更精细的胞内/胞外双空间）： $$\nabla\cdot(\sigma_i\nabla\phi_i) = \chi(C_m\partial_t V_m + I_{ion} - I_{stim,i}) \text{ in } \Omega_i$$ $$\nabla\cdot(\sigma_e\nabla\phi_e) = -\chi(C_m\partial_t V_m + I_{ion} - I_{stim,e}) \text{ in } \Omega_e$$

异质性问题：心脏组织有细胞-纤维-薄片-结缔组织多尺度结构，标准方程 (23)(25) 难以捕捉。

分数阶 monodomain（Bueno-Orovio et al. 2014 [X2] 的核心创新）： $$\partial_t V_m = -D_\alpha(-\Delta)^{\alpha/2}V_m - \frac{1}{C_m}(I_{ion}(V_m, z) - I_{stim}), \quad \alpha \in (1,2]$$ $$\frac{dz}{dt} = f(V_m, z)$$

生物物理基础——Riesz 位势理论：分数阶 Laplacian 在 $\mathbb{R}^N$ 上由 Riesz 位势（[X12] M. Riesz, 1949）定义为 $$(-\Delta)^{\alpha/2}\phi(\mathbf{r}) = C_{N,\alpha}\int_{\mathbb{R}^N}\frac{\phi(\mathbf{r}) - \phi(\mathbf{s})}{\|\mathbf{r}-\mathbf{s}\|^{N+\alpha}}d\mathbf{s}, \quad 1<\alpha\le 2$$ 归一化常数 $C_{N,\alpha} = \pi^{-N/2}2^{2\alpha}\Gamma(\alpha/2)/\Gamma((N-\alpha)/2)$（本章已修正 [X2] 中的笔误）。

关键对应（[X2] eq. 27, 28, 30, 31 推导）： - $N=3, \alpha=2$：$C_{3,2}=4\pi$，$\phi(r) \sim 1/r$（单极子） - $N=3, \alpha=1$：$C_{3,1}=2\pi^2$，$\phi(r) \sim 1/r^2$（偶极子） - $N=2, \alpha=1$：$C_{2,1}=2\pi$，$\phi(r) \sim 1/r$（2D 偶极子） - $N=2, \alpha\to 2$：$C_{2,2}=4\pi$，$\phi(r) \sim \log r$（2D 单极子）

物理解释：心脏中标准方程假设均匀传导率 → 主导单极子场。但实际组织中不同电导率区域的界面产生次级偶极子场 $$\phi(\mathbf{r}) = \phi_m(\mathbf{r}) + \phi_p(\mathbf{r})$$ 即分数阶 Laplacian 正是描述单极子→偶极子过渡的严格数学工具。$\alpha$ 从 1 到 2 的变化对应电场多极展开的截断——这是文章最重要的物理洞察。

数值结果（[X2] 用 Beeler-Reuter 耦合，1D 分数阶 monodomain）：$\alpha$ 减小导致去极化峰值降低、AP 脚部变宽、复极早期更显著、传导速度下降。

3.4 心脏分数阶模型的 DTI 校准

标准 Bloch-Torrey 方程（磁化水分子扩散信号）： $$\partial_t S = -i\gamma\mathbf{r}\cdot G(t)S + \nabla\cdot D\nabla S$$ Stejskal-Tanner 序列下信号为单指数衰减 $S = e^{-b\mathbf{g}^T D\mathbf{g}}S_0$，$b = (\gamma|G|\delta)^2(\Delta - \delta/3)$。

问题：异质组织中信号不是单指数——双指数、拉伸指数经验修正都无法严格推导。

分数阶 Bloch-Torrey 方程： $$\tau_\beta\,_0D_t^\beta S = -i\gamma\mathbf{r}\cdot G(t)S - \tau_\alpha^{\alpha-1}(-\nabla\cdot D\nabla)^\alpha S$$ $0 < \alpha, \beta \le 1$，$\tau_\alpha, \tau_\beta$ 为时间单位常数。

Stejskal-Tanner 解析（$\beta=1$）（[X13]）： $$S = e^{-\tilde b(g^TDg)^\alpha}S_0, \quad \tilde b = \tau_\alpha^{\alpha-1}(\Delta - \frac{\delta(2\alpha-1)}{2\alpha+1})(\gamma\delta G)^{2\alpha}$$ $\beta \ne 1$ 情形仍为开放问题。

Stretched Mittag-Leffler 模型（Calvetti et al. [X14] 的"ansatz"）： $$S = E_\beta(-(b\cdot g^TDg)^\alpha)S_0$$ 不严格解方程 (35)，但能同时拟合超扩散与亚扩散分量。实验（4 只健康大鼠离体心脏，多 $b$ 值扫描）：在左心室拟合 $b$ 高达 10000 s/mm² 数据，结果显示$\alpha \approx 1$（即标准空间扩散），但$\beta$ 呈组织/个体间分布（亚扩散）。这意味着心脏异质性主要来自时间分数阶（亚扩散）而非空间。

方法论要点：不存在捕获心脏异质性的"统一 $\alpha, \beta$"——存在跨组织/跨动物的分布。这与模型群体方法（[X15]）一致。

4. 关键结论

分数阶反应扩散模型中，$\alpha$ 减小 → 模式形成更慢、波长变短、边界核化结构出现——但"分数阶"不是单纯的调控旋钮，而是捕捉了非局部相互作用。
分数阶化学动力学的关键是时间变换（$\alpha$ 远离 1 时，时间尺度被显著放慢）。FCME 解的 Mittag-Leffler 矩阵函数可用 Cauchy 轮廓法快速计算。
心脏分数阶模型的最强论据是 Riesz 位势理论：分数阶 Laplacian 描述从单极子到偶极子的过渡——这恰好是异质组织中电场多极展开的物理结构。
DTI 校准显示：心脏异质性主要表现为时间亚扩散（$\beta < 1$），空间分数阶指数 $\alpha \approx 1$。这与脑组织（$\alpha$ 偏离 1 明显）形成对比。
整个章节传递的元结论：分数阶模型的使用需要严格的"为什么"论证——不能简单地"用分数阶算子替代标准算子"。

5. 挑战和开放性问题

异质界面的处理：当 $\alpha$ 在空间上分段变化（如 [0,L] 上 $\alpha_1$，(L,2L] 上 $\alpha_2$）时，无法用齐次 Neumann 边界条件 + 界面反射来定义正确的非局部算子。Cusimano [X4] 用 $(-\Delta)^{\alpha(x)/2}$ 在整个域上 + 边界反射，但这只是"generalisation of the dynamics"，缺乏严格证明。
$\beta \ne 1$ 的分数阶 Bloch-Torrey 方程目前没有解析解，亦无高效数值方法。这是 DTI 校准的最大瓶颈。
Turing 不稳定性的分数阶判定准则尚未在 [X2] 中给出完整分析——只是说"在 SFRDE 设定下需要重新分析"。
复杂 2D 心脏几何 + 分数阶的数值实现仍非常昂贵，缺乏生产级求解器。
分数阶化学动力学的非线性反应目前只覆盖了单分子反应，双分子及更复杂反应的分数阶 CME 仍未解决。
心脏分数阶模型对复极化离散（dispersion of repolarisation）的预测需要更大的临床数据集验证。
$\alpha$ 的"实际意义"：作者在心脏部分用 Riesz 位势理论给出的物理解释漂亮，但仍是一个 ansatz——缺乏"分数阶 $\alpha$ 直接对应某种具体亚细胞结构"的定量映射。

6. 个人反思与批判性分析

优点： 1. 本章是全书中生物物理基础最扎实的一章。Riesz 位势 → 单极子/偶极子的论证是真正把"分数阶"还原为有物理解释的多极展开——这在其他 9 章中很少见。 2. 章节结构逻辑链非常完整：模式形成 → 化学动力学 → 心脏电传导 → DTI 校准，每一节都用具体生物学现象和真实实验数据锚定，不是空谈理论。 3. 数值方法部分实用性高：Cauchy 轮廓法、矩阵转移技术、Stretched Mittag-Leffler 模型都给出了可操作的具体公式（如 eq. 16, 37, 38），读者可直接复现。 4. 对边界条件和记忆效应两大分数阶 PDE 难题的处理非常坦诚——明确指出 Ilić 矩阵转移技术"无证明"，Cusimano 用主方程 + 反射矩阵补充了证明。

缺点： 1. 章节的心脏部分偏重 Bueno-Orovio 自己组的工作（[X2] 出现频次过高），对其他小组（如 Qu, Garfinkel, Fenton 等在心脏分数阶建模方面的工作）覆盖不足。 2. 对分数阶算子的生物物理等价物的论证虽然漂亮，但 Riesz 位势与细胞异质结构的定量映射未给出。例如，40% 障碍密度 → $\alpha = 0.69$ [X5] 这种具体标定关系仅在细胞质膜上验证，未扩展到组织尺度。 3. DTI 校准结果显示 $\alpha \approx 1$（标准空间扩散），这其实对心脏分数阶模型的核心论据是挑战——如果空间分数阶不显著，那 26 式相对标准方程 (23) 的"改进"到底源自 $\alpha$ 的偏离还是其他因素？作者没有正面回答。 4. 数值方法部分虽然实用，但缺少对实际计算成本的定量分析。例如 Cauchy 轮廓法在 2D/3D 心脏尺度上是否可行？$b$-value 反演的统计鲁棒性如何？ 5. 章节没有给出对所有分数阶模型都成立的一般性局限讨论——例如分数阶 PDE 的适定性（well-posedness）在非线性情形下不平凡，但对应用者至关重要。

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - 本章给出的"分数阶 $\alpha$ 标定异质性"思路，可推广到血管壁的层状结构（内弹力板、中膜、外膜）——不同层的弹性 / 粘弹性 / 渗透性差异本质上是非局部的，分数阶算子是天然选择。 - 异质扩散 → 模式形成（Turing 不稳定性）的思路可应用于动脉粥样硬化斑块的形成——低密度脂蛋白在内膜下层的异常聚集可能就是"分数阶 Fisher-KPP 方程"驱动的 [X3]。 - ASSA 算法的 Mittag-Leffler 等待时间思想可推广到细胞外基质重塑的动力学随机模拟——基质金属蛋白酶（MMP）的激活涉及"长等待 + 短爆发"模式，与重尾分布吻合。 - [X2] 中"分数阶 $\alpha$ 减小 → 模式形成更慢、波长变短"的结论，与本人 SMC 培养实验中观察到的"细胞密度波传播速度对底物刚度敏感"现象思路一致，可能需要类似分数阶反应扩散框架来建模。

7. 重要参考文献（按出现顺序编号）

[X1] L.W. Somathilake, K. Burrage. Cogent Math. Stat. 5, 1426524 (2018). —— 分数阶珊瑚生长模型。
[X2] A. Bueno-Orovio, D. Kay, K. Burrage. BIT Numer. Math. 54, 937 (2014). —— 本章核心：分数阶心脏 monodomain 模型及 Riesz 位势论证。
[X3] H. Engler. Int. J. Differ. Equ. 2010, 315421 (2010). —— 分数阶 Fisher 方程波前速度无界性证明。
[X4] N. Cusimano, K. Burrage, I. Turner, D. Kay. App. Math. Modell. 42, 554 (2017). —— 反射矩阵技术。
[X5] D.V. Nicolau Jr., J. Hancock, K. Burrage. Biophys. J. 92, 1975 (2007). —— 40% 障碍 → $\alpha \approx 0.69$ 的实验。
[X6] J.A.C. Weideman, L.N. Trefethen. Math. Comput. 76, 1341 (2007). —— Mittag-Leffler 矩阵函数的 Cauchy 轮廓法。
[X7] S. MacNamara, B. Henry, W. McLean. SIAM J. Appl. Math. 77, 447 (2017). —— ASSA 算法。
[X8] E. Orsingher, L. Beghin. Ann. Probab. 37, 206 (2009). —— 迭代布朗运动与时间分数阶扩散。
[X9] T. Jahnke, W. Huisinga. J. Math. Biol. 54, 1 (2007). —— 单分子反应 CME 精确解。
[X10] G.W. Beeler, H. Reuter. J. Physiol. 268(1), 177 (1977). —— 心室肌细胞经典模型。
[X11] J. Tomek et al. (人心室 AP 模型). —— 现代人心室细胞模型。
[X12] M. Riesz. Acta Math. 81, 1 (1949). —— Riesz 位势理论原始论文。
[X13] （分数阶 Bloch-Torrey 方程 Stejskal-Tanner 解析）. —— 分数阶 DTI 解析框架。
[X14] D. Calvetti et al. —— Stretched Mittag-Leffler 模型及 4 只大鼠实验数据。
[X15] （模型群体方法）. —— 心脏异质性统计建模。
[X16] D.T. Gillespie. J. Comput. Phys. 22, 403 (1976). —— SSA 原始论文。
[X17] R. FitzHugh. Biophys. J. 1, 445 (1961). —— 神经动作电位 FHN 模型。
[X18] J.D. Murray. Mathematical Biology II (Springer, 3rd ed., 2003). —— 反应扩散 Turing 不稳定性经典教材。
[X19] A.M. Turing. Philos. Trans. Royal Soc. London, Ser. B 237, 37 (1952). —— Turing 不稳定性原始论文。
[X20] R. Metzler, J. Klafter. Phys. Rep. 339, 1 (2000). —— 分数阶反常扩散的综述。

本章由主 agent 亲自从 PDF 提取并撰写，无子任务委托。
撰写日期：2026-06-01