全国大学生数学建模竞赛2004优秀论文：C、D题.doc

C题之一（全国一等奖）
酒精在人体内的分布与排除优化模型
桂林工学院，袁孟强，王哲，张莉
指导教师：数模辅导组

摘要：酒精进入机体后，随血液运输到各个器官和组织，不断的被吸收，分布，代谢，最终排除体外。为了研究酒精在体内吸收，分布和排除的动态过程，以及这些过程与人体反应的定量关系，本文建立了一个酒精在人体内的分布与排除优化模型，在药物动力学的一室模型的基础上，进行优化，改进，分别建立了酒精在人体内分布的房室模型和房室模型，以及酒精在人体内的静态排除模型和动态排除模型，导出模型的体液酒精浓度的状态函数，用常数交叉拟合方法，采用VB编写程序，得到两个重要系数和。根据此模型，计算的体液酒精浓度理论值与实验值十分相符，并很好地解释了给出的所有问题，得到一些有价值的结论。

关键词：房室模型，排除模型，体液酒精浓度，动态和静态的转换

酒精在人体内的分布与排除优化模型

一、 问题的重述
国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。
大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？
参考下面给出的数据（或自己收集资料）建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：
1. 对大李碰到的情况做出解释；
2. 在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：
1） 酒是在很短时间内喝的；
2） 酒是在较长一段时间（比如2小时）内喝的。
3. 怎样估计血液中的酒精含量在什么时间最高。
4. 根据你的模型论证：如果天天喝酒，是否还能开车？
5. 根据你做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机如何驾车提出忠告。
参考数据
1. 人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。
2. 体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克／百毫升），得到数据如下(表—1)：

时间(小时)
0.25
0.5
0.75
1
1.5
2
2.5
3
3.5
4
4.5
5
酒精含量
30
68
75
82
82
77
68
68
58
51
50
41
时间(小时)
6
7
8
9
10
11
12
13
14
15
16

酒精含量
38
35
28
25
18
15
12
10
7
7
4

二、 模型假设
1、 酒精的转移速率，及向体外排除的速率， 与该室的血酒浓度成正比。
2、酒精的转移速率，及向体外的排除速率，与时间有关，与空间（人体的各个部分）无关。
3、中心室与体外有酒精交换，及酒精从体外进入中心室，最后又从中心室排出体外。与转移和排除的数量相比，酒精的吸收可以忽略。

三、模型建立与求解
房室模型Ⅰ（在短时间内喝下酒精量为）
在短时间内喝下酒精量为，酒精进入胃，人体吸收酒精，然后排除出体外。吸收酒精的过程相当于酒精进入体液（中心室）的过程，全过程可以简化为下图：

胃室

中心室

排除体外

建模过程：
——短时间内进入胃的酒精；
——为胃室（吸收室）进入中心室的转移速率系数（由人体机能确定的常数）；
——是t时刻胃室的酒精；其微分方程为：
　　　　 （１）
——是t时刻进入中心室的酒精，其微分方程为：
（2）
酒精进入中心室的速率为： （3）
将方程（1）的解代入（3）得： （4）

房室模型Ⅱ（在较长一段时间内喝酒）
假设在较长的一段时间内喝下的酒是匀速进入胃室，则简化如下图：
常数

胃室

排除
中心室

建模过程：
——为酒精进入胃的速率： ,为喝酒时间。
——为酒精从中心室排除体内的速率
——为酒精进入中心室的速率
——为胃室进入中心室的转移速率（由人体机能确定的常数）
——为是酒从中心室向外排除的速率系数。
——是t时刻胃室的酒精，微分方程为：
（5）
（6） ——是t时刻进入中心室的酒精
将方程（5）的解代入（6）得：
　　　（7）
　　　　　　　　　　（8）
静态排除模型Ⅰ
与房室模型I配套的静态酒精排除模型Ⅰ
——中心室的血酒浓度；
V——人体体液量和人体血液量；
酒精进入中心室的速率：
——中心室的酒精量；微分方程为：
（9）

