基于煙絲分布特征的卷煙硬度預(yù)測模型

何孝強 李泓燊 向 虎 王 龍 張云飛 鄒玉勝 趙 榮 劉戈弋楊 濤 李紹臣

(1. 紅云紅河煙草﹝集團﹞有限責(zé)任公司曲靖卷煙廠，云南 曲靖 655001；2. 紅云紅河煙草﹝集團﹞有限責(zé)任公司，云南 昆明 650000)

卷煙硬度是反映煙支飽滿程度和回彈性能的一項指標(biāo)[1]。硬度過低會造成卷煙輸送過程中外觀容易受損、消費者抽吸時發(fā)生燃燒錐掉落等，過高會造成卷煙吸阻、煙氣含量以及煙絲消耗增加。而且卷煙硬度能直接反映煙絲含量的多少，是制定卷煙重量標(biāo)準(zhǔn)的重要參考指標(biāo)[2]，但是在新品開發(fā)或者配方調(diào)整時，只能通過試生產(chǎn)出成品煙支后進行檢測獲得硬度指標(biāo)，再根據(jù)硬度調(diào)整重量，檢測存在嚴(yán)重的滯后性，不僅會造成較大的原輔料消耗，還存在一定的質(zhì)量隱患，因此在煙支卷制前對卷煙硬度進行預(yù)測對提高生產(chǎn)效率、降低消耗具有重要意義。

目前對卷煙硬度的研究主要包括：煙絲結(jié)構(gòu)對卷煙硬度的影響[3]，卷煙單支重量和圓周對卷煙硬度的影響[4-6]，煙支密度對卷煙硬度的影響[7]，膨脹梗絲填充能力對卷煙硬度的影響[8]等。以上研究內(nèi)容存在以下問題：① 僅對單一規(guī)格的卷煙進行硬度影響因素研究，當(dāng)卷煙規(guī)格發(fā)生變化時研究結(jié)果不再適用；② 分析了卷煙單重對硬度的影響關(guān)系，未將卷煙煙絲和卷煙材料分開，實際卷煙硬度測量的是煙絲部分；③ 卷煙硬度測量的是除去濾棒后的煙支中段位置，未考慮煙支的緊頭和中段煙絲量的區(qū)別，存在一定偏差；④ 僅做了卷煙硬度影響因素研究，未能對卷煙硬度做出準(zhǔn)確預(yù)測。試驗擬從卷煙硬度的檢測原理入手，通過對大量不同規(guī)格卷煙的硬度檢測數(shù)據(jù)的深入挖掘，研究煙支內(nèi)煙絲分布情況與卷煙硬度的關(guān)系，結(jié)合卷煙規(guī)格、煙絲特性和設(shè)備參數(shù)等因素建立卷煙硬度預(yù)測模型，以期根據(jù)模型預(yù)測結(jié)果指導(dǎo)卷煙重量標(biāo)準(zhǔn)的制定及相關(guān)參數(shù)的調(diào)整，提高生產(chǎn)效率。

1 原理分析

卷煙硬度是指在一定時間內(nèi)，試樣的徑向受到一定壓力，試樣受壓后與受壓前直徑的百分比[8]，即：

(1)

式中：

H——硬度，%；

a——壓陷量，mm；

D——壓縮前試樣直徑，mm；

d——壓縮后試樣直徑，mm。

標(biāo)準(zhǔn)[9]規(guī)定，卷煙硬度的測量位置為不包括濾嘴段的煙支中部位置。根據(jù)檢測原理，卷煙規(guī)格一定時，卷煙硬度與煙絲本身的填充力和實際的填充密度存在正相關(guān)關(guān)系[2]，由于在煙支卷制過程中，為了減少煙支空頭概率、提高搓接質(zhì)量，削減煙絲的平準(zhǔn)盤往往設(shè)計有凹槽，使煙支兩端煙絲密度高于中段，因此如果直接用卷煙的整體煙絲平均密度來預(yù)測硬度存在較大誤差，為建立精確的卷煙硬度預(yù)測模型，首先要計算出卷煙的中段煙絲含量。

