体内药物浓度的变化

1、体内药物浓度的变化摘要本文研究体内药物浓度变化的相关问题。在这个问题的研究中，我们采用了典型的房室模型。将人体看作是一个单房室系统，并 且针对几种不同的给药方式分别进行了讨论：针对第一种一次性给药的方式，将人体看作是只输岀不输入的房室模型，第一次给药之 后就不再有输入，从而构建常微分方程来求解：针对第二种等间隔服药的情况，药物以近似脉冲的方式进入人体，从而构建常微分方程 来求解：针对第三种恒速静脉点滴的情况，药物以近似恒速的方式进入人体，从而构建常微分方 程来求解。在各步骤的求解中，都是将人体看做房室模型，通过构造常微分方程来模拟求解关键词：药物浓度变化房室模型常微分方程一、问题重述医生给病人

2、开处方时必须注明两点：服药的剂量和服药的时间间隔.超剂量的药品 会对身体产生严重不良后果，甚至死亡，而剂量不足，则不能达到治病的目的.已知患者 服药后，随时间推移，药品在体内逐渐被吸收，发生生化反应，也就是体内药品的浓度逐 渐降低.药品浓度降低的速率与体内当时药品的浓度成正比.当服药疑为A、服药间隔为 T,试分析体内药的浓度随时间的变化规律.体内药物浓度的变化问题是分析病人进药之后，药物在体内被吸收，药物浓度在血 液之中含量上升，而作用在人身之后时间段内药物浓度在体内变化的情况。讨论问题：1、一次性注射的情况下，药物浓度在人体内如何变化？2、等间隔服药的情况下，药物浓度在人体内如何变化？3、恒

3、速静脉点滴的情况下，药物浓度在人体内如何变化？二、问题分析实际上，体内药物浓度的变化是一个很复杂的问题，涉及到很多因素。在本文讨论的方 法中，主要采用的是房室模型，就是忽略诸多的次要因素，只考虑一些最重要的影响作用最 大的因素，以此来简化问题。药物进入机体之后形成血药浓度（单位体积血液的药物量），在给药方案的设计中，应该 使药物浓度始终保持在一个合适的有效范囤内，既使药物浓度达到一个最有效的水平，又可 以防止浓度过高而对人产生副作用。我们借鉴药物动力学中的有关知识，建立房室模型，模 拟药物在体内吸收、分布和排除过程，我们建立的房室模型属于机体的一部分，药物在一个 房室内均匀分布（血药浓度为常数

4、），在房室间按一泄规律转移,利用这些规律来对问题进行模拟求解。确左了以上这些基本问题之后，就通过高等数学中的常微分方程来模拟药物浓度在体内 的变化。三、模型假设1、假设药物的分解排出速率与药物当前的浓度成正比；2、假设药物从服用到到达体液的过程中没有损失，即被完全吸收;3、假设药物在人体的体液中是均匀分布的：四、符号说明x(t)t时刻人体内的药物浓度值v(t)t时刻人体内药物消耗速率k药物消耗速率与人体内药物浓度的比例系数A每次的注射量v0恒速静脉点滴时体内药物浓度的增量T每次等间隔服药的时间五、模型建立假设注射药物后完全吸收，那么每次增量就是a,随着时间的推移，药物浓度会 下降，若一次性注射

5、则药物递降至无，若等间隔服药则药物浓度会在服药期上升，若恒 速注射则药物浓度可能最终恒速变化。一次注射药物后，体内浓度为a,此时为零时刻则：x(O)=a根据药物在人体内的消耗与人体内药物浓度成正比，列出方程dxv v=kx (t)dt得：=kx(t)dtkt求解上述方程：X (t)二C*幺将x(O)=a代入得C=a因此一次服药后浓度变化方程为x(r) = a * 0六. 模型分析分三种情况来对模型进行分析：1、一次性注射的情况：一次性注射情况即最简单的服药模型，如同上述解答，那么我们将代入实际数值来用 Mat lab模拟作图令a=0. b k=-0. b代入作图执行程序结果如下：上图即一次性服

6、药情况下，人体内药物浓度的变化曲线。从图中可以看出，随着时间 的推移，人体内的药物浓度最终会为0。2、等间隔服药的情况:考虑到等间隔服药的脉冲性得到： 其中T为每次服药之间的间隔时间 dx=一kxj H nTdtx(O) = a、x(nT) = a + 厂)在区间nT, (n+l)T上求解方程得到W) = xnT)ek,-nTt e 时,(” +1)门 在0WtVT内，方程的解为 x(/) = g",OStVT 在TWtV2T内，方程的解为 x(t)(a + aekT)e'k(,Ti在 2T t<Z3r内，方程的解为在心 7")=«+心 *由于方程的

