三、二室模型一级消除动力学
多室模型和前面讨论的单室模型不同，此时存在着由包括血液在内的中央室向周边室分布达到平衡的过程，影响血药浓度的因素更为复杂，下面以单剂静脉注射为例，简介有关药动学知识。
（一）模式图和药-时关系
静脉注射时，不存在吸收过程，中央室中的药量或血药浓度受中央室与周边室间双向分布，以及自中央室消除的影响。其模式如图9-9所示。
x[XB]c[/XB]中央室药量
x[XB]p[/XB]周边室药量
v[XB]c[/XB]中央室分布容积
v[XB]p[/XB]周边室分布容积
k[XB]10[/XB]中央室向周边室转运速率常数
k[XB]21[/XB]周边室向中央室转运速率常数
图9-9 二室分布静脉注射模式图
中央室药量变化的速率微分方程为：
dk[XB]c[/XB]/dt=k[XB]21[/XB]·x[XB]p[/XB]-k[XB]21[/XB]·x[XB]c[/XB]-k[XB]10[/XB]·x[XB]c[/XB]
对上式积分并引入中央室分布容积Vc，可得中央室（血）药物浓度随时间变化的基本表达式：
C＝A·e＋B·e[SB]-βt[/SB]式(25)
式中α为分布速率常数，β为消除速率常数，A、B为经验常数。四者都是由模式参数k[XB]10[/XB]、k[XB]12[/XB]、k[XB]21[/XB]组成的混杂参数（hybridparameters）。存在：
α·β＝k[XB]21[/XB]·k[XB]10[/XB]α＋β＝k[XB]10[/XB]＋k[XB]12[/XB]＋k[XB]21[/XB]
（二）药动学参数计算
二室模型静脉注射药动学参数的求算，类似于单室模型血管外用药，仍采用残数法。即因为α＞β,当t充分大时，A·e[SB]-αt[/SB]→0，则式(25)变为
C＝Be[SB]-βt[/SB]
此即消除相药-时关系表达式。按前述方法可求得B、β和消除半寿期t[XB]1/2β[/XB]。将式(25)移项可得：
C－B·e[SB]-αt[/SB]＝A·e[SB]-αt[/SB]
令C[XB]r[/XB]＝C－B·e[SB]-βt[/SB]，C[XB]r[/XB]即为消除相以前某时点实测血药浓度减去消除相该点外推浓度的残数。代入上式可得分布相药-时关系表达式
C[XB]r[/XB]＝A·e[SB]-αt[/SB]
同样可求得A、α和分布半寿期t[XB]1/2α[/XB]（图9-10）。
图9-10 二室模型静脉注射血药浓度-时间关系示意图
再根据下列公式，可求得各有关药动学参数。
中央室表观分布容积：Vc＝X[XB]0[/XB]/（A＋B）
周边室向中央室的转运速率常数：
自中央室消除的速率常数：
k[XB]10[/XB]＝α·β/k[XB]21[/XB]
中央室向周边室转运的速率常数：
k[XB]12[/XB]＝（α＋β）－（k[XB]10[/XB]＋k[XB]21[/XB]）
曲线下面积：AUC[XB]0→∞[/XB]＝A/α+B/β（积分外推法）
周边室表观分布容积：V[XB]p[/XB]＝V[XB]c[/XB]·k[XB]12[/XB]/k[XB]21[/XB]
总表观分布容积：
稳态表观分布容积：V[XB]ss[/XB]＝V[XB]c[/XB]＋V[XB]p[/XB]
消除率：Cl＝k[XB]10[/XB]·V[XB]c[/XB]＝β·V
血管外用药时，由于存在吸收因素，求算药动学参数的基本方法仍为残数法，但此时需进行两次残数处理。有关二室模型血管外用药及二室以上多室模型的药动学，在TDM的实际工作中很少应用，可参阅有关药动学专著。