酱香型白酒是各大香型白酒中酿造工艺、微生物体系和风味结构最为复杂的,酱香白酒的不同风格特征与不同类型的基酒勾调密切相关,如不同产地基酒中各类风味化合物含量与种类有较大差异。据初步统计,贵州赤水河流域白酒企业有400家,白酒产能达到40万kL,其中酱香型白酒产能约30万kL,白酒中的风味成分会受到产区因素的影响,即使是同种工艺类型的白酒,产区环境和气候条件都会造成白酒风味的差异[1]。贵州各地具有独特的土壤、水分、地势条件等一系列自然条件,因此不同流域酿造出来的白酒具有不同的口感和风味。对赤水河流域各大酱酒产业不同空间分布及特征影响进行深入研究分析,对于建立赤水河流域不同产区酱香白酒特征模型有重要意义。
近年来,随着气相色谱(gas chromatography,GC)分析技术广泛应用于白酒香气物质的定性定量分析使得白酒质量、香型、品牌的鉴别结果更加准确[2-4]。统计学分析已经逐步渗入到食品科学领域,与仪器检测相结合展示了较高的优越性。主成分分析(principal component analysis,PCA)是一种比较常用的多元统计分析方法,常用于酒类、茶叶和食品等领域[5-8],正交偏最小二乘-判别分析(orthogonal partial least-squares-discrimination analysis,OPLS-DA)是在偏最小二乘-判别分析(partial least-squares-discrimination analysis,PLS-DA)的基础上,进行了正交变换的矫正,判别效果及主成分得分图的可视化效果更加明显。李爱兰等[9]使用吹扫捕集-气质联用方法鉴别不同产地、品牌浓香型白酒。基于浓香型白酒的挥发性风味成分的差异,通过OPLS-DA可对其产地进行准确区分,通过变量投影重要性(variable importance for the projection,VIP)值分析筛选出辛酸乙酯、己酸己酯等40种潜在差异性化合物。程平言等[10]对于白酒原产地的鉴定通过偏最小二乘-判别分析法和逐步线性判别分析法筛选出27个重要特征离子,然后用筛选出的27个特征离子构建反向传播(back propagation,BP)神经网络模型和支持向量机(support vector machine,SVM)模型,其原产地预测准确率分别达96.2%和97.7%。
本研究以赤水河不同产区的酱香型白酒为研究对象,采用气相色谱法对酒样中的64种挥发性风味物质进行定量分析及差异分析,采用正交偏最小二乘-判别分析(OPLS-DA)构建模型,并利用变量投影重要性(VIP)值筛选差异挥发性风味化合物,并进行聚类分析(cluster analysis,CA)。以期对酱香型地域划分、企业品牌标准体系建立提供研究基础与理论指导,对酱香产区分类有一定指导意义。
51个酱香白酒样品,其中包括16个MT产区系列酒样(MT1~MT16),9个ZJ产区系列酒样(ZJ1~ZJ9),11个XJ产区系列酒样(XJ1~XJ11),10个JS产区系列(JS1~JS10)酒样,5个LJ产区系列酒样(LJ1~LJ5),均由市场采购获得。
叔戊醇、乙酸正戊酯、2-乙基丁酸(纯度均>98%)、乙酸乙酯、丙酸乙酯、异丁酸乙酯等35个标准品(纯度≥97.0%):上海阿拉丁生化科技股份有限公司;无水乙醇(色谱纯):天津市科密欧化学试剂有限公司;乙醛(纯度≥97.0%):美国Sigma-Aldrich公司;2-甲基丁酸乙酯(纯度≥97.0%):上海贤鼎生物科技有限公司;异戊酸乙酯、丁二酸二乙酯、邻苯二甲酸二丁酯等标准品(纯度均≥97.0%):上海易恩化学技术有限公司;庚酸乙酯、壬酸乙酯、异丁醇、2,3-丁二醇、异戊醛(纯度均≥97.0%):上海麦克林生化科技有限公司。
7890A气相色谱仪:安捷伦科技有限公司;PAL多功能自动进样器:瑞士斯特分析仪器有限公司;Millipore-Q超纯水系统、2 mL进样瓶:佰奥萃(天津)生物科技有限公司。
1.3.1 酒样预处理
移取5 mL酒样过0.22 μL有机滤膜,吸取990 μL过膜后的酒样于2 mL气相进样瓶中,加入体积分数1%的叔戊醇(内标)溶液10 μL,加盖密封,备用,每个样品重复2~3次。
1.3.2 挥发性风味物质检测
采用气相色谱法检测白酒样品中的挥发性风味物质。
GC条件:SH-Rtx-Wax色谱柱(30 m×0.25 mm×0.25 μm);载气为氮气(N2)(99.999%);进样口温度250 ℃;分流比30∶1;检测器温度250 ℃;氢气流量30 mL/min,空气流量300mL/min,升温程序为初始温度30℃,保持3min,以3℃/min升温至180 ℃,再以15 ℃/min升温至210 ℃,保持8 min。
定性、定量分析:采用保留时间定性,根据张晓婕等[11]方法对白酒风味进行定量,以待测物与内标物的含量比(X)为横坐标,峰面积比(Y)为纵坐标,采用内标法定量酒样中各风味物质。绘制各挥发性风味物质的标准曲线,获得其标准曲线回归方程。
香气活度值(odor active value,OAV)分析:通过查阅相关文献确定每个物质阈值计算得到OAV,一般认为OAV>1则代表该物质具有贡献作用,计算公式如下:
式中:C为各挥发性物质质量浓度,mg/L;OT为各挥发性物质阈值,μg/L。
1.3.3 正交偏最小二乘-判别分析
从146个观测值中筛选掉异常值后,从每个产区随机选取3个样本作为测试集,共计15个,其余样本为训练集,利用正交偏最小二乘-判别分析构建不同产区酱香白酒鉴别模型。并利用VIP值筛选不同流域酱香白酒产区的挥发性差异物,根据所选差异物绘制层次聚类热图,直观表达风味化合物的潜在联系。
1.3.4 数据处理
采用SIMCA软件13.0进行主成分分析以及模型构建。采用MetaboAnalyst2.0绘制聚类分析热图。
不同产区各挥发性风味化合物相对含量见表1,不同产区酱香型白酒挥发性风味物质测定结果及OAV见表2。由表1和表2可知,不同酒样中共检出挥发性风味化合物64种,其中,酯类19种,醇类13种,酸类8种,醛酮类7种,吡嗪类2种,呋喃类5种,芳香族类8种,酚类2种,不同产区酒样的挥发性风味物质总含量存在差异,ZJ产区酒样挥发性风味物质总含量略偏低,为9 125.45 mg/L,其余酒样总含量较为接近,为11 010.40~11 716.49 mg/L,这与李利利[14]的研究结果一致。不同产区酱香白酒风味物质主要为酯类物质,其次是醇类、酸类、醛酮类,吡嗪类物质在白酒中含量极低,但赋予白酒坚果和焙烤香气,对其他化合物有明显的烘香作用[11]。酯类物质在白酒中能赋予其花香和果香[15],对白酒空杯后味具有贡献,酸类化合物是重要的呈味物质,能够调和酒体口味、稳定香气[15],醇类物质沸点低易挥发,对白酒风味层次具有贡献作用[16]。XJ和MT产区酒样中的酯类物质和酸类物质含量较高,且酯类物质中的乙酸乙酯和酸类物质中的乙酸在酱香型白酒样品中的含量最高。JS产区酒样的酯类物质含量最低,但醇类、呋喃类化合物含量最高,LJ产区酒样的酚类化合物和醛酮类化合物明显高于其他产区,MT产区酒样与LJ产区酒样相反,ZJ产区酒样中的芳香族类物质相对含量最高。
表1 不同产区酒样各类挥发性风味化合物相对含量
Table 1 Relative contents of various types of volatile flavor compounds in Baijiu samples from different producing areas
酒样编号酯类/%醇类/%酸类/%醛酮类/%吡嗪类/%呋喃类/%芳香族类/%酚类/%XJ ZJ MT JS LJ 40.31 37.42 34.89 25.02 35.28 21.85 21.36 32.62 38.67 21.95 18.02 17.47 15.03 12.73 14.17 15.91 18.30 13.80 17.84 24.12 0.02 0.14 0.21 0.04 0.02 1.74 1.90 1.84 2.26 1.67 1.97 3.17 1.52 3.16 2.31 0.17 0.23 0.08 0.28 0.50
表2 不同产区酱香型白酒挥发性风味物质测定结果及香气活度值
Table 2 Determination results of volatile flavor compounds and odor active value of sauce-flavor Baijiu from different producing areas
化合物酯类乙酸乙酯丙酸乙酯异丁酸乙酯丁酸乙酯2-甲基丁酸乙酯异戊酸乙酯乙酸异戊酯戊酸乙酯己酸乙酯庚酸乙酯乳酸乙酯辛酸乙酯己酸异戊酯壬酸乙酯癸酸乙酯丁二酸二乙酯肉豆蔻酸乙酯棕榈酸乙酯邻苯二甲酸二丁酯总计醇类仲丁醇正丙醇异丁醇2-戊醇正丁醇异戊醇正戊醇2-庚醇2-壬醇含量/(mg·L-1)XJ ZJ MT JS LJ阈值/(μg·L-1)[12-13]OAV XJ ZJ MT JS LJ 2 648.65±301.03 62.01±13.33 43.89±16.26 35.09±6.50 4.05±1.46 9.48±3.87 4.60±0.81 6.14±1.50 26.90±8.36 1.62±0.30 1 545.24±132.00 3.47±0.53 8.28±2.66 13.01±0.85 4.75±0.48 9.15±0.39 7.10±0.42 37.66±10.32 26.88±6.66 4 497.98±342.51 1 981.93±250.64 69.29±20.26 38.82±13.71 55.03±10.93 4.17±1.18 8.74±2.34 4.02±0.32 11.62±3.61 29.32±8.53 2.09±0.31 1 089.84±127.18 2.89±0.22 7.14±1.89 10.99±0.84 4.21±0.14 7.74±2.74 7.67±0.70 37.86±18.25 41.65±8.92 3 415.03±275.16 2 324.21±319.53 71.60±6.58 24.35±6.39 26.13±4.3 9.05±1.88 12.48±3.44 4.37±0.44 3.33±1.00 7.90±3.87 1.09±0.45 1 409.48±94.06 2.18±0.23 8.48±1.75 12.08±0.28 4.09±0.1 8.82±0.15 7.00±0.38 28.29±7.47 56.97±6.32 4 021.9±367.59 1 776.31±320.83 60.64±9.86 23.58±9.18 76.23±12.5 7.51±5.54 11.61±4.3 4.79±1.07 12.78±3.19 47.38±17.13 2.70±0.59 790.65±160.14 4.12±0.83 12.13±7.89 13.46±4.86 4.29±0.64 2.40±2.06 7.85±2.01 27.55±8.74 46.04±22.08 2 932.03±448 1 720.13±248.51 48.92±9.16 27.56±12.73 44.87±3.20 6.87±2.15 19.17±2.24 3.30±0.32 13.56±0.94 88.28±8.39 5.02±0.40 1 801.23±99.66 7.67±0.54 6.73±2.02 14.44±0.82 4.26±0.07 10.25±3.02 13.78±1.12 27.50±9.3 20.46±6.88 3 884±249.08 32 551.6 19 019.33 57.47 81 18 6.89 93.93 268 553 13 153.17 128 083.8 127 1 400 3 150.61 1 122.3 353 193.25 81.37±9.25 3.26±0.70 763.63±282.95 430.57±79.71 224.78±81.00 1 376.51±560.96 49.02±8.58 229.44±56.17 486.12±151.11 0.12±0.02 12.06±1.03 269.74±40.91 5.92±1.90 4.13±0.27 4.23±0.43 0.03±0 60.89±7.70 3.64±1.07 675.52±238.54 675.25±134.06 231.79±65.82 1 268.24±340.22 42.76±3.36 434.01±134.73 529.99±154.1 0.16±0.02 8.51±0.99 224.47±17.38 5.10±1.35 3.49±0.27 3.76±0.12 0.02±0.01 71.40±9.82 3.76±0.35 423.62±111.23 320.64±52.79 502.86±104.31 1 811.68±499.14 46.56±4.73 124.35±37.37 142.76±69.86 0.08±0.03 11.00±0.73 169.04±17.81 6.06±1.25 3.83±0.09 3.65±0.09 0.02±0 54.57±9.86 3.19±0.52 410.38±159.66 935.4±153.4 417.29±307.53 1 684.37±624.30 50.99±11.4 477.32±118.98 856.27±309.65 0.21±0.04 6.17±1.25 320.16±64.21 8.67±5.64 4.27±1.54 3.83±0.57 0.01±0.01 52.84±7.63 2.57±0.48 479.5±221.53 550.61±39.21 381.39±119.51 2 782.70±324.76 35.14±3.43 506.4±35.27 1 595.48±151.65 0.38±0.03 14.06±0.78 596.23±42.28 4.81±1.44 4.58±0.26 3.80±0.06 0.03±0.01-- - --- - --- - --- - --- - --- - -80.01±13.17 1 749.35±536.93 138.00±9.88 2.25±1.58 73.91±11.48 305.73±30.6 3.69±0.92 1.09±0.32 15.89±4.08 51.03±13.98 1 207.77±209.63 158.30±37.46 3.69±3.30 76.54±10.17 351.35±24.66 5.80±1.65 1.38±0.36 19.68±5.02 124.27±28.72 2 971.9±856.63 160.86±14.57 1.54±0.46 51.18±9.22 358.25±19.76 3.22±0.57 0.79±0.14 23.35±2.65 209.99±93.75 3 510.63±774.08 169.34±22.38 4.52±1.63 90.23±18.9 411.26±90.71 8.99±3.36 2.11±0.67 40.87±5.99 35.13±3.19 1 543.50±228.43 164.25±49.42 9.75±2.02 96.12±54.35 451.8±26.42 5.95±0.43 1.64±1.16 31.48±4.15 50 000 53 952.63 28 300-2 733.35 179 190.83 4 000 1 433.94 75 1.60±0.26 32.42±9.95 4.88±0.35-27.04±4.2 1.71±0.17 0.92±0.23 0.76±0.22 211.92±54.33 1.02±0.28 22.39±3.89 5.59±1.32-28.00±3.72 1.96±0.14 1.45±0.41 0.97±0.25 262.39±66.92 2.49±0.57 55.08±15.88 5.68±0.51-18.72±3.37 2.00±0.11 0.81±0.14 0.55±0.10 311.37±35.32 4.20±1.87 65.07±14.35 5.98±0.79-33.01±6.92 2.30±0.51 2.25±0.84 1.47±0.47 544.88±79.91 0.70±0.06 28.61±4.23 5.80±1.75-35.17±19.88 2.52±0.15 1.49±0.11 1.14±0.81 419.79±55.30
续表
化合物2,3-丁二醇辛醇1,2-丙二醇壬醇总计酸类乙酸丙酸异丁酸丁酸异戊酸戊酸己酸癸酸总计醛酮类乙醛异丁醛异戊醛乙缩醛己醛乙偶姻2-壬酮总计吡嗪类2,3,5-三甲基吡嗪2,3,5,6-四甲基吡嗪总计呋喃类2-丁基呋喃糠醛2-乙酰基呋喃5-甲基-2-乙酰基呋喃糠醇总计芳香族类苯甲醛苯乙醛苯乙酮苯甲酸乙酯苯乙酸乙酯乙酸苯乙酯苯丙酸乙酯β-苯乙醇总计酚类愈创木酚4-乙基愈创木酚总计合计含量/(mg·L-1)XJ ZJ MT JS LJ阈值/(μg·L-1)[12-13]OAV XJ ZJ MT JS LJ 43.48±2.13 1.87±0.20 3.84±0.63 19.34±2.55 2 438.44±556.52 49.78±4.92 1.85±0.99 4.79±0.85 17.58±3.40 1 949.53±201.12 35.48±2.21 1.83±0.15 4.11±0.45 23.65±2.99 3 760.43±895.37 59.71±13.71 0.74±0.85 5.29±1.50 17.56±2.48 4 531.24±854.03 51.54±2.94 2.19±0.12 4.57±0.58 18.87±1.42 2 416.80±223.87——--1 100-50-1.70±0.18-386.78±51.00-1.69±0.90-351.52±67.93-1.67±0.14-472.96±59.75-0.67±0.77-351.15±49.59-1.99±0.11-377.42±28.32-1 778.78±125.68 124.44±12.76 15.66±2.09 22.69±2.48 22.67±3.42 5.34±0.80 20.99±5.13 20.04±3.26 2 010.61±130.12 1 348.30±118.41 106.18±16.7 16.88±4.12 32.33±5.99 20.12±4.76 8.86±2.25 24.84±6.99 36.74±6.27 1 594.25±114.54 1 484.28±120.48 127.26±11.23 18.38±2.70 16.69±2.39 28.46±4.21 4.00±0.86 7.74±3.19 45.95±2.91 1 732.74±112.20 1 224.63±80.15 105.81±10.51 13.07±3.48 42.06±7.35 18.22±5.47 10.08±1.17 39.76±10.85 37.39±12.39 1 491.01±73.81 1 292.27±113.87 89.61±6.98 15.52±0.77 28.17±1.80 21.93±2.18 10.45±0.75 75.83±8.81 26.5±7.93 1 560.26±114.02 160 000 18 100 126 964.64 1 045.47 389.11 2 517.16 13 736.77-11.12±0.79 6.88±0.70 124.30±16.58 23.53±2.57 21.68±3.27 13.72±2.07 8.31±2.03 1.46±0.24-8.43±0.74 5.87±0.92 133.99±32.69 33.52±6.21 19.24±4.56 22.76±5.77 9.83±2.77 2.67±0.46-9.28±0.75 7.03±0.62 145.85±21.42 17.30±2.48 27.22±4.02 10.28±2.21 3.06±1.26 3.35±0.21-7.65±0.50 5.85±0.58 103.75±27.63 43.60±7.62 17.43±5.23 25.90±3.01 15.73±4.29 2.72±0.90-8.08±0.71 4.95±0.39 123.15±6.12 29.20±1.87 20.97±2.09 26.85±1.93 30.00±3.49 1.93±0.58-1 151.48±211.93 33.44±4.05 58.10±6.29 468.95±65.99 16.31±1.94 45.02±6.42 2.05±0.08 1 775.35±271.56 1 097.65±227.62 31.22±7.87 51.43±12.34 404.71±74.87 14.58±4.08 67.46±21.81 2.85±0.33 1 669.9±325.94 1 009.57±162.25 38.96±3.47 66.27±5.76 411.78±56.92 18.17±2.28 43.19±7.01 3.16±0.26 1 591.11±225.15 1 359.69±324.05 55.76±16.61 115.80±38.44 481.56±72.42 27.22±11.09 47.70±15.11 2.08±0.21 2 089.82±424.85 2 049.86±177.45 32.3±2.84 40.68±3.29 500.63±247.03 10.58±1.31 19.39±9.03 1.98±0.10 2 655.43±221.98 500-16.51 69 25.48 259 2 302.95±423.87-3 519.16±380.96 6 796.38±956.41 640.02±76.32 173.84±24.80 2 195.31±455.23-3 114.80±747.54 5865.36±1085.06 572.26±159.98 260.45±84.20 2 019.14±324.51-4 013.71±348.78 5 967.87±824.98 713.29±89.37 166.77±27.05 2 719.37±648.10-7013.68±2328.55 6979.18±1049.63 1 068.48±435.13 184.18±58.33 4 099.73±354.90-2 464.07±199.14 7255.51±3580.14 415.39±51.35 74.85±34.87——————0.84±0.81 1.62±1.08 2.45±0.86 2.11±0.89 10.80±7.32 12.91±7.78 2.74±0.31 21.42±5.05 24.16±5.31 0.27±0.23 3.91±2.14 4.18±2.19 1.02±0.57 0.89±1.79 1.91±2.03 729.86 80 073.16 1.15±1.11 0.02±0.01-2.89±1.22 0.13±0.09-3.76±0.42 0.27±0.06-0.37±0.31 0.05±0.03-1.39±0.78 0.01±0.02-4.52±0.53 164.29±40.71 3.60±0.93 3.08±0.16 18.87±4.09 194.36±43.15 1.57±0.68 148.38±43.91 4.33±0.91 3.04±0.25 16.05±2.10 173.37±44.28 3.59±0.67 177.38±18.45 4.46±0.76 3.66±0.3 22.94±2.89 212.03±18.66 6.67±3.78 235.69±74.46 7.36±2.54 3.12±0.42 12.16±3.59 265.01±76 2.09±0.28 168.33±9.34 3.69±0.40 2.62±0.10 6.85±0.33 183.59±9.66-5 44 029.73 58 504.19 40 870.06 2 000-903.36±105.94 3.73±0.92 0.06±0.02 0.08±0 9.44±2.04-314.86±136.69 3.37±1.00 0.07±0.02 0.07±0.01 8.03±1.05-717.89±134.5 4.03±0.42 0.08±0.01 0.09±0.01 11.47±1.44-1 334.36±755.67 5.35±1.69 0.13±0.04 0.08±0.01 6.08±1.8-418.8±56.88 3.82±0.21 0.06±0.01 0.06±0 3.42±0.17-2.41±0.59 166.68±16.92 13.68±1.48 5.13±4.10 6.51±0.75 1.29±0.49 7.30±0.10 16.63±1.31 219.63±17.96 2.84±0.58 232.41±38.96 19.82±3.75 4.62±2.88 5.42±1.29 0.88±0.67 7.29±0.07 16.10±1.41 289.38±40.25 2.98±0.50 125.82±16.27 9.83±1.46 3.56±2.88 6.49±1.63 0.74±0.49 7.44±0.09 18.51±0.98 175.35±20.35 5.15±1.91 306.96±31.74 27.59±3.05 3.29±1.33 5.90±1.94 0.18±0.18 5.85±2.95 15.36±4.91 370.28±37.20 2.46±0.23 198.45±11.74 16.54±0.93 3.65±0.26 10.95±1.12 0.79±0.40 7.25±0.12 13.74±0.52 253.84±13.80 4 203.1 262-1 433.65 406.83 908.83 125.21 28 922.73-0.57±0.14 636.17±64.57-3.58±2.86 16.01±1.85 1.41±0.53 58.33±0.82 0.58±0.05-0.68±0.14 887.06±148.71-3.22±2.01 13.32±3.17 0.97±0.74 58.22±0.58 0.56±0.05-0.71±0.12 480.24±62.11-2.48±2.01 15.94±4.01 0.81±0.54 59.42±0.68 0.64±0.03-1.23±0.46 1 171.61±121.15-2.29±0.93 14.5±4.78 0.20±0.20 46.70±23.6 0.53±0.17-0.59±0.06 757.44±44.82-2.55±0.18 26.92±2.75 0.87±0.44 57.92±0.94 0.47±0.02-15.62±3.71 3.03±0.06 18.66±3.70 11 157.48 18.06±4.87 3.02±0.14 21.08±4.82 9 125.45 6.55±2.16 3.05±0.12 9.60±2.16 11 527.32 29.79±5.08 3.13±0.30 32.92±5.10 11 716.49 51.27±5.34 3.30±0.27 54.57±5.38 11 010.40 13.41 122.74 1 165.08±276.96 24.72±0.50 1 346.45±363.39 24.64±1.15 488.77±161.17 24.85±0.96 2 221.58±379.05 25.48±2.43 3 823.31±398.19 26.85±2.20——————
XJ和MT产区酒样酯类和酸类化合物含量明显高于其他产区的原因主要是乙酸乙酯和乙酸含量高,乙酸乙酯和乙酸是酱香白酒中占比最大的酯类和酸类物质。此外,OAV结果表明,乙酸苯乙酯只对XJ产区酒样有香气贡献,乙酸苯乙酯在酿造过程中合成率较低,XJ产区酒样可能具有更多适宜具有合成乙酸苯乙酯能力的微生物。通过计算OAV可知,正戊醇和2-庚醇对JS和LJ产区酒样有香气贡献,对XJ和MT产区酒样则无香气贡献,这可能是造成它们香气差异的主要物质。辛醇对JS产区酒样无香气贡献,而对其余酒样均有贡献,苯甲醛只对JS产区酒样有香气贡献。
为进一步分析赤水河五个不同产区的酱香白酒的香气特征化合物,以对五种酒样的香气具有贡献的化合物(OAV>1)做载荷分析,结果见图1。由表1和图1可知,碳原子数>4的酯类物质对LJ产区酒样影响较大,包括辛酸乙酯、戊酸乙酯、己酸乙酯,其中己酸乙酯呈水果香和花香,LJ产区酒样果香味较浓,己酸乙酯含量明显高于其他几个酒样,分析原因可能是LJ产区酒样既生产酱香型白酒又生产浓香型白酒,己酸乙酯是浓香型白酒中酯类物质含量最高的物质,在LJ产区酿造环境和酿造工艺中存在一些己酸乙酯的代谢物。乙酸、乙酸乙酯与丙酸、丙酸乙酯之间会发生水解和酯化反应,该四种化合物与MT产区酒样相关性较大,异戊醛、异戊醇、丁酸乙酯、糠醛对JS产区酒样贡献较大,与ZJ和XJ产区酒样相关性较大的多为一些酯类化合物,如异丁酸乙酯、癸酸乙酯、壬酸乙酯、异戊酸乙酯、癸酸乙酯、乳酸乙酯。
图1 基于OAV>1的特征化合物不同产区酒样的主成分分析载荷图
Fig.1 Load diagram of principal component analysis of Baijiu samples from different regions based on characteristic compounds with OAV>1
2.2.1 主成分分析
对赤水河流域五个不同产区的酱香型白酒中的64个挥发性风味物质含量的数据进行处理,PCA散点图见图2。由图2可知,五个产区白酒样品大致被分成三类,LJ产区酒样位于坐标轴左下方,JS产区酒样集中位于坐标轴左上轴。MT产区酒样位于坐标轴右上半轴,XJ和ZJ产区酒样基本位于坐标轴下半轴。ZJ、XJ、MT三大产区酒样则集中分布在原点附近,聚集度较高。JS产区酒样组内分布较为散乱,可能是JS产区酒样中存在部分年份白酒,放置一段时间后风味物质含量发生了改变。据文献报道,金沙地区酱香白酒分布较为集中,茅台镇核心地区、习水地区分布聚拢,这可能是因为金沙酱酒与其他几个地区酱香白酒的矿质含量有差异[17]。由此可见,JS与LJ产区酒样组间差异大,与其余三个产区分布明显区分开,XJ、MT和ZJ三个产区酒样的风味物质分布较一致的原因可能是地理位置接近,环境条件以及微生物种类较相似。
图2 不同产区酒样挥发性风味物质主成分分析散点图
Fig.2 Scatter plots of principal component analysis of volatile flavor compounds in Baijiu samples from different producing regions
2.2.2 模型构建及验证
总体来看,虽然五个产区白酒大致被区分开,但组间和组内之间分布都较为散乱,为建立一个快速准确的赤水河流域产区鉴别模型,采用Hotelling's T2和DmodX法排除异常值,主成分分析的异常值筛选见图3。由图3可知,排除异常值15个,剩下116个训练集、15个测试集基于挥发性风味物质含量结果对不同赤水河流域产区酱香型白酒构建模型。
图3 主成分分析的异常值筛选
Fig.3 Abnormal value screening of principal component analysis
a为Hotelling's T2法;b为DmodX法。
筛选完异常值后,将剩余训练集样品的挥发性风味化合物定量结果构建OPLS-DA模型,结果见图4。由图4可知,MT产区酒样基本与XJ和ZJ产区酒样分离开,但XJ和ZJ产区酒样分离效果依然不佳。总体来说模型分类良好,且模型的解释变量(R2Y)和预测能力(Q2Y)分别为0.966和0.945。
图4 不同酒样挥发性风味化合物正交偏最小二乘-判别分析
Fig.4 Orthogonal partial least squares-discriminant analysis of volatile flavor compounds in different Baijiu samples
此外进行了200次的置换检验,结果见图5。由图5可知,置换检验的R2截距为0.926,Q2截距为-0.454,且图中左侧的点均低于右侧的点,说明模型对数据未过拟合。
图5 置换检验结果
Fig.5 Results of the permutation test
将从每个产区随机抽取三个观测值作为测试集代入上述模型中,以评估模型对于不同产区的识别情况,五个不同产区的训练集和测试集在两个模型中的分布情况见图6。由图6可知,大多数酒样正确分布到相应的区域,个别酒样除外。
图6 测试集酒样分布
Fig.6 Distribution of test set Baijiu samples
VIP值可以量化PLS-DA的每个变量对分类的贡献度[18],VIP值越大,认为变量在不同产地间差异越显著[19],以VIP值>1为筛选标准,确认了28种对赤水河流域不同产区白酒具有贡献的物质。赤水河流域不同产区酱香型白酒风味物质VIP分布图见图7。
图7 基于正交偏最小二乘-判别分析不同酒样挥发性风味物质变量投影重要性值
Fig.7 Variable importance in the projection values of volatile flavor compounds in different Baijiu samples based on orthogonal partial least squares-discriminant analysis
由图7可知,28种对赤水河流域不同产区白酒具有贡献的物质为2-丁基呋喃、壬酸乙酯、癸酸、2-壬酮、乳酸乙酯、仲丁醇、乙酸、乙偶姻、肉豆蔻酸乙酯、2,3,5,6-四甲基吡嗪、癸酸乙酯、异戊醛、乙酸乙酯、己醛、2,3,5-三甲基吡嗪、辛酸乙酯、正丙醇、邻苯二甲酸二丁酯、异戊醇、异丁醛、2-戊醇、苯乙酸乙酯、2-壬醇、庚酸乙酯、异戊酸乙酯、己酸、糠醛、丁酸。其中,2-丁基呋喃(1.60)、壬酸乙酯(1.40)、癸酸(1.30)对区分不同产区酒样贡献度更大。
为进一步探索赤水河流域不同产区酱酒风味差异物,根据VIP值筛选出来的28种潜在差异物绘制聚类热图,直观展示了64种挥发性风味成分在不同产区间的分布规律,结果见图8。由图8可知,不同产区的酒样之间分类效果较好,LJ和JS产区酒样明显区分于其他三个酒样,MT、ZJ和XJ产区酒样风格类似,其中MT和ZJ产区酒样聚类效果更佳。不同产区酱香白酒中的关键差异物是造成不同酒样品质风格不同的主要原因,五个产区被分成不同类别,进一步验证了赤水河流域不同酱香白酒产区风格差异显著性。
图8 不同酒样挥发性风味化合物的聚类分析热图
Fig.8 Heat map of cluster analysis of volatile flavor compounds in different Baijiu samples
从红到蓝颜色越深代表含量越高,颜色越浅代表含量越低。
其中庚酸乙酯、辛酸乙酯、肉豆蔻酸乙酯、己酸在LJ和JS产区酒样含量明显高于其他三个产区,根据定量结果可知,LJ和JS产区酒样的己酸乙酯含量偏高,己酸乙酯是浓香型白酒中酯类物质含量最高的物质,在LJ产区酿造环境和酿造工艺中存在一些己酸乙酯的代谢物,这与唐平等[20-21]研究结果一致。JS产区酒样的异丁醛、异戊醛、糠醛和异戊醇含量较高,糠醛是酱酒中的主要呋喃类物质,高沸点、不易挥发,具有空杯留香的特点,其在JS产区酒样中的含量远远高于其他产区,JS酱酒浓郁的空杯香很有可能就是高含量、高沸点的糠醛与其他香气物质形成的复合香气。这与周靖等[22]对金沙酱酒风味研究研究结果一致,XJ产区酒样中的乙酸和乙酸乙酯含量最高,乙酸是白酒中最丰富的酸,具有酸香和醋香气味,它由酵母菌在发酵过程中产生[23],脂肪酸通过酯类水解来促进XJ产区酒样的香气平衡,这对芳香酯含量有重要影响。MT产区酒样的吡嗪类物质含量相对于其他产区而言较高,大部分吡嗪类物质对人体健康有益,吡嗪类物质是酱香白酒焦香风味的主要来源[24]。异戊醛、异戊醇、丁酸乙酯、糠醛对JS产区酒样贡献较大,这与张情亚等[25]研究结果一致。
本研究围绕了赤水河流域中的JS、LJ、ZJ、MT、ZJ五个不同产区代表性的52个样品酱香型白酒中的64种挥发性物质定量结果进行分析。主成分分析结果显示,JS产区和LJ产区酒样明显与其他产区酒样分开,ZJ、XJ产区酒样聚集度较高,MT产区酒样靠近ZJ和XJ产区酒样,总体来说,五个产区的酱香白酒风味分布有一定的差异性。利用PCA和OPLS-DA建立赤水河流域不同酱香白酒产区分类模型,部分酒样作为验证模型,模型预测率为100%,基本可以实现不同产区的酱香白酒鉴别。根据VIP值筛选出的28种重要差异物做聚类分析热图以直观表达不同流域产区白酒风味成分分布规律,结合载荷图和OAV分析最终确立了19个对贵州不同产区风味物质差异具有贡献的物质,包括6种酯类、3种酸类、2种醛类、3种醇类、2种吡嗪类、1种呋喃类、2种芳香族类物质。本研究可对赤水河不同流域产区白酒风味差异提供一定的数据和理论支持,丰富了不同产区白酒的地理环境存在差异,为了建立更具实用价值、更稳定的白酒判别模型,还需扩大采样范围和样本量进行更深层次的研究。
[1]BOKULICH N A,THORNGATE J H,RICHARDSON P M,et al.Microbial biogeography of wine grapes is conditioned by cultivar,vintage,and climate[J].PNAS,2014,111(1):E139-E148.
[2]孙细珍,杜佳炜,张帆,等.应用顶空固相微萃取-气相色谱质谱技术解析天龙泉米香型白酒风味成分[J].酿酒,2021,48(5):55-60.
[3]胡瑾,马一飞,罗佳雪,等.气相色谱法同时测定白酒中57种风味物质[J].中国酿造,2022,41(5):206-211.
[4]洪泽淳,熊含鸿,李南,等.基于全二维气相色谱-飞行时间质谱对酱香型白酒风味成分的研究[J].酿酒科技,2018(4):121-125.
[5]赵宁,解菲,梁鸡保,等.基于主成分分析的21种粳性黄米品质评价[J].食品与机械,2022,38(8):22-27,33.
[6]ZHANG Z Y,SHA M,WANG Y H,et al.Chemical perturbation two-dimensional correlation ultraviolet visible spectroscopy for quality control of Chinese liquor[J].J Am Soc Brew Chem,2018,76(2):141-146.
[7]雷亚兰,周志梅,李瑾,等.基于主成分分析和聚类分析方法评价宝庆桂丁绿茶品质特性[J].食品工业科技,2022,43(6):269-277.
[8]雷春妮,王波,解迎双,等.化学计量学结合气相色谱法对12种植物油判别模型的构建[J].中国粮油学报,2022,37(6):170-176.
[9]李爱兰,乔文君,汪兴财,等.吹扫捕集-气质联用结合化学计量学用于不同产地和品牌浓香型白酒鉴别[J].化学试剂,2021,43(5):576-583.
[10]程平言,范文来,徐岩.基于质谱与化学计量学的白酒原产地鉴定[J].质谱学报,2014,35(1):32-37.
[11]张晓婕,邱树毅,曾庆军,等.不同工艺酱香型白酒挥发性物质差异分析[J].食品科学,2022,43(18):279-285.
[12]范文来,徐岩.白酒79个风味化合物嗅觉阈值测定[J].酿酒,2011,38(4):80-84.
[13] LIU H L, SUN B G.Effect of fermentation processing on the flavor of Baijiu[J].J Agr Food Chem,2018,66(22):5425-5432.
[14]李利利.基于风味组学对贵州产区酱香典型体基酒风味结构的解析[D].贵阳:贵州大学,2021.
[15]李莉,王秋叶,盛夏,等.白酒中酯类对酒质的影响[J].食品安全导刊,2016(36):124.
[16]朱金玉,解成玉,李玉英.白酒酒体的结构分析与设计原则[J].酿酒科技,2016(11):83-84.
[17]王丽,卢君,山其木格,等.赤水河流域不同地区酱香型白酒矿质元素特征分析及产地辨识初探[J].中国酿造,2021,40(9):58-62.
[18]王梦琪,邵晨阳,朱荫,等.龙井茶香气成分的产区差异分析[J].茶叶科学,2018,38(5):508-517.
[19]张晨霞,王国成,王超,等.普洱市不同产茶区普洱生茶香气成分差异性分析[J].食品研究与开发,2020,41(1):177-184.
[20]唐平,卢君,毕荣宇,等.赤水河流域不同地区酱香型白酒风味化合物分析[J].食品科学,2021,42(6):274-281.
[21] ZHU J C, NIU Y W, XIAO Z B.Characterization of important sulfur and nitrogen compounds in Lang Baijiu by application of gas chromatography-olfactometry,flame photometric detection, nitrogen phosphorus detector and odor activity value[J].Food Res Int,2020,131:109001.
[22]周靖,刘义,李世杰,等.基于GC-MS对金沙酱酒醇柔风格的分析研究[J].中国酿造,2022,41(8):229-234.
[23]ONYENEKWE P C,ERHABOR G O,AKANDE S A.Characterisation of aroma volatiles of indigenous alcoholic beverages:burukutu and pito[J].Nat Prod Res,2016,30(5):558-564.
[24]许佩勤,方毅斐,庄俊钰,等.白酒中吡嗪、呋喃类化合物的分析方法研究[J].食品与发酵科技,2017,53(2):85-93.
[25]张情亚,余德顺,邱树毅,等.利用风味物质进行金沙窖酒原产地溯源[J].食品与发酵工业,2016,42(12):157-162.
Analysis on the difference of volatile flavor compounds in sauce-flavor Baijiu from different production areas in Chishui-river basin