羅彥波 陳小靜 李翔宇 張洪非 朱風鵬 姜興益 龐永強

摘 要 建立了化學衍生結合液相色譜-串聯(lián)質譜（LC-MS/MS）檢測加熱卷煙煙氣中9種醛酮類香味成分的分析方法。樣品和標準溶液中的目標分析物分別用1-（2-肼基-2-氧代乙基）吡啶-1-溴化物（吉拉德P試劑，GP）和d5-GP衍生，以標準溶液的衍生化產物作為穩(wěn)定同位素內標，采用多反應監(jiān)測模式檢測。在最優(yōu)的衍生化反應條件下，9種醛酮類香味成分在各自濃度范圍內線性關系良好，標準曲線相關系數(shù)（R2）≥0.9937，檢出限為0.07～12.66 ng/mL，相對標準偏差為0.2%～11.6%。9種醛酮類香味成分在3個加標水平下的回收率為82.5%～118.0%。本方法可用于加熱卷煙煙氣中醛酮類香味成分的檢測。

關鍵詞 化學衍生;液相色譜-串聯(lián)質譜;加熱卷煙;醛酮類香味成分

1 引 言

煙草及煙氣品質的控制和評價是煙草行業(yè)關注的重要問題[1]，其中，非煙氣的化學組成與煙氣品質密切相關，煙氣的化學組成影響煙氣的香味和吃味，因此，煙氣香味成分的研究一直備受重視。醛酮類香味成分是煙草中數(shù)量最多的一類香味成分[2]，目前，煙氣中已鑒定出的醛酮類香味成分超過500種。許多重要香味成分屬于醛酮類[3]，如3-甲基-2-環(huán)戊烯-1-酮在煙氣中能產生甜的、焦糖-楓槭的氣息;4-庚酮則能產生甜的、水果的氣息;一些高沸點的酮（如紫羅蘭酮）也有明顯的致香作用，不少已用于煙草制品的加香中[2]。有研究表明，醛酮類香味成分的組成和相對含量綜合影響煙氣的品質香味情況[4]。因此，有必要對煙氣中醛酮類香味成分進行分析研究。

目前，醛酮類化合物的檢測方法主要有氣相色譜-質譜聯(lián)用法[5～9]、高效液相色譜-二極管陣列檢測/紫外檢測法[10～13]和高效液相色譜-質譜聯(lián)用法[14～16]等。由于醛酮類化合物缺少發(fā)色團和易電離基團，因此，在進行紫外和質譜分析前，需要通過衍生引入發(fā)色團或易電離基團，以提高紫外檢測的吸光度和質譜檢測的離子化效率[17，18]。質譜檢測常用的衍生化試劑有亞硫酸銨[19]、伯胺類[20，21]、肼類[22～29]和酰肼類[30～34]等。吉拉德P試劑（GP）是常用的酰肼類衍生化試劑，其結構中含有恒帶正電荷的季銨鹽結構，因而能夠顯著提高衍生化產物在質譜中的離子化效率和檢測靈敏度。GP已被用于人和牛體液中類固醇激素[31，35]、血漿和干血點中性激素[36，37]等的分析。結果表明，GP可選擇性地與醛酮類化合物高效反應，在分析上述復雜樣品中具有優(yōu)良的表現(xiàn)。然而，目前尚未見GP用于卷煙煙氣中醛酮類香味成分檢測的報道。

煙氣中的醛酮類化合物主要是碳水化合物在加熱條件下產生的，當加熱溫度高于300℃時，醛酮類化合物是構成煙氣香味的主要成分之一[38]。與傳統(tǒng)卷煙相比，加熱卷煙煙氣中大多數(shù)化學成分的釋放量顯著降低，而以“加熱而非燃燒煙絲”為特征的加熱卷煙（加熱溫度約為350℃）在加熱過程中具備醛酮類香味成分生成的物質基礎和形成條件，因此，加熱卷煙煙氣中醛酮類香味成分的降低比例有限[39，40]。由于制造工藝和添加原料的差異，不同加熱卷煙煙氣中的醛酮類香味成分的種類也不同。目前，煙氣中醛酮類物質的分析方法主要采用外標法或單一內標法進行定量分析，而且目標分析物多為有害的羰基化合物，個別針對香味成分的檢測方法存在目標分析物的種類和結構較為單一的問題。

