基于中药水提液浓缩过程溶液环境特征参数相关性的瞬时能耗计算方法 探索——以玉屏风散水提液为例

摘  要：目的  寻找中药水提液浓缩过程瞬时能耗计算简便方法。方法  以玉屏风散为模型药物，测定其水提液浓缩过程中多种溶液环境特征参数的动态变化数据，构建中药溶液环境特征参数白利度（简便、易测、准确）与浓缩能耗计算需要的相关参数的相关性模型。结果  白利度与固含率、密度、比热容、导热率等具有优越的线性相关，相关系数（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₁。将同一批次取得的物料放入烘干的玻璃培养皿，并准确称定质量，记为m₂。将已称定质量的物料与玻璃培养皿放入电热鼓风干燥机中，设置50 ℃干燥4 h以上确保水份完全蒸发。待干燥的物料与玻璃培养皿放置室温后，准确称定质量，记为m₃。此样品的固含率可用公式（1）计算。

固含率＝(m₃－m₂)/(m₂－m₁)                （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₃－T₁)                          （4）

量取物料被加热前的温度（T₁），被加热一段时间（t）后量取其即时温度记录为T₃；水被加热前的温度为T₂，同理水被加热一段时间t后量取其即时温度记为T₄，则有：

Q_pw＝C_pw(T₃－T₁)                       （5）

联立（3）、（4）和（5）得到公式（6）。

C_p＝C_pw(T₄－T₂)/(T₃－T₁)                 （6）

为控制水和物料都得到相同的热量，实验设计水和物料均被电陶炉（选择70 W）加热240 s（t＝240 s）。

另外，根据传统传热学中导热率的定义，一个体系的导热率（k）如下^[18]。

k＝QΔx/AdT                             （7）

同理可推得，同一组实验已经称得的相同质量的水和物料都得到相同的输入热量（Q）、加热接触面积（A）和热量传递距离（Δx），物料的k_f可用已知的水的导热率（k_w）根据公式（8）计算得到。

k_f＝k_w(T₄－T₂)/(T₃－T₁)                   （8）

k_w选择0.640 6W/(m·K)，水在50 ℃的导热率

2.3.6  渗透压（π_f）  按照YASN OM819型渗透压仪（摩尔浓度渗透压仪，摩尔渗透压仪）操作手册执行，由扬子江药业集团南京海陵中药制药工艺技术国家工程研究中心检测。每次测量样品需10 mL，因此，每组样品收集30mL，共测量3组数据，取平均值，得到渗透浓度b₀（mOsm/kg）。根据公式（9）算得渗透压（摩尔浓度渗透压，摩尔渗透压）（kPa）。

π_f＝b₀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²＝0.998 8，y代表固含率，x代表白利度。因此，在玉屏风散的浓缩过程研究中，可以利用白利度去预测玉屏风散的固含率。除了经典复方玉屏风散在浓缩过程中具有白利度与固含率的上述关系外，此前Criscuoli等^[11]也考察了枣汁的膜蒸馏浓缩过程中白利度与固含率之间的关系，得出类似结果。

需要特别指出的是，中药制药工艺往往涉及澄清、精制工段，淀粉、蛋白质等大分子杂质在浓缩工段之前已被去除。因此中药物料中一些大分子物质可能因为加热等问题引起的析出影响可被忽略。

3.2  浓缩过程中中药物料溶液环境表征参数的相关性研究

3.2.1  白利度与密度、黏度的关系  浓缩过程中同一时间点玉屏风散物料的溶液环境特征参数白利度与密度、黏度关系的实测数据分别如表2、3所示。

根据密度的基本定义，物质的密度可以用质量除以体积算得。中药水提液物料的浓缩过程，可认为是一个去除水，留下固体溶质的过程。一般认为固体溶质的密度比水密度大，因此，除水过程中物料的密度应呈上升趋势。玉屏风散浓缩过程中，物料的白利度与密度呈现优越的线性相关，回归方程为y＝4.266 1 x＋1 072.3，R²＝0.993 2，y代表密度，x代表白利度。