7、解为V(Z)=(a + “L + “严)宀切,2 < (<3ra + 心 +。严 + d = a上上> 丄1 一戶丁 1_£” 由此看岀，在等间隔服药的情况下，药物的浓度在人体中呈上升趋势，且最后会稳定 在一定的水平,当T=8, k=0.1, a=0.1时,利用Mat lab编程得出数值计算结果如下： 运算得到体内药物的浓度应该维持在0. 0817-0. 1816的浓度范围。运算得到的示意图如下：c 0.180.160.140.120.10.080.063、静脉恒速注射的情况:静脉恒速注射时，药物会以一左速率V。持续注射到人体内，则人体内药物浓度的变化率会有变化，列

8、出方程：cLv=X) + V。at解上述方程得：X%+丸)* 才_土kk利用Mat lab编程模拟得岀如下图形：从模拟出的图形中可以看出，随着时间的推移，药物浓度会逐渐上升，但是最终会达到 一个恒泄的浓度，针对这个模型，得岀的最终药物浓度维持在0.1而不再变化，这个浓度显 然是和滴定的速率是有关的。七、模型检验1、一次性注射的情况，查阅文献资料的到如下的数据：此为单次静脉给药PS916血液浓度的观测数据：表4-11大鼠静脉注射625mg/kg PS916后血浆PS916浓度(jig/mL)Table 4-11 Plasma concentration of PS916 after an int

9、ravenous injection of 6.25 mg/kg to rats(pg/mL)Tim 金、123456XSD0.083142.32109.96139.25118.5790.28108.64118.1719.820.16769.8649.4975.2964.4747.4268.5362.5111.440.2525.0230.6619.29363324.3634.472&366.570.510.3913.1514.916.879.0416.2613.443.1817.66338.739.95.466.877.4&1.6323.39414.674.882.433.69

10、3.860.9043.223.373.673.663.263.183.390.2282.213.153.512.982.262.562.780.52121,7&1.912.061.99L851.981.93OJO240.890.950.931.020.881.050.950.07据此绘出如下图形:(二EUOIIeJPU ouoo &U1S&ICIt(h)图4M不同剂PS9I6大鼠静脉给药后乎均药时曲线图比较得，与建立的模型模拟岀的一次性服药的情况类似，模型的模拟性较好。2、等间隔服药模型的检验：査阅文献资料得到等间隔给药后体内药物弄度的变化如下图：01002 00300

11、Time(r/miii)与所建立的数学模型进行比较可以看出模型的符合度较髙。3、恒速静脉滴注的模型检验： 査阅文献资料得到：400肝肿瘤患者静脉恒速滴注高剂虽：(130mg/nr ) L-H0P和GEH后血浆中浓度的测定:将数据进行绘图得到:可以看出在一泄时间以后，药物的浓度就稳左，与上面所建立的数学模型符合较好。八、模型优缺点上文的讨论都是将人体看作是一个较为简单的房室模型，从而将问题简化而忽略诸多的 次要因素，只考虑主要因素对人体内药物浓度变化的影响。这种方法有利于对主要因素的研 究从而较易的观察总结规律。但是实际上人体的实际情况远远要比这种房室模型复杂的多，药物进入血液，通过血液 循环药

12、物被带到身体的各个部位，又通过交换进入各个器官，而并不是单一的简单的输出和 输入。因此，要建立更接近实际情况的数学模型就必须正视机体部位之间的差异及相互之间 的关联关系，这就需要多房室系统模型。九、参考文献1 周义仓等，数学建模实验，西安交通大学出版社，第2版，2007年8月2 朱旭等，MATLAB软件与基础数学实验，西安交通大学岀版社，2008年1月。3 宋晓坤，奥沙利钳腹腔热灌注安全性考察及药代动力学研究，天津医科大学，2008年5 月。4 吕志华，海洋多糖药物PS916的荧光标记及其药代动力学研究，中国海洋大学，2008年 6月。5 陆佳契，韬内动脉灌注甲氨蝶吟在大鼠体内各组织药物浓度监测，复旦大学，2007年4 月。十、附录1、一次性注射情况的程序: t=0:0. 1:80x=0. l*exp(-0. l)*t)plot (t, x, ' linewidth', 2)运行后的图形如下：2、等间隔服药的情况程序:a=0. 1T=8k=0. 1n=0Qap二exp (-k*n*T)t=0；x=ahold onwhile t<100t=t+0. 01;if t<(n+l)*Tx=x*exp (-0. 01*k) elseif t>=(n+l)*Tn=n+l