因此，本研究利用GP和d5-GP分別衍生樣品和標準溶液中的醛酮類香味成分，并以d5-GP的衍生化產物為目標分析物的穩(wěn)定同位素內標。在目標分析物的選擇方面，本研究以煙草制品許可使用的醛酮類化合物和煙草在加熱過程中生成的醛酮類化合物為基礎，同時考慮對卷煙品質的影響和化學結構（不同碳鏈長度、含有其它不同的官能團、不同羰基或醛基數(shù)目），選擇了如圖1所示的9種代表性目標分析物。其中，3-甲基-2-環(huán)戊烯-1-酮、庚醛和4-庚酮為煙草在加熱過程中生成的目標分析物，其余6種化合物為煙草制品許可使用的醛酮類化合物添加劑。將GP應用于加熱卷煙中醛酮類香味成分的檢測，通過衍生化反應向目標分析物分子中引入同位素基團，并作為內標，建立了加熱卷煙煙氣中醛酮類香味成分的分析方法。本方法顯著提高了目標分析物的靈敏度，對煙氣中醛酮類香味成分的分離和分析具有參考意義。

2 實驗部分

2.1 儀器與試劑

Agilent 1200液相色譜儀（美國Agilent公司）和AB SCIEX 4000三重四極桿質譜儀（美國AB公司），配備有G1311A輸液泵、G1322A在線脫氣機、G1329A自動進樣器、G1316A柱溫箱和Analyst 1.5.2工作站。CP224S電子天平（德國Sartorius公司）;SM450型20通道吸煙機（配加熱卷煙分析套件，英國Cerulean公司）; 44 mm劍橋濾片（德國Borgwaldt-KC公司）;SB25-12DTD超聲波清洗機（寧波新芝生物科技股份有限公司）;SHA-B恒溫振蕩器（常州國華電器有限公司）;TurboVap LV氮氣濃縮儀（瑞典Biotage公司）;Heidolph多位渦旋儀（德國Heidolph公司）。

GP和d5-GP（純度≥99%，美國Avanti Polar Lipids Inc公司）;9種醛酮類化合物：3-羥基-2-丁酮（純度≥99.4%）、3-甲基-2-環(huán)戊烯-1-酮（純度≥99.4%）、2，3-丁二酮（純度≥97.0%）、苯甲醛（純度≥995%）、己醛（純度≥95.0%）、4-庚酮（純度≥99. 8%）、庚醛（純度≥97.0%）、壬醛（純度≥99.5%）和紫羅蘭酮（純度≥96.0%），均購自美國Sigma-aldrich公司;甲醇和乙腈（色譜純，德國Merck公司）;甲酸（色譜純，美國Aladdin公司）;冰醋酸（色譜純，美國Fisher Scientific公司）;實驗用水為Q-POD Milli-Q超純水儀（美國Millipore公司）制得的超純水。

2.2 儀器條件

液相色譜條件：色譜柱為Shimadzu VP-ODS（150 mm × 2.0 mm，5 μm），柱溫為40℃，流速為0.2 mL/min，流動相A和B分別為0.1%（V/V）甲酸和乙腈。梯度洗脫：0～5 min，5% B;5～22 min，5%～40% B;22～40 min，40%～90% B; 40～45 min，90% B; 45～48 min，90%～5% B; 48～60 min，5% B。進樣體積為10 μL。

質譜條件：離子源為電噴霧離子源;電噴霧電壓為5500 V，氣簾氣壓力為68.95 kPa，霧化氣壓力為75.84 kPa，碰撞氣（Ar氣）壓力為41.37 kPa;掃描模式為多反應監(jiān)測;其它質譜參數(shù)見表1。

2.3 樣品制備

本研究的目的是建立加熱卷煙中醛酮類香味成分的分析方法，因此，直接采用文獻報道的方法對加熱卷煙進行抽吸，并對其煙氣冷凝物進行提取[41，42]。由于加熱卷煙煙氣冷凝物的基質比傳統(tǒng)卷煙煙氣簡單，主要為水分、丙二醇和丙三醇等。因此，冷凝物提取液直接進行衍生化處理。具體步驟：采用ISO抽吸模式規(guī)定的抽吸參數(shù)對加熱卷煙樣品進行抽吸，再將捕集有煙氣總粒相物的濾片用15 mL甲醇超聲提取30 min，靜置，取1.0 mL上清液進行衍生。衍生時，在上清液中依次加入125 μL冰醋酸和GP溶液（3.0 μmol/L，溶于甲醇），渦旋混合后在40℃下振蕩反應2 h，氮氣吹干，以1.0 mL甲醇-水（1∶1，V/V）溶液復溶。在對標準溶液進行衍生時，稱取適量目標分析物標準品，用甲醇配制成不同質量濃度

