烤煙煙葉物理特性與氣象因子關系研究

郭東鋒，鄒 鵬，劉 非，吳克松，姚忠達，舒俊生，竇玉青，張忠鋒，閆 寧*

(1.安徽中煙工業有限責任公司 技術中心，安徽 合肥 230088；2.中國農業科學院 煙草研究所，山東 青島 266101)

烤煙煙葉物理特性與氣象因子關系研究

郭東鋒1，鄒 鵬1，劉 非1，吳克松1，姚忠達1，舒俊生1，竇玉青2，張忠鋒2，閆 寧2*

(1.安徽中煙工業有限責任公司 技術中心，安徽 合肥 230088；2.中國農業科學院 煙草研究所，山東 青島 266101)

為了探索煙葉物理特性與烤煙生產大田期氣象因素的關系，以我國13個典型植煙區域烤煙的3個部位樣品為研究對象，運用廣義可加模型對煙葉物理特性和氣象因子之間的關系進行了研究分析，結果表明：日照百分率、日降水量≥0.1 mm日數、平均相對濕度、極端最高氣溫、極端最低氣溫對煙葉葉面密度有重要影響。就部位而言，上部葉葉面密度對大田期的降水量比中下部煙葉更為敏感。平均最高氣溫、平均最低氣溫、極端最高氣溫、平均氣壓、平均水汽壓和平均風速對煙葉填充值有重要影響，中部葉的填充值與日降水量≥0.1 mm日數關系顯著；極端最高氣溫、平均最高氣溫和平均溫度對含梗率有重要影響，同時日降水量≥0.1 mm日數對中上部煙葉的含梗率也有顯著的影響。就部位而言，上部葉含梗率受氣象因素影響相對較弱。

物理特性；氣象因素；烤煙

煙草是我國重要的經濟作物，煙葉是卷煙工業企業重要的戰略資源。煙葉的物理特性是煙葉質量的一個重要方面，煙葉物理指標的優劣直接關系到煙葉加工性能，影響卷煙工業的成本控制、倉儲養護、打葉復烤、制絲切絲等重要環節。同時，煙葉的物理特性也是表征煙葉內在質量水平的重要部分。前人關于煙葉物理特性與內在化學成分[1-4]、感官質量[5]、外觀質量[1,6]等方面的關系已有不少研究，也有研究不同海拔高度對煙葉物理特性的影響[7-8]、不同打葉復烤工藝對煙葉物理特性的影響[9-10]、栽培技術措施對煙葉物理特性的影響[11-13]等，也有單項物理指標比如填充值變異研究[14]、物理特性指標綜合評價等[15-16]方面的研究。雖然圍繞植煙區域煙草生產適宜性評價開展了很多研究[17-19]，但是圍繞煙葉物理特性與植煙區域氣象因子之間關系的研究尚不多見。由于我國植煙區域廣泛，生態條件差異較大，氣象因素復雜多變，近年來隨著部分區域干旱、洪澇等自然災害頻發，對煙葉質量產生較大影響。為了探索氣象因素對煙葉物理特性的影響，本文以我國烤煙典型植煙區域標準等級樣品為代表，運用廣義可加模型研究了氣象因子與不同部位煙葉物理特性之間的關系，以期為煙草農業生產、工業企業原料采購調撥、打葉復烤等提供技術參考。

1 材料與方法

1.1 材料

材料來自2012年度宣城、郴州、三門峽、曲靖、文山、涼山、麗江、昆明、玉溪、三明、龍巖、遵義、許昌等13個產煙區域初烤中部煙葉B2F、C3F、X2F標準等級樣品。填充值測定樣品為相應樣品切絲后的煙絲樣品，其它物理特性測定樣品均為相應標準等級樣品。

1.2 方法

1.2.1 物理特性測定

1.2.1.1 葉面密度的測定 在來樣群體中篩選完整葉片，葉片不應有殘傷、破損、砂土率超標等現象。從待測樣品中選出20片完整葉片，放入恒溫恒濕箱48 h，恒溫恒濕箱條件設定為：(22.2±0.5)℃，濕度：(61.0±2.0)%。打孔器一個，測定打孔器面積為S(單位換算為m2)。測定步驟：每個樣品選取具有代表性的煙葉10片，在主脈兩側，盡量避免打孔到主脈或支脈上，自葉片基部至葉尖均勻地在主脈兩側打取10個小葉片；稱取空鋁盒重量，記為M1，將用打孔器取下的100片葉片放入鋁盒，稱重，記為M2；按如下公式計算煙葉葉面密度D(單位：g/m2)：