煙支緊頭密度高于中段密度是由平準(zhǔn)盤凹槽削減過多的煙絲造成的，根據(jù)平準(zhǔn)盤凹槽的弧長及煙支長度可對煙支緊頭和中段進行分解，平準(zhǔn)盤削減瞬間的煙絲條經(jīng)卷煙紙包裹后存在一定的壓縮量，通過計算壓縮后煙支2個緊頭段和中段的體積即可計算出各段的煙絲量，具體計算步驟：

(1) 當(dāng)煙支緊頭位置和中段密度相同時，分別計算無嘴煙支中段、點燃端緊頭和接嘴端緊頭體積，其中煙支中段體積v計算公式：

(2)

式中：

π——圓周率；

v——煙支中段體積，mm3；

r——煙支半徑，mm；

l——無嘴煙支長度，mm；

l1——點燃端平準(zhǔn)盤凹槽弧長，mm；

l2——接嘴端平準(zhǔn)盤凹槽弧長(當(dāng)煙支只有點燃端一個緊頭時，l2為0)，mm。

煙支點燃端緊頭體積v1計算公式：

(3)

式中：

v1——煙支點燃端緊頭體積，mm3；

l1——點燃端平準(zhǔn)盤凹槽弧長，mm；

r——煙支半徑，mm；

h1——點燃端平準(zhǔn)盤凹槽徑深，mm；

d1——點燃端平準(zhǔn)盤凹槽深度，mm。

煙支接嘴端緊頭體積v2計算公式：

(4)

式中：

v2——煙支接嘴端緊頭體積，mm3；

l2——接嘴端平準(zhǔn)盤凹槽弧長，mm；

r——煙支半徑，mm；

h2——接嘴端平準(zhǔn)盤凹槽徑深，mm；

d2——接嘴端平準(zhǔn)盤凹槽深度，mm。

(2) 計算煙支中段煙絲量。首先根據(jù)預(yù)設(shè)卷煙重量和已知的卷煙材料重量計算出煙絲量，然后根據(jù)煙支各段體積計算出中段煙絲量m，具體計算公式：

(5)

式中：

m——煙支中段煙絲量，g；

m1——預(yù)設(shè)單支卷煙重量，g；

m2——平均單支卷煙輔料重量，g；

v——煙支緊頭和中段密度相同時的煙支中段體積，mm3；

v1——點燃端緊頭體積，mm3；

v2——接嘴端緊頭體積，mm3。

2 模型建立與驗證

2.1 煙支中段填充系數(shù)的計算

選擇曲靖卷煙廠2018年1月～12月的15個常規(guī)牌號卷煙物理指標(biāo)檢測數(shù)據(jù)，數(shù)據(jù)包括卷煙重量、圓周、長度、硬度，卷煙規(guī)格包括常規(guī)卷煙和中、細支卷煙，涵蓋了目前主流的卷煙規(guī)格。經(jīng)數(shù)據(jù)清洗后共167 721條數(shù)據(jù)，每個樣本為一支卷煙的物理指標(biāo)檢測數(shù)據(jù)，根據(jù)上述原理在數(shù)據(jù)集中加入計算煙支中段煙絲量所需要的中間變量，包括卷煙材料重量、濾棒長度、煙絲填充值和平準(zhǔn)盤凹槽規(guī)格，計算出每個樣本的煙支中段煙絲量m，結(jié)合煙支規(guī)格和平準(zhǔn)盤凹槽規(guī)格計算出煙支中段煙絲密度，由于卷煙硬度受到煙絲填充值和填充密度的綜合影響，因此用煙支中段的煙絲密度乘以煙絲填充值得到表征煙支中段煙絲填充情況的變量，記作煙支中段填充系數(shù)fdm，計算公式為：

(6)

式中：

fdm——煙支中段填充系數(shù)；

f——煙絲填充值，cm3/g；

m——煙支中段煙絲量，g；

v——煙支中段體積，mm3。

根據(jù)式(6)計算出每個樣本的fdm，加入數(shù)據(jù)集中。

