摘
要：目的
寻找中药水提液浓缩过程瞬时能耗计算简便方法。
方法
以玉屏风散为模型药物，测定其水提液浓缩过程中多种溶液环境特征参数的动态变化数据，构建中药溶液环境特征参数白利度（简便、易测、准确）与浓缩能耗计算需要的相关参数的相关性模型。
结果
白利度与固含率、密度、比热容、导热率等具有优越的线性相关，相关系数（
R
²
）均高达
0.99
左右。
结论
可通过测定白利度及构建与其他物理参数如比热容（
C
p
）、导热率（
k
f
）等关系模型，计算浓缩的动态能耗；为中药制药过程工程原理研究开辟新思路。
传统中药制药过程中，提取、浓缩、干燥等工段都需要使用大量蒸汽
[1]
。其中，中药浓缩工段的能耗惊人，浓缩工段所用蒸汽量占总蒸汽量
60%
以上
[2]
。膜浓缩（反渗透、膜蒸馏等）比一般热法浓缩如单效、双效蒸发节能
50%
[3]
。而近期
Weaver
等
[4]
发现，使用反渗透预浓缩技术可比单独用热效浓缩技术节省能耗
85%
以上。
物料浓缩的原理是实现水（或其他溶剂）与可溶性物质（溶质）的分离，目前常见的产业化浓缩技术有热法蒸发（单效、多效蒸发），膜技术（纳滤、反渗透、膜蒸馏等），均属于物理分离过程。热法、膜法浓缩是海水淡化的常见手段，已有比较成熟的能耗计算与评估方法。但中药制药领域因缺乏相关决策研究，无法采用这些现有的方法。如研究表明中药与相关天然产物的浓缩过程远比海水淡化过程复杂，其主要特征是物料的溶液环境特征参数发生了明显的动态变化，因此海水淡化能耗计算方法对于中药水提液体系存在较大误差
[5-6]
。因此，亟需提出一种将中药溶液环境特征参数变化融入到浓缩能耗计算的新方法。
溶液环境是指溶液体系所具有的电导率、折光指数（实际应用中常换算成白利度）、黏度、
pH
值、离子强度等特征性质。在中药应用领域，精制、浓缩等制药工艺过程所面对的溶液环境体系即为中药液体物料。上述溶液环境特征参数，虽然彼此间存在大量非线性、高噪声、多因子的复杂关系，但借助理论化学对简单物质研究的成果，可从中抽提出若干参数和概念，进而运用人工智能技术，从已知数据和实验事实中抽提规律性，用于简化中药制药过程所采用膜分离等先进技术的工艺优化设计。如张刘红等
[7]
已建立渗透压与电导率相关性模型，借助中药物料浓缩过程中电导率的变化推测其他特征参数的变化，可实现中药溶液环境各相关特征参数的在线实时监测，为中药制药的智能监控提供更简便、经济的手段，并为中药绿色节能膜浓缩技术的推广应用提供有力的计算依据与理论支持
[7-10]
。
中药物料的浓缩过程还涉及到其浓度变化所带来的比热容、固体含量变化以及密度等溶液环境表征参数的变化。其中，折光指数是国内外评判天然产物及中药浓缩物料常用的检测指标
[11-13]
，也是中药制药工业检测物料比重所用传感器不可或缺的重要参数
[14-15]
。利用白利度的变化推测其他物理参数的变化，对中药制药在线检测和浓缩所需能耗计算可达到事半功倍的作用。但过去的研究中并未对特征参数如何影响中药物料浓缩工段的能耗变化做出详细说明（尤其是中药物料的膜浓缩能耗计算）。
多年来中药制药工程设计中，因缺乏基本的中药物性数据（如不同中药物料的密度、黏度、表面张力、导热系数、扩散系数等）
[10]
，往往凭经验或采用经验方式估算，甚至用相近物质的物性代替，其结果必然导致工艺技术选择或设计的失真甚至失败。从而使中药生产的规范化、现代化难以付诸现实。本实验拟以玉屏风散为模型药物，探讨浓缩过程对中药溶液环境特征参数的影响及各特征参数的相关性，指出相关特征参数变化与中药物料浓缩工段能耗的关系，为中药浓缩工艺的瞬时能耗计算提供一种简便的模式，同时也为探讨中药制药过程工程原理开辟新视野。
1
仪器与材料
1.1
药材
防风购于安徽道源堂中药饮片有限公司，批号
200501
，白术购于毫州市永刚饮片厂有限公司，批号
200203
，蜜炙黄芪购于赤峰荣兴堂药业有限责任公司蒙中药饮片厂，批号
2005046
，上述
3
味药材经本课题组郭立玮教授鉴定，分别为伞形科防风属植物防风
Saposhnikovia divaricata
(Turcz.) Schischk.
的干燥根；菊科苍术属植物白术
Atractylodes macrocephala
Koidz.
的干燥根茎；豆科黄芪属植物膜荚黄芪
Astragalus membranaceus
(Fisch.) Bge.
干燥根的炮制（蜜炙）加工品。
1.2
仪器
BSA224S-CW
型电子分析天平，赛多利斯科学仪器（北京）有限公司；
LC-DRT-948
数显折射仪、
NDJ-5S
数显黏度计，上海力辰邦西仪器科技有限公司；
TP101
笔式探针数显温度计，测量范围
−50
～
300
℃
，苏泊尔
C22-IH66E8
电磁炉，
2200 W
，浙江苏泊尔生活电器有限公司；简悠便携式炆火电陶炉
1200 W
，中山市诺洁仕电器有限公司；
DHG-9145A
电热鼓风干燥箱，上海恒勤仪器设备有限公司；
OM819C
渗透压仪，北京雅森博科科学仪器有限公司。
2
方法
2.1
提取方法
参照“十二五”普通高等教育本科国家级规划教材《方剂学》中玉屏风散方各药味的剂量比
[16-17]
，准确称取药材如下：防风
50.0 g
、白术
100.0 g
、蜜炙黄芪
100.0 g
。用
2500 mL
纯净水（经过实验室反渗透滤膜制备的纯净水）煎煮
1 h
，分离药渣与提取液。再在药渣中再加入
2500 mL
纯净水煎煮
1 h
，分离药渣与提取液。将
2
次提取液均匀混合后，
150
目纱布滤过，收集滤液备用。
2.2
中药物料浓缩实验
采用敞口锅（直径为
24 cm
，容积为
6800 mL
的实验室用小型敞口锅）对上述滤液加热蒸发进行浓缩。采集不同时间段浓缩物料，进行标记，按照“
2.3
”项所述方法检测相应的溶液环境特征参数。
2.3
溶液环境特征参数测试
2.3.1
白利度
量筒量取
100 mL
样品，其中
1
～
2 mL
样品用来测量白利度，分
2
组样品测试读数，由机器分别自动读取
15
次数据，直到数据稳定不变再记录数值。取得的
2
组数值记录下来，并且算得平均值。由于实验室采用的是离线折光仪，需要取样读数，而样品数量仅需
1
～
2
滴，因此样品温度很快就被降至室温。为了统一样品折光度与白利度读数温度，此研究的白利度测试均在室温（
30
±
5
）
℃
下进行。
2.3.2
黏度
待同一批次取得的物料温度降至室温（
30
±
5
）
℃
，用黏度计读取物料的电导率和黏度，读取
3
组数据取得平均值。
2.3.3
密度
将干燥好的烧杯放置天平上清零。从上述同一批次取的物料中，用量取精准量取一定体积的物料并记录体积并倒入干燥的烧杯中。密度测量时物料温度为（
45
±
5
）
℃
。将称有物料的烧杯放置天平上，准确称量其质量并记录。同一批次、浓度物料的密度可用传统密度定义公式算得。
2.3.4
固含率
烘干玻璃培养皿并进行准确称定质量，记为
m
1
。将同一批次取得的物料放入烘干的玻璃培养皿，并准确称定质量，记为
m
2
。将已称定质量的物料与玻璃培养皿放入电热鼓风干燥机中，设置
50
℃
干燥
4 h
以上确保水份完全蒸发。待干燥的物料与玻璃培养皿放置室温后，准确称定质量，记为
m
3
。此样品的固含率可用公式（
1
）计算。
固含率＝
(
m
3
－
m
2
)/(
m
2
－
m
1
)
（
1
）
2.3.5
物料比热容（
C
p
）与导热率（
k
f
）
浓缩中的物料比热容与导热性变化由此部分实验操作得到。比热容是由热量（
Q
）、质量（
M
）和温度变化（
Δ
T
）的比值而定义的
[18]
。
C
p
＝
Q
/
M
Δ
T
（
2
）
在同一组测量实验中准确称量相同质量的水和物料。水与物料在相同时间内得到的热量一样（均被同样加热设备加热），根据已知纯水的比热容
［
C
pw
，
4.18 kJ/(kg·K)
］
，利用同样质量的水和物料在相同时间内的温度变化比，用以下推导得到的公式（
6
）可以算得采集到的物料的
C
p
。
Q
p
＝
Q
pw
（
3
）
Q
p
＝
C
p
(
T
3
－
T
1
)
（
4
）
量取物料被加热前的温度（
T
1
），被加热一段时间（
t
）后量取其即时温度记录为
T
3
；
水被加热前的温度为
T
2
，同理水被加热一段时间
t
后量取其即时温度记为
T
4
，则有：
Q
pw
＝
C
pw
(
T
3
－
T
1
)
（
5
）
联立（
3
）、（
4
）和（
5
）得到公式（
6
）
。
C
p
＝
C
pw
(
T
4
－
T
2
)/(
T
3
－
T
1
)
（
6
）
为控制水和物料都得到相同的热量，实验设计水和物料均被电陶炉（选择
70 W
）加热
240 s
（
t
＝
240 s
）。
另外，根据传统传热学中导热率的定义，一个体系的导热率（
k
）如下
[18]
。
k
＝
Q
Δ
x
/
A
d
T
（
7
）
同理可推得，同一组实验已经称得的相同质量的水和物料都得到相同的输入热量（
Q
）、加热接触面积（
A
）和热量传递距离（
Δ
x
），物料的
k
f
可用已知的水的导热率（
k
w
）根据公式（
8
）计算得到。
k
f
＝
k
w
(
T
4
－
T
2
)/(
T
3
－
T
1
)
（
8
）
k
w
选择
0.640 6W/(m·K)
，水在
50
℃的导热率
2.3.6
渗透压（
π
f
）
按照
YASN OM819
型渗透压仪操作手册执行，由扬子江药业集团南京海陵中药制药工艺技术国家工程研究中心检测。每次测量样品需
10 mL
，因此，每组样品收集
30mL
，共测量
3
组数据，取平均值，得到渗透浓度
b
0
（
mOsm/kg
）。根据公式（
9
）算得渗透压（
kPa
）。
π
f
＝
b
0
RT
（
9
）
R
为摩尔气体常数
8.314 J/(K
∙
mol)
，
T
为开氏温度，所有数据在
25
℃（
298.15K
）
时测量得到
3
结果与分析
3.1
白利度与固含率相关性研究
白利度
[19-20]
最早被定义为溶液中蔗糖的质量分数，后被广义定义为与蔗糖质量分数相等的溶液中可溶性固形物质量分数，被广泛运用在浓缩行业（如果汁浓缩和凉茶浓缩）的物料浓度。白利度一般用折射仪测得折射指数，并换算成
°Bx
或
%
。其中，可溶性固形物也包括溶解性物质与胶态物质。
溶液中固含率定义一般为固体成分质量分数，包括了可溶性与不可溶固形物。因此一个系统里的白利度与固含率不一定相同；一般可认为固含率要比白利度高。本实验中，所得到的同一时间点的白利度与固含率相近（表
1
），而测得的白利度基本比固含率低，并且呈现出较好的线性关系，回归方程为
y
＝
0.011 5
x
＋
0.001 7
，
R
2
＝
0.998 8
，
y
代表固含率，
x
代表白利度。因此，在玉屏风散的浓缩过程研究中，可以利用白利度去预测玉屏风散的固含率。除了经典复方玉屏风散在浓缩过程中具有白利度与固含率的上述关系外，此前
Criscuoli
等
[11]
也考察了
枣汁的膜蒸馏浓缩过程中白利度与固含率之间的关系，得出类似结果。
需要特别指出的是，中药制药工艺往往涉及澄清、精制工段，淀粉、蛋白质等大分子杂质在浓缩工段之前已被去除。因此中药物料中一些大分子物质可能因为加热等问题引起的析出影响可被忽略。
3.2
浓缩过程中中药物料溶液环境表征参数的相关性研究
3.2.1
白利度与密度、黏度的关系
浓缩过程中同一时间点玉屏风散物料的溶液环境特征参数白利度与密度、黏度关系的实测数据分别如表
2
、
3
所示。
根据密度的基本定义，物质的密度可以用质量除以体积算得。中药水提液物料的浓缩过程，可认为是一个去除水，留下固体溶质的过程。一般认为固体溶质的密度比水密度大，因此，除水过程中物料的密度应呈上升趋势。玉屏风散浓缩过程中，物料的白利度与密度呈现优越的线性相关，回归方程为
y
＝
4.266 1
x
＋
1 072.3
，
R
2
＝
0.993 2
，
y
代表密度，
x
代表白利度。
黏度是流体的一种特性，客观地表达了流体对流动的阻力，在中药水提液浓缩体系中则反映了水分子在运动过程中的受阻程度。白利度与黏度是
2
个不同的物理概念，没有直接的关联。本实验所得
到的白利度与黏度的关系方程为
y
＝
3.361 9 e
15.593
x
，
R
2
＝
0.933 1
，
y
代表黏度，
x
代表白利度。根据二者的关系方程，虽不能由白利度准确计算物料体系的黏度，但是能粗略地估算玉屏风散物料的黏度。中药液体物料黏度的数值主要是用来计算流体体系的雷诺数，以预测流体的分类（湍流、层流），为膜蒸馏浓缩中药物料的即时能耗动态计算提供一个重要的选取依据。
3.2.2
白利度与比热容、导热率的关系
从物质比热容定义及公式（
2
）可知，某体系的比热容是指在该体系中单位质量的物质改变单位温度（
1 K
）时需要吸收或释放的热量。而在本研究中的比热容，特指定压比热容（在同一压强下）。
一般认为一个复杂体系中，比热容与混合物质质量呈一定比例。因此在中药溶液环境的复杂体系中，水的去除应该使物料比热容的大小有所改变。浓缩过程中同一时间点的溶液环境特征参数白利度与比热容、导热率的实测数据分别见表
4
、
5
。白利度与玉屏风散溶液环境中比热容的关系方程为
y
＝
0.044
x
＋
4.039 7
，
R
2
＝
0.987 9
，
y
代表比热容，
x
代表白利度；在实验初始白利度为
2.6°Bx
时，溶液的比热容与水的比热容非常接近
［
4.18 kJ/(kg
∙K
)
］
，而比热容随着水分的减少而增加。比热容是能耗测算公式的关键参数，但因为中药水体液物料体系是一种高度多维的复杂系统，其比热容难于测定。因此利用白利度去预测不同浓度溶液下的比热容，对浓缩能耗计算具有填补空白的突破性意义。
导热率（又称导热系数）是一个形容系统热传导能力的物理量，由其计算公式（
7
）可知，导热率即为在相同的热量传导条件下，温差为
1 K
的单位厚度的系统在传热为单位面积单位时间条件下的传热量，导热率与系统的成分组成、状态等有关系。根据公式（
8
）可推知，本研究测试的玉屏风散系统中导热率与系统物料温差成正比。利用白利度与导
热率的关系方程（
y
＝
0.006 7
x
＋
0.619 1
，
R
2
＝
0.9879
，
y
代表导热率，
x
代表白利度），可以成功地预测出玉屏风散在不同浓度溶液下的导热率。
3.2.3
白利度与渗透压的关系
浓缩过程中同一时间点的溶液环境特征参数白利度与渗透压实测数据见表
6
；白利度与渗透压的关系方程为
y
＝
160.54
x
－
366.66
，
R
2
＝
0.965 0
，
y
代表渗透压，
x
代表白利度。
渗透压是与单位体积溶剂中溶质微粒的浓度有关，可以理解一个溶液系统的渗透压是由溶解性溶质产生。因此在玉屏风散的物料体系中，可以观察到白利度与渗透压呈现良好的线性关系，而白利度表示一个溶液体系的溶解性固体物质量分数。溶解性固体物质量分数增加，渗透压也随着增加，符合理论期待值。实验中观察到在白利度高于
10°Bx
时，测得的渗透压比线性方程预测的数值稍微低一些，这可能是因为物料中有些物质经过一段时间（样品由广州寄往南京期间）储存后沉淀析出，导致了测得的渗透压比样品收集时的渗透压稍低。利用液体物料白利度可以较好地推算中药物料的瞬时渗透压，对反渗透膜浓缩在中药物料浓缩中的瞬时能耗变化和总能耗计算有重要意义。
3.3
基于本研究结果的热法浓缩及膜蒸馏浓缩的能耗计算
3.3.1
热法浓缩能耗计算相关性
如热法单效蒸发浓缩中药物料，需要大量蒸汽或者电转化成热能，因此，计算单效蒸发最低所需热能（
Q
T
）的公式如下
[3]
。
Q
T
＝
Q
v
＋
Q
L
（
10
）
Q
v
＝
M
v
(Δ
H
v
－
Δ
H
f
)
（
11
）
Q
L
＝
M
L
C
p
Δ
T
（
12
）
式（
10
）中，总能耗
Q
T
由
Q
L
与
Q
v
组成，分别代表物料升温需要的能量以及蒸发需要的能量。
M
v
代表水挥发速率，
Δ
H
v
代表水蒸气在此压强、温度下的焓，
Δ
H
f
代表饱和水蒸气在此压强、温度下的焓。物料守恒定律，可以算得浓缩物料的水挥发速率，见式（
13
）。
M
v
＝
M
L
(1
－
x
)
（
13
）
式（
13
）中，
M
L
代表物料进入单效浓缩器中的速率，
x
代表物料固含率。其中，物料固含率会随着物料浓缩过程变化。此时，根据测得物料的白利度读数，可用“
3.1
”项得到的白利度与固含率关系模型得到某一时刻的固含率，从而算出
M
v
。
另外，式（
12
）中
C
p
（物料的比热容）会随着中药物料的浓缩过程而变化
[3]
。因此
Q
L
也会随浓缩过程中，中药物料浓度升高而变化。根据“
3.2
”项中白利度与热容率的关系，可计算出热法单效浓缩的能耗。
3.3.2
膜蒸馏能耗计算相关性
同理，膜蒸馏浓
缩中药物料的能耗也与中药浓缩物料特征有关系。
膜蒸馏
MD
理论最低所需能耗（
Q
f
）计算公式如下
[21]
。
Q
f
＝
A
f
h
f
(
T
f,b
－
T
f,m
)
（
14
）
式（
14
）中，
A
f
代表膜接触物料端面积，
h
f
代表中药物料端热对流导热系数，
T
f,b
代表中药物料端温度以及
T
f,m
代表中药物料端膜表面温度。其中，膜蒸馏物料端膜表面的温度
T
f,m
不能被直接测出来，因此，需要根据物料衡算和传热传质模型的推导。膜表面温度
T
f,m
可用以下公式进行计
算
[22]
。
T
f,m
＝
[
k
m
(
T
p,b