的系列混合標準工作溶液，再分別取系列混合標準工作溶液進行衍生，衍生步驟同上。用d5-GP溶液（3.0 μmol/L，溶于甲醇）代替GP溶液對混合標準溶液進行衍生，衍生步驟同上。分別取20 μL復溶后的d5-GP衍生化產物添加至GP衍生后的系列混合標準工作溶液和樣品溶液，混勻，待測。

3 結果與討論

3.1 GP標記的醛酮化合物的質譜特征

GP中帶孤電子對的氮原子可以進攻羰基基團上帶正電荷的碳原子發(fā)生親核加成反應，生成的羥胺類化合物經(jīng)脫水后形成穩(wěn)定的腙，其反應示意圖如圖2所示。經(jīng)衍生化反應后目標分析物在結構中引入了帶正電荷的季銨基團，可以提高質譜響應信號。因此，本研究首先考察了GP標記的目標分析物的二級碎裂行為，結果如圖3所示。所有的目標分析物被GP標記后都能產生一個丟失79 Da的子離子和相同的碎片離子m/z 80，這是由衍生化試劑中的吡啶環(huán)引起的[30，35]，說明衍生化試劑GP能夠標記本實驗中的目標分析物。標記后的目標分析物具有較好的碎裂行為，丟失吡啶環(huán)后，離子[M·GP-Br-Py]的質譜響應信號相對較強，而且附近沒有其它碎片離子的干擾，這有利于對標記后的目標分析物進行多反應監(jiān)測模式檢測。

3.2 衍生化條件的優(yōu)化

為了得到較好的標記效率，考察了醋酸和GP的用量、反應溫度和時間等影響衍生化效率的因素。

在衍生化反應中，醋酸可以提供酸性的催化環(huán)境，在冰醋酸占整個反應體系體積的0～10%（V/V）范圍內，考察了冰醋酸用量的影響。結果表明，當冰醋酸的體積為5%（V/V）時，衍生化產物與內標的峰面積比值達到最大，繼續(xù)增大醋酸用量，峰面積的比值基本保持不變。

通過改變衍生化試劑的用量，使GP與目標分析物的物質的量的比值為10～1000。結果表明，隨著衍生化試劑用量的增加，衍生化產物與內標的峰面積比值逐漸增加，當物質的量的比值達到500時，峰面積的比值基本不再變化。標準工作曲線繪制以及實際樣品分析時衍生反應中，為保證衍生化反應完全，選擇GP與待測物的物質的量的比值遠超過500。

文獻[30，35～37]報道的GP試劑標記的反應溫度為室溫～70℃，由于本研究的衍生化反應在甲醇（沸點為65℃）中進行，因此，選擇的反應溫度范圍為30℃～60℃。結果表明，在30℃～60℃范圍內，反應溫度對衍生化產物與內標的峰面積比值影響不大，為了保證穩(wěn)定的衍生化效率，選擇40℃作為后續(xù)的反應溫度。

考察了反應時間對衍生化效率的影響。結果表明，在0.5～6 h范圍內，當反應時間從0.5 h延長到3 h時，衍生化產物與內標的峰面積比值逐漸增加，繼續(xù)延長反應時間，峰面積比值增加不明顯。

綜上，選擇的最優(yōu)的衍生化條件為：冰醋酸體積為5%（V/V），反應溫度為40℃，反應時間為3 h。

3.3 衍生前后靈敏度比較