2.2 數(shù)據(jù)處理

為探索相關(guān)因素對卷煙硬度的影響模式，以原始數(shù)據(jù)集的卷煙牌號為分類變量，以每個類別的均值進行初步探索，得出中段填充系數(shù)fdm與硬度檢測結(jié)果的散點圖(圖1)。

圖1 fdm與煙支硬度散點圖

從圖1可以看出，煙支中段填充系數(shù)與卷煙硬度沒有明顯的關(guān)系，由于硬度是壓縮前后卷煙直徑的比例，相同硬度的不同規(guī)格卷煙壓縮量不同，因此考慮將卷煙硬度H根據(jù)式(1)轉(zhuǎn)換為測量硬度時的徑向壓縮量a進行分析，得出壓縮量與煙支中段填充系數(shù)的散點圖(圖2)。

圖2 fdm與壓縮量散點圖

圖與壓縮量散點圖

2.3 模型建立與驗證

首先將數(shù)據(jù)集分為訓(xùn)練集和測試集，其中訓(xùn)練集為總樣本的70%，采用分層抽樣的方法在每個牌號卷煙的檢測數(shù)據(jù)中隨機抽取70%的樣本量組成訓(xùn)練集建立模型，剩下的30%為作為驗證模型的測試集。為盡可能消除檢測誤差對模型的影響，以fdm和D為變量，利用K-Means聚類算法[10]對測試集進行聚類，首先根據(jù)總的類內(nèi)平方和來選擇最佳的聚類個數(shù)，計算結(jié)果見圖4。

圖4 聚類數(shù)與總的類內(nèi)平方和散點圖

Figure 4 Scatter plot of the number of clusters and total within-cluster sum of squares

(7)

回歸方差的擬合優(yōu)度R2>98%，擬合效果較好。

預(yù)測出煙支徑向壓縮量后，按式(1)計算卷煙硬度。

利用驗證集對模型進行驗證，通過計算，模型預(yù)測的相對誤差平均值為3.7%，測試集的標(biāo)準(zhǔn)化均方誤差[11]為0.6，預(yù)測效果較好。為進一步分析模型的預(yù)測誤差，將各卷煙牌號硬度指標(biāo)在生產(chǎn)中的實際波動和產(chǎn)品技術(shù)標(biāo)準(zhǔn)與模型預(yù)測誤差進行對比，結(jié)果見表1。

表1 模型預(yù)測誤差對比表

從表1可以看出，模型對各不同規(guī)格卷煙的硬度指標(biāo)預(yù)測誤差遠遠小于指標(biāo)正常波動范圍及企業(yè)產(chǎn)品技術(shù)標(biāo)準(zhǔn)規(guī)定的允差范圍，說明模型可應(yīng)用于實際生產(chǎn)過程中的卷煙硬度預(yù)測，且應(yīng)用方便，不受卷煙規(guī)格、設(shè)備型號、煙絲特性等因素限制，適用范圍廣。

3 結(jié)論

試驗研究了煙支中段填充系數(shù)與測量卷煙硬度時的徑向壓縮量的關(guān)系，建立了壓縮量與煙絲填充系數(shù)和煙支直徑的回歸模型，進而建立卷煙硬度預(yù)測模型對卷煙硬度進行預(yù)測，模型擬合優(yōu)度>98%、測試集的標(biāo)準(zhǔn)化均方誤差為0.6，相對誤差均值為3.7%，模型預(yù)測誤差遠遠小于指標(biāo)正常波動范圍及產(chǎn)品技術(shù)標(biāo)準(zhǔn)規(guī)定的允差范圍，說明模型預(yù)測效果較好，同時模型適用于各種規(guī)格的卷煙，可以應(yīng)用于實際生產(chǎn)中的卷煙硬度預(yù)測，為制定卷煙重量標(biāo)準(zhǔn)和產(chǎn)品質(zhì)量改進提供參考。下一步將加入更多規(guī)格牌號的卷煙數(shù)據(jù)對模型進一步優(yōu)化和完善。