——酒精从中心室向体外排除的速率系数（由人体机能确定的常数）
由方程（9）得： （10）
对应的通解为：
微分方程的解为：

令 得：.
.
根据参考数据表——1，已知：短时间内进入胃的酒精,人体体液量V和一批实验数据（，）。用交叉常数拟合原理在VB环境中编写程序，利用该程序算出两个重要系数和，给出模型的状态函数.
若初始值设为，则

动态排除模型
与房室模型Ⅱ配套建立动态酒精排除模型.
——中心室的血酒浓度；
V——人体体液量和人体血液量；
酒精进入中心室的速率：
——中心室的酒精量；微分方程为：

（11）
（12）
——酒精从中心室向体外排除的速率系数（由人体机能确定的常数）.
由方程（11）得： （13）
对应的通解为：.
将及代入得：
（14）
四、 酒过程的描述
1、在短时间内喝下酒精量为
用静态排除模型描述：.
2、在较长一段时间内喝酒
用动态排除模型描述喝酒的过程，用静态排除模型描述酒后的过程.
即先用状态函数

描述喝酒的过程，然后用状态函数

描述酒后的过程.
3、天天喝下酒精量为
用动态排除模型描述第一天喝酒（喝一小时）的过程，用静态排除模型描述酒后23小时内的的过程，用动态排除模型描述第二天喝酒（喝一小时）的过程，用静态排除模型描述第二天酒后23小时内的的过程，再用动态排除模型描述第三天喝酒（喝一小时）的过程，用静态排除模型描述第三天酒后23小时内的的过程……

……

五、参数的选择
一瓶啤酒的酒精量：.
人体体液：
速率系数： 　　
两小时内慢慢喝下两瓶啤酒的输入率：

六、拟合效果.
中心室中酒精含量（毫克／百毫升）浓度状态函数拟合效果图

效果图显示拟合程度极高，说明参数的选择与客观情况相符合.
七、问题分析
问题1的分析：
大李在中午喝了一瓶啤酒下午6点检查时，利用静态排除模型：

由浓度状态函数在时的浓度：,可知大李此时符合新的驾车标准.
紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车

回家.检测时距晚饭喝的那瓶啤酒已过了八小时，其胃中酒精及体液中的酒精含量分别为：

由浓度状态函数可知，检查时被定为饮酒驾车.
之所以被定为饮酒驾车，关键是此时方程中的初始条件而不是第一次喝酒的.
问题2的分析：
在很短时间内喝了3瓶啤酒,问多长时间内驾车就会违反新的驾车标准. 根据静态排除模型Ⅰ：
3瓶啤酒的酒精量：；
浓度状态函数为 .
代入数据得：. （15）
若违反新的驾车标准,则：,根据人体体液酒精浓度曲线图可知，喝酒后约在t=0.07小时（42秒）与t=11.24小时之内驾车会违反新标准。将t=0.07，t=11.24分别代入方程（15）检验得出：

通过验证，证明观测值基本接近实际值。

问题3的分析：
在较长一段时间（比如2小时）内喝3瓶啤酒，多长时间内驾车就会违反新标准.
假设匀速喝酒，则
此过程分两阶段：（1）在喝酒过程中，多少时间后驾车会违反新标准。
（2）喝酒之后，多长时间内驾车会违反新标准。
根据动态模型：

（16）

（1） 由人体体液酒精浓度曲线观测出：当t大于0.62小时（约37.2分钟）时，体液酒精浓度大于毫克，把已知酒精数据，，代入方程（15）检验得：

通过数据验证，证明观测值基本接近实际值。

（2）由阶段（1）的已知数据算出，当t=2小时（即停止喝酒时）人体酒精为：
此时，根据胃里剩余的酒精方程：
得出：
以2小时作为零时刻，设状态初始值为,则取，此时静态排除模型为：

根据此模型，可得曲线：

可观测出：当时，t(t>2)的值约在11.67附近。再将t=11.67代入进行检测得：
（毫克/百毫升）
通过数据验证，证明观测值基本接近实际值。
结论：若连续2小时均匀喝下3瓶啤酒，则在0.62与13.67小时之内驾车会违反新标准。