由于衍生化產物中含有恒帶正電荷的季銨鹽基團，因此，衍生化處理可增加醛酮化合物在電噴霧離子源中的離子化效率和質譜響應，有助于解決醛酮類香味成分因缺少易電離基團而無法用質譜直接檢測的難題。將9種醛酮類香味成分用質譜進行檢測，目標分析物標準溶液未經(jīng)衍生直接進行液相色譜-質譜分析，只檢測到3-甲基-2-環(huán)戊烯-1-酮和紫羅蘭酮的質譜信號;目標分析物標準溶液經(jīng)過GP衍生后進行液相色譜-串聯(lián)質譜分析，所有目標分析物均有質譜響應。經(jīng)過衍生后，3-甲基-2-環(huán)戊烯-1-酮和紫羅蘭酮的靈敏度分別提高了36和27倍，其余7種化合物（直接進行液相色譜-質譜分析沒有信號）在衍生化之后的質譜信號都比3-甲基-2-環(huán)戊烯-1-酮強，可以推測其靈敏度提高的倍數(shù)更大。以上結果說明，GP衍生化處理可顯著提高目標分析物的靈敏度。

3.4 方法分析性能

考察了GP標記后衍生化產物的穩(wěn)定性。將標記后的產物在0～4℃下避光保存，分別在第0、6、12和24 h進行液相色譜-串聯(lián)質譜分析。結果表明，9種衍生化產物的峰面積沒有明顯變化，相對標準偏差在1.0%～11.9%之間，說明此保存條件下衍生化產物在24 h內穩(wěn)定，保證了方法的準確性和重復性。由于經(jīng)過衍生化處理的樣品可在24 h內進行液相色譜-串聯(lián)質譜分析，因此，本研究未對樣品進行更長時間的穩(wěn)定性考察。

在優(yōu)化條件下，考察了方法的檢出限、定量限和線性范圍。檢出限和定量限為衍生化產物信噪比為3和10時對應的濃度。以目標分析物濃度為橫坐標，以目標分析物與內標峰面積之比為縱坐標，繪制工作曲線。如表2所示，9種目標物的檢出限在0.07～12.66 ng/mL之間，定量限在0.22～42.19 ng/mL之間，線性相關系數(shù)R2≥0.9937。

為了考察本方法的重現(xiàn)性，制備3種濃度（各目標分析物線性范圍最低值的2、10和40倍）的樣品，以1天內制備的5個樣品和連續(xù)3天配制的樣品進行分析，計算日內和日間精密度（相對標準偏差）。如表3所示，目標分析物的日內及日間精密度分別不大于10.8%和11.6%，說明本方法重現(xiàn)性良好。

3.5 實際樣品分析

采用本方法對兩種典型加熱卷煙（iQOS和glo）進行分析。在iQOS和glo煙氣中檢測到3-羥基-2-丁酮，其含量均低于本方法的定量限;在glo煙氣中檢測到2，3-丁二酮，其含量為39.18 ng/支。

如圖4所示，目標分析物的定量分析不存在干擾。由于目前沒有分析加熱卷煙煙氣中醛酮類香味成分的標準方法，文獻[39]報道glo煙氣中2，3-丁二酮的含量為33.60～43.80 ng/支。以iQOS和glo進行回收率測定，在超聲提取階段加入9種目標分析物的混合標準溶液，3種加標水平分別為各目標分析物線性范圍最低值的2、10和40倍，每個濃度水平重復測定5次。如表4所示，實際樣品中目標分析物在不同濃度下的回收率為82.5%～118.0%，說明本方法的準確度好，可用于分析加熱卷煙煙氣中的醛酮類香味成分。

4 結 論

采用化學衍生標記技術，建立了液相色譜-串聯(lián)質譜法同時測定醛酮類香味成分含量，并用于實際加熱卷煙煙氣樣品的分析。本方法具有靈敏度高、快速、準確的優(yōu)點，可以滿足加熱卷煙煙氣中醛酮類香味成分高靈敏度檢測的需求，同時也為其它類似樣品中醛酮類化合物的分析提供了參考。

References

1 LI Ang-Ze，WANG Xian-Shuang，XU Xiang-Jun，HE Ya-Ge，GUO Shuai，LIU Yu-Fei，GUO Wei，LIU Rui-Bin. Chin. Opt.，2019，12（5）： 1139-1146

李昂澤，王憲雙，徐向君，何雅格，郭 帥，柳宇飛，郭 偉，劉瑞斌. 中國光學，2019，12（5）： 1139-1146

2 SHI Hong-Zhi，LIU Guo-Shun，YANG Hui-Juan，JI Xiao-Ming. Tobacco Flavor. Beijing： China Agriculture Press，2011

史宏志，劉國順，楊惠娟，姬小明. 煙草香味學. 北京：中國農業(yè)出版社，2011

