第十七章：肾脏生理学（Renal Physiology）
书籍信息：Keener & Sneyd, Mathematical Physiology (II: Systems Physiology), Springer 2009, Chapter 17
学习笔记撰写日期：2026年5月11日
一、概述
肾脏是维持生命体稳态的核心器官，执行两大主要生理功能：第一，排泄体内代谢终产物；第二，调节体液中大多数化学成分的浓度。本章的数学生理学分析聚焦于尿液形成的过程及废物如何从血液中被清除，同时简要讨论体液成分的调控机制。
肾脏的基本功能单位是肾单位（nephron），每个肾脏约含一百万个肾单位。每个肾单位都能独立形成尿液。血液通过入球小动脉（afferent arteriole）进入肾单位，该小动脉位于肾皮质；肾单位的肾小管及周围的管周毛细血管延伸到肾髓质深处。肾单位的主要功能结构包括：肾小球（glomerulus），血液在此被过滤；球旁器（juxtaglomerular apparatus），控制肾小球血流；以及长肾小管，过滤液在此转化为尿液。
在人类正常生理条件下，肾小球滤过率（glomerular filtration rate, GFR）约为125毫升/分钟。滤过分数（filtration fraction）是指成为肾小球滤过液的血浆量占总肾血浆流量的比例，约为20%。超过99%的滤过液在肾小管中被重吸收，仅有小部分成为尿液排出体外。本章的数学模型主要关注三个核心问题：肾小球的过滤机制与自身调节、肾小管液流动的振荡行为、以及尿液浓缩的逆流倍增机制。
二、肾小球：过滤与自身调节
2.1 肾小球过滤模型
肾小球是尿液形成的第一阶段，通过血液血浆产生滤过液。肾小球是由多达50条平行分支且相互吻合的毛细血管组成的网络，被上皮细胞覆盖并包裹在鲍曼囊（Bowman's capsule）中。血液通过入球小动脉进入肾小球，通过出球小动脉（efferent arteriole）离开。肾小球内的血压驱动液体过滤到鲍曼囊中，同时将所有小分子量溶解物质携带出来。肾小球膜对所有血浆蛋白几乎完全不通透（白蛋白的分子量为69,000道尔顿），因此肾小球滤过液与血浆几乎完全相同，只是蛋白质含量极低。
肾小球过滤的数学模型将肾小球毛细血管假设为一维管道，流量为 \(q_1\)，周围鲍曼囊也等效为一维管道，流量为 \(q_2\)。由于跨越肾小球毛细血管壁的过滤流与毛细血管壁两侧的压力差成正比，在稳态下有：
其中 \(P_1\) 和 \(P_2\) 分别为管道1（毛细血管）和管道2（鲍曼囊）的静水压，\(\pi_c\) 是血浆中悬浮蛋白的渗透压，\(K_f\) 为毛细血管过滤系数。渗透压由下式给出：
其中 \(c\) 为浓度（单位：摩尔/升），是 \(x\) 的函数，因为悬浮液在通过肾小球的过程中逐渐浓缩。由于大分子蛋白绕过过滤器，根据物质守恒有：
其中 \(c_i\) 为输入浓度，\(Q_i\) 为输入流量。由此可得：
其中 \(\pi_i = RTc_i\) 为输入渗透压。由于肾小球内的静水压降相对于入球和出球小动脉的压力降很小，可将 \(P_1\) 和 \(P_2\) 视为常数。
方程（17.1）结合（17.4）给出关于 \(q_1\) 的一阶常微分方程，易于求解。设 \(q_1(L) = Q_e\)，得到：
其中 \(Q_e\) 为出球小动脉流量，\(L\) 为过滤器长度，\(\alpha = \pi_i/(P_1 - P_2)\)。
此外，输入和输出小动脉控制压力和流速，通过以下方程描述：
肾小球到近端小管的流量由下式控制：
其中 \(P_a\)、\(P_e\) 和 \(P_d\) 分别为入球小动脉、出球小动脉和下降小管的压力，\(R_a\)、\(R_e\) 和 \(R_d\) 分别为入球和出球小动脉以及近端小管的阻力。典型值为：\(P_1 = 60\) mmHg，\(P_2 = 18\) mmHg，\(P_a = 100\) mmHg，\(P_e = 18\) mmHg，\(P_d = 14-18\) mmHg，\(\pi_i = 25\) mmHg，\(Q_i = 650\) ml/min，\(Q_d = Q_i - Q_e = 125\) ml/min。
2.2 自身调节与管球反馈振荡
肾小球滤过率的自身调节至关重要。若滤过液流速过慢，重吸收将过度，导致肾脏无法清除必要的代谢废物；若流速过快，肾小管无法重吸收那些需要保留的物质，造成有价值物质的尿中丢失。
肾小球通过检测升支粗段顶部 NaCl（很可能是氯离子）浓度来控制过滤速率。亨利袢（loop of Henle）降支深入肾髓质后，在肾小球附近返回，分别靠近入球和出球小动脉。球旁复合体由远曲小管中的致密斑（macula densa）细胞和入球及出球小动脉壁中的球旁细胞组成。
升支粗段的主要功能是将 Na\(^+\) 从肾小管腔主动泵入间质；该段对水几乎完全不通透。因此，流速越慢，从升支粗段主动转运出的 Na\(^+\) 就越多，Cl\(^-\) 则被动跟随。流速过低导致 Na\(^+\) 和 Cl\(^-\) 重吸收过度，使这些离子浓度在亨利袢末端下降过多；流速过高则使致密斑处 Na\(^+\) 和 Cl\(^-\) 浓度较高。
致密斑细胞对 NaCl 浓度下降的响应（机制尚未完全阐明）是释放血管舒张物质，降低入球小动脉阻力。同时球旁细胞释放肾素（renin），该酶促进血管紧张素II的形成，后者收缩出球小动脉。两者的共同作用是增加肾小球滤过液流量。相反，致密斑处 NaCl 浓度升高则降低入球小动脉阻力并增加出球小动脉阻力，从而减少过滤和肾小管液流速。
这种控制系统可以导致亨利袢液流速率的自发振荡。Layton 等人（1991）提出的管球反馈振荡（tubuloglomerular oscillations）简化模型重点研究升支粗段氯离子浓度的作用。该模型将升支粗段模拟为一维管道，Cl\(^-\) 通过液体流动沿管道运输，同时被主动移除（实际上是 Na\(^+\) 被主动移除，但 Cl\(^-\) 被动跟随）。
设 \(y = 0\) 和 \(y = L\) 分别表示长度为 \(L\) 的升支粗段的起点和终点，\(C(y, \tau)\) 表示 \(y\) 处 \(\tau\) 时刻的 Cl\(^-\) 浓度。Cl\(^-\) 的守恒方程为：
其中 \(R(C)\) 表示通过 Na\(^+\) 泵驱动从管道移除 Cl\(^-\) 的速率，\(\phi\) 表示液体沿管道流动的速度。简化模型假设 Cl\(^-\) 通过一级过程被移除：
当致密斑处 NaCl 浓度改变时，入球小动脉阻力的改变存在时间延迟。因此 \(\phi\) 是致密斑处过去某时刻 Cl\(^-\) 浓度的函数，即 \(C(L, \tau - \bar{\tau})\) 的函数。实验研究确定 \(\phi\) 的函数形式为：
注意 \(\phi\) 是 \(C\) 的递减函数，\(\bar{C}\) 是 \(\phi\) 具有拐点的 Cl\(^-\) 浓度。
设 \(y = Lx\)，\(\tau = F_{op}^{-1} t\)，\(\bar{\tau} = F_{op}^{-1} \bar{t}\)，\(C = C_0 c\)，\(K_1 = \tilde{F}/F_{op}\)，\(K_2 = \bar{C}\)，守恒律化为带时滞的偏微分方程：
其中 \(F(c) = 1 + K_1 \tanh(K_2(\bar{c} - c))\)，\(\mu = r_c/F_{op}\)，边界条件为 \(c(0, t) = 1\)。
稳态解 \(s(x)\) 满足：
由此得到：
由于 \(F\) 是递减函数，存在唯一的 \(k\) 值使方程成立，因此稳态解唯一。
稳定性分析通过线性化进行：\(c(x, t) = s(x) + u(x, t)\)，代入方程并保留一阶项，得到关于 \(u\) 的方程。设 \(u(x, t) = f(x)e^{\lambda t}\)，可推导出特征方程：
其中 \(\gamma = -kF'(e^{-k})e^{-k} > 0\)，与反馈调节的敏感性有关。
Hopf 分岔发生在 \(\lambda\) 实部为零且虚部非零时，由此可推导出参数空间中的分岔曲线：
该分析表明，时滞 \(\bar{t}\) 和敏感性反馈参数 \(\gamma\) 是决定振荡发生的关键参数。振荡发生在两个参数中任一足够大时（但两者均非零）。当 \(\mu\) 偏离最优值 \(\mu = -\ln\bar{c}\) 时，\(\gamma\) 减小，振荡发生的可能性降低。
三、尿液浓缩：亨利袢
3.1 亨利袢的结构与功能
尿液浓缩或稀释机制的核心挑战在于理解浓缩和稀释机制如何协同工作以决定尿液成分，并定量解释特定物种的浓缩能力。例如，人类最大尿液浓缩能力为1200 mOsm/升，而某些沙漠动物（如澳大利亚跳跃鼠）可将尿液浓缩至高达10,000 mOsm/升。目前尚不完全清楚如此高浓度尿液是如何产生的。
一个70公斤的成年人每天必须排出约600 mOsm 的溶质（代谢废物和摄入的离子）。运输这些溶质所需的最小尿量（强制性尿量）为：
这解释了为什么饮用海水会导致严重脱水。海洋中盐的浓度平均为3%氯化钠，渗透压在2000至2400 mOsm/升之间。饮用1升浓度为2400 mOsm/升的海水将摄入2400 mOsm 的溶质，这些溶质必须被排出。若最大尿浓缩能力为1200 mOsm/升，则需要2升尿液来清除这些摄入的溶质，产生1升的净失水。这解释了为什么船难受害者饮用海水会迅速脱水，而澳大利亚跳跃鼠却可以毫无顾虑地饮用大量海水。
尿液浓缩或稀释主要在亨利袢中完成。离开鲍曼囊后，肾小球滤过液流入一条由五段组成的肾小管：近端小管、亨利袢降支、亨利袢升支、远曲小管和集合管。这些肾小管被管周毛细血管包围，后者重吸收从肾小管中提取的液体。
近端小管的主要功能是提取大部分水和溶解的化学物质（电解质、葡萄糖、各种氨基酸等），使其重吸收回血液，同时浓缩代谢废物。近端小管上皮细胞含有大量线粒体以支持快速的主动转运过程。实际上，约65%的肾小球滤过液在到达亨利袢降支之前就被重吸收了。葡萄糖、蛋白质、氨基酸、乙酰乙酸离子和维生素几乎完全通过沿近端小管上皮细胞的主动协同转运过程被重吸收。
亨利袢降支衬有薄上皮细胞，线粒体很少，表明代谢活性低；该段对水高度通透，对 Na\(^+\)、尿素和大多数离子为中等通透性。亨利袢升支始于薄壁，但约一半高度处明显增厚。与降支相反，升支对水和尿素高度不通透。升支粗段细胞类似于近端小管细胞，适合将 Na\(^+\) 和 Cl\(^-\) 从肾小管腔主动转运到间质液中。
3.2 逆流倍增机制
亨利袈尿液浓缩的定性的工作原理如下：沿升支，Na\(^+\) 被吸收到间质中（薄升支为被动吸收，升支粗段为主动吸收）。这在间质中创造了高 Na\(^+\) 浓度，随后将水从降支抽出，并使 Na\(^+\) 重新进入降支。因此，进入降支的液体逐渐浓缩，直到在亨利袢转折点处，液体渗透压约为1200 mOsm/升（进入液体约为300 mOsm/升）。显然，由于进入升支的液体已经如此浓缩，从升支提取 Na\(^+\) 得以增强，从而进一步增强从降支提取水，如此循环。这种正反馈过程是逆流倍增机制的核心。
当无抗利尿激素（ADH）存在时，亨利袢形成的稀释尿基本不变地通过集合管，产生大量稀释尿。当存在大量 ADH 时，集合管对水高度通透，因此当滤过液到达亨利袢转折点水平时，基本上与间质处于同一浓度（约1200 mOsm/升），从而产生少量浓缩尿。
四、逆流倍增机制的数学模型
4.1 逆流交换器的基本方程
逆流倍增机制通过跨膜扩散进行溶质交换。由于交换速率受膜两侧浓度差的影响，维持较大浓度差可以提高交换效率。假设两股含有溶质的气体或液体沿长度 \(L\) 的平行管道流动，中间由通透性膜隔开。最简单的一维模型假设溶质转运是浓度差的线性函数，则两个管道中的浓度由以下方程描述：
其中 \(d\) 为交换系数，\(q_1\) 和 \(q_2\) 为两个管道中的体积流率。
在稳态下，假设输入浓度为 \(C_{01}\) 和 \(C_{02}\)，将两个方程相加并积分得到：
消去 \(C_2\) 得到关于 \(C_1\) 的微分方程：
其解为：
其中 \(\kappa = \frac{k}{q_1 + q_2}\)，\(\lambda = \frac{d(q_1 + q_2)}{q_1 q_2}\)。
4.2 同向流与逆向流的比较
需要考虑两种情况：\(q_1\) 和 \(q_2\) 同号（同向流）以及 \(q_1\) 和 \(q_2\) 异号（逆向流）。
当 \(q_1\) 和 \(q_2\) 同向时（均为正），输入端在 \(x = 0\)，有 \(C_1(0) = C_{01}\)，\(C_2(0) = C_{02}\)，输出浓度比为：
其中 \(\gamma = C_{02}/C_{01}\)，\(\rho = q_2/q_1\)，\(\lambda = \frac{d(1 + \rho)}{q_1}\)。
当 \(q_1\) 和 \(q_2\) 异号时（设 \(q_1 > 0\)，\(q_2 < 0\)），管道1的输入在 \(x = 0\)，管道2的输入在 \(x = L\)。输出浓度比为：
其中 \(\gamma = C_{02}/C_{01}\)，\(\rho = -q_2/q_1 > 0\)，\(\lambda = \frac{d(1 - \rho)}{q_1}\)（\(\rho \neq 1\) 时）。当 \(\rho = 1\) 时，有特殊解：
当 \(\gamma < 1\)（即物质从管道1向管道2转移）时，方程（17.36）中 \(C_1(L)/C_{01}\) 的值始终大于方程（17.37）中相应的值。这表明在相同参数条件下，逆流交换的溶质总转运效率始终高于同向流。
五、肾单位中的逆流倍增
5.1 肾单位逆流机制的特殊性
在肾单位中，逆流机制的工作方式与上述简化模型略有不同。两个平行管道——亨利袢的降支和升支——在底端相互连接。因此，降支流出管道的流量和浓度必须与升支流入管道的流量和浓度相匹配。
尿液浓缩机制的数学模型已存在一段时间，所有模型都利用相同的基本物理原理：通过主动转运过程建立化学梯度、通过扩散进行离子运动、以及通过渗透作用进行水的转运。肾单位的独特之处在于：尽管所有脊椎动物都能产生稀释尿，但只有鸟类和哺乳动物能产生高渗尿，而这些动物的肾脏才含有亨利袢。
逆流倍增机制的核心在于：沿升支粗段，Na\(^+\) 和 Cl\(^-\) 被主动泵入肾间质，但该段对水不通透。因此，当滤过液沿升支向上流动时，Na\(^+\) 和 Cl\(^-\) 浓度逐渐降低（因为被不断泵出），而水含量保持恒定。间质中的高 Na\(^+\) 浓度通过渗透作用将水从降支抽出进入间质，而降支对水高度通透。这样形成了一个正反馈循环：升支去除 Na\(^+\) 越多，间质渗透压越高；间质渗透压越高，从降支抽出的水越多；降支液体越浓缩，到达升支的液体 Na\(^+\) 浓度越高，从而支持更多的主动 Na\(^+\) 泵出。这个循环在亨利袢转折点处产生约1200 mOsm/升的高渗环境。
六、肾小管转运的模型
6.1 近端小管的主动与协同转运
近端小管上皮细胞通过三种主要机制将物质转运至血液：第一，Na\(^+\)–K\(^+\) ATP 酶将 Na\(^+\) 从上皮细胞内主动转运至间质，建立 Na\(^+\) 的浓度梯度。第二，次级主动转运蛋白利用 Na\(^+\) 梯度（由 ATP 酶建立）将葡萄糖、氨基酸、磷酸盐、Ca\(^{2+}\) 和 Mg\(^{2+}\) 等从肾小管腔转运到上皮细胞内。第三，水通过细胞膜上的渗透压介导的通道（由水孔蛋白）进行转运。
Na\(^+\) 在上皮细胞膜上与 H\(^+\) 交换体（将 H\(^+\) 从细胞转运到肾小管腔以换取 Na\(^+\)）也是重要的转运机制之一。近端小管中约65%的滤过液在此被重吸收，同时葡萄糖、氨基酸等几乎100%被协同转运重吸收。
6.2 集合管的激素调节
远曲小管和集合管细胞对多种激素敏感，这些激素调节其功能以及尿液最终化学成分。其中最重要的激素是醛固酮（aldosterone）和抗利尿激素（ADH）。醛固酮决定 Na\(^+\) 从肾小管腔转运出去的速率，而 ADH 则决定集合管对水的通透性，从而决定尿液的最终浓度。
当无 ADH 存在时，集合管对水完全不通透，稀释尿的形成不受阻碍；当 ADH 大量存在时，集合管对水高度通透，水分被大量重吸收，使尿液浓缩至高渗状态（最高可达1200 mOsm/升）。
七、公式汇总表
| 编号 | 公式名称/描述 | 公式 |
|---|---|---|
| (17.1) | 肾小球过滤方程 | \(\dfrac{dq_1}{dx} = K_f(P_2 - P_1 + \pi_c)\) |
| (17.2) | 胶体渗透压 | \(\pi_c = RTc\) |
| (17.3) | 蛋白守恒 | \(c_i Q_i = c q_1\) |
| (17.4) | 渗透压与流量关系 | \(\pi_c = \pi_i \dfrac{Q_i}{q_1}\) |
| (17.5) | 肾小球过滤稳态方程 | \(\dfrac{Q_e}{\alpha} + \alpha \ln\dfrac{Q_i - \alpha}{1 - \alpha} = 1 - K_f L \dfrac{\pi_i}{\alpha Q_i}\) |
| (17.6) | 入球小动脉压力-流量 | \(P_a - P_1 = R_a Q_i\) |
| (17.7) | 出球小动脉压力-流量 | \(P_1 - P_e = R_e Q_e\) |
| (17.8) | 肾小球到近端小管 | \(P_2 - P_d = R_d(Q_i - Q_e)\) |
| (17.9) | Cl\(^-\) 守恒方程 | \(\dfrac{\partial C}{\partial \tau} + \phi \dfrac{\partial C}{\partial y} = -R(C)\) |
| (17.10) | 一级 Cl\(^-\) 移除 | \(R(C) = -r_c C\) |
| (17.11) | 流量对浓度的依赖 | \(\phi(C) = F_{op} + \tilde{F} \tanh(\alpha(\bar{C} - C))\) |
| (17.12) | 无量纲化后的带时滞 PDE | \(\dfrac{\partial c}{\partial t} + F(c(1, t - \bar{t})) \dfrac{\partial c}{\partial x} = -\mu c\) |
| (17.14) | 稳态方程 | \(F(s(1)) \dfrac{ds}{dx} = -\mu s\) |
| (17.15) | 稳态解 | \(s(x) = s(0) e^{-kx},\quad k = \dfrac{\mu}{F(s(1))}\) |
| (17.17) | 敏感性定义 | \(\sigma = \dfrac{s(0)}{s(1)} \dfrac{ds(1)}{ds(0)}\) |
| (17.19) | 敏感性表达式 | \(\sigma = \dfrac{1}{1 + \mu k \gamma}\) |
| (17.24) | 特征方程（Hopf 分岔） | \(\lambda = \gamma e^{-\lambda \bar{t}}\left(e^{-\frac{\mu k}{\lambda}} - 1\right)\) |
| (17.28) | Hopf 分岔曲线 | \(\gamma = (-1)^{n+1} \dfrac{\omega}{2\sin\frac{k\omega}{2\mu}},\ \omega = \dfrac{n\pi/k}{\bar{t} + k/(2\mu)}\) |
| (17.29) | 强制性尿量 | \(\text{obligatory volume} = \dfrac{\text{total solute/day}}{\text{maximal urine concentration}}\) |
| (17.30) | 人类强制性尿量计算 | \(\dfrac{600 \text{ mOsm/day}}{1200 \text{ mOsm/L}} = 0.5 \text{ L/day}\) |
| (17.31) | 逆流交换器方程1 | \(\dfrac{\partial C_1}{\partial t} + q_1 \dfrac{\partial C_1}{\partial x} = d(C_2 - C_1)\) |
| (17.32) | 逆流交换器方程2 | \(\dfrac{\partial C_2}{\partial t} + q_2 \dfrac{\partial C_2}{\partial x} = d(C_1 - C_2)\) |
| (17.33) | 逆流稳态不变量 | \(q_1 C_1 + q_2 C_2 = k\)（常数） |
| (17.35) | 交换器稳态解 | \(C_1(x) = \kappa + (C_1(0) - \kappa)e^{-\lambda x}\) |
| (17.36) | 同向流输出浓度比 | \(\dfrac{C_1(L)}{C_{01}} = \dfrac{1 + \gamma\rho}{1 + \rho} + \dfrac{1 - \gamma - \lambda L}{1 + \rho} e^{-\lambda L}\) |
| (17.37) | 逆向流输出浓度比 | \(\dfrac{C_1(L)}{C_{01}} = \dfrac{-\gamma\rho + (1 - \rho + \gamma\rho)e^{-\lambda L}}{e^{-\lambda L} - \rho}\) |
| (17.38) | 逆向流 \(\rho=1\) 特殊解 | \(\dfrac{C_1(L)}{C_{01}} = \dfrac{q_1 + \gamma dL}{q_1 + dL}\) |
八、总结与延伸阅读
本章系统地介绍了肾脏生理学的数学生理学建模方法，涵盖了从肾小球过滤到尿液浓缩的完整生理过程。主要内容包括以下几个方面：
第一，肾小球过滤的物理模型。该模型将肾小球视为一条一维管道，过滤过程由静水压和胶体渗透压共同驱动。通过引入毛细血管过滤系数 \(K_f\) 和血浆蛋白渗透压 \(\pi_c = RTc\)，建立了过滤流量 \(q_1\) 沿肾小球变化的微分方程，并求得稳态解析解。该模型揭示了肾小球滤过率（GFR）随动脉压变化的规律。
第二，管球反馈振荡的动力学分析。肾小球通过致密斑检测亨利袢升支粗段 Cl\(^-\) 浓度，经时滞反馈调节入球和出球小动脉阻力。建立带时滞的偏微分方程模型后，通过线性稳定性分析发现，当时滞 \(\bar{t}\) 和敏感性参数 \(\gamma\) 足够大时，系统通过 Hopf 分岔产生自振荡。这一定性结果与实验观察到的近端小管压力自发性振荡相吻合。
第三，逆流倍增机制的数学描述。将亨利袢的降支和升支模拟为平行管道中的逆流交换系统，通过求解带交换项的对流-扩散方程，证明逆流（反向流动）比同向流具有更高的溶质交换效率。当两管道流量相等（\(\rho=1\)）时，逆流输出浓度与输入浓度之比为 \((q_1 + \gamma dL)/(q_1 + dL)\)，显著优于同向流的结果。
第四，尿液浓缩的生理学机制。通过升支粗段主动泵出 Na\(^+\)（对水不通透）和降支对水的通透性的配合，建立起肾髓质间质的高渗梯度。这一正反馈机制使得亨利袢转折点处的渗透压可达1200 mOsm/升，配合 ADH 对集合管通透性的调节，实现尿液浓度的精细控制。
关于延伸阅读，本章引用的经典文献包括 Leyssac 和 Baumbach（1983）首次描述管球反馈振荡的实验研究、Holstein-Rathlou 和 Marsh（1989, 1990）对管球振荡的详细理论和实验研究、以及 Layton 等人（1991, 1995, 1997, 2000）系列论文中的数学建模工作。Just（2006）的综述文章是了解该领域全貌的良好起点。在更复杂的生理建模层面，Pitman 等人（2004）的工作将自身调节和管球振荡整合到更完整的肾单位模型中。此外，Guyton 和 Hall（1996）的《医学生理学教科书》是理解肾脏生理学背景知识的经典参考书。