全 文 ：中草药 Chinese Traditional and Herbal Drugs 第 42 卷 第 10 期 2011 年 10 月

• 2114 •
红外光谱法结合系统聚类和 SIMCA 模
1, 2 2，

2. 广西大学农
大
目 方法 在
SIM
活的 0%，两者盲
论 应用
025
ntification of Notopterygii Rhiz
SHEN Liang 2, ZHOU
y of
, C
ine Sciences, C 41, China
mental Science and Engineering, Sichuan Uni
Abst izom m infrared spectrometry (FTIR).
0 to 400 c
e se
ce ysis, the
nd 97%, r
testing with the blind sample which the authors picked out from the chosen samples, the accuracy rate reached over 90%. The
established model could be used for samp d FTIR combined with SIMCA provide
an effect s a new
techn ting y
Key wo pterygii Rhizoma et Ra ed identifica TIR; hierar cluster Indepe Model Class
Analog )

羌 terygii Rh et R 传统常用大
宗中药 有散寒、 除湿 痛的功效，
主治风 、骨节酸 痈疽疮 症[1]，入方
成药达 全国关 700 家。
《中国药典》 10 年版规
植物羌活 erygium incisum ex H hang
和宽叶羌 franc 的干燥根茎及
根，为生长于青藏高原东南部高海拔山地的中国特
有属植物 售的羌 材在商 格和利用价
值上历来 羌活优 叶羌活 十余年高强

式识别法快速鉴别羌活种子
孙 辉沈 亮 ，蒋舜媛3*，黄荣韶 ，周2 毅 ，李良波3 4，甘凤琼2
1. 中国医学科学院 北京协和医学院药用植物研究所，北京
学院，广西 南宁
100193
530005
3. 四川省中医药科学院，四川 成都 610041
4. 四川 学 环境科学与工程系，四川 成都 610065
摘 要： 的 采用傅里叶变换红外光谱法鉴别羌活种子。 4 000～400 cm−1内测定羌活种子光谱吸收峰，应用系
运用系统聚类法，当聚类距离为 15 时，两种羌活可明显分
种子识别率分别达 96%和 97%，拒绝率均为 10
统聚类和 CA模式识别法鉴别羌活和宽叶羌活的种子。结果
为 2 类；运用SIMCA模式识别法，所建模型对羌活和宽叶羌
样检测准确率也在 90%以上，所建模型可用于样品检测。结
简单、快速、准确地鉴别羌活种子，该方法为鉴别羌活种子的来源及真伪提供了一种新的技术手段。
关键词：羌活；种子鉴定；红外光谱；系统聚类；SIMCA
中图分类号：R282.5 文献标志码：A 文章编号：
红外光谱法结合系统聚类和SIMCA模式识别法，可以

