跳转至

第1章：多物种行波与实际应用（Multi-Species Waves and Practical Applications）

§1 作者

  • 作者: J. D. Murray (James D. Murray),牛津大学数学生物学教授,Mathematical Biology 跨学科应用数学丛书的创始人之一。Murray 长期致力于将反应扩散方程 (reaction-diffusion equations) 应用于胚胎发育、肿瘤生长、伤口愈合、流行病空间传播等真实生物问题,以"用数学讲故事"的方式让应用数学家可以快速进入生物建模。
  • 章节角色: 本章是《Mathematical Biology II: Spatial Models》第三版第 1 章,定位为"单物种行波 (Vol I Ch 13)"到"多物种空间波现象"两卷之间的桥梁。Vol I 读者已熟悉 Fisher–Kolmogoroff 方程、spruce budworm 模型、单种群 logistic 行波,本章将这些概念扩展到多物种、趋化 (chemotaxis) 与可激发介质 (excitable media) 框架,并展示目标波 (target patterns)、螺旋波 (spiral waves)、追逐波 (pursuit-evasion waves) 等典型现象。
  • 章节风格: 先给出物理直觉 (§1.1),再分若干小专题 (捕食者-猎物、竞争扩散、GMO 扩散、BZ 反应、激发介质、振荡行波、螺旋波、混沌螺旋波) 提供模型与解析,每节都给出方程、线性化、本征值分析或数值结果,适合用作"案例研究"形式的进阶教材。
  • 与本书其他章节的关系: 本章建立的"线性化 → 本征值 → 最小波速"技术贯穿 Ch 2 (反应扩散图灵失稳)、Ch 5 (细菌趋化)、Ch 11 (脑肿瘤)、Ch 13 (传染病);"螺旋波"与"激发介质"是 Ch 8 (血管网络) 和 Ch 9-10 (伤口愈合) 中电生理耦合模型的前置知识。

§2 内容概述 (300-500 中文字符)

本章针对多物种反应扩散与趋化系统 (Vol I Eq 11.18) 给出系统化的行波 (travelling wave) 分析与应用示范。在 §1.1,作者用"局部振荡器作为起搏器 (pacemaker)"的直觉把目标波、螺旋波、阈值脉冲波 (threshold pulse) 三类典型模式统一到反应扩散框架下,并指出多物种情形的解空间远比单物种丰富:稳态之间可形成追逐波、振荡波前,在 2D/3D 中还能自发出现混沌螺旋波。

§1.2 给出修改型 Lotka-Volterra 捕食者-猎物系统,推导 \(c \geq [4a(1-b)]^{1/2}\) 的最小波速条件,并通过特征多项式分析证明存在临界 \(a^*\) 把波前解分为单调趋近和振荡趋近两类;随后讨论"对流型"追逐-逃避模型,把回避/追逐速度作为梯度的函数。§1.3-1.4 用两物种竞争模型 + Allee 效应分析英国灰松鼠扩散 (17.9 km²/yr) 与基因工程生物 (GEMs) 风险。§1.5 在 Oregonator / FKN 模型下推导 BZ 反应的最小波速,与实验 5.4 mm/min 对照。§1.6 转入激发介质 (excitable media) — FitzHugh-Nagumo 模型 — 解释 Hodgkin-Huxley 神经动作电位的传播;§1.7 给出极限环反应动力学下的行波列 (wave train) 速度;§1.8-1.9 详细分析 \(\lambda\)-\(\omega\) 系统的螺旋波与自发混沌螺旋波 (Winfree 1991; Agladze-Krinskii 1982) 的机制。

读本章需要 Vol I 第 1, 3, 11, 13 章作为前置:对非线性常微分方程的相平面分析、Liapunov 函数构造、反应扩散方程的线性化与本征值技术。

§3 核心方程与概念 (500-1000 中文字符,本卷最长的章节)

3.1 多物种反应扩散系统的统一形式

Vol I Eq 11.18 的多物种推广 (1.1):

\[ \frac{\partial \mathbf{u}}{\partial t} = \mathbf{f}(\mathbf{u}) + D \nabla^2 \mathbf{u}, \quad \mathbf{u} = (u_1, u_2, \ldots, u_n)^T \]

  • \(\mathbf{f}(\mathbf{u})\) 是反应动力学向量场;\(D\) 是扩散系数矩阵 (本章取对角常矩阵)。
  • 当 \(n=1\) 时退化为 Vol I Eq 13.4 (Fisher-KPP);\(n=2\) 时即可出现本章所有的复杂现象 (追逐波、螺旋波、混沌波前)。

