• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    調(diào)控主流煙氣中巴豆醛釋放量的制絲工藝參數(shù)模型

    2014-07-02 11:09:00袁月等
    湖北農(nóng)業(yè)科學(xué) 2014年6期
    關(guān)鍵詞:卷煙

    袁月等

    摘要:為了確定制絲工藝對(duì)主流煙氣中巴豆醛釋放量的影響,通過均勻設(shè)計(jì)法,建立了切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度、筒壁溫度、筒體轉(zhuǎn)速6種制絲工藝參數(shù)與卷煙主流煙氣巴豆醛釋放量的二次多項(xiàng)式逐步回歸模型,并對(duì)模型進(jìn)行了因素主效應(yīng)分析、單因素效應(yīng)分析、邊際效應(yīng)分析和雙因素交互作用分析。結(jié)果表明,6種工藝參數(shù)對(duì)巴豆醛釋放量的貢獻(xiàn)率由大到小依次為筒壁溫度、熱風(fēng)風(fēng)門開度、切絲寬度、熱風(fēng)溫度、筒體轉(zhuǎn)速、HT工作蒸氣壓力。切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度的增加會(huì)促進(jìn)巴豆醛的釋放,筒體轉(zhuǎn)速和筒壁溫度的升高會(huì)抑制巴豆醛的釋放;切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、筒壁溫度和筒體轉(zhuǎn)速對(duì)巴豆醛釋放量的影響隨其他工藝參數(shù)水平的升高而增大,熱風(fēng)風(fēng)門開度的影響不變。切絲寬度與HT工作蒸氣壓力、切絲寬度與熱風(fēng)溫度、筒壁溫度與筒體轉(zhuǎn)速之間存在交互作用影響巴豆醛釋放量。

    關(guān)鍵詞:卷煙;制絲工藝參數(shù);巴豆醛;均勻設(shè)計(jì);二次多項(xiàng)式逐步回歸

    中圖分類號(hào):TS452+.3 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):0439-8114(2014)06-1322-05

    Models of Processing Parameters for Controlling the Crotonaldehyde in Cigarette Mainstream Smoke

    YUAN Yue,TAN Lan-lan,DAI Ya,F(xiàn)ENG Guang-lin,YANG Wen-min,CHEN Kun-yan,XIAO Ke-yi,WANG Chang-guo (Technical Research Center of China Tobacco Chuanyu Industrial Limited Liability Corporation/Sichuan key Laboratory for Harmful Components and Tar Reduction in Cigarette, Chengdu 610066, China)

    Abstract: A mathematical model was established by regression techniques from the six processing parameters including tobacco cut width, HT operating steam pressure, hot-air temperature, hot air damper opening, the temperature of cylinder wall, and the revolving speed of casing cylinder to the emission of crotonaldehyde in mainstream smoke. The principal-factor effect, single-factor effect, marginal effect and double-factor effect of the model were analyzed. The results showed that the contribution order for the crotonaldehyde emission was the temperature of cylinder wall>hot air damper opening>tobacco cut width>hot-air temperature>the revolving speed of casing cylinder>HT operating steam pressure. The crotonaldehyde release amount increased with tobacco cut width, HT operating steam pressure, hot-air temperature and hot air damper opening, and decreased with the temperature of cylinder wall. The revolving speed of casing cylinder promoted the crotonaldehyde delivery under low level background, and inhibited it under middle and high level background. The influence of tobacco cut width, HT operating steam pressure, hot-air temperature, temperature of cylinder wall and the revolving speed of casing cylinder on crotonaldehyde emission aggravated with the increasing of other parameters. Hot air damper openings effect didnt change. In addition, interactive influence on crotonaldehyde delivery existed between tobacco cut width and HT operating steam pressure, tobacco cut width and the temperature of cylinder wall, the temperature of cylinder wall and the revolving speed of casing cylinder.

    Key words: cigarette; processing parameters; crotonaldehyde; uniform design; quadratic polyriomial stepwise regression

    巴豆醛即丁烯醛,屬于卷煙主流煙氣中主要揮發(fā)性羰基化合物中的一種,是呼吸道纖毛毒素,對(duì)眼結(jié)膜及上呼吸道黏膜有強(qiáng)烈刺激作用,長(zhǎng)期接觸易引起慢性鼻炎、神經(jīng)系統(tǒng)機(jī)能障礙[1,2],被收錄在Hoffmann等[3]所列的卷煙煙氣危害性成分名單中,謝劍平等[4]將其列為7種卷煙主流煙氣危害性指標(biāo)之一。因此,降低卷煙主流煙氣中巴豆醛的釋放量是當(dāng)前卷煙減害研究的重要組成部分。國(guó)內(nèi)關(guān)于卷煙制絲工藝參數(shù)對(duì)煙氣有害成分釋放量的影響有過一些報(bào)道[5-7],但在試驗(yàn)設(shè)計(jì)上主要采用單因素循環(huán)的方法,雖有組合性的研究,所建立的模型卻存在一定的局限性,不能有效反映單一因素的影響和兩兩因素的交互作用。

    本研究選取卷煙生產(chǎn)中涉及的3個(gè)關(guān)鍵制絲工序的6個(gè)工藝參數(shù),即切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度、筒壁溫度及筒體轉(zhuǎn)速,采用均勻設(shè)計(jì)開展試驗(yàn),以各工藝參數(shù)為自變量,巴豆醛釋放量為因變量,應(yīng)用二次多項(xiàng)式逐步回歸方法,建立了工藝參數(shù)到巴豆醛釋放量的模型,并對(duì)模型進(jìn)行了深入分析,旨在為卷煙實(shí)際生產(chǎn)中工藝參數(shù)的選擇提供指導(dǎo),從而有效調(diào)控卷煙主流煙氣中巴豆醛的含量。

    1 材料與方法

    1.1 材料與設(shè)備

    1.1.1 材料 所用材料為川渝中煙涪陵分廠生產(chǎn)某一品牌配方煙絲。

    1.1.2 設(shè)備 工藝設(shè)備:FT112微波松散機(jī)(江蘇智思機(jī)械制造有限責(zé)任公司)、WQ3313滾筒式葉片回潮機(jī)(昆明船舶制造有限責(zé)任公司)、SQ313切葉絲機(jī)(昆明船舶制造有限責(zé)任公司)、隧道式葉絲回潮機(jī)(昆明船舶制造有限責(zé)任公司)、SH315D滾筒管板式烘絲機(jī)(秦皇島煙草機(jī)械有限責(zé)任公司);檢測(cè)設(shè)備:Agilent 1200高效液相色譜儀和Agilent 1200 FLD檢測(cè)器;20孔道轉(zhuǎn)盤式吸煙機(jī)(德國(guó)BORGWALDT公司)。

    1.2 方法

    1.2.1 樣品制備 按照隨機(jī)化原則選擇試驗(yàn)順序。在設(shè)備運(yùn)行平穩(wěn)及保持烘后葉絲含水率符合工藝要求的前提下,按文獻(xiàn)[8]中規(guī)定的方法在烘絲機(jī)出口處分別取樣5次,每次10 kg?;旌暇鶆蚝笥盟姆址ūA魳悠窡熃z約20 kg。樣品煙絲經(jīng)自然平衡含水率在11.8%~12.5%后,利用同一機(jī)器用相同卷接材料卷制,煙支的圓周、硬度及單支重量均在規(guī)定的允許誤差范圍內(nèi)。卷煙樣品平衡48 h經(jīng)重量和吸阻挑選后,用于巴豆醛釋放量檢測(cè)。

    1.2.2 檢測(cè)方法 按照標(biāo)準(zhǔn)YC/T 255-2008[9]對(duì)樣品卷煙主流煙氣中的巴豆醛含量進(jìn)行測(cè)定。

    2 結(jié)果與分析

    2.1 制絲實(shí)際工藝參數(shù)與試驗(yàn)結(jié)果

    選取卷煙生產(chǎn)中的切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)、熱風(fēng)風(fēng)門開度(x4)、筒壁溫度(x5)、筒體轉(zhuǎn)速(x6)作為主要工藝參數(shù)進(jìn)行均勻設(shè)計(jì)[10]試驗(yàn),共7個(gè)組合,3次重復(fù),各因素范圍見表1。對(duì)試驗(yàn)中各工藝參數(shù)按0~1的范圍進(jìn)行極差歸一化處理,排除量綱的影響,結(jié)果如表2所示。

    對(duì)SDH-1、SDH-2、SDH-3、SDH-4、SDH-5、SDH-6、SDH-7共7組卷煙樣品主流煙氣中的巴豆醛釋放量檢測(cè)結(jié)果進(jìn)行描述性統(tǒng)計(jì)分析,結(jié)果顯示,巴豆醛釋放量為16.96~19.24 μg/支,平均為17.96 μg/支,中位數(shù)為17.70 μg/支,峰度為0.40,偏度為0.58。

    2.2 模型的建立

    選用二次多項(xiàng)式逐步回歸分析方法建立模型,一般公式為:y=b0+bijxixj,以歸一化后的工藝參數(shù)(表2)作為自變量,巴豆醛釋放量作為因變量,得到工藝參數(shù)到巴豆醛釋放量的模型,方程為:y=17.500+1.480x42+0.059x52+0.578x1x2+1.930 x1x3-2.110x5x6。

    同時(shí)得到回歸方程的復(fù)相關(guān)系數(shù)R=0.999 9,P=0.002 5<0.05??梢钥闯鲈撃P托Ч^好,能夠正確反映各工藝參數(shù)與巴豆醛釋放量之間的關(guān)系,可靠性較高。

    2.3 因素主效應(yīng)分析

    采用因素主效應(yīng)分析法[11]計(jì)算模型中各因素的貢獻(xiàn)率,考察各工藝參數(shù)對(duì)巴豆醛釋放量的影響。根據(jù)本試驗(yàn)建立的模型計(jì)算各因素的貢獻(xiàn)率x1為0.992,x2為0.495,x3為0.498,x4為0.998,x5為1.418,x6為0.497,表明6種工藝參數(shù)對(duì)巴豆醛釋放量影響的大小順序?yàn)橥脖跍囟?、熱風(fēng)風(fēng)門開度、切絲寬度、熱風(fēng)溫度、筒體轉(zhuǎn)速、HT工作蒸氣壓力。

    2.4 單因素效應(yīng)分析

    依據(jù)表2歸一化后的工藝參數(shù),分別設(shè)置高、中、低3個(gè)水平,具體見表3。在其他因素置于高、中、低水平時(shí),考察了回歸模型中單一工藝參數(shù)變化對(duì)巴豆醛釋放量的影響,得到的一元函數(shù)方程如表4所示。從表4中可以看出,將其他因素分別置于高水平和中水平時(shí),切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)、筒體轉(zhuǎn)速(x6)4個(gè)因素的效應(yīng)方程都是一元一次方程,說明其與巴豆醛釋放量呈線性相關(guān),其中,切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)為正相關(guān),筒體轉(zhuǎn)速(x6)為負(fù)相關(guān);將其他因素置于低水平時(shí),這4個(gè)因素的效應(yīng)方程均為常數(shù)方程,說明其對(duì)巴豆醛釋放量無影響。當(dāng)熱風(fēng)風(fēng)門開度(x4)水平為0~1,其他因素分別置于高、中、低水平時(shí),隨著熱風(fēng)風(fēng)門開度(x4)的增大,巴豆醛釋放量(y)呈拋物線形式逐漸增大,升高趨勢(shì)先慢后快。當(dāng)筒壁溫度(x5)水平為0~1時(shí),其他因素分別置于高、中、低水平時(shí),巴豆醛釋放量(y)均隨著筒壁溫度(x5)的增加而下降,且各水平下的下降速率變化不明顯,巴豆醛釋放量(y)的下降速率由大到小依次為高水平、中水平、低水平。

    2.5 邊際效應(yīng)分析

    通過分析模型中各因素的邊際效應(yīng),可以進(jìn)一步確定巴豆醛釋放量隨各工藝參數(shù)改變的變化速率和變化方向。根據(jù)建立的模型,求出各因素的一階偏導(dǎo)數(shù),按照因素與水平表(表3),將其他因素分別置于低、中、高水平,求出各因素分別在高、中、低水平下巴豆醛的邊際釋放量(y),結(jié)果依次如表5、表6、表7所示。由表5、表6、表7可知,工藝參數(shù)切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)、筒體轉(zhuǎn)速(x6)分別處于高、中、低水平時(shí),其對(duì)應(yīng)各水平的巴豆醛邊際釋放量不變,即其對(duì)巴豆醛釋放量影響的變化速率不變。且將其他因素置于低水平(表5)時(shí),x1、x2、x3、x6的巴豆醛邊際釋放量均為0,說明這4種工藝參數(shù)對(duì)巴豆醛釋放量無影響。將其他因素置于中水平(表6)、高水平(表7)時(shí),x1、x2、x3的巴豆醛邊際釋放量為正值,說明巴豆醛釋放量會(huì)隨其增加而增加;x6的巴豆醛邊際釋放量為負(fù)值,說明巴豆醛釋放量會(huì)隨其增加而降低。這都與單因素效應(yīng)分析結(jié)果一致。

    熱風(fēng)風(fēng)門開度(x4)的巴豆醛邊際釋放量為正值,且不因其他因素水平的改變而變化,僅隨其自身因素水平的降低而減小,說明該工藝參數(shù)的增加會(huì)促進(jìn)巴豆醛的釋放,且影響程度也會(huì)相應(yīng)增大。筒壁溫度(x5)的邊際效應(yīng)比較復(fù)雜,當(dāng)其他因素置于低水平時(shí),x5的巴豆醛邊際釋放量為正值,并隨其自身水平的提高而變大,說明巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而增加,影響速率也逐漸增大;但當(dāng)其他因素分別置于中水平和高水平時(shí),其巴豆醛邊際釋放量為負(fù)值,且絕對(duì)值隨其自身水平的提高而減小,說明此情況下巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而降低,受影響程度也逐漸減小。

    綜合以上結(jié)果來看,除了熱風(fēng)風(fēng)門開度(x4),其余因素的巴豆醛邊際釋放量絕對(duì)值均隨其他因素水平的升高而增大,說明其他因素的水平越高,巴豆醛釋放量受這5種工藝參數(shù)變化的影響也越大。

    2.6 雙因素交互作用分析

    采用等高線分析法[12]考察建立的回歸模型中任意兩種工藝參數(shù)對(duì)巴豆醛釋放量的交互作用影響結(jié)果。選定模型中的2個(gè)因素,將其他因素按照表4分別置于高、中、低水平,共得到等高線方程42個(gè),結(jié)果只有切絲寬度(x1)與HT工作蒸氣壓力(x2)、切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線方程中含有交互項(xiàng)x1x2、x1x3、x5x6,如表8所示。表明對(duì)于本試驗(yàn)建立的模型,只有x1與x2、x1與x3、x5與x6之間存在交互作用。根據(jù)其等高線方程作圖發(fā)現(xiàn),這3對(duì)因素的等高線圖基本變化趨勢(shì)相同,即隨著2種工藝參數(shù)的減小,巴豆醛釋放量逐漸減小。以工藝參數(shù)切絲寬度(x1)與HT工作蒸汽壓力(x2)的等高線圖為例,當(dāng)其他因素分別置于低、中、高水平時(shí),其等高線圖如圖1所示。切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線圖與圖1a一致,等高線疏密分布均勻,說明2種因素對(duì)巴豆醛釋放量的影響程度沒有明顯差別。從圖1b、1c可以看出,當(dāng)其他工藝參數(shù)置于中水平和高水平時(shí),切絲寬度(x1)與HT工作蒸氣壓力(x2)的交互作用表現(xiàn)為,沿切絲寬度(x1)軸向等高線變化密集,而HT工作蒸氣壓力(x2)軸向等高線變化相對(duì)稀疏,說明切絲寬度(x1)對(duì)巴豆醛釋放量的影響比HT工作蒸氣壓力(x2)大。

    3 結(jié)論

    通過均勻設(shè)計(jì)安排試驗(yàn),采用二次多項(xiàng)式逐步回歸方法建立了包括切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度、筒壁溫度、筒體轉(zhuǎn)速6個(gè)工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的模型。模型的因素主效應(yīng)分析結(jié)果表明,各工藝參數(shù)對(duì)巴豆醛釋放量影響的貢獻(xiàn)率由大到小依次為筒壁溫度、熱風(fēng)風(fēng)門開度、切絲寬度、熱風(fēng)溫度、筒體轉(zhuǎn)速、HT工作蒸氣壓力。單因素效應(yīng)分析發(fā)現(xiàn),切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度與巴豆醛釋放量呈線性正相關(guān),筒體轉(zhuǎn)速則與巴豆醛釋放量呈線性負(fù)相關(guān),而熱風(fēng)風(fēng)門開度和筒壁溫度與巴豆醛釋放量呈較復(fù)雜的拋物線相關(guān)。進(jìn)一步的邊際效應(yīng)分析證實(shí),在0~1水平下,巴豆醛釋放量隨切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度的增大而提高,隨筒體轉(zhuǎn)速和筒壁溫度的升高而降低;且切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、筒壁溫度和筒體轉(zhuǎn)速對(duì)巴豆醛釋放量影響的變化速率隨其他工藝參數(shù)的增大而變大,熱風(fēng)風(fēng)門開度的影響不變。雙因素交互作用的分析結(jié)果表明,只有切絲寬度與HT工作蒸氣壓力、切絲寬度與熱風(fēng)溫度、筒壁溫度與筒體轉(zhuǎn)速之間存在交互作用。本研究建立的模型對(duì)于從工藝參數(shù)的角度控制卷煙主流煙氣中巴豆醛的釋放量具有一定的指導(dǎo)意義,在卷煙實(shí)際生產(chǎn)中還需要考慮到其他因素如卷煙材料、葉組配方等的影響。

    參考文獻(xiàn):

    [1] MANNING D L, ASKERINEC M P, JENKINS R A,et al. High performance liquid chromatographic determinations of selected gas phase carbonyls in tobacco smoke[J]. Journal of Association of Analytical Chemistry,1983,66(1): 8-12.

    [2] WYNDER E L,HOFFMAN D. Tobacco and tobacco smoke[M]. New York: Academic Press,1967.

    [3] HOFFMANN I, HOFFMANN D. The Changing Cigarette: Chemical Studies and Bioassays[M]. New York:Oxford University Press,2001.

    [4] 謝劍平,劉惠民,朱茂祥,等.卷煙煙氣危害性指數(shù)研究[J]. 煙草科技,2009(2):5-15.

    [5] 王 鵬,劉 華,曾 建,等.制絲工藝參數(shù)對(duì)主流煙氣中氫氰酸含量影響研究[J].湖北農(nóng)業(yè)科學(xué),2011,50(2):365-367.

    [6] 譚蘭蘭,施豐成,薛 芳,等.卷煙制絲工藝參數(shù)對(duì)主流煙氣中苯并[a]芘釋放量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(1):165-167.

    [7] 陳昆燕,薛 芳,曾 建,等.重點(diǎn)工序工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的關(guān)系研究[J].江西農(nóng)業(yè)學(xué)報(bào),2009,21(12):40-42.

    [8] 秦前浩.卷煙工藝測(cè)試與分析大綱[M]. 成都:四川大學(xué)出版社,2004.

    [9] YC/T 254-2008,卷煙主流煙氣中主要羰基化合物的測(cè)定 高效液相色譜法[S].

    [10] 方開泰.均勻設(shè)計(jì)與均勻設(shè)計(jì)表[M].北京:科學(xué)出版社,1994.

    [11] 唐啟義,馮明光.實(shí)用統(tǒng)計(jì)分析及其DPS數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [12] 何 為,薛衛(wèi)東,唐 斌.優(yōu)化試驗(yàn)設(shè)計(jì)方法及數(shù)據(jù)分析[M].北京:化學(xué)工業(yè)出版社,2012.

    2.5 邊際效應(yīng)分析

    通過分析模型中各因素的邊際效應(yīng),可以進(jìn)一步確定巴豆醛釋放量隨各工藝參數(shù)改變的變化速率和變化方向。根據(jù)建立的模型,求出各因素的一階偏導(dǎo)數(shù),按照因素與水平表(表3),將其他因素分別置于低、中、高水平,求出各因素分別在高、中、低水平下巴豆醛的邊際釋放量(y),結(jié)果依次如表5、表6、表7所示。由表5、表6、表7可知,工藝參數(shù)切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)、筒體轉(zhuǎn)速(x6)分別處于高、中、低水平時(shí),其對(duì)應(yīng)各水平的巴豆醛邊際釋放量不變,即其對(duì)巴豆醛釋放量影響的變化速率不變。且將其他因素置于低水平(表5)時(shí),x1、x2、x3、x6的巴豆醛邊際釋放量均為0,說明這4種工藝參數(shù)對(duì)巴豆醛釋放量無影響。將其他因素置于中水平(表6)、高水平(表7)時(shí),x1、x2、x3的巴豆醛邊際釋放量為正值,說明巴豆醛釋放量會(huì)隨其增加而增加;x6的巴豆醛邊際釋放量為負(fù)值,說明巴豆醛釋放量會(huì)隨其增加而降低。這都與單因素效應(yīng)分析結(jié)果一致。

    熱風(fēng)風(fēng)門開度(x4)的巴豆醛邊際釋放量為正值,且不因其他因素水平的改變而變化,僅隨其自身因素水平的降低而減小,說明該工藝參數(shù)的增加會(huì)促進(jìn)巴豆醛的釋放,且影響程度也會(huì)相應(yīng)增大。筒壁溫度(x5)的邊際效應(yīng)比較復(fù)雜,當(dāng)其他因素置于低水平時(shí),x5的巴豆醛邊際釋放量為正值,并隨其自身水平的提高而變大,說明巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而增加,影響速率也逐漸增大;但當(dāng)其他因素分別置于中水平和高水平時(shí),其巴豆醛邊際釋放量為負(fù)值,且絕對(duì)值隨其自身水平的提高而減小,說明此情況下巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而降低,受影響程度也逐漸減小。

    綜合以上結(jié)果來看,除了熱風(fēng)風(fēng)門開度(x4),其余因素的巴豆醛邊際釋放量絕對(duì)值均隨其他因素水平的升高而增大,說明其他因素的水平越高,巴豆醛釋放量受這5種工藝參數(shù)變化的影響也越大。

    2.6 雙因素交互作用分析

    采用等高線分析法[12]考察建立的回歸模型中任意兩種工藝參數(shù)對(duì)巴豆醛釋放量的交互作用影響結(jié)果。選定模型中的2個(gè)因素,將其他因素按照表4分別置于高、中、低水平,共得到等高線方程42個(gè),結(jié)果只有切絲寬度(x1)與HT工作蒸氣壓力(x2)、切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線方程中含有交互項(xiàng)x1x2、x1x3、x5x6,如表8所示。表明對(duì)于本試驗(yàn)建立的模型,只有x1與x2、x1與x3、x5與x6之間存在交互作用。根據(jù)其等高線方程作圖發(fā)現(xiàn),這3對(duì)因素的等高線圖基本變化趨勢(shì)相同,即隨著2種工藝參數(shù)的減小,巴豆醛釋放量逐漸減小。以工藝參數(shù)切絲寬度(x1)與HT工作蒸汽壓力(x2)的等高線圖為例,當(dāng)其他因素分別置于低、中、高水平時(shí),其等高線圖如圖1所示。切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線圖與圖1a一致,等高線疏密分布均勻,說明2種因素對(duì)巴豆醛釋放量的影響程度沒有明顯差別。從圖1b、1c可以看出,當(dāng)其他工藝參數(shù)置于中水平和高水平時(shí),切絲寬度(x1)與HT工作蒸氣壓力(x2)的交互作用表現(xiàn)為,沿切絲寬度(x1)軸向等高線變化密集,而HT工作蒸氣壓力(x2)軸向等高線變化相對(duì)稀疏,說明切絲寬度(x1)對(duì)巴豆醛釋放量的影響比HT工作蒸氣壓力(x2)大。

    3 結(jié)論

    通過均勻設(shè)計(jì)安排試驗(yàn),采用二次多項(xiàng)式逐步回歸方法建立了包括切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度、筒壁溫度、筒體轉(zhuǎn)速6個(gè)工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的模型。模型的因素主效應(yīng)分析結(jié)果表明,各工藝參數(shù)對(duì)巴豆醛釋放量影響的貢獻(xiàn)率由大到小依次為筒壁溫度、熱風(fēng)風(fēng)門開度、切絲寬度、熱風(fēng)溫度、筒體轉(zhuǎn)速、HT工作蒸氣壓力。單因素效應(yīng)分析發(fā)現(xiàn),切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度與巴豆醛釋放量呈線性正相關(guān),筒體轉(zhuǎn)速則與巴豆醛釋放量呈線性負(fù)相關(guān),而熱風(fēng)風(fēng)門開度和筒壁溫度與巴豆醛釋放量呈較復(fù)雜的拋物線相關(guān)。進(jìn)一步的邊際效應(yīng)分析證實(shí),在0~1水平下,巴豆醛釋放量隨切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度的增大而提高,隨筒體轉(zhuǎn)速和筒壁溫度的升高而降低;且切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、筒壁溫度和筒體轉(zhuǎn)速對(duì)巴豆醛釋放量影響的變化速率隨其他工藝參數(shù)的增大而變大,熱風(fēng)風(fēng)門開度的影響不變。雙因素交互作用的分析結(jié)果表明,只有切絲寬度與HT工作蒸氣壓力、切絲寬度與熱風(fēng)溫度、筒壁溫度與筒體轉(zhuǎn)速之間存在交互作用。本研究建立的模型對(duì)于從工藝參數(shù)的角度控制卷煙主流煙氣中巴豆醛的釋放量具有一定的指導(dǎo)意義,在卷煙實(shí)際生產(chǎn)中還需要考慮到其他因素如卷煙材料、葉組配方等的影響。

    參考文獻(xiàn):

    [1] MANNING D L, ASKERINEC M P, JENKINS R A,et al. High performance liquid chromatographic determinations of selected gas phase carbonyls in tobacco smoke[J]. Journal of Association of Analytical Chemistry,1983,66(1): 8-12.

    [2] WYNDER E L,HOFFMAN D. Tobacco and tobacco smoke[M]. New York: Academic Press,1967.

    [3] HOFFMANN I, HOFFMANN D. The Changing Cigarette: Chemical Studies and Bioassays[M]. New York:Oxford University Press,2001.

    [4] 謝劍平,劉惠民,朱茂祥,等.卷煙煙氣危害性指數(shù)研究[J]. 煙草科技,2009(2):5-15.

    [5] 王 鵬,劉 華,曾 建,等.制絲工藝參數(shù)對(duì)主流煙氣中氫氰酸含量影響研究[J].湖北農(nóng)業(yè)科學(xué),2011,50(2):365-367.

    [6] 譚蘭蘭,施豐成,薛 芳,等.卷煙制絲工藝參數(shù)對(duì)主流煙氣中苯并[a]芘釋放量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(1):165-167.

    [7] 陳昆燕,薛 芳,曾 建,等.重點(diǎn)工序工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的關(guān)系研究[J].江西農(nóng)業(yè)學(xué)報(bào),2009,21(12):40-42.

    [8] 秦前浩.卷煙工藝測(cè)試與分析大綱[M]. 成都:四川大學(xué)出版社,2004.

    [9] YC/T 254-2008,卷煙主流煙氣中主要羰基化合物的測(cè)定 高效液相色譜法[S].

    [10] 方開泰.均勻設(shè)計(jì)與均勻設(shè)計(jì)表[M].北京:科學(xué)出版社,1994.

    [11] 唐啟義,馮明光.實(shí)用統(tǒng)計(jì)分析及其DPS數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [12] 何 為,薛衛(wèi)東,唐 斌.優(yōu)化試驗(yàn)設(shè)計(jì)方法及數(shù)據(jù)分析[M].北京:化學(xué)工業(yè)出版社,2012.

    2.5 邊際效應(yīng)分析

    通過分析模型中各因素的邊際效應(yīng),可以進(jìn)一步確定巴豆醛釋放量隨各工藝參數(shù)改變的變化速率和變化方向。根據(jù)建立的模型,求出各因素的一階偏導(dǎo)數(shù),按照因素與水平表(表3),將其他因素分別置于低、中、高水平,求出各因素分別在高、中、低水平下巴豆醛的邊際釋放量(y),結(jié)果依次如表5、表6、表7所示。由表5、表6、表7可知,工藝參數(shù)切絲寬度(x1)、HT工作蒸氣壓力(x2)、熱風(fēng)溫度(x3)、筒體轉(zhuǎn)速(x6)分別處于高、中、低水平時(shí),其對(duì)應(yīng)各水平的巴豆醛邊際釋放量不變,即其對(duì)巴豆醛釋放量影響的變化速率不變。且將其他因素置于低水平(表5)時(shí),x1、x2、x3、x6的巴豆醛邊際釋放量均為0,說明這4種工藝參數(shù)對(duì)巴豆醛釋放量無影響。將其他因素置于中水平(表6)、高水平(表7)時(shí),x1、x2、x3的巴豆醛邊際釋放量為正值,說明巴豆醛釋放量會(huì)隨其增加而增加;x6的巴豆醛邊際釋放量為負(fù)值,說明巴豆醛釋放量會(huì)隨其增加而降低。這都與單因素效應(yīng)分析結(jié)果一致。

    熱風(fēng)風(fēng)門開度(x4)的巴豆醛邊際釋放量為正值,且不因其他因素水平的改變而變化,僅隨其自身因素水平的降低而減小,說明該工藝參數(shù)的增加會(huì)促進(jìn)巴豆醛的釋放,且影響程度也會(huì)相應(yīng)增大。筒壁溫度(x5)的邊際效應(yīng)比較復(fù)雜,當(dāng)其他因素置于低水平時(shí),x5的巴豆醛邊際釋放量為正值,并隨其自身水平的提高而變大,說明巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而增加,影響速率也逐漸增大;但當(dāng)其他因素分別置于中水平和高水平時(shí),其巴豆醛邊際釋放量為負(fù)值,且絕對(duì)值隨其自身水平的提高而減小,說明此情況下巴豆醛釋放量會(huì)隨該工藝參數(shù)的增加而降低,受影響程度也逐漸減小。

    綜合以上結(jié)果來看,除了熱風(fēng)風(fēng)門開度(x4),其余因素的巴豆醛邊際釋放量絕對(duì)值均隨其他因素水平的升高而增大,說明其他因素的水平越高,巴豆醛釋放量受這5種工藝參數(shù)變化的影響也越大。

    2.6 雙因素交互作用分析

    采用等高線分析法[12]考察建立的回歸模型中任意兩種工藝參數(shù)對(duì)巴豆醛釋放量的交互作用影響結(jié)果。選定模型中的2個(gè)因素,將其他因素按照表4分別置于高、中、低水平,共得到等高線方程42個(gè),結(jié)果只有切絲寬度(x1)與HT工作蒸氣壓力(x2)、切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線方程中含有交互項(xiàng)x1x2、x1x3、x5x6,如表8所示。表明對(duì)于本試驗(yàn)建立的模型,只有x1與x2、x1與x3、x5與x6之間存在交互作用。根據(jù)其等高線方程作圖發(fā)現(xiàn),這3對(duì)因素的等高線圖基本變化趨勢(shì)相同,即隨著2種工藝參數(shù)的減小,巴豆醛釋放量逐漸減小。以工藝參數(shù)切絲寬度(x1)與HT工作蒸汽壓力(x2)的等高線圖為例,當(dāng)其他因素分別置于低、中、高水平時(shí),其等高線圖如圖1所示。切絲寬度(x1)與熱風(fēng)溫度(x3)、筒壁溫度(x5)與筒體轉(zhuǎn)速(x6)的等高線圖與圖1a一致,等高線疏密分布均勻,說明2種因素對(duì)巴豆醛釋放量的影響程度沒有明顯差別。從圖1b、1c可以看出,當(dāng)其他工藝參數(shù)置于中水平和高水平時(shí),切絲寬度(x1)與HT工作蒸氣壓力(x2)的交互作用表現(xiàn)為,沿切絲寬度(x1)軸向等高線變化密集,而HT工作蒸氣壓力(x2)軸向等高線變化相對(duì)稀疏,說明切絲寬度(x1)對(duì)巴豆醛釋放量的影響比HT工作蒸氣壓力(x2)大。

    3 結(jié)論

    通過均勻設(shè)計(jì)安排試驗(yàn),采用二次多項(xiàng)式逐步回歸方法建立了包括切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度、筒壁溫度、筒體轉(zhuǎn)速6個(gè)工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的模型。模型的因素主效應(yīng)分析結(jié)果表明,各工藝參數(shù)對(duì)巴豆醛釋放量影響的貢獻(xiàn)率由大到小依次為筒壁溫度、熱風(fēng)風(fēng)門開度、切絲寬度、熱風(fēng)溫度、筒體轉(zhuǎn)速、HT工作蒸氣壓力。單因素效應(yīng)分析發(fā)現(xiàn),切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度與巴豆醛釋放量呈線性正相關(guān),筒體轉(zhuǎn)速則與巴豆醛釋放量呈線性負(fù)相關(guān),而熱風(fēng)風(fēng)門開度和筒壁溫度與巴豆醛釋放量呈較復(fù)雜的拋物線相關(guān)。進(jìn)一步的邊際效應(yīng)分析證實(shí),在0~1水平下,巴豆醛釋放量隨切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、熱風(fēng)風(fēng)門開度的增大而提高,隨筒體轉(zhuǎn)速和筒壁溫度的升高而降低;且切絲寬度、HT工作蒸氣壓力、熱風(fēng)溫度、筒壁溫度和筒體轉(zhuǎn)速對(duì)巴豆醛釋放量影響的變化速率隨其他工藝參數(shù)的增大而變大,熱風(fēng)風(fēng)門開度的影響不變。雙因素交互作用的分析結(jié)果表明,只有切絲寬度與HT工作蒸氣壓力、切絲寬度與熱風(fēng)溫度、筒壁溫度與筒體轉(zhuǎn)速之間存在交互作用。本研究建立的模型對(duì)于從工藝參數(shù)的角度控制卷煙主流煙氣中巴豆醛的釋放量具有一定的指導(dǎo)意義,在卷煙實(shí)際生產(chǎn)中還需要考慮到其他因素如卷煙材料、葉組配方等的影響。

    參考文獻(xiàn):

    [1] MANNING D L, ASKERINEC M P, JENKINS R A,et al. High performance liquid chromatographic determinations of selected gas phase carbonyls in tobacco smoke[J]. Journal of Association of Analytical Chemistry,1983,66(1): 8-12.

    [2] WYNDER E L,HOFFMAN D. Tobacco and tobacco smoke[M]. New York: Academic Press,1967.

    [3] HOFFMANN I, HOFFMANN D. The Changing Cigarette: Chemical Studies and Bioassays[M]. New York:Oxford University Press,2001.

    [4] 謝劍平,劉惠民,朱茂祥,等.卷煙煙氣危害性指數(shù)研究[J]. 煙草科技,2009(2):5-15.

    [5] 王 鵬,劉 華,曾 建,等.制絲工藝參數(shù)對(duì)主流煙氣中氫氰酸含量影響研究[J].湖北農(nóng)業(yè)科學(xué),2011,50(2):365-367.

    [6] 譚蘭蘭,施豐成,薛 芳,等.卷煙制絲工藝參數(shù)對(duì)主流煙氣中苯并[a]芘釋放量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(1):165-167.

    [7] 陳昆燕,薛 芳,曾 建,等.重點(diǎn)工序工藝參數(shù)與卷煙主流煙氣中巴豆醛釋放量的關(guān)系研究[J].江西農(nóng)業(yè)學(xué)報(bào),2009,21(12):40-42.

    [8] 秦前浩.卷煙工藝測(cè)試與分析大綱[M]. 成都:四川大學(xué)出版社,2004.

    [9] YC/T 254-2008,卷煙主流煙氣中主要羰基化合物的測(cè)定 高效液相色譜法[S].

    [10] 方開泰.均勻設(shè)計(jì)與均勻設(shè)計(jì)表[M].北京:科學(xué)出版社,1994.

    [11] 唐啟義,馮明光.實(shí)用統(tǒng)計(jì)分析及其DPS數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [12] 何 為,薛衛(wèi)東,唐 斌.優(yōu)化試驗(yàn)設(shè)計(jì)方法及數(shù)據(jù)分析[M].北京:化學(xué)工業(yè)出版社,2012.

    猜你喜歡
    卷煙
    卷煙智能自助迷你倉配送服務(wù)模式探索與應(yīng)用
    基于虛擬視窗的異型卷煙分揀線研制
    卷煙物流節(jié)能新途徑
    卷煙包裝痕跡分析
    卷煙工業(yè)輔料倉儲(chǔ)的建立
    河南科技(2014年3期)2014-02-27 14:05:46
    頂空-氣相色譜-質(zhì)譜技術(shù)應(yīng)用于卷煙真?zhèn)舞b別
    晾曬煙葉模塊在烤煙型卷煙中的應(yīng)用
    降低卷煙主流煙氣NNK釋放量的技術(shù)研究進(jìn)展
    我國(guó)卷煙需求預(yù)測(cè)研究述評(píng)
    卷煙引燃傾向測(cè)試方法的研究與改進(jìn)
    亚洲精品日韩在线中文字幕| 91久久精品国产一区二区成人| 亚洲成人久久爱视频| 九九爱精品视频在线观看| kizo精华| 国产精品无大码| 男女国产视频网站| 亚洲欧美日韩另类电影网站 | 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| 边亲边吃奶的免费视频| 久久久国产一区二区| 精品久久久久久久末码| 国产综合懂色| 日本-黄色视频高清免费观看| 少妇人妻久久综合中文| 色播亚洲综合网| 久热久热在线精品观看| 伊人久久精品亚洲午夜| 99热国产这里只有精品6| 18禁在线播放成人免费| 嘟嘟电影网在线观看| 日韩成人伦理影院| 亚洲精品久久久久久婷婷小说| 麻豆国产97在线/欧美| 色婷婷久久久亚洲欧美| 天堂网av新在线| 亚洲国产高清在线一区二区三| 99热网站在线观看| 嫩草影院入口| 老女人水多毛片| av福利片在线观看| 成人国产麻豆网| 亚洲,一卡二卡三卡| 欧美3d第一页| 久久久久国产网址| 看非洲黑人一级黄片| 久久久久久久精品精品| 内地一区二区视频在线| 亚洲经典国产精华液单| 一本色道久久久久久精品综合| 欧美xxxx黑人xx丫x性爽| 国产精品伦人一区二区| 亚洲av二区三区四区| 最新中文字幕久久久久| 国产免费一级a男人的天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美变态另类bdsm刘玥| 亚洲精品,欧美精品| 国产视频首页在线观看| 亚洲精品中文字幕在线视频 | 狂野欧美激情性bbbbbb| 亚洲国产精品999| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| 亚洲精品中文字幕在线视频 | 国产爽快片一区二区三区| 久久鲁丝午夜福利片| 免费观看a级毛片全部| 免费观看性生交大片5| 一级av片app| 国产探花极品一区二区| 一区二区av电影网| 热99国产精品久久久久久7| 久久久久久久亚洲中文字幕| 亚洲精品日韩在线中文字幕| 久久6这里有精品| 熟妇人妻不卡中文字幕| 你懂的网址亚洲精品在线观看| 国产成人免费观看mmmm| 简卡轻食公司| 特级一级黄色大片| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 直男gayav资源| 色网站视频免费| 亚洲精品乱久久久久久| 亚洲av免费高清在线观看| 男女那种视频在线观看| 男女边摸边吃奶| 在线观看国产h片| 热re99久久精品国产66热6| 在线精品无人区一区二区三 | 色播亚洲综合网| 99久久人妻综合| 色网站视频免费| 国产精品国产三级国产av玫瑰| 午夜福利视频1000在线观看| 国产成人午夜福利电影在线观看| 26uuu在线亚洲综合色| 少妇猛男粗大的猛烈进出视频 | 色哟哟·www| 在线免费观看不下载黄p国产| 午夜免费男女啪啪视频观看| 亚洲精品成人久久久久久| 亚洲aⅴ乱码一区二区在线播放| 免费黄网站久久成人精品| 老司机影院毛片| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 国产毛片在线视频| 美女国产视频在线观看| 国产精品一区www在线观看| 午夜视频国产福利| 精华霜和精华液先用哪个| 男男h啪啪无遮挡| 精品亚洲乱码少妇综合久久| 亚洲av欧美aⅴ国产| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 亚洲欧美日韩东京热| 欧美bdsm另类| 欧美日韩视频高清一区二区三区二| 一级毛片我不卡| 香蕉精品网在线| 国产亚洲一区二区精品| 国产一区二区三区av在线| 日韩精品有码人妻一区| 国产高潮美女av| 久久女婷五月综合色啪小说 | 最近2019中文字幕mv第一页| 亚洲人成网站在线观看播放| 亚洲av中文av极速乱| 国产成人福利小说| a级毛色黄片| 国产毛片在线视频| 国产精品国产三级专区第一集| 一边亲一边摸免费视频| 青春草国产在线视频| 18+在线观看网站| 国产亚洲5aaaaa淫片| 久久人人爽av亚洲精品天堂 | 日日摸夜夜添夜夜添av毛片| 黄片无遮挡物在线观看| 国产成人免费无遮挡视频| 久久精品国产a三级三级三级| 色网站视频免费| 国产精品国产三级国产av玫瑰| 大香蕉97超碰在线| 亚洲精品一二三| 亚洲成人精品中文字幕电影| 男男h啪啪无遮挡| 国产在线一区二区三区精| h日本视频在线播放| 麻豆久久精品国产亚洲av| 国产欧美另类精品又又久久亚洲欧美| 天天躁日日操中文字幕| 99热这里只有精品一区| av在线蜜桃| 一级毛片电影观看| 欧美极品一区二区三区四区| 麻豆久久精品国产亚洲av| 久久ye,这里只有精品| 色视频www国产| 久久久久国产网址| 精品午夜福利在线看| 日本一二三区视频观看| 欧美老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 日韩免费高清中文字幕av| 成人欧美大片| 久久精品国产自在天天线| 亚洲高清免费不卡视频| 99久久精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 亚洲av一区综合| 国产爽快片一区二区三区| 国产高清不卡午夜福利| 亚洲精品国产成人久久av| 欧美丝袜亚洲另类| 人妻制服诱惑在线中文字幕| 黑人高潮一二区| 午夜爱爱视频在线播放| 免费大片黄手机在线观看| 久久久成人免费电影| 成年人午夜在线观看视频| kizo精华| 亚洲av男天堂| 99热这里只有是精品在线观看| 成人高潮视频无遮挡免费网站| 高清毛片免费看| 亚洲,欧美,日韩| 亚洲欧洲国产日韩| 亚洲欧美日韩卡通动漫| av国产精品久久久久影院| 免费电影在线观看免费观看| 亚洲精品,欧美精品| 国产免费一级a男人的天堂| 国产有黄有色有爽视频| 亚洲精品aⅴ在线观看| 蜜臀久久99精品久久宅男| 亚洲自拍偷在线| 久久久久久久久大av| a级毛片免费高清观看在线播放| 久久久久久久精品精品| 最近手机中文字幕大全| 国产黄色视频一区二区在线观看| 婷婷色av中文字幕| 日韩欧美一区视频在线观看 | 午夜福利网站1000一区二区三区| 亚洲色图综合在线观看| 成人欧美大片| 美女国产视频在线观看| 日本黄色片子视频| 制服丝袜香蕉在线| 女人十人毛片免费观看3o分钟| 最近最新中文字幕大全电影3| 日韩一本色道免费dvd| 免费观看性生交大片5| 老女人水多毛片| 亚洲真实伦在线观看| 色播亚洲综合网| 夜夜爽夜夜爽视频| 26uuu在线亚洲综合色| 毛片女人毛片| 国产成人精品福利久久| 久久ye,这里只有精品| 亚洲av在线观看美女高潮| 亚洲精品aⅴ在线观看| 国产av国产精品国产| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件| 国产片特级美女逼逼视频| 麻豆成人av视频| 爱豆传媒免费全集在线观看| 99热这里只有精品一区| 街头女战士在线观看网站| 能在线免费看毛片的网站| 成人午夜精彩视频在线观看| 高清日韩中文字幕在线| 国产欧美日韩一区二区三区在线 | 一级av片app| 插逼视频在线观看| 亚洲欧美日韩卡通动漫| 久久亚洲国产成人精品v| 草草在线视频免费看| 国产美女午夜福利| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜爱| 男人添女人高潮全过程视频| 亚洲国产av新网站| 免费看av在线观看网站| 亚洲精品自拍成人| 国产亚洲av嫩草精品影院| 国产探花极品一区二区| 国产高清国产精品国产三级 | 免费观看av网站的网址| 国产av国产精品国产| 99久久精品一区二区三区| 草草在线视频免费看| 少妇熟女欧美另类| 欧美xxⅹ黑人| 好男人视频免费观看在线| 另类亚洲欧美激情| 中文精品一卡2卡3卡4更新| 91精品国产九色| 交换朋友夫妻互换小说| 精品久久久精品久久久| www.av在线官网国产| 国产男女内射视频| 国产日韩欧美在线精品| 99热网站在线观看| 成年av动漫网址| 男女啪啪激烈高潮av片| 中文乱码字字幕精品一区二区三区| 成人高潮视频无遮挡免费网站| 国产一区二区亚洲精品在线观看| 国产高清三级在线| 一区二区av电影网| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 天天躁日日操中文字幕| 久久鲁丝午夜福利片| 亚洲内射少妇av| 亚洲久久久久久中文字幕| 国产亚洲5aaaaa淫片| 日韩亚洲欧美综合| 欧美区成人在线视频| 一级毛片久久久久久久久女| 80岁老熟妇乱子伦牲交| 亚洲av福利一区| 蜜桃亚洲精品一区二区三区| 亚洲一区二区三区欧美精品 | 亚洲精品久久久久久婷婷小说| 春色校园在线视频观看| 国产黄片美女视频| 久久99热这里只频精品6学生| 国产淫语在线视频| 99热网站在线观看| av线在线观看网站| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 男人舔奶头视频| 日本一二三区视频观看| 全区人妻精品视频| 一区二区三区免费毛片| 伦理电影大哥的女人| 中国国产av一级| 中文字幕久久专区| 国产精品人妻久久久影院| 蜜桃久久精品国产亚洲av| 一区二区三区精品91| 国产片特级美女逼逼视频| 色网站视频免费| 亚洲国产精品成人综合色| 久久亚洲国产成人精品v| 国产精品三级大全| 99久久精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产高清三级在线| 国产成年人精品一区二区| 日本三级黄在线观看| 久久国产乱子免费精品| 成人鲁丝片一二三区免费| 国内少妇人妻偷人精品xxx网站| 久久国内精品自在自线图片| 欧美xxⅹ黑人| 日韩亚洲欧美综合| 永久免费av网站大全| 亚洲va在线va天堂va国产| 少妇高潮的动态图| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放| 69人妻影院| 亚洲欧美成人精品一区二区| 亚洲精品国产色婷婷电影| 日韩视频在线欧美| 国产男人的电影天堂91| 国产淫语在线视频| 国产一区二区亚洲精品在线观看| 国产午夜精品一二区理论片| 成人特级av手机在线观看| 精品一区二区三卡| 啦啦啦在线观看免费高清www| 精品人妻一区二区三区麻豆| 亚洲伊人久久精品综合| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| 中文字幕制服av| 欧美日韩视频精品一区| 国产精品久久久久久精品古装| 婷婷色麻豆天堂久久| 亚洲国产精品成人久久小说| 日韩一本色道免费dvd| 久久久久精品性色| 日韩av在线免费看完整版不卡| 好男人在线观看高清免费视频| 九九在线视频观看精品| 午夜老司机福利剧场| 久久久久久久久久成人| 日韩制服骚丝袜av| 亚洲av一区综合| 亚洲精品乱码久久久v下载方式| 久久女婷五月综合色啪小说 | 毛片女人毛片| 国产精品不卡视频一区二区| 最新中文字幕久久久久| 亚洲成人久久爱视频| 国产一区二区三区综合在线观看 | 亚洲一区二区三区欧美精品 | 免费大片黄手机在线观看| 国产欧美另类精品又又久久亚洲欧美| 人妻系列 视频| 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 91精品一卡2卡3卡4卡| 亚洲精品国产色婷婷电影| 国产精品国产av在线观看| 欧美另类一区| 亚洲欧美成人精品一区二区| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区成人| 大香蕉久久网| 性插视频无遮挡在线免费观看| 伦理电影大哥的女人| 亚洲精品国产色婷婷电影| 国产亚洲av片在线观看秒播厂| 黄色日韩在线| 看免费成人av毛片| 久久国产乱子免费精品| 一本一本综合久久| 另类亚洲欧美激情| 久久精品国产自在天天线| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看| 国产成人精品婷婷| 麻豆精品久久久久久蜜桃| 久久久久久久久久久免费av| 一级毛片电影观看| 欧美成人精品欧美一级黄| 欧美日韩视频精品一区| 亚洲精品乱码久久久v下载方式| 久久99热这里只有精品18| 久久精品国产亚洲av天美| 草草在线视频免费看| av在线app专区| 七月丁香在线播放| 在线观看人妻少妇| 有码 亚洲区| 免费看不卡的av| 噜噜噜噜噜久久久久久91| 亚洲精品一二三| 观看美女的网站| 亚洲最大成人中文| 性插视频无遮挡在线免费观看| 成人亚洲欧美一区二区av| 午夜福利视频1000在线观看| 午夜福利网站1000一区二区三区| 一级黄片播放器| 国产极品天堂在线| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 女的被弄到高潮叫床怎么办| 日本wwww免费看| 国产乱人偷精品视频| 七月丁香在线播放| 国产欧美日韩精品一区二区| 我要看日韩黄色一级片| 国产精品一区二区三区四区免费观看| 亚洲综合精品二区| 边亲边吃奶的免费视频| 麻豆乱淫一区二区| av免费观看日本| 国产久久久一区二区三区| 中国三级夫妇交换| 一级毛片aaaaaa免费看小| 欧美少妇被猛烈插入视频| kizo精华| 在线观看免费高清a一片| 免费av不卡在线播放| 精品久久久精品久久久| 22中文网久久字幕| 久久精品久久久久久久性| 免费观看a级毛片全部| 亚洲图色成人| 日韩大片免费观看网站| 一级毛片aaaaaa免费看小| 久久影院123| 亚洲在久久综合| 大片免费播放器 马上看| 精品熟女少妇av免费看| 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 香蕉精品网在线| 美女国产视频在线观看| 99久久精品一区二区三区| 欧美性感艳星| 亚洲色图av天堂| 免费观看的影片在线观看| 亚洲精品aⅴ在线观看| 麻豆乱淫一区二区| 亚洲欧美成人综合另类久久久| 中文在线观看免费www的网站| 在线观看三级黄色| 日韩在线高清观看一区二区三区| 永久免费av网站大全| 国产 精品1| 免费看a级黄色片| 免费黄网站久久成人精品| 国产一区亚洲一区在线观看| 国产精品国产av在线观看| 亚洲色图综合在线观看| 成年女人看的毛片在线观看| 日韩一区二区三区影片| 男女下面进入的视频免费午夜| av黄色大香蕉| 国产成人精品婷婷| 特级一级黄色大片| 亚洲综合精品二区| 久久久午夜欧美精品| 91精品国产九色| 免费电影在线观看免费观看| 国产精品人妻久久久久久| 丝袜美腿在线中文| 天堂中文最新版在线下载 | 精品一区二区三卡| 综合色丁香网| 大陆偷拍与自拍| 男人爽女人下面视频在线观看| 国产精品偷伦视频观看了| 国产久久久一区二区三区| 大片电影免费在线观看免费| 中国美白少妇内射xxxbb| 丝袜美腿在线中文| 国产免费一级a男人的天堂| 美女高潮的动态| 少妇人妻 视频| 丝袜喷水一区| 亚洲最大成人av| 国产毛片a区久久久久| 日韩免费高清中文字幕av| 午夜日本视频在线| 亚洲图色成人| 成人特级av手机在线观看| 黄色一级大片看看| 男人和女人高潮做爰伦理| 亚洲精品久久久久久婷婷小说| 精品99又大又爽又粗少妇毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 97人妻精品一区二区三区麻豆| 精品亚洲乱码少妇综合久久| 亚洲精品国产成人久久av| 女人久久www免费人成看片| 国产亚洲精品久久久com| 亚洲精品视频女| 午夜精品国产一区二区电影 | 91精品国产九色| 97超视频在线观看视频| 免费av观看视频| 久久韩国三级中文字幕| 纵有疾风起免费观看全集完整版| 特大巨黑吊av在线直播| 国产av不卡久久| 天天一区二区日本电影三级| 禁无遮挡网站| 国产黄频视频在线观看| av黄色大香蕉| 一级二级三级毛片免费看| 黄色日韩在线| 少妇熟女欧美另类| 亚洲av.av天堂| 久热这里只有精品99| 久久精品久久精品一区二区三区| 99re6热这里在线精品视频| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| 毛片一级片免费看久久久久| 女人久久www免费人成看片| av在线播放精品| 欧美+日韩+精品| 免费观看无遮挡的男女| 三级经典国产精品| 色播亚洲综合网| 亚洲成色77777| 18禁裸乳无遮挡免费网站照片| 欧美日韩一区二区视频在线观看视频在线 | 97在线视频观看| 小蜜桃在线观看免费完整版高清| 青春草亚洲视频在线观看| 亚洲在久久综合| 亚洲aⅴ乱码一区二区在线播放| av福利片在线观看| 日本一二三区视频观看| 人妻夜夜爽99麻豆av| 简卡轻食公司| 亚洲精品乱码久久久久久按摩| 少妇熟女欧美另类| 少妇人妻一区二区三区视频| 特大巨黑吊av在线直播| 一级a做视频免费观看| 国产成人精品福利久久| 国产精品一区www在线观看| 日韩中字成人| 成人午夜精彩视频在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲精品自拍成人| 亚洲天堂国产精品一区在线| 国产精品女同一区二区软件| 日本免费在线观看一区| 成人鲁丝片一二三区免费| 777米奇影视久久| 免费看光身美女| 国产人妻一区二区三区在| 99热国产这里只有精品6| 欧美潮喷喷水| 王馨瑶露胸无遮挡在线观看| 国产黄色免费在线视频| 免费黄频网站在线观看国产| 免费大片18禁| 国产免费一级a男人的天堂| 亚洲久久久久久中文字幕| 国产亚洲最大av| 国产有黄有色有爽视频| 91精品国产九色| 精品久久国产蜜桃| 亚洲成人av在线免费| 国产精品国产三级国产专区5o| 爱豆传媒免费全集在线观看| 国产欧美亚洲国产| 久久久久国产精品人妻一区二区| 高清日韩中文字幕在线| 不卡视频在线观看欧美| 亚洲国产精品成人久久小说| 亚洲怡红院男人天堂| 亚洲无线观看免费| 99久久精品国产国产毛片| 亚洲av欧美aⅴ国产| 韩国av在线不卡| 免费av毛片视频| 99re6热这里在线精品视频| 久久国内精品自在自线图片| 亚洲av中文av极速乱| 午夜精品国产一区二区电影 | 国产色爽女视频免费观看| 国产高清有码在线观看视频| 亚洲av欧美aⅴ国产| 又爽又黄a免费视频| 久久久久久久精品精品| 日韩一本色道免费dvd| 国产 一区 欧美 日韩| 国产av国产精品国产| 神马国产精品三级电影在线观看| 中文在线观看免费www的网站| 91在线精品国自产拍蜜月| 亚洲一区二区三区欧美精品 | 欧美性感艳星| 边亲边吃奶的免费视频| 尤物成人国产欧美一区二区三区| av在线天堂中文字幕| 寂寞人妻少妇视频99o| av线在线观看网站| 色视频在线一区二区三区| 中文字幕av成人在线电影| 又黄又爽又刺激的免费视频.|