一、问题分析
小孩误用药,药会先进入胃肠道,部分进入血液系统,血液系统中会对药物做一定速率地排出。但当药物浓度在血液系统中过大时,会有生命危险。(血药浓度达到100 ug/mL会出现严重中毒;血药浓度达到200 ug/mL可能会致命),对于2000mL血液总量的小孩,即200mg与400mg的含量。
二、模型假设
设胃肠道中的药量为 x(t),血液系统中的药量为 y(t),时间 t 以孩子误用药为起点(t=0)。
根据实际资料调查,我们做出下面假设:
1. 胃肠道中药物移向血液系统的转移率与胃药量 x(t) 成正比,比例系数为 λ (>0);
2. 药物总剂量 N 在复用后瞬间进入胃肠道;
3. 血液系统中药物的排除率与血液药量 y(t) 成正比,比例系数为 μ (>0);
4. 血液系统在t=0时无药物(由前面假设知开始时均在胃肠道);
5. 该药物的被吸收的半衰期为 t1 h,排除的半衰期为 t2 h;
6. 孩子的血液总量 L = 2000 mL。
三、模型建立
根据假设,t=0时刻,x(0) = N,
随药物从胃肠道向血液系统移动(由假设与当前胃肠道的药量呈线性关系),
x(t) 满足微分方程: (1)
血液系统的药量 一部分从胃肠道转移而来(+)(相当于上面微分方程输入的),一部分随肾脏排出到体外(-)(由假设与当前血液系统的药量呈线性关系),
y(t) 满足微分方程: (2)
其中,参数 λ,μ由假设中药物的半衰期(t1,t2 )计算确定。
四、模型求解
方程(1)求解得到 (3)
根据吸收的半衰期 t1 求 λ,x(t1) = N/2(药物从胃肠道被血液被吸收了一半)
结合方程(3)即可解得 λ。
方程(2)求解得到 (4)
根据排除的半衰期 t2 求 μ,
但此时需要用的y(t)方程是与x(t)无关,即只与排出的y有关, 有
,与上面解λ同理,结合 y(t2) = N/2求解μ。
到这里,就能得到 x(t) 与 y(t) 的具体方程表达式,只需代入数据得到数值解。
施救方案,若在t = tp时刻,若采取体外血液透析(医学方法加快药物排除),那么在此时,方程(2)改变(μ改变),血液系统中药物浓度设为 z(t),有
μ与具体的医学治疗方法相关。
五、模型分析
六、模型检验
七、模型应用
对相同的药物(及有相同的吸收、排除半衰期 λ,μ),达到血液中药物最大浓度的时间(只与λ,μ相关)是相同的。
模型与建模方程来源于《数学模型》第五版。