3.2 追逐-逃避行波: 修改型 Lotka-Volterra (1.2)-(1.5)

模型 (1.2):

\[ \frac{\partial U}{\partial t} = A U \left(1 - \frac{U}{K}\right) - B U V + D_1 \nabla^2 U \]
\[ \frac{\partial V}{\partial t} = C U V - D V + D_2 \nabla^2 V \]

  • \(U\) 是猎物,\(V\) 是捕食者,\(A, B, C, D, K > 0\) 均为常数,\(D_1, D_2\) 为扩散系数。无量纲化 (1.3):

    \[ \frac{\partial u}{\partial t} = u(1 - u - v) + D \frac{\partial^2 u}{\partial x^2} \]
    \[ \frac{\partial v}{\partial t} = a v (u - b) + \frac{\partial^2 v}{\partial x^2} \]

  • (T) 理论;在 \((u,v) \in \mathbb{R}_+^2\) 中有三个稳态: \((0,0), (1,0), (b, 1-b)\) (假定 \(b < 1\) 即共存稳态存在)。共存态 \((b, 1-b)\) 是全局渐近稳定的,可由 Liapunov 函数证明 (1.4):

    \[ L(u, v) = a \left[ u - b - b \ln \frac{u}{b} \right] + \left[ v - (1-b) - (1-b) \ln \frac{v}{1-b} \right] \]
    \(L(b,1-b) = 0\), \(L > 0\) 其他位置, \(\frac{dL}{dt} \leq 0\)。

行波 ODE (1.5): 令 \(u(x,t) = U(z), v(x,t) = V(z), z = x + ct\) (波向左传播):

