跳转至

第十四章：呼吸读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd。本章把 Vol I Ch 13（血液）和 Ch 16（内分泌反馈）的工具，扩展到呼吸系统——气体在肺泡-毛细血管界面的扩散、CO2/HCO3⁻ 化学对气体运输的放大效应、通气/灌注比、Cheyne-Stokes 周期性呼吸的延迟微分方程模型、呼吸中枢的互抑制神经网络。本章的数学核心是Mackey-Glass 型延迟微分方程（14.57），它给出 Cheyne-Stokes 呼吸（40-60 秒周期）的 Hopf 分岔机制，是临床心衰合并陈-施呼吸的经典数学解释。

内容概述

本章分为 4 大主题：(1) 14.1 肺泡-毛细血管气体交换——Fick 扩散定律、毛细血管内气体的稳态分布、CO2/碳酸酐酶对 CO2 运输的增强（因子 \(1 + K_c\)）、Hb 对 O2 运输的增强（因子可达 200，真实约 32）、CO 中毒数学（0.1% CO 致死）；(2) 14.2 通气与灌注——通气/灌注比 \(\dot{V}/Q\)、BTPS vs STPD 单位转换、O2-CO2 图、呼吸交换比 R；(3) 14.3 通气调节——Mackey-Glass 延迟微分方程 (14.57) 解释 Cheyne-Stokes 呼吸（心衰和脑损伤的标志）、Grodins 多室模型（脑-组织-肺-血液的延迟耦合）；(4) 14.4 呼吸中枢——pre-Bötzinger 复合体（preBötC）、mutual inhibition 神经网络（von Euler 1977/Wyman 1977 互抑制模型，14.81-14.82）、Hering-Breuer 反射。

本章最具数学深度的部分是 14.3 节 —— Mackey-Glass 方程 (14.57) \(dx/dt = \lambda - \alpha x \dot{V}(x(t-\tau))\)，其中 \(\dot{V}(x) = V_m x^n/(\theta^n + x^n)\)（Hill 函数模拟通气对 PCO2 的响应）。线性化分析 (14.60-14.65) 给出 Hopf 分岔条件 \(\sigma > \sigma_c\)，对应 Cheyne-Stokes 呼吸——这是临床心衰合并陈-施呼吸的直接数学解释。

核心方程与概念

14.1 肺泡-毛细血管气体交换

Fick 扩散（14.2）\(q = D_s(P_g - U/\sigma)\)。质量守恒（14.4）\(U_t + (vU)_x = pq/A\)。稳态解（14.7）\(U(x) = \sigma P_g + (U_0 - \sigma P_g) e^{-D_m x/v}\)，其中 \(D_m = \chi D_s/\sigma\) 是膜交换速率倒数，\(\chi = p/A\) 是表面积比。总通量（14.8）\(Q = vA\sigma(P_g - P_0)(1 - e^{-D_m L/v})\)。长毛细血管极限（14.10）\(Q \to vA\sigma(P_g - P_0)\)——流量有限但取决于压差。

14.1.3 CO2 与碳酸酐酶的化学放大

CO2 ⇌ H+ + HCO3⁻ 准稳态（14.14）：\(V = K_c U\)，\(K_c = k_1/(k_{-1}[H+])\)。总通量（14.15-14.16）\(Q = vA(1+K_c)\sigma_{CO_2}(P_0 - P_{CO_2})(1 - e^{-D_{CO_2}L/[v(1+K_c)]})\)——碳酸酐酶把有效流速放大 \(1 + K_c\) 倍。pH = 7.4 时 \(K_c \approx 20\)（\(\log_{10}(k_1/k_{-1}) = -6.1\)），故 CO2 运输被放大 ~21 倍。

14.1.4 血红蛋白对 O2 的放大