黏度是流体的一种特性，客观地表达了流体对流动的阻力，在中药水提液浓缩体系中则反映了水分子在运动过程中的受阻程度。白利度与黏度是2个不同的物理概念，没有直接的关联。本实验所得到的白利度与黏度的关系方程为y＝3.361 9 e^15.593 x，R²＝0.933 1，y代表黏度，x代表白利度。根据二者的关系方程，虽不能由白利度准确计算物料体系的黏度，但是能粗略地估算玉屏风散物料的黏度。中药液体物料黏度的数值主要是用来计算流体体系的雷诺数，以预测流体的分类（湍流、层流），为膜蒸馏浓缩中药物料的即时能耗动态计算提供一个重要的选取依据。

3.2.2  白利度与比热容、导热率的关系  从物质比热容定义及公式（2）可知，某体系的比热容是指在该体系中单位质量的物质改变单位温度（1 K）时需要吸收或释放的热量。而在本研究中的比热容，特指定压比热容（在同一压强下）。

一般认为一个复杂体系中，比热容与混合物质质量呈一定比例。因此在中药溶液环境的复杂体系中，水的去除应该使物料比热容的大小有所改变。浓缩过程中同一时间点的溶液环境特征参数白利度与比热容、导热率的实测数据分别见表4、5。白利度与玉屏风散溶液环境中比热容的关系方程为y＝0.044 x＋4.039 7，R²＝0.987 9，y代表比热容，x代表白利度；在实验初始白利度为2.6°Bx时，溶液的比热容与水的比热容非常接近［4.18 kJ/(kg∙K)］，而比热容随着水分的减少而增加。比热容是能耗测算公式的关键参数，但因为中药水体液物料体系是一种高度多维的复杂系统，其比热容难于测定。因此利用白利度去预测不同浓度溶液下的比热容，对浓缩能耗计算具有填补空白的突破性意义。

导热率（又称导热系数）是一个形容系统热传导能力的物理量，由其计算公式（7）可知，导热率即为在相同的热量传导条件下，温差为1 K的单位厚度的系统在传热为单位面积单位时间条件下的传热量，导热率与系统的成分组成、状态等有关系。根据公式（8）可推知，本研究测试的玉屏风散系统中导热率与系统物料温差成正比。利用白利度与导热率的关系方程（y＝0.006 7x＋0.619 1，R²＝0.9879，y代表导热率，x代表白利度），可以成功地预测出玉屏风散在不同浓度溶液下的导热率。

3.2.3  白利度与渗透压的关系  浓缩过程中同一时间点的溶液环境特征参数白利度与渗透压（摩尔浓度渗透压，摩尔渗透压）实测数据见表6；白利度与渗透压的关系方程为y＝160.54   x－366.66，R²＝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_LC_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_fh_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＋h_fT_f,b/h_p)/δ＋h_fT_f,b－JΔH_latent]/[k_m(1＋h_f/h_p)/δ＋h_f]                                （15）

式（15）中，k_m表示膜的导热系数，δ表示膜的厚度，T_p,b代表纯水透过液端的温度，T_f,b代表中药浓缩物料端的温度，h_f代表物料端的热对流导热系数，h_p代表纯水透过液端的热对流导热系数。J代表膜通量，以及ΔH_latent代表饱和液体水蒸发成水蒸气的汽化潜热。

其中，物料端的膜表面表面温度也需要得到物料的热对流导热系数h_f。然而根据经典流体力学和传热理论，计算热对流的导热系数h_f需要物料的比热容C_p、密度ρ、黏度μ和物料导热率k_f等浓缩过程变化的物理特征，涉及公式（16）～（20）等复杂、繁琐的较多计算过程如下^[22]。

根据（17），首先利用浓缩物料膜端的水力直径d_h、物料端流速ν、物料的黏度μ以及物料的密度ρ算出雷诺数Re。再用（18）与物料的比热容C_p、物料的粘度μ与物料导热率k_f算出普朗特数Pr。根据Re与Pr数，利用（19）或（20）算出新的努塞尔数Nu_f。利用Nu_f代入（16）中，算出新的物料的热对流的导热系数h_f。

h_f＝Nu_fk_f/d_h                            （16）

Re＝d_hνρ/μ                             （17）

Rr＝μC_p/k_f                             （18）

Nu_f＝1.62(RePrd_h/L)^1/3Re＜2300           （19）