3 Perfetti T A，Rodgman A. Contrib. Tob. Res.， 2011，24 （5）： 215-232

4 YE Rong-Fei，ZHAO Rui-Feng. Guangdong Agr. Sci.， 2011，5： 51-53

葉榮飛，趙瑞峰. 廣東農業(yè)科學，2011，5： 51-53

5 Filipowska W，Jaskula-Goiris B，Ditrych M，Schlich J，De Rouck G，Aerts G，de Cooman L. J. Chromatogr. A，2020，1612： 460647

6 Lomonaco T，Romani A，Ghimenti S，Biagini D，Bellagambi F G，Onor M，Salvo P，F(xiàn)uoco R，Di Francesco F. J. Breath Res.，2018，12（4）： 046004

7 CUI Ting-Hui，LI Heng，TANG Shi-Yun，LIU Chun-Bo，GAO Xiao-Hua，LI Zhen-Jie，ZHANG Feng-Mei，YANG Guang-Yu，JIANG Wei，LIU Zhi-Hua，MIAO Ming-Ming，HE Pei，SHEN Qin-Peng. Journal of Instrumental Analysis，2018，37（7）： 766-771

崔婷惠，李 恒，唐石云，劉春波，高曉華，李振杰，張鳳梅，楊光宇，蔣 薇，劉志華，繆明明，何 沛，申欽鵬. 分析測試學報，2018，37（7）： 766-771

8 Moreira N，Araujo A M，Rogerson F，Vasconcelos I，F(xiàn)reitas V D，Pinho P G. Food Chem.，2019，270： 518-526

9 Pang X，Lewis A C，Shaw M D. J. Sep. Sci.，2017，40（3）： 753-766

10 Liu G，Luan B，Xin L，Li M，Tian Q，Liu Y，Xu Y. J. Chromatogr. A，2019，1594： 112-119

11 Prieto-Blanco M C，Penafiel Barba S，Moliner-Martinez Y，Campins-Falco P. J. Chromatogr. A，2019，1596： 241-249

12 FANG Si-Yang，CHEN Shi-Li，HE Jia-Hua，WANG Hao，SHEN Li-Si，LIU Yang-Yang，SUN Hui，LAI Jia-Ping. Chinese Journal of Chromatography，2019，37（6）： 634-643

方思煬，陳詩麗，何嘉鏵，王 浩，沈麗思，劉洋洋，孫 慧，賴家平. 色譜，2019，37（6）： 634-643

13 Liu J F，Yuan B F，F(xiàn)eng Y Q. Talanta，2015，136： 54-59

14 Antonowicz S，Hanna G B，Takats Z，Bodai Z. Talanta，2018，190： 436-442

15 Zhao S，Dawe M，Guo K，Li L. Anal. Chem.，2017，89（12）： 6758-6765

16 CAO Xiao-Ji，ZHOU Ting，XU Jia-Hui，CHAI Yun-Feng，YE Xue-Min. Chinese J. Anal. Chem.，2019，47（3）： 394-402

曹小吉，周 婷，徐家輝，柴云峰，葉學敏. 分析化學，2019，47（3）： 394-402

17 Laghrib F，Lahrich S，El Mhammedi M A. J. Electrochem. Soc.，2019，166（15）： 1543-1551

18 El-Maghrabey M H，Kishikawa N，Kuroda N. Biomed. Chromatogr.，2020，34（3）： e4756

19 Hu Y N，Chen D，Zhang T Y，Ding J，F(xiàn)eng Y Q. Talanta，2020，206： 120172

20 Eggink M，Wijtmans M，Ekkebus R，Lingeman H，de Esch I J P，Kool J，Niessen W M A，Irth H. Anal. Chem.，2008，80（23）： 9042-9051

21 Eggink M，Wijtmans M，Kretschmer A，Kool J，Lingeman H，de Esch I J P，Niessen W M A，Irth H. Anal. Bioanal. Chem.，2010，397（2）： 665-675

22 van Leeuwen S M，Hendriksen L，Karst U. J. Chromatogr. A，2004，1058（1）： 107-112

23 Banos C E，Silva M. J. Chromatogr. A，2009，1216（38）： 6554-6559

24 Shibayama Y，Higashi T，Shimada K，Kashimada K I，Onishi T，Ono M，Miyai K，Mizutani S. J. Chromatogr. B，2008，867（1）： 49-56

