微积分方法建模药物在体内的分布与排除数学建模案例分析.docx

§10 药物在体内的分布与排除 药物进入机体后，在随血液输送到各器官和组织的过程中，不断地被吸收、分布、代谢，最终排出体外.药物在血液中的浓度 ?? g mv?称血药浓度.血药浓度的大小直接影响到药物的疗 效，浓度太低不能达到预期的治疗效果，浓度太高又可能导致中毒、副作用太强或造成浪费 .因此研究药物在体内吸收、分布和排除的动态过程，对于新药研制时剂量的确定、给药方案设计等药理学和临床医学的发展具有重要的指导意义和实用价值. 为了研究目的，将一个机体划分成若干个房室，每个房室是机体的一部分，比如中心室和周边室.在一个房室内药物呈均匀分布，而在不同的房室之间按一定规律进行转移.如果要求的精度 不是太高的情况下，可以只考虑一室模型. 模型假设 1.药物进入机体后，全部进入中心室（血液较丰富的心、肺、肾等器官和组织），中心室的容积在给药过程中保持不变. 药物从中心室排出体外，与排除的数量相比，药物的吸收可以忽略. 药物排除的速率与中心室的血药浓度成正比. 模型构成与求解 f (t)0给药中心室kc( f (t) 0 给药 中心室 k c(t), x(t) 排除 V 0 c?t ?  中心室血药浓度 x?t ? 中心室药量 一、求解各种给药方式下血药浓度变化情况 V 中心室容积 k 排除速率系数 ?上述各量间有关系 x ? f ? 0 ? ?t ? kx ? ? 即 x? kx ? f t 0 又 x?t ?? Vc?t ? ? ? ? ? f ?t? ?? 得方程 c t 1、 快速静脉注射 kc t ? 0 V （1） c t c 0 设给药量 D ,则初始条件c?0?? D ， f V 0 ?t ?? 0 的解为 c?t ?? 2、恒速静脉注射 D e ?k t V 0 t 设持续时间为? ，注射速率为k ，则有 0 f ?t ?? k ，初始条件c?0?? 0 ， ?0 ? t ? ? ? 0 0 c(t) 1 / 3 0 t PAGE PAGE 2 / 3 f ?t ?? 0 ， c?? ?? k0 ?1 ? e?k ? ?， ?t ? ? ? ? ? ? ? ? 0? ? 0 ? ? ? ?（1）的解为 c?t ? ? ? ? 1 ? e?k t 0 ? t ? ?  （3） ? 3、口服或肌肉注射 ? k0 ?Vk 1 ? e?k ? e?k t ?? t ? ? 在药物输入中心室之前先有一个将药物吸收入血液的过程，可以看作为有一个吸收室，药物由吸 收室进入中心室，药物由吸收室进入中心室的转移速率系数记成k ，给药量 D ，吸收室药量 x ?t ?， 1 0 则 f ?t??k x ?t? 吸收室 0 1 0 ?? ? ? ? x x x ?t? 0  ? ?? 中心室 ? 0 t ? ?k x 1 0 ? x t De?k1t ?? ??x ? 0 ?0?? D 0 c t 于是 f 0 ?t ?? k 1 De?k1t ，初始条件c?0?? 0 ，（1）的解为 ? ? k D ? ? ? ? c t ? V ?k1 ? k ? e?k t ? e?k1t 1 ， k ? k 1 （4） 0t 0 t 1、快速静脉注射下估计k 在t ? 0 时刻快速注射剂量为D 的药物后，在一系列时刻t i  ?i ? 1,2, , n? 从中心室取血样获得血药浓度 c?t i 值： ?.由（2）式反解出几个 k?t i ?，取算术平均值就得到 k 的估计 2、恒速静脉滴注下估计k k ? 1 ?n n i?1 k?t ? i 方法与1 类似， t i 取在?0,? ?内，用（3）中第一式反解k ?t ? i i 3、口服或肌肉注射下估计k 和 k1 k D 因为k1 ? k ，记 A ? V ?k1 ? k ?，于是当t 充分大时（4）近似为 1 c?t ?? Ae?k t 或 ln c?t ?? ln A ? kt 对于适当大的t i 和测得的相应c?t i ?，用最小二乘法估计出k 和ln A ，从而再由 k D 就可以估计出 A ? V ?k1 ? k ? 1 ?AVk ? k1 AV ? D 模型分析 当要求精度较高时,可采用二室甚至多室模型,例如二室模型图示如下，这时的机理分析和参数估计都比一室模型难度更大,需要建立微分方程组来进行分析. f f ?t? 中心室 排除 k 周边室 0 给药 c t , x t ? ? ? ? 12 1 1 c t , x t ? ? ? ? 2 2 V 1 k V 2 21 k 13