\[ c U' = U(1 - U - V) + D U'' \]
\[ c V' = a V (U - b) + V'' \]

在 \(D = 0\) 近似下化为三维动力系统 (1.6):

\[ U' = \frac{U(1 - U - V)}{c}, \quad V' = W, \quad W' = c W - a V (U - b) \]

最小波速 (1.10): 对稳态 \((1,0,0)\) 线性化,本征值为

\[ \lambda_1 = -\frac{1}{c}, \quad \lambda_{2,3} = \frac{1}{2}\left\{ c \pm \sqrt{c^2 - 4a(1-b)} \right\} \]
当 \(c < 2\sqrt{a(1-b)}\) 时,本征值是复数,系统在 \((1,0,0)\) 附近振荡发散 — 这意味着只有当 \(c \geq 2\sqrt{a(1-b)}\) 时才能存在连接 \((1,0,0)\) 与 \((b, 1-b, 0)\) 的非负行波前解,即最小波速
\[ c_{\min} = 2\sqrt{a(1-b)} \quad (b < 1) \]

振荡 vs 单调趋近 (1.11): 对 \((b,1-b,0)\) 线性化得特征多项式

\[ p(\lambda) \equiv \lambda^3 - \lambda^2 \left( c - \frac{b}{c} \right) - \frac{\lambda b + a b (1-b)}{c} = 0 \]
存在临界 \(a^*\): 当 \(a > a^*\) 时本征值是含负实部的复数 (即 \(U, V\) 以衰减振荡趋近共存态),\(a < a^*\) 时三个本征值皆为实数 (单调趋近)。这给出 §1.2 图 1.3 两种波前类型的物理解释。

3.3 对流型追逐-逃避 (1.12)-(1.13)

当扩散远小于个体主动运动 (如浮游生物、迁徙兽群) 时,模型写为守恒律形式:

\[ u_t - [(c_1 + h_1 v_x) u]_x = f(u, v) \]
\[ v_t - [(c_2 - h_2 u_x) v]_x = g(v, u) \]

  • \(c_1, c_2\) 是无干扰的游动速度;\(h_1, h_2 > 0\) 是"逃避-追逐"敏感度,\(h_1 v_x\) 意味着猎物背向捕食者密度梯度加速,\(h_2 u_x\) 意味着捕食者向猎物密度梯度加速。
  • 当 \(h_1 = h_2 = 0\) 时退化为平移坐标系下的反应动力学;当 \(h_1, h_2 > 0\) 时,需在特征线分析中考虑梯度项导致的双曲-反应耦合。

3.4 竞争模型与灰松鼠案例 (1.14 附近)

两物种竞争 + 扩散:

\[ \frac{\partial u}{\partial t} = r_u u(1 - u - \alpha v) + D_u \nabla^2 u \]
\[ \frac{\partial v}{\partial t} = r_v v(1 - v - \beta u) + D_v \nabla^2 v \]

  • 当种间竞争 \(\alpha \beta > 1\) 时存在竞争排斥;反之共存 (Vol I Ch 3 已分析)。
  • 行波解连接被排斥方 \((u=1,v=0)\) 与共存态;Skellam (1951) 与 Okubo (1980) 给出的实测波速公式 \(c = 2\sqrt{r_u D_u (1 - u^*)}\),其中 \(u^*\) 是共存态中竞争弱者的平衡密度,本章用英国灰松鼠替代红松鼠给出 \(c \approx 17.9\) km²/yr 的扩散速率。

3.5 BZ 反应的 Oregonator 模型 (§1.5)

FKN 三变量模型简化为 Oregonator 二变量 (无量纲):

\[ \frac{\partial u}{\partial t} = \frac{1}{\epsilon}\left[ u(1-u) - f v \frac{u - q}{u + q} \right] + \nabla^2 u \]
\[ \frac{\partial v}{\partial t} = u - v \]

  • \(u\) 是 HBrO₂ (自催化),\(v\) 是 Br⁻ (抑制剂),\(\epsilon \ll 1\) 是快慢时间尺度分离,\(f, q\) 是反应参数。
  • 通过行波前分析 (Murray §1.5 推导) 给出最小波速 \(c_{\min} = 1.0 \sim 2.0\) (无量纲),与实验测得的 \(5.4\) mm/min 在合理尺度变换下一致。

3.6 激发介质 (Excitable Media) 与 FHN 模型 (§1.6)

FitzHugh-Nagumo 简化模型:

\[ \frac{\partial u}{\partial t} = D_u \nabla^2 u + f(u) - v \]
\[ \frac{\partial v}{\partial t} = \epsilon (u - \gamma v) \]

  • \(u\) 是膜电位(快变量),\(v\) 是恢复变量(慢变量),\(f(u) = u(u-\alpha)(1-u)\) 是三次非线性,\(\epsilon \ll 1\)。
  • 系统具有阈值行为: 弱扰动 → 回到稳态;超过阈值的大扰动 → 长距离行波前 (类神经动作电位, Hodgkin-Huxley 的简化)。
  • 行波前解以 \(c = \sqrt{1/2}\) (无量纲) 的速度传播,这是 Hodgkin-Huxley 经典结果 \(c \approx 18.8\) m/s (枪乌贼大神经) 的简化对应。

3.7 \(\lambda\)-\(\omega\) 螺旋波 (§1.8-1.9)

复变量 Ginzburg-Landau 类模型 (极限环附近):

\[ \frac{\partial W}{\partial t} = W - (1 + i \alpha) |W|^2 W + (1 + i \beta) \nabla^2 W \]

  • \(W(x,t) = R(x,t) e^{i \phi(x,t)}\) 是复振幅;\(\alpha\) 是非线性频率修正,\(\beta\) 是扩散各向异性。
  • 当 \(1 + \alpha \beta < 0\) 时平面波失稳,出现螺旋波;螺旋臂数与 \(\alpha, \beta\) 的关系由 Kuramoto-Sivashinsky 推导给出。
  • 在 2D 模拟中,即使初始条件无混沌源,数值长时间演化也会出现螺旋波湍流 (Winfree 1991 称作 "Spiral Wave Break-up")。

§4 关键结论

  • 多物种反应扩散比单物种拥有显著更丰富的行波解空间: 单物种的 Fisher-KPP 解连接两个稳态;多物种可连接任意一对同维稳态,且可呈现追逐波 (pursuit-evasion)、阈值脉冲波 (threshold pulse) 与振荡波前 (oscillatory wavefront) 三种典型形态。
  • 最小波速公式对所有多物种波前具有普适结构: \(c_{\min} = 2\sqrt{D \cdot (\text{动力学线性增长率})}\)。本卷后续章节 (Ch 5 细菌、Ch 11 脑肿瘤、Ch 13 传染病) 都沿用此公式。
  • 临界参数 \(a^*\) 区分单调与振荡趋近: 在 Lotka-Volterra 追逐波中,\(a > a^*\) 时 \((u, v)\) 以衰减振荡接近共存态,\(a < a^*\) 时单调趋近 — 这一结果预示了"波前尾部的弱阻尼振荡"会导致与后继波前的相互作用 (Ch 11 脑肿瘤中的"波前-波尾耦合"现象的雏形)。
  • 目标波、螺旋波是同一起搏器机制的不同几何表现: 二维上,固定起搏器 → 圆形目标波;移动起搏器 (绕小圆环) → 螺旋波;两者的数学本质都是反应扩散 + 局部振荡器,差别仅在边界条件。
  • BZ 反应的理论最小波速与实验吻合: Oregonator 简化下 \(c_{\min} \sim 1.0\) (无量纲) 变换到物理量 \(\sim 5\) mm/min,与 Winfree (1974) 实验一致,验证了"两步快慢分离"的反应扩散机制可解释 BZ 反应的几何结构。
  • 激发介质的脉冲波速与神经动作电位速度处于同一量级: FHN 模型预测 \(\sim 1\) (无量纲) 对应 \(\sim 20\) m/s,接近 Hodgkin-Huxley 大神经实测值,说明 FHN 是 HH 的合理简化。
  • 螺旋波可在无混沌起搏器下自发湍流化: §1.9 给出数值证据 — 即使所有局部动力学是极限环 (非混沌),通过扩散耦合也可出现螺旋波破碎与时空混沌,这是 Belousov-Zhabotinskii 反应和心肌纤颤 (cardiac fibrillation) 中的关键机制。

§5 挑战和开放性问题

  • 多维最小波速的严格推导: 本章给出的 \(c_{\min}\) 都是 1D 结果;2D/3D 圆周或球面行波的最小波速是否存在类似 Fisher-KPP 公式,目前没有严格证明,是心电生理 (cardiac electrophysiology) 的重要开放问题。
  • 追逐波存在性的全局相图: 本章只在 \(D = 0\) 极限下分析,完整的 \(D_1/D_2\) 空间相图,特别是 \(D_1 \approx D_2\) 时的波前解结构,目前数值上有结果但解析相图尚未完成。
  • GMO 扩散的 Allee 效应下波速公式: §1.4 给出三稳态系统的定性图像,但具体的最小 Allee 阈值与扩散系数的关系缺乏闭式解,工程界 (EPA 风险评估) 急需这种公式。
  • BZ 反应螺旋波破碎的解析触发条件: §1.9 是数值实验,缺乏线性稳定性分析;Kuramoto (1984) 的弱湍流理论仅在 \(1 + \alpha \beta\) 略小于零时有效,强非线性区需要新工具。
  • 混沌螺旋波与心肌纤颤的连接: 临床心律失常中螺旋波破碎被认为是心室纤颤的起因,但 Murray 写本章时(2002 年)这一假说仍处早期,目前仍未被实验直接验证(2020 年代的研究进展包括 Fenton-Karma 模型,见 Qu et al. 2000 之后文献)。
  • 趋化 (chemotaxis) 与对流耦合的稳定性: §1.3 的对流模型假设了"无扩散 + 仅梯度驱动",实际生物系统 (黏菌、细菌) 同时具有扩散和趋化,Keller-Segel (Vol I §11.4) 与本章的耦合缺乏统一的稳定分析。

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

  • "追逐波"的命名其实有误导性: §1.2 给出的 Lotka-Volterra 波前解中,猎物并没有"逃避"动作 (v 始终是空间均匀的),"pursuit and evasion" 的命名暗含了行为学解释,但数学上只是两个标量场的耦合行波。Mimura-Murray 1991 之后的论文中已开始用 "Lotka-Volterra 竞争/捕食波前" 等中性名称。
  • §1.1 的"起搏器"图景在数学上不严谨: 把局部振荡器嵌入反应扩散系统会改变其动力学 (与外部介质耦合),所谓的"pacemaker"实际上是边界条件 (Neumann 边界 + 局部异质性) 的等效,这一区分在 Fenton 等 (2002) 心脏模型中变得关键。
  • \(c_{\min} = 2\sqrt{a(1-b)}\) 与极限环波速混淆的风险: Vol I §13.2 的 Fisher-KPP 最小波速 \(c_{\min} = 2\sqrt{Dr}\) 与本章追逐波最小波速形式相同但物理意义不同 — 后者来自"两个稳态 + 振荡发散"的本征值,前者来自"实正本征值"。Murray 在 §1.2 没有明确指出这一结构对应,读者易误读为"任意反应扩散行波都有相同 \(c_{\min}\) 公式"。
  • GMO 扩散模型过于简单: §1.4 假设了均匀环境 + Allee 效应,实际 GMO 释放需要考虑异质地形、季节性繁殖、基因流强度、监测策略 (如 sterile insect technique),这些因素都让 \(c_{\min}\) 公式失效 — 监管机构使用本章模型做风险评估时,需注意模型局限。
  • 可借鉴/可独立复现的推导:
  • §1.2 的 Lotka-Volterra 行波 ODE 化为三维动力系统 + 本征值分析 (1.6) → (1.9) → (1.10) — 经典且自包含,适合作为反应扩散行波最小波速的"教学例题"。
  • §1.2 末尾的"特征多项式 → 局部极值 → 临界 \(a^*\)"是动力系统分岔的极好范例,可在 FEniCS / Dedalus 中数值复现 (Vol I Ch 13 的复现脚本是好的起点)。
  • 如果要问作者一个问题: 我会问"为什么 Vol I 的 \(c_{\min}\) 公式与本章追逐波的 \(c_{\min}\) 形式相同但对应不同的本征值结构?是反应扩散行波的某种通用性质,还是这两个具体例子的巧合?" — 这个区分对推广到多物种多稳态情形至关重要。
  • 跨章节连结:
  • Ch 2 (Reaction Diffusion Pattern) 的 Turing 失稳是本章"非行波"的另一面 — 本章分析行波沿 1D 传播,Ch 2 分析稳态空间结构。
  • Ch 5 (Bacterial Patterns and Chemotaxis) 给出真实的趋化数据 (Tyson et al. 1998-1999) 验证了 §1.3 的趋化波速公式。
  • Ch 9-10 (Wound Healing) 借用本章"激发介质"框架,把表皮/真皮伤口的电生理-力学耦合作为 FHN 模型的扩展。
  • Ch 13 (Geographic Spread of Epidemics) 借用 §1.3 的竞争波速公式,加上非线性传染率,得到 rabies 的空间波前速度 (与本章追逐波形式上同构)。

§7 重要参考文献

  • [X1] Murray J.D. Mathematical Biology II: Spatial Models and Biomedical Applications (3rd ed.), Springer, 2002 (本书)
  • [X2] Dunbar S.R. "Travelling wave solutions of diffusive Lotka-Volterra equations: a heteroclinic connection in \(\mathbb{R}^4\)", Trans. AMS 1983;269:557-570
  • [X3] Dunbar S.R. "Travelling wave solutions of diffusive Lotka-Volterra equations: a heteroclinic connection in \(\mathbb{R}^4\) (具体 predator-prey pursuit wave 分析)", J. Math. Biol. 1984;17:11-32
  • [X4] Tyson J.J., Fife P.C. "Target patterns in a realistic model of the Belousov-Zhabotinskii reaction", J. Chem. Phys. 1980;73:2224-2237
  • [X5] Winfree A.T. "Rotating chemical reactions", Scientific American 1974;230(6):82-95
  • [X6] Zaikin A.N., Zhabotinskii A.M. "Concentration wave propagation in two-dimensional liquid-phase self-oscillating system", Nature 1970;225:535-537
  • [X7] Agladze K.I., Krinskii V.I. "Spiral waves in active media: chronotopography of rotating waves", Nature 1982;296:424-426
  • [X8] Okubo A., Maini P.K., Williamson M.H., Murray J.D. "On the spatial spread of the grey squirrel in Britain", Proc. R. Soc. Lond. B 1989;238:113-125
  • [X9] Skellam J.G. "Random dispersal in theoretical populations", Biometrika 1951;38:196-218
  • [X10] Keller E.F., Segel L.A. "Model for chemotaxis", J. Theor. Biol. 1971;30:225-234
  • [X11] Sherratt J.A., Lewis M.A., Fowler A.C. "Deterministic chaos in a model of microbial population dynamics", Proc. R. Soc. Lond. B 1995;259:25-32
  • [X12] Kuramoto Y. Chemical Oscillations, Waves, and Turbulence, Springer, 1984 (DOI: 10.1007/978-3-642-69689-3)
  • [X13] Epstein I.R., Showalter K. "Nonlinear chemical dynamics: oscillations, patterns, and chaos", J. Phys. Chem. 1996;100:13132-13147
  • [X14] Merkin J.H., Petrov V., Scott S.K., Showalter K. "Wave-induced chaos in a continuously fed unstirred reactor", J. Chem. Soc. Faraday Trans. 1996;92:2911-2915
  • [X15] Petrov V., Ouyang Q., Swinney H.L. "Resonant generation of spatial patterns in a non-equilibrium system", Phys. Rev. Lett. 1994;83:2893-2896