液体输注动力学 模型的讨论.pptVIP

液体输注动力学模型的讨论 数 学 系 谢晓天 朱一超 内容提要 问题的提出 一室模型的给出 一室模型假设 一室模型建立 一室模型参数估计 改进后的参数估计 对实际问题的模拟 内容提要 模型的改进——二室模型 二室模型假设 二室模型建立 参数估计 改进的参数估计 实验数据测试 附录 问题的提出 液体通过静脉输入机体后，在随血液输运到各个器官和组织的过程中，不断被吸收、分布、代谢，最终排出体外。研究液体在体内吸收、分布和排除的动态过程，对新药研制、剂量确定、输液方案设计等药理学和临床医学问题都具有重要的实用价值。在这样的背景下，瑞典人Lars Stahle在1997年提出了液体输注动力学模型。 问题的提出 液体输注动力学模型用类似药代动力学的方法，研究液体在体内吸收、分布和消除的过程。但它和传统的药代动力学模型不完全一样： 前者假定在输液过程中，液体进入一个可扩张的液体容量空间。依靠测定血浆容量的稀释程度，估算液体容量空间的变化。 后者假定液体容量空间是保持不变的，依靠测定血药浓度，估算药物在血液中的浓度变化 。 一室模型的给出 输液时，人体机体就像一个气球，具有可扩张性，也就是可扩张的液体容量空间（volume-of-fluid-space, VOFS）下面介绍一个模拟其随时间变化的模型——一室模型。 一室模型的假设 输液过程中，液体进入一个可扩张的液体空间容量为v，液体以恒定的速度ki持续输入 可扩张的液体空间有一个靶容量V，这个靶容量是机体极力维持的 V中液体的消除通过两种机制：kb是通过呼吸和排尿实现，大约0.5ml/min，而kr是由肾脏功能决定，与稀释程度有关。 一室模型的假设 基于以上假设，液体的变化规律如下图所示： 一室模型的建立 根据假设我们有如下微分方程： 在输液过程中： 其中：ki(t)表示输液速度，s表示的输液时间 一室模型的建立 根据假设我们有如下微分方程： 在输液过程中： 其中：ki(t)表示输液速度，s表示的输液时间 一室模型的建立 针头拔出后，有如下微分方程： 一室模型的建立意义 我们可以通过对方程中一组（v-V）/V的测量数据，来计算参数kr与靶容量V的大小。 注：测量通过血浆稀释度的测量来完成 的，详情请见《British Journal of Anaesthesia 1997》。 一室模型的解 为简化问题，并从实际出发，我们可以假定注射速度ki为恒定的，我们可以用分段函数来表示注射速度： 注：D为注射量，s为注射时间 一室模型的解 用w=(v-V)/V代入方程（1），（2），得到： 一室模型的解 我们得到以上方程的解析解： 一室模型的参数估计 在没有参数的情况下，我们利用w的表达式，在恒速静脉滴注方式下，根据假定参数来计算一组w的值，并画图。 假定参数如下： kr=50， v=5000, s=30, D=1800。 一室模型的参数估计 通过观察我们发现，在注射过程中w在增加，当针头拔出后，w减少并趋于稳定。 一室模型的参数估计 我们将拟合图象画出来。拟合非常精确。我们基本看不到原来的黑点了。得出：kr=50.0 V=4999.99 一室模型的参数估计 通过这种方法，我们可以得到模型的参数估计，但是这些理论上估计是无扰动的，对于微小的扰动，我们的拟合方法是否合适呢？ 我们知道，根据w的表达式，w在微小扰动下的误差是可控的，但是指数可控，这对于实际问题来说，其实偏差还是挺大的，因此我们可以改进参数估计方法。 改进的的参数估计 我们仍然按前面假设的参数，画出w的图象。 改进的的参数估计 我们对w的值用生成随机数的方法来进行微小扰动，得到右图： 改进的的参数估计 得到： kr=50.1071， V=5074.3。 误差在5%以内。 改进的的参数估计 为了提高精度，我们可以用随机生成50组数据，并取平均的方法来做，下图是50组的kr平均值柱状图表： 改进的的参数估计 让我们再来看一下V的50组数据的柱状图： 改进的的参数估计 通过多次计算取平均值的方法，我们得到： kr=50.0，V=4998.9 精度比以前又有所提高，这有助于我们解决实际问题。 对实际问题的模拟 这是一组临床的数据： 对实际问题的模拟 将数据代入进行拟合得到： kr=156.8 V=4983.6 对实际问题的模拟 我们再看一组临床数据： 对实际问题的模拟 将数据代入进行拟合得到： kr=153.4 V=3354.4 小结 一室模