式中：D：煙葉葉面密度，單位為g/m2；M1：鋁盒質量，單位為g；M2：100片小葉片和鋁盒總質量，單位為g。

1.2.1.2 填充值的測定 參照行業標準《卷煙 煙絲填充值的測定》 YC/T 152─2001進行檢測。

1.2.1.3 含梗率的測定 樣品前處理同葉面密度測度。測定步驟：待樣品在恒溫恒濕箱48 h后取出，用電子天平稱量20片煙葉的總質量，記為M1。用手撕掉葉片主脈兩側的葉片，留下主脈，主脈兩側不應有葉片殘留。用螺旋測微器確定≥1.5 mm煙梗為測量目標，將符合要求的所有20片葉的煙梗稱重，記為M2。

按如下公式計算煙葉含梗率G，以百分數計：

式中：G：煙葉含梗率(%)；M1：煙葉總質量，單位為g；M2：煙梗總質量，單位為g。

1.2.1.4 平衡含水率的測定 將各樣品分別切絲后，在22 ℃、相對濕度(60±2)%條件下平衡48 h，按文獻[20]規定的方法測定其含水率。

1.2.2 烤煙大田期確定 根據各植煙區域生產實際情況確定大田期(表1)，并據此整理氣象數據。

1.2.3 氣象數據來源 17個氣象點數據取自中國氣象數據共享服務網(http://cdc.cma.gov.cn/)。

1.3 數據處理方法

利用Excel 2007軟件進行數據的描述統計分析、相關和回歸分析；作圖在R3.3軟件平臺進行。文中統計描述均在數據清洗、整理、填充缺失值后進行；方差分析采用Turkey法進行；相關分析均采用Pearson相關；廣義可加模型分析(GAM)采用R軟件包Mgcv，進行光滑樣條估計；其它統計分析采用R軟件基礎包進行分析。

2 結果與分析

2.1 大田期不同氣象因子和烤煙物理特性統計描述

13個植煙區域大田期確定見表1，對各植煙區域大田期氣象因子和煙葉物理特性指標的統計描述見表2，由表2可知，各氣象因子中數據變異達到50%以上的極端最低氣溫、降雨量存在廣泛的變異。氣象因子變異系數超過20%的有平均風速、日降水量≥0.1 mm日數、最大風速和最大日降水量，存在中等水平的變異，說明各個植煙區域在氣象因子上存在明顯差異。烤煙葉面密度變異超過20%，主要是上、中、下3個部位煙葉差異引起變異較大。平衡含水率、填充值變異相對較小，僅有5%左右的變異；含梗率變異系數為12.55%，說明煙葉物理特性指標葉面密度和含梗率差異較大。

表1 各植煙區大田期時間跨度確定

表2 大田期不同氣象因子和烤煙物理特性統計描述

2.2 烤煙煙葉不同部位物理特性方差分析結果

由于烤煙不同部位煙葉的物理特性存在較為明顯的差異，為了進一步印證此差異，基于樣品檢測的不同部位物理特性數據，對不同部位數據進行了方差分析，結果見表3。結果可知，葉面密度和填充值在部位間存在顯著性差異，葉面密度自上而下依次降低，填充值反之；含梗率下部葉顯著高于上部葉，下部葉與中部葉差異不顯著；平衡含水率中部葉顯著高于上部葉，中下部差異不顯著、上下兩部位差異不顯著。

表3 烤煙煙葉不同部位物理特性方差分析結果

注：表中不同大寫字母表示在0.01水平下差異顯著，不同小寫字母表示在0.05水平下差異顯著。

2.3 氣象因子與煙葉物理特性的廣義可加模型分析結果

廣義可加模型(General Addative Model，GAM)可對部分或全部的自變量采用平滑函數的方法建立模型，函數可以是非參數的形式，適用于多種分布類型、多種復雜非線性關系的分析，由于此種特性，GAM可以很好地刻畫自變量和依變量之間的線性或非線性的關系。

基于煙葉物理特性部位間的顯著性差異，運用廣義可加模型分析不同部位物理特性與氣象因子的關系，結果見表4～表7。由表4氣象因子與葉面密度GAM分析結果可知，上部葉葉面密度與日降水量≥0.1 mm日數關系達到P<0.001顯著水平，其偏差解釋率達到了96%；與平均相對濕度關系達到P<0.01顯著水平，其偏差解釋率達到了58.5%；與降水量關系達到了P<0.1顯著水平，其偏差解釋率達到了74.7%；與日照因子的關系中，與日照百分率關系達到P<0.05顯著水平；與日照時數關系不顯著，但是兩者偏差解釋率均超過了40%，分別達到了41.8%和48.3%。中部葉葉面密度與氣象因子中的極端最低溫關系達到了P<0.05顯著水平，其偏差解釋率達到了93.2%，與日照百分率關系達到了P<0.05顯著水平，其偏差解釋率達到了28.3%；下部葉葉面密度與氣象因子中的極端最高溫關系達到了P<0.1顯著水平，其偏差解釋率達到了86.8%，另外達到P<0.1顯著水平的氣象因子還有日降水量≥0.1 mm日數、平均相對濕度和日照百分率，三者的偏差解釋率均超過20%。綜合來看，對3個部位有顯著影響的氣象因子是日照百分率、日降水量≥0.1 mm日數、平均相對濕度和極端最高、最低溫。與中上部煙葉相比，大田期的降水量對上部葉影響較為顯著。

表4 氣象因子與煙葉葉面密度關系的廣義可加模型分析結果

注：表中***表示在0.001水平下差異顯著，**表示在0.01水平下差異顯著，*表示在0.05水平下差異顯著，·表示在0.1水平下差異顯著。下同。

由表5可知，上部葉填充值與氣象因子平均氣壓和平均氣溫達到了P<0.01顯著水平，偏差解釋率分別達到了97.9%和66.9%；與平均最低氣溫、平均最高氣溫和平均水汽壓達到了P<0.05顯著水平，偏差解釋率分別達到了81.8%、67.7%和80.9%；與平均風速關系達到了P<0.1顯著水平，其偏差解釋率為23.8%。中部葉填充值與日降水量≥0.1 mm日數關系達到了P<0.05顯著水平，其偏差解釋率為59%；與平均風速關系達到了P<0.1顯著水平，其偏差解釋率為68.8%。下部葉填充值與極端最高氣溫關系達到了P<0.001顯著水平，其偏差解釋率為93.4%；與平均氣壓、平均水汽壓、平均最高氣溫達到了P<0.05顯著水平，三者的偏差解釋率均超過了60%，分別達到了70.1%、64.9%和68.5%；與平均氣溫和平均最低氣溫關系達到了P<0.1顯著水平，其偏差解釋率分別為61.3%和56%。綜合來看，煙葉填充值與氣象因子中的溫度(平均最高氣溫、最低氣溫、極端最高氣溫)、平均氣壓、水汽壓和風速關系密切；此外中部葉的填充值與日降水量≥0.1 mm日數關系密切，達到了顯著水平。

表5 氣象因子與煙葉填充值關系的廣義可加模型分析結果

表6結果表明：上部葉含梗率沒有與氣象因子達到顯著關系的指標，僅有平均最高氣溫偏差解釋率達到了53.5%；中部葉含梗率與日降水量≥0.1 mm日數關系達到了P<0.001顯著水平，其偏差解釋率為98.5%；與平均最高氣溫關系達到了P<0.05顯著水平，其偏差解釋率為76.7%；與極端最高氣溫關系達到了P<0.1顯著水平，其偏差解釋率為66.8%。此外，與日照時數、平均氣溫關系雖未達到顯著水平，但是其偏差解釋率分別達到了65.4%和64.1%，說明中部葉含梗率與這兩個氣象因子也具有較為密切的關系。下部葉含梗率與最大日降水量關系達到了P<0.05顯著水平，其偏差解釋率為32.6%；與極端最高氣溫關系達到了P<0.1顯著水平，其偏差解釋率為69.9%。此外，下部葉含梗率與平均風速、日照時數關系雖未達到顯著水平，但是其偏差解釋率分別達到了60.1%和67.7%，說明下部葉含梗率與這兩個氣象因素也有較為密切的關系。綜合來看，煙葉的含梗率與溫度尤其是高溫有密切的關系，同時日降水量≥0.1 mm日數對其也有較為顯著的影響。與中下部煙葉相比，上部葉含梗率受氣象因素影響的程度相對要弱一些。

表6 氣象因子與煙葉含梗率關系的廣義可加模型分析結果

由表7可知：上部葉平衡含水率與極端最高氣溫關系達到了P<0.05顯著水平，其偏差解釋率為29.7%；與平均氣溫和平均最高氣溫關系達到了P<0.1顯著水平，其偏差解釋率分別為25.1%和24.5%。中部葉平衡含水率與極端最高氣溫、極端最低氣溫、平均氣壓關系達到了P<0.05顯著水平，其偏差解釋率分別達到了57.2%、49.2%和31.1%；與日照百分率關系達到了P<0.1顯著水平，其偏差解釋率為73.1%；與日降水量≥0.1 mm日數關系雖未達到顯著水平，但是其偏差解釋率達到了76%。下部葉平衡含水率與日照百分率關系達到了P<0.01顯著水平，其偏差解釋率達到了53%；與降水量、平均氣壓、平均風速、平均水汽壓關系達到了P<0.05顯著水平，其偏差解釋率分別為48.1%、35.5%、14.2%和31%；與日照時數關系達到了P<0.1顯著水平，其偏差解釋率為66.7%。綜合看來，煙葉的平衡含水率與極端最高溫、平均氣壓、日照百分率、降水量等關系密切。

表7 氣象因子與煙葉平衡含水率關系的廣義可加模型分析結果

3 討論

本研究中，對煙葉葉面密度產生重要影響的氣象因素有日照百分率、日降水量≥0.1 mm日數、平均相對濕度、極端最高氣溫、極端最低氣溫；就部位而言，上部葉葉面密度對大田期的降水量比中下部煙葉更為敏感。說明煙葉的葉面密度主要受氣象因素日照、降水和溫度的影響，尤其是日照百分率對3個部位煙葉均有顯著影響，這也說明有效的光照是煙葉進行光合作用、有效積累干物質的必要條件。

對煙葉填充值產生重要影響的氣象因素有平均最高氣溫、平均最低氣溫、極端最高氣溫、平均氣壓、平均水汽壓和平均風速；就部位而言，中部葉的填充值與日降水量≥0.1 mm日數顯著相關，較上、下部葉對日降水量≥0.1 mm日數更為敏感。對含梗率有重要影響的氣象因素是極端最高氣溫、平均最高氣溫和平均溫度，說明溫度對煙葉含梗率較其它氣象因素影響更大；同時日降水量≥0.1 mm日數對中上部煙葉的含梗率也有顯著的影響；就部位而言，上部葉含梗率受氣象因素影響相對較弱。與煙葉平衡含水率關系密切的氣象因素主要是極端最高溫、平均氣壓、日照百分率、降水量；由于平衡含水率反映的是煙葉持水失水的能力，一方面受煙葉自身物理特性的影響，另一方面也是煙葉內在物質成分的反應，同時也與外在環境條件直接相關，因此，筆者認為煙葉平衡含水率受氣象條件的影響應該是通過影響煙葉內在成分的積累和物理特性的形成，進而間接影響到煙葉的平衡含水率。

雖然本文研究了氣象因素與煙葉物理特性間的關系，但是仍然存在幾個方面的問題：一是煙葉物理特性不僅包含本文研究的4項指標，仍有拉力、伸長率、抗張強度等相關指標，對于煙葉物理特性與氣象因素的研究仍待繼續深入；二是煙葉物理特性指標與氣象因素間的關系，雖然本文用GAM進行了分析，但是由于GAM模型自身特點(GAM模型不能給出具體參數估計模型)，對于氣象因素與煙葉物理特性指標的具體關系模型尚待繼續深入研究；三是基于煙草大田期葉片發育的時序性，大田期氣象因子的復雜性、階段性等因素的影響，針對煙草大田期氣象因素的階段細化研究仍待繼續深入。

[1] 付秋娟,張忠鋒,竇家宇,等.烤煙物理特性與常規化學成分及外觀質量的關系[J].中國煙草科學,2014,35(1):117-122.

[2] 王建民,孫意然,李瑞麗,等.烤煙物理特性與化學成分的典型相關分析[J].鄭州輕工業學院學報:自然科學版,2014,29(5):28-31.

[3] 王冬,趙銘欽,張學杰,等.烤煙物理特性與化學成分的相關及逐步回歸分析[J].中國農業大學學報,2010,15(6):52-58.

[4] 湯朝起,竇玉青,張俊.烤煙物理特性和化學成分與煙氣組分的關系[J].四川農業大學學報,2009,27(4):427-432.

[5] 薛超群,尹啟生,王廣山,等.烤煙煙葉物理特性的變化及其與評吸質量的關系[J].煙草科技,2008(7):52-55.

[6] 陳佳亮,黎娟,周清明,等.湖南濃香型烤煙物理特性與外觀質量關系分析[J].天津農業科學,2016,22(5):80-83.

[7] 周方舟,李志偉,周冀衡,等.海拔高度對陸良烤煙物理特性和化學品質的影響[J].湖南農業科學,2016(5):19-21.

[8] 李永智,楊虹琦,何偉,等.不同海拔高度氣象因素對烤煙物理特性的影響[J].湖南農業科學,2010(7):55-58.

[9] 吳祚友,衛盼盼,竇家宇,等.煙葉物理特性與打葉質量的關系[J].江西農業學報,2014,26(3):73-75.

[10] 衛盼盼,吳祚友,安銀立,等.煙葉物理特性與打葉風分工藝參數的關系[J].煙草科技,2014(8):5-9.

[11] 孟琦,于曉娜,李亞娟,等.河南初烤煙葉物理特性與土壤屬性的典型相關分析[J].河南農業大學學報,2014,48(4):413-419.

[12] 梁紅,張宇亭,楊超,等.不同移栽期對重慶烤煙生長和品質的影響[J].西南大學學報:自然科學版,2014,36(6):22-28.

[13] 徐敏,劉國順,劉小可.磷對不同烤煙品種調制后煙葉物理特性的影響[J].河南農業科學,2007(4):33-36.

[14] 荊聰瑩,杜詠梅,侯小東,等.烤煙填充值與主要理化指標的關系研究[J].中國農學通報,2013,29(30):202-205.

[15] 郭建華,宋紀真,王廣山,等.基于主成分分析和聚類分析的煙葉物理特性區域歸類[J].煙草科技,2014(8):14-17.

[16] 閆鐵軍,吳風光,王海明,等.烤煙物理特性差異分析及核心指標選擇[J].江西農業學報,2014,26(4):76-79.

[17] 劉琰琰,李海燕,陳超,等.攀西地區烤煙氣候適宜性評價指標建立及應用[J].四川農業大學學報,2015,33(3):299-305.

[18] 譚智勇,許龍,徐天養,等.文山煙區烤煙氣候適宜性評價及其變化特征分析[J].中國煙草科學,2015,36(1):78-83.

[19] 胡雪瓊,徐夢瑩,何雨芩,等.未來氣候變化對云南烤煙種植氣候適宜性的影響[J].應用生態學報,2016,27(4):1241-1247.

[20] 劉建福,尹大鋒,尹賓.YC/T 312—1996煙草及煙草制品試樣的制備和水分測定烘箱法[M].北京:中國標準出版社,2010.

(責任編輯：許晶晶)

Study on Relationship between Meteorological Factors and Physical Properties of Flue-cured Tobacco Leaf

GUO Dong-feng1, ZOU Peng1, LIU Fei1, WU Ke-song1, YAO Zhong-da1,SHU Jun-sheng1, DOU Yu-qing2, ZHANG Zhong-feng2, YAN Ning2

(1. Technology Center, China Tobacco Anhui Industrial Company Limited, Hefei 230088, China;2. Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao 266101, China)

Taking 3 different leaf-positions of tobacco leaf samples from 13 typical tobacco-planting areas of China as the studied objects, the author analyzed the relationship between field meteorological factors and physical properties of flue-cured tobacco leaves by using general additive model (GAM). The results indicated that the foliage density of tobacco leaves was influenced significantly by sunshine percentage, number of days with daily rainfall over 0.1 mm, average relative humidity, extreme maximum air temperature, and extreme minimum air temperature. Furthermore, the foliage density of upper leaves was more sensitive to the rainfall in tobacco growth period than that of middle and lower leaves. The average maximum air temperature, average minimum air temperature, extreme maximum air temperature, average air pressure, average vapor pressure, and average wind speed had an important impact on the filling power of tobacco leaves, and there was a significant correlation between the filling power of middle leaves and the number of days with daily rainfall over 0.1 mm. The extreme maximum air temperature, average maximum air temperature, and average air temperature had an important impact on the vein ratio of tobacco leaves, and the number of days with daily rainfall over 0.1 mm also had a significant effect on the vein ratio of upper and middle leaves. The vein ratio of upper leaves was affected more weakly by the meteorological factors than that of middle and lower leaves.

Physical property; Meteorological factor; Flue-cured tobacco

2016-08-29

中國煙草總公司科技重點項目“烤煙生產結構優化效應及關鍵技術研究與應用”(110210402007)；安徽中煙工業有限責 任公司科技計劃項目“皖南煙葉生產GAP管理模式研究”(2014124)、“皖南煙葉生產等級結構優化技術研究”(2014125)。

郭東鋒(1980—)，高級農藝師，博士，主要從事煙葉原料質量評價及生理生態研究。*通訊作者：閆寧。

S572

A

1001-8581(2017)02-0079-06