问题4的分析：
估算出血液中酒精含量何时达到峰值。
最高值的估计要分两种情况讨论：

1、 据静态排除模型（一次喝完的情况）
酒精进入机体之前，胃里的酒精浓度为0，当喝酒后，胃里的酒精浓度逐渐增大，因为存在着酒精进入中心室的过程，根据人体机理可知，。当胃里的酒精浓度等于中心室的酒精浓度时，，此时血液中酒精含量达到最高值。

由静态排除模型得：
若要使的值达到最高，则要求即：

代入已知数据得：

算出：t=1.36小时
即：喝酒1.36小时后，血液中酒精含量最高。
（1） 根据动态排除模型（酒是在很长一段时间内喝完）
喝酒时，由于酒量在不断的增加，根据人体机理可知，胃中的酒精浓度始终是大于中心室的酒精浓度，即，易知血液中的酒精浓度是一个递增的函数。那么在这阶段中，血液的酒精含量的最高值为停止喝酒时血液中的酒精含量。
运用动态排除模型：

可计算出最值。
喝完酒后，胃中的酒精浓度仍然大于中心室的酒精浓度，此时，.随着时间的推移，胃中的酒精浓度在不断的减少，而血液的酒精浓度不断增加，这就使得当时，血液的酒精浓度达到最高值（即），若代入具体数据就可算出最大值。

问题5的分析：
根据我们建立的模型论证：如果天天喝酒，是否还能开车？
天天喝一瓶酒的数学模型见图

一个人天天喝酒，如果每天只喝一瓶啤酒第二天即12个小时后体内的血酒浓度远远小于20毫克/百毫升所以仍能开车；

天天喝两瓶酒的数学模型见图：

一个人天天喝酒，如果每天只喝两瓶啤酒第二天即12个小时后体内的血酒浓度远远小于20毫克/百毫升所以仍能开车；
天天喝三瓶酒的数学模型见图：

如果每天喝三瓶啤酒，12小时后体内的血酒浓度为18.920毫克/百毫升，比新的国家标准低但很接近，此时仍可以开车却比较危险

附表（1）如下：
体内的
血酒浓度

时间
(t)
摄入
的酒
精量

一瓶啤酒的酒精含量
24192
(mg)
两瓶啤酒的酒精含量
48384
(mg)
三瓶啤酒的酒精含量
72576
(mg)
四瓶啤酒的酒精含量
96768
(mg)
五瓶啤酒的酒精含量
120960
(mg)
12
5.7301
12.6
18.9
25.2
28.604
24
0.5198
1.396
1.5595
2.0793
2.5991

注：据资料表明人体酒精浓度达到200(毫克/百毫升)—300(毫克/百毫升) 会使人昏迷,一般人喝五瓶啤酒（含酒精量120960毫克）后，（据静态排出模型）在一小时左右浓度达到213.6729毫克/百毫升已使人昏迷,因而研究五瓶啤酒以上的酒精量已没有意义。
结论：每天喝三瓶以下的啤酒第二天不违反新的国家标准仍可以开车，而喝四瓶到五瓶就不能再开车了，五瓶以上大多数人都会昏迷