协同结合（14.20-14.23）\(v \cdot d/dx(W + 4Z_0 f(W)) = D_{O_2}(\sigma_{O_2} P_{O_2} - W)\)，其中 \(f(W) = W^4/(K_4 + W^4)\) 是 O2 饱和函数。等价形式（14.26）\(v(1 + 4Z_0 f'(W)) dW/dx = D_{O_2}(\sigma_{O_2} P_{O_2} - W)\)——Hb 把有效流速放大 \(1 + 4Z_0 f'(W)\) 倍。

最大放大倍数：\(f'(W)\) 在 \(W = K_4\) 处取最大值 \(1/K_4\)，\(K_4 = 30\sigma_{O_2} \approx 4.2 \times 10^{-2}\) mM，\(Z_0 \approx 2\) mM，\(4Z_0 f'_{max} = 8/(4.2 \times 10^{-2}) \approx 200\)。实际增强：100 ml 血载 20 ml O2（Hb 携带 97%），无 Hb 仅 0.5 ml，实际增强 ~40 倍。

14.1.5 CO 中毒

CO 竞争方程（14.30-14.34）\(\beta = K_{O_2}/K_{CO} \approx 200\)。临界 CO 浓度（14.35）\(u > 4.64\) → 缺氧，等价于 \(P_{CO} > 0.7\) mmHg ≈ 0.1% 大气浓度——这是 CO 中毒致死阈值的数学解释。

治疗原理（14.38）：高压氧（2 atm）使 \(w^*\) 从 3.5 升至 50.7，CO 清除速率增加 ~14 倍。

14.2 通气与灌注

通气/灌注比 \(\dot{V}/Q\) 是肺血液气体含量的决定因素。通气 \(\dot{V} \approx 5250\) ml/min（500 ml × 15 bpm - 150 ml 解剖死腔），灌注 \(Q \approx 5000\) ml/min，正常比值 \(\dot{V}/Q \approx 1\)。

稳态气体交换（14.41）\(\dot{V}/Q = (c_v - c_a)RT/(P_a - P_i)\)。单位换算（14.44）\(V_{STPD}/V_{BTPS} = (P-47)/863\)。O2-CO2 图（Fig. 14.9）：所有可能的 \((P_{a,O_2}, P_{a,CO_2})\) 组合在同一条曲线上，由 \(\dot{V}/Q\) 参数化。呼吸交换比 \(R = \dot{V}_{CO_2}/\dot{V}_{O_2}\)：碳水化合物代谢 RQ = 1.0、脂肪代谢 RQ = 0.7、混合饮食 RQ = 0.825。

通气/血酸碱： - 过度通气（\(\dot{V}/Q\) ↑）→ \(P_{a,CO_2}\) ↓ → 呼吸性碱中毒 - 通气不足（\(\dot{V}/Q\) ↓）→ \(P_{a,CO_2}\) ↑ → 呼吸性酸中毒

14.3 通气调节 — Cheyne-Stokes 呼吸

14.3 Mackey-Glass 模型

延迟微分方程（14.57）\(dx/dt = \lambda - \alpha x \dot{V}(x(t - \tau))\)，其中 \(\dot{V}(x) = V_m x^n/(\theta^n + x^n)\)。Table 14.1 参数：\(\lambda = 6\) mmHg/min、\(V_m = 80\) L/min、\(\tau = 0.25\) min。无量纲化（14.58）\(dy/ds = \beta - y F(y(s - \sigma))\)。

线性化（14.59-14.60）特征方程 \(\mu + F(y^*) + y^* F'(y^*) e^{-\mu \sigma} = 0\)。关键函数（14.61）\(g(y) = F(y) - y F'(y)\)。Hopf 分岔条件（14.64-14.65）\(\tan \omega \sigma = -\omega/F(y^*)\)，\(\sigma = (1/\omega)[\pi + \tan^{-1}(-\omega/F(y^*))]\)。

临床意义：如果 \(\sigma > \sigma_c\)（临界延迟），稳态失稳，出现 Cheyne-Stokes 呼吸——深呼吸与呼吸暂停交替，周期 40-60 秒。Lange & Hecht 1962 9 例心衰患者中 7 例在 2 年内死亡——Cheyne-Stokes 呼吸是心衰预后不良的独立预测因子。

两种诱发机制：(1) 延迟增大（心衰导致脑血流时间延长）；(2) 反馈增益 \(F'\) 增大（脑损伤）。

14.3.1 Grodins 多室模型

四室（脑-组织-肺-血液）耦合延迟 ODE（14.66-14.77）。关键参数（Table 14.2）：\(V_B = 1\) L、\(V_T = 39\) L、\(V_L = 3\) L、\(Q_B = 0.75\) L/min、\(Q = 6\) L/min、\(\tau_{aB} = 0.18\) min、\(\tau_{aT} = 0.32\) min。通气（14.78）\(\dot{V} = \dot{V}_{base} + G_C \max(P_B - I_C, 0)\)。Hopf 分岔 \(G_C = 1.8\)（稳定）→ \(G_C = 9\)（振荡，周期 ~60 s）。

14.4 呼吸中枢

pre-Bötzinger 复合体（Smith 1991 发现）：节律生成的必要充分条件。互抑制网络（von Euler 1977, Wyman 1977）（14.81-14.82）\(\tau dI_1/dt + I_1 = F(E_1 - I_2)\)，\(\tau dI_2/dt + I_2 = F(E_2 - I_1)\)，对称耦合，\(F(x) = 0\)（\(x < 0\)），\(F\) 单调增（\(x > 0\)）。

Hering-Breuer 反射：肺过度充气时，支气管壁牵张感受器激活 → 切断吸气斜坡。

关键结论

  1. CO2 运输被碳酸酐酶放大 20 倍（pH 7.4，\(K_c = 20\)）——这是 CO2 高效从组织运到肺的化学基础。
  2. O2 运输被 Hb 放大 40 倍——20 volumes percent vs 0.5 volumes percent（无 Hb）。
  3. CO 中毒阈值 0.1% 大气浓度（\(P_{CO} = 0.7\) mmHg）——CO 与 Hb 亲和力高 200 倍。
  4. 高压氧治疗使 CO 清除加快 14 倍——\(w^*\) 从 3.5 升至 50.7。
  5. O2 沿毛细血管的摄取不是指数的，而是先线性快速上升（Hb 在陡坡段 \(f'(W)\) 大）后饱和——这是 Fig. 14.5 曲线的解释。
  6. 通气/灌注比 \(\dot{V}/Q \approx 1\)（正常）；过度通气 → 呼吸性碱中毒；通气不足 → 呼吸性酸中毒。
  7. 呼吸交换比 RQ：碳水 1.0、脂肪 0.7、混合饮食 0.825。
  8. Cheyne-Stokes 呼吸 = Hopf 分岔：延迟 \(\sigma > \sigma_c\)（典型心衰时 \(\sigma\) 增大）或增益 \(F'\) 增大（脑损伤）→ 40-60 s 周期的通气振荡。
  9. Lange-Hecht 1962：9 例 Cheyne-Stokes 患者中 7 例 2 年内死亡——Cheyne-Stokes 呼吸是心衰预后不良的独立预测因子。
  10. Grodins 模型的 \(G_C\) 阈值：\(G_C = 1.8\) 稳定，\(G_C = 9\) 振荡——脑 CO2 敏感性增大到一定阈值才会产生 Cheyne-Stokes。

挑战和开放性问题

  1. Cheyne-Stokes 呼吸的精确延迟常数 \(\tau\) 难以测定——脑血流时间受心输出量、血管阻力、肺-脑距离多因素影响，临床不易实时测量。
  2. pre-BötC 是否是真正的 pacemaker——Feldman & Del Negro 2006 主张节律来自网络连接而非细胞内禀振荡，但 14.4 节说"there seems little doubt that such control is exerted by multiple feedbacks"。细胞内禀振荡 vs 网络振荡仍未定论。
  3. Grodins 模型的非线性——(14.78) 中 \(\max(P_B - I_C, 0)\) 是非光滑的，全局稳定性分析困难。
  4. O2 与 Hb 结合的精确动力学——14.1.4 用 Hill 方程 \(f(W) = W^4/(K_4 + W^4)\)，但 Ch 13 显示真实模型是 Adair 中间模型（4 个 \(K_j\) 不相等）。两者预测的 O2 释放曲线在生理范围差异约 5-10%。
  5. 胎儿-新生儿呼吸转换——未在本章展开。第一口呼吸的触发（温度、化学、温度感受器）、新生儿呼吸窘迫综合征（NRDS，表面活性剂缺乏）的模型缺失。
  6. 呼吸暂停综合征（OSA）——上呼吸道阻塞的力学模型未涉及。OSA 与 Cheyne-Stokes 呼吸数学上同构（都是反馈失稳），但物理机制不同。
  7. 高原适应的红细胞增多——14.1 节说"the rate of red blood cell production can be increased as much as tenfold"，但高原长期适应的红细胞过度增多（慢性高原病）的反馈机制未量化。
  8. 运动时通气增加的精确机制——\(\dot{V}\) 在运动时可增加 5 倍，但 \(\dot{V}\) vs \(\dot{V}_{CO_2}\) 的精确斜率（运动刺激 vs 化学刺激的相对贡献）仍存争议。

个人反思与批判性分析

本章的最大优点在于 14.3 节的 Mackey-Glass 模型——它用一个 1 个状态变量 + 1 个时间延迟的简单 DDE，定量解释了 Cheyne-Stokes 呼吸这一临床现象。这是一个 "模型最简性与临床相关性" 的完美平衡：模型有 4 个参数（\(\lambda, \alpha, V_m, \tau\)），预测 40-60 s 周期，与临床实测吻合，并明确区分心衰（延迟增大）和脑损伤（增益增大）两种机制。

本章的主要缺陷在于 14.1 节给出的 CO2/HCO3⁻ 模型过于简化——假设 [H+] 常数、pH 固定为 7.4，但 14.1.3 节末明确说"\(\log_{10}(k_1/k_{-1}) = -6.1\) at pH 7.4" 而不同组织 [H+] 差异很大。真实 pH 变化（7.4 ↔ 7.35）对 \(K_c\) 的影响未在模型中显式追踪。

14.2 节的 O2-CO2 图（Fig. 14.9）是一个临床医生熟悉的工具，本章把它与 \(\dot{V}/Q\) 和 R 的关系清晰展示。这是基础呼吸生理的 "图论"——所有可能状态被压缩到一条曲线上。但本章没有讨论：(a) V/Q 不均匀性（不同肺段的 \(\dot{V}/Q\) 可差 5 倍）；(b) 分流（shunt）（\(\dot{V}/Q = 0\) 但仍有血流）；(c) 弥散障碍（如肺纤维化）。这些扩展对临床肺科医生重要，但超出了"模型最简"的范围。

对我自己研究血管生物力学的启发：(1) CO2 运输的化学放大（碳酸酐酶 + Hb）是一个生物化学缓冲系统的标准案例——血管壁细胞内的 pH 缓冲（HCO3⁻/H+）、细胞外基质的 Ca²⁺ 缓冲（钙结合蛋白）、血液中的 H+ 缓冲（Hb）都遵循同样的"化学放大"逻辑；(2) Mackey-Glass 方程可以推广到血管动力学——血管平滑肌的钙振荡、血管壁的自身调节反馈、心率变异性中的低频振荡都可以用类似 DDE 框架建模；(3) CO2-Hb 相互增强（Bohr 效应）的热力学对称性与血流-代谢耦合的"神经-血管耦合"（neurovascular coupling）在结构上类似——局部 O2 浓度下降触发血管舒张，血管舒张增加 O2 供应，二者满足某种"化学平衡"。

值得指出，本章没有讨论 WKB 近似下的扩散边界层——当 \(D_m L/v\) 很大时，O2 在毛细血管入口附近的薄层内完成大部分摄取，这一"边界层"结构对运动生理（高血流量下 O2 提取率下降）和微循环数学都重要。Tawhai et al. 2004-2006 的全肺模型包含这些细节，但本章没有引用其结果。

重要参考文献

[X1] Glass, L., & Mackey, M. C. (1988). From Clocks to Chaos: The Rhythms of Life. Princeton University Press. — 14.3 节 Mackey-Glass 方程 (14.57) 的来源；Cheyne-Stokes 呼吸的延迟微分方程原型。 [X2] Grodins, F. S., Buell, J., & Bart, A. J. (1967). Mathematical analysis and digital simulation of the respiratory control system. Journal of Applied Physiology, 22, 260-276. — 14.3.1 节多室延迟 ODE 呼吸控制模型的原型。 [X3] Khoo, M. C. K., Kronauer, R. E., Strohl, K. P., & Slutsky, A. S. (1982). Factors inducing periodic breathing in humans: A general model. Journal of Applied Physiology, 53, 644-659. — 14.3.1 节 Grodins 模型的简化版（用于重现 Cheyne-Stokes）。 [X4] Smith, J. C., Ellenberger, H. H., Ballanyi, K., Richter, D. W., & Feldman, J. L. (1991). Pre-Bötzinger complex: A brainstem region that may generate respiratory rhythm in mammals. Science, 254, 726-729. — 14.4 节 pre-Bötzinger 复合体作为节律生成关键区的发现。 [X5] Butera, R. J., Rinzel, J., & Smith, J. C. (1999a). Models of respiratory rhythm generation in the pre-Bötzinger complex. I. Bursting pacemaker neurons. Journal of Neurophysiology, 82, 382-397. — 14.4 节 preBötC 神经元的 HH 模型。 [X6] Butera, R. J., Rinzel, J., & Smith, J. C. (1999b). Models of respiratory rhythm generation in the pre-Bötzinger complex. II. Populations of coupled pacemaker neurons. Journal of Neurophysiology, 82, 398-415. — 14.4 节耦合 preBötC 神经元的网络模型。 [X7] Feldman, J. L., & Del Negro, C. A. (2006). Looking for inspiration: New perspectives on respiratory rhythm. Nature Reviews Neuroscience, 7, 232-242. — 14.4 节双振荡器假说（preBötC + RTN）。 [X8] von Euler, C. (1980). Central pattern generation during breathing. In Neurobiology of the Locomotor System. Elsevier. — 14.4.1 节互抑制网络模型的原型。 [X9] Wyman, R. J. (1977). Neural generation of the breathing rhythm. Annual Review of Physiology, 39, 417-448. — 14.4.1 节互抑制模型的独立提出。 [X10] Lange, R. L., & Hecht, H. H. (1962). The mechanism of Cheyne-Stokes respiration. Journal of Clinical Investigation, 41, 42-52. — 14.3 节末尾引用的临床数据；9 例 Cheyne-Stokes 患者 7 例 2 年内死亡。 [X11] Ben-Tal, A. (2006). Simplified models for gas exchange in the human lungs. Journal of Theoretical Biology, 238, 474-495. — 14.1 节末尾综述：从最简到复杂的肺气体交换模型。 [X12] Whiteley, J. P., Gavaghan, D. J., & Hahn, C. E. W. (2001, 2002, 2003a,b). A series of papers on gas exchange modeling. Mathematical Medicine and Biology. — 14.1 节末尾引用的红细胞形态对气体交换影响的详细模型。 [X13] Tawhai, M. H., Hunter, P., Tschirren, J., Reinhardt, J., McLennan, G., & Hoffman, E. A. (2004). CT-based geometry analysis and finite element models of the human and ovine bronchial tree. Journal of Applied Physiology, 97, 2310-2321. — 14.1 节末尾引用的全肺分支结构计算模型。 [X14] Guyton, A. C., & Hall, J. E. (1996). Textbook of Medical Physiology (9th ed.). W. B. Saunders. — 14.1.1 Fig. 14.1, 14.4, 14.5, 14.11 临床数据出处。 [X15] Roughton, F. J. W. (various years). Series of papers on Hb kinetics and the Bohr effect. Journal of Applied Physiology and Respiration Physiology. — 14.1.3 节碳酸酐酶反应速率常数和 \(K_c\) 值的测定。