乔璐，符佳豪，侯益民，李雪菊，赵博宇，董诚明，*，张园园，吕芳

（1.河南中医药大学药学院，河南郑州450046；2.江苏省中国科学院植物研究所，江苏南京210014）

香附为莎草科植物莎草（Cyperus rotundus L.）的干燥根茎[1]。常用于治疗肝郁气滞、月经不调、脾胃气滞等[2]，同时，香附又有妇科之要药的美名[3]。在临床应用非常广泛，常被制成“香砂养胃片”、“四制香附丸”、“调经益灵片”、“七制香附丸”等多种中成药，香附所治疗的疾病范围非常广泛，因此，它在中医临床理气类中药材中也经常被应用[4]。在临床应用中，香附主要行使解热、镇痛、镇静、解痉、抗菌、保护胃黏膜等作用[5]。正是由于临床应用的需要，目前市场中的香附饮片主要就是生香附和醋炙香附[6-7]。

近红外光谱是指在红外光谱范围中短波长的一类电磁波[8]，依据分子振动吸收的数据不同从而实现物质检测分析的目的[9-10]。在近红外检测范围内，药材中的物质成分都能得到大量有效的信息反馈，从整体上反映中药的化学成分的状况，对差异细微的红外光谱进行鉴别，同时傅里叶近红外光谱（Fourier transform infrared spectroscopy，FTIR）分以方便快捷，对样品无破坏性，对环境无污染，样品损耗少等优点[11-12]。近年来红外光谱分析技术在中药材的应用较为广泛[13-14]。但对于基础指标水分、灰分、浸出物的定量分析较少，本文以傅里叶变换红外光谱结合化学计量法对香附水分、灰分和浸出物进行快速测定，为药材香附质量标准的制定提供科学依据。

1 材料与方法

1.1 材料与试剂

1.1.1 试验材料

分别采购了15批共60个香附样品，经河南中医药大学药学院董诚明教授鉴定为莎草科植物莎草Cyperus rotundus L.的干燥根茎，详细信息见表1。

1.1.2 试剂

乙醇（分析纯）：天津致远化学试剂有限公司；甲苯（分析纯）：山东烟台双双化工有限公司。

表1 香附产地Table 1 Cyperus rotundus place of origin

1.2 仪器与设备

ANTARIS II型傅立叶变换近红外光谱仪：Thermo SCIENTIFIC美国电热公司；BX-2.5-10型数显箱式电炉：上海柏欣仪器设备厂；FA2004N型电子天平：上海民桥精密科学仪器有限公司。

1.3 试验方法

1.3.1 香附的处理

60个香附样品均经过打粉处理，过5号药典筛。

1.3.2 香附的水分、灰分和浸出物测定

依据2015版中国药典规定的甲苯法测定香附样品的水分含量（通则0832第四法），不超过13%。总灰分测定法测定香附样品的灰分含量（通则2302），不超过4%。醇溶性浸出物测定法（通则2201）项下的热浸法测定，含量不低于15%，用稀乙醇作溶剂。

1.4 光谱的采集

本试验运用ANTARIS II FT-NIR Analyzer对60个香附样品进行光谱采集，结果见图1。

图1 香附样品的近红外光谱Fig.1 NIR spectra of Cyperus rotundus

将香附样品粉末均匀放置于5 cm直径采样杯中。方式为旋转采样，光谱范围为12 000 cm-1～4 000 cm-1。背景光谱扫描32次，每个样品扫描64次，扫描时间平均为60 s。试验所得样品光谱为直接去除背景的光谱。

1.5 模型的建立

1.5.1 光谱的预处理

对香附的原始光谱进行预处理，包括光谱的算法、导数、降噪、基线校正等，适宜的预处理方式会对于最终模型的精确度产生重要影响。本试验中采用二阶导数光谱（second derivative），并且选择导数滤波（norris derivative filter，ND）、卷积法（savitzky-golay，SG）和无平滑（no smoothing，NS）方法进行平滑处理。选择的光谱范围是 12 000 cm-1～4 000 cm-1。

图1中60份香附样品的近红外光谱整体趋势相同，全部样品光谱在10 000 cm-1～12 000 cm-1光谱范围中波动较大，应该是噪声所致。60个样品光谱在10 500 cm-1～8 800 cm-1范围内趋势平缓；在 8300、6800、5 600、5 200、4 700、4 400、4 300 cm-1处可见明显波峰，在 5 400、5 000、4500、4200 cm-1处可见明显波谷。

1.5.2 模型的验证与校正

本试验借助于TQ Analyst软件来进行香附快速测定模型的建立，采用定量校正算法，主要使用偏最小二乘法（partial least squares，PLS）与逐步多元线性回归法（sparse multinomial logistic regression，SMLR）。光程选择时采用多元信号修正（multiplicative signal correction，MSC）与内参考标准法（internal reference，IR），以此来减小误差，设置目标组分含量的上下限时，将下限值设置为香附样品目标组分含量最低值的95%，上限值设置为香附样品目标组分含量最高值的105%。光谱范围的选择则依据TQ Analyst软件中的suggest向导进行。