问题6的分析：
通过以上论述，特给爱喝点酒的司机提个醒，希望他们能平安。
科学研究表明，人体内每百毫升血液中酒精含量达到20毫克，会出现头晕脑胀、兴奋健谈；达 60～80毫克时，感情容易冲动，步态不稳；达120～160毫克时，神经进入抑制状态，开始昏睡；达200～400毫克时，意识朦胧，呈木僵状态，如果达400～500毫克时，就可能导致脑损伤和呼吸麻痹，从而死亡。
根据状态函数与静态排除模型的分析方法可得出：
喝完一瓶啤酒经5.74小时后，血液的酒精浓度仍为20.0383，超过新的驾车标准 ，为了安全，建议司机饮酒八小时后驾车。 .
喝完两瓶啤酒经9.2小时后，血液的酒精浓度仍为20.063，饮酒的司机在11小时内不宜驾车。
喝完四瓶啤酒经12.66小时后，血液的酒精浓度仍为20.086，饮酒的司机在14个小时内不宜驾车。
综上所述，为安全起见，对于爱饮酒的司机，每天的饮酒量不宜超过两瓶啤酒。若饮酒量超过四瓶，第二天司机不宜驾车。
八、模型评价
1、 本模型绘出状态函数后，极好的解决了喝酒过程的全程数学描述，定量的解决了所提出的问题。
2、 本模型微分方程的待定系数求解方案没有利用现成的方法。而是自行编程解决，给此类模型提供了一个利用的程序。
3、 本模型通俗易懂，为大学生所接受，但反复运用状态函数，很有特色，尤以初值问题多次循环反复，妙用至极。

4、 有了本模型，人们对喝酒过程有一个较为精确的定量认识。从事有一定危险工作的人们应该如何正确的喝酒的依据。
5、 本模型可推广到有毒物质在人体的分布与排除。这对在有毒物场所工作的人起到一定的提醒作用。

参考文献：
1、[ISBN 7-04-004505-2 ] 姜启源 ，《数学模型》（第二版），北京：高等教育出版社 ，
1993年8月
2、[ISBN 7-04-011943-9] 徐全智 杨晋浩，《数学建模》，北京：高等教育出版社，
2003年七月

附件
常数拟合流程图

开始

输入，D，V，M，
t，N，

i=1

i<N

i=i+1
Y

N
i=1

|

N N

Y

输出结果

i=i+1
i<N
Y
Y

N N

说明：是精度，D是酒精量，V是人体体液量，M是实验的酒精含量数据，t是实验的时间数据，是代定常数初值，是步长值。
源程序
Dim T(23) As Single '实验值，时间（小时）
Dim M(23) As Integer '酒精含量

Private Sub Form_Load
T(1) = 0.25
T(2) = 0.5
T(3) = 0.75
T(4) = 1
T(5) = 1.5
T(6) = 2
T(7) = 2.5
T(8) = 3
T(9) = 3.5
T(10) = 4
T(11) = 4.5

T(12) = 5
T(13) = 6
T(14) = 7
T(15) = 8
T(16) = 9
T(17) = 10
T(18) = 11
T(19) = 12
T(20) = 13
T(21) = 14
T(22) = 15
T(23) = 16

M(1) = 30
M(2) = 68
M(3) = 75
M(4) = 82
M(5) = 82
M(6) = 77
M(7) = 68
M(8) = 68
M(9) = 58
M(10) = 51
M(11) = 50
M(12) = 41
M(13) = 38
M(14) = 35
M(15) = 28
M(16) = 25
M(17) = 18
M(18) = 15
M(19) = 12
M(20) = 10
M(21) = 7
M(22) = 7
M(23) = 4
End Sub

Public Function C1(ByVal D0 As Single, ByVal V As Double, ByVal K01 As Single, ByVal K10 As Single, ByVal T As Single) As Single
C1 = D0 * K01 * (Exp(-K01 * T) - Exp(-K10 * T)) / V / (K10 - K01)
End Function

Private Sub cmdOK_Click

Const N = 23
Const nStep = 0.01 '如果设得太小，计算时间会很长。当然，精确度会更高。
Dim i As Integer

Dim MaxK01 As Single
Dim MaxK10 As Single
Dim initK01 As Single

Dim D0 As Single
Dim V As Double
Dim K01 As Single
Dim K10 As Single
Dim Ei As Single

Dim IsPass As Boolean
Dim nTimer As Single

txtV.Text = Trim(txtV.Text)
If Not IsNumeric(txtV.Text) Then
MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtV.SetFocus
Exit Sub
End If
If txtV.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtV.SetFocus
Exit Sub
End If

txtD0.Text = Trim(txtD0.Text)
If Not IsNumeric(txtD0.Text) Then
MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtD0.SetFocus
Exit Sub
End If
If txtD0.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtD0.SetFocus
Exit Sub
End If

