• EI
  • Scopus
  • 中国科技期刊卓越行动计划项目资助期刊
  • 北大核心期刊
  • DOAJ
  • EBSCO
  • 中国核心学术期刊RCCSE A+
  • 中国精品科技期刊
  • JST China
  • FSTA
  • 中国农林核心期刊
  • 中国科技核心期刊CSTPCD
  • CA
  • WJCI
  • 食品科学与工程领域高质量科技期刊分级目录第一方阵T1
中国精品科技期刊2020

酶解辅助预糊化技术工艺优化及风味物质差异分析

昝学梅, 刘明, 刘艳香, 孟宁, 刘翼翔, 谭斌, 田晓红, 方秀丽

昝学梅,刘明,刘艳香,等. 酶解辅助预糊化技术工艺优化及风味物质差异分析[J]. 食品工业科技,2021,42(8):19−28. doi: 10.13386/j.issn1002-0306.2020060094.
引用本文: 昝学梅,刘明,刘艳香,等. 酶解辅助预糊化技术工艺优化及风味物质差异分析[J]. 食品工业科技,2021,42(8):19−28. doi: 10.13386/j.issn1002-0306.2020060094.
ZAN Xuemei, LIU Ming, LIU Yanxiang, et al. Process Optimization of Enzyme-assisted Pregelatinization Technology and Analysis of Flavor Substance Difference[J]. Science and Technology of Food Industry, 2021, 42(8): 19−28. (in Chinese with English abstract). doi: 10.13386/j.issn1002-0306.2020060094.
Citation: ZAN Xuemei, LIU Ming, LIU Yanxiang, et al. Process Optimization of Enzyme-assisted Pregelatinization Technology and Analysis of Flavor Substance Difference[J]. Science and Technology of Food Industry, 2021, 42(8): 19−28. (in Chinese with English abstract). doi: 10.13386/j.issn1002-0306.2020060094.

酶解辅助预糊化技术工艺优化及风味物质差异分析

基金项目: “十三五”国家重点研发计划(2018YFD0400802)
详细信息
    作者简介:

    昝学梅(1993−),女,硕士研究生,研究方向:食品工程,E-mail:17835396976@163.com

    通讯作者:

    刘翼翔(1982−),男,博士,副教授,研究方向:从事海洋生物活性物质利用与开发研究,E-mail:lyxjmu@jmu.edu.cn

    谭斌(1972−),男,博士,研究员,研究方向:粮食资源深度开发利用,E-mail:tb@chinagrain.org

  • 中图分类号: TS210.1