在进行光谱预处理时，本试验分别采用IR、MSC方法并以内部交叉验证模型相关系数（R2）为指标，其中模型的相关系数（R2）越接近1，说明试验所得样品的化学值与近红外模型的预测值越接近，即相关性越好，而均方根误差校正集预测均方差（root mean square error of calibration，RMSEC）与预测集预测均方差（root mean square error of prediction，RMSEP）的值越小，则说明所建模型的预测性能越好[15-16]。

2 结果与分析

2.1 香附的水分、灰分、浸出物含量

对香附的水分、灰分、浸出物含量测定结果见表2。

表2 香附的各成分含量Table 2 Content of constituents of Cyperus rotundus

由表2可知，香附水分含量6.27%～13.7%，灰分含量1.69%～4.04%，浸出物含量15.04%～25.77%，所得数据均符合2015版药典规定。

2.2 香附水分、灰分、浸出物定量模型的建立

对光谱模型进行拟合，筛选出最佳模型见表3。

表3 光谱处理方法分析结果Table 3 Analysis results of spectral treatment method

续表3 光谱处理方法分析结果Continue table 3 Analysis results of spectral treatment method

对光谱模型进行拟合，筛选出最佳模型。其中，香附的水分模型以SMLR+MSC+2nd Der+ND模型最好，香附的灰分模型以PLS+MSC+2nd Der+NS模型最好，香附的浸出物模型以SMLR+IR+2nd Der+ND模型最好。

2.3 香附水分灰分浸出物模型

分别建立的香附水分、灰分和浸出物的预测模型，结果见表4。

随机选取建模样本为52个，验证样本8个，光谱范围为12 000 cm-1～4 000 cm-1。香附水分的预测模型精度R2=0.902 7，RMSEC=0.774，预测均方差RMSEP=0.521；香附灰分所得模型精度R2=0.988 4，RMSEC=0.061 9，RMSEP=0.058 5；香附浸出物模型中，精度R2=0.890 3，RMSEC=1.25，RMSEP=2.33。所建模型精度高，均方根误差小。

2.4 香附水分、灰分、浸出物模型对未知含量样品的预测

对香附水分、灰分、浸出物模型进行验证，结果见表5。

表4 香附模型Table 4 Cyperus rotundus model

表5 香附水分、灰分、浸出物含量的预测Table 5 Prediction of moisture，ash and extractum content of Cyperus rotundus

从表5中的数据可以看出，在随机挑选出的8个香附样品中，所建立的香附灰分快速测定模型在绝对偏差上是明显要低于香附水分和浸出物快速测定模型，而本身在所得3种模型的精度上就是香附灰分快速测定模型的精度最高，相关系数R2高达0.988 4，也就进一步证明了试验所得模型可以应用于香附水分、灰分、浸出物含量的快速测定，且测定精确度与相关系数R2呈正相关。将未参与建模的8个样本带入模型，回执预测值和实测值的散点图见图2～图4。检验样本集中在趋势线y=x附近，表明模型精度高，稳定性强。

图2 香附水分预测值与实测值Fig.2 Predicted and measured values of Cyperus rotundus moisture content

图3 香附灰分预测值与实测值Fig.3 Predicted and measured values of Cyperus rotundus ash content

图4 香附浸出物预测值与实测值Fig.4 Predicted and measured values of Cyperus rotundus extractum content

3 结论与讨论

本研究利用傅里叶近红外光谱（NIRS）分析技术建立了香附水分、灰分、浸出物的快速预测模型，并应用未参与数据验证，实现了香附样品中相关指标的快速精准的测定，其中香附水分模型中，光谱预处理选择多元信号修正（MSC）+二阶导数光谱+ND，逐步多元线性回归法建模，光谱范围12 000 cm-1～4 000 cm-1。模型 R2为 0.902 7，RMSEC 为 0.774，RMSEP 为 0.521。香附灰分以偏最小二乘法建模，光谱处理包括MSC+2nd Der+NS，模型 R2为 0.988 4，RMSEC 为 0.061 9，RMSEP为0.058 5。香附浸出物模型中，光谱处理包括内参考标准法IR+2nd Der+NS，选择偏最小二乘法建模，模型 R2为 0.890 3，RMSEC为 1.25，RMSEP为2.33。其中香附灰分模型精确度较高，水分和浸出物的模型精确度略低，可能是在香附样品水分、浸出物含量测定时的误差所致。以上所得模型对于香附水分、灰分、浸出物含量的预测均有较高的精确度，可实现近红外光谱分析技术对香附水分、灰分、浸出物含量的快速测定，这为香附质量标准的制定提供理论依据。