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

    微型流化床內(nèi)松木屑和煤泥等溫混合熱解特性

    2017-10-16 08:13:32李天濤郭飛強(qiáng)王巖郭成龍董玉平
    化工學(xué)報(bào) 2017年10期
    關(guān)鍵詞:木屑煤泥流化床

    李天濤,郭飛強(qiáng),王巖,郭成龍,董玉平

    (1中國(guó)礦業(yè)大學(xué)電氣與動(dòng)力工程學(xué)院,江蘇 徐州 221116;2山東大學(xué)機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    微型流化床內(nèi)松木屑和煤泥等溫混合熱解特性

    李天濤1,郭飛強(qiáng)1,王巖1,郭成龍1,董玉平2

    (1中國(guó)礦業(yè)大學(xué)電氣與動(dòng)力工程學(xué)院,江蘇 徐州 221116;2山東大學(xué)機(jī)械工程學(xué)院,山東 濟(jì)南 250061)

    采用微型流化床反應(yīng)器對(duì)松木屑和煤泥的等溫混合熱解氣體釋放行為進(jìn)行實(shí)驗(yàn)研究。探討不同溫度和摻混比例對(duì)CH4、CO、CO2與H2釋放特性的影響,并通過(guò)模型配合法求解其動(dòng)力學(xué)參數(shù),研究松木屑和煤泥混合熱解過(guò)程的相互作用。通過(guò)FT-IR檢測(cè)發(fā)現(xiàn)煤泥的主要成分為含有C—O和CO鍵的芳香化合物,松木屑則以帶—OH鍵的長(zhǎng)鏈脂肪烴為主。在等溫穩(wěn)定反應(yīng)階段,松木屑熱解氣體生成速率高于煤泥,隨著生物質(zhì)摻入比例的不斷提高,混合原料氣體生成反應(yīng)速率亦呈現(xiàn)不同程度的增加。利用模型積分法求解了松木屑、煤泥及其混合物熱解氣體生成動(dòng)力學(xué)參數(shù),并通過(guò)實(shí)驗(yàn)值和計(jì)算值對(duì)比篩選出了最概然機(jī)理函數(shù)。通過(guò)活化能對(duì)比發(fā)現(xiàn),混合熱解對(duì)4種氣體組分生成具有不同的影響作用,其中CO實(shí)驗(yàn)活化能明顯低于計(jì)算值,表現(xiàn)為二者協(xié)同作用利于CO的生成釋放;對(duì) H2而言,在 75%混合比例條件下,混合反應(yīng)導(dǎo)致其生成活化能呈現(xiàn)協(xié)同負(fù)效應(yīng),使得活化能實(shí)驗(yàn)值明顯高于計(jì)算值;相較而言,CH4在混合熱解過(guò)程影響相對(duì)較弱,并呈現(xiàn)小幅度的協(xié)同負(fù)效應(yīng),而CO2的生成特性則受混合比例的影響較為明顯。

    微型流化床;生物質(zhì);煤泥;熱解;等溫反應(yīng);反應(yīng)動(dòng)力學(xué);相互作用

    Abstract:The releasing characteristics of gas species during isothermal co-pyrolysis of pine sawdust (PS) and coal slime (CS) were studied in a micro fluidized bed reactor.Four main pyrolysis gaseous products,CH4,CO,CO2and H2,were investigated under different temperature and blending ratios.The kinetic parameters were calculated employing the universal integral method to examine the interaction of PS and CS during co-pyrolysis.It was observed by FT-IR that CS mainly contained aromatic compounds carrying C—O and CO,while the PS was mainly consist of long-chain aliphatic hydrocarbon with —OH.The gas forming reaction rate of PS was higher than that of CS,and its value increased in a certain extent with the increase of biomass blending ratio.The kinetic parameters of pyrolysis gas formation of PS,CS and their blends were obtained by the mechanism function models,and the most probable mechanism function was determined by comparing the experimental and calculated values.The effect of blending PS and CS on the forming activations of the four gas species varied obviously.The activation energies of CO were significantly lower than the calculated ones,indicating that a positive synergisticeffect existed between PS and CS.The co-pyrolysis reaction led to a negative effect on the formation of H2at a PS blending ratio of 75% as the activation energy was obviously higher than the calculated one.In comparison,biomass blending showed relatively small positive synergistic effect on the formation of CH4,while that of CO2depended upon the blending ratio significantly.

    Key words:micro fluidized bed reactor; biomass,coal slime; pyrolysis; isothermal reaction; reaction kinetics;interaction

    引 言

    煤泥是洗煤過(guò)程中產(chǎn)生的一種副產(chǎn)品,顆粒較細(xì),水分與灰分含量較高,黏度較大,運(yùn)輸十分不便,利用不當(dāng)將造成能量的浪費(fèi)并存在潛在的環(huán)境污染風(fēng)險(xiǎn)[1-5]。目前常見的煤泥燃燒利用容易產(chǎn)生大量的細(xì)灰以及氮和硫的化合物,不利于環(huán)境的保護(hù)[6]。熱解氣化作為一種重要的潔凈煤利用技術(shù),可以實(shí)現(xiàn)煤泥在相對(duì)較低的溫度條件下轉(zhuǎn)化制備清潔燃?xì)?,但煤泥中氫元素含量較低,若在熱解氣化過(guò)程摻入富氫生物質(zhì)則能一定程度上對(duì)煤泥補(bǔ)充氫元素,改善氣體產(chǎn)物分布[7]。同時(shí),生物質(zhì)儲(chǔ)量巨大,來(lái)源廣泛,但密度和熱值較低,若與煤泥混合可提高整體能量密度,實(shí)現(xiàn)與煤泥的規(guī)?;旌限D(zhuǎn)化。再者,生物質(zhì)屬于二氧化碳零排放的新型能源,與煤泥混合熱解氣化轉(zhuǎn)化實(shí)現(xiàn)優(yōu)勢(shì)互補(bǔ),是新能源與傳統(tǒng)能源協(xié)同利用的有效途徑[8]。

    目前國(guó)內(nèi)外對(duì)生物質(zhì)和各類煤樣的混合熱解研究主要采用固定床反應(yīng)器、下降管反應(yīng)器和熱重分析儀等,尤其是熱重分析儀被廣泛用于混合熱解過(guò)程研究,是當(dāng)前研究煤和生物質(zhì)混合熱解過(guò)程相互作用的主要設(shè)備。部分研究結(jié)果證實(shí)煤與生物質(zhì)混合熱解過(guò)程可以降低焦炭產(chǎn)率,并改變氣體生成特性[9-10],視二者混合熱解存在相互作用或協(xié)同效應(yīng)。然而,部分研究結(jié)果認(rèn)為煤和生物質(zhì)混合熱解過(guò)程相互作用關(guān)系微弱,并不能說(shuō)明二者之間存在協(xié)同效應(yīng)[7,11-13]。由此可以看出,煤和生物質(zhì)混合熱解反應(yīng)過(guò)程復(fù)雜,受到多種因素的限制,包括原料種類、H/C比、升溫速率及無(wú)機(jī)金屬離子成分等。因此,研究二者之間混合熱解特性不僅要從整體轉(zhuǎn)化率的角度,還要注意主要產(chǎn)物組分生成的特性,并立足于二者處于同一反應(yīng)條件。

    近期,微型流化床反應(yīng)器研究氣固反應(yīng)動(dòng)力學(xué)特性被研究開發(fā),其可實(shí)現(xiàn)原料等溫快速流化反應(yīng),降低內(nèi)擴(kuò)散影響,在線快速檢測(cè)氣體產(chǎn)物釋放特性,進(jìn)而推測(cè)氣體生成動(dòng)力學(xué)特性[14-15]?;谀P头e分法,微型流化床分析儀已成功應(yīng)用于熱解反應(yīng)等研究中。Liu等[16]采用微型流化床對(duì)稻殼進(jìn)行催化熱解實(shí)驗(yàn),通過(guò)模型積分法求解其動(dòng)力學(xué)參數(shù),進(jìn)而研究了鉀離子對(duì)稻殼熱解的影響。Guo等[17]采用微型流化床研究了焦油催化裂解生成 H2和 CO的特性。Yu等[18]在微型流化床中進(jìn)行石墨粉燃燒反應(yīng),證實(shí)微型流化床是一種簡(jiǎn)單且可靠的動(dòng)力學(xué)研究?jī)x器。

    基于此,本文采用微型流化床分析儀研究煤泥與生物質(zhì)混合熱解動(dòng)力學(xué)特性,分析主要?dú)怏w組分的生成和轉(zhuǎn)化規(guī)律,并通過(guò)模型積分法求算動(dòng)力學(xué)參數(shù),進(jìn)而從氣體生成的角度闡明煤泥和生物質(zhì)混合熱解過(guò)程相互作用特性。

    1 實(shí)驗(yàn)材料和方法

    1.1 實(shí)驗(yàn)原料

    本文實(shí)驗(yàn)原料煤泥(CS)和松木屑(PS)取自徐州本地,其工業(yè)分析及元素分析如表1所示,煤泥有較高的固定碳和灰分含量,松木則是有較高的揮發(fā)分含量。兩種原料分別粉碎、篩分至150~250 μm,并按照生物質(zhì)質(zhì)量分?jǐn)?shù)分別為0%、25%、50%、75%、100% 5種比例進(jìn)行混合,在105℃的鼓風(fēng)干燥箱中干燥24 h,取出密封保存?zhèn)溆谩?/p>

    表1 原料工業(yè)分析及元素分析Table 1 Proximate and ultimate analysis of samples

    1.2 實(shí)驗(yàn)裝置與操作

    實(shí)驗(yàn)裝置和流程如圖1所示。

    實(shí)驗(yàn)在流化床石英反應(yīng)器上進(jìn)行,石英反應(yīng)器出口連接質(zhì)譜儀對(duì)出口氣體成分進(jìn)行連續(xù)監(jiān)測(cè)記錄。反應(yīng)器內(nèi)徑為20 mm,高150 mm,由兩層布風(fēng)板分成3段,兩層布風(fēng)板之間40 mm的區(qū)域?yàn)橹饕磻?yīng)區(qū)。詳盡的實(shí)驗(yàn)裝置介紹可以參考文獻(xiàn)[14-15]。每次實(shí)驗(yàn)前,將3 g石英砂(74~125 μm)放在石英反應(yīng)器底層作為流化介質(zhì)。前期研究表明[19],載氣流速在300 ml·min?1即可實(shí)現(xiàn)良好的流化性,并可以忽略外擴(kuò)散的影響,因而本實(shí)驗(yàn)氬氣的流量設(shè)定為 300 ml·min?1。首先,在流態(tài)化狀態(tài)下將石英反應(yīng)器加熱到設(shè)定溫度,然后用壓縮氣將10 mg樣品迅速注入石英反應(yīng)器,4種主要熱解產(chǎn)氣 H2、CO、CH4、CO2的釋放強(qiáng)度通過(guò)質(zhì)譜儀進(jìn)行連續(xù)檢測(cè)。為了保證實(shí)驗(yàn)結(jié)果的可靠性,每次實(shí)驗(yàn)重復(fù)3次。

    圖1 實(shí)驗(yàn)系統(tǒng)流程圖Fig.1 Schematic diagram of experimental system

    1.3 動(dòng)力學(xué)模型

    快速熱解實(shí)驗(yàn)中,各氣體的釋放特性由在線質(zhì)譜儀實(shí)時(shí)連續(xù)測(cè)得,產(chǎn)氣轉(zhuǎn)化率計(jì)算如下

    式中,t0,t和te分別代表開始反應(yīng)、反應(yīng)進(jìn)行到時(shí)間t和結(jié)束反應(yīng)的時(shí)間,s;分別代表t時(shí)刻和反應(yīng)結(jié)束時(shí)刻的氣體產(chǎn)量,g;φi代表氣體i的體積濃度分?jǐn)?shù),%;qv代表氣體的流量,L·min?1;MR表示氣體摩爾質(zhì)量;22.4 為標(biāo)準(zhǔn)狀態(tài)下氣體摩爾常數(shù);x代表氣體的轉(zhuǎn)化率;dx/dt代表氣體轉(zhuǎn)化速率,s?1。

    其中,f(x)代表反應(yīng)的微分機(jī)理函數(shù);k(T)代表Arrhenius速率常數(shù),由Arrhenius方程所定義,在等溫過(guò)程中是一個(gè)常量。

    本文采用模型積分法對(duì)各氣體生成過(guò)程進(jìn)行動(dòng)力學(xué)求解,計(jì)算方法如下

    表2 用于氣固反應(yīng)的典型機(jī)理函數(shù)Table 2 Typical mechanism functions using in gas-solid reactions

    式中,G(x)是機(jī)理函數(shù)的積分形式;Ea代表表觀活化能,kJ·mol?1;A表示指前因子,s?1;T代表溫度,K;R代表氣體常數(shù),8.314 J·mol?1·K?1。由式(7)可以擬合出一條機(jī)理函數(shù)G(x)關(guān)于時(shí)間t的直線,直線斜率即為k(T)。再根據(jù)式(7),擬合得到速率常數(shù)k(T)與溫度T的直線,即可根據(jù)直線的斜率與截距求得表觀活化能Ea和指前因子A。將求得的動(dòng)力學(xué)參數(shù)進(jìn)行逆向運(yùn)算可以獲得轉(zhuǎn)化率計(jì)算值,進(jìn)而通過(guò)對(duì)實(shí)驗(yàn)值和計(jì)算值的誤差分析,篩選最概然機(jī)理函數(shù)。對(duì)于氣相與固相之間的化學(xué)反應(yīng),常見的非勻相反應(yīng)模型有20種,如表2所示[20],這些函數(shù)已廣泛應(yīng)用于描述熱解反應(yīng)過(guò)程[18,20-23]。

    2 實(shí)驗(yàn)結(jié)果與討論

    2.1 生物質(zhì)和煤的FT-IR分析

    采用傅里葉變換紅外光譜儀(VERTEX 80v)分析松木屑和煤泥的結(jié)構(gòu)特征,結(jié)果如圖2所示。煤泥和松木屑在結(jié)構(gòu)上既有相似之處也有明顯不同,第一個(gè)主要區(qū)別出現(xiàn)在波數(shù) 3568~3730 cm?1范圍內(nèi),煤泥有3個(gè)連續(xù)峰峰值分別在3622 cm?1、3658 cm?1和3695 cm?1,表示煤泥中帶有小分子側(cè)鏈的亞氨基N—H結(jié)構(gòu),而松木屑在這一范圍內(nèi)沒有波峰。吸收峰 3000~3568 cm?1主要是由羥基的振動(dòng)引起,松木屑在此范圍出現(xiàn)一個(gè)較強(qiáng)的吸收峰,而煤泥在此區(qū)域的峰值較弱,其相應(yīng)基團(tuán)含量明顯低于松木屑。在 2765~3000 cm?1波數(shù)范圍內(nèi),煤泥分別在2854 cm?1和2919 cm?1出現(xiàn)兩個(gè)微弱吸收峰,松木屑則只有一個(gè)強(qiáng)峰在2903 cm?1處,表明煤泥中存在少量環(huán)烷或脂肪族的亞甲基—CH2和次甲基—CH,而松木屑中有較多的亞甲基—CH2,即松木屑中有較長(zhǎng)的脂肪鏈。由此可以看出,松木屑中含氫基團(tuán)明顯多于煤泥,在熱解過(guò)程中可以作為主要的氫源,釋放含氫氣體組分。1378~1802 cm?1范圍的吸收峰主要是由芳香結(jié)構(gòu)CC和CO伸縮振動(dòng)產(chǎn)生的,煤泥有1個(gè)主峰和1個(gè)肩峰,峰值分別在1442 cm?1和1602 cm?1處,而松木屑在這一波數(shù)范圍內(nèi)波峰多且雜,表明松木屑中有著豐富的含氧基團(tuán)。在833~1378 cm?1范圍內(nèi),煤泥出現(xiàn)了較強(qiáng)的吸收峰,在 1032 cm?1處,主要為有機(jī)質(zhì)中醇和烷基醚的C—O鍵,可能會(huì)導(dǎo)致較高的CO和CO2產(chǎn)率;松木屑中波峰仍然繁雜,主峰出現(xiàn)在1037 cm?1處,對(duì)應(yīng)于醇和烷基醚的C—O鍵和羧酸中的醚鍵。由此可見,松木屑中有豐富的脂肪族結(jié)構(gòu),且脂肪鏈較長(zhǎng),羧酸基團(tuán)及醚鍵的連接方式多樣,而煤泥中更多的是芳香族結(jié)構(gòu),相較而言更為穩(wěn)定,熱解溫度更高[24]。反應(yīng)過(guò)程中若兩種基團(tuán)發(fā)生交互作用實(shí)現(xiàn)元素互補(bǔ),可以促進(jìn)二者熱解反應(yīng)。

    圖2 松木屑和煤泥的紅外光譜圖Fig.2 FT-IR spectra of pine sawdust and coal slime samples

    2.2 松木屑和煤泥單獨(dú)熱解特性

    松木屑和煤泥在注射入石英反應(yīng)器后迅速進(jìn)入熱解反應(yīng)階段。實(shí)驗(yàn)由在線質(zhì)譜儀實(shí)時(shí)連續(xù)檢測(cè)樣本在750~900℃熱解條件下的 H2、CO、CO2和CH4釋放強(qiáng)度。圖3展示了松木屑和煤泥在750℃和 900℃下的氣體釋放強(qiáng)度隨時(shí)間的變化。可以看到在微型流化床等溫條件下,松木屑和煤泥產(chǎn)氣熱解釋放基本都在反應(yīng)開始后的15~20 s內(nèi)完成且不同氣體釋放的起止時(shí)間并不相同,不同溫度條件下各氣體的釋放強(qiáng)度也有明顯差異。

    在750~900℃熱解條件下,比較煤泥和松木屑熱解氣體釋放強(qiáng)度可以發(fā)現(xiàn),松木屑的CO釋放強(qiáng)度約為煤泥的3倍。同時(shí),松木屑熱解釋放CO2、H2和CH4的釋放強(qiáng)度同樣高于煤泥,這也與常規(guī)煤與生物質(zhì)熱解特性的預(yù)期相符合。生物質(zhì)的主要組分包括纖維素,半纖維素和木質(zhì)素,相互之間以醚鍵(R—O—R)連接,醚鍵的鍵能為 380~420 kJ·mol?1[25],熱解條件下能夠釋放大量的氣體。相比之下,煤泥中固定碳和灰分含量較高,有機(jī)物主要組成為稠多環(huán)芳烴,連接鍵為芳香環(huán)鍵(CC),鍵能達(dá)到1000 kJ·mol?1[26],熱解溫度較高。另外,在相同溫度條件下生物質(zhì)含有更多的氫元素,使得含氫氣體組分(H2和CH4)產(chǎn)率更高,因而如果混合熱解,可為煤泥的熱解提供氫源,使煤熱解產(chǎn)生的大量不穩(wěn)定分子自由基得到配對(duì)電子穩(wěn)定成為低分子揮發(fā)性物質(zhì)[27]。

    在 750℃熱解條件下,松木屑和煤泥最先釋放的組分均為CO2,其次是CO、CH4和H2。隨著溫度的升高,氣體釋放順序發(fā)生變化。900℃下H2的釋放明顯提前,H2的釋放強(qiáng)度也明顯增強(qiáng)。另外CH4和CO的釋放強(qiáng)度也隨著溫度的升高而有所增強(qiáng),而CO2的釋放強(qiáng)度沒有太大變化。各氣體不同的釋放行為意味著其不同的生成路徑和反應(yīng)機(jī)理。在生物質(zhì)熱解中,CO通常由醚鍵(C—O—C)和羰基(CO)的斷裂和揮發(fā)分的二次裂解產(chǎn)生,在高溫下更容易發(fā)生[28]。H2通常由自由基縮聚反應(yīng)和脫氫反應(yīng)產(chǎn)生,CH4的形成通常伴隨—OCH3和—CH2—的斷裂[29],同樣與溫度密切相關(guān)。相較之下,CO2主要由 CO鍵和羧基(—COOH)的裂解和重整產(chǎn)生[28],鍵能相對(duì)較弱,因此CO2的產(chǎn)生受熱解溫度的影響較小。

    圖3 熱解氣體釋放特性Fig.3 Gas releasing characteristics during pyrolysis process

    圖4 不同溫度下的各氣體轉(zhuǎn)化率隨時(shí)間的變化Fig.Conversion profiles of gas species related to time and temperature

    圖4為不同溫度條件下4種氣體相對(duì)生成率隨時(shí)間的變化。對(duì)于松木屑和煤泥,隨著溫度的升高,4種氣體完成轉(zhuǎn)化的時(shí)間都相應(yīng)變短,表明產(chǎn)氣反應(yīng)與溫度高度相關(guān)。除此之外,不同氣體完成轉(zhuǎn)化的時(shí)間也有明顯差別,同樣暗示著其不同的化學(xué)反應(yīng)生成路徑?;谑?5),將圖4中x對(duì)t求微分,可以進(jìn)一步得到轉(zhuǎn)化速率dx/dt與轉(zhuǎn)化率x之間的關(guān)系,如圖5所示。從圖5可以看出在同一轉(zhuǎn)化率下,各氣體的轉(zhuǎn)化速率都隨著溫度的升高而加快。在反應(yīng)初始階段,原料快速升溫并迅速釋放氣體,使得氣體轉(zhuǎn)化速率迅速增長(zhǎng)并在轉(zhuǎn)化率0.25附近達(dá)到最大值。隨后,氣體轉(zhuǎn)化速率隨著轉(zhuǎn)化率的增加逐漸降低,表明熱解反應(yīng)已經(jīng)進(jìn)行到等溫?zé)峤怆A段。反應(yīng)進(jìn)行到后期,氣體組分的生成主要依賴于芳香核的縮合反應(yīng),導(dǎo)致反應(yīng)速率逐漸下降[30]。

    圖5 松木屑和煤泥熱解氣生成反應(yīng)速率隨轉(zhuǎn)化率的變化Fig.5 Reaction rate of gas releasingversusconversion of CS and PS

    圖6 混合樣本產(chǎn)氣釋放速率與轉(zhuǎn)化率的關(guān)系Fig.6 Reaction rate of gas releasingversusconversion of blends

    2.3 混合物的混合熱解特性

    當(dāng)在煤泥中添加松木屑后,以850℃為例,不同混合比例的混合樣品在熱解過(guò)程中各氣體轉(zhuǎn)化速率隨轉(zhuǎn)化率的變化如圖6所示。整體而言,4種氣體在穩(wěn)定的等溫?zé)峤夥磻?yīng)區(qū)(x>0.25),隨著添加松木屑添加比例的提高,反應(yīng)速率均呈現(xiàn)不同程度的升高,表明松木屑熱解反應(yīng)活性高于煤泥。但松木屑的添加對(duì)4種氣體的影響也存在明顯區(qū)別,其中H2受生物質(zhì)摻入比例的影響最為明顯,反應(yīng)速率隨松木屑比例增加呈現(xiàn)明顯增加的趨勢(shì),當(dāng)木屑摻入比例達(dá)到75%時(shí),H2的釋放反應(yīng)速率接近于純松木屑。對(duì)CO而言,松木屑的混合明顯增加了CO的轉(zhuǎn)化速率,特別是當(dāng)生物質(zhì)摻入比例為50%時(shí),其反應(yīng)速率高于75%比例的混合樣品。然而煤泥和松木屑熱解生成CH4和CO2的反應(yīng)速率較為接近,混合熱解對(duì)于兩種氣體轉(zhuǎn)化速率影響較為微弱。由此可見,混合原料熱解過(guò)程中,氣體生成速率一定程度上取決于煤泥和松木屑自身的反應(yīng)速率,另外在部分混合比例條件下反應(yīng)速率出現(xiàn)一定的變化,可能出現(xiàn)了一定的相互作用。

    2.4 模型配合法求解動(dòng)力學(xué)參數(shù)

    本文采用模型積分法對(duì)各氣體進(jìn)行動(dòng)力學(xué)求解。將各機(jī)理模型函數(shù)分別代入式(8)進(jìn)行擬合,得到線性擬合程度較好的機(jī)理模型函數(shù)為G3(三維球形擴(kuò)散模型)、G7(1級(jí)成核生長(zhǎng)模型)和G17(1.5級(jí)化學(xué)反應(yīng)模型)。在松木屑與煤泥 1:1混合比例下,機(jī)理模型函數(shù)G3、G7、G17與反應(yīng)時(shí)間t的關(guān)系如圖7所示,其擬合相關(guān)系數(shù)列于表3,4種氣體相關(guān)系數(shù)均在0.95以上,其他比例下的擬合結(jié)果與之類似,表明3種模型能夠較好地描述4種氣體在微型流化床中的等溫?zé)峤忉尫盘匦浴?種機(jī)理函數(shù)模型獲得的動(dòng)力學(xué)參數(shù)如表4和表5所示。為篩選出最概然機(jī)理函數(shù),本文將模型配合法計(jì)算得到的表觀活化能依次代入式(8)和式(7)求算轉(zhuǎn)化率的計(jì)算值,并與實(shí)驗(yàn)值進(jìn)行比較,采用式(9)計(jì)算實(shí)驗(yàn)值與計(jì)算值的平均相對(duì)誤差[31],式中定義 OF為平均相對(duì)誤差值,OF最小的機(jī)理函數(shù)即被認(rèn)定為最概然機(jī)理函數(shù)。

    圖7 煤泥與松木屑1:1混合樣本的機(jī)理模型函數(shù)Fig.7 Linear goodness of fit mechanism function model determination for 50% blending ratio of PS

    以松木屑與煤泥1:1摻混為例比較實(shí)驗(yàn)值與計(jì)算值的差異如圖8所示。計(jì)算結(jié)果顯示G3模型擬合的OF值范圍為0.03~0.07,G7的OF值集中在0.03~0.05,G17的OF值大部分在0.04~0.08之間,G7可以實(shí)現(xiàn)對(duì)于氣體轉(zhuǎn)化更好的計(jì)算預(yù)測(cè)。由表5可以看出,通過(guò)G7機(jī)理模型函數(shù)進(jìn)行計(jì)算,煤泥熱解生成H2、CO、CO2和CH4的活化能分別為75.0、50.3、37.6、50.8 kJ·mol?1;相應(yīng)的松木屑 4 種氣體生成活化能分別為 70.3、37.9、54.9、68.0 kJ·mol?1。Guo等[15]和Yu等[14]采用類似微型流化床分別研究中藥渣和酒糟熱解動(dòng)力學(xué)特性,獲得的4種氣體熱解活化能數(shù)值與本實(shí)驗(yàn)結(jié)果相當(dāng)。煤泥熱解反應(yīng)過(guò)程 CO2生成活化能最低,H2最高,說(shuō)明熱解反應(yīng)中 CO2最易生成,H2釋放最為困難;而松木屑熱解反應(yīng) CO生成活化能最低,H2最高,因此所獲得的活化能值在一定程度上也與氣體的釋放順序相一致。

    表3 煤泥松木屑1:1混合樣本的機(jī)理模型函數(shù)的線性相關(guān)系數(shù)Table 3 Linear correlation coefficient of mechanism function model determination for 50% blending ratio of PS

    表4 不同混合比例下各氣體指前因子Table 4 Pre-exponential factor of gas species related to blend ratios in different mechanism functions

    表5 不同混合比例下各氣體表觀活化能Table 5 Apparent activation energy of gas species related to blend ratios in different mechanism functions

    2.5 煤泥與生物質(zhì)相互作用特性分析

    在混合熱解反應(yīng)中,煤泥和生物質(zhì)的相互作用可以通過(guò)預(yù)測(cè)其混合樣本的活化能進(jìn)行討論。Goldfarb等[32]認(rèn)為,如果混合樣本在反應(yīng)動(dòng)力學(xué)方面沒有協(xié)同效應(yīng),則混合樣本的活化能應(yīng)為獨(dú)立樣本活化能的加權(quán)平均值?;旌蠘颖镜挠?jì)算活化能可由式(10)得到

    其中,EC代表計(jì)算活化能,kJ·mol?1;EPS和ECS分別表示松木屑和煤泥的實(shí)驗(yàn)活化能;kJ·mol?1;vPS代表松木屑的混合質(zhì)量比,%。

    圖9為采用G7模型計(jì)算得到的各氣體生成活化能的實(shí)驗(yàn)值和計(jì)算值對(duì)比。可以看出混合樣本與單獨(dú)樣本之間并不遵從線性關(guān)系,4種產(chǎn)氣中,混合樣本的H2和CH4實(shí)驗(yàn)活化能較計(jì)算活化能偏高,CO的實(shí)驗(yàn)活化能較計(jì)算活化能偏低。

    為了進(jìn)一步表征松木屑和煤泥間的相互作用,定義計(jì)算活化能與實(shí)驗(yàn)活化能的相對(duì)偏差δ如式(11)所示

    式中,Emeasured代表實(shí)驗(yàn)活化能,kJ·mol?1。不同比例下松木屑與煤泥混合樣本活化能相對(duì)偏差如圖10所示。

    對(duì)于CH4而言,在松木屑混合比例為25%、50%和75%下,相對(duì)偏差均為正值,分別為8.0%、8.6%、7.4%,表明實(shí)驗(yàn)活化能偏高,兩者的混合熱解對(duì)CH4的產(chǎn)生具有一定的不利作用。對(duì)于CO而言,當(dāng)松木屑混合比例為 25%、50%、75%時(shí),活化能實(shí)驗(yàn)值和計(jì)算值呈現(xiàn)明顯偏差,分別為?26.0%、?30.7%、?29.9%,說(shuō)明混合原料熱解過(guò)程CO釋放所需能量壘勢(shì)降低,對(duì)CO的生成有明顯的促進(jìn)作用。在CO2生成過(guò)程,在松木屑與煤泥比例為1:3時(shí),相對(duì)偏差幾乎為 0,兩者的混合熱解沒有協(xié)同效應(yīng)發(fā)生;當(dāng)松木屑與煤泥比例為1:1時(shí),相關(guān)系數(shù)為7.6%,混合熱解不利于CO2的生成;當(dāng)松木屑摻混比例提高到75%時(shí),相關(guān)系數(shù)降到?10.6%,反而促進(jìn)了 CO2的產(chǎn)生。對(duì) H2而言,當(dāng)松木屑混合比例為 50%和 75%時(shí),相關(guān)系數(shù)分別為 10.1%、20.5%,說(shuō)明在較高的松木屑混合比例條件下,H2生成所需能量增加,可能由于在生物質(zhì)比例過(guò)高時(shí),低密度的生物質(zhì)顆粒在熱解中軟化、變形,包裹黏附在煤泥顆粒的表面,堵塞煤泥的孔隙,抑止煤泥中揮發(fā)分的逸出和擴(kuò)散[33]。通過(guò)活化能的變化可以發(fā)現(xiàn),混合熱解過(guò)程對(duì)于不同氣體的活化能影響不同,整體而言混合熱解對(duì)于降低CO生成活化能影響較為明顯,表現(xiàn)為CO生成過(guò)程發(fā)生協(xié)同效應(yīng),降低了其生成的表觀活化能;在高松木屑摻入比例條件下H2的活化能高于理論計(jì)算值,表現(xiàn)出一定的負(fù)協(xié)同作用。

    圖8 煤泥與松木屑1:1混合樣本的機(jī)理模型函數(shù)值與實(shí)驗(yàn)值比較Fig.8 Comparison between mechanism function model value and experimental value of 50% blending ratio of PS

    圖9 G7模型下不同生物質(zhì)比例的各氣體表觀活化能Fig.9 Apparent activation energy of gas species in different biomass blending ratio in G7

    3 結(jié) 論

    (1)煤泥和松木屑在官能團(tuán)結(jié)構(gòu)存在明顯差異,主要體現(xiàn)在含氫官能團(tuán)和含氧官能團(tuán),在熱解過(guò)程中有利于官能團(tuán)相互作用和元素互補(bǔ)。

    (2)松木屑和煤泥等溫?zé)峤膺^(guò)程,松木屑?xì)怏w釋放反應(yīng)速率高于煤泥,混合比例的增加一定程度可以提高氣體生成速率。

    (3)松木屑和煤泥的混合熱解過(guò)程相互作用對(duì)于不同氣體影響存在差異,其中對(duì)于CO生成影響最為明顯,混合反應(yīng)降低了CO生成的活化能;在高生物質(zhì)摻入比例條件下,二者相互作用導(dǎo)致 H2活化能升高,表現(xiàn)出協(xié)同負(fù)效應(yīng);混合反應(yīng)對(duì)于CH4和CO2生成影響相對(duì)較弱,表現(xiàn)為混合反應(yīng)會(huì)導(dǎo)致CH4活化能小幅度增加,而CO2活化能變化受原料混合比例影響明顯。

    [1]SONG Z L,JING C M,YAO L S,et al.Coal slime hot air/microwave combined drying characteristics and energy analysis[J].Fuel Processing Technology,2017,156: 491-499.

    [2]ZHOU K,LIN Q Z,HU H W,et al.The ignition characteristics and combustion processes of the single coal slime particle under different hot-coflow conditions in N2/O2atmosphere[J].Energy,2017,136:173-184.

    [3]張坤.火電廠煤泥型煤制備及工程應(yīng)用方案初步設(shè)計(jì)[D].北京:華北電力大學(xué),2013.ZHANG K.Preliminary design on preparation and engineering applications of slime briquette in thermal power plant[D].Beijing:North China Electric Power University,2013.

    [4]劉明強(qiáng),劉建忠,王睿坤,等.煤泥燃燒動(dòng)力學(xué)機(jī)理研究的新方法探討[J].煤炭學(xué)報(bào),2012,37(a02): 444-448.LIU M Q,LIU J Z,WANG R K,et al.Discussion of a new approach to study kinetic mechanism of coal slurry combustion[J].Journal of China Coal Society,2012,37(a02): 444-448.

    [5]徐曉虹,李法香,吳建鋒,等.利用煤泥制備蓄熱陶瓷及抗熱震性能[J].武漢理工大學(xué)學(xué)報(bào),2016,36(11): 6-10.XU X H,LI F X,WU J F,et al.Preparation and thermal shock resistance of thermal storage ceramic from coal slime[J].Journal of Wuhan University of Technology,2016,36(11): 6-10.

    [6]陳恒寶,周敏,賀國(guó)章.煤泥的二氧化碳?xì)饣瘎?dòng)力學(xué)[J].煤炭轉(zhuǎn)化,2012,35(3): 19-22.CHEN H B,ZHOU M,HE G Z.Carbon dioxide gasification kinetics of coal slime[J].Coal Conversion,2012,35(3): 19-22.

    [7]WEILAND N T,MEANS N C,MORREALE B D.Product distributions from isothermal co-pyrolysis of coal and biomass[J].Fuel,2011,94(1): 563-570.

    [8]WU Z Q,WANG S Z,ZHAO J,et al.Product distribution during co-pyrolysis of bituminous coal and lignocellulosic biomass major components in a drop-tube furnace[J].Energy & Fuels,2015,29(7):4168-4180.

    [9]ZHANG L,XU S P,ZHAO W,et al.Co-pyrolysis of biomass and coal in a free fall reactor[J].Fuel,2007,86(3): 353-359.

    [10]DU Y Y,JIANG X G,LV G,et al.Thermal behavior and kinetics of bio-ferment residue/coal blends during co-pyrolysis[J].Energy Conversion & Management,2014,88: 459-463.

    [11]GIL M V,CASAL D,PEVIDA C,et al.Thermal behaviour and kinetics of coal/biomass blends during co-combustion[J].Bioresource Technology,2010,101(14): 5601-5608.

    [12]ZHOU C C,LIU G J,FANG T,et al.Investigation on thermal and trace element characteristics during co-combustion biomass with coal gangue.[J].Bioresource Technology,2015,175(175C):454-462.

    [13]LU K M,LEE W J,CHEN W H,et al.Thermogravimetric analysis and kinetics of co-pyrolysis of raw/torrefied wood and coal blends[J].Applied Energy,2013,105(2): 57-65.

    [14]YU J,YAO C B,ZENG X,et al.Biomass pyrolysis in a micro-fluidized bed reactor: characterization and kinetics[J].Chemical Engineering Journal,2011,168(2): 839-847.

    [15]GUO F Q,DONG Y P,LV Z C,et al.Pyrolysis kinetics of biomass(herb residue) under isothermal condition in a micro fluidized bed[J].Energy Conversion & Management,2015,93: 367-376.

    [16]LIU Y,WANG Y,GUO F Q,et al.Characterization of the gas releasing behaviors of catalytic pyrolysis of rice husk using potassium over a micro-fluidized bed reactor[J].Energy Conversion and Management,2017,136: 395-403.

    [17]GUO F Q,DONG Y P,FAN P F,et al.Catalytic decomposition of biomass tar compound by calcined coal gangue: a kinetic study[J].International Journal of Hydrogen Energy,2016,41(31):13380-13389.

    [18]YU J,ZENG X,ZHANG J,et al.Isothermal differential characteristics of gas-solid reaction in micro-fluidized bed reactor[J].Fuel,2013,103(1): 29-36.

    [19]GUO F Q,DONG Y P,FAN P F,et al.Detailed kinetic study of phenol decomposition under isothermal conditions to understand tar catalytic cracking process[J].Journal of Analytical & Applied Pyrolysis,2016,118: 155-163.

    [20]MAO Y B,DONG L,DONG Y P,et al.Fast co-pyrolysis of biomass and lignite in a micro fluidized bed reactor analyzer[J].Bioresource Technology,2015,181: 155-162.

    [21]朱劍虹.微型流化床反應(yīng)分析儀的應(yīng)用[D].湘潭: 湘潭大學(xué),2010.ZHU J H.Application of micro fluidized bed reaction analyzer[D].Xiangtan: Xiangtan University,2010.

    [22]余劍,李強(qiáng),段正康,等.微型流化床中的等溫微分反應(yīng)特性[J].中國(guó)科學(xué):化學(xué),2011,41(1): 152-160.YU J,LI Q,DUAN Z K,et al.Isothermal differential characteristics of the reaction in micro fluidized bed[J].Scientia Sinica Chimica,2011,41(1):152-162.

    [23]GAI C,DONG Y P,LV Z C,et al.Pyrolysis behavior and kinetic study of phenol as tar model compound in micro fluidized bed reactor[J].International Journal of Hydrogen Energy,2015,40(25):7956-7964.

    [24]車德勇.基于物化結(jié)構(gòu)特征的生物質(zhì)與煤共氣化特性研究[D].北京: 華北電力大學(xué),2013.CHE D Y.Research on co-gasification mechanism of biomass and coal based on surface physicochemical structure characteristics[D].Beijing: North China Electric Power University,2013.

    [25]BLA?EJ A,KO?íK M.Phytomass : A Raw Material for Chemistry and Biotechnology[M].Ellis Horwood,1993.

    [26]MOGHTADERI B,MEESRI C,WALL T F.Pyrolytic characteristics of blended coal and woody biomass[J].Fuel,2004,83(6): 745-750.

    [27]王鵬,文芳,邊文,等.煤與生物質(zhì)共熱解特性初步研究[J].煤炭轉(zhuǎn)化,2008,31(4): 40-44.WANG P,WEN F,BIAN W,et al.Study on the co-pyrolysis characteristics of coal and biomass[J].Coal Conversion,2008,31(4):40-44.

    [28]YANG H P,YAN R,CHEN H P,et al.Characteristics of hemicellulose,cellulose and lignin pyrolysis[J].Fuel,2007,86(12/13):1781-1788.

    [29]LIU Q,WANG S R,ZHENG Y,et al.Mechanism study of wood lignin pyrolysis by using TG-FTIR analysis[J].Journal of Analytical& Applied Pyrolysis,2008,82(1): 170-177

    [30]STANIS?AW P.The reactions of formation of selected gas products during coal pyrolysis[J].Fuel,2004,83(9): 1191-1196.

    [31]CAPART R,KHEZAMI L,BURNHAM A K.Assessment of various kinetic models for the pyrolysis of a microgranular cellulose[J].Thermochimica Acta,2004,417(1): 79-89.

    [32]GOLDFARB J L,CEYLAN S.Second-generation sustainability:application of the distributed activation energy model to the pyrolysis of locally sourced biomass-coal blends for use in co-firing scenarios[J].Fuel,2015,160: 297-308.

    [33]閻維平,陳吟穎.生物質(zhì)混合物與褐煤共熱解特性的試驗(yàn)研究[J].動(dòng)力工程學(xué)報(bào),2006,26(6): 865-870.YAN W P,CHEN Y Y.Experimental study on co-pyrolysis characteristics of lignite mixed with biomass mixture[J].Journal of Power Engineering,2006,26(6): 865-870.

    Characterization of co-pyrolysis of pine sawdust and coal slime under isothermal conditions in micro fluidized bed reactor

    LI Tiantao1,GUO Feiqiang1,WANG Yan1,GUO Chenglong1,DONG Yuping2
    (1School of Electrical and Power Engineering,China University of Mining and Technology,Xuzhou221116,Jiangsu,China;2School of Mechanical Engineering,Shandong University,Jinan250061,Shandong,China)

    X 752

    A

    0438—1157(2017)10—3923—11

    10.11949/j.issn.0438-1157.20170393

    2017-04-13收到初稿,2017-05-25收到修改稿。

    聯(lián)系人:郭飛強(qiáng)。

    李天濤(1992—),男,碩士研究生。

    國(guó)家自然科學(xué)基金項(xiàng)目(51406226)。

    Received date:2017-04-13.

    Corresponding author:GUO Feiqiang,fqguo@cumt.edu.cn

    Foundation item:supported by the National Natural Science Foundation of China (51406226).

    猜你喜歡
    木屑煤泥流化床
    木屑和鐵屑的分離實(shí)驗(yàn)
    煤泥重介質(zhì)旋流器在動(dòng)力煤選煤廠粗煤泥分選中的應(yīng)用
    MG-9#捕收劑在極細(xì)粒煤泥浮選中的應(yīng)用
    煤泥水中煤與不同礦物相互作用的模擬研究
    木屑
    飛天(2020年4期)2020-04-17 14:48:24
    流化床丙烷脫氫反應(yīng)段的模擬及優(yōu)化
    關(guān)于循環(huán)流化床鍋爐集控運(yùn)行研究
    煤泥分選超凈煤的藥劑研究
    油頁(yè)巖與木屑混合熱解特性研究
    單沉浸管流化床內(nèi)離散顆粒數(shù)值模擬
    亚洲国产精品成人综合色| 伦精品一区二区三区| 亚洲国产高清在线一区二区三| 国产三级在线视频| 国产一区二区在线av高清观看| 国产在视频线精品| 国产久久久一区二区三区| 97超碰精品成人国产| 一区二区三区高清视频在线| 色哟哟·www| 国内精品宾馆在线| 欧美另类亚洲清纯唯美| 尾随美女入室| 国产精品福利在线免费观看| 视频中文字幕在线观看| 国产精品伦人一区二区| 欧美人与善性xxx| 亚洲精品日韩在线中文字幕| 床上黄色一级片| 亚洲内射少妇av| 国产免费视频播放在线视频 | 欧美bdsm另类| 国产成人免费观看mmmm| 亚洲一级一片aⅴ在线观看| 亚洲av成人精品一二三区| 精品一区二区免费观看| 亚州av有码| 国产精品久久视频播放| 真实男女啪啪啪动态图| 寂寞人妻少妇视频99o| 最近最新中文字幕免费大全7| 毛片一级片免费看久久久久| 欧美丝袜亚洲另类| 女的被弄到高潮叫床怎么办| 日韩欧美国产在线观看| 免费黄网站久久成人精品| 国产乱人偷精品视频| 亚洲国产精品国产精品| 黄色一级大片看看| 99久久精品热视频| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 日本免费在线观看一区| 久久人人爽人人爽人人片va| 亚洲人与动物交配视频| 麻豆一二三区av精品| 国产高潮美女av| 中文精品一卡2卡3卡4更新| 男人和女人高潮做爰伦理| 国产亚洲av片在线观看秒播厂 | 国产精品一区二区性色av| 亚洲精品乱久久久久久| 国产真实乱freesex| 欧美一级a爱片免费观看看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av不卡在线观看| 亚洲美女搞黄在线观看| 人人妻人人澡欧美一区二区| 亚洲精品久久久久久婷婷小说 | 偷拍熟女少妇极品色| 成人无遮挡网站| 亚洲欧美一区二区三区国产| 国产大屁股一区二区在线视频| 国产午夜精品论理片| 纵有疾风起免费观看全集完整版 | 久久热精品热| 亚洲自拍偷在线| 欧美日韩精品成人综合77777| 亚洲国产欧美在线一区| 激情 狠狠 欧美| 一级毛片久久久久久久久女| 国产精品久久久久久精品电影小说 | 一级二级三级毛片免费看| 男女视频在线观看网站免费| 成人无遮挡网站| 免费看av在线观看网站| 美女被艹到高潮喷水动态| 欧美变态另类bdsm刘玥| 联通29元200g的流量卡| 国产成人freesex在线| 精品不卡国产一区二区三区| 欧美一区二区精品小视频在线| 国产av在哪里看| 亚洲美女视频黄频| 激情 狠狠 欧美| 一级毛片我不卡| 能在线免费观看的黄片| 国产一级毛片在线| 老女人水多毛片| 亚洲,欧美,日韩| 亚洲欧美精品专区久久| 免费搜索国产男女视频| 国产高清三级在线| 精品一区二区三区人妻视频| 国产精品三级大全| 婷婷色av中文字幕| 精品国产一区二区三区久久久樱花 | 男女边吃奶边做爰视频| 国产亚洲一区二区精品| 黄色日韩在线| 久久久久久国产a免费观看| 日韩 亚洲 欧美在线| 国产成人午夜福利电影在线观看| 五月伊人婷婷丁香| 久久久久九九精品影院| 国产午夜精品一二区理论片| 免费观看a级毛片全部| 日日摸夜夜添夜夜爱| 国产av一区在线观看免费| 亚洲国产高清在线一区二区三| 久久精品影院6| 一本一本综合久久| 国产成人一区二区在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产毛片a区久久久久| 久久这里只有精品中国| 熟女人妻精品中文字幕| 国产精品99久久久久久久久| 免费av毛片视频| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 亚洲真实伦在线观看| 日本av手机在线免费观看| 免费看日本二区| 日本熟妇午夜| 国产成年人精品一区二区| 国产精品一区www在线观看| 国产精品人妻久久久久久| 精品人妻视频免费看| 国产免费视频播放在线视频 | 国产在线男女| 看非洲黑人一级黄片| 能在线免费看毛片的网站| 成人美女网站在线观看视频| 亚洲欧美精品专区久久| 亚洲欧美清纯卡通| 国产在线男女| 婷婷色综合大香蕉| 2021少妇久久久久久久久久久| 亚洲丝袜综合中文字幕| 亚洲精品成人久久久久久| 能在线免费看毛片的网站| 国产探花在线观看一区二区| 免费观看在线日韩| 国产一区二区在线观看日韩| 亚洲内射少妇av| 97热精品久久久久久| 日韩一区二区三区影片| videos熟女内射| 欧美日本视频| 亚洲第一区二区三区不卡| 亚洲美女视频黄频| 一边亲一边摸免费视频| 青春草亚洲视频在线观看| 最近手机中文字幕大全| 成人二区视频| 99热精品在线国产| 午夜福利在线观看免费完整高清在| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧洲国产日韩| 久久久精品欧美日韩精品| 一边摸一边抽搐一进一小说| 国产精品伦人一区二区| 小蜜桃在线观看免费完整版高清| eeuss影院久久| 插阴视频在线观看视频| 小蜜桃在线观看免费完整版高清| 久久这里只有精品中国| 久久久久性生活片| 插逼视频在线观看| 一个人免费在线观看电影| 国产在线一区二区三区精 | 国产免费视频播放在线视频 | 欧美一区二区精品小视频在线| 搡老妇女老女人老熟妇| 中文在线观看免费www的网站| 男的添女的下面高潮视频| 国产高清国产精品国产三级 | 国产精品国产三级国产av玫瑰| 少妇人妻精品综合一区二区| 中文亚洲av片在线观看爽| a级毛片免费高清观看在线播放| 国产精品伦人一区二区| 永久免费av网站大全| 校园人妻丝袜中文字幕| 国产精品野战在线观看| 99热网站在线观看| 少妇熟女欧美另类| 欧美日韩一区二区视频在线观看视频在线 | 精品少妇黑人巨大在线播放 | 小蜜桃在线观看免费完整版高清| 免费电影在线观看免费观看| 色视频www国产| .国产精品久久| 亚洲内射少妇av| 晚上一个人看的免费电影| 亚洲中文字幕一区二区三区有码在线看| 国产精品精品国产色婷婷| 又爽又黄a免费视频| 日日啪夜夜撸| 1024手机看黄色片| 亚洲经典国产精华液单| 深爱激情五月婷婷| 啦啦啦观看免费观看视频高清| av播播在线观看一区| 一本久久精品| www.av在线官网国产| 国产男人的电影天堂91| 欧美成人午夜免费资源| 日本免费一区二区三区高清不卡| 亚洲国产最新在线播放| 欧美高清成人免费视频www| 亚洲怡红院男人天堂| 最近2019中文字幕mv第一页| 亚洲电影在线观看av| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 毛片女人毛片| 网址你懂的国产日韩在线| 国产一区二区亚洲精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 波多野结衣高清无吗| 国产高清国产精品国产三级 | 色吧在线观看| 久久久久九九精品影院| 真实男女啪啪啪动态图| 亚洲av二区三区四区| 91aial.com中文字幕在线观看| 成人特级av手机在线观看| 水蜜桃什么品种好| 欧美+日韩+精品| 国产成年人精品一区二区| 如何舔出高潮| 亚洲乱码一区二区免费版| 最近手机中文字幕大全| 日本熟妇午夜| 99九九线精品视频在线观看视频| 99久久精品一区二区三区| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 日本黄大片高清| 日日干狠狠操夜夜爽| 成人一区二区视频在线观看| 美女国产视频在线观看| 中文字幕免费在线视频6| 日本免费一区二区三区高清不卡| 日本欧美国产在线视频| 欧美日韩综合久久久久久| 一级二级三级毛片免费看| 亚洲在线自拍视频| 国产成人91sexporn| 成人鲁丝片一二三区免费| 亚洲欧美精品自产自拍| 久久午夜福利片| 亚洲无线观看免费| 亚洲欧洲日产国产| 中文字幕久久专区| 在线免费观看的www视频| 国产精品伦人一区二区| 午夜福利在线在线| 日韩高清综合在线| 免费av不卡在线播放| 成人特级av手机在线观看| 观看免费一级毛片| 免费av毛片视频| 国产一级毛片在线| 大香蕉久久网| 精品免费久久久久久久清纯| 欧美日韩国产亚洲二区| 亚洲精品国产成人久久av| 日本一本二区三区精品| 久久久成人免费电影| 日本黄色片子视频| 精品一区二区三区视频在线| 成年女人看的毛片在线观看| 亚洲国产最新在线播放| 国产精品1区2区在线观看.| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 最新中文字幕久久久久| 久久草成人影院| 亚洲精品日韩在线中文字幕| 最近2019中文字幕mv第一页| 国产爱豆传媒在线观看| 国产精品嫩草影院av在线观看| 国产精品.久久久| 中文字幕熟女人妻在线| 中文字幕久久专区| 七月丁香在线播放| videos熟女内射| www日本黄色视频网| 国产乱来视频区| av在线蜜桃| 国语对白做爰xxxⅹ性视频网站| 免费看光身美女| 久久久精品欧美日韩精品| 联通29元200g的流量卡| 国产亚洲91精品色在线| 内地一区二区视频在线| 国产三级中文精品| 免费观看人在逋| 蜜桃久久精品国产亚洲av| 九九在线视频观看精品| 亚洲自偷自拍三级| 国产黄色视频一区二区在线观看 | 日韩制服骚丝袜av| 97热精品久久久久久| 国产日韩欧美在线精品| 22中文网久久字幕| 亚洲精品国产成人久久av| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美精品v在线| 国产精品美女特级片免费视频播放器| 人体艺术视频欧美日本| 免费观看性生交大片5| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 亚洲欧美日韩高清专用| 97超碰精品成人国产| 黄片wwwwww| 美女cb高潮喷水在线观看| 久久久久久久久中文| 亚洲成人中文字幕在线播放| 亚洲综合精品二区| 国产激情偷乱视频一区二区| av国产免费在线观看| 亚洲精品国产av成人精品| 亚洲成av人片在线播放无| 26uuu在线亚洲综合色| 国产真实乱freesex| h日本视频在线播放| 欧美一区二区国产精品久久精品| 国产乱人偷精品视频| 亚洲欧美成人综合另类久久久 | www.色视频.com| 内地一区二区视频在线| 少妇丰满av| 欧美成人免费av一区二区三区| 国产一区二区三区av在线| 一级黄片播放器| 男的添女的下面高潮视频| 久久久久网色| 精品人妻偷拍中文字幕| 亚洲成人久久爱视频| 亚洲电影在线观看av| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| 国产精品久久久久久av不卡| 亚洲成人av在线免费| 一个人观看的视频www高清免费观看| 亚洲在久久综合| 免费观看人在逋| 欧美+日韩+精品| 视频中文字幕在线观看| www.av在线官网国产| 日日干狠狠操夜夜爽| 日韩,欧美,国产一区二区三区 | 在线天堂最新版资源| 18禁裸乳无遮挡免费网站照片| 免费av不卡在线播放| 好男人视频免费观看在线| 黄色日韩在线| 亚洲综合精品二区| 人妻夜夜爽99麻豆av| 国产精品熟女久久久久浪| 亚洲电影在线观看av| 国产亚洲精品久久久com| 亚洲第一区二区三区不卡| 午夜免费激情av| 日本黄色视频三级网站网址| 麻豆成人午夜福利视频| 22中文网久久字幕| 国产精品国产三级国产专区5o | 久久精品国产鲁丝片午夜精品| 九九久久精品国产亚洲av麻豆| 色尼玛亚洲综合影院| 免费av毛片视频| 中文字幕人妻熟人妻熟丝袜美| 七月丁香在线播放| 欧美97在线视频| 少妇的逼好多水| 天堂av国产一区二区熟女人妻| 亚洲内射少妇av| 亚洲av中文字字幕乱码综合| 国产亚洲午夜精品一区二区久久 | 国产黄片视频在线免费观看| 赤兔流量卡办理| 国产高清国产精品国产三级 | 18禁在线无遮挡免费观看视频| av福利片在线观看| 国内精品宾馆在线| 中文字幕av成人在线电影| 久久精品夜色国产| 亚洲欧美日韩无卡精品| 亚洲美女搞黄在线观看| 日韩欧美三级三区| 国产欧美日韩精品一区二区| or卡值多少钱| 久久鲁丝午夜福利片| 91久久精品国产一区二区成人| 国产精品国产高清国产av| 美女内射精品一级片tv| 成年女人看的毛片在线观看| 亚洲欧美一区二区三区国产| 日韩欧美精品v在线| 亚洲av不卡在线观看| 亚洲美女视频黄频| 一级爰片在线观看| 亚洲国产精品专区欧美| 亚洲一区高清亚洲精品| 成年女人看的毛片在线观看| 最近中文字幕高清免费大全6| 国产伦精品一区二区三区四那| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 亚洲真实伦在线观看| 在现免费观看毛片| 色综合色国产| 色综合站精品国产| 级片在线观看| 美女大奶头视频| 简卡轻食公司| 夫妻性生交免费视频一级片| 中文字幕熟女人妻在线| 精品久久久噜噜| 亚洲欧美精品专区久久| 中文乱码字字幕精品一区二区三区 | 伊人久久精品亚洲午夜| 菩萨蛮人人尽说江南好唐韦庄 | 欧美激情国产日韩精品一区| 亚洲18禁久久av| 国产精品久久久久久久久免| 寂寞人妻少妇视频99o| 久久午夜福利片| 91狼人影院| 欧美一区二区国产精品久久精品| 老女人水多毛片| 日本与韩国留学比较| 国产淫片久久久久久久久| 亚洲经典国产精华液单| 深爱激情五月婷婷| 亚洲不卡免费看| 国产极品天堂在线| 成人性生交大片免费视频hd| 国产老妇女一区| 国产黄色小视频在线观看| 成人午夜精彩视频在线观看| 国产大屁股一区二区在线视频| 亚洲精品久久久久久婷婷小说 | 中文天堂在线官网| 欧美3d第一页| 国产成人91sexporn| 国产免费又黄又爽又色| 欧美性感艳星| 亚洲成色77777| 国产精品嫩草影院av在线观看| 男人舔奶头视频| 如何舔出高潮| 亚州av有码| 国产成人午夜福利电影在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久电影| 热99re8久久精品国产| 五月伊人婷婷丁香| 亚洲人成网站高清观看| 舔av片在线| 三级经典国产精品| 99在线视频只有这里精品首页| 国产亚洲5aaaaa淫片| 亚洲人成网站在线播| 亚洲av免费在线观看| 国产午夜精品一二区理论片| 日韩av在线大香蕉| 成人毛片60女人毛片免费| 国产视频内射| 亚洲国产精品专区欧美| 成人午夜精彩视频在线观看| 丝袜美腿在线中文| 永久免费av网站大全| 干丝袜人妻中文字幕| 国产精品av视频在线免费观看| 99热这里只有是精品50| 深爱激情五月婷婷| 国产熟女欧美一区二区| 久久久久性生活片| 亚洲综合色惰| 欧美高清成人免费视频www| 亚洲国产精品成人综合色| 大话2 男鬼变身卡| 欧美3d第一页| 国产精品国产三级国产av玫瑰| 午夜免费男女啪啪视频观看| 99久久无色码亚洲精品果冻| ponron亚洲| 欧美成人免费av一区二区三区| 亚洲av免费在线观看| 国产精品不卡视频一区二区| 国产精品久久视频播放| 国产三级中文精品| 直男gayav资源| 99在线人妻在线中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 久久国内精品自在自线图片| 国产极品天堂在线| 级片在线观看| 青春草亚洲视频在线观看| 在线观看66精品国产| 三级毛片av免费| 亚洲欧美成人综合另类久久久 | 日日干狠狠操夜夜爽| 天堂中文最新版在线下载 | 最近中文字幕高清免费大全6| 欧美日韩一区二区视频在线观看视频在线 | av福利片在线观看| 精品久久久久久久末码| 网址你懂的国产日韩在线| 少妇猛男粗大的猛烈进出视频 | 三级毛片av免费| 国产爱豆传媒在线观看| 国产精品熟女久久久久浪| 天堂中文最新版在线下载 | 黄色一级大片看看| 亚洲av熟女| 亚洲国产精品合色在线| 亚洲av熟女| 亚洲成色77777| 欧美成人一区二区免费高清观看| 亚洲av二区三区四区| 男插女下体视频免费在线播放| 亚洲欧美精品自产自拍| 亚洲欧美日韩高清专用| 久久久久久九九精品二区国产| 日本黄色片子视频| 国产一区亚洲一区在线观看| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看| 久久国产乱子免费精品| 精品国内亚洲2022精品成人| 亚洲国产精品成人综合色| 精品人妻偷拍中文字幕| 亚洲精品一区蜜桃| 中国国产av一级| 在线观看66精品国产| 三级国产精品片| 丰满人妻一区二区三区视频av| 欧美日韩国产亚洲二区| 好男人视频免费观看在线| 寂寞人妻少妇视频99o| 免费在线观看成人毛片| 亚洲精品成人久久久久久| 日本三级黄在线观看| 日韩制服骚丝袜av| 国产日韩欧美在线精品| 久久精品综合一区二区三区| 欧美日韩国产亚洲二区| 国产成人a∨麻豆精品| 熟女电影av网| av女优亚洲男人天堂| 亚洲精品久久久久久婷婷小说 | 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| 麻豆久久精品国产亚洲av| 中文字幕免费在线视频6| 成人特级av手机在线观看| 精品不卡国产一区二区三区| 干丝袜人妻中文字幕| 又粗又硬又长又爽又黄的视频| 久久鲁丝午夜福利片| 国产精品美女特级片免费视频播放器| 亚洲伊人久久精品综合 | 国产亚洲午夜精品一区二区久久 | 免费看光身美女| 国产精品一区www在线观看| 99久久精品热视频| 亚洲性久久影院| 欧美色视频一区免费| 久久精品久久久久久噜噜老黄 | 亚洲精品影视一区二区三区av| 日本五十路高清| 天堂网av新在线| 99久久无色码亚洲精品果冻| 国产精品国产三级专区第一集| 又黄又爽又刺激的免费视频.| av视频在线观看入口| 男人舔奶头视频| 男女那种视频在线观看| 美女xxoo啪啪120秒动态图| 99久久九九国产精品国产免费| 亚洲自偷自拍三级| 精品久久久久久久末码| 国内精品美女久久久久久| 51国产日韩欧美| 欧美又色又爽又黄视频| 内地一区二区视频在线| АⅤ资源中文在线天堂| 麻豆成人av视频| 99久久无色码亚洲精品果冻| 欧美日韩一区二区视频在线观看视频在线 | 成年免费大片在线观看| 美女大奶头视频| 久久久久久久久久久丰满| 97在线视频观看| 欧美最新免费一区二区三区| 成人av在线播放网站| 深夜a级毛片| 国产成人精品一,二区| 天美传媒精品一区二区| 精品一区二区三区人妻视频| 久热久热在线精品观看| 在线免费观看的www视频| 人人妻人人看人人澡| 91精品伊人久久大香线蕉|