-==228 第31卷第9期 2010年9月 食品研究与开发 天然食品添加剂 黑青稞花色素的提取工艺 赵桃,马林,李嘉佳。杜静君,杨娟。杨凡 (上海工程技术大学化学化工学院,上海201620) 摘要:研究青稞花色素的提取工艺,以黑青稞为原料,考察提取溶剂、酸化剂、提取温度、提取时间、料液比、提取次 数和籽粒粉碎度等7个单因素对其花色素提取率的影响,并采用二次回归正交旋转组合设计法建立和优化青稞色素 提取的乙醇浓度、温度、时间的回归方程。试验结果表明,在试验范围内各因子对青稞色素提取率的作用大小依次为: 乙醇浓度>时间>温度。根据优化的回归方程,青稞色素的最佳提取工艺条件为:乙醇浓度5O%.温度53 ,时间3 h, 此时最大提取率为108.33 P ̄g/g。试验所建立的回归模型在试验范围内能较准确的预测青稞色素的提取率。 关键词:青稞;花色素;提取工艺;二次回归正交旋转组合设计 Extraction Technology of Anthocyanins from Black Hulless Barley ZHAO Tao,MA Lin,LI Jia-jia,DU Jing-jon,YANG Juan,YANG Fan (School ofChemistry and Chemical Engineering,Shanghai University ofEngineering Sciences,Shanghai 201620,China) Abstract:Black hulless barley was used to study the optimal extraction technology of anthocyanins from hulless barley.Effects of extraction solvent,acidifying agent,temperature,time,solid-solvent ratio,frequency and particle size of material on extraction yield of anthocyanins from hulless barley were determined individually. Based on the results,a three—factor quadratic regression rotational combinational design was applied to establish a regression model and 1o estimate the optimal conditions for maximal extraction yield.The results showed that the contributions of the three factors were:ethanol concentration>time>temperature.The optimal extracting 基金项目:上海市教委优秀青年教师基金资助项目(B一8500—08—0¨2);上海工程技术大学科研启动基金(校启08—25);上海工程技术大学大学 生创新活动计划项目(CS0904003) 作者简介:赵桃(1982一),女(汉),讲师,博士,研究方向:植物活性成分。 ●『◆rii●●●◆ii◆i[i●i◆i●● ◆ii◆¨●●i1:●●◆1●ii● ◆ !●i◆●li◆◆i●◆●◆◆◆◆●●●i◆● ●i◆●●◆● l● [8】何方奕,张捷莉,程钰杰.菠萝色素的提取及其稳定性的研究fJ]. 鞍山师范学院学报:综合版,1998,19f4):31—34 【9]司新超,黄和.菠萝皮色素的提取及稳定性研究Ⅲ.广西热带农业, 2007.2:1-4 waste bioprocessed by Rhizopus oligosporus[J1.Food Biochemist,'3r, 2004,28:419-434 [1 6]丁道福.菠萝加工和综合利用[M],北京:中国轻工业出版社,1980: 243 【lO】L H McKEE,T A Latner.Underntilized sources of dietary fiber:A [17】杨礼富,谢贵水.菠萝加工废料——果皮渣的综合利用【JJ.热带农 业科学,2002,22(4):67—7 l review[J].Plant Foods for Human Nutrition,2000,55:285—304 [1 1]T Prakongpan,A Nitithamyong,P Luangpituksa.Extraction and Application of Dietary Fiber and Cellulose from Pineapple Cores[J]. journal offood science,2002,67(4):1308—1313 [18]苑艳辉.菠萝皮的综合利用【J].食品与发酵工业,2005,31(2):145— 147 【l2]严浪,田海军,张树明,等.菠萝半纤维素多糖的提取纯化及免疫 活性研究【J】.食品科学,2008,29(2):35—38 [13】C T Tran,L I Sly,D A Mitchel1.Selection of a strain of Aspe llus for the production of citric acid from pineapple waste in solid—state [19】周立国.食用天然色素及其提取应用[M】.山东:山东科学技术出 版社,1993:1 [2O]全国食品添加剂标准化技术委员会.GB 8818—2008食品添加剂 可可壳色素(s].中国标准出版社,2009 [21]中华人民共和国国家标准.GB/T 12291—90水果、蔬菜汁类胡萝 素全量的测定(S】.中国标准出版社,1990 【22]刘良忠.天然红心蛋中红色类胡萝h素的结构表征及其生物活性 fermentation【J】.World Journal of Microbiology&Biotechnology, 1998。14:399—403 【14】Correia,RobertaTP,McCue,et a1.Production ofphenolic antioxidan- ts by the solid—state bioconvemion of pineapple waste mixed with soy 评价[DI.华中农业大学,2004 [23】毛跟年,许牡丹.功能食品生理特性与检测技术【M】.北京:化学_T 业出版社,2005:557 lfourusingRhi-zopus oligosporu"s啪.Process Biochemistry,2004,39 (12):2167—21 72 [15】Roberta T P Correia,Patrick Mccue,Dhi raj A.Vattem.et a1.amylase and helicobacter pylori inhibition by phenolic extracts of pineapple 收稿日期:2009—09—27 天然食品添加剂 食品研究与开发 第32011卷第9期 0年9月 229===一 technology was considered to be extracted with 50%ethanol at 53℃for 3 h and the maximum yield 0f the anthocyanins was 108.33 g儋.At the range of the designed condition the mathematic regression model could forecast the extraction yield of anthocyanins from hulless barley. Key words:hulless barley;anthocyanin;extraction technology;quadratic regression rotational combinational design 青稞(Hordeum vulgare ssp.vulgate)又称裸大麦,是 多功能”的天然色素资源提供参考。 禾本科大麦属作物,在世界许多国家和地区作为主要 粮食作物已有几千年历史。20世纪以后,随着稻米和 小麦的食用量增加,大麦的食用量急剧减少。然而20世 纪80年代以来,由于青稞丰富的营养成分和特有的保 健品质以及在燃料工业中的潜力,在发达国家日趋受 到重视,掀起综合研究利用的热潮【l_5】。我国青藏高原 地区是青稞的起源中心之一,拥有占全世界2/3以上 的裸大麦资源,尤其富藏黑色、紫色、灰色等深色品种, 如能科学合理的利用,将带来巨大的经济效益和社会 效益,且对藏区农业种植结构的调整都将具有潜在意义。 现代研究发现,黑色食品中的天然色素(花色苷 类物质)具有较强的消除体内自由基和抗氧化作用,具 有延缓衰老、预防癌症等功效 】。此外,相对于合成色 素,天然色素主要来源于动植物和微生物,具有安全、 健康的特点 ,利用天然物质提取食用、化妆品和医用 色素,成为当今世界色素研究和发展的主要趋势llq。裸大 麦种皮色素同属花色苷类色素,Deguchi T.等【l1]从紫色 裸大麦麸皮的发酵物中提取的花色素——单宁类紫色 素,可消除机体代谢过程中产生的过多自由基,具有较 强的抗氧化、抗突变和抗肿瘤等功能。目前关于对我国 青稞中天然色素的相关研究还鲜报道,但通过对从西 藏黑青稞种皮中的黑色素的研究发现,该类色素也具 有优异的抗氧化性f结果另文发表)。天然色素的提取 受多个因素的影响,为使提取达到最大效率,本试验采 用单因素试验设计与二次回归正交旋转组合设计相结 合的方法优化青稞色素提取工艺。二次回归正交旋转 组合设计是将回归正交设计与旋转设计组合起来的一 种方法,它不仅保留了回归正交设计的优点,如试验布 点分散均匀,试验次数少,还部分地消除了回归系数间 的相关性,虽然牺牲部分的正交性但是获得旋转性,可 直接找出预测值相对较优的区域,降低了试验成本,缩 短了研究周期。它以回归方程作为函数估算的工具, 克服了正交试验只能给出最佳因素水平组合,而无法 找出整个区域上因素的最佳组合和试验结果的最优值 的缺陷旧。本研究首次对青稞总花色素提取工艺进行 探索,为开发利用深色裸大麦作为新型“天然、营养、 1材料与方法 1.1材料与试剂 试验所用青稞材料由农科所唐亚伟惠赠。 HH一数显恒温水浴锅:上海逸龙科技有限公司; pHB一4便携式pH计:上海精密科技仪器有限公司; HR1727型搅拌机:飞利浦公司;SHZ—D型循环水式真 空泵:巩义市英峪予华仪器厂;AB104一N型分析天平: 梅特勒一托利上海有限公司;UV一7504型单光束紫外 可见分光光度计:上海欣茂仪器有限公司。 试验所用试剂均为国产分析纯。 1.2试验设计 1.2.1单因素试验设计 分别单独考察提取剂(不同浓度甲醇、乙醇和 水)、酸化剂(盐酸、三氯乙酸、酒石酸、甲酸、酒石酸、柠 檬酸)、提取温度(30、40、50、60、70、80、95℃)、提取时 间(1、2、3、4、5 h)、料液比(1:5、1:10、1:20、1:30、1:40、 1:50 g/mL)、提取次数(1、2、3、4)及粉碎度(完整籽粒, 20、40,60、80、100目)对青稞黑色素提取率的影响。 1.2.2二次回归正交旋转组合试验设计 根据单因素试验结果,选择对提取率影响显著的 因子(温度、时间、乙醇浓度)作为二次回归正交旋转组 合试验的因素并确定其零水平和间距。 1.3色素提取方法 称取除杂后的青稞种子3 g置于250 mL圆底烧 瓶中,加入一定体积的提取剂,于水浴锅中回流提取。 提取结束后,趁热过滤,并收集滤液,避光放置待测。 1.4色素提取率的测定 滤液避光冷却至室温后,将pH值调节为2,以提 取剂作为空白,于535 nm处测定吸光度,色素的提取 率按以下公式计算[13】: 色素提取率= 式中:色素提取率,( g); 35为535 nln下的吸 光度, 为稀释后的体积,mL;N为稀释倍数: l 000为 单位换算倍数(mg/g--%zg/g);98.2为花色素的平均消 —===230 1.5数据分析 第卷第期 31 9 ’食品研究与开发 瓦明研九司丌矾 2.1.2 pH值及酸化剂的选择 天然食品添加剂 光系数;m为样品重量,g。 花色苷在中性和弱碱性溶液中不太稳定,因此, 所有试验均设置3次重复,单因素试验数据用 spssl 1.5进行分析,二次回归正交旋转组合试验结果 用DPS v8.01软件进行分析。 2结果与分析 提取过程通常要采用酸性溶剂。酸性溶剂在破坏植物 细胞膜的同时溶解水溶性色素,但过低的pH值能改 变色素的原始状态,为了获得更接近于天然状态的花 色苷一般可采用有机酸作为酸化剂 。此外,根据青稞 色素的稳定性研究,其在pH 3时稳定性较好(另文报 2.1单因素试验 道),因此试验选择在此pH条件下,研究不同酸化剂 2.1.1提取剂的选择 溶剂法提取活性成分时,选择适宜的溶剂非常关 键。天然谷物色素属黄酮类花色苷色素,具有较强极 性。因此,在提取时问为1 h,料液比为l:lO(g/mL)的室 温条件下,对常用的极性溶剂对黑青稞籽粒中花色苷 的提取能力进行测定见图1。 O-15 o.1o 譬o.05 0l00 l_l 。- 甲醇7O%甲醇乙醇70%乙醇水 提取剂 圈1提取剂对黑青稞色素提取率的影响 Fig.1 Effects of solvent Oil extraction yield 结果如图l所示,乙醇一水溶液(70%乙醇)的提 取率最高,甲醇一水溶液(70%甲醇)和水溶液次之,纯 甲醇和乙醇的提取率极低。乙醇安全廉价,提取液黏 度小、易过滤、回收方便,不易发霉变质,所以选择乙 醇一水溶液(70%甲醇)作为提取剂。 进一步考察乙醇浓度对提取率的影响见图2。 立 褂 乙醇浓度/% 图2 乙醇浓度对黑青稞色素提取率的影响 Fig.2 Effects of ethanol concentration on extraction yield 结果发现,色素提取率随乙醇浓度的增加而增加, 当乙醇浓度达到60%时提取率最高,随后又逐渐降 低。青稞色素易溶于乙醇,在水中溶解度相对较低,但 水起着穿透植物细胞的作用,有助于花色苷从细胞中 渗出㈣,因此选用醇一水混合物可提高提取率。根据试 验结果,确定最佳乙醇浓度为60%。 对提取率的影响。 0 0.6 三0.4 甜 聋O.2 O f I.I.-。-.-.I 。. 酸化剂 图3不l司酸化剂对黑青稞色素提取翠的影响 Fig.3 Effect of acidifying agent on extraction yield 结果如图3所示,相同pH值条件下以有机酸为 酸化剂时提取率普遍高于盐酸,其中以柠檬酸提取率 最高。该结果与在黑米色素中的研究结果一致 ,可 能是由于柠檬酸酸性柔和,不会破坏花色素的结构,对 色素无不良影响 1。因此,确定用柠檬酸为酸化剂。 2.1.3提取温度的选择 提高温度是增大花色素溶解度和溶出速率的有 效途径,但温度过高也会加速活性成分的分解,尤其是 对黄酮类具有还原活性物质来说不宜用过高的温度 提取。此外过高的温度在实际生产中会使生产成本增 加。因此,选择合适的提取温度至关重要,结果见图4。 i 铸 温度/(℃ 图4提取温度对黑青稞色素提取率的影响 Fig.4 Effects of extraction temperature on extraction yield 如图4所示,黑青稞提取率在30℃~70℃随温度 升高呈增加趋势,而大于80℃后提取率反而有所下 降。这是由于温度过高情况下,色素发生一定程度的 降解,同时,温度较高促进可溶性蛋白质溶出变性,溶 液黏度增大影响细胞,阻碍了黄酮类物质的溶出, 从而降低了花色素的提取率㈣。综合考虑,确定最佳提 天然食品添加剂 食品研究与开发 第32010年9月 卷第9期 取温度为70℃。 2.1.4提取时间的选择 以pH 3(柠檬酸)的60%乙醇水溶液为提取剂, 提取率随时间的变化规律如图5所示。 瓣 时间/h 图5提取时间对黑青稞色素提取率的影响 Fig.5 Effects of extraction time on extraction yield 在0.5 h~5 h时提取率逐渐增加。天然有机化合物 在提取过程中,时间过短难以有效提取,而时间过长有 效成分可能被破坏,并且造成成本的浪费。花色苷类 物质其本身易被氧化破坏,过长时间提取会导致其受 热降解_191。因此将时间确定为时间为3 h。 2.1.5料液比的选择 料液比的选择见图6。 0. 驰 0. 料液比/(g/mL) 图6料液比对黑青稞色素提取率的影响 Fig.6 Effect of solid-solvent ratio on extraction yield 由图6可知,黑青稞色素提取率随料液比的增大 而升高,但当料液比大于1:lO(g/mL)后,提取率逐渐下 降。这是由于液固两相存在溶解平衡和吸附平衡,低 温有利于吸附,由于料液比的增大,在产生相同热量的 条件下,提取液温度下降,吸附作用增强导致提取率下 降 ,同时考虑到从生产效率,耗能和节约成本等各方 面因素,确定l:10(g/mL)应为最佳的料液比。 2.1.6提取次数的选择 提取次数的选择见图7。 《1¨¨ 2. ==_- 如图7所示,青稞籽粒在60℃提取1 h时,提取 率随提取次数增加而提高。假定提取4次的色素总量 为100%,则1、2、3次提取率分别为85.1%、96.2%、 99.7%,差异不显著(P>0.05 o因此从节约成本,生产 效率等各方面考虑确定进行一次提取。 2.1.7籽粒粉碎度的选择 青稞籽粒含有50%以上的淀粉,在温度较高的条 件下发生糊化【2I】’造成粉碎籽粒的提取液浑浊黏稠,难 以过滤及测定。此外青稞色素主要存在于其:种皮之中, 因此确定以完整籽粒进行色素提取。 2.2三元二次回归正交旋转组合设计的试验结果 对单因素试验结果进行方差分析(ANOVA),发现 提取时间、温度、和乙醇浓度对提取率的影响显著(p< 0.05),因此以这3个因素作为正交试验研究对象,固 定其它试验条件为:以酸化剂为柠檬酸,料液比1:10 (g/mL),完整籽粒提取一次。取试验中效果较好的水平 最为各正交试验研究因素的零水平,并设计因素水平 编码表,见表1。 表1因素水平表 Table 1 FactOrs and levels of the test 规范变量 提取温度/℃ 提取时间/h 乙醇浓度/% 2.2.1正交试验设计表及试验结果 采用DPS软件编制二次正交旋转组合设计表,并 将试验结果列入表2。 表2黑青稞色素的回归旋转组合设计表和色素提取率 Table 2 Quadratic regression rotational combinational design and extraction yield X 提取温廖 提取时 X 乙醇浓度 提取率H 恸 232 第卷第期 31 9 ’食品研究与开发 R明耵九司丌鼠 天然食品添加剂 续表2黑青稞色素的回归旋转组合设计表和色素提取率 Continue table 2 Quadratic regression rotational combinational design and extraction yidd 2.2.2回归方程的建立和检验 利用DPS软件对试验结果进行统计分析,得到影 响多糖提取率的回归方程: Y=92.894 77—7。51 1 87X1—3.053 87X2—5.930 54X3— 1.370 62 —3.668 72X22—3.133 09X32—0.483 75 1X2+ 1.048 75XlX厂4.1 33 75X2X3 式中:y为色素提取率(I.zg/g); 为提取温度,qC; 为提取时间,h; 为乙醇浓度,%。 表3回归旋转组合设计试验方差分析表 Table 3 Variance analysis of quadratic regression rotational combinational design 表3由方差分析结果可知,回归方程的失拟性检 验:F >Fn∞(5,8)-1.64不显著,即所选的回归模型不存 在失拟因子,回归模型是合适的。对回归方程进行显 著性检验:F2>F0.∞f9,13)=8.27极显著,说明得到的回 归方程与试验数据是吻合的,该数学模型成立。对回 归系数进行显著性检验,在 :0.1显著水平剔除不显 著项,得到优化后的回归方程: Y=92.894 77—7.5 1 1 87X。一3.053 87X2—5.930 54X3— 3.668 72 ’2—3.133 09X32-4.133 75X 2.2.3因素重要性分析 为进一步确定各因素色素提取率的影响程度,采 用贡献率法 对二次正交旋转组合设计结果进行分析。 值的标准转换值计算公式如下: 0 若F≤1 1一__1若F≥1 贡献率计算公式如下: m △ + 岛蚂 i≠j 式中:△ 为第. 因素对提取效果的贡献的大小;△, 代表第 .因素F值的标准转换值;△ 代表第i因素与 第_『因素交互作用F值的标准转换值 代表表第 因 素二次项F值的标准转换值。 据表3均方比F值,可知温度、时间、乙醇浓度对 提取率的贡献率分别为△ =1.108 3,△。=2.082 266,△,= 2.186 386。由于△ >△ >△ ,所以各因素对提取率的作用 大小依次为乙醇浓度>时间>温度。 2.2.4提取条件的优化和验证 优化提取方案中x 取值频率分布表见表4。 表4优化提取方案中x。取值频率分布表 Table 4 The probability distribution of xl in the combined application DPS软件分析可得,黑青稞色素提取率最高时提 取方案为xl—1.682,X:,--0,x 1,即温度53 cI二,时间3h, 乙醇浓度50%。此时提取率为108.33 lxg/g,明显高于 原试验设计中所有处理组合的提取率。在此条件下对 结果进行验证,实际提取率为105.27、106.11、 105.54 ,与理论值间误差小于3%,说明该回归方 程与实际情况吻合良好。 本试验为三因素五水平,共有53=125种组合方 式,其中大于88.04的方案有45个,各变量取值的频 率分布如表4所示。在95%的置信区间内提取率大于 88.04的优化提取方案为:温度57.20℃~61.19 cC,时间 2.528 h~3.189 h,乙醇浓度51.53%~57.44%。 天然食品添加剂 食品研究与开发 第32011卷第9期 0年9月 233 取优化后浸提条件的平均值方案,X :一1.080 5, X2=一0.141 5,X3:一0.551 5,即温度为59-2 cc,时间为 川.Journal ofCereal Science.1996,24(2):16l—l f4】Thomas K C,Dhas A,Rossnagel B G.Production of fu el alcohol from hull-less barley by very high gravity technology[J].Cereal chemisty,r 1995,72(4):360-364 2.9 h,乙醇浓度为54.5%组合方案对优化结果验证。 此时理论值Y=103.37 Ixg/g,与实测值100.06 {;/g非常 接近(误差为3.2%),说明该回归模型可以较好的反 映实际规律,也进一步验证了其合理性。 3结论 【5]尼玛扎西.青稞与高原地区的食物保障:优势与作用[J】.农业 科技,1998,20(2):20—25 【6】付湘晋.我国天然色素研究概况[J】.粮食与油脂,21X D5(12):42—45 [7]王桃云,万承永,王飞.我国天然色素研发与初探【 浙江林业科 技,2001,21(2):71—73 【8]邓祥元,王淑军,李富超,等.天然色素的资源和应/- ̄[J1.中国调味 品.2000(10):49—53 【9】Bridle P,Timberlake C F.Anthocyanins as natural food colours— selected aspects[J】.Food Chemistry,1997,58(1/2):103—109 本研究采用经典溶剂提取法提取黑青稞色素,对 仪器设备要求不高。在设计试验时,考虑到众多因素 都对提取效率有着不同程度的影响,因此先进行了单 因素试验,筛选出对提取效率影响显著的因素,再通过 二次正交旋转组合试验确定最佳提取工艺,有效提高 了试验的严谨性和可靠程度。试验结论如下: 黑青稞中花色苷类色素提取条件的优化数 学回归模型为: Y=92.894 77-7.5 1 l 87Xl一3.053 87X2-5.930 54X3— 3.668 72X22—3.133 09X3。4.133 75X2X3 [10】左玉.天然色素的研究与应用[J】.粮食与油脂,2006(19):46—48 [1 l】Deguchi L Ohba R,Ueda S.Radical scavenging activity of a purple pigment,hordeumin,from uncooked barley bran-fern1ented broth[J]. Journal ofAgricultural and Food Chemistry,2000,48(8):3198-3201 【12】杨德.试验设计与分析【M】.北京:中国农业出版社,2002:221 [13]徐金瑞,张名位,刘兴华.黑大豆种皮花色苷的提取及其抗氧化作 用研究【J】.农业工程学报,2005,21(8):t61—164 [14】唐洪波,马冰洁,李艳萍.葡萄籽中提取原花色素工艺的研究[J1. 粮油加工。2008(10):73—76 由回归方程可知,在酸化剂为柠檬酸(pH 3),料 液比1:10 g/mL,完整籽粒提取一次的条件下,黑青稞 【15】卢钰,董现义,杜景平,等.花色苷研究进展【J】l山东农业大学学报: 自然科学版,2004,35(2):315—320 花色素最佳提取工艺条件为:53℃,时间3 h,乙醇浓 度5O%,各因素对提取率的作用大小依次为乙醇浓 度>时间>温度。此数学回归模型在试验范围内能较准 [16】于长顺,杨婷婷,马春,等.黑米色素提取精制工艺的研究[J】_大连 轻工业学院学报,2005,24(2):l28—13l 【l7】张晴.黑米色素提取精制工艺的研究【J1.青岛大学学报:T程技术 版。2000.15(2):24—26 确的预测总花色素的提取率。 参考文献: 【1】Ehrenbergerov&J,Belcrediov6 N,Pryma J,et a1.Effect of cuhivar, year grown,and cropping system on the content of tocopherols and 【l8】严赞开,朱兰花.炮仗花色素的提取及稳定性[J】.食品科技,2007, 32(1):137-139 [19】唐浩国.黄酮类化合物研究[M】.北京:科学出版社,21)09:39—40 【2O]冯年平,郁威.中药提取分离技术原理与应用[MI.北京:中国医药 科技出版社,2005:6—7 tocotrienols in grains of hulled and hulless barley[J1.Plant Foods for Human Nutirtion,2006,61(3):145—150 [2l】邹弈星,潘志芬,邓光兵,等.青藏高原青稞的淀粉特性[J1.麦类作 物学报,2008,28(1):74—79 [2】Liu C,Wesenberg D,Hunt C,et a1.Hulless barley:A new look for barley in Idaho【J/OL].1996.http://WWW.cals.uidaho.edu/edcomm/ pdf/CIS/CIS1050.pdf 【221巩发永,齐桂年,李静.四川边茶中茶多糖提取条件的优化研究 [J】.福建茶叶,2005,25(3):6—9 [3]Oscarsson M,Andersson R,Salomonsson A C P.Chemical composition of barley samples focusing on dietary fibre components 收稿日期:2009—11-07 ●●●●●●●●●●●●…●●●●●●●●●●●●●●●●●●●●●●◆●●●●●●●●●…●一●●●● 《食品研究与开发》编辑部常年办理订阅手续 地址:天津市静海经济开发区(天字园)科技大道9号 邮编:301609 邮汇收款人:《食品研究与开发》编辑部 电话:022—5952567 1