Process Optimization of Enzyme-assisted Pregelatinization Technology and Analysis of Flavor Substance Difference

  • 摘要: 利用响应面试验设计优化糙米的酶解辅助预糊化工艺,以蒸煮时间为指标,得到最佳酶解辅助预糊化工艺条件,用气相色谱–离子迁移谱法(GC-IMS)测定白米、糙米和酶解辅助预糊化糙米及其蒸煮前后的挥发性物质。结果表明,最优工艺条件为酶解时间为4 h、酶的添加量为0.1%、预糊化时间为10 min,酶解温度为39 ℃。在酶解辅助预糊化糙米中检测到56种挥发性物质,其中醛类19种,酮类8种,酯类6种,醇类12种及其它类11种;而蒸煮后的酶解辅助预糊化糙米中检测到44种挥发性物质,其中醛类18种,酮类8种,酯类3种,醇类11种及其它类4种,且2-辛烯醛、苯乙醛、苯乙酮、乙偶姻、γ-丁内酯、戊酸、己酸、丁酸和己二烯二硫化等挥发性物质在蒸煮之后消失。糙米酶解辅助预糊化前、后挥发性风味成分改变最大的为醛类、酮类和醇类;而酶解辅助预糊化糙米蒸煮前、后挥发性风味成分改变最大的为酯类和其它类中的酸类。
    Abstract: The response surface experiment was used to design and optimize the enzyme-assisted pregelatinization process of brown rice. The optimum conditions of enzyme-assisted pregelatinization process were obtained through the cooking time: The enzymatic hydrolysis time was 4 h, the amount of enzyme added was 0.1%, the pregelatinization time was 10 min and the enzymatic hydrolysis temperature was 39 ℃. Gas chromatography-ion migration spectrometry (GC-IMS) was used to determine the volatile substances of white rice, brown rice and enzyme-assisted pregelatinized brown rice before and after cooking. 56 kinds of volatile substances were detected in enzymolysis assisted pregelatinized brown rice, including 19 kinds of aldehydes, 8 kinds of ketones, 6 kinds of esters, 12 kinds of alcohols and 11 kinds of other compounds. Forty-four kinds of volatile substances were detected in the pregelatinized brown rice assisted by enzymatic hydrolysis after cooking, including 18 kinds of aldehydes, 8 kinds of ketones, 3 kinds of esters, 11 kinds of alcohols and 4 kinds of other compounds, and 2-octenal, phenylacetaldehyde, acetophenone, acetoin, gamma-butyralide, pentanoic acid, hexanoic acid, butyric acid and hexadiene disulfide and other volatile substances disappeared after cooking.Aldehydes, ketones and alcohols showed the largest changes in the volatile flavor components of brown rice before and after the enzymatic hydrolysis assisted pregelatinization. The acids in esters and other groups showed the greatest changes in the volatile flavor components of brown rice before and after cooking assisted by enzymatic hydrolysis.
  • 稻谷属于禾本科稻属,是水稻脱粒后得到的带有不可食颖壳的籽粒;经砻谷处理,脱去颖壳,得到的籽粒称为糙米;再经过碾米加工得到大米[1]。但由于全谷物糙米灰暗的外表和坚硬的质地使其蒸煮性、口感和吸收性较差[2],所以不能作为现在的日常主食。而精白米虽然口感好、容易消化[3],但缺乏有益于人体的健康成分,如维生素和纤维素等。

    酶是由生物体活细胞产生的且具有特殊催化作用的一类蛋白质[4],目前用于酶解颗粒态糙米的酶主要是纤维素酶和纤维素酶与木聚糖酶的混合酶等。这主要是因为糙米皮层主要成分为纤维素和阿拉伯木聚糖,木聚糖层与木质素共价连接、与纤维素非共价连接,此种交错联结的形式是赋予糠层抗性的基础[5]。而酶解技术不仅可以破坏全谷物糙米的组织细胞壁,而且还可以改善糙米的蒸煮特性和食用品质等。Das等[6]通过研究木聚糖酶和纤维素酶混合酶溶液对糙米皮层的影响,研究表明酶解糙米的水摄取率和膨胀率增加,糙米的表面结构变化,使糙米的蒸煮时间减少,糙米蒸煮食用品质提高。预糊化是一种物理改性方法[7],淀粉经预糊化处理后具有绿色环保、分散性好、膨润度高等特点,不仅可以作为增稠剂应用于食品行业,还可作为粘结剂、施胶剂等应用在医药化工等领域[8-10]。研究表明,将预糊化淀粉加入面团中,其强吸水性及结构的无序化可以改变面制品的醒发体积及凝胶特性,调控其品质[11-12]。酶解辅助预糊化技术是先利用添加外源酶对糙米表层的纤维素等物质进行降解,再通过预糊化对糙米胚乳进行物理改性的一种复合加工技术。目前,酶解技术和预糊化技术多用于改善糙米的蒸煮特性和食用品质,尚未见酶解辅助预糊化的复合加工技术在全谷物蒸煮品质和风味方面的研究报道。

    风味是全谷物最重要的感官品质之一,挥发性风味物质刺激鼻腔嗅觉感受器直接影响全谷物品质评价和消费者接受能力[13]。迄今为止,在各种水稻中已鉴定出300多种挥发性化合物,包括酮类、醛类、酯类、醇类、杂环类化合物[14]。因此,本试验中采用响应面试验设计研究糙米酶解辅助预糊化工艺参数(酶解时间、酶的添加量、预糊化时间和酶解温度),采用蒸煮时间评价筛选出品质最优工艺条件,并用气相-离子迁移谱(GC-IMS)检测酶解辅助预糊化糙米蒸煮前后的挥发性风味物,旨在为糙米加工过程中挥发性成分变化规律提供参考。

    黑龙江长粒香(HCLX) 五常市金禾米业有限责任公司;纤维素酶(全组分PL) 宁夏夏盛实业集团有限公司提供。

    GC(7890A)-MS(5975C)气质联用分析仪 美国安捷伦公司;50/30 μm DVB/CAR/PDMS固相微萃取头  美国Supelco公司;RDN型分段可编程人工气候箱 宁波东南仪器有限公司;HH-4数显恒温水浴锅 国华电器有限公司。

    稻谷经砻谷机砻出糙米,去除霉变糙米粒、杂质、未成熟粒和碎糙米粒,用聚乙烯自封袋封装保存在4 ℃冰箱中备用。

    每组称量糙米200 g,放入烧杯中,用纯水洗涤3次,然后用0.5%的次氯酸钠溶液消毒15 min,之后再用纯水洗涤3次,按固液比1:3 g/mL添加酶溶液(柠檬酸缓冲液pH5.0),然后将烧杯放入恒温水浴锅中。在设定温度和时间下进行酶解,每15 min 翻料一次,糙米经酶解及预糊化后,90 ℃烘箱烘1 h。晾凉后−4 ℃冷藏保存。

    工艺流程:糙米除杂筛选→酶解→预糊化→干燥→备用。根据预实验结果,以原料糙米为对照,分别在酶解温度为40 ℃,酶的添加量为0.1%,预糊化时间为9 min的固定条件下,考察酶解时间(3、4、5、6 h)对处理糙米蒸煮时间的影响;在酶解时间为4 h,酶解温度为40 ℃,预糊化时间为9 min的固定条件下,考察酶的添加量(0.05%、0.1%、0.15%、0.2%)对处理糙米蒸煮时间的影响;在酶解温度为30 ℃,酶的添加量为0.1%,酶解时间为4 h的固定条件下,考察预糊化时间(3、6、9、12 min)对处理糙米蒸煮时间的影响;在酶解时间为4 h,酶的添加量为0.1%,预糊化时间为9 min的固定条件下,考察酶解温度(20、30、40、50 ℃)对处理糙米蒸煮时间的影响。

    利用Design-Expert 8.0.5软件自带的中心组合试验模型,以酶解时间(A)、酶的添加量(B)、预糊化时间(C)、酶解温度(D)4个因素为自变量,糙米处理后蒸煮时间为响应值,依据单因素实验结果,设计了四因素三水平的响应面分析试验。利用Design-Expert 8.0.6软件进行数据拟合,具体因素水平设置表1

    表  1  响应面试验因素和水平
    Table  1.  Factors and levels of Box-Behnken experiment
    水平因素
    A酶解时间
    (h)
    B酶的添加量
    (%)
    C预糊化时间
    (min)
    D酶解温度
    (℃)
    −130.05630
    040.10940
    150.151250
    下载: 导出CSV 
    | 显示表格

    使用玻璃板-白芯法进行米饭蒸煮时间的测定[15]。称取5 g样品,加入到10 mL沸水中进行蒸煮。蒸煮时,每隔2 min随机取出10粒样品,将其在玻璃板上进行挤压。样品挤压时无不透明白芯出现的时间为蒸煮时间。

    6种待测样品各称取20.00 g于顶空样品瓶中,用锡纸做隔垫进行加塞密闭处理,于80 ℃的恒温水浴中平衡1 h,再将萃取头插入顶空样品瓶中萃取50 min,之后于GC-IMS进样口(250 ℃、不分流模式)解析5 min。

    糙米顶空瓶中孵化温度:75 ℃/95 ℃,孵化时间20 min;IMS:迁移谱的温度45 ℃,GC:柱温40 ℃。其中采集数据程序设置如表2所示。

    表  2  采集数据程序设置
    Table  2.  Collect data application settings
    TimeE1(mL/min)E2(mL/min)
    00:00,001502
    02:00,001502
    10:00,0015010
    20:00,00150100
    25:00,00150130
    30:00,00150130
    下载: 导出CSV 
    | 显示表格

    试验中采用SPSS19.0进行数据分析,以P<0.05作为差异显著性判断标准,每个样品做3个重复。通过GC-IMS自带的IMScontrol TFTP-Server软件进行三维指纹图谱的采集,使用 GC-IMS自带的LAV 2.0.0软件和Matlab 2015b(美国 Mathworks)进行综合分析。使用LAV对谱图进行RIP(Reactive ion peak)标准化处理,依据不同样品对应的特征峰强度变化,选取26个特征信息,以特征峰对应的离子强度作为参数变量,使用Matlab软件对变量进行统计分析,主要采用主成分分析和线性判别分析两种方法进行分析处理。GC-IMS特征峰选取原则:GC-IMS三维谱中每一个特征峰代表一种挥发性有机成分,以不同样品中的特征峰的有无或者特征峰强度的显著变化为原则,手动选择并以长方形标记特征峰所在区域,以该特征峰区域的峰强度(即特征区域离子强度的最大值)作为参数变量进一步做数据分析。

    图1A所示,随着酶解时间的延长,糙米蒸煮时间呈现先下降后上升的趋势,在酶解4 h时蒸煮时间最短。其原因可能是随酶解时间的延长,糙米皮层的致密结构不断被破坏,有利于糙米蒸煮时间降低。但继续延长酶解时间却使糙米的蒸煮时间增加,可能是因为酶解液随着酶解时间的延长,组分越来越复杂,不仅降低酶活性限制酶解反应,而且皮层结构破坏程度越大越有利于酶解液中的小分子与糙米中成分反应,限制糙米糊化,从而使糙米的蒸煮时间增加。这说明继续增加水解时间不能提高酶解效率,故酶解时间选为4 h。

    图  1  酶解时间(A)、酶的添加量(B)、预糊化时间(C)、酶解温度(D)对糙米蒸煮时间的影响
    注:小写字母不同表示差异显著,P<0.05。
    Figure  1.  Effect of processing time(A), enzyme addition(B), steaming time(C) and temperature(D) on the cooking time of brown rice

    图1B显示了酶的添加量对糙米蒸煮时间的影响规律,随着酶的添加量的增加,糙米蒸煮时间呈先下降后上升的趋势。酶的添加量低时,酶解速率较慢,皮层降解程度低,蒸煮时水分进入糙米较慢。随着酶的添加量的增加酶解速率增大,皮层降解程度加大,水分更容易与淀粉作用糊化,使糙米的蒸煮时间缩短。然而,过高的添加量使底物浓度受到制约,致使反应基本达到平衡或降低,使得蒸煮时间增加。

    预糊化时间对糙米蒸煮时间的影响如图1C所示,随着预糊化时间的增加,糙米的蒸煮时间呈显著的下降趋势(P<0.05)。这可能是预糊化过程中的蒸汽处理会对糙米致密的表面糠层结构、内部结构产生较强的破坏作用,且蒸汽处理时间越长,结构破坏越严重,有效促进了糙米蒸煮过程中水分的吸收[16],从而缩短了蒸煮时间。但是,当预糊化时间为12 min时,糙米饭的爆腰率明显增多,所以预糊化时间选为9 min。

    酶解温度对糙米蒸煮时间影响如图1D所示,随着酶解温度的升高,糙米的蒸煮时间不断降低,且降低趋势明显,在40 ℃时达到最低。这是因为随着温度的升高纤维素酶的活性提高,糙米皮层降解吸水速率加快,蒸煮时促进糙米吸水糊化,使糙米蒸煮时间降低。

    根据单因素实验的结果,以酶解时间、酶解温度、酶的添加量和预糊化时间为变量因素,蒸煮时间为响应值,采用 Box-Benhnken 法进行响应面分析,实验结果见表3

    表  3  中心组合实验设计及结果
    Table  3.  Central composite design and results
    实验号A-酶解时间B-酶的添加量C-预糊化时间D-酶解温度Y-蒸煮时间(min)
    1−10−1026.28
    200−1126.20
    3000021.57
    4−101024.57
    50−10127.24
    601−1025.97
    7100126.23
    8010−125.29
    9−1−10027.26
    10000021.02
    111−10024.95
    120−10−126.53
    13100−124.23
    140−11024.80
    15−110026.20
    16000021.26
    17110024.65
    18011025.64
    1900−1−125.69
    20010125.83
    21−100126.25
    22000021.34
    23101024.10
    240−1−1027.62
    25000022.56
    26001124.96
    2710−1025.67
    28001−125.25
    29−100−125.34
    下载: 导出CSV 
    | 显示表格

    对实验方案得出的29组数据进行多元回归拟合,得到以酶解辅助预糊化糙米的蒸煮时间为目标函数的响应面回归方程为:

    Y=21.55−0.51A−0.40B−0.68C+0.36D+0.19AB+0.035AC+0.27AD+0.62BC−0.042BD−0.20CD+ 1.74A2+2.52B2+1.87C2+2.16D2

    表4可以看出,蒸煮时间的回归模型P值小于0.0001,说明回归模型达到极显著水平。酶解时间和预糊化时间两个因素的单次项对蒸煮时间影响有极显著水平且P<0.01;酶的添加量和酶解温度两个因素的单次项对蒸煮时间影响显著(P<0.05)。而酶的添加量、酶解时间、酶解温度和预糊化时间4个因素的各平方项对蒸煮时间的影响均有极显著水平(P<0.01),而交互项中除了酶的添加量和预糊化时间交互项对响应值的影响显著外(P<0.05),其它交互项均不显著。在所选的各因素水平范围内,按照对蒸煮时间结果的影响排序为:C>A>B>D,即预糊化时间>酶解时间>酶的添加量>酶解温度。失拟系数0.6157,达到非显著水平,说明模型合适;蒸煮时间的决定系数R2分别为0.9514。说明了方程拟合度较高,即方程中自变量的变化可以很好地解释因变量Y的变化。

    表  4  回归方程方差分析表
    Table  4.  Analysis of variance (ANOVA) of regression equations
    方差来源平方和df均方 F P值
    模型88.33146.3119.58<0.0001**
    A-酶解时间3.0713.079.530.0080**
    B-酶的添加量1.9311.936.010.0280*
    C-预糊化时间5.4815.4817.010.0010**
    D-酶解温度1.6011.604.960.0428*
    AB0.1410.140.450.5141
    AC4.900E-00314.900E-0030.0150.9037
    AD0.3010.300.920.3533
    BC1.5511.554.810.0457*
    BD7.225E-00317.225E-0030.0220.8831
    CD0.1610.160.500.4926
    A219.71119.7161.07<0.0001**
    B241.34141.34128.31<0.0001**
    C222.70122.7070.46<0.0001**
    D230.18130.1893.68<0.0001**
    残差4.51140.32
    失拟项3.08100.310.860.6157
    纯误差1.4340.36
    总和92.8428
    注:*差异显著,P<0.05;**差异极显著,P<0.01。
    下载: 导出CSV 
    | 显示表格

    酶的添加量、酶解时间、酶解温度和预糊化时间交互作用对蒸煮时间的影响结果见图2。从图2可以看出,在恒定因素条件为酶解温度为40 ℃、酶的添加量0.1%、酶解时间4 h和预糊化时间9 min的条件下,酶的添加量与预糊化时间交互作用均为显著(P<0.05)。随着酶的添加量和预糊化时间的增加,酶解辅助预糊化糙米的蒸煮时间呈均先下降后上升;而酶解时间与酶的添加量、酶解时间与预糊化时间、酶解时间与酶解温度的交互作用、预糊化时间与酶解温度及酶的添加量与酶解温度交互作用均不显著(P>0.05)。这五种交互条件下随着酶解时间、酶的添加量、预糊化时间和酶解温度的增加,糙米的蒸煮时间均先下降后上升。

    图  2  各因素交互作用对糙米蒸煮时间影响的响应面图
    Figure  2.  Response surface plots for the effects of different factors onthe cooking time of brown rice

    分析回归模型可获取响应值蒸煮时间的极值点,即获得优化后的最佳工艺参数:酶解时间4.15 h、酶的添加量为0.1%、预糊化时间为9.5 min和酶解温度为39.14 ℃时,此时模型预测蒸煮时间的值为21.43 min。根据实际可操作性,按酶解时间4 h、酶的添加量为0.1%、预糊化时间为10 min和酶解温度为39 ℃的条件开展三次平行验证试验,所得蒸煮时间平均值为(21.56±0.23)min,与模型预测值的误差在5%以内,说明采用响应面法优化所得工艺参数较为准确可靠。

    根据挥发性成分气相色谱保留时间和离子迁移时间对样品挥发性组分进行定性分析,确定了58种挥发性物质的化学成分见表5。这58种物质主要是醛类、酯类、醇类、酮类、杂环类化合物、酸类及醚类等,与糙米挥发性成分组成的相关报道一致[17-18]。采用GC-IMS Library Search软件内置的2014NIST数据库和IMS数据库对6种样品的风味物质进行定性分析,在酶解辅助预糊化糙米中检测到56种挥发性物质,其中醛类19种,酮类8种,酯类6种,醇类12种及其它类11种;而蒸煮后的酶解辅助预糊化糙米中检测到44种挥发性物质,其中醛类18种,酮类8种,酯类3种,醇类11种及其它类4种。

    表  5  不同样品挥发性物质种类及峰强度信息
    Table  5.  Volatile substance types and peak strength information of different samples
    编号化合物CAS号保留指数气味特征保留时间(s)迁移时间(ms)
    醛类
    1苯甲醛二聚体100-52-7958.8坚果味、苦味434.6551.467
    2苯甲醛单体100-52-7958.8坚果味、苦味434.6551.148
    3己醛二聚体75-07-0793.6青草香251.161.563
    4己醛单体75-07-0797.2青草香254.0851.253
    5糠醛98-01-1958.8甜杏仁香烤味434.6551.148
    6戊醛二聚体110-62-3698.8木香、水果香182.3251.182
    7戊醛单体110-62-3902.9木香、水果香356.6551.331
    8辛醛二聚体124-13-01006青草、柑橘味516.9451.821
    9辛醛单体124-13-01006.6青草、柑橘味518.1151.404
    10庚醛二聚体111-71-7901.3水果香354.7051.697
    11庚醛单体111-71-7797.2水果香254.0851.253
    12正壬醛二聚体124-19-61105.5柑橘、肥皂味729.8851.943
    13正壬醛单体124-19-61107.4柑橘、肥皂味733.981.475
    143-甲基丁醛二聚体590-86-3658.3麦芽味162.4351.404
    153-甲基丁醛单体590-86-3656.3麦芽味161.6551.172
    162-辛烯醛(反式)单体2548-87-01063.7青草、脂肪味636.871.331
    17反式2-庚烯醛18829-55-5958.8辛辣、脂肪味434.6551.666
    18癸醛112-31-21209柑橘香960.961.538
    192-辛烯醛(反式)2548-87-01063.1青草、脂肪味635.5051.818
    20苯乙醛122-78-11043.8青草香味593.7751.253
    酮类
    21苯乙酮98-86-21057.8香辣杏仁味6241.187
    222-3-丁二酮431-03-8625.1清香味149.9551.168
    232-丁酮二聚体78-93-3587.2类似薄荷香味137.4751.243
    242-丁酮单体78-93-3589.6类似薄荷香味138.2551.057
    25丙酮67-64-1504类似薄荷香味110.761.119
    26己酮591-78-6783.4果味、黄油味242.971.188
    272-庚酮二聚体110-43-0893梨香344.9551.631
    282-庚酮单体110-43-0892.9梨香344.761.261
    292-戊酮107-87-9690.4类似薄荷香味177.6451.12
    酯类
    30乙酸乙酯二聚体141-78-6614.2水果香味146.251.335
    31乙酸乙酯单体141-78-6616水果香味146.8351.093
    32二氢草莓酸乙酯108-64-5831.5果香味283.141.257
    33乙酸己酯142-92-7987.4果香,似菠萝香482.6251.414
    34γ-丁内酯96-48-0919.6奶油焦糖味377.911.081
    35乙偶姻513-86-0719.2甜奶油味1951.055
    醇类
    36庚醇单体53535-33-4977.8甜香、坚果味465.8551.397
    37正己醇二聚体111-27-3873.9果香味323.8951.644
    38正己醇单体111-27-3877.7果香味327.991.327
    49正戊醇单体71-41-0768.4水果味231.0751.25
    40正戊醇二聚体71-41-0766.9水果味229.9051.51
    41甲醇二聚体3857-25-8958.9甜焦糖味434.851.565
    42甲醇单体3857-25-8961.1甜焦糖味438.361.254
    43反式2-己烯-1-醇二聚体2305-21-7849.6果香299.7151.516
    44反式2-己烯-1-醇单体2305-21-7849.6果香299.7151.516
    452-己醇626-93-7820.7果香味273.781.278
    46蘑菇醇3391-86-4984.5蘑菇香477.5551.155
    472-辛醇6169-06-8948果香味418.081.443
    其它类
    48正基呋喃3777-69-3993.7果香味494.131.253
    491,8-桉叶素单体470-82-61031.4樟树脑味567.841.295
    50柠檬烯138-86-31029.7柠檬香564.331.218
    51α-蒎烯80-56-8974.1松木香459.6151.216
    521,8-桉叶素二聚体470-82-61031.4樟树脑味567.841.729
    53二甲苯95-47-6923天竺葵382.3951.064
    54戊酸二聚体109-52-4885.9酸味336.961.509
    55戊酸单体109-52-4888.4酸味339.691.233
    56己酸142-62-1888.4干酪酸味339.691.233
    57丁酸107-92-6777.8奶酪酸味238.4851.162
    58己二烯二硫化2179-57-91076.2大蒜味664.3651.196
    下载: 导出CSV 
    | 显示表格

    醛类是本实验样品中检出种类数最多的挥发性物质,其中,未蒸煮样品中检出19种,蒸煮样品中检出18种,结果见图3。醛类化合物阈值一般较低,又具有叠加效应,往往在食品风味中起重要作用[19]。酶解辅助预糊化糙米(CHLXM)中的2-辛烯醛、3-甲基丁醛和糠醛含量相对较高,且可分别产生青草香、麦芽香和甜杏仁味,对CHLXM的整体愉悦气味有较大贡献。

    图  3  蒸煮与未蒸煮样品中醛类物质的气味指纹图谱
    Figure  3.  Major volatile components of aldehydes in cooked and uncooked samples

    CHLXM经蒸煮之后(CHLXM-ZHS),醛类物质的含量有所增加,其中包括苯甲醛、己醛、辛醛、庚醛、2-辛烯醛、2-庚烯醛和癸醛。其中,2-庚烯醛和2-辛烯醛等醛类物质是煮制温度达到一定条件才检测到,预测其原因是蒸煮达到了反应所需条件。大量醛类化合物的产生是因为在风味形成的过程中,当美拉德反应生成含有羰基的戊糖和丁糖时,进一步发生Strecker降解反应也可产生醛类[20],如苯甲醛是苯丙氨酸的降解产物[21]

    图4所示,CHLXM中的酮类化合物主要有苯乙酮、丙酮和2-丁酮,经蒸煮之后形成的CHLXM-ZHS样品中,丙酮和2-丁酮的含量均有所增加,说明蒸煮激活了糙米胚中丙酮酸转氨酶的活力和脂肪酸的氧化作用。酮类阈值一般高于醛类,但性质较稳定、香气持久,又具有独特的清香、水果香、花香和甜味等令人愉悦的风味,不仅是样品中其它杂环类风味物质的中间体,而且对风味物质的构成起到一定的协同作用[14]

    图  4  蒸煮与未蒸煮样品中酮类物质的气味指纹图谱
    Figure  4.  Major volatile components of ketones in cooked and uncooked samples

    通常酯类由低级饱和脂肪酸与醇类化合而成,具有令人愉快的水果香气或酒香味。未蒸煮样品中的酯类物质共检出6种,蒸煮样品中检出3种。如图5所示,酶解辅助预糊化技术对糙米酯类物质的含量影响较大,且蒸煮之后的CHLXM的酯类物质种类减少。酶解辅助预糊化技术和蒸煮处理中都有加热这一工序,所以推测酯类物质受热易分解。总之,酯类物质在CHLXM和CHLXM-ZHS中的贡献较小。

    图  5  蒸煮与未蒸煮样品中酯类物质的气味指纹图谱
    Figure  5.  Major volatile components of esters in cooked and uncooked samples

    大数的支链醇来自氨基酸的Strecker的降解,而多数的直链醇来自不饱和脂肪酸的氧化降解[22]。如图6所示,2-己醇在CHLXM中含量最高,而CHLXM-ZHS中含量较高的醇类物质有甲醇和2-己烯-1-醇。这也说明热处理有利于醇类风味物质的形成[18],其中,不饱和醇比饱和醇的气味强烈,且阈值较低。而酶解辅助预糊化糙米经蒸煮之后2-己烯-1-醇的含量增加,以致果香味增强。饱和醇中:C1~C3香气轻快;C7~C10芳香气味;C10以上,气味逐渐减弱至无味。所以醇类物质对样品整体风味的和谐有一定的贡献,有助于提升整体风味。

    图  6  蒸煮与未蒸煮样品中醇类物质的气味指纹图谱
    Figure  6.  Major volatile components of alcohols in cooked and uncooked samples

    图7所示,与CHLX和CK相比,CHLXM中产生的酸类挥发性物质含量较高,而且CHLXM-ZHS中的酸类物质消失。酸类物质阈值较高,对样品香味的直接贡献不大,但具有调节体系口感和维持酯类物质香气的作用[17]。此外,本实验还检出少量的杂环类和烷烃类等化合物,一般认为其存在不仅对产品整体风味的和谐有一定的贡献,而且也是形成杂环化合物的重要中间体,有助于提升整体风味。

    图  7  蒸煮与未蒸煮样品中其它类物质的气味指纹图谱
    Figure  7.  Major volatile components of other substances in cooked and uncooked samples

    基于以上分析可知,在糙米经酶解辅助预糊化处理过程中,相较于糙米(CHLX)和白米(CK),挥发性物质的种类并没有减少。其中含量明显增加的风味物质有苯乙醛、戊醛、3-甲基丁醛和酸类物质,明显减少的风味物质的含量与白米(CK)相近,如乙酸己酯、蘑菇醇、正戊基呋喃、戊醛、辛醛与庚醛等。而CHLXM-ZHS与CHLXM相比,风味物质的种类减少了13种,包括辛烯醛、苯乙醛、苯乙酮、γ-丁内酯、二甲苯、戊酸、己酸和丁酸等,风味物质含量增加的有正戊基呋喃、2-己烯-1-醇、丙酮、己酮、辛酮和庚醛等。风味物质含量增加的主要原因是前体物质含量随着煮制温度升高而增加直到达到了风味物质产生所需条件而被分解氧化[23]。而含量下降是风味前体物质已基本被反应等多种可能因素导致的。在酶解辅助预糊化处理和煮制过程中,酮类,醇类和酸类物质含量均有所变化,这是因为部分酮、醇、酸类是形成杂环类化合物或其他风味物质的重要中间体[19],并对其它风味化合物起到一定的协同作用。在酶解辅助预糊化技术和煮制过程中,既有美拉德反应及与其相关的Strecker降解反应,也有脂肪的氧化降解反应。如图3图7所示,糙米的关键风味物质主要是醛类物质,其次是醇类物质和酮类物质,崔琳琳等[24]也给出过相似的研究结果。Champagne等[25]究发现,米饭风味主要以醛类、醇类、酮类及呋喃类物质为主。Zeng等[26]的研究证实,在其鉴定出的45种米饭风味化合物中,醛类所占比重最大,其次是醇类及酮类物质。

    用SPSS软件对的挥发性物质进行主成分分析(PCA)。主成分分析是设法将原来众多具有一定相关性的几个指标,重新组合成一组新的互相无关的综合指标来代替原来的指标,以评价原来变量的信息。在PCA主成分分析图中,距离越远的两组差异越大,距离越近的两组差异越小。由表6可知第一主成分的贡献率为91.928%,第二主成分的贡献率为4.998%,累计贡献率为96.926%,累计贡献率大于80%可以进行PCA分析[27]。利用Origin软件制作PCA图,结果如图8所示。可以看出蒸煮前后的白米、糙米和酶解辅助预糊化糙米的主成分有明显差异,可将白米、糙米、酶解辅助预糊化糙米及蒸煮后的白米、糙米、酶解辅助预糊化糙米区分开。

    表  6  解释的总方差
    Table  6.  Total variance explained
    成份初始特征值 提取平方和载入
    合计方差的(%)累积(%)合计方差的(%)累积(%)
    15.51691.92891.928 5.51691.92891.928
    20.34.99896.926
    30.132.1699.086
    40.0550.914100
    51.01E-111.68E-10100
    65.45E-129.08E-11100
    下载: 导出CSV 
    | 显示表格
    图  8  样品挥发性风味物质的PCA分析图
    注:CK:白米;CHLK:糙米;CHLKM:酶解辅助预糊化糙米;CK-ZHS:蒸煮后的白米;CHLK-ZHS:蒸煮后的糙米;CHLKM-ZHS:蒸煮后的酶解辅助预糊化糙米。
    Figure  8.  PCA analysis of the sample with uncooked volatile flavor

    利用响应面试验对糙米的酶解辅助预糊化条件进行优化,得到最佳酶解辅助预糊化工艺条件:酶解时间为4 h、酶的添加量为0.1%、预糊化时间为10 min和酶解温度为39 ℃,在此条件下生产的糙米蒸煮时间短且糠味不明显。用气相色谱–离子迁移谱法(GC-IMS)测定白米、糙米和酶解辅助预糊化糙米及其蒸煮前后的挥发性物质,在酶解辅助预糊化糙米中检测到56种挥发性物质,其中醛类19种,酮类8种,酯类6种,醇类12种及其它类11种;而蒸煮后的酶解辅助预糊化糙米中检测到44种挥发性物质,其中醛类18种,酮类8种,酯类3种,醇类11种及其它类4种。糙米酶解辅助预糊化前、后挥发性风味成分改变最大的为醛类、酮类和醇类;而酶解辅助预糊化糙米蒸煮前、后挥发性风味成分改变最大的为酯类和其它类中的酸类,这是因为酶解辅助预糊化和蒸煮过程中,糙米中的糖类、核苷酸、自由氨基酸、肽类等发生美拉德反应[28-29],从而产生醛类、呋喃衍生物、酮类、吡嗪类等挥发性风味物质[30-31]

  • 图  1   酶解时间(A)、酶的添加量(B)、预糊化时间(C)、酶解温度(D)对糙米蒸煮时间的影响

    注:小写字母不同表示差异显著,P<0.05。

    Figure  1.   Effect of processing time(A), enzyme addition(B), steaming time(C) and temperature(D) on the cooking time of brown rice

    图  2   各因素交互作用对糙米蒸煮时间影响的响应面图

    Figure  2.   Response surface plots for the effects of different factors onthe cooking time of brown rice

    图  3   蒸煮与未蒸煮样品中醛类物质的气味指纹图谱

    Figure  3.   Major volatile components of aldehydes in cooked and uncooked samples

    图  4   蒸煮与未蒸煮样品中酮类物质的气味指纹图谱

    Figure  4.   Major volatile components of ketones in cooked and uncooked samples

    图  5   蒸煮与未蒸煮样品中酯类物质的气味指纹图谱

    Figure  5.   Major volatile components of esters in cooked and uncooked samples

    图  6   蒸煮与未蒸煮样品中醇类物质的气味指纹图谱

    Figure  6.   Major volatile components of alcohols in cooked and uncooked samples

    图  7   蒸煮与未蒸煮样品中其它类物质的气味指纹图谱

    Figure  7.   Major volatile components of other substances in cooked and uncooked samples

    图  8   样品挥发性风味物质的PCA分析图

    注:CK:白米;CHLK:糙米;CHLKM:酶解辅助预糊化糙米;CK-ZHS:蒸煮后的白米;CHLK-ZHS:蒸煮后的糙米;CHLKM-ZHS:蒸煮后的酶解辅助预糊化糙米。

    Figure  8.   PCA analysis of the sample with uncooked volatile flavor

    表  1   响应面试验因素和水平

    Table  1   Factors and levels of Box-Behnken experiment

    水平因素
    A酶解时间
    (h)
    B酶的添加量
    (%)
    C预糊化时间
    (min)
    D酶解温度
    (℃)
    −130.05630
    040.10940
    150.151250
    下载: 导出CSV

    表  2   采集数据程序设置

    Table  2   Collect data application settings

    TimeE1(mL/min)E2(mL/min)
    00:00,001502
    02:00,001502
    10:00,0015010
    20:00,00150100
    25:00,00150130
    30:00,00150130
    下载: 导出CSV

    表  3   中心组合实验设计及结果

    Table  3   Central composite design and results

    实验号A-酶解时间B-酶的添加量C-预糊化时间D-酶解温度Y-蒸煮时间(min)
    1−10−1026.28
    200−1126.20
    3000021.57
    4−101024.57
    50−10127.24
    601−1025.97
    7100126.23
    8010−125.29
    9−1−10027.26
    10000021.02
    111−10024.95
    120−10−126.53
    13100−124.23
    140−11024.80
    15−110026.20
    16000021.26
    17110024.65
    18011025.64
    1900−1−125.69
    20010125.83
    21−100126.25
    22000021.34
    23101024.10
    240−1−1027.62
    25000022.56
    26001124.96
    2710−1025.67
    28001−125.25
    29−100−125.34
    下载: 导出CSV

    表  4   回归方程方差分析表

    Table  4   Analysis of variance (ANOVA) of regression equations

    方差来源平方和df均方 F P值
    模型88.33146.3119.58<0.0001**
    A-酶解时间3.0713.079.530.0080**
    B-酶的添加量1.9311.936.010.0280*
    C-预糊化时间5.4815.4817.010.0010**
    D-酶解温度1.6011.604.960.0428*
    AB0.1410.140.450.5141
    AC4.900E-00314.900E-0030.0150.9037
    AD0.3010.300.920.3533
    BC1.5511.554.810.0457*
    BD7.225E-00317.225E-0030.0220.8831
    CD0.1610.160.500.4926
    A219.71119.7161.07<0.0001**
    B241.34141.34128.31<0.0001**
    C222.70122.7070.46<0.0001**
    D230.18130.1893.68<0.0001**
    残差4.51140.32
    失拟项3.08100.310.860.6157
    纯误差1.4340.36
    总和92.8428
    注:*差异显著,P<0.05;**差异极显著,P<0.01。
    下载: 导出CSV

    表  5   不同样品挥发性物质种类及峰强度信息

    Table  5   Volatile substance types and peak strength information of different samples

    编号化合物CAS号保留指数气味特征保留时间(s)迁移时间(ms)
    醛类
    1苯甲醛二聚体100-52-7958.8坚果味、苦味434.6551.467
    2苯甲醛单体100-52-7958.8坚果味、苦味434.6551.148
    3己醛二聚体75-07-0793.6青草香251.161.563
    4己醛单体75-07-0797.2青草香254.0851.253
    5糠醛98-01-1958.8甜杏仁香烤味434.6551.148
    6戊醛二聚体110-62-3698.8木香、水果香182.3251.182
    7戊醛单体110-62-3902.9木香、水果香356.6551.331
    8辛醛二聚体124-13-01006青草、柑橘味516.9451.821
    9辛醛单体124-13-01006.6青草、柑橘味518.1151.404
    10庚醛二聚体111-71-7901.3水果香354.7051.697
    11庚醛单体111-71-7797.2水果香254.0851.253
    12正壬醛二聚体124-19-61105.5柑橘、肥皂味729.8851.943
    13正壬醛单体124-19-61107.4柑橘、肥皂味733.981.475
    143-甲基丁醛二聚体590-86-3658.3麦芽味162.4351.404
    153-甲基丁醛单体590-86-3656.3麦芽味161.6551.172
    162-辛烯醛(反式)单体2548-87-01063.7青草、脂肪味636.871.331
    17反式2-庚烯醛18829-55-5958.8辛辣、脂肪味434.6551.666
    18癸醛112-31-21209柑橘香960.961.538
    192-辛烯醛(反式)2548-87-01063.1青草、脂肪味635.5051.818
    20苯乙醛122-78-11043.8青草香味593.7751.253
    酮类
    21苯乙酮98-86-21057.8香辣杏仁味6241.187
    222-3-丁二酮431-03-8625.1清香味149.9551.168
    232-丁酮二聚体78-93-3587.2类似薄荷香味137.4751.243
    242-丁酮单体78-93-3589.6类似薄荷香味138.2551.057
    25丙酮67-64-1504类似薄荷香味110.761.119
    26己酮591-78-6783.4果味、黄油味242.971.188
    272-庚酮二聚体110-43-0893梨香344.9551.631
    282-庚酮单体110-43-0892.9梨香344.761.261
    292-戊酮107-87-9690.4类似薄荷香味177.6451.12
    酯类
    30乙酸乙酯二聚体141-78-6614.2水果香味146.251.335
    31乙酸乙酯单体141-78-6616水果香味146.8351.093
    32二氢草莓酸乙酯108-64-5831.5果香味283.141.257
    33乙酸己酯142-92-7987.4果香,似菠萝香482.6251.414
    34γ-丁内酯96-48-0919.6奶油焦糖味377.911.081
    35乙偶姻513-86-0719.2甜奶油味1951.055
    醇类
    36庚醇单体53535-33-4977.8甜香、坚果味465.8551.397
    37正己醇二聚体111-27-3873.9果香味323.8951.644
    38正己醇单体111-27-3877.7果香味327.991.327
    49正戊醇单体71-41-0768.4水果味231.0751.25
    40正戊醇二聚体71-41-0766.9水果味229.9051.51
    41甲醇二聚体3857-25-8958.9甜焦糖味434.851.565
    42甲醇单体3857-25-8961.1甜焦糖味438.361.254
    43反式2-己烯-1-醇二聚体2305-21-7849.6果香299.7151.516
    44反式2-己烯-1-醇单体2305-21-7849.6果香299.7151.516
    452-己醇626-93-7820.7果香味273.781.278
    46蘑菇醇3391-86-4984.5蘑菇香477.5551.155
    472-辛醇6169-06-8948果香味418.081.443
    其它类
    48正基呋喃3777-69-3993.7果香味494.131.253
    491,8-桉叶素单体470-82-61031.4樟树脑味567.841.295
    50柠檬烯138-86-31029.7柠檬香564.331.218
    51α-蒎烯80-56-8974.1松木香459.6151.216
    521,8-桉叶素二聚体470-82-61031.4樟树脑味567.841.729
    53二甲苯95-47-6923天竺葵382.3951.064
    54戊酸二聚体109-52-4885.9酸味336.961.509
    55戊酸单体109-52-4888.4酸味339.691.233
    56己酸142-62-1888.4干酪酸味339.691.233
    57丁酸107-92-6777.8奶酪酸味238.4851.162
    58己二烯二硫化2179-57-91076.2大蒜味664.3651.196
    下载: 导出CSV

    表  6   解释的总方差

    Table  6   Total variance explained

    成份初始特征值 提取平方和载入
    合计方差的(%)累积(%)合计方差的(%)累积(%)
    15.51691.92891.928 5.51691.92891.928
    20.34.99896.926
    30.132.1699.086
    40.0550.914100
    51.01E-111.68E-10100
    65.45E-129.08E-11100
    下载: 导出CSV
  • [1] 谭斌. 全谷物食品: 健康食品新趋势[J]. 农产品加工,2010(4):4−5.
    [2]

    Deng Y, Zhong Y, Yu W, et al. Effect of hydrostatic high pressure pretreatment on flavor volatile profile of cooked rice[J]. Journal of Cereal Science,2013,58(3):479−487.

    [3]

    Wu P, Deng R, Wu X, et al. In vitro gastric digestion of cooked white and brown rice using a dynamic rat stomach model[J]. Food Chemistry,2017,237:1065−1072. doi: 10.1016/j.foodchem.2017.05.081

    [4] 昝学梅, 刘明, 孟宁, 等. 健康谷物及其制品在乳品加工中的应用[J]. 食品与机械,2019,35(11):5−10.
    [5] 蒋龙伟. 糙米发芽前复合酶预处理工艺优化研究[D]. 哈尔滨: 东北农业大学, 2017.
    [6]

    Das M, Banerjee R, Bal S. Evaluation of physicochemical properties of enzyme treated brown rice (Part B)[J]. LWT-Food Science & Technology,2008,41(10):2092−2096.

    [7] 王宏伟, 王新天, 肖乃勇, 等. 预糊化淀粉对糯米粉糊化、流变性能及微观结构的影响[J]. 中国粮油学报,2019,34(4):57−62.
    [8]

    Bortnowska G Y, Balejko J, Schube V, et al. Stability and physicochemical properties of model salad dressings prepared with pregelatinized potato starch[J]. Carbohydr Polym,2002,111:624−632.

    [9]

    Anastasiades A, Thanou S, Loulis D, et al. Rheological and physical characterization of pregelatinized maize starches[J]. Journal of Food Engineering,2016,52(1):57−66.

    [10]

    Hedayati S, Shahidi F, Koocheki A, et al. Comparing the effects of sucrose and glucose on functional properties of pregelatinized maize starch[J]. International Journal of Biological Macromolecules,2016:499−504.

    [11]

    Seetapan N, Fuongfuchat A, Gamonpilas C, et al. Effect of modified tapioca starch and xanthan gum on low temperature texture stability and dough viscoelasticity of a starch-based food gel[J]. Journal of Food Engineering,2013,119(3):446−453. doi: 10.1016/j.jfoodeng.2013.06.010

    [12]

    Ortolan F, Brites L T G, Montenegro F M, et al. Effect of extruded wheat flour and pre-gelatinized cassava starch on process and quality parameters of French-type bread elaborated from frozen dough[J]. Food Research International,2015,76(Pt 3):402−409.

    [13]

    Fang D, Yang W, Kimatu B M, et al. Comparison of flavour qualities of mushrooms (Flammulina velutipes) packed with different packaging materials[J]. Food Chemistry,2017,232(OCT. 1):1−9.

    [14]

    Giri A, Osako K, Ohshima T. Identification and characterisation of headspace volatiles of fish miso, a Japanese fish meat based fermented paste, with special emphasis on effect of fish species and meat washing[J]. Food Chemistry,2009,120(2):621−631.

    [15] 刘明, 孟宁, 朱运恒, 等. 低温等离子体技术改善糙米蒸煮品质工艺优化及热力学特性研究[J]. 粮油食品科技,2020,28(2):49−54.
    [16] 孔祺, 李星骆, 刘庆庆. 蒸汽预糊化处理对发芽糙米结构及吸水特性的影响[J]. 食品科技,2019,44(12):193−198.
    [17]

    Bryant R J, Mcclung A M. Volatile profiles of aromatic and non-aromatic rice cultivars using SPME/GC–MS[J]. Food Chemistry,2011,124(2):501−513. doi: 10.1016/j.foodchem.2010.06.061

    [18]

    Shurong L, Li W, Chunhong Z, et al. Analysis of the key odorants of roasted peanut[J]. Scientia Agricultura Sinica,2010,43(15):3199−3203.

    [19]

    Zhao J, Wang M, Xie J, et al. Volatile flavor constituents in the pork broth of black-pig[J]. Food Chemistry,2017,226(JUL. 1):51−60.

    [20]

    Yaylayan V A, Haffenden L J W. Mechanism of pyrazole formation in [~(13)C-2] labeled glycine model systems: N-N bond formation during Maillard reaction[J]. Food Research International,2003,36(6):571−577. doi: 10.1016/S0963-9969(03)00003-6

    [21] 朱琪, 师希雄, 韩玲, 等. 祁连清泉清炖羊肉汤挥发性化合物的检测[J]. 食品与发酵科技,2017,53(4):109−114.
    [22]

    Garcia-Gonzalez D L, Tena N, Aparicio-Ruiz R, et al. Relationship between sensory attributes and volatile compounds qualifying dry-cured hams[J]. Meat Science,2008,80(2):315−325. doi: 10.1016/j.meatsci.2007.12.015

    [23] 单启梅, 赵晓策, 罗瑞明, 等. 滩羊肌肉在煮制过程中可挥发性化合物的变化[J]. 食品科学: 1-10 [2020-07-20]. http://kns.cnki.net/kcms/detail/11.2206.TS.20200414.1325.002.html.
    [24] 崔琳琳, 赵燊, 周一鸣, 等. 基于GC-MS和电子鼻技术的大米挥发性风味成分分析[J]. 中国粮油学报,2018,33(12):134−141. doi: 10.3969/j.issn.1003-0174.2018.12.021
    [25]

    Champagne E T, Bettgarber K L, Thomson J L, et al. Impact of presoaking on flavor of cooked rice[J]. Cereal Chemistry,2008,85(3):275−276.

    [26]

    Zeng Z, Zhang H, Chen J Y, et al. Direct extraction of volatiles of rice during cooking using solid-phase microextraction[J]. Cereal Chemistry,2007,84(5):423−427. doi: 10.1094/CCHEM-84-5-0423

    [27]

    Lin, Lu, Shaoping, et al. Classification of rice by combining electronic tongue and nose[J]. Food Analytical Methods,2015,8(8):1839−1902.

    [28] 祁岩龙, 冯怀章, 于洋, 等. 美拉德反应研究进展及在食品工业中的应用[J]. 食品工业,2018,39(3):248−252.
    [29] 邓仕彬. 美拉德反应中间体的水相制备及其在卷烟中的应用[D]. 无锡: 江南大学, 2017.
    [30] 苑蘅. Heyns化合物的质谱碎裂规律及其在Maillard反应中含量与变化研究[D]. 北京: 中国农业大学, 2016.
    [31]

    Xie J, Sun B, Zheng F, et al. Volatile flavor constituents in roasted pork of Mini-pig[J]. Food Chemistry,2008,109(3):506−514. doi: 10.1016/j.foodchem.2007.12.074

图(8)  /  表(6)
计量
  • 文章访问数:  284
  • HTML全文浏览量:  74
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-08
  • 网络出版日期:  2021-01-27
  • 刊出日期:  2021-04-14

目录

/

返回文章
返回