txtK01.Text = Trim(txtK01.Text)
If Not IsNumeric(txtK01.Text) Then

MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtK01.SetFocus
Exit Sub
End If
If txtK01.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtK01.SetFocus
Exit Sub
End If

txtK10.Text = Trim(txtK10.Text)
If Not IsNumeric(txtK10.Text) Then
MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtK10.SetFocus
Exit Sub
End If
If txtK10.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtK10.SetFocus
Exit Sub
End If

txtE.Text = Trim(txtE.Text)
If Not IsNumeric(txtE.Text) Then
MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtE.SetFocus
Exit Sub
End If
If txtE.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtE.SetFocus
Exit Sub
End If

txtTime.Text = Trim(txtTime.Text)
If Not IsNumeric(txtTime.Text) Then
MsgBox "请输入数值。", vbOKOnly + vbInformation, "数据拟合"
txtTime.SetFocus
Exit Sub
End If
If txtTime.Text <= 0 Then
MsgBox "请输入大于0的数值。", vbOKOnly + vbInformation, "数据拟合"
txtTime.SetFocus
Exit Sub

End If

On Error GoTo Err1

txtC1.Text = ""
txtC1.Refresh
txtK01v.Text = ""
txtK01v.Refresh
txtK10v.Text = ""
txtK10v.Refresh

D0 = txtD0.Text
V = txtV.Text
MaxK01 = txtK01.Text
MaxK10 = txtK10.Text
Ei = txtE.Text

initK01 = MaxK01
K01 = MaxK01
K10 = MaxK10
nTimer = Timer

Me.MousePointer = 11

While Timer - nTimer < txtTime.Text
If initK01 = MaxK01 Then
While K10 > 0
K01 = MaxK01
While K01 > 0
IsPass = True
For i = 1 To N
If K01 = K10 Then
IsPass = False
Exit For
End If
If Abs((C1(D0, V, K01, K10, T(i)) - M(i)) / M(i)) >= Ei Then
IsPass = False
Exit For
End If
Next i

If IsPass = True Then
GoTo Pass:

Else
K01 = K01 - nStep
End If
Wend
K10 = K10 - nStep
Wend

MaxK01 = MaxK01 + 1
MaxK10 = MaxK10 + 1
K01 = MaxK01
K10 = MaxK10
End If

While K10 > MaxK10 - 1
K01 = MaxK01
While K01 > 0
IsPass = True
For i = 1 To N
If K01 = K10 Then
IsPass = False
Exit For
End If
If Abs((C1(D0, V, K01, K10, T(i)) - M(i)) / M(i)) >= Ei Then
IsPass = False
Exit For
End If
Next i

If IsPass Then
GoTo Pass:
Else
K01 = K01 - nStep
End If
Wend
K10 = K10 - nStep
Wend

While K10 > 0
K01 = MaxK01
While K01 > MaxK01 - 1
IsPass = True
For i = 1 To N

If K01 = K10 Then
IsPass = False
Exit For
End If
If Abs((C1(D0, V, K01, K10, T(i)) - M(i)) / M(i)) >= Ei Then
IsPass = False
Exit For
End If
Next i

If IsPass Then
GoTo Pass:
Else
K01 = K01 - nStep
End If
Wend
K10 = K10 - nStep
Wend

If K10 < 0 Then
MaxK01 = MaxK01 + 1
MaxK10 = MaxK10 + 1
K01 = MaxK01
K10 = MaxK10
End If
Wend
Me.MousePointer = 0
MsgBox " 无结果。 ", vbOKOnly + vbInformation, "数据拟合"

Exit Sub
Pass:
Me.MousePointer = 0
txtC1.Text = CInt(C1(D0, V, K01, K10, T(1)))
For i = 2 To N
txtC1.Text = txtC1.Text & "," & CInt(C1(D0, V, K01, K10, T(i)))
Next i
txtK01v.Text = Format(K01, "0.00")
txtK10v.Text = Format(K10, "0.00")
Ms