Nu_f＝0.023Re^4/5Pr^1/3Re＞2300             （20）

利用“3.2”项得出白利度与其他物理参数如C_p、k_f等关系模型，计算膜蒸馏浓缩中药的能耗。

3.3.3  反渗透能耗相关性  反渗透浓缩中药物料的能耗的计算必须根据物料浓缩过程中渗透压（摩尔浓度渗透压，摩尔渗透压）的变化来计算反渗透技术最低所需能耗（W_min）计算公式如下^[22]。

W_min＝V₀ΔπdR                       （21）

V₀代表中药物料初始体积，R代表浓缩率，Δπ代表物料端与渗透液端渗透压差

Δπ＝π_f－π_p≈π_f                         （22）

π_f为物料渗透压，π_p为渗透液端（清水）渗透压

物料的渗透压（π_f）会随着物料浓缩的浓度上升而上升。根据白利度与渗透压（摩尔浓度渗透压，摩尔渗透压）的关系方程，可利用测得的白利度得到某时反渗透浓缩中药物料的瞬时能耗，从而计算出整个浓缩过程的总能耗。

4 讨论

4.1  有关溶液环境特征参数的讨论

本实验所发现的中药水提液浓缩过程中若干“溶液环境”特征参数的相关性，有关文献亦有过相关报道^[9,17]。这是一种偶然现象，还是具有普遍意义的规律，尚有待进一步深入系统的研究。中药溶液环境的黏度、pH值、离子强度等表征参数，既来源于中药水提液中各种物质的化学组成，又是水提液体系中各种物质不同表现的综合反馈，当然也必定与制药分离过程物料体系中各种物质的热力学和传递性质、多相流、多组分传质等影响分离的因素密切相关^[23]。本团队将在进一步开展系统实验研究的基础上，借助现代分离科学与物理化学等多学科交叉研究的方法，从理论上加以诠释。

4.2  化工过程控制原理是探讨中药制药过程工程原理的重要基石

化工过程控制原理提出：对于复杂的化工过程，不能满足于在现有装置上通过测试获得所需的对象动态特性知识，更重要的是在对象处于设计阶段，就能利用计算方法预估其特性，以改变黑匣无知状态，指导工艺设计中的原理性设计。中药水提液由几十甚至上百种小分子化合物（包括挥发油、氨基酸及生物碱、有机酸、黄酮类、皂苷等化学成分）和生物大分子物质（包括肽、蛋白、糖肽及多糖等）组成，是一种典型的高度复杂的大系统客体，因此，十分有必要引入现代自然科学化工过程控制的方法对其进行研究。节能减排是中药绿色制造的核心目标，智能制造则是中药现代化的重要内容。而传统的浓缩能耗计算公式涉及多达近10个物料理化参数和浓缩过程工艺参数，其中部分参数采集方法繁琐、缺乏精准性，且无法实现在线检测，加上繁琐的推算过程，难以在大生产中推广应用。

本研究将化工过程控制基本原理引入中药生产工艺，结合数学建模的基本方法，将原型的某一部分信息简化、压缩、提炼而构造成的原型替代物，在开展中药制剂、分析化学、现代分离科学及计算机技术等多学科联合攻关的基础上，对中药制药浓缩过程进行深入系统的基础研究。

本研究提出的能耗方法，仅需对某一常年生产的大品种建立以白利度为核心参数，与物料黏度、渗透压（摩尔浓度渗透压，摩尔渗透压）、比热容、密度和导热率等参数的相关性模型，采用有关算法编制相关软件，即可根据实时采集的白利度数据，得出瞬时能耗及同一时刻的物料溶液环境各相关特征参数，为膜浓缩这一绿色制造先进技术在大生产中的推广应用，也为中药生产过程实现在线监测、智能控制提供技术支撑。

利益冲突  所有作者均声明不存在利益冲突

参考文献（略） 

来  源：钟文蔚，黎万钰，丁  菲，李除夕，郑东阳，郭立玮．基于中药水提液浓缩过程溶液环境特征参数相关性的瞬时能耗计算方法探索——以玉屏风散水提液为例  [J]. 中草药, 2021, 52(7): 1937-1944 .