3 - 2670(2011)10 - 2114 - 05
oma et Radix by SIMCA based on FTIR
Yi
Rapid seed ide
1, 2, JIANG Shun-yuan3, HUANG Rong-shao
1. Institute of Medicinal Plant Development, Chinese Academ
100193, China
2. College of Agriculture, Guangxi University, Nanning 530005
3. Sichuan Academy of Traditional Chinese Medic
3, LI Liang-bo2, SUN Hui4, GAN Feng-qiong2
Medical Sciences and Peking Union Medical College, Beijing
hina
hengdu 6100
4. Department of Environ versity, Chengdu 610065, China
a et Radix by Fourier transforract: Objective To identify the seeds of Notopterygii Rh
Methods Based on the fingerprint infrared spectrum from 4 00
Class Analogy (SIMCA) analyses were developed to identify th
seeds could be divided into two categories with 15 Euclidean distan
recognition rate of N. incisum and N. franchetii reached to 96% a
m−1, hierarchical cluster and Soft Independent Modeling of
eds of Notopterygium incisum and N. franchetii. Results The
by using hierarchical cluster analysis. Using SIMCA anal
espectively, and the rejection rate reached to 100%. When
le testing. Conclusion Both clustering analysis an
ive way to evaluate the seeds of Notopterygii Rhizoma et Radix rapidly and accurately. Therefore, this method provide
ique for evalua the source and authenticity of the seeds of Notopter gii Rhizoma et Radix.
rds: Noto dix; se tion; F chical ; Soft ndent ing of
y (SIMCA
活Notop izoma adix为
材，具 祛风、 、止
寒感冒 痛、 毒等
280 余种， 联的制药企业超过
20 定其药材为伞形科羌活属
Notopt Ting . T. C
活N. hetii H. de Boisseu
[2]。市 活药 品规
都是 于宽 ，近

收 1-01-15
“重大新药创制” 大专项（民口） “中药材种子种苗和种植（养 准平台 （2009Z 002）；
099200
（1986 ，山东 在读博士，主 药用植物栽培 保护研
6678824 ail: she 8@126.com
媛 Tel: 72450 E il: jsy007@vip.sin
稿日期：201
基金项目：国家科技部
广西科学研究与技术开发计划项目（桂科攻
科技重 课题；
3A-20）资助
殖）标 ”课题 X09308-
作者简介：沈 亮 —），男 枣庄人， 要从事 育种及 究。
Tel: 13
*通讯作者 蒋舜
55 E-m
189820
nliang0
-ma a.com
中草药 Chinese Traditional and Herbal Drugs 第 42 卷 第 10 期 2011 年 10 月

• 2115 •
鉴别
羌活种子质量控制和人工种植的一个关
合系统聚类分析及SIMCA模式识
有快
新的分
材
1.1
羌活和宽叶羌活种子共 17 份，由四川省中医药
号为 6 份不同来源的宽叶羌活种子
9～1 号为 8份不同来源的羌活种子（定为B种子），
7、8
器内
津光
1.2
mg
明锭片。
备 个溴化钾压片，每个待测样品随机扫描 次，
收集所得光谱图。
原始光谱数据经过 Omnic 8.0 软件进
正，13 点平滑，平均
Unscrambler 10.0
理，以消除样本称量不均带来的光谱峰吸光度的差
异，保证谱图规格
分析使用 SPSS 11.5
Tab
度的采挖已致其资源严重濒危。羌活种子发芽育苗
技术的突破为实现羌活的人工栽培和可持续利用
提供了技术支撑，但两个种的种子在形态上难以分
辨，种子前处理及育苗时间长达 9 个月以上[3]，使
形态法、幼苗法等常用的种子鉴定方法难以实施，
采用电泳法等其他生物技术手段也因实验材料缺
乏稳定性和准确性，羌活种子真伪的快速准确
成为困扰
键问题。
红外光谱结
别法在药材鉴定方面已有相关报道[4]，该方法具
速、准确、经济等特点。本实验尝试运用该
项技术对羌活种子进行鉴别研究，为快速、准确
同来源的羌活属植物种子探索一种鉴定不
析方法。
1 料与方法
材料
科学院蒋舜媛博士提供并鉴定（表 1）。其中 1～6
（定为 A 种子），
6
和 17 号为待测盲样。将各样品清洗后，放
入烘箱，40 ℃低温烘干至恒质量，去种皮，粉碎，
过 200 目筛，将样品包装好后放于恒温烘箱或干燥
保存，待测定。
Nicolet 5700傅氏转换红外线光谱分析仪，检
测器为美国尼高力公司生产；FW—4A型压片机（天
学仪器厂）；分析天平XSI05DU（梅特勒–托
利多仪器公司）；溴化钾（天津光复精细化工研究所，
色谱纯）。
方法
精确称取 1.5 mg样品，将其与完全干燥的 200.0
溴化钾混合研磨均匀，压制成厚度约 1 mm的透
扫描范围 4 000～400 cm−1 分辨率，光谱
4 cm−1，累积次数 32 次。每个产地的羌活种子各制
20 3
行基线校
图谱后，将光谱数据导入
软件中，对谱图进行标准归一化处
的统一性和可比性。系统聚类
软件，聚类采用内平方距离法
表 1 两种羌活种子样品的产地及采收时间
le 1 Habitat and collected time of N.franchetii
and N. incisum
序号 来 源 采收时间 拉丁名
四川小金县家种 2010-09 N. franchetii 1
2 四川壤塘县家种 2010-09 N. franchetii
四川小金县两河乡野生移栽 2010-09 N. franchetii
四川甘孜县野生移
3
4 栽 2010-09 N. franchetii
6 源县家种 2010-08 i
8
9 2010-09 N. incisum
11
13
14 2010-09 N. incisum
15 四川壤塘县家种 1 2010-09 N. incisum
16 四川壤塘县家种 2 2010-09 N. incisum
2010-09
5 四川壤塘县阿斯玛乡野生 2010-09 N. franchetii
甘肃渭 N. francheti
7 待测盲样 1 2010-09
待测盲样 2 2010-09
四川小金县家种
10 四川壤塘县家种 2010-09 N. incisum
四川甘孜县泥柯乡野生移栽 2010-09 N. incisum
12 四川小金县梦笔山野生 2010-09 N. incisum
四川壤塘县阿斯玛乡野生 2010-09 N. incisum
四川壤塘县野生
17 待测盲样 3

[5]
化
归属可知，3 305 cm−1波数附近是多糖、苷类化合物
伸缩振动峰；2 925、2 854 cm−1附近为
白质和油脂类甲基、亚甲基C-H伸缩振动峰；1 741
缩振动
峰；这部分主要体现了羌活种子中蛋白质、挥发油、
cm−1的吸收峰可以推断其种子中含有大量的油脂类
60 cm−1附近为C＝O或C＝C双键的伸缩振
动区，可能为苯环及双键的骨架振动，此外还包括
I带；而在 1 616
-NH2的剪式振动酰胺II
近出现了芳环
活种子
质的光谱
（Word法），采用欧氏距离法 计算距离，SIMCA建
模预测分析在Unscrambler 10.0 软件中进行。
2 结果与分析
2.1 羌活与宽叶羌活种子的红外光谱分析比较
红外光谱的特点是能够解析出样品化学成分的
官能团，确定其基本骨架，从而判断出样品所含
学成分[6]。图l为羌活和宽叶羌活种子的原始红外光
谱图，参考相关文献[7-9]数据对其光谱峰进行指认和
中的羟基O-H
蛋
cm−1附近是羧酸类或酯类物质中C＝O的伸
香豆素等脂类物质的光谱表征，而从 2 852、1 741
物质。1 6
由蛋白质C＝O伸缩振动产生的酰胺
cm−1波数附近出现了氨基酸
带特征吸收峰；1 518、1 456 cm−1附
的骨架振动吸收峰；这部分主要体现了羌
含有氨基酸、香豆素类物
中草药 Chinese Traditional and Herbal Drugs 第 42 卷 第 10 期 2011 年 10 月

• 2116 •
是C-H的弯曲振 分主要体现
有多糖类物质的特征吸收峰。虽然两种羌活种子的
形及峰强的差异还是为鉴别两种羌活种子提供了
理
收
活种
子
吸
cm−1附近也有不同程度的吸收峰，通过原始光谱可
和相似程度。待测盲样 7、8 号与宽
叶羌
2.3 种子的 SIMCA 模式识别分析
为了更加
系统
[10]对

ectra recorded of N.
sum seeds from 17 habitats
以宽叶羌活种子的样品（1～6 号）为纵坐标，羌活
的样品（9～15 号）为横坐标，经过处理得到
子的SIMCA聚类分析图。由图 3
果比较理想，两者


一步的模型诊断，即求出每种羌活种子的识别率和
拒绝率。所谓识别率就是考察某类样品有多少落在
模型的范围中，而拒绝率就是考察某类样品模
样品的拒绝程度，即其他
样品是否都落在该类模型的范围外。当这两个值都

图 1 宽叶羌活（A）和羌活（B）种子的红外原始光谱图
Fig. 1 FTIR spectra of N. franchetii (A) and N. incisum (B) seeds
表征。1 157、1 070 cm−1重叠而成的宽峰，归属为
多糖、苷类等碳水化合物的C-O弯曲振动，而 857 cm−1
动峰，这部 羌活种 franchetii and N. inci子含
红外光谱图较为相似，但是某些吸收峰的峰位、峰
种子的智能识别。
2.3.1 羌活和宽叶羌活种子模型的建立 本研究经
论依据。如羌活种子在 1 741、1 616 cm−1处的吸
过数据预处理后，选择羌活种子光谱信息丰富的
−1
峰峰强明显弱于宽叶羌活种子；在 1 660 cm−1处羌
1 800～400 cm 波段进行SIMCA模式识别分析，
子的酰胺I带吸收峰又明显强于宽叶羌活种
−1
种子
；此外二者在 1 518、1 456 cm 附近的芳环振动 羌活和宽叶羌活种
收也存在明显的差异，在 1 279、1 261、857、673 可见，羌活和宽叶羌活的聚类结
以有效地对两种羌活种子进行鉴别及化学成分分
析。
2.2 羌活和宽叶羌活种子的系统聚类分析
对羌活样品进行归一化处理后，将 1 800～400
cm−1波段的吸光度进行系统聚类分析，结果如图 2
所示。当聚类距离为 15 时，羌活样品明显的分为 2
类，而不出现交叉。其中 1～8 号聚为一类，9～17
号聚为另一类，它们分别代表宽叶羌活和羌活种子
各自聚到一起，没有出现重叠。
的光谱信息，系统聚类方法可以直观地反映出样品
间的亲缘关系之
活类群聚到一起，而待测盲样 17 与羌活类群聚
到一起，这与本课题组采用其他常规方法及分子鉴
定方法所得结果一致[10]。另外系统聚类法还可以推
断出样品亲缘关系的远近，由图 2 可知两种羌活的
差异较大，当聚类距离接近 25 时，两者才出现较为
图 3 宽叶羌活（A）和羌活（B）的 SIMCA 聚类分析图
Fig. 3 SIMCA classification of N. franchetii (A)
and N. incisum (B) seeds
2.3.2 模型的诊断 模型建好后，需要对其进行进
相似的光谱信息而聚为一类。
羌活和宽叶羌活
直观地区分两种羌活种子并验证以上 该类
聚类分析结果，本研究选用SIMCA模式识别法 型对于不属于该类的其他
羌活种子做进一步研究，以实现对两种羌活
图 2 17 个不同产地的羌活种子红外光谱系统聚类图
Fig. 2 Hierarchical cluster of FTIR sp
00 3500 3000 2500 2000 1500 1000 500
A
B
4 000 3 000 2 000 1 000
波数 −1 / cm
0 5 10 15 20 25
15
1
12
16
11
2
3
5
4
7
14
9
13
10
6
7
1
8


0.09
0.08
0.07
0.06
0.05
0.04
0.03
0.02
0.01
0.00
0.01 0.02 0.03 0.04 0.05

A
B
样
品
距
PC
A
模
型
的
距
离

样品距 PCA 模型的距离
中草药 Chinese Traditional and Herbal Drugs 第 42 卷 第 10 期 2011 年 10 月

• 2117 •
为 2 个产地（为待测种子 7、8 号）的
20 个
活种子的鉴别研究。
原始光谱图分析发现，羌活和宽叶羌活
种子
，二者的聚类判别准确率分
别为
外观相似的羌活和宽叶羌活种子的快速
期较长的种子，则需花费大量
的人
的种子 但这
些方法一方面对仪器及人员的技能要求较高，花费
较高[11]， 验材料的要求对某些种类来
说难以满足。
仪收集谱图进行数据处理即可得出鉴定结果，具有
方便、快速、准确、灵敏、经济等特点。其不足之
处在于大型红外光谱仪价格较高，同时由于实验用
材较少，检测灵敏性较高，取样必须有较好的代表
性。目前已有运用红外光谱法对砂仁种子及紫苏子
的真伪进行快速鉴别研究的相关报道 6]，而红外
光 合化 学在种子鉴
如系统聚类法、SIMCA模式识别法具有直观、
可验证性等优点，可以快速客观的区分形态和解剖
手 较 真伪 ，
在种子鉴定领域里，红外光谱鉴定法与其他鉴定方
法相比，具有操作简单、鉴定快速、准确，花费较
少等优点而具有较大的应用价值。随着科学技术的
发展，红外光谱鉴定法可以在种子鉴定领域发挥更
大的作用。
中国特有属—羌活属的
[J]. 云南植物研究 , 1996, 18(4):
[3]
[J]
rimination of the
为 100%时，表明这两种样品之间差异很大，可以
完全分开。经过多次重复试验及数据处理，所建模
型对羌活和宽叶羌活种子的识别率分别达到了 96%
和 97%，拒绝率均为 100%，说明所建模型较为成
功，准确性较高。
2.3.3 模型的检验 为检验所建模型是否具有较全
面的样品代表性，本研究采用测试集样本对模型进
行了测试。测试集是由 40 个非建模集样品构成，其
中宽叶羌活
种子样品，羌活为 2 个产地（16 号和待测种子
17 号）的 20 个种子样品。经过预测得知，宽叶羌
活种子判别准确率为 90%，而羌活种子的判别准确
率为 95%，测试正确率较高，说明建模样品数量合
理，所建模型具有一定的代表性，可以用于两种羌
活及未知羌
3 讨论
通过对
的化学成分存在一定的差异，但尚不能从原始
图谱上进行直接区分。进一步使用系统聚类法进行
分析，两种羌活种子样品各自聚到一起，而不出现
交叉，当聚类距离接近 25 时，二者聚为一类，较好
地实现了不同种羌活种子的区分。采用 SIMCA 模
式识别法对羌活种子和宽叶羌活种子的识别率分别
达到了 96%和 97%，拒绝率均为 100%。用待测盲
样对所建模型检验得出
95%和 90%，可见该方法能够用于对未知样品
的检测。本研究表明，借助羌活种子的红外光谱信
息，从整体上可以解析出其种子内含有的化学成分，
在此基础上采用系统聚类和 SIMCA 模式识别法，
可以实现对
鉴别。
目前有关种子鉴定的方法较多，对于区别性特
征很显著的种子，可采用传统的外观形态鉴定法，
外观相似度较高的种子则需要依赖长期积累的鉴别
经验为基础，且难保证高的准确性[11]。幼苗观察法
可确保鉴定结果准确，对于发芽周期短的植物种子
适用[11]，但对发芽周
力、物力、财力及时间，尤其对从种子前处理
到发芽出苗，并长到形态可鉴别的幼苗所需时间长
达 9 个月以上的羌活种子来说，幼苗法可作为复核
的技术手段，但对种子质量控制和人工种植生产则
缺乏实用价值。近年在种子鉴定方面发展起来的电
泳法、ISSR标记、rDNA标记等方法，对某些类群
能实现快速鉴定，准确性也较高[12-14]，
另一方面对实
红外光谱鉴定法只需少量粉碎样品，使用光谱
[15-1
谱法结 学计量 定领域更具优势，
准确、
段 难鉴别的不同来源及 的种子[17-19]。可见

参考文献
[1] 中国药典 [S]. 一部. 2010.
, 王萍莉, 等. [2] 王幼平, 溥发鼎
系统分类研究
424-430.
史 静, 蒋舜媛, 马小军, 等. 羌活种子发芽及实生苗
生长发育的研究 [J]. 中国中药杂志 2007, 32(18):
1841-1844.
[4] 徐永群, 黄 昊, 周 群, 等. 红外指纹图谱和聚类分
析法在赤芍产域分类鉴别中的应用 [J]. 分析化学 ,
2003, 31(1): 5-9.
[5] 孔德鑫, 黄庶识, 黄荣韶, 等. 基于双指标分析法和聚
类分析法的鸡骨草红外指纹图谱比较研究 [J]. 光谱学
与光谱分析, 2010, 30(1): 45-49.
[6] 翁诗甫. 傅里叶变换红外光谱分析 [M]. 北京: 化学工
业出版社, 2010.
[7] 陈 君, 孙素琴, 徐 荣, 等. 应用红外光谱法鉴别肉
苁蓉及其混淆品草苁蓉和锁阳 [J]. 光谱学与光谱分
析, 2009, 29(6): 1502-1507.
[8] 孙素琴, 周 群, 秦 竹. 中药二维相关红外光谱鉴定
图集 [M]. 北京: 化学工业出版社, 2003.
[9] 黄冬兰, 陈小康, 徐永群, 等. 纹党参与白条党参红外
光谱的 SIMCA 聚类鉴别方法研究 . 分析测试学报,
2009, 28(12): 1440-1443.
[10] Xu K J, Jiang S Y, Zhou Y, et al. Disc
中草药 Chinese Traditional and Herbal Drugs 第 42 卷 第 10 期 2011 年 10 月

• 2118 •
seeds of Notopterygium incisum and Notopterygium
franchetii by validated HPLC-DAD-ESI-MS method and
principal component analysis [J]. J Pharm Biomed Anal,
2011, 56: 1089-1093.
[11] 白 雁, 鲍红娟, 王 东, 等. 红外光谱和聚类分析法
在药用菊花产域分类鉴别中的应用 [J]. 中药材, 2006,
29(7): 663-665.
[12] 王 芳. 种子纯度鉴定方法及其评述 [J]. 中国种业,
2008, 10: 62-63.
[13] 高燕会, 李慧慧, 朱玉球, 等. 基于 ISSR 的栝楼 传
[15] 郝剑瑾, 程 舟, 梁洪卉, 等. 基于 rDNA ITS 序列探
讨我国冬虫夏草的遗传分化及分布格局 [J]. 中草药,
2009, 40(1): 112-116.
[16] 程存归, 阮永明, 李冰岚. 傅里叶变换红外光谱法应用
于中药砂仁真伪鉴别的研究 [J]. 光谱学与光谱分析,
2004, 24(11): 1355-1358.
[17] 程存归. FTIR 直接鉴定紫苏子及其伪品的研究 [J]. 光
谱学与光谱分析, 2003, 23(2): 282-284.
[18] 姜科声, 孔黎春, 余 鹏, 等. 不同花生品种的红外光
谱分析 [J]. 光谱实验室, 2008, 25(4): 550-553.
[19] 周 群, 孙素琴, 梁曦云. 枸杞产地的红外指纹图谱与
聚类分析法研究 [J]. 光谱学与光谱分析, 2003, 23(3):
s（
icin
10/R
所
工程
林研
《
药 ，宣扬我国中草药
传统与现代、东方与西方之间架起一座理解和沟
梁，促进中药现代化、国际化。
文摘、信息和国际动态、人物介绍、来信、书评等栏目。
读者对
邀请相关领域的院士和国内外知名专家加盟，组建一支国际化、高水平、精干的编委会队伍（第一届编辑委员会
9 提高期刊的学术质量；坚持按照国际标准
和英语，加强与国际知名科技期刊的交流合作；充分发挥中医
美国科学引文索引（SCI），把 CHM 办成国际知名期刊之一。
发行，直接与编辑部订阅！
阅！网址: www.tiprpress.com
京编辑部
北京市海淀区
mai
Tel: (0
: (0
号：
遗
多样性分析 [J]. 中草药, 2011, 42(2): 363-366.
[14] 蒲汉丽, 李乃坚, 黄爱兴, 等. RAPD 技术鉴定番茄一代
杂种粤星种子纯度的研究 [J]. 种子, 2000, 111(54): 14-15.


欢迎订阅
Chinese Herbal Medicine
我国第一份中药专业的英文期刊——Chinese Herbal Med
2009 年 10 月正式创刊，国内统一连续出版号为：CN12—14
CHM 由天津药物研究院和中国医学科学院药用植物研究
科学院药用植物研究所名誉所长肖培根教授担任主编；中国
长汤立达研究员，中国医学科学院药用植物研究所所长陈士
中草药》杂志执行主编陈常青研究员担任编辑部主任。
办刊宗旨 以高起点、国际化为特点，继承和发扬祖国医
的传统特色，加强与世界各国在传统药物研究的经验交流，在中医和西医、
通的桥
509-511.
CHM, 中草药英文版）
es（CHM，中草药英文版）经国家新闻出版总署批准，已于
。
主办，天津中草药杂志社出版。中国工程院院士、中国医学
院院士、天津药物研究院刘昌孝研究员，天津药物研究院院
究员共同担任副主编；天津药物研究院医药信息中心主任、
学遗产，报道和反映中草药研究最新进展
主要栏目 综述与述评、论著、快报、简报、
象 国内外从事中医药研究、管理、监督、检验和临床的专业技术人员。
CHM
由 4 位专家组成，其中院士 10 名，国际编委 19 名）。吸引国内外高质量的稿件，
编排，加强刊物规范化和标准化，充分利用计算机、网络技术
药特色，争取在较短时间内进入国际最著名的检索系统——
欢迎广大作者踊跃投稿！ 欢迎广大读者积极订阅！自办
本刊已正式开通网上在线投稿系统。欢迎投稿、欢迎订
Chinese Herbal Medicines（CHM, 中草
天津编辑部 北
地址：天津市南开区鞍山西道 308 号 地址：
邮编：300193 邮编：
E-mail: chm@tiprpress.com E-
Tel: (022)27474913
药英文版）编辑部
马连洼北路 151 号
100193
l: bjchm@tiprpress.com
10)57833035
10)57833035
44114010010081504 户名：天津中草药杂志社
Fax: (022)23006821 Fax
开户银行：兴业银行天津南开支行 账