25 Xu H，Song D，Cui Y，Hu S，Yu Q W，F(xiàn)eng Y Q. Chromatographia，2009，70（5）： 775-781

26 Xia L，Du Y Q，Xiao X H，Li G K. Talanta，2019，202： 580-590

27 CHEN Hui，DENG Qi-Xin，LIU Jiang-Sheng，WU Qing-Hui，KE Wen-Lin，HUANG Sheng-Han.Food Ind.，2019，40（1）： 311-314

陳 輝，鄧其馨，劉江生，吳清輝，柯文林，黃勝翰. 食品工業(yè)，2019，40（1）： 311-314

28 LI Ying，ZHAO Hao-Jun，LIU Fei，ZHAO Fang，YANG Wei-Hua，YANG Chao-Yi，YANG Ting. J. Food Saf. Qua.，2019，10（16）： 5523-5527

李 穎，趙浩軍，劉 飛，趙 芳，楊衛(wèi)花，楊朝義，楊 婷. 食品安全質量檢測學報，2019，10（16）： 5523-5527

29 ZHANG Xin，LI Hong，ZHANG Cheng-Long，ZHANG Yu-Jie，HE Zhen，GAO Rui，WANG Wen-Xing. Res. Environ. Sci.，2019，32（5）： 821-829

張 鑫，李 紅，張成龍，張玉潔，何 振，高 銳，王文興. 環(huán)境科學研究，2019，32（5）： 821-829

30 Guo N，Liu P，Ding J，Zheng S J，Yuan B F，F(xiàn)eng Y Q. Anal. Chim. Acta，2016，905： 106-114

31 Guo N，Peng C Y，Zhu Q F，Yuan B F，F(xiàn)eng Y Q. Anal. Chim. Acta，2017，967： 42-51

32 Cobice D F，Mackay C L，Goodwin R J A，McBride A，Langridge-Smith P R，Webster S P，Walker B R，Andrew R. Anal. Chem.，2013，85（23）： 11576-11584

33 Johnson D W. Rapid Commun. Mass Spectrom.，2007，21（28）： 2926-2932

34 Chen X，Shi Q，Zhou X B，Liu X Z. Roy. Soc. Open Sci.，2019，6（2）： 181832

35 Kirk J M，Tarbin J，Keely B J. Rapid Commun. Mass Spectrom.，2006，20（8）： 1247-1252

36 De Barber A E，Sandlers Y，Pappu A S，Merkens L S，Duell P B，Lear S R，Erickson S K，Steiner R D. J. Chromatogr. B，2011，879（17-18）： 1384-1392

37 Lai C C，Tsai C H，Tsai F J，Lee C C，Lin W D. Rapid Commun. Mass Spectrom.，2001，15（22）： 2145-2151

38 ZHOU Hui-Ming，HUA Qing，TAO Li-Qi，ZHENG Sai-Jing. Tob. Sci. Tech.，2019，52（5）： 67-76

周慧明，華 青，陶立奇，鄭賽晶. 煙草科技，2019，52（5）： 67-76

39 Savareear B，Lizak R，Brokl M，Wright C，Liu C，F(xiàn)ocant J F. J. Chromatogr. A，2017，1520： 135-142

40 Forster M，F(xiàn)iebelkorn S，Yurteri C，Mariner D，Liu C，Wright C，McAdam K，Murphy J，Proctor C. Regul. Toxicol. Pharm.，2018，93： 14-33

41 ZHANG Hong-Fei，JIANG Xing-Yi，PANG Yong-Qiang，LI Xiang-Yu，LUO Yan-Bo，ZHU Feng-Peng，YAN Rui-Bo，MA Tao，CHEN Xiao-Jing. Tob. Sci. Tech.，2018，51（9）： 40-48

張洪非，姜興益，龐永強，李翔宇，羅彥波，朱風鵬，閆瑞波，馬 濤，陳小靜. 煙草科技，2018，51（9）： 40-48

42 ZHOU Quan，YIN Dong-Hong，WEN Jian-Hui，LUO Wei，WANG Liang，LIU You-Nian，ZHONG Ke-Jun. Tob. Sci. Tech.，2018，51（8）： 61-69

周 全，銀董紅，文建輝，羅 瑋，王 亮，劉又年，鐘科軍. 煙草科技，2018，51（8）： 61-69